CN105509946A - 一种辨识飞机升降舵效率的方法 - Google Patents
一种辨识飞机升降舵效率的方法 Download PDFInfo
- Publication number
- CN105509946A CN105509946A CN201510873241.1A CN201510873241A CN105509946A CN 105509946 A CN105509946 A CN 105509946A CN 201510873241 A CN201510873241 A CN 201510873241A CN 105509946 A CN105509946 A CN 105509946A
- Authority
- CN
- China
- Prior art keywords
- test
- aircraft
- moment
- state
- error
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L3/00—Measuring torque, work, mechanical power, or mechanical efficiency, in general
- G01L3/26—Devices for measuring efficiency, i.e. the ratio of power output to power input
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及航空制造领域,涉及航空制造领域的测试领域,具体为一种辨识飞机升降舵效率的方法,其特征在于包括:a.地面滑行试验方法、b.气动参数测定方法和c.误差分析与数据辨识方法,所述地面滑行试验方法包括:a1无动力滑行试验和a2气动参数测定滑行试验,所述气动参数测定方法包括:b1升力系数测试、b2俯仰力矩系数测试和b3升降舵效率测试,所述误差分析与数据辨识方法包括:c1测量误差和c2计算误差,该方法可获得飞机升降舵效率,且简单、数据准确、安全可靠、费用低廉,同时通过测试误差分析方法对结果数据进行分析,解决了无人机首飞前关键气动参数的验证问题,为无人机首飞安全提供保证。
Description
技术领域
本发明涉及航空制造领域,涉及航空制造领域的测试领域,具体为一种辨识飞机升降舵效率的方法。
背景技术
目前,公知的升降舵效率辨识方法是通过进行风洞试验、飞行试验,对试验数据进行一定修正获得。风洞试验是由模型的气动特性间接获得真实飞机的气动特性,包括升降舵效率,这之间要通过相似性、干扰修正等环节,因此要获得真实飞机的气动参数还比较困难,虽然现在对风洞试验数据的洞壁干扰修正和雷诺数修正等方面已有较成熟的方法,但其修正结果的准确性仍然需要进一步提高。飞行试验可以直接获取真实飞机在真实大气环境中的气动特性,包括升降舵效率,是获取真实飞机准确气动参数的重要途径,也是风洞试验与飞行相关性研究的基础,但空中飞行试验费用高昂且风险较大。
发明内容
为了克服现有升降舵效率辨识方法结果精度低、费用高昂和风险较大的不足,本发明提供一种辨识飞机升降舵效率的方法,该方法通过地面滑行试验获得飞机的气动参数,利用升降舵效率测定方法获得飞机的升降舵效率,通过误差分析与数据辨识方法,分析测量误差与计算误差产生机理,可得到合理的滑行试验策略,保证滑行的安全和数据的可靠。
一种辨识飞机升降舵效率的方法,其特征在于包括:
一、地面滑行试验方法
a.无动力滑行试验
无动力滑行试验由牵引车牵引进行,通过无动力滑行试验来检查飞机飞控的部分控制能力(如纠偏),起落架、刹车系统的工作情况和全机各系统的工作匹配情况。
b.气动参数测定滑行试验
气动参数测定滑行试验则是根据制定的具体滑行试验方案从而测得飞机的各项气动参数。气动参数测定滑行试验要求滑行速度尽量大,因为速度越大气动效率越大,风等环境影响的比重越小,测试结果越准确。因此,在地面滑行试验方案是:首先对飞机气动特性(如抬前轮速度)进行分析,然后滑行试验中由低速到高速逐步接近地面滑行的最大速度。
在地面滑行试验中,为了试验安全,一般要避免飞机在滑行过程中前轮抬起,因此,地面滑行试验的滑行终止速度要应尽量接近飞机实际抬前轮速度,同时又要留有足够的余量,这就要求获得准确的飞机加、减速特性。
为了测试飞机升降舵效率,需要偏转舵面和改变飞机重心。对于偏转舵面的滑行试验,要求选择合适的舵面偏度,既保证舵面的气动效率足够大,以便测量,也要保证在滑行过程中飞机不至于抬起前轮,引起危险。
舵面效率测定滑行试验方案如下(具体数据根据飞机不同而调整):
飞机重心位置20.21%bA,分别以升降舵偏度+5°、0°、-2°、-5°,发动机转速95%加速到150km/h,测试升降舵效率;
飞机重心位置30.39%bA,分别以升降舵偏度+5°、0°、-2°、-5°,发动机转速95%加速到150km/h,测试升降舵效率。
对不同重心位置情况分别测定升降舵效率,通过分析重心位置对测试结果的影响,给出测试结果修正量,从而进一步提高升降舵效率测试结果的准确性。
二、气动参数测定方法
对飞机地面滑行过程的受力图如图2:
其中,Fn,Fm分别为前、主起落架支反力,fn、fm分别为前、主起落架摩擦阻力,G为飞机重力,X为飞机气动阻力,Y为飞机气动升力,T为飞机发动机推力。
飞机在地面滑行时,纵向为加速运动,法向则为随起落架支柱伸缩的振动运动,对于大型飞机,此运动的幅值和频率都不大,因此整个过程飞机法向可以看作平衡状态。另外,由于飞机在滑行过程中姿态变化很小,因此,俯仰运动方向也可看作平衡状态。
综上,根据飞机受力情况可知,若已知飞机滑行过程中起落架的载荷和对飞机的力矩,便可以根据平衡方程求出飞机的气动升力和俯仰力矩,从而获得飞机的升降舵效率。
1.升力系数测试
飞机在地面滑跑过程中,根据法向力的平衡关系,有
其中,Y为升力,Nlg为前、主起落架支反力之和,T为发动机推力,G飞机重力。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:
状态2:
若状态1的升力为零,则状态2的升力可由下式计算得到:
在一个架次的滑行过程中,重量的变化主要是油料的消耗。在飞机从滑出到刹车,时间很短,燃油消耗量很小,因此飞机重量的变化较小,而若发动机保持在一定的转速,则推力变化很小,且在飞机开始抬头前,迎角的变化也是小量,因此可以忽略重量和发动机推力法向分量的增量,升力计算则为:
Y=-(Nlg2-Nlg1)
试验测试参数:飞机俯仰角,前起落架行程,主起落架行程。
取状态1为飞机刚滑出的某一状态,飞机速度小,升力近似为零。利用起落架行程载荷曲线,通过对试验中测定的前、主起落架行程插值,计算得到升力以及升力系数。图3为根据某次滑行试验计算得到的升力系数与相同状态下风洞试验结果比较。
2.俯仰力矩系数测试
滑行过程中,由于摩擦力难以估算,因此对主起落架取矩,以避免摩擦力的力矩的计算。对主起的力矩包括:气动力矩,前起支反力力矩,重力力矩,发动机推力力矩,纵向惯性力力矩。
根据力矩平衡:
Mz+Mnlg+MG+Mnx+MT=0
其中,Mz为气动俯仰力矩,Mnlg为前起落架支反力力矩,MG为重力力矩,Mnx为纵向惯性力矩,MT为发动机推力力矩。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:Mz1+Mnlg1+MG1+Mnx1+MT1=0
状态2:Mz2+Mnlg2+MG2+Mnx2+MT2=0
若状态1气动力矩为零,则状态2气动力矩表示为:
Mz2=-[(Mnlg2-Mnlg1)+(MG2-MG1)+(Mnx2-Mnx1)+(MT2-MT1)]
根据上述公式获得的力矩还需变换到参考重心位置,才能与风洞实验数据作比较。变换公式为:
M'z=Mz-Y[lcos(α)+hsin(α)]-X[lsin(α)-hcos(α)]
其中升力根据上一节的方法求得,l为参考重心到主起的水平距离,h为参考重心到主起的垂直距离。由于阻力难以准确求得,且阻力与升力相比小得多,因此采用风洞实验结果。
计算结果及与风洞试验数据比较见图4。
1.升降舵效率测试
俯仰力矩系数表示为:
对于相同重心位置、不同升降舵偏度的两个状态,分别有:
状态1:
状态2:
两个状态的俯仰力矩差量:
若选取的两个状态俯仰角相差不大,则升力系数相差不大,可忽略掉上式右边第一项,则舵面效率为:
不同试验速度下的mz~δz曲线如图5。
三、误差分析与数据辨识方法
1.测量误差
测量误差主要存在于如下几方面:
1)起落架行程载荷曲线测定误差;
2)滑行过程飞机姿态、速度、加速度、起落架行程等测量误差;
对于测量误差,主要体现在如下几方面:
1)传感器测量精度引起的误差;
2)数据采集中噪声引起的误差;
3)测量手段引起的误差。
对于大型无人机,选用的传感器测量精度一般较高,因此,传感器测量误差的影响不大。
飞机在滑行过程中的状态变化是一个慢变化过程,频率不高,而在采集的数据中经常会出现高频噪声,这主要是由于环境(紊流风、跑道面粗糙度等),数据采集、传输过程噪声引起的。因此,在进行试验数据处理时,可以采用低通滤波降低高频噪声引起的误差。如图6和图7为试验俯仰角数据滤波前后比较。
2.计算误差
计算误差主要是由测试方法的相关假设和近似引起的,合理的前提和假设是该方法成立的基础,也是该测试方法误差的主要来源。
计算误差主要体现在如下两方面:
a)平衡假设引起的误差
从图6和图8可见,地面滑行过程中,无人机俯仰角和法向速度都是小量,且变化量不大。因此,平衡状态的假设是合理正确的。
b)小量近似假设引起的误差
为了使计算方法简单,避免难以测量和测量不太准确地量参与计算,通过试验方案设计和计算模型设计,使这些量在模型中成为小量,从而近似忽略,因此也带来了一定计算误差。
本发明的有益效果:
该方法可获得飞机升降舵效率,且简单、数据准确、安全可靠、费用低廉。同时本发明通过地面滑行试验的方法获得直接的飞机气动参数,利用测试气动参数的数学模型及方法计算升降舵效率,通过测试误差分析方法对结果数据进行分析,建立了通过地面滑行测试飞机气动参数的方法,解决了无人机首飞前关键气动参数的验证问题,为无人机首飞安全提供保证;建立了无人机地面滑行气动参数辨识模型,用该模型辨识得到升降舵效率;通过对滑行试验误差机理的分析,得到了合理的滑行试验策略,保证了滑行的安全和数据的可靠。
附图说明:
图1为发明的整体流程图;
图2为飞机滑跑受力图的受力图;
图3为滑行测定的升力系数与风洞试验结果比较图;
图4为测定的俯仰力矩系数与风洞试验结果比较图;
图5为在不同试验速度、升降舵偏度下的俯仰力矩所对应的系数图;
图6为经过计算之后得出的舵效增加百分比列表;
图7为试验俯仰角的滤波之前的数据图;
图8为为试验俯仰角的滤波之后的数据图;
图9为滑行试验法向速度变化图;
图10为通过地面滑行试验辨识飞机升降舵效率方法的结构图。
具体的实施方式:
实施例1:
一种辨识飞机升降舵效率的方法,其特征在于包括:a.地面滑行试验方法、b.气动参数测定方法和c.误差分析与数据辨识方法,所述地面滑行试验方法包括:a1无动力滑行试验和a2气动参数测定滑行试验,所述气动参数测定方法包括:b1升力系数测试、b2俯仰力矩系数测试和b3升降舵效率测试,所述误差分析与数据辨识方法包括:c1测量误差和c2计算误差。
所述a.地面滑行试验方法
a1.无动力滑行试验
无动力滑行试验由牵引车牵引进行,通过无动力滑行试验来检查飞机飞控的部分控制能力(如纠偏),起落架、刹车系统的工作情况和全机各系统的工作匹配情况。
a2.气动参数测定滑行试验
气动参数测定滑行试验则是根据制定的具体滑行试验方案从而测得飞机的各项气动参数;
舵面效率测定滑行试验方案如下(具体数据根据飞机不同而调整):
调整飞机重心位置为30.39%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率;
所述b气动参数测定方法包括:
对飞机地面滑行过程的受力图如图2:
其中,Fn,Fm分别为前、主起落架支反力,fn、fm分别为前、主起落架摩擦阻力,G为飞机重力,X为飞机气动阻力,Y为飞机气动升力,T为飞机发动机推力;
b1.升力系数测试
飞机在地面滑跑过程中,根据法向力的平衡关系,有
其中,Y为升力,Nlg为前、主起落架支反力之和,T为发动机推力,G飞机重力。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:
状态2:
若状态1的升力为零,则状态2的升力可由下式计算得到:
在一个架次的滑行过程中,重量的变化主要是油料的消耗,在飞机从滑出到刹车,时间很短,燃油消耗量很小,因此飞机重量的变化较小,而若发动机保持在一定的转速,则推力变化很小,且在飞机开始抬头前,迎角的变化也是小量,因此可以忽略重量和发动机推力法向分量的增量,升力计算则为:
Y=-(Nlg2-Nlg1)
试验测试参数:飞机俯仰角,前起落架行程,主起落架行程;
取状态1为飞机刚滑出的某一状态,飞机速度小,升力近似为零。利用起落架行程载荷曲线,通过对试验中测定的前、主起落架行程插值,计算得到升力以及升力系数、图3为根据某次滑行试验计算得到的升力系数与相同状态下风洞试验结果比较;
b2.俯仰力矩系数测试
滑行过程中,由于摩擦力难以估算,因此对主起落架取矩,以避免摩擦力的力矩的计算。对主起的力矩包括:气动力矩,前起支反力力矩,重力力矩,发动机推力力矩,纵向惯性力力矩;
根据力矩平衡:
Mz+Mnlg+MG+Mnx+MT=0
其中,Mz为气动俯仰力矩,Mnlg为前起落架支反力力矩,MG为重力力矩,Mnx为纵向惯性力矩,MT为发动机推力力矩。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:Mz1+Mnlg1+MG1+Mnx1+MT1=0
状态2:Mz2+Mnlg2+MG2+Mnx2+MT2=0
若状态1气动力矩为零,则状态2气动力矩表示为:
Mz2=-[(Mnlg2-Mnlg1)+(MG2-MG1)+(Mnx2-Mnx1)+(MT2-MT1)]
根据上述公式获得的力矩还需变换到参考重心位置,才能与风洞实验数据作比较,变换公式为:
M'z=Mz-Y[lcos(α)+hsin(α)]-X[lsin(α)-hcos(α)]
其中升力根据上一节的方法求得,l为参考重心到主起的水平距离,h为参考重心到主起的垂直距离;
计算结果及与风洞试验数据比较见图4;
b3.升降舵效率测试
俯仰力矩系数表示为:
对于相同重心位置、不同升降舵偏度的两个状态,分别有:
状态1:
状态2:
两个状态的俯仰力矩差量:
若选取的两个状态俯仰角相差不大,则升力系数相差不大,可忽略掉上式右边第一项,则舵面效率为:
不同试验速度下的mz~δz曲线如图5;
所述c误差分析与数据辨识方法包括:
c1.测量误差
测量误差主要存在于如下几方面:
1)起落架行程载荷曲线测定误差;
2)滑行过程飞机姿态、速度、加速度、起落架行程等测量误差;
对于测量误差,主要体现在如下几方面:
1)传感器测量精度引起的误差;
2)数据采集中噪声引起的误差;
3)测量手段引起的误差。
对于大型无人机,选用的传感器测量精度一般较高,因此,传感器测量误差的影响不大。
飞机在滑行过程中的状态变化是一个慢变化过程,频率不高,而在采集的数据中经常会出现高频噪声,这主要是由于环境(紊流风、跑道面粗糙度等),数据采集、传输过程噪声引起的。因此,在进行试验数据处理时,可以采用低通滤波降低高频噪声引起的误差。如图6和图7为试验俯仰角数据滤波前后比较。
c2.计算误差
计算误差主要是由测试方法的相关假设和近似引起的,合理的前提和假设是该方法成立的基础,也是该测试方法误差的主要来源。
计算误差主要体现在如下两方面:
c21.平衡假设引起的误差
从图6和图8可见,地面滑行过程中,无人机俯仰角和法向速度都是小量,且变化量不大。因此,平衡状态的假设是合理正确的,
c22.小量近似假设引起的误差
为了使计算方法简单,避免难以测量和测量不太准确地量参与计算,通过试验方案设计和计算模型设计,使这些量在模型中成为小量,从而近似忽略,因此也带来了一定计算误差。
实施例2:
一种辨识飞机升降舵效率的方法,其特征在于包括:a.地面滑行试验方法、b.气动参数测定方法和c.误差分析与数据辨识方法,所述地面滑行试验方法包括:a1无动力滑行试验和a2气动参数测定滑行试验,所述气动参数测定方法包括:b1升力系数测试、b2俯仰力矩系数测试和b3升降舵效率测试,所述误差分析与数据辨识方法包括:c1测量误差和c2计算误差。
所述a.地面滑行试验方法
a1.无动力滑行试验
无动力滑行试验由牵引车牵引进行,通过无动力滑行试验来检查飞机飞控的部分控制能力(如纠偏),起落架、刹车系统的工作情况和全机各系统的工作匹配情况。
a2.气动参数测定滑行试验
气动参数测定滑行试验则是根据制定的具体滑行试验方案从而测得飞机的各项气动参数;
舵面效率测定滑行试验方案如下(具体数据根据飞机不同而调整):
调整飞机重心位置为20.21%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率;
所述b气动参数测定方法包括:
对飞机地面滑行过程的受力图如图2:
其中,Fn,Fm分别为前、主起落架支反力,fn、fm分别为前、主起落架摩擦阻力,G为飞机重力,X为飞机气动阻力,Y为飞机气动升力,T为飞机发动机推力;
b1.升力系数测试
飞机在地面滑跑过程中,根据法向力的平衡关系,有
其中,Y为升力,Nlg为前、主起落架支反力之和,T为发动机推力,G飞机重力。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:
状态2:
若状态1的升力为零,则状态2的升力可由下式计算得到:
在一个架次的滑行过程中,重量的变化主要是油料的消耗,在飞机从滑出到刹车,时间很短,燃油消耗量很小,因此飞机重量的变化较小,而若发动机保持在一定的转速,则推力变化很小,且在飞机开始抬头前,迎角的变化也是小量,因此可以忽略重量和发动机推力法向分量的增量,升力计算则为:
Y=-(Nlg2-Nlg1)
试验测试参数:飞机俯仰角,前起落架行程,主起落架行程;
取状态1为飞机刚滑出的某一状态,飞机速度小,升力近似为零。利用起落架行程载荷曲线,通过对试验中测定的前、主起落架行程插值,计算得到升力以及升力系数、图3为根据某次滑行试验计算得到的升力系数与相同状态下风洞试验结果比较;
b2.俯仰力矩系数测试
滑行过程中,由于摩擦力难以估算,因此对主起落架取矩,以避免摩擦力的力矩的计算。对主起的力矩包括:气动力矩,前起支反力力矩,重力力矩,发动机推力力矩,纵向惯性力力矩;
根据力矩平衡:
Mz+Mnlg+MG+Mnx+MT=0
其中,Mz为气动俯仰力矩,Mnlg为前起落架支反力力矩,MG为重力力矩,Mnx为纵向惯性力矩,MT为发动机推力力矩。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:Mz1+Mnlg1+MG1+Mnx1+MT1=0
状态2:Mz2+Mnlg2+MG2+Mnx2+MT2=0
若状态1气动力矩为零,则状态2气动力矩表示为:
Mz2=-[(Mnlg2-Mnlg1)+(MG2-MG1)+(Mnx2-Mnx1)+(MT2-MT1)]
根据上述公式获得的力矩还需变换到参考重心位置,才能与风洞实验数据作比较,变换公式为:
M'z=Mz-Y[lcos(α)+hsin(α)]-X[lsin(α)-hcos(α)]
其中升力根据上一节的方法求得,l为参考重心到主起的水平距离,h为参考重心到主起的垂直距离;
计算结果及与风洞试验数据比较见图4;
b3.升降舵效率测试
俯仰力矩系数表示为:
对于相同重心位置、不同升降舵偏度的两个状态,分别有:
状态1:
状态2:
两个状态的俯仰力矩差量:
若选取的两个状态俯仰角相差不大,则升力系数相差不大,可忽略掉上式右边第一项,则舵面效率为:
不同试验速度下的mz~δz曲线如图5;
所述c误差分析与数据辨识方法包括:
c1.测量误差
测量误差主要存在于如下几方面:
1)起落架行程载荷曲线测定误差;
2)滑行过程飞机姿态、速度、加速度、起落架行程等测量误差;
对于测量误差,主要体现在如下几方面:
1)传感器测量精度引起的误差;
2)数据采集中噪声引起的误差;
3)测量手段引起的误差。
对于大型无人机,选用的传感器测量精度一般较高,因此,传感器测量误差的影响不大。
飞机在滑行过程中的状态变化是一个慢变化过程,频率不高,而在采集的数据中经常会出现高频噪声,这主要是由于环境(紊流风、跑道面粗糙度等),数据采集、传输过程噪声引起的。因此,在进行试验数据处理时,可以采用低通滤波降低高频噪声引起的误差。如图6和图7为试验俯仰角数据滤波前后比较。
c2.计算误差
计算误差主要是由测试方法的相关假设和近似引起的,合理的前提和假设是该方法成立的基础,也是该测试方法误差的主要来源。
计算误差主要体现在如下两方面:
c21.平衡假设引起的误差
从图6和图8可见,地面滑行过程中,无人机俯仰角和法向速度都是小量,且变化量不大。因此,平衡状态的假设是合理正确的,
c22.小量近似假设引起的误差
为了使计算方法简单,避免难以测量和测量不太准确地量参与计算,通过试验方案设计和计算模型设计,使这些量在模型中成为小量,从而近似忽略,因此也带来了一定计算误差。
实施例3:
一种辨识飞机升降舵效率的方法,其特征在于包括:a.地面滑行试验方法、b.气动参数测定方法和c.误差分析与数据辨识方法,所述地面滑行试验方法包括:a1无动力滑行试验和a2气动参数测定滑行试验,所述气动参数测定方法包括:b1升力系数测试、b2俯仰力矩系数测试和b3升降舵效率测试,所述误差分析与数据辨识方法包括:c1测量误差和c2计算误差。
所述a.地面滑行试验方法
a1.无动力滑行试验
无动力滑行试验由牵引车牵引进行,通过无动力滑行试验来检查飞机飞控的部分控制能力(如纠偏),起落架、刹车系统的工作情况和全机各系统的工作匹配情况。
a2.气动参数测定滑行试验
气动参数测定滑行试验则是根据制定的具体滑行试验方案从而测得飞机的各项气动参数;
舵面效率测定滑行试验方案如下(具体数据根据飞机不同而调整):
调整飞机重心位置为20.21%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率;
调整飞机重心位置为30.39%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率;
对不同重心位置情况分别测定升降舵效率,通过分析重心位置对测试结果的影响,给出测试结果修正量,从而进一步提高升降舵效率测试结果的准确性;
所述b气动参数测定方法包括:
对飞机地面滑行过程的受力图如图2:
其中,Fn,Fm分别为前、主起落架支反力,fn、fm分别为前、主起落架摩擦阻力,G为飞机重力,X为飞机气动阻力,Y为飞机气动升力,T为飞机发动机推力;
b1.升力系数测试
飞机在地面滑跑过程中,根据法向力的平衡关系,有
其中,Y为升力,Nlg为前、主起落架支反力之和,T为发动机推力,G飞机重力。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:
状态2:
若状态1的升力为零,则状态2的升力可由下式计算得到:
在一个架次的滑行过程中,重量的变化主要是油料的消耗,在飞机从滑出到刹车,时间很短,燃油消耗量很小,因此飞机重量的变化较小,而若发动机保持在一定的转速,则推力变化很小,且在飞机开始抬头前,迎角的变化也是小量,因此可以忽略重量和发动机推力法向分量的增量,升力计算则为:
Y=-(Nlg2-Nlg1)
试验测试参数:飞机俯仰角,前起落架行程,主起落架行程;
取状态1为飞机刚滑出的某一状态,飞机速度小,升力近似为零。利用起落架行程载荷曲线,通过对试验中测定的前、主起落架行程插值,计算得到升力以及升力系数、图3为根据某次滑行试验计算得到的升力系数与相同状态下风洞试验结果比较;
b2.俯仰力矩系数测试
滑行过程中,由于摩擦力难以估算,因此对主起落架取矩,以避免摩擦力的力矩的计算。对主起的力矩包括:气动力矩,前起支反力力矩,重力力矩,发动机推力力矩,纵向惯性力力矩;
根据力矩平衡:
Mz+Mnlg+MG+Mnx+MT=0
其中,Mz为气动俯仰力矩,Mnlg为前起落架支反力力矩,MG为重力力矩,Mnx为纵向惯性力矩,MT为发动机推力力矩。
因此对于滑行中两个不同的状态1和2,分别有:
状态1:Mz1+Mnlg1+MG1+Mnx1+MT1=0
状态2:Mz2+Mnlg2+MG2+Mnx2+MT2=0
若状态1气动力矩为零,则状态2气动力矩表示为:
Mz2=-[(Mnlg2-Mnlg1)+(MG2-MG1)+(Mnx2-Mnx1)+(MT2-MT1)]
根据上述公式获得的力矩还需变换到参考重心位置,才能与风洞实验数据作比较,变换公式为:
M'z=Mz-Y[lcos(α)+hsin(α)]-X[lsin(α)-hcos(α)]
其中升力根据上一节的方法求得,l为参考重心到主起的水平距离,h为参考重心到主起的垂直距离;
计算结果及与风洞试验数据比较见图4;
b3.升降舵效率测试
俯仰力矩系数表示为:
对于相同重心位置、不同升降舵偏度的两个状态,分别有:
状态1:
状态2:
两个状态的俯仰力矩差量:
若选取的两个状态俯仰角相差不大,则升力系数相差不大,可忽略掉上式右边第一项,则舵面效率为:
不同试验速度下的mz~δz曲线如图5;
所述c误差分析与数据辨识方法包括:
c1.测量误差
测量误差主要存在于如下几方面:
1)起落架行程载荷曲线测定误差;
2)滑行过程飞机姿态、速度、加速度、起落架行程等测量误差;
对于测量误差,主要体现在如下几方面:
1)传感器测量精度引起的误差;
2)数据采集中噪声引起的误差;
3)测量手段引起的误差。
对于大型无人机,选用的传感器测量精度一般较高,因此,传感器测量误差的影响不大。
飞机在滑行过程中的状态变化是一个慢变化过程,频率不高,而在采集的数据中经常会出现高频噪声,这主要是由于环境(紊流风、跑道面粗糙度等),数据采集、传输过程噪声引起的。因此,在进行试验数据处理时,可以采用低通滤波降低高频噪声引起的误差。如图6和图7为试验俯仰角数据滤波前后比较。
c2.计算误差
计算误差主要是由测试方法的相关假设和近似引起的,合理的前提和假设是该方法成立的基础,也是该测试方法误差的主要来源。
计算误差主要体现在如下两方面:
c21.平衡假设引起的误差
从图6和图8可见,地面滑行过程中,无人机俯仰角和法向速度都是小量,且变化量不大。因此,平衡状态的假设是合理正确的,
c22.小量近似假设引起的误差
为了使计算方法简单,避免难以测量和测量不太准确地量参与计算,通过试验方案设计和计算模型设计,使这些量在模型中成为小量,从而近似忽略,因此也带来了一定计算误差。
Claims (4)
1.一种辨识飞机升降舵效率的方法,其特征在于包括:a.地面滑行试验方法、b.气动参数测定方法和c.误差分析与数据辨识方法,所述地面滑行试验方法包括:a1无动力滑行试验和a2气动参数测定滑行试验,所述气动参数测定方法包括:b1升力系数测试、b2俯仰力矩系数测试和b3升降舵效率测试,所述误差分析与数据辨识方法包括:c1测量误差和c2计算误差。
2.根据权利要求1所述的一种辨识飞机升降舵效率的方法,其特征在于:所述a.地面滑行试验方法
a1.无动力滑行试验
无动力滑行试验由牵引车牵引进行,通过无动力滑行试验来检查飞机飞控的部分控制能力(如纠偏),起落架、刹车系统的工作情况和全机各系统的工作匹配情况;
a2.气动参数测定滑行试验
气动参数测定滑行试验则是根据制定的具体滑行试验方案从而测得飞机的各项气动参数;
舵面效率测定滑行试验方案如下:
调整飞机重心位置为20.21%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率;
调整飞机重心位置为30.39%bA,分别以升降舵偏度+5°、0°、-2°、-5°,将发动机转速95%加速到150km/h,测试升降舵效率。
3.根据权利要求1所述的一种辨识飞机升降舵效率的方法,其特征在于:所述b气动参数测定方法包括:
对飞机地面滑行过程的受力分析:
其中,Fn,Fm分别为前、主起落架支反力,fn、fm分别为前、主起落架摩擦阻力,G为飞机重力,X为飞机气动阻力,Y为飞机气动升力,T为飞机发动机推力;
b1.升力系数测试
飞机在地面滑跑过程中,根据法向力的平衡关系,有
其中,Y为升力,Nlg为前、主起落架支反力之和,T为发动机推力,G飞机重力;
因此对于滑行中两个不同的状态1和2,分别有:
状态1:
状态2:
若状态1的升力为零,则状态2的升力可由下式计算得到:
在一个架次的滑行过程中,重量的变化主要是油料的消耗,在飞机从滑出到刹车,时间很短,燃油消耗量很小,因此飞机重量的变化较小,而若发动机保持在一定的转速,则推力变化很小,且在飞机开始抬头前,迎角的变化也是小量,因此可以忽略重量和发动机推力法向分量的增量,升力计算则为:
Y=-(Nlg2-Nlg1)
试验测试参数:飞机俯仰角,前起落架行程,主起落架行程;
取状态1为飞机刚滑出的某一状态,飞机速度小,升力近似为零,利用起落架行程载荷曲线,通过对试验中测定的前、主起落架行程插值,计算得到升力以及升力系数;
b2.俯仰力矩系数测试
滑行过程中,由于摩擦力难以估算,因此对主起落架取矩,以避免摩擦力的力矩的计算,对主起的力矩包括:气动力矩,前起支反力力矩,重力力矩,发动机推力力矩,纵向惯性力力矩;
根据力矩平衡:
Mz+Mnlg+MG+Mnx+MT=0
其中,Mz为气动俯仰力矩,Mnlg为前起落架支反力力矩,MG为重力力矩,Mnx为纵向惯性力矩,MT为发动机推力力矩,
因此对于滑行中两个不同的状态1和2,分别有:
状态1:Mz1+Mnlg1+MG1+Mnx1+MT1=0
状态2:Mz2+Mnlg2+MG2+Mnx2+MT2=0
若状态1气动力矩为零,则状态2气动力矩表示为:
Mz2=-[(Mnlg2-Mnlg1)+(MG2-MG1)+(Mnx2-Mnx1)+(MT2-MT1)]
根据上述公式获得的力矩还需变换到参考重心位置,才能与风洞实验数据作比较,变换公式为:
M'z=Mz-Y[lcos(α)+hsin(α)]-X[lsin(α)-hcos(α)]
其中升力根据上一节的方法求得,l为参考重心到主起的水平距离,h为参考重心到主起的垂直距离;
b3.升降舵效率测试
俯仰力矩系数表示为:
对于相同重心位置、不同升降舵偏度的两个状态,分别有:
状态1:
状态2:
两个状态的俯仰力矩差量:
若选取的两个状态俯仰角相差不大,则升力系数相差不大,可忽略掉上式右边第一项,则舵面效率为:
4.根据权利要求1所述的一种辨识飞机升降舵效率的方法,其特征在于:所述c误差分析与数据辨识方法包括:
c1.测量误差
测量误差主要存在于如下几方面:
1)起落架行程载荷曲线测定误差;
2)滑行过程飞机姿态、速度、加速度、起落架行程等测量误差;
c2.计算误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510873241.1A CN105509946A (zh) | 2015-12-03 | 2015-12-03 | 一种辨识飞机升降舵效率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510873241.1A CN105509946A (zh) | 2015-12-03 | 2015-12-03 | 一种辨识飞机升降舵效率的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105509946A true CN105509946A (zh) | 2016-04-20 |
Family
ID=55718091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510873241.1A Pending CN105509946A (zh) | 2015-12-03 | 2015-12-03 | 一种辨识飞机升降舵效率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105509946A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106840572A (zh) * | 2016-12-19 | 2017-06-13 | 中国航天空气动力技术研究院 | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 |
CN107066652A (zh) * | 2016-12-15 | 2017-08-18 | 中国航空工业集团公司西安飞机设计研究所 | 一种估算尾撬支反力的方法 |
CN107145646A (zh) * | 2017-04-19 | 2017-09-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种计算飞机尾撬支反力的方法 |
CN109506831A (zh) * | 2018-11-14 | 2019-03-22 | 北京航天发射技术研究所 | 一种火箭起飞底阻检测装置 |
CN110816874A (zh) * | 2019-10-11 | 2020-02-21 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN111241625A (zh) * | 2019-10-17 | 2020-06-05 | 成都飞机工业(集团)有限责任公司 | 飞机发动机特性辨识及平衡极曲线辨识的试飞方法 |
CN111914345A (zh) * | 2020-06-15 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 一种基于参数辨识的飞机剩余推力等效试飞模型 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3988839A (en) * | 1975-07-03 | 1976-11-02 | Coomber Allen J | Model aircraft incidence gauge |
JP2004109067A (ja) * | 2002-09-20 | 2004-04-08 | National Aerospace Laboratory Of Japan | 遷音速フラッタ停止装置 |
CN101267987A (zh) * | 2005-09-22 | 2008-09-17 | 法国空中巴士公司 | 用于在航空器上执行至少一个飞行试验的方法和装置及其应用 |
CN102305699A (zh) * | 2011-05-19 | 2012-01-04 | 北京航空航天大学 | 自由飞模型风洞实验系统 |
CN104101477A (zh) * | 2014-07-14 | 2014-10-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种低速风洞测力试验模型自动舵机 |
CN104344882A (zh) * | 2013-07-24 | 2015-02-11 | 中国国际航空股份有限公司 | 一种飞机抖动检测系统及方法 |
-
2015
- 2015-12-03 CN CN201510873241.1A patent/CN105509946A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3988839A (en) * | 1975-07-03 | 1976-11-02 | Coomber Allen J | Model aircraft incidence gauge |
JP2004109067A (ja) * | 2002-09-20 | 2004-04-08 | National Aerospace Laboratory Of Japan | 遷音速フラッタ停止装置 |
CN101267987A (zh) * | 2005-09-22 | 2008-09-17 | 法国空中巴士公司 | 用于在航空器上执行至少一个飞行试验的方法和装置及其应用 |
CN102305699A (zh) * | 2011-05-19 | 2012-01-04 | 北京航空航天大学 | 自由飞模型风洞实验系统 |
CN104344882A (zh) * | 2013-07-24 | 2015-02-11 | 中国国际航空股份有限公司 | 一种飞机抖动检测系统及方法 |
CN104101477A (zh) * | 2014-07-14 | 2014-10-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种低速风洞测力试验模型自动舵机 |
Non-Patent Citations (2)
Title |
---|
晏传银 等: "基于试飞数据的气动特性分析校核", 《山东工业技术》 * |
汪清: "不稳定飞机气动参数辨识的一种实用方法", 《飞行力学》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066652A (zh) * | 2016-12-15 | 2017-08-18 | 中国航空工业集团公司西安飞机设计研究所 | 一种估算尾撬支反力的方法 |
CN107066652B (zh) * | 2016-12-15 | 2020-09-22 | 中国航空工业集团公司西安飞机设计研究所 | 一种估算尾撬支反力的方法 |
CN106840572B (zh) * | 2016-12-19 | 2019-05-24 | 中国航天空气动力技术研究院 | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 |
CN106840572A (zh) * | 2016-12-19 | 2017-06-13 | 中国航天空气动力技术研究院 | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 |
CN107145646B (zh) * | 2017-04-19 | 2021-07-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种计算飞机尾撬支反力的方法 |
CN107145646A (zh) * | 2017-04-19 | 2017-09-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种计算飞机尾撬支反力的方法 |
CN109506831A (zh) * | 2018-11-14 | 2019-03-22 | 北京航天发射技术研究所 | 一种火箭起飞底阻检测装置 |
CN110816874A (zh) * | 2019-10-11 | 2020-02-21 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN110816874B (zh) * | 2019-10-11 | 2022-04-08 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN111241625A (zh) * | 2019-10-17 | 2020-06-05 | 成都飞机工业(集团)有限责任公司 | 飞机发动机特性辨识及平衡极曲线辨识的试飞方法 |
CN111241625B (zh) * | 2019-10-17 | 2022-01-07 | 成都飞机工业(集团)有限责任公司 | 飞机发动机特性辨识及平衡极曲线辨识的试飞方法 |
CN111914345B (zh) * | 2020-06-15 | 2021-06-08 | 成都飞机工业(集团)有限责任公司 | 一种基于参数辨识的飞机剩余推力等效试飞方法 |
CN111914345A (zh) * | 2020-06-15 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 一种基于参数辨识的飞机剩余推力等效试飞模型 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105509946A (zh) | 一种辨识飞机升降舵效率的方法 | |
CN105373647B (zh) | 一种通过地面滑跑试验辨识气动焦点的方法 | |
US7075457B1 (en) | Energy index for aircraft maneuvers | |
CN103994748B (zh) | 一种采用飞行和风洞试验数据估计无人机配平迎角的方法 | |
US8442705B2 (en) | Method and device for determining aerodynamic characteristics of an aircraft | |
CN109710961B (zh) | 一种基于gps数据的高空无人机升限数据处理方法 | |
CN108132134A (zh) | 基于风洞自由飞试验的气动导数辨识方法和系统 | |
CN108153323B (zh) | 一种高空无人飞行器高精度再入制导方法 | |
CN102680201A (zh) | 基于视频测量的抖振风洞试验方法 | |
CN103612772B (zh) | 一种飞机刹车系统动态载荷谱联试方法 | |
Larrabee et al. | Wind field estimation in UAV formation flight | |
CN104102127B (zh) | 一种机载气动参数辨识系统 | |
CN103995529A (zh) | 一种旋翼飞行器机械安装故障在线诊断方法 | |
CN104504255A (zh) | 一种螺旋翼升力和阻力力矩的确定方法 | |
CN103020456B (zh) | 一种飞行器气动参数辨识误差区间估算方法 | |
CA2745315A1 (en) | Flow determination method | |
CN104462022A (zh) | 飞行器动力学系统参数可辨识性分析方法 | |
CN109094816B (zh) | 一种测试飞机气动升力的方法 | |
CN106372307B (zh) | 一种基于气动模型的民机气流角估计方法 | |
Preisighe Viana | Time-domain system identification of rigid-body multipoint loads model | |
CN103869811A (zh) | 一种四旋翼飞行器的遥控和伺服信号的监控方法 | |
CN107967382B (zh) | 一种吸气式飞行器耦合影响程度评价方法 | |
Ranaudo et al. | Effects of horizontal tail ice on longitudinal aerodynamic derivatives | |
CN106508018B (zh) | 基于飞机运动方程的飞行数据“野值”检测与剔除方法 | |
CN111859540B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160420 |