CN116148350A - 一种基于锤击法的快速获取最佳频响函数的方法 - Google Patents
一种基于锤击法的快速获取最佳频响函数的方法 Download PDFInfo
- Publication number
- CN116148350A CN116148350A CN202310130514.8A CN202310130514A CN116148350A CN 116148350 A CN116148350 A CN 116148350A CN 202310130514 A CN202310130514 A CN 202310130514A CN 116148350 A CN116148350 A CN 116148350A
- Authority
- CN
- China
- Prior art keywords
- signal
- tap
- ith
- frequency response
- function
- 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
Links
- 238000005316 response function Methods 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000004044 response Effects 0.000 claims abstract description 21
- 230000005284 excitation Effects 0.000 claims abstract description 20
- 238000001228 spectrum Methods 0.000 claims abstract description 16
- 230000001133 acceleration Effects 0.000 claims abstract description 15
- 238000013139 quantization Methods 0.000 claims abstract description 10
- 230000001427 coherent effect Effects 0.000 claims description 8
- 238000010079 rubber tapping Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 abstract description 6
- 238000002474 experimental method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- 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
- G01N2291/0235—Plastics; polymers; soft materials, e.g. rubber
-
- 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/26—Scanned objects
- G01N2291/269—Various geometry objects
- G01N2291/2692—Tyres
Abstract
本发明涉及一种基于锤击法的快速获取最佳频响函数的方法,属于轮胎测试领域。对激励信号和加速度响应信号的时域数据进行分段;对各段信号施加窗函数,其中激励信号施加力窗、加速度响应信号施加指数衰减窗;对加窗后的每一段信号进行离散傅里叶变换;计算所述的激励信号和加速度响应信号两信号的互功率谱以及各自的自功率谱,为计算频响函数的H1估计和H2估计以及相干函数做铺垫;依据相干函数为量化指标,择优选择第一次敲击信号;以所找的第一次敲击信号为基础,依据相干函数为量化指标,寻找其余两次敲击信号;依据所找的三次敲击信号确定最终的频响函数。只需要保持同一力度连续敲击轮胎指定位置多次后,能够使用所有的敲击数据获得最优频响函数。
Description
技术领域
本发明涉及轮胎频响函数测试方法领域,详细地讲是基于锤击法的快速获取轮胎胎面至轮辋的最佳频响函数的方法。
背景技术
众所周知,在轮胎频响函数测试模态实验中,激励直接关系到实验的成败,通常采用激振器激励或者锤击法进行试验。其中,采用敲击锤进行结构动态试验的方法我们称之为锤击法。相对于轮胎这种中小型结构,锤击法是一种方便快速的模态试验方法。
锤击法最大的好处在于实验速度快,周期短,但也有其缺点,它对敲击力有较高的需求,不仅要求敲击点位置可重复,敲击方向正确,还要保证敲击力度一致性。除此外,还要确保每次敲击不双击,不过载,力的时域信号呈尖峰状。否则力太小会使得信噪降低比,力太大会使结构出现非线性响应,严重地会导致信号失真。
对于锤击法而言,第一次敲击至关重要,它决定结构的频响函数的基本形态,基于第一次敲击信号,通过比对相干函数的一致性,选取频响函数形态一致的敲击信号,最终确定平均后的频响函数。
通常选取第一次敲击的做法是,连续敲击多次,能够重复出现多次类似的频响函数,并且根据经验判断频响函数的良好性,这样的敲击力度、位置和方向才能确定为第一次敲击。此过程过于繁琐冗余,在把重复性高的敲击选定为第一次敲击后,那么之前与之类似的敲击无法获得,敲击利用率不高。
发明内容
为了克服现有技术锤击法选取第一次敲击费时多、利用率低的不足,本发明提供一种敲击利用率高,且快速便捷得到第一次敲击并能找到后续敲击最终求出频响函数的方法。
本发明解决其技术问题所采用的技术方案是:一种基于锤击法的快速获取最佳频响函数的方法,其特征在于,包括以下步骤:
(1)分别对实际敲击的力的激励信号和加速度响应信号的时域数据进行分段;
(2)分别对各段信号施加窗函数,其中激励信号施加力窗(即可变宽度矩形窗)、加速度响应信号施加指数衰减窗(窗的宽度视信号的长度或衰减程度而定);
(3)分别对加窗后的每一段信号进行离散傅里叶变换;
(4)计算所述的激励信号和加速度响应信号两信号的互功率谱以及各自的自功率谱,为计算频响函数的H1估计和H2估计以及相干函数做铺垫;
(5)依据相干函数为量化指标,择优选择第一次敲击信号;
(6)以所找的第一次敲击信号为基础,仍依据相干函数为量化指标,寻找其余两次敲击信号;
(7)依据所找的三次敲击信号确定最终的频响函数。
本发明还可通过如下措施来实现:
根据所述的实际敲击的力的激励信号,选择合适的力窗的长度;根据加速度响应信号的衰减趋势,选择指数衰减窗合适的衰减率b。
所述的步骤(4)中从两信号的互功率谱以及各自的自功率谱得到频响函数的H1估计和H2估计以及相干函数γ的步骤如下:
所述的步骤(5)中相干函数这一量化指标,也就是后文所定义的相干质量因子,具体是指:
所述的步骤(5)中择优选择第一次敲击信号的步骤如下:
首先,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号频响函数的H1估计、H2估计以及相干函数;
其次,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号的相干质量因子;
接下来,统计相干质量因子高于一定值(自己设定,例如0.8)的次数,次数最多者将作为第一次敲击,记为i1。
所述的步骤(6)中选择第二次敲击信号的步骤如下:
首先,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的频响函数的H1估计、H2估计以及相干函数;
其次,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的相干质量因子;
接下来,选出相干质量因子最大者将作为第二次敲击,记为i2。
所述的步骤(6)中根据选定的第一次敲击信号i1和第二次敲击信号i2,确定出第三次敲击信号的步骤如下:
首先,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的频响函数的H1估计、H2估计以及相干函数;
其次,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的相干质量因子;
接下来,选出相干质量因子最大者将作为第三次敲击,记为i3。
所述的步骤(7)中基于第i1、i2、i3次敲击信号相干平均后,确定整体的频响函数的步骤如下:
本发明的有益效果是,基于所有敲击信号,通过评价相干函数质量选取第一次敲击和后续敲击,得到最终频响函数,敲击利用率高,可快速便捷得到第一次敲击并能找到后续敲击,使得所有有效敲击信号都能用于计算最终的频响函数,则可以大大提高锤击法的敲击利用率。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明的基于锤击法快速获取最佳频响函数方法的流程图;
图2是本发明锤击法获得的力的激励信号和加速度的响应信号;
图3是本发明获取第一次敲击信号的过程;
图4是本发明获取后续敲击信号的过程;
图5是本发明相干平均后的频响函数的对比结果。
具体实施方式
本发明所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明所需的敲击信号数据要求不能出现双击、过载等严重情况,否则应剔除此数据。
结合图1,本发明方法的基本实施步骤为:
(1)分别对力激励信号和加速度响应信号的时域数据进行分段;
(2)分别对各段信号施加窗函数,其中激励信号施加力窗(即可变宽度矩形窗)、加速度响应信号施加指数衰减窗(窗的宽度视信号的长度或衰减程度而定);
(3)分别对加窗后的每一段信号进行离散傅里叶变换;
(4)计算两信号的互功率谱以及各自的自功率谱,为计算频响函数的H1估计和H2估计以及相干函数做铺垫;
(5)依据相干函数为量化指标,择优选择第一次敲击信号;
(6)以所找的第一次敲击信号为基础,仍依据相干函数为量化指标,寻找其余两次敲击信号;
(7)依据所找的三次敲击信号确定最终的频响函数。
下面举例对本发明做出更详细的描述。
假设每次敲击的三角形脉冲激励信号形如:
衰减正弦震荡加速度响应信号形如:
首先将两信号按图1所示分成m段,每段长度为N,得到,
xi(n)=x(n+(i-1)N),n=0,1,...N-1,i=12...m (3)
yi(n)=x(n+(i-1)N),n=0,1,...N-1,i=12...m (4)
为了减少频谱泄露,对分段后的数据进行加窗,
x′i(n)=xi(n)w1(n),y′i(n)=xi(n)w2(n),n=0,1,...N-1. (5)
其中窗函数w1(n)是力窗,窗函数w2(n)是指数窗。
然后对各个分段信号x′i(n)和y′i(n)进行离散傅里叶变换:
这里ω代表离散傅里叶变换所得到的离散频率。
接下来计算各个分段信号在每个频率ω处的自功率谱和互功率谱。其中各段激励信号是x′i(n)的自功率谱为:
各段响应信号y′i(n)的自功率谱为:
各段激励信号x′i(n)和各段响应信号y′i(n)的互功率谱为:
下一步,择优选定第一次敲击信号,详细流程见图3。
这里,我将介绍一种关于相干函数的评价指标,称之为相干质量因子γ,定义如下:
相干函数曲线与两坐标轴所围成的面积与频率最大值ωmax之比(比值在0-1之间),比值越接近1,质量越高。
为选定第一次敲击信号,首先,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号频响函数的H1估计、H2估计以及相干函数,例如下面是第i次和第j次敲击信号的结果。
其次,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号的相干质量因子:
接下来,统计相干质量因子高于一定值(自己设定,例如0.8)的次数,次数最多者将作为第一次敲击。例如以第i1次敲击信号为基础,余下m-1次敲击信号与第i1次敲击信号的相干质量因子高于0.8的次数最高,那么选定第i1次敲击信号作为第一次敲击。
接着,根据选定的第一次敲击信号i1,确定出第二次敲击信号,详细流程见图4。
首先,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的频响函数的H1估计、H2估计以及相干函数,例如下面是第i1次和第j次敲击信号的结果。
其次,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的相干质量因子:
接下来,选出相干质量因子最大者将作为第二次敲击。例如以第i1次敲击信号为基础,第i2次敲击信号与第i1次敲击信号的相干质量因子最大,那么选定第i2次敲击信号作为第二次敲击。
接着,根据选定的第一次敲击信号i1和第二次敲击信号i2,确定出第三次敲击信号。
首先,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的频响函数的H1估计、H2估计以及相干函数,例如下面是第i1次、第i2次和第j次敲击信号的结果。
其次,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的相干质量因子:
接下来,选出相干质量因子最大者将作为第三次敲击。例如以第i1次和第i2次敲击信号为基础,第i3次敲击信号与第i1次和第i2次敲击信号的相干质量因子最大,那么选定第i3次敲击信号作为第三次敲击。
最后,基于第i1、i2、i3次敲击信号相干平均后,确定整体的频响函数Hv估计为:
图5给出了基于不同敲击信号相干平均后频响函数的H1、H2估计。至此,我们已经测量出频响函数,基于频响函数后续可以进行模态测试、固有频率测试等一系列测试。
本发明的优势在于只需要保持同一力度连续敲击轮胎指定位置多次后,能够使用所有的敲击数据,根据上述步骤即可获得最优频响函数。本发明方法既保持了锤击法实验速度快,周期短的特点,又大大提高了敲击信号利用率,避免了数据的不必要浪费。
Claims (8)
1.一种基于锤击法的快速获取最佳频响函数的方法,其特征在于,包括以下步骤:
(1)分别对实际敲击的力的激励信号和加速度响应信号的时域数据进行分段;
(2)分别对各段信号施加窗函数,其中激励信号施加力窗、加速度响应信号施加指数衰减窗;
(3)分别对加窗后的每一段信号进行离散傅里叶变换;
(4)计算所述的激励信号和加速度响应信号两信号的互功率谱以及各自的自功率谱,为计算频响函数的H1估计和H2估计以及相干函数做铺垫;
(5)依据相干函数为量化指标,择优选择第一次敲击信号;
(6)以所找的第一次敲击信号为基础,仍依据相干函数为量化指标,寻找其余两次敲击信号;
(7)依据所找的三次敲击信号确定最终的频响函数。
2.根据权利要求1所述基于锤击法的快速获取最佳频响函数的方法,其特征在于,根据实际敲击的力的激励信号,选择合适的力窗的长度;根据加速度响应信号的衰减趋势,选择指数衰减窗合适的衰减率b。
5.根据权利要求1所述基于锤击法的快速获取最佳频响函数的方法,其特征在于所述的步骤(5)中择优选择第一次敲击信号的步骤如下:
首先,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号频响函数的H1估计、H2估计以及相干函数;
其次,以第i次敲击信号为基础,遍历计算第i次敲击信号和余下m-1次敲击信号的相干质量因子;
接下来,统计相干质量因子高于一定值的次数,次数最多者将作为第一次敲击,记为i1。
6.根据权利要求1所述基于锤击法的快速获取最佳频响函数的方法,其特征在于所述的步骤(6)中选择第二次敲击信号的步骤如下:
首先,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的频响函数的H1估计、H2估计以及相干函数;
其次,以第i1次敲击信号为基础,遍历计算第i1次敲击信号和余下m-1次敲击信号的相干质量因子;
接下来,选出相干质量因子最大者将作为第二次敲击,记为i2。
7.根据权利要求1所述基于锤击法的快速获取最佳频响函数的方法,其特征在于所述的步骤(6)中根据选定的第一次敲击信号i1和第二次敲击信号i2,确定出第三次敲击信号的步骤如下:
首先,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的频响函数的H1估计、H2估计以及相干函数;
其次,以第i1次和第i2次敲击信号为基础,遍历计算第i1次和第i2次敲击信号和余下m-2次敲击信号的相干质量因子;
接下来,选出相干质量因子最大者将作为第三次敲击,记为i3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310130514.8A CN116148350A (zh) | 2023-02-17 | 2023-02-17 | 一种基于锤击法的快速获取最佳频响函数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310130514.8A CN116148350A (zh) | 2023-02-17 | 2023-02-17 | 一种基于锤击法的快速获取最佳频响函数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116148350A true CN116148350A (zh) | 2023-05-23 |
Family
ID=86359714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310130514.8A Pending CN116148350A (zh) | 2023-02-17 | 2023-02-17 | 一种基于锤击法的快速获取最佳频响函数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116148350A (zh) |
-
2023
- 2023-02-17 CN CN202310130514.8A patent/CN116148350A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hudson et al. | Knock measurement for fuel evaluation in spark ignition engines | |
CN102422154A (zh) | 一种结构损伤检测系统、设备以及结构损伤检测方法 | |
Lee et al. | Development of enhanced Wigner–Ville distribution function | |
CN106556647A (zh) | 一种冲击回波数据处理方法 | |
CN108896274B (zh) | 一种基于子集窗长度优化算法的分布式光纤应变解调方法 | |
CN110186556B (zh) | 一种新能源汽车电机台架试验噪声评价方法 | |
CN111639541A (zh) | 基于频率变化率的自适应同步压缩时频分析方法 | |
CN113901580A (zh) | 一种车辆减振器异响预测方法及系统 | |
CN110596247A (zh) | 一种温度变化环境下的超声结构健康监测方法 | |
CN116148350A (zh) | 一种基于锤击法的快速获取最佳频响函数的方法 | |
CN110936977A (zh) | 一种高速无砟轨道结构扣件松脱检测方法 | |
JP4520146B2 (ja) | 突起剛性変化に基づくタイヤピッチ列の分析方法 | |
CN110399683A (zh) | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 | |
Brandt et al. | Impact excitation processing for improved frequency response quality | |
JP6053335B2 (ja) | 波動方程式の反転による内燃機関のノッキング強度の推定方法及び機関の燃焼の制御方法 | |
DE3607913C1 (en) | Method for locating leaks in pipelines | |
KR102067197B1 (ko) | 주파수 영역 분석을 이용한 타이어 패턴소음 예측 방법 | |
JP4431023B2 (ja) | タイヤ振動特性測定方法 | |
CN111272335B (zh) | 一种机器人多维力/力矩传感器动态特性辨识方法 | |
CN111931288B (zh) | 一种火工冲击响应时频域分布相似性测度方法 | |
JP2000346755A (ja) | タイヤ摩耗の測定方法 | |
JP6501049B1 (ja) | 評価方法及び評価装置、並びに記録媒体 | |
US20210087986A1 (en) | Method and device for determining combustion noise parameters | |
Gaidis et al. | New test for determining fundamental frequencies of concrete | |
CN114384156B (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 |