空间目标轨道外热流计算及辐射特性研究.pdf

文章编号2097-1842（2024）01-0187-11空间目标轨道外热流计算及辐射特性研究郑鸿儒1,2，马岩1*，张帅1，王建超2，曲友阳2（1.北京跟踪与通信技术研究所,北京100094；2.长光卫星技术股份有限公司,吉林长春130000）摘要：本文针对空间目标受到的太阳辐射、地球辐射、地球反照辐射，采用蒙特卡洛（MonteCarlo）法，基于非结构四面体网格编写了仿真程序，并对计算结果进行了对比验证。进一步地，对太阳同步轨道卫星受到的轨道外热流，采用带帆板的网格对有无遮挡情况下各表面受到的轨道外热流进行了分析。结果显示，在对地模式下考虑遮挡后，Y 表面平均热流值降低了 53.79W/m2，+YZ 侧帆板表面平均热流值降低了 32.05W/m2。结合表面材料属性，分析了各表面的温度特性，并结合帆板温度的在轨遥测数据，验证了计算的准确性。最后，计算了两种模式下各方向的红外辐射强度。结果表明，不同观测模式下各表面受热流的影响不同，对地模式下各表面温度随时间变化较大，而对日模式下各表面热流较为稳定。两种模式下，太阳能帆板的温度较高，辐射强度较大，具有明显的红外特征，便于开展红外观测。关键词：轨道外热流；空间目标探测；红外辐射；蒙特卡洛方法中图分类号：V416.5文献标志码：Adoi：10.37188/CO.2023-0033CalculationoforbitexternalheatflowandradiationcharacteristicsofspacetargetZHENGHong-ru1,2，MAYan1*，ZHANGShuai1，WANGJian-chao2，QUYou-yang2（1.Beijing Institute of Tracking and Telecommunications Technology,Beijing 100094,China；2.Chang Guang Satellite Technology Co.Ltd.,Changchun 130000,China）*Corresponding author，E-mail:Abstract:Inthispaper,thesolarradiation,theearthradiationandtheearthalbedoradiationreceivedbythespacetargetaresimulatedbyMonteCarlosimulationmethod,andthesimulationprogramiswrittenbasedontheunstructuredtetrahedralgrid,andthecalculationresultsarecomparedandverified.Furthermore,fortheorbitexternalheatflowreceivedbythesun-synchronousorbitsatellite,thegridwithsolarpanelsisusedtoanalyzetheorbitexternalheatflowreceivedbyeachsurfacewithorwithoutocclusion.TheresultsshowthattheaverageheatflowvalueofYsurfacedecreasesby53.79W/m2afterconsideringocclusionintheearth-pointingmode.Theaveragesurfaceheatflowvalueof+YZsidepaneldecreasedby32.05W/m2.Thetem-peraturecharacteristicsofeachsurfacearegivencombinedwiththepropertiesofsurfacematerials,andtheaccuracyofthecalculationisverifiedbycombiningwiththeon-orbittelemetrydataofthesolarpaneltem-perature.Finally,theinfraredradiationintensityineachdirectionofthetwomodesiscalculated.Theresults收稿日期：2023-02-18；修订日期：2023-03-13基金项目：国家自然科学基金（No.61890965）SupportedbyNationalNaturalScienceFoundationofChina(No.61890965)第17卷第1期中国光学（中英文）Vol.17No.12024 年 1 月ChineseOpticsJan.2024showthattheinfluenceofheatflowonthesurfaceisdifferentunderdifferentobservationmodes.Thetem-peratureofeachsurfacevariesgreatlyovertimeintheearth-pointingmode,whiletheheatflowoneachsur-faceisrelativelystableinthesun-pointingmode.Underbothmodes,thetemperatureofthesolarpanelishigher,theradiationintensityislarger,andithasobviousinfraredcharacteristics,whichfacilitatesinfraredobservation.Keywords:orbitexternalheatflow；spacetargetdetection；infraredradiation；MonteCarlomethod1引言目前各国对空间资源的争夺愈演愈烈，空间态势感知技术对国防的战略意义日益凸显，如何对空间目标进行有效监视成为亟需解决的问题之一。空间目标由于距离太远，往往在探测器上成像只有几个像素，可提取的目标信息十分有限1。为了维持载荷正常工作，空间目标不可避免地会向外辐射多余热量，因此，通过目标红外信息获得目标的工作状态是空间态势感知系统的重要组成部分2。然而，空间目标红外探测试验成本较高，通过建立理论模型，开展仿真分析已成为研究目标红外特性的重要手段。其中，目标轨道外热流计算是重要一环，主要包括太阳辐射、地球红外辐射和地球反照辐射。太阳辐射及遮挡的影响相对容易计算，而地球红外辐射和反照辐射则需要关联目标姿态和受照射情况，计算过程十分复杂。国内外的相关学者开展了大量研究，主要分为针对简单结构的解析法或积分法，以及处理复杂结构的蒙特卡洛（MonteCarlo）法。前者在处理常见六面体结构时比较方便，如易桦3等提出了一种针对圆轨道航天器外热流的计算方法，对偏航姿态瞬态外热流进行了分析。GARZN4等人针对3U 卫星建立了简化的热流和温度模型，分析了贝塔角和导热率的影响。李志松5等人针对微纳卫星开展了轨道外热流计算及在轨温度分析。李世俊6、吴愉华7等人将相机简化为六面体结构，对太阳同步轨道和地球同步轨道外热流分别进行了分析，得到变姿态条件下的热流值。积分方法虽然能够快速获得初步结果，但对于处理复杂结构及遮挡问题时则比较困难，且精度不够。蒙特卡洛方法是商业仿真软件常用的计算方法，在处理复杂结构航天器及多次反射等情况具有明显优势。Atra 等人8对卫星热环境的建模和分析进行了综述，并对比分析了几款商业软件的计算效果。韩玉阁等人9利用蒙特卡洛方法计算卫星红外辐射特征，结果表明卫星的散热面是判断卫星是否失效的依据。潘晴10等使用反向蒙特卡洛方法对带天线结构的立方体卫星进行了轨道外热流的计算，并与商业软件进行了对比分析，以验证其精度。刘巨11使用 STK 软件获得太阳矢量关系，使用 IDEAS/TMG 模块获得了轨道周期内的热流变化。到目前为止，国内采用自主编程对空间目标外热流的精细化计算研究的较少。本文采用蒙特卡洛方法，基于非结构化四面体网格，编写空间目标轨道外热流仿真软件，采用 OpenMP 并行加速光线计算，获得目标在任意时间、任意轨道、任意姿态下的轨道外热流，并进行了对比验证。进一步地，对存在遮挡情况下的目标表面热流情况进行了仿真分析，并结合表面材料属性，对表面温度和辐射特性开展了研究。2计算方法本文卫星目标轨道外热流的计算软件主要采用矢量坐标变换法，计算顺序如下：首先读取卫星轨道参数，根据卫星轨道信息，计算得到太阳和卫星在 J2000 坐标系下的位置坐标，获得转换矩阵。然后，读取网格文件，遍历每一个表面网格，计算每个表面受到的太阳辐射、地球辐射和地球反照辐射。判断卫星所在位置是否处于地影区，如果在地影区则将太阳辐射和地球反照太阳辐射值设置为 0，不在地影区则总辐射为 3 种辐射产生的热流之和。最后判断表面是否被遮挡，如果被遮挡，则热流值设置为 0。软件运行时卫星轨道信息输入为 6 个（半长轴、偏心率、倾角、升交点赤经、近地点幅角、真188中国光学（中英文）第17卷近点角），并在内部完成递推。在不进行姿态调整时，本体系与轨道系重合。轨道系定义为，+X 指向飞行方向，+Z 指向地心，+Y 遵循右手定则。卫星轨道外热流计算中，暂不考虑目标自身温度产生的辐射。下面将逐一介绍太阳直接辐射、地球红外辐射和地球反照辐射的计算方法。2.1太阳辐射计算Ssun太阳辐射及其遮挡情况计算与地球辐射和反照辐射计算不同，不需要对大量光线进行随机计算。太阳在 J2000 系的位置计算过程参见文献7。在获得太阳矢量后，目标表面接收到的太阳辐射的热流密度可表示为qs=Sc|cos()|(1),（1）ScSc式中，为太阳辐射照度，为表面网格法向与太阳光向量的夹角，为面元反射率。太阳辐射照度随季节变化，通常取平均太阳辐射照度=1367W/m2进行计算：Sc=Sc(Rs/Rs),（2）RsRsSceee式中为日地距离，为平均日地距离。目标在光照区受到太阳辐射，在地影区则设置为 0。地影区通过计算当前时刻卫星视角下的地球张角判断，当卫星与地球连线和卫星与太阳连线的夹角大于时，卫星受到光照。的计算方法如下：e=arcsin(ReRsat),（3）ReRsat式中为地球半径，为卫星与地心的距离。2.2地球红外辐射计算地球红外辐射及反照外热流示意图如图 1 所示。计算地球红外辐射时认为地球处于热平衡状态，任意表面的辐射强度均匀且相等，则地球辐射热流的表达式3,12为：qe=1e4Sc(1)xcos1max(cos2,0)L2ds.（4）e12dsdALe=scos1max(cos2,0)/L2ds其中，为地球反射率，本文设定为 0.3，为面元红外发射率；和分别为地球表面微元与目标表面微元的连线与二者法向之间的夹角，为目标表面微元与地球表面之间的距离。令，为地球红外辐射角系数，是主要求解对象。Shadow areaEarthSunTarget surfacedsdAnns21图1地球红外辐射及反照外热流示意图Fig.1Schematicdiagramofearthsinfraredradiationandalbedoexternalheatflowek=Re/Rsate通过积分方法计算时，主要判断面元法向与目标-地球连线的夹角 与的余角之间的关系，令，则可由下式得到：0arccoske=k2cos当时，；arccosk (arccosk)当时，e=k2cos+1/(/2(1k2)1/2(k2cos2)1/2arcsin(1k2)1/2/sin)k2cos arccos(1k2)1/2/ktan)；(arccosk)e=0当时，。地球红外角系数公式的推导详见文献 13。在使用蒙特卡洛方法时，则通过统计光线实现地球红外角系数的计算，此时地球辐射热流表达式变换为：qe=2NNabsi=1(1e)Sccosi(1)4,（5）Nabsi式中，N 是面元发出的总光线数量，是被吸收的光线数量，是光线发射位置的天顶角。2.3地球反照辐射计算地球反照辐射计算中将地球设定为漫反射，依据兰贝特余弦定理，则目标表面接收到的地球反照辐射热流3,12可表示为：qa=eSc(1)xcos1max(cos2,0)max(cos,0)L2ds,（6）erer=xcos1max(cos2,0)max(cos,0)/L2ds令为地球反照角系数，。由于涉及光线与地球表面交点，以及和太阳的位置关系，计算较为复杂。第1期郑鸿儒,等:空间目标轨道外热流计算及辐射特性研究189在工程上，地球反照辐射热流可以由地球辐射热流得到:qa=qee(1e)/4cos,（7）其中，为卫星-地球连线与太阳光的夹角。在使用蒙特卡洛方法计算时，地球反照辐射热流变换为：qa=2NNabsi=1eSccosi(1)max(cos,0).（8）2.4温度场计算详细计算空间目标温度场的分布十分复杂，本文对计算条件进行了一定简化。目标在大气层外时，忽略热对流；除帆板外，忽略星体各表面之间的导热和辐射；除太阳和地球外，忽略其他天体辐射的影响。在这种条件下，空间目标接收到的外热流即为太阳辐射、地球辐射和地球反照辐射。星体表面热平衡方程如下：Nsi=1(sqs,i+aqa,i+eqe,i)Ai/A+qin=T4.（9）对于太阳能帆板电池片表面，平衡方程如下：Nsi=1(sqs,i+aqa,i)(1s)+eqe,i)Ai/A(T1T2)/h=T14.（10）对于太阳能帆板背板表面，平衡方程如下：Nsi=1(sqs,i+aqa,i+eqe,i)Ai/A+(T1T2)/h=T24,（11）sea qs,i、qa,i、qe,iAiiNsATs其中，、分别为目标表面对太阳辐射、地球辐射、地球反照的吸收率，太阳辐射和地球反照辐射吸收率取值为材料吸收率，地球辐射吸收系数为表面发射率；分别为第 i 个面元接收到的太阳辐射、地球反照、地球辐射的热流；是第 个面元的面积，是该表面上面元的总数目，是该表面所有面元的总面积，为斯忒芬玻耳兹曼常数，=5.67108W/(m2K4)；为目标表面温度；T1为帆板电池的表面温度；T2为帆板背板的表面温度；为太阳能帆板的光qinh电转换率，设为 0.2；为目标表面的内热源，此处等效为面热源，对于本体表面，其值为 20W/m2。太阳能帆板蜂窝材料厚度 设为 20mm，热导率设为 1.7W/(mK)14，表面材料属性如表 1 所示。表1表面材料热参数15Tab.1Thermalparametersofsurfacematerial表面材料吸收率发射率本体氟460.350.68太阳帆板电池片0.820.81背板0.880.862.5计算结果验证为了验证仿真软件计算结果的准确性，针对800km 太阳同步轨道和地球同步轨道开展轨道外热流计算，并与文献 5 中给出的结果进行对比分析。仿真中采用的光线数为 10000，轨道递推间隔为 90s，仿真日期选择为春分日，对比结果如图 2图 3（彩图见期刊电子版）所示。其中大写字母（线条）代表文献中给出的各表面结果，小写字母（符号）代表本文仿真结果。其中 q 为包括太阳辐射、地球辐射和地球反照辐射的轨道总热流。由图 2图 3 可以看出，两种轨道高度下轨道外热流的计算值与文献中的计算值误差小于 5%，具有较高的一致性。005001 0001 50012True anomaly/radq/(Wm2)3456+XX+YY+ZZ+xx+yy+zz图2800km 太阳同步轨道外热流对比Fig.2Comparison of external heat flow of the 800kmsun-synchronousorbit此外，对地球辐射角系数也进行了验证，如表 2 所示，给出了与文献 16 的对比结果。可以看出，两者一致性较好，可以认为本软件计算精度较高。190中国光学（中英文）第17卷005001 0001 50012True anomaly/radq/(Wm2)34+XX+YY+ZZ+xx+yy+zz56图3地球同步轨道外热流对比Fig.3Comparisonof external heat flow of the geosyn-chronousorbit表2典型位置地球红外角系数随平板俯仰角变化对比Tab.2Comparisonof earth infrared angular coeffi-cientsvaryingwithplatepitchangleattypicalpositions/()本文结果Ref16结果00.91290.9118300.79610.7961600.56730.5639900.31390.31251200.11070.11001500.00810.00771800.00000.00002.6计算域进一步地，建立带帆板的模拟卫星模型网格，计算域如图 4 所示。卫星本体系与轨道系重合（本体系原点为卫星质心）。帆板为固定式帆板，分布在+Y 和-Y 侧，卫星本体尺寸为 30cm20cm30cm，帆板尺寸为 30cm40cm2cm，帆板 边 缘 与 星 体 Y 侧 平 面 4cm，网 格 数 量 为28 万。XYZ图4计算域Fig.4Computationaldomain2.7遮挡计算在计算目标本体结构对太阳辐射的遮挡时，由网格中心点沿着卫星指向太阳的矢量方向发射一条光线。通过给定足够大的步长，保证矢量终点在计算域外。然后遍历目标所有面元，如果光线矢量穿插了其他本体表面，则将该光线所在面元太阳辐射热流值设置为 0。2N在计算对地球红外辐射和地球反照辐射的遮挡时，由表面面元向空间内随机发射条光线，跟踪每条光线的运动。在计算光线运动路径时，采用临近网格检索方法，通过计算光线路径与四面体网格的 4 个三角面元相交情况，当光线穿插本体其他表面时删除光线，以此循环，直到光线到达计算域边界。此时，根据光线、目标-地球连线、地球-太阳连线等的关系，统计地球红外辐射和地球反照辐射量。3分析与讨论3.1不同模式下轨道外热流计算空间目标的辐射特性除与所处位置、结构形状、表面属性因素相关外，还受姿态变化的影响。对于天基目标观测尤为明显。通常，三轴稳定卫星在轨长期模式可分为三轴对地、三轴对日等。本节对固定帆板的卫星在两种模式下的春分日轨道外热流进行了计算。对日坐标系定义为，Z 轴指向太阳，+Y 轴指向黄北极，X 轴遵循右手定则。轨道参数为 535km 太阳同步轨道，降交点地方时为 10:30，计算起始时刻为春分日UTC 时间12:30:00，迭代时间间隔为120s，轨道周期为5720s，其中 2838s 至 4901s 为地影区。图 5图 10（彩图见期刊电子版）给出了卫星本体各表面在一个轨道周期内的 3 种热流变化情况，其中 mode1 代表三轴对地模式，图中以实线表示，mode2 代表三轴对日模式，图中以点划线表示。可以看出，在对地模式下，除+Y 面以外，本体其余各表面均受到太阳辐射，由于太阳辐射在总辐射中占比较大，因此，在卫星设计上一般选择+Y 面为散热面。在各面中，Y 面受到太阳帆板遮挡，太阳辐射在 0550W/m2范围内波动，相较于其他受晒表面幅值较小。+Z 面长期对地，只第1期郑鸿儒,等:空间目标轨道外热流计算及辐射特性研究191有进出地影区时短时间受晒，峰值在 510W/m2左右。除Z 面以外，其他表面均受到地球辐射及地球反照辐射影响。对地模式下，各表面受到的地球辐射恒定，反照辐射随时间变化。002004006008001 0001 2001 4002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图5+X 面外热流曲线Fig.5Externalheatflowof+Xsurface002004006008001 0001 2001 4002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图6X 面外热流曲线Fig.6ExternalheatflowofXsurface00501001502002503003504002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图7+Y 面外热流曲线Fig.7Externalheatflowof+Ysurface在对日模式下，各表面热流情况与对地模式有显著不同。太阳辐射仅存在于Z 面，在光照区恒定，约为 1378W/m2，其他表面只受地球辐射和002004006008001 0002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图8Y 面外热流曲线Fig.8ExternalheatflowofYsurface002004006008001 0002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图9+Z 面外热流曲线Fig.9Externalheatflowof+Zsurface004008001 2001 6002 0002 000Time/s4 0006 000q/(Wm2)qs(mode 1)qe(mode 1)qa(mode 1)qs(mode 2)qe(mode 2)qa(mode 2)图10Z 面外热流曲线Fig.10ExternalheatflowofZsurface192中国光学（中英文）第17卷反照辐射，与太阳辐射相比量级较小，因此受热相对稳定，在卫星热设计上，X、Y、+Z 面均可设计为散热面。对于地球辐射和反照辐射来说，对日模式和对地模式相比，主要不同是由于姿态变化产生的在时间和幅值上的差别。对于载荷在+Z 面的卫星来说，从图 9 可以看出，卫星在两种模式下受到的 3 种辐射值均变化不大，热环境比较温和。3.2遮挡影响分析图 11（彩图见期刊电子版）给出了计算起始时刻对地模式下星体各表面受到的总外热流云图。从图中可以看出，由于卫星降交点地方时选择为 10:30am，太阳从Y 侧照射星体，星体和帆板的Z、+X，Y 侧受光照，+X 面热流值在1000W/m2以上，其他面在 600W/m2左右。受帆板遮挡影响，Y 侧部分位置不被太阳直接照射，+Y 侧帆板靠近星体边缘被遮挡，说明本文编写的软件可以很好地处理复杂结构体遮挡情况。YXZq/(Wm2)1 000800600400200图11外热流分布及遮挡影响示意图Fig.11Schematicdiagramofexternalheatflowdistribu-tionandoccludingeffects表 3 给出了一个轨道周期内各表面的平均热流情况。对于三轴对地模式（mode1），除+Y 面平均热流在 90W/m2左右外，各表面热流分布相对均匀，在 340400W/m2左右变化。而对于对日模式，Z 侧的星体表面和帆板表面热流值在 972W/m2左右，其他表面则小于 150W/m2。针对结构遮挡的影响开展了仿真分析，无遮挡情况下不考虑光线被星体结构或帆板遮挡。计算结果表明，遮挡的影响主要表现在对地模式下，由于帆板的遮挡，Y 表面遮挡后热流值降低了 53.79W/m2，由于星本体的遮挡，+YZ 侧帆板表面遮挡后热流值降低了 32.05W/m2。对本体其他表面影响微弱，对于对日模式则无影响。表3各平面一个轨道周期内的平均外热流Tab.3AverageexternalheatfluxesofeachsurfaceinoneorbitalperiodW/m2SurfaceMode1Mode2+X390.5187.92X373.86116.75+Y89.96100.83Y408.43120.34+Z346.37148.28Z399.60972.33+Y+Zsolarpanel342.52148.17+YZsolarpanel367.55972.22Y+Zsolarpanel346.43148.29YZsolarpanel399.60972.233.3表面温度结合各表面热流变化及表 1 中给出的表面材料热参数，根据式（9）式（11）计算得到各表面温度。图 12图 13（彩图见期刊电子版）分别给出了对地和对日模式下的各表面温度。由于各表面的质量和比热容未知，本节计算结果仅考虑稳态解，相较于实际情况波动较大。可以看出，对地模式下各表面温度在不同时刻变化较大，其中，星本体Z 面温度变化最大，约为 180K，+Y 侧温度变化最小，约为 22K。帆板由于存在导热，电池片侧（Z）表面温度在光照区和地影区的差值相较于本体Z 面较小，约为 120K；随着时间变化，+X 和X 侧表面受到阳光照射，温度变化较大，幅值约在 140K 波动。Y 面在光照区和地影区温度都比较恒定，分别为 283K 和 200K，+Z 表面温度由于在进出地影区时受到太阳照射，存在两个峰值，其余时间波动较小。0100150200250300350400450500+X+Y+Z solar panel+YZ solar panelY+Z solar panelYZ solar panelX+YY+ZZ2 0001 0003 0005 000Time/s4 0006 000Temperature/K图12对地模式下的各表面温度Fig.12Temperatureofeachsurfaceintheearth-pointingmode第1期郑鸿儒,等:空间目标轨道外热流计算及辐射特性研究193在对日模式下，除帆板和本体Z 面在光照区和地影区变化较大外，其余表面变化幅度较小，其中波动最大的是本体+Z 面，波动区间约为 132K。对于帆板和本体Z 面，光照区温度在 340K 左右波动，地影区在 220K 左右波动，整个区间相对恒定。01001502002503003504004505002 0001 0003 0005 000Time/s4 0006 000Temperature/K+X+Y+Z solar panel+YZ solar panelY+Z solar panelYZ solar panelX+YY+ZZ图13对日模式下的各表面温度Fig.13Temperature of each surface in the sun-pointingmode3.4表面温度对比验证为了验证轨道外热流计算的准确性，对“吉林一号”某卫星帆板温度进行仿真与在轨实测值验证分析。卫星运行在 535km 太阳同步轨道，降交点地方时为 12:00，帆板位于Y 侧，帆板运行模式为对日模式。在实际应用中帆板处于非平衡状态，因此表面温度使用以下简化公式进行求解：cmdTdt=AsT4+Q,（12）cJ/(kgK)mAsQ式中，为比热容，取值为 350，是帆板质量，取值为 1.5kg，为帆板面积，此处取值为0.3，为帆板表面总的吸收热流，表面吸收率设为 0.9。采用蒙特卡洛方法获得表面外热流后，代入公式迭代求解帆板温度，并与卫星遥测数据做对比，结果如图 14 所示。可以看出，计算结果在各个轨道周期内的变化幅值和规律与遥测数据符合较好。最高温度约为 95C，最低值约为76C，在光照区升高，在地影区下降，不断循环。本文计算结果可以为太阳帆板温度场分析提供可信参考。同时也可以看出，计算值和遥测值在部分区域存在一定偏差，可能是参数设置与实际情况的差异影响了计算结果，需要开展进一步研究。010080604020020406080100120140Remote measurementThis work10 000Time/s20 00030 000Temperature/图14太阳帆板温度对比Fig.14Temperaturecomparisonofsolarpanel3.5红外辐射强度分布目标自身红外辐射由表面温度和表面材料的发射率决定，在获得温度变化曲线后，目标自身红外光谱辐射出射度可由式（13）表示：Mself,i=i()c15e(c2/(Ti)1,（13）c1c2Tii()其中，为第一辐射常量，值为 3.7421016Wm2，为第二辐射常量，其值为 1.4388102mK，为波长，单位为 m，是面元表面温度。是面元光谱发射率。结合图 12 和图 13 给出的两种模式下的各表面温度，认为各表面为漫反射，根据朗伯余弦定律计算得到两种模式下一个轨道周期内的各方向上的辐射强度，如图 15图 18（彩图见期刊电子版）所示。000.51.01.52.02 0001 0003 0005 000Time/s4 0006 000Radiation intensity/(Wsr1)+XX+YY+ZZ图15对地模式下的各方向辐射强度（35m）Fig.15Radiationintensityineachdirectionintheearth-pointingmode(35m)194中国光学（中英文）第17卷可以看出，在常用的 35m 和 814m 两个探测波段中，两种模式下+Z 和Z 方向的辐射强度较高，35m 波段最大值在 1.5Wsr1左右。814m波段最大值在 22Wsr1左右。主要原因是帆板温度较高且面积较大，是红外信号的主要来源。814m 波段辐射强度明显强于 35m 波段。同一谱段内，对地模式下，不同时间内，+X、X 和Y 方向也存在红外信号较强时刻。但对日模式下，除+Z 和Z 方向的辐射强度较高外，其他方向红外辐射值较小，目标探测存在困难。000.51.01.52.02 0001 0003 0005 000Time/s4 0006 000Radiation intensity/(Wsr1)+XX+YY+ZZ图16对日模式下的各方向辐射强度（35m）Fig.16Radiation intensity in each direction in the sun-pointingmode(35m)00510152025302 0001 0003 0005 000Time/s4 0006 000Radiation intensity/(Wsr1)+XX+YY+ZZ图17对地模式下的各方向辐射强度（814m）Fig.17Radiationintensityineachdirectionintheearth-pointingmode(814m)4结论本文针对绕地卫星轨道外热流，采用蒙特卡洛法、基于非结构四面体网格和 OpenMP 并行算法编写了仿真软件。针对文献中的典型工况进行了对比分析，验证了软件计算结果的准确性。进一步，针对 535km 太阳同步轨道，考虑卫星结构遮挡的影响，对不同姿态控制策略下的卫星轨道外热流进行了仿真分析，并将太阳帆板表面温度与在轨遥测数据进行了对比分析。得出如下结论：（1）不同姿态模式下目标的轨道外热流区别较大。对于 535km 太阳同步，10:30 地方时轨道来说，对地模式下除+Y 面外，各表面热流值随时间变化波动较大，而对日模式下除Z 面本体及帆板表面波动较大外，其他表面变化较小。（2）蒙特卡洛算法对空间目标复杂结构及结构遮挡具有很好的适应性，对地模式下，考虑遮挡后，Y 表面热流值降低了 53.79W/m2，+YZ 侧帆板表面热流值降低了 32.05W/m2。（3）不同模式下目标各表面温度特性不同。对地模式下各表面温度随时间波动较大，使红外观测窗口规划提高了难度。在两种模式下，帆板在光照区温度较高，具有明显的红外特征，便于开展空间目标红外观测。参考文献：谷牧,任栖锋,周金梅,等.基于地基观测的时序卫星红外光谱建模与分析J.物理学报，2019，68（5）：059501.GUM,RENQF,ZHOUJM,et al.Modelingandanalyzingoftime-resolvedsatelliteinfraredspectrumbasedon100510152025302 0001 0003 0005 000Time/s4 0006 000Radiation intensity/(Wsr1)+XX+YY+ZZ图18对日模式下的各方向辐射强度（814m）Fig.18Radiation intensity in each direction in the sun-pointingmode(814m)第1期郑鸿儒,等:空间目标轨道外热流计算及辐射特性研究195ground-baseddetectorJ.Acta Physica Sinica,2019,68（5）:059501.（inChinese）郑鸿儒,马岩,范林东,等.高空离轨发动机流场红外辐射特性研究J.中国光学，2022，15（2）：260-267.ZHENGHR,MAY,FANLD,et al.Infraredradiationcharacteristicsofhigh-altitudeoff-orbitengineplumeJ.Chinese Optics,2022,15（2）:260-267.（inChinese）2易桦,黄兴,江海.一种圆轨道航天器外热流通用计算方法J.航天器工程，2021，30（5）：53-58.YIH,HUANGX,JIANGH.CommoncalculationmethodfororbitalheatfluxofspacecraftoncircularorbitJ.Spacecraft Engineering,2021,30（5）:53-58.（inChinese）3GARZN A,TAMI J A,CAMPOS-JULCA C D,et al.Effect of beta angle and contact conductances on thetemperaturedistributionofa3UCubeSatJ.Thermal Science and Engineering Progress,2022,29:101183.4李志松,马昌健,毛云杰,等.微纳卫星在轨温度场快速分析J.航天器环境工程，2021，38（2）：122-129.LIZHS,MACHJ,MAOYJ,et al.RapidanalysisoftemperaturefieldfororbitingnanosatellitesJ.SpacecraftEnvironment Engineering,2021,38（2）:122-129.（inChinese）5李世俊,陈立恒,冯文田,等.太阳同步轨道二维变姿态空间相机的外热流计算J.红外与激光工程，2018，47（9）：0917008.LISHJ,CHENLH,FENGWT,et al.Calculationofexternalheatfluxesonspacecamerawithtwo-dimensionalchangingattitudesinsun-synchronousorbitJ.Infrared and Laser Engineering,2018,47（9）:0917008.（inChinese）6吴愉华,陈立恒,李行,等.地球静止轨道变姿态空间相机的外热流计算J.红外与激光工程，2019，48（6）：0604001.WUYH,CHENLH,LIH,et al.ComputationofexternalheatfluxesonspacecamerawithattitudechangeingeostationaryorbitJ.Infrared and Laser Engineering,2019,48（6）:0604001.（inChinese）7ATARC,AKTAM.Advancesinthermalmodelinganda