CN105628387B - 采用预测控制法调节试车台进气状态参数的调试方法 - Google Patents
采用预测控制法调节试车台进气状态参数的调试方法 Download PDFInfo
- Publication number
- CN105628387B CN105628387B CN201511017540.1A CN201511017540A CN105628387B CN 105628387 B CN105628387 B CN 105628387B CN 201511017540 A CN201511017540 A CN 201511017540A CN 105628387 B CN105628387 B CN 105628387B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- valve
- pressure
- inlet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种采用预测控制法调节试车台进气状态参数的调试方法,所述调试的类型包括空载调试、带载调试,所述带载调试,包括如下步骤:按预测控制法得出的单个参数的控制量,依次调试进气压力调节阀,加热器燃料比例调节阀,发动机燃料比例节流阀,液氧流量调节阀;与现有技术相比,本发明提供的采用预测控制法调节试车台进气状态参数的调试方法,可以整体掌握试车台的性能,降低试验成本,提高试验效率。
Description
技术领域
本发明涉及一种采用预测控制法调节试车台进气状态参数的调试方法。
背景技术
目前,发动机试验时,要求状态点的一组进气参数,包括进气总压,进气空气流量,进气氧成分,进气总温,发动机供油油量,同时达到设定值及稳定后,记录其性能及参数或考核其性能。但是,在发动机试验前,需要对试车台进行调试,以确定试车台的性能是否能够满足发动机的试验要求,在试车台调试中,进气***的调试是重要组成部分,但是,对于采用预测控制法调节试车台进气状态参数的调试方法,目前尚没有成熟的调试方法,这导致了无法整体掌握试车台的性能,从而增加了模拟时间,提高了试验成本。
发明内容
本发明的目的是为了克服上述现有技术存在的缺陷而提供了一种采用预测控制法调节试车台进气状态参数的调试方法。
本发明的目的可以通过以下技术方案来实现:
一种采用预测控制法调节试车台进气状态参数的调试方法,所述调试的类型包括空载调试、带载调试,所述带载调试,包括如下步骤:
1)、按预测控制法得出的单个参数的控制量,调试进气压力调节阀;
2)、按预测控制法得出的单个参数的控制量,调试加热器燃料比例调节阀
3)、按预测控制法得出的单个参数的控制量,调试发动机燃料比例节流阀
4)、按预测控制法得出的单个参数的控制量,调试液氧流量调节阀。
所述预测控制法,包括如下步骤:
1)、根据当前状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量,与所需模拟状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量的设定值的差值通过进气***状态点参数方程组解算出单个参数的控制量;
2)、根据所述步骤1)中所解算出的单个参数的控制量,调节控制上述参数的各阀门的开度,使得在下一循环开始时,设定各阀门的开度满足所述步骤1)中所解算出的单个参数的控制量;
3)重复所述步骤1)和步骤2)直至模拟精度满足要求;
所述进气***状态点参数方程组为:
qmi=qma+qmo+qmj
pti=pt2+Δpt2
所述进气***状态点参数方程组中,qmi为进气流量,qma为通过进气压力调节阀的质量流量,qmo为补氧流量,qmj为加热器供油流量,pti为进气总压,Ai为进气截面面积,Tti为进气总温,q(λi)为计算截面气体流量系数,Tti0为当前进气总温,η为燃烧效率,Qj为加热器供油汽化潜热,α为过热空气系数,Cp为燃气定压比热,Cv为进气压力调节阀的阀门流量系数,Psa为气源压力,pt2为空气调节阀阀后压力,γ为空气质量比重,T0a为空气源温度,ΔPt2为管道压力损失,Tsa为混合前空气温度,Cpa为空气定压比热,Tso为混合前补氧温度,Cpo为氧气定压比热,Q0为加热器供油汽化潜热,为理论空燃比;
在所述进气***状态点参数方程组中,已知当前态点的进气空气流量,当前态点的补氧流量,当前态点的加热器供油流量,当前态点的进气总压,当前态点的进气总温,当前态点的当前进气总温,当前态点的过热空气系数,和所需模拟状态点的进气总压,所需模拟状态点的进气总温,求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量,所需模拟状态点的补氧流量,所需模拟状态点的加热器供油流量。
所述进气压力调节阀的调试内容包括确定进气压力调节阀的流量-开度特性曲线、临界点、分辨率、重复性精度、负载压力充填时间验证、负载压降曲线;所述进气压力调节阀的分辨率为使流量或进气压力发生有效、最小变化的调节阀门开度变化值,所述负载压降为进气压力调节阀阀后压力与进气模拟总压的差值。
确定进气压力调节阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过进气压力调节阀的质量流量与阀前压力和阀后压力,介质密度,温度、阀门额定面积、阀门开度的关系,由公式
时
时
和qma=Qρa得到:
上述公式中,CV为进气压力调节阀的阀门流量系数,Q为通过进气压力调节阀的体积流量,γ1为介质比重,T0为空气温度,p1为进气压力调节阀阀前压力,p2为进气压力调节阀阀后压力,qma为通过进气压力调节阀的质量流量,ρa为介质密度,Ka为进气压力调节阀的阀门开度,Aa为进气压力调节阀的阀门额定面积,cp1为介质定压比热,R为气体常数,Ta为介质温度,Δp为进气压力调节阀阀后压力与阀前压力的压差;
其中,进气压力调节阀的阀门额定面积,介质密度,介质定压比热,介质比重为常数;
2)、将测量所得的介质温度,进气压力调节阀阀前压力,进气压力调节阀阀后压力与阀前压力的压差作为参变量,将进气压力调节阀的阀门开度作为自变量;
3)以气源设计压力的85%为上限,以所需最小的进气压力调节阀阀前压力为下限,将由所需最小的进气压力调节阀阀前压力和气源设计压力的85%之间的压力分割成若干压力点,在每个压力点下,将进气压力调节阀按10%~90%、间隔10%的开度正反行程开启三遍,记录每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量,将上述记录的每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量的3个数据取平均数,拟合进气压力调节阀流量-开度特性曲线。
确定进气压力调节阀临界点的方法,其步骤为:在进气压力调节阀阀前压力固定和进气压力调节阀的阀门开度变化过程中,在某个进气压力调节阀的阀门开度下出现临界流动时,然后逐步减小进气压力调节阀的阀门开度,直至通过进气压力调节阀的质量流量与进气压力调节阀阀前压力不成正比关系,由此点的进气压力调节阀阀后压力即可确定进气压力调节阀临界点。
发动机燃料比例节流阀的调试方法的调试内容包括确定所述发动机燃料比例节流阀的流量-开度曲线、分辨率和重复性精度;所述发动机燃料比例节流阀的分辨率为使流量发生有效、最小变化的阀门开度变化值。
确定发动机燃料比例节流阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过发动机燃料比例节流阀的质量流量仅与发动机燃料比例节流阀阀后压力与阀前压力的压差、工作介质密度、阀门额定面积、阀门开度有关,由公式
确定发动机燃料比例节流阀的质量流量与发动机燃料比例节流阀的阀门开度的关系;
上述公式中,qm为通过发动机燃料比例节流阀的质量流量,Kb为发动机燃料比例节流阀的阀门开度,CVj为发动机燃料比例节流阀的阀门流量系数,Ab为发动机燃料比例节流阀阀门额定面积,Δpi为发动机燃料比例节流阀阀后压力与阀前压力的压差,ρb为工作介质密度,g为重力加速度,其为已知常数;
2)、将发动机燃料比例节流阀的阀门额定面积、工作介质密度作为常数,将发动机燃料比例节流阀阀后压力与阀前压力的压差作为参变量,将发动机燃料比例节流阀的阀门开度作为自变量;
3)、在试车台进气***的任一工作压力下,均会产生一个发动机燃料比例节流阀阀后压力与阀前压力的压差,对于任一个发动机燃料比例节流阀阀后压力与阀前压力的压差,将发动机燃料比例节流阀的阀门开度从5%、10%、15%、...、90%、95%开启三遍,进行喷液试验,记录通过发动机燃料比例节流阀的质量流量,将每个阀门开度所记录的三个发动机燃料比例节流阀的质量流量取平均值,作为有效的发动机燃料比例节流阀的质量流量,进而计算出发动机燃料比例节流阀流量系数,同时拟合发动机燃料比例节流阀的流量-开度曲线。
所述空载调试和带载调试可以用模拟筒代替。
所述预测控制法用于直连式试车台。
所述预测控制法可用于直连式试车台用于发动机状态点进气模拟。
本发明的有益效果是,与现有技术相比,本发明提供的采用预测控制法调节试车台进气状态参数的调试方法,可以整体掌握试车台的性能,降低试验成本,提高试验效率。
具体实施方式
一种采用预测控制法调节试车台进气状态参数的调试方法,所述调试的类型包括空载调试、带载调试,所述带载调试,包括如下步骤:
1)、按预测控制法得出的单个参数的控制量,调试进气压力调节阀;
2)、按预测控制法得出的单个参数的控制量,调试加热器燃料比例调节阀
3)、按预测控制法得出的单个参数的控制量,调试发动机燃料比例节流阀
4)、按预测控制法得出的单个参数的控制量,调试液氧流量调节阀。
所述预测控制法,包括如下步骤:
1)、根据当前状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量,与所需模拟状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量的设定值的差值通过进气***状态点参数方程组解算出单个参数的控制量;
2)、根据所述步骤1)中所解算出的单个参数的控制量,调节控制上述参数的各阀门的开度,使得在下一循环开始时,设定各阀门的开度满足所述步骤1)中所解算出的单个参数的控制量;
3)重复所述步骤1)和步骤2)直至模拟精度满足要求;
所述进气***状态点参数方程组为:
qmi=qma+qmo+qmj
pti=pt2+Δpt2
所述进气***状态点参数方程组中,qmi为进气流量,qma为通过进气压力调节阀的质量流量,qmo为补氧流量,qmj为加热器供油流量,pti为进气总压,Ai为进气截面面积,Tti为进气总温,q(λi)为计算截面气体流量系数,Tti0为当前进气总温,η为燃烧效率,Qj为加热器供油汽化潜热,α为过热空气系数,Cp为燃气定压比热,Cv为进气压力调节阀的阀门流量系数,Psa为气源压力,pt2为空气调节阀阀后压力,γ为空气质量比重,T0a为空气源温度,ΔPt2为管道压力损失,Tsa为混合前空气温度,Cpa为空气定压比热,Tso为混合前补氧温度,Cpo为氧气定压比热,Q0为加热器供油汽化潜热,为理论空燃比;
在所述进气***状态点参数方程组中,已知当前态点的进气空气流量,当前态点的补氧流量,当前态点的加热器供油流量,当前态点的进气总压,当前态点的进气总温,当前态点的当前进气总温,当前态点的过热空气系数,和所需模拟状态点的进气总压,所需模拟状态点的进气总温,求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量,所需模拟状态点的补氧流量,所需模拟状态点的加热器供油流量。
所述进气压力调节阀的调试内容包括确定进气压力调节阀的流量-开度特性曲线、临界点、分辨率、重复性精度、负载压力充填时间验证、负载压降曲线;所述进气压力调节阀的分辨率为使流量或进气压力发生有效、最小变化的调节阀门开度变化值,所述负载压降为进气压力调节阀阀后压力与进气模拟总压的差值。
确定进气压力调节阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过进气压力调节阀的质量流量与阀前压力和阀后压力,介质密度,温度、阀门额定面积、阀门开度的关系,由公式
时
时
和qma=Qρa得到:
上述公式中,CV为进气压力调节阀的阀门流量系数,Q为通过进气压力调节阀的体积流量,γ1为介质比重,T0为空气温度,p1为进气压力调节阀阀前压力,p2为进气压力调节阀阀后压力,qma为通过进气压力调节阀的质量流量,ρa为介质密度,Ka为进气压力调节阀的阀门开度,Aa为进气压力调节阀的阀门额定面积,cp1为介质定压比热,R为气体常数,Ta为介质温度,Δp为进气压力调节阀阀后压力与阀前压力的压差;
其中,进气压力调节阀的阀门额定面积,介质密度,介质定压比热,介质比重为常数;
2)、将测量所得的介质温度,进气压力调节阀阀前压力,进气压力调节阀阀后压力与阀前压力的压差作为参变量,将进气压力调节阀的阀门开度作为自变量;
3)以气源设计压力的85%为上限,以所需最小的进气压力调节阀阀前压力为下限,将由所需最小的进气压力调节阀阀前压力和气源设计压力的85%之间的压力分割成若干压力点,在每个压力点下,将进气压力调节阀按10%~90%、间隔10%的开度正反行程开启三遍,记录每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量,将上述记录的每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量的3个数据取平均数,拟合进气压力调节阀流量-开度特性曲线。
确定进气压力调节阀临界点的方法,其步骤为:在进气压力调节阀阀前压力固定和进气压力调节阀的阀门开度变化过程中,在某个进气压力调节阀的阀门开度下出现临界流动时,然后逐步减小进气压力调节阀的阀门开度,直至通过进气压力调节阀的质量流量与进气压力调节阀阀前压力不成正比关系,由此点的进气压力调节阀阀后压力即可确定进气压力调节阀临界点。
发动机燃料比例节流阀的调试方法的调试内容包括确定所述发动机燃料比例节流阀的流量-开度曲线、分辨率和重复性精度;所述发动机燃料比例节流阀的分辨率为使流量发生有效、最小变化的阀门开度变化值。
确定发动机燃料比例节流阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过发动机燃料比例节流阀的质量流量仅与发动机燃料比例节流阀阀后压力与阀前压力的压差、工作介质密度、阀门额定面积、阀门开度有关,由公式
确定发动机燃料比例节流阀的质量流量与发动机燃料比例节流阀的阀门开度的关系;
上述公式中,qm为通过发动机燃料比例节流阀的质量流量,Kb为发动机燃料比例节流阀的阀门开度,CVj为发动机燃料比例节流阀的阀门流量系数,Ab为发动机燃料比例节流阀阀门额定面积,Δpj为发动机燃料比例节流阀阀后压力与阀前压力的压差,ρb为工作介质密度,g为重力加速度,其为已知常数;
2)、将发动机燃料比例节流阀的阀门额定面积、工作介质密度作为常数,将发动机燃料比例节流阀阀后压力与阀前压力的压差作为参变量,将发动机燃料比例节流阀的阀门开度作为自变量;
3)、在试车台进气***的任一工作压力下,均会产生一个发动机燃料比例节流阀阀后压力与阀前压力的压差,对于任一个发动机燃料比例节流阀阀后压力与阀前压力的压差,将发动机燃料比例节流阀的阀门开度从5%、10%、15%、...、90%、95%开启三遍,进行喷液试验,记录通过发动机燃料比例节流阀的质量流量,将每个阀门开度所记录的三个发动机燃料比例节流阀的质量流量取平均值,作为有效的发动机燃料比例节流阀的质量流量,进而计算出发动机燃料比例节流阀流量系数,同时拟合发动机燃料比例节流阀的流量-开度曲线。
所述空载调试和带载调试可以用模拟筒代替。
所述预测控制法用于直连式试车台。
所述预测控制法可用于直连式试车台用于发动机状态点进气模拟。
最后说明的是,以上实施例仅用于说明本发明的技术方案而非限制,尽管参照较佳的实施例对本发明进行了详细说明,本领域的技术人员应当理解,可以对本发明的技术方案进行修改和等同替代,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,所述调试的类型包括空载调试、带载调试,所述带载调试,包括如下步骤:
1)、按预测控制法得出的单个参数的控制量,调试进气压力调节阀;
2)、按预测控制法得出的单个参数的控制量,调试加热器燃料比例调节阀;
3)、按预测控制法得出的单个参数的控制量,调试发动机燃料比例节流阀;
4)、按预测控制法得出的单个参数的控制量,调试液氧流量调节阀;
所述预测控制法,包括如下步骤:
a)、根据当前状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量,与所需模拟状态点的进气总压,进气空气流量,补氧流量,进气总温,发动机供油油量的设定值的差值通过进气***状态点参数方程组解算出单个参数的控制量;
b)、根据所述步骤a)中所解算出的单个参数的控制量,调节控制上述参数的各阀门的开度,使得在下一循环开始时,设定各阀门的开度满足所述步骤a)中所解算出的单个参数的控制量;
c)重复所述步骤a)和步骤b)直至模拟精度满足要求;
所述进气***状态点参数方程组为:
qmi=qma+qmo+qmj
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.0404</mn>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
</mrow>
<msqrt>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
</mrow>
</msub>
</msqrt>
</mfrac>
<mi>q</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&eta;q</mi>
<mrow>
<mi>m</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>Q</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&alpha;q</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>287</mn>
<msub>
<mi>C</mi>
<mi>V</mi>
</msub>
<msqrt>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&gamma;T</mi>
<mrow>
<mn>0</mn>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
pti=pt2+△pt2
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>Q</mi>
<mi>o</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mi>a</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>L</mi>
<mi>j</mi>
<mn>0</mn>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>21</mn>
<mi>%</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>21</mn>
<mi>%</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
所述进气***状态点参数方程组中,qmi为进气空气流量,qma为通过进气压力调节阀的质量流量,qmo为补氧流量,qmj为加热器供油流量,pti为进气总压,Ai为进气截面面积,Tti为进气总温,q(λi)为计算截面气体流量系数,Tti0为当前进气总温,η为燃烧效率,Qj为加热器供油汽化潜热,α为过热空气系数,Cp为燃气定压比热,Cv为进气压力调节阀的阀门流量系数,Psa为气源压力,pt2为空气调节阀阀后压力,γ为空气质量比重,T0a为空气源温度,ΔPt2为管道压力损失,Tsa为混合前空气温度,Cpa为空气定压比热,Tso为混合前补氧温度,Cpo为氧气定压比热,Q0为加热器供油汽化潜热,为理论空燃比;
在所述进气***状态点参数方程组中,已知当前态点的进气空气流量,当前态点的补氧流量,当前态点的加热器供油流量,当前态点的进气总压,当前态点的进气总温,当前态点的当前进气总温,当前态点的过热空气系数,和所需模拟状态点的进气总压,所需模拟状态点的进气总温,求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量,所需模拟状态点的补氧流量,所需模拟状态点的加热器供油流量。
2.根据权利要求1所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,所述进气压力调节阀的调试内容包括确定进气压力调节阀的流量-开度特性曲线、临界点、分辨率、重复性精度、负载压力充填时间验证和负载压降曲线;所述进气压力调节阀的分辨率为使流量或进气压力发生有效、最小变化的调节阀门开度变化值,所述负载压降为进气压力调节阀阀后压力与进气模拟总压的差值。
3.根据权利要求2所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,确定进气压力调节阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过进气压力调节阀的质量流量与阀前压力和阀后压力,介质密度,温度、阀门额定面积、阀门开度的关系,由公式
(p1-p2)<p1/2时
(p1-p2)≥p1/2时
和qma=Qρa得到:
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<msub>
<mi>C</mi>
<mi>v</mi>
</msub>
<msub>
<mi>A</mi>
<mi>a</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>a</mi>
</msub>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>p</mi>
</mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>RT</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
</msqrt>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>p</mi>
<mo><</mo>
<msub>
<mi>&beta;</mi>
<mi>p</mi>
</msub>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<msub>
<mi>C</mi>
<mi>v</mi>
</msub>
<msub>
<mi>A</mi>
<mi>a</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>a</mi>
</msub>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>RT</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
</msqrt>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>p</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>&beta;</mi>
<mi>p</mi>
</msub>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
上述公式中,CV为进气压力调节阀的阀门流量系数,Q为通过进气压力调节阀的体积流量,γ1为介质比重,T0为空气温度,p1为进气压力调节阀阀前压力,p2为进气压力调节阀阀后压力,qma为通过进气压力调节阀的质量流量,ρa为介质密度,Ka为进气压力调节阀的阀门开度,Aa为进气压力调节阀的阀门额定面积,cp1为介质定压比热,R为气体常数,Ta为介质温度,△p为进气压力调节阀阀后压力与阀前压力的压差;
其中,进气压力调节阀的阀门额定面积,介质密度,介质定压比热,介质比重为常数;
2)、将测量所得的介质温度,进气压力调节阀阀前压力,进气压力调节阀阀后压力与阀前压力的压差作为参变量,将进气压力调节阀的阀门开度作为自变量;
3)以气源设计压力的85%为上限,以所需最小的进气压力调节阀阀前压力为下限,将由所需最小的进气压力调节阀阀前压力和气源设计压力的85%之间的压力分割成若干压力点,在每个压力点下,将进气压力调节阀按10%~90%、间隔10%的开度正反行程开启三遍,记录每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量,将上述记录的每个压力下每个开度的进气压力调节阀阀前压力、进气压力调节阀阀后压力与阀前压力的压差、以及通过进气压力调节阀的质量流量的3个数据取平均数,拟合进气压力调节阀流量-开度特性曲线。
4.根据权利要求3所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,确定进气压力调节阀临界点的方法,其步骤为:在进气压力调节阀阀前压力固定和进气压力调节阀的阀门开度变化过程中,在某个进气压力调节阀的阀门开度下出现临界流动时,然后逐步减小进气压力调节阀的阀门开度,直至通过进气压力调节阀的质量流量与进气压力调节阀阀前压力不成正比关系,由此点的进气压力调节阀阀后压力即可确定进气压力调节阀临界点。
5.根据权利要求1所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,发动机燃料比例节流阀的调试方法的调试内容包括确定所述发动机燃料比例节流阀的流量-开度曲线、分辨率和重复性精度;所述发动机燃料比例节流阀的分辨率为使流量发生有效、最小变化的阀门开度变化值。
6.根据权利要求5所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,确定发动机燃料比例节流阀流量-开度特性曲线的方法,包括如下步骤:
1)、通过发动机燃料比例节流阀的质量流量仅与发动机燃料比例节流阀阀后压力与阀前压力的压差、工作介质密度、阀门额定面积、阀门开度有关,由公式
<mrow>
<msub>
<mi>q</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>V</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mi>b</mi>
</msub>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>&Delta;p</mi>
<mi>j</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
</mrow>
<mi>g</mi>
</mfrac>
</msqrt>
</mrow>
确定发动机燃料比例节流阀的质量流量与发动机燃料比例节流阀的阀门开度的关系;
上述公式中,qm为通过发动机燃料比例节流阀的质量流量,Kb为发动机燃料比例节流阀的阀门开度,CVj为发动机燃料比例节流阀的阀门流量系数,Ab为发动机燃料比例节流阀阀门额定面积,△pj为发动机燃料比例节流阀阀后压力与阀前压力的压差,ρb为工作介质密度,g为重力加速度,其为已知常数;
2)、将发动机燃料比例节流阀的阀门额定面积、工作介质密度作为常数,将发动机燃料比例节流阀阀后压力与阀前压力的压差作为参变量,将发动机燃料比例节流阀的阀门开度作为自变量;
3)、在试车台进气***的任一工作压力下,均会产生一个发动机燃料比例节流阀阀后压力与阀前压力的压差,对于任一个发动机燃料比例节流阀阀后压力与阀前压力的压差,将发动机燃料比例节流阀的阀门开度从5%、10%、15%、…、90%、95%开启三遍,进行喷液试验,记录通过发动机燃料比例节流阀的质量流量,将每个阀门开度所记录的三个发动机燃料比例节流阀的质量流量取平均值,作为有效的发动机燃料比例节流阀的质量流量,进而计算出发动机燃料比例节流阀流量系数,同时拟合发动机燃料比例节流阀的流量-开度曲线。
7.根据权利要求1所述的采用预测控制法调节试车台进气状态参数的调试方法,其特征在于,所述预测控制法用于直连式试车台。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511017540.1A CN105628387B (zh) | 2015-12-30 | 2015-12-30 | 采用预测控制法调节试车台进气状态参数的调试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511017540.1A CN105628387B (zh) | 2015-12-30 | 2015-12-30 | 采用预测控制法调节试车台进气状态参数的调试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105628387A CN105628387A (zh) | 2016-06-01 |
CN105628387B true CN105628387B (zh) | 2018-01-05 |
Family
ID=56043526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511017540.1A Active CN105628387B (zh) | 2015-12-30 | 2015-12-30 | 采用预测控制法调节试车台进气状态参数的调试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105628387B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106555683B (zh) * | 2016-11-18 | 2019-05-31 | 龙岩学院 | 一种沼气发动机点火控制方法 |
CN109115506B (zh) * | 2018-10-05 | 2020-06-19 | 北京航天三发高科技有限公司 | 一种发动机试车台进气***的使用方法 |
CN109141894A (zh) * | 2018-10-05 | 2019-01-04 | 北京航天三发高科技有限公司 | 一种发动机试车台进气*** |
CN110595790A (zh) * | 2019-08-20 | 2019-12-20 | 南京理工大学 | 加温器补氧流量调节整流装置 |
CN112611567B (zh) * | 2020-11-27 | 2022-05-27 | 中国航发四川燃气涡轮研究院 | 一种航空发动机高空模拟试车台调节阀自动选型方法 |
CN112764437A (zh) * | 2020-12-29 | 2021-05-07 | 北京动力机械研究所 | 基于气动调节阀的试车台空气流量调节的自动控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6131447A (en) * | 1997-08-14 | 2000-10-17 | Bayerische Motoren Werke Aktiengesellschaft | Process and arrangement for checking the operability of a secondary air pump in an internal-combustion engine |
CN101479143A (zh) * | 2007-03-01 | 2009-07-08 | 桂林吉星电子等平衡动力有限公司 | 发动机伺服加裁装置及其动态寻优运行控制方法 |
CN103291464A (zh) * | 2012-02-28 | 2013-09-11 | 通用电气公司 | 基于传感器的追求性能的燃气涡轮发动机控制 |
CN103499448A (zh) * | 2013-10-22 | 2014-01-08 | 北京汽车股份有限公司 | 一种发动机热平衡性能台架测试装置及方法 |
JP2014224786A (ja) * | 2013-05-17 | 2014-12-04 | 株式会社小野測器 | エンジンモデル作成方法およびエンジンモデル作成システム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6460409B1 (en) * | 2000-05-13 | 2002-10-08 | Ford Global Technologies, Inc. | Feed-forward observer-based control for estimating cylinder air charge |
-
2015
- 2015-12-30 CN CN201511017540.1A patent/CN105628387B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6131447A (en) * | 1997-08-14 | 2000-10-17 | Bayerische Motoren Werke Aktiengesellschaft | Process and arrangement for checking the operability of a secondary air pump in an internal-combustion engine |
CN101479143A (zh) * | 2007-03-01 | 2009-07-08 | 桂林吉星电子等平衡动力有限公司 | 发动机伺服加裁装置及其动态寻优运行控制方法 |
CN103291464A (zh) * | 2012-02-28 | 2013-09-11 | 通用电气公司 | 基于传感器的追求性能的燃气涡轮发动机控制 |
JP2014224786A (ja) * | 2013-05-17 | 2014-12-04 | 株式会社小野測器 | エンジンモデル作成方法およびエンジンモデル作成システム |
CN103499448A (zh) * | 2013-10-22 | 2014-01-08 | 北京汽车股份有限公司 | 一种发动机热平衡性能台架测试装置及方法 |
Non-Patent Citations (1)
Title |
---|
某弹用冲压发动机试车台自动测控***设计;牛明;《中国优秀硕士学位论文全文数据库信息科技辑》;20130215;第7-8,41,50-52页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105628387A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105628387B (zh) | 采用预测控制法调节试车台进气状态参数的调试方法 | |
CN105486511B (zh) | 一种试车台状态参数的调试方法 | |
CN105445007B (zh) | 一种燃气轮机用气体燃料控制阀流量特性试验***及方法 | |
CN105628389B (zh) | 一种试车台状态点模拟控制方法 | |
US8700291B2 (en) | System for controlling the exhaust gas return rate by means of virtual NOx sensors with adaptation via a NOx sensor | |
US7281440B2 (en) | Particulate sampling system having flow check device | |
CN111106366B (zh) | 一种燃料电池电堆测试台及其背压控制方法 | |
CN102562329B (zh) | 具有内、外废气再循环的氮氧化物控制 | |
RU2012127271A (ru) | Способ управления двигателем (варианты) и система двигателя | |
CN107269407B (zh) | 用于确定燃烧马达中的新鲜空气质量流量的方法和装置 | |
CN104963780B (zh) | 用于机动车的内燃机,尤其燃气马达 | |
CN203310632U (zh) | 一种发动机增压***模拟装置 | |
CN107063697A (zh) | 一种空气加热***及燃烧室试验台*** | |
EP3040543A1 (en) | Method to control the temperature of the exhaust gases of a supercharged internal combustion engine | |
CN210239874U (zh) | 一种实现柴油单缸机egr的*** | |
CN102434333B (zh) | 柴油机台架试验用辅助增压***及其使用方法 | |
CN105651519B (zh) | 一种试车台进气***进气压力调节阀的调试方法 | |
CN102901632B (zh) | 变缸径变行程中速单缸柴油机试验台的增压***模拟装置 | |
JP5924743B2 (ja) | 内燃エンジンのシリンダの入口における再循環排気ガスの比率を決定する方法、及びこのような方法を実行するエンジン | |
CN105651517B (zh) | 一种并联的试车台进气***压力调节阀的使用方法 | |
CN105628388B (zh) | 一种试车台状态点进气模拟控制方法 | |
CN105547703B (zh) | 一种试车台状态点进气模拟调节方法 | |
CN105628391B (zh) | 一种试车台进气***串联的压力调节阀的使用方法 | |
CN212155023U (zh) | 一种试验平台 | |
CN209027795U (zh) | 一种发动机试车台进气*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |