发明名称 用于自旋恢复的直升机非线性预测控制方法及装置
摘要 本发明公开了一种用于自旋恢复的直升机非线性预测控制方法。在进入自旋之后,利用预先训练好的直升机需求扭矩模型对当前的直升机需求扭矩进行实时在线预测;在进入自旋恢复阶段后,利用预先训练好的发动机动态参数模型对当前的发动机动态参数进行实时在线预测,同时利用直升机需求扭矩模型和发动机动态参数模型的在线预测结果,求解以缩小离合器接通时的直升机需求扭矩和发动机提供的扭矩支持间的差值为目标,并考虑发动机稳定安全的运行条件的滚动优化问题,并将求解得到的控制量序列的第一项作为当前输入的直升机控制量。本发明还公开了直升机非线性预测控制装置。本发明可有效地缩短自旋恢复阶段的时间延迟并减小旋翼转速瞬态下滑。
申请公布号 CN105867121A 申请公布日期 2016.08.17
申请号 CN201610184036.9 申请日期 2016.03.28
申请人 南京航空航天大学 发明人 郑前钢;徐田镇;刘明磊;杜瑶;张海波
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 杨楠
主权项 一种用于自旋恢复的直升机非线性预测控制方法,其特征在于,在进入自旋之后,利用预先训练好的直升机需求扭矩模型对当前的直升机需求扭矩进行实时在线预测;在进入自旋恢复阶段后,利用预先训练好的发动机动态参数模型对当前的发动机动态参数进行实时在线预测,与此同时,利用直升机需求扭矩模型和发动机动态参数模型的在线预测结果,求解以下滚动优化问题,并将求解得到的控制量序列的第一项作为当前输入的直升机控制量:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>x</mi></munder><mi>J</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>&delta;</mi><mi>p</mi></msub></munderover><msub><mi>&omega;</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><mfrac><mrow><mo>(</mo><msub><mi>Q</mi><mi>E</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>-</mo><msub><mi>Q</mi><mi>H</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>)</mo></mrow><msub><mi>Q</mi><mrow><mi>E</mi><mo>,</mo><mi>d</mi><mi>s</mi></mrow></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>&delta;</mi><mi>p</mi></msub></munderover><msub><mi>&omega;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mfrac><mrow><mo>(</mo><msub><mi>N</mi><mi>P</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>-</mo><mn>100</mn><mo>)</mo></mrow><mn>100</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><msub><mi>&omega;</mi><mn>1</mn></msub><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&Element;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&omega;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&omega;</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mi>min</mi></msub><mo>&lt;</mo><mo>=</mo><mi>u</mi><mo>&lt;</mo><mo>=</mo><msub><mi>u</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>|</mo><mi>&Delta;</mi><mi>u</mi><mo>|</mo><mo>&lt;</mo><mo>=</mo><msub><mi>&Delta;u</mi><mi>max</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>N</mi><mrow><mi>p</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>&lt;</mo><mo>=</mo><msub><mi>N</mi><mi>p</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>&lt;</mo><mo>=</mo><msub><mi>N</mi><mrow><mi>p</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>N</mi><mrow><mi>g</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>&lt;</mo><mo>=</mo><msub><mi>N</mi><mi>g</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>&lt;</mo><mo>=</mo><msub><mi>N</mi><mrow><mi>g</mi><mo>,</mo><mi>max</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mn>4</mn></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>&lt;</mo><mo>=</mo><msub><mi>T</mi><mrow><mn>4</mn><mo>,</mo><mi>max</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub><mo>&lsqb;</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>&rsqb;</mo><mo>&gt;</mo><mo>=</mo><msub><mi>S</mi><mrow><mi>m</mi><mi>c</mi><mo>,</mo><mi>min</mi></mrow></msub></mrow></mtd></mtr></mtable><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><msub><mi>&delta;</mi><mi>u</mi></msub></mrow></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000951332710000011.GIF" wi="1213" he="703" /></maths>其中,N<sub>p</sub>、N<sub>g</sub>、T<sub>4</sub>、S<sub>mc</sub>、Q<sub>E</sub>分别为发动机动态参数模型在线预测的动力涡轮转子转速、燃气涡轮转子转速、涡轮前温度、压气机喘振裕度、发动机输出扭矩;Q<sub>H</sub>为直升机需求扭矩模型在线预测的直升机需求扭矩;Q<sub>E,ds</sub>为发动机在设计点所提供扭矩;t<sub>0</sub>为自旋恢复阶段的起始时刻;δ<sub>p</sub>为预测时域;ω<sub>1</sub>和ω<sub>2</sub>为权重系数;F(*)为在[0,1]区间内的单调递减函数;x为优化变量,在k时刻,x={u(k),u(k+2),…,u(k+δ<sub>u</sub>‑1)},其中u为控制变量,δ<sub>u</sub>为控制时域。
地址 210016 江苏省南京市白下区御道街29号