4. 光机结构有限元建模要点

拓扑优化能否给出可信的轻量化方案,取决于有限元模型是否真实反映光机结构在工作状态下的力学行为。若网格、边界条件、连接模拟、材料属性或光学面形评价节点选取存在偏差,优化算法会在错误模型上“高效”地收敛到一个看似减重、实则无法满足光学指标的结果。本节围绕这五个建模要点,说明每一步应如何落地,以及如何判断模型是否具备工程可信度。

一、网格划分:先保证传力路径,再谈单元数量

光机结构通常由镜筒、支撑环、柔性铰链、背板、反射镜镜体等组成,几何尺度差异大。网格划分的核心不是“越密越好”,而是让单元尺寸与结构传力特征匹配。对以弯曲、扭转为主的镜筒和支撑臂,建议采用二阶四面体或六面体单元,并在壁厚方向至少布置 2~3 层单元,否则弯曲刚度会被高估。对薄壁、柔性铰链等局部,需单独加密,但过渡区应设置 1.5~2 倍尺寸渐变,避免相邻单元尺寸突变导致应力集中假象。

以某口径 200 mm 反射镜支撑结构为例,镜体采用 8 mm 厚背部半刚性支撑,镜筒壁厚 5 mm。若整体用 10 mm 单元,镜体一阶固有频率计算值可能比 5 mm 局部加密模型高 8%~12%,原因是背部筋板弯曲变形未被充分捕捉。因此建议先做网格无关性检查:将全局单元尺寸从 10 mm 依次减至 8、6、5 mm,观察一阶频率与镜面 RMS 变化;当连续两次结果差异小于 3% 时,可认为网格收敛。对光学面形敏感区域,评价节点所在表面的单元尺寸应不大于镜面直径的 1/40。

二、边界条件:把“安装状态”翻译成约束

边界条件决定载荷如何传入结构,是建模中最容易被简化的环节。实际光机结构通过法兰、螺钉、销钉或柔性环节与整机连接,不能简单用“全约束”代替。若支撑点位于镜筒底部法兰,应约束法兰安装面的法向位移与两个切向位移,并释放绕光轴的转动,以模拟螺栓连接提供的摩擦与局部刚度。若采用三点柔性支撑,则每个支撑点可简化为沿支撑杆方向的弹簧单元,刚度由柔性杆轴向刚度 k=EA/L 估算,其中 E 为材料弹性模量,A 为柔性杆截面积,L 为有效长度。

重力工况下,需明确重力方向与光轴关系。例如光轴水平时,重力沿垂直光轴方向加载;光轴竖直时,重力沿光轴方向加载。温度工况则通过参考温度与均匀温升施加,若结构存在材料差异,还需考虑热膨胀系数不匹配引起的热应力。判断边界条件是否合理,可检查约束反力是否与外部载荷平衡,以及支撑点附近是否出现非物理的应力尖峰。

三、连接模拟:螺钉、胶接与柔性环节不能“刚性绑定”

光机结构中大量使用螺钉连接、胶接和柔性铰链。若全部用刚性绑定,局部刚度会被显著高估,优化结果可能过度削弱连接区。对螺钉连接,可采用梁单元或弹簧单元模拟:将螺钉简化为一段等效梁,直径取螺纹小径,长度取夹紧长度,两端分别与连接件节点耦合。对胶接面,可用薄层实体单元或界面弹簧,胶层厚度按实际 0.05~0.2 mm 建模,弹性模量取胶粘剂供应商数据。

柔性铰链是光机结构中的关键环节,其转动刚度远低于相邻实体。若用实体单元直接划分,需保证铰链最薄处至少有 3 层单元;若用弹簧单元替代,则需根据铰链几何计算转动刚度。以圆弧柔性铰链为例,其绕转轴的转动刚度可近似为 kθ=Ebt^2.5/(9πR^0.5),其中 b 为宽度,t 为最薄处厚度,R 为圆弧半径。连接模拟完成后,应通过模态试验或静力试验校核:若一阶频率误差超过 10%,需重新检查连接刚度。

四、材料属性:各向异性与温度依赖性不可忽略

光机结构常用材料包括铝合金、钛合金、碳纤维复合材料、微晶玻璃和铟钢。金属材料可视为各向同性,但需注意弹性模量随温度变化;例如铝合金 6061 在 20℃ 时 E≈68.9 GPa,在 100℃ 时降至约 65 GPa。碳纤维复合材料则必须输入正交各向异性参数:E1、E2、E3、ν12、ν13、ν23、G12、G13、G23,并定义铺层方向。若将碳纤维镜筒简化为各向同性,轴向刚度可能被低估 30% 以上。

光学镜体材料如微晶玻璃,其热膨胀系数接近零,但弹性模量约 90 GPa,密度约 2.5 g/cm³。建模时应区分镜体与支撑结构材料,避免统一赋值。对温度工况,还需输入热膨胀系数和热导率;若只做静力分析,至少应输入弹性模量、泊松比、密度和热膨胀系数。材料数据应来自供应商手册或实测,并在模型中记录来源,便于后续复核。

五、光学面形评价节点选取:从“看应力”转向“看面形”

拓扑优化的目标通常不是降低最大应力,而是保证光学面形精度。因此,必须在镜面加工表面上选取评价节点,而不是在镜体内部或支撑点上。评价节点应覆盖镜面有效口径,建议按极坐标网格布置:沿径向取 5~7 圈,沿周向每圈取 12~24 点,总节点数约 100~200 个。节点位置应与干涉仪测量点对应,便于后续对比。

评价指标包括 PV(峰谷值)和 RMS(均方根值)。在有限元中,提取这些节点在重力、温度或装配载荷下的法向位移,去除刚体位移后计算 PV 和 RMS。例如某 200 mm 反射镜要求重力工况下面形 RMS≤λ/40(λ=632.8 nm,即约 15.8 nm)。若优化后 RMS 为 12 nm,则满足;若为 20 nm,则需调整支撑点位置或增加局部刚度。评价节点选取后,应在模型中固定其编号,避免每次优化后重新选取导致结果不可比。

光机结构有限元建模五要点与优化可信度闭环 有限元模型 网格 · 边界 · 连接 材料 · 评价节点 ① 网格划分 壁厚≥2~3层单元 收敛判据:频率差<3% ② 边界条件 法兰约束/弹簧支撑 重力方向与光轴关系 ③ 连接模拟 螺钉梁/胶层弹簧 柔性铰链刚度kθ ④ 材料属性 各向异性/温度依赖 E、ν、ρ、α ⑤ 光学面形评价节点选取 镜面有效口径 · 径向5~7圈 · 周向12~24点 · PV/RMS 优化结果判据 一阶频率误差 ≤10% 面形 RMS ≤ λ/40 约束反力与载荷平衡 模型校核与迭代 频率/面形偏差超限 → 检查连接刚度 → 加密评价区网格 注:图中参数为教学示例,实际取值需依据材料手册、试验数据与光学指标确定

建模检查清单与判断依据

完成上述五步后,建议用一份简短的检查清单确认模型可信度。以下清单可直接用于项目评审或自检:

检查项 合格判据 常见问题
网格收敛性 连续两次加密后一阶频率差<3%,面形RMS差<5% 壁厚方向单元不足,弯曲刚度高估
边界条件 约束反力合力与外部载荷平衡,误差<2% 全约束代替柔性支撑,热应力失真
连接刚度 模态试验一阶频率误差≤10% 螺钉刚性绑定,局部刚度过高
材料属性 各向异性参数完整,温度工况含E(T)、α(T) 碳纤维按各向同性赋值
评价节点 覆盖有效口径,节点数与干涉仪测点对应 取镜体内部节点,面形评价失真
实用提示:若优化后镜面RMS不达标,优先检查支撑点位置与连接刚度,而不是直接增加镜体厚度。多数光机结构的面形误差来自支撑局部变形,而非镜体整体刚度不足。通过调整柔性铰链刚度或增加支撑点数量,往往能在不显著增重的前提下改善面形。

最后需要强调,有限元模型是优化算法的“眼睛”。网格、边界、连接、材料和评价节点五者中任何一项失真,都会让优化结果偏离真实光机性能。建议在优化前先做一次基准工况校核:施加已知重力载荷,提取镜面节点位移,与理论估算或试验数据对比。若偏差在可接受范围内,再启动拓扑优化;否则应先修正模型,再谈轻量化。