3. 热弹性效应：热应力与热应变的基本理论、热弹性方程、材料的热力学本构关系

各位工程师朋友，咱们今天聊聊热弹性效应。说白了，就是温度一变，材料不光会膨胀收缩，还会因此产生内部应力。这个效应，是光学系统无热化设计的根基。你想想看，一个镜头从常温装调到低温环境，镜片和镜筒的膨胀系数不一样，光路就偏了。嗯，这就是热弹性在捣鬼。

3.1 热应力与热应变：从一根杆子说起

先看最简单的模型。一根杆子，两端固定。温度升高，它想伸长，但被约束住了。这时候杆子内部就会产生压应力。反过来，温度降低，它想收缩，但被拉住了，就会产生拉应力。

这个现象，我习惯用两个概念来拆解：

• 自由热应变：如果杆子没被约束，温度变化 ΔT 引起的应变为 ε_th = α · ΔT。α 是线膨胀系数。
• 约束热应力：因为约束的存在，实际应变被“憋”住了，于是产生了应力 σ = -E · α · ΔT。E 是弹性模量。

我在项目中遇到过一个问题：一个铝合金镜筒，里面装着一片玻璃镜片。低温试验时，镜片被镜筒“夹”得太紧，产生了双折射。后来一算，就是热应力超过了玻璃的弹性极限。所以，别小看这个一维公式，它是所有复杂分析的起点。

3.2 热弹性方程：三维世界的扩展

一维太简单了，实际镜片是三维的。我们需要把上面的公式推广到三维空间。这就是热弹性方程，也叫杜哈梅尔-诺伊曼方程。

三维情况下，应变和应力的关系变成了张量形式。我直接给出最常用的形式：

ε_ij = (1+ν)/E · σ_ij - (ν/E) · σ_kk · δ_ij + α · ΔT · δ_ij

其中：

• ε_ij 是应变张量分量
• σ_ij 是应力张量分量
• σ_kk 是应力张量的迹（三个正应力之和）
• ν 是泊松比
• δ_ij 是克罗内克符号（i=j时为1，否则为0）

这个方程看着复杂，其实逻辑很简单：总应变 = 弹性应变 + 热应变。弹性应变由应力引起，热应变由温度引起。两者线性叠加。

3.3 材料的热力学本构关系

上面说的热弹性方程，其实是一种本构关系。本构关系，说白了就是材料“受力-变形-温度”三者之间的内在联系。不同材料，这个关系不一样。

我整理了一个常用光学材料的热力学参数对比表，方便你快速查阅：

材料	弹性模量 E (GPa)	泊松比 ν	线膨胀系数 α (10⁻⁶/℃)	热导率 k (W/m·K)
熔石英 (SiO₂)	72	0.17	0.55	1.4
BK7 玻璃	81	0.21	7.1	1.1
铝合金 6061	69	0.33	23.6	167
钛合金 TC4	110	0.34	8.6	7.2
殷钢 (Invar 36)	141	0.29	1.2	10

你看这个表，熔石英的膨胀系数极低，只有0.55，所以它热稳定性好。但它的热导率也低，只有1.4，意味着温度均匀化很慢。铝合金膨胀系数大，但热导率高，温度响应快。这就是为什么很多无热化设计会采用“铝合金镜筒+殷钢嵌件”的组合——用殷钢的低膨胀来抵消铝合金的高膨胀。

3.4 热弹性效应的知识框架

为了让你更直观地理解本章的知识体系，我画了一张结构图。它把热弹性效应的核心逻辑串了起来：

这张图把整个逻辑串起来了：温度变化是源头，自由膨胀和约束条件是两条路径，最终汇入热弹性本构方程，再指导无热化设计。你把这个框架记在脑子里，后面分析具体案例时就不会乱。

3.5 实际应用中的几点提醒

最后，我结合项目经验，给你几个实用建议：

• 线弹性假设的局限：热弹性方程假设材料是线弹性的。但温度变化很大时，材料可能进入塑性。比如一些塑料镜筒，低温下会变脆，高温下会蠕变。这时候线弹性理论就不够用了。
• 各向异性材料：我遇到过用碳纤维复合材料做镜筒的项目。这种材料在不同方向上的膨胀系数不一样，甚至可以是负值。这时候热弹性方程要写成矩阵形式，不能简单用标量α。
• 瞬态 vs 稳态：热应力分析分两种。稳态分析看最终状态，瞬态分析看温度变化过程中的应力峰值。我记得有一次做空间相机设计，发射阶段的快速温变产生了很大的瞬态热应力，差点把镜片震裂。所以，别只看稳态。