-
液体火箭发动机故障仿真,是通过在描述发动机动态过程的常微分方程组中,嵌入描述发动机故障的特性方程或改变发动机边界条件,然后进行数值求解来实现的。在建立了描述发动机动力学系统的精确数学模型后,故障仿真可以对在故障状态下的发动机工作参数的变化过程进行准确和细致的描述,从而可以比较全面地了解在故障状态下的发动机系统动力学行为[1]。
20世纪70年代,美国针对航天飞机主发动机开展了故障诊断与健康监控的相关研究,随后,美国和欧洲针对液体火箭发动机故障过程的动态特性开展了大量的仿真研究[2]。国内则以国防科技大学为代表[3],主要对常规推进剂发动机、液氧煤油发动机[4]和空间推进系统[5]开展了故障过程的动态特性建模与仿真研究,但对氢氧发动机故障过程的研究则较为缺乏。
为了实现新一代运载火箭芯一级氢氧发动机的故障检测与定位,基于AMESim仿真软件平台,自主编写开发了液体火箭发动机动态仿真模型库,在合理简化系统的基础上,构建了50 t级大推力氢氧发动机动态仿真模型。针对该型发动机某次试验故障,本文开展了多种故障模式的注入与动态仿真,对发动机故障点附近的动态特性进行了分析,与试验实测参数曲线进行详细对比,提出了可能的故障模式,为大推力氢氧发动机故障的准确定位提供了参考。
-
近似认为推进剂在管路内温度恒定,其密度、黏性等物性参数均只与压力相关,由此可将一维可压缩流分段有限元模型进行简化[6]。
设管道分段数为n,则对于状态变量控制单元,采用容性方程形式
$$\begin{aligned}[b] & {\dfrac{{{\rm{d}}{p_i}}}{{{\rm{d}}t}} = \dfrac{1}{{{C_i}}}\left( {{Q_i} - {Q_{i + 1}}} \right) + \varepsilon \dfrac{{{a_i}}}{{{A_i}}} \cdot \rho \left( {\dfrac{{{\rm{d}}{Q_i}}}{{{\rm{d}}t}} - \dfrac{{{\rm{d}}{Q_{i + 1}}}}{{{\rm{d}}t}}} \right),}\\ &{i = 0,1,\cdots,n - 1} \end{aligned}$$ (1) 对于速度变量控制单元
$$\begin{aligned}[b] & \begin{array}{*{20}{l}} {\dfrac{{{\rm{d}}{Q_i}}}{{{\rm{d}}t}} = \dfrac{1}{{{L_i}}}\left( {\dfrac{{{p_{i - 1}} - {p_i}}}{\rho } - f\rho \left| {{Q_i}} \right|{Q_i}} \right),}&{i = 0,1,\cdots,n} \end{array} \\ &{p_{ - 1}} = {p_a},{p_n} = {p_b} \end{aligned}$$ (2) 其中:
${p_i}$ 为第i个状态变量单元的压力;${C_i} = {{{V_i}} / {\left( {\rho {a^2}} \right)}}$ 为第i个状态变量单元的流容;${V_i}$ 为第i个状态变量单元体积;${A_i}$ 为第$i$ 个速度变量单元内的体积流量;${a_i} = \sqrt {\dfrac{{{B_i}}}{{{\rho _i}}}} $ 为流体声速;${Q_i}$ 为第i个速度变量单元内的体积流量;$\varepsilon \dfrac{{{a_i}}}{{{A_i}}}$ 是为避免单元压力与流量耦合产生高频振荡而增加的阻尼项,$\dfrac{{{a_i}}}{{{A_i}}}$ 是容积波动阻力,$\varepsilon $ 为可根据单元数调整的系数,调整范围为0.1~0.5;${L_i} = {{{l_i}} / {{A_i}}}$ 为第i个速度单元的流感;$f$ 为准稳态摩擦系数。 -
如图1所示,将泵入口端面和诱导轮入口端面之间的体积视为集中体积,该体积应包括入口容积以及诱导轮、叶轮、扩散器的过流容积。
1)泵入口处压力方程
泵入口单元压力方程可由集中体积的流容方程计算[6]
$$\dfrac{{{\rm{d}}{p_1}}}{{{\rm{d}}t}} = \dfrac{1}{{C_{{\rm{pump}}}^{{\rm{in}}}}}\left( {{Q_1} - {{Q'}_1}} \right)$$ (3) 其中:
$C_{{\rm{pump}}}^{{\rm{in}}} =( {{V_{{\rm{pump}}}^{{\rm{in}}}}) / {{\rho _1}}}{a^2}$ 为泵入口集中体积等效流容;${\rho _1}$ 、$a$ 分别为泵入口流体密度、声速。2)离心泵效率、稳态扬程和转矩的计算
离心泵效率与扬程、转矩的关系为
$${\eta _p} = \dfrac{{\Delta pQ}}{{{\boldsymbol{\omega}} T}} = \dfrac{{30\Delta pQ}}{{\text{π} nT}}$$ (4) 其中:
${\eta _p}$ 为离心泵效率;${\boldsymbol{\omega}} $ 为泵叶轮角速度;$\Delta p$ 为扬程;$Q$ 为泵体积流量;$T$ 为转矩。泵的稳态扬程和效率计算公式由实验数据拟合得到,具体如下
$$\Delta p = \alpha \bar \rho \left( {{c_1}{Q^2} + {c_2}Q \cdot n + {c_3}{n^2}} \right)$$ (5) $${\eta _p} = \beta \left[ {{d_1}{{\left( {\dfrac{Q}{n}} \right)}^2} + {d_2}\left( {\dfrac{Q}{n}} \right) + {d_3}} \right]$$ (6) 其中:
$\alpha $ 为扬程修正系数;$\bar \rho = {{\left( {{\rho _1} + {\rho _2}} \right)} / 2}$ 为泵内推进剂平均密度;$\beta $ 为效率修正系数。3)泵出口参数计算
离心泵出口压力为
$${p_2} = {p_1} + \Delta p$$ (7) -
将涡轮考虑为一个C-R-C模型,即等效为一个集中体积–阻力–集中体积的组合元件[6-7]。
1)涡轮工质流动的动力学方程
如图2所示,为涡轮等效元件的示意图。对于两端的集中容积,根据质量守恒和能量守恒定律,不难得到
$$\dfrac{{{\rm{d}}{p_1}}}{{{\rm{d}}t}} = \dfrac{\gamma }{{{V_1}}}\left( {p_1^{{\rm{in}}}{Q_1} - {p_1}\dfrac{{{{\dot m}_t}}}{{{\rho _1}}}} \right)$$ (8) $$\dfrac{{{\rm{d}}{\rho _1}}}{{{\rm{d}}t}} = \dfrac{1}{{{V_1}}}\left( {\rho _1^{{\rm{in}}}{Q_1} - {{\dot m}_t}} \right)$$ (9) $$\dfrac{{{\rm{d}}{p_2}}}{{{\rm{d}}t}} = \dfrac{\gamma }{{{V_2}}}\left( {\left( {R{{T'}_2}} \right){{\dot m}_t} - p_2^{{\rm{exit}}}{Q_2}} \right)$$ (10) $$\dfrac{{{\rm{d}}{\rho _2}}}{{{\rm{d}}t}} = \dfrac{1}{{{V_2}}}\left( {{{\dot m}_t} - \rho _2^{{\rm{exit}}}{Q_2}} \right)$$ (11) 其中:
${T'_2}$ 为涡轮动叶出口处气体温度;$R$ 为燃气气体常数;${\dot m_t}$ 为通过涡轮的燃气质量流量,按下式确定$${\dot m_t} = \dfrac{{{C_t}{A_t}p_1^ * }}{{\sqrt {RT_1^ * } }}{\Gamma _t}$$ (12) $${\Gamma _t} = \left\{ {\begin{array}{*{20}{c}} {\sqrt {\gamma {{\left( {\dfrac{2}{{\gamma + 1}}} \right)}^{\frac{{\gamma + 1}}{{\gamma - 1}}}}} ,\dfrac{{{p_{{\rm{bi}}}}}}{{{p_1}}} \leqslant {{\left( {\dfrac{2}{{\gamma + 1}}} \right)}^{\frac{\gamma }{{\gamma - 1}}}}} \\ {\sqrt {\dfrac{{2\gamma }}{{\gamma - 1}}\left[ {{{\left( {\dfrac{{{p_{{\rm{bi}}}}}}{{{p_1}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\dfrac{{{p_{{\rm{bi}}}}}}{{{p_1}}}} \right)}^{\frac{{\gamma + 1}}{\gamma }}}} \right]} ,\dfrac{{{p_{{\rm{bi}}}}}}{{{p_1}}} > {{\left( {\dfrac{2}{{\gamma + 1}}} \right)}^{\frac{\gamma }{{\gamma - 1}}}}} \end{array}} \right.$$ (13) 其中:上标*表示滞止参数;
${C_t}$ 为涡轮喷嘴流量系数;${A_t}$ 为涡轮喷嘴的最小流通截面积;$\gamma $ 为燃气绝热指数;${p_{{\rm{bi}}}}$ 为静叶喷嘴出口压力,按下式确定$${p_{{\rm{bi}}}} = {p_1}{\left[ {{\Omega _T} + \left( {1 - {\Omega _T}} \right){{\left( {\dfrac{{{p_2}}}{{{p_1}}}} \right)}^{\frac{{\gamma - 1}}{\gamma }}}} \right]^{\frac{\gamma }{{\gamma - 1}}}}$$ (14) 其中,
${\Omega _T}$ 为涡轮反力度。2)涡轮功率、效率和转矩方程
涡轮轴功率为
$${N_t} = {\eta _t}{L_{{\rm{ad}}}}{\dot m_t}$$ (15) 其中:
${\eta _t}$ 为涡轮效率;${L_{{\rm{ad}}}}$ 为涡轮单位工质等熵膨胀功,由下式确定$${L_{\rm ad}} = {c_p}\left( {T_1^ * - T_2^ * } \right) = \dfrac{\gamma }{{\gamma - 1}}RT_1^ * \left[ {1 - {{\left( {\dfrac{{p_2^ * }}{{p_1^ * }}} \right)}^{\frac{{\gamma - 1}}{\gamma }}}} \right]$$ (16) 涡轮效率采用实验数据拟合方程。氢、氧涡轮拟合方程分别如下
$${\eta _{t,{H_2}}} = {\alpha _{t,{H_2}}}\left[ {{\xi _0} + {\xi _1}\left( {\dfrac{u}{C}} \right) + {\xi _2}{{\left( {\dfrac{u}{C}} \right)}^2} + {\xi _3}{{\left( {\dfrac{u}{C}} \right)}^3}} \right]$$ (17) $${\eta _{t,{O_2}}} = {\alpha _{t,{O_2}}}\left[ {{\varphi _0} + {\varphi _1}\left( {\dfrac{u}{C}} \right) + {\varphi _2}{{\left( {\dfrac{u}{C}} \right)}^2}} \right]$$ (18) 其中:
${\alpha _{t,{H_2}}}$ 、${\alpha _{t,{O_2}}}$ 分别为氢、氧涡轮效率修正系数;${\xi _i}$ 、${\varphi _i}$ 均为拟合系数;${u / C}$ 为涡轮速比,$u = {{\text{π} {D_t}n} / {60}}$ ,${D_t}$ 为动叶直径,$n$ 为转速;$C = \sqrt {2{L_{{\rm{ad}}}}} $ 。涡轮转矩为
$${T_t} = \dfrac{{{N_t}}}{{\boldsymbol{\omega}} } = \dfrac{{30{N_t}}}{{\text{π} n}}$$ (19) 其中,
${\boldsymbol{\omega}} $ 为动叶角速度。当发动机发生某些故障,如涡轮压比剧烈变化时,涡轮速比增大,但涡轮效率是减小的,如仍采用单一的以
${u / C}$ 描述的公式计算效率,将会得到效率升高的错误结论。为此,必须采用如下速比和压比二维拟合公式的形式计算发动机大范围变工况时的涡轮效率。$$ \begin{aligned} & \quad\quad\quad\quad\quad\quad\quad\quad {\eta _t} = f\left( {{\text{π} _t},\dfrac{u}{C}} \right) = \\ & {\alpha _0} \cdot \left( {{\eta _0} + a \cdot {\text{π} _t} + b \cdot \dfrac{u}{C} + c \cdot {\text{π} _t}^2 + d \cdot {{\left( {\dfrac{u}{C}} \right)}^2} + f \cdot {\text{π} _t} \cdot \dfrac{u}{C}} \right) \end{aligned}$$ (20) -
根据转动的动力学基本定律,涡轮泵转子的动量矩与作用在其上的力矩有如下关系
$$\dfrac{{{\rm{d}}{\boldsymbol{M}}}}{{{\rm{d}}t}} = \sum\limits_i {{{\boldsymbol{T}}_i}} = {{\boldsymbol{T}}_t} - \sum\limits_i {{{\boldsymbol{T}}_{p,i}}} $$ (21) 其中:
${{\boldsymbol{T}}_t}$ 、$\sum\limits_i {{{\boldsymbol{T}}_{p,i}}} $ 分别是作用在转子上的涡轮、泵的合力矩;M为转子动量矩,可由下式得到$${\boldsymbol{M}} = J{\boldsymbol{\omega}} $$ (22) 其中:J为转子转动惯量;
${\boldsymbol{\omega}} $ 为转子角速度。由上述两式及${\boldsymbol{\omega}} = {{2\text{π} n}/ {60}}$ 可得到涡轮泵转子的动力学方程$$\dfrac{{{\rm{d}}n}}{{{\rm{d}}t}} = \dfrac{{30}}{{\text{π} J}}\left( {{T_t} - \sum\limits_i {{T_{p,i}}} } \right)$$ (23) -
发动机系统中燃气发生器、推力室均为热力组件。对于燃气发生器及推力室的燃烧区,建立其数学模型时作如下假设:①忽略燃烧时滞;②燃烧过程是绝热的;③燃烧生成燃气的温度、组分浓度等在燃烧区内瞬时均匀一致;④燃气为理想气体;⑤忽略燃烧区内液相推进剂所占体积[6-8]。
根据以上假设以及质量、能量守恒方程、燃气状态方程可得到燃烧区状态方程如下
$$\dfrac{{{\rm{d}}{p_c}}}{{{\rm{d}}t}} = \dfrac{\gamma }{{{V_c}}}\left[ {{{\left( {RT} \right)}_{\rm{g}}}{{\dot m}_{\rm{g}}} - {{\left( {RT} \right)}_c}{{\dot m}_{{\rm{out}}}}} \right]$$ (24) $$\dfrac{{{\rm{d}}{{\left( {RT} \right)}_c}}}{{{\rm{d}}t}} = \dfrac{{{{\left( {RT} \right)}_c}}}{{{p_c}{V_c}}}\left\{ {\left[ {\gamma {{\left( {RT} \right)}_{\rm{g}}} - {{\left( {RT} \right)}_c}} \right]{{\dot m}_{\rm{g}}} - \left( {\gamma - 1} \right){{\left( {RT} \right)}_c}{{\dot m}_{{\rm{out}}}}} \right\}$$ (25) $$ \dfrac{{\rm{d}}{K}_{c}}{{\rm d}t}=\dfrac{{K}_{c}+1}{{\rho }_{c}{V}_{c}}\left({\dot{m}}_{\rm{in}}^{o}-{K}_{c}\centerdot {\dot{m}}_{\rm{in}}^{f}\right)$$ (26) 其中:
${\dot m_{\rm{g}}}$ 为单位时间生成燃气质量;${\dot m_{{\rm{out}}}}$ 为单位时间流出燃烧区的燃气质量;${p_c}$ 为燃烧区压力;${\left( {RT} \right)_{\rm{g}}}$ 、${\left( {RT} \right)_c}$ 分别为生成燃气、滞留燃气的热值;${V_c}$ 为燃烧区体积;${K_c}$ 为燃烧区平均质量混合比;${\rho _{_c}}$ 为燃烧区燃气平均密度;$\dot m_{{\rm{in}}}^o$ 、$\dot m_{{\rm{in}}}^f$ 为进入燃烧区的氧化剂、燃料的质量流量。 -
阀门、节流孔板等节流元件,可用统一的压力–流量特性方程描述
$$ \dot{m}=\tau \centerdot \overline{{C}_{d}A}\sqrt{2{\rho }_{\rm l}\left({p}_{\rm 1}-{p}_{2}\right)}$$ (27) 其中:
$\tau $ 为阀门相对开度;$\overline {{C_d}A} $ 为阀门额定流通截面积与流量系数的乘积;$\;{\rho _{\rm{l}}}$ 为流体密度;${p_1}$ 为节流元件入口压力;${p_2}$ 为节流元件出口压力。 -
大推力氢氧发动机系统方案如图3所示,搭建的发动机故障模式动态特性仿真模型,如图4所示。液氢与液氧经离心泵增压后,大部分进入燃烧室进行燃烧,然后通过喷管排出产生推力,小部分进入燃气发生器,燃烧产生的燃气驱动涡轮,分别带动氢泵和氧泵做功,然后直接排向大气[8]。在仿真模型中,使用与发动机实际喷管喉部面积相同的音速喷嘴模拟喷管的流量特性,同时为了提高计算效率,对系统模型中不影响故障仿真正确性的因素进行简化,忽略推力室冷却通道的换热,并用节流元件模拟冷却通道压降,这是由于推力室燃烧模型采用绝热模型,因此再生冷却模型中的夹套传热过程并不影响燃烧结果计算的准确性;另一方面将推进剂视为温度不变的弱可压流体[9-12]。
-
为验证上述系统仿真模型的准确性,对该型大推力氢氧发动机地面试车的起动与关机过程进行了动态特性仿真,仿真结果与试车实测参数的对比如图5所示。
通过对比可以看出,本文搭建的发动机系统动态特性仿真模型能够较为准确地反映发动机起动、关机过程中的参数变化趋势,推力室压力曲线与试车基本符合实际试车情况。由于建模过程中对各组件数学模型进行了不同程度的简化,因此仿真结果与试车实测参数存在一定的稳态误差,但是动态趋势的一致性较好。由此验证了上述发动机动态特性仿真模型的准确性。
-
采用前述仿真模型,对大推力氢氧发动机某试验故障多种可能的故障模式[13-18]进行了动态特性仿真分析,并与故障实测参数进行了详细对比。下面给出了几种典型故障模式的仿真结果分析及实测参数对比。
-
仿真时设置发生器氧阀异常关闭故障,阀门按线性关闭,响应时间50 ms,仿真得到的结果与故障实测参数对比如图6所示。
图 6发生器氧阀异常关闭故障仿真与实测曲线对比
Figure 6.Comparison between simulation results and measured data of abnormal closing of generator oxygen valve
从故障仿真结果看,如发生器氧阀异常关闭,则发动机各参数下降趋势与正常关机类似,由于燃气发生器的氧化剂供应被切断,其压力首先迅速下降,但因涡轮泵转子有惯性,因此推力室压力、泵后压力下降相对较为平缓,氢、氧涡轮泵转速下降速率也与之相似,这与故障实测参数迅速下降的现象差异明显。
-
仿真时设置氧涡轮排气流路流通面积减小,相对流通面积在50 ms内按线性方式减小95%。变化曲线如图7所示。
得到的仿真结果与实测参数曲线对比如图8。该故障模式下,仿真得到的发生器压力、推力室压力、氢/氧泵后压力及转速曲线均出现快速下降,趋势与实测参数基本一致,但氢泵转速下降曲线与实测参数曲线之间存在一定的差异,这可能是由于仿真模型与在非正常关机过程中剧烈变化的液氢的实际物性参数之间的差异造成的。
-
仿真时设置燃气喷嘴堵塞模式,对主燃气喷嘴堵塞的情况进行了模拟仿真。仿真时假设主燃气喷嘴流通面积在50 ms内线性减小至0,仿真结果与实测参数曲线对比如图9所示。
图 9主燃气喷嘴堵塞故障仿真与故障实测对比
Figure 9.Comparison between simulation results and measured data of main gas nozzle clogging
仿真结果表明,如果燃气喷嘴堵塞,燃气发生器首先出现一压力峰,随后该压力再迅速下降,氢涡轮泵转速也会有先上升、再迅速下降的趋势。这主要是由于燃气发生器及燃气喷嘴前、氢涡轮前的燃气管路构成了一个大的气体容腔,集中容积的特性较显著,燃气喷嘴堵塞,最先反应的即是上游的燃气发生器压力及氢涡轮入口压力(均出现压力峰),该压力峰未能传播至燃气喷嘴下游的氧涡轮入口。发动机其它参数变化相对滞后。
-
仿真时设置燃气发生器突然熄火,转为气氢膨胀做功状态,假设气氢温度为300 K。仿真结果与实测参数的对比情况如图10所示。
图 10燃气发生器熄火故障仿真与故障实测对比
Figure 10.Comparison between simulation results and measured data of gas generator flameout
仿真结果显示,燃气发生器熄火时,发生器压力迅速下降,比故障实测燃气发生器压力下降速率更快,但与此同时,受到氢氧涡轮泵惯性的影响,推力室压力、泵后压力以及氢氧涡轮泵转速的下降速度相对较缓,与故障实测曲线有明显的差别。
-
仿真时,待发动机工况稳定后,在氧涡轮泵转子上增加一个与泵扭矩方向相同的附加扭矩,大小为1 000
${\rm{N}} \cdot {\rm{m}}$ ,得到的仿真结果如图11所示。图 11氧涡轮泵转子卡滞故障仿真与故障实测对比
Figure 11.Comparison between simulation results and fault measurement data of oxygen turbine pump rotor sticking
可以看到,氧涡轮泵转子负载力矩增大1 000
${\rm{N}} \cdot {\rm{m}}$ 后,发动机参数迅速下降,燃烧装置压力、泵后压力等下降速率与实测参数基本一致,但氢泵转速的下降速率与实测参数差别较大,氧泵转速的下降速率也更大,说明在实际的故障中,氢氧涡轮泵的转动没有受到明显的外力阻碍。 -
通过仿真计算结果与实测参数的对比,发现在不同的故障模式中,氧涡轮排气管路阻塞的故障模式与发动机故障实测参数下降情况符合得最好。
-
通过实测数据分析、历史故障数据分析,结合前述故障模式仿真结果,初步确定导致发动机故障的最大可能原因是氧涡轮排气管路阻塞故障模式,为此开展了验证试验。
为模拟氧涡轮出口排气管阻塞,对大推力氢氧发动机进行了局部改造,在排气管路并联安装两个常开电爆阀。模拟氧涡轮出口阻塞故障时,两路电爆阀同时关闭,其中一路全关,另一路保留部分面积,最终预计阀门关闭后的当量面积与3.3节中仿真所用的管路堵塞95%一致。
通过热试验得到了氧涡轮出口排气管阻塞时,发动机整体参数的变化情况。结果表明,故障注入后(电爆阀起爆),氧泵后压力变化速率最大,整体参数变化趋势与故障实测及仿真结果相似,但下降速率仍略小于实测结果,主要是由于氧涡轮输出功率下降速率较慢所致,可导致该现象的因素包括电爆阀关闭后氧涡轮压比偏差、低压比下涡轮效率偏差等。故障模式验证试验与故障实测、仿真结果的对比如图12所示。
图 12故障模式验证试验与故障实测、仿真结果的对比
Figure 12.Comparison between fault mode verification test with fault measurement and simulation results
从图12可以看出,在氧涡轮出口阻塞故障模式下,发动机主要参数的动态变化趋势与故障实测参数及仿真计算结果基本一致,初步验证了故障仿真结果及故障定位的可信度。
根据仿真与试验结果,分析故障可能的发展过程为氧涡轮排气管路在复杂力热环境下局部结构出现异常,造成氧涡轮排气阻塞,导致氧涡轮后压力升高,涡轮压比减小,氧涡轮失去做功能力;随后发动机失去氧化剂供应,燃气发生器与推力室相继熄火,发动机主要参数快速下降,发动机失去推力。
-
本文建立了液体火箭发动机各组件的动态特性数学模型,在此基础上构建了大推力氢氧发动机系统动态特性仿真模型;将该模型仿真获得的结果与地面试车实测参数进行对比,验证了仿真模型的准确性。
针对大推力氢氧发动机某次试验故障,采用该仿真模型对多种可能的故障模式开展仿真计算,获取了不同故障模式下发动机参数的变化规律,通过对比分析,认为可能性较大的故障模式为氧涡轮排气管路阻塞,随后针对该故障模式开展验证试验,结果吻合良好,证明通过模型仿真实现了故障的准确定位。
综上,本文建立的动态仿真模型可以实现大推力氢氧发动机的故障识别与准确定位,在此基础上,利用故障树形成完整的故障注入模型阵列,通过仿真计算建立完整的故障参数特征库,就能够通过发动机实时工作参数与故障参数特征库的对比,实现发动机故障的实时监测与准确定位,提升发动机的工作可靠性。
Fault Simulation and Experimental Study on High-Thrust LOX/LH2 Rocket Engine
-
摘要:为对大推力氢氧发动机的某次试验开展故障定位,建立了液体火箭发动机动态特性仿真模型库,搭建了大推力氢氧发动机动态特性仿真模型,并通过与试验实测参数对比,验证了仿真模型的准确性。针对发动机可能出现的故障,依次开展仿真计算,将获得的结果与故障实测参数进行对比,得到最有可能的故障模式为氧涡轮排气管路阻塞。通过试验对仿真结果开展验证,结果吻合良好,证明了故障定位的准确性与通过动态仿真实现大推力氢氧发动机故障定位的方法可行性,为未来发展大推力氢氧发动机故障监测与定位技术奠定了基础。Abstract:In order to locate the fault in a test of a high-thrust LOX/LH2 engine, the dynamic characteristic simulation model library of liquid rocket engine is established, and the dynamic characteristic simulation model of the high-thrust LOX/LH2 engine is built. The accuracy of the simulation model is verified by comparing the model with the experimental parameters. According to the possible faults of the engine, the simulation calculation is carried out in turn, and the results are compared with the measured parameters. Through the comparison, it is found that the most likely failure mode is the blockage of the exhaust pipe of the oxygen turbine. The simulation results verified by experiments were in good agreement, which not only proves the accuracy of fault location, but also proves the feasibility of fault location of high-thrust oxyhydrogen engine by dynamic simulation, laying a foundation for further development of fault monitoring and location technology for high-thrust LOX/LH2 rocket engine in the future.Highlights
● The dynamic characteristics simulation model of high thrust LOX/LH2 rocket engine system is established. ● The response characteristics of the parameters of the LOX/LH2 rocket engine system in different failure modes are obtained. ● The simulation results are in good agreement with the experimental results,which verifies the correctness of the method. -
-
[1] 张育林, 刘昆, 程谋森. 液体火箭发动机动力学理论与应用[M]. 北京: 科学出版社, 2005. [2] 熊靖宇. 液氧甲烷发动机故障特征分析与故障预测方法研究[D]. 长沙: 国防科技大学, 2018. [3] 何涛,黄敏超,胡小平,等. 某火箭发动机故障检测及诊断算法设计分析[J]. 南京航空航天大学学报,2019(S1):50-55.HE T,HUANG M C,HU X P,et al. Design analysis of fault detection and diagnosis algorithms for rocket engine[J]. Journal of Nanjing University of Aeronautics & Astronautics,2019(S1):50-55. [4] 黄敏超,刘昆,张育林,等. 分级燃烧循环火箭发动机的动态特性研究[J]. 航空动力学报,2002(4):500-504.doi:10.3969/j.issn.1000-8055.2002.04.023HUANG M C,LIU K,ZHANG Y L,et al. Dynamic characteristic study of staged combustion cycle rocket engine[J]. Journal of Aerospace Power,2002(4):500-504.doi:10.3969/j.issn.1000-8055.2002.04.023 [5] 樊忠泽,黄敏超. 空间推进系统工作过程故障仿真[J]. 国防科技大学学报,2008,30(2):11-15.PAN Z Z,HUANG M C. Fault simulation of space power system in the operation process[J]. Journal of National University of Defense Technology,2008,30(2):11-15. [6] 曹泰岳. 火箭发动机动力学[M]. 长沙: 国防科技大学出版社, 2004. [7] 吴义忠, 陈立平. 多领域物理系统的仿真优化方法[M]. 北京: 科学出版社, 2011. [8] 蔡国飙, 李家文, 田爱梅. 液体火箭发动机设计[M]. 北京: 北京航空航天大学出版社, 2011. [9] 王建波,于达仁,王广雄. 液体火箭发动机泄漏故障实时仿真[J]. 推进技术,1999,20(5):1-5.doi:10.3321/j.issn:1001-4055.1999.05.001WANG J B,YU D R,WANG G X. Real time simulation of leak fault of liquid rocket engine[J]. Journal of Propulsion Technology,1999,20(5):1-5.doi:10.3321/j.issn:1001-4055.1999.05.001 [10] 吴建军. 液体火箭发动机故障检测与诊断研究[D]. 长沙: 国防科学技术大学, 1995. [11] 潘辉,张黎辉. AMESim软件在液体火箭发动机系统动态仿真上的应用[J]. 火箭推进,2011,37(3):6-11.doi:10.3969/j.issn.1672-9374.2011.03.002PAN H,ZHANG L H. Application of AMESim in dynamic characteristic simulation of liquid rocket engine system[J]. Journal of Rocket Propulsion,2011,37(3):6-11.doi:10.3969/j.issn.1672-9374.2011.03.002 [12] 王弘亚,张卫红. 基于AMESim平台的氢氧火箭发动机启动过程仿真研究[J]. 导弹与航天运载技术,2017,357(6):28-31.WANG H Y,ZHANG W H. Research on simulation of start-up process of LOX/LH2 rocket engine based on AMESim[J]. Missiles and Space Vehicles,2017,357(6):28-31. [13] YAMANISHI N, KIMURA T, TAKAHASHI M, et al. Transient analysis of the LE-7A rocket engine using the rocket engine dynamic simulator(REDS)[C]//40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. [S.l]: AIAA, 2004. [14] CHRISTENSEN E, BROWN A, FRADY G. Calculation of dynamic loads due to random vibration environemnts in rocket engine systems[C]//AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, & Materials Conference. [S.l]: AIAA, 2013. [15] JUN S I, TOKUMASU T, KAMIJO K. Analysis of dynamic response of rocket engine LH_2 feed system[C]//The Japan Society of Mechanical Engineers. Japan: JSME, 2002. [16] RUTH E, AHN H, BAKER R, et al. Advanced liquid rocket engine transient model[C]//Joint Propulsion Conference. [S.l]: 26th Joint Propulsion Conference, 2013. [17] MATTEO F D, ROSA M D, ONOFRI M. Start-Up transient simulation of a liquid rocket engine[C]//AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. [S.l]: AIAA, 2011. [18] 王珏, 朱森元. 氢氧火箭发动机(YF-73)动态数学模型及其数值模拟方法[C]//中国宇航学会液体火箭推进专业委员会学术会议. 青岛: 中国宇航学会, 2007. -