CN106092524B - 一种使用振动信号精确提取转速信号的方法 - Google Patents
一种使用振动信号精确提取转速信号的方法 Download PDFInfo
- Publication number
- CN106092524B CN106092524B CN201610322710.5A CN201610322710A CN106092524B CN 106092524 B CN106092524 B CN 106092524B CN 201610322710 A CN201610322710 A CN 201610322710A CN 106092524 B CN106092524 B CN 106092524B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- speed
- vibration signal
- tach
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明的目的在于提供一种使用振动信号精确提取转速信号的方法,用于进行针对性长时间的取证,更加方便和隐蔽。为了实现所述目的,本发明一种使用振动信号精确提取转速信号的方法将可以通过获得的过零点时间信息对振动信号——包括本通道或是其他同步采样的通道——进行信号重采样并平均,以获得转速相关的阶次信息,可用于转子***或齿轮箱的故障分析。在没有转速计模拟或数字输入信号的情况下,本办法是一个非常实用的工具。
Description
技术领域
本发明涉及一种提取转速信号用于数据分析的方法,尤其涉及一种使用振动信号精确提取转速信号的方法。
背景技术
旋转机械构成了机械设备中非常重要的部分,比如发电机、压缩机、飞机引擎等。随着科技的快速发展,旋转机械日益大型化和复杂化,一旦发生故障,损失可能非常严重,因此对于旋转机械的基于振动信号的状态监测变得日益重要。
针对旋转机械的状态监测,机械设备的旋转速度和转速信号(Tachometer或键相信号Key phaser)是极其重要的,除了旋转速度本身是一项重要指标外,转速信号也是后续很多振动参数提取和故障分析方法的基础,比如阶次分析,就需要转速信号作为参考时间序列来进行振动信号的重采样。然而工业现场转速信号的测量并不容易,很多设备上面都没有转速计或转速测量装置,而且后续加装也非常困难;另外,无线振动传感器在工业现场的设备监测也越来越多地得以使用,这些传感器是没有转速信号配合的,即使有,它们和转速信号的同步目前也是行业面临的难题之一。
对此,专利号为201410036882.7的发明专利《一种基于DSP的发动机转速测量方法和装置》中就公开了一种基于DSP的发动机转速测量方法和装置,通过获取发动机振动信号、排气声音信号或发动机点火电压信号,并对获得的信号进行AD转换,根据预设的采样率对AD转换后的信号数据进行采样;对采样后的数据进行限幅数字滤波和加汉宁窗函数处理;对处理后的数据进行快速傅氏变换使时域采样数据变换为频域内频谱值,并计算各频率点幅值;根据频率基波间隔法计算发动机对应的频率值;根据发动机转速与发动机对应的频率值的关系,得到发动机的转速;从而提高了转速测量精度及抗干扰能力,降低了误判率。
但是这样的方法虽然能精确获得具体的某点的数值,并不能用于进行后续的故障分析等操作。
发明内容
本发明的目的在于提供一种使用振动信号精确提取转速信号的方法,用于进行针对性长时间的取证,更加方便和隐蔽。
为了实现所述目的,一种使用振动信号精确提取转速信号的方法,包括
步骤1:获得振动信号;
步骤2:对步骤1中所得振动信号加载汉宁窗;
步骤3:对步骤2中所得信号进行快速傅立叶变换,其特征在于,还包括步骤4:估算已知转速范围;
步骤5:将步骤3所得数据结合步骤4中所估算的转速范围,利用汉宁窗特性获得转速频率位置以及复核转速频率;假定左右两个边频带的幅值分别为a,b,位置为fa,fb,则在a≠b的情况下,准确的转速频率为
步骤6:将步骤3中所得数据结合步骤5中获取的转速频率位置提取转速及周边频率幅值,保留相位信息,合并为复数形式;
步骤7:将步骤6中获得的数据进行反快速傅立叶变换,获得时域信号;
其中,ei是复数的指数形式;
n为反快速傅立叶变换后时域数字信号的索引;
k为频域信号的索引。
步骤8:将步骤7中的时域信号通过最小二乘法正弦曲线拟合获得转速频率分量;
步骤9:生成转速方波信号:使用所述步骤8中获取的正弦信号的过零点时间信息作为转速信号的跳变点,生成转速方波信号,为后续的信号分析备用。
优选的,所述步骤1中获得振动信号的具体方式为:假设设备在采样信号长度的时间区间内转速稳定,获取振动传感器的原始数字信号X(n),截取2^N个采样点,使获得的数字信号时间长度超过但接近1秒钟。这样的取值便于后期操作,也足够确保运算的精确度。
优选的,对获得原始数字信号X(n)加载汉宁窗w(n)
窗函数的数据和原始数字信号的数据作对应点之间的乘法,获得加载窗函数后的新的信号数据点X′(n),
X′(n)=w(n)*X(n)。
优选的,所述步骤3中所述快速傅立叶变换方式为:
优选的,所述步骤6中具体步骤为将频域信号转化为极坐标形式,由幅值谱和相位谱组成,把幅值谱单列出来进行下一步处理,相位谱保持不变。
优选的,所述步骤4中已知转速范围估算方法为:通过设备型号获得恒定转速R,在幅值谱上查询到其所对应的峰值;如果转速处于变动状态当中,则根据前次循环处理的结果,在正负5%的频率范围内查询幅值谱的就近峰值。这样便于计算时的取值。
优选的,所述步骤8中获取的正弦信号的周期的倒数和步骤5中获得的复核转速频率相比较,获取差值,所述差值在1RPM之内则校验为合格,否则返回到步骤1,重置整个过程。通过复核转速频率避免过大的误差产生。
通过这样的设置,本发明一种使用振动信号精确提取转速信号的方法将可以通过获得的过零点时间信息对振动信号——包括本通道或是其他同步采样的通道——进行信号重采样并平均,以获得转速相关的阶次信息,可用于转子***或齿轮箱的故障分析。在没有转速计模拟或数字输入信号的情况下,本办法是一个非常实用的工具。
附图说明
图1为本发明一种使用振动信号精确提取转速信号的方法实施例1中原始时域信号示意图。
图2为对图1原始时域信号加载汉宁窗后获得的示意图。
图3为对图2中信号进行傅立叶变换所得的频域频谱示意图。
图4为对已知转速范围进行放大后,并根据汉宁窗的特性对转速峰值进行比例计算获得准确的转速频率示意图。
图5为对照图4显示的通过几何运算获得精确的转速频率位置,图中加粗的黑线为计算所得位置。
图6为编辑频域幅值谱后的示意图,具体编辑方式为:选取中间峰值和左右各三个边频带幅值线保留不变,把其余的谱线幅值设为10的-5次方。
图7为把幅值谱和分离的相位谱合并,做反傅立叶变换,然后对所获得的时域信号做最小二乘法正弦曲线拟合,获得的正弦信号示意图。其中虚线为拟合后的正弦曲线(显示0.3-0.6秒部分)
图8为利用拟合后的正弦曲线的过零点,生成转速方波信号,可用于后续的分析函数的计算。图中曲线为生成转速方波信号,虚线为生成的转速信号。
图9为本发明一种使用振动信号精确提取转速信号的方法具体流程示意图。
具体实施方式
如图9所示一种使用振动信号精确提取转速信号的方法,一种使用振动信号精确提取转速信号的方法,包括:
步骤1:获得振动信号;
步骤2:对步骤1中所得振动信号加载汉宁窗;
步骤3:对步骤2中所得信号进行快速傅立叶变换,其特征在于,还包括步骤4:估算已知转速范围;
步骤5:将步骤3所得数据结合步骤4中所估算的转速范围,利用汉宁窗特性获得转速频率位置以及复核转速频率;假定左右两个边频带的幅值分别为a,b,位置为fa,fb,则在a≠b的情况下,准确的转速频率为
步骤6:将步骤3中所得数据结合步骤5中获取的转速频率位置提取转速及周边频率幅值,保留相位信息,合并为复数形式;
步骤7:将步骤6中获得的数据进行反快速傅立叶变换,获得时域信号;
其中,ei是复数的指数形式;
n为反快速傅立叶变换后时域数字信号的索引;
k为频域信号的索引。
步骤8:将步骤7中的时域信号通过最小二乘法正弦曲线拟合获得转速频率分量;
步骤9:生成转速方波信号:使用步骤8中所述正弦信号的过零点时间信息作为转速信号的跳变点,生成转速方波信号,为后续的信号分析备用。
所述步骤1中获得振动信号的具体方式为:假设设备在采样信号长度的时间区间内转速稳定,获取振动传感器的原始数字信号X(n),截取2^N个采样点,使获得的数字信号时间长度超过但接近1秒钟。对获得原始数字信号X(n)加载汉宁窗w(n)
窗函数的数据和原始数字信号的数据作对应点之间的乘法,获得加载窗函数后的新的信号数据点X′(n),
X′(n)=w(n)*X(n)。所述步骤3中所述快速傅立叶变换方式为:所述步骤6中具体步骤为将频域信号转化为极坐标形式,由幅值谱和相位谱组成,把幅值谱单列出来进行下一步处理,相位谱保持不变。所述步骤4中已知转速估算方法为:通过设备型号获得恒定转速R,在幅值谱上查询到其所对应的峰值;如果转速处于变动状态当中,则根据前次循环处理的结果,在正负5%的频率范围内查询幅值谱的就近峰值。所述步骤8中获取的正弦信号的周期的倒数和步骤5中获得的复核转速频率相比较,获取差值,所述差值在1RPM之内则校验为合格,否则返回到步骤1,重置整个过程。
如图1所示,某个电机拖动的冷却风扇的振动信号数据采集的采样频率为2560Hz,则选取2^12=4096个采样点,即采样长度为4096。这样信号的采样长度为1.6秒,频域的分辨率为2560/4096=0.625Hz。根据型号查询,风扇的转动速度大约为1000转/分。
将图1中的信号加载汉宁窗,获得图2所示的信号。
对图2中所得加窗信号进行傅立叶变换,得到图3所示的频域频谱。图中仅显示0-80Hz部分。
对已知转速范围进行放大,并根据汉宁窗的特性对转速峰值进行比例计算获得准确的转速频率17.08Hz。
图4中即为放大的频谱区间
通过几何运算获得精确的转速频率位置,即图5中所示的加粗的黑线为计算所得位置。
选取中间峰值和左右各三个边频带幅值线保留不变,把其余的谱线幅值设为10的-5次方。
编辑后的频域幅值谱如图6所示,
此后把幅值谱和分离的相位谱合并,做反傅立叶变换。然后对所获得的时域信号做最小二乘法正弦曲线拟合,获得图7所示的正弦信号。虚线为显示0.3-0.6秒部分的拟合的曲线。利用拟合后的正弦曲线的过零点,生成如图8所示的转速方波信号,虚线为生成的转速信号。该曲线可用于后续的分析函数的计算。该方波信号的频率和频谱计算得到的转速频率进行比较,结果一致。
Claims (7)
1.一种使用振动信号精确提取转速信号的方法,包括
步骤1:获得振动信号;
步骤2:对步骤1中所得振动信号加载汉宁窗;
步骤3:对步骤2中所得信号进行快速傅立叶变换,
其特征在于,还包括步骤4:估算已知转速范围;
步骤5:将步骤3所得数据结合步骤4中所估算的转速范围,利用汉宁窗特性获得转速频率位置以及复核转速频率;假定左右两个边频带的幅值分别为a,b,位置为fa,fb,则在a≠b的情况下,准确的转速频率为
步骤6:将步骤3中所得数据结合步骤5中获取的转速频率位置提取转速及周边频率幅值,保留相位信息,合并为复数形式;
步骤7:将步骤6中获得的数据进行反快速傅立叶变换,获得时域信号;
其中,ei是复数的指数形式;
n为反快速傅立叶变换后时域数字信号的索引;
k为频域信号的索引;
步骤8:将步骤7中的时域信号通过最小二乘法正弦曲线拟合获得转速频率分量;
步骤9:生成转速方波信号:使用所述步骤8中获取的正弦信号的过零点时间信息作为转速信号的跳变点,生成转速方波信号,为后续的信号分析备用。
2.根据权利要求1所述一种使用振动信号精确提取转速信号的方法,其特征在于,所述步骤1中获得振动信号的具体方式为:假设设备在采样信号长度的时间区间内转速稳定,获取振动传感器的原始数字信号X(n),截取2^N个采样点,使获得的数字信号时间长度超过但接近1秒钟。
3.根据权利要求2所述一种使用振动信号精确提取转速信号的方法,其特征在于,对获得原始数字信号X(n)加载汉宁窗w(n)
窗函数的数据和原始数字信号的数据作对应点之间的乘法,获得加载窗函数后的新的信号数据点X′(n),
X′(n)=w (n)*X(n)。
4.根据权利要求1所述一种使用振动信号精确提取转速信号的方法,其特征在于,所述步骤3中所述快速傅立叶变换方式为:
5.根据权利要求1所述一种使用振动信号精确提取转速信号的方法,其特征在于,所述步骤6中具体步骤为将频域信号转化为极坐标形式,由幅值谱和相位谱组成,把幅值谱单列出来进行下一步处理,相位谱保持不变。
6.根据权利要求1所述一种使用振动信号精确提取转速信号的方法,其特征在于,所述步骤4中已知转速范围估算方法为:通过设备型号获得恒定转速R,在幅值谱上查询到其所对应的峰值;如果转速处于变动状态当中,则根据前次循环处理的结果,在正负5%的频率范围内查询幅值谱的就近峰值。
7.根据权利要求1所述一种使用振动信号精确提取转速信号的方法,其特征在于,所述步骤8中获取的正弦信号的周期的倒数和步骤5中获得的复核转速频率相比较,获取差值,所述差值在1RPM之内则校验为合格,否则返回到步骤1,重置整个过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610322710.5A CN106092524B (zh) | 2016-05-13 | 2016-05-13 | 一种使用振动信号精确提取转速信号的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610322710.5A CN106092524B (zh) | 2016-05-13 | 2016-05-13 | 一种使用振动信号精确提取转速信号的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106092524A CN106092524A (zh) | 2016-11-09 |
CN106092524B true CN106092524B (zh) | 2018-11-13 |
Family
ID=57229921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610322710.5A Active CN106092524B (zh) | 2016-05-13 | 2016-05-13 | 一种使用振动信号精确提取转速信号的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106092524B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106769033B (zh) * | 2016-11-30 | 2019-03-26 | 西安交通大学 | 基于阶次包络时频能量谱的变转速滚动轴承故障识别方法 |
CN106872721B (zh) * | 2017-02-14 | 2019-04-30 | 浙江聚励云机械科技有限公司 | 一种基于电瓶正极电压的发动机转速计算方法 |
CN107798298B (zh) * | 2017-09-30 | 2020-06-05 | 安徽容知日新科技股份有限公司 | 一种信号处理方法、装置和计算设备 |
CN107918032B (zh) * | 2017-11-14 | 2020-04-14 | 湖南大学 | 一种空间多源声音信号融合的转速测量方法 |
CN110017957B (zh) * | 2018-01-10 | 2021-08-31 | 国家能源投资集团有限责任公司 | 一种同步分析方法、装置及*** |
CN108709635B (zh) * | 2018-08-17 | 2020-07-28 | 国家电网有限公司 | 确定旋转机械振动信号主要频率成分信息的方法及装置 |
JP7352371B2 (ja) * | 2019-04-04 | 2023-09-28 | 株式会社日立製作所 | 診断装置 |
CN110470465B (zh) * | 2019-09-19 | 2021-09-10 | 武汉市华英电力科技有限公司 | 一种基于振动信号分析的断路器测试方法及*** |
CN110987438B (zh) * | 2019-12-04 | 2021-12-28 | 国网福建省电力有限公司 | 水轮发电机变转速过程周期性振动冲击信号检测的方法 |
CN113188644A (zh) * | 2021-03-31 | 2021-07-30 | 大唐东营发电有限公司 | 一种旋转机械振动测试***谐波振动幅值计算方法 |
CN114264839A (zh) * | 2021-12-30 | 2022-04-01 | 燕山大学 | 一种应用加速度计的旋转体转速测定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095885A (zh) * | 2010-10-28 | 2011-06-15 | 华南理工大学 | 应用四点能量重心校正法的发动机转速的振动测量方法 |
CN102590543A (zh) * | 2012-02-24 | 2012-07-18 | 中国汽车技术研究中心 | 一种转速测量装置及其测量方法 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN103743922A (zh) * | 2014-01-26 | 2014-04-23 | 深圳科瑞技术股份有限公司 | 一种基于dsp的发动机转速测量方法和装置 |
CN104374939A (zh) * | 2014-11-06 | 2015-02-25 | 西安交通大学 | 基于振动信号同步压缩变换的旋转机械瞬时转速估测方法 |
CN104865400A (zh) * | 2015-04-14 | 2015-08-26 | 华北电力大学 | 一种风电机组转速的检测识别方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5293300B2 (ja) * | 2009-03-16 | 2013-09-18 | 富士電機株式会社 | 回転機の振動監視装置および振動監視方法 |
-
2016
- 2016-05-13 CN CN201610322710.5A patent/CN106092524B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095885A (zh) * | 2010-10-28 | 2011-06-15 | 华南理工大学 | 应用四点能量重心校正法的发动机转速的振动测量方法 |
CN102590543A (zh) * | 2012-02-24 | 2012-07-18 | 中国汽车技术研究中心 | 一种转速测量装置及其测量方法 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN103743922A (zh) * | 2014-01-26 | 2014-04-23 | 深圳科瑞技术股份有限公司 | 一种基于dsp的发动机转速测量方法和装置 |
CN104374939A (zh) * | 2014-11-06 | 2015-02-25 | 西安交通大学 | 基于振动信号同步压缩变换的旋转机械瞬时转速估测方法 |
CN104865400A (zh) * | 2015-04-14 | 2015-08-26 | 华北电力大学 | 一种风电机组转速的检测识别方法及*** |
Non-Patent Citations (1)
Title |
---|
从实测振动信号中精确计算叶轮机械主轴转速的方法;雷继祖等;《西南民族学院学报(自然科学版)》;20020831;第28卷(第3期);第348-350页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106092524A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106092524B (zh) | 一种使用振动信号精确提取转速信号的方法 | |
CN110763462B (zh) | 一种基于同步压缩算子的时变振动信号故障诊断方法 | |
Lamraoui et al. | Cyclostationarity approach for monitoring chatter and tool wear in high speed milling | |
EP2199764B1 (en) | Timing analysis | |
CN106769033A (zh) | 基于阶次包络时频能量谱的变转速滚动轴承故障识别方法 | |
Meng et al. | General synchroextracting chirplet transform: Application to the rotor rub-impact fault diagnosis | |
US11927501B2 (en) | Method and device for monitoring a gear system | |
CN109116048A (zh) | 一种汽车发动机转速检测及控制方法 | |
Wang et al. | Size estimation for naturally occurring bearing faults using synchronous averaging of vibration signals | |
Corne et al. | Comparing MCSA with vibration analysis in order to detect bearing faults—A case study | |
Wang et al. | A two-stage method using spline-kernelled chirplet transform and angle synchronous averaging to detect faults at variable speed | |
Wu et al. | A modified tacho-less order tracking method for the surveillance and diagnosis of machine under sharp speed variation | |
JP6354349B2 (ja) | 振動検出装置と工作機械 | |
CN103221791B (zh) | 用于监测风力涡轮机的机械传动装置中的旋转部件的方法 | |
JP2012047707A (ja) | 振動検出装置、振動抑制装置、および、振動情報表示装置 | |
CN106353528A (zh) | 电机转速测量装置和方法 | |
KR101372406B1 (ko) | 회전축 분석 장치 | |
CN112781723A (zh) | 一种基于频谱方差的谐波成分检测方法 | |
CN104236798B (zh) | 旋转机械启停车纯试重奈奎斯特图的单面图解平衡方法 | |
JP7161967B2 (ja) | 回転解析装置 | |
JP2021505872A (ja) | 回転数測定のための勾配決定 | |
Kidar et al. | Control of phases by ESPRIT and WLSE methods for the early detection of gear cracks | |
CN112268688B (zh) | 一种错误数据识别方法及装置 | |
Zhou et al. | A Fast Tacho-Less Order Tracking Method for Gear Fault Diagnosis Under Large Rotational Speed Variation Conditions Based on Multi-stage Generalized Demodulation | |
Zhao et al. | Blind source separation based on EMD and correlation measure for rotating machinery fault diagnosis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |