1. 课程导论与工程问题界定

金属材料在冷变形之后，内部储存了大量位错与空位等缺陷，体系处于热力学亚稳态。当温度升高到再结晶温度以上时，原子获得足够迁移能力，变形基体中的晶界开始移动，新的无畸变晶粒形核并长大，最终取代变形组织。这一过程决定了材料的晶粒尺寸、织构、力学性能与后续服役行为，是热加工、退火、焊接、增材制造等工艺链条中的核心环节。本课程围绕"晶界迁移"与"再结晶"两个相互耦合的现象，建立从物理机制到相场建模、再到工程参数调控的完整认知链条。

1.1 工程背景：为什么必须同时关注晶界迁移与再结晶

晶界迁移是再结晶的动力学基础。再结晶形核需要局部晶界弓出或亚晶界合并，而晶粒长大则完全依赖晶界向变形基体一侧迁移。若只研究再结晶形核而不描述晶界迁移速率，就无法预测再结晶完成时间；若只研究晶界迁移而不考虑变形储能分布，就无法解释为什么某些区域优先形核。二者在工程上表现为同一问题的两个侧面。

典型工程问题包括：冷轧低碳钢在罩式退火中如何控制再结晶晶粒尺寸在 15～25 μm；铝合金挤压件如何避免再结晶导致强度下降；增材制造镍基合金如何利用再结晶消除各向异性。这些问题都要求回答：给定温度历史与初始变形组织，晶界在何处、以多快速度迁移，最终形成怎样的晶粒统计分布。

1.2 相场方法的定位：从界面追踪到界面弥散

传统尖锐界面模型需要显式追踪晶界位置，在三维多晶体系中拓扑变化复杂、计算成本高。相场方法引入连续变化的序参量场 η(r,t)，用 η 从 0 到 1 的过渡区表示晶界，晶界位置由 η 的等值面隐式给出。这样，晶界迁移、合并、消失等拓扑事件无需特殊处理，自然被偏微分方程演化所包含。

相场方法在本课程中的定位是"介观尺度桥梁"：向上承接原子尺度晶界能与迁移率的第一性原理或分子动力学数据，向下为有限元热-力耦合与宏观组织统计提供晶粒演化信息。它不替代位错动力学，也不替代宏观传热计算，而是把变形储能、温度场与晶界迁移率组织成一个可求解的演化方程。

1.3 课程目标与项目任务总览

本课程不满足于复述相场方程，而是要求学习者能够独立完成一条从"工艺参数"到"组织预测"的建模链路。具体目标分为三层：

• 概念层：理解晶界迁移驱动力（变形储能差、曲率、晶界能各向异性）与再结晶形核判据（亚晶长大、晶界弓出、应变诱导晶界迁移）的物理含义。
• 方法层：掌握多序参量相场模型的构建方式，包括自由能密度函数、演化方程、初始组织生成、边界条件与数值求解流程。
• 工程层：能够针对给定材料与工艺，标定迁移率与形核参数，运行模拟并提取再结晶分数-时间曲线、晶粒尺寸分布与织构演化，判断模拟结果是否合理。

贯穿课程的项目任务是：以冷轧纯铜（变形量 60%）为对象，建立二维多晶相场模型，模拟 300 ℃ 等温退火过程中的再结晶与晶粒长大，输出再结晶完成时间与平均晶粒尺寸，并与经典 JMAK 理论及文献实验趋势进行对比。项目将分阶段嵌入各章节，最终形成一份可复现的建模报告。

1.4 关键参数与建模约定

为使后续章节的方程与图示保持一致，此处先约定符号与参照系。空间坐标采用右手系，晶界迁移方向定义为从高储能区指向低储能区，即驱动力 F = -δE/δη。温度 T 为绝对温度，时间 t 单位为秒。相场变量 η_i 表示第 i 个晶粒的取向，在晶粒内部 η_i ≈ 1，在晶界处连续过渡到 0。

示例参数（合理假设，非实测值）：纯铜晶界能 γ = 0.6 J/m²，晶界迁移率 M 在 300 ℃ 下取 1.0×10⁻¹⁴ m⁴/(J·s)，变形储能 E_stored = 5×10⁵ J/m³，网格间距 Δx = 0.5 μm，时间步长 Δt = 0.1 s。这些数值用于演示量级关系，实际建模需根据材料与温度重新标定。

1.5 从工程问题到相场模型的思维路径

面对一个再结晶问题，第一步不是写方程，而是界定"什么在变、什么驱动它变、什么阻碍它变"。以冷轧铜退火为例：变形储能是驱动力，晶界能与曲率是阻力，温度通过迁移率放大动力学。相场模型把这三者分别放入自由能密度、梯度项和演化系数中。若模拟结果显示再结晶分数在 10 s 内从 0 跳到 1，而实验趋势是数百秒，则说明迁移率或形核密度标定偏高，需要回到参数层调整，而不是修改方程结构。