4. 材料本构与界面建模
在CPO光电封装的热机械协同仿真中，材料本构与界面建模是决定仿真可信度的第一道关口。芯片、基板、光纤阵列、透镜、焊点、underfill与封装外壳由多种材料组成，它们的弹性模量、热膨胀系数、玻璃化转变温度、蠕变与黏弹性行为各不相同。若材料参数取值失真，或界面被简化为理想绑定，后续的温度场、应力场与光耦合位移结果都会出现系统性偏差。本节围绕“参数从哪里来、本构如何选、界面层如何建”三个问题展开。
一、材料热力学参数获取的三条路径
材料参数获取通常有三条路径：查手册与供应商数据、实验测试、以及基于分子动力学或微观力学的参数标定。工程仿真中优先使用前两者，第三条用于缺乏数据的新材料或界面层等效参数。
硅作为芯片与光波导主体，单晶硅在<100>、<110>、<111>晶向的弹性模量各向异性明显，室温下杨氏模量约130～188 GPa，热膨胀系数约2.6×10⁻⁶/K。硅在CPO中通常按各向同性或正交各向异性线弹性处理，温度不超过400℃时可不考虑塑性。玻璃（如硼硅玻璃、石英）用于透镜、光纤与玻璃基板，其模量约60～75 GPa，CTE约0.5～5×10⁻⁶/K，需注意玻璃在Tg附近黏弹性显著，若工艺温度接近Tg，应改用黏弹性本构。
聚合物（如光学胶、模塑材料、PI）与underfill是典型的黏弹性材料。其参数包括瞬时模量、平衡模量、松弛时间谱、CTE（Tg前后两段）、固化收缩率与吸湿膨胀系数。焊料（SAC305、SnAgCu、AuSn）在室温至熔点之间经历蠕变与塑性，常用Anand模型、Garofalo双曲正弦蠕变模型或Darveaux模型描述。下表给出典型示例参数，仅用于教学演示，实际项目应以供应商数据或实测为准。
| 材料 | 弹性模量 E | CTE (×10⁻⁶/K) | 泊松比 ν | 本构建议 |
|---|---|---|---|---|
| 单晶硅 | 130～188 GPa | 2.6 | 0.22～0.28 | 各向异性线弹性 |
| 硼硅玻璃 | 64 GPa | 3.3 | 0.20 | 线弹性/黏弹性 |
| Underfill | 8～10 GPa (Tg以下) | 25 (Tg以下)/90 (Tg以上) | 0.30 | 黏弹性+固化收缩 |
| SAC305焊料 | 40 GPa (25℃) | 21 | 0.35 | Anand蠕变+塑性 |
| AuSn焊料 | 68 GPa | 16 | 0.35 | 弹塑性+蠕变 |
参数获取时需注意：CTE应给出参考温度Tref，仿真中热应变按α(T−Tref)计算；黏弹性材料的松弛时间谱需覆盖工艺与服役温度区间；焊料蠕变参数需与温度、应力水平匹配。若供应商只给室温模量，应通过DMA或纳米压痕补充温度扫描数据。
二、本构模型的选择逻辑
本构选择取决于材料在工艺与服役中的实际响应。硅、玻璃在低温段用线弹性即可；underfill在固化冷却与温度循环中会经历应力松弛，应使用Prony级数形式的广义Maxwell模型或时间-温度等效原理（WLF/Arrhenius移位因子）描述。焊料在温度循环中同时发生塑性变形与蠕变，Anand模型用一个流动方程统一描述黏塑性，参数包括A、Q/R、ξ、m、ŝ、n、h₀、a、s₀，适合无铅焊料。
以Anand模型为例，其流动方程为：
ε̇p = A [sinh(ξ·q/s)]^(1/m) · exp(−Q/(R·T))
ṡ = [h₀·(1 − s/s*)^a · sign(1 − s/s*)] · ε̇p
s* = ŝ · [ε̇p/A · exp(Q/(R·T))]^n其中q为等效应力，s为内变量，s*为饱和值。参数标定通常用等温拉伸/剪切蠕变试验，拟合不同温度与应力下的稳态蠕变速率。若缺少试验条件，可先用文献值，但必须做参数敏感性分析，确认关键输出（如焊点等效塑性应变、界面剪应力）对参数的敏感程度。
三、界面层建模方法
CPO中关键界面包括：芯片-焊料、焊料-基板、underfill-芯片、underfill-基板、光纤-玻璃、玻璃-聚合物。界面建模有三种常用方法：
- 共节点绑定：假设界面完美连接，网格共节点。计算量小，适合界面强度远高于体材料、且不关心界面失效的场景。
- 内聚力单元（Cohesive Zone Model, CZM）：在界面插入零厚度或薄层单元，用双线性或指数型牵引-分离法则描述界面损伤起始与演化。参数包括界面强度σmax、断裂能Gc、初始刚度Kn。适合评估分层、开裂风险。
- 等效薄层实体：用一层薄实体单元（厚度1～5 μm）代表界面过渡区，赋予等效模量与CTE。适合界面存在IMC、扩散层或粗糙度影响的场景。
界面层建模的关键是厚度与刚度的匹配。若界面层厚度远小于相邻体网格尺寸，直接建薄层会导致网格畸变与刚度病态，此时应使用CZM或绑定约束。CZM参数可通过双悬臂梁（DCB）试验、四点弯曲试验或微划痕试验获取；若无法试验，可按界面强度约为较软材料屈服强度的0.5～1倍、断裂能按经验值估算，并做区间分析。
四、从参数到模型的落地流程
实际项目中，材料与界面建模可按以下顺序推进：先建立材料卡片，记录每个参数的来源、温度范围与不确定度；再根据工艺温度曲线与服役温度范围选择本构；然后对关键界面做CZM或薄层建模；最后通过网格收敛性与参数敏感性分析确认模型可信度。下面用一张图概括材料参数、本构选择与界面建模之间的对应关系。
图中左侧为材料对象，中间为需要获取的关键参数，右侧为对应的本构与界面建模方法。实际建模时，应先确认每个材料在工艺与服役温度区间内的主导变形机制，再选择本构；界面则根据是否关心分层失效来决定绑定、CZM或薄层实体。参数来源、温度范围与不确定度应记录在材料卡片中，便于后续追溯与敏感性分析。