基于HHT能量和最大Lyapunov指数的蛇行分类方法.pdf

设计与应用计算机测量与控制 ()C o m p u t e r M e a s u r e m e n t&C o n t r o l 收稿日期:;修回日期:.基金项目:国家自然科学基金项目();国家自然科学基金项目();四川省科技计划资助项目(J D T D ).作者简介:王敏(),男,山西大同市人,硕士研究生,主要从事智能化状态监测与故障诊断方向的研究.通讯作者:宁静(),女,四川成都人,博士,副教授,主要从事智能化状态监测与故障诊断、精密测控技术及系统大数据挖掘在高速列车中的应用等方向的研究.引用格式:王敏,宁静,赵飞,等基于HHT能量和最大L y a p u n o v指数的蛇行分类方法J计算机测量与控制,():文章编号:()D O I:/j c n k i /t p 中图分类号:T P 文献标识码:A基于H H T能量和最大L y a p u n o v指数的蛇行分类方法王敏,宁静,赵飞,李艳萍,陈春俊(西南交通大学 机械工程学院,成都 )摘要:现有蛇行监测标准仅考虑构架横向加速度信号的时域幅值信息,通过单一的固定值去评判列车是否蛇行,没有定量反映蛇行失稳对列车的影响程度;针对以上问题,提出基于HHT的能量法,该方法从信号频域主频的大小,频谱的集中性以及频率值等方面来判断高速列车是否发生蛇行运动;在此基础上,提出了快速蛇行收敛的概念,此类情况对车辆系统运行安全影响较小,因此是不影响车辆系统安全的信号;最后,利用最大L y a p u n o v指数表征信号的周期性,从而区分快速蛇行收敛和蛇行运动(小幅蛇行和大幅蛇行);结果表明,通过HHT能量法与最大L y a p u n o v指数相结合可以对高速列车正常运行、快速蛇行收敛、小幅蛇行和大幅蛇行的种运行状态进行定性识别,且HHT能量值能定量反映蛇行程度的大小,并通过实测线路数据验证了方法的可行性.关键词:高速列车;蛇行监测;HHT能量法;快速蛇行收敛;最大L y a p u n o v指数H u n t i n gC l a s s i f i c a t i o nM e t h o dB a s e do nHH TE n e r g ya n dL a r g e s tL y a p u n o vE x p o n e n tWANG M i n,N I NGJ i n g,Z HAOF e i,L IY a n p i n g,CHE NC h u n j u n(S c h o o l o fM e c h a n i c a lE n g i n e e r i n g,S o u t h w e s t J i a o t o n gU n i v e r s i t y,C h e n g d u ,C h i n a)A b s t r a c t:e x i s t i n gh u n t i n gm o n i t o r i n gs t a n d a r do n l yc o n s i d e r s t h e t i m ed o m a i na m p l i t u d e i n f o r m a t i o no f f r a m e l a t e r a l a c c e l e r a t i o ns i g n a l,s i n g l e f i x e dv a l u e i su s e d t o j u d g ew h e t h e r t r a i nh u n t s,w h i c hd o e sn o t r e f l e c t t h e i n f l u e n c eo f h u n t i n g i n s t a b i l i t yo n t r a i nq u a n t i t a t i v e l y A i m e da t a b o v ep r o b l e m,ae n e r g ym e t h o db a s e do nH i l b e r t H u a n g t r a n s f o r m(HHT)i s f i r s t l yp r o p o s e d t o j u d g ew h e t h e rt h eh i g h s p e e dt r a i nh u n t s f r o mt h em a i nf r e q u e n c ys i z e,s p e c t r u mc o n c e n t r a t i o na n df r e q u e n c yv a l u eo f s i g n a l f r e q u e n c yd o m a i n O nt h i sb a s i s,t h e c o n c e p t o f f a s t h u n t i n gc o n v e r g e n c e i sp r o p o s e d,w h i c hh a s l i t t l e i n f l u e n c eo n t h e s a f e t yo f v e h i c l e s y s t e mo p e r a t i o n,s ot h es i g n a l d o e sn o t a f f e c t t h e s a f e t yo f t h ev e h i c l e s y s t e m F i n a l l y,t h e l a r g e s tL y a p u n o ve x p o n e n t i su s e d t o c h a r a c t e r i z e t h ep e r i o d i c i t yo f t h e s i g n a l,t h e f a s t h u n t i n gc o n v e r g e n c e a n dh u n t i n g i n s t a b i l i t y(s m a l l a m p l i t u d eh u n t i n ga n d l a r g e a m p l i t u d eh u n t i n g)a r ed i s t i n g u i s h e d T h e r e s u l t s s h o wt h a t,t h e f o u r o p e r a t i n gs t a t e so f n o r m a l o p e r a t i o n,f a s t h u n t i n gc o n v e r g e n c e,s m a l l a m p l i t u d eh u n t i n ga n dl a r g ea m p l i t u d eh u n t i n go f t h ev e h i c l e s y s t e ma r e i d e n t i f i e dq u a l i t a t i v e l yb yc o m b i n i n g t h eHHTe n e r g yw i t ht h e l a r g e s tL y a p u n o ve x p o n e n t T h eHHTe n e r g yv a l u e sc a nr e f l e c tt h ei n t e n s i t yo fh u n t i n gm o v e m e n tq u a n t i t a t i v e l y A n dt h ef e a s i b i l i t yo ft h em e t h o di sv e r i f i e db yt h em e a s u r e dd a t ao f t h e l i n e K e y w o r d s:h i g hs p e e dt r a i n;h u n t i n gm o n i t o r i n g;HHTe n e r g y;f a s th u n t i n gc o n v e r g e n c e;l a r g e s tL y a p u n o ve x p o n e n t引言蛇行运动是车辆动力学系统的核心问题之一,车辆在运行过程中发生收敛较慢的小幅蛇行或者剧烈的大幅蛇行都会严重地影响车辆运行安全,因此对车辆小幅蛇行和大幅蛇行的在线监测都至关重要.学者们通过大量的理论研究总结出轮轨参数和悬挂系统参数对车辆系统蛇行运动的影响规律.在这些理论研究的成果上,更多的学者致力于探索如何将这些成果应用到高速列车蛇行监测领域.宋兴武提出横向位移峰值方法(LMP),通过计算轮对和构架的横向加速度确定轮对横移量的峰值来识别车辆的蛇行运动.朴明伟等发现当车辆系统服从超临界分岔时,随着运行速度的增加,转向架会发生小幅值蛇行,且随着速度的增加,转向架的蛇行幅值变大,最终会出现大幅蛇行失稳.蔡里军总结了大量的实测数据认为当构架横向加速度峰值连续次达到或超过m/s,车辆发生了蛇行运投稿网址:w w wj s j c l y k z c o m计算机测量与控制第 卷 动.方明宽等考虑到构架横向加速度数据的非平稳性,提出NK J A D E方法对多传感器横向加速度数据进行特征融合,从而对车辆的正常运行,小幅蛇行,大幅蛇行状态进行识别.S u n 发现当车辆出现蛇行运动时,构架与车体的加速度信号存在相位延迟,通过对车体和构架的横向加速度信号进行互相关分析,将互相关指标作为判别车辆蛇行失稳的阀值.曾元辰、赵飞 等发现车辆出现蛇行运动时,轮对、构架和车体的横向加速度都具有较强的周期性,提出利用信号周期性特点来识别车辆的蛇行运动.崔万里 通过样本熵理论以及等距映射算法(I S OMA P)对高速列车构架横向加速度信号进行特征提取,然后通过L S S VM进行车辆小幅蛇行异常识别,达到了较好的识别效果.叶运广 通过ME EMD进行特征提取与使用L S S VM对车辆的运行状态进行特征识别亦取得很好的效果.冉伟 因而提出了一种基于E EMD S V D L T S A的特征提取框架,用于识别高速列车小幅蛇行运动的演变趋势.陈杨 通过建立道岔模型,研究高速列车在小幅蛇行的状态下通过道岔的演变规律.宁云志 提出基于 D C NN和C GAN的预测方法对小幅蛇行进行预警,实现了蛇行失稳的提前预警.王晓东 利用 DC N N自适应地对构架横向加速度数据进行特征提取,使用L S TM对车辆的正常、小幅蛇行和大幅蛇行状态进行分类,取得了很好的分类效果.但是目前,列车实际运行过程中对蛇行运动的监测还存在以下不足:问题一:现有的蛇行监测标准将构架横向加速度信号峰值连续次达到m/s时视为大幅蛇行,但实际过程中即使构架横向加速度没有达到现有的监测标准也会出现大幅的振动.此外,现有标准没有定量反映蛇行失稳对车辆的影响程度,仅用一个固定值来评判是否产生蛇行,具有一定的局限性.问题二:某些车辆受到部分轨道不平顺周期性的干扰,虽然发生了明显的谐波振动,但只是运行中的极少时刻,且会快速收敛,不会形成稳定的周期蛇行运动,对运行安全性没有较大的影响,将此类状态认为快速蛇行收敛.本文针对以上问题,提出HHT(H i l b e r t H u a n gt r a n s f o r m)能量法从信号频域主频的大小与频谱的集中性以及频率值等方面来判断高速列车是否存在蛇行运动,对两个蛇行运动状态(小幅蛇行和大幅蛇行)的蛇行程度的大小进行定量的评估.利用最大L y a p u n o v指数进一步表征信号的周期性从而将快速蛇行收敛和蛇行运动(小幅蛇行和大幅蛇行)区分开来.最后通过HHT能量法与最大L y a p u n o v指数相结合,从高速列车实际运行的需求出发,对车辆系统进行正常运行、小幅蛇行、快速蛇行收敛、大幅蛇行的定性识别和蛇行程度大小的定量分析,以实现对蛇行运动的具体监测.仿真模型建立 高速列车动力模型本文使用S I MP A C K软件建立国内某型号高速列车的整车动力学模型,该车辆的基本参数如表所示.该车辆包括个车体、个构架、个轮对、个轴箱、一系悬挂系统(一系钢簧、一系垂向减振器、转臂式轴箱定位装置)和二系悬挂系统(空气弹簧、二系垂向减振器、二系横向减振器、抗蛇行减振器、抗侧滚扭杆、牵引拉杆、横向止挡).在建立仿真模型时,忽略了车辆系统各部件的弹性变形,将车体、构架、轮对、轴箱等部件假定为刚体.车体、构架和轮对有个自由度(纵向、横向、垂向、侧滚、点头、摇头),轴箱只有个点头自由度.因此,该车辆动力学模型共有 个自由度.表高速列车动力学参数表参数数值转向架中心距 m轴距 m车轮滚动圆横向跨距 m车轮滚动圆直径 m轮对内侧距 m车轮外形S G钢轨外形CHN 车体质量 t车体侧滚转动惯量 tm车体点头转动惯量 tm车体摇头转动惯量 tm构架质量 t构架侧滚转动惯量 tm构架点头转动惯量 tm构架摇头转动惯量 tm轮对质量 t轮对侧滚转动惯量 tm轮对点头转动惯量 tm轮对摇头转动惯量 tm车体重心距轨面高 m构架重心距轨面高 m 建立仿真工况类型关于轮对踏面的选择,踏面类型对高速列车系统极限环分岔影响显著,为了模拟高速列车的不同运行状态,同时使高速列车的运行状态更符合实际,本文分别选用了在一个旋 修 周 期(k m)内种 不 同 行 驶 里 程 的S G实测踏面.踏面至踏面对应行驶的里程依次增大,踏面的磨损程度依次增加.其中,踏面为列车行驶 k m的行驶里程左右时,测量轮对外形轮廓获得的实测磨损踏面;踏面为列车行驶 k m的行驶里程左右时,测量轮对外形轮廓获得的实测磨损踏面;踏面为列车行驶 k m的行驶里程左右时,测量轮对外形轮廓获得的实测磨损踏面.关于抗蛇行减振器的参数设定,抗蛇行减振器的参数对高速列车系统的临界速度和分岔类型也有很大影响.本文为充分模拟高速列车的不同运行状态,使用种不同阻尼特性的抗蛇行减振器.抗蛇行减振器的阻尼特性曲线如图所示,从抗蛇行减振器到抗蛇行减振器对应的卸荷投稿网址:w w wj s j c l y k z c o m第期王敏,等:基于HHT能量和最大L y a p u n o v指数的蛇行分类方法 力和卸荷速度依次增加.使用不同行驶里程的S G实测踏面和不同阻尼特性的抗蛇行减振器模拟高速列车(速度:k m/h)种工况:)工况:使用踏面和抗蛇行减振器仿真,模拟正常运行;)工况:使用踏面和抗蛇行减振器仿真,模拟小幅蛇行;)工况:使用踏面和抗蛇行减振器仿真,模拟大幅蛇行和快速蛇行收敛.图抗蛇行减振器阻尼系数图图车辆系统极限环分岔图根据以上工况,计算得到各工况下的极限环分岔图,如图所示.工况的高速列车动力学模型分岔类型为超临界分岔,在此工况下车辆模型在 k m/h区间运行时,车辆处于正常运行状态,如图所示.工况的高速列车动力学模型分岔类型同样为超临界分岔,在此工况下车辆模型在 k m/h区间运行时,车辆的轮对处于小幅振动状态,此时车辆处于小幅周期蛇行状态,随着速度增加,横向振动幅值也逐渐增大,如图所示.工况下高速列车动力学模型分岔类型为亚临界分岔,在此工况下,车辆模型在 k m/h区间运行时,车辆系统极限环处于不稳定区域,受轨道不平顺的影响,轨道不平顺激扰大时,车辆系统平衡点收敛为稳定极限环值产生大幅蛇行运动,如图所示.当轨道不平顺激扰不足以让车辆系统产生稳定的极限环时,车辆系统可能会出现短暂的蛇行运动随后车辆系统平衡点收敛为,无法形成稳定的周期运动,即快速蛇行收敛,如图所示.图正常运行图小幅蛇行图大幅蛇行图快速蛇行收敛通过建立上述工况,确定了高速列车动力学模型种运行状态的仿真条件.为构建大量的仿真数据集,将车辆运 行 速 度 分 为(k m/h、k m/h、k m/h、k m/h、k m/h、k m/h)个速度级,对上文用到的轨道不平顺分别设置了 、三种幅值比例系数,扩充轨道不平顺数据集.每种运行状态对应个速度投稿网址:w w wj s j c l y k z c o m计算机测量与控制第 卷 级,种幅值比例轨道不平顺,共 组仿真条件,种运行状态共对应 组仿真数据.对每组仿真得到的构架横向加速度信号随机取 个样本,每个样本段取样时长为s,采样频率为 H z,总计 个样本段的仿真数据集,供后续实验使用.理论方法 H H T能量理论大量的研究表明蛇行信号和正常信号的主频不同,基于这一特点,通过研究信号所有频率成分在频域段内的变化,可以更加准确地找到正常信号和蛇行信号在频域内的不同特征.再结合信号的时域特征进行分析,可以准确地得知信号的时频信息,从而更加全面地分析信号特征.孙新建 通过HHT能量方法对非平稳的爆炸震动信号进行分析,对希尔伯特边际谱进行频域积分计算HHT能量,研究结构体的爆炸震动损伤.HHT能量法结合了信号频域主频的大小与频谱的集中性以及频率值和时域幅值的相关信息,因此本文采用此方法研究蛇行运动,首先对蛇行信号和正常信号的主频及频域范围进行分析,确定蛇行频域主要范围.然后通过HHT边际谱统计正常信号和蛇行信号在特定频率段范围内的频域信息,最后结合信号的时域幅值和频域信息,通过HHT能量法全面反映蛇行信号和正常信号不同的时频特征,从而对车辆运行状态准确的识别.HH T具体变换原理如下:对原始信号x(t)进行经验模式分解得到n个I M F(i n t r i n s i cm o d e f u n c t i o n)信号分量:x(t)nicirn()式中,ci为信号的基本模式分量,rn称为残余函数.对ci进行H i l b e r t变换得到:Gi(t)ci()td()式中,为时间延长间隔.构造解析信号zi(t):zi(t)ci(t)j Gi(t)ai(t)ej i(t)()式中,幅值函数为ai(t),相位函数为i(t),ci分量信号的H i l b e r t谱为Hi(,t):Hi(,t)R Pai(t)ej i(t)dt()式中,R P表示取实部.对Hi(,t)进行时域上的积分得到边际谱hi():hi()THi(,t)dt()式中,T为信号时间长度.对hi()(i,n)边际谱分析,对主频在H z以上的边际谱进行叠加求得最终的边际谱h():h()jkhk()()式中,hk()(k,j)为主频在H z以上的边际谱,在频域内对h()进行积分得到最终的HHT能量值.Eh()d()最大李雅普诺夫指数理论蛇行信号和正常信号除了时频上的差异,其周期性特征也有着明显的区别,所以有必要对信号的周期性特征进行分析.李雅普诺夫方法是研究运动稳定性的重要方法,对于非线性动力学系统,李雅普诺夫指数可以从状态变量的时间序列中计算出来.最大李雅普诺夫指数为相空间相邻轨迹的平均指数发散率的数值特征,常用来判定系统的混沌特性,因此其可以用来量化信号的周期性.最大李雅普诺夫指数计算步骤如下.对于时间序列 x(t)构造n维空间Rn,T为延迟时间,Tkt(k,),t为时间间隔即:X(t)x(t),x(tT),x(t(n)T)T()N维相空间中的某一时刻,两条邻近轨迹之间的距离可以分解在n个不同的方向,这个n个不同方向的上的距离增长率是不同的,每一个增长率就是一个李雅普诺夫指数.系统可以写成n个自治一阶微分方程组的形式即:xf(x)xRn()取两条邻近轨迹L和L,起始点分别为x和y两起始点的距离为dyx,经过t时间后分别运动到x和y,此时距离为dyx,如此循环下去经过mt后得到m个di(i,m),diyixi,最大李雅普诺夫指数即:l i mmm tmil ndid()基于仿真数据的蛇行状态研究 研究步骤由于实测数据较少,本文利用仿真数据进行分析,再通过实测数据进行验证,通过HHT能量法和最大L y a p u n o v指数法对蛇行状态进行分类且定量反映蛇行程度的大小,其具体操作如下:)仿真得到不同运行状态下的构架横向加速度信号数据,采用 H z对构架横向加速度进行滤波,提取时间长度为s的信号进行分析.)通过对大量的仿真数据进行EMD(e m p i r i c a lm o d ed e c o m p o s i t i o n)分析,得出蛇行信号主频和能量频带基本在H z以上.)对主频在H z以上I MF(i n t r i n s i cm o d ef u n c t i o n)分量信号的H i l b e r t边际谱进行叠加求和.)在频域内对最终得到的H i l b e r t边际谱进行积分得到最终的HHT能量值.)利 用S VM(s u p p o r tv e c t o r m a c h i n e s)分 类 得 到HHT能量阈值,根据HHT能量值将正常运行、小幅蛇行、大幅蛇行区分开,且通过HHT能量值的大小反映蛇行程度的大小.)利用最大L y a p u n o v指数对高速列车不同运行状态信号的周期性进行分析.)利用S VM分类得到最大李雅普诺夫指数阈值,根据最大L y a p u n o v指数值将快速蛇行收敛与蛇行运动(小幅投稿网址:w w wj s j c l y k z c o m第期王敏,等:基于HHT能量和最大L y a p u n o v指数的蛇行分类方法 蛇行和大幅蛇行)区分开来.H H T能量分析对图、的构架横向加速度信号进行示例分析.如表所示,正常运行信号的主频和能量频带分布比较均匀,主频H z以上I MF信号分量能量占比,H z以下I MF信号分量能量值占比.如表和表所示,蛇行信号(小幅蛇行、大幅蛇行)主频高的I MF分量能量占比高,蛇行频率主频和能量频带基本在H z以上,信号分量能量值均占比达到 以上.如表所示,快速蛇行收敛信号因为包含蛇行特征,I MF分量主频在H z以上的占比也高达.对主频在H z以上I MF分量信号边际谱进行叠加,得到最终的H i l b e r t边际谱,最后在频域内对H i l b e r t边际谱进行积分得到最终的HHT能量值.表正常运行(H z以上HHT能量所占百分比为 )信号分量边际谱主频频带HHT能量百分比I m f I m f I m f I m f I m f 表小幅蛇行(H z以上HHT能量所占百分比为 )信号分量边际谱主频频带HHT能量百分比I m f I m f I m f I m f I m f 表大幅蛇行(H z以上HHT能量所占百分比为 )信号分量边际谱主频频带HHT能量百分比I m f I m f I m f I m f 表快速蛇行收敛(H z以上HHT能量所占百分比为 )信号分量边际谱主频频带HHT能量百分比I m f I m f I m f I m f I m f 对仿真数据集(个样本)的正常运行、小幅蛇行和大幅蛇行类数据的HHT能量值进行分析,利用S VM方法计算得出其分类阈值,分别为 和 ,如图所示.当信号HHT能量值小于 时,将其视为正常运行,当信号HHT能量值大于 小于 时视为小幅蛇行,当信号HHT能量值大于 时视为大幅蛇行.HHT能量值的大小反映了蛇行程度的大小,即HHT能量值越大蛇行程度越剧烈,HHT能量值越小蛇行程度越小.如图所示,在此类数据上加入快速蛇行收敛数据,发现快速蛇行收敛HHT能量值和小幅蛇行、大幅蛇行HHT能量值重叠,表明HHT能量值指标无法区分这种运行状态,因此需要补充新的指标将其区分.图三类仿真数据HHT能量统计图图四类仿真数据HHT能量统计图大幅蛇行和小幅蛇行的最大L y a p u n o v指数值信号周期性较强,正常运行和快速蛇行收敛的信号周期性较弱.因此可通过最大L y a p u n o v指数值将快速蛇行收敛和正常运行归为周期性较弱的一类,将小幅蛇行和大幅蛇行归为周期性较强的一类.对仿真数据集(个样本)的正常运行、小幅蛇行、大幅蛇行、快速蛇行收敛数据的最大L y a p u n o v指数进行分析,如图所示.利用S VM方法计算得出其分类阈值 ,李雅普诺夫指数值低于 时,将其认为属于正常运行和快速蛇行收敛一类,最大L y a p u n o v指数值大于 时将其认为小幅蛇行和大幅蛇行一类.通过最大L y a p u n o v指数法和HHT能量法结合,将正常运行、小幅蛇行、大幅蛇行、快速蛇行收敛区分开来.当信号HHT能量值小于 时,将其视为正常运行,当信号HHT能量值大于 小于 且最大L y a p u n o v指数值小于 时视为小幅蛇行,当信号HHT能量值大于 且最大L y a p u n o v指数值小于 时视为大幅蛇行,当信号HHT能量值大于 且最大L y a p u n o v指数值大于 时视为快速蛇行收敛.此方法既考虑了蛇行的周期性特点,也考虑了蛇行信号的频域特性包括频率主频的大小、投稿网址:w w wj s j c l y k z c o m计算机测量与控制第 卷 频谱的集中性以及频率值等特点,从而可以定性识别车辆系统不同的运行状态,且能够定量表示蛇行信号的能量值从而反映蛇行程度的大小.图仿真数据最大L y a p u n o v指数统计图实测数据验证为验证本文所提方法的正确性,采用实测数据对其进行验证.图 为某高速列车速度在 k m/h区间的部分实测构架横向加速度数据,时长为 s,单个样本数据长度为s,样本总数为 .图 实测构架横向加速度信号用HHT 最大L y a p u n o v指数方法对实测数据进行分析,图 为实测数据的HHT能量值,图 为实测数据的最大L y a p u n o v指数图.图 为局部信号监测放大图.图 实测数据HHT能量图从图 的识别结果中可以观察到,s这一段信号,传统的方法会将其识别为正常运行,但这段信号出现了明显的周期性特征,本文所提方法将其识别为小幅蛇行状态,这段小幅蛇行识别可以起到预警作用.s这一段信号,其加速度峰值未达到目前所使用的蛇行报警标准(连续个周期加速度峰值达到m/s),传统的方法图 实测数据L y a p u n o v指数图图 小幅大幅监测图图 小幅快速蛇行收敛监测图图 正常小幅监测图会将其识别为正常运行,但通过观察可以看出,该信号峰值(m/s之间)接近目前所使用的蛇行报警标准,且出现了明显的周期性,本文所提方法将其识别为大幅蛇行状态,故本文认为其识别结果合理.从图 的识别结果可以看到 s这一段信号,传统的方法无法识别其小幅蛇行状态,但其具有明显的周期性特征,s这一段信号本方法识别为快速蛇行投稿网址:w w wj s j c l y k z c o m第期王敏,等:基于HHT能量和最大L y a p u n o v指数的蛇行分类方法 收敛,从图中可以看出其周期性特征逐渐消失,证明了本方法的合理性.从图 的识别结果看到,s这一段信号从正常运行变为小幅蛇行的状态.其信号周期性特征逐渐变强,幅值不断上升,车辆系统出现小幅蛇行.传统的方法会将这一段信号整体识别为正常运行,本方法可以准确地监测这一变化过程.本文所提方法准确地识别了车辆蛇行运动的不同状态,能够及时采取措施预防列车发生严重的蛇行运动,保证列车运行安全.综上,本文利用实际数据验证了方法具有可行性,能够在车辆运行过程中实现更为准确的在线监测.结束语现有的高速列车蛇行失稳标准仅针对大幅蛇行且单一的从信号时域的角度 分 析.本 文将HHT能 量值 和 最 大L y a p u n o v指数值两个指标进行了结合,综合考虑了信号的时域、频域和周期性特点,以此作为区分高速列车不同蛇行状态的参考指标.本方法重点并非确定不同蛇行运行状态的准确阈值,由于仿真数据的局限性,本文通过仿真数据所确立的阈值标准存在一定的误差.阈值的准确性问题可通过台架试验或线路实测等方式对真实数据进行补充,再结合迁移学习等机器学习算法,提高实际监测中阈值的准确性.参考文献:孙丽霞高速列车横向运动稳定性和动态脱轨理论分析及评价方法研究 D北京:中国铁道科学研究院,P O L A CH O A p p l i c a t i o no fn o n l i n e a rs t a b i l i t ya n a l y s i s i nr a i l w a yv e h i c l e i n d u s t r yMS p r i n g e rB e r l i nH e i d e l b e r g,:曾京车辆系统的蛇行运动分叉及极限环的数值计算 J铁道学报,():董浩铁道车辆运动稳定性及分岔类型研究 D成都:西南交通大学,晏永车辆系统蛇行运动H o p f分岔行为研究 D成都:西南交通大学,宋欣武转向架蛇行失稳监测的横向位移识别方法研究 D成都:西南交通大学,朴明伟,梁树林,孔维刚,等高速转向架非线性稳定性及安全裕度对策 J振动与冲击,():蔡里军基于F P G A的高速列车转向 架 蛇 行 失 稳 检 测 装 置D成都;西南交通大学,N I NGJ,F ANG M K,R AN W,e ta l R a p i dm u l t i s e n s o r f e a t u r ef u s i o nb a s e do nn o n s t a t i o n a r yk e r n e l j a d e f o r t h e s m a l l a m p l i t u d eh u n t i n gm o n i t o r i n go fh i g h s p e e dt r a i n sJS e n s o r s,():S UNJF,ME L IE,C A IW B,e ta l As i g n a la n a l y s i sb a s e dh u n t i n g i n s t a b i l i t yd e t e c t i o nm e t h o d o l o g y f o r h i g h s p e e d r a i l w a yv e h i c l e sJ V e h i c l eS y s t e mD y n a m i c s,():Z E NGYC,Z HANG W H,S O NGDL An e ws t r a t e g yf o rh u n t i n ga l a r m a n d s t a b i l i t y e v a l u a t i o nf o rr a i l w a y v e h i c l e sb a s e do nn o n l i n e a rd y n a m i c sa n a l y s i sJP r o c e e d i n g so f t h eI n s t i t u t i o no f M e c h a n i c a lE n g i n e e r s,P a r tF:J o u r n a lo fR a i la n dR a p i dT r a n s i t,():赵飞,宁静,李艳萍,等基于S I MP A C K仿真分析的高速列车蛇行失稳识别方法 J现代制造工程,():崔万里,宁静,种传杰,等基于EMD I S OMA P高速列车小幅蛇行异常特征提取 J中国测试,():叶运广,宁静,种传杰,等高速列车转向架蛇行失稳的ME EMD L S S VM预测模型 J铁道学报,():冉伟,宁静,陈杨,等基于E EMD S V D L T S A的高速列车蛇行演变特征提取框架 J电子测量技术,():陈杨,宁静,王靖铭,等高速转向架小幅蛇行运动状态下过道岔研究 J计算机测量与控制,():宁云志,宁静,李艳萍,等基于数据不平衡下的高速列车小幅蛇行预测方法 J电子测量技术,():王晓东,宁静,陈春俊 DC NN和L S TM高速列车横向稳定性状态识别研 究 J中 国 测 试,():石怀龙,罗仁,曾京国内外高速列车动力学评价标准综述 J交通运输工程学报,():孙新建基于H i l b e r t能量分析的岩体爆破震动损伤研究D天津:天津大学,(上接第 页)于占东,王世笛,付莹,等感应加热动态负载阻抗匹配及功率因数角控制 J电力电子技术,():傅清并联谐振中频电源逆变角自动移相电路的应用 J电力电子技术,():陈彦桥,郑亚锋,刘吉臻,等基于动态解耦的模糊多模型协调控制 系 统 应 用 研 究 J中 国 电 机 工 程 学 报,():S T E I N GR e s p e c tt h eu n s t a b l eJ I E E EC o n t r o lS y s t e m sM a g a z i n e,():Z HOU KE s s e n t i a l so fr o b u s tc o n t r o lMP r e n t i c H a l l,王 广 雄,何朕工 业 伺 服 控 制 J电 机 与 控 制 学 报,():吴敏,桂卫华,何勇现代鲁棒控制(第二版)M长沙:中南大学出版社,姚郁,傅绍文,贺风华,等基于加权函数选择的伺服系统性能优化设计方法 J自动化学报,():姚方方,张欣雨具有随机增益的电力系统弹性负荷频率H控制 J计算机测量与控制,():投稿网址:w w wj s j c l y k z c o m