发明名称 一种基于零射程线的运载火箭耗尽关机控制方法
摘要 本发明提供一种基于零射程线的运载火箭耗尽关机控制方法,涉及制导控制领域。本发明提出一种闭路制导与零射程线横、法向导引相结合的能量管理及射程精确控制总体方案,具体的包括基于起算点到轨道终端脉冲转移函数的高精度闭路制导和基于零射程线的多余能量耗散控制及快速收敛的横、法向导引控制。本发明通过将姿态调整到零射程线耗散多余的能量,对速度模量的控制精度没有严格要求,还可以利用发动机的能量,继续在零射程线附近进行小量的姿态实时调整,不需要增加额外的硬件,动作简单快速、精度高,可显著简化姿态调制程序,解决了需要速度偏差对射程控制精度影响大的难题以及大偏差条件下导引控制的快速收敛问题。
申请公布号 CN105501467B 申请公布日期 2017.03.08
申请号 CN201510922157.4 申请日期 2015.12.11
申请人 周沁心 发明人 周沁心
分类号 B64G1/26(2006.01)I 主分类号 B64G1/26(2006.01)I
代理机构 武汉楚天专利事务所 42113 代理人 孔敏
主权项 一种基于零射程线的运载火箭耗尽关机控制方法,其特征在于包括如下步骤:步骤一、在第三级点火时,根据当前程序角<img file="FDA0001126001920000011.GIF" wi="82" he="46" />零射程线对应的程序角<img file="FDA0001126001920000012.GIF" wi="51" he="47" />和调姿时间长度Δt<sub>t</sub>,给出调姿程序角<img file="FDA0001126001920000013.GIF" wi="112" he="59" /><img file="FDA0001126001920000014.GIF" wi="438" he="118" />其中t为以调姿起点起算的时间;步骤二、在随后的第三级发动机工作过程中实时进行速度待增量<img file="FDA0001126001920000015.GIF" wi="510" he="77" />及相应发动机工作时间Δt的闭路制导计算:计算过程包括如下三部分:a.以当前轨道参数,即当前轨道的时间T<sub>0</sub>、位置<img file="FDA0001126001920000016.GIF" wi="82" he="70" />速度<img file="FDA0001126001920000017.GIF" wi="51" he="68" />为起点,通过计算得到终端轨道的时间T<sub>e</sub>、位置<img file="FDA0001126001920000018.GIF" wi="83" he="67" />速度<img file="FDA0001126001920000019.GIF" wi="67" he="70" />其中T<sub>e</sub>取轨道y向位置分量等于虚拟目标点的y向位置分量所对应的时间;计算中主要考虑发动机推力<img file="FDA00011260019200000110.GIF" wi="61" he="87" />和引力场<img file="FDA00011260019200000111.GIF" wi="50" he="55" />的影响,即:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><mi>d</mi><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mi>d</mi><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mover><mover><mi>W</mi><mo>&CenterDot;</mo></mover><mo>&RightArrow;</mo></mover><mi>p</mi></msub><mo>+</mo><msub><mover><mi>g</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00011260019200000112.GIF" wi="323" he="262" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>I</mi><mn>0</mn></mrow></msub></mrow>]]></math><img file="FDA00011260019200000113.GIF" wi="230" he="71" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mrow><mi>I</mi><mn>0</mn></mrow></msub></mrow>]]></math><img file="FDA00011260019200000114.GIF" wi="237" he="67" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>e</mi></msub><mo>=</mo><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>I</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>e</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011260019200000115.GIF" wi="238" he="71" /></maths><img file="FDA00011260019200000116.GIF" wi="228" he="71" />其中<img file="FDA00011260019200000117.GIF" wi="147" he="70" />分别为轨道的位置、速度;<img file="FDA00011260019200000118.GIF" wi="276" he="71" />分别为轨道的位置初值、速度初值,其对应的数值分别为<img file="FDA00011260019200000119.GIF" wi="68" he="71" />和<img file="FDA00011260019200000120.GIF" wi="84" he="87" />b.计算T<sub>e</sub>时的虚拟目标点的位置及速度<img file="FDA0001126001920000021.GIF" wi="567" he="167" />其中:T<sub>e0</sub>、<img file="FDA0001126001920000022.GIF" wi="220" he="70" />分别为标准条件下,飞行到虚拟目标点的时间和相应的虚拟目标点的位置及速度;<img file="FDA0001126001920000023.GIF" wi="190" he="83" />分别为虚拟目标点位置及速度的变化率;<img file="FDA0001126001920000024.GIF" wi="186" he="71" />为T<sub>e</sub>时的虚拟目标点的位置及速度;c.计算T<sub>e</sub>时的位置<img file="FDA0001126001920000025.GIF" wi="78" he="70" />及速度偏差<img file="FDA0001126001920000026.GIF" wi="99" he="55" /><maths num="0005"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&Delta;</mi><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>e</mi></msub><mo>=</mo><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mi>e</mi></msub><mo>-</mo><msub><mover><mi>R</mi><mo>&RightArrow;</mo></mover><mrow><mi>m</mi><mi>e</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>&Delta;</mi><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>e</mi></msub><mo>=</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>e</mi></msub><mo>-</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>m</mi><mi>e</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001126001920000027.GIF" wi="332" he="151" /></maths>对导航计算公式进行线性化处理,得到起算点到终端的脉冲转移函数,进而给出需要速度待增量<img file="FDA0001126001920000028.GIF" wi="513" he="77" />及相应发动机工作时间Δt的计算公式:<img file="FDA0001126001920000029.GIF" wi="894" he="447" />其中:ΔX<sub>e</sub>,ΔY<sub>e</sub>,ΔZ<sub>e</sub>为T<sub>e</sub>时的位置偏差<img file="FDA00011260019200000210.GIF" wi="83" he="70" />的分量,ΔV<sub>xe</sub>,ΔV<sub>ye</sub>,ΔV<sub>ze</sub>为T<sub>e</sub>时的速度偏差<img file="FDA00011260019200000211.GIF" wi="82" he="71" />的分量,ΔV<sub>x,i‑1</sub>、ΔV<sub>y,i‑1</sub>、ΔV<sub>z,i‑1</sub>为速度待增量的前点值,ΔV<sub>x,i</sub>、ΔV<sub>y,i</sub>、ΔV<sub>z,i</sub>为当前速度待增量,ΔV为速度待增量<img file="FDA00011260019200000212.GIF" wi="74" he="54" />的模,T<sub>0</sub>为当前时间,按照制导周期,循环上述算法,需要速度待增量及相应发动机工作时间会快速收敛;步骤三、当发动机工作时间Δt=0时,对应的当前时刻为t<sub>0</sub>,按照步骤一中给出的调姿程序角<img file="FDA0001126001920000031.GIF" wi="91" he="55" />公式开始俯仰调姿;步骤四、当俯仰调姿角<img file="FDA0001126001920000032.GIF" wi="187" he="55" />时,开始沿零射程线的横、法向导引控制,横、法向导引控制量根据步骤二得到的速度待增量ΔV<sub>x</sub>、ΔV<sub>y</sub>、ΔV<sub>z</sub>按照如下公式计算得出:<maths num="0006"><math><![CDATA[<mrow><msub><mi>U</mi><mi>z</mi></msub><mo>=</mo><msub><mi>K</mi><mi>z</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;V</mi><mi>z</mi></msub><mo>/</mo><msub><mover><mi>W</mi><mo>&CenterDot;</mo></mover><mi>p</mi></msub></mrow>]]></math><img file="FDA0001126001920000033.GIF" wi="374" he="71" /></maths><img file="FDA0001126001920000034.GIF" wi="766" he="71" />其中U<sub>z</sub>、U<sub>y</sub>为横、法向导引控制量,<img file="FDA0001126001920000035.GIF" wi="60" he="70" />为发动机推力产生的过载,<img file="FDA0001126001920000036.GIF" wi="34" he="47" />为当前俯仰姿态角,K<sub>z</sub>、K<sub>y</sub>为横、法向导引系数;步骤五、当发动机推力耗尽后,全部耗尽关机控制过程结束。
地址 430034 湖北省武汉市解放大道347号武汉市第四中学
您可能感兴趣的专利