亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        PAN纖維炭化過(guò)程中缺陷結(jié)構(gòu)的溫度效應(yīng)

        2015-03-03 07:32:45王夢(mèng)梵蘇世州徐樑華曹維宇
        材料工程 2015年6期
        關(guān)鍵詞:回轉(zhuǎn)半徑炭化微孔

        王夢(mèng)梵,陳 旺,蘇世州,韓 旭,徐樑華,曹維宇

        (1 北京化工大學(xué) 有機(jī)無(wú)機(jī)復(fù)合材料國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100029;2 北京化工大學(xué) 碳纖維及功能高分子教育部重點(diǎn)實(shí)驗(yàn)室,北京 100029)

        ?

        PAN纖維炭化過(guò)程中缺陷結(jié)構(gòu)的溫度效應(yīng)

        王夢(mèng)梵1,陳 旺2,蘇世州2,韓 旭1,徐樑華2,曹維宇1

        (1 北京化工大學(xué) 有機(jī)無(wú)機(jī)復(fù)合材料國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100029;2 北京化工大學(xué) 碳纖維及功能高分子教育部重點(diǎn)實(shí)驗(yàn)室,北京 100029)

        采用小角X射線散射(SAXS)方法研究了PAN纖維在炭化及石墨化過(guò)程中微孔缺陷結(jié)構(gòu)隨熱處理溫度的變化。由微孔結(jié)構(gòu)散射形成的PAN基碳纖維SAXS散射花樣經(jīng)圓積分處理后為光滑曲線。使用一維散射曲線計(jì)算得到微孔結(jié)構(gòu)均方回轉(zhuǎn)半徑、相對(duì)孔體積、孔隙率、比表面積結(jié)構(gòu)參數(shù)隨溫度的變化規(guī)律。結(jié)果表明:隨熱處理溫度提高PAN基碳纖維微孔結(jié)構(gòu)變化分為三個(gè)階段:400~700℃微孔體積變化(孔的融合)為主,形態(tài)變化為輔;700~1800℃微孔體積變化(孔的分裂與縮小)與形態(tài)變化(邊緣復(fù)雜化)基本同步;而1800~2400℃微孔形態(tài)變化為主,體積變化為輔。

        PAN基碳纖維;微孔缺陷結(jié)構(gòu);小角X射線散射(SAXS);溫度效應(yīng)

        纖維內(nèi)部微孔缺陷結(jié)構(gòu)定量分析的常規(guī)方法主要有氣體吸附法、壓汞法及小角X射線散射法(SAXS)等。傳統(tǒng)的氣體吸附法和壓汞法多用于對(duì)活性碳纖維(ACF)孔隙結(jié)構(gòu)的表征, 但無(wú)法對(duì)孔隙含量少、孔隙尺寸小且存在封閉孔的纖維樣品進(jìn)行定量的表征。SAXS法具有測(cè)試時(shí)間短、數(shù)據(jù)易于統(tǒng)計(jì)分析、可靠性高等特點(diǎn), 并且由于X射線可以穿透樣品,所以不受纖維內(nèi)部開(kāi)孔閉孔的影響,其所得到的結(jié)果更具有統(tǒng)計(jì)與實(shí)際意義,所以SAXS方法是研究PAN基碳纖維納米尺度微孔結(jié)構(gòu)最為有效的方法[5]。

        至今國(guó)內(nèi)外使用SAXS法對(duì)PAN基碳纖維的微孔缺陷結(jié)構(gòu)開(kāi)展了研究,主要包括微孔尺寸、長(zhǎng)徑比、界面分形維數(shù),取向等[5-13]。盛毅等[10]使用SAXS法對(duì)T300,T800和自制MHS碳纖維中的微孔結(jié)構(gòu)進(jìn)行對(duì)比,結(jié)果表明,纖維中微孔結(jié)構(gòu)基本沿纖維軸方向取向,MHS碳纖維的界面復(fù)雜程度介于T300和T800之間。Thünemann等[5]使用SAXS研究了預(yù)氧化過(guò)程中微孔結(jié)構(gòu)的變化,結(jié)果表明,預(yù)氧化階段微孔含量較原絲有所增大,并向針狀形態(tài)發(fā)展。

        碳纖維的微孔結(jié)構(gòu)來(lái)源于濕法紡絲過(guò)程中產(chǎn)生的孔洞,在炭化階段隨著熱處理溫度的升高,由于非碳元素的脫除反應(yīng)和類石墨結(jié)構(gòu)高溫重排的發(fā)生,纖維中微孔缺陷結(jié)構(gòu)也會(huì)相應(yīng)發(fā)展,但目前對(duì)此針對(duì)此過(guò)程進(jìn)行的完整研究并不多見(jiàn)。本工作利用SAXS方法研究了PAN基碳纖維炭化(包括石墨化階段)微孔結(jié)構(gòu)隨熱處理溫度的演變規(guī)律,并根據(jù)各結(jié)構(gòu)參數(shù)變化提出了微孔結(jié)構(gòu)的演變模型。

        1 實(shí)驗(yàn)

        1.1 PAN碳纖維的制備

        實(shí)驗(yàn)采用國(guó)產(chǎn)三元共聚PAN原絲, 共聚組成為丙烯腈、衣康酸和丙烯酸甲酯, 單體質(zhì)量比為98∶1∶1。原絲經(jīng)空氣氧化后,在高純氮?dú)鈿夥毡Wo(hù)下,以一定的速率連續(xù)通過(guò)炭化爐及石墨化爐,如圖1所示(圖中溫度表示不同熱處理溫度樣品的收集位置),為了保證研究變量的單一性,走絲過(guò)程中纖維始終處于零張力狀態(tài),收取400~2400℃不同熱處理溫度下的樣品進(jìn)行SAXS測(cè)試。

        圖1 樣品采集示意圖Fig.1 Schematic of the sample collection

        1.2 小角X射線散射(SAXS)

        測(cè)試在上海同步輻射光源(SSRF)BL16B1小角散射實(shí)驗(yàn)站上完成,貯存環(huán)的電子能量為3.5GeV。將0.5mm厚的纖維束平行、緊湊、均勻、垂直地固定在樣品架上進(jìn)行SAXS測(cè)試,檢測(cè)器為Mar165 CCD,曝光時(shí)間為25s。入射X射線的波長(zhǎng)為0.124nm,用標(biāo)樣Collagen(長(zhǎng)周期為66.7nm)標(biāo)定樣品到檢測(cè)器的距離為5220mm。

        2 結(jié)果與討論

        2.1 碳纖維微孔結(jié)構(gòu)的SAXS散射

        根據(jù)SAXS基本原理,散射花樣是纖維內(nèi)部不同微結(jié)構(gòu)的電子密度起伏經(jīng)過(guò)傅里葉變換得到的[14,15],在本文中則主要為微孔結(jié)構(gòu)與纖維基體之間的電子密度差異。一般認(rèn)為碳纖維中微孔結(jié)構(gòu)模型是基于稀粒子體系,具有針狀結(jié)構(gòu),也就是具有較大的長(zhǎng)徑比,且存在沿纖維軸方向的取向結(jié)構(gòu)[5],如圖2所示。根據(jù)倒易關(guān)系,倒易空間與實(shí)空間中的尺寸成反比,在實(shí)空間中微孔的長(zhǎng)軸的在倒易空間中為短軸,短軸則為長(zhǎng)軸。由于纖維中微孔沿著纖維軸有一定的夾角,SAXS散射花樣在赤道方向呈扇形分布,如圖3(a)所示。

        圖2 碳纖維中微孔結(jié)構(gòu)示意圖Fig.2 Model of micro-voids within PAN based carbon fiber

        圖3 碳纖維2D-SAXS散射花樣(圓圈表示不同散射矢量q的方位角掃描)(a)不同溫度碳纖維散射曲線(b)Fig.3 Typical 2D-SAXS pattern of PAN based carbon fiber (where the circle is the schematic of azimuthally scan of the various scattering vector q) (a) and the scattering profiles of carbon fiber at different heat treatment temperatures (b)

        在對(duì)二維散射圖樣進(jìn)行處理時(shí)使用fit-2D軟件,采用沿2θ方向的圓積分獲得纖維中微孔結(jié)構(gòu)的統(tǒng)計(jì)信息。圖3(a)與3(b)為碳纖維SAXS散射圖樣圓積分法示意圖及使用此方法所得不同溫度下的散射曲線。從散射曲線看出,對(duì)于圓積分所得散射曲線,其均為一條I隨q(散射矢量,q=4πsinθ/λ)減小的光滑曲線,結(jié)合碳纖維的特點(diǎn),可判定其均為微孔結(jié)構(gòu)散射所致的信號(hào)。

        2.2 炭化石墨化過(guò)程中孔結(jié)構(gòu)參數(shù)的變化

        根據(jù)圓積分所得到的SAXS數(shù)據(jù),通過(guò)相應(yīng)的計(jì)算方法,本文計(jì)算了炭化及石墨化過(guò)程中微孔結(jié)構(gòu)的均方回轉(zhuǎn)半徑Rg、相對(duì)孔體積Vrel、孔隙率P、比表面積S/V變化,以此來(lái)描述這一過(guò)程中微孔結(jié)構(gòu)的演變規(guī)律。

        2.2.1 相對(duì)孔體積

        在計(jì)算相對(duì)孔體積之前,需要計(jì)算兩個(gè)特殊的散射值:零角散射強(qiáng)度I(0),與積分不變量Q。其中I(0)因?yàn)樵跍y(cè)試時(shí)beamstopper的存在,所以在散射曲線中無(wú)法直接得到q→0nm-1時(shí)的散射強(qiáng)度,需要通過(guò)對(duì)散射曲線中較低散射矢量q的部分使用外推法得到;積分不變量Q值可通過(guò)公式(1)計(jì)算得到。圖4為兩參數(shù)計(jì)算方法示意圖。

        (1)

        圖4 I(0)(a)與Q(b)計(jì)算示意圖Fig.4 Schematic of the calculation of I (0) (a) and Q (b)

        對(duì)于散射曲線,存在以下關(guān)系式:

        (2)

        式中:V為粒子體積;r為相關(guān)長(zhǎng)度;γ(r)為相關(guān)函數(shù);Δρ為體系內(nèi)部?jī)上嗟碾娮用芏炔睢?/p>

        將(1),(2)式聯(lián)立,可得

        (3)

        為了便于比較,本文使用相對(duì)孔體積Vrel代替總孔體積,將所得孔體積與400℃熱處理碳纖維的孔體積作比值處理以表征相對(duì)孔體積Vrel隨溫度的變化規(guī)律。

        2.2.2 均方回轉(zhuǎn)半徑

        均方回轉(zhuǎn)半徑的計(jì)算主要使用Guinier方法[15],如式(4)所示:

        (4)

        圖5 SAXS散射曲線Guinier法擬合示意圖Fig.5 Schematic of SAXS data processing by Guinier method

        2.2.3 孔隙率

        對(duì)于兩相體系來(lái)說(shuō),SAXS散射強(qiáng)度不僅與體系的電子密度漲落有關(guān),也與體系中兩相所占比例有關(guān),如式(5)所示:

        (5)

        (6)

        因?yàn)樘祭w維中微孔結(jié)構(gòu)為稀粒子體系,微孔所占體積比例相對(duì)基體很小,所以可以做以下假設(shè):φv?1,φm≈1。取相關(guān)長(zhǎng)度為r=0時(shí),其自相關(guān)函數(shù)值為1:r=0,γ0(0)=1。由此可得纖維中微孔結(jié)構(gòu)所占體積分?jǐn)?shù)也就是孔隙率:

        (7)

        式中φv即為孔隙率P。

        2.2.4 比表面積

        比表面積計(jì)算使用Porod法[16],Porod法主要研究散射曲線末端所代表得結(jié)構(gòu)信息,因?yàn)閝與r成反比, 當(dāng)q很大時(shí),r很小所以只有微孔結(jié)構(gòu)表面上的r對(duì)γ0(r)有影響,故曲線末端描述的是粒子的表面結(jié)構(gòu)。

        在大散射矢量q區(qū)域,其散射強(qiáng)度與粒子表面積成趨近關(guān)系:

        (8)

        將式(8)與式(4)積分不變量作比值,可得:

        (9)

        可以看出,I(q)/Q→與S/V二者也呈趨近關(guān)系,所以比表面積可用式(10)求得:

        (10)

        2.3 炭化過(guò)程中纖維微孔結(jié)構(gòu)的演變模型

        通過(guò)以上計(jì)算,可得到碳纖維中微孔結(jié)構(gòu)參數(shù)隨溫度的變化,如圖6所示。根據(jù)四條曲線的規(guī)律,可將微孔結(jié)構(gòu)隨溫度的變化分為三個(gè)階段討論:第Ⅰ階段(400~700℃),為低溫炭化階段,此階段相對(duì)孔體積Vrel、孔隙率P、比表面積S/V基本保持不變,而均方回轉(zhuǎn)半徑Rg呈增大趨勢(shì)。說(shuō)明纖維中微孔總體積沒(méi)有變,但是單孔的統(tǒng)計(jì)平均直徑變大,這可能是由于幾個(gè)小孔融并成大孔而導(dǎo)致的,所以只有Rg增大而其他參數(shù)保持不變。雖然在此階段也發(fā)生大量的化學(xué)反應(yīng),并產(chǎn)生小分子氣體,但是由于纖維內(nèi)部結(jié)構(gòu)仍相對(duì)疏松,小分子氣體很容易從纖維中逸出,所以其對(duì)孔總體積的變化幾乎沒(méi)有貢獻(xiàn)。第Ⅱ階段(700~1800℃),為高溫炭化與初級(jí)石墨化階段。此階段相對(duì)孔體積Vrel、均方回轉(zhuǎn)半徑Rg呈減小趨勢(shì),而孔隙率P、比表面積S/V呈增大趨勢(shì)。說(shuō)明在此階段孔的總體積在減小,這是因?yàn)榇穗A段纖維內(nèi)部裂解反應(yīng)已經(jīng)放緩,主要發(fā)生固相縮聚反應(yīng)及石墨結(jié)構(gòu)重排,伴隨著纖維中凝聚態(tài)高次結(jié)構(gòu)的發(fā)展,這些物理結(jié)構(gòu)會(huì)對(duì)微孔產(chǎn)生力的作用,在外力的作用下大孔被擠壓分裂為小孔,孔的總體積也在減小,同時(shí)其邊緣形態(tài)也在復(fù)雜化,在這些因素的綜合作用下,結(jié)構(gòu)參數(shù)發(fā)生如上所述的變化。第Ⅲ階段(1800~2400℃),為石墨化階段。在此階段相對(duì)孔體積Vrel減小,均方回轉(zhuǎn)半徑Rg、比表面積S/V增大,孔隙率P基本保持不變。由于纖維內(nèi)部化學(xué)反應(yīng)基本結(jié)束,并且由于物理結(jié)構(gòu)重排,纖維內(nèi)部結(jié)構(gòu)進(jìn)一步緊密,使得纖維直徑持續(xù)減小,纖維內(nèi)部氣體繼續(xù)溢出,導(dǎo)致纖維中孔總體積減小,所以孔隙率基本不變。而石墨化階段,石墨片層堆砌進(jìn)一步致密,基體內(nèi)部重排產(chǎn)生的應(yīng)力對(duì)微孔結(jié)構(gòu)的壓縮作用也愈發(fā)強(qiáng)烈,導(dǎo)致孔的表面結(jié)構(gòu)進(jìn)一步復(fù)雜化,所以均方回轉(zhuǎn)半徑Rg、比表面積S/V增大。

        圖6 微孔結(jié)構(gòu)參數(shù)隨溫度的變化Fig.6 The change of micro-voids structure parameters during heat treatment process

        根據(jù)以上所得到的結(jié)果,可將PAN基碳纖維微孔缺陷結(jié)構(gòu)隨溫度的變化模型以圖7表示。隨熱處理溫度的提高,纖維中微孔結(jié)構(gòu)變化主要經(jīng)歷從體積變化為主到表面變化為主的過(guò)程,這些過(guò)程的過(guò)渡溫度分別為700,1800℃和2400℃。

        圖7 碳纖維中微孔結(jié)構(gòu)隨溫度的變化模型示意圖Fig.7 The module of the microvoids structure changes during heat treatment process of the PAN based carbon fiber

        3 結(jié)論

        (1)PAN基碳纖維SAXS散射花樣在赤道方向呈扇形分布,將其進(jìn)行全積分處理后,其一維散射曲線為光滑曲線,說(shuō)明其是由微孔結(jié)構(gòu)導(dǎo)致的。

        (2)通過(guò)對(duì)不同熱處理溫度下的一維SAXS散射曲線進(jìn)行理論計(jì)算,得到均方回轉(zhuǎn)半徑Rg、相對(duì)孔體積Vrel、孔隙率P、比表面積S/V四個(gè)孔結(jié)構(gòu)參數(shù)隨溫度的變化,通過(guò)其分析變化規(guī)律,可將PAN基碳纖維微孔結(jié)構(gòu)隨熱處理溫度的變化分為三個(gè)階段:400~700℃以體積變化(孔的融并)為主,表面變化為輔;700~1800℃體積變化(孔的分裂與縮小)與表面變化(邊緣復(fù)雜化)同步;1800~2400℃以分形變化為主,表面變化為輔。

        [1] 陳曉,查劉生,朱育平.高強(qiáng)度、高模量聚丙烯腈基碳纖維的微晶取向研究[J].化工新型材料, 2010,(7): 85-88.

        CHEN Xiao, ZHA Liu-sheng, ZHU Yu-ping. Microstructure of high tensile strength and high tensile modulus PAN-based carbon fibers determined by X-ray diffraction [J]. New Chemical Materials, 2010,(7): 85-88.

        [2] 井敏,譚婷婷,王成國(guó),等.PAN基碳纖維的微觀結(jié)構(gòu)與力學(xué)性能相關(guān)性分析[J].航空材料學(xué)報(bào),2013, 33(1): 78-85.

        JING Min, TAN Ting-ting, WANG Cheng-guo, et al. Correlation between microstructure and mechanical properties of PAN-based carbon fiber [J]. Journal of Aeronautical Materials, 2013, 33(1): 78-85.

        [3] 代志雙,石峰暉,張寶艷,等.碳纖維表面物理特性的表征及其對(duì)浸潤(rùn)性能的影響[J].材料工程,2012,(3): 56-60.

        DAI Zhi-shuang, SHI Feng-hui, ZHANG Bao-yan, et al. Study on the surface physical characteristics of carbon fibers and the effect on wettability [J]. Journal of Materials Engineering, 2012,(3): 56-60.

        [4] 賀福.碳纖維及其應(yīng)用技術(shù)[M].北京:化學(xué)工業(yè)出版社.2004.

        [5] THüNEMANN A F, RULAND W. Microvoids in polyacrylonitrile fibers: a small-angle X-ray scattering study [J]. Macromolecules, 2000, 33(5): 1848-1852.

        [6] SUGIMOTO Y, SHIOYA M, YAMAMOTO K, et al. Relationship between axial compression strength and longitudinal microvoid size for PAN-based carbon fibers [J]. Carbon, 2012, 50 (8):2860-2869.

        [7] LOZANO C D, MACIA J A, CAZORLA A D, et al. Isotropic and anisotropic microporosity development upon chemical activation of carbon fibers, revealed by microbeam small-angle X-ray scattering [J]. Carbon, 2006, 44 (7):1121-1129.

        [8] KABURAGI M, BIN Y, ZHU D, et al. Small angle X-ray scattering from voids within fibers during the stabilization and carbonization stages [J]. Carbon, 2003, 41 (5):915-926.

        [9] LOZANO C D, RAYMUNDO P E, CAZORLA A D, et al. Characterization of pore distribution in activated carbon fibers by microbeam small angle X-ray scattering [J]. Carbon, 2002, 40 (14):2727-2735.

        [10] 盛毅, 張彩紅, 徐耀, 等. 二維小角X射線散射法研究PAN基碳纖維內(nèi)部微孔結(jié)構(gòu)[J]. 新型炭材料, 2009, (3):270-276.

        SHENG Yi, ZHANG Cai-hong, XU Yao, et al. Investigation of PAN-based carbon fiber microstructure by 2D-SAXS [J]. New Carbon Materials, 2009, (3):270-276.

        [11] ZHU C Z, LIU X F, YU X L, et al. A small-angle X-ray scattering study and molecular dynamics simulation of microvoid evolution during the tensile deformation of carbon fibers [J]. Carbon, 2012, 50 (1):235-243.

        [12] ZHU C, YU X, LIU X, et al. 2D SAXS/WAXD analysis of pan carbon fiber microstructure in organic/inorganic transformation-Springer[J]. Chinese Journal of Polymer Science, 2013, 31(5): 823-832.

        [13] 李登華, 吳剛平, 呂春祥, 等. 聚丙烯腈基碳纖維及其原絲中的微孔尺寸分布[J]. 分析測(cè)試學(xué)報(bào), 2010, 29(4):321-326.

        LI Deng-hua, WU Gang-ping, Lü Chun-xiang, et al. Pore size distribution of microvoids in p an-based carbon fibers and the precursors[J]. Journal of Instrumental Analysis, 2010, 29(4):321-326.

        [14] GLATTER O, KRATKY O.Small Angle X-ray Scattering[M]. London: Academic Press, 1982.

        [15] GUINIER A, FOURNET G, YUDOWITCH K L. Small-angle scattering of X-rays [M]. New York: Wiley, 1955.

        [16] 莫志深, 張宏放,張吉東. 晶態(tài)聚合物結(jié)構(gòu)和 X 射線衍射[M]. 北京:科學(xué)出版社, 2010.

        Effect of Temperature on the Micro-void Structure for PAN Fiber During Carbonization

        WANG Meng-fan1,CHEN Wang2,SU Shi-zhou2, HAN Xu1,XU Liang-hua2,CAO Wei-yu1

        (1 State Key Laboratory of Organic-inorganic Composites,Beijing 2 Key Laboratory of Carbon Fiber and Functional Polymer (Ministry of Education),Beijing University of Chemical Technology,Beijing 100029,China)

        Small angle X-ray scattering (SAXS) method was applied to investigate the evolution of microvoids structure during carbonization and graphitization for polyacrylonitrile (PAN) fiber. The SAXS patterns of the PAN based carbon fiber treated by integration are both smooth curves.The scattering signals come from the microvoids structure within fiber. The parameters of the micro-voids structure including mean square radiuses of gyration, relative volume of the microvoids, porosity and specific surface area were calculated. The results indicate that the evolution of the micro-voids structure during carbonization and graphitization can be divided into three stages: 400-700℃, the volume change (mergence of the small microvoids) is faster than the shape change; 700-1800℃, the volume change (split and shrinking of the microvoids) keeps pace with the shape change (complicating of the surface structure); 1800-2400℃, the volume change is slower than the shape change.

        PAN based carbon fiber;micro-void defect structure;small angle X-ray scattering (SAXS);temperature effect

        10.11868/j.issn.1001-4381.2015.06.011

        O631.1

        A

        1001-4381(2015)06-0066-05

        973計(jì)劃項(xiàng)目(2011CB605603)

        2013-12-02;

        2014-02-28

        曹維宇(1972-),男,博士,副研究員,主要從事高性能纖維及高分子物理化學(xué)研究, 聯(lián)系地址:北京市朝陽(yáng)區(qū)北三環(huán)東路15號(hào)北京化工大學(xué)34號(hào)信箱(100029),E-mail:caowy@mail.buct.edu.cn

        猜你喜歡
        回轉(zhuǎn)半徑炭化微孔
        最優(yōu)回轉(zhuǎn)半徑驅(qū)動(dòng)輪總成的設(shè)計(jì)方法研究
        氨基酸帶電量及分布對(duì)蛋白質(zhì)單鏈自組裝行為的影響
        水稻秸稈制備生物制活性碳
        市政污泥炭化時(shí)間與溫度研究
        安徽建筑(2018年4期)2019-01-29 02:16:14
        強(qiáng)疏水性PDMS/PVDF微孔膜的制備及其性能研究
        膜蒸餾用PDMS/PVDF/PTFE三元共混微孔膜制備
        微孔發(fā)泡塑料中成核劑的研究
        炭化米糠經(jīng)臭氧活化制備活性炭及其去除Cr(VI)離子
        新型炭材料(2015年3期)2015-01-01 08:20:20
        芯層微孔發(fā)泡管材的制備
        真空環(huán)保秸稈炭化罐
        精品国产亚洲av久一区二区三区| 亚洲久热无码av中文字幕| 无码少妇一区二区性色av| 国产天堂网站麻豆| 日本岛国大片不卡人妻| 野花视频在线观看免费| 蜜臀av无码人妻精品| 精品国产午夜福利在线观看| 久久久精品国产视频在线| 久久国产精品一区av瑜伽| 日本最新免费二区| 巨大欧美黑人xxxxbbbb| 国产后入内射在线观看| 91精品国产综合久久久密臀九色 | a毛片全部免费播放| 精品国产成人一区二区不卡在线| 亚洲色图专区在线视频| 国产内射爽爽大片| 亚洲av无码之国产精品网址蜜芽| 国产精品原创巨作av无遮| 性色av一区二区三区四区久久 | 国产精品久久久久aaaa| 99福利在线| 亚洲一区二区三区新视频| 日韩av无码社区一区二区三区 | 亚洲a∨无码一区二区| 99久久久精品国产性黑人| 国产实拍日韩精品av在线| 亚洲av日韩av天堂一区二区三区| 亚洲AⅤ精品一区二区三区| 精品女人一区二区三区| 天天躁日日躁狠狠躁欧美老妇小说| 3344永久在线观看视频| 99久久精品国产亚洲av天| 蜜桃tv在线免费观看| ā片在线观看免费观看| 成人日韩av不卡在线观看| 一级一片内射视频网址| 国产喷水1区2区3区咪咪爱av | 91久久精品人妻一区二区| 中文字幕人妻一区二区二区|