李纯飞,龙春伟,董 苑
某单组元发动机系统配套两台贮箱供应推进剂,两台贮箱采用并联方式布局,接入环形管路将推进剂分配至各推力装置,贮箱与推力装置布局结构如图1所示。系统工作时,两贮箱推进剂供应流量不会完全一致,导致贮箱内推进剂排放不均衡,严重时其中一台贮箱推进剂耗尽而另一台贮箱剩余大量推进剂,会影响系统工作性能和可靠性,同时引起飞行器质心偏移,影响控制精度。经梳理,发动机管路流阻特性(连接贮箱与环管管路流阻大小及偏差)、推力装置工作特性(如推力装置喷注压降、床流阻、分解效率的各异性)、发动机结构布局(分布均匀性)等因素影响系统排放均衡性。为定量获得各影响因素进而为系统采取控制措施提供依据,建立了计算模型,对各影响因素进行了计算分析,搭建了试验验证系统。
图1 贮箱并联布局示意图Fig.1 Layout diagram of parallel tank
以两台推力装置同时工作为例,建立了如图2所示分析模型,图2中:A,B为两贮箱;F1,F2为两台推力装置;qA,qB为两贮箱流量;qA1,qA2为A贮箱供应的推进剂经环管分流后供应给两台推力装置的流量;qB1,qB2为B贮箱供应的推进剂经环管分流后供应给两台推力装置的流量;q1,q2为两台推力装置的流量;per为贮箱入口增压压力(采用同一减压阀供应气体,因此假设两贮箱入口压力相等);pA为A贮箱对应环管入口压力;pB为B贮箱对应环管入口压力;p1,p2为两台推力装置入口压力。
图2 模型示意图Fig.2 Diagram of analysis model
系统工作时,推力装置推进剂供应遵循压力平衡和流量平衡,以图2所示系统为例,可建立如下压力平衡和流量平衡方程组:
(1)
(2)
(3)
(4)
qA=qA1+qA2
(5)
qB=qB1+qB2
(6)
q1=qA1+qB1
(7)
q2=qA2+qB2
(8)
对于推进剂供应系统,推进剂增压压力即per为已知,贮箱与管路系统流阻特性通过液流试验或计算模型计算获得,流阻特性可拟合为流量的一元关系式;发动机入口压力与流量关系通过热试车获得,为此,推进剂供应压力平衡方程式(1)~(4)可转化为与流量相关的方程组,见式(9)~(12),结合系统流量平衡方程式(5)~(8),通过迭代计算,即可实现系统流量求取,进而可获得两贮箱供应流量并定量分析对排放均衡性的影响:
per=fq1+fqA1+fqA
(9)
per=fq1+fqA2+fqA
(10)
per=fq1+fqB1+fqB
(11)
per=fq1+fqB2+fqB
(12)
对管路流阻特性、推力装置工作特性、发动机结构布局等影响因素进行了计算。
管路流阻特性影响计算时重点计算了连接贮箱与环管管路流阻大小、流阻偏差对系统排放均衡性的影响,流阻大小变化通过改变管路通径获得,流阻偏差通过调节管路流阻相关系数获得。
推力装置工作特性影响计算时通过改变入口压力与流量拟合关系式,实现不同推力装置工作特性模拟,计算了工作特性无差异以及工作特性存在差异两状态对系统排放均衡性的影响。
发动机结构布局影响计算时对发动机与贮箱对称布局状态和非对称布局两状态下系统供应特性进行了计算。
计算了3种系统状态(系统状态1,2,3),8种工作模式(组合A,B,C,D,E,F,G,H)下系统供应特性。系统状态描述见表1。8种工作模式均为4机同时工作,组合A中发动机推力品种、推力装置工作特性相同,布局相同,贮箱对称布局;组合C,E,G三种工作模式推力品种相同、推力装置工作特性相同,结构布局不同,贮箱为非对称布局;组合B相对组合A,组合D相对组合C,组合F相对组合E以及组合H相对组合G除靠近A贮箱两台推力装置工作特性不同外,其余状态相同;组合B,D,F,H中靠近A贮箱两台发动机在设计入口压力下流量为设计值的1.1倍。计算结果见表1。
.
1.4.1发动机管路流阻特性影响分析
对比系统状态1与2,状态1连接A,B贮箱与环管两管路流阻特性相同,均为设计状态,状态2两管路流阻存在偏差,连接A贮箱管路流阻较设计状态偏大20%,连接B贮箱管路流阻为设计状态。A,C,E,G四种工作组合在状态1下两贮箱排液偏差分别为0.0%,+1.0%,-4.1%,+4.3%,在状态2下两贮箱排液偏差分别为-7.5%,-6.9%,-11.2%,-3.9%,两状态下4种工作组合两贮箱排液偏差值增加7.1%~8.2%,变化较大,表明连接A、B两贮箱与环管管路的流阻偏差对系统排放均衡性影响较大。对比状态1与2四种工作组合计算结果,两管路流阻相等时,系统排放具有较好的均衡性。
对比系统状态1与3,状态1连接A,B贮箱与环管两管路通径为设计状态,状态3两管路通径较设计状态减小了20%,相应的阻尼效应增强。根据计算结果,A,C,E,G四种工作组合在状态3下两贮箱排液偏差分别为0.0%,+0.3%,-1.5%,+1.5%,相比于状态1,除对组合A排液偏差无影响外,组合B,C,D排液偏差均减小,减小0.7%~2.8%,变化较大,表明连接A,B两贮箱与环管管路的流阻大小对系统排放均衡性影响较大。对比状态1与3四种工作组合计算结果,增加两管路设计流阻有利于提高系统排放均衡性。
1.4.2推力装置工作特性影响分析
组合B,D,F及H中靠近A贮箱的两台推力装置工作特性相对设计状态存在变化,在设计入口压力下,其流量取为设计流量的1.1倍。根据计算结果,上述4种组合在状态1下两贮箱排液偏差分别为0.1%,+1.1%,-4.0%及+4.4%,对比A,C,E,G四种工作组合,各组合两贮箱排液偏差值增加约0.1%,相对推力装置工作特性变化10%,两贮箱排液偏差值变化较小,表明推力装置工作特性对系统排放均衡性影响有限。
1.4.3发动机结构布局影响分析
组合C,E,G推力装置品种及数量相同,结构布局不同,对比3种工作组合在3种计算状态下计算结果,各计算状态下,3种工作组合两贮箱排液偏差不同,最大相差3.0%~8.4%不等,表明发动机结构布局对系统排放均衡性影响较大。组合A发动机布局相对两贮箱对称,总体而言,在各计算状态下,组合A下两贮箱排液偏差较小,表明发动机宜设计为对称布局方案,以改善系统排放均衡性。
根据计算分析表明,影响本发动机系统排放均衡性的主导因素有连接贮箱与环管入口管路流阻大小与偏差、发动机结构布局,推力装置工作特性影响有限。
搭建冷态试验系统,将其介质供应流量-压降特性与发动机推进剂供应系统流量-压降特性设计一致,即可实现对发动机系统推进剂供应过程模拟,所建立的试验系统排放均衡性与发动机系统相当。试验系统进行冷态排液试验,获得各工作模式下系统排放特性,进而可对各影响因素影响情况进行试验验证。
试验时各工作模式逐一进行试验,各工作模式单独进行一定时长的排液,通过获得试验排液零秒及排液末秒两贮箱质量,即可获得各工作模式下两贮箱排液量及偏差情况,进而获得各影响因素影响情况。
试验选取了两种试验状态(试验状态1,2)、5种工作模式(组合U,V,W,X,Y)。试验状态描述见表2。组合U,V推力品种、数量、布局相同,推力装置工作特性不同,组合U为设计状态,组合V中靠近B贮箱的两台发动机在设计入口压力下流量为设计状态的1.2倍,另外两台靠近A贮箱的发动机为设计状态;组合W,X推力品种、数量与工作特性相同,布局不同。试验结果见表2。
2.2.1发动机管路流阻特性影响分析
对比试验状态1与2,试验状态2连接贮箱与环管管路在设计状态基础上设置了节流孔板,流阻增加。组合U,V,W及Y均在两状态下进行了试验,根据试验结果,4种组合在试验状态1下两贮箱排液偏差分别为+1.4%,+0.4%,+10.2%及+18.4%,在试验状态2下分别为+0.4%,-0.5%,+3.6%及+7.5%。由于组合U,V排放均衡性相对较好,在两试验状态下变化不明显;组合W,Y排放均衡性相对较差,两组合在试验状态2下排放均衡性提高明显,相对试验状态1偏差分别减小6.6%,10.9%,表明增加连接贮箱与环管管路的流阻特性有利于提高系统排放均衡性。
.
2.2.2推力装置工作特性影响分析
组合U与V除推力装置工作特性不同外,其余状态均相同,其中组合V中两台靠近B贮箱的发动机在设计入口压力下流量为设计状态的1.2倍。根据试验结果,两组合在试验状态1下两贮箱排液偏差分别为+1.4%、+0.4%,相差1%,相对推力装置工作特性存在20%变化,排液偏差变化很小,表明推力装置工作特性变化对系统排放均衡性的影响有限。
2.2.3发动机结构布局影响分析
组合W与X除发动机结构布局不同外,其余状态均相同。组合W在试验状态1下两贮箱排液偏差为+10.2%,组合X在试验状态1下两贮箱排液偏差为-1.2%,两组合下系统排液偏差相差11.4%,表明发动机结构布局对系统排放均衡性具有较大影响。
根据试验结果,影响本发动机系统排放均衡性的主导因素有连接贮箱与环管入口管路流阻大小、发动机结构布局,推力装置工作特性影响有限。
仿真分析和试验验证结果表明,连接贮箱与环管入口管路流阻大小与偏差、发动机结构布局对系统排放均衡性影响较大,推力装置工作特性影响有限。
可采取的控制措施:
1)在连接贮箱与环管入口管路设置一定压降的限流结构或适当减小连接贮箱与环管入口管路的通径,增加管路阻尼效应。
2)对连接贮箱与环管入口管路进行液流试验,通过选配产品进而消除或降低两管路流阻差异。
3)发动机相对两贮箱对称布局;将发动机接口设置于远离两贮箱的环管中部位置,或在环管中部设置转接通集中安装各发动机分机,以降低不对称性影响。
对影响发动机贮箱推进剂排放均衡性的主要因素,进行了计算分析与试验,计算与试验结果表明,连接贮箱至环管的管路流阻大小、两管路流阻偏差、各推力装置相对两贮箱布局是否对称对发动机系统排放均衡性具有较大影响,推力装置工作特性影响有限,在此基础上,提出了可行的控制措施。