基于城轨列车单质点动力学模型的滑模自适应速度跟踪控制 (1).pdf

基于城轨列车单质点动力学模型的滑模自适应速度跟踪控制伍文豪1，李润梅1，熊刚2（1.北京交通大学电子信息工程学院，北京100044；2.中国科学院自动化研究所多模态人工智能系统全国重点实验室，北京100190）摘要：实现城市轨道交通（简称：城轨）列车自动驾驶的关键是实现城轨列车运行速度曲线的跟踪控制。文章通过对城轨列车的运行进行受力分析，建立城轨列车单质点动力学模型，并设计了滑模自适应速度跟踪控制器，实现对城轨列车运行速度曲线的跟踪控制。以广州 22 号线城轨列车数据为研究对象，进行仿真试验，验证了滑模自适应速度跟踪控制器的有效性。关键词：城市轨道交通；单质点动力学模型；滑模自适应控制；速度跟踪控制；受力分析中图分类号：U284.48:U231:TP39文献标识码：ADOI：10.3969/j.issn.1005-8451.2023.08.17Sliding mode adaptive speed tracking control based on single-particledynamic model of urban rail trainWUWenhao1，LIRunmei1，XIONGGang2(1.SchoolofElectronicandInformationEngineering,BeijingJiaotongUniversity,Beijing100044,China；2.StateKeyLaboratoryofMultimodalArtificialIntelligenceSystems,InstituteofAutomation,ChineseAcademyofSciences,Beijing100190,China)Abstract:Thekeytoimplementingautomaticdrivingofurbanrailtransittrainsistotrackandcontrolthespeedcurveofurbanrailtrainoperation.Thispaperanalyzedtheforcesactingontheoperationofurbanrailtrain,establishedasingle-particledynamicmodelofurbanrailtrain,anddesignedaslidingmodeadaptivespeedtrackingcontrollertoimplementtrackingcontrolofthespeedcurveofurbanrailtrainoperation.ThepaperalsotookthedataofurbanrailtrainonGuangzhouLine22astheresearchobjectandconductedsimulationexperimentstoverifytheeffectivenessoftheslidingmodeadaptivespeedtrackingcontroller.Keywords:urbanrailtrain；single-particledynamicsmodel；slidingmodeadaptivecontrol；speedtrackingcontrol；forceanalysis列车自动驾驶（ATO，AutomaticTrainOperation）是指应用人工智能、智能控制等技术，实现列车自动启动、牵引、惰性、制动等基本的列车运行作业，保障列车的准点发车、安全运行和精准停车1。ATO分为 4 个自动化等级（GoA，GradeofAutomatic），分别为：司机监控列车运行（GoA1）、半自动列车运行（GoA2）、无司机驾驶（GoA3）和无人值守列车运行（GoA4）2。目前，在我国已开通的 350 余条城市轨道交通（简称：城轨）线路中，由于地形、气候和人口分布等复杂因素的影响，只有 30 余条线路实现了 GoA4 级自动驾驶，ATO 系统应用普及率较低，因此，对 ATO 的相关研究已成为列车运行控制系统（简称：列控系统）研究的热点。ATO 需要对城轨列车运行速度曲线进行跟踪控制，来实现城轨列车的自动驾驶。实现城轨列车运行速度曲线跟踪控制的关键在于建立准确的城轨列车动力学模型和选择合适的控制算法。城轨列车动力学模型可分为单质点模型和多质点模型。单质点动力学模型较简单，适用于城轨列车的整体控制和调度，但对于城轨列车整体的输入饱和约束及控制时滞问题考虑较少；多质点动力学收稿日期：2023-05-20基金项目：中国国家铁路集团有限公司科技研究开发计划（L2022X002）作者简介：伍文豪，在读硕士研究生；李润梅，教授。智慧城轨Smart Urban Rail Transit第32卷 第8期Vol.32 No.8文章编号：1005-8451（2023）08-0094-06RCA942023.08 总第 317 期模型可体现城轨列车车厢间的相互作用力和相对位移关系，适用于对每一节车厢的控制，但由于其维数增大，导致计算复杂，难以融入复杂的非线性特性和系统时滞。文献 3 对高速列车进行单质点动力学建模，但是未考虑到速度延时及输入饱和约束；文献 4 利用牛顿第二定律分别建立了高速列车的三质点和八质点动力学模型，并根据得到的各车厢间相互作用力进行安全分析；文献 5 利用单质点动力学模型描述高速列车的位移与速度，并基于高速列车多质点模型建立了一种多级牵引与制动的高速列车多质点单位移模型，保留了多质点特性且具有低维数特点，但并未考虑高速列车的速度时滞特性。在控制算法的选择方面，列车在实际运行中常遇到参数突变和扰动的情况，经典控制算法已无法满足控制需求。滑动模态（简称：滑模）控制算法在实际工程中逐渐得到推广应用，具有快速响应、对参数变化及扰动不灵敏、物理实现简单等优点，但对被控对象模型的精准性有较高要求。因为城轨列车实际运行环境复杂多变，其动力学模型中部分参数无法直接测量，具有不确定性，针对城轨列车动力学模型中的不确定性，有学者利用自适应方法对城轨列车动力学模型中存在的不确定、非线性特性的参数进行在线估计，并与期望的城轨列车速度指标进行对比，以此调整对城轨列车牵引力或制动力的控制策略，实现对城轨列车动力学模型不确定参数的实时补偿6。基于上述研究，本文以城轨列车运行速度曲线的跟踪控制为研究目的，建立城轨列车单质点动力学模型。结合滑模控制和自适应控制的特点，设计滑模自适应速度跟踪控制器。1 建立城轨列车单质点动力学模型将整组城轨列车当作一个刚性的无尺寸质点，利用牛顿定律对其进行受力分析，考虑城轨列车的输入饱和约束、控制延时及惯性环节，以微分方程的形式建立了城轨列车单质点动力学模型。1.1 城轨列车受力分析本文采用单质点模型对城轨列车运行的纵向运动特性进行描述，仅考虑城轨列车在平直道上运行的情况，根据牛顿第二定律对城轨列车进行受力分析，如图 1 所示。F城轨列车fNMg城轨列车运行方向图1城轨列车受力分析FfNMg图 1 中，为城轨列车牵引电机/制动装置在运行过程中提供的牵引力/制动力，牵引力和制动力不能同时存在；为城轨列车受到的总运行阻力，包括基本运行阻力和附加阻力；为轨面提供的支持力；为城轨列车的重力。1.1.1城轨列车牵引力和制动力MF1F2F1城轨列车牵引力是驱动列车向前运行的力，其产生的基本原理如图 2 所示，由传动系统对车轮产生旋转力矩而产生牵引力。为与大小相同、方向相反的相互作用力。图2城轨列车牵引力基本原理在 ATO 中，城轨列车牵引力的大小由列控系统计算产生。如要增大城轨列车牵引力，则要增加电动机的牵引功率。制动力是城轨列车运行过程中阻碍其运行、可调节的外力，只能在牵引力为零时产生，且与城轨列车运行速度方向相反。牵引力仅产生在动力单元的动轮与轨道间，而制动力产生在所有车厢的动轮与轨道间。1.1.2基本运行阻力城轨列车受到的基本阻力包括机械阻力和气动阻力，其中，机械阻力与城轨列车实际速度成正比，第32卷 第8期伍文豪等：基于城轨列车单质点动力学模型的滑模自适应速度跟踪控制智慧城轨RCA2023.08 总第 317 期95气动阻力与城轨列车实际速度的平方成正比。但由于城轨列车运行中影响基本阻力的因素较复杂，在实际运用中很难用理论公式进行计算，通常采用Davis 公式7，可表示为fb(v)=a+bv(t)+cv(t)2（1）v(t)km/ha、b、c式中，fb(v)为基本运行阻力；表示城轨列车的运行速度，单位为；为基本阻力系数，其系数值主要与城轨列车类型、结构及运行条件（风速、温度等）有关。1.1.3附加阻力城轨列车在某些特殊场景下（坡道、弯路、隧道、恶劣气候等）运行时，会受到除基本阻力外的阻力，称为附加阻力。附加阻力的计算公式为fa=wi+wr+ws（2）wiwrwsws式中，为坡道附加阻力；为曲线附加阻力；为隧道附加阻力。坡道附加阻力指城轨列车在坡道状况下运行时，自身重力沿坡道斜面产生的分力，上坡时为正值（阻力），下坡时为负值（加速力）；曲线附加阻力指城轨列车在曲线上运行比在直线上运行时额外增加的阻力；隧道附加阻力指城轨列车进入隧道时所受到的空气阻力。为方便计算，本文只考虑城轨列车在平直道上运行的场景，忽略其受到的坡道附加阻力和曲线附加阻力。因此，城轨列车受到的附加阻力即为隧道附加阻力，隧道附加阻力的经验公式8为fa=ws=0.000 13mgLs（3）famgLs式中，的单位为N；0.00013 为隧道附加阻力相关系数的一般取值；为城轨列车的质量，单位为kg；为重力加速度；为隧道长度，单位为m。1.2 列车单质点动力学模型根据牛顿定律，城轨列车单质点动力学模型可用微分方程描述，公式为dS(t)dt=S(t)=v(t)mdv(t)dt=m v(t)=F(t)fa fb(v)（4）S(t)F(t)式中，t 为城轨列车的运行时间，为城轨列车位移；为牵引力或制动力输入。在实际工程中，由于信号传输和机械传导，控制输入需要经过一定的延时才能作用于执行机构，即城轨列车的电机，其输出量（城轨列车加速度）的变化需要一个过程，可用一阶惯性环节来描述。城轨列车电机牵引或制动过程可用模型框图描述，如图 3 所示。()F t-Se1()F tK()a t1+1Ts1()at1s()v t图3列车电机牵引/制动模型框图sF(t)F1(t)esF(t)F1(t)图 3 中，表示对微分方程进行拉普拉斯变换的复变量；为通过计算得到的输入，即期望控制输出；为实际控制输出；为延时环节，为时滞参数；与间的关系可表示为F1(t)=F(t)（5）Ka(t)a1(t)1Ts+1T1sa1(t)v(t)为常数，表示比例环节；为期望加速度，为实际加速度；为一阶惯性环节，为惯性环节时间常数；为积分环节，实际加速度经积分环节得到城轨列车运行速度。在城轨列车实际运行中，应考虑对牵引力或制动力的饱和约束，即存在上下界，公式为F1(t)=FtF1(t)FtF1(t)Fb F1(t)FtFbF1(t)0s(x)0式中，为滑模参数，需满足。取 为 的估计值。（3）选择合适的趋近律。常用的趋近律有等速趋近律和指数趋近律。指数趋近律相较于等速趋近律具有较快的趋近速度和较小的抖振，因此，本文选取指数趋近律进行滑膜自适应速度跟踪控制器的设计，指数趋近律的表达式为 s(x)=sgn(s(x)ks(s(x)（14）k 0且 k 0式中，和 分别为等速趋近项和指数趋近项的参数，满足。sgn(s(x)sat(s(x)（4）为进一步减小抖振，将公式（14）中的符号函数项改进为饱和函数项。其表达式为sat(s(x)=1s(x)ks(x)|s(x)|,1s(x)0（18）0u 0min,maxmaxmin因需要满足，为防止出现控制输入 过大或的情况，需要对自适应律进行修正，使得的变化在范围内，和分别为不确定参数 的上下界。本文采用一种映射自适应算法，对公式（18）所示的自适应律进行修正，修正后的自适应律为=Proj(s)=0 max且 00 min且 0=0 0=0即当 大于最大值时，如有继续增大的趋势，即，则取；当 小于最小值时，如有继续减小的趋势，即，则取；其他情况则保持原自适应律不变。3 仿真分析本文选用广州 22 号线某三站两区间的 ATO 数据，作为列车运行期望速度曲线，在 MATLAB 的Simulink 中12，使用本文设计的滑模自适应速度跟踪控制器对该路线区间内的城轨列车运行速度曲线进行跟踪控制仿真。3.1 仿真条件TATO数据采样总时长为 340s，采样间隔为0.05s，共计6800 组，拟合后的城轨列车运行期望速度曲线如图 5 所示。在 Simulink 仿真中设置固定步长为0.01s；延时环节中时滞参数 为 0.2s；一阶惯性环节中时间常数 为 0.4。城轨列车单质点动力学模型相关参数设定为：m1m2第 1 区间运行时城轨列车的质量为 400t（空载情况下为 378t），第 2 区间运行时城轨列车的质量为 450t；城轨列车最大牵引/制动力均为550kN。ckminmax滑模自适应速度跟踪控制器相关参数设定为：滑模参数 为 0.5；趋近律中等速趋近项参数 为 0.5，指数趋近项参数 为 0.1；自适应律中自适应参数为 0.2；参数 的上下界和分别为 0.01m/s2和0.15m/s2。3.2 仿真结果分析a=9.888,b=0.05,c=0.001 95为进行对比，本文使用滑模控制算法设计了用于对比的列车运行速度跟踪控制器，其基本运行阻力参数设定为，其他参数与 3.1 中设定的参数一致，并将 2 种速度跟踪控制器的仿真结果与期望运行速度进行对比。2 种速度跟踪控制器输出的城轨列车实际运行速度曲线与期望运行速度曲线对比如图 6 所示，城轨列车运行速度跟踪误差如图 7 所示。-11()sat-s(x)s(x)图4饱和函数示意速度/(m s-1)时间/s图5列车运行期望速度曲线期望运行速度曲线滑膜速度跟踪控制滑膜自适应速度跟踪控制速度跟踪误差/(m s-1)时间/s图6城轨列车实际运行速度曲线与期望运行速度曲线对比智慧城轨2023年8月RCA982023.08 总第 317 期由图 6 和图 7 可以看出，在本文设计的滑模自适应速度跟踪控制器的控制下，城轨列车的速度跟踪误差较小，最大跟踪误差仅为 0.7m/s，城轨列车运行速度的最大允许误差在2m/s左右。因此，本文设计的速度跟踪控制器能够满足设计要求。同时，对比 2 种控制器的速度曲线和跟踪误差可发现，相比滑模速度跟踪控制器，本文的滑模自适应速度跟踪控制器的误差收敛速度更快，速度跟踪效果更好。4 结束语本文针对城轨列车自动驾驶的速度跟踪控制问题，考虑了输入饱和约束、系统延时及惯性环节，建立了城轨列车单质点动力学模型，并设计了滑模自适应速度跟踪控制器。仿真结果表明，本文设计的滑模自适应速度跟踪控制器在速度跟踪误差和收敛速度方面都有着较好的控制效果。参考文献杜海宾，王瑞，姜正，等.货运机车自动驾驶系统贯通试验控制研究及应用J.铁路计算机应用，2022，31（5）：60-64.1赵阳，张萍.高速铁路自动驾驶技术研究与展望J.铁道通信信号，2019，55（S1）：75-80.2李中奇，丁俊英，杨辉，等.基于控制器匹配的高速列车广义预测控制方法J.铁道学报，2018，40（9）：82-89.3丁盼.基于多质点模型的高速列车自适应速度跟踪控制D.南昌：华东交通大学，2021.4宋琦.高速列车的鲁棒自适应及容错控制D.北京：北京交通大学，2014.5SongYD,SongQ,CaiWC.Fault-tolerantadaptivecontrolofhigh-speed trains under traction/braking failures:a virtualparameter-basedapproachJ.IEEETransactionsonIntelligentTransportationSystems,2014,15(2):737-748.6Lin X,Dong H R,Yao X M,et al.Neural adaptive fault-tolerantcontrolforhigh-speedtrainswithinputsaturationandunknowndisturbanceJ.Neurocomputing,2017(260):32-42.7袁佳希.列车在线辨识与预测控制研究D.杭州：浙江大学，2017.8任一凡，孙后环，刘陈.基于滑模控制算法的自适应跟车巡航研究J.机械与电子，2023，41（1）：59-64.9李中奇，周靓，杨辉.高速动车组数据驱动无模型自适应控制方法J.自动化学报，2023，49（2）：437-447.10刘豹，唐万生.现代控制理论M.3版.北京：机械工业出版社，2013：164-165.11刘金琨.滑模变结构控制MATLAB仿真：基本理论与设计方法M.3版.北京：清华大学出版社，2015:56-57.12责任编辑李依诺速度跟踪误差/(m s-1)时间/s图7城轨列车运行速度跟踪误差对比第32卷 第8期伍文豪等：基于城轨列车单质点动力学模型的滑模自适应速度跟踪控制智慧城轨RCA2023.08 总第 317 期99