机械与材料工程

基于博弈论组合赋权法的轴承性能退化特征指标选择

  • 王晓燕 , 1, 2 ,
  • 张天晓 2, 3 ,
  • 李景奎 4 ,
  • 付通通 2, 3
展开
  • 1. 沈阳航空航天大学,经济与管理学院,沈阳 110136
  • 2. 沈阳航空航天大学,辽宁省飞机火爆防控及可靠性适航技术 重点实验室,沈阳 110136
  • 3. 沈阳航空航天大学,安全工程学院,沈阳 110136
  • 4. 沈阳航空航天大学,民用航空学院,沈阳 110136

王晓燕(1975—),女,辽宁沈阳人,教授,博士,主要研究方向为机械可靠性工程,E-mail:

收稿日期: 2024-12-23

  修回日期: 2025-02-05

  录用日期: 2025-02-10

  网络出版日期: 2025-08-19

基金资助

辽宁省教育厅基本科研项目(JYTMS20230229)

Selection of performance degradation characteristic indicators of bearings based on game theory combinatorial weighting method

  • Xiaoyan WANG , 1, 2 ,
  • Tianxiao ZHANG 2, 3 ,
  • Jingkui LI 4 ,
  • Tongtong FU 2, 3
Expand
  • 1. College of Economics and Management, Shenyang Aerospace University,Shenyang 110136,China
  • 2. Liaoning Provincial Key Laboratory of Aircraft Fire Control and Reliability Airworthiness Technology, Shenyang Aerospace University,Shenyang 110136,China
  • 3. College of Safety Engineering, Shenyang Aerospace University,Shenyang 110136,China
  • 4. College of Civil Aviation, Shenyang Aerospace University,Shenyang 110136,China

Received date: 2024-12-23

  Revised date: 2025-02-05

  Accepted date: 2025-02-10

  Online published: 2025-08-19

摘要

为了选择更有效的表征空压机轴承性能退化特征指标,提出一种基于博弈论组合赋权法的特征指标选取方法。通过对轴承振动进行时域指标分析及冗余信息预处理,获得轴承时域指标数据集。利用单调性、鲁棒性和趋势性3种特征指标评价方法,分别完成性能特征指标选择。基于此,利用博弈论组合赋权法,对3种特征指标和16种时域指标进行加权,完成轴承性能退化指标的选择,并通过实例验证,说明该方法的有效性。

本文引用格式

王晓燕 , 张天晓 , 李景奎 , 付通通 . 基于博弈论组合赋权法的轴承性能退化特征指标选择[J]. 沈阳航空航天大学学报, 2025 , 42(4) : 37 -44 . DOI: 10.3969/j.issn.2095-1248.2025.04.006

Abstract

In order to select the characteristic indicator to more effectively characterize the performance degradation of the air compressor bearing,a feature indicator selection method based on game theory combinatorial weighting method was proposed.Through the analysis of the time-domain indicator of the bearing vibration and the preprocessing of redundant information,the time-domain indicator dataset of the bearing was obtained.Three characteristic indicator evaluation methods including monotonicity,robustness and trend were used to complete the selection of performance characteristic indicators.Based on this,the game theory combinatorial weighting method was used to weight three characteristic indicator and sixteen time-domain indicator to complete the selection of bearing performance degradation indicator.The effectiveness of this method was illustrated by example verification.

空压机作为主要的动力能源,在石油、汽车和航空航天等许多领域发挥着十分重要的作用。轴承作为空压机中故障率最高的元件,其故障状况影响着整个系统的运行可靠性,因此通常利用监测空压机轴承的振动信号来诊断空压机的运行状态1-3。在空压机运行的过程中,由于受到摩擦、振动、腐蚀等物理化学作用,输出轴轴承性能的退化,功能退化的累积造成零件的故障及停机连锁。因此对空压机输入、输出轴轴承性能的退化情况实时监测是进行空压机及时故障诊断的重要方式。
为了能够选择更有效的评价性能退化的趋势和规律的时域指标,学者们做了大量的研究。马洪涛等4通过时域信号自相关和时域同步平均方法,对振动信号进行降噪处理,得到故障齿轮的峭度指标为最优特征指标。韩毅5针对地铁车辆滚动轴承通过时域分析,得到了方根幅值和绝对平均值可作为诊断轴承性能退化的有效指标。本文利用传感器采集轴承的振动信号数据并进行冗余信息的过滤,获得振动信号的16种时域指标。
为了更加科学地选择特征指标,学者们通过不同的特征指标评价方法进行研究。孙丽等6从时域和频域谱峭度提取信号的特征进行平滑处理,并从单调性的角度出发进行排序,从而选取峰值指标、脉冲指标和谱峭度的均值为最优特征指标。庞世杰等7为了提高轴承预测精度,提取振动信号中的时域特征,对时域特征的趋势性进行研究,得到最优表征特征。然而单一的特征指标评价方法的评价结果不够全面,因此,苏皓南等8分析了时域、频域等特征,捕捉不同的故障信息并根据单调性和趋势性评价方法,得到了均方根和方根幅值等时域指标为优选特征进行研究。周建民等9采用时域方法和集成经验模态分解能量熵提取轴承特征,综合考虑相关性、单调性和鲁棒性,最终得到均方根值、绝对平均幅值等为最优的退化特征指标。
本文根据上述文献轴承性能退化时序性的特点,同时结合已有研究成果的优势,提出对空压机轴承进行时域分析,然后采用不同的特征指标评价方法,并综合这些不同特征的优势,将各个指标进行综合运算,采用博弈论组合赋权法综合主客观赋权方法对特征指标进行加权计算,得到每个特征指标评价方法下每个时域指标的综合权重,由此选出表征空压机轴承性能退化状态的最优时域指标集,来更加精准地对轴承运行状态进行分析。

1 数据采集

本文数据来源于某粮食运输公司空压机轴承传感器的五年的运行数据,包括水平方向振动(VI7051A)、垂直方向振动(VI7051B)和轴承温度(TI7052),轴承相关图例如图1所示。
图1 轴承相关图例

2 数据分析

在实际工况中,通过传感器得到的性能退化数据往往存在着一些冗余信息。为了更好地表征性能退化状态,一般通过时域、频域分析等手段对数据进行预处理。由于性能退化状态信息往往存在严格的时间序列,因此采用时域分析对数据进行预处理。

2.1 空压机轴承振动信号时域分析

时域分析是以时间为自变量,描述振动信号在不同时刻的取值。空压机轴承的时域分析通过振动信号在时间和空间上的变化趋势传递了大量的轴承退化信息。为了更加准确地表达这些信息,采用16种时域指标,有量纲和无量纲时域指标分别如表1表2所示10。表12的公式中,Xi 为传感器收集到的原始数据;N为收集总量。
表1 有量纲时域指标
指标 编号 计算公式 指标 编号 计算公式
均值 1 X ¯ = 1 N i = 1 N x i 最小值 6

X m i n = m i n x i

i = 1,2 , , N

方差 2 σ x 2 = 1 N - 1 i = 1 N x i - X ¯ 2 峰峰值 7 X P = m a x X i - m i n X i
峭度 3 β = 1 N i = 1 N x i 4 歪度 8 α = 1 N i = 1 N x i 3
均方根值 4 X r m s = 1 N i = 1 N x i 2 绝对平均幅值 9 X ¯ = 1 N i = 1 N x i
最大值 5

X m a x = m a x x i ,

i = 1,2 , , N

方根幅值 10 X r = 1 N i = 1 N x i 2
表2 无量纲时域指标
指标 编号 计算公式 指标 编号 计算公式
波动指标 11 S f = X r m s X ¯ 裕度指标 14 C L f = X m a x X r
峰值指标 12 C f = X m a x X r m s 峭度指标 15 K r = β X r m s 4
脉冲指标 13 I f = X m a x X ¯ 倾斜度指标 16 P = α X r m s 3

2.2 特征指标评价方法

为了从16种时域指标中选择更有效表征轴承性能退化的指标,提升故障诊断精度,学者们做了很多探索。但是不同的评价方法选择的特征指标差别较大。

2.2.1 单调性特征指标评价方法

在工程实践中,轴承的退化过程是不可逆的。单调性是一个退化指标的固有属性。一个退化指标的单调性越好,则代表它表征轴承退化轨迹的能力越强11。单调性指标计算方法如公式(1)所示。
M o n X i = 1 N - 1 N o . d i f f x > 0 - N o . d i f f x < 0
式中:Mon(Xi为单调性的有效系数;Mon(Xi取值范围为0到1,Mon(Xi越接近1,则代表其单调性越好;Xi 为时域指标,i=(1,2,…,16);N为退化指标序列的长度;No.diff(x)为一阶差分的数量。

2.2.2 鲁棒性特征指标评价方法

退化指标中包括噪声、随机冲击及异常点等不可控因素,将会降低退化指标的表征退化轨迹的能力。为了量化退化指标的这种能力,将其命名为鲁棒性12,如式(2)所示。
R o b X = 1 N n = 1 N e x p - x n - x n t ¯ x n
式中:Rob(X)为鲁棒性的有效系数,Rob(X)取值范围为0到1,Rob(X)越接近1,则代表鲁棒性越好;xnXtn 时刻的退化指标值; x n t ¯Xtn 时刻的平均趋势值。

2.2.3 趋势性特征指标评价方法

轴承的退化过程是与时间强关联的,随着工作时长的增加,轴承磨损程度会逐渐加重,同时也表现在退化轨迹上,特征指标整体趋势逐渐呈现出递增或递减,故采用趋势性衡量特征与时间的相关性,如公式(3)所示。
T r e X , T = N n = 1 N x n t n - n = 1 N x n · n = 1 N t n N n = 1 N x n 2 - n = 1 N x n 2 N n = 1 N t n 2 - n = 1 N t n 2
式中,xnXtn 时刻的退化指标值;|Tre(X,T)为趋势性的有效系数;Tre(X,T1,其越接近1,趋势性越强13

2.3 博弈论组合赋权法

经上述研究,采用不同的指标评价方法选取的特征指标各不相同,不能确定最优指标。因此为了选取的特征指标能够较为全面地表征轴承的退化性能的三大特性,综合属性最优,拟通过博弈论组合赋权法,对不同的特征指标评价方法和时域特征指标进行加权。

2.3.1 熵值法

熵值法是根据包含的各个变量的变异程度进行赋权的方法,是一种客观赋权方法。在信息论中,熵值法通过计算各特征指标评价方法的熵值来确定它们的权重。主要步骤包括14
1)构建数据源矩阵
Z b = z 11 z 1 m z n 1 z n m
式中:z 1 m 为同一时域指标下不同特征指标评价值;zn 1为同一特征指标评价方法下的时域指标值。
2)数据标准化处理:对构建的数据源矩阵进行正向化、标准化处理。
B i j = z n m i = 1 n z n m
B = b 11 b 1 j b i 1 b i j
3)计算熵值:根据标准化的数据计算各指标熵ej
e j = - 1 l n m · i = 1 n B i j l n B i j
其中:j=1,2,…,m;0≤ ej<1
4)计算指标差异系数:差异系数gj 等于1减去熵值,表示指标的变异程度。
g j = 1 - e j
5)计算权重:对差异系数进行归一化处理,得到各指标权重W j
W j = g j j = 1 m g j

2.3.2 优序图

优序图是通过多个指标重要性程度进行两两比较,最终得到系统指标重要程度的一种主观赋权法。用“0”“0.5”和“1”来表示指标之间的优劣,即“1”表示在两个指标的对比中前者更为重要,“0.5”表示两个指标同等重要,“0”则代表后者更为重要,由此建立各个指标关系表,计算得到各个指标的权重。

2.3.3 基于博弈论方法的组合赋权

该方法将各个评价指标通过构建博弈模型,分析各指标在博弈中的策略选择和行为模式,从而确定各指标的权重,强调在权重确定过程中决策者之间的相互影响和制约关系。具体计算方法如下15
1)利用博弈论思想求取纳什均衡点。
M i n | | w * - w k | | 2 , k = 1,2
式中:||∙|| 2表示矩阵的二范式。
2)根据矩阵的微分性质,将式(10)变换为一阶导数下的线性方程组。
w 1 w 1 T w 1 w 2 T w 2 w 1 T w 2 w 2 T α 1 α 2 = w 1 w 1 T w 2 w 2 T
式中: w 1为熵值法客观赋权法得到的权重值; w 2为利用优序图得到的权重值。
3)求解权重系数α 1α 2,并对其进行归一化处理,得到α 1 *α 2 *
α 1 * = α 1 α 1 + α 2 α 2 * = α 2 α 1 + α 2
4)确定组合权重 w i j
w i j = α 1 * w 1 T + α 2 * w 2 T

2.4 特征指标加权评估

根据2.2中3种特征指标评价方法和式(13)的组合权重,得到不同特征指标评价方法下的不同时域指标的组合数据值。
H M X = M A X w 1 j M O N X + w 2 j T R E X + w 3 j R O B ( X )
HM(X)为不同特征指标评价方法下的不同时域指标中的最优综合评价值。其中, w ij>0 i = 1 3 w i j = 1,数值越大,代表该指标越能精确地表征轴承的性能退化特性。

3 实例验证

以某型空压机轴承运行5年时间,经传感器传输所得的3组运行数据中的水平方向振动数据作为输入数据集对所述方法进行验证。

3.1 数据预处理

为了更好地表征轴承的退化状态,去除原始数据存在的冗余信息,先对空压机轴承水平振动原始数据进行预处理,代入式(1)—(3)中,得到轴承时域指标数据集,如表3所示。
表3 轴承时域指标数据集
时域指标 单调性 位次 鲁棒性 位次 趋势性 位次
均值 0.636 144 578 6 0.735 789 70 14 0.932 007 270 3
方差 0.016 767 460 13 0.990 112 89 12 0.927 302 750 6
峭度 0.012 048 192 14 0.994 918 17 9 0.875 349 346 9
均方根值 0.590 361 440 9 0.997 037 89 3 0.940 076 913 1
MAX 0.010 843 370 15 0.995 480 03 6 0.879 662 330 8
MIN 0.009 638 550 16 0.841 797 58 13 0.451 806 729 14
峰峰值 0.204 819 200 12 0.992 894 52 11 0.883 551 586 7
歪度 0.551 807 220 11 0.994 338 25 10 0.937 960 287 2
绝对平均幅值 0.636 144 570 7 0.997 257 35 2 0.932 007 268 4
方根幅值 0.701 204 810 4 0.715 071 29 15 0.927 434 434 5
波动指标 0.739 759 030 3 0.997 314 31 1 0.731 451 624 10
峰值指标 0.573 493 970 10 0.995 474 90 7 0.678 078 233 11
脉冲指标 0.619 277 100 8 0.995 336 72 8 0.639 755 396 13
裕度指标 0.684 337 340 5 0.586 127 02 16 0.668 321 957 12
峭度指标 0.812 048 190 1 0.995 685 78 5 0.144 893 172 16
倾斜度指标 0.802 409 630 2 0.996 111 89 4 0.318 641 922 15
表3可以看出,同一时域指标在不同评价方法下,计算得到的结果差距较大,而且对性能退化趋势特征表达的强弱程度的位序也存在差异。由于每种方法侧重描述的特征有所不同,因此这种单一特征指标评价方法有可能遗漏掉最佳的表征指标。

3.2 熵值法客观赋权值

采用式(9)对时域指标数据集进行加权计算,求得熵值法权重W j
表4为熵值法变观赋权值。由表4可知,同一时域指标在不同特征指标评价方法下,计算得到的权重值差异较大,并且采用不同的特征指标评价方法下所占权重最大的时域指标也各不相同,因此单一使用熵值法、客观赋权法无法得到最优表征轴承性能退化的时域指标。
表4 熵值法客观赋权值
时域指标 单调性/% 位次 鲁棒性/% 位次 趋势性/% 位次
均值 31.46 10 35.53 3 33.01 11
方差 37.55 3 31.83 13 30.62 16
峭度 37.56 2 31.69 14 30.75 14
均方根值 32.13 7 34.61 5 33.26 8
MAX 37.60 1 31.68 15 30.72 15
MIN 36.46 4 31.31 16 32.23 12
峰峰值 34.98 5 33.02 12 32.00 13
歪度 32.35 6 34.50 7 33.15 9
绝对平均幅值 31.85 8 34.73 4 33.42 7
方根幅值 31.08 12 35.77 1 33.15 9
波动指标 30.97 13 34.61 5 34.42 4
峰值指标 31.75 9 34.06 9 34.19 5
脉冲指标 31.40 11 34.11 8 34.49 3
裕度指标 30.50 14 35.65 2 33.85 6
峭度指标 29.21 16 33.16 11 37.63 1
倾斜度指标 29.77 15 33.72 10 36.51 2

3.3 优序图主观赋权值

经学者对各个指标进行排序对比,根据2.3.2优序图主观赋权法的评价规则,最终得到基于优序图主观赋权值,如表5所示。
表5 优序图主观赋权值
评价方法 单调性 鲁棒性 趋势性 总计 权重/%
单调性 0.5 0 0 0.5 11.11
鲁棒性 1 0.5 0 1.5 33.33
趋势性 1 1 0.5 2.5 55.56

3.4 博弈论组合赋权值

根据式(10)—(13)将主客观赋权法结合,得到空压机轴承性能退化指标的博弈论组合赋权值,如表6所示。
表6 博弈论组合赋权值
时域指标 单调性/% 位次 鲁棒性/% 位次 趋势性/% 位次
均值 30.76 9 35.45 3 33.77 15
方差 33.07 4 32.08 13 34.85 10
峭度 33.12 2 31.97 14 34.91 8
均方根值 31.53 7 34.57 5 33.90 13
MAX 33.12 2 31.96 15 34.93 7
MIN 33.54 1 31.54 16 34.90 9
峰峰值 33.07 4 33.04 12 33.89 14
歪度 31.92 6 34.48 7 33.60 16
绝对平均幅值 31.02 8 34.67 4 34.31 11
方根幅值 30.17 11 35.66 1 34.18 12
波动指标 29.22 13 34.50 6 36.29 3
峰值指标 30.40 10 34.01 9 35.59 5
脉冲指标 29.77 12 34.05 8 36.19 4
裕度指标 28.97 14 35.47 2 35.57 6
峭度指标 25.67 16 33.19 11 41.13 1
倾斜度指标 26.73 15 33.66 10 39.62 2
表6表4对比可知,通过博弈论组合得到的权重值与熵值法单一权重数值与位次均有所不同,因此通过结合优序图赋权法的博弈论组合赋权法得到更加严谨的权重值。

3.5 博弈论组合赋权法评价结果

通过博弈论组合赋权值和特征指标数据集,并根据式(14),最终得到空压机轴承在水平振动数据的博弈论组合赋权法评价结果,如表7所示。
表7 博弈论组合赋权法评价结果
时域指标 单调性 位次 鲁棒性 位次 趋势性 位次 合计 位次
均值 0.195 706 551 7 0.260 872 828 14 0.314 748 703 6 0.771 328 082 6
方差 0.005 544 185 13 0.317 672 133 12 0.323 168 868 1 0.646 385 187 11
峭度 0.003 990 715 14 0.318 026 732 11 0.305 601 486 8 0.627 618 933 13
均方根值 0.186 112 271 8 0.344 707 579 2 0.318 701 985 3 0.849 521 834 2
MAX 0.003 590 942 15 0.318 148 135 10 0.307 288 633 7 0.629 027 710 12
MIN 0.003 233 163 16 0.265 522 765 13 0.157 701 843 14 0.426 457 772 16
峰峰值 0.067 733 955 12 0.328 100 044 9 0.299 392 131 9 0.695 226 130 9
歪度 0.176 146 110 10 0.342 812 091 4 0.315 172 653 5 0.834 130 854 3
绝对平均幅值 0.197 324 892 6 0.345 787 987 1 0.319 745 853 2 0.862 858 732 1
方根幅值 0.211 528 885 3 0.254 982 870 15 0.316 951 873 4 0.783 463 629 5
波动指标 0.216 126 503 1 0.344 042 916 3 0.265 423 771 10 0.825 593 189 4
峰值指标 0.174 342 870 11 0.338 583 484 6 0.241 311 905 11 0.754 238 260 8
脉冲指标 0.184 352 040 9 0.338 881 729 5 0.231 548 837 13 0.754 782 606 7
裕度指标 0.198 225 449 5 0.207 878 523 16 0.237 705 371 12 0.643 809 343 14
峭度指标 0.208 473 982 4 0.330 500 212 8 0.059 600 595 16 0.598 574 789 15
倾斜度指标 0.214 457 206 2 0.335 255 332 7 0.126 236 236 15 0.675 948 774 10
表7可以看出,通过博弈论组合赋权法计算后得到绝对平均幅值时域指标相对于其他时域指标数值最大,说明在同时考虑3种特征指标评价方法时,绝对平均幅值表征轴承性能退化能力更好。通过对比表3表7的数据,使用博弈论组合赋权法可以准确地捕捉到单一特征指标评价方法下不是位于前列位次的时域指标,作为最优表征轴承性能退化的时域指标集,验证了博弈论组合赋权法的准确性和有效性。

4 结论

本文提出一种基于博弈论组合赋权法的特征指标选择方法,并将其应用到空压机轴承振动监测的特征指标选择中,得到如下结论:
1)单一的特征指标评价方法所得到的最优时域指标各不相同,仅能选择出在一种特性下突出优势的指标,并不能综合得到全面表征轴承运行状态的时域指标。
2)将3种特征指标评价方法通过博弈论组合赋权法进行综合研究,相对于单一的特征指标评价方法可以得到更优的表征性能退化状态的指标,为特征指标选择提供了有力的支持。
3)博弈论组合赋权法避免了由于优序图主观赋权法带来的主观偏差,并且结合了熵值法客观赋权法基于数据本身更好的量化能力,保障了时域指标选取的综合性和全面性,使整个特征指标选择更加科学、严谨,评价结果更贴近真实情况。
[1]
Heng A Zhang S Tan A C C,et al.Rotating machinery prognostics:state of the art,challenges and opportunities[J].Mechanical Systems and Signal Processing200923(3):724-739.

[2]
Resta F Ripamonti F Cazzulani G,et al.Independent modal control for nonlinear flexible structures:an experimental test rig[J].Journal of Sound and Vibration2010329(8):961-972.

[3]
Bagordo G Cazzulani G Resta F,et al.A modal disturbance estimator for vibration suppression in nonlinear flexible structures[J].Journal of Sound and Vibration2011330(25):6061-6069.

[4]
马洪涛,李睿,梁雅楠.基于振动时域分析的齿轮故障特征提取[J].汽车实用技术202348(8):162-165.

[5]
韩毅.地铁车辆滚动轴承振动信号的时域分析[J].城市轨道交通研究202124():57-62.

[6]
孙丽,赵俊杰,袁春元,等.基于时域和谱峭度特征融合及指数模型的滚动轴承RUL预测[J].机床与液压202351(10):203-209.

[7]
庞世杰,韩晓明.基于深度特征与随机过程的滚动轴承寿命预测[J/OL].轴承1-9[2024-11-05]

1122.004.html.

[8]
苏皓南,黄倩,胡波,等.基于数据驱动的离心泵轴承特征分析及寿命预测[J].机电工程202441(6):941-955.

[9]
周建民,王发令,张臣臣,等.基于特征优选和GA-SVM的滚动轴承智能评估方法[J].振动与冲击202140(4):227-234.

[10]
王旭.基于维纳过程的分阶段轴承剩余寿命预测研究[D].大庆:东北石油大学,2021.

[11]
Liao L X.Discovering prognostic features using genetic programming in remaining useful life prediction[J].IEEE Transactions on Industrial Electro-nics201461(5):2464-2472.

[12]
Zhang B Zhang L J Xu J W.Degradation feature selection for remaining useful life prediction of rolling element bearings[J].Quality and Reliability Engineering International201632(2):547-554.

[13]
Li N P Lei Y G Liu Z Y,et al.A particle filtering-based approach for remaining useful life predication of rolling element bearings[C]//2014 International Conference on Prognostics and Health Ma-nagement.Cheney:IEEE,2014:1-8.

[14]
林智辉,赖永芳,陈建立,等.基于熵值法-TOPSIS的配电网供电可靠性评价[J].电工技术2024(13):23-26.

[15]
吕安涛,高勇斌,韩汶,等.基于博弈论组合赋权 - TOPSIS法的交通运输企业综合信用评价[J].山东大学学报(理学版)202459(9):88-97.

文章导航

/