日常生活中,壳单元的使用场景极其常见,各种钣金件,汽车、楼板等均由壳单元进行模拟。
壳单元可以看做是薄膜单元+板单元的结合体,既能承受平面内的变形,也能承受弯曲。大概就是下图所示的意思:
一般来说我们都认为壳单元有6个自由度,即3个平动3个转动,然而对于很多壳单元来说其实只有5个自由度,即3个平动,与两个弯曲方向的转动,跟上图画的情况有点类似,就是不考虑壳单元法向方向的转动,这也是最为传统的壳单元理论,而更为“新”的一些单元则会考虑3个方向的转动,这也主要是为了大变形非线性。下面就以Abaqus的命名方式为例,说一下我常用的壳单元。Abaqus中有两种壳单元,分别是传统的壳单元与连续体壳单元(Continuum Shell),我对连续体壳单元应用不多,下面就主要说传统的壳单元:
这里细心的朋友们会发现相比于其他类型的单元这里我并没有介绍2阶单元,其实Abaqus中对于壳来说2阶单元也确实存在,诸如STRI65,S8R,S8R5,S9R5等,但在我真实工作中他们极少被使用(如果各位网友使用较多,请留言),这些单元一般来说既计算量大一些,又不见得更为准确(这也是壳单元的特殊之处)。
那么说了这么多壳单元,我们在真实工作中如何使用呢,以我的经验就是S4/S4R+S3/S3R就行了,他们不会增加过多的计算量,且能考虑各种大变形,如果你确定你做的是小变形,或者说线性的薄壳分析,你可以采用S4R5+STRI3。
上期说过膜单元的截面属性那些,壳单元比膜单元的截面属性更为丰富,但这里也是要强调的一点就是壳单单元,默认的截面泊松比也是为0.5,也就是说壳单元天生也有不可压缩的假设存在,这也同样解释了壳单元的单元算法中看不到Hybird单元的原因,就是说你不会看到有个单元叫S4H之类的。
还有就是壳单元还可以定义一些钢筋排布的内嵌计算,这就对土木行业非常有用,而对于制造业可能大家从来没点过这个地方,想了解的可以尝试一下,这样通过直接钢筋的排布、面积等信息就能计算出加筋的板壳的承载情况。
这里有一点需要注意的就是这种只能在截面“During analysis”计算下的这个选项下才是可用的,其实这也好理解,因为对于一个分布不确定的钢筋层也很难在计算前(Before analysis)确定其各种参数。
最后说一下显式的壳单元往往跟隐式壳单元的算法 会有一些不同,但这个先按下不表,后续聊他们的公式的时候可以详聊。