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

    一種基于二維Helmholtz腔陣列的低頻寬帶隔聲結(jié)構(gòu)實驗研究?

    2017-07-31 01:37:34高東寶劉選俊田章福周澤民曾新吾韓開鋒
    物理學(xué)報 2017年1期
    關(guān)鍵詞:共振頻率禁帶隔聲

    高東寶 劉選俊 田章福 周澤民 曾新吾 韓開鋒

    (國防科學(xué)技術(shù)大學(xué)海洋科學(xué)與工程研究院,長沙 410073)(2016年7月11日收到;2016年9月27日收到修改稿)

    一種基于二維Helmholtz腔陣列的低頻寬帶隔聲結(jié)構(gòu)實驗研究?

    高東寶?劉選俊 田章福 周澤民 曾新吾 韓開鋒

    (國防科學(xué)技術(shù)大學(xué)海洋科學(xué)與工程研究院,長沙 410073)(2016年7月11日收到;2016年9月27日收到修改稿)

    基于圓周排列的Helmholtz共振腔單元,設(shè)計并實現(xiàn)了一種具有低頻寬禁帶的聲人工結(jié)構(gòu),可以在結(jié)構(gòu)中心處實現(xiàn)二維隔聲效果.針對實際模型,搭建了二維聲場測量平臺,進(jìn)行了相應(yīng)的實驗研究,實驗結(jié)果與有限元仿真結(jié)果符合較好.該結(jié)構(gòu)在較寬的頻帶內(nèi)(680—1050 Hz)可以實現(xiàn)較好的隔聲效果,最大隔聲量可達(dá)41dB.實驗中還研究了單元參數(shù)及共振狀態(tài)對隔聲效果的影響.隔聲區(qū)的大小與共振單元的分布形式有直接關(guān)系,而良好的共振狀態(tài)將對提高隔聲量有一定幫助.研究結(jié)果對設(shè)計新型聲防護(hù)結(jié)構(gòu)具有理論與應(yīng)用價值.

    Helmholtz共振腔,低頻隔聲,二維聲場測量,聲防護(hù)

    1引 言

    聲人工結(jié)構(gòu)是聲學(xué)領(lǐng)域近年來發(fā)展起來的一個重要研究方向.一般按照所研究參數(shù)的不同,可以將聲人工結(jié)構(gòu)分為聲子晶體[1](主要關(guān)注結(jié)構(gòu)的聲傳輸特性)和聲學(xué)超材料[2](主要關(guān)注結(jié)構(gòu)的等效參數(shù)特性)兩大類.從功能上講,聲人工結(jié)構(gòu)的作用主要是控制聲波的傳播,從而實現(xiàn)各種超常的聲功能器件[3?7].

    聲人工結(jié)構(gòu)通常由周期排列的人工微結(jié)構(gòu)單元構(gòu)成,微結(jié)構(gòu)單元形式多樣,Helmholtz共振腔(Helmholtz resonator,HR)是其中常見的一種.在聲學(xué)領(lǐng)域,對于Helmholtz腔的研究由來已久[8?10],但隨著聲子晶體研究的深入,人們發(fā)現(xiàn)在周期排列的HR結(jié)構(gòu)共振頻率附近會出現(xiàn)性能較好的局域共振型聲禁帶[10],同時從材料學(xué)的角度研究表明,HR結(jié)構(gòu)在其共振頻率附近會表現(xiàn)出等效質(zhì)量密度和等效體積模量同時為負(fù)的特性[11?13],這也為解釋其局域共振禁帶機(jī)理提供了理論支撐.

    一維HR陣列由于結(jié)構(gòu)簡單,現(xiàn)象明顯且易于分析而受到了更多的關(guān)注.已有文獻(xiàn)中,除了對完美周期結(jié)構(gòu)[11?13]的聲傳輸特性及等效介質(zhì)參數(shù)進(jìn)行了系統(tǒng)研究之外,人們還對含缺陷態(tài)的一維結(jié)構(gòu)[14]和參數(shù)梯度分布的結(jié)構(gòu)[15]都有所研究.Cheng等還對共振單元間存在耦合模式的情況進(jìn)行了分析,結(jié)果表明此時會出現(xiàn)禁帶加寬[16]或產(chǎn)生新禁帶[17]等現(xiàn)象.二維HR結(jié)構(gòu)的研究目前還比較有限[10],但由于同屬局域共振型結(jié)構(gòu),因此其他類似聲子晶體和超材料研究的成果也可以對其提供一些支撐和參考[18].

    理論與實驗研究結(jié)果均證明HR結(jié)構(gòu)是一種具有較好局域共振特性的聲學(xué)超材料,可以對波長遠(yuǎn)大于自身尺寸的聲波(即低頻聲波)進(jìn)行傳播控制.由于是一種局域共振型聲人工結(jié)構(gòu),HR聲子晶體所產(chǎn)生的聲禁帶基本特性主要與單元本身的共振特性相關(guān)[12],即理論上只需要一個單元就可以產(chǎn)生局域共振型聲禁帶,單元間的排列方式不會對禁帶本質(zhì)產(chǎn)生影響.

    目前來看,關(guān)于HR結(jié)構(gòu)的研究還主要處于理論與實驗驗證階段,不過HR結(jié)構(gòu)良好的低頻聲能量局域效果為傳統(tǒng)低頻隔聲隔振問題提供了新的思路[17].由于受到質(zhì)量定律的限制,常規(guī)材料對低于1000 Hz的聲波通常隔聲效果較差.局域共振型超材料通過設(shè)計共振單元的基本參數(shù),可以調(diào)節(jié)聲禁帶的位置,從而實現(xiàn)對波長遠(yuǎn)大于其尺寸的聲波進(jìn)行控制,諸多研究成果有望在此方面獲得新的突破.本文正是基于Helmholtz共振腔結(jié)構(gòu)的局域共振特性設(shè)計了一種具有低頻寬禁帶的二維隔聲結(jié)構(gòu).針對模型特點,搭建了一套二維聲場測量系統(tǒng).通過理論分析與實驗研究,證明了該結(jié)構(gòu)良好的區(qū)域隔聲特性,并討論了隔聲效果與結(jié)構(gòu)參數(shù)之間的基本關(guān)系.由于該結(jié)構(gòu)是“嵌”于水平面之下的,對水平面之上的物體可以實現(xiàn)沒有遮擋的隔聲效果,方便對隔聲區(qū)域內(nèi)的重要設(shè)施或人員進(jìn)行低頻聲防護(hù).

    2模型設(shè)計與仿真結(jié)果

    所設(shè)計的結(jié)構(gòu)中包含了四圈呈圓周排列的HR單元,基本結(jié)構(gòu)如圖1所示.其中R1,R2,R3和R4分別表示各圈HR對應(yīng)的圓環(huán)的半徑,分別為50,100,150和200mm.這里,將圓環(huán)由內(nèi)至外分別編號為1,2,3和4號.共振腔的基本結(jié)構(gòu)如圖1中的附圖所示,其中a表示短管的半徑,d表示短管的長度,L表示腔體的半徑,而H表示腔體的長度,其共振角頻率可以表示為其中c0=343m/s為空氣介質(zhì)的聲速,deffect表示等效短管長度,S=πa2為短管橫截面積,V=πL2H為共振腔體積.在本文中設(shè)a=10mm,d=20mm和L=20mm為固定值,而H可根據(jù)需要進(jìn)行調(diào)節(jié),以形成不同共振頻率的HR單元.為了保證結(jié)構(gòu)的隔聲效果,需要盡量緊密地排列HR單元,所設(shè)計結(jié)構(gòu)中由內(nèi)至外各列HR個數(shù)分別為6,12,18和24個,單元間距約為52mm.

    當(dāng)一列HR單元共振頻率梯度變化時,會產(chǎn)生局域共振型禁帶的重疊和耦合,從而實現(xiàn)寬頻禁帶效果[19].因此,本文所設(shè)計的結(jié)構(gòu)中將腔體高度由內(nèi)至外分別設(shè)置為H1=40mm,H2=35mm,H3=H4=30mm,對應(yīng)共振頻率分別為710,758和819 Hz.之所以將最外側(cè)兩層設(shè)置為相同的參數(shù),是因為最外側(cè)單元所面對的入射波范圍較大,增加單元層數(shù)可以提高隔聲量.本結(jié)構(gòu)中各層單元按照越靠近內(nèi)側(cè)共振頻率越小,而越靠近外側(cè)共振頻率越大的方式設(shè)計.這樣設(shè)計的原因是已有一維結(jié)構(gòu)研究結(jié)果[14]表明,當(dāng)HR單元未共振時(即入射波頻率小于其共振頻率),能量主要集中在共振腔內(nèi);而當(dāng)HR處于已共振狀態(tài)時(即入射波頻率大于其共振頻率),入射波能量則主要分布在波導(dǎo)中.因此,在本結(jié)構(gòu)中,將共振頻率較大的置于外側(cè)可以在內(nèi)部單元處于共振狀態(tài)時也局域一定的能量,從而提高此時的隔聲量.數(shù)值模擬采用有限元軟件COMSOL Multiphysics進(jìn)行,網(wǎng)格劃分嚴(yán)格按照有限元法的基本原理,每個波長約為10個網(wǎng)格.在模擬過程中,聲波沿x正方向入射.

    圖1 基本結(jié)構(gòu)示意圖Fig.1.Schematic structure.

    圖2中繪制出了結(jié)構(gòu)中心點處的隔聲量隨頻率變化曲線.從圖中結(jié)果可以清晰地看到在680—1050 Hz范圍內(nèi)存在一個聲禁帶,最大隔聲量可以達(dá)到近80dB.中心點頻率對應(yīng)的波長約為單元間隔的8.4倍.圖2結(jié)果說明該結(jié)構(gòu)可以對較低頻率的寬頻帶聲波產(chǎn)生良好的隔聲效果.

    圖2 結(jié)構(gòu)中心點處隔聲量Fig.2.Insulation sound pressure level spectrum of the center position.

    這一禁帶是由三種共振頻率的HR產(chǎn)生的局域共振型聲禁帶疊加而成的,所以禁帶特性與單一共振頻率的局域共振型禁帶[12]相比略有差別.從圖2中可以看到,聲禁帶范圍內(nèi)曲線并不是非常光滑,尤其在其左半邊,存在多個數(shù)值波動頻率點.這主要是由兩方面原因造成的.一方面,某些谷值對應(yīng)了不同參數(shù)單元的共振頻率,如778,836和876 Hz處.另外一些波動則是由于多種共振特性的單元間相互耦合形成的,本文不做深入討論.

    通過對比公式計算出的共振頻率和數(shù)值模擬得到的共振頻率可以發(fā)現(xiàn),二者之間存在6%—9%的差異,這主要是因為在模擬過程中采用的是三維模型,存在壁面、空間等相互作用過程,因此與采用集總模型得到的理論公式結(jié)果之間會存在一定的差別.不過,總體來說,差別還是合理的,并且不影響問題的本質(zhì).

    圖3中給出了幾個共振頻點(778,836和876 Hz)處聲壓級分布情況,圖中深顏色代表低聲壓級,即隔聲量大的情況,而淺顏色代表高聲壓級,即隔聲量小的區(qū)域.

    由圖3結(jié)果可見,當(dāng)頻率為778 Hz時,低聲壓級區(qū)域主要集中在第1層圓環(huán)內(nèi)部,圓環(huán)外聲壓級很高,而此時正是第1層圓環(huán)處于共振狀態(tài);當(dāng)頻率為836 Hz時,低聲壓級區(qū)域主要集中在第2層圓環(huán)內(nèi)部,此時第2層圓環(huán)處于共振態(tài);而當(dāng)頻率為876 Hz時,則正好在第4層圓環(huán)內(nèi)部,此時第3,4層圓環(huán)處于共振態(tài).從圖中分布結(jié)果可見,不同共振態(tài)下隔聲區(qū)域內(nèi)聲場分布都較為均勻.結(jié)合圖2中結(jié)果可知,所設(shè)計結(jié)構(gòu)不但具有較大的隔聲量,而且在整個隔聲區(qū)域中聲場均勻分布,隔聲效果較好.實際應(yīng)用過程中,該結(jié)構(gòu)可以為某些重要設(shè)施提供一個“安靜的”的二維隔聲區(qū).

    為更清晰地獲得不同狀態(tài)下的聲場分布特性,從圖3中提取了各共振態(tài)下x和y坐標(biāo)方向上隔聲量的分布情況,如圖4所示.其中圖4(a)表示x方向情況,圖4(b)表示y方向情況.從圖中結(jié)果可見,三種共振態(tài)下,均在各方向上形成了隔聲效果較好的“聲禁區(qū)”.明顯可以發(fā)現(xiàn),不同共振頻率對應(yīng)了不同的禁區(qū)范圍.778 Hz時,“聲禁區(qū)”主要集中在(?0.05m,0.05m)的區(qū)間;836 Hz 時,“聲禁區(qū)”則擴(kuò)展到了(?0.1m,0.1m)的范圍;當(dāng)最外兩層單元共振,即876 Hz時,禁區(qū)范圍進(jìn)一步擴(kuò)大到(?0.15m,0.15m).這說明通過調(diào)整不同位置處單元的共振狀態(tài),可以獲得不同尺寸(面積)的隔聲區(qū).

    另外,在圖4中還可以發(fā)現(xiàn)兩個有趣的現(xiàn)象.兩幅圖中均存在多處非常窄的高隔聲量位置.通過分析它們的位置可以發(fā)現(xiàn),這些非常窄的高隔聲量點正位于兩個坐標(biāo)軸上處于共振態(tài)的單元附近.這主要是因為此時臨近的單元處于共振狀態(tài),較多的能量被局域在共振腔內(nèi)部,造成其臨近位置出現(xiàn)了能量的低點.而距離共振單元較遠(yuǎn)的地方,受到不同方向散射波的影響,聲壓級反而沒那么低.第二個現(xiàn)象則是在圖4(a)中可以明顯地看到聲場分布不對稱的情況.原因在于本文中入射波沿x正方向入射,經(jīng)過x負(fù)半軸上單元的作用和阻擋,自然進(jìn)入x正半軸的能量就小了很多.由于圖4(b)中y方向與入射波方向垂直,所以聲場分布要對稱得多.

    圖3 (網(wǎng)刊彩色)不同共振頻點處聲壓級分布 (a)778 Hz;(b)836 Hz;(c)876 HzFig.3. (color online)Acoustic fi eld distributions at 778 Hz(a),836 Hz(b)and 876 Hz(c).

    上述結(jié)果證明了本文所設(shè)計的二維Helmholtz共振腔結(jié)構(gòu)可以在較低頻率的寬帶內(nèi)(680—1050 Hz)實現(xiàn)區(qū)域性的隔聲效果,而且從隔聲場的特性來看,分布較為均勻,隔聲范圍具有可調(diào)性.

    圖4 共振狀態(tài)下x和y方向上隔聲量的分布Fig.4.(color online)Insulation sound pressure level distributions in x(a)and y(b)directions at di ff erent resonant states.

    3實驗設(shè)計與結(jié)果分析

    根據(jù)仿真結(jié)果設(shè)計了實驗?zāi)P?如圖5所示.模型采用厚度為20mm,半徑為250mm的高密度聚氯乙烯板作為基本材料,材料聲阻抗相對于背景介質(zhì)空氣來說要大得多,可以認(rèn)為壁面為硬邊界.實驗中利用五層薄板疊加,整體尺寸與仿真模型相同.為了形成所需參數(shù)的Helmholtz共振腔,在最上層薄板中開有直徑為20mm的孔,作為Helmholtz腔的短管.下面三層薄板開有直徑為40mm的孔,共同構(gòu)成了Helmholtz腔的共振腔.最下層板不開孔,作為擋板使用.由于腔體采用三層板疊加而成,最大深度可達(dá)到60mm.通過調(diào)整腔體深度,可分別形成不同共振頻率的共振腔單元,最小頻率可低至580 Hz左右.實驗中,每個單元尺寸可單獨調(diào)整,形成任意組合.該模型不但可以滿足本文實驗要求,還可以為其他研究工作提供支持.

    圖5 實驗?zāi)P虵ig.5.Experimental model.

    實驗中所搭建的測量系統(tǒng)如圖6所示.實驗系統(tǒng)中包含信號發(fā)生器、有源音箱、PULSE 3050信號采集系統(tǒng)、B&K 4192聲信號測量傳感器以及用于信號處理與分析的筆記本電腦等.在測量過程中,為了形成較好的二維平面波聲場,在模型上方添加了一個二維波導(dǎo).波導(dǎo)由上下兩塊硬塑料板(材料與模型相同)構(gòu)成,邊界處采用螺釘固定及定位.塑料板厚度為10mm,尺寸為1.5m×1.5m,相對于所測結(jié)構(gòu)要大得多.在波導(dǎo)的三個傳播邊界處都放置了吸聲海綿,避免了界面處反射波對聲場的影響.上層波導(dǎo)表面開有11×11個測量孔,孔間距離為50mm,保證了所測量頻段內(nèi)(400—1200 Hz)一個波長最少有6個測點,滿足測量需求.測量范圍完全覆蓋了隔聲結(jié)構(gòu)表面,可以獲得隔聲區(qū)域范圍內(nèi)的聲場分布情況.測量過程中,波導(dǎo)間距為50mm,根據(jù)矩形聲波導(dǎo)理論[20]可知,其截至頻率可達(dá)到3400 Hz,遠(yuǎn)大于實驗中隔聲結(jié)構(gòu)工作頻段.

    測量過程中,由信號發(fā)生器產(chǎn)生所需單頻信號通過有源音箱發(fā)射,音箱距離波導(dǎo)入口1m,因此可將波導(dǎo)入口處聲場近似認(rèn)為是平面波.波導(dǎo)中聲場通過傳感器進(jìn)行測量,經(jīng)過信號采集系統(tǒng)采集聲信號.同時,波導(dǎo)入口中央也會放置一個傳感器,其測量得到的聲壓級作為入射聲壓級.每一點聲壓級經(jīng)過多次測量取平均值作為最終結(jié)果.各測點隔聲量計算過程可表示為Tn=Σ(Pn?Pni)/N(n=1,2,···,121),其中Pn表示波導(dǎo)中第n號位置處測到的聲壓級,Pni表示第n號位置處波導(dǎo)入口聲壓級,N表示每一點的測量次數(shù).

    圖6 實驗測量系統(tǒng)示意圖Fig.6.Experimental con fi guration.

    圖7 (網(wǎng)刊彩色)中心點處隔聲量隨頻率變化Fig.7.(color online)Spectrum of insulation SPL at center position.

    圖7中繪出了波導(dǎo)中心點處的隔聲量隨頻率變化曲線,其中藍(lán)色實線表示有限元理論計算結(jié)果,而紅色實線表示實驗測量結(jié)果.通過結(jié)果對比可以發(fā)現(xiàn),實驗所得結(jié)果無論在禁帶位置還是在禁帶寬度方面與理論結(jié)果都符合得非常好.實驗中,同樣在大約680—1050 Hz范圍內(nèi)出現(xiàn)了一個聲禁帶,并且可以清晰地發(fā)現(xiàn)其在785,840和890 Hz處存在三處共振態(tài),對應(yīng)了三種幾何參數(shù)的Helmholtz腔共振頻率.實驗結(jié)果與理論結(jié)果之間在Helmholtz腔的共振頻率方面略有差異,這主要是由于實際模型在加工裝配過程中的誤差造成的.從數(shù)值上而言,這種差異還是可以接受的,因為最大的誤差也只有1.6%.不過,在禁帶的深度方面,實驗結(jié)果明顯小于理論結(jié)果.有限元計算結(jié)果中,最大隔聲量可以達(dá)到近80dB,而實驗中最大隔聲量才只有41dB,僅約為其一半.本文中由于實際結(jié)構(gòu)存在各種吸收與散射作用,造成了聲波與共振腔之間耦合減弱,因此并不能如理論計算那樣,使HR共振態(tài)完全激發(fā),所以能量局域程度較為有限,因而禁帶沒有那么深.

    總的來說,圖7結(jié)果依然可以表明本文中所設(shè)計結(jié)構(gòu)具有較好的寬頻帶隔聲效果,而且最大隔聲量達(dá)到了41dB.

    圖8中給出了三種共振頻率下波導(dǎo)中聲場的分布情況,其中圓形表示不同層單元所處的位置.由圖可見,不同共振頻率下,波導(dǎo)中隔聲區(qū)大小不同.通過對比可知,圖8(a)與圖3(a)中結(jié)果分布非常相近,低聲壓級區(qū)域只分布在第1層單元范圍內(nèi).圖8(c)中聲場分布情況與圖3(c)中情況也比較類似,隔聲區(qū)域分布較廣,幾乎在整個模型區(qū)域隔聲量均為25dB以上,隔聲效果較好.

    對比圖8(b)與圖3(b)可以發(fā)現(xiàn),二者之間存在一些明顯的差異.在圖3(b)中,低聲壓級區(qū)域主要集中在第2層單元圍繞的范圍內(nèi);而在圖8(b)中,這一區(qū)域則擴(kuò)展到了第3、第4層單元區(qū)域,甚至某些位置的聲壓級還要低一些.在圖8(b)中某些位置,如(0.1,0),(0.05,?0.1),(?0.05,0.1),(0.15,0),(0.05,0.2)和(0.15,0.15)等都存在明顯的低聲壓級狀態(tài).上述位置正對應(yīng)了不同層上共振單元的位置,如(0.1,0)和(0.05,?0.1)對應(yīng)于第2層的某些單元,(?0.05,0.1)和(0.15,0)對應(yīng)于第3層的某些單元,而(0.05,0.2)和(0.15,0.15)恰好對應(yīng)了第4層上的某些單元.結(jié)合圖4的結(jié)果可知,此時這些單元處于共振態(tài).然而此時入射波頻率為840 Hz,只等于第2層單元的共振頻率,而不等于第3、第4層單元的共振頻率,理論上這兩層的單元不應(yīng)該共振才對.不過,由于最外側(cè)兩層單元的共振頻率(890 Hz)與第2層單元的共振頻率(840 Hz)相差較小,同時某些單元在裝配過程中可能也造成了尺寸的偏差,從而導(dǎo)致這些單元共振頻率較為相近,所以當(dāng)入射波頻率在840—890 Hz頻段內(nèi),都會造成這些單元不同程度的起振.相似地,這一現(xiàn)象也可以在圖8(c)中發(fā)現(xiàn),此時不但在第3、第4層位置處的單元(?0.15,0),(?0.1,0.1)和(0.05,0.15)產(chǎn)生了共振,而且在第2層(0.05,?0.05)位置附近的某些單元也出現(xiàn)了共振現(xiàn)象.

    圖8 (網(wǎng)刊彩色)三種共振頻率對應(yīng)的聲場分布情況(a)785 Hz;(b)840 Hz;(c)890 HzFig.8.(color online)Acoustic fi eld distributions at frequencies of 785 Hz(a),840 Hz(b)and 890 Hz(c)for experimental results.

    裝配誤差所導(dǎo)致的共振頻率的差別不但增加了隔聲區(qū)域控制上的難度,而且也是造成圖7中隔聲量較低的另外一個重要原因.由于同一層的單元不能同時共振,在某些共振態(tài)較弱的位置便出現(xiàn)了“漏洞”,使入射聲波從這些位置傳播到了隔聲區(qū)域內(nèi)部,以至于降低了隔聲量.

    為了定量地說明該結(jié)構(gòu)的實際隔聲效果,圖9中分別給出了x和y軸上隔聲量的分布結(jié)果.整體來看,圖9中結(jié)果與圖4中結(jié)果相類似,“聲禁區(qū)”現(xiàn)象比較明顯.對于785 Hz的情況,“聲禁區(qū)”相對較小,明顯存在漏斗形分布情況,其他兩種頻率下“聲禁區(qū)”分布較廣,隔聲量分布曲線更為平緩.類似地,在圖9中,也能發(fā)現(xiàn)明顯的高隔聲量位置點,如圖9(a)中840 Hz時?0.1m和0.15m兩點,890 Hz時的?0.15m處等.結(jié)合圖8結(jié)果可知,此時該位置單元處于共振態(tài),這也直接證明了該結(jié)構(gòu)中隔聲區(qū)是由于共振單元出現(xiàn)了共振現(xiàn)象從而阻止了聲波的傳播而產(chǎn)生的.對比發(fā)現(xiàn),實驗結(jié)果(圖9)與理論結(jié)果(圖4)仍存在一定的差異,除了數(shù)值方面,實驗結(jié)果的分布曲線的對稱性明顯不如理論結(jié)果的好.這一原因與上文中的分析結(jié)果相似,主要還是因為某些單元的共振態(tài)出現(xiàn)了偏差,從而使結(jié)構(gòu)中聲場分布不均勻?qū)е碌?不過總體來說,實驗現(xiàn)象符合理論預(yù)期,證明了該結(jié)構(gòu)具有較好的二維隔聲效果.

    圖9 (網(wǎng)刊彩色)共振狀態(tài)下x和y方向上隔聲量的分布(實驗結(jié)果)Fig.9.(color online)Insulation sound pressure level distributions in x(a)and y(b)directions at di ff erent resonant states(experimental results).

    4結(jié) 論

    本文基于二維分布的Helmholtz共振腔陣列設(shè)計了一種具有低頻寬禁帶隔聲效果的人工結(jié)構(gòu).該結(jié)構(gòu)具有單元共振頻率可調(diào)的特性,可以通過幾何參數(shù)的優(yōu)化設(shè)計實現(xiàn)針對不同頻帶的二維區(qū)域隔聲.針對該模型,搭建了一套二維聲場測量系統(tǒng),進(jìn)行了理論與實驗研究,結(jié)果表明該結(jié)構(gòu)可在680—1050 Hz范圍內(nèi)實現(xiàn)有效隔聲,實驗中隔聲量最大可達(dá)到41dB.隔聲區(qū)內(nèi)聲場分布較為均勻,隔聲區(qū)大小與共振態(tài)HR的排列情況有關(guān),可通過調(diào)節(jié)不同位置處HR的共振頻率,改變隔聲區(qū)大小.由于實驗中裝配誤差的存在,導(dǎo)致隔聲量等參數(shù)受到了一定的影響,需要在后續(xù)的研究中進(jìn)一步增強(qiáng)單元的一致性改善隔聲效果.該結(jié)構(gòu)的另一個優(yōu)點是無需采用具有遮擋效果的隔聲材料便能實現(xiàn)區(qū)域性的隔聲效果,可以有效地對器件、裝備或人員進(jìn)行聲防護(hù).

    本文所分析的結(jié)構(gòu)只是一種原理性模型,在實際應(yīng)用中還需要根據(jù)具體問題來恰當(dāng)?shù)卦O(shè)計相應(yīng)參數(shù).如本文模型中最內(nèi)層隔聲區(qū)域直徑只有0.1m,相對較小,在實際應(yīng)用中只要隔聲區(qū)域的尺寸,即最內(nèi)層HR陣列所圍成的區(qū)域面積大于所需隔聲的目標(biāo)物體的尺寸即可.另外,隨著聲人工結(jié)構(gòu)技術(shù)的發(fā)展,研究人員也發(fā)現(xiàn)了工作在深亞波長尺度的超材料和單向傳輸超材料等[21?23],為低頻隔聲問題提供了更廣闊的思路.

    [1]Kushwaha M S,Halevi P,Dobrzynski L,Djafari-Rouhani B 1993 Phys.Rev.Lett.71 2022

    [2]Liu Z Y,Zhang X X,Mao Y W,Zhu Y Y,Yang Z Y,Chan C T,Sheng P 2000 Science 289 1734

    [3]Zhu R,Liu X N,Hu G K,Sun C T,Huang G L 2014 Nat.Commun.5 5510

    [4]Zhang S,Xia C G,Fang N 2011 Phys.Rev.Lett.106 024301

    [5]Zhu X F,Liang B,Kan W W,Zou X Y,Cheng J C 2011 Phys.Rev.Lett.106 014301

    [6]Kan W W,Garcia-Chocano V M,Cervera F,Liang B,Zou X Y,Yin L L,Cheng J C,Sanchez-Dehesa J 2015 Phys.Rev.Appl.3 064019

    [7]Wu L Y,Chen L W,Liu C M 2009 Appl.Phys.Lett.95 013506

    [8]Norris A N,Wickham G 1993 J.Acoust.Soc.Am.93 617

    [9]Sugimoto N,Horioka T 1995 J.Acoust.Soc.Am.97 1446

    [10]Hu X H,Chan C T 2005 Phys.Rev.E 71 055601

    [11]Fang N,Xi D J,Xu J Y,Ambati M,Srituravanich W,Sun C,Zhang X 2006 Nat.Mater.5 452

    [12]Cheng Y,Xu J Y,Liu X 2008 J.Phys.Rev.B 77 045134

    [13]Ding C L,Zhao X P 2009 Acta Phys.Sin.58 6351(in Chinese)[丁昌林,趙曉鵬 2009物理學(xué)報 58 6351]

    [14]Gao D B,Zeng X W,Zhou Z M,Tian Z F 2013 Acta Phys.Sin.62 094304(in Chinese)[高東寶,曾新吾,周澤民,田章福2013物理學(xué)報62 094304]

    [15]Gao D B,Zeng X W,Zhou Z M,Tian Z F 2014 Proceedings of the 21th Inter.Congress on Sound and Vib.Beijing,China,July 13–17,2014p2679

    [16]Cheng Y,Xu J Y,Liu X J 2008 Appl.Phys.Lett.92 051913

    [17]Cheng Y,Liu X J 2012 Appl.Phys.A 109 805

    [18]Liu J,Hou Z L,Fu X J 2015 Acta Phys.Sin.64 154302(in Chinese)[劉嬌,侯志林,傅秀軍 2015物理學(xué)報 64 154302]

    [19]Fey J,Robertson M 2011 J.Appl.Phys.109 114903

    [20]Du G H,Zhu Z M,Gong X F 2001 Fundamentals of Acoustics(Nanjing:Nanjing University Press)pp279–283(in Chinese)[杜功煥,朱哲民,龔秀芬 2001聲學(xué)基礎(chǔ)(南京:南京大學(xué)出版社)第279—283頁]

    [21]Zhu X F,Liang B,Kan W W,Peng Y G,Cheng J C 2016 Phys.Rev.Appl.5 054015

    [22]Zhu X F,Li K,Zhang P,Zhu J,Zhang J T,Tian C,Liu S C 2016 Nat.Commun.7 11731

    [23]Zhu X F,Ramezani H,Shi C Z,Zhu J,Zhang X 2014 Phys.Rev.X 4 031042

    PACS:43.40.+s,43.35.Gk,43.20.+gDOI:10.7498/aps.66.014307

    *Project supported by the National Natural Science Foundation of China(Grant Nos.11504425,41374005).

    ?Corresponding author.E-mail:gaodongbao@nudt.edu.cn

    A broadband low-frequency sound insulation structure based on two-dimensionally inbuilt Helmholtz resonator?

    Gao Dong-Bao?Liu Xuan-Jun Tian Zhang-Fu Zhou Ze-Min Zeng Xin-Wu Han Kai-Feng

    (Academy of Ocean Science and Engineering,National University of Defense Technology,Changsha 410073,China)(Received 11 July 2016;revised manuscript received 27 September 2016)

    Helmholtz resonator(HR)has already been demonstrated both theoretically and experimentally to be a metamaterial with negative mass density and negative bulk modulus simultaneously.The HR can resonate at a frequency corresponding to a wavelength much longer than its geometrical parameters.At this time,the incident acoustic energy can be located.Therefore,the HR structures are considered to be good choices for controlling low-frequency sound waves.Furthermore,existing results indicate that the wide forbidden band could be formed by a one-dimensional structure shunted with detuned HRs.Based on these aforementioned theories,a man-made acoustical structure with broadband low-frequency sound insulation e ff ect is designed by circularly inbuilt HRs.Beyond this structure’s surface,a two-dimensional quiet zone can be created.With the same simulated model,an experimental structure is fabricated based on PVC plastic material.The structure consists of fi ve layerd circular plates.In the top four plates,two kinds of holes are drilled.The smaller holes in the top plate act as shot necks of the HR,while the bigger holes in the middle three plates serve as the cavities of the HR.They can construct 60resonators with di ff erent resonant frequencies.Experiments are carried out to study its sound insulation properties.In the experiments,three kinds of HRs with resonant frequencies 785,840and 890 Hz from inner loop to outer loop,respectively,are formed.The experimental results are very coincident with the simulation results from the software of COMSOL Multiphysics based on fi nite element method,which shows that this structure has an excellent sound insulation e ff ect in a frequency band of 680–1050 Hz,and the maximum insulation sound pressure level can reach 41dB.Meanwhile,the distribution of the two-dimensional sound fi eld is measured.The results point out that the range of the insulation area can be changed with the incident frequency.In addition,the sound insulation e ff ect is sensitive to the resonant state of the HRs.When all of the resonators at the same loop resonate simultaneously,the insulation sound pressure level will be higher.On the contrary,the insulation sound pressure level will be lower because of the energy leaking through the positions where the HRs do not resonate with the others.This work will be of help for designing new sound protection devices for low-frequency sound waves.

    Helmholtz resonator,low-frequency insulation,two-dimensional acoustic fi eld distribution,sound protection

    10.7498/aps.66.014307

    ?國家自然科學(xué)基金(批準(zhǔn)號:11504425,41374005)資助的課題.

    ?通信作者.E-mail:gaodongbao@nudt.edu.cn

    猜你喜歡
    共振頻率禁帶隔聲
    多激勵下結(jié)構(gòu)反共振頻率的計算和配置
    壓電周期板中耦合禁帶影響規(guī)律分析
    常用樓板隔聲材料的撞擊聲隔聲性能測試分析
    綠色建筑(2021年1期)2021-07-19 02:15:04
    一種隔聲軟簾的研制及應(yīng)用
    Multi-functional roles of TaSSI2 involved in Fusarium head blight and powdery mildew resistance and drought tolerance
    寬絕對禁帶的一維磁性光子晶體結(jié)構(gòu)?
    三層軟磁鎳超晶格薄膜的共振頻率*
    基于介質(zhì)圓柱結(jié)構(gòu)的復(fù)式二維光子晶體禁帶研究
    寬禁帶固態(tài)功放在導(dǎo)航衛(wèi)星中的應(yīng)用
    基于混響室—半消聲室的V型聲屏障隔聲測試與分析
    亚洲一区高清亚洲精品| 免费观看精品视频网站| 小蜜桃在线观看免费完整版高清| 中文字幕av在线有码专区| 国产精品乱码一区二三区的特点| 亚洲av二区三区四区| 国产高清不卡午夜福利| 亚洲欧美成人精品一区二区| 日韩 亚洲 欧美在线| 亚洲中文日韩欧美视频| 99久国产av精品国产电影| 国产亚洲av嫩草精品影院| 久久6这里有精品| 国产一区二区亚洲精品在线观看| 久久这里只有精品中国| 三级经典国产精品| av专区在线播放| 中文字幕免费在线视频6| 三级毛片av免费| 久久6这里有精品| 亚洲精品国产成人久久av| 国产精品一区二区三区四区久久| 我的女老师完整版在线观看| 在线观看一区二区三区| 国产av在哪里看| ponron亚洲| 亚洲一区二区三区色噜噜| av福利片在线观看| 亚洲成人av在线免费| 国产成人freesex在线 | 在线观看av片永久免费下载| 日韩欧美精品v在线| 国内揄拍国产精品人妻在线| 波多野结衣高清作品| 在线观看av片永久免费下载| 一级av片app| 久久精品国产自在天天线| 免费看av在线观看网站| 国产真实乱freesex| 欧美xxxx黑人xx丫x性爽| 久久久久久国产a免费观看| 国产精品,欧美在线| 亚洲精品日韩av片在线观看| 亚洲第一区二区三区不卡| 亚洲经典国产精华液单| 日韩,欧美,国产一区二区三区 | 国产伦一二天堂av在线观看| 成人精品一区二区免费| 国产精品电影一区二区三区| 日韩大尺度精品在线看网址| 国产在线精品亚洲第一网站| 国产高清三级在线| 又黄又爽又刺激的免费视频.| 欧美色欧美亚洲另类二区| 久久午夜亚洲精品久久| 亚洲精品国产av成人精品 | 九九爱精品视频在线观看| 直男gayav资源| 午夜日韩欧美国产| 日本黄色视频三级网站网址| 国产一区二区三区av在线 | 亚洲熟妇中文字幕五十中出| 一夜夜www| 国产真实乱freesex| 国产久久久一区二区三区| 久久久欧美国产精品| 99视频精品全部免费 在线| 自拍偷自拍亚洲精品老妇| 亚洲18禁久久av| 狠狠狠狠99中文字幕| 1024手机看黄色片| 免费黄网站久久成人精品| 在线国产一区二区在线| 精品熟女少妇av免费看| 大型黄色视频在线免费观看| 免费电影在线观看免费观看| av在线天堂中文字幕| aaaaa片日本免费| 久久久欧美国产精品| 久久久午夜欧美精品| 日韩强制内射视频| 狂野欧美白嫩少妇大欣赏| 国产在线男女| 嫩草影视91久久| 国产高清视频在线观看网站| 超碰av人人做人人爽久久| 成年女人看的毛片在线观看| 啦啦啦韩国在线观看视频| 精品久久久久久久久久免费视频| 国产大屁股一区二区在线视频| 蜜桃久久精品国产亚洲av| 日本在线视频免费播放| 亚洲av电影不卡..在线观看| 亚洲国产色片| 欧美丝袜亚洲另类| 男女视频在线观看网站免费| 夜夜看夜夜爽夜夜摸| 免费一级毛片在线播放高清视频| 波多野结衣巨乳人妻| 在线播放无遮挡| 内射极品少妇av片p| 三级经典国产精品| 看片在线看免费视频| 一本久久中文字幕| 国产日本99.免费观看| 特大巨黑吊av在线直播| 天堂网av新在线| 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 我的老师免费观看完整版| 深夜a级毛片| 精品国产三级普通话版| 亚洲国产欧洲综合997久久,| 久久综合国产亚洲精品| 18禁在线无遮挡免费观看视频 | 男女之事视频高清在线观看| 免费观看的影片在线观看| 日韩欧美精品免费久久| 一区二区三区四区激情视频 | 亚洲国产精品成人久久小说 | 日本在线视频免费播放| 色综合色国产| 婷婷色综合大香蕉| 观看美女的网站| 夜夜夜夜夜久久久久| 亚洲av成人精品一区久久| 一进一出抽搐动态| 十八禁国产超污无遮挡网站| 国产黄色视频一区二区在线观看 | 国产成人福利小说| 卡戴珊不雅视频在线播放| 老司机福利观看| 午夜日韩欧美国产| 亚洲电影在线观看av| 最新中文字幕久久久久| 深爱激情五月婷婷| 美女内射精品一级片tv| 欧美在线一区亚洲| 最近2019中文字幕mv第一页| 精品99又大又爽又粗少妇毛片| 熟妇人妻久久中文字幕3abv| 男人舔奶头视频| 色综合站精品国产| 六月丁香七月| 久99久视频精品免费| 国产又黄又爽又无遮挡在线| 亚洲美女搞黄在线观看 | 精品熟女少妇av免费看| 国产午夜福利久久久久久| 久久热精品热| 日韩一本色道免费dvd| 99久久九九国产精品国产免费| 欧美人与善性xxx| 日本成人三级电影网站| 国产一级毛片七仙女欲春2| 日韩欧美精品v在线| 最近视频中文字幕2019在线8| 1024手机看黄色片| 国产成人a区在线观看| 亚洲欧美精品自产自拍| 黄色一级大片看看| 日韩精品青青久久久久久| 久久久午夜欧美精品| 99视频精品全部免费 在线| 亚洲欧美成人综合另类久久久 | videossex国产| 中出人妻视频一区二区| 伦精品一区二区三区| 人人妻,人人澡人人爽秒播| 国产黄a三级三级三级人| 久久精品人妻少妇| 精华霜和精华液先用哪个| 国产精品精品国产色婷婷| 久久天躁狠狠躁夜夜2o2o| 美女cb高潮喷水在线观看| 欧美3d第一页| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 99国产极品粉嫩在线观看| 蜜桃久久精品国产亚洲av| 精品久久久久久久久av| 变态另类丝袜制服| 狠狠狠狠99中文字幕| 18禁裸乳无遮挡免费网站照片| 午夜老司机福利剧场| 99久国产av精品国产电影| 69人妻影院| 深爱激情五月婷婷| 亚洲无线在线观看| а√天堂www在线а√下载| 久久6这里有精品| 亚洲av第一区精品v没综合| 亚洲va在线va天堂va国产| 嫩草影院入口| 淫秽高清视频在线观看| 最近视频中文字幕2019在线8| 日本爱情动作片www.在线观看 | 狂野欧美激情性xxxx在线观看| 久久国内精品自在自线图片| 成人美女网站在线观看视频| 九九在线视频观看精品| 丰满的人妻完整版| 99久国产av精品国产电影| 欧美一级a爱片免费观看看| 午夜福利在线观看免费完整高清在 | 好男人在线观看高清免费视频| 男女那种视频在线观看| 国产免费一级a男人的天堂| 国产高清视频在线观看网站| 一本一本综合久久| 亚洲人成网站高清观看| 免费看光身美女| www日本黄色视频网| 成人午夜高清在线视频| 亚洲欧美精品综合久久99| 婷婷亚洲欧美| 国产成人91sexporn| 国产精品三级大全| 国产精品久久视频播放| 真实男女啪啪啪动态图| 国语自产精品视频在线第100页| 亚洲美女视频黄频| 国模一区二区三区四区视频| 精品人妻一区二区三区麻豆 | 国产精品免费一区二区三区在线| 国产黄色视频一区二区在线观看 | 村上凉子中文字幕在线| 成人午夜高清在线视频| 日韩一区二区视频免费看| 美女免费视频网站| 老司机影院成人| 欧美激情久久久久久爽电影| 国产成人福利小说| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕 | 午夜久久久久精精品| 亚洲精品在线观看二区| 18禁裸乳无遮挡免费网站照片| av在线亚洲专区| 我的女老师完整版在线观看| 国产v大片淫在线免费观看| 成人av一区二区三区在线看| 色播亚洲综合网| 国产精品野战在线观看| 一进一出好大好爽视频| 老司机午夜福利在线观看视频| 免费观看在线日韩| 一区二区三区四区激情视频 | 亚洲精品在线观看二区| 一本一本综合久久| 在线观看66精品国产| 久久久精品欧美日韩精品| 男女那种视频在线观看| 日本 av在线| 午夜激情福利司机影院| 亚洲欧美精品自产自拍| 亚洲精品影视一区二区三区av| 人人妻人人看人人澡| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类| 国产免费一级a男人的天堂| 一个人观看的视频www高清免费观看| 久久久久久久久中文| 97超视频在线观看视频| 晚上一个人看的免费电影| 久久久久久九九精品二区国产| 国模一区二区三区四区视频| 国内精品一区二区在线观看| 亚洲精品国产成人久久av| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 亚洲av免费高清在线观看| 国产白丝娇喘喷水9色精品| 亚洲成人中文字幕在线播放| 最近的中文字幕免费完整| 在线观看免费视频日本深夜| 91久久精品国产一区二区三区| 久久精品国产99精品国产亚洲性色| 12—13女人毛片做爰片一| 91久久精品国产一区二区成人| 欧美潮喷喷水| 国产精品国产三级国产av玫瑰| 午夜福利成人在线免费观看| 色在线成人网| 亚洲国产精品成人久久小说 | 亚洲综合色惰| 女的被弄到高潮叫床怎么办| 人人妻人人澡欧美一区二区| 亚洲国产色片| 午夜福利在线观看免费完整高清在 | 欧美成人一区二区免费高清观看| 黄色配什么色好看| 日本免费a在线| 真人做人爱边吃奶动态| 免费不卡的大黄色大毛片视频在线观看 | 全区人妻精品视频| 有码 亚洲区| 又爽又黄无遮挡网站| 人妻制服诱惑在线中文字幕| 精品国产三级普通话版| 亚洲七黄色美女视频| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲最大成人手机在线| 欧美潮喷喷水| 亚洲精品日韩av片在线观看| 国产精品久久视频播放| 精品国内亚洲2022精品成人| 免费看a级黄色片| 亚洲av熟女| 日韩欧美国产在线观看| 麻豆久久精品国产亚洲av| 舔av片在线| 别揉我奶头 嗯啊视频| 精品熟女少妇av免费看| 国产精品爽爽va在线观看网站| 亚洲av不卡在线观看| 午夜福利在线在线| 亚洲人与动物交配视频| 老熟妇仑乱视频hdxx| 最近手机中文字幕大全| 中文亚洲av片在线观看爽| 91久久精品国产一区二区成人| 欧美激情国产日韩精品一区| 一级毛片久久久久久久久女| 国产精品人妻久久久久久| 女生性感内裤真人,穿戴方法视频| 亚洲在线观看片| 免费av毛片视频| 亚洲不卡免费看| 看非洲黑人一级黄片| 欧美色欧美亚洲另类二区| 亚洲激情五月婷婷啪啪| 人妻少妇偷人精品九色| 1000部很黄的大片| 亚洲国产精品成人久久小说 | 国产精品人妻久久久影院| 国产乱人偷精品视频| 草草在线视频免费看| 精品午夜福利视频在线观看一区| 国产一区二区三区av在线 | 特级一级黄色大片| 日韩欧美免费精品| 国产精品99久久久久久久久| av在线天堂中文字幕| 亚洲精华国产精华液的使用体验 | 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 一边摸一边抽搐一进一小说| 欧美区成人在线视频| 91久久精品国产一区二区成人| 亚洲欧美日韩东京热| 我要搜黄色片| 国产精品乱码一区二三区的特点| 狠狠狠狠99中文字幕| 伊人久久精品亚洲午夜| 一夜夜www| 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 天天一区二区日本电影三级| 日本爱情动作片www.在线观看 | 天天躁日日操中文字幕| 成年版毛片免费区| 免费搜索国产男女视频| 日本熟妇午夜| 少妇猛男粗大的猛烈进出视频 | 国产成人91sexporn| 色尼玛亚洲综合影院| 22中文网久久字幕| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一区二区视频在线观看视频在线 | 少妇裸体淫交视频免费看高清| 欧美高清成人免费视频www| 天堂√8在线中文| 成人一区二区视频在线观看| 丰满乱子伦码专区| 激情 狠狠 欧美| 小蜜桃在线观看免费完整版高清| 中国国产av一级| 看黄色毛片网站| 久久久久久国产a免费观看| 看片在线看免费视频| 天堂动漫精品| 可以在线观看毛片的网站| 国产精品三级大全| 小蜜桃在线观看免费完整版高清| 波多野结衣高清作品| 十八禁国产超污无遮挡网站| 99在线人妻在线中文字幕| 久久久精品大字幕| 少妇的逼好多水| 啦啦啦观看免费观看视频高清| 国产精品一区二区三区四区久久| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| 中文资源天堂在线| 日本免费a在线| 成人高潮视频无遮挡免费网站| 69av精品久久久久久| 精品福利观看| 日本一二三区视频观看| 我的老师免费观看完整版| 午夜爱爱视频在线播放| 日韩精品青青久久久久久| 99热这里只有是精品在线观看| 国产高潮美女av| 午夜激情欧美在线| 熟女电影av网| 六月丁香七月| 国产伦精品一区二区三区四那| 国产精品99久久久久久久久| 少妇的逼好多水| 天堂影院成人在线观看| 久久精品久久久久久噜噜老黄 | 女生性感内裤真人,穿戴方法视频| 最近最新中文字幕大全电影3| 欧美一区二区国产精品久久精品| 久久天躁狠狠躁夜夜2o2o| 国内精品一区二区在线观看| 成人漫画全彩无遮挡| 一级毛片aaaaaa免费看小| 黑人高潮一二区| 国产精品一区二区免费欧美| 久久精品国产清高在天天线| 内地一区二区视频在线| 五月伊人婷婷丁香| 一进一出抽搐动态| 亚洲电影在线观看av| 亚洲五月天丁香| 亚洲最大成人中文| 亚洲av一区综合| 久久亚洲国产成人精品v| a级毛片免费高清观看在线播放| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 特大巨黑吊av在线直播| av在线观看视频网站免费| 乱系列少妇在线播放| 熟女电影av网| 国产成人a∨麻豆精品| 少妇熟女aⅴ在线视频| 少妇裸体淫交视频免费看高清| 在线播放无遮挡| 国产一区二区在线av高清观看| 中文在线观看免费www的网站| 日韩欧美一区二区三区在线观看| 久久精品国产99精品国产亚洲性色| 亚洲激情五月婷婷啪啪| 国产亚洲精品久久久com| 亚洲精品一卡2卡三卡4卡5卡| 69av精品久久久久久| 看黄色毛片网站| a级毛片免费高清观看在线播放| 人妻丰满熟妇av一区二区三区| 一级a爱片免费观看的视频| 国产日本99.免费观看| 亚洲精品日韩在线中文字幕 | 性欧美人与动物交配| 男人舔奶头视频| 能在线免费观看的黄片| 国产精品久久久久久久电影| 少妇裸体淫交视频免费看高清| 午夜视频国产福利| 国产黄a三级三级三级人| 国产精品一区二区免费欧美| 日韩精品有码人妻一区| 特大巨黑吊av在线直播| 精品人妻视频免费看| 国产单亲对白刺激| 人人妻人人澡欧美一区二区| 国产人妻一区二区三区在| 亚洲欧美精品自产自拍| 一边摸一边抽搐一进一小说| 国产精品一区二区三区四区免费观看 | 男插女下体视频免费在线播放| 国产精品一区二区三区四区免费观看 | 在线看三级毛片| 国产男靠女视频免费网站| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 亚洲精品久久国产高清桃花| 黄色日韩在线| 91在线精品国自产拍蜜月| 久久亚洲国产成人精品v| 久久久久久久午夜电影| 欧美潮喷喷水| 亚洲欧美成人精品一区二区| 99热网站在线观看| 色尼玛亚洲综合影院| 少妇丰满av| 成年女人毛片免费观看观看9| 日本免费一区二区三区高清不卡| 久久久久九九精品影院| 看黄色毛片网站| 变态另类丝袜制服| 欧美中文日本在线观看视频| 国产免费一级a男人的天堂| 久久午夜福利片| 黄色一级大片看看| 日本免费a在线| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 免费观看精品视频网站| 国产黄a三级三级三级人| 亚洲av中文字字幕乱码综合| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 变态另类丝袜制服| 可以在线观看毛片的网站| АⅤ资源中文在线天堂| 国产色婷婷99| 免费高清视频大片| 亚洲丝袜综合中文字幕| 久久精品综合一区二区三区| 少妇的逼水好多| 国产精品精品国产色婷婷| 日日摸夜夜添夜夜添av毛片| 男女之事视频高清在线观看| 国产白丝娇喘喷水9色精品| 国产精华一区二区三区| 麻豆久久精品国产亚洲av| 美女 人体艺术 gogo| av天堂在线播放| 亚洲精品456在线播放app| 久久久久久久久大av| 国产在线男女| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 欧美高清性xxxxhd video| 久久精品综合一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx性猛交bbbb| 国产高清激情床上av| 成人精品一区二区免费| videossex国产| 国产男靠女视频免费网站| 国产aⅴ精品一区二区三区波| 偷拍熟女少妇极品色| 亚洲电影在线观看av| 国产单亲对白刺激| 欧美精品国产亚洲| 亚洲欧美日韩东京热| 亚洲最大成人中文| 亚洲av免费高清在线观看| 美女免费视频网站| 欧美潮喷喷水| 天天躁日日操中文字幕| 婷婷精品国产亚洲av在线| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品电影| 亚洲欧美成人精品一区二区| 男插女下体视频免费在线播放| 亚洲欧美日韩无卡精品| 中文字幕久久专区| 免费搜索国产男女视频| 亚洲精品日韩在线中文字幕 | 久久久久久久久中文| 日韩欧美精品免费久久| 22中文网久久字幕| 精品久久久久久久久久久久久| a级毛片a级免费在线| 在线观看免费视频日本深夜| 性欧美人与动物交配| 国产av不卡久久| h日本视频在线播放| 在线看三级毛片| 免费电影在线观看免费观看| 伊人久久精品亚洲午夜| 乱系列少妇在线播放| 色综合亚洲欧美另类图片| 色5月婷婷丁香| 欧美色视频一区免费| 一区福利在线观看| 日韩在线高清观看一区二区三区| 禁无遮挡网站| 久久精品久久久久久噜噜老黄 | 欧美最黄视频在线播放免费| 日本a在线网址| 国产爱豆传媒在线观看| 尤物成人国产欧美一区二区三区| 综合色丁香网| 22中文网久久字幕| 中文在线观看免费www的网站| 日本免费a在线| 22中文网久久字幕| 乱码一卡2卡4卡精品| 免费大片18禁| 99久久久亚洲精品蜜臀av| 日本三级黄在线观看| 国产精品精品国产色婷婷| 麻豆久久精品国产亚洲av| 色尼玛亚洲综合影院| 亚洲av一区综合| 国产爱豆传媒在线观看| 老师上课跳d突然被开到最大视频| 中国国产av一级| 人妻夜夜爽99麻豆av| 亚洲第一电影网av| 国产精品女同一区二区软件| 成人av一区二区三区在线看| 亚洲精品乱码久久久v下载方式| 久久久精品94久久精品| 成年女人看的毛片在线观看| 蜜桃久久精品国产亚洲av| 久久综合国产亚洲精品| ponron亚洲| 最近最新中文字幕大全电影3| 女人十人毛片免费观看3o分钟| 国产精品久久久久久亚洲av鲁大| 精品一区二区三区av网在线观看| av卡一久久| 国产 一区精品|