CN106248300B - 基于成对推力器连续工作的卫星质心位置测量方法 - Google Patents
基于成对推力器连续工作的卫星质心位置测量方法 Download PDFInfo
- Publication number
- CN106248300B CN106248300B CN201610601089.6A CN201610601089A CN106248300B CN 106248300 B CN106248300 B CN 106248300B CN 201610601089 A CN201610601089 A CN 201610601089A CN 106248300 B CN106248300 B CN 106248300B
- Authority
- CN
- China
- Prior art keywords
- thruster
- satellite
- group
- centroid
- moment
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/12—Static balancing; Determining position of centre of gravity
- G01M1/122—Determining position of centre of gravity
- G01M1/125—Determining position of centre of gravity of aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供一种基于成对推力器连续工作的卫星质心位置测量方法,对推力器安装布局进行设计;选择成对推力器进行连续喷气工作;地面处理遥测下传的陀螺测量数据,对卫星的质心位置进行计算。本发明通过选定推力器进行固定时间长度连续工作,可避免推力器的输出特性变化带来的测量误差;在数据处理中通过对推力器产生的力矩取比值的方式,避免了推力大小偏差引起的测量误差,提高了质心测量精度。
Description
技术领域
本发明涉及卫星姿态与轨道控制技术,尤其是对于储箱平铺安装结构的地球静止轨道卫星在转移轨道段期间的姿轨控技术,利用本发明对卫星的质心位置进行测量,为储箱的均衡排放、混合比调节提供参考。
背景技术
受运载能力所限,地球静止轨道卫星需要星上的变轨发动机进行远地点变轨,将卫星由大椭圆轨道送入地球静止轨道。变轨过程中随着燃料消耗,由于储箱的不均衡排放、混合比偏差等因素的影响,卫星的质心位置会发生较大变化,可能会导致发动机推力偏心产生较大的干扰,增加姿控消耗的燃料。由于不均衡排放、混合比偏差等不确定因素,地面难以准确计算出卫星的质心位置。为减小变轨时的姿态扰动力矩,实现对卫星质心位置在轨调节,有必要对卫星的质心位置进行在轨测量。
当前大多数方案均基于推力器随机离散喷气的方式,推力器工作在小脉冲宽度条件下,利用递推算法对质心位置进行在轨计算,星上处理数据量大,计算复杂;另外,在小脉冲宽度条件下,推力器实际产生的推力大小变化较大,并且与标称值之间存在较大偏差,导致质心测量误差大。
发明内容
针对现有技术存在的缺点,本发明提供一种基于成对推力器连续工作的方式,利用陀螺测量数据实现对卫星质心位置测量。该发明通过选择相应的推力器在固定时间内连续工作,将陀螺测量数据遥测下传,结合推力器的安装位置及角度,通过对陀螺测量数据进行处理计算得到卫星的质心位置。
为了达到上述目的,本发明的技术方案是提供一种基于成对推力器连续工作的卫星质心位置测量方法:
布置至少一组推力器,每一组中的两个推力器产生相反控制力矩;
选取被布置的其中任意一组推力器,控制这一组的两个推力器在固定时间内同时连续喷气;
将这一组的两个推力器同时连续喷气期间的陀螺测量数据遥测下传进行处理,得到卫星三轴惯性角速度;
根据卫星三轴惯性角速度与卫星转动惯量,基于姿态动力学方程计算这一组中两个推力器在同时连续喷气期间各自产生的力矩及力矩的比值;
被布置的其中一组推力器是第一组推力器(2A、3A),所述第一组推力器(2A、3A)安装位置的y向坐标,相对于卫星y向平均质心位置对称;计算第一组推力器(2A、3A)同时连续喷气产生的俯仰力矩与偏航力矩,根据俯仰力矩与偏航力矩的比值以及第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标,计算卫星x向质心位置。
优选地,被布置的其中一组推力器是第二组推力器(4A、5A),所述第二组推力器(4A、5A)安装位置的x向坐标,相对于卫星x向平均质心位置对称;计算第二组推力器(4A、5A)同时连续喷气产生的滚动力矩与偏航力矩,根据滚动力矩与偏航力矩的比值以及第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标,计算卫星y向质心位置。
优选地,被布置的其中一组推力器是第三组推力器(6A、7A),所述第三组推力器(6A、7A)安装位置的z向坐标,相对于卫星z向平均质心位置对称;计算第三组推力器(6A、7A)同时连续喷气产生的滚动力矩与俯仰力矩,根据滚动力矩与俯仰力矩的比值以及第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标,计算卫星z向质心位置。
优选地,每一组中的各个推力器均倾斜安装,保证每个推力器产生的推力与卫星本体的轴不平行。
优选地,第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标,第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标,第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标,是在地面安装相应推力器时分别通过测量得到。
本发明提出了一种基于成对推力器连续工作的卫星质心位置测量方法。依据此方法可实现对卫星的质心位置进行测量。本发明中对推力器安装布局进行设计;选择推力器进行连续喷气工作;地面处理遥测下传的陀螺测量数据,对卫星的质心位置进行计算。
本发明采用的方法与现有技术相比,其优点在于:
本发明采用产生相反控制力矩的两个推力器同时连续喷气的方式,保证推力器产生的推力近似为恒值,从而使得喷气过程中卫星的角加速度近似为常值,避免了推力器的推力大小变化所引起的测量误差;此外,通过对力矩取比值的方式,可消除推力器稳态推力大小偏差带来的误差。
附图说明
图1为星上姿控推力器的安装布局示意图;
图2为质心测量过程流程图。
具体实施方式
本发明提供一种基于成对推力器连续工作的卫星质心位置测量方法,以下将结合附图和实施例对本发明作进一步说明。
1.推力器安装布局
为了保证能够利用该发明进行质心位置计算,对推力器的安装布局要求如下:
(1)推力器倾斜安装,保证每个推力器产生的推力与卫星本体轴不平行;
(2)滚动推力器的推力与星体的YOZ平面相互平行;俯仰推力器的推力与星体的XOZ平面相互平行;偏航推力器的推力与星体的XOY平面相互平行;
(3)产生相反控制力矩的两个推力器按照平均质心位置对称安装。
本发明中典型的推力器安装布局如图1所示,推力器2A、3A安装位置的y向坐标相对于卫星y向平均质心位置对称;推力器4A、5A安装位置的x向坐标相对于卫星x向平均质心位置对称;推力器6A、7A安装位置的z向坐标相对于卫星z向平均质心位置对称。
2.选择推力器及工作时间长度
在质心测量过程中,选择产生方向相反控制力矩的两个推力器在固定时间内连续喷气的方式,可保证测量过程中推力器产生的推力近似为常值,同时可保证卫星的角速度不会过大;此外,延长推力器的工作时间,可以获取足够多的测量信息,减小陀螺测量噪声对测量精度的影响。
3.对遥测下传的陀螺测量的角速度信息进行处理
将推力器喷气期间陀螺测量得到的三轴惯性角速度遥测下传,通过拟合的方式可估算出推力器喷气引起的三轴惯性角加速度;由姿态动力学方程可计算出推力器产生的作用力矩然后根据推力器的安装位置及角度(地面安装时可测量得到),由力矩公式可计算出卫星的质心位置。
按照图1中推力器的安装布局图,选择成对工作的推力器2A与3A、4A与5A、6A与7A分别“喷气”,实现对卫星X、Y、Z向质心位置的测量。如图2所示,本发明的方法包含以下步骤:
①t0时刻注数使系统进入到质心测量模式后,星上计算机停止发送控制指令,通过地面注数控制推力器2A、3A在t0时刻开始同时工作,喷气时间长度为10s(暂定,可进行修改),将喷气期间卫星三轴惯性角速度遥测下传。
②在t0+100s时刻,地面注数使系统进入惯性保持模式。
③当卫星三轴角速度均小于0.01°/s时,通过地面注数使系统再次进入到质心测量模式。
④t1时刻注数使系统进入到质心测量模式后,通过地面注数控制推力器4A、5A在t1时刻开始同时工作,喷气时间长度为10s,将喷气期间卫星三轴惯性角速度遥测下传。
⑤在t1+100s时刻,地面注数使系统进入惯性保持模式。
⑥t2时刻注数使系统进入到质心测量模式后,通过地面注数控制推力器6A、7A在t2时刻开始同时工作,喷气时间长度为10s,将喷气期间卫星三轴惯性角速度遥测下传。
⑦在t2+100s时刻,地面注数使系统退出质心测量模式。
⑧对遥测下传的陀螺测量数据进行处理,计算卫星质心位置。
对上述①、④、⑥中遥测下传的数据进行处理的算法如下:
假设每个控制周期内的惯性角加速度不变,将遥测下传的陀螺测量数据按照形式进行拟合,其中,将为卫星三轴惯性角速度,由陀螺组合测量遥测下传得到,t为遥测下传的时间,t=t0+k*T(k=1,2,…),t0为第一拍数据对应的时间,常选择t0=0,T为遥测下传周期;c为常数;为推力器产生的卫星三轴惯性角加速度。通过利用最小二乘法估算得到推力器工作所引起的惯性角加速度考虑到卫星的三轴惯性角速度均为小量,卫星的姿态动力学方程可简化为下式所示:
式中,I为卫星转动惯量矩阵;为卫星三轴惯性角加速度;为卫星三轴惯性角速度;[Tx Ty Tz]T分别为卫星受到的外部力矩。
结合卫星转动惯量I及陀螺测量数据处理得到的三轴惯性角加速度由式(1)可估算出推力器产生的作用力矩。
令推力器2A、3A同时工作时在星体上产生的力矩为Tx_23、Ty_23、Tz_23,根据推力器的安装布局,标称条件下推力器2A、3A在星体坐标系下产生的推力分别为:
F2a=F2a_0·[0 cosβ2 cosγ2]
F3a=F3a_0·[0 cosβ3 cosγ3]
其中,F2a_0,F3a_0分别为推力器2A、3A产生推力的大小;
β2,γ2分别为推力器2A推力的俯仰、偏航方向余弦角;
β3,γ3分别为推力器3A推力的俯仰、偏航方向余弦角。
推力器2A、3A在卫星布局坐标系下的安装位置分别为:
p2a=[p2a_x p2a_y p2a_z],p3a=[p3a_x p3a_y p3a_z]
由力矩公式可知,推力器2A、3A同时工作时产生的俯仰、偏航力矩如下式所示:
Ty_23=-F2a_0·cosγ2·(p2ax-px)-F3a_0·cosγ3·(p3ax-px)
Tz_23=F2a_0.cosβ2·(p2ax-px)+F3a_0.cosβ3·(p3ax-px)
推力器2A、3A产生的推力大小偏差较小时,即F2a_0/F3a_0≈1,则有:
由式(2)可计算得到卫星x向质心位置:
式中,k1=Ty_23/Tz_23,
p2a_x,p3a_x分别为推力器2A、3A在卫星布局坐标系下的x向位置坐标,由推力器安装时测量得到。
同理,通过利用推力器4A、5A同时工作,由陀螺测量数据处理得到的三轴惯性角加速度及式(1)可估算出推力器4A、5A同时工作产生的作用力矩Tx_45Ty_45Tz_45,可计算出卫星y向质心位置为:
其中,k2=Tx_45/Tz_45,;
α4,γ4分别为推力器4A推力矢量的滚动、偏航方向余弦角;
α5,γ5分别为推力器5A推力矢量的滚动、偏航方向余弦角;
p4a_y,p5a_y分别为推力器4A、5A在卫星布局坐标系下位置的y向坐标,由安装时测量得到。
通过利用推力器6A、7A同时工作,由陀螺测量数据处理得到的三轴惯性角加速度及式(1)可估算出推力器6A、7A同时工作产生的作用力矩Tx_67Ty_67Tz_67,可计算得到卫星z向质心位置:
式中,α6,β6分别为推力器6A推力矢量的滚动、俯仰余弦角,由安装时测量得到;
α7,β7分别为推力器7A推力矢量的滚动、俯仰余弦角,由安装时测量得到;
p6a_z,p7a_z分别为推力器6A、7A在卫星布局坐标系下的z向位置坐标,由安装时测量得到;
k3=Tx_67/Ty_67,其中[Tx_67 Ty_67 Tz_67]T为推力器6A、7A同时工作产生的作用力矩,可由遥测下传的陀螺测量数据处理得到。
由px,py,pz表达式可知,卫星的质心位置与推力器的推力大小无关,避免推力大小变化引起的测量误差。
综上所述,与以往测量方式相比,本发明的方法通过选定推力器进行固定时间长度连续工作,可避免推力器的输出特性变化带来的测量误差;此外,在数据处理中通过对推力器产生的力矩取比值的方式,避免了推力大小偏差引起的测量误差,提高了质心测量精度。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (3)
1.一种基于成对推力器连续工作的卫星质心位置测量方法,其特征在于,
布置至少一组推力器,每一组中的两个推力器产生相反控制力矩;
选取被布置的其中任意一组推力器,控制这一组的两个推力器在固定时间内同时连续喷气;
将这一组的两个推力器同时连续喷气期间的陀螺测量数据遥测下传进行处理,得到卫星三轴惯性角速度;
根据卫星三轴惯性角速度与卫星转动惯量,基于姿态动力学方程计算这一组中两个推力器在同时连续喷气期间各自产生的力矩及力矩的比值;
被布置的其中一组推力器是第一组推力器(2A、3A),所述第一组推力器(2A、3A)安装位置的y向坐标,相对于卫星y向平均质心位置对称;计算第一组推力器(2A、3A)同时连续喷气产生的俯仰力矩与偏航力矩,根据俯仰力矩与偏航力矩的比值以及第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标,计算卫星x向质心位置;
被布置的其中一组推力器是第二组推力器(4A、5A),所述第二组推力器(4A、5A)安装位置的x向坐标,相对于卫星x向平均质心位置对称;计算第二组推力器(4A、5A)同时连续喷气产生的滚动力矩与偏航力矩,根据滚动力矩与偏航力矩的比值以及第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标,计算卫星y向质心位置;
被布置的其中一组推力器是第三组推力器(6A、7A),所述第三组推力器(6A、7A)安装位置的z向坐标,相对于卫星z向平均质心位置对称;计算第三组推力器(6A、7A)同时连续喷气产生的滚动力矩与俯仰力矩,根据滚动力矩与俯仰力矩的比值以及第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标,计算卫星z向质心位置。
2.如权利要求1所述的卫星质心位置测量方法,其特征在于,
每一组中的各个推力器均倾斜安装,保证每个推力器产生的推力与卫星本体的轴不平行。
3.如权利要求1所述的卫星质心位置测量方法,其特征在于,
第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标,第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标,第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标,是在地面安装相应推力器时分别通过测量得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610601089.6A CN106248300B (zh) | 2016-07-27 | 2016-07-27 | 基于成对推力器连续工作的卫星质心位置测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610601089.6A CN106248300B (zh) | 2016-07-27 | 2016-07-27 | 基于成对推力器连续工作的卫星质心位置测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106248300A CN106248300A (zh) | 2016-12-21 |
CN106248300B true CN106248300B (zh) | 2018-10-26 |
Family
ID=57603347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610601089.6A Active CN106248300B (zh) | 2016-07-27 | 2016-07-27 | 基于成对推力器连续工作的卫星质心位置测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106248300B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103152B (zh) * | 2017-05-16 | 2021-02-02 | 上海航天控制技术研究所 | 一种有利于抑制液体晃动的推力器安装布局方法 |
CN107933967A (zh) * | 2017-11-23 | 2018-04-20 | 北京控制工程研究所 | 一种卫星转动惯量的在轨辨识方法 |
CN111324137B (zh) * | 2020-03-05 | 2022-08-05 | 中国工程物理研究院总体工程研究所 | 一种基于常值连续推力的区域悬停轨道控制方法 |
CN112461324B (zh) * | 2020-10-30 | 2022-04-12 | 中国空间技术研究院 | 一种基于发动机干扰力矩确定并联贮箱推进剂剩余量方法 |
CN112577669B (zh) * | 2020-12-15 | 2022-08-12 | 上海卫星工程研究所 | 分体式卫星在轨两舱质心辨识方法及系统 |
CN115752897B (zh) * | 2022-11-25 | 2023-08-08 | 哈尔滨工业大学 | 姿态控制发动机的动态质心测量方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5459669A (en) * | 1994-02-14 | 1995-10-17 | Space Systems/Loral, Inc. | Control system and method for spacecraft attitude control |
CN101499220B (zh) * | 2009-01-24 | 2010-09-29 | 哈尔滨工业大学 | 一种模拟航天器上大型推力器的方法及其装置 |
CN103235598B (zh) * | 2013-05-14 | 2016-08-10 | 北京理工大学 | 一种调整推力器方向指向组合体航天器质心的方法 |
JP2015168315A (ja) * | 2014-03-06 | 2015-09-28 | 三菱電機株式会社 | 誘導装置、誘導装置搭載宇宙機 |
CN105620792B (zh) * | 2016-02-05 | 2017-12-05 | 上海微小卫星工程中心 | 一种采用斜装推力器进行卫星姿态和轨道控制的方法 |
-
2016
- 2016-07-27 CN CN201610601089.6A patent/CN106248300B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106248300A (zh) | 2016-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106248300B (zh) | 基于成对推力器连续工作的卫星质心位置测量方法 | |
US9643740B2 (en) | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter | |
US3866025A (en) | Spacecraft attitude control system | |
US8352101B2 (en) | Algorithm for simultaneous attitude maneuver and momentum dumping | |
US11325726B2 (en) | Method and apparatus for spacecraft gyroscope scale factor calibration | |
JP3027734B2 (ja) | 衛星の位置を機内で自主的に求める方法と装置 | |
US20080315039A1 (en) | System and methods for space vehicle torque balancing | |
US8321076B2 (en) | On-line inertia estimation for use in controlling an aerospace vehicle | |
CN103235509B (zh) | 一种基于动量轮的转动部件干扰补偿方法 | |
CN104527994A (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN108181916B (zh) | 小卫星相对姿态的控制方法及装置 | |
How et al. | Analysis of the impact of sensor noise on formation flying control | |
CN103955223A (zh) | 一种深空探测软着陆过程的姿轨耦合控制方法 | |
CN104960674A (zh) | 一种运动目标的指向跟踪控制方法 | |
US11787569B2 (en) | System and method for optimizing a low-thrust trajectory of a spacecraft trajectory | |
CN103019247A (zh) | 一种火星探测器无陀螺自主空间姿态机动控制方法 | |
RU2318188C1 (ru) | Способ автономной навигации и ориентации космических аппаратов | |
CN110803304A (zh) | 一种卫星姿态控制系统 | |
CN111638643A (zh) | 一种位移模式无拖曳控制动力学协调条件确定方法 | |
CN110955255B (zh) | 基于cmg的高精度轨控姿态维持方法、系统及介质 | |
CN106494643B (zh) | 一种基于地心矢量信息的卫星姿态异常恢复控制方法 | |
CN104932520B (zh) | 一种存在测量误差下的航天器雨滴形状悬停构型控制方法 | |
CN107132850B (zh) | 基于角速度跟踪的变轨姿态保持控制方法 | |
Fan et al. | Attitude optimization control method of agile optical small satellite for nonparallel ground track imaging | |
CN114802818A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |