3. 热应力理论基础：CTE失配与热应力产生机理
各位同学，今天我们来聊聊热应力。说实话，做TGV封装，热应力是绕不开的坎。我刚开始接触玻璃通孔时，总觉得玻璃这么硬，应该没问题吧？结果第一次热循环测试，样品直接裂了。嗯，从那以后，我再也不敢小看热应力了。
3.1 热膨胀系数（CTE）失配——问题的根源
热膨胀系数，简称CTE。说白了，就是材料受热时膨胀的"脾气"。每种材料都有自己的CTE值。比如玻璃大约是3-5 ppm/°C，铜大约是17 ppm/°C，硅大约是2.6 ppm/°C。
你想想看，当温度变化时，不同材料膨胀或收缩的程度不一样。它们又紧紧粘在一起，谁也不让谁。这就产生了内应力。我管这叫"材料间的内讧"。
关键概念：CTE失配是指两种或多种材料在温度变化时，因热膨胀系数不同而产生的变形不协调。这是热应力的根本来源。
我在项目中遇到过最典型的例子：铜填充的TGV，升温时铜膨胀得比玻璃快，玻璃就被撑开了。降温时铜收缩得比玻璃快，玻璃又被拉住了。一来二去，界面处就积累了很大的应力。
3.2 热应力产生机理——应力是怎么来的？
热应力的产生，其实分三步走：
- 温度变化：工艺或使用过程中，温度从T1变到T2
- 自由膨胀受阻：每种材料都想按自己的CTE自由膨胀/收缩
- 约束产生应力：但材料之间互相约束，谁也不能自由变形
为什么会这样？我打个比方。你坐地铁，车厢里挤满了人。到站了你想下车，但前后左右都有人挡着。你使劲往外挤，别人也使劲往里挤。这种"互相挤压"就是应力。材料也是一样，温度变了，它们想动，但邻居不让。
我的经验：在TGV中，最危险的界面是铜-玻璃界面。铜的CTE是玻璃的4-5倍。这个失配量，足够让器件在热循环中失效。我曾经见过一个样品，只做了50次热循环，界面就出现了微裂纹。
3.3 热应变与热应力计算公式
好，理论讲完了，我们来点干货。热应变和热应力的计算，其实公式很简单。
热应变公式
热应变只与温度变化和CTE有关：
ε_th = α × ΔT其中：
- ε_th：热应变（无量纲）
- α：热膨胀系数（1/°C）
- ΔT：温度变化量（°C）
举个例子。铜的CTE是17×10⁻⁶ /°C，温度从25°C升到125°C，ΔT=100°C。那么铜的热应变就是：
ε_th = 17×10⁻⁶ × 100 = 0.0017 = 0.17%也就是说，100毫米长的铜柱，会膨胀0.17毫米。看起来不大，但别忘了，玻璃只膨胀0.03%左右。这个差异，就是应力的来源。
热应力公式
当材料被完全约束时，热应力由胡克定律给出：
σ_th = E × ε_th = E × α × ΔT其中：
- σ_th：热应力（Pa）
- E：弹性模量（Pa）
- 其他符号同上
但实际情况更复杂。因为材料不是完全自由的，也不是完全约束的。所以实际应力会小于这个理论值。我一般用这个公式做初步估算，然后靠仿真来精确计算。
注意：这个公式假设材料是线弹性的。如果温度变化很大，材料可能进入塑性区。那时就要用弹塑性模型了。我吃过这个亏，第一次仿真用了线弹性模型，结果应力算出来大得离谱，后来才发现材料早就屈服了。
多层结构的应力计算
对于TGV这种多层结构，情况更复杂。不同层之间互相约束，应力分布不均匀。我常用的方法是：
σ_i = E_i × (α_avg - α_i) × ΔT其中α_avg是各层CTE的加权平均值。这个公式虽然简化了，但能快速判断哪层受拉、哪层受压。
| 材料 | CTE (ppm/°C) | 弹性模量 (GPa) | 热应力方向 |
|---|---|---|---|
| 玻璃 | 3.2 | 70 | 受压（升温时） |
| 铜 | 17.0 | 110 | 受拉（升温时） |
| 硅 | 2.6 | 130 | 受压（升温时） |
你看这个表，升温时铜受拉，玻璃受压。铜的屈服强度大约200-300 MPa，如果热应力超过这个值，铜就会发生塑性变形。我见过不少案例，铜柱在热循环后出现了"挤出"现象，就是塑性变形的结果。
3.4 知识体系总览
为了让大家更直观地理解本章的知识结构，我画了一张图：
这张图把本章的核心逻辑串起来了。从CTE失配出发，到产生机理，再到计算公式，最后落到实际后果。你记住这个链条，热应力的问题就理解了一半。
实用技巧：做仿真前，先用这个公式估算一下应力水平。如果估算值已经接近材料强度，那仿真结果大概率也会出问题。提前心里有数，省得白跑仿真。
好了，热应力的理论基础就讲到这里。记住一句话：CTE失配是根，温度变化是引子，约束条件是放大器。三者缺一，热应力都不会产生。下一节我们聊聊怎么在Ansys里把这些理论变成仿真模型。
公众号:蓝海资料掘金营，微信deep3321