CN109376445B - 燃气涡轮发动机起动建模方法 - Google Patents
燃气涡轮发动机起动建模方法 Download PDFInfo
- Publication number
- CN109376445B CN109376445B CN201811317089.9A CN201811317089A CN109376445B CN 109376445 B CN109376445 B CN 109376445B CN 201811317089 A CN201811317089 A CN 201811317089A CN 109376445 B CN109376445 B CN 109376445B
- Authority
- CN
- China
- Prior art keywords
- calculating
- speed
- engine
- moment
- rotation
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Engines (AREA)
Abstract
本发明涉及一种燃气涡轮发动机起动建模方法,涉及燃气涡轮发动机建模仿真技术领域。本发明提供一种基于相似理论的燃气涡轮发动机起动建模方法,该方法包括:基于相似理论的高转转加速度、高转转速、低转转速和高压压气机出口压力计算方法;以及使用上述计算结果根据核心机气动热力学过程计算涡轮级间温度的方法。该方法能够以较高精度模拟发动机起动过程的主要参数,仅基于发动机试验数据建立发动机起动模型,具有方法简便,计算量小,对建模前提条件要求低的优点,模型精度可满足工程使用的需要。
Description
技术领域
本发明涉及燃气涡轮发动机建模仿真技术领域,具体涉及一种燃气涡轮发动机起动建模方法。
背景技术
在航空燃气涡轮发动机起动控制规律和控制算法的设计开发过程中,需要较为准确的发动机起动模型。如果采用传统的机理建模方法,需要通过试验精确测定出起动过程中的气动部件特性、燃烧效率特性、风车特性和内外涵匹配特性等,需要专门的试验设备、较高的试验成本和较长的周期,对于许多小型发动机而言根本不具备以上试验条件。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何设计一种简便易行,计算量小,对建模前提条件要求低,模型精度可满足工程使用的需要的燃气涡轮发动机起动建模方法。
(二)技术方案
为了解决上述技术问题,本发明提供了一种燃气涡轮发动机起动建模方法,包括以下步骤:
高转折合转加速度是高转折合转速N2c和折合燃油流量Wfc的函数,即有对于i=1,2,3...n次试验,其中第i次试验中记录了j=1,2,3...mi个采样时刻的数据,mi为第i次试验中采样时刻总数,采样间隔为dt,对于第i次试验中的第j个采样时刻,进行以下计算:
根据发动机风扇入口总温Tt1i,j计算温度折合系数θi,j
θi,j=Tt1i,j/Tstd (1)
其中Tstd为海平面标准大气温度288.15K;
根据发动机风扇入口总压Pt1i,j计算压力折合系数δi,j
δi,j=Pt1i,j/Pstd (2)
其中Pstd为海平面标准大气压力10132.5Pa;
使用温度折合系数θi,j、高转转速N2i,j计算高转折合转速N2ci,j
使用温度折合系数θi,j、压力折合系数δi,j、燃油流量Wfi,j计算折合燃油流量Wfci,j为:
Ni,j+1为第i次试验中的第j+1个采样时刻的发动机转速,Ni,j-1为第i次试验中的第j-1个采样时刻的发动机转速;
步骤2、根据发动机试验数据提取低转折合转速N1c变化规律:
低转折合转速N1c是高转折合转速N2c的函数,即有N1c=g(N2c);
使用第i次试验中的第j时刻的低转转速N1i,j计算高转折合转速N1ci,j
使用发动机地面起动试验数据计算得到的N2ci,j和N1ci,j,进行多项式拟合,即可得到N1c=g(N2c);
步骤3、根据发动机试验数据提取高压压气机出口折合压力Pt3c变化规律:
高压压气机出口折合压力Pt3c是高转折合转速N2c的函数,即有Pt3c=h(N2c);
使用第i次试验中的第j时刻的高压压气机出口压力Pt3i,j计算高压压气机出口折合压力Pt3ci,j
Pt3ci,j=Pt3i,j/δi,j (11)
使用发动机地面起动试验数据计算得到的N2ci,j和Pt3ci,j,进行多项式拟合,即可得到Pt3c=h(N2c);
然后利用步骤1至3的规律建立发动机起动过程建模。
优选地,慢车过程结束时仿真结束。
优选地,发动机转速达到慢车转速时仿真结束。
优选地,然后利用步骤1至3的规律建立发动机起动过程建模的计算过程包括以下步骤:
步骤4、模型初始化:高转转速的初始值N20=N2ignit,其中,N2ignit为发动机点火转速;
对于第j时刻根据该时刻模型的输入参数:发动机风扇入口总温Tt1j、发动机风扇入口总压Pt1j、燃油流量Wfj和起动机有效输出功率Psj进行步骤5的计算;
步骤5、计算j时刻温度折合系数θj和压力折合系数δj
根据Tt1j计算θj为:
θj=Tt1j/Tstd (12)
根据Pt1j计算δj为:
δj=Pt1j/Pstd (13)
根据j时刻高转转速N2j计算j时刻高转折合转速N2cj为:
利用拟合得到的低转折合转速特性N1c=g(N2c)计算j时刻低转折合转速为:
N1cj=g(N2cj) (15)
计算j时刻低转转速N1j为:
利用拟合得到的高压压气机出口折合压力特性Pt3c=h(N2c)计算j时刻高压压气机折合压力Pt3cj为:
Pt3cj=h(N2cj) (17)
计算j时刻高压压气机压力Pt3j为:
Pt3j=Pt3cjδj (18)
计算j时刻折合燃油流量Wfcj为:
利用拟合得到的高转折合转加速度变化规律计算j时刻高转折合转加速度,在风车起动时为:
其中Psj为j时刻的起动机有效输出功率;
计算j+1时刻高转转速N2j+1
之后重复步骤5的计算过程,直至仿真结束。
优选地,所述燃气涡轮发动机为航空燃气涡轮发动机。
优选地,所述燃气涡轮发动机为涡扇发动机。
(三)有益效果
本发明提供一种基于相似理论的燃气涡轮发动机起动建模方法,该方法包括:基于相似理论的高转转加速度、高转转速、低转转速和高压压气机出口压力计算方法;以及使用上述计算结果根据核心机气动热力学过程计算涡轮级间温度的方法。该方法能够以较高精度模拟发动机起动过程的主要参数,仅基于发动机试验数据建立发动机起动模型,具有方法简便,计算量小,对建模前提条件要求低的优点,模型精度可满足工程使用的需要。
附图说明
图1是采用本发明的方法得到的3km高度风车起动过程仿真结果;
图2是采用本发明的方法得到的4.5km高度风车起动过程仿真结果。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本发明实施例以涡扇发动机为例,提出一种基于相似理论的燃气涡轮起动建模方法,包括以下步骤:
高转折合转加速度是高转折合转速N2c和折合燃油流量Wfc的函数,即有对于i=1,2,3...n次试验,其中第i次试验中记录了j=1,2,3...mi个采样时刻的数据(mi为第i次试验中采样时刻总数),采样间隔为dt。对于第i次试验中的第j个采样时刻,进行以下计算:
根据发动机风扇入口总温Tt1i,j计算温度折合系数θi,j
θi,j=Tt1i,j/Tstd (1)
其中Tstd为海平面标准大气温度288.15K。
根据发动机风扇入口总压Pt1i,j计算压力折合系数δi,j
δi,j=Pt1i,j/Pstd (2)
其中Pstd为海平面标准大气压力10132.5Pa。
使用温度折合系数θi,j、高转转速N2i,j计算高转折合转速N2ci,j
使用温度折合系数θi,j、压力折合系数δi,j、燃油流量Wfi,j计算折合燃油流量Wfci,j为:
Ni,j+1为第i次试验中的第j+1个采样时刻的发动机转速,Ni,j-1为第i次试验中的第j-1个采样时刻的发动机转速;
步骤2、根据发动机试验数据提取低转折合转速N1c变化规律
低转折合转速N1c是高转折合转速N2c的函数,即有N1c=g(N2c)。
使用第i次试验中的第j时刻的低转转速N1i,j计算高转折合转速N1ci,j
使用发动机地面起动试验数据计算得到的N2ci,j和N1ci,j,进行多项式拟合,即可得到N1c=g(N2c)。
步骤3、根据发动机试验数据提取高压压气机出口折合压力Pt3c变化规律。
高压压气机出口折合压力Pt3c是高转折合转速N2c的函数,即有Pt3c=h(N2c)。
使用第i次试验中的第j时刻的高压压气机出口压力Pt3i,j计算高压压气机出口折合压力Pt3ci,j
Pt3ci,j=Pt3i,j/δi,j (11)
使用发动机地面起动试验数据计算得到的N2ci,j和Pt3ci,j,进行多项式拟合,即可得到Pt3c=h(N2c)。
在完成上述规律的提取后,使用上述规律建立发动机起动过程建模,模型计算过程包括以下步骤:
步骤4、模型初始化:高转转速的初始值N20=N2ignit,其中,N2ignit为发动机点火转速。
对于第j时刻根据该时刻模型的输入参数:发动机风扇入口总温Tt1j、发动机风扇入口总压Pt1j、燃油流量Wfj和起动机有效输出功率Psj进行步骤五的计算。
步骤5、计算j时刻温度折合系数θj和压力折合系数δj
根据Tt1j计算θj为:
θj=Tt1j/Tstd (12)
根据Pt1j计算δj为:
δj=Pt1j/Pstd (13)
根据j时刻高转转速N2j计算j时刻高转折合转速N2cj为:
利用拟合得到的低转折合转速特性N1c=g(N2c)计算j时刻低转折合转速为:
N1cj=g(N2cj) (15)
计算j时刻低转转速N1j为:
利用拟合得到的高压压气机出口折合压力特性Pt3c=h(N2c)计算j时刻高压压气机折合压力Pt3cj为:
Pt3cj=h(N2cj) (17)
计算j时刻高压压气机压力Pt3j为:
Pt3j=Pt3cjδj (18)
计算j时刻折合燃油流量Wfcj为:
利用拟合得到的高转折合转加速度变化规律计算j时刻高转折合转加速度。在风车起动时为:
其中Psj为j时刻的起动机有效输出功率。
计算j+1时刻高转转速N2j+1
之后重复步骤5的计算过程,直至仿真结束,本实施例中,发动机转速达到慢车转速时,即慢车过程结束时仿真结束。
采用本发明的方法进行3km高度风车起动的仿真,仿真与试验对比结果如图1所示。采用本发明的方法还进行了4.5km高度地面起动的仿真,仿真与试验对比结果如图2所示。仿真从点火时刻开始运行,运行周期间隔为0.01s。可以看出,N2和Pt3的模型示值误差可以控制在10%以内;N1起动初期误差较大,在第8s后示值误差也可以降低到10%以下。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.一种燃气涡轮发动机起动建模方法,其特征在于,包括以下步骤:
高转折合转加速度是高转折合转速N2c和折合燃油流量Wfc的函数,即有对于i=1,2,3...n次试验,其中第i次试验中记录了j=1,2,3...mi个采样时刻的数据,mi为第i次试验中采样时刻总数,采样间隔为dt,对于第i次试验中的第j个采样时刻,进行以下计算:
根据发动机风扇入口总温Tt1i,j计算温度折合系数θi,j
θi,j=Tt1i,j/Tstd (1)
其中Tstd为海平面标准大气温度288.15K;
根据发动机风扇入口总压Pt1i,j计算压力折合系数δi,j
δi,j=Pt1i,j/Pstd (2)
其中Pstd为海平面标准大气压力10132.5Pa;
使用温度折合系数θi,j、高转转速N2i,j计算高转折合转速N2ci,j
使用温度折合系数θi,j、压力折合系数δi,j、燃油流量Wfi,j计算折合燃油流量Wfci,j为:
Ni,j+1为第i次试验中的第j+1个采样时刻的发动机转速,Ni,j-1为第i次试验中的第j-1个采样时刻的发动机转速;
步骤2、根据发动机试验数据提取低转折合转速N1c变化规律:
低转折合转速N1c是高转折合转速N2c的函数,即有N1c=g(N2c);
使用第i次试验中的第j时刻的低转转速N1i,j计算高转折合转速N1ci,j
使用发动机地面起动试验数据计算得到的N2ci,j和N1ci,j,进行多项式拟合,即得到N1c=g(N2c);
步骤3、根据发动机试验数据提取高压压气机出口折合压力Pt3c变化规律:
高压压气机出口折合压力Pt3c是高转折合转速N2c的函数,即有Pt3c=h(N2c);
使用第i次试验中的第j时刻的高压压气机出口压力Pt3i,j计算高压压气机出口折合压力Pt3ci,j
Pt3ci,j=Pt3i,j/δi,j (11)
使用发动机地面起动试验数据计算得到的N2ci,j和Pt3ci,j,进行多项式拟合,即得到Pt3c=h(N2c);
然后利用步骤1至3的规律建立发动机起动过程建模。
2.如权利要求1所述的方法,其特征在于,然后利用步骤1至3的规律建立发动机起动过程建模的计算过程包括以下步骤:
步骤4、模型初始化:高转转速的初始值N20=N2ignit,其中,N2ignit为发动机点火转速;
对于第j时刻根据该时刻模型的输入参数:发动机风扇入口总温Tt1j、发动机风扇入口总压Pt1j、燃油流量Wfj和起动机有效输出功率Psj进行步骤5的计算;
步骤5、计算j时刻温度折合系数θj和压力折合系数δj
根据Tt1j计算θj为:
θj=Tt1j/Tstd (12)
根据Pt1j计算δj为:
δj=Pt1j/Pstd (13)
根据j时刻高转转速N2j计算j时刻高转折合转速N2cj为:
利用拟合得到的低转折合转速特性N1c=g(N2c)计算j时刻低转折合转速为:
N1cj=g(N2cj) (15)
计算j时刻低转转速N1j为:
利用拟合得到的高压压气机出口折合压力特性Pt3c=h(N2c)计算j时刻高压压气机折合压力Pt3cj为:
Pt3cj=h(N2cj) (17)
计算j时刻高压压气机压力Pt3j为:
Pt3j=Pt3cjδj (18)
计算j时刻折合燃油流量Wfcj为:
利用拟合得到的高转折合转加速度变化规律计算j时刻高转折合转加速度,在风车起动时为:
在地面起动时,先计算未考虑起动机作用的高转净折合转加速度
其中Psj为j时刻的起动机有效输出功率;
计算j+1时刻高转转速N2j+1
之后重复步骤5的计算过程,直至仿真结束。
3.如权利要求1所述的方法,其特征在于,慢车过程结束时仿真结束。
4.如权利要求3所述的方法,其特征在于,发动机转速达到慢车转速时仿真结束。
5.如权利要求1所述的方法,其特征在于,所述燃气涡轮发动机为航空燃气涡轮发动机。
6.如权利要求5所述的方法,其特征在于,所述燃气涡轮发动机为涡扇发动机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811317089.9A CN109376445B (zh) | 2018-11-07 | 2018-11-07 | 燃气涡轮发动机起动建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811317089.9A CN109376445B (zh) | 2018-11-07 | 2018-11-07 | 燃气涡轮发动机起动建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109376445A CN109376445A (zh) | 2019-02-22 |
CN109376445B true CN109376445B (zh) | 2022-10-14 |
Family
ID=65384116
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811317089.9A Active CN109376445B (zh) | 2018-11-07 | 2018-11-07 | 燃气涡轮发动机起动建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109376445B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110728052B (zh) * | 2019-10-11 | 2023-09-22 | 中国航发沈阳发动机研究所 | 一种确定旋转盘腔相似试验边界条件的方法 |
CN111058984B (zh) * | 2019-11-20 | 2021-04-30 | 中国人民武装警察部队工程大学 | 一种基于关联有向图的燃油供给系统测试性建模方法 |
CN112539941B (zh) * | 2020-12-02 | 2023-01-20 | 西安航天动力研究所 | 考虑真实气体效应的液体火箭发动机热试验参数设置方法 |
CN112613119B (zh) * | 2020-12-17 | 2024-04-12 | 大连理工大学 | 一种多几何参数可调的进/排/发一体化航空推进系统建模方法 |
CN113886983B (zh) * | 2021-09-27 | 2022-09-20 | 中国航发沈阳发动机研究所 | 一种基于试验数据的航空发动机起动仿真模型修正方法 |
CN114217722B (zh) * | 2021-12-22 | 2024-05-17 | 沈阳东睿科技有限公司 | 一种基于动态绑定的发动机卷宗数据自动采集方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103838920A (zh) * | 2014-02-13 | 2014-06-04 | 南京航空航天大学 | 一种涡轴发动机动态实时模型建模方法及系统 |
CN106323640A (zh) * | 2015-06-30 | 2017-01-11 | 中航空天发动机研究院有限公司 | 一种航空发动机加减速供油试验方法 |
CN106569977A (zh) * | 2016-09-26 | 2017-04-19 | 南京航空航天大学 | 涡轴发动机起动过程建模方法 |
CN107315875A (zh) * | 2017-06-26 | 2017-11-03 | 南京航空航天大学 | 分开排气三涵道涡扇发动机仿真模型 |
CN108647428A (zh) * | 2018-05-08 | 2018-10-12 | 南京航空航天大学 | 一种涡扇发动机自适应部件级仿真模型构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7233884B2 (en) * | 2002-10-31 | 2007-06-19 | United Technologies Corporation | Methodology for temporal fault event isolation and identification |
-
2018
- 2018-11-07 CN CN201811317089.9A patent/CN109376445B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103838920A (zh) * | 2014-02-13 | 2014-06-04 | 南京航空航天大学 | 一种涡轴发动机动态实时模型建模方法及系统 |
CN106323640A (zh) * | 2015-06-30 | 2017-01-11 | 中航空天发动机研究院有限公司 | 一种航空发动机加减速供油试验方法 |
CN106569977A (zh) * | 2016-09-26 | 2017-04-19 | 南京航空航天大学 | 涡轴发动机起动过程建模方法 |
CN107315875A (zh) * | 2017-06-26 | 2017-11-03 | 南京航空航天大学 | 分开排气三涵道涡扇发动机仿真模型 |
CN108647428A (zh) * | 2018-05-08 | 2018-10-12 | 南京航空航天大学 | 一种涡扇发动机自适应部件级仿真模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109376445A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109376445B (zh) | 燃气涡轮发动机起动建模方法 | |
CN106569977B (zh) | 涡轴发动机起动过程建模方法 | |
CN110717219A (zh) | 航空发动机整机状态下压气机进口流量的获取方法及装置 | |
CN110348078B (zh) | 一种涡轴发动机容积动力学结合热惯性效应的建模方法 | |
CN108108528A (zh) | 一种分轴式航改燃机动力涡轮一维匹配设计方法 | |
Benson et al. | Comparison of experimental and simulated transient responses of a turbocharged diesel engine | |
CN109829238B (zh) | 基于折合平衡流形展开模型的航空发动机系统辨识方法 | |
Gao et al. | Modeling and simulation of an aero turbojet engine with gasturb | |
CN109635512A (zh) | 一种基于修正控制方程的离心叶轮进口设计方法 | |
Aghaei-Togh et al. | Effects of nozzle arrangement angle on the performance of partially admitted turbines | |
Reitz et al. | Full High Pressure Compressor Investigations to Determine Aerodynamic Changes due to Deterioration | |
CN107451331A (zh) | 一种增压柴油机工作过程三维仿真方法 | |
Mataczynski et al. | Design, Simulation, and Testing of a Pressure Wave Supercharger for a Small Internal Combustion Engine | |
CN115419509A (zh) | 一种燃气涡轮发动机动态模型 | |
Riegler et al. | Validation of a mixed flow turbofan performance model in the sub-idle operating range | |
Ma et al. | Modeling and simulation of gas turbine starter and fuel control system | |
Schobeiri | Active aerodynamic control of multi-stage axial compressor instability and surge by dynamically adjusting the stator blades | |
Tan et al. | Turbo engine starting control law design and process simulation | |
Pakanati et al. | High fidelity engine performance models for windmill relight predictions | |
Agrawal et al. | A generalized mathematical model to estimate gas turbine starting characteristics | |
Castner et al. | An engine research program focused on low pressure turbine aerodynamic performance | |
Li et al. | Calculation Model Based Design-Point Gas Generator Performance Adaptation Method | |
Hao et al. | Modeling of aero-engine starting process based on rising rate of rotating speed | |
Krivosheev et al. | Analysis of Options for Converting Aviation Two Spool Turbojet Engines with Afterburner when Developing Gas-Turbine-Driven Compressor Plant for Gas-Compressor Unit | |
Biraud et al. | Simulation of the WR-21 Advanced Cycle Engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |