• <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法
    麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 亚洲精品自拍成人| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 中出人妻视频一区二区| 国产探花极品一区二区| 亚洲精品久久久久久婷婷小说 | 国产女主播在线喷水免费视频网站 | 亚洲欧美精品综合久久99| 成人午夜精彩视频在线观看| 亚洲一区二区三区色噜噜| 草草在线视频免费看| 日本免费a在线| 日韩欧美精品免费久久| 亚洲国产精品国产精品| 亚洲成人av在线免费| 久久精品久久久久久噜噜老黄 | 成人性生交大片免费视频hd| 韩国av在线不卡| 校园春色视频在线观看| 国产探花在线观看一区二区| 看免费成人av毛片| 成人毛片60女人毛片免费| 国产高清激情床上av| 我要搜黄色片| 在线免费十八禁| 一卡2卡三卡四卡精品乱码亚洲| 黄色欧美视频在线观看| 97超碰精品成人国产| 直男gayav资源| 欧美潮喷喷水| 可以在线观看毛片的网站| 亚洲欧美日韩东京热| 亚洲中文字幕一区二区三区有码在线看| 午夜福利在线在线| 色综合亚洲欧美另类图片| 丰满乱子伦码专区| av在线播放精品| 性色avwww在线观看| 69人妻影院| 亚洲高清免费不卡视频| videossex国产| 老司机影院成人| 青春草国产在线视频 | 日韩国内少妇激情av| 日本成人三级电影网站| 天天躁夜夜躁狠狠久久av| 在线观看一区二区三区| 爱豆传媒免费全集在线观看| 狂野欧美激情性xxxx在线观看| 亚洲高清免费不卡视频| 99热6这里只有精品| 级片在线观看| 看黄色毛片网站| 久久热精品热| 99九九线精品视频在线观看视频| 91精品国产九色| 免费看光身美女| 美女大奶头视频| 我要搜黄色片| 日本在线视频免费播放| 夜夜看夜夜爽夜夜摸| 91麻豆精品激情在线观看国产| 欧美xxxx黑人xx丫x性爽| 欧美高清性xxxxhd video| 国产毛片a区久久久久| 哪个播放器可以免费观看大片| 观看美女的网站| 干丝袜人妻中文字幕| 欧美日韩国产亚洲二区| 午夜a级毛片| 久久韩国三级中文字幕| 老师上课跳d突然被开到最大视频| 国产精品一区二区在线观看99 | 你懂的网址亚洲精品在线观看 | 男人和女人高潮做爰伦理| 国产色婷婷99| 国产精品久久久久久亚洲av鲁大| 搡老妇女老女人老熟妇| 欧美高清性xxxxhd video| 亚洲国产精品合色在线| 国内少妇人妻偷人精品xxx网站| 久久精品人妻少妇| 99久国产av精品国产电影| 中文字幕精品亚洲无线码一区| 久久久精品大字幕| 夫妻性生交免费视频一级片| 久久久午夜欧美精品| 日韩欧美三级三区| 欧美zozozo另类| 亚洲人与动物交配视频| 中文字幕人妻熟人妻熟丝袜美| 国产黄片美女视频| 色噜噜av男人的天堂激情| 午夜亚洲福利在线播放| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩东京热| av在线观看视频网站免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲成人av在线免费| 九九在线视频观看精品| 国产一区二区在线观看日韩| 最近视频中文字幕2019在线8| 国产精品美女特级片免费视频播放器| 黄色欧美视频在线观看| 亚洲第一区二区三区不卡| 精品一区二区免费观看| 成人午夜高清在线视频| 女人被狂操c到高潮| 中文字幕av在线有码专区| ponron亚洲| 久久精品影院6| 99riav亚洲国产免费| 欧美成人免费av一区二区三区| 国产成人一区二区在线| 国产精品久久久久久亚洲av鲁大| 国产一区二区激情短视频| 午夜老司机福利剧场| 国产69精品久久久久777片| 免费一级毛片在线播放高清视频| 插逼视频在线观看| 久久亚洲精品不卡| 青青草视频在线视频观看| 国产精品一区www在线观看| 尤物成人国产欧美一区二区三区| 国产精品一区二区三区四区免费观看| 亚洲精品久久久久久婷婷小说 | 村上凉子中文字幕在线| 国产不卡一卡二| 国产精品三级大全| 一边摸一边抽搐一进一小说| 麻豆国产av国片精品| 日韩一区二区三区影片| 国产精品电影一区二区三区| 老熟妇乱子伦视频在线观看| 国产三级中文精品| 99热网站在线观看| 午夜福利在线观看吧| 欧美色视频一区免费| 十八禁国产超污无遮挡网站| 成年版毛片免费区| 国产一区二区在线观看日韩| 激情 狠狠 欧美| 亚洲欧美中文字幕日韩二区| 国产精品麻豆人妻色哟哟久久 | 成人无遮挡网站| 97热精品久久久久久| 成年免费大片在线观看| 国产精品麻豆人妻色哟哟久久 | 久久韩国三级中文字幕| 国产精品嫩草影院av在线观看| 国产单亲对白刺激| 赤兔流量卡办理| 悠悠久久av| 日韩欧美一区二区三区在线观看| 舔av片在线| 国产精品一区www在线观看| 国产乱人偷精品视频| 欧美不卡视频在线免费观看| 国产真实伦视频高清在线观看| 亚洲国产欧美人成| 午夜精品在线福利| 日本三级黄在线观看| 亚洲国产精品合色在线| 国产精品无大码| 综合色丁香网| 久久草成人影院| kizo精华| 男女下面进入的视频免费午夜| 一边摸一边抽搐一进一小说| 久久久久免费精品人妻一区二区| 国产高清视频在线观看网站| 精品久久久久久成人av| 久久精品人妻少妇| 一本久久中文字幕| 国产免费男女视频| 九色成人免费人妻av| 久久精品久久久久久久性| h日本视频在线播放| 欧美精品国产亚洲| 国产单亲对白刺激| 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 国产伦在线观看视频一区| 精品一区二区免费观看| 三级毛片av免费| 一个人看的www免费观看视频| 国产一区亚洲一区在线观看| 亚洲av电影不卡..在线观看| 日韩一本色道免费dvd| 狠狠狠狠99中文字幕| 久久久精品大字幕| 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 国产伦理片在线播放av一区 | 欧美zozozo另类| 亚洲真实伦在线观看| 看非洲黑人一级黄片| 亚洲国产欧洲综合997久久,| 99久久无色码亚洲精品果冻| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 日本黄大片高清| 97超视频在线观看视频| 好男人视频免费观看在线| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 黄片无遮挡物在线观看| 我要看日韩黄色一级片| 国产人妻一区二区三区在| 国产欧美日韩精品一区二区| 悠悠久久av| 国产男人的电影天堂91| 麻豆成人午夜福利视频| 99久久精品一区二区三区| 日本免费一区二区三区高清不卡| 国产一区二区在线av高清观看| 桃色一区二区三区在线观看| 两个人的视频大全免费| 青春草亚洲视频在线观看| 亚洲国产精品成人久久小说 | 久久中文看片网| 97人妻精品一区二区三区麻豆| 久久午夜亚洲精品久久| 色综合色国产| 亚洲国产欧美人成| 国产在线精品亚洲第一网站| 18禁黄网站禁片免费观看直播| 久久午夜亚洲精品久久| 欧美日韩综合久久久久久| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 欧美精品一区二区大全| av视频在线观看入口| 国产av麻豆久久久久久久| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区在线观看99 | 国产极品天堂在线| 男人和女人高潮做爰伦理| 亚洲国产精品合色在线| 我要搜黄色片| 欧美日本视频| 亚洲精品粉嫩美女一区| 老师上课跳d突然被开到最大视频| 男人舔奶头视频| 校园人妻丝袜中文字幕| 天堂√8在线中文| 亚洲欧美精品专区久久| 一级毛片电影观看 | 亚洲精品久久久久久婷婷小说 | 亚洲三级黄色毛片| 久久久精品大字幕| 国产伦精品一区二区三区四那| 国产成人午夜福利电影在线观看| 午夜免费男女啪啪视频观看| 91狼人影院| av福利片在线观看| 国产精品久久久久久亚洲av鲁大| 色哟哟·www| 三级毛片av免费| 日日摸夜夜添夜夜添av毛片| 18禁黄网站禁片免费观看直播| 夜夜爽天天搞| 级片在线观看| 精品一区二区三区人妻视频| 国产精品日韩av在线免费观看| 午夜激情福利司机影院| kizo精华| 26uuu在线亚洲综合色| 精品久久久久久久末码| 男女那种视频在线观看| 又粗又爽又猛毛片免费看| 青青草视频在线视频观看| 一卡2卡三卡四卡精品乱码亚洲| 美女脱内裤让男人舔精品视频 | 麻豆国产av国片精品| 欧美激情久久久久久爽电影| 一夜夜www| 天堂av国产一区二区熟女人妻| 久久这里只有精品中国| 99热这里只有精品一区| 久久久午夜欧美精品| 久久久久网色| 国产极品精品免费视频能看的| 色哟哟·www| 麻豆精品久久久久久蜜桃| 美女被艹到高潮喷水动态| 日韩一本色道免费dvd| 国产精品久久久久久亚洲av鲁大| 久久久午夜欧美精品| 欧美日韩精品成人综合77777| 亚洲av男天堂| videossex国产| 非洲黑人性xxxx精品又粗又长| 国产亚洲欧美98| 久久精品久久久久久久性| av在线亚洲专区| 精品一区二区三区人妻视频| АⅤ资源中文在线天堂| 少妇熟女欧美另类| 国产精品国产高清国产av| 国产成人91sexporn| 美女cb高潮喷水在线观看| 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| 夜夜爽天天搞| 欧美最黄视频在线播放免费| 亚洲真实伦在线观看| av卡一久久| 日韩视频在线欧美| 自拍偷自拍亚洲精品老妇| 女人十人毛片免费观看3o分钟| 久久久精品大字幕| 国产午夜精品久久久久久一区二区三区| 国产成人a∨麻豆精品| 欧美潮喷喷水| 精品久久国产蜜桃| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 久久精品国产99精品国产亚洲性色| 男女啪啪激烈高潮av片| 99在线视频只有这里精品首页| 小说图片视频综合网站| 好男人视频免费观看在线| 成人永久免费在线观看视频| 国产精品免费一区二区三区在线| 狂野欧美白嫩少妇大欣赏| 国产午夜精品一二区理论片| 神马国产精品三级电影在线观看| 哪个播放器可以免费观看大片| 色哟哟·www| 日本三级黄在线观看| 不卡一级毛片| 欧美日韩乱码在线| 午夜精品一区二区三区免费看| 亚洲国产高清在线一区二区三| 久久99精品国语久久久| 神马国产精品三级电影在线观看| 欧美bdsm另类| 亚洲av第一区精品v没综合| 亚洲va在线va天堂va国产| 男人的好看免费观看在线视频| 一本精品99久久精品77| 伦理电影大哥的女人| 国产精品国产三级国产av玫瑰| 99热6这里只有精品| 久久久久久久久久成人| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久精品电影小说 | 国产三级在线视频| 一级毛片我不卡| 午夜免费男女啪啪视频观看| 美女大奶头视频| 毛片女人毛片| 婷婷亚洲欧美| av又黄又爽大尺度在线免费看 | 午夜精品国产一区二区电影 | 日本黄色片子视频| 亚洲乱码一区二区免费版| 午夜激情欧美在线| 国产亚洲5aaaaa淫片| 美女 人体艺术 gogo| 又粗又爽又猛毛片免费看| 欧美xxxx性猛交bbbb| 男的添女的下面高潮视频| 国产精品一区www在线观看| ponron亚洲| 男插女下体视频免费在线播放| 欧美高清成人免费视频www| av天堂在线播放| av视频在线观看入口| 99热这里只有精品一区| 日韩大尺度精品在线看网址| 欧美高清成人免费视频www| 成人毛片60女人毛片免费| 99久久精品热视频| 国产一区二区三区在线臀色熟女| 18禁在线无遮挡免费观看视频| 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 亚洲精品成人久久久久久| 高清毛片免费观看视频网站| 一本精品99久久精品77| 日本与韩国留学比较| 高清毛片免费观看视频网站| 国产伦一二天堂av在线观看| 日韩三级伦理在线观看| 女的被弄到高潮叫床怎么办| 亚洲成人久久性| 嫩草影院入口| 尾随美女入室| 精品午夜福利在线看| 中文字幕免费在线视频6| 久久精品人妻少妇| 美女cb高潮喷水在线观看| 免费看光身美女| 久久久精品94久久精品| 精品人妻一区二区三区麻豆| 18禁黄网站禁片免费观看直播| 亚洲成人久久爱视频| 精品熟女少妇av免费看| 成人性生交大片免费视频hd| 自拍偷自拍亚洲精品老妇| 久久精品久久久久久久性| av在线观看视频网站免费| 国产亚洲欧美98| 婷婷色综合大香蕉| 日韩欧美国产在线观看| 亚洲内射少妇av| 国产精品野战在线观看| 偷拍熟女少妇极品色| 久久鲁丝午夜福利片| 男女视频在线观看网站免费| 熟女电影av网| 熟女人妻精品中文字幕| 2022亚洲国产成人精品| 毛片女人毛片| 看免费成人av毛片| 男女边吃奶边做爰视频| 欧洲精品卡2卡3卡4卡5卡区| 天天躁夜夜躁狠狠久久av| 欧美xxxx黑人xx丫x性爽| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| 高清毛片免费观看视频网站| 国产免费男女视频| 99热只有精品国产| 亚洲成a人片在线一区二区| 色综合色国产| 欧美性感艳星| a级一级毛片免费在线观看| 我的女老师完整版在线观看| 可以在线观看的亚洲视频| 看十八女毛片水多多多| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 在线天堂最新版资源| 观看美女的网站| 国内精品一区二区在线观看| 能在线免费看毛片的网站| av在线蜜桃| 秋霞在线观看毛片| 亚洲成人中文字幕在线播放| 国产乱人偷精品视频| 长腿黑丝高跟| 国产亚洲精品av在线| 国产高清有码在线观看视频| 亚洲人成网站在线播放欧美日韩| 亚洲成人久久爱视频| 亚洲乱码一区二区免费版| 一卡2卡三卡四卡精品乱码亚洲| av免费在线看不卡| 国产av麻豆久久久久久久| 一级毛片久久久久久久久女| 三级男女做爰猛烈吃奶摸视频| 午夜精品在线福利| 国产三级在线视频| 99久久成人亚洲精品观看| 国产高清视频在线观看网站| 久久久久免费精品人妻一区二区| 美女内射精品一级片tv| 麻豆久久精品国产亚洲av| 人妻系列 视频| 26uuu在线亚洲综合色| 丝袜喷水一区| 91久久精品国产一区二区三区| 久久鲁丝午夜福利片| 一本久久精品| 亚洲激情五月婷婷啪啪| 国产91av在线免费观看| 五月玫瑰六月丁香| 国产伦在线观看视频一区| 好男人视频免费观看在线| 两个人的视频大全免费| 免费看日本二区| 精品久久久久久久久亚洲| 如何舔出高潮| 日韩精品青青久久久久久| 男人的好看免费观看在线视频| .国产精品久久| 一级毛片电影观看 | 亚洲国产欧洲综合997久久,| 美女 人体艺术 gogo| 免费看光身美女| 国产高清视频在线观看网站| 一级黄片播放器| 在线观看66精品国产| 亚洲七黄色美女视频| 97超碰精品成人国产| 国产av不卡久久| 亚洲国产精品久久男人天堂| 亚洲精品日韩av片在线观看| 一进一出抽搐动态| 日本免费一区二区三区高清不卡| 国产精品嫩草影院av在线观看| 亚洲中文字幕日韩| 国产精品久久视频播放| av福利片在线观看| 免费人成视频x8x8入口观看| 黄色配什么色好看| 久久综合国产亚洲精品| 国产亚洲精品av在线| 国产精品精品国产色婷婷| a级一级毛片免费在线观看| 麻豆国产97在线/欧美| 国产精品电影一区二区三区| 国产成人精品婷婷| 国产探花在线观看一区二区| 波多野结衣高清作品| 日韩一区二区视频免费看| 国产精品99久久久久久久久| 男人的好看免费观看在线视频| 亚洲精品久久久久久婷婷小说 | 欧美最黄视频在线播放免费| 91在线精品国自产拍蜜月| 亚洲中文字幕一区二区三区有码在线看| 亚洲成a人片在线一区二区| 久久婷婷人人爽人人干人人爱| a级一级毛片免费在线观看| 91久久精品国产一区二区成人| 边亲边吃奶的免费视频| 一卡2卡三卡四卡精品乱码亚洲| 少妇的逼好多水| 岛国毛片在线播放| 亚洲成人久久爱视频| 亚洲欧洲国产日韩| 中文字幕av在线有码专区| 亚洲五月天丁香| 欧美+亚洲+日韩+国产| 校园春色视频在线观看| 国产精品日韩av在线免费观看| 白带黄色成豆腐渣| 人妻少妇偷人精品九色| 国产不卡一卡二| 级片在线观看| 国产精品一区二区三区四区久久| 日本一本二区三区精品| 国产私拍福利视频在线观看| 在线播放国产精品三级| 精品一区二区免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲av中文av极速乱| 国产淫片久久久久久久久| 亚洲成人av在线免费| 男女视频在线观看网站免费| 精品少妇黑人巨大在线播放 | 久久精品影院6| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 最近最新中文字幕大全电影3| 国产精品久久久久久精品电影| 国产久久久一区二区三区| 最新中文字幕久久久久| 日本在线视频免费播放| 色综合亚洲欧美另类图片| 观看免费一级毛片| 哪里可以看免费的av片| 在线国产一区二区在线| 免费观看人在逋| 国产成人freesex在线| 国产成人a∨麻豆精品| 高清午夜精品一区二区三区 | 悠悠久久av| 国产中年淑女户外野战色| 精品人妻视频免费看| 国产91av在线免费观看| 国产精品一及| 啦啦啦观看免费观看视频高清| 综合色av麻豆| 搞女人的毛片| 美女内射精品一级片tv| 亚洲中文字幕一区二区三区有码在线看| 午夜爱爱视频在线播放| 久久99蜜桃精品久久| 99热这里只有精品一区| 免费电影在线观看免费观看| 亚洲性久久影院| 国产片特级美女逼逼视频| 最新中文字幕久久久久| 91精品一卡2卡3卡4卡| 男人的好看免费观看在线视频| 免费观看a级毛片全部| 精品国内亚洲2022精品成人| 美女内射精品一级片tv| 黄片无遮挡物在线观看| 一级毛片aaaaaa免费看小| 毛片一级片免费看久久久久| 人人妻人人澡欧美一区二区| 欧美一级a爱片免费观看看| 色播亚洲综合网| 三级国产精品欧美在线观看| 国产日本99.免费观看| 美女大奶头视频| 校园春色视频在线观看| 高清日韩中文字幕在线| 国产伦精品一区二区三区四那| 久久久精品欧美日韩精品| 亚洲性久久影院| 美女内射精品一级片tv| 在线天堂最新版资源| 老司机影院成人| 亚洲精品亚洲一区二区| 99热这里只有是精品在线观看| 91aial.com中文字幕在线观看| 别揉我奶头 嗯啊视频| 丰满人妻一区二区三区视频av| 国产久久久一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲图色成人| 欧美精品国产亚洲| kizo精华|