1. 课程导入与多物理场认知

在工程材料与结构失效分析中，有一类问题很难用单一学科解释：一块高温服役的合金叶片，表面氧化层不断增厚，内部却出现裂纹；一块锂离子电池电极，在多次充放电后容量骤降，极片发生褶皱；一块焊接接头，在腐蚀介质中放置数月后，焊缝附近出现延迟开裂。这些现象背后往往同时涉及原子扩散、应力演化与温度、电场、化学势等多个物理场的相互作用。本课程要建立的，正是把扩散、应力与多物理场耦合放在同一分析框架下理解与计算的工程能力。

1.1 三个基本概念：扩散、应力、多物理场耦合

扩散是原子或分子在浓度梯度、化学势梯度或电势梯度驱动下发生的定向迁移。工程中最常用的是Fick第一定律与第二定律：

J = -D · ∂C/∂x          （Fick第一定律，稳态）
∂C/∂t = ∂/∂x ( D · ∂C/∂x )   （Fick第二定律，非稳态）

其中J为扩散通量，D为扩散系数，C为浓度。D对温度极其敏感，通常服从Arrhenius关系 D = D₀·exp(−Q/RT)。以钢中碳扩散为例，在727℃附近D约为10⁻¹¹ m²/s量级，而降到200℃时可低至10⁻²⁰ m²/s量级，相差近九个数量级。这意味着温度场哪怕只变化几十摄氏度，扩散进程就可能发生数量级改变。

应力是材料内部单位面积上的内力。扩散与应力之间存在双向耦合：一方面，原子半径差异、浓度梯度会产生化学膨胀应变，例如锂嵌入石墨时体积膨胀可达约10%，碳在奥氏体中富集也会引起晶格畸变；另一方面，应力又会反过来改变原子迁移的驱动力，即扩散的驱动力不只是浓度梯度，还包括应力梯度对应的势能项。因此扩散方程中常出现应力耦合项，应力平衡方程中又包含浓度引起的本征应变。

多物理场耦合指两个或以上物理场通过共享变量、材料参数或边界条件相互影响。材料扩散应力问题中，典型耦合链为：温度场→扩散系数与浓度场→化学应变→应力场→应力梯度→扩散通量→浓度场再分布。这是一个闭环，不能按单向顺序简单求解。

1.2 为什么必须做耦合分析

如果忽略耦合，只做单向计算，会得到什么？以锂电极为例，若只算扩散浓度场，可以知道锂浓度分布；若只算应力，可以知道给定膨胀下的应力。但真实问题是：锂浓度分布决定膨胀应变分布，膨胀应变分布决定应力分布，应力分布又通过应力梯度项改变锂的扩散路径，最终浓度分布与初始假设不同。单向计算会低估峰值应力，甚至给出错误的裂纹位置。工程上，这类误差可能导致寿命预测偏差数倍。

耦合分析的核心难点在于时间尺度差异与空间尺度差异。扩散特征时间 t_d ~ L²/D，应力波传播时间 t_s ~ L/c。以L=1 mm、D=10⁻¹⁴ m²/s为例，t_d约为10⁶ s量级；而应力波传播时间约为10⁻⁷ s量级，相差十三个数量级。因此工程中通常采用准静态力学+瞬态扩散的弱耦合策略：在每个扩散时间步内求解力学平衡，再把应力结果代入下一步扩散计算。

1.3 典型失效案例与工程背景

案例一：高温合金氧化层剥落。NiCrAlY涂层在1100℃氧化时，Al向外扩散形成Al₂O₃层，同时涂层内Al浓度下降。氧化层与基体热膨胀系数不匹配，冷却时产生压应力，当氧化层厚度超过约5～8 μm时，界面附近应力超过结合强度，发生剥落。剥落后新鲜表面暴露，氧化加速，形成循环失效。

案例二：锂离子电池硅负极颗粒断裂。硅嵌锂时体积膨胀约300%，颗粒内部产生巨大应力。若充电倍率过高，锂浓度梯度大，颗粒表面与中心应变差大，表面应力超过硅的断裂强度（约1～2 GPa量级），颗粒开裂，SEI膜反复生长，锂库存损失，容量快速衰减。工程上通过纳米化、碳包覆、预锂化等手段降低浓度梯度和应力。

案例三：焊接接头氢致开裂。焊接时高温区溶解大量氢，冷却后氢在应力梯度驱动下向焊缝热影响区扩散富集。当局部氢浓度超过临界值且应力强度因子超过门槛值时，发生延迟开裂。这类失效可能在焊后数小时至数天才出现，具有隐蔽性。

1.4 本课程路线图

本课程共30章，按“基础→耦合→专题→综合”四段推进。前段建立扩散方程、应力平衡、热传导与电场的基本描述；中段重点讲耦合项的来源、弱耦合与强耦合求解策略、有限元实现要点；后段进入氧化、电池、焊接、互连等专题，每类问题都给出可复现的参数设置、边界条件与结果判读方法；最后以综合案例串联多场耦合分析全流程。