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

    艦船空氣凈化材料吸附CO2的分子模擬研究

    2012-11-12 08:03:14李躍華徐筱欣
    中國艦船研究 2012年1期
    關(guān)鍵詞:力場混合氣構(gòu)型

    李躍華 徐筱欣

    上海交通大學(xué) 動(dòng)力裝置及其自動(dòng)化研究所,上海 200240

    1 引 言

    在艦艇上,CO2是一種有害氣體。大氣中的CO2濃度通常為0.03%,艦艇上該氣體的安全濃度通常為0.5%~1.0%[1]。在其濃度高于這一范圍時(shí),艇員就會(huì)感到不適,甚至是難以完成工作。因此,當(dāng)艦艇上的CO2濃度超標(biāo)時(shí),必須采用一定的去除技術(shù)以降低濃度。

    目前,艦艇上采用的空氣凈化裝置對CO2的吸收主要是通過吸收劑來完成。曾有人進(jìn)行過一項(xiàng)統(tǒng)計(jì):美國和英國的核潛艇主要使用一乙醇胺(MEA)吸收CO2;法國海軍核潛艇采用分子篩吸附CO2;東歐和前蘇聯(lián)國家的常規(guī)潛艇一般采用超氧化物(超氧化鉀或超氧化鈉)吸收CO2并釋放O2,也有采用MEA的。

    對于動(dòng)力系統(tǒng)在某一封閉的箱裝體內(nèi)運(yùn)行的情況,由于CO2的成分也比較高,其內(nèi)部氣體在排出到艙室之前必須經(jīng)過凈化處理,以達(dá)到艙室排放要求。孟帥等[2]研制了用于該場合的空氣凈化器,吸附CO2的材料主要也是分子篩。

    沸石是一種重要的分子篩,是結(jié)晶態(tài)的硅酸鹽或硅鋁酸鹽,由硅氧四面體或鋁氧四面體通過氧橋鍵相連而形成分子尺寸大?。ㄍǔ?.3~2.0 nm)的孔道和空腔體系,從而具有篩分分子的特性。

    生產(chǎn)廠商提供的各種分子篩對CO2的吸附量通常是在25℃下測得。而對于艦船,無論是用于凈化艙室空氣還是凈化箱裝體之類的密閉環(huán)境,其壓力約為1個(gè)大氣壓,但溫度有可能會(huì)不同。因此,有必要研究分子篩在1個(gè)大氣壓左右、不同溫度下的CO2吸附性能。另外,孟帥等人開發(fā)的空氣凈化裝置在吸附CO2時(shí),往往還要面臨空氣中大量的 N2、O2以及少量 CO、CxHy等的競爭。因此,研究其競爭吸附性能同樣非常重要。

    吸附CO2最常用的分子篩為13X或13XAPG分子篩,后者為應(yīng)用于空分領(lǐng)域的專門分子篩,但二者的結(jié)構(gòu)分子式一樣。本文將著重分析13X-APG分子篩在不同溫度下的CO2吸附性能和競爭吸附性能。

    分析采用的方法是計(jì)算機(jī)分子模擬下的蒙特卡洛方法[3-4]。所謂分子模擬[5],是指利用量子力學(xué)或經(jīng)典牛頓力學(xué)的知識,從微觀層面模擬某些分子或粒子,研究其微觀活動(dòng),進(jìn)而體現(xiàn)其宏觀性質(zhì)。量子力學(xué)法適用于簡單的分子或電子數(shù)目較少的體系,在實(shí)際運(yùn)用中,一般使用經(jīng)典力學(xué)方法,借助于分子的立場來模擬。本文應(yīng)用的蒙特卡洛方法就是在經(jīng)典的分子力學(xué)基礎(chǔ)上,結(jié)合統(tǒng)計(jì)力學(xué)中概率分配原理后得出。

    模擬中采用的軟件為美國Accelrys公司的Materials Studio(MS)軟件。

    2 13X-APG分子篩建模

    13X-APG分子篩屬于NaX分子篩,Si為+4價(jià),硅氧四面體沒有骨架外陽離子,Al為+3價(jià),所以,1個(gè)鋁氧四面體外得有1個(gè)骨架外陽離子Na+來平衡。根據(jù)文獻(xiàn)[6],13X-APG分子篩的模型分子式為Na86Al86Si106O384·265H2O。去掉結(jié)晶水后,在MS中建立含有86個(gè)AlO2,106個(gè)SiO2的晶胞,Al和Si骨架外的86個(gè)Na使用MS自帶的Sorption模塊下的Locate子模塊吸附,并用Forcite模塊優(yōu)化幾何模型。最后得到13X-APG分子篩的一個(gè)晶胞,晶胞參數(shù)為:a=b=c=25.028 ?,α=β=γ=90°(a、b、c 為晶胞邊長,90°說明晶胞為四方體)。

    由于分子篩骨架原子數(shù)目較大,需要從模型中切取3簇小模型,使用MS中的Dmol3模塊計(jì)算 3簇模型的靜電勢 (Electro-static Potential,ESP)電荷值[7]。最后分別取三者的平均值,得到的計(jì)算結(jié)果如表1所示。

    表1 13X-APG分子篩原子電荷Tab.1 Atomic charge in 13X-APG zeolite

    CO2模型比較常用,根據(jù)文獻(xiàn)[8],CO2的模型參數(shù)如表2所示。

    表2 CO2的模型參數(shù)Tab.2 Parameters of Molecular Model of CO2

    3 巨正則系綜蒙特卡洛模擬

    蒙特卡洛(Monte Carlo,MC)算法是一種統(tǒng)計(jì)意義上的算法,具有較大的隨機(jī)性。在重要抽樣的分子模擬中,MC通過改變當(dāng)前體系里粒子的位置、朝向和構(gòu)象,例如旋轉(zhuǎn)、平移和交換粒子,得到粒子的新構(gòu)型。新的構(gòu)型會(huì)對應(yīng)粒子新的能量,通過比較新舊構(gòu)型能量的變化情況來決定這個(gè)新的構(gòu)型是否被接受,其是以概率1完全接受,還是以某一小于1的概率接受。重復(fù)這樣的步驟,收集每一個(gè)被接受的構(gòu)型。對于每一個(gè)被接受的構(gòu)型,計(jì)算其某些物理量(如總能量、均方根位移等)的平均值,如果這些物理量的平均值最后是收斂的,就說明產(chǎn)生了可靠的平衡系統(tǒng)。

    在本文的模擬中,假定CO2為剛性,分子篩和CO2之間以及CO2彼此之間通過原子間非鍵能發(fā)生作用,也就是范德華力和靜電力作用下的勢能,分別用 Lennard-Jones(LJ)勢和 Coulombic勢來描述。二者組合起來便是:

    式中,“+”前為LJ勢(范德華力產(chǎn)生的勢能),“+”后為Coulombic勢(靜電力產(chǎn)生的勢能);下標(biāo)i、j均代表2個(gè)粒子,它們可能相同,也可能不同;rij為粒子之間的距離;ε0為真空介電常數(shù),具體數(shù)值為ε0=8.8542 ×10-12C2N-1m-2;εij和 σij為粒子之間相互作用的能量參數(shù)和尺寸參數(shù) (它們也是軟件進(jìn)行分子模擬時(shí)力場的主要參數(shù)),特別是在i和j不相同時(shí),ε和σ往往要采用某種方式來估計(jì),2 種常用的估算公式如式(2)、式(3)所示。

    另外,值得一提的是,本文使用的MS模擬軟件的力場參數(shù)通常以σ和ε/kB表示。其中kB為玻爾茲曼常數(shù),單位為J/K;ε/kB的單位為K。

    13X-APG分子篩采用MS自帶的UFF和Dreiding力場,UFF不同原子間σ和ε的混合規(guī)則為:

    Dreiding力場的混合規(guī)則為:

    本文介紹的模擬是以軟件給定的概率旋轉(zhuǎn)、平移和交換吸附質(zhì)分子(CO2、N2、O2),然后再與原構(gòu)型相比,計(jì)算二者勢能的差值Δu,即利用式(1)計(jì)算新舊構(gòu)型的勢能差。若Δu≤0,則完全接受這個(gè)新構(gòu)型;若Δu>0,則并不是一味地拒絕,而是讓系統(tǒng)另外產(chǎn)生一個(gè)小于1的隨機(jī)數(shù)p,用玻爾茲曼因子 exp(-Δu/kBT)與這個(gè) p 比較,若其大于等于 p,則這個(gè)新構(gòu)型以概率 exp(-Δu/kBT)被接受,若小于p,則拒絕,構(gòu)型不做改變(其中T為系統(tǒng)溫度)。如此往復(fù)計(jì)算300萬步,直至能量值收斂、系統(tǒng)達(dá)到平衡。其中,前150萬步用于使系統(tǒng)達(dá)到平衡,后150萬步用于取點(diǎn),統(tǒng)計(jì)各種熱力學(xué)物理量。計(jì)算Δu時(shí),會(huì)涉及到計(jì)算式(1)中加號后的靜電作用,模擬中采用Ewald于1921年提出的Ewald求和法進(jìn)行,其截?cái)喟霃綖?2.5×10-10m,計(jì)算精度為4.186 J/mol,“+”前面的 LJ勢主要是力場參數(shù)的改變。

    4 結(jié)果與討論

    4.1 驗(yàn)證力場

    本文CO2采用的力場為UFF力場,具體的力場參數(shù)如表3所示。分子篩采用UFF和Dreiding力場時(shí),最后的CO2吸附量是不一樣的,所以有必要事先驗(yàn)證合適的分子篩力場,其力場參數(shù)如表4所示。

    根據(jù)這些參數(shù),在25℃(298 K)下,分別做出分子篩使用UFF力場和Dreiding力場時(shí)其CO2吸附量隨壓力變化的情況,即其等溫吸附線,如圖1所示。 圖中包括了在 33.320 kPa、60 kPa、101 kPa和298 K下13X-APG分子篩吸附CO2的實(shí)測值,具體數(shù)據(jù)如表5所示(數(shù)據(jù)來源于環(huán)球分子篩有限公司)。

    表3 CO2的力場參數(shù)Tab.3 Force field parameters of CO2

    表4 分子篩所使用的力場參數(shù)Tab.4 Force field parameters of zeolite

    表5 13X-APG分子篩吸附CO2的實(shí)際值Tab.5 Actual CO2value absorbed by 13X-APG zeolite

    不難看出,分子篩采用Dreiding力場時(shí),模擬的CO2吸附量要比實(shí)際值小很多,而采用UFF力場則基本吻合。因此,在后面的模擬當(dāng)中,選擇的分子篩力場為UFF力場(表4),選擇的CO2力場也為UFF(表 3)。

    4.2 空氣中CO2的吸附模擬

    由于空氣凈化器使用環(huán)境不同,從陸上試驗(yàn)到裝船實(shí)用,一般溫度不同,因此溫度對吸附量的影響還有待確認(rèn)。另外,空氣中N2和O2的含量約占99%,CO2只是極小一部分,吸附的時(shí)候,CO2還要面臨N2和O2的競爭。本節(jié)將通過使用上節(jié)驗(yàn)證過的力場,研究分子篩在20~40℃范圍內(nèi)與N2、O2的競爭吸附性能。

    事實(shí)上,艦船密閉環(huán)境的空氣中除以上3種氣體外,還有CO以及柴油揮發(fā)出的烴類氣體,但這些氣體所占的比重與CO2相比都極小,因此模擬時(shí)忽略了。根據(jù)文獻(xiàn)[9],CO2的濃度通常在1.2%~1.4%之間,本次模擬以 1.4%計(jì)算。 另外,考慮到密閉環(huán)境中機(jī)組運(yùn)行時(shí)需要很小的真空度,本次模擬的總壓以100 kPa計(jì)算。表6所示為艦船密閉環(huán)境中各氣體的組成。

    表6 艦船密閉環(huán)境中氣體的組成Tab.6 Air composition in warship closed environment

    模擬時(shí),N2和 O2采用聯(lián)合原子 TraPPE力場[10],用球形粒子表示,其LJ參數(shù)如表7所示。

    表7 N2和O2的力場參數(shù)Tab.7 Force field parameters of N2and O2

    表8所示為N2、O2、CO2在各自分壓下,以混合空氣形式在13X-APG分子篩中的競爭吸附情況,它用2種單位顯示氣體的吸附量:一種是每個(gè)分子篩晶胞吸附的氣體分子數(shù) (也是模擬軟件直接顯示的單位);另一種是每克分子篩吸附的氣體毫摩爾數(shù)。表9所示為在各自分壓下,N2、O2、CO2在13X-APG分子篩中的單組份吸附情況,同樣也用2種單位表示氣體的吸附量。

    表8 混合氣體在分子篩中的競爭吸附量Tab.8 Competition adsorption amount of gas mixture in zeolite

    表9 單組份氣體在分子篩中的吸附量Tab.9 Adsorption amount of pure gas in zeolite

    圖2所示為混合氣中CO2和單組份CO2在分子篩中的吸附情況。圖3所示為混合氣中N2、O2和其單組份時(shí)的吸附情況。

    由表8可見,混合氣在經(jīng)過分子篩后,與CO2相比,N2與O2的吸附量只是極小的一部分,幾乎可以認(rèn)為N2和O2不被分子篩吸附。這也為分子篩作為艦船空氣凈化器的吸附材料提供了可能。

    結(jié)合表8、表9和圖2、圖3還可發(fā)現(xiàn),不管是CO2還是N2、O2,不管是混合氣體還是單組份氣體,它們的吸附量隨著溫度的升高都在減小,這是因?yàn)闅怏w在分子篩上的吸附為物理吸附,溫度升高吸附量便會(huì)下降。具體到混合氣CO2,隨著溫度的升高,其吸附量從20℃時(shí)的2.894 mmol/g下降到了 40℃時(shí)的 2.465 mmol/g,也就是從 127.3 mg/g下降到了108.5 mg/g,降幅達(dá)14.8%。應(yīng)該說這個(gè)降幅還是非常大的,對于在常溫下做實(shí)驗(yàn)的樣機(jī),這樣的降幅必須予以考慮。

    另外,混合氣中CO2的分壓是1.4 kPa,單獨(dú)測1.4 kPa下單組份CO2的吸附情況時(shí) (圖2),CO2在混合氣中要比其單組份情況下吸附的量多,但不是特別明顯。而N2和O2在單組份時(shí)則要比在混合氣中被吸附的更多(圖3)。對于N2和O2而言,雖然它們的分壓高,但因其本身就屬于不被吸附的一類,因而在混合氣情況下,分子篩中大部分的吸附位都被CO2所占據(jù),留給它們的吸附位顯然就少了。而對于CO2來說,分子篩的大部分活性位都被其占據(jù),O2和N2只占據(jù)極少的活性位。但隨著CO2吸附的增加,一部分O2和N2分子也會(huì)吸附在CO2分子上,根據(jù)文獻(xiàn)[11]的觀點(diǎn),此時(shí)的吸附更多的是由于O2和N2在分壓上占據(jù)的絕對優(yōu)勢所帶來的一種堆積效應(yīng),從而為CO2的進(jìn)一步吸附提供了少量的活性位。也就是說,CO2不僅受分子篩骨架和自身的作用,還受已經(jīng)吸附上來的N2和O2分子的影響。因此,CO2在混合氣中的吸附量反而會(huì)有少量的上升,這一點(diǎn)也符合該文獻(xiàn)所提到的低壓下流體之間相互作用影響吸附量的觀點(diǎn)。當(dāng)然,這些都只是理論的分析,實(shí)際結(jié)果還有待于實(shí)驗(yàn)驗(yàn)證。

    5 結(jié) 論

    從某種程度上說,艦船上應(yīng)用的空氣凈化裝置材料的選擇要比凈化裝置本身結(jié)構(gòu)的設(shè)計(jì)更加重要。即使是同一種吸附劑,在不同溫度下,其吸附CO2的性能也會(huì)有所差異。本文選擇了一種微觀層面上的模擬,得到了一些重要結(jié)果:

    1)13X-APG分子篩對CO2有著良好的吸附效果。對于混有1.4%CO2的空氣而言,在總壓為100 kPa的情況下,CO2的吸附量能達(dá)到108.5~127.3 mg/g。

    2)該分子篩對空氣中N2和O2的吸附量極小,可以認(rèn)為是不吸附N2和O2。

    3)溫度對13X-APG分子篩的CO2吸附效果有一定的影響,在20~40℃溫度范圍內(nèi),吸附量減小了15%。因此,在設(shè)計(jì),尤其是選擇吸附材料時(shí),要充分考慮到溫度的影響。

    [1]彭光明.AIP潛艇艙室大氣環(huán)境控制系統(tǒng)研究[J].中國艦船研究,2006,1(2):62-65.PENG G M.A Study on the cabin atmosphere control system of AIP submarines[J].Chinese Journal of Ship Research,2006,1(2):62-65.

    [2]孟帥,徐筱欣.船用裝置三維流場的模擬研究[J].風(fēng)機(jī)技術(shù),2010(3):22-25.MENG S,XU X X.Numerical research on 3-D flow field for marine machinery[J].Compressor,Blower & Fan Technology,2010(3):22-25.

    [3]陳正龍,徐為人,湯立達(dá).分子模擬的理論與實(shí)踐[M].北京:化學(xué)工業(yè)出版社,2007.

    [4]姚敏強(qiáng),程智斌,董紅.武器任意攻擊角下艦艇的不沉性概率計(jì)算[J].中國艦船研究,2009,4(1):47-51.YAO M Q,CHENG Z B,DONG H.Calculation of the insinkability probability for ships in any attack angles[J].Chinese Journal of Ship Research,2009,4(1):47-51.

    [5]XIONG Y Z,SHEN X Z,ZHOU H,et al.Molecular simulation of pyridine derivatives sorption in faujasite zeolite[J].Computers and Applied Chemistry,2008,25 (12):1553-1556.

    [6]黃艷芳,馬正飛,劉曉勤,等.LTA和FAU型分子篩吸附CO2力場的研究[J].計(jì)算機(jī)與應(yīng)用化學(xué),2010,27(6):759-763.HUANG Y F,MA Z F,LIU X Q,et al.Force field study of adsorption of CO2in LTA and FAU zeolites[J].Computers and Applied Chemistry,2010,27(6):759-763.

    [7]朱玉俊,黃永明,胡軍,等.三維共價(jià)有機(jī)框架材料吸附分離CO2-CH4的計(jì)算機(jī)模擬——CO2電荷模型的影響[J].華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,37(1):8-15.ZHU Y J,HUANG Y M,HU J,et al.Computer simulation of CO2-CH4adsorption and separation in 3D-COFs—effect of CO2charge models[J].Journal of East China U-niversity of Science and Technology(Natural Science Edition),2011,37(1): 8-15.

    [8]FURUKAWA H,YAGHI O M.Storage of hydrogen,methane,and carbon dioxide in highly porous covalent organic frameworks for clean energy applications[J].Journal of the American Chemical Sociely,2009,131(25):8875-8883.

    [9]徐筱欣,施錫鉅.某空氣凈化器凈化材料篩選和阻力特性試驗(yàn)研究[J].實(shí)驗(yàn)室研究與探索,2004,23(7):14-16.XU X X,SHI X J.Experimental study on depuration materials selection and resistance characteristic of a certain air depuration device[J].Research and Exploration in Laboratory,2004,23(7):14-16.

    [10]TAKABA H,MATSUDA E,NAIR B N,et al.Molecular modeling of gas permeation through an amorphous microporous silica membrane[J].Journal of Chemical Engineering of Japan,2002,35(12):1312-1321.

    [11]GARBEROGLIO.Computer simulation of the adsorption of light gases in covalent organic frameworks[J].Langmuir,2007,23(24):12154-12158.

    猜你喜歡
    力場混合氣構(gòu)型
    基于物理模型的BaZrO3鈣鈦礦機(jī)器學(xué)習(xí)力場
    力場防護(hù)罩:并非只存在于科幻故事中
    調(diào)性的結(jié)構(gòu)力場、意義表征與聽覺感性先驗(yàn)問題——以貝多芬《合唱幻想曲》為例
    分子和離子立體構(gòu)型的判定
    脫氧核糖核酸柔性的分子動(dòng)力學(xué)模擬:Amber bsc1和bsc0力場的對比研究?
    航天器受迫繞飛構(gòu)型設(shè)計(jì)與控制
    ER308L焊絲硫含量對Ar-He混合氣TIG焊焊縫成形的影響
    焊接(2015年6期)2015-07-18 11:02:24
    遙感衛(wèi)星平臺與載荷一體化構(gòu)型
    兩個(gè)具stp三維拓?fù)錁?gòu)型的稀土配位聚合物{[Ln2(pda)3(H2O)2]·2H2O}n(Ln=Nd,La)
    Audi公司新一代1.8L增壓燃油分層噴射汽油機(jī)(第2部分)——混合氣形成、燃燒過程和增壓
    中文字幕高清在线视频| 亚洲在久久综合| 欧美亚洲日本最大视频资源| 日本欧美国产在线视频| 好男人视频免费观看在线| 精品国产国语对白av| 国产av国产精品国产| 国产精品一区二区精品视频观看| 精品一区二区三卡| 国产精品av久久久久免费| 一区在线观看完整版| 久久ye,这里只有精品| 啦啦啦视频在线资源免费观看| av一本久久久久| 天美传媒精品一区二区| 亚洲色图综合在线观看| 男男h啪啪无遮挡| av卡一久久| 丰满少妇做爰视频| 国产在视频线精品| 色播在线永久视频| 亚洲婷婷狠狠爱综合网| 午夜日韩欧美国产| 婷婷色综合www| 国产片特级美女逼逼视频| 下体分泌物呈黄色| 亚洲精品第二区| 天堂俺去俺来也www色官网| www.自偷自拍.com| 一区福利在线观看| 捣出白浆h1v1| 国产精品久久久久久人妻精品电影 | 高清av免费在线| 色综合欧美亚洲国产小说| 免费av中文字幕在线| 欧美中文综合在线视频| 丰满乱子伦码专区| 久久ye,这里只有精品| 精品少妇久久久久久888优播| 日日爽夜夜爽网站| 毛片一级片免费看久久久久| 国产 精品1| 免费黄网站久久成人精品| 亚洲精品在线美女| 久久久久精品性色| 亚洲一区中文字幕在线| 亚洲国产欧美网| 国产视频首页在线观看| 可以免费在线观看a视频的电影网站 | 久久久久久久久免费视频了| 亚洲精品久久午夜乱码| 啦啦啦在线免费观看视频4| 99久久综合免费| 国精品久久久久久国模美| 老司机深夜福利视频在线观看 | 日韩中文字幕欧美一区二区 | 少妇人妻 视频| 亚洲免费av在线视频| 色吧在线观看| 亚洲成国产人片在线观看| 亚洲欧洲日产国产| 欧美国产精品一级二级三级| 亚洲国产欧美在线一区| 国产亚洲午夜精品一区二区久久| 国产精品一区二区在线不卡| 校园人妻丝袜中文字幕| 亚洲色图综合在线观看| 国产免费一区二区三区四区乱码| 欧美乱码精品一区二区三区| 老司机影院成人| 最近2019中文字幕mv第一页| 一区二区三区激情视频| 国产探花极品一区二区| 欧美人与性动交α欧美精品济南到| 人人澡人人妻人| 色婷婷久久久亚洲欧美| 超碰97精品在线观看| www.自偷自拍.com| 国产亚洲av片在线观看秒播厂| 午夜福利乱码中文字幕| 中文字幕人妻丝袜一区二区 | 大码成人一级视频| 亚洲国产欧美在线一区| 国产日韩一区二区三区精品不卡| 青青草视频在线视频观看| 国产亚洲最大av| 久久久久久久久免费视频了| 一级毛片我不卡| 国产精品99久久99久久久不卡 | 伊人亚洲综合成人网| 日本wwww免费看| 国产一区亚洲一区在线观看| 丰满迷人的少妇在线观看| 人妻人人澡人人爽人人| 下体分泌物呈黄色| 777久久人妻少妇嫩草av网站| 亚洲五月色婷婷综合| 妹子高潮喷水视频| 男女午夜视频在线观看| 国产不卡av网站在线观看| 51午夜福利影视在线观看| 久久天躁狠狠躁夜夜2o2o | 男女下面插进去视频免费观看| a级毛片在线看网站| 成人国产麻豆网| 欧美日韩亚洲国产一区二区在线观看 | 国产精品欧美亚洲77777| 欧美日韩国产mv在线观看视频| 亚洲精品日本国产第一区| 最近2019中文字幕mv第一页| 亚洲av在线观看美女高潮| 看十八女毛片水多多多| 日韩大片免费观看网站| 亚洲欧美精品综合一区二区三区| 国产xxxxx性猛交| 肉色欧美久久久久久久蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲熟女精品中文字幕| 七月丁香在线播放| av天堂久久9| 一区二区三区四区激情视频| 久久久久久人人人人人| 免费观看人在逋| 十分钟在线观看高清视频www| 欧美国产精品va在线观看不卡| 精品一品国产午夜福利视频| 制服诱惑二区| 亚洲综合精品二区| 久久久欧美国产精品| 超色免费av| 看免费av毛片| 欧美成人精品欧美一级黄| 一区福利在线观看| 日韩熟女老妇一区二区性免费视频| av网站免费在线观看视频| 欧美人与性动交α欧美精品济南到| 国产男女内射视频| 久久精品久久久久久久性| 国产精品二区激情视频| 免费不卡黄色视频| 最近的中文字幕免费完整| 亚洲第一青青草原| 日本wwww免费看| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 最黄视频免费看| 亚洲国产欧美网| 精品亚洲成国产av| 国产成人精品福利久久| 天天躁夜夜躁狠狠久久av| 日韩制服骚丝袜av| 最近的中文字幕免费完整| 亚洲久久久国产精品| 日本一区二区免费在线视频| 曰老女人黄片| 久久久久精品人妻al黑| 亚洲av福利一区| 久热爱精品视频在线9| 精品一品国产午夜福利视频| 视频区图区小说| 亚洲精品国产av成人精品| 日韩欧美一区视频在线观看| 在线看a的网站| 国产成人免费观看mmmm| 又粗又硬又长又爽又黄的视频| 亚洲 欧美一区二区三区| 水蜜桃什么品种好| 欧美黑人欧美精品刺激| av卡一久久| 老汉色∧v一级毛片| 国产成人系列免费观看| 19禁男女啪啪无遮挡网站| 最新在线观看一区二区三区 | 色婷婷久久久亚洲欧美| 午夜福利视频在线观看免费| 国产xxxxx性猛交| 国产成人精品在线电影| 亚洲美女黄色视频免费看| 18禁裸乳无遮挡动漫免费视频| 色视频在线一区二区三区| 热re99久久精品国产66热6| 亚洲国产日韩一区二区| 熟女少妇亚洲综合色aaa.| 国产乱来视频区| 19禁男女啪啪无遮挡网站| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 久久国产精品大桥未久av| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av综合色区一区| 19禁男女啪啪无遮挡网站| 日韩中文字幕视频在线看片| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| 最近中文字幕高清免费大全6| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩亚洲综合一区二区三区_| 99精品久久久久人妻精品| 桃花免费在线播放| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 国产又色又爽无遮挡免| 国产成人91sexporn| 下体分泌物呈黄色| 欧美日韩综合久久久久久| 一级毛片电影观看| 韩国精品一区二区三区| 汤姆久久久久久久影院中文字幕| 日韩免费高清中文字幕av| 巨乳人妻的诱惑在线观看| 大香蕉久久成人网| 王馨瑶露胸无遮挡在线观看| avwww免费| 哪个播放器可以免费观看大片| 99久久99久久久精品蜜桃| 久久婷婷青草| 最新在线观看一区二区三区 | 男人舔女人的私密视频| 亚洲欧美一区二区三区黑人| 国产av国产精品国产| 国产一区有黄有色的免费视频| 成人手机av| 午夜福利视频精品| 熟女av电影| 十分钟在线观看高清视频www| 一边摸一边做爽爽视频免费| 免费女性裸体啪啪无遮挡网站| 久久久久久久久免费视频了| 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄| 久久久久精品性色| 香蕉丝袜av| 在线观看国产h片| 老熟女久久久| 日本av免费视频播放| 国产在视频线精品| 久久久国产精品麻豆| 一边亲一边摸免费视频| 亚洲精品国产av蜜桃| 一级片免费观看大全| 亚洲 欧美一区二区三区| 亚洲精品一区蜜桃| 成年美女黄网站色视频大全免费| 女人精品久久久久毛片| 中文字幕av电影在线播放| 国产欧美日韩一区二区三区在线| 久久久久久久久久久免费av| 免费av中文字幕在线| 久久久久久免费高清国产稀缺| 午夜福利,免费看| 视频在线观看一区二区三区| 成年女人毛片免费观看观看9 | 色视频在线一区二区三区| 亚洲成人免费av在线播放| 国产极品天堂在线| 国产精品嫩草影院av在线观看| 国产探花极品一区二区| 亚洲一级一片aⅴ在线观看| 少妇人妻精品综合一区二区| 欧美最新免费一区二区三区| 日本av手机在线免费观看| videos熟女内射| 国产精品无大码| 国产精品嫩草影院av在线观看| 黄色一级大片看看| 综合色丁香网| 精品人妻在线不人妻| 成人国产麻豆网| 水蜜桃什么品种好| 大话2 男鬼变身卡| 久久精品国产综合久久久| 高清视频免费观看一区二区| 黑人猛操日本美女一级片| 久久精品国产a三级三级三级| 美女大奶头黄色视频| 午夜福利视频精品| 国产欧美亚洲国产| 久久久久久久久免费视频了| 赤兔流量卡办理| 国产成人免费观看mmmm| 亚洲av成人不卡在线观看播放网 | 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 97在线人人人人妻| 成人国产麻豆网| 激情视频va一区二区三区| 欧美日韩综合久久久久久| 久久久久久久精品精品| 1024香蕉在线观看| 激情五月婷婷亚洲| 久久精品熟女亚洲av麻豆精品| 视频区图区小说| 国产精品一二三区在线看| 亚洲一码二码三码区别大吗| 日日爽夜夜爽网站| 午夜免费男女啪啪视频观看| 国产精品欧美亚洲77777| 波多野结衣av一区二区av| 观看av在线不卡| 国产日韩欧美视频二区| 纵有疾风起免费观看全集完整版| 免费日韩欧美在线观看| av.在线天堂| 欧美久久黑人一区二区| 国产乱人偷精品视频| av在线观看视频网站免费| 亚洲欧美清纯卡通| 国产有黄有色有爽视频| 在线观看免费高清a一片| 老熟女久久久| 丝袜人妻中文字幕| 国产成人精品在线电影| 丰满乱子伦码专区| 午夜福利一区二区在线看| 七月丁香在线播放| av卡一久久| 日韩欧美一区视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 777久久人妻少妇嫩草av网站| 丝袜美腿诱惑在线| 好男人视频免费观看在线| 国产不卡av网站在线观看| 美女扒开内裤让男人捅视频| 亚洲第一av免费看| 亚洲国产欧美在线一区| 日本欧美视频一区| 777米奇影视久久| 99久久人妻综合| 国产精品一国产av| 性色av一级| 男女下面插进去视频免费观看| 亚洲色图综合在线观看| 午夜精品国产一区二区电影| 叶爱在线成人免费视频播放| 狂野欧美激情性xxxx| 亚洲精品成人av观看孕妇| 女人久久www免费人成看片| 精品久久久久久电影网| www日本在线高清视频| 熟女av电影| e午夜精品久久久久久久| 妹子高潮喷水视频| 91aial.com中文字幕在线观看| 最近的中文字幕免费完整| 亚洲精品aⅴ在线观看| 91精品国产国语对白视频| 久久久国产欧美日韩av| 街头女战士在线观看网站| 中文字幕高清在线视频| 精品一区二区三卡| 国产 一区精品| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 久久天躁狠狠躁夜夜2o2o | netflix在线观看网站| 日本黄色日本黄色录像| www.熟女人妻精品国产| 伊人久久大香线蕉亚洲五| 亚洲一码二码三码区别大吗| 两个人看的免费小视频| 亚洲av在线观看美女高潮| 伊人亚洲综合成人网| 一二三四中文在线观看免费高清| 天天躁狠狠躁夜夜躁狠狠躁| av在线观看视频网站免费| 热re99久久国产66热| 久久久精品免费免费高清| 亚洲欧美成人综合另类久久久| 老司机靠b影院| 91国产中文字幕| 久久人人爽人人片av| 男女边摸边吃奶| 黄色视频在线播放观看不卡| 麻豆av在线久日| 中文乱码字字幕精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 午夜91福利影院| 亚洲,欧美,日韩| 亚洲国产中文字幕在线视频| www.自偷自拍.com| 亚洲av国产av综合av卡| av卡一久久| 国产不卡av网站在线观看| 国产成人91sexporn| 免费黄色在线免费观看| 久久精品人人爽人人爽视色| 国产高清国产精品国产三级| 新久久久久国产一级毛片| 少妇人妻久久综合中文| 男女午夜视频在线观看| xxx大片免费视频| av在线老鸭窝| 国产男人的电影天堂91| 久久精品aⅴ一区二区三区四区| 久久久国产一区二区| 久久狼人影院| 国产视频首页在线观看| 久久久欧美国产精品| 搡老乐熟女国产| 久久毛片免费看一区二区三区| 日韩中文字幕视频在线看片| 熟妇人妻不卡中文字幕| 亚洲国产欧美一区二区综合| 国产av国产精品国产| 欧美av亚洲av综合av国产av | 精品少妇一区二区三区视频日本电影 | 午夜福利免费观看在线| 熟女少妇亚洲综合色aaa.| 一区福利在线观看| 国产精品欧美亚洲77777| 另类亚洲欧美激情| 色视频在线一区二区三区| 捣出白浆h1v1| 国产免费一区二区三区四区乱码| 国产精品av久久久久免费| 又黄又粗又硬又大视频| 久久精品国产亚洲av高清一级| 不卡av一区二区三区| 亚洲精品aⅴ在线观看| 一边摸一边抽搐一进一出视频| 中文天堂在线官网| 精品少妇一区二区三区视频日本电影 | √禁漫天堂资源中文www| 欧美97在线视频| 国产精品欧美亚洲77777| kizo精华| 日本av手机在线免费观看| 国产欧美亚洲国产| 午夜91福利影院| 男女之事视频高清在线观看 | 国产日韩欧美视频二区| 日本午夜av视频| 中文字幕亚洲精品专区| 国产一区二区在线观看av| 国产乱人偷精品视频| 成年女人毛片免费观看观看9 | 国产av一区二区精品久久| 国产片内射在线| 18禁国产床啪视频网站| 丰满迷人的少妇在线观看| 久久久久精品人妻al黑| av卡一久久| 大香蕉久久网| 最近中文字幕高清免费大全6| 国产毛片在线视频| 高清不卡的av网站| 美女脱内裤让男人舔精品视频| 中文字幕精品免费在线观看视频| 大话2 男鬼变身卡| 国产乱来视频区| av免费观看日本| 亚洲av在线观看美女高潮| 精品久久蜜臀av无| h视频一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 一级,二级,三级黄色视频| 亚洲精品在线美女| 亚洲欧美精品自产自拍| 综合色丁香网| 欧美变态另类bdsm刘玥| 一级片'在线观看视频| 久久久久久久国产电影| 又大又爽又粗| 老司机在亚洲福利影院| 免费观看人在逋| 亚洲欧美一区二区三区久久| 亚洲男人天堂网一区| 熟女少妇亚洲综合色aaa.| 国产精品久久久av美女十八| 看非洲黑人一级黄片| 老汉色∧v一级毛片| 1024视频免费在线观看| av不卡在线播放| 五月开心婷婷网| 午夜免费鲁丝| 欧美黄色片欧美黄色片| 久久久久视频综合| 亚洲图色成人| 国产熟女午夜一区二区三区| 国产精品一二三区在线看| 十八禁人妻一区二区| 国产亚洲av高清不卡| 亚洲av欧美aⅴ国产| 国产亚洲一区二区精品| 日本vs欧美在线观看视频| 精品视频人人做人人爽| 校园人妻丝袜中文字幕| 欧美激情高清一区二区三区 | 中文字幕人妻熟女乱码| 日本91视频免费播放| 日本猛色少妇xxxxx猛交久久| 欧美黄色片欧美黄色片| 精品卡一卡二卡四卡免费| 国产女主播在线喷水免费视频网站| 日韩中文字幕欧美一区二区 | 午夜91福利影院| 精品国产乱码久久久久久男人| 99久久99久久久精品蜜桃| 亚洲人成77777在线视频| 亚洲成人手机| 精品酒店卫生间| 精品一区二区三区av网在线观看 | 成人影院久久| 亚洲中文av在线| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 久久精品国产a三级三级三级| 国产激情久久老熟女| 人人妻人人爽人人添夜夜欢视频| 国语对白做爰xxxⅹ性视频网站| 久久天躁狠狠躁夜夜2o2o | 国产一区二区激情短视频 | 午夜日韩欧美国产| 天天躁夜夜躁狠狠久久av| 久久久久久久大尺度免费视频| 国产精品人妻久久久影院| 欧美精品高潮呻吟av久久| 亚洲综合色网址| 日本欧美视频一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产区一区二| 一级,二级,三级黄色视频| 亚洲国产精品成人久久小说| 日韩熟女老妇一区二区性免费视频| 青青草视频在线视频观看| 宅男免费午夜| 精品一区在线观看国产| 久久毛片免费看一区二区三区| 男女边吃奶边做爰视频| 80岁老熟妇乱子伦牲交| 久久ye,这里只有精品| 色94色欧美一区二区| 国产精品99久久99久久久不卡 | 亚洲欧美色中文字幕在线| 久久ye,这里只有精品| 午夜福利在线免费观看网站| 久热爱精品视频在线9| 一区二区av电影网| 午夜免费男女啪啪视频观看| netflix在线观看网站| 亚洲av中文av极速乱| 伊人久久国产一区二区| 国产精品一区二区在线不卡| 国产精品av久久久久免费| 亚洲三区欧美一区| 国产精品亚洲av一区麻豆 | 亚洲国产最新在线播放| 欧美日韩成人在线一区二区| 日本91视频免费播放| 国产一区二区 视频在线| 中文字幕另类日韩欧美亚洲嫩草| 妹子高潮喷水视频| 大码成人一级视频| 啦啦啦 在线观看视频| 在现免费观看毛片| 少妇被粗大猛烈的视频| 精品一区二区三卡| 卡戴珊不雅视频在线播放| 久久免费观看电影| 人人妻,人人澡人人爽秒播 | 丁香六月欧美| 人成视频在线观看免费观看| 久久久久久人人人人人| 国产精品欧美亚洲77777| 亚洲av电影在线进入| 最近2019中文字幕mv第一页| 日韩一卡2卡3卡4卡2021年| a级片在线免费高清观看视频| 99九九在线精品视频| 这个男人来自地球电影免费观看 | 欧美精品av麻豆av| 99精国产麻豆久久婷婷| 啦啦啦啦在线视频资源| 久久久久视频综合| 成年人免费黄色播放视频| 成年人午夜在线观看视频| 国产欧美日韩综合在线一区二区| 欧美变态另类bdsm刘玥| 美国免费a级毛片| 国产 一区精品| avwww免费| 亚洲国产欧美网| 国产精品偷伦视频观看了| 国产日韩欧美在线精品| 高清不卡的av网站| 视频在线观看一区二区三区| 日本91视频免费播放| xxxhd国产人妻xxx| 永久免费av网站大全| 国产欧美日韩一区二区三区在线| www.熟女人妻精品国产| 你懂的网址亚洲精品在线观看| 国产成人欧美| 美女扒开内裤让男人捅视频| 亚洲人成77777在线视频| 啦啦啦啦在线视频资源| 一边亲一边摸免费视频| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 国产男女内射视频| 色网站视频免费| 性色av一级| 欧美日韩亚洲国产一区二区在线观看 | 成年动漫av网址| 国产免费视频播放在线视频| 色精品久久人妻99蜜桃| 久久久国产欧美日韩av| 啦啦啦在线免费观看视频4| 亚洲av中文av极速乱| 午夜影院在线不卡| 久久韩国三级中文字幕| 亚洲精品,欧美精品| 国产精品麻豆人妻色哟哟久久| 日本爱情动作片www.在线观看| 久久久亚洲精品成人影院|