回弹模拟器 返回
塑性加工模拟器

回弹模拟器 — 板金弯曲的弹性回复

基于板金纯弯曲后的回弹比 K = 4X³ - 3X + 1(弹性参数 X = σ_y·R_i/(E·t))实时计算卸载后的内侧半径 R_f = R_i/K、90° 弯曲时的回弹角 Δθ = 90·(1-K) 和弹性参数 X,输入参数为降伏应力 σ_y、杨氏模量 E、内侧弯曲半径 R_i 和板厚 t。通过重叠绘制弯曲后(蓝色)和卸载后(红色)的板形状,以及 K-X 曲线的动作点可视化,直观理解回弹的缩放规律。

参数设置
降伏应力 σ_y
MPa
杨氏模量 E
GPa
内侧弯曲半径 R_i
mm
板厚 t
mm

默认值:σ_y = 300 MPa、E = 210 GPa(软钢)、R_i = 50 mm、t = 2.0 mm。X 越大(σ_y 或 R_i 增加,E 或 t 减小),K 越小,回弹越大。薄厚板区域 X 小时 K ≈ 1,回弹接近零。

计算结果
回弹比 K
复原后半径 R_f
90° 回弹角 Δθ
弹性参数 X
板弯曲示意图(弯曲后 → 卸载后)

上部冲头(黄色)与下部凹模(灰色)之间板被弯曲的过程。蓝线 = 加载下的板(R_i),红线 = 卸载后因回弹而恢复的板(R_f = R_i/K)。红线总是开口更大,半径更大。

K-X 曲线(弹性参数与回弹比)

横轴 = X = σ_y·R_i/(E·t),纵轴 = K = 4X³ - 3X + 1。X = 0 时 K = 1(无回弹),X 增大时 K 单调递减。黄色标记 = 当前动作点。K = 0.9、0.7、0.5 辅助线帮助估算"回弹量"的量级。

理论与主要公式

板材纯弯曲后的回弹比通过 Boothroyd/Marciniak 简化式给出:

$$K = \frac{R_i}{R_f} = 4X^3 - 3X + 1, \quad X = \frac{\sigma_y\, R_i}{E\, t}$$

$K$ 是卸载后与卸载前的曲率比,$X$ 是弹性参数(无量纲),$\sigma_y$ 是降伏应力,$R_i$ 是内侧弯曲半径,$E$ 是杨氏模量,$t$ 是板厚。X 越小(厚板、刚性高、降伏应力低),K 越接近 1,回弹越小。

$$R_f = \frac{R_i}{K}, \quad \Delta\theta_{90} = 90^\circ\,(1 - K)$$

卸载后的内侧半径 $R_f$ 按 K 的倒数倍扩大,90° 弯曲时的回弹角 $\Delta\theta$ 与 K 呈线性关系。实际补正中,弯曲角设为 $\theta/K$ 倍进行过度弯曲。

$$\sigma = E\,\varepsilon \;(\varepsilon \le \varepsilon_y),\quad \sigma = \sigma_y \;(\varepsilon \gt \varepsilon_y)$$

所采用的完全弹塑性应力-应变关系:降伏应变 $\varepsilon_y = \sigma_y/E$ 之前为线性,之后为塑性流动。实际材料有加工硬化,K 倾向于接近 1。

回弹模拟器简介

🙋
我用压力机的弯曲工具把板金弯成 90°,但卸下工具后它弯回到约 88°。这是怎么发生的?材料"记得"之前的形状吗?
🎓
很好的观察!这正是回弹现象。板被弯曲时,板厚方向的外侧受拉,内侧受压,产生应力。塑性变形的部分留下永久变形,但接近中性轴的部分仍在弹性范围内。卸下工具时,这些弹性应变释放,弯曲角就恢复了一些。本工具在默认值下(σ_y=300 MPa、E=210 GPa、R_i=50 mm、t=2.0 mm 的软钢),K=0.893,即约 11% 的回弹,90° 弯曲会变成约 80.4°。
🙋
本工具中叫做"弹性参数"的 X = σ_y·R_i/(E·t) 是什么意思?
🎓
X 可以理解为"最外层纤维的弹性应变与板的总应变的比"。当厚度为 t 的板被弯成半径 R_i 时,最外层的全应变约为 ε ≈ t/(2R_i),降伏应变为 ε_y = σ_y/E。X = σ_y·R_i/(E·t) 接近这些比值,X 越大表示"弹性区越广" = 回弹越大,X 越小表示"基本完全塑性化" = 回弹越小。例如,本工具中把板厚从 2 改为 8 mm,X 变为 1/4,K 接近 0.99,回弹角基本消失。
🙋
我听说高强度钢板的回弹特别大。为什么?
🎓
高强度钢的降伏应力 σ_y 很大(590、780、1180 MPa 等),但杨氏模量 E 与普通钢相同(约 210 GPa),所以 X 会跳到普通钢的 2~4 倍。本工具中把 σ_y 从 300 改为 780 MPa,X=0.0929、K=0.722、Δθ=25°,回弹角立即增大。这就是为什么汽车车身用高强度钢时,需要精细计算"过度弯曲 + 金型补正"。相比之下,铝合金 σ_y 中等(约 250 MPa)但 E 很小(70 GPa,仅为钢的 1/3),所以 X 也很大,回弹约为钢的 3 倍。
🙋
实际生产中怎样补正呢?每次都要计算过度弯曲多少吗?
🎓
实务中主要有 3 种方法。(1) 角度补正:目标 90°、K=0.893 时,冲头角设为 90/0.893 ≈ 100.8°。本工具看一眼 K 就能算出来。(2) 半径补正:同时缩小内侧半径,R_punch = R_target × K。(3) 底部冲压/硬币冲压:把冲头压到凹模底部,用高压(3~10 倍)把塑性区扩展到整个板厚,从物理上抑制回弹。本工具主要用于 (1)(2) 的理论计算,底部冲压需要单独建模。但即使这样,初期金型设计和换批次时的试打次数也能大大减少。

常见问题

此公式源自受纯弯曲的线性弹塑性板的弯矩-曲率关系。当板厚为 t、弹塑性应力分布围绕中性轴时,截面弯矩为 M = σ_y·b·(t²/4 - c²/3)/2,其中 c = E·t/(2·σ_y·R_i) 定义弹性核的宽度。卸载时的曲率变化 Δκ 等效于 M·6/(E·b·t³),整理后得到回弹比 K = R_i/R_f = 4·(σ_y·R_i/(E·t))³ - 3·(σ_y·R_i/(E·t)) + 1。Boothroyd 和 Marciniak 等教科书对此有详细推导。
在理论上,纯弯曲、线性弹塑性、无加工硬化的理想条件下,K ≤ 1(K=1 时无回弹)。K > 1 只在 X 为负时出现,但在真实材料参数(σ_y > 0、E > 0、R_i > 0、t > 0)下不可能。本工具出于数值稳定性,在显示时将 K 限制在 [0.01, 1] 范围内。"过度弯曲回弹"的报告通常源于加工硬化或 Bauschinger 效应显著的情况,或弯曲与拉伸并存的深拉工序,这些超出了本公式的适用范围,需要额外建模。
K = R_i/R_f 是独立于弯曲角的比率,所以任意弯曲角都适用。角度的回弹 Δθ 与弯曲角 θ 成正比,公式为 Δθ = θ·(1-K)。例如,本工具默认值 K=0.893 时,30° 弯曲的回弹为 Δθ = 30·0.107 = 3.2°,120° 弯曲为 12.8°。本工具统计栏显示的是 90° 基准的 Δθ,其他角度可用 Δθ_target = θ_target·(1-K) 换算。需要注意的是,极小的角度(10°以下)或极大的角度(170°以上)会使纯弯曲假设失效,摩擦和接触几何影响变大。
通常理论值偏大(K_calc > K_real,即计算对回弹的估计不足),原因有:(1) 加工硬化(σ_y 随变形增加)→ 实际 K 更小,(2) Bauschinger 效应(压缩侧降伏应力降低)→ 同样 K 更小,(3) 凹模肩部摩擦提供拉伸力矩,使 K 接近 1,(4) 板厚在弯曲中变薄,(5) 弹性系数的各向异性,等。反之,在薄厚板的薄板弯曲中,偏离平面应变状态时,理论 K 可能偏小。实际补正的做法是以本工具计算值为初值,通过试打调试微调。

实际应用

汽车车身冲压成形:现代车体为了轻量化采用了大量 590~1180 MPa 级高强度钢板和 6000 系铝合金,这些材料的回弹比普通钢大 2~5 倍,导致车身外板、中柱等的形状精度难以保证。本工具中把 σ_y 改为 590 MPa,K 约为 0.79,Δθ=18.5°,补正量立即增大。奥迪 A8、特斯拉 Model 3 的中柱采用底部冲压和多段弯曲组合,最终精度控制在 ±0.5° 以内。

压力机折弯加工:板金加工厂日常的 V 形折、L 形折中,目标角与成品角的偏差是品质投诉的主要原因。对于 t=1.0 mm、σ_y=250 MPa、R_i=1.0 mm 的薄板 V 折,本工具计算得 X=0.00119、K=0.996,基本无回弹。但改为 t=3.0 mm、R_i=8 mm(典型的自由弯条件)时,X=0.00317、K=0.9905,回弹约 0.85°。进一步把高强度 σ_y=780 MPa 代入,K=0.973,回弹 2.4°,这个值被输入到 CNC 压力机的角度补正表中使用。

轧辊成形的形状管理:把钢卷连续折弯成形的太阳能面板框架或建筑用夹芯板,各轧辊站的折弯角累积起来决定最终截面精度。各站的弯曲角通过 K 计算,轧辊弯曲角设为 1/K 倍进行过度弯曲设计。对于 σ_y=350 MPa(建筑用一般结构钢)、t=2.3 mm、R_i=5 mm,本工具计算 K=0.978,回弹 2°,设计中每个轧辊弯曲 92°,最后成品 90°。

航空器零件的精密弯曲:铝合金(2024-T3、7075-T6)航空外皮肋骨和加强筋,E=70 GPa(低弹性)且 σ_y=300~500 MPa(中等),所以 X 也很大。本工具中设 σ_y=400 MPa、E=70 GPa、R_i=10 mm、t=1.6 mm,得 X=0.0357(与软钢同样数量级!),K=0.893。再加上热处理状态保留和残留应力管理,采用张力成形(边拉伸边弯曲)将 K 逼近 0.95 以上。本工具不涵盖张力效应,但可作为纯弯曲极限参考。

常见误区与注意事项

最常见的误解是"K 由材料特性唯一确定"。但看公式就明白,K 取决于 σ_y、E、R_i、t 的组合 X,即便同一材料,改变弯曲半径和板厚也会大幅改变 K。比如 SPHC 钢(σ_y=300 MPa、E=210 GPa)在 R_i=10 mm × t=4 mm 厚板弯曲时,X=0.00357、K=0.989,基本无回弹;但改为 R_i=100 mm × t=1 mm 薄板大半径弯曲,X=0.143、K=0.583,回弹剧增。建议用本工具拖动滑块,体验"相同材料,不同几何形状,K 完全不同"的效果。数据表上的回弹值只是标准工况的一例。

其次常见的简化错误是"只用过度弯曲能补正"。实际操作中分为 (a) 角度补正:冲头角乘以 1/K;(b) 半径补正:内侧半径乘以 K;(c) 底部冲压:冲头强压到凹模底,用 3~10 倍压力扩大塑性区抑制回弹,三者需视情况搭配。只补正角度而不补正半径,最终内侧 R 可能严重偏差,导致尺寸不合格。本工具的 R_f 值也应参考,金型内半径设计应遵循 R_punch = R_target × K。再者,底部冲压条件下本理论式不适用,需用 PAM-STAMP、AutoForm、LS-DYNA 等有限元工具。

第三个拓展误用是"平面应变板弯曲公式能直接套用 3D 弯曲"。本工具基于平板纯弯曲(无限长、板宽方向应变为零),仅在板厚相对于拉深深度充分小,或长条法兰折弯时成立。对于深拉(Drawing)、法兰化(Flanging)、卷边(Hemming)等,板面内的拉压应力复杂耦合,K 会严重偏离理论值。本工具的适用范围是"压力机 V 形折的角度补正""轧辊成形站点设计"等简单纯弯场景,第一近似可用。复杂成形必须配合实机试模和 CAE 详细分析。

使用指南

  1. 用滑块设置降伏应力 σ_y(MPa)。软钢标准值 250 MPa,硬化铝合金 300 MPa
  2. 调整杨氏模量 E(GPa)。软钢 200 GPa,铝合金 70 GPa
  3. 输入内侧弯曲半径 R_i(mm)和板厚 t(mm),模拟器自动计算弹性参数 X = σ_y·R_i/(E·t)
  4. 回弹比 K = 4X³ - 3X + 1 得出,复原后半径 R_f = K·R_i 和 90° 回弹角 Δθ 实时显示

具体计算案例

SPCC 软钢板(σ_y = 270 MPa、E = 205 GPa、t = 1.5 mm)内侧半径 R_i = 5 mm 弯 90°:X = 270×5/(205000×1.5) = 0.0439,K = 4(0.0439)³ - 3(0.0439) + 1 = 0.869,复原后半径 R_f = 0.869×5 = 4.35 mm,回弹角 Δθ ≈ 4.2°。实务中压力机上死点需调整 -4.2° 偏移。

实操注意事项