3. 光机热结构协同建模基础
光机热协同仿真的精度,往往不取决于求解器本身,而取决于建模阶段是否把“几何、材料、连接、规模”四件事处理一致。本节围绕统一坐标系与单位制、几何简化、材料温度依赖性、接触与连接、模型降阶与子模型五个基础问题展开,目标是让光学、结构、热三个物理场在同一套模型语言下可交换、可追溯。
一、统一坐标系与单位制:先定参照,再谈耦合
协同仿真的第一步不是画网格,而是约定坐标系。工程上通常定义三类坐标系:全局坐标系(GCS)用于装配与结果汇总;局部坐标系(LCS)附着于光学元件或镜座,用于描述面形与法向;热分析中还会用到材料主方向坐标系。三者之间必须通过明确的平移向量与旋转矩阵建立映射关系,且全文、全模型保持一致。
以某离轴三反相机为例,主镜顶点在 GCS 下为 (0, 0, 0),次镜顶点为 (0, −320, 180) mm,三镜为 (0, 260, −90) mm。若光学面形用局部坐标描述,则需给出每个镜面的局部原点与法向。热变形结果从结构网格映射到光学面时,必须先把结构节点位移从 GCS 变换到对应 LCS,再沿法向提取面形误差,否则会出现“位移正确、面形错误”的假象。
单位制同样需要统一。推荐采用 mm–N–MPa–s–℃ 的封闭单位制:长度 mm、力 N、应力 MPa、质量 t、密度 t/mm³、热膨胀系数 1/℃、热流密度 mW/mm²。若热分析使用 W 与 m,而结构使用 N 与 mm,则热应变与热应力会在界面处出现 10³ 量级错误。建议在模型说明中固定一张单位对照表,并在脚本入口处做一次量纲检查。
二、几何模型简化原则:保刚度路径,舍细节噪声
几何简化不是“能删就删”,而是判断哪些特征参与刚度传递、哪些只影响局部应力。可遵循三条原则:第一,保留主传力路径上的筋、壁、法兰与安装耳;第二,删除倒角、小孔、螺纹、铭牌槽等对整体刚度贡献小于约 1% 的特征;第三,对光学有效口径内的面形特征必须保留,口径外可合并。
例如镜筒上均布的 24 个 M4 螺纹孔,若用于连接,可简化为等效刚性连接或梁单元;若仅用于工艺,可直接删除。蜂窝夹层板可等效为正交各向异性壳,但需给出等效厚度与等效弯曲刚度。简化后应做一次模态或静刚度对比:若前 5 阶模态频率偏差小于 5%、关键点位移偏差小于 3%,可认为简化可接受。这一判断依据比“看起来差不多”更可靠。
三、材料属性温度依赖性:把 CTE 和 E 写成温度函数
光机结构常用材料在 −20 ℃~+60 ℃ 范围内性能会变化。以铝合金 6061-T6 为例,20 ℃ 时弹性模量约 68.9 GPa,200 ℃ 时降至约 60 GPa;热膨胀系数从 23.6×10⁻⁶/℃ 增至约 25.5×10⁻⁶/℃。若仿真只输入室温常数,热变形会被系统性低估或高估。
正确做法是:在材料卡片中把 E、CTE、导热系数 k、比热容 cp 都写成温度表格,求解器在每步按当前温度插值。对于光学玻璃如熔石英,CTE 在 20 ℃ 附近接近 0.5×10⁻⁶/℃,但到 200 ℃ 会明显上升,且存在“零膨胀点”附近非线性,必须用实测或手册曲线而非单点值。若缺少数据,可用线性近似,但需在报告中标注适用范围与不确定度。
四、接触与连接建模:刚度与热阻要成对考虑
光学组件中大量使用螺栓、胶接、销钉与柔性铰。结构上,螺栓可简化为梁单元或刚性连接;胶层可简化为薄层单元或弹簧;柔性铰则需保留其转动刚度。热分析中,接触面存在接触热阻,不能默认理想导热。典型螺栓连接面接触热阻约为 1×10⁻⁴~5×10⁻⁴ m²·K/W,胶层热阻由厚度除以导热系数得到。
建议对每个连接点同时定义力学刚度与热阻,并在同一位置建立“结构–热”耦合节点。若只做结构分析,可先忽略热阻;一旦进入热–结构耦合,必须把热阻作为参数扫描项,观察其对镜面温度梯度与面形的影响。判断依据是:当接触热阻变化 2 倍时,若关键面形 PV 变化超过 λ/20,则该连接必须精细建模。
五、模型降阶与子模型技术:先全局粗算,再局部精算
整机模型动辄数百万单元,直接热–结构耦合代价高。常用策略是:全局模型用粗网格或超单元降阶,获得整体温度场与位移边界;再在光学元件或关键连接处切出子模型,把全局结果作为边界条件插值到子模型上,进行精细网格计算。子模型边界应远离关注区,一般取特征尺寸的 3~5 倍。
降阶方法包括模态综合、Craig–Bampton 与 Guyan 缩聚。以某支撑结构为例,全局模型 120 万单元,缩聚后保留 80 个主自由度,前 20 阶模态频率误差小于 2%,计算时间从 4 h 降至 12 min。子模型再取镜座局部 15 万单元,面形计算精度可恢复到与全精细模型相差 5% 以内。判断是否可接受,应同时看频率误差、关键位移误差与面形 PV/RMS 误差三项。
六、建模检查清单(可直接使用)
| 检查项 | 关键输入 | 判断依据 |
|---|---|---|
| 坐标系 | GCS/LCS 原点、法向、变换矩阵 | 面形提取方向与光学法向一致 |
| 单位制 | mm-N-MPa-℃ 封闭单位 | 热应变与热应力量级无 10³ 偏差 |
| 几何简化 | 保留筋/壁/法兰,删除倒角小孔 | 前 5 阶模态偏差 <5%,位移偏差 <3% |
| 材料属性 | E(T)、CTE(T)、k(T)、cp(T) 表格 | 工作温区内插值,标注适用范围 |
| 连接建模 | 螺栓/胶层刚度与接触热阻 | 热阻×2 时 PV 变化 <λ/20 |
| 降阶与子模型 | 主自由度、子模型边界距离 | 频率误差 <2%,PV/RMS 误差 <5% |
把上述五件事在建模初期一次性对齐,后续的热–结构–光学数据传递才有稳定的基准。统一坐标系解决“往哪映射”,单位制解决“量级对不对”,几何简化解决“算得动”,材料温度依赖解决“算得准”,接触与连接解决“界面不丢信息”,降阶与子模型解决“算得起”。六者共同构成光机热协同建模的基础层。