• <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)和拖曳式誘餌頻域分離
    亚洲国产精品合色在线| 波多野结衣高清无吗| 色老头精品视频在线观看| 国产成+人综合+亚洲专区| 757午夜福利合集在线观看| 国产视频一区二区在线看| 亚洲精品一卡2卡三卡4卡5卡| 婷婷亚洲欧美| 脱女人内裤的视频| 啦啦啦观看免费观看视频高清| 免费一级毛片在线播放高清视频| 成人欧美大片| 91九色精品人成在线观看| 搡老熟女国产l中国老女人| 成人永久免费在线观看视频| 日本熟妇午夜| 亚洲色图av天堂| 熟女电影av网| 一区二区三区高清视频在线| 亚洲中文字幕一区二区三区有码在线看| 精品一区二区三区视频在线观看免费| ponron亚洲| 亚洲精品一区av在线观看| 欧美在线黄色| 99久久精品热视频| 美女黄网站色视频| 国产亚洲欧美在线一区二区| 国内毛片毛片毛片毛片毛片| 岛国在线免费视频观看| 91久久精品国产一区二区成人 | 国产在视频线在精品| 又黄又爽又免费观看的视频| 九九久久精品国产亚洲av麻豆| 香蕉丝袜av| 亚洲成a人片在线一区二区| 一区二区三区国产精品乱码| 老熟妇仑乱视频hdxx| 在线观看av片永久免费下载| 美女大奶头视频| 中文字幕人妻丝袜一区二区| 久久性视频一级片| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区色噜噜| 宅男免费午夜| 高清毛片免费观看视频网站| 国产色婷婷99| 中文字幕av成人在线电影| 床上黄色一级片| 真人一进一出gif抽搐免费| 真实男女啪啪啪动态图| 一区二区三区高清视频在线| 亚洲性夜色夜夜综合| 午夜激情福利司机影院| 欧美乱色亚洲激情| 日本 av在线| 免费搜索国产男女视频| 午夜久久久久精精品| 激情在线观看视频在线高清| 香蕉av资源在线| 天美传媒精品一区二区| 性色avwww在线观看| 亚洲七黄色美女视频| 99riav亚洲国产免费| 美女 人体艺术 gogo| 欧美一级毛片孕妇| 最近最新中文字幕大全电影3| 18+在线观看网站| 国内精品一区二区在线观看| 国产精品一区二区三区四区免费观看 | 法律面前人人平等表现在哪些方面| 在线视频色国产色| 丰满人妻一区二区三区视频av | 最新在线观看一区二区三区| 国产精品 欧美亚洲| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 日韩人妻高清精品专区| 他把我摸到了高潮在线观看| 国产高清有码在线观看视频| 欧美日韩福利视频一区二区| 午夜福利成人在线免费观看| 国产激情偷乱视频一区二区| 在线十欧美十亚洲十日本专区| 婷婷精品国产亚洲av| 少妇熟女aⅴ在线视频| 中文字幕人妻丝袜一区二区| 91在线精品国自产拍蜜月 | 午夜日韩欧美国产| 欧美三级亚洲精品| 搡老岳熟女国产| 亚洲成人中文字幕在线播放| 午夜亚洲福利在线播放| 99久久精品国产亚洲精品| 少妇的逼水好多| 动漫黄色视频在线观看| 日韩欧美国产一区二区入口| 午夜福利在线在线| 亚洲国产精品sss在线观看| 2021天堂中文幕一二区在线观| 欧美日韩黄片免| 99久久久亚洲精品蜜臀av| 18禁黄网站禁片免费观看直播| 国产黄片美女视频| 欧美精品啪啪一区二区三区| bbb黄色大片| 精品无人区乱码1区二区| 欧美黄色淫秽网站| xxx96com| 天堂影院成人在线观看| 九色国产91popny在线| 成人一区二区视频在线观看| 色综合婷婷激情| 99久久精品一区二区三区| 免费av观看视频| 熟妇人妻久久中文字幕3abv| 在线观看66精品国产| 国产午夜精品久久久久久一区二区三区 | 亚洲av五月六月丁香网| 国产三级中文精品| 美女大奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩av在线大香蕉| 久久久色成人| 国产精品 欧美亚洲| 好看av亚洲va欧美ⅴa在| 国产午夜福利久久久久久| 村上凉子中文字幕在线| 亚洲成人中文字幕在线播放| 草草在线视频免费看| 女人高潮潮喷娇喘18禁视频| 欧美bdsm另类| 国产精品三级大全| 天天躁日日操中文字幕| 久久性视频一级片| 中出人妻视频一区二区| 国产av一区在线观看免费| 真实男女啪啪啪动态图| 黄色片一级片一级黄色片| 欧美中文综合在线视频| 中文字幕av在线有码专区| 婷婷精品国产亚洲av| 亚洲精品日韩av片在线观看 | 女人被狂操c到高潮| 国产色婷婷99| 欧美精品啪啪一区二区三区| 国产高潮美女av| 午夜免费观看网址| 两个人视频免费观看高清| 免费在线观看亚洲国产| 亚洲无线观看免费| 午夜福利在线在线| 观看美女的网站| 色综合欧美亚洲国产小说| 欧美一区二区精品小视频在线| 亚洲人成网站在线播| 国产免费男女视频| 97超级碰碰碰精品色视频在线观看| 欧美成人性av电影在线观看| 国产男靠女视频免费网站| 无限看片的www在线观看| 日韩av在线大香蕉| 亚洲av中文字字幕乱码综合| 久久亚洲真实| 中亚洲国语对白在线视频| 性色av乱码一区二区三区2| 亚洲av免费在线观看| 国产一区二区在线av高清观看| 九色成人免费人妻av| 少妇丰满av| 精品电影一区二区在线| 国产一区二区三区视频了| 日本黄大片高清| 变态另类成人亚洲欧美熟女| 1000部很黄的大片| 九九在线视频观看精品| 国产av在哪里看| 国产伦精品一区二区三区视频9 | 欧美性猛交╳xxx乱大交人| 亚洲狠狠婷婷综合久久图片| 中文字幕熟女人妻在线| 免费看a级黄色片| 国产男靠女视频免费网站| 午夜福利免费观看在线| 欧美黑人欧美精品刺激| 久久这里只有精品中国| 亚洲精品日韩av片在线观看 | 国产一区二区三区在线臀色熟女| 少妇的逼好多水| 毛片女人毛片| 日本在线视频免费播放| 人人妻人人看人人澡| 国产爱豆传媒在线观看| 亚洲精品一区av在线观看| 国产三级黄色录像| 午夜福利在线观看吧| 国产亚洲欧美98| 真人做人爱边吃奶动态| 夜夜看夜夜爽夜夜摸| 黄色女人牲交| 国产高清有码在线观看视频| 欧美3d第一页| 18+在线观看网站| 日韩免费av在线播放| 日韩欧美三级三区| 黄片大片在线免费观看| 国产久久久一区二区三区| 国内精品一区二区在线观看| 在线视频色国产色| 国产精品精品国产色婷婷| 一进一出抽搐动态| 国产毛片a区久久久久| 波多野结衣高清作品| 老司机午夜十八禁免费视频| 身体一侧抽搐| 亚洲国产精品成人综合色| 精品久久久久久久久久免费视频| 日本五十路高清| 亚洲五月天丁香| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 国产精品嫩草影院av在线观看 | 久久久久久久久大av| 亚洲av熟女| 很黄的视频免费| 精品久久久久久成人av| 欧美日韩亚洲国产一区二区在线观看| 精品熟女少妇八av免费久了| 88av欧美| ponron亚洲| 三级男女做爰猛烈吃奶摸视频| 国产在视频线在精品| 韩国av一区二区三区四区| 狂野欧美激情性xxxx| 国产成人啪精品午夜网站| 韩国av一区二区三区四区| 真人做人爱边吃奶动态| 国产男靠女视频免费网站| 免费观看人在逋| 蜜桃亚洲精品一区二区三区| 波多野结衣巨乳人妻| 国产真实伦视频高清在线观看 | 精品一区二区三区人妻视频| 国产午夜精品久久久久久一区二区三区 | 制服人妻中文乱码| 欧美绝顶高潮抽搐喷水| 我的老师免费观看完整版| 久久久成人免费电影| 欧美zozozo另类| 欧美丝袜亚洲另类 | 男女之事视频高清在线观看| 欧美国产日韩亚洲一区| 老熟妇仑乱视频hdxx| 舔av片在线| 99热这里只有精品一区| 在线观看av片永久免费下载| 欧美bdsm另类| 俺也久久电影网| 国产色婷婷99| 美女免费视频网站| 午夜老司机福利剧场| 国产激情欧美一区二区| 国产男靠女视频免费网站| 亚洲精华国产精华精| 日本一本二区三区精品| 一进一出抽搐gif免费好疼| 欧美日韩乱码在线| 黑人欧美特级aaaaaa片| 人人妻人人澡欧美一区二区| 欧美性猛交╳xxx乱大交人| 亚洲人与动物交配视频| 国产乱人伦免费视频| 99热只有精品国产| 色视频www国产| 亚洲av五月六月丁香网| www.999成人在线观看| 俺也久久电影网| 久久精品国产综合久久久| 中文字幕人成人乱码亚洲影| 亚洲最大成人手机在线| 中文资源天堂在线| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 夜夜看夜夜爽夜夜摸| 亚洲精品色激情综合| 在线观看日韩欧美| 激情在线观看视频在线高清| 欧美黑人欧美精品刺激| 亚洲欧美激情综合另类| 熟妇人妻久久中文字幕3abv| 丰满人妻一区二区三区视频av | 9191精品国产免费久久| 一级a爱片免费观看的视频| 亚洲精华国产精华精| 淫妇啪啪啪对白视频| 97碰自拍视频| 51午夜福利影视在线观看| 热99re8久久精品国产| 窝窝影院91人妻| 午夜福利在线在线| 欧美成人a在线观看| 久久久久久久亚洲中文字幕 | 人人妻,人人澡人人爽秒播| 观看美女的网站| 色视频www国产| 成人高潮视频无遮挡免费网站| 亚洲熟妇熟女久久| 三级毛片av免费| 婷婷亚洲欧美| 国产午夜精品久久久久久一区二区三区 | 午夜免费成人在线视频| 亚洲欧美日韩卡通动漫| 啦啦啦观看免费观看视频高清| 久久人人精品亚洲av| 露出奶头的视频| 老司机午夜福利在线观看视频| h日本视频在线播放| 欧美黑人欧美精品刺激| 亚洲欧美激情综合另类| 国产黄片美女视频| 天堂√8在线中文| 夜夜夜夜夜久久久久| 亚洲专区国产一区二区| 九色成人免费人妻av| www国产在线视频色| 九九在线视频观看精品| 中文字幕熟女人妻在线| 国产精品永久免费网站| 国产一级毛片七仙女欲春2| 女生性感内裤真人,穿戴方法视频| av片东京热男人的天堂| 成人高潮视频无遮挡免费网站| 亚洲美女视频黄频| 午夜免费激情av| 麻豆一二三区av精品| 97人妻精品一区二区三区麻豆| 色精品久久人妻99蜜桃| 国产探花在线观看一区二区| 神马国产精品三级电影在线观看| а√天堂www在线а√下载| 中文在线观看免费www的网站| 午夜福利高清视频| 欧美一区二区亚洲| 亚洲av电影不卡..在线观看| 99久久99久久久精品蜜桃| www.色视频.com| 国产精品嫩草影院av在线观看 | 色播亚洲综合网| 99国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 国产午夜精品论理片| 淫秽高清视频在线观看| 国产熟女xx| 国产精品精品国产色婷婷| 国产亚洲精品综合一区在线观看| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添小说| 19禁男女啪啪无遮挡网站| 日韩大尺度精品在线看网址| 亚洲精品粉嫩美女一区| 小蜜桃在线观看免费完整版高清| 最近视频中文字幕2019在线8| 日韩欧美在线乱码| 偷拍熟女少妇极品色| 国产亚洲精品综合一区在线观看| 国产精品女同一区二区软件 | 日日摸夜夜添夜夜添小说| 19禁男女啪啪无遮挡网站| 99精品在免费线老司机午夜| 欧美日韩中文字幕国产精品一区二区三区| 亚洲午夜理论影院| 欧美日韩乱码在线| 手机成人av网站| 国产精品久久久久久人妻精品电影| 成人国产综合亚洲| 中亚洲国语对白在线视频| 亚洲专区国产一区二区| 午夜福利视频1000在线观看| av在线蜜桃| 久9热在线精品视频| 12—13女人毛片做爰片一| 欧美日韩一级在线毛片| 精品电影一区二区在线| 99在线视频只有这里精品首页| 国产黄片美女视频| 久久久久亚洲av毛片大全| 色哟哟哟哟哟哟| 在线观看美女被高潮喷水网站 | 超碰av人人做人人爽久久 | 不卡一级毛片| 观看美女的网站| 男女那种视频在线观看| 女人高潮潮喷娇喘18禁视频| 免费看十八禁软件| 久久精品国产综合久久久| 国产精品久久久人人做人人爽| 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 18禁在线播放成人免费| 婷婷亚洲欧美| 久久精品影院6| 精品99又大又爽又粗少妇毛片 | 三级国产精品欧美在线观看| 亚洲人成电影免费在线| 午夜福利在线在线| 婷婷丁香在线五月| 午夜免费观看网址| 看免费av毛片| 国产精品精品国产色婷婷| 内射极品少妇av片p| 岛国在线免费视频观看| 又紧又爽又黄一区二区| 最近在线观看免费完整版| 精品电影一区二区在线| 成人鲁丝片一二三区免费| 3wmmmm亚洲av在线观看| 欧美大码av| 男女床上黄色一级片免费看| 麻豆一二三区av精品| 99国产综合亚洲精品| 亚洲最大成人中文| 性色avwww在线观看| 色哟哟哟哟哟哟| 国产美女午夜福利| 亚洲国产精品999在线| 国产精品久久电影中文字幕| 97超视频在线观看视频| 成人精品一区二区免费| 免费观看人在逋| 成年版毛片免费区| 国产精品久久视频播放| 内射极品少妇av片p| 国产精品日韩av在线免费观看| 亚洲av不卡在线观看| 久久久久亚洲av毛片大全| 青草久久国产| 国产毛片a区久久久久| 国产一区二区三区在线臀色熟女| 午夜精品久久久久久毛片777| 午夜福利欧美成人| 18禁裸乳无遮挡免费网站照片| 日日摸夜夜添夜夜添小说| 久久精品国产自在天天线| 国产97色在线日韩免费| 国产一区二区三区视频了| 老鸭窝网址在线观看| 国产免费男女视频| 90打野战视频偷拍视频| 国产探花在线观看一区二区| 深夜精品福利| 性色av乱码一区二区三区2| a级一级毛片免费在线观看| 美女黄网站色视频| 制服人妻中文乱码| 国产蜜桃级精品一区二区三区| 99久国产av精品| 午夜福利高清视频| 久久精品91无色码中文字幕| 亚洲国产日韩欧美精品在线观看 | 三级毛片av免费| 在线观看免费视频日本深夜| 午夜福利视频1000在线观看| 欧美又色又爽又黄视频| 色av中文字幕| 国产精品香港三级国产av潘金莲| 夜夜夜夜夜久久久久| 久久久国产成人精品二区| 色老头精品视频在线观看| 免费一级毛片在线播放高清视频| 国产一级毛片七仙女欲春2| 成人无遮挡网站| 国产精品国产高清国产av| 日韩欧美免费精品| 首页视频小说图片口味搜索| 又爽又黄无遮挡网站| 999久久久精品免费观看国产| 午夜免费激情av| 日韩中文字幕欧美一区二区| 国模一区二区三区四区视频| 成年版毛片免费区| 色播亚洲综合网| 久久久久国产精品人妻aⅴ院| 欧美性感艳星| 欧美高清成人免费视频www| 国产黄a三级三级三级人| 一级毛片女人18水好多| АⅤ资源中文在线天堂| 男女那种视频在线观看| 国产私拍福利视频在线观看| 午夜久久久久精精品| 亚洲国产欧洲综合997久久,| 国产精品野战在线观看| 国产一区二区三区在线臀色熟女| 成人精品一区二区免费| 国产黄a三级三级三级人| 中文在线观看免费www的网站| 香蕉av资源在线| 日韩欧美在线二视频| 国产高清激情床上av| 91麻豆av在线| 久久国产乱子伦精品免费另类| 欧美黄色淫秽网站| 色在线成人网| 欧美3d第一页| 国产真实伦视频高清在线观看 | 91字幕亚洲| 国产真实乱freesex| 日韩欧美 国产精品| 欧美日韩一级在线毛片| 午夜福利高清视频| 在线免费观看不下载黄p国产 | 丰满的人妻完整版| 免费看日本二区| 亚洲成a人片在线一区二区| 日本五十路高清| 国产单亲对白刺激| 婷婷六月久久综合丁香| 最新美女视频免费是黄的| 国产av不卡久久| a级一级毛片免费在线观看| 亚洲真实伦在线观看| 欧美日韩福利视频一区二区| 1000部很黄的大片| 热99re8久久精品国产| 国产野战对白在线观看| 女人被狂操c到高潮| 午夜a级毛片| 天堂影院成人在线观看| 99精品欧美一区二区三区四区| 国产伦精品一区二区三区四那| 18禁在线播放成人免费| 99国产精品一区二区三区| 欧美黑人巨大hd| 国产乱人伦免费视频| 久久久久久久精品吃奶| 熟妇人妻久久中文字幕3abv| 日本a在线网址| 久久精品国产亚洲av涩爱 | 每晚都被弄得嗷嗷叫到高潮| 国产精品美女特级片免费视频播放器| 亚洲成人久久性| 狂野欧美白嫩少妇大欣赏| 国产男靠女视频免费网站| 国产日本99.免费观看| 久久久久亚洲av毛片大全| 97超级碰碰碰精品色视频在线观看| 国产野战对白在线观看| 国产精品亚洲av一区麻豆| 精华霜和精华液先用哪个| 亚洲久久久久久中文字幕| 香蕉av资源在线| 香蕉丝袜av| 成人精品一区二区免费| 国产精品 欧美亚洲| 国产精品久久久久久精品电影| 亚洲在线自拍视频| 国产精品一区二区三区四区久久| 国产精品野战在线观看| 香蕉丝袜av| 高潮久久久久久久久久久不卡| 免费电影在线观看免费观看| 成熟少妇高潮喷水视频| 美女大奶头视频| 亚洲人成网站在线播| 18禁黄网站禁片免费观看直播| 老司机福利观看| 色综合亚洲欧美另类图片| 99久久综合精品五月天人人| 神马国产精品三级电影在线观看| 最近在线观看免费完整版| 日韩精品青青久久久久久| 精品熟女少妇八av免费久了| 亚洲人成网站在线播放欧美日韩| 日韩欧美精品免费久久 | 少妇熟女aⅴ在线视频| 18禁在线播放成人免费| 伊人久久精品亚洲午夜| 亚洲七黄色美女视频| 757午夜福利合集在线观看| 桃红色精品国产亚洲av| 日韩欧美 国产精品| 午夜激情福利司机影院| 免费看a级黄色片| 成人亚洲精品av一区二区| 黄色女人牲交| 人妻丰满熟妇av一区二区三区| 丰满乱子伦码专区| 熟妇人妻久久中文字幕3abv| 国产精品久久视频播放| 麻豆久久精品国产亚洲av| 美女 人体艺术 gogo| 欧美在线黄色| 最近最新中文字幕大全电影3| 久久久久久国产a免费观看| 好看av亚洲va欧美ⅴa在| 久久亚洲真实| 国产免费一级a男人的天堂| 精品99又大又爽又粗少妇毛片 | 国产精品一区二区三区四区久久| 国产精品一区二区三区四区免费观看 | 日日摸夜夜添夜夜添小说| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| 欧美成狂野欧美在线观看| av中文乱码字幕在线| 亚洲国产中文字幕在线视频| 久久香蕉国产精品| 国产亚洲精品久久久久久毛片| 久久久国产成人精品二区| 久久香蕉国产精品| 亚洲成人免费电影在线观看| 欧美成人免费av一区二区三区| 淫秽高清视频在线观看| 欧美不卡视频在线免费观看| 俺也久久电影网| 人人妻人人看人人澡|