4. 热-力耦合基础：热传导方程、热应力产生机理、热膨胀系数的影响

各位同行，咱们今天聊聊热-力耦合。说实话，做复合材料叶片，最让我头疼的往往不是单纯的力学问题，而是温度变化带来的那一堆麻烦。你想想看，叶片从模具里脱模、冷却，再到实际运行中经历高温燃气冲刷，这中间的温度变化少说也有几百度。温度一变，应力就跟着来了，搞不好就开裂、变形。

我个人习惯把热-力耦合拆成三块来理解：热量怎么传进去、温度怎么产生应力、材料本身怎么响应。咱们一个一个说。

4.1 热传导方程：热量是怎么跑的？

先说热传导。说白了，就是热量在材料内部怎么扩散。对于复合材料叶片这种各向异性材料，热传导比金属复杂得多。我记得刚入行时，用各向同性的公式去算碳纤维预浸料的温度场，结果偏差很大——因为纤维方向导热快，垂直方向导热慢，差了好几倍。

热传导方程的一般形式是这样的：

ρ·Cp·(∂T/∂t) = ∂/∂x(kx·∂T/∂x) + ∂/∂y(ky·∂T/∂y) + ∂/∂z(kz·∂T/∂z) + Q

这里：

• ρ — 密度，单位 kg/m³
• Cp — 比热容，单位 J/(kg·K)
• kx, ky, kz — 三个方向的热导率，单位 W/(m·K)
• Q — 内热源，比如固化反应放热

对于复合材料，kx、ky、kz 往往不一样。比如单向碳纤维预浸料，沿纤维方向（0°）的热导率可能是垂直方向的 5~10 倍。这一点在有限元建模时必须单独定义材料坐标系，否则算出来的温度场就是错的。

我在做风电叶片固化仿真时遇到过一个问题：模具加热从底部传上来，但叶片厚度方向导热慢，导致内部温度滞后。结果固化度不均匀，产生了内应力。后来我调整了升温速率，才把这个问题压下去。

4.2 热应力产生机理：温度变化为什么会产生力？

热应力的本质其实很简单：材料受热会膨胀，受冷会收缩。但如果膨胀或收缩被约束住了，应力就产生了。

举个例子。你把一根金属棒两端固定，然后加热。金属棒想伸长，但两端被卡住，伸不出去。于是棒内部就产生了压应力。反过来，冷却时产生拉应力。这就是热应力的基本来源。

对于复合材料叶片，情况更复杂：

• 不同材料的热膨胀系数不同 — 比如碳纤维是负膨胀（受热收缩），树脂是正膨胀（受热膨胀）。两者粘在一起，温度一变，界面处就会产生剪切应力。
• 不同铺层方向的热膨胀不同 — 0°层和90°层的热变形方向不一样，层间会产生应力。
• 厚度方向温度梯度 — 叶片表面温度高、内部温度低，内外膨胀不一致，产生弯曲应力。

我曾经做过一个案例：某型风机叶片在脱模后出现表面裂纹。排查了很久，发现是冷却速率太快，外层已经硬化收缩，内层还在高温膨胀状态，结果外层被拉裂。后来把冷却速率从 5°C/min 降到 1.5°C/min，问题就解决了。

4.3 热膨胀系数的影响：材料自己的“脾气”

热膨胀系数（CTE，Coefficient of Thermal Expansion）是衡量材料受热后尺寸变化程度的参数。单位是 1/°C 或 ppm/°C。

对于复合材料，CTE 同样具有各向异性：

材料	纤维方向 CTE (×10⁻⁶/°C)	垂直方向 CTE (×10⁻⁶/°C)
碳纤维/环氧（T300/5208）	-0.3 ~ 0.5	25 ~ 35
玻璃纤维/环氧	5 ~ 8	20 ~ 30
纯环氧树脂	50 ~ 70	50 ~ 70

看到没？碳纤维在纤维方向甚至是负膨胀——温度升高反而收缩。这听起来反直觉，但确实如此。原因在于碳纤维的晶体结构在受热时沿轴向收缩。

这就带来一个有意思的现象：碳纤维复合材料叶片在升温时，纤维方向收缩，树脂方向膨胀。两者互相拉扯，界面处产生很大的剪切应力。如果界面强度不够，就会脱粘。

我个人习惯在有限元软件中这样设置：

*MATERIAL, NAME=CFRP
*ELASTIC, TYPE=ENGINEERING CONSTANTS
155000, 8500, 8500, 0.3, 0.3, 0.4, 4700, 4700, 3200
*EXPANSION, TYPE=ORTHO
-0.3e-6, 28e-6, 28e-6, 0, 0, 0
*CONDUCTIVITY, TYPE=ORTHO
6.5, 0.8, 0.8
*SPECIFIC HEAT
850
*DENSITY
1600

这段代码定义了一个正交各向异性的碳纤维复合材料，其中 *EXPANSION 那行就是三个方向的热膨胀系数。注意第一个值是负的，代表纤维方向收缩。

4.4 知识体系总览

下面这张图是我自己整理的，把热-力耦合的核心逻辑串起来了。你看一遍应该就能明白整个分析流程：

这张图把热-力耦合的核心逻辑串起来了。你看，左边是热传导，中间是热膨胀，右边是热应力。三者顺序耦合：先算温度场，再算热应变，最后算热应力。下面三个框是我总结的关键因素、常见失效和工程对策，做分析时对照着看，不容易漏项。

好了，关于热-力耦合的基础就聊到这儿。说白了，就是搞清楚热量怎么进去、材料怎么变形、约束怎么产生应力。这三件事搞明白了，叶片制造过程中的热应力问题就解决了一大半。

---

公众号:蓝海资料掘金营，微信deep3321