CN1828253A - 一种基于多流速传感器的飞行器飞行参数测量方法 - Google Patents

一种基于多流速传感器的飞行器飞行参数测量方法 Download PDF

Info

Publication number
CN1828253A
CN1828253A CN 200610065681 CN200610065681A CN1828253A CN 1828253 A CN1828253 A CN 1828253A CN 200610065681 CN200610065681 CN 200610065681 CN 200610065681 A CN200610065681 A CN 200610065681A CN 1828253 A CN1828253 A CN 1828253A
Authority
CN
China
Prior art keywords
aircraft
angle
speed
attack
sensor
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
Application number
CN 200610065681
Other languages
English (en)
Other versions
CN100491949C (zh
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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CNB2006100656815A priority Critical patent/CN100491949C/zh
Publication of CN1828253A publication Critical patent/CN1828253A/zh
Application granted granted Critical
Publication of CN100491949C publication Critical patent/CN100491949C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明公开了一种用多个流速传感器组合测量飞行器飞行参数的方法。该方法,包括如下步骤:1)在飞行器机翼前缘上下左右对称位置布置至少一组流速传感器,流速传感器以四个为一组;2)对于一组流速传感器,四个流速传感器测量得到的流速值分别为:s1、s2、s3和s4,飞行器参数按照如下公式计算:空速:u1=1/4(s1+s2+s3+s4),迎角:α1=k1(s1+s2-s3-s4-k2β1),侧滑角:β1=k3(s1-s2-s3+s4-k4α1),滚转角速率:ω1=k5(s1-s2+s3-s4);3)对于一组以上流速传感器,对每组传感器测得的数值取平均值。采用本发明可增加飞行器的有效载荷。

Description

一种基于多流速传感器的飞行器飞行参数测量方法
技术领域
本发明涉及飞行器的飞行参数如空速、迎角、侧滑角及滚转角速率等的测量方法,特别涉及用多个流速传感器组合测量上述信息量的方法。
背景技术
空速、迎角、侧滑角、滚转角速率等飞行参数的准确测量对于飞行器的有效控制和操纵有着重要的作用,通常,已知技术对上述信息量的测量是独立的。空速测量一般用空速管,即皮托管,通过测量空速管处的总压与静压的压差,间接得出空速。迎角与侧滑角的测量方法类似,主要有两类典型传感器,即风标式迎角/侧滑角传感器和差压式迎角/侧滑角传感器,如文献1:“戚宗旭,迎角测量方法讨论,飞行试验,1992年第一期,33-37页”中公开的技术,其中,风标式传感器由对称剖面的翼型叶片(即风标)、转轴、角度变换器、配重等部分组成,当叶片与远前方气流形成一定角度时,气动力的作用将使叶片向气流方向转动,直到偏角趋于零,由此测出飞行器的迎角或侧滑角。差压式传感器通常由差压管和开口膜盒式压力传感器组成,差压管与空速管类似,其端头孔感受总压,侧面的小孔感受静压,在端头与差压管轴线成一定角度有两个对称孔,在一定安装方式下,通过测这两个孔的压差便可求得相应迎角或侧滑角。滚转角速率一般用速率陀螺直接测量,比如采用美国Analog Device公司的ADXRS系列速率陀螺芯片,文献2:“曾庆华,张为华,iMEMS速率陀螺芯片在MAV飞行控制系统中应用研究,测控技术,2004年23卷2期,68-70页”中公开的技术。
上述方法中,四个飞行参数的测量相对较独立,传感器体积和重量都比较大,而且像风标、空速管等器件若受到较大冲击则很容易损坏,因而适用性较差。
由于飞行器的上述飞行参数与其周围的流场情况有直接的联系,因此,人们一直希望有改进的飞行器飞行参数测量方法,比如通过测量飞行器机翼表面流场来推算相关信息量。随着新材料、新工艺的不断发展,流速传感器的微型化、低功耗已成为现实,采用微型流速传感器可以实现飞行参数的测量。
发明内容
本发明的目的是克服现有技术的不足,提供一种基于多流速传感器的飞行器飞行参数的测量方法。
为了达到上述目的,本发明采取如下技术方案:
一种基于多流速传感器的飞行器飞行参数的测量方法,包括如下步骤:
1)在飞行器机翼前缘布置至少一组流速传感器,流速传感器以四个为一组;
2)对于一组流速传感器,四个流速传感器测量得到的流速值分别为:s1、s2、s3和s4,飞行器参数按照如下公式计算:
飞行器空速: u 1 = 1 4 ( s 1 + s 2 + s 3 + s 4 )
飞行器迎角:α1=k1(s1+s2-s3-s4-k2β1)
飞行器侧滑角:β1=k3(s1-s2-s3+s4-k4α1)
飞行器滚转角速率:ω1=k5(s1-s2+s3-s4)
式中,k1、k2、k3、k4、k5为空速u1的函数,由预先标定确定。
3)对于组数为n的一组以上流速传感器,对每组传感器测得的数值取平均值,飞行器参数按照如下公式计算:
飞行器空速: u ‾ = 1 n Σ i = 1 n u i
飞行器迎角: α ‾ = 1 n Σ i = 1 n α i
飞行器侧滑角: β ‾ = 1 n Σ i = 1 n β i
飞行器滚转角速率: ω ‾ = 1 n Σ i = 1 n ω i
在上述技术方案中,进一步地,如图1所示,所述步骤1)中每组的四个传感器分别布置在机翼前缘上、下、左、右四个对称位置。
在上述技术方案中,进一步地,所述步骤1)中根据飞行器的负载情况、功耗情况及所需测量的飞行参数的精确程度确定所需要安装的流速传感器的组数。当飞行器负载能力和所能提供的功耗都充裕的情况下,多组传感器不会影响飞行器的飞行状态,可在机翼前缘布置多组传感器,以使测量更精确。本领域技术人员知道,每组传感器之间的距离视飞行器翼型及所使用的传感器而定,在流速传感器之间互不干扰的前提下。
在上述技术方案中,进一步地,所述步骤1)多组流速传感器之间是互不干扰的。
在上述技术方案中,进一步地,所述流速传感器包括压力流速传感器、光纤流速传感器和热式流速传感器等流速传感器。
在上述技术方案中,进一步地,k1、k2、k3、k4、k5为空速u1的函数,空速可由四个传感器测得的流速值直接得到;迎角和侧滑角之间有耦合关系,根据步骤2)和步骤3)中的表达式进行联合标定;滚转角速率按照步骤2)和步骤3)中的表达式单独标定。标定参数k1、k2、k3、k4、k5具体包括如下步骤:
(1)在一定空速u下,改变迎角和侧滑角,由四个传感器分别得到不同迎角αi(i=1,2…n)和侧滑角βj(j=1,2…m)组合下的相应位置的流速值s1ij,s2ij,s3ij,s4ij,并记:
αtij=s1ij+s2ij-s3ij-s4ij
βtij=s1ij-s2ij-s3ij+s4ij
然后,根据上述数据,采用一般的参数估计方法,如最小二乘法、极大似然估计法等,估计参数k1、k2、k3、k4,由此,可得到迎角和侧滑角的表达式如下:
飞行器迎角:α=k1(s1+s2-s3-s4-k2β)
飞行器侧滑角:β=k3(s1-s2-s3+s4-k4α)
(2)在同样空速下,对滚转角速率进行标定,不同滚转角速率ωg(g=1,2...l)下四个传感器测得的相应位置流速分别为s1g,s2g,s3g,s4g,记:
ωtg=s1g-s2g+s3g-s4g
同样,采用一般的参数估计方法估计参数k5,由此,得到滚转角速率的表达式如下:
飞行器滚转角速率:ω=k5(s1-s2+s3-s4)
(3)在其它空速下,做同样标定,得到不同空速下相关的系数k,飞行器实际飞行时根据测得的空速选取相应系数k,用以求解当时的其它三个飞行参数;或者,预先标定拟和各系数k关于空速u的函数关系式,用以求解飞行器飞行时的其它三个飞行参数。
本发明的基于多流速传感器的飞行器飞行参数测量方法主要包括流速传感器在飞行器机翼前缘的布置方式和通过多个流速传感器的测量值推算飞行器空速、迎角、侧滑角及滚转角速率的方法。
与现有技术相比,本发明具有如下有益效果:
1、采用组合测量方式测量飞行器的四个飞行参数,相对于传统各飞行参数独立测量方法,系统更简洁;而且,可选用小型化的流速传感器,从而减少了测量系统的质量,增加了飞行器的有效载荷。
2、由于流速传感器形式多样,使用时可根据飞行器的实际情况进行选择。特别是对于微型飞行器,由于其负载能力极为有限,而且在着陆时往往与地面发生冲撞,所以传统的相关传感器很难应用于其上,而本发明则克服了这个缺陷。
附图说明
图1表示本发明飞行器飞行参数的测量方法使用1组(4个)流速传感器测量飞行器飞行参数的立体示意图;
图2表示本发明一实施例中使用1组(4个)流速传感器测量飞行器飞行参数的俯视示意图;
图3表示本发明一实施例采用的恒压式热式流速传感器的电路原理图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细描述:
参照图2,图中所示为一微型飞行器的机翼20,机翼20左右关于轴线00’对称,翼展300mm,蒙皮厚0.1mm。
四个流速传感器,为第一流速传感器1、第二流速传感器2、第三流速传感器3、第四流速传感器4分别按顺时针表贴在飞行器机翼前缘左右对称位置的上下表面。从正视飞行器机翼前缘的方向看,这4个传感器按顺时针方向依次为1、2、3、4。本实施例中选用热式流速传感器,其敏感元件为由陶瓷封装的薄膜铂热敏电阻,在市场上买到。铂热敏电阻本身性能稳定,从而使整个传感器网络有稳定的性能。铂薄膜通过激光喷溅在陶瓷表层,因而元件具有优良的防振和防冲击性能,比较适合用于飞行器。
铂热敏电阻通电发热后,其阻值会随着流经传感器的空速的改变而改变,因而可以通过测量热敏电阻的阻值来感受空气流速的变化。每个热式流速传感器均采用恒压式工作方式,其电路原理图如附图3所示。接点5处接恒压源(图中未示出),数字6表示定值电阻,数字7表示第一铂热敏电阻,作流速敏感元件,数字8为第二铂热敏电阻,放置在飞行器内部,对第一铂热敏电阻7起环境温度补偿作用,数字9表示调零电阻。第一铂热敏电阻7的工作电阻Rw和其所在位置的空气流速v之间有如下关系:
U 2 ( R c + R w ) 2 R w = ( R w - R g ) ( A + B v )
式中,Rc——定值电阻6的阻值;
      U——电桥桥顶电压,为定值;
      Rg——环境温度下铂热敏电阻7的阻值;
      A,B——常数,由预先标定确定。
式中,常数A、B的标定在零侧滑角、零迎角及零滚转角速率下进行。测量不同空气流速下传感器的输出,用得到的数据,采用最小二乘法估计参数A、B。
由第一、第二、第三和第四流速传感器测得的相应位置的流速值s1、s2、s3和s4,可直接得到飞行器的空速,即:
飞行器空速: u = 1 4 ( s 1 + s 2 + s 3 + s 4 )
由于飞行器空速对迎角、侧滑角和滚转角速率的测量值有影响,所以需要在飞行器实际飞行时可能的空速范围内对迎角、侧滑角和滚转角速率进行标定,标定时还要考虑迎角与侧滑角之间的耦合关系。本实施例中,迎角的标定范围为0~8°,侧滑角的标定范围为-10°~10°,均比较小,在标定过程中线性拟和相关曲线。在一定空速u下,改变迎角和侧滑角,则由四个传感器分别得到不同迎角αi(i=1,2…n)和侧滑角βj(j=1,2…m)组合下的相应位置的流速值s1ij,s2ij,s3ij,s4ij,并记:
αtij=s1ij+s2ij-s3ij-s4ij
βtij=s1ij-s2ij-s3ij+s4ij
然后,根据上述数据,采用最小二乘法估计参数k1、k2、k3、k4,由此,得到迎角和侧滑角的表达式如下:
飞行器迎角:α=k1(s1+s2-s3-s4-k2β)
飞行器侧滑角:β=k3(s1-s2-s3+s4-k4α)
在同样空速下,对滚转角速率进行标定,不同滚转角速率ωg(g=1,2...l)下四个传感器测得的相应位置流速分别为s1g,s2g,s3g,s4g,记:
ωtg=s1g-s2g+s3g-s4g
同样,采用最小二乘法估计参数k5,由此,得到滚转角速率的表达式如下:
飞行器滚转角速率:ω=k5(s1-s2+s3-s4)
在其它空速下,做同样标定,得到不同空速下相关的系数k,飞行器实际飞行时根据测得的空速选取相应系数k,用以求解当时的其它三个飞行参数。或者,也可预先标定拟和各系数k关于空速u的函数关系式,用以求解飞行器飞行时的其它三个飞行参数。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1、一种基于多流速传感器的飞行器飞行参数的测量方法,包括如下步骤:
1)在飞行器机翼前缘布置至少一组流速传感器,流速传感器以四个为一组;
2)对于一组流速传感器,四个流速传感器测量得到的流速值分别为:s1、s2、s3和s4,飞行器参数按照如下公式计算:
飞行器空速: u 1 = 1 4 ( s 1 + s 2 + s 3 + s 4 )
飞行器迎角:α1=k1(s1+s2-s3-s4-k2β1)
飞行器侧滑角:β1=k3(s1-s2-s3+s4-k4α1)
飞行器滚转角速率:ω1=k5(s1-s2+s3-s4)
式中,k1、k2、k3、k4、k5为空速u1的函数,由预先标定确定;
3)对于组数为n的一组以上流速传感器,对每组传感器测得的数值取平均值,飞行器参数按照如下公式计算:
飞行器空速: u ‾ = 1 n Σ i = 1 n u i
飞行器迎角: α ‾ = 1 n Σ i = 1 n α i
飞行器侧滑角: β ‾ = 1 n Σ i = 1 n β i
飞行器滚转角速率: ω ‾ = 1 n Σ i = 1 n ω i
2、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述步骤1)中每组的四个传感器分别布置在机翼前缘上、下、左、右四个对称位置。
3、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述步骤1)中在机翼前缘位置布置多组传感器,且流速传感器之间互不干扰。
4、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述k1、k2、k3、k4、k5为空速u1的函数,空速由四个传感器测得的流速值直接得到;迎角和侧滑角之间有耦合关系,根据步骤2)和步骤3)中的公式进行联合标定;滚转角速率按照步骤2)和步骤3)中的公式单独标定。
5、根据权利要求4所述基于多流速传感器的飞行器飞行参数的测量方法,其特征在于,标定参数k1、k2、k3、k4、k5具体包括如下步骤:
(1)在一定空速u下,改变迎角和侧滑角,由四个传感器分别得到不同迎角αi(i=1,2…n)和侧滑角βj(j=1,2…m)组合下的相应位置的流速值s1ij,s2ij,s3ij,s4ij,并记:
αtij=s1ij+s2ij-s3ij-s4ij
βtij=s1ij-s2ij-s3ij+s4ij
采用参数估计方法估计参数k1、k2、k3、k4,由此,得到迎角和侧滑角的表达式如下:
飞行器迎角:α=k1(s1+s2-s3-s4-k2β)
飞行器侧滑角:β=k3(s1-s2-s3+s4-k4α)
(2)在同样空速下,对滚转角速率进行标定,不同滚转角速率ωg(g=1,2…l)下四个传感器测得的相应位置流速分别为s1g,s2g,s3g,s4g,记:
ωtg=s1g-s2g+s3g-s4g
采用参数估计方法估计参数k5,由此,得到滚转角速率的表达式如下:
飞行器滚转角速率:ω=k5(s1-s2+s3-s4)
(3)在其它空速下,重复步骤(1)-(2)做同样标定,得到不同空速下相关的系数k,飞行器实际飞行时根据测得的空速选取相应系数k;或者,预先标定拟和各系数k关于空速u的函数关系式,用以求解飞行器飞行时的其它三个飞行参数。
6、根据权利要求5所述基于多流速传感器的飞行器飞行参数的测量方法,其特征在于,所述参数估计方法包括最小二乘法、极大似然估计法。
CNB2006100656815A 2006-03-21 2006-03-21 一种基于多流速传感器的飞行器飞行参数测量方法 Expired - Fee Related CN100491949C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100656815A CN100491949C (zh) 2006-03-21 2006-03-21 一种基于多流速传感器的飞行器飞行参数测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100656815A CN100491949C (zh) 2006-03-21 2006-03-21 一种基于多流速传感器的飞行器飞行参数测量方法

Publications (2)

Publication Number Publication Date
CN1828253A true CN1828253A (zh) 2006-09-06
CN100491949C CN100491949C (zh) 2009-05-27

Family

ID=36946750

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100656815A Expired - Fee Related CN100491949C (zh) 2006-03-21 2006-03-21 一种基于多流速传感器的飞行器飞行参数测量方法

Country Status (1)

Country Link
CN (1) CN100491949C (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101750514B (zh) * 2010-01-27 2011-12-28 中国科学院光电研究院 高空驻空飞行器风速风向实时测量方法与装置
CN102323440A (zh) * 2011-06-03 2012-01-18 清华大学 直升机二轴式空速检测系统
CN101316478B (zh) * 2007-05-29 2012-02-08 清华大学 微型飞行器
CN102662078A (zh) * 2012-05-07 2012-09-12 北京航空航天大学 一种直升机多轴式全向空速传感器
CN106372307A (zh) * 2016-08-30 2017-02-01 中国航空工业集团公司西安飞行自动控制研究所 一种基于气动模型的民机气流角估计方法
CN108139425A (zh) * 2015-10-28 2018-06-08 欧姆龙株式会社 空速测量系统
CN109781107A (zh) * 2017-11-15 2019-05-21 北京自动化控制设备研究所 一种低精度惯导滚动角确定方法
CN113654707A (zh) * 2021-08-12 2021-11-16 北京航空航天大学 基于压力和流速的大气参数传感系统和需求参数计算方法
CN114353739A (zh) * 2021-12-10 2022-04-15 武汉航空仪表有限责任公司 一种小型化多功能探头传感器

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101316478B (zh) * 2007-05-29 2012-02-08 清华大学 微型飞行器
CN101750514B (zh) * 2010-01-27 2011-12-28 中国科学院光电研究院 高空驻空飞行器风速风向实时测量方法与装置
CN102323440A (zh) * 2011-06-03 2012-01-18 清华大学 直升机二轴式空速检测系统
CN102662078A (zh) * 2012-05-07 2012-09-12 北京航空航天大学 一种直升机多轴式全向空速传感器
CN102662078B (zh) * 2012-05-07 2014-06-25 北京航空航天大学 一种直升机多轴式全向空速传感器
CN108139425A (zh) * 2015-10-28 2018-06-08 欧姆龙株式会社 空速测量系统
US11008118B2 (en) 2015-10-28 2021-05-18 Omron Corporation Airspeed measurement system
CN106372307A (zh) * 2016-08-30 2017-02-01 中国航空工业集团公司西安飞行自动控制研究所 一种基于气动模型的民机气流角估计方法
CN106372307B (zh) * 2016-08-30 2020-04-07 中国航空工业集团公司西安飞行自动控制研究所 一种基于气动模型的民机气流角估计方法
CN109781107A (zh) * 2017-11-15 2019-05-21 北京自动化控制设备研究所 一种低精度惯导滚动角确定方法
CN113654707A (zh) * 2021-08-12 2021-11-16 北京航空航天大学 基于压力和流速的大气参数传感系统和需求参数计算方法
CN114353739A (zh) * 2021-12-10 2022-04-15 武汉航空仪表有限责任公司 一种小型化多功能探头传感器

Also Published As

Publication number Publication date
CN100491949C (zh) 2009-05-27

Similar Documents

Publication Publication Date Title
CN1828253A (zh) 一种基于多流速传感器的飞行器飞行参数测量方法
US8408871B2 (en) Method and apparatus for measuring air flow condition at a wind turbine blade
CN201662459U (zh) 流场压力和速度的组合测量工具
CN203053472U (zh) 测量旋流场的压力和速度的工具
Gardner et al. Review of measurement techniques for unsteady helicopter rotor flows
Lueck Calculating the rate of dissipation of turbulent kinetic energy
CN112197932B (zh) 一种飞机防除冰试验模型表面温度修正方法及测量方法
JP2006343179A (ja) 境界層乱流遷移現象を利用したレイノルズ数等の計測
CN105675901A (zh) 用于估计飞行器空速的方法和设备
CN101050974A (zh) 上下游温度分布二次差动流量测量方法、传感器、及应用
Brezina et al. Measurement of static and dynamic performance characteristics of electric propulsion systems
JP2005522700A (ja) 風角度計測用センサ
CN106052668B (zh) 一种大量程硅微陀螺仪非线性数字补偿方法
CN105301275A (zh) 估算飞行器的马赫数的方法和装置
CN106679932B (zh) 基于微型柔性热膜传感器阵列的攻角测量方法
US7185538B2 (en) Methods and systems for sensing air vehicle airspeed
CN113340473B (zh) 一种转子表面热膜测试修正方法
Liu et al. A directional cylindrical anemometer with four sets of differential pressure sensors
Dellwik et al. Flux–profile relationships over a fetch limited beech forest
US10073921B1 (en) Thermal air data (TAD) system
Hannon et al. Trapezoidal wing experimental repeatability and velocity profiles in the 14-by 22-foot subsonic tunnel
EP3685170B1 (en) An airflow measurement device
CN109238738B (zh) 基于可测试性格栅的特种车辆冷却风量测试装置及方法
WO2014036701A1 (zh) 测量旋流场的压力和速度的工具
Halstead Boundary layer development in multi-stage low-pressure turbines

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090527

Termination date: 20170321

CF01 Termination of patent right due to non-payment of annual fee