[go: up one dir, main page]

CN110488853B - 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法 - Google Patents

一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法 Download PDF

Info

Publication number
CN110488853B
CN110488853B CN201910807545.6A CN201910807545A CN110488853B CN 110488853 B CN110488853 B CN 110488853B CN 201910807545 A CN201910807545 A CN 201910807545A CN 110488853 B CN110488853 B CN 110488853B
Authority
CN
China
Prior art keywords
stable
angular velocity
inner frame
inertial navigation
motor
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.)
Active
Application number
CN201910807545.6A
Other languages
English (en)
Other versions
CN110488853A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201910807545.6A priority Critical patent/CN110488853B/zh
Publication of CN110488853A publication Critical patent/CN110488853A/zh
Application granted granted Critical
Publication of CN110488853B publication Critical patent/CN110488853B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Manufacturing & Machinery (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法,通过该方法可将混合式惯导系统的台体稳定在地理坐标系,同时可对惯性器件进行旋转调制。该方法的优点,在于计算内框电机稳定控制指令时,针对转轴涡动与旋转调制角速度耦合对系统姿态误差的影响,将旋转角速度施加于台体坐标系z轴方向,提高了台体的稳定精度,降低了系统姿态误差。

Description

一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计 算方法
技术领域
本发明涉及一种提高台体稳定精度、降低系统姿态误差的稳定控制指令计算方法,适用于对系统姿态精度有较高要求、需要稳定在地理坐标系的惯性导航系统,属于伺服控制领域。
技术背景
惯导系统被广泛应用于军用、民用领域的多种运载器。传统的捷联式惯导系统和平台式惯导系统,导航精度的提升主要依靠提高惯性器件的性能。混合式惯导系统不仅可以隔离载体角运动,也可以对惯性器件进行旋转调制,从而在使用同等精度的惯性器件的条件下,显著提高导航精度。
为了达到隔离载体运动的目的,一般以陀螺输出为反馈量,计算稳定指令。但是对于双轴或三轴混合式惯导系统,旋转调制会造成陀螺敏感轴与电机轴不重合。因此,需要根据载体的不同运动、不同旋转策略,将陀螺输出折算到电机轴方向,得到沿电机轴方向的稳定指令。
同时,应该注意旋转调制引入的误差对混合式惯导系统姿态输出精度的影响。例如,对于绕内框轴长时间旋转调制的混合式惯导系统,由于内框轴的涡动,平台坐标系的z轴指向不固定。此时如果将旋转调制角速度施加到地理坐标系的天向,会导致平台坐标系的水平陀螺敏感到旋转调制角速度的投影,这会导致稳定指令计算错误,系统俯仰角和横滚角出现震荡。
因此,需要为混合式惯导系统设计一套稳定指令计算方法,使得旋转机构快速、平稳、准确地跟踪导航坐标系,尽量减小旋转调制角速度与转轴涡动的耦合误差。
发明内容
本发明解决的技术问题是:针对混合式惯导系统台体需稳定于地理坐标系以及隔离载体角运动的需求,减少涡动误差对系统姿态的影响,提出一种稳定控制指令的计算方法,在隔离载体角运动的同时提高系统姿态的精度。
本发明的技术解决方案是:一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法,实现步骤如下:
(1)进行姿态解算,计算平台坐标系(p系)相对于导航坐标系(n系)的旋转角速度
Figure GDA0002959092020000021
(2)在平台坐标系内,
Figure GDA0002959092020000022
扣除旋转调制角速度分量,得到稳定指令角速度
Figure GDA0002959092020000023
(3)对稳定指令角速度积分,得到稳定指令角度
Figure GDA0002959092020000024
(4)将稳定指令角度折算到电机轴上,得到电机稳定指令。
所述步骤(2)稳定指令角速度的计算公式如下:
Figure GDA0002959092020000025
式中,ΔRrx、ΔRry、ΔRrz分别为当前周期的中框、外框、内框的电机旋转指令。dt为控制周期。
所述步骤(3)稳定指令角度的计算公式如下:
Figure GDA0002959092020000026
本发明的原理是:
从陀螺输出中可以提取台体相对地理坐标系的角运动信息,以此为反馈,可以计算出电机轴的稳定指令,电机带动框架旋转,实现台体快速、准确地跟踪地理坐标系。
本发明与现有技术相比的优点在于:
(1)本发明适用于不同载体角运动、不同旋转策略下的稳定指令计算。由于载体运动、框架旋转会影响电机轴的正交关系以及陀螺敏感轴与电机轴的空间角度关系,需要将陀螺输出折算到电机轴方向,才能得到正确的稳定指令计算方法。本发明利用陀螺输出、系统姿态角、IMU姿态角等数据,给出了稳定指令计算方法。
(2)本发明减小了绕内框轴旋转调制时的系统姿态误差。绕内框轴旋转调制时,如果旋转调制角速度施加在地理坐标系天向,内框转轴涡动与旋转调制角速度耦合,会造成系统姿态误差。针对这一问题,本发明更改了稳定指令计算方法,将旋转调制角速度施加在平台坐标系z轴方向,可以减小系统姿态误差。
附图说明
图1为本发明所指的稳定指令计算方法流程图。图2为摘要附图。
具体实施方式
下面以一种三轴混合式惯性导航系统为例来阐述本发明的具体实施过程。
该惯导系统可以隔离载体角运动,三轴旋转机构的内框轴、中框轴、外框轴分别指向系统的方位、俯仰、横滚,三个框架各配一个增量式光电码盘。IMU安装于该旋转机构上。
根据图1所示的流程,电机稳定指令计算方法如下:
1.姿态解算,得到姿态矩阵
Figure GDA0002959092020000031
2.计算
Figure GDA0002959092020000032
更新地球自转以及地理系相对于地球系的角速率:
Figure GDA0002959092020000033
Figure GDA0002959092020000034
Figure GDA0002959092020000035
式中,ωie为地球自转角速度;L为当地地理纬度;vE、vN为上一周期计算出的东向速度、北向速度;RM、RN为当前周期计算出的地球曲率半径。
Figure GDA0002959092020000036
经过上一时刻的
Figure GDA0002959092020000037
转化到地理系得到
Figure GDA0002959092020000038
Figure GDA0002959092020000039
设陀螺输出为
Figure GDA00029590920200000310
扣除
Figure GDA00029590920200000311
得到
Figure GDA00029590920200000312
Figure GDA00029590920200000313
Figure GDA00029590920200000314
补偿圆锥误差后,可以用来计算稳定指令。
3.计算电机稳定指令。
光栅角为零时,x、y、z陀螺敏感轴基本指向中框轴、外框轴、内框轴。但是当载体存在俯仰运动时,三个电机轴之间不正交;当IMU绕某个电机轴进行旋转时,陀螺敏感轴与电机轴不重合;当载体俯仰角达到90°或者惯导系统绕中框轴旋转调制时,将无法隔离载体的航向运动。因此需要根据系统姿态以及不同旋转策略计算电机稳定指令。
(1)电机绕内框轴旋转时的稳定指令的计算公式如下:
Figure GDA00029590920200000315
Figure GDA0002959092020000041
式中,Rsx、Rsy、Rsz分别为中框、外框、内框的电机稳定指令;ψb为系统航向角,R(-ψb)是由-ψb构成的旋转矩阵;ΔRrx、ΔRry、ΔRrz分别为当前周期的中框、外框、内框的电机旋转指令;
Figure GDA0002959092020000042
为当前时刻的姿态矩阵。θb为系统俯仰角。
式(8)、式(9)适用于内框轴指天情况。
对于绕中框轴或者外框轴旋转180°,使内框轴指地旋转的情况,稳定指令计算公式如下:
Figure GDA0002959092020000043
式中,
Figure GDA0002959092020000044
由式(8)计算得到。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
最后所应说明的是:以上实施实例仅用以说明而非限制本发明的技术方案,所有的不脱离本发明的精神和范围的修改或局部替换,均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种混合式惯导系统稳定控制指令的计算方法,其特征在于实现步骤如下,
(1)进行姿态解算,计算平台坐标系相对于导航坐标系的旋转角速度
Figure FDA0002959092010000012
Figure FDA0002959092010000011
所述步骤1旋转角速度
Figure FDA0002959092010000012
的计算公式如下,
Figure FDA0002959092010000013
Figure FDA0002959092010000014
Figure FDA0002959092010000015
式中,ωie为地球自转角速度;L为当地地理纬度;vE、vN为上一周期计算出的东向速度、北向速度;RM、RN为当前周期计算出的地球曲率半径;
Figure FDA0002959092010000016
为上一时刻的姿态矩阵;
Figure FDA0002959092010000017
为当前时刻的陀螺输出,
(2)在平台坐标系内,
Figure FDA0002959092010000018
扣除旋转调制角速度分量,得到稳定指令角速度
Figure DEST_PATH_IMAGE002
Figure FDA0002959092010000019
所述步骤2稳定指令角速度的计算公式如下,
Figure FDA00029590920100000110
式中,ΔRrx、ΔRry、ΔRrz分别为当前周期的中框、外框、内框的电机旋转指令,dt为控制周期,
(3)对稳定指令角速度积分,得到稳定指令角度
Figure DEST_PATH_IMAGE004
Figure FDA00029590920100000111
Figure FDA00029590920100000112
(4)将稳定指令角度折算到电机轴上,得到电机稳定指令,
当内框轴指天,电机绕内框轴旋转时的稳定指令的计算公式如下,
Figure FDA0002959092010000021
Figure FDA0002959092010000022
式中,Rsx、Rsy、Rsz分别为中框、外框、内框的电机稳定指令;ψb为系统航向角,R(-ψb)是由-ψb构成的旋转矩阵;
Figure FDA0002959092010000023
为当前时刻的姿态矩阵, θb为系统俯仰角,
对于绕中框轴或者外框轴旋转180°,使内框轴指地旋转的情况,稳定指令计算公式如下,
Figure FDA0002959092010000024
CN201910807545.6A 2019-08-29 2019-08-29 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法 Active CN110488853B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910807545.6A CN110488853B (zh) 2019-08-29 2019-08-29 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910807545.6A CN110488853B (zh) 2019-08-29 2019-08-29 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法

Publications (2)

Publication Number Publication Date
CN110488853A CN110488853A (zh) 2019-11-22
CN110488853B true CN110488853B (zh) 2021-06-08

Family

ID=68553682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910807545.6A Active CN110488853B (zh) 2019-08-29 2019-08-29 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法

Country Status (1)

Country Link
CN (1) CN110488853B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114877890A (zh) * 2022-06-13 2022-08-09 武汉华中天易星惯科技有限公司 一种三轴混合式惯导系统的复合控制方法及控制装置
CN115077521B (zh) * 2022-08-22 2022-11-01 中国船舶重工集团公司第七0七研究所 一种基于虚拟框架载体坐标系惯性导航系统姿态解耦方法
CN115127552B (zh) * 2022-08-31 2022-11-18 中国船舶重工集团公司第七0七研究所 旋转调制方法、装置、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0619134A1 (de) * 1993-04-08 1994-10-12 ABB Management AG Mischkammer
WO2001046647A2 (en) * 1999-12-22 2001-06-28 Honeywell International Inc. Method, apparatus and computer program product for estimating airplane attitude with reduced sensor set
KR20150010318A (ko) * 2013-07-19 2015-01-28 현대모비스 주식회사 요 레이트 센서의 오프셋 보정 장치와 방법 및 상기 장치를 구비하는 차량 속도 제어 시스템

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05141308A (ja) * 1991-11-22 1993-06-08 Honda Motor Co Ltd 内燃機関のウオータジヤケツト構造
JPH09273938A (ja) * 1996-04-04 1997-10-21 Takao Yamaguchi ハイブリッド慣性航法装置
BRPI0511385A (pt) * 2004-06-15 2007-12-04 Aly El-Shafei métodos para controlar a instabilidade em rolamentos de pelìcula de fluido
FR2878954B1 (fr) * 2004-12-07 2007-03-30 Sagem Systeme de navigation inertielle hybride base sur un modele cinematique
CN101858353B (zh) * 2010-06-03 2012-09-19 浙江大学 一种用于离心泵的可控涡动装置
CN202452059U (zh) * 2012-01-12 2012-09-26 西安市瑞特测控技术有限责任公司 陀螺稳定云台
CN104064869B (zh) * 2014-06-13 2016-10-05 北京航天万达高科技有限公司 基于mems惯导的双四元数动中通天线控制方法及系统
CN106052682B (zh) * 2016-05-13 2018-12-11 北京航空航天大学 一种混合式惯性导航系统及导航方法
CN106289322B (zh) * 2016-08-23 2019-01-04 北京航空航天大学 一种基于陀螺数据的带旋转机构的惯导系统涡动补偿方法
CN106767806B (zh) * 2017-04-01 2019-10-22 北京航空航天大学 一种用于混合式惯性导航系统的物理平台
CN108873928A (zh) * 2017-05-10 2018-11-23 中国航空工业集团公司西安飞行自动控制研究所 一种基于传感器状态评估结果的无人直升机高度保持方法
CN108006410B (zh) * 2017-11-15 2019-11-01 湖北三江航天红峰控制有限公司 一种用于捷联惯导系统位置标定的双轴转位机构
CN207570556U (zh) * 2017-12-31 2018-07-03 宁波迪泰电子科技有限公司 捷联式惯性导航与校准系统
CN108869540B (zh) * 2018-01-12 2024-06-25 刘慕华 一种推力轴承、转子系统及推力轴承的控制方法
CN108318052B (zh) * 2018-01-24 2020-08-14 北京航天控制仪器研究所 一种基于双轴连续旋转的混合式平台惯导系统标定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0619134A1 (de) * 1993-04-08 1994-10-12 ABB Management AG Mischkammer
WO2001046647A2 (en) * 1999-12-22 2001-06-28 Honeywell International Inc. Method, apparatus and computer program product for estimating airplane attitude with reduced sensor set
KR20150010318A (ko) * 2013-07-19 2015-01-28 현대모비스 주식회사 요 레이트 센서의 오프셋 보정 장치와 방법 및 상기 장치를 구비하는 차량 속도 제어 시스템

Also Published As

Publication number Publication date
CN110488853A (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
CN109211269B (zh) 一种双轴旋转惯导系统姿态角误差标定方法
CN108318052B (zh) 一种基于双轴连续旋转的混合式平台惯导系统标定方法
CN106052682B (zh) 一种混合式惯性导航系统及导航方法
CN102749079B (zh) 一种光纤捷联惯导双轴旋转调制方法及双轴旋转机构
CN110488853B (zh) 一种降低转轴涡动影响的混合式惯导系统稳定控制指令的计算方法
CN104596546B (zh) 一种单轴旋转惯导系统的姿态输出补偿方法
CN103743413B (zh) 倾斜状态下调制寻北仪安装误差在线估计与寻北误差补偿方法
CN103630146B (zh) 一种离散解析与Kalman滤波结合的激光陀螺IMU标定方法
CN102288133B (zh) 一种陀螺间接稳定系统安装偏角标定方法
CN103292809A (zh) 一种单轴旋转式惯导系统及其专用误差自补偿方法
CN115574817B (zh) 一种基于三轴旋转式惯导系统的导航方法及导航系统
CN112179340B (zh) 一种冗余配置惯性测量单元双轴旋转调制方法
CN102798399A (zh) 基于双轴转动方案的捷联惯导系统误差抑制方法
CN108731674A (zh) 一种基于单轴旋转调制的惯性天文组合导航系统及计算方法
CN114061572A (zh) 一种用于旋转式惯导系统的双轴旋转调制方法
CN116481564B (zh) 基于Psi角误差修正模型的极地双惯导协同标定方法
Dong MEMS inertial navigation systems for aircraft
CN110501028A (zh) 一种用于双轴旋转mems-sins的十六位置旋转调制方法
CN111551164B (zh) 速率偏频激光陀螺寻北仪航向效应误差的补偿方法
CN110986934B (zh) 一体化双轴旋转惯导天文组合导航系统的导航方法及系统
CN116242397A (zh) 一种速度误差修正模型下的双惯导协同标定方法
CN115790668B (zh) 一种高精度惯导标度因数误差计算方法
CN114777777B (zh) 一种捷联化平台式惯导系统及导航方法
CN116519011B (zh) 基于Psi角误差修正模型的长航时双惯导协同标定方法
CN113108786B (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