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

    Cu、Nb對(duì)Sanicro 25鋼結(jié)構(gòu)穩(wěn)定性及力學(xué)性能影響的第一性原理研究

    2021-08-16 07:27:04賈瑞瑞韓培德
    關(guān)鍵詞:延展性泊松比晶面

    賈瑞瑞,董 楠,王 劍,王 蘭,劉 思,賈 偉,韓培德

    (1.太原理工大學(xué)材料科學(xué)與工程學(xué)院,太原 030024;2.太原理工大學(xué)新材料界面科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室,太原 030024)

    1 引 言

    高效清潔的超超臨界發(fā)電技術(shù)是燃煤火力電站的重點(diǎn)發(fā)展方向.為了進(jìn)一步提高熱效率,降低煤耗,630-650℃等級(jí)超超臨界機(jī)組是火電建設(shè)的下一個(gè)重要目標(biāo)[1].然而蒸汽參數(shù)的提高會(huì)加快機(jī)組中過熱器和再熱器等核心部件老化的速率,縮短其使用壽命.因此發(fā)展超超臨界機(jī)組,保證其持續(xù)、穩(wěn)定、安全運(yùn)行的關(guān)鍵是研發(fā)成本合理、結(jié)構(gòu)穩(wěn)定性良好、高溫力學(xué)性能優(yōu)異的結(jié)構(gòu)材料[2].Sanicro 25鋼具有良好的抗氧化性能、抗蠕變性能、抗蒸汽腐蝕性能和組織穩(wěn)定性,是更高參數(shù)下電站鍋爐中制造過熱器和再熱器的理想材料,有望用于630-650℃等級(jí)超超臨界火電機(jī)組中[3].

    對(duì)比HR3C和Sanicro 25鋼的化學(xué)成分可以發(fā)現(xiàn)合金化元素Cu、Nb等在結(jié)構(gòu)穩(wěn)定性和力學(xué)性能方面起到重要作用.趙博等人[4]研究表明提高Sanicro 25鋼中的Nb含量會(huì)促進(jìn)Z相和MX相的析出,有利于提高Sanicro 25鋼的蠕變強(qiáng)度.Zhou等人[5]的研究結(jié)果表明次生NbCrN,富Cu相,M23C6和MX的沉淀強(qiáng)化極大地提高了Sanicro 25鋼的蠕變斷裂強(qiáng)度.這些實(shí)驗(yàn)研究均表明Cu、Nb可以提高Sanicro 25鋼的力學(xué)性能.

    目前主要通過實(shí)驗(yàn)方法從析出相角度研究Cu、Nb對(duì)Sanicro 25鋼高溫力學(xué)性能的影響,然而僅通過實(shí)驗(yàn)方法研究其強(qiáng)化機(jī)理是非常有限的,仍需要從原子尺度進(jìn)行更詳細(xì)的分析.本文利用第一性原理從原子層次上研究了Cu、Nb對(duì)Sanicro 25鋼結(jié)構(gòu)穩(wěn)定性和力學(xué)性能的影響.首先通過形成能,結(jié)合能和吉布斯自由能分析Cu、Nb對(duì)Sanicro 25鋼結(jié)構(gòu)穩(wěn)定性和熱力學(xué)穩(wěn)定性的影響.之后通過計(jì)算彈性模量,泊松比,B/G和派納力等探討Cu、Nb對(duì)Sanicro 25鋼力學(xué)性能的影響,并從電荷密度圖和態(tài)密度圖進(jìn)一步分析其微觀作用機(jī)理.

    2 計(jì)算模型與方法

    Fe在自然狀態(tài)下存在體心立方(α-Fe和δ-Fe)和面心立方(γ-Fe)兩種晶體結(jié)構(gòu).本文的研究體系是奧氏體耐熱鋼Sanicro 25,因此選擇γ-Fe作為初始晶體結(jié)構(gòu).根據(jù)Sanicro 25鋼的化學(xué)成分,考慮合金原子數(shù)和實(shí)際計(jì)算效率,選擇含32個(gè)原子的2×2×2 Fe超胞作為基體,分別用8個(gè)Cr原子和8個(gè)Ni原子置換16個(gè)Fe原子搭建了Fe16Cr8Ni8模型,其對(duì)應(yīng)化學(xué)成分為:22 wt.%Cr,25 wt.% Ni和53 wt.% Fe,與Sanicro 25鋼中Fe、Cr、Ni的質(zhì)量分?jǐn)?shù)吻合.隨后在Fe16Cr8Ni8模型基礎(chǔ)上通過置換Fe引入不同含量的Cu和Nb,得到Fe15Cr8Ni8Cu、Fe14Cr8Ni8Cu2、Fe15Cr8Ni8Nb和Fe14Cr8Ni8Nb2模型如圖1所示.

    圖1 計(jì)算模型示意圖:(a)Fe16 Cr8 Ni8,(b)Fe15 Cr8 Ni8 X(c)Fe14 Cr8 Ni8 X2(X=Cu,Nb)Fig.1 Schematic diagrams of models for calculations:(a)Fe16 Cr8 Ni8,(b)Fe15 Cr8 Ni8 X(c)Fe14 Cr8 Ni8 X2(X=Cu,Nb)

    本文的計(jì)算任務(wù)都是在Materials Studio 8.0軟件的CASTEP模塊中完成的.計(jì)算參數(shù)設(shè)置如下[6]:交換-相關(guān)能選用廣義梯度近似GGA下的PBE泛函描述,并采用超軟贗勢(shì)描述離子與價(jià)電子之間的相互作用,平面波截?cái)嗄苤翟O(shè)定為500 eV.在結(jié)構(gòu)優(yōu)化期間,采用4×4×7 Monkhorst-Pack k點(diǎn)網(wǎng)格進(jìn)行布里淵區(qū)采樣.總能計(jì)算時(shí)采用自洽迭代的方法,自洽場(chǎng)的收斂精度取1.0×10-5eV/atom,每個(gè)原子上的力低于0.3 eV/nm,公差偏移低于1×10-4nm,應(yīng)力偏差低于0.03 GPa.

    3 結(jié)果與討論

    3.1 結(jié)構(gòu)穩(wěn)定性分析

    對(duì)上述模型進(jìn)行結(jié)構(gòu)優(yōu)化后得到的晶格參數(shù)如表1所示.在Fe16Cr8Ni8體系中加入Cu、Nb后形成的體系仍為四方晶系,但都發(fā)生了不同程度的體積膨脹.由于Cu、Nb的原子半徑相差較大,含Cu體系的體積變化相對(duì)于含Nb體系的體積變化較小.此外,隨著體系中Cu和Nb含量的增加,體系的體積也隨之增加.

    表1 不同模型結(jié)構(gòu)優(yōu)化后的晶格參數(shù)Table 1 Calculated lattice constants of different models after optimization

    形成能是指不同原子從其單質(zhì)狀態(tài)下生成化合物所釋放或吸收的能量.在Fe16Cr8Ni8體系中加入Cu、Nb后形成新體系的難易程度可以通過形成能來判斷,其計(jì)算公式如下所示[7]:

    其中Etotal表示每個(gè)體系經(jīng)過結(jié)構(gòu)優(yōu)化后在平衡晶格常數(shù)下的晶胞總能量,Ni是晶胞中原子i的個(gè)數(shù)(i=Fe,Cr,Ni,Cu或Nb),Eiatom是原子i處于單質(zhì)狀態(tài)下的單原子能量,在計(jì)算單原子能量時(shí)需使用與結(jié)構(gòu)優(yōu)化相同的計(jì)算條件.當(dāng)體系的形成能為負(fù)值時(shí),說明該體系易于形成,且形成能的絕對(duì)值越大,表明體系越容易形成.相反地,當(dāng)體系的形成能為正值時(shí),表明該體系難以形成,在體系形成過程中需要吸收能量[8].計(jì)算得到不同體系的形成能如圖2(a)示.從圖2(a)中可以發(fā)現(xiàn)所有體系的形成能值均在0.125~0.175 eV/atom之間,這說明所有體系都相對(duì)容易形成.在Fe16Cr8Ni8體系中加入Cu后形成新體系的形成能隨著Cu含量的增加而增加,這表明新體系相對(duì)于Fe16Cr8Ni8體系不易形成.相反地,在Fe16Cr8Ni8體系中加入Nb后形成新體系的形成能相對(duì)于Fe16Cr8Ni8體系有所降低,這表明Fe14Cr8Ni8Nb2體系和Fe15Cr8Ni8Nb體系相對(duì)于Fe16Cr8Ni8體系容易形成.但是隨著Nb含量的增加,體系的形成能值也隨之增加,即Fe14Cr8Ni8Nb2體系相對(duì)于Fe15Cr8Ni8Nb體系不易形成.

    為了進(jìn)一步分析Cu、Nb在Fe16Cr8Ni8體系中的結(jié)構(gòu)穩(wěn)定性,還計(jì)算了各個(gè)體系的結(jié)合能.結(jié)合能是將自由原子結(jié)合為晶體所釋放的能量,也就是將晶體分解成單個(gè)原子時(shí)外界所做的功.結(jié)合能為負(fù)值時(shí)其絕對(duì)值越大,體系的結(jié)構(gòu)穩(wěn)定性就越好.結(jié)合能的計(jì)算公式如下[9]:

    吉布斯自由能通常用于評(píng)估材料在高溫下的結(jié)構(gòu)穩(wěn)定性.當(dāng)吉布斯自由能為負(fù)值時(shí)表明體系可以穩(wěn)定存在,且吉布斯自由能越小,其熱力學(xué)穩(wěn)定性越好[10].為了研究在Fe16Cr8Ni8體系中加入Cu和Nb后體系在高溫下的結(jié)構(gòu)穩(wěn)定性,我們計(jì)算了不同體系吉布斯自由能隨溫度變化的曲線如圖3所示.可以看出所有體系的吉布斯自由能均隨著溫度的升高而逐漸減小,這說明所有體系隨著溫度的升高而更加穩(wěn)定,且在630-650℃下所有體系的結(jié)構(gòu)穩(wěn)定性更好.630-650℃時(shí)在Fe16Cr8Ni8體系中加入Cu和Nb后體系的吉布斯自由能相較于Fe16Cr8Ni8體系明顯下降,這說明在體系中加入Cu和Nb會(huì)提高其熱力學(xué)穩(wěn)定性,Nb效果更明顯.圖2和圖3一致說明Cu和Nb可以在Fe16Cr8Ni8體系中穩(wěn)定存在.

    圖2 不同體系的形成能Ef和結(jié)合能Eb(a)Ef;(b)EbFig.2 Formation energies Ef(a)and binding energies Eb(b)of different systems

    圖3 不同模型中吉布斯自由能與溫度的關(guān)系Fig.3 The relationships between Gibbs free energy and temperature in different models

    3.2 彈性性質(zhì)計(jì)算

    除了結(jié)構(gòu)穩(wěn)定性和熱力學(xué)穩(wěn)定性之外,力學(xué)穩(wěn)定性也是判斷晶體結(jié)構(gòu)是否能夠穩(wěn)定存在的一個(gè)重要依據(jù).彈性常數(shù)反映了晶體材料應(yīng)對(duì)外部應(yīng)變的能力,并且與晶體材料的力學(xué)性能密切相關(guān)[11].本文研究的四方晶系具有6個(gè)獨(dú)立彈性常數(shù)(C11,C12,C13,C33,C44和C66),力學(xué)穩(wěn)定的條件是:C11>0,C33>0,C44>0,C66>0,C11–C12>0,C11+C33–2C12>0且2(C11+C12)+C33+4C13>0[12].不同體系計(jì)算的彈性常數(shù)如表2所示,通過計(jì)算發(fā)現(xiàn)含Cu和Nb的體系均滿足力學(xué)穩(wěn)定條件,這說明所有體系的力學(xué)性能都很穩(wěn)定,均能穩(wěn)定存在.

    表2 計(jì)算不同體系的彈性常數(shù)Table 2 Calculated elastic constants for different systems

    材料的力學(xué)性質(zhì)在宏觀上體現(xiàn)為材料的延展性、脆性、硬度、強(qiáng)度等方面.與之相對(duì)應(yīng),在微觀上,通過一系列的物理量如彈性常數(shù)、體積模量、剪切模量、楊氏模量、泊松比、體模量和剪切模量的比值等來體現(xiàn).

    多晶的彈性模量可由單晶的彈性常數(shù)進(jìn)一步計(jì)算得到,本文采用Voight-Reuss-Hill的計(jì)算方法,體積模量B,剪切模量G和楊氏模量E的計(jì)算公式如下[13]:

    不同體系的B,G和E如表3所示.B表征材料在施加的應(yīng)力下抵抗變形的能力,B越大,材料抗壓縮能力越強(qiáng).通常B越大材料的強(qiáng)度越高.G表征材料抵抗剪切應(yīng)變的能力,G越大,材料的抗剪切應(yīng)變能力越強(qiáng).E是描述固體材料抵抗形變能力的物理量,E越小材料的剛度就越小.從表3中可以發(fā)現(xiàn)隨著Fe16Cr8Ni8體系中Cu和Nb含量的增加,體系的B、G和E均逐漸減小,這說明隨著Cu和Nb含量的增加,體系的強(qiáng)度、抗壓縮能力、抗剪切應(yīng)變能力和剛度都會(huì)隨之下降.

    表3 計(jì)算不同體系的彈性模量Table 3 Calculated elastic moduli of different systems

    此外,泊松比是反映材料橫向變形的彈性常數(shù),通常用于估量材料晶體抗剪切應(yīng)變的穩(wěn)定性,泊松比的值越大,材料的塑性就越好.泊松比的計(jì)算公式如下所示[14]:

    不同體系的泊松比值如圖4所示.可以發(fā)現(xiàn)在Fe16Cr8Ni8體系中加入Cu后體系的泊松比降低,而隨著Cu含量的增加,體系的泊松比沒有明顯變化,這說明在Fe16Cr8Ni8體系中加入Cu會(huì)降低體系的塑性,而隨著Cu含量的增加,體系的塑性沒有明顯的變化.相反在Fe16Cr8Ni8體系中加入Nb后體系的泊松比值增加,且隨著Nb含量的增加,體系的泊松比逐漸增加,這說明在Fe16Cr8Ni8體系中加入Nb會(huì)提高體系的塑性,隨著Nb含量的增加,體系的塑性也會(huì)明顯增加.

    圖4 計(jì)算不同體系的泊松比值Fig.4 Calculated Poisson's ratio values of different systems

    B/G通常用于評(píng)估材料的延展性,當(dāng)B/G≥1.75時(shí),材料的延展性較好,B/G值越大,材料的延展性就越好[6].圖5為不同體系的B/G值,可以發(fā)現(xiàn)只有Fe14Cr8Ni8Nb2體系的B/G值大于1.75,這說明該體系的延展性較好.在Fe16Cr8Ni8體系中加入Cu后體系的B/G降低,說明在Fe16Cr8Ni8體系中加入Cu會(huì)降低體系的延展性.相反地,在Fe16Cr8Ni8體系中加入Nb后體系的B/G值隨著Nb含量的增加而增加,這說明Fe16Cr8Ni8體系中加入Nb原子會(huì)提高體系的延展性,且隨著Nb含量的增加,體系的延展性也會(huì)隨之提高.

    圖5 計(jì)算不同體系的B/G值Fig.5 Calculated B/G ratio values of different systems

    3.3 派納力計(jì)算

    派納力也可以判斷體系的塑性變形能力,從而進(jìn)一步研究其力學(xué)性能.在位錯(cuò)移動(dòng)的過程中(到達(dá)平衡位置之前),位錯(cuò)中心將偏離平衡位置使晶體能量增加構(gòu)成能壘,這就是由晶體點(diǎn)陣造成的位錯(cuò)運(yùn)動(dòng)阻力,也稱為派納力.派納力越小,位錯(cuò)越容易發(fā)生滑移,其計(jì)算公式如下[15]:

    其中p是滑動(dòng)面的面間距,q是滑動(dòng)方向上的原子間距,W=p/(1-ν)是位錯(cuò)寬度,ν是泊松比.本文選擇計(jì)算不同體系的(111)晶面和(001)晶面的派納力,如圖6所示.其中(111)晶面和(001)晶面的滑移方向分別為[11 2-]和[010].通過對(duì)比發(fā)現(xiàn),(111)晶面的派納力明顯小于(001)晶面的派納力,其計(jì)算結(jié)果符合面心立方結(jié)構(gòu)優(yōu)先沿其密排面(111)面滑移的規(guī)律.此外,在Fe16Cr8Ni8體系中加入Cu和Nb對(duì)(111)晶面的派納力均無明顯影響,Cu對(duì)(001)晶面的派納力也沒有明顯影響.而在Fe16Cr8Ni8體系中加入Nb會(huì)明顯降低(001)晶面的派納力,促進(jìn)位錯(cuò)沿[010]方向發(fā)生滑移.一般來說,材料中位錯(cuò)越容易發(fā)生滑移,其塑性變形能力就越好.泊松比和派納力的計(jì)算結(jié)果一致說明了Nb會(huì)提高體系的塑性變形能力.

    圖6 計(jì)算不同體系的派納力:(a)(111)(b)(001)晶面Fig.6 Calculated Peierls stresses of different systems:(a)(111)crystal plane(b)(001)crystal plane

    3.4 電子性質(zhì)分析

    圖7為不同體系的電荷密度圖,觀察Fe16Cr8Ni8體系發(fā)現(xiàn)原子周圍電荷分布均勻,近似球狀的電荷區(qū)域?qū)?yīng)體系中金屬鍵的存在.對(duì)比發(fā)現(xiàn)在Fe16Cr8Ni8體系中加入Cu后周圍電荷分布減少,Cu與周圍原子之間的結(jié)合略微變差,因此Cu會(huì)降低體系的結(jié)構(gòu)穩(wěn)定性.而在Fe16Cr8Ni8體系中加入Nb后周邊電荷分布增加,Nb與周圍原子之間的結(jié)合更加緊密,因此Nb提高了體系的結(jié)構(gòu)穩(wěn)定性.電荷密度分布結(jié)果與結(jié)構(gòu)穩(wěn)定性結(jié)果一致.

    圖7 不同體系沿(100)面的電荷密度圖:(a)Fe16 Cr8 Ni8;(b)Fe15 Cr8 Ni8 Cu;(c)Fe14 Cr8 Ni8 Cu2;(d)Fe15 Cr8 Ni8 Nb;(e)Fe1 4 Cr8 Ni8 Nb2Fig.7 Charge densities distribution maps of different systems in(100)plane:(a)Fe1 6 Cr8 Ni8;(b)Fe15 Cr8 Ni8 Cu;(c)Fe14 Cr8 Ni8 Cu2;(d)Fe15 Cr8 Ni8 Nb;(e)Fe1 4 Cr8 Ni8 Nb2

    繪制不同體系的總態(tài)密度圖如圖8所示.通過觀察發(fā)現(xiàn)所有體系的總態(tài)密度曲線的形狀類似,這說明Cu、Nb沒有改變體系的能級(jí)結(jié)構(gòu).所有體系在費(fèi)米能級(jí)處的電子數(shù)均不為零,因此所有體系均表現(xiàn)出金屬特性.此外,電子數(shù)密度越大,費(fèi)米能級(jí)處的電子數(shù)越高,該體系的電化學(xué)穩(wěn)定性就越差[16].在Fe16Cr8Ni8體系中加入Cu、Nb后體系費(fèi)米能級(jí)處的電子數(shù)均減少,這表明含Cu、Nb的體系的電化學(xué)穩(wěn)定性均有所提高,且Nb效果更為明顯.

    圖8 不同體系的總態(tài)密度曲線Fig.8 Total densities of states of different systems

    4 結(jié) 論

    本文在Fe-Cr-Ni體系的基礎(chǔ)上,利用第一性原理研究了Cu、Nb對(duì)Sanicro 25鋼結(jié)構(gòu)穩(wěn)定性和力學(xué)性能的影響.結(jié)果表明:在Fe-Cr-Ni體系中加入Cu、Nb后體系仍具有較好的結(jié)構(gòu)穩(wěn)定性,且在630-650℃下結(jié)構(gòu)穩(wěn)定性更好;Cu會(huì)降低體系的抗壓縮能力、抗剪切應(yīng)變能力、剛度、強(qiáng)度、塑性和延展性,而Nb會(huì)降低體系的抗壓縮能力、抗剪切應(yīng)變能力、剛度和強(qiáng)度,提高塑性和延展性.

    猜你喜歡
    延展性泊松比晶面
    乙酸乙酯與ε-CL-20不同晶面的微觀作用機(jī)制
    突發(fā)重大公共衛(wèi)生事件響應(yīng)中的出版偏倚、科學(xué)客觀性和科學(xué)倫理——“奧司他韋事件”及其延展性思考
    具有負(fù)泊松比效應(yīng)的紗線研發(fā)
    負(fù)泊松比功能的結(jié)構(gòu)復(fù)合紡紗技術(shù)進(jìn)展
    毛紡科技(2020年6期)2021-01-06 03:41:48
    考慮粘彈性泊松比的固體推進(jìn)劑蠕變型本構(gòu)模型①
    NaCl單晶非切割面晶面的X射線衍射
    固體推進(jìn)劑粘彈性泊松比應(yīng)變率-溫度等效關(guān)系
    (100)/(111)面金剛石膜抗氧等離子刻蝕能力
    不同硅晶面指數(shù)上的類倒金字塔結(jié)構(gòu)研究與分析?
    課堂有效教學(xué)之大學(xué)課堂
    只有这里有精品99| 亚洲精品日本国产第一区| 插阴视频在线观看视频| 欧美国产精品一级二级三级| 91国产中文字幕| 啪啪无遮挡十八禁网站| 国产在线视频一区二区| 99在线人妻在线中文字幕 | 女同久久另类99精品国产91| 一边摸一边做爽爽视频免费| 水蜜桃什么品种好| 欧美激情久久久久久爽电影 | 国产成人精品久久二区二区免费| 国产在视频线精品| 亚洲色图av天堂| 亚洲国产成人一精品久久久| 久久久久久人人人人人| h视频一区二区三区| 国产成人精品在线电影| 精品久久久久久电影网| 美国免费a级毛片| 久久久久久人人人人人| 性高湖久久久久久久久免费观看| 国产成人av教育| 脱女人内裤的视频| a级毛片黄视频| 日韩成人在线观看一区二区三区| 日本a在线网址| 免费不卡黄色视频| 又黄又粗又硬又大视频| 男女午夜视频在线观看| 久久国产精品大桥未久av| 日本一区二区免费在线视频| 啦啦啦视频在线资源免费观看| 丰满饥渴人妻一区二区三| 国产精品久久久人人做人人爽| 大型av网站在线播放| 中文字幕色久视频| 五月天丁香电影| 精品亚洲成国产av| 久久99一区二区三区| 国产区一区二久久| 国产精品一区二区免费欧美| 人人妻人人添人人爽欧美一区卜| 老熟妇乱子伦视频在线观看| 精品视频人人做人人爽| 欧美在线一区亚洲| 多毛熟女@视频| 亚洲自偷自拍图片 自拍| 十八禁高潮呻吟视频| 高清视频免费观看一区二区| 精品国产乱码久久久久久小说| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久 | 首页视频小说图片口味搜索| 国产精品国产高清国产av | 在线看a的网站| 日韩大片免费观看网站| 精品久久久精品久久久| 国产精品麻豆人妻色哟哟久久| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 免费在线观看视频国产中文字幕亚洲| 91成人精品电影| 99精国产麻豆久久婷婷| 777米奇影视久久| 国产精品熟女久久久久浪| tube8黄色片| 免费观看av网站的网址| √禁漫天堂资源中文www| 一边摸一边抽搐一进一出视频| 亚洲欧美色中文字幕在线| 色94色欧美一区二区| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 超碰97精品在线观看| 18禁国产床啪视频网站| 中文欧美无线码| 亚洲欧美日韩高清在线视频 | 一本—道久久a久久精品蜜桃钙片| 国产成人精品久久二区二区免费| 在线观看一区二区三区激情| 啪啪无遮挡十八禁网站| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 精品久久久久久电影网| 老熟妇乱子伦视频在线观看| 高清视频免费观看一区二区| 日本撒尿小便嘘嘘汇集6| 欧美乱码精品一区二区三区| 国产精品香港三级国产av潘金莲| 免费看十八禁软件| 日日夜夜操网爽| 久久久久久久国产电影| 51午夜福利影视在线观看| 91麻豆av在线| 亚洲精品在线美女| 国产在线免费精品| 大码成人一级视频| 亚洲五月色婷婷综合| 亚洲精品一二三| 久久精品亚洲av国产电影网| 脱女人内裤的视频| 国产一区有黄有色的免费视频| 高清黄色对白视频在线免费看| 欧美中文综合在线视频| 丝袜在线中文字幕| 国产av精品麻豆| 国产一区二区在线观看av| 怎么达到女性高潮| 久久久久久久大尺度免费视频| 丝袜在线中文字幕| 俄罗斯特黄特色一大片| 人人妻人人澡人人看| 国精品久久久久久国模美| 又紧又爽又黄一区二区| 国产精品国产av在线观看| 丝袜美腿诱惑在线| 欧美另类亚洲清纯唯美| 精品福利观看| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 少妇被粗大的猛进出69影院| 亚洲专区中文字幕在线| 高清欧美精品videossex| 狂野欧美激情性xxxx| 大码成人一级视频| 国产高清videossex| 免费av中文字幕在线| 亚洲情色 制服丝袜| 色播在线永久视频| 不卡av一区二区三区| 久久 成人 亚洲| 亚洲中文日韩欧美视频| 91精品国产国语对白视频| 9191精品国产免费久久| 在线十欧美十亚洲十日本专区| 在线观看免费视频日本深夜| 三上悠亚av全集在线观看| 久久久精品免费免费高清| 亚洲精品美女久久av网站| 一级片免费观看大全| 波多野结衣一区麻豆| 欧美亚洲 丝袜 人妻 在线| 欧美乱妇无乱码| 欧美日韩视频精品一区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产一区二区精华液| 淫妇啪啪啪对白视频| 久久久久网色| 一级黄色大片毛片| 国产精品亚洲一级av第二区| 老司机靠b影院| 在线观看一区二区三区激情| 久久精品aⅴ一区二区三区四区| 女人久久www免费人成看片| 人成视频在线观看免费观看| 男人操女人黄网站| 最新在线观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 中国美女看黄片| 欧美在线一区亚洲| 亚洲av日韩在线播放| 这个男人来自地球电影免费观看| 在线观看免费午夜福利视频| av天堂久久9| 欧美精品av麻豆av| 夫妻午夜视频| 免费在线观看视频国产中文字幕亚洲| 精品一区二区三区视频在线观看免费 | 亚洲国产中文字幕在线视频| 欧美大码av| 免费日韩欧美在线观看| 欧美亚洲日本最大视频资源| 免费黄频网站在线观看国产| 捣出白浆h1v1| 精品免费久久久久久久清纯 | svipshipincom国产片| 午夜激情av网站| 日韩成人在线观看一区二区三区| 久久久久久久久久久久大奶| 欧美日韩国产mv在线观看视频| 老熟女久久久| 亚洲精品中文字幕一二三四区 | 午夜视频精品福利| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 免费在线观看日本一区| 欧美精品亚洲一区二区| 天堂8中文在线网| 国产高清国产精品国产三级| 国产精品一区二区免费欧美| 王馨瑶露胸无遮挡在线观看| 国产视频一区二区在线看| 老熟妇仑乱视频hdxx| 天堂动漫精品| 一边摸一边抽搐一进一小说 | 黄频高清免费视频| 欧美精品一区二区免费开放| 黄色片一级片一级黄色片| 最黄视频免费看| 在线观看免费高清a一片| 黄色 视频免费看| 蜜桃在线观看..| 亚洲av日韩精品久久久久久密| 亚洲一区二区三区欧美精品| 国产精品免费大片| 亚洲成人手机| 老司机亚洲免费影院| 丁香六月欧美| 超碰成人久久| 黄网站色视频无遮挡免费观看| 精品一区二区三区视频在线观看免费 | 亚洲一区中文字幕在线| 男人舔女人的私密视频| av在线播放免费不卡| 19禁男女啪啪无遮挡网站| 18在线观看网站| 色婷婷av一区二区三区视频| av片东京热男人的天堂| 亚洲国产欧美在线一区| 美女福利国产在线| 亚洲色图av天堂| 啦啦啦免费观看视频1| 岛国毛片在线播放| 中文欧美无线码| 精品第一国产精品| 久久人人爽av亚洲精品天堂| 我的亚洲天堂| 老司机深夜福利视频在线观看| 纯流量卡能插随身wifi吗| 男人舔女人的私密视频| 1024香蕉在线观看| 欧美精品人与动牲交sv欧美| e午夜精品久久久久久久| 精品久久蜜臀av无| 亚洲精品国产色婷婷电影| 99国产精品99久久久久| 动漫黄色视频在线观看| 国产91精品成人一区二区三区 | 制服人妻中文乱码| 国产人伦9x9x在线观看| 欧美精品啪啪一区二区三区| 老熟妇仑乱视频hdxx| a级片在线免费高清观看视频| 高清黄色对白视频在线免费看| 亚洲第一欧美日韩一区二区三区 | 成人国产一区最新在线观看| 欧美老熟妇乱子伦牲交| 精品福利永久在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 丰满迷人的少妇在线观看| 91九色精品人成在线观看| 国产精品久久久久成人av| 色在线成人网| 久久午夜亚洲精品久久| 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美日韩高清在线视频 | www.精华液| 久久婷婷成人综合色麻豆| 一二三四社区在线视频社区8| 女性生殖器流出的白浆| 交换朋友夫妻互换小说| 久久久精品国产亚洲av高清涩受| 正在播放国产对白刺激| 午夜91福利影院| 一级黄色大片毛片| 亚洲精品国产区一区二| 国产高清激情床上av| 亚洲精品国产精品久久久不卡| 亚洲精品国产一区二区精华液| 亚洲欧美激情在线| 高清黄色对白视频在线免费看| 天堂8中文在线网| 国产野战对白在线观看| 人人妻人人澡人人看| 一区二区三区国产精品乱码| 国产一区有黄有色的免费视频| 99国产综合亚洲精品| 性少妇av在线| 97人妻天天添夜夜摸| 久久人妻福利社区极品人妻图片| 青青草视频在线视频观看| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 人人澡人人妻人| 亚洲精品国产区一区二| 久久久水蜜桃国产精品网| 国产国语露脸激情在线看| 国产亚洲精品一区二区www | av在线播放免费不卡| 两性夫妻黄色片| 91九色精品人成在线观看| 两个人免费观看高清视频| 两性夫妻黄色片| 黄片播放在线免费| 国产xxxxx性猛交| a级毛片在线看网站| 自拍欧美九色日韩亚洲蝌蚪91| 欧美 日韩 精品 国产| 久久人妻福利社区极品人妻图片| 国产精品电影一区二区三区 | 亚洲精品国产精品久久久不卡| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 亚洲性夜色夜夜综合| 久久人妻av系列| 人人妻,人人澡人人爽秒播| 欧美精品一区二区大全| 丝袜美足系列| 亚洲视频免费观看视频| 久久精品成人免费网站| 国产97色在线日韩免费| 美女扒开内裤让男人捅视频| 日本黄色视频三级网站网址 | 性色av乱码一区二区三区2| 99香蕉大伊视频| 精品午夜福利视频在线观看一区 | 国产在线视频一区二区| 成人国产av品久久久| 啦啦啦视频在线资源免费观看| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线 | 国产日韩一区二区三区精品不卡| 国产一区二区三区视频了| 狠狠婷婷综合久久久久久88av| 黄色丝袜av网址大全| 丁香六月天网| 99热国产这里只有精品6| 一级a爱视频在线免费观看| 国产成人精品在线电影| 国产在线精品亚洲第一网站| 侵犯人妻中文字幕一二三四区| 动漫黄色视频在线观看| 香蕉丝袜av| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 国产又爽黄色视频| 亚洲黑人精品在线| 精品少妇久久久久久888优播| 建设人人有责人人尽责人人享有的| 午夜福利欧美成人| 麻豆国产av国片精品| 久热爱精品视频在线9| 十八禁网站网址无遮挡| 国产一区二区 视频在线| 国产男靠女视频免费网站| 亚洲熟女精品中文字幕| 麻豆成人av在线观看| 午夜福利免费观看在线| 操美女的视频在线观看| 建设人人有责人人尽责人人享有的| 免费观看人在逋| 人人妻人人添人人爽欧美一区卜| 亚洲视频免费观看视频| 美女扒开内裤让男人捅视频| 国产在线观看jvid| 国产亚洲欧美精品永久| 国产精品影院久久| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 欧美精品一区二区免费开放| 女同久久另类99精品国产91| 亚洲一区二区三区欧美精品| 欧美日韩精品网址| 宅男免费午夜| 91老司机精品| 欧美激情 高清一区二区三区| 免费在线观看黄色视频的| 精品一区二区三区视频在线观看免费 | 日本一区二区免费在线视频| 亚洲一区中文字幕在线| 51午夜福利影视在线观看| 99精品在免费线老司机午夜| 日韩制服丝袜自拍偷拍| 美女主播在线视频| 每晚都被弄得嗷嗷叫到高潮| tocl精华| 欧美性长视频在线观看| 少妇裸体淫交视频免费看高清 | 欧美日韩黄片免| 久久亚洲真实| 变态另类成人亚洲欧美熟女 | 91精品国产国语对白视频| 久久精品aⅴ一区二区三区四区| 免费久久久久久久精品成人欧美视频| 国产真人三级小视频在线观看| 色婷婷av一区二区三区视频| 在线观看人妻少妇| 性少妇av在线| 成人国产一区最新在线观看| 亚洲成a人片在线一区二区| 自线自在国产av| 国产麻豆69| 老司机靠b影院| 丝袜人妻中文字幕| 亚洲精品中文字幕在线视频| 国产成人影院久久av| 天堂俺去俺来也www色官网| 久久久久久久久久久久大奶| 久久精品亚洲熟妇少妇任你| 叶爱在线成人免费视频播放| 一本久久精品| 亚洲九九香蕉| www日本在线高清视频| 大陆偷拍与自拍| 老司机午夜十八禁免费视频| 国产精品98久久久久久宅男小说| 在线永久观看黄色视频| 欧美成狂野欧美在线观看| 亚洲色图av天堂| 久久中文字幕一级| 不卡av一区二区三区| 中文字幕人妻丝袜制服| 一本大道久久a久久精品| 国产一区二区激情短视频| 正在播放国产对白刺激| 他把我摸到了高潮在线观看 | 国产精品一区二区精品视频观看| 国产欧美亚洲国产| av网站在线播放免费| 一个人免费看片子| 亚洲精华国产精华精| 青青草视频在线视频观看| 国产高清激情床上av| 成人国语在线视频| 97在线人人人人妻| 超色免费av| 五月天丁香电影| 国产精品免费一区二区三区在线 | 免费看十八禁软件| 美女国产高潮福利片在线看| 亚洲欧美一区二区三区黑人| 精品国内亚洲2022精品成人 | 午夜福利影视在线免费观看| 久久中文看片网| 夫妻午夜视频| 精品亚洲乱码少妇综合久久| 性少妇av在线| 亚洲av成人一区二区三| 日本黄色日本黄色录像| 国产成人一区二区三区免费视频网站| 色视频在线一区二区三区| 在线观看人妻少妇| bbb黄色大片| 久久久久久久大尺度免费视频| 最近最新中文字幕大全免费视频| 91老司机精品| 亚洲第一青青草原| 午夜成年电影在线免费观看| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 12—13女人毛片做爰片一| 丝袜人妻中文字幕| 亚洲精品国产区一区二| 99香蕉大伊视频| 自线自在国产av| 男女边摸边吃奶| 天天躁日日躁夜夜躁夜夜| 叶爱在线成人免费视频播放| 在线av久久热| 国产99久久九九免费精品| 在线观看人妻少妇| 亚洲免费av在线视频| 亚洲全国av大片| 丝袜在线中文字幕| 一级,二级,三级黄色视频| 日日摸夜夜添夜夜添小说| 老汉色av国产亚洲站长工具| 国产高清videossex| 无遮挡黄片免费观看| 亚洲国产欧美一区二区综合| 成人精品一区二区免费| 日韩欧美一区视频在线观看| 三级毛片av免费| 精品人妻1区二区| 精品一区二区三区四区五区乱码| 久久久久网色| 亚洲性夜色夜夜综合| 青青草视频在线视频观看| av线在线观看网站| 极品教师在线免费播放| 丁香六月欧美| 久久国产精品大桥未久av| 一区二区三区乱码不卡18| 亚洲第一av免费看| 国产成人影院久久av| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸| 精品国产国语对白av| 亚洲成人免费av在线播放| 午夜福利在线免费观看网站| av国产精品久久久久影院| 国产精品国产av在线观看| 电影成人av| 欧美中文综合在线视频| 女人被躁到高潮嗷嗷叫费观| 国产av又大| 午夜精品国产一区二区电影| 国产精品秋霞免费鲁丝片| 老司机在亚洲福利影院| 成年人免费黄色播放视频| 两个人看的免费小视频| 99re6热这里在线精品视频| 成人影院久久| 国产亚洲精品第一综合不卡| 制服人妻中文乱码| 国产真人三级小视频在线观看| 亚洲色图 男人天堂 中文字幕| 日韩一卡2卡3卡4卡2021年| 悠悠久久av| 欧美精品人与动牲交sv欧美| 国产免费视频播放在线视频| 国产精品香港三级国产av潘金莲| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 亚洲精品国产一区二区精华液| h视频一区二区三区| 国产一区二区在线观看av| 国产一区二区三区视频了| 久久久国产成人免费| 亚洲国产成人一精品久久久| 一区二区三区精品91| 老司机亚洲免费影院| 免费高清在线观看日韩| 日韩一卡2卡3卡4卡2021年| 日韩 欧美 亚洲 中文字幕| 午夜激情久久久久久久| 亚洲视频免费观看视频| 咕卡用的链子| 欧美激情高清一区二区三区| 动漫黄色视频在线观看| 久久久久国内视频| 高清欧美精品videossex| 久久毛片免费看一区二区三区| 亚洲国产中文字幕在线视频| 亚洲少妇的诱惑av| 中国美女看黄片| 国产精品一区二区精品视频观看| 人妻 亚洲 视频| 少妇的丰满在线观看| 亚洲欧美激情在线| 亚洲国产精品一区二区三区在线| 又黄又粗又硬又大视频| av视频免费观看在线观看| aaaaa片日本免费| 黑人巨大精品欧美一区二区mp4| 久久精品国产亚洲av高清一级| 亚洲va日本ⅴa欧美va伊人久久| 麻豆av在线久日| 一区二区三区乱码不卡18| 91av网站免费观看| 又紧又爽又黄一区二区| 日韩欧美免费精品| 国产精品一区二区精品视频观看| 国产精品 国内视频| 成人精品一区二区免费| 久久久久精品人妻al黑| 电影成人av| www日本在线高清视频| 国产精品 国内视频| 久久精品亚洲精品国产色婷小说| 91精品三级在线观看| a在线观看视频网站| 老司机在亚洲福利影院| 1024视频免费在线观看| 午夜激情av网站| 久久人妻熟女aⅴ| 久久 成人 亚洲| 国产一区有黄有色的免费视频| 国产精品久久久久久人妻精品电影 | 老汉色∧v一级毛片| 久久香蕉激情| 老司机靠b影院| 免费看十八禁软件| 免费在线观看视频国产中文字幕亚洲| 亚洲精品国产精品久久久不卡| 亚洲情色 制服丝袜| 一进一出好大好爽视频| 色视频在线一区二区三区| 国产91精品成人一区二区三区 | 丰满少妇做爰视频| 捣出白浆h1v1| 免费av中文字幕在线| 国产成人系列免费观看| 成人三级做爰电影| 日本欧美视频一区| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 欧美国产精品一级二级三级| 国产午夜精品久久久久久| 少妇 在线观看| 国产精品成人在线| 丰满饥渴人妻一区二区三| 免费看十八禁软件| 变态另类成人亚洲欧美熟女 | 十八禁人妻一区二区| 免费在线观看视频国产中文字幕亚洲| 国产免费现黄频在线看| 国产aⅴ精品一区二区三区波| 久久久欧美国产精品| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 人人妻,人人澡人人爽秒播| 国产成人影院久久av| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人| 国产老妇伦熟女老妇高清| 久久人妻熟女aⅴ| 免费日韩欧美在线观看| 黄色视频,在线免费观看|