CN100491949C - 一种基于多流速传感器的飞行器飞行参数测量方法 - Google Patents
一种基于多流速传感器的飞行器飞行参数测量方法 Download PDFInfo
- Publication number
- CN100491949C CN100491949C CNB2006100656815A CN200610065681A CN100491949C CN 100491949 C CN100491949 C CN 100491949C CN B2006100656815 A CNB2006100656815 A CN B2006100656815A CN 200610065681 A CN200610065681 A CN 200610065681A CN 100491949 C CN100491949 C CN 100491949C
- Authority
- CN
- China
- Prior art keywords
- aircraft
- angle
- attack
- speed
- group
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000014509 gene expression Effects 0.000 claims description 10
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 18
- 229910052697 platinum Inorganic materials 0.000 description 9
- 238000005259 measurement Methods 0.000 description 8
- 238000000691 measurement method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000000919 ceramic Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000010409 thin film Substances 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000000703 anti-shock Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002775 capsule Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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,飞行器参数按照如下公式计算:
飞行器空速:
飞行器迎角:α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的一组以上流速传感器,对每组传感器测得的数值取平均值,飞行器参数按照如下公式计算:
飞行器空速:
飞行器迎角:
飞行器侧滑角:
飞行器滚转角速率:
在上述技术方案中,进一步地,如图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左右关于轴线OO’对称,翼展300mm,蒙皮厚0.1mm。
四个流速传感器,为第一流速传感器1、第二流速传感器2、第三流速传感器3、第四流速传感器4分别按顺时针表贴在飞行器机翼前缘左右对称位置的上下表面。从正视飞行器机翼前缘的方向看,这4个传感器按顺时针方向依次为1、2、3、4。本实施例中选用热式流速传感器,其敏感元件为由陶瓷封装的薄膜铂热敏电阻,在市场上买到。铂热敏电阻本身性能稳定,从而使整个传感器网络有稳定的性能。铂薄膜通过激光喷溅在陶瓷表层,因而元件具有优良的防振和防冲击性能,比较适合用于飞行器。
铂热敏电阻通电发热后,其阻值会随着流经传感器的空速的改变而改变,因而可以通过测量热敏电阻的阻值来感受空气流速的变化。每个热式流速传感器均采用恒压式工作方式,其电路原理图如附图3所示。接点5处接恒压源(图中未示出),数字6表示定值电阻,数字7表示第一铂热敏电阻,作流速敏感元件,数字8为第二铂热敏电阻,放置在飞行器内部,对第一铂热敏电阻7起环境温度补偿作用,数字9表示调零电阻。第一铂热敏电阻7的工作电阻Rw和其所在位置的空气流速v之间有如下关系:
式中,Rc——定值电阻6的阻值;
U——电桥桥顶电压,为定值;
Rg——环境温度下铂热敏电阻7的阻值;
A,B——常数,由预先标定确定。
式中,常数A、B的标定在零侧滑角、零迎角及零滚转角速率下进行。测量不同空气流速下传感器的输出,用得到的数据,采用最小二乘法估计参数A、B。
由第一、第二、第三和第四流速传感器测得的相应位置的流速值s1、s2、s3和s4,可直接得到飞行器的空速,即:
飞行器空速:
由于飞行器空速对迎角、侧滑角和滚转角速率的测量值有影响,所以需要在飞行器实际飞行时可能的空速范围内对迎角、侧滑角和滚转角速率进行标定,标定时还要考虑迎角与侧滑角之间的耦合关系。本实施例中,迎角的标定范围为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,飞行器参数按照如下公式计算:
飞行器空速;
飞行器迎角:α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的一组以上流速传感器,对每组传感器测得的数值取平均值,飞行器参数按照如下公式计算:
飞行器空速
飞行器迎角
飞行器侧滑角
飞行器滚转角速率
2、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述步骤1)中每组的四个传感器分别布置在机翼前缘上、下、左、右四个对称位置。
3、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述步骤1)中在机翼前缘位置布置多组传感器,且流速传感器之间互不干扰。
4、根据权利要求1所述基于多流速传感器的飞行器飞行参数的测量方法,其特征是,所述空速由四个传感器测得的流速值直接得到;迎角和侧滑角之间有耦合关系,根据步骤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,记:
ωig=s1g-s2g+s3g-s4g
采用参数估计方法估计参数k5,由此,得到滚转角速率的表达式如下:
飞行器滚转角速率:ω=k5(s1-s2+s3-s4)
(3)在其它空速下,重复步骤(1)—(2)做同样标定,得到不同空速下相关的系数k,飞行器实际飞行时根据测得的空速选取相应系数k;或者,预先标定拟和各系数k关于空速u的函数关系式,用以求解飞行器飞行时的其它三个飞行参数。
6、根据权利要求5所述基于多流速传感器的飞行器飞行参数的测量方法,其特征在于,所述参数估计方法包括最小二乘法、极大似然估计法。
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 CN1828253A (zh) | 2006-09-06 |
CN100491949C true 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) |
Families Citing this family (9)
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 | 清华大学 | 直升机二轴式空速检测系统 |
CN102662078B (zh) * | 2012-05-07 | 2014-06-25 | 北京航空航天大学 | 一种直升机多轴式全向空速传感器 |
JP6673667B2 (ja) * | 2015-10-28 | 2020-03-25 | オムロン株式会社 | 対気速度計測システム |
CN106372307B (zh) * | 2016-08-30 | 2020-04-07 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种基于气动模型的民机气流角估计方法 |
CN109781107A (zh) * | 2017-11-15 | 2019-05-21 | 北京自动化控制设备研究所 | 一种低精度惯导滚动角确定方法 |
CN113654707B (zh) * | 2021-08-12 | 2022-08-16 | 北京航空航天大学 | 基于压力和流速的大气参数传感系统和需求参数计算方法 |
CN114353739A (zh) * | 2021-12-10 | 2022-04-15 | 武汉航空仪表有限责任公司 | 一种小型化多功能探头传感器 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5639964A (en) * | 1994-10-24 | 1997-06-17 | Djorup; Robert S. | Thermal anemometer airstream turbulent energy detector |
US6531967B2 (en) * | 2001-05-16 | 2003-03-11 | Robert S. Djorup | Aircraft disturbed air direct sensing system |
CN1664506A (zh) * | 2004-03-05 | 2005-09-07 | 清华大学 | 一种载体姿态测量方法及其系统 |
-
2006
- 2006-03-21 CN CNB2006100656815A patent/CN100491949C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5639964A (en) * | 1994-10-24 | 1997-06-17 | Djorup; Robert S. | Thermal anemometer airstream turbulent energy detector |
US6531967B2 (en) * | 2001-05-16 | 2003-03-11 | Robert S. Djorup | Aircraft disturbed air direct sensing system |
CN1664506A (zh) * | 2004-03-05 | 2005-09-07 | 清华大学 | 一种载体姿态测量方法及其系统 |
Non-Patent Citations (1)
Title |
---|
iMEMS 速率陀螺芯片在MAV 飞行控制系统中应用研究. 曾庆华,张为华.测控技术,第23卷第2期. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN1828253A (zh) | 2006-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100491949C (zh) | 一种基于多流速传感器的飞行器飞行参数测量方法 | |
US7908096B2 (en) | Integrated micromachined thermal mass flow sensor and methods of making the same | |
CN105548606B (zh) | 基于mems的柔性流速传感器的流速测量方法 | |
US6490510B1 (en) | Fixed multifunction probe for aircraft | |
US6609421B2 (en) | Sideslip correction for a multi-function three probe air data system | |
JPH0257853B2 (zh) | ||
CN104048808A (zh) | 一种动态熵探针 | |
CN108051610B (zh) | 一种智能双检测模式的风速风向传感器及测量方法 | |
Svedin et al. | A new silicon gas-flow sensor based on lift force | |
Zhu et al. | Sensitivity improvement of a 2D MEMS thermal wind sensor for low-power applications | |
JP2005522700A (ja) | 風角度計測用センサ | |
Xu et al. | High sensitivity and wide dynamic range thermoresistive micro calorimetric flow sensor with CMOS MEMS technology | |
CN108195510A (zh) | 一种热膜式剪应力传感器的热风洞校准方法 | |
Mansoor et al. | SOI CMOS multi-sensors MEMS chip for aerospace applications | |
US7185538B2 (en) | Methods and systems for sensing air vehicle airspeed | |
Löfdahl et al. | Characteristics of a hot-wire microsensor for time-dependent wall shear stress measurements | |
CN106679932A (zh) | 基于微型柔性热膜传感器阵列的攻角测量方法 | |
Wang et al. | Temperature effects of a ceramic MEMS thermal wind sensor based on a temperature-balanced mode | |
EP3685170B1 (en) | An airflow measurement device | |
Svedin et al. | A lift force sensor with integrated hot-chips for wide range flow measurements | |
Zhu et al. | A robust and low-power 2-D thermal wind sensor based on a glass-in-silicon reflow process | |
CN112556917B (zh) | 一种利用测压装置进行测压的方法 | |
CN108519189B (zh) | 一种恒温式热膜剪应力微传感器闭环反馈控制系统 | |
Beutel et al. | Design and evaluation process of a robust pressure sensor for measurements in boundary layers of liquid fluids | |
Wei et al. | Hybrid system for airspeed measurement using dual MEMS sensors |
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 |