作者caron0225 (淯仔)
看板Cad_Cae
标题[请益] Abaqus 非元素节点 接触 插值
时间Thu Feb 22 14:03:25 2018
各位大大好:
小弟在abaqus上有些问题想请教~
我想要在接触面上，输出一些点的CSTRESS
但是这些点并不在元素的节点上
在abaqus中是否有内建的功能能这样作呢??
起初，小弟有尝试过建立与那些点配对的离散刚体并输出其RF来代替，
但是如果我的点数量要是很多且密，相对离散刚体的元素量也很多且小，
导致计算成本增加，且因为过小的主面造成数值上的误差增加。不是好办法...
而解析刚体只能耦合一个点，只能输出整体的RF。
剩下办法应该就是利用元素节点的cstress值去内插指定点的值了。
虽然我使用的是一阶元素，但是从ODB的结果来看非节点的值不是单纯的线性内插??
还是应该要用积分点去外插??但是好像不能输出积分点的CSTRESS??
还是有其他更好办法??
希望能给小弟一些建议或参考!!感谢!!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 60.248.190.51
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Cad_Cae/M.1519279407.A.026.html
1F：推 gamer: 可以输出surface上的contact force再平均到元素表面上。 02/22 17:53
2F：→ caron0225: 如果一个元素的面有数个点，这样值都一样，我的模型的 02/24 15:51
3F：→ caron0225: 网格相对於点的密度大多了，我也不希望细化模型网格 02/24 15:52
4F：推 gamer: 因为contact force的解是用牛顿法从节点上的力平衡求解来的 02/25 04:10
5F：→ gamer: 所以我猜CSTRESS只是把这个力平均到两节点的Surface上而已 02/25 04:11
6F：推 gamer: 那ABAQUS本身应该没有输出任意点上这个平均值解的预设功能 02/25 04:13
7F：→ gamer: 除非你自己写plug-in做 02/25 04:13
8F：→ gamer: 至於积分点没值是因为contact求解不是从K矩阵来的，所以解 02/25 04:15
9F：→ gamer: 非来自积分点上。 02/25 04:15
的确如G大所说的，积分点是在元素内部，并不在表面，无法计算CANTACT FORCE，用积分
点做插值并不合理。
其实我想计算的是一个面压在一个压力毯上，压力毯是由数个单位矩形组成的，并输出每
个单位矩形上的反力来表示压力。
但是用离散刚体作类似这样的压力毯相对於我的模型来说，元素数量实在太多且密，计算
成本高且收敛精度低。
而仅仅输出各个模型节点上的CSTRESS又不太够，似乎只能自己写plug-in作插值
※ 编辑: caron0225 (60.248.190.51), 02/26/2018 08:30:45
10F：推 gamer: 好奇问一下你要怎麽内插，如果是用相邻节点的值做线性内插 02/27 23:59
11F：→ gamer: 那可以利用X-Y DATA建立分布图，然後勾选interpolate betwe 02/28 00:02
12F：→ gamer: en point选项，就可以用probe选出你要的位置，ABAQUS会帮你 02/28 00:03
13F：→ gamer: 计算两点之间的线性内插值 02/28 00:03
14F：→ caron0225: 原本是打算用element的形状函数去内插~ 02/28 12:43
15F：→ caron0225: G大的这个方法似乎更方便!!我尝试一下!! 02/28 12:43