本发明属于故障诊断领域，具体涉及一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法。
背景技术：
通风机是利用叶轮旋转来提升气体压力并输送气体的机械，担负着更换空气，保障工作人员健康的重要任务，广泛应用于矿山、农业、冶金、石油、制药、化工、航天、航海、能源及效能等工程领域。叶片是通风机转子系统的重要零件，运行中工作条件恶劣，载荷复杂。故而，长期运行的风机叶片会产生不同程度的裂纹，甚至引起叶片断裂，对机组的安全运行带来威胁，进而引发重大安全事故。
现有通风机叶片裂纹故障诊断的方法均存在各自的局限性：首先，基于信号处理的诊断方法——现场振动情况复杂且不稳定，分离振动信号与噪声信号仍是技术难题，小尺寸裂纹的信号常会被噪声信号或整体振动信号湮灭，且目前的故障监测系统主要是依据故障特征是否超出设定的阈值来判断是否有故障，尚不能做到定位检测；其次，基于人工智能的故障诊断方法——虽然基于人工智能的故障诊断方法快速方便，但是无法融入经验性的知识，需要经过庞大数量的样本学习才能保障诊断的可靠性，且存在由于规模庞大导致诊断模型学习和训练时间过长的问题，对叶片裂纹的非线性问题没有通用解决方案，难以用于通风机叶片裂纹位置与深度的诊断。
对损伤识别来说，结构局部的损伤往往对应着信号的奇异点。小波分析作为信号处理的有效工具，在时、频两域都具有表征信号局部特征的能力，越来越多地被应用到损伤识别领域。但是，小波分析在通风机叶片裂纹故障诊断中存在一个十分重要的问题就是小波基的选择问题：不同的小波基分析同一个问题会产生不同的结果。因此，小波分析的结果受到所选小波基的限制，通过小波分析得到的小波分量和小波谱只相对于所选择的小波基有意义，而不是直接对应于原始数据。
hilbert—huang变换(hilberthuangtransform，简称hht)是以hilbert变换为基础的一种新的信号处理技术。该方法首先利用经验模式分解(emd)方法将给定信号自适应地逐级分解成若干基本的本征模函数(imf)，通过将这种对信号局部特征的分解方法的引入，使得信号瞬时频率这一概念具有了实际的物理意义；然后对分解得到的本征模函数分量进行hilbert变换，得到相应的时频域上的能量分布即hilbert谱。hilbert—huang变换方法没有固定的先验基底，是自适应的，不同的序列得出不同组的本征模函数，同时，瞬时频率被定义为相位函数的导数，不需要整个波来定义局部频率，因而可以实现从低频信号中分辨出奇异信号，非常适合非线性信号的处理。
技术实现要素：
针对叶片裂纹局部的损伤往往对应着信号奇异点的特点，提出一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法。
本发明采用的技术方案为一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，该方法包括以下步骤：
步骤s1：通过仿真建立通风机叶片裂纹故障的扩充样本数据库。
步骤s2：通过测试系统获取通风机叶片在正常以及不同程度裂纹故障状态下的振型位移参数(即振型位移值)，建立通风机叶片裂纹故障的样本数据库。
步骤s3：利用经验模式分解(emd)方法将样本数据库中的振型位移参数自适应地逐级分解成若干基本的本征模函数(imf)。
步骤s4：然后对分解得到的本征模函数(imf)分量进行空间hilbert变换，得到相应的能量分布即hilbert谱，从图中可以明显看出裂纹截面处奇异点突变,由此便可判定叶片裂纹的位置。
步骤s5：依据振型位移参数中的奇异性大小定义叶片裂纹深度的特征频率,空间峰值采样频率的极大值和分析频率带宽的线性关系得到有效拟合直线的斜率。有效拟合直线斜率与裂纹深度之间存在比例关系，由此实现叶片裂纹深度的定量判定。
步骤s6：将扩充本数据库中的振型位移参数输入到裂纹故障诊断模型中，对叶片裂纹故障进行诊断。
本发明提出一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，该方法将hht变换的优越性与通风机叶片结构损伤识别(位置、深度)的要求联系起来，对hilbert—huang变换方法进行相应改动，将其从时域转化为空间域，识别出损伤位置，以及定义归一化的损伤频率来表征叶片裂纹的深度，实现了通风机叶片裂纹位置与深度的故障诊断。本方法能够得到裂纹位置及深度的定量数据，操作起来更加简单,且保证一定的准确率，实现了对风机叶片裂纹位置的定量分析,诊断速度快,并且保证了一定的故障诊断率，可靠性比较高。
附图说明
图1是本发明的整体步骤流程图。
具体实施
为了详细说明本发明的主要特点及具体操作步骤，现配合附图进行说明。
步骤s1：通过仿真建立通风机叶片裂纹故障的扩充样本数据库。
具体进行如下处理：
步骤s1.1：通过仿真建立通风机叶片裂纹模型，设置多种程度的叶片裂纹故障。
步骤s1.2：通过有限元仿真对叶片裂纹故障模型进行模态分析，得到仿真数据，建立通风机叶片裂纹故障的扩充样本数据库。
步骤s2：通过测试系统获取通风机叶片在正常以及不同程度裂纹故障状态下的振型位移参数(即振型位移值)，建立通风机叶片裂纹故障的数据库。
具体进行如下处理：
步骤s2.1：测试过程为离线测试模式，测试系统由传感器部分、激励部分、数据采集与分析部分组成,传感器部分通过激励部分与数据采集与分析部分进行连接。其中,传感器部分为压电式力传感器和压电式加速度传感器。
步骤s2.2：激励部分为激振器或冲击锤装置；数据采集与分析部分要针对锤击、环境激励识别方法的不同,采用相应的数据采集与分析应用软件。
步骤s2.3：包括传函分析、参数选择、输入几何结构及约束信息、模态拟合和振型编辑,最后获得位移模态参数，依照通风机叶片的状态——正常、裂纹位置、裂纹深度，适当地进行数据分组，建立通风机叶片裂纹故障的数据库。
步骤s3：利用经验模式分解(emd)方法将振型位移参数自适应地逐级分解成若干基本的本征模函数(imf)。
具体进行如下处理：
步骤s3.1：分别获通风机叶片在正常、裂纹位置、裂纹深度不同故障状态下各n(n≥10)组振型位移参数信号，这三种状态分别得到n组信号样本集，分别记为xmn(t)，其中m＝1，2，3；n＝1，2，3，……，n。
步骤s3.2：对得到的振型位移参数样本进行经验模式分解，得到m个本征模函数(简称imf)分量，将每个imfmn分量记为cmni(t)，则
其中j为imfmn分量的个数，rj(t)为趋势余量。
步骤s4：然后对分解得到的本征模函数(imf)分量进行空间hilbert变换，得到相应的位频能量分布谱即横坐标为x的hilbert谱，从图中可以明显看出裂纹截面处奇异点突变,由此便可判定叶片裂纹的位置。
具体进行如下处理：
步骤s4.1：对于空间hilbert-huang变换，需要把hilbert-huang变换中的t替换为x，将其从时域转化为空间域，其空间采样频率可以类似的按照实际的空间坐标的分辨率来取，即取为空间采样间隔的倒数。
步骤s4.2：对振型数据进行空间hilbert-huang变换,得到hilbert谱；图中横坐标的值为0至叶片总长度为l，纵坐标为由于奇异性所带来的峰值变化，从图中能够看出上裂纹截面的处奇异点的突变,由此判定叶片裂纹的位置。
步骤s5：依据振型位移参数中的奇异性大小定义叶片裂纹深度的特征频率,空间峰值采样频率的极大值和分析频率带宽的线性关系得到有效拟合直线的斜率。有效拟合直线斜率与裂纹深度之间存在比例关系，由此实现叶片裂纹深度的定量判定。
具体进行如下处理：
步骤s5.1：对于叶片损伤故障诊断来说，空间频率峰值的大小可以很好的反映信号的奇异性。定义损伤特征频率，即：
步骤s5.2：式中fx即为空间的采样频率，fx反映了分析频率的带宽max[ωi(x)]表示空问频率的极大值，f反映了空间频率的极大值在分析频率带宽中所占的比例，以此来表征损伤程度——叶片裂纹深度。
步骤s6：将扩充本数据库中的振型位移参数输入到裂纹故障诊断模型中，对叶片裂纹故障进行诊断。
1.一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：该方法包括以下步骤：
步骤s1：通过仿真建立通风机叶片裂纹故障的扩充样本数据库；
步骤s2：通过测试系统获取通风机叶片在正常以及不同程度裂纹故障状态下的振型位移参数即振型位移值，建立通风机叶片裂纹故障的数据库；
步骤s3：利用经验模式分解emd方法将振型位移参数自适应地逐级分解成若干基本的本征模函数imf；
步骤s4：然后对分解得到的本征模函数imf分量进行空间hilbert变换，得到相应的位频能量分布谱即横坐标为x的hilbert谱，从hilbert谱中得出裂纹截面处奇异点突变，判定叶片裂纹的位置；
步骤s5：依据振型位移参数中的奇异性大小定义叶片裂纹深度的特征频率,空间峰值采样频率的极大值和分析频率带宽的线性关系得到有效拟合直线的斜率；拟合直线斜率与裂纹深度之间存在比例关系，由此实现叶片裂纹深度的定量判定；
步骤s6：将扩充本数据库中的振型位移参数输入到裂纹故障诊断模型中，对叶片裂纹故障进行诊断。
2.根据权利要求1所述的一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：步骤s1的具体实施过程如下，
步骤s1.1：通过仿真建立通风机叶片裂纹模型，设置多种程度的叶片裂纹故障；
步骤s1.2：通过有限元仿真对叶片裂纹故障模型进行模态分析，得到仿真数据，建立通风机叶片裂纹故障的扩充样本数据库。
3.根据权利要求1所述的一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：步骤s2的具体实施过程如下，
步骤s2.1：测试过程为离线测试模式，测试系统由传感器部分、激励部分、数据采集与分析部分组成，传感器部分通过激励部分与数据采集与分析部分进行连接；传感器部分为压电式力传感器和压电式加速度传感器；
步骤s2.2：激励部分为激振器或冲击锤装置；数据采集与分析部分要针对锤击、环境激励识别方法的不同,采用相应的数据采集与分析应用软件；
步骤s2.3：包括传函分析、参数选择、输入几何结构及约束信息、模态拟合和振型编辑,最后获得位移模态参数，依照通风机叶片的状态——正常、裂纹位置、裂纹深度，进行数据分组，建立通风机叶片裂纹故障的数据库。
4.根据权利要求1所述的一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：步骤s3的具体实施过程如下，
步骤s3.1：分别获通风机叶片在正常、裂纹位置、裂纹深度不同故障状态下各n组振型位移参数信号，这三种状态分别得到n组信号样本集，分别记为xmn(t)，其中m＝1，2，3；n＝1，2，3，……，n，n≥10；
步骤s3.2：对得到的振型位移参数样本进行经验模式分解，得到m个本征模函数分量，将每个imfmn分量记为cmni(t)，则
其中j为imfmn分量的个数，rj(t)为趋势余量。
5.根据权利要求1所述的一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：步骤s4的具体实施过程如下，
步骤s4.1：对于空间hilbert-huang变换，把hilbert-huang变换中的t替换为x，将其从时域转化为空间域，其空间采样频率按照实际的空间坐标的分辨率来取，即取为空间采样间隔的倒数；
步骤s4.2：对振型数据进行空间hilbert-huang变换,得到hilbert谱；横坐标的值为0至叶片总长度为l，纵坐标为由于奇异性所带来的峰值变化，得到裂纹截面的处奇异点的突变,由此判定叶片裂纹的位置。
6.根据权利要求1所述的一种基于hilbert—huang变换的通风机叶片裂纹故障诊断方法，其特征在于：步骤s5的具体实施过程如下，
步骤s5.1：对于叶片损伤故障诊断来说，空间频率峰值的大小反映信号的奇异性；定义损伤特征频率，即：
步骤s5.2：式中fx即为空间的采样频率，fx反映了分析频率的带宽max[ωi(x)]表示空问频率的极大值，f反映了空间频率的极大值在分析频率带宽中所占的比例，以此来表征损伤程度——叶片裂纹深度。
技术总结