发明名称 一种飞行器非线性姿态控制系统稳定性分析方法
摘要 本发明公开了一种飞行器非线性姿态控制系统稳定性分析方法,属于飞行器姿态控制领域。该方法首先计算飞行器飞行状态对应的动力系数,用死区阈值h表达其中的对应项;之后绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断非线性系统稳定,设计俯仰通道死区阈值为h的取值,否则判断非线性系统不稳定;该方法给出了不同阈值条件下闭环系统的稳定性判别方法,从理论上提出了死区阈值对应的闭环回路稳定性判定依据。
申请公布号 CN106383523A 申请公布日期 2017.02.08
申请号 CN201610874908.4 申请日期 2016.09.30
申请人 湖北航天技术研究院总体设计所 发明人 鲍晓强;朱伟;陈兴褔;陈红艳;杨登峰;田群方;蒋金龙;郭一江;李焰;崔跃军;叶昌
分类号 G05D1/08(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 方放
主权项 一种飞行器非线性姿态控制系统稳定性分析方法,其特征在于,该方法包括以下步骤:(1)采集飞行器飞行状态下的动压q、参考面积s、参考长度l、质量m、马赫数M<sub>a</sub>、速度V、侧滑角β、攻角α、弹道倾角θ;(2)计算飞行器飞行状态对应的俯仰通道动力系数b<sub>1</sub>、b<sub>2</sub>、c<sub>1</sub>和c<sub>2</sub>:<maths num="0001"><math><![CDATA[<mrow><msub><mi>b</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msubsup><mi>m</mi><mi>z</mi><mover><mi>&omega;</mi><mo>&OverBar;</mo></mover></msubsup><msup><mi>qsl</mi><mn>2</mn></msup></mrow><mrow><msub><mi>J</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><mi>V</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001125581780000011.GIF" wi="299" he="135" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>b</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mn>57.3</mn><msubsup><mi>m</mi><mi>z</mi><mi>&alpha;</mi></msubsup><mi>q</mi><mi>s</mi><mi>l</mi></mrow><msub><mi>J</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub></mfrac></mrow>]]></math><img file="FDA0001125581780000012.GIF" wi="361" he="127" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>57.3</mn><msubsup><mi>C</mi><mi>L</mi><mi>&alpha;</mi></msubsup><mi>qs</mi></mrow><mi>mV</mi></mfrac></mrow>]]></math><img file="FDA0001125581780000019.GIF" wi="299" he="118" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mi>g</mi><mi>sin</mi><mi>&theta;</mi></mrow><mi>V</mi></mfrac></mrow>]]></math><img file="FDA0001125581780000013.GIF" wi="246" he="117" /></maths>式中,J<sub>z1</sub>为弹体绕OZ1轴的转动惯量;<img file="FDA0001125581780000014.GIF" wi="77" he="71" />为俯仰通道阻尼力矩系数;<img file="FDA0001125581780000015.GIF" wi="62" he="70" />为俯仰力矩系数对攻角导数;<img file="FDA0001125581780000016.GIF" wi="62" he="63" />为升力系数对攻角导数;g为重力加速度;b<sub>3</sub>和c<sub>3</sub>可表示为含有死区阈值h的表达式:<maths num="0005"><math><![CDATA[<mrow><msub><mi>b</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>F</mi><mrow><mi>c</mi><mi>t</mi><mi>r</mi><mo>_</mo><mi>p</mi></mrow></msub><msub><mi>L</mi><mrow><mi>c</mi><mi>t</mi><mi>r</mi><mo>_</mo><mi>p</mi></mrow></msub></mrow><mrow><msub><mi>J</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><mi>h</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001125581780000017.GIF" wi="331" he="135" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>c</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>F</mi><mrow><mi>c</mi><mi>t</mi><mi>r</mi><mo>_</mo><mi>p</mi></mrow></msub><mi>c</mi><mi>o</mi><mi>s</mi><mi>&alpha;</mi></mrow><mrow><mi>m</mi><mi>V</mi><mi>h</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001125581780000018.GIF" wi="340" he="142" /></maths>式中,F<sub>ctr_p</sub>为俯仰通道单个姿控喷管推力大小;L<sub>ctr_p</sub>为俯仰通道单个姿控喷管推力作用力臂;(3)预选一个死区阈值h,代入到动力系数中,绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断对应非线性系统稳定,设计俯仰通道死区阈值为h的预选值;否则判断对应非线性系统不稳定。
地址 430040 湖北省武汉市金山大道9号