法律状态公告日
法律状态信息
法律状态
2015-09-23
授权
授权
2013-12-18
实质审查的生效 IPC(主分类):G05D1/00 申请日:20130725
实质审查的生效
2013-11-27
公开
公开
技术领域
本发明涉及一种预估轨控综合效率和推进剂消耗量的方法,适用于中高 轨道卫星利用推力器进行轨道机动领域。
背景技术
中高轨道卫星在完成定点捕获、相位捕获及位置保持时,均利用推力器 进行轨道机动。导航卫星由于是组网星座卫星,其极限环内漂移天数要求多, 多颗星处于多星共轨状态,因此需要保证倾斜轨道各卫星之间的相位关系, 同时保证卫星轨道调整定点后的剩余推进剂的充沛。
现有导航卫星共安装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%。证明本次轨道机动控制精度高,轨控综合效率和推进剂余量估计 准确,降低了控制误差。记录轨控综合效率和推进剂余量用于以后的轨 道机动参数制定。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技 术。
机译: 质量预估算设备,质量预估算方法和质量预估算程序
机译: 预估距离计算器,预估距离计算方法,预估距离计算程序和自动计划器
机译: 推进剂在轨航天器补给的推进剂转移系统和方法