基于微型空间相机的红外加热笼仿真与设计

黄春玮，张旭升，郭 亮

（1.中国科学院长春光学精密机械与物理研究所，吉林 长春 130033；2.中国科学院大学，北京 100049）

0 引言

对于航天产品，无论是卫星平台还是有效载荷，都需要进行充分而适度的环境模拟试验，以达到发现隐患、优化设计，提高在轨可靠性的目的。在热控分系统中，包含一项极为重要、且贯穿整个研制阶段的环境模拟试验，热平衡试验。通过该试验，达到验证热设计的正确性和合理性、修正热分析数学模型的目的。在热平衡试验中，能否准确地模拟外热流直接决定试验的成败，是影响试验结果的关键热边界条件。对于近地轨道航天器，空间外热流主要由太阳直射、地球反照和地球红外3 部分辐射能量组成[1]。

在热平衡试验中，模拟外热流的方式主要包括太阳模拟器、表面粘贴电阻加热片、红外加热笼、红外灯阵等。其中，太阳模拟器主要模拟与太阳光相接近的光谱分布、辐照度等，但是这种模拟方式成本高、技术复杂，需配套复杂的运动模拟器模拟航天器在轨的姿态变化。红外灯阵模拟的外热流辐照强度可以根据电压来调节，对瞬态的外热流模拟较容易实现，但是红外灯的辐照强度具有明显的方向性，不容易控制模拟表面的热流均匀程度，并且红外灯阵模拟外热流，需要考察模拟表面对红外灯的辐射光谱吸收特性。表面粘贴电阻加热片模拟外热流，模拟的是航天器表面的吸收外热流，要求加热器粘贴后不能对模拟表面造成损伤，并且要在加热器表面喷涂与模拟表面相同的涂层。针对红外加热笼模拟外热流，其以结构灵活、模拟简单等特点而被广泛应用在热试验中。红外加热笼模拟的目标外热流即为空间相机在热分析过程中所得到的轨道周期平均外热流。通常配合热流计，计算该热流密度下热流计敏感片的温度，并以热流计的温度作为反馈来控制红外加热笼的功率[2-5]。

杨晓宁[6-7]等通过优化红外灯阵不同的参数来降低模拟的热流密度的不均匀性。同时，还研究了红外加热笼带条的覆盖系数对到达卫星表面外热流均匀性的影响，仿真结果显示加热笼距离卫星表面越远，覆盖系数越大模拟出的外热流均匀性增长的越快。季琨等[8]针对热试验中红外加热笼之间的干扰进行了优化设计，通过防辐射屏和补偿边设计，可以达到减小热干扰的目的。韩继广等[9]研制出一种运动装置，可以改变红外加热笼与航天器表面之间的相对距离，来实现瞬态外热流的模拟。孙玉玮等[10]对红外加热笼的边缘效应进行了研究，通过增添边缘加热笼以及将边缘加热笼翻转一定角度的方式来改善航天器表面到达热流密度的均匀程度。此外，对圆台型航天器表面和与之相匹配的红外加热笼进行了研究，得出了圆台型加热笼的优化设计参数[11]。谢吉慧等[12]针对红外加热笼与航天器表面结构匹配问题，建立一套数字化结构匹配方法，可以使红外加热笼与航天器表面的匹配成功率达到100%。房红军等[13]对加热笼的布点进行了优化设计，通过对热流计布置位置的优化使外热流的误差均方根由9.2%降低至3.2%。但是调整布点位置改变的仅为实际到达卫星表面的平均外热流的大小，并不能改善到达卫星表面上外热流的均匀程度。

而对于特征尺寸较小的表面，采用红外加热笼模拟外热流时，可能存在总到达能量偏低、热流密度均匀性差等问题。因此，本文将基于某微型空间相机，以对外热流最为敏感的散热面为研究对象，开展红外加热笼仿真分析与优化设计的研究。采用有限元法建立红外加热笼-热流计的系统仿真模型，分析加热笼尺寸和热流计粘贴位置对模拟表面总到达能量和热流密度均匀性的影响。

1 计算模型

在某微型空间相机的热平衡试验中，通过红外加热笼模拟散热面的到达外热流，在一种工况下，模拟出的外热流是恒定的。使用热流计粘贴在散热面的中心来测量模拟的到达外热流。实际试验外热流模拟装置示意图如图1所示，相机整机包络尺寸为980 mm×200 mm×253 mm，散热面宽和高的尺寸为160 mm×167 mm；加热笼由框架、弹簧、钢带、底座组成，其中钢带通过弹簧与框架相连，形成闭合回路。加热笼的宽和高的尺寸为233 mm×255 mm，分别为加热笼左右边框的最大宽度和上边框至底板之间的距离，加热笼钢带的带条宽度为6 mm，厚度0.1 mm，带条间隔为9 mm；黑片热流计的外径为40 mm；红外加热笼与散热面之间的距离为50 mm。

图1 外热流模拟装置示意图Fig.1 Schematic diagram of the external heat flow simulator

为了分析该装置模拟出的外热流情况，采用UG/SST 模块建立红外加热笼-黑片热流计系统的有限元模型，如图2所示。其中，壳单元为15055 个、节点18408 个，共计33463 个。对模型进行不同程度的简化，由于加热笼电阻丝与加热笼框架之间采用弹簧相连接，实际热阻很大，所以在建立有限元模型过程中忽略其耦合关系，朝向被加热表面的一侧喷涂黑漆，红外半球发射率为0.9，背向被加热表面的一侧为材料表面自然状态。结合散热面、红外加热笼的轮廓尺寸和黑片热流计的外径，将散热面划分为9 个区域，分别建立9 个黑片热流计模型，均匀分布在散热面的相应区域，经测量黑片热流计的敏感片及补偿片的红外半球发射率为0.9，并对对应区域的热流计进行编号；热沉温度设置为100 K；散热面喷涂白漆，红外半球发射率和太阳吸收率分别为0.92 和0.12。

图2 有限元网格示意图Fig.2 Finite element mesh diagram

根据热分析结果得出的散热面高温工况的轨道周期平均外热流，以及散热面参数可以计算出热流计的目标温度：

式中：αs为散热面的太阳吸收率；ε为散热面的红外半球发射率；Q1、Q2、Q3分别为太阳直射、地球反照、地球红外的辐射能量；σ为斯特藩-玻尔兹曼常数。

在热试验中，红外加热笼和被加热表面均为平面，并且加热笼的尺寸与被加热表面尺寸基本一致，加热笼与被加热表面之间存在一定距离，这就会出现被加热表面呈现中间受到的热流多，周围受到的热流少的现象[3]。仅采用单一点处的热流计温度值反映出的整个面的热流密度值与实际到达外热流存在误差，因此需要对试验中的误差进行分析。

目前对于均匀程度的定义还不完善，本文采用9个热流计测得的热流密度的均方偏差S来表征到达散热面的均匀程度。

式中：t1～t9为9 个敏感片达到平衡时的温度；t为目标热流计温度。

2 仿真分析

由于红外加热笼主要模拟高温工况下的外热流，这里仅对高温工况的模拟情况进行分析。经热分析得到的轨道周期平均外热流为423.9 W/m2，并以此作为红外加热笼模拟的目标外热流，代入公式(1)计算出高温工况下的目标热流计温度为20.9℃，以此作为控制条件，改变加热笼的加热功率，反复进行稳态计算，当目标热流计达到目标温度且加热笼的占空比为75%时，以此时的功率作为热试验中加热笼需要调节的功率。

9 个热流计敏感片达到平衡时的温度如表1所示，温度云图如图3所示。

表1 传统加热笼控制方法下敏感片平衡温度Table 1 Equilibrium temperature of sensitive sheet under conventional heating cage control method

图3 敏感片平衡时的温度云图Fig.3 Temperature cloud when the sensitive sheet reaches steady state

由表1、图3 可以看到，在热试验中利用红外加热笼模拟外热流，并使用中心位置的RLJ-5 热流计温度作为反馈，来控制红外加热笼的开关，达到平衡时的平均到达热流密度为341.2 W/m2，较目标到达外热流偏低19.6%，不满足总到达能量的保守设计原则。所引起的各点温度偏差最大值可达－34.0℃，9个区域的热流均方偏差为102.0。由图3 的仿真结果可知，该种传统加热笼设计带来的外热流模拟与热分析所得的外热流存在较大误差，由于中心的RLJ-5 敏感片温度与周围的温度相比较高，所以利用中心RLJ-5 热流计的温度来控制加热笼，会导致模拟的平均外热流明显低于在轨期间的外热流，这会使试验工况无法覆盖在轨期间的极端高温工况。并且加热笼模拟出的外热流的不均匀程度很高，因此若要改善这一现象，需对加热笼进行优化设计。

3 优化设计

本文针对上述缺点提出改进方法：

①针对模拟的平均到达热流密度较在轨期间偏低的问题，可以采用周围温度较低的热流计对加热笼进行控制，使得模拟的平均热流密度略大于在轨期间平均热流密度；由于9 个热流计具有对称性，分别采用RLJ-1、RLJ-2 和RLJ-4 热流计的温度作为反馈控制加热笼，进行仿真分析，仿真结果如下，9 个热流计敏感片达到平衡时的温度如表2所示。

表2 不同热流计控制加热笼时敏感片的平衡温度Table 2 Equilibrium temperature of the sensitive sheet when controlling the heating cage with different heat flow meters

经计算，采用RLJ-2 和RLJ-4 热流计的温度来控制加热笼得到的模拟平均热流密度分别为375.7 W/m2和385.5 W/m2，依然小于在轨期间的平均热流密度，而采用RLJ-1 热流计作为控制加热笼的温度数据时，平均热流密度为493.1 W/m2，高于在轨期间的平均热流密度。因此，在试验中可以使用RLJ-1 热流计来测量到达外热流并采用其温度数据控制加热笼。

②针对模拟的平均外热流均匀性差的问题，可以在空间条件允许的条件下，尽可能扩大加热笼尺寸，并使散热面处在加热笼的中心位置，以提高均匀性。另外，扩大加热笼的尺寸势必会导致除被模拟表面的其他表面受到影响，可以采用防辐射屏来遮挡热流。改进后结构示意图如图4所示。

图4 改进后装置示意图Fig.4 Schematic diagram of the improved device

在原有的加热笼尺寸的基础上，分别将宽和高依次增加30 cm 和20 cm，并利用RLJ-1 热流计的温度反馈对加热笼进行控制，进行仿真分析，得到的结果如下，9 个热流计敏感片达到平衡时的温度如表3所示。各个尺寸下，平均热流密度曲线如图5所示，热流密度均方偏差曲线如图6所示。

表3 改进结构的敏感片的平衡温度Table 3 Equilibrium temperature of sensitive sheet with improved structure/℃

图5 各个尺寸下平均热流密度曲线Fig.5 Average heat flow density curves for different sizes

图6 各个尺寸下热流密度均方偏差曲线Fig.6 Mean square deviation curve for different sizes

由表3 和图5、图6 可以看出利用RLJ-1 热流计温度作为反馈，可以提升模拟的平均外热流至高温工况下的轨道周期平均外热流以上；并且扩大加热笼可以有效改善模拟外热流的均匀程度。但是当尺寸扩大到353 mm×335 mm 时，模拟的平均外热流与在轨期间平均外热流的差值为69.2 W/m2，热流密度均方偏差下降至27.0，继续扩大尺寸不能明显减小模拟的平均外热流与在轨期间平均外热流的差值，且无法有效地降低热流密度的均方偏差，也就是说继续扩大尺寸已经无法有效地改善到达散热面的热流均匀程度，由于扩大加热笼的尺寸势必会提升加热笼所需的功率，因此在实际的热试验中可以选用尺寸为353 mm×335 mm 的加热笼模拟外热流，既可以满足模拟的平均外热流的均匀性要求，又可以做到空间和功率方面的平衡。

4 结论

本文根据某微型空间相机在热试验中的外热流模拟装置，采用有限元法建立了系统仿真模型，分析了在传统加热笼控制方法下，到达模拟表面的总能量及热流密度均匀性，并对其结构进行了改进，通过分析改进装置的仿真结果，得出以下结论：

1）对于采用传统加热笼控制方法出现的到达模拟表面的总能量偏低的问题，通过调整热流计粘贴位置，可以使到达模拟表面的热流密度高于在轨期间到达散热面的热流密度，进而满足总到达能量的保守设计原则。

2）对于红外加热笼在该热试验中模拟外热流时出现的热流密度不均匀的问题，通过扩大加热笼能够实现提升均匀度的效果，热流密度的均方偏差从最初的102.0 降低至27.0，并且平均热流密度与目标热流密度的差值从69.2 W/m2下降至13 W/m2，从而使试验工况更接近空间相机在轨期间的工作状态。