打开APP
userphoto
未登录

开通VIP,畅享免费电子书等14项超值服

开通VIP
连载(17):剪切闭锁、体积闭锁现象的原因和解决方法

剪切闭锁现象一般发生在出现弯曲变形的线性完全积分单元中(例如 CPS4、CPE4、C3D8)。由于线性单元的直边无法承受弯曲载荷作用,分析过程中可能出现本来不存在的虚假剪应变,使单元的弯曲刚度过大,计算的位移值偏小,即单元的位移场不能模拟由于弯曲而引起的剪切变形和弯曲变形,这就是所谓的“剪切闭锁”现象。当单元长度与厚度的数量级相同或长度大于厚度时,此现象会更严重。

需要注意的是,剪切闭锁现象仅影响受弯曲载荷作用的完全积分线性单元。如果模型受到剪力载荷作用,完全积分线性单元的计算结果则非常好。对于二次单元,由于单元的边上可以发生弯曲变形,一般不会出现剪切闭锁现象。但是如果单元畸变非常严重,或者应力状态非常复杂,存在弯曲应力梯度,二次单元也会出现某种程度的闭锁现象。

只有当载荷所产生的弯矩非常小时,才可以考虑采用完全积分线性单元。如果不能肯定弯曲载荷作用的影响大小,可以尝试使用不同的单元类型,并比较其分析结果,然后选定合适的单元类型。

如果怀疑模型中出现了剪切闭锁现象,可以考虑采用非协调单元或减缩积分单元,但要注意前面几节中介绍过的注意事项。如果模型中网格扭曲非常厉害,仅仅改变单元类型往往不会使计算结果得到很大的改进,划分网格时应尽可能保证单元形状是规则的。

体积闭锁是完全积分单元受到过度约束(overconstraint)时的一种闭锁现象。如果材料为不可压缩或近似不可压缩,完全积分单元可能变得特别刚硬而不会产生体积变形,即所谓的“体积闭锁”。

体积闭锁的一个显著特征是:各个积分点之间或各个单元之间的静水压应力出现急剧变化。在 Visualization 功能模块中绘制静水压应力的云纹图时,如果看到静水压应力从一个积分点到另一个积分点的变化很大,呈棋盘形分布,就有可能出现了体积闭锁的数值计算问题。

如果模型中出现了体积闭锁的数值问题,可以采用下列方法解决:

(1)选取适当的单元类型

如果 Abaqus / Standard 分析中包含不可压缩材料,当塑性应变与弹性应变在同一个数量级上时,二次完全积分单元容易出现体积闭锁现象,往往还伴随着沙漏模式的数值问题。因此不能用于弹塑性分析中。如果必须采用完全积分的二次实体单元,则需要选择这种单元类型的杂交单元形式,但其计算费用将大大增加。

如果使用二次减缩积分单元(例如 C3D20R),当应变大于20% ~ 40%时,需要划分足够密的网格才不会产生体积自锁。

建议使用的单元类型为:非协调单元、线性减缩积分单元或修正的二次四面体单元(C3D10M)。

(2)细化网格

在塑性应变较大的区域应划分足够细化的网格。

(3)引入少量的可压缩性

对于不可压缩材料(泊松比 ),适当引入少量的可压缩性可以减轻体积闭锁现象。由于几乎不可压缩材料和完全不可压缩材料的计算结果很接近,因此可以将不可压缩材料的泊松比取为0.475~0.5之间的值。

本站仅提供存储服务,所有内容均由用户发布,如发现有害或侵权内容,请点击举报
打开APP,阅读全文并永久保存 查看更多类似文章
猜你喜欢
类似文章
【热】打开小程序,算一算2024你的财运
Abaqus基本使用方法
Abaqus单元小结
积分点和单元节点的区别
有限元分析中的剪切锁定(shear lock),ANSYS结果不合理也许就是这货引起的
弹簧设计仿真及疲劳计算
ANSYS/LSDYNA 平时积累手册
更多类似文章 >>
生活服务
热点新闻
分享 收藏 导长图 关注 下载文章
绑定账号成功
后续可登录账号畅享VIP特权!
如果VIP功能使用有故障,
可点击这里联系客服!

联系客服