本发明涉及振动故障诊断，特别是一种用于齿轮箱电信号残差信号的提取方法。
背景技术：
齿轮箱具有重量轻、体积小、传动比大、承载能力强、传动效率高等诸多优点，被广泛应用于风力发电、直升机、船舶、冶金、石化、矿山、起重运输等行业的复杂机电装备系统。齿轮箱某一个部件出现故障，可能导致机电系统的停机，产生昂贵的修理和维修费用，甚至造成灾难性的事故，造成巨大的经济损失和恶劣的社会影响。齿轮箱早期故障检测与识别对于优化维修周期、减小经济损失有着十分重要的意义。振动信号具有非平稳性和调频特点，用于振动故障特征提取的占统治地位的关键技术主要是通过傅里叶变换进行频域特征提取。虽然可以获得满意的特征提取结果，但是以傅里叶变换为基础的频域特征提取方法也有一定的局限性，导致频谱估计发生非常明显的错误，造成微弱故障信号不能准确的检测和诊断。动态时间规整方法是一种有效的匹配校准定位时域信号通过匹配比较待测波形与参考波形时间特征量相似度的大小获得诊断结果的方法，具有概念简单、算法鲁棒的优点，是一种非线性的规划算法，它将时间和距离测度结合起来，能对存在全局或局部扩展、压缩或变形的模式进行匹配，解决动态模式的相似度量和分类问题，并在参考信号和测试信号之间找到一条优化的时间校准匹配路径。利用不同的动态信号处理方法提取特征信息以诊断齿轮箱的故障，得到了国内外学者的关注，已开展了大量的研究工作，并取得了一定的研究成果。
技术实现要素：
本发明的目的是为了解决上述问题，设计了一种用于齿轮箱电信号残差信号的提取方法。
实现上述目的本发明的技术方案为，一种用于齿轮箱电信号残差信号的提取方法，包括以下步骤：建立齿轮箱传动系统的振动信号和驱动电机的电流信号的耦合模型，计算得出故障边带特征的电信号残差向量信号。
所述建立齿轮箱传动系统的振动信号和驱动电机的电流信号的耦合模型步骤，其具体为：
(1)建立齿轮箱传动系统的振动信号的非线性集总参数模型。
(2)建立驱动电机的电流信号的等效二相电机。
(3)建立齿轮箱齿轮传动系统的振动信号和驱动电机的电流信号的耦合模型。
(4)使用时间步长积分法和单侧正常接触算法对耦合模型进行迭代求解，得到耦合电信号。
所述计算得出故障边带特征的电信号残差向量信号步骤，其具体为：
(1)耦合电信号带通滤波。
(2)耦合电信号初次相位补偿。
(3)savitzky-golay时间序列法对耦合信号平滑滤波。
(4)逻辑加权衍生动态时间规整算法对耦合电信号进行处理。
(5)得出故障边带特征的电信号残差向量信号。
所述建立齿轮箱的驱动电机电流信号和传动系统振动信号的耦合模型的步骤(1)中，齿轮箱传动系统的振动信号的非线性集总参数模型，该模型通过将齿轮简化为随时间变化的非线性弹簧连接的刚性气缸，转轴模拟为二节点单元，轴承模拟为附加集总刚度单元，从而建立振动信号的非线性集总参数模型。
所述建立齿轮箱的驱动电机电流信号和传动系统振动信号的耦合模型的步骤(2)中，驱动电机的电流信号的等效二相电机是采用克朗变换将静子电流信号系转换为动子d-q模型建立。
所述建立齿轮箱的驱动电机电流信号和传动系统振动信号的耦合模型的步骤(3)中，齿轮箱齿轮传动系统的振动信号和驱动电机的电流信号的耦合模型是非线性一阶差分动态耦合关系。
所述计算得出故障边带特征的电信号残差向量信号步骤(4)中，逻辑加权衍生动态时间规整算法是将逻辑权重函数和坡度加权函数应用于动态时间规整算法之中。
所述计算得出故障边带特征的电信号残差向量信号步骤(4)中，逻辑加权衍生动态时间规整算法计算公式：
其中：
一种用于齿轮箱电信号故障特征提取的装置，包括驱动的电机电信号测量模块、信号耦合模块和电信号残差向量信号计算模块。
所述驱动电机的电流信号模块具有电流信号输出端口，电流信号输出端口将电流信号传递到信号耦合模块。
所述信号耦合模块具有信号耦合功能，输出耦合信号到电信号残差向量信号计算模块。
所述电信号残差向量信号计算模块内置有齿轮箱电信号故障特征提取算法程序，将耦合模块的耦合信号计算处理得到电信号残差向量信号。
利用本发明的技术方案制作的用于齿轮箱电信号残差信号的提取方法，通过该方法得到齿轮箱传动系统的振动信号与驱动电机的电流信号的动态非线性耦合作用规律与机理，从而得到其耦合信号，通过理论分析和仿真实验选择统计特征参数进行故障检测和定位，由于不同的电信号残差信号向量信号带有不同的故障信息，因此，可以准确定位故障来源，该方法计算简单，应用范围广泛，能够精准预测故障趋势，满足齿轮箱齿轮、架、轴承等不同故障诊断的需要。
附图说明
图1是本发明所述用于齿轮箱电信号残差信号的提取方法的流程图；
图2是本发明所述用于齿轮箱电信号故障特征提取的装置的工作流程图；
图中，1、驱动的电机电信号测量模块；2、耦合模块；3、电信号残差向量信号计算模块。
具体实施方式
下面结合附图对本发明进行具体描述，如图1-2所示，一种用于齿轮箱电信号残差信号的提取方法，一种用于齿轮箱电信号残差信号的提取方法，包括以下步骤：建立齿轮箱传动系统的振动信号和驱动电机的电流信号的耦合模型，计算得出故障边带特征的电信号残差向量信号。
所述建立齿轮箱传动系统的振动信号和驱动电机的电流信号的耦合模型步骤，其具体为：
(1)建立齿轮箱传动系统的振动信号的非线性集总参数模型，该模型通过将齿轮简化为随时间变化的非线性弹簧连接的刚性气缸，转轴模拟为二节点单元，轴承模拟为附加集总刚度单元，从而建立振动信号的非线性集总参数模型。
(2)建立驱动电机的电流信号的等效二相电机，采用克朗变换将静子电流信号系转换为动子d-q模型建立。
(3)建立齿轮箱齿轮传动系统的振动信号和驱动电机的电流信号的耦合模型，是非线性一阶差分动态耦合关系。
(4)使用时间步长积分法和单侧正常接触算法对耦合模型进行迭代求解，得到耦合电信号。
所述计算得出故障边带特征的电信号残差向量信号步骤，其具体为：
(1)耦合电信号带通滤波。
(2)耦合电信号初次相位补偿。
(3)savitzky-golay时间序列法对耦合信号平滑滤波。
(4)逻辑加权衍生动态时间规整算法对耦合电信号进行处理，将逻辑权重函数和坡度加权函数应用于动态时间规整算法之中，最终逻辑加权衍生动态时间规整算法计算公式：
其中：
(5)得出故障边带特征的电信号残差向量信号。
一种用于齿轮箱电信号故障特征提取的装置，包括驱动的电机电信号测量模块、信号耦合模块和电信号残差向量信号计算模块，驱动电机的电流信号模块具有电流信号输出端口，电流信号输出端口将电流信号传递到信号耦合模块，信号耦合模块具有信号耦合功能，输出耦合信号到电信号残差向量信号计算模块，电信号残差向量信号计算模块内置有齿轮箱电信号故障特征提取算法程序，将耦合模块的耦合信号计算处理得到电信号残差向量信号。
在本实施方案中，在求解耦合模型时，通过使用时间步长积分法和单侧正常接触算法进行迭代求解，得到耦合电信号初值，利用savitzky-golay时间序列平滑滤波方法对相位补偿后的耦合电信号进行平滑滤波处理，增强算法的平滑性和鲁棒性，在时域内基于多项式，通过移动窗口利用最小二乘法进行最佳拟合，对处于l×m滑动窗口(l，m为窗口大小，且均为奇数)内各点的值d，进行p次二元多项式拟合，多项式对应的系数根据最小二乘法准则使拟合误差最小来确定，由此得出l×m滑动窗口内中心点(i，j)的最佳拟合值p(i，j)，即为d(i，j)去噪处理后的新值，滑动数据窗口依次沿着三维空间每一点滑动，从而实现了平滑处理。
在本实施方案中，逻辑加权衍生动态时间规整优化方法，将逻辑权重函数和坡度加权函数的方法，建立基于相位偏差的修正逻辑权重函数，最优分配时间序列各点数据的权重，使用衍生技术时间序列数据处理为带有波形特征的数据，建立规整矩阵和优化规整路径，对耦合电信号进行逻辑加权衍生动态时间规整处理，修正逻辑权重函数：
其中：i＝1，...，m；m为时间序列长度；
mc为时间序列中间点；
wmax为权重系数最大值；
g为惩罚因子；
逻辑加权衍生动态时间规整处理方法：
通过解方程得到具有故障边带特征的电信号残差向量信号，通过理论分析和仿真实验选择统计特征参数进行故障检测和定位，得到故障的残差信号向量故障特征与定位算法，预测故障趋势，满足齿轮箱齿轮、架、轴承等不同故障诊断的需要。
上述技术方案仅体现了本发明技术方案的优选技术方案，本技术领域的技术人员对其中某些部分所可能做出的一些变动均体现了本发明的原理，属于本发明的保护范围之内。