4. 热应力产生机理：温度梯度引起的热膨胀不均、热传导方程、热应力计算基础
各位，咱们今天聊热应力。说实话，这玩意儿在精密光学制造里，是个绕不开的坎儿。我做了这么多年光学加工，见过太多因为热应力导致镜片变形、甚至炸裂的案例。说白了，热应力的根源就一句话：温度不均匀，材料又不想被束缚，结果内部就打架了。
4.1 温度梯度：热应力的“始作俑者”
你想想看，一块光学玻璃放在炉子里加热。表面先热起来，内部还凉着。表面想膨胀，内部不让它膨胀——这不就产生应力了吗？反过来，冷却的时候也一样，表面先收缩，内部还热着，又得打架。
我个人习惯把温度梯度想象成“热膨胀的错配”。就像两个人跳舞，一个快一个慢，步子肯定要踩到对方脚上。这个“踩脚”的力，就是热应力。
核心要点：温度梯度越大，热应力越严重。尤其是在光学玻璃的退火过程中，温度梯度控制不好，轻则折射率不均匀，重则直接报废。
我记得有一次，客户送来一块直径300mm的K9玻璃，要求面形精度λ/10。结果退火时升温速率快了那么一点点，中心温度比边缘低了15°C。抛光后一测，面形直接变成了“马鞍形”。嗯，从那以后，我对温度梯度的控制就再也不敢马虎了。
4.2 热传导方程：温度场的“数学画像”
要定量分析热应力，首先得知道温度是怎么分布的。这就得请出热传导方程了。
对于各向同性的均匀材料，三维热传导方程长这样：
∂T/∂t = α · (∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²) + Q/(ρ·Cp)其中：
- T：温度（K）
- t：时间（s）
- α：热扩散系数（m²/s），α = k/(ρ·Cp)
- k：导热系数（W/(m·K)）
- ρ：密度（kg/m³）
- Cp：比热容（J/(kg·K)）
- Q：内热源项（W/m³）
说白了，这个方程就是在描述热量怎么在材料里“跑”。左边是温度随时间的变化，右边是热量扩散的贡献加上内部热源。
实战经验：在光学加工中，我们通常关注稳态温度场（∂T/∂t = 0）。因为很多工艺过程，比如精密退火，都是缓慢进行的，可以近似看作稳态。这时候方程简化为拉普拉斯方程：∇²T = 0。
我曾经用有限元软件算过一个透镜在磨削过程中的温度场。磨削区温度高达200°C，但距离磨削点仅5mm的地方，温度就降到了60°C。这个梯度，你想想看，有多恐怖？
4.3 热应力计算基础：从温度到应力的“桥梁”
有了温度场，怎么算热应力？这里需要引入热弹性力学的基本方程。
首先，热应变由两部分组成：
- 弹性应变：由应力引起，服从胡克定律
- 热应变：由温度变化引起，ε_th = α · ΔT
总应变 = 弹性应变 + 热应变。写成张量形式：
ε_ij = (1+ν)/E · σ_ij - (ν/E) · σ_kk · δ_ij + α · ΔT · δ_ij其中：
- E：弹性模量（Pa）
- ν：泊松比
- α：线膨胀系数（1/K）
- ΔT：温度变化量（K）
- δ_ij：克罗内克符号
这个公式看着复杂，其实核心思想很简单：如果材料能自由膨胀，就不会有热应力。但实际中，材料内部各部分相互约束，或者外部有边界约束，热膨胀被“憋”住了，应力就出来了。
注意：对于光学玻璃，α值通常很小（5~10×10⁻⁶/K），但E值很大（70~90GPa）。所以即使ΔT只有几十度，产生的热应力也可能达到几十兆帕。而光学玻璃的抗拉强度通常只有30~50MPa。嗯，你懂的，一不小心就裂了。
4.4 知识体系：一张图说清楚
下面我用一张SVG图，把热应力的产生逻辑串起来。这样你看起来更直观。
4.5 一个简单的算例
咱们来算个简单的。假设一块BK7玻璃板，长100mm，厚10mm。从20°C均匀加热到120°C，但玻璃板两端被刚性夹具固定，不能自由伸长。
BK7的参数：
- E = 81 GPa
- α = 7.1 × 10⁻⁶ /K
- ν = 0.208
如果完全约束，热应力为：
σ = -E · α · ΔT
= -81×10⁹ × 7.1×10⁻⁶ × 100
= -57.51 MPa负号表示压应力。57.5 MPa，已经接近BK7的抗压强度极限了。如果温度梯度再大一点，或者玻璃有微裂纹，那就危险了。
我的建议：在实际工艺中，尽量让温度变化平缓。比如退火时，升温速率控制在0.5°C/min以内。另外，夹具设计要留出膨胀余量，别把玻璃“卡死”。
4.6 避坑指南
我曾经犯过一个低级错误。有一次做胶合透镜，用紫外胶把两个镜片粘在一起。固化时用了大功率紫外灯，结果镜片局部温度飙升到80°C，而边缘还是室温。胶层还没固化，镜片就因为热应力裂了。嗯，从那以后，我固化胶水时都会先用小功率预热，或者用风扇吹着降温。
所以，各位记住：热应力不是算出来的，是设计出来的。你设计工艺时，就得把温度梯度考虑进去。
公众号:蓝海资料掘金营，微信deep3321