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-4]。NASA喷气发动机风扇机匣、波音787飞机风扇机匣都采用的是编织复合材料。复合材料结构是由纤维铺层叠加而成,在面临碰撞和冲击时,纤维和纤维之间都是由基体传递载荷,冲击容易造成材料基体开裂、分层、脱粘等缺陷,编织层压板冲击损伤模式主要是层间分层,不会出现纤维劈裂等损伤模式[5]。在结构表面看不出来损伤缺陷,目视检测也无法检测出缺陷信息,但其内部可能因瞬时冲击能量的吸收,产生了严重的损伤,使得材料结构失效、材料性能降低,对结构承载能力造成大幅度下降[6-7]。因此,有效及时检测编织复合材料内部损伤或缺陷情况,对于保证复合材料结构件可靠性具有重要的意义。
超声红外热成像技术以超声作为热激励源,激励过程中仅界面类缺陷或闭合裂纹发热,正常区域不发热,属于选择性激励方式,更适合于不规则形状材料界面缺陷或闭合裂纹检测。因此目前超声红外热成像技术主要应用在对金属裂纹的检测[8-11]。韩梦[12]等采用超声红外热成像对铝金属的疲劳微裂纹进行检测,微裂纹最小宽度为20 μm;习小文[13]等采用超声红外热成像技术对航空发动机叶片裂纹进行检测,最小检测裂纹宽度为0.5 μm。G. Busse等[14]最早将超声红外锁相热成像检测技术引入到复合材料界面类缺陷检测。田干等[15]采用超声红外热成像技术对航空复合材料进行数值仿真研究,建立含裂纹缺陷复合材料的有限元模型;金国锋、张炜等[16-17]通过数值计算和试验研究了超声红外热成像技术对复合材料冲击损伤检测的适用性;李胤等[18]采用超声红外热像检测对复合材料层合板24 J和29 J冲击损伤进行检测研究,并与超声C扫进行对比,结果表明超声红外热像检测可检测出复合材料不同损伤形式,并且比超声C扫检测速度快、检测精度高。
目前对冲击能量损伤研究大部分基于1个或2个冲击能量对复合材料损伤的研究,为了解不同冲击能量对编织复合材料损伤及扩展情况,采用10 J、20 J、30 J、40 J、50 J冲击能量对不同编织复合材料进行冲击,并对冲击过程中的冲击能量-时间曲线进行了记录,采用超声红外热成像技术对不同冲击能量损伤试件进行检测与分析。
1. 基本原理
1.1 超声热成像技术的检测原理
超声红外热成像技术是将超声波激励和红外热成像两者结合的新型无损检测方法,属于主动式红外无损检测方法的一种。超声红外热成像技术的检测原理如图 1所示,是通过将高频超声波(20~40 kHz)耦合注入待测试件,超声波经过界面耦合在试件中传播,在超声波的作用下,试件缺陷由于塑性变形、摩擦作用、粘弹性效应等产生热量,使得裂纹、损伤区域温度升高,在热传导作用下,其裂纹或损伤位置对应的试件表面温度也相应的升高,利用红外热像仪采集试件表面的红外热序列图像,对红外序列图像进行处理便可得到其内部缺陷信息。相较于其他激励方式的红外热成像技术,例如:闪光灯激励、红外灯激励、激光激励等,这些激励方式都属于面激励,不论是缺陷区域还是无缺陷区域,试件表面被整体加热,表面吸收的能量通过热传导到试件内部,遇到隔热性缺陷时,热波反射回表面。超声红外热成像技术是对试件局部加热,其仅对缺陷的区域进行加热,无缺陷区域不进行加热,属于选择性加热方式,在超声激励下产生相对位移的缺陷区域温度都将升高,相当于是体加热方式,并且从缺陷区域通过热传导的方式传播到试件表面,仅传播了单次的距离。
1.2 超声热成像技术的检测原理
超声波在试件内界面缺陷处衰减产生热量,缺陷处则成为热源。由此可知,超声波激励产生的热流率q的传播满足热传导方程:
$$ {k_x}\frac{{{\partial ^2}T}}{{\partial {x^2}}} + {k_y}\frac{{{\partial ^2}T}}{{\partial {y^2}}} + {k_z}\frac{{{\partial ^2}T}}{{\partial {z^2}}}{\text{ + }}q{\text{ = }}\rho c\frac{{\partial T}}{{\partial t}} $$ (1) 式(1)中:T为温度,℃;kx、ky、kz分别为在三个维度上的导热率,W/(m·℃);ρ为密度,kg/m3,c为比热容,J/(kg·℃)。
热传导方程(1)的初始条件为:
$$ T(x, y, z, t)=T_{0}\qquad t=0$$ (2) 式(2)中:T0为环境温度,℃。
热传导方程支配热传导及其扩散过程,试件检测过程时间比较短,表面温度变化比较小,忽视试件表面对流和热辐射,则:
$$ \left\{ \begin{array}{l} \frac{{\partial T}}{{\partial t}} - \alpha {\nabla ^2} = q(t) \hfill \\ T(0) = {T_0} \hfill \\ \end{array} \right. $$ (3) 式(3)中:α为热扩散率,m2/s;∇2为拉普拉斯算子,
$$ {{\mathit{\nabla}} ^2}{\text{ = }}\frac{{{\partial ^2}}}{{\partial {x^2}}} + \frac{{{\partial ^2}}}{{\partial {y^2}}} + \frac{{{\partial ^2}}}{{\partial {z^2}}} $$ 分别对上式进行关于空间(x, y, z)的傅里叶变换得到:
$$ \left\{ \begin{array}{l}\frac{{\rm{d}}{T}^{f}}{{\rm{d}}t}-\alpha {ϵ}^{2}{T}^{f}={q}^{f}(ϵ, t)\\ {T}^{f}(ϵ, 0)={T}_{0}{}^{f}\end{array} \right. $$ (4) 式(4)中:qf为对生热率q作温度T的傅里叶变换;$ ϵ $为傅里叶变换过程中的中间变量,没有特殊的物理含义;Tf为对温度T的傅里叶变换。
常微分方程的解析解为:
$$ {T}^{f}(ϵ, t)={T}_{0}{}^{f}{{\rm{e}}}^{-\alpha {ϵ}^{2}t}+{\displaystyle {\int }_{0}^{t}{q}^{f}}(ϵ, t){{\rm{e}}}^{-\alpha {ϵ}^{2}(t-\tau )}{\rm{d}}\tau $$ (5) 对Tf(∈, t)进行傅里叶反变换得到:
$$ \begin{array}{l}T\left(X, t\right)={T}_{0}{{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K(X-ϵ, t){\rm{d}}ϵ+\\ {{\displaystyle \int }}_{0}^{t}{\rm{d}}\tau {{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K\left(X-ϵ, t-\tau \right){q}^{f}(ϵ, \tau ){\rm{d}}ϵ\end{array} $$ (6) 式(6)为超声波在试件内界面缺陷热源公式,式子表示,热源公式是一个复杂函数表达式。
2. 试验与结果分析
2.1 碳纤维复合材料冲击损伤
试验用的碳纤维复合材料是T700编织复合材料层压板,编织技术制造网格状堆叠而成,铺层方式为-45/90/45/0,经过浸渍固化直接形成复合材料结构,试件共5件,试件尺寸均为150 mm×100 mm×4 mm,对5件试件进行冲击试验,冲击能量分别是10 J、20 J、30 J、40 J、50 J,冲击能量-时间曲线如图 2所示,冲击后试件光学图像如图 3所示,图 3(a)为试件正面,为冲击面,图 3(b)为试件背面,光学图像上很难看到试件表面有损伤。
2.2 超声热成像检测系统
超声红外热成像系统如图 4所示,系统主要包括超声激励电源、超声枪、红外热像仪、夹具、图像处理和控制系统。超声枪包括换能器、调幅器、工具杆,最大功率为1200 W,激励频率为20 kHz,超声激励时间可设置为0.2~20 s。通过超声激励电源箱控制器调节输出功率大小,预紧力大小通过机械控制,预紧力大小最大可设置1000 N。红外热像仪型号为FLIR A655sc,属于非制冷型红外热像仪,工作波段为7.5~14 μm,室温下等效噪声温差为0.04 K,图像分辨率为640 pixel×480 pixel。
试验过程中超声枪不直接作用于试件表面,试件固定在夹具上,夹具固定在平台上,超声枪顶在夹具上面,在试验过程中,需要确保超声波有效注入到夹具中,夹具和试件要接触良好,超声枪在20 kHz的激励下,振幅约为几十微米,如果超声枪与夹具之间接触不好,则超声能量无法传递到夹具上,因此设置超声枪与夹具的预紧力为200 N,预紧力是通过旋转预紧力手柄迫使预紧力系统单元向前移动,内部弹簧受到挤压,从而对超声枪头产生压力,超声枪头压紧夹具。校正好红外热像仪,超声激励时间为1.5 s,采集频率为50 Hz,采集时间为10 s。
2.3 检测结果分析
试验中共有5种不同冲击能量的试件,选取30 J能量冲击能量的试件3#进行分析,图 5所示为超声红外热成像系统检测试件3#所示的红外序列热图,从图 5中可以看出,超声激励前,看不到任何热斑信息,随着超声持续激励,图像中热斑信息逐渐清晰,0.08 s时稍微可以看到零星散点,到0.5 s时逐渐形成了较为清晰的热斑,超声激励时间为1.5 s,红外序列中时间为1.6 s时热斑最为清晰,在超声激励结束后,冲击损伤区域不能持续发热,仅能靠自身热传导把热斑区域的能量传递到低温区域,逐渐达到温度平衡,从红外序列图中表现为图像越来越模糊,对比度越来越低,到8.66 s时已不能分辨出热斑轮廓信息,随着热传导的进行,将不能检测到任何热斑信息,所检测到的红外图像与0.02 s时图像一致。
图 6为不同位置处A点、B点、C点、D点、E点的温升变化曲线,由图 6曲线可以看出,超声激励后,损伤区域温度迅速上升,且各点温度上升不尽相同,A点为试件3#的温升最高点,最高温升可达到8℃,B点和C点为其他区域温升点,这3个点温升曲线有一个特点,温升曲线都是呈现先上升再下降的过程,这3个点代表着试件内部有持续发热的热源,这些热源就是在超声激励下,损伤处由于摩擦生热等原因产生的热源。D点温升曲线只有温升的过程,没有下降的过程,从图像上看,C点位热斑区中间位置,D点处于C点上方,D点处于热斑区域的边缘,边缘位置损伤比较轻微,温升没有那么高,但依然会在持续激励过程中温度不断上升,等热激励结束以后,随着热扩散的进行温度将下降,由于其处于热斑的边缘区域,热斑处的热量会由于热扩散向低温区域传播,C点和D点的温度由于热扩散会逐渐接近,因此导致D点温度下降曲线基本是平的。E点代表着无热斑区域,即无损伤区域,温升曲线基本是平坦的,E点对应位置的内部没有发热源,与上文分析一致。
为了更形象描述不同帧之间空间温度变化曲线,分别取图 5所示前6个时间点的对应图像的空间温度变化曲线,直线所对应的位置为图 7所示红色空间直线对应的位置,左边是起始点,对应的距离为0,右边对应的距离为141.5 mm。曲线横坐标为距离,纵坐标为温度,环境温度约为22℃,t=0.02 s对应的空间曲线是超声激励之前的空间温度曲线,温度基本相同,没有太明显变化,t=0.08 s红外图像上可以看到零星散点,红色空间直线正好位于散点上,在空间温度曲线上正好可以看到损伤处对应位置处有凸起,随着超声激励的加载,损伤区域空间温度曲线温度逐渐升高,并且其对应的宽度稍微有所增加,但温度升高的点仅限于损伤区域,非损伤区域对应的曲线基本是不变的,重叠在一起;损伤区域t=0.24 s对应的空间温度曲线就比t=0.08 s对应的空间温度曲线要高,同样的损伤区域t=0.50 s对应的空间温度曲线比t=0.24 s对应的空间温度曲线高,到超声激励结束时间附近t=1.6 s,空间温度曲线升高最高值,图 7对应的最高曲线。超声激励结束后,损伤区域不再主动发热,内部热源消失,温度开始下降,空间温度曲线也随之下降,t=2.28 s损伤区域对应的曲线的值就比t=1.6 s对应曲线值有所下降。图像中红色箭头指示的空间距离约60 mm时间为t=0.50 s和t=1.60 s损伤区域的边缘点,温度比其旁边两个尖峰都要低,当超声激励结束以后,随着热扩散的进行,t=2.28 s对应的温度反而要比t=1.60 s对应的温度要高,此时该区域对应的2个尖峰也变成了一个大而平滑的尖峰曲线,在红外序列图像中就表现为图像变模糊,与图 5所示图像相对应。
2.4 不同冲击能量试件结果与分析
图 8为不同冲击能量冲击试件的检测图像,从光学图像上都看不到表面的损伤信息;从检测图上看,冲击能量为10 J的1#试件检测图上有3个小的热斑区域,冲击能量为20 J的2#试件看出热斑区域连成了一个整体,冲击能量为30 J的3#试件可以看到热斑已延展开,冲击能量为40 J的4#试件检测图上看出形成的热斑延展成一条长的缝,冲击能量为50 J的5#试件检测热图上热斑几乎贯穿对角线,并且还有额外2条分支。试件是二维编织技术制作的二维编织复合材料,编织使得纤维之间连成了有机的整体,使得二维编织复合材料展现出了较强的整体性。在冲击能量的作用下,外表面几乎看不出损伤信息,当编织复合材料受到小的外部冲击后,不足以破坏有机的整体,损伤区域仅限于冲击点周围;当编织复合材料受到大的外部冲击后,其内部将会产生大的面积分层,损伤展现出延展性。
2.5 损伤区域识别
衡量超声红外热成像技术对编织复合材料检测能力关键是对损伤区域进行提取,采用双阈值分割,对图 5中时间为t=0.50 s、t=1.60 s、t=2.28 s图像进行分割,分割效果如图 9所示,无法达到理想的分割效果。底部区域是超声激励中夹具位置,夹具在超声激励作用下温度升高,通过热传导会把试件底部的温度升高,采用双阈值分割方法无法把该区域给剔除,对比3张图分割效果可以看出,选用不同时间的图像,分割的结果不尽相同。
由于红外序列图中单帧图像由于热传导一直在发生变化,采用单帧图像分割算法会由于选取的帧不同而导致分割区域的差异,为了准确提取损伤区域,本文采用曲线分类算法,利用每一个像素点的温升曲线进行分类,定义损伤区域强弱程度,温升曲线有上升沿和下降沿的点为热斑中心区域,该区域损伤比较严重,图 6所示对应的A、B、C三个点,该像素点置红色;温升曲线有上升沿但无下降沿的点为损伤比较轻微,图 6所示对应的D点,该像素点置绿色;温升曲线基本不变化的为无损伤区域,图 6所示的E点,该像素点置黑色。不同冲击能量试件分割效果如图 10所示。对比于双阈值分割,采用曲线分类算法分割效果比较好,且图像中有两种颜色表示,红色区域基本对应于检测图像中的亮热斑区域。
表 1为不同冲击能量试件损伤面积统计,从表格中可以看出,冲击能量越大,损伤面积越大,但是试件4#损伤面积反而比试件3#损伤面积小。注意到在冲击试验中,图 2所示冲击能量-时间曲线,试件4#峰值冲击能量是40 J,正常冲击结束处的能量是30 J左右,但冲击结束处的能量只有20 J,反而低于试件3#冲击结束处的能量24 J,所以造成这种对应关系不匹配的情况。
表 1 不同冲击能量试件损伤面积Table 1. Damage area of specimens with different impact energySpecimen Red area/mm2 Total area/mm2 1#-10J 630 926 2#-20J 1475 2525 3#-30J 2909 3929 4#-40J 2896 3383 5#-50J 4517 6291 由于试件4#冲击结束处能量偏低,因此把其他4个试件的冲击能量和损伤面积作为损伤面积关系图,如图 11所示,可以发现损伤面积和冲击能量成线性关系,与碳纤维层合板冲击损伤面积和冲击能量关系曲线相同[19]。
3. 结论
编织复合材料冲击损伤主要涉及复合材料的分层损伤,这些损伤无法采用目视检测,且损伤会使得结构承载能力严重降低,从而威胁复合材料构件的安全使用。制作了编织复合材料层压板试件共5件,并对不同冲击能量下冲击损伤开展了超声红外热成像技术研究。超声激励试件时,冲击损伤区域由于摩擦等因素产生热量,相当于内部有热源的存在;无损伤区域不生热,使用红外热像仪采集红外图像序列便可实现对损伤区域检测,本文以试件3#为例对超声激励过程的温升曲线、空间温度曲线进行了详细的分析。并对10~50 J冲击能量试件进行了对比分析,最后采用曲线分类算法对损伤区域进行了定量识别,发现编织复合材料损伤面积和冲击能量成线性关系。后续的工作是采用其他技术手段(如超声C扫描无损检测技术)确认试件损伤区域面积,以及对损伤区域的评估,将对损伤定量化及曲线分类算法得到更进一步提高。
-
表 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