小推力长工作时间固体火箭发动机喷管温度仿真分析.pdf

第 卷第 期 固 体 火 箭 技 术 小推力长工作时间固体火箭发动机喷管温度仿真分析邓康清，王鹍鹏，余小波，向 进，朱雯娟，王相宇，杨育文，汤 亮，张琪敏（航天化学动力技术重点实验室，襄阳；湖北航天化学技术研究所，襄阳）摘要：利用 发展了一种瞬态流固热一体化一步耦合三维有限元分析方法，并用该方法分析了一种采用无铝 推进剂的小推力长工作时间（）固体火箭发动机喷管温度场分布情况。分析结果表明，从喷管中心到喷管表面温度逐渐降低，喷管喉部中心温度达到 以上，喉衬温度为 ，隔热层温度为 ，壳体温度为 ；喷管中心轴线上各点的温度开始时急剧增加，后基本保持不变，从喷管出口到喷管入口，沿喷管中心轴线温度逐渐增加；喷管壳体表面温度随工作时间增加而增加，随离喷管出口的距离增加先增加而后逐渐降低；喷管喉部壳体表面温度最高，达 左右；推进剂燃温增加，喷管金属壳体表面温度呈线性增加；换热系数增加，金属壳体表面温度近似线性降低；燃烧室压强变化对金属壳体表面温度影响不大；隔热材料、和 降低金属壳体表面温度的效果突出。选用合适的隔热材料，做好喷管喉部位置的热防护是降低喷管壳体表面温度的有效方法。发动机试车结果与仿真结果吻合，验证了仿真分析方法的正确性。该瞬态流固热耦合一体化一步分析法，克服了过去流固热耦合方法步骤多、绝热假设不符合实际情况的局限性，可为合理设计固体发动机喷管、控制喷管温度提供依据。关键词：固体火箭发动机；喷管；长工作时间；流热固耦合；温度场分布中图分类号：文献标识码：文章编号：（）：，（，；，）：，（）（），收稿日期：；修回日期：。通讯作者：邓康清，男，博士 研究员，研究方向为宇航推进理论与工程。，：；引言固体火箭发动机喷管在发动机工作过程中承受高温、高压和高冲刷的恶劣环境。因此，耐烧蚀和热防护技术一直是发动机研究的重要内容。这就离不开发动机喷管的流场仿真和热力分析，国内外对此开展了大量研究。等利用二维轴对称可压缩 雷诺时均方程数值模拟发动机内流场温度分布，并通过试车试验验证。和 等基于有限元的方法对潜入式喷管的热响应进行了数值模拟分析。付鹏等采用有限元法研究了固体火箭发动机喷管喉衬结构的烧蚀量、瞬态温度场和应力场；张晓光等建立了基于 流体计算软件和 结构分析软件的流场热结构耦合分析二维模型，实现了流场与烧蚀传热的双向耦合及流场、热到结构的单向耦合；王寅虎利用 软件一维数值模拟了发动机火药启动器工作期间壳体瞬态温度场。韩珺礼等基于三维有限元法模拟分析了喷管工作时的温度场。孙林等研究了旋转固体火箭发动机喷管热结构，得到不同转速下的流场、温度场和应力场。白俊华等建立了喷管喉衬流固耦合换热模型，分析了喷管粗糙度、燃气成分、燃烧室压强和推进剂燃温等因素对喷管喉衬热结构及换热规律的影响。熊永亮等建立起轴对称的有限元计算模型，计算了喉衬组件与扩张段结构的瞬态温度场，研究了轴对称载荷下热结构耦合场的应力，得到喷管各层材料间的膨胀系数差异对应力的影响较大。陈国光等系统说明了固体火箭发动机喷管石墨喉衬的瞬态温度场和热应力场有限元计算方法，导出了横向同性材料轴对称瞬态温度场、热应力场的有限元公式，并对两种喉径的喉衬进行了数值计算，计算结果与发动机试车试验时 喉衬安全、喉衬发生裂纹结果是一致的。吴川等采用流固耦合方法，对金属及多种非金属材料组成的长尾喷管工作过程进行的非稳态传热数值计算表明，燃气对长尾喷管由内向外进行传热，但外壁面各区域温度场因内部材料不同差别较大。张小英等计算了美国 弹道导弹一级发动机的流场、壁面温度场和壁面烧蚀，发现喷管收敛段和喉部主要受到高温喷流的辐射加热，内壁辐射热流约为对流热流的 倍，喉部下游因喷流温度下降，速度激增，内壁对流热流超过辐射热流，在扩张段尾部，内壁的辐射热流再次超过对流热流；沿喷管壁面厚度自内向外，壁面温度急剧下降，发动机工作后 时，喉部截面处内壁温度达到 ，而外壁温度仅为 。目前，对固体火箭发动机喷管流场研究主要集中在短时间工作的喷管（小于 ）及喷管内表面温度场的模拟，而对长时间工作的固体火箭发动机喷管温度场的模拟研究不多，特别是对降低喷管外表面温度场的技术还少有研究。本文通过瞬态流固热一体化一步耦合分析法对一种采用无铝 推进剂的小推力长时间工作的固体火箭发动机喷管进行了传热分析，得到固体火箭发动机喷管燃气流场情况、喷管内外表面温度场分布及其随工作时间的变化规律，并用发动机试验验证了计算结果的正确性，从而为合理设计小推力长时间工作固体火箭发动机喷管，控制发动机温度提供指导和依据。计算模型 有限元模型由于本文研究的喷管具有对称性，取出一个子结构进行分析，从而极大简化了计算规模。取 进行分析，子结构三维几何结构见图。发动机喷管由喉衬、隔热层和壳体组成，外径 ，长度为 。瞬态流固热耦合计算区域及边界条件通常的流固热耦合方法是：假设流体与外界是绝热的，对发动机喷管中的流体进行静态流场仿真分析，得到喷管中燃气的热流场分布；再通过流体与喉衬界面耦合，得到喉衬内表面的热流密度情况；最后，通过瞬态热仿真得到各部件在各个时刻的温度场分布。这种方法步骤多且繁琐，并且绝热的假设与实际不符，将导致后续仿真结果与实际结果存在偏差。本文的瞬态流固热一体化一步耦合法将整个流体、喷管喉衬、隔热层和壳体作为一个整体考虑，只需一步瞬态流场仿真，即可得到流体流场分布和所有部件的温度场分布，年 月固体火箭技术第 卷而不再需要瞬态热仿真这一步。瞬态流固热耦合一体化一步仿真法是选取喷管入口为燃气进口，采用商用软件 对喷管流场进行计算。整个流场区域如图 所示。图 中，为喷管入口边界，给定燃气总压和总温，分别为 和 （）；为喷管壁面，采用无滑移壁面边界；是出口边界，采用压力场边界，压力为 。流体和固体间、固体和固体间采用耦合壁面，外壁面采用对流边界条件。整个流场和喷管主要采用六面体网格，流场壁面网格逐步加密。典型计算规模为 个节点，个单元。对该喷管同时进行瞬态流固热耦合，得到整个喷管温度场分布和燃气流的温度分布和压强分布，结果可作为降低固体发动机壁面温度的依据。NozzlethroatThermalinsulatedlayerCase（）Combustion gasInletNozzlethroatThermalinsulatedlayerCase（）图 发动机喷管局部模型 材料性能通过热力计算得到采用无金属 推进剂发动机喷管气相入口燃气种类和质量分数如表 所示。用于瞬态热分析的金属壳体、绝热层、耐烧蚀喷管材料的性能参数见表。模型验证 网格无关性验证为确定合理的计算网格分布，兼顾计算准确性和效率，对于发动机分别采用 （个节点）和（个节点）不同数量的网格进行仿真，对比不同网格尺寸的计算结果，差异性很小，最高温度仅相差 （密网格 ，疏网格 ）。综合考虑计算精度及资源，使用 网格尺度进行计算。表 气相入口燃气种类和质量分数 表 材料的性能参数 （）（）（）实验验证采用无金属 型固体推进剂，燃烧室理论燃温，制成相同喷管结构的固体发动机，该发动机长 ，直径 ，在 下进行了 试车试验，采用喷管表面埋置热电偶方法测试了发动机喷管喉径处壳体表面温度，并与表面温度仿真结果进行对比，结果如图 所示。10080604020-200 20 40 60 80 100 120 140Predicted=20 GFM2TsExperimentt/sT/图 发动机喷管段壳体表面温度实测结果与仿真结果 年 月邓康清，等：小推力长工作时间固体火箭发动机喷管温度仿真分析第 期 由图 可知，发动机喷管喉径处壳体表面温度时间仿真曲线与实验测试曲线趋势相同，但喷管喉径处壳体表面温度预示结果比实测结果高 左右，偏差约，喷管喉径处壳体表面预示温度结果与实测结果较吻合，验证了模型的可靠性。计算结果及分析 喷管流场和温度场典型计算结果利用 软件，对喷管中的气相组分，采用组分输运模型处理，燃气设为理想气体，在三维 方程的基础上，加入组分输运方程，采用标准 两方程模型对流场进行瞬态分析求解。其中，喷管材料为，推进剂燃温为 ，工作压强为 ，换热系数取 ，每个时间点迭代计算 次，直至计算到工作时间 结束。喷管内燃气温度场和喷管本身温度场分析结果如图 所示。从喷管中的压强分析结果最大压强 和温度场分析结果最高温度 可以看出，流场分析结果是可信的。p/Pat=130 s（）T/t=130 s（）图 燃烧 时发动机流场分析结果（）（）图 是燃烧 时发动机各部件传热分析结果。可知，喉衬温度最高，最低也有 ；隔热层温度次之，在 之间；壳体表面温度最低，在 之间，且从喷口到入口，壳体表面温度逐渐增加。（）（）（）图 燃烧 时发动机各部件传热分析结果（）（）25002000150010005000 20 40 60 80 100 120 140 160 1805 mm10 mm20 mm25 mm30 mm35 mm40 mm45 mmT/t/s（）22002000180016001400120010008000 10 20 30 40 50t=130 sT/L/mm（）图 工作时间和喷口距离对中心轴线上各点温度的影响 年 月固体火箭技术第 卷 图 是发动机工作过程中喷管中心轴线上燃气温度时间曲线和 时中心轴线上各点燃气温度距离曲线，图（）中的 指离喷管出口的距离。可知，随工作时间增加，喷管中心轴线上各点温度在点火后急剧增加，后基本保持不变；从喷管出口到喷管入口，喷管中心轴线上各点温度逐渐增加。图 为喷管中几个典型特征点的温度随工作时间的变化情况。个点的位置见图。点 在喷管喉部中心轴线上，点 在喷管喉衬直线段内表面上，点 在喷管喉衬直线段外表面上，点 在喷管喉部隔热层外表面上，点 在喷管喉部金属壳体外表面上。从中可知，随工作时间增加，喷管中心轴线上点的温度点火后急剧增加，后基本保持不变；喷管喉部中心轴线上温度最高，喉衬内表面温度次之，喉衬外表面温度比内表面温度低，但两者很接近；由于隔热层的隔热作用，隔热层外表面温度大大降低，金属壳体外表面温度与隔热层外表面温度相近。20001500100050000 20 40 60 80 100 120 140T/t/spoint 1point 2point 3point 4point 5图 喷管中几个典型特征点的温度变化情况 发动机喷管表面温度变化规律 发动机喷管表面位置的影响图 为发动机喷管工作 过程中壳体表面温度工作时间曲线和壳体表面最高温度距离曲线，图（）中的 指离喷管出口的距离。可知，喷管壳体表面温度随工作时间增加而增加，随离喷管出口的距离增加先增加，喷管喉部壳体表面温度最高，而后逐步降低。经过 工作后，这种发动机喷管表面温度最高达到 左右，这对舱内的设备不利，因此需要降低喷管表面温度。推进剂燃温和燃烧室压强的影响推进剂燃温和燃烧室压强对金属壳体表面温度的影响见图 和图。推进剂燃温 变化是在压强 不变下，由改变推进剂配方组分含量达到的，再由全流场计算得到 金属壳体表面最高温度。同时，由于表 中热力计算得到，在燃烧室压强 范 围 内，推 进 剂 燃 温 仅 变 化 （），推进剂燃气成分最大变化不到，可以认为，燃烧室压强对推进剂燃温和燃气成分影响可以忽略，因此在流场仿真时可以直接改变燃烧室压强进行计算。从图 和图 可知，推进剂燃温增加，喷管金属壳体表面温度呈线性增加；燃烧室压强变化对金属壳体表面温度影响不大。2502001501005000 20 40 60 80 100 120 140 160 1805 mm10 mm15 mm20 mm25 mm30 mm35 mm40 mm45 mmT/t/s（）2402202001801601401201000 10 20 30 40 50t=130 sT/L/mm（）图 离喷管出口距离对喷管金属壳体表面温度的影响 2802602402202001801601400 1600 1800 2000 2200 2400 2600Tc/T/t=130 s=2.2+0.1098=0.9999TTRc图 推进剂燃温对金属壳体表面温度的影响 换热系数对喉径处金属壳体表面温度的影响换热系数对喉径处金属壳体表面温度的影响见图。可知，喷管壳体表面温度随工作时间增加而增 年 月邓康清，等：小推力长工作时间固体火箭发动机喷管温度仿真分析第 期加；换热系数增加，金属壳体表面温度几乎呈线性降低。这是因为换热系数越大，通过换热导出的热量越多，金属壳体表面温度降低越大。2502001501005000 20 40 60 80 100 120 140pc=3.4 MPapc=4.4 MPa=5.4 MPapc=6.4 MPapc=7.4 MPapcT/t/s图 燃烧室压强对金属壳体表面温度的影响 2502001501005000 20 40 60 80 100 120 140t/sT/h=25 Wm-2-1K-2-1h=20 Wm K-2-1h=100 Wm K-2-1h=150 Wm K-2-1h=200 Wm K-2-1h=250 Wm K（）2602402202001801601401200 50 100 150 200 250T/t=130 s-2-1h/(Wm K)（）图 换热系数对喉径处外金属壳体表面温度的影响 喉衬材料的影响喉衬材料和隔热材料对金属壳体表面温度的影响见图 和图。可知，、和 的壳体表面温度相同，只有采用 有一定的降低金属壳体表面温度的作用；隔热材料具有显著的降低金属壳体表面温度的作用，其作用大小为：该顺序与上述几种材料的导热系数大小顺序正好相反。说明隔热材料导热系数是影响壳体表面温度的主要因素：隔热材料导热越慢，隔热材料阻止热量往外导出到壳体表面的作用越大，隔热材料降低壳体表面温度的作用越大。2502001501005000 20 40 60 80 100 120 140t/sCSWCu10WCCT/图 喉衬材料对金属壳体表面温度的影响 40030020010000 20 40 60 80 100 120 140t/sTFMGFM5-IIIGFMBEPDMGFM2T/图 隔热材料对金属壳体表面温度的影响 结论（）得到了喷管的温度变化规律，喷管中心轴线上各点的温度开始时急剧增加，后基本保持不变；从喷管出口到喷管入口，喷管中心轴线上各点的温度逐渐增加。喷管喉部中心温度最高（以上），喉衬内表面温度次之，喉衬外表面温度比内表面温度低，但两者很接近（）；由于隔热作用，隔热层外表面温度大大降低（）；金属壳体外表面温度最低（），和隔热层外表面温度接近。说明隔热层是喷管各部组件中降低壳体外表面温度最有效部件。（）揭示了喷管壳体表面温度的变化规律，喷管壳体表面温度随工作时间增加而增加，随离喷管出口的距离增加先增加而后逐步降低。喷管喉部壳体表面温度最高，时最高达 左右。说明喷管喉部位置的绝热防护是降低壳体表面温度最重要部位。（）推进剂燃温增加，喷管金属壳体表面温度呈线性增加；换热系数增加，金属壳体表面温度近似线性降低；燃烧室压强变化对金属壳体表面温度影响不大。（）隔热材料具有显著降低金属壳体表面温度的 年 月固体火箭技术第 卷作用，其作用大小为 。（）本文采用的瞬态流固热一体化一步耦合分析法，克服了以往先假设绝热进行流场分析，再将结果耦合进行热分析方法的步骤多、绝热假设不符合实际情况的局限性，具有开拓性，对类似研究工作具有指导性。可为合理设计固体火箭发动机喷管，控制发动机喷管表面温度提供设计依据。本方法预估的喷管表面温度与实测结果仍有约的偏差，后续可进一步开展提高喷管表面温度预估精度的方法研究。参考文献：，（）：付鹏，蹇泽群，张钢锤，等发动机喷管喉衬烧蚀及热结构工程计算固体火箭技术，（）：，（）：张晓光，王长辉，刘宇，等固体火箭发动机喉衬流场及热结构耦合分析固体火箭技术，（）：，（）：王寅虎火药启动器壳体瞬态温度场的数值模拟分析导弹与航天运载技术，（）：，（）：韩珺礼，蒙上阳，杨军辉，等固体火箭发动机喷管热防护设计计算机辅助工程，（），（），孙林，鲍福廷，时伟旋转固体火箭发动机喷管热结构分析航空计算技术，（）：，（）：白俊华，胡春波，李佳明固体发动机喷管喉衬流固耦合换热规律数值研究西北工业大学学报，（）：，（）：熊永亮，郜冶喷管温度与应力场的数值研究哈尔滨工程大学学报，（）：，（）：陈国光，孙南生，薛献新固体火箭发动机喷管石墨喉衬的温度场和应力场研究推进技术，（）：，（）：吴川，邢国强，门们固体火箭发动机长尾喷管传热数值模拟弹箭与制导学报，（）：，（）：张小英，向红军固体火箭发动机喷管传热与壁面烧蚀的耦合计算分析固体火箭技术，（）：，（）：郑亚，陈军，鞠玉涛，等固体火箭发动机传热学北京：北京航空航天大学出版社，王元有，赵伯华，余世方，等固体火箭发动机设计北京：国防工业出版社，（编辑：吕耀辉）年 月邓康清，等：小推力长工作时间固体火箭发动机喷管温度仿真分析第 期