深入解析ABAQUS单元选择:剪切自锁与沙漏现象

在进行网格划分时,用户需根据分析步的类型,网格位移插值函数的阶数选择网格的类型(如图1所示)。通常ABAQUS会默认采用缩减积分单元(Reduced Integration),本文主要介绍相关概念。以下内容参考了《塑性非线性分析原理》等书籍。

ABAQUS中的单元选择-理解剪切自锁和沙漏的图1

图1

一、完全积分Full Integration和缩减积分Reduced Integration的区别

根据有限元理论,对平面4节点单元,如果在一个方向上采用2个高斯积分点,其所计算的刚度矩阵是精确的。因此,4节点4边形的2*2=4个高斯积分点方案为完全积分方案,而缩减积分方案(4节点4边形采用1*1个高斯积分点)不能对刚度矩阵进行精确积分。(图2)

ABAQUS中的单元选择-理解剪切自锁和沙漏的图2

图2

二、完全积分方案中的剪切自锁shear locking现象

当模拟具有弯曲特征的问题时,如果采用完全积分方案单元(如平面4节点应变单元CPE4,平面4节点应力单元CPS4,三维9节点单元C3D8等)进行网格划分,会出现剪切自锁的现象,此时加密网格也不能解决问题。为了说明问题,考察如图3所示的纯弯梁。由于平面4节点单元采用线性的位移插值函数,因此变形后的网格必为一梯形。变形后单元的内角不再保持为直角,也即产生了剪应变,这与实际情况不同。换句话说,这类单元在剪切时过于刚硬,与实际情况不同。由于采用了完全积分方案,刚度矩阵被精确积分,剪切自锁现象将对计算结果产生影响。即使加密网格,计算结果也不会有明显的改善。读者可以自行对一悬臂梁算例进行验证。

需要指出如果采用平面8节点单元,如CPE8,位移在边上2次分布,可以反映弯曲特征,不会产生剪切自锁现象。

ABAQUS中的单元选择-理解剪切自锁和沙漏的图3

图3


三、缩减积分方案中的沙漏现象Hourglass

如果采用缩减积分方案,平面4节点单元只在单元中心点设置一个积分点,该点变形后的剪切应变为0,不会出现剪切自锁现象。未能对刚度矩阵进行精确积分的误差可以通过加密网格来弥补,也即加密缩减积分网格的数量可以有效改善计算结果。


但是同时注意到,单元中心点(积分点)x和y方向上的长度也没有变化(图4),基于该点所计算的整个单元的应变能将等于0,单元表现为零刚度,单元过软,称为沙漏hourglass现象。在网格较粗时,零刚度现象将通过网格进一步扩散,从而产生无意义的结果。为了限制沙漏现象的扩展,ABAQUS引进了“防沙漏刚度”Hourglass stiffness,一般情况下采用默认值即可,如果确有需要可在图1中的Hourglass control选项中设置。

ABAQUS中的单元选择-理解剪切自锁和沙漏的图4
图4


四、小结

如果模型中有比较明显的弯曲现象,为避免出现剪切自锁现象,优先选择二阶单元,或者采用缩减积分方案(网格需要更细,通常厚度方向4层以上)。

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空