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

    Boussinesq類方程描述有限沙壩地形Bragg反射的性能比較

    2010-09-05 03:59:48張俊張日向張永剛
    海洋通報(bào) 2010年3期
    關(guān)鍵詞:沙壩入射波反射系數(shù)

    張俊,張日向,張永剛

    (1. 大連理工大學(xué)海岸及近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116023;2. 海軍大連艦艇學(xué)院軍事海洋系,遼寧 大連 116018)

    Boussinesq類方程描述有限沙壩地形Bragg反射的性能比較

    張俊1,張日向1,張永剛2

    (1. 大連理工大學(xué)海岸及近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116023;2. 海軍大連艦艇學(xué)院軍事海洋系,遼寧 大連 116018)

    利用Boussinesq方程,采用線性攝動(dòng)展開(kāi)法,求解波浪正向通過(guò)有限沙壩地形的一階反射波解,研究比較Boussinesq類方程描述沙壩地形對(duì)波浪的反射作用的性能。通過(guò)研究反射系數(shù)并與勢(shì)流理論結(jié)果及實(shí)驗(yàn)結(jié)果對(duì)比,發(fā)現(xiàn):鄒志利的高階方程、張永剛及Madsen的方程適用的水深范圍較廣,而Nwogu和Peregrine的方程僅在共振點(diǎn)附近有效;當(dāng)入射波的波長(zhǎng)為沙壩波長(zhǎng)兩倍時(shí),反射波產(chǎn)生共振效應(yīng)(即Bragg 反射),反射系數(shù)與沙壩振幅和水深的比值以及地形中沙壩的條數(shù)成正比;相對(duì)勢(shì)流理論的共振時(shí)的反射系數(shù),以張永剛為代表的一系列Boussinesq方程色散精度越高,適用水深范圍越廣,而高階方程適用的水深很淺。

    Boussinesq類方程;有限沙壩;線性攝動(dòng)展開(kāi);一階反射波解; Bragg 反射

    Boussinesq方程能夠用于模擬表面重力波傳播過(guò)程中的折射、繞射、反射及淺化,非線性作用等現(xiàn)象。為改善經(jīng)典Boussinesq方程的色散性和非線性,許多科學(xué)工作者提出了一系列改善色散性的技術(shù)和高階方程,形成Boussinesq類方程[1]。致力于改善線性色散性具有代表性的是張永剛等[2]采用不同的兩個(gè)水深層水平速度變量組合導(dǎo)出的新型Boussinesq方程,選擇方程中的速度變量,該方程代表一系列不同型式的Boussinesq方程,例如Peregrine[3]、Madsen[4]及Nwogu[5]的方程。對(duì)于高階方程,有代表性的有Madsen等[6]精確到的方程及鄒志利[7]適應(yīng)任意海底地形的非線性和色散性精確到的方程。

    近海岸附近通常存在數(shù)條周期變化的呈一定間距的規(guī)則沙壩。關(guān)于沙壩的形成機(jī)理及其與波浪的相互作用,最先是二戰(zhàn)盟軍諾曼底登陸時(shí)提出該課題,Keulegan通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)平底沙灘上,在長(zhǎng)期的正弦波浪入射下,會(huì)沿著波浪破碎線形成一些周期的沙壩[8]。Davies[9]基于勢(shì)流理論,利用線性攝動(dòng)展開(kāi)法求解得波浪通過(guò)有限正弦沙壩地形的遠(yuǎn)場(chǎng)理論波解,發(fā)現(xiàn)當(dāng)入射波的波長(zhǎng)為沙壩波長(zhǎng)的兩倍時(shí),反射波將會(huì)發(fā)生共振,即Bragg 反射,并給出反射波共振時(shí)的反射系數(shù)公式。該結(jié)論隨后由Heathershaw[10]及Davies[11]的實(shí)驗(yàn)所證實(shí),并引發(fā)了人們利用人造沙壩來(lái)保護(hù)岸線的思路。

    由于近岸通常是淺水條件,Boussinesq方程具有較好的適應(yīng)性。本文利用Bossinesq類方程中的幾個(gè)代表方程,忽略底部摩擦及波浪傳播過(guò)程中的能量損失,采用線性攝動(dòng)展開(kāi)法求解出沙壩地形引起的反射波解,給出Bragg反射發(fā)生時(shí)的反射系數(shù)公式。將上述Boussinesq類方程描述Bragg反射的結(jié)果與Davies 勢(shì)流理論結(jié)果及Heathershaw的實(shí)驗(yàn)結(jié)果作對(duì)比,發(fā)現(xiàn)不同Boussinesq方程描述沙壩對(duì)波浪的反射效應(yīng)的性能有差異,為選用Bossinesq方程研究Bragg反射提供理論依據(jù)。

    1 有限沙壩地形

    圖1 有限沙壩的輪廓圖Fig. 1 Profile of finite sandbar topography

    2 Boussinesq類方程中的代表方程

    張永剛等[2]采用兩個(gè)水深層水平速度變量組合導(dǎo)出的新型Boussinesq方程的一階線性形式如下:

    鄒志利[7]適用于任意海底地形的新型高階方程,方程的色散性近似到O(μ4)階,其一維線性形式如下:

    3 Bragg反射問(wèn)題求解

    以張永剛等的Boussinesq方程為例,采用攝動(dòng)法求解波浪通過(guò)有限沙壩地形的反射波解。

    零階方程組:

    一階方程組:

    零階方程組(3.1)描述平底時(shí)波浪的傳播;一階方程組(3.2)反映緩慢變化地形對(duì)波浪傳播的作用。

    3.1 零階方程的波解

    同時(shí)有:

    3.2 一階方程的波解

    3.2.1 一階方程傅立葉變換(變量x) 由傅立葉變換定義,方程組(3.2)變換為:

    聯(lián)立方程組,消去uα1、uβ1,有:

    式(3.9?3.10)中:

    對(duì)方程(3.10)進(jìn)行傅立葉逆變換,得到一階波面升高:

    3.2.2 無(wú)窮積分 無(wú)窮積分式(3.13)可記為:

    式(3.14)中各項(xiàng)關(guān)于變量λ的無(wú)窮積分被積表達(dá)式可表示為(C=x或者和Qiλ分別5次和7次多項(xiàng)式。該類積分可以通過(guò)圍道積分來(lái)實(shí)現(xiàn)。的三個(gè)實(shí)根將對(duì)應(yīng)沙壩地形上的周期波解,而虛根對(duì)非傳播態(tài)的波解有貢獻(xiàn)。由于非傳播態(tài)波解隨傳播距離快速衰減,這里不予考慮。利用留數(shù)定理有[13],C>0時(shí),有:

    可以證明,C<0時(shí),有:

    式(3.16)證明從略。

    根據(jù)留數(shù)定義,有:

    給出x<0側(cè)的反射波波面升高如下:

    定義反射系數(shù)為反射波波幅與入射波波幅比值,得到反射系數(shù)為:

    同理,利用鄒志利的高階方程求解出描述有限沙壩Bragg反射的反射系數(shù)公式,形式同上,結(jié)果略。

    4 結(jié)果分析與驗(yàn)證

    由式(3.20)及(3.21)可見(jiàn):反射系數(shù)和沙壩振幅同水深的比值成正比;入射波與沙壩波長(zhǎng)滿足時(shí),反射波發(fā)生共振,此時(shí),反射系數(shù)與沙壩地形包含單位沙壩條數(shù)成正比。

    Davies[9]給出水平海床上含有數(shù)條小波幅的沙壩時(shí)的反射系數(shù)公式:

    圖2給出本文由不同形式的Boussinesq方程求解的有限沙壩共振系數(shù)(式(3.21))與相對(duì)水深的變化關(guān)系,并與式(4.2)的Davies結(jié)果對(duì)比。圖中假設(shè),地形中的沙壩條數(shù)m=10,入射波長(zhǎng)等于兩倍沙壩波長(zhǎng)。

    由圖2可見(jiàn),Davies勢(shì)流理論求解的反射波發(fā)生共振時(shí)的反射系數(shù)隨相對(duì)水深增大而減小,相對(duì)水深達(dá)到0.6時(shí),反射系數(shù)接近零。相對(duì)勢(shì)流理論的結(jié)果,本文選用的以張永剛等人的方程為代表的Boussinesq方程,線性頻散精度越高,在描述Bragg反射時(shí)適用相對(duì)水深范圍越大:Peregine方程在0.15倍波長(zhǎng)水深范圍內(nèi)與勢(shì)流理論結(jié)果誤差控制在5%,隨相對(duì)水深增大,反射系數(shù)相對(duì)勢(shì)流理論結(jié)果偏?。徽`差控制在5%時(shí),Madsen方程適用相對(duì)水深范圍為0.25,Nwogu方程達(dá)到0.33,張永剛等的方程能達(dá)到0.47,且三者的反射系數(shù)均隨相對(duì)水深增大而偏大;而本文選用的高階方程雖然色散精度達(dá)到0(μ)4,但描述Bragg 反射的反射系數(shù)時(shí),適用水深很淺,相對(duì)水深在0.1以內(nèi)時(shí)誤差控制在3%,隨相對(duì)水深增大,相對(duì)勢(shì)流理論結(jié)果偏大。

    圖2 發(fā)生共振時(shí),本文選用的五個(gè)Boussinseq方程計(jì)算的反射系數(shù)與Davies勢(shì)流理論的結(jié)果的比較Fig. 2 Comparison of reflection coefficient when oscillation

    為驗(yàn)證Davies[9]勢(shì)流理論結(jié)果,Heathershaw[10]在Coastal Engineering Research Center, Virginia進(jìn)行了物理模型試驗(yàn)。在水槽中設(shè)置振幅b=0.05m,條數(shù)m=10,波長(zhǎng)的剛性正弦沙壩,沙壩與水深比值。通過(guò)調(diào)整入射波波長(zhǎng)L,在范圍內(nèi)給出實(shí)測(cè)反射系數(shù)。

    圖3給出在該實(shí)驗(yàn)參數(shù)下,本文選用的Boussinesq方程計(jì)算結(jié)果(式(3.20))與由式(4.1)計(jì)算的Davies[9]理論結(jié)果以及Heathershaw[10]實(shí)驗(yàn)結(jié)果的關(guān)于反射系數(shù)隨的變化關(guān)系的對(duì)比。

    由圖3可見(jiàn),在該實(shí)驗(yàn)參數(shù)下,2k/ kd在(0,1.2)范圍內(nèi)(即入射波長(zhǎng)與沙壩波長(zhǎng)滿足L>1.67Ld)時(shí),由鄒志利、張永剛等及Madsen方程計(jì)算的結(jié)果與Davies勢(shì)流理論得到的結(jié)果一致吻合,但隨著入射波長(zhǎng)減小,即相對(duì)水深增大,Boussinesq方程的結(jié)果呈周期振蕩并逐漸增大,與實(shí)際經(jīng)驗(yàn)背離;Peregrine和Nwogu的方程僅在共振點(diǎn)附近與勢(shì)流理論比較吻合,淺水時(shí)結(jié)果偏小,而在水深較大時(shí)結(jié)果振蕩并偏大,與實(shí)際不符。需要指出的是,該實(shí)驗(yàn)中,共振點(diǎn)處的相對(duì)水深為,高階方程的計(jì)算結(jié)果偏大,而Peregrine的結(jié)果偏小,這與圖2的結(jié)論是吻合的。另外,入射波長(zhǎng)與沙壩波長(zhǎng)滿足時(shí),鄒志利、張永剛等及Madsen的方程的結(jié)果與試驗(yàn)結(jié)果趨勢(shì)一致,基本吻合,但是在共振點(diǎn)處Boussinesq方程的結(jié)果偏大。

    圖3 本文選用的五個(gè)Boussinesq方程的計(jì)算結(jié)果與Davies理論結(jié)果及Heathershaw實(shí)驗(yàn)結(jié)果的對(duì)比Fig. 3 Comparison of Davies’ Potential Theory result and Heathershaw’s experiment result about Reflection Coefficient with several Boussinesq equations.

    5 結(jié) 論

    本文采用線性攝動(dòng)展開(kāi)法,導(dǎo)出基于一組Boussinesq類方程的波浪通過(guò)含有數(shù)條沙壩地形時(shí)的反射波面升高,通過(guò)研究其反射系數(shù)并與勢(shì)流理論結(jié)果及實(shí)驗(yàn)結(jié)果對(duì)比發(fā)現(xiàn):

    1、當(dāng)入射波的波長(zhǎng)等于沙壩波長(zhǎng)的兩倍時(shí),反射波發(fā)生Bragg反射;此時(shí),反射系數(shù)與沙壩振幅和水深的比值、沙壩地形中沙壩的條數(shù)成正比。

    2、Boussinesq方程描述沙壩地形對(duì)波浪的反射作用時(shí),本文所選取的五個(gè)方程中鄒志利的高階方程、張永剛及Madsen的方程適用的水深范圍較廣,而Nwogu和Peregrine的方程僅在共振點(diǎn)附近有效;當(dāng)時(shí),高階方程適用水深范圍最淺(0.1倍波長(zhǎng)水深范圍內(nèi)相對(duì)勢(shì)流理論的反射系數(shù)誤差控制在2%),其他方程色散精度越高適用水深范圍越廣(反射系數(shù)誤差控制在5%時(shí),Peregrine、Madsen、Nwogu及張永剛等的方程適用相對(duì)水深范圍分別為0.15,0.25,0.33,0.47)。

    附錄A

    致謝:感謝威爾士大學(xué)的A G Davies教授提供相關(guān)文獻(xiàn)資料及關(guān)于圍道積分的建議及大連理工大學(xué)倪漢根教授的關(guān)于圍道積分的幫助。

    [1] 陶建華. 水波的數(shù)值模擬 [M]. 天津:天津大學(xué)出版社. 2005:104-108.

    [2] 張永剛,李玉成. 一種新型式的Bossineseq方程 [J]. 科學(xué)通報(bào),1997, 42(21).

    [3] Peregrine D H. Long wave on a beach [J]. Fluid Mech. 1967,27:815-827.

    [4] Madsen P A, Murry R, Sorensen O R. A new form of the Boussinesq equations with improved linear dispersion characteristics. Coast. Engrg. 1991, 15:371-388.

    [5] Nwogu O. An alternative form of the Boussinesq equations for nearshore wave propagation [J]. waterway, port, coastal and ocean Engineering. 1993, 119: 618-638.

    [6] Madsen P A, Schaffer H A. Higher-order Boussinesq–type equations for surface gravity waves: derivation and analysis [J].Phil. Trans. R Soc. Lond. A printed in Great Britain,1998, 356:3123-3184.

    [7] 鄒志利. 水波理論及其應(yīng)用 [M]. 北京: 科學(xué)出版社. 2005:416-418.

    [8] C C Mei, T Hara, J Yu. Longshore Bars and Bragg Resonance [J].

    [9] Some interaction between surface waves and ripples and dunes on the seabed. Inst. Oeanogr. Sci. Rep 108(1980).

    [10] Heathershaw A D. Seadbed-wave resonance and sand bar growth[J]. Nature, 1982, 296: 343-345.

    [11] DAVIES A G, HEATHERSHAW A D. Surface-wave propagation over sinusoidally varying topography [J]. Fluid Mech, 1984, 144:419-443.

    [12] 張永剛, 李玉成. 促進(jìn)其線性頻散特征另一種形式的Boussinesq方程 [J]. 力學(xué)學(xué)報(bào), 1997, 29(2): 142-150.

    [13] Dragoslav S Mitrinovic, Jowan D Keckic. The Cauchy Method of Residues Theory and Applications [M]. Holland: Dordrecht:D.Reidel Publishing Company. 1984: 108-114.

    Comparison among several Boussinesq-type equations about the character of describing Bragg resonance

    ZHANG Jun1, ZHANG Ri-xiang1, ZHANG Yong-gang2
    (1. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116023, China;2. Department of Military Oceanography, Dalian Naval Academy, Dalian 116018, China)

    In order to make a comparison of the character of the Boussinesq-type equations when these equations are used to describe the Bragg resonance, this paper utilized Linear Perturbation Method on the basis of several types of Boussinesq equations to derive the first order reflection wave solution. Through focusing on the reflection coefficient and making comparisons with both Potential flow theory and experiment data, it demonstrated that Zou Zhili’s higher order equation, the equations of Zhang Yonggang and Madsen are appropriate for a wide range of relative water depth,by contrast, Nwogu’s and Peregrine’s are only available in the vicinity of L=2Ld. For the special cases, where the Bragg Resonance occurs, the reflection coefficient is proportional to both the ratio of amplitude of sandbar to the water depth and to the number of bars; in addition, for these Boussinesq-type equations interpreted by Zhang Yonggang, the one with higher precision of the linear dispersion is applicable to wider range of relative water depth, while the higher order equations of Zou Zhili are only appropriate for the relatively shallow water.

    Boussinesq-type equations; finite sandbar; linear perturbation method; first order reflection wave;Bragg resonance

    TV139.2

    A

    1001-6932(2010)02-0295-07

    2009-09-14;

    2009-10-19

    張俊(1984-),男,碩士,主要從事近岸水動(dòng)力研究。電子郵箱:zhangjun.dlut@gmail.com

    猜你喜歡
    沙壩入射波反射系數(shù)
    1977—2017年芝罘連島沙壩地貌演變
    SHPB入射波相似律與整形技術(shù)的試驗(yàn)與數(shù)值研究
    Contribution to the taxonomy of the genus Lycodon H.Boie in Fitzinger,1827 (Reptilia:Squamata:Colubridae) in China,with description of two new species and resurrection and elevation of Dinodon septentrionale chapaense Angel,Bourret,1933
    多道隨機(jī)稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    瞬態(tài)激勵(lì)狀態(tài)下樁身速度以及樁身內(nèi)力計(jì)算
    球面波PP反射系數(shù)的頻變特征研究
    海岸沙壩剖面和灘肩剖面特征研究
    對(duì)機(jī)械波半波損失現(xiàn)象的物理解釋
    電子科技(2015年11期)2015-03-06 01:32:24
    沙質(zhì)沉積物反射系數(shù)的寬帶測(cè)量方法
    后彎管式波力發(fā)電裝置氣室結(jié)構(gòu)的試驗(yàn)研究*
    国产又色又爽无遮挡免| 国产成人午夜福利电影在线观看| 亚洲电影在线观看av| 久久久久精品久久久久真实原创| 日本爱情动作片www.在线观看| 简卡轻食公司| 亚洲av中文字字幕乱码综合| 亚洲电影在线观看av| 久久久成人免费电影| 国产男人的电影天堂91| 国产精品一区www在线观看| 亚洲精品456在线播放app| 日韩欧美在线乱码| 国产三级中文精品| 成人特级av手机在线观看| 日韩亚洲欧美综合| 免费观看精品视频网站| 日韩一本色道免费dvd| 国产男人的电影天堂91| 国产黄色视频一区二区在线观看 | 亚洲第一区二区三区不卡| 久久精品国产自在天天线| 一夜夜www| 建设人人有责人人尽责人人享有的 | 亚洲国产成人一精品久久久| 天堂影院成人在线观看| 亚洲国产精品久久男人天堂| 国产91av在线免费观看| 亚洲美女搞黄在线观看| 99热这里只有是精品50| 天堂网av新在线| 国产色婷婷99| 亚洲人成网站在线观看播放| 欧美变态另类bdsm刘玥| 一个人看视频在线观看www免费| 久久久久久久久久久免费av| 亚洲经典国产精华液单| 精品一区二区三区人妻视频| 国产av在哪里看| 免费av毛片视频| 嫩草影院精品99| 搡老妇女老女人老熟妇| 欧美性感艳星| 成人亚洲精品av一区二区| 18禁动态无遮挡网站| 国产亚洲av片在线观看秒播厂 | 久久人人爽人人片av| 又爽又黄无遮挡网站| 观看免费一级毛片| 内射极品少妇av片p| 国产乱人视频| 少妇的逼好多水| 国产在视频线精品| 成人亚洲精品av一区二区| 久久久久免费精品人妻一区二区| www.av在线官网国产| 日韩av在线免费看完整版不卡| 国产欧美另类精品又又久久亚洲欧美| 91精品伊人久久大香线蕉| 水蜜桃什么品种好| 别揉我奶头 嗯啊视频| 久久亚洲国产成人精品v| 人体艺术视频欧美日本| 欧美一区二区精品小视频在线| 日韩欧美精品v在线| 国产男人的电影天堂91| 日日摸夜夜添夜夜爱| 1024手机看黄色片| 精品久久久久久成人av| 久久精品熟女亚洲av麻豆精品 | 日韩视频在线欧美| www.av在线官网国产| 国产免费一级a男人的天堂| 日韩欧美三级三区| 五月玫瑰六月丁香| 日日摸夜夜添夜夜爱| 麻豆乱淫一区二区| 亚洲在线观看片| 人妻制服诱惑在线中文字幕| 汤姆久久久久久久影院中文字幕 | 久久国产乱子免费精品| 亚洲精品乱久久久久久| 69人妻影院| 最近2019中文字幕mv第一页| 国产久久久一区二区三区| 在线播放无遮挡| 两个人的视频大全免费| 少妇被粗大猛烈的视频| 欧美一区二区国产精品久久精品| 日韩av在线大香蕉| 三级毛片av免费| 日韩欧美国产在线观看| 日韩国内少妇激情av| 国产高潮美女av| 精品一区二区免费观看| 日韩成人av中文字幕在线观看| 国产美女午夜福利| 国产成人a区在线观看| 欧美激情久久久久久爽电影| 午夜福利网站1000一区二区三区| 久久这里只有精品中国| 午夜亚洲福利在线播放| 国产真实乱freesex| 国产精品久久电影中文字幕| 在线播放无遮挡| 成年女人看的毛片在线观看| av国产免费在线观看| 国产成人91sexporn| 91久久精品国产一区二区三区| 欧美激情在线99| 人妻少妇偷人精品九色| 久久久久久久久久久丰满| 国产成人午夜福利电影在线观看| 国内精品美女久久久久久| 亚洲欧美成人综合另类久久久 | 91狼人影院| 精品人妻偷拍中文字幕| 最近最新中文字幕免费大全7| 卡戴珊不雅视频在线播放| 国产精品国产三级国产专区5o | 男女国产视频网站| 国内精品宾馆在线| 亚洲美女搞黄在线观看| 欧美+日韩+精品| 婷婷色麻豆天堂久久 | 天堂网av新在线| 一级黄色大片毛片| 免费播放大片免费观看视频在线观看 | 偷拍熟女少妇极品色| 在线免费十八禁| 国产在视频线在精品| 色综合色国产| 日本与韩国留学比较| 国产成人91sexporn| 偷拍熟女少妇极品色| 成人av在线播放网站| av线在线观看网站| 亚洲,欧美,日韩| 午夜日本视频在线| 麻豆国产97在线/欧美| 搡老妇女老女人老熟妇| 最近的中文字幕免费完整| av黄色大香蕉| 精品少妇黑人巨大在线播放 | 亚洲av免费高清在线观看| 国产精品伦人一区二区| 成人午夜精彩视频在线观看| 色网站视频免费| 国产真实伦视频高清在线观看| 亚洲av免费在线观看| 丰满乱子伦码专区| av在线播放精品| av卡一久久| 久久久色成人| 精品无人区乱码1区二区| videos熟女内射| 大香蕉久久网| 91aial.com中文字幕在线观看| 久久人人爽人人爽人人片va| 国产在视频线精品| 别揉我奶头 嗯啊视频| 日韩,欧美,国产一区二区三区 | 午夜爱爱视频在线播放| 美女被艹到高潮喷水动态| 午夜福利在线观看免费完整高清在| 亚洲欧洲国产日韩| 欧美zozozo另类| 成人鲁丝片一二三区免费| 久久这里有精品视频免费| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| av在线播放精品| 干丝袜人妻中文字幕| 91狼人影院| av在线观看视频网站免费| 狂野欧美激情性xxxx在线观看| 国产美女午夜福利| 国产又色又爽无遮挡免| 亚洲精品亚洲一区二区| АⅤ资源中文在线天堂| 国产亚洲精品久久久com| 精品久久久久久久久久久久久| 波野结衣二区三区在线| 日韩制服骚丝袜av| av卡一久久| 亚洲美女搞黄在线观看| 国产黄片美女视频| 国产不卡一卡二| av国产久精品久网站免费入址| 欧美另类亚洲清纯唯美| 亚洲丝袜综合中文字幕| 在线观看66精品国产| 2021天堂中文幕一二区在线观| 欧美精品国产亚洲| 午夜免费激情av| av线在线观看网站| 国产精品.久久久| 亚洲精品日韩av片在线观看| 中文字幕熟女人妻在线| 非洲黑人性xxxx精品又粗又长| 美女被艹到高潮喷水动态| 日本免费在线观看一区| 国产精品久久久久久精品电影| av在线播放精品| 特大巨黑吊av在线直播| 国产免费男女视频| 中文天堂在线官网| 亚洲av成人av| 国产黄色视频一区二区在线观看 | 国产精品三级大全| 午夜福利在线观看吧| 内地一区二区视频在线| 少妇熟女aⅴ在线视频| 免费观看a级毛片全部| 免费大片18禁| 青春草国产在线视频| 午夜免费男女啪啪视频观看| 亚洲激情五月婷婷啪啪| 色哟哟·www| 久久久国产成人免费| 丝袜喷水一区| 免费av不卡在线播放| 97人妻精品一区二区三区麻豆| 欧美性猛交黑人性爽| 国产黄色视频一区二区在线观看 | 高清av免费在线| 国产在视频线精品| 久久久精品大字幕| 男的添女的下面高潮视频| 日日撸夜夜添| 丰满少妇做爰视频| 亚洲精品乱码久久久v下载方式| 久久鲁丝午夜福利片| av国产免费在线观看| 我要看日韩黄色一级片| 18禁裸乳无遮挡免费网站照片| 久久精品国产99精品国产亚洲性色| 国产伦一二天堂av在线观看| 亚洲精品,欧美精品| 嫩草影院新地址| 如何舔出高潮| 亚洲成色77777| 亚洲精品色激情综合| 亚洲国产精品成人久久小说| 国产又色又爽无遮挡免| 黄片无遮挡物在线观看| 国产高清有码在线观看视频| 亚洲久久久久久中文字幕| 中文亚洲av片在线观看爽| 亚洲av免费在线观看| 亚洲人与动物交配视频| 日本熟妇午夜| 亚洲怡红院男人天堂| 男插女下体视频免费在线播放| 一级毛片aaaaaa免费看小| 午夜日本视频在线| 免费观看性生交大片5| 内地一区二区视频在线| 国产午夜福利久久久久久| 国产高清视频在线观看网站| 99久久精品国产国产毛片| 九九爱精品视频在线观看| 99热精品在线国产| 看十八女毛片水多多多| 免费看a级黄色片| 九九热线精品视视频播放| 色视频www国产| www.av在线官网国产| 亚洲综合色惰| 国产精品美女特级片免费视频播放器| 色5月婷婷丁香| 国产伦精品一区二区三区视频9| 男插女下体视频免费在线播放| 国产高清三级在线| 欧美又色又爽又黄视频| 欧美zozozo另类| 久久人妻av系列| 欧美精品一区二区大全| 在线观看美女被高潮喷水网站| 欧美极品一区二区三区四区| 插逼视频在线观看| 精品久久久久久久久av| 干丝袜人妻中文字幕| 波野结衣二区三区在线| 国产亚洲91精品色在线| 美女高潮的动态| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 午夜激情福利司机影院| 亚洲国产精品合色在线| 99久久中文字幕三级久久日本| 国产淫片久久久久久久久| 欧美人与善性xxx| 99久久精品热视频| 国产淫片久久久久久久久| 欧美日本视频| 搡女人真爽免费视频火全软件| 午夜福利视频1000在线观看| 国产成人午夜福利电影在线观看| www日本黄色视频网| 一边摸一边抽搐一进一小说| av专区在线播放| 国产av在哪里看| 在线观看美女被高潮喷水网站| 国产大屁股一区二区在线视频| 国产熟女欧美一区二区| 亚洲av成人av| 免费在线观看成人毛片| 九草在线视频观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美97在线视频| 在线观看av片永久免费下载| 丝袜美腿在线中文| 99热精品在线国产| 国产一区二区三区av在线| 久久99热这里只频精品6学生 | 欧美日韩一区二区视频在线观看视频在线 | 亚洲四区av| 国产伦理片在线播放av一区| 久久久久免费精品人妻一区二区| 夜夜看夜夜爽夜夜摸| 黄片无遮挡物在线观看| 日日干狠狠操夜夜爽| 观看免费一级毛片| 日日啪夜夜撸| 国产成人a∨麻豆精品| 久久久久久久亚洲中文字幕| 亚洲精品,欧美精品| 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 亚洲综合色惰| a级毛色黄片| 夫妻性生交免费视频一级片| 久久久久久国产a免费观看| 中文字幕制服av| 午夜福利高清视频| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 日韩精品有码人妻一区| 国产黄色视频一区二区在线观看 | 免费搜索国产男女视频| 天堂网av新在线| 婷婷色综合大香蕉| 国产老妇伦熟女老妇高清| 亚洲最大成人中文| 高清av免费在线| 亚洲精品一区蜜桃| 亚洲精品自拍成人| 好男人视频免费观看在线| 麻豆一二三区av精品| 国产伦精品一区二区三区视频9| 亚洲成人久久爱视频| 色网站视频免费| 国产精品日韩av在线免费观看| 中文字幕制服av| 婷婷色麻豆天堂久久 | 国产av码专区亚洲av| 淫秽高清视频在线观看| 丝袜喷水一区| 成人三级黄色视频| 男人和女人高潮做爰伦理| 国产成人一区二区在线| 国产一区亚洲一区在线观看| 成人二区视频| 一区二区三区高清视频在线| 91久久精品国产一区二区成人| 看片在线看免费视频| 热99在线观看视频| 麻豆精品久久久久久蜜桃| 在线天堂最新版资源| 超碰av人人做人人爽久久| 国产日韩欧美在线精品| 日韩精品有码人妻一区| 在线免费观看不下载黄p国产| 三级男女做爰猛烈吃奶摸视频| 欧美一区二区国产精品久久精品| 内地一区二区视频在线| 免费看a级黄色片| 2021少妇久久久久久久久久久| 日本熟妇午夜| 国产亚洲午夜精品一区二区久久 | 欧美一区二区亚洲| 又黄又爽又刺激的免费视频.| 亚洲精品成人久久久久久| 成人性生交大片免费视频hd| 黄色一级大片看看| 国产亚洲午夜精品一区二区久久 | 搡老妇女老女人老熟妇| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 成人鲁丝片一二三区免费| 日本免费一区二区三区高清不卡| 丰满少妇做爰视频| 欧美激情国产日韩精品一区| 欧美精品国产亚洲| 乱码一卡2卡4卡精品| 少妇人妻一区二区三区视频| 国产麻豆成人av免费视频| 最近中文字幕2019免费版| or卡值多少钱| 久久人妻av系列| 日本午夜av视频| 免费在线观看成人毛片| 三级国产精品欧美在线观看| 好男人视频免费观看在线| 中国美白少妇内射xxxbb| 国产一区二区在线av高清观看| 别揉我奶头 嗯啊视频| 久久99精品国语久久久| 午夜免费激情av| 三级经典国产精品| 乱人视频在线观看| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 在现免费观看毛片| 欧美区成人在线视频| 日日干狠狠操夜夜爽| 波多野结衣巨乳人妻| 精品一区二区免费观看| 久久久国产成人精品二区| 91午夜精品亚洲一区二区三区| 久久韩国三级中文字幕| 十八禁国产超污无遮挡网站| 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| 国产成人freesex在线| 丰满乱子伦码专区| 欧美97在线视频| 亚洲伊人久久精品综合 | 国产伦一二天堂av在线观看| 久久久久久国产a免费观看| 精品人妻一区二区三区麻豆| 日韩欧美精品v在线| 亚洲av熟女| 看片在线看免费视频| 五月玫瑰六月丁香| 又爽又黄a免费视频| 久久国产乱子免费精品| 最近2019中文字幕mv第一页| 精品国产一区二区三区久久久樱花 | 神马国产精品三级电影在线观看| 高清视频免费观看一区二区 | 国产免费视频播放在线视频 | 特大巨黑吊av在线直播| 岛国毛片在线播放| 日韩三级伦理在线观看| 欧美成人免费av一区二区三区| 中文欧美无线码| 国产美女午夜福利| 午夜福利网站1000一区二区三区| 久久国内精品自在自线图片| 国产精品久久久久久久电影| 精品久久久久久久末码| 不卡视频在线观看欧美| 欧美一区二区亚洲| 国产精品野战在线观看| 免费av毛片视频| 色尼玛亚洲综合影院| 九九久久精品国产亚洲av麻豆| 午夜福利成人在线免费观看| 2021天堂中文幕一二区在线观| 亚洲欧美一区二区三区国产| 一区二区三区高清视频在线| 国产精品久久久久久精品电影小说 | 少妇高潮的动态图| 国产黄色小视频在线观看| 久久亚洲国产成人精品v| 国产亚洲最大av| 国产精品野战在线观看| 有码 亚洲区| 啦啦啦观看免费观看视频高清| 亚洲一级一片aⅴ在线观看| 丝袜喷水一区| 黑人高潮一二区| 久久久a久久爽久久v久久| 色综合亚洲欧美另类图片| 一卡2卡三卡四卡精品乱码亚洲| 美女大奶头视频| 欧美3d第一页| 国产精品国产三级专区第一集| 最近的中文字幕免费完整| 欧美精品一区二区大全| 国产真实乱freesex| 国产爱豆传媒在线观看| av卡一久久| 国产探花在线观看一区二区| 国产高清视频在线观看网站| 99热这里只有是精品50| 亚洲av成人av| 联通29元200g的流量卡| 又粗又硬又长又爽又黄的视频| 免费人成在线观看视频色| www.色视频.com| 99久久成人亚洲精品观看| 我的女老师完整版在线观看| 日本一本二区三区精品| 亚洲国产欧洲综合997久久,| 亚洲av成人精品一区久久| 久久久久久久久久黄片| 国产成人精品一,二区| 色吧在线观看| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线播| 国产不卡一卡二| 一卡2卡三卡四卡精品乱码亚洲| 26uuu在线亚洲综合色| 午夜精品国产一区二区电影 | 成年免费大片在线观看| av免费观看日本| 国产91av在线免费观看| 99热网站在线观看| 男女那种视频在线观看| 亚洲18禁久久av| 青春草视频在线免费观看| 久久久久久久午夜电影| 国产伦理片在线播放av一区| 综合色丁香网| 看免费成人av毛片| 亚洲婷婷狠狠爱综合网| 美女脱内裤让男人舔精品视频| av国产久精品久网站免费入址| 亚洲电影在线观看av| 精品一区二区三区视频在线| 精品人妻偷拍中文字幕| 午夜激情欧美在线| 国产真实伦视频高清在线观看| 免费看日本二区| 男女下面进入的视频免费午夜| av线在线观看网站| 日韩三级伦理在线观看| 联通29元200g的流量卡| 综合色av麻豆| 伦精品一区二区三区| 91久久精品国产一区二区成人| 免费av不卡在线播放| 国产精华一区二区三区| 寂寞人妻少妇视频99o| 最近的中文字幕免费完整| 男女啪啪激烈高潮av片| 老司机影院成人| 91午夜精品亚洲一区二区三区| 韩国av在线不卡| 亚洲精品乱码久久久久久按摩| 麻豆乱淫一区二区| 日韩欧美精品v在线| 18禁动态无遮挡网站| 国产成人91sexporn| 精品久久久噜噜| or卡值多少钱| 国产成人一区二区在线| 国产精品女同一区二区软件| 欧美高清成人免费视频www| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 长腿黑丝高跟| 免费观看精品视频网站| 久久久久久久午夜电影| 国产免费男女视频| 一二三四中文在线观看免费高清| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品成人久久小说| 免费看a级黄色片| 九九爱精品视频在线观看| 美女黄网站色视频| 看片在线看免费视频| av在线天堂中文字幕| 国语自产精品视频在线第100页| 国内精品宾馆在线| 中文字幕av成人在线电影| 免费看av在线观看网站| 国产成人91sexporn| 国产精品野战在线观看| 国产成人福利小说| 欧美人与善性xxx| 好男人视频免费观看在线| 国产高潮美女av| 免费看日本二区| 国产精品av视频在线免费观看| 啦啦啦韩国在线观看视频| 噜噜噜噜噜久久久久久91| 国产精品乱码一区二三区的特点| av视频在线观看入口| 全区人妻精品视频| 少妇人妻精品综合一区二区| 国产 一区精品| 尤物成人国产欧美一区二区三区| 亚洲av成人av| 国产成人一区二区在线| 午夜福利在线观看吧| ponron亚洲| 又爽又黄无遮挡网站| 日本wwww免费看| 成人欧美大片| 午夜亚洲福利在线播放| 熟女电影av网| 少妇人妻精品综合一区二区| 深夜a级毛片| av在线老鸭窝| 乱人视频在线观看| 国产精品1区2区在线观看.| 亚洲国产精品成人久久小说| 只有这里有精品99| 国产久久久一区二区三区| 九色成人免费人妻av| 久久久久久久久久成人| 国模一区二区三区四区视频| 国产熟女欧美一区二区| 成年版毛片免费区| 久久久国产成人免费| 97人妻精品一区二区三区麻豆| 一级爰片在线观看| 免费大片18禁| 精品国产露脸久久av麻豆 | 国产黄色小视频在线观看|