为了提高分析精度以及解决在大变形分析中网格畸形的问题,Abaqus提供了三种网格自适应的方法,需注意各自的目标以及适用范围。
![22a57b11bc8ce8a91007919a764cba03.png](https://img-blog.csdnimg.cn/img_convert/22a57b11bc8ce8a91007919a764cba03.png)
本文以一个挤压成型的实例,着重介绍ALE自适应网格和自适应网格重划(adaptive remeshing)。网格间的求解变换(mapsolution)将在下文中介绍。
ALE自适应网格
1. 建立模型
分别建立压头(2D轴对称,解析刚体)和坯料(2D轴对称,变形体)的模型,如下图所示。
![3ae2418b69cffb57a2e3bea12cdfc1cc.png](https://img-blog.csdnimg.cn/img_convert/3ae2418b69cffb57a2e3bea12cdfc1cc.png)
压头和坯料模型
2. 创建坯料材料
创建坯料的材料参数,密度为7.8e-9 t/m3,弹性模量为2e5 MPa,泊松比为0.3,屈服强度为100 MPa,切线模量为300 MPa。
![adbbae3728c3dd6555ec9356bc4f0ebb.png](https://img-blog.csdnimg.cn/img_convert/adbbae3728c3dd6555ec9356bc4f0ebb.png)
塑性参数定义
3.生成装配体(略)
4.设置分析类型
创建静力分析载荷步,打开大变形,并定义合适的迭代增量。
![f31f3315dc378e1c7b802aae8edd68da.png](https://img-blog.csdnimg.cn/img_convert/f31f3315dc378e1c7b802aae8edd68da.png)
step定义
5.定义接触关系
定义压头和坯料上方及右侧的无摩擦接触关系。
![5854830f306a1ffdd985765a582a185e.png](https://img-blog.csdnimg.cn/img_convert/5854830f306a1ffdd985765a582a185e.png)
接触关系定义
6.定义边界条件
胚料底部约束Y方向的位移,左侧约束X方向的位移(对称边界),压头施加向下7.6mm的位移。
![15d1660f6a813175fd9895964258b729.png](https://img-blog.csdnimg.cn/img_convert/15d1660f6a813175fd9895964258b729.png)
约束条件
7.划分网格
定义坯料的单元类型为CAX4R,并在长度方向划分较密的网格,高度方向的网格可以较为稀疏,最终划分的网格如下图所示。
![8c3bbdcb211b7fd53208c0e2578f1017.png](https://img-blog.csdnimg.cn/img_convert/8c3bbdcb211b7fd53208c0e2578f1017.png)
网格划分图
8.建立Job并计算
计算过程中出现error “Time increment required isless than the minimum specified”,计算不收敛。因此下面将改为显式动力学进行分析。
9.显式动力学分析
将step改为“dynamic,explicit”,计算时间为3.8e-4s。同时给予压头向下20000mm/s的恒定速度(替换向下的位移)。重新建立任务提交计算。
最终计算完成,没有出现error,但是网格出现了畸变,这对于计算精度有着较大的影响。
![d794b38cdb044914ab21238202da3fe5.png](https://img-blog.csdnimg.cn/img_convert/d794b38cdb044914ab21238202da3fe5.png)
压缩过程动图(显示动力学)
10. 显式动力学分析+ALE
在step模块中,点击菜单栏other : ALE Adaptive Mesh Domain,选择需要进行ALE网格自适应的区域,并将frequency设为5,remeshing sweeps per increment设为3。(每5个增量步进行一次remesh过程,每一次remesh进行3次sweep)
![07fcee621b137d0960749f9ce97491ac.png](https://img-blog.csdnimg.cn/img_convert/07fcee621b137d0960749f9ce97491ac.png)
ALE设置
进行ALE设置后,重新建立任务提交计算。计算成功完成,并且在整个压缩过程中,网格没有出现畸变。这说明ALE对于控制网格变形有着明显的效果。
![4da939efe6e82e046a45c2c945c8c9cc.png](https://img-blog.csdnimg.cn/img_convert/4da939efe6e82e046a45c2c945c8c9cc.png)
自适应网格重画
如果要采用自适应网格重画功能(adaptive remeshing),需要注意以下两点:
●对于三维实体模型,必须使用四面体单元网格;对于二维模型,必须使用三角形单元或以进阶算法(advancing front)生成的四边形单元网格,否则在提交分析时将会提示错误;
●自适应网格重画根据计算结果误差因子,对误差较大的区域进行网格细化,其根本目的是提高分析精度,而并非解决网格畸形问题。
11. 静力学分析+adaptive remeshing
由于adaptive remeshing只能在ABAQUS/standard求解器中使用,因此分析类型仍设置为static,general。
在mesh模块中,将胚料划分为四边形网格。在菜单栏中选择adaptivity:remeshing rule,,选择胚料区域,进行误差监测变量的设置。
![4fe9c5ee9466008625e68d54c986bea2.png](https://img-blog.csdnimg.cn/img_convert/4fe9c5ee9466008625e68d54c986bea2.png)
Adaptive remeshing 设置
为了节省计算时间,将压头向下移动的距离改为3.8mm。
在job模块中,在菜单栏中选择adaptivity:,创建自适应网格计算任务,保持默认迭代次数3次,提交计算。
迭代次数=1 迭代次数=2 迭代次数=3
![96877142ffc0aa94f9d6ff7f679fc3ea.png](https://img-blog.csdnimg.cn/img_convert/96877142ffc0aa94f9d6ff7f679fc3ea.png)
网格remeshing过程
![9773c50c45c15758b23543b1a0c4f0c9.png](https://img-blog.csdnimg.cn/img_convert/9773c50c45c15758b23543b1a0c4f0c9.png)
Deformation plot
由最终的变形图可以看出,随着网格的细化迭代,模型的变形更加贴近ALE的计算结果,可以认为计算结果更为精确。
本例主要介绍了ABAQUS中ALE和adaptive remeshing这两种自适应网格方法。ALE对于控制网格畸变有着很好的效果,而adaptive remeshing更着重于网格细化,个人认为ALE方法实用性更强。