基于時(shí)間反轉(zhuǎn)法的汽車車門Rattle異響噪聲源定位研究*

趙衛(wèi)東，徐鑫蔚，宋 睿，楊明亮，齊 潮

（1.江蘇大學(xué)汽車與交通工程學(xué)院，鎮(zhèn)江 212013；2.中國(guó)汽車技術(shù)研究中心，天津 300300；3.西南交通大學(xué)機(jī)械工程學(xué)院，成都 610031；4.安徽理工大學(xué)礦業(yè)工程學(xué)院，淮南 232001）

前言

汽車異響噪聲（buzz，squeak and rattle，BSR）是指汽車運(yùn)行過(guò)程中，由于結(jié)構(gòu)缺陷、材料對(duì)的不相容和不良的幾何匹配［1-2］等所導(dǎo)致的非正常、沒(méi)有規(guī)律的聲音。近年來(lái)，汽車輕量化使異響出現(xiàn)的概率大大增加［1］，在傳統(tǒng)噪聲振動(dòng)源和傳遞路徑得到很好控制的背景下，汽車異響噪聲已成為決定汽車聲品質(zhì)和消費(fèi)者滿意度的關(guān)鍵因素［3-5］。異響噪聲源的準(zhǔn)確定位是異響控制的重要前提。由于BSR的特征和持續(xù)時(shí)間均無(wú)規(guī)律，是一種非線性很強(qiáng)的物理現(xiàn)象，因此難以用成熟的理論和方法分析處理。目前工程應(yīng)用中異響噪聲源定位仍以主觀評(píng)價(jià)為主［6-7］。主觀評(píng)價(jià)法對(duì)工程師的經(jīng)驗(yàn)要求高，試驗(yàn)工作量大，且易引起爭(zhēng)議。發(fā)展異響噪聲源客觀定位技術(shù)與方法已成為行業(yè)的迫切需求。

目前較為成熟的客觀聲源定位方法主要有聲強(qiáng)測(cè)量法［8-9］、波束形成法［10-11］和聲全息法［12-13］。聲強(qiáng)測(cè)量法只適用于穩(wěn)態(tài)聲源，無(wú)法應(yīng)用于非穩(wěn)態(tài)異響噪聲源定位；波束形成法在應(yīng)用于中高頻遠(yuǎn)場(chǎng)信號(hào)源識(shí)別定位時(shí)精度較高，但對(duì)低頻異響噪聲源分辨率較低［14］；聲全息法雖適用于低頻噪聲信號(hào)的定位，但實(shí)施時(shí)需獲取整個(gè)全息面的聲信號(hào)，布置陣元較多，系統(tǒng)復(fù)雜，測(cè)試設(shè)備和成本較高，阻礙了該技術(shù)在異響噪聲源識(shí)別定位技術(shù)領(lǐng)域的推廣和應(yīng)用。探索操作簡(jiǎn)單、成本低、定位精度高的異響噪聲源客觀定位技術(shù)已成為領(lǐng)域前沿。

時(shí)間反轉(zhuǎn)法（time reversal，TR），作為一種客觀信號(hào)處理方法，源于光學(xué)中的相位共軛技術(shù)［15-18］。信號(hào)在頻域中的相位共軛對(duì)應(yīng)時(shí)域中的時(shí)間反轉(zhuǎn)，利用收發(fā)信道互易性原理［15-18］，在信號(hào)接收處發(fā)射時(shí)間反轉(zhuǎn)信號(hào)，能夠克服多徑效應(yīng)，在信號(hào)源位置實(shí)現(xiàn)信號(hào)在時(shí)間和空間上的聚焦，從而達(dá)到定位信號(hào)源的目的?；诼暡ê凸獾膫鬏斚嗨菩?，F(xiàn)ink和Wu等［19-21］最早將時(shí)間反轉(zhuǎn)法引入到聲學(xué)領(lǐng)域；在水聲通訊［22-24］、材料結(jié)構(gòu)損傷識(shí)別［25-26］等方面的成功應(yīng)用中，體現(xiàn)出適用于非穩(wěn)態(tài)低頻聲振信號(hào)，所需定位信號(hào)少等技術(shù)優(yōu)勢(shì)，為汽車異響噪聲源的定位提供了一種新的途徑。

Kwak等［27］基于時(shí)間反轉(zhuǎn)定位原理，成功實(shí)現(xiàn)了白車身沖擊異響噪聲源的定位。雖然Kwak等人的開(kāi)創(chuàng)性工作充分證實(shí)了TR定位異響噪聲源的可行性，但其建立的異響噪聲源定位模型在實(shí)施時(shí)，須采用譜元法對(duì)被測(cè)對(duì)象進(jìn)行精準(zhǔn)的力學(xué)分析，以求得振動(dòng)信號(hào)的傳遞函數(shù)，存在建模工作量大、測(cè)試工作周期長(zhǎng)等不足。

汽車鈑金件易發(fā)生變形，是Rattle異響噪聲的高發(fā)部件。在鈑金件中，振動(dòng)波的傳播可近似為L(zhǎng)amb波，鑒于Lamb波在介質(zhì)中的傳遞函數(shù)已有精確的解析解［28］，故理論上以該解析解為基礎(chǔ)，對(duì)發(fā)生在鈑金件中的異響噪聲源進(jìn)行時(shí)間反轉(zhuǎn)定位，可規(guī)避Kwak等人研究中復(fù)雜的譜元法建模過(guò)程，使TR異響噪聲定位模型的構(gòu)建及應(yīng)用大為簡(jiǎn)化。

但具體分析Lamb波傳遞函數(shù)可知［28］，該函數(shù)包含n階對(duì)稱與反對(duì)稱模式的幅值項(xiàng)和相位項(xiàng)，無(wú)法實(shí)際應(yīng)用于振動(dòng)時(shí)間反轉(zhuǎn)信號(hào)的疊加計(jì)算；以實(shí)際異響振動(dòng)信號(hào)做時(shí)間反轉(zhuǎn)處理，計(jì)算量大，且寬頻Lamb波在傳播時(shí)的頻散現(xiàn)象將導(dǎo)致波速計(jì)算產(chǎn)生較大誤差，進(jìn)而影響定位精度。

基于上述分析，為應(yīng)用TR定位原理，實(shí)現(xiàn)汽車鈑金件中異響噪聲源的快速、精準(zhǔn)定位，本文中根據(jù)汽車鈑金件中異響振動(dòng)信號(hào)的特征，將無(wú)限階次的Lamb波傳遞函數(shù)簡(jiǎn)化為有限階次；采用Morlet小波變換提取異響振動(dòng)信號(hào)中的窄帶成分，進(jìn)而依據(jù)時(shí)間反轉(zhuǎn)定位原理，建立Rattle異響噪聲源定位模型。在此基礎(chǔ)上，以某型汽車車門為例，對(duì)模型定位精度進(jìn)行仿真和試驗(yàn)。期望通過(guò)研究，為汽車鈑金件異響噪聲源的定位提供一種操作簡(jiǎn)單、實(shí)施方便的新方法。

1 Rattle異響噪聲源TR聚焦定位模型

1.1 TR信號(hào)聚焦原理

異響振動(dòng)時(shí)間反轉(zhuǎn)信號(hào)（下簡(jiǎn)稱時(shí)反信號(hào)）聚焦定位系統(tǒng)示意圖如圖1所示。設(shè)待測(cè)件P上某處存在異響振動(dòng)源S，激發(fā)頻譜為s(ω)的振動(dòng)信號(hào)，振動(dòng)信號(hào)在P上任意兩點(diǎn)之間的傳遞函數(shù)為h(r，ω)，其中r為兩點(diǎn)之間的距離，ω為振動(dòng)信號(hào)的頻率。為定位異響振動(dòng)源S，在P表面布置n(n≥2)個(gè)振動(dòng)信號(hào)傳感器。

圖1 時(shí)間反轉(zhuǎn)信號(hào)聚焦定位系統(tǒng)示意圖

設(shè)S與i(i=1，2，…n)傳感器之間的距離為ri，則兩者之間振動(dòng)信號(hào)傳遞函數(shù)為h(ri，ω)。

i傳感器接收到的振動(dòng)信號(hào)為

將i傳感器接收到的信號(hào)yi(ω)時(shí)間反轉(zhuǎn)處理得到其共軛信號(hào)(ω)：

設(shè)待測(cè)件P上任意位置U到i傳感器之間的距離為，若在相應(yīng)傳感器位置處將n個(gè)時(shí)間反轉(zhuǎn)信號(hào)同時(shí)發(fā)射，則在U位置處，n個(gè)時(shí)間反轉(zhuǎn)信號(hào)疊加為

當(dāng)U與異響振動(dòng)源S位置重疊時(shí)，r′i=ri，式（3）中h*(ri，ω)?h(r′i，ω)=h*(ri，ω)?h(ri，ω)，表 現(xiàn) 自 相關(guān)，疊加信號(hào)模值增強(qiáng)。當(dāng)U與S位置不重疊時(shí)，因r′i≠ri，h*(ri，ω)?h(r′i，ω)表現(xiàn)弱相關(guān)，疊加信號(hào)模值減小。

通過(guò)上述分析可知，時(shí)間反轉(zhuǎn)信號(hào)在待測(cè)件中傳播具有空間聚焦特性，表現(xiàn)為在異響振動(dòng)源處，疊加信號(hào)模值出現(xiàn)極值［29］。通過(guò)計(jì)算由待測(cè)件幾何形狀所限定的空間內(nèi)各點(diǎn)處的合成信號(hào)模值，即可確定異響振動(dòng)源位置。

1.2 車門Rattle異響振動(dòng)信號(hào)傳播特性

薄板結(jié)構(gòu)中，應(yīng)力波在傳播過(guò)程中以反射與折射的方式與邊界發(fā)生作用，縱波和橫波不斷發(fā)生反射和模式轉(zhuǎn)換，干涉疊加形成Lamb波［30］。單頻或窄帶Lamb波傳播的傳遞函數(shù)［28］可以表述為

式中：aAi(ri，ω)、asi(ri，ω)、kAi、ksi分別表示信號(hào)以對(duì)稱和反對(duì)稱模式進(jìn)行傳播時(shí)的幅值和波數(shù)；kAi=ω/cAi，ksi=ω/csi，cAi、csi分別為信號(hào)以對(duì)稱和反對(duì)稱模式傳播時(shí)的波速，i=0，1，…，n。

由于式（4）所示的函數(shù)包含有n階對(duì)稱和反對(duì)稱模式的幅值項(xiàng)和相位項(xiàng)，因此無(wú)法直接代入式（3）進(jìn)行時(shí)間反轉(zhuǎn)定位計(jì)算。實(shí)際應(yīng)用時(shí)，尚須根據(jù)異響振動(dòng)噪聲信號(hào)特征，對(duì)式（4）進(jìn)行簡(jiǎn)化，將其近似表達(dá)為有限階。

已有研究表明，當(dāng)信號(hào)頻率低于100 kHz時(shí)，Lamb波主要存在A0和S0模式，且A0模式占主導(dǎo)地位［31］。汽車車門Rattle異響振動(dòng)信號(hào)實(shí)測(cè)頻譜如圖2所示。由圖2可見(jiàn)，信號(hào)能量主要集中在0-50 kHz頻率范圍內(nèi)，滿足文獻(xiàn)［31］給出的判據(jù)。

圖2 車門Rattle信號(hào)頻譜

因此，車門Rattle異響振動(dòng)信號(hào)在鈑金件中傳播的頻響函數(shù)可近似為

i傳感器接收到的信號(hào)可表達(dá)為

將n個(gè)Rattle異響振動(dòng)時(shí)間反轉(zhuǎn)信號(hào)在相應(yīng)傳感器處同時(shí)發(fā)射，則在被測(cè)件上距i傳感器為r′i的任意一點(diǎn)U處合成信號(hào)為

式中：c為單頻或窄帶信號(hào)傳播的波速；y*(ω)可通過(guò)對(duì)傳感器采集到的異響振動(dòng)信號(hào)時(shí)間反轉(zhuǎn)處理獲得。

相位項(xiàng)是時(shí)反信號(hào)能否聚焦的決定因素，因此可忽略幅值項(xiàng)的影響，僅保留相位項(xiàng)［28］。此時(shí)式（7）可表述為

對(duì)式（8）進(jìn)行傅里葉反變換，可得待測(cè)件上任意一點(diǎn)U處合成信號(hào)的時(shí)域表達(dá)式為

式中：Const為常數(shù)項(xiàng)；τ為所取固定時(shí)長(zhǎng)；yi(τ-t+r′i/c)表示對(duì)傳感器接收信號(hào)yi(t)進(jìn)行時(shí)間反轉(zhuǎn)處理。

式（9）即為簡(jiǎn)化的窄帶低頻Lamb波時(shí)間反轉(zhuǎn)定位數(shù)學(xué)模型。將檢測(cè)區(qū)域離散為l×w個(gè)單元，求得各單元處合成信號(hào)的模值，由時(shí)間反轉(zhuǎn)信號(hào)的聚焦特性可知，疊加信號(hào)模值出現(xiàn)極值處的單元，即為異響噪聲源位置［31］。由式（9）可知，對(duì)于窄帶低頻Lamb波信號(hào)，TR信號(hào)的疊加計(jì)算可簡(jiǎn)化為對(duì)信號(hào)的相位調(diào)制。

1.3 異響振動(dòng)信號(hào)窄帶成分提取

式（9）所示的異響噪聲源TR定位模型，建立在信號(hào)為窄帶Lamb波基礎(chǔ)之上，因此提取Rattle振動(dòng)信號(hào)中的窄帶成分，是式（9）應(yīng)用的重要前提。此外，以提取的窄帶信號(hào)代替實(shí)際異響振動(dòng)信號(hào)，進(jìn)行信號(hào)時(shí)間反轉(zhuǎn)、疊加計(jì)算，還可消除Lamb波頻散對(duì)定位精度的影響［32］，有效提高定位精度。

使用Morlet小波變換，可提取異響振動(dòng)信號(hào)中所需的窄帶信號(hào)成分［32］，判斷Lamb波信號(hào)到達(dá)時(shí)刻，進(jìn)而計(jì)算其傳播的相速度。

Morlet小波函數(shù)［33-34］定義為

式中ω0為小波中心頻率，其傅里葉變換為

式中Φ(ω)為高斯窗函數(shù)，可用于提取中心頻率為ω0的窄帶信號(hào)。

由于Rattle振動(dòng)信號(hào)可看作由兩個(gè)振幅相同但頻率不同的諧波組成的發(fā)散波［31］，則沿x方向傳播，單位幅值的Rattle信號(hào)可表示為

信號(hào)u(x，t)的Morlet小波變換為

式中：*表示復(fù)函數(shù)的共軛變換；a為尺度因子；b為時(shí)間因子。設(shè)

使用Morlet復(fù)數(shù)小波對(duì)u(x，t)做小波變換，并計(jì)算其模值：

式中b=?kx/?ω=x/c時(shí)模值最大，即Morlet小波變換系數(shù)模值最大值對(duì)應(yīng)的時(shí)刻，即為提取的窄帶信號(hào)到達(dá)時(shí)刻［35］。

1.4 TR異響噪聲源定位實(shí)施流程

將檢測(cè)區(qū)域離散為l×w個(gè)單元，在待測(cè)件表面布置n(n>2)個(gè)加速度傳感器。設(shè)異響振動(dòng)信號(hào)發(fā)出時(shí)刻為T0，傳播到各傳感器的時(shí)間分別為t1，t2，…，tn，t1

根據(jù)式（16），可計(jì)算得到檢測(cè)區(qū)域內(nèi)各單元疊加信號(hào)模值的大小。為了更加直觀地定位異響噪聲源，可對(duì)檢測(cè)區(qū)域進(jìn)行信號(hào)成像，并對(duì)成像圖閾值化處理。TR異響噪聲源定位實(shí)施流程如圖3所示。

圖3 Rattle異響噪聲源TR定位實(shí)施流程

2 仿真分析

為驗(yàn)證上述基于TR的汽車鈑金件異響噪聲源定位方法的可行性，并探究定位精度的影響因素，進(jìn)行了薄板集中力沖擊點(diǎn)定位的仿真。

2.1 仿真模型設(shè)置

使用ABAQUS建立尺寸為100 cm×100 cm×0.08 cm的仿真模型，如圖4所示。以模型左下角為原點(diǎn)建立直角坐標(biāo)系，在坐標(biāo)（40，30，0）處沿Z軸方向施加10 N的沖擊載荷。為在模型全區(qū)域范圍內(nèi)均能獲得理想的信號(hào)聚焦效果，在模型四周選?。?0，90，0）、（90，90，0）、（10，10，0）、（90，10，0）4個(gè)點(diǎn)，設(shè)置為1-4號(hào)振動(dòng)信號(hào)接收點(diǎn)，提取振動(dòng)信號(hào)。

圖4 薄板沖擊點(diǎn)定位仿真模型

薄板材料屬性如表1所示。

表1 薄板材料屬性

2.2 沖擊仿真與信號(hào)處理

采用動(dòng)力學(xué)顯式step，設(shè)置分析步時(shí)長(zhǎng)為4 ms，采樣頻率為96 kHz，進(jìn)行仿真分析。2.56×10-4s時(shí)刻加速度云圖如圖5所示。

圖5 沖擊仿真加速度云圖

提取信號(hào)接收點(diǎn)處振動(dòng)加速度信號(hào)，其中1號(hào)點(diǎn)處Z向振動(dòng)加速度時(shí)域信號(hào)如圖6所示。

圖6 沖擊振動(dòng)加速度信號(hào)

對(duì)1號(hào)點(diǎn)處時(shí)域信號(hào)進(jìn)行Morlet小波變換，得到信號(hào)的時(shí)頻圖，如圖7所示。

由圖7可見(jiàn)，隨著時(shí)間的推移，沖擊信號(hào)波包混疊嚴(yán)重，不利于信號(hào)到達(dá)時(shí)刻的準(zhǔn)確判斷，從而影響異響噪聲源定位。

圖7 沖擊振動(dòng)信號(hào)時(shí)頻圖

為探究波包混疊成因，設(shè)定XZ、YZ面及其對(duì)稱面為低反射邊界，進(jìn)行沖擊仿真對(duì)比研究，其中1號(hào)點(diǎn)接收到的振動(dòng)加速度信號(hào)的時(shí)域和時(shí)頻圖如圖8所示。由圖8（b）可見(jiàn)，低反射邊界條件下沖擊振動(dòng)信號(hào)沒(méi)有波包混疊現(xiàn)象出現(xiàn)。因此可推斷，沖擊信號(hào)波包混疊是由端面反射所致。

圖8 低反射邊界條件下薄板沖擊振動(dòng)加速度信號(hào)

但在實(shí)際工程應(yīng)用中，測(cè)試件均為有限介質(zhì)，故端面反射必然存在。對(duì)比圖7與圖8可知，首達(dá)波包主要成分為沖擊振動(dòng)源直達(dá)信號(hào)，且無(wú)端面反射波疊加。故提取窄帶信號(hào)時(shí)，為減少信號(hào)端面反射對(duì)定位精度的影響，應(yīng)截取首達(dá)波包，并以首達(dá)波包能量峰值頻率為小波中心頻率。

由圖7可見(jiàn)，首達(dá)波包信號(hào)能量主要集中在15 kHz左右，故采用Morlet小波變換提取中心頻率為15 kHz的窄帶信號(hào)。為減少信號(hào)衰減對(duì)定位精度的影響，將提取的信號(hào)模值進(jìn)行歸一化處理，結(jié)果如圖9所示。

圖9 15 kHz窄帶信號(hào)歸一化模值

由1.3節(jié)分析可知，圖9中信號(hào)模值極大值對(duì)應(yīng)時(shí)刻即為信號(hào)到達(dá)時(shí)刻。計(jì)算傳感器與沖擊點(diǎn)位置距離差與信號(hào)到達(dá)時(shí)刻差之間的比值，得到15 kHz窄帶Lamb波傳播平均波速為1 310.91 m/s。

2.3 沖擊點(diǎn)定位仿真

在確定窄帶信號(hào)頻率及其波速的基礎(chǔ)上，任意選取一點(diǎn)（40，70，0）進(jìn)行沖擊點(diǎn)定位仿真。以1 cm×1 cm的單元格尺寸，將模型X-Y平面劃分為100×100個(gè)單元，將計(jì)算所得波速c=1310.91 m/s代入式（16）計(jì)算檢測(cè)區(qū)域合成信號(hào)模值，并對(duì)計(jì)算結(jié)果成像，結(jié)果如圖10和圖11所示。由圖11閾值化處理結(jié)果可知，定位沖擊點(diǎn)位置為（38，72，0），與實(shí)際沖擊點(diǎn)位置（40，70，0）相比，定位誤差為2.83 cm。實(shí)際工程應(yīng)用中，定位誤差不超過(guò)5 cm，即可判定異響發(fā)生部件。仿真結(jié)果表明，本文中提出的基于TR的異響噪聲源定位方法，滿足定位精度要求。

圖10 定位成像圖

圖11 成像圖閾值化結(jié)果

3 車門Rattle異響噪聲源定位試驗(yàn)

在通過(guò)仿真確定窄帶信號(hào)提取原則，并驗(yàn)證基于TR定位異響噪聲源方法可行性的基礎(chǔ)上，進(jìn)行了車門鈑金件Rattle異響噪聲源定位試驗(yàn)。

3.1 試驗(yàn)系統(tǒng)與方法

試驗(yàn)系統(tǒng)主要由國(guó)產(chǎn)某品牌汽車車門、力錘（型號(hào)：086C03，PCB）、LMS信號(hào)采集系統(tǒng)（型號(hào)：LMS SCADAS Mobile，SIEMENS）、單向振動(dòng)加速度傳感器（型號(hào)：4507B，BK）和計(jì)算機(jī)（型號(hào)：X13，Thinkpad）組成，系統(tǒng)結(jié)構(gòu)如圖12所示，試驗(yàn)現(xiàn)場(chǎng)如圖13所示。

圖12 基于TR的汽車車門Rattle異響噪聲源定位試驗(yàn)系統(tǒng)示意圖

圖13 Rattle異響源定位試驗(yàn)現(xiàn)場(chǎng)圖

試驗(yàn)時(shí)，使用力錘沿-Z方向敲擊車門激發(fā)Rattle異響，LMS信號(hào)采集系統(tǒng)獲取振動(dòng)加速度信號(hào)并輸入TR定位模型，進(jìn)行窄帶信號(hào)提取、時(shí)間反轉(zhuǎn)、時(shí)反信號(hào)疊加計(jì)算，進(jìn)而確定異響噪聲源位置。

試驗(yàn)時(shí)，將車門固定在MB零部件異響臺(tái)架上。理論上，增加傳感器數(shù)量可提高噪聲源定位信號(hào)強(qiáng)度；但傳感器過(guò)多，將使試驗(yàn)操作變得復(fù)雜，增加試驗(yàn)和計(jì)算時(shí)長(zhǎng)，且一定程度上會(huì)影響待測(cè)件整體模態(tài)。在進(jìn)行試驗(yàn)時(shí)，綜合考慮試驗(yàn)車門形狀和信號(hào)衰減程度，在車門表面沿待檢測(cè)區(qū)域周邊，安裝6個(gè)單向振動(dòng)加速度傳感器。以左下傳感器所在位置為坐標(biāo)原點(diǎn)，建立直角坐標(biāo)系。檢測(cè)區(qū)域尺寸約為76 cm×35 cm，1-6號(hào)傳感器坐標(biāo)分別為（76，0，0）、（76，35，0）、（0，35，0）、（0，0，0）、（36，35，0）、（36，0，0）。

3.2 試驗(yàn)結(jié)果與分析

3.2.1 信號(hào)采集與處理

試驗(yàn)時(shí)，首先用力錘在確定點(diǎn)（30，10，0）處敲擊，獲得振動(dòng)加速度信號(hào)，從而確定信號(hào)提取的小波中心頻率為20 kHz，計(jì)算得到平均波速為595.10 m/s。在此基礎(chǔ)上，在隨機(jī)選取點(diǎn)（35，15，0）處敲擊，進(jìn)行定位試驗(yàn)。截取1號(hào)加速度傳感器1.33-1.34 s時(shí)段的振動(dòng)加速度信號(hào)，如圖14所示。

圖14 1號(hào)加速度傳感器信號(hào)

對(duì)圖14所示振動(dòng)加速度信號(hào)進(jìn)行Morlet小波變換，提取信號(hào)首達(dá)波包中頻率為20 kHz的窄帶信號(hào)，將提取的信號(hào)模值進(jìn)行歸一化處理，結(jié)果如圖15所示。

圖15 20 kHz窄帶信號(hào)模值歸一化結(jié)果

3.2.2 異響振動(dòng)信號(hào)發(fā)出時(shí)刻T0的判定

由1.4節(jié)的論述可知，在進(jìn)行異響噪聲源定位時(shí)，應(yīng)以tn-T0為時(shí)窗，截取各傳感器信號(hào)，時(shí)間反轉(zhuǎn)處理后進(jìn)行疊加計(jì)算。當(dāng)時(shí)窗長(zhǎng)度小于或大于tn-T0時(shí)，時(shí)反信號(hào)無(wú)法在源處聚焦，從而無(wú)法精確定位異響噪聲源位置；當(dāng)時(shí)窗長(zhǎng)度等于tn-T0時(shí)，時(shí)反信號(hào)恰好在異響噪聲源處疊加聚焦，且異響噪聲源附近區(qū)域疊加信號(hào)模值較小，從而準(zhǔn)確定位異響噪聲源位置。

在進(jìn)行仿真時(shí)，仿真沖擊開(kāi)始時(shí)刻即為異響振動(dòng)信號(hào)發(fā)出時(shí)刻，T0=0。不同于仿真，實(shí)際Rattle異響是偶發(fā)的，異響振動(dòng)信號(hào)發(fā)出時(shí)刻T0未知，這將導(dǎo)致時(shí)窗長(zhǎng)度tn-T0確定困難，進(jìn)而影響異響源位置的準(zhǔn)確定位。

工程應(yīng)用中，信號(hào)發(fā)出時(shí)刻T0可通過(guò)四點(diǎn)圓弧定位法計(jì)算求得［36］。但該方法計(jì)算工作量大，且因測(cè)量誤差的存在，計(jì)算精度較低。根據(jù)上述時(shí)窗長(zhǎng)度對(duì)定位精度影響的分析可知，當(dāng)所取時(shí)窗長(zhǎng)度恰好等于tn-T0時(shí)，疊加信號(hào)模值較大區(qū)域連續(xù)且面積最小；而所取時(shí)窗長(zhǎng)度大于或小于tn-T0時(shí)，疊加信號(hào)模值較大區(qū)域離散且面積較大。因此，可通過(guò)疊加信號(hào)成像區(qū)域連續(xù)性和面積大小，判定信號(hào)發(fā)出時(shí)刻T0。

具體實(shí)施時(shí)，以檢測(cè)區(qū)域信號(hào)可能傳播的最遠(yuǎn)距離與波速的比值為上限Tu，n個(gè)傳感器最先接收到信號(hào)的時(shí)刻為下限Td，在區(qū)間[tn-Tu，tn-Td]范圍內(nèi)，遞增或遞減賦值T0，比較成像圖，從而判定T0，并定位異響源位置。

不同T0時(shí)刻，車門Rattle異響源定位成像圖和閾值化處理結(jié)果如圖16所示。由圖16可見(jiàn)：T0=1.3345 s，T0=1.3347 s時(shí)，疊加信號(hào)模值較大區(qū)域離散且面積更大；T0=1.3346 s時(shí)，疊加信號(hào)模值較大區(qū)域連續(xù)且面積最小；因此判定，實(shí)際振動(dòng)信號(hào)發(fā)出時(shí)刻T0應(yīng)為1.334 6 s。對(duì)圖16（c）進(jìn)行閾值化處理，結(jié)果如圖16（d）所示，由圖16（d）定位異響噪聲源位置為（34，16，0）。

圖16 不同T0時(shí)刻車門Rattle異響噪聲源定位成像圖

3.2.3 定位精度與誤差

在試驗(yàn)車門上隨機(jī)選取6點(diǎn)，使用力錘敲擊，定位結(jié)果與誤差如表2所示。由表2可見(jiàn)，6次試驗(yàn)定位誤差均不超過(guò)3.16 cm，定位平均誤差為2.01 cm，滿足工程應(yīng)用要求。

表2 時(shí)間反轉(zhuǎn)法定位結(jié)果與誤差

車門Rattle異響噪聲源定位試驗(yàn)結(jié)果表明，本文所建立的基于TR的鈑金件Rattle異響噪聲源定位方法，具備較強(qiáng)的工程實(shí)際應(yīng)用價(jià)值。

4 結(jié)論

為滿足工程應(yīng)用中對(duì)汽車車門等鈑金件Rattle異響噪聲源客觀定位方法的需求，本文從Lamb波傳遞函數(shù)出發(fā)，建立了基于TR的異響噪聲源定位模型，通過(guò)仿真與試驗(yàn)探究了信號(hào)波包混疊的成因，確定了窄帶信號(hào)提取原則，提出了異響振動(dòng)信號(hào)發(fā)出時(shí)刻T0的疊加信號(hào)成像圖判定方法，并以某型汽車車門為例，驗(yàn)證了基于TR的汽車鈑金件Rattle異響噪聲源定位方法的可行性。主要結(jié)論如下：

（1）Rattle異響振動(dòng)信號(hào)分析結(jié)果表明，汽車鈑金件Rattle異響振動(dòng)信號(hào)為低頻Lamb波，其傳遞函數(shù)可簡(jiǎn)化為僅包含0階對(duì)稱模式的有限階表達(dá)式。

（2）仿真對(duì)比結(jié)果表明，端面反射是Lamb波傳播過(guò)程中波包混疊的原因，首達(dá)波包主要成分為振動(dòng)源直達(dá)信號(hào)，無(wú)端面反射波疊加。為避免波包混疊影響定位精度，應(yīng)截取首達(dá)波包，并以首達(dá)波包能量峰值頻率為小波中心頻率，提取窄帶信號(hào)進(jìn)行TR處理。

（3）分析了Rattle異響振動(dòng)信號(hào)發(fā)出時(shí)刻T0對(duì)TR定位成像的影響規(guī)律，提出了T0時(shí)刻的信號(hào)成像圖判定法。當(dāng)所取時(shí)窗長(zhǎng)度恰好等于tn-T0時(shí)，疊加信號(hào)模值較大區(qū)域連續(xù)且面積最小；而所取時(shí)窗長(zhǎng)度大于或小于tn-T0時(shí)，疊加信號(hào)模值較大區(qū)域離散且面積較大。通過(guò)對(duì)比不同T0時(shí)刻的Rattle信號(hào)TR定位成像圖，即可確定異響振動(dòng)信號(hào)發(fā)出時(shí)刻T0，并定位異響噪聲源位置。

（4）車門Rattle異響噪聲源定位試驗(yàn)結(jié)果表明，本文所建立的基于TR的汽車鈑金件異響噪聲定位新方法，定位最大誤差為3.16 cm，平均誤差為2.01 cm，定位精度滿足工程應(yīng)用要求，具有較強(qiáng)的工程實(shí)用價(jià)值。