4、材料属性定义：硅、铜、二氧化硅、焊料等材料的热、电、力学参数、温度依赖性设置
各位工程师，咱们今天聊聊材料属性定义。说实话，这是多物理场仿真里最容易被忽视，却又最容易翻车的一环。我见过太多仿真结果漂漂亮亮，一到实测就崩盘，最后查来查去，问题出在材料参数上。
你想想看，TSV 结构里就那么几种材料：硅衬底、铜导体、二氧化硅绝缘层、还有焊料。但每种材料在不同温度、不同工艺条件下，表现完全不一样。今天我就把我在项目中踩过的坑、积累的经验，一次性说清楚。
核心观点：材料属性不是常数，是变量。温度、频率、应力状态都会改变它。忽略这一点，仿真就是自欺欺人。
4.1 硅（Si）—— 衬底材料的双面性
硅是半导体，不是纯导体也不是纯绝缘体。它的电导率随温度变化非常明显。我个人习惯用以下模型来描述硅的电阻率：
ρ(T) = ρ₀ * [1 + α*(T - T₀)]其中 ρ₀ 是室温下的电阻率，α 是温度系数。对于轻掺杂硅，α 大约是正值，温度越高电阻越大。但重掺杂硅就不一样了，温度系数可能变成负的。
热导率方面，硅其实是个不错的导热体，室温下约 130 W/(m·K)。但温度升高到 200°C 以上时，热导率会下降到 80 左右。我在做功率器件仿真时，就因为这个变化没考虑进去，导致结温预测偏差了 15°C 以上。
| 参数 | 室温值 (25°C) | 高温值 (200°C) | 温度依赖性 |
|---|---|---|---|
| 热导率 (W/m·K) | 130 | ~80 | ∝ T⁻¹·³ |
| 热膨胀系数 (ppm/K) | 2.6 | 3.5 | ∝ T⁰·⁵ |
| 杨氏模量 (GPa) | 130 | 125 | 弱依赖 |
| 电阻率 (Ω·cm) | 10 (轻掺杂) | ~5 | ∝ T⁻¹·⁵ |
我的经验：硅的各向异性不能忽略。不同晶向的杨氏模量能差 30% 以上。如果你做的是高精度应力仿真，记得用单晶硅的弹性矩阵，别用各向同性近似。
4.2 铜（Cu）—— 导电导热都靠它
铜是 TSV 的灵魂材料。它的电导率和热导率都很高，但问题在于——它和硅的热膨胀系数不匹配。铜的 CTE 大约是 17 ppm/K，硅只有 2.6 ppm/K。温度一变化，应力就来了。
电导率温度模型：
σ(T) = σ₀ / [1 + α_cu*(T - 20°C)]α_cu 大约是 0.0039 /°C。什么意思呢？温度每升高 100°C，铜的电阻率增加约 39%。我在做电流密度仿真时，如果不考虑这个，算出来的焦耳热会偏小很多。
力学参数要注意：铜在高温下会软化。室温下杨氏模量约 110 GPa，到 300°C 时只剩 80 GPa 左右。而且铜有塑性变形行为，应力超过屈服强度后会发生蠕变。我曾经在热循环仿真中忽略了蠕变，结果预测的寿命比实测长了 3 倍。
避坑指南：铜的电迁移效应和温度强相关。电流密度超过 1e5 A/cm² 时，铜原子会沿着电子流方向迁移，形成空洞。这个效应在 150°C 以上会急剧加速。我建议在仿真中至少设置两个温度点来评估电迁移风险。
4.3 二氧化硅（SiO₂）—— 绝缘层的隐形杀手
二氧化硅是 TSV 侧壁的绝缘层，厚度通常只有 0.5-2 μm。别小看这层薄薄的氧化物，它的热导率极低，只有 1.4 W/(m·K)，是硅的百分之一。这意味着什么？
热量从铜柱传到硅衬底时，二氧化硅层会成为主要热阻。我做过一个对比仿真：不考虑 SiO₂ 热阻时，TSV 结温是 85°C；考虑之后，结温直接飙到 105°C。差了 20°C！
介电常数方面：SiO₂ 的介电常数约 3.9，频率依赖性很小。但温度升高到 200°C 以上时，介电常数会略微增加，大约到 4.2 左右。如果你做的是高频信号完整性仿真，这个变化不能忽略。
| 参数 | 典型值 | 温度影响 | 频率影响 |
|---|---|---|---|
| 热导率 (W/m·K) | 1.4 | 随温度升高略有增加 | 无 |
| 介电常数 | 3.9 | 200°C 时约 4.2 | 弱 |
| 击穿场强 (MV/cm) | 10 | 高温下降至 6-7 | 无 |
| 热膨胀系数 (ppm/K) | 0.5 | 基本恒定 | 无 |
我的建议：在热-电耦合仿真中，一定要把 SiO₂ 层单独建模，不要用等效热阻代替。因为它的厚度很薄，但温度梯度很大，等效模型会丢失很多细节。
4.4 焊料（Solder）—— 最不稳定的环节
焊料是 TSV 互连中温度依赖性最强的材料。常用的 SAC305（Sn-3.0Ag-0.5Cu）焊料，熔点约 217°C。在服役温度范围内（-55°C 到 125°C），它的力学性能变化非常大。
蠕变行为：焊料在高温下会发生明显的蠕变。我用的是 Anand 粘塑性模型来描述：
ε̇ = A * exp(-Q/RT) * [sinh(ξσ/s)]^(1/m)这个模型有 9 个参数，每个参数都有物理意义。A 是常数，Q 是激活能，R 是气体常数，T 是绝对温度。说白了，温度越高，蠕变速率越快。
电导率方面：焊料的电导率大约是铜的 1/10，约 7.5e6 S/m。温度系数和铜类似，约 0.003 /°C。但要注意，焊料在服役过程中会发生金属间化合物（IMC）生长，这会增加电阻。
| 参数 | SAC305 (25°C) | SAC305 (125°C) | 变化趋势 |
|---|---|---|---|
| 杨氏模量 (GPa) | 50 | 25 | 下降 50% |
| 屈服强度 (MPa) | 35 | 10 | 下降 70% |
| 蠕变激活能 (eV) | 0.8 | 0.8 | 恒定 |
| 热导率 (W/m·K) | 58 | 55 | 弱下降 |
我曾经踩过的坑：焊料的疲劳寿命和温度循环范围直接相关。ΔT 从 100°C 增加到 150°C，寿命可能下降一个数量级。做可靠性仿真时，一定要用实测的 Coffin-Manson 参数，别用默认值。
4.5 温度依赖性设置实操
在 ANSYS 或 COMSOL 中设置温度依赖参数时，我一般用两种方法：
- 多项式拟合：对于硅的热导率，用二次多项式拟合效果不错。比如：k(T) = 130 - 0.3*(T-25) + 0.0005*(T-25)²
- 分段线性插值：对于焊料这种非线性很强的材料，用 5-10 个温度点的实测值做插值更准确。
我个人更推荐第二种方法。为什么？因为多项式在高温区容易发散，而插值法只要数据点够密，精度就有保障。
关键提醒：所有温度依赖参数都要在仿真前做一次验证。拿一个简单结构（比如单根 TSV）跑一下，看结果是否和文献或实测数据吻合。这一步花不了多少时间，但能避免后面返工。
4.6 知识体系总览
下面这张图总结了材料属性定义的核心逻辑。你看，从四种基础材料出发，每个材料都有热、电、力学三个维度的参数，而温度依赖性贯穿始终。做仿真时，脑子里要有这张图。
好了，材料属性这块就聊到这儿。记住一句话：参数不准，仿真白费。下次做仿真前，先花半小时把材料参数理清楚，你会发现后面省下的时间远不止半小时。
公众号:蓝海资料掘金营，微信deep3321