主权项 |
一种电动绳系航天器渐近稳定释放控制方法,其特征在于包括以下步骤:步骤A、构建描述电动绳索释放过程中系统俯仰角、滚转角摆动的无量纲动力学方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mover><mi>θ</mi><mo>··</mo></mover><mo>+</mo><mn>2</mn><mrow><mo>(</mo><mover><mi>θ</mi><mo>·</mo></mover><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mover><mi>ξ</mi><mo>·</mo></mover><mi>ξ</mi></mfrac><mo>-</mo><mi>φ</mi><mover><mi>φ</mi><mo>·</mo></mover><mo>)</mo></mrow><mo>+</mo><mn>3</mn><mi>θ</mi><mo>=</mo><msub><mi>H</mi><mi>θ</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>φ</mi><mo>··</mo></mover><mo>+</mo><mn>2</mn><mover><mi>φ</mi><mo>·</mo></mover><mrow><mo>(</mo><mfrac><mi>ξ</mi><mi>ξ</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mo>[</mo><msup><mrow><mo>(</mo><mover><mi>θ</mi><mo>·</mo></mover><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>3</mn><mo>]</mo><mi>φ</mi><mo>=</mo><msub><mi>H</mi><mi>φ</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000783497770000011.GIF" wi="1150" he="350" /></maths>式中,θ表示电动绳索摆动的面内俯仰角,<img file="FDA0000783497770000012.GIF" wi="49" he="70" />表示面内俯仰角对轨道真近点角ν求导,即绳索摆动的面内俯仰角速度,<img file="FDA0000783497770000013.GIF" wi="46" he="71" />表示绳索摆动的面内俯仰角的二阶导数,φ表示电动绳索摆动的面外滚转角,<img file="FDA0000783497770000014.GIF" wi="45" he="86" />表示面外滚转角对轨道真近点角ν求导,即绳索摆动的面外滚转角速度,<img file="FDA0000783497770000015.GIF" wi="56" he="87" />表示绳索摆动的面外滚转角的二阶导数,ξ表示绳索无量纲长度,<img file="FDA0000783497770000016.GIF" wi="46" he="86" />表示绳索无量纲长度的变化速率,参数H<sub>θ</sub>=‑[μ<sub>m</sub>I(m<sub>M</sub>‑m<sub>S</sub>)cosδ]/[2μ<sub>E</sub>m<sub>M</sub>m<sub>S</sub>]、H<sub>φ</sub>=0,其中μ<sub>m</sub>表示地球偶极子磁矩,I表示绳索电流强度,m<sub>M</sub>表示航天器质量,m<sub>S</sub>表示末端探测器质量,δ表示轨道倾角,μ<sub>E</sub>表示地球引力常数,若电流I恒定则参数H<sub>θ</sub>也为恒定值;步骤B、根据公式(1)确定系统的平衡点:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>θ</mi><mn>10</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>ξH</mi><mi>θ</mi></msub><mo>-</mo><mn>2</mn><mover><mi>ξ</mi><mo>·</mo></mover></mrow><mrow><mn>3</mn><mi>ξ</mi></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>θ</mi><mn>20</mn></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>φ</mi><mn>10</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>H</mi><mi>φ</mi></msub><mn>4</mn></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>φ</mi><mn>20</mn></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000783497770000017.GIF" wi="917" he="519" /></maths>式中,θ<sub>10</sub>为平衡位置俯仰角,θ<sub>20</sub>为对应的俯仰角速度,φ<sub>10</sub>为平衡位置滚转角,φ<sub>20</sub>为对应的滚转角速度,由于H<sub>φ</sub>=0故平衡位置滚转角φ<sub>10</sub>=0;步骤C、任意设置一个期望俯仰角θ<sub>e</sub>,并令平衡位置俯仰角θ<sub>10</sub>=θ<sub>e</sub>,确定θ<sub>e</sub>与参数H<sub>θ</sub>、绳索无量纲长度ξ及其变化速率<img file="FDA0000783497770000021.GIF" wi="44" he="87" />的关系,如下<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>ξ</mi><mo>·</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>ξH</mi><mi>θ</mi></msub><mo>-</mo><mn>3</mn><msub><mi>ξθ</mi><mi>e</mi></msub></mrow><mn>2</mn></mfrac><mo>;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000783497770000022.GIF" wi="948" he="143" /></maths>步骤D、根据系统线性化矩阵的特征根及释放特性来确定期望俯仰角θ<sub>e</sub>取值范围,使得在对电动绳索释放时,保证绳索释放过程渐近稳定;步骤E、在期望俯仰角θ<sub>e</sub>满足取值范围的条件下,根据式(3)绳索释放的无量纲释放速率对电动绳索进行释放控制。 |