可重复使用运载器热防护系统热/力耦合数值计算研究

论文摘要

热防护系统(Thermal Protection System—TPS)设计与分析是可重复使用航天器(Reusable Launch Vehicle—RLV)的重要关键技术之一，开展热防护系统研究，确定其结构参数对RLV质量与温度场分布的影响作用、再入时热冲击导致的热力耦合效应等是热防护系统设计分析的重要研究课题。 基于有限元Galerkin方法原理，本文细致研究了线性瞬态热传导方程和非线性瞬态热传导方程的数值计算方法，进行了数值分析计算，主要研究工作包括，推导了正交各向异性的复合材料等参元有限元列式，给出了正交各向异性材料的热传导系数从材料主轴、整体坐标系到参数正则坐标系下的有限元等参元计算变换式；为计及热载荷具有热冲击性质，温度场在单元间采用二次插值函数，并推导出不同边界条件下的单元算法和有限元数值计算列式，获得到各向异性层状材料热分析超参数壳体单元的有限元完整计算方法；在此基础上，突破传统单元构造方法，提出了一种中面层8节点高精度等参三维退化热壳单元计算列式，使壳结构的热特性数值分析变量减少了一倍，通过数值算例验证了该算法的正确性和高效性；根据相对自由度概念，构造了相对自由度壳体单元的热分析算法，通过具体算例与ANSYS软件计算结果的比较，吻合很好，验证了本文算法的正确性；与此同时，把使用超参数温度壳元计算结果和相对自由度温度壳元计算结果相比较，表明了超参数温度壳元和相对自由度温度壳元有良好的一致性。 针对热防护系统金属蜂窝典型单元，本文构造了蜂窝核细观导热有限元数值计算分析模型，建立了空间二维薄板热分析有限元算法，推导了一般性计算公式，建立了空间四节点薄板等参元计算格式；在此基础上，讨论了典型蜂窝单元的数值计算模型和边界条件处理算法，确定了具体结构尺寸、边界条件以及热流函数曲线。鉴于典型蜂窝单元材料物性参数的复杂性，在没有具体的物性参数可用的情况下，参考国外技术资料给出了计算热物性经验公式。在大量数值计算基础上，研究分析了结构尺寸和不同辐射条件对典型蜂窝单元的温度场分布影响作用，获得了结构构形参数的敏感性作用。在典型蜂窝单元的传热分析基础上，本文分析推导了准静态热力耦合有限元列式，采用ANSYS软件

论文目录

• 摘要
• Abstract
• 目录
• 第一章 绪论
• 1.1 可重复使用热防护系统概述
• 1.1.1 国外RLV的热防护系统发展特点以及形式
• 1.1.2 国内外RLV的热防护系统结构和材料形式
• 1.2 可重复使用热防护系统热分析问题的研究现状
• 1.2.1 国外金属热防护系统热分析现状与新进展
• 1.2.2 国内热防护系统热分析现状与新进展
• 1.3 可重复使用热防护系统热力耦合问题的研究现状
• 1.4 瞬态系统的时程精细积分方法
• 1.5 本文主要工作
• 第二章 热传导方程的自适应算法研究
• 2.1 引言
• 2.2 热传导方程有限元计算的Galerkin格式
• 2.2.1 热传导微分控制方程
• 2.2.2 有限元计算的Galerkin格式
• 2.3 二维热传导方程的自适应算法研究
• 2.3.1 坐标变换
• 2.3.2 温度插值函数
• 2.3.3 温度有限元列式
• 2.3.4 边界条件的应用
• 2.3.5 二维四边形阶谱单元
• 2.4 三维八节点各向异性体单元
• 2.5 线性有限元瞬态导热微分方程组的数值解法
• 2.6 非线性有限元瞬态导热微分方程组的数值解法
• 2.7 算例分析
• 2.8 小结
• 第三章 各向异性层状材料的热壳元分析
• 3.1 引言
• 3.2 超参数温度壳体单元
• 3.2.1 壳体几何形状的规定
• 3.2.2 参数坐标系下的温度函数插值
• 3.2.3 有限元列式
• 3.3 层壳复合材料单元的积分简化汁算
• 3.4 相对自由度温度壳体单元
• 3.5 算例分析
• 3.6 小结
• 第四章 金属热防护系统蜂窝结构的热分析
• 4.1 引言
• 4.2 蜂窝核微观导热有限元分析模型
• 4.2.1 蜂窝结构示意图
• 4.2.2 三维金属薄壳的有限元数值分析方法
• 4.2.3 蜂窝平面四节点有限元列式
• 4.3 蜂窝型热防护板数值计算模型以及边界条件处理
• 4.3.1 建立模型
• 4.3.2 物质热属性的确定
• 4.3.3 确定边界条件和载荷
• 4.4 蜂窝型金属热防护板热分析
• 4.4.1 结构尺度对于温度场分布的影响
• 4.4.2 辐射对于温度场分布的影响
• 4.5 小结
• 第五章 热力耦合有限元数值算法研究
• 5.1 引言
• 5.2 热力耦合有限元列式
• 5.2.1 热应力与热应变
• 5.2.2 热力平衡微分方程
• 5.2.3 热力耦合有限元列式
• 5.3 蜂窝型金属热防护板二维热力耦合分析
• 5.3.1 单层蜂窝热防护板二维热力耦合分
• 5.3.2 双层蜂窝热防护板二维热力耦合分析
• 5.3.3 两块热防护板的二维热力耦合分析
• 5.4 温度场与应变场的耦合数值计算列式
• 5.5 小结
• 第六章 热传导方程的时程精细积分方法
• 6.1 引言
• 6.2 线性问题的时程精细积分法
• 6.3 非线性问题的时程精细积分法
• 6.4 瞬态热传导分析的时程精细积分法
• 6.5 非线性瞬态热传导分析的时程精细积分法
• 6.6 算例分析
• 6.7 小结
• 第七章 热力耦合系统的精细积分研究
• 7.1 引言
• 7.2 能量方程
• 7.3 增维精细积分方法
• 7.3.1 非齐次项为方程组的解
• 7.3.2 增加常数维
• 7.4 准静态热力耦合
• 7.5 动态热力耦合
• 7.6 算例分析
• 7.7 小结
• 第八章 总结与展望
• 8.1 工作总结
• 8.2 工作创新点
• 8.3 工作展望
• 参考文献
• 攻读博士学位期间发表的相关学术论文
• 致谢