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

    球面配流副二維穩(wěn)態(tài)壓力場的數(shù)值求解方法研究

    2010-12-03 09:46:46李小金苑士華胡紀濱
    中國機械工程 2010年2期
    關(guān)鍵詞:配流油口球面

    李小金 苑士華 胡紀濱

    北京理工大學車輛傳動國家重點實驗室,北京,100081

    0 引言

    缸體與配流盤之間的平衡分析是軸向柱塞液壓元件設計中不可缺少的環(huán)節(jié),而精確求得配流副間油膜的壓力場則是平衡分析的關(guān)鍵。國內(nèi)外學者針對平面配流副進行了大量的研究工作。Shute等[1]通過電模擬法計算了平面配流副上高低壓油槽間過渡區(qū)域的壓力分布,并給出了等效延伸角。Yamaguchi等[2-3]、Pan等[4]、陳卓如等[5]、王有榮等[6]利用不同形式的數(shù)值計算方法給出了平面配流副壓力場分布的數(shù)值解。而針對目前應用較多的球面配流副的研究則極少,且多數(shù)為以測量油膜形狀為目的的試驗研究[7-8]。在壓力場研究方面,基本上采用忽略過渡區(qū)壓力場的方法,只研究一維徑向壓力場分布[9],這與實際的情況通常有較大的偏差。

    油膜的二維壓力場可由雷諾方程求得。雷諾方程是非線性二階偏微分方程,通常很難求得解析解,需要使用計算流體力學方法進行數(shù)值求解。有限差分法是其中發(fā)展最為成熟、應用最為廣泛的數(shù)值求解方法之一。有限差分法采用結(jié)構(gòu)化網(wǎng)格進行計算,由于配流副油膜形狀復雜,因此,為保證網(wǎng)格質(zhì)量及求解精度,需要在求解的物理區(qū)域建立貼體網(wǎng)格,貼體網(wǎng)格在計算域中將被映射為矩形網(wǎng)格。貼體網(wǎng)格的生成方法主要有三種:代數(shù)方法、微分方程方法和保角變換方法[10],其中微分方程方法可以處理各種不規(guī)則邊界,應用最為廣泛。但微分方程法需要以平面直角坐標系作為物理區(qū)域的坐標系,而對于球面配流副,其物理區(qū)域為球面,因此,無法將物理區(qū)域通過平面直角坐標系表示。本文探討球面坐標系下貼體網(wǎng)格的生成方法,解決球面配流副二維壓力場的求解問題。

    1 貼體坐標系下球面配流副的雷諾方程

    本文在文獻[11]中提到的雷諾假設條件下,利用不可壓縮黏性流體的雷諾方程求解球面配流副壓力場。張量形式的雷諾方程表達式為:

    由雷諾假設,在油膜厚度方向不計壓力的變化,因此球面配流副的壓力場可以簡化為二維穩(wěn)態(tài)壓力場。式(1)是可以直接在球坐標系下進行數(shù)值求解的,但由于配流副的形狀復雜,求解時會給邊界條件的描述帶來很大困難,同時數(shù)值求解的精度也不高。因此通常選擇在與油膜邊界貼合的貼體坐標系下對雷諾方程進行求解[12-13]。設貼體坐標系為(ξ,η),球面油膜所在球面坐標系為(θ,φ)。應用張量分析的相關(guān)知識[14],并將貼體坐標系下的速度分量以球坐標系下的缸體的轉(zhuǎn)速 ω表示。最終,式(1)在貼體坐標系下可展開為

    2 貼體網(wǎng)格的生成方程

    球面配流副屬于空間曲面,因此已有的網(wǎng)格生成方程并不適用。本文將推導適用于球面配流副的貼體網(wǎng)格生成方程。由微分幾何學可知,空間曲面上的坐標系均應滿足方程[15]:

    式中,r為曲面上任意一點的矢徑;rα、rβ為坐標曲線在該點的切矢量;ξα、ξβ為坐標系的坐標;L為微分算子;g為度規(guī)張量分量gαβ的行列式的值;Δ2為Beltrami二階微分算子;n為空間曲面的單位法矢量為空間曲面的平均曲率。

    式中,P、Q為控制函數(shù)。

    式(4)、式(5)即為油膜貼體網(wǎng)格的生成方程。通過生成方程求解貼體坐標系中每個網(wǎng)格點在球坐標系中的坐標值即可自動生成油膜貼體網(wǎng)格。當不對貼體網(wǎng)格進行特定控制時,可令P=Q=0,而使用不同形式的控制函數(shù)可以對貼體網(wǎng)格的局部密度產(chǎn)生不同的控制效果。常用的控制函數(shù)為[13]

    球面配流副上存在多個油口(對于液壓泵/液壓馬達為2個,對于液壓變壓器為3個),而且各個油口的包角、壓力邊界條件也不盡相同。為了簡化雷諾方程的求解,形成一套統(tǒng)一的建模方法,本文采用在各個油口分別建立貼體坐標系生成貼體網(wǎng)格,然后采用將其兩兩耦合的方法建立整個球面配流副上的貼體網(wǎng)格。這種建模方法的優(yōu)點在于貼體網(wǎng)格的生成與油口數(shù)量及包角無關(guān),而且對于各個貼體坐標系,雷諾方程都有相同的形式,僅僅邊界條件的數(shù)值不同而已。如圖1所示,斜軸式液壓變壓器的球面配流副共有A、T、B 3個油口,因此作切口7—8、9—10、11—12,將物理域分為3塊,分別建立3個貼體坐標系。可以看出,在物理域中邊界復雜的空間曲面在計算域中被映射為3個矩形平面。A、T、B口復雜的壓力邊界在計算域中分別被映射為橫坐標ξ的等值線1—1′、3 —3′和5—5′;回油壓力邊界在計算域中分別被映射為橫坐標ξ的等值線2—2′、4—4′和6—6′,其中 ,切口邊界 7—8 、9 —10、11—12為油口間的耦合邊界,并不是壓力邊界,該處的壓力值需在計算過程中確定。綜上所述,物理域中原本復雜的油膜邊界條件在計算域中變得十分簡單,而求解區(qū)域在計算域中也轉(zhuǎn)化為標準的矩形區(qū)域,從而大大降低了雷諾方程的求解難度,提高了求解精度。

    3 網(wǎng)格生成方程與雷諾方程的求解

    網(wǎng)格生成方程(式(4)、式(5))與雷諾方程(式(2))均為二階偏微分方程,可以用常規(guī)的有限差分迭代法在圖1b所示的計算域中求解。圖2所示為完整的求解流程。

    這里需要說明的是,在迭代求解網(wǎng)格生成方程(式4)時,會發(fā)現(xiàn)在該方程的若干系數(shù)中還包含待求變量θ的三角函數(shù)項,這樣迭代方程將是隱式方程,直接求解將極大影響計算速度。為此,本文提出:在迭代過程中使用待求變量θi,j周圍的4個點的均值替換三角函數(shù)項中出現(xiàn)的θi,j,進而將迭代方程轉(zhuǎn)化為顯式方程的形式,從而大大提高計算速度。表1所示為網(wǎng)格密度不同、但精度相同條件下使用顯式方程與隱式方程求解所花費的時間。從表1可以看出,采用近似平均的方法將隱式方程轉(zhuǎn)化為顯式方程求解明顯縮短了求解時間。

    表1 隱式方程與顯式方程求解時間對比

    求解雷諾方程前,還需確定配流副上各點的油膜厚度。球面配流副的間隙可認為是由兩個球面半徑相同的球面沿軸向偏移形成的偏心圓球縫隙。球面配流副上各點油膜厚度可表示[9]為

    4 算例及討論

    本文在MAT LAB環(huán)境下以斜軸式液壓變壓器為例,對球面配流副的二維穩(wěn)態(tài)壓力場進行了編程求解,并針對不同油膜狀態(tài)下的壓力分布特性進行比較分析。計算時液壓變壓器球面配流副的相關(guān)參數(shù)為:球面半徑R=112mm,密封帶半徑 r1=6.75mm、r2 = 10.35mm、r3=20.75mm、r4=24mm,油口包角 120°均布,配流盤控制角度60°,A口壓力20MPa,T口壓力0.5MPa,B口壓力20MPa,回油壓力0.1MPa。

    圖3所示為求解網(wǎng)格生成方程生成的球面配流副貼體網(wǎng)格。為顯示清晰,圖中所示的網(wǎng)格數(shù)量為30×300(3個計算域,每個計算域劃分為10×100個網(wǎng)格)。

    圖4a所示為e=5μ m 、ω=0時的油膜壓力場分布??梢钥闯鰤毫鲫P(guān)于x軸呈對稱分布。當兩相鄰油口壓力水平較為接近時,油口間過渡區(qū)域大部分的壓力水平接近油口壓力。顯然,在這種情況下過渡區(qū)域的壓力場是不能假設為零的。圖4b為e=10μ m、ω=2000 r/min時的油膜壓力場分布??梢钥闯?壓力場分布與圖4a完全相同。

    綜上,缸體與配流盤處于偏心狀態(tài)(對配流副進行理論分析時常假設為這一狀態(tài))時壓力場分布與偏心距、缸體轉(zhuǎn)速無關(guān)。這一點與常規(guī)的理論計算結(jié)果是一致的,這也是應用剩余壓緊力理論對配流副進行設計以及計算時無法直接確定油膜厚度的原因。此外,本文完全按照真實油膜區(qū)域求解,因而可以獲得包括過渡區(qū)域在內(nèi)的二維穩(wěn)態(tài)壓力場分布。值得注意的是,過渡區(qū)域的壓力場在相鄰油口壓力接近時不能假設為零,其產(chǎn)生的分離力是不能被忽略的,否則將給分離力的計算造成較大的誤差。

    李小寧[9]給出了忽略過渡區(qū)壓力場(認為只在徑向上存在壓力梯度)的情況下,同軸偏心時球面配流副密封帶徑向壓力場分布的解析解。這一解析解雖然不能描述過渡區(qū)的壓力分布,但是在描述油口區(qū)域密封帶的壓力分布時是準確的。圖5所示為配流盤A口區(qū)域內(nèi)外密封帶壓力分布的數(shù)值解與解析解的對比圖。可以看出,在貼體坐標系下求出的數(shù)值解與解析解非常接近。

    由于配流副尺寸通常較小,壓力場也很難通過試驗測得。有關(guān)球面配流副過渡區(qū)域的壓力場分布的試驗數(shù)據(jù)幾乎沒有,因此無法直接對計算結(jié)果進行驗證。而平面配流副是進行理論與試驗研究時最常使用的配流副,對其過渡區(qū)域壓力場的研究則相對較多一些。當球面半徑很大時,可以認為球面配流副即近似為平面配流副。因此,作為對比和驗證,本文使用球面配流副的貼體網(wǎng)格生成方程及相應的雷諾方程進行平面配流副的近似計算,計算時取球面半徑R=100m、配流盤半徑r4=30mm。由簡單的幾何計算可知,此時配流盤的球面高度僅為4.5μ m,完全滿足通常平面配流盤平面度不大于5μ m的要求,這樣的近似是可行的。圖6所示為取量綱一寬度=0.4、量綱一槽寬=0.4、高低壓槽端部區(qū)夾角ε=32°時,計算所得量綱一等壓力分布線與文獻[1]中電模擬法結(jié)果的對比。可以看出,兩者非常接近。圖7所示為取=0.4時,計算得出的等效延伸角與文獻[1]中結(jié)果的對比??梢钥闯?兩者誤差很小,絕大部分結(jié)果誤差小于0.5°。在實際應用中等效延伸角變化1°大約相當于分離力變化1%,這樣的誤差完全可以接受。

    5 結(jié)論

    (1)推導出了基于球面坐標系的貼體網(wǎng)格生成方程,在貼體坐標系下求解雷諾方程得到了球面配流副的二維穩(wěn)態(tài)壓力場,為今后進一步開展球面配流副二維壓力場特性的理論研究提供了條件。

    (2)針對迭代求解貼體網(wǎng)格生成方程時迭代方程為隱式方程的問題,提出了通過近似平均將隱式方程轉(zhuǎn)化為顯式方程的方法,可在不降低求解精度的條件下大大縮短求解時間。

    (3)對液壓變壓器球面配流副不同工況下的壓力場分布進行了計算求解。計算結(jié)果表明壓力場分布與油膜厚度、缸體轉(zhuǎn)速無關(guān)。油口區(qū)域密封帶徑向壓力分布的計算結(jié)果與解析解相吻合。此外,當兩相鄰油口壓力接近時,過渡區(qū)域的壓力場不能忽略。

    (4)用球面配流副近似平面配流副進行了壓力場計算。結(jié)果表明本文得出的計算結(jié)果與采用電模擬法得到的壓力場分布、等效延伸角等結(jié)果非常接近,證明本文提出的二維壓力場計算方法正確、可行。

    [1]Shute N A,Turnbull D E.The Thrust Balancing of Axial Piston Machines[R].Bedford:BHRA,1963.

    [2]Yamaguchi A.Pressure Distribution on Valve Plate of Axial Plunger Pumps and Motors[J].Sensankenkyu,1966,18(7):18-20.

    [3]Yamaguchi A,Fujitani Y,Isoda Y.Characteristics of Fluid Film between a Valve Plate and a Cylinder Block of Axial Piston Pumps and Motors[J].Japan Hydraulics and Pneumatics Society,1984,15(4):314-322.

    [4]Pan H C,Sheng J C,Lu Y X.Finite Difference Computation of Valve Plate Fluid Film in Axial Piston Machines[J].Int.J.Mech.Sci.,1989,31(10):779-791.

    [5]陳卓如,范莉,金朝銘.液壓馬達新型端面配流副液壓分離力的數(shù)值求解及分析[J].機床與液壓,2000(2):7-9.

    [6]王有榮,劉恒龍,柯堅,等.水壓軸向柱塞泵配流副流場的CFD解析[J].液壓與氣動,2007(7):21-24.

    [7]Kim J K,Kim H E,Lee Y B.Measurement of Fluid Film Thickness on the Valve Plate in Oil Hydraulic Axial Piston Pumps(PartII:Spherical Design Effects)[J].Journal of M echanical Science and Technology,2005,19(2):655-663.

    [8]王彬,周華,楊華勇,等.基于電液反饋的軸向柱塞機械潤滑試驗系統(tǒng)[J].哈爾濱工程大學學報,2008,29(4):348-353.

    [9]李小寧.斜軸式軸向柱塞泵主體部件的計算機輔助設計系統(tǒng)[D].哈爾濱:哈爾濱工業(yè)大學,1989.

    [10]周正貴.計算流體力學基礎理論與實際應用[M].南京:東南大學出版社,2008.

    [11]溫詩鑄,黃平.摩擦學原理[M].2版.北京:清華大學出版社,2002.

    [12]Thompson J F,Thames F C,Mastin C W.Automatic Numerical Generation of Body—fitted Curvilinear Coordinates System for Field Containing Any Number of Arbitrary Two—Dimensional Bodies[J].Journal of Computation Physics.1974,15:299-319.

    [13]Thompson J F.Numerical Grid Generation:Foundation and Application[M].North Holland:Elsevier Science Publishing Co.,Inc.,1985.

    [14]呂盤明.張量算法簡明教程[M].合肥:中國科學技術(shù)大學出版社,2005.

    [15]Warsi Z U A.A Note on the Mathematical Formulation of the Problem of Numerical Coordinate Generation[J].Quarterly of Applied M athematics.1983,41:221-236.

    猜你喜歡
    配流油口球面
    重載柱塞泵球面配流副承載特性研究*
    機電工程(2022年12期)2022-12-26 08:21:14
    一種全液壓助力制動叉車液壓系統(tǒng)
    安徽科技(2021年3期)2021-04-06 08:50:16
    微觀織構(gòu)配流副熱-流-固耦合潤滑特性
    球面檢測量具的開發(fā)
    Clinical outcomes of endoscopic management of pancreatic fluid collections in cirrhotics vs non-cirrhotics: Α
    KD424:小型挖掘機用泵
    Heisenberg群上移動球面法的應用——一類半線性方程的Liouville型定理
    帶有轉(zhuǎn)閥換向機構(gòu)的手動油泵
    科技資訊(2016年23期)2016-05-30 18:00:34
    鐵路編組站動態(tài)配流分層模型
    球面穩(wěn)定同倫群中的ξn-相關(guān)元素的非平凡性
    99九九线精品视频在线观看视频| 亚洲国产高清在线一区二区三| 看十八女毛片水多多多| 婷婷色麻豆天堂久久| 国产真实伦视频高清在线观看| 久久这里只有精品中国| 亚洲最大成人av| 搡老妇女老女人老熟妇| 秋霞在线观看毛片| 美女被艹到高潮喷水动态| 国产成人一区二区在线| 午夜福利高清视频| 日韩三级伦理在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美xxxx黑人xx丫x性爽| 男女国产视频网站| 国产美女午夜福利| av在线老鸭窝| 国产精品.久久久| 边亲边吃奶的免费视频| 亚洲怡红院男人天堂| 国国产精品蜜臀av免费| 亚洲综合精品二区| 欧美 日韩 精品 国产| 精品久久久噜噜| 天天一区二区日本电影三级| 中文字幕久久专区| 美女脱内裤让男人舔精品视频| 亚洲欧美精品专区久久| 97在线视频观看| av国产久精品久网站免费入址| 国产中年淑女户外野战色| 欧美成人午夜免费资源| 日韩,欧美,国产一区二区三区| 三级国产精品欧美在线观看| 免费电影在线观看免费观看| 91久久精品国产一区二区成人| 简卡轻食公司| 国产探花在线观看一区二区| 一个人看视频在线观看www免费| 亚洲av二区三区四区| 天天躁夜夜躁狠狠久久av| 国产精品国产三级专区第一集| 亚州av有码| 欧美激情在线99| 亚洲成人久久爱视频| 亚洲国产精品国产精品| 亚洲av二区三区四区| 国产极品天堂在线| 国产亚洲最大av| 美女cb高潮喷水在线观看| 极品教师在线视频| 97超碰精品成人国产| av.在线天堂| 婷婷色综合大香蕉| av一本久久久久| 午夜免费激情av| 国产单亲对白刺激| 麻豆成人av视频| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| 亚洲va在线va天堂va国产| 色视频www国产| 真实男女啪啪啪动态图| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 国产亚洲91精品色在线| 毛片女人毛片| 亚洲精品乱码久久久v下载方式| 国产亚洲精品av在线| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 国产黄片视频在线免费观看| 天堂俺去俺来也www色官网 | 欧美+日韩+精品| 两个人的视频大全免费| 免费av毛片视频| 国产精品女同一区二区软件| 在线免费十八禁| 99久久人妻综合| 亚洲在久久综合| 久久久亚洲精品成人影院| 国产精品国产三级国产av玫瑰| 最近中文字幕高清免费大全6| 亚洲国产色片| 亚洲天堂国产精品一区在线| 国产高清不卡午夜福利| 日本wwww免费看| 国产 亚洲一区二区三区 | 亚洲av中文av极速乱| 久久精品夜夜夜夜夜久久蜜豆| 亚洲18禁久久av| 91精品伊人久久大香线蕉| 嘟嘟电影网在线观看| 免费看不卡的av| h日本视频在线播放| 国产黄频视频在线观看| 亚洲最大成人av| 最近2019中文字幕mv第一页| 欧美变态另类bdsm刘玥| 亚洲真实伦在线观看| 男女边吃奶边做爰视频| 精品人妻熟女av久视频| 精品久久久久久久人妻蜜臀av| 亚洲第一区二区三区不卡| 26uuu在线亚洲综合色| 青青草视频在线视频观看| 少妇人妻精品综合一区二区| av专区在线播放| 国产女主播在线喷水免费视频网站 | 久久这里只有精品中国| 天天躁日日操中文字幕| 色尼玛亚洲综合影院| 亚洲欧美日韩东京热| 国产极品天堂在线| 一个人免费在线观看电影| 91久久精品国产一区二区成人| a级毛色黄片| 啦啦啦韩国在线观看视频| 成年女人看的毛片在线观看| 亚洲自偷自拍三级| 久久久色成人| 成年女人在线观看亚洲视频 | 黄色日韩在线| 亚洲一级一片aⅴ在线观看| 六月丁香七月| 亚洲国产精品国产精品| 网址你懂的国产日韩在线| av免费观看日本| 深爱激情五月婷婷| 秋霞在线观看毛片| 特级一级黄色大片| 街头女战士在线观看网站| 久久综合国产亚洲精品| 亚洲成人一二三区av| 欧美日本视频| 中文欧美无线码| 亚洲最大成人av| 日韩一区二区三区影片| 熟妇人妻久久中文字幕3abv| 免费观看a级毛片全部| 超碰av人人做人人爽久久| 日韩一区二区三区影片| 久久久久国产网址| 欧美精品一区二区大全| 色综合站精品国产| 人妻系列 视频| 亚洲国产成人一精品久久久| 又爽又黄无遮挡网站| 美女高潮的动态| 国产精品蜜桃在线观看| 国产在线一区二区三区精| 欧美成人午夜免费资源| .国产精品久久| 免费播放大片免费观看视频在线观看| 六月丁香七月| 蜜桃久久精品国产亚洲av| 成年版毛片免费区| 七月丁香在线播放| 精品欧美国产一区二区三| 国产高清三级在线| 一区二区三区高清视频在线| 中文精品一卡2卡3卡4更新| 99久国产av精品国产电影| 国产一区有黄有色的免费视频 | 国内精品一区二区在线观看| 永久免费av网站大全| 亚洲av成人av| 亚洲欧美中文字幕日韩二区| 国产综合懂色| 纵有疾风起免费观看全集完整版 | 午夜福利视频1000在线观看| 亚洲va在线va天堂va国产| 少妇人妻一区二区三区视频| 亚洲av成人精品一二三区| 777米奇影视久久| 久久6这里有精品| 99久久精品国产国产毛片| 水蜜桃什么品种好| 国产v大片淫在线免费观看| 观看免费一级毛片| 中国美白少妇内射xxxbb| 欧美最新免费一区二区三区| 免费观看av网站的网址| 免费av不卡在线播放| 久久久久国产网址| 国产老妇女一区| 日韩人妻高清精品专区| 日韩中字成人| 草草在线视频免费看| 亚洲精品乱久久久久久| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| 国产成人a∨麻豆精品| 最近手机中文字幕大全| 久久久久久久久久久丰满| 在线a可以看的网站| av免费在线看不卡| 国产美女午夜福利| 亚洲性久久影院| 日韩伦理黄色片| 午夜福利在线观看吧| 在线观看人妻少妇| 日本免费在线观看一区| 91在线精品国自产拍蜜月| 亚洲精品乱码久久久v下载方式| 欧美最新免费一区二区三区| 国产一区亚洲一区在线观看| 水蜜桃什么品种好| 18+在线观看网站| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 免费看不卡的av| 日韩欧美精品v在线| 国产av国产精品国产| 99热这里只有精品一区| 人人妻人人澡人人爽人人夜夜 | 亚洲国产欧美人成| 久热久热在线精品观看| 1000部很黄的大片| 舔av片在线| 99热6这里只有精品| 国产日韩欧美在线精品| 联通29元200g的流量卡| 国产麻豆成人av免费视频| 乱人视频在线观看| 亚洲va在线va天堂va国产| 亚洲精品456在线播放app| 亚洲国产成人一精品久久久| 欧美潮喷喷水| 国产一区有黄有色的免费视频 | 美女内射精品一级片tv| 水蜜桃什么品种好| 日本黄色片子视频| 联通29元200g的流量卡| 欧美激情久久久久久爽电影| 国产一级毛片七仙女欲春2| 26uuu在线亚洲综合色| 亚洲内射少妇av| 91久久精品国产一区二区成人| 久久99热6这里只有精品| 日韩成人伦理影院| av播播在线观看一区| 一个人观看的视频www高清免费观看| 亚洲国产精品国产精品| 午夜福利在线观看免费完整高清在| 夫妻午夜视频| 午夜精品国产一区二区电影 | 亚洲精品456在线播放app| 久久久久久久久久久免费av| 麻豆国产97在线/欧美| 亚洲国产色片| 欧美+日韩+精品| 精品久久久久久电影网| 人妻少妇偷人精品九色| 伦精品一区二区三区| 国产精品av视频在线免费观看| 亚洲va在线va天堂va国产| 国产激情偷乱视频一区二区| 只有这里有精品99| 成人无遮挡网站| 欧美97在线视频| 亚洲av中文av极速乱| 深爱激情五月婷婷| 久久久久久久亚洲中文字幕| 欧美日韩精品成人综合77777| 六月丁香七月| 国产亚洲精品av在线| 菩萨蛮人人尽说江南好唐韦庄| 男人爽女人下面视频在线观看| 日韩成人伦理影院| 爱豆传媒免费全集在线观看| 婷婷色av中文字幕| 熟女电影av网| 听说在线观看完整版免费高清| 欧美高清性xxxxhd video| 亚洲成色77777| 久久久久久国产a免费观看| 国产精品伦人一区二区| 九九在线视频观看精品| 国产高清国产精品国产三级 | 国产三级在线视频| 一本一本综合久久| 一个人看的www免费观看视频| 深爱激情五月婷婷| 99久久精品一区二区三区| 国产精品久久久久久精品电影| 免费看日本二区| 亚洲精品自拍成人| 成人av在线播放网站| 亚洲成人av在线免费| 少妇裸体淫交视频免费看高清| 久久精品熟女亚洲av麻豆精品 | 日韩欧美三级三区| 99久久精品国产国产毛片| 日本一本二区三区精品| 永久免费av网站大全| 国产91av在线免费观看| 亚洲精品一二三| 肉色欧美久久久久久久蜜桃 | 久久久久久久久久久免费av| 91av网一区二区| 久久久久久伊人网av| 18禁在线无遮挡免费观看视频| 中文在线观看免费www的网站| 国产白丝娇喘喷水9色精品| 日韩一本色道免费dvd| 综合色av麻豆| 少妇人妻精品综合一区二区| 亚洲久久久久久中文字幕| 午夜免费观看性视频| 肉色欧美久久久久久久蜜桃 | 亚洲国产高清在线一区二区三| 男女下面进入的视频免费午夜| 床上黄色一级片| 狂野欧美白嫩少妇大欣赏| 久久久久久久久中文| 国产免费一级a男人的天堂| 亚洲欧美成人综合另类久久久| 国产精品精品国产色婷婷| 一级毛片我不卡| 久久精品国产亚洲av天美| 一个人看的www免费观看视频| av又黄又爽大尺度在线免费看| 久久久久久久久久久丰满| 免费高清在线观看视频在线观看| 国产在视频线精品| 黄色欧美视频在线观看| 91午夜精品亚洲一区二区三区| 成人漫画全彩无遮挡| 嘟嘟电影网在线观看| .国产精品久久| av免费观看日本| 国产精品三级大全| 一区二区三区乱码不卡18| 99久国产av精品| 亚洲丝袜综合中文字幕| 99视频精品全部免费 在线| 久久精品国产自在天天线| 国产国拍精品亚洲av在线观看| 又粗又硬又长又爽又黄的视频| 一级a做视频免费观看| 91在线精品国自产拍蜜月| 观看免费一级毛片| 亚洲熟女精品中文字幕| 国产精品久久久久久精品电影小说 | 人体艺术视频欧美日本| 亚洲最大成人av| 免费黄频网站在线观看国产| 日韩伦理黄色片| 久久精品国产亚洲网站| 精品熟女少妇av免费看| 亚洲18禁久久av| 中文字幕制服av| 亚洲精品乱码久久久久久按摩| 一二三四中文在线观看免费高清| 国产激情偷乱视频一区二区| 免费观看a级毛片全部| 国产淫片久久久久久久久| 成人高潮视频无遮挡免费网站| 久久久久久久久久黄片| 日本三级黄在线观看| 精品人妻偷拍中文字幕| 男女视频在线观看网站免费| 最近中文字幕高清免费大全6| 麻豆av噜噜一区二区三区| 国产毛片a区久久久久| 久久久久久国产a免费观看| 久久久久精品性色| 一级av片app| 亚洲精品色激情综合| 天天躁日日操中文字幕| 国产成人aa在线观看| 免费黄频网站在线观看国产| 精品一区二区三卡| 色综合色国产| 国产精品1区2区在线观看.| ponron亚洲| 97超视频在线观看视频| 蜜桃久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说| 午夜福利视频1000在线观看| 内射极品少妇av片p| 九色成人免费人妻av| 欧美成人午夜免费资源| 婷婷色综合大香蕉| 国产麻豆成人av免费视频| 国产av码专区亚洲av| 国产av不卡久久| 久久精品国产亚洲av天美| 人人妻人人澡欧美一区二区| 久久国内精品自在自线图片| 欧美最新免费一区二区三区| 少妇裸体淫交视频免费看高清| 久久这里有精品视频免费| 久久人人爽人人爽人人片va| 美女xxoo啪啪120秒动态图| 十八禁国产超污无遮挡网站| 在线免费观看的www视频| 免费黄频网站在线观看国产| 久久久国产一区二区| 青春草视频在线免费观看| 国产亚洲午夜精品一区二区久久 | 欧美变态另类bdsm刘玥| 嫩草影院精品99| 欧美高清成人免费视频www| 国产视频首页在线观看| 亚洲性久久影院| 国产黄片视频在线免费观看| .国产精品久久| 亚洲精品国产av成人精品| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 麻豆久久精品国产亚洲av| 国产片特级美女逼逼视频| 国产av码专区亚洲av| 一个人观看的视频www高清免费观看| 亚洲经典国产精华液单| 免费观看无遮挡的男女| 九九在线视频观看精品| 简卡轻食公司| 亚洲乱码一区二区免费版| 婷婷色av中文字幕| 97人妻精品一区二区三区麻豆| 简卡轻食公司| 欧美日韩在线观看h| 97精品久久久久久久久久精品| 久久久午夜欧美精品| 亚洲熟女精品中文字幕| 国产久久久一区二区三区| 亚洲伊人久久精品综合| 亚洲人成网站高清观看| 淫秽高清视频在线观看| 十八禁网站网址无遮挡 | 午夜福利网站1000一区二区三区| 最近2019中文字幕mv第一页| 日本与韩国留学比较| 亚洲一区高清亚洲精品| 精品国产三级普通话版| 国产成人精品一,二区| 高清av免费在线| 国产在线男女| 校园人妻丝袜中文字幕| 亚洲av二区三区四区| 久久国内精品自在自线图片| 精品人妻偷拍中文字幕| 嘟嘟电影网在线观看| 国产三级在线视频| 大香蕉97超碰在线| 久久99热6这里只有精品| 狂野欧美激情性xxxx在线观看| 男人狂女人下面高潮的视频| 日韩,欧美,国产一区二区三区| 亚洲怡红院男人天堂| 亚洲精品久久午夜乱码| 亚洲精品国产av蜜桃| 国内少妇人妻偷人精品xxx网站| 日本av手机在线免费观看| 国产成人精品一,二区| 美女大奶头视频| 好男人在线观看高清免费视频| 日日摸夜夜添夜夜添av毛片| 日韩成人av中文字幕在线观看| 日本熟妇午夜| 国产精品久久久久久av不卡| 免费观看av网站的网址| 午夜激情福利司机影院| 最新中文字幕久久久久| 十八禁国产超污无遮挡网站| av免费在线看不卡| 久久这里有精品视频免费| 免费看美女性在线毛片视频| 岛国毛片在线播放| 久久久色成人| av在线观看视频网站免费| 晚上一个人看的免费电影| 久久99热这里只有精品18| 国产综合懂色| 亚洲最大成人av| 九草在线视频观看| 大香蕉97超碰在线| 只有这里有精品99| 别揉我奶头 嗯啊视频| xxx大片免费视频| kizo精华| 国产免费视频播放在线视频 | 日本av手机在线免费观看| 人人妻人人澡人人爽人人夜夜 | 少妇熟女aⅴ在线视频| 精品久久久久久久末码| 永久网站在线| 国产中年淑女户外野战色| 精华霜和精华液先用哪个| 十八禁国产超污无遮挡网站| 91久久精品国产一区二区三区| 九草在线视频观看| 午夜老司机福利剧场| 天堂中文最新版在线下载 | 精品一区二区三区视频在线| 2021天堂中文幕一二区在线观| 免费黄网站久久成人精品| 七月丁香在线播放| 欧美日本视频| 免费观看无遮挡的男女| av在线蜜桃| 男人狂女人下面高潮的视频| 色吧在线观看| 国产在线男女| 91久久精品国产一区二区成人| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 少妇丰满av| 成年女人在线观看亚洲视频 | 建设人人有责人人尽责人人享有的 | 亚洲av电影不卡..在线观看| 97超视频在线观看视频| 中文资源天堂在线| 亚洲国产最新在线播放| 男女啪啪激烈高潮av片| 丝瓜视频免费看黄片| 久久人人爽人人爽人人片va| 国内精品一区二区在线观看| av女优亚洲男人天堂| 黄片无遮挡物在线观看| 日韩欧美国产在线观看| 麻豆国产97在线/欧美| 国产精品一区二区三区四区免费观看| 能在线免费看毛片的网站| 亚洲精品乱码久久久久久按摩| 欧美高清性xxxxhd video| 日韩av免费高清视频| 看免费成人av毛片| 亚洲精品日本国产第一区| 亚洲乱码一区二区免费版| 老师上课跳d突然被开到最大视频| 亚洲精品国产成人久久av| 免费av不卡在线播放| 国产中年淑女户外野战色| 久久这里只有精品中国| 18+在线观看网站| 午夜老司机福利剧场| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲最大av| 看十八女毛片水多多多| 亚洲精品第二区| 少妇丰满av| 80岁老熟妇乱子伦牲交| 国产精品久久视频播放| freevideosex欧美| 国产 一区 欧美 日韩| 欧美性猛交╳xxx乱大交人| 肉色欧美久久久久久久蜜桃 | 亚洲最大成人中文| 精品一区二区三区视频在线| 成人国产麻豆网| 国产成人午夜福利电影在线观看| videossex国产| 国产三级在线视频| 国产免费又黄又爽又色| 国产高清有码在线观看视频| 欧美精品一区二区大全| 成年av动漫网址| 色综合站精品国产| 国产高清不卡午夜福利| 国产69精品久久久久777片| 美女国产视频在线观看| 欧美激情国产日韩精品一区| 国产淫片久久久久久久久| 3wmmmm亚洲av在线观看| 少妇的逼水好多| av国产久精品久网站免费入址| 三级国产精品欧美在线观看| 亚洲美女搞黄在线观看| 日韩强制内射视频| 欧美日韩综合久久久久久| a级毛色黄片| 亚洲精品乱久久久久久| 亚洲最大成人中文| 午夜激情欧美在线| 久久99热这里只有精品18| 亚洲无线观看免费| 2021少妇久久久久久久久久久| 中文欧美无线码| 午夜精品一区二区三区免费看| 嫩草影院新地址| 国产精品综合久久久久久久免费| 国产乱人偷精品视频| 亚洲av电影在线观看一区二区三区 | 亚洲第一区二区三区不卡| 日韩一区二区三区影片| 爱豆传媒免费全集在线观看| 精品酒店卫生间| 久久人人爽人人片av| 97超视频在线观看视频| 精品酒店卫生间| 麻豆精品久久久久久蜜桃| 国产 亚洲一区二区三区 | 国产伦一二天堂av在线观看| 久久久色成人| 久久99热这里只有精品18| 欧美+日韩+精品| av女优亚洲男人天堂| 美女国产视频在线观看| 亚洲最大成人中文| 51国产日韩欧美| 色吧在线观看| 少妇高潮的动态图| 免费在线观看成人毛片| 久久久久久久久久黄片| 成年版毛片免费区| 国产乱来视频区| 亚洲内射少妇av| 国产熟女欧美一区二区| 欧美 日韩 精品 国产|