发明名称 一种动车组运行过程建模与自适应控制方法
摘要 本发明公开了一种动车组运行过程建模与自适应控制方法。所述方法针对动车组运行过程复杂、信息不完备、非线性特征明显的特性,采用数据驱动的建模方法,提出了T-S双线性模型辨识方法;依据动车组实际运行图、前方线路情况、限速条件,牵引力/制动力饱和非线性特性等约束条件建立动车组约束模型;根据上述模型设计即时学模型预测控制算法来改善多目标优化控制性能指标。本发明为动车组乘务员进行优化操纵,保证安全正点运行,改善运行舒适性,降低能源消耗提供了一套可靠的依据。本发明适用于动车组运行状态在线辨识与多目标优化控制。
申请公布号 CN102981408B 申请公布日期 2015.05.27
申请号 CN201210524440.8 申请日期 2012.12.10
申请人 华东交通大学 发明人 杨辉;张坤鹏;李中奇;付雅婷;刘杰民
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 南昌市平凡知识产权代理事务所 36122 代理人 姚伯川
主权项 一种动车组运行过程建模与自适应控制方法,其特征是,所述方法依据动车组实际运行图、前方线路情况、限速条件,牵引力/制动力饱和非线性特性约束条件建立动车组运行多目标约束优化模型,据此设计双线性模型预测控制器实现动车组优化运行;并根据上述模型设计即时学习模型预测控制算法来改善动车组多目标约束优化控制性能指标;当模型输出误差e在系统允许范围内时,T‑S双线性模型参数无需优化;当模型输出误差超过阈值时,采用即时学习算法对模型参数进行在线校正;据此动态调整双线性模型预测控制器的参数,实现模型参数和控制器参数的同时优化,减少了双线性模型预测控制器在线计算量;所述动车组运行多目标约束优化模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mi>F</mi><mo>=</mo><msub><mi>W</mi><mi>T</mi></msub><mo>|</mo><mfrac><mrow><mi>t</mi><mo>-</mo><msup><mi>t</mi><mo>*</mo></msup></mrow><msup><mi>t</mi><mo>*</mo></msup></mfrac><mo>|</mo><mo>+</mo><msub><mi>W</mi><mi>E</mi></msub><mo>|</mo><mfrac><mrow><mi>E</mi><mo>-</mo><msup><mi>E</mi><mo>*</mo></msup></mrow><msup><mi>E</mi><mo>*</mo></msup></mfrac><mo>|</mo></mrow>]]></math><img file="FDA0000666289080000011.GIF" wi="657" he="156" /></maths>(ν<sub>r</sub>‑δ)≤ν≤(ν<sub>r</sub>+δ)t≤t<sup>*</sup>U<sub>min</sub>≤U≤U<sub>max</sub>Δu<sub>min</sub>≤Δu≤Δu<sub>max</sub>式中:F表示目标函数;W<sub>T</sub>为运行时分的权重;W<sub>E</sub>为能耗权重;v<sub>r</sub>为目标速度;δ为速度跟踪误差范围;t<sup>*</sup>为运行图的运行时分;E<sup>*</sup>为区间运行期望的能耗;U<sub>min</sub>代表最大制动力;Δu为控制量增量;U<sub>max</sub>代表最大牵引力;Δu<sub>min</sub>为控制量增量的下限值;Δu<sub>max</sub>为控制量增量的上限值;所述双线性模型为:N辆编组的动车组动态特性、能量消耗和运行时分的双线性模型为:<img file="FDA0000666289080000021.GIF" wi="1406" he="347" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub></msubsup><mfrac><mn>1</mn><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac><mi>dx</mi></mrow>]]></math><img file="FDA0000666289080000022.GIF" wi="418" he="168" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub></munderover><mi>U</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>dx</mi></mrow>]]></math><img file="FDA0000666289080000023.GIF" wi="692" he="183" /></maths>式中:M为动车组等价总质量;x是动车组的运行距离;系统输入是作用于动车组上的不同手柄级位下的控制力U(x);系统输出是速度v(x);符号<img file="FDA0000666289080000024.GIF" wi="310" he="73" />算子使得U(x)和v(x)满足乘积关系;G为线路参数;x<sub>1</sub>表示起点的位置;x<sub>2</sub>表示终点的位置;t为运行时分;E(x)为动车组运行过程中的能量消耗;C<sub>1</sub>为机械阻力系数,大小为7.4×10<sup>‑3</sup>,C<sub>2</sub>为空气阻力系数,大小为1.24×10<sup>‑4</sup>。
地址 330013 江西省南昌市昌北区双港路