《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于改进粒子群算法的红外灯阵自动优化设计

张筱娴 林士峰 张旭 林宝军 蒋桂忠

张筱娴, 林士峰, 张旭, 等. 基于改进粒子群算法的红外灯阵自动优化设计[J]. 航天器环境工程, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
引用本文: 张筱娴, 林士峰, 张旭, 等. 基于改进粒子群算法的红外灯阵自动优化设计[J]. 航天器环境工程, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
ZHANG X X, LIN S F, ZHANG X, et al. Automatic optimization design of infrared lamp array based on enhanced particle swarm algorithm[J]. Spacecraft Environment Engineering, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
Citation: ZHANG X X, LIN S F, ZHANG X, et al. Automatic optimization design of infrared lamp array based on enhanced particle swarm algorithm[J]. Spacecraft Environment Engineering, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009

基于改进粒子群算法的红外灯阵自动优化设计

doi: 10.12126/see.2022.05.009
详细信息
    作者简介:

    张筱娴(1992—),女,硕士学位,研究方向为航天器热控设计及在轨热防护。E-mail: zhxxxy@126.com

  • 中图分类号: TP301.6; V416.5

Automatic optimization design of infrared lamp array based on enhanced particle swarm algorithm

图(17)
计量
  • PDF下载量:  8
  • 文章访问数:  73
  • HTML全文浏览量:  6
文章相关
  • 中图分类号:  TP301.6; V416.5
  • 收稿日期:  2022-05-05
  • 修回日期:  2022-09-18
  • 刊出日期:  2022-10-26

基于改进粒子群算法的红外灯阵自动优化设计

doi: 10.12126/see.2022.05.009
    作者简介:

    张筱娴(1992—),女,硕士学位,研究方向为航天器热控设计及在轨热防护。E-mail: zhxxxy@126.com

  • 中图分类号: TP301.6; V416.5

摘要: 为提高航天器热试验用红外灯阵的设计效率、提升灯阵热流分布均匀性,提出一种基于粒子群算法的红外灯阵自动优化设计方法:建立红外灯单灯的热流分布数据库与灯阵位置参数的自动优化设计方法;并通过优化参数的重新选择和优化参数间约束条件的建立,改进传统的粒子群算法以提高设计效率。计算结果表明,相比遗传算法,使用改进粒子群算法进行灯阵优化设计用时缩减82%,而优化所得的热流分布均匀度仅低了0.77%。

English Abstract

张筱娴, 林士峰, 张旭, 等. 基于改进粒子群算法的红外灯阵自动优化设计[J]. 航天器环境工程, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
引用本文: 张筱娴, 林士峰, 张旭, 等. 基于改进粒子群算法的红外灯阵自动优化设计[J]. 航天器环境工程, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
ZHANG X X, LIN S F, ZHANG X, et al. Automatic optimization design of infrared lamp array based on enhanced particle swarm algorithm[J]. Spacecraft Environment Engineering, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
Citation: ZHANG X X, LIN S F, ZHANG X, et al. Automatic optimization design of infrared lamp array based on enhanced particle swarm algorithm[J]. Spacecraft Environment Engineering, 2022, 39(5): 501-508 doi:  10.12126/see.2022.05.009
    • 红外灯是常用的航天器外热流模拟装置,因其使用灵活,对航天器表面遮挡较小,且能够重复使用,被广泛应用于航天器真空热试验中。对于红外灯阵的优化设计,国外航天机构的相关研究工作较早就有报道:1986年,Sanger[1]研究了一种新型的红外灯,并计算得到其较精确的热流分布;1994年,Turner等[2]采用蒙特卡罗方法计算红外灯的热流分布,并以此作为红外灯阵优化的基础;2004年,Ziemke[3]使用有限元法分析了红外灯的热流分布并优化灯阵设计。

      国内这方面的研究起步较晚:2010年之前,主要采用试验方法测试红外灯阵热流的不均匀度;2010年,张坤等[4]对带抛物面型反射器的红外灯泡进行了辐射仿真分析;2010年—2011年,北京卫星环境工程研究所[5-8]对管状红外灯进行了单灯热流仿真,并使用遗传算法优化了灯阵设计,此后又提出曲面红外灯阵热流仿真方法,结合遗传算法和复合形算法进行灯阵优化设计[9],并在单灯热流仿真的基础上分析影响灯阵热流均匀度的参数,以指导灯阵的优化设计[10];2012年之后,国内对提高红外灯阵热流均匀度的研究较少,鲜有的研究[11-13]仅基于实验和数据拟合,误差较大,通用性也有限。

      单纯使用数学计算方法对已有的红外灯阵进行优化设计,需要输入初始的阵列参数,且遗传算法虽然擅长找到全局最优解,但收敛较慢、运行时间长,在实际工程应用上不具优势。目前国内各大科研院所在进行航天器热试验时,仍以传统的灯阵设计手段为主,即先根据经验设计灯阵的间距和高度,再根据测得的热流均匀度来调整灯阵布局,不仅费时费力,且热流均匀度较差,影响真空热试验结果的准确性,进而影响航天器的热设计效率。

      为能快速地针对不同尺寸的卫星散热面设计出适用的灯阵布局,本文提出一种基于改进粒子群算法的红外灯阵自动优化设计方法,结合工程经验对优化参数进行约束,以求在提高收敛速度的同时能较准确地获得全局最优解,更加适用于红外灯阵优化设计的工程应用。

    • 红外灯阵自动优化设计的基础是较为准确地计算掌握红外灯阵的热流分布,由于各类商业热分析软件难以分析灯管对灯丝不同波段能量的折射、反射等情况,所以当前设计中普遍采用蒙特卡罗法计算红外灯的辐射热流。蒙特卡罗法通过追踪光线轨迹来获得热流分布,能够有效解决复杂的热流计算问题,但获得精确结果的代价是较长的计算时间。为了简化计算,通常先计算不同安装高度的单灯热流分布,获得红外单灯热流分布数据库,然后忽略各个红外灯之间的互相影响,通过坐标变换、叠加网格热流来计算红外灯阵的热流分布。

    • 以国内常用的飞利浦13169X/98型号红外灯为例,灯丝为钨丝,发射率接近1,绕制成单螺旋圆柱体样式,可等效为一个圆柱体。由于红外灯的灯管是封闭体,灯丝发出的能量一部分被灯管吸收,一部分由灯管透射出去,其余经反射等过程最终被灯管吸收或透射。

      灯管的材质为石英玻璃,对波长0.2~4.0 μm的粒子只透射或反射,对此波段外的粒子只吸收或反射,即有3.8 μm宽的透射谱带。代表能量的粒子在灯管上的反射/透射/吸收与粒子的波长和入射方向均有关。光谱反射率的计算方法[2]可表示为

      $$\left\{ \begin{gathered} {\rho _{\text{n}}}\left( \theta \right) = \frac{{{a^2} + {b^2} - 2a\cos \theta + {{\cos }^2}\theta }}{{{a^2} + {b^2} + 2a\cos \theta + {{\cos }^2}\theta }} \\ {\rho _{\text{p}}}\left( \theta \right) = {\rho _{\text{n}}}\left( \theta \right)\frac{{{a^2} + {b^2} - 2a\sin \theta \tan \theta + {{\sin }^2}\theta {{\tan }^2}\theta }}{{{a^2} + {b^2} + 2a\sin \theta \tan \theta + {{\sin }^2}\theta {{\tan }^2}\theta }} \\ \end{gathered} \right. \text{,} $$ (1)

      其中,

      $$\left\{ \begin{gathered} 2{a^2} = {\left[ {{{\left( {{\eta ^2} - {\kappa ^2} - {{\sin }^2}\theta } \right)}^2} + 4{\eta ^2}{\kappa ^2}} \right]^{1/2}} + \left( {{\eta ^2} - {\kappa ^2} - {{\sin }^2}\theta } \right) \\ 2{b^2} = {\left[ {{{\left( {{\eta ^2} - {\kappa ^2} - {{\sin }^2}\theta } \right)}^2} + 4{\eta ^2}{\kappa ^2}} \right]^{1/2}} - \left( {{\eta ^2} - {\kappa ^2} - {{\sin }^2}\theta } \right) \\ \end{gathered} \right.。 $$ (2)

      式(1)和式(2)中:ρn(θ)和ρp(θ)分别为反射率的法向分量和平行分量;θ为粒子的入射角;ηκ分别为石英玻璃的折射率和吸收系数,可查表得到。

      灯管吸收部分灯丝的能量后温度升高,成为不可忽略的二次辐射源,可视其为灰体,发射率取0.94。叠加由灯丝和灯管发出的辐射热流,可以计算出被照平面上的热流分布,求解流程如图1所示。

      图  1  单灯热流分布求解流程

      Figure 1.  Solution process for calculating heat flux distributions for a single lamp

      根据工程经验,红外灯的安装高度通常在0.3~0.6 m之间,故本文计算此安装高度区间内的红外单灯的热流分布,计算间隔为0.1 m,xy方向分别垂直和平行于红外灯的轴线方向,h为灯的安装高度,图2给出部分安装高度的单灯热流分布。

      图  2  不同安装高度的单灯热流分布

      Figure 2.  Heat flux distributions for a single lamp at different installation heights

    • 图3为某卫星热试验用红外灯阵。灯阵平行于卫星散热面排布,同一平面的红外灯可安装在相同或不同高度,为了提高热流均匀性,一般在灯阵的周围加设反光板。计算灯阵热流时需考虑多个红外灯直射热流与四周反光板反射热流的叠加。

      图  3  某卫星热试验用红外灯阵

      Figure 3.  The infrared lamp array for a satellite’s thermal test

      图4所示,被照面上的每一个网格点,不仅接收来自每只红外灯的直接辐射热流,还接收来自每块反射板对每只灯的反射热流。因此,被照面上任一网格点的热流计算式为

      图  4  被照面热流示意

      Figure 4.  The heat flux on the illuminated surface

      $$ q = \sum\limits_{i = 1}^N {{{\left( {{q_i}} \right)}_{\rm{dir}}}} + \sum\limits_{i = 1}^N {\sum\limits_{j = 1}^M {{{\left( {{q_{i j}}} \right)}_{\rm{ref}}}} } \text{,} $$ (3)

      式中:q为网格点的总热流;qdirqref分别为红外灯的直接热流和反射热流;N为红外灯的总数量;M为反光板的总数量。为简化示意,图4只展示了1块反光板的反射热流和4只红外灯的直接热流,其中反光板对红外灯的反射可视为镜面反射。

      工程上通常使用热流不均匀度来衡量被照面红外灯阵的热流分布均匀性[14],其计算公式为

      $$ \delta =\frac{{q}_{\mathrm{m}\mathrm{a}\mathrm{x}}-{q}_{\mathrm{m}\mathrm{i}\mathrm{n}}}{{q}_{\mathrm{m}\mathrm{a}\mathrm{x}}+{q}_{\mathrm{m}\mathrm{i}\mathrm{n}}}\times 100\% \text{,} $$ (4)

      其中,qmaxqmin分别为红外灯阵在被照面的最大和最小热流密度。

    • 1)单灯热流计算验证

      计算输入功率为71 W、红外灯高度为0.3 m条件下的单灯热流分布,计算结果由灯丝热流和灯管热流叠加而得。由于灯管的热流计算不涉及石英玻璃的光谱特性,使用商业软件同样能够得到较为精确的结果,所以将本文程序计算得到的灯管辐射热流与商业软件仿真得到的结果进行比对验证,如图5所示,可以看出两者一致性良好,证明本文计算方法基本无误。

      图  5  单灯灯管热流的本文计算结果与商业软件仿真结果对比

      Figure 5.  Comparison of the heat flux for a single lamp tube between the calculation results obtained with this paper’s method and the simulation results obtained with commercial software

      以本文程序计算方法,采用文献[2]中的红外灯参数,计算输入功率1032 W、安装高度0.152 4 m条件下的单灯热流分布;并与文献[2]的计算和实测结果进行比对验证,如图6所示,可以看出,本文方法与文献[2]的数据在中心区域一致性良好,证明本文计算方法基本无误。

      图  6  单灯热流的本文计算结果与文献数据对比

      Figure 6.  Comparison of the heat flux for a single lamp between the calculation results obtained with this paper’s method and the reference data

      2)灯阵计算验证

      为了验证本文灯阵热流分布计算的正确性,借用文献[7]中的灯阵尺寸进行计算,并与文献[7]的灯阵热流数据进行比对验证。灯阵为7行×6列均匀分布,每只灯安装高度均为0.45 m,本文计算所得灯阵热流分布如图7所示,不均匀度为23.68%,与文献[7]中的24.61%高度符合,且数值范围基本一致,证明了本文采用简化方法计算灯阵热流分布是可行的。

      图  7  本文的红外灯阵热流分布计算结果

      Figure 7.  Calculated heat flux for lamp array on the work-plane

    • 红外灯阵的优化是个多参数优化问题,参数均连续变化且为实数,本文针对此特点引入粒子群算法进行优化计算。相比于遗传算法,粒子群算法具有记忆性,在搜索过程中仅跟随当前最优解,计算不需要编码,也没有交叉和变异操作,只通过内部速度进行更新,因此原理更简单、参数更少,多数情况下收敛速度更快,适用于实数问题求解。

    • 粒子群算法,又称粒子群优化算法(PSO),是通过模拟鸟群觅食行为而发展起来的一种基于群体协作的随机搜索算法[15-16]。假设搜索空间是L维,群体中有N个粒子,则群体中第i个粒子在L维搜索空间的位置为XiXi=(xi1, xi2,…, xiL),i=1, 2,…, NXi经过的最佳位置为Pbesti=(Pi1, Pi2,…, PiL),i=1, 2,…, N。粒子的每个位置代表要求的一个潜在解,代入目标函数即可得到它的适应度值,用于评价粒子的“优”“劣”。整个群体的最佳位置为Gbestg=(Pg1, Pg2,…, PgL),g是最优粒子位置的索引。粒子飞行速度与位置的关系可表示为

      $$ {V_{i}}^{{t + 1}} = \omega {V_{i}}^{{t}} + {c_1}{r_1}\left( {P^{t}_{{\rm{best}}{i} }- {X_{i}^{t}}} \right) + {c_2}{r_2}\left( {G^{t}_{{\rm{best}}g}- {X_i^t}} \right) \text{,} $$ (5)
      $$ X_i^{t + 1} = X_i^t + V_i^{t + 1}。 $$ (6)

      式(5)和式(6)中:ω为惯性权重;$P^{t}_{{\rm{best}}{i} } $ 为第i个粒子第t代为止搜索到的历史最优解;$G^{t}_{{\rm{best}}{g} } $ 为整个粒子群目前为止搜索到的最优解;XitVit分别为第i个粒子当前的位置和飞行速度;c1c2为加速因子,均为非负常数;r1r2是[0, 1]之间的随机数。公式(5)的第1部分是粒子当前速度,表明粒子状态;第2部分是认知部分,表示粒子本身的思考;第3部分是社会认知部分,表示粒子之间的信息共享。

      从式(5)和式(6)可以看出,粒子群算法是一种非常简便的优化算法,单次计算需要调整的参数个数少,各粒子“自我学习”的同时“相互学习”,计算速度非常快。

    • 1)优化参数的确定

      红外灯阵的热流分布其实是每只灯辐射热流在空间上的叠加结果,因此,实质上影响红外灯阵热流分布均匀度的是各灯在xy方向的间距以及h向的高度。传统方法直接使用各红外灯的位置而不是间距作为优化参数,算法每次更新迭代时,一旦灯的位置改变,灯阵间距参数也会发生改变,而且这种改变的影响因素有2个,因此产生了很多计算冗余。

      本文直接使用两灯间距作为优化参数,能够减少这种计算冗余,加速优化。相邻两只灯x向的间距为(∆x1∆x2、…),y向的间距为(∆y1∆y2、…),每只灯距被照面的高度为(h1h2、…)。为简化计算,将灯阵分为4个象限,不同象限的灯阵位置对x轴、y轴对称,以减少待优化的参数总数。

      选取右上角1/4灯阵的高度参数为优化参数。当灯阵的行数m或列数n为偶数时,选取的1/4灯阵的行数为m/2、列数为n/2;当灯阵的行数m或列数n为奇数时,选取的1/4灯阵的行数为(m+1)/2、列数为(n+1)/2。以4行×4列的灯阵和3行×3列的灯阵为例,如图8图9所示。

      图  8  4行×4列灯阵的1/4灯阵选取示意

      Figure 8.  A quarter selection of the 4×4 lamp array

      图  9  3行×3列灯阵的1/4灯阵选取示意

      Figure 9.  A quarter selection of the 3×3 lamp array

      对于灯阵的间距参数,采取同样的方法选取右上角1/4的灯阵间距参数作为优化参数。当灯阵的行数为m、列数为n时,有m‒1个行间距参数、n‒1个列间距参数,因此,当灯阵的行数m或列数n为偶数时,则选取的1/4灯阵间距参数矩阵的行数为m/2、列数为n/2;当灯阵的行数m或列数n为奇数时,则选取的1/4灯阵间距参数矩阵的行数为(m‒1)/2、列数为(n‒1)/2。

      2)优化参数间的约束

      优化计算时,根据工程经验,遵从灯阵间距由中心至外围逐渐减小、灯阵安装高度由中心至外围逐渐降低的原则约束优化参数,以加快收敛速度。以6行×6列的灯阵为例,选取的间距参数矩阵为3×3,选取的高度参数矩阵为3×3,x向的间距由左至右分别为∆x1、∆x2、∆x3y向的间距由下至上分别为∆y1、∆y2、∆y3,高度参数矩阵如图10所示。

      图  10  6行×6列灯阵优化计算时选取的高度参数矩阵

      Figure 10.  Altitude parameter matrix selected for the optimization of 6×6 lamp array

      则有

      $$ \left\{ \begin{gathered} \Delta {x_1} {\text{>}} \Delta {x_2} {\text{>}} \Delta {x_3} \\ \Delta {y_1} {\text{>}} \Delta {y_2} {\text{>}} \Delta {y_3} \\ {h_{11}} {\text{>}} {h_{12}} {\text{>}} {h_{13}} \\ {h_{11}} {\text{<}} {h_{21}} {\text{<}} {h_{31}} \\ {h_{i,j}} {\text{<}} {h_{i,j - 1}}{\text{,}}i,j {\text{>}} 1 \\ {h_{i,j}} {\text{>}} {h_{i - 1,j}}{\text{,}}i,j {\text{>}} 1 \\ \end{gathered} \right. \text{。} $$ (7)

      由于随机迭代过程中各参数之间存在约束关系,能够很大程度上减少计算量,提高迭代速度,且能减弱陷入局部最优的倾向。

    • 为缩短计算收敛时间,做以下处理:

      1)由工程经验,xy向间距一般在0.3~0.5 m之间,高度h一般在0.3~0.6 m之间;

      2)每行灯的y坐标相同、每列灯的x坐标相同;

      3)xy向的优化步长为5 mm,h向的优化步长为10 mm;

      4)对于迭代过程中出现越界的参数,使用图11所示的方式处理。

      图  11  优化参数的越界处理

      Figure 11.  Handling of overrange parameters in the optimization process

      针对典型的矩形航天器散热面,红外灯阵自动优化设计的基本思路为:

      1)使用蒙特卡罗法计算高度h在0.3~0.6 m之间、固定输入电流下的单灯辐射热流分布,计算步长为0.1 m,得到红外灯单灯热流分布数据库。

      2)基于工程经验设定灯阵x向、y向的初始间距和h向的初始高度,根据散热面的长l与宽w计算灯阵的初始行数m与列数n,并计算灯阵的初始热流不均匀度。

      3)划分象限,选取待优化的参数。

      4)采用粒子群算法对m行×n列红外灯阵进行优化,目标函数为灯阵的热流不均匀度δ(∆x1,…, ∆y1,…, h11,…)。

      5)在设定的最大循环次数之内,当δ小于目标值(通常为10%)时停止迭代;当循环次数超出且δ仍未达到目标值时,将灯阵的行数或列数增加1,重复步骤2)~5),直至δ小于目标值。

      6)输出红外灯阵设计参数。

    • 为了显示本文的改进策略对红外灯阵优化速度的促进,以2.1 m×2.0 m的卫星散热面为例,指定灯阵为7行×6列,分别使用遗传算法、粒子群算法和改进的粒子群算法进行灯阵的优化设计,图12为各算法的计算速度对比。

      图  12  各算法的计算速度对比

      Figure 12.  Comparison of calculation speed among different methods

      图12中可以看出:使用遗传算法优化灯阵,热流不均匀度降至10%以下耗时75.8 min,计算305.4 min后的最终热流不均匀度为6.75%;使用粒子群算法进行优化,热流不均匀度降至10%以下耗时56 min,比遗传算法用时缩减26.1%,但计算100 min后的最终热流不均匀度稳定在9.88%,无法继续下降;使用改进的粒子群算法进行优化,计算仅8 min后热流不均匀度即降至9.25%,比遗传算法用时缩减89.4%,且计算55 min后的最终热流不均匀度为7.52%。可见与遗传算法相比,改进的粒子群算法的优化总时长可缩减82%,虽然最终优化结果相对差了0.77%,但加速效果显著,更贴合实际工程应用的要求。

      为更全面评估改进的粒子群算法的优势,将其优化结果与传统粒子群算法的优化结果进行对比,并将二者坐标范围设为一致以便于比较,如图13所示。图中可明显看出算法改进后的优化灯阵热流分布更加均匀,而算法改进前优化结果中灯阵的边缘热流密度往往较大。

      图  13  粒子群算法改进前后的计算结果对比

      Figure 13.  Comparison of calculation results before and after improvement of PSO

    • 对2.1 m×2.0 m的矩形散热面进行红外灯阵自动优化设计,将热流不均匀度的目标值设为8%。自动计算的初始灯阵为6行×6列,迭代计算150次后的热流不均匀度为12.05%;因其未达到目标值,程序自动将灯阵行数增加1行,再对7行×6列的红外灯阵进行优化,热流不均匀度为7.52%时计算完成。同理,对1 m×1 m的散热面进行灯阵优化设计,设计所得的灯阵为4行×4列,热流不均匀度为6.8%。上述灯阵的位置参数和热流分布如图14图17所示。

      图  14  7行×6列灯阵位置参数

      Figure 14.  Optimized result of position parameters of 7×6 lamp array

      图  15  7行×6列灯阵热流分布

      Figure 15.  Optimized result of heat flux distributions for 7×6 lamp array

      图  16  4行×4列灯阵位置参数

      Figure 16.  Optimized result of position parameters of 4×4 lamp array

      图  17  4行×4列灯阵热流分布

      Figure 17.  Optimized result of heat flux distributions for 4×4 lamp array

    • 本文基于粒子群算法建立了红外灯阵自动优化设计的方法,不需要预设灯阵的初始参数,即可自动为不同尺寸的矩形散热面设计出符合热流均匀度要求的红外灯阵,直接供工程使用。为提高优化设计效率,对粒子群算法做了针对性的改进:使用两灯间距作为优化参数,减少计算冗余;遵从灯阵间距由中心至外围逐渐减小、灯阵安装高度由中心至外围逐渐降低的原则约束优化参数,以加快收敛速度,减弱陷入局部最优的倾向。

      改进的粒子群算法与传统遗传算法的计算对比结果证明本文改进策略能有效加快优化速度,且方法简便,可降低采用红外灯阵进行热真空试验的时间成本,对卫星批量化生产具有一定实用价值。

参考文献 (16)

目录

    /

    返回文章
    返回