CN101819284A - 一种利用陀螺仪测定地球极移的方法 - Google Patents

一种利用陀螺仪测定地球极移的方法 Download PDF

Info

Publication number
CN101819284A
CN101819284A CN 201010174339 CN201010174339A CN101819284A CN 101819284 A CN101819284 A CN 101819284A CN 201010174339 CN201010174339 CN 201010174339 CN 201010174339 A CN201010174339 A CN 201010174339A CN 101819284 A CN101819284 A CN 101819284A
Authority
CN
China
Prior art keywords
formula
moment
gyroscope
meridian
latitude
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
CN 201010174339
Other languages
English (en)
Other versions
CN101819284B (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.)
Changan University
Original Assignee
Changan 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 Changan University filed Critical Changan University
Priority to CN2010101743395A priority Critical patent/CN101819284B/zh
Publication of CN101819284A publication Critical patent/CN101819284A/zh
Application granted granted Critical
Publication of CN101819284B publication Critical patent/CN101819284B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Gyroscopes (AREA)

Abstract

本发明公开了一种利用陀螺仪测定地球极移的方法,包括以下步骤:将陀螺仪安置于已知点,使陀螺旋转轴的方向指向北方向,偏差小于10°;启动陀螺仪测量此时指向力矩Mα0,计算此时陀螺旋转轴与子午线方向的夹角α0;启动陀螺仪测量此时的指向力矩
Figure 201010174339.5_AB_0
,将α0作为已知值,计算此时陀螺仪所在位置的纬度
Figure 201010174339.5_AB_1
;将
Figure 201010174339.5_AB_1
作为已知值,再次启动陀螺仪测量指向力矩Mα1,计算此时陀螺仪所在位置的纬度;如此反复进行测量指向力矩M,并计算相应的α,
Figure 201010174339.5_AB_2
的值,然后将每一次得到的α,
Figure 201010174339.5_AB_2
值作为下一次计算的已知值进行迭代,计算得到αi:选取任意两个时刻的结果
Figure 201010174339.5_AB_4

Description

一种利用陀螺仪测定地球极移的方法
技术领域
本发明属于地球物理、大地测量学、惯性导航技术领域,具体涉及一种测定地球极移参数的方法,特别是一种利用陀螺仪测定地球极移参数的方法。本发明可以广泛应用于航空、航天、全球板块运动监测、深空探测研究等领域。
背景技术
地球本身由于其内部物质流分布的不均匀造成地球自转轴的位置不停的发生着微小的变化。这种变化被称为“极移”。准确的测量地球的极移变化对于航空、航天、地球物理学的研究都有着重要的实际意义和应用价值。
自20世纪70年代以来,人们采用多种观测技术对地球的极移进行测量,如天文观测方法、VLBI技术以及GPS技术等,而这些方法和技术大多需要大范围的地空联测,所需的设备造价昂贵,且获取观测数据的周期较长,数据后处理过程复杂。
发明内容
针对目前采用的测量地球极移的多种技术的缺陷或不足,本发明的目的在于,提供一种利用陀螺仪测定地球极移的方法,本发明的方法所需设备简单易行,不需要进行大范围的地空联测,可以独立的进行地球极移参数的测量,后期数据处理速度快,可实时监测地球极移参数的变化。
为了实现上述任务,本发明采用如下的技术解决方案:
一种利用陀螺仪测定地球极移的方法,包括以下步骤:
第一步,将陀螺仪安置于已知点,使陀螺旋转轴的方向指向北方向,偏差小于10°;
第二步,启动陀螺仪测量此时的指向力矩Mα0,将已知点的纬度
Figure GDA0000021456440000011
作为已知值,根据公式8计算得到此时陀螺旋转轴与子午线方向的夹角α0
Figure GDA0000021456440000021
(公式8)
其中,Mα0为指向力矩;H为陀螺角动量,H为已知量;ωe为地球自转角速度,ωe为已知量;
Figure GDA0000021456440000022
为陀螺安置点纬度;α0为陀螺旋转轴与子午线的夹角。
第三步,再次启动陀螺仪测量此时的指向力矩
Figure GDA0000021456440000023
将第二步得到的α0作为已知值,根据公式9计算得到此时陀螺仪所在位置的纬度
Figure GDA0000021456440000025
(公式9);
第四步,再将
Figure GDA0000021456440000026
作为已知值,再次启动陀螺仪测量指向力矩Mα1,根据公式10计算得到陀螺仪所在位置的纬度α1
Figure GDA0000021456440000027
(公式10);
第五步,如此反复进行测量指向力矩M,并根据公式2、公式3计算相应的α,
Figure GDA0000021456440000028
的值,然后将每一次得到的α,
Figure GDA0000021456440000029
值作为下一次计算的已知值进行迭代,根据公式7计算得到αi
Figure GDA00000214564400000210
Figure GDA00000214564400000211
(公式7)
即可得到一组测量成果
Figure GDA00000214564400000212
第六步,选取任意两个时刻的结果
Figure GDA00000214564400000213
根据公式11、公式12和公式6即可得到i=m时刻到i=n时刻之间地球的极移量
Figure GDA00000214564400000215
Figure GDA00000214564400000216
(公式11),
Figure GDA00000214564400000217
(公式12),
AD → = AB → + AC → (公式6)。
本发明的理论根据如下:
高速旋转的陀螺,在地球自转与自身重力矩的共同影响下会产生向子午线方向进动的指向力矩,指向力矩的大小为:
Figure GDA0000021456440000031
(公式1)
其中,M为指向力矩;H为陀螺角动量,为已知量;ωe为地球自转角速度,为已知量;为陀螺仪所在位置的纬度;α为陀螺旋转轴与子午线的夹角。
在指向力矩的作用下,陀螺旋转轴会绕着子午线往复摆动,在摆动的过程中,随着陀螺旋转轴与子午线夹角α的变化,指向力矩M的大小也在不断的发生变化。如果通过施加与指向力矩大小相等,方向相反的外加力矩使陀螺旋转轴处于静止状态,并通过测量外加力矩的大小测定指向力矩的大小。再根据公式1,可得计算陀螺仪所在位置的纬度公式为:
(公式2)
计算陀螺旋转轴与子午线夹角α的公式为:
Figure GDA0000021456440000035
(公式3)
如图1所示,地球半径为R时,假设地轴由AA′移动到DD′,现将陀螺仪安置在地面上的任意点T,当地轴位于AA′时,利用陀螺安置点的纬度为
Figure GDA0000021456440000036
陀螺旋转轴与子午线的夹角为αA,当地轴位于DD′时,利用陀螺安置点的纬度为陀螺旋转轴与子午线的夹角为αD
现将地轴由A→D的位移分解为
Figure GDA0000021456440000038
两个互相垂直方向上位移的合成,其中,
Figure GDA0000021456440000039
位移方向与地轴位于AA′时的T点所在子午面垂直,位移方向在地轴位于AA′时T点所在的子午面内。根据球面三角形的几何关系可得,
(公式4)
Figure GDA00000214564400000312
(公式5)
则地轴的极移量为:
AD → = AB → + AC → (公式6)
如果利用陀螺仪测定
Figure GDA0000021456440000041
αAαD,便可根据公式4、公式5和公式6计算地轴在两个位置的极移量。
由于地球的极移是一个缓慢微小的变化过程(每年极移量变化约20米),而且指向力矩的测量可以在短时间内(15分钟)完成,那么纬度
Figure GDA0000021456440000043
与夹角α的测量计算就可以通过相互迭代的方法来得到。并且可以通过天文测量或GPS测量的方法,测定陀螺仪所在位置纬度作为迭代计算的初始值。
通过常年累月的测量即可获取一系列
Figure GDA0000021456440000044
α0
Figure GDA0000021456440000045
α1,…αi,…,通过这些数据即可模拟出地球自转轴的变化过程。假设如果在2006年7月1日12时获得的测量的成果为在2007年3月1日12时获得的测量成果为根据公式4、公式5、公式6即可计算出2006年7月1日到2007年3月1日之间地球自转轴的极移量。
附图说明
图1为地球极移分解示意图;
图2为本发明的实施例所采用的陀螺仪的结构示意图;
以下结合附图与具体实施例对本发明做进一步的解释说明。
具体实施方式
本实施例所用陀螺仪结构如图2所示:包括外部壳体6以及位于外部壳体6内部的磁悬浮系统,外部壳体6底部中心位置设有对中标识23,外部壳体6上方设有与外部计算机系统相连接的数据接口26和北向标识27。
磁悬浮系统包括回转壳体7、力矩器壳体19、第一电感线圈8、第二电感线圈9、弹簧10、压片11、反射棱镜组15、力矩器定子21和陀螺灵敏部;该陀螺灵敏部包括:磁浮球12、连接杆14、光电传感器16、力矩器转子21、陀螺马达房17以及位于陀螺马达房17内的陀螺马达18和陀螺旋转轴28;回转壳体7为一门形;力矩器壳体19为一倒门形,回转壳体7与力矩器壳体19固连并位于力矩器壳体19上方;回转壳体7内壁顶部中心设有第一电感线圈8;回转壳体7内壁顶部围绕第一电感线圈8的八个方向均匀设有8个第二电感线圈9;每个第二电感线圈9下方固连一个弹簧10,弹簧10下方固连压片11;力矩器壳体19内壁上方设有两个对称的反射棱镜组15;力矩器壳体19内壁下方设有力矩器定子21;磁浮球12位于第一电感线圈8下方,连接杆14位于磁浮球12下方与磁浮球12固连,连接杆14包括主杆体140和垂直于主杆体140的平面141,平面141上设有与8个与第二电感线圈9位置相对应的触头142,平面141下方设有两个光电传感器16与连接杆14固连,连接杆14位于陀螺房17上方,与陀螺房17固连,陀螺房17内设有陀螺马达18,陀螺旋转轴28穿过陀螺马达18中心,且与陀螺房17固连;陀螺房17下方设有力矩器转子21。
上述陀螺仪的工作过程如下:启动陀螺马达,待陀螺马达转速达到额定转速后,第二电感线圈9通电,压片11由于磁力作用被吸附起来,弹簧10处于压缩状态,然后第一电感线圈8通电,由于磁力作用,磁浮球12被向上拉起,在连接杆14的作用下,陀螺灵敏部处于悬浮状态,力矩器定子20与力矩器转子21之间形成水平磁场,在该磁场力的作用下陀螺灵敏部保持静止,此时光电传感器16的发射/接收面与反射棱镜组15正对,光电传感器16向反射棱镜组15发射光束,光束到达反射棱镜组15后被反射回来,当发射光束与反射光束稳定重合时,力矩器定子20与力矩器转子21开始采集测量数据,并通过数据接口26传输给外部计算机,外部计算机计算力矩值。
本发明的利用陀螺仪测定地球极移的方法,包括以下步骤:
1)将陀螺仪安置于已知点,已知点的经纬度为已知值,使陀螺仪的北向标识27指向北方向±0°;
2)启动陀螺仪测量此时的指向力矩Mα0,将已知点的纬度
Figure GDA0000021456440000051
作为已知值,根据公式8计算陀螺旋转轴28与子午线方向的夹角α0
Figure GDA0000021456440000052
(公式8)
其中,Mα0为指向力矩;H为陀螺角动量,H为已知量;ωe为地球自转角速度,ωe为已知量;为陀螺安置点纬度;α0为陀螺旋转轴28与子午线的夹角。
3)再次启动陀螺仪测量此时的指向力矩
Figure GDA0000021456440000062
将第二步得到的α0作为已知值,根据公式9计算此时陀螺仪所在位置的纬度
Figure GDA0000021456440000063
Figure GDA0000021456440000064
(公式9);
第四步,再将
Figure GDA0000021456440000065
作为已知值,再次启动陀螺仪测量指向力矩Mα1,根据公式10计算此时陀螺旋转轴28与子午线的夹角α1
Figure GDA0000021456440000066
(公式10);
第五步,如此反复进行测量指向力矩M,并根据公式2、公式3计算相应的α,
Figure GDA0000021456440000067
的值,然后将每一次得到的α,
Figure GDA0000021456440000068
值作为下一次计算的已知值进行迭代,根据公式7计算得到αi
Figure GDA0000021456440000069
Figure GDA00000214564400000610
(公式7),
即可得到一组测量成果
Figure GDA00000214564400000611
α0
Figure GDA00000214564400000612
α1,…
Figure GDA00000214564400000613
αi,…;
第六步,选取任意两个时刻的结果
Figure GDA00000214564400000615
根据公式11、公式12和公式6即可得到i=m时刻到i=n时刻之间地球的极移量
Figure GDA00000214564400000616
Figure GDA00000214564400000617
(公式11),
Figure GDA00000214564400000618
(公式12),
AD → = AB → + AC → (公式6)。

Claims (1)

1.一种利用陀螺仪测定地球极移的方法,其特征在于,包括以下步骤:
第一步,将陀螺仪安置于已知点,使陀螺旋转轴的方向指向北方向,偏差小于10°;
第二步,启动陀螺仪测量此时的指向力矩Mα0,将已知点的纬度
Figure FDA0000021456430000011
作为已知值,根据公式8计算得到此时陀螺旋转轴与子午线方向的夹角α0
(公式8)
其中,Mα0为指向力矩;H为陀螺角动量,H为已知量;ωα为地球自转角速度,ωα为已知量;
Figure FDA0000021456430000013
为陀螺安置点纬度;α0为陀螺旋转轴与子午线的夹角。
第三步,再次启动陀螺仪测量此时的指向力矩
Figure FDA0000021456430000014
将第二步得到的α0作为已知值,根据公式9计算得到此时陀螺仪所在位置的纬度
Figure FDA0000021456430000015
Figure FDA0000021456430000016
(公式9);
第四步,再将
Figure FDA0000021456430000017
作为已知值,再次启动陀螺仪测量指向力矩Mα1,根据公式10计算得到此时陀螺仪所在位置的纬度α1
Figure FDA0000021456430000018
(公式10);
第五步,如此反复进行测量指向力矩M,并根据公式2、公式3计算相应的α,
Figure FDA0000021456430000019
的值,然后将每一次得到的α,
Figure FDA00000214564300000110
值作为下一次计算的已知值进行迭代,根据公式7计算得到αi
Figure FDA00000214564300000111
Figure FDA00000214564300000112
(公式7),
即可得到一组测量成果
Figure FDA00000214564300000113
α0
Figure FDA00000214564300000114
α1,…
Figure FDA00000214564300000115
αi,…;
第六步,选取任意两个时刻的结果
Figure FDA00000214564300000116
Figure FDA00000214564300000117
根据公式11、公式12和公式6即可得到i=m时刻到i=n时刻之间地球的极移量
Figure FDA00000214564300000118
Figure FDA00000214564300000119
(公式11),
Figure FDA0000021456430000021
(公式12),
AD ‾ = AB ‾ + AC ‾ (公式6)。
CN2010101743395A 2010-05-18 2010-05-18 一种利用陀螺仪测定地球极移的方法 Expired - Fee Related CN101819284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101743395A CN101819284B (zh) 2010-05-18 2010-05-18 一种利用陀螺仪测定地球极移的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101743395A CN101819284B (zh) 2010-05-18 2010-05-18 一种利用陀螺仪测定地球极移的方法

Publications (2)

Publication Number Publication Date
CN101819284A true CN101819284A (zh) 2010-09-01
CN101819284B CN101819284B (zh) 2012-07-04

Family

ID=42654452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101743395A Expired - Fee Related CN101819284B (zh) 2010-05-18 2010-05-18 一种利用陀螺仪测定地球极移的方法

Country Status (1)

Country Link
CN (1) CN101819284B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635200A (zh) * 2015-01-30 2015-05-20 长安大学 基于陀螺定向的标校相控阵天线北方向的测量方法
CN110108265A (zh) * 2019-04-25 2019-08-09 长安大学 一种自主获取地理纬度并自动寻北的陀螺测量仪器

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4449817A (en) * 1981-03-02 1984-05-22 Abler William L Celestial navigator
US4691798A (en) * 1985-10-10 1987-09-08 Engelbach Thomas S System for automatically preventing turnover of two wheeled vehicles
CN1052371A (zh) * 1990-05-23 1991-06-19 航空航天部第一设计研究院十五所 陀螺寻北仪
FR2659061A1 (fr) * 1990-03-05 1991-09-06 Gen Electric Procede de stabilisation d'un satellite.
WO2006115041A1 (ja) * 2005-04-25 2006-11-02 Alps Electric Co., Ltd. 傾斜センサおよびこれを用いた方位計測装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4449817A (en) * 1981-03-02 1984-05-22 Abler William L Celestial navigator
US4691798A (en) * 1985-10-10 1987-09-08 Engelbach Thomas S System for automatically preventing turnover of two wheeled vehicles
FR2659061A1 (fr) * 1990-03-05 1991-09-06 Gen Electric Procede de stabilisation d'un satellite.
CN1052371A (zh) * 1990-05-23 1991-06-19 航空航天部第一设计研究院十五所 陀螺寻北仪
WO2006115041A1 (ja) * 2005-04-25 2006-11-02 Alps Electric Co., Ltd. 傾斜センサおよびこれを用いた方位計測装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《2009全国时间频率学术会议论文集》 20091231 何战科;杨旭海等 利用GPS资料解算地球自转参数 全文 1 , *
《矿山测量》 19841231 周维生 极移及其对地面方位的影响极移及其对地面方位的影响 全文 1 , *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635200A (zh) * 2015-01-30 2015-05-20 长安大学 基于陀螺定向的标校相控阵天线北方向的测量方法
CN104635200B (zh) * 2015-01-30 2017-04-26 长安大学 基于陀螺定向的标校相控阵天线北方向的测量方法
CN110108265A (zh) * 2019-04-25 2019-08-09 长安大学 一种自主获取地理纬度并自动寻北的陀螺测量仪器
CN110108265B (zh) * 2019-04-25 2022-11-18 长安大学 一种自主获取地理纬度并自动寻北的陀螺测量仪器

Also Published As

Publication number Publication date
CN101819284B (zh) 2012-07-04

Similar Documents

Publication Publication Date Title
CN101344391B (zh) 基于全功能太阳罗盘的月球车位姿自主确定方法
CN101887068B (zh) 三轴矢量传感器及两轴矢量传感器的标定补偿方法
CN104374388B (zh) 一种基于偏振光传感器的航姿测定方法
CN104880191B (zh) 一种基于太阳矢量的偏振辅助导航方法
CN103852085B (zh) 一种基于最小二乘拟合的光纤捷联惯导系统现场标定方法
CN102506898A (zh) 一种基于遗传算法的惯性/地磁传感器标定方法
CN102278108B (zh) 小口径定向陀螺测斜仪连续测量模式标定方法
CN103090870A (zh) 一种基于mems传感器的航天器姿态测量方法
CN102071924B (zh) 一种陀螺测斜仪全方位连续测量方法
CN103308045A (zh) 移动物体在三维空间中的实时位置定位系统和定位方法
CN102288201B (zh) 用于星敏感器的精度测量方法
CN110108300A (zh) 一种基于卧式三轴转台的imu正六面体标定方法
CN105737848B (zh) 一种系统级星敏感器观星系统及观星方法
CN103743413A (zh) 倾斜状态下调制寻北仪安装误差在线估计与寻北误差补偿方法
CN103115623A (zh) 基于偏振光仿生导航的定位系统及其定位方法
CN106017452A (zh) 双陀螺抗扰动寻北方法
CN110285810A (zh) 一种基于惯性导航数据的采煤机自主定位方法及装置
CN105674982A (zh) 一种六参数量子惯性传感器及其测量方法
CN106842080A (zh) 一种磁场测量装置姿态摆动干扰去除方法
CN103091723B (zh) 降低重力卫星质心调整误差对地球重力场精度影响的方法
CN101819284B (zh) 一种利用陀螺仪测定地球极移的方法
CN110986961A (zh) 磁光矩阵标定装置及磁光矩阵标定方法
CN104655123B (zh) 一种利用光纤陀螺测定地球自转角速度的方法
CN103983259A (zh) 基于摆式陀螺最大进动速度粗寻北原理的全方位寻北方法
CN102494679B (zh) 一种超导磁悬浮的寻北装置

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120704