Simdroid有限元分析，壳单元的非线性屈曲
一、问题描述几何https://pic2.zhimg.com/80/v2-3f9e15fe7e1905f29a07b60041febb1d_1440w.jpg网格分别采用TMQ单元和MITC4单元求解，网格剖分如下图：https://pic3.zhimg.com/80/v2-57d197cc79d9cfbc024fdb807444232e_1440w.jpg
材料和截面弹性模量7.124e+10 Pa，泊松比 0.3，壳厚度0.0006 m。
分析与荷载壳的 x = 0 m的边固定约束，分为两个分析步：（1）第一步采用线性屈曲分析，在壳的y = -0.24 m 处的边上施加总力为 Fx = 1 N的荷载，计算5阶屈曲特征值和模态。（2）第二步采用非线性屈曲分析，在壳的y = -0.24 m 处的边上施加总力为 Fx = 0.01 N的荷载，注意任何线性分析步的荷载均不会向下传递，停止条件为最大荷载乘子LPF大于300，打开几何非线性。注意本例需要限制最大弧长增量以避免卸载时沿原路径返回。非线性屈曲的控制面板如下图：https://pic1.zhimg.com/80/v2-572947847c541c6a8f84a72f9ace0658_1440w.jpg其中初始缺陷采用了线性屈曲步计算得到的第二个屈曲模态，因为由其结果可知，第一阶屈曲特征值为负数，第二阶才为正值。缺陷值具有长度量纲，意味着采用的初始缺陷最大位移点的位移为1e-4 m。二、计算结果（1）分析步1，线性屈曲的屈曲特征值为：
模态号TMQMITC4
1-0.802722-0.804231
2，临界荷载1.358971.36548
3-2.93052-2.93919
43.79823.80684
5-5.31526-5.35355
可知线性屈曲求得的屈曲临界荷载约为 1.36 N。各阶屈曲模态如下图：https://pic3.zhimg.com/80/v2-95ac1571feba21834d6af3f30628abca_1440w.jpg
（2）分析步2， 非线性屈曲分析，两种单元得到的荷载位移曲线如下图：https://pic4.zhimg.com/80/v2-34952ab41e9508d579cdf30327dd9f17_1440w.jpg可知非线性屈曲求得的临界荷载也为约 1.36 N左右作为参考，此处贴出Abaqus手册中本例的计算结果：https://pic3.zhimg.com/80/v2-6ec9135cc92d71b07d58f723634ff316_1440w.jpg
LPF达到300时的位移云图如下：https://pic4.zhimg.com/80/v2-75898ce1982ac75a6980e7809d9af68f_1440w.jpg三、结果讨论与上一篇 snap through问题不同，本例展示了另一种典型的屈曲类型，即分叉屈曲，在荷载达到临界值时发生的分叉现象。Simdroid软件的几何非线性壳单元可以用来分析这类问题，需要注意的是，必须引入初始缺陷以触发这种分叉屈曲，软件中提供了将线性屈曲模态、振动模态或线性静力的某个工况的结果作为初始缺陷的方法。在Abaqus手册中，荷载被考虑成中间节点上的集中力，这与本例的荷载加在边上略有不同，使得Abaqus的线性屈曲和Riks两个分析步得到的极限荷载均为1.2。如果将集中力均匀分散布置到边的整个节点上，则其线性屈曲分析步的第二阶模态结果为1.3654，两个分析步的结果与本例的结果是一致的。 Simdroid是北京云道智造科技有限公司采用“PaaS大平台+APP小应用”的仿真软件系统架构自主研发的新一代通用CAE平台。该平台能够帮助用户对真实世界中的设备或过程进行仿真分析，实现对产品设计的快速评估和优化；同时，为广大用户提供统一图形交互式的建模仿真和开发环境，支持工程师便捷完成建模仿真和APP开发封装工作，帮助工程师将仿真知识和经验转化为仿真APP，实现自我价值。
[*]自主可控的仿真内核
[*]首创“PaaS大平台+APP小应用”仿真软件系统架构
[*]图形交互式的仿真开发环境
[*]无代码化仿真APP开发、编译工具
扫码试用 基于自主仿真内核的第三代通用CAE及仿真APP开发平台
http://www.simwe.com/uploadfile/2021/1103/20211103041317467.png
页:
[1]