空间望远镜焦平面基板结构热-力耦合双目标拓扑优化设计
Dual-objective topology optimization design for focal plane substrate structure of space telescope with thermo-mechanical coupling
通讯作者:
收稿日期: 2025-05-07 修回日期: 2025-07-02
| 基金资助: |
|
Received: 2025-05-07 Revised: 2025-07-02
作者简介 About authors
曹雪虎(1997—),男,硕士生,从事结构轻量化研究,E-mail:
关键词:
Keywords:
本文引用格式
曹雪虎, 纪敬虎.
CAO Xuehu, JI Jinghu.
焦平面是空间望远镜的关键组件,其表面安装了大量光学元件,承担着信息转换的重要功能[3]。焦平面基板结构为这些光学元件提供稳定支撑,其稳定性直接影响最终的成像质量。因此,焦平面基板结构需满足多维度性能要求:首先,应具备卓越的刚度特性,既能维持光学元件的位形精度,又能确保发射阶段的结构基频满足动力学要求;其次,需具备有效的热控体系,可将温度梯度控制在较小范围内,避免热应力过大而引发变形;最后,必须实现轻量化设计,以降低发射成本[4]。传统的经验设计方法在小口径望远镜焦平面基板结构的设计中尚可适用,但对于大口径望远镜,由于结构惯性载荷急剧增长,常规设计方法难以兼顾刚度、热稳定性与轻量化的多维需求。因此,需要引入一种新的设计方法,以实现大口径望远镜焦平面基板结构的多目标协同优化设计。
拓扑优化理论起源于Michell的桁架理论[5],经过数十年的发展,先后演变出均匀化法[6]、变密度法[7]、水平集法[8]和渐进结构进化法[9]等方法。在结构设计初期应用拓扑优化,优化空间大,且更易获得形状合理、布局紧凑的轻量化结构。例如:刘伟等[10]分析了齿轮箱体的静动态特性并对其进行了多目标拓扑优化,改进后的箱体在质量减小11.2%的同时,最大应力下降了39.6%,位移减小了5.1%。韩博文等[11]提出了基于折衷规划法与熵权-层次分析的多目标拓扑优化方法,并将其应用于机载雷达支架的轻量化设计与结构优化,优化后支架的质量减小了53.8%。沈保山等[12]采用变密度拓扑优化法对转向摇臂支架进行了轻量化设计,减重比例高达54.8%。Zeng等[13]受铸造约束的拓扑优化方法的启发,提出了一种适用于大口径空间望远镜反射镜结构设计的优化模型,该模型可同时优化加强筋的分布和高度,优化后反射镜的光学性能提升了62.5%。Mao等[14]针对自主式水下航行器(autonomous underwater vehicle, AUV)的电池架在工作过程中面临的热-力耦合挑战(包括电池发热和载荷冲击),建立了热-固耦合拓扑优化模型,为AUV电池架的创新设计提供了理论支撑和优化方法。张小强等[15]针对聚变堆偏滤器的W/Cu模块进行了拓扑优化设计,在10 MW/m²热负荷下,W/Cu模块的最高温度降低了108.5 K,界面热应力下降了264.2 MPa,冷却性能显著提升。田立勇等[16]基于煤矿井下的复杂环境,通过对更换托辊机器人的履带式底盘进行动力学仿真和对其姿态调整机构进行拓扑优化,使横移平台的应力降低了13.71 MPa,质量减小了36.92%,从而显著提升了机器人的行驶稳定性与结构性能。
近年来,面向航空航天领域复杂热力环境下的多物理场耦合多目标拓扑优化研究不断深入,并取得了一系列重要进展。郭伟超等[17]提出了一种宏微观并行的拓扑优化方法,并结合聚类算法划分微结构类型,解决了热-力耦合载荷下导弹导引头的轻量化与热防护一体化设计问题,优化后导引头的质量减小了50%。叶帅等[18]结合拓扑优化与尺寸灵敏度分析,对星载临边光谱仪反射镜背部筋的构型进行了设计,并引入综合评价因子优化了柔性支撑结构,在重力载荷与热-力耦合工况下,实现了在保证反射镜面型精度满足使用要求的同时减重60%。朱兆基等[19]提出了一种静-动-随机振动多工况优化流程,并结合柔度约束和固有频率最大化目标,对火箭发动机的承力结构进行了优化,使其质量减小了30%,一阶模态频率提升了35.4%,进而增强了发动机对振动环境的适应性。
尽管上述研究在热-力耦合多目标优化方面取得了显著成果,并逐渐应用于航空航天结构,但当前的拓扑优化研究主要集中在单物理场下多工况或者双物理场下单目标的结构优化上,无法满足空间望远镜焦平面基板结构对高刚度、低温度梯度和轻量化的多维要求。为此,本文针对大口径空间望远镜焦平面基板结构,构建其在热-力耦合作用下的双目标拓扑优化模型,并以柔度最小化与温度梯度最小化为目标进行协同优化,致力于开发出能够适应复杂工况的高性能轻量化基板结构。
1 空间望远镜焦平面基板结构的设计要求
初始的大口径空间望远镜焦平面基板结构如图1所示,其安装方式为:四周由螺栓固定,光学元件通过夹紧装置、定位装置等安装在中间凹槽内。针对航空航天领域的轻量化和高刚度要求,基板结构的材料选择TC4钛合金。TC4钛合金的密度为4.43 g/cm3,弹性模量为114 GPa,导热系数为8 W/(m·K),屈服强度为860 MPa。
图1
图1
初始空间望远镜焦平面基板结构示意
Fig.1
Schematic of original substrate structure of space telescope focal plane
空间望远镜焦平面基板结构所受的载荷除了自身重力外,还包括安装在凹槽内的光学元件以及其余装置所产生的压力;热量主要由光学元件工作时所产生。所设计的基板结构需满足以下要求:1)最大变形量不超过0.05 mm;2)最大von Mises应力不超过所用材料的屈服极限和抗拉极限;3)温度梯度小于0.05 K/mm;4)基频不低于120 Hz。
2 基板结构热-力耦合双目标拓扑优化设计
2.1 拓扑优化方法
连续体结构拓扑优化通常基于有限元离散化方法,通过引入材料分布表征函数来确定单元域的材料分布状态。其中:函数值为1表示实体材料,0表示空隙区域。现阶段,最常用的拓扑优化方法为固体各向同性材料惩罚模型(solid isotropic material with penalization, SIMP)。在SIMP中,将每个单元的材料视为多孔材料,以相对密度作为设计变量,并将单元的材料特性(如弹性模量)表示为相对密度的函数,数学表达式如下:
式中:Ei 为单元i中材料的弹性模量;
由
2.2 热-力耦合双目标拓扑优化数学模型
本文针对空间望远镜焦平面基板结构的热-力耦合拓扑优化问题,构建了综合考虑柔度最小化与温度梯度最小化的双目标优化模型。结合结构力学与热传导理论,建立如下目标函数:
式中:П为双目标函数;φ为基板结构的柔度;ψ为基板结构的温度梯度;w为目标的权重系数;φ0和ψ0分别为初始柔度和初始温度梯度,引入初始参数的目的是实现目标函数的无量纲化。
基板结构柔度的目标函数可表示为:
式中: U 为位移矩阵; K 为刚度矩阵; F 为外载荷,包括机械载荷 Fm和热载荷 Ft。
对于基板结构的温度梯度,为了量化结构的导热均匀性,以温度梯度的积分作为目标函数,表示为:
其中:
式中: Kt为热刚度矩阵,
由于基板结构的材料是均匀的,因此可使用体积约束替代质量约束,体积约束定义如下:
式中:V为设计域的体积,f为优化体积比,V0为整个设计域的初始体积,vi 为优化后的单元体积。
结合
式中:q为热源强度,
2.3 双目标权重系数
在解决热-力耦合双目标拓扑优化问题时,目标权重系数的确定是关键环节。针对热力学与结构力学双场耦合的复杂场景,优化问题本质上属于典型的多属性决策问题,各目标函数之间既存在涉及能量传递的物理耦合关系,又表现出不同量纲与数量级的数学特性差异。为此,本文通过引入层次分析法来构建系统化的权重决策模型[20]。
在层次分析法中,层次结构模型的科学构建是确定权重系数的基础。本文基于热-力耦合双目标拓扑优化的特性,构建了聚焦于结构柔度和温度梯度这2个关键因素的层次结构模型。
在完成层次结构模型构建后,着手构造判断矩阵。判断矩阵的构造是基于各因素间的相对重要性来实现的。在比较各因素重要性的过程中,采用1~9标度法进行赋值操作,具体赋值规则详见表1。
表1 标度法的赋值规则
Table 1
| 相对重要性① | 赋值 |
|---|---|
| 同等重要 | 1 |
| 稍微重要 | 3 |
| 明显重要 | 5 |
| 强烈重要 | 7 |
| 极端重要 | 9 |
① 相对重要性是指因素j相对于因素t的重要程度。
标度法为衡量各因素之间的相对重要性提供了一套标准化、可量化的方法。当对因素j和因素t进行比较时,若认为因素j相对于因素t的重要程度为a(a的取值依据1~9标度法确定),则判断矩阵中因素j与因素t对应的位置赋值a;同时,由于比较是相互的,则判断矩阵中因素t与因素j对应的位置赋值1/a。基于上述规则,对所有因素进行两两比较并完成赋值操作,即可构造完整的判断矩阵。
对于本文所研究的焦平面基板结构,柔度占主导地位,故柔度因素相对于温度梯度因素为极端重要,由此可得到表2所示的判断矩阵。
在层次分析法的理论架构内,当决策系统的维度n=2时,其判断矩阵具有内在的完全一致性。凭借相关的数学推导,可知一致性比率CR=0恒成立,这是由正互反矩阵的特定结构属性决定的。设二维决策系统的判断矩阵 A 为:
判断矩阵 A 的特征方程可表示为:
式中:λ为特征值。
根据
定义一致性指数CI为:
由
通过几何平均法计算优化目标的权重系数,计算式如下:
将极端重要对应的赋值代入
3 基板结构拓扑优化与有限元分析
3.1 优化流程
空间望远镜焦平面基板结构热-力耦合双目标拓扑优化设计流程遵循系统化实施路径:在预分析阶段,利用COMSOL Multiphysics软件建立热-力耦合有限元模型并定义设计域边界;随后,分别进行柔度最小化与温度梯度最小化的单目标拓扑优化,以获取各自的最优值作为后续双目标函数归一化处理的基准。接着,通过层次分析法确定柔度和温度梯度的权重系数,并采用SIMP变密度法进行迭代求解,以获得最优的材料分布;最后,对基板结构的拓扑构型进行几何重构,以提升后续的工艺可行性,并通过有限元仿真来验证其刚度、强度及温度均匀性等核心指标。具体优化流程如图2所示。
图2
图2
基板结构热-力耦合双目标拓扑优化流程
Fig.2
Dual-objective topology optimization flow of substrate structure with thermo-mechanical coupling
3.2 优化结果与几何重构
在开展热-力耦合双目标拓扑优化时,鉴于目标函数已进行归一化处理,为有效防止优化过程过早达到收敛状态,将收敛容差设为0.000 01,其余各项参数的设置均与单目标拓扑优化时保持一致。
图3所示为基板结构拓扑优化目标函数随迭代次数的变化曲线。结果表明:在初始优化阶段(0~20次迭代),目标函数值呈急剧下降的趋势,反映出基板结构具有显著的可优化空间;从第20次迭代开始一直到第60次迭代期间,目标函数值虽仍持续下降,但下降速度明显减缓,这意味着优化过程正在逐步逼近最优解,相应的优化空间也逐渐变小;在第60次迭代后,目标函数值基本上保持不变,趋于稳定状态,这表明已十分接近最优解。
图3
图3
基板结构拓扑优化目标函数的变化曲线
Fig.3
Variation curve of objective function for substrate structure topology optimization
图4
图4
优化后基板结构多切面单元密度分布云图
Fig.4
Cloud map of multi-slice element density distribution of substrate structure after optimization
综上,所构建的热-力耦合双目标拓扑优化模型具有优良的收敛特性与稳定性,收敛过程符合优化问题的典型特征,验证了该拓扑优化方法的工程适用性,实现了预期优化目标。
图5
图5
双目标优化重构后的基板结构
Fig.5
Substrate structure after dual-objective optimization and reconstruction
3.3 有限元仿真验证
图6
图6
优化前后基板结构的位移分布云图
Fig.6
Cloud maps of displacement distribution of substrate structure before and after optimization
图7
图7
优化前后基板结构的von Mises应力分布云图
Fig.7
Cloud maps of von Mises stress distribution of substrate structure before and after optimization
图8
图8
优化前后基板结构的热应变分布云图
Fig.8
Cloud maps of thermal strain distribution of substrate structure before and after optimization
根据有限元仿真分析结果,导出优化前后基板结构在热-力耦合作用下各个性能指标的最大值,如表3所示。
表3 优化前后基板结构的性能对比
Table 3
| 对比项 | 最大位移/10-3 mm | 最大应力/MPa | 最大热应变/10-5 | 最大温度梯度/ (K/mm) | 基频/Hz | 质量/kg |
|---|---|---|---|---|---|---|
| 变化幅度/% | -4.92 | +367.89 | +10.22 | +22.67 | -21.54 | -39.94 |
| 优化前 | 1.22 | 21.8 | 1.86 | 0.034 4 | 1 021.1 | 3.220 |
| 优化后 | 1.16 | 102.0 | 2.05 | 0.042 2 | 801.2 | 1.934 |
为进一步量化分析优化前后基板结构的性能,设定以下阈值:以最大位移的1/2为高位移阈值,最大应力的1/2为高应力阈值,最大热应变的7/10为高热应变阈值。通过计算可得,相比于优化前,优化后基板结构的高位移区域占比从27.0%降至19.5%,高应力区域占比从0.003 0%降至0.000 5%,高热应变区域占比从87.7%降至76.9%。定量分析结果表明,优化后基板结构的高值区域显著减小,且材料分布更均匀,这表明拓扑优化技术的应用显著改变了结构的力传递路径和散热路径,同时使变形量分布和应力分布更加合理,从而提高了结构的整体性能。
优化后基板结构的基频为801.2 Hz,相较于优化前的1 021.1 Hz有所降低,但远高于设计要求的120.0 Hz,且与其他组件的基频相差较大,可有效避免因共振而导致的结构破坏。
综上,优化后的基板结构满足设计要求,同时其质量减小了39.94%,实现了轻量化设计。采用3D打印技术制造优化后的新基板结构,如图9所示。
图9
4 总 结
本文对大口径空间望远镜焦平面基板结构进行了热-力耦合双目标拓扑优化,并通过有限元仿真对优化结果进行了验证,得到以下结论。
1)针对刚度、热稳定性与轻量化的协同设计难题,提出了一种热-力耦合双目标拓扑优化模型,并将其成功应用于空间望远镜焦平面基板结构的优化设计,获得了满足设计要求的拓扑构型。
2)根据拓扑优化结果,考虑加工性和装配性,对基板结构进行了几何重构,并对重构后的基板结构进行了有限元仿真。结果表明:优化后基板结构的质量为1.934 kg,相比于优化前减小了39.94%;最大位移为1.16×10-3 mm,相比于优化前降低了4.92%;最大应力与最大热应变均满足设计要求。新基板结构在刚度性能、热稳定性满足设计要求的前提下,质量大幅减小,实现了轻量化目标。
3)通过3D打印技术制作了优化后的基板结构,验证了热-力耦合双目标拓扑优化设计方法的工程可行性,这可为其他类似航空航天结构的轻量化设计提供参考。
参考文献
Questions of constructing lightened primary mirrors of space telescopes
[J].
Space telescope design considerations
[J].
Hubble space telescope: servicing mission 3A media reference guide
[EB/OL]. [
The limits of economy of material in frame-structures
[J].
Generating optimal topologies in structural design using a homogenization method
[J].
Design of material structures using topology optimization
[D].
Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations
[J].
A simple evolutionary procedure for structural optimization
[J].
齿轮箱体静动态特性分析与多目标拓扑优化
[J].
Static and dynamic characteristic analysis and multi-objective topology optimization of gearbox
[J].
某机载雷达支架轻量化与结构拓扑优化研究
[J].
Research on lightweight and structural topology optimization of an airborne radar support
[J].
双前桥液压助力转向摇臂支架轻量化设计
[J].
Lightweight design of double front axle hydraulic power steering rocker arm bracket
[J].
Optimal design of rectangular mirror based on topology and size optimization
[C]//
Design and analysis of the thermal-stress coupled topology optimization of the battery rack in an AUV
[J].
聚变堆偏滤器拓扑优化设计与稳态热分析
[J].
Topology optimization design and steady-state thermal analysis of fusion reactor divertor
[J].
更换托辊机器人履带式底盘的仿真与优化
[J].
Simulation and optimization of crawler chassis of idler replacement robot
[J].
热力耦合情况下的导弹导引头多尺度并行拓扑优化设计方法
[J].
Multi-scale parallel topology optimization design method for missile seeker with thermo-dynamic coupling loads
[J].
星载临边光谱仪反射镜组件优化设计分析
[J].
Optimization design of the spaceborne edge spectral imager reflective mirror assembly
[J].
液体火箭发动机典型承力结构动态拓扑优化设计
[J].
Dynamic topology optimization design of typical load-bearing structures in liquid rocket engines
[J].
基于加权系数分配的火炮座圈优化设计
[J].
Optimal design of artillery seat based on weighting coefficient distribution
[J].
汽车悬架下节叉多工况结构拓扑优化设计
[J].
Structural topology optimization design of automotive suspension lower knuckle fork under multiple conditions
[J].
/
| 〈 |
|
〉 |


