【JY】有限单元分析的常见问题及单元选择

        因你精彩

即刻关注

    再次整理了下笔记，在看本文前，可以先看下：【JY】有限元分析的单元类型分享一波~

    我们常用的有限元方法有以下非常需要注意的要点（特别是实体单元的应用）：剪切锁死、体积锁死、沙漏模式、零能模式，对于单元选择又需要注意：完全积分、减缩积分、强化应变、杂交分析的概念。

    简而言之，正确的实体单元的模拟是最贴近真实状态的，但是往往由于使用者自身水平受限，包括考虑各种复杂非线性等因素(如材料非线性、几何非线性、接触非线性等)，实体单元是比较难以驾驭的，模型在一定条件下若能简化成杆系单元（【JY】从一根悬臂梁说起）、板壳单元（【JY】板壳单元的分析详解）进行分析是一种非常有利简便、精准快捷的数值计算方法。

☞

剪切锁死（Shear locking）

    概念：在理论上没有剪切变形的单元中发生了剪切变形，该剪切变形也常称伴生剪切。以结构弯曲变形为主的问题中，单元整体出现“刚度过大”的情况，结构形变明显小于预期。

    根据梁的基本理论，对于纯弯变形，轴向应变在厚度方向呈线性变化，而剪应变为0。但如果在计算中采用了低阶完全积分的四节点或四边形单元，线性单元在没有中间节点的情况下无法模拟纯弯。

    在每一个积分点，初始时竖直方向的虚线与水平方向的虚线夹角为90°，变形后却改变了，说明这些点上的剪应力不为零，这与纯弯曲时剪应力为零不一致。产生这种伪剪应力的原因是因为单元的边不能弯曲，他的出现意味着应变能正在产生剪切变形，而不是所希望的弯曲变形，因此总的挠度变小，即单元过硬。剪力自锁仅影响受弯曲载荷完全积分的线性单元行为，而二次单元的边界可以弯曲，故它不存在剪力自锁的问题。

    采用了四节点的线性单元，就只能以上下缘节点相对位移变化来体现“弯曲”但是，纯弯问题的特点是只存在沿高度方向的纤维长短变化，纯弯构件的每个截面与中线总是垂直的。当出现四个节点模拟纯弯的时候，无法体现“中心线的弯曲"，即在数学描述上形成了单元水平线与垂直线之间的夹角变化，即“产生了”名义上的剪应力。从而提高了单元的刚度，然而这部分刚度显然是不存在的，无形中使总的变形量减小，使得分析失真。

受纯弯曲的减缩积分线形单元的位移

     总结发生的条件：1.一阶、全积分单元；2.受纯弯状态(或弯比占大的情况)；

产生的结果：使得弯曲变形偏小，即计算的弯曲刚度太刚，模拟失真。

    剪切锁死的解决方法：

1、采用减缩积分；

2、细化网格；

3、强化应变单元模式，采用非协调单元；

4、假定剪切应变法；

5、采用高阶单元。

    建议采用方法：采用减缩积分，适当细化网格~

    注意：采用高阶单元即增加单元的形函数阶次，就可以在一定程度上消除剪切锁死现象，但由于单元节点增加，计算工作量就提高不少。虽然减缩积分可以在一定程度上消除剪切锁死现象，但有时候效果不佳，有时候还会出现沙漏现象。

☞

体积锁死（Volumetric locking）

    概念：应该有单元的体积变化的时候体积却没发生变化，该原因是受到了伪围压应力（Spurious pressure stresses）。分析结果显示体积几乎不可压缩，体积应变表现为无穷小，体现为结构过硬，甚至导致非线性分析的不收敛。

    材料力学中有弹性模量、剪切模量等等诸多“模量”，其中，有一个与体积有关的模量，即“体积模量"，表达式如下：

    体积模量可描述均质各向同性固体的弹性，可表示为单位面积的力，表示不可压缩性。当泊松比v接近0.5时，上式中分母趋近于零，导致体积模量无穷大、体积应变无穷小。材料表现为不可压缩，在超弹性材料、塑性流动时出现这种不可压缩性的时候，会导致计算困难，产生单元伪应力。（注意：特别橡胶材料）

    选择二阶单元对于弹塑性材料（塑性部分几乎属于不可压缩），二阶全积分四边形和六面体单元在塑性应变和弹性应变在一个数量级时会发生体积锁死，二次减缩积分单元发生大应变时体积锁死也伴随出现。但值得注意的是，一阶全积分单元当采用选择性减缩积分时在变形较小时可以避免出现体积锁死。

    叠层橡胶支座模拟详见：