首页/文章/ 详情

我所常用的单元:壳单元

2月前浏览737

日常生活中,壳单元的使用场景极其常见,各种钣金件,汽车、楼板等均由壳单元进行模拟。

壳单元可以看做是薄膜单元+板单元的结合体,既能承受平面内的变形,也能承受弯曲。大概就是下图所示的意思:

一般来说我们都认为壳单元有6个自由度,即3个平动3个转动,然而对于很多壳单元来说其实只有5个自由度,即3个平动,与两个弯曲方向的转动,跟上图画的情况有点类似,就是不考虑壳单元法向方向的转动,这也是最为传统的壳单元理论,而更为“新”的一些单元则会考虑3个方向的转动,这也主要是为了大变形非线性。下面就以Abaqus的命名方式为例,说一下我常用的壳单元。Abaqus中有两种壳单元,分别是传统的壳单元与连续体壳单元(Continuum Shell),我对连续体壳单元应用不多,下面就主要说传统的壳单元:

  • S4:4边形壳单元,支持大转角、大应变
  • S4R:带缩减积分的以上单元
  • S3:3角形壳单元,支持大转角、大应变,注意这里他其实是由S4退化而来
  • S3R:带缩减积分的以上单元,而且正是因为S3为S4退化而来,所以这里S3才也有缩减积分的选项
  • STRI3:传统意义上的3角形壳单元,采用薄壳理论,对于壳比较厚的时候可能计算不准确,支持大转角、但不支持大应变
  • S4R5:传统意义上的4边型壳单元,采用薄壳理论,对于壳比较厚的时候可能计算不准确,支持大转角、但不支持大应变

这里细心的朋友们会发现相比于其他类型的单元这里我并没有介绍2阶单元,其实Abaqus中对于壳来说2阶单元也确实存在,诸如STRI65,S8R,S8R5,S9R5等,但在我真实工作中他们极少被使用(如果各位网友使用较多,请留言),这些单元一般来说既计算量大一些,又不见得更为准确(这也是壳单元的特殊之处)。

那么说了这么多壳单元,我们在真实工作中如何使用呢,以我的经验就是S4/S4R+S3/S3R就行了,他们不会增加过多的计算量,且能考虑各种大变形,如果你确定你做的是小变形,或者说线性的薄壳分析,你可以采用S4R5+STRI3

上期说过膜单元的截面属性那些,壳单元比膜单元的截面属性更为丰富,但这里也是要强调的一点就是壳单单元,默认的截面泊松比也是为0.5,也就是说壳单元天生也有不可压缩的假设存在,这也同样解释了壳单元的单元算法中看不到Hybird单元的原因,就是说你不会看到有个单元叫S4H之类的。

还有就是壳单元还可以定义一些钢筋排布的内嵌计算,这就对土木行业非常有用,而对于制造业可能大家从来没点过这个地方,想了解的可以尝试一下,这样通过直接钢筋的排布、面积等信息就能计算出加筋的板壳的承载情况。

这里有一点需要注意的就是这种只能在截面“During analysis”计算下的这个选项下才是可用的,其实这也好理解,因为对于一个分布不确定的钢筋层也很难在计算前(Before analysis)确定其各种参数。

最后说一下显式的壳单元往往跟隐式壳单元的算法 会有一些不同,但这个先按下不表,后续聊他们的公式的时候可以详聊。



来源:大狗子说数值模拟
Abaqus非线性汽车UM理论钣金
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-12-12
最近编辑:2月前
大狗子说数值模拟
博士 传播国际一流的数值模拟算法
获赞 6粉丝 6文章 45课程 0
点赞
收藏
作者推荐

我所理解的张量:为啥需要张量

从事数值仿真、计算物理、计算力学等等相关专业的大家,可能从本科或者研究生开始学习张量分析,我当年学这些的时候就感觉极为痛苦,特别是若干年后看到一本书上,将张量列为“高等数学启蒙小丛书“,再回想当年学习的时候的痛苦经历,我觉得结合我这些年跟张量日夜斗法过程中一些自己对它的理解,可以帮助大家更好的也更轻松的理解张量。今天就以此为开端,给大家说一说我所认为的张量,后续会持续更新的。首先我认为为啥会有张量呢,原因就在于人类本性都很懒,所以尽可能的减轻自己的工作量,科学家也懒,所以在写公式的时候想少写一些,所以为了方便与满足自己懒惰的意愿,张量就诞生了。我们看一下我们通常用的矩阵形式来表示应力:而如果用张量的形式来表示应力张量:这样就很简单,ei与ej分别是两个基准向量,这里可以理解为两个"轴",两个方向向量,比如直角坐标系下的x,y,z轴。这里向量也叫一阶张量,其实也很好理解,因为向量可以表示成一个数(标量)与基准的乘积的形式。那么,我们就可以很天真的认为,向量这种有一个下标的就是个1阶张量,应力这种有两个下标的就是2阶张量,那么有3个下标的是不是就是3阶张量,依此类推,也确实是这样的。那么到现在,科学家为了懒惰而作的努力达到了,可以用一种统一方式表示很多不同的内容。那这个东西未来能方便做研究呢,原因就在于刚刚提到的“轴”,由于科学家在研究中“轴”是变的,基向量是变的。我每次学到这我就想起秦始皇统一度量衡的故事。这里我们可以一起看一下,假如小明身高1.8m,那么他的身高是多少英尺呢,讲到这里大家可能已经明白了,小明的身高自始至终就是那么高,但是在不同的度量体系下那个数值是不一样的。那么推广到高维的几何世界,我们有直角坐标系x,y,z三个方向,我们同样也可以用三个不互相垂直的e1,e2,e3方向向量来表示这个宇宙,当然带来的与之匹配的也就是前边那个带有下标的(标量、矩阵)里面数值不一样,但这个宇宙还是这个宇宙,小明还是小明。这里提一句,张量还跟爱因斯坦有关系,大家都听过广义相对论,“时空“的概念,就是说世界和时间也是有关的,浅显的说你面前的东西,有x,y,z还有个时间t,x,y,z还和时间t有关系,那么也就引入了通常我们所说的“时空”有四个“轴”x,y,z,t,然后张量的表示形式就更为重要了,这里先按下不表,后续详述。那么一个n阶张量就可以这么表示:无论我们怎么改变这里的ei,ej....,这个张量本身是不变的,也就是很多人说的客观的,那么其实我们就可以用不同的形式表达这个自始至终一直在这的物理量(位移、应力、应变、温度梯度、曲率),我们都熟悉的坐标变换,其实也可以认为是这个其中的一种方式,当然这里我们不要求“轴”互相垂直,也不要求“轴”都是单位长度。以下是计算几何和计算力学中一些比较常见的张量,后面我们会一个接着一个用到他们的,并会说到张量表达的好处:1阶:节点位移2阶:应力、应变3阶:位移-应变张量,也就是有限元书中经常用的那个B矩阵4阶:本构,也就是表达应变与应力关系那东西我平时工作中好像没用到更高阶的张量,所以就到这,之后会有张量计算等内容,敬请期待。来源:大狗子说数值模拟

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈