航天器在轨全周期热变形分析.docx

航天器在轨全周期热变形分析提出一种适应于在轨全周期热变形的分析方法，采用基于热传导算法进行“热分析模型构造分析模型全周期温度场映射，利用数学拟合算法开展对各类结果数据的分析，通过相关程序实现全周期多工况温度场映射、计算、数据分析的自动化。对某遥感卫星进行全周期热变形分析，结果表明：全周期温度场映射时间由天缩短至小时量级，温度场映射精度可控制在以内，相对于以往基于极端工况的热变形分析方法，可显著地提升分析精度与验证覆盖性，获得在轨热变形量级、全周期变化规律。文章的研究结果可为航天器热稳定设计提供参考。关键词：航天器；在轨全周期；热变形；稳定性随着对地、对天观测航天器指标要求的日益提高，高图像定位精度成为高性能遥感航天器的典型特征。在轨构造变形直接影响相机、星敏感器、陀螺等关键部件本身空间指向及相互间的几何关系，甚至影响相机内部各镜片间的空间位置关系，是决定图像定位精度、相机成像质量的重要因素之一。构造在轨热变形在相机安装处引起的位移一般为微米级，对于低分辨率观测航天器，这些扰动能够忽略；但对分辨率优于的航天器，则必须考虑热变形扰动影响，由于这些影响可能直接决定了成像质量能否到达设计指标。一般情况下，微米级的构造变形可能导致角秒级的设备安装面法向指向变化，进而出现米级的成像误差。同时，与其他因素相比，航天器构造在轨热变形具有一定的随机性，很难通过后期在轨处理消除其影响。因而，在轨稳定性对高分辨率航天器的性能指标至关重要，在地面研制阶段就应结合航天器系统需求开展航天器高稳定构造设计、验证工作。仿真分析是航天器机械系统研制的有力支撑，对于设计工况复杂、影响因素诸多的构造热稳定设计而言，在轨热变形分析直接支撑了系统指标分配、构造研制、热变形地面试验方案设计等各研制环节，国内已经开展的构造热变形分析工作多数是针对构造热膨胀、吸湿性的理论研究，以及零膨胀铺层设计的研究等，对航天器构造的在轨热变形分析及试验验证主要立足于某特定温度场或模拟温度场，如模拟在轨工况的最高温工况、最低温工况、最大温差工况等，或通过施加最大包络载荷实现对在轨热变形的预估及构造低膨胀设计。随着构造热稳定性要求的日益提高，国外对构造热稳定性的研究已经深化到在轨微裂纹、微蠕变等领域，并开展了构造热稳定高保真仿真方法研究、影响因素灵敏度研究等，这些研究的基础在于对在轨全周期热变形的高效、高精度仿真。在轨全周期热变形分析还能为在轨成像标定策略的制定提供参考，进而实如今传统图像修正方法的基础上进一步提升修正精度。因而，新一代遥感航天器研制对在轨全周期热变形分析需求日益迫切。由于在轨全周期热变形分析涉及的工况数量较大，一般可达数百、甚至上千个温度工况，且面临机热耦合效应复杂、温度边界及力边界模拟难度大等技术瓶颈，因而对热控设计温度场数据与构造分析温度场数据高效、高精度映射，以及多工况下热变形高保真分析、结果高效处理及判读等，均提出了更高要求。在的部分设计资料里，固然提到了热变形分析方法，但没有阐述怎样开展在轨全周期各个时刻的热变形分析。目前，国内关于全周期热变形分析的案例较少，支撑全周期热变形分析、实现海量数据快速映射的方法更是鲜有提及；固然有学者根据热分析模型重新划分构造分析模型，然后从数据文件中读取相应节点的各个时刻的温度数值，根据不同载荷工况的形式写入计算文件进行温度场分析，但是数据转换经过较多，影响了全周期热变形分析的效率。本文首先对适应于在轨全周期热变形分析方法进行讨论，重点介绍了基于热传导算法的“热分析模型构造分析模型全周期温度场映射方法，以及基于数学拟合算法开展结果数据处理及判读方法，阐述了全周期热变形分析流程。基于上述方法对某遥感卫星进行全周期热变形算例分析，获取了在轨热变形量级、全周期变化规律等，并与传统分析方法进行了比照，结果表明本文提出的方法可显著地提升分析精度与验证覆盖性。全周期热变形分析方法总体思路用于热变形分析输入的温度场，通常是基于、等软件开展的热分析得到，用于热分析的数学模型与用于构造分析的数学模型一般为通过、等有限元分析软件建立的构造分析有限元模型在节点位置、网格离散程度、建模简化方式等方面均存在差异性。例如，某航天器高稳定载荷适配支撑构造热分析模型约有个节点，构造分析模型有个节点。因热分析与构造分析所采用的软件差异性，以及二者分析模型的差异性，在开展热变形分析前，首先要将热分析温度场映射至构造分析模型上，进而实现温度数据从热分析模型传递至构造分析模型，并作为构造分析输入载荷。本文提出的方法是首先生成无温度场的构造分析计算文件，然后进行全周期工况判读及分析，基于热传导算法实现热分析模型向构造分析模型的温度场映射，并对可能存在的不能映射节点和奇异节点温度远高于或远低于在轨实际温度的节点进行二次映射，继而生成映射后的温度场及构造分析计算文件，判定无误后进行计算，并对全周期热变形分析结果进行数据拟合、生成报告。详细流程如图所示。基于热传导算法的全周期温度场映射热传导算法概述从映射算法上讲，以往热变形分析主要采用基于几何差值的映射算法，该算法仅与空间位置相关，是目前商业软件中广为采用的映射方法。其局限性在于，对于非连续构造、不同组件连接构造的映射工况，易产生映射奇异的现象，无法识别各部位或各组件之间温度的差异性。从温度场映射实现工具来讲，目前广泛采用基于有限元商业软件进行温度场映射，此类方法对于温度场单次映射较为通用，但对于全周期热变形分析则具有一定的局限性，详细表现为：温度场导入及映射功能主要基于手动实现，很难知足全周期成百上千个工况温度场高效映射分析。商业软件一般仅内嵌基于几何算法的映射方法。因单位不同、设计师不同，热分析经过可采用、等不同软件，进而导致热分析结果数据格式存在显著差异，商业软件在数据导入、映射方式上对于各类热分析软件适应性较差。本文提出的基于热传导算法的全周期温度场映射方法，能以既有节点温度场为基础，根据构造热传导特性进行映射计算，避免因个组件空间距离较近、但并不属于同一温度范围的节点发生映射关系。同时，采用二次开发程序对映射方法进行封装，分别建立适应于不同热分析软件的温度场映射模块。通过开发与有限元商业软件前后处理工具的接口，实现温度场映射批量处理、结果数据批量处理、结果数据批量判读等，可显著地提升映射精度和分析效率。通过图一个简单的构造组件分析模型，能够证实基于热传导算法进行温度场映射的优势。算例中的构造组件由高温构造、低温构造、室温构造部分组成，相互间存在隔热层，能够阻断热传递。基于几何插值映射方法无法考虑低热导率层的影响，见图；而基于热传导算法得到的温度场映射结果，可更为真实地反映实际温度情况，见图。温度场映射经过基于热传导算法的全周期温度场映射经过包括个步骤。构建热分析模型单元节点和构造分析模型单元节点之间的对应关系，如图所示。基于数学方法实现热分析模型节点与构造分析模型节点的关联，一般准则为：覆盖热分析模型节点的构造分析模型节点温度，根据热分析模型节点温度取值。使用有限元形函数获取加权系数。式中：为第个热分析模型节点温度；为第个构造分析模型节点温度。写成矩阵形式为式中：为热分析模型节点温度矩阵；为权重系数矩阵；为构造分析模型节点温度矩阵。通过用热控材料替代构造分析模型材料例如用材料卡片替换材料卡片，基于构造分析模型可计算得到热传导矩阵。求解如下的插值方程。式中：为拉格朗日乘子。通过式即可求解构造分析模型节点温度矩阵。数据拟合算法航天器在轨热变形分析的目的，是获取关键设备指向变化或各设备间的夹角变化，由此引申出采用何种方式来表征设备指向及其夹角的问题。目前，国内外广泛采用的表征方式有种。对于光学相机、星敏感器等设备，主镜、次镜等关键部件均沿设备轴向且近似在一条直线上，可选取此线上的多个关键点，应用“多点拟合线的方式获取设备指向。选取设备，安装面上的多个关键点，采用多点拟合面的形式获取设备安装面矢量，以此模拟设备安装指向，此种方法对于各类设备均具有通用性。由于构造分析结果一般为有限元模型节点位移，因而要借助其他程序并选取相应数学算法，对有限元分析结果进行二次处理和判读。对于直线矢量计算，设待拟合直线矢量个节点的坐标为，写成如下矩阵形式。熿燀燄燅计算式的协方差矩阵如下。式中：式中：为平均值矩阵；，；，。计算协方差矩阵的特征值为，则由特征值，组成的向量即为待拟合直线的向量。对于平面法线向量计算，设待拟合平面法线矢量的个节点的坐标为，写成如下矩阵形式。熿燀燄燅由式各列减去各自的均值，得到矩阵如下。计算矩阵的特征值为，则由特征值，组成的向量即为待拟合平面法线矢量。算例分析温度场映射某遥感卫星个星敏感器通过支架安装于相机承力框上，相机安装在卫星构造平台上，其中整星机械坐标系向为卫星纵向相机对地观测方向，整星机械坐标系向、向为卫星横向相机承力框及载荷适配构造面内方向。该卫星天天运行个轨道周期，热控设计时要对卫星全生命周期所有极端工况取一个最大包络，即卫星在轨运行天天承受的温度工况均不会超过目前给定的个轨道周期状态，以载荷适配构造为例，其中个典型位置的个轨道周期热分析节点温度见图。根据高定位精度设计需求，须开展个轨道周期不同姿态下相机成像、数传记录等关键时刻点的热变形分析，由此获取相机安装面法向转角、星敏感器安装面法向转角，以及相机安装面法向与星敏感器安装面法向间的夹角变化。同时，分析不同的姿态、星敏感器和工作形式下上述各项分析结果的变化规律，为高定位精度指标分析提供支撑。将采用软件得到的全周期热分析温度场作为输入，采用图仿真流程、基于热传导算法进行全周期近千余时间点温度场映射，并对映射奇异节点进行二次修正，生成可用于软件进行有限元分析的批处理求解文件。此外，通过程序实现“热分析温度场输入温度场映射温度场修正有限元计算高度集成化与自动化，进而实现全周期温度场映射时间由天缩短至小时量级。为验证映射准确性，选取第一轨道周期温度梯度较大个时刻点，分别对应第一轨道周期温度场“正弦曲线的起点时刻、波谷时刻、波峰时刻、终点时刻。将映射前的热分析温度场与映射后的构造分析有限元分析温度场进行比照，由图热分析温度场与构造分析温度场比照结果能够看出，基于热传导算法可实现热分析与构造分析温度场的准确匹配，由热分析模型到构造分析模型的温度场映射精度可控制在以内详见表。全周期热变形分析结果在以往热变形评估工作中，主要选取模拟在轨工况的极端工况进行分析，尤其是以某个高温或低温状态的均匀温度作为输入载荷，基于上述思路，针对图所示第一轨道周期典型时刻温度场，选取个极端温度点，作为均匀温度载荷输入模拟工况，开展热变形分析，获取星敏感器安装面法向与相机安装面法向夹角变化，并与全周期分析结果真实工况进行比照，见图。从图能够看出：个模拟工况固然能够反映一个周期内变形的平均值，但不能覆盖全周期各时刻点可能出现的真实变形情况，讲明采用全周期热变形分析的必要性和优势。在轨运行期间，卫星有效载荷如光学遥感卫星的对地或对天观测相机部分会采取精细控温、相对常温变化仅为几摄氏度，而卫星平台部分相对常温存在几十摄氏度的温度波动。选取第一轨道周期变形量最大时刻点进行变形分析，由图整星、平台、主承力立柱、相机适配支撑构造等部分变形云图能够看出，在轨温度交变引起的平台部分变形在几百微米，接近毫米级。而由图、图相机、相机主承力框变形云图能够看出，相对于平台部分而言，相机主承力框变形相对较小。通过全周期热变形分析，获取了轨道周期个星敏感器安装面与相机安装面法向夹角变化，见图。由分析结果能够看出：个星敏感器安装面与相机安装面法向夹角变化呈现正弦周期性变化，每一圈对应一个完好正弦波，此种变化趋势与图所示的热分析温度场周期性变化情况相对应。星敏感器安装面与相机安装面法向夹角变化存在显著差异性，如星敏感器安装面与相机安装面法向全周期最大夹角变化超过，而另外个星敏感器安装面与相机安装面法向最大夹角变化未超过。个星敏感器支架同一变形形态出现时刻存在差异，即一个星敏感器支架变形位于“波峰之时，另一个可能位于“波谷。由此，获得了不同星敏感器安装面与相机安装面法向最大夹角变化量级及相互差异性，可为整星在轨稳定性评估、卫星定位精度评估等提供重要参考。本文提出了在轨全周期热变形分析方法，基于热传导算法进行“热分析模型构造分析模型温度场映射，基于数学拟合算法开展结果数据处理，获得了航天器在轨全周期热变形量级、全周期变化规律等，可得出如下结论。基于热传导算法并通过相关程序实现全周期千余时间点映射经过的自动化，可将全周期温度场映射时间由天缩短至小时，温度场映射精度可控制在以内，相对于传统分析方法显著地提升了分析精度和全周期覆盖性。获取了整星在轨全周期热变形情况，从变形云图能够得到整星平台、相机、星敏感器支架等关键部位变形状态、量级及宏观变形传递趋势。通过各星敏感器安装面与相机安装面夹角变化全周期分析，能够看出个星敏感器安装面与相机安装面夹角变化呈现正弦周期性变化，且变形量级、同一变形形态出现时刻均存在差异。本文的热变形分析方法及所获取的变形规律具有一定的普适性，可用于航天器机械系统的热稳定设计，对于高精度、高稳定性航天器的研制具有参考价值。