2. 锻造过程力学基础

锻造仿真的可信度，首先取决于本构关系是否真实反映了材料在高温、高应变率下的力学行为。有限元求解器每一步都在做同一件事：根据当前应变、应变率和温度，计算应力，再判断材料是否进入塑性。理解这条计算链，才能正确设置材料卡片、读懂等效应力云图，并判断结果是否可信。

一、应力应变：从三维张量到仿真中的等效应力

锻造是典型的三维大变形问题，应力必须用张量描述。以拉应力为正、压应力为负的约定，应力张量可分解为球张量与偏张量：

σij = σm·δij + sij
σm = (σx + σy + σz)/3      （静水压力，影响体积变化）
sij = σij − σm·δij          （偏应力，驱动塑性流动）

金属塑性变形由偏应力驱动，静水压力主要影响压实与损伤。仿真后处理中常见的 von Mises 等效应力，就是把三维应力状态折算为一个标量：

σ̄ = √( ( (σx−σy)² + (σy−σz)² + (σz−σx)² + 6(τxy²+τyz²+τzx²) ) / 2 )

应变同样分弹性与塑性两部分。锻造中塑性应变远大于弹性应变，因此仿真通常采用增量形式累积塑性应变，并用等效应变 ε̄ 表征变形程度。需要特别注意：有限元输出的是对数应变（真应变），而工程应变仅在小变形下近似成立。例如坯料高度由 100 mm 压缩到 60 mm，真应变为 ln(100/60)≈0.511，而工程应变为 0.40，两者相差近 28%，用错会直接导致流动应力取值偏低。

二、屈服准则：判断"何时开始塑性变形"

屈服准则是弹性到塑性的开关。金属材料普遍采用 von Mises 准则：当等效应力达到当前屈服应力时进入塑性。

f = σ̄ − σs(ε̄, ε̇, T) = 0

其中 σs 不是常数，而是随等效应变、应变率和温度变化的流动应力。这一句是整章的核心：屈服面不是固定的，它随变形历史不断扩张。仿真中求解器每步都要求解 f=0，若 f<0 则仍为弹性加载或卸载，若 f=0 且继续加载则发生塑性流动，并按关联流动法则确定应变增量方向。对大多数锻件，von Mises 准则精度足够；对有明显各向异性或织构的板材，则需考虑 Hill 或 Barlat 类各向异性屈服函数。

三、硬化模型：流动应力如何随变形增长

冷锻时加工硬化显著，流动应力随应变上升；热锻时硬化与软化（动态回复、动态再结晶）同时存在，曲线可能先升后降。仿真中常用的硬化表达有三类：

• 幂硬化模型：σ = K·ε̄ⁿ，适合冷锻、室温附近成形，参数少、易拟合。
• Swift / Ludwik 模型：σ = K(ε₀+ε̄)ⁿ，能描述初始屈服段，常用于冷锻与温锻。
• Johnson-Cook 模型：σ = (A + B·ε̄ⁿ)(1 + C·ln ε̇*)(1 − T*m)，把应变硬化、应变率强化和温度软化三项相乘，是热锻与高应变率成形最常用的形式。

以某低碳钢示例参数（合理假设，非实测）：A=350 MPa，B=600 MPa，n=0.45，C=0.015，m=0.8，熔点 1500 ℃。在 ε̄=0.5、ε̇=10 s⁻¹、T=1000 ℃ 时，三项分别约为 350+600×0.5^0.45≈790 MPa、1+0.015×ln(10/1)≈1.035、1−((1000−25)/(1500−25))^0.8≈0.30，相乘得流动应力约 245 MPa。可以看到温度项把应力压低了近 70%，这正是热锻载荷远低于冷锻的原因。

四、温度与应变率效应：热锻仿真的关键变量

温度升高使原子扩散与位错攀移更容易，流动应力下降、延性提高；应变率升高则使位错运动来不及热激活，流动应力上升，即"率敏感性"。二者耦合，常用 Zener-Hollomon 参数统一描述：

Z = ε̇ · exp(Q / (R·T))

Z 越大，流动应力越高。仿真中必须把温度场与变形场耦合：变形生热、接触传热、环境散热同时发生。若只做等温分析，锻件心部与模具接触表面的温差可达 100～200 ℃，流动应力差异会明显改变充填与载荷预测。

五、在仿真中的表达：材料卡片与结果判读

把上述理论落到软件里，通常需要填写：弹性模量与泊松比、密度、热导率与比热、流动应力曲线或本构方程参数、以及摩擦与传热系数。以 Johnson-Cook 为例，材料卡片中要给出 A、B、n、C、m、参考应变率 ε̇₀ 和熔点 Tm。若软件只接受表格形式的流动应力曲线，则需按实际温区（如 900、1000、1100、1200 ℃）和应变率（0.1、1、10 s⁻¹）分别生成曲线族，插值使用。

结果判读时，重点看三处：一是等效应力是否超过材料在该温度下的流动应力上限，若远超则可能是网格畸变或接触穿透导致；二是载荷—行程曲线是否与设备吨位匹配，若峰值载荷接近压机额定吨位，需检查摩擦或温度设置；三是锻件心部与表面的等效应变差，若差异过大，说明变形不均匀，可能产生内部缺陷。只有本构参数、温度场和摩擦条件三者一致，仿真结果才具备工程参考价值。