喷管喉衬温度场及热应力计算分析研究.pdf

喷管喉衬温度场及热应力计算分析 李耿周为民张钢锤 (航天科技集团公司第四研究院4l所。西安。710025) 摘要：在已知喉材烧蚀率的每件下，运用有限元分析软件^NS丫s的生死单元技术。考虑工 作过程中喉村的型面变化。对莱固体火箭发动机喉种几个不同时刻的温度场和热应力进行了耦合计 算，并矸结果进行分析． 主题词：固体火箭发动机，喉衬．生死单元 l引言 喉衬是固体火箭发动机唼管的关键部件，它构成了气流通道的最小截面，控制燃烧室压强，维持 药柱正常燃烧，使燃气从亚音速加速到超音速。它的性能直接影响发动机性能及其工作可靠性。由于 喉村较长时间工作在3000k以上的高温环境中，既承受者高温热载．困此对其热环境的分析一直是人 们比较关注的。由于当前结合烧蚀和传热计算喉衬热结构的手段还不成熟。计算中喉衬始终保持初始 状态，不能模拟出烧蚀对几何构型以及传热的影响．为了克服这一缺陷，本文在己知喉衬烧蚀率的条 件下，考虑工作过程中烧蚀引起的喉衬型面变化，对几个不同时刻喉挂的温度场和热应力进行了计算， ’ 并对结果进行分析。 2计算模型 2．1模型的简化 喷管喉衬的几何形状以及工作中的载荷都是轴对称的，可以认为热传导只发生在轴向和径向．而 无周向热流，因而可将实际上的三维空间问题简化为轴对称问题．并在分析求解时作如下假设： (1)忽略两层材科之间的接触热阻； (2)燃烧室温度视为常数，初始温度为室温i (3)喷管内流场的计算采用一维等熵流； (4)喉衬的热学性能参数是随温度而变化的，力学性能参数为常温下的数据： 经过上述的简化处理．实际的结构被简化为如图1所示的模型 2．2边界条件和初始条件 边界条件为： 对漉边界 绝热边界 竺锄栏a 初始条件为： 矿卜，鱼忭霉 每《一 2．3喷管内流场的计算 假设燃烧产物在喷管中的流动属于一维等熵过程，咀速度系数九为自变量．其控制方程如下： r广 m…数椰，每(州斟丑[t一丽k-I 啪刺以，=；小等∥ +I 砌，=寺=(t一等kz2尸 p ＼ ／ (3)压强函数 五= 肘：兰 (4)速度系数 口 V ． ％= ^=一 ％ (5)临界声速 式中 V，P，丁，吖，口分别为半径R的计算截面内的速度。压强，温度．马赫数和音速 如为最小截面半径ir，P‘为滞止温度和压强；K，R为燃气比热比和气体常效． 根据以上公式即可算出喷管内任一截面的速度．压强，温度和马赫数。 2．4对流换热系数的计算 喷管壁面与燃气的换热系数采用半经验的Bartz公式， 公式如下； 矗=击c等等x等)o。c钞睁”盯 式中，Ml，嗍．026；M=l，b=0．0245，MI：嘲．023 Pr： 墨 (1)普朗特数 1．77k一0．45； 式中k为燃气比热比。 0·87×1