一种基于数据驱动的高速列车悬挂系统早期故障检测方法与流程

1.本发明公开了一种基于数据驱动的高速列车悬挂系统早期故障检测方法，属于高速列车悬挂系统故障检测技术领域。


背景技术：

2.高速列车悬挂系统起到支撑车体和转向架，同时隔离由轨道不平顺所引起的轮轨作用力，保证列车在高速运行时的稳定性和安全性，因此对其有高可靠性的要求。从2008年在主要干线上开通了时速200公里以上的crh(china railway high
‑
speed)动车组列车到现在，经过10多年来的发展，我国已经建成全世界规模最大、运营速度最快的高速铁路网。对提高我国高铁运行的安全性，针对悬挂系统早期故障检测的研究具有重要的意义。
3.高速列车悬挂系统分为主动悬挂和半主动悬挂，采用闭环控制结构。悬挂系统采用两系悬挂装置，一系悬挂装置设在轴箱和转向架构架之间，二系悬挂装置设在转向架构架和车体之间。包含有大量的部件，包括螺旋弹簧、横向/垂向阻尼器、空气弹簧、主动式作动器以及传感器。其中主动式作动器作为重要的执行器部件，对于高速列车的安全运行和乘坐舒适性至关重要；主动式作动器是根据传感器测得的车辆输出信号计算出主动悬挂系统所需的主动控制力，因此，传感器同样具有很高的重要性。
4.悬挂系统故障可分为执行器故障、传感器故障和机械部件故障。随着列车在轨运行时间的增长，悬挂系统中的一些部件，例如螺旋弹簧、避震器、空气弹簧、主动/半主动式作动器、传感器都会产生一定程度的性能衰退，诱发诸如螺旋弹簧细微开裂、避震器轻微漏油/漏液、空气弹簧轻度漏气、作动器小幅度作动效能丧失、传感器偏差漂移等早期故障，给列车的行车安全带来潜在的危险。


技术实现要素：

5.本发明的目的是为了解决上述现有技术存在的技术问题和缺陷，提供一种基于数据驱动的高速列车悬挂系统早期故障检测方法。
6.本方法由离线设计和在线检测两部分组成。在离线设计部分中，主要包括利用采集到的无故障数据构造残差信号、建立评价指标体系、以及根据训练数据估计出阈值。在线检测部分利用在线测量值对正交矩阵进行更新，然后根据计算出的评价指标与阈值进行比较，做出故障判定。
7.为达到上述目的，本发明采用以下技术方案予以实现：
8.一种基于数据驱动的高速列车悬挂系统早期故障检测方法，包括如下步骤：
9.步骤一、利用位移传感器、陀螺仪获取高速列车运行时列车车厢质心的垂向位移、角速度信号及转向架构架质心的垂向位移信号，所述传感器及陀螺仪的信号即为系统输出；
10.步骤二、依据高速列车悬挂系统的离散空间状态方程进行数据建模，得到隐含悬挂系统的输入输出关系的悬挂系统输入输出数据模型；步骤一中所得数据按数据模型构造
输入输出数据矩阵和y
k,l
；
11.步骤三、利用步骤二中所得输入输出数据模型，设计基于数据驱动的故障检测残差，并构造评价指标j；
12.步骤四、实时获取步骤一中的系统输出，更新步骤二中的输入输出数据矩阵，然后依据新的输入输出数据矩阵更新步骤三中构造的评价指标；
13.步骤五、实时判断步骤四中更新的评价指标是否超出设定阈值，当评价指标j超过阈值，则给出故障报警：
14.评价指标j的报警阈值为j
th
，由历史数据计算获得，故障检测机制为：
[0015][0016]
进一步优选方案，步骤二中，所述的输入输出数据模型的建模方法，具体步骤和内容如下：
[0017]
步骤1、根据列车悬挂系统运行机理，其离散的系统模型可表示为：
[0018][0019]
其中，a,b,d,c,e,f分别为空间状态方程相应的系数矩阵；x为系统状态，包括车体和转向架的垂向位移以及车体和转向架的俯仰角；u为系统输入，包括主动悬挂中作动器的控制输出量；d为轨道扰动激励；f表征所有可能的故障；y是系统的输出信号，包括车体和转向架的垂向位移以及车体和转向架的俯仰角；w,v分别是过程噪声和测量噪声。
[0020]
步骤2、由步骤1中建立的离散状态空间方程，当系统无故障时，可得如下隐含悬挂系统的输入输出关系的高速列车悬挂系统的数据模型：
[0021][0022]
其中，输入数据矩阵与输出数据矩阵y
k,l
在实际应用中可知。
[0023]
进一步优选方案，步骤三中，所述的基于数据驱动的故障检测残差的设计方法，其设计步骤和内容如下：
[0024]
步骤1、引入正交投影矩阵
[0025]
步骤2、步骤二中建立高速列车悬挂系统的输入输出数据模型两边同时右乘正交投影矩阵由于与相互独立，可得如下基于输入输出数据的残差：
[0026][0027]
进一步优选方案，步骤三中所述的评价指标j的构造方法，步骤如下：
[0028]
步骤1、为消除系统噪音对残差量的影响，引入与噪音独立的辅助变量z
k
，有
[0029][0030]
步骤2、构造对称矩阵
[0031][0032]
步骤3、选取评价指标j＝tr(r)。
[0033]
进一步优选方案，步骤四中，所述的更新输入输出数据矩阵的方法，其更新步骤如下：
[0034]
步骤1、收集在线数据
[0035]
步骤2、通过在数据矩阵最后一列增加新数据列、将删除数据矩阵第一列的方式更新数据矩阵和z
k
：
[0036][0037]
作为本发明的进一步优选方案，步骤四中所述更新评价指标的方法，其步骤为：
[0038]
步骤1、将更新后的输入数据矩阵做如下计算：
[0039][0040]
步骤2、更新正交投影矩阵：
[0041][0042]
步骤3、利用更新后的数据矩阵、正交投影矩阵，依据步骤三中的评价指标构造方法更新评价指标。
[0043]
进一步优选方案，步骤五中，所述的故障检测阈值的计算方法为：
[0044]
j
th
＝max{j
i
}
[0045]
其中j
i
为悬挂系统无故障时，不同时刻下的评价指标。
[0046]
本发明具有以下技术效果：
[0047]
本发明针对高速列车闭环控制结构的悬挂系统，以主动式作动器、传感器、机械部件为故障检测对象，该检测方法对悬挂系统早期故障十分敏感，在悬挂系统发生较小故障的时候就可以实时检测到并准确估计出故障的动态特性。本发明可以用于基于simpack
‑
matlab/simulink仿真架构的执行机构故障分析和系统可靠性分析。
附图说明
[0048]
图1是本发明的故障检测方法流程图。
[0049]
图2是高速列车悬挂系统故障检测结构示意框图。
[0050]
图3是高速列车悬挂系统传感器、陀螺仪安装位置示意图。
[0051]
图4是悬挂执行器早期故障检测仿真曲线图。
[0052]
图5是悬挂传感器早期故障检测仿真曲线图。
[0053]
图6是悬挂机械部件早期故障检测仿真曲线图。
具体实施方式
[0054]
下面结合附图对本发明的技术方案做进一步的详细说明：
[0055]
如附图1
‑
3所示，本发明的一种基于数据驱动的高速列车悬挂系统早期故障检测方法，具体步骤如下：
[0056]
步骤一：建立高速列车悬挂系统的输入输出数据模型，具体步骤和内容是：
[0057]
步骤1.根据列车悬挂系统运行机理，其离散的系统模型可表示为：
[0058][0059]
其中，a,b,d,c,e,f分别为空间状态方程相应的系数矩阵；x为系统状态，包括车体和转向架的垂向位移以及车体和转向架的俯仰角；u为系统输入，包括主动悬挂中作动器的控制输出量；d为轨道扰动激励；f表征所有可能的故障；y是系统的输出信号，包括车体和转向架的垂向位移以及车体和转向架的俯仰角；w,v分别是过程噪声和测量噪声。
[0060]
步骤2.由步骤一中建立的离散状态空间方程，当系统无故障时，有：
[0061][0062]
其中
[0063][0064][0065]
步骤3.步骤二所得表达式含有状态变量x(k)，为消除状态变量x(k)，可引入堆栈矩阵由步骤一中建立的离散状态空间方程可得到
[0066][0067]
步骤4.将步骤三所得表达式代入步骤二所得表达式，有
[0068]
[0069]
化简得
[0070][0071]
其中
[0072][0073][0074][0075]
步骤5.考虑时间间隔n，定义堆栈矩阵步骤5.考虑时间间隔n，定义堆栈矩阵与由步骤四中所得化简后表达式，可得如下隐含悬挂系统的输入输出关系的高速列车悬挂系统的数据模型：
[0076][0077]
其中，输入数据矩阵与输出数据矩阵y
k,l
在实际应用中可知。
[0078]
步骤二：设计基于数据驱动的故障检测残差，具体步骤和内容是：
[0079]
步骤1.引入正交投影矩阵
[0080]
步骤2.步骤五中建立高速列车悬挂系统的输入输出数据模型两边同时右乘正交投影矩阵由于与相互独立，可得如下基于输入输出数据的残差：
[0081][0082]
步骤三：构造评价指标，具体步骤和内容是：
[0083]
步骤1.为消除系统噪音对残差量的影响，引入与噪音独立的辅助变量z
k
，有
[0084][0085]
步骤2.令构造对称矩阵
[0086][0087]
步骤3.选取评价指标j＝tr(r)。
[0088]
步骤四：更新输入输出数据矩阵，具体步骤和内容是：
[0089]
步骤1.收集在线数据
[0090]
步骤2.通过在数据矩阵最后一列增加新数据列、将删除数据矩阵第一列的方式更新数据矩阵和z
k
：
[0091][0092]
步骤五：更新评价指标，具体步骤和内容是：
[0093]
步骤1、将更新后的输入数据矩阵做如下计算：
[0094][0095]
步骤2、更新正交投影矩阵：
[0096][0097]
步骤3、利用更新后的数据矩阵、正交投影矩阵，依据步骤三中的评价指标构造方法更新评价指标。
[0098]
步骤六：故障检测阈值计算，具体步骤和内容是：
[0099]
步骤1.计算n组无故障时的评价指标j
i
,i＝1,2,l n，n为足够大的整数
[0100]
步骤2.计算故障检测阈值：j
th
＝max{tr(r
i
)}
[0101]
步骤七：根据步骤五所得的实时评价指标，和步骤六中的故障检测阈值，对故障进行在线检测，判断检测指标是否达到报警阈值，在检测指标超出阈值时进行报警。
[0102]
当悬挂系统发生故障时，
[0103][0104]
评价指标
[0105][0106]
故障项不为零时，j＞j
th
，即指示故障的发生。
[0107]
下面对本发明的方法进行仿真验证，
[0108]
步骤1.设定执行器故障的特征信息，包括：故障类型为失效1％的执行器早期故障、故障开始时间为第30秒、故障结束时间为仿真结束时间以及故障大小，所述故障通过故障注入模块进行软件注入，建立故障模型；
[0109]
步骤2.设定传感器故障的特征信息，包括：故障类型为传感器发生缓变漂移早期故障，漂移量为θ
s
(t)＝0.001
×
(0.15 0.02sin(t) 0.01sin(0.2t))，故障开始时间为第30秒、故障结束时间为仿真结束时间以及故障大小，所述故障通过故障注入模块进行软件注入，建立故障模型；
[0110]
步骤3.设定机械部件故障的特征信息，包括：拖车二系悬挂系统弹簧失效10％、故障开始时间为第30秒、故障结束时间为仿真结束时间以及故障大小，所述故障通过故障注入模块进行软件注入，建立故障模型；
[0111]
步骤4.采用simpack和matlab/simulink联合仿真，将simpack中建立的轨道
‑
列车耦合的模型导入到matlab/simulink中，并在simulink中建立列车垂向悬挂系统控制仿真模型。设置车辆的运行速度为250km/h，仿真时常为34秒，总共生成1700组样本数据，故障发生时间30s为第1000个采样点对应时刻。
[0112]
如图4所示，当高速列车悬挂系统执行器在第30s时发生失效1％的早期故障，专利提出的方法能够很好的检测出执行器早期故障的发生。
[0113]
如图5所示，当高速列车悬挂系统传感器在第30s时发生缓变漂移早期故障，专利提出的方法能够很好的检测出传感器缓变漂移早期故障的发生。
[0114]
如图6所示，当高速列车悬挂系统机械部件(以二系弹簧为例)在第30s时失效10％的早期故障，专利提出的方法能够很好的检测出列车悬挂系统机械部件早期故障的发生。
[0115]
本发明的故障检测方法对高速列车悬挂系统早期故障具有较高的灵敏度，可以有效的实现高速列车悬挂系统早期故障的检测，有效地解决了闭环控制结构下早期故障的检测及其工程实用问题，这对于高速列车悬挂故障的早期预警和实时监控具有重要的意义。
[0116]
上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。