有限元理论基础及Abaqus实现方式研究:S4壳单元剪切自锁与沙漏控制

==概述==


在CAE领域,从学校、实验室的自研算法到实现真正的商业化软件是一条无比漫长的道路。我们不研究有限元的新方法、新理论,只是研究商用有限元软件的实现方式。有限元的理论发展了几十年已经相当成熟,商用有限元软件同样也是采用这些成熟的有限元理论,只是在实际应用过程中,商用软件在这些传统的理论基础上会做相应的修正以解决工程中遇到的不同问题,且各家软件的修正方法都不一样,每个主流商用软件手册中都会注明各个单元的理论采用了哪种理论公式,但都只是提一下用什么方法修正,很多没有具体的实现公式。

FEMath.png

商用软件对外就是一个黑盒子,除了开发人员,使用人员只能在黑盒子外猜测内部实现方式,一方面我们查阅Abaqus软件手册得到修正方法的说明,另一方面我们自己编程实现简单的结构有限元求解器,通过自研求解器和Abaqus的结果比较结合理论手册如同管中窥豹一般来研究Abaqus的修正方法,从而猜测商用有限元软件的内部计算方法。为了理解方便有很多问题在数学上其实并不严谨,同时由于水平有限也可能有许多的理论错误,欢迎交流讨论。




==S4壳单元剪切自锁和沙漏控制==

商用有限元软件的健壮性体现在对各种特殊情况,求解过程和解的正确性依然能得到保证,而这些特殊情况在自编程序中如果没有考虑到,那么结果就可能相差极大。其中剪切自锁和沙漏现象是最常见的会影响正确性的两个特殊情况。这两者具有相似性,所以我们在本文中一起研究Abaqus中线性壳单元S4针对这两种情况下的内部实现方式。剪切自锁和沙漏现象影响的是刚度矩阵和应力,我们研究方式是在自编程序iSolver中根据成熟的消除剪切自锁和沙漏控制的理论实现刚度矩阵的修正,通过比较同一模型的Abaqus的刚度矩阵结果,结合帮助文档猜测Abaqus软件单元消除剪切自锁和控制沙漏的内部实现方法。

MeshLock.png

图1:剪切自锁


Hourglass.png

图2:沙漏



===S4壳单元剪切自锁和沙漏控制研究总结===

 完全积分单元才有剪切自锁,虽然Abaqus的S4单元是完全积分,但内部已经做了修正完全消除了剪切自锁,所以不需要用户做任何设置。

 减缩积分单元才有沙漏现象,Abaqus的S4R默认增加一个人工的沙漏刚度来控制沙漏现象,如果发现结果还是不理想,那么需要采用其它建模方法才能控制沙漏了。

 Abaqus针对剪切自锁和沙漏控制做的修正如下:

单元类型元素修正情况
剪切自锁S4薄膜刚度假设应变方法修正
面外弯曲刚度
面外横向剪切刚度假设应变方法修正
S4R所有项
沙漏控制S4薄膜刚度
面外弯曲刚度
面外横向剪切刚度沙漏控制,和Belytchko公式不一致,暂时没研究
S4R薄膜刚度和Belytchko公式一致,因子取0.005G
面外弯曲刚度和Belytchko公式一致,因子取0.00375G*4
面外横向剪切刚度沙漏控制,和Belytchko公式不一致,暂时没研究


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

相关推荐
技术文档
QR Code
微信扫一扫,欢迎咨询~

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

* 公司名称:

姓名不为空

手机不正确

公司不为空