紀 琳,黃震宇
(1.山東大學機械工程學院,濟南 250061;2.上海交通大學電子信息與電氣工程學院,上海 200240)
連接界面不定性對于加筋板結(jié)構振動能量傳遞特性的影響分析
紀 琳1,黃震宇2
(1.山東大學機械工程學院,濟南 250061;2.上海交通大學電子信息與電氣工程學院,上海 200240)
加筋板結(jié)構有著極為廣泛的工程應用,準確預測振動能量(結(jié)構聲)在加筋板中的傳遞特性是對于工程結(jié)構進行有效、合理聲學設計的一個極為重要的環(huán)節(jié)。本文以板-筋-板組合系統(tǒng)為研究對象,將加強筋視為板件結(jié)構的中間連接件,通過仿真分析計算,研究了加筋板的參數(shù)不確定性對筋板組合系統(tǒng)的功率流傳遞特性影響。考慮到加強筋與所連接板件通常具有較大的動態(tài)特性差異,仿真模擬以中頻混合分析法(Hybrid FE-SEAmethod)為主。研究表明:在長短波并存的中頻域,加筋板中的結(jié)構聲傳遞特性受加強筋剛度和邊界條件變化的影響較大,而受其材料屬性和連接自由度數(shù)量變化的影響較小。
振動傳遞;中頻混合模型法;統(tǒng)計能量法;參數(shù)不確定性;加筋板結(jié)構
加筋板結(jié)構是構成飛機、航天器、汽車、火車、船舶等大型交通運載工具的重要結(jié)構類型,準確預測加筋板結(jié)構的振動噪聲能量傳遞特性對于這些運載工具的艙室聲學設計有著非常重要的工程意義[1]。
加筋板結(jié)構的動態(tài)特性表現(xiàn)在結(jié)構本身的振動響應及振動能量在加筋板各子結(jié)構間的傳遞[2-3]。加筋板結(jié)構的理論建模方法,一般可以按照結(jié)構尺寸與結(jié)構波波長的相對長短關系分為低頻、中頻和高頻三個頻域的分析方法[4](注:與可測頻率的高低可能不一致)。其中,低頻域是指結(jié)構波波長遠大于結(jié)構尺寸的情況,高頻域指結(jié)構波波長遠小于結(jié)構尺寸的情況,而中頻域是指結(jié)構中既包含有長波子結(jié)構,又包含有短波子結(jié)構的情況。對于低頻域,分析方法以傳統(tǒng)的有限元、邊界元法為主,對于高頻域,以統(tǒng)計能量法分析為主,而對于中頻域則以混合模型法為主[5],即:將長波子系統(tǒng)用有限元模型構建,而對于短波子系統(tǒng)用統(tǒng)計能量法模擬[6]。在過去幾十年中,針對加筋板結(jié)構振聲響應的研究取得了長足的發(fā)展,尤其是對于加強筋與板之間存在較大動態(tài)特性反差的中頻域,取得了一系列里程碑式的研究成果[7]。
然而,在現(xiàn)有中頻研究方法和理論模型中,對于加強筋的動態(tài)模擬多以確定性分析方法為主[5-7],即:假定加強筋的材料和尺寸、連接邊界條件均為確定已知。顯然,這一假定條件與實際工程中加強筋自身必然具有一定程度的參數(shù)不確定性這一事實不相符合。
為了揭示連接邊界的參數(shù)不確定性對于加筋板整體結(jié)構振動能量傳遞特性的影響規(guī)律,本文選取具有典型中頻振動特征的加筋板結(jié)構,針對加強筋的結(jié)構屬性、材料屬性、連接邊界自由度數(shù)量及邊界條件的不確定性對于振動能量的傳遞特性的影響進行了仿真分析研究,分別將加強筋的結(jié)構尺寸(如長度、高度等)和材料屬性(如密度、阻尼損耗因子)視為圍繞其名義參數(shù)值上下波動的隨機參數(shù),同時也探討了當筋板之間的連接自由度數(shù)目存在某些不確定性以及加強筋自身的邊界條件存在一定的不確定性條件下的影響特性。
由于實際工程中,加強筋一般可視為薄板件間的中間連接件;同時,考慮到加筋板結(jié)構通常存在著加強筋的波長遠大于板件中的波長的關系,加筋板結(jié)構中的振動響應將在相當大的頻域內(nèi)呈現(xiàn)顯著的中頻振動特征,因此,本項研究中,對于加筋板振動能量傳遞分析主要是針對中頻域,所用建模方法為當前中頻振動領域影響力最為廣泛的Hybrid FE-SEA方法[6-8]。
關于Hybrid FE-SEA的建模原理和方法的詳細介紹請參見文獻[7]。本文中,為了便于讀者理解,作者先就一個簡單的梁-板組合結(jié)構(如圖1所示,其中f表示系統(tǒng)所受到的外部激振力向量,點力和/或分布力均可),對于Hybrid FE-SEA理論的振動能量傳遞計算過程簡述如下。

圖1 簡單加筋板模型Fig.1 A simple rib-stiffened platemodel
首先,假定板具有足夠的隨機不定性[6],那么,對于板可建立基于統(tǒng)計能量法(SEA)理論的能量平衡方程:

式中,左側(cè)Ep表示板的統(tǒng)計能量響應,ηp和ηb,p則分別表示板的阻尼損耗因子和梁板間的能量傳遞系數(shù)(power transmission coefficient);而式(1)右側(cè)的則表示在外部激振力矢量f的作用下,由梁輸入到板的功率流,分別可表達為[6]:



式(5)中,E[·]表示群體平均。
將式(4)~(5)代入到式(2)~(3),進而代入到式(1),就能夠得出板的統(tǒng)計能量響應Ep,再根據(jù)統(tǒng)計能量法理論可得出由外部激振力傳遞到板中的平均功率流為:

當加筋板結(jié)構中包含多個板子結(jié)構時(如圖2所示的板-梁-板組合結(jié)構),可以將式(1)~(4)拓展為矩陣形式并進行相應的求解,從而得出由外部激振力f分別傳遞到各個板子結(jié)構中的平均功率流及統(tǒng)計能量響應情況[7-8]。
圖2所示的加筋板結(jié)構為由板1和板2通過梁(加強筋)聯(lián)接而構成,其中,假定梁的邊界條件為兩端自由,而兩個板的邊界條件為:與加強筋平行的邊界為簡支,而另外兩對邊為自由邊界條件,各子結(jié)構的材料和尺寸參數(shù)分別示于表1和表2。

圖2 板-筋-板組合模型Fig.2 A plate-beam-plate couplingmodel
當板1在任意位置處受到一個點簡諧激振力的作用時,振動能量輸入到板1中,進而通過中間連接件(梁)傳遞到板2中。
圖3所示為板1和板2通過中間加強筋相聯(lián)和無加強筋直接相聯(lián)的平均功率流傳遞的對比情況。這里所采用的平均功率流的計算方法如下:首先,激勵點位于板1上的位置用隨機選取的方式選用了7個不同的位置;其次,對于每一個激勵位置,采用子系統(tǒng)模態(tài)法[7]得出時間平均傳遞功率流的精確計算值;最后,對于所有激勵位置處所對應的時間平均傳遞功率流計算結(jié)果進行空間平均,得出結(jié)構在統(tǒng)計能量法(SEA)意義上的平均傳遞功率流。
這里應該指出:在子系統(tǒng)模態(tài)法的計算過程中,梁和板之間的線連接形式采用間隔為1/4波長(板的彎曲波)的均布離散點連接形式來模擬,顯然,由于板的彎曲波長隨著頻率的升高而減小,對于相同尺寸的線連接來說,所需連接點的數(shù)量將隨著頻率的提高而增大,計算耗時也逐漸增加。

表1 材料特性參數(shù)(鋼)Tab.1 Parameters ofmaterial properties

表2 尺寸特性參數(shù)(板1-梁-板2)Tab.2 Parameters of geometry properties

圖3 板1與板2間的平均傳遞功率流Fig.3 Power transmission between plate 1 and plate 2
由圖3可見:中間加強筋的存在能夠顯著降低板結(jié)構間的能量傳遞,顯然,梁的動態(tài)特性對于系統(tǒng)的功率流傳遞特性能夠產(chǎn)生很大的影響,這就使得研究加強筋的參數(shù)不確定性對于加筋板結(jié)構功率流傳遞的影響具有十分重要的意義。
對于圖2所示的加筋板結(jié)構,由其參數(shù)特性(表1、表2)可知:筋的模態(tài)密度遠低于板的模態(tài)密度,這就意味著加筋板結(jié)構的振動特性將以中頻振動為主,因此,本節(jié)對于結(jié)構中功率流傳遞特性的計算中采用了本文第一節(jié)所述的Hybrid FE-SEA理論。
圖4中對比了梁的長度和高度與其名義尺寸相比分別變化5%時對于結(jié)構中傳遞功率流的影響情況。可見:隨著頻率的增大,即使是梁的結(jié)構尺寸的微小變化亦能夠?qū)τ趥鬟f功率流的峰值位置產(chǎn)生較大的影響,這就意味著:連接界面剛度的微小變化對于加筋板振動能量的傳遞的影響程度隨著頻率的升高而增大;同時,圖4亦表明:增大梁的剛度,能夠在一定程度上減少加筋板結(jié)構振動能量的傳遞。
圖5對比了梁的邊界條件分別為兩端自由和兩端簡支條件下,結(jié)構中振動能量傳遞的計算結(jié)果。由圖5可見:加強筋自身的邊界條件的改變能夠?qū)τ诮Y(jié)構的振動能量傳遞的峰值和谷值的位置產(chǎn)生顯著影響,這就使得對于一個特定的頻域而言,減小加筋板結(jié)構的結(jié)構聲傳遞性能有可能通過改變加強筋的邊界條件來實現(xiàn),如:圖5中,對于頻率范圍(1 200,1 600)Hz和(1 900,2 400)Hz而言,梁的邊界條件為兩端簡支時的傳遞功率流顯著低于兩端自由時的傳遞功率流。

圖4 梁的結(jié)構尺寸參數(shù)對于結(jié)構傳遞功率流的影響對比Fig.4 Comparison of transmitted power when the dimensional parameters of the beam are varied

圖5 梁的邊界條件對于結(jié)構傳遞功率流的影響對比Fig.5 Comparison of transmitted power when the boundary conditions of the beam are varied
然而值得注意的是,圖5亦同時表明:對于寬頻激勵,加強筋邊界條件的變化對于結(jié)構的功率流傳遞的平均水平的影響隨著頻率的增大而減小,也就意味著:連接邊界處兩端的具體邊界條件對于板1與板2之間的有效耦合因子的影響隨著頻率的增加而減小。
圖6對比了當作為連接件的梁的質(zhì)量密度在-5%~+5%變化時,對于板1與板2間的傳遞功率流的影響。由圖可見:連接邊界處質(zhì)量密度的變化對于加筋板振動能量傳遞的影響程度弱于連接邊界剛度變化的影響。
圖7對比了加強筋的阻尼特性的變化對于傳遞功率流的影響,表明:即使是在梁的材料阻尼值呈現(xiàn)20%的大幅度變化時,其對于加筋板結(jié)構的振動能量傳遞的影響也幾乎可以忽略不計。這就意味著:增大連接界面處的阻尼值,并不能夠有效降低板之間的結(jié)構聲傳遞水平。
圖8對比了連接界面自由度數(shù)量(即模擬梁板間連續(xù)線連接所采用的連接點數(shù)量)對于加筋板結(jié)構聲傳遞預測結(jié)果的影響情況,這里,梁與板之間的實際連續(xù)線連接用一組均勻分布的離散連接點模擬,而連接點數(shù)量的變化通過改變相鄰連接點之間的間距來實現(xiàn)。圖8中連接點的間距分別采用了板的彎曲波長的一半、三分之一和四分之一的情況,由圖可見:連接邊界自由度數(shù)量的變化對于加筋板中的結(jié)構聲傳遞特性影響主要體現(xiàn)于谷值域,而對于峰值處的影響相對較小。

圖6 梁的質(zhì)量密度對于結(jié)構傳遞功率流的影響對比Fig.6 Comparison of transmitted power when the mass density of the beam material is varied

圖7 梁的阻尼損耗因子對于結(jié)構傳遞功率流的影響對比Fig.7 Comparison of transmitted power when the damping loss factor of the beam is varied

圖8 連接界面自由度數(shù)量對于結(jié)構傳遞功率流的影響對比Fig.8 Comparison of transmitted power when the numbers of interface degrees of freedom(DOFs)are varied
加筋板作為若干工程結(jié)構中的重要組成部分,其結(jié)構聲傳遞特性對于實現(xiàn)工程結(jié)構的減振降噪設計有著非常重要的影響作用。由于加強筋通常可視為不同板件間的連接邊界,且加強筋與所連接板件通常存在著顯著的動態(tài)特性差異(如加強筋的振動主要表現(xiàn)為長波變形,而板件振動以短波局部變形為主),因此,本文用加筋板的參數(shù)變化來模擬連接邊界處的參數(shù)不確定,并采用中頻混合模型法(Hybrid FE-SEA Method)作為主要分析方法,通過一個板-筋-板的簡單加筋板模型,對于連接邊界的參數(shù)不確定性的結(jié)構聲傳遞影響規(guī)律進行了仿真計算分析。
研究結(jié)果表明:在長、短波并存的中高頻域,連接邊界對于加筋板的振動能量傳遞特性(即:結(jié)構聲傳遞)能夠產(chǎn)生較大的影響,尤為表現(xiàn)為連接界面處剛度和邊界條件的不確定性,而連接界面處的材料特性(如質(zhì)量密度和阻尼特性)等參數(shù)不確定性對于加筋板中結(jié)構聲傳遞的影響較小。
以上研究結(jié)論將對采用商用軟件(如VA-One[6]等)建立諸如飛機、高鐵、潛艇等的艙室筋板結(jié)構的中頻聲振響應預測分析模型時,如何引入連接邊界的參數(shù)不確定性的影響提供有效借鑒。
[1]Galbrun L.Vibration transmission through plate/beam structures typical of lightweight buildings:Applicability and limitations of fundamental theories[J].Applied Acoustics,2010,71:587-596.
[2]喬百杰,趙彤,陳雪峰.功率流測量以及振動能量參數(shù)估計試驗研究[J].振動與沖擊,2014,33(7):194-198.
QIAO Bai-jie,ZHAO Tong,CHEN Xue-feng.Tests for power flow measurement and vibrational energy estimation[J].Journal of Vibration and Shock,2014,33(7):194-198.
[3]薛開,王久法,李秋紅,等.中厚矩形板的振動功率流特性分析[J].振動與沖擊,2013,32(21):178-181.
XUE Kai,WANG Jiu-fa,LIQiu-hong,etal.Vibration power flow analysis of moderately thick rectangular plates[J].Journal of Vibration and Shock,2013,32(21):178-181.
[4]Soize C.A model and numerical method in the medium frequency range for vibroacoustic predictions using the theory of structural fuzzy[J].Journal of the Acoustical Society of America,1993,94:849-865.
[5]DesmetW.Mid-frequency vibro-acoustic modelling:challenges and potential solutions[C]//Proceedings of International Conference on Noise and Vibration Engineering.Leuven:Leuven University,2002:835-862.
[6]Shorter P J,Langley R S.Vibro-acoustic analysis of complex systems[J].Journal of Sound and Vibration,2005,288:669 -700.
[7]紀琳.中頻分析方法[M].北京:機械工業(yè)出版社,2013.
[8]PAM VA-One用戶手冊.ESI集團,2013.
Effects of ribs'parameter uncertainty on vibration energy transm ission of Rib-stiffened p late structures
JILin1,HUANG Zhen-yu2
(1.School of Mechanical Engineering,Shandong University,Jinan 250061,China;2.School of Electronic Information and Electrical Engineering,Shanghai Jiao Tong University,Shanghai200240,China)
It is of great importance to predict the vibration energy transmission(structure-borne sound transmission)within rib-stiffened plate structures due to their large applications in practical engineering.Since ribs usually were taken as junctions between plates,here parameter uncertainty effects of stiffening beams(ribs)on the power transmission between coupled thin plates were studied based on a plate-beam-plate model.Considering the significant dynamic features'mismatches between the ribs and plates,the hybrid FE-SEA method was adopted for simulations.Numerical investigations showed that the variations of stiffness and boundaries of ribs have relatively bigger effects on the vibration energy transmission of ribs-stiffened plate structures than the variations of material properties and interface coupling DOFs do.
vibration transmission;mid-frequency hybrid FE-SEA theory;statistical energy analysis(SEA);parameter uncertainty;rib-stiffened plate structures
TB53;TH212
A
10.13465/j.cnki.jvs.2015.21.032
國家國際科技合作專項資助項目(2012DFG01790);國家自然科學基金(51175300)
2014-05-08 修改稿收到日期:2014-06-17
紀琳女,博士,教授,1967年9月生
黃震宇男,博士,副教授,1975年1月生