CN112762960A - 一种飞行器所处风场的在线计算方法 - Google Patents
一种飞行器所处风场的在线计算方法 Download PDFInfo
- Publication number
- CN112762960A CN112762960A CN202011598737.XA CN202011598737A CN112762960A CN 112762960 A CN112762960 A CN 112762960A CN 202011598737 A CN202011598737 A CN 202011598737A CN 112762960 A CN112762960 A CN 112762960A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- wind field
- axis component
- speed
- ground
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 32
- 239000013598 vector Substances 0.000 claims description 71
- 230000009466 transformation Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 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 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 17
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C23/00—Combined instruments indicating more than one navigational value, e.g. for aircraft; Combined measuring devices for measuring two or more variables of movement, e.g. distance, speed or acceleration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
本发明属于飞行控制与飞行导航技术领域,提供一种飞行器所处风场的在线计算方法,通过飞行器地速向量与空速向量的向量合成和投影分解构造飞行器所处风场的信息,提高对飞行器所处风场测量的准确性。本发明能够为飞行器当地风场提供高精度、快动态、低成本的实时计算,有助于提升飞行测试的分析精度和飞行性能的计算精度,在航空、航天行业中有较高的应用价值。
Description
技术领域
本发明属于飞行控制与飞行导航技术领域,涉及一种可对飞行器所处风场进行在线、高精度计算的方法。
背景技术
风场会直接改变飞行器受到的气动力和力矩,进而影响飞行器的飞行姿态和飞行轨迹。在飞行试验中,大气中的自然风始终存在,并对飞行器产生持续干扰,使飞行试验的测试结果产生偏差。对飞行器所处风场的精准测量是从测试结果中消除偏差的技术基础。
此外,风场的存在使飞行器地速与空速产生差异,进而影响航程等飞行性能。例如,当飞行器以同样的空速飞行时,若为顺风,则地速大于空速,使飞行器的航程加长;若为逆风,则地速小于空速,使飞行器的航程缩短。因而对飞行器所处风场的精确测量也是飞行导航中精确计算飞行性能的技术基础。
为了对飞行器所处风场进行测量,通常的做法是从地面释放探空气球,在气球上升过程中使用气球搭载的仪器对各个高度的风速、风向进行测量,作为这一时段内飞行器在该区域飞行的大气风场数据。这一方法能获得测量点当地的瞬时风场数据,但当飞行器距离测量点较远时,使用测量点处的风场数据作为飞行器当地的风场数据将会带来大的偏差。
为提升对飞行器所处风场的测量精度和实时性,需要通过机载传感器信号对飞行器当地的风场数据进行实时、高精度的在线构造。
发明内容
本发明的目的是:提供一种飞行器所处风场的在线计算方法,提高对飞行器所处风场测量的准确性。
本发明的技术方案是:
一种飞行器所处风场的在线计算方法,通过飞行器地速向量与空速向量的向量合成和投影分解构造飞行器所处风场的信息。
包括飞行器地速的体轴分量解算、飞行器空速的体轴分量解算、风场向量的体轴分量解算、风场向量的地面坐标系投影、风场向量的航迹坐标系投影。
飞行器地速的体轴分量解算,使用飞行器北向地速、东向地速和地向地速,通过飞行器俯仰角、滚转角、偏航角的坐标变换,获得飞行器地速向量在机体坐标系的三轴分量,包括地速纵轴分量、地速侧轴分量和地速竖轴分量。
飞行器空速的体轴分量解算,使用飞行器真空速,通过迎角和侧滑角的坐标变换,获得飞行器空速向量在机体坐标系的三轴分量,包括空速纵轴分量、空速侧轴分量和空速竖轴分量。
风场向量的体轴分量解算,使用权利要求3中求得的地速体轴分量和权利要求4中求得的空速体轴分量,通过向量合成获得飞行器所处风场向量在机体坐标系中的三轴分量,包括风场纵轴分量、风场侧轴分量和风场竖轴分量。
风场向量的地面坐标系投影,使用权利要求5中求得的风场向量的体轴分量,通过飞行器俯仰角、滚转角、偏航角的坐标变换获得风场向量在地面坐标系中的三轴分量,包括北向风速、东向风速、地向风速;
风场向量的航迹坐标系投影,使用权利要求6中求得的风场向量的地面坐标系投影,通过航迹偏角、航迹倾角的坐标变换获得风场向量在航迹坐标系中的三轴分量,包括前向风速、侧向风速和垂向风速。
飞行器北向地速、东向地速、地向地速、飞行器真空速、迎角、侧滑角、飞行器俯仰角、滚转角、偏航角、航迹偏角、航迹倾角都是通过机载系统数据测量获得。
本发明所产生的有益效果:本发明一种飞行器所处风场的在线计算方法,使用飞行器地速、空速、迎角、侧滑角等机载传感器信号构造风场向量,能够获得飞行器当地的实时风场信息,避免传统探空气球方法使用测量点风场代替飞机当地风场带来的测量偏差。同时,本发明中,飞行器地速的体轴分量解算、飞行器空速的体轴分量解算、风场向量的体轴分量解算等算法均基于严格的数学推导建立,无估计、模拟等近似的计算方法,构造出的风场向量与飞行器当地的真实风场信息相比具有较高的精度。此外,本发明中,所有算法均为代数方程,无微分方程,无需考虑方程的收敛性问题,同时能避免微分方程动力学解算带来的信号延迟问题。最后,构造飞行器当地风场使用的飞行器地速、姿态角信号通常由惯性导航设备产生,空速、迎角、侧滑角通常由大气数据系统产生,惯性导航设备和大气数据系统均为现代飞行控制系统的常用设备,对风场向量的构造通过软件和算法即可实现,无需增加额外的硬件设备,具有良好的经济性。综合以上特点,本发明能够为飞行器当地风场提供高精度、快动态、低成本的实时计算,有助于提升飞行测试的分析精度和飞行性能的计算精度,在航空、航天行业中有较高的应用价值。
附图说明
图1是本发明的方法原理图。
具体实施方式
下面结合附图和实施方式对一种飞行器所处风场的在线计算方法做具体说明。
通过飞行器地速向量与空速向量的向量合成和投影分解构造飞行器所处风场的信息,如图1所示。具体包括:飞行器地速的体轴分量解算、飞行器空速的体轴分量解算、风场向量的体轴分量解算、风场向量的地面坐标系投影、风场向量的航迹坐标系投影。
(1)飞行器地速的体轴分量解算。使用飞行器北向地速、东向地速和地向地速构成飞行器相对于地面运动的速度(地速)向量,使用俯仰角、滚转角、偏航角表征飞行器相对于地面的姿态,通过俯仰角、滚转角、偏航角构造地面坐标系到飞行器机体坐标系的转换矩阵,通过转换矩阵与飞行器地速向量的矩阵乘法运算,获得飞行器地速在机体坐标系的三轴分量,包括地速纵轴分量、地速侧轴分量和地速竖轴分量。
(2)飞行器空速的体轴分量解算。使用飞行器真空速表征飞行器相对于大气运动的速度(空速)向量,使用迎角、侧滑角表征飞行器相对于大气来流的姿态,通过迎角、侧滑角构造空速坐标系到飞行器机体坐标系的转换矩阵,通过转化矩阵与飞行器空速向量的矩阵乘法运算,获得飞行器空速在机体坐标系的三轴分量,包括空速纵轴分量、空速侧轴分量和空速竖轴分量。
(3)风场向量的体轴分量解算。使用步骤(1)中求得的地速体轴分量和步骤(2)中求得的空速体轴分量,在机体坐标系内通过向量合成获得飞行器所处风场向量在机体坐标系中的三轴分量,包括风场纵轴分量、风场侧轴分量和风场竖轴分量。
(4)风场向量的地面坐标系投影。使用步骤(3)中求得的风场向量的体轴分量,通过俯仰角、滚转角、偏航角构造飞行器机体坐标系到地面坐标系的转换矩阵,通过转换矩阵与风场向量体轴分量的矩阵乘法运算,获得风场向量在地面坐标系的三轴分量,包括:北向风速、东向风速、地向风速。
(5)风场向量的航迹坐标系投影。使用步骤(4)中求得的风场向量的地面坐标系投影,通过航迹偏角、航迹倾角构造地面坐标系到航迹坐标系的转换矩阵,通过转换矩阵与风场向量地面坐标系分量的矩阵乘法运算,获得风场向量在航迹坐标系中的三轴分量,包括前向风速、侧向风速和垂向风速。
具体实现方式如下:
(1)飞行器地速的体轴分量解算
设飞行器相对于地面运动的北向速度为VN、东向速度为VE、地向速度为VD,则飞行器相对于地面运动的速度(地速)向量可表达为:[VN VE VD]T。
设飞行器相对于地面坐标系的俯仰角为θ、滚转角为φ、偏航角为ψ,则地面坐标系i到飞行器机体坐标系b的坐标变换矩阵为:
设飞行器地速在机体坐标系的纵轴分量为u,侧轴分量为v,竖轴分量为w,则地速在机体坐标系的三轴分量计算方法为:
(2)飞行器空速的体轴分量解算
设大气数据系统测得的飞行器真空速为VT,则飞行器相对于大气的速度(空速)向量在空速坐标系a中可表达为:[VT 0 0]T。
设飞行器相对于大气来流的迎角为α、侧滑角为β,则空速坐标系a到飞行器机体坐标系b的坐标变换矩阵为:
设飞行器空速在机体坐标系的纵轴分量为ua,侧轴分量为va,竖轴分量为wa,则空速在机体坐标系的三轴分量计算方法为:
(3)风场向量的体轴分量解算
设风场向量在机体坐标系的纵轴分量为uw,侧轴分量为vw,竖轴分量为ww,则风场向量在机体坐标系的三轴分量计算方法为:
(4)风场向量的地面坐标系投影
设飞行器相对于地面坐标系的俯仰角为θ、滚转角为φ、偏航角为ψ,则飞行器机体坐标系b到地面坐标系i的坐标变换矩阵为:
设风场向量在地面坐标系的三轴分量为:北向风速Vw,N,东向风速Vw,E,地向风速为Vw,D,则风场向量在地面坐标系的三轴分量计算方法为:
(5)风场向量的航迹坐标系投影
设飞行器地速向量相对于地面坐标系的航迹倾角为γ,航迹偏角为χ,则从地面坐标系i到飞行器航迹坐标系k的坐标变换矩阵为:
设风场向量在航迹坐标系的前向分量为Vw,X,侧向分量为Vw,Y,垂向分量为Vw,Z,则风场向量在航迹坐标系中的三轴分量计算方法为:
Claims (8)
1.一种飞行器所处风场的在线计算方法,其特征在于:通过飞行器地速向量与空速向量的向量合成和投影分解构造飞行器所处风场的信息。
2.如权利要求1所述的一种飞行器所处风场的在线计算方法,其特征在于:包括飞行器地速的体轴分量解算、飞行器空速的体轴分量解算、风场向量的体轴分量解算、风场向量的地面坐标系投影、风场向量的航迹坐标系投影。
3.如权利要求2所述的一种飞行器所处风场的在线计算方法,其特征在于:飞行器地速的体轴分量解算,使用飞行器北向地速、东向地速和地向地速,通过飞行器俯仰角、滚转角、偏航角的坐标变换,获得飞行器地速向量在机体坐标系的三轴分量,包括地速纵轴分量、地速侧轴分量和地速竖轴分量。
4.如权利要求2所述的一种飞行器所处风场的在线计算方法,其特征在于:飞行器空速的体轴分量解算,使用飞行器真空速,通过迎角和侧滑角的坐标变换,获得飞行器空速向量在机体坐标系的三轴分量,包括空速纵轴分量、空速侧轴分量和空速竖轴分量。
5.如权利要求2所述的一种飞行器所处风场的在线计算方法,其特征在于:风场向量的体轴分量解算,使用权利要求3中求得的地速体轴分量和权利要求4中求得的空速体轴分量,通过向量合成获得飞行器所处风场向量在机体坐标系中的三轴分量,包括风场纵轴分量、风场侧轴分量和风场竖轴分量。
6.如权利要求2所述的一种飞行器所处风场的在线计算方法,其特征在于:风场向量的地面坐标系投影,使用权利要求5中求得的风场向量的体轴分量,通过飞行器俯仰角、滚转角、偏航角的坐标变换获得风场向量在地面坐标系中的三轴分量,包括北向风速、东向风速、地向风速。
7.如权利要求2所述的一种飞行器所处风场的在线计算方法,其特征在于:风场向量的航迹坐标系投影,使用权利要求6中求得的风场向量的地面坐标系投影,通过航迹偏角、航迹倾角的坐标变换获得风场向量在航迹坐标系中的三轴分量,包括前向风速、侧向风速和垂向风速。
8.如权利要求3-7中的任意一项所述的一种飞行器所处风场的在线计算方法,其特征在于:
飞行器北向地速、东向地速、地向地速、飞行器真空速、迎角、侧滑角、飞行器俯仰角、滚转角、偏航角、航迹偏角、航迹倾角都是通过机载系统数据测量获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011598737.XA CN112762960A (zh) | 2020-12-29 | 2020-12-29 | 一种飞行器所处风场的在线计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011598737.XA CN112762960A (zh) | 2020-12-29 | 2020-12-29 | 一种飞行器所处风场的在线计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112762960A true CN112762960A (zh) | 2021-05-07 |
Family
ID=75697045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011598737.XA Pending CN112762960A (zh) | 2020-12-29 | 2020-12-29 | 一种飞行器所处风场的在线计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112762960A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281531A (zh) * | 2021-05-20 | 2021-08-20 | 北京科技大学 | 一种无人机当前风速风向测量方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2461801C1 (ru) * | 2011-09-20 | 2012-09-20 | Открытое акционерное общество "Раменское приборостроительное конструкторское бюро" (ОАО "РПКБ") | Способ определения скорости ветра на борту летательного аппарата и комплексная навигационная система для его реализации |
WO2013174559A1 (en) * | 2012-05-25 | 2013-11-28 | Prox Dynamics As | Method and device for estimating a wind field |
CN106324643A (zh) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | 一种无人机空速估计和空速管故障检测方法 |
CN108152529A (zh) * | 2017-11-02 | 2018-06-12 | 成都飞机工业(集团)有限责任公司 | 一种基于飞行参数计算风速及风向的方法 |
CN109374924A (zh) * | 2018-10-25 | 2019-02-22 | 西北工业大学 | 一种基于四旋翼无人机倾斜角的横侧向风场估计方法 |
CN109541963A (zh) * | 2018-11-12 | 2019-03-29 | 北京应用气象研究所 | 一种基于侧滑角信息的无人机测风建模技术 |
CN109781138A (zh) * | 2019-01-25 | 2019-05-21 | 北京电子工程总体研究所 | 一种基于浮空平台的风速风向测量方法及装置 |
-
2020
- 2020-12-29 CN CN202011598737.XA patent/CN112762960A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2461801C1 (ru) * | 2011-09-20 | 2012-09-20 | Открытое акционерное общество "Раменское приборостроительное конструкторское бюро" (ОАО "РПКБ") | Способ определения скорости ветра на борту летательного аппарата и комплексная навигационная система для его реализации |
WO2013174559A1 (en) * | 2012-05-25 | 2013-11-28 | Prox Dynamics As | Method and device for estimating a wind field |
CN106324643A (zh) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | 一种无人机空速估计和空速管故障检测方法 |
CN108152529A (zh) * | 2017-11-02 | 2018-06-12 | 成都飞机工业(集团)有限责任公司 | 一种基于飞行参数计算风速及风向的方法 |
CN109374924A (zh) * | 2018-10-25 | 2019-02-22 | 西北工业大学 | 一种基于四旋翼无人机倾斜角的横侧向风场估计方法 |
CN109541963A (zh) * | 2018-11-12 | 2019-03-29 | 北京应用气象研究所 | 一种基于侧滑角信息的无人机测风建模技术 |
CN109781138A (zh) * | 2019-01-25 | 2019-05-21 | 北京电子工程总体研究所 | 一种基于浮空平台的风速风向测量方法及装置 |
Non-Patent Citations (2)
Title |
---|
刘骏,范彦铭,郑伟: "一种通过ADS/INS 系统估算风速的方法", 《飞机设计》 * |
袁燎原: "基于信息融合的飞行风场估计技术", 《机电信息》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281531A (zh) * | 2021-05-20 | 2021-08-20 | 北京科技大学 | 一种无人机当前风速风向测量方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108152529A (zh) | 一种基于飞行参数计算风速及风向的方法 | |
CN107544067B (zh) | 一种基于高斯混合近似的高超声速再入飞行器跟踪方法 | |
CN110672092B (zh) | 一种降低固定翼无人机平台磁干扰的航迹生成方法 | |
CN102607639A (zh) | 基于bp神经网络的大攻角飞行状态下大气数据测量方法 | |
CN109635494A (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN111766397B (zh) | 一种基于惯性/卫星/大气组合的气象风测量方法 | |
CN112683446B (zh) | 一种飞机实时重心位置估计方法 | |
CN106248082B (zh) | 一种飞行器自主导航系统及导航方法 | |
CN106885573A (zh) | 面向四旋翼飞行器的运动捕捉系统实时测姿方法 | |
CN114636842B (zh) | 一种高超声速飞行器的大气数据估计方法及装置 | |
CN112556724A (zh) | 动态环境下的微型飞行器低成本导航系统初始粗对准方法 | |
CN109141802A (zh) | 用于捕获轨迹试验中的外挂物控制律的模拟方法 | |
RU2564379C1 (ru) | Бесплатформенная инерциальная курсовертикаль | |
CN109599674B (zh) | 一种基于解耦的相控阵天线稳定角跟踪方法 | |
CN109541963B (zh) | 一种基于侧滑角信息的无人机测风建模方法 | |
CN112762960A (zh) | 一种飞行器所处风场的在线计算方法 | |
CN111412887B (zh) | 一种基于卡尔曼滤波的攻角、侧滑角辨识方法 | |
RU2373562C2 (ru) | Способ и устройство контроля горизонтальной ориентации аппарата | |
CN109445283A (zh) | 一种用于欠驱动浮空器在平面上定点跟踪的控制方法 | |
CN114325666A (zh) | 高速运动下的相干外差测风雷达修正方法 | |
CN113188540B (zh) | 基于恒星数目和构型的惯性/天文自适应滤波方法 | |
CN110017809B (zh) | 利用地磁信息和光流传感器解算飞行器姿态的方法 | |
CN110017831B (zh) | 利用地磁信息和声呐传感器解算飞行器姿态的方法 | |
Cho et al. | Air data system calibration using GPS velocity information | |
Sun et al. | A gnss/imu-based 5-hole pitot tube calibration algorithm |
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: 20210507 |