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

    納米級(jí)缺口起裂臨界條件的晶體相場(chǎng)模擬*

    2017-01-03 02:45:38溫振川寧歲婷高英俊
    廣西科學(xué) 2016年5期
    關(guān)鍵詞:相場(chǎng)缺口步長(zhǎng)

    溫振川,寧歲婷,張 睿,高英俊

    (廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院,廣西南寧 530004)

    ?

    納米級(jí)缺口起裂臨界條件的晶體相場(chǎng)模擬*

    溫振川,寧歲婷,張睿,高英俊**

    (廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院,廣西南寧530004)

    (College of Physical Science and Technology,Guangxi University,Nanning,Guangxi,530004,China)

    摘要:【目的】研究初始缺口大小對(duì)裂紋起裂的影響?!痉椒ā坎捎镁w相場(chǎng)(phase-field-crystal,PFC)模型研究單軸拉伸應(yīng)力作用下,材料裂紋起裂的動(dòng)態(tài)演化過(guò)程,并結(jié)合自由能分析裂紋起裂與擴(kuò)展過(guò)程中,不同缺口尺寸對(duì)裂紋起裂的影響。【結(jié)果】沿x方向施加恒應(yīng)變,初始缺口y方向尺寸較小時(shí),改變?nèi)笨趚方向尺寸,缺口并不起裂;當(dāng)y方向尺寸達(dá)到一定值時(shí),初始缺口會(huì)起裂。當(dāng)y方向尺寸接近起裂臨界值時(shí),x方向尺寸的增大會(huì)延緩缺口的起裂。當(dāng)x方向尺寸改變量較小時(shí),裂紋并不起裂,當(dāng)增大到x方向尺寸臨界值時(shí),缺口開(kāi)始起裂并擴(kuò)展?!窘Y(jié)論】在應(yīng)力大小一定的情況下,缺口起裂與應(yīng)力作用方向垂直的缺口尺寸密切相關(guān),只有與應(yīng)變作用方向垂直的尺寸達(dá)到臨界值,缺口才會(huì)發(fā)生起裂。

    關(guān)鍵詞:晶體相場(chǎng)模型臨界缺口尺寸缺口起裂

    0 引言

    【研究意義】材料在生產(chǎn)制造和實(shí)際應(yīng)用的過(guò)程中,受加載作用下容易產(chǎn)生微裂紋,這些微裂紋在一定的條件下會(huì)發(fā)生擴(kuò)展,從而影響材料的壽命[1]。早期的研究是從材料斷裂的宏觀表象出發(fā),建立了斷裂力學(xué)理論和模型,來(lái)解釋裂紋萌生和擴(kuò)展的相關(guān)行為[2]。但實(shí)際上,裂紋的萌生與擴(kuò)展[3]在納米尺度就開(kāi)始了。因此,深入研究材料的斷裂特征,就不能只停留在給出表象的宏觀尺度上,而應(yīng)當(dāng)從納觀尺度分析裂紋萌生與擴(kuò)展機(jī)理?!厩叭搜芯窟M(jìn)展】材料在介觀和宏觀尺度上的性能,在很大程度上由其內(nèi)部的微觀缺陷所決定[4],例如:空洞、空位、位錯(cuò)、晶界以及微裂紋等。這些缺陷起因于在原子尺度發(fā)生的復(fù)雜非平衡動(dòng)力學(xué)過(guò)程[5]。就目前所采用的多尺度模擬研究方法來(lái)說(shuō),對(duì)這些微觀特征結(jié)構(gòu)的實(shí)時(shí)演化過(guò)程的建模和模擬是一個(gè)重大的挑戰(zhàn)。最近,Elder等[5-7]基于密度泛函理論提出了晶體相場(chǎng)(phase-field-crystal,PFC)模型。該模型給出了新的自由能函數(shù)形式,引入周期序參量作為局域密度場(chǎng),它將無(wú)序相(液相等)的密度場(chǎng)定義為一常數(shù),將固相的密度場(chǎng)表示成周期性函數(shù)(波)的形式,進(jìn)而通過(guò)周期的原子密度函數(shù)表現(xiàn)晶體的周期結(jié)構(gòu)。這樣的密度場(chǎng)就很自然地與彈性效應(yīng)、晶粒取向和位錯(cuò)的運(yùn)動(dòng)等由周期性結(jié)構(gòu)產(chǎn)生的物理特性緊密的關(guān)聯(lián)起來(lái)[8]。目前,PFC方法已有許多具體的應(yīng)用,例如:模擬位錯(cuò)滑移和攀移[9]、相結(jié)構(gòu)轉(zhuǎn)變[10]、異質(zhì)外延生長(zhǎng)[11]、納米晶粒長(zhǎng)大[12]等?!颈狙芯壳腥朦c(diǎn)】材料微裂紋擴(kuò)展的研究目前主要集中在材料受到外部載荷累計(jì)改變作用下的疲勞裂紋尖端擴(kuò)展行為以及裂紋分叉特性[13-14],而應(yīng)用先進(jìn)的PFC模型,研究初始缺口尺寸大小對(duì)裂紋起裂影響的工作還鮮有報(bào)道?!緮M解決的關(guān)鍵問(wèn)題】應(yīng)用PFC方法模擬拉伸應(yīng)變作用下[15],納米級(jí)裂紋起裂的動(dòng)態(tài)演化過(guò)程,詳細(xì)分析預(yù)制中心缺口尺寸大小對(duì)裂紋起裂的影響。

    1 PFC模型與方法

    1.1體系的自由能密度函數(shù)

    與傳統(tǒng)的相場(chǎng)模型不同,PFC模型采用具有周期性的局域原子密度場(chǎng)作為相場(chǎng)變量[6]。因此,PFC模型能夠反映晶體周期結(jié)構(gòu)特性以及原子尺度的行為。對(duì)于固態(tài)金屬材料,其原子的位置不依賴(lài)于時(shí)間,呈現(xiàn)規(guī)則排列。文獻(xiàn)[6]給出了既能夠描述晶體周期結(jié)構(gòu),又能描述體系均勻結(jié)構(gòu)的密度函數(shù),其表達(dá)式可寫(xiě)成

    (1)

    式中,G為倒格子,r為空間位置矢量,an,m為Fourier系數(shù)。式(1)中等號(hào)右邊第一項(xiàng)反映晶格原子的周期性排列特征;第二項(xiàng)ρ0反映均勻無(wú)序相的原子密度分布平均值,是一個(gè)常量。此時(shí)體系無(wú)量綱的自由能函數(shù)F可以寫(xiě)成

    (2)

    (3)

    式中:A0是一個(gè)特定的常數(shù),反映固相原子密度周期結(jié)構(gòu)的振幅。由能量函數(shù)取極小值得到:波矢q=2π/a,A0表達(dá)式為

    (4)

    式中,ρ0和γ為體系能量函數(shù)的兩個(gè)重要參數(shù)。由體系的能量密度函數(shù)可以計(jì)算并畫(huà)出體系的相圖。對(duì)于二維體系,相區(qū)有均勻無(wú)序相和晶體相,且晶體相又有六角格子相和條狀相兩種。文獻(xiàn)[6]已給出均勻相和條狀相的能量密度函數(shù)表達(dá)式,按照平衡相圖的計(jì)算方法可得到二維PFC相圖如圖1所示。

    S、T、L分別表示條狀相、三角相、液相;圖(b)為圖(a)中虛線(xiàn)框的放大圖

    S,T and L represent stripe phase,triangular phase and liquid phase,respectively;(b)is enlarged image of the dashed box of (a)

    圖1PFC模型的二維相圖[6]

    Fig.1Two-dimensional phase diagrams of PFC model[6]

    1.2動(dòng)力學(xué)方程與數(shù)值化處理

    保守的原子密度場(chǎng)變量的演化可用與時(shí)間相關(guān)的Cahn-Hilliard動(dòng)力學(xué)方程[6]描述

    (5)

    該動(dòng)力學(xué)方程是六階非線(xiàn)性高階偏微分方程,在坐標(biāo)空間求解十分復(fù)雜和繁瑣。因此,為了提高求解動(dòng)力學(xué)方程(5)的效率,需要轉(zhuǎn)換到波矢空間求解。這里采用半隱式傅里葉法[6,17]求解,對(duì)式(5)變換到傅里葉空間[18],整理可得

    (6)

    1.3拉伸應(yīng)變的施加以及參數(shù)設(shè)置

    對(duì)樣品采用等應(yīng)變(等位移)單軸拉伸加載。在受到拉伸作用后,樣品中所有格點(diǎn)的坐標(biāo)位置都會(huì)發(fā)生相應(yīng)的變化。因此,通過(guò)格點(diǎn)坐標(biāo)位置的變化來(lái)反映應(yīng)變作用所引起的樣品變形。設(shè)變形前的樣品中格點(diǎn)位置的坐標(biāo)用(x,y)表示,施加等應(yīng)變作用引起變形后,格點(diǎn)位置的坐標(biāo)用(x′,y′)表示。樣品在y方向受到應(yīng)變?chǔ)舮的拉伸作用。變形前后坐標(biāo)變換關(guān)系為x′=x和y′=y/(1+εy)。由于等應(yīng)變作用引起樣品變形,其內(nèi)部原子密度函數(shù)變?yōu)棣?x′,y′)=ρ(x,y/(1+εy)),因此,在考慮等應(yīng)變作用引起變形的情況下,原子密度函數(shù)表達(dá)式(4)將變換為ρ(x′,y′),體系的能量也隨之發(fā)生相應(yīng)的變化,但控制演化的動(dòng)力學(xué)方程的形式保持不變。將樣品變形后的原子密度ρ(x′,y′)代入動(dòng)力學(xué)方程(5)中,采用半隱式傅里葉變換方法整理得到(6),對(duì)其求解并進(jìn)行可視化處理,即可得到裂紋擴(kuò)展的演化圖。

    本研究模擬不涉及具體材料的物性參數(shù),所用參數(shù)均已無(wú)量綱化處理,并將連續(xù)空間離散為四方格子網(wǎng)格,采用周期性邊界條件。計(jì)算網(wǎng)格為512Δx×512Δy,空間步長(zhǎng)Δx=Δy=π/3,時(shí)間步長(zhǎng)Δt=0.05。根據(jù)單模近似下的二維相圖(圖1),本研究選用無(wú)量綱密度為ρ0=0.49的三角相表征晶體相,對(duì)應(yīng)的溫度參量取為γ=-1.0,如圖1b中A點(diǎn)所示。為模擬中心裂紋,在樣品中心設(shè)置一定大小的缺口作為初始裂口,其參數(shù)設(shè)置為ρ0=0.79,γ=-1.0,如圖1b中B點(diǎn)所示。之所以選取該點(diǎn)是因?yàn)槠鋵?duì)應(yīng)液相與三角相共存區(qū)域附近的原子密度,有利于三角相向裂紋結(jié)構(gòu)轉(zhuǎn)變。原子排列取向角θ=1°。

    2 結(jié)果與分析

    2.1缺口x方向起裂的臨界尺寸

    在施加相同應(yīng)變情況下,不同的初始缺口尺寸會(huì)影響裂紋的萌生與擴(kuò)展。如圖2所示,本研究施加拉應(yīng)變方向?yàn)閱蝬方向,應(yīng)變量為ε=10%。固定缺口y方向長(zhǎng)度為10個(gè)格子,改變?nèi)笨趚方向的尺寸,分別取10,12,14,16,18,20,22個(gè)格子(如表1樣品Ai所示)。圖3a~d為實(shí)驗(yàn)條件下,樣品A6缺口形貌圖,從圖3可以看出,整個(gè)過(guò)程中初始缺口沒(méi)有起裂。圖4給出了不同尺寸初始缺口樣品的自由能曲線(xiàn)。由于采取恒應(yīng)變施加的方法,從該圖中可觀察到,在整個(gè)過(guò)程,自由能曲線(xiàn)呈水平直線(xiàn),體系能量保持常量,這表明裂紋沒(méi)有起裂。由此可以推斷,缺口y方向長(zhǎng)度固定且較小時(shí),改變?nèi)笨趚方向的尺寸,不產(chǎn)生缺口起裂;但是隨著x方向尺寸的增大,同等應(yīng)變作用下,整個(gè)體系的自由能升高。

    圖2施加單軸拉伸應(yīng)變的示意圖

    Fig.2Effect diagram of single biaxial stress

    表1樣品參數(shù)及起裂情況

    Table 1Sample parameters and crack initiation

    樣品Sample缺口尺寸NotchsizeΔx×Δy起裂CrackinitiationA110×10NA212×10NA314×10NA416×10NA518×10NA620×10NA722×10NB110×10NB210×12NB310×14NB410×16YB510×18YB610×20YB710×22YC110×14NC212×14YC314×14YC416×14YC518×14YC620×14YC722×14Y

    注:N表示沒(méi)有起裂,Y表示發(fā)生起裂

    Note:N means no crack initiation,Y means crack initiation

    t表示時(shí)間步長(zhǎng)t is time step size圖3 樣品A6缺口形貌Fig 3 The morphology of crack of sample A6

    圖4樣品Ai體系的自由能變化曲線(xiàn)

    Fig.4Free energy change curve of the evolution of sample Ai

    2.2缺口y方向起裂的臨界尺寸

    固定缺口x方向尺寸為10個(gè)格子,改變初始缺口y方向的尺寸,分別取10,12,14,16,18,20,22個(gè)格子(如表1樣品Bi所示)。圖5給出了不同尺寸缺口樣品的自由能曲線(xiàn),從圖中可以觀察到,當(dāng)y方向缺口尺寸小于等于14個(gè)格子時(shí),自由能幾乎是一條水平的直線(xiàn),反映了整個(gè)體系能量基本沒(méi)有變化,也表明缺口沒(méi)有起裂。但隨著y方向缺口尺寸增大,不同樣品能量升高;當(dāng)y方向尺寸大于14個(gè)格子時(shí),自由能曲線(xiàn)在時(shí)間步長(zhǎng)大約為25萬(wàn)的時(shí)候開(kāi)始降低,這意味著缺口在該時(shí)刻開(kāi)始起裂,體系自由能得到釋放。對(duì)于y方向缺口尺寸不同的起裂樣品,起裂時(shí)間與曲線(xiàn)降低速率變化相一致。圖6為樣品B6在x方向的應(yīng)變10%作用下的裂口起裂擴(kuò)展演化圖。從圖6a~b中能看出,缺口沒(méi)有發(fā)生明顯的變化,圖6c中缺口出現(xiàn)明顯的細(xì)小裂縫,對(duì)應(yīng)圖4自由能曲線(xiàn)能量開(kāi)始降低的時(shí)刻,說(shuō)明缺口已經(jīng)開(kāi)始起裂,是由于起裂擴(kuò)展過(guò)程中釋放應(yīng)變能,導(dǎo)致自由能曲線(xiàn)降低。圖7顯示了樣品B1~B7在時(shí)間步長(zhǎng)40萬(wàn)步時(shí)的裂紋形貌,圖7a~c缺口隨著選取缺口參數(shù)的不同而有細(xì)微的變化,但無(wú)明顯起裂情況,符合圖4中自由能曲線(xiàn)1,2和3對(duì)應(yīng)的能量變化規(guī)律;圖7d~g這4個(gè)形貌圖中,裂紋的形貌基本一致,符合圖4中自由能曲線(xiàn)4,5,6和7對(duì)應(yīng)的能量變化規(guī)律。該結(jié)果與格里菲斯強(qiáng)度理論[19]基本一致,即裂紋擴(kuò)展的應(yīng)變大小與垂直應(yīng)變作用方向的缺口尺寸密切相關(guān)[20]。

    由2.1,2.2的結(jié)果可推斷,在單x方向施加固定的應(yīng)變量作用下,裂紋擴(kuò)展主要是受初始缺口y方向尺寸變化的影響。只有y方向缺口尺寸達(dá)到臨界值,缺口才會(huì)起裂。

    圖5 樣品Bi裂紋演化過(guò)程的體系自由能曲線(xiàn)Fig.5 Free energy curve of sample Bi

    t表示時(shí)間步長(zhǎng)

    t is time step size

    圖6x方向應(yīng)變10%作用下樣品B6的裂紋擴(kuò)展演化

    Fig.6The morphology evolution of crack propagation of sample B6under 10% of strain at x direction

    圖7時(shí)間步長(zhǎng)40萬(wàn)步時(shí)樣品B1~B7的裂紋形貌

    Fig.7The morphology of the sample B1~B7at the time step of 400 thousands

    2.3缺口y方向臨界尺寸與起裂的關(guān)系

    設(shè)置初始缺口y方向尺寸為14個(gè)格子,x方向尺寸分別取10,12,14,16,18,20,22個(gè)格子(如表1樣品Ci所示),研究在單x方向施加固定的應(yīng)變量作用下,缺口x方向尺寸變化是否對(duì)樣品裂紋擴(kuò)展情況產(chǎn)生影響。圖8給出了不同x方向尺寸樣品裂紋演化的自由能曲線(xiàn)。從圖中可以觀察到,當(dāng)x方向尺寸取10個(gè)格子時(shí),曲線(xiàn)是一條比較平滑的直線(xiàn),這意味著整個(gè)樣品體系的能量幾乎沒(méi)有變化,裂紋沒(méi)有起裂;當(dāng)x方向尺寸取大于10個(gè)格子時(shí),自由能曲線(xiàn)均有明顯的下滑過(guò)程,能量顯著降低,但曲線(xiàn)走向并不是在同一時(shí)刻開(kāi)始降低,這表明不同的x尺寸對(duì)應(yīng)的能量開(kāi)始降低的起始步長(zhǎng)不同。隨著缺口x方向尺寸增大,能量降低時(shí)間不斷延緩,缺口起裂時(shí)間不斷滯后。由此可知,當(dāng)缺口y方向尺寸達(dá)到裂紋擴(kuò)展臨界尺寸值時(shí),x方向尺寸增加,會(huì)延緩缺口的起裂。

    圖8樣品Ci裂紋演化過(guò)程的體系自由能曲線(xiàn)

    Fig.8Free energy curve of evolution of crack propagation of sample Ci

    從圖9中可以看出,當(dāng)時(shí)間步長(zhǎng)t達(dá)到15萬(wàn)步(圖9b)時(shí),初始缺口還沒(méi)有出現(xiàn)明顯的擴(kuò)展,缺口邊緣已經(jīng)出現(xiàn)了原子層局部滑移。當(dāng)時(shí)間步長(zhǎng)t達(dá)到25萬(wàn)步(圖9c)時(shí),缺口已經(jīng)出現(xiàn)了明顯的裂紋擴(kuò)展,并且隨著時(shí)間步長(zhǎng)的增加,裂紋擴(kuò)展越明顯。圖10顯示了樣品C1~C7在時(shí)間步長(zhǎng)40萬(wàn)步時(shí)的形貌,圖10a對(duì)應(yīng)圖8中的自由能曲線(xiàn)1,整個(gè)過(guò)程初始裂口沒(méi)有起裂;圖10b和圖10c中樣品的裂紋形態(tài)基本一致,分別對(duì)應(yīng)圖8中的自由能曲線(xiàn)2和3;圖10d~f中3個(gè)樣品的裂紋形態(tài)一致,但是不同于圖10b~c樣品裂紋擴(kuò)展情況,圖10d,10e和10f中的樣品分別對(duì)應(yīng)圖8中的自由能曲線(xiàn)4,5和6;圖10g裂紋擴(kuò)展情況落后于圖10其他已經(jīng)起裂樣品,與圖8中其對(duì)應(yīng)的自由能曲線(xiàn)7降低的拐點(diǎn)落后于其他自由能曲線(xiàn)的情況相符。由此可推斷,在施加單方向恒應(yīng)變及缺口y方向尺寸達(dá)到裂紋擴(kuò)展臨界尺寸值的條件下,x方向尺寸改變量較小時(shí),裂紋并不起裂;當(dāng)增大到臨界值時(shí),缺口開(kāi)始起裂并擴(kuò)展。

    圖9單x方向10%應(yīng)變作用下樣品C6的裂紋擴(kuò)展演化

    Fig.9The morphology evolution of crack propagation of sample C6under 10% of strain in uniaxial x direction

    圖10時(shí)間步長(zhǎng)為40萬(wàn)步時(shí)樣品C1~C7的形貌

    Fig.10The morphology of the sample C1~C7at the time step of 400 thousands

    3 結(jié)論

    在沿x方向施加恒定10%應(yīng)變量條件下,缺口y方向尺寸遠(yuǎn)小于臨界值時(shí),改變?nèi)笨趚方向尺寸并使其增大,缺口不起裂;固定缺口x方向尺寸,改變?nèi)笨趛方向尺寸并使其增大,當(dāng)缺口y方向尺寸增大到臨界值時(shí),開(kāi)始起裂并擴(kuò)展;當(dāng)缺口y方向尺寸接近臨界值時(shí),則改變?nèi)笨趚方向尺寸時(shí),當(dāng)x方向尺寸較小時(shí),裂紋并不起裂,隨著x增大,當(dāng)達(dá)到x尺寸臨界值時(shí),缺口開(kāi)始起裂并擴(kuò)展。表明在施加恒應(yīng)變作用下,缺口起裂與應(yīng)力作用方向垂直的缺口尺寸密切相關(guān),只有與應(yīng)變作用方向垂直的尺寸達(dá)到臨界值,缺口才會(huì)發(fā)生起裂。

    參考文獻(xiàn):

    [1]衣林,陳躍良,徐麗,等.金屬材料疲勞微裂紋的萌生與擴(kuò)展研究[J].飛機(jī)設(shè)計(jì),2012,32(2):63-67. YI L,CHEN Y L,XU L,et al.Research on initiation and propagation of fatigue micro-crack for metal materials[J].Aircraft Design,2012,32(2):63-67.

    [2]鐘群鵬,周煜,張崢.裂紋學(xué)[M].北京:高等教育出版社,2015. ZHONG Q P,ZHOU Y,ZHANG Z.Crack Science[M].Beijing:Higher Education Press,2015.

    [3]高英俊,羅志榮,黃禮琳,等.韌性材料的微裂紋擴(kuò)展和連通的晶體相場(chǎng)模擬[J].中國(guó)有色金屬學(xué)報(bào),2013,23(07):1892-1899. GAO Y J,LUO Z R,HUANG L L,et al.Phase field crastal modeling for microcrack propagation and connecting of ductile materials[J].The Chinese Journal of Nonferrous Metals,2013,23(07):1892-1899.

    [4]劉曉波,徐慶軍,劉劍.鋁裂紋擴(kuò)展行為的分子動(dòng)力學(xué)模擬[J].中國(guó)有色金屬學(xué)報(bào),2014,24(6):08. LIU X B,XU Q B,LIU J.Molecular dynamics simulation of crack propagation behavior of aluminum[J].The Chinese Journal of Nonferrous Metals,2014,24(6):08.

    [5]ELDER K R,KATAKOWSKI M,HAATAJA M,et al.Modeling elasticity in crystal growth[J].Physical Review Letters,2002,88(24):245701.

    [6]ELDER K R,GRANT M.Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals[J].Physical Review E,2004,70(5Pt1):051605.

    [7]ELDER K R,PROVATAS N,BERRY J,et al.Phase-field crystal modeling and classical density functional theory of freezing[J].Physical Review B,2007,75(6):064107.

    [8]BERRY J,GRANT M,ELDER K R.Diffusive atomistic dynamics of edge dislocations in two dimensions[J].Physical Review E,2006,73(3Pt1):031609.

    [9]楊濤,陳錚,董衛(wèi)平.應(yīng)力誘發(fā)雙位錯(cuò)組亞晶界湮沒(méi)的晶體相場(chǎng)模擬[J].金屬學(xué)報(bào),2011,47(10):1301-1306. YANG T,CHEN Z,DONG W P,et al.Phase field crystal simulation of stress-induced annihilation of sub-grain boundary with double-array dislocation[J].Acta Metallurgica Sinica,2011,47(10):1301-1306.

    [10]羅志榮,黃世葉,茹謝辛,等.晶體相場(chǎng)法模擬大角度晶界的變形過(guò)程[J].廣西科學(xué),2013,20(4):311-315. LUO Z R,HUANG S Y,RU X X,et al.Phase-field-crystal modeling for deformation process of high-angle grain boundary[J].Guangxi Sciences,2013,20(4):311-315.

    [11]高英俊,羅志榮,黃創(chuàng)高,等.晶體相場(chǎng)方法研究二維六角相向正方相結(jié)構(gòu)轉(zhuǎn)變[J].物理學(xué)報(bào),2013,62(5):050507-050510. GAO Y J,LUO Z R,HUANG C G,et al.Phase-field-crystal modeling for two-dimensional transformation from hexagonal to square structure[J].Acta Physica Sinica,2013,62(5):050507.

    [12]YU Y M,BACKOFEN R,VOIGT A.Morphological instability of heteroepitaxial growth on vicinal substrates:A phase-field crystal study[J].Journal of Crystal Growth,2011,318(1):18-22.

    [13]高英俊,王江帆,羅志榮,等.晶體相場(chǎng)方法模擬納米孿晶結(jié)構(gòu)[J].計(jì)算物理,2013,30(4):577-581. GAO Y J,WANG J F,LUO Z R,et al.Nano-twin structure simulation with phase field crystal method[J].Chinese Journal of Computational Physics,2013,30(4):577-581.

    [14]HU S,CHEN Z,PENG Y Y,et al.Modeling and simulation of microcrack propagation behavior under shear stress using phase-field-crystal[J].Computational Materials Science,2016,121:143-150

    [15]GOMEZ H,NOGUEIRA X.An unconditionally energy-stable method for the phase field crystal equation[J].Computer Methods in Applied Mechanics and Engineering,2012,249-252:52-61.

    [16]ELDER K R,ROSSI G,KANERVA P,et al.Patterning of heteroepitaxial overlayers from nano to micron scales[J].Physical Review Letters,2012,108(22):226102.

    [17]盧成健,蔣麗婷,王玉玲,等.晶體相場(chǎng)法模擬小角度晶界的位錯(cuò)結(jié)構(gòu)及其演化[J].廣西科學(xué),2013,20(4):316-320. LU C J,JIANG L T,WANG Y L,et al.Dislocation structure evolution in low angle grain boundary[J].Guangxi Sciences,2013,20(4):316-320.

    [18]毛鴻,羅志榮,黃世葉,等.材料裂紋擴(kuò)展分叉機(jī)理的晶體相場(chǎng)法研究[J].廣西科學(xué),2015(5):499-505. MAO H,LUO Z R,HUANG S Y,et al.Phase-field-crystal modeling for crack propagation and branch of materials[J].Guangxi Sciences,2015(5):499-505.

    [19]張俊善.材料強(qiáng)度學(xué)[M].黑龍江:哈爾濱工業(yè)大學(xué)出版社,2004. ZHANG S J.Strength of Materials[M].Heilongjiang:Harbin Institute of Technology Press,2004.

    [20]師俊平,胡義鋒,羅峰.裂紋起裂方向及擴(kuò)展路徑的數(shù)值模擬[J].力學(xué)季刊,2013,34(02):233-239. SHI J P,HU Y F,LUO F.Numerical simulation of crack initiation directions and its propagation path[J].Chinese Quarterly of Mechanics,2013,34(2):233-239.

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

    Phase-field-crystal Simulation of Critical Conditions for the Crack Initiation of Nanoscale Notch

    WEN Zhenchuan,NING Suiting,ZHANG Rui,GAO Yingjun

    Key words:phase-field-crystal model,critical notch size,notch initiation

    Abstract:【Objective】To study the effect of initial crack size on crack initiation.【Methods】The dynamic of process crack propagation under the uni-axial tension deformation in single crystals is simulated by the phase-field-crystal (PFC) model.The free energy is used to analyze the influence of different initial crack size on the crack initiation in the process of crack initiation and crack propagation.【Results】When the initial crack of y direction size is smaller,there is no initiation of crack under the uniaxial tension in the x axis,although the notch size in x direction changes.when the y direction size reaches a certain value,the initial cracks will expand.The time of crack initiation is delayed with the increase of x direction size.When the initial crack of x direction size changes less,there is no initiation of crack,and the initial cracks will expand when the x direction size reaches a certain value.【Conclusion】The initiation of crack is closely related to the crack size of vertical stress direction.The crack extension occurs,if only the size of notch reaches a certain value.

    收稿日期:2016-07-25

    作者簡(jiǎn)介:溫振川(1991-),男,碩士研究生,主要從事納米材料設(shè)計(jì)與模擬實(shí)驗(yàn)研究。 **通信作者:高英俊(1962-),男,教授,博士生導(dǎo)師,主要從事材料納微結(jié)構(gòu)的設(shè)計(jì)與模擬實(shí)驗(yàn)研究,E-mail:gaoyj@gxu.edu.cn。

    中圖分類(lèi)號(hào):TG111

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

    文章編號(hào):1005-9164(2016)05-0448-06

    *國(guó)家自然科學(xué)基金項(xiàng)目(51161003,50661001),廣西自然科學(xué)基金重點(diǎn)項(xiàng)目(2012GXNSFDA053001)和廣西大學(xué)大創(chuàng)項(xiàng)目(201610593220,201610593218)資助。

    廣西科學(xué)Guangxi Sciences 2016,23(5):448~453

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

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

    猜你喜歡
    相場(chǎng)缺口步長(zhǎng)
    基于Armijo搜索步長(zhǎng)的BFGS與DFP擬牛頓法的比較研究
    基于子單元光滑有限元的混凝土相場(chǎng)損傷模型研究
    必須堵上尾款欠薪“缺口”
    堵缺口
    鑄件凝固微觀組織仿真程序開(kāi)發(fā)
    基于相場(chǎng)理論的瀝青自愈合微觀進(jìn)程與機(jī)理研究進(jìn)展
    石油瀝青(2018年1期)2018-04-12 07:31:51
    基于COMSOL的相場(chǎng)模擬研究
    科技視界(2017年8期)2017-07-31 10:31:17
    我國(guó)醫(yī)學(xué)物理師缺口巨大
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    一種新型光伏系統(tǒng)MPPT變步長(zhǎng)滯環(huán)比較P&O法
    精品日产1卡2卡| 国产日本99.免费观看| 国产一区二区激情短视频| 搡老岳熟女国产| 国产午夜福利久久久久久| 亚洲国产日韩欧美精品在线观看 | 日日爽夜夜爽网站| 又爽又黄无遮挡网站| 毛片女人毛片| 久久精品综合一区二区三区| 亚洲成人精品中文字幕电影| 我的老师免费观看完整版| 国产成+人综合+亚洲专区| 91大片在线观看| 久久精品成人免费网站| 国产91精品成人一区二区三区| av天堂在线播放| 欧美一区二区精品小视频在线| 亚洲国产日韩欧美精品在线观看 | 欧美 亚洲 国产 日韩一| 正在播放国产对白刺激| 狂野欧美白嫩少妇大欣赏| 久久久久九九精品影院| 亚洲七黄色美女视频| 亚洲国产欧美人成| 国产精品香港三级国产av潘金莲| 久久久国产欧美日韩av| 日韩精品青青久久久久久| 不卡一级毛片| 久久香蕉精品热| 成人国产综合亚洲| 国产亚洲精品综合一区在线观看 | 久久99热这里只有精品18| 婷婷六月久久综合丁香| 精华霜和精华液先用哪个| 亚洲欧美日韩无卡精品| 国产爱豆传媒在线观看 | 又大又爽又粗| 日韩av在线大香蕉| 国产伦人伦偷精品视频| 九色成人免费人妻av| 国产成年人精品一区二区| 免费电影在线观看免费观看| 99热只有精品国产| 波多野结衣高清作品| 久久这里只有精品中国| 亚洲,欧美精品.| 久久天堂一区二区三区四区| 久久精品国产亚洲av香蕉五月| 国产人伦9x9x在线观看| 精品一区二区三区四区五区乱码| 亚洲一卡2卡3卡4卡5卡精品中文| 青草久久国产| 给我免费播放毛片高清在线观看| 久久久久久久精品吃奶| 欧美黑人巨大hd| av超薄肉色丝袜交足视频| 最近最新免费中文字幕在线| 18美女黄网站色大片免费观看| 中国美女看黄片| 禁无遮挡网站| 午夜福利在线观看吧| 免费在线观看完整版高清| 五月伊人婷婷丁香| 成人国产一区最新在线观看| 亚洲精品在线观看二区| 亚洲成人中文字幕在线播放| 国产主播在线观看一区二区| 亚洲精品久久成人aⅴ小说| 久久这里只有精品中国| 久久中文看片网| www.999成人在线观看| 亚洲精品色激情综合| 国产欧美日韩精品亚洲av| 麻豆国产97在线/欧美 | xxxwww97欧美| 成人18禁在线播放| 亚洲国产中文字幕在线视频| 精品午夜福利视频在线观看一区| 成熟少妇高潮喷水视频| 亚洲国产中文字幕在线视频| 亚洲人成网站高清观看| 中文字幕久久专区| 窝窝影院91人妻| 国产亚洲精品av在线| www日本在线高清视频| 黄色成人免费大全| 亚洲熟女毛片儿| 亚洲精品中文字幕在线视频| 大型黄色视频在线免费观看| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩无卡精品| 久久天堂一区二区三区四区| 99精品欧美一区二区三区四区| 国产黄a三级三级三级人| 色综合站精品国产| 国产精品av久久久久免费| 国产成人精品久久二区二区91| 母亲3免费完整高清在线观看| 亚洲精品久久国产高清桃花| 久久午夜亚洲精品久久| 亚洲人成伊人成综合网2020| 亚洲男人的天堂狠狠| 日本在线视频免费播放| 脱女人内裤的视频| 精品高清国产在线一区| 757午夜福利合集在线观看| 激情在线观看视频在线高清| 亚洲熟女毛片儿| 999精品在线视频| 啦啦啦免费观看视频1| 美女黄网站色视频| 变态另类丝袜制服| 成人国产一区最新在线观看| 97碰自拍视频| 亚洲中文字幕一区二区三区有码在线看 | 男人的好看免费观看在线视频 | 欧美三级亚洲精品| 欧美日本亚洲视频在线播放| 午夜a级毛片| 免费在线观看完整版高清| 久久久水蜜桃国产精品网| 久久久久久久午夜电影| 国产伦人伦偷精品视频| 久久精品91无色码中文字幕| 两性夫妻黄色片| 亚洲人成伊人成综合网2020| 精品久久久久久成人av| 欧美大码av| 黑人欧美特级aaaaaa片| 日日爽夜夜爽网站| 99国产精品一区二区三区| 日本成人三级电影网站| 又粗又爽又猛毛片免费看| 中文字幕人成人乱码亚洲影| 国产精品一及| 国产精品九九99| 亚洲av熟女| 99精品久久久久人妻精品| 可以免费在线观看a视频的电影网站| 中出人妻视频一区二区| 色尼玛亚洲综合影院| 精品第一国产精品| 看片在线看免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产99白浆流出| av中文乱码字幕在线| 巨乳人妻的诱惑在线观看| 日韩欧美一区二区三区在线观看| 丁香欧美五月| 亚洲男人天堂网一区| 国产亚洲精品久久久久5区| 男女视频在线观看网站免费 | 日本熟妇午夜| 亚洲中文字幕日韩| 此物有八面人人有两片| 啦啦啦免费观看视频1| 九色成人免费人妻av| 免费看十八禁软件| 很黄的视频免费| 波多野结衣巨乳人妻| 我的老师免费观看完整版| 搡老岳熟女国产| 亚洲av日韩精品久久久久久密| 国产成人精品无人区| 国产成人一区二区三区免费视频网站| 国产精品 欧美亚洲| 啦啦啦观看免费观看视频高清| 久久久久国产精品人妻aⅴ院| 老汉色∧v一级毛片| 99re在线观看精品视频| 亚洲人成77777在线视频| a在线观看视频网站| 久久久久久亚洲精品国产蜜桃av| 亚洲天堂国产精品一区在线| 在线国产一区二区在线| 女同久久另类99精品国产91| 国产精品影院久久| 亚洲熟妇中文字幕五十中出| 露出奶头的视频| 又粗又爽又猛毛片免费看| 又粗又爽又猛毛片免费看| 精品久久久久久,| 九色成人免费人妻av| 国产伦人伦偷精品视频| 国产欧美日韩一区二区精品| 成人18禁在线播放| 9191精品国产免费久久| 亚洲专区国产一区二区| 欧美又色又爽又黄视频| 美女高潮喷水抽搐中文字幕| 性色av乱码一区二区三区2| 欧美日韩一级在线毛片| 成熟少妇高潮喷水视频| 亚洲 欧美一区二区三区| 日韩国内少妇激情av| ponron亚洲| 久久中文看片网| 国产成人精品无人区| 国产成人aa在线观看| 伦理电影免费视频| 国产伦一二天堂av在线观看| 麻豆久久精品国产亚洲av| 在线观看一区二区三区| 欧美性猛交黑人性爽| 1024视频免费在线观看| 女人被狂操c到高潮| 久久久久久久久免费视频了| 欧美乱色亚洲激情| 日本免费一区二区三区高清不卡| 色综合亚洲欧美另类图片| АⅤ资源中文在线天堂| 中亚洲国语对白在线视频| 欧美3d第一页| 99精品在免费线老司机午夜| 国产伦一二天堂av在线观看| 欧美日韩黄片免| 国产在线观看jvid| 亚洲无线在线观看| 给我免费播放毛片高清在线观看| 欧美成人免费av一区二区三区| 精品久久久久久久久久免费视频| 成人国产综合亚洲| 亚洲精品一区av在线观看| 久久精品国产清高在天天线| 18禁美女被吸乳视频| 美女 人体艺术 gogo| 精品国内亚洲2022精品成人| 久热爱精品视频在线9| 999精品在线视频| 精华霜和精华液先用哪个| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 最新在线观看一区二区三区| 淫妇啪啪啪对白视频| 超碰成人久久| 亚洲专区中文字幕在线| 99热只有精品国产| 69av精品久久久久久| 桃色一区二区三区在线观看| 国产av麻豆久久久久久久| 午夜两性在线视频| 国产三级在线视频| 天堂动漫精品| 18美女黄网站色大片免费观看| 成年女人毛片免费观看观看9| 国产精品电影一区二区三区| 黄色成人免费大全| 国产在线观看jvid| 日韩欧美免费精品| 精品一区二区三区四区五区乱码| 久久久久久久久免费视频了| 变态另类成人亚洲欧美熟女| 美女大奶头视频| 国产午夜精品久久久久久| tocl精华| 在线播放国产精品三级| 亚洲最大成人中文| 欧美一区二区精品小视频在线| 亚洲激情在线av| 久久久久久久精品吃奶| 桃色一区二区三区在线观看| 国产精品野战在线观看| 亚洲成人久久性| 十八禁网站免费在线| 国产成人精品无人区| 成人国语在线视频| netflix在线观看网站| 欧美日韩黄片免| 国产又色又爽无遮挡免费看| 中文字幕熟女人妻在线| 久久这里只有精品19| 色综合婷婷激情| 12—13女人毛片做爰片一| 欧美午夜高清在线| 最近最新免费中文字幕在线| 欧美中文日本在线观看视频| 久久精品国产综合久久久| 久久精品人妻少妇| 亚洲av成人av| 亚洲欧美日韩高清专用| 国产高清videossex| 激情在线观看视频在线高清| 美女午夜性视频免费| 午夜亚洲福利在线播放| 老司机午夜十八禁免费视频| 国产精品野战在线观看| 国产精品综合久久久久久久免费| 亚洲18禁久久av| 久久精品成人免费网站| 午夜成年电影在线免费观看| 国产精品一区二区三区四区免费观看 | 特级一级黄色大片| 国产av在哪里看| 亚洲av日韩精品久久久久久密| 国产熟女午夜一区二区三区| 精品一区二区三区av网在线观看| 国产亚洲精品久久久久久毛片| 黄色女人牲交| 久久久水蜜桃国产精品网| 99热只有精品国产| 一级片免费观看大全| 久久草成人影院| 激情在线观看视频在线高清| 亚洲美女视频黄频| 久久九九热精品免费| 国产成+人综合+亚洲专区| 男人舔奶头视频| 十八禁网站免费在线| 给我免费播放毛片高清在线观看| 亚洲天堂国产精品一区在线| 老司机在亚洲福利影院| 亚洲国产看品久久| 亚洲性夜色夜夜综合| 日本一区二区免费在线视频| 黄色丝袜av网址大全| 免费无遮挡裸体视频| av天堂在线播放| 国产黄色小视频在线观看| 国产视频内射| 在线观看舔阴道视频| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 久久久久国产一级毛片高清牌| 两性午夜刺激爽爽歪歪视频在线观看 | 无遮挡黄片免费观看| 黄频高清免费视频| 久久精品国产99精品国产亚洲性色| 亚洲自偷自拍图片 自拍| 午夜福利18| а√天堂www在线а√下载| 午夜激情福利司机影院| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 波多野结衣高清作品| 欧美一区二区国产精品久久精品 | 窝窝影院91人妻| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器 | 亚洲电影在线观看av| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 婷婷丁香在线五月| 国产黄片美女视频| 香蕉av资源在线| tocl精华| 国产精品野战在线观看| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 国产区一区二久久| 无限看片的www在线观看| 亚洲国产看品久久| 亚洲aⅴ乱码一区二区在线播放 | 人人妻,人人澡人人爽秒播| 在线观看美女被高潮喷水网站 | 国产精品 国内视频| 午夜福利成人在线免费观看| 婷婷精品国产亚洲av| 1024香蕉在线观看| 成人av在线播放网站| 亚洲中文日韩欧美视频| 亚洲全国av大片| 国产精品自产拍在线观看55亚洲| 中文字幕精品亚洲无线码一区| 日日摸夜夜添夜夜添小说| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 熟女少妇亚洲综合色aaa.| av福利片在线观看| 91麻豆精品激情在线观看国产| 给我免费播放毛片高清在线观看| 人妻丰满熟妇av一区二区三区| 99久久久亚洲精品蜜臀av| av有码第一页| 亚洲免费av在线视频| 久久久久国产精品人妻aⅴ院| 日韩欧美一区二区三区在线观看| 禁无遮挡网站| 99国产精品99久久久久| 美女大奶头视频| 亚洲精华国产精华精| 天堂影院成人在线观看| 悠悠久久av| 亚洲av电影在线进入| 成人三级做爰电影| 欧美不卡视频在线免费观看 | 国产精品亚洲av一区麻豆| av免费在线观看网站| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 亚洲欧美日韩东京热| 黄片小视频在线播放| 久99久视频精品免费| 久久精品亚洲精品国产色婷小说| avwww免费| 国产熟女午夜一区二区三区| 韩国av一区二区三区四区| avwww免费| 制服诱惑二区| 国产伦一二天堂av在线观看| 国产精品av视频在线免费观看| 1024手机看黄色片| 久久精品亚洲精品国产色婷小说| 国产主播在线观看一区二区| 青草久久国产| 久久久久免费精品人妻一区二区| 久久亚洲精品不卡| 久久久久九九精品影院| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 中亚洲国语对白在线视频| 免费在线观看日本一区| 2021天堂中文幕一二区在线观| 最近最新中文字幕大全免费视频| 丰满人妻熟妇乱又伦精品不卡| 国产爱豆传媒在线观看 | 真人做人爱边吃奶动态| 波多野结衣高清无吗| 69av精品久久久久久| 国产成人av激情在线播放| 欧美又色又爽又黄视频| 99久久精品国产亚洲精品| 日韩欧美 国产精品| 亚洲黑人精品在线| 丁香六月欧美| 人妻夜夜爽99麻豆av| 亚洲av第一区精品v没综合| 亚洲第一电影网av| 不卡一级毛片| 制服丝袜大香蕉在线| 精品久久久久久成人av| 一区二区三区高清视频在线| 51午夜福利影视在线观看| 亚洲av成人精品一区久久| 国产精品久久久久久亚洲av鲁大| 母亲3免费完整高清在线观看| 国产av麻豆久久久久久久| 日韩精品免费视频一区二区三区| 99在线视频只有这里精品首页| 成人永久免费在线观看视频| 久久中文看片网| 欧美3d第一页| 黄色 视频免费看| 国产成人精品无人区| 99国产精品99久久久久| 久久婷婷人人爽人人干人人爱| 国产成人一区二区三区免费视频网站| 午夜激情av网站| 欧美日韩中文字幕国产精品一区二区三区| 国产激情久久老熟女| 久久精品人妻少妇| 久久天躁狠狠躁夜夜2o2o| 欧美精品啪啪一区二区三区| 日韩欧美国产一区二区入口| 在线观看午夜福利视频| 无限看片的www在线观看| 在线永久观看黄色视频| 国产69精品久久久久777片 | 激情在线观看视频在线高清| 国产三级黄色录像| 最新美女视频免费是黄的| 亚洲成人精品中文字幕电影| www.www免费av| 欧美日本视频| 久久久久久久久久黄片| 欧美色欧美亚洲另类二区| 久久久久久免费高清国产稀缺| 叶爱在线成人免费视频播放| 国产人伦9x9x在线观看| 午夜免费成人在线视频| 又大又爽又粗| 欧美乱色亚洲激情| 十八禁网站免费在线| 美女扒开内裤让男人捅视频| 天天躁夜夜躁狠狠躁躁| 又大又爽又粗| 天天一区二区日本电影三级| 亚洲精品粉嫩美女一区| 蜜桃久久精品国产亚洲av| 97碰自拍视频| 高清毛片免费观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲七黄色美女视频| 丁香六月欧美| 亚洲激情在线av| 国产亚洲欧美98| 亚洲成人免费电影在线观看| 国产亚洲精品久久久久5区| 观看免费一级毛片| www日本黄色视频网| 非洲黑人性xxxx精品又粗又长| 小说图片视频综合网站| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 亚洲 国产 在线| 午夜免费激情av| 岛国在线观看网站| 午夜福利免费观看在线| 国产在线精品亚洲第一网站| 中文字幕人成人乱码亚洲影| x7x7x7水蜜桃| 国内少妇人妻偷人精品xxx网站 | 日韩av在线大香蕉| 久久久久免费精品人妻一区二区| www.999成人在线观看| 国产高清videossex| 美女免费视频网站| 亚洲欧美激情综合另类| 欧美日韩黄片免| 午夜福利在线观看吧| 日韩精品中文字幕看吧| 制服丝袜大香蕉在线| 国内精品久久久久精免费| 久久精品国产清高在天天线| 校园春色视频在线观看| 亚洲国产精品sss在线观看| 99久久国产精品久久久| 久热爱精品视频在线9| 成人午夜高清在线视频| 9191精品国产免费久久| www.999成人在线观看| 国内久久婷婷六月综合欲色啪| 无遮挡黄片免费观看| 操出白浆在线播放| 不卡一级毛片| 十八禁人妻一区二区| 白带黄色成豆腐渣| 一个人免费在线观看电影 | 色尼玛亚洲综合影院| 欧美在线黄色| 性欧美人与动物交配| 人人妻,人人澡人人爽秒播| 亚洲成人久久爱视频| 一区二区三区激情视频| 丝袜人妻中文字幕| 午夜激情av网站| 亚洲国产欧美一区二区综合| 99国产精品一区二区三区| 九色成人免费人妻av| 日韩大尺度精品在线看网址| 88av欧美| 久久久久精品国产欧美久久久| 99re在线观看精品视频| e午夜精品久久久久久久| 精品国产亚洲在线| 99riav亚洲国产免费| 一级黄色大片毛片| 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 国产精品自产拍在线观看55亚洲| 9191精品国产免费久久| 天天添夜夜摸| 免费一级毛片在线播放高清视频| 日本精品一区二区三区蜜桃| 91大片在线观看| 97人妻精品一区二区三区麻豆| aaaaa片日本免费| 嫩草影视91久久| 欧美日韩福利视频一区二区| 欧美三级亚洲精品| 高清毛片免费观看视频网站| 麻豆成人av在线观看| 最近最新免费中文字幕在线| 亚洲欧美精品综合一区二区三区| 亚洲av第一区精品v没综合| 无限看片的www在线观看| 天堂av国产一区二区熟女人妻 | 黄色视频不卡| 男男h啪啪无遮挡| 国产精品日韩av在线免费观看| 精品国产乱子伦一区二区三区| 美女高潮喷水抽搐中文字幕| 91老司机精品| 亚洲精品色激情综合| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 亚洲 欧美一区二区三区| 九色成人免费人妻av| 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区三区四区久久| 高潮久久久久久久久久久不卡| 亚洲熟妇中文字幕五十中出| 一级毛片女人18水好多| 九色成人免费人妻av| 天堂av国产一区二区熟女人妻 | 给我免费播放毛片高清在线观看| 欧美三级亚洲精品| tocl精华| 精品电影一区二区在线| 非洲黑人性xxxx精品又粗又长| 免费看a级黄色片| 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| 欧美极品一区二区三区四区| 99re在线观看精品视频| 日韩欧美在线二视频| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久久久久久| 美女高潮喷水抽搐中文字幕| 啦啦啦韩国在线观看视频| 床上黄色一级片| 99久久99久久久精品蜜桃| 国产亚洲精品久久久久5区| 亚洲一区中文字幕在线| 这个男人来自地球电影免费观看| 亚洲成人免费电影在线观看| 亚洲欧美日韩高清专用| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 五月伊人婷婷丁香| 国产精品 欧美亚洲| 国产精品免费一区二区三区在线| 久久热在线av|