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

    靜電微泵的3D流固耦合動態(tài)特性分析

    2012-12-03 14:50:30董金新林謝昭
    中國機(jī)械工程 2012年5期
    關(guān)鍵詞:微泵瞬態(tài)流體

    董金新 林謝昭 應(yīng) 濟(jì) 李 俊

    1.溫州職業(yè)技術(shù)學(xué)院,溫州,325035 2.浙江大學(xué),杭州,310027

    0 引言

    微流體系統(tǒng)中大都要求有微泵。其中,薄膜往復(fù)振動式無閥微泵是以擴(kuò)散/收縮單元為閥門,通過泵膜的振動驅(qū)動流體,其制造簡單,可以驅(qū)動一些非均相的流體,是眾多研究的焦點。根據(jù)驅(qū)動原理不同,無閥微泵可劃分為電磁、靜電、壓電、形狀記憶合金等多種類型。其中,靜電無閥微泵具有較低功耗,易于與IC工藝集成兼容等優(yōu)點,而得到廣泛關(guān)注。

    由于微泵的尺度很小,薄膜往復(fù)振動式無閥微泵的絕大部分物理量難以通過實驗測量,特別是瞬態(tài)量。因此,在微泵的研制階段,利用數(shù)值方法對動態(tài)性能進(jìn)行仿真、預(yù)測顯得很重要。它不僅可降低費(fèi)用,而且能更好地了解微泵的工作原理以及可能的潛在缺陷。但因微泵中柔性泵膜、電驅(qū)動和流體相互耦合,加之流體在擴(kuò)散/收縮單元中的不同流動方向表現(xiàn)出的不同壓力特性,增加了對整個泵的仿真難度和復(fù)雜性。早期的研究采用低階集總參數(shù)模型[1-2]和等效電網(wǎng)絡(luò)模型[3]等簡化方法,忽略了參數(shù)的空間分布特性。目前大部分的研究沒有考慮泵的流固耦合效應(yīng),也無法建立電驅(qū)動信號量與流體動力學(xué)量(比如流量)之間的關(guān)系。或是研究驅(qū)動器的動態(tài)特性[4-7];或是簡化流體場,根據(jù)擴(kuò)散/收縮單元的壓力損失系數(shù),將泵膜的運(yùn)動與泵腔內(nèi)外壓力以及擴(kuò)散/收縮單元流量聯(lián)系起來,研究泵膜的動態(tài)特性[8-11];或是將驅(qū)動器的位移輸出,作為流體場的移動壁面問題來研究泵流體的動力學(xué)特性[12-14]。實際上,泵膜的動態(tài)特性與流體黏滯損失的非線性和不穩(wěn)定性等是相互影響的。為了理解整個微泵的動態(tài)特性,模型必須能有效地描述這些耦合效應(yīng)[15-16]。本文以靜電無閥微泵為對象,建立靜電-結(jié)構(gòu)-流體全耦合的3維模型,利用數(shù)值方法,仿真并分析各個變量的非線性動態(tài)特性。

    1 模型的建立

    1.1 控制方程、邊界條件以及初始條件

    圖1所示為靜電無閥微泵的結(jié)構(gòu)示意圖。對于泵送流體,采用任意拉格朗日-歐拉(ALE)描述比較方便。某一時刻,在參考系中^x位置流動的流體質(zhì)點速度u的時間導(dǎo)數(shù)可以表示為

    圖1 靜電無閥微泵的結(jié)構(gòu)示意圖

    流體域Ωf的不可壓縮流動動量守恒定律以及連續(xù)方程可以表示為

    式中,ρf、f、σ分別為流體密度、體積力向量和Cauchy應(yīng)力張量;[0,T]為所考察的時間間隔。

    牛頓流體的本構(gòu)方程為

    在泵膜流固耦合界面Γf-s處,邊界條件如下:

    式(5)為流固界面Γf-s處的無滑移條件;在流固界面Γf-s上,參考系的速度^v及其位置^x需要滿足連續(xù)條件式(6),以保證參考系的邊界準(zhǔn)確地表達(dá)結(jié)構(gòu)的當(dāng)前構(gòu)型。沿著界面的應(yīng)力平衡由式(7)表達(dá)。

    對于泵膜,可以利用標(biāo)準(zhǔn)的Lagrangian描述建立其運(yùn)動學(xué)方程,即

    式中,ρs為固體的密度。

    泵膜流固耦合邊界Γf-s滿足:

    邊界條件式(9)、式(10)分別對應(yīng)于式(5)和式(7)。

    如果沒有包含自由電荷,則描述靜電場的Poisson方程可以表示為

    式中,φ為靜電勢;ε為介電常數(shù)張量;E為靜電場電場強(qiáng)度。

    Maxwell靜電應(yīng)力張量TE由下式計算:

    式中,F(xiàn)E為靜電載荷。

    各個場變量的初始條件都設(shè)置為零。為了計算易于收斂,將一個很小的初值輸入到模型中。

    1.2 模型參數(shù)及材料屬性

    圖1所示靜電微泵的泵膜厚度為10μm,泵腔深度為100μm,泵膜半徑為2000μm;擴(kuò)散/收縮單元的幾何結(jié)構(gòu)相同,長度為1000μm,最小端面寬度尺寸為80μm,擴(kuò)張角為10°;連接的進(jìn)出口管道長度皆為2000μm,寬度為200μm。流體介質(zhì)為去離子水,物性參數(shù)如表1所示。

    表1 與微泵仿真相關(guān)的材料屬性參數(shù)

    假設(shè)層流流動,考慮對稱性,在CFD-ACE中建立微泵的一半模型。利用結(jié)構(gòu)化網(wǎng)格劃分流體,單元數(shù)為52 880個,泵膜劃分為4974個殼單元。使用基于壓力修正的SIMPLE-C 算法對流動以及能量進(jìn)行順序積分,獲得流場的解。流體域的對流和擴(kuò)散采用一階迎風(fēng)格式,時間積分采用了Crank-Nicolson格式。使用基于有限體積法的計算格式求解靜電場,這樣可以處理不同介質(zhì)介電常數(shù)問題。利用一階殼單元來表達(dá)泵膜,并進(jìn)行大變形和接觸等幾何非線性分析??紤]到靜電、結(jié)構(gòu)以及流動耦合求解的需要,利用網(wǎng)格重劃分工具,使用預(yù)測-校正方法對局部變形進(jìn)行連續(xù)修正,將流動速度和結(jié)構(gòu)速度聯(lián)系起來,實現(xiàn)流固耦合分析。流體和結(jié)構(gòu)求解器的耦合迭代頻率設(shè)置為fc=1,使得流體場和結(jié)構(gòu)場信息能及時得到交換和更新。當(dāng)量綱一殘差小于10-4時,認(rèn)為計算是收斂的。施加的周期性驅(qū)動電壓V=200(1-cos(20πt))(泵腔底電壓為零),在Pentium-3.0GHz的PC機(jī)中,完成3.5個周期的求解,共350個子步數(shù),求解時間約為59h。

    2 結(jié)果及討論

    2.1 泵膜的動力學(xué)特性

    圖2所示為計算得到的泵膜中心點垂直方向上的位移隨時間變化曲線,圖中最大變形量為28.433μm??紤]流固耦合效應(yīng)后,在流體阻尼作用下,泵膜無法回復(fù)至初始零位置,0.653μm為新的平衡位置。穩(wěn)態(tài)情況下,泵膜位移的變化滯后于輸入電壓信號的變化,一個周期的滯后時間約為0.01s。

    圖2 泵膜中心垂直方向上的位移隨時間變化曲線

    2.2 流場的瞬態(tài)特性

    實際上,泵流體動態(tài)特性更為設(shè)計者所關(guān)注,為了更好地詮釋泵膜振動與流體特性之間的耦合關(guān)系,用圖3所示的泵膜流固耦合界面上中心點處的垂直方向位移與泵腔內(nèi)流體壓力來顯示其對應(yīng)關(guān)系,從圖中可以明顯地看出,流體壓力的變化與泵膜的位移存在著密切關(guān)系。在靜電力的作用下,泵膜下拉變形的速度大于回彈的速度,所以在流體壓出泵腔的過程中其壓力變化較大,最高壓力為945.1Pa;而在流體吸入泵腔過程中,其壓力變化相對較小,最大負(fù)壓為-558.0Pa;整個泵送過程中,流體壓力呈現(xiàn)非對稱性變化。圖3還表明,泵初始工作點從零點開始,經(jīng)過第一個周期的瞬態(tài)過程后,第二個周期的軌跡與第三個周期的軌跡已經(jīng)重合,表明泵已經(jīng)進(jìn)入穩(wěn)定工作狀態(tài)。流體壓力的波動周期同樣滯后于驅(qū)動電壓的變化(約0.01s)。

    圖3 泵膜中心點處垂直方向上位移與流體壓力的關(guān)系

    圖4所示為進(jìn)出口單元處的瞬態(tài)體積流量變化曲線,從圖中可以看出,在泵膜下拉(往固定電極)運(yùn)動階段(泵送階段),出口單元的瞬態(tài)流量大于進(jìn)口單元的瞬態(tài)流量;而在泵膜回彈階段(泵吸入階段),進(jìn)口單元的流量則大于出口單元的流量。最大的瞬態(tài)流量差發(fā)生在泵腔內(nèi)壓力達(dá)到極大值時刻,而在泵腔內(nèi)壓力與進(jìn)出口壓力相差不大的時間段,進(jìn)出口單元的流量相差不大。表明要使擴(kuò)散/收縮單元發(fā)揮“整流”作用,其兩端的壓力差必須達(dá)到一定值,并且壓力差要盡可能地大。對穩(wěn)態(tài)情況下的泵凈流量進(jìn)行積分計算,得到的凈輸出流量為10.764nL/s。

    圖4 進(jìn)出口瞬態(tài)流量隨時間的變化

    圖5所示為微泵橫截面方向上流體壓力的分布情況,由圖可知泵腔內(nèi)的壓力幾乎分布一致,說明在平面布置的無閥微泵中,采用Reynolds方程來描述泵腔內(nèi)流體的瞬態(tài)動力學(xué)特性具有一定的可行性。另外,微泵所連接的進(jìn)出口管道對微泵的工作性能也有影響。

    圖5 不同時刻下的微泵橫截面方向上的壓力分布情況

    考慮了流固耦合特性后,模型不僅能夠仿真泵流場的動態(tài)特性,還能夠仿真泵膜的應(yīng)力特性變化。圖6所示為t=28ms時刻的泵腔內(nèi)流體壓力分布與泵膜應(yīng)力分布情況,以及這一時刻的擴(kuò)散/收縮單元內(nèi)流體的速度場分布。此時,泵腔內(nèi)流體壓力達(dá)到最大值945.3Pa,泵膜中心撓度為13.20μm;泵 膜 邊 緣 處 VonMises 應(yīng) 力 為28.39MPa。微泵腔內(nèi)的流速較小,擴(kuò)張管/收縮管頸部的平均流速達(dá)到最大,計算得到的平均流速為0.0867m/s,雷諾數(shù)為7.71,遠(yuǎn)小于宏觀條件下通常認(rèn)定的臨界雷諾數(shù)(2000左右),這也證實了前面假設(shè)層流流動模式的正確性。另外,圖6b和圖6c表明,收縮/擴(kuò)散管內(nèi)未出現(xiàn)流體與固體壁面分離流動現(xiàn)象。

    圖6 t=28ms時刻的泵膜應(yīng)力和微泵的流體場瞬態(tài)結(jié)果

    3 結(jié)論

    (1)擴(kuò)散/收縮單元的最大流量差都發(fā)生在泵腔內(nèi)壓力與進(jìn)出口壓力相差最大的時刻,表明要使得擴(kuò)散單元發(fā)揮“整流”作用,其單元兩端的壓力差必須達(dá)到一定值。

    (2)在仿真的每個時刻,泵腔內(nèi)流體壓力分布幾乎一致,表明用Reynolds方程來描述泵腔內(nèi)流體動力學(xué)特性具有可行性。

    (3)仿真得到的最大雷諾數(shù)遠(yuǎn)小于通常認(rèn)定的臨界雷諾數(shù),表明微泵流體具有層流特性。

    (4)泵腔內(nèi)流體動態(tài)特性與泵膜的運(yùn)動存在著密切關(guān)系,如果忽略了流固耦合效應(yīng),簡單地將微泵的驅(qū)動器與腔內(nèi)的流場分離處理,必然歪曲了泵流體場的動態(tài)特性。泵膜變形、流體壓力和進(jìn)出口流速等響應(yīng)與輸入電壓信號之間存在著滯后現(xiàn)象。對微泵的3D全耦合仿真,能夠獲得驅(qū)動電信號與流體流量、泵送壓力等輸出變量之間的直接關(guān)系,有利于從整體意義上實現(xiàn)微泵的優(yōu)化設(shè)計。

    [1]Morris C J,F(xiàn)orster F K.Low-order Modeling of Resonance for Fixed-valve Micropumps Based on First Principles[J].Journal Microelectromech.Syst.,2003,12:325-334.

    [2]Olsson A,Stemme G,Stemme E.A Numerical Design Study of the Valveless Diffuser Pump Using a Lumped-mass Model[J].Journal Micromech.Microeng.,1999,9:34-44.

    [3]Bourouina T,Grandchamp J P.Modeling Micropumps with Electrical Equivalent Networks[J].Journal Micromech.Microeng.,1996,6:398-404.

    [4]Lin X Z,Ying J,Chen Zichen.Research on Terminal Behavior of Electrostatically Actuated Micropump Membrane Based on Modal Analysis[C]//Wei Y L,Chong K T,Takahashi T,et al.ICMIT 2005:Mechatronics,MEMS and Smart Materials,SPIE.Bellingham,2005,6040:P60400U.

    [5]Faris W F,Abdel-Rahman E M,Nayfeh A H.Mechanical Behavior of an Electrostatically Actuated Micropump[C]//Agnes G,Chamis C,Noor A,et al.43rd AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference.Denver,2002:AIAA2002-1303.

    [6]Morris C J,F(xiàn)orster F K .Optimization of a Circular Piezoelectric Bimorph for a Micropump Driver[J].Journal Micromech.Microeng.,2000,10:65-459.

    [7]王立文,高殿榮,楊林杰,等.壓電驅(qū)動微泵泵膜振動有限元分析[J].機(jī)械工程學(xué)報,2006,42(4):230-235.

    [8]Ullman A,F(xiàn)ono I.The Piezoelectric Valve-less Pump—improved Dynamic Model[J].Journal Microelectromech.Syst.,2002,11:655-664.

    [9]Pan L S,Ng T Y,Wu X H,et al.Analysis of Valveless Micropumps with Inertial Effects[J].Journal Micromech.Microeng.,2003,13:390-399.

    [10]MacHauf A,Nemirovsky Y,Dinnar U.A Membrane Micropump Electrostatically Actuated Across the Working Fluid[J].Journal Micromech.Microeng.,2005,15(12):2309.

    [11]白蘭,馮志慶,吳一輝.無閥微泵動態(tài)特性的固液耦合分析[J].機(jī)械工程學(xué)報,2008,44(7):69-74.

    [12]謝海波,傅新,楊華勇.微型無閥泵流動特征仿真與試驗研究[J].機(jī)械工程學(xué)報,2002,38(7):54-57.

    [13]楊愷祥.壓電無閥式微泵浦制造與測量分析[D].臺灣云林縣:臺灣云林科技大學(xué),2004.

    [14]Tsui Y Y,Lu S L.Evaluation of the Performance of a Valveless Micropump by CFD and Lumpedsystem Analyses[J].Sensors and Actuators A:Physical,2008,148(1):138.

    [15]Yao Q,Xu D,Pan L S,et al.CFD Simulation of Flows in Valveless Micropumps[J].Eng.Appl.Comput.Fluid Mech.,2007,1(3):181-188.

    [16]林謝昭,應(yīng)濟(jì),陳子辰.靜電無閥微泵兩種仿真模型的比較研究[J].哈爾濱工業(yè)大學(xué)學(xué)報,2010,42(12):1943-1947.

    猜你喜歡
    微泵瞬態(tài)流體
    流體壓強(qiáng)知多少
    基于微流控芯片的五腔室壓電蠕動微泵的設(shè)計
    分娩鎮(zhèn)痛聯(lián)合微泵靜注縮宮素的臨床效果觀察
    高壓感應(yīng)電動機(jī)斷電重啟時的瞬態(tài)仿真
    山雨欲來風(fēng)滿樓之流體壓強(qiáng)與流速
    等效流體體積模量直接反演的流體識別方法
    一種微泵懸掛裝置的設(shè)計和應(yīng)用
    利用體塊PZT制備膜片式壓電微泵
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    欧美一级a爱片免费观看看| 热99在线观看视频| 久久久精品94久久精品| 99久久精品国产国产毛片| 亚洲人成网站在线播| 久久久久免费精品人妻一区二区| 国产三级在线视频| 亚洲无线在线观看| 简卡轻食公司| 日本三级黄在线观看| 免费看光身美女| 欧美一区二区精品小视频在线| 女人被狂操c到高潮| 久久这里有精品视频免费| 女同久久另类99精品国产91| 亚州av有码| 欧美成人一区二区免费高清观看| 午夜久久久久精精品| 黄色一级大片看看| 91aial.com中文字幕在线观看| 国产一区二区三区在线臀色熟女| 久久精品国产99精品国产亚洲性色| kizo精华| 久久精品人妻少妇| 亚洲欧洲日产国产| 久久精品国产亚洲网站| 国产探花在线观看一区二区| 99久久久亚洲精品蜜臀av| 国产精品国产三级国产av玫瑰| 日韩成人av中文字幕在线观看| 国产视频内射| 日韩成人伦理影院| 国产精品久久久久久亚洲av鲁大| 一级毛片久久久久久久久女| 欧美极品一区二区三区四区| 波野结衣二区三区在线| 2021天堂中文幕一二区在线观| 97在线视频观看| 亚洲欧美中文字幕日韩二区| 亚洲精华国产精华液的使用体验 | 久久久久性生活片| 99久久精品一区二区三区| 日本黄大片高清| 日韩三级伦理在线观看| 禁无遮挡网站| 亚洲丝袜综合中文字幕| 日日啪夜夜撸| 岛国在线免费视频观看| 精品日产1卡2卡| 国产精品人妻久久久久久| 又粗又爽又猛毛片免费看| 久久人人爽人人爽人人片va| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 欧美日本亚洲视频在线播放| 一区二区三区四区激情视频 | 国产精品,欧美在线| 亚洲精品国产av成人精品| 免费黄网站久久成人精品| 色综合色国产| 久久亚洲精品不卡| 成人午夜高清在线视频| 九九爱精品视频在线观看| 伦理电影大哥的女人| 一区福利在线观看| av在线播放精品| 亚洲av不卡在线观看| 免费看光身美女| 精品久久久久久久久亚洲| 国产黄色小视频在线观看| 久久精品国产99精品国产亚洲性色| 亚洲欧美成人精品一区二区| 国产不卡一卡二| 日韩av在线大香蕉| 别揉我奶头 嗯啊视频| 三级男女做爰猛烈吃奶摸视频| 婷婷六月久久综合丁香| 国产精华一区二区三区| 联通29元200g的流量卡| 精品人妻偷拍中文字幕| 亚洲精品国产成人久久av| 一级黄片播放器| 91精品一卡2卡3卡4卡| 久久99精品国语久久久| 淫秽高清视频在线观看| 国产亚洲欧美98| 国产免费男女视频| 欧美bdsm另类| 亚洲四区av| 成人综合一区亚洲| 天天躁夜夜躁狠狠久久av| 亚洲精品国产成人久久av| 久久精品国产亚洲网站| 一个人看视频在线观看www免费| 日本色播在线视频| 日本一二三区视频观看| 国产综合懂色| 寂寞人妻少妇视频99o| 国产黄a三级三级三级人| 最近的中文字幕免费完整| 国产成人午夜福利电影在线观看| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 夜夜夜夜夜久久久久| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 一本久久精品| 大又大粗又爽又黄少妇毛片口| 午夜福利高清视频| 国产视频首页在线观看| 国产人妻一区二区三区在| 午夜免费男女啪啪视频观看| 久久精品国产自在天天线| 好男人视频免费观看在线| av在线老鸭窝| 日产精品乱码卡一卡2卡三| 观看美女的网站| 婷婷色av中文字幕| 国产精品一区二区在线观看99 | 夜夜看夜夜爽夜夜摸| 久久久久久久久大av| 久久午夜亚洲精品久久| 国产在线男女| 一区福利在线观看| 99久久精品国产国产毛片| a级一级毛片免费在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲精品国产av成人精品| 联通29元200g的流量卡| 亚洲真实伦在线观看| 日本av手机在线免费观看| 亚洲不卡免费看| 亚洲国产高清在线一区二区三| 哪里可以看免费的av片| 国产精品美女特级片免费视频播放器| 伦理电影大哥的女人| 美女黄网站色视频| 久久精品夜色国产| 国内精品一区二区在线观看| 九色成人免费人妻av| 国产精品一二三区在线看| 黄色一级大片看看| 精品午夜福利在线看| 人妻少妇偷人精品九色| 欧美+日韩+精品| 男女做爰动态图高潮gif福利片| 99riav亚洲国产免费| 99精品在免费线老司机午夜| 黄片wwwwww| 午夜免费激情av| 亚洲欧美日韩高清在线视频| 色噜噜av男人的天堂激情| 禁无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 久久99热这里只有精品18| 亚洲av免费高清在线观看| 小蜜桃在线观看免费完整版高清| 欧美一区二区精品小视频在线| 国产精品人妻久久久影院| 日韩成人av中文字幕在线观看| 美女高潮的动态| 伦理电影大哥的女人| 亚洲av免费在线观看| 午夜激情福利司机影院| 成人亚洲欧美一区二区av| 久久人人爽人人爽人人片va| 内地一区二区视频在线| 国产黄片美女视频| eeuss影院久久| av视频在线观看入口| 久久精品国产亚洲网站| 小蜜桃在线观看免费完整版高清| 日韩人妻高清精品专区| 综合色丁香网| 99热这里只有是精品50| 亚洲真实伦在线观看| 亚洲高清免费不卡视频| 看片在线看免费视频| 国产不卡一卡二| 一本一本综合久久| 久久久久免费精品人妻一区二区| 亚洲国产欧美在线一区| 夜夜夜夜夜久久久久| 小说图片视频综合网站| 综合色av麻豆| 国产又黄又爽又无遮挡在线| 国产一区二区三区av在线 | 久久精品人妻少妇| 免费搜索国产男女视频| 联通29元200g的流量卡| 国产极品天堂在线| 激情 狠狠 欧美| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 国产成人91sexporn| 午夜激情欧美在线| 国产精品.久久久| 国产男人的电影天堂91| 欧美丝袜亚洲另类| 久久久久久久久中文| 2021天堂中文幕一二区在线观| 欧美xxxx黑人xx丫x性爽| 国产在线男女| 深夜精品福利| 国产精品久久久久久久久免| 身体一侧抽搐| 国产精品女同一区二区软件| 乱码一卡2卡4卡精品| 麻豆久久精品国产亚洲av| 午夜亚洲福利在线播放| 久久久精品大字幕| 欧美日韩精品成人综合77777| 久久久色成人| 国产色婷婷99| 热99在线观看视频| 99热6这里只有精品| 婷婷精品国产亚洲av| 综合色av麻豆| 精品国内亚洲2022精品成人| 高清日韩中文字幕在线| 又爽又黄无遮挡网站| 亚洲欧美精品自产自拍| 69av精品久久久久久| 久久久a久久爽久久v久久| 91aial.com中文字幕在线观看| 免费观看的影片在线观看| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| 成人性生交大片免费视频hd| 日本成人三级电影网站| 国产一区二区亚洲精品在线观看| 亚洲图色成人| 免费看a级黄色片| 久久国内精品自在自线图片| 免费无遮挡裸体视频| 色综合站精品国产| 中文字幕人妻熟人妻熟丝袜美| 国产精品综合久久久久久久免费| 一级毛片我不卡| 国产伦在线观看视频一区| 国产精品人妻久久久影院| 亚洲av不卡在线观看| 国产精品综合久久久久久久免费| 国产亚洲精品久久久com| 搞女人的毛片| av在线天堂中文字幕| 久久人妻av系列| 国产 一区精品| 天美传媒精品一区二区| 国产激情偷乱视频一区二区| 亚洲国产欧美在线一区| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 国产亚洲精品久久久com| 亚洲国产精品合色在线| 九九热线精品视视频播放| av黄色大香蕉| 婷婷六月久久综合丁香| 色视频www国产| 麻豆久久精品国产亚洲av| 天天躁日日操中文字幕| 久久久久久久久中文| 99热只有精品国产| 九九爱精品视频在线观看| 只有这里有精品99| 久久久久网色| 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 18+在线观看网站| 精品午夜福利在线看| 免费黄网站久久成人精品| 人人妻人人澡欧美一区二区| 麻豆成人午夜福利视频| 九九在线视频观看精品| 亚洲七黄色美女视频| 看黄色毛片网站| 久久精品国产清高在天天线| 日本在线视频免费播放| 黄片无遮挡物在线观看| 国产黄色视频一区二区在线观看 | 我的女老师完整版在线观看| 久久精品国产亚洲av天美| 亚洲自拍偷在线| 最近2019中文字幕mv第一页| 欧美精品一区二区大全| 日本一二三区视频观看| 免费av不卡在线播放| 免费看美女性在线毛片视频| 中文资源天堂在线| 国产精品日韩av在线免费观看| 日本一本二区三区精品| eeuss影院久久| 国产69精品久久久久777片| 日韩一本色道免费dvd| 禁无遮挡网站| 人体艺术视频欧美日本| 亚洲av熟女| 国产成人午夜福利电影在线观看| 成人特级av手机在线观看| 亚洲在线自拍视频| 久久久久久久久大av| 午夜精品国产一区二区电影 | 好男人在线观看高清免费视频| 免费观看精品视频网站| 午夜a级毛片| 国产日本99.免费观看| 亚洲一区二区三区色噜噜| 精品久久久久久久末码| 国产精品.久久久| 一个人免费在线观看电影| 看黄色毛片网站| 变态另类丝袜制服| 99热这里只有是精品50| 美女黄网站色视频| 中文字幕人妻熟人妻熟丝袜美| 最新中文字幕久久久久| 嘟嘟电影网在线观看| av在线观看视频网站免费| 一进一出抽搐动态| 国产爱豆传媒在线观看| 久久久色成人| 边亲边吃奶的免费视频| 午夜精品在线福利| 你懂的网址亚洲精品在线观看 | 一本精品99久久精品77| 久久99热这里只有精品18| 欧美另类亚洲清纯唯美| 精品人妻视频免费看| av女优亚洲男人天堂| av.在线天堂| 久久6这里有精品| 国产色婷婷99| 国产午夜福利久久久久久| 熟女人妻精品中文字幕| 国产成人精品一,二区 | 免费观看a级毛片全部| 嫩草影院新地址| 亚洲图色成人| 久久亚洲精品不卡| 欧洲精品卡2卡3卡4卡5卡区| 亚洲美女视频黄频| 欧美另类亚洲清纯唯美| eeuss影院久久| 亚洲精品456在线播放app| 中文欧美无线码| 久久人人爽人人片av| 在线观看av片永久免费下载| 熟女人妻精品中文字幕| 国产乱人视频| 国产黄色视频一区二区在线观看 | 久久久精品94久久精品| 午夜视频国产福利| 淫秽高清视频在线观看| 国产精品不卡视频一区二区| 最近视频中文字幕2019在线8| 久久久久久久久中文| 九九热线精品视视频播放| 成人午夜精彩视频在线观看| 久久久久久大精品| 免费av观看视频| 久久久久久大精品| 日韩成人伦理影院| 观看美女的网站| 一个人看视频在线观看www免费| 男女边吃奶边做爰视频| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线 | 中文欧美无线码| 美女cb高潮喷水在线观看| 在线天堂最新版资源| 99久久无色码亚洲精品果冻| 美女脱内裤让男人舔精品视频 | 国产精品一区二区三区四区免费观看| 国产精品久久久久久av不卡| 国产精品国产高清国产av| а√天堂www在线а√下载| 久久鲁丝午夜福利片| 欧美zozozo另类| 国产白丝娇喘喷水9色精品| 夜夜爽天天搞| 亚洲成人av在线免费| 国产精品久久视频播放| 成人亚洲欧美一区二区av| 插阴视频在线观看视频| 成年女人看的毛片在线观看| 久99久视频精品免费| 日韩av不卡免费在线播放| 中文字幕精品亚洲无线码一区| 2022亚洲国产成人精品| 成人鲁丝片一二三区免费| 国产精品无大码| 又爽又黄无遮挡网站| 亚洲成人中文字幕在线播放| 亚洲不卡免费看| 日日撸夜夜添| 啦啦啦韩国在线观看视频| 有码 亚洲区| 看黄色毛片网站| 99九九线精品视频在线观看视频| 不卡视频在线观看欧美| 国产高清有码在线观看视频| 在线a可以看的网站| 亚洲欧美清纯卡通| 不卡一级毛片| 一级黄片播放器| 日本与韩国留学比较| 国产v大片淫在线免费观看| 一区福利在线观看| 久久久a久久爽久久v久久| 亚洲av成人av| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 日韩,欧美,国产一区二区三区 | 日本三级黄在线观看| 国产私拍福利视频在线观看| 一级毛片aaaaaa免费看小| a级毛片a级免费在线| 日韩成人伦理影院| 亚洲七黄色美女视频| 亚洲av二区三区四区| 插阴视频在线观看视频| 国产女主播在线喷水免费视频网站 | 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 插逼视频在线观看| 久99久视频精品免费| 2022亚洲国产成人精品| 午夜福利成人在线免费观看| 丰满的人妻完整版| 中文字幕制服av| 亚洲av男天堂| 天堂av国产一区二区熟女人妻| 青春草亚洲视频在线观看| 人妻制服诱惑在线中文字幕| 久久久久九九精品影院| 在线观看av片永久免费下载| 亚洲精品456在线播放app| 国产精品av视频在线免费观看| 真实男女啪啪啪动态图| 久久久久久国产a免费观看| 国产成人一区二区在线| 性插视频无遮挡在线免费观看| av.在线天堂| 一进一出抽搐动态| 国产精品福利在线免费观看| 校园春色视频在线观看| 少妇丰满av| 欧美另类亚洲清纯唯美| 国产精品人妻久久久影院| 欧美变态另类bdsm刘玥| 91在线精品国自产拍蜜月| 欧美日韩综合久久久久久| 亚洲av免费高清在线观看| 亚洲精品亚洲一区二区| 亚洲av电影不卡..在线观看| 我的女老师完整版在线观看| 午夜视频国产福利| 女人十人毛片免费观看3o分钟| www.av在线官网国产| 久久久久九九精品影院| 亚洲精品日韩在线中文字幕 | 深夜精品福利| 国产精品不卡视频一区二区| 毛片一级片免费看久久久久| 99久久久亚洲精品蜜臀av| 亚洲高清免费不卡视频| 欧美日韩综合久久久久久| 99riav亚洲国产免费| 国内少妇人妻偷人精品xxx网站| 欧美3d第一页| 久久久久久久久大av| 99久久精品一区二区三区| 噜噜噜噜噜久久久久久91| 国内精品一区二区在线观看| 国内精品宾馆在线| 蜜臀久久99精品久久宅男| 亚洲人成网站在线播| 国产高清视频在线观看网站| av在线亚洲专区| 美女内射精品一级片tv| 男的添女的下面高潮视频| 成年免费大片在线观看| 日韩欧美一区二区三区在线观看| 亚洲经典国产精华液单| 在线免费十八禁| 深夜a级毛片| 久久精品国产鲁丝片午夜精品| 天堂√8在线中文| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 亚洲欧美精品自产自拍| 国产日本99.免费观看| 国产一区二区在线观看日韩| 国产色婷婷99| 日本av手机在线免费观看| 中文亚洲av片在线观看爽| 热99re8久久精品国产| av国产免费在线观看| 久久久久久久亚洲中文字幕| 听说在线观看完整版免费高清| 一个人看的www免费观看视频| 国产淫片久久久久久久久| 亚洲欧美精品专区久久| 尤物成人国产欧美一区二区三区| 淫秽高清视频在线观看| 欧美最新免费一区二区三区| 国产人妻一区二区三区在| www日本黄色视频网| 亚洲在线观看片| 免费观看人在逋| 亚洲国产精品合色在线| 久久久国产成人免费| 久久国产乱子免费精品| 精品久久久噜噜| 少妇裸体淫交视频免费看高清| 欧美日韩综合久久久久久| 99热6这里只有精品| 国产伦精品一区二区三区视频9| 国产成人a区在线观看| 一级毛片久久久久久久久女| 亚洲电影在线观看av| 欧美丝袜亚洲另类| 91麻豆精品激情在线观看国产| 久久久欧美国产精品| 嫩草影院入口| 亚洲精品日韩在线中文字幕 | 深夜精品福利| 久久精品国产亚洲av涩爱 | 国产亚洲5aaaaa淫片| 免费黄网站久久成人精品| 成年女人永久免费观看视频| h日本视频在线播放| 国产美女午夜福利| 国产伦精品一区二区三区视频9| 精品久久久久久久人妻蜜臀av| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 亚洲第一电影网av| 噜噜噜噜噜久久久久久91| 日本成人三级电影网站| 国产精品一二三区在线看| 久久中文看片网| 中国美白少妇内射xxxbb| 18禁裸乳无遮挡免费网站照片| av天堂在线播放| 久久人妻av系列| 综合色丁香网| 悠悠久久av| 欧美成人免费av一区二区三区| 亚洲国产欧美在线一区| www.av在线官网国产| 熟女人妻精品中文字幕| 九九久久精品国产亚洲av麻豆| 欧美高清成人免费视频www| 一级毛片我不卡| 黄片无遮挡物在线观看| 啦啦啦观看免费观看视频高清| 两个人视频免费观看高清| 午夜福利高清视频| 国产乱人偷精品视频| 欧美性感艳星| 日本黄大片高清| 中文字幕免费在线视频6| 成人特级av手机在线观看| 国产老妇女一区| 亚洲一区高清亚洲精品| 国产国拍精品亚洲av在线观看| 欧美不卡视频在线免费观看| 少妇人妻精品综合一区二区 | 欧美成人免费av一区二区三区| 床上黄色一级片| h日本视频在线播放| 麻豆精品久久久久久蜜桃| 黑人高潮一二区| 黄色配什么色好看| 51国产日韩欧美| www.色视频.com| 久久亚洲国产成人精品v| 久久国内精品自在自线图片| 日韩中字成人| 午夜免费男女啪啪视频观看| 美女cb高潮喷水在线观看| www日本黄色视频网| 午夜视频国产福利| 插阴视频在线观看视频| av免费观看日本| 国产av麻豆久久久久久久| 亚洲精品久久国产高清桃花| 国产真实伦视频高清在线观看| 欧美一级a爱片免费观看看| 免费搜索国产男女视频| 亚洲真实伦在线观看| 精品久久久噜噜| 男女视频在线观看网站免费| 只有这里有精品99| 免费观看精品视频网站| 午夜福利在线观看吧| 中文欧美无线码| 亚洲丝袜综合中文字幕| 最后的刺客免费高清国语| 少妇人妻一区二区三区视频| a级一级毛片免费在线观看| 18禁裸乳无遮挡免费网站照片| 97人妻精品一区二区三区麻豆| 国产伦精品一区二区三区四那| 高清毛片免费观看视频网站| 免费观看a级毛片全部| 看黄色毛片网站| 日本欧美国产在线视频| 少妇丰满av| 色5月婷婷丁香| 久99久视频精品免费| 国产精品国产三级国产av玫瑰| 大香蕉久久网| 色视频www国产|