• <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)研究*
    中文字幕熟女人妻在线| 美女高潮喷水抽搐中文字幕| www日本在线高清视频| 99久久国产精品久久久| 亚洲一码二码三码区别大吗| avwww免费| 一区二区三区高清视频在线| 中文字幕久久专区| 非洲黑人性xxxx精品又粗又长| 日韩高清综合在线| 俄罗斯特黄特色一大片| 日韩 欧美 亚洲 中文字幕| 久久久久久人人人人人| 婷婷亚洲欧美| 午夜久久久久精精品| 国产精品美女特级片免费视频播放器 | 精品一区二区三区av网在线观看| 日本一区二区免费在线视频| 亚洲国产精品成人综合色| 日韩欧美三级三区| 日韩成人在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 国产熟女xx| 国产在线精品亚洲第一网站| 欧美3d第一页| 好男人在线观看高清免费视频| 岛国在线免费视频观看| 久久精品综合一区二区三区| 午夜福利欧美成人| 久久久久久免费高清国产稀缺| 日韩国内少妇激情av| 制服人妻中文乱码| 日本撒尿小便嘘嘘汇集6| 亚洲五月婷婷丁香| 欧美日韩亚洲综合一区二区三区_| 精品电影一区二区在线| 国产免费男女视频| 精品久久久久久久末码| 色综合站精品国产| 1024香蕉在线观看| 亚洲乱码一区二区免费版| 女同久久另类99精品国产91| 两个人看的免费小视频| 国产视频内射| tocl精华| 午夜亚洲福利在线播放| 九九热线精品视视频播放| 男女床上黄色一级片免费看| 丰满的人妻完整版| 国产高清视频在线播放一区| 亚洲va日本ⅴa欧美va伊人久久| 99精品在免费线老司机午夜| 色av中文字幕| 好男人在线观看高清免费视频| 欧美日韩亚洲综合一区二区三区_| 丰满的人妻完整版| 精品久久久久久,| 999精品在线视频| 欧美不卡视频在线免费观看 | 亚洲精品在线美女| 国产伦在线观看视频一区| 国产久久久一区二区三区| 国产麻豆成人av免费视频| 亚洲五月婷婷丁香| 99久久无色码亚洲精品果冻| 久久香蕉国产精品| 在线永久观看黄色视频| 欧美+亚洲+日韩+国产| 亚洲中文字幕日韩| 日韩国内少妇激情av| 亚洲自偷自拍图片 自拍| 真人做人爱边吃奶动态| 欧美大码av| 亚洲第一电影网av| 一级作爱视频免费观看| 免费人成视频x8x8入口观看| 久久中文字幕人妻熟女| 久久久久久国产a免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 最近视频中文字幕2019在线8| 黄片小视频在线播放| 亚洲美女视频黄频| 国产精品香港三级国产av潘金莲| 一a级毛片在线观看| 日韩欧美精品v在线| 两性夫妻黄色片| 国产私拍福利视频在线观看| 人人妻人人澡欧美一区二区| 12—13女人毛片做爰片一| 亚洲中文字幕一区二区三区有码在线看 | 亚洲九九香蕉| av在线天堂中文字幕| 亚洲熟妇熟女久久| 亚洲片人在线观看| 美女免费视频网站| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲欧美98| 1024视频免费在线观看| 久久精品国产亚洲av高清一级| 91成年电影在线观看| 视频区欧美日本亚洲| 国产一区二区激情短视频| 欧美午夜高清在线| 国产av一区二区精品久久| 国产真人三级小视频在线观看| 一区二区三区高清视频在线| 国产精品久久久久久久电影 | 岛国在线免费视频观看| 一个人观看的视频www高清免费观看 | 国内精品久久久久精免费| 国产人伦9x9x在线观看| 午夜福利成人在线免费观看| e午夜精品久久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 俺也久久电影网| 国产亚洲av嫩草精品影院| 给我免费播放毛片高清在线观看| 久久天躁狠狠躁夜夜2o2o| 麻豆av在线久日| 久久亚洲真实| 这个男人来自地球电影免费观看| 久久亚洲真实| 嫁个100分男人电影在线观看| 真人一进一出gif抽搐免费| av福利片在线| 久久久精品国产亚洲av高清涩受| 操出白浆在线播放| 女人被狂操c到高潮| 首页视频小说图片口味搜索| 中出人妻视频一区二区| 中出人妻视频一区二区| 国产v大片淫在线免费观看| 午夜福利成人在线免费观看| 天堂影院成人在线观看| 亚洲av成人精品一区久久| 久久久久久久久免费视频了| 精品一区二区三区av网在线观看| 久久 成人 亚洲| 一级毛片女人18水好多| 国产激情偷乱视频一区二区| 韩国av一区二区三区四区| 一本大道久久a久久精品| 精品高清国产在线一区| 欧美日本视频| 免费在线观看视频国产中文字幕亚洲| 白带黄色成豆腐渣| 男人舔女人的私密视频| 亚洲 国产 在线| 亚洲真实伦在线观看| 亚洲av成人不卡在线观看播放网| 美女午夜性视频免费| 亚洲片人在线观看| 丁香欧美五月| 国产精品,欧美在线| 老汉色av国产亚洲站长工具| 欧美性猛交╳xxx乱大交人| 亚洲精品中文字幕一二三四区| a在线观看视频网站| 欧美日韩一级在线毛片| 可以免费在线观看a视频的电影网站| 欧美另类亚洲清纯唯美| 亚洲人成网站高清观看| 一本综合久久免费| 一进一出抽搐gif免费好疼| 一a级毛片在线观看| 哪里可以看免费的av片| 亚洲av成人av| 中文字幕熟女人妻在线| 精品久久久久久久末码| 国产成人aa在线观看| 99久久国产精品久久久| 欧美另类亚洲清纯唯美| 两个人看的免费小视频| 精品高清国产在线一区| 欧美高清成人免费视频www| 色av中文字幕| 精品熟女少妇八av免费久了| 亚洲五月天丁香| 男插女下体视频免费在线播放| 99在线视频只有这里精品首页| 男人舔奶头视频| 国产精品影院久久| 亚洲专区中文字幕在线| 国产成人系列免费观看| 亚洲精品色激情综合| 欧美日韩乱码在线| 人成视频在线观看免费观看| 成年人黄色毛片网站| 女人被狂操c到高潮| 国产av麻豆久久久久久久| 在线十欧美十亚洲十日本专区| 激情在线观看视频在线高清| 不卡一级毛片| 少妇人妻一区二区三区视频| ponron亚洲| 精品人妻1区二区| 国产三级在线视频| 久久这里只有精品中国| 午夜免费观看网址| 午夜福利免费观看在线| 91av网站免费观看| 99在线视频只有这里精品首页| 在线观看免费日韩欧美大片| 国产视频内射| 国产又黄又爽又无遮挡在线| 欧美乱码精品一区二区三区| 女同久久另类99精品国产91| 欧美+亚洲+日韩+国产| 免费在线观看完整版高清| 亚洲av第一区精品v没综合| 日韩欧美免费精品| cao死你这个sao货| 在线观看免费日韩欧美大片| 午夜精品久久久久久毛片777| 亚洲一区高清亚洲精品| 级片在线观看| 一级毛片高清免费大全| 亚洲va日本ⅴa欧美va伊人久久| 欧美成人午夜精品| 一区二区三区高清视频在线| 天天躁狠狠躁夜夜躁狠狠躁| 精华霜和精华液先用哪个| 中文字幕人成人乱码亚洲影| 男女做爰动态图高潮gif福利片| 精品乱码久久久久久99久播| 麻豆国产97在线/欧美 | 成人午夜高清在线视频| 一本一本综合久久| 久久精品国产99精品国产亚洲性色| 国产久久久一区二区三区| 亚洲av成人精品一区久久| 在线国产一区二区在线| 国产精品1区2区在线观看.| ponron亚洲| 国产野战对白在线观看| 午夜免费观看网址| 俄罗斯特黄特色一大片| 免费看十八禁软件| 老司机午夜十八禁免费视频| 精品久久蜜臀av无| 日本在线视频免费播放| 99久久99久久久精品蜜桃| 国内精品一区二区在线观看| 亚洲 欧美 日韩 在线 免费| 18禁黄网站禁片免费观看直播| 日韩三级视频一区二区三区| 亚洲中文av在线| 国产成人精品久久二区二区91| 久9热在线精品视频| 中文在线观看免费www的网站 | 免费高清视频大片| 青草久久国产| 亚洲国产中文字幕在线视频| 久久精品亚洲精品国产色婷小说| av中文乱码字幕在线| 19禁男女啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 欧美久久黑人一区二区| 三级国产精品欧美在线观看 | 麻豆久久精品国产亚洲av| 欧美一区二区精品小视频在线| 亚洲成人中文字幕在线播放| 午夜福利免费观看在线| 99在线视频只有这里精品首页| 成人18禁在线播放| 欧美日本视频| 老熟妇仑乱视频hdxx| 国产精品一区二区三区四区久久| 欧美成人一区二区免费高清观看 | 制服诱惑二区| 又爽又黄无遮挡网站| 欧美日韩亚洲国产一区二区在线观看| xxx96com| 亚洲国产欧美网| 老司机靠b影院| 色噜噜av男人的天堂激情| 国产午夜精品久久久久久| 亚洲精华国产精华精| 国产日本99.免费观看| 日本三级黄在线观看| 一区二区三区激情视频| 999久久久国产精品视频| 亚洲精品久久成人aⅴ小说| 久久久久久大精品| 日韩高清综合在线| 精品久久蜜臀av无| 亚洲第一电影网av| 91在线观看av| 国产精品一区二区精品视频观看| 久久热在线av| 精品国产乱子伦一区二区三区| 少妇熟女aⅴ在线视频| 日韩欧美 国产精品| 亚洲av成人不卡在线观看播放网| 岛国视频午夜一区免费看| 免费搜索国产男女视频| 亚洲国产精品久久男人天堂| 久久久久国产精品人妻aⅴ院| 午夜福利18| 天天添夜夜摸| 97超级碰碰碰精品色视频在线观看| 色综合婷婷激情| 亚洲 国产 在线| 亚洲av日韩精品久久久久久密| 好男人在线观看高清免费视频| 日日夜夜操网爽| 一区二区三区激情视频| 日韩精品免费视频一区二区三区| 久久中文字幕一级| 夜夜夜夜夜久久久久| 两个人看的免费小视频| 久久久久免费精品人妻一区二区| 日韩免费av在线播放| 宅男免费午夜| 99在线视频只有这里精品首页| 50天的宝宝边吃奶边哭怎么回事| 日本成人三级电影网站| 国产精品一及| 国产成人影院久久av| √禁漫天堂资源中文www| 国产成人精品久久二区二区免费| 亚洲 欧美 日韩 在线 免费| 婷婷六月久久综合丁香| 国产高清有码在线观看视频 | 精品欧美一区二区三区在线| 欧美日本视频| 亚洲av成人不卡在线观看播放网| 亚洲国产欧美人成| 免费在线观看亚洲国产| 9191精品国产免费久久| 久久人人精品亚洲av| 我要搜黄色片| x7x7x7水蜜桃| 精品国产乱子伦一区二区三区| 国产精品99久久99久久久不卡| 久久精品夜夜夜夜夜久久蜜豆 | 听说在线观看完整版免费高清| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 精品一区二区三区视频在线观看免费| 欧美黄色片欧美黄色片| 男插女下体视频免费在线播放| 久久久精品大字幕| 国产精品综合久久久久久久免费| 亚洲男人天堂网一区| 黄片大片在线免费观看| 日本一本二区三区精品| 亚洲成人久久性| 亚洲国产欧美人成| 精品久久久久久成人av| 9191精品国产免费久久| 午夜福利在线在线| 少妇裸体淫交视频免费看高清 | 两个人看的免费小视频| 亚洲熟女毛片儿| 一进一出好大好爽视频| 亚洲精品久久国产高清桃花| 最近最新中文字幕大全电影3| 亚洲电影在线观看av| 国产成年人精品一区二区| 国产欧美日韩一区二区精品| 成人欧美大片| 亚洲乱码一区二区免费版| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 国产精品一区二区三区四区免费观看 | 亚洲最大成人中文| 午夜福利在线在线| 99国产综合亚洲精品| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 久久精品国产清高在天天线| 欧美乱妇无乱码| 国产精品电影一区二区三区| 欧美一区二区精品小视频在线| 天天添夜夜摸| 久久草成人影院| 久热爱精品视频在线9| 岛国视频午夜一区免费看| av福利片在线| 亚洲专区字幕在线| 久久亚洲精品不卡| 一本久久中文字幕| 真人做人爱边吃奶动态| 欧美av亚洲av综合av国产av| 51午夜福利影视在线观看| √禁漫天堂资源中文www| 动漫黄色视频在线观看| 国产亚洲欧美在线一区二区| 97超级碰碰碰精品色视频在线观看| 在线十欧美十亚洲十日本专区| 国产蜜桃级精品一区二区三区| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 成人国产一区最新在线观看| 制服丝袜大香蕉在线| svipshipincom国产片| 精品国产亚洲在线| 久久久久久久午夜电影| e午夜精品久久久久久久| 岛国在线观看网站| 99热只有精品国产| 国产精品美女特级片免费视频播放器 | 可以在线观看的亚洲视频| 成人精品一区二区免费| 亚洲成人精品中文字幕电影| 亚洲五月天丁香| 脱女人内裤的视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久天躁狠狠躁夜夜2o2o| 怎么达到女性高潮| 2021天堂中文幕一二区在线观| 亚洲国产欧美一区二区综合| 午夜精品一区二区三区免费看| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 国产av麻豆久久久久久久| 日韩欧美精品v在线| 一级片免费观看大全| 亚洲人成网站高清观看| 午夜免费成人在线视频| 一级片免费观看大全| 国产激情偷乱视频一区二区| 国产成人系列免费观看| a在线观看视频网站| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 十八禁人妻一区二区| 精品久久久久久久人妻蜜臀av| 国产一区在线观看成人免费| www.www免费av| 真人做人爱边吃奶动态| 亚洲成人中文字幕在线播放| 欧美在线一区亚洲| 久久久久久久久免费视频了| 亚洲av成人不卡在线观看播放网| 免费看日本二区| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 久久久精品国产亚洲av高清涩受| 久久久久国产精品人妻aⅴ院| 88av欧美| 首页视频小说图片口味搜索| 在线观看免费视频日本深夜| 视频区欧美日本亚洲| 国产精品 欧美亚洲| 每晚都被弄得嗷嗷叫到高潮| 叶爱在线成人免费视频播放| 成人18禁在线播放| 波多野结衣高清无吗| av超薄肉色丝袜交足视频| 黄色毛片三级朝国网站| 欧美在线黄色| 久久香蕉激情| 手机成人av网站| 亚洲国产欧美人成| 日本三级黄在线观看| 91字幕亚洲| 欧美成人一区二区免费高清观看 | 国产在线精品亚洲第一网站| 最好的美女福利视频网| 久久久久九九精品影院| www日本在线高清视频| 天堂√8在线中文| 成熟少妇高潮喷水视频| 午夜视频精品福利| 69av精品久久久久久| 国产激情偷乱视频一区二区| 国产免费av片在线观看野外av| 久久久精品国产亚洲av高清涩受| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产精品麻豆| www.熟女人妻精品国产| 波多野结衣高清无吗| 午夜免费激情av| 欧美绝顶高潮抽搐喷水| 99在线人妻在线中文字幕| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 亚洲精品久久国产高清桃花| 超碰成人久久| 亚洲精品粉嫩美女一区| svipshipincom国产片| 午夜福利在线在线| 亚洲成av人片免费观看| www日本在线高清视频| 欧美一区二区国产精品久久精品 | 成人国产综合亚洲| videosex国产| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 女警被强在线播放| av福利片在线观看| 欧美日韩亚洲综合一区二区三区_| 国内精品久久久久精免费| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 国产午夜福利久久久久久| 国产精品久久久久久人妻精品电影| 中文字幕高清在线视频| 国产精品免费一区二区三区在线| 久久久久久亚洲精品国产蜜桃av| 老汉色∧v一级毛片| 91国产中文字幕| 亚洲成人免费电影在线观看| 18禁美女被吸乳视频| 免费一级毛片在线播放高清视频| 免费av毛片视频| 色综合欧美亚洲国产小说| 亚洲av成人不卡在线观看播放网| 久久性视频一级片| 精品福利观看| 91av网站免费观看| 一区二区三区国产精品乱码| 久久久久久久久久黄片| 欧美中文综合在线视频| 亚洲成人精品中文字幕电影| 亚洲成人久久爱视频| 999久久久精品免费观看国产| 夜夜爽天天搞| 亚洲最大成人中文| 亚洲九九香蕉| 看片在线看免费视频| 男女床上黄色一级片免费看| 日韩欧美在线二视频| 美女大奶头视频| 99久久精品国产亚洲精品| 日韩欧美 国产精品| 国产精华一区二区三区| 亚洲av五月六月丁香网| 性色av乱码一区二区三区2| 波多野结衣巨乳人妻| 天堂动漫精品| 级片在线观看| av欧美777| 特大巨黑吊av在线直播| av天堂在线播放| 一级作爱视频免费观看| 12—13女人毛片做爰片一| 嫁个100分男人电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产亚洲av高清一级| 中文资源天堂在线| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 欧美午夜高清在线| 欧美中文综合在线视频| 国产一区二区在线观看日韩 | 91在线观看av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站高清观看| 成人高潮视频无遮挡免费网站| 一个人免费在线观看的高清视频| 日日夜夜操网爽| 十八禁网站免费在线| 午夜精品一区二区三区免费看| 国产精品免费一区二区三区在线| 国产精品永久免费网站| 久久亚洲精品不卡| 精品久久久久久成人av| 亚洲国产看品久久| 老司机午夜福利在线观看视频| 国产三级在线视频| 桃色一区二区三区在线观看| 国产精品影院久久| aaaaa片日本免费| 国产三级中文精品| 午夜激情福利司机影院| 又粗又爽又猛毛片免费看| 免费在线观看黄色视频的| 少妇粗大呻吟视频| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 午夜精品在线福利| 亚洲中文字幕一区二区三区有码在线看 | 国产成人系列免费观看| 国产激情久久老熟女| 免费av毛片视频| 久久人妻av系列| 五月玫瑰六月丁香| 国产激情欧美一区二区| 视频区欧美日本亚洲| 亚洲成av人片在线播放无| 99国产综合亚洲精品| 国产亚洲精品一区二区www| 叶爱在线成人免费视频播放| 免费在线观看视频国产中文字幕亚洲| 国产黄色小视频在线观看| 757午夜福利合集在线观看| 大型av网站在线播放| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 日韩中文字幕欧美一区二区| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播放欧美日韩| 国产成人影院久久av| 亚洲 欧美 日韩 在线 免费| 久久午夜亚洲精品久久| 我的老师免费观看完整版| 精品熟女少妇八av免费久了| 久久精品人妻少妇| 国产视频内射| 99re在线观看精品视频| 日本三级黄在线观看| 亚洲激情在线av| 国产又黄又爽又无遮挡在线| 天堂动漫精品| 色综合婷婷激情| 日韩欧美免费精品| 日韩精品中文字幕看吧| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 给我免费播放毛片高清在线观看| 免费在线观看影片大全网站| 黄色片一级片一级黄色片|