Thermal Design and Simulation Verification of High-precision Integrated Star Sensor in Near-Earth Orbit
-
摘要: 某型高精度一体式星敏感器指向精度高,对温度变化非常敏感。其所处近地轨道外热流复杂多变,一体式结构和内热源集中的综合因素不仅导致散热设计困难,而且镜头直接受到内热源发热影响难以保障指向精度。首先,结合轨道参数,安装布局获得星敏感器平均吸收外热流。然后,通过分析外热流与内热源工作情况,采用被动热控和主动热控相结合的热设计方法,并对星敏感器散热面的位置与大小进行设计与计算。最后,根据轨道环境和热控措施并利用热仿真软件进行热分析验证。仿真结果表明,安装法兰温度为19.82℃~20.10℃,镜头轴向温差小于2.23℃,周向温差小于0.48℃,电路盒温度为19.10℃~23.49℃,满足热控指标。通过合理的热控设计保证了极高精度星敏感器的稳定工作条件,星敏感器的热设计合理有效。Abstract: High-precision integrated star sensors have a high pointing accuracy and are highly sensitive to temperature changes. The external heat flow in near-Earth orbits is complex and variable. The comprehensive factors of the integrated structure and the concentration of the internal heat source not only lead to difficulty in the heat dissipation design, but also make it difficult to guarantee the pointing accuracy of the lens directly affected by the internal heat source. First, combined with the orbital parameters, the installation layout provides the average absorbed external heat flux of the star sensor. Subsequently, by analyzing the working conditions of the external and internal heat flows, a thermal design method combining passive and active thermal control is adopted, and the position and size of the heat dissipation surface of the star sensor are designed and calculated. Finally, thermal analysis and verification are performed using thermal simulation software according to the orbital environment and thermal control measures. The simulation results show that the installation flange temperature is 19.82-20.10℃, the axial temperature difference of the lens is less than 2.23℃, the circumferential temperature difference is less than 0.48℃, and the circuit box temperature is 19.10-23.49℃, which meets the thermal control index. The stable working conditions of the extremely high-precision star sensor are ensured by a reasonable thermal control design, and the thermal design of the star sensor is reasonable and effective.
-
Keywords:
- high-precision star sensor /
- thermal design /
- thermal analysis
-
0. 引言
混凝土材料往往因为预埋钢筋的腐蚀、水的侵蚀、连续的冻融循环、荷载等多种因素共同作用,产生裂缝、孔洞和分层等缺陷[1]。这些内部缺陷随着时间的推移会逐渐积聚、扩展,严重威胁着结构的耐久性和安全性[2]。主动红外热像检测方法作为一种无损检测手段,由于具有检测面积大、非接触式和应用操作简单等优点在混凝土结构检测领域得到了广泛的关注,并逐渐被应用到了隧道、桥梁现场检测中。
主动红外热像检测技术是一种利用红外热像仪检测材料内部缺陷的无损检测方法,其检测的基本原理是:材料被热激励后其表面红外辐射量因内部缺陷的存在而表现出差异,红外热像仪通过记录该差异,以获得材料内部的缺陷信息[3]。红外图像的信息提取与分析是红外热像检测技术的核心内容。理论上,根据热图像缺陷区域和正常区域之间的红外辐射差异可获得缺陷的信息。但实际上,红外图像本质是根据缺陷与正常区域之间的温度和辐射率生成的灰度图,而且红外光辐射的能量远比可见光低,因此到达光学镜头的红外辐射能很小,导致了红外图像与可见光图像相比只有较低的对比度,分辨弱小目标和细节能力差[4];另外,红外光子的随机性、红外成像系统的固有特性以及热激励不均匀等给红外图像引入了许多噪声,在一定程度上改变了图像的原始信息,增加了图像分割及后续图像处理的难度[5-6]。传统的基于空域的红外图像处理方法主要分为两个大的步骤:首先,是对红外图像进行降噪和增强,常用的方法有高斯滤波、中值滤波、引导滤波、小波变换、Contourlet变换、三维块匹配算法等[7-10];其次,对红外图像进行边缘特征的提取,常用的方法有:基于边缘的检测算法、基于数学形态的检测算法和最近发展较快的基于网络模型的检测算等[11-12]。传统的基于单张空域的红外图像处理方法虽然在一定程度上可以消除噪声的影响、提高图像的对比度,但是仍存在一些问题。第一,采集的数据集包含大量的热图像,每一帧红外图像都对应着某一时刻,依靠视觉手动选择缺陷信息丰富的红外图像进行缺陷信息的提取的识别费时费力。第二,在图像降噪、增强和图像分割过程中设定阈值会引入主观成分,干扰红外图像的解释。第三,仅仅分析单张红外图像,忽略采集过程中的时序信息可能会导致忽略一些边缘的缺陷信息[13-14]。
针对上述问题,本文提出了一种基于时序信息的红外图像缺陷信息提取方法。首先,通过室内实验制作含缺陷分层的混凝土试块,其次利用主动红外热像检测技术进行红外数据的采集,提取每个像素点的温度特征曲线,然后采用基于时序信息的提取方法进行含分层缺陷混凝土的缺陷提取,并与传统的基于空域的处理方法进行对比分析。
1. 基本原理
1.1 主动红外热像检测技术原理和系统
主动红外热像检测技术中,当以一定的距离对混凝土材料表面持续加热时,物体会吸收入射波的能量并将其转化为热能,并以温度场的形式表现出来。如果试样存在缺陷(本研究采用的是隔热型缺陷,模拟混凝土内部的空气分层缺陷),到达缺陷的大部分能量将反射到试样表面,引起表面正常区域和有缺陷区域温度场的差异。
利用红外热像仪记录整个检测过程中物体表面温度场的空间和时间分布信息,通过对红外图像数据进行分析和处理进而可以提取材料近表层的分层缺陷信息。主动红外热像检测技术采集到的一系列红外图像组成了一个三维的图像数据集[15](尺寸:Nx×Ny×Nt,见图 1(a)),在空间维度上,每一帧对应着某一时刻红外热像仪所采集到的物体表面温度场信息,其像素个数为Nx×Ny;在时间维度上,每个像素点可以看作是一条随时间变化的温度特征曲线(Nt时间点组成),如图 1(b)所示,其中标签NT11 DEFECT表示缺陷区域的温度特征曲线,标签NT11 NORMAL表示正常区域的温度特征曲线。
1.2 基于时序信息的缺陷信息提取原理
主动红外热像技术采集到的三维红外图像数据,其正常区域和含缺陷的区域每个像素点的温度特征曲线在冷却阶段会存在差异,如图 1(b)所示,利用这个特点我们可以采用许多能够分辨这个特征曲线差异值的相关算法进行缺陷信息的提取,本文采用相对简单、容易实施的K-means方法来提取混凝土的缺陷信息。K-means聚类方法的原理是同一类内的实体是相似的,一个类是测试空间中点的集合,同一类内任意一个点到其类中心的距离小于其到其它类中心的距离,我们可以用此算法分辨正常区域和缺陷区域的温度特征曲线[16]。K-means方法的关键是核函数的选取,核函数决定了其分辨能力,本文采用常用的马氏距离进行度量。其具体算法原理如下[17]:
假设给定的数据集为Y={y1, y2, …, ym},其中yi=[X1, X2, …, Xn],将数据对象划分为k类C={c1, c2, …, ck},每个类有一个类中心U={u1, u2, …, uk}。选取欧式距离作为相似性和距离判断准则,计算ck类内各点到聚类中心uk的距离平方和,见公式(1):
$$ J\left( {{c_k}} \right) = \sum\limits_{{x_i} \in {C_k}} {{{\left\| {{x_i} - {u_k}} \right\|}^2}} $$ (1) 聚类的目标是使各类总的距离平方和最小,见公式(2):
$$ J(C) = \sum\limits_{k = 1}^K J \left( {{c_k}} \right) = \sum\limits_{k = 1}^K {\sum\limits_{{x_i} \in {C_k}} {{{\left\| {{x_i} - {u_k}} \right\|}^2}} } $$ (2) 根据最小二乘法和拉格朗日原理,聚类中心uk应该取为类别ck类各数据点的平均值。
K-means算法是一个反复迭代过程,目的是使聚类域中所有的点到聚类中心距离的平方和最小。
2. 实验与数据分析
2.1 实验过程
本次实验混凝土试件的设计强度为C50,尺寸为50 cm×50 cm×20 cm,在混凝土中嵌入4块尺寸为10 cm×10 cm的聚苯乙烯材料(隔热型材料),模拟混凝土中的分层缺陷,混凝土试块见图 2(a)。试验采用美国FLIR公司的A655SC非制冷型红外热像仪,其热灵敏度为30 mK,图像分辨率为640×480,标准测温范围为-40℃~150℃,波长范围为7.5~14 μm,见图 2(b)。主动热激励系统采用自制的配有自动控制加热时间的碳化硅远红外加热板,共2块,每块加热板的加热功率为500 W,见图 2(c)。整个实验过程中加热时间为5 min,冷却时间为10 min;加热完成后,立即用红外热像仪采集试样表面的冷却温度场,整个采集系统见图 3所示,本文选用的是持续加热法。
2.2 红外图像数据
实验采集到的红外图像序列如图 4所示,根据红外图像可以发现一共存在4块分层缺陷,各个缺陷之间存在温度场相互影响的干扰区,而且各个缺陷的轮廓比较模糊,存在被隐藏的缺陷信息。整个采集的红外图像数据的第一帧(0 s)最清晰,其正常区域和缺陷区域存在较大的温差;之后随着冷却的进行,正常区域和缺陷区域的温差逐渐减小,缺陷信息逐渐模糊,至最后一帧很难用肉眼分辨。
3. 红外图像处理和对比分析
3.1 红外图像缺陷提取
主动红外热像检测技术所采集的图像每一帧都含有大量的噪声,为了测试基于时序信息红外图像缺陷信息提取的可行性以及性能,本次在图像特征提取前不对图像进行预处理。首先,将红外图像数据按像素点逐个提取时间序列,组成聚类数据集;然后将时序数据序列输入到K-means程序中,进而得到每一帧各个像素点的聚类信息。为了提高运算效率,本次在图像序列处理中采样间隔为1,其序列长度减小为原来的一半。聚类完成后得到图像中每一个像素点所属的类别,进而完成对图像缺陷信息的提取。
图 5为经过聚类后的提取结果,通过观察发现4个缺陷信息都被提取出来,与原始红外图像相比(图 4),其缺陷信息有了明显的增强,缺陷信息完全肉眼可见,每个识别出的缺陷区域呈近似正方形,完整性相当好,接近于预埋缺陷的形状,而且缺陷提取后的图像序列的时序信息得以保留。
3.2 对比分析
为了测试基于时序信息红外图像分层缺陷的信息提取效果,本文采用基于空域的二维K-means图像缺陷信息提取算法与其进行对比。由于采集的红外图像数据集第一帧(0 s)缺陷信息最清晰,因此选用第一帧红外图像进行基于空域的分层缺陷信息提取。两种方法的分层缺陷信息提取结果见图 6,由图 6(b)所示,基于时序信息的分层缺陷提取方法4个分层缺陷都被完全提取出来,图 6(a)中模糊、隐藏的信息也被提取出来,缺陷细节分辨能力较好。由图 6(c)所示,基于空域信息的分层缺陷信息提取方法4个分层缺陷被识别为一个缺陷,缺陷细节分辨能力较差。通过分析可能是受缺陷之间温度场叠加的相互影响,缺陷之间区域(图 6(a)中虚线框内部所示的无缺陷区域)的温度信息和分层缺陷的温度场信息较为接近,而和图 6(a)的无缺陷边缘温度场信息差异较大,导致基于空域信息的提取方法效果较差;根据含缺陷混凝土试块的物理特性可知,存在分层缺陷和正常区域的温度特征曲线是存在差异的,基于时序信息的提取方法主要是根据每个像素点的温度特征曲线差异来提取特征,其提取缺陷的分辨能力和效果更好。
4. 结论
主动红外热像检测技术中,传统的红外图像处理方法在一定程度上可以消除噪声、提高图像的对比度,但是仍存在一些问题。针对上述问题,本文根据主动红外热成像的数据特点提出了一种基于时序信息的红外图像缺陷信息提取方法。结果表明,基于时序信息的缺陷提取方法是可行的,其可以提取到隐藏的缺陷信息,分层缺陷信息提取效果优于基于空域的K-means方法。
-
表 1 重要结构热控指标
Table 1 Important structure thermal control index
Structure Index content Thermal control index Installation flange Temperature 20℃±0.5℃ Lens cone Axial temperature difference <3℃ Circumferential temperature difference <1℃ Power supply chassis Temperature <25℃ 表 2 极端工况
Table 2 Extreme conditions
Solar position β Heat sink surface coating surface properties High temperature
working conditionWinter
solstice60° End period Low temperature
working conditionSummer
solstice0° Initial period -
[1] 吕建伟, 王领华, 苏生, 等. 微型星敏热控设计及仿真[J]. 红外与激光工程, 2022, 51(11): 265-270. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202211025.htm LV Jianwei, WANG Linghua, SU Sheng, et al. Thermal control design and simulation of micro star sensor[J]. Infrared and Laser Engineering, 2022, 51(11): 265-270 https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202211025.htm
[2] 余成武, 隋杰, 陈超, 等. 高分七号卫星多探头甚高精度星敏感器热设计与验证[J]. 航天器工程, 2020, 29(3): 144-150. https://www.cnki.com.cn/Article/CJFDTOTAL-HTGC202003024.htm YU C W, SUN J, CHEN C, et al. Thermal design and verification of multi-head very high accuracy star sensor onboard GF-7 satellite[J]. Spacecr Eng., 2020, 29(3): 144-150. https://www.cnki.com.cn/Article/CJFDTOTAL-HTGC202003024.htm
[3] 江帆, 王忠素, 陈立恒, 等. 星敏感器组件的热设计[J]. 红外与激光工程, 2014, 43(11): 3740-3745. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201411041.htm JIANG FAN, WANG Zhongsu, CHEN Liheng, et al. Thermal design of star sensor assembly[J]. Infrared and Laser Engineering, 2014, 43(11): 3740-3745. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201411041.htm
[4] 江帆, 吴清文, 王忠素, 等. 星敏感器支架的结构/热稳定性分析及验证[J]. 红外与激光工程, 2015, 44(11): 6. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201511049.htm JIANG Fan, WU Qingwen, WANG Zhongsu., et al. Analysis and verification of structure stability and thermal stability of a bracket of star sensors[J]. Infrared and Laser Engineering, 2015, 44(11): 6. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201511049.htm
[5] 张彧, 赵吉喆, 张翔, 等. 环月轨道一体式星敏感器热设计及仿真验证[J]. 航天器环境工程, 2016, 33(5): 516-520. https://www.cnki.com.cn/Article/CJFDTOTAL-HTHJ201605011.htm ZHANG Yu, ZHAO Jizhe, ZHANG Xiang, et al. Thermal design of an integrated star sensor on lunar orbit with simulation validation[J]. Spacecraft Environment Engineering, 2016, 33(5): 516-520. https://www.cnki.com.cn/Article/CJFDTOTAL-HTHJ201605011.htm
[6] 訾克明, 吴清文, 李泽学, 等. 空间光学遥感器的热设计实例及其仿真分析[J]. 计算机仿真, 2008, 25(12): 4. DOI: CNKI:SUN:JSJZ.0.2008-12-023. ZI K M, WU Q W, LI Z X. Simulation analysis of a space optical remote-sensor's thermal design[J]. Computer Simulation, 2008, 25(12): 4. DOI: CNKI:SUN:JSJZ.0.2008-12-023.
[7] 苗建印. 航天器热控制技术[M]. 北京: 北京理工大学出版社有限责任公司, 2018. MIAO J Y. Spacecraft Thermal Control Technology[M]. Beijing: Beijing Institute of Technology Press, 2018.
[8] 李强, 陈立恒. 复杂外热流条件下红外探测器组件热设计[J]. 红外与激光工程, 2016, 45(9): 73-79. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201609011.htm LI Qiang, CHEN Liheng. Thermal design of infrared detector components in complex heat flux[J]. Infrared and Laser Engineering, 2016, 45(9): 73-79. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201609011.htm
[9] 董奎辰, 郭亮, 黄美娇, 等. 空间相机大功率热源散热设计[J]. 红外技术, 2023, 45(5): 521-526. http://hwjs.nvir.cn/article/id/f7193fd5-4039-46a6-b862-00a13bc28417 DONG Kuichen, GUO Liang, HUANG Meijiao, et al. Heat-dissipation design for space camera high-power heat source[J]. Infrared Technology, 2023, 45(5): 521-526. http://hwjs.nvir.cn/article/id/f7193fd5-4039-46a6-b862-00a13bc28417
-
期刊类型引用(3)
1. 邱欣欣,温强,何婧. 基于参考帧的数字媒体视频图像信息隐藏算法. 吉林大学学报(信息科学版). 2025(02): 377-383 . 百度学术
2. 李环宇,冯国会,刘馨,蒲毅,王涵. 基于红外热成像的围护结构热工性能定量检测方法研究进展. 太阳能学报. 2024(07): 427-437 . 百度学术
3. 张迁,王剑,楚瑞博,陈欢欢. 基于BRISK-BEBLID特征的无人机图像快速配准方法. 激光杂志. 2023(06): 92-98 . 百度学术
其他类型引用(1)