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

    相場法研究初始微結(jié)構(gòu)對晶粒長大的影響*

    2017-01-03 02:45:22羅志榮盧成健高英俊
    廣西科學(xué) 2016年5期
    關(guān)鍵詞:相場柱狀微結(jié)構(gòu)

    羅志榮,盧成健,高英俊

    (1.玉林師范學(xué)院物理科學(xué)與工程技術(shù)學(xué)院,廣西玉林 537000;2.廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院,廣西南寧 530004)

    ?

    相場法研究初始微結(jié)構(gòu)對晶粒長大的影響*

    羅志榮1,盧成健1,高英俊2**

    (1.玉林師范學(xué)院物理科學(xué)與工程技術(shù)學(xué)院,廣西玉林537000;2.廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院,廣西南寧530004)

    (1.College of Physics Science and Engineering,Yulin Normal University,Yulin,Guangxi,537000,China;2.College of Physical Science and Technology,Guangxi University,Nanning,Guangxi,530004,China)

    摘要:【目的】研究不同初始微結(jié)構(gòu)對晶粒長大過程及生長動力學(xué)的影響?!痉椒ā坎捎孟鄨龇?Phase Field)模擬二維多晶材料中正常晶粒長大及初態(tài)分別為柱狀和梯度微結(jié)構(gòu)的晶粒長大過程?!窘Y(jié)果】正常晶粒長大的相對晶粒尺寸分布具有時間不變性的特點;柱狀微結(jié)構(gòu)的長寬比和梯度微結(jié)構(gòu)的梯度指數(shù)均直接影響晶粒長大動力學(xué),這兩種初始微結(jié)構(gòu)在演化過程中均向均勻等軸微結(jié)構(gòu)轉(zhuǎn)變?!窘Y(jié)論】晶粒長大是大晶粒不斷吞噬相鄰小晶粒的過程;晶界曲率對晶粒長大過程有顯著影響,曲率越大,晶粒長大越快。

    關(guān)鍵詞:微結(jié)構(gòu)晶粒長大曲率晶粒尺寸相場法

    0 引言

    【研究意義】多晶材料的晶粒尺寸直接影響材料的強度、塑性、韌性以及耐腐蝕性等性能,因此研究多晶材料晶粒尺寸的演化規(guī)律具有重要意義[1-2]?!厩叭搜芯窟M(jìn)展】近年來,隨著材料科學(xué)與計算機(jī)技術(shù)的快速發(fā)展,利用計算機(jī)模擬研究材料微結(jié)構(gòu)演化過程已成為大趨勢。當(dāng)前,利用計算機(jī)模擬晶粒長大的主要模擬方法有蒙特卡洛法(Monte Carlo)[3-5],元胞自動機(jī)法(Cellular Automaton)[6-7]以及相場法(Phase Field)[8-12]。在這些模擬方法中,相場法因具有深刻的物理思想,不需要追蹤復(fù)雜界面的運動等優(yōu)勢而成為強有力的模擬方法。1994年,Chen和Yang[8]構(gòu)建了單相晶粒長大的相場模型,研究二維晶粒長大過程及晶粒尺寸分布特征。1996年,Chen和Fan[9]構(gòu)建了兩相晶粒長大的相場模型,研究Al2O3-ZrO2兩相晶粒長大過程及微結(jié)構(gòu)特征。2001年,Kazaryan等[10]利用相場法研究晶界遷移率各向異性對單相晶粒長大的形貌、動力學(xué)和拓?fù)浣Y(jié)構(gòu)的影響。2002年,Krill和Chen[11]利用相場法研究三維晶粒長大過程及其動力學(xué)和拓?fù)涮匦浴?011年,魏承煬和李賽毅[12]采用相場法研究單相多晶材料退火過程中,由溫度梯度引起的非均勻二維晶粒長大和定向晶粒長大行為?!颈狙芯壳腥朦c】目前,利用相場法研究初始微結(jié)構(gòu)對晶粒長大過程的影響還鮮有報到?!緮M解決的關(guān)鍵問題】采用相場法模擬單相多晶材料中正常晶粒及初態(tài)分別為柱狀和梯度微結(jié)構(gòu)的晶粒長大過程,從晶粒形貌和晶粒尺寸等方面進(jìn)行分析討論,揭示晶粒長大的規(guī)律。

    1 模型與方法

    1.1相場模型

    在相場模型中,多晶材料的微結(jié)構(gòu)可用一系列非保守的取向場變量ηi(r,t),(i=1,2,…,p)來表征,其中p為晶粒取向個數(shù)。這些場變量是空間r和時間t的連續(xù)函數(shù),其值為0~1或-1~0。例如,在第q個取向的晶粒中,ηi=q=1(或-1),ηi≠q=0;在晶界處,ηi在0和1(或0和-1)之間連續(xù)變化。材料微觀結(jié)構(gòu)演化的驅(qū)動力來自系統(tǒng)總自由能的降低。系統(tǒng)總自由能F可表示為[8]

    (1)

    式中,ki為能量梯度系數(shù);f0為自由能密度函數(shù),其表達(dá)式為[8]

    (2)

    式中,α、β和γ為正實數(shù),要求γ>β。

    由于取向場變量是非保守場變量,它們隨時間的演化由Ginzburg-Landau方程組描述[8]

    (3)

    式中,Li為晶界遷移率。

    為求解動力學(xué)方程(3),需要在時間和空間上對其進(jìn)行離散化處理。在時間上,采用顯式Euler迭代公式[8]

    (4)

    式中,Δt為時間步長。

    在空間上,對式(3)中的Laplace算子采用九點差分格式離散[8]

    (5)

    式中,Δx為空間步長,j和n分別代表i的最近鄰與次近鄰格點。

    1.2參數(shù)設(shè)置

    本研究不針對具體材料,所采用的參數(shù)均已無量綱化。計算時,將連續(xù)空間離散為四方格子,計算區(qū)域為512×512 gp(gp表示格子點數(shù))。對于正常晶粒長大,其取向場變量ηi的初始值設(shè)置為在[-0.001,0.001]之間的隨機(jī)數(shù)。其他參數(shù)設(shè)置為p=36,α=β=1.0,γ=2.0,Li=1.0,ki=2.0,Δx=2.0,Δt=0.25。對所有格點遍歷處理一次后計為一個模擬時間步(ts),并最終模擬出整個晶粒的長大過程。

    2 結(jié)果與分析

    2.1正常晶粒長大

    圖1為晶粒長大過程中不同時刻的晶粒組織,其中不同顏色區(qū)域表征不同取向晶粒。可見,隨著演化的進(jìn)行,大晶粒逐漸長大,小晶粒逐漸縮小直至消失,使得晶粒數(shù)量減小,平均晶粒尺寸增大。由圖1還可見,在演化初期,晶界較彎曲,曲率較大,邊數(shù)大于6的晶粒晶界呈現(xiàn)內(nèi)凹形狀,而邊數(shù)少于6的晶粒具有凸出的特征。隨著演化的進(jìn)行,晶界在向其曲率中心遷移的過程中逐漸趨于平直,三叉晶界夾角趨于120°,與文獻(xiàn)[13]的實驗觀察結(jié)果及晶粒長大理論[14]吻合。

    圖2a和b分別為不同演化時間的晶粒尺寸分布和相對晶粒尺寸分布曲線??梢姡S著演化時間的增加,晶粒尺寸分布沿大尺寸方向轉(zhuǎn)移(圖2a),而相對晶粒尺寸分布基本保持不變,具有自相似性(圖2b),與文獻(xiàn)[15]的實驗結(jié)果一致。

    (a)1×103ts;(b)6×103ts;(c)2×104ts圖1 正常晶粒長大過程中的晶粒形貌Fig.1 Grain morphology during normal grain growth

    圖2晶粒長大過程中晶粒尺寸分布(a)和相對晶粒尺寸分布(b)

    Fig.2Grain size distribution (a) and relative grain size distribution (b) during grain growth

    為研究多晶長大過程中單個晶粒的尺寸隨時間的變化情況,在晶粒演化到500 ts時隨機(jī)選取9個大小各異的晶粒,并監(jiān)測其面積(晶粒的格點數(shù)量)。圖3為這9個晶粒的面積隨時間的變化曲線??梢姡S著演化的進(jìn)行,在演化初期尺寸較小的晶粒,其面積一直減小,并在較短時間內(nèi)消失;而在演化初期尺寸較大的晶粒,一部分持續(xù)長大,另一部分在長大到某一時刻面積開始減小,直至完全消失。這說明小晶粒在演化初期就被周圍的大晶粒吞噬,而初始晶粒中尺寸較大的晶粒在長大過程中其相鄰晶粒不斷發(fā)生變化,當(dāng)其被更大尺寸的晶粒包圍時,其面積就會減小,直至晶粒消失。只有在長大過程中一直處于優(yōu)勢地位的大晶粒,才能持續(xù)長大。

    單相多晶系統(tǒng)晶粒演化的驅(qū)動力來源于系統(tǒng)界面能的減少,而界面能與晶界長度密切相關(guān)[6]。晶粒長大過程中,晶界總長度隨演化時間的變化曲線如圖4所示??梢姡谘莼^程中晶界長度不斷減小,減小幅度由快變慢。說明在晶粒長大過程中系統(tǒng)界面能不斷降低,遵循能量降低原則。

    圖3隨機(jī)選取的9個不同尺寸晶粒的面積與演化時間的關(guān)系

    Fig.3Temporal evolution of the areas of 9 randomly selected grains with different sizes

    圖4晶粒長大過程中晶界總長度與演化時間的關(guān)系

    Fig.4Total length of grain boundary with time during grain growth

    2.2柱狀微結(jié)構(gòu)的晶粒長大

    為研究初始狀態(tài)為柱狀晶粒的生長動力學(xué),設(shè)置了長寬比不同且平均晶粒尺寸相同(初始平均晶粒尺寸R0為13 gp)的3種初始柱狀微結(jié)構(gòu),如圖5所示。然后對系統(tǒng)進(jìn)行均勻退火,即進(jìn)行正常晶粒演化,所得結(jié)果見圖6。

    圖6顯示了晶粒長大過程中柱狀微結(jié)構(gòu)向等軸微結(jié)構(gòu)轉(zhuǎn)變過程。為定量證實該轉(zhuǎn)變現(xiàn)象,統(tǒng)計了3種不同初始長寬比的柱狀微結(jié)構(gòu)的平均晶粒尺寸和長寬比與演化時間的關(guān)系,如圖7所示。由圖7a可見,晶粒尺寸隨初始長寬比的增大而增大。這是因為圖5中晶粒為扁平狀或竹節(jié)狀,其晶界上下方向較平直,而左右方向具有較大曲率,且晶粒長寬比越大,晶粒左右方向的晶界曲率也越大,在演化過程中晶界曲率驅(qū)使晶界更快地遷移,晶粒生長就越快。從圖7b可以看出,在晶粒演化過程中,晶粒長寬比不斷減小,且減小幅度隨初始長寬比的增加而增大。該結(jié)果與冷軋鋼在退火過程中晶粒長寬比隨退火時間變化的規(guī)律[16]一致。

    (a)Aspect ratio=2.2;(b)Aspect ratio=4.4;(c)Aspect ratio=7.6

    圖53種不同長寬比的初始柱狀微結(jié)構(gòu)

    Fig.5Initial columnar microstructures with three different grain aspect ratios

    (a1)~(c1)1×103ts;(a2)~(c2)1×104ts;(a1),(a2) Initial aspect ratio=2.2;(b1),(b2) Initial aspect ratio=4.4;(c1),(c2) Initial aspect ratio=7.6

    圖63種不同初始長寬比的柱狀微結(jié)構(gòu)的演化過程

    Fig.6Evolution process of columnar microstructures with three different initial aspect ratios

    2.3梯度微結(jié)構(gòu)的晶粒長大

    利用局部退火條件下線性溫度梯度方法[4],獲得梯度指數(shù)不同且平均晶粒尺寸相同(平均晶粒半徑R0=13 gp)的3種梯度微結(jié)構(gòu),如圖8所示。其中,梯度指數(shù)是系統(tǒng)底部與頂部區(qū)域平均晶粒尺寸的比值。由圖8可見,系統(tǒng)頂部的晶粒細(xì)小均勻,而系統(tǒng)底部的晶粒較粗大,且梯度指數(shù)越大,系統(tǒng)底部晶粒的平直晶界就越多,夾角趨于穩(wěn)定值120°的三叉晶界也越多。

    為研究梯度微結(jié)構(gòu)的晶粒長大行為,對具有不同梯度指數(shù)的微結(jié)構(gòu)進(jìn)行均勻退火處理,即進(jìn)行正常晶粒演化,所得結(jié)果如圖9所示??梢?,晶粒發(fā)生了不同程度的長大,其中越靠近系統(tǒng)頂部的晶粒,長大現(xiàn)象就越明顯。這是因為系統(tǒng)頂部的晶粒細(xì)小,晶界曲率大,晶界多,使得晶粒長大的驅(qū)動力大,晶粒長大速度快;而系統(tǒng)底部的晶粒粗大,大部分晶界已趨于直線,晶粒長大的驅(qū)動力較小,因此晶粒長大速度緩慢。

    圖7柱狀微結(jié)構(gòu)的平均晶粒尺寸(a)和長寬比(b)與演化時間的關(guān)系

    Fig.7Changes of average grain size (a) and aspect ratio (b) of columnar microstructure with time

    (a)Gradation index=3.2;(b)Gradation index=5.3;(c)Gradation index=8.4

    圖83種不同梯度指數(shù)的初始微結(jié)構(gòu)

    Fig.8Initial microstructures with three different gradation indexes

    圖10為3種梯度微結(jié)構(gòu)的平均晶粒尺寸和梯度指數(shù)隨模擬時間的變化曲線。由圖10a可見,平均晶粒尺寸隨演化時間的增加而增大,且初始梯度指數(shù)越大,晶粒尺寸增長越快。從圖10b可以看出,隨著演化時間的增加,3種梯度微結(jié)構(gòu)的梯度指數(shù)均逐漸減小,且減小幅度隨初始梯度指數(shù)的增大而增大。這歸因于晶粒長大過程遵循拋物線長大規(guī)律,晶粒生長率隨晶粒尺寸的增加而減小,導(dǎo)致系統(tǒng)底部的小晶粒比頂部的大尺寸晶粒長得更快,使得梯度指數(shù)減小。

    (a1)~(c1) 2×103ts;(a2)~(c2) 1×104ts;(a1),(a2) Initial gradation index=3.2;(b1),(b2) Initial gradation index=5.3;(c1),(c2) Initial gradation index=8.4

    圖93種初始梯度指數(shù)的梯度微結(jié)構(gòu)的演化過程

    Fig.9Evolution process of graded microstructures with three different initial gradation indexes

    圖10梯度微結(jié)構(gòu)的平均晶粒尺寸(a)和梯度指數(shù)(b)隨演化時間的變化

    Fig.10Changes of average grain size (a) and gradation index (b) of graded microstructure with time

    3 結(jié)論

    本研究采用相場法研究二維多晶材料中正常晶粒長大及初態(tài)分別為柱狀和梯度微結(jié)構(gòu)的晶粒長大過程。研究結(jié)果表明,晶粒長大是大晶粒吞噬相鄰小晶粒的過程,在此過程中晶界長度不斷減小,遵循能量降低原則;正常晶粒長大過程中相對晶粒尺寸分布具有時間不變性特點;柱狀微結(jié)構(gòu)的長寬比和梯度微結(jié)構(gòu)的梯度指數(shù)均直接影響晶粒長大動力學(xué),這兩種初始微結(jié)構(gòu)在演化過程中均向均勻等軸微結(jié)構(gòu)轉(zhuǎn)變;晶界曲率顯著影響晶粒長大過程,曲率越大,晶粒長大就越快。

    參考文獻(xiàn):

    [1]MEYERS M A,CHAWLA K K.Mechanical Behavior of Materials[M].Cambridge:Cambridge University Press,2009.

    [2]GEERS M G D,YVONNET J.Multiscale modeling of microstructure-property relations[J].MRS Bulletin,2016,41(8):610-616.

    [3]HOLM E A,FOILES S M.How grain growth stops:A mechanism for grain-growth stagnation in pure materials[J].Science,2010,328(5982):1138-1141.

    [4]GODFREY A W,MARTIN J W.Some Monte Carlo studies of grain growth in a temperature gradient[J].Philosophical Magazine A,1995,72(3):737-749.

    [5]WILLIAMSON A,DELPLANQUE J P.Strain energy limitations in Monte Carlo Potts modeling of grain growth[J].Computational Materials Science,2015,108:48-55.

    [6]HE Y Z,DING H L,LIU L F,et al.Computer simulation of 2D grain growth using a cellular automata model based on the lowest energy principle[J].Materials Science and Engineering:A,2006,429(1/2):236-246.

    [7]RAGHAVAN S,SAHAY S S.Modeling the topological features during grain growth by cellular automaton[J].Computational Materials Science,2009,46(1):92-99.

    [8]CHEN L Q,YANG W.Computer simulation of the domain dynamics of a quenched system with a large number of nonconserved order parameters:The grain-growth kinetics[J].Physical Review B,1994,50(21):15752-15756.

    [9]CHEN L Q,FAN D N.Computer simulation model for coupled grain growth and Ostwald Ripening-application to Al2O3-ZrO2two-phase systems[J].Journal of the American Ceramic Society,1996,79(5):1163-1168.

    [10]KAZARYAN A,WANG Y,DREGIA S A,et al.Grain growth in systems with anisotropic boundary mobility:Analytical model and computer simulation[J].Physical Review B,2001,63(18):184102.

    (Continue on page 442)(Continue from page 436)

    [11]KRILL Ⅲ C E,CHEN L Q.Computer simulation of 3-D grain growth using a phase-field model[J].Acta Materialia,2002,50(12):3059-3075.

    [12]魏承煬,李賽毅.溫度梯度對晶粒生長行為影響的相場模擬[J].物理學(xué)報,2011,60(10):132-139. WEI C Y,LI S Y.Effect of temperature gradient on grain growth behavior from phase field simulations[J].Acta Physica Sinica,2011,60(10):132-139.

    [13]KURTZ S K,CARPAY F M A.Microstructure and normal grain growth in metals and ceramics.Part Ⅱ.Experiment[J].Journal of Applied Physics,1980,51(11):5745-5754.

    [14]LOUAT N P.On the theory of normal grain growth[J].Acta Metallurgica,1974,22(6):721-724.

    [15]SAHAY S S,MALHOTRA C P,KOLKHEDE A M.Accelerated grain growth behavior during cyclic annealing[J].Acta Materialia,2003,51(2):339-346.

    [16]SAHAY S S,KUMAR B V H,KRISHNAN S J.Microstructure evolution during batch annealing[J].Journal of Materials Engineering and Performance,2003,12(6):701-707.

    (責(zé)任編輯:陸雁)

    Phase Field Study on Effect of Initial Microstructure on Grain Growth

    LUO Zhirong1,LU Chengjian1,GAO Yingjun2

    Key words:microstructure,grain growth,curvature,grain size,phase field method

    Abstract:【Objective】The effects of different initial microstructures on grain growth process and growth dynamics were studied.【Methods】The growth processes of normal grain,columnar and graded microstructures in two-dimensional polycrystalline materials were studied by phase field method.【Results】The relative normalized grain size distribution of normal grain growth has time-invariant nature,and the aspect ratio of columnar microstructure and gradient index of graded microstructure directly affect grain growth dynamics,while these two kinds of initial microstructures evolve towards uniform and equiaxed microstructures during grain growth.【Conclusion】Grain growth is a process in which large grains grow at the expense of their smaller neighbours.Grain boundary curvature has a significant influence on the grain growth process,i.e.,the larger the boundary curvature is,the faster the grain grows.

    收稿日期:2016-09-09

    作者簡介:羅志榮(1981-),男,博士,副教授,主要從事材料微結(jié)構(gòu)演化及物理性能的研究。 **通信作者:高英俊(1962-),男,教授,博士生導(dǎo)師,主要從事材料納微結(jié)構(gòu)的設(shè)計與模擬實驗研究,E-mail:gaoyj@gxu.edu.cn。

    中圖分類號:TG111

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1005-9164(2016)05-0432-05

    *國家自然科學(xué)基金項目(51561031),廣西自然科學(xué)基金項目(2015GXNSFBA139240)和廣西高??蒲许椖?YB2014318)資助。

    廣西科學(xué)Guangxi Sciences 2016,23(5):432~436,442

    網(wǎng)絡(luò)優(yōu)先數(shù)字出版時間:2016-11-21【DOI】10.13656/j.cnki.gxkx.20161121.004

    網(wǎng)絡(luò)優(yōu)先數(shù)字出版地址:http://www.cnki.net/kcms/detail/45.1206.G3.20161121.1520.008.html

    猜你喜歡
    相場柱狀微結(jié)構(gòu)
    硅片上集成高介電調(diào)諧率的柱狀納米晶BaTiO3鐵電薄膜
    基于子單元光滑有限元的混凝土相場損傷模型研究
    鑄件凝固微觀組織仿真程序開發(fā)
    基于相場理論的瀝青自愈合微觀進(jìn)程與機(jī)理研究進(jìn)展
    石油瀝青(2018年1期)2018-04-12 07:31:51
    基于COMSOL的相場模擬研究
    科技視界(2017年8期)2017-07-31 10:31:17
    在K2O—Al2O3—SiO2系致密瓷中原位合成莫來石
    江蘇陶瓷(2017年2期)2017-05-30 10:48:04
    一種柱狀金屬物質(zhì)量檢測器的研究
    電子制作(2017年10期)2017-04-18 07:23:14
    金屬微結(jié)構(gòu)電鑄裝置設(shè)計
    用于視角偏轉(zhuǎn)的光學(xué)膜表面微結(jié)構(gòu)設(shè)計
    粘結(jié)型La0.8Sr0.2MnO3/石墨復(fù)合材料的微結(jié)構(gòu)與電輸運性質(zhì)
    久久精品亚洲精品国产色婷小说| 伦理电影免费视频| 中文字幕精品亚洲无线码一区| 可以在线观看的亚洲视频| 最近在线观看免费完整版| 可以在线观看的亚洲视频| 午夜福利在线观看吧| 亚洲熟妇中文字幕五十中出| 欧美黄色片欧美黄色片| 一级a爱片免费观看的视频| 老司机在亚洲福利影院| 91老司机精品| 天堂av国产一区二区熟女人妻| netflix在线观看网站| 伦理电影免费视频| 国产av一区在线观看免费| 校园春色视频在线观看| 国产日本99.免费观看| 久久久久久大精品| 午夜免费观看网址| 国产av不卡久久| 草草在线视频免费看| 欧美又色又爽又黄视频| 日韩人妻高清精品专区| 日本五十路高清| 午夜免费成人在线视频| 欧美绝顶高潮抽搐喷水| 丁香六月欧美| 免费观看精品视频网站| 香蕉国产在线看| 亚洲色图 男人天堂 中文字幕| 最近最新免费中文字幕在线| 男女做爰动态图高潮gif福利片| 窝窝影院91人妻| 中文字幕av在线有码专区| 久久人妻av系列| 特大巨黑吊av在线直播| 亚洲成人久久爱视频| tocl精华| 国产亚洲精品久久久com| 很黄的视频免费| 麻豆成人午夜福利视频| 日韩 欧美 亚洲 中文字幕| 精品久久久久久,| 久久国产精品影院| 亚洲精品一卡2卡三卡4卡5卡| 国产综合懂色| 黑人欧美特级aaaaaa片| 老司机午夜十八禁免费视频| 男插女下体视频免费在线播放| 好男人在线观看高清免费视频| 亚洲中文av在线| 国产伦精品一区二区三区视频9 | 亚洲国产中文字幕在线视频| 老司机福利观看| 99热这里只有精品一区 | 日韩欧美在线乱码| 一级毛片精品| 窝窝影院91人妻| 欧美高清成人免费视频www| 又黄又粗又硬又大视频| 亚洲欧美精品综合一区二区三区| 国产精品乱码一区二三区的特点| 国产又黄又爽又无遮挡在线| 欧美国产日韩亚洲一区| 欧美色视频一区免费| 久久久国产成人精品二区| 青草久久国产| 99久国产av精品| 国产精品亚洲av一区麻豆| 日韩有码中文字幕| 69av精品久久久久久| 午夜福利视频1000在线观看| 狠狠狠狠99中文字幕| 亚洲精华国产精华精| 久久中文看片网| 日本精品一区二区三区蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 久久国产乱子伦精品免费另类| 国产爱豆传媒在线观看| 国产精品久久久久久久电影 | 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| 美女午夜性视频免费| 在线十欧美十亚洲十日本专区| 亚洲专区字幕在线| 曰老女人黄片| 成人鲁丝片一二三区免费| 欧美日韩乱码在线| 午夜两性在线视频| 日韩欧美国产一区二区入口| 久久久久久久久久黄片| 精品国产乱码久久久久久男人| 一本综合久久免费| 搡老熟女国产l中国老女人| 夜夜看夜夜爽夜夜摸| 村上凉子中文字幕在线| 午夜免费观看网址| 国产伦人伦偷精品视频| 午夜精品在线福利| 国产成人欧美在线观看| 免费观看精品视频网站| 黄色成人免费大全| 久久久久亚洲av毛片大全| 国产成年人精品一区二区| 国产一区二区在线av高清观看| 麻豆av在线久日| 精品欧美国产一区二区三| av片东京热男人的天堂| 男女那种视频在线观看| 一个人看视频在线观看www免费 | 真人做人爱边吃奶动态| 日本成人三级电影网站| 亚洲欧美日韩无卡精品| 一二三四社区在线视频社区8| 麻豆av在线久日| 亚洲欧美精品综合久久99| 99热精品在线国产| 在线a可以看的网站| 神马国产精品三级电影在线观看| 美女高潮喷水抽搐中文字幕| 欧美激情在线99| 无人区码免费观看不卡| 制服人妻中文乱码| 亚洲国产精品合色在线| 国产精品综合久久久久久久免费| e午夜精品久久久久久久| 欧美中文日本在线观看视频| 国产av麻豆久久久久久久| 国产黄片美女视频| 亚洲精品一区av在线观看| 床上黄色一级片| 18禁裸乳无遮挡免费网站照片| 99久国产av精品| netflix在线观看网站| 色哟哟哟哟哟哟| 99久久精品热视频| 我要搜黄色片| 欧美乱妇无乱码| 亚洲男人的天堂狠狠| 日韩国内少妇激情av| 国产爱豆传媒在线观看| 久久久久久久久久黄片| 最近最新中文字幕大全免费视频| 国产精品亚洲av一区麻豆| 免费在线观看日本一区| 小说图片视频综合网站| 五月伊人婷婷丁香| 露出奶头的视频| 欧美日韩一级在线毛片| 亚洲第一电影网av| 日本熟妇午夜| 日韩精品青青久久久久久| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影| 国产成人av教育| 夜夜躁狠狠躁天天躁| 曰老女人黄片| 在线视频色国产色| 视频区欧美日本亚洲| 俺也久久电影网| 亚洲男人的天堂狠狠| 欧美日韩国产亚洲二区| www.自偷自拍.com| 美女高潮喷水抽搐中文字幕| 午夜福利在线观看吧| 熟女少妇亚洲综合色aaa.| 久久久久久久久中文| 国产高潮美女av| 天堂动漫精品| 又大又爽又粗| 中出人妻视频一区二区| 久久久久免费精品人妻一区二区| 国产精品综合久久久久久久免费| 精品国产美女av久久久久小说| 老鸭窝网址在线观看| 91麻豆av在线| 韩国av一区二区三区四区| 日韩欧美 国产精品| 岛国视频午夜一区免费看| 国产精品 国内视频| 欧美日韩亚洲国产一区二区在线观看| 欧美高清成人免费视频www| 国产乱人视频| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| av黄色大香蕉| 亚洲精品一卡2卡三卡4卡5卡| 亚洲片人在线观看| 国产成人精品久久二区二区91| 成人三级做爰电影| 亚洲欧美日韩东京热| 美女免费视频网站| 亚洲av电影在线进入| 国产成人一区二区三区免费视频网站| 淫妇啪啪啪对白视频| 黄色成人免费大全| 精品国产美女av久久久久小说| 母亲3免费完整高清在线观看| 亚洲av电影不卡..在线观看| 老司机深夜福利视频在线观看| 一级a爱片免费观看的视频| 欧美成人性av电影在线观看| 女警被强在线播放| 美女被艹到高潮喷水动态| 精品国产乱码久久久久久男人| tocl精华| 欧美一级a爱片免费观看看| 女同久久另类99精品国产91| 久久国产乱子伦精品免费另类| 亚洲午夜精品一区,二区,三区| 啦啦啦免费观看视频1| 九九在线视频观看精品| 日本 av在线| 日本 欧美在线| 亚洲第一电影网av| 麻豆成人av在线观看| 日韩欧美三级三区| 免费观看的影片在线观看| 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 精品国产三级普通话版| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 成人国产综合亚洲| av欧美777| 国产成人精品久久二区二区91| 美女扒开内裤让男人捅视频| 亚洲专区国产一区二区| 国产私拍福利视频在线观看| 嫁个100分男人电影在线观看| 国内精品久久久久精免费| 麻豆久久精品国产亚洲av| 床上黄色一级片| 免费在线观看亚洲国产| 国产精品亚洲美女久久久| 美女午夜性视频免费| a级毛片在线看网站| 亚洲天堂国产精品一区在线| 999久久久精品免费观看国产| 亚洲国产欧洲综合997久久,| 99精品久久久久人妻精品| 99久久99久久久精品蜜桃| 日韩欧美精品v在线| 69av精品久久久久久| 国产人伦9x9x在线观看| 亚洲国产高清在线一区二区三| 国产av不卡久久| 99久久国产精品久久久| 久久精品亚洲精品国产色婷小说| 亚洲熟女毛片儿| 美女免费视频网站| 精品久久久久久久末码| xxx96com| 国内少妇人妻偷人精品xxx网站 | 亚洲成av人片在线播放无| 村上凉子中文字幕在线| 亚洲av熟女| 亚洲国产精品国产精品| 国产精品一区二区性色av| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 国产老妇女一区| 亚洲最大成人av| 久久精品国产99精品国产亚洲性色| 99视频精品全部免费 在线| 国产高清国产精品国产三级 | 如何舔出高潮| 亚洲中文字幕一区二区三区有码在线看| 日韩国内少妇激情av| eeuss影院久久| 插阴视频在线观看视频| 99久久人妻综合| 精品国产一区二区三区久久久樱花 | 欧美激情国产日韩精品一区| 黄片wwwwww| 亚洲国产最新在线播放| 欧美潮喷喷水| 成人三级黄色视频| 午夜免费激情av| 亚洲av中文字字幕乱码综合| 中文欧美无线码| 日韩欧美精品免费久久| 国产黄色视频一区二区在线观看 | 亚洲精品色激情综合| 久99久视频精品免费| 尾随美女入室| 亚洲伊人久久精品综合 | 日韩在线高清观看一区二区三区| 97人妻精品一区二区三区麻豆| 一区二区三区免费毛片| av专区在线播放| 国产真实伦视频高清在线观看| 美女xxoo啪啪120秒动态图| 国产高清有码在线观看视频| 亚洲国产精品国产精品| 久久欧美精品欧美久久欧美| 亚洲高清免费不卡视频| 成人美女网站在线观看视频| 国产精品伦人一区二区| 一区二区三区免费毛片| 国产成人a∨麻豆精品| av.在线天堂| 欧美潮喷喷水| 亚洲成人av在线免费| 女人久久www免费人成看片 | 乱码一卡2卡4卡精品| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 亚洲不卡免费看| 我的老师免费观看完整版| av在线天堂中文字幕| 内地一区二区视频在线| 成人毛片60女人毛片免费| 一级二级三级毛片免费看| 国产久久久一区二区三区| 我的老师免费观看完整版| 在线观看美女被高潮喷水网站| 国产乱人偷精品视频| 91久久精品电影网| 日本av手机在线免费观看| 免费黄色在线免费观看| 看片在线看免费视频| 欧美最新免费一区二区三区| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| .国产精品久久| 久久久久久久午夜电影| 99热全是精品| av天堂中文字幕网| 久久国产乱子免费精品| 精品久久久久久久久亚洲| 18禁在线无遮挡免费观看视频| 日韩一区二区视频免费看| 夜夜爽夜夜爽视频| 汤姆久久久久久久影院中文字幕 | 大话2 男鬼变身卡| 午夜a级毛片| 亚洲欧美中文字幕日韩二区| 精品酒店卫生间| 午夜福利在线观看免费完整高清在| 草草在线视频免费看| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 欧美不卡视频在线免费观看| 国产精品无大码| 成人无遮挡网站| 欧美成人精品欧美一级黄| 欧美三级亚洲精品| 少妇丰满av| 人妻少妇偷人精品九色| 国产高清视频在线观看网站| 中文字幕精品亚洲无线码一区| 一个人看的www免费观看视频| 色综合站精品国产| 国产单亲对白刺激| 欧美最新免费一区二区三区| 国产av不卡久久| 精品无人区乱码1区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产成人免费观看mmmm| 国产一区二区亚洲精品在线观看| 亚洲精品aⅴ在线观看| 精品人妻偷拍中文字幕| 色哟哟·www| 高清毛片免费看| 波多野结衣高清无吗| 一本久久精品| 精品少妇黑人巨大在线播放 | 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 国产精品美女特级片免费视频播放器| 久久久亚洲精品成人影院| 男女视频在线观看网站免费| 免费黄网站久久成人精品| 亚洲欧美精品自产自拍| 99九九线精品视频在线观看视频| 最近手机中文字幕大全| 国产色爽女视频免费观看| 亚洲自拍偷在线| 99久久精品热视频| 99热这里只有是精品在线观看| 国产精品蜜桃在线观看| 日韩人妻高清精品专区| 18禁在线播放成人免费| 午夜精品国产一区二区电影 | 国产精品人妻久久久影院| 熟女人妻精品中文字幕| 国产精品久久视频播放| 国产伦理片在线播放av一区| 亚洲国产精品专区欧美| 成年免费大片在线观看| 九草在线视频观看| 美女黄网站色视频| 亚洲欧美精品综合久久99| 在线a可以看的网站| 欧美一区二区亚洲| 啦啦啦观看免费观看视频高清| 免费人成在线观看视频色| 精品久久久久久电影网 | 又爽又黄a免费视频| 精品人妻视频免费看| 一区二区三区四区激情视频| 日韩欧美在线乱码| 九九久久精品国产亚洲av麻豆| 亚洲国产精品成人综合色| 久热久热在线精品观看| 精品人妻偷拍中文字幕| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 在线免费观看的www视频| 欧美成人免费av一区二区三区| 亚洲经典国产精华液单| 黑人高潮一二区| www日本黄色视频网| 三级毛片av免费| 天堂中文最新版在线下载 | 亚洲精品国产av成人精品| 白带黄色成豆腐渣| 有码 亚洲区| 亚洲五月天丁香| 成人高潮视频无遮挡免费网站| 一级毛片电影观看 | 中文字幕制服av| 国产极品天堂在线| 一个人观看的视频www高清免费观看| 国产精品麻豆人妻色哟哟久久 | 人妻夜夜爽99麻豆av| 久久精品影院6| 欧美色视频一区免费| 精品国产一区二区三区久久久樱花 | 久久久久久久久久久丰满| 高清午夜精品一区二区三区| 久久99精品国语久久久| 能在线免费观看的黄片| 黄色日韩在线| 一夜夜www| kizo精华| 永久网站在线| 精品久久久久久久久av| 久久精品国产99精品国产亚洲性色| 国产日韩欧美在线精品| 高清日韩中文字幕在线| 一边摸一边抽搐一进一小说| 久久久久久久久久久免费av| or卡值多少钱| 久久精品夜色国产| 精品午夜福利在线看| 国模一区二区三区四区视频| 日韩,欧美,国产一区二区三区 | 毛片一级片免费看久久久久| 亚洲aⅴ乱码一区二区在线播放| 成人午夜精彩视频在线观看| 又粗又硬又长又爽又黄的视频| 国产91av在线免费观看| 亚洲精品久久久久久婷婷小说 | 国产精品电影一区二区三区| 26uuu在线亚洲综合色| 天堂√8在线中文| 男女那种视频在线观看| 久久99蜜桃精品久久| 久久久精品欧美日韩精品| av专区在线播放| 又粗又爽又猛毛片免费看| 国产人妻一区二区三区在| 国产高清不卡午夜福利| 久久精品人妻少妇| 两个人视频免费观看高清| 熟妇人妻久久中文字幕3abv| 黄片wwwwww| 欧美成人免费av一区二区三区| 国产白丝娇喘喷水9色精品| 九色成人免费人妻av| 中文字幕久久专区| 亚洲精品日韩在线中文字幕| 简卡轻食公司| 欧美精品一区二区大全| 亚洲欧洲日产国产| 国产精品电影一区二区三区| ponron亚洲| 最近的中文字幕免费完整| 小说图片视频综合网站| 丝袜美腿在线中文| 欧美成人午夜免费资源| 国产一区二区在线观看日韩| 国产69精品久久久久777片| 不卡视频在线观看欧美| 日本av手机在线免费观看| 国产成人a区在线观看| 亚洲av成人精品一二三区| 国产精品一区www在线观看| 永久免费av网站大全| 午夜亚洲福利在线播放| 国产伦在线观看视频一区| 久久这里只有精品中国| 69av精品久久久久久| 亚洲av二区三区四区| 精品久久久久久久末码| 久久久久久久久久成人| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 国产视频首页在线观看| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 最近手机中文字幕大全| 日本免费在线观看一区| 狂野欧美激情性xxxx在线观看| 国产色爽女视频免费观看| 五月伊人婷婷丁香| 女人被狂操c到高潮| 夜夜看夜夜爽夜夜摸| 一夜夜www| 国产单亲对白刺激| 黄片无遮挡物在线观看| 好男人在线观看高清免费视频| 亚洲精品aⅴ在线观看| 久99久视频精品免费| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 中文字幕亚洲精品专区| 欧美最新免费一区二区三区| 丰满少妇做爰视频| 久久精品91蜜桃| 高清日韩中文字幕在线| 18禁动态无遮挡网站| 亚洲精华国产精华液的使用体验| 国产成人免费观看mmmm| 午夜精品国产一区二区电影 | 久久久久久大精品| 天堂网av新在线| 久久久久久久午夜电影| 久久人人爽人人片av| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影小说 | 女人被狂操c到高潮| 国产色婷婷99| 亚洲国产欧美人成| 午夜精品一区二区三区免费看| 国产国拍精品亚洲av在线观看| 中文字幕制服av| 亚洲av成人精品一二三区| 色播亚洲综合网| 人体艺术视频欧美日本| 18禁在线无遮挡免费观看视频| 男人狂女人下面高潮的视频| 久久久久精品久久久久真实原创| 亚洲一区高清亚洲精品| 人体艺术视频欧美日本| 国产老妇伦熟女老妇高清| av天堂中文字幕网| 日韩欧美精品免费久久| 亚洲第一区二区三区不卡| 久99久视频精品免费| 亚洲三级黄色毛片| 日韩强制内射视频| 成人综合一区亚洲| 亚洲精品国产成人久久av| 成人性生交大片免费视频hd| 国产亚洲精品av在线| 久久综合国产亚洲精品| 熟妇人妻久久中文字幕3abv| 免费大片18禁| 免费黄网站久久成人精品| 日本午夜av视频| 别揉我奶头 嗯啊视频| av在线蜜桃| 久久精品国产亚洲网站| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| 内射极品少妇av片p| 好男人视频免费观看在线| 纵有疾风起免费观看全集完整版 | 看十八女毛片水多多多| 禁无遮挡网站| 国产免费一级a男人的天堂| 久久99热6这里只有精品| 日韩中字成人| 久久久国产成人精品二区| 国产精品无大码| 成年av动漫网址| 五月玫瑰六月丁香| 黄色欧美视频在线观看| 亚洲成人精品中文字幕电影| av线在线观看网站| 国产精品一区二区三区四区免费观看| 国语对白做爰xxxⅹ性视频网站| av天堂中文字幕网| 亚洲自偷自拍三级| kizo精华| 中国国产av一级| 日本熟妇午夜| 久久6这里有精品| 国产极品精品免费视频能看的| 七月丁香在线播放| 精品久久久久久久人妻蜜臀av| 亚洲熟妇中文字幕五十中出| 免费搜索国产男女视频| 夫妻性生交免费视频一级片| 18禁在线无遮挡免费观看视频| 永久免费av网站大全| 校园人妻丝袜中文字幕| 99视频精品全部免费 在线| 色播亚洲综合网| 夜夜爽夜夜爽视频| 综合色丁香网| 国产成人一区二区在线| 亚洲高清免费不卡视频| 三级经典国产精品| 国产精品久久久久久精品电影小说 | 1024手机看黄色片| 久久99蜜桃精品久久| 亚洲欧美精品综合久久99| 欧美最新免费一区二区三区| 久99久视频精品免费|