一种振动数据驱动的直升机旋翼异常检测方法
【专利摘要】本发明公开一种振动数据驱动的直升机旋翼异常检测方法，包括：原始振动数据分析和处理；提取直升机旋翼状态特征；用正常状态的旋翼数据特征训练支持向量数据模型；最后根据建立的模型对旋翼进行异常检测。本发明的优点是：针对非平稳、非线性的振动数据，提取了有效的故障特征，在只有正常状态数据可用的情况下，本发明给出了一种直升机旋翼异常检测方法，克服了常规方法必须有异常特征数据的需求。
【专利说明】
-种振动数据驱动的直升机旋翼异常检测方法
技术领域
[0001] 本发明设及振动数据驱动的直升机旋翼异常检测方法，是一种专口针对直升机旋 翼振动数据的异常检测方法，属于航空工程与信息科学的交叉领域。
【背景技术】
[0002] 直升机旋翼作为升力面和操纵面，不仅为前进提供动力，还是各种姿态变换的动 力来源。旋翼系统工作在周期性变化的气动环境中，承受了复杂的气弹载荷，如前飞过程中 的不平衡周期的交变载荷，运些载荷容易造成轴承、减摆器等出现磨损现象，更有甚者振动 载荷会导致直升机的重要的飞行组件出现疲劳损伤，使旋翼产生故障进而导致事故的发 生。在直升机故障中旋翼故障率高达20%，其中直升机一级事故中有35%与旋翼故障相关。 因此，对直升机旋翼系统的检测和故障诊断就显得尤为重要。
[0003] 由于直升机旋翼系统的工作特性，直接对旋翼系统进行检测和诊断技术难度较 大。旋翼发生故障时，机体的振动频率中包含了旋翼的故障信息。通过分析采集的机体振动 数据为旋翼系统故障诊断开辟了新途径，经理论分析和实践证明，该方法行之有效且易于 操作。直升机旋翼复杂的结构和工作环境决定了直升机振动数据是非平稳和非线性的，选 用时域或者频域的分析方法很难获取有效的旋翼故障特征。时-频分析方法在两维空间中 同时展示了振动数据的时域和频域信息，常用于提取旋转器械的故障特征。常用的时-频分 析方法有短时傅里叶变化，小波分析，S变换、Cohen类双线性变换方法W及希尔伯特黄变换 等。
[0004] 相比于国外的研究，我国在直升机诊断方面的认识不够全面，研究起步较晚，因此 积极深入直升机故障诊断方面的研究将具有实际意义。目前已经从理论和实践两个方面证 实了旋翼频率与故障之间具有映射关系，利用异常模式识别，可进行故障检测。直升机故障 诊断研究主要集中于旋翼的不平衡故障，有关旋翼异常检测方面的研究很少。传统的故障 诊断一般是二分类问题，主要借助分类算法进行处理，其中人工神经网络和支持向量机已 经取得了不错的诊断效果。然而，运些故障诊断方法的应用前提是有大量的故障数据。在真 实场景中，考虑到直升机旋翼的未知故障模式多种多样，难W预料，无法事先大量获取各种 状态下的故障数据;正常状态的数据很容易获取，故障状态时旋翼可能已经损坏;又因为故 障模拟的代价非常高，在实践中往往不允许，导致故障状态样本很难捕获而致出现严重匿 乏。因此，基于正常状态的振动数据，及时检测出旋翼异常信息显得尤为重要。
【发明内容】
[0005] 本发明的目的是基于振动试验数据，发明一个针对直升机旋翼异常检测的方法， 基于该方法，可检测直升机旋翼状态的变化，从而提供直升机旋翼的早期故障预警，降低旋 翼的维修成本的同时，提高旋翼飞行的安全性和可靠性。
[0006] 本发明的具体技术方案包括W下几个步骤：
[0007] 步骤一:对旋翼振动数据进行预处理，选用小波阔值去噪对原始数据进行去噪处 理。
[0008] 步骤二:对降噪后的振动数据进行经验模态分解，获取有限个内蕴模式函数。
[0009] 步骤=:用获取的内蕴模式函数构建矩阵，并进行奇异值分解，将得到奇异值归一 化到区间[-1,1]，最终组成旋翼故障特征向量。
[0010] 步骤四：利用支持向量数据描述，基于旋翼正常特征构建超球体，计算出超球体的 球屯、和半径，并根据此超球体检测旋翼是否出现异常。
[0011] 本发明的有益效果是：
[0012] 本发明基于直升机旋翼振动数据进行分析，选用时-频分析方法提取有效的旋翼 状态特征，利用支持向量数据描述针对旋翼正常状态特征构建状态模型，解决了在只有正 常状态样本数据可用的情况下，如何检测出旋翼是否出现故障的问题。
【附图说明】
[OOK]图1是本发明方法的总体流程图。
[0014] 图2是经验模态分解的流程图。
[0015] 图3旋翼状态特征提取的流程图。
【具体实施方式】
[0016] 下面结合附图和相关算法，对本发明做进一步的说明。
[0017] 本发明的总体流程如图1所示。
[0018] 本发明基于直升机振动数据进行异常检测，采用小波阔值去噪法对振动数据进行 去噪处理。利用经验模态分解获取有限个内蕴模式函数，对内蕴模式函数矩阵进行奇异值 分解，提取直升机旋翼状态特征。采用支持向量数据描述，针对直升机正常状态特征向量构 建超球体，检测旋翼异常。其具体实施步骤如下：
[0019] 1.小波阔值去噪
[0020] 本发明应用于直升机旋翼振动数据，通过传感器测得的振动数据不可避免地含有 噪声，在特征提取之前选用小波阔值去噪法对振动数据进行去噪处理。小波阔值去噪具体 过程如下：
[0021] 对含噪声的振动数据x(t)做小波变换，得到一组小波系数W;
[0022] 通过对W进行阔值处理，得到估计系数，使得W与两者的差值尽可能小；
[0023] 利用进行小波重构，得到去躁后的振动数据。
[0024] 2.经验模态分解
[0025] 本发明针对去噪后的振动数据进行经验模态分解。经验模态分解能把一个复杂的 振动数据分解为若干不同的内蕴模式函数之和。每一个内蕴模式函数反应了信号内部固有 的波动特性，且只包含一个波动模态，避免了多波动模态的混叠。
[0026] 内蕴模式函数必须满足如下两个条件：
[0027] 在整个数据段内，极值点的个数(包括极大值和极小值)和过零点的个数相等或者 最多只相差一个；
[0028] 在任意时间点上，由局部极大值点形成的上包络线和局部极值点形成的下包络线 的平均值为零。
[0029] 针对振动数据x(t)，经验模态分解具体步骤如下：
[0030] 确定x(t)上所有的极值点，然后用=次样条曲线分别把极大值点和极小值点连起 来，形成上、下包络线；
[003。 求出上下包络线的平均值，记为mi，将mi从x(t)中分离出来
[0032] hi = x(t)-mi (8)
[0033] 如果hi满足内蕴模式函数的两个条件，那么hi就是x(t)分解出来的第一个内蕴模 式函数，若hi不满足内蕴模式函数的条件，则把hi作为原始数据，重复前面的过程直至满足 相应的条件或者预先设置的停止准则，记第一个内蕴模式函数为C1。
[0034] 把Cl从X (t)中分离出来，得到
[0035] ri = x(t)-ci (9)
[0036] 将ri作为原始数据，根据上述方法继续对ri进行分解，直至rn变成一个单调函数， 不能再从中选出满足经验模态函数条件的分量时，循环结束，最终得到
[0037]
(10)
[0038] 其中rn是残余函数，代表振动数据的平均趋势。利用经验模态分解获取的内蕴模 式函数为C1，C2，. . .，Cn。
[0039] 3.奇异值分解
[0040] 经验模态分解利用数据中的特征时间尺度识别出包含其内的内蕴模式函数Cl~ Cn,构建内蕴模式函数矩阵T= [Cl C2…CnL对矩阵T进行奇异值分解。奇异值分解的步骤 如下：
[0041 ] 假设T是M*N的矩阵，U是M*M的矩阵，其中U的列为TTT的正交特征向量，V为N*N的矩 阵，其中V的列为TTT的正交特征向量，m为矩阵T的秩，则存在奇异值分解：
[0042] T = UEyT (11)
[0043] 其中TTT和tTt的特征值相同，记为Al, . . .，Am，E为M*N的对角矩阵，对角元素为
，其余位置数值为0,则奇异值为
I由奇异值构成的奇异值谱为U= [01， 曰2，? ? ?，。皿]O
[0044] 4.数据归一化
[0045] 本发明针对获取的旋翼特征向量U= [01，02, ...，0m]，采用线性函数对各属性进行 归一化处理，其过程可用算法1来描述：
[0046] 算法1:数据归一化
[0047] 输入:旋翼特征向量集{ui}
[0048] 输出：归一化后的数据{V i}
[0049]
[(K)加]5.旋翼正常状态建模
[0051] 本发明选取旋翼正常状态下的特征向量集，训练得到旋翼健康状态模型。
[0052] 其本质是构建一个超球体包含所有或者几乎所有旋翼正常状态下的特征向量集X = {xl，X2，…XN}，此处的x即为提取的特征u，超球体的球屯、为a半径为R。训练的目标是最小 化超球体的体积，降低接受异常样本的概率。考虑到训练特征向量集中可能有少数正常的 特征偏离目标类，通过引入松弛因子Cl和惩罚参数C，允许少数样本点分布在超球体之外。 [0化3] 优化巧敬化下所示：
[0化4] (12)
[0055]为了避免过拟合，引入了惩罚参数C，表示对错分样本惩罚的程度，实现了在错分 样本的比例和区域容量之间的折衷。将上述的述优化问题，转换为其拉格朗日对偶问题进 行求解： 闺
(13)
[0057]其中，Qi为拉格朗日乘子。如果参数设置合理，高斯径向基函数KU响）可W逼近 大多数核函数。本发明采用K(xi，xj)来代替函数内积（(Hxi) ? (Hx如，实现由低维空间到 高维空间的邮射-
[005引 （14)
[0059] 假定支持向量的数量是化v，Xm表示一个支持向量，则超球体球屯、和半径的计算公 式如下：
[0060]
[0061] (16)
[0062] 6.异常检测
[0063] 针对直升机旋翼测试状态特征向量集，根据正常状态建立的模型，计算测试特征 的决策函数值。本发明将决策函数值定义为状态特征向量偏离正常状态模型中屯、a的程度。 若测试特征的决策函数值大于20%，则认为有异常发生。
[0064] 对于测试样本X可计算出其距离球屯、的距离为
[00化]
(17)
[0066] 化笛品激(吿皆f;n 了.
[0067] (18)
[0068] 决策函数的具体计算过程可用算法2来描述：
[0069] 算法2:计算决策函数值
[0070] 输入:待测试的样本X，基于健康状态特征建立的模型mode 1 [0071 ]输出：样本X的决策函数值
[0072]
【主权项】
1. 一种振动数据驱动的直升机旋翼异常检测的方法，其特征包括如下步骤： (1) 针对直升机旋翼振动数据的特点，按异常检测过程的特殊需要对原始数据进行分 析和处理； (2) 选用信号处理领域的时-频分析方法，对振动数据进行经验模态分解获取有限个内 蕴模式函数； (3) 用获取的内蕴模式函数构建矩阵，并进行奇异值分解，将得到奇异值归一化到区间 [_1，1]，最终组成旋翼状态特征向量； (4) 利用支持向量数据描述，选用旋翼正常特征构建超球体，检测旋翼是否出现故障。2. 根据权利要求1所述的振动数据驱动的直升机旋翼异常检测的方法，其特征在于，步 骤(1)针对直升机振动数据的分析和处理，其实现方法如下： (21)选用小波阈值去噪方法对原始的振动数据进行去噪处理。3. 根据权利要求1所述的方法，其特征在于，步骤(2)对去噪后的数据进行经验模态分 解，得到有限个相互独立的内蕴模式函数，获取旋翼振动数据的时频特性，其实现方法包 括： 针对振动数据x(t)，确定x(t)上所有的极值点，然后用三次样条曲线分别把极大值点 和极小值点连起来，形成上、下包络线； 求出上下包络线的平均值，记为m，将m从X (t)中分离出来： hi = x(t)-mi (1) 如果In满足内蕴模式函数的两个条件，那么In就是x(t)分解出来的第一个内蕴模式函 数，若hi不满足内蕴模式函数的条件，则把In看作原始数据，重复前面的过程直至满足相应 的条件或者预先设置的停止准则，记第一个内蕴模式函数为c 1; 把cdAx (t)中分离出来，得到 ri = x(t)-ci (2) 将^作为原始数据，根据以上过程继续对^进行分解，直至rn变成一个单调函数，不能 再从中选出满足内蕴模式函数条件的分量时，循环结束，最终满足如下关系：(3) 其中rn表示残余函数，代表振动数据的平均趋势。利用经验模态分解获取的内蕴模式函 数为C1，Cy，Cn，获取的内蕴模式函数之间两两相互独立。4. 根据权利要求1所述方法，其特征在于，步骤(3)用获取的内蕴模式函数构建矩阵，并 进行奇异值分解，将得到的奇异值归一化到区间[_1，1]，组成旋翼状态特征向量 (41) 构建内蕴模式函数矩阵： 振动数据经经验模态分解后变成了η个相互独立的内蕴模式函数，记做C1~Cn，以内蕴 模式函数为列向量构建矩阵T; T=[ci C2 ··· Cn] (4) (42) 奇异值分解： 对内蕴模式函数矩阵T进行奇异值分解，获取m个奇异值，将奇异值按从大到小的顺序 进行排列，则〇1彡〇2彡...，构成的奇异值谱...，om]代表了旋翼振动数据的 状态特征； (43)归一化处理： 为了消除量纲属性对异常检测的影响，对获取的奇异值进行归一化处理，将数值规约 到[_1，1]之间。归一化函数如下所示：⑶ 式中，σ/表示归一化后的奇异值，分别表示第i个属性的下限和上限，归一化 处理后下限和下限分别变为-1和1。5.根据权利要求1所述方法，其特征在于，步骤(4)利用支持向量数据描述，选用旋翼正 常状态的特征构建超球体，检测直升机旋翼是否出现故障，实现方法如下： 基于旋翼正常状态下提取的特征，利用支持向量数据描述，依据结构风险最小化原理， 对正常状态的特征向量集X ={ XI，X2，…XN}构建超球体，该超球体尽可能多的包含正常状态 的样本特征同时尽可能降低接受异常特征的概率。考虑到训练特征向量集中可能有少数正 常的特征偏离目标类，通过引入松弛因子L和惩罚参数C，允许少数样本点分布在超球体区 域之外。构建超球体可转化为以下优化问题的求解：(6) 通过求解上述优化问题，计算得到超球体的球心是a，半径是R，对于新的样本特征X，决 策函数定义如下：(7) 若样本特征的决策函数值f (X )> 20 %，则认为有异常发生。
【文档编号】G06K9/62GK105956514SQ201610246471
【公开日】2016年9月21日
【申请日】2016年4月15日
【发明人】皮德常, 何云
【申请人】南京航空航天大学