CN1865886A - 一种用于振动时效的频谱分析判峰方法 - Google Patents

一种用于振动时效的频谱分析判峰方法 Download PDF

Info

Publication number
CN1865886A
CN1865886A CN 200610009897 CN200610009897A CN1865886A CN 1865886 A CN1865886 A CN 1865886A CN 200610009897 CN200610009897 CN 200610009897 CN 200610009897 A CN200610009897 A CN 200610009897A CN 1865886 A CN1865886 A CN 1865886A
Authority
CN
China
Prior art keywords
frequency
acceleration
workpiece
peak
oscillating aging
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.)
Pending
Application number
CN 200610009897
Other languages
English (en)
Inventor
王鸿鹏
李会滨
韩殿文
董德权
Original Assignee
韩殿文
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 韩殿文 filed Critical 韩殿文
Priority to CN 200610009897 priority Critical patent/CN1865886A/zh
Publication of CN1865886A publication Critical patent/CN1865886A/zh
Pending legal-status Critical Current

Links

Images

Abstract

本发明公开了一种用于振动时效的频谱分析判峰方法。该方法的主要步骤为:数据采集设备对工件受迫振动加速度值进行采样,对采样后的加速度值进行量化、低通滤波、加窗处理、离散傅立叶变换处理之后找到共振频率及其它各阶谐振频率。该方法主要优点是能快速检测出工件的共振频率和谐振频率,提高了小型工件的共振频率和谐振频率检测的准确度,并且能够减少时效处理时间、降低噪声、节约电能、减小振动时效设备的损失。

Description

一种用于振动时效的频谱分析判峰方法
技术领域
本发明涉及对机械工件进行振动时效过程中对采集到的数据进行频谱分析以及根据分析结果对工件进行振动时效的技术领域。
背景技术
在机械加工行业,振动时效工艺是通过振动时效设备的控制系统控制联结于工件上的激振器的转动使工件被迫振动,工件受迫振动的动应力与工件内残余应力叠加的和超过某一极限,使工件发生微观塑性变形,来降低或均化工件内部的残余应力,使工件尺寸趋于稳定,达到要求精度。
通常振动时效设备由以下主要部件组成:控制器、激振器、加速度(振幅)传感器、电缆等,如附图1所示。
判峰方法是振动时效设备寻找工件固有频率(共振频率)及谐振频率的方法。常规振动时效设备的判峰方法是用控制器控制激振器转速从激振器的最低转速缓慢升到最高转速(这一过程称为扫频),同时利用加速度传感器采集工件受迫振动的加速度值(振幅值),通过分析加速度-转速曲线的极值点来判断工件的共振峰和谐振峰(工件的固有频率和谐振频率对应的加速度-转速曲线上的极值点分别称为共振峰和谐振峰),这是传统的判峰方法,是在时域进行的。共振峰与谐振峰在加速度-转速曲线中的位置如附图2所示。使用这种时域判峰方法的时效设备一般具有以下缺点:
一、常规振动时效设备扫频时一般从激振器的最低转速缓慢升到最高转速,进行扫频耗费时间较长(一般需要3-6分钟,如果再进行二次扫频以判断时效的效果时,时间加倍),这样在进行大批量工件的振动时效时,就会严重影响工作效率,并且振动时效设备工作时间较长,会导致电力资源的浪费,振动时效设备本身的使用寿命缩短等问题;另外,操作人员也要忍受长时间的噪声。
二、对于大型工件,采用常规的判峰方法的振动时效设备通常可以分析出工件的共振频率和谐振频率。而对于小型工件,情况则不然。一方面,由于小型工件的刚性一般都比较大,扫频的时候,在不同的激振器转速下工件的受迫振动加速度值差异不大,导致用普通的判峰方法很难找到可以进行时效的共振峰或者谐振峰;另一方面,联接在小型工件上的激振器在扫频时容易导致工件整体随着激振器进行振动,这样加速度传感器采集回来的加速度值就不准确,有可能误判共振频率,致使振动时效不能达到预期的效果。小型工件加速度-转速曲线如附图3所示,共振峰与谐振峰差异不大,很难发现共振峰。
三、对于固有频率很高的一些工件,在激振器的额定转速范围内,用常规的扫频方式,不能得到可以进行时效的共振峰或者谐振峰。
发明内容
本发明的目的是提供一种能够减少扫频时间、对小型工件判峰准确的用于振动时效的频谱分析判峰方法。
本发明的目的是通过以下技术方案实现的:被处理工件定位后,振动时效设备中的激振器和加速度传感器分别被安置于工件上,激振器和加速度传感器分别通过导线与振动时效设备的控制器相连接,通过加速度传感器对加速度信号进行采样、量化得到加速度值,通过控制器对经过采样和量化处理的加速度信号进行低通滤波、加窗处理,使用离散傅立叶变换DFT对数据进行分析,得到各个频率下的加速度值,从中选出加速度最大的值所对应的频率即为工件的共振频率,其余振幅较大的频率为谐振频率。如附图4所示。
具体步骤是:
第一步:通过数据采集设备对加速度信号进行采样和量化。设采样频率为FS,采样频率是根据采样定理来设定的,采样频率要大于两倍的激振器最高频率,对采样数据采用若干比特进行量化。通常情况下,振动时效设备的激振器带动工件振动,工件的振动频率与转速成正比。由于激振器本身具有一个最高转速限制,即具有一个最高振动频率,那么我们选择激振器能够产生的最高振动频率作为信号的最高频率,信号中高于激振器能够产生的最高振动频率的部分将在后续的低通滤波过程中被滤除。量化过程是对经过采样的加速度信号用设定数量的比特进行编码,通常采用8位或者16位比特。
第二步:对经过采样和量化处理的加速度信号进行低通滤波处理。滤除频率高于激振器最高振动频率的采样信号。低通滤波器有多种设计方法,可根据实际情况选择。
第三步:对通过低通滤波的加速度采样信号进行加窗处理,加窗处理对加速度采样信号进行了截断,确定了采用多少个采样点参与后续处理。设这个点数为N,同时加窗过程中不可避免要产生频谱泄露,所以加窗时可以根据实际情况采用不同的窗函数来抑制频谱泄露。
例如:可以采用海明窗,它的公式如下:
w ( n ) = [ 0.54 - 0.46 cos ( 2 πn N - 1 ) ] R N ( n )
RN(n)为矩形窗,即在加速度采样信号中截取N个采样信号。
第四步:对经过加窗处理的采样信号进行离散傅立叶变换,也可通过离散傅立叶变换(DFT)的其它改进方法,例如:快速傅立叶变换FFT来完成这个转换。
离散傅立叶变换将时域信号转化为频域信号,它的公式如下:
X [ k ] = Σ n = 0 N - 1 x ( n ) e - j 2 π N nk = Σ n = 0 N - 1 x ( n ) W N nk , 0 ≤ k ≤ N - 1
其中:X[k]为经过离散傅立叶变换处理过的频域值;
x(n)为加速度采样信号;
W N nk = e - j 2 π N nk
N为经过加窗处理后截取的加速度采样点数量;
第五步:加速度采样信号经过DFT或者FFT处理过以后,产生加速度-数字频率曲线,设加速度-数字频率曲线中的数字频率为w,频率为F,数字频率与频率之间存在着换算关系: w = 2 π F F S
FS为采样频率,通过对加速度-数字频率曲线上加速度值的比较,找到加速度值最大点对应的数字频率,根据这个数字频率计算得到工件的固有频率。然后根据需要的谐振峰数量找出等数量的次高点,根据次高点对应数字频率利用上述公式计算出次高点对应的频率,即得到需要的工件的谐振峰频率。由于频率与转速成正比,所以可以找到工件共振峰和谐振峰对应的转速。
本发明能够产生的有益效果:
一、减少了扫频时间,节省了能源消耗,减少了振动时效设备本身的损耗。通过激振器快速升速低速振动或者采用敲击方式使工件产生受迫振动,节省了扫频时间,在很短时间内达到了采集加速度数据,判断工件共振频率和谐振频率的目的,达到了省时、省电、减少设备损耗、降低时效时的噪音的目的。
二、解决了采用传统判峰方法的振动时效设备无法找出部分小工件和刚性很大的工件的谐振峰的问题。而使用本发明中的基于频谱分析判峰方法的振动时效设备,可以计算出工件的共振频率和各阶谐振频率,控制器中的软件可以在激振器的额定转速范围内选择几个谐振峰进行时效。
三、提高对小型工件判峰的准确性。频谱分析判峰方法通过对信号的各个频率上的加速度进行分析,能够清晰分辨出工件的共振振动频率,减少了时域判峰方法对小型工件共振频率误判的问题。
四、对于共振频率超过激振器振动频率范围的工件,可以找到工件在激振器振动频率内的谐振峰,解决了共振频率超过激振器振动频率范围工件的振动时效问题。
附图说明
图1是振动时效设备组成示意图;
图2是工件的加速度-转速曲线示意图;
图3是小型工件的加速度-转速曲线示意图;
图4是本发明频谱分析判峰方法示意图。
具体实施方式
如图1所示;1.1控制器、1.2工件、1.3激振器、1.4激振器连接电缆、1.5加速度传感器、1.6加速度传感器连接导线。激振器和加速度传感器分别通过导线与工件连接,激振器和加速度传感器分别通过导线与振动时效设备的控制器相连接。
如图2所示,2.1为共振峰,2.2、2.3、2.4为谐振峰。
将激振器和加速度传感器安装在工件上,启动具备频谱分析判峰方法的振动时效设备,在频谱分析的功能菜单下,显示(一)激振式(二)敲击式,用户可以选择其中之一以便进行基于频谱分析判峰方法的振动时效。如果选择激振式,激振器在低速区快速扫频,同时数据采集设备对加速度信号进行采样、量化得到加速度值。采集过程结束时,控制器中的软件系统中的频谱分析判峰方法对加速度信号低通滤波、加窗处理、进行频谱分析,得到最大加速度对应的频率,即得到工件的共振频率,同时得到各阶谐振频率,然后控制器控制激振器在共振频率或者几个选定的谐振频率对应的转速附近进行振动时效,如果工件的共振频率对应的转速超过激振器的转速范围,则选择工件在激振器转速范围内的谐振峰附近进行时效,同时通过数据采集设备采集加速度信号值,当满足预定的判断标准时,在相应的振动频率处的时效结束。如果选择敲击式,则立即用锤子在工件上敲击一下,这时数据采集设备采集工件的加速度信号,以下处理过程同激振式相同。

Claims (6)

1、一种用于振动时效的频谱分析判峰方法,被处理工件定位后,振动时效设备中的激振器和加速度传感器分别被安置于工件上,激振器和加速度传感器分别通过导线与振动时效设备的控制器相连接,其特征是:通过加速度传感器对加速度信号进行采样、量化得到加速度值,通过控制器对经过采样和量化处理的加速度信号进行低通滤波、加窗处理,使用离散傅立叶变换DFT对数据进行分析,得到各个频率下的加速度值,从中选出加速度最大的值所对应的频率即为工件的共振频率,其余振幅较大的频率为谐振频率。
2、根据权利要求1所述的一种用于振动时效的频谱分析判峰方法,其特征是:所述通过加速度传感器对加速度信号进行采样和量化,采样频率要大于两倍的激振器最高频率,对采样数据采用若干比特进行量化。
3、根据权利要求1所述的一种用于振动时效的频谱分析判峰方法,其特征是:所述对经过采样和量化处理的加速度信号进行低通滤波处理时,滤除频率高于激振器最高振动频率的采样信号。
4、根据权利要求1所述的一种用于振动时效的频谱分析判峰方法,其特征是:所述对通过低通滤波的加速度采样信号进行加窗处理,加窗处理对加速度采样信号进行了截断,确定了采用多少个采样点参与后续处理。
5、根据权利要求1所述的一种用于振动时效的频谱分析判峰方法,其特征是:对经过加窗处理的采样信号进行离散傅立叶变换,也可通过快速傅立叶变换FFT来完成这个转换。
6、根据权利要求1所述的一种用于振动时效的频谱分析判峰方法,其特征是:所述的使用离散傅立叶变换DFT对数据进行分析,是指加速度采样信号经过DFT或者FFT处理过以后,产生加速度-数字频率曲线,通过对加速度-数字频率曲线上加速度值的比较,找到加速度值最大点对应的数字频率,根据这个数字频率计算得到工件的固有频率,然后根据需要的谐振峰数量找出等数量的次高点,根据次高点对应数字频率计算出次高点对应的频率,即得到需要的工件的谐振峰频率。
CN 200610009897 2006-04-01 2006-04-01 一种用于振动时效的频谱分析判峰方法 Pending CN1865886A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200610009897 CN1865886A (zh) 2006-04-01 2006-04-01 一种用于振动时效的频谱分析判峰方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200610009897 CN1865886A (zh) 2006-04-01 2006-04-01 一种用于振动时效的频谱分析判峰方法

Publications (1)

Publication Number Publication Date
CN1865886A true CN1865886A (zh) 2006-11-22

Family

ID=37424969

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200610009897 Pending CN1865886A (zh) 2006-04-01 2006-04-01 一种用于振动时效的频谱分析判峰方法

Country Status (1)

Country Link
CN (1) CN1865886A (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101303575B (zh) * 2007-05-10 2010-05-19 肖维鑫 具有十五个操作键的智能型振动时效控制器
CN101672692B (zh) * 2009-10-14 2011-02-09 大连理工大学 一种基于虚拟仪器的音叉谐振频率快速测量方法
CN102416419A (zh) * 2011-07-27 2012-04-18 中南大学 一种金属构件振动蠕变成形的方法和装置
CN101691632B (zh) * 2009-06-26 2012-11-14 哈尔滨工业大学深圳研究生院 快速判峰的方法
CN102816920A (zh) * 2011-06-10 2012-12-12 安徽省科捷再生能源利用有限公司 频谱谐波时效设备
CN101674520B (zh) * 2009-03-30 2013-01-16 东莞泉声电子有限公司 振动式拾音器参数分析测试方法和装置
CN101539113B (zh) * 2008-03-17 2013-04-24 西门子公司 确定风力涡轮机塔架的谐振频率的装置和方法
CN103335802A (zh) * 2013-07-01 2013-10-02 柳州五菱汽车有限责任公司 机动车组合开关共振点检测方法
CN103464536A (zh) * 2013-08-12 2013-12-25 西北工业大学 一种弹性变形条件下的振动时效成形方法及装置
CN104250686A (zh) * 2013-06-27 2014-12-31 李义根 小波频谱分析振动时效装置
CN105928677A (zh) * 2016-04-14 2016-09-07 北京工业大学 异形水箱结构液体晃动频率测试方法
CN105935857A (zh) * 2015-03-02 2016-09-14 通用汽车环球科技运作有限责任公司 用于减少裂纹倾向的机械糙化气缸孔的应力释放
CN106834657A (zh) * 2016-09-20 2017-06-13 上海海事大学 一种多维高频微观振动时效系统及方法
CN107144388A (zh) * 2017-05-17 2017-09-08 苏交科集团股份有限公司 一种柔性绳索振动频率的全域搜峰法
CN108020688A (zh) * 2017-10-30 2018-05-11 中国科学院上海微系统与信息技术研究所 自谱分时段傅里叶变换高量程加速度计共振频率提取方法
CN109632020A (zh) * 2019-01-30 2019-04-16 天津大学 能够排除共振干扰的涡街流量计测量方法

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101303575B (zh) * 2007-05-10 2010-05-19 肖维鑫 具有十五个操作键的智能型振动时效控制器
CN101539113B (zh) * 2008-03-17 2013-04-24 西门子公司 确定风力涡轮机塔架的谐振频率的装置和方法
CN101674520B (zh) * 2009-03-30 2013-01-16 东莞泉声电子有限公司 振动式拾音器参数分析测试方法和装置
CN101691632B (zh) * 2009-06-26 2012-11-14 哈尔滨工业大学深圳研究生院 快速判峰的方法
CN101672692B (zh) * 2009-10-14 2011-02-09 大连理工大学 一种基于虚拟仪器的音叉谐振频率快速测量方法
CN102816920A (zh) * 2011-06-10 2012-12-12 安徽省科捷再生能源利用有限公司 频谱谐波时效设备
CN102416419B (zh) * 2011-07-27 2014-04-23 中南大学 一种金属构件振动蠕变成形的方法
CN102416419A (zh) * 2011-07-27 2012-04-18 中南大学 一种金属构件振动蠕变成形的方法和装置
CN104250686A (zh) * 2013-06-27 2014-12-31 李义根 小波频谱分析振动时效装置
CN103335802A (zh) * 2013-07-01 2013-10-02 柳州五菱汽车有限责任公司 机动车组合开关共振点检测方法
CN103335802B (zh) * 2013-07-01 2016-03-23 广西汽车集团有限公司 机动车组合开关共振点检测方法
CN103464536A (zh) * 2013-08-12 2013-12-25 西北工业大学 一种弹性变形条件下的振动时效成形方法及装置
CN103464536B (zh) * 2013-08-12 2015-06-10 西北工业大学 一种弹性变形条件下的振动时效成形方法及装置
CN105935857A (zh) * 2015-03-02 2016-09-14 通用汽车环球科技运作有限责任公司 用于减少裂纹倾向的机械糙化气缸孔的应力释放
CN105928677A (zh) * 2016-04-14 2016-09-07 北京工业大学 异形水箱结构液体晃动频率测试方法
CN106834657A (zh) * 2016-09-20 2017-06-13 上海海事大学 一种多维高频微观振动时效系统及方法
CN106834657B (zh) * 2016-09-20 2020-02-07 上海海事大学 一种多维高频微观振动时效系统及方法
CN107144388A (zh) * 2017-05-17 2017-09-08 苏交科集团股份有限公司 一种柔性绳索振动频率的全域搜峰法
CN108020688A (zh) * 2017-10-30 2018-05-11 中国科学院上海微系统与信息技术研究所 自谱分时段傅里叶变换高量程加速度计共振频率提取方法
CN108020688B (zh) * 2017-10-30 2019-11-12 中国科学院上海微系统与信息技术研究所 自谱分时段傅里叶变换高量程加速度计共振频率提取方法
CN109632020A (zh) * 2019-01-30 2019-04-16 天津大学 能够排除共振干扰的涡街流量计测量方法

Similar Documents

Publication Publication Date Title
CN1865886A (zh) 一种用于振动时效的频谱分析判峰方法
CN110696051B (zh) 一种基于多分量信号分解的机械臂关节振动识别方法
CN1825082B (zh) 滚动轴承故障自动诊断系统
CN1884990A (zh) 大型发电机定子槽楔松动振动检测方法及装置
CN108844733B (zh) 一种基于kl散度与均方根值的齿轮状态监测指标提取方法
CN105058165A (zh) 基于振动信号的刀具磨损量监测系统
CN1844908A (zh) 一种超声波检测方法及其检测装置
CN105500115A (zh) 一种铣削加工中刀具颤振的检测系统及其检测方法
CN102252852A (zh) 柴油发动机信号采集装置及其燃烧质量检测和无负荷测功方法
CN1955710A (zh) 发动机燃油系统故障的检测装置及方法
CN109605128B (zh) 一种基于功率谱熵差的铣削颤振在线检测方法
CN1595188A (zh) 基于振动分析的真空断路器触头关合时刻的在线检测方法
CN109894925B (zh) 基于内嵌式压电传感器的薄壁件铣削加工振动监测方法
CN1873406A (zh) 声学材料性能综合测试仪
CN111175378A (zh) 一种螺栓拉力快速检测方法
CN1166934C (zh) 具有波纹板固定结构的旋转电机槽楔松动检测方法
CN2864705Y (zh) 利用扫频电源激振检测变压器绕组状态的装置
CN101691632B (zh) 快速判峰的方法
CN103630604B (zh) 离心式压缩机半开式叶轮裂纹故障识别方法
CN112539933A (zh) 一种基于共振基带宽傅立叶分解的齿轮箱故障诊断方法
CN110243400B (zh) 基于主动激励信号获取共振信号的触滑觉传感器
CN102680080A (zh) 一种基于改进的自适应形态滤波的非稳态信号检测方法
CN201438134U (zh) 一种快速判峰的振动时效系统
CN1629608A (zh) 设备状态数据采集器
CN2700889Y (zh) 实时电能质量综合分析仪

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication