硅基材料烧蚀模型研究

对冲压发动机燃烧室所使用的硅基热防护材料，用以往的单纯因二氧化硅液态层被气动吹除减薄的模型进行烧蚀计算，计算结果和实际偏差很大。我们根据硅基热防护材料的实际特性，在原有的液态层模型的基础上，考虑了夹杂碳化层对熔融态硅基的补强作用，提出了既有液态层吹除也有化学反应烧蚀的新模型，即考虑了热解、熔融和碳化、气动吹除结合化学反应烧蚀的多因素过程。并根据硅基热防护材料的烧蚀特点，根据发动机模型的试验结果，初步总结出了液态层粘附力的经验公式。从物理概念上分析可以看出该模型基本符合实际情况；从多种工况计算结果可看出，新模型和经验公式是比较合理的。     

热容燃烧室基于单点温度的热流测量方法研究 *

为了研究推力室内壁面与燃气间的传热特性,获得多喷嘴燃烧室内壁热载分布,从而为推力室设计及热结构分析提供参考,开展了热容燃烧室壁面热流测量方法研究。以位于距内壁一定距离位置的单点测量温度作为输入条件,文章提出了三种热流计算方法,并进行火箭发动机试验研究,得到了对应位置燃气对内壁的热流密度。可以看出,运用单点方法得到的热流密度与传统两点法结果符合较好,两者计算结果误差在10%以内。结果显示:不同点火时长工况下在同一位置相同时间点上重复性很好,在燃烧室圆柱段热流密度随着点火时长的增加而减小,而在喷管段变化规律则相反,热流密度随时间增加而增加。在提出的计算平均热流的三种方法中,热积累方法的计算结果最高,热平衡方法次之,瞬态法最低。在火箭发动机热容试验中,本文提出的方法可以用于测量燃烧室内壁面热流密度。     

再生冷却推力室热障涂层系统的热结构分析

为了设计和优化适用于液体火箭发动机推力室的热障涂层,应用ANSYS的热-结构分析功能,对再生冷却推力室-热障涂层系统进行了热结构有限元分析,得到在不同涂层覆盖下,推力室壁中的温度场和应变场,并通过对热障涂层中应变场的分析,研究不同涂层发生分层剥落的关键位置以及主要驱动力。结果表明,陶瓷层厚度较大的YSZ+Ni Cr Al Y涂层拥有更优异的性能,使推力室壁在热试阶段的最大应变量减少约36.1%;工作循环中,涂层与推力室壁的接触面上会产生较大的应变量,最终有可能导致涂层剥落失效;粘结层能缓解涂层与推力室壁间的热膨胀系数不匹配,使陶瓷层在热试阶段的最大应变量减少约80%。     

基于细观颗粒夹杂模型的复合固体推进剂松弛模量预测

为更准确地预测不同固体颗粒体积分数的复合固体推进剂的松弛模量,采用了分子动力学方法对不同体积分数的复合固体推进剂细观模型进行建模.根据有限元理论及细观力学均匀化方法,计算在定应变工况下复合固体推进剂细观模型的平均应力随时间的变化,从而有效地预测复合固体推进剂的松弛模量.该方法有效地体现了随填充颗粒体积分数的增大,复合固体推进剂瞬时模量逐渐增大的变化规律及颗粒随机分布对复合固体推进剂瞬时模量的影响.将其应用到复合固体推进剂的设计过程中,可有效降低设计成本,缩短设计周期.      

原子氧和紫外线对聚酰亚胺综合作用数值模拟

在建立数学物理模型的基础上,对低地球轨道环境和地面试验环境下有无保护涂层的聚酰亚胺所受原子氧冲蚀及紫外线的综合作用进行了数值模拟,获得了具有工程应用价值的计算结果,并讨论了数学物理模型中各参数对基蚀曲线形状的影响.从数值模拟结果与美国太空试验结果的比较可以看出,得出的数值模拟的结果是正确的,对航天器设计具有重要的指导意义.     

固体火箭发动机复合结构喷管温度场与应力场理论计算

对某长时间工作的固体火箭发动机复合结构喷管的二维轴对称瞬态温度场及二维轴对称准静态应力场进行了理论预估。温度场的计算采用时间域上的有限差分及空间域上的有限元相结合的方法。应力场的计算采用以位移作为未知量的有限元法。计算中考虑了材料的方向性及物性参数随温度的变化，采用变带宽一维压缩存贮刚度矩阵，节约了大量的内存。计算结果与试验结果吻合较好。     

基于Stereo Lithography文件和有限元法的在轨航天器关键部件热计算

针对STL(Stereo Lithography)文件在传递复杂几何实体模型信息方面具有精度高的特点,提出了一种基于STL文件计算复杂结构角系数和外热流的方法。给出了有限元网格离散和边界识别方法。根据单元和节点的拓扑关系进行辐射换热计算。利用有限元法计算辐射换热角系数和外热流。最后利用有限元法计算航天器关键部件在轨的三维瞬态温度场,并进行了计算效率和误差分析。     

液膜再生复合冷却中液膜传热特性

对液体火箭发动机液膜再生复合冷却进行了算法研究.综合考虑了发动机内部化学反应、蒸发、卷吸、对流、导热、辐射等因素,将冷却液膜分为显热区、潜热区及气膜区三个区域进行了计算.推导了液膜长度和厚度的计算方法,分析了液膜再生复合冷却效率及各因素对液膜传热特性的影响.计算结果表明:①液膜入口质量流量越大,液膜区长度越长,冷却效率越高,复合冷却效率可维持在0.57以上.②高温燃烧室内膜的液体段长度很短,在液膜存在区域内冷却效率高达0.9.③液膜消失后,头部冷却液膜的设计仍对室壁起了很好的冷却保护作用,低温边区一直延伸至出口.④液膜吸收的显热和液膜蒸发吸收的热量及高温燃气与膜间的对流在液膜区内起了主要作用,而卷吸造成的质量损失及传热不可忽略.      

固体火箭发动机药柱三维结构非线性分析

非线性分析是药柱结构分析中的难点，针对其对结构分析的影响，基于不可压缩材料的粘弹性本构关系，应用完全拉格朗日（T．L）法的虚功方程，综合考虑药柱的近似不可压缩性和几何非线性，推导了三维粘弹性几何非线性有限元增量方程，编写了有限元程序对星型药柱在受压力载荷以及固化降温载荷作用下的结构进行了分析，并与线性计算的结果进行了对比，结果表明，在大载荷作用下，非线性对药柱结构分析的影响比较显著，计算时考虑非线性的影响是必要的。     

三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。