CN114397906B - 地球静止卫星电推进转移快速高精度计算方法 - Google Patents
地球静止卫星电推进转移快速高精度计算方法 Download PDFInfo
- Publication number
- CN114397906B CN114397906B CN202111487363.9A CN202111487363A CN114397906B CN 114397906 B CN114397906 B CN 114397906B CN 202111487363 A CN202111487363 A CN 202111487363A CN 114397906 B CN114397906 B CN 114397906B
- Authority
- CN
- China
- Prior art keywords
- control law
- stage
- inclination angle
- angle
- electric propulsion
- 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
- 238000012546 transfer Methods 0.000 title claims abstract description 23
- 238000004364 calculation method Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 14
- 230000010354 integration Effects 0.000 claims abstract description 9
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- 238000013461 design Methods 0.000 claims abstract description 4
- 230000003247 decreasing effect Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开的地球静止卫星电推进转移快速高精度计算方法，包括确定控制设计变量偏航角，建立法向控制律；降低倾角，建立第一阶段面内控制律来改变半长轴，直到半长轴达到同步轨道高度，本阶段数值积分终止转入下一步骤；继续降低倾角，同时建立第二阶段面内控制律来改变偏心率，直到偏心率达到0，本阶段数值积分终止；改变偏航角来更新法向控制律，来调整改变倾角的变化速率，直到倾角在终止时刻等于0，则表示迭代收敛，即半长轴、偏心率和倾角在数值积分终端时刻同时满足同步轨道条件，最后输出结果。本发明的地球静止卫星电推进转移快速高精度计算方法，控制律具有简明的解析解，控制过程简单，便于工程实施，和最优解相差较小。
Description
技术领域
本发明属于航天动力学与控制技术领域，具体涉及一种地球静止卫星电推进转移快速高精度计算方法。
背景技术
电推进技术的成熟发展，电推进技术的应用逐渐从位置保持向完成整个轨道转移控制任务转变。以往传统的做法是火箭将静止卫星发射到地球同步转移轨道(GTO)以后，通常采用化学推进发动机在远地点进行3-5次变轨可以将GTO轨道变成地球静止轨道(GEO)。如果用电推进代替化学推进，由于电推进高比冲的特点，可以节省大量燃料。但电推进系统推力很小，一般只有几百毫牛，远远小于化学推进的几百牛，因而采用电推进从GTO变轨至GEO的时间将长达数月，且变轨策略与化学推进的策略有很大区别。
地球静止卫星轨道转移是一类经典问题，但以往的研究大多针对的情况是推力比较大，轨道转移圈次只有几十圈的情况。以目前电推进系统的水平，推力还不能达到几牛的量级。由于初值估计偏差随时间的传播和累积，一般来说轨道转移时间越长、圈次越多，优化难度越大。Graham和Rao采用伪谱法优化多圈静止卫星转移轨道，伪谱法本质上是配点法的一种，解的最优性只能事后验证。很多学者将转移轨道根数进行均化处理，这样得到的解对估计时间和燃料消耗固然有意义，但也会一定程度上降低轨道动力学模型的准确性。国内学者田百义等人基于李雅普诺夫优化原理进行了GTO-GEO转移轨道的优化设计，但和最优解相差较大。
发明内容
本发明的目的在于提供一种地球静止卫星电推进转移快速高精度计算方法，解决了同步卫星全电推进入轨的多圈轨道优化难度大的问题。
本发明所采用的技术方案是：地球静止卫星电推进转移快速高精度计算方法，包括以下步骤：
步骤1、确定控制设计变量偏航角，建立法向控制律；
步骤2、采用步骤1的法向控制律降低倾角，建立第一阶段面内控制律来改变半长轴，直到半长轴达到同步轨道高度，本阶段数值积分终止转入下一步骤；
步骤3、采用步骤1的法向控制律继续降低倾角，同时建立第二阶段面内控制律来改变偏心率，直到偏心率达到0，本阶段数值积分终止；
步骤4、改变偏航角来更新步骤1的法向控制律，来调整改变倾角的变化速率，直到倾角在步骤3终止时刻等于0，则表示迭代收敛，即半长轴、偏心率和倾角在数值积分终端时刻同时满足同步轨道条件，最后输出结果。
本发明的特点还在于，
步骤1中建立的法向控制律表示为：
式(1)中，F表示推力大小，m表示卫星质量，β表示偏航角，取值范围0到90度，u＝ω+ν表示纬度幅角，ω表示近地点幅角，ν表示真近点角，sign表示符号函数。
步骤2中建立的第一阶段面内控制律中径向和横向控制律分别表示为：
fr＝0 (2)
步骤3中建立的第二阶段面内控制律中径向和横向控制律分别表示为：
步骤4具体包括一维搜索β满足终端时刻倾角为0的条件。
本发明的有益效果是：本发明的地球静止卫星电推进转移快速高精度计算方法，控制律具有简明的解析解，控制过程简单，便于工程实施，和最优解相差较小。
附图说明
图1是本发明的地球静止卫星电推进转移快速高精度计算方法的流程图。
具体实施方式
下面结合附图以及具体实施方式对本发明进行详细说明。
本发明提供了一种地球静止卫星电推进转移快速高精度计算方法，如图1所示，具体按照以下步骤实施：
步骤1、建立法向控制律
其中F表示推力大小，m表示卫星质量，β表示偏航角(取值范围0到90度)，u＝ω+ν表示纬度幅角，ω表示近地点幅角，ν表示真近点角，sign表示符号函数。
步骤2、建立面内控制律改变半长轴，径向和横向控制律分别为
fr＝0 (2)
步骤3、建立面内控制律改变偏心率，径向和横向控制律分别为
步骤4、一维搜索β满足终端时刻倾角为0的条件。
实施例
假设同步卫星入轨的初始条件如下：近地点高度hp0＝200km、远地点高度ha0＝10000km、倾角i0＝20.8°、升交点赤经Ω0＝0°、近地点幅角ω＝180°、平近点角M＝2°。可见，偏心率大约是0.73，星箭分离后卫星入轨远地点高度等于同步轨道高度35786km。卫星初始质量是4000kg，比冲3500s，推力350mN。因此初始加速度是0.875×10-4m/s2。
采用本发明的地球静止卫星电推进转移快速高精度计算方法，获得最优343.7天的次优解，包含463圈的轨道转移，采用偏航角β＝32.8°。采用同样的初始条件和最省时间的性能指标，李雅普诺夫方法获得指标是338.9天，本发明的方法和李雅普诺夫方法指标仅相差1.4％，但是本发明方法的控制律更便于工程实施。
通过上述方式，本发明的地球静止卫星电推进转移快速高精度计算方法，控制律具有简明的解析解，控制过程简单，便于工程实施，和最优解相差较小。
Claims (2)
1.地球静止卫星电推进转移快速高精度计算方法，其特征在于，包括以下步骤：
步骤1、确定控制设计变量偏航角，建立法向控制律，表示为：
式(1)中，F表示推力大小，m表示卫星质量，β表示偏航角，取值范围0到90度，u＝ω+ν表示纬度幅角，ω表示近地点幅角，ν表示真近点角，sign表示符号函数；
步骤2、采用步骤1的法向控制律降低倾角，建立第一阶段面内控制律来改变半长轴，直到半长轴达到同步轨道高度，本阶段数值积分终止转入下一步骤；建立的第一阶段面内控制律中径向和横向控制律分别表示为：
fr＝0 (2)
步骤3、采用步骤1的法向控制律继续降低倾角，同时建立第二阶段面内控制律来改变偏心率，直到偏心率达到0，本阶段数值积分终止；建立的第二阶段面内控制律中径向和横向控制律分别表示为：
步骤4、改变偏航角来更新步骤1的法向控制律，来调整改变倾角的变化速率，直到倾角在步骤3终止时刻等于0，则表示迭代收敛，即半长轴、偏心率和倾角在数值积分终端时刻同时满足同步轨道条件，最后输出结果。
2.如权利要求1所述的地球静止卫星电推进转移快速高精度计算方法，其特征在于，所述步骤4具体包括一维搜索β满足终端时刻倾角为0的条件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111487363.9A CN114397906B (zh) | 2021-12-06 | 2021-12-06 | 地球静止卫星电推进转移快速高精度计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111487363.9A CN114397906B (zh) | 2021-12-06 | 2021-12-06 | 地球静止卫星电推进转移快速高精度计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114397906A CN114397906A (zh) | 2022-04-26 |
| CN114397906B true CN114397906B (zh) | 2023-06-06 |
Family
ID=81227025
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111487363.9A Active CN114397906B (zh) | 2021-12-06 | 2021-12-06 | 地球静止卫星电推进转移快速高精度计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114397906B (zh) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106628261B (zh) * | 2016-10-20 | 2018-11-20 | 上海航天控制技术研究所 | 一种电推位保过程中的卫星高精度高稳度姿态控制方法 |
| CN106802125B (zh) * | 2017-01-24 | 2020-02-11 | 上海空间推进研究所 | 航天器用电推力器推力矢量偏心的确定方法 |
| CN110304279B (zh) * | 2019-05-31 | 2021-03-26 | 北京控制工程研究所 | 一种电推进卫星的质心在轨标定补偿方法 |
| CN112506211B (zh) * | 2020-12-07 | 2022-08-12 | 上海卫星工程研究所 | 面向未来重力场测量的分离式卫星平台无拖曳控制方法及系统 |
-
2021
- 2021-12-06 CN CN202111487363.9A patent/CN114397906B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114397906A (zh) | 2022-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110276159B (zh) | 一种基于多模型融合的卫星系统多学科优化方法 | |
| CN105631095A (zh) | 一种等间隔发射的多约束地月转移轨道簇搜索方法 | |
| CN108491668A (zh) | 一种基于动态多模型融合的飞行器系统优化方法 | |
| CN114933028A (zh) | 双星轨控策略控制方法、装置、电子设备及存储介质 | |
| CN111366984B (zh) | 一种基于重力卫星星间激光测距系统确定引力场模型的方法 | |
| CN108762285A (zh) | 一种航天器多级复合控制的目标姿态协同规划方法及系统 | |
| CN112607065B (zh) | 一种基于电推进系统的高精度相位控制方法 | |
| CN112016187B (zh) | 一种基于混合动力的近地小行星交会任务轨道优化方法 | |
| EP0438229A2 (en) | Method for controlling east/west motion of a geostationary satellite | |
| CN111301715A (zh) | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 | |
| CN103226631A (zh) | 一种小推力转移轨道快速设计与优化方法 | |
| CN110378012A (zh) | 一种考虑高阶重力场的严格回归轨道设计方法 | |
| CN114715435A (zh) | 一种解析的小推力圆轨道异面交会优化方法 | |
| CN115373423A (zh) | 一种用于商业卫星的编队捕获方法 | |
| CN114397906B (zh) | 地球静止卫星电推进转移快速高精度计算方法 | |
| CN103235870B (zh) | 兼顾多任务高度的太阳同步轨道倾角偏置方法 | |
| CN109325288A (zh) | 一种基于不确定性优化的固体运载器总体参数确定方法及系统 | |
| CN113093246A (zh) | 地面多目标点成像快速判定及任务参数计算方法 | |
| CN105912020A (zh) | 一种弱稳定边界区域内航天器编队飞行方法 | |
| Barbosa et al. | Multidisciplinary design optimization of sounding rocket fins shape using a tool called MDO-SONDA | |
| CN104950668A (zh) | 卫星编队解析式燃料优化控制方法及系统 | |
| CN107844462A (zh) | 一种行星际连续推力转移轨道评估方法 | |
| CN108984903B (zh) | 一种制导参数的优选/优化设计方法 | |
| CN116466576A (zh) | 一种运载火箭末级在线规划动力学约束的伪谱凸化方法 | |
| CN115993777A (zh) | 基于轨道摄动模型反演的径切联控解耦迭代标定方法 |
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 |