第四章：应力与裂纹仿真——热应力与冷裂纹的形成机理
各位工程师朋友，大家好。今天我们来聊聊铸造里最让人头疼的问题之一——应力与裂纹。说实话，我干铸造这行二十多年，见过太多因为裂纹导致废品的情况。有时候一个铸件，前面所有工序都好好的，偏偏在冷却阶段裂了，那叫一个心疼。
这一章，我们就来掰扯清楚：热应力是怎么产生的？冷裂纹又是怎么形成的？以及，我们怎么用仿真手段提前预判这些问题。
4.1 热应力的形成机理
先问大家一个问题：为什么铸件冷却会产生应力？
说白了，就是温度不均匀。铸件不同部位，冷却速度不一样。厚大部位冷得慢，薄壁部位冷得快。冷得快的地方想收缩，但被冷得慢的地方拽着，这就产生了拉应力或压应力。
我个人习惯把热应力分成三个阶段来看：
- 液态阶段：金属还是液体，应力基本为零。液体嘛，你拉它也没用。
- 凝固阶段：这时候开始形成枝晶骨架，应力开始出现。但此时金属还有一定的塑性，能"扛"一下。
- 固态冷却阶段：这是应力积累最严重的阶段。温度越低，材料越脆，应力越大。裂纹往往就发生在这里。
核心观点：热应力的本质是"热胀冷缩不协调"。你想想看，一个铸件里，有的地方想缩，有的地方不让缩，这不就打架了吗？
我在项目中遇到过一种典型情况：某大型机床床身，壁厚差异很大。厚大部位中心温度还很高时，表面已经冷透了。结果表面产生拉应力，直接开裂。后来我们调整了浇注温度和冷却工艺，才把问题解决。
4.2 冷裂纹的形成机理
冷裂纹和热裂纹不一样。热裂纹是在凝固过程中产生的，而冷裂纹是在固态冷却到较低温度时出现的。嗯，这里要注意：冷裂纹往往发生在马氏体转变温度以下，或者是在材料的脆性温度区间。
冷裂纹的形成条件，我总结为三个：
- 足够的拉应力：应力值超过材料的抗拉强度。
- 脆性组织：比如高碳钢中的马氏体，或者铸铁中的渗碳体。
- 氢的存在：氢脆会大大降低材料的断裂韧性。
为什么会这样？因为低温下材料的塑性很差，一旦应力集中，裂纹就会迅速扩展。我曾经处理过一个案例：某高锰钢铸件，热处理后出现大量微裂纹。查了半天，发现是冷却速度太快，产生了大量马氏体，加上残余应力没释放，结果一碰就裂。
避坑指南：我曾经在调试一个大型铸钢件的工艺时，忽略了氢的影响。结果铸件在放置24小时后，突然出现延迟裂纹。从那以后，我每次做高强钢铸件，都会严格控制原材料的氢含量，并增加去氢退火工序。
4.3 应力场仿真设置
好了，理论讲完了，咱们来点实际的。应力场仿真怎么做？我用的是ProCAST和Abaqus联合仿真的方式。当然，其他软件如Magma、AnyCasting也类似。
我个人习惯的仿真流程是这样的：
- 第一步：先做温度场仿真。应力场是温度场的结果，温度算不准，应力就别想了。
- 第二步：设置材料属性。这里要注意，高温下的力学性能数据很难找。我一般用JMatPro计算，或者查文献。
- 第三步：定义边界条件。包括铸件与铸型的接触、冷却条件、约束条件等。
- 第四步：设置求解参数。时间步长要足够小，尤其是相变阶段。
下面是一个简单的ProCAST应力仿真设置示例：
// 应力分析设置
$ STRESS ANALYSIS
STRESS ON
MECHANICAL MODEL ELASTIC_PLASTIC
YIELD CRITERION VON_MISES
HARDENING ISOTROPIC
// 材料属性
MATERIAL STEEL_1045
YOUNG_MODULUS 210000 MPa
POISSON_RATIO 0.3
THERMAL_EXPANSION 1.2e-5 1/K
YIELD_STRESS 350 MPa
TANGENT_MODULUS 2000 MPa
END
// 约束条件
FIXED NODE 1001 X Y Z
FIXED NODE 1002 Y Z
// 冷却条件
HEAT_TRANSFER 500 W/m2K
你可能会问：这些参数怎么取？我的经验是，弹性模量和泊松比相对好办，但屈服强度和硬化模量随温度变化很大。我建议至少取5个温度点的数据：室温、200°C、400°C、600°C、800°C。
4.4 仿真预测与实际裂纹对照
仿真做完了，怎么验证？最好的办法就是跟实际铸件的裂纹图谱对照。
我整理了一个对照表，大家可以参考：
| 裂纹类型 | 仿真特征 | 实际图谱特征 | 对应关系 |
|---|---|---|---|
| 热裂纹 | 应力集中区在凝固前沿 | 沿晶界分布，呈树枝状 | 高度吻合 |
| 冷裂纹 | 高应力区在厚大部位或尖角 | 穿晶断裂，呈直线或弧线 | 位置一致 |
| 应力腐蚀裂纹 | 拉应力+腐蚀介质 | 龟裂状，有腐蚀产物 | 需结合工况 |
| 疲劳裂纹 | 循环应力作用区 | 贝壳纹，有疲劳源 | 需结合载荷 |
我记得有一次，仿真显示某铸件的内腔拐角处应力高达450MPa，而材料的抗拉强度只有400MPa。我建议客户增加圆角半径，从R5改成R15。结果实际生产后，裂纹果然消失了。你看，仿真不是万能的，但没有仿真万万不能。
小技巧：仿真预测裂纹位置时，不要只看应力绝对值。我习惯用"应力安全系数"来判断：安全系数 = 材料抗拉强度 / 仿真应力。如果安全系数小于1.2，就要警惕了。如果小于1.0，那基本可以断定会裂。
4.5 知识体系总览
为了让大家更直观地理解本章的知识结构，我画了一张流程图：
这张图把本章的核心逻辑串起来了。从热应力和冷裂纹的机理出发，到仿真设置，再到实际对照验证，环环相扣。你想想看，如果每一步都做到位了，裂纹问题至少能减少80%。
好了，这一章就到这里。记住：应力与裂纹不是玄学，是科学。用仿真工具去预判、去优化，比事后补救强一百倍。
公众号:蓝海资料掘金营，微信deep3321