CN101738314A - 一种基于抗体浓度的齿轮故障诊断方法 - Google Patents

一种基于抗体浓度的齿轮故障诊断方法 Download PDF

Info

Publication number
CN101738314A
CN101738314A CN200910262984A CN200910262984A CN101738314A CN 101738314 A CN101738314 A CN 101738314A CN 200910262984 A CN200910262984 A CN 200910262984A CN 200910262984 A CN200910262984 A CN 200910262984A CN 101738314 A CN101738314 A CN 101738314A
Authority
CN
China
Prior art keywords
antibody
antigen
gear
fault
clone
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
Application number
CN200910262984A
Other languages
English (en)
Other versions
CN101738314B (zh
Inventor
刘韬
尚丽
周燕
赵志强
陈杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Modern Enterprise Informatization Application And Support Software Engineering Technology Research And Development Center
Original Assignee
Jiangsu Modern Enterprise Informatization Application And Support Software Engineering Technology Research And Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu Modern Enterprise Informatization Application And Support Software Engineering Technology Research And Development Center filed Critical Jiangsu Modern Enterprise Informatization Application And Support Software Engineering Technology Research And Development Center
Priority to CN200910262984XA priority Critical patent/CN101738314B/zh
Publication of CN101738314A publication Critical patent/CN101738314A/zh
Application granted granted Critical
Publication of CN101738314B publication Critical patent/CN101738314B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Peptides Or Proteins (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了基于抗体浓度的齿轮故障诊断方法,基于齿轮故障问题与生物免疫系统所遇到的问题的相似性,齿轮工作过程中出现了故障,类似于免疫系统受到了抗原的入侵,利用自适应免疫网络模型,将齿轮故障作为抗原,建立抗体种类及其浓度与齿轮故障种类(对应于抗原)及其强度关系模型,通过对模型中抗体种类和抗体浓度的统计,判断齿轮故障是否存在和故障出现时故障的种类。

Description

一种基于抗体浓度的齿轮故障诊断方法
技术领域
本发明涉及一种齿轮故障诊断方法,特别涉及一种利用人工免疫系统模型,建立抗体种类及浓度与齿轮故障种类及强度关系的齿轮故障诊断方法。
背景技术
齿轮传动是机械设备中最常用的传动方式之一。齿轮由于经常工作于高速、高负荷、高冲击的环境下,容易产生磨损、断齿、裂纹等多种故障,并进一步诱发其他机械故障。因此,对齿轮的故障诊断具有重要的意义。
当齿轮发生故障时,其振动信号大多为多分量的调幅-调频信号,因此解调分析是业界常用的齿轮故障诊断方法,从振动信号中提取调制信号、分析其强度和频次就可以判断零件损伤的程度和部位。现有的齿轮故障诊断方法,如基于支持向量机的诊断方法和基于神经网络的诊断方法等,仅只能对判断出齿轮有无故障和故障类型,但不能判断出故障发生的严重程度,本发明提供的方法不仅可以判断出齿轮有无故障和故障类型,还可以计算出故障发生的严重程度,这对于评价齿轮的工作状态是非常重要的。
发明内容
本发明的发明目的是为了解决现有齿轮故障诊断信号存在的上述问题,提供一种基于抗体浓度的齿轮故障诊断方法。
本发明的技术方案包括以下步骤:
1)获取齿轮振动信号,通过振动传感设备,测试提取诊断齿轮正常工作和发生故障时的有关振动信号;
2)提取特征向量,采用小波包方法,提取正常工作和发生故障时的特征向量值,对每一信号经小波包分解,求各频带的能量,得到一组与信号相对应的能量序列{Ej|j=1,2,...,2N},由此确定与此能量序列对应的特征向量 P = ( E 1 , E 2 , . . . , E 2 N ) . 选择齿轮磨损、断齿、裂纹和正常状态进行故障特征提取,将选择的振动信号时域波形,用小波包方法进行分解,分别提取不同频段的特征信号,即由各频段的能量谱求得齿轮发生故障时其振动信号的特征向量值;
3)建立诊断系统,选择齿轮正常工作和不同故障类型发生时的特征向量作为训练样本,生成诊断系统,系统采用自适应免疫网络模型产生抗体群体,采用故障分类及其强度评价模型计算各种故障类型的故障强度系数,计算所得的齿轮不同故障状态下的故障强度系数作为具体诊断时的故障判断标准;
4)输入诊断数据,将第一步检测并通过第二步特征提取所获得的待诊断齿轮故障特征向量值作为抗原,输入到建立的故障诊断系统,系统发生免疫应答,产生相应的故障特征抗体群体,然后,计算不同故障类型的强度系数;
5)输出诊断结果,根据强度系数判断出可能存在的故障类型并输出相应结果。
上述步骤3)中的所述的自适应免疫网络模型的算法描述如下:
①将待处理的数据定义为抗原集合AG={ag1,ag2,...,agI},抗体识别抗原的过程就是免疫应答识别待处理数据模式的过程。通过免疫应答识别过程,获得的记忆细胞集合M即为需要求解的模式。设系统中的抗原具体表示为:
AG = ( ag i . . . ag l ) T = ag 11 K ag 1 m M O M ag I 1 L ag Im - - - ( 1 )
agij∈[0,1]m,i=1,...,I;j=1,...,m。式(1)中:agi表示第i个抗原;agij表示第i个抗原的第j个向量。通过随机产生抗体:
AB = ( ab 1 . . . ab k ) T = ab 11 K ab 1 m M O M ab k 1 L ab km - - - ( 2 )
abij∈[0,1]m,i=1,...,k;j=1,...,m。作为新生成的抗体。式中:abi表示第i个抗体,abij表示第i个抗体的第j个向量。
②把抗原与抗体间的欧氏距离作为衡量二者之间匹配程度的一个指标。将抗原中的每个元素分别与抗体中的每个元素相匹配,对任意抗原agi和任意抗体abj,二者之间的距离d(i,j)可通过下式计算:
d ( i , j ) = ( Σ t = 1 m ( ag it - ab jt ) 2 ) 1 / 2 - - - ( 3 )
对任意抗体abi和任意抗体abj,二者之间的距离s(i,j)可通过下式计算:
s ( i , j ) = ( Σ t = 1 m ( ab it - ab jt ) 2 ) 1 / 2 - - - ( 4 )
对任意abj(abj∈AB,j=1,L,k),所受到的来自抗原agj(agi∈AG,i=1,L,I)的激励sl(i,j)为:
sl ( i , j ) = 1 / d ( i , j ) d ( i , j ) < T sl 0 other - - - ( 5 )
i=1,2,...,I;j=1,2,..,k。在式(5)中:
T sl = k sl &Sigma; i = 1 I &Sigma; j = 1 k d ( i , j ) I &times; k 2 - - - ( 6 )
Tsl表示抗原对抗体激励的阈值,其中:ksl∈(0,1]且为常数。
③当抗体abj所受到来自抗原的刺激达到一定程度(d(i,j)<Tsl时,免疫系统被激活,系统开始对抗体进行克隆增殖。增殖是指免疫细胞发生无性、有丝分裂。自我复制(没有交叉)的过程。抗体克隆的数目与其所受到的激励水平成正比,受激励越大的抗体克隆的数目越多;受激励越小的抗体克隆的数目越少。由于受抗原agi激励所导致的抗体abj的克隆数目为:
式(7)中:
Figure G200910262984XD00044
表示取与sl(i,j)最接近的整数。抗体abj经agi刺激产生克隆的集合为Abji,抗体abj受抗原群体刺激所产生的克隆集合Cj为:
C j = Ab j 1 &cup; Ab j 2 &cup; L &cup; Ab jc 1 ( i , j ) = U t = 1 c 1 ( i , j ) Ab jt - - - ( 8 )
相应的抗体集合AB经克隆过程后新生成的抗体集合C为:
C = C 1 &cup; C 2 &cup; L &cup; C t = U t = 1 k C t - - - ( 9 )
克隆细胞要经历突变(伴随着高变异率)过程,从而使得新产生的细胞与所选择的抗原以更高的亲和力相匹配。在变异过程中每种免疫细胞变异的速率与细胞受体受抗原的激励程度的大小成反比:其受抗原的激励越大,变异越小;其受抗原的激励越小,变异越大。定义:
Mut(i,j)=d(i,j)    (10)
为抗体abj受到抗原agi刺激后新产生的抗体的变异率。由抗原agi入侵引起的abj的克隆所生成的cj变异为:
c′j=cj+mut(i,j)×(cj-agi)
(11)
抗体abj被抗愿群刺激产生的变异集合为C′j,变异产生的全部抗体集合为:
C &prime; = C 1 &prime; &cup; C 2 &prime; &cup; L &cup; C t &prime; = U t = 1 k C t &prime; - - - ( 12 )
变异保证了新产生的抗体具有充足的多样性,这是免疫系统强化学习机制的内在的特性,系统通过重复进行上述工作来增强自己识别抗原的性能。
④克隆变异之后首先进行克隆抑制(Clone suppression)。对集合C′,再次计算抗原与克隆变异后的抗体细胞的亲和力,计算克隆变异后的抗体细胞之间的亲和力。同时计算抗原-抗体抑制阈值Td和抗体-抗体抑制阈值Ts
T d = k d 1 p &Sigma; j = 1 p min i = 1 I d ( i , j ) - - - ( 13 )
式(13)中:kd∈(0,1]且为常数
T s = k s &Sigma; i = 1 p &Sigma; j = 1 p s ( i , j ) p ( p - 1 ) 2 - - - ( 14 )
式(14)中:ks∈(0,1]且为常数
在克隆抑制阶段,首先删除那些与抗原亲和力大(于阈值Td)的抗体细胞,然后再根据抗体细胞之间的亲和力,清除亲和力小(于阈值Ts)的记忆细胞,此阶段完成后形成的抗体集合为C″,其元素满足以下条件:
ab j = ab j &Exists; ag i d ( i , j ) < T d 0 other - - - ( 15 )
ab j = ab j &Exists; ab i s ( i , j ) < T s 0 other - - - ( 16 )
通过克隆抑制阶段,形成的抗体集合C″比C′的规模缩减。随后再进行克隆选择,对集合C″的抗体,按照其与抗原的亲和力大小,选择ξ%克隆变异后的细胞作为该抗体的记忆细胞,ξ%可以预先设定。将选择出的变异克隆细胞放入记忆细胞集合M中。
⑤经克隆变异后的抗体中那些具有较高抗原刺激度的抗体转化为具有较长寿命的记忆细胞,而那些具有较低抗原刺激度的抗体则随着抗原的消亡而消亡。由于抗体细胞之间的互相识别易引起自身免疫应答反应,因此,为了消除自身免疫,在每次加入新的抗体记忆细胞后,需要对那些能够自我识别的抗体进行抑制。
设定记忆细胞集合中元素的个数为p1,定义记忆细胞亲和力成熟阈值为:
T s &prime; = k s &prime; &Sigma; i = 1 p 1 &Sigma; j = 1 p 1 s ( i , j ) p 1 &times; ( p 1 - 1 ) 2 - - - ( 17 )
式(17)中:k′s∈(0,1]且为常数
那么抗体经亲和力成熟过程后转化后,记忆细胞集合M中的元素应该满足以下条件:
ab j = ab j s ( i , j ) > T s &prime; 0 other - - - ( 18 )
式(18)中:i=1,2,...,p1;j=1,2,...,p1
⑥网络训练终止条件判断,任给一个抗体agi,在记忆细胞集合里都可以找到一个抗体abj,判断是否满足:d(i,j)≤T′s。如果满足,则训练结束,否则转步骤③继续。
上述步骤3)中所述的故障分类及其强度评价模型为:基于齿轮故障问题与生物免疫系统所遇到的问题的相似性,齿轮工作过程中出现了故障,类似于免疫系统受到了抗原的入侵,利用自适应免疫网络模型,将齿轮故障作为抗原,建立抗体种类及其浓度与齿轮故障种类(对应于抗原)及其强度关系模型,通过对模型中抗体种类和抗体浓度的统计,判断齿轮故障是否存在和故障出现时故障的种类;
设定ni为对应于齿轮的第i类故障发生时所产生的抗体数量,ci为齿轮正常状态下系统所产生的抗体数量,当齿轮的第i类故障发生时,定义其故障强度系数为:
r i = 1 - 1 1 + lg ( | n i - c i | + 1 ) - - - ( 19 )
由本公式的定义可知:0≤ri≤1,当ri=0时,表明齿轮当前工作状态下没有第i类故障,当ri=1时,表明第i类故障存在,当0<ri<1时,存在第i类故障的可能性,且ri值越大,表明当前齿轮面临的第i类故障可能性越高,在模型中,当没有故障存在时,系统各种抗体数目相对稳定;当故障存在时各抗体的数目急剧增加,因此,可以根据故障存在时刺激所产生的抗体种类及其浓度来识别故障种类及强度。
本发明采用自适应免疫网络模型产生抗体群体,采用故障分类及其强度评价模型计算各种故障类型的故障强度系数,计算所得的齿轮不同故障状态下的故障强度系数作为具体诊断时的故障判断标准;将检测到得并通过特征提取所获得的待诊断齿轮故障特征向量值作为抗原,输入到建立的故障诊断系统,系统发生免疫应答,产生相应的故障特征抗体群体,然后,计算不同故障类型的强度系数;根据强度系数判断出可能存在的故障类型并输出相应结果。本发明提供的方法不仅可以判断出齿轮有无故障和故障类型,还可以计算出故障发生的严重程度。
附图说明
图1为本发明的流程图;
图2为齿轮振动信号波形图;
具体实施方式
附图非限制的公开了本发明的一个具体实施例,下面结合附图为本发明做进一步描述。
如图1所示本发明的诊断流程如下:
(1)获得齿轮振动信号,齿轮故障包括齿面磨损、裂纹、轮齿折断、塑性变形、传动轴的弯曲变形以及轴承故障等多个方面,检测故障最有效的方法是检测箱体或轴的振动信号,其振动信号中包含了复杂的多种频率成份及噪声,其中包含齿轮故障产生的附加振动信号,使用测振传感器从3个处于工作状态下的C618型车床主轴传动齿轮箱轴端部测取振动波形,经过低通滤波对检测到的振动信号进行采样,采样频率为2500Hz,以每20个采样值为一组,检测到的2幅振动信号图附图2所示;
(2)特征向量提取,根据诊断实现步骤中的第二步对获取的齿轮振动信号进行特征提取,利用建立的小波包工具分解检测到的振动信号,得到一组与信号相对应的能量序列{Ej|j=1,2,...,2N},并将其作为特征向量 P = ( E 1 , E 2 , . . . , E 2 N ) , 这里将测试获得的10组振动信号时域波形,用小波包方法进行分解,分解层数为5层,分别提取第5层从低频到高频的8个频率段的特征信号,然后求各频段的能量值作为变速箱齿轮振动信号的特征向量值,并对这些数据进行归一化处理,使它们都落入[0,1]上,其中一台设备的20个特征向量值举例如下(见附表1);
表1部分齿轮特征向量值
Figure G200910262984XD00101
(3)建立诊断系统,建立诊断系统包括构建基于自适应免疫网络模型和构建故障分类及其强度评价模型,并对构建的系统进行训练,系统训练时,选取了C618型车床主轴传动齿轮箱在磨损、断齿、裂纹和正常状态下的典型振动信号数据200组作为训练样本,以建立特征向量与故障种类及强度之间的关系,通过训练所获得的对应关系见附表2。
表2特征向量与故障种类及强度之间的关系
Figure G200910262984XD00102
(4)输入诊断数据,将提取的待诊断齿轮故障特征向量值作为抗原,输入到建立的故障诊断系统,系统发生免疫应答,产生相应的故障特征抗体群体,然后,计算不同故障类型的强度系数。
(5)输出诊断结果,根据强度系数判断出可能存在的故障类型并输出相应结果,本例诊断的三个C618型车床主轴传动齿轮箱故障类型及强度见附表3。
表3诊断结果输出
Figure G200910262984XD00111

Claims (3)

1. 一种基于抗体浓度的齿轮故障诊断方法,其特征是包括以下步骤:
1)获取齿轮振动信号,通过振动传感设备,测试提取诊断齿轮正常工作和发生故障时的有关振动信号;
2)提取特征向量,采用小波包方法,提取正常工作和发生故障时的特征向量值,对每一信号经小波包分解,求各频带的能量,得到一组与信号相对应的能量序列{Ej|j=1,2,...,2N},由此确定与此能量序列对应的特征向量 P = ( E 1 , E 2 , . . . , E 2 N ) 。选择齿轮磨损、断齿、裂纹和正常状态进行故障特征提取,将选择的振动信号时域波形,用小波包方法进行分解,分别提取不同频段的特征信号,即由各频段的能量谱求得齿轮发生故障时其振动信号的特征向量值;
3)建立诊断系统,选择齿轮正常工作和不同故障类型发生时的特征向量作为训练样本,生成诊断系统,系统采用自适应免疫网络模型产生抗体群体,采用故障分类及其强度评价模型计算各种故障类型的故障强度系数,计算所得的齿轮不同故障状态下的故障强度系数作为具体诊断时的故障判断标准;
4)输入诊断数据,将第一步检测并通过第二步特征提取所获得的待诊断齿轮故障特征向量值作为抗原,输入到建立的故障诊断系统,系统发生免疫应答,产生相应的故障特征抗体群体,然后,计算不同故障类型的强度系数;
5)输出诊断结果,根据强度系数判断出可能存在的故障类型并输出相应结果。
2.如权利要求1所述的基于抗体浓度的齿轮故障诊断方法,其特征是所述步骤3)中描述的自适应免疫网络模型的算法描述如下:
①将待处理的数据定义为抗原集合AG={ag1,ag2,...,agI},抗体识别抗原的过程就是免疫应答识别待处理数据模式的过程。通过免疫应答识别过程,获得的记忆细胞集合M即为需要求解的模式。设系统中的抗原具体表示为:
AG = ( a g i . . . a g I ) T = a g 11 K a g 1 m M O M a g I 1 M a g Im - - - ( 1 )
agij∈[0,1]m,i=1,...,I;j=1,...,m。式(1)中:agi表示第i个抗原;agij表示第i个抗原的第j个向量。通过随机产生抗体:
AB = ( a b 1 . . . a b k ) T = a b 11 K a b 1 m M O M a b k 1 L a b km - - - ( 2 )
abij∈[0,1]m,i=1,...,k;j=1,...,m。作为新生成的抗体。式中:abi表示第i个抗体,abij表示第i个抗体的第j个向量。
②把抗原与抗体间的欧氏距离作为衡量二者之间匹配程度的一个指标。将抗原中的每个元素分别与抗体中的每个元素相匹配,对任意抗原agi和任意抗体abj,二者之间的距离d(i,j)可通过下式计算:
( i , j ) = ( &Sigma; t = 1 m ( a g it - a b jt ) 2 ) 1 / 2 - - - ( 3 )
对任意抗体abi和任意抗体abj,二者之间的距离s(i,j)可通过下式计算:
s ( i , j ) = ( &Sigma; t = 1 m ( a b it - a b jt ) 2 ) 1 / 2 - - - ( 4 )
对任意abj(abj∈AB,j=1,L,k),所受到的来自抗原agi(agi∈AG,i=1,L,I)的激励sl(i,j)为:
sl ( i , j ) = 1 / d ( i , j ) d ( i , j ) < T sl 0 other - - - ( 5 )
i=1,2,...,I;j=1,2,..,k。在式(5)中:
T sl = k sl &Sigma; i = 1 I &Sigma; j = 1 k d ( i , j ) I &times; k 2 - - - ( 6 )
Tsl表示抗原对抗体激励的阈值,其中:ksl∈(0,1]且为常数。
③当抗体abj所受到来自抗原的刺激达到一定程度(d(i,j)<Tsl时,免疫系统被激活,系统开始对抗体进行克隆增殖。增殖是指免疫细胞发生无性、有丝分裂。自我复制(没有交叉)的过程。抗体克隆的数目与其所受到的激励水平成正比,受激励越大的抗体克隆的数目越多;受激励越小的抗体克隆的数目越少。由于受抗原agi激励所导致的抗体abj的克隆数目为:
Figure F200910262984XC00033
式(7)中:
Figure F200910262984XC00034
表示取与sl(i,j)最接近的整数。抗体abj经agi刺激产生克隆的集合为Abji,抗体abj受抗原群体刺激所产生的克隆集合Cj为:
C j = A b j 1 &cup; A b j 2 &cup; L &cup; A b j c 1 ( i , j ) = U t = 1 c 1 ( i , j ) A b jt - - - ( 8 )
相应的抗体集合AB经克隆过程后新生成的抗体集合C为:
C = C 1 &cup; C 2 &cup; L &cup; C t = U t = 1 k C t - - - ( 9 )
克隆细胞要经历突变(伴随着高变异率)过程,从而使得新产生的细胞与所选择的抗原以更高的亲和力相匹配。在变异过程中每种免疫细胞变异的速率与细胞受体受抗原的激励程度的大小成反比:其受抗原的激励越大,变异越小;其受抗原的激励越小,变异越大。定义:
Mut(i,j)=d(i,j)                                (10)
为抗体abj受到抗原agi刺激后新产生的抗体的变异率。由抗原agi入侵引起的abj的克隆所生成的cj变异为:
c′j=cj+mut(i,j)×(cj-agi)                      (11)
抗体abj被抗愿群刺激产生的变异集合为C′j,变异产生的全部抗体集合为:
C &prime; = C 1 &prime; &cup; C 2 &prime; &cup; L &cup; C t &prime; = U t = 1 k C t &prime; - - - ( 12 )
变异保证了新产生的抗体具有充足的多样性,这是免疫系统强化学习机制的内在的特性,系统通过重复进行上述工作来增强自己识别抗原的性能。
④克隆变异之后首先进行克隆抑制(Clone suppression)。对集合C′,再次计算抗原与克隆变异后的抗体细胞的亲和力,计算克隆变异后的抗体细胞之间的亲和力。同时计算抗原-抗体抑制阈值Td和抗体-抗体抑制阈值Ts
T d = k d 1 p &Sigma; j = 1 p min i = 1 I d ( i , j ) - - - ( 13 )
式(13)中:kd∈(0,1]且为常数
T s = k s &Sigma; i = 1 p &Sigma; j = 1 p s ( i , j ) p &times; ( p - 1 ) 2 - - - ( 14 )
式(14)中:ks∈(0,1]且为常数
在克隆抑制阶段,首先删除那些与抗原亲和力大(于阈值Td)的抗体细胞,然后再根据抗体细胞之间的亲和力,清除亲和力小(于阈值Ts)的记忆细胞,此阶段完成后形成的抗体集合为C″,其元素满足以下条件:
a b j = a b j &Exists; a g i d ( i , j ) < T d 0 other - - - ( 15 )
a b j = a b j &Exists; a b i s ( i , j ) > T s 0 other - - ( 16 )
通过克隆抑制阶段,形成的抗体集合C″比C′的规模缩减。随后再进行克隆选择,对集合C″的抗体,按照其与抗原的亲和力大小,选择ξ%克隆变异后的细胞作为该抗体的记忆细胞,ξ%可以预先设定。将选择出的变异克隆细胞放入记忆细胞集合M中。
⑤经克隆变异后的抗体中那些具有较高抗原刺激度的抗体转化为具有较长寿命的记忆细胞,而那些具有较低抗原刺激度的抗体则随着抗原的消亡而消亡。由于抗体细胞之间的互相识别易引起自身免疫应答反应,因此,为了消除自身免疫,在每次加入新的抗体记忆细胞后,需要对那些能够自我识别的抗体进行抑制。
设定记忆细胞集合M中元素的个数为p1,定义记忆细胞亲和力成熟阈值为:
T s &prime; = k s &prime; &Sigma; i = 1 p 1 &Sigma; j = 1 p 1 s ( i , j ) p 1 &times; ( p 1 - 1 ) 2 - - - ( 17 )
式(17)中:k′s∈(0,1]且为常数
那么抗体经亲和力成熟过程后转化后,记忆细胞集合M中的元素应该满足以下条件:
a b j = a b j s ( i , j ) > T s &prime; 0 other - - - ( 18 )
式(18)中:i=1,2,...,p1;j=1,2,...,p1
⑥网络训练终止条件判断,任给一个抗体agi,在记忆细胞集合里都可以找到一个抗体abj,判断是否满足:d(i,j)≤T′s。如果满足,则训练结束,否则转步骤③继续。
3.如权利要求1所述的基于抗体浓度的齿轮故障诊断方法,其特征是上述步骤3)中所述的故障分类及其强度评价模型为:利用自适应免疫网络模型,将齿轮故障作为抗原,建立抗体种类及其浓度与齿轮故障种类及其强度关系模型,通过对模型中抗体种类和抗体浓度的统计,判断齿轮故障是否存在和故障出现时故障的种类;
设定ni为对应于齿轮的第i类故障发生时所产生的抗体数量,ci为齿轮正常状态下系统所产生的抗体数量,当齿轮的第i类故障发生时,定义其故障强度系数为:
r i = 1 - 1 1 + lg ( | n i - c i | + 1 ) - - - ( 19 )
由本公式的定义可知:0≤ri≤1,当ri=0时,表明齿轮当前工作状态下没有第i类故障,当ri=1时,表明第i类故障存在,当0<ri<1时,存在第i类故障的可能性,且ri值越大,表明当前齿轮面临的第i类故障可能性越高,在模型中,当没有故障存在时,系统各种抗体数目相对稳定;当故障存在时各抗体的数目急剧增加,因此,可以根据故障存在时刺激所产生的抗体种类及其浓度来识别故障种类及强度。
CN200910262984XA 2009-12-14 2009-12-14 一种基于抗体浓度的齿轮故障诊断方法 Expired - Fee Related CN101738314B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910262984XA CN101738314B (zh) 2009-12-14 2009-12-14 一种基于抗体浓度的齿轮故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910262984XA CN101738314B (zh) 2009-12-14 2009-12-14 一种基于抗体浓度的齿轮故障诊断方法

Publications (2)

Publication Number Publication Date
CN101738314A true CN101738314A (zh) 2010-06-16
CN101738314B CN101738314B (zh) 2011-08-17

Family

ID=42462074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910262984XA Expired - Fee Related CN101738314B (zh) 2009-12-14 2009-12-14 一种基于抗体浓度的齿轮故障诊断方法

Country Status (1)

Country Link
CN (1) CN101738314B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937522A (zh) * 2012-08-30 2013-02-20 桂林电子科技大学 一种齿轮箱复合故障诊断方法及系统
CN103310224A (zh) * 2013-06-08 2013-09-18 上海电机学院 一种无监督的人工免疫分类方法
CN103940608A (zh) * 2014-04-29 2014-07-23 中能电力科技开发有限公司 一种提高风电机组齿轮箱故障等级判断精度的方法
CN104122086A (zh) * 2014-08-11 2014-10-29 长沙理工大学 基于Kriging模型的风力机齿轮箱故障诊断方法
CN106092564A (zh) * 2016-06-06 2016-11-09 电子科技大学 基于esmd和能量算子解调的齿轮故障诊断方法
CN106198000A (zh) * 2016-07-11 2016-12-07 太原理工大学 一种采煤机摇臂齿轮故障诊断方法
CN107168054A (zh) * 2017-05-10 2017-09-15 沈阳工业大学 多机器人任务分配及路径规划方法
CN108896932A (zh) * 2018-05-29 2018-11-27 南京理工大学 基于小波包分解和csa的开关电源电路故障诊断方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19919604B4 (de) * 1999-04-29 2009-08-13 Olaf Unbehaun Verfahren und Vorrichtung zur Ermittlung von im Betrieb autretenden Fehlern an Rädern von Eisenbahnfahrzeugen
CN101266659B (zh) * 2008-05-08 2010-06-09 山东大学 基于免疫自适应遗传算法的机器人栅格子地图融合方法
CN101299007A (zh) * 2008-06-25 2008-11-05 茂名学院 组合式旋转机械复合故障诊断试验装置
CN100585366C (zh) * 2008-07-10 2010-01-27 武汉科技大学 一种同步齿轮箱初始故障的在线监测与诊断装置

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937522A (zh) * 2012-08-30 2013-02-20 桂林电子科技大学 一种齿轮箱复合故障诊断方法及系统
CN103310224A (zh) * 2013-06-08 2013-09-18 上海电机学院 一种无监督的人工免疫分类方法
CN103940608A (zh) * 2014-04-29 2014-07-23 中能电力科技开发有限公司 一种提高风电机组齿轮箱故障等级判断精度的方法
CN103940608B (zh) * 2014-04-29 2016-10-19 中能电力科技开发有限公司 一种提高风电机组齿轮箱故障等级判断精度的方法
CN104122086A (zh) * 2014-08-11 2014-10-29 长沙理工大学 基于Kriging模型的风力机齿轮箱故障诊断方法
CN104122086B (zh) * 2014-08-11 2016-05-04 长沙理工大学 基于Kriging模型的风力机齿轮箱故障诊断方法
CN106092564A (zh) * 2016-06-06 2016-11-09 电子科技大学 基于esmd和能量算子解调的齿轮故障诊断方法
CN106198000A (zh) * 2016-07-11 2016-12-07 太原理工大学 一种采煤机摇臂齿轮故障诊断方法
CN107168054A (zh) * 2017-05-10 2017-09-15 沈阳工业大学 多机器人任务分配及路径规划方法
CN108896932A (zh) * 2018-05-29 2018-11-27 南京理工大学 基于小波包分解和csa的开关电源电路故障诊断方法

Also Published As

Publication number Publication date
CN101738314B (zh) 2011-08-17

Similar Documents

Publication Publication Date Title
CN101738314B (zh) 一种基于抗体浓度的齿轮故障诊断方法
CN106323636A (zh) 栈式稀疏自动编码深度神经网络的机械故障程度特征自适应提取与诊断方法
CN108805059B (zh) 稀疏正则化滤波与自适应稀疏分解的齿轮箱故障诊断方法
CN109100143A (zh) 基于ceemdan和cfsfdp的滚动轴承故障诊断方法及设备
CN103439110B (zh) 滚动轴承早期微弱故障诊断方法
CN103955750B (zh) 基于特征融合和粒子滤波的滚动轴承剩余寿命预测方法
CN100485342C (zh) 机械故障的集成支持向量机混合智能诊断方法
CN108343599B (zh) 一种基于多粒度级联森林的水泵机组智能故障诊断方法
CN110674892A (zh) 一种基于加权多特征融合和svm分类的故障特征筛选方法
CN106092578A (zh) 一种基于小波包和支持向量机的机床主轴轴承装配状态在线检测方法
CN102998118B (zh) 一种基于形态学滤波和复杂度测度的轴承定量诊断方法
CN104330255A (zh) 一种基于多传感器信息融合的齿轮故障诊断方法
CN105547465A (zh) 一种变压器振动信号绕组状态特征提取方法
CN107817106A (zh) 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法
CN104502126A (zh) 一种基于模态区间的高速列车转向架故障诊断方法
CN106934421A (zh) 基于2dpca及svm的变电站电力变压器故障检测系统及检测方法
CN105424366A (zh) 基于eemd自适应消噪的轴承故障诊断方法
CN110175369A (zh) 一种基于二维卷积神经网络的齿轮剩余寿命预测方法
CN103968933A (zh) 基于模糊域特征的光纤振动信号识别方法
CN105862935A (zh) 一种用于挡土墙结构系统的损伤识别方法
CN107003665A (zh) 用于评估连接到电马达的旋转机械的状态的方法
CN111581762A (zh) 早期故障诊断方法及系统
CN107655690A (zh) 基于人工神经网络的电机轴承故障诊断方法
CN114112398A (zh) 一种变速工况下滚动轴承故障诊断方法
CN103983692B (zh) 一种用于挡土墙结构系统的损伤预警方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110817

Termination date: 20161214