CN109444263B - 一种基于频响函数的装配质量检测系统及方法 - Google Patents
一种基于频响函数的装配质量检测系统及方法 Download PDFInfo
- Publication number
- CN109444263B CN109444263B CN201811438317.8A CN201811438317A CN109444263B CN 109444263 B CN109444263 B CN 109444263B CN 201811438317 A CN201811438317 A CN 201811438317A CN 109444263 B CN109444263 B CN 109444263B
- Authority
- CN
- China
- Prior art keywords
- frequency
- unit
- signal
- output
- frequency response
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/045—Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/12—Analysing solids by measuring frequency or resonance of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/014—Resonance or resonant frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
Abstract
本发明公开了一种基于频响函数的装配质量检测系统及方法,包括:冲击激励单元、测量单元、分析单元以及输出显示单元;在被测衬板上设置至少一个测点,冲击激励单元在每一个测点位置产生振动信号;测量单元包括加速度传感器,在衬板上与冲击激励单元作用面相对的一面的每一个测点位置设置加速度传感器,采集加速度信号并传送至分析单元,所述分析单元对采集到的数据进行分析,判断装配质量是否合格;分析单元将分析结果传送至输出显示单元进行显示。本方案在检测衬板装配质量时不需要视觉系统,有效解决作业环境视觉遮挡、光线恶劣等问题;采用自动化方式实现检测分析,节约人力,减少工作人员工作量。
Description
技术领域
本发明涉及球磨机筒体内衬板的装配质量检测技术领域,具体而言,涉及到一种基于频响函数的装配质量检测系统及方法。
背景技术
球磨机筒体内侧安装着各种各样的衬板,用来保护磨机,使用一段时间后就会产生不同程度的磨损,需要对其不定期更换。更换的衬板的装配质量,影响衬板的使用寿命。在磨机衬板装配过程中,由于装配力或其他不确定性因素,会导致装配质量不稳定。目前装配质量的检测方法大部分基于机器视觉的方法,但在视觉遮挡等情况下,不能很好的进行装配质量判别。
发明内容
为了解决上述问题,本发明提出了一种基于频响函数(频率响应函数)的装配质量检测系统及方法,该系统及方法不需要视觉系统,有效解决作业环境视觉遮挡、光线恶劣等问题。
为了实现上述目的,本发明采用如下技术方案:
在一个或多个实施方式中公开的一种基于频响函数的装配质量检测系统,包括:冲击激励单元、测量单元、分析单元以及输出显示单元;
在被测衬板上设置至少一个测点,冲击激励单元在每一个测点位置产生振动信号;测量单元包括加速度传感器,在衬板上与冲击激励单元作用面相对的一面的每一个测点位置设置加速度传感器,采集加速度信号并传送至分析单元,所述分析单元对采集到的数据进行分析,判断装配质量是否合格;分析单元将分析结果传送至输出显示单元进行显示。
进一步地,所述冲击激励单元包括:用于在测点位置产生振动信号的力锤或者振动台。
进一步地,所述测量单元还包括:依次串联连接的放大电路和AD转换电路,放大电路与加速度传感器连接,AD转换电路与分析单元连接。
进一步地,所述分析单元对各测点进行频响函数估计,根据估计的频响函数,求出幅频特性和相频特性,并根据幅频特性和相频特性计算相干函数,通过相干函数评定频响函数估计的可信度;将幅频特性、相频特性及相干函数曲线输入装配质量库,采用模板匹配法与装配装量对应的特性曲线进行匹配,如果匹配成功则装配质量合格,否则不合格。
在一个或多个实施方式中公开的一种基于频响函数的装配质量检测方法,包括:
以反映被测衬板的构型为原则,在被测衬板上布置合适的测点;
采集设定采样频率,采用多点冲击激励法,采集每一个测点的加速度信号;
根据采集到的加速度信号计算出真正作用于被测衬板上的作用力;
对各测点进行频响函数估计;
根据估计的频响函数,求出幅频特性和相频特性;
根据幅频特性和相频特性计算相干函数,通过相干函数评定频响函数估计的可信度;
将幅频特性、相频特性及相干函数曲线输入装配质量库,采用模板匹配法与装配装量对应的特性曲线进行匹配,如果匹配成功则装配质量合格,否则不合格。
进一步地,所述根据采集到的加速度信号计算出真正作用于被测衬板上的作用力ft,具体为:
进一步地,所述对各测试的响应点进行频响函数估计,具体为:
设输入端激励力矩阵f(t)N×1,输出端响应矩阵x(t)O×1,对应频率域的激励力和响应力矩阵分别为F(f)N×1和X(f)O×1;
假设输入端的噪声矩阵m(t)N×1,对应测量单元输出端矩阵n(t)O×1,对应频率域输入噪声和输出噪声矩阵为M(f)N×1和N(f)O×1;
测量单元测得时域输入信号a(t)N×1,输出信号b(t)O×1,频域为A(f)N×1、B(f)O×1,则输入输出之间的关系为:
a(t)=f(t)+m(t)
b(t)=x(t)+n(t)
或者
A(f)=F(f)+M(f)
B(f)=X(f)+N(f)
则理想情况下,系统频响函数估计式为:X(f)=H(f)F(f);
进一步地,系统的频响函数估计模型根据H3估计法得到,即
其中,
H1(f)=Gfx(f)Gff(f)-1
H2(f)=Gxx(f)Gxf(f)-1
其中,Gfx(f)实测输入信号和实测输出信号的互谱均值,Gff(f)实测输入信号的的自谱均值,Gab(f)实测输入噪声信号的自谱均值,GaQ(f)系统输入噪声的自谱均值,Gxx(f)实测输出信号的自谱均值,Gxf(f)实测输出信号和实测输入信号的互谱均值,Gbb(f)实测输出噪声的自谱均值,GbQ(f)系统输出信号的自谱均值。
进一步地,根据估计的频响函数,求出幅频特性和相频特性,具体为:
θ(w)=∠G(jw);
其中,Φ(jw)为幅频特性、G(jw)为频率响应、X是信号的输入、Y是信号的输出、θ(w)信号的相频特性。
进一步地,根据幅频特性和相频特性计算相干函数,具体为:
与现有技术相比,本发明的有益效果是:
本方案在检测衬板装配质量时不需要视觉系统,有效解决作业环境视觉遮挡、光线恶劣等问题;采用自动化方式实现检测分析,节约人力,整个过程可在测量单元的提醒下完成,减少工作人员工作量。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1是基于频响函数的装配质量检测系统示意图;
图2(a)-(b)分别是衬板测点布置示意图;
图3是冲击激励单元示意图;
图4是基于频响函数的装配质量检测方法流程图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本发明使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
在一个或多个实施方式中公开的一种基于频响函数的装配质量检测系统,如图1所示,包括:冲击激励单元、被测结构件、测量单元、分析单元及输出显示单元。
其中,冲击激励单元主要通过力锤或振动台等产生振动信号;
被测结构件是磨机衬板,包括长方形和三角形两种型号,如图2(a)和图2(b)所示,根据衬板的型号,在衬板上选择合适的测点;测点的选择以反映被测结构件的构型为原则。冲击激励单元在每一个测点产生冲击激励信号。
测量单元主要由传感器和调理电路组成;传感器为三轴加速度传感器,设置在每一个测点位置;调理电路包括信号放大电路和AD转换电路;将传感器采集到的信号放大、转换后送入分析单元;
分析单元主要通过计算机对采集来的信号进行分析处理;
输出显示单元主要包括显示器和打印机等,将分析单元输出的检测结果进行显示和打印。
在另外一些实施方式中,公开了一种基于频响函数的装配质量检测方法,如图4所示,包括以下步骤:
(1)以反映被测结构件的构型为原则,根据检测结构件布置合适的测点,见图二;
(2)安装力锤、三轴加速度传感器等,采用冲击激励法。由于衬板结构巨大,需要大的能量激励,采用多点激励法,如图3,设实际施加在结构上的激励力fm,真正作用于结构件上的作用力为ft,是扣除三维力传感器附加质量和激振系统引入的附加质量所消耗的惯性加速度a,即
(3)设定采样频率,采集加速度等信号,经放大器、AD卡等组成的调理电路转换后,上传到计算机系统;
(4)将采集来的加速度信号利用公式(1)计算出作用于结构件上的作用力,采用频响函数的方法对采集来的信号进行分析处理:
(4-1)对各测试的响应点进行频响函数(FRF)估计;
设系统输入端激励力矩阵f(t)N×1,输出端响应矩阵x(t)O×1,对应频率域的激励力和响应力矩阵分别为F(f)N×1和X(f)O×1。假设系统输入端的噪声矩阵m(t)N×1,对应测量系统输出端矩阵n(t)O×1,对应频率域输入噪声和输出噪声矩阵为M(f)N×1和N(f)O×1。测量系统测得时域输入信号a(t)N×1,输出信号b(t)O×1,频域为A(f)N×1、B(f)O×1,则输入输出之间的关系为:
a(t)=f(t)+m(t) (2)
b(t)=x(t)+n(t) (3)
或者是:
A(f)=F(f)+M(f) (4)
B(f)=X(f)+N(f) (5)
于是,理想情况下,系统频响函数估计式为
X(f)=H(f)F(f) (6)
当有噪声污染时,频响函数估计式变为
将公式(4)(5)代入可得:
B(f)-N(f)=H(f)(A(f)-M(f)) (8)
系统的频响函数估计模型可根据H3估计法得到,即
其中,
H1(f)=Gfx(f)Gff(f)-1 (11)
H2(f)=Gxx(f)Gxf(f)-1 (13)
G为激励信号和响应信号的各态历经随机信号。
(4-2)对频响函数进行幅频、相频及相干特性的分析;
(4-2-1)根据估计的频响函数,求出幅频特性和相频特性,公式(15)(16),并根据幅频和相频曲线进行;
θ(w)=∠G(jw) (16)
(4-2-2)计算相干函数,通过相干函数评定频响函数估计的可信度,可估算出在频域内总输出中真正输入信号产生的输出所占的比例。
(4-2-3)将幅频特性、相频特性及相干特性曲线输入装配质量库,采用模板匹配法与装配装量对应的特性曲线进行判别。
(5)根据判别结果输出显示打印装配质量结果。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (2)
1.一种基于频响函数的装配质量检测方法,采用基于频响函数的装配质量检测系统,包括:冲击激励单元、测量单元、分析单元以及输出显示单元;在被测衬板上设置至少一个测点,冲击激励单元在每一个测点位置产生振动信号;测量单元包括加速度传感器,在衬板上与冲击激励单元作用面相对的一面的每一个测点位置设置加速度传感器,采集加速度信号并传送至分析单元,所述分析单元对采集到的数据进行分析,判断装配质量是否合格;分析单元将分析结果传送至输出显示单元进行显示;其特征在于,包括:
以反映被测衬板的构型为原则,在被测衬板上布置合适的测点;其中,所述衬板为球磨机筒体内的衬板;
设定采集采样频率,采用多点冲击激励法,采集每一个测点的加速度信号;
根据采集到的加速度信号计算出真正作用于被测衬板上的作用力;
对各测点进行频响函数估计;
根据估计的频响函数,求出幅频特性和相频特性;
根据幅频特性和相频特性计算相干函数,通过相干函数评定频响函数估计的可信度;
将幅频特性、相频特性及相干函数曲线输入装配质量库,采用模板匹配法与装配质量对应的特性曲线进行匹配,如果匹配成功则装配质量合格,否则不合格;
根据幅频特性和相频特性计算相干函数,具体为:
S x (f)输入谱密度、S y (f)输出谱密度、S xy (f)条件谱密度;
所述对各测点进行频响函数估计,具体为:
设输入端激励力矩阵f(t) N×1 ,输出端响应矩阵x(t) O×1 ,对应频率域的激励力和响应力矩阵分别为F(f) N×1 和X(f) O×1 ;
假设输入端的噪声矩阵m(t) N×1 ,对应测量单元输出端矩阵n(t) O×1 ,对应频率域输入噪声和输出噪声矩阵为M(f) N×1 和N(f) O×1 ;
测量单元测得时域输入信号a(t) N×1 ,输出信号b(t) O×1 ,频域为A(f) N×1 、B(f) O×1 ,则输入输出之间的关系为:
a(t)=f(t)+m(t)
b(t)=x(t)+n(t)
或者
A(f)=F(f)+M(f)
B(f)=X(f)+N(f)
则理想情况下,系统频响函数估计式为:X(f)=H(f)F(f);
系统的频响函数估计模型根据H3估计法得到,即
其中,
H
1
(f)=G
fx
(f)G
ff
(f)
-1
H
2
(f)=G
xx
(f)G
xf
(f)
-1
其中,G fx (f)实测输入信号和实测输出信号的互谱均值,G ff (f)实测输入信号的的自谱均值,G ab (f)实测输入噪声信号的自谱均值,G aQ (f)系统输入噪声的自谱均值,G xx (f)实测输出信号的自谱均值,G xf (f)实测输出信号和实测输入信号的互谱均值,G bb (f)实测输出噪声的自谱均值,G bQ (f)系统输出信号的自谱均值;
根据估计的频响函数,求出幅频特性和相频特性,具体为:
θ(w)=∠G(jw)
其中,Φ(jw)为幅频特性、G(jw)为频率响应、X是信号的输入、Y是信号的输出、θ(w)是信号的相频特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811438317.8A CN109444263B (zh) | 2018-11-27 | 2018-11-27 | 一种基于频响函数的装配质量检测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811438317.8A CN109444263B (zh) | 2018-11-27 | 2018-11-27 | 一种基于频响函数的装配质量检测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444263A CN109444263A (zh) | 2019-03-08 |
CN109444263B true CN109444263B (zh) | 2021-07-20 |
Family
ID=65554443
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811438317.8A Expired - Fee Related CN109444263B (zh) | 2018-11-27 | 2018-11-27 | 一种基于频响函数的装配质量检测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444263B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110243941A (zh) * | 2019-07-08 | 2019-09-17 | 西安航天化学动力有限公司 | 一种固体火箭发动机界面脱粘缺陷的无损检测方法 |
CN111122085B (zh) * | 2019-12-30 | 2021-11-16 | 武昌船舶重工集团有限公司 | 一种基于功率分布特性的结构装配质量的评价方法 |
CN114242434B (zh) * | 2021-12-14 | 2024-03-01 | 天航长鹰(江苏)科技有限公司 | 一种具备振动频响测试的变压器组装设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107024141A (zh) * | 2017-05-23 | 2017-08-08 | 中北大学 | 供输弹机构装配质量的声振监测与缺陷定位方法 |
CN108151870A (zh) * | 2017-11-30 | 2018-06-12 | 武昌船舶重工集团有限公司 | 一种基于频响函数的施工质量问题检测方法 |
CN108956122A (zh) * | 2018-08-22 | 2018-12-07 | 清华大学 | 一种基于结构动力学特性的装配质量检测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100465609C (zh) * | 2006-12-30 | 2009-03-04 | 北京航空航天大学 | 数控机床加工动力学特性测试分析系统 |
CN102288374A (zh) * | 2011-07-22 | 2011-12-21 | 哈尔滨工业大学 | 同时识别多点随机载荷的试验平台及试验方法 |
CN104502450B (zh) * | 2014-10-21 | 2018-03-13 | 南京航空航天大学 | 简易管道安装应力识别方法 |
CN107767975B (zh) * | 2017-10-17 | 2020-01-21 | 中北大学 | 一种核动力装置关重部件质量性能的监测与故障诊断方法 |
-
2018
- 2018-11-27 CN CN201811438317.8A patent/CN109444263B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107024141A (zh) * | 2017-05-23 | 2017-08-08 | 中北大学 | 供输弹机构装配质量的声振监测与缺陷定位方法 |
CN108151870A (zh) * | 2017-11-30 | 2018-06-12 | 武昌船舶重工集团有限公司 | 一种基于频响函数的施工质量问题检测方法 |
CN108956122A (zh) * | 2018-08-22 | 2018-12-07 | 清华大学 | 一种基于结构动力学特性的装配质量检测方法 |
Non-Patent Citations (2)
Title |
---|
一种基于冲击响应的装配式板梁桥铰接缝病害动力评估方法;战家旺 等;《土木工程学报》;20180630;第51卷(第06期);第103-110页 * |
机械螺栓连接状态监测和辨识方法研究进展;徐超 等;《强度与环境》;20090430;第36卷(第02期);第28-36页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109444263A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444263B (zh) | 一种基于频响函数的装配质量检测系统及方法 | |
CN110389170B (zh) | 基于Lamb波成像的列车部件裂纹损伤检测方法和系统 | |
Rivola et al. | Bispectral analysis of the bilinear oscillator with application to the detection of fatigue cracks | |
Avitabile et al. | Prediction of full field dynamic strain from limited sets of measured data | |
US8401820B2 (en) | In situ health monitoring of piezoelectric sensors | |
CN106568563A (zh) | 一种可定量激励的主轴固有频率多点测试系统 | |
CN106596723A (zh) | 一种多层复合材料结构力学参数声学检测方法 | |
CN103234742A (zh) | 振动筛减振弹簧故障诊断方法 | |
US6920790B2 (en) | Apparatus for in-situ nondestructive measurement of Young's modulus of plate structures | |
Zanarini | On the role of spatial resolution in advanced vibration measurements for operational modal analysis and model updating | |
CN202362104U (zh) | 一种液压破碎锤冲击能检测装置 | |
CN106324097A (zh) | 一种基于声波传感测定材料碰撞恢复系数的方法 | |
CN106932162A (zh) | 轨道动刚度测试方法及系统 | |
CN102901611B (zh) | 一种基于模态试验的星载天线结构损伤快速排查方法 | |
CN106645399B (zh) | 一种复合材料损伤检测评估方法和系统 | |
CN201993345U (zh) | 一种混凝土结构无损测试仪 | |
CN108593769B (zh) | 利用多通道控制弹性波的归一化信号激振无损检测方法 | |
CN114084764B (zh) | 一种电梯乘运质量检测方法及检测系统 | |
Pasquali et al. | Advanced System Identification of Plates Using a Higher-Order-Spectral Approach: Theory and Experiment | |
CN103759951A (zh) | 可同时测量多台农机车的制动性能测试仪及测试方法 | |
CN106289961A (zh) | 一种基于声波传感测定材料碰撞恢复系数的装置 | |
CN202362112U (zh) | 一种爆震传感器测试装置 | |
US6575036B1 (en) | Method for in-situ nondestructive measurement of Young's modulus of plate structures | |
Zhao et al. | A method for structural damage identification using residual force vector and mode shape expansion | |
CN220040202U (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 |
Granted publication date: 20210720 Termination date: 20211127 |
|
CF01 | Termination of patent right due to non-payment of annual fee |