› 金属AM 热变形·残留应力预测模拟器 返回
AM·金属粉末造形

金属AM 热变形·残留应力预测模拟器 — Inherent Strain 法

用于 LPBF·EBM·DED 造形金属部品的反翘（warping）和残留应力的 Inherent Strain 法概略估计工具。通过改变激光功率·扫描速度·层厚来确认能量密度的工艺窗口，可在构建前评价 Ti-6Al-4V 和 Inconel 718 等的热变形风险。

参数设置
材料
杨氏模量 E·热膨胀系数 α·屈服应力 σ_y 自动设定
工艺
LPBF：激光粉末床／EBM：电子束／DED：定向能量
部品代表寸法 L
mm
部品高度 h
mm
激光功率 P
W
扫描速度 v
mm/s
层厚 t
μm
行距以 0.1 mm 固定计算
计算结果
—
能量密度 (J/mm³)
—
工艺品质
—
残留应力 σ_res (MPa)
—
屈服应力比 σ_res/σ_y
—
热变形 δ (mm)
—
构建时间 (h)
—
粉末床·激光扫描·残留应力地图

粉末层上激光扫描，已造形部品积累热应力的过程可视化。色彩表示局部应力（蓝→绿→橙→红）。

工艺地图 — 激光功率 P 和扫描速度 v
材料别残留应力·热变形比较
理论·主要公式

$$E = \frac{P}{v \cdot h \cdot t}, \qquad \epsilon^{*} = \alpha\,(T_m - T_0)\cdot 0.5, \qquad \delta \approx \frac{\epsilon^{*} L^{2}}{h}$$

E：体积能量密度 [J/mm³]、ε*：本征应变（Inherent Strain）、δ：反翘变位、α：热膨胀系数、L：部品代表寸法、h：部品高度。系数 0.5 是简化弹性回复的试验补正。

$$\sigma_{\text{res}} = \min\!\bigl(\sigma_y,\; E\cdot\epsilon^{*}\bigr)$$

残留应力是弹性拘束 E·ε* 和材料屈服应力 σ_y 中较小者。发生屈服时塑性变形截断，应力不再增加。

金属AM 热变形·残留应力预测 — Inherent Strain 法

🙋
金属 3D 打印机是用激光把粉末熔化再凝固，对吧？那为什么部品会反翘、甚至裂纹呢？
🎓
问得好。LPBF（激光粉末床熔融）中，直径 0.1 mm 左右的熔融池以接近 1m/s 的速度移动，同时以每秒 100 万度（10⁶ K/s）的速率冷却。这样的事会重复数百层。每层冷却时，下层的凝固部品会约束它不能自由收缩，那部分「残留应力」就积累起来。最后如果应力超过屈服，去掉支撑的瞬间部品就会裂开，最坏情况下造形中间就会裂。
🙋
这么复杂的现象能用模拟预测吗？如果要逐层逐遍计算，得要好几个月吧…
🎓
完全同意，如果要完全物理地求解熔融池到凝固，会有数百万网格×数十万时间步，对实际零件来说根本不现实。所以就用 Inherent Strain 法（Bugatti 2001）。「每层产生的塑性应变，条件相同的话基本是常数」这个假设，用实验测出「等效本征应变 ε*」，把它作为初始条件给线性弹性 CAE，就能算出全体变形。速度快 100~1000 倍，Ansys Workbench Additive、Autodesk Netfabb 都采用这个。
🙋
所以左边改「激光功率」和「扫描速度」，能量密度就变了，变形量也就变了，对吧？
🎓
完全正确。E = P/(v·h·t) 算出体积能量密度，Ti-6Al-4V 的 LPBF 里 50~80 J/mm³ 是最优。30 以下会产生气孔（lack of fusion），150 以上会蒸发产生键孔气孔。默认 P=250W、v=1200mm/s、t=30μm 的话 E=69.4 J/mm³，评价是「最优」。试试把扫描速度改到 3000 吧，能量会变薄，应该变成「融合不足」。反过来 P=1000W 就是 277 J/mm³，「键孔风险」了。
🙋
改成 AlSi10Mg 的话反翘变大了。为什么？
🎓
热膨胀系数 α 的差异。Ti-6Al-4V 是 α≈9.4×10⁻⁶ /K，但 AlSi10Mg 大约 21×10⁻⁶ /K，2 倍以上。ε* = α·ΔT·0.5，本征应变直接就翻倍，反翘 δ 也成比例增大。AlSi10Mg 在布加迪 Chiron 的发动机脚和宝马 i8 的车顶支架实用化了，但薄壁结构必须加密支撑才能控制反翘。Inconel 718 热膨胀系数中等（13×10⁻⁶），但屈服应力高达 1100 MPa，应力不容易被屈服截断，SpaceX Raptor 喷射器也要必做应力释放热处理（HIP）。
🙋
除了模拟，还有其他办法在实际中减少热变形吗？
🎓
办法很多。最有效的是「构建方向优化」，长边竖起来让 L 变短，δ ∝ L² 就能大幅下降。其次是「加密支撑」，约束增强能压制层间自由收缩。EBM 一开始就 700°C 预热，T_0 高了 ε* 就小——这就是 EBM 在大型 Ti-6Al-4V 医疗植入物上强的原因。还有扫描策略，用岛分割或棋盘式分散热。Magics Simulation 可以生成「逆向变形的补正形状」提前织进去抵消，这是标准做法。

常见问题

Inherent Strain 法是 Bugatti（2001）提出的金属 AM 全体变形预测手法，通过将各层发生的塑性应变作为「等效初始应变」按实验确定，以其为初始条件进行线性弹性 CAE 解析，可高速预测部品整体的反翘和残留应力。由于不直接求解微小尺度的熔融池·凝固，与完全物理热应力解析相比快 100~1000 倍，已被 Ansys Workbench Additive·Autodesk Netfabb·Materialise Magics Simulation·Simufact Additive 等采用。本工具用 ε* = α(T_m - T_0)·0.5 的简易式估计本征应变，用悬臂梁代理 δ ≈ ε*·L²/h 概略估计反翘量。
体积能量密度 E = P/(v·h·t) [J/mm³] 由激光功率 P、扫描速度 v、行距 h、层厚 t 确定，表示「单位体积投入能量」。Ti-6Al-4V 的 LPBF 中一般 50~80 J/mm³ 为最优，30 以下会产生融合不足（lack of fusion）气孔，150 以上会产生蒸发导致的键孔气孔。本工具以 0.1 mm 固定行距计算 E，用 30/80/150 J/mm³ 的阈值将工艺品质分为 4 类：融合不足 / 最优 / 可接受 / 键孔风险。
当残留应力 σ_res 超过屈服应力 σ_y 的 80% 时，在去除支撑或热处理前的搬运过程中部品会反翘，从基板剥离，最坏情况下造形中会产生裂纹。航空发动机部品（GE LEAP 燃油喷嘴）和火箭发动机（SpaceX Raptor 喷射器）在造形后必须进行应力释放热处理（HIP·应力缓解）。本工具返回的 σ_res 是弹性拘束 ε*·E 的上限用 σ_y 截断的值，因此屈服应力比 σ_res/σ_y 超过 0.95 的设计需增加支撑密度或重新检视热处理工程。
根据 δ ≈ ε*·L²/h，首先效果最大的是「缩短代表寸法 L、增加高度 h」。水平面（悬挑）越长越容易反翘，基本原则是通过改变构建方向使 L 变短。其次是通过构建板加热（EBM 的 700°C 预热为代表）提高 T_0 降低 ε*，增加支撑密度提高拘束，或采用岛分割·棋盘式扫描策略分散热积累。AlSi10Mg 这样的热膨胀系数大的材料本质上容易反翘，而 Ti-6Al-4V 和 Inconel 718 这样高温强度高的材料应力不易被屈服截断，易出现高残留应力，需考虑材料特性。

实际应用

航空发动机部品：GE Aviation 的 LEAP 发动机用燃油喷嘴采用 Ti-6Al-4V / Co-Cr LPBF 造形，原来需要 20 个部品组装的结构现在用 1 个部品实现。包含复杂内冷流道的涡轮叶片，支撑除去后如果翘曲超过 0.5 mm 就无法装入发动机，所以造形前需要用 Inherent Strain 解析生成「逆向变形补正形状」，对造形文件进行变形抵消。

火箭发动机·宇宙部品：SpaceX Raptor 甲烷发动机喷射器、Rocket Lab Rutherford 电动泵壳体都采用 Inconel 718 LPBF 制造。由于要承受高温高压循环，造形后必做 HIP（热间等向压加）释放残留应力、压实内部气孔。本工具的「屈服应力比」可作为 HIP 需求的早期筛查指标。

汽车·机动性：布加迪 Chiron 的钛制制动卡钳、宝马 i8 敞篷的 AlSi10Mg 车顶支架、布加迪发动机脚等轻量化和复杂形状并行的部品实用化在推进。AlSi10Mg 热膨胀系数大易反翘，所以通常用拓扑优化确保刚性、用支撑设计释放热应力。

医疗植入物：股骨头、脊椎笼等 Ti-6Al-4V 植入物采用 AM 一体造形多孔结构促进骨长入。EBM（Arcam、GE Additive）的 700°C 预热使残留应力只有 LPBF 的 1/5，特别适合大型植入物。本工具改成 EBM 时，同样条件下变形风险会明显降低。

常见误区和注意事项

首先最大的陷阱是「能量密度 E 相同就造形结果相同」这个误解。E 是 P、v、h、t 乘积的倒数，比如功率翻倍、速度也翻倍的话 E 不变。但实际熔融池长度、深度、热积累、冷速都完全不同。高 P·高 v 容易产生键孔气孔，低 P·低 v 热影响域宽粒子粗化。E 只是「设计起点」，最优点要通过 P 和 v 独立变化的 2D 工艺地图才能找到（参考本工具的工艺地图图表）。

其次是「Inherent Strain 法的应力 σ_res 就是实际应力」的思考固化。本方法是为整体变形预测优化的模型，ε* 必须实验标定。Ansys Additive 默认值中 Ti-6Al-4V LPBF 的 ε* ≈ 0.7~1.5%，但同样机器（EOS M290 vs SLM 280 vs Renishaw RenAM 500）、粉末批次、扫描策略会轻易变化 30%。本工具的简易式适合「相对比较」和「初期检讨」，用于零件合格/不合格判定的绝对值时必须做实机标定。

最后是「支撑加密就能消除反翘」的过度简化。支撑确实约束变形，但那部分应力会在去支撑时转移到部品本体。去掉后会出现反向回弹，支撑痕迹处应力集中容易成为疲劳裂纹起点。实务中是「支撑量最小化」+「逆向变形补正形状织入」组合应对。Materialise Magics Simulation 和 Autodesk Netfabb 标准配置有这个补正功能，Inherent Strain 解析的价值其实在补正形状生成上。

使用指南

  1. 输入部品寸法（mm）、部品高度（mm），从 Ti-6Al-4V、AlSi10Mg、SS316L、Inconel 718、Cu 中选择材料。
  2. 设定激光输出（W）和扫描速度（mm/s），计算能量密度（J/mm³）。一般来说，Ti-6Al-4V 的最优范围是 80~120 J/mm³，AlSi10Mg 是 50~90 J/mm³。
  3. 模拟器运行后，残留应力 σ_res（MPa）、屈服应力比 σ_res/σ_y、热变形 δ（mm）、构建时间（h）实时输出。

具体计算例

Ti-6Al-4V 造形例：部品寸法 50mm、部品高度 30mm、激光输出 200W、扫描速度 1000mm/s 时，能量密度约 100 J/mm³。用 Inherent Strain 法推定残留应力约 250~300MPa，屈服应力比 0.35~0.42，热变形 2.1~2.8mm。构建时间约 4.2 小时。而同条件下 AlSi10Mg（激光输出 150W、扫描速度 1200mm/s）的情况，能量密度 60 J/mm³、残留应力 160~180MPa、热变形 1.4mm，构建时间 3.8 小时。

实务中的注意事项

  1. Inconel 718 在能量密度 110~140 J/mm³ 范围内 HIP 工艺后残留应力可降至 85~120MPa。过度功率设置会导致粗大晶粒，屈服应力比超过 0.5 时裂纹风险急增。
  2. SS316L（奥氏体）热膨胀系数 15.9×10⁻⁶/K，反翘灵敏度高，当 δ>3mm 时需追加支撑结构设计。
  3. Cu 在造形中容易出现 200MPa 以上局部过热残留应力，建议激光功率限制在 150W 以下，扫描速度设在 1500mm/s 以上。
  4. 本模拟器预测值基于平均 LPBF/EBM 条件，实际工艺品质（缺陷率、表面粗糙度）依赖气体控制、喷嘴温度、粉末品质。