基于行星齿轮结构等效轴承模型的振动频谱特征分析方法技术

本发明专利技术公开了一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法；其具体包括以下步骤：输入原始振动信号并利用提升小波包变换方法进行分解，计算分解得到的每层分解频带信号的峭度值、得到峭度分析图并选取最大峭度值，获取对应的分解频带信号并利用希尔伯特包络解调分析方法进行分析得到振动频谱图，建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型，对行星齿轮传动系统各部件故障振动频率进行精确定量分析。本发明专利技术具有精度高、简单可靠的优点，实现了理论计算的故障频率峰值和实际故障频率的峰值直接对应，从而精确表征了行星齿轮箱传动系统各部件的故障特征频率。

【技术实现步骤摘要】
基于行星齿轮结构等效轴承模型的振动频谱特征分析方法
本专利技术属于旋转机械装备故障诊断与预示
，尤其涉及一种基于行星齿轮结构等效轴承模型的行星齿轮箱传动系统中各部件故障振动频谱特征分析方法。
技术介绍
行星齿轮传动系统是机械传动系统中的关键部分，其广泛应用于航空航天、船舶、风力发电、工程机械等领域。行星齿轮传动系统通常工作在大功率、高转速、高负载的工作环境下，行星齿轮系极容易发生齿面磨损、齿面接触疲劳、轮齿弯曲疲劳、乃至断齿或轴断裂等失效现象，最终导致整个系统的完全失效，造成严重的经济损失，甚至人员伤亡和灾难性的损失。列如：行星齿轮传动系统是风力发电系统的关键装置之一，长期在无规律变向变载荷的风力作用以及强阵风的瞬时冲击下工作，导致其故障频率极高。但由于传统的监测手段无法及时发现故障，造成了严重的经济损失，成为了制约风力发电技术发展和推广的棘手难题。2007年5月，美国军方在对UH-60A“黑鹰”直升机的一次常规检测中，意外发现了未被监测装置探测出的行星齿轮传动系统中行星架的严重裂纹，而这种故障将导致机毁人亡，为此美军下令停飞和检测所有在役直升机。目前行星齿轮传动系统由于特殊的物理结构和复合运动形式导致了故障振动机理比与普通的平行齿轮传动系统更为复杂，其振动信号具有复杂时变调制特点。而现有的行星齿轮振动频谱特征研究模型较复杂，主要对啮合频率的边频带分析，理论计算和实际故障频率的峰值之间存在偏差，而且观测到的故障边频带太多，对于故障位置判断可能产生误判。
技术实现思路
本专利技术的专利技术目的是：为了解决现有技术中行星齿轮振动频谱特征研究模型复杂且对故障位置可能产生误判等问题，本专利技术提出了一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法。本专利技术的技术方案是：一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法，包括以下步骤：A、输入行星齿轮的原始振动信号，利用提升小波包变换方法对原始振动信号进行分解；B、对步骤A中分解得到的每层分解频带信号计算其峭度值，得到峭度分析图，选取最大峭度值；C、根据步骤B中得到的最大峭度值获取其对应的分解频带信号，利用希尔伯特包络解调分析方法进行分析，得到振动频谱图；D、建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型；E、通过步骤D中建立的基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型计算振动频谱特征，并对步骤C中得到的振动频谱图中的振动频谱特征进行精确定量分析。进一步地，所述步骤A中利用提升小波包变换方法对原始振动信号进行分解，具体为：设定行星齿轮的原始振动信号为X(t)，利用提升小波包变换方法将原始振动信号进行m层分解，得到每层2m分解频带。进一步地，所述m的取值为3或4。进一步地，所述步骤B中峭度的表达式具体为：其中，E[·]为求期望符号，x为采集样本，μ为样本均值，σ为样本标准差。进一步地，所述步骤B中得到峭度分析图的方法具体为：对计算得到的峭度值进行归一化处理，以各层对应的各节点信号频带为横坐标，以分解信号层数为纵坐标，得到峭度分析图。进一步地，所述步骤D中建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型，具体包括以下分步骤：D1、获取行星齿轮箱中太阳轮齿数Zs、行星轮齿数Zp、内齿圈齿数Zr、行星轮个数zp、太阳轮转动频率fns及轴承中由局部缺陷引起冲击的间隔振动频谱特征：内圈故障振动频率fbi为：外圈故障振动频率fbo为：滚动体故障振动频率fbb为：其中，f为轴的转频，α为接触角，z为滚动体数目，d为滚动体直径，D为节径、即滚动体中心到滚动轴承中心距离的两倍；D2、利用齿轮分度圆直径定义：df＝mZ对行星齿轮箱中太阳轮和行星轮进行处理，其中，df为分度圆直径，m为模数，Z为齿轮齿数；D3、将行星齿轮箱中太阳轮故障振动频率、行星轮故障振动频率、内齿圈故障振动频率、分别对应轴承中内圈故障振动频率、滚动体故障振动频率、外圈故障振动频率，得到行星齿轮传动系统各部件的故障振动频谱特征：太阳轮故障振动频率fs为：内齿圈故障振动频率fr为：行星轮故障振动频率fp为：进一步地，所述步骤E中通过步骤D建立的基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型计算振动频谱特征，具体为：先计算基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型中一个部件的故障振动频谱特征，再通过其它部件故障振动频谱特征与该部件故障振动频谱特征的比值关系得到其它部件的故障振动频谱特征。本专利技术具有以下有益效果：(1)本专利技术建立的基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型，实现了理论计算的故障频率峰值和实际故障频率的峰值直接对应，从而精确表征了行星齿轮箱传动系统各部件的故障特征频率。(2)本专利技术实现了基于提升小波包变换、峭度、希尔伯特包络解调分析方法和行星齿轮结构等效轴承模型的振动频谱特征理论计算模型的算法上混合，具有精度高、简单可靠的优点，为精确定量提取行星齿轮传动系统故障诊断提供了一种新思路和新方法，具有广阔的工程应用前景。附图说明图1为本专利技术的基于行星齿轮结构等效轴承模型的振动频谱特征分析方法流程示意图。图2为本专利技术的实施例中得到的峭度分析图示意图。图3为本专利技术的行星齿轮结构等效轴承模型的振动频谱特征计算模型示意图。图4为基于希尔伯特包络解调的太阳轮磨损故障振动频谱特征分析图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，为本专利技术的基于行星齿轮结构等效轴承模型的振动频谱特征分析方法流程示意图。一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法，包括以下步骤：A、输入行星齿轮的原始振动信号，利用提升小波包变换方法对原始振动信号进行分解；B、对步骤A中分解得到的每层分解频带信号计算其峭度值，得到峭度分析图，选取最大峭度值；C、根据步骤B中得到的最大峭度值获取其对应的分解频带信号，利用希尔伯特包络解调分析方法进行分析，得到振动频谱图；D、建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型；E、通过步骤D中建立的基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型计算振动频谱特征，并对步骤C中得到的振动频谱图中的振动频谱特征进行精确定量分析。在步骤A中，输入行星齿轮的原始振动信号，设为X(t)，再利用提升小波包变换方法将原始振动信号进行m层分解，得到每层2m分解频带；优选的，m的取值为3或4。这里的提升小波包变换方法为本领域技术人员常用技术手段，本专利技术不作赘述。在步骤B中，对步骤A分解得到的每层分解频带信号分别计算其峭度值，从而得到峭度分析图，再根据峭度分析图选取得到最大峭度值。这里得到峭度分析图的方法具体为：对计算得到的峭度值进行归一化处理，以各层对应的各节点信号频带为横坐标，以分解信号层数为纵坐标，从而得到峭度分析图。如图2所示，为本专利技术的实施例中得到的峭度分析图示意图。为了更直观的对峭度进行分析，本专利技术采用不同的色度对应不同的峭度值；色度越接近红色，表示峭度值越大；色度越接近蓝色，表示峭度值越小。在步骤C中，根据步骤B得到的最大峭度值获取其对应的分解频带信号，特别的，若最大峭度值对应2个本文档来自技高网...

【技术保护点】
一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法，其特征在于，包括以下步骤：A、输入行星齿轮的原始振动信号，利用提升小波包变换方法对原始振动信号进行分解；B、对步骤A中分解得到的每层分解频带信号计算其峭度值，得到峭度分析图，选取最大峭度值；C、根据步骤B中得到的最大峭度值获取其对应的分解频带信号，利用希尔伯特包络解调分析方法进行分析，得到振动频谱图；D、建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型；E、通过步骤D中建立的基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型计算振动频谱特征，并对步骤C中得到的振动频谱图中的振动频谱特征进行精确定量分析。

【技术特征摘要】
1.一种基于行星齿轮结构等效轴承模型的振动频谱特征分析方法，其特征在于，包括以下步骤：A、输入行星齿轮的原始振动信号，利用提升小波包变换方法对原始振动信号进行分解；B、对步骤A中分解得到的每层分解频带信号计算其峭度值，得到峭度分析图，选取最大峭度值；C、根据步骤B中得到的最大峭度值获取其对应的分解频带信号，利用希尔伯特包络解调分析方法进行分析，得到振动频谱图；D、建立基于行星齿轮结构等效轴承模型的振动频谱特征理论计算模型，具体包括以下分步骤：D1、获取行星齿轮箱中太阳轮齿数Zs、行星轮齿数Zp、内齿圈齿数Zr、行星轮个数zp、太阳轮转动频率fns及轴承中由局部缺陷引起冲击的间隔振动频谱特征：内圈故障振动频率fbi为：外圈故障振动频率fbo为：滚动体故障振动频率fbb为：其中，f为轴的转频，α为接触角，z为滚动体个数，d为滚动体直径，D为节径、即滚动体中心到滚动轴承中心距离的两倍；D2、利用齿轮分度圆直径定义：df＝mZ对行星齿轮箱中太阳轮和行星轮进行处理，其中，df为分度圆直径，m为模数，Z为齿轮齿数；D3、将行星齿轮箱中太阳轮故障振动频率、行星轮故障振动频率、内齿圈故障振动频率分别对应轴承中内圈故障振动频率、滚动体故障振动频率、外圈故障振动频率，得到行星齿轮传动系统各部件的故障振动频谱特征：太阳轮故障振动频率fs为：内齿圈故障振动频率fr为：行星轮故障振动频率fp为：

【专利技术属性】
技术研发人员：刘治汶，刘剑锋，左明健，郭伟，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51