解决ABAQUS中橡胶大变形问题的有效方法

密封橡胶的数值仿真是一类典型的非线性问题,牵涉到材料非线性(超弹性)、边界非线性(接触)和几何非线性(大变形)问题的集合,如果设置不当,极容易导致求解困难。特别是在密封橡胶的变形复杂,比如和多个不规则边界接触、变形很大等情况,需要更谨慎的设置相关参数,以求得到合适的解答。



模型的适当简化

对薄板问题可忽略厚度方向的应力,作为平面应力(plane stress)问题;对长柱体可忽略第三方向的应变,作为平面应变(plane strain)问题;对O型圈等可作为轴对称问题。平面应力和平面应变在建立part时需选中2D Planar,轴对称问题需选中Axisymmetric;在选择单元时也应注意三者的区别(CPS* ,CPE*, CAX* )。


求解器的选择

因为问题复杂,使用Standard求解容易不收敛,在精度允许的情况下,可选用Explicit求解器。只是多数时候Explicit求解时间较长。应当知道的是,对于橡胶这种典型的不可压缩材料,使用杂交单元(含字母H)是恰当的,但Explicit中没有杂交单元(庄茁书中的例子选用减缩单元)。并且在Explicit中,橡胶材料默认泊松比为0.475。


材料模型的选择

我只用过其中三个,Neo-hookean, 简单易用,就一个参数。对于初学者和简单的模拟比较方便。但是当变形增加到一定范围就不能得到准确的结果了,因为它的参数是来自小变形部分的应力-应变关系。Mooney-Rivlin 是比较常用的本构模型。对于没有加碳黑的橡胶来说,这模型能得到比较准确的结果。但是用它来模拟加了碳黑的橡胶就不太精确了。Yeoh 是用来模拟加碳黑后的橡胶,三个参数都比较容易得到。 可是这个模型在小变形 extension ration=1.5时结果不准确。另外,在橡胶材料模型中引入少量的可压缩性能够减少体积自锁。  


分析步的设置

多设置一个小载荷的分析步,让接触关系平稳的建立起来以后,在下一个分析步中施加更大的真实的载荷,可以减小收敛的困难。  


接触的设置

把刚体面或刚度较大的面作为主面,橡胶一般作为从面。刚性面上的尖角会引起收敛问题,可定义倒角光滑解析刚性面上的尖角。定义接触时应当注意接触面的法向方向应该相反。两个接触面之间滑动较大,故选用有限滑动算法。在ABAQUS/Standard中,对“Constraint enforcement method”的选择,平面应力问题建议选用“Node to surface”,平面应变问题建议选用“Surface to surface”。在结果出来后往往会发现出现穿透现象,这时需要重点检查两个方面:

1.接触定义是否正确,包括主从面、法线方向等;2.从面网格需要细化。


约束的设置

可能的情况下,应当增加约束,可以增大求解的成功的可能,而约束不足容易产生奇异解,比如橡胶“到处乱跑”。比如本模型中,可以对橡胶下方的接触点全约束,而橡胶上方的接触点只允许竖直方向的移动。


网格

网格的疏密、网格划分技术(Mesh controls)和单元类型的选择,对计算的成败影响很大。只能具体问题具体分析,可参考石亦平的《ABAQUS有限元分析实例详解》。划分好网格后最好使用verify来检查一下网格的质量,发生扭曲的地方要改善网格质量。变形严重的地方需要细化网格。优先选用线性单元,高阶单元容易发生单元畸变,易于导致计算失败和不准确的结果。

20.jpg

30.jpg


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

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

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

* 公司名称:

姓名不为空

手机不正确

公司不为空