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

    六圓弧蛋形斷面共軛水深計算方法的研究

    2015-02-21 13:15:34張志昌李若冰
    關(guān)鍵詞:形心水躍蛋形

    張志昌,賈 斌,李若冰

    (西安理工大學(xué) 水利水電學(xué)院,陜西 西安 710048)

    六圓弧蛋形斷面共軛水深計算方法的研究

    張志昌,賈 斌,李若冰

    (西安理工大學(xué) 水利水電學(xué)院,陜西 西安 710048)

    【目的】 建立明渠六圓弧蛋形斷面水躍共軛水深的計算方法,為該斷面的設(shè)計提供支持。【方法】 通過分塊計算六圓弧蛋形斷面不同水深時的面積、分塊形心位置和總形心位置,并以此為依據(jù)分析相對面積、相對形心位置與相對水深的關(guān)系,根據(jù)動量方程研究并建立水躍共軛水深的計算方法?!窘Y(jié)果】 給出了不同水深時六圓弧蛋形斷面形心和面積的計算公式以及水躍共軛水深的試算方法;擬合了相對斷面形心和相對面積與相對水深的關(guān)系,給出了水躍共軛水深的簡化迭代計算公式,并驗證了公式的收斂性。與采用理論公式的試算法相比,簡化計算方法對同一算例的躍前、躍后水深計算誤差分別為0.04%和0.16%?!窘Y(jié)論】 推導(dǎo)的水躍共軛水深的試算法和簡化計算方法,其計算精度完全可以滿足工程實際要求,可為工程設(shè)計提供參考。

    六圓弧蛋形斷面;水躍;共軛水深;簡化計算

    水躍是明渠水流從急流向緩流過渡時發(fā)生的一種水面突然躍起的局部水流現(xiàn)象。對于矩形斷面,水躍共軛水深已有理論計算公式[1],而對于其他斷面,由于形狀復(fù)雜,水躍共軛水深不易計算。近年來,學(xué)者們對復(fù)雜斷面的水躍共軛水深進(jìn)行了很多研究。王學(xué)斌等[2]通過求解一元四次方程,得到了梯形斷面水躍共軛水深的精確解,但其計算過程仍比較繁瑣。為了簡化計算,劉計良等[3]根據(jù)水躍函數(shù)曲線的性質(zhì)和幾何意義,提出了梯形斷面水躍共軛水深的簡化計算公式;孫道宗[4]給出了梯形斷面水躍共軛水深的迭代方法;趙延風(fēng)等[5]通過引入單位水面寬度的概念,提出了梯形斷面共軛水深的直接計算公式;馬吉明等[6]研究了城門洞型及平底馬蹄型隧洞內(nèi)的水躍計算方法,由于計算比較復(fù)雜,僅給出了水躍共軛水深的計算曲線。張志昌等[7]和李若冰等[8]研究了U型渠道、標(biāo)準(zhǔn)Ⅰ型馬蹄型斷面水躍的共軛水深,并給出了詳細(xì)的推導(dǎo)過程和計算步驟。

    常見的蛋形斷面分為四圓弧和六圓弧2種形式,被認(rèn)為是最優(yōu)的受力斷面[9]。四圓弧蛋形斷面常用于城市排水,關(guān)于其過流能力和臨界水深已有一些研究成果[10-12]。六圓弧蛋形斷面常用于大型的水利工程,例如湖南省臨澧縣青山水輪泵站灌區(qū)的冉鋪灣隧洞[13]、湖南衡東縣的白蓮灌區(qū)武家坳隧洞[14]均采用六圓弧蛋形斷面,但目前對于其水力計算的研究甚少,僅有文獻(xiàn)[15]給出了臨界水深的計算方法,對于水躍共軛水深的研究尚未見相關(guān)文獻(xiàn)報道。為此,本研究詳細(xì)推導(dǎo)了六圓弧蛋形斷面不同水深時斷面面積和斷面形心的計算過程,并從水躍的基本方程出發(fā),給出了水躍共軛水深的試算計算方法;為了簡化計算,還擬合了相對斷面形心和相對面積與相對水深的關(guān)系,給出了共軛水深的簡化計算公式,以期為六圓弧蛋形斷面的設(shè)計與計算提供參考。

    1 六圓弧蛋形斷面的形狀及參數(shù)

    六圓弧蛋形斷面如圖1所示,該斷面由6段圓弧組成,其中弧ab、ce、df半徑最大,均為r1,弧ef的半徑r2=0.426 042r1,弧ac、bd的半徑r3=(17/56)r1,弧ab的半角α1=25.375°,弧ef的角度α2=121.183 8°,弧ce、df的角度α3=29.408 1°,弧ac、bd的角度α4=40.128 11°,渠道最寬處寬度cd與最大半徑r1相同。由于六圓弧蛋形斷面的復(fù)雜性,現(xiàn)以圖1中的ab、cd、ef處作為分界線,對應(yīng)的水深分別為hab、hcd、hef,因此可以將斷面水深分為4個部分,即h≤hab、hab≤h≤hcd、hcd≤h≤hef和hef≤h。

    由圖1可以看出,底部弓形的高度hab和兩側(cè)扇形斷面的高度h2分別為:

    hab=r1-r1cosα1=r1(1-cos 25.375)=0.096 478r1,

    h2=r3sinα4=r3sin 40.128 11=0.644 499r3=0.195 651 4r1。

    2 六圓弧蛋形斷面不同水深時面積和形心位置的計算

    2.1 水深處于底部弓形斷面(h≤hab)

    如圖2所示,φ為水深處于底部弓形斷面時對應(yīng)的半圓心角,0<φ≤25.375°,A0為計算區(qū)域斷面面積,yc0為計算區(qū)域形心距水面的距離。當(dāng)h≤hab時,斷面面積A和形心距水面的距離yc為:

    (1)

    yc=yc0=4r1sin3φ/[3(2φ-sin 2φ)]-r1cosφ。

    (2)

    水深計算公式為:h=r1-r1cosφ。

    (3)

    2.2 水深處于最大寬度(含最大寬度)以下(hab≤h≤hcd)斷面

    水深處于最大寬度以下斷面時,斷面面積和形心距水面的距離如圖3所示。將計算區(qū)域分為底部、兩側(cè)弓形和中央梯形4個部分計算,其面積分別為A1、A2、A3、A4,對應(yīng)形心距水面的距離分別為yc1、yc2、yc3、yc4,則由圖3的幾何關(guān)系可得:

    ab=2r1sinα1=2r1sin 25.375=0.857 082r1,

    c1d1=r1-2r3(1-cosβ)=[11/28+(17/28)cosβ]r1,

    A4=0.5(ab+c1d1)(h2-r3sinβ)=

    [0.624 969 5+(17/56)cosβ][0.195 651 4-

    h3=(r1/2)tanα3=(r1/2)tan 29.408 1=0.281 829r1。

    將A1、h3、r2=0.426 042r1、r3=(17/56)r1代入yc1得:

    則斷面面積A和形心距水面的距離yc為:

    (4)

    yc=(A1yc1+A2yc2+A3yc3+A4yc4)/A=(A1yc1+2A2yc2+A4yc4)/A。

    (5)

    水深計算公式為:

    h=hab+h2-r3sinβ=0.292 129 4r1-(17/56)r1sinβ。

    (6)

    式中:0≤β≤40.128 11°。

    2.3 水深處于ef線(含ef線)以下(hcd≤h≤hef)斷面

    將r3=17r1/56,h2=0.195 651 4r1,α1=25.375°,α4=40.128 11°代入以上各式,整理得斷面面積A為:

    (7)

    則斷面形心距水面的距離為:

    (8)

    水深計算公式為:

    h=hab+h2+r1sinβ1=0.292 129 4r1+r1sinβ1。

    (9)

    式中:0≤β1≤α3。

    圖4 水深處于ef線(含ef線)以下斷面
    Fig.4 Water depth at theefline of cross section

    圖5 水深處于ef線以上斷面內(nèi)
    Fig.5 Water depth in the upper bow of cross section

    2.4 水深處于頂部弓形ef線以上(hef≤h)斷面

    則斷面面積A為:

    (10)

    形心距水面的距離分別為:

    斷面形心距水面的距離yc為:

    (11)

    水深的計算公式為

    h=hab+h2+r1sinα3+r2[cos (α2/2-β2)-

    cos (α2/2)]=0.292 129 4r1+r1sinα3+

    0.426 042r1[cos (α2/2-β2)-cos (α2/2)]。

    (12)

    式中:0≤β2<α2/2。

    3 共軛水深的計算

    水躍共軛水深的一般計算公式[1]為:

    (13)

    式(13)可以寫成:J(h′)=J(h″)。

    (14)

    式中:h′、h″分別為躍前、躍后斷面水深。

    水躍方程的計算一般采用試算法,根據(jù)躍前和躍后水深位置的不同,其組合形式見表1。

    4 斷面面積和形心的簡化計算方法

    由以上推導(dǎo)過程可以看出,除了水深處于底部弓形斷面內(nèi)時,斷面面積和形心計算比較簡單外,當(dāng)水深處于底部弓形斷面以上時,水躍共軛水深的計算十分復(fù)雜,計算工作量大。為了簡化計算,本研究分析了水深大于底部弓形斷面時的相對形心、相對面積與相對水深的關(guān)系,結(jié)果見圖6和圖7。由圖6和圖7擬合可得:

    yc/r1=0.126 3(h/r1)3-0.052 7(h/r1)2+
    0.460 7h/r1-0.005 5,

    (15)

    (16)

    公式(15)和(16)的相關(guān)系數(shù)均為0.999 9。在h/r1=0.095~0.960時,公式(15)的平均誤差為0.585%,最大誤差為1.7%;公式(16)的平均誤差為0.337%,最大誤差為0.71%。

    (17)

    式中:y″c和A″用公式(15)和公式(16)計算,計算時公式(15)、(16)中的h用h″代替。

    如果已知躍后斷面水深h″,則可由公式(13)計算出J(h″),躍前斷面相對水深的迭代公式為:

    (18)

    由式(17)、(18)可以看出,h″/r1和h′/r1均小于1,所以迭代的初值在0到1之間選取。

    下面證明公式的收斂性。先證明公式(17),公式(17)可以寫成下面的函數(shù)關(guān)系

    (19)

    對式(19)求導(dǎo)數(shù),得:

    式中:φ′(h″/r1)為函數(shù)φ(h″/r1)對h″/r1的求導(dǎo)。

    將上式寫成:

    由水躍方程可知:

    由上面的水躍方程可得:

    由此得:

    (20)

    對公式(15)和公式(16)求導(dǎo)數(shù),得:

    (y″c/r1)′=0.378 9(h″/r1)2-0.105 4(h″/r1)+0.460 7,

    (21)

    (22)

    式中:Ak為臨界水深對應(yīng)的面積,Bk為臨界水深對應(yīng)的水面寬度。

    將上式變形為:

    (23)

    將公式(23)代入公式(20),得:

    (24)

    現(xiàn)證明公式(18)的收斂性,與公式(17)的類似,設(shè):

    (25)

    對式(25)求導(dǎo)數(shù),得:

    (26)

    由水躍方程可知:

    Q2/(gA′)+A′y′=J(h″)。

    上式可以變形為:

    (27)

    (28)

    將公式(16)、(27)和(28)代入公式(26),可得:

    (29)

    需要強(qiáng)調(diào)的是,公式(17)、(18)適應(yīng)于躍前和躍后斷面水深在底部弓形斷面以上(圖1中ab線以上)的情況。躍后水深大多數(shù)情況下均在底部弓形斷面以上(圖1中ab線以上),在用公式(18)計算時,如果求得的躍前斷面水深h′

    算例 某泄洪隧洞采用六圓弧蛋形斷面,已知六圓弧蛋形斷面的最大半徑r1=10 m,消力池以上總水頭E0=50 m,流速系數(shù)φ=0.86,渠道通過的流量Q=150 m3/s,試判斷是否會發(fā)生水躍,若發(fā)生水躍,試計算躍后水深。

    已知六圓弧蛋形斷面參數(shù)如下:

    r2=0.426 042r1=4.260 42 m,

    r3=17/56r1=3.035 71 m,

    α1=25.375°,α2=121.183 8°,

    α3=29.408 1°,α4=40.128 11°,

    hab=r1(1-cosα1)=0.964 78 m,

    h2=0.195 651 4r1=1.956 51 m。

    解法一:用理論公式計算水躍共軛水深。壩趾收縮斷面以上的總水頭E0為:

    臨界水深用公式(23)計算,假定臨界水深處于最大寬度以上且頂拱以下的斷面內(nèi),臨界水深對應(yīng)的面積Ak和臨界水深hk分別用公式(7)和公式(9)計算,計算時將公式中的β1用β1k表示,其中下標(biāo)k表示臨界關(guān)系,則:

    hk=0.292 129 4r1+r1sinβ1k。

    臨界水深對應(yīng)的水面寬度Bk由圖4的幾何關(guān)系獲得:

    Bk=r1(2cosβ1k-1)。

    由以上公式試算得β1k=2.403 6°,臨界水深為3.341 m,大于躍前斷面水深,因此會發(fā)生水躍。

    即當(dāng)躍后水深為9.384 1 m時,計算的J(h″)與J(h′)近似相等,此時的水深即為躍后水深。

    解法二:用簡化公式計算水躍共軛水深。仍假定躍前水深大于底部弓形斷面,躍前斷面面積用公式(16)計算,則躍前斷面相對水深的迭代公式為:

    將躍前斷面的水躍函數(shù)J(h′)、流量Q和r1代入公式(17),迭代得h″/r1=0.939 916 7,即躍后水深為h″=9.399 167 m,驗算得J(h″)=410.269 m3。求得的躍后水深與理論公式計算的躍后水深相差0.16%。

    5 結(jié) 語

    本研究根據(jù)水躍的一般方程和斷面分割法,研究了六圓弧蛋形斷面水躍共軛水深的計算方法,給出了詳細(xì)的推導(dǎo)過程,通過優(yōu)化擬合,得出了六圓弧蛋形斷面底部弓形以上的相對形心和相對斷面面積的簡化計算公式。與采用理論公式的試算法相比,簡化計算方法對同一算例得到的躍前水深相差0.04%,躍后水深相差0.16%,表明該計算結(jié)果與傳統(tǒng)試算法一樣精確,但該公式形式簡單,計算方便,計算精度完全可以滿足工程實際要求。

    [1] 張志昌.水力學(xué) [M].北京:中國水利水電出版社,2011.

    Zhang Z C.Hydraulics [M].Beijing:China Water Power Pre-ss,2011.(in Chinese)

    [2] 王學(xué)斌,張 毅.梯形明渠水躍共軛水深的精確解 [J].電網(wǎng)與清潔能源,2013,29(11):118-122.

    Wang X B,Zhang Y.The exact solution of conjugate water depth of trapezoidal cross-section [J].Power System and Clean Energy,2013,29(11):118-122.(in Chinese)

    [3] 劉計良,王正中,楊曉松,等.梯形渠道水躍共軛水深理論計算方法初探 [J].水力發(fā)電學(xué)報,2010,29(5):216-219.

    Liu J L,Wang Z Z,Yang X S,et al.Preliminary study on theoretical method for calculating conjugate depth of trapezoidal channel [J].Journal of Hydroelectric Engineering,2010,29(5):216-219.(in Chinese)

    [4] 孫道宗.梯形斷面渠道中水躍共軛水深計算 [J].江西水利科技,2003,29(3):133-137.

    Sun D Z.The calculation of the conjugate depth about the hydraulic jump of the canal of the trapezidal section [J].Jiangxi Hydraulic Science & Technology,2003,29(3):133-137.(in Chinese)

    [5] 趙延風(fēng),王正中,蘆 琴,等.梯形明渠水躍共軛水深的直接計算方法 [J].山東大學(xué)學(xué)報:工學(xué)版,2009,39(2):131-136.

    Zhao Y F,Wang Z Z,Lu Q,et al.Direct calculation method for conjugate water depth of the trapezoidal open channel [J].Journal of Shandong University:Engineering Science,2009,39(2):131-136.(in Chinese)

    [6] 馬吉明,謝省宗,梁元博.城門洞形及馬蹄型輸入隧洞內(nèi)的水躍 [J].水利學(xué)報,2000(7):20-24.

    Ma J M,Xie S Z,Liang Y B.Hydraulic jumps in rectangular conduit with circular upper wall and horseshoe tunnel [J].Journal of Hydraulic Engineering,2000(7):20-24.(in Chinese)

    [7] 張志昌,李若冰.U形渠道臨界水深、弗勞德數(shù)和水躍的研究 [J].西安理工大學(xué)學(xué)報,2012,28(2):198-203.

    Zhang Z C,Li R B.Research on critical water depth,froude number and hydraulic jump of U-shaped channel [J].Journal of Xi’an University of Technology,2012,28(2):198-203.(in Chinese)

    [8] 李若冰,張志昌.標(biāo)準(zhǔn)Ⅰ型馬蹄型斷面水躍共軛水深的計算 [J].西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2012,40(8):230-234.

    Li R B,Zhang Z C.Calculation of conjugate water depth in standard typeⅠof horseshoe cross section [J].Journal of Northwest A&F University:Natural Science Edition,2012,40(8):230-234.(in Chinese)

    [9] Hilbino S,Motojima M.Stability of a large scale egg-shaped cavern [C]//Kobe.Proc 2nd International Symposium on Field Measurements in Geomechanics.Publ Rotterdam:A A Balkema,1988:683-692.

    [10] Bijankhan M,Kouchakzadeh S.Egg-shaped cross section:Uniform flow direct solution and stability identification [J].Flow Measurement and Instrumentation,2011,22(6):511-516.

    [11] Mohammad A G.Hydraulics of partlally filled egg sewers [J].Journal of Environmental Engineering,1987,113(2):407-425.

    [12] Raikar R V,Shiva Reddy M S,Vishwanadh G K.Normal and critical depth computations for egg-shaped conduit sections [J].Flow Measurement and Instrumentation,2010,21(3):367-372.

    [13] 傅功年,唐耿紅.冉鋪灣隧洞襯砌斷面型式選擇 [J].人民長江,2002,33(5):14-16.

    Fu G N,Tang G H.Type selection of lined cross-section for Ranpuwan tunnel [J].Yangtze River,2002,33(5):14-16.(in Chinese)

    [14] 王子宜.蛋形襯砌在白蓮灌區(qū)武家坳隧洞中的應(yīng)用 [J].湖南水利水電,2006(6):70-71.

    Wang Z Y.Application of egg-shaped lining in Bailian irrigation area Wujiaao tunnel [J].Hunan Hydro & Power,2006(6):70-71.(in Chinese)

    [15] 李若冰,張志昌.明渠六圓弧蛋形斷面臨界水深和收縮斷面水深的計算 [J].武漢大學(xué)學(xué)報,2012,45(4):463-467.

    Li R B,Zhang Z C.Calculation of critical water depth and contraction section water depth in egg-shaped cross-section with six arcs of open channel [J].Engineering Journal of Wuhan University,2012,45(4):463-467.(in Chinese)

    [16] 鄧建中,葛仁杰,程正興.計算方法 [M].西安:西安交通大學(xué)出版社,1985.

    Deng J Z,Ge R X,Cheng Z X.Computing method [M].Xi’an:Xi’an Jiaotong University Press,1985.(in Chinese)

    Calculation method of conjugate water depth of egg-shaped cross section with six arcs in open channel

    ZHANG Zhi-chang,JIA Bin,LI Ruo-bing

    (InstituteofWaterResourcesandHydro-electricEngineering,Xi’anUniversityofTechnology,Xi’an,Shaanxi710048,China)

    【Objective】 The calculation method of conjugate water depth of egg-shaped cross section with six arcs in open channel was established to support the design of this kind of section.【Method】 Blocking method was used to calculate areas and centroids under different water depths,and the relationships between relative area,relative centroid position and relative depth were analyzed.Then the calculation method of conjugate water depth of egg-shaped cross section with six arcs was established according to momentum equation.【Result】 The computational formulas of centroids and areas of egg-shaped cross section with six arcs under different water depths and trial method of conjugate water depth were established.The relationships between the relative area,relative centroid position and relative depth were fitted,the simplified calculation method and processes of conjugate water depth were supplied and the convergence was verified.Compared with the theoretical formula of trial method,the errors of depths before and after water jump by simplified calculation method were 0.04% and 0.16%,respectively.【Conclusion】 The accuracy of simplified calculation method for conjugate water depth met the requirement of actual projects.

    egg-shaped cross section with six arcs;hydraulic jump;conjugate water depth;simplified calculation

    2013-09-22 [作者簡介] 張志昌(1954-),男,陜西西安人,教授級高級工程師,主要從事明渠測流、高速水流等研究。 E-mail:zhangzhichang1954@163.com

    時間:2014-12-12 09:30

    10.13207/j.cnki.jnwafu.2015.01.025

    TV133+.1

    A

    1671-9387(2015)01-0220-09

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20141212.0930.025.html

    猜你喜歡
    形心水躍蛋形
    Heisenberg李代數(shù)的形心
    那吉航運(yùn)樞紐大壩消力池?fù)p毀成因和搶修方案分析
    基于MATLAB圖像特征提取的零件位置識別
    突然擴(kuò)散水躍方程的改進(jìn)與比較
    基于MATLAB圖像特征提取的零件位置識別
    漸擴(kuò)式消力池水躍長度的分析與計算
    蛋形輥道窯結(jié)構(gòu)的窯體內(nèi)部溫度場分析
    蛋形斷面無壓隧洞水面線解析算法
    無壓流六圓弧蛋形斷面臨界水深近似算法
    基于空間混淆位置隱私保護(hù)的位置隱私區(qū)域生成算法*
    国产成人啪精品午夜网站| 亚洲国产毛片av蜜桃av| 免费久久久久久久精品成人欧美视频| 国产精华一区二区三区| 久久国产乱子伦精品免费另类| 精品高清国产在线一区| 免费在线观看黄色视频的| 午夜免费观看网址| 亚洲 国产 在线| 看黄色毛片网站| 日本vs欧美在线观看视频| 日本免费a在线| 亚洲第一av免费看| 在线十欧美十亚洲十日本专区| 日韩欧美在线二视频| 日韩视频一区二区在线观看| 波多野结衣高清无吗| 欧美人与性动交α欧美软件| 少妇的丰满在线观看| 日韩大码丰满熟妇| 国产精品国产av在线观看| 国产一区二区三区视频了| 色尼玛亚洲综合影院| 99香蕉大伊视频| a级片在线免费高清观看视频| 久久久久国产精品人妻aⅴ院| 久久久久精品国产欧美久久久| 亚洲国产欧美一区二区综合| 国产精品二区激情视频| 日韩免费高清中文字幕av| 99香蕉大伊视频| 交换朋友夫妻互换小说| 精品国产亚洲在线| 欧美久久黑人一区二区| 黄片大片在线免费观看| 国产91精品成人一区二区三区| 999精品在线视频| 亚洲精品成人av观看孕妇| 国产无遮挡羞羞视频在线观看| 黄色成人免费大全| 色老头精品视频在线观看| 中文字幕精品免费在线观看视频| 男男h啪啪无遮挡| 精品高清国产在线一区| 免费一级毛片在线播放高清视频 | 亚洲av日韩精品久久久久久密| 精品乱码久久久久久99久播| 久热这里只有精品99| 国产精品久久久av美女十八| 人人妻人人爽人人添夜夜欢视频| 1024香蕉在线观看| 欧美+亚洲+日韩+国产| 少妇 在线观看| 亚洲欧美激情在线| 亚洲成人久久性| 亚洲欧美精品综合一区二区三区| 在线十欧美十亚洲十日本专区| 激情在线观看视频在线高清| 99久久99久久久精品蜜桃| 午夜91福利影院| 成人永久免费在线观看视频| 窝窝影院91人妻| 午夜福利一区二区在线看| 亚洲专区国产一区二区| 岛国在线观看网站| 在线观看66精品国产| 一级毛片高清免费大全| 高清在线国产一区| 成年女人毛片免费观看观看9| 精品无人区乱码1区二区| 欧美黄色淫秽网站| 老司机午夜十八禁免费视频| 91麻豆av在线| 成人18禁高潮啪啪吃奶动态图| 校园春色视频在线观看| 多毛熟女@视频| 女人被狂操c到高潮| svipshipincom国产片| 亚洲成a人片在线一区二区| 99香蕉大伊视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲激情在线av| 久久久精品欧美日韩精品| 亚洲人成77777在线视频| 黑丝袜美女国产一区| 91av网站免费观看| 大陆偷拍与自拍| 国产精品久久久人人做人人爽| 免费少妇av软件| 天天躁夜夜躁狠狠躁躁| 免费在线观看影片大全网站| 在线视频色国产色| 亚洲成人精品中文字幕电影 | 巨乳人妻的诱惑在线观看| 日本免费一区二区三区高清不卡 | 欧美日韩av久久| 又黄又粗又硬又大视频| 国产成人av激情在线播放| 久久草成人影院| 97超级碰碰碰精品色视频在线观看| 亚洲人成伊人成综合网2020| 国产97色在线日韩免费| 欧美日韩视频精品一区| 日韩av在线大香蕉| 黄频高清免费视频| 亚洲熟妇中文字幕五十中出 | 操出白浆在线播放| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久久5区| 久久久水蜜桃国产精品网| 天堂动漫精品| 久久天堂一区二区三区四区| 一进一出抽搐gif免费好疼 | 另类亚洲欧美激情| 国产成人欧美| 欧美日韩亚洲高清精品| 新久久久久国产一级毛片| 曰老女人黄片| 9191精品国产免费久久| 精品国内亚洲2022精品成人| 人妻久久中文字幕网| 中文字幕色久视频| 十分钟在线观看高清视频www| 极品人妻少妇av视频| 亚洲人成电影免费在线| 久久久久国内视频| 亚洲av电影在线进入| 在线天堂中文资源库| 国产不卡一卡二| 亚洲成人免费av在线播放| 纯流量卡能插随身wifi吗| 悠悠久久av| 免费日韩欧美在线观看| 国产成人av激情在线播放| svipshipincom国产片| 免费久久久久久久精品成人欧美视频| 在线观看免费午夜福利视频| 国产色视频综合| 后天国语完整版免费观看| e午夜精品久久久久久久| 国产极品粉嫩免费观看在线| 国产精华一区二区三区| 真人一进一出gif抽搐免费| 国产精品美女特级片免费视频播放器 | 国产精品香港三级国产av潘金莲| 一级毛片高清免费大全| 久久人妻av系列| 欧美一级毛片孕妇| 在线观看免费午夜福利视频| 丰满饥渴人妻一区二区三| 一级毛片高清免费大全| 中文字幕高清在线视频| 黄色毛片三级朝国网站| 两个人看的免费小视频| 国产成人精品久久二区二区免费| 成年版毛片免费区| 午夜视频精品福利| 国产男靠女视频免费网站| 在线免费观看的www视频| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久毛片微露脸| 国产精品一区二区三区四区久久 | 成人亚洲精品一区在线观看| 在线天堂中文资源库| 俄罗斯特黄特色一大片| 叶爱在线成人免费视频播放| 在线观看免费日韩欧美大片| 精品一区二区三区四区五区乱码| 成人18禁高潮啪啪吃奶动态图| 不卡av一区二区三区| 欧美成人午夜精品| 99久久久亚洲精品蜜臀av| 女人高潮潮喷娇喘18禁视频| 精品久久久久久,| 精品人妻1区二区| 美女大奶头视频| 婷婷六月久久综合丁香| 亚洲中文av在线| 日本a在线网址| 国产黄色免费在线视频| 国产精品久久电影中文字幕| 亚洲欧美精品综合一区二区三区| 日韩欧美一区二区三区在线观看| 国产精品 国内视频| 国产精品亚洲av一区麻豆| 看片在线看免费视频| 97碰自拍视频| 大陆偷拍与自拍| 成人三级做爰电影| 老司机在亚洲福利影院| 久久亚洲精品不卡| 极品教师在线免费播放| 国产成人啪精品午夜网站| 9191精品国产免费久久| 不卡一级毛片| 日本欧美视频一区| 亚洲熟妇中文字幕五十中出 | svipshipincom国产片| 91在线观看av| 夜夜看夜夜爽夜夜摸 | 51午夜福利影视在线观看| 热99re8久久精品国产| 夜夜躁狠狠躁天天躁| 亚洲专区中文字幕在线| 天天添夜夜摸| 在线av久久热| 高清欧美精品videossex| 精品一品国产午夜福利视频| 久久精品91无色码中文字幕| 久久热在线av| 免费在线观看影片大全网站| 99热只有精品国产| 最近最新中文字幕大全免费视频| 亚洲人成电影免费在线| 91麻豆av在线| 亚洲男人的天堂狠狠| √禁漫天堂资源中文www| 最新美女视频免费是黄的| 成人精品一区二区免费| 又黄又爽又免费观看的视频| 日韩欧美国产一区二区入口| 久久精品国产清高在天天线| 亚洲成av片中文字幕在线观看| 久久中文字幕人妻熟女| 啪啪无遮挡十八禁网站| 亚洲一区二区三区欧美精品| 人人妻人人澡人人看| 成熟少妇高潮喷水视频| 看片在线看免费视频| 成年版毛片免费区| 国产精品一区二区精品视频观看| 老汉色∧v一级毛片| 欧美日本亚洲视频在线播放| 天天影视国产精品| 岛国视频午夜一区免费看| 久久精品国产亚洲av高清一级| 满18在线观看网站| 一区二区三区精品91| 亚洲一区高清亚洲精品| 精品久久久久久成人av| 久久国产精品人妻蜜桃| 高清欧美精品videossex| 97人妻天天添夜夜摸| 999精品在线视频| 韩国精品一区二区三区| 波多野结衣高清无吗| 日韩 欧美 亚洲 中文字幕| 亚洲熟妇中文字幕五十中出 | 午夜精品久久久久久毛片777| 国产精品亚洲av一区麻豆| 欧美另类亚洲清纯唯美| 怎么达到女性高潮| 香蕉久久夜色| 亚洲熟女毛片儿| 天天影视国产精品| 亚洲成av片中文字幕在线观看| 免费看a级黄色片| 999精品在线视频| 在线观看一区二区三区激情| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品久久久久久人妻精品电影| cao死你这个sao货| 国产成人影院久久av| 性欧美人与动物交配| 9色porny在线观看| 水蜜桃什么品种好| 欧美av亚洲av综合av国产av| 男女床上黄色一级片免费看| 欧美日韩中文字幕国产精品一区二区三区 | 老司机午夜福利在线观看视频| 国产成人免费无遮挡视频| 日本黄色视频三级网站网址| 最好的美女福利视频网| 日本 av在线| av在线天堂中文字幕 | 久久久久久亚洲精品国产蜜桃av| 亚洲欧洲精品一区二区精品久久久| 中文字幕av电影在线播放| 在线十欧美十亚洲十日本专区| 一进一出抽搐gif免费好疼 | 精品久久蜜臀av无| 最近最新中文字幕大全电影3 | 九色亚洲精品在线播放| 好男人电影高清在线观看| 日韩 欧美 亚洲 中文字幕| 精品国产美女av久久久久小说| 中文字幕人妻丝袜一区二区| 在线天堂中文资源库| 亚洲狠狠婷婷综合久久图片| 如日韩欧美国产精品一区二区三区| 法律面前人人平等表现在哪些方面| 淫妇啪啪啪对白视频| 在线观看免费高清a一片| 97人妻天天添夜夜摸| 最好的美女福利视频网| 亚洲激情在线av| 久久精品成人免费网站| 日韩欧美在线二视频| 99精品在免费线老司机午夜| 一边摸一边抽搐一进一出视频| 神马国产精品三级电影在线观看 | 日韩中文字幕欧美一区二区| 成人av一区二区三区在线看| 一a级毛片在线观看| 久久人妻福利社区极品人妻图片| 后天国语完整版免费观看| 久久 成人 亚洲| av免费在线观看网站| 亚洲国产精品sss在线观看 | 99国产精品一区二区蜜桃av| 亚洲成a人片在线一区二区| 黄片播放在线免费| 日韩精品中文字幕看吧| 精品久久久久久,| 欧美日韩乱码在线| 亚洲第一欧美日韩一区二区三区| 精品国产乱子伦一区二区三区| 美女国产高潮福利片在线看| 精品电影一区二区在线| 国产一卡二卡三卡精品| 久久久久久久精品吃奶| 久久久久国产一级毛片高清牌| 一区二区三区国产精品乱码| 欧美+亚洲+日韩+国产| 五月开心婷婷网| 51午夜福利影视在线观看| aaaaa片日本免费| 国产精品野战在线观看 | 天天添夜夜摸| 99香蕉大伊视频| 国产高清国产精品国产三级| 一本综合久久免费| 亚洲国产精品999在线| 国产一区二区三区在线臀色熟女 | 国产亚洲精品一区二区www| 天天躁夜夜躁狠狠躁躁| 亚洲熟妇中文字幕五十中出 | 国产精品久久久久成人av| 悠悠久久av| 长腿黑丝高跟| 日韩大码丰满熟妇| 两个人看的免费小视频| 日韩精品中文字幕看吧| 天天影视国产精品| 亚洲视频免费观看视频| 亚洲精华国产精华精| 在线播放国产精品三级| 在线观看66精品国产| 男女午夜视频在线观看| 亚洲av日韩精品久久久久久密| 午夜福利在线观看吧| 人妻久久中文字幕网| 91在线观看av| 精品人妻1区二区| 国产一区二区激情短视频| 制服诱惑二区| 精品欧美一区二区三区在线| 看黄色毛片网站| 999精品在线视频| 欧美黄色片欧美黄色片| av福利片在线| 亚洲色图综合在线观看| 亚洲成人精品中文字幕电影 | 亚洲视频免费观看视频| 满18在线观看网站| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 日韩 欧美 亚洲 中文字幕| 免费搜索国产男女视频| √禁漫天堂资源中文www| 老司机深夜福利视频在线观看| 91字幕亚洲| 大陆偷拍与自拍| 欧美中文日本在线观看视频| 岛国在线观看网站| 亚洲国产精品sss在线观看 | 午夜精品在线福利| 天堂影院成人在线观看| 激情视频va一区二区三区| 精品乱码久久久久久99久播| 国产精品成人在线| 一级毛片女人18水好多| 老汉色∧v一级毛片| 午夜免费观看网址| 制服诱惑二区| 国产精品一区二区免费欧美| 纯流量卡能插随身wifi吗| 十八禁网站免费在线| 免费在线观看黄色视频的| 国产激情欧美一区二区| 自线自在国产av| 国产有黄有色有爽视频| 狂野欧美激情性xxxx| 老司机深夜福利视频在线观看| 亚洲精品一区av在线观看| 国产精品自产拍在线观看55亚洲| x7x7x7水蜜桃| 国产精品99久久99久久久不卡| 国产不卡一卡二| 一区在线观看完整版| 99国产精品99久久久久| 丁香欧美五月| 十八禁人妻一区二区| 丝袜美足系列| 日韩 欧美 亚洲 中文字幕| 极品教师在线免费播放| 成人亚洲精品av一区二区 | 久久中文看片网| 一级作爱视频免费观看| 在线观看舔阴道视频| 男人舔女人下体高潮全视频| 久久九九热精品免费| 久久久久国产精品人妻aⅴ院| 国内毛片毛片毛片毛片毛片| 人人妻人人添人人爽欧美一区卜| 最好的美女福利视频网| cao死你这个sao货| 欧美一级毛片孕妇| 亚洲人成网站在线播放欧美日韩| www.999成人在线观看| 这个男人来自地球电影免费观看| 法律面前人人平等表现在哪些方面| 动漫黄色视频在线观看| videosex国产| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 久久中文字幕人妻熟女| 日韩欧美在线二视频| 亚洲成人久久性| 日韩av在线大香蕉| 黄色怎么调成土黄色| 午夜免费成人在线视频| 国产精品 国内视频| xxxhd国产人妻xxx| 亚洲伊人色综图| 超色免费av| 亚洲av成人一区二区三| 午夜福利一区二区在线看| av片东京热男人的天堂| 在线观看www视频免费| 极品教师在线免费播放| 黄色丝袜av网址大全| 久久久久久人人人人人| 男女之事视频高清在线观看| 高清在线国产一区| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 麻豆av在线久日| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 日本黄色日本黄色录像| 美女扒开内裤让男人捅视频| 国产精品综合久久久久久久免费 | 老熟妇乱子伦视频在线观看| 国产91精品成人一区二区三区| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 亚洲av美国av| 国产成人欧美在线观看| 精品国产超薄肉色丝袜足j| 国内毛片毛片毛片毛片毛片| 久久人人爽av亚洲精品天堂| 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 国产乱人伦免费视频| 99在线人妻在线中文字幕| 法律面前人人平等表现在哪些方面| 国产欧美日韩一区二区三区在线| 麻豆av在线久日| 香蕉久久夜色| 亚洲熟妇中文字幕五十中出 | 国产一区二区三区视频了| 99精国产麻豆久久婷婷| www.熟女人妻精品国产| 国产精品综合久久久久久久免费 | 日韩精品青青久久久久久| 母亲3免费完整高清在线观看| 99精品在免费线老司机午夜| 午夜福利一区二区在线看| 久久这里只有精品19| 少妇 在线观看| 久久精品人人爽人人爽视色| 美女扒开内裤让男人捅视频| 亚洲 国产 在线| 激情视频va一区二区三区| 香蕉丝袜av| 丁香六月欧美| 日韩欧美三级三区| 成人亚洲精品av一区二区 | 欧美激情高清一区二区三区| 级片在线观看| 美女大奶头视频| x7x7x7水蜜桃| 丝袜在线中文字幕| 午夜免费激情av| 97人妻天天添夜夜摸| 日本a在线网址| 久久精品国产综合久久久| 99国产精品免费福利视频| 久久九九热精品免费| 国产成人av激情在线播放| 欧美在线黄色| 三级毛片av免费| 亚洲国产欧美日韩在线播放| 一级,二级,三级黄色视频| 中文字幕精品免费在线观看视频| 19禁男女啪啪无遮挡网站| 高清在线国产一区| 久久国产精品男人的天堂亚洲| 操出白浆在线播放| 不卡一级毛片| 日韩大码丰满熟妇| 日韩精品中文字幕看吧| 极品教师在线免费播放| 美女福利国产在线| 深夜精品福利| 中文字幕最新亚洲高清| 国产xxxxx性猛交| 两人在一起打扑克的视频| 嫩草影视91久久| 国产国语露脸激情在线看| ponron亚洲| 欧美乱妇无乱码| 中文亚洲av片在线观看爽| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 中文字幕精品免费在线观看视频| 国产免费现黄频在线看| 丰满迷人的少妇在线观看| 国产高清激情床上av| 亚洲国产中文字幕在线视频| 男人舔女人下体高潮全视频| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影 | www国产在线视频色| 日本五十路高清| 久久久精品国产亚洲av高清涩受| 大型黄色视频在线免费观看| 制服诱惑二区| ponron亚洲| 一级作爱视频免费观看| 激情视频va一区二区三区| 18美女黄网站色大片免费观看| 亚洲专区国产一区二区| 电影成人av| 757午夜福利合集在线观看| 亚洲国产精品一区二区三区在线| 99久久99久久久精品蜜桃| 国产一区二区三区在线臀色熟女 | 亚洲成av片中文字幕在线观看| 欧美久久黑人一区二区| 宅男免费午夜| 欧美大码av| 国产成年人精品一区二区 | 欧美黑人精品巨大| 国产亚洲欧美精品永久| 搡老岳熟女国产| 亚洲欧美精品综合久久99| 天天躁狠狠躁夜夜躁狠狠躁| 国产单亲对白刺激| 在线观看一区二区三区| 国产高清国产精品国产三级| 91麻豆精品激情在线观看国产 | av欧美777| 国产片内射在线| 日韩欧美三级三区| 亚洲欧美日韩高清在线视频| 亚洲精品粉嫩美女一区| 精品高清国产在线一区| 精品国产乱码久久久久久男人| 免费日韩欧美在线观看| 男女下面插进去视频免费观看| 亚洲色图av天堂| 亚洲av熟女| 亚洲欧美精品综合一区二区三区| 日韩人妻精品一区2区三区| 国产亚洲精品久久久久久毛片| 亚洲五月婷婷丁香| 身体一侧抽搐| 69精品国产乱码久久久| 精品电影一区二区在线| 女警被强在线播放| 大型黄色视频在线免费观看| 精品国产一区二区久久| 国产男靠女视频免费网站| 国产区一区二久久| 久久国产乱子伦精品免费另类| 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产一区二区三区四区第35| www.精华液| 美女高潮喷水抽搐中文字幕| 黄色毛片三级朝国网站| 日韩欧美免费精品| 美女国产高潮福利片在线看| 中亚洲国语对白在线视频| 国产99白浆流出| 国产国语露脸激情在线看| 久久人人精品亚洲av| 女人精品久久久久毛片| 美女午夜性视频免费| 日本撒尿小便嘘嘘汇集6| 丁香欧美五月| 精品久久久久久成人av| av视频免费观看在线观看| 伦理电影免费视频| 99国产精品免费福利视频| 看片在线看免费视频| 久久久精品欧美日韩精品| 中文字幕精品免费在线观看视频| 国产成人精品久久二区二区91| av免费在线观看网站| 久久九九热精品免费|