📚 课程目录第4章成形极限预测的理论模型

4. 成形极限预测的理论模型

成形极限曲线（FLC）是判断板材能否安全成形的核心依据，而预测 FLC 的理论模型决定了仿真中失稳判据的物理基础。本节围绕 Swift 分散性失稳、Hill 局部颈缩、M-K 模型以及 Keeler-Brazier 经验公式展开，重点说明各模型的力学假设、参数输入、适用条件与精度差异，并给出在仿真中选用模型的判断依据。

一、从分散性失稳到局部颈缩：两条基本路线

板材在双向拉伸下，失稳并非瞬间发生。Swift 模型认为，当板材的硬化率不足以补偿承载截面的减小，即载荷达到极值点时，发生分散性失稳。其判据可写为：

dσ/dε = σ   （Swift 分散性失稳条件）
σ = K(ε₀ + ε)^n  （幂硬化本构）
→ ε_u = n - ε₀

其中 K 为硬化系数，n 为硬化指数，ε₀ 为预应变。Swift 模型给出的极限应变与 n 值直接相关，形式简洁，适合描述均匀变形阶段的终点。但它假设材料各向同性、变形均匀，无法解释颈缩为何在某一局部突然出现。

Hill 局部颈缩理论则从另一角度出发：当板材在某一方向上的应变硬化率与厚度减薄率满足特定关系时，局部颈缩开始。对于幂硬化材料，Hill 给出的局部颈缩条件为：

dσ/dε = σ / (1 + ρ)   （Hill 局部颈缩条件）
ρ = dε₂/dε₁  （应变路径比）

当 ρ = 0（单向拉伸）时，Hill 条件退化为 dσ/dε = σ，与 Swift 一致；当 ρ = 1（等双向拉伸）时，Hill 条件变为 dσ/dε = σ/2，对应的极限应变明显低于 Swift 预测。这说明 Hill 模型对双向拉伸路径更敏感，预测的 FLC 左半段与右半段呈现不同斜率。

二、M-K 模型：引入初始缺陷的局部颈缩判据

M-K（Marciniak-Kuczyński）模型的核心思想是：板材内部存在一个与主应力方向垂直的初始厚度不均带，其厚度比周围区域薄一个微小量 f₀ = t_b/t_a。在拉伸过程中，凹槽区域的应变增长快于均匀区，当凹槽内外应变比达到临界值时，视为局部颈缩发生。

M-K 模型的关键输入包括：初始不均匀度 f₀（通常取 0.99～0.995）、硬化指数 n、各向异性系数 r、应变路径比 ρ。其计算流程为：

给定 f₀ 和应变路径 ρ，假设凹槽区应变增量 dε_b；
由力平衡条件求解均匀区应变增量 dε_a；
迭代至凹槽内外应变比 dε_b/dε_a 超过临界值（通常取 10）；
记录此时的主次应变，作为该路径下的极限点。

M-K 模型的优势在于能自然预测 FLC 的右半段下降趋势，且对 f₀ 敏感——f₀ 越小，预测的极限应变越高。实际应用中，f₀ 常作为标定参数，用少量试验点反推确定。其局限是计算量较大，且 f₀ 缺乏明确的物理测量手段。

三、Keeler-Brazier 经验公式：快速估算的工程工具

对于低碳钢等常用板材，Keeler 和 Brazier 基于大量试验数据拟合出 FLC 的经验公式。该公式以平面应变状态下的极限应变 FLD₀ 为基准，通过 n 值和厚度 t 进行修正：

FLD₀ = (23.3 + 14.13t) / (0.21 + n)   （t 单位：mm，n 为硬化指数）
FLC 右半段：ε₁ = FLD₀ + ε₂ · (1 - 0.5ε₂)   （近似抛物线）

例如，取 n = 0.22、t = 1.0 mm，则 FLD₀ ≈ (23.3 + 14.13)/(0.21 + 0.22) ≈ 87.0/0.43 ≈ 202%，即平面应变极限应变约为 0.202。该公式计算极快，适合在工艺初筛阶段快速判断成形风险，但它仅适用于 n 值在 0.15～0.25、厚度 0.5～3 mm 的低碳钢，对铝合金、高强钢或非线性应变路径的预测误差较大。

四、模型适用条件与精度对比

四类模型的选用并非互斥，而是根据材料类型、应变路径和精度要求分层使用。下表汇总了关键对比：

模型	核心判据	主要输入	适用条件	精度特点
Swift 分散性失稳	dσ/dε = σ	n、ε₀	均匀变形、单向或近似单向拉伸	左半段尚可，右半段偏高
Hill 局部颈缩	dσ/dε = σ/(1+ρ)	n、ρ	比例加载、各向同性材料	右半段趋势合理，左半段偏保守
M-K 模型	凹槽应变比临界值	f₀、n、r、ρ	比例或非比例加载，可标定	整体精度最高，依赖 f₀ 标定
Keeler-Brazier	经验拟合	n、t	低碳钢、比例加载、厚度 0.5～3 mm	快速估算，外推能力弱

在仿真实践中，若材料为低碳钢且仅需快速评估，可优先用 Keeler-Brazier 公式生成初始 FLC；若材料为铝合金或高强钢，建议采用 M-K 模型，并用单向拉伸与等双向拉伸试验点标定 f₀。对于非线性应变路径（如先拉伸后剪切），M-K 模型可通过增量形式更新应变路径，而 Swift 和 Hill 的解析式则需重新推导。

选用建议：仿真前先确认材料类型与应变路径。比例加载 + 低碳钢 → Keeler-Brazier 或 Hill；比例加载 + 铝合金/高强钢 → M-K；非比例加载 → M-K 增量形式。若仅需左半段安全裕度，Swift 可作为快速下限估计。

需要强调的是，上述模型均建立在比例加载或近似比例加载假设上。当应变路径发生显著偏转（如先单向拉伸再等双向拉伸），Swift 和 Hill 的解析式不再适用，必须采用 M-K 的增量迭代形式，逐段更新凹槽内外应变。此外，M-K 模型中的 f₀ 并非材料常数，而是与板厚、轧制工艺和表面质量相关的标定参数，同一材料不同批次可能需重新标定。

常见误区：直接将 Keeler-Brazier 公式用于铝合金 FLC 预测，会导致右半段极限应变被高估 30% 以上。铝合金的 n 值通常低于 0.2，且应变速率敏感性显著，应改用 M-K 模型并考虑速率相关修正。

综合来看，Swift 和 Hill 提供了失稳判据的力学基础，M-K 模型将初始缺陷引入局部颈缩分析，Keeler-Brazier 则给出了工程快速估算的捷径。在板材成形仿真中，建议以 M-K 模型为主、Keeler-Brazier 为辅进行交叉验证，并根据材料类型和应变路径灵活调整 f₀ 与判据阈值，从而获得与试验更一致的 FLC 预测结果。