第一章 热变形基础：热膨胀系数、热应力与热应变、热弹性力学基本方程
各位工程师朋友，咱们今天聊聊热变形。说实话，我在光机系统这行干了十几年，最头疼的就是温度变化带来的麻烦。你辛辛苦苦设计好的光学系统，温度一变，像差全跑出来了。嗯，这背后的根源，就是热变形。
1.1 热膨胀系数——材料的热敏感度
热膨胀系数，说白了就是材料对温度变化的敏感程度。我习惯用线膨胀系数 α 来表征，单位是 1/℃ 或 ppm/℃。
核心公式：
ΔL = α × L₀ × ΔT其中：ΔL 为长度变化量，L₀ 为初始长度，ΔT 为温度变化量。
举个例子，铝合金的 α 大约 23 ppm/℃，而殷钢只有 1.2 ppm/℃。我在项目中遇到过，一个 200mm 长的铝合金镜筒，温度升高 10℃，长度就变了 0.046mm。你想想看，这对光学系统来说，可不是小数目。
| 材料 | 线膨胀系数 (ppm/℃) | 典型应用 |
|---|---|---|
| 铝合金 6061 | 23.6 | 镜筒、结构件 |
| 殷钢 4J36 | 1.2 | 精密光学基准 |
| 熔石英 | 0.55 | 反射镜基底 |
| 碳纤维复合材料 | ≈0 (可设计) | 轻量化结构 |
我的经验：选材料时别只看 α 值，还要看各向同性。有些复合材料，不同方向的热膨胀系数差好几倍，装调时特别容易出问题。
1.2 热应力与热应变——约束带来的麻烦
如果材料能自由膨胀，那就只有应变，没有应力。但实际结构中，约束无处不在。比如镜片被镜框压着，镜筒两端被固定着。这时候，热膨胀被限制，应力就来了。
热应变 ε_th 的计算很简单：
ε_th = α × ΔT但实际总应变 ε_total 是热应变和弹性应变之和：
ε_total = ε_th + ε_elastic而热应力 σ 则通过胡克定律与弹性应变挂钩：
σ = E × ε_elastic其中 E 是弹性模量。我曾经吃过一次亏，一个铝合金镜框和熔石英镜片，温度变化 20℃ 后，镜片边缘应力直接超过了 10MPa。嗯，后来镜片裂了。从那以后，我每次做设计都会先算一下热应力。
避坑指南：我曾经以为只要材料匹配好就没事，结果忽略了胶粘剂的膨胀系数。光学胶的 α 通常在 50-100 ppm/℃，比金属大得多。温度变化时，胶层会传递很大的应力到镜片表面。
1.3 热弹性力学基本方程——仿真的数学基础
做有限元仿真，本质上就是在解热弹性力学方程。我个人习惯把方程分成三组来理解：
1.3.1 平衡方程
∂σ_x/∂x + ∂τ_xy/∂y + ∂τ_xz/∂z + f_x = 0
∂τ_yx/∂x + ∂σ_y/∂y + ∂τ_yz/∂z + f_y = 0
∂τ_zx/∂x + ∂τ_zy/∂y + ∂σ_z/∂z + f_z = 0说白了，就是力的平衡。每个微元体上的应力梯度加上体力，总和为零。
1.3.2 几何方程
ε_x = ∂u/∂x
ε_y = ∂v/∂y
ε_z = ∂w/∂z
γ_xy = ∂u/∂y + ∂v/∂x
γ_yz = ∂v/∂z + ∂w/∂y
γ_zx = ∂w/∂x + ∂u/∂z这是应变和位移的关系。你想想看，位移场 u、v、w 一旦确定，应变场就全出来了。
1.3.3 本构方程（含热效应）
ε_x = (1/E)[σ_x - ν(σ_y + σ_z)] + αΔT
ε_y = (1/E)[σ_y - ν(σ_x + σ_z)] + αΔT
ε_z = (1/E)[σ_z - ν(σ_x + σ_y)] + αΔT
γ_xy = τ_xy / G
γ_yz = τ_yz / G
γ_zx = τ_zx / G这里多了 αΔT 这一项，就是热应变。G 是剪切模量，ν 是泊松比。
关键点：这三组方程联立起来，加上边界条件，就能求解出温度变化后的位移场和应力场。有限元软件做的就是这件事——把连续方程离散化，然后求解大型线性方程组。
1.4 知识体系框架
为了让大家更直观地理解本章的知识结构，我画了一张图：
1.5 小结
这一章我们聊了热变形的三个核心概念：热膨胀系数是材料的固有属性，热应力来自约束，热弹性力学方程是仿真的理论基础。我个人觉得，理解这三者的关系，比死记公式更重要。你想想看，温度变化引起膨胀，膨胀被约束产生应力，应力又反过来影响变形——这是一个耦合的过程。
做光机系统设计，说白了就是在和热变形做斗争。选对材料、算准应力、用好仿真工具，这三板斧缺一不可。嗯，今天就到这里，下一章我们聊聊有限元仿真的具体操作流程。