CN110686813A - 一种通过地面滑行试验辨识双发飞机发动机推力的方法 - Google Patents

一种通过地面滑行试验辨识双发飞机发动机推力的方法 Download PDF

Info

Publication number
CN110686813A
CN110686813A CN201910962534.5A CN201910962534A CN110686813A CN 110686813 A CN110686813 A CN 110686813A CN 201910962534 A CN201910962534 A CN 201910962534A CN 110686813 A CN110686813 A CN 110686813A
Authority
CN
China
Prior art keywords
engine
thrust
deceleration
acceleration
speed
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.)
Pending
Application number
CN201910962534.5A
Other languages
English (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.)
Chengdu Aircraft Industrial Group Co Ltd
Original Assignee
Chengdu Aircraft Industrial Group 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 Chengdu Aircraft Industrial Group Co Ltd filed Critical Chengdu Aircraft Industrial Group Co Ltd
Priority to CN201910962534.5A priority Critical patent/CN110686813A/zh
Publication of CN110686813A publication Critical patent/CN110686813A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0028Force sensors associated with force applying means
    • G01L5/0038Force sensors associated with force applying means applying a pushing force
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/14Testing gas-turbine engines or jet-propulsion engines

Abstract

本发明公开了一种通过地面滑行试验辨识双发飞机发动机推力的方法,忽略发动机转速变化而引起的进气畸变变化量;在重量m、速度V下分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验;通过加/减速运动曲线计算得到加/减速度,进而计算得到不同转速下的推力。本发明克服了发动机台架试验推力数据的差异性、空中飞行试验费用的高昂性,降低首飞的试飞风险。

Description

一种通过地面滑行试验辨识双发飞机发动机推力的方法
技术领域
本发明属于飞机发动机测试的技术领域,具体涉及一种通过地面滑行试验辨识双发飞机发动机推力的方法。
背景技术
目前,公知的发动机推力数据是通过台架试验得到的,但由于发动机的装配以及在机上的安装误差影响,台架推力与机上发动机的实际推力存在差异。为了明确机上发动机的实际推力特性,合理制定试飞规划和准确预估飞机性能指标,国内一些科研人员基于飞行试验数据对发动机模型进行辨识。通过在飞机上加装传感器,利用传感器辨识发动机的工作参数。如雷晓波等通过在推力销上加装剪力应变桥路,建立载荷标定方程和温度修正方程获取发动机安装节推力,同时,利用进气道测量耙测试参数,计算飞行中进气道冲压阻力和压差阻力。从而,辨识出发动机的推力特性(雷晓波,李密,张强,高扬,文敏.航空发动机推力直接测量飞行试验[J].航空动力学报.2018,33(7):1631-1638.)。
目前,在飞行试验阶段进行发动机推力辨识方面存在以下缺陷:
1)发动机台架推力与机上发动机推力的存在差异;
2)空中飞行试验费用高昂且风险较大;
3)特别是,新研飞机在首飞之前,就需要充分的掌握机上发动机的推力特性。
发明内容
本发明的目的在于提供一种通过地面滑行试验辨识双发飞机发动机推力的方法,旨在克服发动机台架试验推力数据的差异性、空中飞行试验费用的高昂性,降低首飞的试飞风险。
本发明主要通过以下技术方案实现:一种通过地面滑行试验辨识双发飞机发动机推力的方法,忽略发动机转速变化而引起的进气畸变变化量;在重量m、速度V下分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验;通过加/减速运动曲线计算得到加/减速度,进而计算得到不同转速下的推力。
为了更好的实现本发明,进一步的,双发同步加/减速的双发转速为n1;所述双发异步加/减速的一发转速为n1,且另一发关车。
为了更好的实现本发明,进一步的,所述推力的计算公式如下:
2Tn1-D-f=ma1 (7)
Tn1+T关车–D-f=ma2 (8)
f=F (9)
其中,a1,a2为滑行过程中的加速度,
Tn1为左发动机转速n1对应的推力,
T关车为右发动机关车后的推力,发动机关车后存在溢流阻力,
D为该滑行状态下的气动阻力,
f为该滑行状态下的摩擦阻力,
F为该滑行状态下的阻力即气动阻力和摩擦阻力的合力;
在相同重量和相同速度时,2种转速状态下的气动阻力和摩擦阻力相同,即:
2Tn1-F=ma1 (10)
Tn1+T关车–F=ma2 (11)
通过方程(7)-(11),可得:
Tn1–T关车=m(a1–a2) (12)
其中,T关车为发动机关车后的推力,即溢流阻力,此时,溢流阻力相对于发动机的推力Tn1仅占据小量,认为T关车的推力为0,方程(12)化简为
Tn1≈m(a1–a2) (13)
此时,得到了转速n1对应的推力,即Tn1
为了更好的实现本发明,进一步的,滑行试验中加/减速运动的力学模型如下:
L+N=G (1)
T-D-f=ma (2)
Figure BDA0002229419870000021
Figure BDA0002229419870000022
f=μN (5)
其中,L为飞机的升力,N为飞机轮胎的支撑力,G为飞机的重力,T为两个发动机的总推力,D为飞机的阻力,f为飞机受到的摩擦力,m为飞机的质量,a为飞机的加/减速度,ρ为密度,V为速度,S为参考面积,CL为升力系数,CD为阻力系数。
为了更好的实现本发明,进一步的,所述加/减速度计算公式如下:
Figure BDA0002229419870000023
其中,选取速度为V的滑行状态点作为不同发动机状态的基准速度,进而得到基准速度区间的加/减速度。
为了更好的实现本发明,进一步的,在重量m、速度V下改变转速分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验,测试得到相同速度下不同转速的推力。
为了更好的实现本发明,进一步的,在重量m下改变速度分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验,测试得到不同速度下不同转速的推力。
本发明的有益效果:
(1)本发明克服了发动机台架试验推力数据的差异性、空中飞行试验费用的高昂性,降低首飞的试飞风险。
(2)本发明无需加装传感器和复杂的推导,能够直接准确的测量出发动机的实际推力,帮助设计人员修正台架试验数据,掌握该发动机的性能。
(3)本发明通过辨识出的发动机推力值,能够准确掌握飞机性能,制定合理的试飞规划,避免因发动机数据掌握不准确而导致的首飞风险。
(4)本发明利用地面滑行试验进行,极大的节约了人力和财力。同时,为后续试飞规划的合理制定提供了数据支撑,保证试飞架次的顺利进行,缩短了新机试飞的周期,有利于产品更快的进入市场和占领市场。
附图说明
图1为本发明的流程图;
图2为滑行加/减速的力学模型;
图3双发同步在转速n1下的加/减速运动曲线图;
图4为双发异步在一发转速n1、另一发关车下的加/减速运动曲线图;
图5为不同状态下的加速度曲线图;
图6为发动机最大状态推力曲线图。
具体实施方式
实施例1:
一种通过地面滑行试验辨识双发飞机发动机推力的方法,忽略发动机转速变化而引起的进气畸变变化量;在重量m、速度V下分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验;通过加/减速运动曲线计算得到加/减速度,进而计算得到不同转速下的推力。
本发明克服了发动机台架试验推力数据的差异性、空中飞行试验费用的高昂性,降低首飞的试飞风险。本发明无需加装传感器和复杂的推导,能够直接准确的测量出发动机的实际推力,帮助设计人员修正台架试验数据,掌握该发动机的性能。本发明利用地面滑行试验进行,极大的节约了人力和财力。同时,为后续试飞规划的合理制定提供了数据支撑,保证试飞架次的顺利进行,缩短了新机试飞的周期,有利于产品更快的进入市场和占领市场。
实施例2:
一种通过地面滑行试验辨识双发飞机发动机推力的方法,假定条件:忽略发动机转速变化而引起的进气畸变变化量,认为不同转速时,进气畸变变化量相同。如图1所示,主要包括以下步骤:
1)选取工况,即重量和速度点;
2)分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速3种不同状态的滑行试验;
3)计算3种不同的滑行过程中,相同速度区间的加/减速度;
4)通过3种不同状态的滑行过程,分别求解不同转速对应的推力值;
5)改变转速,重复2)~4)的步骤,辨识出相同速度,不同转速对应的发动机推力值;
6)通过调整速度点,重复2)~5)的步骤,辨识出不同速度和转速对应的发动机推力值。
实施例3:
本实施例是在实施例1或2的基础上进行优化,为了使飞机的发动机推力方便求解,选取重量m,速度V的滑行状态点进行测试。如表1所示,双发同步加/减速的双发转速为n1;所述双发异步加/减速的一发转速为n1,且另一发关车。进行2种不同发动机转速状态下的飞机滑行试验。
表1 3种不同转速状态点
依据表1的转速状态点,进行滑行试验:1)双发同步加/减速(双发转速n1);2)双发异步加/减速(一发转速n1,另一发关车)。
本实施例的其他部分与上述实施例1或2相同,故不再赘述。
实施例4:
本实施例是在实施例1-3任一个的基础上进行优化,滑行加/减速运动的力学模型见图2,具体的方程组如下:
L+N=G (1)
T-D-f=ma (2)
Figure BDA0002229419870000051
Figure BDA0002229419870000052
f=μN (5)
其中,L为飞机的升力,N为飞机轮胎的支撑力,G为飞机的重力,T为两个发动机的总推力,D为飞机的阻力,f为飞机受到的摩擦力,m为飞机的质量,a为飞机滑行的加/减速度,ρ为密度,V为速度,S为参考面积,CL为升力系数,CD为阻力系数。
本实施例的其他部分与上述实施例1-3任一项相同,故不再赘述。
实施例5:
本实施例是在实施例1-4任一个的基础上进行优化,2种不同的双发转速状态下的加/减速运动过程见图3和图4。通过该加/减速运动图线,选取速度为V的滑行状态点作为2种状态的基准速度,得到该基准速度区间的加/减速度,具体的方程如下:
Figure BDA0002229419870000053
运用上述方法,得到2种不同情况的飞机加/减速度,见表2。
表2 3种不同转速状态下的加/减速度
本实施例的其他部分与上述实施例1-4任一项相同,故不再赘述。
实施例6:
本实施例是在实施例1-5任一个的基础上进行优化,在上述滑行试验状态下,根据2种不同的滑行试验方法以及加/减速度的数据分析,可列方程组:
2Tn1–D-f=ma1 (7)
Tn1+T关车–D-f=ma2 (8)
D-f=F (9)
其中,a1,a2为2.2.1节中基准速度下的加/减速度,Tn1为左发动机转速n1对应的推力,T关车为右发动机关车后的推力(发动机关车后存在溢流阻力),D为该滑行状态下的气动阻力,f为该滑行状态下的摩擦阻力,F为该滑行状态下的阻力即气动阻力和摩擦阻力的合力。在相同重量和相同速度时,2种转速状态下的气动阻力和摩擦阻力是相同,可合为一项。即:
2Tn1–F=ma1 (10)
Tn1+T关车–F=ma2 (11)
通过上述方程组,可得
Tn1–T关车=m(a1–a2) (12)
其中,T关车为发动机关车后的推力(即溢流阻力),此时,溢流阻力相对于发动机的推力Tn1仅占据小量,仅约占1%~2%。因此,认为T关车的推力为0,方程(12)化简为
Tn1≈m(a1–a2) (13)
此时,得到了转速n1对应的推力,即Tn1
本实施例的其他部分与上述实施例1-5任一项相同,故不再赘述。
实施例7:
本实施例是在实施例1-6任一个的基础上进行优化,在相同速度下,为了辨识出不同转速下的发动机推力值。可调节n1对应的转速,按照1节(滑行试验方法)和2节(滑行试验数据分析方法)的方法进行反复滑行测试和分析计算,从而,得到相同速度下对应的不同转速推力值,即n-T,见表3。
表3相同速度下的n-T
Figure BDA0002229419870000061
为了辨识出不同速度下的发动机推力值。可改变速度,反复滑行测试,从而,得到不同速度下对应的不同转速推力值,即n-T,见表4。
表4不同速度下的n-T
Figure BDA0002229419870000062
本实施例的其他部分与上述实施例1-6任一项相同,故不再赘述。
实施例8:
一种通过地面滑行试验辨识双发飞机发动机推力的方法,以某型机的多次滑行试验结果为例,机场高度约500m,飞机的重量为2108kg,滑行速度(IAS)为100km/h,进行n1对应最大状态下的多次双发同步和异步的滑行试验,加速度结果见图5。以式(13)进行分析计算,得到的T最大的推力见图6。同时,相同状态下的台架推力也同样给出,可以看出:试验推力多次结果略有波动,与相同状态的台架推力差仅为4%~7%。
本发明克服了发动机台架试验推力数据的差异性、空中飞行试验费用的高昂性,降低首飞的试飞风险。本发明无需加装传感器和复杂的推导,能够直接准确的测量出发动机的实际推力,帮助设计人员修正台架试验数据,掌握该发动机的性能。本发明利用地面滑行试验进行,极大的节约了人力和财力。同时,为后续试飞规划的合理制定提供了数据支撑,保证试飞架次的顺利进行,缩短了新机试飞的周期,有利于产品更快的进入市场和占领市场。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。

Claims (7)

1.一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,忽略发动机转速变化而引起的进气畸变变化量;在重量m、速度V下分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验;通过加/减速运动曲线计算得到加/减速度,进而计算得到不同转速下的推力。
2.根据权利要求1所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,双发同步加/减速的双发转速为n1;所述双发异步加/减速的一发转速为n1,且另一发关车。
3.根据权利要求1或2所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,所述推力的计算公式如下:
2Tn1-D-f=ma1 (7)
Tn1+T关车–D-f=ma2 (8)
D-f=F (9)
其中,a1,a2为滑行过程中的加速度,
Tn1为左发动机转速n1对应的推力,
T关车为右发动机关车后的推力,发动机关车后存在溢流阻力,
D为该滑行状态下的气动阻力,
f为该滑行状态下的摩擦阻力,
F为该滑行状态下的阻力即气动阻力和摩擦阻力的合力;
在相同重量和相同速度时,2种转速状态下的气动阻力和摩擦阻力相同,即:
2Tn1-F=ma1 (10)
Tn1+T关车–F=ma2 (11)
通过方程(7)-(11),可得:
Tn1–T关车=m(a1–a2) (12)
其中,T关车为发动机关车后的推力,即溢流阻力,此时,溢流阻力相对于发动机的推力Tn1仅占据小量,认为T关车的推力为0,方程(12)化简为
Tn1≈m(a1–a2) (13)
此时,得到了转速n1对应的推力,即Tn1
4.根据权利要求3所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,滑行试验中加/减速运动的力学模型如下:
L+N=G (1)
T-D-f=ma (2)
Figure FDA0002229419860000021
Figure FDA0002229419860000022
f=μN (5)
其中,L为飞机的升力,N为飞机轮胎的支撑力,G为飞机的重力,T为两个发动机的总推力,D为飞机的阻力,f为飞机受到的摩擦力,m为飞机的质量,a为飞机的加/减速度,ρ为密度,V为速度,S为参考面积,CL为升力系数,CD为阻力系数。
5.根据权利要求1所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,所述加/减速度计算公式如下:
Figure FDA0002229419860000023
其中,选取速度为V的滑行状态点作为不同发动机状态的基准速度,进而得到基准速度区间的加/减速度。
6.根据权利要求1所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,在重量m、速度V下改变转速分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验,测试得到相同速度下不同转速的推力。
7.根据权利要求1所述的一种通过地面滑行试验辨识双发飞机发动机推力的方法,其特征在于,在重量m下改变速度分别进行双发在不同转速下的同步加/减速以及双发异步的加/减速的滑行试验,测试得到不同速度下不同转速的推力。
CN201910962534.5A 2019-10-11 2019-10-11 一种通过地面滑行试验辨识双发飞机发动机推力的方法 Pending CN110686813A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910962534.5A CN110686813A (zh) 2019-10-11 2019-10-11 一种通过地面滑行试验辨识双发飞机发动机推力的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910962534.5A CN110686813A (zh) 2019-10-11 2019-10-11 一种通过地面滑行试验辨识双发飞机发动机推力的方法

Publications (1)

Publication Number Publication Date
CN110686813A true CN110686813A (zh) 2020-01-14

Family

ID=69112060

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910962534.5A Pending CN110686813A (zh) 2019-10-11 2019-10-11 一种通过地面滑行试验辨识双发飞机发动机推力的方法

Country Status (1)

Country Link
CN (1) CN110686813A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114638181A (zh) * 2022-05-18 2022-06-17 中国航空工业集团公司西安飞机设计研究所 一种校正飞机起飞性能模型的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090320445A1 (en) * 2008-06-30 2009-12-31 Airbus France Method and device for the detection of a thrust dissymmetry of an aircraft during a braking
CN202511930U (zh) * 2012-02-23 2012-10-31 北京理工大学 一种用于测量小型螺旋桨发动机推力和扭矩的装置
EP2930341A1 (en) * 2014-04-11 2015-10-14 Rolls-Royce Corporation Strain gauge and accelerometer measurement for thrust estimation
CN105973516A (zh) * 2015-12-11 2016-09-28 北京强度环境研究所 一种用于识别固体火箭发动机的脉动推力的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090320445A1 (en) * 2008-06-30 2009-12-31 Airbus France Method and device for the detection of a thrust dissymmetry of an aircraft during a braking
CN202511930U (zh) * 2012-02-23 2012-10-31 北京理工大学 一种用于测量小型螺旋桨发动机推力和扭矩的装置
EP2930341A1 (en) * 2014-04-11 2015-10-14 Rolls-Royce Corporation Strain gauge and accelerometer measurement for thrust estimation
CN105973516A (zh) * 2015-12-11 2016-09-28 北京强度环境研究所 一种用于识别固体火箭发动机的脉动推力的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
范卫平 等: "航空器地面滑行状态和发动机推力研究", 《航空计算技术》 *
蔡良才 等: "高原机场飞机起飞滑跑距离计算方法", 《交通运输工程学报》 *
雷晓波 等: "航空发动机推力直接测量飞行试验", 《航空动力学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114638181A (zh) * 2022-05-18 2022-06-17 中国航空工业集团公司西安飞机设计研究所 一种校正飞机起飞性能模型的方法

Similar Documents

Publication Publication Date Title
CN110816874B (zh) 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法
CN108871788B (zh) 一种自动变速器换挡品质测试台架的校验方法
CN108152050B (zh) 一种整车参数校准方法
CN111382522B (zh) 一种基于起飞滑跑数据的航空发动机安装推力评估方法
CN103308325B (zh) 电动汽车驱动系统半实物仿真平台
CN102323206B (zh) 一种车辆滚动阻力系数的测量方法
CN108332937A (zh) 一种风洞连续变攻角测力试验数据修正方法
CN109977448A (zh) 一种基于实测结构载荷的旋翼气动载荷分析方法
CN111709195B (zh) 一种基于盘旋试飞数据巡航耗油特性的修正方法
CN103612772B (zh) 一种飞机刹车系统动态载荷谱联试方法
CN101813554A (zh) 可在同一模型上进行测量的进气道实验装置及工作方法
CN102305715A (zh) 一种用于汽车动力系统测试的动态负载模拟装置及方法
CN103359293A (zh) 飞机制动控制方法和相应的制动控制装置及飞机
CN111767609A (zh) 一种基于试飞数据标准重量下爬升率的修正方法
CN107002578A (zh) 用于确定推进力矩的方法和装置
CN104298804A (zh) 一种飞行载荷设计方法
CN110686813A (zh) 一种通过地面滑行试验辨识双发飞机发动机推力的方法
CN116086756A (zh) 一种前飞状态下直升机全机配平的地面模拟试验方法
CN115655642A (zh) 一种有效评估直升机旋翼气动性能的试验方法
CN109094816B (zh) 一种测试飞机气动升力的方法
CN111717411B (zh) 一种基于试飞数据标准重量下巡航推力增量的修正方法
CN100419402C (zh) 汽车底盘测功机系统
CN112149299B (zh) 用于飞行模拟机数据包开发的试飞任务的优化分析方法
CN109543297A (zh) 一种飞机机身对弹射座椅气动干扰的修正方法
CN102506890A (zh) 一种车辆滑行距离的测量方法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200114

RJ01 Rejection of invention patent application after publication