双面共体反射镜结构优化设计与分析

2024-02-18谭双龙马思宇吴洪波

应用光学 2024年1期

马霖，谭双龙，闫磊，马思宇，吴洪波，张新

（1.中国科学院长春光学精密机械与物理研究所，吉林长春 130033；2.应用光学国家重点实验室，中国科学院光学系统先进制造技术重点实验室，吉林长春 130033）

引言

红外成像具有全天候工作、穿透力强、抗干扰能力强、受天气影响小、制导精度高等优点，在侦察探测、导引等领域发挥了重要作用[1-2]。反射镜作为红外光学系统中的重要部件，面形稳定性、轻量化率是衡量其性能的关键指标。通过拓扑优化技术可以在保证刚度最大化的同时，提高反射镜的轻量化率，国内外学者也对此做出了诸多研究。在较早的时候，国外Park 等以自重、抛光压力作为载荷条件，获得高轻量化的镜体为目标，利用拓扑优化技术优化单面反射镜背部的材料分布形式，并加工制造反射镜体进行验证[3]。David K 通过尺寸优化的方式以最大刚度为目标，优化镜面厚度、筋板厚度等参数，得到了一种背板开放的高轻量化铍反射镜[4]。韩国LEE 等人通过尺寸优化的方式，对某大口径镜体背部肋的厚度和安装位置进行优化，使其应力应变分布更加均匀[5]。如图1 所示，从左到右分别为Park，David，LEE 优化所获的单面反射镜。

图1 国外反射镜优化示例Fig.1 Examples of mirror optimization overseas

国内方面，大连理工大学Liu Shutian 等人通过拓扑优化技术，以最大刚度为目标，对大型空间反射镜背部的结构形式进行尺寸参数优化和形状优化，完成高刚度高轻量化反射镜的优化设计[6]。长春光机所Liu Guang 等人将拓扑优化结合尺寸优化，利用连续体拓扑优化技术获得了材料的最优分布，根据拓扑结果对轻质肋分组进行尺寸优化，得到最优解[7]。如图2 所示，从左到右依次为Liu Shutian，Liu Guang 优化的单面反射镜。

图2 国内反射镜优化示例Fig.2 Examples of mirror optimization in domestic

此外，西安光机所、安徽光机所的研究人员也基于拓扑优化方法对反射镜展开很多研究[8-10]。但上述研究均为针对单面反射镜的优化工作。随着光学系统产品高集成化、轻量化的发展，多镜面共体化、集成化方案的紧凑优势逐渐凸显出来，受到了越来越多研究人员的关注[11-13]。相比于单面反射镜，双面共体反射镜集成度更高，更适合应用在对重量要求苛刻的轻巧红外光学系统中。但双镜面对支撑刚度要求更高，如果刚度不足，2 个反射镜面变形相互耦合，对整个光学系统的像质影响会叠加放大，故设计难度更高。需确保2 个反射镜面自身刚度足够大且相互之间影响最小，同时为了满足使用需求，还需要有足够的轻量化率，这就极大地增加了设计难度。基于上文所述国内外研究人员对单面反射镜的优化方法，将其应用在双面共体反射镜的设计中是否可行，目前此方面的研究甚少。在此基础上，本文结合拓扑优化和尺寸优化方法，对某同轴四反红外光学系统中主四镜共体结构展开优化设计，并对最终设计方案进行分析和测试，验证优化设计方法的有效性。

1 优化方法

1.1 拓扑优化

由于结构的多样性，传统设计优化时需要大量的迭代，而拓扑优化可以设定载荷和约束条件，给出既定目标和优化变量，通过算法计算出设计域中相对最优布局，给设计提供参考和支撑，极大地提高了设计效率。拓扑优化流程如图3 所示。

图3 拓扑优化设计流程Fig.3 Flow chart of topology optimization design

拓扑优化中常用变密度法，它源于微观结构等效的思路，以元素的相对密度为设计变量，通过密度从0 到1 的空间配置变化来实现拓扑结构。0 表示该区域无密度，存在的必要性低，应尽量优化掉该部分；1 表示该区域存在的必要性大，需要保留[14]。反射镜结构拓扑优化过程中，通常关注反射镜的面形精度变化RMS 最小，而现有的拓扑优化商业软件不能直接使用RMS 作为优化目标[15]。为了最小化RMS，需要保证整个反射镜镜体具有足够大的结构刚度来抵抗变形，避免外载荷影响下面形精度变化过大，故本文以刚度最大即柔度最小作为目标函数，约束整体镜体体积分数，保证轻量化设计条件下，设计刚度最大的双面共体反射镜镜体。数学模型如下：

式中：ρi为相对密度；U为全局位移矢量；K为全局刚度矩阵；F为全局负载向量；V为结构的总体积；υi为单位体积；N为变量数[16]。

1.2 尺寸优化

通过拓扑优化可以确定材料的合理分布形式，但最优分布形势下的最佳尺寸还需要进一步确认，这就用到尺寸优化[17]。尺寸优化是在形状不变的基础上，对参数进行优化，其设计变量可以定义为壳单元厚度、梁单元截面积等。本文以最小柔度为目标，一阶频率和筋板质量为约束条件，寻找最佳尺寸参数。数学模型如下：

式中：ti为筋板厚度；f1为一阶频率；为设定的约束频率；Mass为求解域质量；为设定的约束质量。

2 优化设计

2.1 光学系统参数

传统同轴四反光学系统设有四块反射镜，本文采用主、四共体的思路，使布局更加紧凑，但对主、四共体反射镜镜体的设计提出了更高的要求。系统光路如图4 所示。红外系统方案中，反射镜口径为186 mm，工作波段8 μm～12 μm，设计要求总质量低于1 kg，主、四镜面形RMS 要求低于λ/5（λ=632.8 nm），自重对主、四镜面形RMS 影响要求低于λ/20，一阶频率要求大于400 Hz。

图4 同轴四反光路图Fig.4 Optical path diagram of coaxial four-mirror system

2.2 材料选择

随着超精密加工技术的突飞猛进[18]，铝合金材料因重量较轻、加工性能优良，广泛应用在光学领域。铝合金既可以做结构件，也可以做光学件，可以设计成光、机同种材料，更好地实现整机被动无热化，且制造成本相对玻璃陶瓷更低[19-21]。本文主四镜材料选择铝合金6061-T6，材料属性见表1 所示。

表1 铝合金6061-T6 属性Table 1 Properties of aluminum alloy 6061-T6

2.3 初始模型设计

首先需要确定支撑形式，反射镜支撑种类一般包括中心支撑、背部支撑、侧面支撑、周边支撑四种方式。根据平板理论[22]，反射镜的支撑点数目可由下式计算：

式中：R为反射镜半口径，单位为in；E为材料的弹性模量，单位为lb/in2；δ为允许的PV 值，单位为in，为保证裕度充足取δ=λ/20（λ=632.8 nm）；t为反射镜厚度，单位为in；ρ为材料密度，单位kg/m2。主镜口径大于四镜，将表2 所示主镜各参数代入公式（1）进行计算。

表2 反射镜结构参数Table 2 Structural parameters of mirror mm

计算得到N=1.54，表明采用3 点支撑形式已足够保证镜体的面形不受支撑点数量多少的影响。为了避免过约束带来的非必要的镜面变形，不再引入过多支撑，最终采用3 点支撑的形式。

本文拓扑优化设计的双面共体反射镜有限元模型如图5 所示。其中蓝色部分为非设计区域，分别包括3 个支撑耳，主镜面和四镜面；黄色区域代表设计区域，F1、F4 分别代表主、四镜上模拟加工抛光压力设置的负载。优化目标时，选择主、四镜面柔度最小为优化目标。柔度是变形与负载的比值，反射镜柔度最小可以等效负载条件下的面形变化最小[23]。约束设计区域时优化体积分数在0.3 以内，确保足够的轻量化率。固定安装耳的圆孔，按水平放置的自重载荷和加工抛光压力设置负载。本文双面共体反射镜无表面改性及镀膜工序，可以被动消热且无需考虑双金属效应等影响，设置去除材料的方向关于光轴对称。

图5 拓扑优化有限元初始模型Fig.5 Initial finite element model of topology optimization

2.4 拓扑优化

经过51 轮迭代，计算结果逐渐收敛到一个较稳定的状态。拓扑优化后隐去主镜面非设计区域及隐去四镜面非设计区域的材料密度分布图如图6 所示。图6 中红色表示材料密度为1，即对刚度贡献最大的地方（主、四镜镜面区域为非设计区域，也显示为红色）；蓝色表示材料密度为0，为可以去除的材料。

图6 拓扑材料密度分布图Fig.6 Density distribution diagram of materials by topology

为了便于观察，隐掉材料密度小于0.6 kg/m2的部分，如图7 所示。拓扑优化结果表明，在上述设置条件下，主镜面与四镜面通过以光轴方向为中心呈辐射形排列的筋板固定时刚度最高。

图7 材料密度大于0.6 kg/m2 分布结果Fig.7 Distribution results of density greater than 0.6 kg/m2

2.5 尺寸优化

根据最佳材料分布密度，考虑实际加工情况，对筋板分布进行微调，建立主、四镜及筋板面网格模型，如图8 所示，镜面和筋板均赋予为壳单元属性。尺寸优化方式上，有尺寸参数化优化和自由尺寸优化供选择。参数优化结果为等厚结果，本文中需要对2 个镜面厚度尺寸优化，参数化优化对镜面尺寸限制较大，采用自由尺寸优化可以获得镜面不同半径位置处的厚度，故采用自由尺寸优化，将筋板、主镜面、四镜面设为3 个独立的设计域。在分析软件中将主镜面柔度变量和四镜面柔度变量通过光学分配的权重系数整合成一个变量，以加权后的变量柔度最小为优化目标，负载设置同上述2.3 节，约束整体一阶频率及筋板的总质量。

图8 尺寸优化网格模型Fig.8 Mesh model for size optimization

两镜面和筋板自由尺寸优化后的结果如图9所示。为确保足够的刚度，在镜面内孔与筋板交界处，9 根筋板中，与安装耳（绿色）连接的尺寸值较大，为4.5 mm，其余筋板为1.5 mm，2 个镜面靠近内径部分厚度为7.5 mm，边缘为4 mm。

图9 自由尺寸优化结果Fig.9 Results of free size optimization

参考优化结果，综合考虑镜体刚度、轻量化率、安装适应性、传统加工铣削进给刀方向、抛光镀膜等加工工艺制造过程，对部分结构进行细化，最终设计的主四镜模型如图10 所示。

图10 最终设计的主四镜模型Fig.10 Final design of primary-fourth mirror model

最终模型3 个支撑耳分别连接2 个筋板，每2 个支撑耳中间再设置1 个筋板，共9 个筋板。该模型原始镜坯质量为4.53 kg，轻量化后为0.8 kg，符合质量低于1 kg 的要求，且双面反射镜的轻量化率达到82.4%。

3 设计结果分析

3.1 模态分析

对最终设计的主四镜模态进行分析，考量所设计反射镜的刚度特性，模态分析结果如表3 所示，对应的前三阶振型如图11～图13 所示。X轴、Y轴、Z轴方向见图11 所示。

表3 主四反射镜前3 阶模态Table 3 The first three modes of primary-fourth mirror

图11 第一阶振型Fig.11 The first-order vibration mode

图12 第二阶振型Fig.12 The second-order vibration mode

图13 第三阶振型Fig.13 The third-order vibration mode

模态分析结果表明，镜体一阶频率已达到417 Hz，具有足够高的刚度，符合大于400 Hz 的要求。振型主要以安装耳的变形为主，主四镜镜面及筋板构成的整体刚度更高，更加稳定，可以更好地保证面形精度要求。

3.2 自重变形分析

为了进一步验证反射镜的刚度，分析径向1g重力、轴向1g 重力做用下反射镜的面形变化。在有限元分析结果中导出反射镜面形变化前后的各节点坐标进行面形拟合，主、四镜面形云图如图14所示。

图14 主、四镜面形云图Fig.14 Cloud atlas of surface shape for primary-fourth mirror

主、四镜在X，Y，Z3 个方向1g 重力下的RMS值和PV 值结果如表4 所示。从表4 可以看出，面形RMS 变化最大发生在主镜Z方向上，换算成波长约为λ/22（λ=632.8 nm），符合低于λ/20 的要求，其余方向自重对主、四镜面形RMS 的影响均低于λ/37，说明其抵抗重力变形的能力出色，侧面证明其刚度充足。

表4 主、四镜1g 重力变形Table 4 Primary-fourth mirror deformation under 1 g gravity nm

4 面形测试

在完成优化设计与分析验证后，采用最终设计方案对镜体进行加工，在粗加工、半精加工、精加工后对镜体进行时效处理，消除镜面的残余应力，以确保面形的稳定。采用干涉法中的零位检测法测试主四镜面形RMS 是否满足要求。零位检测法通过计算机生成的全息图（computer-generated holograms，CGH），将干涉仪发出的球面波转化为主四镜面形所对应的非球面波，从而消除非球面的法线像差，可将测试精度提高至nm 级[24]。面形测试过程如图15 所示。