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

    納米ZnO改性瀝青分子動力學(xué)模擬研究

    2021-12-17 10:44:30蘇曼曼司春棣張洪亮
    關(guān)鍵詞:回轉(zhuǎn)半徑擴散系數(shù)模量

    蘇曼曼, 司春棣, 張洪亮

    (1. 煙臺大學(xué) 土木工程學(xué)院, 山東 煙臺 264005; 2. 石家莊鐵道大學(xué) 交通運輸學(xué)院, 河北 石家莊 050010; 3. 長安大學(xué) 特殊地區(qū)公路工程教育部重點實驗室, 陜西 西安 710064)

    0 引 言

    為延長瀝青道路的使用壽命并減緩瀝青病害的發(fā)生,改性瀝青在高速公路建設(shè)中被廣泛使用。納米材料擁有顆粒粒徑較小、比表面積大等特性,能改變共混材料的微觀結(jié)構(gòu)。因此,納米材料改性瀝青受到越來越多道路研究者的關(guān)注[1-2]。G.SHAFABAKHSH等[3]、D.WANG等[4]研究了納米TiO2對AH-70基質(zhì)瀝青的影響,研究結(jié)果表明,納米TiO2增大了瀝青的軟化點,并提高瀝青的高溫抗車轍能力,同時納米TiO2具有較好的耐光氧化能力,可提高瀝青的抗氧化和抗老化能力;XU Xu等[5]對納米ZnO熱拌瀝青混合料進行了室內(nèi)試驗研究,研究表明,納米ZnO可改善瀝青的抗老化及抗疲勞性能;H.Y.LIU等[6]則通過室內(nèi)試驗驗證了經(jīng)過表面修飾后的納米ZnO能均勻分散于瀝青中,并提高了瀝青的抗老化性能;ZHANG Hongliang等[7]則通過納米ZnO改性瀝青的三大指標(biāo)試驗及老化試驗等,驗證了納米ZnO能有效改善瀝青的抗老化性能,并通過紅外光譜試驗驗證了納米ZnO與瀝青之間發(fā)生了化學(xué)反應(yīng)。

    隨著對納米改性瀝青研究的不斷深入,越來越多的學(xué)者從不同角度采用不同的手段對納米改性瀝青的改性機理進行研究。ZHU Juncai等[8]、LI Rui等[9]分別通過原子力顯微鏡及掃描電鏡試驗觀察納米黏土、納米ZnO在瀝青中的微觀結(jié)構(gòu),并通過對比基質(zhì)瀝青和改性瀝青的微觀形貌觀察改性劑在瀝青中的分散性;M.ABDELRAHMAN[10]通過紅外光譜試驗發(fā)現(xiàn),納米黏土與瀝青之間產(chǎn)生了化學(xué)反應(yīng),打破了基質(zhì)瀝青的結(jié)構(gòu),同時具有較大比表面積的納米顆粒與瀝青之間產(chǎn)生較強的作用力,導(dǎo)致瀝青具有較好的物理性能;張慶等[11]提出,可以通過差式掃描量熱法研究橡膠改性瀝青在老化過程中的相態(tài)轉(zhuǎn)變溫度,采用熱重分析技術(shù)分析改性瀝青組分,揭示改性瀝青的改性機理;陳淵召等[12]采用組分分離試驗分析了納米ZnO改性瀝青的改性機理,試驗表明,納米ZnO加入瀝青中后,瀝青質(zhì)處于膠束中心,增大了瀝青質(zhì)的含量,進而提高了瀝青的高溫性能。目前對納米改性瀝青的研究大多均停留在宏觀試驗層面,主要研究改性劑對瀝青物理性能改善效果,對改性瀝青的改性機理研究僅限于紅外光譜和掃描電鏡等表象試驗。掃描電鏡和熒光顯微鏡僅能通過影像圖片研究改性劑在瀝青中的分散形貌,但僅根據(jù)分散情況推測瀝青形貌的改變及對瀝青性能的影響。紅外光譜可以判定改性劑與瀝青之間是否出現(xiàn)新的官能團,是否發(fā)生化學(xué)反應(yīng),但不能測定瀝青分子與改性劑之間的相互作用力。因此,紅外光譜和掃描電鏡試驗結(jié)果并不能全面、深入解釋改性劑對瀝青性能的改善原因。

    分子動力學(xué)作為研究材料性質(zhì)、揭示材料力學(xué)行為機理的一種新途徑,在很大程度上起到溝通微觀機理與宏觀力學(xué)特性的橋梁作用。近幾年,越來越多的學(xué)者將分子動力學(xué)研究方法應(yīng)用于改性瀝青的研究中。WANG Peng 等[13]構(gòu)建了瀝青分子模型、碳納米管分子模型并進行分子動力學(xué)模擬,采用鍵接能和非鍵接能表征了各分子體系中的相互作用;SU Manman等[14-15]采用分子動力學(xué)方法研究了SBS改性劑與瀝青之間的相互作用及瀝青物理性能改善效果,進而,對聚合物/納米ZnO共混材料的物理性能及共混體系中瀝青分子結(jié)構(gòu)變化進行了分子力學(xué)計算。筆者在前述研究的基礎(chǔ)上進一步探究納米ZnO與瀝青的相互作用,揭示納米ZnO對瀝青物理性能及分子結(jié)構(gòu)的影響。

    筆者以SK-70瀝青為研究對象,采用MS軟件構(gòu)建瀝青分子模型、納米ZnO模型及納米ZnO/瀝青共混體系,對各體系進行分子動力學(xué)模擬,然后研究納米ZnO分子與瀝青之間的相互作用、納米ZnO在基質(zhì)瀝青中的擴散、納米ZnO對瀝青物理性能及分子結(jié)構(gòu)的影響,通過分子動力學(xué)模擬結(jié)果解釋納米ZnO改性瀝青的改性機理。

    1 計算參數(shù)

    1.1 相互作用能

    采用分子間的非鍵接、范德華和靜電相互作用能作為各體系穩(wěn)定性的評價指標(biāo)。以a、b兩體系為例,相互作用能計算公式如式(1)~式(3)[16]:

    En=Eabn-Ean-Ebn

    (1)

    EV=EabV-EaV-EbV

    (2)

    Eε=Eabε-Eaε-Ebε

    (3)

    式中:En為a、b體系非鍵接相互作用能,kcal/mol;EV為a、b體系范德華相互作用能,kcal/mol;Eε為a、b體系靜電相互作用能,kcal/mol;Eabn、Ean、Ebn分別為ab共混體系、a體系、b體系的非鍵接能,kcal/mol;EabV、EaV、EVb分別為ab共混體系、a體系、b體系的范德華勢能,kcal/mol;Eabε、Eaε、Ebε分別為ab共混體系、a體系、b體系的靜電勢能,kcal/mol。

    1.2 擴散系數(shù)

    納米材料在瀝青中的分散、遷移能力采用擴散系數(shù)進行表征,計算公式如式(4)[17]:

    (4)

    式中:D為擴散系數(shù);t為時間,ps;r(t)為t時刻分子的坐標(biāo);s(t)為分子均方位移,?2;m為均方位移隨時間變化曲線的斜率。

    1.3 物理模量參數(shù)

    在分子模擬計算中,任一受到外力作用的體系均處于應(yīng)力狀態(tài)下,會引起體系內(nèi)粒子相對位置的改變。對于各向同性的材料,其應(yīng)力應(yīng)變行為僅由拉梅常數(shù)便可完全描述,此時體系的剛度矩陣[c]則可通過拉梅常數(shù)建立其與應(yīng)力應(yīng)變之間的關(guān)系式,進而可計算各體系的體積模量K、剪切模量G,如式(5)~(10):

    (5)

    (6)

    (7)

    (8)

    (9)

    (10)

    式中:KH為Voigt-Reuss-Hill法體積模量近似均值;KV為Voigt-Reuss-Hill法體積模量近似上限值;KR為Voigt-Reuss-Hill法體積模量近似下限值;cij(i=1,2,…,6;j=1,2,…,6)為剛度矩陣中的分量值;Sij(i=1,2,…,6;j=1,2,…,6)為柔度矩陣分量值;GH為Hill法剪切模量近似均值;GV為Voigt法剪切模量近似上限值;GR為Reuss法剪切模量近似下限值。

    1.4 結(jié)構(gòu)參數(shù)

    1.4.1 徑向分布函數(shù)

    徑向分布函數(shù)反映分子間相互作用的本質(zhì),可用于分析共混體系微觀結(jié)構(gòu)[18],其表征公式如式(11):

    (11)

    式中:NAB為體系中A、B原子的個數(shù);ΔNAB為對于A原子或者B原子,在r~r+δr范圍內(nèi)出現(xiàn)B原子或者A原子的個數(shù);K為時步,ps;δr為間隔寬度,?;ρAB為體系密度,g/cm3。

    1.4.2 回轉(zhuǎn)半徑

    回轉(zhuǎn)半徑能夠反映分子體積與形狀的動態(tài)變化規(guī)律,可用于評價分子體系的緊密程度和狀態(tài)[19],其計算公式如式(12):

    (12)

    式中:Rg為回轉(zhuǎn)半徑;ri為分子質(zhì)心到第i個鏈單元的距離;mi為第i個鏈單元的分子質(zhì)量。

    2 模型構(gòu)建與模擬方法

    2.1 瀝青分子模型

    SK-70基質(zhì)瀝青廣泛應(yīng)用于中國瀝青路面,其基本性能參數(shù)如表1。

    表1 SK-70基質(zhì)瀝青基本性能Table 1 Basic properties of SK-70 matrix asphalt

    目前,瀝青分子建模的方法有2種:一種為平均分子法,采用核磁共振、紅外光譜等技術(shù)獲取瀝青的平均分子結(jié)構(gòu)并建立瀝青分子模型,此種方法無法獲得瀝青中不同組分化合物與其他物質(zhì)的相互作用;另一種為組裝法,可選定代表性分子代表瀝青各組分,進而組裝成瀝青分子模型,組裝法可更好地了解組分與其他物質(zhì)的相互作用。筆者采用組裝法構(gòu)建瀝青分子模型。

    2.1.1 代表性分子

    參照現(xiàn)有研究成果,選擇瀝青四組分代表性分子(圖1),其中:圖1(a)~(c)為瀝青質(zhì)分子模型;圖1(d)~(h)為膠質(zhì)分子模型;圖1(i)~(k)為芳香分分子模型;圖1(l)~(m)為飽和分分子模型[20-25]。

    圖1 瀝青各組分代表性化合物Fig. 1 Representative compounds of asphalt components

    2.1.2 瀝青分子模型組裝與驗證

    SK-70瀝青分子模型組成如下:瀝青質(zhì)代表性分子a~c的分子個數(shù)分別為2、2、1;膠質(zhì)代表性分子d~h的分子個數(shù)分別為3、6、6、4、9;芳香分代表性分子i~k的分子個數(shù)分別為10、16、15;飽和分代表性分子l~m的分子個數(shù)分別為10、13。各組分與元素含量的計算值與實驗值如表2。顯然,瀝青分子模型各組分與各元素含量的計算值與實驗值均較為接近。

    表2 瀝青組分相對含量及元素試驗值與計算值對比Table 2 Comparison between test and calculated values of relative content of asphalt components and elements

    2.2 納米ZnO簇團模型

    根據(jù)表3中納米ZnO的晶格常數(shù)及空間群號構(gòu)建納米ZnO簇團模型(圖2),簇團直徑分別設(shè)為4、6、8、10 ?。

    圖2 不同粒徑納米ZnO簇團模型Fig. 2 Nano-ZnO cluster model with different particle sizes

    表3 納米ZnO晶格常數(shù)及坐標(biāo)[26]Table 3 Lattice constants and coordinates of ZnO

    2.3 共混體系模型

    采用MS軟件的Amorphous模塊構(gòu)建ZnO/瀝青共混體系模型,各共混體系組成信息如表4,其中粒徑為6 ?的ZnO/瀝青混體系三維模型如圖3。

    表4 各體系組成信息Table 4 Composition information of each system

    圖3 納米ZnO/瀝青共混體系Fig. 3 Nano-ZnO/asphalt blends system

    2.4 模擬方法

    筆者主要研究了納米ZnO與瀝青分子間的相互作用、納米ZnO在瀝青中的擴散、納米ZnO對瀝青分子結(jié)構(gòu)及性能的影響等,分子模擬時選用力場為COMPASS力場,系綜為等溫等壓系綜,具體模擬過程如下:

    1)首先,運用Focite模塊對各體系進行能量優(yōu)化和幾何結(jié)構(gòu)優(yōu)化。

    2)然后,進行退火處理,退火過程采用Amorphous模塊的Protocols程序,溫度為200~450 K,間隔為50 K。

    3)體系結(jié)構(gòu)穩(wěn)定后利用Focite模塊進行分子動力學(xué)模擬,體系分子動力學(xué)計算前后對比見圖4。顯然,分子動力學(xué)模擬后瀝青四組分的排布基本符合膠體模型,即以瀝青質(zhì)為核心、膠質(zhì)包裹著瀝青質(zhì),隨后是芳香分和飽和分,納米ZnO則填充與瀝青分子間隙。

    圖4 ZnO/瀝青共混體系分子動力學(xué)模擬前后對比Fig. 4 Comparison before and after molecular dynamics simulation of ZnO/asphalt blend system

    4)對完成分子動力學(xué)計算的體系進行模擬,并計算物理性能、均方位移、回轉(zhuǎn)半徑等參數(shù)。

    3 結(jié)果與討論

    3.1 ZnO對共混體系相互作用能的影響

    圖5為各共混體系在不同溫度下的非鍵接相互作用能、范德華相互作用能和靜電相互作用能模擬計算結(jié)果。顯然,對任一共混體系,靜電相互作用能幾乎不受溫度的影響,但范德華、非鍵接相互作用能則隨溫度的增長變化較大。對粒徑為4、6 ?納米ZnO/瀝青體系,模擬計算溫度小于410 K時,范德華相互作用能和非鍵接相互作用能隨溫度的增長波動較??;溫度大于410 K時,二者隨溫度的增長波動較為劇烈。簇團粒徑為8 ?時,共混體系范德華相互作用能和非鍵接相互作用能隨溫度的波動幅度最大。簇團粒徑為4、6、8、10 ?的納米/瀝青共混體系的范德華相互作用能和非鍵接相互作用能絕對值分別在溫度為423.3、422.5、418.3、423.2 K時達到最大??梢?,粒徑對納米ZnO/瀝青體系結(jié)構(gòu)最穩(wěn)定對應(yīng)的溫度影響不大,約在150 ℃時各體系分子間相互作用能最大。通常,分子間相互作用能越大,分子結(jié)構(gòu)越穩(wěn)定,從熱力學(xué)觀點可以認為,當(dāng)納米ZnO/瀝青體系結(jié)構(gòu)破壞時需要更多的能量,可以理解為加入ZnO后瀝青結(jié)構(gòu)體系變得更穩(wěn)定了。

    圖5 不同粒徑納米ZnO/瀝青體系相互作用能隨溫度的變化Fig. 5 Variation of interaction between nano-ZnO with different particle sizes chaning with temperature

    3.2 ZnO粒徑對擴散系數(shù)的影響

    圖6為不同納米粒徑ZnO在瀝青中時間-均方位移的變化規(guī)律。

    圖6 ZnO粒徑均方位移隨模擬時間的變化Fig. 6 Variation of mean square displacement of nano-ZnO particle size changing with simulation time

    由圖6可知,不同納米粒徑ZnO/瀝青體系中納米ZnO均方位移隨模擬時間的延長逐漸增大,且納米ZnO粒徑對納米ZnO均方位移有一定影響。

    對各曲線進行一次線性擬合的擬合方程見表5。根據(jù)式(4)可計算粒徑為4、6、8、10 ?時納米ZnO在瀝青中的擴散系數(shù),分別為0.336 8×10-4、0.328 4×10-4、0.293 1×10-4、0.258 0×10-4。粒徑為6、8、10 ?納米ZnO簇團在瀝青中的擴散系數(shù)分別比粒徑為4 ?納米ZnO在瀝青中的擴散系數(shù)降低了2.5%、12.9%和23.4%??梢?,隨著納米ZnO粒徑的增大,納米ZnO粒子在瀝青中的擴散能力減弱。因此,僅從納米顆粒擴散能力角度考慮,在實際工程中應(yīng)選取粒徑較小的納米材料。

    表5 不同粒徑ZnO均方位移擬合式Table 5 MSD formula of nano-ZnO with different particle size

    3.3 ZnO對瀝青物理模量的影響

    圖7為各共混體系剪切模量、體積模量和彈性模量的模擬計算結(jié)果。向瀝青體系中加入納米ZnO顆粒后,瀝青的彈性模量(E)、體積模量(K)和剪切模量(G)均發(fā)生了改變。4、6、8、10 ?納米ZnO/瀝青體系的彈性模量比基質(zhì)瀝青體系的彈性模量分別增長了2.03%、6.27%、6.5%和5.85%;體積模量則較瀝青體系分別增長了15.09%、12.46%、10.06%和8.51%;剪切模量較瀝青體系分別增長了1.33%、1.71%、5.33%和2.21%??梢?,粒徑大于4 ?后的納米ZnO對彈性模量的影響較小;粒徑為8 ?時,彈性模量增長幅度最大,此時剪切模量增長幅度亦達到最大。納米ZnO對瀝青物理性能的改善原因在于納米ZnO顆粒粒徑較小,可以在瀝青分子孔隙中穿越,在一定程度上起到了填充作用,增大了瀝青的體積模量。同時,也在一定程度上提高了瀝青的剪切模量和彈性模量。而剪切模量的提升則意味著瀝青在高溫下的抗剪能力得到提升,從而改善了瀝青的高溫性能。該模擬計算結(jié)論與室內(nèi)試驗結(jié)果一致[7]。

    圖7 瀝青體系與共混體系力學(xué)參數(shù)Fig. 7 Mechanical parameters of asphalt system and blend system

    3.4 ZnO對瀝青分子結(jié)構(gòu)的影響

    取粒徑為8 ?的ZnO簇團與瀝青分子構(gòu)建共混體系,研究ZnO對瀝青分子結(jié)構(gòu)的影響。

    3.4.1 ZnO對瀝青各組分芳環(huán)質(zhì)心徑向分布函數(shù)的影響

    圖8為納米ZnO對瀝青各組分芳環(huán)質(zhì)心徑向分布函數(shù)分子動力學(xué)模擬結(jié)果。不同體系的g(r)隨著原子間距的增加均趨近于1,此為典型的非晶結(jié)構(gòu)特點。未加入ZnO時,瀝青質(zhì)、膠質(zhì)、芳香分和飽和分體系出現(xiàn)第一個峰值的位置分別為1.08、1.12、1.11、1.08 ?;加入ZnO后,瀝青質(zhì)和膠質(zhì)體系第一個峰位置分別右移了0.02、0.03 ?,飽和分體系出現(xiàn)第一個峰值的位置則左移了0.07 ?,芳香分體系第一峰的位置未改變。 未加入納米ZnO簇團時,瀝青質(zhì)、膠質(zhì)、芳香分和飽和分體系的首個峰值分別為41.61、29.03、24.18和41.71;加入納米ZnO后,各體系峰值分別改變?yōu)?5.77、56.51、45.77和57.23,提高量分別為34%、49%、48%和37%。峰值強度提高說明芳環(huán)在該范圍內(nèi)堆積密度提高,而且各體系徑向分布函數(shù)在不同位置的峰值均呈現(xiàn)出高而尖的特點,這表明分子的有序性增強,原子之間聯(lián)系較為緊密。

    圖8 各組分在添加納米ZnO前后的峰值強度比較Fig. 8 Peak intensity for asphalt components before and after adding nano-ZnO

    3.4.2 ZnO對瀝青各組分支鏈回轉(zhuǎn)半徑的影響

    圖9為加入納米ZnO前后瀝青各組分代表性分子支鏈回轉(zhuǎn)半徑變化情況。未加入納米ZnO時,瀝青質(zhì)、膠質(zhì)、芳香分分子側(cè)鏈和飽和分回轉(zhuǎn)半徑值分別在4.32、5.31、5.75、4.81 ?處。加入納米ZnO后,瀝青質(zhì)、飽和分、膠質(zhì)和芳香分支鏈回轉(zhuǎn)半徑峰位分別左移了0.66、0.71、0.21、0.17 ?。未加入納米ZnO時,瀝青質(zhì)、膠質(zhì)、芳香分和飽和分分子支鏈回轉(zhuǎn)半徑峰寬為1.5、1.9、2.6、3.0 ?,加入納米ZnO后,各組分支鏈回轉(zhuǎn)半徑峰寬則變?yōu)?.6、1.9、1.62、2.2 ???梢?,納米ZnO對瀝青質(zhì)峰寬影響較小,對膠質(zhì)峰寬沒有影響,對芳香分和飽和分峰寬影響較大,二者峰寬分別減小了0.58、0.80 ?。各組分的回轉(zhuǎn)半徑可以反映其致密程度,回轉(zhuǎn)半徑峰值左移表明體系發(fā)生了塌縮,體系致密程度增大。峰寬越小則表明支鏈在空間的延展性越強,越容易包裹周圍的分子。

    圖9 不同組分分子支鏈回轉(zhuǎn)半徑計算結(jié)果Fig. 9 Calculation results of molecular branch chain radius of gyration of different components

    由納米ZnO對瀝青分子結(jié)構(gòu)影響的模擬結(jié)果可知,納米ZnO增大了瀝青質(zhì)與膠質(zhì)體系分子間的芳環(huán)質(zhì)心距離,減緩了強極性組分的堆積,同時縮小了芳香分和飽和分分子間芳環(huán)質(zhì)心的距離,加強了支鏈在分子間的延展性。可見,納米ZnO從整體上加強了瀝青各組分之間的交錯,增加了瀝青結(jié)構(gòu)的致密性,促使瀝青形成更穩(wěn)定的膠體結(jié)構(gòu),從而提高了瀝青的物理性能。

    4 結(jié) 論

    1)ZnO粒徑對ZnO與瀝青分子間相互作用能影響較小,溫度約150 ℃時,ZnO與瀝青之間范德華相互作用能和非鍵接相互作用能最大,納米ZnO/瀝青體系結(jié)構(gòu)最穩(wěn)定。

    2)隨著納米ZnO粒徑的增大,納米ZnO在瀝青中的擴散系數(shù)降低。粒徑為6、8、10 ?納米ZnO簇團在瀝青中的擴散系數(shù)分別比粒徑為4 ?納米ZnO在瀝青中的擴散系數(shù)降低了2.5%、12.9%和23.4%。

    3)瀝青體系中加入納米ZnO顆粒后,瀝青的彈性模量、體積模量和剪切模量均出現(xiàn)了不同程度的增長。粒徑為8 ?時,彈性模量和剪切模量增長幅度達到最大。

    4)納米ZnO改變了各組分芳環(huán)質(zhì)心徑向分布函數(shù)首個峰位和峰值強度,增強了分子的有序性,提高了芳環(huán)之間的緊密型,加強了各組分支鏈的延展性,使得各組分之間聯(lián)系較為緊密,促使瀝青形成更穩(wěn)定的膠體結(jié)構(gòu),從而改善了瀝青的物理性能。

    猜你喜歡
    回轉(zhuǎn)半徑擴散系數(shù)模量
    最優(yōu)回轉(zhuǎn)半徑驅(qū)動輪總成的設(shè)計方法研究
    機電信息(2025年2期)2025-01-27 00:00:00
    氨基酸帶電量及分布對蛋白質(zhì)單鏈自組裝行為的影響
    高勁度模量瀝青混合料在京臺高速車轍維修段的應(yīng)用
    室內(nèi)回彈模量和回彈再壓縮模量試驗參數(shù)探討
    山西建筑(2020年11期)2020-06-04 00:09:48
    關(guān)于現(xiàn)行規(guī)范路基頂面回彈模量的理解和應(yīng)用
    上海公路(2018年4期)2018-03-21 05:57:24
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴散系數(shù)測定
    上海金屬(2015年6期)2015-11-29 01:09:09
    非時齊擴散模型中擴散系數(shù)的局部估計
    稻谷堆的壓縮密度與體變模量的測定與分析
    Ni-Te 系統(tǒng)的擴散激活能和擴散系數(shù)研究
    上海金屬(2013年4期)2013-12-20 07:57:07
    欧美日韩国产mv在线观看视频| 夜夜看夜夜爽夜夜摸| 国产精品偷伦视频观看了| 十分钟在线观看高清视频www | 国产在视频线精品| 乱系列少妇在线播放| 大片免费播放器 马上看| 老司机影院毛片| 人妻制服诱惑在线中文字幕| 夫妻性生交免费视频一级片| 国产免费一级a男人的天堂| 街头女战士在线观看网站| 亚洲怡红院男人天堂| 免费av中文字幕在线| 国产免费一区二区三区四区乱码| 日韩 亚洲 欧美在线| 欧美日韩亚洲高清精品| 涩涩av久久男人的天堂| 国产在线男女| 不卡视频在线观看欧美| 看十八女毛片水多多多| 亚洲国产精品一区二区三区在线| 免费播放大片免费观看视频在线观看| 免费久久久久久久精品成人欧美视频 | kizo精华| 精品国产一区二区三区久久久樱花| 日韩av在线免费看完整版不卡| 99久久精品热视频| 亚洲内射少妇av| 啦啦啦视频在线资源免费观看| 中文乱码字字幕精品一区二区三区| 亚洲人与动物交配视频| 简卡轻食公司| 亚洲精品456在线播放app| 国产一区有黄有色的免费视频| av线在线观看网站| 久久狼人影院| 久久av网站| 精品一品国产午夜福利视频| 中国美白少妇内射xxxbb| 欧美日本中文国产一区发布| 边亲边吃奶的免费视频| 一级毛片 在线播放| 91精品国产国语对白视频| 最近2019中文字幕mv第一页| 又黄又爽又刺激的免费视频.| 女人久久www免费人成看片| 国产黄频视频在线观看| 97在线视频观看| 美女cb高潮喷水在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品国产av成人精品| 国产精品久久久久久久久免| 国产午夜精品一二区理论片| 九九在线视频观看精品| 少妇丰满av| 99九九线精品视频在线观看视频| 五月天丁香电影| 精品亚洲成a人片在线观看| a级片在线免费高清观看视频| 插逼视频在线观看| 国产永久视频网站| 一区在线观看完整版| 日韩三级伦理在线观看| 国产黄片视频在线免费观看| 丰满少妇做爰视频| 天堂俺去俺来也www色官网| 精品亚洲成国产av| 亚洲精品一二三| 久久人妻熟女aⅴ| 国产一区二区三区av在线| 18禁裸乳无遮挡动漫免费视频| 精品一区在线观看国产| 欧美日韩视频高清一区二区三区二| 91久久精品国产一区二区成人| 下体分泌物呈黄色| 色5月婷婷丁香| av播播在线观看一区| 综合色丁香网| 国产深夜福利视频在线观看| 在线看a的网站| 九色成人免费人妻av| 国产成人freesex在线| 欧美日韩视频精品一区| 精品少妇内射三级| 一级爰片在线观看| 亚洲精品亚洲一区二区| 久久97久久精品| 久久精品国产亚洲av天美| 国产成人精品福利久久| 精品一区二区免费观看| 国产精品嫩草影院av在线观看| 免费少妇av软件| 日韩人妻高清精品专区| 亚洲国产精品一区三区| 九九爱精品视频在线观看| 黄色日韩在线| 国产男人的电影天堂91| 三级国产精品片| 亚洲,欧美,日韩| 精品久久久久久久久亚洲| 国产成人a∨麻豆精品| 乱人伦中国视频| 九九在线视频观看精品| 国产亚洲91精品色在线| 青春草国产在线视频| 国产精品99久久久久久久久| 国产一区二区在线观看av| 最后的刺客免费高清国语| 国产高清三级在线| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 久热久热在线精品观看| 全区人妻精品视频| 丰满乱子伦码专区| 香蕉精品网在线| 人妻制服诱惑在线中文字幕| 女人精品久久久久毛片| 日韩 亚洲 欧美在线| 中国美白少妇内射xxxbb| 少妇人妻一区二区三区视频| 国产精品久久久久久av不卡| 日韩视频在线欧美| 国产 一区精品| 日日爽夜夜爽网站| 观看免费一级毛片| 亚洲精品日韩在线中文字幕| 美女内射精品一级片tv| 日韩一区二区视频免费看| 久久99精品国语久久久| 丰满少妇做爰视频| 亚洲国产精品999| 哪个播放器可以免费观看大片| 亚洲,一卡二卡三卡| 欧美日韩av久久| 丝瓜视频免费看黄片| 十八禁高潮呻吟视频 | 亚洲国产欧美在线一区| 青春草亚洲视频在线观看| 国产片特级美女逼逼视频| 日本91视频免费播放| 久久久国产精品麻豆| 少妇被粗大猛烈的视频| 人人妻人人添人人爽欧美一区卜| 国产精品熟女久久久久浪| 狠狠精品人妻久久久久久综合| av黄色大香蕉| 国产精品99久久久久久久久| 夫妻性生交免费视频一级片| 日韩一本色道免费dvd| 日韩三级伦理在线观看| 91久久精品国产一区二区成人| 狂野欧美激情性bbbbbb| 国产色婷婷99| 永久网站在线| 久久久久久伊人网av| 丰满迷人的少妇在线观看| 人人妻人人添人人爽欧美一区卜| 中文乱码字字幕精品一区二区三区| 久久久久久久国产电影| 青春草视频在线免费观看| 你懂的网址亚洲精品在线观看| 亚洲美女视频黄频| 久久综合国产亚洲精品| 久久精品国产亚洲网站| 尾随美女入室| 2022亚洲国产成人精品| www.色视频.com| 久久国产亚洲av麻豆专区| 国产美女午夜福利| 国产色爽女视频免费观看| 国产成人freesex在线| 国内少妇人妻偷人精品xxx网站| 国产熟女午夜一区二区三区 | 少妇人妻 视频| 黑人猛操日本美女一级片| 老司机亚洲免费影院| 中文字幕人妻熟人妻熟丝袜美| 高清黄色对白视频在线免费看 | 夫妻性生交免费视频一级片| 亚洲内射少妇av| 亚洲精品456在线播放app| 国产av码专区亚洲av| 欧美另类一区| 老女人水多毛片| 亚洲激情五月婷婷啪啪| 日韩免费高清中文字幕av| 国产精品人妻久久久久久| 看非洲黑人一级黄片| 99久久精品国产国产毛片| 国产精品不卡视频一区二区| 亚洲精品色激情综合| 久久精品久久久久久久性| 人体艺术视频欧美日本| 色视频在线一区二区三区| 国产亚洲av片在线观看秒播厂| 2022亚洲国产成人精品| 一本一本综合久久| 性高湖久久久久久久久免费观看| 99九九线精品视频在线观看视频| 精品视频人人做人人爽| 精品一区二区免费观看| 嘟嘟电影网在线观看| 精品熟女少妇av免费看| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 在线观看三级黄色| a级毛色黄片| 在线精品无人区一区二区三| 丝袜喷水一区| 久热这里只有精品99| 日韩熟女老妇一区二区性免费视频| 在线观看免费日韩欧美大片 | 亚洲欧美精品专区久久| 免费av中文字幕在线| 伊人久久国产一区二区| www.av在线官网国产| 久久精品久久久久久噜噜老黄| 亚洲国产精品一区二区三区在线| 国产色爽女视频免费观看| 成人国产av品久久久| 能在线免费看毛片的网站| 亚州av有码| 噜噜噜噜噜久久久久久91| 久久97久久精品| 大片免费播放器 马上看| 成人18禁高潮啪啪吃奶动态图 | .国产精品久久| 日本免费在线观看一区| 精品熟女少妇av免费看| 精品熟女少妇av免费看| 男人舔奶头视频| 97超视频在线观看视频| 91久久精品国产一区二区三区| 久久99蜜桃精品久久| 成年美女黄网站色视频大全免费 | 亚洲精品,欧美精品| 99热网站在线观看| 色哟哟·www| 日日摸夜夜添夜夜爱| 国产精品嫩草影院av在线观看| av国产精品久久久久影院| 国产av国产精品国产| 久久久久精品性色| 亚洲人成网站在线播| 国产在线男女| 日韩欧美 国产精品| 日本爱情动作片www.在线观看| 又黄又爽又刺激的免费视频.| 国产黄频视频在线观看| 成年人午夜在线观看视频| 久久精品国产亚洲av涩爱| 各种免费的搞黄视频| 国产亚洲一区二区精品| 日本黄大片高清| 我的老师免费观看完整版| 午夜久久久在线观看| 亚洲精品456在线播放app| 天堂中文最新版在线下载| 天堂8中文在线网| 日本欧美视频一区| 日韩制服骚丝袜av| 国产精品久久久久久久久免| 亚洲精品视频女| 国产精品久久久久久av不卡| 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 欧美xxxx性猛交bbbb| 色94色欧美一区二区| 欧美另类一区| 久久久a久久爽久久v久久| 国产熟女午夜一区二区三区 | 欧美精品高潮呻吟av久久| 国产亚洲5aaaaa淫片| 春色校园在线视频观看| 国产白丝娇喘喷水9色精品| 五月玫瑰六月丁香| 我要看黄色一级片免费的| 亚洲欧洲日产国产| 欧美日韩在线观看h| 国产日韩一区二区三区精品不卡 | 日韩av免费高清视频| 51国产日韩欧美| 国产伦在线观看视频一区| av免费观看日本| 一二三四中文在线观看免费高清| 国产精品一二三区在线看| 欧美97在线视频| 天天操日日干夜夜撸| av天堂久久9| 亚洲人成网站在线播| 欧美xxⅹ黑人| 色哟哟·www| 久久午夜福利片| 国产无遮挡羞羞视频在线观看| 久久毛片免费看一区二区三区| 久久人妻熟女aⅴ| av网站免费在线观看视频| 99热国产这里只有精品6| 午夜激情久久久久久久| 最后的刺客免费高清国语| 国产爽快片一区二区三区| 不卡视频在线观看欧美| 免费人妻精品一区二区三区视频| 一边亲一边摸免费视频| 久热久热在线精品观看| 免费高清在线观看视频在线观看| 国产精品99久久99久久久不卡 | 中国三级夫妇交换| av线在线观看网站| 少妇丰满av| 99热这里只有是精品在线观看| 在线看a的网站| 精品卡一卡二卡四卡免费| 人妻一区二区av| 男女国产视频网站| 精品一区二区三卡| 国产精品人妻久久久影院| 内地一区二区视频在线| 精品久久久久久久久亚洲| 建设人人有责人人尽责人人享有的| 寂寞人妻少妇视频99o| 美女内射精品一级片tv| 亚洲中文av在线| 亚洲情色 制服丝袜| 久久这里有精品视频免费| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产毛片av蜜桃av| 亚洲自偷自拍三级| 国产色婷婷99| 九色成人免费人妻av| 精品亚洲成国产av| 99热网站在线观看| 国产亚洲av片在线观看秒播厂| 久久精品久久久久久噜噜老黄| 亚洲精品中文字幕在线视频 | av不卡在线播放| 波野结衣二区三区在线| 毛片一级片免费看久久久久| 五月开心婷婷网| 欧美日韩视频高清一区二区三区二| 两个人免费观看高清视频 | 18+在线观看网站| 2018国产大陆天天弄谢| 精品一区二区三卡| 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 成人国产av品久久久| 在线精品无人区一区二区三| 美女xxoo啪啪120秒动态图| 人妻 亚洲 视频| 亚洲精品色激情综合| 日日啪夜夜爽| 日日撸夜夜添| 国产一区有黄有色的免费视频| 99久久精品热视频| 丁香六月天网| 亚洲精品色激情综合| 国产色爽女视频免费观看| 我的女老师完整版在线观看| 丁香六月天网| 国产男人的电影天堂91| 日韩视频在线欧美| 国产亚洲欧美精品永久| 免费人成在线观看视频色| 亚州av有码| 国产色爽女视频免费观看| 桃花免费在线播放| 成人美女网站在线观看视频| 九色成人免费人妻av| 久久精品熟女亚洲av麻豆精品| 黑人高潮一二区| 极品少妇高潮喷水抽搐| 七月丁香在线播放| tube8黄色片| 女人精品久久久久毛片| 精品人妻一区二区三区麻豆| av免费在线看不卡| 亚洲精品一区蜜桃| 久久99精品国语久久久| 国产淫片久久久久久久久| 寂寞人妻少妇视频99o| av国产精品久久久久影院| 中文字幕久久专区| 2018国产大陆天天弄谢| 人人妻人人看人人澡| 男女国产视频网站| 人人妻人人爽人人添夜夜欢视频 | 久热这里只有精品99| 国产亚洲欧美精品永久| 国产午夜精品一二区理论片| 噜噜噜噜噜久久久久久91| 国产精品不卡视频一区二区| 国产综合精华液| 国产成人freesex在线| 午夜免费鲁丝| 国产成人精品福利久久| av卡一久久| 黄色日韩在线| 婷婷色综合www| 亚洲av日韩在线播放| 丝袜脚勾引网站| 永久网站在线| 简卡轻食公司| freevideosex欧美| 成人二区视频| 亚洲国产av新网站| 免费少妇av软件| 男女边吃奶边做爰视频| 成人午夜精彩视频在线观看| 青春草亚洲视频在线观看| 久久亚洲国产成人精品v| 天美传媒精品一区二区| 伦理电影大哥的女人| 国产91av在线免费观看| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 久久久久视频综合| 精品视频人人做人人爽| 美女xxoo啪啪120秒动态图| 国产在线视频一区二区| 91精品伊人久久大香线蕉| 欧美最新免费一区二区三区| 中文字幕精品免费在线观看视频 | 久久久久久久久久人人人人人人| 春色校园在线视频观看| 日韩一本色道免费dvd| 国产视频内射| 哪个播放器可以免费观看大片| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 少妇精品久久久久久久| 男女国产视频网站| 国产欧美亚洲国产| 最后的刺客免费高清国语| 我的老师免费观看完整版| 免费黄色在线免费观看| 少妇人妻一区二区三区视频| av黄色大香蕉| 高清av免费在线| 久久国产精品男人的天堂亚洲 | 我的女老师完整版在线观看| 久久这里有精品视频免费| 成人漫画全彩无遮挡| www.av在线官网国产| 热re99久久国产66热| 欧美激情极品国产一区二区三区 | 亚洲欧洲精品一区二区精品久久久 | 亚洲av成人精品一二三区| 日本黄色片子视频| 日本av免费视频播放| 熟女人妻精品中文字幕| 精品熟女少妇av免费看| 亚洲va在线va天堂va国产| 五月开心婷婷网| 国产av国产精品国产| 80岁老熟妇乱子伦牲交| 成人毛片60女人毛片免费| 一级片'在线观看视频| 日本91视频免费播放| 日本午夜av视频| 三级国产精品片| 两个人的视频大全免费| 男男h啪啪无遮挡| 高清视频免费观看一区二区| 成人特级av手机在线观看| h视频一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲美女搞黄在线观看| 欧美精品一区二区免费开放| 精品一区二区三区视频在线| 国产免费福利视频在线观看| 成人无遮挡网站| 精品酒店卫生间| 高清不卡的av网站| 久久久久国产网址| 亚洲精品久久午夜乱码| 人妻一区二区av| 国产日韩一区二区三区精品不卡 | 亚洲av不卡在线观看| 国产午夜精品久久久久久一区二区三区| 欧美xxxx性猛交bbbb| 天天操日日干夜夜撸| 国产黄色视频一区二区在线观看| 少妇人妻一区二区三区视频| 久久久a久久爽久久v久久| av一本久久久久| 日本黄色日本黄色录像| 日韩免费高清中文字幕av| 亚洲精品自拍成人| 中文字幕制服av| 五月玫瑰六月丁香| 久久久久久久久久成人| 我要看黄色一级片免费的| 蜜臀久久99精品久久宅男| 久久久久久久国产电影| 大片电影免费在线观看免费| 人人妻人人爽人人添夜夜欢视频 | 国产精品一区二区三区四区免费观看| 夜夜骑夜夜射夜夜干| 日韩av在线免费看完整版不卡| 亚洲精品456在线播放app| 国产精品国产三级国产av玫瑰| kizo精华| 精品久久久精品久久久| 亚洲精品乱码久久久久久按摩| 香蕉精品网在线| 18禁裸乳无遮挡动漫免费视频| 国产午夜精品一二区理论片| av卡一久久| 欧美精品人与动牲交sv欧美| 极品教师在线视频| 久久久久视频综合| 99re6热这里在线精品视频| 最近手机中文字幕大全| 日产精品乱码卡一卡2卡三| 成年av动漫网址| 午夜久久久在线观看| www.av在线官网国产| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 在线观看一区二区三区激情| 国产精品女同一区二区软件| 99精国产麻豆久久婷婷| 99热这里只有是精品在线观看| 久久久久久久久大av| kizo精华| 午夜福利影视在线免费观看| 久久99蜜桃精品久久| 丝袜喷水一区| 18禁在线播放成人免费| 欧美日韩一区二区视频在线观看视频在线| 日本黄色片子视频| 亚洲精品色激情综合| 99国产精品免费福利视频| 日韩免费高清中文字幕av| 在线观看免费日韩欧美大片 | 久久人妻熟女aⅴ| 亚洲av综合色区一区| 大码成人一级视频| 国产精品99久久久久久久久| 我要看黄色一级片免费的| 国产精品无大码| 久久国产精品男人的天堂亚洲 | 国产日韩一区二区三区精品不卡 | 在线观看免费日韩欧美大片 | 纵有疾风起免费观看全集完整版| 国产一级毛片在线| 免费久久久久久久精品成人欧美视频 | 国产黄片视频在线免费观看| 亚洲欧美精品自产自拍| 纯流量卡能插随身wifi吗| 国产亚洲午夜精品一区二区久久| 女人久久www免费人成看片| av女优亚洲男人天堂| 日韩中字成人| 一级毛片aaaaaa免费看小| 一区二区三区四区激情视频| 永久免费av网站大全| 久久久久久人妻| 夫妻午夜视频| 成人亚洲精品一区在线观看| 99热6这里只有精品| 色吧在线观看| 麻豆成人av视频| 久久久久网色| 少妇人妻 视频| 国产午夜精品一二区理论片| 插阴视频在线观看视频| 亚洲欧美精品专区久久| 国产色婷婷99| 精品少妇内射三级| 久久久午夜欧美精品| 国产色爽女视频免费观看| 国产黄片视频在线免费观看| 久久热精品热| 丁香六月天网| 久久青草综合色| 永久免费av网站大全| 亚洲图色成人| 青春草亚洲视频在线观看| 午夜福利在线观看免费完整高清在| 女的被弄到高潮叫床怎么办| 丰满乱子伦码专区| 色网站视频免费| 男人狂女人下面高潮的视频| 色哟哟·www| 欧美区成人在线视频| 欧美丝袜亚洲另类| 免费久久久久久久精品成人欧美视频 | 国产片特级美女逼逼视频| 草草在线视频免费看| 日韩在线高清观看一区二区三区| 2021少妇久久久久久久久久久| 制服丝袜香蕉在线| 99热这里只有精品一区| 精品一区在线观看国产| 九九在线视频观看精品| 婷婷色av中文字幕| 一级a做视频免费观看| 美女脱内裤让男人舔精品视频| 丝袜在线中文字幕| 丰满乱子伦码专区| 久久人妻熟女aⅴ| 内地一区二区视频在线| 久久综合国产亚洲精品| 伊人亚洲综合成人网| 丁香六月天网| 久久6这里有精品| 欧美亚洲 丝袜 人妻 在线| 99热这里只有是精品50| 日日爽夜夜爽网站| 亚洲一区二区三区欧美精品| 观看免费一级毛片| 91精品国产国语对白视频|