首页/文章/ 详情

我所常用的单元:桁架+梁

1月前浏览167

今天大雨飞机延误,在机场继续这个系列。

梁与桁架单元相比实体单元有更多的假设,单元类型上也有很多的迷惑性。下面简述一下我常用的桁架与梁单元。

桁架:

T3D2:空间2节点桁架单元,3D是空间的意思,2代表两节点,采用线性插值,对于桁架(只承受轴向力)足够了,虽然也有T3D3,但我从来没用过。

梁单元:

1. B31, 空间(3代表3D)一阶Timosenko梁,考虑了截面的剪切变形,梁节点的位移转角均采用线性(1阶)插值,一般用来模拟短粗梁,也可以用来模拟细长梁,当模拟细长梁的时候,精度会比B32与B33低。

2. B32,空间二阶Timosenko梁,考虑了截面的剪切变形,梁节点的位移转角均采用二次插值,可以有效的解决Timosenko梁的剪切自锁问题,一般模拟短粗梁。

3. B33, Euler梁单元,不考虑剪切变形,一般用来模拟细长的梁,由于其不考虑剪切变形,对于短粗梁则得不到准确的结果,由于Euler梁单元采用3次的Hermite插值,对于细长梁往往可以通过较少的单元得到更高的精度。

4. B31OS,B32OS:开放截面的一阶二阶Timosenko梁单元,那么什么叫开放截面呢,开放截面对应的是闭合截面,闭合截面就是我们最常用的那些圆形、矩形、箱型、环形,而对于C型、L型、工字型、T型则属于开口截面梁,理论上需要用带OS的单元类型才能仿真的更为准确,但是现实中很多都直接用B31与B32,并不会过多的关注截面是否为开口截面。不同于正常梁单元的每个节点6个自由度(3个平动、3个转动),开口截面梁每个节点有7个自由度(3个平动、3个转动、1个截面翘曲)。

以上就是我常用的一些梁、桁架单元,之后会详述Timosenko梁单元的剪切自锁与Euler梁的Hermite插值。

对于梁单元的大变形、大转角实现则有非常多的说法,如果大家感兴趣,未来也可以讲一讲这方面的内容。


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

我所理解的张量:基本计算

这次来说一下张量的基本计算。对于加减,标量、向量、张量其实都是用传统的方式进行,就无须赘述,但需要注意的是张量加减的时候需要注意两个张量的“基准”必须是相同的,也就是说张量的各个分量只有两个张量在统一的基向量的情况下才能相加,类比于两组坐标如果想相加,则必须要统一坐标系下。而对于张量乘,则有很多不同的地方,正确的叫法也不应该是乘,而是“积”(product),而又有多种“积”,包括内积,叉积、外积,这里就跟标量区别有些大,在说这些之前,先说一下爱因斯坦标记法,我个人认为其实就是因为科学家懒弄出来的标记方法。爱因斯坦在推导相对论的时候,为了记号方便,就发明了张量表示法。他又发现经常会有相同下标求和的情况发生,于是规定:凡是在一项内有相同下标出现的,一律默认求和,而重复的下标就称为dummyvariable(哑标),不重复的就称为freevariable(自由标)。所以以下公式的意思就是这里熟悉矩阵计算的大家,就可以看到,其实这就非常类似于一个矩阵与向量的乘积。那么下边分别介绍一下内积、叉积、外积。内积是降阶的,这是至关重要的一点,类比两个向量点乘得到一个标量,同理在张量计算中,一个张量为m阶,一个为n阶,且,那么他们内积的结果为阶张量。如下所示,一个2阶张量与2阶张量的内积得到了一个0阶张量,也就是个标量(应变能的计算),另外一个例子则是一个四阶张量与2阶张量的内积得到了一个2阶张量(应力-应变的关系)。叉积不会升降阶,一般叉积其实在计算力学中就是用于两个向量的运算(一阶张量),两个一阶张量的叉积仍得到一个一阶张量那么用张量的技法怎么写呢,两个向量的叉积其实可以表述为下边的形式:其中是Levi-Civita符号,被定义为:外积(outterproduct)也叫kroneckerproduct,tensorproduct,directproduct,就问你晕不晕,不同人有不同的叫法。外积是升阶的,在张量计算中,一个张量为m阶,一个为n阶,那么他们外积的结果为阶,比如如下两个一阶张量的外积有很多人喜欢写外积的时候,直接省略,虽然我非常不喜欢,但是我们得知道这个,所以用这种写法,上边的公式就写成:好了,今天就到这里,下次我们讲一下,求导、拉普拉斯、散度、旋度等。来源:大狗子说数值模拟

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