4. Ansys Mechanical建模准备

在仿真与试验相关性验证中,一个常被低估的事实是:相关性误差的根源往往不在求解器,而在建模准备阶段。几何清理、材料模型选择、单元类型与网格策略这四项准备工作,直接决定了仿真模型能否在物理上"代表"试验件。本节围绕这四项工作,说明它们各自如何影响相关性,以及在实际操作中应如何取舍。

几何清理:先决定"保留什么",再决定"怎么修"

试验件上的倒角、小孔、螺纹、圆角、凸台等特征,在有限元中会带来大量小尺寸单元和应力奇异点。几何清理的核心不是"把模型修干净",而是判断哪些特征对关心的响应量有实质贡献。

以一块带四个安装孔的支架为例,若关注的是整体一阶弯曲频率,安装孔附近的局部应力集中对模态频率影响很小,可将孔简化为等效质量或直接忽略;但若关注孔边裂纹萌生位置,则孔必须保留,且孔边需要局部加密。清理时常用的操作包括:缝合自由边与缝隙、去除微小面与短边、合并重合面、将小圆角改为直角或倒角、抽取中面(壳)或梁线(梁)。

需要特别注意的是,几何清理引入的简化必须在相关性报告中记录。例如把 R2 圆角简化为直角,会使局部刚度略微提高,若试验测得的应变片恰好贴在圆角附近,仿真与试验的偏差就可能主要来自这一简化,而非材料或边界条件。

材料模型选择:线性还是非线性,取决于应变水平

材料模型的选择应基于试验工况下的实际应变水平与温度。对于大多数金属结构在弹性范围内工作的情况,线弹性各向同性模型(弹性模量 E、泊松比 ν、密度 ρ)已足够。但当试验中出现明显塑性变形、大变形、蠕变或橡胶类超弹性行为时,线性模型会系统性高估刚度、低估位移。

下表给出常见材料模型与适用场景的对应关系,可作为选择时的参考:

材料模型关键参数适用场景对相关性的典型影响
线弹性各向同性E、ν、ρ金属弹性范围、小变形频率与刚度相关性好
双线性等向强化塑性屈服强度、切线模量局部屈服、过载试验改善残余应变与极限载荷相关性
超弹性(Mooney-Rivlin 等)C10、C01橡胶减振件显著改善大变形刚度曲线
正交各向异性E1、E2、G12、ν12复合材料层合板决定铺层方向相关性的关键

参数来源同样影响相关性。若 E 取自材料手册而非同批次试样拉伸试验,偏差可达 5%~10%,直接反映为频率偏差。因此在相关性要求较高的项目中,建议用试验件同批次材料做标准拉伸试验,反推 E 与 ν,再代入模型。

单元类型:壳、实体与梁的取舍

单元类型决定了模型对结构力学行为的"表达方式"。薄壁件用壳单元(Shell181)通常比实体单元更高效,但壳单元默认忽略厚度方向应力,若试验关注厚度方向弯曲或层间剪切,则必须用实体单元(Solid186/187)。梁单元(Beam188)适合长细比大的杆件,但连接处的局部柔度需要额外处理。

一个实用的判断原则是:当结构某一方向尺寸远小于其他方向(如厚度/跨度 < 1/10),且不关心该方向梯度时,优先用壳或梁;否则用实体。混合建模时,壳与实体的连接区域需要设置耦合或约束方程,否则会引入虚假的局部刚度。

网格策略:相关性验证中的"收敛"与"局部加密"

网格对相关性的影响体现在两个层面:一是整体网格密度是否足够收敛,二是关心区域是否有局部加密。对于模态分析,通常要求一阶频率在网格加密后变化小于 1%~2%;对于应力热点,热点附近至少布置 3~5 层单元,并避免使用退化单元(如四面体在接触面附近)。

下面这张图概括了建模准备四项工作与相关性验证之间的信号链关系,以及每项工作最需要关注的判断点:

建模准备四项工作 → 相关性验证信号链 几何清理 保留/简化特征 缝合·去小面·抽中面 材料模型 线弹性/塑性/超弹性 E·ν·ρ·屈服·C10 单元类型 壳/实体/梁 厚度方向·连接柔度 网格策略 收敛性·局部加密 热点3~5层单元 有限元模型(刚度矩阵 K、质量矩阵 M、阻尼 C) 仿真响应 频率·振型·位移·应变·应力 试验响应 模态·应变片·位移计 相关性判据 频率误差·MAC·应变偏差 偏差超限 → 回溯建模准备环节 关键判断点 几何简化是否影响关心响应 · 材料参数是否同批次 · 单元是否匹配变形模式 · 网格是否收敛且热点加密

从图中可以看到,四项建模准备工作共同决定了刚度矩阵、质量矩阵与阻尼矩阵,进而决定仿真响应。当仿真与试验出现偏差时,应沿着这条链回溯,而不是直接调整求解设置。

一个可复用的建模准备检查清单

为便于在实际项目中执行,下面给出一份可直接使用的检查清单,覆盖本节四项工作:

  • 几何:是否记录了所有简化特征及其对关心响应的影响?自由边是否缝合?小面/短边是否清理?中面抽取后厚度是否与图纸一致?
  • 材料:E、ν、ρ 是否来自同批次试样?是否根据应变水平选择了线性/非线性模型?温度相关参数是否输入?
  • 单元:壳/实体/梁的选择是否与结构尺寸比匹配?壳与实体连接处是否设置耦合?是否避免了退化单元?
  • 网格:整体网格加密后一阶频率变化是否小于 2%?热点附近是否有 3~5 层单元?接触面网格是否匹配?
建议在相关性报告中为每一项建模准备决策附上简短理由,例如"R2 圆角简化为直角,因应变片距圆角 15 mm,局部应力集中不影响测量值"。这样当后续出现偏差时,可以快速定位是建模简化还是其他因素导致。

建模准备不是一次性的"前处理",而是相关性验证的物理基础。几何清理决定模型代表什么,材料模型决定本构是否真实,单元类型决定变形模式能否被表达,网格策略决定数值解是否收敛。四者中任何一项与试验条件不匹配,都会在相关性对比中表现为难以解释的偏差。因此,在进入求解之前,花时间把这四项工作做扎实,往往比事后反复调参更有效。