• <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
    荒銑加工軸向切深識別方法
    玻璃漿料鍵合中的孔洞抑制和微復合調控
    石英玻璃的熱輔助高效塑性域干磨削
    微小型薄底零件的軸向車銑實驗研究
    欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 黄频高清免费视频| 午夜老司机福利片| 日日夜夜操网爽| 国产精品1区2区在线观看.| 日本五十路高清| 黄色视频不卡| av福利片在线| 99国产精品99久久久久| 国产爱豆传媒在线观看 | a级毛片在线看网站| 在线av久久热| 女性被躁到高潮视频| 欧美日本亚洲视频在线播放| 香蕉丝袜av| 久久久精品国产亚洲av高清涩受| 男男h啪啪无遮挡| 一区二区三区高清视频在线| 日韩三级视频一区二区三区| 好男人在线观看高清免费视频 | 久久久久久大精品| 国产亚洲欧美98| 国产精品电影一区二区三区| 自线自在国产av| 变态另类丝袜制服| 99热6这里只有精品| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 亚洲人成77777在线视频| 亚洲人成电影免费在线| 免费一级毛片在线播放高清视频| 精品电影一区二区在线| 露出奶头的视频| 免费在线观看黄色视频的| 中文字幕久久专区| 波多野结衣高清无吗| 亚洲中文av在线| 国产欧美日韩一区二区精品| 国产1区2区3区精品| 亚洲一区二区三区色噜噜| 欧美黄色片欧美黄色片| 在线观看日韩欧美| 人人妻人人看人人澡| 国内少妇人妻偷人精品xxx网站 | 亚洲精品美女久久av网站| 制服诱惑二区| 免费在线观看日本一区| 神马国产精品三级电影在线观看 | 亚洲人成伊人成综合网2020| 老司机深夜福利视频在线观看| 国产野战对白在线观看| 满18在线观看网站| 99热这里只有精品一区 | 亚洲熟妇熟女久久| 国产极品粉嫩免费观看在线| 一级毛片高清免费大全| 国产野战对白在线观看| 俄罗斯特黄特色一大片| 欧美成人免费av一区二区三区| 757午夜福利合集在线观看| 精品国产一区二区三区四区第35| 国产精品久久久久久人妻精品电影| 99久久综合精品五月天人人| 一二三四在线观看免费中文在| 久久精品aⅴ一区二区三区四区| 夜夜爽天天搞| 久久热在线av| 欧美日韩精品网址| 日韩欧美一区二区三区在线观看| 十八禁人妻一区二区| 成人18禁高潮啪啪吃奶动态图| 又紧又爽又黄一区二区| 亚洲第一青青草原| 别揉我奶头~嗯~啊~动态视频| 脱女人内裤的视频| 国产三级在线视频| 老鸭窝网址在线观看| 国产亚洲欧美精品永久| 亚洲精品一卡2卡三卡4卡5卡| 男人舔女人下体高潮全视频| 91国产中文字幕| 国产aⅴ精品一区二区三区波| 精品不卡国产一区二区三区| 美女高潮喷水抽搐中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产人伦9x9x在线观看| 亚洲精品中文字幕在线视频| 国产精品九九99| 热99re8久久精品国产| 在线国产一区二区在线| 波多野结衣av一区二区av| 国产伦人伦偷精品视频| 亚洲欧洲精品一区二区精品久久久| 又紧又爽又黄一区二区| 长腿黑丝高跟| 欧美日韩福利视频一区二区| 老司机深夜福利视频在线观看| 中文字幕人成人乱码亚洲影| 成人手机av| 丁香六月欧美| 免费看a级黄色片| 18禁美女被吸乳视频| 美女扒开内裤让男人捅视频| 精品少妇一区二区三区视频日本电影| 国产三级黄色录像| 老司机福利观看| av免费在线观看网站| 91麻豆精品激情在线观看国产| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 美女免费视频网站| 桃红色精品国产亚洲av| 色综合婷婷激情| 亚洲国产精品sss在线观看| cao死你这个sao货| 国内精品久久久久精免费| 国内少妇人妻偷人精品xxx网站 | 国产一区二区三区在线臀色熟女| 成人亚洲精品av一区二区| 国产黄片美女视频| 亚洲精品美女久久av网站| 免费在线观看亚洲国产| 女同久久另类99精品国产91| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲欧美精品永久| 久久久久久人人人人人| 色播亚洲综合网| 久久午夜综合久久蜜桃| 大型黄色视频在线免费观看| 久久久久久久精品吃奶| 一级毛片女人18水好多| 久久国产精品男人的天堂亚洲| 99久久无色码亚洲精品果冻| 丁香六月欧美| 草草在线视频免费看| 久久九九热精品免费| 日韩大尺度精品在线看网址| 美女免费视频网站| 久久久国产成人免费| 真人做人爱边吃奶动态| 日韩精品免费视频一区二区三区| 最近在线观看免费完整版| 黄色毛片三级朝国网站| 美女大奶头视频| 99国产精品一区二区蜜桃av| 两个人看的免费小视频| 国产野战对白在线观看| 午夜精品在线福利| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人精品中文字幕电影| 大香蕉久久成人网| 搡老熟女国产l中国老女人| 国产主播在线观看一区二区| 非洲黑人性xxxx精品又粗又长| 99国产精品一区二区三区| 午夜a级毛片| 后天国语完整版免费观看| 久久久久久久久中文| 久久国产精品男人的天堂亚洲| 久久久久九九精品影院| 久久精品亚洲精品国产色婷小说| 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区| 午夜视频精品福利| 麻豆av在线久日| 久久青草综合色| 亚洲美女黄片视频| 91国产中文字幕| 亚洲av成人一区二区三| 午夜免费观看网址| 国产亚洲欧美98| 亚洲五月色婷婷综合| 视频区欧美日本亚洲| 精品欧美一区二区三区在线| 麻豆久久精品国产亚洲av| 久久久精品欧美日韩精品| 国产激情久久老熟女| 亚洲一区二区三区色噜噜| 国内毛片毛片毛片毛片毛片| 精品一区二区三区av网在线观看| 欧美在线黄色| x7x7x7水蜜桃| 亚洲电影在线观看av| 午夜福利免费观看在线| 久久热在线av| 美女午夜性视频免费| 制服诱惑二区| 黄色成人免费大全| 俄罗斯特黄特色一大片| xxx96com| 亚洲片人在线观看| 免费在线观看亚洲国产| 麻豆久久精品国产亚洲av| 精品免费久久久久久久清纯| 午夜免费鲁丝| 黄色女人牲交| 国产成人一区二区三区免费视频网站| 亚洲人成网站在线播放欧美日韩| 91成人精品电影| 久久亚洲真实| 国产亚洲精品一区二区www| 国产欧美日韩一区二区三| bbb黄色大片| 美女国产高潮福利片在线看| 成人亚洲精品一区在线观看| 国产在线精品亚洲第一网站| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品中文字幕看吧| 成人免费观看视频高清| 九色国产91popny在线| www.自偷自拍.com| 这个男人来自地球电影免费观看| 欧美国产精品va在线观看不卡| 十八禁网站免费在线| 中出人妻视频一区二区| 亚洲国产高清在线一区二区三 | 免费在线观看日本一区| 91麻豆精品激情在线观看国产| xxxwww97欧美| 亚洲熟女毛片儿| 99国产综合亚洲精品| 日本三级黄在线观看| 成人av一区二区三区在线看| а√天堂www在线а√下载| 国内精品久久久久精免费| 99re在线观看精品视频| 国产极品粉嫩免费观看在线| 不卡一级毛片| 久久精品成人免费网站| 一级毛片女人18水好多| 他把我摸到了高潮在线观看| 欧美一级毛片孕妇| 国产精品久久久久久亚洲av鲁大| 国产av又大| 成人国产综合亚洲| 国产精品精品国产色婷婷| www.精华液| 好男人电影高清在线观看| 国产精品 国内视频| 亚洲国产欧美网| 91麻豆精品激情在线观看国产| 在线av久久热| av天堂在线播放| 午夜福利视频1000在线观看| 免费无遮挡裸体视频| 久久久久久大精品| 人人澡人人妻人| 亚洲国产精品999在线| 一级作爱视频免费观看| 国产视频一区二区在线看| 18禁美女被吸乳视频| 日日摸夜夜添夜夜添小说| 久久精品aⅴ一区二区三区四区| 很黄的视频免费| 夜夜看夜夜爽夜夜摸| av有码第一页| 91大片在线观看| 丁香欧美五月| 中文字幕人妻熟女乱码| 妹子高潮喷水视频| 久久久久九九精品影院| 一级黄色大片毛片| 99国产精品一区二区蜜桃av| 搞女人的毛片| 久久中文字幕人妻熟女| 亚洲av成人一区二区三| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 精品久久久久久,| 麻豆成人午夜福利视频| 欧美三级亚洲精品| 国产精品久久视频播放| 成人国产综合亚洲| а√天堂www在线а√下载| 精品久久久久久久人妻蜜臀av| 少妇被粗大的猛进出69影院| 精品久久久久久成人av| 97人妻精品一区二区三区麻豆 | 亚洲精品在线美女| 色播亚洲综合网| 一本综合久久免费| 国产精品av久久久久免费| 免费在线观看日本一区| 免费高清视频大片| 我的亚洲天堂| 在线免费观看的www视频| 成人亚洲精品av一区二区| 精品人妻1区二区| 男人操女人黄网站| 身体一侧抽搐| 欧美成狂野欧美在线观看| 亚洲午夜精品一区,二区,三区| 视频区欧美日本亚洲| 夜夜躁狠狠躁天天躁| 亚洲精品久久成人aⅴ小说| 国产日本99.免费观看| 每晚都被弄得嗷嗷叫到高潮| 制服诱惑二区| 妹子高潮喷水视频| 成年人黄色毛片网站| 老司机在亚洲福利影院| 免费人成视频x8x8入口观看| 在线观看舔阴道视频| 在线免费观看的www视频| 日本三级黄在线观看| 18禁黄网站禁片免费观看直播| 非洲黑人性xxxx精品又粗又长| 丝袜在线中文字幕| 一本久久中文字幕| 久久国产精品影院| 国产精品久久久久久人妻精品电影| 亚洲国产日韩欧美精品在线观看 | av电影中文网址| 哪里可以看免费的av片| 夜夜躁狠狠躁天天躁| 久久天堂一区二区三区四区| www.熟女人妻精品国产| 国产蜜桃级精品一区二区三区| 久久久久国产一级毛片高清牌| 神马国产精品三级电影在线观看 | 国产成人精品久久二区二区免费| av有码第一页| 看黄色毛片网站| 久久久水蜜桃国产精品网| 精品第一国产精品| 亚洲精品色激情综合| 日韩欧美在线二视频| 亚洲国产看品久久| 99热这里只有精品一区 | 国产精品98久久久久久宅男小说| 国产精品 国内视频| 日韩精品中文字幕看吧| 夜夜爽天天搞| 麻豆成人av在线观看| cao死你这个sao货| 婷婷精品国产亚洲av| 日韩欧美免费精品| 一本大道久久a久久精品| 露出奶头的视频| 中文字幕高清在线视频| 日韩视频一区二区在线观看| netflix在线观看网站| 国产三级黄色录像| 色播在线永久视频| 欧美性猛交黑人性爽| 中文字幕av电影在线播放| 男女床上黄色一级片免费看| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产欧洲综合997久久, | 免费高清在线观看日韩| 99久久99久久久精品蜜桃| 91成年电影在线观看| 久久 成人 亚洲| 亚洲熟妇中文字幕五十中出| 99精品在免费线老司机午夜| 久久久久久久午夜电影| 欧美日本视频| 国产黄色小视频在线观看| 欧美乱码精品一区二区三区| 少妇裸体淫交视频免费看高清 | 一二三四在线观看免费中文在| 一级片免费观看大全| 老熟妇乱子伦视频在线观看| 欧美在线一区亚洲| 日本成人三级电影网站| 在线观看免费日韩欧美大片| 999久久久精品免费观看国产| 最好的美女福利视频网| 制服诱惑二区| 好男人电影高清在线观看| 午夜精品在线福利| 日韩免费av在线播放| 国产国语露脸激情在线看| 12—13女人毛片做爰片一| 1024手机看黄色片| 国产爱豆传媒在线观看 | 欧美日本视频| а√天堂www在线а√下载| 欧美 亚洲 国产 日韩一| 制服丝袜大香蕉在线| 国产伦一二天堂av在线观看| 亚洲av熟女| 久久久久久久久中文| or卡值多少钱| 中文资源天堂在线| 搡老妇女老女人老熟妇| 亚洲av五月六月丁香网| 俺也久久电影网| 国产午夜福利久久久久久| 男人操女人黄网站| 日韩欧美国产一区二区入口| 欧美丝袜亚洲另类 | 日韩免费av在线播放| 黄色片一级片一级黄色片| 好男人电影高清在线观看| 精品一区二区三区四区五区乱码| 国产三级在线视频| 熟女少妇亚洲综合色aaa.| 国产主播在线观看一区二区| xxx96com| 久久久久久久午夜电影| 国产精品亚洲美女久久久| 国产精品精品国产色婷婷| 国产免费男女视频| 亚洲av熟女| 国产黄片美女视频| 精品久久久久久久末码| 久久香蕉精品热| 桃红色精品国产亚洲av| 12—13女人毛片做爰片一| 人人妻人人澡人人看| 久久久久亚洲av毛片大全| 国产日本99.免费观看| 最好的美女福利视频网| 亚洲欧美精品综合久久99| 日韩 欧美 亚洲 中文字幕| 美女扒开内裤让男人捅视频| 91麻豆精品激情在线观看国产| 精品一区二区三区av网在线观看| 欧美性猛交黑人性爽| 国产成+人综合+亚洲专区| 亚洲一码二码三码区别大吗| 亚洲三区欧美一区| 18禁黄网站禁片午夜丰满| 人人妻人人澡欧美一区二区| 成年免费大片在线观看| 国内揄拍国产精品人妻在线 | 欧美又色又爽又黄视频| 免费在线观看成人毛片| 午夜福利欧美成人| 九色国产91popny在线| 变态另类丝袜制服| 一本综合久久免费| 最近最新免费中文字幕在线| 999久久久精品免费观看国产| 女警被强在线播放| 国产精品久久久久久亚洲av鲁大| 欧美成人免费av一区二区三区| 色在线成人网| 高清在线国产一区| 日韩三级视频一区二区三区| 欧美乱色亚洲激情| 一级片免费观看大全| 在线观看www视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻1区二区| 色av中文字幕| 亚洲,欧美精品.| 国产激情偷乱视频一区二区| 欧美日韩亚洲综合一区二区三区_| 免费在线观看影片大全网站| 国产精品久久视频播放| 一级a爱视频在线免费观看| 久久中文字幕一级| 亚洲 欧美一区二区三区| 欧美日韩精品网址| 精品国产美女av久久久久小说| 亚洲av电影不卡..在线观看| 久久人人精品亚洲av| 国产真实乱freesex| 麻豆成人av在线观看| 国产精品久久久久久人妻精品电影| 侵犯人妻中文字幕一二三四区| 91国产中文字幕| 日韩精品青青久久久久久| 久久久久久亚洲精品国产蜜桃av| 久久狼人影院| 国产成人av激情在线播放| 亚洲中文av在线| 搡老熟女国产l中国老女人| 精品久久久久久久毛片微露脸| 中文字幕av电影在线播放| 精品国产乱码久久久久久男人| 久热爱精品视频在线9| 亚洲成av人片免费观看| av天堂在线播放| 欧美中文综合在线视频| 亚洲熟妇中文字幕五十中出| 中文字幕精品免费在线观看视频| 国产成人影院久久av| 国产爱豆传媒在线观看 | 成人18禁在线播放| 99在线人妻在线中文字幕| 99久久99久久久精品蜜桃| 两个人视频免费观看高清| 久久国产精品影院| 人人妻人人澡人人看| 精华霜和精华液先用哪个| 一本精品99久久精品77| 久久久国产精品麻豆| 激情在线观看视频在线高清| 国产精华一区二区三区| 婷婷精品国产亚洲av在线| 亚洲男人的天堂狠狠| 久久精品影院6| 亚洲国产日韩欧美精品在线观看 | 欧美一区二区精品小视频在线| 国产单亲对白刺激| 88av欧美| 国产免费男女视频| 给我免费播放毛片高清在线观看| av天堂在线播放| 成人18禁在线播放| 淫秽高清视频在线观看| 婷婷精品国产亚洲av| 亚洲人成电影免费在线| 在线观看一区二区三区| 可以在线观看毛片的网站| 亚洲欧美一区二区三区黑人| 国产99白浆流出| 国产97色在线日韩免费| 精品国产美女av久久久久小说| 亚洲中文日韩欧美视频| 国产精品 国内视频| 亚洲全国av大片| 国产精品久久电影中文字幕| 哪里可以看免费的av片| 亚洲一区中文字幕在线| 欧美色视频一区免费| 久久国产乱子伦精品免费另类| 久久人妻av系列| 久久久精品欧美日韩精品| 少妇 在线观看| 很黄的视频免费| 精品国产美女av久久久久小说| 热re99久久国产66热| 老鸭窝网址在线观看| 大型黄色视频在线免费观看| svipshipincom国产片| videosex国产| 在线观看66精品国产| 很黄的视频免费| 国产一卡二卡三卡精品| 国产伦在线观看视频一区| 国产精华一区二区三区| 美女大奶头视频| 怎么达到女性高潮| cao死你这个sao货| 久久精品人妻少妇| 久久伊人香网站| 国产精品爽爽va在线观看网站 | 久久精品国产亚洲av高清一级| 老汉色∧v一级毛片| 亚洲欧美精品综合久久99| 亚洲熟妇熟女久久| 午夜激情av网站| 无人区码免费观看不卡| 亚洲激情在线av| 久久精品国产综合久久久| 免费看十八禁软件| 午夜福利18| 日本一本二区三区精品| 狠狠狠狠99中文字幕| 欧美一区二区精品小视频在线| 法律面前人人平等表现在哪些方面| 欧美激情极品国产一区二区三区| 欧美不卡视频在线免费观看 | 欧美日本视频| 久久九九热精品免费| 啦啦啦免费观看视频1| 脱女人内裤的视频| 国产精品 国内视频| 婷婷精品国产亚洲av在线| 久久国产精品人妻蜜桃| 制服丝袜大香蕉在线| 悠悠久久av| √禁漫天堂资源中文www| 视频区欧美日本亚洲| 一边摸一边做爽爽视频免费| 国产单亲对白刺激| 国产一区二区在线av高清观看| 国产精品1区2区在线观看.| 99久久精品国产亚洲精品| 亚洲精品色激情综合| www日本在线高清视频| 嫁个100分男人电影在线观看| 一二三四社区在线视频社区8| 欧美不卡视频在线免费观看 | 国产亚洲精品一区二区www| 久久午夜亚洲精品久久| av欧美777| 热re99久久国产66热| 亚洲五月婷婷丁香| 中文字幕人成人乱码亚洲影| 久久久久久久久中文| 国产精华一区二区三区| 国产不卡一卡二| 人人妻人人看人人澡| 久久九九热精品免费| av欧美777| 9191精品国产免费久久| 欧美成狂野欧美在线观看| 日韩精品免费视频一区二区三区| 9191精品国产免费久久| 亚洲色图 男人天堂 中文字幕| 国产高清激情床上av| av电影中文网址| 麻豆成人午夜福利视频| 在线永久观看黄色视频| 超碰成人久久| 国产成人一区二区三区免费视频网站| 欧美激情 高清一区二区三区| 又黄又爽又免费观看的视频| 免费在线观看亚洲国产| 亚洲精品美女久久久久99蜜臀| 曰老女人黄片| 国产亚洲欧美98| 欧美激情极品国产一区二区三区| 欧美成人午夜精品| 91麻豆av在线| or卡值多少钱| 国产一级毛片七仙女欲春2 | 日韩三级视频一区二区三区|