编者按：

熟悉ANSYS Mechanical的朋友知道，早年ANSYS经典界面风行一时，后来随着2000年后ANSYS Workbench平台的推出，经过十多年的发展完善，其易用性、功能性进入了一个非常强健稳定的状态，现在用Workbench平台进行分析的工程师越来越多，毋容置疑的是其易用性远超ANSYS经典界面，在功能角度也实现了相当的水平。早年学习ANSYS的朋友会使用一些经典的练习案例进行学习，熟悉软件的操作及基本特性，那这些经典案例是非常有学习意义的，不过这些官方的经典案例并没有Workbench的版本，所以南京安世亚太集中资源对一些经典的ANSYS学习算例进行了梳理，在workbench中进行了一些复现的尝试，并将以连载的方式与爱好者们分享，希望能对大家的学习工作有所帮助。

基于网格重划分的金属成型仿真分析

案例背景

对复杂的金属成型问题采用有限元分析，可以提高产品质量，同时提高生产效率。但是，

由于结构的大变形、金属的非线性以及非线性的接触，往往大幅度增加了仿真计算的难度，尤其是因为网格变形过大所导致的网格畸形，极易导致计算的失败。网格重划分可以很好地解决由于网格畸形所导致的收敛困难。由于有了网格重划分，即使在计算过程中出现了网格地严重扭曲，程序也会对问题网格进行修复，从而保证计算顺利进行。

采用二维轴对称模型进行该案例地分析。弹性圆柱体为工件，刚性面为模具。其中，动模位于工件上端，以一个微小地速度向下挤压工件，最终，工件成为一个齿圈地形状。

图1 模型介绍

材料输入

工件地材料为AISI8620。杨氏弹性模量为200GPa；泊松比为0.3。

该材料可以近似看成是线弹性材料，且满足等向强化准则。此案例中，将使用幂强化准则(Power-hardening)。

其中为当前屈服强度；为初始屈服强度(此处为385.4MPa)；G为剪切模量；N为幂系数(此处为0.13)；为等效塑性应变。

接触设置

案例中采用无摩擦接触，增强拉格朗日算法，并且每次迭代都更新法向刚度。具体接触对地选择请参照下图：

图2 接触设置

边界条件

如下图所示，在左侧施加无摩擦约束(这里近似可以理解为对称约束)；下端的几何(静模)固定，上端的几何(动模)向下移动135mm。整个工件添加Nonlinear Adaptive Region，表示该区域在计算过程中会进行网格重划分计算。

图3 边界条件

NonlinearAdaptive Region是Workbench中网格重划分功能。该功能完全自动，无需用户在计算过程中进行操作。用户可以在计算前处理时定义网格重划分的准则，重划分的区域，甚至在计算的那些时间区域进行网格重划分。

对于该选项，Workbench中存在一些限制，需要用户了解：

1 该功能目前只能用于结构静力学分析；

2 计算需要打开大变形功能，并且在Output Control中将结果保存设置为All Time Points；

3 不能与一些接触设置(如法向拉格朗日算法，MPC算法……)和边界设置(周期对称边界条件，运动副连接，弹簧连接，远场载荷……)同时使用，具体请参见帮助文档；

4 不能用于混凝土材料、铸铁材料、记忆合金材料以及部分断裂力学材料的计算，具体请参见帮助文档；

5 在二维仿真中，网格重划分仅能用于低阶单元，且三维仿真中，网格重划分仅支持四面体单元，不支持壳单元、梁单元和杆单元，且不支持自接触计算……

NonlinearAdaptive Region的具体设置如下图所示

图4 网格重划分设置

求解设置如下图所示

图5 求解选项设置

结果分析

用户在查看收敛曲线时，可以阅读网格重划分的相关信息。下图橙色实线所对应的时间点，即为网格重划分的时刻。

图6 收敛曲线

此外，在阅读应力变形等结果时，最右侧的Changed Mesh一列也显示了网格重划分的相关信息。

图7 结果列表中的网格重划分信息

齿圈的成型过程如下组图所示：

T=0.0175s

T=0.245s

T=0.50457s

T=0.80234s

T=0.90094s

T=1s

图8 各时间点的工件变形

我们 | 南京安世亚太

▼