2. 金属冲压成形力学基础

冲压成形 CAE 仿真的本质，是把一块金属板料在模具作用下的连续变形过程，离散成有限个单元、有限个时间步上的力学平衡问题来求解。要让求解结果可信，必须先回答三个问题：板料内部每一点的受力状态如何描述？它什么时候开始塑性变形？变形过程中它的抵抗能力如何变化？这三个问题分别对应应力应变张量、屈服准则与硬化模型。再叠加轧制板料的方向性（各向异性）和破裂判据（成形极限图），就构成了冲压仿真从输入到输出的完整力学链条。

一、应力应变张量与坐标约定

板料内任一点的应力状态由二阶对称张量 σ_ij 描述，共 6 个独立分量：3 个正应力 σ_x、σ_y、σ_z 和 3 个剪应力 τ_xy、τ_yz、τ_zx。冲压板料通常处于平面应力状态，即厚度方向应力 σ_z≈0，因此实际需要跟踪的主要是 σ_x、σ_y、τ_xy。应变张量 ε_ij 与之对应，塑性变形中体积近似不可压缩，故 ε_x+ε_y+ε_z=0，厚度应变 ε_z 可由面内两向应变直接推出，这也是仿真中判断减薄率的依据。

坐标约定必须全文统一：取板料轧制方向为 0°（RD），垂直轧制方向为 90°（TD），与轧制方向成 45° 为对角方向。拉应力为正、压应力为负，拉应变为正、压应变为负。仿真后处理中读取的“主应变”即由 ε_x、ε_y、γ_xy 求得的两个主值 ε₁≥ε₂，成形极限图正是以 (ε₂, ε₁) 为坐标绘制。

二、屈服准则：从等效应力到塑性流动

屈服准则回答“应力达到什么组合时开始塑性变形”。金属板料常用两种：

• von Mises 准则：以等效应力 σ̄ 为判据，σ̄ = √(σ_x² − σ_xσ_y + σ_y² + 3τ_xy²)。当 σ̄ 达到当前屈服应力时进入塑性。它假设各向同性，适合铝合金、部分退火钢的近似分析。
• Hill'48 各向异性屈服准则：引入 R₀、R₄₅、R₉₀ 三个塑性应变比，描述板料在不同方向屈服强度的差异。对深冲用低碳钢板，R 值差异显著，必须使用 Hill'48 或 Barlat Yld2000 等准则，否则无法预测制耳和拉深高度。

塑性流动方向由关联流动法则确定：塑性应变增量方向垂直于屈服面法向。这意味着屈服面形状直接决定材料往哪个方向流动，是仿真中“材料往哪跑”的力学根源。

三、硬化模型：屈服面如何长大

进入塑性后，屈服应力随塑性应变增大而提高，即硬化。仿真中常用三种形式：

• 各向同性硬化：屈服面中心不动、均匀扩大，σ̄ = K(ε₀+ε̄_p)ⁿ，即 Hollomon 幂律。参数 K 为强度系数、n 为硬化指数、ε₀ 为预应变。n 越大，均匀变形能力越强，越不易早期颈缩。
• 随动硬化：屈服面中心平移，用于描述包辛格效应，适合反复加载的回弹分析。
• 混合硬化：两者叠加，兼顾循环加载与单调加载，是回弹仿真精度较高的选择。

以某深冲钢板示例（合理假设值）：K=520 MPa，n=0.22，ε₀=0.005，屈服强度约 160 MPa。当等效塑性应变 ε̄_p=0.20 时，σ̄=520×(0.205)^0.22≈520×0.706≈367 MPa。可见从屈服到 20% 塑性应变，流动应力提高了约 2.3 倍，这正是仿真中必须用真实硬化曲线而非单一屈服强度的原因。

四、各向异性与成形极限图

轧制板料的晶粒取向使 R 值随方向变化。R 值定义为宽向塑性应变与厚向塑性应变之比，R>1 表示板料更难减薄，拉深性能更好。R₀、R₄₅、R₉₀ 的差异还会造成圆筒拉深件口部高低不平，即制耳。仿真中若忽略各向异性，制耳高度和拉深力都会被低估。

成形极限图（FLD）是判断破裂的核心工具：横轴为次应变 ε₂，纵轴为主应变 ε₁，曲线为成形极限曲线 FLC。ε₂<0 的左侧为拉-压区，极限应变较高；ε₂>0 的右侧为拉-拉区，极限应变迅速下降。仿真输出的每个单元 (ε₂, ε₁) 点若落在 FLC 上方，即判定破裂风险。工程上常留 10%~20% 安全裕度，即用 FLC 下移后的曲线作为判据。

五、仿真输入输出关系与结果判读

把上述力学要素串起来，冲压仿真的输入输出关系可以概括为：材料参数 + 工艺参数 + 模具几何 → 有限元求解 → 应变/应力/厚度/回弹场。其中材料参数包括 E、ν、硬化曲线、R 值和 FLC；工艺参数包括压边力、摩擦系数、冲压速度、拉深筋阻力；模具几何包括凸凹模圆角、间隙、拉深深度。

结果判读时，先看减薄率：一般冷轧钢板减薄率超过 20%~25% 即接近颈缩，超过 30% 极可能破裂。再看 FLD：把每个单元的 (ε₂, ε₁) 点投到 FLD 上，统计落在安全裕度线以上的单元比例和位置。最后看回弹：回弹量由卸载后的弹性恢复决定，与屈服强度、硬化模型和残余应力分布直接相关，随动硬化模型通常比各向同性硬化预测的回弹更接近实际。

理解这条力学链条，才能在仿真结果异常时快速定位原因：是材料硬化曲线给低了，还是 R 值方向搞反了，还是 FLC 用错了厚度版本。下一阶段将进入网格划分与单元类型选择，但所有网格与求解设置都建立在本节的力学基础之上。