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

    碳納米管受介電泳作用三維運(yùn)動(dòng)仿真

    2016-12-06 07:59:02王小沖安立寶
    關(guān)鍵詞:電液電泳碳納米管

    王小沖,安立寶,龔 亮,陳 琰

    (華北理工大學(xué)機(jī)械工程學(xué)院,河北省無(wú)機(jī)非金屬材料重點(diǎn)實(shí)驗(yàn)室,河北唐山 063009)

    碳納米管受介電泳作用三維運(yùn)動(dòng)仿真

    王小沖,安立寶,龔 亮,陳 琰

    (華北理工大學(xué)機(jī)械工程學(xué)院,河北省無(wú)機(jī)非金屬材料重點(diǎn)實(shí)驗(yàn)室,河北唐山 063009)

    利用溶于介電液中的碳納米管在外加交變電場(chǎng)作用下極化產(chǎn)生介電泳力和液體粘滯阻力的共同作用,建立了碳納米管受力和運(yùn)動(dòng)模型;對(duì)碳納米管的運(yùn)動(dòng)過(guò)程進(jìn)行了仿真,得到了碳納米管從特定初始點(diǎn)開(kāi)始的三維運(yùn)動(dòng)軌跡和可以實(shí)現(xiàn)組裝的碳納米管初始點(diǎn)分布區(qū)域.計(jì)算了碳納米管運(yùn)動(dòng)過(guò)程中所受介電泳力和速度的變化規(guī)律,距離電極間隙越近,碳納米管受力和速度越大,最大分別能達(dá)到10-9N和105μm/s的數(shù)量級(jí).仿真結(jié)果為碳納米管介電電泳組裝提供指導(dǎo).

    碳納米管;介電泳;數(shù)學(xué)模型;仿真;三維運(yùn)動(dòng)軌跡

    碳納米管(Carbon Nano Tube,CNT)具備優(yōu)異的電學(xué)、力學(xué)、熱學(xué)性能以及良好的化學(xué)穩(wěn)定性,因此,可被應(yīng)用于眾多領(lǐng)域,如場(chǎng)效應(yīng)晶體管、場(chǎng)發(fā)射元件、互連線、傳感器、催化劑載體以及新能源電極材料等[1-3].組裝是上述大多應(yīng)用得以實(shí)現(xiàn)的重要前提條件[4-5].碳納米管組裝常用的方法包括化學(xué)氣相沉積[6]、流體排布法[7]、原子力顯微鏡直接組裝法[8]和介電電泳組裝[9]等.相對(duì)于其他組裝方法,介電電泳組裝具有操作簡(jiǎn)單、速度快和成本低的優(yōu)點(diǎn).因此,介電電泳法成為目前廣泛使用的碳納米管組裝方法.

    在外加交變電場(chǎng)作用下,介電液中的碳納米管將被極化而產(chǎn)生介電泳力[10].此外,碳納米管在介電液中運(yùn)動(dòng),還會(huì)受到液體的粘滯阻力.為了更好地分析碳納米管的介電電泳組裝過(guò)程,筆者將建立碳納米管的受力和運(yùn)動(dòng)數(shù)學(xué)模型,并依據(jù)所建模型對(duì)組裝過(guò)程中碳納米管的運(yùn)動(dòng)進(jìn)行仿真,分析碳納米管在運(yùn)動(dòng)過(guò)程中的受力和速度變化規(guī)律,得到碳納米管三維運(yùn)動(dòng)軌跡,獲取可以成功組裝碳納米管的初始點(diǎn)分布區(qū)域.

    1 介電電泳模型

    現(xiàn)有的介電電泳模型大都是基于有效偶極矩理論[11]建立的.該理論只有在偶極子長(zhǎng)度小于電場(chǎng)的非均勻度,即“偶極子近似”條件成立時(shí)才適用.碳納米管的直徑為納米級(jí),長(zhǎng)度一般在微米級(jí)[12],為長(zhǎng)徑比較大的一維材料.在距離電極間隙較近的地方,碳納米管上各點(diǎn)所處位置的電場(chǎng)強(qiáng)度變化很大,即作為偶極子的碳納米管的長(zhǎng)度大于電場(chǎng)的非均勻度,“偶極子近似”條件不成立使得有效偶極矩理論不再適用.因此,需要建立新的適用于碳納米管等一維納米材料組裝的介電電泳模型.

    碳納米管是由碳原子組成的一維納米材料,可以假設(shè)碳納米管為連續(xù)介質(zhì),對(duì)碳納米管進(jìn)行分段處理,劃分后的每段滿足“偶極子近似”條件.根據(jù)每段碳納米管所在位置的局部電場(chǎng),依據(jù)有效偶極矩理論計(jì)算出各段的介電泳力,然后再把各段介電泳力相加得到碳納米管所受總的介電泳力[13].考慮到碳納米管在電場(chǎng)中兩端部分電場(chǎng)強(qiáng)度變化較大,中間部分場(chǎng)強(qiáng)變化較小[14],因此,兩端部分的劃分要比中間部分精細(xì).

    碳納米管中間電場(chǎng)強(qiáng)度變化較小的部分劃分成N1個(gè)圓柱體.對(duì)場(chǎng)強(qiáng)變化較大的碳納米管兩端進(jìn)行精細(xì)劃分,兩端部分不再劃分成圓柱體,而是劃分成球體如圖1所示,兩端分別劃分成N2/2個(gè)球體.

    圖1 碳納米管分段示意圖

    第i(i=1,2,3,…,N1)個(gè)圓柱體所受的介電泳力為[15]

    其中,r為碳納米管的半徑;li=(L-N2×2×r)/N1,為第i個(gè)圓柱體的軸向長(zhǎng)度,其中L為碳納米管的長(zhǎng)度;Ei為第i個(gè)圓柱體的局部電場(chǎng)強(qiáng)度;和,分別為碳納米管和介電液的復(fù)合介電常數(shù),εCNT和εm分別為碳納米管和介電液的介電常數(shù),σCNT和σm分別為碳納米管和介電液的電導(dǎo)率,ω=2πfe,fe為電場(chǎng)頻率;,為去偏極系數(shù).

    第n(n=1,2,3,…,N2)個(gè)球體所受的介電泳力為[16]

    其中,En為第n個(gè)球體的局部電場(chǎng)強(qiáng)度.

    碳納米管所受到的總介電泳力為所有圓柱體段和球體段所受介電泳力之和FDEP.碳納米管在介電液中受介電泳力作用運(yùn)動(dòng)時(shí)還會(huì)受到液體粘滯阻力的影響,碳納米管的運(yùn)動(dòng)速度為[17]

    在碳納米管組裝過(guò)程中,碳納米管的軸向方向總是沿著電場(chǎng)線的方向[18].碳納米管組裝到電極上后由于靜電力[19]和范德華力[20]的作用,碳納米管被固定在電極上.碳納米管自身重力大約為10-19N[17],相對(duì)于介電泳力(10-14~10-9N)而言,可以忽略不計(jì).當(dāng)粒子直徑大于1μm時(shí),布朗運(yùn)動(dòng)對(duì)粒子運(yùn)動(dòng)的影響非常小,布朗運(yùn)動(dòng)作用同樣可以忽略不計(jì)[11].

    2 仿真研究

    仿真研究容易重復(fù)操作,研究成本低,還可以排除復(fù)雜的實(shí)驗(yàn)環(huán)境帶來(lái)的影響.根據(jù)上述的介電電泳模型建立仿真模型,設(shè)定仿真參數(shù),利用COMSOL Multiphysics軟件進(jìn)行仿真研究.在仿真結(jié)果中對(duì)電場(chǎng)和電場(chǎng)平方梯度的分布、碳納米管運(yùn)動(dòng)過(guò)程中的受力和速度變化情況、碳納米管的運(yùn)動(dòng)軌跡情況和可以成功組裝碳納米管的初始位置分布區(qū)域進(jìn)行分析,得到了碳納米管的受力和速度與碳納米管到電極間隙距離之間的關(guān)系,以及形象直觀的三維運(yùn)動(dòng)軌跡圖.為了加快運(yùn)算速度,根據(jù)電場(chǎng)的分布情況,在距離電極間隙遠(yuǎn)的地方仿真步長(zhǎng)設(shè)置的較大,在距離近的地方步長(zhǎng)設(shè)置的較小.

    2.1仿真模型

    圖2為所建立的仿真模型圖,模型中硅基片長(zhǎng)寬均為60μm,厚度為10μm;二氧化硅覆蓋在硅基片上,厚度為1μm;電極為金屬鋁,電極寬度為2μm,厚度為0.1μm,電極間隙為2μm.在電極上施加幅值為5 V, fe=100 k Hz的交流電,溶液完全覆蓋二氧化硅基片,高度設(shè)置為30μm.

    圖2 三維運(yùn)動(dòng)軌跡仿真模型圖

    2.2仿真參數(shù)

    仿真過(guò)程中碳納米管設(shè)定為L(zhǎng)=4μm,r=20 nm的多壁碳納米管, N1=20,N2=6.介電溶液物理參數(shù)對(duì)仿真結(jié)果有著重要影響,文中綜合文獻(xiàn)選取異丙醇為介電液.異丙醇的粘度系數(shù)η=1×10-3Pa·s[14],真空絕對(duì)介電常數(shù)ε0=8.85×10-12F/m[10],各種材料電學(xué)參數(shù)如表1所示[14].

    表1 模型中各種材料的電學(xué)參數(shù)[14]

    2.3仿真結(jié)果

    X-Y平面(Z=0)上電場(chǎng)強(qiáng)度的分布情況如圖3(a)所示,從圖中可以看出,在電極間隙處電場(chǎng)強(qiáng)度最強(qiáng).在文中碳納米管受正介電泳力作用,向電場(chǎng)強(qiáng)度大的區(qū)域運(yùn)動(dòng)[13],即碳納米管向電極間隙運(yùn)動(dòng).圖3(b)仿真Y軸(X=0,Z=0)方向上Y=[-20μm,20μm]范圍內(nèi)的電場(chǎng)平方梯度分布情況,發(fā)現(xiàn)在電極間隙靠近電極尖端處電場(chǎng)平方梯度的值最大.由式(1)和式(2)可知,當(dāng)碳納米管到達(dá)該區(qū)域時(shí)受到的介電泳力也最大,同時(shí)在電極間隙的中間區(qū)域電場(chǎng)平方梯度變化較平穩(wěn),在此區(qū)域更容易實(shí)現(xiàn)碳納米管組裝,與文獻(xiàn)[21]得出的結(jié)果一致.

    圖3 電場(chǎng)強(qiáng)度分布和電場(chǎng)平方梯度分布

    以碳納米管從初始點(diǎn)(14.44μm,15μm,25μm)開(kāi)始運(yùn)動(dòng)為例,對(duì)組裝過(guò)程中碳納米管受力和速度變化情況進(jìn)行分析.從圖4(a)中可以看出,隨著碳納米管到電極間隙距離的減小,碳納米管所受的介電泳力逐漸增大,最大能達(dá)到10-9N的數(shù)量級(jí).圖4(b)描述的是隨著碳納米管到電極間隙距離的減小,碳納米管的速度也在逐漸增大,最大速度能達(dá)到105μm/s的數(shù)量級(jí),與文獻(xiàn)[14]得出的結(jié)果近似.比較圖4(a)和(b)發(fā)現(xiàn),碳納米管受力和速度隨距離的變化趨勢(shì)是一致的,即當(dāng)碳納米管受力較大時(shí),運(yùn)動(dòng)速度也會(huì)相應(yīng)地很大.

    圖5(a)所示為碳納米管從初始點(diǎn)(14.44μm,15μm,25μm)開(kāi)始運(yùn)動(dòng)的三維運(yùn)動(dòng)軌跡.從圖中可以看出,碳納米管在組裝過(guò)程中不斷向電極間隙位置靠近,其軸線方向和電場(chǎng)線方向始終保持一致.組裝完成時(shí)碳納米管左端和左側(cè)電極剛好接觸,同時(shí)順著電場(chǎng)線方向排列,如圖5(b)所示,與文獻(xiàn)[18]實(shí)驗(yàn)觀察到的結(jié)果相同.點(diǎn)(14.44μm,15μm,25μm)就是可以成功組裝碳納米管到電極間隙上的空間區(qū)域的一個(gè)臨界位置.

    圖4 介電泳力和速度隨碳納米管到電極間隙距離的變化

    圖5(c)描繪了碳納米管從初始點(diǎn)(0μm,0μm,25μm)開(kāi)始運(yùn)動(dòng)的三維運(yùn)動(dòng)軌跡.從圖中可以看出,碳納米管在組裝過(guò)程中一直處在電極間隙正上方的位置,并且軸線方向和電場(chǎng)線的方向保持一致.碳納米管組裝完成后兩端與電極接觸長(zhǎng)度相等,同時(shí)順著電場(chǎng)線方向排列,如圖5(d)所示.由圖5可知,碳納米管如果從比臨界位置更靠近電極間隙的初始點(diǎn)開(kāi)始運(yùn)動(dòng),則碳納米管能夠組裝到電極上.

    圖5 碳納米管三維運(yùn)動(dòng)軌跡及最后組裝位置示意圖

    在坐標(biāo)系統(tǒng)的第一象限內(nèi)找到高度Z=25μm的6個(gè)初始點(diǎn),從這6個(gè)初始點(diǎn)開(kāi)始運(yùn)動(dòng)的碳納米管恰好可以組裝到電極上,如圖6所示,這6個(gè)點(diǎn)所在位置就是碳納米管恰好能夠組裝到電極上的臨界位置.從圖6碳納米管的6條運(yùn)動(dòng)軌跡可以看出,如果初始點(diǎn)在6條運(yùn)動(dòng)軌跡的包絡(luò)面與YZ面、XZ面所圍成的區(qū)域內(nèi),碳納米管就可以組裝到電極間隙上.由于模型中電場(chǎng)是對(duì)稱分布的,其余3個(gè)象限內(nèi)的碳納米管恰好可以組裝到電極上的臨界點(diǎn)與第一象限內(nèi)的臨界點(diǎn)是對(duì)稱分布的.4個(gè)象限內(nèi)的臨界點(diǎn)及其運(yùn)動(dòng)軌跡可以圍成一個(gè)三維空間初始點(diǎn)分布區(qū)域,如果碳納米管的初始點(diǎn)在這個(gè)區(qū)域內(nèi),則碳納米管就可以實(shí)現(xiàn)組裝.

    圖6 第一象限內(nèi)的6個(gè)臨界初始點(diǎn)及碳納米管運(yùn)動(dòng)軌跡

    3 結(jié) 論

    筆者建立了適合碳納米管等一維納米材料的介電電泳組裝數(shù)學(xué)模型.通過(guò)對(duì)碳納米管分段減小了偶極子長(zhǎng)度,使每段碳納米管滿足“偶極子近似”條件.采用每段碳納米管所在處的局部電場(chǎng)計(jì)算各段所受介電泳力,然后累加求得碳納米管所受總的介電泳力,提高了模型的計(jì)算精度.仿真中考慮了介電泳力和液體粘滯阻力的共同作用,碳納米管所受重力和布朗運(yùn)動(dòng)作用很小而忽略不計(jì).

    隨著碳納米管到電極間隙距離的減小,碳納米管所受的介電泳力和運(yùn)動(dòng)速度增大,兩者隨距離變化曲線的趨勢(shì)一致,最大分別達(dá)到10-9N和105μm/s的數(shù)量級(jí),且碳納米管在運(yùn)動(dòng)過(guò)程中其軸線方向與電場(chǎng)線方向始終保持一致.文中得到了碳納米管三維運(yùn)動(dòng)軌跡和可以成功組裝碳納米管的初始點(diǎn)分布區(qū)域,揭示了碳納米管在介電泳力和液體粘滯阻力共同作用下的運(yùn)動(dòng)規(guī)律,將為進(jìn)一步分析碳納米管介電電泳組裝過(guò)程和參數(shù)優(yōu)化提供幫助.

    [1]李玉魁,李德昌,朱長(zhǎng)純.新型陰極粘貼技術(shù)制備三極CNT-FED的研究[J].西安電子科技大學(xué)學(xué)報(bào),2005,32(2): 314-317. LI Yukui,LI Dechang,ZHU Changchun.The Fabrication of the Triode CNT-FED with New Cathode-pasting Technology[J].Journal of Xidian University,2005,32(2):314-317.

    [2]RAKOV E G.Carbon Nanotubes in New Materials[J].Russian Chemical Reviews,2013,82(1):27-41.

    [3]劉紅霞,宋久旭,張鶴鳴.(4,4)碳納米管/碳化硅納米管異質(zhì)結(jié)的輸運(yùn)特性[J].西安電子科技大學(xué)學(xué)報(bào),2010,37 (3):520-523. LIU Hongxia,SONG Jiuxu,ZHANG Heming.Electronic Transport Properties of a(4,4)Carbon Nanotube/Silicon Carbide Nanotube Heterojunction[J].Journal of Xidian University,2010,37(3):520-523.

    [4]NAIENI A K,NOJEH A.Dielectrophoretic Deposition of Carbon Nanotubes:The Role of Field Frequency and Its Dependence on Solution Conductivity[J].Microelectronic Engineering,2014,114(2):26-30.

    [5]YAN Y H,CHAN-PARK M B,ZHANG Q.Advances in Carbon-Nanotube Assembly[J].Small,2007,3(1):24-42.

    [6]CHEN Y B,ZHANG J.Chemical Vapor Deposition Growth of Single-Walled Carbon Nanotubes with Controlled Structures for Nanodevice Applications[J].Accounts of Chemical Research,2014,47(8):2273-2281.

    [7]譚苗苗,葉雄英,王曉浩,等.流體排布與介電電泳結(jié)合排列組裝碳納米管[J].納米技術(shù)與精密工程,2010,8(5): 389-393. TAN Miaomiao,YE Xiongying,WANG Xiaohao,et al.Combining Fluidic Alignment with Dielectrophoresis to Assemble Carbon Nanotubes[J].Nanotechnology and Precision Engineering,2010,8(5):389-393.

    [8]POSTMA H W C,SELLMEIJER A,DEKKER C.Manipulation and Imaging of Individual Single-walled Carbon Nanotubes with an Atomic Force Microscope[J].Advanced Materials,2000,12(17):1299-1302.

    [9]AN L,FRIEDRICH C R.Dielectrophoretic Assembly of Carbon Nanotubes and Stability Analysis[J].Progress in Natural Science:Materials International,2013,23(4):367-373.

    [10]OLIVA-AVILéS A I,AVILéS F,SOSA V,et al.Dielectrophoretic Modeling of the Dynamic Carbon Nanotube Network Formation in Viscous Media under Alternating Current Electric Fields[J].Carbon,2014,69(4):342-354.

    [11]MORGAN H,GREEN N G.AC Electrokinetics:Colloids and Nanoparticles[M].Philadelphia:Research Studies Press,2003.

    [12]李德昌,李玉魁,李昕,等.碳納米管薄膜陰極的制備及分析[J].西安電子科技大學(xué)學(xué)報(bào),2005,32(6):877-880. LI Dechang,LI Yukui,LI Xin,et al.The Preparation and Analysis of the Carbon Nanotube Film Cathode[J].Journal of Xidian University,2005,32(6):877-880.

    [13]TAO Q,JIANG M,LI G.Simulation and Experimental Study of Nanowire Assembly by Dielectrophoresis[J].IEEE Transactions on Nanotechnology,2014,13(3):517-526.

    [14]AN L,FRIEDRICH C R.Process Parameters and Their Relations for the Dielectrophoretic Assembly of CarbonNanotubes[J].Journal of Applied Physics,2009,105(7):074314.

    [15]BERGER S D,MCGRUER N E,ADAMS G G.Simulation of Dielectrophoretic Assembly of Carbon Nanotubes Using 3D Finite Element Analysis[J].Nanotechnology,2015,26(15):155602.

    [16]BANERJEE S,WHITE B,HUANG L,et al.Precise Positioning of Carbon Nanotubes by Ac Dielectrophoresis Using Floating Posts[J].Applied Physics A,2007,86(4):415-419.

    [17]KIM J E,HAN C S.Use of Dielectrophoresis in the Fabrication of an Atomic Force Microscope Tip with a Carbon Nanotube:a Numerical Analysis[J].Nanotechnology,2005,16(10):2245-2250.

    [18]YANG B,YANG Z,ZHAO Z,et al.The Assembly of Carbon Nanotubes by Dielectrophoresis:Insights into the Dielectrophoretic Nanotube--Nanotube Interactions[J].Physica E,Low-dimentional Systems and Nanostructures,2014, 56(2):117-122.

    [19]MAYER A,LAMBIN P.Calculation of the Electrostatic Forces that Act on Carbon Nanotubes Placed in the Vicinity of Metallic Protrusions[J].Nanotechnology,2005,16(11):2685-2695.

    [20]HIROTANI J,IKUTA T,NISHIYAMA T,et al.Measuring the Thermal Boundary Resistance of Van Der Waals Contacts Using an Individual Carbon Nanotube[J].Journal of Physics:Condensed Matter,2013,25(2):025301.

    [21]KWON S G,KIM S H,KIM K H,et al.Distribution of Electric Field for Carbon Nanotube Assembly:Simulation(I) [J].The Chinese Journal of Nonferrous Metals,2011,21(11):117-120.

    (編輯:李恩科)

    Three-dimensional simulation of kinetics of carbon nanotubes under dielectrophoresis

    WANG Xiaochong,AN Libao,GONG Liang,CHEN Yan
    (College of Mechanical Engineering,Hebei Provincial Key Lab.of Inorganic Nonmetallic Materials, North China Univ.of Science and Technology,Tangshan 063009,China)

    Based on the dielectrophoretic force arising from the polarization of carbon nanotubes(CNTs)in the dielectric fluid under an alternating current electric field and the viscous resistance from the fluid,force and motion models for CNTs under dielectrophoresis(DEP)are established.Simulation is conducted on the translational and rotational motion of CNTs.Three-dimensional trajectories and a distribution region of initial points from where a CNT can be successfully assembled onto the electrode gap are obtained.By calculating the DEP force and translational velocity of a CNT during assembly,we know that the nearer the CNT to the electrode gap,the greater the DEP force and CNT velocity,and that the maximal DEP force and CNT velocity can be respectively of orders 10-9N and 105μm/s.The results of simulation provide guidance for DEP assembly of CNTs.

    carbon nanotube;dielectrophoresis;mathematical model;simulation;three-dimensional trajectory

    TN4

    A

    1001-2400(2016)04-0154-06

    10.3969/j.issn.1001-2400.2016.04.027

    2015-04-09 網(wǎng)絡(luò)出版時(shí)間:2015-10-21

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51172062,51472074);河北省引進(jìn)海外高層次人才“百人計(jì)劃”資助項(xiàng)目(E2012100005)

    王小沖(1990-),男,華北理工大學(xué)碩士研究生,E-mail:purplean@163.com.

    安立寶(1965-),男,教授,博士,E-mail:lan@ncst.edu.cn.

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1076.TN.20151021.1046.054.html

    猜你喜歡
    電液電泳碳納米管
    北京航空航天大學(xué)學(xué)報(bào)(2022年6期)2022-07-02 01:59:46
    電液復(fù)合控制系統(tǒng)SEM工作原理及測(cè)試方法
    化工管理(2021年7期)2021-05-13 00:46:16
    挖掘機(jī)屬具電液控制系統(tǒng)的研究開(kāi)發(fā)
    輔助陽(yáng)極在輕微型廂式車身電泳涂裝中的應(yīng)用
    PPG第四屆電泳涂料研討會(huì)在長(zhǎng)沙成功舉辦
    上海建材(2017年4期)2017-04-06 07:32:03
    改良的Tricine-SDS-PAGE電泳檢測(cè)胸腺肽分子量
    碳納米管陣列/環(huán)氧樹(shù)脂的導(dǎo)熱導(dǎo)電性能
    基于AEMSim的某重型AT電液調(diào)壓閥仿真分析
    專用汽車(2016年4期)2016-03-01 04:14:27
    聚賴氨酸/多壁碳納米管修飾電極測(cè)定大米中的鉛
    拓?fù)淙毕輰?duì)Armchair型小管徑多壁碳納米管輸運(yùn)性質(zhì)的影響
    亚洲国产精品一区三区| 岛国在线观看网站| 丰满少妇做爰视频| 日日夜夜操网爽| 高清在线国产一区| 亚洲国产欧美在线一区| 久久精品久久久久久噜噜老黄| 欧美成狂野欧美在线观看| 亚洲五月婷婷丁香| 久久久国产成人免费| 午夜福利在线观看吧| 国产成人免费观看mmmm| 一区二区三区乱码不卡18| 久久久久久人人人人人| 捣出白浆h1v1| 99精品久久久久人妻精品| 成人影院久久| 国产精品一二三区在线看| 国产主播在线观看一区二区| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 在线 av 中文字幕| 国产又爽黄色视频| 51午夜福利影视在线观看| 日韩制服骚丝袜av| 亚洲欧美一区二区三区久久| av在线老鸭窝| 高清黄色对白视频在线免费看| 亚洲熟女精品中文字幕| 亚洲伊人久久精品综合| 成年女人毛片免费观看观看9 | 成人手机av| 99香蕉大伊视频| 免费在线观看黄色视频的| 黄片播放在线免费| 久久久久久人人人人人| 国产精品自产拍在线观看55亚洲 | 这个男人来自地球电影免费观看| 欧美精品av麻豆av| 国产无遮挡羞羞视频在线观看| 欧美变态另类bdsm刘玥| 欧美亚洲 丝袜 人妻 在线| 夫妻午夜视频| 亚洲欧美激情在线| 欧美日本中文国产一区发布| 午夜成年电影在线免费观看| 99热网站在线观看| 午夜老司机福利片| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 肉色欧美久久久久久久蜜桃| 日本五十路高清| 12—13女人毛片做爰片一| 国产一区二区 视频在线| 91大片在线观看| 97精品久久久久久久久久精品| 精品人妻1区二区| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 女警被强在线播放| 狠狠狠狠99中文字幕| 国产av一区二区精品久久| 久久国产精品影院| 丰满饥渴人妻一区二区三| a级毛片黄视频| 午夜福利,免费看| 亚洲国产av新网站| netflix在线观看网站| 久久国产精品男人的天堂亚洲| 久久久久国产精品人妻一区二区| 久久精品国产综合久久久| 色视频在线一区二区三区| 中文字幕色久视频| 在线观看免费视频网站a站| av电影中文网址| 热re99久久国产66热| 大码成人一级视频| 久久久精品94久久精品| 丰满人妻熟妇乱又伦精品不卡| 久久影院123| 99久久综合免费| 日本av免费视频播放| 日本一区二区免费在线视频| 18禁黄网站禁片午夜丰满| av片东京热男人的天堂| 爱豆传媒免费全集在线观看| 国产精品 欧美亚洲| 秋霞在线观看毛片| 天堂俺去俺来也www色官网| 丰满饥渴人妻一区二区三| 久久久久久久国产电影| 午夜视频精品福利| av线在线观看网站| av欧美777| 亚洲中文av在线| 精品少妇黑人巨大在线播放| 欧美老熟妇乱子伦牲交| 欧美精品高潮呻吟av久久| 欧美亚洲日本最大视频资源| 久久精品国产a三级三级三级| 韩国精品一区二区三区| 亚洲 国产 在线| 亚洲全国av大片| 涩涩av久久男人的天堂| 国产日韩欧美在线精品| 国产精品久久久久久精品古装| 99久久综合免费| 大型av网站在线播放| 999久久久精品免费观看国产| 久久久久久久精品精品| 久久综合国产亚洲精品| 免费不卡黄色视频| 日韩视频一区二区在线观看| 十八禁人妻一区二区| 午夜福利免费观看在线| 欧美老熟妇乱子伦牲交| 国产野战对白在线观看| 亚洲精品中文字幕在线视频| 国产亚洲av高清不卡| www日本在线高清视频| 午夜日韩欧美国产| 国产99久久九九免费精品| 精品欧美一区二区三区在线| 精品久久久久久电影网| 免费女性裸体啪啪无遮挡网站| 久久久久精品人妻al黑| 亚洲国产av新网站| a 毛片基地| 久久精品aⅴ一区二区三区四区| 一级,二级,三级黄色视频| 一级毛片电影观看| 黄色 视频免费看| 日韩,欧美,国产一区二区三区| 久9热在线精品视频| 一区二区三区四区激情视频| av在线app专区| 少妇裸体淫交视频免费看高清 | 精品福利观看| 亚洲av美国av| 日韩制服骚丝袜av| av不卡在线播放| 极品少妇高潮喷水抽搐| 在线看a的网站| 久久久精品94久久精品| 悠悠久久av| 男女床上黄色一级片免费看| 超色免费av| 黄色a级毛片大全视频| 美女视频免费永久观看网站| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 欧美 亚洲 国产 日韩一| 美女午夜性视频免费| 精品福利永久在线观看| 久久青草综合色| 91字幕亚洲| 久久中文字幕一级| 多毛熟女@视频| 精品国产超薄肉色丝袜足j| 老汉色av国产亚洲站长工具| 国产一区二区在线观看av| 91精品国产国语对白视频| 精品亚洲乱码少妇综合久久| 中亚洲国语对白在线视频| 制服人妻中文乱码| 乱人伦中国视频| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 午夜视频精品福利| 交换朋友夫妻互换小说| 脱女人内裤的视频| 日本欧美视频一区| 国产av又大| 午夜福利视频精品| 日韩精品免费视频一区二区三区| 超色免费av| 在线观看免费高清a一片| 精品亚洲成a人片在线观看| 欧美在线一区亚洲| 啦啦啦在线免费观看视频4| 97在线人人人人妻| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频| 欧美av亚洲av综合av国产av| 亚洲色图 男人天堂 中文字幕| 国产免费视频播放在线视频| 一个人免费看片子| 丝袜人妻中文字幕| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 久久久久久久久久久久大奶| 久久 成人 亚洲| 国产一区有黄有色的免费视频| 免费看十八禁软件| 一本综合久久免费| 欧美变态另类bdsm刘玥| 亚洲男人天堂网一区| 久久久久视频综合| 亚洲精品中文字幕一二三四区 | 高清在线国产一区| 久久久国产一区二区| 精品熟女少妇八av免费久了| 最新的欧美精品一区二区| 免费看十八禁软件| 国产精品成人在线| 成人黄色视频免费在线看| 99精品久久久久人妻精品| 国产av又大| 麻豆国产av国片精品| 欧美亚洲 丝袜 人妻 在线| 国产在线免费精品| 亚洲第一青青草原| 精品一区二区三卡| 老司机靠b影院| 亚洲精品国产av成人精品| 国精品久久久久久国模美| 精品少妇内射三级| 人妻人人澡人人爽人人| 高清欧美精品videossex| 一二三四在线观看免费中文在| 天堂俺去俺来也www色官网| 天天躁夜夜躁狠狠躁躁| 日本欧美视频一区| 18禁国产床啪视频网站| 一个人免费在线观看的高清视频 | 国产麻豆69| av天堂久久9| 最黄视频免费看| 真人做人爱边吃奶动态| 视频区图区小说| 熟女少妇亚洲综合色aaa.| 欧美日韩成人在线一区二区| 欧美激情高清一区二区三区| 岛国毛片在线播放| 久久久国产精品麻豆| 91av网站免费观看| 黄色片一级片一级黄色片| 国产精品秋霞免费鲁丝片| 国产97色在线日韩免费| 黄色怎么调成土黄色| 国产亚洲av高清不卡| 国产精品久久久人人做人人爽| 捣出白浆h1v1| 免费高清在线观看视频在线观看| 精品欧美一区二区三区在线| 真人做人爱边吃奶动态| 精品福利观看| 久久久国产欧美日韩av| 亚洲中文av在线| 中国美女看黄片| 午夜激情久久久久久久| 中文字幕av电影在线播放| 亚洲国产成人一精品久久久| kizo精华| 中文字幕人妻丝袜制服| 一区在线观看完整版| 精品一区在线观看国产| 欧美激情久久久久久爽电影 | 国产野战对白在线观看| 人人妻人人澡人人看| 老司机在亚洲福利影院| 一个人免费在线观看的高清视频 | 日韩 欧美 亚洲 中文字幕| 18禁黄网站禁片午夜丰满| 国产亚洲av高清不卡| 成年美女黄网站色视频大全免费| 丝袜美足系列| 成年人免费黄色播放视频| www.熟女人妻精品国产| 日韩视频在线欧美| 精品熟女少妇八av免费久了| 亚洲精品成人av观看孕妇| 国产精品av久久久久免费| 亚洲情色 制服丝袜| 久久国产精品男人的天堂亚洲| 国产国语露脸激情在线看| 亚洲欧美日韩另类电影网站| 欧美xxⅹ黑人| 亚洲欧美一区二区三区久久| 99国产精品免费福利视频| 国产激情久久老熟女| 人人妻,人人澡人人爽秒播| 欧美亚洲 丝袜 人妻 在线| 亚洲va日本ⅴa欧美va伊人久久 | 久久精品成人免费网站| 欧美97在线视频| 亚洲九九香蕉| 中亚洲国语对白在线视频| 嫁个100分男人电影在线观看| 亚洲欧洲日产国产| 一本大道久久a久久精品| 欧美成人午夜精品| 亚洲中文字幕日韩| 99精品欧美一区二区三区四区| 精品人妻熟女毛片av久久网站| 亚洲欧美一区二区三区黑人| 国产极品粉嫩免费观看在线| 真人做人爱边吃奶动态| 黄色 视频免费看| 亚洲精品粉嫩美女一区| 国产熟女午夜一区二区三区| 精品高清国产在线一区| 国产亚洲av片在线观看秒播厂| 肉色欧美久久久久久久蜜桃| 午夜久久久在线观看| 国产精品免费视频内射| 国产在线视频一区二区| 中文字幕精品免费在线观看视频| 亚洲情色 制服丝袜| 欧美亚洲日本最大视频资源| 色94色欧美一区二区| 日本一区二区免费在线视频| 久久久久久久精品精品| 蜜桃国产av成人99| 女人久久www免费人成看片| 女警被强在线播放| 人妻人人澡人人爽人人| 亚洲伊人色综图| 精品免费久久久久久久清纯 | 巨乳人妻的诱惑在线观看| 日韩一卡2卡3卡4卡2021年| 999精品在线视频| 丝袜人妻中文字幕| 亚洲欧美激情在线| 又大又爽又粗| 97精品久久久久久久久久精品| 最近最新中文字幕大全免费视频| 免费日韩欧美在线观看| 中文字幕人妻丝袜制服| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产精品999| 国产激情久久老熟女| 午夜影院在线不卡| 亚洲成国产人片在线观看| 嫩草影视91久久| 精品熟女少妇八av免费久了| 久久久久久久国产电影| 99久久国产精品久久久| 亚洲精品中文字幕在线视频| 国产熟女午夜一区二区三区| 在线亚洲精品国产二区图片欧美| 18禁黄网站禁片午夜丰满| 亚洲国产欧美日韩在线播放| 最黄视频免费看| 每晚都被弄得嗷嗷叫到高潮| 男人操女人黄网站| 国产成人欧美| 国产黄色免费在线视频| 天堂俺去俺来也www色官网| 亚洲一区二区三区欧美精品| a级毛片黄视频| 免费观看人在逋| 精品卡一卡二卡四卡免费| 日韩三级视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 久久影院123| 国产不卡av网站在线观看| av网站在线播放免费| 在线观看一区二区三区激情| 少妇的丰满在线观看| 乱人伦中国视频| 国产精品秋霞免费鲁丝片| 成年av动漫网址| 欧美人与性动交α欧美精品济南到| 久久久久久免费高清国产稀缺| 亚洲人成电影观看| 亚洲精品国产区一区二| 亚洲熟女精品中文字幕| 18禁观看日本| 99热网站在线观看| 国产成人系列免费观看| 成人国产av品久久久| 波多野结衣av一区二区av| 日本vs欧美在线观看视频| 乱人伦中国视频| 老司机午夜十八禁免费视频| 人人妻,人人澡人人爽秒播| 久久这里只有精品19| 午夜精品国产一区二区电影| 亚洲国产av影院在线观看| 999久久久国产精品视频| 日韩有码中文字幕| 欧美亚洲日本最大视频资源| 啪啪无遮挡十八禁网站| 少妇人妻久久综合中文| 香蕉国产在线看| 黄色怎么调成土黄色| 美女中出高潮动态图| 一级,二级,三级黄色视频| 精品国产乱子伦一区二区三区 | 国产一区二区激情短视频 | 一级黄色大片毛片| 老汉色av国产亚洲站长工具| 国产黄频视频在线观看| 精品国产国语对白av| 国产精品一区二区精品视频观看| av在线播放精品| 少妇被粗大的猛进出69影院| av片东京热男人的天堂| 一区二区三区激情视频| 久久人妻福利社区极品人妻图片| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 国产精品秋霞免费鲁丝片| 黄色视频在线播放观看不卡| 啦啦啦 在线观看视频| 在线观看免费午夜福利视频| 婷婷色av中文字幕| 成年美女黄网站色视频大全免费| 脱女人内裤的视频| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 国产一级毛片在线| 嫁个100分男人电影在线观看| 9191精品国产免费久久| 一区福利在线观看| 日韩有码中文字幕| 老司机在亚洲福利影院| 成年人免费黄色播放视频| 黑人猛操日本美女一级片| 亚洲 国产 在线| 久久久久久久久久久久大奶| 男男h啪啪无遮挡| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| 亚洲精品av麻豆狂野| 黄片小视频在线播放| 久久人妻福利社区极品人妻图片| 免费观看a级毛片全部| 超碰成人久久| 欧美日韩亚洲综合一区二区三区_| 丰满少妇做爰视频| 日韩熟女老妇一区二区性免费视频| 黄色怎么调成土黄色| 久久亚洲精品不卡| 人妻久久中文字幕网| 午夜福利视频在线观看免费| 精品国产乱码久久久久久小说| 秋霞在线观看毛片| 中亚洲国语对白在线视频| 欧美日韩精品网址| 69av精品久久久久久 | 欧美大码av| 热99久久久久精品小说推荐| 一级片'在线观看视频| 精品一区二区三卡| 狂野欧美激情性xxxx| 欧美av亚洲av综合av国产av| 欧美黄色淫秽网站| 丁香六月天网| 久久久精品国产亚洲av高清涩受| 老司机靠b影院| 大香蕉久久成人网| 一级a爱视频在线免费观看| 韩国精品一区二区三区| 波多野结衣av一区二区av| 午夜老司机福利片| 久久久久久久大尺度免费视频| 男女边摸边吃奶| 欧美国产精品va在线观看不卡| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久5区| 亚洲 欧美一区二区三区| 精品人妻1区二区| 午夜福利乱码中文字幕| 国产主播在线观看一区二区| av片东京热男人的天堂| 极品少妇高潮喷水抽搐| 操美女的视频在线观看| av在线app专区| 精品熟女少妇八av免费久了| 丝袜在线中文字幕| 日韩人妻精品一区2区三区| 五月天丁香电影| 亚洲国产精品一区二区三区在线| 亚洲自偷自拍图片 自拍| 国产精品一区二区在线观看99| 精品一区在线观看国产| 80岁老熟妇乱子伦牲交| 成人国产av品久久久| 久久这里只有精品19| 国产精品久久久久成人av| 久久av网站| 亚洲国产欧美在线一区| 丁香六月天网| 一进一出抽搐动态| 国产在视频线精品| 亚洲精品一二三| 国产高清videossex| 免费少妇av软件| 日日摸夜夜添夜夜添小说| 激情视频va一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人一区二区三| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 18禁国产床啪视频网站| 欧美另类亚洲清纯唯美| 精品亚洲成a人片在线观看| av片东京热男人的天堂| 国产日韩欧美亚洲二区| 国产成人一区二区三区免费视频网站| 18禁观看日本| bbb黄色大片| www.自偷自拍.com| 热99re8久久精品国产| 男人舔女人的私密视频| 波多野结衣一区麻豆| 另类精品久久| 成人手机av| 欧美亚洲日本最大视频资源| 国产欧美日韩一区二区精品| 国产成人欧美| 最黄视频免费看| 好男人电影高清在线观看| 69精品国产乱码久久久| 婷婷色av中文字幕| 一区在线观看完整版| 亚洲专区中文字幕在线| 麻豆国产av国片精品| 涩涩av久久男人的天堂| 国产日韩一区二区三区精品不卡| 俄罗斯特黄特色一大片| 亚洲欧美日韩另类电影网站| 老司机午夜十八禁免费视频| 蜜桃在线观看..| 午夜91福利影院| 国产精品 国内视频| 好男人电影高清在线观看| 免费久久久久久久精品成人欧美视频| 亚洲精品粉嫩美女一区| 下体分泌物呈黄色| 美女主播在线视频| 亚洲精品久久成人aⅴ小说| 久久综合国产亚洲精品| 91av网站免费观看| 不卡一级毛片| 大码成人一级视频| 久久久水蜜桃国产精品网| 免费在线观看视频国产中文字幕亚洲 | 一级毛片女人18水好多| 性色av乱码一区二区三区2| 亚洲精品第二区| 日本撒尿小便嘘嘘汇集6| 欧美变态另类bdsm刘玥| 亚洲欧美一区二区三区久久| 大香蕉久久成人网| 亚洲欧美色中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美免费精品| 老汉色∧v一级毛片| 午夜激情久久久久久久| 国产精品九九99| 一边摸一边抽搐一进一出视频| 操出白浆在线播放| 一区福利在线观看| 精品卡一卡二卡四卡免费| 亚洲精品久久久久久婷婷小说| 国产野战对白在线观看| 人人妻人人添人人爽欧美一区卜| 国产福利在线免费观看视频| 亚洲全国av大片| 夫妻午夜视频| 久久精品人人爽人人爽视色| 国产精品一区二区精品视频观看| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产亚洲av香蕉五月 | 久久久久久久久免费视频了| 黄频高清免费视频| 亚洲国产成人一精品久久久| 国产一区二区在线观看av| 午夜影院在线不卡| 伦理电影免费视频| 成年人免费黄色播放视频| 亚洲欧美清纯卡通| 久久精品成人免费网站| 天天躁夜夜躁狠狠躁躁| 午夜激情久久久久久久| 高清视频免费观看一区二区| 99国产精品99久久久久| 51午夜福利影视在线观看| 日韩中文字幕欧美一区二区| 精品国产一区二区久久| 久久人人97超碰香蕉20202| 亚洲色图 男人天堂 中文字幕| 亚洲一区中文字幕在线| 一级片免费观看大全| www日本在线高清视频| 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 亚洲国产av影院在线观看| 丁香六月欧美| 日本a在线网址| 黄色视频在线播放观看不卡| 亚洲av国产av综合av卡| 欧美日韩黄片免| 黄色片一级片一级黄色片| videos熟女内射| svipshipincom国产片| 国产免费av片在线观看野外av| 国产精品 欧美亚洲| 欧美午夜高清在线| 日日夜夜操网爽| 欧美激情高清一区二区三区| 欧美日韩精品网址| 国产色视频综合| av线在线观看网站| 亚洲avbb在线观看| 国产欧美日韩精品亚洲av| 又黄又粗又硬又大视频| 久久精品成人免费网站| 18禁观看日本| 久久久久国产一级毛片高清牌| 在线观看免费视频网站a站| 久久久久久久久久久久大奶|