一种车轮多边形状态识别方法和系统技术方案

本发明专利技术公开一种车轮多边形状态识别方法和系统，包括S1：获取车轮的环境信号，并对环境信号进行处理得到角域信号；S2：对角域信号进行FFT得到阶次域振动信号，再根据阶次域振动信号构建非线性双稳态系统，增强第k阶信号，从而提取因车轮多边形激励的加强振动信号频率；S3：将加强振动信号频率范围内的阶次域振动信号截取后进行傅里叶反变换得到振动加速度阶次域信号，再对一个周期内的振动加速度阶次域信号进行二次积分得到周期位移时域信号，进而提取车轮的波长；S4：将S3中车轮的波长与预设最大波长进行比较，当车轮的波长不小于预设最大波长时，提出预警、报警信息。报警信息。报警信息。

【技术实现步骤摘要】
一种车轮多边形状态识别方法和系统


[0001]本专利技术涉及轨道交通
，特别涉及一种车轮多边形状态识别方法和系统。

技术介绍

[0002]城市轨道交通中的地铁车辆在运行过程中，车轮是链接列车与路轨、承载列车重量以及受力的最重要构件，因此车轮和路轨的接触状态与列车行驶是否平稳、乘客是否舒适、安全息息相关。
[0003]车轮在运行过程中往往会因为侵蚀、轮轴共振、磨损和行车特性等原因产生不同形式的磨损，而其中表现车轮不圆顺的主要特征为车轮多边形现象。而车轮上不同阶数的多边形会对车辆行驶造成不同程度的影响，因此建立一种有效的方法对车轮多边形情况进行分析具有很高的重要性。
[0004]轨道交通中用于检测车轮磨损的方法主要分为两种。一是轨旁式动态检测法，其主要是利用各种类型的传感器对路轨状态进行监测，主要采用轮轨力检测法、激光与机器视觉检测法和声面状态，但抗干扰能力差，精度低，不能很好反映车轮多边形情况而且成本高，技术难度大很难成规模进行推广；另外一种车轮磨损检测方法为车载式动态检测法，通过将传感器直接安装在车辆的轴箱和转向架上进行车轮的不间断监测，实时获取车辆行驶数据进而对车轮故障进行定量定位分析，但是车载式动态监测法采集通道多，数据量大，实用性仍有不足，并且安装在轴箱处会受到轴承振动以及其他结构振动的影响，导致识别精度较低。

技术实现思路

[0005]针对现有技术中车轮多边形状态识别精度低的技术问题，本专利技术提出一一种车轮多边形状态识别方法和系统，通过对地铁车辆轴箱处的振动状况进行信号采集和故障监测，从而判别车轮多边形所处状态，提高了识别精度。
[0006]为了实现上述目的，本专利技术提供以下技术方案：
[0007]一种车轮多边形状态识别方法，包括以下步骤：
[0008]S1：通过增量型编码器获取车轮的角域信号；
[0009]S2：对角域信号进行FFT得到阶次域振动信号，再根据阶次域振动信号构建非线性双稳态系统，增强第k阶信号，从而提取因车轮多边形激励的加强振动信号频率；
[0010]S3：将加强振动信号频率范围内的阶次域振动信号截取后进行傅里叶反变换得到振动加速度阶次域信号，再对一个周期内的振动加速度阶次域信号进行二次积分得到周期位移时域信号，进而提取车轮在周期内的最大波长；
[0011]S4：将S3中车轮在周期内的最大波长与预设最大波长进行比较，当车轮在周期内的最大波长不小于预设最大波长时，车轮处于故障状态，提出预警、报警信息。
[0012]优选地，所述S1中，环境信号包括振动信号、轴箱机械构件所产生的有色噪声以及环境中的白噪声。
[0013]优选地，所述S1中，对环境信号进行处理的方式为：
[0014]在车辆轴箱处安装增量型编码器，车轴旋转时增量型编码器的光栅发生变化，从而产生采样脉冲，再通过采样脉冲触发采集信号；由于采样脉冲频率对应增量型编码器分辨率且由于增量型编码器分辨率转角确定，从而得到对应的角域信号。
[0015]优选地，所述S2包括：
[0016]S2
‑
1：对角域信号进行FFT可以得到阶次域振动信号：
[0017][0018]公式(2)中，S(k)为阶次域振动信号；s(θ
n
)为角域信号；N表示信号长度；k表示阶次变量；表示复指数常量；n表示阶次域振动信号的阶次；j表示复数常量；
[0019]S2
‑
2：构建非线性双稳态系统，对振动信号进行增强，从而提取多边形激励的振动信号频率；
[0020]非线性双稳态系统的Langevin方程表示为：
[0021][0022]公式(3)中，a、b表示状态参数；x为非线性双稳态系统变量；k为阶次变量；S(k)表示阶次域振动信号；Γ(k)表示阶次域噪声信号，D表示噪声强度，ξ(k)为零均值的高斯白噪声，满足数学期望E(Γ(k))＝0，为零均值的高斯白噪声，满足数学期望E(Γ(k))＝0，表示相位移动。
[0023]优选地，所述S2
‑
2中，a、b的求取方法为：
[0024]A、设阶次域振动信号的阶次为1
‑
12，计算得到阶次域振动信号的阶次功率：
[0025][0026]公式(5)中，P
s
表示阶次域振动信号的阶次功率；n表示阶次域振动信号的阶次；k表示阶次变量；S(k)用参数a、b表达，S(k)表示阶次域振动信号；Γ(k)表示噪声信号；
[0027]B、根据阶次功率得到信噪比：
[0028]Snr＝10log(Ps/Pn)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0029]公式(5)中，Snr表示信噪比；P
s
表示阶次域振动信号的阶次功率；Pn表示整个阶次域中信号功率，计算公式为：
[0030]C、令a、b∈[k1，k2]，K1、K2分别表示取值下限、上限，a和b均以Δ＝(k2
‑
k1)/n为梯度计算S(K)，并计算对应Snr
(i,j)
值，从而生成信噪比粒子群矩阵，并给出对应初始速度矩阵其中V
(i,j)
＝1：
[0031][0032]公式(6)中，POP表示信噪比粒子群矩阵；Snr
(i,j)
表示通过第i个a，第j个b计算得到的信号信噪比，V
(i,j)
表示对应信噪比粒子的初始速度；
[0033]D、设迭代次数为m，适应度阈值为ω，最大允许迭代步数为t，适度函数选择为p
i
＝max{Snr
(i,j)
|Snr
(i,j)
∈POP}，则迭代速度计算公式为：
[0034][0035]公式(7)中，表示迭代速度；ω
i
表示适应度阈值，表示粒子移动速度，c1、c2表示两个正常数，rand1，rand2表示两个独立的、介于[0，1]之间的随机数，p
besti
表示优化解；表示第i个粒子，g
besti
表示全局最好适应度函数值，i表示粒子序号，d表示维度参数，m表示迭代次数；
[0036]E、计算每个粒子的适应函数值p
i
，并与优化解p
besti
进行比较：若p
i
<p
besti
，用新的粒子替代前一轮粒子，即令p
besti
＝Snr
i
，再与全局最好适应度函数值g
besti
进行比较，若p
i
≥p
besti
继续寻找下一个粒子进行比较；若p
i
<g
besti
则取代每个粒子，p
besti
＝g
besti
＝Snr
i
，并保存当前状态进行下本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种车轮多边形状态识别方法，其特征在于，包括以下步骤：S1：通过增量型编码器获取车轮的角域信号；S2：对角域信号进行FFT得到阶次域振动信号，再根据阶次域振动信号构建非线性双稳态系统，增强第k阶信号，从而提取因车轮多边形激励的加强振动信号频率；S3：将加强振动信号频率范围内的阶次域振动信号截取后进行傅里叶反变换得到振动加速度阶次域信号，再对一个周期内的振动加速度阶次域信号进行二次积分得到周期位移时域信号，进而提取车轮在周期内的最大波长；S4：将S3中车轮在周期内的最大波长与预设最大波长进行比较，当车轮在周期内的最大波长不小于预设最大波长时，车轮处于故障状态，提出预警、报警信息。2.如权利要求1所述的一种车轮多边形状态识别方法，其特征在于，所述S1中，环境信号包括振动信号、轴箱机械构件所产生的有色噪声以及环境中的白噪声。3.如权利要求1所述的一种车轮多边形状态识别方法，其特征在于，所述S1中，对环境信号进行处理的方式为：在车辆轴箱处安装增量型编码器，车轴旋转时增量型编码器的光栅发生变化，从而产生采样脉冲，再通过采样脉冲触发采集信号；由于采样脉冲频率对应增量型编码器分辨率且由于增量型编码器分辨率转角确定，从而得到对应的角域信号。4.如权利要求1所述的一种车轮多边形状态识别方法，其特征在于，所述S2包括：S2
‑
1：对角域信号进行FFT可以得到阶次域振动信号：公式(2)中，S(k)为阶次域振动信号；s(θ
n
)为角域信号；N表示信号长度；k表示阶次变量；表示复指数常量；n表示阶次域振动信号的阶次；j表示复数常量；S2
‑
2：构建非线性双稳态系统，对振动信号进行增强，从而提取多边形激励的振动信号频率；非线性双稳态系统的Langevin方程表示为：公式(3)中，a、b表示状态参数；x为非线性双稳态系统变量；k为阶次变量；S(k)表示阶次域振动信号；Γ(k)表示阶次域噪声信号，D表示噪声强度，ξ(k)为零均值的高斯白噪声，满足数学期望E(Γ(k))＝0，均值的高斯白噪声，满足数学期望E(Γ(k))＝0，表示相位移动。5.如权利要求4所述的一种车轮多边形状态识别方法，其特征在于，所述S2
‑
2中，a、b的求取方法为：A、设阶次域振动信号的阶次为1
‑
12，计算得到阶次域振动信号的阶次功率：
公式(5)中，P
s
表示阶次域振动信号的阶次功率；n表示阶次域振动信号的阶次；k表示阶次变量；S(k)用参数a、b表达，S(k)表示阶次域振动信号；Γ(k)表示噪声信号；B、根据阶次功率得到信噪比：Snr＝10log(Ps/Pn)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)公式(5)中，Snr表示信噪比；P
s
表示阶次域振动信号的阶次功率；Pn表示整个阶次域中信号功率，计算公式为：C、令a、b∈[k1，k2]，K1、K2分别表示取值下限、上限，a和b均以Δ＝(k2
‑
k1)/n为梯度计算S(K)，并计算对应Snr
(i,j)
值，从而生成信噪比粒子群矩阵，并给出对应初始速度矩阵其中V
(i,j)
＝1：公式(6)中，POP表示信噪比粒子群矩阵；Snr
(i,j)
表示通过第i个a，第j个b计算得到的信号信噪比，V
(i,j)
表示对应信噪比粒子的初始速度；D、设迭代次数为m，适应度阈值为ω，最大允许迭代步数为t，适度函数选择为pi＝max{Snr
(i,j)
||Snr
(i,j)
∈POP}，则迭代速度计算公式为：公式(7)中，表示迭代速度；ω
i
表示适应度阈值，表示粒子...

【专利技术属性】
技术研发人员：王利明，赵林华，孙辰，李林，黄文彬，丁晓喜，
申请(专利权)人：重庆大学，
类型：发明
国别省市：