温度梯度对结构的影响 — 热应力、变形机制与FEM耦合分析
理论与物理
为何温度梯度会产生应力？
老师，存在温度梯度会对结构产生什么影响？温度升高不就只是膨胀而已吗？
问得好。确实，如果均匀加热且能自由膨胀，应力就是零。但是存在温度梯度时，情况就完全不同了。高温部分想要大幅膨胀，而低温部分“还不想膨胀那么多”——这种内部约束会即使没有外部约束也会产生应力。
诶，即使没有外部压制也会产生应力吗？这和直觉不太一样呢…
嗯，这就是温度梯度问题的核心。我们用一个具体例子来思考。想象一下涡轮增压器的涡轮壳体。排气侧的内表面大约800°C，被冷却的外表面大约500°C。内表面想要大幅膨胀，但外表面却把它拉住。结果，内表面产生压缩应力，外表面产生拉伸应力。这种应力分布在发动机每次启动和停止时重复发生，成为热疲劳裂纹的原因。
原来如此…内表面和外表面在互相拉扯对吧。那么温度差越大，应力也越大吗？
没错。应力与温度差 $\Delta T$ 成比例增加。此外，还与杨氏模量 $E$ 和线膨胀系数 $\alpha$ 成比例。所以铝（$\alpha \approx 23 \times 10^{-6}$/K）产生的热应变大约是钢（$\alpha \approx 12 \times 10^{-6}$/K）的2倍，但杨氏模量是钢的1/3，因此实际热应力的大小取决于材料的组合。
基本公式：约束体的热应力
定量上可以用什么公式表示呢？
首先考虑最基本的情况——完全约束的均匀材料的平面应力。对于温度变化 $\Delta T$，自由热应变是 $\varepsilon_{th} = \alpha \Delta T$，但由于完全约束无法膨胀，这部分直接作为弹性应变产生应力：
这里 $E$ 是杨氏模量，$\alpha$ 是线膨胀系数，$\nu$ 是泊松比。除以 $(1-\nu)$ 是双轴约束（平面应力面内两个方向约束）的效果。如果是单轴约束，分母就是1；三轴约束（体积变化受抑制）则是 $(1-2\nu)$。负号表示“温度上升→压缩应力”。
我想代入具体数字试试。比如钢在温度差100°C时呢？
以钢为例，设 $E = 200$ GPa，$\alpha = 12 \times 10^{-6}$ /K，$\nu = 0.3$，计算如下：
$$\sigma = \frac{200 \times 10^3 \times 12 \times 10^{-6} \times 100}{1 - 0.3} = \frac{240}{0.7} \approx 343 \;\text{MPa}$$
仅仅100°C的温度差就达到343 MPa——超过了软钢的屈服应力（约250 MPa）。这下明白温度梯度的影响有多大了吧？
100°C就超过屈服极限！这在设计中可不能忽视啊…
板厚方向梯度：膜力与弯矩
像刚才涡轮增压器那样，板的正面和背面温度不同时，情况会更复杂吗？
是的。对于厚度为 $h$ 的板，存在厚度方向的温度分布 $T(y)$ 时，温度效应可以分解为膜力分量和弯矩分量。设板的中立面为 $y=0$：
直观地说，$N_{th}$ 是“板整体均匀伸长的效果”，$M_{th}$ 是“板发生翘曲的效果”。如果温度分布关于中立面对称，则 $M_{th} = 0$（无翘曲）；如果反对称，则 $N_{th} = 0$（无伸长，只有弯曲）。
例如，如果只加热单面，就会同时产生膜力和弯矩，对吧？
没错。例如，被阳光直射的桥梁上表面为50°C，背阴的下表面为30°C，假设为线性梯度 $T(y) = 40 + 20y/h$，则 $N_{th}$ 和 $M_{th}$ 都会产生。这种热梯度弯曲会使长大桥梁在白天呈弓形变形数毫米——这在桥梁设计中作为温度荷载已在规范中予以考虑。
双金属条的曲率公式
恒温器中使用的双金属也是温度梯度的应用吗？
你注意到了关键点。双金属是将两种金属粘合在一起的结构，利用线膨胀系数的差异将温度变化转换为曲率变化的装置。Timoshenko（1925年）推导出了经典的曲率公式：
这里 $m = t_1/t_2$（两层厚度比），$n = E_1/E_2$（杨氏模量比），$h = t_1 + t_2$（总板厚）。当 $m = 1$（等厚）、$n = 1$（等刚度）时曲率最大——也就是说，厚度相同、硬度相同，只有线膨胀系数不同的组合弯曲得最厉害。
原来如此！这个公式似乎也能用于预测电子电路板的翘曲呢。
正是如此。在电子封装中，将硅芯片（$\alpha \approx 2.6 \times 10^{-6}$/K）用焊料连接到有机基板（$\alpha \approx 15 \times 10^{-6}$/K）上时，从回流炉冷却下来会产生与双金属相同原理的翘曲。这直接关系到焊点接合部的可靠性。
面内温度梯度与平面应力
不仅有板厚方向的温度梯度，也有面内方向的温度梯度吧？
是的。对于存在面内温度梯度 $T(x, y)$ 的平板，热弹性控制方程可以用应力函数 $\phi$ 写成如下形式：
右边的 $\nabla^2 T$ 是关键。如果温度分布是调和函数（拉普拉斯算子 = 0）——也就是说，在稳态热传导下没有内部热源，温度分布平滑时——右边为零，仅凭面内温度梯度不会产生内部应力。不过，这仅限于边界自由的情况。如果边界受到约束，当然会产生应力。
温度分布的形状本身会影响应力呢。稳态和非稳态结果会不同，是这个意思吗？
没错。在非稳态（瞬态）温度场中，$\nabla^2 T = (\rho c_p / k) \, \partial T / \partial t \neq 0$，因此即使边界自由也会产生内部应力。发动机启动时或喷气发动机循环等伴随急剧温度变化的瞬态最为危险——因为此时会产生最大的温度梯度。
各项的物理意义
- $\sigma = E\alpha\Delta T/(1-\nu)$：双轴约束平板的热应力。$E\alpha\Delta T$ 是单轴热应力，用 $(1-\nu)$ 修正面内两个方向的泊松效应。直接关系到钢轨的屈曲、受约束螺栓接合部的评估。
- $N_{th}$（热膜力）：相当于板厚方向温度分布的“平均值”分量。使板整体均匀伸缩的效果。在受约束的板中会产生面内应力。
- $M_{th}$（热弯矩）：相当于温度分布的“非对称分量”。引起板的翘曲/弯曲。太阳热引起的桥梁变形、PCB基板翘曲的主要原因。
- Timoshenko的双金属曲率 $\kappa$：预测异种材料接合体在温度变化时的曲率。不仅应用于恒温器、双金属开关，还广泛应用于电子封装的翘曲、MEMS执行器的设计。
假设条件与适用范围
- 假设为线弹性体。高温下发生屈服/蠕变时，需要非线性材料模型
- 假设材料特性与温度无关（$E$, $\alpha$, $\nu$ 恒定）。实际上高温下变化很大，因此500°C以上必须使用温度相关数据
- 小变形理论。双金属大曲率变形时需要考虑几何非线性
- 假设单向耦合（结构变形不影响温度场）。存在摩擦发热或接触变化时，需要双向耦合
- 双金属公式假设完全接合（界面无滑移）。未考虑粘接剂层的剪切变形
量纲分析与单位制
| 变量 | SI单位 | 典型值・注意事项 |
|---|---|---|
| 杨氏模量 $E$ | Pa (N/m²) | 钢: 200 GPa, Al: 70 GPa, Cu: 120 GPa |
| 线膨胀系数 $\alpha$ | 1/K | 钢: 12×10⁻⁶, Al: 23×10⁻⁶, Si: 2.6×10⁻⁶ |
| 泊松比 $\nu$ | 无量纲 | 金属: 0.25〜0.35, 橡胶: ≈0.5 |
| 温度差 $\Delta T$ | K（与℃等同） | 注意单位制混淆（°F→K转换遗忘频发） |
| 曲率 $\kappa$ | 1/m | 曲率半径 $R = 1/\kappa$ |
弯曲新干线轨道的“无形之力”
夏日的午后，新干线轨道上表面因日照可达60°C以上。另一方面，与道砟（碎石）接触的下表面则停留在40°C左右。这约20°C的温度梯度会在钢轨上产生弯矩，使其向上弯曲。在长钢轨区间，这种效应会累积，引发称为“轨道不平顺”的路基微小变形。新干线的安全运行要求以0.5mm为单位管理轨道，因此这种温度梯度效应是养护计划中的重要参数。类似的现象也发生在混凝土桥上，白天和夜间桥面高度会变化数毫米。“温度不同，结构必动”——这个简单的原理支配着基础设施设计。
数值解法与实现
顺序耦合分析的流程
用FEM分析温度梯度的影响需要什么步骤？
标准方法是使用顺序耦合（单向耦合）。步骤分为两步：
- Step 1 — 热传导分析：设定边界条件（对流、辐射、接触热阻等），求解温度场 $T(\mathbf{x}, t)$
- Step 2 — 结构分析：将Step 1得到的温度场作为“温度荷载”映射到结构模型，计算位移和应力
温度→结构是单向传递，假设结构变形不影响温度场，所以是“单向”。这在许多工业应用中足以提供足够的精度。
两个分析是分开做的啊。网格使用相同的吗？
原则上使用相同网格，温度→位移的映射会更准确。但在实际工作中，用于热分析的网格（流体边界层细化）和用于结构分析的网格（应力集中部位细化）常常不同。此时，需要通过插值映射（最近邻、投影法、RBF插值等）来传递温度场。在ANSYS Workbench中，如果网格相同，只需设置一个标志即可联动；即使网格不同，也能自动进行映射。
热传导方程的离散化
热传导分析的方程是什么形式？
傅里叶热传导方程为：
$\rho$ 是密度，$c_p$ 是比热容，$k$ 是热导率，$\dot{q}$ 是内部发热（焦耳热等）。用FEM离散化后：
$$[C]\{\dot{T}\} + [K_T]\{T\} = \{Q\}$$
$[C]$ 是热容矩阵，$[K_T]$ 是热传导矩阵，$\{Q\}$ 是热流向量。稳态分析时设 $\{\dot{T}\} = 0$，只需求解 $[K_T]\{T\} = \{Q\}$ 即可。
热荷载向量的构成
Step 1得到温度后，怎么传递给结构分析呢？温度作为“荷载”有点难以想象。
问得好。结构的FEM方程是 $[K]\{u\} = \{F\}$，而温度荷载是包含在右边的力向量 $\{F\}$ 中的。单元级别的热荷载向量为：
$$\{F_{th}\}_e = \int_{\Omega_e} [B]^T [D] \{\varepsilon_{th}\} \, d\Omega$$
这里 $[B]$ 是应变-位移矩阵，$[D]$ 是弹性矩阵，$\{\varepsilon_{th}\} = \alpha \Delta T \{1, 1, 1, 0, 0, 0\}^T$（3D）是热应变向量。也就是说，将根据各节点温度差 $\Delta T$ 计算出的热应变，转换为等效的节点力。
啊，原来是“将温差引起的膨胀转换为力学荷载”的意思啊！
稳态分析 vs 瞬态分析的判断
温度场用稳态分析就可以吗？还是需要瞬态分析？
判断标准是比较结构的热时间常数和温度变化的速度。热时间常数是 $\tau = \rho c_p L^2 / k$，$L$ 是特征尺寸。如果温度变化相对于 $\tau$ 足够缓慢（即总能几乎跟随稳态），则稳态分析即可。如果是快速加热/冷却，则必须进行瞬态分析。
| 判断标准 | 稳态分析 | 瞬态分析 |
|---|---|---|
| 温度变化速度 | 缓慢（分钟〜小时） | 急剧（秒〜分钟） |
| 典型例子 | 管道的稳态运行 | 发动机启动・停止 |
| 主要关注点 | 最大变形量 | 应力历程・热疲劳 |
| 计算成本 | 低 | 高（与时间步数成正比） |
像涡轮增压器或燃气轮机这样，启动时瞬态温度梯度达到最大的情况，需要通过瞬态分析将各时刻的温度分布传递给结构。最危险的时刻往往不在稳态，而在瞬态，因此“稳态是否足够？”需要慎重判断。
网格
なった
詳しく
報告