CN107562979A - 一种基于foa‑wsvdd的滚动轴承性能退化评估方法 - Google Patents
一种基于foa‑wsvdd的滚动轴承性能退化评估方法 Download PDFInfo
- Publication number
- CN107562979A CN107562979A CN201710559685.7A CN201710559685A CN107562979A CN 107562979 A CN107562979 A CN 107562979A CN 201710559685 A CN201710559685 A CN 201710559685A CN 107562979 A CN107562979 A CN 107562979A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- munderover
- msup
- mfrac
- 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.)
- Granted
Links
Classifications
-
- 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
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于FOA‑WSVDD的滚动轴承性能退化评估方法,其特征是:先提取轴承振动数据的时域及时频域的特征向量,并基于单调性进行特征选择;然后利用正常状态特征向量建立性能退化评估模型;针对高斯核函数对早期故障不敏感问题,将小波核函数引入到SVDD算法中;针对SVDD算法参数选择困难问题,以支持向量个数与总样本数的比值为适应度函数,采用改进的FOA对其核参数进行优化,建立FOA‑WSVDD评估模型。最后计算轴承全生命周期振动数据的评估指标,以此来评估轴承的性能退化状况。相比基于高斯核函数建立的SVDD性能退化评估模型,基于FOA‑WSVDD的滚动轴承性能退化评估模型对轴承早期故障更为敏感,与故障程度之间的单调性更好,准确的反映出了轴承的健康状态。
Description
技术领域
本发明属于轴承性能退化评估技术领域,具体涉及一种基于FOA-WSVDD(Fruitfly Optimization Algorithm-Wavelet support vector data description)的滚动轴承性能退化评估方法。
背景技术
随着工业及科技水平的不断提高,机械设备不断在复杂、高效、轻型等方面进行改进,同时也面临更加苛刻的工作环境。一旦设备的关键部件出现故障,就可能会影响整个生产过程,造成巨大的经济损失,甚至会导致人员伤亡等问题。因此,设备维修正由传统的事后维修和计划维修向基于状态的视情维修转变,而作为建立合理维修策略的前提,设备性能退化评估也开始备受关注。
滚动轴承作为旋转机械中的关键零部件之一,其性能状态的好坏直接影响整台设备的运行可靠性。一般来说,滚动轴承在使用过程中都会经历从正常到退化直至失效的过程,而这期间通常要经历一系列不同的性能退化状态。如果能够在滚动轴承性能退化的过程中监测到轴承性能退化的程度,那么就可以有针对性地组织生产和制定合理的维护计划,防止设备异常失效的发生。
目前,对于滚动轴承运行状态的评估可以分为:特征指标法和评估模型法。特征指标法包括:时域指标、频域指标、时频指标。特征指标法能在一定程度上反映出轴承的运行状态,例如时域指标中的稳定性指标(如均方根值、方根幅值等)会随着故障发展而逐渐增大,但无法判断初始损伤的位置。而敏感性指标(如峭度指标)虽可以识别初始损伤的位置,但随着故障的发展会呈现先升后降的趋势,并不符合轴承故障程度的发展趋势。评估模型法就是综合选取各种特征指标,利用智能学习算法,建立性能退化评估模型,从而得到能够全面地反映滚动轴承性能退化过程的特征指标。常用的智能学习算法有逻辑回归(logistic regression,LR),支持向量数据描述(support vector data description,SVDD)等。LR算法在应用过程中,需要使用不同故障模式下的数据建立模型,限制了它的应用范围。SVDD算法仅需要少量正常状态的数据建模,在应用中具有较高的精度和较强的鲁棒性,克服了LR算法的缺点,得到了广泛的应用。但当数据具有多模态分布时,其精度会受到很大的影响,而且不能反映出信号的局部信息。
发明内容
本发明为了更及时准确地反映出轴承的早期退化状态,提出一种基于FOA-WSVDD的滚动轴承性能退化评估方法。
为了实现上述目的,本发明通过如下技术方案实现:
步骤(1):获取轴承正常状态下的振动数据,进行去噪和归一化预处理;提取振动数据的时域、时频域特征向量,并基于单调性进行特征选择;
步骤(2):将FOA算法的寻优空间扩展到三维空间,并改变其迭代步长,得到改进后的FOA算法;
步骤(3):将小波核函数引入到SVDD算法中,利用步骤(1)提取的特征向量建立WSVDD性能退化评估模型,并利用改进的FOA算法对WSVDD模型的核参数寻优;
步骤(4):对轴承全生命周期振动数据同样根据步骤(1)的方法提取特征向量,并作为WSVDD模型的输入,计算特征向量距离超球中心的距离,以此作为轴承的性能退化评估指标;
步骤(5):设置自适应报警阈值曲线,对轴承的早期退化状态做出及时准确的预警。
步骤(6):利用基于经验模态分解和Hilbert包络解调的故障诊断方法验证评估结果进行的正确性。
根据以上的技术方案,可以实现以下的有益效果:
(1)本发明提取了轴承振动数据的时域及基于小波包的时频特征,并基于单调性进行特征选择,使最终的评估指标与故障程度之间具有良好的单调性;
(2)本发明将小波核函数引入到SVDD算法中,构造了基于WSVDD算法的性能退化评估模型,相比于LR算法,本算法仅需要少量正常状态下的数据建模,克服了LR算法对故障模式下的数据的依赖;相比于SVDD算法,本算法对早期的故障更为敏感,更能及时的对早期故障作出预警;
(3)本发明采用改进的FOA,以支持向量个数与总样本数的比值作为适应度函数,对SVDD算法的核参数进行优化,既避免了标准FOA寻优空间不足,易陷入局部最优的缺陷,又消除人为选择核参数的盲目性。
附图说明
图1本发明的实验装置示意图;
图2本发明的具体操作流程。
具体实施方式
为使本发明的目的、技术方案和优点等更加清楚明白,以下结合实施例,并参照附图,对本发明进一步详细说明。
一种基于FOA-WSVDD的滚动轴承性能退化评估方法,操作流程图如图2所示,所述方法包括以下步骤:
步骤(1):提取振动数据的时域、时频域特征向量,并基于单调性进行特征选择;
利用轴承1正常状态下的数据,提取其时域特征RMS(均方根)、AM(绝对值)、SMR(方根幅值)、Kurtosis(峭度)、Skewness(偏斜度)、Peak(峰值),利用db8小波对数据进行三层小波包分解,得到八个节点能量的归一化值,作为时频特征。由于轴承在使用过程中,性能逐渐退化,并且是不可恢复的,所以其性能变化应该是随着时间单调改变的。为了更准确高效的对退化性能进行评估,以单调性作为特征向量质量的评价指标进行特征选择,特征向量的单调性定义为:
其中,x为去噪后的特征向量,K为特征向量的长度,为单位阶跃函数。
计算出每个特征向量的单调性后,从中选取单调性较高的五个特征向量作为最终的特征向量。
步骤(2):将FOA算法的寻优空间扩展到三维空间,并改变其迭代步长,得到改进后的FOA算法;
改进的FOA算法是用来对核参数寻优的,步骤如下:
1)参数初始化:初始化种群规模sizepop;最大迭代次数maxgen,果蝇群体初始位置X_axis、Y_axis、Z_axis,随机飞行范围[a,b];
2)赋予果蝇个体搜索食物的飞行随机方向与距离:
其中,β、η为调节因子,gen当前迭代次数。在寻优的初始阶段,为了增加果蝇种群的多样性,避免陷入局部最优,应增大随机飞行的距离,此时调节因子β大于1,记为β1。在寻优的第二阶段,为了增加寻优精度,随机飞行的距离应减少,此时调节因子β小于1,记为β2。
3)计算每个个体与原点的距离Di以及味道浓度判断值Si:
4)将味道浓度判断值Si代入味道浓度判断函数(适应度函数),求得果蝇个体的味道浓度,找出果蝇种群中味道浓度最佳的果蝇:
Smelli=f(Si)
[bestSmell bestIndex]=min(Smell)
5)保留最佳味道浓度值S与个体对应坐标X、Y、Z,果蝇群体将向该位置飞去:
Smellbest=bestSmell
X_axis=X(bestIndex)
Y_axis=Y(bestIndex)
Z_axis=Z(bestIndex)
6)迭代寻优:重复执行步骤2)~4),并判断当前最佳味道浓度值是否由于上次迭代的最佳味道浓度值,若是执行步骤5),当达到最大迭代次数时,寻优结束。
步骤(3):将小波核函数引入到SVDD算法中,利用步骤(1)提取的特征向量建立WSVDD性能退化评估模型,并利用改进的FOA算法对WSVDD模型的核参数寻优,步骤如下:
1)首先,定义一个训练样本集X={x1,...,xn},其中xi∈Rd(1≤i≤n)为列向量,通过非线性函数把样本特征空间o∈Rd映射到高维超球空间Φ,即φ:x∈O→φ(x)∈Φ,则WSVDD优化问题描述为:
s.t||φ(xi)-a||2≤r2+ξi,ξi≥0,1≤i≤n
其中,a为超球的圆心,r为半径,ξi是引入的松弛项,C>0是控制参数,调节错分训练样本数(球外样本数)和r的大小。
2)优化问题的拉格朗日方程为:
其中α=(α1,...,αn)T≥0,β=(β1,...,βn)T≥0是拉格朗日乘子向量.上式分别对变量r,a,ξi求偏导数并置为0,得:
3)原始问题的对偶形式:
k(xi,xj)为核函数,用它代替内积计算,即这里小波函数取Morlet小波,即此时核函数为:
4)利用改进FOA算法,以最小化支持向量个数Nsv与总样本数N的比值Fsv作为适应度函数,选择最优核参数;原始优化的对偶形式为凸规划问题,利用凸规划问题求解方法求取拉格朗日乘子α,对于新样本z到球心的距离为:
置f(x)=cv(Confidence Value),以CV作为轴承性能退化评估的指标。
步骤(4):对轴承全生命周期振动数据同样根据步骤(1)的方法提取特征向量,并作为WSVDD模型的输入,计算特征向量距离超球中心的距离,以此作为轴承的性能退化评估指标,得到轴承全生命周期的CV值;
步骤(5):设置自适应报警阈值曲线,对轴承的早期退化状态做出及时准确的的预警。
CV是一个连续变化的参数,表示设备偏离正常状态的程度。设定报警阈值,能够对轴承的健康状态进行实时的监控,报警阈值的计算公式如下:
式中:T(t)表示t时刻的CV值,mean、std分别表示求均值和标准偏差。阈值Th的求解分为3个阶段,第1阶段数据来源于早期无故障状态,阈值Th是一个固定值,计算公式如式(1)。第2阶段将t时刻的T(t)与t-1时刻的Th(t-1)进行比较,若T(t)在Th(t-1)范围内,则用式(2)计算Th值。如果其后连续Nu个CV值均超限,则定义t=te为性能退化状态发生变化的时刻,进入第三阶段,计算公式如式(3)。
步骤(6):利用基于经验模态分解和Hilbert包络解调的故障诊断方法验证评估结果进行的正确性。
根据轴承的包络知识可知,当轴承的某个部件发生故障时,其冲击信号的包络谱会在故障特征频率及其倍频处出现幅值逐渐衰减的离散谱线。所以利用基于经验模态分解和Hilbert包络解调的方法,提取轴承振动信号的包络谱,对评估结果的正确性进行验证。
Claims (3)
1.一种基于FOA-WSVDD的滚动轴承性能退化评估方法,具体包括以下几个步骤:
步骤(1):获取轴承正常状态下的振动数据,进行去噪和归一化预处理;提取振动数据的时域、时频域特征向量,并基于单调性进行特征选择;
步骤(2):将FOA算法的寻优空间扩展到三维空间,并改变其迭代步长,得到改进后的FOA算法;
步骤(3):将小波核函数引入到SVDD算法中,利用步骤(1)提取的特征向量建立WSVDD性能退化评估模型,并利用改进的FOA算法对WSVDD模型的核参数寻优;
步骤(4):对轴承全生命周期振动数据根据步骤(1)的方法提取特征向量,并作为WSVDD模型的输入,计算特征向量距离超球中心的距离,以此作为轴承的性能退化评估指标;
步骤(5):设置自适应报警阈值曲线,对轴承的早期退化状态做出及时准确的预警;
步骤(6):利用基于经验模态分解和Hilbert包络解调的故障诊断方法验证评估结果进行的正确性。
2.根据权利要求1所述一种基于FOA-WSVDD的滚动轴承性能退化评估方法,其特征是:所述步骤(2)中对FOA算法寻优空间的扩展和迭代步长的自适应改变避免其陷入局部最优,包括以下步骤:
1)参数初始化:初始化种群规模sizepop;最大迭代次数maxgen,果蝇群体初始位置X_axis、Y_axis、Z_axis,随机飞行范围[a,b];
2)赋予果蝇个体搜索食物的飞行随机方向与距离:
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>X</mi>
<mo>_</mo>
<mi>a</mi>
<mi>x</mi>
<mi>i</mi>
<mi>s</mi>
<mo>+</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mi>o</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>a</mi>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>&beta;</mi>
<mfrac>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
<mi>&eta;</mi>
</mfrac>
</msup>
</mrow>
<mrow>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>Y</mi>
<mo>_</mo>
<mi>a</mi>
<mi>x</mi>
<mi>i</mi>
<mi>s</mi>
<mo>+</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mi>o</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>a</mi>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>&beta;</mi>
<mfrac>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
<mi>&eta;</mi>
</mfrac>
</msup>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>Z</mi>
<mo>_</mo>
<mi>a</mi>
<mi>x</mi>
<mi>i</mi>
<mi>s</mi>
<mo>+</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mi>o</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>a</mi>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>&beta;</mi>
<mfrac>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
<mi>&eta;</mi>
</mfrac>
</msup>
</mrow>
其中,β、η为调节因子,gen当前迭代次数;在寻优的初始阶段,为了增加果蝇种群的多样性,避免陷入局部最优,应增大随机飞行的距离,此时调节因子β大于1,记为β1;在寻优的第二阶段,为了增加寻优精度,随机飞行的距离应减少,此时调节因子β小于1,记为β2;
3)计算每个个体与原点的距离Di以及味道浓度判断值Si:
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>X</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Y</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
4)将味道浓度判断值Si代入味道浓度判断函数(适应度函数),求得果蝇个体的味道浓度,找出果蝇种群中味道浓度最佳的果蝇:
Smelli=f(Si)
[bestSmellbestIndex]=min(Smell)
5)保留最佳味道浓度值S与个体对应坐标X、Y、Z,果蝇群体将向该位置飞去:
Smellbest=bestSmell
X_axis=X(bestIndex)
Y_axis=Y(bestIndex)
Z_axis=Z(bestIndex)
6)迭代寻优:重复执行步骤2)~4),并判断当前最佳味道浓度值是否由于上次迭代的最佳味道浓度值,若是执行步骤5),当达到最大迭代次数时,寻优结束。
3.根据权利要求1所述一种基于FOA-WSVDD的滚动轴承性能退化评估方法,其特征是:所述步骤(3)中基于WSVDD的性能退化评估模型的建立,并利用改进的FOA算法对WSVDD模型的核参数寻优,包括以下步骤:
1)首先,定义一个训练样本集X={x1,...,xn},其中xi∈Rd(1≤i≤n)为列向量,通过非线性函数把样本特征空间o∈Rd映射到高维超球空间Φ,即φ:x∈O→φ(x)∈Φ,则WSVDD优化问题描述为:
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>c</mi>
<mo>,</mo>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
</mrow>
</munder>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>c</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
</mrow>
s.t ||φ(xi)-a||2≤r2+ξi,ξi≥0,1≤i≤n
其中,a为超球的圆心,r为半径,ξi是引入的松弛项,C>0是控制参数,调节错分训练样本数(球外样本数)和r的大小;
2)优化问题的拉格朗日方程为:
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>C</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>|</mo>
<mo>|</mo>
<mi>&phi;</mi>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>-</mo>
<mi>a</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
</mrow>
其中α=(α1,...,αn)T≥0,β=(β1,...,βn)T≥0是拉格朗日乘子向量,上式分别对变量r,a,ξi求偏导数并置为0,得:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>L</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>a</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>&DoubleRightArrow;</mo>
<mi>a</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>L</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>&DoubleRightArrow;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>C</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>L</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>&DoubleRightArrow;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
3)原始问题的对偶形式:
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>&alpha;</mi>
</munder>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<mi>C</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
k(xi,xj)为核函数,用它代替内积计算,即这里小波函数取Morlet小波,即此时核函数为:
<mrow>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>1.75</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mi>a</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
4)利用改进FOA算法,以最小化支持向量个数Nsv与总样本数N的比值Fsv作为适应度函数,选择最优核参数;原始优化的对偶形式为凸规划问题,利用凸规划问题求解方法求取拉格朗日乘子α,对于新样本z到球心的距离为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>1.75</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>a</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&alpha;</mi>
<mi>j</mi>
</msub>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>1.75</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>a</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
置f(x)=cv(Confidence Value),以CV作为轴承性能退化评估的指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710559685.7A CN107562979B (zh) | 2017-07-11 | 2017-07-11 | 一种基于foa-wsvdd的滚动轴承性能退化评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710559685.7A CN107562979B (zh) | 2017-07-11 | 2017-07-11 | 一种基于foa-wsvdd的滚动轴承性能退化评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107562979A true CN107562979A (zh) | 2018-01-09 |
CN107562979B CN107562979B (zh) | 2023-07-11 |
Family
ID=60973146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710559685.7A Active CN107562979B (zh) | 2017-07-11 | 2017-07-11 | 一种基于foa-wsvdd的滚动轴承性能退化评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107562979B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108733899A (zh) * | 2018-05-02 | 2018-11-02 | 上海理工大学 | 频域响应计算的精密机床动态性能优化方法 |
CN109241915A (zh) * | 2018-09-11 | 2019-01-18 | 浙江大学 | 基于振动信号平稳非平稳性判别及特征甄别的智能电厂泵机故障诊断方法 |
CZ307847B6 (cs) * | 2018-02-14 | 2019-06-19 | ZKL - Výzkum a vývoj a.s. | Způsob konstruování valivých ložisek s využitím integrovaných počítačových nástrojů |
CN109974782A (zh) * | 2019-04-10 | 2019-07-05 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及*** |
CN110046423A (zh) * | 2019-04-12 | 2019-07-23 | 哈工大机器人(合肥)国际创新研究院 | 一种上楼辅助装置的故障预警方法及*** |
CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
CN111723449A (zh) * | 2020-06-30 | 2020-09-29 | 华北科技学院 | 一种矿井提升机恒减速制动***性能退化评价方法 |
CN112765703A (zh) * | 2020-12-23 | 2021-05-07 | 桂林电子科技大学 | 一种危险品码头结构性能退化评估方法 |
CN113139251A (zh) * | 2021-04-23 | 2021-07-20 | 桂林电子科技大学 | 一种优化主题相关性分析的变工况滚动轴承故障诊断方法 |
CN113627088A (zh) * | 2021-08-23 | 2021-11-09 | 上海交通大学 | 基于基因编程和数据融合的机器性能退化评估方法及*** |
CN115901249A (zh) * | 2022-11-07 | 2023-04-04 | 昆明理工大学 | 结合特征优选与多策略优化svdd的滚动轴承性能退化评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854015A (zh) * | 2012-10-15 | 2013-01-02 | 哈尔滨理工大学 | 一种滚动轴承故障位置及性能退化程度诊断方法 |
CN105528504A (zh) * | 2016-03-01 | 2016-04-27 | 哈尔滨理工大学 | 基于cfoa-mkhsvm的滚动轴承健康状态评估方法 |
US20160282229A1 (en) * | 2014-08-28 | 2016-09-29 | Beijing Jiaotong University | Fault Prediction and Condition-based Repair Method of Urban Rail Train Bogie |
CN106644481A (zh) * | 2016-12-27 | 2017-05-10 | 哈尔滨理工大学 | 基于数学形态学和ifoa‑svr的滚动轴承可靠度预测方法 |
-
2017
- 2017-07-11 CN CN201710559685.7A patent/CN107562979B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854015A (zh) * | 2012-10-15 | 2013-01-02 | 哈尔滨理工大学 | 一种滚动轴承故障位置及性能退化程度诊断方法 |
US20160282229A1 (en) * | 2014-08-28 | 2016-09-29 | Beijing Jiaotong University | Fault Prediction and Condition-based Repair Method of Urban Rail Train Bogie |
CN105528504A (zh) * | 2016-03-01 | 2016-04-27 | 哈尔滨理工大学 | 基于cfoa-mkhsvm的滚动轴承健康状态评估方法 |
CN106644481A (zh) * | 2016-12-27 | 2017-05-10 | 哈尔滨理工大学 | 基于数学形态学和ifoa‑svr的滚动轴承可靠度预测方法 |
Non-Patent Citations (2)
Title |
---|
BIN ZHANG ET AL: "Degradation Feature Selection for Remaining Useful Life Prediction of Rolling Element Bearings" * |
何青 等: "基于EEMD 和MFFOA-SVM 滚动轴承故障诊断" * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CZ307847B6 (cs) * | 2018-02-14 | 2019-06-19 | ZKL - Výzkum a vývoj a.s. | Způsob konstruování valivých ložisek s využitím integrovaných počítačových nástrojů |
CN108733899A (zh) * | 2018-05-02 | 2018-11-02 | 上海理工大学 | 频域响应计算的精密机床动态性能优化方法 |
CN109241915A (zh) * | 2018-09-11 | 2019-01-18 | 浙江大学 | 基于振动信号平稳非平稳性判别及特征甄别的智能电厂泵机故障诊断方法 |
CN109241915B (zh) * | 2018-09-11 | 2022-03-25 | 浙江大学 | 基于振动信号平稳非平稳性判别及特征甄别的智能电厂泵机故障诊断方法 |
CN109974782A (zh) * | 2019-04-10 | 2019-07-05 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及*** |
CN109974782B (zh) * | 2019-04-10 | 2021-03-02 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及*** |
CN110046423B (zh) * | 2019-04-12 | 2023-01-17 | 合肥哈工慈健智能科技有限公司 | 一种上楼辅助装置的故障预警方法及*** |
CN110046423A (zh) * | 2019-04-12 | 2019-07-23 | 哈工大机器人(合肥)国际创新研究院 | 一种上楼辅助装置的故障预警方法及*** |
CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
CN111597651B (zh) * | 2020-04-30 | 2023-05-02 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
CN111723449A (zh) * | 2020-06-30 | 2020-09-29 | 华北科技学院 | 一种矿井提升机恒减速制动***性能退化评价方法 |
CN111723449B (zh) * | 2020-06-30 | 2024-02-20 | 华北科技学院 | 一种矿井提升机恒减速制动***性能退化评价方法 |
CN112765703A (zh) * | 2020-12-23 | 2021-05-07 | 桂林电子科技大学 | 一种危险品码头结构性能退化评估方法 |
CN112765703B (zh) * | 2020-12-23 | 2022-06-21 | 桂林电子科技大学 | 一种危险品码头结构性能退化评估方法 |
CN113139251B (zh) * | 2021-04-23 | 2023-03-14 | 桂林电子科技大学 | 一种优化主题相关性分析的变工况滚动轴承故障诊断方法 |
CN113139251A (zh) * | 2021-04-23 | 2021-07-20 | 桂林电子科技大学 | 一种优化主题相关性分析的变工况滚动轴承故障诊断方法 |
CN113627088A (zh) * | 2021-08-23 | 2021-11-09 | 上海交通大学 | 基于基因编程和数据融合的机器性能退化评估方法及*** |
CN113627088B (zh) * | 2021-08-23 | 2024-04-09 | 上海交通大学 | 基于基因编程和数据融合的机器性能退化评估方法及*** |
CN115901249A (zh) * | 2022-11-07 | 2023-04-04 | 昆明理工大学 | 结合特征优选与多策略优化svdd的滚动轴承性能退化评估方法 |
CN115901249B (zh) * | 2022-11-07 | 2024-02-27 | 昆明理工大学 | 结合特征优选与多策略优化svdd的滚动轴承性能退化评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107562979B (zh) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107562979A (zh) | 一种基于foa‑wsvdd的滚动轴承性能退化评估方法 | |
CN106682814B (zh) | 一种基于故障知识库的风电机组故障智能诊断方法 | |
CN107341349B (zh) | 风机健康评估的方法、***、存储器及控制器 | |
TWI669617B (zh) | 設備健康狀態監控方法及其系統 | |
Kusiak et al. | Analyzing bearing faults in wind turbines: A data-mining approach | |
An et al. | Self-learning transferable neural network for intelligent fault diagnosis of rotating machinery with unlabeled and imbalanced data | |
CN111597651B (zh) | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 | |
CN104632521B (zh) | 一种基于偏航校正的风电功率优化***及其方法 | |
CN106885697A (zh) | 基于fcm‑hmm的滚动轴承的性能退化评估方法 | |
CN113570138B (zh) | 一种时间卷积网络的设备剩余使用寿命预测方法及装置 | |
Masseran | Integrated approach for the determination of an accurate wind-speed distribution model | |
CN108898251A (zh) | 考虑气象相似性和功率波动的海上风电场功率预测方法 | |
CN109858104A (zh) | 一种滚动轴承健康评估与故障诊断方法及监测*** | |
CN112016251B (zh) | 一种核动力装置故障的诊断方法及*** | |
CN105129109A (zh) | 一种基于多重分形理论和自组织映射网络的飞机副翼作动器***健康评估方法 | |
CN102135021B (zh) | 一种工业抽汽冷凝式蒸汽轮机的轴功率预测方法 | |
CN103413174A (zh) | 基于深度学习方法的短期风速多步预测方法 | |
CN112033463A (zh) | 一种核动力设备状态评估与预测一体化方法和*** | |
Baraldi et al. | An ensemble approach to sensor fault detection and signal reconstruction for nuclear system control | |
CN109062180A (zh) | 一种基于ifoa优化svm模型的油浸式电抗器故障诊断方法 | |
CN104218571B (zh) | 一种风力发电设备的运行状态评估方法 | |
CN111191832A (zh) | 一种台风灾害配电网杆塔故障预测方法及*** | |
Yang et al. | Twin broad learning system for fault diagnosis of rotating machinery | |
CN104268662B (zh) | 一种基于分步优化分位数回归的沉降预测方法 | |
CN116976227B (zh) | 一种基于lstm机器学习的风暴增水预报方法及*** |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |