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

    基于壓電材料的薄膜聲學(xué)超材料隔聲性能研究*

    2019-08-27 06:56:26賀子厚趙靜波姚宏蔣娟娜陳鑫
    物理學(xué)報 2019年13期
    關(guān)鍵詞:特征頻率隔聲壓電

    賀子厚 趙靜波 姚宏 蔣娟娜 陳鑫

    1)(空軍工程大學(xué)研究生院,西安 710051)

    2)(空軍工程大學(xué)基礎(chǔ)部,西安 710051)

    1 引 言

    隨著科技的進(jìn)步,各類裝備不斷向著高速、輕質(zhì)、重型方面發(fā)展,隨之而來的噪聲污染也日益嚴(yán)重,影響到了人們?nèi)粘I詈脱b備的性能[1].尤其對于1000 Hz以下的低頻噪聲,由于其波長較長,難以在傳播中衰減,同時由于質(zhì)量定律[2]的限制,傳統(tǒng)的材料難以對其進(jìn)行有效地隔離.

    近年來,有許多學(xué)者對低頻噪聲的隔離問題進(jìn)行了研究.2000年,Liu等[3]提出了局域共振機制,基于這一機制的聲學(xué)超材料可以實現(xiàn)“小尺寸控制大波長”的功能.張思文和吳九匯[4]以及張帥等[5]進(jìn)一步討論了其低頻帶隙特性.同時,聲學(xué)超材料也在其他方面有廣泛的應(yīng)用.趙甜甜等[6]借助聲子晶體原理對超聲塑料焊接工具的橫向振動進(jìn)行了抑制.王莎和林書玉[7]基于二維聲子晶體原理對大尺寸夾心換能器進(jìn)行了優(yōu)化設(shè)計.張振方等[8]研究了內(nèi)插擴(kuò)張室聲子晶體管路帶隙特性.杜春陽等[9]對X形超阻尼局域共振聲子晶體梁彎曲振動帶隙特性進(jìn)行了研究.

    薄膜型聲學(xué)超材料由于輕質(zhì)、性能好等原因,在低頻噪聲隔離方面有著廣闊的應(yīng)用前景.2008年,Yang等[10]首先提出了“薄膜聲學(xué)超材料”的概念.2012年,Mei等[11]設(shè)計的“聲學(xué)暗材料”,由于薄膜結(jié)構(gòu)中軟硬介質(zhì)的相互作用所導(dǎo)致的局域共振與聲波振動的耦合,產(chǎn)生了一種“拍動”[12]的振動模式,實現(xiàn)了在一些頻率上對聲波的完全吸收.Chen等[13]對薄膜附加質(zhì)量塊結(jié)構(gòu)進(jìn)行了解析建模.Langfeldt等[14]優(yōu)化了薄膜聲學(xué)超材料傳輸損失的計算方法,大大減少了計算量.還有許多學(xué)者設(shè)計了諸多性能良好的薄膜聲學(xué)超材料[15?19].

    但是目前大部分的薄膜聲學(xué)超材料,結(jié)構(gòu)設(shè)計一旦固定,其性能便難以修改,滿足不了實際應(yīng)用的需求.為了實現(xiàn)材料性能的可調(diào)性,壓電分流阻尼技術(shù)被運用到了聲學(xué)超材料的結(jié)構(gòu)設(shè)計之中.壓電分流阻尼技術(shù)是一種起源于20世紀(jì)的振動控制技術(shù),其原理是通過壓電片與外接分流電路互相作用產(chǎn)生阻尼衰減效果[20].Chen等[21]通過在梁上附加壓電片,實現(xiàn)了一維梁結(jié)構(gòu)帶隙的可調(diào)性.Zhang等[22]則在鋁板上附加了壓電片,并實現(xiàn)了其隔聲性能的可調(diào)性.董亞科等[23]在外接電路中加入負(fù)電容,拓寬了壓電局域共振帶隙的頻帶.廖濤等[24]對二維壓電聲子晶體板帶隙可調(diào)性進(jìn)行了研究.孫煒海等[25]對含磁電彈夾層的壓電/壓磁聲子晶體帶隙特性進(jìn)行了研究.

    本文設(shè)計了一種附加壓電質(zhì)量塊的薄膜聲學(xué)超材料,使用有限元法計算了材料的特征頻率與20—1200 Hz頻段內(nèi)材料的傳輸損失曲線,并通過實驗驗證了有限元計算的可靠性,發(fā)現(xiàn)這種材料在此頻段內(nèi)隔聲性能優(yōu)良,并具有一個可調(diào)隔聲峰.

    2 結(jié)構(gòu)設(shè)計與計算方法

    通過以往的研究[26]可知,在聲學(xué)超材料中加入偏心質(zhì)量單元,并進(jìn)行合理的參數(shù)設(shè)置,可以顯著提升材料的性能.故此,在這里提出一種附加偏心質(zhì)量塊的薄膜聲學(xué)超材料,如圖1所示,該材料由壓電質(zhì)量塊嵌入緊張彈性膜,并固定在金屬外框上構(gòu)成,壓電材料的上下表面接有電極,與外部電路連接.其中A為壓電材料,其長為b=8 mm,寬為a=5 mm,厚度為e=2 mm; B為方形膜,其邊長為c=20 mm厚度為f=0.2 mm; C為正方形邊框,其外邊長為d=25 mm,厚度為f=0.2 mm,m=2 mm,n=7.5 mm分別為壓電質(zhì)量塊距離邊框上邊界與左邊界的距離.壓電質(zhì)量塊材料為PZT-5H,方形膜材料為硅橡膠,金屬外框材料為鋼.表1和表2為所涉及到的材料參數(shù),其中壓電材料參數(shù)來源于文獻(xiàn)[27],其中為彈性柔順系數(shù),為垂直膜面方向的自由介電常數(shù),d31為壓電應(yīng)變常數(shù).

    表1 壓電材料參數(shù)Table 1.Piezoelectric material parameters.

    表2 材料參數(shù)Table 2.Material parameters.

    在材料中,由于金屬框厚度較小而且受到的激勵垂直于版面,在計算中可將其看作薄板發(fā)生橫向振動,振動方程為

    其中,w為板的橫向位移,D0為金屬板的抗彎剛度,ρ為板的密度,h為板厚度;p(x,y,t)為薄板受到的橫向激勵,x,y和t分別為位置坐標(biāo)和時間.

    材料中膜的振動方程為

    對于壓電材料[28]來說,其壓電本構(gòu)方程,動力學(xué)方程和準(zhǔn)靜態(tài)電荷方程分別為:

    得到壓電晶體中的振動方程為

    本文通過Comsol有限元計算軟件計算結(jié)構(gòu)的傳輸損失曲線與特征頻率,將結(jié)構(gòu)四周的邊界條件設(shè)置為固定,將壓電片上邊界設(shè)置為懸浮電位,與電路模塊耦合,將壓電片下邊界設(shè)置為接地.在有限元計算中,設(shè)置D型壓電方程作為壓電片的本構(gòu)方程,D型壓電方程為

    其中:T1,S1分別為壓電片垂直于膜面方向的應(yīng)力與應(yīng)變;D3,E3分別為壓電材料上下表面的電場密度與電勢.

    外界電路結(jié)構(gòu)如圖1(a)所示,其中R為電阻,L為電感,二者串聯(lián).在計算結(jié)構(gòu)傳輸損失曲線時,構(gòu)建如圖2所示腔體結(jié)構(gòu).整個腔體長為400 mm,材料設(shè)置為空氣,空氣中的聲速為p0=343 m/s,空氣密度為c0=1.25 kg/m3.在腔體結(jié)構(gòu)的左側(cè)邊界垂直入射平面波激勵,并在結(jié)構(gòu)的右側(cè)邊界上拾取響應(yīng),計算兩者的差值得到結(jié)構(gòu)的傳輸損失(transmission loss,TL),其單位為分貝.

    其中Win與Wout分別為入射聲能與出射聲能,即為

    S1與S2分別為圖2腔體結(jié)構(gòu)的左側(cè)邊界與右側(cè)邊界.

    圖2 腔體結(jié)構(gòu)Fig.2.Cavity structure.

    3 結(jié)果分析

    在設(shè)置串聯(lián)電阻R=0 W,串聯(lián)電感L=200 H的情況下,使用上節(jié)所述方法計算結(jié)構(gòu)的特征頻率與傳輸損失曲線.表3所示為結(jié)構(gòu)前14階特征頻率與其模態(tài),圖3為20—1000 Hz頻率范圍的傳輸損失曲線.

    通過分析模態(tài)圖可知,包含質(zhì)量塊振動的模態(tài)都集中于500 Hz以下,500 Hz以上則展現(xiàn)了豐富的薄膜振動的模態(tài).材料的傳輸損失曲線如圖3實線所示.作為對照組,計算了未嵌入壓電質(zhì)量塊材料的傳輸損失曲線,如圖3虛線所示.

    材料在200—1200 Hz頻率范圍內(nèi)具有較為良好的隔聲性能.雖然在嵌入壓電質(zhì)量塊后,傳輸損失在100 Hz附近出現(xiàn)衰減,但總體上來看,材料的隔聲性能有了較大提升,并在458 Hz處出現(xiàn)壓電隔聲峰.

    結(jié)構(gòu)的共振點都已在圖3中表示處,尤其值得注意的是在216,1251 Hz頻率處出現(xiàn)非對稱的共振峰,即為傳輸損失的“突變”.在第4節(jié)中,將結(jié)合特征頻率與振動模式圖,進(jìn)一步討論結(jié)構(gòu)的隔聲機理.

    表3 模態(tài)圖Table 3.Modal diagram.

    圖3 傳輸損失曲線Fig.3.Transmission loss curve.

    觀察傳輸損失曲線可以發(fā)現(xiàn),一些特征頻率處并未出現(xiàn)共振點,這是由于在計算傳輸損失曲線時,僅僅考慮了平面波垂直膜面入射的情況,而一些模態(tài)難以與此方向上的行波發(fā)生耦合[29]作用,故此不出現(xiàn)共振點.此外,大量的工程實踐也證明,并不是所有的共振模態(tài)都可以被激起.

    4 實驗驗證

    為了驗證有限元計算的可靠性,這里使用AWA6290T型吸聲測量系統(tǒng)測試了聲學(xué)超材料結(jié)構(gòu)的隔聲性能,由于此測量系統(tǒng)的阻抗管呈圓柱形,故設(shè)計圓形結(jié)構(gòu)如圖4所示.

    圖4 實驗示意圖(a)樣件結(jié)構(gòu);(b)實驗裝置;(c)樣件實物圖Fig.4.Experimental schematic diagram:(a)Sample structure;(b)experimental facility;(c)physical samples.

    圖4中環(huán)形外框材質(zhì)為類 丙烯腈-丁二烯-苯乙烯(ABS)白色SLA(stereo lithography apparatus)樹脂,其密度為1.13 g/cm2,楊氏模量為2600 MPa,泊松比為0.37,薄膜材質(zhì)為硅橡膠,張力為1.5 MPa,所貼壓電材料為鍍銀電極PZT-5H壓電陶瓷附加銅質(zhì)基板.壓電陶瓷與硅橡膠間使用氰基丙烯酸酯粘合劑連接.壓電片外部串聯(lián)0—1 H可調(diào)電感與10 W電阻.由于在電路中所需電感值較大,這里使用Antoniou’s電路來模擬大電感(參見圖5).由于采樣精度有限,通過調(diào)節(jié)電感可以更方便地監(jiān)測到壓電隔聲峰的存在.

    圖5 大電感(a)電路圖;(b)實物圖Fig.5.Large inductance:(a)Circuit diagram;(b)physical diagram.

    分別使用實驗方法與有限元計算方法得到結(jié)構(gòu)的傳輸損失曲線,如圖6所示.實驗數(shù)據(jù)由500次采樣求均值得到.

    圖6 傳輸損失曲線Fig.6.Transmission loss curve.

    可以看出,實驗雖然存在一些誤差,但兩支曲線基本吻合,驗證了有限計算的真實性.在實驗中,通過調(diào)節(jié)電感,可以在630 Hz檢測到尖銳隔聲峰,與有限元計算得到的583 Hz壓電隔聲峰相印證,驗證了壓電隔聲峰的真實性.

    5 隔聲機理分析

    5.1 等效模型

    一般來說,薄膜聲學(xué)超材料的首階特征頻率對其低頻隔聲性能有極大的影響,更高的首階特征頻率,往往意味著更好的低頻隔聲性能.雖然本文所設(shè)計的結(jié)構(gòu)性能較好,但是由于其結(jié)構(gòu)形狀較為不規(guī)則,難以進(jìn)行分析.為進(jìn)一步揭示薄膜附質(zhì)量塊結(jié)構(gòu)的隔聲機理,這里對實驗中的結(jié)構(gòu)進(jìn)行了解析建模,最后構(gòu)建了“彈簧-振子”等效模型來估算其首階特征頻率.同理,由于此簡單與前文所設(shè)計的結(jié)構(gòu)相似,有理由相信二者隔聲機理相同.

    由于首階特征頻率不涉及壓電材料.可以將結(jié)構(gòu)簡化為如圖7所示,其中小圓部分為薄膜,材質(zhì)為硅橡膠,其余部分為附加質(zhì)量塊,薄膜張力為44 N/m,邊緣固定.其結(jié)構(gòu)尺寸如圖7所示.

    圖7 簡化結(jié)構(gòu)示意圖Fig.7.Simplified structure sketch.

    設(shè)一不附加質(zhì)量塊的薄膜,其振動方程為

    其中ξ為薄膜的縱向位移;為二維坐標(biāo)拉普拉斯算符;其中,T為膜表面張力,s為膜的面密度.

    已知首階共振模態(tài)為圓對稱情形,故可將振動方程的解寫為徑向距離g與時間t的函數(shù):

    其中 J0為零階貝塞爾函數(shù),

    設(shè)圓膜的半徑為R,則邊界條件為ξ(γ=R)=0.代入解中可得,J0(κR)=0.

    在MATLAB調(diào)用Besselj函數(shù),可以畫出貝塞爾函數(shù)在[0,1]區(qū)間上的圖線,如圖8所示.

    圖8 貝塞爾函數(shù)曲線Fig.8.Bessel function curve.

    通過最小二乘法搜索得到J0(μ)=0的解,精確到三位小數(shù)得到μ=2.405.根據(jù)κR=μ,可得首階特征頻率f為

    在圓膜上取微元,其徑向距離為(r,r+dr),其動能為dE:

    將運動方程代入(11)式并取其實部可得此模態(tài)下的平均動能為

    其中J1為一階貝塞爾函數(shù).

    若將此薄膜的振動等效為“彈簧—振子”結(jié)構(gòu),設(shè)振子的位移為薄膜中心的位移.可以得到振子速度v的值為

    設(shè)M為系統(tǒng)等效質(zhì)量,則系統(tǒng)在一個周期內(nèi)的平均動能為

    其中m為薄膜的質(zhì)量,從而得到如圖9(a)所示“彈簧—振子”系統(tǒng),其中K為等效彈簧的質(zhì)量.此系統(tǒng)的特征頻率令可得

    現(xiàn)考慮在原膜中心附加質(zhì)量塊的情況,其等效模型如圖9(b)所示,即在原先質(zhì)量M之上在附加質(zhì)量Madd.在附加質(zhì)量塊之后,可知系統(tǒng)的等效彈性模量也會發(fā)生變化,但可以近似地認(rèn)為K1=K.

    圖9 等效模型示意圖Fig.9.Schematic diagram of equivalent model.

    故此,簡化結(jié)構(gòu)的首階特征頻率為

    分別使用有限元法與等效模型方法,計算了簡化模型的首階特征頻率,取附加質(zhì)量塊的密度為1000—5000 kg/m2,結(jié)果如圖10所示,兩種方法可以基本吻合,但是仍存在一定誤差,這是由于簡化模型所附加的質(zhì)量塊底面積較大,對薄膜的等效彈性模量產(chǎn)生了影響.

    圖10 首階特征頻率Fig.10.First natural frequency.

    5.2 模態(tài)分析

    為了進(jìn)一步研究結(jié)構(gòu)的隔聲機理,這里計算了各隔聲峰、隔聲谷和傳輸損失“突變”處的振動模式圖.

    根據(jù)如圖11所示的消聲原理,首先對隔聲峰出現(xiàn)的機理進(jìn)行分析.圖12(a)所示為結(jié)構(gòu)在185 Hz處隔聲峰處的振動模式圖,此處的傳輸損失高達(dá)69.5 dB.當(dāng)材料受到此頻率聲波激勵時,壓電質(zhì)量塊與薄膜反向振動,呈現(xiàn)出“拍動”模式,由于二者振幅相等,相位相反,使得在遠(yuǎn)場處,薄膜振動產(chǎn)生的聲波與壓電片振動產(chǎn)生的聲波干涉相消,從而使得入射聲波得到極大的衰減.同時材料本身也發(fā)生了極大的形變,將聲波的能量局域在材料之中.

    圖11 消聲原理圖Fig.11.Anechoic schematic diagram.

    這一振動模式未體現(xiàn)在共振模態(tài)中,之所以出現(xiàn)這樣的現(xiàn)象,一是由于在計算傳輸損失時考慮了空氣阻尼,二是由于聲波入射的方向垂直與膜面,致使多個平行于膜面的模態(tài)難以展現(xiàn).但這種振動模式仍可以看作是多個共振模態(tài)的疊加.

    位于485.62 Hz的壓電隔聲峰的帶寬只有0.59 Hz,這一頻率下的傳輸損失為23.65 dB,相較于其左右1 Hz頻率范圍內(nèi),此處傳輸損失上升了12 dB,其振動模式圖如圖12(b)所示.從圖中可以看出,此模態(tài)下壓電質(zhì)量塊橫向壓縮形變,行波的能量被壓電片所吸收.

    圖12 隔聲峰處的振動模式圖(a)185 Hz;(b)485.6 Hz;(c)896 HzFig.12.Vibration mode diagram at sound insulation peak:(a)185 Hz;(b)485.6 Hz;(c)969 Hz.

    在材料中,壓電片可看作電容(電容值可通過Comsol有限元軟件可求得)連接外接電路中的電感與電阻,共同組成了LC振蕩電路.由D型壓電方程可知,當(dāng)壓電片受到一定頻率聲波激勵時,其上下表面的電荷密度發(fā)生同頻率的改變.當(dāng)壓電片兩端電荷密度變化頻率接近LC振蕩電路的共振頻率時,電路就會發(fā)生強烈的振蕩,聲波的能量被轉(zhuǎn)化為電能并消耗.這一隔聲峰的存在,是本結(jié)構(gòu)的重要特性.

    圖12(c)為896 Hz隔聲峰處結(jié)構(gòu)的振動模式圖,對應(yīng)于第13階特征頻率,此振動模式表現(xiàn)為薄膜分兩部分發(fā)生“振幅相等,相位相反”的振動,形成“拍動”模式,這一隔聲峰形成機理與185 Hz隔聲峰類似.

    圖13(a)所示為結(jié)構(gòu)在115 Hz隔聲谷處的振動模式圖,對應(yīng)于首階特征頻率,圖13(b)為對照組在457 Hz隔聲谷處的振動模式圖.從以往的研究可知,結(jié)構(gòu)的首階特征頻率會對其低頻隔聲性能產(chǎn)生極大的影響.在附加質(zhì)量塊的結(jié)構(gòu)中,115 Hz隔聲谷處的傳輸損失幾乎下降至0,此時材料發(fā)生強烈的橫向共振,聲波的大部分能量得以透過薄膜.未附加質(zhì)量塊的薄膜結(jié)構(gòu)出現(xiàn)隔聲谷的機理與之類似,但是由于未附加質(zhì)量塊薄膜的等效質(zhì)量更小,根據(jù)公式其首階特征頻率會更高,根據(jù)有限元計算可知,其首階特征頻率為456 Hz,因此,未附加質(zhì)量塊薄膜在低頻范圍內(nèi)隔聲量衰減得更慢.但也由于未附加壓電材料,其隔聲性能較為平庸.

    圖13 隔聲谷處的振動模式圖(a)115 Hz;(b)457HzFig.13.Vibration mode diagram at sound insulation peak:(a)115 Hz;(b)457Hz.

    圖14(a)所示為687 Hz隔聲谷處結(jié)構(gòu)的振動模式圖,對應(yīng)于第8階特征頻率,可以看出,其與圖13(b)的振動模式類似,隔聲谷產(chǎn)生機理也類似,但由于壓電片的存在,膜的振動面積較小,使其頻率更升高.圖14(b)所示為969 Hz隔聲谷處結(jié)構(gòu)的振動模式圖,可以看出,其振動模式與896 Hz隔聲峰處的振動模式及其類似,表現(xiàn)為壓電片長邊側(cè)兩邊與短邊一側(cè)的振動,但是在896 Hz隔聲峰處,此兩部分的振幅與膜的振動面積大致相等,故此在遠(yuǎn)場處其波動可以相互抵消,但是在969 Hz隔聲谷處,長邊兩側(cè)振動的振幅與膜的振動面積要遠(yuǎn)大于其短邊側(cè)的振幅,使其在遠(yuǎn)場處其波動難以干涉相消,故此在傳輸損失曲線上表現(xiàn)為隔聲谷.與之機理類似的還有處于1129—1136 Hz頻率范圍內(nèi)的突變,圖14(c)為1129 Hz處的振動模式圖,圖14(d)為1136 Hz處的振動模式圖,在1129 Hz處結(jié)構(gòu)表現(xiàn)高達(dá)33 dB的傳輸損失,緊接著在1136 Hz急劇下降,造成這樣突變的原因也是由于膜上不同部位振動的細(xì)微差別導(dǎo)致,當(dāng)不同振動相位相反的兩部分的振動面積與振幅大致相等時,則表現(xiàn)為傳輸損失上升,反之,則表現(xiàn)為傳輸損失的下降.

    圖14 隔聲谷與傳輸損失突變處的振動模式圖(a)687 Hz;(b)969 Hz;(c)1129 Hz;(d)1136 HzFig.14.Vibration mode diagram at TL peak and TL sudden change:(a)687 Hz;(b)969 Hz;(c)1129 Hz;(d)1136 Hz.

    在229—235 Hz附近也出現(xiàn)了一次突變,其中229 Hz的振動模式圖如圖15(a)所示,235 Hz的振動模式如圖15(b)所示.從傳輸損失曲線中可以看出在229 Hz時,結(jié)構(gòu)傳輸損失迅速衰減,但在235 Hz時急劇升高.從圖中可以看到,兩者的振動模式圖極其類似,都表現(xiàn)為壓電片以平行于短邊的線為轉(zhuǎn)軸的轉(zhuǎn)動振動,但是也可以明顯看出,在235 Hz時,結(jié)構(gòu)的轉(zhuǎn)軸更加靠近壓電片的中軸線,這意味著,在235 Hz時,結(jié)構(gòu)的振動更加接近于“拍動”的模式,而在228 Hz時,壓電片兩端的振幅不同,使得在遠(yuǎn)場處二者難以相互抵消,傳輸損失下降.

    5.3 Fano共振對結(jié)構(gòu)隔聲性能的影響

    圖15 傳輸損失突變處的振動模式圖(a)229 Hz;(b)235 HzFig.15.Vibration mode diagram at TL sudden change:(a)229 Hz;(b)235 Hz.

    綜合分析特征頻率與傳輸損失曲線可以發(fā)現(xiàn),結(jié)構(gòu)在各特征頻率處的傳輸損失表現(xiàn)出不同的特征,例如出現(xiàn)“隔聲峰”、“隔聲谷”,以及“突變”等.5.2節(jié)已經(jīng)從模態(tài)分析的角度討論了其產(chǎn)生機理,這里從Fano共振的角度,解釋在不同特征頻率處傳輸損失曲線表現(xiàn)出不同特征的原因.

    根據(jù)共振機理不同,共振可以分為Lorentz共振與Fano共振.其中Lorentz共振只涉及一種模態(tài),在頻譜中呈現(xiàn)對稱的譜線形狀,如前文所述的“隔聲峰”、“隔聲谷”等一類情況.當(dāng)涉及多種模態(tài)的耦合,則會出現(xiàn)Fano共振,即在頻譜中表現(xiàn)出非對稱的譜線形狀,例如前文所述的“突變”.

    Fano resonance(FR)概念起源于量子物理領(lǐng)域,于1961年被Ugo Fano提出[30],并得到了Fano諧振的公式:

    其中I(ω)為頻率響應(yīng);q為Fano參數(shù)(Fano parameter),其定義為兩種互相干涉模態(tài)的強度比值,影響著Fano共振頻譜的非對稱性;e為歸一化的調(diào)諧頻率.圖16所示為Fano參數(shù)不同時,特征頻率附近的頻率響應(yīng)曲線.從圖中觀察可得,當(dāng)q=0時,在特征頻率附近的頻率響應(yīng)呈現(xiàn)Lorentz共振模式,這時由于q=0時,兩種模態(tài)有一種的強度為0,故此不存在干涉的情況; 當(dāng)q=± 1時,此時兩種模態(tài)的能量最為接近,頻率響應(yīng)呈現(xiàn)標(biāo)準(zhǔn)的Fano共振的非對稱特性.同理,當(dāng)q值逐漸趨向于無窮時,頻率響應(yīng)又會變?yōu)長orentz共振模式.

    后來Fano共振的概念拓展到光學(xué)領(lǐng)域[31],廣泛應(yīng)用于光子晶體[32]等方面的研究中,被用于描述兩個模態(tài)耦合的情況下,材料對入射波的響應(yīng)情況.由于兩個模態(tài)的相位響應(yīng)不同,當(dāng)二者互相影響時,其譜線并不是簡單的強度疊加,當(dāng)兩種模態(tài)存在相位差時,會出現(xiàn)干涉相消的情況,在譜線中表現(xiàn)為谷值,或者當(dāng)兩種振動模態(tài)的相位響應(yīng)隨入射波頻率變化的速率差別較大時,則會出現(xiàn)類似于上文所述的隔聲量“突變”的情況.

    圖16 Fano共振Fig.16.Fano resonance.

    在對薄膜聲學(xué)超材料的研究過程中發(fā)現(xiàn),當(dāng)薄膜受到較高頻率聲波激勵時,薄膜的振動往往被“分割”開來,例如第10階模態(tài)、第14階模態(tài),薄膜的各個部分表現(xiàn)為相互獨立的振動,這種獨立的振動可以看作分立的模態(tài),在薄膜后的傳播過程中,各分立模態(tài)發(fā)生干涉作用.

    類比于光學(xué)領(lǐng)域,由于各振動模態(tài)的Fano參數(shù)不同,相位響應(yīng)也不同,導(dǎo)致在不同特征頻率處展現(xiàn)出不同的共振模式.例如在各隔聲峰與隔聲谷處表現(xiàn)為Lorentz共振模式,在216,1251 Hz表現(xiàn)為非對稱的Fano共振模式,同時896 Hz隔聲峰與969 Hz隔聲谷也可以看作Fano參數(shù)適中的Fano共振.

    6 結(jié)構(gòu)參數(shù)對結(jié)構(gòu)隔聲性能的影響

    6.1 壓電片偏心量對材料性能的影響

    圖17所示為壓電片偏心量m不同時,20—1200 Hz頻段內(nèi)結(jié)構(gòu)的傳輸損失曲線.從圖中可知,壓電片的偏心量對結(jié)構(gòu)的低頻隔聲性能影響較小,但是對結(jié)構(gòu)的高頻隔聲性能影響較大,并且,無論m值如何變化,在458 Hz處,結(jié)構(gòu)的隔聲量都會有約12 dB的上升.

    圖17 傳輸損失曲線m=0.001,0.004,0.006 mFig.17.TL curve m=0.001,0.004,0.006 m.

    之所以傳輸損失曲線出現(xiàn)這樣的變化,是由于低頻特征頻率對m值較敏感,高頻特征頻率對m值較為不敏感.圖18所示為m=0.001,0.004,0.006 m時材料的各階特征頻率.

    圖18 特征頻率Fig.18.Eigen frequencies.

    如圖17所示,壓電片的偏心量對300 Hz以下的特征頻率影響甚小,這是因為300 Hz以下的共振模態(tài)表現(xiàn)為“彈簧-振子”模式,由于膜的面積較大,壓電片的偏心量對“彈簧-振子”系統(tǒng)的等效剛度影響不大.例如圖19所示為m值不同時第五階特征頻率的模態(tài)圖,其振動表現(xiàn)為壓電片帶動周圍一部分膜進(jìn)行振動,壓電片的位置對其影響較小.

    但是嵌入壓電片的位置對500 Hz以上的特征頻率影響較大.從前文的分析可知,500 Hz以上的共振模態(tài)表現(xiàn)為薄膜的振動,嵌入薄膜的壓電片則會在一定程度上 “分割”薄膜,改變薄膜振動的位置與面積,從而影響材料的特征頻率.

    從圖18中也可以看到,m值改變對結(jié)構(gòu)的第七階特征頻率影響很小.這是由于第7階特征頻率是由壓電片外接的LC電路所決定的,與壓電片的位置無關(guān),使得無論m值如何變化,在458 Hz附近總會出現(xiàn)隔聲峰.

    6.2 電路參數(shù)對隔聲量的影響

    圖19 第五階共振模態(tài)(a)m=0.002 m;(b)m=0.004 m;(c)m=0.006 mFig.19.Fifth order vibration:(a)m=0.002 m;(b)m=0.004 m;(c)m=0.006 m.

    圖20 電路參數(shù)不同時隔聲量的變化(a)不同電阻;(b)不同電感Fig.20.TL with different circuit parameters:(a)Different resistors;(b)different inductances.

    對于一般的薄膜聲學(xué)超材料而言,一旦其結(jié)構(gòu)固定,其性能便也同時固定.但是本材料中加入了壓電材料,可以通過改變外接電路的系數(shù)來調(diào)節(jié)材料的性能,使其能夠更好地滿足需求.圖20(a)和圖20(b)分別為外接電阻為1,10,100,1000 W時與外接電感為200,201,202,203 H時結(jié)構(gòu)的隔聲量.

    從圖20(a)中可以看出,隨著電阻的增大,隔聲峰的峰值逐漸減小(此時設(shè)定L=200 H),而隔聲峰的位置不發(fā)生改變.峰值減小的原因主要是電阻的增大使電磁振蕩作用減弱.從圖20(b)中可以看出,隨著電感增大(此時R=1 W),隔聲峰逐漸向低頻移動.這是由于LC振蕩電路的諧振頻率fr為

    隨著電感值的增大,諧振頻率fr逐漸變小.

    7 結(jié) 論

    本文設(shè)計了一種基于壓電材料的薄膜聲學(xué)超材料,得到結(jié)論如下:

    1)使用有限元方法計算了結(jié)構(gòu)的特征頻率與20—1200 Hz頻段的隔聲量,發(fā)現(xiàn)材料此頻率范圍內(nèi)隔聲性能良好,存在兩個50 dB以上的隔聲峰和一個可調(diào)隔聲峰,并在實驗中驗證了有限元計算的真實性;

    2)建立了簡化結(jié)構(gòu)首階共振模態(tài)等效模型,并且使用有限元法驗證了其合理性; 綜合分析材料的傳輸損失曲線與共振模態(tài),討論了材料的隔聲機理,并通過Fano共振理論解釋了傳輸損失“突變”出現(xiàn)的機理;

    3)探究了壓電片偏心量對結(jié)構(gòu)性能的影響,通過調(diào)整電路參數(shù),實現(xiàn)了材料性能的可調(diào)性.

    猜你喜歡
    特征頻率隔聲壓電
    常用樓板隔聲材料的撞擊聲隔聲性能測試分析
    綠色建筑(2021年1期)2021-07-19 02:15:04
    瓷磚檢測機器人的聲音信號處理
    一種隔聲軟簾的研制及應(yīng)用
    《壓電與聲光》征稿啟事
    壓電與聲光(2019年1期)2019-02-22 09:46:06
    光學(xué)波前參數(shù)的分析評價方法研究
    基于振動信號特征頻率的數(shù)控車床故障辨識方法
    新型壓電疊堆泵設(shè)計及仿真
    基于小波去噪和EMD算法在齒輪故障檢測中的應(yīng)用
    基于混響室—半消聲室的V型聲屏障隔聲測試與分析
    基于壓電激振的彈性模量測量方法
    在线天堂最新版资源| 一个人看视频在线观看www免费| 国产深夜福利视频在线观看| 中文乱码字字幕精品一区二区三区| 男女边吃奶边做爰视频| 欧美日韩视频精品一区| 国产男女超爽视频在线观看| 日韩欧美一区视频在线观看 | 亚洲欧洲国产日韩| 日本免费在线观看一区| 日韩成人伦理影院| 精品亚洲乱码少妇综合久久| 国产精品嫩草影院av在线观看| 超碰av人人做人人爽久久| 啦啦啦在线观看免费高清www| 亚洲精品乱码久久久久久按摩| 女性生殖器流出的白浆| 国产av一区二区精品久久 | 久热这里只有精品99| 干丝袜人妻中文字幕| 午夜福利视频精品| 国产精品精品国产色婷婷| 亚州av有码| 国产精品国产三级专区第一集| 乱码一卡2卡4卡精品| 亚洲va在线va天堂va国产| 国产精品精品国产色婷婷| 美女国产视频在线观看| 边亲边吃奶的免费视频| 大香蕉久久网| 中文资源天堂在线| 欧美97在线视频| 伊人久久国产一区二区| 如何舔出高潮| 国产黄色视频一区二区在线观看| 青青草视频在线视频观看| 丰满人妻一区二区三区视频av| 国产精品国产三级专区第一集| 免费观看av网站的网址| 老司机影院毛片| 熟女av电影| 国产日韩欧美亚洲二区| 精品99又大又爽又粗少妇毛片| 99热6这里只有精品| 亚洲精品国产av成人精品| 天美传媒精品一区二区| 一级毛片 在线播放| 哪个播放器可以免费观看大片| 日韩av不卡免费在线播放| 男女免费视频国产| 成人亚洲欧美一区二区av| 99国产精品免费福利视频| 国产又色又爽无遮挡免| 成人午夜精彩视频在线观看| 免费黄频网站在线观看国产| 青春草国产在线视频| 国产中年淑女户外野战色| 毛片女人毛片| 亚洲综合精品二区| 午夜福利视频精品| 亚洲欧美日韩卡通动漫| 日本爱情动作片www.在线观看| 搡老乐熟女国产| 高清日韩中文字幕在线| 搡女人真爽免费视频火全软件| 精品久久久久久久久亚洲| 国产极品天堂在线| 边亲边吃奶的免费视频| 精品国产一区二区三区久久久樱花 | 精品久久国产蜜桃| 一区二区av电影网| 日韩精品有码人妻一区| 亚洲欧美日韩无卡精品| 一区二区av电影网| 国产一区亚洲一区在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产av新网站| 99热网站在线观看| 日韩大片免费观看网站| 国产精品av视频在线免费观看| 久久99精品国语久久久| 男女啪啪激烈高潮av片| 日韩欧美 国产精品| 日本-黄色视频高清免费观看| 性高湖久久久久久久久免费观看| 亚洲国产精品999| 97超碰精品成人国产| 男人添女人高潮全过程视频| 国产日韩欧美亚洲二区| 欧美高清性xxxxhd video| 成年av动漫网址| 欧美成人午夜免费资源| 国产精品蜜桃在线观看| 如何舔出高潮| 女的被弄到高潮叫床怎么办| 乱系列少妇在线播放| 亚洲欧美日韩卡通动漫| 男女国产视频网站| 我要看日韩黄色一级片| 免费av中文字幕在线| 欧美另类一区| 在线观看美女被高潮喷水网站| 免费观看无遮挡的男女| 国产综合精华液| 久久久久久久精品精品| 国产综合精华液| 色视频在线一区二区三区| 一级毛片久久久久久久久女| 综合色丁香网| 日本爱情动作片www.在线观看| 我要看日韩黄色一级片| 亚洲aⅴ乱码一区二区在线播放| av在线播放精品| 在线观看免费高清a一片| 啦啦啦视频在线资源免费观看| 亚洲,一卡二卡三卡| 国产免费视频播放在线视频| 91精品伊人久久大香线蕉| 日本爱情动作片www.在线观看| 少妇被粗大猛烈的视频| 久久99精品国语久久久| 免费高清在线观看视频在线观看| av播播在线观看一区| 欧美日韩在线观看h| 亚洲一区二区三区欧美精品| 久久99精品国语久久久| 久久99蜜桃精品久久| 精品国产乱码久久久久久小说| 3wmmmm亚洲av在线观看| 视频区图区小说| 国产一区亚洲一区在线观看| 中文精品一卡2卡3卡4更新| 看非洲黑人一级黄片| 亚洲美女视频黄频| 亚洲av.av天堂| 精品午夜福利在线看| 中国三级夫妇交换| 日本vs欧美在线观看视频 | 欧美亚洲 丝袜 人妻 在线| 免费高清在线观看视频在线观看| 国产精品av视频在线免费观看| 亚洲精品色激情综合| 多毛熟女@视频| 国国产精品蜜臀av免费| 欧美一区二区亚洲| 欧美zozozo另类| 99久国产av精品国产电影| 国产精品久久久久久精品电影小说 | 日韩一区二区三区影片| 三级国产精品片| 国产91av在线免费观看| 尤物成人国产欧美一区二区三区| 免费人成在线观看视频色| 国产精品国产av在线观看| 国产成人freesex在线| 看十八女毛片水多多多| .国产精品久久| 久久人妻熟女aⅴ| 亚洲精品视频女| 五月伊人婷婷丁香| 一本—道久久a久久精品蜜桃钙片| 搡老乐熟女国产| 亚洲婷婷狠狠爱综合网| 国产淫片久久久久久久久| 国产成人aa在线观看| 中文字幕免费在线视频6| 亚洲人成网站高清观看| 久久av网站| av.在线天堂| 男女啪啪激烈高潮av片| 美女主播在线视频| 亚洲自偷自拍三级| 一级毛片久久久久久久久女| 高清不卡的av网站| 国产在视频线精品| 日韩一区二区三区影片| 久久精品国产鲁丝片午夜精品| 久热这里只有精品99| 波野结衣二区三区在线| 国产高潮美女av| 国产精品国产三级专区第一集| 国产一区二区在线观看日韩| 亚洲av福利一区| 久久影院123| 六月丁香七月| 一本一本综合久久| 久久精品久久久久久久性| 新久久久久国产一级毛片| 日韩一区二区视频免费看| 五月开心婷婷网| 亚洲av日韩在线播放| 这个男人来自地球电影免费观看 | 26uuu在线亚洲综合色| 啦啦啦在线观看免费高清www| 亚洲精品久久午夜乱码| 国内精品宾馆在线| 亚洲精品乱码久久久v下载方式| 久久精品熟女亚洲av麻豆精品| 夜夜看夜夜爽夜夜摸| 国产女主播在线喷水免费视频网站| 少妇人妻一区二区三区视频| 国产男人的电影天堂91| 亚洲色图综合在线观看| 日本一二三区视频观看| 热re99久久精品国产66热6| 精品一区二区三卡| 高清视频免费观看一区二区| 免费观看在线日韩| av播播在线观看一区| 一本色道久久久久久精品综合| 免费av中文字幕在线| 日日摸夜夜添夜夜爱| 亚洲欧美清纯卡通| 国产精品久久久久久久久免| 黄色欧美视频在线观看| 麻豆成人av视频| 欧美xxⅹ黑人| 久久国产乱子免费精品| 国产成人免费观看mmmm| 亚洲av电影在线观看一区二区三区| 免费观看无遮挡的男女| 亚洲国产精品成人久久小说| 全区人妻精品视频| 91久久精品电影网| 亚洲精品日韩av片在线观看| 少妇的逼好多水| 国产免费视频播放在线视频| 日本欧美国产在线视频| 嘟嘟电影网在线观看| av专区在线播放| 国产在线免费精品| 亚洲精品aⅴ在线观看| 欧美日韩视频精品一区| 国产国拍精品亚洲av在线观看| av卡一久久| 国产成人aa在线观看| 又黄又爽又刺激的免费视频.| 国产探花极品一区二区| 午夜日本视频在线| 国产成人免费观看mmmm| 校园人妻丝袜中文字幕| 在现免费观看毛片| 在线观看人妻少妇| 欧美激情极品国产一区二区三区 | av天堂中文字幕网| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 超碰av人人做人人爽久久| 男人和女人高潮做爰伦理| 成年免费大片在线观看| 青春草视频在线免费观看| 成人亚洲精品一区在线观看 | 亚洲国产欧美人成| 一区二区三区乱码不卡18| 国产v大片淫在线免费观看| 国产精品人妻久久久久久| av.在线天堂| 少妇被粗大猛烈的视频| 国产亚洲最大av| 大香蕉97超碰在线| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 自拍偷自拍亚洲精品老妇| 2018国产大陆天天弄谢| 欧美成人一区二区免费高清观看| 国产精品国产三级国产专区5o| 国产成人精品一,二区| 欧美精品一区二区大全| 有码 亚洲区| a级毛色黄片| 女人久久www免费人成看片| 自拍欧美九色日韩亚洲蝌蚪91 | 成年免费大片在线观看| 老司机影院毛片| 一级片'在线观看视频| 最近2019中文字幕mv第一页| 亚洲在久久综合| 精品久久久久久电影网| 少妇裸体淫交视频免费看高清| 五月天丁香电影| 极品少妇高潮喷水抽搐| 九草在线视频观看| 日日啪夜夜撸| 少妇的逼好多水| 午夜福利影视在线免费观看| 狠狠精品人妻久久久久久综合| 91精品一卡2卡3卡4卡| 午夜福利高清视频| 最黄视频免费看| 少妇的逼好多水| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 国产免费福利视频在线观看| 亚洲国产精品专区欧美| 天堂8中文在线网| 少妇熟女欧美另类| 久久国产精品男人的天堂亚洲 | 日本免费在线观看一区| 免费播放大片免费观看视频在线观看| 视频中文字幕在线观看| 成年人午夜在线观看视频| 亚洲人成网站高清观看| 精品熟女少妇av免费看| 久久99热这里只有精品18| 日本av手机在线免费观看| 久久99热6这里只有精品| 中文资源天堂在线| 亚洲激情五月婷婷啪啪| 成人影院久久| 久久99热6这里只有精品| 欧美日本视频| 国产黄色免费在线视频| 国国产精品蜜臀av免费| 国产片特级美女逼逼视频| 又黄又爽又刺激的免费视频.| 久久久久久久久久人人人人人人| 久久99热这里只有精品18| 六月丁香七月| 美女视频免费永久观看网站| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 久久 成人 亚洲| 午夜精品国产一区二区电影| 国产成人精品福利久久| 男女下面进入的视频免费午夜| 亚洲色图综合在线观看| 91精品国产九色| 亚洲欧美日韩东京热| 亚洲精品,欧美精品| 国产精品女同一区二区软件| 97精品久久久久久久久久精品| 美女中出高潮动态图| 欧美精品一区二区大全| 精品久久久久久久久av| 一本一本综合久久| 99久久人妻综合| 精品熟女少妇av免费看| 丰满人妻一区二区三区视频av| 2021少妇久久久久久久久久久| 免费人成在线观看视频色| 日韩一本色道免费dvd| 国产男女超爽视频在线观看| 中文字幕av成人在线电影| 色哟哟·www| 51国产日韩欧美| 熟女av电影| 亚洲精品久久久久久婷婷小说| 成人影院久久| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 久久久久精品久久久久真实原创| 啦啦啦视频在线资源免费观看| 18禁动态无遮挡网站| 久久久久久久久久久免费av| 少妇人妻久久综合中文| 精品久久久精品久久久| 国产日韩欧美在线精品| 22中文网久久字幕| 丰满迷人的少妇在线观看| 男人狂女人下面高潮的视频| 成人高潮视频无遮挡免费网站| 日韩中字成人| 嫩草影院入口| 丝袜喷水一区| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 欧美成人午夜免费资源| 亚洲欧美精品自产自拍| 多毛熟女@视频| 欧美成人午夜免费资源| 久久人妻熟女aⅴ| 国产国拍精品亚洲av在线观看| 国产中年淑女户外野战色| 亚洲久久久国产精品| 看非洲黑人一级黄片| 日韩免费高清中文字幕av| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片 | 欧美高清成人免费视频www| 欧美日韩精品成人综合77777| 一二三四中文在线观看免费高清| 精品99又大又爽又粗少妇毛片| av天堂中文字幕网| 黑人高潮一二区| 国产片特级美女逼逼视频| 国产又色又爽无遮挡免| 亚洲精品日韩av片在线观看| 久久国产亚洲av麻豆专区| 又黄又爽又刺激的免费视频.| 国产永久视频网站| 国产综合精华液| 久久久久久人妻| 欧美精品一区二区大全| 一级二级三级毛片免费看| 特大巨黑吊av在线直播| 老司机影院毛片| 成人特级av手机在线观看| 亚洲精品中文字幕在线视频 | 内地一区二区视频在线| 成人亚洲精品一区在线观看 | 欧美日韩国产mv在线观看视频 | 九色成人免费人妻av| 超碰97精品在线观看| 国产白丝娇喘喷水9色精品| 黑人高潮一二区| 国产亚洲一区二区精品| 五月开心婷婷网| xxx大片免费视频| 啦啦啦啦在线视频资源| 人体艺术视频欧美日本| 欧美国产精品一级二级三级 | 免费黄网站久久成人精品| 欧美精品一区二区大全| 观看美女的网站| 伊人久久国产一区二区| 91久久精品电影网| 人人妻人人添人人爽欧美一区卜 | 国产黄片美女视频| 国产亚洲欧美精品永久| 国产一级毛片在线| 免费观看av网站的网址| 国产亚洲最大av| 国产大屁股一区二区在线视频| 亚洲精品视频女| 少妇裸体淫交视频免费看高清| 国产乱人偷精品视频| 日产精品乱码卡一卡2卡三| 日韩欧美精品免费久久| 亚洲av中文字字幕乱码综合| 欧美bdsm另类| av女优亚洲男人天堂| 亚洲最大成人中文| 亚洲精品乱久久久久久| 狂野欧美白嫩少妇大欣赏| 岛国毛片在线播放| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三区在线 | 美女视频免费永久观看网站| 91久久精品电影网| 久久97久久精品| 国产免费又黄又爽又色| 3wmmmm亚洲av在线观看| 久久久久国产精品人妻一区二区| 亚洲成人av在线免费| 成人18禁高潮啪啪吃奶动态图 | 国产女主播在线喷水免费视频网站| 人人妻人人添人人爽欧美一区卜 | 制服丝袜香蕉在线| 高清不卡的av网站| 联通29元200g的流量卡| 哪个播放器可以免费观看大片| av不卡在线播放| 亚洲欧美成人精品一区二区| 亚洲国产欧美在线一区| 最近中文字幕高清免费大全6| 日韩一区二区视频免费看| 男女下面进入的视频免费午夜| 男人爽女人下面视频在线观看| 岛国毛片在线播放| 一级毛片电影观看| 国产日韩欧美亚洲二区| 亚洲美女黄色视频免费看| av不卡在线播放| 国产精品人妻久久久久久| 一本色道久久久久久精品综合| 99热这里只有是精品在线观看| 一级毛片黄色毛片免费观看视频| 超碰av人人做人人爽久久| 精品久久久噜噜| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美人成| 亚洲图色成人| 人体艺术视频欧美日本| 夫妻午夜视频| 久久国产精品男人的天堂亚洲 | 一区二区三区精品91| 男人和女人高潮做爰伦理| 国产在线免费精品| 日韩强制内射视频| 亚洲人成网站高清观看| 精品久久久久久久久亚洲| 在线免费观看不下载黄p国产| 免费观看av网站的网址| 六月丁香七月| 91狼人影院| 国产男女内射视频| 国产深夜福利视频在线观看| 国产黄频视频在线观看| 亚洲精品日本国产第一区| .国产精品久久| 视频中文字幕在线观看| 人人妻人人看人人澡| 日韩三级伦理在线观看| 免费黄网站久久成人精品| 亚洲精品,欧美精品| 午夜老司机福利剧场| 久久久久国产网址| 成人18禁高潮啪啪吃奶动态图 | 1000部很黄的大片| 美女福利国产在线 | 22中文网久久字幕| 日本欧美国产在线视频| 黄色日韩在线| 国产综合精华液| 小蜜桃在线观看免费完整版高清| 黑人高潮一二区| 高清在线视频一区二区三区| 成人黄色视频免费在线看| 国产黄色视频一区二区在线观看| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品999| 亚洲国产欧美人成| 国产免费又黄又爽又色| 七月丁香在线播放| 内地一区二区视频在线| 成人黄色视频免费在线看| 超碰97精品在线观看| 蜜臀久久99精品久久宅男| 亚洲精品成人av观看孕妇| 国产高潮美女av| 午夜福利影视在线免费观看| 黄色一级大片看看| 国产av国产精品国产| 十八禁网站网址无遮挡 | 久久人人爽人人爽人人片va| 极品少妇高潮喷水抽搐| 成人国产麻豆网| 国产欧美日韩一区二区三区在线 | 国产精品伦人一区二区| 中文字幕av成人在线电影| 亚洲国产色片| 亚洲欧美一区二区三区黑人 | 亚洲欧美精品自产自拍| 2021少妇久久久久久久久久久| 精品国产露脸久久av麻豆| 亚洲高清免费不卡视频| 国产成人freesex在线| 成人漫画全彩无遮挡| 少妇的逼水好多| 国产精品爽爽va在线观看网站| 99久久精品一区二区三区| 国产亚洲91精品色在线| av在线老鸭窝| 国产成人91sexporn| 亚洲精品乱码久久久v下载方式| 亚洲欧洲国产日韩| 国产爽快片一区二区三区| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 亚洲精品视频女| 久久久久国产精品人妻一区二区| 午夜福利高清视频| av女优亚洲男人天堂| 身体一侧抽搐| 色婷婷久久久亚洲欧美| 一级毛片黄色毛片免费观看视频| 午夜老司机福利剧场| 欧美日韩一区二区视频在线观看视频在线| 男男h啪啪无遮挡| h日本视频在线播放| 国产成人免费无遮挡视频| 精品酒店卫生间| 少妇猛男粗大的猛烈进出视频| 免费大片黄手机在线观看| 观看美女的网站| 亚洲久久久国产精品| 亚洲色图综合在线观看| 日本黄色日本黄色录像| 亚洲精品久久午夜乱码| 亚洲久久久国产精品| 综合色丁香网| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 久久久久国产网址| 一二三四中文在线观看免费高清| 大香蕉久久网| videossex国产| 国产熟女欧美一区二区| 春色校园在线视频观看| 亚洲电影在线观看av| 日韩亚洲欧美综合| 五月开心婷婷网| 国产精品久久久久久av不卡| 三级经典国产精品| 美女内射精品一级片tv| 九九在线视频观看精品| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区黑人 | 一区二区三区精品91| 91精品国产九色| 男女下面进入的视频免费午夜| 内地一区二区视频在线| 99久久精品热视频| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 少妇的逼水好多| 久久人人爽av亚洲精品天堂 | 国产高清国产精品国产三级 | 免费av中文字幕在线| 国产伦在线观看视频一区| 一级毛片久久久久久久久女| 免费看不卡的av| 黄色怎么调成土黄色| 看十八女毛片水多多多| 插阴视频在线观看视频| 亚洲av日韩在线播放| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 亚洲国产成人一精品久久久| 99热国产这里只有精品6| 中文字幕亚洲精品专区| 精品久久久久久电影网| 亚洲国产毛片av蜜桃av| 麻豆国产97在线/欧美|