1.一种基于稀疏低秩协同正则框架的齿轮箱局部故障诊断方法,包括如下步骤:S100:采集齿轮箱的振动信号并构建分块算子;
S200:根据所述分块算子建立稀疏低秩协同优化框架;
所述稀疏低秩协同优化框架表示为:
其中,表示振动信号矩阵, 表示故障特征信号矩阵, 表示从振动信号 中分离出来的故障特征信号, 表示从振动信号 中分离出来的谐波干扰信号, 表示求解得到的故障特征信号矩阵, 表示分块算子, 表示分块逆算子,为分块算子 的逆过程,为谐波分量的稀疏表示字典,为谐波分量的稀疏表示系数,表示求解得到的稀疏表示系数, 为正则化参数, 表示矩阵的 范数, 表示稀疏正则范数且 ,为矩阵 的低秩正则约束,所述稀疏低秩协同优化框架简化为:
其中,稀疏范数 中令 ,得到 , 为加权核范数,为权系数向量;
S300:求解所述稀疏低秩协同优化框架,获得振动信号中所含的故障特征信号;
所述步骤S300包括如下步骤:
S301:将所述稀疏低秩协同优化框架进行简化,将对简化后的稀疏低秩协同优化框架的求解转化为对增广拉格朗日方程的求解:其中,为代理中间变量,且 , 为罚参数,为拉格朗日乘子,是由已知的振动信号经分块算子处理得到的振动信号矩阵; 为谐波分量的稀疏表示字典, 均为拉格朗日方程中的待求解变量;
S302:基于所述增广拉格朗日方程,获得所述待求解变量 的闭式解;
S303:令迭代数 ,重复步骤S302直到满足停止准则,停止迭代,并获得 ;
S304:根据所获得的 ,由分块逆算子对 进行分块逆操作可获得从振动信号中分离出来的故障特征信号 ,即 ;
S400:通过对所述故障特征信号进行包络分析,诊断齿轮箱局部故障并辨识故障类型。
2.根据权利要求1所述的方法,其中,所述齿轮箱的振动信号表示为:其中, 为一维向量形式的齿轮箱局部故障特征信号, 为齿轮箱谐波干扰信号, 为噪声干扰信号, 为采集到的振动信号的点数, 表示实数域。
3.根据权利要求1所述的方法,其中,所述分块算子表示为:且
其中, 为信号块的长度, 为信号块的个数, 表示对采集的振动信号 经过分块算子 处理后得到的二维振动信号矩阵 , 中每一列表示一个信号块,表示实数域;
相应的,
表示齿轮箱局部故障特征信号经过分块算子处理后得到的二维故障特征信号矩阵;
表示齿轮箱谐波干扰信号经过分块算子处理后得到的二维谐波干扰信号矩阵;
表示噪声干扰信号经过分块算子处理后得到的二维噪声干扰信号矩阵。
4.根据权利要求3所述的方法,其中,信号块的长度 和信号块的个数 分别表示为:其中, 为采样频率, 为潜在发生局部故障部件的故障特征频率, 为相邻两个信号块之间的冗余度,其取值范围为 , 表示向上取整符号。
5.根据权利要求1所述的方法,其中,所述待求解变量 的闭式解分别为:,且
,且
,对 进行奇异值分解可得
其中, 是符号函数,提取变量的正负号; 是相对值比较函数,取较大值;
表示向量或矩阵逐点运算操作,是函数变量, 为迭代数, 表示矩阵 的奇异值分解, 为奇异矩阵, 为右奇异矩阵, 为相应奇异值序列构成的对角阵, 表示转置运算; 表示逆运算, 中上标 表示第 次迭代。
6.根据权利要求1所述的方法,其中,所述停止准则包括如下任一:(1)当前迭代数 等于最大迭代数 ;
(2)目标函数 的函数值在相
邻两次迭代中的相对误差小于给定误差 。
7.根据权利要求1所述的方法,其中,所述故障类型包括如下任一:齿轮箱中轴承外圈局部故障、轴承内圈局部故障、轴承滚动体局部故障、齿轮齿面局部故障。