4. 弹塑性有限元理论

锻造仿真与普通线弹性分析的根本差别在于：材料会屈服、会硬化、会经历大位移与大转动，坯料还会与模具发生不断变化的接触。这四类效应分别对应弹塑性有限元中的增量本构、屈服与返回映射、几何非线性和接触非线性。理解它们各自解决什么问题、在锻造中何时必须开启，是判断仿真结果是否可信的前提。

一、为什么必须用增量本构

锻造过程中应力路径不是简单加载，而是“加载—卸载—再加载”的复杂历史。塑性变形不可逆，应力不能由当前应变唯一确定，只能由应变增量逐步积分得到。因此弹塑性本构写成增量形式：

Δσ = D_ep · Δε
D_ep = D_e − D_e·(∂f/∂σ)(∂f/∂σ)ᵀ·D_e / [ H + (∂f/∂σ)ᵀ·D_e·(∂f/∂σ) ]

其中 D_e 为弹性矩阵，D_ep 为弹塑性切线矩阵，f 为屈服函数，H 为硬化模量。锻造中常用等向硬化假设，H 由流动应力曲线斜率给出。若忽略增量形式而用全量本构，就无法正确反映卸载回弹和残余应力，对锻后尺寸与回弹预测影响显著。

二、屈服准则与返回映射

金属锻造多采用 von Mises 屈服准则，屈服函数为 f = q − σ_y(ε̄ᵖ)，q 为等效应力，σ_y 为当前屈服应力，ε̄ᵖ 为等效塑性应变。当试探应力超出屈服面时，需要把应力“拉回”屈服面，这就是返回映射。

常用径向返回算法分两步：先按弹性预测计算试探应力 σ_tr；再判断 f(σ_tr) 是否大于零。若大于零，则求解塑性乘子 Δλ，使更新后的应力满足屈服条件，并同步更新等效塑性应变与硬化变量。对 von Mises 与等向硬化，Δλ 可由下式迭代求得：

Δλ = f(σ_tr) / (3G + H)
σ = σ_tr − 2G·Δλ·(∂f/∂σ)
ε̄ᵖ ← ε̄ᵖ + Δλ

返回映射的精度直接决定大变形下应力是否“漂移”。锻造单道次压下量常达 20%～60%，若增量步过大，返回映射误差会累积，表现为载荷—行程曲线偏高或偏低。因此实际设置中应控制每步等效应变增量在 0.01～0.05 量级，并在屈服附近适当加密。

三、几何非线性：大变形锻造绕不开

锻造中坯料高度可压缩 50% 以上，网格发生大位移、大转动，应变与位移不再线性相关。几何非线性通过更新拉格朗日或全拉格朗日格式处理，核心是采用 Green-Lagrange 应变或变形率张量，并在刚度矩阵中加入几何刚度项与初应力项。若忽略几何非线性，镦粗时会出现“体积不守恒”的假象，鼓肚形状和接触压力分布都会失真。

实际判断标准是：当最大等效应变超过 0.05～0.1，或单元转动角度超过 5°～10° 时，应开启大变形选项（NLGEOM）。锻造几乎总是满足这一条件，因此几何非线性是默认开启项，而非可选项。

四、接触非线性：锻造仿真的“硬骨头”

坯料与模具的接触面随变形不断变化，接触状态在“分离—黏着—滑动”之间切换，属于典型的边界非线性。常用处理方法有罚函数法和拉格朗日乘子法。罚函数法通过接触刚度 k_n 施加法向约束，k_n 过小会穿透，过大则收敛困难；拉格朗日乘子法精度高但增加自由度。

锻造中还需处理摩擦，常用剪切摩擦模型 τ = m·k，m 为摩擦因子，k 为剪切屈服强度。热锻 m 常取 0.3～0.7，冷锻取 0.05～0.15。摩擦因子对载荷影响很大，例如某轴对称镦粗示例中，m 从 0.3 增至 0.6，相同压下量下成形载荷可上升约 15%～25%。

五、在锻造中的适用场景对照

四类非线性并非彼此独立：几何非线性改变接触搜索的几何位置，接触非线性又影响应力增量路径，而返回映射在每个积分点内完成应力更新。只有把它们放在同一个增量求解框架中协同处理，锻造仿真的载荷、充填与缺陷预测才具备工程参考价值。