CN109186971B

CN109186971B - 基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法

Info

Publication number: CN109186971B
Application number: CN201810884359.8A
Authority: CN
Inventors: 薛红涛; 陈震宇; 江洪
Original assignee: Jiangsu University
Current assignee: Jiangyin Intellectual Property Operation Co., Ltd
Priority date: 2018-08-06
Filing date: 2018-08-06
Publication date: 2020-05-05
Anticipated expiration: 2038-08-06
Also published as: CN109186971A

Abstract

本发明公开一种基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，分别计算每段加速度信号时域和频域高敏感特征参数，构建当前时间片段k对应的速度片v_k的贝叶斯网络结构，得到各网络节点的条件概率分布，确定前一时间片段k‑1对应的速度片v_k‑1与当前时间片段k对应的速度片v_k之间的二速度片v_k‑1‑v_k状态转移概率分布，建立动态贝叶斯网络模型，根据速度片v_k‑1和速度片v_k之间的多种状态转移概率分布建立故障诊断模型群，在线采集轮毂电机运行信息，从故障诊断模型群中选择相应的二速度片v_k‑1‑v_k诊断模型，计算后验概率分布，判断轮毂电机正常或故障，提高了识别轮毂电机机械故障的准确性与时效性。

Description

基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法

技术领域

本发明涉及汽车轮毂电机机械故障的状态监测与智能诊断领域，特别是一种基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法。

背景技术

采用轮毂电机驱动的电动汽车具有结构简单紧凑、传动效率高等突出优点。轮毂电机是电动汽车驱动系统的核心，由于轮毂电机安装在狭小的轮毂空间内，磁场饱和、转矩波动、负载突变等因素对其性能影响明显，且多变的汽车行驶工况和复杂的路况极易诱发轮毂电机机械故障，从而引起振动增强、效率下降、温升加剧，持续长时间运行还会导致绝缘材料性能下降、定转子摩擦，进而引发绕组受损、匝间断路、相间短路等次生故障，危及车辆运行安全，因此有必要对轮毂电机运行状态进行在线监测和诊断，提高识别轮毂电机机械故障的准确性与时效性。

轮毂电机在运行过程中产生的振动及其特征信息反映了其运行状态的变化，因此有效地分析振动信号、提炼出对故障敏感度高的特征参数是进行精确诊断的关键。目前采用较多的故障信号提取方法包括：

(1)2013年Li K等在杂志《Sensors》(第13卷第6期)的“Sequential fuzzydiagnosis method for motor roller bearing in variable operating conditionsbased on vibration analysis”论文中提到了一种区分度指标法(Distinguish index,DI)，它由DI值和两状态区分率(Discrimination Rate,DR)组成，但是该方法仅对两种状态下单一高敏感特征参数快捷有效，无法同时选取多个高敏感特征参数。

(2)2012年Li K等在杂志《Sensors》(第12卷第5期)的“An IntelligentDiagnosis Method for Rotating Machinery Using Least Squares Mapping and aFuzzy Neural Network”论文中公布了一种综合区分度指标法(Synthetic detectionindex,SDI)，用于提取多个敏感特征参数，它由多个状态下的DI值组成，但是该方法忽略了通常情况下机械设备均有多种工况，不能够选取对多种工况均高度敏感的特征参数。

贝叶斯网络(Bayesian Networks,BNs)使用了一种方便的框架结构来表示因果关系，使不确定性推理在逻辑上变得更为清晰。已有不少将贝叶斯网络应用于故障诊断技术中，具体包括：

(1)2011年Li Z等在杂志《Procedia Engineering》(第16卷第16期)的“FaultDiagnosis of Motor Bearing Based on the Bayesian Network”论文中公布了一种将贝叶斯网络应用于电机轴承故障的诊断方法，它基于振动信号和贝叶斯诊断模型进行诊断，但是该方法没有把前一时间段对当前时间段的状态影响考虑进去，不能够动态地调整诊断结果。

(2)2016年赵月南等在杂志《机电工程》(第33卷第2期)的“采用布谷鸟算法的贝叶斯网络在异步机故障诊断中的应用”论文中公布了一种采用布谷鸟算法的贝叶斯网络，并将其应用于异步电机的故障诊断中，它基于异步电机故障电流信号和贝叶斯诊断模型进行诊断，但是由于轮毂电机运行环境复杂，该方法对轮毂电机故障状态的识别比较低，并且该方法是基于离线数据进行的故障诊断，不能够实现快速有效的在线识别诊断，无法满足轮毂电机的运行安全要求。

发明内容

针对上述问题，本发明提供一种基于动态贝叶斯网络(Dynamic BayesianNetworks,DBNs)的轮毂电机机械故障在线诊断方法，考虑到相邻时间片之间的状态影响，因此能动态地调整诊断结果，并且能在线识别诊断，快速有效地识别诊断。

本发明所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法采用的技术方案是具有以下步骤：

步骤一：采集轮毂电机在正常和故障运行状态下的运行信息，包括车速信号、加速度信号，将每种加速度信号等分为15段，分别计算每段加速度信号时域和频域高敏感特征参数；

步骤二：根据运行状态、加速度、高敏感特征参数建立训练数据集，将运行状态、加速度作为隐藏节点，高敏感特征参数作为可观测的高斯节点，对训练数据集进行参数学习，构建当前时间片段k对应的速度片v_k的贝叶斯网络结构，得到各网络节点的条件概率分布，确定前一时间片段k-1对应的速度片v_k-1与当前时间片段k对应的速度片v_k之间的二速度片v_k-1-v_k状态转移概率分布；

步骤三：结合当前时间片段k对应的速度片v_k的贝叶斯网络结构和二速度片v_k-1-v_k状态转移概率分布，建立动态贝叶斯网络模型；

步骤四：根据速度片v_k-1和速度片v_k之间的多种状态转移概率分布建立二速度片v_k-1-v_k轮毂电机机械故障诊断模型群；

步骤五：在线采集轮毂电机的车速V_k、加速度a，计算出当前时间片段k中加速度信号a的高敏感特征参数；

步骤六：根据前一时间片段k-1的车速V_k-1和当前时间片段k的车速V_k从所述的轮毂电机故障诊断模型群中选择相应的二速度片v_k-1-v_k诊断模型。

步骤七：将步骤五中所述的高敏感特征参数作为选择的二速度片v_k-1-v_k诊断模型的输入，计算轮毂电机运行状态的后验概率分布；

步骤八：根据所述的后验概率判断轮毂电机正常或故障。

本发明的有益效果是：

1、本发明通过提出的综合权重诊断指标法(Synthetic weight detectionindex，SWDI)可以提炼出对多种工况下多种状态均具有高敏感度的多个特征参数，更有利于所建模型的推理，提高诊断结果的精度。

2、本发明提出的速度片，解决轮毂电机运行状态在相邻时间片之间无法构建转移概率分布的问题，从而建立轮毂电机故障诊断模型群。

3、本发明基于DBNs的轮毂电机故障诊断模型群能够实现轮毂电机机械故障的在线诊断，提高了识别轮毂电机机械故障的准确性与时效性。

4、本发明能有效降低故障诊断错误率，提高轮毂电机运行的安全性。

附图说明

图1是本发明在第一阶段轮毂电机机械故障诊断模型群的建立流程图：

图2是图1中获得加速度信号的时域、频域高敏感特征参数的流程图；

图3是图1中当前时间片段k对应的速度片v_k内的轮毂电机贝叶斯网络诊断模型；

图4是图1中基于二速度片的动态贝叶斯网络的轮毂电机机械故障诊断模型；

图5是本发明中在第二阶段基于诊断模型群的轮毂电机机械故障在线诊采用的模块结构图；

图6是本发明在第二阶段实施轮毂电机机械故障在线诊断方法的流程图。

具体实施方式

下面结合附图和具体实施方式对本发明技术方案作进一步的详细描述。

本发明包括两个阶段，第一阶段为基于离线数据的轮毂电机机械故障诊断模型群的建立阶段；第二阶段是基于诊断模型群的机械故障在线诊断阶段。

将轮毂电机运行状态的改变过程理解为一系列随速度改变的快照，每个快照描述了轮毂电机在相应时间片段内特定速度下的状态，这样的快照被称为“速度片”，则可先构建单个“速度片”内的贝叶斯网络模型，再通过确定不同“速度片”之间的状态转移概率分布构建“二速度片”动态贝叶斯网络，即构建前一时间片段k-1对应的速度片v_k-1与当前时间片段k对应的速度片v_k的轮毂电机机械故障诊断模型群。如图1所示的第一阶段轮毂电机机械故障诊断模型群的建立流程图，其具体步骤如下：

步骤一：收集并整理轮毂电机在正常和故障状态下的运行信息作为建立故障诊断模型群的离线数据。具体包括在车速5km/h，10km/h，15km/h，20km/h和25km/h下的正常和故障时轮毂电机轴承处的垂向的加速度信号，这样，正常时轮毂电机轴承处的垂向加速信号有5种，故障时轮毂电机轴承处的垂向加速度信号有5种，共10种加速度信号，每种信号采样频率为12800Hz，采样时间为45秒，所采集的所有加速度信号均为时域信号；而轮毂电机运行状态信息S分为正常状态和故障状态。

在10种加速度信号中，将加速度信号类型信息C共计分为10种类型，分别为：加速度信号类型c₁：车速5km/h下轮毂电机的正常运行状态加速度信号类型；加速度信号类型c₂：车速10km/h下轮毂电机的正常运行状态加速度信号类型；加速度信号类型c₃：车速15km/h下轮毂电机的正常运行状态加速度信号类型；加速度信号类型c₄：车速20km/h下轮毂电机的正常运行状态加速度信号类型；加速度信号类型c₅：车速25km/h下轮毂电机的正常运行状态加速度信号类型；加速度信号类型c₆：车速5km/h下轮毂电机的故障运行状态加速度信号类型；加速度信号类型c₇：车速10km/h下轮毂电机的故障运行状态加速度信号类型；加速度信号类型c₈：车速15km/h下轮毂电机的故障运行状态加速度信号类型；加速度信号类型c₉：车速20km/h下轮毂电机的故障运行状态加速度信号类型；加速度信号类型c₁₀：车速25km/h下轮毂电机的故障运行状态加速度信号类型。

步骤二：10种加速信号中，将每种加速度信号等分为15段，再分别计算每段加速度信号时域和频域高敏感特征参数。特征参数的敏感性是指特征参数在区分机械设备不同运行状态时的难易程度，高敏感特征参数越容易识别机械设备的不同运行状态，本发明设计了以提取加速度信号在时域和频域中的多个高敏感特征参数为目标的综合权重诊断指标法(SWDI)，如图2所示是该方法的流程图，具体步骤如下：

(1)针对轮毂电机运行工况，从时域和频域中预选取出8个反映轮毂电机运行状态的特征参数，在将每种加速度信号分段的基础上分别计算10种加速度信号下的8个特征参数值，分别是4个时域信号特征参数和4个频域信号特征参数。

具体的8个特征参数如下所示：

加速度信号的时域信号表示为x_i，i＝1～N，N为采样点数，在本发明中因为加速度信号等分为15段，因此N＝38400，计算出其平均值

和标准差σ分别为：

计算出时域信号特征参数SP₁为：

定义{x_i}中极大值为{x_pi}，其中：j＝1～N_p，N_p为{x_i}中极大值的总数，则可计算出极大值平均值

和极大值标准差

由此再计算出时域信号特征参数SP₂、SP₃、SP₄分别为：

式中，N_v为{x_i}中极小值的总数。

将加速度x_i信号的时域信号经过快速傅里叶变换后转化为频域信号，表示为f_i，其中，i＝1～I，I是采样频率的一半，即I＝6400，F(f_i)是第f_i的频谱值，频域信号特征参数SP₅、SP₆、SP₇、SP₈、根据以下公式进行计算：

计算完成后，得到每种加速度信号下的15组特征参数组，每组特征参数组包含8个特征参数值，例如下表1所示为5km/h车速下某轮毂电机正常运转状态下加速度信号的特征参数组：

表1 5km/h车速下某轮毂电机正常运转状态下加速度信号的特征参数组

组号	SP<sub>1</sub>	SP<sub>2</sub>	SP<sub>3</sub>	SP<sub>4</sub>	SP<sub>5</sub>	SP<sub>6</sub>	SP<sub>7</sub>	SP<sub>8</sub>
									1	0.167	0.394	2.481	0.981	48410.63	0.413	0.013	0.0032
2	0.199	0.426	2.667	0.964	48180.15	0.403	0.019	0.0059
									3	0.169	0.384	2.585	0.919	48308.41	0.412	0.014	0.0034
4	0.144	0.334	2.343	0.967	48903.75	0.433	0.010	0.0030
									5	0.141	0.317	2.246	0.926	48891.71	0.434	0.009	0.0032
6	0.182	0.422	2.534	0.920	47994.67	0.408	0.016	0.0035
									7	0.166	0.372	2.364	0.951	48293.72	0.414	0.013	0.0035
8	0.145	0.335	2.435	0.940	48835.91	0.435	0.010	0.0031
									9	0.155	0.368	2.688	0.950	48616.67	0.422	0.011	0.0031
10	0.225	0.471	2.495	1.02	47279.22	0.399	0.025	0.0062
									11	0.171	0.386	2.522	0.971	48198.96	0.407	0.014	0.0032
12	0.146	0.329	2.256	0.953	48601.9	0.433	0.010	0.0027
									13	0.175	0.397	2.563	0.904	47958.63	0.412	0.015	0.0030
14	0.200	0.430	2.605	0.969	47863.55	0.411	0.019	0.0052
									15	0.142	0.324	2.341	0.951	48500.77	0.431	0.010	0.0025

(2)根据得到的特征参数组，分别计算5种车速下、轮毂电机正常和故障运行状态的8个信号特征参数的平均值μ_SP和标准差σ_SP。

某一车速下某一特征参数的平均值μ_SP和标准差σ_SP用下式计算：

式中(SP)_i表示为该特征参数的第i组的值，因分段为15，所以n＝15。

计算具体结果如下表2、表3所示：

表2

5种车速下、轮毂电机正常和故障运行状态的8个特征参数数值的平均值μ_SP

表3

5种车速下、轮毂电机正常和故障运行状态的8个特征参数数值的标准差σ_SP

(3)根据计算得到的不同加速度信号下各特征参数的平均值μ_SP和标准差σ_SP，再计算各特征参数在各车速工况下正常与故障状态之间的区分度指标值DI(Distinguishindex,DI)如下：

式中，μ_SP1是某一特征参数在某一车速下正常状态对应的平均值，μ_SP2是该特征参数在同一车速下故障状态时对应的平均值，σ_SP1、σ_SP2是对应的标准差。

(4)根据区分度指标值DI计算各特征参数对两种状态的区分率DR(Discrimination Rate,DR)：

x是积分变量，dx是x的变化量。

各特征参数DI值与DR值如表4所示。

表4 DI值与DR值的对应关系

DI	DR
		＜0.85	＜80％
0.85-1.30	80％-90％
		1.30-1.65	90％-95％
1.65-2.33	95％-99％
		＞2.33	＞99％

(5)根据区分度指标值DI、区分率DR计算出综合权重指标值(Synthetic weightdetection index，SWDI)，定义如下：

式中，M代表不同状态(正常状态，故障1状态，故障2状态等)的数目，本发明涉及正常与故障两种状态，因此M取值为2；R代表所要选取高敏感特征参数的数量，需从8个特征参数中优选出4个高敏感特征参数，因此R取值为4；L代表运行工况类别的数目，涉及5种不同车速的工况，因此L取值为5。将大于2.33的DI值全部近似取为2.33，用f(DI)表示；w是权重系数，式中α是高敏感特征参数阈值，取α为95％；y是积分变量，dy是y的变化量。

(6)取综合权重指标值SWDI的最大值所对应的一组特征参数作为优选出的高敏感特征参数组。

例如表5所示，通过计算，第11组SWDI值最大，为45.97，因此选取的高敏感特征参数为SP₁、SP₂、SP₅和SP₇。

表5不同特征参数组的SWDI值

因此，得到了加速度信号的时域高敏感特征参数SP₁、SP₂和频域高敏感特征参数SP₅、SP₇。

步骤三：基于步骤一中整理的轮毂电机在正常和故障状态下的运行状态信息S、加速度信号类型信息C和步骤二计算得出的高敏感特征参数SP₁、SP₂、SP₅、SP₇建立训练数据集。

步骤四：利用专家知识构建当前时间片段k对应的速度片v_k内的贝叶斯网络结构，如图3所示，将具体的将轮毂电机机运行状态S、加速度信号类型C作为隐藏节点，高敏感特征参数SP₁、SP₂、SP₅、SP₇作为可观测的多个高斯节点，则在速度片v_k轮毂电机运行状态节点用S^k表示，加速度信号类型节点用C^k表示，各高敏感特征参数节点分别用SP₁ ^k、SP₂ ^k、SP₅ ^k、SP₇ ^k表示。根据2013年朱金林在杂志《信息与控制》(第42卷第4期)的“基于动态贝叶斯网络的缺失数据系统故障识别”论文中提到参数学习方法，基于步骤三中的训练数据集进行参数学习，得到各网络节点的条件概率分布，具体包括速度片v_k中高敏感特征参数节点SP₁ ^k和加速度信号类型节点C^k的条件概率分布P(SP₁ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₂ ^k和加速度信号类型节点C^k的条件概率分布P(SP₂ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₅ ^k和加速度信号类型节点C^k的条件概率分布P(SP₅ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₇ ^k和加速度信号类型节点C^k的条件概率分布P(SP₇ ^k|C^k)，速度片v_k中轮毂电机运行状态节点S^k的先验概率分布P_f(S^k)，速度片v_k中加速度信号节点C^k的先验概率分布P_f(C^k)。

步骤五：利用专家知识制定如下表6所示的不同“速度片”之间的状态转移概率分布选择表，确定前一时间片段k-1对应的速度片v_k-1与当前时间片段k对应的速度片v_k之间的状态转移概率分布P(S^k-1|S^k)，S^k-1是轮毂电机在前一时间片段k-1对应的速度片v_k-1中的运行状态，S^k是轮毂电机在当前时间片段k对应的速度片v_k中的运行状态。特别的，速度片对应的车速为5km/h，10km/h，15km/h，20km/h和25km/h。

表6轮毂电机二速度片之间的运行状态转移概率分布选择表

步骤六：结合当前时间片段k对应的速度片v_k的贝叶斯网络和二速度片v_k-1-v_k间的状态转移概率分布，建立如图4所示的二速度片v_k-1-v_k轮毂电机机械故障诊断的动态贝叶斯网络模型。将步骤四和步骤五中得到的各网络节点间的条件概率分布和状态转移概率分布设置为所建动态贝叶斯网络的网络参数，具体包括速度片v_k中高敏感特征参数节点SP₁ ^k和加速度信号类型节点C^k的条件概率分布P(SP₁ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₂ ^k和加速度信号类型节点C^k的条件概率分布P(SP₂ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₅ ^k和加速度信号类型节点C^k的条件概率分布P(SP₅ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₇ ^k和加速度信号类型节点C^k的条件概率分布P(SP₇ ^k|C^k)，速度片v_k中轮毂电机运行状态节点S^k的先验概率分布P_f(S^k)，速度片v_k中加速度信号节点C^k的先验概率分布P_f(C^k)以及速度片v_k-1和速度片v_k之间的状态转移概率分布P(S^k-1|S^k)。根据速度片v_k-1和速度片v_k之间的多种状态转移概率分布P(S^k-1|S^k)，建立二速度片轮毂电机机械故障诊断模型群，每个模型的参数中仅状态转移概率分布不同，其他参数一样。

如图5所示，为第二阶段基于机械故障诊断模型群的机械故障的在线诊阶段。将簧下质量加速度采集模块1安装在轮毂电机轴承处，用于在线采集车辆行驶时轮毂电机的振动信号信息，即垂向加速度a信号。采集当前时间片段k内的加速度信号，采样时长为3s，采样频率为12800Hz。簧下质量加速度采集模块1的输出与信号特征参数提取模块3的输入连接，通过信号传输线将垂向加速度a信息传递至信号特征参数提取模块3。

再采用车速采集模块2采集车辆行驶的车速V_k信息。具体为当前时间片段k中间时刻的车速，车速V_k信息从汽车CAN总线通讯系统中读取。车速采集模块2的输出连接诊断模型选择模块4的输入，通过信号传输线将车速V_k信息传递至诊断模型选择模块4。

信号特征参数提取模块3用于对所述垂向加速度信号a的处理，具体为计算加速度信号在时域和频域的高敏感特征参数SP₁、SP₂、SP₅、SP₇，并通过信号传输路线将得到的高敏感特征参数提供给集成诊断推理算法的ECU5。

诊断模型选取模块4中内置了第一阶段离线得到的轮毂电机机械故障诊断群，并储存有前一时间片段k-1的车速信息V_k-1，则根据前一时间片段和当前时间片段的车速V_k-1和V_k从轮毂电机机械故障诊断模型群中选择相应的“二速度片”轮毂电机诊断模型。选择诊断模型的具体步骤如下：

(1)将车速V_k-1和V_k划分为相应的速度片v_k-1和v_k，划分标准如下：

(2)根据得到的v_k-1-v_k从轮毂电机机械故障诊断模型群中选择相应的诊断模型。诊断模型选择模块4的输出与集成诊断推理算法的ECU5的输入连接，将选择的轮毂电机机械故障诊断模型传输至集成诊断推理算法的ECU5，同时将前一时间片段k-1的车速信息V_k-1删除，并储存当前时间片段k的车速信息V_k。

集成诊断推理算法的ECU5集成了诊断推理算法，并储存有前一时间片段k-1对应的速度片v_k-1中轮毂电机运行状态S^k-1的诊断结果。ECU5用于对当前时间片段k对应的速度片v_k中轮毂电机运行状态S^k的在线推理，将高敏感特征参数SP₁、SP₂、SP₅、SP₇作为轮毂电机机械故障诊断模型的输入，从而完成对当前时间片段k对应的速度片v_k中轮毂电机运行状态S^k的诊断，并将轮毂电机运行状态S^k的诊断结果提供给诊断结果显示模块6，同时将前一时间片段k-1对应的速度片v_k-1中轮毂电机运行状态S^k-1诊断结果删除，并储存当前时间片段k对应的速度片v_k中轮毂电机运行状态S^k的诊断结果。

诊断结果显示模块6用于将接收到的轮毂电机运行状态S^k的诊断结果对应到当前的时间片段k，并将结果显示出来提供给驾驶员。

如图6所示，为第二阶段基于诊断模型群的机械故障在线诊断流程图，在线诊断流程是循环进行的，具体步骤如下：

步骤一：在线采集当前时间片段k中的加速度信号a，采集时长为3s，采样频率12800Hz，以及该时间片段中间时刻的车速V_k。

步骤二：计算当前时间片段k中加速度信号a的高敏感特征参数SP₁、SP₂、SP₅、SP₇。

步骤三：根据前一时间片段的车速信息V_k-1和当前时间片段的车速信息V_k从轮毂电机故障诊断模型群中选择相应的二速度片v_k-1-v_k诊断模型。进一步的，所述选择诊断模型具体步骤如下：

(2)根据得到的v_k-1-v_k从轮毂电机机械故障诊断模型群中选择相应的诊断模型。

步骤四：将步骤二中的计算得到的高敏感特征参数作为二速度片诊断模型的输入，对当前时间片段k对应的速度片v_k中轮毂电机运行状态S^k进行推理诊断，即计算轮毂电机运行状态S^k的后验概率分布P_u(S^k)。进一步的，基于二速度片诊断模型的诊断过程如下：

(1)ECU5中储存有前一时间片段k-1对应的速度片v_k-1的轮毂电机运行状态S^k-1的诊断结果即轮毂电机运行状态S^k-1的后验概率P_u(S^k-1)；根据步骤二中计算的高敏感特征参数，输入模型可获得在速度片v_k内各高敏感特征参数节点的概率分布，包括：高敏感特征参数节点SP₁ ^k概率分布P(SP₁ ^k)，高敏感特征参数节点SP₂ ^k概率分布P(SP₂ ^k)，高敏感特征参数节点SP₅ ^k概率分布P(SP₅ ^k)，高敏感特征参数节点SP₇ ^k概率分布P(SP₇ ^k)。

结合速度片v_k中加速度信号类型节点C^k的先验概率分布P_f(C^k)推理得到加速度信号类型节点C^k的后验概率分布：

其中，下标u表示所属概率分布为后验概率分布，下标f表示所属概率分布为先验概率分布。

(2)结合速度片v_k-1和速度片v_k的轮毂电机状态转移概率分布P(S^k-1|S^k)，速度片v_k中加速度信号类型节点C^k与轮毂电机运行状态节点S^k的条件概率分布P(C^k|S^k)，速度片v_k中轮毂电机运行状态节点S^k的先验概率P_f(S^k)推理得到速度片v_k中轮毂电机运行状态节点S^k的后验概率：

(3)根据马尔可夫性，更新得到时间片段k+1对应的速度片v_k+1中加速度信号类型节点C^k+1和轮毂电机运行状态节点S^k+1的先验概率分布：

P_f(C^k+1)＝P_u(C^k)，

P_f(S^k+1)＝P_u(S^k)，

根据后验概率P_u(S^k)判断电机运行状态，当P_u(S^k)≤50％时判断轮毂电机运行状态为正常，当P_u(S^k)>50％时判断轮毂电机运行状态为故障。

步骤五：将速度片v_k中轮毂电机运行状态S^k的诊断结果对应到当前时间片段k中，并将诊断结果提供给驾驶员。

步骤六：采集下一时间片段k+1的加速度信号和该时间片段中间时刻的车速V_k+1，重复步骤二～步骤六。

Claims

1.一种基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是具有以下步骤：

步骤六：根据前一时间片段k-1的车速V_k-1和当前时间片段k的车速V_k从所述的轮毂电机故障诊断模型群中选择相应的二速度片v_k-1-v_k诊断模型；

步骤八：根据所述的后验概率判断轮毂电机正常或故障。

2.根据权利要求1所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是：步骤一中：先采集5km/h，10km/h，15km/h，20km/h和25km/h这5车速种下的正常和故障时轮毂电机轴承处的加速度x_i时域信号，计算出4个时域信号特征参数分别为：

i＝1～N，N为采样点数，x_pi是x_i中的极大值，j＝1～N_p，N_p为x_i中极大值的总数，极大值平均值

极大值标准差

N_v为{x_i}中极小值的总数；

将加速度x_i时域信号变换为频域信号f_i，计算出4个频域信号特征参数分别为：

i＝1～I，I是采样频率的一半；F(f_i)是第f_i的频谱值；

再分别计算5种车速下、轮毂电机正常和故障运行状态的8个信号特征参数的平均值μ_SP、标准差σ_SP、各信号特征参数在各车速工况下正常与故障状态之间的区分度指标值

和区分率

μ_SP1是正常状态对应的平均值，μ_SP2是故障状态对应的平均值，σ_SP1、σ_SP2是对应的标准差；x是积分变量，dx是x的变化量；

最后根据区分度指标值DI、区分率DR计算出综合权重指标值

M＝2，R＝4，L＝5，大于2.33的DI值近似取为2.33，用f(DI)表示；w是权重系数，α是高敏感特征参数阈值，α＝95％，y是积分变量，dy是y的变化量；

取综合权重指标值SWDI的最大值所对应的一组特征参数为高敏感特征参数组，得到高敏感特征参数。

3.根据权利要求2所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是：步骤二中：所述的各网络节点的条件概率分布是速度片v_k中高敏感特征参数节点SP₁ ^k和加速度信号类型节点C^k的条件概率分布P(SP₁ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₂ ^k和加速度信号类型节点C^k的条件概率分布P(SP₂ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₅ ^k和加速度信号类型节点C^k的条件概率分布P(SP₅ ^k|C^k)，速度片v_k中高敏感特征参数节点SP₇ ^k和加速度信号类型节点C^k的条件概率分布P(SP₇ ^k|C^k)，速度片v_k中轮毂电机运行状态节点S^k的先验概率分布P_f(S^k)，速度片v_k中加速度信号节点C^k的先验概率分布P_f(C^k)。

4.根据权利要求3所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是：步骤七中，根据式

计算出轮毂电机运行状态S^k的后验概率分布P_u(S^k)。

5.根据权利要求4所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是：步骤八中，当后验概率分布P_u(S^k)≤50％时判断轮毂电机运行状态为正常，当后验概率分布P_u(S^k)>50％时判断轮毂电机运行状态为故障。

6.根据权利要求1所述的基于动态贝叶斯网络的轮毂电机机械故障在线诊断方法，其特征是：步骤五中，将车速V_k-1和V_k划分为相应的速度片v_k-1和v_k是：