基金项目:863计划项目(2015AA7023021)
作者简介:巩岩博(1994—),男,硕士,研究领域为液体火箭发动机系统设计
(Beijing Aerospace Propulsion Institute,Beijing 100076,China)
cryogenic rocket engine; steady-state characteristic simulation; performance reliability; performance sensitivity
为了提高发动机工程研制的效率,以低温火箭发动机为研究对象,使用通用性好、支持模块化建模与仿真的Modelica语言,建立了具有较高仿真精度的低温火箭发动机稳态特性仿真模型库与仿真分析平台,并针对工程应用的需要开展了氢氧发动机故障仿真、液氧/甲烷发动机性能可靠性评估和性能敏感性分析等工作,为某型氢氧发动机的故障定位提供了依据,同时得到了某型液氧/甲烷发动机的性能可靠性指标,以及对其性能影响最大的干扰因素,为该型发动机的研制工作提供了保障。
In order to improve the efficiency of engine engineering development, the cryogenic rocket engine was used as the research object, a simulation model library and simulation analysis platform with high simulation accuracy for the steady-state characteristics of the cryogenic rocket engine were established by using Modelica language with good generality, supporting the modular modeling and simulation.In addition, aiming at the requirements of engineering application, the fault simulation of LOX/LH2 engine, the performance evaluation and performance reliability sensitivity analysis of LOX/LCH4 engine were carried out.It provides a basis for the fault location of a certain LOX/LH2 engine, and obtains the performance reliability index of a certain LOX/LCH4 engine, as well as the interference factors that have the greatest impact on its performance, which provides a guarantee for the development of this type of engine.
以液氧、液氢、液甲烷等低温介质作推进剂的液体火箭发动机被称作低温火箭发动机; 低温火箭发动机具有低温、高能、无毒、环保等特点[1],是各航天大国的重点发展方向。稳态特性仿真分析是低温火箭发动机研发过程中必不可少的一个重要环节,从发动机方案优化设计、参数选择与预估到系统级试车方案的制定、偏差分析与故障分析,都需要进行发动机稳态特性仿真分析,它与试验研究相辅相成,是实现发动机设计目标、降低发动机试验成本、提高工程研制效率的重要手段。
大推力低温火箭发动机结构复杂,由大量组件构成,是一个非线性很强的复杂热力学系统。对于大推力低温火箭发动机来说,发动机包含的零组件众多,系统复杂,描述发动机状态的非线性方程众多且具有很强的刚性,如算法不当将出现振荡而无法收敛,数值求解难度较大。相比于小推力低温火箭发动机,大推力低温火箭发动机泵后温度和密度变化很大,特别是液氢具有一定的可压缩性且泵后处于超临界状态,若不考虑推进剂在涡轮泵及管路中的温度和密度变化情况,采用传统模型计算出来的推进剂流量将产生很大偏差[2],严重影响仿真精度。
本文开展了低温火箭发动机稳态特性仿真分析研究,建立了组件方程合理、实用化、通用化良好,具有较高工程精度的可视化仿真分析平台,并成功应用于发动机故障诊断与分析、性能可靠性评估、性能参数敏感性分析等工作,使用效果良好,提升发动机系统稳态特性认知水平,指导发动机工程研制。
发动机在稳态工况下,各组件的流量、压力功率处于平衡状态,因此可根据发动机系统组成,按照发动机各流路的压力平衡、流量平衡和功率平衡,结合发动机各部组件的特性方程,建立发动机系统的非线性方程组[3]
Fi(D,X)=0 i=1,2,…,n(1)
式中:D表示发动机各部件的特性参数; X表示发动机性能参数。本文针对低温推进剂温度与密度变化较大的特点,引入泵后温升模型; 针对现有的夹套传热计算模型在计算氢氧火箭发动机时精度不高的问题,通过理论分析与试验数据归纳提出了相应的改进方案,针对基于不可压缩流伯努利方程得到的汽蚀文氏管模型不适用于液氢的问题,引入了等熵算法。
推进剂的物性会对泵的特性产生较大的影响。对于低温火箭发动机,推进剂的温度和密度在流经泵的过程中会发生较大的改变,如果忽略这种变化将会造成较大的计算误差,因此为了准确计算离心泵的特性,必须对泵后密度和温度进行计算。本文引入基于泵后压力和比焓计算推进剂温升和密度变化的方法进行计算[4]。
已知泵前推进剂的比焓h1时,泵后推进剂的比焓
h2≈h1+(1-ζ+ζ/(ηp))(Δp)/(p-)(2)
式中ζ为泵功率总损失中转化为热量的比例。根据试验数据统计,一般取0.88~0.95。
泵前比焓可以根据泵前推进剂的温度与压力得到,即
h1=fh(p1,T1)(3)
根据热力学原理,有
T2=fT(p2,ρ2)(4)
ρ2=fρ(h2,p2)(5)
式中:下标1为泵入口参数; 下标2为泵出口参数。式(4)和式(5)的计算需要介质热力参数之间的函数关系,可以查热力参数表获得。
将式(2)~(5)与泵特性方程联立,即可迭代求解出泵出口的介质温度与介质密度。
使用改进的涡轮泵模型对某型发动机的主泵进行计算,并与试验数据进行对比,得到表1和表2。
由对比可知,考虑泵后推进剂温升与密度变化的涡轮泵模型能够较好地模拟低温涡轮泵的工作特性,具有较高的计算精度。
本文采用巴兹公式[5]对冷却夹套开展传热计算,并针对氢氧推力室传热试验数据对其进行了修正,引入修正系数
Cre=7.55[lg((pc)/(pe))-1]·((dt/dl)0.1)/(rc)(6)
式中:pc为推力室室压; pe为大气压强; rc为推力室混合比; dt为推力室喉部直径; dl为圆筒段直径。所有参数均采用国际制基本单位。则修正后的燃气侧传热系数
h'g=Crehg(7)
式中hg为通过巴兹公式计算得到的传热系数。
使用修正后的传热计算模型分别对某型50 t级氢氧火箭发动机的低、中、高工况状态进行计算,得到夹套温升与压降,并与试车试验数据对比,发现改进后的模型具有更高的计算精度。如表3~表5所示。
表3 低工况状态(pc=9.937 MPa, r=6.06)
Tab.3 Low operating condition(pc=9.937 MPa, r=6.06)
表4 额定工况状态(pc=10.195 MPa, r=6.517)
Tab.4 Rated working condition(pc=10.195 MPa, r=6.517)
表5 高工况状态(pc=10.391 MPa, r=6.791)
Tab.5 High operating condition(pc=10.391 MPa, r=6.791)
常规的汽蚀文氏管计算公式由不可压缩流伯努利方程[6]导出,并近似认为喉部温度等于入口温度,只适用于不可压缩或压缩性较弱的介质。液氢与甲烷都具有较强的压缩性,不可压缩流伯努利方程不适用。特别是液氢,随着温度、压力的变化,其密度会发生显著改变。在汽蚀管喉部截面,由于流速的增大,液氢的温度会明显低于入口温度,因此必须根据汽蚀管喉部的实际参数计算流量。
一般认为液氢在汽蚀管中流动是等熵膨胀过程,因此有
S1=St(8)
式中:S1为汽蚀文氏管入口比熵; St为汽蚀文氏管喉部比熵。其中入口比熵可以根据入口压力和温度查表得到,即
S1=fs(p1,T1)(9)
喉部处于饱和态,为了方便计算,近似认为喉部处于汽蚀临界点,即此处的气体体积分数为0,有
St=fS(Tt,Q=0)(10)
联立式(8)~(10)即可迭代求解出喉部温度βt。
根据入口温度和入口压力查表可得入口比焓
h1=fh(p1,T1)(11)
在喉部处于饱和态的前提下,根据喉部温度和喉部气体体积分数为零,查表可得喉部比焓和密度
ht=fh(Tt,Q=0)(12)
ρt=fρ(Tt,Q=0)(13)
根据热力学原理[7],假设汽蚀管入口速度为0时,喉部处的流速
式中:Vt为喉部流速,m/s; h1为入口处推进剂比焓,J/kg; ht为喉部处推进剂比焓,J/kg。最后,根据连续方程计算喉部流量
qm=ρtVtAt(15)
式中qm为喉部质量流率,kg/s。
表6所示为使用等熵算法计算得到的汽蚀文氏管流量与试验结果的对比,可以看出,等熵算法具有较高的计算精度。
Modelica是一种面向对象的非因果建模语言,可以实现多领域物理系统建模和工程仿真,支持模块化建模与仿真,在使用Modelica进行建模时,不需要对数学模型进行额外的推导与变形,可以直接基于原始的非线性方程进行建模[8]。本文针对低温火箭发动机,基于MWorks仿真平台,在Modelica标准库的基础上进行了二次开发,得到低温火箭发动机稳态特性模块化仿真模型库,如图1与图2所示。
该模型库包含低温火箭发动机的主要部件模块,支持所有基于Modelica语言的建模与仿真平台,能够进行可视化与模块化操作,具备对各类低温火箭发动机的分系统和全系统开展模块化建模与仿真的能力,可以大大提高建模与仿真的效率。
某型氢氧发动机为燃气发生器循环、并联涡轮泵布局,推力室身部使用液氢进行再生冷却,喷管延伸段以液氢为介质进行排放冷却[9]。发动机系统如图3所示。
在发动机稳定工作阶段,当故障出现并发展时,发动机将出现不同程度的性能下降、参数不协调及元组件损坏等问题。利用发动机稳态特性仿真分析平台,可以在不分解发动机的情况下,通过稳态特性仿真的方法对故障模式开展初步的定位与评估。如某次发动机试车过程中主要性能参数出现了约3%的波动如图4所示。根据试验数据,针对可能的故障模式建立发动机故障树(见图5),利用发动机稳态特性分析平台对试验数据进行了仿真分析,找到仿真结果与试验参数相吻合的故障模式,在未对产品进行分解检查的情况下,实现对故障进行了定性和定量分析,帮助故障快速成功定位,结果与产品分解检查情况一致,与试验参数相吻合的仿真计算结果如图6所示,直方图表示氧泵扬程与效率发生变化时发动机各主要参数的偏差情况,可以看出,仿真结果与试验参数的实际变化非常接近。
某型液氧/甲烷发动机采用燃气发生器循环方式,推进剂供应系统采用单台富甲烷燃气发生器、同轴涡轮泵方案,再生冷却推力室。发动机系统如图7所示。
图7 液氧/甲烷发动机系统示意图
Fig.7 System schematic of LOX/LCH4 engine
在发动机研制过程中的某个阶段,虽然发动机整机级试验样本较少,但却拥有相对数量可观的子系统级与部件级试验数据样本,这些数据包含了丰富的部件特性信息,由这些部件组成的发动机整机性能可靠性与部件特性可靠性密切相关,通过建立合理的映射关系与数学模型,就能够在没有或者仅有很少整机性能试车样本的情况下获得发动机的整机性能可靠性[10-11]。
基于发动机组件及子系统的试验数据,通过随机方法[12-13]可以生成大量满足实际分布规律的部组件特性参数,运用低温火箭发动机稳态特性仿真分析模型库建立发动机部件特性与整机性能参数之间的映射关系,就可以通过仿真获得大样本条件下的发动机整机性能参数,然后通过统计学方法,就可以得到发动机的性能可靠性指标。本文通过上述步骤,经过统计分析,获得了某型液氧/甲烷火箭发动机整机性能可靠性指标,其中,各部件的特性参数X均服从正态分布,以各参数的设计值为基准,对参数特性进行无量纲化,通过统计得到各参数的数学期望μ和标准差σ如表7所示。通过大量部组件特性参数开展随机仿真的结果如图8和图9 所示,发动机性能参数落在规定范围内的可靠性点估计结果如表8所示,这为在研发动机的性能可靠性评估提供了一种可行的方法。
火箭发动机的主要性能参数既受到自身部组件特性的影响,同时也受到外部环境因素的影响。通常把性能参数受自身部组件特性和外部因素的影响而产生的变化量称为性能参数对影响因素的敏感度[14]。利用发动机稳态特性仿真分析平台对影响图7所示液氧/甲烷发动机性能的因素开展了仿真与分析,得到了对发动机性能的影响较大的因素以及这些因素的影响程度,有利于设计人员聚焦于这些关键的影响因素,帮助设计人员加强敏感因素的质量控制与性能分析[15]。
能够对发动机性能产生影响的因素可分为两类,即内部影响因素和外部影响因素。外部影响因素是指由外部环境或发动机推进剂带来的干扰,如泵入口压力的变化、推进剂温度的变化和喷管出口的背压等,这种干扰一般可以准确测量并采取相应的措施进行控制; 内部因素则是指发动机自身带来的干扰,如在制造过程中的加工误差、零部件的安装误差或部件性能试验的测量系统误差等,这种误差具有随机性,难以准确测量与控制[10]。因此在进行性能敏感性分析时,主要关注内部干扰因素。
从仿真准确性和计算规模两个方面综合考虑,在单因素敏感性分析中,每个因素的水平数均取5,即基准值和分别偏离基准值-20%、-10%、10%、20%; 在多因素敏感性分析中,每个因素的水平数均取3,即基准值和分别偏离基准值-10%、10%。从计算结果可以看出,在本文所选定的参数变化范围之内,各内外因素对发动机的推力和混合比的影响有着较为相似的敏感度,可以得到基本相同的变化趋势,即随着影响因素变化率的增大,发动机推力和混合比的变化增大,二者近似呈线性关系。对于发动机推力,其敏感度大小为βt>ηpf>ηpo>Cqgov>Cqgfv>Cqgo>μt>Cqgf,其中涡轮、甲烷泵及氧泵的效率对发动机推力的影响程度最大; 对于发动机混合比,其敏感度大小为βt>ηpf>ηpo>Cqmov>Cqco>Cqcf>Cqgfv>μt,其中对发动机混合比影响程度最大的同样是涡轮效率、甲烷泵效率和氧泵效率。
推力和混合比对不同影响因素的平均敏感度直方图如图 10所示。
图 10 单因素下发动机主要性能敏感度分析结果
Fig.10 Sensitivity analysis results of engine main performance under single factor
由分析结果可知,氧泵和甲烷泵入口压力,以及推进剂供应管路的流阻系数对发动机性能的影响相对较小,为提高分析效率,在开展多因素敏感性分析时可以忽略这些因素的影响,仅选取影响较大的因素开展分析。根据仿真与分析结果,在多因素共同作用时,发动机性能参数对各内外影响因素所表现出的敏感性与单因素作用时的结果基本一致。其中,影响发动机推力的因素重要性排序依次为:βt>ηpf>ηpo>Cqgov>Cqgfv>Cqgo>μt>Cqgf,影响发动机混合比的因素重要性依次为: βt>ηpf>ηpo>Cqmov>Cqco>Cqcf>Cqgfv>μt>Cqgf>Cqgov。根据极差分析的结果绘制发动机推力和混合比的极差直方图,分别如图 11和图 12所示。
可以看出,不同影响因素在多因素共同作用下对发动机性能的影响程度有所不同,如发动机推力对发生器氧阀的流阻系数较为敏感,而发动机混合比则对推力室氧喷注器的流阻系数较为敏感。通过分析得到对发动机性能影响较大的因素后,就可以采取措施,如加强零部件生产过程中的质量控制,或针对相应的部组件开展大量的地面试验,掌握部组件特性,降低因部组件特性的不确定性而对发动机整体性能造成的影响。此外,通过敏感性分析还可以给出发动机主要调节元件与发动机性能之间的定量关系,对提高调整精度和分析效率有一定帮助。以低温火箭发动机为研究对象,使用通用性良好、具备模块化建模与仿真能力的Modelica语言[16],建立组件方程合理、通用化、实用化良好,能够更为准确地反映组件实际工作特性的低温火箭发动机可视化稳态特性仿真分析平台,具有界面友好、可模块化建模与仿真、通用化良好的特点[17-19]。
在此基础上,通过发动机稳态特性仿真分析的方法[20],开展了氢氧火箭发动机故障诊断分析、液氧/甲烷发动机的性能可靠性评估与性能参数敏感性分析,效果良好,促进了发动机工程研制,提升了发动机稳态特性认知水平。在该平台基础上进行适应性改进与调整,补充相应的物性参数后,还可以用于其它低温火箭发动机稳态特性仿真计算与分析,具有良好的应用前景。