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

    超聲導(dǎo)波任意波形激勵(lì)技術(shù)研究

    2016-10-14 06:44:22鄭甜甜潘曉宇
    電子設(shè)計(jì)工程 2016年12期
    關(guān)鍵詞:有限元信號(hào)檢測(cè)

    鄭甜甜,潘曉宇

    (北京交通大學(xué)機(jī)械與電子控制工程學(xué)院,北京100044)

    超聲導(dǎo)波任意波形激勵(lì)技術(shù)研究

    鄭甜甜,潘曉宇

    (北京交通大學(xué)機(jī)械與電子控制工程學(xué)院,北京100044)

    高效準(zhǔn)確的無(wú)縫鋼軌溫度應(yīng)力檢測(cè)與探傷是確保高速鐵路安全運(yùn)營(yíng)的重要技術(shù)手段,近年來(lái),基于超聲導(dǎo)波的鋼軌應(yīng)力檢測(cè)技術(shù)受到越來(lái)越大的關(guān)注,精確控制鋼軌中導(dǎo)波的波形是該技術(shù)的關(guān)鍵,但目前得到的接收波形精確度卻不高。本文利用超聲導(dǎo)波換能器和時(shí)間反轉(zhuǎn)法,通過(guò)ANSYS有限元仿真分析,獲得了任意波形的激勵(lì)信號(hào),提高了接收波形與理想波形的相關(guān)度。仿真結(jié)果表明,由時(shí)間反轉(zhuǎn)法得到的接收波形精確度提高了約12%。

    超聲導(dǎo)波換能器;時(shí)間反轉(zhuǎn)法;波形;激勵(lì)信號(hào)

    在鐵路系統(tǒng)中,無(wú)縫線路(CWR——continuouslyweldedrails)已經(jīng)成為世界鐵路領(lǐng)域的主要發(fā)展潮流。已有研究表明:軌溫變化1℃,無(wú)縫鋼軌固定區(qū)內(nèi)縱向應(yīng)力變化18.82 kN,若軌溫變化50℃,則鋼軌內(nèi)應(yīng)力變化為941 kN[1]??梢?jiàn)無(wú)縫線路長(zhǎng)鋼軌所承受的溫度應(yīng)力要比普通鋼軌大得多。因?yàn)檫@種使鋼軌變形的溫度應(yīng)力會(huì)導(dǎo)致脹軌或斷軌的發(fā)生[2],嚴(yán)重威脅到了列車的運(yùn)行安全,因此準(zhǔn)確及時(shí)地掌握鋼軌內(nèi)部溫度應(yīng)力具有重要意義。

    目前,根據(jù)原理的不同,鋼軌溫度應(yīng)力的檢測(cè)方法可以分為應(yīng)變法,應(yīng)力法和能量釋放法,具體方法又主要分為8種。其中,觀測(cè)樁法、應(yīng)變計(jì)法、光纖光柵法、巴克豪森法受內(nèi)部結(jié)構(gòu)和外部因素影響較大,應(yīng)力檢測(cè)結(jié)果的準(zhǔn)確率較低,標(biāo)定軌長(zhǎng)法實(shí)時(shí)性較差,效率較低[3-7];橫向加力法操作繁瑣,使用不方便;X射線法測(cè)量范圍有限[8];超聲波雖然有很強(qiáng)的穿透能力,在傳播時(shí)有很好的方向性,但超聲波法測(cè)量結(jié)果反映的是超聲波傳播過(guò)程中經(jīng)過(guò)的結(jié)構(gòu)路徑內(nèi)應(yīng)力情況,不能反映結(jié)構(gòu)內(nèi)部整體的平均應(yīng)力[9-10]。由于超聲導(dǎo)波具有傳輸距離遠(yuǎn)、可以引起鋼軌截面全部質(zhì)點(diǎn)振動(dòng)等優(yōu)點(diǎn),基于超聲導(dǎo)波的鋼軌應(yīng)力狀態(tài)在線監(jiān)測(cè)技術(shù)在國(guó)內(nèi)外有了長(zhǎng)足發(fā)展[11-13]。

    應(yīng)用超聲導(dǎo)波的一個(gè)關(guān)鍵是要精確控制鋼軌中導(dǎo)波的波形,這樣才能精確定位和匹配探測(cè)波包,測(cè)量其形狀、速度等的變化,進(jìn)而感知溫度應(yīng)力、缺陷等鋼軌關(guān)鍵狀態(tài)。目前對(duì)于超聲導(dǎo)波的研究基本都采用將需要的波形直接作為激勵(lì)波形的方法來(lái)獲得接收波形,但最終的接收波形與需要的波形的相關(guān)度有待提高,這就影響了最終應(yīng)力檢測(cè)的精確度。本文提出利用時(shí)間反轉(zhuǎn)法來(lái)控制激勵(lì)波形的策略,借助ANSYS有限元仿真,成功獲得了任意波形激勵(lì),提高了接收波形的精確度。

    1 超聲導(dǎo)波檢測(cè)原理

    超聲導(dǎo)波是一種在桿、板、管等結(jié)構(gòu)中傳播的波,是波不斷與介質(zhì)的上邊界和下邊界發(fā)生反射、折射以及橫波和縱波之間的相互轉(zhuǎn)換所產(chǎn)生的波[14]。超聲導(dǎo)波有兩個(gè)重要特性:群速度和相速度。相速度cp指的是導(dǎo)波上固定的相位點(diǎn)在介質(zhì)中傳播的速度;群速度cg是多個(gè)相近頻率的波在同一個(gè)介質(zhì)中傳播時(shí)共同合成包絡(luò)線的傳播速度,是波包的能量傳播速度,也就是導(dǎo)波的傳播速度。根據(jù)超聲導(dǎo)波的聲彈性原理,鋼軌內(nèi)的應(yīng)力與超聲導(dǎo)波在鋼軌內(nèi)的傳播速度成正比關(guān)系,因此,監(jiān)測(cè)到超聲導(dǎo)波的速度就能監(jiān)測(cè)鋼軌應(yīng)力,而在導(dǎo)波的速度檢測(cè)過(guò)程中波形控制是關(guān)鍵。

    超聲導(dǎo)波一般需要壓電換能器來(lái)進(jìn)行信號(hào)的發(fā)射和接收。超聲導(dǎo)波換能器固定于鋼軌上,發(fā)射器向相應(yīng)的壓電換能器發(fā)出電壓信號(hào),通過(guò)逆壓電效應(yīng),換能器中的壓電晶體發(fā)生機(jī)械變形,產(chǎn)生振動(dòng),即向鋼軌中發(fā)出導(dǎo)波信號(hào),導(dǎo)波經(jīng)過(guò)鋼軌傳播到接收端,同理,經(jīng)過(guò)正壓電效應(yīng),換能器產(chǎn)生電壓信號(hào),并由導(dǎo)波探頭接收。采用經(jīng)典的一發(fā)兩收裝置,令兩個(gè)接收器位于發(fā)射器的同側(cè),并記錄兩個(gè)接收器間的距離和接收電壓信號(hào)的時(shí)間差,就可以得到超聲導(dǎo)波在鋼軌中的傳播速度,進(jìn)而感知鋼軌的溫度應(yīng)力狀態(tài),該過(guò)程的關(guān)鍵就是對(duì)接收端接收到的導(dǎo)波波形的精確控制。經(jīng)過(guò)前期的有限元仿真對(duì)比分析,文中最后選擇壓電片直接耦合的導(dǎo)波換能器。

    2 時(shí)間反轉(zhuǎn)法方案設(shè)計(jì)

    時(shí)間反轉(zhuǎn)法是指將換能器接收到的時(shí)域信號(hào)進(jìn)行時(shí)間反轉(zhuǎn)放大后,再在接收換能器處重新發(fā)射,使信號(hào)后到先發(fā),先到后發(fā),根據(jù)聲的互易性原理,最終聲波能量將自適應(yīng)聚焦于原激發(fā)換能器處,實(shí)現(xiàn)聲源信號(hào)的重構(gòu)。不用掌握換能器特性及信號(hào)傳播的介質(zhì)屬性,時(shí)間反轉(zhuǎn)法就可以實(shí)現(xiàn)聲波的自適應(yīng)聚焦和檢測(cè),克服了相控聚焦法需要預(yù)先計(jì)算各通道延遲參數(shù)的缺點(diǎn),可有效提高檢測(cè)效率。

    時(shí)間反轉(zhuǎn)過(guò)程中,信號(hào)反轉(zhuǎn)聚焦的前提條件是超聲導(dǎo)波傳播的結(jié)構(gòu)必須是線性的,如果不存在結(jié)構(gòu)損傷或缺陷,最終得到的信號(hào)是各種模態(tài)的超聲導(dǎo)波信號(hào)的聚焦,與原激勵(lì)信號(hào)只有幅值上的差異,波形基本一致;如果存在結(jié)構(gòu)損傷或缺陷,結(jié)構(gòu)內(nèi)部的幾何邊界條件會(huì)發(fā)生變化,當(dāng)導(dǎo)波通過(guò)損傷和缺陷時(shí)會(huì)發(fā)生透射和散射等現(xiàn)象,而各模態(tài)導(dǎo)波對(duì)結(jié)構(gòu)損傷和缺陷的響應(yīng)則會(huì)體現(xiàn)在反轉(zhuǎn)聚焦信號(hào)的主波包中,因此,超聲導(dǎo)波和時(shí)間反轉(zhuǎn)法常用來(lái)檢測(cè)結(jié)構(gòu)內(nèi)部的損傷和缺陷。而本文則利用時(shí)間反轉(zhuǎn)法中聚焦信號(hào)與激勵(lì)信號(hào)波形一致的特點(diǎn)來(lái)獲得鋼軌中超聲導(dǎo)波任意波形的激勵(lì)。

    當(dāng)需要在接收器處得到某接收波形時(shí),目前一般選擇將該種波形信號(hào)直接作為激勵(lì)波形,但由于現(xiàn)實(shí)條件并不理想,再加上超聲導(dǎo)波本來(lái)的多模態(tài)基頻散特性,實(shí)際接收到的波形與需要的波形還有一定差距,會(huì)影響后續(xù)研究的精確度。本文利用時(shí)間反轉(zhuǎn)法獲得任意波形激勵(lì)的示意圖如圖1所示。當(dāng)需要在接收器B處得到接收信號(hào)S1時(shí),將信號(hào)S1作為激勵(lì)施加于接收器B處,在激發(fā)器A處得接收信號(hào)U1,將U1的主波包經(jīng)過(guò)反轉(zhuǎn)放大得信號(hào)U2,再將U2作為激勵(lì)加載到激發(fā)器A上,在接收器B處得到反轉(zhuǎn)聚焦信號(hào)S2,將最初的激勵(lì)信號(hào)S1和聚焦信號(hào)S2的主波包作歸一化處理,進(jìn)行對(duì)比,按照時(shí)間反轉(zhuǎn)法原理,兩個(gè)信號(hào)之間差別應(yīng)該不大,則信號(hào)U2即為需要的激勵(lì)信號(hào)。在時(shí)反過(guò)程中,需求信號(hào)S1可以改變,由不同的S1可以得到不同的U2,從而實(shí)現(xiàn)任意波形的激勵(lì)控制。

    3 有限元仿真過(guò)程參數(shù)設(shè)置

    有限元分析軟件ANSYS能夠有效地分析本文需要的結(jié)構(gòu)-電場(chǎng)的耦合,選擇的鋼軌模型為中國(guó)鐵路采用的CHN60軌,壓電片直接耦合作為壓電換能器,構(gòu)建的模型如圖2所示。

    圖1 利用時(shí)間反轉(zhuǎn)法獲得激勵(lì)波形的示意圖

    圖2 鋼軌和換能器模型

    3.1激勵(lì)頻率與激勵(lì)位置

    激勵(lì)頻率是一個(gè)非常重要的參數(shù),當(dāng)激勵(lì)頻率過(guò)低時(shí)信噪比很低,距離稍遠(yuǎn)就無(wú)法采集到有用信號(hào);當(dāng)激勵(lì)頻率過(guò)大時(shí),會(huì)使導(dǎo)波模態(tài)增多,不利于后續(xù)信號(hào)處理[15],因此,選擇合適的激勵(lì)頻率是必要的。將0~50 kHz分為10組,5 kHz為一個(gè)間隔,進(jìn)行有限元仿真分析,對(duì)比接收波形,最終選擇出35 kHz作為鋼軌中超聲導(dǎo)波的最佳激勵(lì)頻率。

    由于鐵路行車需要,鋼軌的軌頭不能安裝換能器,軌底由于有鋼軌扣件的存在,振動(dòng)在軌底的傳播會(huì)受限導(dǎo)致衰減過(guò)快,致使檢測(cè)距離極劇減小,因此超聲導(dǎo)波換能器只能安裝于軌腰處。通過(guò)對(duì)鋼軌進(jìn)行模態(tài)分析發(fā)現(xiàn)鋼軌的振動(dòng)集中在鋼軌軌腰處,軌頭和軌底的振動(dòng)都很小,這也證明了超聲導(dǎo)波安裝在鋼軌軌腰的必要性。在激勵(lì)超聲導(dǎo)波時(shí),應(yīng)盡量使能量沿著軌腰傳播,使能量足夠集中于軌腰,這樣可以使產(chǎn)生的能量被最大程度地利用。

    3.2壓電耦合

    超聲導(dǎo)波在壓電耦合分析中遵循的壓電方程為:

    因此,在有限元分析計(jì)算中,式(1)中需要設(shè)定的壓電材料參數(shù)為彈性常數(shù)矩陣[c],介電常數(shù)矩陣[ε],壓電應(yīng)力常數(shù)矩陣[e],這也是壓電材料最重要及必不可少的3個(gè)參數(shù),各參數(shù)的值分別如表1,式(2),表2所示。

    表1 彈性常數(shù)矩陣1010Pa

    值得一提的是,設(shè)置材料屬性時(shí),一定要保持單位的一致性,網(wǎng)格劃分以及有限元分析都是在mm-kg-s的單位制環(huán)境中進(jìn)行的,在參數(shù)設(shè)置之前,必須將標(biāo)準(zhǔn)單位制m-kg-s下的參數(shù)轉(zhuǎn)換成單位制mm-kg-s下的參數(shù)。

    3.3網(wǎng)格劃分及計(jì)算步長(zhǎng)

    為保證計(jì)算精度,網(wǎng)格劃分的單元大小Le和有限元瞬態(tài)分析的時(shí)間步長(zhǎng)Δt都必須精確,但二者過(guò)小又會(huì)大大降低仿真效率,故選擇合適的網(wǎng)格大小和計(jì)算時(shí)間步長(zhǎng)尤為重要。根據(jù)經(jīng)驗(yàn)可知,網(wǎng)格劃分單元大小Le應(yīng)滿足以下關(guān)系式:

    而有限元的瞬態(tài)分析時(shí)間步長(zhǎng)Δt則應(yīng)滿足:

    由前期得到的實(shí)驗(yàn)數(shù)據(jù)可知,超聲導(dǎo)波在鋼軌中的傳播速度v≈3 000 m/s,又知λ=v/f,而f=35 kHz,因此,本文在鋼軌超聲導(dǎo)波換能器的仿真中網(wǎng)格單元大小取Le=5 mm,同理,計(jì)算時(shí)間步長(zhǎng)取Δt=1×10-6s。

    表2 壓電矩陣常數(shù)c/m2

    4 有限元仿真與分析

    鋼軌模型長(zhǎng)5 m,激發(fā)換能器A(左)和接收換能器B(右)均設(shè)于鋼軌軌腰處,且兩個(gè)換能器相距0.5 m,激發(fā)器A距離鋼軌左端面2.3 m。假設(shè)需要在接收器B處接收電壓信號(hào)a,本文的信號(hào)a以5個(gè)周期,幅值為1,漢寧窗調(diào)制的正弦信號(hào)為例,如圖3所示。

    圖3 激勵(lì)信號(hào)

    將需要的接收信號(hào)a作為激勵(lì)信號(hào)施加在接收器B上,激勵(lì)頻率為35 kHz,進(jìn)行ANSYS有限元仿真后在激發(fā)器A處得到正向電壓接收信號(hào)b如圖4(a),取正向接收信號(hào)b的主波包進(jìn)行反轉(zhuǎn)放大后得反轉(zhuǎn)信號(hào)c,如圖4(b),將反轉(zhuǎn)信號(hào)c作為新的電壓激勵(lì)信號(hào)施加在激發(fā)器A上進(jìn)行反轉(zhuǎn)仿真。最終在接收器B處得到的反轉(zhuǎn)聚焦信號(hào)d如圖4(c),將激勵(lì)信號(hào)a與反轉(zhuǎn)聚焦信號(hào)d的主波包進(jìn)行歸一化處理之后得到的對(duì)比圖如圖4(d)所示。

    圖4 時(shí)間反轉(zhuǎn)的有限元仿真過(guò)程及結(jié)果

    由圖4(d)中激勵(lì)信號(hào)a和反轉(zhuǎn)聚焦信號(hào)d歸一化處理后的對(duì)比圖可知,聚焦信號(hào)d的波形曲線平滑,不論是從主波包的周期數(shù),還是波形走向方面,聚焦信號(hào)d對(duì)激勵(lì)信號(hào)a的還原度都很高。

    對(duì)比圖4(a)中的正向接收波形b和圖4(c)中的反轉(zhuǎn)聚焦信號(hào)d可見(jiàn),當(dāng)需要獲得接收信號(hào)a時(shí),直接選擇a作為激勵(lì)信號(hào)時(shí)得到的接收波形b也有明顯主波包,并與信號(hào)a很相似,但本文并沒(méi)有采用這種方法,而是選擇利用時(shí)反法先求得激勵(lì)波形,進(jìn)而再獲得接收波形d,是因?yàn)樾盘?hào)波形b與d相比,d的波形與需要得到的波形a更相似。通過(guò)matlab分別求得波形b和a、波形d和a的互相關(guān)曲線及互相關(guān)系數(shù),如圖5所示。

    圖5 波形b和a、d和a的互相關(guān)曲線及互相關(guān)系數(shù)

    互相關(guān)曲線的峰值越高,互相關(guān)系數(shù)越大代表相關(guān)度越高,即波形相似度越高。由圖5可以很明顯看出波形b和a與波形d和a相比,后者的互相關(guān)曲線峰值更高,互相關(guān)系數(shù)更大,因此波形d與a的相似度更高,證明利用時(shí)間反轉(zhuǎn)法獲得的接收波形d精確度更高。

    當(dāng)需要在接收器B處得到接收信號(hào)a時(shí),則可以利用時(shí)反法獲得反轉(zhuǎn)信號(hào)c作為激勵(lì)波形施加在激發(fā)器A處,當(dāng)信號(hào)a改變時(shí),則可以得到不同的激勵(lì)信號(hào)c,因此可以通過(guò)時(shí)反法獲得任意波形的激勵(lì)。

    5 結(jié)論

    本文基于鋼軌超聲導(dǎo)波檢測(cè)技術(shù),借助ANSYS有限元仿真,確定了激勵(lì)頻率、激勵(lì)位置、壓電耦合、仿真步長(zhǎng)等參數(shù),最終利用時(shí)間反轉(zhuǎn)法成功獲得了任意波形的激勵(lì)。仿真結(jié)果表明,和直接將需要的接收波形作為激勵(lì)的方法相比,利用時(shí)間反轉(zhuǎn)法獲得的接收波形與需要的波形相關(guān)度更高,提高了下一步導(dǎo)波速度檢測(cè)的精確度,進(jìn)而提高了鋼軌應(yīng)力監(jiān)測(cè)系統(tǒng)的精確度。

    [1]李旭.用電容式傳感器測(cè)量鋼軌溫度應(yīng)力[D].成都:電子科技大學(xué),2005.

    [2]劉興漢.鋼軌溫度和溫度應(yīng)力[J].西鐵科技,2002(3):30-32.

    [3]閆進(jìn)學(xué).測(cè)標(biāo)測(cè)量無(wú)縫線路鎖定軌溫方法在新線的應(yīng)用[J].鐵道勘查,2007(6):81-83.

    [4]高宏偉.無(wú)縫鋼軌溫度應(yīng)力檢測(cè)技術(shù)研究[D].北京交通大學(xué),2012.

    [5]張兆亭,閆連山,王平,等.基于光纖光柵的鋼軌應(yīng)變測(cè)量關(guān)鍵技術(shù)研究[J].鐵道學(xué)報(bào),2012,34(5):65-69.

    [6]王驍,劉輝,祁欣,等.巴克豪森噪訊無(wú)縫線路應(yīng)力檢測(cè)儀的研制及應(yīng)用[J].北京化工大學(xué)學(xué)報(bào):自然科學(xué)版,2010,37(3):123-126.

    [7]田浩,于石生,趙小瑩,等.利用巴克豪森效應(yīng)測(cè)定鋼軌縱向應(yīng)力[J].材料科學(xué)與工藝,2004,12(2):196-198.

    [8]劉金艷.X射線殘余應(yīng)力的測(cè)量技術(shù)與應(yīng)用研究[D].北京:北京工業(yè)大學(xué),2009.

    [9]Zhu Z H,Post M A,Xu P C.Stress evaluation using ultrasonic interference spectrum of leaky lamb waves[J].Experimental Mechanics,2011,51:971-980.

    [10]Kleitsa D,Kawai K,Shiotani T,et al.Assessment of metal strand wire pre-stress in anchor head by ultrasonics[J].NDT &EInternational,2010,43:547-554.

    [11]Ivan Bartoli,Alessandro Marzani,et al.Modeling wave propagation in damped waveguides of arbitrary cross-section[J]. Journal of Sound and Vibration,2006,295:685-707.

    [12]Alessandro Marzani.Time-transient response for ultrasonic guided waves propagating in damped cylinders[J].International Journal of Solids and Structures,2008,45:6347-6368.

    [13]Mazzotti M,Marzani A,Bartoli I,et al.Guided waves dispersion analysis for prestressed viscoelastic waveguides by means of the SAFE method[J].International Journal of Solids and Structures,2012,49:2359-2372.

    [14]Steve Dixon,Rachel S.Edwards.Detection and depth gauging of transverse cracks in the rail head[J].Ironmaking and Steelmaking,2005,32(5):376-378.

    [15]Joseph L.Rose,Michael J.Avioli,Peter Mudge,et al.Guided wave inspection potential of defects in rail[J].NDT&E International,2004,37:153-161.

    The technique of ultrasonic guided wave with arbitrary waveform excitation

    ZHENG Tian-tian,PAN Xiao-yu
    (School of Mechanical,Electronic and Control Engineering,Beijing Jiao-tong University,Beijing 100044,China)

    High efficiency and accurate seamless rail temperature stress detection and flaw detection is an important technology to ensure the safe operation of high speed railway.In recent years,the rail stress detection technology based on Ultrasonic guided wave has caused more and more attention,in which precise control of guided wave waveform in rail is the key.But at present the accuracy of received waveform is not high.Using ultrasonic guided wave transducer and time reversal method,the excitation signal of arbitrary waveform is gotand the correlation between received waveform andideal waveform is improvedby ANSYS finite element simulation.The simulation results show that the accuracy of the received waveform obtained by time reversal method is improved by about 12%.

    ultrasonic guided wave transducer;time reversal method;waveform;excitation signal

    TN98

    A

    1674-6236(2016)12-0001-04

    2016-02-28稿件編號(hào):201602170

    國(guó)家自然科學(xué)基金(61134003)

    鄭甜甜(1990—),女,河北徐水人,碩士研究生。研究方向:超聲導(dǎo)波理論研究及有限元仿真。

    猜你喜歡
    有限元信號(hào)檢測(cè)
    “不等式”檢測(cè)題
    “一元一次不等式”檢測(cè)題
    “一元一次不等式組”檢測(cè)題
    信號(hào)
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
    電子制作(2018年11期)2018-08-04 03:25:42
    小波變換在PCB缺陷檢測(cè)中的應(yīng)用
    基于LabVIEW的力加載信號(hào)采集與PID控制
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    麻豆国产av国片精品| 国产av在哪里看| 国产亚洲精品久久久久久毛片| 国产精品久久久久成人av| 9191精品国产免费久久| 国产高清视频在线播放一区| 日本黄色日本黄色录像| 久久久精品欧美日韩精品| 亚洲av日韩精品久久久久久密| 国产成年人精品一区二区 | 99久久精品国产亚洲精品| 两人在一起打扑克的视频| 一级毛片女人18水好多| 亚洲中文日韩欧美视频| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看| 人成视频在线观看免费观看| 国产成人啪精品午夜网站| 咕卡用的链子| 国产男靠女视频免费网站| 99在线人妻在线中文字幕| 国产精品偷伦视频观看了| 亚洲精品国产一区二区精华液| 国产成人精品久久二区二区免费| 久久香蕉精品热| 国产一区在线观看成人免费| 亚洲av第一区精品v没综合| 黄色片一级片一级黄色片| 免费av中文字幕在线| 欧美人与性动交α欧美软件| 欧美av亚洲av综合av国产av| 国产精品日韩av在线免费观看 | 久久伊人香网站| 欧美 亚洲 国产 日韩一| av中文乱码字幕在线| 19禁男女啪啪无遮挡网站| 18禁裸乳无遮挡免费网站照片 | 黄频高清免费视频| 69av精品久久久久久| 成人国语在线视频| 香蕉久久夜色| 成人国产一区最新在线观看| 黄色视频不卡| 18禁黄网站禁片午夜丰满| 高潮久久久久久久久久久不卡| 亚洲五月色婷婷综合| 日本 av在线| 男女之事视频高清在线观看| 亚洲国产中文字幕在线视频| 热99国产精品久久久久久7| 午夜免费成人在线视频| 90打野战视频偷拍视频| 桃色一区二区三区在线观看| 色播在线永久视频| 成人国产一区最新在线观看| 五月开心婷婷网| 欧美人与性动交α欧美精品济南到| 制服人妻中文乱码| 一边摸一边做爽爽视频免费| www.999成人在线观看| 成人手机av| 黑人欧美特级aaaaaa片| 12—13女人毛片做爰片一| 一二三四社区在线视频社区8| 在线观看免费日韩欧美大片| 看免费av毛片| 久久人人97超碰香蕉20202| 亚洲精品粉嫩美女一区| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添小说| 久久婷婷成人综合色麻豆| 99精品欧美一区二区三区四区| a在线观看视频网站| 欧美日韩亚洲国产一区二区在线观看| 亚洲午夜精品一区,二区,三区| 天天影视国产精品| 高清黄色对白视频在线免费看| 国产成人欧美在线观看| 美国免费a级毛片| 国产男靠女视频免费网站| 国产av精品麻豆| 亚洲欧美激情综合另类| 91九色精品人成在线观看| 久久久久亚洲av毛片大全| 久久久久久久午夜电影 | 巨乳人妻的诱惑在线观看| 免费女性裸体啪啪无遮挡网站| 精品无人区乱码1区二区| 黄色a级毛片大全视频| 黑人巨大精品欧美一区二区mp4| 啦啦啦免费观看视频1| 亚洲九九香蕉| 欧洲精品卡2卡3卡4卡5卡区| 欧美精品亚洲一区二区| 欧美黑人欧美精品刺激| 制服诱惑二区| 午夜激情av网站| avwww免费| 欧美在线一区亚洲| 丰满迷人的少妇在线观看| 亚洲中文字幕日韩| 99在线人妻在线中文字幕| 日韩有码中文字幕| 亚洲专区字幕在线| 国产真人三级小视频在线观看| 51午夜福利影视在线观看| 中文字幕色久视频| 国产精品一区二区精品视频观看| 99精品欧美一区二区三区四区| 69精品国产乱码久久久| 好看av亚洲va欧美ⅴa在| 午夜影院日韩av| 免费观看人在逋| 中出人妻视频一区二区| 国产精品亚洲一级av第二区| 亚洲av成人av| 久久久精品国产亚洲av高清涩受| 欧美激情久久久久久爽电影 | a级毛片在线看网站| www.熟女人妻精品国产| 欧美成人午夜精品| 黄色视频不卡| 成人三级做爰电影| 变态另类成人亚洲欧美熟女 | 高清毛片免费观看视频网站 | 两个人免费观看高清视频| 国产黄色免费在线视频| 真人做人爱边吃奶动态| 国产不卡一卡二| 麻豆av在线久日| 淫妇啪啪啪对白视频| 久久久久久人人人人人| 日本免费a在线| av网站在线播放免费| 777久久人妻少妇嫩草av网站| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 成人18禁高潮啪啪吃奶动态图| 超碰97精品在线观看| 在线天堂中文资源库| 热99re8久久精品国产| 91九色精品人成在线观看| 国产真人三级小视频在线观看| 一级,二级,三级黄色视频| 搡老乐熟女国产| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 成人三级做爰电影| 国产成人欧美| 日本a在线网址| 午夜久久久在线观看| 日本黄色视频三级网站网址| 国产av精品麻豆| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 88av欧美| 久久久久国产一级毛片高清牌| 欧美不卡视频在线免费观看 | 亚洲av片天天在线观看| 在线视频色国产色| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 国产男靠女视频免费网站| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲欧美在线一区二区| 夫妻午夜视频| 久久婷婷成人综合色麻豆| 大陆偷拍与自拍| 午夜久久久在线观看| 欧美日韩黄片免| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 热re99久久精品国产66热6| av有码第一页| 国产人伦9x9x在线观看| 亚洲美女黄片视频| 一区在线观看完整版| 亚洲午夜理论影院| 1024香蕉在线观看| 国产单亲对白刺激| 亚洲一区二区三区色噜噜 | 久久青草综合色| 91成人精品电影| 18美女黄网站色大片免费观看| 国产99久久九九免费精品| 亚洲少妇的诱惑av| 91国产中文字幕| 色综合婷婷激情| 国产精品香港三级国产av潘金莲| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 午夜亚洲福利在线播放| 亚洲精品av麻豆狂野| 亚洲自拍偷在线| 麻豆av在线久日| 一级毛片高清免费大全| 另类亚洲欧美激情| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 男女下面插进去视频免费观看| 国产日韩一区二区三区精品不卡| 免费少妇av软件| 桃红色精品国产亚洲av| 欧美日韩瑟瑟在线播放| 美女高潮到喷水免费观看| 狂野欧美激情性xxxx| 国产亚洲av高清不卡| 国产一区二区在线av高清观看| 狂野欧美激情性xxxx| 久久精品91蜜桃| 国产精品av久久久久免费| 国产人伦9x9x在线观看| 日本欧美视频一区| 国产精品秋霞免费鲁丝片| a级毛片黄视频| 黄色视频不卡| 欧美中文综合在线视频| 高清欧美精品videossex| 一a级毛片在线观看| 视频在线观看一区二区三区| 搡老乐熟女国产| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| 日日爽夜夜爽网站| 国产精品免费一区二区三区在线| 精品欧美一区二区三区在线| 亚洲专区中文字幕在线| www.www免费av| 咕卡用的链子| 99国产精品免费福利视频| 日本vs欧美在线观看视频| 一级黄色大片毛片| 香蕉丝袜av| 国产人伦9x9x在线观看| 国产av在哪里看| 久久久久国内视频| 久久国产亚洲av麻豆专区| 一级a爱片免费观看的视频| 桃红色精品国产亚洲av| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品美女久久久久99蜜臀| 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 美女扒开内裤让男人捅视频| 日韩欧美国产一区二区入口| 一区在线观看完整版| 老熟妇乱子伦视频在线观看| 国产有黄有色有爽视频| 男男h啪啪无遮挡| 国产蜜桃级精品一区二区三区| 国产成人一区二区三区免费视频网站| 亚洲三区欧美一区| 国产亚洲精品久久久久5区| 最新美女视频免费是黄的| 精品久久久久久,| 欧美精品一区二区免费开放| 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 亚洲欧洲精品一区二区精品久久久| 乱人伦中国视频| 国产黄色免费在线视频| a级毛片黄视频| 国产免费现黄频在线看| 精品久久久久久久久久免费视频 | 国产精品美女特级片免费视频播放器 | 亚洲欧美日韩另类电影网站| 最近最新中文字幕大全免费视频| 精品国产国语对白av| 91麻豆精品激情在线观看国产 | 黄色女人牲交| 一级片'在线观看视频| 国产欧美日韩一区二区三区在线| 亚洲av熟女| 久久国产亚洲av麻豆专区| 亚洲激情在线av| 老司机午夜福利在线观看视频| 国产精品综合久久久久久久免费 | 国产欧美日韩一区二区精品| 午夜两性在线视频| 中亚洲国语对白在线视频| 午夜91福利影院| www.999成人在线观看| 欧美中文综合在线视频| 国产真人三级小视频在线观看| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲高清精品| 两性夫妻黄色片| 成年版毛片免费区| 精品免费久久久久久久清纯| 亚洲中文日韩欧美视频| 国产一区二区三区综合在线观看| 国产aⅴ精品一区二区三区波| 国产在线观看jvid| 性欧美人与动物交配| 少妇粗大呻吟视频| 午夜精品在线福利| 动漫黄色视频在线观看| 十八禁人妻一区二区| 日韩国内少妇激情av| 亚洲av日韩精品久久久久久密| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三卡| 亚洲成国产人片在线观看| 好男人电影高清在线观看| 超色免费av| 国产亚洲精品一区二区www| 亚洲精品粉嫩美女一区| av国产精品久久久久影院| av视频免费观看在线观看| 成人国语在线视频| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| 女人被躁到高潮嗷嗷叫费观| 亚洲色图 男人天堂 中文字幕| 九色亚洲精品在线播放| 免费观看人在逋| 丝袜美足系列| 美女高潮到喷水免费观看| 国产精品一区二区免费欧美| 久久中文看片网| 精品久久久久久成人av| 啦啦啦免费观看视频1| 免费不卡黄色视频| 免费看a级黄色片| bbb黄色大片| 精品国产亚洲在线| 嫁个100分男人电影在线观看| 日韩精品中文字幕看吧| 大型av网站在线播放| 国产成+人综合+亚洲专区| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 嫩草影院精品99| 看免费av毛片| 黄色视频不卡| 国产成人免费无遮挡视频| 99精品欧美一区二区三区四区| www.熟女人妻精品国产| 亚洲视频免费观看视频| 黄色丝袜av网址大全| av欧美777| 国产单亲对白刺激| 在线观看一区二区三区激情| 99久久综合精品五月天人人| 亚洲国产毛片av蜜桃av| 少妇的丰满在线观看| 神马国产精品三级电影在线观看 | 无遮挡黄片免费观看| 国产国语露脸激情在线看| 身体一侧抽搐| 欧美日韩视频精品一区| 老司机靠b影院| 成年人黄色毛片网站| 黄网站色视频无遮挡免费观看| 久久天躁狠狠躁夜夜2o2o| 9热在线视频观看99| 久久99一区二区三区| 日韩免费高清中文字幕av| 成人免费观看视频高清| 国产亚洲精品综合一区在线观看 | 欧美乱色亚洲激情| 欧美日韩亚洲综合一区二区三区_| 成人精品一区二区免费| 亚洲av电影在线进入| 一级毛片高清免费大全| 国产精品自产拍在线观看55亚洲| 精品久久久精品久久久| 久久人妻福利社区极品人妻图片| 亚洲av片天天在线观看| 欧美日韩亚洲综合一区二区三区_| 在线观看舔阴道视频| 久久久国产精品麻豆| 麻豆一二三区av精品| 成年人免费黄色播放视频| 午夜免费鲁丝| 中文字幕人妻丝袜一区二区| 人妻丰满熟妇av一区二区三区| 亚洲成av片中文字幕在线观看| 久久影院123| 亚洲成人国产一区在线观看| 欧美日韩亚洲国产一区二区在线观看| 成人国产一区最新在线观看| 在线观看日韩欧美| 亚洲全国av大片| 69精品国产乱码久久久| 成人国产一区最新在线观看| 免费在线观看亚洲国产| 亚洲熟女毛片儿| 亚洲av电影在线进入| 丝袜美足系列| 欧美激情极品国产一区二区三区| 国产高清videossex| 在线观看一区二区三区激情| 免费久久久久久久精品成人欧美视频| 少妇被粗大的猛进出69影院| 免费av中文字幕在线| 国产av一区二区精品久久| 亚洲国产中文字幕在线视频| 91麻豆精品激情在线观看国产 | 亚洲人成77777在线视频| 久久久久久久久中文| 午夜福利免费观看在线| 韩国av一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 亚洲一区二区三区不卡视频| 18禁美女被吸乳视频| 女性被躁到高潮视频| 日韩av在线大香蕉| 香蕉久久夜色| 狠狠狠狠99中文字幕| 国产精品国产高清国产av| 久久99一区二区三区| 曰老女人黄片| 1024视频免费在线观看| 不卡一级毛片| 美女 人体艺术 gogo| 99国产极品粉嫩在线观看| 国产97色在线日韩免费| 亚洲第一av免费看| 免费久久久久久久精品成人欧美视频| 91麻豆精品激情在线观看国产 | 国产av精品麻豆| 无限看片的www在线观看| 欧美激情高清一区二区三区| 欧美日韩av久久| 亚洲一区二区三区欧美精品| 动漫黄色视频在线观看| 一级作爱视频免费观看| 最好的美女福利视频网| 大码成人一级视频| 久久精品aⅴ一区二区三区四区| 亚洲精品国产色婷婷电影| 在线国产一区二区在线| 亚洲,欧美精品.| 国产在线精品亚洲第一网站| 91精品三级在线观看| 亚洲 欧美一区二区三区| 成人亚洲精品一区在线观看| 免费在线观看视频国产中文字幕亚洲| 国产成人精品久久二区二区91| 一二三四在线观看免费中文在| 桃红色精品国产亚洲av| 老熟妇仑乱视频hdxx| 热99re8久久精品国产| 亚洲成人免费av在线播放| 99久久综合精品五月天人人| 午夜免费观看网址| 一级毛片女人18水好多| 黄色视频,在线免费观看| 99国产精品一区二区三区| 亚洲国产欧美一区二区综合| 亚洲 欧美 日韩 在线 免费| 三级毛片av免费| 波多野结衣av一区二区av| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 精品第一国产精品| 国产一区二区激情短视频| 国产亚洲精品一区二区www| 国产成人精品无人区| 一边摸一边做爽爽视频免费| 99热国产这里只有精品6| 午夜影院日韩av| 熟女少妇亚洲综合色aaa.| 宅男免费午夜| 男男h啪啪无遮挡| 免费女性裸体啪啪无遮挡网站| 777久久人妻少妇嫩草av网站| 成人手机av| 久久人人97超碰香蕉20202| 国产成人欧美在线观看| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人| 91av网站免费观看| 国产精品99久久99久久久不卡| 中文欧美无线码| 国产深夜福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美黑人欧美精品刺激| 在线观看舔阴道视频| 国产精品久久久久久人妻精品电影| 日韩中文字幕欧美一区二区| 免费在线观看黄色视频的| 最好的美女福利视频网| 男女高潮啪啪啪动态图| 国产午夜精品久久久久久| 一本大道久久a久久精品| av网站免费在线观看视频| 国产无遮挡羞羞视频在线观看| 久久这里只有精品19| 99re在线观看精品视频| 91老司机精品| 欧美日韩亚洲综合一区二区三区_| 国产乱人伦免费视频| 国产欧美日韩一区二区精品| 久久久国产欧美日韩av| 欧洲精品卡2卡3卡4卡5卡区| 国产黄a三级三级三级人| 无限看片的www在线观看| 亚洲一区二区三区欧美精品| 男人舔女人的私密视频| av网站在线播放免费| 99精品欧美一区二区三区四区| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 大型黄色视频在线免费观看| 校园春色视频在线观看| 91精品三级在线观看| 亚洲一区二区三区不卡视频| www.www免费av| 嫩草影视91久久| 校园春色视频在线观看| 99在线人妻在线中文字幕| 亚洲自偷自拍图片 自拍| 国产精品久久视频播放| 自线自在国产av| 久久久国产一区二区| 亚洲精品美女久久久久99蜜臀| 国产成+人综合+亚洲专区| 日韩高清综合在线| 大型黄色视频在线免费观看| 欧美另类亚洲清纯唯美| 天堂√8在线中文| 精品卡一卡二卡四卡免费| av超薄肉色丝袜交足视频| 日韩免费av在线播放| 亚洲欧美精品综合久久99| 欧美激情 高清一区二区三区| 亚洲国产欧美网| 亚洲五月天丁香| 亚洲 国产 在线| 午夜久久久在线观看| 国产午夜精品久久久久久| 欧美国产精品va在线观看不卡| 亚洲五月婷婷丁香| 亚洲美女黄片视频| 欧美久久黑人一区二区| 国内毛片毛片毛片毛片毛片| 丰满饥渴人妻一区二区三| 国产区一区二久久| 亚洲,欧美精品.| 欧美乱码精品一区二区三区| 久久久久九九精品影院| 久久人人精品亚洲av| 91大片在线观看| av在线天堂中文字幕 | 啦啦啦 在线观看视频| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 国产在线观看jvid| xxxhd国产人妻xxx| 亚洲一码二码三码区别大吗| 国产成人欧美| 欧美日韩亚洲国产一区二区在线观看| 19禁男女啪啪无遮挡网站| 手机成人av网站| 别揉我奶头~嗯~啊~动态视频| 亚洲国产看品久久| 国产精品一区二区免费欧美| 久久人妻熟女aⅴ| 激情视频va一区二区三区| 夜夜夜夜夜久久久久| 99热只有精品国产| 午夜福利影视在线免费观看| 成年人黄色毛片网站| 亚洲第一av免费看| 欧美最黄视频在线播放免费 | 咕卡用的链子| 大香蕉久久成人网| 亚洲成人精品中文字幕电影 | 一区福利在线观看| 看黄色毛片网站| 日韩大码丰满熟妇| 精品久久久久久,| 精品国产乱子伦一区二区三区| 国产成人精品久久二区二区91| 三上悠亚av全集在线观看| svipshipincom国产片| 久久久国产精品麻豆| 亚洲精品成人av观看孕妇| aaaaa片日本免费| 精品乱码久久久久久99久播| 国产伦一二天堂av在线观看| 十八禁人妻一区二区| 亚洲国产毛片av蜜桃av| 脱女人内裤的视频| av超薄肉色丝袜交足视频| 午夜免费成人在线视频| 精品久久久久久久久久免费视频 | 亚洲av片天天在线观看| 黄片大片在线免费观看| 日本欧美视频一区| 91精品国产国语对白视频| 麻豆av在线久日| 亚洲美女黄片视频| 脱女人内裤的视频| 国产精品久久视频播放| 久久久国产欧美日韩av| 免费在线观看日本一区| 99国产精品一区二区蜜桃av| 少妇粗大呻吟视频| 在线播放国产精品三级| 黄色视频,在线免费观看| 国产99白浆流出| 美女扒开内裤让男人捅视频| 亚洲人成77777在线视频| 久久久久久免费高清国产稀缺| 欧美日韩亚洲国产一区二区在线观看| videosex国产| 动漫黄色视频在线观看| 欧美另类亚洲清纯唯美| 日韩视频一区二区在线观看| 色哟哟哟哟哟哟| 亚洲国产中文字幕在线视频|