CN103678814A - 临界倾角近圆轨道的偏心率预偏置设计方法 - Google Patents
临界倾角近圆轨道的偏心率预偏置设计方法 Download PDFInfo
- Publication number
- CN103678814A CN103678814A CN201310699927.4A CN201310699927A CN103678814A CN 103678814 A CN103678814 A CN 103678814A CN 201310699927 A CN201310699927 A CN 201310699927A CN 103678814 A CN103678814 A CN 103678814A
- Authority
- CN
- China
- Prior art keywords
- mean
- orbit
- excentricity
- taken
- prebias
- 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.)
- Granted
Links
Images
Landscapes
- Braking Arrangements (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明提供了一种临界倾角近圆轨道的偏心率预偏置设计方法,属于卫星轨道技术领域。本方法首先将偏心率、近地点幅角、平近点角、升交点赤经、轨道半长轴和轨道倾角六个参数的平根数转换为瞬根数,然后将瞬根数带入STK卫星仿真软件中进行轨道计算,得到寿命期末时六个参数的瞬根数,再将瞬根数转换为平根数,最后将偏心率的平根数预偏置量设计为寿命期末时偏心率平根数的一半,近地点幅角的平根数设计为180°。本发明可以将寿命期末的偏心率减少一半,以克服临界倾角近圆轨道偏心率逐渐变大的不足。
Description
技术领域
本发明涉及卫星轨道技术领域,具体涉及一种临界倾角近圆轨道的偏心率预偏置设计方法。
背景技术
卫星在轨飞行期间,地球非球型摄动会带来轨道的长期项变化;在卫星寿命期内,为避免轨道根数超出阈值,需要施加必要的控制予以抑制。常用的控制策略包括:1)根据根数的实时测量值进行多次控制,对于寿命比较长或轨道稳定性要求比较高的任务可以采用这种策略;2)针对某根数进行预偏置设计而在整个寿命期间不再进行调整,对于寿命不太长或允许变化范围较大的任务可以采用这种策略。
对于第2种策略,常应用于太阳同步轨道抑制降交点地方时漂移而采取倾角预偏置[参见:杨维廉.太阳同步回归轨道的长期演变与控制,航天器工程,2008年3月,17(2):26-30]。半长轴用以调整星下点轨迹或相位,一般不进行预偏置设计。而对于非临界倾角轨道,偏心率会以冻结偏心率为中心进行往复变化(以645km高度太阳同步轨道为例,往复周期约为100天),一般不进行预偏置设计。以往发射的临界倾角轨道以大椭圆轨道为主,例如前苏联Molniya轨道偏心率达到0.5以上,且保持不变,无须进行预偏置设计;而以往发射的临界倾角近圆轨道均以返回式卫星为主,其在轨寿命较短,也不进行预偏置设计。
在非球型摄动下,临界倾角大椭圆轨道可以保持偏心率和近地点幅角不变;而对于临界倾角近圆轨道,非球型摄动将造成偏心率缓慢增加,甚至超出近圆轨道所容忍的范围。考察高度为1100km、倾角为63.4°轨道的偏心率和近地点幅角(或偏心率矢量)的长期演化规律,如图1所示。偏心率矢量(ex,ey)将在100年内完成长周期变化;且在(0,0)附近偏心率矢量的变化接近线性。偏心率矢量定义如下:ex=e·cosω,ey=e·sinω,其中,e为偏心率,ω为近地点幅角。
以高度为1100km、倾角为63.4°圆轨道为初始轨道,在卫星寿命期3年内偏心率由0变化为0.0113326,且偏心率变化以近似线性方式增加,如图2所示。但地球遥感或侦察任务不希望卫星在飞行一圈内轨道高度变化过大,且偏心率越小越好。
发明内容
本发明为了克服临界倾角近圆轨道偏心率逐渐变大的不足,提供了一种临界倾角近圆轨道的偏心率预偏置设计方法,本发明方法可减小寿命期内的偏心率。
本发明的临界倾角近圆轨道的偏心率预偏置设计方法,包括如下步骤:
步骤1:将偏心率、近地点幅角、平近点角、升交点赤经、轨道半长轴和轨道倾角六个参数的平根数转换为瞬根数;偏心率的平根数取为0,近地点幅角的平根数取为0°,平近点角的平根数取为0°,升交点赤经取为任意值,轨道倾角的平根数为63.4°;
步骤2:把六个参数的瞬根数带入STK卫星仿真软件中进行轨道计算,得到寿命期末时六个参数的瞬根数;
步骤3:将寿命期末时六个参数的瞬根数转换为平根数,得到偏心率的平根数ef;
步骤4:将偏心率的平根数预偏置量设计为ef/2,近地点幅角的平根数设计为180°。
本发明与现有技术相比,本发明通过预偏置偏心率和近地点幅角,可以将寿命期末的偏心率减少一半,以克服临界倾角近圆轨道偏心率逐渐变大的不足。
附图说明
图1为偏心率矢量的长期演化曲线图;
图2为预偏置前偏心率的长期演化曲线示意图;
图3为本发明的偏心率预偏置设计方法的流程图;
图4为预偏置后偏心率的长期演化曲线示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
以3年寿命期、高度为1100km的临界倾角轨道为例,结合图3来说明本发明的偏心率预偏置设计方法。
步骤1:将偏心率、近地点幅角、平近点角、升交点赤经、轨道半长轴和轨道倾角的平根数转换为瞬根数。
初始六个参数的平根数设置为:偏心率取为0,近地点幅角取为0°,平近点角的平根数取为0°,升交点赤经取为任意值,半长轴取为7471.004km,轨道倾角取为63.4°,即为临界倾角(Critical Inclination)。将六个参数的平根数代入由平根到瞬根的转换公式,得到六个参数的瞬根数,本发明实施例得到:瞬半长轴为7478.074km,瞬偏心率为0.00055262,瞬倾角为63.4236°,瞬近地点幅角为0°,瞬平近点角为0°,瞬升交点赤经为任意值。
步骤2:把转换所得到的六个参数的瞬根数带入STK软件进行轨道计算,得到寿命期末六个参数的瞬根数。
在STK卫星仿真软件中的具体操作是:将轨道预报法“Propagator”取为“LOP”,开始时间“Start Time”取为任意,轨道历元时刻“Orbit Epoch”与“Start Time”相同,结束时间“Stop Time”为“Start Time”加上卫星寿命期,参数“Cood Type”取为“Classical”,参数“Cood System”取为“TEMEOfDate”。LOP表示长期轨道预报,Classical表示经典型轨道根数(包括半长轴、偏心率、倾角、近地点幅角、平近点角和升交点赤经),TEMEOfDate表示真赤道平春分点惯性坐标系。Cood Type表示所选取的参数类型,Cood System表示所选取的坐标系类型。
本发明实施例中得到寿命期末六个参数的瞬根数为:瞬半长轴为7469.759km,瞬偏心率为0.011805,瞬倾角为63.386°,瞬近地点幅角为1.516°,瞬平近点角为127.149°,瞬升交点赤经为任意值。
步骤3:将计算得到的寿命期末时六个参数的瞬根数,代入瞬根到平根的转换公式,得到六个参数的平根数。本发明中需要获取寿命期末时偏心率的平根数ef。
本发明实施例中获得寿命期末偏心率的平根数ef为0.0113326。
步骤4:根据步骤3得到的寿命期末偏心率设计偏心率预偏置量。
本发明实施例中,将偏心率的平根数预偏置量设计为0.005666,近地点幅角的平根数设计为180°。
本发明实施例中,利用本发明偏心率预偏置设计方法进行偏心率预偏置量设置后的作用下,高度为1100km临界倾角轨道在3年寿命期内的偏心率变化如图3所示;从图3中可以看出,寿命期内,偏心率的最大值不超过0.005666,克服了临界倾角近圆轨道偏心率逐渐变大的不足。
Claims (2)
1.一种临界倾角近圆轨道的偏心率预偏置设计方法,其特征在于,包括如下步骤:
步骤1:将偏心率、近地点幅角、平近点角、升交点赤经、轨道半长轴和轨道倾角六个参数的平根数转换为瞬根数;偏心率的平根数取为0,近地点幅角的平根数取为0°,平近点角的平根数取为0°,升交点赤经取为任意值,轨道倾角的平根数为63.4°;
步骤2:把六个参数的瞬根数带入STK卫星仿真软件中进行轨道计算,得到寿命期末时六个参数的瞬根数;
步骤3:将寿命期末时六个参数的瞬根数转换为平根数,得到偏心率的平根数ef;
步骤4:将偏心率的平根数预偏置量设计为ef/2,近地点幅角的平根数设计为180°。
2.根据权利要求1所述的一种临界倾角近圆轨道的偏心率预偏置设计方法,其特征在于,所述的步骤2中,在STK卫星仿真软件中的具体操作是:将轨道预报法“Propagator”取为“LOP”,开始时间“Start Time”取为任意,轨道历元时刻“Orbit Epoch”与“Start Time”相同,结束时间“Stop Time”为“Start Time”加上卫星寿命期,参数类型“Cood Type”取为“Classical”,坐标系类型“Cood System”取为“TEMEOfDate”;LOP表示长期轨道预报,Classical表示经典型轨道根数,TEMEOfDate表示真赤道平春分点惯性坐标系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310699927.4A CN103678814B (zh) | 2013-12-18 | 2013-12-18 | 临界倾角近圆轨道的偏心率预偏置设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310699927.4A CN103678814B (zh) | 2013-12-18 | 2013-12-18 | 临界倾角近圆轨道的偏心率预偏置设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103678814A true CN103678814A (zh) | 2014-03-26 |
CN103678814B CN103678814B (zh) | 2016-06-29 |
Family
ID=50316350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310699927.4A Active CN103678814B (zh) | 2013-12-18 | 2013-12-18 | 临界倾角近圆轨道的偏心率预偏置设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103678814B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106092105A (zh) * | 2016-06-03 | 2016-11-09 | 上海航天控制技术研究所 | 一种近地卫星严格回归轨道的确定方法 |
CN106250623A (zh) * | 2016-08-01 | 2016-12-21 | 上海航天控制技术研究所 | 一种基于状态平稳切换的半物理快速仿真方法 |
CN109211225A (zh) * | 2017-06-29 | 2019-01-15 | 中国科学院国家天文台 | 获取大椭圆轨道空间物体剩余轨道寿命的方法、系统及设备 |
CN109613574A (zh) * | 2018-11-13 | 2019-04-12 | 中国人民解放军战略支援部队航天工程大学 | 计算北斗中轨卫星坟墓轨道最早穿越其他全球卫星导航系统轨道时间的方法 |
CN110060552A (zh) * | 2019-03-26 | 2019-07-26 | 南京航空航天大学 | 一种卫星星座轨道物理仿真实验平台 |
CN110068845A (zh) * | 2019-04-30 | 2019-07-30 | 上海微小卫星工程中心 | 一种基于平根数理论确定卫星理论轨道的方法 |
CN110068846A (zh) * | 2019-04-30 | 2019-07-30 | 上海微小卫星工程中心 | 一种基于星载gnss接收机在星上自主确定轨道平根数的方法 |
CN113778112A (zh) * | 2021-08-13 | 2021-12-10 | 北京九天微星科技发展有限公司 | 一种严格回归轨道平面外管道控制的优化方法及装置 |
CN115149995A (zh) * | 2022-05-16 | 2022-10-04 | 亚太卫星通信(深圳)有限公司 | 一种heo星座轨道设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3914086B2 (ja) * | 2002-04-17 | 2007-05-16 | 株式会社国際電気通信基礎技術研究所 | 生体分子計算システム設計システムおよび設計方法 |
CN101010675A (zh) * | 2005-04-12 | 2007-08-01 | 松下电器产业株式会社 | 电子元件设计方法 |
CN102495950A (zh) * | 2011-11-25 | 2012-06-13 | 北京航空航天大学 | 一种适用于太阳同步轨道的倾角偏置量获取方法 |
EP2665007A1 (en) * | 2012-05-16 | 2013-11-20 | Sumitomo Rubber Industries, Ltd. | Simulation method for polymer material comprising filler |
-
2013
- 2013-12-18 CN CN201310699927.4A patent/CN103678814B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3914086B2 (ja) * | 2002-04-17 | 2007-05-16 | 株式会社国際電気通信基礎技術研究所 | 生体分子計算システム設計システムおよび設計方法 |
CN101010675A (zh) * | 2005-04-12 | 2007-08-01 | 松下电器产业株式会社 | 电子元件设计方法 |
CN102495950A (zh) * | 2011-11-25 | 2012-06-13 | 北京航空航天大学 | 一种适用于太阳同步轨道的倾角偏置量获取方法 |
EP2665007A1 (en) * | 2012-05-16 | 2013-11-20 | Sumitomo Rubber Industries, Ltd. | Simulation method for polymer material comprising filler |
Non-Patent Citations (3)
Title |
---|
侯育卓: "近圆轨道编队飞行星座相对构型稳定性分析", 《航天控制》, no. 1, 29 February 2004 (2004-02-29) * |
安雪滢: "椭圆轨道航天器编队飞行动力学及应用研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》, no. 5, 15 November 2007 (2007-11-15) * |
徐明 等: "绕月飞行的大幅值逆行轨道研究", 《宇航学报》, no. 5, 30 September 2009 (2009-09-30) * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106092105A (zh) * | 2016-06-03 | 2016-11-09 | 上海航天控制技术研究所 | 一种近地卫星严格回归轨道的确定方法 |
CN106250623A (zh) * | 2016-08-01 | 2016-12-21 | 上海航天控制技术研究所 | 一种基于状态平稳切换的半物理快速仿真方法 |
CN106250623B (zh) * | 2016-08-01 | 2019-05-24 | 上海航天控制技术研究所 | 一种基于状态平稳切换的半物理快速仿真方法 |
CN109211225A (zh) * | 2017-06-29 | 2019-01-15 | 中国科学院国家天文台 | 获取大椭圆轨道空间物体剩余轨道寿命的方法、系统及设备 |
CN109613574B (zh) * | 2018-11-13 | 2023-05-12 | 中国人民解放军战略支援部队航天工程大学 | 计算北斗中轨卫星坟墓轨道最早穿越其他全球卫星导航系统轨道时间的方法 |
CN109613574A (zh) * | 2018-11-13 | 2019-04-12 | 中国人民解放军战略支援部队航天工程大学 | 计算北斗中轨卫星坟墓轨道最早穿越其他全球卫星导航系统轨道时间的方法 |
CN110060552A (zh) * | 2019-03-26 | 2019-07-26 | 南京航空航天大学 | 一种卫星星座轨道物理仿真实验平台 |
CN110068845A (zh) * | 2019-04-30 | 2019-07-30 | 上海微小卫星工程中心 | 一种基于平根数理论确定卫星理论轨道的方法 |
CN110068845B (zh) * | 2019-04-30 | 2021-07-23 | 上海微小卫星工程中心 | 一种基于平根数理论确定卫星理论轨道的方法 |
CN110068846B (zh) * | 2019-04-30 | 2022-01-07 | 上海微小卫星工程中心 | 一种基于星载gnss接收机在星上自主确定轨道平根数的方法 |
CN110068846A (zh) * | 2019-04-30 | 2019-07-30 | 上海微小卫星工程中心 | 一种基于星载gnss接收机在星上自主确定轨道平根数的方法 |
CN113778112A (zh) * | 2021-08-13 | 2021-12-10 | 北京九天微星科技发展有限公司 | 一种严格回归轨道平面外管道控制的优化方法及装置 |
CN113778112B (zh) * | 2021-08-13 | 2024-05-14 | 北京九天微星科技发展有限公司 | 一种严格回归轨道平面外管道控制的优化方法及装置 |
CN115149995A (zh) * | 2022-05-16 | 2022-10-04 | 亚太卫星通信(深圳)有限公司 | 一种heo星座轨道设计方法 |
CN115149995B (zh) * | 2022-05-16 | 2023-11-24 | 亚太卫星通信(深圳)有限公司 | 一种heo星座轨道设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103678814B (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103678814A (zh) | 临界倾角近圆轨道的偏心率预偏置设计方法 | |
CN102424116B (zh) | 一种静止轨道卫星变轨策略优化方法 | |
CN105883008B (zh) | 卫星推力器布局方法 | |
Lu et al. | Highly constrained optimal launch ascent guidance | |
CN103257653B (zh) | 一种基于燃料消耗优化的卫星编队构形控制方法 | |
CN104828261B (zh) | 一种小型动量轮的角动量卸载方法 | |
Rasheed | On Solving Hyperbolic Trajectory Using New Predictor-Corrector Quadrature Algorithms | |
CN107168372B (zh) | 基于在轨参数辨识和偏置的卫星跟飞长期摄动补偿方法 | |
CN103678787B (zh) | 一种星下点圆迹地球同步轨道设计方法 | |
CN104076819B (zh) | 一种圆参考轨道下卫星有界伴飞的边界控制方法 | |
Vtipil et al. | Determining an Earth observation repeat ground track orbit for an optimization methodology | |
CN105353621A (zh) | 一种地球静止轨道卫星电推力器故障模式推力分配方法 | |
CN103853887B (zh) | 一种冻结轨道的偏心率的卫星轨道确定方法 | |
CN109539903A (zh) | 一种固体运载火箭椭圆转移轨道迭代制导控制方法 | |
CN103112600A (zh) | 一种星际转移轨道设计方法 | |
RU2608186C2 (ru) | Способ и система для управления группой, по меньшей мере, из двух спутников, выполненных с возможностью обеспечения обслуживания | |
CN104309822A (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN106840164B (zh) | 一种多碎片主动清除在线重规划方法 | |
JP2015000703A (ja) | 宇宙機の軌道制御方法およびその装置 | |
CN106815400B (zh) | 一种调轨方案自动化设计方法 | |
Xia et al. | Adaptive control for spacecraft rendezvous subject to time-varying inertial parameters and actuator faults | |
CN103235870B (zh) | 兼顾多任务高度的太阳同步轨道倾角偏置方法 | |
RU2013131329A (ru) | Способ автономной коллокации на геостационарной орбите | |
Kos et al. | Altair descent and ascent reference trajectory design and initial dispersion analyses | |
CN103853047A (zh) | 一种基于状态量反馈的小推力跟踪制导方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |