2. 光机结构设计约束与轻量化边界

轻量化不是把材料“挖掉”那么简单。对光机结构而言，减重的每一步都在与刚度、强度、热稳定性、模态、面形和指向精度博弈。拓扑优化的价值，在于把这些工程语言翻译成数学语言——目标函数与约束条件，让算法在可行域内寻找最优材料分布。本节的核心任务，就是完成这场“翻译”。

一、从工程指标到数学约束的映射逻辑

一个可计算的优化问题通常写成：在满足约束 g_i(x) ≤ 0 的前提下，最小化目标 f(x)。对光机结构，设计变量 x 常取单元相对密度 ρ_e∈[0,1]，目标多为柔度最小（即刚度最大）或质量最小。关键在于约束如何建立。

刚度约束最直接：以柔度 C = U^TKU 表征，要求 C ≤ C₀，其中 C₀ 由许用位移反推。例如某反射镜支撑结构要求镜面在自重下最大位移不超过 50 nm，则先做一次基准分析得到参考柔度，再将其设为上限。强度约束则用应力判据：σ_max ≤ σ_s/n，安全系数 n 通常取 1.5～2.5，对脆性材料如微晶玻璃还需考虑拉压不对称。

二、热稳定性：从温度场到面形误差

热稳定性约束比力学约束更隐蔽。温度变化 ΔT 通过热膨胀系数 α 引起结构变形，进而传递到光学面形。对反射镜，面形误差常用 PV（峰谷值）和 RMS（均方根）描述，二者与 Zernike 系数直接相关。可计算约束可写为：在给定 ΔT 工况下，镜面节点法向位移的 RMS ≤ λ/20（λ 为工作波长）。

这里有一个实用技巧：把热载荷等效为热应变 ε_th = αΔT，在有限元中作为初应变施加，再提取镜面节点位移。若结构存在热致刚体位移（如支撑点不对称），还需分离刚体位移与面形畸变，否则 RMS 会被刚体项“污染”。

三、模态与指向精度约束

模态约束用于保证结构在动态环境下不产生共振。对空间光学载荷，一阶固有频率 f₁ 通常要求高于运载器激励频率的 2～3 倍，例如要求 f₁ ≥ 120 Hz。在优化中，这转化为特征值问题：Kφ = ω²Mφ，约束 ω₁² ≥ ω_min²。由于特征值对设计变量的灵敏度计算量较大，工程上常采用模态置信因子或近似代理模型来降低计算成本。

指向精度约束则关注光学元件之间的相对角位移。例如，主镜与次镜之间的光轴夹角变化 θ 需控制在 10 μrad 以内。这可以通过提取两个镜面中心节点的位移差，除以间距得到。若结构存在柔性铰链或可调机构，还需将刚体转动与弹性变形分开处理。

四、约束的归一化与尺度处理

多约束优化中，各约束量纲不同：位移是纳米级，频率是百赫兹级，应力是兆帕级。若直接代入算法，数值差异会导致灵敏度矩阵病态。标准做法是归一化：将每个约束除以其许用值，使所有约束在初始设计点附近量级为 1。例如：

g1 = C / C_allow - 1 ≤ 0
g2 = σ_max / (σ_s / n) - 1 ≤ 0
g3 = RMS / (λ/20) - 1 ≤ 0
g4 = f_min / f1 - 1 ≤ 0
g5 = θ_tilt / θ_allow - 1 ≤ 0

这样，算法对每个约束的“敏感程度”趋于一致，避免某一约束主导迭代方向。

五、轻量化边界的确定

轻量化边界不是单一数值，而是多约束共同围成的可行域边界。实际工作中，可先做单约束优化，得到各指标下的最小质量，再取最大值作为初始质量上限。例如：

• 仅刚度约束下，最小质量 3.2 kg；
• 仅模态约束下，最小质量 4.1 kg；
• 仅热稳定性约束下，最小质量 3.8 kg；
• 仅面形约束下，最小质量 4.5 kg。

则综合轻量化边界至少为 4.5 kg，且最终设计需同时满足全部约束。若某约束在优化中始终无法满足，说明该约束是“瓶颈”，需回到设计空间或材料选择层面调整，而非继续减重。

六、一个简化的约束转化示例

假设某反射镜支撑结构，材料为碳化硅（E = 410 GPa，α = 2.5×10^-6/K），镜面直径 200 mm，要求：自重下镜面 RMS ≤ 15 nm，ΔT = 20 K 时 RMS ≤ 25 nm，一阶频率 ≥ 150 Hz，总质量 ≤ 2.0 kg。转化后的约束为：

g1: RMS_gravity / 15nm - 1 ≤ 0
g2: RMS_thermal / 25nm - 1 ≤ 0
g3: 150Hz / f1 - 1 ≤ 0
g4: mass / 2.0kg - 1 ≤ 0

在有限元中，RMS 通过镜面节点法向位移的均方根计算：RMS = sqrt( (1/N) Σ δ_i² )。注意去除刚体位移分量，否则结果偏大。频率通过模态分析提取，质量通过单元体积乘以密度求和。所有约束归一化后代入 MMA 算法，即可开始迭代。