劉豆豆, 賀西平, 李家星

(陜西師范大學 物理學與信息技術學院， 陜西 西安710119)

?

縱振激振面積對條紋振動矩形板聲場指向性的影響

劉豆豆, 賀西平*, 李家星

(陜西師范大學 物理學與信息技術學院， 陜西 西安710119)

摘要：在大功率超聲應用中，縱振激振換能器與階梯板間有一定接觸面積，導致振動節(jié)線與凸臺邊緣不重合。利用數(shù)值計算和有限元分析相結(jié)合的方法，研究了(8,0)條紋振動模式下激振源面積大小對矩形板聲場指向性的影響。結(jié)果表明，對尺寸為180 mm×120 mm×5 mm(長×寬×厚)的矩形薄板，隨著激振源半徑增大，條紋振動節(jié)線逐漸彎曲，(8,0)模態(tài)共振頻率增大，指向性尖銳程度減小。當觀察平面垂直于矩形板表面時，沿y軸方向指向性更為尖銳。當半徑增大到12 mm時，副瓣已明顯出現(xiàn)，故對該矩形板，應選擇激振源半徑在8 mm以內(nèi)。

關鍵詞:有限元法；矩形板；激振源面積；節(jié)線；指向性

輻射體的聲場指向性在大功率超聲應用中是一個重要問題，研究者對此進行了大量研究。文獻[1]對楔形換能器遠場作了分析和實驗，理論上給出了介質(zhì)表面的瑞利波指向性圖案。何正耀等利用邊界元理論計算了水聲換能器陣的指向性，并進行了實驗驗證[2]。文獻[3-5]計算了圓柱殼類聲場的指向性。賀西平等計算了節(jié)圓與階梯凸臺位置重合時，階梯圓盤聲場的指向性[6]。文獻[7]測量了激勵面積較小時階梯型矩形板的聲場指向性。文獻[8]利用聲場迭加原理推導出了階梯圓盤輻射聲壓，并計算了面激振時(即激勵換能器與輻射盤之間有一定面積大小相聯(lián)接)節(jié)圓偏離階梯位置對聲場指向性的影響。文獻[9]計算了縱振激振面積不同對矩形薄板振動特性的影響。

縱振-彎曲輻射體振動和輻射系統(tǒng)由Gallego等提出，此類換能器綜合了縱振換能器的高效率和彎曲振動平圓盤大輻射面積的特點，能有效地與空氣介質(zhì)進行阻抗匹配并輻射強功率超聲[10]。平圓盤彎振時節(jié)線兩側(cè)相位相反，導致聲場指向性變差，文獻[11-13]在該平圓盤節(jié)線兩側(cè)反相位區(qū)域增加輻射介質(zhì)中聲波波長1/2的高度變?yōu)殡A梯狀。此類階梯狀輻射體類似于活塞振動盤，振動面上的各部分振動相位相同。激振換能器與該階梯圓盤接觸面積很小(可視為點接觸)，激振頻率為該階梯圓盤本征頻率時，階梯圓盤節(jié)線位置與階梯位置相重合，輻射聲場指向性尖銳[14-15]。在振動面內(nèi)，矩形板比圓形板振幅分布更加均勻，應力分布也較均勻，可以承受更大的機械強度。同理，可以在矩形板本征條紋彎曲振動模式的基礎上，改進成為階梯板?？v振換能器在其本征頻率點接觸激振該板中心，換能器的振動幾乎不影響階梯板彎曲振動特性，此時振動節(jié)線將與階梯邊緣位置重合，輻射聲場具有尖銳指向性。實際應用中，尤其是大功率狀態(tài)下，激振換能器一般與輻射盤或板肯定有一定的接觸面積， 即為面接觸激振。面接觸激振時，就會導致板的節(jié)線與預設的階梯位置不重合。

本文以矩形板為例，研究在自由邊界條件下， (8,0)模態(tài)節(jié)線和預設階梯位置不重合時的輻射聲場指向性。將瑞利積分公式離散化，結(jié)合有限元方法計算矩形板的振動模態(tài)，提取其各有限單元上的振動位移，再利用聲場迭加原理，計算研究激振源面積變化對矩形板輻射聲場指向性的影響。發(fā)現(xiàn)面激振的矩形板振動節(jié)線變形而不再保持為直線，且接觸面積越大，節(jié)線彎曲也越明顯[16]。計算結(jié)果為分析面激振下階梯形矩形板節(jié)線偏離預設階梯位置對指向性的影響提供了參考。

1矩形板的振動特性

取材料為45#鋼，尺寸為(180×120×5)mm(長×寬×厚)的矩形板,其楊氏模量E=2.16×1011N2/m，泊松比δ=0.28，密度ρ=7 800 kg/m3?？諝庵新曀贋?40 m/s，對應(8,0)模態(tài)本征振動頻率為20 216 Hz，振動節(jié)線基本為直線。按照此時的直節(jié)線分布將矩形板改進成階梯形，若仍為點源激振，階梯形矩形板的振動節(jié)線將與階梯邊緣的位置相一致，但若面激振改進后的階梯形矩形板，則板的振動節(jié)線將彎曲變形不再與預設的階梯邊緣重合。

圖1為激振源半徑為r=16 mm時，矩形板彎曲振動振型圖，在該振型圖上，依照(8,0)本征模態(tài)直振動節(jié)線位置在圖1上標注，從左至右共八條，其中第1、2條，3、4條，5、6條以及7、8條直線之間為預設的加階梯位置，所加階梯的高度為介質(zhì)中輻射聲波波長的1/2。換能器面激振矩形板，條紋振動模式時，節(jié)線扭曲變形偏離原直節(jié)線位置，陰影部分表示不加階梯時矩形板面激振下相位為負的區(qū)域，空白部分表示相位為正的區(qū)域，陰影部分的邊緣線即為扭曲后的節(jié)線。也就是說，面激振時，振動節(jié)線不再保持為直線而呈現(xiàn)彎曲狀。從左至右共有8條，從小到大依次編號，原第1、2、7、8條節(jié)線中間部分已經(jīng)斷開，第3條、第6條節(jié)線彎曲程度較大，中間部分分別與第2條和第7條節(jié)線相交，第4條、第5條節(jié)線彎曲程度較大，整個矩形板振型復雜。

2聲場指向性計算方法

圖2為用于計算彎曲振動矩形板聲場指向性的幾何示意圖，假設該彎曲振動矩形板嵌在無限大障板上，彎振時會向z軸正向的無限大半空間內(nèi)輻射聲波，其中A為聲場中任意觀察點，直角坐標和球坐標分別為(x0,y0,z0)和(r,θ,φ)，ds為矩形薄板上一微小面元，ds=dxdy,r為觀察點A與坐標原點O之間距離，h為觀察點到面元距離。

由圖2可知,觀察點與面元ds間距為

(1)

通過瑞利積分公式，有面元ds在觀察點A處產(chǎn)生的聲壓為

(2)

其中u(x,y)表示面元處的振速幅值。將所有點源輻射聲波進行疊加可得

(3)

對于r?ω和l的遠場區(qū)域，不考慮矩形板自身尺寸，上式中振幅部分的h可近似用板中心到觀察點距離r來替代。結(jié)合遠場條件，近似表示為

(4)

觀察點A所處的方向角為

cosα=x0/r=sinθcosφ，

cosβ=y0/r=sinθsinφ。

(5)

將(4)、(5)兩式代入聲壓表達式(3)可得

[(cos(kxsinθcosφ)cos(kysinθsinφ)-

sin(kxsinθcosφ)sin(kysinθsinφ))+

j(cos(kxsinθcosφ)sin(kysinθsinφ)+

cos(kysinθsinφ)sin(kxsinθcosφ))]dxdy。

(6)

令

(7)

離散上述兩式得

A點的聲壓幅值表達式為

(9)

其中：N為離散后組成振動面的單元總數(shù)，Δsi為第i個單元的面積，vi表示第i個單元的中心點速度，xi、yi為第i個單元的橫坐標和縱坐標，這些量可以從有限元數(shù)值計算中獲取。

在有限元模態(tài)計算模塊，提取板輻射表面各單元撓度Li。假設矩形板為小撓度簡諧振動，則板的橫向位移可表示為

η=Liejωt。

(10)

由此可得矩形板表面各單元振速為

(11)

將(11)式帶入(8)式計算彎曲振動矩形板輻射聲場指向性表達式為

(12)

結(jié)合(8)、(12)式，提取振動位移后，對位移值進行處理，如圖1中所示，將預設加階梯部分的負相位位移值全部變?yōu)檎?，正相位位移值全部變?yōu)樨撝?，未加階梯部分位移值相位不變，將處理后的數(shù)據(jù)匯總，編制程序計算得到面激振情況下矩形板彎曲振動輻射聲場指向性。

3計算結(jié)果及分析

3.1縱振激振源點激勵矩形板

利用有限元方法得到本征頻率下彎曲振動矩形板(8,0)模態(tài)振型圖，如圖3a所示，矩形板各彎振節(jié)線基本為直線，陰影部分兩側(cè)相鄰節(jié)線位置即為預設階梯的邊緣位置。提取矩形板前表面各面元位移值代入計算指向性程序中，得到觀察平面垂直于矩形板表面并分別沿x軸和y軸方向的指向性如圖3b、3c所示。

可以看出，不加階梯情況下，矩形板(8,0)模態(tài)聲場指向性較差，如3b所示，觀察平面垂直于矩形板表面并沿x軸方向時，指向性圖案中在θ=0方向上沒有主瓣，即軸線方向上沒有聲能量輻射。3c中觀察平面沿y軸時，指向性圖案中在θ=0方向上出現(xiàn)主瓣，但在其他方向有較大副瓣，在指向性要求高的場合無法使用。若以陰影部分的邊緣直線為界給矩形板增加高度為空氣中聲波波長1/2的階梯成為活塞板，則圖3a中矩形板陰影部分的負相位位移反相，整個輻射表面振動相位基本都為正值，本文對預設加階梯部分的位移值進行反向處理后計算得加階梯矩形板的指向性如圖4所示。

圖3彎曲振動矩形板的振型和聲場指向性

Fig.3The directivity diagram of a rectangular plate in flexural vibration

圖4彎曲振動階梯矩形板聲場指向性

Fig.4The directivity diagram of a rectangular

plate in flexural vibration

圖4a中，當觀察平面沿x軸方向時，指向性圖案在θ=0方向上有較尖銳的主瓣，在其他方向有較小副瓣；當觀察平面垂直于矩形板表面并沿y軸方向時，指向性圖案只在θ=0方向上也有尖銳的主瓣。上述方法計算得到階梯形矩形板的輻射聲場指向性。

3.2縱振激振源面激振矩形板

分別取縱振換能器激振源半徑為4、8、12和16 mm，建立面激勵振動模型。利用有限元計算，得到各振動模型的(8,0)模態(tài)振型圖，如圖5所示。

隨著激振源半徑的增加，(8,0)模態(tài)的共振頻率逐漸增大，節(jié)線彎曲程度逐漸增大。圖5a中，當激振源半徑r=4 mm時，矩形薄板(8,0)模態(tài)振動節(jié)線基本呈直線，在預設階梯位置加階梯后圖中陰影部分負相位位移基本上全都落在階梯上而變?yōu)檎?，只有小部分負相位位移落在階梯外。當激振半徑增大至r=8 mm時(圖5b)，第2、3、6、7條節(jié)線彎曲程度不明顯，基本為直線，第1、4、5、8條節(jié)線向外彎曲，彎曲程度較小,在預設階梯位置加階梯后落在階梯外的負相位位移增多。圖5c為r=12 mm時矩形板(8,0)模態(tài)振型圖，可以看出，節(jié)線彎曲程度較大，第1、3、4、5、6、8條節(jié)線向外彎曲，第2、7條節(jié)線向內(nèi)彎曲，負相位位移部分偏離加階梯位置嚴重。圖5d中，當激振半徑增大到r=16 mm時，振動節(jié)線彎曲程度更加明顯，第1、2、7、8條節(jié)線中間部分已經(jīng)斷開，第3、6條節(jié)線中間部分分別與第2、7條節(jié)線相交。隨著激振源半徑逐漸增大，條紋彎曲振動狀態(tài)下節(jié)線偏離預設階梯邊緣位置程度增大，將處理后的位移值匯總，帶入指向性程序，計算得到不同激振半徑下加階梯矩形板的指向性(如圖6、圖7所示)。

圖6不同激振半徑下觀察平面沿x軸指向性圖

Fig.6Directivity diagram with different excitation radius at observed surface along thexaxis

圖6為觀察平面垂直于矩形板并沿x軸，當激振源半徑較小時，指向性在θ=0方向上有尖銳的主瓣，其他方向上副瓣較小，隨著激振源半徑逐漸增大，副瓣越來越大，聲能量逐漸分散，矩形板指向性越來越差。圖7為觀察平面沿y軸方向，指向性尖銳程度較沿x軸時好，r≤8 mm時，觀察平面沿y軸方向，指向性只有主瓣，基本沒有副瓣，指向性尖銳。隨著激振源半徑增大，副瓣增多，指向性尖銳程度逐漸減小，當半徑增大到12 mm時，副瓣已十分明顯，聲能量分散，指向性尖銳程度較差。對該矩形板，應選擇縱振激振源半徑在8 mm以內(nèi)。

4結(jié)論

對于縱振換能器、變幅桿及矩形板組成的復合振動系統(tǒng)，本文計算了激振源變幅桿面積不同對矩形板振動特性的影響，得到如下結(jié)論：

(1)對于尺寸為(180×120×5) mm的矩形板，隨著激振源面積的增加，(8,0)模態(tài)的共振頻率逐漸增大。

(2)隨著縱振激振源面積的逐漸增加，矩形板條紋振動節(jié)線與本征振動時的節(jié)線相比，不再為直條紋狀，彎曲程度逐漸明顯。激振源半徑增大，節(jié)線偏離預設階梯邊緣位置程度增大，反相區(qū)域增多。

(3)面激振階梯形矩形板，觀察平面垂直于矩形板表面且沿x軸方向的指向性尖銳程度較y軸好，激振源半徑增大，尖銳程度減小，當半徑增大到12 mm時，副瓣較大。對該矩形板，應盡量選擇激振源半徑在8 mm以內(nèi)。這一結(jié)果為分析面激振下階梯形矩形板節(jié)線偏離預設階梯位置對指向性的影響提供了參考。

參考文獻:

[1] 嚴仁博.超聲楔形換能器的體波和瑞利表面波指向性圖案[J].物理學報,1974,23(6):419-428.

[2] 何正耀，馬遠良.水聲共形陣輻射指向性計算方法及其實驗驗證[J].聲學學報,2007,32(3):270-274.

[3] TETSURO O, BORIS A, DAVID A B. Broadband multimode baffled piezoelectric cylindrical shell transducers[J].The Journal of the Acoustical Society of America,2007,121(6):3465-3471.

[4] BORIS A,DAVID A B, COREY L B. Effects of coupledvibrations on the acoustical performance of underwater cylindrical shell transducers[J]. The Journal of the Acoustical Society of America,2007,122(6):3419-3427.

[5] 商德江，何祚鏞.加肋雙層圓柱殼振動聲輻射數(shù)值計算分析[J].聲學學報,2001,26(3):193-201.

[6] 賀西平，張頻，宋旭霞．單臺階同心圓盤輻射聲場的指向性計算[J].聲學技術，2009，28(5):673-677．

[7] GALLEGO-JUAREZ J A, RODRIGUEZ G, ACOSTA V, et al. Power ultrasonic transducers with extensive radiators for industrial processing[J]. Ultrasonics Sonochemistry, 2010, 17: 953-964.

[8] 蘭正康，賀西平，馬煥培,等.階梯圓盤節(jié)圓偏離對輻射聲場指向性的影響[J].振動與沖擊,2013,32(12):105-115.

[9] HE X P,YAO J,ZHANG H P,et al.Effect of excitation area of longitudinal transducer on the flexural vibration characteristics of a rectangular plate in stripe mode[J]. Ultrasonics,2015,58:104-105.

[10] CAMPOS-POZUELO C, LAVIE A, DUBUS B G.Numerical study of air-borne acoustic field of stepped-plate high-power ultrasonic transducer[J]. Acta Acustica United With Acustica,1998,84:1042-1046.

[11] 賀西平.彎曲振動階梯圓盤輻射阻抗的計算方法[J].物理學報，2010,59(5):3290-3293.

[12] 賀西平，宋旭霞，李偉.階梯圓盤的設計及聲參數(shù)計算[J].聲學學報,2010,35(5): 502-507.

[13] XIANG Y,ZHANG L.Free vibration analysis of stepped circular midline plates[J].Joural of Sound and Vibration,2005,280: 633-655.

[14] JE Y，LEE H，PARK J，et al． A stepped plate bi-frequency source for generating a difference frequency sound with aparametric array[J]．The Journal of Acoustical society of America，2010，127(6): 3494-3502．

[15] GALLEGO-JUAREZ J A，RODRIGUEZ G，RIERA-FRANCO DE S E，et al.Recent developments in vibrating-plate microsonic transducers[J]． Ultrasonics，2002，40(1/8):889-893．

[16] 姚景，賀西平.激振面積對彎振矩形板振動特性的影響[J].陜西師范大學學報(自然科學版)2014,42(2):46-49.

〔責任編輯 李博〕

Effect of longitudinal vibration excitation source area on radiated acoustic field directivity pattern of stripe vibration mode rectangular plate

LIU Doudou, HE Xiping*, LI Jiaxing

(School of Physics and Information Technology,Shaanxi Normal University, Xi′an 710119, Shaanxi, China)

Abstract:The connection between the rectangular plat and the longitudinal vibration transducer always has a certain contact area in high-power applications.Thus the deformed nodal line of stripe vibration mode don′t coincide with the expected step edge positions.Using combination of the method of numerical calculation and finite element, the influence of the excitation area on the directivity pattern in (8,0) stripe mode flexural vibration is investigated. The result shows that, for the 180×120×5 mm (length×width×thickness) rectangular plate, the nodal lines deform gradually with the increasing of excitation area.The (8,0) modal resonance frequency increases.The directivity pattern becomes less sharply and the acoustic field directivity pattern along theyaxis is sharper than that along thexaxis.When the radius of the excitation area increases to 12 mm，the side lobe is obvious.The radius of the excitation area should not be greater than 8mm for the stepped plate.Keywords: finite element method; rectangular plate; excitation area; nodal line; directivity patternPACS: 43.35.+d

文章編號：1672-4291(2016)03-0038-05

doi:10.15983/j.cnki.jsnu.2016.03.232

收稿日期：2015-03-12

基金項目：國家自然科學基金(11374201)

*通信作者：賀西平，男，教授，博士生導師。E-mail:hexiping@snnu.edu.cn

中圖分類號:O426.1

文獻標志碼:A