CN104504240A - 航天器总装精度测量计算方法 - Google Patents

航天器总装精度测量计算方法 Download PDF

Info

Publication number
CN104504240A
CN104504240A CN201410710525.4A CN201410710525A CN104504240A CN 104504240 A CN104504240 A CN 104504240A CN 201410710525 A CN201410710525 A CN 201410710525A CN 104504240 A CN104504240 A CN 104504240A
Authority
CN
China
Prior art keywords
coordinate system
transit
prism square
computing method
spacecraft
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
CN201410710525.4A
Other languages
English (en)
Other versions
CN104504240B (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.)
Shanghai Institute of Satellite Equipment
Original Assignee
Shanghai Institute of Satellite Equipment
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 Shanghai Institute of Satellite Equipment filed Critical Shanghai Institute of Satellite Equipment
Priority to CN201410710525.4A priority Critical patent/CN104504240B/zh
Publication of CN104504240A publication Critical patent/CN104504240A/zh
Application granted granted Critical
Publication of CN104504240B publication Critical patent/CN104504240B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提供了一种航天器总装精度测量计算方法,其包括以下步骤:步骤一:定义经纬仪测量坐标系、各立方镜坐标系;步骤二:确定航天器中立方镜方位的测量过程;步骤三:使用单台经纬仪坐标系变换原理进行准直镜面矢量的方位计算;步骤四:使用两台经纬仪坐标系变换原理进行互瞄方位传递的计算;步骤五:利用坐标系变换计算立方镜之间的方位关系。本发明根据经纬仪的测量原理构建测站坐标系,并将立方镜的法线矢量转移到同一测站坐标系下,最终完成各个矢量关系的解算,简便且不容易出错。

Description

航天器总装精度测量计算方法
技术领域
本发明涉及一种计算方法,具体地,涉及一种航天器总装精度测量计算方法。
背景技术
在航天器总装过程中,通常需要测量有效载荷、单机(如星敏、陀螺、飞轮等)与航天器基准之间以及有关载荷与载荷、单机之间的方位关系。由于航天器需要测量的待测部件往往空间狭小,结果复杂,因此普遍的做法是在这些单机或载荷上安置固定一些立方镜,以此将其姿态或光学路径转换到立方镜上。
立方镜是由6个平面度较高的镜面组成的正立方体,其中5个为工作面,一搬采用石英晶体或金属制作,表面镀有反射膜,如图1所示。其相邻两个面有着较高的垂直度,垂直度加工误差约为1″~3″。利用立方镜的三个相邻面的法线就可以构成一个纯角度坐标系,因此对部件的姿态测量就可以转换到对立方镜姿态的测量。
目前立方镜之间的方位关系都是通过经纬仪建站测量实施的。实际测量时,使用3台及以上经纬仪进行准直、互瞄等过程,将被测对象上光学立方镜反射面的法线联系起来,以此计算它们之间的角度关系,如解放军信息工程大学测绘学院的杨振等刊登的《准直测量进行高精度立方镜间关系标定》。但目前采用的计算方法繁杂且容易出错,应用范围很局限。因此有必要应用一种新的航天器总装精度测量计算方法。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种航天器总装精度测量计算方法,其根据经纬仪的测量原理构建测站坐标系,并将立方镜的法线矢量转移到同一测站坐标系下,最终完成各个矢量关系的解算,简便且不容易出错。
根据本发明的一个方面,提供一种航天器总装精度测量计算方法,其特征在于,所述航天器总装精度测量计算方法包括以下步骤:
步骤一:定义经纬仪测量坐标系、各立方镜坐标系;
步骤二:确定航天器中立方镜方位的测量过程;
步骤三:使用单台经纬仪坐标系变换原理进行准直镜面矢量的方位计算;
步骤四:使用两台经纬仪坐标系变换原理进行互瞄方位传递的计算;
步骤五:利用坐标系变换计算立方镜之间的方位关系。
优选地,所述航天器总装精度测量计算方法结合经纬仪测量设备的特点建立测量坐标系;利用立方镜的形状特点建立立方镜坐标系;经纬仪准直测量立方镜时镜面法线矢量的表示;经纬仪互瞄测量时其各自测量坐标系的传递;最后将各经纬仪测量坐标系下的镜面矢量统一到一个测量坐标系下描述,实现立方镜的方位关系计算。
优选地,所述步骤三的单台经纬仪坐标系变换原理如下:经纬仪在准直时的坐标系定义为准直测站坐标系{A},当其旋转θ角时将产生一个新的坐标系{B},两个坐标系之间的变换关系为如下式:
R B A = f ( θ ) = cos ( θ ) - sin ( θ ) 0 sin ( θ ) cos ( θ ) 0 0 0 1 .
优选地,所述步骤四的两台经纬仪坐标系变换原理如下:两台经纬仪互瞄时将要进行测站坐标系的传递,按坐标系转换原理可以知道变换关系为如下式:
优选地,所述经纬仪是一种基于大地水平和自身调平进行测量的仪器。
优选地,所述步骤五采用四台经纬仪分别准直两个立方镜的各自相邻两个镜面。
优选地,所述两个立方镜为基准立方镜和单机立方镜。
与现有技术相比,本发明具有如下的有益效果:本发明根据经纬仪的测量原理构建测站坐标系,并将立方镜的法线矢量转移到同一测站坐标系下,最终完成各个矢量关系的解算,简便且不容易出错,应用范围广。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明中经纬仪(以Leica TM5100A为例)测量坐标系定义的示意图。
图2为本发明中单台经纬仪测站坐标系旋转示意图。
图3为本发明中两台经纬仪互瞄时测站坐标系关系示意图;
图4为本发明中经纬仪准直立方镜时镜面矢量关系示意图;
图5为本发明中航天器总装精度测量示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,本发明航天器总装精度测量计算方法包括以下步骤:
步骤一:定义经纬仪测量坐标系、各立方镜坐标系;首先定义经纬仪测站坐标系(本发明涉及到的坐标系只有角度关系,不涉及位置关系):经纬仪是一种基于大地水平和自身调平进行测量的仪器。测量时,所有仪器的竖直轴互相平行,并垂直于大地水平面,可定义其测站坐标系为如图1所示,竖直向上为+Z轴,水平度盘的任意位置为+X轴(即视准轴在水平面的投影方向),第三个方向由右手法则确定。
步骤二:确定航天器中立方镜方位的测量过程;
步骤三:使用单台经纬仪坐标系变换原理进行准直镜面矢量的方位计算;单台经纬仪坐标系变换原理如下:如图2所示,经纬仪在准直时的坐标系定义为准直测站坐标系{A},当其旋转θ角(按右手法则逆时针为正)时将产生一个新的坐标系{B},两个坐标系之间的变换关系为如下式(1):
R B A = f ( θ ) = cos ( θ ) - sin ( θ ) 0 sin ( θ ) cos ( θ ) 0 0 0 1 . . . ( 1 )
其中R表示坐标系{B}到{A}的角度旋转变换矩阵,f定义为θ的函数。
步骤四:使用两台经纬仪坐标系变换原理进行互瞄方位传递的计算;两台经纬仪坐标系变换原理如下:如图3所示,两台经纬仪互瞄时将要进行测站坐标系的传递,按坐标系转换原理可以知道变换关系为如下式(2):
单台经纬仪准直测量时立方镜镜面法线矢量的表示:如图4所示,当经纬仪准直立方镜镜面时,按经纬仪测站坐标系的定义,此时立方镜镜面法线矢量可以表示为如下式(3):
v i Ti = - sin ( 180 - σ i ) 0 cos ( 180 - σ i ) ′ . . . ( 3 )
其中,i表示测站号,σi为经纬仪俯仰角的度数,也即其望远镜视准轴与竖直向上方向的夹角,v表示立方镜镜面法线矢量。
步骤五:利用坐标系变换计算立方镜之间的方位关系,也即各立方镜坐标系之间的坐标转换。
两立方镜方位关系测量过程如图5所示,采用四台经纬仪(第一经纬仪T1、第二经纬仪T2、第三经纬仪T3、第四经纬仪T4)分别准直两个立方镜(基准立方镜Q1和单机立方镜Q2)的各自相邻两个镜面;然后瞄准基准立方镜的两个经纬仪(第一经纬仪T1、第二经纬仪T2)互瞄,建立第一经纬仪T1和第二经纬仪T2之间的关系;最后瞄准单机立方镜的第三经纬仪T3、第四经纬仪T4与第一经纬仪T1、第二经纬仪T2中之一进行互瞄,以便将单机立方镜镜面矢量传递到基准立方镜坐标系中。
按照图5所示测量过程,经纬仪之间坐标系转换关系计算步骤,其中表示第一经纬仪T1和第二经纬仪T2互瞄时,第一经纬仪T1所在坐标系到基准立方镜时的坐标系的转换关系,其余以此类推:a1,第一经纬仪T1和第二经纬仪T2准直互瞄: a2,第三经纬仪T3准直及其与第一经纬仪T1互瞄:a3,第四经纬仪T4准直及其与第二经纬仪T2互瞄:
两立方镜方位关系计算步骤如下,其中表示镜面矢量v1在第一经纬仪T1坐标系下的矢量表示:将各镜面矢量转换到测站坐标系T11下,将镜面2矢量转换到坐标系T11中:此时可以计算立方镜两个镜面的夹角关系,用于现场测量的数据检核:将镜面3矢量转换到坐标系T11中:将镜面4矢量转换到坐标系T11中: v 4 T 1 = R T 12 T 1 · R T 21 T 12 · R T 2 T 21 · R T 24 T 2 · R T 42 T 24 R T 4 T 42 · v 4 T 4 .
构建两立方镜坐标系并计算它们之间的关系:构建基准坐标系A,镜面1保持不变: v x A = v 1 T 1 , v y ′ A = v 2 T 1 , v z A = v x A × v y ′ A , v y A = v z A × v x A ; 计算坐标系A与测站坐标系T11之间的关系: R T 1 A = v x A v y A v z A - 1 ; 计算镜面矢量在A坐标系下的表达: v 3 A = R T 1 A · v 3 T 1 , v 4 A = R T 1 A · v 4 T 1 ; 构建单机坐标系B,镜面3保持不变: v x B = v 3 A , v y ′ B = v 4 A , v z B = v x B × v y ′ B , v y B = v z B × v x B ; 计算单机坐标系与基准坐标系之间的变换关系: R B A = v x B v y B v z B .
本发明结合经纬仪测量设备的特点建立测量坐标系;利用立方镜的形状特点建立立方镜坐标系;经纬仪准直测量立方镜时镜面法线矢量的表示;经纬仪互瞄测量时其各自测量坐标系的传递;最后将各经纬仪测量坐标系下的镜面矢量统一到一个测量坐标系下描述,实现立方镜的方位关系计算。本发明可用于航天器中利用经纬仪测量安装在航天器上的立方镜方位的计算,也能用于经纬仪建站测量点坐标的过程中。立方镜的坐标系定义根据实际需要进行确定,一般按右手法则,确定相邻三个垂直面的法线方向分布为+X、+Y、+Z轴,并按一定的规则进行正交化修正。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (7)

1.一种航天器总装精度测量计算方法,其特征在于,所述航天器总装精度测量计算方法包括以下步骤:
步骤一:定义经纬仪测量坐标系、各立方镜坐标系;
步骤二:确定航天器中立方镜方位的测量过程;
步骤三:使用单台经纬仪坐标系变换原理进行准直镜面矢量的方位计算;
步骤四:使用两台经纬仪坐标系变换原理进行互瞄方位传递的计算;
步骤五:利用坐标系变换计算立方镜之间的方位关系。
2.根据权利要求1所述的航天器总装精度测量计算方法,其特征在于,所述航天器总装精度测量计算方法结合经纬仪测量设备的特点建立测量坐标系;利用立方镜的形状特点建立立方镜坐标系;经纬仪准直测量立方镜时镜面法线矢量的表示;经纬仪互瞄测量时其各自测量坐标系的传递;最后将各经纬仪测量坐标系下的镜面矢量统一到一个测量坐标系下描述,实现立方镜的方位关系计算。
3.根据权利要求1所述的航天器总装精度测量计算方法,其特征在于,所述步骤三的单台经纬仪坐标系变换原理如下:经纬仪在准直时的坐标系定义为准直测站坐标系{A},当其旋转θ角时将产生一个新的坐标系{B},两个坐标系之间的变换关系为如下式:
R B A = f ( θ ) = cos ( θ ) - sin ( θ ) 0 sin ( θ ) cos ( θ ) 0 0 0 1 .
4.根据权利要求1所述的航天器总装精度测量计算方法,其特征在于,所述步骤四的两台经纬仪坐标系变换原理如下:两台经纬仪互瞄时将要进行测站坐标系的传递,按坐标系转换原理可以知道变换关系为如下式:
5.根据权利要求1所述的航天器总装精度测量计算方法,其特征在于,所述经纬仪是一种基于大地水平和自身调平进行测量的仪器。
6.根据权利要求1所述的航天器总装精度测量计算方法,其特征在于,所述步骤五采用四台经纬仪分别准直两个立方镜的各自相邻两个镜面。
7.根据权利要求6所述的航天器总装精度测量计算方法,其特征在于,所述两个立方镜为基准立方镜和单机立方镜。
CN201410710525.4A 2014-11-27 2014-11-27 航天器总装精度测量计算方法 Active CN104504240B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410710525.4A CN104504240B (zh) 2014-11-27 2014-11-27 航天器总装精度测量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410710525.4A CN104504240B (zh) 2014-11-27 2014-11-27 航天器总装精度测量计算方法

Publications (2)

Publication Number Publication Date
CN104504240A true CN104504240A (zh) 2015-04-08
CN104504240B CN104504240B (zh) 2017-10-10

Family

ID=52945637

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410710525.4A Active CN104504240B (zh) 2014-11-27 2014-11-27 航天器总装精度测量计算方法

Country Status (1)

Country Link
CN (1) CN104504240B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105953771A (zh) * 2016-06-14 2016-09-21 西安中科光电精密工程有限公司 一种主动式经纬仪系统及测量方法
CN107121123A (zh) * 2017-05-18 2017-09-01 上海卫星工程研究所 卫星精度单机测量方法
CN108426523A (zh) * 2018-03-02 2018-08-21 北京空间技术研制试验中心 一种精度基准转移方法
CN109470272A (zh) * 2018-12-05 2019-03-15 中国科学院长春光学精密机械与物理研究所 一种imu测量基准的标定方法
CN109489642A (zh) * 2018-11-09 2019-03-19 北京航天计量测试技术研究所 一种空间任意姿态下两立方镜相对姿态的动态测量方法
CN109631828A (zh) * 2019-01-23 2019-04-16 中国科学院长春光学精密机械与物理研究所 基于立方棱镜的相邻空间坐标系的坐标轴夹角的检测方法
CN111238458A (zh) * 2020-03-17 2020-06-05 散裂中子源科学中心 一种设备基准快速调平方法
CN112648956A (zh) * 2021-03-01 2021-04-13 北京卫星制造厂有限公司 一种基于联合标定的空间位姿实时测调方法
CN113358024A (zh) * 2021-05-08 2021-09-07 上海卫星工程研究所 一种卫星仪器相对精测数据处理方法及系统
CN114646261A (zh) * 2022-03-14 2022-06-21 北京卫星环境工程研究所 一种基于斜向观测镜面法方向的测量方法和系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102032918A (zh) * 2010-10-20 2011-04-27 郑州辰维科技股份有限公司 三探头星敏指向标定方法
CN103363949A (zh) * 2013-07-19 2013-10-23 北京卫星制造厂 一种卫星天线混合测量分析的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102032918A (zh) * 2010-10-20 2011-04-27 郑州辰维科技股份有限公司 三探头星敏指向标定方法
CN103363949A (zh) * 2013-07-19 2013-10-23 北京卫星制造厂 一种卫星天线混合测量分析的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
中国测绘地理信息学会仪器装备专业委员会: "《测绘地理信息仪器装备发展研究》", 1 November 2014, 测绘出版社 *
孙方金 等: "《定向原理与方位角的传递》", 30 June 2014, 中国宇航出版社 *
李广云 等: "《工业测量系统原理与应用》", 31 January 2011, 测绘出版社 *
杨再华: "航天器总装精度测量方法分析", 《航天器环境工程》 *
杨振 等: "多种互瞄条件下立方镜准直传递方法研究", 《宇航计测技术》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105953771B (zh) * 2016-06-14 2018-02-13 西安中科光电精密工程有限公司 一种主动式经纬仪系统及测量方法
CN105953771A (zh) * 2016-06-14 2016-09-21 西安中科光电精密工程有限公司 一种主动式经纬仪系统及测量方法
CN107121123A (zh) * 2017-05-18 2017-09-01 上海卫星工程研究所 卫星精度单机测量方法
CN107121123B (zh) * 2017-05-18 2019-06-28 上海卫星工程研究所 卫星精度单机测量方法
CN108426523A (zh) * 2018-03-02 2018-08-21 北京空间技术研制试验中心 一种精度基准转移方法
CN109489642A (zh) * 2018-11-09 2019-03-19 北京航天计量测试技术研究所 一种空间任意姿态下两立方镜相对姿态的动态测量方法
CN109489642B (zh) * 2018-11-09 2020-10-23 北京航天计量测试技术研究所 一种空间任意姿态下两立方镜相对姿态的动态测量方法
CN109470272A (zh) * 2018-12-05 2019-03-15 中国科学院长春光学精密机械与物理研究所 一种imu测量基准的标定方法
CN109631828B (zh) * 2019-01-23 2020-11-13 中国科学院长春光学精密机械与物理研究所 基于立方棱镜的相邻空间坐标系的坐标轴夹角的检测方法
CN109631828A (zh) * 2019-01-23 2019-04-16 中国科学院长春光学精密机械与物理研究所 基于立方棱镜的相邻空间坐标系的坐标轴夹角的检测方法
CN111238458A (zh) * 2020-03-17 2020-06-05 散裂中子源科学中心 一种设备基准快速调平方法
CN111238458B (zh) * 2020-03-17 2022-03-25 散裂中子源科学中心 一种设备基准快速调平方法
CN112648956A (zh) * 2021-03-01 2021-04-13 北京卫星制造厂有限公司 一种基于联合标定的空间位姿实时测调方法
WO2022183761A1 (zh) * 2021-03-01 2022-09-09 北京卫星制造厂有限公司 一种基于联合标定的的空间位姿实时测调方法
CN113358024A (zh) * 2021-05-08 2021-09-07 上海卫星工程研究所 一种卫星仪器相对精测数据处理方法及系统
CN113358024B (zh) * 2021-05-08 2022-08-12 上海卫星工程研究所 一种卫星仪器相对精测数据处理方法及系统
CN114646261A (zh) * 2022-03-14 2022-06-21 北京卫星环境工程研究所 一种基于斜向观测镜面法方向的测量方法和系统
CN114646261B (zh) * 2022-03-14 2023-09-19 北京卫星环境工程研究所 一种基于斜向观测镜面法方向的测量方法和系统

Also Published As

Publication number Publication date
CN104504240B (zh) 2017-10-10

Similar Documents

Publication Publication Date Title
CN104504240A (zh) 航天器总装精度测量计算方法
CN103363949B (zh) 一种卫星天线混合测量分析的方法
CN104266649B (zh) 基于陀螺经纬仪测量基准立方镜姿态角度的方法
CN105068065B (zh) 星载激光测高仪在轨检校方法及系统
CN103591949B (zh) 三轴姿态测量系统非正交性误差的正交补偿方法
CN104848833A (zh) 一种电子经纬仪和激光跟踪仪联合测量系统的建立方法
CN103323026A (zh) 星敏感器和有效载荷的姿态基准偏差估计与修正方法
CN109520486B (zh) 一种基于水平张量重力梯度的垂线偏差实时计算方法
CN104457684A (zh) 免固定设站点的全站仪三维变形监测方法
CN106125069A (zh) 一种基于指向角残差的星载激光测高仪指向角系统误差标定方法
CN103217159A (zh) 一种sins/gps/偏振光组合导航系统建模及动基座初始对准方法
CN102087117B (zh) 飞船交会对接用测距敏感器精度的地面测量方法
CN104049269B (zh) 一种基于激光测距和mems/gps组合导航系统的目标导航测绘方法
Gao et al. Development and calibration of an accurate 6-degree-of-freedom measurement system with total station
CN103737433B (zh) 大尺寸构件复合式精度测量方法
CN102207380B (zh) 一种高精度的横轴倾斜误差补偿方法
CN103913179A (zh) 舰载单轴旋转调制捷联惯导姿态角修正方法
CN102564417A (zh) 非接触式动态主动定位方法
Guo et al. Application of a self-compensation mechanism to a rotary-laser scanning measurement system
Lu et al. Improved calibration of IMU biases in analytic coarse alignment for AHRS
CN104330078B (zh) 一种基于三点后方交会模型的联合测量方法
CN105652333A (zh) 一种低低跟踪重力测量卫星四点三线模型及其建立方法
Gao et al. Integrated calibration of a 3D attitude sensor in large-scale metrology
CN106338261A (zh) 一种两束干涉仪出射平面波光束间角度偏差标定方法
Lu et al. Self-calibration of gyro asymmetry for single-axis forward–reverse rotating inertial navigation system under arbitrary attitude

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Li Qiang

Inventor after: Ling Min

Inventor after: Liu Siwei

Inventor after: Che Lamei

Inventor before: Li Qiang

Inventor before: Ling Min

Inventor before: Liu Siwei

Inventor before: Che Lamei

COR Change of bibliographic data
CB03 Change of inventor or designer information

Inventor after: Li Qiang

Inventor after: Jiang Jian

Inventor after: Ling Min

Inventor after: Liu Siwei

Inventor after: Che Lamei

Inventor before: Li Qiang

Inventor before: Ling Min

Inventor before: Liu Siwei

Inventor before: Che Lamei

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant