发明名称 一种航天器编队维持或绕飞撤离的控制方法
摘要 本发明公开了一种航天器编队维持或绕飞撤离的控制方法,在目标航天器相邻的两个轨道周期中,分别对主控航天器相对于目标航天器轨道坐标系x坐标分量进行积分,对积分结果关于该周期的平方做差商,得到周期平均漂移速度;根据需要设定编队保持控制完成之后或绕飞撤离之后期望的周期平均漂移速度,该期望漂移速度减去控前周期平均漂移速度得到周期平均漂移速度增量,编队保持控制的双脉冲及绕飞撤离控制单脉冲的水平分量均由该增量计算得到,从而使编队保持控制具有自主及低频度等特征,使绕飞撤离自主且安全快捷。
申请公布号 CN104317303A 申请公布日期 2015.01.28
申请号 CN201410532652.X 申请日期 2014.10.10
申请人 北京控制工程研究所 发明人 苟兴宇;韩冬;李克行;张斌;董筠;赵键;张洪华;刘其睿;曾春平;谢晓兵;李鹤;张欣;徐子荔;刘一薇;赵志明
分类号 G05D1/10(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 褚鹏蛟
主权项 一种航天器编队维持控制方法,其特征在于步骤如下:(1)在预期后续一段时间内主控航天器与目标航天器均不会发生轨道控制的情况下,指定一个进行周期平均漂移速度计算的采样初始时刻t<sub>00</sub>;(2)按以下公式进行周期平均漂移速度<img file="FDA0000584414990000011.GIF" wi="140" he="81" />的计算:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>drift</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mo>&Integral;</mo><mrow><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mi>T</mi></mrow><mrow><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mn>2</mn><mi>T</mi></mrow></msubsup><msub><mi>x</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi><mo>-</mo><msubsup><mo>&Integral;</mo><msub><mi>t</mi><mn>0</mn></msub><mrow><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mi>T</mi></mrow></msubsup><msub><mi>x</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi></mrow><msup><mi>T</mi><mn>2</mn></msup></mfrac></mrow>]]></math><img file="FDA0000584414990000012.GIF" wi="750" he="174" /></maths>式中x<sub>r</sub>(t)为采样时刻t对应的主控航天器在目标航天器的轨道坐标系中的x轴位置分量,根据测量得到;T为目标航天器轨道周期;t<sub>0</sub>为相邻两个轨道周期中前一个轨道周期的积分初始时刻,t<sub>0</sub>取值为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>t</mi><mn>0</mn></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>t</mi><mn>00</mn></msub></mtd><mtd><mi>t</mi><mo>&le;</mo><msub><mi>t</mi><mn>00</mn></msub><mo>+</mo><mn>2</mn><mi>T</mi></mtd></mtr><mtr><mtd><mi>t</mi><mo>-</mo><mn>2</mn><mi>T</mi></mtd><mtd><mi>t</mi><mo>></mo><msub><mi>t</mi><mn>00</mn></msub><mo>+</mo><mn>2</mn><mi>T</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000584414990000013.GIF" wi="536" he="157" /></maths>(3)判断时刻t是否满足t&gt;t<sub>00</sub>+2T,当满足t&gt;t<sub>00</sub>+2T时,获得第一个收敛可用的周期平均漂移速度结果,并转入步骤(4);否则,返回步骤(2);(4)判断时刻t是否到达编队维持控制启动时刻t<sub>beforeKeep</sub>,如果到达时刻t<sub>beforeKeep</sub>,获得该时刻t<sub>beforeKeep</sub>对应的周期平均漂移速度<img file="FDA0000584414990000014.GIF" wi="314" he="94" />转入步骤(5);如果没有到达时刻t<sub>beforeKeep</sub>,则返回步骤(2);(5)根据周期平均漂移速度<img file="FDA0000584414990000015.GIF" wi="296" he="97" />和期望的维持控制完成后的周期平均漂移速度<img file="FDA0000584414990000016.GIF" wi="274" he="100" />计算周期平均漂移速度增量<img file="FDA0000584414990000017.GIF" wi="221" he="87" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Delta;</mi><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mrow><mi>drift</mi><mo>,</mo><mi>Keep</mi></mrow></msub><mo>=</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>drift</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>afterKeep</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>drift</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>beforeKeep</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000584414990000018.GIF" wi="847" he="99" /></maths>(6)根据周期平均漂移速度增量<img file="FDA0000584414990000019.GIF" wi="189" he="85" />计算编队维持控制水平双脉冲中实现相对长半轴长调整的构成分量Δv<sub>a</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&Delta;v</mi><mi>a</mi></msub><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mi>&Delta;</mi><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mrow><mi>drift</mi><mo>,</mo><mi>Keep</mi></mrow></msub></mrow></mrow>]]></math><img file="FDA00005844149900000110.GIF" wi="392" he="131" /></maths>(7)根据编队维持控制实施之前的相对偏心率矢量<img file="FDA00005844149900000111.GIF" wi="70" he="60" />和期望的维持控制完成后的目标相对偏心率矢量<img file="FDA0000584414990000021.GIF" wi="74" he="72" />计算相对偏心率矢量增量<img file="FDA0000584414990000022.GIF" wi="127" he="60" />其中<img file="FDA0000584414990000023.GIF" wi="353" he="76" />并确定相对偏心率矢量增量<img file="FDA0000584414990000024.GIF" wi="92" he="61" />在目标航天器节点坐标系Ox<sub>n</sub>y<sub>n</sub>z<sub>n</sub>中x<sub>n</sub>轴与y<sub>n</sub>轴的坐标分量δΔe<sub>x</sub>与δΔe<sub>y</sub>;(8)根据公式<img file="FDA0000584414990000025.GIF" wi="517" he="103" />计算得到相对偏心率矢量增量的幅值δΔe;根据相对偏心率矢量增量的幅值δΔe计算出编队维持控制水平双脉冲中实现相对偏心率调整的构成分量Δv<sub>e</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&Delta;v</mi><mi>e</mi></msub><mo>=</mo><mfrac><mrow><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mrow><mi>t</mi><mi>arg</mi><mi>et</mi></mrow></msub><msub><mi>n</mi><mrow><mi>t</mi><mi>arg</mi><mi>et</mi></mrow></msub></mrow><mn>4</mn></mfrac><mi>&delta;&Delta;e</mi></mrow>]]></math><img file="FDA0000584414990000026.GIF" wi="422" he="131" /></maths>式中<img file="FDA0000584414990000027.GIF" wi="115" he="75" />为目标航天器的平长半轴长,n<sub>target</sub>为目标航天器的轨道角速率;(9)根据Δv<sub>a</sub>和Δv<sub>e</sub>确定编队维持控制的水平双脉冲大小,计算公式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&Delta;</mi><msub><mi>v</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&Delta;v</mi><mi>a</mi></msub><mo>+</mo><msub><mi>&Delta;v</mi><mi>e</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000584414990000028.GIF" wi="402" he="129" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Delta;v</mi><mrow><mi>x</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&Delta;v</mi><mi>a</mi></msub><mo>-</mo><msub><mi>&Delta;v</mi><mi>e</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000584414990000029.GIF" wi="406" he="130" /></maths>其中,Δv<sub>x1</sub>为第一脉冲的大小,Δv<sub>x2</sub>为第二脉冲的大小;(10)根据δΔe<sub>x</sub>、δΔe<sub>y</sub>确定水平双脉冲在主控航天器上作用点的轨道幅角:u<sub>x1</sub>=atan2(δΔe<sub>y</sub>,δΔe<sub>x</sub>)u<sub>x2</sub>=atan2(‑δΔe<sub>y</sub>,‑δΔe<sub>x</sub>)其中,u<sub>x1</sub>为所述第一脉冲在主控航天器上作用点的轨道幅角,u<sub>x2</sub>为所述第二脉冲在主控航天器上作用点的轨道幅角,atan2(*,*)为工程适用的拓展值域到[‑π,π]范围的反正切函数;(11)在轨道幅角u<sub>x1</sub>处对主控航天器施加第一脉冲Δv<sub>x1</sub>,在轨道幅角u<sub>x2</sub>处对主控航天器施加第二脉冲Δv<sub>x2</sub>,完成编队维持控制。
地址 100080 北京市2729信箱