3. 热稳定性的工程定义与量化
做光机热分析，第一步就是搞清楚——热从哪来。
很多人一上来就建模、加边界条件，结果算出来温度场跟实测差一大截。为什么？热源类型搞错了，或者载荷施加方式不对。我早年吃过这个亏，一个激光投影系统，仿真温升比实测低了8℃，后来发现是忽略了腔体内的辐射换热。
这一节，咱们就把热源和热载荷这件事彻底捋清楚。
3.1 光学系统中的三大热源
光学系统里的热源，说白了就三类：主动光源发热、环境热辐射、结构自身耗散。其中前两类是主要矛盾。
| 热源类型 | 典型场景 | 热流密度范围 | 时间特性 |
|---|---|---|---|
| 激光光源 | 激光雷达、激光投影、光纤耦合 | 10⁴ ~ 10⁶ W/m² | 连续/脉冲 |
| LED光源 | 照明、投影、显示背光 | 10² ~ 10⁴ W/m² | 连续/调制 |
| 环境辐射 | 户外设备、卫星、高温车间 | 10¹ ~ 10³ W/m² | 准稳态/变化 |
3.1.1 激光光源——高密度点热源
激光的特点是能量高度集中。一个10W的半导体激光器，发光面可能只有几百微米，热流密度轻松上百万W/m²。这种热源对光学镜片的局部加热效应非常明显。
我在做激光雷达项目时遇到过一个问题：激光通过分光棱镜后，棱镜表面温度在几秒内上升了15℃，导致透射波前畸变。后来分析发现，棱镜镀膜对808nm激光的吸收率虽然只有0.5%，但架不住功率高、光斑小。
3.1.2 LED光源——面热源与热沉耦合
LED跟激光不一样，它是大面积发光，热流密度相对低，但总功率不小。一个100W的COB LED，发光面积约1~2cm²，热流密度在10⁴ W/m²量级。
LED热分析的核心难点在于热沉设计。LED芯片产生的热量，90%以上要通过热沉传导出去。如果热沉与散热器之间的接触热阻没处理好，结温会飙升，光效下降，寿命缩短。
3.1.3 环境辐射——容易被忽略的“隐形热源”
环境辐射这东西，很多人做仿真时会忽略。但实际项目中，它往往是导致系统温升的主要贡献者。
举个例子：一个户外使用的光学监控设备，外壳被太阳直射，表面温度能达到60~70℃。这个热量通过外壳传导到内部光学件，导致镜片温度升高、像质下降。如果不考虑太阳辐射，仿真结果跟实测能差10℃以上。
3.2 热载荷的施加方式
热源确定了，接下来就是怎么把热载荷加到有限元模型上。我见过不少工程师，直接把热流密度加在节点上，结果算出来的温度场奇奇怪怪。其实，热载荷施加有四种标准方式，每种对应不同的物理场景。
3.2.1 热流密度载荷（Heat Flux）
适用于已知单位面积热功率的场景。比如激光光斑照射在镜面上，或者LED芯片贴在热沉上。
施加时要注意：热流密度是矢量，方向垂直于加载面。如果光斑是斜入射，需要分解法向分量。
! ANSYS APDL 示例：在镜面施加高斯热流
! 光斑半径 2mm，中心功率密度 5e5 W/m²
NSEL, S, LOC, X, -0.002, 0.002
NSEL, R, LOC, Y, -0.002, 0.002
*GET, NMIN, NODE, 0, NUM, MIN
*GET, NMAX, NODE, 0, NUM, MAX
*DO, I, NMIN, NMAX
*IF, NSEL(I), EQ, 1, THEN
DIST = SQRT(NX(I)**2 + NY(I)**2)
Q = 5E5 * EXP(-DIST**2 / (2*0.001**2))
F, I, HEAT, Q * AREA_NODE(I)
*ENDIF
*ENDDO
3.2.2 体热源载荷（Heat Generation Rate）
适用于体积发热的场景，比如激光晶体内部吸收泵浦光、LED芯片的PN结发热、电阻发热等。
体热源的单位是W/m³。施加时，需要知道发热体的体积和总发热功率。如果发热体形状不规则，可以用平均体热源近似，但局部精度会下降。
| 场景 | 总功率 | 体积 | 体热源密度 |
|---|---|---|---|
| 激光晶体（Nd:YAG） | 50 W | 3×3×20 mm³ | 2.78×10⁸ W/m³ |
| LED芯片（单颗） | 1 W | 1×1×0.1 mm³ | 1×10¹⁰ W/m³ |
| 电阻加热片 | 10 W | 10×10×0.5 mm³ | 2×10⁷ W/m³ |
3.2.3 温度边界条件（Temperature BC）
适用于已知某部分结构温度的场景。比如光学系统安装在恒温平台上，或者某些部件与散热器接触面温度已知。
温度边界条件是最强约束，它会强制节点温度等于给定值。使用时要谨慎——如果给错了位置，整个温度场都会偏。
3.2.4 对流与辐射边界条件
这是最灵活也最容易出错的载荷方式。对流边界需要给定换热系数和环境温度；辐射边界需要给定发射率和环境辐射温度。
我见过一个案例：工程师把自然对流的换热系数设为10 W/m²·K，但实际系统有强制风冷，换热系数应该在50~100 W/m²·K。结果仿真温升比实测高了20℃。
- 自然对流（空气）：5~15 W/m²·K
- 强制风冷（2~5 m/s）：20~100 W/m²·K
- 强制风冷（10~20 m/s）：100~500 W/m²·K
- 水冷：500~5000 W/m²·K
3.3 热载荷施加的工程要点
说了这么多，最后总结几条我在项目中反复踩过的坑：
- 热源位置要精确——激光光斑位置偏差0.5mm，镜片温度分布就完全不一样。我习惯在CAD模型里用坐标系标记光斑中心，确保网格划分后节点位置准确。
- 吸收率不是猜的——镀膜吸收率可以从镀膜厂商拿到，一般在0.1%~1%之间。如果没有数据，可以用分光光度计实测。不要随便用0.5%代替，不同膜系差异很大。
- 多物理场耦合——热源往往不是独立的。比如激光加热导致镜片温度升高，镜片热膨胀改变光路，光斑位置偏移，又反过来影响热源分布。这种热-结构-光耦合问题，需要迭代求解。
- 瞬态 vs 稳态——激光脉冲加热要用瞬态分析，时间步长要小于脉冲宽度的1/10。连续激光或LED可以用稳态分析，但要注意开机瞬间的温升过程。