999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

原位橢圓偏振光譜研究多孔陽極氧化鋁初期生長過程

2010-11-06 07:01:26雷驚雷張李娜李凌杰謝昭明張勝濤
物理化學學報 2010年9期
關(guān)鍵詞:實驗模型

雷驚雷 張李娜 李凌杰 鄭 莎 王 超 謝昭明 張勝濤

(重慶大學化學化工學院,重慶 400044)

原位橢圓偏振光譜研究多孔陽極氧化鋁初期生長過程

雷驚雷 張李娜?李凌杰*鄭 莎 王 超 謝昭明 張勝濤

(重慶大學化學化工學院,重慶 400044)

根據(jù)鋁陽極氧化初期不同階段Al-H2SO4體系的界面特征建立多個光學模型,并采用有效介質(zhì)近似解析了高時間分辨率的原位橢圓偏振光譜數(shù)據(jù),詳細獲得鋁陽極氧化初期的Al2O3-Al界面層、Al2O3阻擋層以及多孔層組成、厚度變化的動態(tài)信息.根據(jù)光學模型及計算結(jié)果,可以清楚地分辨鋁陽極氧化過程中的阻擋層形成、孔洞萌生、孔洞發(fā)展以及多孔層穩(wěn)定生長等4個階段.在多孔層穩(wěn)定生長階段,阻擋層厚度及多孔層孔隙率幾乎不變,而多孔層厚度隨時間線性增長,其速率為5.8 nm·s-1.

鋁;陽極氧化;模板;橢圓偏振光譜;有效介質(zhì)近似

多孔陽極氧化鋁膜由于具有整齊規(guī)則的納米結(jié)構(gòu),是制備低維光學、電學或磁學納米功能材料的重要模板,近年來隨著納米技術(shù)的迅速發(fā)展而受到廣泛關(guān)注[1].然而,目前的相關(guān)研究大多集中于多孔陽極氧化鋁膜的制備及應用[2-10],探討其多孔結(jié)構(gòu)的生長機制及生長動力學的理論研究相對較少[11-12].這一方面是由于鋁陽極氧化過程涉及復雜的化學、電學、力學等多種因素的耦合而導致理論研究困難,另一方面也是由于鋁化學性質(zhì)活潑而致使很難應用常規(guī)研究方法獲取反應中間過程動態(tài)信息并推測其反應機理.

橢圓偏振光譜法作為一種高靈敏度、非破壞性的原位光譜電化學研究方法,可在進行電化學研究的同時獲取電極表面膜層厚度、組成等信息,為推測反應機理提供有力證據(jù)[13-14].將橢圓偏振光譜法應用于鋁陽極氧化過程的研究已有少量報道[15-18],然而這些研究均是單波長測量或非原位測量,且時間分辨率較低,不能準確、及時地提供反應過程中,特別是鋁陽極氧化初期的電極表面膜層的實時動態(tài)變化信息.由于鋁陽極氧化初期涉及阻擋層和多孔層的生成、發(fā)展等階段[19-20],對于探明鋁多孔陽極氧化膜的生長機制非常重要,已引起人們的關(guān)注[21].隨著儀器制造技術(shù)和計算機技術(shù)的發(fā)展,單次采集橢圓偏振光譜數(shù)據(jù)時間縮短到數(shù)十毫秒,時間分辨率大大提高[22].本文通過解析多孔陽極氧化鋁膜生長初期的高時間分辨率原位橢圓偏振光譜數(shù)據(jù),獲取膜層厚度、組成隨時間的變化信息.

1 實驗

實驗采用兩電極體系,工作電極為環(huán)氧樹脂封裝的純度為99.9%的金屬Al(西南鋁業(yè)集團有限公司),用400#、800#、1200#和2000#水砂紙逐級打磨至鏡面光滑后進行電化學拋光(在體積比為1∶4的HClO4與 C2H5OH溶液中以 15.0 V電壓拋光 1 min);對電極為大Pt片;電解液為1.1 mol·L-1的H2SO4(重慶川東化工有限公司化學試劑廠)溶液;采用直流穩(wěn)壓電源(北京大華電子集團DH1716A-13)控制電壓為16.0 V;實驗溫度為23.0℃.橢圓偏振光譜(美國J.A.Woollam公司M2000U)測量的波長λ范圍為600-800 nm,入射角為70°;儀器單次采集光譜時間需38 ms,為提高信噪比,以40次采集的平均值作為輸出信號,即光譜的實際時間分辨率為1.5 s(在實驗總時間90 s內(nèi),共記錄63次光譜數(shù)據(jù));橢圓偏振數(shù)據(jù)解析采用儀器自帶軟件WVASE 32. Al2O3的光學常數(shù)采用軟件提供的數(shù)據(jù)庫[23];在進行陽極氧化之前,采用橢圓偏振光譜法在空氣中測量并擬合Al基底的光學常數(shù).

2 結(jié)果與討論

圖1是氧化電流隨時間的變化曲線,其形狀和文獻報道[20,24]的非常相似,表現(xiàn)出先急劇減小,越過最低點后逐漸回升直至幾乎穩(wěn)定的特征.其中插圖為前28 s的電流-時間曲線放大圖,可以發(fā)現(xiàn)約在2.6 s時,電流出現(xiàn)最小值.

記錄電流的同時記錄橢圓偏振光譜數(shù)據(jù).圖2分別繪出波長為600和800 nm的橢圓偏振光譜參量△-ψ實驗曲線.借鑒解析橢圓偏振光譜的“環(huán)境-膜層-基底”經(jīng)典模型[13],可以建立如圖3所示的“H2SO4溶液-Al2O3層-Al基底”3層膜光學模型.由于Al2O3和H2SO4的消光系數(shù)均為0,根據(jù)此模型計算,隨Al2O3層厚度增加,△-ψ理論曲線應為一條封閉曲線[25-26].但從圖2可以看出,△-ψ實驗曲線并不封閉,表明進行實驗的90 s內(nèi),多孔層已經(jīng)生成并顯著影響橢圓偏振光譜測量結(jié)果.逐一對各時刻的ψ-λ和△-λ曲線進行擬合,可以得知膜層厚度、組成隨時間的變化信息.

圖1 Al在1.1 mol·L-1H2SO4溶液中的陽極氧化電流-時間曲線Fig.1 Anodic current transient of Al in 1.1 mol·L-1H2SO4solutionThe inset shows the current from 0 to 28 s.The letters A,B,C,and D correspond to the time of the first 4 ellipsometric measurements,i.e., 0.6,2.2,3.5,5.0 s,respectively.The letter E corresponds to the time of porous layer starts to grow at a stable rate(26.2 s),which is determined from the ellipsometric measurements.

圖2 波長λ分別為600和800 nm的橢圓偏振光譜參量△-ψ實驗曲線Fig.2 Experimental curves of△-ψ at the wavelength of 600 and 800 nm

圖3 解析阻擋層形成階段的橢圓偏振光譜實驗數(shù)據(jù)的“H2SO4溶液-Al2O3層-Al基底”3層膜模型Fig.3 Three-layer model of“H2SO4solution-Al2O3-Al substrate”used to fit the ellipsometric data(a)physical model;(b)optical model

采用圖3所示的3層膜光學模型(僅擬合Al2O3層厚度)可以很好地擬合對應于通電后0.6 s(圖1中標記為A點)的第1次原位橢圓偏振光譜數(shù)據(jù),擬合得到的理論曲線與實驗曲線符合很好(如圖4所示).此時,在Al陽極氧化初期,將按照如下反應生成Al2O3[20]:

圖4 在0.6 s時的橢圓偏振光譜實驗曲線和由3層膜模型計算得到的理論曲線的比較Fig.4 Comparison of the experimental ellipsometric spectra at 0.6 s with the theoretic spectra generated with the three-layer optical model

Al2O3形成致密的阻擋層,由于其電阻很大,隨阻擋層厚度增加,通過體系的電流迅速減小,此即文獻中通常認為的“阻擋層形成”階段[19-20].采用同樣模型,也可以很好地擬合對應于通電后2.2 s(圖1中標記為B點)的第2次原位橢圓偏振光譜數(shù)據(jù).比較發(fā)現(xiàn),阻擋層厚度隨時間的延長而增大,這與電流隨時間減小的實驗事實相符.

隨著反應的進行,Al2O3層靠近溶液一側(cè)部分開始發(fā)生如下的化學溶解反應[17]:

溶解將形成孔洞,此即“孔洞萌生”階段[19-20].采用如圖5所示的“溶液-多孔層-Al2O3層-Al基底”4層膜光學模型可以比前述3層膜模型更好地擬合對應于通電后3.5 s(圖1中標記為C點)的第3次原位橢圓偏振光譜數(shù)據(jù),其中多孔層可看作由Al2O3和溶液按照一定比例組成,其等效復折射率n可以按照有效介質(zhì)近似(EMA)計算而得[27-28]:

式中nsolution和nAl2O3分別表示H2SO4溶液和Al2O3的復折射率,φsolution表示溶液在多孔層中的體積分數(shù).由于溶液填充孔洞,因此可以認為溶液的體積分數(shù)等于多孔層中孔洞的體積分數(shù),即孔隙率,它與孔洞的密度、直徑有關(guān).值得說明的是,盡管在上述4層膜光學模型中,我們同時擬合Al2O3層厚度(dAl2O3)、多孔層的厚度(dPL)及其中溶液的體積分數(shù)(φsolution),但并不能很好地擬合對應于圖1中A、B點的原位橢圓偏振光譜數(shù)據(jù),意味著在B點對應時刻(2.2 s)前還沒有孔洞萌生.由于剛萌生的孔洞非常細小,很難應用常規(guī)方法進行研究,通常是和僅生成阻擋型氧化鋁膜的電流-時間曲線相比較以確定孔洞的萌生時間:當僅生成阻擋型氧化鋁膜時,電流單調(diào)減小;而生成多孔型氧化鋁膜時電流會出現(xiàn)先減小再增大最后逐漸穩(wěn)定的過程,故往往將電流出現(xiàn)差別的時刻作為孔洞萌生的開始[26],也有文獻[29]認為,孔洞萌生于電流-時間曲線上最低點對應時間的1/2.然而,由于阻擋型和多孔型膜在不同電解質(zhì)、電壓等條件下生成,因此上述方法并不準確.我們高靈敏度的原位橢圓偏振光譜實驗為辨別孔洞萌生提供了實驗依據(jù),結(jié)果表明孔洞的產(chǎn)生非常接近于電流-時間曲線上的最低點.

圖5 解析孔洞萌生階段橢圓偏振光譜實驗數(shù)據(jù)的“H2SO4溶液-多孔層-Al2O3層-Al基底”4層膜模型Fig.5 Four-layer model of“H2SO4solution-porous layer-Al2O3-Al substrate”used to fit the ellipsometric data(a)physical model;(b)optical model

圖6 解析孔洞生長階段橢圓偏振光譜實驗數(shù)據(jù)的“H2SO4溶液-多孔層-Al2O3層-界面層-Al基底”5層膜模型Fig.6 Five-layer model of“H2SO4solution-porous layer-Al2O3-interphase layer-Al substrate”used to fit the ellipsometric data(a)physical model;(b)optical model

然而,前述的4層膜模型不能擬合對應于通電后5.0 s(圖1中標記為D點)的第4次原位橢圓偏振光譜數(shù)據(jù).考慮到隨反應進行,Al2O3阻擋層和Al基底界面逐漸呈球弧狀而出現(xiàn)由Al2O3阻擋層和Al基底組成的“界面層”,光學模型演變?yōu)槿鐖D6所示的“溶液-多孔層-Al2O3層-界面層-Al基底”5層膜模型[15-16],其中界面層的等效光學常數(shù)可以由Al2O3和Al的光學常數(shù)按有效介質(zhì)近似計算而得.用此5層膜模型,同時擬合多孔層的厚度及其中溶液的體積分數(shù)、Al2O3層厚度、界面層厚度(dIL)及其中Al2O3的體積分數(shù),可以很好地解析剩余所有60次橢圓偏振光譜數(shù)據(jù).

圖7 原位橢圓光譜擬合結(jié)果Fig.7 Results of deconvoluting of the in situ ellipsometric spectraI,II,III,and IV correspond to the stages of the barrier layer grows, pores form,pores enlarge,and pores grow at a stable rate,respectively. The lines B,C,and E correspond to the time of 2.2,3.5,and 26.2 s, which are marked as the same letters in Fig.1. φAl2O3:volume fraction of Al2O3;φsolution:volume fraction of solution; dIL:interface layer thickness;dPL:porous layer thickness

所有擬合結(jié)果如圖7(a-e)所示.圖中字母B、C、E對應于圖1中相應字母表示的時刻.由圖7(a,b)可以發(fā)現(xiàn),26.2 s(圖1中標記為E點)之后,多孔層孔隙率φsolution(等同于溶液的體積分數(shù))幾乎不變,而僅僅厚度dPL隨時間改變,表明此時體系進入“多孔層穩(wěn)定生長”階段.文獻中通常以電流達到穩(wěn)定(本實驗中21.0 s)作為此階段的標志[29],但原位橢圓偏振光譜實驗結(jié)果表明多孔層穩(wěn)定生長應稍晚于電流穩(wěn)定.在孔洞萌生與多孔層穩(wěn)定生長階段之間,則是“孔洞發(fā)展”階段.此階段由原位橢圓偏振光譜計算得到的各膜層厚度、組成改變的規(guī)律不明顯,這可能與舊孔洞發(fā)展的同時有新孔洞萌生(即孔洞萌生與發(fā)展階段重疊)[20]以及孔洞同時橫向、縱向發(fā)展[30]有關(guān).

綜上所述,根據(jù)解析橢圓偏振光譜實驗數(shù)據(jù)的光學模型及計算結(jié)果,可以清楚地分辨鋁陽極氧化過程中阻擋層形成、孔洞萌生、孔洞發(fā)展以及多孔層穩(wěn)定生長等階段(在圖7中分別以I、II、III、IV表示).文獻中通常采用非原位測量方法對上述4個階段進行研究[16,21,30-31],由于其測量的不連續(xù)性而可能遺漏各階段演變的部分信息.本文采用高時間分辨率、高靈敏度的原位橢圓偏振光譜電化學方法詳細地獲得了各階段的界面層、阻擋層以及多孔層的組成、厚度變化的動態(tài)信息,為各階段的劃分提供了實驗依據(jù).

進一步研究發(fā)現(xiàn),可以用線性關(guān)系擬合“多孔層穩(wěn)定生長”階段多孔層厚度隨時間t的改變,其方程及相關(guān)系數(shù)R2為:

上式表明多孔層的生長速度約為5.8 nm·s-1,由此可估算制備一定厚度的多孔陽極氧化鋁納米模板所需要的時間.此生長速度與文獻報道相近條件(21℃,3 mol·L-1H2SO4溶液中,10mA·cm-2恒電流陽極氧化)下采用非原位方法測定的多孔層生長速度5.9 nm·s-1非常接近[16].另外,通常認為阻擋層厚度dAl2O3與制備電壓的關(guān)系為1.0 nm·V-1[24],對應于本實驗16.0 V電壓的阻擋層厚度為16.0 nm,此數(shù)值與根據(jù)橢圓偏振光譜實驗計算得到的多孔層穩(wěn)定生長階段的阻擋層厚度平均值18.2 nm非常接近.上述結(jié)果表明我們的模型和計算結(jié)果的合理性.

3 結(jié)論

借助于合理的理論模型,原位橢圓偏振光譜法可以很好地用于研究組成、厚度同時改變的多層膜復雜體系,為深入理解電化學反應機制提供更豐富的信息.本文利用高時間分辨率、高靈敏度的原位橢圓偏振光譜法研究了鋁陽極氧化過程,結(jié)果表明: (1)利用恰當?shù)墓鈱W模型以及有效介質(zhì)近似可以很好地解析鋁陽極氧化過程中的原位橢圓偏振光譜實驗數(shù)據(jù),從而獲得界面層、阻擋層以及多孔層的組成、厚度變化的動態(tài)信息;根據(jù)所用光學模型及計算結(jié)果,可清楚地分辨鋁陽極氧化過程中的各個階段; (2)多孔層穩(wěn)定生長階段中,阻擋層厚度及多孔層的孔隙率幾乎不變而多孔層厚度隨時間線性增長,其速度為5.8 nm·s-1.

1 Chik,H.;Xu,J.M.Mater.Sci.Eng.R,2004,434:103

2 Losic,D.;Lillo,M.;Losic,J.D.Small,2009,5:1392

3 Rao,C.N.R.;Govindaraj,A.Adv.Mater.,2009,21:4208

4 Tao,F.F.;Xu,Z.Acta Phys.-Chim.Sin.,2009,25:977 [陶菲菲,徐 正.物理化學學報,2009,25:977]

5 Chen,W.;Wu,J.S.;Xia,X.H.ACS nano,2008,2:959

6 Chen,W.;Wang,F.B.;Ruan,M.H.;Xu,J.J.;Xia,X.H.Sci.Adv. Mater.,2009,1:25

7 Chen,W.;Wu,J.S.;Yuan,J.H.;Xia,X.H.;Lin,X.H. J.Electroanal.Chem.,2007,600:257

8 Lee,W.;Ji,R.;Goesele,U.;Nielsch,K.Nat.Mater.,2006,5:741

9 Chen,S.S.;Ling,Z.Y.;Hu,X.;Yang,H.;Li,Y.J.Mater.Chem., 2010,20:1794

10 Sabzi,R.;Kant,K.;Losic,D.Electrochemica Acta,2010,55: 1829

11 Gong,Y.L.;Wang,W.;Wang,H.;Guo,H.T.Acta Phys.-Chim. Sin.,2004,20:199 [鞏運蘭,王 為,王 惠,郭鶴桐.物理化學學報,2004,20:199]

12 Stanton,L.G.;Golovin,A.A.Phys.Rev.B,2009,79:035414

13 Azzam,R.M.A.;Bashara,N.M.Ellipsometry and polarized light. Amsterdam:North-Holland Publishing Company,1977:283,452

14 Lei,J.L.;Li,L.J.;Cai,S.M.;Zhang,S.T.;Li,D.;Yang,M.Z. Acta Phys.-Chim.Sin.,2001,17:1107 [雷驚雷,李凌杰,蔡生民,張勝濤,李 荻,楊邁之.物理化學學報,2001,17:1107]

15 Laet,J.D.;Vanhellemont,J.;Terryn,H.;Vereccken,J.Appl.Phys. A,1992,54:72

16 Schnyder,B.;K?tz,R.J.Electroanal.Chem.,1992,339:167

17 Laet,J.D.;Terryn,H.;Vereccken,J.Thin Solid Films,1998,320: 241

18 Kooij,E.S.;Wormeester,H.;Ga?lcǎ,A.C.;Poelsema,B. Electrochem.Solid-State Lett.,2003,6:B52

19 Parkhutik,V.P.;Shershulsky,V.I.J.Phys.D-Appl.Phys.,1992, 25:1258

20 Zhu,Z.F.Techniques of anodic oxidizing and surface treating for aluminum alloys.Beijing:Chemical Industry Press,2004:94 [朱祖芳.鋁合金陽極氧化與表面處理技術(shù).北京:化學工業(yè)出版社,2004:94]

21 Wu,Z.;Richter,C.;Menon,L.J.Electrochem.Soc.,2007,154:E8

22 Okabe,H.;Hayakawa,M.;Matoba,J.;Naito,H.;Oka,K.Rev.Sci. Instrum.,2009,80:083104

23 Palik,E.D.Handbook of Optical Constants of Solids(III).San Diego:Academic Press,1998:653

24 Hoar,T.P.;Yahalom,J.J.Electrochem.Soc.,1963,110:614

25 Chen,T.X.Experimental electrochemistry.Xiamen:Xiamen University Press,1993:226 [陳體銜.實驗電化學.廈門:廈門大學出版社,1993:226]

26 Diggle,J.W.;Downie,T.C.;Goulding,C.W.Chem.Rev.,1969, 69:365

27 Aspnes,D.E.;Kinsbron,E.;Bacon,D.D.Phys.Rev.B,1980,21: 3290

28 Fujiwara,H.Spectroscopic dllipsometry,principles and applications.West Sussex:John Willy&Sons Ltd.,2007:177

29 Siejka,J.;Ortega,C.J.Electrochem.Soc.,1977,124:883

30 Xu,Y.;Thompson,G.E.;Wood,G.C.J.Chin.Soc.Corr.Prot., 1989,9:1 [徐 源,Thompson,G.E.,Wood,G.C.中國腐蝕與防護學報,1989,9:1]

31 Thompson,G.E.;Furneaux,R.C.;Wood,G.C.;Richardson,J.A.; Goode,J.S.Nature,1978,272:433

In situ Ellipsometric Study on the Initial Stages of Porous Anodization of Aluminum

LEI Jing-Lei ZHANG Li-Na?LI Ling-Jie*ZHENG Sha WANG Chao XIE Zhao-Ming ZHANG Sheng-Tao
(College of Chemistry and Chemical Engineering,Chongqing University,Chongqing 400044,P.R.China)

The initial stages of the anodization of aluminum are influential during the preparation of anodic aluminum oxide nanotemplates and provide an insight into their formation mechanism.The formation and development of both the barrier layer and the porous layer are involved.In this paper,in situ ellipsometric spectra with a high time-resolution were collected.To deconvolute the ellipsometric spectra,several optical models were extracted from the physical models of the Al-H2SO4interface.Using these models and the effective medium approximation, detailed information about the composition and thickness of the Al2O3-Al interphase,Al2O3barrier layer,and porous layer was acquired.Based on the deconvoluted results and the optical models,4 stages of the porous anodization of aluminum,i.e.,the barrier layer grows,pores form,pores enlarge,and pores grow at a stable rate,were clearly distinguished.Moreover,during the last stage the thickness of the porous layer changes linearly with time at a rate of 5.8 nm·s-1while both the thickness of the barrier layer and the porosity of the porous layer change very little.

Aluminum;Anodization;Template;Spectroscopic ellipsometry;Effective medium approximation

O646.54

Received:February 19,2010;Revised:April 18,2010;Published on Web:July 7,2010.

*Corresponding author.Email:LJLi@cqu.edu.cn.

The project was supported by the National Natural Science Foundation of China(20603049,20803097),Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education,China([2007]1108-4)and Natural Science Foundation Project of CQ CSTC

(2008BB4174).

國家自然科學基金(20603049,20803097)和教育部留學回國人員科研啟動基金(教外司留[2007]1108-4),重慶市自然科學基金

(CSTS2008BB4174)資助項目

?現(xiàn)在地址:四川天華股份有限公司,郵編646207.

?Editorial office of Acta Physico-Chimica Sinica

猜你喜歡
實驗模型
一半模型
記一次有趣的實驗
微型實驗里看“燃燒”
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
做個怪怪長實驗
3D打印中的模型分割與打包
NO與NO2相互轉(zhuǎn)化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
主站蜘蛛池模板: 午夜天堂视频| 在线观看国产黄色| 国产成人8x视频一区二区| 人妻精品久久久无码区色视| 找国产毛片看| 亚洲午夜片| 青青青国产视频| 亚洲视频色图| 亚洲va欧美va国产综合下载| 精品日韩亚洲欧美高清a| 国产网站一区二区三区| 欧美成人二区| 伊人国产无码高清视频| 国产人成在线视频| 亚洲欧美精品日韩欧美| 亚洲欧美精品一中文字幕| 中文字幕亚洲无线码一区女同| 91小视频在线播放| 波多野结衣第一页| 国产亚洲精久久久久久久91| 亚洲精品777| 米奇精品一区二区三区| 91av成人日本不卡三区| jizz亚洲高清在线观看| 无码网站免费观看| 在线观看视频99| 亚洲自偷自拍另类小说| A级毛片高清免费视频就| 国产精品亚洲五月天高清| 国产制服丝袜无码视频| 人妻无码中文字幕一区二区三区| 潮喷在线无码白浆| 亚洲国产成熟视频在线多多| 999精品色在线观看| 97国产在线视频| 亚洲天堂免费| 亚洲丝袜第一页| 国产一级精品毛片基地| 日韩成人在线网站| 国产精品成人一区二区不卡| 日本免费福利视频| 毛片在线区| 不卡无码h在线观看| 高清无码一本到东京热| 亚洲AV无码乱码在线观看裸奔| 亚洲一区精品视频在线| 欧美激情视频在线观看一区| 在线视频97| 国产午夜福利片在线观看| 99这里精品| 久久a级片| 中文国产成人久久精品小说| 欧美日韩综合网| 最新国产成人剧情在线播放| 精品丝袜美腿国产一区| 亚洲全网成人资源在线观看| 中文字幕亚洲无线码一区女同| 国产乱码精品一区二区三区中文 | 日本午夜精品一本在线观看| 在线欧美日韩国产| 精品人妻无码中字系列| 国产在线自乱拍播放| 欧美va亚洲va香蕉在线| 亚洲欧美激情小说另类| A级全黄试看30分钟小视频| 女人毛片a级大学毛片免费| 婷婷六月在线| 久久久久国产精品免费免费不卡| 国产全黄a一级毛片| 国产精品无码制服丝袜| 视频二区欧美| 99re热精品视频国产免费| 亚洲天堂网在线视频| 国产精品综合色区在线观看| 亚洲资源站av无码网址| 日韩毛片免费| 国产精品亚洲专区一区| 亚洲视频色图| 最新国语自产精品视频在| 精品视频福利| 4虎影视国产在线观看精品| 国产97视频在线观看|