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

    電磁超聲換能器線圈設(shè)計(jì)與提高換能效率研究*

    2016-03-22 02:26:25范吉志吳運(yùn)新石文澤譚良辰中南大學(xué)高性能復(fù)雜制造國家重點(diǎn)實(shí)驗(yàn)室長沙410083中南大學(xué)機(jī)電工程學(xué)院長沙410083
    傳感技術(shù)學(xué)報(bào) 2016年1期
    關(guān)鍵詞:有限元

    范吉志,吳運(yùn)新*,石文澤,龔 海,譚良辰(1.中南大學(xué)高性能復(fù)雜制造國家重點(diǎn)實(shí)驗(yàn)室,長沙410083;2.中南大學(xué)機(jī)電工程學(xué)院,長沙410083)

    ?

    電磁超聲換能器線圈設(shè)計(jì)與提高換能效率研究*

    范吉志1,2,吳運(yùn)新1,2*,石文澤1,2,龔海1,2,譚良辰1,2
    (1.中南大學(xué)高性能復(fù)雜制造國家重點(diǎn)實(shí)驗(yàn)室,長沙410083;2.中南大學(xué)機(jī)電工程學(xué)院,長沙410083)

    摘要:提出了一種提高電磁超聲換能器EMAT(Electromagnetic Acoustic Transducer)換能效率的方法,首先在COMSOL Multiphysics軟件中建立了電磁超聲無損檢測激發(fā)探頭有限元模型,然后研究了電磁超聲換能器中雙層疊加螺旋線圈的基板厚度、銅箔厚度以及提離距離對換能效率的影響。結(jié)果表明:減小基板厚度以及提離距離,增加線圈銅箔厚度可以提高EMAT換能效率。最后通過實(shí)驗(yàn)進(jìn)行驗(yàn)證,仿真結(jié)果與實(shí)驗(yàn)基本一致。關(guān)鍵詞:電磁超聲換能器;有限元;螺旋線圈;換能效率

    電磁超聲無損檢測技術(shù)由于具有不需要耦合介質(zhì),檢測溫度最高可達(dá)1 000℃,對一般的被探工件表面不需經(jīng)過特殊清理,可直接進(jìn)行無損探傷等優(yōu)點(diǎn),已經(jīng)成為國際上快速發(fā)展的一項(xiàng)新的檢測手段,是無損檢測領(lǐng)域發(fā)展較快且較前沿的技術(shù)[1-2]。近年來,國內(nèi)外學(xué)者主要的研究方向已經(jīng)集中在將電磁超聲技術(shù)推向工業(yè)應(yīng)用階段。JafariShapoorabadi采用考慮到EMAT線圈之間的集膚效應(yīng)和鄰近效應(yīng)源電流密度的完整方程,并利用有限元方法分析了EMAT換能過程中的電流和洛倫茲力特性,改進(jìn)之后的方法具有更好的穩(wěn)定性和收斂性[3-4]。英國華威大學(xué)Steve Dixon團(tuán)隊(duì)對高溫電磁超聲無損檢測技術(shù)進(jìn)行了大量的研究,利用電磁鐵及陶瓷封裝間距線圈實(shí)現(xiàn)高溫環(huán)境下的應(yīng)用[5-8]。張廣純教授在國內(nèi)最早對EMAT的應(yīng)用技術(shù)進(jìn)行了研究,可以實(shí)現(xiàn)對500℃以下的鋼板進(jìn)行在線的缺陷檢測,并設(shè)計(jì)了大功率脈沖功率放大器來驅(qū)動(dòng)激發(fā)線圈,以便可以得到更高的激發(fā)能量。哈爾濱工業(yè)大學(xué)王淑娟教授科研團(tuán)隊(duì)設(shè)計(jì)了對列車輪對踏面的探傷裝置[9]以及利用Lamb波對鋁合金板材進(jìn)行無損檢測的裝置[10]。

    電磁超聲無損檢測探頭是一種新型UT發(fā)射和接收的直接裝置,一般稱之為換能器,是電磁超聲探傷系統(tǒng)中的核心部件[11]。EMAT的局限性在于換能效率相比其它超聲激發(fā)形式低很多,并且換能機(jī)理較復(fù)雜,通過對EMAT的結(jié)構(gòu)參數(shù)以及尺寸參數(shù)進(jìn)行合理優(yōu)化設(shè)計(jì),提高其換能效率[12]。上述研究已經(jīng)成為提升EMAT無損檢測能力中重要的手段[13]。EMAT的換能效率跟線圈的形狀密切相關(guān)[14-15],為此,利用COMSOL Multiphysics軟件研究了改變雙層疊加線圈結(jié)構(gòu)尺寸對EMAT換能效率的影響,并與實(shí)驗(yàn)結(jié)果進(jìn)行比較。

    1 電磁超聲原理

    電磁超聲的換能機(jī)理包括洛倫茲力機(jī)理、磁致伸縮機(jī)理和磁化力機(jī)理。目前的研究表明在非鐵磁性材料中主要是洛倫茲力機(jī)理的存在,而在鐵磁性材料中不僅包含洛倫茲力機(jī)理還包含了洛倫茲力機(jī)理,而且還包括了磁致伸縮機(jī)理和磁化力機(jī)理,但是一般情況下磁化力都比較小,在研究過程中可以忽略。

    在非鐵磁性材料中僅有洛倫茲力來激發(fā)超聲波,其激發(fā)過程為:EMAT線圈中通過高頻、大功率的電流,其在被測試樣內(nèi)部產(chǎn)生與之流向相反的感生渦流,此時(shí)被測試樣的集膚層內(nèi)相當(dāng)于有電流流過,根據(jù)安培定律可知,永磁鐵的靜態(tài)磁場與感生渦流相互作用產(chǎn)生洛倫茲力,由于感生渦流為高頻變化的,為此產(chǎn)生的洛倫茲力在試樣的表面產(chǎn)生,并向內(nèi)部傳播。洛倫茲力機(jī)理的接收過程與發(fā)射過程互為逆過程。換能過程可表述如下:

    式中:H為磁場強(qiáng)度;JC為激發(fā)電流密度;B為磁通密度;EE為交變磁場產(chǎn)生的交變電場強(qiáng)度;JE為交變電流密度;μr為導(dǎo)體相對磁導(dǎo)率;fL為洛倫茲力密度;σ為導(dǎo)體電導(dǎo)率。

    2 EMAT探頭的建模

    本文主要研究基于洛倫茲力機(jī)理的EMAT所激發(fā)的超聲橫波信號(hào),其探頭主要是由提供偏置磁場的圓柱形磁體、線圈以及非鐵磁性材料的被測試樣組成的。根據(jù)計(jì)算建立如圖1所示的EMAT發(fā)射探頭的二維軸對稱模型,圖中顯示的是隱藏了空氣域之后的EMAT有限元模型剩余部分。為了能夠使線圈的尺寸更加精確,更容易控制其結(jié)構(gòu)參數(shù),本文采用印刷電路板PCB(Printed Circuit Board)作為EMAT的螺旋線圈,為了增加信號(hào)的強(qiáng)度,本文采用雙層PCB電路板的結(jié)構(gòu)作為研究對象。EMAT結(jié)構(gòu)參數(shù)如下:EMAT線圈的半徑(外徑)r21,EMAT線圈的內(nèi)圈半徑r11,EMAT線圈中銅箔的厚度h1,PCB電路板基板厚度D1,EMAT的提離距離l1以及永磁體與線圈的距離l2。

    圖1 電磁超聲探頭結(jié)構(gòu)參數(shù)示意圖

    EMAT中探頭的雙層螺旋線圈和被測鋁板試樣的尺寸參數(shù)等如表1所示。

    表1 雙層螺旋線圈有限元模型的結(jié)構(gòu)參數(shù)

    在EMAT二維軸對稱有限元計(jì)算過程中設(shè)定的螺旋線圈和被測鋁板試樣的材料參數(shù)如下:銅箔磁導(dǎo)率4π×10-7H/m,銅箔電導(dǎo)率2.667×10-7S/m,鋁板磁導(dǎo)率4π×10-7H/m,鋁板電導(dǎo)率2.667×10-7S/m,試樣彈性模量71 GPa,試樣泊松比0.33,剩磁密度1.21 T,磁鐵相對磁導(dǎo)率5.31,試樣相對磁導(dǎo)率1,空氣相對磁導(dǎo)率1。

    電磁超聲無損檢測有限元仿真模擬中,網(wǎng)格大小以及計(jì)算時(shí)間步長的選擇,對仿真結(jié)果有很大的影響。經(jīng)過計(jì)算發(fā)現(xiàn)當(dāng)試樣內(nèi)部的網(wǎng)格大小為超聲波波長的1/20時(shí),超聲波在試樣內(nèi)部的傳播有限元計(jì)算結(jié)果已經(jīng)完全收斂。為了計(jì)算的準(zhǔn)確性,被測試樣內(nèi)部的集膚深度內(nèi)應(yīng)至少有6個(gè)~7個(gè)網(wǎng)格,為此要對集膚層內(nèi)的網(wǎng)格進(jìn)行細(xì)化,圖2為對二維軸對稱模型以及線圈和集膚層進(jìn)行網(wǎng)格劃分后的圖示。

    圖2 網(wǎng)格劃分后的有限元模型

    對于上述多物理場有限元模型的求解,永磁體的靜態(tài)偏置磁場為穩(wěn)定求解,螺旋線圈的動(dòng)態(tài)磁場以及機(jī)械力場均為瞬態(tài)求解。瞬態(tài)求解過程中需要對于求解時(shí)間步長以及其相對誤差、絕對誤差進(jìn)行設(shè)定。本文選用信號(hào)的頻率為1 MHz,根據(jù)對有限元結(jié)果的分析,可知當(dāng)設(shè)定的最大時(shí)間步長為1/1 000f時(shí)有限元仿真結(jié)果的收斂性較好,當(dāng)設(shè)定的絕對誤差為0.01時(shí)有限元模型才能順利的完成計(jì)算。

    一般的EMAT的激勵(lì)信號(hào)為窄帶或調(diào)制tone burst信號(hào),該信號(hào)由信號(hào)發(fā)生器產(chǎn)生并經(jīng)脈沖功率放大器放大之后得到。調(diào)制的tone burst信號(hào)利用下式進(jìn)行表征:

    電流激勵(lì)信號(hào)的波形如圖3所示。

    圖3 螺旋線圈所施加的tone burst電流信號(hào)的波形

    由于電磁超聲激發(fā)過程的有限元模型中涉及到靜態(tài)磁場、螺旋線圈的脈沖感生渦流場以及機(jī)械力場,多場之間需要設(shè)置耦合變量。將靜態(tài)偏置磁場中求解所得到的靜態(tài)磁感應(yīng)強(qiáng)度和脈沖渦流場所得到的渦流密度的乘積,設(shè)置為機(jī)械力場中被測試樣的集膚層內(nèi)所受到的體力。

    3 激發(fā)探頭結(jié)構(gòu)尺寸參數(shù)對換能效率影響的仿真分析

    3.1基板厚度對換能效率的影響

    選取PCB的基板厚度D1為0.3 mm至2 mm之間的7組數(shù)據(jù),為了消除由于基板厚度的變化導(dǎo)致永磁體與被測試樣距離的增加,設(shè)定永磁體與被測試樣之間的距離即基板厚度D1與磁鐵線圈距離l2之和保持3 mm不變,磁鐵、空氣域的尺寸參數(shù)以及線圈的內(nèi)外徑具體參照表1。通過有限元模型得到EMAT的等效阻抗,計(jì)算在不同的PCB厚度的參數(shù)影響下,通過螺旋線圈中的激勵(lì)的強(qiáng)度。

    經(jīng)過有限元模型的計(jì)算,通過提取被測試樣的集膚層距離被測試樣表面0.01 mm處與線圈平行的直線上的渦流密度,求取在每個(gè)時(shí)刻上的平均值,作為EMAT線圈在此刻所產(chǎn)生的感生渦流密度。同樣的方法得到r軸與z軸方向振動(dòng)幅值的平均值,作為此刻所激發(fā)的超聲波信號(hào)強(qiáng)度,并提取不同PCB基板厚度下所激發(fā)超聲波的峰值,以此來評(píng)判電磁超聲激發(fā)探頭的換能效率。如圖4所示為不同PCB基板厚度的EMAT所產(chǎn)生的感生渦流密度與橫縱波峰值。

    圖4 不同PCB基板厚度感生渦流與超聲波幅值

    由上可知當(dāng)PCB板厚由0.5 mm增加至1.5 mm的過程中,螺旋線圈在被測試樣內(nèi)部產(chǎn)生的感生渦流密度不斷減小,超聲橫波與縱波的峰值分別下降了約23%和20%。

    3.2銅箔厚度對換能效率的影響

    選取銅箔厚度h1為0.035 mm、0.05 mm、0.07 mm和0.105 mm 4組數(shù)據(jù),磁鐵、空氣域的尺寸參數(shù)以及線圈的內(nèi)外徑參照表1。通過有限元模型計(jì)算EMAT的等效阻抗,計(jì)算得到在不同銅箔厚度的影響下,螺旋線圈中激勵(lì)電流的強(qiáng)度。圖5所示為不同銅箔厚度的EMAT所產(chǎn)生的感生渦流密度與橫縱波峰值。

    圖5 不同銅箔厚度下激發(fā)的超聲波信號(hào)

    由上可知,當(dāng)銅箔厚度由0.035 mm增加至0.07 mm的過程中,感生渦流密度不斷的增加。超聲橫波與縱波的峰值均增加了17%。

    3.3提離距離對換能效率的影響

    提離距離由0.5 mm增加至2 mm的過程中,感生渦流密度急劇減小。超聲橫波與縱波的峰值均下降了約73%。

    選取提離距離l1為0.5 mm、1 mm、1.5 mm和2 mm 4組數(shù)據(jù),磁鐵、空氣域的尺寸參數(shù)以及線圈的內(nèi)外徑具體參照表1。通過有限元模型計(jì)算EMAT的等效阻抗,計(jì)算得到不同提離距離的EMAT,螺旋線圈中激勵(lì)電流的強(qiáng)度。圖6所示為不同提離距離的EMAT所產(chǎn)生的感生渦流密度與橫縱波峰值。

    圖6 不同提離距離下激發(fā)的超聲波信號(hào)

    4 實(shí)驗(yàn)驗(yàn)證

    為了驗(yàn)證電磁超聲激發(fā)探頭模型和仿真結(jié)果的正確性,將以實(shí)驗(yàn)的方法對螺旋線圈電磁超聲探頭所激發(fā)的橫波進(jìn)行檢測。通過對采集信號(hào)的處理來衡量超聲波激發(fā)信號(hào)的幅值,用來驗(yàn)證有限元模型中螺旋線圈的參數(shù)對激發(fā)信號(hào)的影響規(guī)律。

    4.1實(shí)驗(yàn)方案

    利用同一個(gè)電磁超聲接收探頭來接收不同螺旋線圈所激發(fā)的超聲波信號(hào),通過對接收探頭所接收的超聲信號(hào)的幅值進(jìn)行對比,實(shí)現(xiàn)不同探頭的結(jié)構(gòu)參數(shù)對電磁超聲激發(fā)探頭所激發(fā)的超聲波的信號(hào)的影響規(guī)律的研究和對有限元模型正確性的驗(yàn)證。

    4.2實(shí)驗(yàn)結(jié)果與分析

    按照實(shí)驗(yàn)方案,分別對利用不同的線圈激發(fā)超聲波,利用同一線圈接收,確保單一變量,圖7為實(shí)驗(yàn)測試現(xiàn)場的圖片以及各組別線圈的實(shí)物圖。

    圖8為在其它參數(shù)不變時(shí),PCB板的厚度由0.5 mm增加至1.5 mm,接收到的超聲波的橫波幅值由145.35 mV減小至103.75 mV,減小了約26%,縱波幅值由43.15 mV減小至32.75 mV,減小了約24%,而在相同情況下有限元計(jì)算的結(jié)果分別為23%和20%,實(shí)驗(yàn)與仿真的結(jié)果有很好的吻合度。

    圖9為銅箔的厚度分別為0.035 mm和0.07 mm,提離距離為0.5 mm的螺旋線圈組成的EMAT所激發(fā)的超聲波的對比,由圖可知隨著激發(fā)線圈銅箔厚度由0.035 mm增加至0.07 mm的過程中,所激發(fā)的超聲波的橫波幅值由101.45 mV增加至114.36 mV,增加了約13%,縱波幅值由29.25 mV增加至32.35 mV,增加了約11%,而有限元仿真的結(jié)果均增加了17%,說明了仿真結(jié)果的正確性。

    圖10為提離距離由0.5 mm增加至2 mm的EMAT所激發(fā)的超聲波信號(hào),對比分析可知,隨著提離距離的不斷增加,EMAT所激勵(lì)的超聲波信號(hào)的幅值逐漸減小,當(dāng)提離距離由0.5 mm增加至2 mm時(shí),超聲信號(hào)的橫波幅值由160.04 mV下降至32.21 mV,下降的比例達(dá)到了80%,縱波幅值由49.55 mV下降至10.65 mV,下降的比例達(dá)到了79%,這與之前的仿真結(jié)果73%比較吻合。

    圖7 實(shí)驗(yàn)測試系統(tǒng)及所用線圈

    圖8 相同線圈在不同基板厚度所激發(fā)的超聲波

    圖9 不同銅箔厚度的螺旋線圈所激發(fā)的超聲波

    圖10 相同線圈在不同提離距離所激發(fā)的超聲波

    5 結(jié)束語

    在建立的電磁超聲有限元模型的基礎(chǔ)上,通過仿真以及實(shí)驗(yàn)研究了EMAT雙層螺旋線圈基板的厚度、銅箔的厚度以及提離距離對激發(fā)的超聲橫波和縱波信號(hào)的影響規(guī)律。分析結(jié)果表明:①雙層螺旋線圈探頭的設(shè)計(jì)中應(yīng)該盡量的減小基板的厚度,以提升EMAT的激發(fā)效率,確保激發(fā)超聲橫波信號(hào)的強(qiáng)度。②螺旋線圈的設(shè)計(jì)中應(yīng)該適量的增加線圈中銅箔的厚度,以提升EMAT的激發(fā)效率。③盡量減小提離距離,嚴(yán)格控制提離距離的大小,避免因?yàn)樘犭x距離的變化導(dǎo)致超聲回波信號(hào)幅值的下降,從而引起對缺陷定量估計(jì)的誤差。

    參考文獻(xiàn):

    [1]Kazys R,Voleisis A,Voleisiene B. High Temperature Ultrasonic Transducers:Review[J]. Ultragarsas(Ultrasound),2008,63(2):7-17.

    [2]江念,王召巴,陳友興,等.電磁超聲檢測鋼板厚度實(shí)驗(yàn)的參數(shù)優(yōu)化[J].傳感技術(shù)學(xué)報(bào),2015,28(4):498-502.

    [3]Shapoorabadi R J,Konrad A,Sinclair A N. Computation of Cur?rent Densities in the Receiving Mode of EMATs[J]. Journal of Ap?plied Physics,2005,97(10Q106):1-3.

    [4]Shapoorabadi R J,Konrad A,Sinclair A N. The Governing Electro?dynamic Equations of EMATs[J]. Journal of Applied Physics, 2005,97(10E102):1-3.

    [5]Hernandez-Valle F,Dixon S. Initial Tests for Designing a High Temperature EMAT with Pulsed Electromagnet[J]. NDT&E Inter?national,2010,43(2):171-175.

    [6]Baillie I,Griffith P,Jian X,et al. Implementing an Ultrasonic In?spection System to Find Surface and Internal Defects in Hot,Mov?ing Steel Using EMATs[J]. Insight,2007,49(2):87-92,.

    [7]Jian X,Dixon S. Enhancement of EMAT and Eddy Current Using a Ferrite Back- Plate[J]. Sensors and Actuators A(Physical),2007,136(1):132-136.

    [8]Hernandez-Valle F,Dixon S. Pulsed Electromagnet EMAT for Ul?trasonic Measurements at Elevated Temperatures[J]. Insight,2011,53(2):96-99.

    [9]王淑娟,趙再新,翟國富.基于電磁超聲的火車車輪裂紋檢測系統(tǒng)[J].儀表技術(shù)與傳感器,2005(11):27-29.

    [10]宮佳鵬,許霽,邱玉,等.基于電磁超聲導(dǎo)波的鋁合金板材缺陷自動(dòng)檢測裝置[J].儀表技術(shù)與傳感器,2011(5):78-81.

    [11]夏志敏.電磁超聲檢測技術(shù)初探[D].武漢:華中科技大學(xué),2007:25-29.

    [12]黃鳳英,周正干.靜態(tài)偏置磁場對電磁超聲換能器靈敏度的影響[J].機(jī)械工程學(xué)報(bào),2011,4(10):1-6.

    [13]陳鵬,韓德來,蔡強(qiáng)富,等.電磁超聲檢測技術(shù)的研究進(jìn)展[J].國外電子測量技術(shù),2012,31(11):18-21.

    [14]Kaltenbacher M,Lerch R,Landes H,et al. Computer Optimiza?tion of Electromagnetic a Coustic Transducers[J]. Ultrasonics Symposium,1998,1(2):1029-1034.

    [15]張端,張帥.閉合磁路電磁式低頻振動(dòng)能量收集裝置[J].傳感技術(shù)學(xué)報(bào),2015,28(7):1091-1096.

    范吉志(1989-),男,安徽阜陽人,碩士研究生,主要從事電磁超聲無損檢測技術(shù)的研究,fjzcsu@163.com;

    吳運(yùn)新(1963-),男,博士,博士生導(dǎo)師、教授,主要從事機(jī)械結(jié)構(gòu)動(dòng)力學(xué)、無損檢測、冶金機(jī)械研究,wuyunxin@csu. edu.cn。

    DSP Implementation of Relative Position and Attitude Calculation*

    FENG Yu,WANG Xiangjun*,CHEN Wenliang
    (MOEMS Education Ministry Key Lab,Tianjin University,Tianjin 300072,China)

    Abstract:To improve the speed and reduce the volume in position and attitude measurement,a real-time image ac?quisition and processing system based on TMS320DM6437 digital signal processor was designed.The Video Process SubSystem of TMS320DM6437 and ping-pong operation are used in processing the raw data acquired by CMOS im?age sensor to reduce the volume of system and improve the frame rate.The software and hardware of the system were debugged and the experiments of position and attitude measurement were carried out.When the frame rate of camera is 40fps,system can complete the real-time solution.Experiment results show that the designed system is reasonable and feasible in applications.

    Key words:pose estimation;DSP;CMOS;ping-pong operation;image processing

    doi:EEACC:6140C10.3969/j.issn.1004-1699.2016.01.007

    收稿日期:2015-08-13修改日期:2015-09-24,

    中圖分類號(hào):TB552

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1004-1699(2016)01-0029-06

    項(xiàng)目來源:國家科技支撐計(jì)劃項(xiàng)目(2014BAF12B01);高性能復(fù)雜制造國家重點(diǎn)實(shí)驗(yàn)室自主研究課題項(xiàng)目(zzyjkt2013-06B)

    猜你喜歡
    有限元
    基于擴(kuò)展有限元的疲勞裂紋擴(kuò)展分析
    非線性感應(yīng)加熱問題的全離散有限元方法
    Extended Fisher-Kolmogorov方程的間斷有限元分析
    TDDH型停車器制動(dòng)過程有限元分析
    基于I-DEAS的履帶起重機(jī)主機(jī)有限元計(jì)算
    基于有限元模型對踝模擬扭傷機(jī)制的探討
    10MN快鍛液壓機(jī)有限元分析
    基于HyperWorks的某重型鑄造橋殼有限元分析及改進(jìn)
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    亚洲国产av新网站| 性色av一级| 国产在线一区二区三区精| 午夜视频国产福利| 天天操日日干夜夜撸| 日本色播在线视频| 我要看日韩黄色一级片| 亚洲精品日本国产第一区| 亚洲美女搞黄在线观看| 国产精品99久久久久久久久| 伊人亚洲综合成人网| 日韩大片免费观看网站| 日韩成人av中文字幕在线观看| 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 伦精品一区二区三区| 交换朋友夫妻互换小说| 爱豆传媒免费全集在线观看| 亚洲欧美日韩卡通动漫| 欧美丝袜亚洲另类| a级片在线免费高清观看视频| 三级国产精品欧美在线观看| 日韩大片免费观看网站| 2021少妇久久久久久久久久久| 国产爽快片一区二区三区| 伦精品一区二区三区| 亚洲精品亚洲一区二区| 成人特级av手机在线观看| 久久国内精品自在自线图片| 久久国产精品男人的天堂亚洲 | 超碰97精品在线观看| 久久这里有精品视频免费| 插逼视频在线观看| 女人精品久久久久毛片| 热re99久久国产66热| 免费看日本二区| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区| 男女边吃奶边做爰视频| 51国产日韩欧美| 男人添女人高潮全过程视频| 精品人妻熟女av久视频| 一级毛片aaaaaa免费看小| 全区人妻精品视频| 成人国产av品久久久| 一区二区三区乱码不卡18| 中文字幕免费在线视频6| 中文资源天堂在线| 国产极品天堂在线| 久久久精品94久久精品| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 国产精品一区二区性色av| 亚洲综合色惰| 99久国产av精品国产电影| 十八禁高潮呻吟视频 | 热99国产精品久久久久久7| 亚洲精品一区蜜桃| 免费观看av网站的网址| 一级片'在线观看视频| 最近手机中文字幕大全| av专区在线播放| 国产女主播在线喷水免费视频网站| 97超视频在线观看视频| 男人爽女人下面视频在线观看| 国产色爽女视频免费观看| 黑丝袜美女国产一区| 免费播放大片免费观看视频在线观看| 校园人妻丝袜中文字幕| 国产爽快片一区二区三区| 看免费成人av毛片| 国产成人91sexporn| 亚洲成人手机| 久久久a久久爽久久v久久| av天堂中文字幕网| 亚洲精品色激情综合| 亚洲精品国产av成人精品| 日日啪夜夜爽| 亚洲欧美精品自产自拍| 黑人高潮一二区| 最后的刺客免费高清国语| 精品99又大又爽又粗少妇毛片| 国产成人一区二区在线| 免费人妻精品一区二区三区视频| 国产综合精华液| 99热这里只有是精品在线观看| av视频免费观看在线观看| 亚洲欧美精品自产自拍| 免费av中文字幕在线| 日韩亚洲欧美综合| 精品视频人人做人人爽| 亚洲精品日韩在线中文字幕| 免费大片18禁| 极品教师在线视频| 成人亚洲精品一区在线观看| 又大又黄又爽视频免费| 亚洲成人手机| av黄色大香蕉| 夫妻性生交免费视频一级片| 日日啪夜夜撸| 97在线视频观看| 91在线精品国自产拍蜜月| 最近2019中文字幕mv第一页| 美女中出高潮动态图| 国产免费一级a男人的天堂| 亚洲精品乱久久久久久| 久久久久网色| 另类亚洲欧美激情| 国产黄频视频在线观看| 啦啦啦中文免费视频观看日本| 一级av片app| 亚洲久久久国产精品| 亚洲精品,欧美精品| 9色porny在线观看| 狠狠精品人妻久久久久久综合| 黄色一级大片看看| 我的女老师完整版在线观看| 一区二区三区乱码不卡18| 在线观看一区二区三区激情| 午夜日本视频在线| av国产精品久久久久影院| 自线自在国产av| 夫妻性生交免费视频一级片| 国产精品三级大全| 晚上一个人看的免费电影| 亚洲,欧美,日韩| 久久99精品国语久久久| 国产成人免费无遮挡视频| 99热这里只有精品一区| 久久鲁丝午夜福利片| 国内揄拍国产精品人妻在线| 一本大道久久a久久精品| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 高清av免费在线| 春色校园在线视频观看| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 国产精品久久久久久av不卡| 日本欧美视频一区| 日韩av不卡免费在线播放| 狂野欧美激情性xxxx在线观看| 久久人人爽人人片av| 午夜福利,免费看| 五月伊人婷婷丁香| 中文字幕精品免费在线观看视频 | 亚洲国产欧美日韩在线播放 | 99久久精品国产国产毛片| 日韩视频在线欧美| 久久99蜜桃精品久久| 青青草视频在线视频观看| 国产欧美日韩一区二区三区在线 | 亚洲国产av新网站| 下体分泌物呈黄色| 国产毛片在线视频| 99热网站在线观看| 最近2019中文字幕mv第一页| 免费大片黄手机在线观看| av福利片在线| 国模一区二区三区四区视频| 精品一品国产午夜福利视频| 日韩一本色道免费dvd| 亚洲欧洲国产日韩| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 男女无遮挡免费网站观看| 男的添女的下面高潮视频| 激情五月婷婷亚洲| 亚洲国产欧美日韩在线播放 | 精品一区二区三区视频在线| 另类精品久久| 欧美日韩综合久久久久久| 少妇人妻 视频| 最后的刺客免费高清国语| 亚洲丝袜综合中文字幕| 亚洲熟女精品中文字幕| 一二三四中文在线观看免费高清| 黄色欧美视频在线观看| 这个男人来自地球电影免费观看 | 免费观看的影片在线观看| videos熟女内射| 在线看a的网站| 伦理电影免费视频| 亚洲不卡免费看| 全区人妻精品视频| 国产精品女同一区二区软件| 97在线人人人人妻| freevideosex欧美| 免费观看性生交大片5| 久久精品国产亚洲av天美| 制服丝袜香蕉在线| 青春草亚洲视频在线观看| 丝袜在线中文字幕| 国产成人精品久久久久久| 国产国拍精品亚洲av在线观看| 乱系列少妇在线播放| 中文资源天堂在线| 国产极品粉嫩免费观看在线 | 狂野欧美白嫩少妇大欣赏| 亚洲国产精品999| 少妇被粗大猛烈的视频| 免费观看的影片在线观看| 午夜免费观看性视频| 亚洲精品久久久久久婷婷小说| 久久97久久精品| 亚洲国产毛片av蜜桃av| 亚洲美女黄色视频免费看| 午夜福利视频精品| 成年女人在线观看亚洲视频| 国产精品成人在线| 99久国产av精品国产电影| 另类亚洲欧美激情| 欧美精品一区二区免费开放| 综合色丁香网| 日韩视频在线欧美| 美女脱内裤让男人舔精品视频| 精品国产一区二区久久| 日韩中文字幕视频在线看片| 日日摸夜夜添夜夜爱| 在线观看人妻少妇| 成年美女黄网站色视频大全免费 | 国产亚洲午夜精品一区二区久久| 大香蕉97超碰在线| h视频一区二区三区| 午夜福利在线观看免费完整高清在| 最近最新中文字幕免费大全7| 午夜激情福利司机影院| 欧美xxⅹ黑人| 亚洲成人手机| 精品人妻熟女毛片av久久网站| 亚洲精品色激情综合| 嘟嘟电影网在线观看| 制服丝袜香蕉在线| 国模一区二区三区四区视频| 黄色配什么色好看| 欧美性感艳星| 国产高清三级在线| 视频区图区小说| 最近2019中文字幕mv第一页| 成年美女黄网站色视频大全免费 | 欧美 日韩 精品 国产| 天堂中文最新版在线下载| 欧美日韩视频高清一区二区三区二| 99热这里只有是精品50| 国产欧美日韩综合在线一区二区 | 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频| 色5月婷婷丁香| 久久久久视频综合| 亚洲精品乱码久久久v下载方式| 亚洲精品视频女| 最近手机中文字幕大全| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产成人久久av| 丁香六月天网| 免费播放大片免费观看视频在线观看| 成人二区视频| 国内少妇人妻偷人精品xxx网站| 久久国产亚洲av麻豆专区| 99久久精品一区二区三区| 亚洲精品,欧美精品| 亚洲精品久久午夜乱码| 狂野欧美白嫩少妇大欣赏| 新久久久久国产一级毛片| 国产精品福利在线免费观看| 黄色毛片三级朝国网站 | av.在线天堂| 天美传媒精品一区二区| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| 能在线免费看毛片的网站| 秋霞伦理黄片| 观看av在线不卡| 偷拍熟女少妇极品色| 99久久人妻综合| 天堂中文最新版在线下载| 日本猛色少妇xxxxx猛交久久| 欧美xxxx性猛交bbbb| 十八禁网站网址无遮挡 | 日韩电影二区| 婷婷色麻豆天堂久久| 91精品国产国语对白视频| 亚洲av在线观看美女高潮| av在线观看视频网站免费| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 成人黄色视频免费在线看| 亚洲婷婷狠狠爱综合网| 国产精品国产av在线观看| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 国产欧美日韩精品一区二区| 日本黄大片高清| 成人无遮挡网站| 日韩一区二区视频免费看| 美女国产视频在线观看| 久久精品国产a三级三级三级| 国产真实伦视频高清在线观看| 日韩 亚洲 欧美在线| 成年美女黄网站色视频大全免费 | 丰满迷人的少妇在线观看| 特大巨黑吊av在线直播| 街头女战士在线观看网站| 毛片一级片免费看久久久久| 亚洲丝袜综合中文字幕| 少妇的逼好多水| 18禁动态无遮挡网站| 久久人妻熟女aⅴ| 一级毛片电影观看| 免费人成在线观看视频色| 最新的欧美精品一区二区| 国产免费视频播放在线视频| 97精品久久久久久久久久精品| 99久久精品一区二区三区| 国产黄色免费在线视频| av女优亚洲男人天堂| 中国美白少妇内射xxxbb| 大片电影免费在线观看免费| 晚上一个人看的免费电影| 国产免费福利视频在线观看| 国语对白做爰xxxⅹ性视频网站| kizo精华| 黑人高潮一二区| 97在线人人人人妻| 国产精品成人在线| 韩国高清视频一区二区三区| 交换朋友夫妻互换小说| 国产精品一区www在线观看| 在线天堂最新版资源| 成人特级av手机在线观看| 久久狼人影院| 少妇丰满av| 丝袜在线中文字幕| 国产成人精品无人区| 午夜免费鲁丝| 亚洲久久久国产精品| 亚洲av男天堂| 一级av片app| 国产精品成人在线| 欧美日韩av久久| 国产成人freesex在线| 久久午夜综合久久蜜桃| 亚洲av成人精品一区久久| 免费观看a级毛片全部| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 精品亚洲成a人片在线观看| 久久热精品热| 丝袜脚勾引网站| 成人免费观看视频高清| 欧美另类一区| 午夜视频国产福利| 国产一区亚洲一区在线观看| 亚洲国产欧美日韩在线播放 | 国产日韩欧美视频二区| 黄色一级大片看看| 日本爱情动作片www.在线观看| 国产色爽女视频免费观看| 日产精品乱码卡一卡2卡三| 国产成人aa在线观看| 成人亚洲欧美一区二区av| 国产高清国产精品国产三级| 狠狠精品人妻久久久久久综合| a级毛片在线看网站| 亚洲欧洲精品一区二区精品久久久 | 国产日韩欧美在线精品| 你懂的网址亚洲精品在线观看| 一本—道久久a久久精品蜜桃钙片| 久久午夜福利片| 久久精品久久久久久久性| 另类精品久久| 国产午夜精品一二区理论片| av免费观看日本| 色视频在线一区二区三区| 中文字幕制服av| 国产黄片美女视频| 国产亚洲午夜精品一区二区久久| 国产精品蜜桃在线观看| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 久久久久久伊人网av| 插阴视频在线观看视频| 熟女av电影| xxx大片免费视频| 久久久久久久大尺度免费视频| 久久国产乱子免费精品| a 毛片基地| 十八禁网站网址无遮挡 | 夜夜骑夜夜射夜夜干| av有码第一页| 青春草亚洲视频在线观看| av在线app专区| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 日韩av在线免费看完整版不卡| 最近的中文字幕免费完整| 全区人妻精品视频| 午夜福利视频精品| 欧美日韩一区二区视频在线观看视频在线| av在线播放精品| 日韩欧美一区视频在线观看 | 三级国产精品片| 亚洲欧美日韩东京热| 欧美少妇被猛烈插入视频| 国产精品蜜桃在线观看| 我的老师免费观看完整版| 男男h啪啪无遮挡| h视频一区二区三区| 日韩大片免费观看网站| 亚洲av成人精品一区久久| 视频区图区小说| 精品国产国语对白av| 亚洲av二区三区四区| 一级毛片久久久久久久久女| 久久这里有精品视频免费| 日韩强制内射视频| 日韩视频在线欧美| 边亲边吃奶的免费视频| 久久影院123| 久久久久久久精品精品| 99久久精品国产国产毛片| 精品一区二区三卡| 深夜a级毛片| 成人影院久久| 青春草国产在线视频| 国产欧美亚洲国产| 2022亚洲国产成人精品| 亚洲怡红院男人天堂| a级毛片免费高清观看在线播放| 热99国产精品久久久久久7| 另类亚洲欧美激情| 国产黄频视频在线观看| 多毛熟女@视频| 丝袜脚勾引网站| 久久久国产欧美日韩av| 在线天堂最新版资源| 多毛熟女@视频| 日本黄色片子视频| 免费人成在线观看视频色| 男女边摸边吃奶| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品日本国产第一区| 大又大粗又爽又黄少妇毛片口| 亚洲av日韩在线播放| 亚洲综合色惰| 久久97久久精品| 国产无遮挡羞羞视频在线观看| 亚洲自偷自拍三级| 偷拍熟女少妇极品色| 国产黄片视频在线免费观看| 亚洲欧美成人精品一区二区| 日日爽夜夜爽网站| 精品少妇黑人巨大在线播放| 成人无遮挡网站| 不卡视频在线观看欧美| 嫩草影院入口| 国产成人午夜福利电影在线观看| 我要看日韩黄色一级片| 99久久人妻综合| 国产熟女欧美一区二区| 妹子高潮喷水视频| 这个男人来自地球电影免费观看 | 午夜激情久久久久久久| 精品亚洲成国产av| 一本—道久久a久久精品蜜桃钙片| 亚洲丝袜综合中文字幕| 97超视频在线观看视频| kizo精华| 男男h啪啪无遮挡| 国产伦精品一区二区三区四那| 91久久精品国产一区二区三区| 日韩成人伦理影院| 国产精品女同一区二区软件| 久久久久久久久久久免费av| 久久国产亚洲av麻豆专区| 免费观看在线日韩| 在现免费观看毛片| 亚洲欧美精品专区久久| 色视频www国产| 99热国产这里只有精品6| 韩国高清视频一区二区三区| 极品人妻少妇av视频| 免费黄网站久久成人精品| 亚洲一级一片aⅴ在线观看| 乱码一卡2卡4卡精品| 欧美另类一区| 亚洲丝袜综合中文字幕| av在线播放精品| 久久久久国产网址| 我要看黄色一级片免费的| 波野结衣二区三区在线| 日韩成人av中文字幕在线观看| 亚洲av二区三区四区| 天堂8中文在线网| av有码第一页| 在线观看免费高清a一片| 免费观看在线日韩| 国产片特级美女逼逼视频| 亚洲美女黄色视频免费看| 99热全是精品| 亚洲美女视频黄频| 成年av动漫网址| 国产亚洲欧美精品永久| 中文资源天堂在线| 久久精品久久精品一区二区三区| 亚洲电影在线观看av| 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| av福利片在线| 超碰97精品在线观看| 中国三级夫妇交换| 亚洲三级黄色毛片| 嫩草影院新地址| 国产成人freesex在线| 久久国内精品自在自线图片| 边亲边吃奶的免费视频| 亚洲四区av| 国产一区二区三区av在线| 亚洲国产欧美日韩在线播放 | 特大巨黑吊av在线直播| 99热国产这里只有精品6| 大又大粗又爽又黄少妇毛片口| 最近最新中文字幕免费大全7| 涩涩av久久男人的天堂| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 久久国内精品自在自线图片| 我的老师免费观看完整版| 下体分泌物呈黄色| 亚洲精品日韩在线中文字幕| 国产乱来视频区| 久久久久久伊人网av| 免费黄频网站在线观看国产| av又黄又爽大尺度在线免费看| 国产精品不卡视频一区二区| av一本久久久久| 亚洲人成网站在线观看播放| 大片免费播放器 马上看| 高清午夜精品一区二区三区| 免费黄色在线免费观看| freevideosex欧美| 老司机亚洲免费影院| 亚洲欧美精品专区久久| 国产91av在线免费观看| 欧美国产精品一级二级三级 | 国产在线视频一区二区| 国产亚洲精品久久久com| av天堂中文字幕网| 国模一区二区三区四区视频| 国产亚洲5aaaaa淫片| 天堂俺去俺来也www色官网| 美女国产视频在线观看| 中文字幕久久专区| 久久人人爽人人爽人人片va| av有码第一页| 亚洲精品国产av成人精品| 中文欧美无线码| 99热这里只有精品一区| 精品久久国产蜜桃| 日韩在线高清观看一区二区三区| 一级毛片 在线播放| 午夜免费观看性视频| av国产精品久久久久影院| 性高湖久久久久久久久免费观看| 久久久久国产精品人妻一区二区| 精品人妻偷拍中文字幕| 99久久人妻综合| 亚洲欧洲国产日韩| 日韩在线高清观看一区二区三区| 亚洲精品国产av成人精品| 成人毛片a级毛片在线播放| 纯流量卡能插随身wifi吗| 国产极品粉嫩免费观看在线 | 免费看av在线观看网站| 黄片无遮挡物在线观看| 日韩视频在线欧美| 国产视频首页在线观看| 99精国产麻豆久久婷婷| 人妻人人澡人人爽人人| 十八禁网站网址无遮挡 | 看十八女毛片水多多多| 最近的中文字幕免费完整| 久久久久久久久大av| 日日摸夜夜添夜夜添av毛片| 丰满少妇做爰视频| 国产 一区精品| 男人添女人高潮全过程视频| 啦啦啦视频在线资源免费观看| 三级经典国产精品| 久久99热6这里只有精品| 日日啪夜夜撸| 在线观看一区二区三区激情| 噜噜噜噜噜久久久久久91| 国产熟女欧美一区二区| 欧美xxxx性猛交bbbb| 欧美日韩视频精品一区| 亚洲精品乱码久久久v下载方式| 色婷婷av一区二区三区视频| 欧美激情国产日韩精品一区| 免费观看无遮挡的男女| 伦理电影大哥的女人| 久久6这里有精品| 亚洲第一区二区三区不卡| 中文字幕免费在线视频6| 国产成人免费观看mmmm| 日韩,欧美,国产一区二区三区| 国语对白做爰xxxⅹ性视频网站| 最新的欧美精品一区二区| 老熟女久久久| 草草在线视频免费看| 中文资源天堂在线| 国产精品久久久久久av不卡| 中文欧美无线码|