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

    θ-Al2Cu金屬間化合物中氫原子占位的第一性原理研究

    2021-12-23 06:13:46宋莉莉傅高升李龍澤王火生
    機(jī)械工程材料 2021年12期
    關(guān)鍵詞:原子間八面體氫原子

    宋莉莉,傅高升,李龍澤,王火生

    (1.福州大學(xué)材料科學(xué)與工程學(xué)院,福州 350108;2.寧德職業(yè)技術(shù)學(xué)院,福安 355000;3.閩江學(xué)院,福州 350108;4.福建工程學(xué)院材料科學(xué)與工程學(xué)院,福州 350118)

    0 引 言

    鋁合金由于具有高強(qiáng)度、低密度等優(yōu)異的綜合性能而廣泛應(yīng)用于航空航天、交通運(yùn)輸、汽車(chē)工業(yè)等領(lǐng)域。一般來(lái)說(shuō),鋁合金的制備都要經(jīng)過(guò)熔煉過(guò)程,在熔煉過(guò)程中鋁熔體大量吸氣而引入氫,氫在熔體凝固過(guò)程中以微氣孔滯留合金中,對(duì)合金的力學(xué)性能和耐腐蝕性能造成不利影響[1];并且氫還會(huì)通過(guò)擴(kuò)散進(jìn)入合金內(nèi)部而形成固溶體、氫化物或與金屬中第二相交互作用引起氫脆[2-3]。由于氫原子半徑極小,難以通過(guò)試驗(yàn)方法研究氫對(duì)材料的作用機(jī)理,因此,目前主要通過(guò)第一性原理來(lái)研究氫對(duì)金屬及其化合物的影響機(jī)理[4-6]。金屬間化合物θ-Al2Cu是2XXX、7XXX系鋁合金的重要強(qiáng)化相,是在合金鑄造凝固和均勻化過(guò)程中形成的,然而其脆性特征明顯,在較低的應(yīng)力下即可導(dǎo)致微裂紋萌生,進(jìn)而形成孔洞。Al2Cu相對(duì)合金的性能有十分重要的影響[7],然而對(duì)其性能研究的相關(guān)報(bào)道較少。ZHOU等[8]通過(guò)第一性原理對(duì)金屬間化合物Al2Cu晶體的彈性模量進(jìn)行了計(jì)算,發(fā)現(xiàn)Al2Cu在彈性模量上表現(xiàn)為明顯的各向異性。WANG等[9]采用第一性原理研究了金屬摻雜對(duì)Al2Cu力學(xué)性能的影響,發(fā)現(xiàn):鈣和鍶原子易取代Al2Cu化合物中鋁位點(diǎn),提高該化合物強(qiáng)度,降低其延展性;而砷、鎵、鈷和釩原子更容易占據(jù)銅位點(diǎn),提高化合物延展性并降低強(qiáng)度。氫原子的擴(kuò)散勢(shì)必對(duì)θ-Al2Cu相的成核、穩(wěn)定性和性能產(chǎn)生影響,但影響機(jī)制還不清楚。已有研究[9]中的金屬摻雜原子以置換方式進(jìn)入θ-Al2Cu晶體結(jié)構(gòu)而產(chǎn)生作用,而氫原子半徑很小,一般通過(guò)占據(jù)晶胞的間隙位置溶入晶體[10]。作者采用第一性原理方法研究了氫在θ-Al2Cu結(jié)構(gòu)中的占位行為,結(jié)合鍵長(zhǎng)、集居數(shù)的變化,討論了氫占位傾向的電子結(jié)構(gòu)機(jī)制。

    1 模型建立與計(jì)算方法

    1.1 結(jié)構(gòu)模型

    θ-Al2Cu金屬間化合物的晶胞為C16四方結(jié)構(gòu),空間群號(hào)為第140號(hào)空間群I4/MCM,晶格常數(shù)[10]a和b均為0.601 9 nm,c為0.490 8 nm,棱間夾角α,β和γ均為90°。Al2Cu晶胞有12個(gè)原子,其中鋁原子8個(gè),銅原子4個(gè);晶胞中有2種間隙,分別為八面體間隙和四面體間隙。當(dāng)氫原子溶入Al2Cu晶體時(shí),由于原子半徑較小,一般會(huì)占據(jù)晶體間隙位置。采用vesta軟件[11]構(gòu)建無(wú)缺陷以及引入1個(gè)間隙氫原子后Al2Cu晶胞的計(jì)算模型,并對(duì)各原子進(jìn)行標(biāo)識(shí),如圖1所示。將Al2Cu晶胞中八面體間隙結(jié)構(gòu)記為oct間隙,氫在oct間隙中的近鄰原子分別為Al0、Al2、Al3、Al5、Cu8和Cu11。晶胞中的四面體間隙有2種,一種由3個(gè)鋁原子和1個(gè)銅原子構(gòu)成,記為tet1間隙;另一種由2個(gè)鋁原子和2個(gè)銅原子構(gòu)成,記為tet2間隙。氫在tet1間隙中的近鄰原子分別為Al0、Al5、Al6和Cu11,在tet2間隙中的近鄰原子分別為Al2、Al3、Cu8和Cu11。

    圖1 無(wú)缺陷和含1個(gè)間隙氫原子Al2Cu晶胞的計(jì)算模型Fig.1 Calculation model of Al2Cu crystal cells without defects (a)and containing one interstitial hydrogen atom (b-d):(b)with hydrogen at oct interstitial site;(c)with hydrogen at tet1 interstitial site and (d)with hydrogen at tet2 interstitial site

    1.2 計(jì)算方法

    采用CASTEP(Cambridge Sequential Total Energy Package)程序軟件包對(duì)氫摻雜前后的θ-Al2Cu晶胞結(jié)構(gòu)進(jìn)行計(jì)算,該程序使用基于密度反函數(shù)理論(DFT)和平面波贗勢(shì)的從頭算量子力學(xué)方法,系綜能量包含動(dòng)能、靜電能和交換關(guān)聯(lián)能,其中交換關(guān)聯(lián)能采用局域密度近似(LDA)中的CA-PZ形式,勢(shì)函數(shù)選用倒易空間中表述的OTFG超軟贗勢(shì)。采用自洽迭代(SCF)方法進(jìn)行計(jì)算,其中平面波階段能和K點(diǎn)網(wǎng)格對(duì)計(jì)算結(jié)果精確度影響最大,為平衡計(jì)算效率和計(jì)算精度,經(jīng)收斂測(cè)試,設(shè)定平面波截?cái)嗄転?00 eV,第一布里淵區(qū)k-point矢量分別取8×8×10。采用Broyden-Fletcher-Goldfarb-Shanno (BFGS)方案首先對(duì)Al2Cu結(jié)構(gòu)進(jìn)行幾何優(yōu)化,然后摻入1個(gè)氫原子,再對(duì)摻雜后的結(jié)構(gòu)進(jìn)行優(yōu)化,計(jì)算其結(jié)構(gòu)參數(shù)與能量。幾何結(jié)構(gòu)優(yōu)化的收斂指標(biāo):自洽場(chǎng)能量收斂標(biāo)準(zhǔn)為每原子10-6eV,最大內(nèi)應(yīng)力為0.05 GPa,最大原子位移為10-4nm,每個(gè)原子的受力小于0.3 eV·nm-1。

    為了分析氫原子在Al2Cu中的占位傾向,計(jì)算氫原子在θ-Al2Cu不同間隙位置的形成能。無(wú)缺陷Al2Cu和含間隙氫原子Al2Cu晶胞的形成能計(jì)算公式分別為

    Ef=Et(Al2Cu)-{x(Al)Ea(Al)+

    [1-x(Al)]Ea(Cu)}

    (1)

    EFh=Estru-E(H)-E(Al2Cu)

    (2)

    式中:Ef,EFh分別為無(wú)缺陷和含間隙氫原子Al2Cu晶胞的形成能;Et(Al2Cu)為Al2Cu晶胞的能量;Ea(Al),Ea(Cu)分別為穩(wěn)定自由原子鋁和銅的原子能量;x(Al),x(Cu)分別為Al2Cu中鋁和銅的原子分?jǐn)?shù);Estru為氫摻雜Al2Cu晶胞的總能量;E(H)為引入氫原子的單點(diǎn)能;E(Al2Cu)為無(wú)缺陷Al2Cu晶胞的能量。

    理想晶體中原子排列十分規(guī)則且具有周期性,而單個(gè)氫原子不是長(zhǎng)程有序的。為計(jì)算氫原子單點(diǎn)能,需人為地使氫原子有序排列,為此構(gòu)建了一個(gè)尺寸為1 nm×1 nm×1 nm的盒子,并在中心位置放置一個(gè)氫原子以忽略氫原子之間的相互作用,然后對(duì)盒子進(jìn)行充分馳豫,使計(jì)算結(jié)果真實(shí)可信。

    2 計(jì)算結(jié)果及討論

    2.1 晶體結(jié)構(gòu)

    采用第一性原理計(jì)算晶體在基態(tài)時(shí)的能量極小值,從而得到穩(wěn)定的晶體結(jié)構(gòu)。在晶格參數(shù)附近對(duì)晶胞進(jìn)行結(jié)構(gòu)優(yōu)化,得到不同結(jié)構(gòu)的晶胞參數(shù)。

    表1中:V為晶胞體積;ΔV為晶胞體積膨脹率;E為結(jié)構(gòu)優(yōu)化后晶胞的能量。由表1可以看出,氫原子進(jìn)入Al2Cu晶胞間隙后,無(wú)論是占據(jù)四面體間隙還是占據(jù)八面體間隙,均會(huì)引起晶胞結(jié)構(gòu)不同程度的畸變。氫原子溶入導(dǎo)致的體積膨脹率越大,晶胞的畸變程度越大,畸變能越高,結(jié)構(gòu)的穩(wěn)定性越低,溶解度也越小。氫原子占據(jù)oct間隙導(dǎo)致的晶胞體積膨脹率為1.127%,小于氫原子占據(jù)tet1和tet2間隙導(dǎo)致的晶胞體積膨脹率,說(shuō)明氫原子溶入八面體間隙時(shí)晶胞結(jié)構(gòu)更穩(wěn)定。氫原子溶于Al2Cu晶體后,晶胞能量降低,說(shuō)明氫可以穩(wěn)定地存在于Al2Cu晶體中;當(dāng)氫占據(jù)八面體間隙時(shí),晶胞能量最小,說(shuō)明氫更容易占據(jù)Al2Cu晶胞中的八面體間隙。

    表1 無(wú)缺陷和含1個(gè)間隙氫原子Al2Cu的晶胞參數(shù)和能量Table 1 Crystal cell parameters and energy of Al2Cu without defects and containing one interstitial hydrogen atom

    2.2 形成能

    結(jié)構(gòu)充分馳豫后,得到氫原子的單點(diǎn)能為-12.999 9 eV,無(wú)缺陷Al2Cu晶胞的形成能為-359.828 kJ·mol-1,Al2Cu晶胞的形成能與文獻(xiàn)[12]中交換關(guān)聯(lián)勢(shì)采用廣義梯度近似(GGA)下PBE方法處理得到的結(jié)果(-353.181 kJ·mol-1)接近,說(shuō)明結(jié)果可靠。由式(2)計(jì)算得到氫原子占據(jù)oct間隙、tet1間隙和tet2間隙時(shí)的形成能分別為-70.665 6,-3.465 6,10.550 4 kJ·mol-1,換算成電子伏特分別為-0.736 1,-0.036 1,0.109 9 eV。一般而言,如果形成能為正,則需要輸入能量,如果為負(fù)則說(shuō)明反應(yīng)會(huì)自發(fā)進(jìn)行。因此,可采用形成能判斷晶胞的穩(wěn)定性:間隙原子進(jìn)入晶胞的形成能越負(fù),間隙原子越容易進(jìn)入晶胞,晶胞越穩(wěn)定。占據(jù)tet2間隙時(shí)的形成能為正值,說(shuō)明氫原子無(wú)法自發(fā)進(jìn)入tet2間隙,即該摻雜結(jié)構(gòu)不能穩(wěn)定存在;占據(jù)oct間隙時(shí)的形成能為負(fù)值且最小,說(shuō)明氫原子會(huì)優(yōu)先占據(jù)八面體間隙。

    2.3 氫對(duì)θ-Al2Cu晶體電子結(jié)構(gòu)的影響

    為了從電子層面進(jìn)一步解釋氫與近鄰原子之間的相互作用,分別計(jì)算了氫占據(jù)Al2Cu晶體結(jié)構(gòu)不同間隙位置時(shí)的化學(xué)鍵鍵長(zhǎng)、電荷數(shù)和集居數(shù)?;瘜W(xué)鍵的鍵長(zhǎng)是判斷原子間作用力強(qiáng)弱的重要依據(jù),鍵長(zhǎng)越短,鍵能越大,原子間的作用力越強(qiáng)。氫鍵鍵能的大小與成鍵原子的電負(fù)性大小也有關(guān),電負(fù)性越大,氫鍵越強(qiáng),鍵能也越大。在氫原子摻雜的θ-Al2Cu晶胞中,與氫原子近鄰的有鋁和銅原子,鋁原子的電負(fù)性是1.61,銅原子的電負(fù)性是1.9,因此Cu—H鍵的鍵能較Al—H鍵大,Cu—H之間的鍵長(zhǎng)對(duì)鍵能的貢獻(xiàn)起主要作用。由表2可以看出:當(dāng)氫原子位于八面體間隙時(shí),與相鄰鋁原子之間的鍵長(zhǎng)為0.227 8 nm,略大于占據(jù)四面體間隙時(shí)與對(duì)應(yīng)鋁原子之間的鍵長(zhǎng);八面體間隙中氫原子與相鄰銅原子之間的鍵長(zhǎng)為0.145 4 nm,小于四面體間隙中氫原子與對(duì)應(yīng)銅原子之間的鍵長(zhǎng)。因此,氫進(jìn)入八面體間隙更穩(wěn)定。

    表2 含1個(gè)間隙氫原子Al2Cu晶胞中的化學(xué)鍵鍵長(zhǎng)Table 2 Chemical bond length of Al2Cu crystal cells containing one interstitial hydrogen atom nm

    由表3可知,當(dāng)氫原子溶入晶體中,由于異類(lèi)原子之間相互作用,晶胞中的鋁、銅電荷數(shù)發(fā)生變化。與無(wú)缺陷Al2Cu晶胞相比,當(dāng)氫原子占據(jù)oct間隙時(shí),氫原子得到電子,電荷數(shù)下降,與氫原子相鄰的鋁原子電荷數(shù)不變,相鄰的銅原子(Cu8和Cu11)均失去電子,電荷數(shù)增加(-0.76>-0.93),氫與銅原子之間的吸引力增強(qiáng),鍵合作用較強(qiáng);當(dāng)氫原子占據(jù)tet1間隙時(shí),相鄰的Al5原子得到電子,電荷數(shù)減少,氫與鋁原子間吸引力降低,Cu11原子失去電子,電荷數(shù)增加,氫與銅原子間吸引力增大但較氫原子占據(jù)八面體結(jié)構(gòu)時(shí)略弱;當(dāng)氫原子占據(jù)tet2間隙時(shí),氫的電荷數(shù)為-0.25,相鄰的Al3原子得到電子,氫與鋁原子間吸引力降低,Cu8和Cu11原子失去電子,氫與銅原子間的吸引力增強(qiáng)。當(dāng)溶入不同間隙位置時(shí),氫原子均轉(zhuǎn)移了一定數(shù)量的電子。一般而言,電子轉(zhuǎn)移的數(shù)量越多,離子鍵的相互作用就越強(qiáng)。氫原子占據(jù)oct間隙和tet1間隙時(shí)的電子轉(zhuǎn)移數(shù)量相同(電荷數(shù)均為-0.28),高于占據(jù)tet2間隙時(shí);此外,氫占據(jù)tet1間隙時(shí),銅電荷數(shù)增加得比占據(jù)oct間隙時(shí)多,而銅電荷數(shù)的增加會(huì)減弱鋁向銅轉(zhuǎn)移電子的傾向。因此,氫占據(jù)八面體間隙的傾向最大。

    表3 無(wú)缺陷和含1個(gè)間隙氫原子Al2Cu晶胞中相鄰原子的電荷數(shù)Table 3 Neighboring atomic charge number of Al2Cu crystal cell without defects and containing one interstitial hydrogen atom

    可以通過(guò)原子間的Mulliken集居數(shù)來(lái)評(píng)價(jià)兩原子間的成鍵情況。Mulliken集居數(shù)可以對(duì)電荷的得失情況進(jìn)行量化:集居數(shù)為正值,則表示存在共價(jià)鍵,正值越大,則共價(jià)鍵作用越強(qiáng);集居數(shù)為0,則表示存在離子鍵;集居數(shù)為負(fù)值,則表示存在反鍵,即原子間相互排斥,負(fù)值越小,則原子間的排斥越強(qiáng),原子間鍵會(huì)作用就越弱,晶胞的穩(wěn)定性越差。由圖2可以看出,氫占據(jù)oct間隙時(shí),近鄰原子鍵集居數(shù)為正值,鍵長(zhǎng)最短,晶胞穩(wěn)定性最好,占據(jù)tet2間隙時(shí)晶胞穩(wěn)定性次之,占據(jù)tet1間隙時(shí)晶胞穩(wěn)定性最差。

    圖2 無(wú)缺陷及含1個(gè)間隙氫原子Al2Cu不同間隙結(jié)構(gòu)的近鄰原子鍵集居數(shù)和鍵長(zhǎng)Fig.2 Bond population (a,c,e)and bond length (b,d,f)of neighboring atoms in different interstitial structures of Al2Cu without defects and containing one interstitial hydrogen atom:(a-b)oct interstitial structure;(c-d)tet1 interstitial structure and (e-f)tet2 interstitial structure

    3 結(jié) 論

    (1)無(wú)缺陷以及氫占據(jù)八面體間隙、tet1四面體間隙和tet2四面體間隙后,Al2Cu晶胞總能量分別為-7 585.662,-7 599.655,-7 598.955,-7 598.809 eV,氫占據(jù)八面體間隙時(shí)系統(tǒng)能量最低,且引起的晶格畸變最小,故氫進(jìn)入Al2Cu時(shí)優(yōu)先占據(jù)八面體間隙位置。

    (2)氫占據(jù)八面體間隙、tet1四面體間隙和tet2四面體間隙的形成能分別為-0.736 1,-0.036 1,0.109 9 eV,氫占據(jù)八面體間隙需要的形成能最小,故氫更易占據(jù)八面體間隙。

    (3)氫占據(jù)八面體間隙后,近鄰原子鍵集居數(shù)為正值,鍵長(zhǎng)減小,鍵合力增強(qiáng);氫占據(jù)四面體間隙后,近鄰原子之間鍵長(zhǎng)增加,鍵合力降低。故氫處于Al2Cu八面體間隙位置更穩(wěn)定。

    猜你喜歡
    原子間八面體氫原子
    納米八面體二氧化鈦的制備及光催化性能研究
    一維模型氫原子在啁啾激光場(chǎng)中的閾上電離
    數(shù)學(xué)文化原創(chuàng)題(一)
    原子間相互作用勢(shì)對(duì)中Al濃度Ni75AlxV25?x合金沉淀序列的影響?
    與熱庫(kù)耦合的光學(xué)腔內(nèi)三原子間的糾纏動(dòng)力學(xué)?
    關(guān)于氫原子躍遷問(wèn)題的難點(diǎn)突破
    團(tuán)簇Mn3BP的電子自旋密度
    當(dāng)鈣鈦礦八面體成為孤寡老人
    對(duì)氫原子能量吸收的再認(rèn)識(shí)お
    氫原子n=2能級(jí)的二級(jí)斯塔克效應(yīng)
    ponron亚洲| 丰满人妻一区二区三区视频av| 麻豆国产97在线/欧美| 精品久久久久久久末码| 欧美国产日韩亚洲一区| 亚洲专区国产一区二区| 精品国产三级普通话版| 亚洲av第一区精品v没综合| 美女 人体艺术 gogo| 淫妇啪啪啪对白视频| 欧美不卡视频在线免费观看| 中文亚洲av片在线观看爽| 一个人看视频在线观看www免费| 国产成人一区二区在线| 在线播放无遮挡| 黄色丝袜av网址大全| 18禁黄网站禁片午夜丰满| 久久九九热精品免费| 亚洲va在线va天堂va国产| 久久久久久久精品吃奶| 精品久久久久久久末码| www.www免费av| 国模一区二区三区四区视频| 嫁个100分男人电影在线观看| 日韩欧美精品免费久久| 最近最新免费中文字幕在线| 99在线视频只有这里精品首页| 别揉我奶头~嗯~啊~动态视频| 舔av片在线| 91久久精品国产一区二区成人| 麻豆成人午夜福利视频| 亚洲人成网站在线播| 久久精品国产亚洲av天美| 午夜免费男女啪啪视频观看 | 成人亚洲精品av一区二区| 老师上课跳d突然被开到最大视频| av专区在线播放| 国内精品一区二区在线观看| 国产女主播在线喷水免费视频网站 | 精品一区二区三区视频在线| 日本三级黄在线观看| 可以在线观看毛片的网站| 白带黄色成豆腐渣| 日日啪夜夜撸| 亚洲自拍偷在线| 天堂动漫精品| 毛片一级片免费看久久久久 | 天堂网av新在线| 看黄色毛片网站| 天天一区二区日本电影三级| 永久网站在线| 日本与韩国留学比较| 亚州av有码| 亚洲熟妇中文字幕五十中出| 欧美zozozo另类| 国内精品久久久久久久电影| 欧美3d第一页| 99久久成人亚洲精品观看| 欧美最黄视频在线播放免费| 久久久久免费精品人妻一区二区| 日韩国内少妇激情av| 国产高清不卡午夜福利| 成人特级黄色片久久久久久久| 人人妻人人澡欧美一区二区| 国产高清有码在线观看视频| 女的被弄到高潮叫床怎么办 | 校园人妻丝袜中文字幕| 毛片一级片免费看久久久久 | 成年免费大片在线观看| 如何舔出高潮| 亚洲人成网站高清观看| 日韩精品中文字幕看吧| 我的老师免费观看完整版| 五月伊人婷婷丁香| 国产乱人伦免费视频| 久久热精品热| 国产精品无大码| 久久婷婷人人爽人人干人人爱| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 欧美成人免费av一区二区三区| 看免费成人av毛片| 黄色女人牲交| 成人av在线播放网站| av专区在线播放| 91午夜精品亚洲一区二区三区 | 日韩在线高清观看一区二区三区 | 亚洲四区av| 中文字幕久久专区| 国产午夜精品久久久久久一区二区三区 | 久久久久久伊人网av| 18禁裸乳无遮挡免费网站照片| 国产 一区 欧美 日韩| 国产午夜福利久久久久久| 波多野结衣高清作品| 久9热在线精品视频| a级一级毛片免费在线观看| 人人妻,人人澡人人爽秒播| 美女xxoo啪啪120秒动态图| 干丝袜人妻中文字幕| 成人高潮视频无遮挡免费网站| avwww免费| 日本 欧美在线| 中文字幕人妻熟人妻熟丝袜美| 在线天堂最新版资源| 国产成年人精品一区二区| 真人一进一出gif抽搐免费| 国产精品野战在线观看| 亚洲欧美日韩无卡精品| 一个人观看的视频www高清免费观看| 国产精品久久久久久久久免| 国产精品永久免费网站| 在线播放国产精品三级| 成人美女网站在线观看视频| 国产高清三级在线| 91狼人影院| 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 3wmmmm亚洲av在线观看| 少妇的逼好多水| 一个人观看的视频www高清免费观看| 国产美女午夜福利| 亚洲不卡免费看| 免费在线观看日本一区| 亚洲欧美日韩东京热| 12—13女人毛片做爰片一| 免费av毛片视频| 亚洲国产精品成人综合色| 亚洲欧美激情综合另类| 国产精品1区2区在线观看.| 深夜精品福利| 可以在线观看的亚洲视频| 亚洲午夜理论影院| 免费av毛片视频| 欧美色视频一区免费| 国产精品不卡视频一区二区| 国内揄拍国产精品人妻在线| bbb黄色大片| 欧美日韩国产亚洲二区| 亚洲av美国av| 色视频www国产| 国产精品一区二区三区四区免费观看 | 国产精品免费一区二区三区在线| aaaaa片日本免费| 少妇熟女aⅴ在线视频| 99热网站在线观看| 国产麻豆成人av免费视频| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 日本一二三区视频观看| 欧美潮喷喷水| 男人舔女人下体高潮全视频| 国产成人av教育| 特级一级黄色大片| 免费电影在线观看免费观看| 国产69精品久久久久777片| 一个人看的www免费观看视频| 精品乱码久久久久久99久播| 成人鲁丝片一二三区免费| 亚洲av二区三区四区| 国产精品人妻久久久影院| 午夜激情福利司机影院| 国产老妇女一区| 国产精品三级大全| 99热这里只有精品一区| 久久精品人妻少妇| 午夜免费激情av| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 国产视频内射| 欧美三级亚洲精品| 69人妻影院| 久久国产精品人妻蜜桃| 91av网一区二区| 春色校园在线视频观看| 黄色配什么色好看| 亚洲精品国产成人久久av| 欧美成人免费av一区二区三区| 亚洲国产高清在线一区二区三| 亚洲无线在线观看| 久久人妻av系列| 日韩精品青青久久久久久| 欧美激情在线99| 亚洲在线观看片| 最好的美女福利视频网| 18禁裸乳无遮挡免费网站照片| 色综合婷婷激情| avwww免费| 国产真实乱freesex| 国产探花极品一区二区| 精品久久国产蜜桃| 级片在线观看| 欧美日韩精品成人综合77777| 一本一本综合久久| 黄色日韩在线| 黄色视频,在线免费观看| 在线播放国产精品三级| 美女大奶头视频| 国产亚洲av嫩草精品影院| 欧美日本视频| 亚洲熟妇中文字幕五十中出| 男人狂女人下面高潮的视频| 大型黄色视频在线免费观看| 午夜福利视频1000在线观看| 人妻制服诱惑在线中文字幕| 亚洲av成人av| 亚洲不卡免费看| 又粗又爽又猛毛片免费看| 天天躁日日操中文字幕| 色5月婷婷丁香| 亚洲最大成人中文| 变态另类丝袜制服| 老师上课跳d突然被开到最大视频| 人人妻,人人澡人人爽秒播| 国产精品久久视频播放| 亚洲最大成人中文| 黄色视频,在线免费观看| 国产又黄又爽又无遮挡在线| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| 亚洲人成网站在线播| 两个人的视频大全免费| videossex国产| 国产亚洲精品综合一区在线观看| 久久人人爽人人爽人人片va| 国产在视频线在精品| 精品人妻偷拍中文字幕| 日本一二三区视频观看| h日本视频在线播放| 亚洲最大成人av| 国产综合懂色| avwww免费| 少妇人妻精品综合一区二区 | 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 亚洲av免费高清在线观看| 搞女人的毛片| 91久久精品国产一区二区三区| 日本三级黄在线观看| 国产精品久久久久久亚洲av鲁大| 午夜精品在线福利| 一级黄色大片毛片| 欧美zozozo另类| 久久久久久久久中文| 日日撸夜夜添| 日本三级黄在线观看| 国产男靠女视频免费网站| 桃红色精品国产亚洲av| av福利片在线观看| 丝袜美腿在线中文| 看片在线看免费视频| 色综合亚洲欧美另类图片| 日本黄色片子视频| 波多野结衣高清无吗| a在线观看视频网站| 日本熟妇午夜| 18禁在线播放成人免费| 免费在线观看日本一区| 深夜精品福利| 老女人水多毛片| 亚洲精品成人久久久久久| 精品人妻熟女av久视频| 国产精品亚洲一级av第二区| 精品免费久久久久久久清纯| 国产精品不卡视频一区二区| 国产成人福利小说| 国产大屁股一区二区在线视频| 97人妻精品一区二区三区麻豆| 成人av在线播放网站| 三级男女做爰猛烈吃奶摸视频| 乱人视频在线观看| 日韩欧美在线乱码| a级毛片免费高清观看在线播放| 亚洲天堂国产精品一区在线| 日韩精品有码人妻一区| 人妻久久中文字幕网| 久久欧美精品欧美久久欧美| 麻豆成人午夜福利视频| 色哟哟哟哟哟哟| 夜夜爽天天搞| 日本熟妇午夜| 精品免费久久久久久久清纯| 欧美人与善性xxx| 亚洲精品粉嫩美女一区| 人妻丰满熟妇av一区二区三区| 精品国产三级普通话版| 精品久久久久久成人av| 五月伊人婷婷丁香| 中出人妻视频一区二区| 午夜福利18| 51国产日韩欧美| 亚洲国产高清在线一区二区三| 欧美在线一区亚洲| 不卡视频在线观看欧美| 一夜夜www| 日韩欧美在线乱码| 综合色av麻豆| 国产精品久久久久久久久免| 夜夜看夜夜爽夜夜摸| 国产精品人妻久久久影院| 国产成人aa在线观看| 国产亚洲精品综合一区在线观看| 亚洲乱码一区二区免费版| 成人午夜高清在线视频| 亚洲人成网站高清观看| 国产精品,欧美在线| 亚洲美女搞黄在线观看 | 丰满乱子伦码专区| 人妻丰满熟妇av一区二区三区| 久久久久久久久大av| 一卡2卡三卡四卡精品乱码亚洲| 免费看a级黄色片| 欧美色视频一区免费| 日日摸夜夜添夜夜添av毛片 | 在线观看午夜福利视频| 少妇人妻精品综合一区二区 | 久久99热这里只有精品18| 久久九九热精品免费| 在线播放无遮挡| 久久精品国产亚洲av香蕉五月| 99精品久久久久人妻精品| 国产精品亚洲美女久久久| 日韩欧美一区二区三区在线观看| 久9热在线精品视频| 内射极品少妇av片p| 色吧在线观看| 免费在线观看成人毛片| 网址你懂的国产日韩在线| 国产av不卡久久| 国产视频内射| 两性午夜刺激爽爽歪歪视频在线观看| 美女免费视频网站| 一区二区三区免费毛片| 国产大屁股一区二区在线视频| 成人av一区二区三区在线看| 乱人视频在线观看| 日本在线视频免费播放| 日韩人妻高清精品专区| 岛国在线免费视频观看| 九色国产91popny在线| 国产精品久久久久久亚洲av鲁大| 免费人成视频x8x8入口观看| 日韩大尺度精品在线看网址| 亚洲国产精品sss在线观看| 一本一本综合久久| 久久精品国产亚洲av涩爱 | 欧美日韩瑟瑟在线播放| 有码 亚洲区| 精品久久久久久成人av| 老熟妇仑乱视频hdxx| 深爱激情五月婷婷| 国产精品伦人一区二区| 看免费成人av毛片| 97超视频在线观看视频| 国产中年淑女户外野战色| 欧美不卡视频在线免费观看| АⅤ资源中文在线天堂| 极品教师在线视频| 国产极品精品免费视频能看的| 免费大片18禁| 真实男女啪啪啪动态图| 亚洲自拍偷在线| 午夜精品在线福利| 搞女人的毛片| 免费高清视频大片| 校园春色视频在线观看| 97超级碰碰碰精品色视频在线观看| 精品不卡国产一区二区三区| 亚洲av熟女| 国产精品野战在线观看| 久久婷婷人人爽人人干人人爱| 久99久视频精品免费| 国产精品乱码一区二三区的特点| 国产免费av片在线观看野外av| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 久久人人精品亚洲av| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区 | av在线天堂中文字幕| 国产精品一区二区免费欧美| 国产 一区精品| 国产精品美女特级片免费视频播放器| 国产毛片a区久久久久| 国产一区二区激情短视频| 亚洲不卡免费看| 自拍偷自拍亚洲精品老妇| 久久久色成人| 波野结衣二区三区在线| 舔av片在线| 亚洲久久久久久中文字幕| 国产精品综合久久久久久久免费| 观看美女的网站| 久久人妻av系列| 美女免费视频网站| 一进一出抽搐动态| 波野结衣二区三区在线| 美女xxoo啪啪120秒动态图| 黄色欧美视频在线观看| 亚洲精品成人久久久久久| 国产成人aa在线观看| 熟女电影av网| 天美传媒精品一区二区| 可以在线观看毛片的网站| 日韩欧美一区二区三区在线观看| 最新在线观看一区二区三区| 身体一侧抽搐| 最近最新免费中文字幕在线| 国产伦在线观看视频一区| 亚洲熟妇熟女久久| 亚洲图色成人| 久久国产乱子免费精品| 亚洲国产色片| 特级一级黄色大片| 美女黄网站色视频| 成年免费大片在线观看| 国产极品精品免费视频能看的| 欧美在线一区亚洲| 91在线观看av| 波多野结衣高清作品| 亚洲不卡免费看| av国产免费在线观看| 又爽又黄无遮挡网站| 人人妻人人澡欧美一区二区| 亚洲国产色片| 九九在线视频观看精品| 亚洲av熟女| 国产亚洲精品av在线| 国产黄a三级三级三级人| 国产综合懂色| 国产伦精品一区二区三区视频9| 又爽又黄无遮挡网站| 波多野结衣高清无吗| 在现免费观看毛片| 日本成人三级电影网站| 岛国在线免费视频观看| 国产精品不卡视频一区二区| 天堂√8在线中文| 99久久成人亚洲精品观看| 亚洲av中文av极速乱 | 又紧又爽又黄一区二区| 搡老妇女老女人老熟妇| 欧美最黄视频在线播放免费| 亚洲美女黄片视频| 99久久精品国产国产毛片| 久久久久国产精品人妻aⅴ院| a级毛片免费高清观看在线播放| 精品久久久久久久久久免费视频| 国产精品av视频在线免费观看| 国产精品伦人一区二区| 午夜福利在线在线| 乱系列少妇在线播放| 国产探花极品一区二区| 美女高潮喷水抽搐中文字幕| 亚洲图色成人| 婷婷色综合大香蕉| 久久99热6这里只有精品| 免费看a级黄色片| 中文字幕免费在线视频6| 美女黄网站色视频| 麻豆一二三区av精品| 精品免费久久久久久久清纯| 中文亚洲av片在线观看爽| 可以在线观看的亚洲视频| 国产真实伦视频高清在线观看 | 国产国拍精品亚洲av在线观看| 91在线观看av| 国产精品久久视频播放| 在线观看av片永久免费下载| 99久国产av精品| 精品一区二区免费观看| 免费一级毛片在线播放高清视频| 狠狠狠狠99中文字幕| 国产亚洲精品久久久久久毛片| 18禁裸乳无遮挡免费网站照片| 免费av不卡在线播放| 国产私拍福利视频在线观看| 精品久久久久久久末码| 色综合婷婷激情| 在线免费观看的www视频| 免费在线观看成人毛片| 亚洲精品亚洲一区二区| 美女cb高潮喷水在线观看| 少妇高潮的动态图| 亚洲性夜色夜夜综合| 亚洲精华国产精华液的使用体验 | 又爽又黄a免费视频| 午夜a级毛片| 人妻久久中文字幕网| 精品久久久久久久久亚洲 | 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | x7x7x7水蜜桃| 又粗又爽又猛毛片免费看| 国产精品美女特级片免费视频播放器| 老熟妇仑乱视频hdxx| 久99久视频精品免费| 老熟妇仑乱视频hdxx| 成人av在线播放网站| 国产爱豆传媒在线观看| 能在线免费观看的黄片| 成人毛片a级毛片在线播放| 国产真实乱freesex| 日韩人妻高清精品专区| 人妻少妇偷人精品九色| 亚洲美女黄片视频| 九色国产91popny在线| 黄色视频,在线免费观看| 变态另类成人亚洲欧美熟女| 非洲黑人性xxxx精品又粗又长| 久久亚洲精品不卡| 赤兔流量卡办理| 夜夜爽天天搞| 嫩草影视91久久| 日韩欧美在线二视频| 国产一区二区激情短视频| 国产精品国产高清国产av| 麻豆国产av国片精品| 91av网一区二区| 欧美一级a爱片免费观看看| 一个人免费在线观看电影| 亚洲电影在线观看av| 麻豆av噜噜一区二区三区| 国产私拍福利视频在线观看| 在线播放无遮挡| 亚洲精品久久国产高清桃花| 欧美又色又爽又黄视频| 天堂影院成人在线观看| 俄罗斯特黄特色一大片| 日本黄色视频三级网站网址| 最新在线观看一区二区三区| 88av欧美| 国产精品乱码一区二三区的特点| 麻豆国产av国片精品| 亚洲第一电影网av| 亚洲av美国av| 不卡视频在线观看欧美| 最近在线观看免费完整版| 深夜a级毛片| 亚洲在线观看片| 日本a在线网址| 日本黄大片高清| 99久国产av精品| АⅤ资源中文在线天堂| 伦理电影大哥的女人| 国产亚洲精品久久久com| 国内精品美女久久久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲午夜理论影院| 精品福利观看| or卡值多少钱| 日本一本二区三区精品| 国产伦精品一区二区三区视频9| 中文字幕人妻熟人妻熟丝袜美| 精品乱码久久久久久99久播| 亚洲av成人精品一区久久| 精品久久久久久久久av| 尾随美女入室| 免费看美女性在线毛片视频| 国产亚洲精品av在线| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 色精品久久人妻99蜜桃| 亚洲电影在线观看av| 国产精品美女特级片免费视频播放器| 欧美又色又爽又黄视频| 久久草成人影院| 国产高清激情床上av| 十八禁国产超污无遮挡网站| 精品久久国产蜜桃| 搡老熟女国产l中国老女人| 欧美激情在线99| 国产视频一区二区在线看| 在现免费观看毛片| 久久久午夜欧美精品| 精品午夜福利视频在线观看一区| 九九久久精品国产亚洲av麻豆| 男女那种视频在线观看| 亚洲av日韩精品久久久久久密| 久久久久久久久久成人| 色精品久久人妻99蜜桃| 身体一侧抽搐| 久久精品国产亚洲网站| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 麻豆国产97在线/欧美| 欧美+日韩+精品| 最近最新免费中文字幕在线| 最好的美女福利视频网| 大又大粗又爽又黄少妇毛片口| 色播亚洲综合网| av在线观看视频网站免费| 精品一区二区三区av网在线观看| 亚洲黑人精品在线| 日本五十路高清| 99热这里只有是精品50| 久久久久久久精品吃奶| 婷婷精品国产亚洲av在线| 成人美女网站在线观看视频| 免费搜索国产男女视频| 亚洲精品一区av在线观看| 亚洲三级黄色毛片| 欧美成人一区二区免费高清观看| 亚洲欧美日韩无卡精品| a级一级毛片免费在线观看| 俺也久久电影网| 搡老岳熟女国产| 久久久久久久久中文| 国内精品宾馆在线| 婷婷六月久久综合丁香| 内射极品少妇av片p|