• <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的吸嘴支撐臂有限元分析
    亚洲国产精品久久男人天堂| 欧美日韩精品网址| 51午夜福利影视在线观看| 亚洲激情在线av| 欧美日本亚洲视频在线播放| 久久久久久免费高清国产稀缺| 香蕉久久夜色| 亚洲精品av麻豆狂野| 男女那种视频在线观看| 99久久综合精品五月天人人| 少妇的丰满在线观看| 久9热在线精品视频| 两个人视频免费观看高清| 淫秽高清视频在线观看| 88av欧美| 国产精品免费一区二区三区在线| 国产av一区二区精品久久| 999精品在线视频| 在线观看66精品国产| 熟妇人妻久久中文字幕3abv| 亚洲国产日韩欧美精品在线观看 | 超碰成人久久| 看片在线看免费视频| xxx96com| 欧美国产日韩亚洲一区| 免费在线观看成人毛片| 在线观看免费日韩欧美大片| 日本 欧美在线| 天堂动漫精品| 草草在线视频免费看| 国产精品久久久久久亚洲av鲁大| 午夜两性在线视频| 精品福利观看| 亚洲成av人片在线播放无| 中文字幕久久专区| √禁漫天堂资源中文www| 国产精品亚洲av一区麻豆| 欧美日韩福利视频一区二区| 最近最新免费中文字幕在线| 欧美在线一区亚洲| 欧美成人性av电影在线观看| 欧美性猛交╳xxx乱大交人| 久久香蕉国产精品| 免费av毛片视频| 亚洲精品久久成人aⅴ小说| 精品人妻1区二区| 一区二区三区高清视频在线| 国产99白浆流出| 国产视频一区二区在线看| 激情在线观看视频在线高清| 岛国在线观看网站| 欧美日韩精品网址| 成人高潮视频无遮挡免费网站| 在线十欧美十亚洲十日本专区| 免费在线观看视频国产中文字幕亚洲| 国产一级毛片七仙女欲春2| 久久 成人 亚洲| 成人一区二区视频在线观看| 99精品在免费线老司机午夜| 久久人妻av系列| 91老司机精品| 成人欧美大片| 国产探花在线观看一区二区| 欧美一级a爱片免费观看看 | 午夜福利在线在线| 神马国产精品三级电影在线观看 | 一本久久中文字幕| 18禁裸乳无遮挡免费网站照片| 最近最新免费中文字幕在线| 男人的好看免费观看在线视频 | 国产精品一区二区精品视频观看| 婷婷亚洲欧美| 老熟妇乱子伦视频在线观看| 免费无遮挡裸体视频| 国产亚洲av高清不卡| 999久久久精品免费观看国产| 国产精品久久久久久久电影 | 午夜精品久久久久久毛片777| 精品熟女少妇八av免费久了| 国产熟女午夜一区二区三区| 热99re8久久精品国产| 99热这里只有精品一区 | 草草在线视频免费看| 美女黄网站色视频| 99在线人妻在线中文字幕| АⅤ资源中文在线天堂| 午夜福利在线在线| 中文字幕人成人乱码亚洲影| 两个人免费观看高清视频| 久久精品国产清高在天天线| 国产成年人精品一区二区| 亚洲中文av在线| 一个人免费在线观看电影 | 午夜福利视频1000在线观看| 最近在线观看免费完整版| 无人区码免费观看不卡| 久久久久久亚洲精品国产蜜桃av| 亚洲男人的天堂狠狠| 男女那种视频在线观看| 国产精品野战在线观看| 欧美日本视频| 午夜激情av网站| 国产精品99久久99久久久不卡| 成人av一区二区三区在线看| 欧美+亚洲+日韩+国产| 天天躁狠狠躁夜夜躁狠狠躁| 国产主播在线观看一区二区| 久久久久久久久久黄片| 亚洲av电影在线进入| 欧美zozozo另类| 国产亚洲av嫩草精品影院| 非洲黑人性xxxx精品又粗又长| 精品久久蜜臀av无| 一本精品99久久精品77| 可以在线观看毛片的网站| 熟妇人妻久久中文字幕3abv| 老熟妇乱子伦视频在线观看| 全区人妻精品视频| 免费看a级黄色片| 亚洲片人在线观看| 国产高清视频在线观看网站| 国产黄a三级三级三级人| 国内毛片毛片毛片毛片毛片| 日韩高清综合在线| bbb黄色大片| 久久久久国产一级毛片高清牌| 国内久久婷婷六月综合欲色啪| 国产99久久九九免费精品| 国产午夜精品论理片| 小说图片视频综合网站| 又粗又爽又猛毛片免费看| 巨乳人妻的诱惑在线观看| 亚洲av美国av| 床上黄色一级片| 欧美日韩乱码在线| 日韩高清综合在线| 夜夜爽天天搞| 在线观看免费午夜福利视频| 亚洲免费av在线视频| 高清在线国产一区| 国产亚洲精品av在线| 成人特级黄色片久久久久久久| 国产三级在线视频| 岛国视频午夜一区免费看| 搡老妇女老女人老熟妇| 久久午夜亚洲精品久久| xxxwww97欧美| 久久人妻av系列| 精品高清国产在线一区| 天天一区二区日本电影三级| 在线a可以看的网站| 999久久久国产精品视频| 一级a爱片免费观看的视频| 亚洲精品久久成人aⅴ小说| 欧美黑人欧美精品刺激| 午夜福利免费观看在线| 色播亚洲综合网| 午夜福利在线观看吧| 国产精品美女特级片免费视频播放器 | 国产视频一区二区在线看| 男人舔奶头视频| 男人舔女人下体高潮全视频| 亚洲精品国产精品久久久不卡| 国产黄色小视频在线观看| 18美女黄网站色大片免费观看| 别揉我奶头~嗯~啊~动态视频| 99热6这里只有精品| 久久伊人香网站| 国产精品野战在线观看| 色老头精品视频在线观看| 少妇粗大呻吟视频| 中文字幕人成人乱码亚洲影| 99精品欧美一区二区三区四区| 国产精品自产拍在线观看55亚洲| 国产成人精品无人区| 午夜福利视频1000在线观看| 国产单亲对白刺激| 法律面前人人平等表现在哪些方面| 久久久精品大字幕| 亚洲最大成人中文| 午夜两性在线视频| www.自偷自拍.com| 草草在线视频免费看| 精品久久久久久久毛片微露脸| 国产伦人伦偷精品视频| 成人av在线播放网站| 在线看三级毛片| 免费看日本二区| 国内精品久久久久精免费| 国内精品一区二区在线观看| ponron亚洲| 亚洲国产精品sss在线观看| 欧美黑人欧美精品刺激| 久久草成人影院| 午夜日韩欧美国产| 老熟妇乱子伦视频在线观看| 少妇粗大呻吟视频| 啦啦啦免费观看视频1| 十八禁网站免费在线| 精品免费久久久久久久清纯| 在线观看日韩欧美| 欧美中文综合在线视频| 国产1区2区3区精品| 小说图片视频综合网站| 午夜免费成人在线视频| 亚洲激情在线av| 欧美性猛交╳xxx乱大交人| tocl精华| 国产v大片淫在线免费观看| 国产探花在线观看一区二区| 精品国产亚洲在线| 国内揄拍国产精品人妻在线| 久久久久久大精品| 1024手机看黄色片| 丰满人妻熟妇乱又伦精品不卡| 好看av亚洲va欧美ⅴa在| 亚洲一区二区三区色噜噜| 国内揄拍国产精品人妻在线| 精品久久久久久久毛片微露脸| 日本一二三区视频观看| 中出人妻视频一区二区| 一区二区三区国产精品乱码| 亚洲一区二区三区色噜噜| 久久精品国产清高在天天线| 69av精品久久久久久| 老司机午夜福利在线观看视频| 国产亚洲精品综合一区在线观看 | 最近最新中文字幕大全免费视频| 亚洲成人久久爱视频| 黑人欧美特级aaaaaa片| 色老头精品视频在线观看| 久久久久亚洲av毛片大全| 国产欧美日韩精品亚洲av| 国产成人精品久久二区二区91| 国产成人av教育| 九色国产91popny在线| 老熟妇乱子伦视频在线观看| av中文乱码字幕在线| 九色成人免费人妻av| 在线观看午夜福利视频| 91麻豆av在线| 中文字幕av在线有码专区| 亚洲国产高清在线一区二区三| 国产一区二区在线观看日韩 | 午夜老司机福利片| 观看免费一级毛片| 成人国产综合亚洲| 午夜免费成人在线视频| 麻豆国产97在线/欧美 | 日本三级黄在线观看| 韩国av一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 美女午夜性视频免费| 极品教师在线免费播放| 在线观看日韩欧美| 久久热在线av| 国产精品1区2区在线观看.| 日韩欧美国产一区二区入口| 久久久水蜜桃国产精品网| 亚洲精品久久国产高清桃花| 老司机午夜十八禁免费视频| 啦啦啦韩国在线观看视频| 久久香蕉精品热| 国产亚洲精品第一综合不卡| 一二三四社区在线视频社区8| 中亚洲国语对白在线视频| 国产日本99.免费观看| 99精品在免费线老司机午夜| 亚洲电影在线观看av| 国内揄拍国产精品人妻在线| 男女之事视频高清在线观看| x7x7x7水蜜桃| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品久久久久久毛片| 国产主播在线观看一区二区| 极品教师在线免费播放| 制服诱惑二区| 精品日产1卡2卡| 久久伊人香网站| 国产精品一区二区三区四区免费观看 | 色播亚洲综合网| 国产69精品久久久久777片 | 麻豆一二三区av精品| 动漫黄色视频在线观看| 草草在线视频免费看| 搡老妇女老女人老熟妇| 亚洲欧美激情综合另类| 中文字幕av在线有码专区| 妹子高潮喷水视频| 久久这里只有精品中国| 国产亚洲精品久久久久5区| 后天国语完整版免费观看| 天堂影院成人在线观看| 国产成+人综合+亚洲专区| 日韩精品青青久久久久久| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 欧美日本视频| 美女大奶头视频| 久久久久久久精品吃奶| 亚洲av电影不卡..在线观看| 又黄又爽又免费观看的视频| 亚洲国产看品久久| 久久久久久久精品吃奶| 蜜桃久久精品国产亚洲av| 亚洲av五月六月丁香网| 午夜激情av网站| 午夜福利欧美成人| 人成视频在线观看免费观看| 成人三级做爰电影| 曰老女人黄片| 一区福利在线观看| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 欧美精品亚洲一区二区| 他把我摸到了高潮在线观看| 亚洲av第一区精品v没综合| 精品欧美一区二区三区在线| 搞女人的毛片| 美女大奶头视频| 51午夜福利影视在线观看| 校园春色视频在线观看| 欧美乱码精品一区二区三区| 欧美又色又爽又黄视频| 露出奶头的视频| 国产69精品久久久久777片 | 91在线观看av| 天天躁夜夜躁狠狠躁躁| 国产精品日韩av在线免费观看| 很黄的视频免费| 看黄色毛片网站| 午夜福利在线在线| 制服丝袜大香蕉在线| 国产精品亚洲av一区麻豆| 少妇熟女aⅴ在线视频| 国产精品久久久av美女十八| 国产高清有码在线观看视频 | 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 又大又爽又粗| 人妻丰满熟妇av一区二区三区| 国产精品综合久久久久久久免费| 欧美一区二区国产精品久久精品 | 日本熟妇午夜| 国产蜜桃级精品一区二区三区| 久久午夜亚洲精品久久| 9191精品国产免费久久| 女人高潮潮喷娇喘18禁视频| 亚洲aⅴ乱码一区二区在线播放 | 国产av一区二区精品久久| tocl精华| 久久久久免费精品人妻一区二区| 村上凉子中文字幕在线| 香蕉丝袜av| 99热这里只有精品一区 | 成人欧美大片| 精品久久久久久,| 精品久久久久久久末码| 欧美一级a爱片免费观看看 | 亚洲人与动物交配视频| 欧美日韩国产亚洲二区| 婷婷六月久久综合丁香| 久久午夜亚洲精品久久| 看黄色毛片网站| 久久亚洲精品不卡| 精品熟女少妇八av免费久了| 欧美乱色亚洲激情| 久久精品国产清高在天天线| 一本久久中文字幕| 精品久久久久久久末码| 女生性感内裤真人,穿戴方法视频| av片东京热男人的天堂| 窝窝影院91人妻| 两个人的视频大全免费| 18禁美女被吸乳视频| 久久午夜综合久久蜜桃| 男人舔女人的私密视频| 国产成人精品久久二区二区91| 欧美人与性动交α欧美精品济南到| 国产午夜精品论理片| 亚洲av成人不卡在线观看播放网| 免费电影在线观看免费观看| 亚洲七黄色美女视频| 国产精品亚洲av一区麻豆| 日本成人三级电影网站| 精华霜和精华液先用哪个| 日本免费a在线| 国产69精品久久久久777片 | 人人妻,人人澡人人爽秒播| 亚洲av美国av| 麻豆av在线久日| www.熟女人妻精品国产| 日韩精品青青久久久久久| 五月玫瑰六月丁香| xxxwww97欧美| 男男h啪啪无遮挡| 欧美日韩精品网址| 男女做爰动态图高潮gif福利片| 亚洲五月天丁香| 婷婷丁香在线五月| 少妇熟女aⅴ在线视频| 精品久久久久久,| 日本一二三区视频观看| 久久人妻av系列| 亚洲成人久久性| 一个人免费在线观看电影 | 啪啪无遮挡十八禁网站| 欧美丝袜亚洲另类 | avwww免费| 亚洲成人久久爱视频| 老熟妇乱子伦视频在线观看| www日本黄色视频网| 两个人视频免费观看高清| 脱女人内裤的视频| 亚洲精品在线观看二区| 18美女黄网站色大片免费观看| 欧美久久黑人一区二区| 一级a爱片免费观看的视频| 一本一本综合久久| 国产成人av激情在线播放| 亚洲国产精品999在线| 国产成人精品久久二区二区91| 欧美3d第一页| 丁香六月欧美| 日韩欧美 国产精品| 51午夜福利影视在线观看| 婷婷丁香在线五月| 国产成+人综合+亚洲专区| 一进一出好大好爽视频| 免费在线观看完整版高清| 精品久久久久久成人av| 男人舔奶头视频| 久久精品国产亚洲av高清一级| 精品久久久久久久久久免费视频| 人妻夜夜爽99麻豆av| 五月玫瑰六月丁香| 国产v大片淫在线免费观看| 免费无遮挡裸体视频| 99re在线观看精品视频| 精品熟女少妇八av免费久了| 欧洲精品卡2卡3卡4卡5卡区| 女人被狂操c到高潮| 99精品久久久久人妻精品| 亚洲国产精品成人综合色| 一级a爱片免费观看的视频| 久久精品人妻少妇| 国产成+人综合+亚洲专区| 欧美成人一区二区免费高清观看 | 国产野战对白在线观看| 亚洲狠狠婷婷综合久久图片| 国产视频内射| 欧美绝顶高潮抽搐喷水| 国产精品免费一区二区三区在线| 在线观看免费日韩欧美大片| 久久久精品大字幕| 日韩 欧美 亚洲 中文字幕| 一本久久中文字幕| 搞女人的毛片| 婷婷六月久久综合丁香| 中国美女看黄片| 黑人操中国人逼视频| 在线视频色国产色| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 久久欧美精品欧美久久欧美| 黄片大片在线免费观看| 窝窝影院91人妻| 亚洲精品在线观看二区| 在线观看午夜福利视频| 久久 成人 亚洲| 国产91精品成人一区二区三区| 男插女下体视频免费在线播放| 欧美日韩瑟瑟在线播放| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利18| 黄色片一级片一级黄色片| 成人永久免费在线观看视频| 国产不卡一卡二| 久久久精品欧美日韩精品| 校园春色视频在线观看| 人人妻人人澡欧美一区二区| 日本 欧美在线| 欧美性长视频在线观看| 草草在线视频免费看| 亚洲狠狠婷婷综合久久图片| 久久久久久久久免费视频了| 午夜福利在线在线| 久久久国产成人免费| 男男h啪啪无遮挡| 草草在线视频免费看| 亚洲第一欧美日韩一区二区三区| 亚洲欧美精品综合一区二区三区| 国产精品,欧美在线| 99热这里只有精品一区 | 亚洲成a人片在线一区二区| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 国产69精品久久久久777片 | 可以在线观看毛片的网站| 好男人电影高清在线观看| 嫁个100分男人电影在线观看| 国产亚洲精品综合一区在线观看 | av视频在线观看入口| 黑人操中国人逼视频| 国产高清激情床上av| 少妇的丰满在线观看| 久久久久久免费高清国产稀缺| www.www免费av| 男男h啪啪无遮挡| 草草在线视频免费看| 校园春色视频在线观看| 黄色 视频免费看| 十八禁人妻一区二区| 欧美乱码精品一区二区三区| 国产av一区二区精品久久| 校园春色视频在线观看| 日韩免费av在线播放| 少妇裸体淫交视频免费看高清 | 国产69精品久久久久777片 | 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 色在线成人网| 亚洲色图 男人天堂 中文字幕| 日韩国内少妇激情av| 88av欧美| 久久精品国产99精品国产亚洲性色| 亚洲专区中文字幕在线| 久久中文看片网| 伦理电影免费视频| 精华霜和精华液先用哪个| 桃红色精品国产亚洲av| 久久久精品大字幕| 美女免费视频网站| 俺也久久电影网| 熟女电影av网| 黄片小视频在线播放| 一本精品99久久精品77| 国产三级中文精品| 久99久视频精品免费| 国产精品98久久久久久宅男小说| 欧美成人性av电影在线观看| 国产片内射在线| 亚洲全国av大片| 一级作爱视频免费观看| 欧美精品啪啪一区二区三区| 制服丝袜大香蕉在线| 欧美久久黑人一区二区| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 在线a可以看的网站| 18美女黄网站色大片免费观看| 欧美日本视频| 十八禁人妻一区二区| 91九色精品人成在线观看| 欧美日韩乱码在线| www日本在线高清视频| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 日韩 欧美 亚洲 中文字幕| 色av中文字幕| 亚洲全国av大片| 三级毛片av免费| 国产亚洲av高清不卡| 十八禁网站免费在线| 欧美性猛交╳xxx乱大交人| 亚洲精品国产一区二区精华液| 亚洲av片天天在线观看| 在线免费观看的www视频| 国产精品亚洲av一区麻豆| 国产精品亚洲一级av第二区| 18禁裸乳无遮挡免费网站照片| 亚洲自偷自拍图片 自拍| 99国产综合亚洲精品| 啦啦啦韩国在线观看视频| 午夜福利免费观看在线| 欧美另类亚洲清纯唯美| 一进一出抽搐gif免费好疼| 欧美日韩一级在线毛片| 一级黄色大片毛片| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦免费观看视频1| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 欧美一级毛片孕妇| 男人舔奶头视频| 亚洲欧美精品综合久久99| 欧美午夜高清在线| 精品久久蜜臀av无| 亚洲片人在线观看| 亚洲天堂国产精品一区在线| 老汉色av国产亚洲站长工具| 成人精品一区二区免费| 亚洲精品色激情综合| 高潮久久久久久久久久久不卡| 国产一区二区在线观看日韩 | 欧美最黄视频在线播放免费| 精华霜和精华液先用哪个| 99热这里只有精品一区 | 90打野战视频偷拍视频| 亚洲国产欧美网| 波多野结衣高清作品| 成人精品一区二区免费| 欧美黑人欧美精品刺激| 啦啦啦韩国在线观看视频| 亚洲成人国产一区在线观看| 中文在线观看免费www的网站 | 妹子高潮喷水视频| 亚洲欧美日韩高清专用| 在线观看免费视频日本深夜| 最近在线观看免费完整版| 亚洲专区字幕在线| 69av精品久久久久久|