引言
自适应网格法,简称AMR,抛弃了等距均匀的差分网格,代之以能够自动地适应所研究问题中解的特征的疏密程度不均的曲线网格。
ANSYS中可设置非线性自适应网格,对于解决超弹性材料大变形不收敛问题具有很好的效果。一般出现的errror:网格扭曲也能较好地处理。
详细参数设置及介绍,可参考博客网格重划分技术
超弹性橡胶材料压缩
该案例为典型的超弹性材料压缩,类似于钣金冲压,其中取1/4
模型进行分析。
设置压缩体为刚体,故网格只存在刚体表面,这是为了减少运算成本。
此外,刚体的边界条件不能通过位移(displacement)直接施加,可通过远端位移或运动副的形式。之后会详细介绍。
网格自适应结果展示
网格划分
由于网格自适应不支持3D六面体网格,故插入网格控制,设置网格为四面体划分,采取默认尺寸。上面也提到了,网格只存在刚体表面。
此外,应设置低阶单元,如下图所示。一方面是网格自适应法的要求,另一方面低阶单元多数情况下更适合超弹性体大变形。
边界条件
首先,插入运动副,选择物体-地面,平动副,再选择刚体上表面。
插入运动副边界条件,如下图设置。刚体向下位移200mm,注意图中数据有误。
因为采取1/4模型,故对弹性体对称面插入无摩擦支撑,固定约束弹性体底面。
求解设置
选择求解器为direct,打开大变形,关闭弱弹簧,设置时间子步为100,50,1000.
直接求解
中途发生不收敛,程序报错原因为局部网格扭曲。
混合U-P单元技术
对超弹性体插入一条命令流:
ketopt,matid,6,1
该命令流为处理体积自锁,采取混合U-P单元技术
此外,还可以在求解栏插入体积相容公差值,与混合U-P单元技术对应使用。取值范围为0<vtol<1,默认为1e-5,推荐为1e-5~1e-2
/solu
solc,,,,vtol
由于系统默认该低阶单元不存在该关键字,故需调整为高阶单元,再使用混合U-P单元技术。结果仍不收敛。
网格自适应
插入非线性自适应区域,并进行如下设置。
结果收敛,其仿真结果如上面展示动图。
参考文献
ANSYS命令流学习教程及书籍参考:ANSYS Workbench命令流