CN104408243A - 一种发动机模型安装推力修正方法 - Google Patents
一种发动机模型安装推力修正方法 Download PDFInfo
- Publication number
- CN104408243A CN104408243A CN201410665841.4A CN201410665841A CN104408243A CN 104408243 A CN104408243 A CN 104408243A CN 201410665841 A CN201410665841 A CN 201410665841A CN 104408243 A CN104408243 A CN 104408243A
- Authority
- CN
- China
- Prior art keywords
- model
- engine
- thrust
- loss
- pla
- 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.)
- Granted
Links
Landscapes
- Testing Of Engines (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于航空技术领域动力装置系统,涉及一种发动机模型安装推力修正方法。所述方法包括:1.1)得到发动机台架推力Fi数据表格;1.2)采用三维插值算法构建发动机台架推力模型Ft;1.3)开展进气道吹风试验;1.4)使用插值方法建模;1.5)建立进气道损失修正模型K1;1.6)建立功率提取损失系数数据库;1.7)构建功率提取损失修正模型K2;1.8)建立喷管后体阻力损失系数数据库;1.9)构建喷管后体阻力损失修正模型K3;1.10)完成发动机安装推力模型的建模。本发明计算精度高,仿真具有实时性,易于实现,建模周期短,工程性应用性好。
Description
技术领域
本发明属于航空技术领域动力装置系统,涉及一种发动机模型安装推力修正方法。
背景技术
飞机设计中,经常用发动机仿真模型,根据性能仿真结果进行飞机仿真验证。发动机仿真模型为发动机装机后的包含进排气损失的发动机安装推力性能模型。
现有发动机台架推力模型采用部件法建模,其特点是已有各部件特性,选定控制方案,按照各部件共同工作条件确定共同工作点,再计算发动机推力和耗油率。由于发动机各部件之间相互制约关系比较复杂,在求共同工作点时要对多个参数的数值进行反复迭代,仿真周期长。
现有进气道损失模型,通常计算进气道的溢流阻力、附面层抽吸阻力、旁路放气阻力和辅助进气阻力,通过迭代求解方程组得到进气道损失,仿真周期长。
现有功率提取损失模型,通常采用经验公式修正得到发动机功率提取损失,仿真精度低。
现有喷管后体阻力损失模型,通常计算喷管外流附面层分离引起压力阻力,结合发动机部件级模型通过迭代求解方程组得到发动机喷管后体阻力损失,仿真周期同样很长。
这种方法建立的发动机安装推力仿真模型,仿真精度一般,仿真周期长达数秒,无法满足实时仿真的要求。
发明内容
本发明所要解决的技术问题是:提供一种仿真精度高、周期短的发动机模型安装推力修正方法。
本发明采取的技术方案为:一种发动机模型安装推力方法,其特征为所述方法包括以下步骤:
1.1)开展发动机地面台架试验和高空台试验,给定飞机飞行状态和油门杆角度PLA,得到发动机台架推力Fi数据表格;所述飞机飞行状态包括飞行高度H和飞行马赫数Ma;
1.2)使用插值方法建模,基于发动机推力数据表格,采用三维插值算法构建发动机台架推力模型Ft:
Ft=f(H,Ma,PLA)
式中:
飞行高度H;
飞行马赫数Ma;
油门杆角度PLA;
1.3)开展进气道吹风试验,测得不同飞行马赫数Ma下进气道总压恢复系数σin:
σin=P1/P0
式中:
进气道出口截面总压P1;
发动机远前方总压P0;
1.4)使用插值方法建模,基于试验测得的进气道总压恢复系数,采用一维插值算法构建进气道总压恢复系数模型σ:
σ=f(Ma)
1.5)依据进气道总压恢复系数,建立进气道损失修正模型K1:
K1=[1-(1-σ)(1+PH·Ac/Ft)]
式中:
PH—给定高度的大气压力(通过查标准大气数据表得到);
Ac—给定发动机状态的尾喷口喉道截面面积;
1.6)通过功率提取试验测得不同功率W提取条件下发动机推力Fw,建立功率提取损失系数数据库η:
Fw=f(W,H,Ma,PLA)
η=(Ft-Fw)/Ft
1.7)基于功率提取损失系数数据库,构建功率提取损失修正模型K2:
K2=1-η
1.8)通过飞机后体缩比模型吹风试验测得不同飞行状态(飞行高度H和飞行马赫数Ma)和油门杆角度PLA条件下喷管后体阻力损失系数,建立喷管后体阻力损失系数数据库ΔPc:
ΔPc=f(H,Ma,PLA)
1.9)基于喷管后体阻力损失系数数据库,构建喷管后体阻力损失修正模型K3:
K3=1-ΔPc
1.10)结合发动机台架推力模型Ft、进气道损失模型K1、功率提取损失模型K2和喷管后体阻力损失模型K3完成发动机安装推力模型Fa的建模:
Fa=Ft·K1·K2·K3。
本发明的优点是:
1)考虑影响因素全面,可以修正发动机装机后的进气、功率提取、排气导致的推力损失;
2)计算精度高,误差不超过5%,可以精确模拟发动机性能参数;
3)仿真具有实时性,仿真周期小于10ms,可以支持飞机进行实时飞行仿真;
4)软件开发简单,易于实现,建模周期短,工程性应用性好;
5)通用性好,可以应用于各类喷气式发动机仿真建模。
本发明主要应用于涡轮风扇发动机仿真模型安装推力修正。发动机安装推力修正具有重要的工程价值,可为飞/发一体化性能设计、飞机性能仿真设计等工作提供重要依据,从而降低风险。各类喷气式发动机仿真模型安装推力修正均可按照本发明提供的方法实施,因此具有广泛的应用性。
附图说明
图1是发动机台架推力模型建模流程示意图。
图2是进气道损失修正模型建模流程示意图。
图3是功率提取损失修正模型建模流程示意图。
图4是喷管后体阻力损失修正模型建模流程示意图。
图5是发动机安装推力模型建模流程示意图。
具体实施方式
下面通过具体的实施例并结合附图对本发明作进一步详细的描述。
实施例:本例是美国F110涡轮风扇发动机配装F16飞机,采用试验与插值方法相结合的修正方法,对发动机性能仿真中使用的发动机性能模型进行装机后的安装推力性能修正。具体步骤如下:
1)开展发动机地面台架试验和高空台试验,给定飞机飞行状态(飞行高度H和飞行马赫数Ma)和油门杆角度PLA,得到发动机台架推力Fi数据表格;设计点H=0km、Ma=0、PLA=110°,试验测得发动机台架推力Fi=105.9kN:
2)使用插值方法建模,基于发动机推力数据表格,采用三维插值算法构建发动机台架推力模型Ft,请参阅图1:
Ft=f(H,Ma,PLA)
式中:
飞行高度H;
飞行马赫数Ma;
油门杆角度PLA;
3)开展进气道吹风试验,测得不同飞行马赫数Ma下进气道总压恢复系数σin:
σin=P1/P0
式中:
进气道出口截面总压P1;
发动机远前方总压P0;
设计点H=0km、Ma=0,试验测得进气道总压恢复系数σin=1;
4)使用插值方法建模,基于试验测得的进气道总压恢复系数,采用一维插值算法构建进气道总压恢复系数模型σ:
σ=f(Ma)
5)依据进气道总压恢复系数,建立进气道损失修正模型K1,请参阅图2:
K1=[1-(1-σ)(1+PH·Ac/Ft)]
式中:
PH—给定高度的大气压力(通过查标准大气数据表得到);
Ac—给定发动机状态的尾喷口喉道截面面积;
设计点H=0km、Ma=0,K1=1;
6)通过功率提取试验测得不同功率W提取条件下发动机推力Fw,建立功率提取损失系数数据库η:
Fw=f(W,H,Ma,PLA)
η=(Ft-Fw)/Ft
设计点H=0km、Ma=0、功率160kW,Fw=104.8kN,η=0.0104;
7)基于功率提取损失系数数据库,构建功率提取损失修正模型K2,请参阅图3:
K2=1-η
设计点H=0km、Ma=0,K2=0.9896;
8)通过飞机后体缩比模型吹风试验测得不同飞行状态(飞行高度H和飞行马赫数Ma)和油门杆角度PLA条件下喷管后体阻力损失系数,建立喷管后体阻力损失系数数据库ΔPc:
ΔPc=f(H,Ma,PLA)
设计点H=0km、Ma=0、PLA=110°,ΔPc=0.11;
9)基于喷管后体阻力损失系数数据库,构建喷管后体阻力损失修正模型K3,请参阅图4:
K3=1-ΔPc
设计点H=0km、Ma=0,K3=0.89;
10)结合发动机台架推力模型Ft、进气道损失模型K1、功率提取损失模型K2和喷管后体阻力损失模型K3完成发动机安装推力模型Fa的建模,请参阅图5:
Fa=Ft·K1·K2·K3。
设计点H=0km、Ma=0、PLA=110°,安装推力Fa=93.2kN。
Claims (1)
1.一种发动机模型安装推力方法,其特征为所述方法包括以下步骤:
1.1)开展发动机地面台架试验和高空台试验,给定飞机飞行状态和油门杆角度PLA,得到发动机台架推力Fi数据表格;所述飞机飞行状态包括飞行高度H和飞行马赫数Ma;
1.2)使用插值方法建模,基于发动机推力数据表格,采用三维插值算法构建发动机台架推力模型Ft:
Ft=f(H,Ma,PLA)
式中:
飞行高度H;
飞行马赫数Ma;
油门杆角度PLA;
1.3)开展进气道吹风试验,测得不同飞行马赫数Ma下进气道总压恢复系数σin:
σin=P1/P0
式中:
进气道出口截面总压P1;
发动机远前方总压P0;
1.4)使用插值方法建模,基于试验测得的进气道总压恢复系数,采用一维插值算法构建进气道总压恢复系数模型σ:
σ=f(Ma)
1.5)依据进气道总压恢复系数,建立进气道损失修正模型K1:
K1=[1-(1-σ)(1+PH·Ac/Ft)]
式中:
PH—给定高度的大气压力(通过查标准大气数据表得到);
Ac—给定发动机状态的尾喷口喉道截面面积;
1.6)通过功率提取试验测得不同功率W提取条件下发动机推力Fw,建立功率提取损失系数数据库η:
Fw=f(W,H,Ma,PLA)
η=(Ft-Fw)/Ft
1.7)基于功率提取损失系数数据库,构建功率提取损失修正模型K2:
K2=1-η
1.8)通过飞机后体缩比模型吹风试验测得不同飞行状态(飞行高度H和飞行马赫数Ma)和油门杆角度PLA条件下喷管后体阻力损失系数,建立喷管后体阻力损失系数数据库ΔPc:
ΔPc=f(H,Ma,PLA)
1.9)基于喷管后体阻力损失系数数据库,构建喷管后体阻力损失修正模型K3:
K3=1-ΔPc
1.10)结合发动机台架推力模型Ft、进气道损失模型K1、功率提取损失模型K2和喷管后体阻力损失模型K3完成发动机安装推力模型Fa的建模:
Fa=Ft·K1·K2·K3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410665841.4A CN104408243B (zh) | 2014-11-19 | 2014-11-19 | 一种发动机模型安装推力修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410665841.4A CN104408243B (zh) | 2014-11-19 | 2014-11-19 | 一种发动机模型安装推力修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408243A true CN104408243A (zh) | 2015-03-11 |
CN104408243B CN104408243B (zh) | 2019-02-12 |
Family
ID=52645874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410665841.4A Active CN104408243B (zh) | 2014-11-19 | 2014-11-19 | 一种发动机模型安装推力修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408243B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107991010A (zh) * | 2017-11-24 | 2018-05-04 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种发动机上多线性约束拉紧力测量装置 |
CN110259583A (zh) * | 2019-07-19 | 2019-09-20 | 中国航发沈阳发动机研究所 | 双外涵变循环航空发动机及其功率的控制方法、控制器 |
CN111717411A (zh) * | 2020-05-22 | 2020-09-29 | 成都飞机工业(集团)有限责任公司 | 一种基于试飞数据标准重量下巡航推力增量的修正方法 |
CN112133156A (zh) * | 2020-09-25 | 2020-12-25 | 中国直升机设计研究所 | 一种三发直升机单发失效训练方法 |
CN113188799A (zh) * | 2021-04-27 | 2021-07-30 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于速度差极值法的航空发动机推力修正方法 |
US11667392B2 (en) | 2019-06-20 | 2023-06-06 | Pratt & Whitney Canada Corp. | Method and system for operating a rotorcraft engine |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080073609A1 (en) * | 2006-09-23 | 2008-03-27 | Eldert Akkermann | Compressed-air needle valve for controlling an air flow for driving engine simulators in aircraft models for wind tunnel experiments |
CN101445106A (zh) * | 2007-11-07 | 2009-06-03 | 通用汽车环球科技运作公司 | 控制混合动力系排气后处理系统温度的方法和装置 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
CN103192996A (zh) * | 2013-04-24 | 2013-07-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种大型飞机飘降航迹的确定方法 |
CN103258105A (zh) * | 2012-02-17 | 2013-08-21 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机起落架寿命消耗计算方法 |
CN103593501A (zh) * | 2013-10-11 | 2014-02-19 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种航空发动机滑油系统设计方法 |
CN103984848A (zh) * | 2014-04-17 | 2014-08-13 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机拦阻装置载荷特性曲线的修正方法 |
-
2014
- 2014-11-19 CN CN201410665841.4A patent/CN104408243B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080073609A1 (en) * | 2006-09-23 | 2008-03-27 | Eldert Akkermann | Compressed-air needle valve for controlling an air flow for driving engine simulators in aircraft models for wind tunnel experiments |
CN101445106A (zh) * | 2007-11-07 | 2009-06-03 | 通用汽车环球科技运作公司 | 控制混合动力系排气后处理系统温度的方法和装置 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
CN103258105A (zh) * | 2012-02-17 | 2013-08-21 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机起落架寿命消耗计算方法 |
CN103192996A (zh) * | 2013-04-24 | 2013-07-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种大型飞机飘降航迹的确定方法 |
CN103593501A (zh) * | 2013-10-11 | 2014-02-19 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种航空发动机滑油系统设计方法 |
CN103984848A (zh) * | 2014-04-17 | 2014-08-13 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机拦阻装置载荷特性曲线的修正方法 |
Non-Patent Citations (2)
Title |
---|
刘鹏超等: "某型飞机发动机的安装推力计算", 《航空工程进展》 * |
赵肃等: "航空发动机推力变化对飞机基本性能影响的敏感性分析", 《航空发动机》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107991010A (zh) * | 2017-11-24 | 2018-05-04 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种发动机上多线性约束拉紧力测量装置 |
CN107991010B (zh) * | 2017-11-24 | 2020-09-25 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种发动机上多线性约束拉紧力测量装置 |
US11667392B2 (en) | 2019-06-20 | 2023-06-06 | Pratt & Whitney Canada Corp. | Method and system for operating a rotorcraft engine |
CN110259583A (zh) * | 2019-07-19 | 2019-09-20 | 中国航发沈阳发动机研究所 | 双外涵变循环航空发动机及其功率的控制方法、控制器 |
CN111717411A (zh) * | 2020-05-22 | 2020-09-29 | 成都飞机工业(集团)有限责任公司 | 一种基于试飞数据标准重量下巡航推力增量的修正方法 |
CN112133156A (zh) * | 2020-09-25 | 2020-12-25 | 中国直升机设计研究所 | 一种三发直升机单发失效训练方法 |
CN113188799A (zh) * | 2021-04-27 | 2021-07-30 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于速度差极值法的航空发动机推力修正方法 |
CN113188799B (zh) * | 2021-04-27 | 2022-09-30 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于速度差极值法的航空发动机推力修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104408243B (zh) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408243A (zh) | 一种发动机模型安装推力修正方法 | |
Otoguro et al. | Space–time VMS flow analysis of a turbocharger turbine with isogeometric discretization: computations with time-dependent and steady-inflow representations of the intake/exhaust cycle | |
CN104392039A (zh) | 一种过渡态发动机仿真建模方法 | |
CN110937129B (zh) | 一种用于涡桨动力无人机飞行性能计算的推阻体系划分方法 | |
CN104298826A (zh) | 一种反推力状态下的航空发动机气动稳定性预测与评估方法 | |
Reitenbach et al. | Optimization of compressor variable geometry settings using multi-fidelity simulation | |
CN109460628B (zh) | 一种进气道与发动机共同工作的流量匹配评估方法 | |
CN102650565A (zh) | 风洞模拟实验中涡轮动力模拟器短舱唇口及其设计方法 | |
CN104298804A (zh) | 一种飞行载荷设计方法 | |
CN112380795A (zh) | 一种基于预置模型库的试车台流场数值仿真系统 | |
CN105354400A (zh) | 一种发动机反推力装置设计方法 | |
CN103544053B (zh) | 一种考虑织物透气性的降落伞气动性能分析方法 | |
CN105785791A (zh) | 一种超声速状态下机载推进系统的建模方法 | |
CN112257256A (zh) | 一种基于稳态数据的发动机简化动态模型设计方法 | |
CN104200107A (zh) | 不同飞行迎角对航空发动机失速/喘振特性影响的适航审定方法 | |
Smith | An Assessment of CREATE-AV Kestrel for F-35 Aero/Performance Applications | |
McWaters | F-35 conventional mode jet-effects testing methodology | |
Xie et al. | Numerical investigation of inlet distortion for different rear mounted engine installations at taking-off conditions | |
Klepper et al. | A-10 Analysis Using HPCMP CREATETM-AV Kestrel Product Utilizing the Firebolt Propulsion Component | |
Liu et al. | Simulating effects of high angle of attack on turbofan engine performance | |
CN209535509U (zh) | 一种航天推进系统的地面试验系统 | |
Hale et al. | A numerical simulation capability for analysis of aircraft inlet–engine compatibility | |
CN109900486B (zh) | 一种带飞飞行器热态气动性能获取方法 | |
Deere et al. | CFD predictions for transonic performance of the ERA hybrid wing-body configuration | |
Liu et al. | Numerical investigation on engine inlet distortion under crosswind for a commercial transport aircraft |
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 |