CN103412563B - 一种预估轨控综合效率和推进剂消耗量的方法 - Google Patents
一种预估轨控综合效率和推进剂消耗量的方法 Download PDFInfo
- Publication number
- CN103412563B CN103412563B CN201310316818.XA CN201310316818A CN103412563B CN 103412563 B CN103412563 B CN 103412563B CN 201310316818 A CN201310316818 A CN 201310316818A CN 103412563 B CN103412563 B CN 103412563B
- Authority
- CN
- China
- Prior art keywords
- thruster
- delta
- rail control
- control
- satellite
- 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
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种预估轨控综合效率和推进剂消耗量的方法,该方法根据之前在轨标定出的轨控推力器工作期间对卫星的扰动力矩,计算出其它推力器对扰动力矩进行补偿控制时对轨道位置产生的不期望影响,根据对轨道的不期望影响,修正轨控推力器工作时间,并计算出这颗卫星轨控推力器的轨控综合效率,根据修正过的所有推力器工作时间,估算本次轨道机动的推进剂消耗量及推进剂余量。本发明使卫星在沿轨道坐标系X方向进行轨道机动时的控制参数确定和推进剂消耗量的估算可以更加准确,有效的提高了中高轨道卫星轨道机动的精度,延长卫星在轨使用寿命。
Description
技术领域
本发明涉及一种预估轨控综合效率和推进剂消耗量的方法,适用于中高轨道卫星利用推力器进行轨道机动领域。
背景技术
中高轨道卫星在完成定点捕获、相位捕获及位置保持时,均利用推力器进行轨道机动。导航卫星由于是组网星座卫星,其极限环内漂移天数要求多,多颗星处于多星共轨状态,因此需要保证倾斜轨道各卫星之间的相位关系,同时保证卫星轨道调整定点后的剩余推进剂的充沛。
现有导航卫星共安装12个10N推力器,分成A、B两个分支,安装布局如图2所示,推力方向和寿命初期的力矩特性见表1:
表1某卫星推力器推力方向和寿命初期三轴力矩分量的工程参数
每个分支中推力器2和3是一对,分别产生-Z和+Z方向控制力矩,两者成对同时工作可产生+X的轨控推力;推力器4和5分别产生+Y和-Y力矩,成对工作可产生-X推力;推力器6和7分别产生+X和-X力矩,成对工作可产生+Y推力。在进行轨道机动的过程中,使用一个分支某2个推力器进行轨道控制。在实施轨道机动之前,需根据轨道机动目标预估推进剂消耗量并制定轨道机动控制参数,其中包括确定使用哪2个10N推力器作为轨控推力器及轨控推力器点火时间长度。
轨控综合效率是轨道机动控制参数制定的必要输入参数,需通过轨控综合效率计算出轨控推力器点火时间长度和推进剂消耗量。根据10N推力器的布局设计,在轨控推力器产生轨控推力的同时也会产生对卫星的扰动力矩和对卫星轨道的不期望影响,为实现力矩平衡和精确的轨道控制,需用其它推力器进行主动控制和补偿控制。轨控综合效率将这些因素综合考虑,制定出轨控推力器点火时间长度。轨控综合效率的准确性直接影响轨控效果和推进剂余量的计算,进而影响卫星寿命的估计。
卫星轨控综合效率和推进剂消耗量现有估算方法都是根据经验估算出来的,每颗卫星的质心位置、10N推力器安装、推力大小均存在差异,经验值不能保证准确体现每颗星的具体情况和准确性,所以现在需要研制一种沿轨道坐标系X轴(切向)轨道机动时高精度估算轨控综合效率和推进剂消耗量的方法。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种预估轨控综合效率和推进剂消耗量的方法,使卫星在沿轨道坐标系X方向进行轨道机动时的控制参数确定、轨控综合效率和推进剂消耗量的估算更加准确,有效地提高了中高轨道卫星轨道机动的精度,延长了卫星在轨使用寿命。
本发明的技术解决方案是:
一种预估轨控综合效率和推进剂消耗量的方法,其步骤如下:
(1)根据轨道机动选择的轨控推力器和控制速度增量ΔV,计算得到轨控推力器点火时间长度Δtc。
其中Md为干星质量;Mp为轨道机动前推进剂质量;ΔV为控制速度增量;Ke为推力器控制效率;Isp为推力器比冲;W0为氧化剂流率,Wf为燃烧剂流率;
(2)根据轨控推力器点火时间长度Δtc、轨控推力器工作时对卫星三轴的干扰力矩Txi、Tyi、Tzi和各推力器产生的力矩Mij,其中i为轨控推力器编号,取23A、45A、67A、23B、45B、67B,j取x、y、或z轴,计算出轨道机动过程中用于补偿控制的喷气时间ti:
使用A或B分支:
若Txi为正,则:
若Txi为负,则:
若Tyi为正,则:
若Tyi为负,则:
若Tzi为正,则:
若Tzi为负,则:
(3)根据步骤(2)中得到的各推力器用于补偿控制的喷气时间ti,计算其对轨道产生冲量Pi:
推力器2A、3A或2B、3B在+X方向的冲量为:
推力器4A、5A或4B、5B在+X方向的冲量为:
推力器6A、7A或6B、7B在+X方向的冲量为:
其中,偏航姿态偏置量为ψb,α2/3、α4/5、α6/7分别为推力器2/3、4/5、6/7的推力方向与卫星本体X轴的夹角;
(4)根据冲量Pi修正轨控推力器工作时间为预估出推力器的总工作时间t总,其中补偿控制引起的对沿轨道坐标系+X轴切向轨控量的补偿量Δtcom计算公式为:
补偿轨控推力器工作时间:
计算所有推力器总工作时间:
或
(5)计算轨控综合效率为
(6)依据推力器总工作时间t总计算推进剂消耗量ΔM*和推进剂余量
ΔM*=2(W0+Wf)×t总
(7)记录轨控综合效率和推进剂余量
本发明与现有技术相比的有益效果是:
(1)本发明提高了中高轨道卫星实施沿轨道坐标系X轴轨道机动的轨道控制精度和剩余推进剂估计的准确度,采用本发明提出的预估轨控综合效率和推进剂消耗量的方法,可给出适合每颗卫星具体在轨情况的轨控综合效率,提高了卫星轨道机动控制参数制定的准确性,从而延长了卫星轨道机动周期和在轨工作寿命。
(2)本发明设计过程简单,易于实现,是十分经济可行的估算轨控综合效率和推进剂消耗量方法,可推广到任何需要进行轨道机动的卫星上。
附图说明
图1为本发明流程图;
图2为本发明推力器安装示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
如图1所示,一种预估轨控综合效率和推进剂消耗量的方法,其具体步骤如下:
(1)根据此次轨道机动选择的轨控推力器和控制速度增量ΔV,计算得到轨控推力器点火时间长度Δtc。
其中ΔM为中间量,Md为干星(为不含任何燃料的卫星)质量;Mp为轨道机动前推进剂质量;ΔV为控制速度增量;Ke为推力器控制效率;Isp为推力器比冲;W0为氧化剂流率,Wf为燃烧剂流率;
(2)根据轨控推力器点火时间长度Δtc、轨控推力器工作时对卫星三轴的干扰力矩Txi、Tyi、Tzi和各推力器产生的力矩Mij,其中i为轨控推力器编号,取23A、45A、67A、23B、45B、67B,j取x、y、或z轴,计算出轨道机动过程中用于补偿控制的喷气时间ti:
使用A或B分支:
若Txi为正,则:
若Txi为负,则:
若Tyi为正,则:
若Tyi为负,则:
若Tzi为正,则:
若Tzi为负,则:
(3)根据步骤(2)中得到的各推力器用于补偿控制的喷气时间ti,计算其对轨道产生的影响,以冲量形式Pi描述:
设偏航姿态偏置量为ψb,根据推力器安装参数,选取α2/3、α4/5、α6/7值,其中α2/3、α4/5、α6/7分别为推力器2/3、4/5的推力方向与卫星本体X轴的夹角以及推力器6/7与卫星本体Y轴的夹角。
推力器2A、3A或2B、3B在+X方向的冲量为:
推力器4A、5A或4B、5B在+X方向的冲量为:
推力器6A、7A或6B、7B在+X方向的冲量为:
(4)根据冲量形式Pi修正轨控推力器工作时间为预估出推力器的总工作时间t总,其中补偿控制引起的对沿轨道坐标系+X轴切向轨控量的补偿量Δtcom计算公式为:
补偿轨控推力器工作时间:
计算所有推力器总工作时间:
或
(5)计算轨控综合效率为
(6)依据推力器总工作时间t总计算推进剂消耗量ΔM*和推进剂余量
ΔM*=2(W0+Wf)×t总
(7)记录轨控综合效率和推进剂余量
下面以一个具体实施例来进一步说明本发明的方法。
以一次定点捕获前制定其轨道机动控制参数操作为例,说明采用本发明实现计算轨控综合效率和推进剂消耗量的实施方法。
(1)根据定点捕获控制策略,计划进行卫星减速控制将轨道半长轴a降低-14.52km,控制量ΔV=-0.529m/s。由某高轨卫星工程参数,Md=1354.8kg,Mp=363.8kg轨道控制前卫星质量Md+MP=1714.8kg。因此,确定选用4A、5B推力器,推力方向为-X。并确定参数值:Ke=0.766;Isp=2696N·s/kg;W0=2.29×10-3kg/s,Wf=1.41×10-3kg/s。
计算得到理论轨控推力器点火时间长度Δtc。
(2)4A、5B推力器工作时对卫星三轴的干扰力矩Tx45A=-0.08Nm、Ty45A=0.2Nm、Tz45A=-0.04Nm,M6Ax=6.00359Nm,M5Ay=-14.56311Nm,M3Az=-11.80633Nm。
计算出轨道机动过程中用于补偿控制的喷气时间:
其它推力器喷气时间为0。
(3)根据步骤(2)计算得到的t5A、t6A和t3A,偏航姿态偏置量为ψb=0°,,α4/5=40°,α2/3=15°,α6/7=-30°(根据表1中推力方向)情况确定。
Px45=-10cosψbcosα4/5×t5A=-6.24Ns
Px67=-10sinψbcosα6/7×t6A=0Ns
Px23=10cos(ψb+α2/3)×t3A=1.94Ns
(4)修正轨控推力器工作时间为
补偿轨控推力器工作时间:
计算所有推力器总工作时间:
(5)计算轨控综合效率为
(6)依据推力器总工作时间t总计算推进剂消耗量ΔM*和推进剂余量
ΔM*=2(W0+Wf)×t总=0.454kg
(7)通过轨控综合效率的修正,轨控推力器的点火时间变化了-0.5%。经过测轨,此次轨道机动的效果为半长轴a降至42165.222km,比控制目标值42165.208km略低0.014km,计划控制量为-14.518km。控制误差0.0096%。证明本次轨道机动控制精度高,轨控综合效率和推进剂余量估计准确,降低了控制误差。记录轨控综合效率和推进剂余量用于以后的轨道机动参数制定。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (1)
1.一种预估轨控综合效率和推进剂消耗量的方法,其特征在于步骤如下:
(1)根据轨道机动选择的轨控推力器和控制速度增量ΔV,计算得到轨控推力器点火时间长度Δtc;
其中Md为干星质量;Mp为轨道机动前推进剂质量;ΔV为控制速度增量;Ke为推力器控制效率;Isp为推力器比冲;W0为氧化剂流率,Wf为燃烧剂流率;
(2)根据轨控推力器点火时间长度Δtc、轨控推力器工作时对卫星三轴的干扰力矩Txi、Tyi、Tzi和各推力器产生的力矩Mij,其中i为轨控推力器编号,取2A、3A、4A、5A、6A、7A、2B、3B、4B、5B、6B、7B,j取x、y、或z轴,计算出轨道机动过程中用于补偿控制的喷气时间ti:
使用A或B分支:
若Txi为正,则:
若Txi为负,则:
若Tyi为正,则:
若Tyi为负,则:
若Ti为正,则:
若Tzi为负,则:
(3)根据步骤(2)中得到的各推力器用于补偿控制的喷气时间ti,计算其对轨道产生冲量Pi:
推力器2A、3A或2B、3B在+X方向的冲量为:
推力器4A、5A或4B、5B在+X方向的冲量为:
推力器6A、7A或6B、7B在+X方向的冲量为:
其中,偏航姿态偏置量为ψb,α2/3、α4/5、α6/7分别为推力器2/3、4/5、6/7的推力方向与卫星本体X轴的夹角;
(4)根据冲量Pi修正轨控推力器工作时间为预估出推力器的总工作时间t总,其中补偿控制引起的对沿轨道坐标系+X轴切向轨控量的补偿量Δtcom计算公式为:
补偿轨控推力器工作时间:
计算所有推力器总工作时间:
或
(5)计算轨控综合效率为
(6)依据推力器总工作时间t总计算推进剂消耗量ΔM*和推进剂余量
ΔM*=2(W0+Wf)×t总
(7)记录轨控综合效率和推进剂余量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310316818.XA CN103412563B (zh) | 2013-07-25 | 2013-07-25 | 一种预估轨控综合效率和推进剂消耗量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310316818.XA CN103412563B (zh) | 2013-07-25 | 2013-07-25 | 一种预估轨控综合效率和推进剂消耗量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103412563A CN103412563A (zh) | 2013-11-27 |
CN103412563B true CN103412563B (zh) | 2015-09-23 |
Family
ID=49605588
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310316818.XA Active CN103412563B (zh) | 2013-07-25 | 2013-07-25 | 一种预估轨控综合效率和推进剂消耗量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103412563B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104590588B (zh) * | 2014-12-04 | 2016-06-08 | 哈尔滨工业大学 | 一种基于隔离余量方法与脉宽融合策略的挠性卫星姿态轨道耦合控制方法 |
CN105197257B (zh) * | 2015-08-31 | 2017-03-15 | 中国空间技术研究院 | 一种分舱优化设计的桁架式geo卫星推力器布局方法 |
CN106114909B (zh) * | 2016-06-23 | 2018-08-31 | 中国空间技术研究院 | 一种卫星变轨推进剂消耗量计算方法 |
CN106114911B (zh) * | 2016-06-23 | 2018-08-31 | 中国空间技术研究院 | 一种卫星推进系统在轨飞行混合比调整方法 |
CN107063567B (zh) * | 2016-11-09 | 2019-06-28 | 上海卫星工程研究所 | 高轨并联平铺贮箱卫星转移段质量特性计算方法 |
CN107352050B (zh) * | 2017-06-07 | 2019-06-18 | 北京空间飞行器总体设计部 | 大型纵轴对地微波天线载荷卫星的推力器配置布局方法 |
CN108516107B (zh) * | 2018-02-27 | 2020-11-20 | 北京控制工程研究所 | 一种在线发动机推力和比冲估计方法及系统 |
CN109515758B (zh) * | 2018-11-05 | 2020-09-18 | 中国空间技术研究院 | 一种利用残余推进剂和氦气的geo卫星离轨方法 |
CN109649689B (zh) * | 2018-12-07 | 2021-10-01 | 北京空间飞行器总体设计部 | 一种有限推力变轨重力损耗计算方法、推力计算装置 |
CN110239744B (zh) * | 2019-06-28 | 2020-12-22 | 北京理工大学 | 一种小天体着陆定推力轨迹跟踪控制方法 |
CN111688954B (zh) * | 2020-07-02 | 2021-10-22 | 北京空间技术研制试验中心 | 一种航天器贮箱系统的排放不平衡系数在轨估计方法 |
CN112078832B (zh) * | 2020-08-04 | 2021-12-07 | 上海航天控制技术研究所 | 一种在轨剩余燃料的确定方法 |
CN112149224B (zh) * | 2020-09-12 | 2022-05-31 | 中国运载火箭技术研究院 | 飞行器推进剂消耗量计算方法、设备、飞行器及存储介质 |
CN112607063B (zh) * | 2020-12-29 | 2021-10-08 | 中国人民解放军32039部队 | 一种卫星位置保持的控制方法和装置 |
CN115636109B (zh) * | 2022-11-17 | 2023-04-18 | 北京九天微星科技发展有限公司 | 一种使用双推力系统的轨控方法、装置及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6089507A (en) * | 1996-12-05 | 2000-07-18 | Parvez; Shabbir Ahmed | Autonomous orbit control with position and velocity feedback using modern control theory |
CN102424116A (zh) * | 2011-12-08 | 2012-04-25 | 中国空间技术研究院 | 一种静止轨道卫星变轨策略优化方法 |
CN103192997A (zh) * | 2013-03-20 | 2013-07-10 | 上海空间推进研究所 | 一种航天器推进剂混合比主动调节方法 |
-
2013
- 2013-07-25 CN CN201310316818.XA patent/CN103412563B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6089507A (en) * | 1996-12-05 | 2000-07-18 | Parvez; Shabbir Ahmed | Autonomous orbit control with position and velocity feedback using modern control theory |
CN102424116A (zh) * | 2011-12-08 | 2012-04-25 | 中国空间技术研究院 | 一种静止轨道卫星变轨策略优化方法 |
CN103192997A (zh) * | 2013-03-20 | 2013-07-10 | 上海空间推进研究所 | 一种航天器推进剂混合比主动调节方法 |
Non-Patent Citations (2)
Title |
---|
一种用于导弹姿态控制系统的直接力非线性控制法;奚勇,张巍巍;《上海航天》;20070806;第24卷(第3期);第42-45页 * |
近圆低轨道控制量简易计算方法;李英良;《载人航天》;20130331;第19卷(第2期);第52-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103412563A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103412563B (zh) | 一种预估轨控综合效率和推进剂消耗量的方法 | |
CN103072702B (zh) | 卫星轨道和姿态控制方法 | |
CN104015938B (zh) | 一种电推进静止轨道卫星的位置保持方法 | |
CN104142686A (zh) | 一种卫星自主编队飞行控制方法 | |
CN105629732B (zh) | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 | |
CN102591349B (zh) | 高轨道卫星大初始角速率情况的无陀螺太阳捕获控制方法 | |
CN103136444B (zh) | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 | |
CN104176275B (zh) | 一种使用动量轮与磁力矩器联合的速率阻尼方法 | |
CN101794154A (zh) | 一种编队卫星相对轨道与姿态解耦控制方法 | |
CN104058104B (zh) | 无加速度计情况下一种基于关调制的高精度轨控方法 | |
CN103424116A (zh) | 一种适应轨道机动的地球同步卫星精密定轨方法 | |
CN103587681A (zh) | 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法 | |
CN103076809B (zh) | 一种利用干扰积累角动量自平衡的轨道控制方法 | |
CN102004491B (zh) | 一种卫星初入轨段的初始太阳捕获方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
CN104015191B (zh) | 基于基座卫星角速度的空间机械臂工具坐标下的运动补偿方法 | |
CN107992062A (zh) | 一种基于混合执行机构的空间高动态目标高精度姿态跟踪控制方法 | |
CN104085539A (zh) | 成像定标的姿态控制方法 | |
CN108516107B (zh) | 一种在线发动机推力和比冲估计方法及系统 | |
CN103171775A (zh) | 一种基于轨控干扰的角动量耦合卸载方法 | |
CN109484675A (zh) | 一种利用空间矢量匹配的航天器轨道入轨控制方法 | |
CN103116361B (zh) | 一种卫星动量轮控制下的变轨间隔确定方法 | |
CN104035333B (zh) | 一种深空探测器动力下降初始关键参数优化方法 | |
CN104848857B (zh) | 弹道导弹惯性测量系统精度指标自动分配方法 | |
CN104504188B (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 |