CN110940517A - 一种直升机传动系统的振动信号分离方法 - Google Patents

一种直升机传动系统的振动信号分离方法 Download PDF

Info

Publication number
CN110940517A
CN110940517A CN201910963096.4A CN201910963096A CN110940517A CN 110940517 A CN110940517 A CN 110940517A CN 201910963096 A CN201910963096 A CN 201910963096A CN 110940517 A CN110940517 A CN 110940517A
Authority
CN
China
Prior art keywords
vibration
data
gear
planet
sensor
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.)
Pending
Application number
CN201910963096.4A
Other languages
English (en)
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.)
China Helicopter Research and Development Institute
Original Assignee
China Helicopter Research and Development Institute
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 China Helicopter Research and Development Institute filed Critical China Helicopter Research and Development Institute
Priority to CN201910963096.4A priority Critical patent/CN110940517A/zh
Publication of CN110940517A publication Critical patent/CN110940517A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明实施例提供一种直升机传动系统的振动信号分离方法,传动系统包括齿圈,多个行星轮,振动传感器,转速传感器,方法包括:根据振动传感器和转速传感器的安装位置,获取传动系统中行星轮系的位置信息;对齿圈的齿轮啮合时序进行分析,得到振动传感器与行星轮的啮合顺序;对振动数据进行数据点的采样,并对行星轮架的每个转动周期内的数据点进行重采样;对重采样得到的振动数据的数据点进行数据提取,得到四维组合集;对四维组合集进行振动分离向量重构,得到齿圈的振动特征信号。本发明实施例可以从复杂的行星传动机构的振动信号中,针对多个振动传感器采集的数据,提取出完整的齿圈振动信号,为行星传动机构的齿圈故障诊断奠定基础。

Description

一种直升机传动系统的振动信号分离方法
技术领域
本申请涉及但不限于直升机检测技术领域,尤指一种直升机传动系统的振动信号分离方法。
背景技术
目前,国内外关于直升机传动系统中,行星传动机构中齿圈故障诊断的研究主要集中在谱分析方法、时频分析法、粒子滤波方法等方面,虽然提出了一些新的状态指标,并且针对行星轮系信号特征复杂、相互淹没的特点进行了信号增强检测方法研究。
但是总体来看,这些研究仍然将行星传动机构中的齿圈当做普通的轮齿啮合形式进行分析,并未深入到各个部件,在实际中仍然无法有效解决信号成分相互淹没的问题。
发明内容
本发明的目的:本发明实施例提供一种直升机传动系统的振动信号分离方法,以转速脉冲信号为依据,实现从振动信号中分离重构齿圈的振动特征信号,从而有效的提取出齿圈啮合时冲击产生的振动信号。
本发明的技术方案:
本发明实施例提供一种直升机传动系统的振动信号分离方法,所述传动系统包括齿圈,设置于所述齿圈内、且与所述齿圈啮合转动的多个行星轮,用于支撑所述行星轮的行星轮架,设置于所述齿圈周围的至少一个振动传感器,以及设置于齿圈周围、用于测量行星轮架转动的转速传感器,所述振动信号分离方法包括:
根据所述振动传感器和所述转速传感器的安装位置,获取所述传动系统中行星轮系的位置信息,所述行星轮系的位置信息包括每个振动传感器的位置信息和每个行星轮的位置信息;
根据所述行星轮系的位置信息,对所述齿圈的齿轮啮合时序进行分析,得到每个振动传感器与每个行星轮的啮合顺序;
在所述传动系统的传动过程中对振动数据进行数据点的采样,并对所述行星轮架的每个转动周期内的数据点进行重采样,所述振动数据为行星轮与振动传感器对齐时的振动信号;
对重采样得到的振动数据的数据点进行数据提取,通过提取出每个行星轮与每个振动传感器对齐时对应的振动数据得到四维组合集,所述四维组合集包括振动传感器编号,行星轮编号,行星架循环编号,以及振动数据;
对所述四维组合集进行振动分离向量重构,得到所述齿圈的振动特征信号,所述齿圈的振动特征信号包括每个行星轮与每个振动传感器对齐时的振动信号。
可选地,如上所述的直升机传动系统的振动信号分离方法中,
所述每个振动传感器的位置信息为:
Figure RE-GDA0002368192240000021
所述每个行星轮的位置信息为:
Figure RE-GDA0002368192240000022
其中,所述
Figure RE-GDA0002368192240000023
为转速传感器与第j个振动传感器的夹角,所述
Figure RE-GDA0002368192240000024
为转速传感器与第i个行星轮的夹角。
可选地,如上所述的直升机传动系统的振动信号分离方法中,所述根据所述行星轮系的位置信息,对所述齿圈的齿轮啮合时序进行分析,得到每个振动传感器与每个行星轮的啮合顺序,包括:
根据所述行星轮系的位置信息,获取每个行星轮与每个传感器的分离角矩阵为:
Figure RE-GDA0002368192240000025
其中,所述D中的其中一个分离角为:
Figure RE-GDA0002368192240000026
所述
Figure RE-GDA0002368192240000027
为第i个行星轮从初始位置转动到第j个振动传感器时,行星轮架的转动角度;
根据所述分离角度矩阵,获取每个振动传感器与每个行星轮的啮合顺序。
可选地,如上所述的直升机传动系统的振动信号分离方法中,所述对所述行星轮架的每个转动周期内的数据点进行重采样,包括:
对所述行星轮架的每个转动周期内的数据点进行重采样,使得所述行星轮架在每次循环周期中采集的数据点数相同,重采样的数据点数为:
Figure RE-GDA0002368192240000031
其中,Nr为所述齿圈的轮齿数量,也为行星轮架在每次循环周期中行星轮与齿圈的啮合周期,所述NTP为齿圈的每个轮齿对应的数据点数。
可选地,如上所述的直升机传动系统的振动信号分离方法中,对重采样得到的振动信号数据点进行数据提取,通过提取出每个行星轮与每个振动传感器对齐时对应的振动信号得到四维组合集,包括:
确定数据点的提取位置,并根据所述提取位置得到振动数据的提取长度;
对每个振动传感器编号,每个行星轮编号和每个行星架循环编号的组合,对应一个振动数据,得到所述四维组合集,所述振动数据的长度为根据所述提取长度提取的。
可选地,如上所述的直升机传动系统的振动信号分离方法中,所述确定数据点的提取位置,并根据所述提取位置得到振动数据的提取长度,包括:
所述提取位置为所采集的数据点的提取索引:
Figure RE-GDA0002368192240000032
其中,
Figure RE-GDA0002368192240000033
所述Δθc为重采样得到的两个相邻数据点之间的间隔,所述提取索引为第i个行星轮与第j个振动传感器对齐时刻的振动数据;
根据所述提取索引确定所述振动传感器所提取的振动数据的中间位置、起始位置和终止位置,得到所述振动数据的提取长度为:
Figure RE-GDA0002368192240000034
可选地,如上所述的直升机传动系统的振动信号分离方法中,所述对所述四维组合集进行振动分离向量重构,得到所述齿圈的振动特征信号,包括:
对所述四维组合集进行时间同步平均、加窗映射和求和重构,得到所述齿圈的振动特征信号。
可选地,如上所述的直升机传动系统的振动信号分离方法中,对所述四维组合集进行时间同步平均、加窗映射和求和重构,得到所述齿圈的振动特征信号,包括:
所述时间同步平均包括:对行星架在多次循环中每个固定位置提取的振动数据进行平均处理,得到振动传感器所采集到的振动数据的时域同步平均TSA 信号集,所述TSA信号集包括多个固定位置对应的TSA信号,所述固定位置为行星轮与传感器的对齐位置;
所述加窗映射包括:对所述TSA信号集中的每个TSA信号分配匹配索引
Figure RE-GDA0002368192240000041
并通过所述匹配索引确定振动数据的数据位置;所述
Figure RE-GDA0002368192240000042
所述振动数据的数据位置为RANGE:
Figure RE-GDA0002368192240000043
其中,
Figure RE-GDA0002368192240000044
k为行星架的循环次数;
通过所述匹配索引对所述TSA信号进行加窗映射,得到组合保持矩阵,所述组合保持矩阵包括:振动传感器编号,行星轮编号,以及每个振动传感器编号和每个行星轮编号的组合所对应的振动数据;
所述求和重构包括:对组合保持矩阵进行求和,得到振动分离向量重构的结果。
本发明的有益效果:本发明实施例提供一种直升机传动系统的振动信号分离方法,该方法采用传动系统和相应的传感器采集针对信号,并对振动信号进行采集、提取和重构,具体实施方式为:根据振动传感器和转速传感器的安装位置,获取传动系统中行星轮系的位置信息;对齿圈的齿轮啮合时序进行分析,得到振动传感器与行星轮的啮合顺序;对振动数据进行数据点的采样,并对行星轮架的每个转动周期内的数据点进行重采样;对重采样得到的振动数据的数据点进行数据提取,得到四维组合集;对四维组合集进行振动分离向量重构,得到齿圈的振动特征信号。采用本发明实施例提供的直升机传动系统的振动信号分离方法,可以从复杂的行星传动机构的振动信号中,针对多个振动传感器采集的数据,提取出完整的齿圈振动信号,为行星传动机构的齿圈故障诊断奠定基础。
附图说明
图1为本发明实施例提供的一种直升机传动系统的振动信号分离方法的流程图;
图2为传动系统的结构示意图;
图3为振动传感器A1到A5采集的原始振动信号的示意图;
图4为振动传感器A1到A5重采样得到的振动信号的示意图;
图5为在行星架的第二至第五次循环中振动传感器A1数据提取与组合的过程;
图6为振动传感器A1到A5在行星架的一次循环中提取的振动数据;
图7为对行星轮在多次循环中提取的信号进行信号重构的重构信号图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
从上世纪80年代开始,英美等国一直积极发展直升机健康与使用监测系统(简称为:HUMS),对齿轮箱、旋翼和发动机等关键部件进行健康监测、故障诊断与预测。MarcosE.Orchard等人使用粒子滤波方法提出一种可以用于在线故障诊断和预测的方法,他们所提出的方法在UH-60直升机行星齿轮上进行故障实验进行验证,结果表明该方法能够成功指出UH-60的异常情况。
但是,上述分析方法仍然将行星传动机构中的齿圈当做普通的轮齿啮合形式进行分析,并未深入到各个部件,在实际中仍然无法有效解决信号成分相互淹没的问题。
本发明实施例提供的直升机传动系统的振动信号分离方法,针对行星传动机构中齿圈故障诊断提出的信号处理方法,目的是以转速脉冲信号为依据,从振动信号中分离重构齿圈的振动特征信号,有效的提取出齿圈啮合时冲击产生的振动信号。
图1为本发明实施例提供的一种直升机传动系统的振动信号分离方法的流程图,图2为传动系统的结构示意图,该传动系统包括:齿圈,设置于齿圈内、且与齿圈啮合转动的多个行星轮Pi,用于支撑行星轮的行星轮架(图2中未示出),设置于齿圈周围的至少一个振动传感器Aj,以及设置于齿圈周围、用于测量行星轮架转动的转速传感器OPR。
基于传动系统的结构基础,针对该传动系统的振动信号分离方法可以包括如下步骤:
S110,根据振动传感器和转速传感器的安装位置,获取传动系统中行星轮系的位置信息,行星轮系的位置信息包括每个振动传感器的位置信息和每个行星轮的位置信息;
S120,根据行星轮系的位置信息,对齿圈的齿轮啮合时序进行分析,得到每个振动传感器与每个行星轮的啮合顺序;
S130,在传动系统的传动过程中对振动数据进行数据点的采样,并对行星轮架的每个转动周期内的数据点进行重采样,振动数据为行星轮与振动传感器对齐时的振动信号;
S140,对重采样得到的振动数据的数据点进行数据提取,通过提取出每个行星轮与每个振动传感器对齐时对应的振动数据得到四维组合集,该四维组合集包括振动传感器编号,行星轮编号,行星架循环编号,以及振动数据;
S150,对四维组合集进行振动分离向量重构,得到齿圈的振动特征信号,齿圈的振动特征信号包括每个行星轮与每个振动传感器对齐时的振动信号。
本发明实施例提供振动信号分离方法,主要是针对齿圈振动信号的分离,的目的是在原始振动信号中提取出齿圈的动态行为。齿圈振动信号分离的基本前提是,当行星轮Pi与振动传感器Aj对齐时,振动传感器Aj获取的振动信号由行星轮Pi对应的行星轮-齿圈啮合副所主导。由于行星轮围绕太阳轮公转运动,引起行星轮Pi啮合区振动信号在振动传感器Aj中周期性强弱变化,当行星轮Pi与振动传感器Aj对齐时,两者距离最近,振动信号由于传递路径引起的衰减最小。
由于振动传感器Aj固定安装在齿圈上,每个振动传感器Aj只能有效感知到安装位置附近的轮齿啮合产生的振动信号。采用多个振动传感器布置在齿圈外围,每个振动传感器感知一定数量齿圈轮齿的信号,然后从每个振动传感器 Aj中提取相对应的轮齿的信号,剔除远离振动传感器Aj的轮齿啮合振动信号,重构出整个齿圈轮齿的振动信号,有效的实现了齿圈振动信号的分离。
振动信号分离方法的主要思想是求解在行星架的第n次循环中,行星轮Pi 与齿圈上安装的振动传感器Aj的对准时间和对齐时参与啮合的齿圈轮齿编号,根据对准时间在振动传感器Aj信号进行振动信号提取,然后以参与啮合的齿圈轮齿编号为依据,使用提取出的数据重构齿圈的振动信号。
采用本发明实施例提供的直升机传动系统的振动信号分离方法,可以从复杂的行星传动机构的振动信号中,针对多个振动传感器采集的数据,提取出完整的齿圈振动信号,为行星传动机构的齿圈故障诊断奠定基础。
本发明实施例的具体实现方式为:
1,参数设置,目的为了得到上述S110中的位置信息;
振动信号分离技术的参数主要包含两类:位置参数与结构参数,其中位置参数是指初始时刻行星轮相对于转速传感器的角位置以及齿圈上振动传感器安装处相对转速传感器的角位置,结构参数是指传动机构中行星轮的个数和齿数,以及齿圈的齿数。
位置参数确定必须在一定的坐标系下进行,如图2所示。定义齿圈和行星轮的齿数分别为Nr和Np。行星架上有NP个行星轮,沿着齿圈共安装NA个加速度计传感器,转速传感器OPR安装在齿圈的外侧,位置为图2中OPR (once-per-revolution)所示。从转速传感器开始,将行星轮按顺时针方向依次编号为Pi,(i=1,2…NP),振动传感器逆时针依次编号为Aj,(j=1,2…NA)。定义转速传感器OPR和振动传感器Aj的夹角为
Figure RE-GDA0002368192240000071
转速传感器OPR和行星轮 Pi的夹角为
Figure RE-GDA0002368192240000072
则行星轮系的位置参数可用向量θA与θP表示。
根据上述参数设定可以得到:
每个振动传感器的位置信息为:
Figure RE-GDA0002368192240000073
每个行星轮的位置信息为:
Figure RE-GDA0002368192240000074
其中,
Figure RE-GDA0002368192240000075
为转速传感器与第j个振动传感器的夹角,
Figure RE-GDA0002368192240000076
为转速传感器与第 i个行星轮的夹角。
2,轮齿啮合时序分析
为了实现提取后振动信号的精确重构,需要对齿圈与行星轮的啮合时序进行分析,上述S120可以包括:
S121,根据行星轮系的位置信息,获取每个行星轮与每个传感器的分离角矩阵;
S122,根据分离角度矩阵,获取每个振动传感器与每个行星轮的啮合顺序。
在具体实现中,如图2所示的行星轮系坐标图,转速传感器和传感器Aj 的夹角为
Figure RE-GDA0002368192240000077
转速传感器和行星轮Pi的夹角为
Figure RE-GDA0002368192240000078
定义行星轮Pi与传感器Aj 分离角度为:
Figure RE-GDA0002368192240000079
式中mod为求余函数,由于行星轮可能循环了多圈;
Figure RE-GDA00023681922400000710
表示行星轮Pi从初始位置运动到传感器Aj过程中,行星架转过的角度,即相遇时刻在一个行星架循环中的相位,定义分离角矩阵为:
Figure RE-GDA0002368192240000081
其中,D中的其中一个分离角为:
Figure RE-GDA0002368192240000082
为第i个行星轮从初始位置转动到第j个振动传感器时,行星轮架的转动角度;通过上述矩阵,可以确定出齿圈与行星轮的啮合顺序。
3,数据采集与重采样
采集就是直接振动传感器直接采集的振动数据,数据重采样的目的是确保行星架在每次循环中采集的振动数据的数据点数相同。
上述S130可以包括:对行星轮架的每个转动周期内的数据点进行重采样,使得行星轮架在每次循环周期中采集的数据点数相同,重采样的数据点数为: N=Nr*NTP;
其中,Nr为所述齿圈的轮齿数量,也为行星轮架在每次循环周期中行星轮与齿圈的啮合周期,NTP为齿圈的每个轮齿对应的数据点数。
在行星架的每一次循环中,有Nr个行星轮-齿圈的轮齿啮合周期。由于传动系统输入轴转速的轻微变化,行星架的每次循环中采集到的振动信号点数不一致,需要进行重采样,以使得每个行星架循环过程中数据点数相同。重采样点数
Figure RE-GDA0002368192240000083
尽量接近原始数据点数,
Figure RE-GDA0002368192240000084
可由上述公式
Figure RE-GDA0002368192240000085
得出,NTP为每个齿圈的每个轮齿对应的数据点数,且为正整数。
4,数据提取
数据提取是从预处理后的振动信号中提取出每个行星轮与每个振动传感器对齐时对应的振动数据并保存在四维组合集中,四维组合集是一个四维矩阵,用于保存提取数据以及数据的相关信息。数据提取与组合需要明确数据提取位置、提取长度、提取数据对应的行星轮编号。
上述S140可以包括:
S141,确定数据点的提取位置,并根据提取位置得到振动数据的提取长度;
S142,对每个振动传感器编号,每个行星轮编号和每个行星架循环编号的组合,对应一个振动数据,得到四维组合集,振动数据的长度为根据提取长度提取的。
1)确定提取位置,上述S141中,确定的数据点的提取位置为所采集的数据点的提取索引:
Figure RE-GDA0002368192240000091
其中,
Figure RE-GDA0002368192240000092
Δθc为重采样得到的两个相邻数据点之间的间隔,提取索引
Figure RE-GDA0002368192240000093
为第i个行星轮与第j个振动传感器对齐时刻的振动数据。
由上述实施例已知:行星轮角Pi的初始角度为
Figure RE-GDA0002368192240000094
振动传感器Aj的安装位置角度为
Figure RE-GDA0002368192240000095
行星轮Pi与传感器Aj的分离角度
Figure RE-GDA0002368192240000096
定义为
Figure RE-GDA0002368192240000097
式中mod为模函数。
行星架的每个周期经过重采样,有
Figure RE-GDA0002368192240000098
个数据点,每两个点之间的位置间隔为
Figure RE-GDA0002368192240000099
数据处理中取与Pi与Aj对准时刻距离的点,该数据点的索引是:
Figure RE-GDA00023681922400000910
式中round函数是最近取整函数。
2)确定提取长度,上述S141中,根据提取索引确定述振动传感器所提取的振动数据的中间位置、起始位置和终止位置,得到振动数据的提取长度为:
Figure RE-GDA00023681922400000911
由于齿圈相对于振动传感器位置不变,因此,在齿圈振动信号分离过程中,数据长度lE参数确定方法为:使用提取索引
Figure RE-GDA00023681922400000912
确定振动传感器Aj中数据提取的中间位置,数据起始位置由
Figure RE-GDA00023681922400000913
决定,数据终止位置由
Figure RE-GDA00023681922400000914
决定。即从相邻两个振动传感器的几何位置中点截取。因此,对于振动传感器Aj,齿圈振动信号分离的数据提取长度lE为:
Figure RE-GDA00023681922400000915
3)四维组合集,上述S142,对每个振动传感器编号,每个行星轮编号和每个行星架循环编号的组合,对应一个振动数据,得到四维组合集,振动数据的长度为根据提取长度提取的。
思维组合集是用于存放提取数据及相关信息的四维矩阵,该四维矩阵的前三个维度分别表示为:振动传感器编号、行星轮编号和行星架循环编号,这三个维度唯一确定的列向量用于存放对应的提取数据。
5,振动分离向量重构
振动分离向量重构的目的是在组合集的基础上,重构得到齿圈的振动特征信号,重构得到的振动特征信号称为振动分离向量。
S150可以包括:对四维组合集进行时间同步平均、加窗映射和求和重构,得到齿圈的振动特征信号。
即振动分离向量重构的步骤为:时间同步平均、加窗映射,以及求和重构。
1)时间同步平均
时间同步平均的步骤包括:对行星架在多次循环中每个固定位置提取的振动数据进行平均处理,得到振动传感器所采集到的振动数据的时域同步平均(简称为:TSA)信号集,该TSA信号集包括多个固定位置对应的TSA信号,该固定位置为行星轮与传感器的对齐位置。
由于在振动传感器Aj采集的数据子集中,该数据子集中Pi层对应的是通过行星轮Pi在行星架的多次循环中提取的数据。由于齿圈相对于振动传感器是固定的,因此,行星架的每次循环中提取的数据对应相同的齿圈轮齿。对行星架的多次循环的提取数据进行平均处理,得到振动传感器Aj采集的数据子集的 TSA信号。
2)加窗映射
加窗映射的步骤包括:Stp1,对TSA信号集中的每个TSA信号分配匹配索引
Figure RE-GDA0002368192240000101
并通过匹配索引确定振动数据的数据位置;该
Figure RE-GDA0002368192240000102
振动数据的数据位置为:RANGE:
Figure RE-GDA0002368192240000103
其中,
Figure RE-GDA0002368192240000104
k为行星架的循环次数。
Stp2,通过匹配索引对TSA信号进行加窗映射,得到组合保持矩阵,该组合保持矩阵包括:振动传感器编号,行星轮编号,以及每个振动传感器编号和每个行星轮编号的组合所对应的振动数据。
实际上,加窗映射过程是数据提取的逆过程,目的是将平均处理后的TSA 信号映射到振动分离向量中的相应位置,映射通过对TSA信号分配匹配索引
Figure RE-GDA0002368192240000105
实现,由于齿圈相对于传感器是固定的,提取信号的相位即映射目标位置,因此,在齿圈振动分离向量重构中,振动传感器Aj数据子集的TSA信号在对应列向量中的位置由匹配索引
Figure RE-GDA0002368192240000106
唯一确定,且匹配索引为
Figure RE-GDA0002368192240000107
数据映射范围(即振动数据的数据位置)为:
RANGE:
Figure RE-GDA0002368192240000108
式中,
Figure RE-GDA0002368192240000109
k为行星架的循环次数。
式中floor为向下取整函数。匹配索引将各传感器中Pi层得到的TSA信号加窗后映射到组合保持矩阵中。
3)求和重构
求和重构的步骤包括:对组合保持矩阵进行求和,得到振动分离向量重构的结果。
得到组合保持矩阵后,对组合保持矩阵进行行求和,实现振动分离向量重构。
以下通过一个应用实例对本发明上述实施例提供的直升机传动系统的振动信号分离方法进行详细说明:
在工程应用中,针对我们某型号常用的直升机主减速器,采用的是行星齿轮传动系统。其中行星轮3个,齿数23,太阳轮1个,齿数46,齿圈1个,齿数92。在主减速器上安装有5个振动传感器,1个转速传感器。主减速器的信号采集频率为20k/s,对主减速器采集的信号进行振动分离。
第一步:确定参数设置
根据传感器的安装位置和行星轮的安装位置,可以得出行星轮系的位置参数为:
传感器位置参数:θA=[45°,90°,180°,225°,270°];
行星轮位置参数为:θP=[180°,300°,60°];
第二步:啮合时序分析
通过计算可得,分离矩阵为:
Figure RE-GDA0002368192240000111
在行星架的第一次循环内,与5个传感器对应的齿圈轮齿啮合的行星轮编号为:
Figure RE-GDA0002368192240000112
基于上述分析过程,行星架的其他循环周期的啮合顺序可以此类推。
第三步:数据采集与重采样
原始信号长度例如为400000点。行星架的每个周期的重采样点数
Figure RE-GDA0002368192240000122
例如为 3936点。图3所示为振动传感器A1到A5采集的原始振动信号的示意图,图4 所示为振动传感器A1到A5重采样得到的振动信号的示意图。
第四步:数据提取
以振动传感器A1的振动信号为例,图5所示为在行星架的第二至第五次循环中振动传感器A1数据提取与组合的过程,每个循环中四个截取位置分别对应行星轮每次循环的数据提取区域,第二次循环为cycle2,第三次循环为 cycle3,第4次循环为cycle4,第五次循环为cycle5;其中心表示提取索引
Figure RE-GDA0002368192240000121
宽度表示数据提取长度lE。图6所示为振动传感器A1到A5在行星架的一次循环中提取的振动数据。
第五步:振动分离向量重构
对每个振动传感器中提取的信号进行重构,可以重构行星轮每次循环的信号,如图7所示,为对行星轮在多次循环中提取的信号进行信号重构的重构信号图。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种直升机传动系统的振动信号分离方法,其特征在于,所述传动系统包括齿圈,设置于所述齿圈内、且与所述齿圈啮合转动的多个行星轮,用于支撑所述行星轮的行星轮架,设置于所述齿圈周围的至少一个振动传感器,以及设置于齿圈周围、用于测量行星轮架转动的转速传感器,所述振动信号分离方法包括:
根据所述振动传感器和所述转速传感器的安装位置,获取所述传动系统中行星轮系的位置信息,所述行星轮系的位置信息包括每个振动传感器的位置信息和每个行星轮的位置信息;
根据所述行星轮系的位置信息,对所述齿圈的齿轮啮合时序进行分析,得到每个振动传感器与每个行星轮的啮合顺序;
在所述传动系统的传动过程中对振动数据进行数据点的采样,并对所述行星轮架的每个转动周期内的数据点进行重采样,所述振动数据为行星轮与振动传感器对齐时的振动信号;
对重采样得到的振动数据的数据点进行数据提取,通过提取出每个行星轮与每个振动传感器对齐时对应的振动数据得到四维组合集,所述四维组合集包括振动传感器编号,行星轮编号,行星架循环编号,以及振动数据;
对所述四维组合集进行振动分离向量重构,得到所述齿圈的振动特征信号,所述齿圈的振动特征信号包括每个行星轮与每个振动传感器对齐时的振动信号。
2.根据权利要求1所述的直升机传动系统的振动信号分离方法,其特征在于,
所述每个振动传感器的位置信息为:
Figure RE-FDA0002368192230000011
所述每个行星轮的位置信息为:
Figure RE-FDA0002368192230000012
其中,所述
Figure RE-FDA0002368192230000013
为转速传感器与第j个振动传感器的夹角,所述
Figure RE-FDA0002368192230000014
为转速传感器与第i个行星轮的夹角。
3.根据权利要求2所述的直升机传动系统的振动信号分离方法,其特征在于,所述根据所述行星轮系的位置信息,对所述齿圈的齿轮啮合时序进行分析,得到每个振动传感器与每个行星轮的啮合顺序,包括:
根据所述行星轮系的位置信息,获取每个行星轮与每个传感器的分离角矩阵为:
Figure RE-FDA0002368192230000021
其中,所述D中的其中一个分离角为:
Figure RE-FDA0002368192230000022
所述
Figure RE-FDA0002368192230000023
为第i个行星轮从初始位置转动到第j个振动传感器时,行星轮架的转动角度;
根据所述分离角度矩阵,获取每个振动传感器与每个行星轮的啮合顺序。
4.根据权利要求3所述的直升机传动系统的振动信号分离方法,其特征在于,所述对所述行星轮架的每个转动周期内的数据点进行重采样,包括:
对所述行星轮架的每个转动周期内的数据点进行重采样,使得所述行星轮架在每次循环周期中采集的数据点数相同,重采样的数据点数为:
Figure RE-FDA0002368192230000024
其中,Nr为所述齿圈的轮齿数量,也为行星轮架在每次循环周期中行星轮与齿圈的啮合周期,所述NTP为齿圈的每个轮齿对应的数据点数。
5.根据权利要求4所述的直升机传动系统的振动信号分离方法,其特征在于,对重采样得到的振动信号数据点进行数据提取,通过提取出每个行星轮与每个振动传感器对齐时对应的振动信号得到四维组合集,包括:
确定数据点的提取位置,并根据所述提取位置得到振动数据的提取长度;
对每个振动传感器编号,每个行星轮编号和每个行星架循环编号的组合,对应一个振动数据,得到所述四维组合集,所述振动数据的长度为根据所述提取长度提取的。
6.根据权利要求5所述的直升机传动系统的振动信号分离方法,其特征在于,所述确定数据点的提取位置,并根据所述提取位置得到振动数据的提取长度,包括:
所述提取位置为所采集的数据点的提取索引:
Figure RE-FDA0002368192230000025
其中,
Figure RE-FDA0002368192230000026
所述Δθc为重采样得到的两个相邻数据点之间的间隔,所述提取索引为第i个行星轮与第j个振动传感器对齐时刻的振动数据;
根据所述提取索引确定所述振动传感器所提取的振动数据的中间位置、起始位置和终止位置,得到所述振动数据的提取长度为:
Figure RE-FDA0002368192230000027
7.根据权利要求6所述的直升机传动系统的振动信号分离方法,其特征在于,所述对所述四维组合集进行振动分离向量重构,得到所述齿圈的振动特征信号,包括:
对所述四维组合集进行时间同步平均、加窗映射和求和重构,得到所述齿圈的振动特征信号。
8.根据权利要求7所述的直升机传动系统的振动信号分离方法,其特征在于,对所述四维组合集进行时间同步平均、加窗映射和求和重构,得到所述齿圈的振动特征信号,包括:
所述时间同步平均包括:对行星架在多次循环中每个固定位置提取的振动数据进行平均处理,得到振动传感器所采集到的振动数据的时域同步平均TSA信号集,所述TSA信号集包括多个固定位置对应的TSA信号,所述固定位置为行星轮与传感器的对齐位置;
所述加窗映射包括:对所述TSA信号集中的每个TSA信号分配匹配索引
Figure RE-FDA0002368192230000031
并通过所述匹配索引确定振动数据的数据位置;所述
Figure RE-FDA0002368192230000032
所述振动数据的数据位置为
Figure RE-FDA0002368192230000033
其中,
Figure RE-FDA0002368192230000034
k为行星架的循环次数;
通过所述匹配索引对所述TSA信号进行加窗映射,得到组合保持矩阵,所述组合保持矩阵包括:振动传感器编号,行星轮编号,以及每个振动传感器编号和每个行星轮编号的组合所对应的振动数据;
所述求和重构包括:对组合保持矩阵进行求和,得到振动分离向量重构的结果。
CN201910963096.4A 2019-10-11 2019-10-11 一种直升机传动系统的振动信号分离方法 Pending CN110940517A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910963096.4A CN110940517A (zh) 2019-10-11 2019-10-11 一种直升机传动系统的振动信号分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910963096.4A CN110940517A (zh) 2019-10-11 2019-10-11 一种直升机传动系统的振动信号分离方法

Publications (1)

Publication Number Publication Date
CN110940517A true CN110940517A (zh) 2020-03-31

Family

ID=69906040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910963096.4A Pending CN110940517A (zh) 2019-10-11 2019-10-11 一种直升机传动系统的振动信号分离方法

Country Status (1)

Country Link
CN (1) CN110940517A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108362492A (zh) * 2018-01-11 2018-08-03 中国人民解放军国防科技大学 一种适用于低转速下行星轮系故障诊断的振动分离方法
CN108896303A (zh) * 2018-06-19 2018-11-27 昆明理工大学 一种行星齿轮或太阳轮齿根裂纹故障特征的检测方法
CN109443768A (zh) * 2018-11-14 2019-03-08 中国直升机设计研究所 一种直升机主减速器行星齿轮振动信号分离方法
US20190078975A1 (en) * 2017-09-12 2019-03-14 General Electric Company Gearbox sensor arrangement
CN109632295A (zh) * 2019-01-08 2019-04-16 中国人民解放军国防科技大学 基于连续振动分离和最小熵解卷积的行星轮故障检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190078975A1 (en) * 2017-09-12 2019-03-14 General Electric Company Gearbox sensor arrangement
CN108362492A (zh) * 2018-01-11 2018-08-03 中国人民解放军国防科技大学 一种适用于低转速下行星轮系故障诊断的振动分离方法
CN108896303A (zh) * 2018-06-19 2018-11-27 昆明理工大学 一种行星齿轮或太阳轮齿根裂纹故障特征的检测方法
CN109443768A (zh) * 2018-11-14 2019-03-08 中国直升机设计研究所 一种直升机主减速器行星齿轮振动信号分离方法
CN109632295A (zh) * 2019-01-08 2019-04-16 中国人民解放军国防科技大学 基于连续振动分离和最小熵解卷积的行星轮故障检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张伦: "直升机传动系统行星轮系故障仿真与振动分离方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Similar Documents

Publication Publication Date Title
CN108362492B (zh) 一种适用于低转速下行星轮系故障诊断的振动分离方法
CN108195587B (zh) 一种电机滚动轴承故障诊断方法及其诊断系统
CN109443768A (zh) 一种直升机主减速器行星齿轮振动信号分离方法
JP5757936B2 (ja) 回転部材を含んでいる機械的なアッセンブリにおける構造的な異常の検出方法
US20170067860A1 (en) Control Area Network Machine Diagostic
CN110987438B (zh) 水轮发电机变转速过程周期性振动冲击信号检测的方法
CN108896303B (zh) 一种行星齿轮或太阳轮齿根裂纹故障特征的检测方法
CN109187023B (zh) 一种汽车发电机轴承故障诊断方法
CN109632295A (zh) 基于连续振动分离和最小熵解卷积的行星轮故障检测方法
Guo et al. Tooth root crack detection of planet and sun gears based on resonance demodulation and vibration separation
CN105675113A (zh) 基于微型传感器的旋转机械角域振动信号获取装置和方法
CN109682597A (zh) 一种齿轮箱振动信号处理和分析方法
CN111413095A (zh) 基于瞬时角速度的行星轴承分布式故障诊断分析方法
KR101622264B1 (ko) 유성기어의 진동신호 추출 장치 및 추출 방법
CN111855211A (zh) 一种用于故障特征增强的最优解调频带确定方法
CN109238717A (zh) 一种基于vmd-omeda的齿轮箱复合故障诊断方法
CN109839272B (zh) 基于故障冲击提取与自相关总体平均的轴承故障诊断方法
Hu et al. Extraction of the largest amplitude impact transients for diagnosing rolling element defects in bearings
CN108106838B (zh) 一种机械故障精密诊断的整周期时域回归及报警方法
EP0889314B1 (en) Method of monitoring a transmission assembly of a vehicle equipped with acceleration sensors, in particular a helicopter
CN110937135A (zh) 一种直升机减速器的振动信号分离方法
CN110940517A (zh) 一种直升机传动系统的振动信号分离方法
Zeng et al. Planetary gear fault diagnosis based on an instantaneous angular speed measurement system with a dual detector setup
CN112345247B (zh) 一种滚动轴承的故障诊断方法及装置
CN116150585A (zh) 基于乘积包络谱的旋转机械故障诊断方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200331