ANSYS结构仿真：从网格划分到强度校核的工程实践

ANSYS结构仿真覆盖了从线性静力分析到非线性大变形、疲劳寿命预测的完整链条。对于工程师而言，掌握ANSYS Mechanical的工作流不只是学会操作界面，更关键的是理解每一步设置对最终结果的影响权重——网格质量、边界条件、材料模型，这三者决定了仿真的可信度。

结构仿真的前处理流程

几何模型处理

ANSYS Workbench中的几何模型可以直接从CAD软件导入（支持SolidWorks、Creo、CATIA等），也可以用SpaceClaim或DesignModeler创建。导入后的第一步是特征简化：

• 去除不影响结构强度的小孔、倒角、圆角（这些特征会大幅增加网格数量但不影响整体应力分布）
• 将装配体中的螺栓连接简化为绑定接触或圆柱坐标系下的约束
• 抽取中面将薄壁件转化为壳单元（厚度<1/10整体尺寸时）

过度简化会丢失应力集中区域，简化不足则浪费计算资源。一个实用原则：在关心的应力区域保留几何细节，远离关键区域的特征可以抑制。

网格划分策略

ANSYS Mechanical提供了多种网格类型，选择依据是几何特征和精度需求：

关键指标是网格质量Skewness（歪斜度）和Orthogonal Quality（正交质量）：

• Skewness < 0.8（可接受），< 0.5（良好），< 0.3（优秀）
• Orthogonal Quality > 0.2（可接受），> 0.5（良好），> 0.8（优秀）

应力集中区域（孔边、圆角根部、截面突变处）需要局部加密。ANSYS提供Refinement工具，也可以用Sphere of Influence定义加密球区。加密区域的网格尺寸通常是全局尺寸的1/3到1/5。

接触设置

装配体仿真中接触定义是最容易出错的环节。ANSYS的接触类型：

• Bonded（绑定）：面完全粘接，无相对滑动和分离。适用于焊缝、胶接、螺栓预紧后的连接面。
• No Separation（不分离）：面可以滑动但不会分开。适用于滑轨、导轨。
• Frictionless（无摩擦）：面可自由滑动和分离。适用于间隙配合。
• Rough（粗糙）：面不滑动但可分离。μ=∞。
• Frictional（摩擦）：最真实的接触类型，需定义摩擦系数。适用于螺栓连接、轴承配合。

摩擦接触是非线性的，会显著增加计算时间。如果装配体在工作载荷下接触面不发生相对滑动，用Bonded替代Frictional可以将计算时间缩短3-5倍。

材料模型选择

线性弹性

最基础的材料模型，只需要弹性模量E和泊松比ν。适用于小变形（应变<0.5%）、材料在弹性范围内工作的场景。多数金属结构件的常规分析用这个模型就够了。

双线性等向强化（BISO）

屈服后用切线模量E_t描述硬化行为：

• 屈服强度σ_y
• 切线模量E_t（通常为E的1/100到1/10）
• 需要定义屈服准则：von Mises（金属）、Tresca（保守估计）、Mohr-Coulomb（岩土）

适用于金属的弹塑性分析。BISO模型简单但精度有限——真实金属的硬化曲线是幂律的（Hollomon: σ=K*ε^n），双线性只是分段线性近似。

多线性等向强化（MISO）

用真实应力-应变曲线的多个数据点描述硬化行为，精度比BISO高得多。输入数据需要从工程应力-应变转换为真实应力-应变：

σ_true = σ_eng * (1 + ε_eng) ε_true = ln(1 + ε_eng)

超弹性模型

橡胶、弹性体等大变形材料需要超弹性模型。ANSYS提供Mooney-Rivlin、Ogden、Yeoh、Arruda-Boyce等本构。选择依据是试验数据的类型和变形范围：

• 单轴拉伸数据：Mooney-Rivlin 2参数足够（应变<100%）
• 双轴+剪切数据：Ogden 3阶（应变>100%）
• 大变形（>200%）：Arruda-Boyce或Gent模型

边界条件与载荷施加

边界条件的选择直接影响应力分布的合理性。几个工程经验：

约束过度 vs 约束不足

全约束（固定所有自由度）是最安全的做法，但可能引入不真实的应力集中。更合理的做法是：

• 用远程位移（Remote Displacement）代替固定约束，允许面内旋转
• 用对称边界条件利用对称性（减少模型规模50-75%）
• 用弹簧连接模拟弹性支撑（比刚性约束更真实）

载荷步设置

对于非线性分析（接触、塑性、大变形），载荷需要分步施加：

• 初始子步数：10-20（帮助收敛）
• 最小子步数：5
• 最大子步数：100
• 打开自动时间步长（Auto Time Stepping = On）

如果非线性不收敛，先检查接触设置和网格质量，再考虑减小载荷步长。盲目增加迭代次数通常无效。

后处理与强度校核

应力结果解读

ANSYS输出的应力有多种：

• 等效应力（Equivalent Stress / von Mises）：用于金属的屈服校核，σ_vM < σ_y / n（n为安全系数）
• 主应力（Principal Stresses）：σ_1 > σ_2 > σ_3，用于脆性材料（铸铁、陶瓷）的最大拉应力准则
• 应力强度（Stress Intensity）：σ_1 – σ_3，用于压力容器分析（ASME规范）

网格收敛性验证

结果可信的必要条件是网格收敛——加密网格后应力变化<5%。操作方法：

1. 初始网格计算，记录关心点的应力σ_1
2. 将该区域网格尺寸减半，重新计算得σ_2
3. 再减半得σ_3
4. 如果|σ_3 – σ_2| / σ_2 < 5%，认为收敛

应力集中区域的网格收敛尤其重要——粗网格下圆角根部的应力可能偏低30-50%，直接导致强度评估不安全。

实际工程案例：压力容器封头分析

以DN800椭圆封头为例：

• 材料：Q345R（σ_y=345 MPa，E=206 GPa，ν=0.3）
• 设计压力：2.5 MPa
• 封头壁厚：10 mm
• 网格：Hex Dominant，过渡区加密至2mm
• 接触：封头与筒体焊接面Bonded
• 约束：筒体底部端面全约束
• 载荷：内表面施加2.5 MPa压力

结果：封头过渡区最大von Mises应力约178 MPa，远低于屈服强度345 MPa，安全系数1.94。封头中心区域薄膜应力约125 MPa，与理论公式计算值128 MPa偏差2.3%，验证了仿真精度。