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

    壓電材料全矩陣材料常數(shù)超聲諧振譜反演技術(shù)中的變溫模式識別?

    2017-08-01 00:35:54湯立國
    物理學(xué)報 2017年2期
    關(guān)鍵詞:變溫模式識別壓電

    湯立國

    (廈門大學(xué),水聲通信與海洋信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,廈門 361102)

    壓電材料全矩陣材料常數(shù)超聲諧振譜反演技術(shù)中的變溫模式識別?

    湯立國?

    (廈門大學(xué),水聲通信與海洋信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,廈門 361102)

    (2016年9月20日收到;2016年10月18日收到修改稿)

    利用傳統(tǒng)的超聲脈沖-回波與電諧振技術(shù)定征壓電材料全矩陣材料參數(shù),必須采用多塊尺寸差異顯著的樣品,故很可能導(dǎo)致定征結(jié)果不自洽.超聲諧振譜(RUS)技術(shù)僅需一塊樣品即可對壓電材料全矩陣材料參數(shù)進(jìn)行定征,故可確保定征結(jié)果的自洽.由于實(shí)際測量諧振譜中模式混疊與遺漏現(xiàn)象不可避免,使得諧振譜中諧振模式的準(zhǔn)確識別成為RUS技術(shù)順利實(shí)施的最大難點(diǎn).本文提出一種諧振模式的變溫識別技術(shù).溫度變化可導(dǎo)致壓電體材料參數(shù)發(fā)生變化,材料參數(shù)的改變可影響各諧振模式的振動頻率,且對不同諧振模式影響不一致,因此改變測量環(huán)境溫度,有可能使得所測量超聲諧振譜中某些原本混疊的模式分開或使得某些原本遺漏的模式出現(xiàn).壓電陶瓷(PZT-8)的實(shí)驗(yàn)結(jié)果表明,該技術(shù)可有效提高諧振譜中諧振模式識別準(zhǔn)確率,從而保證了RUS反演的可靠性.

    壓電材料,材料常數(shù),超聲諧振譜,變溫模式識別

    1 引 言

    壓電材料無論在民用還是軍用領(lǐng)域,皆有著廣泛應(yīng)用[1,2],如工業(yè)檢測超聲與醫(yī)療超聲設(shè)備中的超聲壓電換能器及水聲工程中水聲壓電換能器的核心材料皆為壓電材料.壓電材料在制備完成之后,首先必須對其材料參數(shù)進(jìn)行定征,方能用于換能器的設(shè)計(jì).若利用錯誤的材料參數(shù)進(jìn)行壓電換能器的模擬設(shè)計(jì),很可能導(dǎo)致錯誤的結(jié)果,因此對壓電材料全矩陣材料參數(shù)進(jìn)行精確定征尤為重要.

    目前,定征壓電材料全矩陣材料參數(shù)最普遍的方法是超聲脈沖-回波與電諧振技術(shù)[3-5].利用上述方法定征壓電材料的全矩陣材料參數(shù),需多塊尺寸差異顯著的樣品,從不同樣品可分別定征出若干個材料參數(shù),如從d31振子可得恒電場順度系數(shù)橫向機(jī)電耦合系數(shù)k31,壓電應(yīng)變系數(shù)d31及自由介電常數(shù)而從d33振子可得恒電位移順度系數(shù)橫向機(jī)電耦合系數(shù)k33,恒電場順度系數(shù)及壓電應(yīng)變系數(shù)d33,最后將從不同樣品得到的參數(shù)進(jìn)行綜合可得全矩陣材料參數(shù)[6].尺寸差異顯著的不同壓電振子極化程度往往不一樣,如薄片狀的d31振子顯然比桿狀的d33振子更易極化,而極化程度不同將導(dǎo)致樣品材料參數(shù)存在差異.因此,利用多個樣品定征出的壓電材料全矩陣材料參數(shù)很有可能不自洽[7-9].

    欲保證所定征全矩陣材料參數(shù)自洽,采用樣品數(shù)目必須盡可能少,最理想定征是所有參數(shù)來自于同一塊樣品.迄今為止,能夠?qū)崿F(xiàn)僅用一塊樣品即可定征壓電材料全矩陣材料參數(shù)的方法有兩種:超聲諧振譜(RUS)技術(shù)[9-13]與電阻抗譜技術(shù)[14].與RUS技術(shù)相比,利用電阻抗譜技術(shù)定征材料參數(shù)效率較低,因?yàn)樵诜囱莸^程中需不斷調(diào)用有限元軟件計(jì)算樣品電阻抗譜,而計(jì)算頻率分辨率足夠高的電阻抗譜非常耗時.RUS技術(shù)的基本原理如下:利用兩個超聲換能器頂住長方體樣品的兩個對角或圓柱體樣品的邊沿,掃頻信號經(jīng)功率放大器放大后驅(qū)動發(fā)射換能器產(chǎn)生振動,接收換能器所接收樣品振動時域信號經(jīng)傅里葉變換后即可得到樣品的超聲諧振譜,對該譜進(jìn)行模式識別可得到樣品的若干諧振頻率;樣品的諧振頻率決定于樣品的幾何及材料參數(shù),其中幾何參數(shù)可直接測量得到,反之,由所測若干諧振頻率可對樣品的材料參數(shù)進(jìn)行反演.RUS技術(shù)最初由Frazer與LeCraw[15]引入材料學(xué)領(lǐng)域進(jìn)行固體彈性常數(shù)的反演.該技術(shù)曾被成功應(yīng)用于月球巖石樣品材料參數(shù)的精確反演.Ohno[10]最先將RUS技術(shù)引入壓電體材料參數(shù)的反演.Ogi等[16,17]與Nakamura等[18]利用RUS技術(shù)對鈮酸鋰(LiNbO3),硅酸鎵鑭(La3Ga5SiO14),α-石英(α-SiO2)等機(jī)械品質(zhì)因子QM值非常高的壓電材料的彈性及壓電常數(shù)進(jìn)行了反演.Tang與Cao[9,19]利用變溫RUS技術(shù)獲得了自洽的PZT-4壓電陶瓷及錳摻雜0.24Pb(In1/2Nb1/2)O3-0.46Pb(Mg1/3Nb2/3)O3-0.30PbTiO3單晶隨溫度變化的全矩陣材料參數(shù).

    RUS技術(shù)實(shí)施的最大難點(diǎn)在于對所測超聲諧振譜進(jìn)行分析并識別出足夠多諧振模式,此為樣品材料參數(shù)正確反演的基礎(chǔ).一般,材料獨(dú)立參數(shù)越多,模式識別與反演難度越高.然而,在超聲諧振譜測量過程中,諧振模式的遺漏與混疊幾乎不可避免,其對模式識別造成極大干擾.在RUS技術(shù)實(shí)施過程中,允許模式遺漏與模式混疊現(xiàn)象的發(fā)生,進(jìn)行參數(shù)反演時,舍棄這些模式即可,但是在進(jìn)行參數(shù)反演時,各模式的次序必須明確,因此必須明確何處發(fā)生了模式遺漏及模式混疊,否則將無法對諧振譜中各模式的次序作準(zhǔn)確識別.

    本文提出一種諧振模式的變溫識別方法.溫度變化可導(dǎo)致樣品材料參數(shù)發(fā)生變化,而材料參數(shù)的變化往往導(dǎo)致各模式諧振頻率發(fā)生改變,然而各材料參數(shù)的變化對不同諧振模式的影響并不一致,因此在溫度T1下混疊的模式在溫度T2下可能分開,原本遺漏的模式也有可能出現(xiàn).諧振模式的變溫識別方法正是基于此.壓電陶瓷(PZT-8)樣品的變溫超聲諧振譜測試實(shí)驗(yàn)結(jié)果表明,通過變溫可有效提高模式識別效率,從而確保了后續(xù)材料參數(shù)的正確反演.

    2 壓電體本征振動理論簡介

    壓電體的Lagrangian量可表示為

    其中,Sij是應(yīng)變張量分量,是恒電場彈性剛度系數(shù),φ是電勢,是自由介電常數(shù),emkl是壓電系數(shù),ρ是密度,ω是角頻率,ui是位移分量.位移ui與電勢φ可展開為

    其中,{vp}與{Ψr}為正交基函數(shù),本文數(shù)值計(jì)算中采用Legendre函數(shù)作為正交基函數(shù),與br為展開系數(shù).利用Rayleigh-Ritz方法,可得壓電體本征振動方程:

    3 模式識別障礙與變溫模式識別

    對于PZT-8長方體樣品,其振動模式根據(jù)對稱特征可分為Ag,Bg,Au及Bu四大類[20].圖1為PZT-8長方體樣品諧振模式Ag-8-Ag-10,Au-8-Au-10,Bg-8-Bg-10,及Bu-8-Bu-10的頻率隨恒電場彈性常數(shù)與及壓電應(yīng)力常數(shù)e15與e33的變化關(guān)系.樣品尺寸為4.8 mm×5.4 mm×5.0 mm,密度為7600 kg/m3.進(jìn)行數(shù)值求解時,改變某個參數(shù)的值,其他材料參數(shù)如表1所示.

    圖1 材料常數(shù)變化對諧振頻率的影響Fig.1.Resonance frequencies as functions of elastic and piezoelectric constants:(a)

    表1 PZT-8材料常數(shù)Table 1.Material constants of PZT-8.

    3.1 模式識別障礙分析

    在RUS技術(shù)實(shí)施過程中,為保證反演結(jié)果的正確,首先必須對所測超聲諧振譜中諧振模式進(jìn)行準(zhǔn)確識別,此為RUS技術(shù)實(shí)施的最大障礙與難點(diǎn).一般,材料獨(dú)立參數(shù)越多,模式識別與反演難度越高,因?yàn)楠?dú)立參數(shù)越多,所需反演模式數(shù)目越多.由圖1(a)可知,當(dāng)在圖示的范圍內(nèi)變化時,Ag-9與Au-9模式,Bg-10與Au-10模式始終非常接近,當(dāng)約為13.6×1010N/m2時,Bu-9,Bu-10,Au-8模式非常接近.不同諧振模式頻率非常接近的現(xiàn)象由圖1(b)-圖1(d)亦可觀察到.進(jìn)行超聲諧振譜測量時,上述相鄰且非常接近的模式易混疊在一起而難以識別.此外,實(shí)際測量過程中,有些模式振動非常弱,接收信號的信噪比很低,在諧振譜中無法觀察到該模式的諧振峰,因此雖然該模式未與相鄰模式發(fā)生混疊,但依然無法識別.利用RUS技術(shù)進(jìn)行材料參數(shù)反演時,允許模式的遺漏,但是所識別各模式的次序必須明確,即必須明確何處發(fā)生了模式混疊及遺漏,否則將無法對諧振譜中所識別各諧振頻率進(jìn)行準(zhǔn)確排序.

    3.2 變溫模式識別

    一般,溫度變化可導(dǎo)致壓電樣品材料參數(shù)發(fā)生變化,而材料參數(shù)的變化往往導(dǎo)致各諧振模式頻率發(fā)生改變,且各材料參數(shù)的變化對不同諧振模式的影響并不一致.例如,由圖1(a)可知,與模式Bg-8相比,模式Ag-10對的變化更加敏感;由圖1(c)可知,與模式Au-9或Ag-9相比,模式Au-10對e15的變化更加敏感.正是由于材料參數(shù)變化對不同諧振模式的影響不一致,使得超聲諧振譜的變溫模式識別成為可能,因?yàn)闇囟萒1下所測諧振譜中混疊的模式,在溫度T2下所測諧振譜中有可能分開而變得易于識別.此外,溫度改變導(dǎo)致材料常數(shù)發(fā)生改變,亦有可能導(dǎo)致實(shí)際測量中,不同模式的振動強(qiáng)度發(fā)生改變,從而使得原先振動很弱以至于在諧振譜中無法識別的模式,振動得到加強(qiáng)而易于識別.

    超聲諧振譜測試系統(tǒng)如文獻(xiàn)[9]中圖1所示.系統(tǒng)采用DRS Q9000超聲諧振譜測試儀測試樣品的諧振譜.圖2為PZT-8長方體樣品在20與80°C環(huán)境下所測得超聲諧振譜圖.樣品尺寸為5.195 mm×4.818 mm×5.318 mm,密度為7550 kg/m3.由圖2(a)可知,當(dāng)溫度T=20°C時,諧振譜中Bg-4與Ag-6模式諧振頻率非常接近,彼此混疊而易被誤判成僅有一個模式,當(dāng)T=80°C時,材料參數(shù)發(fā)生了改變,且材料參數(shù)變化對該兩個模式影響不一致,此二模式逐漸分離,而變得易于識別.由圖2(b)可知,當(dāng)溫度T=20°C時,諧振譜中模式Au-9雖然未與相鄰模式發(fā)生混疊,但是由于振動非常弱而無法識別,但當(dāng)溫度T=80°C時,模式Au-9的振動得到了加強(qiáng)而非常容易識別,此外,隨著溫度升高,Bg-8模式振動亦得到了加強(qiáng).由圖2(c)可知,隨著溫度升高,Bg-11模式振動得到了加強(qiáng),諧振譜中該模式的諧振峰越顯尖銳,而越易識別.

    圖2 PZT-8長方體樣品在20與80?C環(huán)境下所測得超聲諧振譜圖 (a)290-350 kHz;(b)420-470 kHz;(c)530-570 kHzFig.2.Resonant ultrasound spectra of a PZT-8 sample measured at 20 and 80?C:(a)290-350 kHz;(b)420-470 kHz;(c)530-570 kHz.

    4 PZT-8樣品諧振譜模式識別與參數(shù)反演

    由于PZT-8樣品諧振模式頻率對介電常數(shù)的變化非常不敏感,如圖3所示,故介電常數(shù)無法利用RUS技術(shù)反演,但是,可利用阻抗分析儀測得的電容直接計(jì)算得出.本文所采用PZT-8樣品的夾持介電常數(shù)測量結(jié)果為3.2小節(jié)中所述樣品的諧振模式識別結(jié)果如表2所列,其中,fmeas為從所測超聲諧振譜中識別出的各諧振模式頻率.表3為彈性剛度系數(shù)及壓電應(yīng)力系數(shù)反演結(jié)果.根據(jù)反演結(jié)果所計(jì)算恒電位移彈性剛度系數(shù)和的值分別為15.02×1010及3.850×1010N/m2,利用超聲脈沖-回波技術(shù)測量所得和的值分別為15.08×1010及3.832×1010N/m2,相對誤差皆小于0.5%,表明了反演結(jié)果的自洽.表2中fcal為根據(jù)反演結(jié)果所計(jì)算各諧振模式頻率,其與測量結(jié)果非常一致,相對誤差大多小于0.3%.

    圖3 介電常數(shù)變化對諧振頻率的影響Fig.3.Resonance frequencies as functions of dielectric constants:

    表2 實(shí)測及計(jì)算所得諧振頻率Table 2.Measured and calculated resonance frequencies.

    表2 實(shí)測及計(jì)算所得諧振頻率(續(xù))Table 2.Measured and calculated resonance frequencies(continued).

    若相鄰模式混疊在一起,則難以從諧振譜中將各自的諧振頻率準(zhǔn)確讀出,故反演時,舍棄了如Bg-4與Ag-6之類混疊在一起的模式.若將圖2(a)中混疊在一起的Bg-4與Ag-6兩個模式誤判為一個模式,則導(dǎo)致后續(xù)模式的排序錯誤,表4為該錯誤排序下的反演結(jié)果.雖然僅僅排序時漏了一個模式,但是e15的反演結(jié)果與正確結(jié)果相對誤差達(dá)到了38%,e33的反演結(jié)果與正確結(jié)果的相對誤差達(dá)到了27%.總而言之,識別模式排序的錯誤將直接導(dǎo)致反演的失敗,而諧振模式的變溫識別技術(shù)將有效提高模式的識別準(zhǔn)確率.

    表3 PZT-8材料常數(shù)反演結(jié)果Table 3.Inversion results of PZT-8.

    表4 排序有誤時的反演結(jié)果Table 4.Inversion results from the resonance frequencies with wrong order.

    5 結(jié) 論

    與傳統(tǒng)的超聲脈沖-回波與電諧振技術(shù)相比,RUS技術(shù)最大優(yōu)點(diǎn)在于僅需一塊樣品即可定征壓電材料全矩陣材料參數(shù),從而確保了定征結(jié)果的自洽.然而,實(shí)際測量所得超聲諧振譜中諧振模式混疊與遺漏現(xiàn)象不可避免,使得諧振模式的準(zhǔn)確識別成為RUS技術(shù)順利實(shí)施的最大難點(diǎn),且材料獨(dú)立參數(shù)越多,模式識別與反演難度越高.提高模式識別效率與準(zhǔn)確率一直是該領(lǐng)域研究重點(diǎn)之一.本文提出一種諧振模式的變溫識別技術(shù).PZT-8樣品的實(shí)驗(yàn)結(jié)果表明,該技術(shù)可有效提高諧振譜中諧振模式識別準(zhǔn)確率.需指出,變溫識別技術(shù)并不能徹底解決模式識別問題,將不同的識別技術(shù)結(jié)合起來,如多點(diǎn)測振與變溫識別,以提高模式識別的準(zhǔn)確度將是進(jìn)一步需要展開的研究.

    感謝美國賓夕法尼亞州立大學(xué)曹文武教授為本文測試及分析所提供的幫助.

    [1]Muralt P 2000J.Micromech.Microeng.10 136

    [2]Zhou Q F,Lam K H,Zheng H R,Qiu W B,Shung K K 2014Prog.Mater.Sci.66 87

    [3]Zhang S J,Li F 2012J.Appl.Phys.111 031301

    [4]Zhang S J,Lee S M,Kim D H,Lee H Y,Shrout T R 2008J.Am.Ceram.Soc.91 683

    [5]Sun E W,Zhang R,Wu F M,Cao W W 2013J.Alloys Compd.553 267

    [6]Sun E W 2011Ph.D.Dissertation(Harbin:Harbin Institute of Technology)(in Chinese)[孫恩偉 2011博士學(xué)位論文(哈爾濱:哈爾濱工業(yè)大學(xué))]

    [7]Topolov V Y 2010Appl.Phys.Lett.96 196101

    [8]Topolov V Y,Bowen C R 2011J.Appl.Phys.109 094107

    [9]Tang L G,Cao W W 2015Appl.Phys.Lett.106 052902

    [10]Ohno I 1990Phys.Chem.Miner.17 371

    [11]Leisure R G,Willis F A 1997J.Phys.:Condens.Matter9 6001

    [12]Zadler B J,Rousseau J H L,Scales J A,Smith M L 2004Geophys.J.Int.156 154

    [13]Migliori A,Maynard J D 2005Rev.Sci.Instrum.76 121301

    [14]Li S Y,Zheng L M,Jiang W H,Sahul R,Gopalan V,Cao W W 2013J.Appl.Phys.114 104505

    [15]Frazer D B,LeCraw R C 1964Rev.Sci.Instrum.35 1113

    [16]Ogi H,Kawasaki Y,Hirao M,Ledbetter H 2002J.Appl.Phys.92 2451

    [17]Ogi H,Ohmori T,Nakamura N,Hirao M 2006J.Appl.Phys.100 053511

    [18]Nakamura N,Ogi H,Hirao M 2012J.Appl.Phys.111 013509

    [19]Tang L G,Tian H,Zhang Y,Cao W W 2016Appl.Phys.Lett.108 082901

    [20]Mochizuki E 1987J.Phys.Earth35 159

    PACS:77.84.-s,43.35.+d,43.35.Yb DOI:10.7498/aps.66.027703

    Mode identification via temperature variation in resonant ultrasonic spectroscopy technique for piezoelectric material?

    Tang Li-Guo?
    (Key Laboratory of Underwater Acoustic Communication and Marine Information Technology,Ministry of Education,Xiamen University,Xiamen 361102,China)

    20 September 2016;revised manuscript

    18 October 2016)

    The full matrix material constants of piezoelectric materials should be characterized first before they have been used to make actuators or sensors.Up to now,they are usually determined by the ultrasonic pulse-echo and electric impedance resonance techniques through using multiple samples with drastically different sizes.However,the constants determined by the aforementioned techniques are probably inconsistent because the sample-to-sample variation cannot be eliminated.The technique of resonant ultrasonic spectroscopy(RUS)only needs one sample to determine the full matrix constants of piezoelectric material.Therefore,the consistency of the constants is guaranteed.During the implementation of the RUS technique,the elastic stiffnessand piezoelectric constantseijcan be determined from the resonance modes identified from the resonant ultrasonic spectrum.The free and clamped dielectric constants cannot be determined by the RUS technique because they have very weak influence on resonance frequency.However,they can be directly measured from the same sample by using an impedance analyzer.To ensure the reliable inversion of material constants,enough resonance modes should be identified from the measured resonant ultrasonic spectrum.However,there are many missing and overlapped modes in the spectrum,which makes mode identification become a biggest obstacle to the implementation of the RUS technique.The adjacent modes may overlap if the resonance frequencies corresponding to them have a very small difference.In addition,the lower the mechanical quality factorQM,the more likely to overlap the adjacent modes are.During the RUS measurement,the rectangular parallelepiped sample is placed between the transmitting and receiving transducers with contacts only at the opposite corners of the sample.Resonance modes would not be detected if the receiving point,i.e.,one corner of the sample,is the node of these modes.Therefore,there are missing modes in the resonant ultrasonic spectrum.To overcome the difficulty in identifying the modes,caused by modes missing and overlapping,the mode identifying method via temperature variation is presented in this study.Note that a change of temperature may change the material properties of a piezoelectric sample.The material properties have a great influence on the resonance frequency of the sample.Moreover,the influences corresponding to resonance modes are different.Therefore,the variation of temperature may make the overlapped modes separated from each other and the missing modes appear,namely,the missing and overlapped modes may be identified by comparing the resonant ultrasonic spectra measured at different temperatures.The experimental results of piezoelectric ceramics(PZT-8)show that this method can effectively improve the accuracy of mode identification and guarantee the reliability of inversion in the RUS technique.

    piezoelectric materials,material constants,resonant ultrasonic spectroscopy,mode identi fication via temperature variation

    :77.84.-s,43.35.+d,43.35.Yb

    10.7498/aps.66.027703

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

    ?通信作者.E-mail:liguotang@xmu.edu.cn

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

    ?Corresponding author.E-mail:liguotang@xmu.edu.cn

    猜你喜歡
    變溫模式識別壓電
    《壓電與聲光》征稿啟事
    壓電與聲光(2019年1期)2019-02-22 09:46:06
    氯乙烯生產(chǎn)中變溫吸附脫水工藝的使用及改進(jìn)
    新型壓電疊堆泵設(shè)計(jì)及仿真
    淺談模式識別在圖像識別中的應(yīng)用
    電子測試(2017年23期)2017-04-04 05:06:50
    第四屆亞洲模式識別會議
    凍融處理對甘薯變溫壓差膨化干燥動力學(xué)的影響
    非共面四頻激光陀螺變溫零偏周期性波動
    Mn摻雜ZnSe量子點(diǎn)變溫發(fā)光性質(zhì)研究
    第3屆亞洲模式識別會議
    基于壓電激振的彈性模量測量方法
    高清午夜精品一区二区三区| 在线观看一区二区三区激情| 亚洲精品久久成人aⅴ小说 | 七月丁香在线播放| 黄色怎么调成土黄色| 久久久久久久精品精品| 一区二区三区免费毛片| 丁香六月天网| 在线精品无人区一区二区三| 精品久久久久久久久亚洲| 久久影院123| 中文字幕av电影在线播放| 在线观看www视频免费| av视频免费观看在线观看| 在线精品无人区一区二区三| 26uuu在线亚洲综合色| 精品久久国产蜜桃| 99热6这里只有精品| 一区二区日韩欧美中文字幕 | 国产精品国产三级专区第一集| 极品人妻少妇av视频| 青春草亚洲视频在线观看| 亚洲欧美一区二区三区国产| 欧美老熟妇乱子伦牲交| 日本av免费视频播放| 中文字幕最新亚洲高清| 亚洲欧美日韩另类电影网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜福利视频在线观看免费| 你懂的网址亚洲精品在线观看| 日本wwww免费看| 少妇猛男粗大的猛烈进出视频| 美女视频免费永久观看网站| 国产免费现黄频在线看| 成人国语在线视频| 亚洲国产色片| 国产又色又爽无遮挡免| 欧美日韩在线观看h| 国产高清国产精品国产三级| 欧美亚洲日本最大视频资源| 久久精品人人爽人人爽视色| 少妇人妻 视频| 免费观看无遮挡的男女| av在线app专区| 国产成人精品在线电影| 亚洲综合色惰| 精品一区二区三区视频在线| 久热这里只有精品99| av网站免费在线观看视频| 久久人人爽人人爽人人片va| 免费黄色在线免费观看| 99re6热这里在线精品视频| 街头女战士在线观看网站| 亚洲国产精品一区三区| 性色av一级| 久久av网站| 久久99精品国语久久久| 日韩成人av中文字幕在线观看| 久久精品久久久久久久性| 国产精品不卡视频一区二区| 欧美激情 高清一区二区三区| 交换朋友夫妻互换小说| 夫妻午夜视频| 国产无遮挡羞羞视频在线观看| 99久久精品一区二区三区| 母亲3免费完整高清在线观看 | 美女大奶头黄色视频| 国产亚洲欧美精品永久| 五月天丁香电影| 日韩视频在线欧美| av在线观看视频网站免费| 成人漫画全彩无遮挡| 青春草国产在线视频| 91精品一卡2卡3卡4卡| 一级,二级,三级黄色视频| a级毛片免费高清观看在线播放| 欧美xxxx性猛交bbbb| 少妇被粗大的猛进出69影院 | 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 国产免费又黄又爽又色| 国产极品天堂在线| 少妇猛男粗大的猛烈进出视频| 久久精品国产亚洲av天美| 国产成人精品在线电影| 日本猛色少妇xxxxx猛交久久| 曰老女人黄片| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 91在线精品国自产拍蜜月| 精品99又大又爽又粗少妇毛片| av有码第一页| 亚洲综合色网址| 熟女人妻精品中文字幕| 九草在线视频观看| 大码成人一级视频| 丝袜脚勾引网站| 少妇被粗大的猛进出69影院 | 七月丁香在线播放| 少妇人妻精品综合一区二区| 久久久久国产网址| 亚洲国产av新网站| 熟女人妻精品中文字幕| 黄色配什么色好看| 成人国产麻豆网| 一区二区三区乱码不卡18| 亚洲av成人精品一区久久| 成人毛片60女人毛片免费| 黑人猛操日本美女一级片| 麻豆乱淫一区二区| 国产免费现黄频在线看| 午夜久久久在线观看| 秋霞在线观看毛片| 欧美最新免费一区二区三区| 午夜久久久在线观看| 麻豆乱淫一区二区| 大陆偷拍与自拍| 777米奇影视久久| 亚洲欧美一区二区三区国产| 五月开心婷婷网| 久久久久久伊人网av| 99久久精品国产国产毛片| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| 黄色怎么调成土黄色| 亚洲av在线观看美女高潮| 欧美三级亚洲精品| 9色porny在线观看| 黑人巨大精品欧美一区二区蜜桃 | 成人国产av品久久久| 久久免费观看电影| 丝袜美足系列| 一区二区日韩欧美中文字幕 | 女人精品久久久久毛片| 国产日韩欧美视频二区| 夫妻性生交免费视频一级片| 男女国产视频网站| 中文字幕亚洲精品专区| 人妻夜夜爽99麻豆av| 亚洲欧洲精品一区二区精品久久久 | 99九九在线精品视频| tube8黄色片| 亚洲av欧美aⅴ国产| 在线观看国产h片| 在线亚洲精品国产二区图片欧美 | 热99国产精品久久久久久7| 老熟女久久久| 精品一品国产午夜福利视频| 久久97久久精品| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区国产| av线在线观看网站| 国产免费福利视频在线观看| xxxhd国产人妻xxx| 精品一区二区三卡| 欧美日韩精品成人综合77777| av卡一久久| 在线观看一区二区三区激情| 哪个播放器可以免费观看大片| 午夜老司机福利剧场| av专区在线播放| 日本av免费视频播放| 最后的刺客免费高清国语| 成人午夜精彩视频在线观看| 精品久久久久久电影网| 久久免费观看电影| 最后的刺客免费高清国语| av专区在线播放| 观看av在线不卡| 午夜福利影视在线免费观看| 一个人看视频在线观看www免费| 国产成人免费观看mmmm| 水蜜桃什么品种好| 色网站视频免费| 久久久精品94久久精品| 中文字幕人妻熟人妻熟丝袜美| 一级,二级,三级黄色视频| 久久 成人 亚洲| 中文字幕久久专区| 成人黄色视频免费在线看| 搡老乐熟女国产| 免费日韩欧美在线观看| 亚洲少妇的诱惑av| 欧美三级亚洲精品| 中文字幕av电影在线播放| 亚洲中文av在线| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 国产精品一区二区在线观看99| 亚洲内射少妇av| 亚洲成人一二三区av| 插逼视频在线观看| 日本爱情动作片www.在线观看| 国产色爽女视频免费观看| 久久99精品国语久久久| 亚洲国产色片| 美女脱内裤让男人舔精品视频| 免费av中文字幕在线| 免费黄网站久久成人精品| 最近中文字幕2019免费版| 日本wwww免费看| 国产日韩一区二区三区精品不卡 | 男人爽女人下面视频在线观看| 午夜福利视频在线观看免费| 日本午夜av视频| 欧美 亚洲 国产 日韩一| 午夜福利,免费看| 我的女老师完整版在线观看| 性色avwww在线观看| 九九爱精品视频在线观看| 国产午夜精品一二区理论片| 美女大奶头黄色视频| 精品午夜福利在线看| 日韩 亚洲 欧美在线| 精品酒店卫生间| 一边亲一边摸免费视频| 黑人巨大精品欧美一区二区蜜桃 | 91久久精品国产一区二区三区| 99久久人妻综合| freevideosex欧美| 国产精品久久久久久精品电影小说| www.av在线官网国产| 九色成人免费人妻av| 女的被弄到高潮叫床怎么办| 国产精品偷伦视频观看了| 亚洲欧美清纯卡通| videosex国产| 久久青草综合色| 亚洲精品久久成人aⅴ小说 | av线在线观看网站| 观看av在线不卡| 黄色视频在线播放观看不卡| 精品熟女少妇av免费看| 视频中文字幕在线观看| 两个人免费观看高清视频| 亚洲综合色惰| 蜜桃久久精品国产亚洲av| 亚洲av中文av极速乱| 午夜福利影视在线免费观看| 久久久国产精品麻豆| 国产亚洲最大av| 又黄又爽又刺激的免费视频.| 国产精品成人在线| 少妇 在线观看| 美女脱内裤让男人舔精品视频| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 视频区图区小说| 亚洲精品成人av观看孕妇| 欧美日韩一区二区视频在线观看视频在线| 91在线精品国自产拍蜜月| 午夜视频国产福利| xxx大片免费视频| 18禁观看日本| 久久免费观看电影| 男人爽女人下面视频在线观看| 成人免费观看视频高清| 黄色怎么调成土黄色| 久热这里只有精品99| 七月丁香在线播放| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人看| 三级国产精品片| 午夜福利在线观看免费完整高清在| 中文字幕av电影在线播放| kizo精华| 国产片内射在线| 日本91视频免费播放| 日本午夜av视频| 亚洲五月色婷婷综合| 欧美xxⅹ黑人| 高清欧美精品videossex| 免费大片黄手机在线观看| 免费av中文字幕在线| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 麻豆精品久久久久久蜜桃| 国产女主播在线喷水免费视频网站| 国产精品久久久久久av不卡| 国产精品一区二区三区四区免费观看| 久久午夜综合久久蜜桃| 国产精品国产av在线观看| 成人国产av品久久久| 久久精品国产自在天天线| 秋霞在线观看毛片| 亚洲av国产av综合av卡| 美女国产高潮福利片在线看| 在线看a的网站| 亚洲高清免费不卡视频| 成人亚洲精品一区在线观看| 国精品久久久久久国模美| 在线播放无遮挡| 日本猛色少妇xxxxx猛交久久| 极品少妇高潮喷水抽搐| 国产伦精品一区二区三区视频9| 999精品在线视频| 亚洲激情五月婷婷啪啪| 51国产日韩欧美| 制服丝袜香蕉在线| 国内精品宾馆在线| 你懂的网址亚洲精品在线观看| 中文字幕av电影在线播放| 亚洲无线观看免费| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 日韩av免费高清视频| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 精品亚洲乱码少妇综合久久| 国产成人91sexporn| 国产精品蜜桃在线观看| 免费观看a级毛片全部| 麻豆乱淫一区二区| av在线观看视频网站免费| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠久久av| 色婷婷久久久亚洲欧美| 尾随美女入室| 女人精品久久久久毛片| 最近2019中文字幕mv第一页| 两个人免费观看高清视频| 久久热精品热| 91久久精品国产一区二区三区| 亚洲三级黄色毛片| 国产亚洲欧美精品永久| 亚洲精品久久成人aⅴ小说 | 精品亚洲成a人片在线观看| 免费播放大片免费观看视频在线观看| 美女视频免费永久观看网站| 五月伊人婷婷丁香| 欧美日韩亚洲高清精品| 国产欧美日韩综合在线一区二区| 三上悠亚av全集在线观看| 国产精品一国产av| 亚洲精品中文字幕在线视频| 免费人妻精品一区二区三区视频| 99热这里只有是精品在线观看| 亚洲精品自拍成人| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 韩国高清视频一区二区三区| 国产精品一国产av| 蜜桃国产av成人99| 69精品国产乱码久久久| 久久久久久久大尺度免费视频| 91aial.com中文字幕在线观看| 99热全是精品| 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 亚洲精品日韩av片在线观看| 国内精品宾馆在线| 肉色欧美久久久久久久蜜桃| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 国产 精品1| 人成视频在线观看免费观看| 亚洲欧美一区二区三区国产| 日产精品乱码卡一卡2卡三| 久久精品夜色国产| av在线观看视频网站免费| 最近的中文字幕免费完整| 国产永久视频网站| 久久精品夜色国产| 国产精品成人在线| 欧美 日韩 精品 国产| 五月伊人婷婷丁香| 自线自在国产av| 搡女人真爽免费视频火全软件| 熟妇人妻不卡中文字幕| 亚洲国产av影院在线观看| 老熟女久久久| 爱豆传媒免费全集在线观看| 精品视频人人做人人爽| 国产视频首页在线观看| a级片在线免费高清观看视频| 91久久精品国产一区二区三区| 赤兔流量卡办理| 黑人猛操日本美女一级片| 考比视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 国产国拍精品亚洲av在线观看| 成人综合一区亚洲| 日本猛色少妇xxxxx猛交久久| 99久久精品一区二区三区| 黑人猛操日本美女一级片| 成人午夜精彩视频在线观看| 欧美性感艳星| 欧美日本中文国产一区发布| 热re99久久国产66热| 欧美最新免费一区二区三区| 2021少妇久久久久久久久久久| 天堂中文最新版在线下载| 亚洲精品日本国产第一区| 亚洲av欧美aⅴ国产| 亚洲人与动物交配视频| 亚洲精品,欧美精品| 亚洲成人一二三区av| 3wmmmm亚洲av在线观看| 伊人久久国产一区二区| 国产男人的电影天堂91| 国产又色又爽无遮挡免| 最近中文字幕高清免费大全6| 日韩,欧美,国产一区二区三区| av国产精品久久久久影院| 岛国毛片在线播放| 久久午夜福利片| 九九在线视频观看精品| 下体分泌物呈黄色| 91国产中文字幕| 蜜桃国产av成人99| 亚洲怡红院男人天堂| 国产亚洲av片在线观看秒播厂| 爱豆传媒免费全集在线观看| 两个人免费观看高清视频| 制服人妻中文乱码| 人妻制服诱惑在线中文字幕| 精品久久久精品久久久| 欧美精品一区二区大全| 国产亚洲一区二区精品| 99热国产这里只有精品6| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| h视频一区二区三区| 97在线视频观看| 亚洲无线观看免费| 国产成人a∨麻豆精品| 精品国产露脸久久av麻豆| 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 国产片特级美女逼逼视频| 亚洲激情五月婷婷啪啪| 如何舔出高潮| 女人精品久久久久毛片| 一区二区日韩欧美中文字幕 | av国产久精品久网站免费入址| 亚洲成人av在线免费| 免费日韩欧美在线观看| 国产极品天堂在线| 熟妇人妻不卡中文字幕| 在线观看三级黄色| 久久亚洲国产成人精品v| 国产av国产精品国产| 国产成人精品在线电影| 亚洲国产色片| 制服诱惑二区| 精品亚洲成a人片在线观看| 大片电影免费在线观看免费| 亚洲无线观看免费| 亚洲婷婷狠狠爱综合网| 亚洲精品aⅴ在线观看| 欧美3d第一页| 免费播放大片免费观看视频在线观看| 又黄又爽又刺激的免费视频.| 日本猛色少妇xxxxx猛交久久| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 国产免费一级a男人的天堂| 亚洲精品成人av观看孕妇| 哪个播放器可以免费观看大片| 国产在线一区二区三区精| 啦啦啦在线观看免费高清www| 国产成人午夜福利电影在线观看| 免费观看无遮挡的男女| 寂寞人妻少妇视频99o| 韩国av在线不卡| 五月伊人婷婷丁香| 亚洲av成人精品一区久久| 伦理电影免费视频| 国产精品一区二区三区四区免费观看| 欧美 亚洲 国产 日韩一| 视频在线观看一区二区三区| 18禁观看日本| 亚洲综合色网址| 欧美激情 高清一区二区三区| 啦啦啦视频在线资源免费观看| 在线天堂最新版资源| 简卡轻食公司| 亚洲精品日韩av片在线观看| 老女人水多毛片| 欧美少妇被猛烈插入视频| 99re6热这里在线精品视频| 成人无遮挡网站| 久久久久久久久久久久大奶| 嘟嘟电影网在线观看| av免费在线看不卡| 亚洲人成网站在线观看播放| 91精品伊人久久大香线蕉| 满18在线观看网站| 欧美亚洲日本最大视频资源| 一本色道久久久久久精品综合| 国产又色又爽无遮挡免| 嫩草影院入口| 日韩一本色道免费dvd| 国产一级毛片在线| 各种免费的搞黄视频| 热99久久久久精品小说推荐| 99久久综合免费| 两个人的视频大全免费| 欧美成人精品欧美一级黄| 久久久久视频综合| 人妻一区二区av| videossex国产| 精品人妻偷拍中文字幕| 欧美日韩亚洲高清精品| 一级片'在线观看视频| 亚洲在久久综合| 中文字幕精品免费在线观看视频 | 9色porny在线观看| 熟女电影av网| 国产精品一区www在线观看| 婷婷色综合大香蕉| 国产高清三级在线| 视频中文字幕在线观看| 午夜福利在线观看免费完整高清在| 国产精品一二三区在线看| 久久久久久伊人网av| 免费看不卡的av| 王馨瑶露胸无遮挡在线观看| 亚洲激情五月婷婷啪啪| 乱人伦中国视频| 欧美少妇被猛烈插入视频| 女性生殖器流出的白浆| 国产成人aa在线观看| 国产国语露脸激情在线看| 日韩精品免费视频一区二区三区 | 热99国产精品久久久久久7| 中文欧美无线码| 又黄又爽又刺激的免费视频.| 久久综合国产亚洲精品| 成年美女黄网站色视频大全免费 | 桃花免费在线播放| a级片在线免费高清观看视频| 国模一区二区三区四区视频| 免费高清在线观看视频在线观看| 涩涩av久久男人的天堂| 国产日韩一区二区三区精品不卡 | 少妇猛男粗大的猛烈进出视频| 91久久精品电影网| 自线自在国产av| 激情五月婷婷亚洲| 国产av精品麻豆| 99热全是精品| 亚洲一级一片aⅴ在线观看| 亚洲国产精品999| 一本大道久久a久久精品| 国产精品成人在线| 精品久久久噜噜| xxxhd国产人妻xxx| 永久免费av网站大全| 亚洲av不卡在线观看| 亚洲内射少妇av| 久久久久久久久久久免费av| 亚洲五月色婷婷综合| 99热网站在线观看| 伦精品一区二区三区| 老司机亚洲免费影院| 少妇被粗大猛烈的视频| 久久这里有精品视频免费| 99久久综合免费| 免费高清在线观看视频在线观看| tube8黄色片| 精品少妇久久久久久888优播| 国产免费一级a男人的天堂| 中文字幕亚洲精品专区| 欧美少妇被猛烈插入视频| 成人亚洲欧美一区二区av| 青春草国产在线视频| 精品久久蜜臀av无| 国产精品一区二区三区四区免费观看| 最近2019中文字幕mv第一页| 亚洲av男天堂| 国产精品99久久久久久久久| 插阴视频在线观看视频| 街头女战士在线观看网站| 国产色婷婷99| 热re99久久国产66热| 日产精品乱码卡一卡2卡三| 亚洲欧美色中文字幕在线| 高清不卡的av网站| 午夜日本视频在线| videos熟女内射| 久久久精品94久久精品| 国产有黄有色有爽视频| 婷婷色av中文字幕| 亚洲欧洲日产国产| 99视频精品全部免费 在线| 国产69精品久久久久777片| 欧美精品国产亚洲| 国产成人精品在线电影| 在线观看www视频免费| 久久人人爽人人爽人人片va| 久久女婷五月综合色啪小说| 欧美丝袜亚洲另类| 精品久久久久久久久av| 久久久久国产网址| 能在线免费看毛片的网站| 在线看a的网站| 一边摸一边做爽爽视频免费| 最黄视频免费看| 亚洲少妇的诱惑av| 2018国产大陆天天弄谢| 尾随美女入室| 黄色毛片三级朝国网站| 99国产综合亚洲精品| 女性被躁到高潮视频| 精品人妻偷拍中文字幕| 亚洲成人一二三区av| 草草在线视频免费看| 免费少妇av软件| 久久久久人妻精品一区果冻| 伊人久久精品亚洲午夜| 亚洲综合精品二区| 欧美日韩视频精品一区| 精品一品国产午夜福利视频| 一区二区三区乱码不卡18|