4. 热弹性力学基础：热应力与热应变的基本概念、热传导方程、热力耦合问题的控制方程

各位同学，欢迎来到《复合材料叶片热力耦合分析实战》的第四讲。

说实话，热力耦合这块内容，我当年刚接触时也觉得挺绕的。又是温度场，又是应力场，两个东西搅在一起，脑子容易打结。但你别怕，咱们今天就把这层窗户纸捅破。

这一章，我带你从最基础的概念开始，一步步搭建起热弹性力学的知识框架。说白了，就是搞清楚三件事：温度变化怎么让材料变形？热量在材料里怎么跑？这两个过程怎么互相影响？

4.1 热应力与热应变的基本概念

先问大家一个问题：一根铁棍，你把它加热，它会变长对吧？但如果我把它两头死死固定住，不让它伸长，会发生什么？

嗯，它内部会产生很大的力，这就是热应力。

热应变其实很好理解。材料受热后，分子运动加剧，原子间距变大，宏观上就表现为膨胀。这个膨胀量就是热应变。对于各向同性材料，热应变只与温度变化量有关：

ε_th = α · ΔT

其中：

• α 是线膨胀系数（单位：1/℃）
• ΔT 是温度变化量（单位：℃）

但注意了，复合材料可不是各向同性的。我在做风机叶片项目时，就吃过这个亏。碳纤维增强树脂基复合材料，沿着纤维方向的热膨胀系数可能只有 1×10⁻⁶ /℃，而垂直纤维方向可能高达 30×10⁻⁶ /℃。你想想看，这差别有多大？

所以对于复合材料，热应变张量必须写成矩阵形式：

{ε_th} = [α] · ΔT

这里的 [α] 是一个 6×6 的矩阵，包含了各个方向的热膨胀系数。说白了，复合材料的热变形是各向异性的，你加热它，它可能不仅会膨胀，还会发生剪切变形。

4.2 热传导方程

搞清楚了热应变，接下来我们看看热量是怎么在材料内部传递的。这就是热传导方程要解决的问题。

对于各向异性材料，傅里叶热传导定律可以写成：

q = -[k] · ∇T

其中：

• q 是热流密度向量
• [k] 是导热系数矩阵（3×3）
• ∇T 是温度梯度

对于复合材料叶片，导热系数同样具有方向性。纤维方向的导热系数可能比基体方向高一个数量级。这意味着热量会优先沿着纤维方向传递。

结合能量守恒定律，我们可以得到瞬态热传导方程：

ρ·c·∂T/∂t = ∇·([k]·∇T) + Q

这里：

• ρ 是密度
• c 是比热容
• Q 是内部热源（比如固化反应放热）

这个方程看着复杂，其实说白了就是一句话：单位时间内，材料内部积蓄的热量 = 传入的热量 + 内部产生的热量。

4.3 热力耦合问题的控制方程

好了，现在我们把热和力这两条线串起来。

热力耦合问题的核心，就是同时求解两个场：温度场和位移场。它们之间通过热应变和热传导相互影响。

完整的控制方程包括三部分：

1. 平衡方程（力的平衡）：

∇·σ + f = 0

2. 几何方程（应变与位移的关系）：

ε = (∇u + (∇u)ᵀ)/2

3. 本构方程（应力与应变、温度的关系）：

σ = [C] · (ε - ε_th)

其中 [C] 是材料的刚度矩阵。对于复合材料，这个矩阵也是各向异性的，需要根据铺层角度和材料属性计算。

把这三个方程联立起来，再加上热传导方程，就构成了热力耦合问题的完整控制方程组。在有限元分析中，我们通常采用顺序耦合或完全耦合两种方式求解。

我个人习惯，如果温度变化对材料属性影响不大，就用顺序耦合——先算温度场，再把温度结果作为载荷施加到结构分析中。但如果温度变化剧烈，或者材料属性随温度变化明显，那就必须用完全耦合了。

4.4 知识体系框架图

为了让你更直观地理解这一章的知识结构，我画了一张图：

4.5 小结与思考

这一章的内容，说白了就是三个关键词：热应变、热传导、耦合方程。你只要记住：

• 温度变化 → 产生热应变 → 如果约束住 → 产生热应力
• 热量在材料中传递，遵循傅里叶定律，复合材料有方向性
• 热和力是双向耦合的，分析时根据情况选择顺序或完全耦合

嗯，到这里，热弹性力学的基础就讲完了。下一章我们会把这些理论落实到有限元实现中，到时候我会带你手写一段 Python 代码，把今天讲的方程离散化，真正跑起来。

---

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