3. 刚塑性/刚粘塑性有限元理论

锻造属于典型的大变形塑性成形过程，弹性变形量相对塑性变形量通常可以忽略，因此刚塑性/刚粘塑性有限元成为锻造仿真中最常用的力学模型。与弹塑性有限元同时求解位移和应力不同，刚塑性有限元以速度场为基本未知量，通过变分原理把平衡问题转化为能量泛函的极值问题，再对速度场离散求解。理解这一条主线，是读懂锻造仿真求解器行为、判断结果是否可信的关键。

一、从虚功原理到能量泛函

刚塑性材料在塑性变形阶段满足 Levy-Mises 流动法则，忽略弹性应变后，等效应力与等效应变速率之间通过流动应力关联。对变形体建立虚功原理：外力功率等于内部塑性变形功率与摩擦、惯性等消耗功率之和。在体积不可压缩条件下，引入 Lagrange 乘子 λ（其物理意义为静水压力 p）构造泛函：

Π = ∫_V σ̄·ε̄̇ dV  −  ∫_Sf f·v dS  +  ∫_V λ·ε̇_v dV
其中：σ̄ 等效应力，ε̄̇ 等效应变速率
      ε̇_v = ε̇_x + ε̇_y + ε̇_z 为体积应变速率
      λ = p（静水压力），f 为表面摩擦力

第一项是塑性变形功率，第二项是外力（含摩擦）功率，第三项是体积不可压缩约束。对泛函取驻值 δΠ=0，即可得到同时满足平衡方程、屈服条件和不可压缩条件的解。刚粘塑性模型在此基础上把 σ̄ 写成等效应变速率 ε̄̇ 的函数（如 σ̄ = σ̄(ε̄̇, ε̄, T)），泛函形式不变，只是材料本构中多了速率敏感性。

二、速度场插值与单元离散

求解时把变形体划分为有限个单元，在每个单元内用节点速度插值出内部速度场。以常用的四节点四边形单元为例，单元内任意点的速度为：

v(x,y) = Σ N_i(ξ,η) · v_i
应变速率：ε̇ = B · v_e
其中 B 为应变速率-速度矩阵，v_e 为单元节点速度向量

插值函数 N_i 通常取双线性形函数。由于刚塑性模型不含位移未知量，单元刚度矩阵由 B 矩阵和材料切线模量组装，形成关于节点速度的非线性方程组。等效应变速率 ε̄̇ 由 ε̇ 的分量按 Mises 定义计算，它是连接几何变形与材料流动应力的桥梁。

三、体积不可压缩的两种处理方式

金属塑性变形近似体积不变，即 ε̇_v = ε̇_x + ε̇_y + ε̇_z = 0。若不加约束，纯速度场求解会出现体积漂移，导致结果失真。常用处理方式有两类：

• Lagrange 乘子法：把 λ 作为独立未知量引入泛函，与节点速度同时求解。优点是约束严格满足，λ 直接给出静水压力；缺点是方程组规模增大，刚度矩阵出现零对角块，对求解器要求较高。
• 罚函数法：用大正数 K 乘以 ε̇_v² 作为惩罚项加入泛函，即 Π_罚 = Π + (K/2)∫ε̇_v² dV。当 K 足够大时，ε̇_v 被逼近于零。优点是方程规模不变、实现简单；缺点是 K 取值需权衡：过小则体积不闭合，过大则刚度矩阵病态、收敛困难。工程中常取 K 为材料剪切模量的 10³～10⁶ 倍量级，并通过试算检查体积应变速率是否控制在 10⁻³～10⁻⁴ 量级。

四、求解流程与增量控制

锻造仿真通常按时间增量步推进。每个增量步内，模具按给定速度或位移移动一小段，程序在当前构形上组装并求解速度场，收敛后更新节点坐标与状态变量，再进入下一步。典型流程如下：

1. 读入网格、材料流动应力曲线、边界条件与模具运动规律；
2. 对当前增量步，以上一收敛速度场作为初值；
3. 组装泛函的一阶变分（残差 R）与切线刚度矩阵 K；
4. 求解线性化方程组得到速度修正量 Δv，沿搜索方向更新速度；
5. 检查残差范数与体积应变速率是否满足容差，不满足则回到第 3 步；
6. 收敛后更新节点坐标、等效应变、温度等，输出该步结果；
7. 判断网格是否畸变，必要时触发重划分并插值场变量；
8. 若未达到总行程，进入下一增量步。

增量步长直接影响收敛性。步长过大，速度场变化剧烈，Newton 迭代容易发散；步长过小，计算量增加且可能积累舍入误差。实际中常采用自适应步长：当迭代次数少于 4 次时适当放大步长，超过 8 次仍不收敛则减半重试。对于锻造这类接触状态频繁变化的问题，还需在模具刚接触、填充圆角、飞边形成等阶段主动减小步长。

五、结果判读要点

刚塑性有限元输出的是速度场和应力场，没有弹性回弹信息，因此不能直接用于回弹预测，但非常适合分析金属流动、填充情况和载荷。判读时重点关注：等效应变分布是否与锻造流线一致；静水压力 p 的正负是否合理（压缩区为负，拉应力区为正，后者易引发开裂）；载荷-行程曲线是否出现异常陡升（可能对应飞边充满或网格畸变）；以及体积应变速率是否始终受控。若发现局部 ε̄̇ 异常大或单元出现负 Jacobian，应优先检查网格质量与不可压缩约束参数，而不是直接采信应力结果。