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

    基于有限元數(shù)值方法的金屬材料內(nèi)壁缺陷超聲頻譜特性研究

    2017-01-03 07:33:31孫繼華趙揚南鋼洋馬健劉帥巨陽
    山東科學(xué) 2016年5期
    關(guān)鍵詞:鋁板內(nèi)壁頻域

    孫繼華,趙揚,南鋼洋,馬健,劉帥,巨陽

    (山東省科學(xué)院激光研究所,山東省無損檢測工程技術(shù)研究中心,山東 濟南 250103)

    ?

    【新材料】

    基于有限元數(shù)值方法的金屬材料內(nèi)壁缺陷超聲頻譜特性研究

    孫繼華,趙揚,南鋼洋,馬健,劉帥,巨陽

    (山東省科學(xué)院激光研究所,山東省無損檢測工程技術(shù)研究中心,山東 濟南 250103)

    根據(jù)金屬材料內(nèi)壁缺陷的特點,建立了吸收邊界和自由邊界相結(jié)合的模擬模型,對含有不同寬度內(nèi)壁凹陷形平底槽缺陷的模型進(jìn)行了超聲頻譜特性模擬研究。通過對波場快照、時域信號和頻譜分析得到了不同尺寸平底槽對聲波傳播頻譜特性的影響以及中心頻率、帶寬、半波高寬度和峰值等呈規(guī)律性變化的頻域特征量。研究結(jié)果對超聲檢測頻譜分析技術(shù)的理論和實驗提供了參考。

    金屬材料;超聲檢測;模擬;內(nèi)壁缺陷;頻域特征量

    隨著現(xiàn)代超聲檢測技術(shù)的發(fā)展,模擬研究越來越多地應(yīng)用到相關(guān)領(lǐng)域中。板材、管道和壓力容器等金屬材料的內(nèi)壁由于腐蝕、外力以及摩擦等作用,易形成凹陷形的內(nèi)壁缺陷,這種缺陷對金屬材料的性能影響極大,不及時發(fā)現(xiàn)就會帶來安全隱患。超聲無損檢測具有無損、快捷、成本低及對人體無害等特點[1],是對該類型缺陷較為可行、有效的檢測方法,但傳統(tǒng)超聲檢測方法通常在時域中對檢測結(jié)果進(jìn)行定性分析,而超聲檢測信號中通常攜帶了大量的結(jié)構(gòu)信息,不同結(jié)構(gòu)的時頻特性不同[2],所以目前對檢測信號的分析正在從時域法向頻域法發(fā)展,頻域分析法為超聲無損檢測技術(shù)注入了活力,一定程度上起到了變革作用[3-4]。通過對超聲檢測信號進(jìn)行頻譜分析可以獲得缺陷在頻域中的各種特征,進(jìn)而對缺陷進(jìn)行進(jìn)一步的識別。

    數(shù)值模擬方法可以應(yīng)用于超聲無損檢測技術(shù)的研究,是理論和實驗重要的驗證手段,其中有限元模擬方法理論基礎(chǔ)比較成熟,具有通用性和靈活性,在模擬聲場方面具有獨特的優(yōu)點,可以正確地模擬金屬材料中超聲傳播情況,是一種有效的數(shù)值模擬工具[5-6]。國內(nèi)外學(xué)者近十幾年來應(yīng)用有限元數(shù)值模擬方法對超聲檢測進(jìn)行了一系列的研究[7-11], 但針對金屬材料內(nèi)壁缺陷超聲頻譜特性方面的研究還比較少。為了深入地認(rèn)識超聲波在有凹陷形的內(nèi)壁缺陷金屬中的傳播情況,本文采用有限元數(shù)值模擬方法研究超聲波在金屬材料中的傳播,通過對所得信號進(jìn)行頻譜分析,得到了不同尺寸內(nèi)壁平底槽對超聲波頻譜特性的影響,模擬結(jié)果可以為超聲檢測頻譜分析技術(shù)的理論和實驗研究提供參考。

    1 建立有限元數(shù)值模型

    1.1 金屬材料中超聲檢測有限元數(shù)值模擬原理[12]

    在各向同性彈性固體中,聲波的Navier控制方程的矢量形式為:

    (1)

    式中,λ和μ為Lamé常數(shù),ρ為介質(zhì)密度,u為位移。

    按照有限元數(shù)值方法計算步驟,首先選取合適網(wǎng)格對計算區(qū)域離散化,并在單元內(nèi)選擇合適的插值函數(shù)得到單元結(jié)點位移、速度和加速度;利用伽遼金(Galerkin)法,由式(1)可導(dǎo)出整體有限元運動方程:

    (2)

    1.2 建立模型

    為了模擬金屬材料中的超聲波傳播,選用在民用、軍工等行業(yè)比較常用的鋁板為母材模型,鋁板的縱波速度vP為6 200 m/s。本文從有限元數(shù)值模擬方法基本理論出發(fā),采用邊長為0.1 mm的三角形單元對一個長22 mm,厚度為16 mm的鋁板材料進(jìn)行有限元剖分[13],建立的有限元模擬模型如圖1所示。圖1中鋁板的上表面和下表面采用自由邊界用于模擬鋁-空氣界面,兩側(cè)采用吸收邊界用于模擬鋁板無限延伸[11]。在模型底部設(shè)置了3種深度均為0.8 mm,寬度則分別為4、5、6 mm的平底槽缺陷,用于模擬材料的內(nèi)壁缺陷。

    圖1 有限元模型Fig.1 Finite element model

    為了使模擬結(jié)果與實際超聲檢測時Φ10 mm超聲探頭產(chǎn)生的平面縱波聲場一致,本文在模擬模型上表面中心10 mm區(qū)域內(nèi)加載高斯分布的應(yīng)力邊界條件以激發(fā)出脈沖超聲縱波[14],以模擬超聲波激勵源,該激勵源同時可以作為接收源,用于模擬接收超聲波信號并以A-掃描形式顯示出來。

    2 模擬結(jié)果與分析

    2.1 波場快照圖分析

    模擬過程中超聲波由上表面激發(fā)出向下傳播,對超聲波分別傳播到0.51、2.58、3.39、4.38 μs時刻的全波場y方向位移波場快照進(jìn)行分析,不同時刻波場快照如圖2所示,從圖上能明顯地看出模擬激發(fā)的超聲波在鋁板中的傳播情況。

    圖2 不同時刻波場快照圖Fig.2 Wave field snapshots at different moments

    本文模擬的是垂直入射的超聲縱波,由圖2可以清晰地看到在0.51 μs時,模擬激發(fā)的超聲縱波P已經(jīng)開始向下傳播,在2.58 μs時超聲波已經(jīng)傳播到底面的平底槽位置,在3.39 μs和4.38 μs時可以看到超聲波碰到平底槽表面和底面后發(fā)生的反射和散射現(xiàn)象。從圖上可以清晰地觀察到很多類型的超聲波,其中信號比較明顯的波形有反射縱波R-P和模式轉(zhuǎn)換的橫波M-S等,此后時間里這些不同類型的超聲波會按自己的速度和方向繼續(xù)傳播。從圖2上也可以看出,本文建立的有限元模型中加了吸收邊界條件的兩個側(cè)面沒有明顯的超聲反射波,說明模型加的吸收邊界條件效果很好,達(dá)到了模擬鋁板的無限延伸的目的。所以通過以上波場快照圖分析,直觀地顯示了超聲波在帶有底面凹槽鋁中傳播的情況,模擬的結(jié)果與實際相符。

    2.2 時域信號分析

    為了獲得時域A-掃描信號,我們在上表面激發(fā)出超聲波的區(qū)域內(nèi)對y方向位移進(jìn)行掃描,可以得到圖3所示的不同缺陷尺寸下A-掃描曲線圖,其中圖3a是底面沒有缺陷時的A-掃描曲線圖,圖3b、c、d是底面有不同尺寸的平底槽時的A-掃描曲線圖。

    從鋁板上表面激發(fā)出的超聲波向下傳播過程中遇到鋁板底面和平底槽時就會發(fā)生反射和散射,從圖3a可以看出上表面激發(fā)的始波(P)和底面回波(R-P)信號,這和實際的超聲檢測儀顯示的A-掃描形式一致。為了確定獲得信號的正確性,我們把圖3a中無缺陷底面回波R-P的接收時間t=5.115 μs,代入公式d=tvP/2,可以計算出鋁板厚度為15.9 mm,與實際模型厚度16 mm誤差僅為0.6%,采用同樣方法計算出的各種平底槽深度也與實際模型數(shù)值相符,說明模擬的超聲波傳播過程是和實際相符的。

    2.3 頻譜分析

    為了獲得超聲檢測信號頻域特征,對底面回波信號進(jìn)行頻譜分析。如圖3b中所示,我們選取合適長度的時間閘門分別截取圖3a、b、c、d中A-掃描曲線上的底面回波,經(jīng)頻譜分析得到幅度譜圖如圖4所示,同時計算得到頻域特征量如表1所示。

    由圖4可以看出,當(dāng)鋁板底面有不同尺寸的平底槽時接收到的回波幅度比沒有平底槽時的幅度大,而且不同尺寸平底槽的回波幅度也不同。由表1可以看出,當(dāng)鋁板底面有寬度為4 mm和5 mm的平底槽時,頻域特征量中的中心頻率和峰值頻率比沒有平底槽時要大;而當(dāng)有寬度為6 mm的平底槽時,其中心頻率和峰值頻率比沒有平底槽時要小。由表1還可以看到當(dāng)有不同寬度的平底槽時,頻域特征量中的半高寬度、帶寬等特征量比沒有平底槽時都要小,這也反映了平底槽對不同頻域特征量的影響。

    圖3 A-掃描曲線Fig.3 The A-scanning curve

    圖4 回波幅度譜Fig.4 Amplitude spectrum of echo

    表1 頻域特征量Table 1 Frequency domain characteristics

    為了得到鋁板中不同尺寸的平底槽對超聲頻域特征量的影響,我們把平底槽寬度d以1 mm為間隔從1 mm增加到12 mm分別進(jìn)行模擬,通過對結(jié)果進(jìn)行計算分析,可以得到峰值、中心頻率、半高寬度和帶寬等的變化趨勢如圖5所示。

    圖5 頻域特征量變化圖Fig.5 Variation chart of frequency domain characteristics

    從圖5可以看出,峰值變化的曲線是隨著平底槽寬度的增加首先呈上升趨勢,當(dāng)平底槽寬度在6 mm左右時出現(xiàn)極大值,隨后就開始變小。頻域特征量中心頻率在平底槽寬度為2 mm左右時有一個極大值,然后開始變??;在8 mm左右出現(xiàn)最小值,此后又開始變大。半高寬度變化曲線在平底槽寬度為4 mm左右有極小值,另外在11 mm左右有極大值,隨后開始變小。帶寬和半高寬度的變化趨勢基本一致,只是實際數(shù)值比半高寬度的要小。

    通過對頻域特征量變化規(guī)律的分析可知,隨著鋁板中平底槽尺寸的變化,會引起各種頻域特征量規(guī)律性的變化,這可以為缺陷識別提供一定的參考。

    3 結(jié)論

    本文通過有限元數(shù)值方法模擬計算了超聲波在有平底槽缺陷鋁板中的傳播情況,通過波場快照、時域A-掃描圖和頻譜等分析了鋁板中不同尺寸內(nèi)壁平底槽缺陷對超聲傳播的影響,主要得出以下結(jié)論:

    (1)通過模擬過程中的波場快照圖,能清晰地觀察到超聲波在鋁板中的傳播以及遇到底面和平底槽時的反射、散射和模式轉(zhuǎn)換等各種現(xiàn)象,模擬結(jié)果與實際情況相符。

    (2)由A-掃描圖上無缺陷底面回波時間計算出的模擬區(qū)域厚度與模型厚度一致,說明采用有限元數(shù)值模擬方法可以正確地模擬超聲在鋁板中的傳播情況。

    (3)通過對鋁板不同尺寸平底槽回波進(jìn)行頻譜分析,得到各種頻域特征量峰值、中心頻率、半高寬和帶寬等呈規(guī)律性的變化,可以為缺陷識別提供一定的參考。

    通過以上分析可知,采用有限元數(shù)值方法模擬超聲波在金屬材料中的傳播是可行的,通過頻譜分析獲得的結(jié)果可以為實際中的內(nèi)壁凹陷形缺陷識別提供依據(jù),并為超聲檢測頻譜分析技術(shù)的理論和實驗研究提供參考。但是本文只對金屬內(nèi)壁較規(guī)則的平底槽進(jìn)行了模擬計算,下一步將開展針對金屬內(nèi)壁不同形狀缺陷的模擬計算研究,以期獲得不同類型缺陷對超聲頻譜特性的影響。

    [1]齊子誠,鄭子軍,唐盛明,等.滾動式鋼板超聲檢測設(shè)備研制[J].兵器材料科學(xué)與工程,2016,39(2):76-79.

    [2]艾春安,蔡笑風(fēng),劉繼方,等.短時傅里葉變換聲-超聲檢測信號缺陷識別[J].中國測試,2015,41(4):29-31.

    [3] 劉立群,趙振南.頻譜分析技術(shù)在超聲檢測中的應(yīng)用[J].一重技術(shù),2013(5):57-61.

    [4]蔣志峰,吳作倫,吳瑞明,等.材料缺陷超聲檢測中譜分析技術(shù)的研究[J].兵器材料科學(xué)與工程,2009,32(4):117-120.

    [5]FELICE M V, VELICHKO A, WILCOX P D,e tal. Obtaining geometries of real cracks and using an efficient finite element method to simulate their ultrasonic array response[J].Insight: Non-Destructive Testing and Condition Monitoring, 2014, 56(9): 492-498 .

    [6] 趙軍輝,魏勤,袁媛,等. Lamb波檢測板中裂紋的有限元模擬[J].壓電與聲光,2013,35(3):320-324.

    [7]XUE T,LORD W,UDPA S. Numerical analysis of the radiated fields of ultrasonic transducers [J]. Journal of Nondestructive Evaluation, 1995,14(3): 137-146.

    [8]DATTA D,KISHORE N N. Features of ultrasonic wave propagation to identify defects in composite materials modelled by finite element method [J]. NDT & E International,1996,29(4): 213-223.

    [9]KISHORE N N,SRIDHAR I,IYENGR N G R.Finite element modeling of the scattering ultrasonic waves by isolated flaws[J]. NDT & E International,2000, 33(5):297-305.

    [10]BLANLOEUIL P, MEZIANE A, BACON C. 2D finite element modeling of the non-collinear mixing method for detection and characterization of closed cracks [J]. NDT & E International,2015, 76: 43-51.

    [11]孫朝明.超聲橫波檢測的有限元模擬計算與分析[J].壓電與聲光,2014,36(6):1018-1021.

    [12] 李樹榜,李書光,劉學(xué)鋒.裂紋超聲散射的有限元模擬[J].無損檢測,2007,29(1):27-31.

    [13]LIU T,LIU K,ZHANG J. Triangular grid method for stress-wave propagation in 2-D orthotropic materials[J]. Archive of Applied Mechanics, 2005, 74(7):477-488.

    [14] 韓秀甫,滕永平,吳迪,等.超聲換能器暫態(tài)聲場的理論分析與計算[J].安慶師范學(xué)院學(xué)報(自然科學(xué)版),2004,10(3):25-28.

    Finite element based research on ultrasonic spectral characteristic of inner wall defect of metal materials

    SUN Ji-hua, ZHAO Yang, NAN Gang-yang, MA Jian, LIU Shuai,JU Yang

    (Shandong Provincial Research Center of Non-destructive Testing Engineering Technology,Laser Institute, Shandong Academy of Sciences, Jinan 250103, China)

    ∶We construct an absorption boundary and free boundary combined model for inner wall defect of metal materials. We further simulate ultrasonic spectral characteristic of concave groove of different width inner wall. We also acquire the impact of different sized groove on spectral characteristic of ultrasonic propagation by wave field snapshots, time-domain signal and spectral analysis, and such regularly variant frequency domain characteristics as central frequency, bandwidth, half width and peaks. The results provide a reference for the theory and experiment of ultrasonic testing spectral analysis.

    ∶ metal materials; ultrasonic testing; simulation; inner wall defect; frequency domain characteristic quantity

    10.3976/j.issn.1002-4026.2016.05.012

    2016-05-05

    山東省科學(xué)院杰出青年計劃( 201505);濟南市高校院所自主創(chuàng)新基金( 201401220,201401209)

    孫繼華(1979—),男,助理研究員,研究方向為材料的無損檢測與評價研究。

    TP391.9;TB52+6

    A

    猜你喜歡
    鋁板內(nèi)壁頻域
    哪塊邊角料的面積大
    大型鋁板拉伸機液壓底座的設(shè)計計算
    垣曲北白鵝墓地出土的青銅匽姬甗(M3:10)
    文物季刊(2021年1期)2021-03-23 08:12:58
    膀胱內(nèi)壁子宮內(nèi)膜異位1例
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設(shè)計
    雙曲弧形純鋁板內(nèi)幕墻的施工探討
    智能城市(2018年7期)2018-07-10 08:30:24
    集流管內(nèi)壁沖壓模具設(shè)計
    長脈沖激光與連續(xù)激光對鋁板熱破壞仿真對比
    核電反應(yīng)堆壓力容器頂蓋J型接頭內(nèi)壁殘余應(yīng)力
    焊接(2016年1期)2016-02-27 12:54:45
    基于改進(jìn)Radon-Wigner變換的目標(biāo)和拖曳式誘餌頻域分離
    最新中文字幕久久久久| 亚洲av电影在线观看一区二区三区| 亚洲综合色网址| 亚洲内射少妇av| 亚洲一区二区三区欧美精品| 欧美+日韩+精品| 国产精品久久久av美女十八| 久久久久久久大尺度免费视频| 国产亚洲精品久久久com| 亚洲色图综合在线观看| 亚洲精品av麻豆狂野| 看免费av毛片| 天堂8中文在线网| 亚洲av电影在线观看一区二区三区| 亚洲情色 制服丝袜| 伦理电影免费视频| 精品一区在线观看国产| 精品一区在线观看国产| 99久久中文字幕三级久久日本| 精品一区二区免费观看| 午夜av观看不卡| 最黄视频免费看| av福利片在线| 中文欧美无线码| av黄色大香蕉| 男人舔女人的私密视频| 纯流量卡能插随身wifi吗| 人妻少妇偷人精品九色| 欧美人与性动交α欧美软件 | 乱码一卡2卡4卡精品| 丝袜喷水一区| 韩国精品一区二区三区 | 午夜av观看不卡| 制服诱惑二区| 国产精品一国产av| 香蕉国产在线看| 一区二区三区精品91| 精品人妻偷拍中文字幕| 精品人妻偷拍中文字幕| 欧美日韩一区二区视频在线观看视频在线| 捣出白浆h1v1| 免费在线观看完整版高清| 亚洲性久久影院| 宅男免费午夜| 亚洲,欧美,日韩| 国产精品人妻久久久久久| 午夜免费观看性视频| 日韩视频在线欧美| 在线观看一区二区三区激情| 国产xxxxx性猛交| 国产淫语在线视频| 美女中出高潮动态图| 亚洲欧美成人精品一区二区| 免费看光身美女| 99视频精品全部免费 在线| 国产免费一级a男人的天堂| 天天躁夜夜躁狠狠久久av| 亚洲精品,欧美精品| 欧美日韩精品成人综合77777| 成人黄色视频免费在线看| 极品人妻少妇av视频| 久久久久精品久久久久真实原创| 亚洲天堂av无毛| 哪个播放器可以免费观看大片| 欧美日韩精品成人综合77777| 十分钟在线观看高清视频www| 国产精品秋霞免费鲁丝片| 日韩制服骚丝袜av| 国产亚洲欧美精品永久| 国产高清国产精品国产三级| 久久韩国三级中文字幕| 亚洲欧美清纯卡通| 少妇 在线观看| 女的被弄到高潮叫床怎么办| 久久久a久久爽久久v久久| 久久精品aⅴ一区二区三区四区 | 水蜜桃什么品种好| 丝袜在线中文字幕| 久久精品国产a三级三级三级| 久久精品国产鲁丝片午夜精品| 午夜影院在线不卡| 另类亚洲欧美激情| 亚洲中文av在线| 国产女主播在线喷水免费视频网站| av.在线天堂| 男的添女的下面高潮视频| 久久99精品国语久久久| 制服诱惑二区| 国产精品秋霞免费鲁丝片| 日韩中文字幕视频在线看片| 精品久久国产蜜桃| 成人影院久久| 日日爽夜夜爽网站| 色婷婷久久久亚洲欧美| 人妻系列 视频| 久久精品国产鲁丝片午夜精品| 久久久精品免费免费高清| 两个人免费观看高清视频| 国产毛片在线视频| 久久国产精品大桥未久av| 母亲3免费完整高清在线观看 | 亚洲成人av在线免费| 国产成人精品久久久久久| 成人手机av| 桃花免费在线播放| 国产亚洲欧美精品永久| 中国国产av一级| 咕卡用的链子| 水蜜桃什么品种好| 黄色 视频免费看| 国产色爽女视频免费观看| 两个人免费观看高清视频| 成人亚洲精品一区在线观看| 永久网站在线| 日产精品乱码卡一卡2卡三| 亚洲天堂av无毛| 777米奇影视久久| 天堂8中文在线网| 欧美 亚洲 国产 日韩一| 99热全是精品| 晚上一个人看的免费电影| 免费黄频网站在线观看国产| 人妻 亚洲 视频| 女人久久www免费人成看片| 在线天堂最新版资源| 精品一区二区三区视频在线| 另类精品久久| 精品一品国产午夜福利视频| 日韩av在线免费看完整版不卡| 桃花免费在线播放| 亚洲国产最新在线播放| 精品99又大又爽又粗少妇毛片| 亚洲国产看品久久| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 国产亚洲欧美精品永久| videossex国产| 97在线人人人人妻| 夫妻性生交免费视频一级片| 免费播放大片免费观看视频在线观看| 亚洲精品久久午夜乱码| 狂野欧美激情性bbbbbb| 亚洲少妇的诱惑av| 蜜桃在线观看..| 久久狼人影院| 大香蕉久久网| 午夜福利,免费看| 一级片免费观看大全| 中文精品一卡2卡3卡4更新| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 欧美日本中文国产一区发布| 日本av免费视频播放| 在线观看国产h片| 久久影院123| av在线播放精品| 日日摸夜夜添夜夜爱| 水蜜桃什么品种好| 秋霞伦理黄片| 国产福利在线免费观看视频| xxx大片免费视频| 一级毛片我不卡| 国产欧美日韩一区二区三区在线| 国产在线视频一区二区| 日本午夜av视频| 一区二区三区四区激情视频| 婷婷色av中文字幕| 午夜福利,免费看| 制服诱惑二区| 国产精品女同一区二区软件| 五月天丁香电影| 国产又爽黄色视频| 欧美xxⅹ黑人| 99热网站在线观看| av卡一久久| 亚洲国产成人一精品久久久| 国产淫语在线视频| 亚洲精品aⅴ在线观看| 九色成人免费人妻av| 最后的刺客免费高清国语| av线在线观看网站| 亚洲人成77777在线视频| 日韩精品有码人妻一区| 高清不卡的av网站| 99热全是精品| 少妇精品久久久久久久| 国产高清不卡午夜福利| 精品熟女少妇av免费看| 日韩中字成人| 久久午夜综合久久蜜桃| av不卡在线播放| 亚洲欧美中文字幕日韩二区| 一级爰片在线观看| 免费观看在线日韩| 18禁裸乳无遮挡动漫免费视频| 亚洲精品久久午夜乱码| 免费黄色在线免费观看| 巨乳人妻的诱惑在线观看| 免费av不卡在线播放| 亚洲精品自拍成人| 一区二区av电影网| 久久久久久久精品精品| 国产老妇伦熟女老妇高清| 免费黄频网站在线观看国产| 久久女婷五月综合色啪小说| 最黄视频免费看| 免费观看性生交大片5| 欧美丝袜亚洲另类| 国产毛片在线视频| 亚洲av国产av综合av卡| 国产国拍精品亚洲av在线观看| 亚洲精品日本国产第一区| 久久国产精品男人的天堂亚洲 | 婷婷色麻豆天堂久久| 国产成人精品无人区| 狠狠婷婷综合久久久久久88av| 成人手机av| 人人妻人人添人人爽欧美一区卜| 国产精品国产三级国产专区5o| 涩涩av久久男人的天堂| 亚洲av电影在线进入| 午夜av观看不卡| 亚洲精品美女久久久久99蜜臀 | 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 成年av动漫网址| 国产熟女欧美一区二区| 久久久久视频综合| av电影中文网址| 日韩精品免费视频一区二区三区 | 观看美女的网站| 热re99久久精品国产66热6| 狠狠精品人妻久久久久久综合| 欧美激情国产日韩精品一区| 亚洲色图 男人天堂 中文字幕 | 美女视频免费永久观看网站| 免费在线观看完整版高清| 亚洲美女黄色视频免费看| 国产亚洲最大av| 少妇熟女欧美另类| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久精品94久久精品| 自拍欧美九色日韩亚洲蝌蚪91| 免费观看性生交大片5| 美女视频免费永久观看网站| 性高湖久久久久久久久免费观看| 亚洲美女视频黄频| 久久午夜福利片| 美女中出高潮动态图| 亚洲av国产av综合av卡| 久久99热这里只频精品6学生| 欧美3d第一页| 欧美日本中文国产一区发布| 成人黄色视频免费在线看| 综合色丁香网| 蜜桃国产av成人99| 纯流量卡能插随身wifi吗| 男女午夜视频在线观看 | 少妇人妻久久综合中文| 一二三四中文在线观看免费高清| 热99国产精品久久久久久7| 在线亚洲精品国产二区图片欧美| 久久毛片免费看一区二区三区| 久久久久国产网址| 亚洲国产精品一区二区三区在线| 日韩av不卡免费在线播放| 精品国产乱码久久久久久小说| 又粗又硬又长又爽又黄的视频| 久久综合国产亚洲精品| 日韩欧美精品免费久久| 一级黄片播放器| 99久久人妻综合| 亚洲人成网站在线观看播放| 日韩熟女老妇一区二区性免费视频| 国产免费福利视频在线观看| 国产一区二区三区av在线| 七月丁香在线播放| 青春草亚洲视频在线观看| 久久午夜福利片| 免费黄色在线免费观看| 韩国精品一区二区三区 | www.色视频.com| 男女国产视频网站| 美国免费a级毛片| 啦啦啦啦在线视频资源| videos熟女内射| 中国美白少妇内射xxxbb| 精品国产一区二区三区久久久樱花| 91午夜精品亚洲一区二区三区| 99热国产这里只有精品6| 日本猛色少妇xxxxx猛交久久| 午夜福利视频精品| 99热6这里只有精品| 欧美日韩亚洲高清精品| 秋霞伦理黄片| 亚洲av电影在线进入| 蜜臀久久99精品久久宅男| 菩萨蛮人人尽说江南好唐韦庄| 丰满乱子伦码专区| 欧美人与性动交α欧美软件 | 亚洲欧美精品自产自拍| 成人黄色视频免费在线看| 伦理电影免费视频| 天堂俺去俺来也www色官网| 捣出白浆h1v1| 精品国产乱码久久久久久小说| av国产久精品久网站免费入址| 亚洲精品一区蜜桃| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲国产日韩| 国产1区2区3区精品| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 黄片无遮挡物在线观看| 在线观看国产h片| 天堂俺去俺来也www色官网| 熟女av电影| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 国产爽快片一区二区三区| 男女免费视频国产| 有码 亚洲区| 成人无遮挡网站| 人人妻人人添人人爽欧美一区卜| 久久精品国产鲁丝片午夜精品| 伦精品一区二区三区| 国产精品国产三级国产av玫瑰| 在线免费观看不下载黄p国产| 国产精品.久久久| 中文字幕最新亚洲高清| 亚洲少妇的诱惑av| 国产 一区精品| 欧美日韩av久久| 欧美最新免费一区二区三区| 国产精品免费大片| 久久99热6这里只有精品| 最近中文字幕2019免费版| 成人影院久久| 丝袜喷水一区| 日韩精品有码人妻一区| 99热这里只有是精品在线观看| 亚洲av日韩在线播放| 一二三四中文在线观看免费高清| kizo精华| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产av蜜桃| 亚洲精品视频女| av电影中文网址| 精品少妇黑人巨大在线播放| 我的女老师完整版在线观看| 久久人人爽人人片av| 精品国产一区二区三区久久久樱花| 亚洲情色 制服丝袜| 久久精品国产a三级三级三级| 日韩av在线免费看完整版不卡| 久久午夜综合久久蜜桃| 欧美国产精品va在线观看不卡| 精品久久蜜臀av无| 色哟哟·www| 国产精品99久久99久久久不卡 | 成人手机av| 欧美国产精品va在线观看不卡| 女性生殖器流出的白浆| 亚洲av电影在线观看一区二区三区| 考比视频在线观看| 中文字幕人妻熟女乱码| 午夜福利视频精品| 精品一品国产午夜福利视频| 少妇被粗大猛烈的视频| 五月天丁香电影| 亚洲国产精品999| 菩萨蛮人人尽说江南好唐韦庄| 久久午夜福利片| av在线app专区| 大陆偷拍与自拍| 中文字幕人妻熟女乱码| 亚洲,一卡二卡三卡| 999精品在线视频| 国产在线视频一区二区| 欧美日韩精品成人综合77777| 日本爱情动作片www.在线观看| 欧美亚洲日本最大视频资源| 午夜精品国产一区二区电影| 国产av一区二区精品久久| 亚洲成色77777| 男女啪啪激烈高潮av片| 免费久久久久久久精品成人欧美视频 | 国产成人精品无人区| 午夜激情久久久久久久| 国语对白做爰xxxⅹ性视频网站| 欧美成人午夜免费资源| 色94色欧美一区二区| 亚洲av.av天堂| 大香蕉97超碰在线| 久久国产精品男人的天堂亚洲 | 久久午夜综合久久蜜桃| 国产乱来视频区| 少妇的逼水好多| 18在线观看网站| 久久ye,这里只有精品| 亚洲精品第二区| 在线精品无人区一区二区三| 99香蕉大伊视频| 国产1区2区3区精品| 国产色爽女视频免费观看| 欧美精品一区二区大全| 亚洲 欧美一区二区三区| 青春草视频在线免费观看| 亚洲色图 男人天堂 中文字幕 | 亚洲人成网站在线观看播放| 日日撸夜夜添| 热re99久久精品国产66热6| 大话2 男鬼变身卡| 国产免费视频播放在线视频| 少妇被粗大的猛进出69影院 | 美女国产高潮福利片在线看| 久久久久精品人妻al黑| 久久人人爽人人爽人人片va| 亚洲精品中文字幕在线视频| 亚洲人成77777在线视频| www.熟女人妻精品国产 | 国产永久视频网站| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频在线播放观看不卡| 久久久欧美国产精品| av播播在线观看一区| 高清不卡的av网站| av黄色大香蕉| 亚洲欧美一区二区三区黑人 | 日韩精品有码人妻一区| 国产 精品1| 国产精品不卡视频一区二区| 精品一区二区三区四区五区乱码 | 人人妻人人添人人爽欧美一区卜| 97在线人人人人妻| 天堂中文最新版在线下载| 制服人妻中文乱码| 乱码一卡2卡4卡精品| 国产精品成人在线| 9热在线视频观看99| 黑丝袜美女国产一区| 亚洲精品国产色婷婷电影| 乱人伦中国视频| 亚洲精品自拍成人| 99久久综合免费| 欧美+日韩+精品| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的| 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 在线观看www视频免费| 国产免费一级a男人的天堂| 啦啦啦中文免费视频观看日本| 日韩一区二区视频免费看| av免费观看日本| 51国产日韩欧美| 永久网站在线| 久久国产亚洲av麻豆专区| 亚洲精华国产精华液的使用体验| 久久久久视频综合| 国产精品久久久久久久久免| 午夜老司机福利剧场| 大码成人一级视频| 亚洲四区av| 免费大片黄手机在线观看| 欧美精品国产亚洲| 18禁观看日本| 日韩一本色道免费dvd| av卡一久久| 国产探花极品一区二区| 日韩av免费高清视频| 18禁观看日本| 1024视频免费在线观看| 18禁在线无遮挡免费观看视频| a级片在线免费高清观看视频| 国产男人的电影天堂91| 免费黄频网站在线观看国产| 亚洲国产av影院在线观看| 久久久久久久国产电影| 久久人人爽人人片av| 久久精品人人爽人人爽视色| av天堂久久9| 日本91视频免费播放| 久久99蜜桃精品久久| 侵犯人妻中文字幕一二三四区| 天天操日日干夜夜撸| 一二三四中文在线观看免费高清| 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www| 亚洲精品国产av蜜桃| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片| 精品国产乱码久久久久久小说| 18禁动态无遮挡网站| 国产精品一国产av| 99热全是精品| 1024视频免费在线观看| 黄色视频在线播放观看不卡| 久久99一区二区三区| 黄色毛片三级朝国网站| 国产精品人妻久久久久久| 日韩伦理黄色片| 最后的刺客免费高清国语| 国产成人一区二区在线| 国产激情久久老熟女| 欧美 日韩 精品 国产| 中文天堂在线官网| 久久久久精品性色| 国产女主播在线喷水免费视频网站| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| 久久综合国产亚洲精品| 少妇的逼水好多| 国产亚洲午夜精品一区二区久久| 久久精品熟女亚洲av麻豆精品| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 亚洲精品自拍成人| 99久国产av精品国产电影| 99视频精品全部免费 在线| 国产成人免费无遮挡视频| 尾随美女入室| 久久影院123| 考比视频在线观看| 久久久国产一区二区| 国产又爽黄色视频| 日日撸夜夜添| 自线自在国产av| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 我的女老师完整版在线观看| 丝袜喷水一区| 精品国产一区二区久久| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 黄网站色视频无遮挡免费观看| 午夜日本视频在线| 免费日韩欧美在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产色片| 亚洲欧美精品自产自拍| 久久 成人 亚洲| 国产探花极品一区二区| 国国产精品蜜臀av免费| 最后的刺客免费高清国语| 国产一区二区激情短视频 | 精品酒店卫生间| 久久国产亚洲av麻豆专区| 国产在线免费精品| 亚洲欧美日韩卡通动漫| 色5月婷婷丁香| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 色哟哟·www| 秋霞伦理黄片| 亚洲精品国产色婷婷电影| 十八禁网站网址无遮挡| 午夜福利视频精品| 天堂中文最新版在线下载| 亚洲国产精品成人久久小说| 男女免费视频国产| 精品久久久精品久久久| 久久久久久久久久人人人人人人| 男男h啪啪无遮挡| 免费高清在线观看日韩| 欧美人与性动交α欧美软件 | 亚洲精品av麻豆狂野| 另类亚洲欧美激情| 免费av中文字幕在线| 考比视频在线观看| 欧美少妇被猛烈插入视频| 精品人妻在线不人妻| 国产成人精品婷婷| 丝袜喷水一区| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 中文字幕人妻熟女乱码| 自线自在国产av| 少妇的丰满在线观看| 国产乱人偷精品视频| 七月丁香在线播放| 欧美精品高潮呻吟av久久| 亚洲一级一片aⅴ在线观看| 又黄又爽又刺激的免费视频.| 各种免费的搞黄视频| 国产日韩欧美在线精品| 啦啦啦在线观看免费高清www| 国产色爽女视频免费观看| 一本色道久久久久久精品综合| 亚洲欧美色中文字幕在线| 国产男人的电影天堂91| 午夜视频国产福利| 免费大片18禁| 久久久国产精品麻豆| 久久国内精品自在自线图片| 亚洲av欧美aⅴ国产| 最近的中文字幕免费完整| 亚洲精品视频女| 激情视频va一区二区三区| 亚洲熟女精品中文字幕| 久久久久久人人人人人| 一区在线观看完整版| 三上悠亚av全集在线观看| 精品一区二区免费观看| 国产高清三级在线| 欧美另类一区| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 99国产综合亚洲精品| 欧美 亚洲 国产 日韩一| 亚洲中文av在线| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说|