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

    貼附型局域共振聲子晶體雙層板的帶隙特性*

    2019-06-25 06:32:56錢登輝史治宇吳靜紅
    振動(dòng)、測試與診斷 2019年3期
    關(guān)鍵詞:帶隙板結(jié)構(gòu)聲子

    錢登輝, 史治宇, 吳靜紅

    (1.蘇州科技大學(xué)土木工程學(xué)院 蘇州,215011) (2. 南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國家重點(diǎn)實(shí)驗(yàn)室 南京,210016)

    引 言

    雙層板以其質(zhì)量輕、剛度大以及抗沖擊性能好等優(yōu)點(diǎn)而被廣泛應(yīng)用于潛艇、飛機(jī)、汽車以及大型發(fā)電機(jī)組等工業(yè)產(chǎn)品的外殼結(jié)構(gòu)[1-2]。振動(dòng)的傳播大多以固體結(jié)構(gòu)為主要介質(zhì),其傳播距離遠(yuǎn)、衰減小,且在傳播過程中不斷向外輻射噪聲。以飛機(jī)艙室為例,發(fā)動(dòng)機(jī)振動(dòng)和空氣動(dòng)力等引起的結(jié)構(gòu)振動(dòng)傳播至艙室周邊的結(jié)構(gòu)引起其振動(dòng),繼而向艙室內(nèi)輻射噪聲,從而給乘客的舒適性帶來一定的影響。聲子晶體概念的提出為結(jié)構(gòu)減振降噪的理論研究和結(jié)構(gòu)設(shè)計(jì)開辟了新的研究方向。聲子晶體是一種具有彈性波帶隙特性的人工周期材料,在帶隙頻率范圍內(nèi)的振動(dòng)無法通過聲子晶體傳播,所以在減振降噪領(lǐng)域具備廣闊的應(yīng)用前景。過去幾十年來,聲子晶體的研究大多集中在帶隙計(jì)算方法和特性研究。Bragg散射[3-5]和局域共振[3,6-8]是帶隙能夠被打開的兩種主要機(jī)理,前者的帶隙所對(duì)應(yīng)頻段比后者高兩個(gè)數(shù)量級(jí)[6],這就意味著只有局域共振型聲子晶體才可以通過較小的周期尺寸來抑制較低頻率的振動(dòng)傳播。因此,通過將傳統(tǒng)局域共振型聲子晶體的設(shè)計(jì)思路引入到雙層板中構(gòu)成局域共振型聲子晶體雙層板結(jié)構(gòu),并對(duì)其展現(xiàn)出的帶隙特性進(jìn)行詳細(xì)研究,為處理工程上低頻域減振降噪設(shè)計(jì)具有理論意義和應(yīng)用價(jià)值。

    近年來,通過在一些如桿、梁、單板等基本彈性結(jié)構(gòu)中引入局域共振型聲子晶體設(shè)計(jì)思路所構(gòu)成的周期性結(jié)構(gòu)獲得了廣泛的研究。對(duì)于桿梁結(jié)構(gòu),Wang等[9]研究了縱向彈性波在由周期貼附諧振子的細(xì)長梁所構(gòu)成的準(zhǔn)一維結(jié)構(gòu)中的傳播特性。Yu等[10]通過實(shí)驗(yàn)和理論相結(jié)合研究了在周期布置諧振子的鐵摩辛柯梁中的彎曲振動(dòng)特性,理論和實(shí)驗(yàn)結(jié)果具備很好的一致性。對(duì)于局域共振型聲子晶體單板結(jié)構(gòu),在研究過程中形成了兩種結(jié)構(gòu)設(shè)計(jì)思路:a.填充型系統(tǒng),主要是通過將單板周期挖孔并填充軟材料而形成;b.貼附型系統(tǒng),通過在單板表面周期貼附共振單元而形成。文獻(xiàn)[11-12]分別以橡膠和包覆橡膠的硬材料作為共振單元研究了二元和三元填充型局域共振聲子晶體單板結(jié)構(gòu)的帶隙特性,結(jié)果表明,通過一定的參數(shù)調(diào)節(jié),在低頻域均可獲得一條完整的帶隙,并且通過對(duì)半徑、厚度等相關(guān)參數(shù)的調(diào)節(jié)可以實(shí)現(xiàn)帶隙的調(diào)控。與填充型系統(tǒng)類似,Oudich等[13]分別研究了在單板表面周期貼附柱狀橡膠共振單元和橡膠-鉛共振單元所構(gòu)成的二元和三元貼附型系統(tǒng)的帶隙特性,并指出低頻帶隙是所貼附共振單元的局域共振與基板Lamb 波模態(tài)相互耦合的結(jié)果。此外,以彈簧振子作為貼附型系統(tǒng)共振單元的簡化模型, Xiao等[14]研究了彎曲波在周期布置彈簧-質(zhì)量塊的單板結(jié)構(gòu)上的傳播特性,并且發(fā)現(xiàn)彈簧-質(zhì)量塊的共振頻率與帶隙存在很大的對(duì)應(yīng)關(guān)系。上述所有關(guān)于貼附型系統(tǒng)的研究均展示了能帶結(jié)構(gòu)可以在低頻域打開一條帶隙,但是由于對(duì)應(yīng)于沿著板面方向振動(dòng)模態(tài)和垂直于板面方向振動(dòng)模態(tài)的能帶錯(cuò)綜交織在一起,因而帶隙寬度很窄。基于此, Ma等[15]提出了一個(gè)由3層球狀共振單元構(gòu)成的新的局域共振型聲子晶體單板結(jié)構(gòu),發(fā)現(xiàn)由該種結(jié)構(gòu)形成的能帶結(jié)構(gòu)圖中出現(xiàn)了一個(gè)較寬的低頻段帶隙。在填充型系統(tǒng)和貼附型系統(tǒng)的研究基礎(chǔ)上, Li等[16]將這兩種系統(tǒng)的共振單元結(jié)合起來,研究了Lamb波在由該種綜合型共振單元構(gòu)成的聲子晶體單板結(jié)構(gòu)中的傳播特性,發(fā)現(xiàn)該種聲子晶體板的帶隙較原有兩類聲子晶體板的帶隙均有了很大改善。Qian等[17]通過在彈簧振子周圍環(huán)繞附加彈簧,研究了由該類共振單元構(gòu)成的局域共振型聲子晶體雙層板結(jié)構(gòu)的彎曲振動(dòng)傳輸特性。

    筆者在現(xiàn)有對(duì)形式各異的局域共振型聲子晶體單板結(jié)構(gòu)帶隙特性研究的基礎(chǔ)上,通過在雙層板間周期貼附由軟-硬-軟材料構(gòu)成的共振單元,從而構(gòu)造貼附型局域共振聲子晶體雙層板結(jié)構(gòu)。采用有限元法研究其帶隙的形成機(jī)制和調(diào)節(jié)規(guī)律,并通過與相應(yīng)單板結(jié)構(gòu)的帶隙特性進(jìn)行對(duì)比,研究了單/雙板的通性以及雙板的特性。此外,還研究了空腔與雙層板間聲振耦合效應(yīng)以及軟材料黏性對(duì)帶隙的影響規(guī)律和機(jī)理。

    1 模型和方法

    如圖1(a)所示,貼附型局域共振聲子晶體雙層板結(jié)構(gòu)是通過在雙層板結(jié)構(gòu)的上下板面之間周期貼附2組分圓柱狀局域共振柱體而形成。對(duì)于該柱體,其由上中下3層粘貼而成,其中上層和下層完全一致,由軟材料橡膠構(gòu)成,而中間層由硬材料鉛構(gòu)成。此外,雙層板的上下板完全一致,且由鋁材料構(gòu)成。在晶胞中,晶格常數(shù)、上下板的厚度、圓柱狀柱體的半徑、橡膠層柱體的高度以及鉛層柱體的高度分別為a,e,r,h1和h2,如圖1(b)所示。表1給出了計(jì)算所用到的所有材料參數(shù)。筆者在研究材料為線彈性和各向同性時(shí)該結(jié)構(gòu)的帶隙特性基礎(chǔ)上,繼而研究了軟材料黏性對(duì)帶隙的影響規(guī)律。

    對(duì)于貼附型局域共振聲子晶體雙層板的能帶結(jié)構(gòu),采用有限元法來計(jì)算,并且借助于商業(yè)軟件COMSOL Multiphysics來實(shí)現(xiàn)。對(duì)于有限元網(wǎng)格的劃分,均采用四面體網(wǎng)格且晶胞的有限元模型如圖1(c)所示。從圖中可以看出,在計(jì)算能帶結(jié)構(gòu)時(shí)只需考慮單個(gè)晶胞,這可以歸因于結(jié)構(gòu)的周期性。對(duì)于邊界條件,非接觸面均采用無壓力自由邊界條件,而相鄰晶胞的交界面采用由Bloch-Floquet定理導(dǎo)出的周期邊界條件[12-13],即

    ui(x+a,y+a)=ui(x,y)e-i(kxa+kya)

    (i=x,y,z)

    (1)

    其中:當(dāng)i取x,y和z時(shí),ui分別表示沿x方向、y方向和z方向的位移;kx和ky為不可約Brillouin區(qū)邊界上的二維Bloch波矢k沿x方向和y方向的分量(見圖1(d))。

    對(duì)于有限元網(wǎng)格劃分,網(wǎng)格尺寸越小,計(jì)算收斂性越好,但相應(yīng)計(jì)算時(shí)間越長。所以,選擇合適的網(wǎng)格劃分很重要。模型中由于橡膠比環(huán)氧樹脂和鉛軟,變形更大,所以橡膠層的網(wǎng)格劃分必須要更細(xì)化。

    將周期邊界條件式(1)代入到自由振動(dòng)的有限元特征方程中,可以得到

    (K-ω2M)u=0

    (2)

    在式(2)中,剛度矩陣K和質(zhì)量矩陣M均含有與Bloch波矢相耦合的項(xiàng),與傳統(tǒng)有限元特征方程中的剛度矩陣和質(zhì)量矩陣有所區(qū)別。

    式(2)即為關(guān)于ω2的典型廣義特征值問題。對(duì)于每一給定的Bloch波矢k,通過求解特征值可得到相應(yīng)的一系列特征頻率。通過遍歷所有的不可約Brillouin區(qū)邊界上的波矢,最終可得到該貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的能帶結(jié)構(gòu)。

    2 數(shù)值結(jié)果和分析

    2.1 能帶結(jié)構(gòu)、位移場及傳輸曲線

    圖2(b)給出了貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的能帶結(jié)構(gòu),其中計(jì)算所用到的材料參數(shù)和幾何參數(shù)分別如表1和表2所示。為了驗(yàn)證計(jì)算所得能帶結(jié)構(gòu)的準(zhǔn)確性,相應(yīng)有限雙層板結(jié)構(gòu)的彎曲振動(dòng)和縱向振動(dòng)的傳輸曲線分別如圖2(a)和(c)所示。該有限雙層板結(jié)構(gòu)由8×8個(gè)晶胞周期排列而成,并且激勵(lì)點(diǎn)和響應(yīng)點(diǎn)分別放置在下板的一端和上板的另一端,如圖3所示。

    表1 計(jì)算所用到的材料參數(shù)

    表2 計(jì)算所用到的幾何參數(shù)

    圖2 貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的能帶結(jié)構(gòu)和相應(yīng)有限8×8雙層板結(jié)構(gòu)的傳輸曲線Fig.2 Band structure of the proposed double panel structure and transmission power spectrums of the corresponding finite 8×8 system

    圖3 由8×8個(gè)晶胞周期排列而成的有限局域共振雙層板結(jié)構(gòu)的有限元模型Fig.3 Meshing of the finite locally resonant double panel structure made of 8×8 unit cells

    由圖2(b)可以看到,能帶結(jié)構(gòu)在72~83 Hz之間打開了一條狹窄的帶隙;但是由圖2(a)和(c)可以看到,不論是彎曲振動(dòng)還是縱向振動(dòng)在有限結(jié)構(gòu)中傳遞,傳輸曲線的衰減頻域均非常寬,這與已有研究中所得的無限結(jié)構(gòu)的帶隙頻段和有限結(jié)構(gòu)的衰減頻段相吻合這一結(jié)論[17-19]相違背。為了發(fā)現(xiàn)更多雙層板結(jié)構(gòu)傳輸曲線的衰減特性,當(dāng)激勵(lì)點(diǎn)和響應(yīng)點(diǎn)均放置在下板時(shí),縱向振動(dòng)傳輸曲線如圖2(d)所示。由圖可以看出,激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在不同板時(shí)的縱向振動(dòng)傳輸曲線所存在的很寬的衰減頻域并未出現(xiàn)在這里,并且與之對(duì)應(yīng)的頻域內(nèi)的振動(dòng)還有所加強(qiáng)。為了揭示該有限局域共振雙層板結(jié)構(gòu)的傳輸曲線展示出來的典型現(xiàn)象的機(jī)理,作出了與圖2(b)中標(biāo)注的特征模態(tài)相對(duì)應(yīng)的位移場,如圖4所示。

    圖4 與圖2(b)中標(biāo)注的特征模態(tài)相對(duì)應(yīng)的位移場Fig.4 The displacement fields of eigenmodes labeled in Fig. 2(b)

    對(duì)于模態(tài)B3,B8和B9,位移場為柱體共振單元沿z方向拉伸振動(dòng)的模態(tài)與雙層板垂直板面向外彎曲振動(dòng)的模態(tài)共同作用的結(jié)果。其中:模態(tài)B3中振動(dòng)能量集中在柱體中而雙層板保持靜止;模態(tài)B8和B9則是柱體中的鉛層充當(dāng)靜止層,但不同的是,模態(tài)B8中上下板以相反的相位達(dá)到彎曲振動(dòng)動(dòng)態(tài)平衡,而模態(tài)B9中上下板以相同的相位達(dá)到彎曲振動(dòng)動(dòng)態(tài)平衡?;诖?,模態(tài)B8和B9被分別稱之為對(duì)稱彎曲振動(dòng)模態(tài)和反對(duì)稱彎曲振動(dòng)模態(tài)。這些模態(tài)耦合作用的結(jié)果,一條在模態(tài)B3,B8和B9之間的頻率區(qū)間為72~230 Hz的局部彎曲振動(dòng)帶隙被打開,該帶隙頻段與圖2(a)所示的彎曲振動(dòng)傳輸曲線中的衰減頻段基本吻合。

    對(duì)于模態(tài)B1和B2,位移場為柱體共振單元繞xy平面內(nèi)軸線旋轉(zhuǎn)振動(dòng)的模態(tài)與雙層板垂直板面向外剪切振動(dòng)的模態(tài)相互耦合的結(jié)果。其中,模態(tài)B1和B2中振動(dòng)能量集中在柱體中而雙層板保持靜止。對(duì)于模態(tài)B4~B7,位移場為柱體共振單元沿xy平面內(nèi)軸線平移振動(dòng)的模態(tài)與雙層板在板面內(nèi)剪切振動(dòng)的模態(tài)共同作用的結(jié)果。在這4種振動(dòng)模態(tài)中,柱體的中間鉛層均保持靜止。不同的是,模態(tài)B4和B5中上下板以相反的相位達(dá)到縱向振動(dòng)動(dòng)態(tài)平衡,而模態(tài)B6和B7中上下板以相同的相位達(dá)到縱向振動(dòng)動(dòng)態(tài)平衡。相似地,模態(tài)B4~B5和模態(tài)B6~B7被分別稱之為反對(duì)稱縱向振動(dòng)模態(tài)和對(duì)稱縱向振動(dòng)模態(tài)。作為這些模態(tài)耦合作用的結(jié)果,一條頻率區(qū)間為66~83Hz的局部縱向振動(dòng)帶隙被打開,但在圖2(c)和(d)所示的相應(yīng)頻段的縱向振動(dòng)傳輸曲線中的衰減并不明顯,這與文獻(xiàn)[20]所描述的現(xiàn)象一致。

    此外,由圖4可以看出,反對(duì)稱縱向振動(dòng)模態(tài)B4中上板的振動(dòng)相位與對(duì)稱縱向振動(dòng)模態(tài)B6中上板的振動(dòng)相位相反。同樣地,反對(duì)稱縱向振動(dòng)模態(tài)B5中上板的振動(dòng)相位與對(duì)稱縱向振動(dòng)模態(tài)B7中上板的振動(dòng)相位相反。于是,當(dāng)振動(dòng)模態(tài)B4~B7相疊加時(shí),上層板的振動(dòng)便被削弱,這就是圖2(c)所示的激勵(lì)點(diǎn)和響應(yīng)點(diǎn)分別放置在上下板時(shí)縱向振動(dòng)傳輸曲線中存在很寬頻段衰減域的原因。與此同時(shí),在疊加過程中,下層板的振動(dòng)會(huì)加強(qiáng),這同樣是因?yàn)閳D2(d)所示的激勵(lì)點(diǎn)和響應(yīng)點(diǎn)均放置在下板時(shí)縱向振動(dòng)傳輸曲線中存在很寬頻段加強(qiáng)域。為了進(jìn)一步說明該衰減特性,給出了頻率落在衰減域的有限雙層板結(jié)構(gòu)的振型圖,如圖5所示。這里,計(jì)算所用到的模型如圖3所示,激勵(lì)點(diǎn)施加三個(gè)方向的位移激勵(lì)并且頻率點(diǎn)選為f=140 Hz。由圖5可以看出,當(dāng)彎曲激勵(lì)和縱向激勵(lì)均施加于下板時(shí),彎曲振動(dòng)和縱向振動(dòng)均無法通過上層板傳播,而僅僅縱向振動(dòng)可以通過下層板傳播??偟膩碚f,當(dāng)振源和響應(yīng)域分別在雙層板結(jié)構(gòu)的異側(cè)時(shí),以模態(tài)B2和B3對(duì)應(yīng)的更高頻率為起始頻率以及模態(tài)B8對(duì)應(yīng)的頻率為終止頻率的帶隙可以被認(rèn)為是一條在相應(yīng)有限結(jié)構(gòu)中彎曲和縱向振動(dòng)均有衰減的完全帶隙。

    圖5 頻率落在衰減域的有限雙層板結(jié)構(gòu)的振型圖(f=140 Hz)Fig.5 The vibration mode of the frequency located inside the frequency range of attenuation(f=140 Hz)

    2.2 聲振耦合對(duì)能帶結(jié)構(gòu)的影響

    圖6 晶胞內(nèi)結(jié)構(gòu)有限元模型及空腔有限元模型Fig.6 Meshing of the structure and cavity in a unit cell

    如果不忽略雙板間的空氣,則可進(jìn)一步考慮聲振耦合效應(yīng)對(duì)貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的影響。這里分別對(duì)晶胞內(nèi)結(jié)構(gòu)部分和空腔部分進(jìn)行有限元建模,如圖6所示,并分別在結(jié)構(gòu)和空腔部分的邊界上布置周期邊界條件。計(jì)算用到的結(jié)構(gòu)部分材料參數(shù)和幾何參數(shù)均與表1和表2所示相同。此外,空腔部分空氣的參數(shù)為:密度ρ0=1.29 kg/m3;聲速c0=340 m/s。

    與圖2類似,圖7給出了考慮聲振耦合效應(yīng)時(shí)的貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的能帶結(jié)構(gòu)、激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在不同板時(shí)的彎曲振動(dòng)傳輸曲線、激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在不同板時(shí)的縱向振動(dòng)傳輸曲線,以及激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在相同板時(shí)的縱向振動(dòng)傳輸曲線。通過對(duì)比7(b)和2(b)可以發(fā)現(xiàn),聲振耦合對(duì)帶隙形成的主要能帶B1~B7和反對(duì)稱彎曲振動(dòng)能帶B9基本無影響,但卻對(duì)對(duì)稱彎曲振動(dòng)能帶B8影響較大,使得原有帶隙將不再被打開。但是對(duì)比圖7(a),(c),(d)和圖2(a),(c),(d)可以發(fā)現(xiàn),雖然聲振耦合對(duì)傳輸曲線有影響,但是并不影響激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在不同板時(shí)的彎曲振動(dòng)和縱向振動(dòng)的傳輸衰減頻域,且同樣激勵(lì)點(diǎn)和響應(yīng)點(diǎn)放置在相同板時(shí)的縱向振動(dòng)傳輸曲線不存在衰減區(qū)域。由此可見,傳輸曲線展現(xiàn)出的現(xiàn)象與能帶結(jié)構(gòu)展現(xiàn)的現(xiàn)象不一致。為了進(jìn)一步揭示該特性,筆者作出了頻率落在衰減域的有限雙層板結(jié)構(gòu)的振型和空腔內(nèi)的聲壓場。

    圖7 考慮聲振耦合效應(yīng)的貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的能帶結(jié)構(gòu)和相應(yīng)有限8×8雙層板結(jié)構(gòu)的傳輸曲線Fig.7 Band structure of the proposed double panel structure and transmission power spectrums of the corresponding finite 8×8 system with the coupling between sound and vibration

    圖8 頻率落在衰減域的有限雙層板結(jié)構(gòu)的振型圖(f=160 Hz)Fig.8 The vibration mode of the frequency located inside the frequency range of attenuation(f=160 Hz)

    圖9 頻率落在衰減域的有限雙層板結(jié)構(gòu)的聲壓場(f=160 Hz)Fig.9 The sound pressure field of the frequency located inside the frequency range of attenuation(f=160 Hz)

    圖8和圖9分別給出了頻率落在衰減域的有限雙層板結(jié)構(gòu)的振型和空腔內(nèi)的聲壓場。這里,激勵(lì)點(diǎn)施加三個(gè)方向的位移激勵(lì)并且頻率點(diǎn)選為f=160 Hz。由圖可以看出,在衰減頻域內(nèi),當(dāng)彎曲激勵(lì)和縱向激勵(lì)均施加于下板時(shí),彎曲振動(dòng)和縱向振動(dòng)均無法通過上層板傳播,而僅僅縱向振動(dòng)可以通過下層板傳播;但是聲壓則可以沿著空腔無限傳播。這就是圖7(b)所示能帶結(jié)構(gòu)中帶隙不能被打開的原因,即能帶B8表征了雙板對(duì)稱彎曲振動(dòng)模態(tài)與空腔內(nèi)聲場相耦合,雖然結(jié)構(gòu)部分振動(dòng)被抑制,但聲場部分卻一直在傳播。對(duì)于外場來講,結(jié)構(gòu)的振動(dòng)和經(jīng)由結(jié)構(gòu)振動(dòng)帶來的輻射噪聲可以認(rèn)為被很好地抑制。因此,基于減振降噪的應(yīng)用背景,在研究包含空氣聲腔的雙層板能帶結(jié)構(gòu)的帶隙特性時(shí),聲振耦合效應(yīng)可以完全被忽略。

    2.3 同參數(shù)下雙板結(jié)構(gòu)與單板結(jié)構(gòu)的能帶結(jié)構(gòu)對(duì)比

    若圖1(a)所示的貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的上層板和上層柱體橡膠層被取走,則可形成貼附型局域共振聲子晶體單板結(jié)構(gòu)。該單板結(jié)構(gòu)的能帶結(jié)構(gòu)如圖10(a)所示。其中,計(jì)算用到的參數(shù)與圖2(b)所示的參數(shù)一致,具體如表1和2所示。為了便于對(duì)比,再次給出了圖2(b)所示同參數(shù)下的雙層板結(jié)構(gòu)的能帶結(jié)構(gòu),如圖10(b)所示。此外,給出了與圖10(a)中標(biāo)注的特征模態(tài)相對(duì)應(yīng)的位移場,如圖11所示。

    圖10 同參數(shù)下局域共振型聲子晶體Fig.10 Band structures of the locally resonant single plate and the locally resonant double panel with the same parameters

    圖11 與圖10(a)中標(biāo)注的特征模態(tài)相對(duì)應(yīng)的位移場Fig.11 The displacement fields of eigenmodes labeled in Fig. 10(a)

    與雙層板結(jié)構(gòu)的能帶結(jié)構(gòu)相似,位移場為柱體共振單元沿z方向拉伸振動(dòng)的模態(tài)與單板垂直板面向外彎曲振動(dòng)的模態(tài)耦合作用的振動(dòng)模態(tài)b3和b6在51~235 Hz之間打開了一條局部彎曲振動(dòng)帶隙。此外,位移場為柱體共振單元繞xy平面內(nèi)軸線旋轉(zhuǎn)振動(dòng)的模態(tài)與單板垂直板面向外剪切振動(dòng)的模態(tài)耦合作用的振動(dòng)模態(tài)b1和b2以及位移場為柱體共振單元沿xy平面內(nèi)軸線平移振動(dòng)的模態(tài)與雙層板在板面內(nèi)剪切振動(dòng)的模態(tài)耦合作用的振動(dòng)模態(tài)b4和b5的共同作用下,在47.2~-83 Hz之間打開了一條局部縱向振動(dòng)帶隙。因此,在彎曲振動(dòng)和縱向振動(dòng)的共同作用下一條較為狹窄的帶隙在51~83Hz之間被打開。

    圖12 用于幫助理解振動(dòng)模態(tài)的“基底-彈簧-質(zhì)量塊”簡化模型Fig.12 The simplified base-spring-mass models applied to help understand the vibration modes

    總的來說,同參數(shù)下局域共振型聲子晶體雙層板結(jié)構(gòu)和單板結(jié)構(gòu)的能帶結(jié)構(gòu)之間有很多相似之處,且相應(yīng)能帶的整體趨勢(shì)大致相同。此外,一些雙層板結(jié)構(gòu)的能帶可以被認(rèn)為是單板結(jié)構(gòu)的能帶分開而形成的。但是從減振降噪的角度來說,當(dāng)激勵(lì)和響應(yīng)域在結(jié)構(gòu)的不同側(cè)時(shí),雙層板結(jié)構(gòu)所擁有的很寬的完整帶隙這一性質(zhì)是單板結(jié)構(gòu)所遠(yuǎn)遠(yuǎn)不能比擬的。

    2.4 各參數(shù)對(duì)雙層板結(jié)構(gòu)帶隙的影響規(guī)律

    振源和接收者常常處于雙層板結(jié)構(gòu)的不同板的一側(cè)。因此,在下面的研究中由能帶B2,B3和B8所打開的帶隙被用來作為研究對(duì)象。這里柱體橡膠層的高度h1、柱體鉛層的高度h2、柱體半徑r、上下板的厚度e以及晶格常數(shù)a被選取為影響參數(shù),用來研究幾何參數(shù)對(duì)帶隙的起始頻率f2(或f3)、終止頻率f8和帶隙寬度fw的影響規(guī)律。在研究過程中,除了影響參數(shù)外別的參數(shù)均與圖2(a)所示算例的參數(shù)一致,具體如表1和表2所示。

    圖13 帶隙臨界頻率f2,f3,f8和帶隙寬度fw與柱體橡膠層高度h1之間的關(guān)系Fig.13 The influences of rubber layer height in the pillar h1 on f2,f3,f8 and fw

    圖14給出了柱體鉛層高度h2對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系,fw同樣由起始頻率f3和終止頻率f8決定。由圖可以看出,f2和f3均隨著h2的增大而減小,這是因?yàn)橹w鉛層的等效質(zhì)量m隨著h2的增大而增大。而f8基本不受h2的影響,這是因?yàn)樵趫D4所對(duì)應(yīng)的B8振動(dòng)模態(tài)中鉛層是靜止不動(dòng)的。于是,h2對(duì)f3和f8的影響關(guān)系使得fw隨h2的增大而增大,如圖14所示。

    圖14 帶隙臨界頻率f2,f3,f8和帶隙寬度fw與柱體鉛層高度h2之間的關(guān)系Fig.14 The influences of Pb layer height in the pillar h2 on f2, f3, f8 and fw

    圖15 帶隙臨界頻率f2,f3,f8和帶隙寬度fw與柱體半徑r之間的關(guān)系Fig.15 The influences of pillar radius r on f2,f3,f8 and fw

    柱體半徑r對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系如圖15所示。由圖可以看出,隨著r的增大,f2和f3緩慢地增大,這是因?yàn)榈刃偠萲和等效質(zhì)量塊質(zhì)量m均隨著r的增大而增大,并且k的增大幅度略大于m的增大幅度。此外,隨著r的增大,f8急速地增大,這是因?yàn)槌洚?dāng)質(zhì)量塊的鉛層在這里保持靜止,僅僅k的增大使得f8的增長幅度比f2和f3大得多。于是,終止頻率f8的急速增大和起始頻率f3的緩慢增大使得fw隨著r的增大而增大。

    圖16 帶隙臨界頻率f2,f3,f8和帶隙寬度fw與上下板厚度e之間的關(guān)系Fig.16 The influences of base plate thickness e on f2,f3,f8 and fw

    圖17 當(dāng)e分別取3和7 mm時(shí)振動(dòng)模態(tài)B2 和B3所對(duì)應(yīng)的位移場Fig.17 The displacement fields of modes B2 and B3 when e takes the value of 3 and 7mm respectively

    圖16給出了上下板的厚度e對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系。由圖可以看出,e=3.42 mm是臨界值。當(dāng)板的厚度小于該值時(shí),fw由起始頻率f2和終止頻率f8決定;而當(dāng)板的厚度大于該值時(shí),fw由起始頻率f3和終止頻率f8決定。此外,f2和f3均隨著e的增大而增大,這可以借助于圖17來幫助理解。該圖分別展示了當(dāng)e取3和7 mm時(shí)振動(dòng)模態(tài)B2和B3所對(duì)應(yīng)的位移場。由圖可以看出,振動(dòng)模態(tài)B2和B3中上下板并非絕對(duì)靜止,并且當(dāng)e取3 mm時(shí)上下板的振動(dòng)更容易展現(xiàn)出來,這是因?yàn)閷?shí)際振動(dòng)中板存在一定的剛度,并且該剛度隨著板的厚度的減小而減小。因此,在這里等效彈簧剛度k可以進(jìn)一步認(rèn)為是板與柱體橡膠層共同作用的結(jié)果。由于隨著e的增大,板的彎曲剛度增大,繼而k也增大,這就導(dǎo)致f2和f3也隨之增大。此外,由于隨著e的增大,板的等效質(zhì)量M也增大,因此f8隨之減小。于是,終止頻率f8的減小和起始頻率f2(或f3)的增大使得fw隨著e的增大而減小。

    圖18 帶隙臨界頻率f2,f3,f8和帶隙寬度fw與晶格常量a之間的關(guān)系Fig.18 The influences of lattice constant a on f2,f3,f8 and fw

    圖19 當(dāng)a分別取85和200 mm時(shí)振動(dòng)模態(tài)B2和B3所對(duì)應(yīng)的位移場Fig.19 The displacement fields of modes B2 and B3 when a takes the value of 85and 200 mm respectively

    圖18給出了晶格常量a對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系。從圖中可以看出,a=0.132m是臨界值。當(dāng)晶格常量小于該值時(shí),fw由起始頻率f3和終止頻率f8決定;而當(dāng)晶格常量大于該值時(shí),fw由起始頻率f2和終止頻率f8決定。此外,隨著a的增大,f2基本保持不變,而f3減小,這同樣可以借助于圖19所示的當(dāng)a分別取0.085和0.2 m時(shí),振動(dòng)模態(tài)B2和B3所對(duì)應(yīng)的位移場來幫助理解。由圖可以看出,板的等效剛度隨著a的增大而減小,但是振動(dòng)模態(tài)B2中板的等效剛度的變化不是很明顯,而振動(dòng)模態(tài)B3中板的等效剛度的變化則很大,這就是f2和f3隨著a的增大分別保持不變和減小的原因。另外,由于隨著a的增大,板的等效質(zhì)量M也增大,所以f8隨之減小。于是,終止頻率f8的快速減小和起始頻率f2的基本不變(或f3的緩慢減小)使得fw隨著e的增大而減小。

    2.5 柱體橡膠層黏性對(duì)雙層板結(jié)構(gòu)帶隙的影響規(guī)律

    當(dāng)考慮復(fù)合板軟材料的黏性作用時(shí),在頻域內(nèi),軟材料的各模量是與頻率相關(guān)的,方程(2)的特征值問題需要通過迭代過程來求解,具體計(jì)算流程可參照文獻(xiàn)[21-23]。

    圖20(a),(b)和(c)顯示了橡膠層松弛時(shí)間為τ=3×10-5s,初始-終止值比α分別取1.5,2.0和3.0時(shí)的能帶結(jié)構(gòu)圖。其中,計(jì)算用到的參數(shù)如表1和表2所示。由圖可以看出,α對(duì)所有能帶均有顯著影響。為了進(jìn)一步說明該影響規(guī)律,作出了帶隙臨界頻率f2,f3,f8和帶隙寬度fw隨α的變化曲線圖,如圖20(d)所示。由圖可以看出,隨著α的增大,頻率f2,f3和f8均減小。此外,帶隙終止頻率f8的減小速率大于起始頻率f2的減小速率,這是因?yàn)樵谘芯康念l域內(nèi)(300 Hz以下),軟材料的材料參數(shù)的變化率隨著頻率的增大而增大[23],導(dǎo)致帶隙寬度fw也隨著α的增大而減小。由此可見,通過增大黏彈性軟材料的初始-終止值比對(duì)降低貼附型局域共振聲子晶體雙層板結(jié)構(gòu)的帶隙起始頻率起到積極作用,但同時(shí)也會(huì)縮減其帶隙寬度。

    圖20 黏彈性材料松弛時(shí)間為τ=3×10-5 s時(shí)不同初始-終止值比下的能帶結(jié)構(gòu)對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系Fig.20 Band structures with the viscidity of soft material considered when relax time τ=3×10-5 s, the influences of initial-final value ratio α on critical frequencies f2,f3,f8 and bandgap width fw

    圖21 黏彈性材料初始-終止值比為α=2.0時(shí)不同松弛時(shí)間下的能帶結(jié)構(gòu)對(duì)帶隙臨界頻率f2,f3,f8和帶隙寬度fw的影響關(guān)系Fig.21 Band structures with the viscidity of soft material considered when initial-final value ratio α=2.0, the influences of relax time τ on critical frequencies f2,f3,f8 and bandgap width fw

    圖21(a),(b)和(c)顯示了軟材料初始-終止值比為α=2.0,松弛時(shí)間τ分別取3×10-5,1×10-4,5×10-4s時(shí)的能帶結(jié)構(gòu)圖。其中,計(jì)算用到的參數(shù)同樣如表1和表2所示。由圖可以看出,τ對(duì)高頻段能帶影響顯著,但對(duì)低頻段能帶幾乎沒影響。為了進(jìn)一步說明該影響規(guī)律,作出了帶隙臨界頻率f2,f3,f8和帶隙寬度fw隨τ的變化曲線圖,如圖21(d)所示。由圖可以看出,隨著τ的增大,頻率f8增大,而f2和f3幾乎不變,這可以歸因于在研究的頻域內(nèi)(300 Hz以下),軟材料的材料參數(shù)的變化率隨著頻率的增大而增大,并且在低頻段不同的松弛時(shí)間所對(duì)應(yīng)的軟材料參數(shù)均趨于同一定值,為黏彈性軟材料參數(shù)的終止值[23]。由于終止頻率f8的增大以及起始頻率f2幾乎不變,所以帶隙寬度fw也隨著τ的增大而增大。由此可見,通過增大黏彈性軟材料的松弛時(shí)間對(duì)增大貼附型局域共振聲子晶體雙層板結(jié)構(gòu)帶隙寬度起到積極作用,且不會(huì)提高帶隙起始頻率。

    3 結(jié) 論

    1) 當(dāng)激勵(lì)點(diǎn)和響應(yīng)點(diǎn)位于上下異側(cè)板時(shí),通過對(duì)能帶結(jié)構(gòu)和傳輸曲線進(jìn)行分析,一條起始頻率低并且?guī)秾挼耐耆珟犊梢员淮蜷_。通過對(duì)帶隙邊界能帶的特征模態(tài)相對(duì)應(yīng)的位移場進(jìn)行分析,不論是彎曲振動(dòng)還是縱向振動(dòng),都展現(xiàn)出對(duì)稱與反對(duì)稱振動(dòng)模式。帶隙的打開可以認(rèn)為是共振單元的振動(dòng)模態(tài)和上下板的板振動(dòng)模態(tài)相互耦合作用的結(jié)果。通過對(duì)比分析,雙板結(jié)構(gòu)的能帶可以被考慮為在上下板作用下單板相對(duì)應(yīng)能帶分開而導(dǎo)致的。研究還表明,由空氣聲腔所帶來的聲振耦合效應(yīng)可以不考慮。

    2) 分析相應(yīng)幾何參數(shù)對(duì)能帶結(jié)構(gòu)的影響,可以歸納為:a.通過增加鉛層柱體的高度或者減小上下基板的厚度,可以降低帶隙起始頻率且增大帶隙寬度;b.通過增加橡膠層柱體高度,帶隙起始頻率和帶隙寬度同時(shí)獲得增長;c.通過增加柱體半徑或者減小晶格常數(shù),可以增大帶隙寬度且對(duì)帶隙起始頻率影響不大。所有上述影響規(guī)律均可借助于“基底-彈簧-質(zhì)量塊”簡化模型獲得相應(yīng)解釋。此外,當(dāng)考慮軟材料黏性時(shí),通過增大黏彈性軟材料的初始-終止值比或者減小黏彈性軟材料的松弛時(shí)間可以降低帶隙起始頻率,但同時(shí)也會(huì)縮減其帶隙寬度。

    猜你喜歡
    帶隙板結(jié)構(gòu)聲子
    半無限板類聲子晶體帶隙仿真的PWE/NS-FEM方法
    密度泛函理論計(jì)算半導(dǎo)體材料的帶隙誤差研究
    環(huán)氧樹脂基復(fù)合材料加筋板結(jié)構(gòu)吸濕行為研究
    納米表面聲子 首次實(shí)現(xiàn)三維成像
    聲子晶體覆蓋層吸聲機(jī)理研究
    一種基于BJT工藝的無運(yùn)放低溫度系數(shù)的帶隙基準(zhǔn)源
    間距比對(duì)雙振子局域共振軸縱振帶隙的影響
    一款高PSRR低溫度系數(shù)的帶隙基準(zhǔn)電壓源的設(shè)計(jì)
    電子制作(2018年1期)2018-04-04 01:48:38
    基于聲子晶體理論的導(dǎo)線防舞方法及數(shù)值驗(yàn)證
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    夜夜爽天天搞| 老司机深夜福利视频在线观看| 人人妻人人澡欧美一区二区| 天堂动漫精品| 国产精品1区2区在线观看.| 免费看日本二区| 国产乱人伦免费视频| 两人在一起打扑克的视频| 丰满人妻一区二区三区视频av | 97超级碰碰碰精品色视频在线观看| 久久精品影院6| 丁香欧美五月| 在线观看www视频免费| 久久 成人 亚洲| 国产片内射在线| 日本在线视频免费播放| 女警被强在线播放| 精品国产乱码久久久久久男人| 99热只有精品国产| 日本五十路高清| 757午夜福利合集在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品无人区乱码1区二区| 亚洲黑人精品在线| 哪里可以看免费的av片| 国产精品免费视频内射| 一区二区三区国产精品乱码| xxx96com| 亚洲一区高清亚洲精品| tocl精华| 久久久久久久精品吃奶| 国产午夜精品久久久久久| 日韩 欧美 亚洲 中文字幕| 99久久久亚洲精品蜜臀av| 国产伦一二天堂av在线观看| 观看免费一级毛片| 日韩大码丰满熟妇| 久久中文字幕一级| 18禁国产床啪视频网站| 村上凉子中文字幕在线| 久久国产精品影院| 亚洲 欧美一区二区三区| 欧美中文综合在线视频| 午夜老司机福利片| 免费搜索国产男女视频| 一二三四社区在线视频社区8| 丁香欧美五月| 亚洲中文字幕一区二区三区有码在线看 | 看片在线看免费视频| 国产精品美女特级片免费视频播放器 | 可以免费在线观看a视频的电影网站| 日韩欧美三级三区| 在线观看一区二区三区| 日韩免费av在线播放| 久久久久免费精品人妻一区二区| 啦啦啦观看免费观看视频高清| 久久 成人 亚洲| 亚洲成av人片在线播放无| 啦啦啦韩国在线观看视频| 亚洲国产中文字幕在线视频| 国产精品免费视频内射| 成人一区二区视频在线观看| 蜜桃久久精品国产亚洲av| 欧美国产日韩亚洲一区| www.精华液| 在线播放国产精品三级| 蜜桃久久精品国产亚洲av| www.999成人在线观看| 18禁美女被吸乳视频| 色综合站精品国产| 亚洲精品国产精品久久久不卡| 午夜福利免费观看在线| 叶爱在线成人免费视频播放| 九色成人免费人妻av| 黄片小视频在线播放| 亚洲av美国av| 免费av毛片视频| 丰满的人妻完整版| 1024香蕉在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产男靠女视频免费网站| 亚洲国产精品sss在线观看| 久久亚洲精品不卡| 精品乱码久久久久久99久播| 亚洲精品色激情综合| 国产一区二区在线av高清观看| 九色成人免费人妻av| 九九热线精品视视频播放| 国产精品自产拍在线观看55亚洲| 欧美性猛交黑人性爽| 亚洲无线在线观看| 午夜福利免费观看在线| 国产精品九九99| 亚洲av成人av| 亚洲av电影在线进入| 中文字幕熟女人妻在线| 久久中文字幕一级| 日日干狠狠操夜夜爽| 香蕉丝袜av| 夜夜爽天天搞| 桃红色精品国产亚洲av| 国产av又大| av福利片在线| 亚洲人成伊人成综合网2020| 久久久久久亚洲精品国产蜜桃av| 欧美色欧美亚洲另类二区| 熟女少妇亚洲综合色aaa.| 黑人欧美特级aaaaaa片| 精品久久久久久成人av| 91麻豆av在线| 久久伊人香网站| 天堂av国产一区二区熟女人妻 | 又黄又粗又硬又大视频| 在线国产一区二区在线| 日韩精品中文字幕看吧| 变态另类成人亚洲欧美熟女| 久久精品91蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久成人aⅴ小说| 午夜福利免费观看在线| 麻豆av在线久日| 免费在线观看黄色视频的| 精品久久蜜臀av无| 一本精品99久久精品77| 啦啦啦观看免费观看视频高清| 亚洲av成人一区二区三| 亚洲av中文字字幕乱码综合| 国产精品乱码一区二三区的特点| 国产精品久久久久久久电影 | 老汉色av国产亚洲站长工具| 日本五十路高清| 色在线成人网| 亚洲av成人av| 日韩欧美免费精品| 男插女下体视频免费在线播放| 久久久久九九精品影院| 日韩有码中文字幕| 国产精品久久电影中文字幕| 亚洲国产精品久久男人天堂| 50天的宝宝边吃奶边哭怎么回事| 美女 人体艺术 gogo| 一个人观看的视频www高清免费观看 | 性欧美人与动物交配| 蜜桃久久精品国产亚洲av| 国产在线精品亚洲第一网站| 国产成人精品久久二区二区免费| 成人高潮视频无遮挡免费网站| 在线观看66精品国产| 国产精品综合久久久久久久免费| 免费看a级黄色片| 国产高清激情床上av| 亚洲美女视频黄频| 99精品欧美一区二区三区四区| 国产精品自产拍在线观看55亚洲| 欧美黑人巨大hd| 精品无人区乱码1区二区| 亚洲国产精品久久男人天堂| 精品久久久久久久人妻蜜臀av| 真人做人爱边吃奶动态| 曰老女人黄片| 欧美一区二区精品小视频在线| 精品人妻1区二区| 国产真人三级小视频在线观看| 两个人免费观看高清视频| 日韩三级视频一区二区三区| 国产一区二区三区视频了| 嫁个100分男人电影在线观看| 1024手机看黄色片| e午夜精品久久久久久久| 男女那种视频在线观看| 亚洲自偷自拍图片 自拍| 999精品在线视频| 婷婷六月久久综合丁香| 国产精品免费一区二区三区在线| 国产精品九九99| 最近最新免费中文字幕在线| 国产精品美女特级片免费视频播放器 | 在线永久观看黄色视频| 香蕉丝袜av| 最近最新免费中文字幕在线| 午夜精品在线福利| 熟女少妇亚洲综合色aaa.| 亚洲中文日韩欧美视频| 国产亚洲av高清不卡| 午夜福利18| 欧美av亚洲av综合av国产av| 国产三级黄色录像| 日韩三级视频一区二区三区| 18美女黄网站色大片免费观看| 欧美色欧美亚洲另类二区| 亚洲 欧美一区二区三区| 日韩国内少妇激情av| 在线观看www视频免费| 他把我摸到了高潮在线观看| 两个人视频免费观看高清| 国产成人影院久久av| 精品一区二区三区四区五区乱码| 制服人妻中文乱码| 欧美激情久久久久久爽电影| 欧美+亚洲+日韩+国产| 国产成年人精品一区二区| 狠狠狠狠99中文字幕| 亚洲成人久久爱视频| 欧美成人免费av一区二区三区| 国产99久久九九免费精品| 中文亚洲av片在线观看爽| 中文字幕人成人乱码亚洲影| 亚洲中文av在线| 久久香蕉激情| 精品欧美一区二区三区在线| 欧美最黄视频在线播放免费| 亚洲欧美精品综合一区二区三区| 国产男靠女视频免费网站| 99在线人妻在线中文字幕| 亚洲精华国产精华精| 两人在一起打扑克的视频| 国产99久久九九免费精品| 麻豆成人av在线观看| 欧美av亚洲av综合av国产av| 国产精品九九99| 麻豆国产av国片精品| 亚洲人成网站高清观看| 777久久人妻少妇嫩草av网站| 午夜影院日韩av| 久久久水蜜桃国产精品网| 午夜免费激情av| 亚洲精品国产一区二区精华液| 亚洲欧美日韩东京热| 亚洲国产精品sss在线观看| 狠狠狠狠99中文字幕| 国内精品久久久久久久电影| 免费在线观看日本一区| 激情在线观看视频在线高清| 人人妻人人澡欧美一区二区| 亚洲av美国av| 母亲3免费完整高清在线观看| 色综合欧美亚洲国产小说| 在线国产一区二区在线| 亚洲中文字幕日韩| 欧美av亚洲av综合av国产av| 99在线人妻在线中文字幕| 性色av乱码一区二区三区2| www.www免费av| 一本一本综合久久| 成人三级黄色视频| 国产精品日韩av在线免费观看| 1024视频免费在线观看| av片东京热男人的天堂| 啦啦啦观看免费观看视频高清| 麻豆成人av在线观看| ponron亚洲| 午夜福利欧美成人| 制服丝袜大香蕉在线| 欧美黑人巨大hd| 国产成人啪精品午夜网站| 欧美日韩瑟瑟在线播放| 国产野战对白在线观看| АⅤ资源中文在线天堂| 欧美在线一区亚洲| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情久久久久久爽电影| 丁香欧美五月| 免费看美女性在线毛片视频| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看| 男女那种视频在线观看| 精品少妇一区二区三区视频日本电影| 亚洲中文字幕一区二区三区有码在线看 | 亚洲性夜色夜夜综合| 亚洲一码二码三码区别大吗| 禁无遮挡网站| 欧美色视频一区免费| 人妻丰满熟妇av一区二区三区| 亚洲av片天天在线观看| 搞女人的毛片| 日本免费一区二区三区高清不卡| 中文字幕人妻丝袜一区二区| 在线观看美女被高潮喷水网站 | 色噜噜av男人的天堂激情| av国产免费在线观看| 国产私拍福利视频在线观看| 国产欧美日韩精品亚洲av| 国产熟女午夜一区二区三区| 美女 人体艺术 gogo| 制服诱惑二区| 精品一区二区三区av网在线观看| 欧美最黄视频在线播放免费| 搡老熟女国产l中国老女人| 日本黄色视频三级网站网址| 亚洲18禁久久av| 51午夜福利影视在线观看| 99精品久久久久人妻精品| 国产视频内射| 好看av亚洲va欧美ⅴa在| 精品不卡国产一区二区三区| 精品电影一区二区在线| 日韩欧美国产一区二区入口| 午夜a级毛片| cao死你这个sao货| 亚洲人成电影免费在线| 很黄的视频免费| 日本一本二区三区精品| 黄色a级毛片大全视频| 中文字幕av在线有码专区| 日韩中文字幕欧美一区二区| 亚洲一区二区三区色噜噜| 亚洲乱码一区二区免费版| 18禁美女被吸乳视频| 久99久视频精品免费| 伦理电影免费视频| 午夜久久久久精精品| 精品福利观看| 国内少妇人妻偷人精品xxx网站 | 最近最新中文字幕大全电影3| 麻豆一二三区av精品| 免费无遮挡裸体视频| 亚洲乱码一区二区免费版| 黄频高清免费视频| 麻豆国产av国片精品| 国产v大片淫在线免费观看| 久久国产精品人妻蜜桃| 在线观看66精品国产| 亚洲 国产 在线| 亚洲成人中文字幕在线播放| 国产精品九九99| 久久这里只有精品19| 亚洲欧洲精品一区二区精品久久久| 日本三级黄在线观看| 黄色丝袜av网址大全| 91九色精品人成在线观看| 精品国产超薄肉色丝袜足j| 精品久久久久久,| 欧美黄色片欧美黄色片| 一级作爱视频免费观看| 法律面前人人平等表现在哪些方面| 欧美一级a爱片免费观看看 | 亚洲va日本ⅴa欧美va伊人久久| 每晚都被弄得嗷嗷叫到高潮| 久久久久久九九精品二区国产 | 丁香欧美五月| 男女做爰动态图高潮gif福利片| 90打野战视频偷拍视频| 最近最新免费中文字幕在线| 精品一区二区三区av网在线观看| 在线十欧美十亚洲十日本专区| 日本免费一区二区三区高清不卡| 国产精品香港三级国产av潘金莲| 51午夜福利影视在线观看| 男女视频在线观看网站免费 | 日韩大尺度精品在线看网址| 国产精品综合久久久久久久免费| 18禁黄网站禁片免费观看直播| 国产一区二区在线av高清观看| 欧美日韩乱码在线| 国产真实乱freesex| 欧美中文日本在线观看视频| 亚洲天堂国产精品一区在线| 国产伦一二天堂av在线观看| 免费无遮挡裸体视频| 香蕉av资源在线| 亚洲aⅴ乱码一区二区在线播放 | 精品午夜福利视频在线观看一区| 亚洲欧美日韩东京热| 18禁观看日本| 久久精品aⅴ一区二区三区四区| 日本免费a在线| 久久精品91无色码中文字幕| 久久人妻av系列| 亚洲国产精品合色在线| 亚洲av成人不卡在线观看播放网| 欧美黑人精品巨大| 一个人观看的视频www高清免费观看 | 欧美中文综合在线视频| www日本在线高清视频| av中文乱码字幕在线| 老鸭窝网址在线观看| 国内精品久久久久精免费| 国产亚洲av高清不卡| 亚洲熟妇熟女久久| 亚洲成av人片免费观看| 国产片内射在线| 午夜两性在线视频| 日韩欧美精品v在线| 琪琪午夜伦伦电影理论片6080| 人成视频在线观看免费观看| 色综合欧美亚洲国产小说| 波多野结衣巨乳人妻| 国产爱豆传媒在线观看 | 久久久精品国产亚洲av高清涩受| 欧美一级毛片孕妇| 一本大道久久a久久精品| 丰满的人妻完整版| 国产成人精品无人区| 精品无人区乱码1区二区| 人妻丰满熟妇av一区二区三区| 又粗又爽又猛毛片免费看| 美女黄网站色视频| 久久久久久久精品吃奶| 怎么达到女性高潮| 手机成人av网站| 国产真人三级小视频在线观看| 午夜福利视频1000在线观看| 精品第一国产精品| 老司机靠b影院| 国产主播在线观看一区二区| 国产精品日韩av在线免费观看| 精品久久蜜臀av无| 无限看片的www在线观看| 成人av一区二区三区在线看| АⅤ资源中文在线天堂| x7x7x7水蜜桃| 日日爽夜夜爽网站| 亚洲国产欧美网| 亚洲激情在线av| 国产精品 国内视频| 老司机福利观看| 日日摸夜夜添夜夜添小说| 成人三级黄色视频| 国产免费av片在线观看野外av| 精品久久久久久久末码| 国产69精品久久久久777片 | 精品高清国产在线一区| 神马国产精品三级电影在线观看 | 午夜成年电影在线免费观看| 波多野结衣巨乳人妻| 草草在线视频免费看| 母亲3免费完整高清在线观看| 国产高清videossex| 国产精品久久久久久亚洲av鲁大| 国产精品九九99| 在线观看免费午夜福利视频| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 亚洲人成伊人成综合网2020| 久久久久久久久免费视频了| 在线观看美女被高潮喷水网站 | 国产精品久久久人人做人人爽| 亚洲九九香蕉| 99久久精品国产亚洲精品| 可以在线观看毛片的网站| 国内少妇人妻偷人精品xxx网站 | 亚洲 欧美 日韩 在线 免费| 久久人人精品亚洲av| 舔av片在线| 午夜精品久久久久久毛片777| 国产人伦9x9x在线观看| 婷婷精品国产亚洲av在线| 亚洲av日韩精品久久久久久密| 99久久国产精品久久久| 免费一级毛片在线播放高清视频| 国产主播在线观看一区二区| 亚洲一码二码三码区别大吗| 特大巨黑吊av在线直播| 国内少妇人妻偷人精品xxx网站 | 亚洲精品一卡2卡三卡4卡5卡| www日本黄色视频网| 日本 欧美在线| 成人18禁在线播放| 婷婷精品国产亚洲av在线| 欧美性猛交╳xxx乱大交人| 黄色片一级片一级黄色片| 黑人巨大精品欧美一区二区mp4| 久久亚洲精品不卡| 国产精华一区二区三区| 日本熟妇午夜| 国产在线观看jvid| 最近最新中文字幕大全免费视频| 免费看日本二区| 91字幕亚洲| 色综合站精品国产| 美女 人体艺术 gogo| 亚洲欧美日韩高清在线视频| 精品一区二区三区视频在线观看免费| 天堂动漫精品| 亚洲人与动物交配视频| 国产一级毛片七仙女欲春2| 国产av不卡久久| 国产亚洲欧美98| 18禁观看日本| 欧美一区二区国产精品久久精品 | 久久久久国内视频| 欧美 亚洲 国产 日韩一| 亚洲精品一区av在线观看| 国产精品久久久久久人妻精品电影| 国内精品久久久久精免费| 大型黄色视频在线免费观看| 国产精品九九99| 99久久国产精品久久久| 亚洲av中文字字幕乱码综合| 小说图片视频综合网站| 深夜精品福利| 黄片大片在线免费观看| 女生性感内裤真人,穿戴方法视频| 精品无人区乱码1区二区| 每晚都被弄得嗷嗷叫到高潮| 九色成人免费人妻av| 在线观看免费午夜福利视频| 国产精品98久久久久久宅男小说| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 欧美国产日韩亚洲一区| 亚洲自拍偷在线| 1024视频免费在线观看| 久久久精品大字幕| 香蕉av资源在线| 日韩有码中文字幕| 国产精品国产高清国产av| 国产精品av久久久久免费| 99久久精品国产亚洲精品| 欧美日韩精品网址| 我要搜黄色片| 最新美女视频免费是黄的| 亚洲 国产 在线| 午夜福利视频1000在线观看| 国产成人av激情在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看日韩欧美| 亚洲欧美日韩东京热| 亚洲18禁久久av| 最新美女视频免费是黄的| 日本黄色视频三级网站网址| or卡值多少钱| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 18禁美女被吸乳视频| netflix在线观看网站| 美女午夜性视频免费| 成人手机av| 欧美成人午夜精品| or卡值多少钱| 国产精品野战在线观看| 最好的美女福利视频网| 两个人视频免费观看高清| 国产视频内射| 制服人妻中文乱码| 老汉色∧v一级毛片| 视频区欧美日本亚洲| 亚洲欧美日韩高清在线视频| 神马国产精品三级电影在线观看 | a级毛片在线看网站| 欧美日韩瑟瑟在线播放| 一进一出好大好爽视频| 日韩有码中文字幕| 少妇粗大呻吟视频| cao死你这个sao货| 日本一本二区三区精品| 男女做爰动态图高潮gif福利片| 欧美日韩一级在线毛片| 国产精品98久久久久久宅男小说| 岛国在线免费视频观看| 日本一本二区三区精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成a人片在线一区二区| 在线国产一区二区在线| 又粗又爽又猛毛片免费看| av有码第一页| 久久久久精品国产欧美久久久| 老汉色∧v一级毛片| 精品国产乱码久久久久久男人| 真人做人爱边吃奶动态| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 国产精品久久久久久久电影 | 国产成人精品久久二区二区91| 精品欧美一区二区三区在线| 婷婷精品国产亚洲av在线| 一个人观看的视频www高清免费观看 | 精品熟女少妇八av免费久了| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 夜夜躁狠狠躁天天躁| 久久精品人妻少妇| 热99re8久久精品国产| 国产精品国产高清国产av| 欧美黑人欧美精品刺激| 久久久久久久久免费视频了| 欧美av亚洲av综合av国产av| 一边摸一边抽搐一进一小说| 亚洲人成网站在线播放欧美日韩| a级毛片a级免费在线| 国产成人av激情在线播放| 美女大奶头视频| 国产成人欧美在线观看| 久久久久免费精品人妻一区二区| 国产精品香港三级国产av潘金莲| 欧美人与性动交α欧美精品济南到| 欧美绝顶高潮抽搐喷水| 国产欧美日韩一区二区三| 日韩欧美在线二视频| 老司机靠b影院| 大型av网站在线播放| 白带黄色成豆腐渣| 制服丝袜大香蕉在线| www日本在线高清视频| 黄片小视频在线播放| 成年免费大片在线观看| 男女午夜视频在线观看| 一级作爱视频免费观看| 好男人电影高清在线观看| 久久亚洲真实| 国产激情偷乱视频一区二区| 色综合站精品国产| 国产精品一区二区免费欧美| 床上黄色一级片| 亚洲男人的天堂狠狠| 久久久国产欧美日韩av| 校园春色视频在线观看| 国产男靠女视频免费网站| 91麻豆av在线| 欧美又色又爽又黄视频| 日日干狠狠操夜夜爽| 精品福利观看| 香蕉丝袜av| 国内精品一区二区在线观看| 久久久久久大精品|