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-℃ 几何模型简化 保留:主传力路径 / 光学口径 删除:倒角 / 小孔 / 螺纹 判据:模态偏差 <5% 位移偏差 <3% 材料温度依赖性 E(T)、CTE(T)、k(T)、cp(T) 6061-T6:E 68.9→60 GPa CTE 23.6→25.5 ×10⁻⁶/℃ 熔石英:零膨胀点附近非线性 接触与连接建模 螺栓 → 梁/刚性连接 胶层 → 薄层单元 / 弹簧 接触热阻 1~5×10⁻⁴ m²·K/W 判据:热阻×2 → PV 变化 <λ/20 模型降阶与子模型 全局粗算 → 温度/位移边界 子模型边界距关注区 3~5 倍 Craig-Bampton / Guyan 缩聚 判据:频率 <2%,PV/RMS <5% 五要素共享同一参照系与单位制;任一环节不一致,都会在耦合界面放大为面形误差

六、建模检查清单(可直接使用)

检查项 关键输入 判断依据
坐标系 GCS/LCS 原点、法向、变换矩阵 面形提取方向与光学法向一致
单位制 mm-N-MPa-℃ 封闭单位 热应变与热应力量级无 10³ 偏差
几何简化 保留筋/壁/法兰,删除倒角小孔 前 5 阶模态偏差 <5%,位移偏差 <3%
材料属性 E(T)、CTE(T)、k(T)、cp(T) 表格 工作温区内插值,标注适用范围
连接建模 螺栓/胶层刚度与接触热阻 热阻×2 时 PV 变化 <λ/20
降阶与子模型 主自由度、子模型边界距离 频率误差 <2%,PV/RMS 误差 <5%
建议在项目启动时先写一份“建模约定”文档,固定坐标系、单位制、材料曲线来源与简化判据。后续所有光学、结构、热模型都引用同一份约定,可显著减少耦合界面上的返工。
常见错误:热分析用 m 与 W,结构用 mm 与 N,导致热应变被放大 1000 倍;或只给室温 CTE,使高温工况下的热变形被低估。两者都会让面形结果失去参考价值。

把上述五件事在建模初期一次性对齐,后续的热–结构–光学数据传递才有稳定的基准。统一坐标系解决“往哪映射”,单位制解决“量级对不对”,几何简化解决“算得动”,材料温度依赖解决“算得准”,接触与连接解决“界面不丢信息”,降阶与子模型解决“算得起”。六者共同构成光机热协同建模的基础层。