4.1.3 壳单元质量分配和转动惯量

无论是4节点壳单元，还是3节点壳单元，都是将单元质量m=ρ·A·t 均分到各个节点上，如图4-2所示。

那么壳单元的转动惯量计算公式为

图4-2 4节点壳单元质量分配

用式 （4-17） 计算转动惯量时， 在某些转动速度很大的运动中会出现计算的不稳定 （不收敛）。 所以Belytsch-ko单元 （I_shell=1） 中假定了I_xx=I_yy， 也就是假定单元是个正方形。 当然这和实际网格会产生一些误差， 大多数情况下长方形的单元网格居多， 但是单元长宽比不大的情况下这样的误差也是在可接受范围内的， 尤其这种方法解决了大转动速度下的不稳定问题。 所以划分网格时通常推荐尽量规则划分。