4. 材料本构模型与温度相关参数
TSV与RDL结构在热循环、回流焊和功率循环中，会同时经历弹性变形、塑性屈服、蠕变松弛和聚合物粘弹性响应。若把材料简化为线弹性，仿真往往在升降温阶段就低估应力，在保温阶段又完全丢失时间相关变形。因此，本节的核心任务是：为铜、焊料、介质和聚合物四类材料分别选择匹配的本构模型，并给出温度相关参数的获取与拟合路径。
四类材料的本构选择
铜互连与RDL在室温下强度高，但在200～400℃退火或回流过程中会发生明显塑性流动和蠕变。铜的常用模型是弹塑性+蠕变：塑性部分用等向强化描述屈服面扩张，蠕变部分用Norton幂律或双曲正弦形式描述稳态蠕变速率。焊料（SAC305、SnAg、SnPb等）在同系温度高于0.5时蠕变主导，通常采用弹塑性+蠕变或Anand统一粘塑性模型；Anand模型用一个内变量同时描述塑性和蠕变，适合温度跨度大的焊点分析。介质层（SiO₂、SiN、低k介质）以脆性弹性为主，温度相关参数主要是弹性模量和热膨胀系数，一般不需要塑性或蠕变。聚合物（PI、PBO、BCB、环氧模塑料）则表现为粘弹性，常用Prony级数形式的广义Maxwell模型描述应力松弛和蠕变柔量。
温度相关参数的获取方法
弹性模量E(T)和泊松比ν(T)通常来自动态力学分析（DMA）或高温拉伸试验。对于铜，可在-55℃、25℃、100℃、150℃、200℃、250℃六个温度点做拉伸，得到不同温度下的应力-应变曲线，取0.2%偏移屈服强度作为σy(T)。对于焊料，DMA或蠕变试验机在恒温恒载下记录应变-时间曲线，通过不同应力水平下的稳态蠕变速率拟合Norton或双曲正弦参数。对于聚合物，DMA在1Hz、多温度下扫描，得到储能模量E'和损耗模量E''，再通过时温等效原理（WLF方程或Arrhenius方程）平移成主曲线，最后用Prony级数拟合松弛模量。
热膨胀系数CTE(T)一般用热机械分析仪（TMA）测量，升温速率控制在5℃/min，记录尺寸变化。铜的CTE从室温约17 ppm/℃升至250℃约20 ppm/℃；SAC305从21 ppm/℃升至约25 ppm/℃；PI在玻璃化转变前后CTE差异显著，需分段拟合。介质材料CTE较低，SiO₂约0.5 ppm/℃，SiN约2.5 ppm/℃，温度依赖性弱，可取常数或线性近似。
拟合方法与实操要点
以SAC305的Norton蠕变拟合为例，假设在25℃、75℃、125℃三个温度下分别施加10 MPa、20 MPa、30 MPa应力，记录稳态蠕变速率。将双曲正弦模型两边取对数：
ln(ε̇) = ln(A) + n·ln[sinh(ασ)] - Q/(R·T)
步骤：
1. 固定α初值（如0.05 MPa⁻¹），对每个温度做ln(ε̇) vs ln[sinh(ασ)]线性回归，斜率得n，截距得ln(A)-Q/(RT)。
2. 用三个温度的截距对1/T线性回归，斜率得-Q/R，截距得ln(A)。
3. 以R²最大为目标，用遗传算法或fminsearch优化α。
4. 验证：将拟合参数回代，计算预测蠕变速率与实验值对比，误差应小于15%。对于聚合物Prony级数，通常取5～8个松弛时间τᵢ，在时间轴上对数均匀分布（如10⁻⁴、10⁻³、10⁻²、10⁻¹、1、10 s）。用最小二乘拟合各τᵢ对应的权重gᵢ，约束gᵢ≥0且Σgᵢ≤1。拟合后检查主曲线在玻璃化转变区是否平滑，避免过拟合导致外推失真。
参数表与检查清单
| 材料 | 本构模型 | 关键温度相关参数 | 典型获取手段 | 拟合验证指标 |
|---|---|---|---|---|
| 铜 | 弹塑性+蠕变 | E(T)、σy(T)、A、n、Q | 高温拉伸+蠕变试验 | R²≥0.95，蠕变速率误差≤15% |
| SAC305 | Anand/双曲正弦 | E(T)、s₀、Q/R、A、ξ、m、n | 恒温蠕变+DMA | R²≥0.95，多温度外推合理 |
| SiO₂/SiN | 线弹性 | E(T)、ν、CTE(T) | 纳米压痕+TMA | E随T单调下降，CTE线性 |
| PI/EMC | 粘弹性Prony | E(T)、ν(T)、CTE(T)、gᵢ、τᵢ | DMA多温度扫描 | 主曲线平滑，gᵢ非负 |
完成参数拟合后，建议在单单元模型上做一次验证：施加与实验相同的温度和载荷历史，比较仿真与实验的应力-应变响应。若偏差超过20%，应回到数据预处理或模型选择环节排查。只有通过单单元验证的参数，才适合放入TSV/RDL整体模型进行热机电可靠性仿真。