1. 课程导入与工程问题定义
热处理是决定金属零件最终性能的关键工序，而淬火又是其中最容易“出事”的一步。零件在加热、保温、急冷的过程中，内部同时发生两件事：一是温度快速变化引起的热应力，二是奥氏体向马氏体转变引起的组织应力。这两类应力叠加，轻则让零件尺寸超差、圆度失稳，重则直接开裂报废。本课程要解决的，就是如何用数值仿真的方法，把这两种应力“算清楚、看明白、控得住”。
1.1 为什么必须做组织应力耦合仿真
如果只做纯热分析，软件会告诉你零件表面冷得快、心部冷得慢，于是表面受压、心部受拉。但真实淬火远不止如此：当表面温度降到马氏体开始转变温度 Ms 以下时，表层奥氏体开始转变为马氏体，体积膨胀约 4%，这种“相变膨胀”会在表层引入压应力，而心部尚未转变的区域则承受拉应力。更麻烦的是，马氏体转变是随温度连续进行的，转变顺序、转变比例、材料的高温屈服强度都随温度变化，热应力和组织应力在时间轴上此消彼长，甚至发生符号反转。
举一个典型场景：直径 40 mm 的 42CrMo 圆柱试样，奥氏体化 850 ℃ 保温 30 min，随后油淬。若只考虑热应力，表面在冷却初期承受拉应力，最大约 300 MPa；但耦合相变后，表层马氏体膨胀会把表面应力翻转为压应力，峰值可达 −500 MPa 以上，而心部拉应力反而升高。这种“表面压、心部拉”的分布，正是淬火开裂常从心部或过渡区起源的原因。仅靠经验公式或试错，很难判断某一时刻、某一位置的应力是否超过材料当时温度下的抗拉强度。
1.2 两类典型失效：淬火开裂与畸变
淬火开裂通常发生在冷却后期或回火前。裂纹多起源于应力集中处，如键槽、尖角、截面突变区。其判据不是“应力大于屈服”这么简单，而是看最大主应力是否超过该点当前温度、当前相组成下的抗拉强度。仿真中我们关注三个量：最大主应力峰值、出现时刻、以及该时刻该点的马氏体体积分数。若峰值出现在马氏体转变中后期，且位置在心部或过渡区，开裂风险极高。
淬火畸变则表现为尺寸变化、椭圆度、弯曲或翘曲。它由不均匀塑性应变累积造成。仿真中通过提取冷却结束后的等效塑性应变分布和节点位移，可以判断畸变趋势。例如，一个薄壁套筒淬火后内孔可能胀大或缩小，取决于热应力与组织应力谁占主导。若表层先转变马氏体并发生塑性流动，内孔往往胀大；若心部后转变并受表层约束，则可能缩小。这种方向性判断，正是仿真能给出的关键信息。
1.3 仿真在工艺开发中的实际位置
在真实的工艺开发流程里，仿真不是替代试验，而是把试验次数从“几十次”压缩到“几次”。典型做法是：先建立材料数据库，包括 42CrMo 的 CCT 曲线、TTT 曲线、各相的热物性参数（导热系数、比热、密度）、力学参数（弹性模量、泊松比、屈服强度随温度变化）以及相变塑性系数。然后用有限元软件建立零件网格，施加热边界条件（表面换热系数随温度变化），求解温度场；再把温度场结果作为载荷，调用相变动力学模型计算马氏体体积分数；最后把温度场和相变场同时作为初始应变来源，求解应力场。
这里有一个容易被忽略的细节：相变塑性。当材料在相变过程中同时承受应力时，即使应力低于屈服强度，也会产生不可逆变形。这一项如果不加，畸变预测会明显偏小。在 42CrMo 的示例中，相变塑性系数约取 5×10⁻⁵ MPa⁻¹，对薄壁件畸变贡献可达 30% 以上。
1.4 课程项目任务说明
本课程以一个完整的工程案例贯穿始终：42CrMo 圆柱试样与薄壁套筒的淬火过程组织应力耦合仿真。你需要完成以下任务链：
- 材料数据准备：整理 42CrMo 的 CCT 曲线关键点（Ms、Mf、临界冷速），建立各相热物性与力学性能随温度变化的表格。
- 几何与网格：建立圆柱试样（Φ40×80 mm）和薄壁套筒（外径 Φ60，壁厚 8 mm，高 50 mm）的三维模型，采用六面体网格，表面边界层加密。
- 热分析：施加 850 ℃ 奥氏体化保温 30 min，随后油淬，表面换热系数按温度分段设置（示例：850 ℃ 时 1200 W/(m²·K)，400 ℃ 时 2500 W/(m²·K)，200 ℃ 时 800 W/(m²·K)）。
- 相变计算：采用 Koistinen-Marburger 模型计算马氏体体积分数，f_M = 1 − exp[−α(Ms − T)]，其中 α 取 0.011 K⁻¹。
- 应力耦合：将温度场、相变场作为热应变和相变应变输入，考虑相变塑性，求解应力与变形。
- 结果判读：提取最大主应力时程曲线、马氏体分布云图、冷却后位移场，判断开裂风险与畸变趋势，并提出工艺改进建议（如调整冷速、改变装炉方式、增加回火前保温）。