CN112330609A - 一种磁悬浮转子系统的故障智能诊断方法 - Google Patents
一种磁悬浮转子系统的故障智能诊断方法 Download PDFInfo
- Publication number
- CN112330609A CN112330609A CN202011130376.6A CN202011130376A CN112330609A CN 112330609 A CN112330609 A CN 112330609A CN 202011130376 A CN202011130376 A CN 202011130376A CN 112330609 A CN112330609 A CN 112330609A
- Authority
- CN
- China
- Prior art keywords
- rotor system
- magnetic suspension
- axis
- suspension rotor
- shape
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000000725 suspension Substances 0.000 title claims abstract description 35
- 238000003745 diagnosis Methods 0.000 title claims abstract description 10
- 238000012706 support-vector machine Methods 0.000 claims abstract description 55
- 238000010586 diagram Methods 0.000 claims abstract description 43
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 34
- 238000012549 training Methods 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims abstract description 18
- 238000012216 screening Methods 0.000 claims abstract description 3
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 230000001788 irregular Effects 0.000 claims 1
- 238000005339 levitation Methods 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000003066 decision tree Methods 0.000 description 8
- 238000013528 artificial neural network Methods 0.000 description 6
- 238000002790 cross-validation Methods 0.000 description 6
- 238000013145 classification model Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 239000002245 particle Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000007635 classification algorithm Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012418 validation experiment Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30164—Workpiece; Machine component
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Quality & Reliability (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种磁悬浮转子系统的故障智能诊断方法,包括以下步骤:采集不同故障情况下的磁悬浮转子系统的轴心轨迹图;对各轴心轨迹图进行提纯去噪处理;通过Hu不变矩理论将各轴心轨迹图转为对应的7个不变矩,7个不变矩形成一个特征向量;通过帝国竞争算法对各轴心轨迹图的特征向量进行筛选精简,剔除不合格轴心轨迹图;用精简后多个不同故障情况下轴心轨迹图的特征向量对有向无环图支持向量机模型进行分类训练;对磁悬浮转子系统的轨迹实时监测,提取故障中磁悬浮转子系统的轴心轨迹图,将提取的轴心轨迹图输入训练后的支持向量机模型中进行分类识别,通过分类识别结果判断出对应的故障类型。对故障中磁悬浮转子系统的轴心轨迹的识别更准确。
Description
技术领域
本发明涉及轴承保护领域,具体涉及一种磁悬浮转子系统的故障智能诊断方法。
背景技术
主动磁悬浮轴承-转子系统(Active Magnetic Bearing,AMB)因其与传统轴承-转子系统相比具有无摩擦、能耗低等多项优点而备受关注。AMB系统具有高度非线性,转子在临界转速下出现的混沌现象是一前沿问题,因此,如何为磁悬浮-转子系统提供一个稳定的控制算法及精准的混沌分析方法成为了该系统的热点研究问题之一。
对于非线性轴承-转子系统的控制,可以通过传统的比例微分控制算法、滑模控制算法等实现,但时滞可能带来的分岔等不稳定因素会对系统控制的稳定性造成不利影响。同时,高速转动下的AMB系统呈现柔性特征,通常用有限元法进行建模分析,但其具有较大的计算复杂度。
AMB系统正常的轴心轨迹应为稳定,长短轴相差不大的椭圆轨迹,由于一些故障情况下便形成其他形状,在混沌过程中表现尤为明显,在跨阶运动后可能会一直保持故障轨迹转动。轨迹形状通常可以直接代表故障类型,通过机器学习中的分类算法对轨迹进行识别可以帮助分析故障原因。
而在实验过程中所采集的轨迹数据会出现大量的冗余,会影响到机器学习的分类效果。因此需要精简数据并通过算法实现智能识别。
发明内容
本发明要解决的技术问题是,针对现有技术存在的上述缺陷,提供了一种磁悬浮转子系统的故障智能诊断方法。
本发明为解决上述技术问题所采用的技术方案是:
一种磁悬浮转子系统的故障智能诊断方法,其特征在于,包括以下步骤:
1)采集不同故障情况下的磁悬浮转子系统的轴心轨迹图;
2)对各轴心轨迹图进行提纯去噪处理;
3)通过Hu不变矩理论将各轴心轨迹图转为对应的7个不变矩,7个不变矩形成一个特征向量;
4)通过帝国竞争算法(ICA)对各轴心轨迹图的特征向量进行筛选精简,剔除不合格轴心轨迹图;
5)用精简后多个不同故障情况下轴心轨迹图的特征向量对有向无环图支持向量机(SVM)模型进行分类训练;
6)对磁悬浮转子系统的轨迹实时监测,提取故障中磁悬浮转子系统的轴心轨迹图,将提取的轴心轨迹图输入训练后的支持向量机(SVM)模型中进行分类识别,通过分类识别结果判断出对应的故障类型。
按上述技术方案,所述步骤1)中提取磁悬浮转子系统轴心轨迹的方法为:
1.1)通过位移传感器捕捉磁悬浮轴承转子系统的轴心位移信号;
1.2)使用数据采集卡将位移信号采样保存,并转存到数据库;
1.3)通过MATLAB对位移信号进行小波变换去噪处理,得到转子的运动轨迹。
按上述技术方案,所述步骤2)中,对各轴心轨迹图进行提纯去噪处理的具体过程包括以下步骤:采用Canny算子算法高斯滤波去除轴心轨迹图噪声,再对轴心轨迹图依次进行边缘提取、二值化、填充处理;最后得到二值化边缘轮廓。
按上述技术方案,在所述的步骤4)中,不合格轴心轨迹图为梅花形或紊乱形干扰轨迹。
按上述技术方案,如果轴心在前后周期中的运动轨迹,其形状和大小不断变化,不能重合,即可判断为轴心轨迹紊乱,将导致故障分类识别的准确性和判断效率减弱。紊乱形轨迹其Hu不变矩特征不收敛,经本智能诊断方法寻优及映射,可从最优解集中予以剔除。因此通过智能诊断方法,精简轴心轨迹不变矩特征数据的冗余数据,即完成不合格的轴心轨迹图像剔除工作。
按上述技术方案,在所述的步骤5)中,先对各种轨迹进行两两组合,分别进行二分类,根据识别的准确率设置最佳的分类节点。
按上述技术方案,所述步骤1)中,轴心轨迹图的形状为4种,分别为外八字形、椭圆形、内八字形和月牙形;
在所述的步骤5)中,有向无环图支持向量机的底层节点分别为外八形、椭圆形、内八形和月牙形。
按上述技术方案,在所述的步骤1)中,所述的不合格轴心轨迹图是指具有不属于用于分类的4种轴心轨迹图的不变矩特征;即其他紊乱型无法得到可用于分类的不变矩特征,均为为不合格轴心轨迹。
按上述技术方案,非故障的轴心轨迹为一个较为稳定的、长短轴近似的椭圆,其他为故障轴心轨迹,故障轨迹分别为月牙状、外八字形、内八字型四种,具有明显的不变矩特征。
按上述技术方案,在所述的步骤5)中,有向无环图SVM是将多个二分类支持向量机组合在一起,对于一个n元的分类问题,一共有n(n-1)/2个二分类器与内部节点,最后产生n个叶子节点,即n个最终分类.训练分类开始时,由第一个节点,即根节点出发一次得到下一层的分类导向,最后得到所需分类。
按上述技术方案,在所述的步骤5)中,将外八字形和内八字形的轨迹分别与椭圆形和月牙形轨迹进行组合,形成外八字形+椭圆形,外八字形+月牙形,内八字形+椭圆形,内八字形+月牙形的两两组合的子集,利用之前寻优后得到的训练样本输入到已设置的有向无环图SVM支持向量机进行训练,然后将测试样本输入到训练的SVM模型中进行测试。
按上述技术方案,将样本集数据被等分为4个子集,每个子集中的样本数为40,每个子集是测试集。
按上述技术方案,有向无环图支持向量机算法结合了图论中有向无环图的思想,将多个二分类支持向量机组合在一起,对于一个n元的分类问题,一共有n个二分类器与内部节点,最后产生n个叶子节点,即n个最终分类。训练分类开始时,由第一个节点,即根节点出发一次得到下一层的分类导向,最后得到所需分类。
由于样本的数量不够多,根据k-fold交叉验证方法,样本集数据被等分为4个子集,每个子集中的样本数为40,每个子集是测试集,任何不重复的3个子集构成4个训练集;k-fold交叉验证方法没有测试实验,而是由验证实验代替。
与C4.5决策树和BP神经网络进行对比,在小样本的情况下,C4.5决策树和BP神经网络的准确性不如SVM的精度.对于C4.5决策树,增加训练样本对信息熵增益率的影响不明显,识别率的增加并不显著.对于BP神经网络,训练样本数量不足以训练模型,并且准确性明显低于SVM和C4.5决策树.与另外两种方法相比,SVM对小样本分类更有效,更适用于磁轴承-转子系统的轴心轨迹识别。
通过表3中各阶不变矩的数据可以看出,同类型的轴心轨迹的Hu不变矩具有一定的相似性,相似的形状其不变矩在相同位置呈现近似的极大值与极小值,同时,不同形状的不变矩也存在一定的差异性,在相应位置不变矩的极大值和极小值数值呈现较大的差异。
根据上述建立的有向无环图SVM模型可知,若每一层的节点设置不佳,会导致特征值相似的两类轨迹不易区分,对最后的识别分类造成的影响较大,这里要先对每两种轨迹组合,分别进行二分类.根据识别的准确率设置最佳的分类节点,二分类的结果如表4所示。
本发明具有以下有益效果:
1、本发明使用帝国竞争算法剔除同种类轴心轨迹图特征向量中的冗余数据和误差数据,使得支持向量机(SVM)学习分类的效率大大提升,支持向量机(SVM)对故障中磁悬浮转子系统的轴心轨迹的识别更准确,为轴心轨迹识别及故障预警提供了新的思路与方法。
2、本发明围绕非线性AMB系统的稳定控制与柔性转子的混沌现象进行深入分析,提出相应的解决算法,构建了仿真计算理论分析平台,对实验获得的混沌状态下AMB系统分别处于正常与故障工况下的轴心轨迹数据进行辨识与分类。
附图说明
图1是本发明实施例中所述的转子四种故障下的轴心轨迹图;
图2是本发明实施例中所述的帝国竞争算法的流程图;
图3是本发明实施例中所述的帝国竞争算法筛选出的3种不合格轴心轨迹图;
图4是本发明实施例中所述的帝国竞争算法与粒子群算法的优化过程对比图;
图5是本发明实施例中所述的有向无环图算法的算法结构图;
图6是本发明实施例中所述的k-fold交叉验证支持向量机的识别结果的线性图。
图7是本发明实施例中所述的相同训练样本的支持向量机、C4.5决策树和BP神经网络分类模型的分类结果对比线性图。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
参照图1~图6所示,本发明提供的一种实施例中磁悬浮转子系统的故障智能诊断方法,包括以下步骤:
1)磁悬浮转子位移传感器捕捉不同故障情况下的磁悬浮轴承转子系统的20组位移信号;
2)使用数据采集卡将位移信号采样保存,并转存到数据库;
3)通过MATLAB对位移信号进行小波变换等去噪处理,得到转子的运动轨迹;
4)采用Canny算子算法高斯滤波,去除轴心轨迹图噪声,边缘提取轴心轨迹图,进行二值化、填充等处理;
5)根据轴心轨迹图的形状将其分为4类,如图1所示,分别为外八字形(图1中a图)、椭圆形(图1中b图)、内八字形(图1中c图)和月牙形(图1中d图);
6)通过Hu不变矩理论将单个轴心轨迹图对应转化为7个不变矩,计算7个不变矩的特征值,7个特征值组成一个特征向量;
7)使用帝国竞争算法(ICA)剔除同种类轴心轨迹图特征向量中的冗余数据和误差数据,从而筛选;
8)使用支持向量机(SVM)中的有向无环图算法对精简后的轴心轨迹图多分类进行训练;
9)提取故障中磁悬浮转子系统的轴心轨迹,将轴心轨迹图输入训练后的支持向量机(SVM)中进行识别,对识别结果进行分析处理即可得到磁悬浮转子系统的故障类型。
进一步地,步骤6)中的Hu不变矩理论,将一维信号经归一化处理后拓展到二维空间,从而得到较为敏感的二维图形的各阶矩。该理论能够提取轴心轨迹真实边缘的特征值,且不变矩具有平移、旋转、比例不变性。不变矩并不是在原图形的各种几何变换下看都是严格一致的,而是在一定范围内波动,每种轴心轨迹则对应在某个特定范围内波动的7个不变矩。
7个不变矩h1—h7的函数式为
其中定义(p+q)阶归一化中心矩ηpq为
中心矩μpq为
图像重心为
原点矩mpq为
通过上述不变矩函数式将单个轴心轨迹图对应转化为7个不变矩,计算7个不变矩的特征值,7个特征值组成一个特征向量,得到20组特征向量,如表1所示,作为训练样本。
表1
通过表1中各阶不变矩的数据可以看出,同类型的轴心轨迹的Hu不变矩具有一定的相似性,相似的形状其不变矩在相同位置呈现近似的极大值与极小值,同时,不同形状的不变矩也存在一定的差异性,在相应位置不变矩的极大值和极小值数值呈现较大的差异。
进一步地,步骤7)中的帝国竞争算法(ICA)其使用原理为:在ICA中,一个帝国的总权力为帝国主义者的权力加上其殖民地权力平均值。一个拥有更大权力的强大帝国将以较少的权力剥夺较弱帝国的殖民地。帝国的力量越弱,就会失去越多的殖民地,直至消失。期间还将有新的帝国,继续重复帝国同化,直至找到全局最优点。
帝国的总权力定义如下
其中impn是第n个帝国的帝国主义国家,TCn是第n个帝国的总权力,NCn是第n帝国的初始殖民地数量,ξ是帝国主义的影响因素,规定了殖民地成本对帝国总权力的影响,并且0<ξ<1。
根据帝国的总权力,选择最弱的帝国的最弱殖民地作为帝国竞争的对象。帝国越强大,就越有可能赢得殖民地。
帝国的标准化总权力是
每个帝国在竞争中赢得最弱殖民地的概率为
较弱的帝国将根据上述概率Pn失去其殖民地,帝国的之间的优势差异变得更加明显,而较弱的帝国最终会因大量殖民地的丧失而消亡。每次迭代后选择权力最大的帝国的Imp为最优解,若满足全局最优解,则停止,否则返回帝国同化,重复帝国竞争。
帝国竞争算法(ICA)的具体步骤如表2所示。
表2
帝国竞争算法流程图如图2所示。
在帝国竞争算法(ICA)对上述实验得到的轨迹数据优化过程中,轨迹在第60次迭代是得到了最优解Y,即样本与标准值间的差值在3.4×10-6附近。经过帝国竞争算法(ICA)寻优后的不变矩特征值如表3所示。
表3
寻优后的图像经过与原图像对比可知,帝国竞争算法将如图3所示的不合格的图像剔除,由此可知,不合格的紊乱形被成功剔除,完整保留了需要识别的4种图形样本,很好的精简了冗余数据。
进一步地,将步骤6)得到20组特征向量,采用粒子群算法(PSO)对水电机组轴心轨迹的冗余数据精简,相同的数据集通过帝国竞争算法(ICA)与粒子群算法(PSO)进行优化得到的对比优化过程如图4所示。
由图4可知,与粒子群算法(PSO)相比,帝国竞争算法(ICA)的收敛速度更快,适合于快速地对得到的轨迹数据进行优化。由此可知帝国竞争算法(ICA)在主动磁悬浮轴承转子系统冗余轨迹数据的精简上具有更优异的成效。
进一步地,步骤8)中的支持向量机最初主要用于二分类情况,但是在实际工程应用中,多类分类问题更为广泛。支持向量机的多分类方法主要分为两类:第一类方法是针对整体中寻找一个多分类模型,同时进行多分类,但这类方法中目标函数十分复杂,计算量巨大,不易实现;第二类方法是通过构造多个二分类模型进行多分类,此类方法应用更为广泛,如二叉树算法,“一对一”、“一对多”多分类算法,有向无环图算法等。
有向无环图算法的算法结构图如图5所示,根据有向无环图算法的算法结构图可知,该算法结合了图论中有向无环图的思想,将多个二分类支持向量机组合在一起,对于一个n元的分类问题,一共有n(n-1)/2个二分类器与内部节点,最后产生n个叶子节点,即n个最终分类。训练分类开始时,由第一个节点,即根节点出发一次得到下一层的分类导向,最后得到所需分类。
根据上述建立的有向无环图SVM模型可知,若每一层的节点设置不佳,会导致特征值相似的两类轨迹不易区分,对最后的识别分类造成的影响较大。这里要先对每两种轨迹组合,分别进行二分类,根据识别的准确率设置最佳的分类节点,二分类的结果如表4所示。
表4
由表4的二分类结果显示,内八与外八,环形与月牙这两个组合的识别率稍低,故分类节点按其它组合设置即可,这里将节点设置为:1、外八形;2、椭圆形;3、内八形;4、月牙形。利用之前寻优后得到的训练样本输入到已设置好节点的有向无环图SVM支持向量机进行训练。
进一步地,使用k-fold交叉验证方法,将测试样本输入到经步骤8)训练后的支持向量机(SVM)中进行测试。根据k-fold交叉验证方法,样本集数据被等分为4个子集,每个子集中的样本数为40,每个子集是测试集,任何不重复的3个子集构成4个训练集。k-fold交叉验证方法没有测试实验,而是由验证实验代替。试验结果如图6所示。
由图6可知,平均正确率达到93.13%,表明训练后的多分类SVM对验证集数据的识别效果比较理想。
进一步地,为了研究样本数量对支持向量机(SVM)分类识别率的影响,采用样本数为96、128、160、192和240共5组不同的训练样本进行训练,对比相同样本数量下的平均识别率,分类结果如表5所示。
表5
从表5的分类结果可以看出,SVM的准确率都会随着样本数量的增加而增加,在128个样本的情况下就已经可以达到很高的识别率。
进一步地,为了验证支持向量机(SVM)的准确性,基于k-fold交叉验证方法对相同训练样本的支持向量机(SVM)、C4.5决策树和BP神经网络分类模型进行对比。分类识别结果如图7所示。
从图7的分类结果可以看出,三种方法的准确度都会随着样本数量的增加而增加。在小样本的情况下,C4.5决策树和BP神经网络的准确性不如支持向量机(SVM)的精度。
对于C4.5决策树,增加训练样本对信息熵增益率的影响不明显,识别率的增加并不显著。
对于BP神经网络,训练样本数量不足以训练模型,并且准确性明显低于支持向量机(SVM)和C4.5决策树。
与另外两种方法相比,支持向量机(SVM)对小样本分类更有效,更适用于磁轴承-转子系统的轴心轨迹识别。
以上的仅为本发明的较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明申请专利范围所作的等效变化,仍属本发明的保护范围。
Claims (10)
1.一种磁悬浮转子系统的故障智能诊断方法,其特征在于,包括以下步骤:
1)采集不同故障情况下的磁悬浮转子系统的轴心轨迹图;
2)对各轴心轨迹图进行提纯去噪处理;
3)通过Hu不变矩理论将各轴心轨迹图转为对应的7个不变矩,7个不变矩形成一个特征向量;
4)通过帝国竞争算法(ICA)对各轴心轨迹图的特征向量进行筛选精简,剔除不合格轴心轨迹图;
5)用精简后多个不同故障情况下轴心轨迹图的特征向量对有向无环图支持向量机(SVM)模型进行分类训练;
6)提取故障中磁悬浮转子系统的轴心轨迹图,将提取的轴心轨迹图输入训练后的支持向量机(SVM)模型中进行分类识别,通过分类识别结果判断出对应的故障类型。
2.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,所述步骤1)中提取磁悬浮转子系统轴心轨迹的方法为:
1.1)通过位移传感器捕捉磁悬浮轴承转子系统的轴心位移信号;
1.2)使用数据采集卡将位移信号采样保存,并转存到数据库;
1.3)通过MATLAB对位移信号进行小波变换去噪处理,得到转子的运动轨迹。
3.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,所述步骤2)中,对各轴心轨迹图进行提纯去噪处理的具体过程包括以下步骤:采用Canny算子算法高斯滤波去除轴心轨迹图噪声,再对轴心轨迹图依次进行边缘提取、二值化、填充处理。
4.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,在所述的步骤4)中,不合格轴心轨迹图为梅花形或紊乱形干扰轨迹。
5.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,在所述的步骤5)中,先对各种轨迹进行两两组合,分别进行二分类,根据识别的准确率设置最佳的分类节点。
6.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,所述步骤1)中,轴心轨迹图的形状为4种,分别为外八字形、椭圆形、内八字形和月牙形;
在所述的步骤5)中,有向无环图支持向量机的底层节点分别为外八形、椭圆形、内八形和月牙形。
7.根据权利要求6所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,在所述的步骤1)中,所述的不合格轴心轨迹图是指具有不属于用于分类的4种轴心轨迹图的不变矩特征。
8.根据权利要求1所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,在所述的步骤5)中,有向无环图SVM是将多个二分类支持向量机组合在一起,对于一个n元的分类问题,一共有n(n-1)/2个二分类器与内部节点,最后产生n个叶子节点,即n个最终分类.训练分类开始时,由第一个节点,即根节点出发一次得到下一层的分类导向,最后得到所需分类。
9.根据权利要求8所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,轴心轨迹图的形状为4种,分别为外八字形、椭圆形、内八字形和月牙形;
在所述的步骤5)中,将外八字形和内八字形的轨迹分别与椭圆形和月牙形轨迹进行组合,形成外八字形+椭圆形,外八字形+月牙形,内八字形+椭圆形,内八字形+月牙形的两两组合的子集,利用之前寻优后得到的训练样本输入到已设置的有向无环图SVM支持向量机进行训练,然后将测试样本输入到训练的SVM模型中进行测试。
10.根据权利要求8所述的磁悬浮转子系统的故障智能诊断方法,其特征在于,将样本集数据被等分为4个子集,每个子集中的样本数为40,每个子集是测试集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011130376.6A CN112330609A (zh) | 2020-10-21 | 2020-10-21 | 一种磁悬浮转子系统的故障智能诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011130376.6A CN112330609A (zh) | 2020-10-21 | 2020-10-21 | 一种磁悬浮转子系统的故障智能诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112330609A true CN112330609A (zh) | 2021-02-05 |
Family
ID=74311777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011130376.6A Pending CN112330609A (zh) | 2020-10-21 | 2020-10-21 | 一种磁悬浮转子系统的故障智能诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112330609A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115143193A (zh) * | 2022-09-01 | 2022-10-04 | 江苏明磁动力科技有限公司 | 基于灵敏度检测结果的磁悬浮电机失稳预诊断方法、系统 |
CN116050030A (zh) * | 2023-04-03 | 2023-05-02 | 亿昇(天津)科技有限公司 | 一种鼓风机转子轴心位置确定方法、装置及设备 |
CN117192371A (zh) * | 2023-11-03 | 2023-12-08 | 南通清浪智能科技有限公司 | 一种新能源汽车电机驱动器的测试方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292243A (zh) * | 2017-06-02 | 2017-10-24 | 太原理工大学 | 一种基于图像处理的转子系统轴心轨迹识别方法 |
CN109030791A (zh) * | 2018-09-07 | 2018-12-18 | 广西电网有限责任公司电力科学研究院 | 一种基于帝国殖民竞争算法的优化svm变压器故障诊断方法 |
CN110569888A (zh) * | 2019-08-21 | 2019-12-13 | 国网江苏省电力有限公司电力科学研究院 | 基于有向无环图支持向量机的变压器故障诊断方法及装置 |
CN110866502A (zh) * | 2019-11-19 | 2020-03-06 | 安徽工业大学 | 基于线性判别分析与粒子群优化支持向量机的故障诊断方法 |
CN111291783A (zh) * | 2020-01-15 | 2020-06-16 | 北京市燃气集团有限责任公司 | 一种燃气调压设备的智能化故障诊断方法、系统、终端及存储介质 |
-
2020
- 2020-10-21 CN CN202011130376.6A patent/CN112330609A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292243A (zh) * | 2017-06-02 | 2017-10-24 | 太原理工大学 | 一种基于图像处理的转子系统轴心轨迹识别方法 |
CN109030791A (zh) * | 2018-09-07 | 2018-12-18 | 广西电网有限责任公司电力科学研究院 | 一种基于帝国殖民竞争算法的优化svm变压器故障诊断方法 |
CN110569888A (zh) * | 2019-08-21 | 2019-12-13 | 国网江苏省电力有限公司电力科学研究院 | 基于有向无环图支持向量机的变压器故障诊断方法及装置 |
CN110866502A (zh) * | 2019-11-19 | 2020-03-06 | 安徽工业大学 | 基于线性判别分析与粒子群优化支持向量机的故障诊断方法 |
CN111291783A (zh) * | 2020-01-15 | 2020-06-16 | 北京市燃气集团有限责任公司 | 一种燃气调压设备的智能化故障诊断方法、系统、终端及存储介质 |
Non-Patent Citations (1)
Title |
---|
张东波: "AMB的轴心轨迹非线性建模与仿真", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》, no. 07, 15 July 2020 (2020-07-15), pages 4 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115143193A (zh) * | 2022-09-01 | 2022-10-04 | 江苏明磁动力科技有限公司 | 基于灵敏度检测结果的磁悬浮电机失稳预诊断方法、系统 |
CN115143193B (zh) * | 2022-09-01 | 2023-01-31 | 江苏明磁动力科技有限公司 | 基于灵敏度检测结果的磁悬浮电机失稳预诊断方法、系统 |
CN116050030A (zh) * | 2023-04-03 | 2023-05-02 | 亿昇(天津)科技有限公司 | 一种鼓风机转子轴心位置确定方法、装置及设备 |
CN116050030B (zh) * | 2023-04-03 | 2023-07-28 | 亿昇(天津)科技有限公司 | 一种鼓风机转子轴心位置确定方法、装置及设备 |
CN117192371A (zh) * | 2023-11-03 | 2023-12-08 | 南通清浪智能科技有限公司 | 一种新能源汽车电机驱动器的测试方法及系统 |
CN117192371B (zh) * | 2023-11-03 | 2024-01-30 | 南通清浪智能科技有限公司 | 一种新能源汽车电机驱动器的测试方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112330609A (zh) | 一种磁悬浮转子系统的故障智能诊断方法 | |
CN107830996B (zh) | 一种飞行器舵面系统故障诊断方法 | |
CN104931263B (zh) | 一种基于符号化概率有限状态机的轴承故障诊断方法 | |
Tang et al. | Latent fingerprint minutia extraction using fully convolutional network | |
Kozerawski et al. | Clear: Cumulative learning for one-shot one-class image recognition | |
CN103366180A (zh) | 一种基于自动特征学习的细胞图像分割方法 | |
CN105930792A (zh) | 一种基于视频局部特征字典的人体动作分类方法 | |
Gumusbas et al. | Offline signature identification and verification using capsule network | |
CN103646255A (zh) | 一种基于Gabor特征和极限学习机的人脸检测方法 | |
CN109344856B (zh) | 一种基于多层判别式特征学习的脱机签名鉴别方法 | |
CN105469080A (zh) | 一种人脸表情识别方法 | |
CN113609569B (zh) | 一种判别式的广义零样本学习故障诊断方法 | |
Oh et al. | Convolutional neural network and 2-D image based fault diagnosis of bearing without retraining | |
Sun et al. | Brushstroke based sparse hybrid convolutional neural networks for author classification of Chinese ink-wash paintings | |
Dahiya et al. | Refining of image using self-organizing map with clustering | |
Park | A feature selection method using hierarchical clustering | |
Lawgali | Recognition of handwritten digits using histogram of oriented gradients | |
Sun et al. | Prototypical inception network with cross branch attention for time series classification | |
CN101877065A (zh) | 小样本条件下的人脸图像非线性鉴别特征抽取和识别方法 | |
CN116756225B (zh) | 一种基于计算机网络安全的态势数据信息处理方法 | |
Eskander et al. | Adaptation of writer-independent systems for offline signature verification | |
CN102332087A (zh) | 一种基于稀疏表示的人脸识别方法 | |
Kumar et al. | k-NN based writer independent offline signature verification system | |
CN108414228A (zh) | 基于平均多粒度决策粗糙集和nnbc轴承故障诊断方法 | |
Yang et al. | Convolutional autoencoder-based sensor fault classification |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210205 |