3. 热机械协同仿真基础
CPO（共封装光学）器件的可靠性问题，本质上很少是单一物理场造成的。激光器发热导致基板翘曲，翘曲又改变光纤耦合对准，耦合损耗上升进一步加剧局部发热——热与力在这里是互相咬合的两条链。要预测这种咬合，就必须把热传导方程与力学平衡方程放在同一个求解框架里，这就是热机械协同仿真的起点。
一、热-力耦合的控制方程
热机械耦合的数学描述由两组方程构成。第一组是瞬态热传导方程：
ρ·c_p·∂T/∂t = ∇·(k∇T) + Q其中 ρ 为密度，c_p 为比热容，k 为导热系数张量，Q 为体热源密度（如激光器有源区的焦耳热与光吸收热）。第二组是力学平衡方程（准静态，忽略惯性项）：
∇·σ + f = 0两者通过热应变耦合：σ = C : (ε − ε_th)，其中热应变 ε_th = α·(T − T_ref)，α 为热膨胀系数（CTE），T_ref 为无应力参考温度。这个式子说明：温度场是力学场的"载荷源"，而力学场一般不会显著反过来改变温度场（除接触热阻随压力变化等弱耦合情形），因此CPO热机械问题通常按单向耦合处理——先算温度，再把温度作为体载荷传给力学求解。
需要特别强调参考温度 T_ref 的约定：它必须与工艺装配温度一致。若芯片在 250 ℃ 回流焊后冷却到 25 ℃，则 T_ref 应取 250 ℃，室温下的残余应力才不会被漏算。全文与图示均采用"以装配温度为无应力参考态、升温为正、压应变为负"的约定。
二、有限元法基础：从弱形式到刚度矩阵
上述偏微分方程在任意几何上无解析解，有限元法（FEM）的思路是：把求解域离散为单元，用形函数 N 近似场变量 T ≈ ΣN_i·T_i、u ≈ ΣN_i·u_i，再通过加权残差（Galerkin 法）把强形式转化为弱形式，最终得到代数方程组。
热问题离散后得到：C·Ṫ + K_th·T = Q，其中 C 为热容矩阵，K_th 为热导矩阵。力学问题离散后得到：K_u·u = F_th + F_ext，其中 K_u 为刚度矩阵，F_th = ∫Bᵀ·C·ε_th dV 是由温度场产生的等效热载荷。求解流程是：先解热方程得到节点温度 T，再逐单元计算 ε_th，组装 F_th，最后解力学方程得到位移 u，进而由 ε = B·u 得到应变、由 σ = C:(ε − ε_th) 得到应力。
对CPO封装而言，网格策略直接决定结果可信度。焊点、微凸点、光纤阵列单元（FAU）胶层这些尺寸跨度大的区域，宜采用二阶四面体或六面体单元，并在界面处加密；而硅基板、PCB 等大块体可用较粗网格。一个实用的经验是：在焊点最细颈处至少保证 3 层单元，否则塑性应变与疲劳寿命会被严重低估。
三、多物理场耦合策略
耦合策略按信息传递方向可分为三类。顺序（弱）耦合：热→力单向传递，计算量小，适合热源稳定、力学反馈不影响温度的场合，是CPO封装应力分析的主力方法。强耦合：把温度与位移放在同一方程组中同时求解，能捕捉接触热阻随压力变化、大变形改变散热路径等效应，但矩阵规模大、收敛困难。迭代耦合：交替求解两个物理场，直到界面量（温度、热流、位移）收敛，介于两者之间。
在CPO中还有一类容易被忽视的耦合：热-光-力。温度不仅产生应力，还通过热光效应改变硅波导折射率（dn/dT ≈ 1.8×10⁻⁴ /K），应力又通过光弹效应改变折射率。因此严格的光学对准分析需要把温度场、位移场作为输入传给光学求解器，计算模式有效折射率偏移与耦合损耗。工程上常先做热-力，再把位移结果映射到光学模型做一次灵敏度评估。
四、工具链概览与协同工作流
CPO 热机械协同仿真的工具链通常由三部分构成。几何与网格前处理：SpaceClaim、ANSYS Discovery 或 COMSOL 的 CAD 导入模块，负责清理封装模型、抽取流体域（若含风冷/液冷）、生成共形网格。多物理场求解器：ANSYS Mechanical + Icepak（热-力顺序耦合）、COMSOL Multiphysics（内置热膨胀与热光模块）、Abaqus（擅长焊点非线性与疲劳）。后处理与光学映射：把位移场导出为 CSV 或 VTK，再导入 Lumerical FDTD/MODE 或 COMSOL Wave Optics 计算模式偏移与耦合效率。
一个可复现的最小工作流如下（以 ANSYS 为例）：
- 在 Icepak 中建立封装热模型，设定芯片功耗 0.5 W、TEC 功耗 0.3 W、环境 25 ℃、对流系数 10 W/(m²·K)，求解稳态温度场。
- 将温度结果作为体载荷导入 Mechanical，设置 T_ref = 250 ℃（回流焊峰值），材料 CTE：硅 2.6 ppm/K、铜 17 ppm/K、焊料 21 ppm/K、Underfill 30 ppm/K（玻璃态）。
- 施加位移约束：PCB 底面四角简支，避免刚体位移。
- 求解后提取 FAU 与硅光芯片界面的相对位移 Δx、Δy、Δz，以及焊点等效塑性应变。
- 把 Δx、Δy、Δz 映射到光学模型，计算耦合损耗变化；若损耗超过 0.5 dB，则回到步骤 1 调整 TEC 控制策略或 Underfill 材料。