CN106768260A - 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 - Google Patents
可抑制直流扰动的振动信号最大功率频率分量实时检测算法 Download PDFInfo
- Publication number
- CN106768260A CN106768260A CN201611151180.9A CN201611151180A CN106768260A CN 106768260 A CN106768260 A CN 106768260A CN 201611151180 A CN201611151180 A CN 201611151180A CN 106768260 A CN106768260 A CN 106768260A
- Authority
- CN
- China
- Prior art keywords
- vibration signal
- maximum power
- frequency component
- power frequency
- direct current
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
- G01H1/14—Frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
- G01H11/06—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法。首先,给出了迭代处理振动信号获取最大功率频率分量的四阶系统状态方程、自适应律及四阶系统阻尼比的自动计算公式;其次,将振动信号乘一增益系数后,迭代执行四阶系统状态方程、自适应律及阻尼比自动计算公式,可动态跟踪振动信号角频率。最后,在频率跟踪进入稳态后,根据角频率收敛值与状态变量计算振动信号频率、幅值与相位。算法实时性强,可避免直流扰动对检测结果的影响,有效提取振动信号中最大功率频率分量信号。
Description
技术领域
本发明涉及振动控制与振动测试技术领域,尤其涉及一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法。
背景技术
随着计算机技术及控制理论的发展,利用控制器产生反向力来抑制机械振动的反馈补偿控制方法在精密仪器、机加工及化工过程机械等领域受到广泛重视。其中,振动信号最大功率对应频率分量的频率、幅值、相位检测精度及检测算法的实时性能是制约控制系统带宽、稳定性的关键因素。本申请人曾在专利申请201510375357.2中公开过一种振动信号幅值、相位实时计算方法,但该方法需预先测量信号中的直流分量,并在软件做偏移处理将振动信号转换成双极性归一化信号。如果直流分量处理不当,将导致频率、幅值、信号估计结果不正确。
发明内容
本发明目的是提出一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法,克服已有技术的缺点,避免直流分量的影响,快速准确地检测出具有功率谱最大分量振动信号的频率、幅值和相位。
一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于具体步骤包括:
采用四阶系统状态方程及自适应律迭代处理振动信号获取最大功率频率分量,其中,状态方程为:
设为振动信号,为增益,式中,为四阶系统输入信号,并且有,为四阶系统状态变量,为阻尼比,为振动信号最大功率频率分量对应的角频率估计值;与之相应的自适应律为,式中,为自适应因子;对于正弦输入信号,迭代执行状态方程及自适应律,角频率估计值收敛于,并且有,四阶系统微分方程的周期解为:
可见,状态变量,准确跟踪上了输入信号。当输入信号包含多种频率成分时,,跟踪的是输入信号中的最大功率频率分量。
预处理。设置、,通带带宽,角频率估计初值,采样频率;计算极半径,计算采样时间。对ADC采样振动信号u 做乘增益k 处理以加快频率跟踪的动态过程。预处理参数设置规则为:假设为最大功率频率分量,预设,时,需保证;通常采样频率越高相位估计越准确,并且有,。
迭代执行算法
①计算阻尼比,,,;
②根据输入及该瞬时的系统状态,计算四阶系统下一时刻的状态,算法如下:
③根据当前时刻的频率估计值、状态变量及自适应律公式计算下一个时刻的频率估计值,算法如下:
④迭代计算①、②、③。
根据角频率估计进入稳态后的收敛值及状态变量,得到振动信号最大功率频率分量的频率、幅值及相位,算法如下:
。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
本发明通过将拾取到的转子振动信号乘一增益后送入一四阶系统,并通过角频率自适应律动态跟踪振动信号角频率。在频率跟踪进入稳态后,根据角频率收敛值及状态变量的稳态值计算振动信号相位、幅值。本发明的主要优点有:
可消除直流偏置对振动信号频率、幅值、相位计算值的影响,无需预先测量振动信号中的直流量,再将振动信号转换成双极性归一化信号。
只需要设置自适应因子,角频率预估值及通带带宽三个初值,并给出了阻尼比的计算算法。
振动最大功率频率分量相位、幅值的计算值不受阻尼比影响。
无需整周期截取、整数倍采样、抗混叠滤波处理输入信号。
算法实时性强,可内嵌至电机驱动器内实现振动控制或振动测试。
附图说明
图1是本发明所述的算法原理框图。
图2是带直流偏置的正弦波曲线。
图3是直流偏置正弦波频率跟踪曲线。
图4是原始机械振动信号时域曲线。
图5是最大功率频率分量振动信号频率估计曲线。
图6是12~13s期间振动信号的FFT幅度谱。
图7是12~12.09s期间振动信号的FFT幅度谱。
图8是12~12.09s期间的四阶系统状态变量曲线。
图9是最大功率频率分量振动信号相位曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合实施例及附图对本发明作进一步详细的描述。根据本发明的技术方案,在不变更本发明的实质精神下,本领域一般技术人员可以提出本发明的多种相关设计。因此以下具体实施方式以及附图仅是本发明技术方案的具体说明,而不应当视为本发明的全部或者视为对本发明技术方案的限定或限制。
参照图1,一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于具体步骤包括:
采用四阶系统状态方程及自适应律迭代处理振动信号获取最大功率频率分量,其中,状态方程为:
设为振动信号,为增益,式中,为四阶系统输入信号,并且有,为四阶系统状态变量,为阻尼比,为振动信号最大功率频率分量对应的角频率估计值;与之相应的自适应律为,式中,为自适应因子;对于正弦输入信号,迭代执行状态方程及自适应律,角频率估计值收敛于,并且有,四阶系统微分方程的周期解为:
可见,状态变量,准确跟踪上了输入信号。当输入信号包含多种频率成分时,,跟踪的是输入信号中的最大功率频率分量。
预处理。设置、,通带带宽,角频率估计初值,采样频率;计算极半径,计算采样时间。对ADC采样振动信号u 做乘增益k 处理以加快频率跟踪的动态过程。预处理参数设置规则为:假设为最大功率频率分量,预设,时,需保证;通常采样频率越高相位估计越准确,并且有,。
迭代执行算法
①计算阻尼比,,,;
②根据输入及该瞬时的系统状态,计算四阶系统下一时刻的状态,算法如下:
③根据当前时刻的频率估计值、状态变量及自适应律公式计算下一个时刻的频率估计值,算法如下:
④迭代计算①、②、③。
根据角频率估计进入稳态后的收敛值及状态变量,得到振动信号最大功率频率分量的频率、幅值及相位,算法如下:
。
下面采用具体的实例验证算法性能:
例1:跟踪频率为150Hz的纯正弦波,设输入信号为,采样频率,,角频率预估初值 rad/s,通带带宽Hz,自适应因子。可计算出,发明算法迭代执行步骤的离散化形式为:
。
图2为带直流偏置的原始正弦波曲线,对正弦波进行定标处理,正弦波幅值为1PU,直流偏置为3PU。图3为增益分别为1、5时的频率跟踪曲线,可见算法很好抑制了直流干扰。同时,增加可加快频率跟踪速度,当时,频率跟踪可在0.8s内进入稳态,频率估计最大相对误差为0.25%。
例2:实验装置采用单级渐开线圆柱齿轮减速器,减速比2.24,减速箱所用齿轮模数均为1.5,齿数分别为91、41,故障轴承座滚动轴承型号为6313ZZ(内含8个人为磨损滚珠)。通过变频器设定三相异步电机转速为1200rpm,采用AC192振动传感器拾取轴承座径向机械振动信号,AC192输出信号经过调理电路调整电平范围,再经过ADS1205实现Δ-Σ过采样,AMC1210实现sinc滤波, TMS320LF2808以10KHz的采样频率读取AMC1210的寄存器获取振动信号。设置角频率预估初值rad/s,通带带宽Hz,自适应因子。发明算法迭代执行步骤的离散化形式与例1相同。
图4为带直流偏置的原始机械振动时域曲线图,图5为最大功率频率分量振动信号频率估计曲线,进入稳态后,估计频率在152Hz上下波动。图6为12~13s期间振动信号的FFT幅度谱,幅值谱峰值对应的频率分量为152.58Hz。图7为12~12.09s期间振动信号的FFT幅度谱,幅值谱峰值对应的频率分量为152.39Hz。图5、图6、图7说明算法可动态跟踪转速变化,频率估计最大相对误差0.31%。图8为12~12.09s期间的四阶系统状态变量曲线,其中为最大功率频率分量振动信号的时域曲线。图9为最大功率频率分量振动信号相位曲线图,可验证所估计的相位跟信号频率一致,但相位估计精度受信号频率及采样频率的影响。在输入信号频率不变的前提下,采样频率越高,相位估算精度越高。
Claims (5)
1.一种可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于,采用四阶系统状态方程及自适应律迭代处理振动信号获取最大功率频率分量,其中,状态方程为:
设为振动信号,为增益,式中,为四阶系统输入信号,并且有,为四阶系统状态变量,为阻尼比,为振动信号最大功率频率分量对应的角频率估计值;
与之相应的自适应律为:,式中,为自适应因子。
2.根据权利要求1所述的可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于,对输入信号做预处理,预处理方法为:设置、,通带带宽,角频率估计初值,采样频率;计算极半径,计算采样时间。
3.根据权利要求1、2所述的可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于,预处理参数设置:假设为最大功率频率分量,预设,时,需保证;通常采样频率越高相位估计越准确,并且有,。
4.根据权利要求1所述的可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于迭代执行算法,具体步骤为:
①计算阻尼比,,,;
②根据输入及该瞬时的系统状态,计算四阶系统下一时刻的状态;
③根据当前时刻的频率估计值、状态变量及自适应律公式计算下一个时刻的频率估计值;
④迭代计算①、②、③。
5.根据权利要求1~4所述的可抑制直流扰动的振动信号最大功率频率分量实时检测算法,其特征在于,根据角频率估计进入稳态后的收敛值及状态变量,得到振动信号最大功率频率分量的频率、幅值及相位,算法如下:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611151180.9A CN106768260B (zh) | 2016-12-14 | 2016-12-14 | 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611151180.9A CN106768260B (zh) | 2016-12-14 | 2016-12-14 | 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106768260A true CN106768260A (zh) | 2017-05-31 |
CN106768260B CN106768260B (zh) | 2019-05-31 |
Family
ID=58887532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611151180.9A Expired - Fee Related CN106768260B (zh) | 2016-12-14 | 2016-12-14 | 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106768260B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107733313A (zh) * | 2017-12-04 | 2018-02-23 | 广东石油化工学院 | 旋转机械谐振自适应抑制方法 |
CN108037351A (zh) * | 2017-12-04 | 2018-05-15 | 广东石油化工学院 | 一种变频电源功率因素测量方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101451898A (zh) * | 2009-01-13 | 2009-06-10 | 华北电力大学 | 汽轮发电机组转子热弯曲不平衡故障实时诊断方法 |
CN101576752A (zh) * | 2009-06-15 | 2009-11-11 | 南京航空航天大学 | 柔性结构主动减振装置及其控制方法 |
US20110036166A1 (en) * | 2009-08-14 | 2011-02-17 | Lenz Michael A W | Method and apparatus for in situ unbalance and corrective balance determination for a non-vertical axis rotating assembly |
DE102011012477A1 (de) * | 2011-02-25 | 2012-08-30 | Schenck Rotec Gmbh | Verfahren zur Bestimmung von Größe und Lage einer Unwucht |
CN102967760A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 指定频率的信号分析方法 |
CN104111109A (zh) * | 2014-07-21 | 2014-10-22 | 石家庄铁道大学 | 一种基于不同阶次统计量及支持向量机的机械振动状态识别方法 |
CN105005548A (zh) * | 2015-06-01 | 2015-10-28 | 莆田学院 | 基于arma模型获取四阶振动系统动态特性参数的方法及装置 |
CN105043667A (zh) * | 2015-07-01 | 2015-11-11 | 广东石油化工学院 | 转子不平衡振动信号幅值、相位实时计算方法 |
-
2016
- 2016-12-14 CN CN201611151180.9A patent/CN106768260B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101451898A (zh) * | 2009-01-13 | 2009-06-10 | 华北电力大学 | 汽轮发电机组转子热弯曲不平衡故障实时诊断方法 |
CN101576752A (zh) * | 2009-06-15 | 2009-11-11 | 南京航空航天大学 | 柔性结构主动减振装置及其控制方法 |
US20110036166A1 (en) * | 2009-08-14 | 2011-02-17 | Lenz Michael A W | Method and apparatus for in situ unbalance and corrective balance determination for a non-vertical axis rotating assembly |
DE102011012477A1 (de) * | 2011-02-25 | 2012-08-30 | Schenck Rotec Gmbh | Verfahren zur Bestimmung von Größe und Lage einer Unwucht |
CN102967760A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 指定频率的信号分析方法 |
CN104111109A (zh) * | 2014-07-21 | 2014-10-22 | 石家庄铁道大学 | 一种基于不同阶次统计量及支持向量机的机械振动状态识别方法 |
CN105005548A (zh) * | 2015-06-01 | 2015-10-28 | 莆田学院 | 基于arma模型获取四阶振动系统动态特性参数的方法及装置 |
CN105043667A (zh) * | 2015-07-01 | 2015-11-11 | 广东石油化工学院 | 转子不平衡振动信号幅值、相位实时计算方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107733313A (zh) * | 2017-12-04 | 2018-02-23 | 广东石油化工学院 | 旋转机械谐振自适应抑制方法 |
CN108037351A (zh) * | 2017-12-04 | 2018-05-15 | 广东石油化工学院 | 一种变频电源功率因素测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106768260B (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xue et al. | An adaptively fast ensemble empirical mode decomposition method and its applications to rolling element bearing fault diagnosis | |
Cao et al. | Vibration signal correction of unbalanced rotor due to angular speed fluctuation | |
CN105247432B (zh) | 频率响应测定装置 | |
CN106546400B (zh) | 一种多输入多输出非高斯随机振动试验系统及试验方法 | |
CN103038620A (zh) | 滚动轴承的异常诊断装置 | |
Wu et al. | Order-domain-based harmonic injection method for multiple speed harmonics suppression of PMSM | |
CN107576851B (zh) | 基于旋转电流相量的系统谐波阻抗测量方法 | |
CN107490722A (zh) | 一种低信噪比实信号的频率估计方法 | |
CN101832849B (zh) | 基于三参量控制的振动台软启动控制方法 | |
CN106768260A (zh) | 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 | |
CN115575707A (zh) | 基于改进fft算法与小波变换结合的谐波检测装置及方法 | |
CN105137174A (zh) | 应用于apf的变步长lms自适应谐波检测方法 | |
CN105043667B (zh) | 转子不平衡振动信号幅值、相位实时计算方法 | |
CN112255457B (zh) | 适用于自动准同期装置的相角差测量方法 | |
CN107733313A (zh) | 旋转机械谐振自适应抑制方法 | |
CN112346381A (zh) | 一种基于fpga的数字相敏检波方法及系统 | |
CN112595893B (zh) | 一种基于自适应滤波器的电网相位同步信号检测方法 | |
CN111353415A (zh) | 一种脉冲响应中谐波成分的检测方法 | |
CN106324342A (zh) | 一种基于查表的谐波检测方法 | |
Nguyen et al. | Auto-calibration and noise reduction for the sinusoidal signals of magnetic encoders | |
CN110263482A (zh) | 一种基于互相关算法的涡流阻抗求解方法及装置 | |
Yang et al. | A novel phase locked loop for grid-connected converters under non-ideal grid conditions | |
CN113804949A (zh) | 一种适用于宽变频交流发电系统的幅值确定方法及装置 | |
CN108631679A (zh) | 一种位置信号滤波方法及系统 | |
CN107576842B (zh) | 一种宽频带同步采样方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190531 Termination date: 20191214 |