CN105823632B - 行星齿轮箱故障信息获取方法 - Google Patents
行星齿轮箱故障信息获取方法 Download PDFInfo
- Publication number
- CN105823632B CN105823632B CN201610019138.5A CN201610019138A CN105823632B CN 105823632 B CN105823632 B CN 105823632B CN 201610019138 A CN201610019138 A CN 201610019138A CN 105823632 B CN105823632 B CN 105823632B
- Authority
- CN
- China
- Prior art keywords
- tooth root
- tooth
- variable period
- signal
- sta
- 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.)
- Active
Links
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
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种行星齿轮箱故障信息获取方法，包括组成行星齿轮箱的内齿圈、行星轮和太阳轮，其中内齿圈固定不动，太阳轮绕自身中心轴线转动，行星轮除自转外还绕太阳轮公转；所述内齿圈的若干连续轮齿作为检测齿，检测齿的齿根部位粘贴有光纤光栅阵列，宽带光源发出的光通过光纤a、耦合器和光纤c传输到光纤光栅阵列，反射光通过光纤c、耦合器和光纤b进入波长解调部分，解调完成之后的电压信号通过数据采集部分进入计算机进行数据处理，计算机用于完成齿根应变序列构建和行星齿轮箱故障特征提取，本发明能够克服现有技术由于振动信号强烈的非线性和非平稳性特点所造成的行星齿轮箱故障诊断的难题。
Description
技术领域
本发明涉及机械传动装置的动态测试与故障诊断技术，具体涉及行星齿轮箱故障信息获取方法。
背景技术
行星齿轮箱由于体积小、重量轻、传动效率高、承载能力强等优点被广泛用于航空、风力发电、起重运输、能源等领域的机械传动***中，然而受制造、安装误差和低速重载的恶劣工作环境影响，行星齿轮箱发生故障的概率很高，由于行星齿轮箱位于机械***的关键环节，当其出现故障后，若不及时采取有效的补救措施，轻则会影响***性能，重则会使***瘫痪，造成严重事故，因此，对行星齿轮箱进行故障诊断具有重要意义。
目前国内外研究人员绝大部分是通过对振动信号的分析处理来判断行星齿轮箱的健康状况，振动传感器安装在与内齿圈相连的箱体上，将振动信号作为故障载体信号的主要优势在于其测量简单、实时性强且包含有丰富的故障信息，然而行星齿轮箱的振动信号呈现出比一般机械设备更复杂的特点：多个啮合振动之间相互耦合且振动信号被传输路径调制，导致传感器测到的振动信号中包含有多种频率成分和丰富的边频带信息，表现出强烈的非线性特点；变速、变载工况导致振动信号中的齿轮啮合频率和故障特征频率发生变化，表现出强烈的非平稳性特点，由于以上原因，尽管国内外学者在动力学建模和振动信号处理等方面进行了大量的工作，但现在行星齿轮箱的故障诊断仍是一个技术难题。
实际上，行星齿轮箱中最关键的运动是齿轮之间的啮合运动，啮合力是能够反映齿轮啮合状态和特性的最直接参数，然而啮合力由轮齿相互接触、挤压产生，不可能在轮齿的接触区域布置传感器以测量该力，因此研究者选择振动信号来间接反映啮合力的变化情况，但在行星齿轮箱中啮合力信息与振动信息的映射关系非常复杂，导致出现了行星齿轮箱故障诊断的困境，考虑到轮齿的齿根区域和啮合区域距离非常近，啮合力的变化会直接导致齿根应变的变化，所以与振动信号相比，齿根应变信号更能够反映齿轮的啮合状态，是一种更加优质的故障特征载体信号，然而，到目前为止，基于齿根应变信号的行星齿轮箱故障诊断方法仍未见报道。
发明内容
本发明的目的是提供一种行星齿轮箱故障信息获取方法，该方法可解决背景技术中由于振动信号强烈的非线性和非平稳性特点所造成的行星齿轮箱故障诊断的难题。
为了达到以上目的，本发明采取如下技术方案予以实现：行星齿轮箱故障信息获取***，包括组成行星齿轮箱的内齿圈、行星轮和太阳轮，其中内齿圈固定不动，太阳轮绕自身中心轴线转动，行星轮除自转外还绕太阳轮公转；所述内齿圈的若干连续轮齿作为检测齿，检测齿的齿根部位粘贴有光纤光栅阵列，宽带光源发出的光通过光纤a、耦合器和光纤c传输到光纤光栅阵列，在光纤光栅阵列发生光反射，反射光通过光纤c、耦合器和光纤b进入波长解调部分，解调完成之后的电压信号通过数据采集部分进入计算机进行数据处理，计算机用于完成齿根应变序列构建和行星齿轮箱故障特征提取。
所述光纤光栅阵列为在一根光纤上刻写多段光纤光栅，光纤光栅的数量与检测齿的数量相同。
所述光纤光栅的长度小于检测齿的齿宽，光纤光栅的轴向与内齿圈的轴向平行。
所述波长解调部分用于将光纤光栅阵列的反射波长信号转换为包含有应变信息的电压信号，包括光电转换模块、信号调理模块、电源模块、驱动电路模块、控制电路模块以及光学元器件。
所述耦合器采用3dB耦合器。
本发明还提出一种行星齿轮箱故障信息获取方法，包括以下步骤：
步骤一、选定内齿圈的连续n个轮齿作为检测齿，并在检测齿的齿根部位粘贴光纤光栅阵列，获取检测齿的齿根应变信号；
步骤二、对步骤一中获取的齿根应变信号进行序列构建；
步骤三、对步骤二构建出的齿根应变序列进行行星齿轮箱故障特征提取：
在每一组齿根应变序列中，将相邻两个轮齿的齿根应变信号相加，形成N+1组时域上连续的齿根应变信号，设反映行星轮j和检测齿直接啮合的齿根应变连续时域信号为fpj(t)，反映太阳轮和检测齿间接啮合的齿根应变连续时域信号为fs(t)，行星齿轮箱的啮合频率为fm，行星轮齿数、太阳轮齿数、内齿圈齿数分别为Zp、Zs、Zr，行星轮数量为N；
3.1提取行星齿轮箱故障的频域特征
若在fpj(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明行星轮j出现故障；
若在fs(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明太阳轮出现故障；
3.2提取行星齿轮箱故障的时域特征
取fpj(t)区间(k为自然数)内的信号fpj(k,t)，设gpj(k,t)＝fpj(0,t)-fpj(k,t)；若行星轮j发生故障，则gpj(k,t)在某一区间(t1,t2)内的能量水平(有效值)便会偏高，(t1,t2)区间在gpj(k,t)中的位置反映故障在行星轮j上出现的位置，(t1,t2)区间的长度反映行星轮j故障的范围，(t1,t2)区间内gpj(k,t)的形式反映行星轮j故障的类型，(t1,t2)区间内gpj(k,t)有效值的大小反映行星轮j故障的剧烈程度；
取fs(t)区间内的信号fs(k,t)，设gs(k,t)＝fs(0,t)-fs(k,t)；若太阳轮发生故障，则gs(k,t)在某一区间(t1,t2)内的能量水平(有效值)便会偏高，(t1,t2)区间在gs(k,t)中的位置反映故障在太阳轮上出现的位置，(t1,t2)区间的长度反映太阳轮故障的范围，(t1,t2)区间内gs(k,t)的形式反映太阳轮故障的类型，(t1,t2)区间内gs(k,t)有效值的大小反映太阳轮故障的剧烈程度。
所述步骤一具体包括以下步骤：
1.1、确定内齿圈上所需连续检测齿的数量n：设在行星齿轮箱中，行星轮齿数、太阳轮齿数、内齿圈齿数分别为Zp、Zs、Zr；行星轮数量为N；rem(w,v)表示w与v相除的余数；i为正整数；令其中a、c为整数，b、d为真分数；H＝bZp；K＝dZs；(若rem(b·2π·i,2π)＝0，则anglep(i)＝2π；若rem(b·2π·i,2π)≠0，则anglep(i)＝rem(b·2π·i,2π))；(若rem(d·2π·i,2π)＝0，则angles(i)＝2π；若rem(d·2π·i,2π)≠0，则angles(i)＝rem(d·2π·i,2π))；
分以下四种情况讨论：
(1)当b＝0且d≠0时，则n为Zp和K中的最大值；
(2)当b≠0且d＝0时，则n为Zs和H中的最大值；
(3)当b＝0且d＝0时，则n为Zp和Zs中的最大值；
(4)当b≠0且d≠0时，设nmax(i)为np(i)和ns(i)中的最大值，当i＝is时，nmax(i)最小，则n＝nmax(is)；
1.2、安装光纤光栅阵列
选择在一根光纤上刻写有n段光栅的光纤光栅阵列，在每个检测齿的齿根部位粘贴一段光纤光栅。
所述步骤二具体包括以下步骤：
每次行星轮与所有检测齿完成一次啮合时，光纤光栅阵列检测到一组齿根应变信号，称为一个齿根应变周期；一个齿根应变周期的持续时间为T1，相邻两个齿根应变周期的时间间隔为T2，每个齿根应变周期的起始端标记为STA，距离终止端(向前)时间间隔为te处标记为END，te为相邻两个齿根应变信号重叠部分的时间长度，在以下的构建过程中，需要对齿根应变周期进行截取，截取信号的起始端标记为STAJ(与STA重合)，距离终止端(向前)时间间隔为te处标记为ENDJ，齿根应变周期编号为p，选定任一组齿根应变周期为构建序列的起点，从该起点开始，p值依次为1，2，3……；
2.1、构建反映行星轮与检测齿直接啮合的齿根应变序列
设所构建的序列针对编号为j(j＝1,2,3......N)的行星轮，若rem(p,N)≠0，则j＝rem(p,N)，若rem(p,N)＝0，则j＝N，按照下述方法，共构建N组反映不同行星轮和检测齿直接啮合的齿根应变序列；
(1)对于齿根应变信号获取方法第一步中的第(1)种情况
若n＝Zp，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝K，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(2)对于齿根应变信号获取方法第一步中的第(2)种情况
若n＝H，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝Zs，首先截取每个齿根应变周期中从STA端开始H个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(3)对于齿根应变信号获取方法第一步中的第(3)种情况
若n＝Zp，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝Zs，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(4)对于齿根应变信号获取方法第一步中的第(4)种情况
若n＝np(is)，将时间间隔为isNT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝ns(is)，首先截取每个齿根应变周期中从STA端开始np(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
2.2、构建反映太阳轮与检测齿间接啮合的齿根应变序列
(1)对于齿根应变信号获取方法第一步中的第(1)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始K个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝K，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(2)对于齿根应变信号获取方法第一步中的第(2)种情况
若n＝H，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝Zs，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(3)对于齿根应变信号获取方法第一步中的第(3)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝Zs，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(4)对于齿根应变信号获取方法第一步中的第(4)种情况
若n＝np(is)，首先截取一个齿根应变周期中从STA端开始ns(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝ns(is)，将时间间隔为isNT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
与现有技术相比较，本发明至少具有以下有益效果：
(1)本发明使用齿根应变信号作为行星齿轮箱故障特征载体信号，与传统的振动信号相比，该信号受故障的影响更直接，信噪比更高，更能反映行星齿轮箱的健康状况；
(2)本发明利用光纤光栅阵列获取内齿圈的多点齿根应变信息，传感器安装方便，测试***结构简单，测量精度高；
(3)本发明根据行星齿轮箱的参数优化应变测点的数量，可有效降低测试***成本，提高本方法的工程应用价值；
(4)本发明对内齿圈齿根应变信号进行时间序列构建，并从时域和频域两个方面提取行星齿轮箱的故障特征，与现有技术相比，能够更加直观、准确地获得行星齿轮箱的故障位置、故障范围、故障剧烈程度等信息。
附图说明
图1为本发明一种光纤获取行星齿轮箱故障信息的方法框图；
图2为本发明在确定检测齿数目时用到的曲线图；
图3为本发明光纤光栅阵列在齿根上的粘贴方式示意图；
图4为本发明光纤光栅阵列测到的齿根应变信号示意图；
图5为本发明将齿根应变信号时间序列转化为连续时域信号的示意图；
图6为本发明行星齿轮箱故障的时域特征提取方法示意图。
附图中：1、宽带光源，2、耦合器，3、内齿圈，4、行星轮，5、太阳轮，6、光纤光栅阵列，7、波长解调部分，8、数据采集部分，9、计算机，10、光纤光栅。
具体实施方式
以下结合附图和具体实施方式对本发明作进一步的详细说明。
如图1所示，本发明包括组成行星齿轮箱的内齿圈3、行星轮4和太阳轮5，在行星齿轮箱中，内齿圈3固定不动，太阳轮5绕自身中心轴线转动，行星轮4除自转外还绕太阳轮公转，内齿圈3的若干连续轮齿作为检测齿，光纤光栅阵列6粘贴在内齿圈3的齿根部位，宽带光源1发出的光通过光纤a、3dB耦合器2和光纤c传输到光纤光栅阵列6，在光纤光栅阵列6发生光反射，反射光通过光纤c、3dB耦合器2和光纤b进入波长解调部分7，解调完成之后的电压信号通过数据采集部分8进入计算机9内进行数据处理，完成齿根应变序列构建和行星齿轮箱故障特征提取。
光纤光栅阵列部分：为了简化测试***结构，提高本发明的工程应用价值，采用光纤光栅阵列6作为应变传感器，即在一根光纤上刻写多段光纤光栅10，各光纤光栅10的结构参数不同，当宽带光入射到光纤光栅阵列6中时，会在多个波长上形成反射光，反射光的波长值与光纤光栅的折射率和周期有关，当某个光纤光栅10发生轴向变形时，其对应的反射波长值会在原始位置附近发生变化，而波长的偏移量与应变量之间为正比关系，所以通过检测光纤光栅阵列6中各段光纤光栅10反射波长值的变化，即可得到测点处的应变值，实现分布式应变测量。
波长解调部分：主要包括光电转换模块、信号调理模块、电源模块、驱动电路模块、控制电路模块以及光学元器件。该部分的目的是将光纤光栅阵列6的反射波长信号转换为包含有应变信息的电压信号。
数据采集部分：通过数据采集部分8对波长解调部分7传输过来的电压信号进行A/D转换、采集并向计算机9传输数据。
数据处理部分：数据处理工作在计算机9中完成，通过软件对数据采集部分8得到的信号进行数据处理、数据显示等操作，并完成***的参数设计、控制指令的下达、***工作状态的监测，该部分的主要工作是进行齿根应变信号时间序列构建和行星齿轮箱故障特征提取。
本发明的方法包括以下步骤：
步骤一、利用光纤光栅阵列获取齿根应变信号
选定内齿圈3的连续n个轮齿作为检测齿，在这n个检测齿的齿根部位粘贴光纤光栅10。
1.1、确定检测齿的数量n。设在行星齿轮箱中，行星轮齿数、太阳轮齿数、内齿圈齿数分别为Zp、Zs、Zr；行星轮数量为N；rem(w,v)表示w与v相除的余数；i为正整数。令其中a、c为整数，b、d为真分数；H＝bZp；K＝dZs。分以下四种情况讨论：
(1)若b＝0且d≠0，则n为Zp和K中的最大值；
(2)若b≠0且d＝0，则n为Zs和H中的最大值；
(3)若b＝0且d＝0，则n为Zp和Zs中的最大值；
(4)若b≠0且d≠0，首先计算行星轮公转i周时，行星轮、太阳轮上的某一点在各自的自转方向上相对于初始位置的角位移anglep(i)和angles(i)。其中，若rem(b·2π·i,2π)＝0，则anglep(i)＝2π，否则anglep(i)＝rem(b·2π·i,2π)；若rem(d·2π·i,2π)＝0，则angles(i)＝2π，否则angles(i)＝rem(d·2π·i,2π)。然后计算以i周为检测周期时，检测行星轮所有齿和检测太阳轮所有齿所需要的检测齿数量np(i)和ns(i)。其中，np(i)＝anglep(i)·Zp/(2π)；ns(i)＝angles(i)·Zs/(2π)。nmax(i)为np(i)和ns(i)中的最大值，设当i＝is时，nmax(i)最小，则n＝nmax(is)。
对于第(4)种情况，举例说明如下：
设在某一行星齿轮箱中，Zs＝19，Zp＝31，Zr＝81，N＝4。按照上述方法做出np(i)和ns(i)的曲线如图2所示。在图2中，可以看到，当i＝5时，nmax(i)＝6，为最小值，所以，检测齿数量n＝6。
1.2、安装光纤光栅阵列。如图3所示，在每个检测齿的齿根部位粘贴一段光纤光栅。光纤光栅的长度Wf小于齿宽Wg，光纤光栅的轴向11与内齿圈的轴向12平行。
步骤二、构建齿根应变序列
如图4所示，每次行星轮与所有检测齿完成一次啮合时，光纤光栅阵列会检测到一组齿根应变信号，该组信号由n个检测齿的齿根应变信号组成，包含有行星轮与检测齿直接啮合的信息和太阳轮与检测齿间接啮合的信息，称为一个齿根应变周期。一个齿根应变周期的持续时间为T1，相邻两个齿根应变周期的时间间隔为T2，每个齿根应变周期的起始端标记为STA，距离终止端(向前)时间间隔为te处标记为END，te为相邻两个齿根应变信号重叠部分的时间长度。在以下的构建过程中，需要对齿根应变周期进行截取，截取信号的起始端标记为STAJ(与STA重合)，距离终止端(向前)时间间隔为te处标记为ENDJ。齿根应变周期编号为p，选定任一组齿根应变周期为构建序列的起点，从该起点开始，p值依次为1，2，3……；对图4所示信号进行序列构建的目的是得到行星轮上连续的所有轮齿在与检测齿发生直接啮合时或太阳轮上连续的所有轮齿在与检测齿发生间接啮合时，检测齿的齿根应变信号。具体将得到以下两类信号：
2.1、反映行星轮与检测齿直接啮合的齿根应变序列
设所构建的序列针对编号为j(j＝1,2,3......N)的行星轮，若rem(p,N)≠0，则j＝rem(p,N)，若rem(p,N)＝0，则j＝N。按照下述方法，共构建N组反映不同行星轮和检测齿直接啮合的齿根应变序列。
(1)对于步骤1.2中的第(1)种情况
若n＝Zp，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
若n＝K，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
(2)对于步骤1.2中的第(2)种情况
若n＝H，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
若n＝Zs，首先截取每个齿根应变周期中从STA端开始H个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
(3)对于步骤1.2中的第(3)种情况
若n＝Zp，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
若n＝Zs，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
(4)对于步骤1.2中的第(4)种情况
若n＝np(is)，将时间间隔为isNT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
若n＝ns(is)，首先截取每个齿根应变周期中从STA端开始np(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
2.2、反映太阳轮与检测齿间接啮合的齿根应变序列
(1)对于步骤1.2中的第(1)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始K个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
若n＝K，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
(2)对于步骤1.2中的第(2)种情况
若n＝H，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
若n＝Zs，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
(3)对于步骤1.2中的第(3)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
若n＝Zs，将时间间隔为NT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
(4)对于步骤1.2中的第(4)种情况
若n＝np(is)，首先截取一个齿根应变周期中从STA端开始ns(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的截取信号进行构建，即将后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端。
若n＝ns(is)，将时间间隔为isNT2的齿根应变周期进行构建，即将后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
步骤三、基于齿根应变序列的行星齿轮箱故障特征的提取
按照以上齿根应变序列的构建方法，构建N+1组齿根应变序列。在每一组序列中，将相邻两个齿根应变信号相加，形成N+1组时域上连续的齿根应变信号，如图5所示。设反映行星轮j和检测齿直接啮合的齿根应变连续时域信号为fpj(t)，反映太阳轮和检测齿间接啮合的齿根应变连续时域信号为fs(t)，行星齿轮箱的啮合频率为fm。
3.1、行星齿轮箱故障的频域特征提取
一个行星轮的所有轮齿与内齿圈完成一次啮合所需要的时间为即行星轮的自转周期为当行星轮出现故障时，会出现啮合频率被故障频率(即行星轮自转频率)调制的现象，所以若在fpj(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明行星轮j出现故障。
太阳轮的所有轮齿与内齿圈完成一次间接啮合所需要的时间为即太阳轮相对于行星架的旋转周期为当太阳轮出现故障时，会出现啮合频率被故障频率(即太阳轮相对于行星架的旋转频率)调制的现象，所以若在fs(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明太阳轮出现故障。
3.2、行星齿轮箱故障的时域特征提取
如图6所示，取fpj(t)区间(k为自然数)内的信号fpj(k,t)，fpj(k,t)表示在行星轮j所有轮齿与检测齿完成一次啮合的过程中，检测齿的齿根应变连续时域信号，区间长度为行星轮j的一个自转周期，不同的k值代表行星轮j的不同自转周期。设gpj(k,t)＝fpj(0,t)-fpj(k,t)，gpj(k,t)表示在行星轮j不同自转周期下，检测齿齿根应变连续时域信号的差异。当gpj(k,t)＝0时，说明行星轮j在不同自转周期下，检测齿齿根应变连续时域信号没有差异，则行星轮j没有出现故障，若行星轮j发生故障，则gpj(k,t)在某一区间(t1,t2)内的能量水平(有效值)便会偏高。(t1,t2)区间在gpj(k,t)中的位置反映故障在行星轮j上出现的位置，(t1,t2)区间的长度反映行星轮j故障的范围，(t1,t2)区间内gpj(k,t)的形式反映行星轮j故障的类型，(t1,t2)区间内gpj(k,t)有效值的大小反映行星轮j故障的剧烈程度。
如图6所示，取fs(t)区间内的信号fs(k,t)，fs(k,t)表示在太阳轮所有轮齿与检测齿完成一次间接啮合的过程中，检测齿的齿根应变连续时域信号，区间长度为太阳轮相对于行星架的一个旋转周期，不同的k值代表太阳轮的不同相对旋转周期。设gs(k,t)＝fs(0,t)-fs(k,t)，gs(k,t)表示在太阳轮不同相对旋转周期下，检测齿齿根应变连续时域信号的差异。当gs(k,t)＝0时，说明太阳轮在不同相对旋转周期下，检测齿齿根应变连续时域信号没有差异，则太阳轮没有出现故障，若太阳轮发生故障，则gs(k,t)在某一区间(t1,t2)内的能量水平(有效值)便会偏高。(t1,t2)区间在gs(k,t)中的位置反映故障在太阳轮上出现的位置，(t1,t2)区间的长度反映太阳轮故障的范围，(t1,t2)区间内gs(k,t)的形式反映太阳轮故障的类型，(t1,t2)区间内gs(k,t)有效值的大小反映太阳轮故障的剧烈程度。
Claims (3)
1.行星齿轮箱故障信息获取方法，其特征在于，包括以下步骤：
步骤一、选定内齿圈的连续n个轮齿作为检测齿，并在检测齿的齿根部位粘贴光纤光栅阵列，获取检测齿的齿根应变信号；
步骤二、对步骤一中获取的齿根应变信号进行序列构建；
步骤三、对步骤二构建出的齿根应变序列进行行星齿轮箱故障特征提取：
在每一组齿根应变序列中，将相邻两个轮齿的齿根应变信号相加，形成N+1组时域上连续的齿根应变信号，设反映行星轮j和检测齿直接啮合的齿根应变连续时域信号为fpj(t)，反映太阳轮和检测齿间接啮合的齿根应变连续时域信号为fs(t)，行星齿轮箱的啮合频率为fm，行星轮齿数、太阳轮齿数、内齿圈齿数分别为Zp、Zs、Zr，行星轮数量为N；
3.1提取行星齿轮箱故障的频域特征
若在fpj(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明行星轮j出现故障；
若在fs(t)的功率谱中，啮合频率及其谐波附近出现间隔为的边频带，则说明太阳轮出现故障；
3.2提取行星齿轮箱故障的时域特征
取fpj(t)区间(k为自然数)内的信号fpj(k,t)，设gpj(k,t)＝fpj(0,t)-fpj(k,t)；若行星轮j发生故障，则gpj(k,t)在某一区间(t1,t2)内的能量水平有效值便会偏高，(t1,t2)区间在gpj(k,t)中的位置反映故障在行星轮j上出现的位置，(t1,t2)区间的长度反映行星轮j故障的范围，(t1,t2)区间内gpj(k,t)的形式反映行星轮j故障的类型，(t1,t2)区间内gpj(k,t)有效值的大小反映行星轮j故障的剧烈程度；
取fs(t)区间内的信号fs(k,t)，设gs(k,t)＝fs(0,t)-fs(k,t)；若太阳轮发生故障，则gs(k,t)在某一区间(t1,t2)内的能量水平(有效值)便会偏高，(t1,t2)区间在gs(k,t)中的位置反映故障在太阳轮上出现的位置，(t1,t2)区间的长度反映太阳轮故障的范围，(t1,t2)区间内gs(k,t)的形式反映太阳轮故障的类型，(t1,t2)区间内gs(k,t)有效值的大小反映太阳轮故障的剧烈程度。
2.根据权利要求1所述的行星齿轮箱故障信息获取方法，其特征在于，步骤一具体包括以下步骤：
1.1、确定内齿圈上所需连续检测齿的数量n：设在行星齿轮箱中，行星轮齿数、太阳轮齿数、内齿圈齿数分别为Zp、Zs、Zr；行星轮数量为N；rem(w,v)表示w与v相除的余数；i为正整数；令其中a、c为整数，b、d为真分数；H＝bZp；K＝dZs；(若rem(b·2π·i,2π)＝0，则anglep(i)＝2π；若rem(b·2π·i,2π)≠0，则anglep(i)＝rem(b·2π·i,2π))；(若rem(d·2π·i,2π)＝0，则angles(i)＝2π；若rem(d·2π·i,2π)≠0，则angles(i)＝rem(d·2π·i,2π))；
分以下四种情况讨论：
(1)当b＝0且d≠0时，则n为Zp和K中的最大值；
(2)当b≠0且d＝0时，则n为Zs和H中的最大值；
(3)当b＝0且d＝0时，则n为Zp和Zs中的最大值；
(4)当b≠0且d≠0时，设nmax(i)为np(i)和ns(i)中的最大值，当i＝is时，nmax(i)最小，则n＝nmax(is)；
1.2、安装光纤光栅阵列
选择在一根光纤上刻写有n段光栅的光纤光栅阵列，在每个检测齿的齿根部位粘贴一段光纤光栅。
3.根据权利要求2所述的行星齿轮箱故障信息获取方法，其特征在于，所述步骤二具体包括以下步骤：
每次行星轮与所有检测齿完成一次啮合时，光纤光栅阵列检测到一组齿根应变信号，称为一个齿根应变周期；一个齿根应变周期的持续时间为T1，相邻两个齿根应变周期的时间间隔为T2，每个齿根应变周期的起始端标记为STA，距离终止端向前的时间间隔为te处标记为END，te为相邻两个齿根应变信号重叠部分的时间长度，在以下的构建过程中，需要对齿根应变周期进行截取，截取信号的起始端标记为STAJ，STAJ与STA重合，距离终止端向前的时间间隔为te处标记为ENDJ，齿根应变周期编号为p，选定任一组齿根应变周期为构建序列的起点，从该起点开始，p值依次为1，2，3……；
2.1，构建反映行星轮与检测齿直接啮合的齿根应变序列
设所构建的序列针对编号为j(j＝1,2,3......N)的行星轮，若rem(p,N)≠0，则j＝rem(p,N)，若rem(p,N)＝0，则j＝N，按照下述方法，共构建N组反映不同行星轮和检测齿直接啮合的齿根应变序列；
(1)对于齿根应变信号获取方法第一步中的第(1)种情况
若n＝Zp，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝K，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(2)对于齿根应变信号获取方法第一步中的第(2)种情况
若n＝H，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝Zs，首先截取每个齿根应变周期中从STA端开始H个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(3)对于齿根应变信号获取方法第一步中的第(3)种情况
若n＝Zp，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝Zs，首先截取每个齿根应变周期中从STA端开始Zp个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
(4)对于齿根应变信号获取方法第一步中的第(4)种情况
若n＝np(is)，将时间间隔为isNT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
若n＝ns(is)，首先截取每个齿根应变周期中从STA端开始np(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
2.2，构建反映太阳轮与检测齿间接啮合的齿根应变序列
(1)对于齿根应变信号获取方法第一步中的第(1)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始K个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝K，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(2)对于齿根应变信号获取方法第一步中的第(2)种情况
若n＝H，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝Zs，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(3)对于齿根应变信号获取方法第一步中的第(3)种情况
若n＝Zp，首先截取每个齿根应变周期中从STA端开始Zs个检测齿的齿根应变信号，然后将时间间隔为NT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝Zs，将时间间隔为NT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端；
(4)对于齿根应变信号获取方法第一步中的第(4)种情况
若n＝np(is)，首先截取一个齿根应变周期中从STA端开始ns(is)个检测齿的齿根应变信号，然后将时间间隔为isNT2的后一个截取信号的STAJ端平移至前一个截取信号的ENDJ端；
若n＝ns(is)，将时间间隔为isNT2的后一个齿根应变周期的STA端平移至前一个齿根应变周期的END端。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610019138.5A CN105823632B (zh) | 2016-01-12 | 2016-01-12 | 行星齿轮箱故障信息获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610019138.5A CN105823632B (zh) | 2016-01-12 | 2016-01-12 | 行星齿轮箱故障信息获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105823632A CN105823632A (zh) | 2016-08-03 |
| CN105823632B true CN105823632B (zh) | 2019-02-05 |
Family
ID=56986894
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610019138.5A Active CN105823632B (zh) | 2016-01-12 | 2016-01-12 | 行星齿轮箱故障信息获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105823632B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106225910B (zh) * | 2016-08-23 | 2019-01-04 | 武汉理工大学 | 基于光纤光栅的转轮振动测量方法及装置 |
| CN108613806A (zh) * | 2018-07-23 | 2018-10-02 | 潍柴动力股份有限公司 | 一种齿轮检测方法及装置 |
| CN111307057B (zh) * | 2020-03-13 | 2022-04-12 | 西安工程大学 | 一种利用光纤光栅检测轴系齿轮齿根应变的方法 |
| CN112231859B (zh) * | 2020-10-11 | 2022-05-17 | 昆明理工大学 | 一种行星齿轮箱行星轮局部故障振动仿真方法 |
| CN112378343B (zh) * | 2020-11-12 | 2022-06-07 | 新疆大学 | 风力发电机行星齿轮箱多路应变动态检测***与方法 |
| CN113761675B (zh) * | 2021-07-23 | 2023-09-22 | 东北大学 | 基于边频分布规律的行星轮轮齿裂纹故障特征判定方法 |
| CN114297801A (zh) * | 2021-12-31 | 2022-04-08 | 重庆大学 | 一种行星轮系内部动态啮合力和动态摩擦力测试方法 |
| CN114266125A (zh) * | 2021-12-31 | 2022-04-01 | 重庆大学 | 一种行星轮系均载分析方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419252A (zh) * | 2011-08-18 | 2012-04-18 | 黑龙江大学 | 光纤在线式高速列车齿轮箱检测装置 |
| CN104101380A (zh) * | 2014-07-02 | 2014-10-15 | 西安交通大学 | 行星轮系固定中心齿轮齿根应力应变的测量方法 |
| DE102014208800A1 (de) * | 2014-05-09 | 2015-11-12 | Schaeffler Technologies AG & Co. KG | Kunststoffplanetenrad mit eingespritztem Lager sowie Planetengetriebe |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9927578B2 (en) * | 2010-09-14 | 2018-03-27 | Trex Enterprises Corp. | Fiber optic rotary connector |
-
2016
- 2016-01-12 CN CN201610019138.5A patent/CN105823632B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419252A (zh) * | 2011-08-18 | 2012-04-18 | 黑龙江大学 | 光纤在线式高速列车齿轮箱检测装置 |
| DE102014208800A1 (de) * | 2014-05-09 | 2015-11-12 | Schaeffler Technologies AG & Co. KG | Kunststoffplanetenrad mit eingespritztem Lager sowie Planetengetriebe |
| CN104101380A (zh) * | 2014-07-02 | 2014-10-15 | 西安交通大学 | 行星轮系固定中心齿轮齿根应力应变的测量方法 |
Non-Patent Citations (2)
| Title |
|---|
| 光纤动态检测技术的研究与进展;张小栋等;《振动、测试与诊断》;20150630;第35卷(第3期);第409-416页 |
| 基于光纤光栅的直齿圆柱齿轮应力修正系数优化;李晗;《机械工程师》;20151010;全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105823632A (zh) | 2016-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105823632B (zh) | 行星齿轮箱故障信息获取方法 | |
| US11220999B1 (en) | Deep hybrid convolutional neural network for fault diagnosis of wind turbine gearboxes | |
| CN104101380B (zh) | 行星轮系固定中心齿轮齿根应力应变的测量方法 | |
| CN202661241U (zh) | 一种模拟风电齿轮箱的振动测试装置 | |
| CN102419252B (zh) | 光纤在线式高速列车齿轮箱检测装置 | |
| CN106368908A (zh) | 一种基于scada***风电机组功率曲线测试方法 | |
| CN105510023A (zh) | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 | |
| CN103604601A (zh) | 基于风电齿轮箱工况模拟的故障诊断试验台 | |
| CN203616120U (zh) | 基于风电齿轮箱工况模拟的故障诊断试验台 | |
| CN101833014A (zh) | 光纤光栅式风电场风速远程检测仪 | |
| CN104792523A (zh) | 基于行星齿轮结构等效轴承模型的振动频谱特征分析方法 | |
| CN103712792A (zh) | 一种针对风电齿轮箱的故障诊断方法 | |
| CN106568597A (zh) | 圆柱齿轮齿面综合啮合刚度的高精度测量方法 | |
| TW201120310A (en) | The state telemetry technology and fault diagnosing system in large-scale wind power farms | |
| CN106055734A (zh) | 一种风电齿轮箱故障诊断方法 | |
| CN103759891A (zh) | 一种双馈风力发电机叶片不平衡在线故障诊断方法 | |
| CN103411727A (zh) | 用于压气机压力测量的光纤压力传感器及其测量方法 | |
| CN101762247A (zh) | 一种光纤光栅位移传感器 | |
| CN104564542A (zh) | 一种基于大数据技术的故障诊断***及其诊断方法 | |
| CN104501772B (zh) | 水工结构土石结合区域不均匀沉降组合监测装置及方法 | |
| Wang et al. | Tacholess order-tracking approach for wind turbine gearbox fault detection | |
| Gutierrez Santiago et al. | Input torque measurements for wind turbine gearboxes using fiber-optic strain sensors | |
| CN104093976A (zh) | 用于确定驱动链的至少一个旋转部件的总受损程度的方法和计算单元 | |
| CN104501773B (zh) | 一种水工建筑物竖向变形监测装置及方法 | |
| CN102592172A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |