CN107633142A - 相对运动轨道构型的模拟方法 - Google Patents
相对运动轨道构型的模拟方法 Download PDFInfo
- Publication number
- CN107633142A CN107633142A CN201710868984.9A CN201710868984A CN107633142A CN 107633142 A CN107633142 A CN 107633142A CN 201710868984 A CN201710868984 A CN 201710868984A CN 107633142 A CN107633142 A CN 107633142A
- Authority
- CN
- China
- Prior art keywords
- star
- relative motion
- perigee
- speed
- argument
- 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
Abstract
本发明公开了一种相对运动轨道构型的模拟方法,其包括如下步骤:步骤一,设计低轨轨道的近地点高度和远地点高度,满足不同的交会速度如100m/s,200m/s,及300m/s;步骤二,设计合理的近地点幅角,使得双星椭圆轨道交会点在我国上空,满足测控弧段要求;步骤三,设计两星距离,使得第一星在相对运动椭圆半长轴的合适位置,可调整交会时刻的两星距离;步骤四,拟定试验步骤,从变轨、跟踪、指向、试验的先后顺序,分解试验内容等。本发明遵行先易后难的原则,双星运行在低轨轨道,创新性的利用近地点幅角相差半周的椭圆轨道,调节双椭圆不同的远地点高度获得不同的相对运动速度,调节双星在两轨道交会处的距离获得不同的相对运动角速度。
Description
技术领域
本发明涉及一种模拟方法,特别是涉及一种相对运动轨道构型的模拟方法。
背景技术
高轨轨道在轨碎片清除任务包括探测跟踪和识别太空碎片、监测太空环境、对碎片进行抓捕清除,为了高效的抵近到碎片目标附近,设想利用零倾角椭圆轨道对高轨碎片等目标进行交会访问,在椭圆轨道的远地点接近高轨目标,先期在低轨开展模拟高轨交会的验证试验,利用低轨共面双椭圆交会轨道模拟椭圆轨道对高轨轨道的相对交会运动工况,验证对目标的大动态测量与轨道预报等技术。
发明内容
本发明所要解决的技术问题是提供一种相对运动轨道构型的模拟方法,其能通过对椭圆轨道对高轨高轨轨道碎片的相对运动情况分析,计算相对运动的周期、角速度、速度、距离变化情况,分析构建一种双星轨道来模拟这种交会运动,以较高性价比验证椭圆轨道对高轨轨道的相对运动情况,遵行先易后难的原则,双星运行在低轨轨道,创新性的利用近地点幅角相差半周的椭圆轨道,调节双椭圆不同的远地点高度获得不同的相对运动速度,调节双星在两轨道交会处的距离获得不同的相对角速度,通过这样一种模拟未来零倾角椭圆轨道对高轨轨道的交会工况,不论从相对运动的速度和角速度都满足模拟的要求,相当于以较低代价在低轨验证了高轨的交会工况,该方法具有重要的应用价值。
本发明是通过下述技术方案来解决上述技术问题的:一种相对运动轨道构型的模拟方法,其包括如下步骤:
步骤一,设计低轨轨道的近地点高度和远地点高度,满足不同的交会速度;
步骤二,设计合理的近地点幅角,使得双星椭圆轨道交会点在我国上空,满足测控弧段要求;
步骤三,设计两星距离,使得第一星在相对运动椭圆半长轴的合适位置,可调整交会时刻的两星距离;
步骤四,拟定试验步骤,从变轨、跟踪、指向、试验的先后顺序,分解试验内容;
步骤五,在STK里进行相对运动的仿真,验证设计正确性。
优选地,所述步骤一中计算交会点的相对运动速度dv=vLE0*de,因此对应100m/s速度差时,近地点幅角相差180°的双星,偏心率应为100/7635/2=0.007。
优选地,所述步骤二中计算交会点的纬度幅角与第一星或者第二星近地点幅角的对应关系。
优选地,所述步骤三中调整交会时刻的相对距离关系,设置第一星在相对运动椭圆半长轴上的不同位置来实现。
优选地,所述步骤五通过建立STK场景,输入卫星轨道参数,拉出相对运动构型图、数据及曲线。
本发明的积极进步效果在于:本发明通过对椭圆轨道对高轨高轨轨道碎片的相对运动情况分析,计算相对运动的周期、角速度、速度、距离变化情况,分析构建一种双星轨道来模拟这种交会运动,以较高性价比验证椭圆轨道对高轨轨道的相对运动情况,遵行先易后难的原则,双星运行在低轨轨道,创新性的利用近地点幅角相差半周的椭圆轨道,调节双椭圆不同的远地点高度获得不同的相对运动速度,调节双星在两轨道交会处的距离获得不同的相对角速度,通过这样一种模拟未来零倾角椭圆轨道对高轨轨道的交会工况,不论从相对运动的速度和角速度都满足模拟的要求,相当于以较低代价在低轨验证了高轨的交会工况。
附图说明
图1为本发明的椭圆轨道对高轨轨道交会实际情况示意图。
图2为本发明的低轨双星椭圆模拟交会示意图。
图3为本发明对高轨轨道接近过程相对距离随时间变化图。
图4为本发明对高轨轨道接近过程视线转率随时间变化图。
图5为本发明对高轨轨道接近过程视线转率变化率随时间变化图。
具体实施方式
下面结合附图给出本发明较佳实施例,以详细说明本发明的技术方案。
如图1至图5所示,本发明相对运动轨道构型的模拟方法包括如下步骤:
步骤一,设计低轨轨道的近地点高度和远地点高度,满足不同的交会速度,如100m/s(米/秒),200m/s(米/秒),及300m/s(米/秒);
步骤二,设计合理的近地点幅角,使得双星椭圆轨道B交会点A在我国上空,满足测控弧段要求;
步骤三,设计两星距离,使得第一星(D星)在相对运动椭圆半长轴的合适位置,可调整交会时刻的两星距离;H代表第一星轨道,G代表第二星轨道。
步骤四,拟定试验步骤,从变轨、跟踪、指向、试验的先后顺序,分解试验内容;
步骤五,在STK(卫星工具包)里进行相对运动的仿真,验证设计正确性。
所述步骤一中计算交会点的相对运动速度dv=vLE0*de(dv:相对运动速度Vleo:低轨轨道运动速度de:双星偏心率矢量差),因此对应100m/s速度差时,近地点幅角相差180°的双星,偏心率应为100/7635/2=0.007。
所述步骤二中计算交会点的纬度幅角与第一星(D星)或者第二星(E星)近地点幅角的对应关系。
所述步骤三中调整交会时刻的相对距离关系,设置第一星在相对运动椭圆半长轴上的不同位置来实现。
所述在STK里进行相对运动的仿真,可通过建立STK场景,输入卫星轨道参数,拉出相对运动构型图、数据及曲线。
本发明的工作原理如下:建立零倾角椭圆轨道对高轨轨道C目标的访问场景,椭圆轨道近地点高度约4500公里;取得椭圆轨道对高轨轨道的相对运动关系;以零倾角椭圆轨道访问高轨轨道目标时,轨道面内切向速度差1090m/s(米/秒),轨道面外速度差取决于倾角差,20°(度)倾角差时,相对速度1000m/s(米/秒),因此,合成速度范围在1-1.5km/s(千米/秒),在10公里处,角速度6-8°/s(度/秒);设计低轨双椭圆交会轨道的轨道高度和卫星相位;取得低轨双椭圆轨道的相对运动数据,与高轨交会运动参数进行对照,验证模拟效果,结果表明发明的低轨双椭圆模拟交会轨道运动参数可等效模拟高轨交会工况;分析试验步骤;
以形成相对运动速度300m/s为例,试验步骤如下:
一、在轨道相位差180°的对称约束下,通过星地联合控制,先后将两星远地点高度抬升,形成相对运动速度300m/s的轨道面内交会关系,第一星在前,轨道飞行时先通过交会过捷点;
二、在上述“双椭圆”300m/s相对速度轨道面内交会条件下,下一次轨道交会点前,调整两星相位与轨道高度,使得两星交会时相距10km;
三、运控评估系统完成两星状态检查,进行试验规划,生成控制数据;
四、空间碎片感知与清除载荷在作用距离满足时适时工作。
本发明的实施例采用如表1所示的不同交会速度下的双星轨道参数设计表:
表1不同交会速度下的双星轨道参数设计表
综上所述,本发明通过对椭圆轨道对高轨高轨轨道碎片的相对运动情况分析,计算相对运动的周期、角速度、速度、距离变化情况,分析构建一种双星轨道来模拟这种交会运动,以较高性价比验证椭圆轨道对高轨轨道的相对运动情况,遵行先易后难的原则,双星运行在低轨轨道,创新性的利用近地点幅角相差半周的椭圆轨道,调节双椭圆不同的远地点高度获得不同的相对运动速度,调节双星在两轨道交会处的距离获得不同的相对角速度,通过这样一种模拟未来零倾角椭圆轨道对高轨轨道的交会工况,不论从相对运动的速度和角速度都满足模拟的要求,相当于以较低代价在低轨验证了高轨的交会工况,该方法具有重要的应用价值。
以上所述的具体实施例,对本发明的解决的技术问题、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种相对运动轨道构型的模拟方法,其特征在于,其包括如下步骤:
步骤一,设计低轨轨道的近地点高度和远地点高度,满足不同的交会速度;
步骤二,设计合理的近地点幅角,使得双星椭圆轨道交会点在我国上空,满足测控弧段要求;
步骤三,设计两星距离,使得第一星在相对运动椭圆半长轴的合适位置,可调整交会时刻的两星距离;
步骤四,拟定试验步骤,从变轨、跟踪、指向、试验的先后顺序,分解试验内容;
步骤五,在STK里进行相对运动的仿真,验证设计正确性。
2.如权利要求1所述的相对运动轨道构型的模拟方法,其特征在于,所述步骤一中计算交会点的相对运动速度dv=vLEO*de,因此对应100m/s速度差时,近地点幅角相差180°的双星,偏心率应为100/7635/2=0.007。
3.如权利要求1所述的相对运动轨道构型的模拟方法,其特征在于,所述步骤二中计算交会点的纬度幅角与第一星或者第二星近地点幅角的对应关系。
4.如权利要求1所述的相对运动轨道构型的模拟方法,其特征在于,所述步骤三中调整交会时刻的相对距离关系,设置第一星在相对运动椭圆半长轴上的不同位置来实现。
5.如权利要求1所述的相对运动轨道构型的模拟方法,其特征在于,所述步骤五通过建立STK场景,输入卫星轨道参数,拉出相对运动构型图、数据及曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710868984.9A CN107633142B (zh) | 2017-09-22 | 2017-09-22 | 相对运动轨道构型的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710868984.9A CN107633142B (zh) | 2017-09-22 | 2017-09-22 | 相对运动轨道构型的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107633142A true CN107633142A (zh) | 2018-01-26 |
CN107633142B CN107633142B (zh) | 2021-03-16 |
Family
ID=61103583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710868984.9A Active CN107633142B (zh) | 2017-09-22 | 2017-09-22 | 相对运动轨道构型的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107633142B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110471432A (zh) * | 2019-07-04 | 2019-11-19 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354218A (zh) * | 2011-06-24 | 2012-02-15 | 哈尔滨工业大学 | 一种航天器相对运动的采样控制方法 |
CN102841966A (zh) * | 2012-08-28 | 2012-12-26 | 曾安里 | Vpp-STK卫星仿真开发与运行平台系统 |
CN103699129A (zh) * | 2013-12-06 | 2014-04-02 | 上海卫星工程研究所 | 基于低轨观测卫星的地球同步轨道目标快速遍历方法 |
CN104614987A (zh) * | 2014-12-19 | 2015-05-13 | 中国空间技术研究院 | 一种用于高轨卫星锥杆式抓捕机构的实时控制系统及方法 |
CN106325099A (zh) * | 2016-07-28 | 2017-01-11 | 中国科学院国家空间科学中心 | 一种基于伪相对运动的航天器实时轨道改进方法 |
-
2017
- 2017-09-22 CN CN201710868984.9A patent/CN107633142B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354218A (zh) * | 2011-06-24 | 2012-02-15 | 哈尔滨工业大学 | 一种航天器相对运动的采样控制方法 |
CN102841966A (zh) * | 2012-08-28 | 2012-12-26 | 曾安里 | Vpp-STK卫星仿真开发与运行平台系统 |
CN103699129A (zh) * | 2013-12-06 | 2014-04-02 | 上海卫星工程研究所 | 基于低轨观测卫星的地球同步轨道目标快速遍历方法 |
CN104614987A (zh) * | 2014-12-19 | 2015-05-13 | 中国空间技术研究院 | 一种用于高轨卫星锥杆式抓捕机构的实时控制系统及方法 |
CN106325099A (zh) * | 2016-07-28 | 2017-01-11 | 中国科学院国家空间科学中心 | 一种基于伪相对运动的航天器实时轨道改进方法 |
Non-Patent Citations (1)
Title |
---|
吴昊 等: "基于STK的小卫星轨道交会设计研究", 《飞行器测控学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110471432A (zh) * | 2019-07-04 | 2019-11-19 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
CN110471432B (zh) * | 2019-07-04 | 2020-09-08 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107633142B (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105792135B (zh) | 一种定位车辆所在车道的方法及装置 | |
WO2020192149A1 (zh) | 轨迹跟踪控制器的测试方法、装置、介质及设备 | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN105404296A (zh) | 一种两轮自平衡智能车姿态控制方法 | |
CN108152529A (zh) | 一种基于飞行参数计算风速及风向的方法 | |
Bolívar-Vincenty et al. | Modelling the ball-and-beam system from newtonian mechanics and from lagrange methods | |
US20120212622A1 (en) | Moving object image tracking apparatus and method | |
Zhang et al. | Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions | |
CN103900576A (zh) | 一种深空探测自主导航的信息融合方法 | |
CN108917764A (zh) | 一种双星编队仅测距相对导航方法 | |
CN108981750A (zh) | X射线脉冲双星光子序列仿真方法 | |
CN102393200A (zh) | 基于飞行仿真的通用惯导测试方法 | |
CN102620748A (zh) | 捷联惯导系统晃动基座条件下杆臂效应的估计和补偿方法 | |
CN103955223A (zh) | 一种深空探测软着陆过程的姿轨耦合控制方法 | |
CN102538819A (zh) | 基于双圆锥红外和星敏感器的自主导航半物理仿真试验系统 | |
CN102915036A (zh) | 一种参数不确定性飞行器倾斜角控制系统极限环抑制方法 | |
Evans et al. | DragonFly: a versatile UAV platform for the advancement of aircraft navigation and control | |
CN103983954A (zh) | 用于雷达跟瞄高精度地面测试的误差补偿系统及方法 | |
CN103712623A (zh) | 基于角速率输入的光纤陀螺惯导系统姿态优化方法 | |
CN108957509A (zh) | 一种双星编队周期相对运动仅测距相对导航解析方法 | |
CN106672265A (zh) | 一种基于光流信息的小天体定点着陆制导控制方法 | |
CN107633142A (zh) | 相对运动轨道构型的模拟方法 | |
CN107024702A (zh) | 一种基于惯性/卫星导航的半实物控制系统闭环仿真平台 | |
Liu et al. | Mass and mass center identification of target satellite after rendezvous and docking | |
Briales et al. | Track frame approach for heading and attitude estimation in operating railways using on-board MEMS sensor and encoder |
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 |