金属结构航天器陨落过程三维瞬态传热有限元算法研究

来源：计算力学学报杂志2019年第2期北京时间：

作者：石卫波,孙海浩,唐小伟,马强,李志辉

单位：中国空气动力研究与发展中心, 绵阳 621000,中国空气动力研究与发展中心, 绵阳 621000,中国空气动力研究与发展中心, 绵阳 621000,四川大学 数学学院, 成都 610043,中国空气动力研究与发展中心, 绵阳 621000,国家计算流体力学实验室, 北京 100191

摘要：构建金属桁架结构航天器陨落再入气动热环境有限元传热模型，是准确预测在轨服役期满大型航天器陨落再入解体过程温度分布的关键。本文采用四节点四面体单元对空间进行离散，依据泛函理论，将传热控制方程离散为代数方程组；利用有限单元法总体合成得到具有对称正定、高度稀疏和非0元素分布的规则性刚度矩阵，发展一维变带宽压缩存贮技术，有效解决大型稀疏矩阵的数据存贮问题；为有效抑制求解过程出现的温度在时间和空间上的振荡问题，发展集中热容矩阵系数处理方法，将热容矩阵的同行或同列元素相加代替对角线元素，使非对角线元素化为0，构造求解三维瞬态温度场的两点向后差分格式、Crank-Nicolson格式和Galerkin格式。通过对正方体瞬态传热计算验证分析，在相同条件下，采用以上三种格式均可获得一致稳定的温度解，并得到与现有ANSYS有限元软件较为吻合的计算结果，验证了所建立三维瞬态传热有限元计算模型的准确可靠性。在此基础上，对铝合金低轨航天器薄壳结构进行了传热计算，给出了类天宫飞行器两舱体陨落飞行107.5 km~90 km不同高度的瞬态温度分布，为这类寿命末期航天器陨落再入解体预报提供理论支撑与可计算模型。

关键词：金属桁架结构航天器,陨落再入,三维瞬态传热,有限元计算模型,

基金资助：国家重点基础研究发展计划（2014CB744100）；国家自然科学基金（11325212；91530319）资助项目.