CN109033582A - 一种高速列车多参数智能阈值准则 - Google Patents
一种高速列车多参数智能阈值准则 Download PDFInfo
- Publication number
- CN109033582A CN109033582A CN201810761712.3A CN201810761712A CN109033582A CN 109033582 A CN109033582 A CN 109033582A CN 201810761712 A CN201810761712 A CN 201810761712A CN 109033582 A CN109033582 A CN 109033582A
- Authority
- CN
- China
- Prior art keywords
- parameter
- monitoring data
- bullet train
- inter
- condition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Emergency Alarm Devices (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种高速列车多参数智能阈值准则,首先采集高速列车的工况参数与状态监测数据,按照工况参数对状态监测数据进行划分;随后利用模糊核聚类对状态监测数据进行分析,得到各类中心坐标、类内距离;最后计算中心坐标、类内距离与工况多参数之间回归关系,进而智能计算各种工况下的阈值,对高速列车进行智能报警。该方法简单易行,适用于变工况下的高速列车监测数据智能阈值的确定。
Description
技术领域
本发明属监测诊断领域,具体涉及一种高速列车多参数智能阈值准则。
背景技术
随着我国铁路的高速发展,高速列车已经成为中国高端制造业崛起的重要标志。然而,高速列车属于典型复杂机电***,以分布式、网络化方式集成了机、电、气、热等多个物理域的部件,部件之间以多种物理作用复杂交互,导致故障表现方式高度复杂化。高速列车的维护普遍沿用不计成本保安全的劳动力密集型计划维修体制,已成为制约我国高铁发展乃至出口的瓶颈,为此,铁路维修保障部门特别提出未来维修方式应该在精确掌握列车状态的前提下,逐渐向状态维修体制过渡,从而保障运行安全、提高维修效率,满足国内和海外维修保障需求。
在列车的实时状态监测与故障诊断过程中,如何制定阈值是行业内的一大难题。主要原因在于:高速通常在高转速、强振动、大应力的恶劣环境下工作,单一的阈值指标难以适应工况的变化。如果阈值过高,只能在最大转速的极限情况下报警,会产生较大的漏报率,起不到实时预警、报警的作用。如果阈值过低,虽能适应低转速工况,一旦高转速工作,正常状态误认为故障,会导致很大的维修资源浪费,并带来人员恐慌,影响正常工作。因此,高速列车监测数据的阈值确定需从实时监测数据出发,并适应各种工况的变化。
发明内容
本发明的目的在于,提供一种高速列车多参数智能阈值准则,从监测数据出发,融入工况多参数,智能确定多种工况下的监测阈值,实现高速列车自适应工况的实时报警。
为了实现上述目的,本发明采取的技术方案是:
1)采集高速列车的工况参数与状态监测数据,按照工况参数对状态监测数据进行划分;
2)利用模糊核聚类对状态监测数据进行分析,得到各类中心坐标、类内距离;
3)计算中心坐标、类内距离与工况多参数之间回归关系,进而智能计算各种工况下的阈值,对高速列车进行智能报警。
所述的步骤1),具体包含以下内容:
首先,采集各种工况下的状态监测数据,其对应的工况参数一并采集,监测状态数据主要为振动加速度信号,工况参数主要包括转速、载荷、温度等。
随后,计算状态监测数据的统计特征:如均方根值、方根幅值、偏斜度指标、峭度指标、频域幅值特征、频率集散特征等,并按照工况参数对状态监测数据进行简单划分。
所述的步骤2),具体包含以下内容:
首先,利用权利要求2中求得的监测数据统计特征构建N×d的样本集X={x1,x2,…,xN},其中N为监测数据时间序列的个数;d为统计特征的种类数;xi,i=1,2,…,N是维数为d的样本。按照工况参数可将样本集分成J类,定义隶属度μij表示分类结果,其含义为样本xi属于第j类的可能性,j=1,2,…,J。
随后,给出模糊核聚类分析的目标函数:
其中U=(μij)J×N为隶属度矩阵,V={v1,v2,…,vJ},vj为第j类的中心坐标,m≥1为加权指数。
为解决原始样本非线性的问题,引入非线性函数φ(x)把原始样本映射至高维空间中,构建核函数K(xj,vi)=φ(xj)T·φ(xj)代替xi·vj,得到下式
对上式进行求解,得到隶属度μij和类中心坐标vj
基于上述公式,模糊核聚类算法的步骤为:
①设定聚类数J、模糊聚类指数m,以及停止阈值ε;
②根据以上参数初始化隶属度矩阵U;
③计算目标函数,判断是否达到中止条件,如果满足,聚类结束,否则继续。
④根据式(3)和式(4)更新模糊聚类中心和模糊隶属度矩阵,并转向步骤③。
最后,第j类聚类的类内距离为
dj=max||xl-vj||,l=1,2,…,L (11)
所述的步骤3),具体包含以下内容:
首先,将类心坐标V={v1,v2,…,vJ}、类内距离D={d1,d2,…,dJ}与工况多参数E={e1,e2,…,es}构建多元非线性回归模型:
其中ek为J维向量,包含各种工况。利用数据自动拟合求解函数f1、f2,得到类心坐标V、类内距离D与工况多参数E之间的函数关系,继而可得到各种工况对应类内距离D之间的关系。
然后规定额定工况J0条件下,类内距离超出正常其50%视为报警门限,然后计算出各种工况对应的报警门限d'j,f1、f2的逆函数求得,进而实现对高速列车的智能报警。
由于本发明一种高速列车多参数智能阈值准则,具有下列区别于传统方法的显著优势:
1)基于模糊核聚类分析方法实现了监测数据的信息融合,并对各工况下的监测数据进行了自适应聚类,将各工况下的监测数据在同一个维度下展开,便于制定监测阈值。
2)将工况多参数融入阈值的确定过程中,实现阈值随工况的自动更新,智能识别高速列车故障。
3)计算简单,运算速度快,便于在工程实际中推广。
附图说明
图1所示为高速列车轴承信号原始特征;
图2所示为高速列车轴承信号模糊核聚类结果;
图3所示为高速列车轴承信号到类心的距离;
具体实施方式
参照图1所示的高速列车轴承信号原始特征,在模拟试验台获取的振动加速度信号计算获得,横坐标为时间单位min,由于各个统计特征含义不同,纵坐标无单位。各个统计特征差别很大,很难再一个尺度下展示,且随工况的影响很大。
参考图2所示的高速列车轴承信号模糊核聚类结果,自适应按照工况对特征进行分类,为在各种工况下制定阈值奠定基础,横坐标、纵坐标均为合成参数,没有明显的物理意义。
参考图3所示的高速列车轴承信号到类心的距离,摆脱的工况的影响,在同一个维度下展开信号特征。横坐标为时间单位min,纵坐标是距离。
附图是本发明的具体实施例;
下面结合附图对本发明的内容作进一步详细说明:
1)采集高速列车的工况参数与状态监测数据,按照工况参数对状态监测数据进行划分;
2)利用模糊核聚类对状态监测数据进行分析,得到各类中心坐标、类内距离;
3)计算中心坐标、类内距离与工况多参数之间回归关系,进而智能计算各种工况下的阈值,对高速列车进行智能报警。
所述的步骤1),具体包含以下内容:
首先,采集各种工况下的状态监测数据,其对应的工况参数一并采集,监测状态数据主要为振动加速度信号,工况参数主要包括转速、载荷、温度等。
随后,计算状态监测数据的统计特征:如均方根值、方根幅值、偏斜度指标、峭度指标、频域幅值特征、频率集散特征等,并按照工况参数对状态监测数据进行简单划分。
表1.统计特征
所述的步骤2),具体包含以下内容:
首先,利用权利要求2中求得的监测数据统计特征构建N×d的样本集X={x1,x2,…,xN},其中N为监测数据时间序列的个数;d为统计特征的种类数;xi=(T1,T2,…,T6),i=1,2,…,N是维数为d=6的样本。按照工况参数可将样本集分成J类,定义隶属度μij表示分类结果,其含义为样本xi属于第j类的可能性,j=1,2,…,J。
随后,给出模糊核聚类分析的目标函数:
其中U=(μij)J×N为隶属度矩阵,V={v1,v2,…,vJ},vj为第j类的中心坐标,m≥1为加权指数。
为解决原始样本非线性的问题,引入非线性函数φ(x)把原始样本映射至高维空间中,构建核函数K(xj,vi)=φ(xj)T·φ(xj)代替xi·vj,得到下式
对上式进行求解,得到隶属度μij和类中心坐标vj
基于上述公式,模糊核聚类算法的步骤为:
①设定聚类数J、模糊聚类指数m,以及停止阈值ε;
②根据以上参数初始化隶属度矩阵U;
③计算目标函数,判断是否达到中止条件,如果满足,聚类结束,否则继续。
④根据式(3)和式(4)更新模糊聚类中心和模糊隶属度矩阵,并转向步骤③。
第j类聚类的类内距离为
dj=max||xl-vj||,l=1,2,…,L (17)
所述的步骤3),具体包含以下内容:
首先,将类心坐标V={v1,v2,…,vJ}、类内距离D={d1,d2,…,dJ}与工况多参数E={e1,e2,...,es}构建多元非线性回归模型:
其中ek为J维向量,包含各种工况。利用数据自动拟合求解函数f1、f2,得到类心坐标V、类内距离D与工况多参数E之间的函数关系,继而可得到各种工况对应类内距离D之间的关系。
然后规定额定工况J0条件下,类内距离超出正常其50%视为报警门限,然后计算出各种工况对应的报警门限d'j,f1、f2的逆函数求得,进而实现对高速列车的智能报警。
实施例:
该实施例给出了本发明在高速列车轴承试验中的具体实施过程,同时验证了该发明的有效性。试验按照四种转速、两种载荷分成以下8中工况,如表2所示,其中工况6为额定工况。
表2高速列车工况
采集高速列车轴承振动加速度计算统计特征值,如图1所示,主要分为两部分,频域幅值指标数据较大,在图形的上部展现,其余指标均较小,在图形下部展现。为进一步显示下部指标的趋势,在左下部出现了一个局部放大图,可以更好地看到均方根值、方根幅值及峭度指标随工况的变化情况。
利用模糊核聚类对统计特征进行信息融合,可以按照工况进行聚类,在高维特征空间中展开,结果如图2所示。计算各个样本到类心的距离,如图3所示。故各种工况下类内距离即为图3所示的各种工况下的最大值,如表3所示。
表3各工况类内距离
工况1 | 工况2 | 工况3 | 工况4 | 工况5 | 工况6 | 工况7 | 工况8 |
5.08 | 7.96 | 11.99 | 14.05 | 9.22 | 13.01 | 8.75 | 12.56 |
其中,工况6为额定工况,在额定工况下,其类内距离超出50%视为报警门限,通过求解类心坐标V、类内距离D与工况多参数E的多元非线性回归模型,最终计算出各种工况下的报警内类距离如表4所示,实现对高速列车的智能报警。
表4各工况报警类内距离
工况1 | 工况2 | 工况3 | 工况4 | 工况5 | 工况6 | 工况7 | 工况8 |
7.36 | 11.8 | 16.75 | 21.22 | 12.18 | 19.65 | 12.21 | 18.83 |
Claims (4)
1.一种高速列车多参数智能阈值准则,其特征在于:
1)采集高速列车的工况参数与状态监测数据,按照工况参数对状态监测数据进行划分;
2)利用模糊核聚类对状态监测数据进行分析,得到各类中心坐标、类内距离;
3)计算中心坐标、类内距离与工况多参数之间回归关系,进而智能计算各种工况下的阈值,对高速列车进行智能报警。
2.根据权利要求1所述的一种高速列车多参数智能阈值准则,其特征在于,所说的采集高速列车的工况参数与状态监测数据,按照工况参数对状态监测数据进行划分,包括以下步骤:
首先,采集各种工况下的状态监测数据,其对应的工况参数一并采集,监测状态数据主要为振动加速度信号,工况参数主要包括转速、载荷、温度等。
随后,计算状态监测数据的统计特征:如均方根值、方根幅值、偏斜度指标、峭度指标、频域幅值特征、频率集散特征等,并按照工况参数对状态监测数据进行简单划分。
3.根据权利要求1所述的一种高速列车多参数智能阈值准则,其特征在于,所说的利用模糊核聚类对状态监测数据进行分析,得到各类中心坐标、类内距离,包括以下步骤:
首先,利用权利要求2中求得的监测数据统计特征构建N×d的样本集X={x1,x2,…,xN},其中N为监测数据时间序列的个数;d为统计特征的种类数;xi,i=1,2,…,N是维数为d的样本。按照工况参数可将样本集分成J类,定义隶属度μij表示分类结果,其含义为样本xi属于第j类的可能性,j=1,2,…,J。
随后,给出模糊核聚类分析的目标函数:
其中U=(μij)J×N为隶属度矩阵,V={v1,v2,…,vJ},vj为第j类的中心坐标,m≥1为加权指数。
为解决原始样本非线性的问题,引入非线性函数φ(x)把原始样本映射至高维空间中,构建核函数K(xj,vi)=φ(xj)T·φ(xj)代替xi·vj,得到下式
对上式进行求解,得到隶属度μij和类中心坐标vj
基于上述公式,模糊核聚类算法的步骤为:
①设定聚类数J、模糊聚类指数m,以及停止阈值ε;
②根据以上参数初始化隶属度矩阵U;
③计算目标函数,判断是否达到中止条件,如果满足,聚类结束,否则继续。
④根据式(3)和式(4)更新模糊聚类中心和模糊隶属度矩阵,并转向步骤③。
最后,第j类聚类的类内距离为
dj=max||xl-vj||,l=1,2,…,L (5) 。
4.根据权利要求1所述的一种高速列车多参数智能阈值准则,其特征在于,所说的计算中心坐标、类内距离与工况多参数之间回归关系,进而智能计算各种工况下的阈值,对高速列车进行智能报警,包括以下步骤:
首先,将类心坐标V={v1,v2,…,vJ}、类内距离D={d1,d2,…,dJ}与工况多参数E={e1,e2,…,es}构建多元非线性回归模型:
其中ek为J维向量,包含各种工况。利用数据自动拟合求解函数f1、f2,得到类心坐标V、类内距离D与工况多参数E之间的函数关系,继而可得到各种工况对应类内距离D之间的关系。
然后规定额定工况J0条件下,类内距离超出正常其50%视为报警门限,然后计算出各种工况对应的报警门限d'j,f1、f2的逆函数求得,进而实现对高速列车的智能报警。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810761712.3A CN109033582A (zh) | 2018-07-12 | 2018-07-12 | 一种高速列车多参数智能阈值准则 |
CN201910510051.1A CN110781553B (zh) | 2018-07-12 | 2019-06-12 | 一种高速列车多参数智能阈值监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810761712.3A CN109033582A (zh) | 2018-07-12 | 2018-07-12 | 一种高速列车多参数智能阈值准则 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109033582A true CN109033582A (zh) | 2018-12-18 |
Family
ID=64642350
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810761712.3A Pending CN109033582A (zh) | 2018-07-12 | 2018-07-12 | 一种高速列车多参数智能阈值准则 |
CN201910510051.1A Active CN110781553B (zh) | 2018-07-12 | 2019-06-12 | 一种高速列车多参数智能阈值监测方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910510051.1A Active CN110781553B (zh) | 2018-07-12 | 2019-06-12 | 一种高速列车多参数智能阈值监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN109033582A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539374A (zh) * | 2020-05-07 | 2020-08-14 | 上海工程技术大学 | 基于多维数据空间的轨道列车轴承故障诊断***及方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113297291A (zh) * | 2021-05-08 | 2021-08-24 | 上海电气风电集团股份有限公司 | 监测方法、***、可读存储介质和风力发电机 |
CN114323707B (zh) * | 2022-01-04 | 2023-07-07 | 中车株洲电力机车有限公司 | 磁浮列车及其振动信号计算方法、模拟发生方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0318339D0 (en) * | 2003-08-05 | 2003-09-10 | Oxford Biosignals Ltd | Installation condition monitoring system |
CN102636991A (zh) * | 2012-04-18 | 2012-08-15 | 国电科学技术研究院 | 一种基于模糊集关联规则的火电机组运行参数优化方法 |
CN103455635A (zh) * | 2013-09-24 | 2013-12-18 | 华北电力大学 | 基于最小二乘支持向量机集成的热工过程软测量建模方法 |
CN104392071B (zh) * | 2014-12-12 | 2017-09-29 | 北京交通大学 | 一种基于复杂网络的高速列车***安全评估方法 |
CN106203856A (zh) * | 2016-07-18 | 2016-12-07 | 交通运输部公路科学研究所 | 一种组合主成分分析和模糊c均值聚类的车辆行驶工况制定方法 |
-
2018
- 2018-07-12 CN CN201810761712.3A patent/CN109033582A/zh active Pending
-
2019
- 2019-06-12 CN CN201910510051.1A patent/CN110781553B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539374A (zh) * | 2020-05-07 | 2020-08-14 | 上海工程技术大学 | 基于多维数据空间的轨道列车轴承故障诊断***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110781553B (zh) | 2023-08-25 |
CN110781553A (zh) | 2020-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106124175B (zh) | 一种基于贝叶斯网络的压缩机气阀故障诊断方法 | |
CN109033582A (zh) | 一种高速列车多参数智能阈值准则 | |
CN110110768B (zh) | 基于并行特征学习和多分类器的滚动轴承故障诊断方法 | |
CN109934206A (zh) | 一种非平稳工况下的旋转机械故障诊断方法 | |
CN104832418B (zh) | 一种基于局部均值变换和Softmax的液压泵故障诊断方法 | |
CN104596780B (zh) | 一种动车组制动***传感器故障的诊断方法 | |
CN112819059B (zh) | 一种基于流行保持迁移学习的滚动轴承故障诊断方法 | |
CN101718634A (zh) | 基于多元概率模型的设备状态综合动态报警方法 | |
CN109708907B (zh) | 一种基于包络信息的设备故障特征提取方法 | |
CN111539374B (zh) | 基于多维数据空间的轨道列车轴承故障诊断方法 | |
CN113542241B (zh) | 一种基于CNN-BiGRU混合模型的入侵检测方法及装置 | |
CN105894024A (zh) | 多重核的可能性模糊聚类算法 | |
CN103291600A (zh) | 一种基于emd-ar和mts的液压泵故障诊断方法 | |
CN101178703A (zh) | 基于网络分割的故障诊断谱聚类方法 | |
CN105225523A (zh) | 一种车位状态检测方法及装置 | |
CN108267312A (zh) | 一种基于快速搜索算法的地铁列车轴承智能诊断方法 | |
CN102156873A (zh) | 一种基于混沌的机械零部件早期单点故障检测与分类方法 | |
CN104622446A (zh) | 一种基于khm聚类算法的心率变异信号优化方法 | |
CN109933040A (zh) | 基于层次密度峰值聚类和最相似模态的故障监测方法 | |
CN105956318B (zh) | 基于改进***h-k聚类方法的风电场机群划分方法 | |
CN114548295A (zh) | 基于多尺度领域自适应网络的轴承故障分类***及方法 | |
CN104239411B (zh) | 一种基于颜色、位置聚类和角点检测的栅格状雷达的检测方法 | |
CN107392979B (zh) | 时间序列的二维可视状态构图及定量分析指标方法 | |
CN105590167A (zh) | 电场多元运行数据分析方法及装置 | |
Khalid et al. | Brain abnormalities segmentation performances contrasting: adaptive network-based fuzzy inference system (ANFIS) vs K-nearest neighbors (k-NN) vs fuzzy c-means (FCM) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181218 |