4. 材料模型与工程数据

在 ANSYS 静强度与刚度分析中,几何决定“力往哪里走”,网格决定“计算精度”,而材料模型决定“结构到底有多硬、多强、能撑多久”。如果材料数据输入错误,即使网格再细、边界条件再合理,位移、应力和安全系数也会整体失真。本节围绕线弹性各向同性材料、温度相关属性、塑性基础以及材料数据的输入与校验展开,目标是让读者能够独立完成一套可用于静强度与刚度评估的材料定义。

一、线弹性各向同性材料:静力分析的基本假设

线弹性各向同性材料是静强度分析中最常用的材料模型。它假设应力与应变成正比,卸载后变形完全恢复,且材料在各个方向上的力学性能相同。对于大多数金属材料(如结构钢、铝合金、钛合金)在小变形、常温、远低于屈服点的工况下,这一假设是合理且高效的。

在 ANSYS Workbench 的 Engineering Data 中,线弹性各向同性材料至少需要输入三个参数:

  • 弹性模量 E(Young's Modulus):描述材料抵抗弹性变形的能力,单位常用 MPa 或 GPa。例如结构钢 E ≈ 2.0×10⁵ MPa,铝合金 E ≈ 7.0×10⁴ MPa。
  • 泊松比 ν(Poisson's Ratio):描述横向收缩与纵向伸长的比值,金属通常取 0.27~0.33。结构钢常取 0.3。
  • 密度 ρ(Density):用于考虑自重、惯性载荷或重力加速度工况,单位常用 kg/m³ 或 t/mm³。结构钢 ρ ≈ 7850 kg/m³。

这三个参数直接决定刚度矩阵的组装。E 越大,结构越刚;ν 影响剪切模量 G = E / [2(1+ν)],进而影响扭转与剪切响应;ρ 只在有体积力或惯性载荷时参与计算。若只做静力线性分析且无重力,密度可以暂不填写,但建议一并录入,便于后续扩展。

二、温度相关属性:热-结构耦合与高温工况

当结构工作温度明显偏离室温时,材料属性会随温度变化。例如碳钢在 400 ℃ 时弹性模量约下降 10%~15%,屈服强度下降更明显。若仍按室温数据计算,会高估刚度、低估热应力风险。ANSYS 允许将 E、ν、热膨胀系数 α、屈服强度等定义为温度的函数。

在 Engineering Data 中,温度相关属性的输入方式是:在属性行右侧选择“Tabular”,逐行填入温度与对应数值。例如某低合金钢的弹性模量:

Temperature (℃)    Young's Modulus (MPa)
20                  2.06e5
200                 1.93e5
400                 1.78e5
600                 1.55e5

求解时,ANSYS 会根据节点温度在表格中线性插值。若温度超出表格范围,程序默认采用端点值外推,这可能导致非保守结果,因此表格温度范围应覆盖实际工况。热膨胀系数 α 用于计算热应变 ε_th = α·ΔT,是热应力分析的关键参数;若只做常温静力分析,可不定义 α。

温度相关属性不是“越多越好”。若工况温度变化小于 20 ℃,且材料属性变化在 5% 以内,使用室温常数即可,避免引入不必要的插值误差和计算量。

三、塑性基础:从屈服到强化的静强度评估

线弹性模型在应力超过屈服点后不再适用。若关注的是极限承载能力、局部塑性扩展或残余应力,就需要引入塑性模型。静强度分析中最常用的是双线性等向强化(Bilinear Isotropic Hardening)和多线性等向强化(Multilinear Isotropic Hardening)。

双线性模型只需两个参数:屈服强度 σ_y 和切线模量 E_t。例如某结构钢 σ_y = 345 MPa,E_t = 2.0×10³ MPa。多线性模型则输入一组“塑性应变—应力”数据点,更真实地描述硬化曲线。ANSYS 中塑性数据通常以真实应力—真实塑性应变形式输入,若只有工程应力—应变曲线,需先转换:

真实应力 σ_true = σ_eng × (1 + ε_eng)
真实应变 ε_true = ln(1 + ε_eng)
真实塑性应变 ε_pl = ε_true - σ_true / E

塑性模型直接影响静强度评估结论:线弹性模型只能给出“是否超过屈服”的判断,而塑性模型可以给出“屈服区扩展到哪里、极限载荷是多少”。但塑性分析对网格、载荷步和收敛控制更敏感,若仅做常规刚度校核,线弹性模型通常已足够。

四、材料数据输入与校验:从录入到可信

材料数据输入不是“填完即可”,而应经过校验。建议按以下顺序操作:

  1. 来源确认:优先采用标准(如 GB/T 700、ASTM A36)或供应商实测数据,避免使用来源不明的“经验值”。
  2. 单位统一:ANSYS 不强制单位制,但必须自洽。推荐 mm-t-N-MPa 或 m-kg-N-Pa 体系。若 E 用 MPa、密度用 kg/m³,而几何用 mm,则密度需换算为 t/mm³(7850 kg/m³ = 7.85×10⁻⁹ t/mm³)。
  3. 量级检查:金属 E 通常在 0.7×10⁵~2.1×10⁵ MPa,ν 在 0.2~0.35,ρ 在 2700~8000 kg/m³。若输入值偏离一个数量级,应立即排查。
  4. 曲线合理性:塑性数据中,应力应随塑性应变单调递增,且斜率逐渐减小;若出现下降段,需确认是否为软化模型,否则会导致收敛困难。
  5. 试算验证:用一个简单拉伸试件做单轴验证,检查反力—位移曲线斜率是否与 E·A/L 一致,屈服点是否与 σ_y 对应。
常见错误:将工程应力—应变曲线直接填入多线性塑性表,未转换为真实应力—塑性应变,导致屈服点偏高、硬化段偏硬,静强度评估偏于危险。

五、材料模型选择与静强度评估的对应关系

材料模型的选择应与评估目标匹配:若只关心刚度(位移、变形),线弹性各向同性模型足够;若关心是否屈服、安全系数,需至少输入屈服强度并采用线弹性+屈服判据;若关心极限承载与塑性重分布,则需完整塑性数据。温度相关属性在高温或大温差工况下必须启用。材料数据一旦录入,应通过单位、量级、曲线和试算四道校验,才能进入正式分析。

材料模型选择与数据校验流程 评估目标 刚度 / 位移 变形是否超限 屈服 / 安全系数 是否进入塑性 极限承载 塑性重分布 高温 / 大温差 热应力主导 材料模型 线弹性各向同性 E、ν、ρ 线弹性 + 屈服强度 E、ν、σ_y 多线性塑性 σ_true - ε_pl 曲线 温度相关属性 E(T)、σ_y(T)、α(T) 数据校验 ① 来源确认 标准 / 供应商实测 ② 单位统一 mm-t-N-MPa 自洽 ③ 量级检查 E、ν、ρ 合理范围 ④ 曲线合理性 单调递增、斜率递减 ⑤ 试算验证 单轴拉伸对比 E·A/L 校验不通过 → 返回修正材料数据

完成材料定义后,建议在 Workbench 中通过“材料指派”将模型与几何体关联,并在求解前用“材料数据摘要”核对关键参数。只有材料模型与评估目标匹配、数据经过校验,静强度与刚度分析的结果才具备工程可信度。