• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    結(jié)構(gòu)-聲強(qiáng)耦合腔振聲響應(yīng)預(yù)報(bào)研究

    2022-03-18 00:46:20廖金龍朱海潮侯九霄
    振動(dòng)與沖擊 2022年5期
    關(guān)鍵詞:均方聲壓聲場

    廖金龍, 朱海潮, 侯九霄

    (海軍工程大學(xué) 振動(dòng)與噪聲研究所,船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室,武漢 430033)

    由彈性結(jié)構(gòu)和封閉腔組成的結(jié)構(gòu)-聲耦合系統(tǒng)廣泛地存在于生活和實(shí)際工程中,如起居臥室,船舶艙室,聲吶探測腔等,彈性結(jié)構(gòu)在外激勵(lì)下產(chǎn)生的聲輻射是艙室內(nèi)噪聲的主要來源。近年來,通過采用各種控制方案抑制艙內(nèi)聲場水平,提高艙室舒適性已成為聲學(xué)領(lǐng)域的研究熱點(diǎn)[1-5],因此如何準(zhǔn)確地預(yù)報(bào)外激勵(lì)源下結(jié)構(gòu)-聲耦合系統(tǒng)的響應(yīng)是進(jìn)行聲學(xué)設(shè)計(jì)和噪聲控制的關(guān)鍵,具有重要的工程應(yīng)用價(jià)值。

    多年來,國內(nèi)外學(xué)者圍繞此問題開展了大量的理論和實(shí)驗(yàn)研究,Pan[6]最早提出采用模態(tài)耦合理論建立結(jié)構(gòu)-聲耦合系統(tǒng)模型,實(shí)現(xiàn)了系統(tǒng)響應(yīng)的預(yù)測。Kim等[7]以模態(tài)耦合理論為基礎(chǔ),從阻抗和導(dǎo)納角度進(jìn)行了響應(yīng)分析,并通過實(shí)驗(yàn)驗(yàn)證了預(yù)測的準(zhǔn)確性。Geng等[8]利用有限元工具獲取復(fù)雜封閉空間的聲模態(tài),結(jié)合模態(tài)耦合理論研究了復(fù)雜封閉空間耦合系統(tǒng)在外力激勵(lì)和聲激勵(lì)作用下響應(yīng)和有源消聲問題,理論上可適用于任意封閉空間,在分析頻率較高時(shí),由于網(wǎng)格質(zhì)量的要求,會(huì)耗費(fèi)大量的計(jì)算資源。張肅等[9]采用小波迦遼金法分析了復(fù)雜結(jié)構(gòu)-聲耦合系統(tǒng)的響應(yīng),并通過實(shí)驗(yàn)驗(yàn)證了方法的準(zhǔn)確性。上述有關(guān)結(jié)構(gòu)-聲耦合系統(tǒng)響應(yīng)預(yù)報(bào)的研究主要適用于空腔深度較深、腔內(nèi)介質(zhì)為非致密的弱耦合情形,當(dāng)空腔深度較淺、腔內(nèi)介質(zhì)具有較大的質(zhì)量密度時(shí),結(jié)構(gòu)與聲場形成強(qiáng)耦合,由于采用剛性壁面聲模態(tài)函數(shù)和真空中結(jié)構(gòu)模態(tài)函數(shù)作為基函數(shù),在彈性邊界處速度不連續(xù),導(dǎo)致不能準(zhǔn)確地預(yù)報(bào)系統(tǒng)響應(yīng)。為此,Du等[10-11]通過引入輔助函數(shù)的方式對(duì)聲場聲壓和結(jié)構(gòu)位移傅里葉級(jí)數(shù)進(jìn)行改進(jìn),使其滿足邊界連續(xù)性條件,實(shí)現(xiàn)了結(jié)構(gòu)-聲強(qiáng)耦合系統(tǒng)響應(yīng)預(yù)報(bào),但是輔助函數(shù)參數(shù)的引入增加了公式復(fù)雜性和系統(tǒng)矩陣的尺寸。陳躍華等[12-13]在此基礎(chǔ)上,采用Chebyshev多項(xiàng)式級(jí)數(shù)替代改進(jìn)傅里葉級(jí)數(shù)對(duì)耦合系統(tǒng)進(jìn)行了建模,形式更加簡潔,不需要借助有限元工具就可以獲取耦合系統(tǒng)的固有頻率和模態(tài),但是由于Chebyshev多項(xiàng)式是加權(quán)正交,在計(jì)算中包含大量高階積分項(xiàng),計(jì)算速度較慢。Kim等[14]提出耦合系統(tǒng)建模降階方法,實(shí)現(xiàn)了強(qiáng)耦合系統(tǒng)響應(yīng)的準(zhǔn)確預(yù)報(bào)。

    在現(xiàn)有的文獻(xiàn)中,對(duì)于強(qiáng)耦合腔體的響應(yīng)預(yù)報(bào)研究并不多見,因此開展結(jié)構(gòu)-聲強(qiáng)耦合腔的振聲響應(yīng)分析具有重要意義。本文從能量角度出發(fā),考慮結(jié)構(gòu)-聲耦合以及邊界力激勵(lì)的作用,對(duì)結(jié)構(gòu)-聲強(qiáng)耦合腔響應(yīng)進(jìn)行了分析。建立了彈性板-充水矩形背腔耦合模型,首先對(duì)系統(tǒng)結(jié)構(gòu)和聲場及其滿足的邊界條件進(jìn)行了說明,根據(jù)能量原理用拉格朗日函數(shù)對(duì)耦合系統(tǒng)動(dòng)力學(xué)方程進(jìn)行了描述,然后用權(quán)函數(shù)為1的 Legendre多項(xiàng)式級(jí)數(shù)分別將結(jié)構(gòu)位移函數(shù)和聲場聲壓函數(shù)進(jìn)行展開,采用Rayleigh-Ritz法對(duì)動(dòng)力學(xué)方程求解得到了耦合系統(tǒng)的振聲響應(yīng),由于Legendre多項(xiàng)式滿足L2內(nèi)積正交性,使得方程中高重積分項(xiàng)簡化為低重積分或者直接乘積項(xiàng),大大地提升了計(jì)算效率。通過與文獻(xiàn)和數(shù)值仿真結(jié)果對(duì)比,驗(yàn)證了本文方法的正確性。

    1 理論推導(dǎo)

    1.1 耦合系統(tǒng)模型

    考慮如圖1的彈性板-矩形背腔耦合系統(tǒng),背腔體積為V,腔內(nèi)充滿聲速為cc,密度為ρc的重介質(zhì)流體,其邊界由彈性板和剛性壁面組成,彈性板厚度為h,面積為S,邊界為Γ,材料密度為ρp,在簡諧外力激勵(lì)Fp的作用下,彈性壁面發(fā)生彎曲振動(dòng)并向腔體內(nèi)部輻射噪聲。由于腔內(nèi)介質(zhì)為重流體(比如水介質(zhì)),彈性邊界與封閉腔形成強(qiáng)耦合系統(tǒng)。

    圖1 彈性板-矩形耦合腔

    假設(shè)腔內(nèi)介質(zhì)在受到激勵(lì)前處于靜止?fàn)顟B(tài),則聲壓分布p(x,y,z)滿足Helmholtz波動(dòng)方程及邊界條件

    (1)

    (2)

    式中:?為拉普拉斯算子;ω為激勵(lì)頻率;w(x,y)表示彈性板的位移分布;n為彈性板外法線方向。

    對(duì)于彈性板,在外激勵(lì)的作用下其位移振動(dòng)方程為

    (D?4-ρphω2)w(x,y)=p(x,y)-

    Fpδ(x-xext)δ(y-yext)

    (3)

    式中,D為結(jié)構(gòu)彎曲剛度。對(duì)于彈性壁面約束邊界條件可以用結(jié)構(gòu)彎曲振動(dòng)的彈簧剛度系數(shù)來描述,以邊界x=0為例,其邊界條件滿足

    (4)

    (5)

    式中:η為材料泊松比;ax0和Ax0分別為x=0處平移彈簧剛度和旋轉(zhuǎn)彈簧剛度,通過設(shè)置平移彈簧剛度為無限大,旋轉(zhuǎn)彈簧剛度為零來表示簡支邊界條件,設(shè)置兩者均為無限大來表示固支邊界條件。

    直接由Helmholtz波動(dòng)方程、振動(dòng)方程及其邊界條件求解聲壓分布和位移分布的理論解是困難的,本文從能量角度出發(fā),采用拉格朗日函數(shù)對(duì)系統(tǒng)動(dòng)力學(xué)方程進(jìn)行描述,該耦合系統(tǒng)分為彈性板系統(tǒng)和封閉腔系統(tǒng),其拉格朗日函數(shù)分別為

    Lp=Up-Tp-Wext-Wctp

    (6)

    Lc=Uc-Tc-Wptc

    (7)

    式中:Up為彈性結(jié)構(gòu)彎曲振動(dòng)時(shí)的總勢能;Tp為總動(dòng)能;Wext為外激勵(lì)對(duì)彈性結(jié)構(gòu)做的功;Wctp為彈性結(jié)構(gòu)與聲場交界處聲壓對(duì)結(jié)構(gòu)做的功;Uc為聲場的總勢能;Tc為聲場的總動(dòng)能;Wptc為彈性結(jié)構(gòu)與聲場交界處結(jié)構(gòu)振動(dòng)對(duì)聲場做的功,由牛頓第三定理可知Wctp=-Wptc。具體的表達(dá)式為

    (8)

    (9)

    Wext=Fpw(xext,yext)

    (10)

    (11)

    (12)

    (13)

    1.2 動(dòng)力學(xué)方程求解

    利用Rayleigh-Ritz法對(duì)耦合系統(tǒng)的拉格朗日函數(shù)方程進(jìn)行求解,在求解時(shí)需要選擇一組試函數(shù)對(duì)聲場聲壓和結(jié)構(gòu)位移函數(shù)進(jìn)行展開。通常,要求選擇的試函數(shù)在對(duì)聲場聲壓和結(jié)構(gòu)位移展開時(shí)能滿足任何可能的邊界條件??紤]試函數(shù)條件及聲模態(tài)函數(shù)和結(jié)構(gòu)模態(tài)函數(shù)的正交特性,選擇正交多項(xiàng)式簇中Legendre多項(xiàng)式對(duì)聲場聲壓函數(shù)和結(jié)構(gòu)位移函數(shù)進(jìn)行展開,由于Legendre多項(xiàng)式的定義區(qū)間為[-1,1],在對(duì)聲壓函數(shù)和位移函數(shù)進(jìn)行展開時(shí)需要進(jìn)行坐標(biāo)變換,假設(shè)笛卡爾坐標(biāo)系下x,y,z的定義區(qū)間分別為[0,Lx],[0,Ly],[0,Lz],則坐標(biāo)變換公式為

    α=(2x-Lx)/Lx

    β=(2y-Ly)/Ly

    γ=(2z-Lz)/Lz

    (14)

    聲壓函數(shù)和位移函數(shù)表達(dá)式分別為

    (15)

    (16)

    式中,Blmn和Cjq分別為(l,m,n)階聲壓展開系數(shù)和(j,q)階位移展開系數(shù);Pl(α)為Legendre多項(xiàng)式。Pl(α)的表達(dá)式為

    (17)

    應(yīng)當(dāng)說明的是,相比于其他正交多項(xiàng)式,除了有良好的完備性和遞推性之外,選擇使用Legendre多項(xiàng)式作為試函數(shù)還有一個(gè)明顯的優(yōu)點(diǎn)。Legendre多項(xiàng)式在定義區(qū)間內(nèi)L2內(nèi)積滿足權(quán)函數(shù)為1的正交條件

    (18)

    式中,δll′為克羅內(nèi)克函數(shù),當(dāng)l=l′時(shí)為1,否則為0。這將簡化耦合系統(tǒng)動(dòng)力學(xué)方程中矩陣元素的運(yùn)算。

    將聲壓函數(shù)展開式(15)和位移函數(shù)展開式(16)代入式(8)~(13),聯(lián)立式(6)和(7),按照Rayleigh-Ritz法,使彈性板系統(tǒng)和封閉腔系統(tǒng)的拉格朗日函數(shù)分別對(duì)聲壓展開系數(shù)和位移展開系數(shù)取極值

    (19)

    整理即可得到耦合系統(tǒng)動(dòng)力學(xué)方程

    (20)

    式中:B為聲壓展開系數(shù)組成的列向量;C為位移展開系數(shù)組成的列向量;W為作用在彈性結(jié)構(gòu)上的外激勵(lì)向量;Kp和Mp分別為彈性結(jié)構(gòu)的剛度矩陣和質(zhì)量矩陣;Kc和Mc分別為聲腔聲場的剛度矩陣和質(zhì)量矩陣;Kptc為彈性結(jié)構(gòu)對(duì)聲腔聲場的耦合作用矩陣。由于Legendre多項(xiàng)式L2內(nèi)積滿足正交性,各矩陣元素中高重積分項(xiàng)可簡化為低重積分或者直接乘積項(xiàng),具體表達(dá)式為

    (21)

    (22)

    (23)

    (24)

    (25)

    通過求解上述動(dòng)力學(xué)方程,可以得到聲壓展開系數(shù)和位移展開系數(shù),分別代入式(15)和(16),即可求得外激勵(lì)下結(jié)構(gòu)的位移振動(dòng)響應(yīng)和聲場的聲壓分布。由于結(jié)構(gòu)均方振速和聲場均方聲壓是從系統(tǒng)能量層面描述響應(yīng)水平的參量,在聲學(xué)設(shè)計(jì)和全局噪聲控制時(shí)更具有代表性,本文選擇其作為響應(yīng)分析目標(biāo)。由振速和位移的關(guān)系v(x,y)=jωw(x,y),結(jié)合Legendre多項(xiàng)式的正交性,均方振速和均方聲壓分別可表示為

    式中:Λp和Λc分別為結(jié)構(gòu)和聲場歸一化對(duì)角矩陣,由式(18)可得,其對(duì)角元素分別為

    (28)

    (29)

    2 理論驗(yàn)證

    2.1 收斂性和準(zhǔn)確性分析

    考慮一彈性板-矩形腔強(qiáng)耦合腔模型,幾何參數(shù)為:Lx×Ly×Lz=0.35 m×0.29 m×0.14 m,四邊簡支,腔內(nèi)充滿水以使板腔形成強(qiáng)耦合。彈性板結(jié)構(gòu)參數(shù)為:密度ρp=2 700 kg/m3,泊松比η=0.3,板厚h=0.015 m,楊氏模量E=72 GPa,結(jié)構(gòu)阻尼ζp=0.01;水介質(zhì)參數(shù)為:密度ρc=1 000 kg/m3,聲速cc=1 500 m/s,介質(zhì)阻尼ζc=0.01,在板表面(0.039,0.272)m位置設(shè)置一幅值為1 N的法向簡諧點(diǎn)力。在實(shí)際的計(jì)算過程中,分別截取前Nj,Nq階和前Nl,Nm,Nn階Legendre多項(xiàng)式對(duì)聲壓函數(shù)和位移函數(shù)進(jìn)行展開。為確保求解結(jié)果的準(zhǔn)確性,需要對(duì)截?cái)嚯A數(shù)的收斂性進(jìn)行分析。根據(jù)板腔幾何參數(shù),取Nj=Nq=Nl=Nm=2Nn=2N,以有限元計(jì)算結(jié)果為準(zhǔn),分別構(gòu)建均方振速和均方聲壓截?cái)嗾`差函數(shù)為

    (30)

    (31)

    式中:[fL,fH]為分析頻率范圍,由于0~350 Hz已經(jīng)包含足夠多的模態(tài)頻率,選擇該頻段進(jìn)行分析。截?cái)嗾`差隨截?cái)嚯A數(shù)的變化曲線如圖2所示,隨著截?cái)嚯A數(shù)增加,截?cái)嗾`差迅速減小,當(dāng)N=6時(shí),均方振速和均方聲壓截?cái)嗾`差分別為2.50%和1.62%??紤]計(jì)算資源,選擇截?cái)嚯A數(shù)為:Nj=Nq=Nl=Nm=2Nn=2N=12。均方振速和均方聲壓的理論計(jì)算結(jié)果與文獻(xiàn)[13]的方法以及有限元計(jì)算結(jié)果的對(duì)比如圖3所示,可以發(fā)現(xiàn),在分析頻段內(nèi),本文理論計(jì)算結(jié)果和文獻(xiàn)[13]以及有限元結(jié)果吻合良好。

    圖2 截?cái)嗾`差曲線

    (a) 均方聲壓

    在選擇COMSOL-Multiphysics中聲-殼耦合模塊進(jìn)行有限元分析時(shí),理論上要獲得可信的結(jié)果,最大網(wǎng)格尺寸應(yīng)小于分析頻段最短波長的1/6。本文中設(shè)置最大網(wǎng)格尺寸為0.008 m,對(duì)應(yīng)的上限分析頻率為31 250 Hz,對(duì)于本文分析的頻段是足夠的。為了進(jìn)一步分析網(wǎng)格尺寸對(duì)響應(yīng)預(yù)報(bào)的影響,考慮使用更小的網(wǎng)格尺寸0.006 m。圖4給出了兩種網(wǎng)格尺寸下系統(tǒng)均方振速響應(yīng)對(duì)比圖,可以看出兩者結(jié)果基本一致,驗(yàn)證了使用0.008 m的網(wǎng)格尺寸有限元計(jì)算結(jié)果作為對(duì)比的可靠性。

    圖4 不同網(wǎng)格尺寸下均方振速有限元結(jié)果

    2.2 計(jì)算效率分析

    文獻(xiàn)[13]提出使用切比雪夫多項(xiàng)式級(jí)數(shù)對(duì)結(jié)構(gòu)位移函數(shù)和聲場聲壓函數(shù)進(jìn)行展開,與之相比,本文使用Legendre多項(xiàng)式級(jí)數(shù)描述聲壓和位移,由于其滿足L2內(nèi)積正交性,簡化了計(jì)算。在同等算力的計(jì)算機(jī)條件下使用代碼實(shí)現(xiàn)時(shí),影響計(jì)算效率的主要是積分項(xiàng)的數(shù)量,下面從一重積分項(xiàng)的維數(shù)對(duì)兩種方法計(jì)算效率進(jìn)行分析,附錄給出了文獻(xiàn)[13]中剛度矩陣,質(zhì)量矩陣及耦合作用矩陣各元素的表達(dá)式。對(duì)于相同的分析模型,相同的截?cái)嚯A數(shù),根據(jù)式(21)~(25)和附錄可得本文方法和文獻(xiàn)[13]的方法的積分項(xiàng)的總維數(shù)分別為

    (32)

    (33)

    比較式(32)、(33)可以看出,對(duì)于相同的截?cái)?,在采用Legendre多項(xiàng)式級(jí)數(shù)描述聲壓和位移時(shí),由于積分項(xiàng)得到了簡化,極大地減小了積分項(xiàng)的總維數(shù),使得代碼處理更加容易,計(jì)算效率得到提升。

    3 振聲響應(yīng)參數(shù)分析

    彈性板-矩形背腔是典型的結(jié)構(gòu)-聲耦合系統(tǒng),雖然其結(jié)構(gòu)簡單,但是當(dāng)腔內(nèi)介質(zhì)為重流體時(shí),系統(tǒng)耦合程度增加,結(jié)構(gòu)和聲場的響應(yīng)趨于復(fù)雜,更難以預(yù)測,尋找系統(tǒng)參數(shù)的影響規(guī)律是進(jìn)行聲學(xué)設(shè)計(jì)的關(guān)鍵。選擇2.1節(jié)的彈性板-矩形腔模型,不改變激勵(lì)信息,分析彈性板邊界條件和背腔深度的變化對(duì)響應(yīng)的影響。

    3.1 彈性板邊界條件

    假設(shè)彈性板四邊具有一致的邊界條件,不同的邊界條件可通過設(shè)置平移彈簧剛度a和旋轉(zhuǎn)彈簧剛度A的不同組合來表示:平移彈簧剛度為無窮大,旋轉(zhuǎn)彈簧剛度為零代表簡支邊界,兩者均為無窮大代表固支邊界,這里計(jì)算時(shí)取1010表示無窮大。

    圖5給出了耦合系統(tǒng)隨邊界條件和頻率變化的振聲響應(yīng)等高線圖,保持平移彈簧剛度為無窮大,旋轉(zhuǎn)彈簧剛度逐漸由零增加到無窮大,代表簡支向固支約束邊界的過渡。一般來講,相比于真空條件下,存在耦合作用的系統(tǒng)峰值頻率會(huì)產(chǎn)生一定的偏差,這也反映了耦合作用對(duì)系統(tǒng)特性的影響,可用系統(tǒng)第一階非零峰值頻率的偏差量與真空條件下第一階非零峰值頻率的比值來表征耦合作用的強(qiáng)度,即:

    Δ=|fcoupled,1-fin-vacuo,1|/fin-vacuo,1

    (34)

    (a) 均方聲壓

    耦合強(qiáng)度系數(shù)隨邊界旋轉(zhuǎn)彈簧剛度的變化如圖6所示,結(jié)合圖5和圖6可以看出,聲場均方聲壓與結(jié)構(gòu)均方振速響應(yīng)峰值頻率的變化趨勢是一樣的。邊界約束可按旋轉(zhuǎn)彈簧剛度大小分為三個(gè)階段:當(dāng)A=100~102趨近于簡支,當(dāng)A=102~105邊界約束由簡支向固支過渡,當(dāng)A=105~1010趨近于固支。當(dāng)邊界約束處于趨近于簡支和趨近于固支階段時(shí),系統(tǒng)耦合強(qiáng)度系數(shù)基本不變,具體表現(xiàn)為耦合系統(tǒng)的響應(yīng)峰值幅值和頻率不隨邊界旋轉(zhuǎn)彈簧剛度的變化而變化;處于過渡階段時(shí),隨著邊界旋轉(zhuǎn)彈簧剛度的增大,邊界約束由簡支逐漸轉(zhuǎn)變?yōu)楣讨В瑫r(shí)系統(tǒng)耦合強(qiáng)度逐漸減弱,具體表現(xiàn)為響應(yīng)峰值頻率向高頻偏移,但由于耦合作用變化量較小,響應(yīng)峰值幅值基本不變。

    圖6 邊界旋轉(zhuǎn)彈簧剛度變化對(duì)耦合強(qiáng)度系數(shù)的影響

    3.2 背腔深度

    改變分析模型的背腔深度,保持模型其他參數(shù)不變,分析背腔深度對(duì)耦合系統(tǒng)響應(yīng)的影響。需要說明的是,為了保證求解精度,在利用Legendre多項(xiàng)式對(duì)聲壓函數(shù)進(jìn)行展開時(shí),腔深方向上的截?cái)嚯A數(shù)應(yīng)當(dāng)隨腔深變化進(jìn)行調(diào)整,在計(jì)算中,當(dāng)0

    圖7給出了耦合系統(tǒng)隨背腔深度和頻率變化的振聲響應(yīng)等高線圖。為了清楚的說明峰值響應(yīng)幅值的變化,圖8給出了前六階均方振速和均方聲壓峰值響應(yīng)幅值隨腔深的變化圖。同時(shí)根據(jù)式(34),可以得到耦合強(qiáng)度系數(shù)隨邊界背腔深度的變化,如圖9所示,結(jié)合圖7~圖9可以看出,隨著背腔深度的增加,耦合強(qiáng)度系數(shù)迅速減小,具體表現(xiàn)為耦合系統(tǒng)的響應(yīng)峰值頻率隨背腔深度的增加迅速向高頻偏移,由于耦合作用迅速減小,由板傳遞到腔內(nèi)的能量減少,均方聲壓峰值幅值逐漸降低,均方振速峰值幅值逐漸增加,當(dāng)腔深進(jìn)一步增大,系統(tǒng)的耦合強(qiáng)度系數(shù)不再變化,振聲響應(yīng)峰值頻率和幅值也不再隨腔深的變化而變化。

    圖7 不同背腔深度下的振聲響應(yīng)

    4 結(jié) 論

    本文提出一種基于能量原理的強(qiáng)耦合腔系統(tǒng)響應(yīng)計(jì)算方法,該方法利用Legendre多項(xiàng)式級(jí)數(shù)分別將聲場聲壓函數(shù)和結(jié)構(gòu)位移函數(shù)展開,結(jié)合Rayleigh-Ritz法對(duì)能量原理下拉格朗日函數(shù)形式的耦合系統(tǒng)動(dòng)力學(xué)方程進(jìn)行了求解,得到聲場聲壓響應(yīng)和結(jié)構(gòu)振速響應(yīng)。分析了彈性板邊界條件和背腔深度對(duì)振聲響應(yīng)的影響,結(jié)果表明:

    (1) 本方法能快速準(zhǔn)確的預(yù)測強(qiáng)耦合腔的振聲響應(yīng),由于Legendre多項(xiàng)式L2內(nèi)積滿足正交性,簡化了動(dòng)力學(xué)方程中高重積分項(xiàng)的計(jì)算,大大提升了計(jì)算效率。

    (2) 結(jié)構(gòu)邊界由強(qiáng)約束過渡到弱約束時(shí)會(huì)導(dǎo)致峰值頻率向低頻偏移,邊界條件的變化對(duì)峰值幅值的影響不明顯。

    (3) 隨著背腔深度的增大,耦合系統(tǒng)峰值頻率向高頻偏移,同時(shí)聲場均方聲壓峰值幅值降低,結(jié)構(gòu)均方振速峰值幅值增加,由結(jié)構(gòu)傳遞到聲場的能量減少。背腔深度進(jìn)一步增大,耦合系統(tǒng)的振聲響應(yīng)不再改變。

    附錄A

    (A.1)

    (A.2)

    (A.3)

    (A.4)

    (A.5)

    式中:Tl(α)為切比雪夫多項(xiàng)式,其表達(dá)式為:Tl(α)=cos(larccosα)

    猜你喜歡
    均方聲壓聲場
    一類隨機(jī)積分微分方程的均方漸近概周期解
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    Beidou, le système de navigation par satellite compatible et interopérable
    基于BIM的鐵路車站聲場仿真分析研究
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    探尋360°全聲場發(fā)聲門道
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    電測與儀表(2016年9期)2016-04-12 00:30:02
    基于抗差最小均方估計(jì)的輸電線路參數(shù)辨識(shí)
    基于隨機(jī)牽制控制的復(fù)雜網(wǎng)絡(luò)均方簇同步
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    中文欧美无线码| 97热精品久久久久久| 成年免费大片在线观看| 国产av一区二区精品久久 | freevideosex欧美| 日韩三级伦理在线观看| 日韩欧美精品免费久久| 亚洲精品国产av蜜桃| 国产伦精品一区二区三区四那| 国产av码专区亚洲av| 少妇被粗大猛烈的视频| 嫩草影院入口| av线在线观看网站| 国产有黄有色有爽视频| 日韩大片免费观看网站| 我要看黄色一级片免费的| 尤物成人国产欧美一区二区三区| 你懂的网址亚洲精品在线观看| 国产一区二区在线观看日韩| 久久 成人 亚洲| 日本黄色片子视频| 国产探花极品一区二区| 日韩一本色道免费dvd| 精品午夜福利在线看| 午夜老司机福利剧场| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 国产精品偷伦视频观看了| 国产视频首页在线观看| 国产精品伦人一区二区| 久久久精品免费免费高清| 国产免费福利视频在线观看| 久久久久视频综合| 国产片特级美女逼逼视频| 国产色爽女视频免费观看| 久久久久久久久久人人人人人人| 日韩伦理黄色片| 最近中文字幕高清免费大全6| 美女高潮的动态| 免费看不卡的av| 国产永久视频网站| 国产精品爽爽va在线观看网站| 国产成人aa在线观看| 青春草视频在线免费观看| 国产免费福利视频在线观看| a级毛片免费高清观看在线播放| 亚洲欧美中文字幕日韩二区| av专区在线播放| 边亲边吃奶的免费视频| 欧美精品国产亚洲| 最近2019中文字幕mv第一页| 国产av一区二区精品久久 | 最后的刺客免费高清国语| 久久精品国产a三级三级三级| h日本视频在线播放| 免费播放大片免费观看视频在线观看| 插阴视频在线观看视频| 国产 一区 欧美 日韩| 国产精品不卡视频一区二区| 久久精品熟女亚洲av麻豆精品| a 毛片基地| 欧美成人精品欧美一级黄| 国产91av在线免费观看| 欧美97在线视频| 黄色日韩在线| 国产精品福利在线免费观看| 国产成人免费观看mmmm| 日韩av免费高清视频| 18禁在线无遮挡免费观看视频| 99热全是精品| 亚洲国产毛片av蜜桃av| 三级经典国产精品| 1000部很黄的大片| 免费观看的影片在线观看| 天天躁夜夜躁狠狠久久av| 亚洲天堂av无毛| 一区二区三区精品91| 视频区图区小说| 久久久精品94久久精品| 亚洲精品国产色婷婷电影| 久久久国产一区二区| 最近中文字幕高清免费大全6| 热99国产精品久久久久久7| 毛片女人毛片| 国产乱来视频区| 国产精品蜜桃在线观看| 精品人妻一区二区三区麻豆| 久久久久久久久久成人| 只有这里有精品99| 一边亲一边摸免费视频| 91精品国产国语对白视频| 精品人妻熟女av久视频| 一区二区三区乱码不卡18| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 尤物成人国产欧美一区二区三区| 国产色婷婷99| 久久亚洲国产成人精品v| 在线天堂最新版资源| 久久久久久九九精品二区国产| 搡女人真爽免费视频火全软件| 免费久久久久久久精品成人欧美视频 | 在线免费观看不下载黄p国产| 久久久久久久国产电影| 在线观看免费视频网站a站| 亚洲经典国产精华液单| 1000部很黄的大片| 亚洲av成人精品一二三区| 国产亚洲午夜精品一区二区久久| 丰满少妇做爰视频| 老女人水多毛片| 欧美xxⅹ黑人| 国产综合精华液| 蜜桃亚洲精品一区二区三区| 全区人妻精品视频| 搡老乐熟女国产| 国产精品久久久久久精品古装| 欧美激情极品国产一区二区三区 | 最近中文字幕高清免费大全6| 婷婷色麻豆天堂久久| 大片电影免费在线观看免费| 天堂8中文在线网| 少妇人妻精品综合一区二区| 国产亚洲av片在线观看秒播厂| 中国美白少妇内射xxxbb| 你懂的网址亚洲精品在线观看| 日韩中字成人| 青青草视频在线视频观看| 三级国产精品片| h视频一区二区三区| 日本黄大片高清| 18禁裸乳无遮挡动漫免费视频| 精品一区二区三区视频在线| 狂野欧美激情性bbbbbb| 一级av片app| 成年女人在线观看亚洲视频| 欧美日韩综合久久久久久| 免费观看a级毛片全部| 国产亚洲最大av| 亚洲国产精品成人久久小说| 久久久久久久亚洲中文字幕| 人人妻人人澡人人爽人人夜夜| 日本av免费视频播放| 国产视频首页在线观看| 乱系列少妇在线播放| 美女中出高潮动态图| 免费观看在线日韩| 色吧在线观看| 午夜老司机福利剧场| 亚洲不卡免费看| 亚洲欧美日韩东京热| 午夜日本视频在线| 久久久久视频综合| 在线观看免费高清a一片| 国产美女午夜福利| 久久国产精品男人的天堂亚洲 | 久久国产精品大桥未久av | 七月丁香在线播放| 综合色丁香网| 人妻少妇偷人精品九色| 日日啪夜夜爽| 91精品一卡2卡3卡4卡| 在线看a的网站| 99久久精品热视频| 欧美一区二区亚洲| 国产精品久久久久久精品电影小说 | 色5月婷婷丁香| 国产乱人视频| 久久久国产一区二区| 99久久精品国产国产毛片| 亚洲图色成人| 人妻制服诱惑在线中文字幕| 汤姆久久久久久久影院中文字幕| 777米奇影视久久| 国产在视频线精品| 老熟女久久久| 黄色日韩在线| 日产精品乱码卡一卡2卡三| 十分钟在线观看高清视频www | 精华霜和精华液先用哪个| 激情 狠狠 欧美| 人妻少妇偷人精品九色| 美女中出高潮动态图| 肉色欧美久久久久久久蜜桃| 亚洲欧洲国产日韩| 久久久久久久精品精品| 日本黄色片子视频| 国产又色又爽无遮挡免| 欧美日韩一区二区视频在线观看视频在线| 99热国产这里只有精品6| 亚洲国产精品一区三区| 国产中年淑女户外野战色| 欧美一级a爱片免费观看看| 美女脱内裤让男人舔精品视频| 少妇人妻一区二区三区视频| 国产视频首页在线观看| 秋霞在线观看毛片| 80岁老熟妇乱子伦牲交| 两个人的视频大全免费| 97超视频在线观看视频| 久久午夜福利片| 久久精品国产亚洲av天美| 五月开心婷婷网| 久久99精品国语久久久| 亚洲av电影在线观看一区二区三区| 精品人妻熟女av久视频| 免费人成在线观看视频色| 亚洲电影在线观看av| 日日撸夜夜添| 亚洲高清免费不卡视频| 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 自拍偷自拍亚洲精品老妇| 麻豆国产97在线/欧美| 在线 av 中文字幕| 精品一区二区三卡| 大话2 男鬼变身卡| 在线观看三级黄色| 2021少妇久久久久久久久久久| 国产色爽女视频免费观看| 国内少妇人妻偷人精品xxx网站| 丰满人妻一区二区三区视频av| 97热精品久久久久久| 99久久人妻综合| 亚洲色图av天堂| 在线观看一区二区三区| 亚洲,一卡二卡三卡| 亚洲精品中文字幕在线视频 | 晚上一个人看的免费电影| 亚洲无线观看免费| 久久精品人妻少妇| a级毛片免费高清观看在线播放| 精品久久久精品久久久| 亚洲av在线观看美女高潮| h日本视频在线播放| 亚洲国产精品专区欧美| 免费观看无遮挡的男女| 狠狠精品人妻久久久久久综合| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 日本爱情动作片www.在线观看| 极品教师在线视频| 国产精品久久久久久精品电影小说 | 国产精品麻豆人妻色哟哟久久| av天堂中文字幕网| 一区二区av电影网| 久久国产亚洲av麻豆专区| 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 精华霜和精华液先用哪个| 欧美性感艳星| 久久久久精品性色| 在线观看免费视频网站a站| 亚洲精品aⅴ在线观看| 国国产精品蜜臀av免费| 亚洲欧美一区二区三区国产| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 免费久久久久久久精品成人欧美视频 | 久久久久久久久大av| 欧美成人一区二区免费高清观看| 熟妇人妻不卡中文字幕| 伊人久久精品亚洲午夜| 久久99蜜桃精品久久| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 日韩,欧美,国产一区二区三区| 免费人妻精品一区二区三区视频| 黑人猛操日本美女一级片| 久久久久性生活片| 国产伦在线观看视频一区| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 色综合色国产| 久久99蜜桃精品久久| av在线播放精品| 亚洲av综合色区一区| 91狼人影院| 日本av免费视频播放| 高清不卡的av网站| 亚洲av不卡在线观看| 三级经典国产精品| 18禁裸乳无遮挡动漫免费视频| 99久久精品国产国产毛片| 久热这里只有精品99| 国产极品天堂在线| 久久久久久久精品精品| 午夜精品国产一区二区电影| 亚洲中文av在线| 色婷婷av一区二区三区视频| 性色avwww在线观看| av专区在线播放| 午夜日本视频在线| 综合色丁香网| 狠狠精品人妻久久久久久综合| 大陆偷拍与自拍| 亚洲国产欧美在线一区| 亚洲欧美中文字幕日韩二区| 亚洲精品一二三| 黄片wwwwww| 中文字幕精品免费在线观看视频 | 国产精品一区二区三区四区免费观看| 搡女人真爽免费视频火全软件| 日韩免费高清中文字幕av| 欧美97在线视频| 这个男人来自地球电影免费观看 | 久久久国产一区二区| 视频中文字幕在线观看| 亚洲欧美中文字幕日韩二区| 能在线免费看毛片的网站| 青青草视频在线视频观看| 久久国产乱子免费精品| 日日啪夜夜撸| 日韩精品有码人妻一区| 日日撸夜夜添| 观看美女的网站| 黄色欧美视频在线观看| 黄片wwwwww| 欧美成人一区二区免费高清观看| 免费少妇av软件| 亚洲欧美精品专区久久| 性色avwww在线观看| 国产在线男女| 欧美日韩精品成人综合77777| 国语对白做爰xxxⅹ性视频网站| 国产伦精品一区二区三区四那| 欧美精品一区二区免费开放| 卡戴珊不雅视频在线播放| 国产亚洲av片在线观看秒播厂| 国内少妇人妻偷人精品xxx网站| 超碰av人人做人人爽久久| 欧美极品一区二区三区四区| 最近中文字幕2019免费版| 久久人人爽人人爽人人片va| 久久精品国产亚洲网站| 国产探花极品一区二区| 国产成人免费观看mmmm| 狂野欧美激情性bbbbbb| 国产探花极品一区二区| 97超视频在线观看视频| 国产乱来视频区| 99精国产麻豆久久婷婷| 色哟哟·www| 欧美日韩亚洲高清精品| 大话2 男鬼变身卡| 狂野欧美激情性xxxx在线观看| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 午夜免费鲁丝| tube8黄色片| 国产精品无大码| 国产在视频线精品| 人人妻人人爽人人添夜夜欢视频 | 国产乱人视频| av在线播放精品| 夜夜爽夜夜爽视频| 免费观看a级毛片全部| 在线 av 中文字幕| 18+在线观看网站| 国产一区二区在线观看日韩| 成人18禁高潮啪啪吃奶动态图 | 你懂的网址亚洲精品在线观看| 亚洲国产欧美人成| 特大巨黑吊av在线直播| 免费观看av网站的网址| 黑人猛操日本美女一级片| 亚洲精品国产av蜜桃| 欧美精品一区二区大全| 免费观看性生交大片5| av线在线观看网站| 久久精品人妻少妇| 欧美成人a在线观看| 久久久欧美国产精品| 亚洲美女视频黄频| 黄色视频在线播放观看不卡| 久久久久国产网址| 爱豆传媒免费全集在线观看| 日本黄色片子视频| 波野结衣二区三区在线| 亚洲国产欧美人成| 欧美成人精品欧美一级黄| 国产免费福利视频在线观看| 久久青草综合色| 成年女人在线观看亚洲视频| 日韩中字成人| 国产淫片久久久久久久久| 国产探花极品一区二区| 亚洲自偷自拍三级| 嘟嘟电影网在线观看| 在线看a的网站| 最近中文字幕高清免费大全6| 久久精品熟女亚洲av麻豆精品| 又大又黄又爽视频免费| 乱系列少妇在线播放| av天堂中文字幕网| 久久影院123| 亚洲第一av免费看| 欧美bdsm另类| 最近2019中文字幕mv第一页| 精品酒店卫生间| 97在线视频观看| 国产乱人偷精品视频| 久久久午夜欧美精品| 国产精品久久久久成人av| 久久久久久久久久成人| 国产国拍精品亚洲av在线观看| 国产91av在线免费观看| 午夜日本视频在线| 亚洲va在线va天堂va国产| 99国产精品免费福利视频| 一本—道久久a久久精品蜜桃钙片| 男人狂女人下面高潮的视频| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 赤兔流量卡办理| www.色视频.com| 亚洲成人一二三区av| 男女啪啪激烈高潮av片| 午夜福利影视在线免费观看| 免费大片黄手机在线观看| 精品一区在线观看国产| 涩涩av久久男人的天堂| 亚洲欧美精品专区久久| 国产精品久久久久久精品电影小说 | 这个男人来自地球电影免费观看 | 欧美变态另类bdsm刘玥| 久久久久久久久大av| 欧美 日韩 精品 国产| 在线观看av片永久免费下载| 亚洲精品国产av成人精品| 看免费成人av毛片| 国产视频首页在线观看| tube8黄色片| 婷婷色综合大香蕉| 少妇 在线观看| 寂寞人妻少妇视频99o| 美女福利国产在线 | 日本与韩国留学比较| 久久久久网色| 黑丝袜美女国产一区| 久久久精品94久久精品| 中文字幕人妻熟人妻熟丝袜美| 蜜臀久久99精品久久宅男| a 毛片基地| 全区人妻精品视频| 高清不卡的av网站| 亚洲无线观看免费| 2018国产大陆天天弄谢| 国产高清三级在线| 国产一区二区在线观看日韩| 久久人人爽人人片av| 成人毛片a级毛片在线播放| 99精国产麻豆久久婷婷| 七月丁香在线播放| 极品教师在线视频| 亚洲精品自拍成人| 麻豆成人午夜福利视频| 大陆偷拍与自拍| 午夜福利视频精品| 亚洲欧美日韩东京热| 欧美一区二区亚洲| 国产伦精品一区二区三区四那| 97超视频在线观看视频| 久久精品国产a三级三级三级| 欧美bdsm另类| 青春草国产在线视频| 国产高潮美女av| 成年人午夜在线观看视频| 欧美最新免费一区二区三区| 天美传媒精品一区二区| 久久久欧美国产精品| 国产精品欧美亚洲77777| 成人综合一区亚洲| 久久久久视频综合| 97在线人人人人妻| 日韩不卡一区二区三区视频在线| 超碰97精品在线观看| 蜜臀久久99精品久久宅男| 国产精品99久久99久久久不卡 | 亚洲av男天堂| 女的被弄到高潮叫床怎么办| 中文欧美无线码| 日韩一区二区三区影片| 在线看a的网站| 男男h啪啪无遮挡| 亚洲久久久国产精品| 在线观看美女被高潮喷水网站| 80岁老熟妇乱子伦牲交| 高清av免费在线| 涩涩av久久男人的天堂| 又大又黄又爽视频免费| 女的被弄到高潮叫床怎么办| 国产真实伦视频高清在线观看| 亚洲欧美日韩另类电影网站 | 精品久久国产蜜桃| 一级爰片在线观看| 99久久中文字幕三级久久日本| 久热久热在线精品观看| 国产精品偷伦视频观看了| 国产精品久久久久久久电影| 久久 成人 亚洲| 91狼人影院| 欧美国产精品一级二级三级 | 国产淫语在线视频| 亚洲av.av天堂| 亚洲国产精品国产精品| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 热re99久久精品国产66热6| 久久99精品国语久久久| 久久久久性生活片| 亚洲一区二区三区欧美精品| 人妻少妇偷人精品九色| 亚洲欧美一区二区三区黑人 | 久久久久精品久久久久真实原创| 蜜桃亚洲精品一区二区三区| 少妇人妻一区二区三区视频| 秋霞在线观看毛片| 亚洲av中文av极速乱| 国产黄色视频一区二区在线观看| 91久久精品电影网| 亚洲国产精品成人久久小说| 国产精品一区二区性色av| 亚洲精品日本国产第一区| 高清黄色对白视频在线免费看 | 国产精品久久久久久精品电影小说 | 久久久久久久久大av| 精品99又大又爽又粗少妇毛片| av国产免费在线观看| 干丝袜人妻中文字幕| 91aial.com中文字幕在线观看| 热re99久久精品国产66热6| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的 | 成人黄色视频免费在线看| 如何舔出高潮| 久久精品久久久久久噜噜老黄| av黄色大香蕉| 美女国产视频在线观看| av在线app专区| 日本vs欧美在线观看视频 | 最近最新中文字幕大全电影3| 成人国产麻豆网| 97在线人人人人妻| 97超视频在线观看视频| 日韩电影二区| 国产精品一区二区在线不卡| 99re6热这里在线精品视频| 91在线精品国自产拍蜜月| 五月开心婷婷网| 好男人视频免费观看在线| 波野结衣二区三区在线| 校园人妻丝袜中文字幕| 97超碰精品成人国产| 天堂8中文在线网| 欧美激情极品国产一区二区三区 | 丝袜脚勾引网站| 亚洲第一区二区三区不卡| 亚洲欧美清纯卡通| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 在线看a的网站| 国产成人aa在线观看| 国产极品天堂在线| 成人二区视频| 男女下面进入的视频免费午夜| 高清毛片免费看| 最近最新中文字幕免费大全7| 国产日韩欧美在线精品| 有码 亚洲区| 色5月婷婷丁香| 久久6这里有精品| 久久久久久久久久人人人人人人| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 天天躁夜夜躁狠狠久久av| 亚洲精品成人av观看孕妇| 精品人妻视频免费看| 欧美极品一区二区三区四区| 99热网站在线观看| 日韩视频在线欧美| 欧美国产精品一级二级三级 | 欧美人与善性xxx| 午夜福利视频精品| 天堂中文最新版在线下载| 国产 精品1| 亚洲欧美精品专区久久| 中文字幕亚洲精品专区| 久久精品久久久久久噜噜老黄| 极品少妇高潮喷水抽搐| 少妇的逼好多水| 91狼人影院| 好男人视频免费观看在线| 国产免费福利视频在线观看| 欧美xxxx性猛交bbbb| 中文欧美无线码| av线在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 人妻 亚洲 视频| 观看美女的网站| 成人免费观看视频高清| 婷婷色综合www| 亚洲精品自拍成人| 欧美 日韩 精品 国产| 丰满人妻一区二区三区视频av| 赤兔流量卡办理| 免费观看性生交大片5| 国产黄片美女视频| 成年人午夜在线观看视频| 日本与韩国留学比较| 国产在线一区二区三区精| 九色成人免费人妻av|