Abaqus最终计算出的是积分点上的应力分量,但是在后处理中显示的时候,单元之间会出现应力梯度,为了使contour光滑需要将积分点应力插值到节点处
最终可以以等效应力的方式显示的(通常采用mises应力)或者是其他变量显示
节点的mises应力是它相邻的单元积分点的应力插值过来的,但一个节点通常有多个相邻单元,比如六面体节点有8个相邻单元,但是每个单元插值到这个节点的应力并不相同(通常情况下单元越密
,不同单元插值到同一个节点的应力越接近),
这就需要一种方法来平均一下在节点处得到唯一的应力值,使单元之间应力连续
Abaqus提供了两种插值方式,compute
scalars before averaging和compute scalars after
averaging
compute scalars before
averaging:首先计算出积分点上的等效应力,然后分别插值到节点处,再平均
compute scalars after
averaging:现将积分点处的应力分量插值到节点处,再平均出节点处的应力分量,最后通过应力分量求出等效应力
其中第二种方法不需要设置阀值
第一种方法将积分点处的应力插值到节点处时,不同的单元插值过来的等效应力由于单元等原因可能相差很大的,这样就出现了阀值
averaging
计算公式
relative nodal variation =
(maximum at node - minimum at node)/( maximum over active regions-
minimum over active regions)
relative nodal variation =
(maximum at node - minimum at node)/(maximum within region- minimum
within region)
第一个公式和第二个公式分母部分是选取的不同的参考区域
如果 relative nodal
variation
如果 relative nodal
variation>averaging 则不做平均处理(会在节点处出现应力不连续)
Averaging处于0到100%
具体可参考 abaqus帮助文档CAE部分
我的不严谨的补充(未经严格证实):
(1)采用compute scalars before averaging
当阀值选择为100%的时候,等价于compute
scalars after averaging
(2)一般情况下,两者之间的差距较小,如果两者差距很大,则说明在差距大的地方,你的网格太粗糙。
上述两条未经严格证实,后果自负
|
|