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

    高應變率下延性金屬中微孔洞貫通行為的數值分析

    2012-02-23 06:43:16王永剛徐東明
    兵工學報 2012年9期
    關鍵詞:孔洞塑性軸向

    王永剛,徐東明

    (寧波大學 力學與材料科學研究中心,浙江 寧波315211)

    0 引言

    延性金屬材料在高應變率加載下發(fā)生拉伸斷裂主要包括3 個物理過程:微孔洞的成核、長大和貫通[1]。迄今,人們對于微孔洞的成核和長大規(guī)律進行了大量的實驗和理論分析工作[2-4],但微孔洞相互貫通行為的研究工作相對很少。已有學者研究指出微孔洞的貫通在動態(tài)拉伸斷裂過程中起著很關鍵的作用[5],因此,研究微孔洞貫通機理對揭示材料動態(tài)拉伸斷裂的整個物理過程和本質有重要意義。

    Curran 等[6]基于實驗觀測結果認為有兩類微孔洞貫通機理,即微孔洞之間直接貫穿和微孔洞之間的變形局部化。事實上,微孔洞相互作用是通過它們周圍的塑性應變場來實現的。當兩個微孔洞之間距離達到一個臨界值時,它們周圍的塑性應變場發(fā)生交接,從而改變微孔洞的長大速率。一些學者曾對微孔洞之間開始貫通臨界條件進行了理論研究。Needleman 和Tvergarrd[7]認為材料內部微孔洞體積分數f 達到一個臨界值fc時微孔洞貫通開始。Brock和Sun[8]研究指出當材料中初始微孔洞體積分數較小時,臨界微孔洞體積分數fc對應力三軸度和加載應變率不敏感。而Pardoen 和Hutchinson[9]認為fc依賴于微孔洞初始體積分數、加載應力三軸度及微孔洞形狀。Brown 和Embury[10]提出了一個微孔洞開始貫通的臨界韌帶距離(ILDc)準則,指出當微孔洞間韌帶長度等于微孔洞的直徑時,微孔洞從獨立長大開始向貫通轉變。Horstemeyer 等[11]采用微力學有限元計算方法研究準靜態(tài)加載條件下304L不繡鋼和6061T6 鋁中兩個微孔洞的貫通行為,結果表明:ILDc的值依賴于邊界條件和材料的特性,約為2~8 倍微孔洞直徑。最近,Sepp?l? 等[12-13]采用分子動力學方法對金屬銅動態(tài)拉伸斷裂的微孔洞長大和聚集行為進行了數值模擬,計算結果顯示:當ILDc等于0.5 倍微孔洞直徑時,微孔洞周邊塑性形變區(qū)相遇,微孔洞長大的速率突然增加,最終導致了微孔洞之間發(fā)生了直接貫穿。

    無論是采用fc還是ILDc作為微孔洞貫通開始的準則,人們對微孔洞貫通行為的認識還是比較初步的,存在一些爭議。本文擬采用動力學有限元計算模型,對高應變率加載條件下純鋁材料中雙孔洞貫通行為進行數值分析,重點討論初始孔洞間韌帶距離和加載應變速率對雙孔洞貫通行為的影響規(guī)律,為合理建立孔洞貫通的理論模型奠定一定的基礎,從而進一步揭示延性金屬拉伸型斷裂的物理機制。

    1 計算模型建立

    高應變率拉伸載荷下延性金屬中出現的損傷主要是大量的隨機分布的球形微孔洞。因此,在一定假設條件下,可以把含損傷材料離散成大量含孔洞的圓柱體胞的集合體。這里,為了研究孔洞之間的貫通行為,在數值模擬中設定每一個圓柱體胞含有兩個大小相同的球形孔洞,二分之一模型見圖1(a)所示。在對稱載荷條件下,我們采用軸對稱的二維單元來進行計算分析,有限元模型見圖1(b)所示。微孔洞的直徑d0為50 μm,圓柱體胞的直徑D 和高度H 相等,均取為1 mm.圓柱體胞的直徑和高度是球形微孔洞直徑的20 倍。

    圖1 二分之一體胞模型及軸對稱有限元分析模型Fig.1 Half of the axisymmetric unit cell and FEM model

    計算中,邊界條件采用速度邊界,如圖1所示,vz、vr分別表示在軸向和徑向施加的速度函數。速度從0 開始經歷0.1 μs 線性增加到恒值,隨后保持不變。選取的模擬材料為純鋁,在高應變率加載下,采用考慮塑性硬化的流體彈塑性本構關系和Grüneisen 狀態(tài)方程進行數值分析。

    式中:σij為應力張量;p 為壓力;sij為偏應力張量。

    式中:σy為屈服應力考慮塑性硬化時,

    式中:σ0為初始屈服應力;為等效塑性應變;Eh塑性硬化模量。

    在材料處于壓縮狀態(tài)下,Grüneisen 狀態(tài)方程的壓力計算:

    當材料處于拉伸狀態(tài)下,

    計算參數見表1.

    表1 純鋁的本構和狀態(tài)方程參數Tab.1 Constitutive and equation of state parameters for pure aluminium

    2 計算結果討論

    在軸向速度場和徑向速度場相同條件下,圖2顯示軸對稱體胞模型邊界單元上宏觀軸向應力和徑向應力幾乎完全相等,應力開始時基本呈線性增加,隨后由于微孔洞長大的影響,出現非線性增長,很快達到極值點,隨后快速下降。由此可見,體胞模型中微孔洞承受的是一個準靜水拉伸載荷條件,這與一維應變層裂實驗中樣品中微孔洞所承受了應力狀態(tài)基本一致。因此,后續(xù)的有關微孔洞之間貫通行為的研究結果可以直接應用于層裂實驗的理論分析和數值模擬中。另外,在準靜水拉伸載荷下,圖3顯示不同時刻兩個微孔洞周圍的塑性應變場的演化云圖。圖3(a)顯示兩個微孔洞各自形成幾乎均勻的球對稱塑性應變場,微孔洞呈現相互獨立球形快速長大,而隨著微孔洞不斷長大,兩個微孔洞之間韌帶區(qū)域塑性應變場相互交疊,如圖3(b)所示,當此區(qū)域的應變場達到一個臨界值時,兩個微孔洞即發(fā)生貫通行為,直至直接貫穿,如圖3(c)和圖3(d)所示。

    圖3給出了微孔洞貫通行為的一些定性描述,為了定量討論微孔洞之間的貫通行為,這里采用模擬計算中實時跟蹤微孔洞長大過程中形狀改變來判斷,即通過比較圖1中標識的微孔洞上頂點A 和下頂點B 的軸向位移。如果微孔洞處于獨立長大階段,A 點和B 點軸向位移是一致的,如果A 和B 點軸向位移發(fā)生明顯的改變,表明微孔洞之間韌帶區(qū)域塑性應變場交疊改變了B 點的軸向位移增長速率,此時刻即可判斷微孔洞貫通開始,所對應的圖1中標識的B 點和C 點之間距離即為微孔洞開始貫通的臨界韌帶距離(ILDc).下面首先討論初始孔洞間韌帶距離ILD0對微孔洞貫通動力學行為的影響。

    圖2 宏觀軸向應力和徑向應力的比較Fig.2 Comparison of global axial stress and radial stress

    圖3 不同時刻微孔洞周圍等效塑性應變云圖Fig.3 Effective plastic strain contours around void

    2.1 ILD0的影響

    在加載應變率104s-1條件下,考察不同初始微孔洞間韌帶距離對孔洞貫通動力學行為的影響,選擇ILD0/d0為2、3、4 和5,即ILD0分別為100 μm、150 μm、200 μm、250 μm(上面圖2和圖3的數據對應的ILD0/d0=3).圖4(a)~圖4(d)左邊縱軸給出了不同ILD0/d0條件下實時的微孔洞間韌帶距離ILD 與實時微孔洞直徑d 之比隨時間演化過程,而其右邊縱軸同時給出了A 點和B 點的軸向位移時程。這里需要說明一下,為什么沒有直接來討論臨界ILD,而是討論它與實時微孔洞直徑d 之比,這主要是因為微孔洞貫通是通過微孔洞周圍的塑性應變場來實現的,而微孔洞的直徑直接影響著其周圍的塑性應變場。對比A 點和B 點的軸向位移時程,發(fā)現:微孔洞貫通開始前A 點和B 點軸向位移是相同的,兩個微孔洞之間沒有發(fā)生相互作用,微孔洞處于獨立長大階段,隨著微孔洞不斷長大,它們周圍的塑性應變場發(fā)生相互影響,導致隨后A 點和B 點軸向位移明顯發(fā)生偏離,利用此點恰好可以判斷微孔洞開始貫通時刻,因此,這里的計算結果也檢驗了上面提出的通過微孔洞形狀改變來判斷微孔洞貫通開始的有效性。而從ILD/d 的時程曲線上來看,在0.3 μs 之前,ILD/d 減小速率緩慢,此階段微孔洞基本還處于彈性變形階段,隨后進入微孔洞獨立快速長大階段,這時ILD/d 迅速減小,緊接著微孔洞間發(fā)生貫通行為,導致ILD/d 減小速率再次變緩,最后慢慢趨于0,即微孔洞發(fā)生了直接貫穿,如上面圖3(d)所示。對比圖4(a)~圖4(d),由A 點和B 點軸向位移明顯發(fā)生偏離點容易確定微孔洞間貫通開始的臨界韌帶距離,在圖中用黑點標識,結果顯示:微孔洞從獨立長大向微孔洞貫通轉變開始的臨界韌帶距離約為0.5 倍(ILDc=0.5 ±0.1d)的微孔洞直徑,并且它對初始的ILD0/d0不敏感。

    臨界韌帶距離作為微孔洞貫通開始的準則,比較直觀,但相對缺乏物理本質的認識和提煉。分子動力學研究成果[12-13]已經在微觀層次上揭示微孔洞間韌帶區(qū)域位錯密度迅速增長導致微孔洞發(fā)生貫通。位錯運動宏觀上對應于塑性應變,因此為了進一步理解微孔洞貫通動力學行為,圖5給出了不同ILD0/d0條件下微孔洞間韌帶距離中點處典型單元的等效塑性應變時程曲線。利用微孔洞開始貫通的時間對應條件,可以在圖5中找出微孔洞貫通時刻所對應的塑性應變,圖中用黑點標識,結果顯示:臨界塑性應變閾值約為0.5,隨著ILD0/d0增大,臨界塑性應變閾值緩慢減小。

    2.2 加載應變速率的影響

    通過改變邊界速度函數,可以調整加載應變速率。這里以ILD0/d0=3 模型為例,討論加載應變速率對微孔洞貫通行為的影響。圖6給出了不同加載應變速率下ILD/d 時程曲線和A 和B 點軸向位移時程曲線,結果顯示:隨著加載應變速率增大,微孔洞長大速率提高,微孔洞間貫通時間提前,但對微孔洞間開始貫通臨界韌帶距離沒有影響,還保持在0.5 倍微孔洞直徑。圖7給出了不同加載應變速率條件下微孔洞間韌帶距離中點處典型單元的等效塑性應變時程曲線,結果也同樣顯示微孔洞開始貫通的臨界塑性應變閾值為0.5,圖中黑點標識。圖8給出了不同應變速率條件下體胞模型邊界單元上的宏觀平均應力隨時間的變化曲線,從中可見:隨著加載應變速率增大,在微孔洞開始貫通時宏觀平均應力(黑點標識處)在增大,不對應于應力峰值,也就是說微孔洞間貫通開始后,宏觀平均應力場還在增大,材料仍然可以繼續(xù)承載。以前有很多學者為了回避對微孔洞貫通過程的理論描述,常常假定微孔洞貫通開始后材料迅速失去承載能力,這個觀點是值得商榷的。以后在建立合理的微孔洞貫通理論模型時必須考慮微孔洞貫通也是需要一個過程的,它不是瞬時事件。

    圖4 ILD/d 演化時程曲線及A 點和B 點軸向位移時程曲線Fig.4 Evolution of ILD/d profiles and A and B points’axial displacement profiles

    圖5 微孔洞間韌帶中心處典型單元的等效塑性應變時程Fig.5 Effective plastic strain profiles of representative element in the center of intervoid ligament distance

    圖6 ILD/d 時程曲線和A 和B 點軸向位移時程曲線Fig.6 Evolution of ILD/d profiles and A and B points’axial displacement profiles

    圖7 微孔洞間韌帶中心處典型單元的等效塑性應變時程Fig.7 Effective plastic strain profiles of representative element in the center of intervoid ligament distance

    圖8 微孔洞間韌帶中心處典型單元的等效塑性應變時程Fig.8 Effective plastic strain profiles of representative element in the center of intervoid ligament distance

    3 關于ILDc和fc的討論

    目前,關于微孔洞開始貫通的準則主要有兩種觀點,有些學者采用微孔洞間臨界韌帶距離ILDc,還有些學者常采用臨界微孔洞體積分數fc.如果基于微孔洞在材料中統(tǒng)計均勻分布的假定,設用邊長為L 的立方體來代表受損傷材料的一個特征體積,假設有8 個直徑為d 的微孔洞均勻分布于立方體的8 個角上,如圖9所示。每個微孔洞只有八分之一屬于此立方體,因此微孔洞開始貫通時的臨界微孔洞體積分數:

    由此可見,ILDc和fc是可以相互轉換的,在物理本質也是一致的。這里近似假定圖9中微孔洞兩兩之間的相互貫通行為是相互獨立的,當然這個假定具有一定的合理性,因為微孔洞之間相互貫通是通過塑性應變場交疊來實現的,圖9中其他6 個微孔洞由于距離相對比較遠,對兩個微孔洞之間的塑性應變場的影響相對較弱。因此,在微孔洞兩兩之間的相互貫通行為是相互獨立的假定下,可以通過含兩個微孔洞圓柱體胞組合排列形成空間的8 個微孔洞,同時體胞模型中微孔洞間臨界韌帶距離ILDc在物理本質上可以近似替代含8 個微孔洞的立方體損傷單元中ILDc.對于純鋁,上述數值模擬得到ILDc為0.5 倍微孔洞直徑,通過(7)式,計算得到的fc約為0.15.文獻[14]曾對純鋁開展過系列層裂實驗研究,基于對回收樣品中含有的微孔洞宏觀統(tǒng)計分析,獲得了微孔洞開始貫通的微孔洞體積分數fc約為0.11,與數值分析結果非常接近,這表明通過數值分析獲得微孔洞開始貫通臨界韌帶距離或臨界孔洞體積分數具有一定的合理性。

    圖9 臨界孔洞間韌帶距離和臨界孔洞體積分數的幾何說明Fig.9 Geometry relationship between ILDc and fc

    4 結論

    采用LS-DYNA 動力學有限元程序,對軸對稱體胞模型中雙孔洞的貫通行為進行了數值分析,討論了初始孔洞間韌帶距離和加載應變速率的影響,得到了如下幾點認識:

    1)提出了一個微孔洞開始貫通的判斷標準,即通過比較微孔洞上頂點和下頂點的軸向位移曲線的明顯偏離來判斷微孔洞開始貫通,數值結果顯示該方法具有一定的合理性和有效性。

    2)微孔洞從獨立長大向微孔洞貫通轉變開始的臨界韌帶距離約為0.5 倍微孔洞實時直徑,并且不依賴于初始微孔洞間韌帶距離和加載應變速率。

    3)微孔洞間貫通行為是通過塑性應變場相互作用實現的,數值計算結果顯示:在孔洞間韌帶區(qū)域的塑性應變場需要達到一個臨界閾值才能驅動微孔洞開始貫通,該臨界塑性應變閾值基本保持在0.5左右。

    4)微孔洞在材料中統(tǒng)計均勻分布和微孔洞之間的相互貫通行為是相互獨立的假定下,臨界孔洞間韌帶距離ILDc判據和臨界孔洞體積分數判據在物理本質也是一致的,可以相互轉換。

    References)

    [1] Antoun T,Seaman L,Curran D R,et al.Spall fracture[M].Berlin:Springer,2003.

    [2] Johnson J N.Dynamic fracture and spallation in ductile solids[J].Journal of Applied Physics,1981,52:2812 -2825.

    [3] Molinari A,Wright T W.A physical model for nucleation and early growth of voids in ductile materials under dynamic loading[J].Journal of the Mechanics and Physics of Solids,2005,53:1476-1504.

    [4] Chen D,Tan H,Yu Y,et al.A void coalescence-based spall model[J].International Journal of Impact Engineering,2006,32(11):1752 -1767.

    [5] Czarnota C,Mercier S,Molinari A.Modeling of nucleation and void growth in dynamic pressure loading,application to spall test in tantalum[J].International Journal of Fracture,2006,141(1 -2):177 -194.

    [6] Curran D R,Seaman L,Shockey D A.Dynamic failure of solids[J].Physics Reports,1987,147:253 -388.

    [7] Tvergarrd V,Needleman A.Analysis of the cup-cone fracture in a round tensile bar[J].Acta Metal.,1984,27:1 -15.

    [8] Brock W,Sun D Z,Honign A.Verification of the transferability of micromechanical parameters by cell model calculations with visco-plastic materials[J].International Journal Plasticity,1995,11(8):971 -989.

    [9] Pardoen T,Hutchinson J W.An extended model for void growth and coalescence[J].Journal of the Mechanics and Physics of Solids,2000,48:2467 -2512.

    [10] Brown L M,Embury J D.The initiation and growth of void at second phase particles[C]∥Proceedings 3rd International Conference on Strength of Metals and Alloys.London:Institute of Metals,1973.

    [11] Horstemeyer M F,Matalanis M M,Sieber A M,et al.Micromechanical finite element calculations of temperature and void configuration effects on void growth and coalescence[J].International Journal Plasticity,2000,16:979 -1015.

    [12] Sepp?l? E T,Belak J,Rudd R E.Onset of void coalescence during dynamic fracture of ductile metals[J].Physics Review Letters,2004,93(24):245503.

    [13] Sepp?l? E T,Belak J,Rudd R E.Three-dimensional molecular dynamic simulation of void coalescence during dynamic fracture of ductile metals[J].Physics Review B,2005,71:064112.

    [14] Qi M L,He H L.Simulation of critical behavior on damage evolution[J].Chinese Physics B,2010,19(3):036201-5.

    猜你喜歡
    孔洞塑性軸向
    基于應變梯度的微尺度金屬塑性行為研究
    大型立式單級引黃離心泵軸向力平衡的研究
    水泵技術(2021年5期)2021-12-31 05:26:40
    一種面向孔洞修復的三角網格復雜孔洞分割方法
    硬脆材料的塑性域加工
    鈹材料塑性域加工可行性研究
    孔洞加工工藝的概述及鑒定要點簡析
    收藏界(2019年3期)2019-10-10 03:16:22
    荒銑加工軸向切深識別方法
    玻璃漿料鍵合中的孔洞抑制和微復合調控
    石英玻璃的熱輔助高效塑性域干磨削
    微小型薄底零件的軸向車銑實驗研究
    日韩成人伦理影院| av免费观看日本| 最近最新中文字幕免费大全7| 亚洲,一卡二卡三卡| 国产成人aa在线观看| 国产在线一区二区三区精| 日日摸夜夜添夜夜添av毛片| 这个男人来自地球电影免费观看 | 亚洲久久久国产精品| 国产av码专区亚洲av| 亚洲美女黄色视频免费看| 成人国产av品久久久| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| a级毛色黄片| 精品一区二区三卡| 国产精品一区二区性色av| 99久久中文字幕三级久久日本| 高清不卡的av网站| 亚洲精品456在线播放app| 人妻一区二区av| 男女边摸边吃奶| 一级片'在线观看视频| 成人无遮挡网站| 中文字幕免费在线视频6| 日韩国内少妇激情av| 亚洲精品亚洲一区二区| 成人特级av手机在线观看| 女性被躁到高潮视频| 精品亚洲成国产av| 免费观看在线日韩| 高清av免费在线| 在现免费观看毛片| 国产成人a区在线观看| 人体艺术视频欧美日本| 久久久久久伊人网av| 国产成人a∨麻豆精品| 91精品国产九色| 久久av网站| 日日啪夜夜撸| 舔av片在线| 狂野欧美白嫩少妇大欣赏| 青春草视频在线免费观看| 日本免费在线观看一区| 久久人人爽人人爽人人片va| 国产精品久久久久久精品电影小说 | 天天躁夜夜躁狠狠久久av| 777米奇影视久久| 国产亚洲精品久久久com| 一个人免费看片子| 国产精品一区二区性色av| 国产男人的电影天堂91| 欧美xxxx性猛交bbbb| 久久婷婷青草| 久久99精品国语久久久| av国产精品久久久久影院| 嘟嘟电影网在线观看| 亚洲性久久影院| 国产精品偷伦视频观看了| 久久久久人妻精品一区果冻| 日本黄色片子视频| 日韩国内少妇激情av| 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 极品少妇高潮喷水抽搐| av福利片在线观看| 色吧在线观看| 欧美精品亚洲一区二区| 日韩在线高清观看一区二区三区| 在线观看免费日韩欧美大片 | 日韩av免费高清视频| 亚洲国产精品成人久久小说| 成年女人在线观看亚洲视频| av在线观看视频网站免费| 精品久久久久久电影网| 最后的刺客免费高清国语| 国产91av在线免费观看| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡动漫免费视频| 99久久精品国产国产毛片| 又黄又爽又刺激的免费视频.| 日本-黄色视频高清免费观看| 亚洲精品乱久久久久久| 99热全是精品| 亚洲av二区三区四区| 久久精品久久久久久噜噜老黄| 少妇猛男粗大的猛烈进出视频| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 搡女人真爽免费视频火全软件| 欧美xxxx黑人xx丫x性爽| 97在线人人人人妻| 国产免费福利视频在线观看| 搡老乐熟女国产| 久久久久性生活片| 亚洲最大成人中文| 国产精品秋霞免费鲁丝片| 97热精品久久久久久| 老师上课跳d突然被开到最大视频| 最近中文字幕2019免费版| 久久久久久伊人网av| 99久国产av精品国产电影| 男人添女人高潮全过程视频| 街头女战士在线观看网站| 色视频在线一区二区三区| 如何舔出高潮| 精品熟女少妇av免费看| 狂野欧美激情性xxxx在线观看| 免费在线观看成人毛片| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 国产精品伦人一区二区| 欧美成人精品欧美一级黄| 午夜精品国产一区二区电影| 国产免费又黄又爽又色| 国产成人午夜福利电影在线观看| 看十八女毛片水多多多| 两个人的视频大全免费| 最新中文字幕久久久久| 91久久精品国产一区二区三区| 最后的刺客免费高清国语| 干丝袜人妻中文字幕| 又黄又爽又刺激的免费视频.| 人妻制服诱惑在线中文字幕| 免费观看av网站的网址| 五月伊人婷婷丁香| 亚洲av电影在线观看一区二区三区| 久久久久网色| 在线观看人妻少妇| 色综合色国产| 国产午夜精品久久久久久一区二区三区| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 国产永久视频网站| 国产黄频视频在线观看| www.色视频.com| 一级二级三级毛片免费看| 亚洲四区av| 欧美 日韩 精品 国产| 久久这里有精品视频免费| 大话2 男鬼变身卡| 亚洲国产av新网站| 色视频在线一区二区三区| 久久久久国产精品人妻一区二区| 国产大屁股一区二区在线视频| 日本黄色日本黄色录像| 日韩人妻高清精品专区| 18禁裸乳无遮挡动漫免费视频| 大片免费播放器 马上看| 亚洲精品456在线播放app| 国产精品99久久久久久久久| 男女边摸边吃奶| 日本av手机在线免费观看| 又爽又黄a免费视频| av在线蜜桃| 免费久久久久久久精品成人欧美视频 | 免费av中文字幕在线| 91在线精品国自产拍蜜月| 能在线免费看毛片的网站| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩东京热| 亚洲精品久久久久久婷婷小说| 美女国产视频在线观看| 国产精品三级大全| 插阴视频在线观看视频| 在线亚洲精品国产二区图片欧美 | 中文字幕精品免费在线观看视频 | 欧美激情极品国产一区二区三区 | 超碰av人人做人人爽久久| 亚洲国产精品成人久久小说| 新久久久久国产一级毛片| 亚洲电影在线观看av| 免费av中文字幕在线| 欧美成人精品欧美一级黄| 久久人妻熟女aⅴ| 亚洲精品国产成人久久av| 晚上一个人看的免费电影| 人体艺术视频欧美日本| 亚洲精品国产色婷婷电影| 熟妇人妻不卡中文字幕| 亚洲av欧美aⅴ国产| 少妇裸体淫交视频免费看高清| 99久久人妻综合| 1000部很黄的大片| 少妇猛男粗大的猛烈进出视频| 美女中出高潮动态图| 国产精品福利在线免费观看| 一本久久精品| 天堂中文最新版在线下载| 男女国产视频网站| 涩涩av久久男人的天堂| 少妇人妻一区二区三区视频| 中文乱码字字幕精品一区二区三区| 欧美+日韩+精品| 亚洲av男天堂| 如何舔出高潮| freevideosex欧美| 狂野欧美激情性xxxx在线观看| 久久久久久久国产电影| 2022亚洲国产成人精品| 久久午夜福利片| 国产欧美日韩精品一区二区| www.色视频.com| 亚洲精品成人av观看孕妇| 免费av不卡在线播放| 成人特级av手机在线观看| 熟妇人妻不卡中文字幕| 街头女战士在线观看网站| 中文字幕精品免费在线观看视频 | 王馨瑶露胸无遮挡在线观看| 亚洲国产精品999| 美女国产视频在线观看| 秋霞在线观看毛片| 亚洲av二区三区四区| 欧美变态另类bdsm刘玥| 久久精品久久久久久久性| 国产av码专区亚洲av| 97在线视频观看| 国产爽快片一区二区三区| 日韩欧美 国产精品| 日本爱情动作片www.在线观看| 丝袜喷水一区| 少妇被粗大猛烈的视频| 一级黄片播放器| 老司机影院毛片| 在线观看国产h片| 涩涩av久久男人的天堂| 最新中文字幕久久久久| 干丝袜人妻中文字幕| 午夜福利在线观看免费完整高清在| 国产视频首页在线观看| 色视频www国产| 99国产精品免费福利视频| 亚洲国产色片| 狠狠精品人妻久久久久久综合| h视频一区二区三区| 三级国产精品片| 日韩电影二区| 插逼视频在线观看| 成人国产麻豆网| 欧美成人a在线观看| 亚洲精品乱码久久久v下载方式| .国产精品久久| 国产成人免费观看mmmm| 高清欧美精品videossex| 国产成人精品一,二区| 边亲边吃奶的免费视频| 伦理电影免费视频| 免费在线观看成人毛片| 人妻一区二区av| 一区在线观看完整版| 久久国产亚洲av麻豆专区| 高清毛片免费看| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 秋霞在线观看毛片| 伦精品一区二区三区| 亚洲自偷自拍三级| 久久精品国产鲁丝片午夜精品| 亚洲精品,欧美精品| 中文字幕精品免费在线观看视频 | 国产亚洲午夜精品一区二区久久| 免费少妇av软件| 欧美bdsm另类| 国产无遮挡羞羞视频在线观看| 日本黄大片高清| 久久国产亚洲av麻豆专区| 欧美老熟妇乱子伦牲交| 男女国产视频网站| 99久久精品一区二区三区| 久久久久精品久久久久真实原创| 青春草国产在线视频| 制服丝袜香蕉在线| 美女cb高潮喷水在线观看| 亚洲国产日韩一区二区| 综合色丁香网| 在线播放无遮挡| 国产成人精品婷婷| 尤物成人国产欧美一区二区三区| 日本av免费视频播放| 黄色视频在线播放观看不卡| 插逼视频在线观看| 久久 成人 亚洲| av视频免费观看在线观看| 成人国产av品久久久| 国产精品熟女久久久久浪| 免费大片黄手机在线观看| 亚洲精品国产色婷婷电影| 国产成人精品一,二区| 国产伦精品一区二区三区视频9| 成年免费大片在线观看| 午夜日本视频在线| 97精品久久久久久久久久精品| 在线天堂最新版资源| 午夜福利视频精品| 人妻制服诱惑在线中文字幕| 嫩草影院新地址| 亚洲人成网站在线观看播放| 亚洲国产欧美人成| 99久久精品一区二区三区| 亚洲精品色激情综合| 亚洲精品一区蜜桃| 精品一区二区免费观看| av线在线观看网站| 亚洲欧洲国产日韩| 国产精品久久久久久精品古装| 日韩一本色道免费dvd| 热99国产精品久久久久久7| 久久久久久久久久久免费av| 日日摸夜夜添夜夜爱| 久久久久久久久久成人| 欧美性感艳星| 欧美人与善性xxx| 黑人猛操日本美女一级片| 夜夜爽夜夜爽视频| 在线观看三级黄色| 99热这里只有是精品50| 十分钟在线观看高清视频www | 亚洲三级黄色毛片| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 国产一区二区三区综合在线观看 | 久久av网站| 色哟哟·www| av在线老鸭窝| 亚洲精品色激情综合| 尤物成人国产欧美一区二区三区| 亚洲av男天堂| 亚洲国产精品国产精品| 国产伦精品一区二区三区视频9| 久久久久久久大尺度免费视频| 自拍偷自拍亚洲精品老妇| 久久久久精品性色| 纵有疾风起免费观看全集完整版| 国产人妻一区二区三区在| 亚洲第一区二区三区不卡| 国产免费视频播放在线视频| 亚洲色图av天堂| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频 | 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 日本av手机在线免费观看| 亚洲,欧美,日韩| 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 国产精品爽爽va在线观看网站| 久久这里有精品视频免费| 九色成人免费人妻av| 国产日韩欧美在线精品| xxx大片免费视频| 亚洲人成网站高清观看| 久久久国产一区二区| 久久国产精品男人的天堂亚洲 | 色5月婷婷丁香| 激情五月婷婷亚洲| 欧美高清成人免费视频www| 亚洲欧洲日产国产| 国产在线男女| 久久ye,这里只有精品| 国产又色又爽无遮挡免| 久久人人爽av亚洲精品天堂 | 午夜老司机福利剧场| 伊人久久精品亚洲午夜| av在线播放精品| 亚洲国产日韩一区二区| 欧美成人a在线观看| 女性生殖器流出的白浆| 日韩欧美精品免费久久| 国产免费福利视频在线观看| 乱系列少妇在线播放| 久久热精品热| 99热这里只有是精品在线观看| 男女边摸边吃奶| 少妇人妻精品综合一区二区| 亚洲一区二区三区欧美精品| 国产视频内射| 美女主播在线视频| 久久久久久久国产电影| 亚洲综合色惰| 欧美97在线视频| 亚洲内射少妇av| 国产伦精品一区二区三区四那| 欧美另类一区| 欧美老熟妇乱子伦牲交| 在线观看国产h片| 成人综合一区亚洲| 国产精品无大码| 在线观看免费高清a一片| 最近中文字幕高清免费大全6| 国产在线一区二区三区精| 亚洲精品日韩在线中文字幕| 亚洲国产最新在线播放| 久久久精品免费免费高清| 夫妻午夜视频| 在线播放无遮挡| 天堂8中文在线网| 欧美+日韩+精品| av国产精品久久久久影院| 国产成人精品福利久久| 性色avwww在线观看| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 欧美精品人与动牲交sv欧美| 亚洲精品成人av观看孕妇| 亚洲欧美精品专区久久| 熟妇人妻不卡中文字幕| 18禁裸乳无遮挡动漫免费视频| 夜夜爽夜夜爽视频| 成人影院久久| 国产探花极品一区二区| 六月丁香七月| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| 亚洲真实伦在线观看| 一级av片app| 久久精品国产自在天天线| 久久久久久人妻| 日日撸夜夜添| 久久这里有精品视频免费| 久久久久人妻精品一区果冻| 黄色一级大片看看| 午夜老司机福利剧场| 国产免费福利视频在线观看| 秋霞在线观看毛片| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久| 久久人人爽人人爽人人片va| 中文资源天堂在线| 亚洲欧美日韩卡通动漫| 国产精品.久久久| 好男人视频免费观看在线| 99视频精品全部免费 在线| 亚洲av免费高清在线观看| 日韩人妻高清精品专区| 蜜桃亚洲精品一区二区三区| 九九久久精品国产亚洲av麻豆| 又黄又爽又刺激的免费视频.| 日本猛色少妇xxxxx猛交久久| 在线播放无遮挡| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 成年人午夜在线观看视频| 男女啪啪激烈高潮av片| 极品教师在线视频| 国产精品国产三级国产av玫瑰| 国产乱人视频| 日韩av不卡免费在线播放| 女人十人毛片免费观看3o分钟| 国产伦精品一区二区三区四那| 欧美精品国产亚洲| 国产成人精品一,二区| 草草在线视频免费看| 嘟嘟电影网在线观看| 只有这里有精品99| 免费在线观看成人毛片| 欧美精品一区二区大全| 久久韩国三级中文字幕| 久久这里有精品视频免费| 18禁裸乳无遮挡动漫免费视频| 国产乱来视频区| 久久热精品热| 精品一区二区三区视频在线| 大陆偷拍与自拍| 高清不卡的av网站| 不卡视频在线观看欧美| a级毛片免费高清观看在线播放| 青春草视频在线免费观看| 日韩大片免费观看网站| 久久久久久人妻| 少妇丰满av| 欧美极品一区二区三区四区| 亚洲精品日韩av片在线观看| 国产乱来视频区| 这个男人来自地球电影免费观看 | 欧美激情极品国产一区二区三区 | 国产精品麻豆人妻色哟哟久久| 黑人猛操日本美女一级片| 久久精品国产亚洲av涩爱| 午夜老司机福利剧场| 精品人妻偷拍中文字幕| 亚洲精品国产av蜜桃| 黑人高潮一二区| 国产成人午夜福利电影在线观看| 男男h啪啪无遮挡| 身体一侧抽搐| 美女高潮的动态| av天堂中文字幕网| 亚洲四区av| 一区二区三区精品91| 精品99又大又爽又粗少妇毛片| av又黄又爽大尺度在线免费看| 成人综合一区亚洲| 在线观看美女被高潮喷水网站| 色网站视频免费| 乱系列少妇在线播放| 观看av在线不卡| 免费黄频网站在线观看国产| 亚洲国产精品一区三区| 欧美区成人在线视频| 水蜜桃什么品种好| 美女主播在线视频| 中文字幕久久专区| 夫妻午夜视频| 国产又色又爽无遮挡免| av免费在线看不卡| 多毛熟女@视频| 日韩视频在线欧美| 亚洲怡红院男人天堂| av线在线观看网站| 日韩大片免费观看网站| 日韩视频在线欧美| 国产亚洲午夜精品一区二区久久| 精品一品国产午夜福利视频| 国产国拍精品亚洲av在线观看| 联通29元200g的流量卡| 在线观看三级黄色| 成人一区二区视频在线观看| 国产男女内射视频| 秋霞伦理黄片| 黑人猛操日本美女一级片| 国产精品人妻久久久久久| 亚洲av电影在线观看一区二区三区| 色5月婷婷丁香| 久久青草综合色| 九草在线视频观看| 美女内射精品一级片tv| 九九在线视频观看精品| 婷婷色综合www| av在线蜜桃| 最近2019中文字幕mv第一页| 2022亚洲国产成人精品| 超碰97精品在线观看| 26uuu在线亚洲综合色| 久久久久久久久久久丰满| 美女国产视频在线观看| 国产成人精品一,二区| 一级a做视频免费观看| 成人午夜精彩视频在线观看| 噜噜噜噜噜久久久久久91| 中文资源天堂在线| 亚洲精品久久午夜乱码| 伦理电影大哥的女人| 午夜老司机福利剧场| 1000部很黄的大片| 在线免费观看不下载黄p国产| 国产精品福利在线免费观看| 国产一级毛片在线| 亚洲美女视频黄频| 久久国产精品男人的天堂亚洲 | 日本av手机在线免费观看| 色综合色国产| 好男人视频免费观看在线| 一级毛片黄色毛片免费观看视频| 在现免费观看毛片| 久久久亚洲精品成人影院| 国产中年淑女户外野战色| 91aial.com中文字幕在线观看| 成年免费大片在线观看| 少妇被粗大猛烈的视频| 精品人妻一区二区三区麻豆| 国产成人精品婷婷| 18+在线观看网站| 色视频在线一区二区三区| av天堂中文字幕网| 国产老妇伦熟女老妇高清| 婷婷色麻豆天堂久久| 国产v大片淫在线免费观看| 亚洲国产av新网站| 亚洲第一av免费看| 大片电影免费在线观看免费| 亚洲国产精品999| 午夜福利视频精品| av.在线天堂| 国产成人aa在线观看| 久久国产乱子免费精品| 熟妇人妻不卡中文字幕| 免费不卡的大黄色大毛片视频在线观看| 日韩电影二区| 美女脱内裤让男人舔精品视频| 久久青草综合色| 免费av不卡在线播放| 简卡轻食公司| av播播在线观看一区| 女的被弄到高潮叫床怎么办| 免费播放大片免费观看视频在线观看| 中文字幕免费在线视频6| 一级毛片 在线播放| 99久久精品一区二区三区| 青青草视频在线视频观看| 亚洲中文av在线| 2021少妇久久久久久久久久久| 在线观看av片永久免费下载| 国产无遮挡羞羞视频在线观看| 草草在线视频免费看| 国产淫语在线视频| 在线观看美女被高潮喷水网站| 亚洲成色77777| 国产片特级美女逼逼视频| 亚洲精品,欧美精品| 国产乱人偷精品视频| 免费av中文字幕在线| 高清视频免费观看一区二区| 亚洲欧美清纯卡通| 在线观看免费日韩欧美大片 | 色哟哟·www| 少妇高潮的动态图| 下体分泌物呈黄色| 少妇被粗大猛烈的视频| 美女中出高潮动态图| 高清不卡的av网站| 大话2 男鬼变身卡| 国产精品久久久久久久久免| 中文精品一卡2卡3卡4更新| 99国产精品免费福利视频|