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

    流場(chǎng)中填充吸聲材料夾層板結(jié)構(gòu)的聲振耦合特性

    2016-04-07 07:49:32寧少武史治宇李曉松
    振動(dòng)與沖擊 2016年3期

    寧少武, 史治宇, 李曉松

    (1.南京航空航天大學(xué) 機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京 210016; 2.中國(guó)運(yùn)載火箭研究院,北京 100076)

    ?

    流場(chǎng)中填充吸聲材料夾層板結(jié)構(gòu)的聲振耦合特性

    寧少武1, 史治宇1, 李曉松2

    (1.南京航空航天大學(xué) 機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京210016; 2.中國(guó)運(yùn)載火箭研究院,北京100076)

    摘要:采用等效流體模擬吸聲材料,建立了外部流場(chǎng)作用下填充吸聲材料夾層板結(jié)構(gòu)的聲振耦合模型,應(yīng)用波動(dòng)分析方法研究結(jié)構(gòu)中聲的透射特性,分析了入射聲波入射角和方位角、流場(chǎng)流速和流向、夾層結(jié)構(gòu)幾何尺寸等參數(shù)對(duì)填充吸聲材料夾層板結(jié)構(gòu)聲振耦合特性的影響。仿真計(jì)算表明吸聲材料提高了雙層板結(jié)構(gòu)的隔聲性能;隔聲性能隨著面板厚度和夾層厚度的增加而提高,隨著入射角和方位角的增大而減??;在計(jì)算頻段內(nèi)(0~5 000 Hz),逆流入射時(shí)傳聲損失隨著馬赫數(shù)的增大而減小,順流入射時(shí)卻隨著馬赫數(shù)的增大而增大。

    關(guān)鍵詞:等效流體模型;夾層板結(jié)構(gòu);吸聲材料;聲振耦合模型;駐波共振頻率

    艙內(nèi)噪聲嚴(yán)重影響飛機(jī)機(jī)艙的舒適性。艙內(nèi)噪聲主要源于艙外的湍流邊界層激勵(lì)噪聲和發(fā)動(dòng)機(jī)噴氣噪聲透過(guò)機(jī)身結(jié)構(gòu)進(jìn)入機(jī)艙內(nèi)部[1-2]。機(jī)艙外部有高速流體,飛機(jī)外蒙皮和艙內(nèi)裝飾板構(gòu)成類(lèi)似雙層板結(jié)構(gòu),在內(nèi)外結(jié)構(gòu)之間填充隔聲隔熱材料,研究這種更加接近實(shí)際工程應(yīng)用的復(fù)雜結(jié)構(gòu)系統(tǒng)的聲振耦合特性對(duì)飛機(jī)結(jié)構(gòu)低噪聲設(shè)計(jì)和艙內(nèi)噪聲的主被動(dòng)控制至關(guān)重要[3-6]。

    對(duì)雙板空腔結(jié)構(gòu)聲振耦合特性的研究通常簡(jiǎn)化為無(wú)限大結(jié)構(gòu),沒(méi)有考慮結(jié)構(gòu)邊界條件的影響[7-8],在文獻(xiàn)[9]中采用了平面波空間窗截?cái)嗟姆椒▉?lái)近似模擬有限大結(jié)構(gòu)的傳聲特性;文獻(xiàn)[10]研究了有限大雙板空腔結(jié)構(gòu)聲振耦合特性,充分考慮了邊界(簡(jiǎn)支和固支)條件對(duì)結(jié)構(gòu)聲振耦合特性的影響,并進(jìn)行了相應(yīng)的實(shí)驗(yàn)測(cè)量研究。文獻(xiàn)[11-13]系統(tǒng)研究了在外部流動(dòng)流體作用下結(jié)構(gòu)的聲振耦合問(wèn)題,建立了基本的理論模型,分析研究了湍流邊界層的影響;文獻(xiàn)[14]中則研究了外部平均流體對(duì)聲波穿過(guò)雙層彈性板結(jié)構(gòu)的影響,研究揭示外部流體對(duì)聲波透射的影響機(jī)理。

    對(duì)多孔吸聲材料主要有兩類(lèi)處理方法[15-16]:其一就是采用嚴(yán)格的Biot理論模型[17],該理論將連續(xù)介質(zhì)力學(xué)應(yīng)用于流體飽和的多孔介質(zhì)體系,考慮了流體和固體骨架的應(yīng)力應(yīng)變關(guān)系及運(yùn)動(dòng)特性,其求解通常借助于有限元數(shù)值方法,計(jì)算量大;其二就是采用等效流體模型[18],即將多孔材料等效為具有同樣動(dòng)態(tài)密度和動(dòng)態(tài)體積模量的流體,其假設(shè)多孔材料是剛性骨架,在飽和的空氣介質(zhì)中存在壓縮波,且滿足Helmholtz方程。文獻(xiàn)[19]采用等效流體模型分別研究了填充吸聲材料夾層板結(jié)構(gòu)的聲輻射和聲透射特性,分析了填充吸聲材料對(duì)正交加筋板結(jié)構(gòu)聲輻射和聲透射的影響,在本文中則研究了外部流場(chǎng)作用下填充吸聲材料夾層板結(jié)構(gòu)的聲振耦合特性。

    夾層板結(jié)構(gòu)在航空結(jié)構(gòu)中廣泛采用,其由外蒙皮-中間隔熱隔聲材料-內(nèi)裝飾板組成,艙外有高速流體,而上述已有的計(jì)算模型只考慮部分因素對(duì)結(jié)構(gòu)聲振耦合特性的影響;本文以飛機(jī)機(jī)身結(jié)構(gòu)的聲振特性研究為背景,綜合考慮各方面因素,采用等效流體模擬吸聲材料,建立外部流場(chǎng)作用下外蒙皮-中間隔熱隔聲材料-內(nèi)裝飾板構(gòu)成的復(fù)雜結(jié)構(gòu)-聲學(xué)系統(tǒng)的聲振耦合模型,應(yīng)用波動(dòng)分析方法研究結(jié)構(gòu)中聲的透射特性。

    1等效流體模型

    本文應(yīng)用等效流體模型模擬聲波在填充吸聲材料中的傳播。按照等效流體模型[16,18-19],聲波在多孔纖維吸聲材料中的傳播滿足波動(dòng)方程

    (1)

    式中:pm為多孔纖維吸聲材料中的聲壓;km為相應(yīng)的復(fù)波數(shù),

    K(ω)=

    G1(ρ0f/R)=[1+iπ(ρ0f/R)]1/2

    G2(ρ0f/R)=G1[(ρ0f/R)4Npr]

    (2)

    式中:ρ(ω)為多孔纖維吸聲材料的動(dòng)態(tài)密度;K(ω)為多孔纖維吸聲材料的動(dòng)態(tài)體積模量;R為多孔纖維材料的靜態(tài)流阻抗;γ和ρ0分別為空氣的比熱容率和密度;P0為大氣壓;Npr為普朗特常數(shù);k0為空氣中聲波波數(shù),k0=ω/c0,c0為空氣中聲波傳播速度;ρm為填充多孔材料的復(fù)密度;σ為填充多孔材料的孔隙率。從上式可以看出km與多孔纖維吸聲材料的動(dòng)態(tài)密度ρ(ω)和動(dòng)態(tài)體積模量K(ω)相關(guān),這兩個(gè)變量都強(qiáng)烈地依賴(lài)于ρ0f/R,反映了吸聲過(guò)程本身隨頻率變化的動(dòng)態(tài)特性和吸聲機(jī)理主要源于材料的流阻抗作用[16]。

    2外部流場(chǎng)作用下填充吸聲材料夾層板結(jié)構(gòu)的傳聲特性

    如圖1所示,平行的雙層彈性薄板簡(jiǎn)支安裝在剛性聲障上,長(zhǎng)寬分別為a和b,將聲場(chǎng)分為入射聲場(chǎng)1、中間聲場(chǎng)2和輻射聲場(chǎng)3。在中間聲場(chǎng)2中填充多孔纖維吸聲材料,厚度為H。在入聲聲場(chǎng)1中存在平行于板平面沿x軸方向的理想平均流,流速為υ。入射聲波的方位角為α,入射角為θ1。平面簡(jiǎn)諧入射聲波的聲壓速度勢(shì)為

    φ=Ie-j(k1xx+k1yy+k1zz-ωt)

    k1x=k1sinθ1cosα

    k1y=k1sinθ1sinα

    k1z=k1cosθ1

    圖1 耦合系統(tǒng)示意圖Fig.1 The sketch of the coupling system

    2.1板的振動(dòng)控制方程

    在簡(jiǎn)諧平面入射聲波的激勵(lì)下,雙層彈性板的振動(dòng)控制方程分別為

    (3)

    (4)

    式中:D1和D2分別為入射板和輻射板的彎曲剛度;w1和w2分別為入射板和輻射板的橫向振動(dòng)位移;m1和m2分別為入射板和輻射板的面密度;Φ1、Φ2和Φ3分別為入射聲場(chǎng)1、中間聲場(chǎng)2和輻射聲場(chǎng)3的速度勢(shì)函數(shù),其中中間聲場(chǎng)2采用上述等效流體模型,模擬聲在多孔吸聲材料中的傳播。

    對(duì)于簡(jiǎn)支邊界條件,則板在邊界處的橫向位移和彎矩都應(yīng)該等于零,即

    2.2聲場(chǎng)的聲學(xué)波動(dòng)方程和邊界條件

    理想流體平行于板平面流動(dòng),入射聲場(chǎng)1滿足對(duì)流波動(dòng)方程

    (5)

    式中:Ux為理想流體流速。

    在入射聲場(chǎng)1中流固界面應(yīng)滿足位移連續(xù)性條件,即

    (6)

    式中:p1為入射聲場(chǎng)聲壓,

    中間聲場(chǎng)2滿足波動(dòng)方程式(1),聲壓pm采用速度勢(shì)函數(shù)表示為

    pm=jωρmΦ2

    輻射聲場(chǎng)3滿足聲場(chǎng)波動(dòng)方程

    p3=jωρ0Φ3

    式中:p3為輻射聲場(chǎng)聲壓。在中間聲場(chǎng)2和輻射聲場(chǎng)3的流固界面上滿足法向速度連續(xù)條件,即速度連續(xù)條件

    (7)

    (8)

    3外部流場(chǎng)作用下填充吸聲材料夾層板結(jié)構(gòu)傳聲特性求解

    雙層彈性薄板為簡(jiǎn)支支撐,則其振動(dòng)位移可以表示成簡(jiǎn)支模態(tài)函數(shù)的形式

    (9)

    (10)

    式中:α1,mn和α2,mn分別為上下兩板的振動(dòng)模態(tài)系數(shù),φmn(x,y)為簡(jiǎn)支模態(tài)函數(shù),

    簡(jiǎn)支雙板的剛性聲障約束密封腔內(nèi)的空氣,使得空腔內(nèi)平行于板平面的聲壓分布與板簡(jiǎn)支模態(tài)函數(shù)的級(jí)數(shù)形式相同[16],因此聲場(chǎng)聲壓速度勢(shì)函數(shù)可以表示為

    (11)

    (12)

    (13)

    式中:kmz和k3z分別為中間聲場(chǎng)2和輻射聲場(chǎng)3中z向波數(shù)分量;Imn、εmn和ξmn分別為入射聲場(chǎng)1、中間聲場(chǎng)2和輻射聲場(chǎng)3中正行波的幅值;βmn和ζmn分別為入射聲場(chǎng)1和中間聲場(chǎng)2中負(fù)行波的幅值。入射波幅值Imn由入射聲波的聲壓速度勢(shì)φ得到

    中間聲場(chǎng)聲壓和煤質(zhì)質(zhì)點(diǎn)速度分別為

    pm=jωρmΦ2

    式中:vmx和vmy分別為x方向和y方向的質(zhì)點(diǎn)速度。從上式可以看出:在中間聲場(chǎng)四周邊界處聲壓為零,在聲腔中形成邊界聲壓為零的駐波聲場(chǎng),因此Φ2形式的速度勢(shì)函數(shù)可以理解為在中間聲場(chǎng)的四周邊界上采用了絕對(duì)軟邊界的假設(shè);而剛性邊界假設(shè)形成邊界質(zhì)點(diǎn)速度為零的駐波聲場(chǎng)。

    首先,將式(11)代入方程(5)可以得到

    (14)

    式中:Ma=Ux/c0為入射聲場(chǎng)中平均流在x向的馬赫數(shù)。

    其次,在與板相鄰的流體中傳播的聲波波長(zhǎng)應(yīng)該與在板中傳播的結(jié)構(gòu)彎曲波波長(zhǎng)相一致[20],即

    k1x=kmx=k3x=ktcosα=k1sinθ1cosα

    k1y=kmy=k3y=ktsinα=k1sinθ1sinα

    式中:kmx和k3x分別為中間聲場(chǎng)2和輻射聲場(chǎng)3中x向波數(shù)分量;kmy和k3y分別為中間聲場(chǎng)2和輻射聲場(chǎng)3中y向波數(shù)分量,則有

    式中km由上述等效流體模型確定,為中間聲場(chǎng)2中的聲波復(fù)波數(shù);k0為空氣中聲波波數(shù)。將式(11)~式(13)代入邊界條件式(6)~式(8),聯(lián)立求解,可得

    將式(9)~式(10),式(11)~式(13)代入方程(3)和(4)中,利用簡(jiǎn)支模態(tài)函數(shù)的正交性,整理可得

    式中:

    通過(guò)上述計(jì)算求解,可以得到上下兩板的振動(dòng)模態(tài)系數(shù)α1,mn和α2,mn,從而求得εmn、ξmn、βmn和ζmn的值。

    聲場(chǎng)的聲功率定義為

    (15)

    式中:ΔSi為第i個(gè)單元面積,Nxy為單元的總數(shù),Nxy=Nx×Ny,pi為第i個(gè)結(jié)構(gòu)單元表面附近聲場(chǎng)聲壓;vi為第i個(gè)結(jié)構(gòu)單元表面附近聲場(chǎng)的質(zhì)點(diǎn)速度。

    4仿真計(jì)算

    為了驗(yàn)證本文理論研究的正確性,同時(shí)計(jì)算了空氣夾層和填充吸聲材料夾層雙層板結(jié)構(gòu)的聲振耦合特性,對(duì)比分析入射角、方位角、馬赫數(shù)、面板尺寸、面板厚度以及夾層厚度等對(duì)結(jié)構(gòu)聲振耦合特性的影響。

    圖1中,上下兩層鋁合金板的結(jié)構(gòu)尺寸為:長(zhǎng)度a=1.0 m,寬度b=1.0 m,厚度t=0.005 m,中間聲腔厚度H=0.06 m??紤]材料阻尼的影響,采用復(fù)彈性模量E(1+iχ),其中,彈性模量E=70 GPa,阻尼損耗因子χ=0.01,泊松比υ=0.33,則彎曲剛度

    密度ρ=2 700 kg/m3。多孔纖維材料的靜態(tài)流阻抗R=20 000 N·m/s4,空氣的比熱容率γ=cp/cv=1.4,密度ρ0=1.21 kg/m3,空氣中的聲速c0=343 m/s,大氣壓P0=101 920 Pa,普朗特常數(shù)Npr=0.702,多孔纖維材料的孔隙率σ=0.62。入射聲場(chǎng)中,流場(chǎng)流速的馬赫數(shù)Ma=0.6,入射角θ1和方位角α都為30°。

    本文計(jì)算隔聲結(jié)構(gòu)在頻率區(qū)間0~5 000 Hz內(nèi)的隔聲性能。為了保證式(15)的計(jì)算精度,單元長(zhǎng)度Δ0=c0/(4fmax),Nx=Ny=[max{a/Δ0,b/Δ0}],fmax為最大計(jì)算頻率,[·]表示向上取整;在5 000 Hz時(shí),式(9)和式(10)中模態(tài)階數(shù)M和N對(duì)傳聲損失收斂性的影響如圖2所示,對(duì)不同尺寸的結(jié)構(gòu)按照?qǐng)D中所示選取M和N階數(shù),a=b取0.5 m,1.0 m,和1.5 m時(shí),M=N分別取20,40和60。

    圖2 耦合系統(tǒng)的傳聲損失級(jí)數(shù)解收斂性檢驗(yàn)Fig.2Convergencecheckofsoundtransmissionlossseriessolutionforthecouplingsystem圖3 入射角和方位角對(duì)有限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.3EffectsofincidentanglesandazimuchanglesonSTLoffinitesandwichstructurewithabsorptivematerial圖4 不同入射角和方位角對(duì)無(wú)限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.4EffectsofincidentanglesandazimuchanglesonSTLofinfinitesandwichstructurewithabsorptivematerial

    圖3和圖4分別對(duì)應(yīng)有限尺寸和無(wú)限尺寸系統(tǒng),對(duì)比了入射角和方位角對(duì)填充吸聲材料夾層板結(jié)構(gòu)傳聲損失的影響,夾層包括空氣(Air)和吸聲材料(AM)兩種情況。有限尺寸雙層板結(jié)構(gòu)上下面板尺寸為1.0×1.0m2,無(wú)限大雙層板結(jié)構(gòu)上下面板尺寸取108×108m2來(lái)代替無(wú)限尺寸平板。從圖3和圖4中可以看出,填充吸聲材料提高了結(jié)構(gòu)的隔聲性能。

    當(dāng)為空氣夾層時(shí),“板-空氣-板”的共振頻率和駐波共振頻率[14]分別為

    (16)

    (17)

    將中間吸聲材料等效為流體,同理可得“板-等效流體-板”共振頻率為

    (18)

    式中:cm為等效流體中聲速,

    (19)

    對(duì)于有吸聲材料的結(jié)構(gòu),多孔纖維材料的等效密度和等效體積模量是隨頻率變化的動(dòng)態(tài)函數(shù),因此中間夾層的折射角θ2和中間夾層的聲速cm都隨頻率變化;特別地,當(dāng)θ1=0°時(shí),cosθ2=1。代入計(jì)算數(shù)據(jù),可得:θ1=0°,α=0°,fα=94.4 Hz,fm,α=101.3 Hz,fd,1=2 858.0 Hz;θ1=30°,α=30°,fα=102.8 Hz,fd,1=3 114.1 Hz;θ1=30°,α=45°,fα=103.58 Hz,fd,1=3 137.7Hz;θ1=45°,α=30°,fα=110.2 Hz,fd,1=3 339.5 Hz。從圖3和圖4中可以看出,空氣中“板-空氣-板”共振頻率、駐波共振頻率以及“板-等效流體-板”共振頻率的理論和數(shù)值計(jì)算結(jié)果相符很好。對(duì)于空氣夾層,隨著入射角的增大,“板-空氣-板”共振頻率和駐波共振頻率向高頻移動(dòng);隨著方位角的增大,“板-空氣-板”共振頻率和駐波共振頻率也向高頻移動(dòng),但入射角的變化對(duì)其影響更明顯。對(duì)于填充吸聲材料夾層,“板-等效流體-板”共振頻率對(duì)入射角和方位角的變化并不敏感,傳聲損失隨著入射角和方位角的增大而減小。

    對(duì)比圖3和圖4也可以看出:當(dāng)夾層為空氣和吸聲材料時(shí),無(wú)限大尺寸結(jié)構(gòu)的傳聲損失曲線遠(yuǎn)比有限尺寸結(jié)構(gòu)的傳聲曲線光滑,這是由于有限大尺寸結(jié)構(gòu)比無(wú)限大尺寸結(jié)構(gòu)有更密集的模態(tài),也反映出采用無(wú)限大結(jié)構(gòu)代替有限尺寸結(jié)構(gòu)計(jì)算傳聲損失是不可取的[16]。

    圖5~圖8對(duì)比了在不同馬赫數(shù)流體中順流和逆流入射時(shí)對(duì)結(jié)構(gòu)傳聲損失的影響,其中θ1=α=30°。

    圖5 在逆流入射情況下馬赫數(shù)對(duì)有限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.5EffectsofMachnumberonSTLoffinitesandwichstructurewithabsorptivematerialintheupstream圖6 在逆流入射情況下馬赫數(shù)對(duì)無(wú)限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.6EffectsofMachnumberonSTLofinfinitesandwichstructurewithabsorptivematerialintheupstream圖7 在順流入射情況下馬赫數(shù)對(duì)有限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.7EffectsofMachnumberonSTLoffinitesandwichstructurewithabsorptivematerialinthedownstream

    由式(16~19)可得知當(dāng)cosθ2→0時(shí),fα,fm,α和fd,n趨于無(wú)窮而消失。

    當(dāng)

    則有

    當(dāng)為空氣夾層時(shí),km=k0,代入相關(guān)計(jì)算數(shù)據(jù),則有Ma>-1.154 7,即對(duì)于逆流入射,當(dāng)馬赫數(shù)大于1.154 7時(shí),fα,fm,α和fd,n趨于無(wú)窮,即“板-空氣-板”共振和駐波共振消失。從圖6可以看出,當(dāng)Ma=-1.2時(shí),空氣夾層雙層板結(jié)構(gòu)的“板-空氣-板”共振和駐波共振消失。對(duì)于填充吸聲材料的夾層板結(jié)構(gòu),km/k0與頻率有關(guān),傳聲損失曲線與空氣夾層的結(jié)構(gòu)完全不同;從仿真計(jì)算過(guò)程中得到(km/k0)>1.1,取km/k0=1.1,計(jì)算可以得到Ma>-1.26,當(dāng)Ma=-1.2時(shí),“板-等效流體-板”共振頻率依然存在,從圖6完全證明這一點(diǎn)。代入計(jì)算數(shù)據(jù),空氣夾層板結(jié)構(gòu)的“板-空氣-板”共振頻率和駐波共振頻率分別為:Ma=-0.6,fα=128.0 Hz,fd,1=3 876.4 Hz;Ma=0,fα=108.0 Hz,fd,1=3 300.5 Hz;Ma=0.6,fα=102.8 Hz,fd,1=3 114.1 Hz;Ma=1.2,fα=99.9 Hz,fd,1=3 027.0 Hz。

    圖8 在順流入射情況下馬赫數(shù)對(duì)無(wú)限尺寸填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.8 Effects of Mach number on STL of infinite sandwich structure with absorptive material in the downstream

    從圖6~圖8和上面的計(jì)算結(jié)果可以看出,對(duì)于空氣夾層的雙層板結(jié)構(gòu),當(dāng)逆流入射時(shí),隨著Ma的增大,“板-等效流體-板”共振頻率和駐波共振頻率向高頻

    移動(dòng),但當(dāng)Ma大于某一馬赫數(shù)以后,因發(fā)生全反射現(xiàn)象而消失;在計(jì)算頻段內(nèi)(0~5 000 Hz),傳聲損失隨著Ma的增大而減小。當(dāng)順流入射時(shí),隨著Ma的增大,“板-等效流體-板”共振頻率和駐波共振頻率向低頻移動(dòng);在計(jì)算頻段內(nèi)(0~5 000 Hz),傳聲損失隨著Ma的增大而減大。

    對(duì)于填充吸聲材料的雙層板結(jié)構(gòu),等效密度和等效體積模量是隨頻率變化的動(dòng)態(tài)函數(shù),反映了吸聲過(guò)程隨頻率變化的動(dòng)態(tài)特性,沒(méi)有駐波共振頻率的出現(xiàn),且延遲了全反射現(xiàn)象的出現(xiàn);當(dāng)逆流入射時(shí),在計(jì)算頻段內(nèi)(0~5 000 Hz)傳聲損失隨著Ma的增大而減??;當(dāng)順流入射時(shí),在計(jì)算頻段內(nèi)(0~5 000 Hz)傳聲損失隨著Ma的增大而增大。

    圖9~圖11分別對(duì)比了不同尺寸大小、不同面板厚度和不同夾層厚度對(duì)雙層板結(jié)構(gòu)傳聲損失的影響。

    圖9中雙層板結(jié)構(gòu)具有相同的夾層厚度和面板面密度,因此具有相同的“板-等效流體-板”共振頻率。由于支撐邊界的影響,在頻率小于“板-等效流體-板”共振頻率的低頻段,板的尺寸越小,隔聲性能越好,即低頻隔聲性能越優(yōu),其中無(wú)限大尺寸板對(duì)應(yīng)的隔聲性能曲線為其隔聲性能下限;在頻率大于“板-等效流體-板”共振頻率的高頻段,板的尺寸越大,隔聲性能越好,即高頻隔聲性能越優(yōu),其中無(wú)限大尺寸板對(duì)應(yīng)的隔聲性能曲線為其隔聲性能上限。

    圖9 結(jié)構(gòu)尺寸對(duì)填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.9EffectsofpanelsizeonSTLofsandwichstructurewithabsorptivematerial圖10 面板厚度對(duì)填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.10EffectsofpanelthicknessonSTLofsandwichstructurewithabsorptivematerial圖11 吸聲材料厚度對(duì)填充吸聲材料夾層板結(jié)構(gòu)的傳聲損失的影響Fig.11EffectsofabsorptionmaterialthicknessonSTLofsandwichstructurewithabsorptivematerial

    由式(18)可知“板-等效流體-板”共振頻率與雙層板夾層厚度以及上下面板面密度密切相關(guān),雙層板結(jié)構(gòu)的夾層厚度和上下面板厚度越大,相應(yīng)的“板-等效流體-板”共振頻率越低。從圖10和圖11可以看出:隨著面板厚度和中間夾層厚度的增加,結(jié)構(gòu)的隔聲性能都提高,而“板-等效流體-板”共振頻率向低頻移動(dòng)。

    5結(jié)論

    本文采用等效流體模型模擬了聲波在吸聲材料中的傳播,討論了外部流場(chǎng)作用下填充吸聲材料夾層板結(jié)構(gòu)的聲振耦合特性,為噪聲控制策略選取提供理論依據(jù),通過(guò)計(jì)算分析得出以下幾點(diǎn)結(jié)論:

    (1)填充吸聲材料提高了雙層板結(jié)構(gòu)的隔聲性能,消除了采用空氣夾層時(shí)出現(xiàn)的由駐波共振引起的隔聲低谷。

    (2)傳聲損失隨著入射角和方位角的增大而減小。對(duì)于空氣夾層,“板-空氣-板”共振頻率和駐波共振頻率向高頻移動(dòng),且對(duì)入射角的變化更敏感;對(duì)于填充吸聲材料夾層,結(jié)構(gòu)的聲學(xué)性能對(duì)入射角和方位角的變化不敏感;

    (3)對(duì)于空氣和填充吸聲材料的夾層板結(jié)構(gòu),當(dāng)逆流入射時(shí),在計(jì)算頻段內(nèi)(0~5 000 Hz)傳聲損失隨著Ma的增大而減??;當(dāng)順流入射時(shí),在計(jì)算頻段內(nèi)(0~5 000 Hz)傳聲損失隨著Ma的增大而增大。對(duì)于逆流入射情況,當(dāng)馬赫數(shù)大于一定值以后會(huì)出現(xiàn)全反射現(xiàn)象,“板-空氣-板”共振頻率、駐波共振頻率和“板-等效流體-板”共振頻率消失,但在填充吸聲材料結(jié)構(gòu)中出現(xiàn)全反射現(xiàn)象需要更大馬赫數(shù);

    (4)由于支撐邊界的影響,夾層板結(jié)構(gòu)上下面板尺寸越小,對(duì)低頻噪聲隔聲性能越好,而對(duì)中高頻噪聲隔聲性能越差;反之,對(duì)低頻噪聲隔聲性能越差,而對(duì)中高頻噪聲隔聲性能越好;

    (5)夾層板結(jié)構(gòu)上下面板和中間夾層厚度越厚,噪聲隔聲性能越好,且“板-等效流體-板”共振頻率向低頻移動(dòng)。

    參 考 文 獻(xiàn)

    [1] Howe M S, Shah P L.Influence of mean flow on boundary layer generated interior noise[J].Journal of the Acoustical Society of America, 1996, 99(6): 3401-3411.

    [2] Graham W R.Boundary layer induced noise in aircraft, part I & II [J]. Journal of Sound and vibration, 1996, 192(1): 101-138.

    [3] Gardonio P, Elliott S J. Active control of structure-borne and airborne sound transmission through doublepanel[J].Journal of Aircraft, 1999, 36(6): 1023-1032.

    [4] Maury C, Gardonio P, Elliott S J. Model for active control of flow-induced noise transmitted through doublepartitions[J].AIAA Journal, 2002, 40(6): 1113-1121.

    [5] Pietrzko S J, Mao Q. New results in active and passive control of sound transmission through double wallstructures[J].Aerospace Science and Technology, 2008, 12(1): 42-53.

    [6] Carneal J P, Fuller C R.An analytical and experimental investigation of active structural acoustic control of noise transmission through double panel systems[J].Journal of Sound and Vibration, 2004, 272: 749-771.

    [7] Antonio J M P, Tadeu A, Godinho L. Analytical evaluation of the acoustic insulation provided by double infinite walls[J].Journal of Sound and Vibration, 2003, 263(1): 113-129.

    [8] Brunskog J. The influence of finite cavities on the sound insulation of double-plate structures[J].Acta Astronautica, 2005, 117(6): 3727-3739.

    [9] Villot M, Guigou C, Gagliardini L. Predicting the acoustical radiation of finite size multi-layered structures by applying spatial windowing on infinite structures[J]. Journal of Sound and Vibration, 2001, 245(3): 433-455.

    [10] Xin F X, Lu T J. Analytical and experimental investigation on transmission loss of clamped double panels: Implication of boundary effects[J].Journal of the Acoustical Society of America, 2009, 125(3): 1506-1517.

    [11] Frampton K D. Radiation efficiency of convected fluid-loaded plates[J].Journal of the Acoustical Society of America, 2003, 113(5): 2663-2673.

    [12] Frampton K D. The effect of flow-induced coupling on sound radiation from convected fluid loaded plates[J].Journal of the Acoustical Society of America, 2005, 117(3): 1129-1137.

    [13] Schmidt P L, Frampton K D.Effect of in-plane forces on sound radiation from convected, fluid loaded plates[J]. Journal of Vibration and Acoustics, 2009, 131(021001): 1-8.

    [14] Xin F X, Lu T J, Chen C Q. External mean flow influence on noise transmission through double-leaf aeroelastic plates[J].AIAA Journal, 2009, 47(8): 1939-1951.

    [15] Panneton R, Atalla N. Numerical prediction of sound transmission through finite multilayer systems with poroelastic materials[J].Journal of the Acoustical Society of America, 1996, 100(1): 346-354.

    [16] 盧天健,辛鋒先. 輕質(zhì)板殼結(jié)構(gòu)設(shè)計(jì)的振動(dòng)和聲學(xué)基礎(chǔ)[M].北京:科學(xué)出版社,2012: 21.

    [17] Biot M A.Theory of propagation of elastic waves in a fluid-saturated porous solid.Ⅰ & Ⅱ[J].Journal of the Acoustical society of America, 1956, 28(2): 168-191.

    [18] Allard J F, Champoux Y.New empirical equations for sound propagation in rigid frame fibrous materials[J].Journal of the Acoustical Society of America,1992, 91(6): 3346-3353.

    [19] Xin F X, Lu T J. Transmission loss of orthogonally rib-stiffened double-panel structures with cavity absorption[J].Journal of the Acoustical Society of America, 2011, 12 (4): 1919-1934.

    [20] Koval L R. Effect of air flow, panel curvature, and internal pressurization on field-incidence transmission loss[J].Journal of the Acoustical Society of America, 1976, 59(6): 1379-1385.

    Vibro-acoustic coupling characteristics of the sandwich structure with absorptive material in external convected fluids

    NINGShao-wu1,SHIZhi-yu1,LIXiao-song2

    (1. State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China;2. China Academy of Launch Vehicle Technology, Beijing 100076, China)

    Abstract:An equivalent fluid model was employed to characterize the absorption of sound in the sound absorptive material. A vibro-acoustic coupling model was developed for the sound insulation of an sandwich structure filled with sound absorptive material in convected fluids. The performance of sound transmission was analysed by employing the wave method. The analysed influential factors of vibro-acoustic responses include the incident angles and azimuch angles, the velocity and direction of convected flow and the geometrical dimensions of the double panels. The studies show that the insulation of the structure filled with absorptive material instead of air is improved; the larger the thicknesses of the up and low panels and the gap are, the larger the sound transmission loss is; the larger the incident elevation angles and azimuch angles are, the smaller the sound transmission loss is; in the calculation frequency band(0~5 000 Hz), the sound transmission loss decreases with the increase of Mach number when the sound is incident in the upstream but increases with the increase of Mach number when the sound is incident in the downstream.

    Key words:equivalent fluid model; sandwich structure; absorptive material; vibro-acoustic coupling model; standing-wave resonance frequency

    中圖分類(lèi)號(hào):TB535

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.13465/j.cnki.jvs.2016.03.017

    通信作者史治宇 男,教授,博士生導(dǎo)師,1967年生

    收稿日期:2014-11-27修改稿收到日期:2015-01-30

    基金項(xiàng)目:江蘇省普通高校研究生科研創(chuàng)新計(jì)劃資助項(xiàng)目(CXZZ13_0147); 中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金資助; 江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目; 國(guó)家自然基金(11172131; 11232007)

    第一作者 寧少武 男,博士生,1985年生

    郵箱: zyshi@nuaa.edu.cn

    久久久色成人| 永久免费av网站大全| 精品一区在线观看国产| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久电影| 成人综合一区亚洲| 狠狠精品人妻久久久久久综合| 麻豆精品久久久久久蜜桃| 亚洲高清免费不卡视频| 97热精品久久久久久| 啦啦啦中文免费视频观看日本| 有码 亚洲区| 一边亲一边摸免费视频| 免费观看a级毛片全部| 寂寞人妻少妇视频99o| 国产美女午夜福利| 天天躁日日操中文字幕| 亚洲精品国产成人久久av| 久久久久精品久久久久真实原创| 国产精品国产三级专区第一集| 国产精品一区二区性色av| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 波野结衣二区三区在线| 伊人久久国产一区二区| 美女视频免费永久观看网站| 国产人妻一区二区三区在| 午夜福利视频1000在线观看| 亚洲精品影视一区二区三区av| 欧美精品国产亚洲| 成人毛片60女人毛片免费| 成年版毛片免费区| 亚洲人成网站在线观看播放| 观看免费一级毛片| 中文乱码字字幕精品一区二区三区| 精品少妇久久久久久888优播| 又粗又硬又长又爽又黄的视频| 国产精品偷伦视频观看了| 欧美亚洲 丝袜 人妻 在线| 两个人的视频大全免费| 一级毛片黄色毛片免费观看视频| 国产探花极品一区二区| 日韩电影二区| 国产日韩欧美亚洲二区| 如何舔出高潮| 午夜福利视频1000在线观看| 亚洲aⅴ乱码一区二区在线播放| 2018国产大陆天天弄谢| 日韩精品有码人妻一区| 男女啪啪激烈高潮av片| 午夜日本视频在线| 永久免费av网站大全| 麻豆成人av视频| 欧美潮喷喷水| 青春草视频在线免费观看| av国产久精品久网站免费入址| av免费在线看不卡| 国产乱人视频| 99热这里只有精品一区| av线在线观看网站| 大片免费播放器 马上看| 啦啦啦中文免费视频观看日本| 久久久久久久久大av| 国产午夜精品久久久久久一区二区三区| 亚洲国产日韩一区二区| 男女无遮挡免费网站观看| 午夜精品一区二区三区免费看| 免费在线观看成人毛片| 97热精品久久久久久| 男女边吃奶边做爰视频| 黄片wwwwww| 国产爱豆传媒在线观看| av福利片在线观看| 日本熟妇午夜| 高清在线视频一区二区三区| 少妇熟女欧美另类| 国产亚洲精品久久久com| 久久人人爽av亚洲精品天堂 | 伦精品一区二区三区| 亚洲av在线观看美女高潮| 久热久热在线精品观看| 99热全是精品| 人妻少妇偷人精品九色| 秋霞在线观看毛片| 亚洲精品日本国产第一区| 最后的刺客免费高清国语| 深爱激情五月婷婷| 精品一区二区免费观看| 国内精品美女久久久久久| 边亲边吃奶的免费视频| 午夜福利网站1000一区二区三区| 三级经典国产精品| 99久久精品一区二区三区| 涩涩av久久男人的天堂| 蜜桃亚洲精品一区二区三区| 一级黄片播放器| 狂野欧美激情性bbbbbb| 国产中年淑女户外野战色| 日韩免费高清中文字幕av| 少妇猛男粗大的猛烈进出视频 | 成人亚洲精品av一区二区| 国产 精品1| 一级黄片播放器| 久久人人爽人人片av| 国产中年淑女户外野战色| 亚洲精品日本国产第一区| 国产淫片久久久久久久久| 好男人视频免费观看在线| 日韩成人av中文字幕在线观看| 欧美国产精品一级二级三级 | 欧美3d第一页| 少妇人妻精品综合一区二区| 日韩免费高清中文字幕av| 女人十人毛片免费观看3o分钟| 舔av片在线| 美女内射精品一级片tv| 婷婷色av中文字幕| 高清毛片免费看| a级一级毛片免费在线观看| 九九在线视频观看精品| 成人黄色视频免费在线看| 国产一区亚洲一区在线观看| 我要看日韩黄色一级片| 亚洲四区av| 亚洲美女搞黄在线观看| 十八禁网站网址无遮挡 | 国产精品久久久久久久久免| 亚洲三级黄色毛片| 三级男女做爰猛烈吃奶摸视频| 午夜日本视频在线| 国产成人精品婷婷| 久久久久久久午夜电影| 日韩强制内射视频| 亚洲,一卡二卡三卡| 日韩欧美精品免费久久| 性色avwww在线观看| 国产美女午夜福利| 国产综合懂色| 极品少妇高潮喷水抽搐| 国产探花在线观看一区二区| 在线观看一区二区三区激情| 成年女人在线观看亚洲视频 | 超碰av人人做人人爽久久| 精品久久久精品久久久| 欧美丝袜亚洲另类| 亚洲精品国产av成人精品| 一区二区三区四区激情视频| 日韩av免费高清视频| 一级毛片 在线播放| 一级a做视频免费观看| 又大又黄又爽视频免费| 亚洲av免费在线观看| 少妇的逼好多水| 亚洲国产精品成人综合色| 啦啦啦中文免费视频观看日本| 一个人观看的视频www高清免费观看| 国产成人免费无遮挡视频| 欧美zozozo另类| 国模一区二区三区四区视频| 久久久色成人| 99热6这里只有精品| 国产 一区 欧美 日韩| 国产在线男女| 亚洲在久久综合| 十八禁网站网址无遮挡 | 久久精品久久久久久噜噜老黄| 街头女战士在线观看网站| 街头女战士在线观看网站| 久久精品久久久久久噜噜老黄| 国产精品一及| 午夜激情久久久久久久| 国产亚洲午夜精品一区二区久久 | 久久久a久久爽久久v久久| 午夜免费男女啪啪视频观看| 日韩在线高清观看一区二区三区| 欧美高清成人免费视频www| 午夜视频国产福利| 国产精品一区二区在线观看99| 午夜福利高清视频| 五月伊人婷婷丁香| 一本色道久久久久久精品综合| 亚洲av免费在线观看| 深爱激情五月婷婷| 能在线免费看毛片的网站| 亚洲人成网站在线播| 欧美一级a爱片免费观看看| 在线a可以看的网站| 亚洲精品久久久久久婷婷小说| 国产成人福利小说| 亚洲国产最新在线播放| 麻豆精品久久久久久蜜桃| 久久这里有精品视频免费| 成人高潮视频无遮挡免费网站| 三级国产精品欧美在线观看| 国产精品秋霞免费鲁丝片| 中国国产av一级| 嫩草影院入口| 国产精品嫩草影院av在线观看| 国产v大片淫在线免费观看| 国产在线一区二区三区精| 国产综合懂色| 人人妻人人看人人澡| 中文资源天堂在线| www.av在线官网国产| 一级毛片黄色毛片免费观看视频| 国产精品福利在线免费观看| 亚洲欧美成人综合另类久久久| 亚洲美女视频黄频| 亚洲怡红院男人天堂| 18禁裸乳无遮挡动漫免费视频 | av国产精品久久久久影院| 国精品久久久久久国模美| 卡戴珊不雅视频在线播放| 91精品伊人久久大香线蕉| 六月丁香七月| 成人无遮挡网站| 精品一区二区三区视频在线| 成人亚洲精品av一区二区| 国产探花在线观看一区二区| 日本熟妇午夜| 免费高清在线观看视频在线观看| 99热网站在线观看| 在线观看国产h片| 热99国产精品久久久久久7| a级毛色黄片| 另类亚洲欧美激情| 久久久久久九九精品二区国产| 亚洲国产精品专区欧美| 久久精品综合一区二区三区| 国产免费视频播放在线视频| 久久久久网色| 国产欧美日韩精品一区二区| 看十八女毛片水多多多| 日本午夜av视频| 国产av不卡久久| 国产免费视频播放在线视频| 久久久色成人| 亚洲精品中文字幕在线视频 | 99热这里只有是精品50| 国内精品宾馆在线| 国产精品一区www在线观看| 晚上一个人看的免费电影| 国产精品久久久久久精品古装| 欧美高清性xxxxhd video| 亚洲天堂国产精品一区在线| 精品一区二区三卡| 国产精品人妻久久久久久| 日韩成人伦理影院| 人妻制服诱惑在线中文字幕| 免费av不卡在线播放| 99精国产麻豆久久婷婷| 亚洲电影在线观看av| 欧美潮喷喷水| 成人二区视频| 欧美97在线视频| 全区人妻精品视频| av免费观看日本| 国产亚洲午夜精品一区二区久久 | 少妇的逼水好多| 亚洲经典国产精华液单| 欧美潮喷喷水| av福利片在线观看| 少妇被粗大猛烈的视频| 成人特级av手机在线观看| 亚洲国产最新在线播放| 国产精品国产三级国产专区5o| 九九爱精品视频在线观看| 亚洲欧美日韩东京热| 在线观看av片永久免费下载| 亚洲国产精品专区欧美| 亚洲精品影视一区二区三区av| 色视频www国产| 国产精品国产av在线观看| 午夜福利网站1000一区二区三区| 亚洲欧美成人综合另类久久久| 国产综合精华液| 国产片特级美女逼逼视频| 欧美高清性xxxxhd video| 91精品一卡2卡3卡4卡| 国产av码专区亚洲av| 只有这里有精品99| 久久久久久伊人网av| 午夜亚洲福利在线播放| 亚洲精品第二区| 久久这里有精品视频免费| 亚洲av二区三区四区| 午夜福利高清视频| 又大又黄又爽视频免费| 亚洲国产欧美人成| 亚洲精品久久午夜乱码| 亚洲精品第二区| 久久久久久久亚洲中文字幕| 九九在线视频观看精品| 久久国产乱子免费精品| 在线观看美女被高潮喷水网站| 国产精品99久久99久久久不卡 | 亚洲性久久影院| av网站免费在线观看视频| www.色视频.com| 日韩成人av中文字幕在线观看| 熟女人妻精品中文字幕| 又爽又黄a免费视频| 久久久久久久久久久免费av| 亚洲精品成人久久久久久| 国产成年人精品一区二区| 日韩成人av中文字幕在线观看| 亚洲欧美清纯卡通| 国产69精品久久久久777片| 乱码一卡2卡4卡精品| 在线播放无遮挡| 免费观看av网站的网址| 欧美日韩精品成人综合77777| 国产国拍精品亚洲av在线观看| 国产高清有码在线观看视频| 国产精品久久久久久精品电影小说 | 国产精品蜜桃在线观看| 日日啪夜夜爽| 人人妻人人爽人人添夜夜欢视频 | 你懂的网址亚洲精品在线观看| 国产黄色免费在线视频| 高清日韩中文字幕在线| 有码 亚洲区| 亚洲成人久久爱视频| 波多野结衣巨乳人妻| 国产免费福利视频在线观看| 久久久成人免费电影| 久久6这里有精品| 婷婷色av中文字幕| 高清视频免费观看一区二区| 国产精品国产三级专区第一集| 嫩草影院精品99| 亚洲欧美一区二区三区国产| 成年女人在线观看亚洲视频 | 欧美日韩视频高清一区二区三区二| 人妻少妇偷人精品九色| 日韩电影二区| 天堂中文最新版在线下载 | 综合色丁香网| 国产精品久久久久久精品电影小说 | 九色成人免费人妻av| 一本久久精品| 尤物成人国产欧美一区二区三区| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区黑人 | 国产午夜精品一二区理论片| 中文字幕久久专区| 精品亚洲乱码少妇综合久久| 亚洲人成网站高清观看| 亚洲三级黄色毛片| 午夜精品一区二区三区免费看| 国产高清三级在线| 精品视频人人做人人爽| 国产69精品久久久久777片| 日韩 亚洲 欧美在线| 熟女人妻精品中文字幕| 免费在线观看成人毛片| 久久久久性生活片| 晚上一个人看的免费电影| 老师上课跳d突然被开到最大视频| 成人午夜精彩视频在线观看| 国产成人午夜福利电影在线观看| 人人妻人人看人人澡| 夜夜爽夜夜爽视频| 亚洲av中文字字幕乱码综合| 国产高清不卡午夜福利| 国产精品无大码| 欧美少妇被猛烈插入视频| 亚洲va在线va天堂va国产| 亚洲国产精品国产精品| 插逼视频在线观看| 免费电影在线观看免费观看| 午夜亚洲福利在线播放| 水蜜桃什么品种好| 亚洲av成人精品一区久久| 久久久久久久久久久免费av| 国产午夜精品久久久久久一区二区三区| 丰满人妻一区二区三区视频av| 卡戴珊不雅视频在线播放| 欧美丝袜亚洲另类| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 在线观看三级黄色| 亚洲精品成人av观看孕妇| 亚洲欧美日韩东京热| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 精品久久久久久久人妻蜜臀av| 婷婷色麻豆天堂久久| a级毛片免费高清观看在线播放| 亚洲精品成人av观看孕妇| 成人高潮视频无遮挡免费网站| 欧美日韩一区二区视频在线观看视频在线 | 91狼人影院| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频| 国产乱人偷精品视频| 全区人妻精品视频| 五月玫瑰六月丁香| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 看非洲黑人一级黄片| 丰满少妇做爰视频| 国产乱人偷精品视频| 亚洲av国产av综合av卡| 亚洲国产精品成人综合色| 中文字幕制服av| 精品酒店卫生间| 精品一区二区免费观看| 18禁在线无遮挡免费观看视频| 丝袜喷水一区| 国产高潮美女av| 偷拍熟女少妇极品色| 国产精品伦人一区二区| www.av在线官网国产| 一区二区av电影网| 午夜福利视频1000在线观看| 欧美97在线视频| 亚洲欧美精品专区久久| 午夜福利视频精品| 男女无遮挡免费网站观看| 日本与韩国留学比较| 丰满乱子伦码专区| 亚洲精品成人av观看孕妇| 天堂网av新在线| 午夜福利网站1000一区二区三区| 国产男人的电影天堂91| 亚洲成人久久爱视频| 免费看不卡的av| 汤姆久久久久久久影院中文字幕| 97在线视频观看| 99热6这里只有精品| 国产精品不卡视频一区二区| 2021天堂中文幕一二区在线观| 久久久精品欧美日韩精品| 久久精品国产亚洲网站| 久久亚洲国产成人精品v| 亚洲欧美精品专区久久| 久久综合国产亚洲精品| 日本黄色片子视频| 亚洲国产日韩一区二区| 成人国产麻豆网| 18禁在线无遮挡免费观看视频| 新久久久久国产一级毛片| 精品国产一区二区三区久久久樱花 | 观看美女的网站| 国产一级毛片在线| 成年版毛片免费区| 精品久久久久久电影网| av福利片在线观看| 色网站视频免费| 亚洲av.av天堂| 亚洲精品日本国产第一区| 免费av观看视频| av免费观看日本| videossex国产| 久久久久久久亚洲中文字幕| av福利片在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 91狼人影院| 亚洲欧美成人综合另类久久久| 男插女下体视频免费在线播放| 日韩视频在线欧美| 2022亚洲国产成人精品| 偷拍熟女少妇极品色| 熟妇人妻不卡中文字幕| 亚洲电影在线观看av| 大码成人一级视频| 久久精品久久久久久久性| 国产精品.久久久| 高清日韩中文字幕在线| 日韩 亚洲 欧美在线| 免费看光身美女| 国产成人免费观看mmmm| 男人爽女人下面视频在线观看| 精品酒店卫生间| 中国国产av一级| 成年女人在线观看亚洲视频 | 少妇人妻精品综合一区二区| 亚洲,欧美,日韩| 免费电影在线观看免费观看| 午夜激情福利司机影院| 成人免费观看视频高清| 国产亚洲91精品色在线| 少妇熟女欧美另类| 女人久久www免费人成看片| 一本色道久久久久久精品综合| 亚洲欧美日韩另类电影网站 | 国产亚洲一区二区精品| 汤姆久久久久久久影院中文字幕| 51国产日韩欧美| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区三区| 蜜桃亚洲精品一区二区三区| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| 亚洲国产日韩一区二区| 亚洲国产欧美人成| videos熟女内射| 国产 一区精品| 麻豆成人午夜福利视频| 免费av不卡在线播放| 日本爱情动作片www.在线观看| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 99热这里只有精品一区| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 你懂的网址亚洲精品在线观看| 国产成人a∨麻豆精品| 中国美白少妇内射xxxbb| 一级片'在线观看视频| 亚洲久久久久久中文字幕| 蜜臀久久99精品久久宅男| 国产老妇伦熟女老妇高清| 波野结衣二区三区在线| 少妇裸体淫交视频免费看高清| 80岁老熟妇乱子伦牲交| 色吧在线观看| 欧美日韩亚洲高清精品| 国内少妇人妻偷人精品xxx网站| 久久午夜福利片| 成人免费观看视频高清| 丝袜美腿在线中文| 麻豆成人av视频| 国产亚洲91精品色在线| av在线老鸭窝| 男女下面进入的视频免费午夜| 国产成人aa在线观看| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 午夜精品一区二区三区免费看| 最新中文字幕久久久久| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 欧美激情在线99| 亚洲人成网站在线观看播放| 亚洲成人中文字幕在线播放| 男人爽女人下面视频在线观看| 伊人久久精品亚洲午夜| 亚洲国产欧美人成| 久久久久久久亚洲中文字幕| 精品人妻视频免费看| 春色校园在线视频观看| 波野结衣二区三区在线| 日日摸夜夜添夜夜爱| 人妻制服诱惑在线中文字幕| av在线天堂中文字幕| 成人一区二区视频在线观看| 国产精品一及| 国产v大片淫在线免费观看| 2022亚洲国产成人精品| 麻豆成人午夜福利视频| 身体一侧抽搐| 性插视频无遮挡在线免费观看| 秋霞伦理黄片| 亚洲av男天堂| 尤物成人国产欧美一区二区三区| 少妇人妻一区二区三区视频| 热99国产精品久久久久久7| 国产毛片a区久久久久| 综合色av麻豆| 国产欧美日韩精品一区二区| 777米奇影视久久| 免费少妇av软件| 啦啦啦在线观看免费高清www| 中文在线观看免费www的网站| 久久精品久久久久久噜噜老黄| 特大巨黑吊av在线直播| 亚洲欧洲国产日韩| 国产91av在线免费观看| 又大又黄又爽视频免费| 久久精品国产a三级三级三级| www.色视频.com| 欧美高清性xxxxhd video| 精品国产乱码久久久久久小说| 晚上一个人看的免费电影| 亚洲av成人精品一二三区| 99热这里只有精品一区| 少妇被粗大猛烈的视频| 香蕉精品网在线| 搡女人真爽免费视频火全软件| 97在线人人人人妻| 各种免费的搞黄视频| 久久久久久久亚洲中文字幕| 2018国产大陆天天弄谢| 午夜免费观看性视频| 嫩草影院精品99| 日韩,欧美,国产一区二区三区| 国产免费视频播放在线视频| 成人亚洲精品一区在线观看 | 久久热精品热| 精品一区二区三区视频在线| 国内精品美女久久久久久| av女优亚洲男人天堂| 交换朋友夫妻互换小说| 国产精品熟女久久久久浪| 成人亚洲精品一区在线观看 | 久久久成人免费电影| av线在线观看网站| 26uuu在线亚洲综合色| 蜜桃久久精品国产亚洲av| 高清日韩中文字幕在线| 欧美一级a爱片免费观看看| 亚洲国产欧美人成| 午夜免费男女啪啪视频观看| 精品久久久久久久久av| 国产午夜精品一二区理论片| 欧美日韩视频高清一区二区三区二| 黄色配什么色好看| 久久99热这里只频精品6学生| 人人妻人人澡人人爽人人夜夜| 国产亚洲精品久久久com| 在线看a的网站| 国产精品国产三级国产专区5o| 免费播放大片免费观看视频在线观看|