CN110849294B - 一种基于光纤陀螺的转台非正交度测试方法 - Google Patents

一种基于光纤陀螺的转台非正交度测试方法 Download PDF

Info

Publication number
CN110849294B
CN110849294B CN201911145605.9A CN201911145605A CN110849294B CN 110849294 B CN110849294 B CN 110849294B CN 201911145605 A CN201911145605 A CN 201911145605A CN 110849294 B CN110849294 B CN 110849294B
Authority
CN
China
Prior art keywords
turntable
outer frame
frame
inner frame
axis
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
CN201911145605.9A
Other languages
English (en)
Other versions
CN110849294A (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.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
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 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN201911145605.9A priority Critical patent/CN110849294B/zh
Publication of CN110849294A publication Critical patent/CN110849294A/zh
Application granted granted Critical
Publication of CN110849294B publication Critical patent/CN110849294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Gyroscopes (AREA)

Abstract

本发明涉及一种基于光纤陀螺的转台非正交度测试方法,其特征在于:包括1)转台结构与陀螺组件建模;2)转台旋转次序与模型;3)转台非正交度计算。本发明设计科学合理,在不使用光电自准直仪的条件下,使用三只光纤陀螺相互正交安装构建陀螺组件,基于陀螺组件实现转台转轴之间非正交度的测试,该技术通用性强,可以简单方便实现测试,而且不需拆装转台。

Description

一种基于光纤陀螺的转台非正交度测试方法
技术领域
本发明属于惯性器件测试领域,涉及脱离光电自准直仪的测试设备,特别涉及一种基于光纤陀螺的转台非正交度测试方法。
背景技术
随着导航技术的发展,转台作为一种测试仪器在惯性领域的使用越来越广泛,各类精密仪器的标定和性能测试等试验的准确性都依赖于转台的精度。转台各转轴之间的非正交度是转台精度一个非常重要的指标。
通常对转台的各转轴之间非正度的测试是通过架设光电自准直仪,在转轴上安装反射镜,使用光学的方法进行测试标校。在实际使用中,需要对转台外壳进行拆装,架设光电自准直仪也需要在转台周围具有合适的工作空间,整个过程操作复杂,所需时间也较长。
因此寻找一种测试精度高,简单方便而且易于操作的转台非正交度测试方法,具有十分重要的意义。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于光纤陀螺的转台非正交度测试方法,在不使用光电自准直仪的条件下,使用三只光纤陀螺相互正交安装构建陀螺组件,基于陀螺组件实现转台转轴之间非正交度的测试,通用性强,可以简单方便实现测试,而且不需拆装转台。
本发明解决其技术问题是通过以下技术方案实现的:
一种基于光纤陀螺的转台非正交度测试方法,其特征在于:所述测试方法的步骤为:
1)转台结构与陀螺组件建模:
1、转台模型:建立转台模型并将转台结构分为外框方位、中框俯仰及内框横滚,并根据上述结构分别建立转台外框坐标系T1,中框坐标系T2,内框坐标系T3,陀螺三轴坐标系g系,其中g系是非正交坐标系;
从T1向T2转换,存在存在中框转轴与理想中框轴的夹角θy1和θz1,以及中框本身构建的水平面与外框旋转平面关于中框转轴的旋转角度,设为α0,那么从T1系向T2系转换如下:
Figure BDA0002282091040000021
从T2向T3转换,存在存在内框轴与理想内框轴的夹角θx2和θz2,以及内框本身构建的水平面与中框平面关于内框转轴的旋转角度,设为β0,那么从T2系向T3系转换如下:
Figure BDA0002282091040000022
2、陀螺组件模型:
单只陀螺的参数包括标度因数、零位偏置及敏感轴偏置角度,由三只陀螺构成的陀螺组件模型参数为12个,包括陀螺标度因数Kgx,Kgy,Kgz,陀螺安装偏角Ugxy,Ugxz,Ugyx,Ugyz,Ugzx,Ugzy,陀螺零偏Bgx0,Bgy0,Bgz0;
陀螺组件的数学模型如下:
Figure BDA0002282091040000023
其中:ωxyz分别为三个相互正交坐标轴上的角速度输入;
Ngx,Ngy,Ngz分别为三个陀螺脉冲数输出;
2)转台旋转次序与模型
将步骤1)中的陀螺组件安装到三轴转台上,预先定义的X,Y和Z轴分别和转台的中框、内框和外框平行,旋转顺序依次为:
1、转台外框、中框、内框归零,Z轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000031
2、中框转至90°,外框、内框归零,Y轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000032
3、内框转至270°,外框、中框归零,X轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000033
4、中框转至180°,外框、内框归零,Z轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000034
5、中框转至270°,外框、内框归零,Y轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000041
6、内框转至90°,外框、中框归零,X轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000042
3)转台非正交度计算
结合步骤1)、2),将求解出的三个陀螺标度因数带入可求出如下变量,见表1:
表1
各轴指天 各轴指地
U<sub>gxy</sub>+θ<sub>z2</sub>-θ<sub>y1</sub> U<sub>gxy</sub>+θ<sub>z2</sub>+θ<sub>y1</sub>
U<sub>gxz</sub>-β<sub>0</sub>-θ<sub>y1</sub> U<sub>gxz</sub>-β<sub>0</sub>+θ<sub>y1</sub>
U<sub>gyx</sub>+θ<sub>x2</sub>+α<sub>0</sub> U<sub>gyx</sub>-θ<sub>x2</sub>-α<sub>0</sub>
U<sub>gyz</sub>+θ<sub>x2</sub>+α<sub>0</sub> U<sub>gyz</sub>+θ<sub>x2</sub>+α<sub>0</sub>
U<sub>gzx</sub>+β<sub>0</sub>+θ<sub>y1</sub> U<sub>gzx</sub>+β<sub>0</sub>+θ<sub>y1</sub>
U<sub>gzy</sub>-θ<sub>x2</sub>-α<sub>0</sub> U<sub>gzy</sub>-θ<sub>x2</sub>-α<sub>0</sub>
由此,可以进一步确定以下变量的值:
Ugxyz2,θy1,Ugxz0,Ugyx,θx20,Ugyz,Ugzx0,Ugzy
同时,由于三个坐标轴的两两之间不正交度的关系,只有三个自由度,因此还具备如下约束条件:
Ugxy+Ugyx=A
Ugxz+Ugzx=B
Ugyz+Ugzy=C
其中:A、B、C对于同一套陀螺组件是固定值,可以事先标定精确得到,进一步可以确认变量θz2的值,因此可计算出转台转轴之间的不正交角为:θy1,θx20,θz2,其中θx2与α0耦合。
本发明的优点和有益效果为:
1、本发明的基于光纤陀螺的转台非正交度测试方法,在不使用光电自准直仪的条件下,使用三只光纤陀螺相互正交安装构建陀螺组件,基于陀螺组件实现转台转轴之间非正交度的测试,该技术通用性强,可以简单方便实现测试,而且不需拆装转台。
附图说明
图1为本发明的转台结构示意图。
具体实施方式
下面通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
如图1所示,三只光纤陀螺的敏感轴两两之间相互正交安装在结构体上,组成一个可以敏感三个轴向角速率信息的陀螺组件,将陀螺组件安装到待测试转台上。建立陀螺组件和转台的数学模型,根据转台非正交度误差数学模型,设计旋转顺序,控制转台旋转,采集三只陀螺输出的角速率信息,利用采集到的陀螺角速率信息,计算出转台转轴之间的非正交度。
一种基于光纤陀螺的转台非正交度测试方法,其特征在于:所述测试方法的步骤为:
1)转台结构与陀螺组件建模:
1、转台模型:建立转台模型并将转台结构分为外框方位、中框俯仰及内框横滚,并根据上述结构分别建立转台外框坐标系T1,中框坐标系T2,内框坐标系T3,陀螺三轴坐标系g系,其中g系是非正交坐标系;
从T1向T2转换,存在存在中框转轴与理想中框轴的夹角θy1和θz1,以及中框本身构建的水平面与外框旋转平面关于中框转轴的旋转角度,设为α0,那么从T1系向T2系转换如下:
Figure BDA0002282091040000061
从T2向T3转换,存在存在内框轴与理想内框轴的夹角θx2和θz2,以及内框本身构建的水平面与中框平面关于内框转轴的旋转角度,设为β0,那么从T2系向T3系转换如下:
Figure BDA0002282091040000062
2、陀螺组件模型:
单只陀螺的参数包括标度因数、零位偏置及敏感轴偏置角度,由三只陀螺构成的陀螺组件模型参数为12个,包括陀螺标度因数Kgx,Kgy,Kgz,陀螺安装偏角Ugxy,Ugxz,Ugyx,Ugyz,Ugzx,Ugzy,陀螺零偏Bgx0,Bgy0,Bgz0;
陀螺组件的数学模型如下:
Figure BDA0002282091040000063
其中:ωxyz分别为三个相互正交坐标轴上的角速度输入;
Ngx,Ngy,Ngz分别为三个陀螺脉冲数输出;
2)转台旋转次序与模型
将步骤1)中的陀螺组件安装到三轴转台上,预先定义的X,Y和Z轴分别和转台的中框、内框和外框平行,旋转顺序依次为:
1、转台外框、中框、内框归零,Z轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000071
2、中框转至90°,外框、内框归零,Y轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000072
3、内框转至270°,外框、中框归零,X轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000073
4、中框转至180°,外框、内框归零,Z轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000074
5、中框转至270°,外框、内框归零,Y轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000081
6、内框转至90°,外框、中框归零,X轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure BDA0002282091040000082
3)转台非正交度计算
结合步骤1)、2),将求解出的三个陀螺标度因数带入可求出如下变量,见表1:
表1
各轴指天 各轴指地
U<sub>gxy</sub>+θ<sub>z2</sub>-θ<sub>y1</sub> U<sub>gxy</sub>+θ<sub>z2</sub>+θ<sub>y1</sub>
U<sub>gxz</sub>-β<sub>0</sub>-θ<sub>y1</sub> U<sub>gxz</sub>-β<sub>0</sub>+θ<sub>y1</sub>
U<sub>gyx</sub>+θ<sub>x2</sub>+α<sub>0</sub> U<sub>gyx</sub>-θ<sub>x2</sub>-α<sub>0</sub>
U<sub>gyz</sub>+θ<sub>x2</sub>+α<sub>0</sub> U<sub>gyz</sub>+θ<sub>x2</sub>+α<sub>0</sub>
U<sub>gzx</sub>+β<sub>0</sub>+θ<sub>y1</sub> U<sub>gzx</sub>+β<sub>0</sub>+θ<sub>y1</sub>
U<sub>gzy</sub>-θ<sub>x2</sub>-α<sub>0</sub> U<sub>gzy</sub>-θ<sub>x2</sub>-α<sub>0</sub>
由此,可以进一步确定以下变量的值:
Ugxyz2,θy1,Ugxz0,Ugyx,θx20,Ugyz,Ugzx0,Ugzy
同时,由于三个坐标轴的两两之间不正交度的关系,只有三个自由度,因此还具备如下约束条件:
Ugxy+Ugyx=A
Ugxz+Ugzx=B
Ugyz+Ugzy=C
其中:A、B、C对于同一套陀螺组件是固定值,可以事先标定精确得到,进一步可以确认变量θz2的值,因此可计算出转台转轴之间的不正交角为:θy1,θx20,θz2,其中θx2与α0耦合。
尽管为说明目的公开了本发明的实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此,本发明的范围不局限于实施例和附图所公开的内容。

Claims (1)

1.一种基于光纤陀螺的转台非正交度测试方法,其特征在于:所述测试方法的步骤为:
1)转台结构与陀螺组件建模:
1、转台模型:建立转台模型并将转台结构分为外框方位、中框俯仰及内框横滚,并根据上述结构分别建立转台外框坐标系T1,中框坐标系T2,内框坐标系T3,陀螺三轴坐标系g系,其中g系是非正交坐标系;
从T1向T2转换,存在中框转轴与理想中框轴的夹角θy1和θz1,以及中框本身构建的水平面与外框旋转平面关于中框转轴的旋转角度,设为α0,那么从T1系向T2系转换如下:
Figure FDA0003014914520000011
从T2向T3转换,存在内框轴与理想内框轴的夹角θx2和θz2,以及内框本身构建的水平面与中框平面关于内框转轴的旋转角度,设为β0,那么从T2系向T3系转换如下:
Figure FDA0003014914520000012
2、陀螺组件模型:
单只陀螺的参数包括标度因数、零位偏置及敏感轴偏置角度,由三只陀螺构成的陀螺组件模型参数为12个,包括陀螺标度因数Kgx,Kgy,Kgz,陀螺安装偏角Ugxy,Ugxz,Ugyx,Ugyz,Ugzx,Ugzy,陀螺零偏Bgx0,Bgy0,Bgz0;
陀螺组件的数学模型如下:
Figure FDA0003014914520000013
其中:ωxyz分别为三个相互正交坐标轴上的角速度输入;
Ngx,Ngy,Ngz分别为三个陀螺脉冲数输出;
2)转台旋转次序与模型
将步骤1)中的陀螺组件安装到三轴转台上,预先定义的X,Y和Z轴分别和转台的中框、内框和外框平行,旋转顺序依次为:
1、转台外框、中框、内框归零,Z轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000021
2、中框转至90°,外框、内框归零,Y轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000022
3、内框转至270°,外框、中框归零,X轴指天,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000023
4、中框转至180°,外框、内框归零,Z轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000031
5、中框转至270°,外框、内框归零,Y轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000032
6、内框转至90°,外框、中框归零,X轴指地,外框按照角速率±ω转动整圈,对应的陀螺组件输出为:
Figure FDA0003014914520000033
3)转台非正交度计算
结合步骤1)、2),将求解出的三个陀螺标度因数带入可求出如下变量,见表1:
表1
Figure FDA0003014914520000034
Figure FDA0003014914520000041
由此,可以进一步确定以下变量的值:
Ugxyz2,θy1,Ugxz0,Ugyx,θx20,Ugyz,Ugzx0,Ugzy
同时,由于三个坐标轴的两两之间不正交度的关系,只有三个自由度,因此还具备如下约束条件:
Ugxy+Ugyx=A
Ugxz+Ugzx=B
Ugyz+Ugzy=C
其中:A、B、C对于同一套陀螺组件是固定值,可以事先标定精确得到,进一步可以确认变量θz2的值,因此可计算出转台转轴之间的不正交角为:θy1,θx20,θz2,其中θx2与α0耦合。
CN201911145605.9A 2019-11-21 2019-11-21 一种基于光纤陀螺的转台非正交度测试方法 Active CN110849294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911145605.9A CN110849294B (zh) 2019-11-21 2019-11-21 一种基于光纤陀螺的转台非正交度测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911145605.9A CN110849294B (zh) 2019-11-21 2019-11-21 一种基于光纤陀螺的转台非正交度测试方法

Publications (2)

Publication Number Publication Date
CN110849294A CN110849294A (zh) 2020-02-28
CN110849294B true CN110849294B (zh) 2021-07-27

Family

ID=69603070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911145605.9A Active CN110849294B (zh) 2019-11-21 2019-11-21 一种基于光纤陀螺的转台非正交度测试方法

Country Status (1)

Country Link
CN (1) CN110849294B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111998776B (zh) * 2020-08-27 2022-02-22 中国科学院长春光学精密机械与物理研究所 二维平台正交度检测装置及其检测方法
CN112344964B (zh) * 2020-11-13 2023-03-07 中国船舶重工集团公司第七0七研究所 一种捷联惯导系统载体航迹仿真设计方法
CN116026370B (zh) * 2023-03-30 2023-06-09 中国船舶集团有限公司第七〇七研究所 基于矩阵等价转换的光纤陀螺误差校准方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103196462A (zh) * 2013-02-28 2013-07-10 南京航空航天大学 一种mimu中mems陀螺仪的误差标定补偿方法
CN103234560A (zh) * 2013-05-17 2013-08-07 哈尔滨工程大学 一种捷联惯导系统零位标定的方法
CN103591949A (zh) * 2013-10-25 2014-02-19 西安石油大学 三轴姿态测量系统非正交性误差的正交补偿方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7726188B2 (en) * 2006-04-24 2010-06-01 Millisensor Systems + Actuators Scale factor measurement for mems gyroscopes and accelerometers

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103196462A (zh) * 2013-02-28 2013-07-10 南京航空航天大学 一种mimu中mems陀螺仪的误差标定补偿方法
CN103234560A (zh) * 2013-05-17 2013-08-07 哈尔滨工程大学 一种捷联惯导系统零位标定的方法
CN103591949A (zh) * 2013-10-25 2014-02-19 西安石油大学 三轴姿态测量系统非正交性误差的正交补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A new multi-position calibration method for MEMS intertial navigation systems;Z F Syed 等;《MEASUREMENT SCIENCE AND TECHNOLOGY》;20071231;第1897-1907页 *
低成本三轴微数字陀螺仪标定方法及实验研究;陈瑜 等;《仪表技术与传感器》;20141231;第48-51页 *

Also Published As

Publication number Publication date
CN110849294A (zh) 2020-02-28

Similar Documents

Publication Publication Date Title
CN110849294B (zh) 一种基于光纤陀螺的转台非正交度测试方法
JP6613236B2 (ja) 高精度な光ファイバジャイロ慣性測量装置の標定方法
CN103808331B (zh) 一种mems三轴陀螺仪误差标定方法
CN109974749B (zh) 一种三轴转台综合指向误差的评定方法
CN102735232B (zh) 光纤陀螺组合体惯性测量装置及其标定方法
CN109459054B (zh) 一种基于自准直跟踪的动基座姿态校准方法
CN106969783B (zh) 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术
CN101975872B (zh) 石英挠性加速度计组件零位偏置的标定方法
CN110361031B (zh) 一种基于回溯理论的imu全参数误差快速标定方法
CN108168575A (zh) 一种五轴冗余配置十表光纤惯组的标定方法和系统
CN101629830A (zh) 一种三轴一体高精度光纤陀螺的标定方法及装置
CN102636183A (zh) 基于光纤监测和双轴转台离心机的挠性陀螺二次过载项测试方法
CN105509768A (zh) 一种单轴寻北仪的误差标定方法
CN103940445B (zh) 一种单轴旋转惯导系统惯性器件误差补偿方法
CN108318051B (zh) 一种闭环光纤陀螺阈值测试方法
CN104006828B (zh) 一种惯性器件轴系误差标定方法
CN104034347B (zh) 一种星用半球谐振陀螺组合指标体系测量方法
CN111220817B (zh) 一种捷联惯组三轴加速度计的标定方法
CN110954131B (zh) 一种光纤陀螺输入轴失准角标定工装
CN104792339B (zh) 光纤陀螺测斜仪的校准方法
CN113899324B (zh) 基于单轴激光陀螺测角仪的多轴转台垂直度误差检测方法
CN115574797A (zh) 一种基于转台位置模式的角速度测量误差的测试方法
Shun-qing et al. Impacts of installation errors on the calibration accuracy of gyro accelerometer tested on centrifuge
CN114858187A (zh) 一种mems惯导系统标定方法
CN116989826B (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