CN106681138B - 一种编队卫星系统燃料消耗均衡在轨实现方法 - Google Patents
一种编队卫星系统燃料消耗均衡在轨实现方法 Download PDFInfo
- Publication number
- CN106681138B CN106681138B CN201611097493.0A CN201611097493A CN106681138B CN 106681138 B CN106681138 B CN 106681138B CN 201611097493 A CN201611097493 A CN 201611097493A CN 106681138 B CN106681138 B CN 106681138B
- Authority
- CN
- China
- Prior art keywords
- star
- formation control
- orbit
- satellite
- formation
- 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
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明一种编队卫星系统燃料消耗均衡在轨实现方法,步骤如下:1)第k次编队控制后,计算出该次编队控制燃料消耗量Ek;2)第k次编队控制后,确定2个卫星轨道周期编队控制参数Δlk确定;3)将编队控制燃料累计消耗量之和SEk以及该次编队控制残差Δlk代入公式SEk≥k1&&Δlk≥k2进行判断;4)若公式成立,星上自主完成主辅星标志切换,先由地面上注指令或星上自主进行主辅星切换控制,使得新主星运行在参考轨道上,然后参数置初始值,接着新的辅星执行编队控制任务,重复步骤1)~3);若公式不成立,继续由原辅星进行编队控制。该方法工程可实现性强,节省燃料,同时兼顾编队卫星任务需求,能够提升提高卫星自主化水平。
Description
技术领域
本发明涉及一种编队卫星系统燃料均衡在轨实现方法。
背景技术
随着卫星技术的不断发展,越来越多的国家和组织都越来越倾向于利用多颗卫星组成编队卫星系统来实现空间任务目标。因此编队卫星系统控制任务越来越复杂,同时也面临着高控制精度、自主化以及长寿命等方面压力。
目前影响卫星寿命的决定性因素主要包括燃料以及电能,其中燃料消耗具有不可逆性。因此,由于编队卫星之间的相对关系需要通过编队控制消耗燃料实现,所以编队任务的顺利开展以及寿命、控制能力不仅取绝对编队控制精度,同时取决于编队卫星系统的燃料均衡管理。目前,编队卫星燃料均衡管理的相关研究成果一般与编队控制路径规划或者优化耦合在一起,同时燃料均衡管理方法多采用建立优化目标函数与复杂算法(如遗传算法、粒子群算法、匈牙利算法)相结合。根据实际在轨编队任务,存在约束条件较理想、算法复杂和缺乏在轨自主可实现性等不足。
发明内容
本发明解决的技术问题是:针对现有技术存在的不足,提供一种编队卫星燃料均衡管理工程可实现的方法,能够实现星上自主燃料均衡管理。
本发明的技术方案是:一种编队卫星系统燃料均衡在轨实现方法,包括如下步骤:
1)第k次编队控制后,卫星自主根据气瓶和推进剂贮箱内的压力和温度测量数据,利用气体状态方程计算得出推进剂贮箱内氦气体积,再由贮箱总容积和推进剂密度计算出贮箱内推进剂体积和质量,即PVT算法计算得出该次编队控制燃料消耗量Ek,k≥1;
2)第k次编队控制后,根据星载GNSS接收机星间测量数据,利用EKF/UKF滤波算法确定两个卫星轨道周期编队控制参数Δlk;
3)将编队控制燃料累计消耗量之和SEk以及该次编队控制残差Δlk代入下式:
SEk≥k1&&Δlk≥k2 (1)
其中,k1、k2为阈值,k1由星间面质比差异约束条件确定,k2由编队控制精度指标确定;SEk=E1+.......+Ek,“&&”为逻辑或关系;
4)若上式(1)成立,则星上自主完成主辅星标志切换,先由地面上注指令或星上自主进行主辅星切换控制,即主星、辅星绝对运行轨道切换控制,使得新主星运行在参考轨道上,然后参数置初始值,即k重新从1开始,SEk重新从0开始,接着新的辅星执行编队控制任务,重复步骤1)~步骤3);若上式(1)不成立,继续由原辅星进行编队控制;所述主星为参考星,辅星为跟随星。
本发明采用的方法与现有技术相比:
1、工程可实现性强,节省燃料,同时兼顾编队卫星任务需求。目前现有的编队任务燃料均衡方法仅是为了均衡而均衡,通过每次控制规划使得不同卫星扮演辅星,这样将导致参考轨道的地面轨迹出现跳变,不利于编队任务有效载荷的工作以及任务长期规划;另外,编队任务燃料均衡的根本原因不单纯是影响任务寿命,更重要的是由于燃料不均衡导致星间面质比差异变大,从而使得编队构形不稳定,因此燃料均衡的判断条件应该是依据构形稳定引出阈值判断,而非现在与单次控制相结合的均衡判断。实际上,工程上燃料均衡是存在一定裕度边界的,这样是现有技术的燃料均衡所忽略缺失的地方,因此更加消耗燃料;
2、提高卫星自主化水平卫星自主化运行是未来发展的趋势,能够克服传统的依赖地面系统控制方式需要大量的人力物力,运营成本高,处理故障和突发事件的能力差,不能胜任复杂的控制任务等不足,该方法能够在轨自主运行,是在轨自主精密编队的支撑技术与前提条件。
3、算法可实现、适用和可靠性强编队卫星系统燃料均衡管理全过程清晰,从相关判断量计算到通过主辅星切换控制实现燃料均衡,意义明确,无诸如遗传算法、粒子群算法等复杂算法,能够在轨自主实现;特别的,现有技术手段仅针对控制过程燃料消耗建议约束,本方法不仅考虑燃料消耗,同时考虑编队控制精度,因此当无法预估编队燃料消耗或者预估不准确的情况下,仍能进行燃料均衡管理,方法可靠性和适用性强,保证了编队任务运行的品质。
附图说明
以下将结合附图和实施例对本发明作进一步说明。
图1为本发明涉及的编队卫星系统燃料均衡管理工程化流程图;
图2为本发明涉及的Hill坐标系下编队参数示意图;
图3为本发明涉及的编队卫星系统主辅星切换示意图;
图4为本发明涉及的编队卫星系统全寿命燃料均衡实现过程示意图。
图5为本发明涉及的编队卫星系统燃料均衡管理双层触发条件示意图。
具体实施方式
星间燃料消耗的差异不仅影响编队卫星的寿命,更通过加剧大气阻力的摄动影响从而影响编队构形(编队卫星之间的空间相对关系)的稳定性以及编队控制的精度。因此,编队卫星系统需要进行星间的燃料均衡管理。针对现有技术存在的不足,本发明解决的技术问题是提供一种编队卫星燃料均衡在轨实现的方法,能够实现星上自主燃料均衡管理。
下面进一步对本发明的工作过程,整个方法流程如如图1所示,具体详见以下描述:
步骤1第k次(k≥1)编队控制后,卫星自主根据气瓶和推进剂贮箱内的压力和温度测量数据,利用气体状态方程计算出推进剂贮箱内氦气体积,再由贮箱总容积和推进剂密度计算出贮箱内推进剂体积和质量,即压力—体积—温度(PVT法)计算出该次编队控制燃料消耗量Ek;
其中,PVT法基本原理过程参考《合多种剩余推进剂测量方法的应用研究》。
步骤2第k次(k≥1)编队控制后,根据星载GNSS接收机星间测量数据利用EKF/UKF滤波算法确定运行2个卫星轨道周期编队控制参数Δlk确定;
其中,为了描述目辅星相对于主星的相对运动,需要定义基于开普勒轨道根数相对轨道根数矢量:
式中,下标1表示主星绝对轨道根数,下标2表示辅星绝对轨道根数。则Δl=a1(coti1Δiy+Δu)。
步骤3将编队控制燃料累计消耗量之和SEk以及该次编队控制残差Δlk代入下式:
SEk≥k1&&|Δlk|≥k2 (1)
其中,k1、k2为阈值,k1由星间面质比差异约束条件确定,一般在8kg~20kg之间,k2由编队控制精度指标约束确定,一般在20m~40m之间,具体值需要根据具体卫星状态确定,SEk=E1+.......+Ek,“&&”表示逻辑或关系。
步骤4若上式成立,星上自主完成主辅星(主星为参考星,辅星为跟随星)标志切换,先由地面上注指令或星上自主进行主辅星切换控制(即主星、辅星绝对运行轨道切换控制),使得新主星运行在参考轨道上,然后参数置初始值(包括k重新从1开始,SEk重新从0开始),接着新的辅星执行编队控制任务,重复1~3;若上式不成立,继续由原辅星进行编队控制。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种编队卫星系统燃料均衡在轨实现方法,其特征在于包括如下步骤:
1)第k次编队控制后,卫星自主根据气瓶和推进剂贮箱内的压力和温度测量数据,利用气体状态方程计算得出推进剂贮箱内氦气体积,再由贮箱总容积和推进剂密度计算出贮箱内推进剂体积和质量,即PVT算法计算得出该次编队控制燃料消耗量Ek,k≥1;
2)第k次编队控制后,根据星载GNSS接收机星间测量数据,利用EKF/UKF滤波算法确定两个卫星轨道周期编队控制参数Δlk;
其中,定义基于开普勒轨道根数相对轨道根数矢量:
式中,下标1表示主星绝对轨道根数,下标2表示辅星绝对轨道根数;则Δlk=a1(coti1Δiy+Δu);
3)将编队控制燃料累计消耗量之和SEk以及该次编队控制参数Δlk代入下式:
SEk≥k1&&Δlk≥k2 (1)
其中,k1、k2为阈值,k1由星间面质比差异约束条件确定,k2由编队控制精度指标确定;SEk=E1+…….+Ek,“&&”为逻辑或关系;
4)若上式(1)成立,则星上自主完成主辅星标志切换,先由地面上注指令或星上自主进行主辅星切换控制,即主星、辅星绝对运行轨道切换控制,使得新主星运行在参考轨道上,然后参数置初始值,即k重新从1开始,SEk重新从0开始,接着新的辅星执行编队控制任务,重复步骤1)~步骤3);若上式(1)不成立,继续由原辅星进行编队控制;所述主星为参考星,辅星为跟随星。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611097493.0A CN106681138B (zh) | 2016-12-02 | 2016-12-02 | 一种编队卫星系统燃料消耗均衡在轨实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611097493.0A CN106681138B (zh) | 2016-12-02 | 2016-12-02 | 一种编队卫星系统燃料消耗均衡在轨实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106681138A CN106681138A (zh) | 2017-05-17 |
CN106681138B true CN106681138B (zh) | 2019-06-21 |
Family
ID=58867478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611097493.0A Active CN106681138B (zh) | 2016-12-02 | 2016-12-02 | 一种编队卫星系统燃料消耗均衡在轨实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106681138B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110096069B (zh) * | 2019-04-25 | 2020-07-28 | 南京航空航天大学 | 基于nsga ii超幅宽成像卫星编队构型的优化方法 |
CN110377047B (zh) * | 2019-06-03 | 2022-08-02 | 上海航天控制技术研究所 | 一种双星卫星编队防碰撞方法 |
CN111443722B (zh) * | 2020-03-23 | 2022-09-23 | 上海航天控制技术研究所 | 一种编队卫星定时段自主保持方法 |
CN114237040A (zh) * | 2021-11-17 | 2022-03-25 | 航天东方红卫星有限公司 | 一种卫星燃料余量在轨实时自主计算与控制系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2674047A1 (fr) * | 1991-03-15 | 1992-09-18 | Thomson Csf | Calculateur de navigation multisenseurs a filtre de kalman modulaire. |
FR2937735B1 (fr) * | 2008-10-24 | 2010-12-24 | Thales Sa | Systeme et procede de detection de bulles presentes dans un fluide circulant a l'interieur d'un conduit |
CN101445160A (zh) * | 2008-12-23 | 2009-06-03 | 中国科学院广州能源研究所 | 用于微纳卫星的微型液体燃料推进系统 |
CN102122171B (zh) * | 2010-12-28 | 2013-04-17 | 北京航空航天大学 | 一种基于智能移动机器人的多微纳探测器组网联合演示验证系统 |
CN102354123B (zh) * | 2011-07-18 | 2013-05-08 | 北京航空航天大学 | 一种跨平台可扩展的卫星动态仿真测试系统 |
CN102541070A (zh) * | 2012-01-20 | 2012-07-04 | 哈尔滨工业大学 | 一种卫星编队飞行地面试验系统的碰撞规避方法 |
CN103257653B (zh) * | 2013-05-22 | 2016-04-06 | 上海新跃仪表厂 | 一种基于燃料消耗优化的卫星编队构形控制方法 |
CN104330971B (zh) * | 2014-10-28 | 2018-01-23 | 蔡远文 | 微小卫星群编队消耗量优化方法 |
CN104317306B (zh) * | 2014-10-28 | 2017-07-21 | 蔡远文 | 卫星编队重构方法 |
CN105737834B (zh) * | 2014-12-09 | 2018-06-26 | 上海新跃仪表厂 | 一种基于轨道平根数的相对导航鲁棒滤波方法 |
CN105891859B (zh) * | 2016-05-23 | 2019-01-15 | 成都市精准时空科技有限公司 | 一种卫星导航pvt解算方法以及相应地芯片、模块 |
-
2016
- 2016-12-02 CN CN201611097493.0A patent/CN106681138B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106681138A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106681138B (zh) | 一种编队卫星系统燃料消耗均衡在轨实现方法 | |
Chen et al. | Optimal positioning of flying relays for wireless networks: A LOS map approach | |
CN105883008B (zh) | 卫星推力器布局方法 | |
Liu et al. | Active disturbance rejection control for delayed electromagnetic docking of spacecraft in elliptical orbits | |
CN110426029A (zh) | 用于无人机蜂群协同导航的动态互观测在线建模方法 | |
CN103731221B (zh) | 一种天地一体化网络系统可用性确定方法 | |
CN110276159A (zh) | 一种基于多模型融合的卫星系统多学科优化方法 | |
CN103112601B (zh) | 导航geo卫星与静止通信卫星共位方法 | |
CN101101650B (zh) | 以机动过载为控制量的低空突防导弹三维航迹规划方法 | |
Fang et al. | Emergency landing trajectory optimization for fixed-wing UAV under engine failure | |
CN107566026A (zh) | 多层次星座组网的卫星信息网络 | |
CN107820260A (zh) | 平动点基干网络 | |
Luo et al. | Skip entry guidance using numerical predictor–corrector and patched corridor | |
Gao et al. | Fast path planning for firefighting UAV based on A-Star algorithm | |
Dell’Elce et al. | Comparison between analytical and optimal control techniques in the differential drag based rendez-vous | |
CN106672269A (zh) | 纳卫星电磁对接重构设备的设计方法 | |
Chasset et al. | 3-axis magnetic control with multiple attitude profile capabilities in the PRISMA mission | |
Thyagarajan et al. | University small satellite program—ANUSAT | |
Park et al. | A new approach to on-board stationkeeping of geo-satellites | |
Luo et al. | Blended skip entry guidance for low-lifting lunar return vehicles | |
Li et al. | Research on Cooperative Navigation Algorithm for UAV Cluster Based on Belief Propagation | |
Qingzhen et al. | Reentry trajectory planning optimization based on ant colony algorithm | |
CN108120993A (zh) | 一种基于对流层先验信息约束的卫星钟差实时测定方法 | |
Li et al. | An optimal energy absorption and utilization design for deformable/reconfigurable modular solar-powered uavs in near space | |
Du et al. | Dynamic coverage performance of wind-assisted balloons mesh based on Voronoi partition and energy constraint |
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 |