CN106802125B - 航天器用电推力器推力矢量偏心的确定方法 - Google Patents

航天器用电推力器推力矢量偏心的确定方法 Download PDF

Info

Publication number
CN106802125B
CN106802125B CN201710062745.4A CN201710062745A CN106802125B CN 106802125 B CN106802125 B CN 106802125B CN 201710062745 A CN201710062745 A CN 201710062745A CN 106802125 B CN106802125 B CN 106802125B
Authority
CN
China
Prior art keywords
electric thruster
ion current
plane
plume
current peak
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
CN201710062745.4A
Other languages
English (en)
Other versions
CN106802125A (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 Space Propulsion
Original Assignee
Shanghai Institute of Space Propulsion
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 Space Propulsion filed Critical Shanghai Institute of Space Propulsion
Priority to CN201710062745.4A priority Critical patent/CN106802125B/zh
Publication of CN106802125A publication Critical patent/CN106802125A/zh
Application granted granted Critical
Publication of CN106802125B publication Critical patent/CN106802125B/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
    • G01B7/00Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/30Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/30Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes
    • G01B7/31Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes
    • G01B7/312Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes for measuring eccentricity, i.e. lateral shift between two parallel axes

Abstract

本发明公开了一种航天器用电推力器推力矢量偏心的确定方法,所述航天器用电推力器推力矢量偏心的确定方法采用羽流测量装置、电推力器,使用羽流测量装置测量电推力器的离子电流分布,由离子电流峰的平面位置得到离子电流峰的空间位置,由离子电流峰的空间位置最终确定电推力器的推力矢量偏心。本发明相比多探针阵列确定电推力器推力偏心的传统方法,对测试设备要求简单,测量精度高、操作便利。

Description

航天器用电推力器推力矢量偏心的确定方法
技术领域
本发明涉及一种确定方法,特别是涉及一种航天器用电推力器推力矢量偏心的确定方法。
背景技术
卫星、深空探测器等各类航天器已广泛使用电推力器。但装配、加工尺寸偏差等非对称因素造成电推力器推力方向偏离其中心轴线,推力起始点不从推力器中心,形成推力矢量绕卫星质心的力矩,使卫星的推力效率降低,需及时启动飞轮或姿控化学推进系统产生反力矩,以防卫星姿态发生偏离。
目前电推力器推力矢量偏心采用安装在移动机构上的多探针阵列进行直接测量,多探针阵列面上密排大量探针,多探针阵列通过移动机构移动到被测羽流区域内某一位置,阵列面上的每个探针同时测量羽流区平面上各点的离子电流,得到离子电流峰的空间坐标和距离推力器中心的距离,计算推力矢量的偏角。
这种多探针阵列测量方式,存在如下问题:
(1)需要采用大量探针(数量从几十到几百个不等)有序排列进行测量,对探针供电、数据采集及机电接口要求大幅提升,系统复杂性明显增加;
(2)空间分辨率不高,只能测量羽流面内各个离散点的数据,特别在羽流中心区附近可能会遗漏掉离子电流峰的特征;
(3)多探针阵列各探针间的距离较近,各探针独立测量就会产生相互干扰,影响测量准确性,但间隔过大又会带来测量关键测量参数的遗漏,因此,需要在测量准确性和空间分辨率之间权衡利弊,确定合理的探针间隔;
(4)引入探针的数量越多,对羽流的畸变干扰越大,不利于测量的真实性;
(5)探针阵列板面积较大,正对羽流区,会产生较严重的返流,使测量失真;
(6)受探针本身物理尺寸的影响,探针视角会产生干涉,测量精度受限。
发明内容
本发明所要解决的技术问题是提供一种航天器用电推力器推力矢量偏心的确定方法,其系统复杂性低;空间分辨率高,从而测量精度高;探针较小,对羽流干扰忽略不计,保证测量不失真,有利于测量的真实性。
本发明是通过下述技术方案来解决上述技术问题的:一种航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述航天器用电推力器推力矢量偏心的确定方法采用羽流测量装置、电推力器,使用羽流测量装置测量电推力器的离子电流分布,由离子电流峰的平面位置得到离子电流峰的空间位置,由离子电流峰的空间位置最终确定电推力器的推力矢量偏心。
优选地,所述羽流测量装置包括法拉第探针、弧度运动机构、激光瞄准器、数据采集处理和控制机构,弧度运动机构包括可调摆杆、安装板、中心转轴、高精度步进电机,电推力器位于羽流测量装置正前方,通过旋转电推力器的安装方向,羽流测量装置上弧度运动机构的摆杆以两个不同半径沿电推力器的中心水平面、中心垂直面进行扫描测量,羽流测量装置上的法拉第探针获取离子电流分布,得到离子电流峰的平面位置,再由离子电流峰的平面位置得到离子电流峰空间位置,由离子电流峰空间位置最终确定电推力器的推力矢量偏心。
优选地,所述航天器用电推力器推力矢量偏心的确定方法包括以下步骤:
步骤一:将电推力器正立安装;
步骤二:弧度运动机构的摆杆半径设定为R1;
步骤三:确保电推力器喷口面的中心垂直轴和羽流测量装置的中心转轴共轴;
步骤四:确保羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心水平面;
步骤五:安装法拉第探针的弧度运动机构对电推力器的中心水平面扫描以获取其上离子电流分布;
步骤六:弧度运动机构的摆杆半径设定为R2,重复步骤五的操作;
步骤七:将电推力器侧立安装;
步骤八:弧度运动机构的摆杆半径设定为R1;
步骤九:确保电推力器喷口面的中心水平轴和羽流测量装置的中心转轴共轴;
步骤十:确保羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心垂直面;
步骤十一:安装法拉第探针的弧度运动机构对电推力器的中心垂直面扫描以获取其上离子电流分布;
步骤十二:弧度运动机构的摆杆半径设定为R2,重复步骤十一的操作;
步骤十三:由步骤五获取的电推力器中心水平面二维离子电流峰的平面位置和由步骤十一获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R1的球面上离子电流峰空间位置;
步骤十四:由步骤六获取的电推力器中心水平面二维离子电流峰的平面位置和步骤十二获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R2的球面上离子电流峰空间位置;
步骤十五:连接半径R1、R2球面上的离子电流峰的空间位置点,形成电推力器推力矢量方向;
步骤十六:对电推力器推力矢量偏角精确求解;
步骤十七:确定电推力器推力矢量偏心距。
优选地,所述步骤三、步骤九采用小锤确定羽流测量装置中心转轴相对于电推力器喷口面中心轴的位置。
优选地,所述步骤四、步骤十采用激光瞄准器对电推力器进行中心定位。
优选地,所述步骤五、步骤六以电推力器喷口面的中心垂直轴为圆心,弧度运动机构分别以摆杆半径R1、R2对电推力器的中心水平面从0°至180°进行扫描,获得相应的离子电流分布,得到离子电流峰在该面上的位置。
优选地,所述步骤十一、步骤十二以电推力器喷口面的中心水平轴为圆心,弧度运动机构分别以摆杆半径R1、R2对电推力器的中心垂直面从0°至180°进行扫描,获得相应的离子电流分布,得到离子电流峰在该面上的位置。
本发明的积极进步效果在于:本发明具有以下优点:
一、所采用的装置简单易用,记录形式直观,直接获得二维离子电流分布的形状;
二、所采用装置对羽流干扰较小,确保测试数据的真实性;
三、所采用装置对羽流返流较小,不影响电推力器正常工作;
四、采用单探针连续测量,避免了多探针测量时的相互干扰和测量数据的离散化;
五、空间分辨率、测量精度均较高;
六、所采用数学变换方法基于立体几何,简单、快速、实用。
附图说明
图1为离子速度分布示意图。
图2(a)、(b)为本发明的示意图。
图3为本发明的装置图。
图4为本发明安装法拉第探针的弧度运动机构对电推力器中心水平面进行扫描图。
图5为本发明安装法拉第探针的弧度运动机构对电推力器中心垂直面进行扫描图。
图6为本发明确定空间离子电流峰位置并取得推力矢量作用方向的示意图。
图7(a)、(b)、(c)为本发明中心垂直面上离子电流峰所在位置的张角示意图。
图8(a)、(b)、(c)为本发明中心水平面上离子电流峰所在位置的张角示意图。
图9(a)、(b)、(c)为本发明偏角精确求解图。
图10(a)、(b)、(c)为本发明偏心距精确求解图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1至图10所示,本发明航天器用电推力器推力矢量偏心的确定方法采用羽流测量装置1、电推力器2,使用羽流测量装置测量电推力器的离子电流分布,由离子电流峰的平面位置得到离子电流峰的空间位置,由离子电流峰的空间位置最终确定电推力器的推力矢量偏心。
其中,羽流测量装置1包括法拉第探针3、弧度运动机构4、激光瞄准器5、数据采集处理和控制机构6,弧度运动机构4包括可调摆杆7、安装板8、中心转轴9、高精度步进电机10,电推力器2位于羽流测量装置1正前方,通过旋转电推力器的安装方向,羽流测量装置上弧度运动机构的摆杆以两个不同半径沿电推力器的中心水平面、中心垂直面进行扫描测量,羽流测量装置上的法拉第探针获取离子电流分布,得到离子电流峰的平面位置,再由离子电流峰的平面位置得到离子电流峰空间位置,由离子电流峰空间位置最终确定电推力器的推力矢量偏心。
电推力器2的定位基准包括:电推力器中心水平面11、电推力器中心垂直面12、电推力器喷口面中心水平轴13、电推力器喷口面中心垂直轴14。
由悬挂在电推力器喷口中心的小锤15,使电推力器喷口面中心水平轴13(中心垂直轴14)和弧度运动机构的中心转轴9对准。
激光瞄准器5由螺纹结构16安装在法拉第探针3的壳体上。
电推力器是一种利用电离产生的离子在电场下加速产生推力的装置。电推力器的推力T为:
Figure BDA0001218835790000051
其中离子质量为mi、离子电流为Ii、离子速度为vi、θ为离子速度相对于推力器轴线的夹角、<vicosθ>为对离子电流分布取平均得到的离子轴向平均速度、e为电子电荷,可见,推力T只和离子电流Ii、离子轴向平均速度<vicosθ>有关,而<vi>可用Ii来反映,因此,可用离子电流及其轴向分量来表征推力矢量,如图1所示。
离子出射速度、密度、能量和方向反映了推力作用的大小、方向。电推力器的离子从圆环放电室内发出,每个离子出射形成单位推力矢量,每个离子的推力合成总推力,出射离子以辐射状离散地向外喷射,在电推力器外的空间形成离子电流分布。各离子能量不同,近似地,能量高的离子主要集中在推力器中心轴线附近,边缘能量较低,可以按最大离子电流对应空间点作为推力矢量作用点。选取距推力器喷口中心不同的两个半径,将其在两个半球面上的最大离子电流所对应的空间点连接起来,形成推力矢量,它和推力器中心轴的夹角即为推力矢量的偏角α,将它反向延长和推力器喷口面相交,其交点和推力器喷口面中心点之间的距离叫推力矢量偏心距δ。
电推力器分别处于正立、侧立两种状态,弧度运动机构分别以两个不同摆杆半径R1、R2对电推力器羽流进行扫描,得到其中心水平面、垂直面上的四组离子电流分布。由半径R1、在电推力器中心水平、垂直两个面上的最高离子电流所在的坐标可以确定离子电流峰对应的空间坐标点a1,由半径R2、在电推力器中心水平、垂直两个面上的最高离子电流所在的坐标可以确定离子电流峰对应的空间坐标点a2,a1a2连线即为推力矢量作用方向。a1a2连线和推力器中心轴的夹角即为推力矢量的偏角α;将a1a2连线反向延伸至推力器喷口面,形成交点a’,该点和推力器中心点的距离即为推力矢量的偏心距δ。
本发明航天器用电推力器推力矢量偏心的确定方法包括以下步骤:
如图2(a)所示,电推力器推力矢量偏心由电推力器推力矢量作用方向偏离电推力器中心轴线的角度α和在电推力器喷口面上推力矢量作用点偏离中心轴线的距离δ来反映。
步骤一:将电推力器正立安装;
步骤二:弧度运动机构的摆杆半径设定为R1,先将羽流测量装置的弧度运动机构的臂杆调节至R1
步骤三:由电推力器喷口中心悬挂的小锤尖锋对准弧度运动机构的中心转轴,确保电推力器喷口面中心垂直轴和弧度运动机构的中心转轴共轴;
步骤四:由安装在法拉第探针的激光瞄准器的光心对准电推力器的喷口中心,使羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心水平面;
如图3所示,采用激光瞄准器对电推力器进行中心定位,激光瞄准器和法拉第探针的外壳以螺纹拧紧,由螺纹精度确保法拉第探针的探测收集孔中心线和激光瞄准器的激光准心同轴,使激光瞄准器发出一束激光对准电推力器的中心。若不对中,适当调整电推力器的安装位置使其对中。
步骤五:定位完毕后,将小锤取下,将激光瞄准器从法拉第探针上拆卸下来,法拉第探针进入待命测试状态;安装法拉第探针的弧度运动机构对电推力器的中心水平面扫描以获取其上离子电流分布;
步骤六:弧度运动机构的摆杆半径设定为R2,安装法拉第探针的弧度运动机构对电推力器的中心水平面扫描以获取其上离子电流分布;
如图4所示,对电推力器的中心水平面进行扫描,其上安装的法拉第探针测得0°~180°的离子电流分布;再将弧度运动机构的臂杆调节至R2,得到另一组0°~180°的离子电流分布。
步骤七:如图5所示,将电推力器旋转90°,将电推力器侧立安装;
步骤八:先将羽流测量装置的弧度运动机构的臂杆调节至R1
步骤九:由电推力器喷口中心悬挂的小锤尖锋对准弧度运动机构的中心转轴,确保电推力器喷口面中心水平轴和弧度运动机构的中心转轴共轴;
步骤十:确保羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心垂直面,由安装在法拉第探针的激光瞄准器的光心对准电推力器的喷口中心,使羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心垂直面;
步骤十一:定位完毕后,将小锤取下,将激光瞄准器从法拉第探针上拆卸下来,法拉第探针进入待命测试状态;安装法拉第探针的弧度运动机构对电推力器的中心垂直面扫描以获取其上离子电流分布;
步骤十二:弧度运动机构的摆杆半径设定为R2,安装法拉第探针的弧度运动机构对电推力器的中心垂直面扫描以获取其上离子电流分布;
如图5所示,对电推力器的中心垂直面进行扫描,其上安装的法拉第探针测得0°~180°的离子电流分布;再将弧度运动机构的臂杆调节至R2,得到另一组0°~180°的离子电流分布。
步骤十三:由步骤五获取的电推力器中心水平面二维离子电流峰的平面位置和由步骤十一获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R1的球面上离子电流峰空间位置,由电推力器二维离子电流峰位置确定空间离子电流峰位置;
步骤十四:由步骤六获取的电推力器中心水平面二维离子电流峰的平面位置和步骤十二获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R2的球面上离子电流峰空间位置;
如图6所示,由电推力器中心垂直面、中心水平面上的二维离子电流峰所在的位置以及臂杆半径R1这三个维度参数,可以确定半径R1的球面上空间离子电流峰P1所对应的空间位置点a1。同理,可以确定半径R2的球面上空间离子电流峰P2的位置点a2。由弧度运动机构的摆杆半径R1所确定空间离子电流峰的位置点a1和摆杆半径R2所确定空间离子电流峰的位置点a2连接成线a1a2,该连线就是电推力器推力矢量作用方向,该连线和推力器中心轴的夹角即为电推力器推力矢量偏角α。该连线反向延长,与推力器喷口面形成交点,该交点和推力器中心轴线的距离就是电推力器推力矢量的偏心距δ。
步骤十五:由R1、R2球面上的离子电流峰的空间位置点,分别确定中心垂直面上离子电流峰所在位置的垂直张角;
由图7(a)、(b)、(c)可见,当半径R1时,垂直张角为θ1;当半径为R2,垂直张角为θ2;由图7(a)可见,∠a1O’f1=∠e1O’O1=θ1,由图7(b)可见,∠a2O’f2=∠e2O’O2=θ2
步骤十六:由R1、R2球面上的离子电流峰的空间位置点,分别确定中心水平面上离子电流峰所在位置的水平张角;
由图8(a)、(b)、(c)可见,当半径R1时,水平张角为φ1;当半径为R2,水平张角为φ2;由图8(a)可见,∠a1O’e1=∠f1O’O1=φ1,由图8(b)可见,∠a2O’e2=∠f2O’O2=φ2
步骤十七:电推力器推力矢量偏角精确求解。
如图9(a)所示,以电推力器喷口中心为原点O’点建立第一个笛卡尔坐标系XYZ,电推力器喷口面位于该坐标系的XO’Y平面上,以O’点为球心,半径为R1的球面上离子电流峰所对应的空间位置点为a1,点a1在XO’Z面上的投影为点b1,点b1在YO’Z面上的投影为点O1,将点a1、b1、O1连接构成一个直角三角形Δa1b1O1。以O1点为原点可建立第二个笛卡尔坐标系X1Y1Z。
如图9(b)所示,以O’点为球心,半径为R2的球面上离子电流峰所对应的空间位置点为a2。点a2在XO’Z面上的投影为点b2,点b2在在YO’Z面上的投影为点O2,将点a2、b2、O2连接构成一个直角三角形Δa2b2O2。以O2点为原点建立第三个笛卡尔坐标系X2Y2Z。
如图9(c)所示,a’点是a1a2的连线在推力器喷口面XO’Y上的交点,点a’在XO’Z面上的投影为点b’,因此,点a’、b’、O’连接构成一个直角三角形Δa’b’O’。
由示XYZ、X1Y1Z、X2Y2Z这三个坐标系的特点为:
(1)O’Z、O1Z、O2Z轴共轴;
(2)XO’Z、X1O1Z、X2O2Z面共面;
(3)XO’Y、X1O1Y1、X2O2Y2面相互平行。
电推力器推力矢量偏角α的确定过程如公式(1)至(16):
如见图9(a),用如下公式(1)至(6):
a1b1=R1·sinθ1……………………(1)
O'b1=R1·cosθ1……………………(2)
O1b1=O'b1·sinφ1=R1·cosθ1·sinφ1…………………(3)
O'O1=O'b1·cosφ1=R1·cosθ1·cosφ1…………………(4)
Figure BDA0001218835790000082
同理,如见图9(b),用如下公式(7)至(12):
a2b2=R2·sinθ2………………(7)
O'b2=R2·cosθ2…………………(8)
O2b2=O'b2·sinφ2=R2·cosθ2·sinφ2……………………(9)
O'O2=O'b2·cosφ2=R2·cosθ2·cosφ2…………………(10)
Figure BDA0001218835790000083
如见图9(c),从O2点引O2e2,使O1a1//O2e2,O1a1=O2e2,得到:a1e2//O1O2,a1e2=O1O2,O1O2⊥X2O2Y2面,有a1e2⊥X2O2Y2面,a2e2位于X2O2Y2面内,因此,Δa1a2e2为直角三角形。
O1O2=O'O2-O'O1………………(13)
将(4)、(10)式代入(13)式,得
O1O2=R2·cosθ2·cosφ2-R1·cosθ1·cosφ1………………(14)
在Δa2e2O2中,由余弦定理可得:
Figure BDA0001218835790000091
其中,O2e2=O1a1,而O1a1由(6)式给出,O2a2由(12)式给出。
推力矢量偏角即为:
Figure BDA0001218835790000092
其中,a1e2=O1O2,O1O2由(14)式给出。
步骤十七:确定电推力器推力矢量偏心距。
由图10(a)、(b)可见,直角Δa’b’O’、直角Δa1b1O1、直角Δa2b2O2相互平行。由图10(c)可见,从b’点引线段b’c2,使b'c2//O'O2,线段b’c2和b1O1的交点为c1,因此,得到O'b'=O1c1=O2c2,b'c1=O'O1,b'c2=O'O2
在Δb2b’c2中,b1c1//b2c2,由平行定理有如下公式(17):
Figure BDA0001218835790000093
在式中,O1b1、O2b2、O’O1、O’O2已分别由(3)、(9)、(4)、(10)式给出,因此,由(17)式可求出O’b’。
从a’点引线段a’d2,使a'd2//b'b2,a’d2和a1b1的交点为d1,因此有a'b'=d1b1=d2b2,b'b1=a'd1,b'b2=a'd2
在Δa2a’d2中,a1d1//a2d2,根据平行传递关系,得如下公式(18):
Figure BDA0001218835790000094
在式中,a1b1、a2b2、O’O1、O’O2已分别由(1)、(7)、(4)、(10)式给出,因此,由(18)式可求出a’b’。
在直角Δa’b’O’中,a’b’和O’b’已分别由(17)、(16)式给出,偏心距为如下公式(19):
Figure BDA0001218835790000095
本发明为一种航天器用电推力器推力矢量偏心的确定方法,它通过先测量而后经数学变换,简单直观地反映出电推力器推力矢量偏心:包括电推力器推力矢量作用方向偏离其中心轴线的角度和在电推力器喷口面推力矢量作用点偏离中心轴线的距离。这种方法相比传统方法具有操作便利、测试接口简单、测试数据准确、干扰小、可信度高的优点。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (5)

1.一种航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述航天器用电推力器推力矢量偏心的确定方法采用羽流测量装置、电推力器,使用羽流测量装置测量电推力器的离子电流分布,由离子电流峰的平面位置得到离子电流峰的空间位置,由离子电流峰的空间位置最终确定电推力器的推力矢量偏心;
所述羽流测量装置包括法拉第探针、弧度运动机构、激光瞄准器、数据采集处理和控制机构,弧度运动机构包括可调摆杆、安装板、中心转轴、高精度步进电机,电推力器位于羽流测量装置正前方,通过旋转电推力器的安装方向,羽流测量装置上弧度运动机构的摆杆以两个不同半径沿电推力器的中心水平面、中心垂直面进行扫描测量,羽流测量装置上的法拉第探针获取离子电流分布,得到离子电流峰的平面位置,再由离子电流峰的平面位置得到离子电流峰空间位置,由离子电流峰空间位置最终确定电推力器的推力矢量偏心;电推力器分别处于正立、侧立两种状态,弧度运动机构分别以两个不同摆杆半径R1、R2对电推力器羽流进行扫描,得到其中心水平面、垂直面上的四组离子电流分布;
所述航天器用电推力器推力矢量偏心的确定方法包括以下步骤:
步骤一:将电推力器正立安装;
步骤二:弧度运动机构的摆杆半径设定为R1;
步骤三:确保电推力器喷口面的中心垂直轴和羽流测量装置的中心转轴共轴;
步骤四:确保羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心水平面;
步骤五:安装法拉第探针的弧度运动机构对电推力器的中心水平面扫描以获取其上离子电流分布;
步骤六:弧度运动机构的摆杆半径设定为R2,重复步骤五的操作;
步骤七:将电推力器侧立安装;
步骤八:弧度运动机构的摆杆半径设定为R1;
步骤九:确保电推力器喷口面的中心水平轴和羽流测量装置的中心转轴共轴;
步骤十:确保羽流测量装置的弧度运动机构所扫描的平面即为电推力器的中心垂直面;
步骤十一:安装法拉第探针的弧度运动机构对电推力器的中心垂直面扫描以获取其上离子电流分布;
步骤十二:弧度运动机构的摆杆半径设定为R2,重复步骤十一的操作;
步骤十三:由步骤五获取的电推力器中心水平面二维离子电流峰的平面位置和由步骤十一获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R1的球面上离子电流峰空间位置;
步骤十四:由步骤六获取的电推力器中心水平面二维离子电流峰的平面位置和步骤十二获取的电推力器中心垂直面的二维离子电流峰平面位置,联合确定半径R2的球面上离子电流峰空间位置;
步骤十五:连接半径R1、R2球面上的离子电流峰的空间位置点,形成电推力器推力矢量方向;
步骤十六:对电推力器推力矢量偏角精确求解;
步骤十七:确定电推力器推力矢量偏心距。
2.如权利要求1所述的航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述步骤三、步骤九采用小锤确定羽流测量装置中心转轴相对于电推力器喷口面中心轴的位置。
3.如权利要求1所述的航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述步骤四、步骤十采用激光瞄准器对电推力器进行中心定位。
4.如权利要求1所述的航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述步骤五、步骤六以电推力器喷口面的中心垂直轴为圆心,弧度运动机构分别以摆杆半径R1、R2对电推力器的中心水平面从0°至180°进行扫描,获得相应的离子电流分布,得到离子电流峰在该面上的位置。
5.如权利要求1所述的航天器用电推力器推力矢量偏心的确定方法,其特征在于,所述步骤十一、步骤十二以电推力器喷口面的中心水平轴为圆心,弧度运动机构分别以摆杆半径R1、R2对电推力器的中心垂直面从0°至180°进行扫描,获得相应的离子电流分布,得到离子电流峰在该面上的位置。
CN201710062745.4A 2017-01-24 2017-01-24 航天器用电推力器推力矢量偏心的确定方法 Active CN106802125B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710062745.4A CN106802125B (zh) 2017-01-24 2017-01-24 航天器用电推力器推力矢量偏心的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710062745.4A CN106802125B (zh) 2017-01-24 2017-01-24 航天器用电推力器推力矢量偏心的确定方法

Publications (2)

Publication Number Publication Date
CN106802125A CN106802125A (zh) 2017-06-06
CN106802125B true CN106802125B (zh) 2020-02-11

Family

ID=58987244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710062745.4A Active CN106802125B (zh) 2017-01-24 2017-01-24 航天器用电推力器推力矢量偏心的确定方法

Country Status (1)

Country Link
CN (1) CN106802125B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109580261A (zh) * 2017-09-28 2019-04-05 南京淳泰控制设备有限公司 一种离子推力器矢量调节机构驱动测试装置
CN112326253B (zh) * 2020-10-28 2022-08-16 哈尔滨工业大学 一种推力矢量偏心诊断装置
CN113465494B (zh) * 2021-06-29 2022-04-08 哈尔滨工业大学 一种霍尔推力器推力矢量偏心计算方法
CN114397906B (zh) * 2021-12-06 2023-06-06 中国西安卫星测控中心 地球静止卫星电推进转移快速高精度计算方法
CN115682919B (zh) * 2022-10-28 2023-08-11 哈尔滨工业大学 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614606B (zh) * 2009-07-30 2011-02-16 中国科学院力学研究所 一种检测空间等离子体推力器推力矢量的测量装置和方法
CN103287590A (zh) * 2013-05-07 2013-09-11 中国科学院力学研究所 一种用于空间环境模拟真空设备的复合装置
CN104202894A (zh) * 2014-07-29 2014-12-10 北京航空航天大学 一种用于离子推力器测量的法拉第探针
CN104406761A (zh) * 2014-10-25 2015-03-11 哈尔滨工业大学 霍尔推力器低频振荡时间尺度内羽流发散角测量方法
CN105004529A (zh) * 2015-07-20 2015-10-28 兰州空间技术物理研究所 一种电推进航天器羽流参数检测装置
CN105067274A (zh) * 2015-07-20 2015-11-18 兰州空间技术物理研究所 一种电推进航天器羽流参数获取方法
CN105116435A (zh) * 2015-07-13 2015-12-02 兰州空间技术物理研究所 一种基于法拉第探针阵列的离子推力器束流测试方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614606B (zh) * 2009-07-30 2011-02-16 中国科学院力学研究所 一种检测空间等离子体推力器推力矢量的测量装置和方法
CN103287590A (zh) * 2013-05-07 2013-09-11 中国科学院力学研究所 一种用于空间环境模拟真空设备的复合装置
CN104202894A (zh) * 2014-07-29 2014-12-10 北京航空航天大学 一种用于离子推力器测量的法拉第探针
CN104406761A (zh) * 2014-10-25 2015-03-11 哈尔滨工业大学 霍尔推力器低频振荡时间尺度内羽流发散角测量方法
CN105116435A (zh) * 2015-07-13 2015-12-02 兰州空间技术物理研究所 一种基于法拉第探针阵列的离子推力器束流测试方法
CN105004529A (zh) * 2015-07-20 2015-10-28 兰州空间技术物理研究所 一种电推进航天器羽流参数检测装置
CN105067274A (zh) * 2015-07-20 2015-11-18 兰州空间技术物理研究所 一种电推进航天器羽流参数获取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于LabVIEW 的霍尔推力器羽流诊断系统;李雪松;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20160630(第06期);第19-40页 *
离子推力器推力矢量偏角测试;唐福俊;《真空与低温》;20091231;第15卷(第2期);第81-85页 *
霍尔推力器束流分布特性实验研究;张乾鹏,康小录,施晨毅,余水淋;《上海航天》;20121231;第29卷(第4期);第49-53页 *

Also Published As

Publication number Publication date
CN106802125A (zh) 2017-06-06

Similar Documents

Publication Publication Date Title
CN106802125B (zh) 航天器用电推力器推力矢量偏心的确定方法
CN113739700B (zh) 坐标测量仪器的正交轴系统的调校方法
CN107703499B (zh) 一种基于自制地基激光雷达对准误差的点云误差校正方法
CN107290734B (zh) 一种基于自制地基激光雷达垂直度误差的点云误差校正方法
GB2516528A (en) Automatic measurement of dimensional data with a laser tracker
CN111811395B (zh) 基于单目视觉的平面位姿动态测量方法
CN101992407A (zh) 设备的误差辨识方法和误差辨识程序
CN101169323A (zh) 在室内用旋转靶标测量电视经纬仪动态测角精度的方法
CN111366908B (zh) 一种激光雷达转台及其测量装置和测量方法
CN106425691B (zh) 基于激光干涉原理的精密主轴回转精度检测装置及方法
CN208968469U (zh) 工业机器人重复定位精度分析系统
CN104197831B (zh) 一种六关节工业机器人的精度标定装置
CN107478195A (zh) 一种基于光学的空间物体姿态测量装置及其测量方法
CN110940296A (zh) 一种高超声速飞行器舵偏角测量方法
US20100325863A1 (en) Surface sensing device
CN112697074B (zh) 动态待测物角度测量仪及测量方法
CN114440790A (zh) 同时检测薄壁回转体内外壁面形与厚度分布的方法与装置
CN109945824B (zh) 一种载机航向基准测量与传递的方法和系统
CN110631523B (zh) 一种精密二维转台轴间位置误差的测量装置及方法
CN108489396A (zh) 一种二维转顶精度检测方法
CN114562962A (zh) 一种基于激光跟踪仪的设备同轴度测量方法
CN112432622A (zh) 基于单测距传感器的曲面法向的测量装置及方法
CN107514990B (zh) 一种基于位移传感器的火炮弹丸卡膛姿态测试方法
CN218455443U (zh) 激光跟踪仪姿态测量校准装置
RU2785949C1 (ru) Устройство для измерения параметров паза на наружной сферической поверхности детали

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