翘曲壳单元的算法公式选择 原创

ELFORM=2

Belytschko-Tsay 壳单元，缺省的壳单元公式，面内单点积分，计算速度很快，通常对于大变形问题是最稳定有效的公式。采用 Co-rotational 应力更新，单元坐标系统置于单元中心，基于平面单元假定，所以对于翘曲的几何体不适用（容易负体积），参考 BWC 壳公式。建议在大多数的分析中使用。

★ 对于几何翘曲问题——此时通过对*CONTROL_SHELL关键字设置——参数 BWC=1 施加翘曲刚度公式，同时参数PROJ=1，以及设置*CONTROL_ACCURACY中参数 INN=2 使节点编号不变，保证计算精度。

★ 对于几何翘曲问题，也可以使用 10 号单元公式：

ELFORM=10

Belytschko-Wong-Chiang 壳单元，面内单点积分，用于单元过度翘曲情况，对于大变形问题没有 B-T 壳稳定。

★ 对于几何翘曲问题，也可以全积分单元替代缩减积分

ELFORM=16

16号公式为全积分壳单元公式（非常快），同样采用Belytschko-Tsay 壳单元公式，采用 Co-rotational 应力更新，采用 2*2 积分点方式，计算精度较高，但时间开销增加并不多，只比 B-T 壳时间多开销 2-3 倍，可以处理翘曲的几何体问题，此时需要激活第 8 种沙漏控制公式【是*CONTROL_HOURGLASS，不是*HOURGLASS】，为了提高计算精度还可以适当增加积分点（NIP）的数量。