CN105628388B - 一种试车台状态点进气模拟控制方法 - Google Patents
一种试车台状态点进气模拟控制方法 Download PDFInfo
- Publication number
- CN105628388B CN105628388B CN201511017589.7A CN201511017589A CN105628388B CN 105628388 B CN105628388 B CN 105628388B CN 201511017589 A CN201511017589 A CN 201511017589A CN 105628388 B CN105628388 B CN 105628388B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- point
- inlet total
- state point
- 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
Abstract
本发明公开了一种试车台状态点进气模拟控制方法，包括如下步骤：1)、根据当前状态点的进气总压，进气空气流量，补氧流量，进气总温，发动机供油油量，与所需模拟状态点的进气总压，进气空气流量，补氧流量，进气总温，发动机供油油量的设定值的差值通过进气系统状态点参数方程组解算出单个参数的控制量；2)、根据所述步骤1)中所解算出的单个参数的控制量，调节控制上述参数的各阀门开度，使得在下一循环开始时，设定各阀门开度满足所述步骤1)中所解算出的单个参数的控制量；3)重复所述步骤1)和步骤2)直至模拟精度满足要求；与现有技术相比，本发明提供的试车台状态点进气模拟控制方法，缩短了发动机试验时间，降低试验成本，提高试验效率。
Description
技术领域
本发明涉及一种试车台状态点进气模拟控制方法，特别是涉及一种直连式试车台状态点进气模拟控制方法。
背景技术
目前，发动机试验时，要求状态点的一组进气参数，包括进气总压，进气空气流量，进气氧成分，进气总温，发动机供油油量，同时达到设定值及稳定后，记录其性能及参数或考核其性能。这组状态点的进气参数一般采用单个参数自动闭环控制方法进行控制，但是，由于状态点的进气参数间有耦合关系，并且气动负载和温度测量包含惯性环节，单个参数自动闭环控制方法会引起参数的反复震荡，造成模拟时间很难控制。
发明内容
本发明的目的是为了克服上述现有技术存在的缺陷而提供了一种试验时间短，试验成本低，试验效率高的试车台状态点进气模拟控制方法。
本发明的目的可以通过以下技术方案来实现：
一种试车台状态点进气模拟控制方法，包括如下步骤：
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为燃气定压比热，Cva为空气调节阀流量系数，Psa为气源压力，pt2为空气调节阀阀后压力，γ为空气质量比重，T0a为空气源温度，ΔPt2为管道压力损失，Tsa为混合前空气温度，Cpa为空气定压比热，Tso为混合前补氧温度，Cpo为氧气定压比热，Q0为加热器供油汽化潜热，为理论空燃比；
在所述进气系统状态点参数方程组中，已知当前态点的进气空气流量，当前态点的补氧流量，当前态点的加热器供油流量，当前态点的进气总压，当前态点的进气总温，当前态点的当前进气总温，当前态点的过热空气系数，和所需模拟状态点的进气总压，所需模拟状态点的进气总温，求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量，所需模拟状态点的补氧流量，所需模拟状态点的加热器供油流量。
所述进气系统状态点参数方程组中，可将式替换为
所述进气系统状态点参数方程组中，
所述进气系统状态点参数方程组中，所述进气总温为燃烧后的温度，所述当前进气总温为进一步燃烧前温度。
所述试车台状态点进气模拟控制方法可用于直连式试车台。
所述试车台状态点进气模拟控制方法可用于发动机状态点进气模拟。
本发明的有益效果是，与现有技术相比，本发明提供的试车台状态点模拟控制方法，缩短了发动机试验时间，降低试验成本，提高试验效率。
具体实施方式
一种试车台状态点进气模拟控制方法，包括如下步骤：
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为燃气定压比热，Cva为空气调节阀流量系数，Psa为气源压力，pt2为空气调节阀阀后压力，γ为空气质量比重，T0a为空气源温度，ΔPt2为管道压力损失，Tsa为混合前空气温度，Cpa为空气定压比热，Tso为混合前补氧温度，Cpo为氧气定压比热，Q0为加热器供油汽化潜热，为理论空燃比；
在所述进气系统状态点参数方程组中，已知当前态点的进气空气流量，当前态点的补氧流量，当前态点的加热器供油流量，当前态点的进气总压，当前态点的进气总温，当前态点的当前进气总温，当前态点的过热空气系数，和所需模拟状态点的进气总压，所需模拟状态点的进气总温，求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量，所需模拟状态点的补氧流量，所需模拟状态点的加热器供油流量。
所述进气系统状态点参数方程组中，可将式替换为
所述进气系统状态点参数方程组中，
所述进气系统状态点参数方程组中，所述进气总温为燃烧后的温度，所述当前进气总温为进一步燃烧前温度。
所述试车台状态点进气模拟控制方法可用于直连式试车台。
所述试车台状态点进气模拟控制方法可用于发动机状态点进气模拟。
最后说明的是，以上实施例仅用于说明本发明的技术方案而非限制，尽管参照较佳的实施例对本发明进行了详细说明，本领域的技术人员应当理解，可以对本发明的技术方案进行修改和等同替代，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种试车台状态点进气模拟控制方法，其特征在于，包括如下步骤：
1)、根据当前状态点的进气总压，进气空气流量，补氧流量，进气总温，发动机供油油量，与所需模拟状态点的进气总压，进气空气流量，补氧流量，进气总温，发动机供油油量的设定值的差值通过进气系统状态点参数方程组解算出单个参数的控制量；
2)、根据所述步骤1)中所解算出的单个参数的控制量，调节控制上述参数的各阀门的开度，使得在下一循环开始时，设定各阀门的开度满足所述步骤1)中所解算出的单个参数的控制量；
3)重复所述步骤1)和步骤2)直至模拟精度满足要求；
所述进气系统状态点参数方程组为：
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>
<mrow>
<mi>V</mi>
<mi>a</mi>
</mrow>
</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为燃气定压比热，Cva为空气调节阀流量系数，Psa为气源压力，pt2为空气调节阀阀后压力，γ为空气质量比重，T0a为空气源温度，ΔPt2为管道压力损失，Tsa为混合前空气温度，Cpa为空气定压比热，Tso为混合前补氧温度，Cpo为氧气定压比热，Q0为加热器供油汽化潜热，为理论空燃比；
在所述进气系统状态点参数方程组中，已知当前态点的进气空气流量，当前态点的补氧流量，当前态点的加热器供油流量，当前态点的进气总压，当前态点的进气总温，当前态点的当前进气总温，当前态点的过热空气系数，和所需模拟状态点的进气总压，所需模拟状态点的进气总温，求解与所需模拟状态点的进气总压和所需模拟状态点的进气总温对应的所需模拟状态点的进气空气流量，所需模拟状态点的补氧流量，所需模拟状态点的加热器供油流量。
2.根据权利要求1所述的试车台状态点进气模拟控制方法，其特征在于，所述进气系统状态点参数方程组中，将式替换为
3.根据权利要求1或2所述的试车台状态点进气模拟控制方法，其特征在于，所述进气系统状态点参数方程组中，
4.根据权利要求1所述的试车台状态点进气模拟控制方法，所述进气系统状态点参数方程组中，所述进气总温为燃烧后的温度，所述当前进气总温为进一步燃烧前温度。
5.根据权利要求1所述的试车台状态点进气模拟控制方法，所述试车台状态点进气模拟控制方法可用于直连式试车台。
6.根据权利要求1所述的试车台状态点进气模拟控制方法，所述试车台状态点进气模拟控制方法可用于发动机状态点进气模拟。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511017589.7A CN105628388B (zh) | 2015-12-30 | 2015-12-30 | 一种试车台状态点进气模拟控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511017589.7A CN105628388B (zh) | 2015-12-30 | 2015-12-30 | 一种试车台状态点进气模拟控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105628388A CN105628388A (zh) | 2016-06-01 |
| CN105628388B true CN105628388B (zh) | 2018-01-05 |
Family
ID=56043527
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511017589.7A Active CN105628388B (zh) | 2015-12-30 | 2015-12-30 | 一种试车台状态点进气模拟控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105628388B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108458852B (zh) * | 2018-05-24 | 2024-03-29 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高温风洞快速变温变压装置及变温变压方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19735318A1 (de) * | 1997-08-14 | 1999-02-18 | Bayerische Motoren Werke Ag | Verfahren und Vorrichtung zur Überprüfung der Funktionsfähigkeit einer Sekundärluftpumpe bei einer Brennkraftmaschine |
| GB2362226B (en) * | 2000-05-13 | 2004-06-16 | Ford Global Tech Inc | Feed-forward observer-based control for estimating cylinder air charge |
| CN101479143B (zh) * | 2007-03-01 | 2012-04-25 | 桂林吉星电子等平衡动力有限公司 | 发动机伺服加载装置及其动态寻优运行控制方法 |
| EP2634396A2 (en) * | 2012-02-28 | 2013-09-04 | General Electric Company | Sensor-based performance-seeking gas turbine engine control |
| JP2014224786A (ja) * | 2013-05-17 | 2014-12-04 | 株式会社小野測器 | エンジンモデル作成方法およびエンジンモデル作成システム |
-
2015
- 2015-12-30 CN CN201511017589.7A patent/CN105628388B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19735318A1 (de) * | 1997-08-14 | 1999-02-18 | Bayerische Motoren Werke Ag | Verfahren und Vorrichtung zur Überprüfung der Funktionsfähigkeit einer Sekundärluftpumpe bei einer Brennkraftmaschine |
| GB2362226B (en) * | 2000-05-13 | 2004-06-16 | Ford Global Tech Inc | Feed-forward observer-based control for estimating cylinder air charge |
| CN101479143B (zh) * | 2007-03-01 | 2012-04-25 | 桂林吉星电子等平衡动力有限公司 | 发动机伺服加载装置及其动态寻优运行控制方法 |
| EP2634396A2 (en) * | 2012-02-28 | 2013-09-04 | General Electric Company | Sensor-based performance-seeking gas turbine engine control |
| JP2014224786A (ja) * | 2013-05-17 | 2014-12-04 | 株式会社小野測器 | エンジンモデル作成方法およびエンジンモデル作成システム |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105628388A (zh) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105628389B (zh) | 一种试车台状态点模拟控制方法 | |
| CN101907043B (zh) | 高频燃烧不稳定性全过程模拟试验自动调节系统及方法 | |
| CN105628387B (zh) | 采用预测控制法调节试车台进气状态参数的调试方法 | |
| CN105486511B (zh) | 一种试车台状态参数的调试方法 | |
| CN102023096B (zh) | 一种航空活塞发动机内流高空模拟试验装置及其试验方法 | |
| JP5858885B2 (ja) | ガスタービン制御装置及びガスタービンを制御するための方法 | |
| CN104948304B (zh) | 一种航空燃气涡轮发动机加速供油方法 | |
| EP2795084B1 (en) | Control of the gas composition in a gas turbine power plant with flue gas recirculation | |
| CN105445007B (zh) | 一种燃气轮机用气体燃料控制阀流量特性试验系统及方法 | |
| CN108843451B (zh) | 燃气轮机燃烧室出口温度计算方法 | |
| CN107166361B (zh) | 火电机组引风机故障减负荷过程中送风机自动控制系统及方法 | |
| RU2015139836A (ru) | Способ работы газотурбинной установкой со ступенчатым и/или последовательным сгоранием | |
| EP2940275A1 (en) | Method and apparatus for controlling egr | |
| CN103195730B (zh) | 使用压气机进口可转导叶控制燃气轮机排气温度的方法 | |
| BR102015016886B1 (pt) | Motor com conjunto regulador | |
| CN106839746B (zh) | 一种高炉喷煤系统烟气炉全自动优化烧炉控制方法 | |
| CN103674463A (zh) | 空气炮发射控制系统 | |
| CN103939938B (zh) | 脉冲燃烧式燃气及助燃气压力前馈优化控制方法 | |
| CN105628388B (zh) | 一种试车台状态点进气模拟控制方法 | |
| CN203310632U (zh) | 一种发动机增压系统模拟装置 | |
| CN209027800U (zh) | 一种发动机试车台 | |
| CN105651517B (zh) | 一种并联的试车台进气系统压力调节阀的使用方法 | |
| CN105547703B (zh) | 一种试车台状态点进气模拟调节方法 | |
| CN105628391B (zh) | 一种试车台进气系统串联的压力调节阀的使用方法 | |
| CN109141906A (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 |