CN113188794A - 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 - Google Patents
一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 Download PDFInfo
- Publication number
- CN113188794A CN113188794A CN202110472851.6A CN202110472851A CN113188794A CN 113188794 A CN113188794 A CN 113188794A CN 202110472851 A CN202110472851 A CN 202110472851A CN 113188794 A CN113188794 A CN 113188794A
- Authority
- CN
- China
- Prior art keywords
- data
- neural network
- formula
- fault diagnosis
- value
- 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
Images
Classifications
-
- 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/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- 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/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- 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/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种基于改进PSO‑BP神经网络的齿轮箱故障诊断方法及装置,故障诊断方法具体包括以下步骤:(1)采集齿轮箱正常状态的振动信号以及不同故障状态下的振动信号;(2)对数据样本进行打标签;(3)通过归一化方法对数据进行归一化处理;(4)采用主成分分析法对数据进行降维处理;(5)采用改进后的PSO算法优化BP神经网络权值和阈值,基于优化后的BP神经网络建立齿轮箱故障诊断模型;(6)利用所述优化后的BP神经网络故障诊断模型进行训练,得到最终的BP神经网络故障诊断模型。本发明提供的故障诊断方法及装置有效的提高了齿轮箱故障诊断的准确率,对于提高安全生产效率、降低经济损失,具有重要的实际工程意义。
Description
技术领域
本发明涉及旋转机械故障诊断技术领域,具体涉及一种基于改进PSO-BP神经网络齿轮箱故障诊断方法及装置。
背景技术
齿轮箱作为机械设备的重要部件,在机械领域中有着广泛的应用,其运行状态正常与否直接决定了整台机械设备的运行情况。齿轮箱是一个结构复杂的系统,其中不仅包含起整体支撑和密封作用的箱体,还有传递动力和运动的齿轮、轴、轴承,以及其它必要的构件等。如果齿轮箱在运行过程中产生故障问题,不仅会给生产效益带来一定的损失,而且还会对人身安全造成伤害。因此,为了保障设备正常运行和保护生命财产安全,对齿轮箱的运行状态进行及时的评价具有重要的应用价值和意义。
传统齿轮箱故障诊断方法是分析收集到的振动信号的频谱,然后通过人工或机器模式识别判断是否有故障。然而,由于齿轮箱故障的原因与征兆之间关系存在不确定性、由于季节和环境等因素也会对正确诊断产生较大影响,导致齿轮箱故障错判、漏判率较高,对设备的正常运行造成较大安全隐患。因此,研究齿轮箱故障诊断方法,对于保障设备的安全运行至关重要。
近年来,由于人工智能技术的发展步伐加快,进而促进故障诊断技术也进入一个崭新的阶段。基于人工智能的齿轮箱故障诊断技术不需要故障原因与征兆之间精确的物理模型,可以很好地解决不确定性问题。因此,进一步深入研究基于人工智能的齿轮箱故障诊断方法,提高诊断效率和准确率,对于提高安全生产效率、降低经济损失,具有重要的实际工程意义。
发明内容
本发明提供了一种基于改进PSO-BP神经网络齿轮箱故障诊断方法及装置,将改进的粒子群算法(Particle Swarm Optimization,PSO)和BP神经网络相结合,提高齿轮箱故障诊断的准确率,弥补神经网络算法在故障诊断领域存在的不足。
为实现上述目的,本发明提供了如下的技术方案:
一种基于改进PSO-BP神经网络齿轮箱故障诊断方法,具体步骤如下:
步骤1000:采集齿轮箱正常状态的振动信号以及不同故障状态下的振动信号,形成第一数据集;
步骤2000:对所述第一数据集中的数据样本打标签,形成第二数据集;
步骤3000:通过归一化方法对所述第二数据集中的数据进行归一化处理,形成第三数据集;
步骤4000:采用主成分分析法对所述第三数据集中的数据进行降维处理,形成第四数据集;
步骤5000:采用改进后的PSO算法优化BP神经网络权值和阈值,基于优化后的BP神经网络建立齿轮箱故障诊断模型;
步骤6000:利用所述优化后的BP神经网络齿轮箱故障诊断模型进行训练,得到最终的BP神经网络故障诊断模型。
优选的,所述步骤2000中,采用0、1组合来表示正常工作状态和故障状态;
优选的,所述步骤3000中,数据归一化采用Min-Max归一化方法,归一化公式为:
公式(1)中,MaxValue表示样本数据的最大值;MinValue表示样本数据的最小值;x表示样本原始数据;y表示归一化之后的数据;
所述步骤4000中,采用主成分分析法对第三数据集中的数据进行降维处理,进一步包括以下步骤:
步骤4100:齿轮箱第三数据集中的n个样本,每个样本测得m个属性,将这n×m个数据构成矩阵X:
可知Xij(i=1,2,…,n;j=1,2,…,m)为第i个样本的第j个特征属性的取值;
步骤4200:根据公式(2)和(3)计算每个样本各个特征属性的均值和标准差:
步骤4300:根据公式(4)求解标准化矩阵:
步骤4400:根据公式(5)求解相关系数矩阵R=(rjk)m×m:
步骤4500:求解特征值及特征向量:由特征方程|λE-R|=0(其中E为单位矩阵)求出相关系数矩阵R的特征值由大到小排序λ1>λ2>...>λm以及对应特征向量为l1,l2,...lm;
公式(6)中zij是标准化矩阵Z的向量;
所述步骤5000中,进一步包括以下步骤:
步骤5010:确定BP神经网络的拓扑结构,该结构包括:神经网络输入层节点数m、隐含层节点数h和输出层节点数d;输入层节点数m取决于样本的特征个数;输出层节点数d取决于故障种类的数量;隐含层节点数h采用经验法和试凑法来确定合适的隐含层节点数;隐含层节点数h与输入层节点m、输出层节点数d满足公式(7)所示的函数关系:
公式(7)中,m为输入层节点数(由实际的输入个数决定);h为隐含层节点数;d为输出层节点数;a为随机数,取值范围为[1,10];
步骤5020:计算粒子个体编码长度,并把粒子向量映射为神经网络的权值和阈值,其中粒子个体编码长度计算公式如公式(8)所示:
Length=h×m+h+d×h+d (8)
公式(8)中,m为输入层节点数;h为隐含层节点数;d为输出层节点数;
步骤5030:粒子群种群的初始化,包括初始化粒子的速度、位置、学习因子和迭代次数等;
步骤5040:神经网络输入训练样本前向传播通过公式(9)计算各个粒子适应度值:
f(xi)=-MSE (9)
公式(9)中,训练误差MSE为神经网络的实际输出和期望输出的均方误差,MSE表达式如公式(10)所示:
步骤5050:如果该粒子当前的适应度值优于其历史最优值,用当前位置替代历史最优值,作为个体极值;
步骤5060:如果该粒子的历史最优值优于全局最优值,用该粒子的历史最优值替代全局最优值,作为全局极值;
步骤5070:迭代更新,根据公式(11)、公式(12)更新所有粒子位置和速度,并检查粒子速度和位置是否超出设置范围,如果超出范围,则用边界值作为粒子的速度和位置;
xi+1(t+1)=xi(t)+vi+1(t+1) (11)
公式(11)中,t表示当前迭代的次数;xi(t)表示第t次迭代第i个粒子位置;vi(t)表示第t次迭代第i个粒子速度;
vi+1(t+1)=wvi(t)+c1r1(pbesti(t)-xi(t))+c2r2(gbest-xi(t)) (12)
公式(12)中,vi(t)表示第t次迭代第i个粒子速度;pbesti(t)表示第i个粒子第t次迭代搜索到的最优位置;gbest表示整个粒子群搜索到的最优位置;c1表示学习因子,影响着粒子的“自我学习”能力;c2表示学习因子,影响着粒子的“社会学习”能力;r1、r2是介于[0,1]之间的随机数;w见公式(13);
公式(13)中,w为粒子群算法的惯性权重;wmax为惯性权重最大值;wmin为惯性权重最小值;T为最大迭代次数;t为当前迭代次数,取值范围在wmin和wmax之间;
步骤5080:根据各个粒子的适应度值确定粒子的个体历史最优和群体历史最优;
步骤5090:根据群体历史最优的变化对算法进行“早熟”判断,如果出现“早熟”转到步骤5100,如果没有则转到步骤5110;“早熟”判断是通过引入一个是否对粒子施加扰动的评判指标ε,如果连续多次迭代中,邻近两次粒子历史最优适应度值差的绝对值都小于ε,则判断该粒子陷入“早熟”,然后对粒子施加一次扰动;
步骤5100:利用公式(14)给陷入“早熟”的粒子施加一次扰动,并转到步骤5090,其中公式(14)如下:
其中,
公式(14)中,X为粒子位置;v为粒子速度;i为粒子编号;k为当前迭代次数;N为粒子群的空间维度;μ和α为常数;m为粒子个数,abs()为求绝对值函数;公式(15)中,rand()为[0,1]均匀分布的随机数,normrnd()为生成服从正态分布的随机数;
步骤5110:判断最小误差或最大迭代次数是否达到限定值;如果最小误差和最大迭代次数都没达到限定值,则转到步骤5070;如果最小误差没有达到要求,而粒子群迭代次数到达限定值,则转到步骤5120;
步骤5120:把粒子群优化的最优权值和阈值作为BP神经网络的初始权值和阈值,当实际输出与期望输出不符时,进入误差的反向传播阶段,误差通过输出层,以网络误差最小为目标函数,按照误差梯度下降的方式不断调整各层权值,当达到最小误差或神经网络的最大训练次数,选取此时的权值和阈值,算法结束。
一种基于改进PSO-BP神经网络齿轮箱故障诊断装置,具体包括:
数据采集模块:用于齿轮箱数据振动信号数据的采集,获得第一数据集;
数据处理模块:用于根据所述数据采集模块获得的第一数据信息进行数据预处理操作,具体包括对数据进行标签处理、数据归一化处理、主成分分析降维处理,最终获得第四数据集;
训练模块:用于建立基于改进PSO-BP神经网络故障诊断模型,通过所述第四数据集中的训练集训练基于改进PSO-BP神经网络模型,再通过所述第四数据集中的验证集验证所述基于改进PSO-BP神经网络故障诊断模型的准确性;
故障诊断模块:对齿轮箱振动实时检测数据进行预处理后得到数据样本,将所述数据样本输入训练后的基于改进PSO-BP神经网络故障诊断模型,最终得到齿轮箱故障诊断结果。
优选的,一种基于改进PSO-BP神经网络齿轮箱故障诊断装置,其特征在于,所述数据采集模块,通过在齿轮箱内安装多个传感器,通过多个传感器对齿轮箱齿轮不同的部位进行振动信号进行采样,每个传感器测得一种故障模式的一个特征,经过多次信号采集,获得第一数据集。
与现有技术相比,本发明的有益效果是:
本发明采用一种基于权重指数变化和粒子扰动策略的粒子群算法优化BP神经网络权值和阈值,解决了因神经网络算法全局寻找最优权值和阈值能力不强,而导致收敛速度慢、故障诊断准确率不高的问题,为齿轮箱故障诊断提供了一种基于人工智能的新方法,提高安全生产效率和降低经济损失。
附图说明
图1是基于改进PSO-BP神经网络齿轮箱故障诊断方法流程图;
具体实施方式
为了更清楚地理解本发明的上述方案,下面结合附图对本发明进行进一步的详细描述。需要说明的是,此处所描述的具体实施仅用于解释本申请,并不用于限定本申请。
图1是本发明一种基于改进PSO-BP神经网络齿轮箱故障诊断方法流程图,具体步骤如下:
步骤1000:采集齿轮箱正常状态的振动信号以及不同故障状态下的振动信号,形成第一数据集;
步骤2000:对所述第一数据集中的数据样本打标签,形成第二数据集;
步骤3000:通过归一化方法对所述第二数据集中的数据进行归一化处理,形成第三数据集;
步骤4000:采用主成分分析法对所述第三数据集中的数据进行降维处理,形成第四数据集;
步骤5000:采用改进后的PSO算法优化BP神经网络权值和阈值,基于优化后的BP神经网络建立齿轮箱故障诊断模型;
步骤6000:利用所述优化后的BP神经网络齿轮箱故障诊断模型进行训练,得到最终的BP神经网络故障诊断模型。
具体的,所述步骤1000中,所得第一数据集是一个n×m的矩阵,其中,n为所测数据的样本数,m是传感器个数;
具体的,所述步骤2000中,采用0、1组合来表示正常工作状态和故障状态;
具体的,所述步骤3000中,数据归一化采用Min-Max归一化方法,归一化公式为:
公式(1)中,MaxValue表示样本数据的最大值;MinValue表示样本数据的最小值;x表示样本原始数据;y表示归一化之后的数据;
所述步骤4000中,采用主成分分析法对第三数据集中的数据进行降维处理,进一步包括以下步骤:
步骤4100:齿轮箱第三数据集中的n个样本,每个样本测得m个属性,将这n×m个数据构成矩阵X:
可知Xij(i=1,2,…,n;j=1,2,…,m)为第i个样本的第j个特征属性的取值;
步骤4200:根据公式(2)和(3)计算每个样本各个特征属性的均值和标准差:
步骤4300:根据公式(4)求解标准化矩阵:
步骤4400:根据公式(5)求解相关系数矩阵R=(rjk)m×m:
步骤4500:求解特征值及特征向量:由特征方程|λE-R|=0(其中E为单位矩阵)求出相关系数矩阵R的特征值由大到小排序λ1>λ2>...>λm以及对应特征向量为l1,l2,...lm;
公式(6)中zij是标准化矩阵Z的向量;
具体的,所述步骤5000中,用改进之后的PSO优化BP神经网络就是用改进后的PSO算法对BP神经网络的权值和阈值进行寻优,进一步包括以下步骤:
步骤5010:确定BP神经网络的拓扑结构,该结构包括:神经网络输入层节点数m、隐含层节点数h和输出层节点数d;输入层节点数m取决于样本的特征个数;输出层节点数d取决于故障种类的数量;隐含层节点数h采用经验法和试凑法来确定合适的隐含层节点数;隐含层节点数h与输入层节点m、输出层节点数d满足公式(7)所示的函数关系:
公式(7)中,m为输入层节点数(由实际的输入个数决定);h为隐含层节点数;d为输出层节点数;a为随机数,取值范围为[1,10];
步骤5020:计算粒子个体编码长度,并把粒子向量映射为神经网络的权值和阈值,其中粒子个体编码长度计算公式如公式(8)所示:
Length=h×m+h+d×h+d (8)
公式(8)中,m为输入层节点数;h为隐含层节点数;d为输出层节点数;
步骤5030:粒子群种群的初始化,包括初始化粒子的速度、位置、学习因子和迭代次数等;
步骤5040:神经网络输入训练样本前向传播通过公式(9)计算各个粒子适应度值:
f(xi)=-MSE (9)
公式(9)中,训练误差MSE为神经网络的实际输出和期望输出的均方误差,MSE表达式如公式(10)所示:
步骤5050:如果该粒子当前的适应度值优于其历史最优值,用当前位置替代历史最优值,作为个体极值;
步骤5060:如果该粒子的历史最优值优于全局最优值,用该粒子的历史最优值替代全局最优值,作为全局极值;
步骤5070:迭代更新,根据公式(11)、公式(12)更新所有粒子位置和速度,并检查粒子速度和位置是否超出设置范围,如果超出范围,则用边界值作为粒子的速度和位置;
xi+1(t+1)=xi(t)+vi+1(t+1) (11)
公式(11)中,t表示当前迭代的次数;xi(t)表示第t次迭代第i个粒子位置;vi(t)表示第t次迭代第i个粒子速度;
vi+1(t+1)=wvi(t)+c1r1(pbesti(t)-xi(t))+c2r2(gbest-xi(t)) (12)
公式(12)中,vi(t)表示第t次迭代第i个粒子速度;pbesti(t)表示第i个粒子第t次迭代搜索到的最优位置;gbest表示整个粒子群搜索到的最优位置;c1表示学习因子,影响着粒子的“自我学习”能力;c2表示学习因子,影响着粒子的“社会学习”能力;r1、r2是介于[0,1]之间的随机数;w见公式(13);
公式(13)中,w为粒子群算法的惯性权重;wmax为惯性权重最大值;wmin为惯性权重最小值;T为最大迭代次数;t为当前迭代次数,取值范围在wmin和wmax之间;
步骤5080:根据各个粒子的适应度值确定粒子的个体历史最优和群体历史最优;
步骤5090:根据群体历史最优的变化对算法进行“早熟”判断,如果出现“早熟”转到步骤5100,如果没有则转到步骤5110;“早熟”判断是通过引入一个是否对粒子施加扰动的评判指标ε,如果连续多次迭代中,邻近两次粒子历史最优适应度值差的绝对值都小于ε,则判断该粒子陷入“早熟”,然后对粒子施加一次扰动;
步骤5100:利用公式(14)给陷入“早熟”的粒子施加一次扰动,并转到步骤5090,其中公式(14)如下:
其中,
公式(14)中,X为粒子位置;v为粒子速度;i为粒子编号;k为当前迭代次数;N为粒子群的空间维度;μ和α为常数;m为粒子个数,abs()为求绝对值函数;公式(15)中,rand()为[0,1]均匀分布的随机数,normrnd()为生成服从正态分布的随机数;
步骤5110:判断最小误差或最大迭代次数是否达到限定值;如果最小误差和最大迭代次数都没达到限定值,则转到步骤5070;如果最小误差没有达到要求,而粒子群迭代次数到达限定值,则转到步骤5120;
步骤5120:把粒子群优化的最优权值和阈值作为BP神经网络的初始权值和阈值,当实际输出与期望输出不符时,进入误差的反向传播阶段,误差通过输出层,以网络误差最小为目标函数,按照误差梯度下降的方式不断调整各层权值,当达到最小误差或神经网络的最大训练次数,选取此时的权值和阈值,算法结束。
一种基于改进PSO-BP神经网络齿轮箱故障诊断装置,具体包括:
数据采集模块:用于齿轮箱数据振动信号数据的采集,获得第一数据集;
数据处理模块:用于根据所述数据采集模块获得的第一数据信息进行数据预处理操作,具体包括对数据进行标签处理、数据归一化处理、主成分分析降维处理,最终获得第四数据集;
训练模块:用于建立基于改进PSO-BP神经网络故障诊断模型,通过所述第四数据集中的训练集训练基于改进PSO-BP神经网络模型,再通过所述第四数据集中的验证集验证所述基于改进PSO-BP神经网络故障诊断模型的准确性;
故障诊断模块:对齿轮箱振动实时检测数据进行预处理后得到数据样本,将所述数据样本输入训练后的基于改进PSO-BP神经网络故障诊断模型,最终得到齿轮箱故障诊断结果。
优选的,所述数据采集模块,通过在齿轮箱内安装多个传感器,通过多个传感器对齿轮箱齿轮不同的部位进行振动信号进行采样,每个传感器测得一种故障模式的一个特征,经过多次信号采集,获得第一数据集。
本发明提供的另一个实施例,采用QPZZ-II.旋转机械振动分析及故障诊断试验平台系统对齿轮箱数据进行采集,此系统可以对工作中此轮箱齿轮的振动信号进行提取和分析,进而得到反应齿轮状态的数据集合。数据的采集方法是用多个传感器对齿轮箱齿轮不同的部位进行采样,采样频率为2000×2.56hz,传感器的数量是9个。
9个传感器编号分别为CH1到CH8、TACH1,9个传感器后面的文字说明是每一个传感器测得一个故障模式下一个特征的含义,即一个故障数据集里的一个列数据所表示的含义:
1---TACH1 转速,光电,每周1个脉冲
2---CH1 输入轴X方向位移
3---CH2 Y方向位移
4---CH3 加速度输入轴电机侧轴承Y
5---CH4 加速度输出轴电机侧轴承Y
6---CH5 加速度输入轴负载侧轴承Y
7---CH6 加速度输出轴负载侧轴承Y
8---CH7 加速度输出轴负载侧轴承X
9---CH8 磁电式速度输出轴负载侧轴承X
所述的9个传感器测得齿轮的四种状态模式数据,每个传感器测得一种故障模式的一个特征,因此每个故障模式数据由9个特征数据组成,原始数据集共有53248个样本数据,因此每个故障模式的原始数据用53248×9的数据矩阵表示。
本实施例诊断的四种齿轮故障(包含正常模式)包括:正常模式、磨损故障、断齿故障、点蚀故障。
对故障诊断系统进行训练时,需要把正常、磨损、点蚀、断齿四种状态的数据集矩阵合在一起,样本总数共有212992个,用数据矩阵表示即212992行9列,形成第一数据集;
根据所述步骤2000对第一数据集打标签,用于识别数据所述的故障类型,形成第二数据集;在本实施例中,采用4位0、1组合来表示四种不同的故障,故障类型所对应的标签见表1。
表1故障类型对应的标签
根据所述步骤3000对第二数据集进行归一化处理,获得第三数据集。归一化处理之后的第三数据集中所有数据分布在[0,1]区间内;
根据所述步骤4000对第三数据集进行降维处理,有5个特征累计方差贡献率超过90%,即降维后的数据集由原来的9维降到5维,形成第四数据集;
根据所述步骤5000对第四数据集进行训练。各模型参数表见表2。
表2参数表
对比了BP网络模型、PSO-BP网络模型以及本发明公开的改进PSO-BP网络模型这三种故障诊断模型的诊断效果,从训练误差、迭代次数、诊断正确率这三个方面作对比,这三种模型的性能见表3。
表3三种诊断模型性能对比
从表3中三种故障诊断模型的结果对比可知,在系统学习收敛效率方面,改进PSO-BP神经网络诊断模型在第146次迭代之后,诊断模型已经收敛到设置的精度10-3。PSO-BP神经网络诊断模型用了251次迭代,系统的收敛精度只到10-2.5,在这之后系统的收敛精度基本不变。BP神经网络诊断模型收敛效果最差,经过401次迭代收敛精度为10-1,之后精度也不再变化。最终,BP模型和PSO-BP模型都未达到设置的精度10-3。在诊断效果方面,改进PSO-BP网络模型的诊断正确率为77.0%,PSO-BP网络模型的诊断正确率为63.5%,而BP网络模型的诊断正确率为52.0%。由以上可知,本发明提出的故障诊断模型在系统本身的学习效率上和故障分类的精度上与改进之前的算法相比都具有很大的提高。
本发明公开了一种基于改进PSO-BP神经网络齿轮箱故障诊断模型,解决了因神经网络算法全局寻找最优权值和阈值能力不强,而导致收敛速度慢和故障诊断准确率不高的问题,为齿轮箱故障诊断提供了一种基于人工智能的新方法,有效提高生产效率和降低经济损失。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所做的任何修改、同等替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于改进PSO-BP神经网络齿轮箱故障诊断方法,其特征在于,包括以下步骤:
步骤1000:采集齿轮箱正常状态的振动信号以及不同故障状态下的振动信号,形成第一数据集;
步骤2000:对所述第一数据集中的数据样本打标签,形成第二数据集;
步骤3000:通过归一化方法对所述第二数据集中的数据进行归一化处理,形成第三数据集;
步骤4000:采用主成分分析法对所述第三数据集中的数据进行降维处理,形成第四数据集;
步骤5000:采用改进后的PSO算法优化BP神经网络权值和阈值,基于优化后的BP神经网络建立齿轮箱故障诊断模型;
步骤6000:利用所述优化后的BP神经网络齿轮箱故障诊断模型进行训练,得到最终的BP神经网络故障诊断模型。
2.根据权利要求1所述的基于改进PSO-BP神经网络齿轮箱故障诊断方法,其特征在于,
所述步骤1000中,所得第一数据集是一个n×m的矩阵,其中,n为所测数据的样本数,m是传感器个数;
所述步骤2000中,采用0、1组合来表示正常工作状态和故障状态;
所述步骤3000中,数据归一化采用Min-Max归一化方法,归一化公式为:
公式(1)中,MaxValue表示样本数据的最大值;MinValue表示样本数据的最小值;x表示样本原始数据;y表示归一化之后的数据;
所述步骤4000中,采用主成分分析法对第三数据集中的数据进行降维处理,进一步包括以下步骤:
步骤4100:齿轮箱第三数据集中的n个样本,每个样本测得m个属性,将这n×m个数据构成矩阵X:
可知Xij(i=1,2,…,n;j=1,2,…,m)为第i个样本的第j个特征属性的取值;
步骤4200:根据公式(2)和(3)计算每个样本各个特征属性的均值和标准差:
步骤4300:根据公式(4)求解标准化矩阵:
步骤4400:根据公式(5)求解相关系数矩阵R=(rjk)m×m:
步骤4500:求解特征值及特征向量:由特征方程|λE-R|=0(其中E为单位矩阵)求出相关系数矩阵R的特征值由大到小排序λ1>λ2>...>λm以及对应特征向量为l1,l2,...lm;
公式(6)中zij是标准化矩阵Z的向量;
所述步骤5000中,用改进PSO算法优化BP神经网络的权值和阈值,进一步包括以下步骤:
步骤5010:确定BP神经网络的拓扑结构,该结构包括:神经网络输入层节点数m、隐含层节点数h和输出层节点数d;输入层节点数m取决于样本的特征个数;输出层节点数d取决于故障种类的数量;隐含层节点数h采用经验法和试凑法来确定合适的隐含层节点数;隐含层节点数h与输入层节点m、输出层节点数d满足公式(7)所示的函数关系:
公式(7)中,m为输入层节点数(由实际的输入个数决定);h为隐含层节点数;d为输出层节点数;a为随机数,取值范围为[1,10];
步骤5020:计算粒子个体编码长度,并把粒子向量映射为神经网络的权值和阈值,其中粒子个体编码长度计算公式如公式(8)所示:
Length=h×m+h+d×h+d (8)
公式(8)中,m为输入层节点数;h为隐含层节点数;d为输出层节点数;
步骤5030:粒子群种群的初始化,包括初始化粒子的速度、位置、学习因子和迭代次数等;
步骤5040:神经网络输入训练样本前向传播通过公式(9)计算各个粒子适应度值:
f(xi)=-MSE (9)
公式(9)中,训练误差MSE为神经网络的实际输出和期望输出的均方误差,MSE表达式如公式(10)所示:
步骤5050:如果该粒子当前的适应度值优于其历史最优值,用当前位置替代历史最优值,作为个体极值;
步骤5060:如果该粒子的历史最优值优于全局最优值,用该粒子的历史最优值替代全局最优值,作为全局极值;
步骤5070:迭代更新,根据公式(11)、公式(12)更新所有粒子位置和速度,并检查粒子速度和位置是否超出设置范围,如果超出范围,则用边界值作为粒子的速度和位置;
xi+1(t+1)=xi(t)+vi+1(t+1) (11)
公式(11)中,t表示当前迭代的次数;xi(t)表示第t次迭代第i个粒子位置;vi(t)表示第t次迭代第i个粒子速度;
vi+1(t+1)=wvi(t)+c1r1(pbesti(t)-xi(t))+c2r2(gbest-xi(t)) (12)
公式(12)中,vi(t)表示第t次迭代第i个粒子速度;pbesti(t)表示第i个粒子第t次迭代搜索到的最优位置;gbest表示整个粒子群搜索到的最优位置;c1表示学习因子,影响着粒子的“自我学习”能力;c2表示学习因子,影响着粒子的“社会学习”能力;r1、r2是介于[0,1]之间的随机数;w见公式(13);
公式(13)中,w为粒子群算法的惯性权重;wmax为惯性权重最大值;wmin为惯性权重最小值;T为最大迭代次数;t为当前迭代次数,取值范围在wmin和wmax之间;
步骤5080:根据各个粒子的适应度值确定粒子的个体历史最优和群体历史最优;
步骤5090:根据群体历史最优的变化对算法进行“早熟”判断,如果出现“早熟”转到步骤5100,如果没有则转到步骤5110;“早熟”判断是通过引入一个是否对粒子施加扰动的评判指标ε,如果连续多次迭代中,邻近两次粒子历史最优适应度值差的绝对值都小于ε,则判断该粒子陷入“早熟”,然后对粒子施加一次扰动;
步骤5100:利用公式(14)给陷入“早熟”的粒子施加一次扰动,并转到步骤5090,其中公式(14)如下:
其中,
公式(14)中,X为粒子位置;v为粒子速度;i为粒子编号;k为当前迭代次数;N为粒子群的空间维度;μ和α为常数;m为粒子个数,abs()为求绝对值函数;公式(15)中,rand()为[0,1]均匀分布的随机数,normrnd()为生成服从正态分布的随机数;
步骤5110:判断最小误差或最大迭代次数是否达到限定值;如果最小误差和最大迭代次数都没达到限定值,则转到步骤5070;如果最小误差没有达到要求,而粒子群迭代次数到达限定值,则转到步骤5120;
步骤5120:把粒子群优化的最优权值和阈值作为BP神经网络的初始权值和阈值,当实际输出与期望输出不符时,进入误差的反向传播阶段,误差通过输出层,以网络误差最小为目标函数,按照误差梯度下降的方式不断调整各层权值,当达到最小误差或神经网络的最大训练次数,选取此时的权值和阈值,算法结束。
3.一种采用如权利要求书1-2任一项所述的一种基于改进PSO-BP神经网络齿轮箱故障诊断方法的故障诊断装置,其特征在于,所述装置包括:
数据采集模块:用于齿轮箱数据振动信号数据的采集,获得第一数据集;
数据处理模块:用于根据所述数据采集模块获得的第一数据信息进行数据预处理操作,具体包括对数据进行标签、数据归一化处理、主成分分析降维处理,最终获得第四数据集;
训练模块:用于建立基于改进PSO-BP神经网络故障诊断模型,通过所述第四数据集中的训练集训练基于改进PSO-BP神经网络模型,再通过所述第四数据集中的验证集验证所述基于改进PSO-BP神经网络故障诊断模型的准确性;
故障诊断模块:对齿轮箱振动实时检测数据进行预处理后得到数据样本,将所述数据样本输入训练后的基于改进PSO-BP神经网络故障诊断模型,最终得到齿轮箱故障诊断结果。
4.根据权利要求3所述的一种基于改进PSO-BP神经网络齿轮箱故障诊断装置,其特征在于,所述数据采集模块,通过在齿轮箱内安装多个传感器,通过多个传感器对齿轮箱齿轮不同的部位进行振动信号进行采样,每个传感器测得一种故障模式的一个特征,经过多次信号采集,获得第一数据集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110472851.6A CN113188794B (zh) | 2021-04-29 | 2021-04-29 | 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110472851.6A CN113188794B (zh) | 2021-04-29 | 2021-04-29 | 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113188794A true CN113188794A (zh) | 2021-07-30 |
CN113188794B CN113188794B (zh) | 2021-11-23 |
Family
ID=76980498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110472851.6A Active CN113188794B (zh) | 2021-04-29 | 2021-04-29 | 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113188794B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113567117A (zh) * | 2021-09-24 | 2021-10-29 | 南京农业大学 | 基于psoobp-cs算法的齿轮箱故障诊断方法 |
CN113742638A (zh) * | 2021-08-30 | 2021-12-03 | 南通大学 | 基于峭度的FastICA和逼近求解域的STLBO电机轴承故障诊断方法 |
CN113923104A (zh) * | 2021-12-07 | 2022-01-11 | 南京信息工程大学 | 基于小波神经网络的网络故障诊断方法、设备及存储介质 |
CN117171659A (zh) * | 2023-10-31 | 2023-12-05 | 国网浙江省电力有限公司宁波供电公司 | 一种智能化配变高压侧故障诊断方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221213A (zh) * | 2008-01-25 | 2008-07-16 | 湖南大学 | 基于粒子群算法的模拟电路故障诊断神经网络方法 |
CN103410662A (zh) * | 2013-08-06 | 2013-11-27 | 江苏科技大学 | 风电系统最大风能捕获神经网络补偿控制方法 |
CN104638671A (zh) * | 2015-01-16 | 2015-05-20 | 合肥工业大学 | 基于引入粒子群算法及相位扰动的孤岛检测方法 |
US20180089566A1 (en) * | 2015-05-22 | 2018-03-29 | Huawei Technologies Co., Ltd. | Method and apparatus for positioning of artificial neural network |
CN108182500A (zh) * | 2018-02-02 | 2018-06-19 | 沈阳理工大学 | 基于加速寿命试验的弹药贮存可靠性预测方法 |
CN109344517A (zh) * | 2018-10-14 | 2019-02-15 | 毛述春 | 一种新能源汽车的高压绝缘故障诊断方法 |
CN110533484A (zh) * | 2019-09-05 | 2019-12-03 | 四川长虹电器股份有限公司 | 一种基于pca和改进bp神经网络的产品销量预测方法 |
US20200241487A1 (en) * | 2019-01-24 | 2020-07-30 | Qingdao university of technology | Networked control system time-delay compensation method based on predictive control |
CN112504682A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军63963部队 | 基于粒子群优化算法的底盘发动机故障诊断方法及系统 |
CN112529285A (zh) * | 2020-12-08 | 2021-03-19 | 湖南工业大学 | 一种基于相似日理论和pca-pso-bp的光伏发电功率预测方法 |
-
2021
- 2021-04-29 CN CN202110472851.6A patent/CN113188794B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221213A (zh) * | 2008-01-25 | 2008-07-16 | 湖南大学 | 基于粒子群算法的模拟电路故障诊断神经网络方法 |
CN103410662A (zh) * | 2013-08-06 | 2013-11-27 | 江苏科技大学 | 风电系统最大风能捕获神经网络补偿控制方法 |
CN104638671A (zh) * | 2015-01-16 | 2015-05-20 | 合肥工业大学 | 基于引入粒子群算法及相位扰动的孤岛检测方法 |
US20180089566A1 (en) * | 2015-05-22 | 2018-03-29 | Huawei Technologies Co., Ltd. | Method and apparatus for positioning of artificial neural network |
CN108182500A (zh) * | 2018-02-02 | 2018-06-19 | 沈阳理工大学 | 基于加速寿命试验的弹药贮存可靠性预测方法 |
CN109344517A (zh) * | 2018-10-14 | 2019-02-15 | 毛述春 | 一种新能源汽车的高压绝缘故障诊断方法 |
US20200241487A1 (en) * | 2019-01-24 | 2020-07-30 | Qingdao university of technology | Networked control system time-delay compensation method based on predictive control |
CN110533484A (zh) * | 2019-09-05 | 2019-12-03 | 四川长虹电器股份有限公司 | 一种基于pca和改进bp神经网络的产品销量预测方法 |
CN112529285A (zh) * | 2020-12-08 | 2021-03-19 | 湖南工业大学 | 一种基于相似日理论和pca-pso-bp的光伏发电功率预测方法 |
CN112504682A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军63963部队 | 基于粒子群优化算法的底盘发动机故障诊断方法及系统 |
Non-Patent Citations (6)
Title |
---|
LAN SHI 等: "PCA-PSO-BP Neural Network Application in IDS", 《INTERNATIONAL POWER, ELECTRONICS AND MATERIALS ENGINEERING CONFERENCE》 * |
ZHAO LIU 等: "Midterm Power Load Forecasting Model Based on Kernel Principal Component Analysis and Back Propagation Neural Network with Particle Swarm Optimization", 《BIG DATA》 * |
任少伟: "基于PCA的PSO-DE混合算法优化BP神经网络在煤与瓦斯突出预测中的应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 * |
张义宏: "基于PCA的BP神经网络优化的研究与应用", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 * |
曹伟: "基于改进型PSO-BP神经网络算法的水环境质量评价", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
洪文鹏、陈重: "基于自适应粒子群优化BP神经网络的氨法烟气脱硫效率预测", 《动力工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113742638A (zh) * | 2021-08-30 | 2021-12-03 | 南通大学 | 基于峭度的FastICA和逼近求解域的STLBO电机轴承故障诊断方法 |
CN113567117A (zh) * | 2021-09-24 | 2021-10-29 | 南京农业大学 | 基于psoobp-cs算法的齿轮箱故障诊断方法 |
CN113923104A (zh) * | 2021-12-07 | 2022-01-11 | 南京信息工程大学 | 基于小波神经网络的网络故障诊断方法、设备及存储介质 |
CN113923104B (zh) * | 2021-12-07 | 2022-03-11 | 南京信息工程大学 | 基于小波神经网络的网络故障诊断方法、设备及存储介质 |
CN117171659A (zh) * | 2023-10-31 | 2023-12-05 | 国网浙江省电力有限公司宁波供电公司 | 一种智能化配变高压侧故障诊断方法 |
CN117171659B (zh) * | 2023-10-31 | 2024-01-26 | 国网浙江省电力有限公司宁波供电公司 | 一种智能化配变高压侧故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113188794B (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113188794B (zh) | 一种基于改进pso-bp神经网络齿轮箱故障诊断方法及装置 | |
CN109918752B (zh) | 基于迁移卷积神经网络的机械故障诊断方法、设备及介质 | |
CN110070060B (zh) | 一种轴承设备的故障诊断方法 | |
CN110132554B (zh) | 一种深度拉普拉斯自编码的旋转机械故障诊断方法 | |
CN111562108A (zh) | 一种基于cnn和fcmc的滚动轴承智能故障诊断方法 | |
CN112257530B (zh) | 基于盲信号分离和支持向量机的滚动轴承故障诊断方法 | |
Di et al. | Ensemble deep transfer learning driven by multisensor signals for the fault diagnosis of bevel-gear cross-operation conditions | |
Zhang et al. | A novel data-driven method based on sample reliability assessment and improved CNN for machinery fault diagnosis with non-ideal data | |
CN112257767B (zh) | 针对类不均衡数据的产品关键零部件状态分类方法 | |
CN113822139A (zh) | 一种基于改进1DCNN-BiLSTM的设备故障诊断方法 | |
CN111898686A (zh) | 基于门控循环单元网络的轴承故障识别方法 | |
CN114297918A (zh) | 基于全注意力深度网络和动态集成学习的航空发动机剩余寿命预测方法 | |
CN110837851A (zh) | 一种电静液作动器液压泵故障诊断方法 | |
CN115859077A (zh) | 一种变工况下多特征融合的电机小样本故障诊断方法 | |
CN116793682A (zh) | 基于iCORAL-MMD和对抗迁移学习的轴承故障诊断方法 | |
CN114091525A (zh) | 一种滚动轴承退化趋势预测方法 | |
CN116358871A (zh) | 基于图卷积网络的滚动轴承微弱信号复合故障诊断方法 | |
CN114757365A (zh) | 一种基于深度学习的高速铁路路基沉降预测及预警方法 | |
CN112763215B (zh) | 一种基于模块化联邦深度学习的多工况在线故障诊断方法 | |
CN117540285A (zh) | 基于能量熵和回归型支持向量机的轴承运行状态评估方法 | |
CN113269221A (zh) | 一种基于多尺度卷积神经网络的滚动轴承故障诊断方法 | |
CN112163630A (zh) | 一种基于不平衡学习的复合故障诊断方法及装置 | |
CN116625686A (zh) | 一种航空发动机轴承故障在线诊断方法 | |
CN116644348A (zh) | 基于传递式对抗迁移的跨机械部件故障诊断方法及装置 | |
Zhang et al. | Fault Diagnosis with Bidirectional Guided Convolutional Neural Networks Under Noisy Labels |
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 |