CN112810834A - 一种同时考核惯性导航和模拟飞行的地面试验方法 - Google Patents

一种同时考核惯性导航和模拟飞行的地面试验方法 Download PDF

Info

Publication number
CN112810834A
CN112810834A CN202011542965.5A CN202011542965A CN112810834A CN 112810834 A CN112810834 A CN 112810834A CN 202011542965 A CN202011542965 A CN 202011542965A CN 112810834 A CN112810834 A CN 112810834A
Authority
CN
China
Prior art keywords
navigation
aircraft
flight
initial
inertial
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
CN202011542965.5A
Other languages
English (en)
Other versions
CN112810834B (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.)
Beijing Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Beijing Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN202011542965.5A priority Critical patent/CN112810834B/zh
Publication of CN112810834A publication Critical patent/CN112810834A/zh
Application granted granted Critical
Publication of CN112810834B publication Critical patent/CN112810834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种同时考核惯性导航和模拟飞行的地面试验方法,步骤包括:确定模拟飞行初始参数;获取初始时刻测试地点导航初始参数;通过推力曲线模拟轨道飞行;根据实测主从惯性测量组合实时输出的数据和测试点的初始速度、位置和姿态信息为初值以及主从惯性导航设备实时测量的陀螺数据和加表数据,主从惯组分别使用一套导航进行惯性导航解算,第一套导航获得主惯性测量组合导航结果,第二套导航获得从惯性测量组合导航结果,用于考核主从惯组的导航精度;由此,获得三套导航结果,第三套导航完成的飞行器推力曲线飞行模拟用于飞行控制。本发明采用推力曲线模拟,既可以模拟真实轨道,又反应了真实的惯组特性,使地面试验更加真实可靠。

Description

一种同时考核惯性导航和模拟飞行的地面试验方法
技术领域
本发明涉及一种同时考核惯性导航和模拟飞行的地面试验方法,属于运载火箭控制技术领域。
背景技术
飞行器研制过程中需要在飞行前对飞行过程中的各项控制过程和惯性测量组合产品进行地面模拟飞行的验证,因此,需要采用模拟飞行试验对飞行器制导控制流程和惯性测量组合硬件系统的全面考核,现有模拟飞行方法主要采用推力曲线或1g模拟飞行方法,而推力曲线模拟飞行主要采用软件计算模拟控制流程,没有考虑真实的惯组信息不能考核惯组精度,1g模拟飞行方法引入了真实的惯组信息,但不能完全模拟制导控制流程,两种方法各有优缺点,无法同时实现高效并采用惯组实测数据考核的目的。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提出一种同时考核惯性导航和模拟飞行的地面试验方法,建立三套导航,主从惯组分别使用一套导航,第三套导航方程参与模拟飞行控制,实现一次同时考核惯组和完成制导控制。
本发明解决技术的方案是:
一种同时考核惯性导航和模拟飞行的地面试验方法,具体步骤如下:
S1、确定模拟飞行初始参数
获取轨道起点的初始速度V0、位置S0和姿态信息θ0;并且获取标准轨道所有点的程序角
Figure BDA0002854448550000011
ψcx、γcx和发动机推力Pmf及飞行器质量Mmf0
S2、获取初始时刻测试地点导航初始参数
根据测试地点的纬度B0、射向A0,获得惯性导航初值:初始速度VG、位置SG和姿态信息θG
S3、通过推力曲线模拟轨道飞行
采用推力曲线模拟得到飞行器角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下的视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航完成飞行器推力曲线飞行模拟;
S4、根据实测主从惯性测量组合实时输出的数据和测试点的初始速度、位置和姿态信息为初值以及主从惯性导航设备实时测量的陀螺数据和加表数据,主从惯组分别使用一套导航进行惯性导航解算,第一套导航获得主惯性测量组合导航结果,第二套导航获得从惯性测量组合导航结果,用于考核主从惯组的导航精度;
S5、由此,获得三套导航结果,第三套导航完成的飞行器推力曲线飞行模拟用于飞行控制,第一套导航获得的主惯性测量组合导航结果和第二套导航获得从惯性测量组合导航结果用于考核主从惯组的导航精度。
进一步的,S3中,飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf通过如下方法得到:
程序角角增量=本周期程序角-上周期程序角
通过如下公式得到飞行器的程序角角增量:
Figure BDA0002854448550000021
Δψcx_mf=ψcxcx,-1
Δγcx_mf=γcxcx,-1
其中,
Figure BDA0002854448550000022
ψcx、γcx分别为导航坐标系下本周期俯仰、偏航、滚动方向的程序角,
Figure BDA0002854448550000023
ψcx,-1、γcx,-1为导航坐标系下上周期俯仰、偏航、滚动方向的程序角,
Figure BDA0002854448550000024
Δψcx_mf、Δγcx_mf为导航坐标系下俯仰、偏航、滚动方向的程序角角增量。
进一步的,由飞行器导航坐标系下的程序角角增量
Figure BDA0002854448550000025
Δψcx_mf、Δγcx_mf,通过如下公式得到飞行器体坐标系下的三方向角增量δθx_mf、δθy_mf、δθz_mf
Figure BDA0002854448550000031
进一步的,S3中,飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf通过如下方法得到:
飞行器质量Mmf通过如下公式获得:
Mmf=Mmf-1-dmCD_mf*T,其中Mmf-1为上一制导周期的飞行器质量,dmCD_mf为飞行器质量消耗率,T为制导周期;质量初值为Mmf0
进一步的,根据飞行器推力Pmf与飞行器质量Mmf,通过如下公式得到飞行器视速度增量δWx1_mf、δWy1_mf、δWz1_mf
Figure BDA0002854448550000032
进一步的,S3中,以Vt、St、θt作为推力曲线段模拟使用的初始惯性导航速度、位置和姿态信息,以及飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航进行惯性导航解算,直至飞行结束,获得推力曲线段轨道数据,完成接入推力曲线后飞行模拟。
本发明与现有技术相比的有益效果是:
(1)本发明采用推力曲线模拟,既可以模拟真实轨道,又反应了真实的惯组特性,使地面试验更加真实可靠;
(2)本发明两套惯组的惯性导航可以测试惯性测量组合硬件精度,相比于传统的模拟方法可实现软件复用且完成惯组考核;
(3)本发明的模拟方法简单、效率高。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
一种同时考核惯性导航和模拟飞行的地面试验方法,如图1所示,具体步骤如下:
S1、确定模拟飞行初始参数
获取轨道起点的初始速度V0、位置S0和姿态信息θ0;并且获取标准轨道所有点的程序角
Figure BDA0002854448550000045
ψcx、γcx和发动机推力Pmf及飞行器质量Mmf0
S2、获取初始时刻测试地点导航初始参数
根据测试地点的纬度B0、射向A0,获得惯性导航初值:初始速度VG、位置SG和姿态信息θG
S3、通过推力曲线模拟轨道飞行
采用推力曲线模拟得到飞行器角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下的视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航完成飞行器推力曲线飞行模拟,其中:
S3.1飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf通过如下方法得到:
程序角角增量=本周期程序角-上周期程序角
通过如下公式得到飞行器的程序角角增量:
Figure BDA0002854448550000041
Δψcx_mf=ψcxcx,-1
Δγcx_mf=γcxcx,-1
其中,
Figure BDA0002854448550000042
ψcx、γcx分别为导航坐标系下本周期俯仰、偏航、滚动方向的程序角,
Figure BDA0002854448550000043
ψcx,-1、γcx,-1为导航坐标系下上周期俯仰、偏航、滚动方向的程序角,
Figure BDA0002854448550000044
Δψcx_mf、Δγcx_mf为导航坐标系下俯仰、偏航、滚动方向的程序角角增量;
由飞行器导航坐标系下的程序角角增量
Figure BDA0002854448550000051
Δψcx_mf、Δγcx_mf通过如下公式得到飞行器体坐标系下的三方向角增量δθx_mf、δθy_mf、δθz_mf
Figure BDA0002854448550000052
S3.2飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf通过如下方法得到:
飞行器质量Mmf通过如下公式获得:
Mmf=Mmf-1-dmCD_mf*T,其中Mmf-1为上一制导周期的飞行器质量,dmCD_mf为飞行器质量消耗率,T为制导周期;质量初值为Mmf0
根据飞行器推力Pmf与飞行器质量Mmf通过如下公式得到飞行器视速度增量δWx1_mf、δWy1_mf、δWz1_mf
Figure BDA0002854448550000053
S3.3以Vt、St、θt作为推力曲线段模拟使用的初始惯性导航速度、位置和姿态信息,以及飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航进行惯性导航解算,直至飞行结束,获得推力曲线段轨道数据,完成接入推力曲线后飞行模拟;
S4、采用实测主从惯性测量组合输出的飞行模拟
采用实测主从惯性测量组合实时输出的数据和测试点的初始速度、位置和姿态信息为初值,采用主从惯性导航设备实时测量的陀螺数据和加表数据,主从惯组分别使用一套导航进行惯性导航解算,第一套导航获得主惯性测量组合导航结果,第二套导航获得从惯性测量组合导航结果,用于考核主从惯组的导航精度;
S5、由此,获得三套导航结果,第三套导航完成的飞行器推力曲线飞行模拟用于飞行控制,第一套导航获得的主惯性测量组合导航结果和第二套导航获得从惯性测量组合导航结果用于考核主从惯组的导航精度。
本发明采用推力曲线模拟,既可以模拟真实轨道,又反应了真实的惯组特性,使地面试验更加真实可靠;
本发明两套惯组的惯性导航可以测试惯性测量组合硬件精度,相比于传统的模拟方法可实现软件复用且完成惯组考核;本发明的模拟方法简单、效率高。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。

Claims (6)

1.一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,具体步骤如下:
S1、确定模拟飞行初始参数
获取轨道起点的初始速度V0、位置S0和姿态信息θ0;并且获取标准轨道所有点的程序角
Figure FDA0002854448540000011
ψcx、γcx和发动机推力Pmf及飞行器质量Mmf0
S2、获取初始时刻测试地点导航初始参数
根据测试地点的纬度B0、射向A0,获得惯性导航初值:初始速度VG、位置SG和姿态信息θG
S3、通过推力曲线模拟轨道飞行
采用推力曲线模拟得到飞行器角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下的视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航完成飞行器推力曲线飞行模拟;
S4、根据实测主从惯性测量组合实时输出的数据和测试点的初始速度、位置和姿态信息为初值以及主从惯性导航设备实时测量的陀螺数据和加表数据,主从惯组分别使用一套导航进行惯性导航解算,第一套导航获得主惯性测量组合导航结果,第二套导航获得从惯性测量组合导航结果,用于考核主从惯组的导航精度;
S5、由此,获得三套导航结果,第三套导航完成的飞行器推力曲线飞行模拟用于飞行控制,第一套导航获得的主惯性测量组合导航结果和第二套导航获得从惯性测量组合导航结果用于考核主从惯组的导航精度。
2.根据权利要求1所述的一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,S3中,飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf通过如下方法得到:
程序角角增量=本周期程序角-上周期程序角
通过如下公式得到飞行器的程序角角增量:
Figure FDA0002854448540000021
Δψcx_mf=ψcxcx,-1
Δγcx_mf=γcxcx,-1
其中,
Figure FDA0002854448540000022
ψcx、γcx分别为导航坐标系下本周期俯仰、偏航、滚动方向的程序角,
Figure FDA0002854448540000023
ψcx,-1、γcx,-1为导航坐标系下上周期俯仰、偏航、滚动方向的程序角,
Figure FDA0002854448540000024
Δψcx_mf、Δγcx_mf为导航坐标系下俯仰、偏航、滚动方向的程序角角增量。
3.根据权利要求2所述的一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,由飞行器导航坐标系下的程序角角增量
Figure FDA0002854448540000025
Δψcx_mf、Δγcx_mf,通过如下公式得到飞行器体坐标系下的三方向角增量δθx_mf、δθy_mf、δθz_mf
Figure FDA0002854448540000026
4.根据权利要求1所述的一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,S3中,飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf通过如下方法得到:
飞行器质量Mmf通过如下公式获得:
Mmf=Mmf-1-dmCD_mf*T,其中Mmf-1为上一制导周期的飞行器质量,dmCD_mf为飞行器质量消耗率,T为制导周期;质量初值为Mmf0
5.根据权利要求4所述的一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,根据飞行器推力Pmf与飞行器质量Mmf,通过如下公式得到飞行器视速度增量δWx1_mf、δWy1_mf、δWz1_mf
Figure FDA0002854448540000031
6.根据权利要求1所述的一种同时考核惯性导航和模拟飞行的地面试验方法,其特征在于,S3中,以Vt、St、θt作为推力曲线段模拟使用的初始惯性导航速度、位置和姿态信息,以及飞行器体坐标系下角增量δθx_mf、δθy_mf、δθz_mf和飞行器体坐标系下视速度增量δWx1_mf、δWy1_mf、δWz1_mf,使用第三套导航进行惯性导航解算,直至飞行结束,获得推力曲线段轨道数据,完成接入推力曲线后飞行模拟。
CN202011542965.5A 2020-12-23 2020-12-23 一种同时考核惯性导航和模拟飞行的地面试验方法 Active CN112810834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011542965.5A CN112810834B (zh) 2020-12-23 2020-12-23 一种同时考核惯性导航和模拟飞行的地面试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011542965.5A CN112810834B (zh) 2020-12-23 2020-12-23 一种同时考核惯性导航和模拟飞行的地面试验方法

Publications (2)

Publication Number Publication Date
CN112810834A true CN112810834A (zh) 2021-05-18
CN112810834B CN112810834B (zh) 2022-11-11

Family

ID=75853664

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011542965.5A Active CN112810834B (zh) 2020-12-23 2020-12-23 一种同时考核惯性导航和模拟飞行的地面试验方法

Country Status (1)

Country Link
CN (1) CN112810834B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100070178A1 (en) * 2008-09-15 2010-03-18 Cnh America Llc Method and system for vehicle orientation measurement
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法
CN106895855A (zh) * 2017-04-13 2017-06-27 北京航天自动控制研究所 一种惯性导航初始基准的估计与补偿方法
CN107037739A (zh) * 2016-12-02 2017-08-11 上海航天控制技术研究所 一种运载火箭半实物仿真试验惯组模拟方法
CN107063244A (zh) * 2017-04-14 2017-08-18 北京航天自动控制研究所 一种飞行器飞行过程模拟方法
CN112034703A (zh) * 2020-11-03 2020-12-04 蓝箭航天空间科技股份有限公司 航天运载器的自适应迭代制导方法及制导装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100070178A1 (en) * 2008-09-15 2010-03-18 Cnh America Llc Method and system for vehicle orientation measurement
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法
CN107037739A (zh) * 2016-12-02 2017-08-11 上海航天控制技术研究所 一种运载火箭半实物仿真试验惯组模拟方法
CN106895855A (zh) * 2017-04-13 2017-06-27 北京航天自动控制研究所 一种惯性导航初始基准的估计与补偿方法
CN107063244A (zh) * 2017-04-14 2017-08-18 北京航天自动控制研究所 一种飞行器飞行过程模拟方法
CN112034703A (zh) * 2020-11-03 2020-12-04 蓝箭航天空间科技股份有限公司 航天运载器的自适应迭代制导方法及制导装置

Also Published As

Publication number Publication date
CN112810834B (zh) 2022-11-11

Similar Documents

Publication Publication Date Title
CN108413887B (zh) 光纤光栅辅助分布式pos的机翼形变测量方法、装置和平台
CN108562288A (zh) 一种激光捷联惯组系统级在线自标定系统及方法
CN107421534B (zh) 一种冗余式捷联惯导系统多故障隔离方法
CN106969783A (zh) 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术
CN109635494A (zh) 一种飞行试验与地面仿真气动力数据综合建模方法
CN107037739B (zh) 一种运载火箭半实物仿真试验惯组模拟方法
CN105136166B (zh) 一种指定惯导位置精度的捷联惯导系统误差模型仿真方法
CN103984237B (zh) 基于运动状态综合识别的轴对称飞行器三通道自适应控制系统设计方法
CN107063244B (zh) 一种飞行器飞行过程模拟方法
CN101963513A (zh) 消除水下运载体捷联惯导系统杆臂效应误差的对准方法
CN103940442A (zh) 一种采用加速收敛算法的定位方法及装置
CN109059914B (zh) 一种基于gps和最小二乘滤波的炮弹滚转角估计方法
CN111537002B (zh) 一种激光捷联惯组安装误差的标定方法和定向方法
CN109297665B (zh) 使用高精准度气动力模型的捕获轨迹在线试验方法
CN109211231B (zh) 一种基于牛顿迭代法的炮弹姿态估计方法
CN106527122A (zh) 平流层飞艇定高飞行非线性pid控制方法
CN104121927A (zh) 一种适用于低精度无方位基准单轴转位设备的惯性测量单元标定方法
CN108680901A (zh) 一种新型的声源方位定位方法
CN103674068B (zh) 一种基于激光跟踪仪的传递对准的验证方法
CN101694390B (zh) 基于光纤惯性测量系统的舰船升沉运动测量方法
CN109211232B (zh) 一种基于最小二乘滤波的炮弹姿态估计方法
CN109029499B (zh) 一种基于重力视运动模型的加速度计零偏迭代寻优估计方法
CN112810834B (zh) 一种同时考核惯性导航和模拟飞行的地面试验方法
CN114326440A (zh) 一种用于运载火箭的仿真试验系统及方法
CN103983274A (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
GR01 Patent grant
GR01 Patent grant