基于叶端定时传感器提取叶片间固有频率差值方法与流程

1.本发明属于叶片非接触无损检测领域，特别是一种基于叶端定时传感器提取叶片间固有频率差值方法。


背景技术：

2.叶片被广泛应用于压气机、燃气轮机、航空发动机等旋转式流体机械设备中，而这些叶片往往在高速、高振动、大冲击等极端恶劣的环境下工作，叶片的故障大都源于异物损伤(foreign object damage，fod)，这些初始缺陷会随着时间的推移而加深，最终导致设备故障甚至发生安全事故。所以对关键设备的旋转叶片进行在线诊断十分必要。叶端定时技术(blade tip timing，btt)是一种非接触在线测量旋转叶片振动的方法，但叶端定时采样速率和转速和传感器数量有关，由于实际情况下传感器安装位置有限，所以叶端定时数据具有严重欠采样特性。叶端定时测量方法在实际使用中，往往选用使用多个叶端定时传感器来减弱欠采样所造成的混叠影响，但多个叶端定时传感器安装在实际有限的空间中会造成较大的困扰，并且也增加了测量成本，传统方法对叶端定时传感器均布和任意布置的数据进行模态参数识别时，如压缩感知、子空间方法、最小二乘迭代。一方面这类算法都涉及大量的运算，无法实现在线实时检测和诊断，另一方面这类算法是直接对单个叶片的固有频率等参数进行识别，存在较大的误差。所以对叶端定时测量方法和模态参数识别方法进行改进十分重要。
3.在背景技术部分中公开的上述信息仅仅用于增强对本发明背景的理解，因此可能包含不构成在本国中本领域普通技术人员公知的现有技术的信息。


技术实现要素：

4.针对现有技术中存在的问题，本发明提出一种基于叶端定时传感器提取叶片间固有频率差值方法，本发明对叶片的健康状态给出更快速和更准确的识别。
5.本发明的目的是通过以下技术方案予以实现，一种基于叶端定时传感器提取叶片间固有频率差值方法包括以下步骤：
6.第一步骤中，利用多个叶端定时传感器获取旋转叶片的实际达到时间，并根据旋转叶片的转速和叶片长度，理论到达时间和实际达到时间之差转换为叶端的位移数据；
7.第二步骤中，对多个叶端定时传感器编号，基于所述编号区分每个叶片的位移数据以分别分析每个叶端定时传感器的位移数据；
8.第三步骤中，选择并截取同转速下的两个旋转叶片的叶端的位移数据；
9.第四步骤中，截取的所述位移数据分别离散傅里叶变换，采样频率近似为平均转速以得到频谱数据，
10.第五步骤中，对每个叶端定时传感器采集到的所述两个旋转叶片的位移数据重复第三步骤和第四步骤，将得到的频谱数据进行线性叠加绘制出总幅频图，
11.第六步骤中，从所述总幅频图中提取两个旋转叶片固有频率混叠后对应的频率余
数，将两个频率余数作差得到两个叶片间的固有频率差。
12.所述的方法中，第一步骤中，单个叶端定时传感器收取每个叶片的达到时刻t，并根据叶片的转速f
r
和叶片长度r将理论到达和实际达到时间差转换为叶端位移，表达式如下：其中表示第i个叶片在第j圈的实际到达第k个传感器的时间，表示第i个叶片在第j圈的理论到达第k个传感器的时间，即叶片无振动时的到达时间，其中θ
i
表示以转速传感器安装位置为基准，第i个叶片的角度。表示以转速传感器安装位置为基准，第k个传感器的角度，n
j
为第j圈时的转速，表示第i个叶片在第时刻的位移。
13.所述的方法中，叶片的旋转过程为预定加速度的升速、减速过程或者匀速过程，在旋转过程中使用周向均布的气嘴喷气模拟气体激励。
14.所述的方法中，对于升速或减速过程，选取叶端定时传感器整数圈记录下的位移数据，利用相同的位置区间[n，m]对数据进行截取，获得两个叶片的同转速下的位移数据，采样频率f
s
近似等于平均转速，其中表示第k圈时的转速，由于以上截取的是索引序号为[n，m]，对于单传感器而言，就是截取了第n到第m圈的数据，这里是在求平均转频。这里的截取是任意的：比如两个长度均为m的位移向量l1，l2，选定索引范围[n1，n2]，即截取l1，l2索引范围为[n1，n2]内的数据。但是l1，l2需要满足是同一次实验中，测量得到的两个叶片的位移，另外需要说明的是[n1，n2]范围是任意确定，如果是对于变速情况下测量测到的数据，该范围不宜过长。
[0015]
所述的方法中，第四步骤中，对所截取两段离散位移数据进行离散傅里叶变换得到频谱数据，
[0016]
其中x(n)为采样得到的信号，i是虚数符号，n为采集到的信号的长度，x中元素的个数，n是一个迭代数，从0遍历到n
‑
1，即取遍x中的所有元素，k是一个0到n
‑
1整数，x(k)表示离散傅里叶变换后的第k个数据。所述的方法中，第五步骤中，得到的n
p
组频谱数据线性叠加绘制出频谱数据am
i，j
，并将频谱数据绘制成幅频图，其中n
p
为传感器数量，
[0017]
所述的方法中，第六步骤中，从总幅频图中提取两个叶片对应的固有频率混叠后的在欠采样频谱中的频率作为频率余数f
sub
，将两个频率余数f
sub
相减得到2个叶片间的固有频率差：将固有频率发生混叠后，在欠采样频谱中呈现的频率
称为频率余数，f
1sub
分别指叶片2和叶片1的频率余数。通过对原本的信号序列进行离散傅里叶变换就可以得到欠采样频谱。(对信号进行傅里叶变换可以得到频谱，只不过btt信号是欠采样信号，所以以下都称为欠采样频谱)。这个离散傅里叶变换就是上面那个式子。在欠采样频谱中提取幅值最高频率成分分量，将该频率作为频率余数f
sub
。
[0018]
本发明方法只需要任意布置的叶端定时传感器即可实现从严重欠采样的数据中提取不同叶片之间的固有频率差，对传感器分布没有特殊的要求，传感器可以任意分布，不需要进行额外的信号重构和更多的叶端定时传感器，运算快速稳定，简单可行，可实现旋转叶片的实时健康监测。
[0019]
上述说明仅是本发明技术方案的概述，为了能够使得本发明的技术手段更加清楚明白，达到本领域技术人员可依照说明书的内容予以实施的程度，并且为了能够让本发明的上述和其它目的、特征和优点能够更明显易懂，下面以本发明的具体实施方式进行举例说明。
附图说明
[0020]
通过阅读下文优选的具体实施方式中的详细描述，本发明各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。说明书附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。显而易见地，下面描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。而且在整个附图中，用相同的附图标记表示相同的部件。
[0021]
在附图中：
[0022]
图1为基于频率余数的任意布置叶端定时叶片间固有频率差值提取步骤示意图；
[0023]
图2为所截取的传感器1提取到的1号和2号叶片位移图；
[0024]
图3为所截取的传感器2提取到的1号和2号叶片位移图；
[0025]
图4为所截取的传感器3提取到的1号和2号叶片位移图；
[0026]
图5为所截取的传感器4提取到的1号和2号叶片位移图；
[0027]
图6为所截取传感器1提取到的1号和2号叶片位移的dft幅频图；
[0028]
图7为所截取传感器2提取到的1号和2号叶片位移的dft幅频图；
[0029]
图8为所截取传感器3提取到的1号和2号叶片位移的dft幅频图；
[0030]
图9为所截取传感器4提取到的1号和2号叶片位移的dft幅频图；
[0031]
图10为4个传感器提取的1号叶片位移频谱数据叠加后的总幅频图；
[0032]
图11为4个传感器提取的2号叶片位移频谱数据叠加后的总幅频图。
[0033]
以下结合附图和实施例对本发明作进一步的解释。
具体实施方式
[0034]
下面将参照附图图1至图11更详细地描述本发明的具体实施例。虽然附图中显示了本发明的具体实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
[0035]
需要说明的是，在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可以理解，技术人员可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名词的差异来作为区分组件的方式，而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语，故应解释成“包含但不限定于”。说明书后续描述为实施本发明的较佳实施方式，然所述描述乃以说明书的一般原则为目的，并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。
[0036]
为便于对本发明实施例的理解，下面将结合附图以具体实施例为例做进一步的解释说明，且各个附图并不构成对本发明实施例的限定。
[0037]
基于叶端定时传感器提取叶片间固有频率差值方法包括，
[0038]
器叶片间固有频率差值提取方法，包括下述步骤：
[0039]
(1)利用任意布置的叶端定时传感器获取旋转叶片的达到时间，并根据转速和叶片长度，将理论到达和实际达到时间差转换为叶端位移。
[0040]
在本示例性实例中，具体为将4个光纤型叶端定时传感器固定在机匣上的随机位置，事后测量可知传感器的安装角度为：48
°
、108
°
、158
°
、168
°
，由于本发明专利所设计的方法并没有使用具体的安装角度信息，所以对传感器安装位置没有要求，将初始转速设定为60hz，转速加速度为0.5hz/s，转速变化范围内60hz
‑
100hz
‑
60hz，其中100hz匀速段时间为20s。叶盘采用8叶片的整体式铝合金叶盘，叶盘半径为r＝68mm，叶片厚度d＝1mm，叶片宽度w＝20mm。在机匣上均布4个喷嘴，喷射0.5mpa的高压气体，利用任意布置的4个叶端定时传感器获取旋转叶片的达到时间，并根据转速和叶片长度，将理论到达和实际达到时间差转换为叶端位移。
[0041]
(2)根据传感器编号将每个叶片的位移数据分开，对每个传感器的位移数据进行分析；
[0042]
本实例使用的叶端定时系统可以直接输出每个传感器通道的叶片实际达到时间，利用转速和叶片半径信息将实际达到时间与理论达到时间之差转换成位移，经过该过程得到了每个传感器通道测量到的所有叶片末端振动位移。
[0043]
(3)选择所要分析的两个叶片的近似同转速下的两段位移数据。若选取的是缓慢升速或降速数据，截取的数据长度不宜过长，以达到近似恒定采样率频的要求。
[0044]
在本示例性实例中，具体为选取叶片1和叶片2的位移数据相同序号范围的240个数据，序号范围为[4719，4958]，每个传感器通道的位移数据选取范围相同，均为[4719，4958]。如图2所示，对应的转速变化范围为：84.66hz～85.59hz，得到单个传感器采集到的位移信号的近似采样频率f
s
＝85.125hz。
[0045]
(4)对截取的两段数据分别进行离散傅里叶变换，采样频率近似为平均转速，得到频谱数据；
[0046]
在本示例性实例中，离散傅里叶变换的计算公式为：
[0047][0048]
其中x(n)为采样得到的信号，i是虚数符号，n为采集到的信号的长度，x
中元素的个数，n是一个迭代数，从0遍历到n
‑
1，即取遍x中的所有元素，k是一个0到n
‑
1整数，x(k)表示离散傅里叶变换后的第k个数据。
[0049]
(5)对每个传感器采集到的两个叶片的数据重复(3)和(4)，将得到的频谱数据进行线性叠加，绘制出总幅频图；
[0050]
在本示例性实例中，具体为对1号、2号、3号、4号叶端定时传感器采集到的1号叶片和2号叶片末端位移数据按照步骤(3)和步骤(4)进行分析，得到叶片1的4个频谱数据向量和叶片2的4个频谱数据向量如图2至图9所示。分别将这4个频谱数据向量进行线性叠加，得到叶片1和2的总幅频数据am
1，2
和am
2，1
：
[0051][0052][0053]
根据总幅频数据am
1，2
和am
2，1
分别绘制4个传感器提取的1号和2号叶片位移频谱数据叠加后总幅频图，如图10和图11所示。
[0054]
(6)从总幅频图中提取两个叶片固有频率混叠后对应的频率余数，将两个频率余数作差即可得到2个叶片间的固有频率差。
[0055]
在本示例性实例中，具体包括以下步骤：
[0056]
a)对上述的总频谱数据am
1，2
和am
2，1
进行分析，绘制出幅频曲线，
[0057]
b)排除掉转速频率分量，由于信号中除了叶片振动频率成分，还存在转频及其倍频，由于对每个传感器通道的位移数据而言，采样频率等于转频，对转频及其倍频的混叠情况进行分析可知，该频率在欠采样频谱中出现在0频率附近，排除掉转速频率分量后，提取叶片振动频率混叠后的在欠采样频谱中的频率，以下称为频率余数f
sub
。
[0058][0059]
对于满足条件的两个叶片，也可以直接得到δf＝δf
sub
。
[0060]
【应用实例】
[0061]
如图1所示的叶端定时试验台，具体为将4个光纤型叶端定时传感器固定在机匣上的随机位置，事后测量可知传感器的安装角度为：48
°
、108
°
、158
°
、168
°
，由于本发明专利所设计的方法并没有使用具体的安装角度信息，所以对传感器安装位置没有要求，将初始转速设定为60hz，转速加速度为0.5hz/s，转速变化范围内60hz
‑
100hz
‑
60hz，其中100hz匀速段时间为20s。叶盘采用8叶片的整体式铝合金叶盘，叶盘半径为r＝68mm，叶片厚度d＝1mm，叶片宽度w＝20mm。在机匣上均布4个喷嘴，喷射0.5mpa的高压气体，利用任意布置的4个叶端定时传感器获取旋转叶片的达到时间，并根据转速和叶片长度，将理论到达和实际达到时间差转换为叶端位移。
[0062]
具体为选取叶片1和叶片2的位移数据相同序号范围的240个数据，序号范围为[4719，4958]，每个传感器通道的位移数据选取范围相同，均为[4719，4958]。如图2所示，对应的转速变化范围为：84.66hz～85.59hz，得到单个传感器采集到的位移信号的近似采样频率f
s
＝85.125hz。
[0063]
对1号、2号、3号、4号叶端定时传感器采集到的1号叶片和2号叶片末端位移数据按照步骤(3)和步骤(4)进行分析，得到叶片1的4个频谱数据向量和叶片2的4个频谱数据向量和叶片2的4个频谱数据向量如图2至图9所示。将这4个频谱数据向量进行线性叠加，得到叶片1和2的总幅频数据am
1，2
和am
2，1
。
[0064]
根据总幅频数据am
1，2
和am
2，1
分别绘制4个传感器提取的1号和2号叶片位移频谱数据叠加后总幅频图，如图10和图11所示。
[0065]
其中分析4个传感器提取的1号叶片位移频谱数据叠加后的总幅频图，即图10可知，叶片1的固有频率混叠后的频率成分为0.35hz，即叶片1的频率余数f
1sub
＝0.35hz。同时可看到图11，即4个传感器提取的2号叶片位移频谱数据叠加后的总幅频图中也含有频率成分0.35hz。观察所截取的1号叶片的位移图可知，所截取的1号叶片振动数据接近共振区，转频的倍频接近于固有频率，可以知道转频及其倍频混叠后的频率也为0.35hz，所以在图11上也会有0.35hz这个转频及其倍频混叠后的频率分量。图10中的0.35hz的频率成分既包括了转频及其倍频混叠后的频率分量，也包含了叶片1的固有频率混叠后的频率分量，所以其幅值远大于图11中的对应频率的幅值。
[0066]
从图11可知，叶片2的固有频率混叠后的频率成分为19.09hz，即叶片2的频率余数
[0067]
通过模态实验分析，可知该铝合金叶盘的1号叶片和2号叶片的一阶固有频率分别为：341hz、361hz，固有频率差为20hz，实际叶片间的固有频率差不会过大，而选择的数据段的采样频率为f
s
＝85.125hz，很明显满足的条件，所以此时通过发明专利提出的频率余数的任意布置叶端定时传感器叶片间固有频率差值提取方法可以得到叶片1和叶片2的频率差值δf＝δf
sub
，进一步可计算由此可得δf＝18.74hz。这与实验测量的20hz非常接近，相差仅为1.26hz，并且该方法不涉及复杂运算，仅需要对2组少量数据进行离散傅里叶变换分析，所以运算快速稳定，简单可行，可实现旋转叶片的实时健康监测。
[0068]
尽管以上结合附图对本发明的实施方案进行了描述，但本发明并不局限于上述的具体实施方案和应用领域，上述的具体实施方案仅仅是示意性的、指导性的，而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下，还可以做出很多种的形式，这些均属于本发明保护之列。