CN105628387B - Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter - Google Patents

Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter Download PDF

Info

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
Application number
CN201511017540.1A
Other languages
Chinese (zh)
Other versions
CN105628387A (en
Inventor
周培好
徐元元
陆小平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Aerospace Sanfa High Tech Co Ltd
Original Assignee
Beijing Aerospace Sanfa High Tech Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Aerospace Sanfa High Tech Co Ltd filed Critical Beijing Aerospace Sanfa High Tech Co Ltd
Priority to CN201511017540.1A priority Critical patent/CN105628387B/en
Publication of CN105628387A publication Critical patent/CN105628387A/en
Application granted granted Critical
Publication of CN105628387B publication Critical patent/CN105628387B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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)
  • Feedback Control In General (AREA)
  • Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)

Abstract

The invention discloses a kind of adjustment method that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, the type of the debugging includes unloaded debugging, with debugging is carried, and the band carries debugging, comprised the following steps：The controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, inlet pressure adjusting valve, heater fuel ratio adjusting valve, engine fuel proportional throttle valve, liquid oxygen flow control valve are debugged successively；Compared with prior art, the adjustment method provided by the invention that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, the performance of test bay can be integrally grasped, reduces experimentation cost, improves test efficiency.

Description

Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter
Technical field
The present invention relates to a kind of adjustment method that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method.
Background technology
At present, during engine test, it is desirable to one group of inlet condition of state point, including inlet total pressure, Induction air flow, Air inlet oxygen composition, inlet total temperature, engine oil oil mass, at the same reach setting value and it is stable after, record its performance and parameter or Examine its performance.But, it is necessary to be debugged to test bay, to determine that the performance of test bay whether can before engine test Enough meet the test requirements document of engine, in test bay debugging, the debugging of gas handling system is important component, still, for Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter, at present still without the adjustment method of maturation, this leads Cause integrally to grasp the performance of test bay, so as to add simulated time, improve experimentation cost.
The content of the invention
The defects of the invention aims to overcome above-mentioned prior art to exist and provide and a kind of use PREDICTIVE CONTROL Method adjusts the adjustment method of test bay air-suction state parameter.
The purpose of the present invention can be achieved through the following technical solutions：
A kind of adjustment method that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, the type of the debugging are included Zero load debugging, band carry debugging, and the band carries debugging, comprised the following steps：
1), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, inlet pressure adjusting valve is debugged；
2), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, heater fuel ratio adjusting valve is debugged
3), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, engine fuel proportional throttle valve is debugged
4), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, liquid oxygen flow control valve is debugged.
The PREDICTIVE CONTROL method, comprises the following steps：
1), according to the inlet total pressure of current state point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil Oil mass, the inlet total pressure with required emulation mode point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil oil mass The difference of setting value the controlled quentity controlled variable of single parameter is calculated by gas handling system status point parameter solution of equations；
2), according to the controlled quentity controlled variable of the single parameter calculated in the step 1), each valve of adjustment control above-mentioned parameter The aperture of door so that when subsequent cycle starts, set each valve aperture meet to be calculated in the step 1) it is single The controlled quentity controlled variable of parameter；
3) repeating said steps 1) and step 2) until simulation precision meet require；
The gas handling system status point parameter equation group is：
qmi=qma+qmo+qmj
pti=pt2+Δpt2
In the gas handling system status point parameter equation group, qmiFor charge flow rate, qmaTo pass through inlet pressure adjusting valve Mass flow, qmoFor oxygen supplying flow, qmjFor heater fuel supply flow rate, ptiFor inlet total pressure, AiFor intake area, TtiFor Inlet total temperature, q (λi) it is to calculate section Ventialtion rate, Tti0For current inlet total temperature, η is efficiency of combustion, QjFor heater The fuel feeding latent heat of vaporization, α are overfire air coefficient, CpFor combustion gas specific heat at constant pressure, CvFor the valve flow system of inlet pressure adjusting valve Number, PsaFor bleed pressure, pt2For air control valve downstream pressure, γ is air quality proportion, T0aFor air source temperature, Δ Pt2 For duct pressure loss, TsaFor air themperature before mixing, CpaFor pressurization by compressed air specific heat, TsoFor oxygenating temperature, C before mixingpoFor oxygen Gas specific heat at constant pressure, Q0For the heater fuel feeding latent heat of vaporization,For chemically correct fuel；
In the gas handling system status point parameter equation group, it is known that the Induction air flow of current state point, current state point Oxygen supplying flow, the heater fuel supply flow rate of current state point, the inlet total pressure of current state point, the inlet total temperature of current state point, when The current inlet total temperature of preceding state point, the overfire air coefficient of current state point, and the inlet total pressure of required emulation mode point, required mould Intend the inlet total temperature of state point, solve and the inlet total pressure of required emulation mode point and the inlet total temperature pair of required emulation mode point The Induction air flow for the required emulation mode point answered, the oxygen supplying flow of required emulation mode point, required emulation mode point add Hot device fuel supply flow rate.
The commissioning content of the inlet pressure adjusting valve includes determining that flow-open degree characteristic of inlet pressure adjusting valve is bent Line, critical point, resolution ratio, repeatability precision, the checking of load pressure loading time, loaded pressure drop curve；The admission pressure is adjusted The resolution ratio of section valve is to make flow or admission pressure that effective, the controlling opening of valve changing value of minimum change, the load occur Pressure drop is inlet pressure adjusting valve downstream pressure and the difference of intake simulation stagnation pressure.
The method for determining inlet pressure adjusting valve flow-open degree characteristic curve, comprises the following steps：
1) mass flow and upstream pressure and downstream pressure of inlet pressure adjusting valve, Media density, temperature, valve, are passed through Door nominal area, the relation of valve opening, by formula
When
When
With qma=Q ρaObtain：
In above-mentioned formula, CVFor the valve flow coefficient of inlet pressure adjusting valve, Q is the body by inlet pressure adjusting valve Product flow, γ1For medium proportion, T0For air themperature, p1For inlet pressure adjusting valve upstream pressure, p2Adjusted for admission pressure Valve downstream pressure, qma are by the mass flow of inlet pressure adjusting valve, ρaFor Media density, KaFor inlet pressure adjusting valve Valve opening, AaFor the valve nominal area of inlet pressure adjusting valve, cp1For medium specific heat at constant pressure, R is gas constant, TaTo be situated between Matter temperature, Δ p are the pressure difference of inlet pressure adjusting valve downstream pressure and upstream pressure；
Wherein, the valve nominal area of inlet pressure adjusting valve, Media density, medium specific heat at constant pressure, medium proportion are normal Number；
2) medium temperature of gained, inlet pressure adjusting valve upstream pressure, inlet pressure adjusting valve downstream pressure will, be measured Pressure difference with upstream pressure is as parameter, using the valve opening of inlet pressure adjusting valve as independent variable；
3) it is the upper limit with the 85% of design of gas source pressure, using required minimum inlet pressure adjusting valve upstream pressure under Limit, if will be divided into by the pressure between required minimum inlet pressure adjusting valve upstream pressure and the 85% of source of the gas design pressure Dry-pressing force, under each pressure spot, inlet pressure adjusting valve is opened by the positive revesal of aperture of 10%~90%, interval 10% Open three times, record the inlet pressure adjusting valve upstream pressure of each aperture, inlet pressure adjusting valve downstream pressure under each pressure Pressure difference with upstream pressure and the mass flow by inlet pressure adjusting valve, will be each under each pressure of above-mentioned record The inlet pressure adjusting valve upstream pressure of aperture, the pressure difference of inlet pressure adjusting valve downstream pressure and upstream pressure and pass through 3 data of the mass flow of inlet pressure adjusting valve are taken the mean, and fitting inlet pressure adjusting valve flow-open degree characteristic is bent Line.
The method for determining inlet pressure adjusting valve critical point, its step are：Fixed in inlet pressure adjusting valve upstream pressure In the valve opening change procedure of inlet pressure adjusting valve, occur under the valve opening of some inlet pressure adjusting valve critical During flowing, then progressively reduce the valve opening of inlet pressure adjusting valve, until the mass flow for passing through inlet pressure adjusting valve It is not proportional with inlet pressure adjusting valve upstream pressure, the inlet pressure adjusting valve downstream pressure thus put can determine that into Air pressure force regulating valve critical point.
The commissioning content of the adjustment method of engine fuel proportional throttle valve includes determining the engine fuel ratio section Flow flow-opening curve, resolution ratio and the repeatability precision of valve；The resolution ratio of the engine fuel proportional throttle valve is to make stream Effective, the valve opening changing value of minimum change occurs for amount.
The method for determining engine fuel proportional throttle valve flow-open degree characteristic curve, comprises the following steps：
1), by the mass flow of engine fuel proportional throttle valve only with engine fuel proportional throttle valve downstream pressure Pressure difference, working media density with upstream pressure, valve nominal area, valve opening are relevant, by formula
Determine the mass flow of engine fuel proportional throttle valve and the valve opening of engine fuel proportional throttle valve Relation；
In above-mentioned formula, qm is by the mass flow of engine fuel proportional throttle valve, KbFor engine fuel ratio The valve opening of choke valve, CVjFor the valve flow coefficient of engine fuel proportional throttle valve, AbFor engine fuel ratio section Flow valve valve nominal area, Δ piFor engine fuel proportional throttle valve downstream pressure and the pressure difference of upstream pressure, ρbIt is situated between for work Matter density, g are acceleration of gravity, and it is known constant；
2), using the valve nominal area of engine fuel proportional throttle valve, working media density as constant, by engine The pressure difference of fuel ratio choke valve downstream pressure and upstream pressure is as parameter, by the valve of engine fuel proportional throttle valve Aperture is as independent variable；
3), under any operative pressure of test bay gas handling system, an engine fuel proportional throttle valve can be produced The pressure difference of downstream pressure and upstream pressure, for any one engine fuel proportional throttle valve downstream pressure and the pressure of upstream pressure Difference, by the valve opening of engine fuel proportional throttle valve from 5%, 10%, 15% ..., 90%, 95% open three times, carry out Hydrojet is tested, and record is by the mass flow of engine fuel proportional throttle valve, three hairs that each valve opening is recorded The mass flow of motivation fuel ratio choke valve is averaged, the quality stream as effective engine fuel proportional throttle valve Amount, and then calculates engine fuel proportional throttle valve discharge coefficient, at the same be fitted the flow of engine fuel proportional throttle valve- Opening curve.
The zero load is debugged and can replaced with debugging is carried with simulation cylinder.
The PREDICTIVE CONTROL method is used for direct-connected test bay.
The PREDICTIVE CONTROL method is used for engine condition point intake simulation available for direct-connected test bay.
The invention has the advantages that compared with prior art, it is provided by the invention using the regulation test run of PREDICTIVE CONTROL method The adjustment method of platform air-suction state parameter, the performance of test bay can be integrally grasped, reduce experimentation cost, improve test efficiency.
Embodiment
A kind of adjustment method that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, the type of the debugging are included Zero load debugging, band carry debugging, and the band carries debugging, comprised the following steps：
1), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, inlet pressure adjusting valve is debugged；
2), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, heater fuel ratio adjusting valve is debugged
3), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, engine fuel proportional throttle valve is debugged
4), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, liquid oxygen flow control valve is debugged.
The PREDICTIVE CONTROL method, comprises the following steps：
1), according to the inlet total pressure of current state point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil Oil mass, the inlet total pressure with required emulation mode point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil oil mass The difference of setting value the controlled quentity controlled variable of single parameter is calculated by gas handling system status point parameter solution of equations；
2), according to the controlled quentity controlled variable of the single parameter calculated in the step 1), each valve of adjustment control above-mentioned parameter The aperture of door so that when subsequent cycle starts, set each valve aperture meet to be calculated in the step 1) it is single The controlled quentity controlled variable of parameter；
3) repeating said steps 1) and step 2) until simulation precision meet require；
The gas handling system status point parameter equation group is：
qmi=qma+qmo+qmj
pti=pt2+Δpt2
In the gas handling system status point parameter equation group, qmiFor charge flow rate, qmaTo pass through inlet pressure adjusting valve Mass flow, qmoFor oxygen supplying flow, qmjFor heater fuel supply flow rate, ptiFor inlet total pressure, AiFor intake area, TtiFor Inlet total temperature, q (λi) it is to calculate section Ventialtion rate, Tti0For current inlet total temperature, η is efficiency of combustion, QjFor heater The fuel feeding latent heat of vaporization, α are overfire air coefficient, CpFor combustion gas specific heat at constant pressure, CvFor the valve flow system of inlet pressure adjusting valve Number, PsaFor bleed pressure, pt2For air control valve downstream pressure, γ is air quality proportion, T0aFor air source temperature, Δ Pt2 For duct pressure loss, TsaFor air themperature before mixing, CpaFor pressurization by compressed air specific heat, TsoFor oxygenating temperature, C before mixingpoFor oxygen Gas specific heat at constant pressure, Q0For the heater fuel feeding latent heat of vaporization,For chemically correct fuel；
In the gas handling system status point parameter equation group, it is known that the Induction air flow of current state point, current state point Oxygen supplying flow, the heater fuel supply flow rate of current state point, the inlet total pressure of current state point, the inlet total temperature of current state point, when The current inlet total temperature of preceding state point, the overfire air coefficient of current state point, and the inlet total pressure of required emulation mode point, required mould Intend the inlet total temperature of state point, solve and the inlet total pressure of required emulation mode point and the inlet total temperature pair of required emulation mode point The Induction air flow for the required emulation mode point answered, the oxygen supplying flow of required emulation mode point, required emulation mode point add Hot device fuel supply flow rate.
The commissioning content of the inlet pressure adjusting valve includes determining that flow-open degree characteristic of inlet pressure adjusting valve is bent Line, critical point, resolution ratio, repeatability precision, the checking of load pressure loading time, loaded pressure drop curve；The admission pressure is adjusted The resolution ratio of section valve is to make flow or admission pressure that effective, the controlling opening of valve changing value of minimum change, the load occur Pressure drop is inlet pressure adjusting valve downstream pressure and the difference of intake simulation stagnation pressure.
The method for determining inlet pressure adjusting valve flow-open degree characteristic curve, comprises the following steps：
1) mass flow and upstream pressure and downstream pressure of inlet pressure adjusting valve, Media density, temperature, valve, are passed through Door nominal area, the relation of valve opening, by formula
When
When
With qma=Q ρaObtain：
In above-mentioned formula, CVFor the valve flow coefficient of inlet pressure adjusting valve, Q is the body by inlet pressure adjusting valve Product flow, γ1For medium proportion, T0For air themperature, p1For inlet pressure adjusting valve upstream pressure, p2Adjusted for admission pressure Valve downstream pressure, qma are by the mass flow of inlet pressure adjusting valve, ρaFor Media density, KaFor inlet pressure adjusting valve Valve opening, AaFor the valve nominal area of inlet pressure adjusting valve, cp1For medium specific heat at constant pressure, R is gas constant, TaTo be situated between Matter temperature, Δ p are the pressure difference of inlet pressure adjusting valve downstream pressure and upstream pressure；
Wherein, the valve nominal area of inlet pressure adjusting valve, Media density, medium specific heat at constant pressure, medium proportion are normal Number；
2) medium temperature of gained, inlet pressure adjusting valve upstream pressure, inlet pressure adjusting valve downstream pressure will, be measured Pressure difference with upstream pressure is as parameter, using the valve opening of inlet pressure adjusting valve as independent variable；
3) it is the upper limit with the 85% of design of gas source pressure, using required minimum inlet pressure adjusting valve upstream pressure under Limit, if will be divided into by the pressure between required minimum inlet pressure adjusting valve upstream pressure and the 85% of source of the gas design pressure Dry-pressing force, under each pressure spot, inlet pressure adjusting valve is opened by the positive revesal of aperture of 10%~90%, interval 10% Open three times, record the inlet pressure adjusting valve upstream pressure of each aperture, inlet pressure adjusting valve downstream pressure under each pressure Pressure difference with upstream pressure and the mass flow by inlet pressure adjusting valve, will be each under each pressure of above-mentioned record The inlet pressure adjusting valve upstream pressure of aperture, the pressure difference of inlet pressure adjusting valve downstream pressure and upstream pressure and pass through 3 data of the mass flow of inlet pressure adjusting valve are taken the mean, and fitting inlet pressure adjusting valve flow-open degree characteristic is bent Line.
The method for determining inlet pressure adjusting valve critical point, its step are：Fixed in inlet pressure adjusting valve upstream pressure In the valve opening change procedure of inlet pressure adjusting valve, occur under the valve opening of some inlet pressure adjusting valve critical During flowing, then progressively reduce the valve opening of inlet pressure adjusting valve, until the mass flow for passing through inlet pressure adjusting valve It is not proportional with inlet pressure adjusting valve upstream pressure, the inlet pressure adjusting valve downstream pressure thus put can determine that into Air pressure force regulating valve critical point.
The commissioning content of the adjustment method of engine fuel proportional throttle valve includes determining the engine fuel ratio section Flow flow-opening curve, resolution ratio and the repeatability precision of valve；The resolution ratio of the engine fuel proportional throttle valve is to make stream Effective, the valve opening changing value of minimum change occurs for amount.
The method for determining engine fuel proportional throttle valve flow-open degree characteristic curve, comprises the following steps：
1), by the mass flow of engine fuel proportional throttle valve only with engine fuel proportional throttle valve downstream pressure Pressure difference, working media density with upstream pressure, valve nominal area, valve opening are relevant, by formula
Determine the mass flow of engine fuel proportional throttle valve and the valve opening of engine fuel proportional throttle valve Relation；
In above-mentioned formula, qm is by the mass flow of engine fuel proportional throttle valve, KbFor engine fuel ratio The valve opening of choke valve, CVjFor the valve flow coefficient of engine fuel proportional throttle valve, AbFor engine fuel ratio section Flow valve valve nominal area, Δ pjFor engine fuel proportional throttle valve downstream pressure and the pressure difference of upstream pressure, ρbIt is situated between for work Matter density, g are acceleration of gravity, and it is known constant；
2), using the valve nominal area of engine fuel proportional throttle valve, working media density as constant, by engine The pressure difference of fuel ratio choke valve downstream pressure and upstream pressure is as parameter, by the valve of engine fuel proportional throttle valve Aperture is as independent variable；
3), under any operative pressure of test bay gas handling system, an engine fuel proportional throttle valve can be produced The pressure difference of downstream pressure and upstream pressure, for any one engine fuel proportional throttle valve downstream pressure and the pressure of upstream pressure Difference, by the valve opening of engine fuel proportional throttle valve from 5%, 10%, 15% ..., 90%, 95% open three times, carry out Hydrojet is tested, and record is by the mass flow of engine fuel proportional throttle valve, three hairs that each valve opening is recorded The mass flow of motivation fuel ratio choke valve is averaged, the quality stream as effective engine fuel proportional throttle valve Amount, and then calculates engine fuel proportional throttle valve discharge coefficient, at the same be fitted the flow of engine fuel proportional throttle valve- Opening curve.
The zero load is debugged and can replaced with debugging is carried with simulation cylinder.
The PREDICTIVE CONTROL method is used for direct-connected test bay.
The PREDICTIVE CONTROL method is used for engine condition point intake simulation available for direct-connected test bay.
Finally illustrate, above example is merely to illustrate technical scheme and unrestricted, although with reference to compared with The present invention is described in detail good embodiment, it will be understood by those of skill in the art that can be to the technology of the present invention Scheme is modified and equivalent substitute, and without departing from the objective and scope of technical solution of the present invention, it all should cover in the present invention Right among.

Claims (7)

  1. A kind of 1. adjustment method that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is characterised in that the debugging Type include unloaded debugging, band carries debugging, the band carries debugging, comprised the following steps：
    1), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, inlet pressure adjusting valve is debugged；
    2), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, heater fuel ratio adjusting valve is debugged；
    3), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, engine fuel proportional throttle valve is debugged；
    4), the controlled quentity controlled variable of the single parameter drawn by PREDICTIVE CONTROL method, liquid oxygen flow control valve is debugged；
    The PREDICTIVE CONTROL method, comprises the following steps：
    A), according to the inlet total pressure of current state point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil oil Amount, the inlet total pressure with required emulation mode point, Induction air flow, oxygen supplying flow, inlet total temperature, engine oil oil mass The difference of setting value calculates the controlled quentity controlled variable of single parameter by gas handling system status point parameter solution of equations；
    B), according to the controlled quentity controlled variable of the single parameter calculated in the step a), each valve of adjustment control above-mentioned parameter Aperture so that when subsequent cycle starts, the aperture for setting each valve meets the single parameter calculated in the step a) Controlled quentity controlled variable；
    C) repeating said steps a) and step b) is until simulation precision meets to require；
    The gas handling system status point parameter equation group is：
    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>&amp;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>&amp;eta;q</mi> <mrow> <mi>m</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>Q</mi> <mi>j</mi> </msub> </mrow> <mrow> <msub> <mi>&amp;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>&amp;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>&amp;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>&amp;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>&amp;alpha;</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>q</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow>
    In the gas handling system status point parameter equation group, qmiFor Induction air flow, qmaTo pass through inlet pressure adjusting valve Mass flow, qmoFor oxygen supplying flow, qmjFor heater fuel supply flow rate, ptiFor inlet total pressure, AiFor intake area, TtiFor Inlet total temperature, q (λi) it is to calculate section Ventialtion rate, Tti0For current inlet total temperature, η is efficiency of combustion, QjFor heater The fuel feeding latent heat of vaporization, α are overfire air coefficient, CpFor combustion gas specific heat at constant pressure, CvFor the valve flow system of inlet pressure adjusting valve Number, PsaFor bleed pressure, pt2For air control valve downstream pressure, γ is air quality proportion, T0aFor air source temperature, Δ Pt2 For duct pressure loss, TsaFor air themperature before mixing, CpaFor pressurization by compressed air specific heat, TsoFor oxygenating temperature, C before mixingpoFor oxygen Gas specific heat at constant pressure, Q0For the heater fuel feeding latent heat of vaporization,For chemically correct fuel；
    In the gas handling system status point parameter equation group, it is known that the Induction air flow of current state point, the benefit of current state point Oxygen flow, the heater fuel supply flow rate of current state point, the inlet total pressure of current state point, the inlet total temperature of current state point, current state The current inlet total temperature of point, the overfire air coefficient of current state point, and the inlet total pressure of required emulation mode point, required simulation shape The inlet total temperature of state point, solve corresponding with the inlet total pressure of required emulation mode point and the inlet total temperature of required emulation mode point The Induction air flow of required emulation mode point, the oxygen supplying flow of required emulation mode point, the heater of required emulation mode point Fuel supply flow rate.
  2. 2. the adjustment method according to claim 1 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is, the commissioning content of the inlet pressure adjusting valve include determining flow-open degree characteristic curve of inlet pressure adjusting valve, Critical point, resolution ratio, repeatability precision, the checking of load pressure loading time and loaded pressure drop curve；The admission pressure regulation The resolution ratio of valve is to make flow or admission pressure that effective, the controlling opening of valve changing value of minimum change, the load pressure occur It is reduced to the difference of inlet pressure adjusting valve downstream pressure and intake simulation stagnation pressure.
  3. 3. the adjustment method according to claim 2 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is, the method for determining inlet pressure adjusting valve flow-open degree characteristic curve, comprises the following steps：
    1) mass flow and upstream pressure and downstream pressure of inlet pressure adjusting valve, Media density, temperature, valve volume, are passed through Area, the relation of valve opening are determined, by formula
    (p1-p2) ＜ p1When/2
    (p1-p2)≥p1When/2
    And qma=Q ρaObtain：
    <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>&amp;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>&amp;gamma;</mi> <mn>1</mn> </msub> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>p</mi> </mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>2</mn> <msub> <mi>RT</mi> <mi>a</mi> </msub> </mrow> </mfrac> </msqrt> <mrow> <mo>(</mo> <mi>&amp;Delta;</mi> <mi>p</mi> <mo>&lt;</mo> <msub> <mi>&amp;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>&amp;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>&amp;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>&amp;Delta;</mi> <mi>p</mi> <mo>&amp;GreaterEqual;</mo> <msub> <mi>&amp;beta;</mi> <mi>p</mi> </msub> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    In above-mentioned formula, CVFor the valve flow coefficient of inlet pressure adjusting valve, Q is the volume flow by inlet pressure adjusting valve Amount, γ1For medium proportion, T0For air themperature, p1For inlet pressure adjusting valve upstream pressure, p2For inlet pressure adjusting valve valve Pressure afterwards, qmaTo pass through the mass flow of inlet pressure adjusting valve, ρaFor Media density, KaFor the valve of inlet pressure adjusting valve Aperture, AaFor the valve nominal area of inlet pressure adjusting valve, cp1For medium specific heat at constant pressure, R is gas constant, TaFor medium temperature Degree, △ p are the pressure difference of inlet pressure adjusting valve downstream pressure and upstream pressure；
    Wherein, the valve nominal area of inlet pressure adjusting valve, Media density, medium specific heat at constant pressure, medium proportion are constant；
    2) medium temperature of gained, inlet pressure adjusting valve upstream pressure, inlet pressure adjusting valve downstream pressure and valve will, be measured The pressure difference of preceding pressure is as parameter, using the valve opening of inlet pressure adjusting valve as independent variable；
    3) it is the upper limit with the 85% of design of gas source pressure, will using required minimum inlet pressure adjusting valve upstream pressure as lower limit If dry-pressing is divided into by the pressure between required minimum inlet pressure adjusting valve upstream pressure and the 85% of source of the gas design pressure Force, under each pressure spot, inlet pressure adjusting valve is opened three by the positive revesal of aperture of 10%~90%, interval 10% Time, record inlet pressure adjusting valve upstream pressure, inlet pressure adjusting valve downstream pressure and the valve of each aperture under each pressure The pressure difference of preceding pressure and the mass flow by inlet pressure adjusting valve, by each aperture under each pressure of above-mentioned record Inlet pressure adjusting valve upstream pressure, inlet pressure adjusting valve downstream pressure and upstream pressure pressure difference and pass through air inlet 3 data of the mass flow of pressure-regulating valve are taken the mean, fitting inlet pressure adjusting valve flow-open degree characteristic curve.
  4. 4. the adjustment method according to claim 3 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is, the method for determining inlet pressure adjusting valve critical point, and its step is：Inlet pressure adjusting valve upstream pressure fix and In the valve opening change procedure of inlet pressure adjusting valve, there is critical flow under the valve opening of some inlet pressure adjusting valve When dynamic, then progressively reduce the valve opening of inlet pressure adjusting valve, until by the mass flow of inlet pressure adjusting valve with Inlet pressure adjusting valve upstream pressure is not proportional, and the inlet pressure adjusting valve downstream pressure thus put can determine that air inlet Pressure-regulating valve critical point.
  5. 5. the adjustment method according to claim 1 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is that the commissioning content of the adjustment method of engine fuel proportional throttle valve includes determining the engine fuel ratio throttling Flow-opening curve, resolution ratio and the repeatability precision of valve；The resolution ratio of the engine fuel proportional throttle valve is to make flow Effective, the valve opening changing value of minimum change occurs.
  6. 6. the adjustment method according to claim 5 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is, the method for determining engine fuel proportional throttle valve flow-open degree characteristic curve, comprises the following steps：
    1), by the mass flow of engine fuel proportional throttle valve only with engine fuel proportional throttle valve downstream pressure and valve The pressure difference of preceding pressure, working media density, valve nominal area, valve opening are relevant, by formula
    <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>&amp;Delta;p</mi> <mi>j</mi> </msub> <msub> <mi>&amp;rho;</mi> <mi>b</mi> </msub> </mrow> <mi>g</mi> </mfrac> </msqrt> </mrow>
    Determine the relation of the mass flow of engine fuel proportional throttle valve and the valve opening of engine fuel proportional throttle valve；
    In above-mentioned formula, qmTo pass through the mass flow of engine fuel proportional throttle valve, KbFor engine fuel proportional throttle valve Valve opening, CVjFor the valve flow coefficient of engine fuel proportional throttle valve, AbFor engine fuel proportional throttle valve valve Door nominal area, △ pjFor engine fuel proportional throttle valve downstream pressure and the pressure difference of upstream pressure, ρbIt is close for working media Degree, g is acceleration of gravity, and it is known constant；
    2), using the valve nominal area of engine fuel proportional throttle valve, working media density as constant, by engine fuel The pressure difference of proportional throttle valve downstream pressure and upstream pressure is as parameter, by the valve opening of engine fuel proportional throttle valve As independent variable；
    3), under any operative pressure of test bay gas handling system, after producing an engine fuel proportional throttle valve valve The pressure difference of pressure and upstream pressure, for any one engine fuel proportional throttle valve downstream pressure and the pressure difference of upstream pressure, By the valve opening of engine fuel proportional throttle valve from 5%, 10%, 15% ..., 90%, 95% open three times, carry out hydrojet Experiment, record is by the mass flow of engine fuel proportional throttle valve, three engines that each valve opening is recorded The mass flow of fuel ratio choke valve is averaged, and as the mass flow of effective engine fuel proportional throttle valve, is entered And engine fuel proportional throttle valve discharge coefficient is calculated, while it is fitted flow-aperture of engine fuel proportional throttle valve Curve.
  7. 7. the adjustment method according to claim 1 that test bay air-suction state parameter is adjusted using PREDICTIVE CONTROL method, it is special Sign is that the PREDICTIVE CONTROL method is used for direct-connected test bay.
CN201511017540.1A 2015-12-30 2015-12-30 Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter Active CN105628387B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511017540.1A CN105628387B (en) 2015-12-30 2015-12-30 Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511017540.1A CN105628387B (en) 2015-12-30 2015-12-30 Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter

Publications (2)

Publication Number Publication Date
CN105628387A CN105628387A (en) 2016-06-01
CN105628387B true CN105628387B (en) 2018-01-05

Family

ID=56043526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511017540.1A Active CN105628387B (en) 2015-12-30 2015-12-30 Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter

Country Status (1)

Country Link
CN (1) CN105628387B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106555683B (en) * 2016-11-18 2019-05-31 龙岩学院 A kind of methane fuelled engine ignition control method
CN109115506B (en) * 2018-10-05 2020-06-19 北京航天三发高科技有限公司 Use method of air inlet system of engine test bed
CN109141894A (en) * 2018-10-05 2019-01-04 北京航天三发高科技有限公司 A kind of engine testsand gas handling system
CN110595790A (en) * 2019-08-20 2019-12-20 南京理工大学 Heater oxygenating flow adjusting and rectifying device
CN112611567B (en) * 2020-11-27 2022-05-27 中国航发四川燃气涡轮研究院 Automatic model selection method for regulating valve of high-altitude simulation test bed of aero-engine
CN112764437A (en) * 2020-12-29 2021-05-07 北京动力机械研究所 Automatic control method for air flow regulation of test bed based on pneumatic regulating valve

Citations (5)

* Cited by examiner, † Cited by third party
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 (en) * 2007-03-01 2009-07-08 桂林吉星电子等平衡动力有限公司 Engine load servo device and dynamic optimize control method thereof
CN103291464A (en) * 2012-02-28 2013-09-11 通用电气公司 Sensor-based performance-seeking gas turbine engine control
CN103499448A (en) * 2013-10-22 2014-01-08 北京汽车股份有限公司 Engine heat balancing performance bed tester and engine heat balancing performance bed testing method
JP2014224786A (en) * 2013-05-17 2014-12-04 株式会社小野測器 Engine model creation method and engine model creation system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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 (en) * 2007-03-01 2009-07-08 桂林吉星电子等平衡动力有限公司 Engine load servo device and dynamic optimize control method thereof
CN103291464A (en) * 2012-02-28 2013-09-11 通用电气公司 Sensor-based performance-seeking gas turbine engine control
JP2014224786A (en) * 2013-05-17 2014-12-04 株式会社小野測器 Engine model creation method and engine model creation system
CN103499448A (en) * 2013-10-22 2014-01-08 北京汽车股份有限公司 Engine heat balancing performance bed tester and engine heat balancing performance bed testing method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
某弹用冲压发动机试车台自动测控***设计;牛明;《中国优秀硕士学位论文全文数据库信息科技辑》;20130215;第7-8，41，50-52页 *

Also Published As

Publication number Publication date
CN105628387A (en) 2016-06-01

Similar Documents

Publication Publication Date Title
CN105628387B (en) Using the adjustment method of PREDICTIVE CONTROL method regulation test bay air-suction state parameter
CN105486511B (en) A kind of adjustment method of test bay state parameter
CN105445007B (en) A kind of gas turbine gaseous fuel control valve Flow characteristic test system and method
CN105628389B (en) A kind of test bay state point analog control method
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 (en) Fuel cell stack test board and back pressure control method thereof
CN102562329B (en) Oxynitride control having inner and external exhaust gas recirculation
RU2012127271A (en) ENGINE CONTROL METHOD (OPTIONS) AND ENGINE SYSTEM
CN107269407B (en) Method and device for determining a fresh air mass flow in a combustion motor
CN104963780B (en) For the internal combustion engine of motor vehicle, especially gas motor
CN203310632U (en) Device for simulating engine boost system
CN107063697A (en) A kind of air heating system and combustor test platform system
EP3040543A1 (en) Method to control the temperature of the exhaust gases of a supercharged internal combustion engine
CN102434333B (en) Auxiliary pressurization system for diesel engine bench test and using method of auxiliary pressurization system
CN105651519B (en) A kind of adjustment method of test bay gas handling system inlet pressure adjusting valve
CN102901632B (en) Pressurization system simulating device for test bed of cylinder diameter-variable stroke-variable intermediate-speed single cylinder diesel
JP5924743B2 (en) Method for determining the ratio of recirculated exhaust gas at the inlet of a cylinder of an internal combustion engine and engine implementing such a method
CN105651517B (en) A kind of application method of test bay gas handling system pressure-regulating valve in parallel
CN105628388B (en) A kind of test bay state point intake simulation control method
CN105547703B (en) A kind of test bay state point intake simulation adjusting method
CN105628391B (en) A kind of application method of the pressure-regulating valve of test bay gas handling system series connection
CN212155023U (en) Test platform
CN209027795U (en) A kind of engine testsand gas handling system
KR102593928B1 (en) Pressure-controlled mass flow system for multipoint gaseous fuel injection

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