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

    孤立子內(nèi)波存在下的聲傳播仿真研究

    2019-07-31 09:28:16邢傳璽劉文博劉佳琪姜思源岳露露
    關(guān)鍵詞:內(nèi)波聲線聲壓

    邢傳璽,宋 揚(yáng),劉文博,劉佳琪,姜思源,崔 琳,岳露露

    (1.云南民族大學(xué) 電氣信息工程學(xué)院, 云南 昆明 650500; 2.哈爾濱工程大學(xué) 水聲技術(shù)重點(diǎn)實(shí)驗(yàn)室,黑龍江 哈爾濱 150001;3.哈爾濱工程大學(xué) 水聲工程學(xué)院, 黑龍江 哈爾濱 150001)

    海洋中由于海水密度垂向分布不均勻會(huì)導(dǎo)致的海水密度垂直梯度突變,產(chǎn)生密度躍層.此密度躍層內(nèi)界面,由于水流擾動(dòng)所引起的波動(dòng)被稱作海洋內(nèi)波,孤立子內(nèi)波通常是指其中的非線性大振幅波動(dòng).孤立子內(nèi)波一般是由強(qiáng)度比較大的流通過海底陡坡、海脊等形成的,其形狀且和波速在相互碰撞后只有微弱的變化或者保持不變[1-5].

    孤立子內(nèi)波又叫做非線性內(nèi)波[2],通常具有非線性、大振幅、短周期的特點(diǎn),經(jīng)常以3種方式出現(xiàn),分別為單個(gè)孤立波、多個(gè)孤立波形成波包、多個(gè)波包形成波群,其特征寬度從幾百米到幾千米、周期為10到30分鐘.孤立子內(nèi)波最常以波群方式傳播,其中多個(gè)波包間距離可達(dá)到幾十至一百千米,特征波長(zhǎng)最大可達(dá)幾百千米.由于其在界面處振幅較大,孤立子內(nèi)波在傳播過程中不但會(huì)引起海水強(qiáng)烈的垂向混合、影響海洋動(dòng)植物的生存環(huán)境,而且在傳播過程中引起海面劇烈輻聚和突發(fā)性的強(qiáng)流.孤立子內(nèi)波又?jǐn)y帶巨大的能量,會(huì)對(duì)水下聲吶、海洋工程、海洋捕撈、石油鉆井等海面、海底設(shè)施構(gòu)成重大威脅.有研究表明,孤立子內(nèi)波的傳播過程中會(huì)產(chǎn)生假潮,對(duì)海岸港口會(huì)產(chǎn)生極大的破壞性[3-5].

    由于孤立子內(nèi)波在海洋上帶來的威脅僅次于臺(tái)風(fēng),國(guó)內(nèi)外學(xué)者紛紛開始展開了對(duì)孤立子內(nèi)波的研究.2005年Sugisoto[6]研究了孤立子內(nèi)波對(duì)質(zhì)量、動(dòng)量、能量的傳遞,指出了孤立子內(nèi)波對(duì)聲場(chǎng)能量的作用.LUO[7]在2008年,利用射線理論的方法,研究了孤立子內(nèi)波對(duì)聲線的聚焦和散焦作用.2009年,Katsnelson Boris[8]研究了淺水孤立子內(nèi)波存在下聲場(chǎng)的時(shí)空波動(dòng).國(guó)內(nèi)的學(xué)者主要研究集中在中國(guó)南部海域,2004年,陳守虎[9]等對(duì)中國(guó)南部海域內(nèi)波的特點(diǎn)和內(nèi)波的時(shí)空相關(guān)性進(jìn)行了研究,并給出了內(nèi)波的簡(jiǎn)正波模態(tài)分析結(jié)果.劉進(jìn)忠[10]詳細(xì)分析了內(nèi)波對(duì)聲場(chǎng)傳播損失、傳播時(shí)間、簡(jiǎn)正波系數(shù)等的影響,而且還提出了一種利用簡(jiǎn)正波群速度來反演孤立子內(nèi)波的方法,這十分有望應(yīng)用到孤立子內(nèi)波檢測(cè)技術(shù)當(dāng)中.2008年,季桂花[11]等研究了淺海各種內(nèi)波對(duì)匹配場(chǎng)時(shí)間相關(guān)的影響,提出了內(nèi)波會(huì)使溫躍層局部大幅度下降,影響水下設(shè)備的檢測(cè)定位性能.2009年,哈爾濱工程大學(xué)的馬樹青[12]進(jìn)一步進(jìn)行了孤立子內(nèi)波對(duì)聲場(chǎng)聲壓傳播損失及對(duì)其對(duì)聲源定位的研究.2011年,李整林[13]等進(jìn)行了孤立子內(nèi)波存在下簡(jiǎn)正波到達(dá)時(shí)間的研究.2015年,周江濤[14]等亞研究了孤立子內(nèi)波對(duì)聲場(chǎng)結(jié)構(gòu)的影響.2016年,秦繼興[15]等對(duì)淺海能量起伏進(jìn)行了研究,并對(duì)有孤立子內(nèi)波和沒有孤立子內(nèi)波時(shí)進(jìn)行了對(duì)比.

    我國(guó)擁有廣闊的海域和海岸線,在我國(guó)沿海有很多淺海地區(qū)都存在著海洋內(nèi)波,隨著海洋工程的日益開展,海洋捕撈業(yè)的蓬勃發(fā)展,海底資源的大力開發(fā),內(nèi)波尤其是對(duì)水面、水底威脅最大的孤立子內(nèi)波(非線性內(nèi)波)就越需要去探究.聲波是探知海洋、水下定位導(dǎo)航和水下通信的有效載體,建設(shè)“數(shù)字海洋”,必須要了解、掌握孤立子內(nèi)波對(duì)聲傳播的影響,因此,本文基于COMSOL有限元方法,對(duì)孤立子內(nèi)波對(duì)聲傳播的影響進(jìn)行了仿真研究,為孤立子內(nèi)部的研究提供理論基礎(chǔ)和研究依據(jù).

    1 孤立子內(nèi)波建模

    1.1 線性內(nèi)波模型

    因?yàn)榫€性內(nèi)波具有很大的隨機(jī)性,所以一般用隨機(jī)場(chǎng)在描述.在不考慮背景流場(chǎng)的情況下,設(shè)水體滿足密度層化的條件,可知內(nèi)波的定解方程為[2]:

    (1)

    其中W為內(nèi)波水質(zhì)點(diǎn)垂向運(yùn)動(dòng)速度ω的振幅,表達(dá)式為ω=W(z)exp[i(σt-kx-ly)];k和l分別為水平坐標(biāo)(x,y)方向上的波數(shù),且χ2=k2+l2;σ表示頻率,t為時(shí)間,g為重力加速度,f為科氏參數(shù),β=df/dy;H為不變值海洋深度,N為浮力頻率表達(dá)式為:

    (2)

    引入f-平面近似,即取β=0,且設(shè)l=0,則有:

    (3)

    海洋中常見的內(nèi)波波長(zhǎng)一般都不超過1 km,因此f-平面近似可以成立[2].

    然后,設(shè)海面是沒有波動(dòng)和起伏的,W(0)=0, 引入“剛蓋近似”:

    (4)

    W(-H)=0,W(0)=0.

    再進(jìn)一步引入Boussinesq近似,內(nèi)波方程可進(jìn)一步的簡(jiǎn)化為:

    (5)

    W(-H)=0,W(0)=0.

    這就是Fjelsted內(nèi)波方程.

    若忽略地轉(zhuǎn)效應(yīng)使f=0,且滿足大尺度波動(dòng)的情況下,或者當(dāng)f2<<σ2<

    (6)

    W(-H)=0,W(0)=0,

    其中c=σ/k為波動(dòng)的相速度.

    1.2 孤立子內(nèi)波模型

    通常情況下,可通過KdV方程來求解方程的孤立子解[3].

    如果KdV方程的形式如下:

    θt-6θθx+θxxx=0,

    (7)

    其中θ為KdV方程的孤立子解,θt為其對(duì)時(shí)間t的一階偏導(dǎo),θx為其對(duì)位移x的一階偏導(dǎo),θxxx為其對(duì)位移x的三階偏導(dǎo).

    設(shè)θ(t,x)=f(u),其中u=x-ct,代入公式(7)后求解可得孤立子解:

    (8)

    如果不考慮地轉(zhuǎn)影響,引入Boussinesq近似可得到KdV方程的標(biāo)準(zhǔn)形式,求解之得到孤立子解.不同孤立子內(nèi)波模型的KdV方程可得到不同的孤立子內(nèi)波解.本文使用的模型為如圖1所示,

    此模型下,孤立子內(nèi)波的垂向位移為:

    (9)

    公式9的模型適用于實(shí)際的海洋環(huán)境,由n個(gè)孤立子內(nèi)波組成的模型,其中的W(kh,j)為簡(jiǎn)正波的模態(tài)函數(shù),Λm為孤立子內(nèi)波的振幅,Vm為其移動(dòng)速度,Δm為孤立子內(nèi)波的波包寬度.

    2 孤立子內(nèi)波對(duì)聲線軌跡的影響仿真

    2.1 孤立子內(nèi)波對(duì)海水聲速梯度的影響仿真

    孤立子內(nèi)波會(huì)改變海水的聲速,海水聲速的變化會(huì)使聲線軌跡發(fā)生變化,最后可通過聲線軌跡的改變來反映孤立子內(nèi)波對(duì)二維聲場(chǎng)的影響.

    設(shè)孤立子內(nèi)波的振幅15 m,波包寬度為50 m,在距聲源水平距離200 m處出現(xiàn)其振幅的極大值:

    (10)

    為了便于觀察設(shè)置海水為溫躍層淺海,海水深度為100 m,海面視為絕對(duì)軟邊界,海底視為絕對(duì)硬邊界,海水中存在負(fù)聲速梯度,圖2為存在孤立子內(nèi)波時(shí)的海水聲速分布.

    2.2 孤立子內(nèi)波對(duì)聲線軌跡的影響仿真研究

    利用存在孤立子內(nèi)波時(shí)的聲速分布進(jìn)行聲線軌跡仿真,設(shè)聲場(chǎng)中存在一點(diǎn)聲源,聲源位于10 m處,聲源頻率為100 Hz.在COMSOL中把物理場(chǎng)設(shè)為射線聲學(xué),研究設(shè)為聲線軌跡,仿真得到有無(wú)孤立子內(nèi)波時(shí)聲線軌跡圖如圖3所示,為便于討論孤立子內(nèi)波對(duì)聲線軌跡的影響,仿真中省略了海底反射等因素.

    從圖3中可以看出,孤立子內(nèi)波可以明顯的改變聲線的傳播軌跡,水平200~400 m區(qū)域是孤立子內(nèi)波影響的主要區(qū)域,其聲線軌跡變得稀疏,400 m后的聲線軌跡則變得十分密集.可近似的認(rèn)為孤立子內(nèi)波使200~400 m區(qū)域的能量減少,使400 m后區(qū)域能量增加.

    3 孤立子內(nèi)波對(duì)聲場(chǎng)的影響仿真

    3.1 聲源頻率不同時(shí)孤立子內(nèi)波對(duì)聲傳播的影響仿真

    設(shè)聲源深度10 m,接收深度10 m,分別計(jì)算當(dāng)聲源發(fā)射頻率為50 Hz、100 Hz、150 Hz、200 Hz時(shí)的傳播損失,其中聲源頻率為50 Hz和250 Hz時(shí)的計(jì)算結(jié)果如圖4所示.從圖4可知,聲源頻率為50 Hz時(shí),孤立子內(nèi)波對(duì)聲信號(hào)的傳播沒有特別大的影響,這是因?yàn)榈皖l聲波激發(fā)的簡(jiǎn)正波模態(tài)相對(duì)較少,孤立子內(nèi)波對(duì)各階簡(jiǎn)正波的耦合沒有產(chǎn)生太大的影響,聲波對(duì)孤立子內(nèi)波產(chǎn)生了饒射,所以內(nèi)波只對(duì)聲傳播產(chǎn)生了微小的影響,并沒有對(duì)聲場(chǎng)干涉結(jié)構(gòu)有很大的改變.隨著聲源頻率從50 Hz升到250 Hz,孤立子內(nèi)波對(duì)聲壓傳播損失的影響逐漸加大.當(dāng)聲源頻率f=250 Hz時(shí),聲壓傳播損失會(huì)因?yàn)楣铝⒆觾?nèi)波的存在急劇增加,對(duì)水平距離2 700 m后的傳播損失都產(chǎn)生了較大的影響.

    孤立子內(nèi)波會(huì)改變聲場(chǎng)的干涉結(jié)構(gòu),進(jìn)而引起聲壓傳播損失的起伏幅度和位置的變化,在不同距離處的聲壓傳播損失也會(huì)有起伏.仿真得到在不同聲源發(fā)射頻率下,存在孤立子內(nèi)波時(shí)的聲壓傳播損失與無(wú)內(nèi)波時(shí)的傳播損失差值如圖5所示:

    在聲壓傳播損失差值較大的聲源發(fā)射頻率,可以近似認(rèn)為孤立子內(nèi)波與聲波產(chǎn)生了共鳴現(xiàn)象.通過仿真分析可知,一般在聲源頻率低于70 Hz時(shí),孤立子內(nèi)波對(duì)聲傳播影響較小,隨著頻率的升高其影響逐漸加大并波動(dòng)起伏.

    3.2 聲源深度與接收深度不同時(shí)孤立子內(nèi)波對(duì)聲傳播的影響仿真

    為了更好研究孤立子內(nèi)波對(duì)聲傳播的影響,將聲源深度和接收深度分別設(shè)置為10 m、20 m、30 m,其他環(huán)境參數(shù)不變,發(fā)射信號(hào)頻率為250 Hz,仿真計(jì)算得到聲源深度10 m,分別計(jì)算不同聲源和接收深度下有孤立子內(nèi)波時(shí)聲壓傳播損失減去無(wú)內(nèi)波時(shí)聲壓傳播損失差值如表1、2、3所示.

    表1 聲源深度10 m時(shí)不同接收深度和水平距離處聲壓傳播損失差值 dB

    表2 聲源深度20 m時(shí)不同接收深度和水平距離處聲壓傳播損失差值 dB

    表3 聲源深度30 m時(shí)不同接收深度和水平距離處聲壓傳播損失差值 dB

    對(duì)仿真結(jié)果進(jìn)行分析可知:

    1)當(dāng)聲源位于海底深度10 m時(shí),有內(nèi)波影響的情況下,在接收深度10 m處的傳播損失增加了8 dB以上,在躍層處的傳播損失增加了2 dB,而躍層的下方,深度30 m處的傳播損失減少了11 dB.這種情況產(chǎn)生的原因是因?yàn)槁曉丛谲S層的上方,能量集中在這一部分,孤立子內(nèi)波使聲波產(chǎn)生散射和折射,造成其傳播損失增加.

    2)當(dāng)聲源位于海底深度30 m時(shí),躍層上方的傳播損失減少了6 dB以上,躍層處的傳播損失變化不是很明顯,躍層下方增加了1.7 dB左右.這種情況產(chǎn)生的原因是此時(shí)聲源能量集中在下層,孤立子內(nèi)波使躍層下方的傳播損失增加.

    3)當(dāng)聲源位于海底深度20 m時(shí),躍層上方的傳播損失增加了12 dB,下方減少了13 dB,躍層處的傳播損失增加了2 dB.這是因?yàn)楣铝⒆觾?nèi)波就位于躍層附近,此時(shí)孤立子內(nèi)波對(duì)傳播損失影響較小.

    由仿真分析可知,孤立子內(nèi)波使聲傳播產(chǎn)生散射,并改變聲場(chǎng)的能量分布.孤立子內(nèi)波會(huì)使聲場(chǎng)的干涉結(jié)構(gòu)更復(fù)雜.

    3.3 孤立子內(nèi)波對(duì)聲矢量場(chǎng)的影響仿真研究

    在低頻時(shí),孤立子內(nèi)波對(duì)聲場(chǎng)的作用較小,設(shè)聲源發(fā)射信號(hào)頻率為250 Hz,環(huán)境參數(shù)不變,聲源深度為30 m,接收深度為10 m,仿真計(jì)算得到孤立子內(nèi)波影響下聲壓、質(zhì)點(diǎn)振速和聲能流傳播損失對(duì)比如圖6所示.

    通過圖6(b)、(c)看出,孤立子內(nèi)波會(huì)使水平質(zhì)點(diǎn)振速傳播損失減小,使垂向質(zhì)點(diǎn)振速略微增加.從圖6(d)、(e)可看出孤立子內(nèi)波使水平、垂向聲能流密度傳播損失減少.

    4 結(jié)語(yǔ)

    基于COMSOL軟件進(jìn)行仿真,建立孤立子內(nèi)波模型,通過仿真研究分析了孤立子內(nèi)波對(duì)聲傳播產(chǎn)生的影響.研究結(jié)果表示,孤立子內(nèi)波對(duì)聲線具有折射作用,聲線發(fā)生了偏折,使聲線密集程度、聲波到達(dá)時(shí)間、聲場(chǎng)能量分布等都發(fā)生了改變.孤立子內(nèi)波對(duì)聲壓傳播損失和聲場(chǎng)干涉影響的研究發(fā)現(xiàn),孤立子內(nèi)波對(duì)聲場(chǎng)能量有明顯的改變.淺海的情況下,孤立子內(nèi)波對(duì)聲傳播的影響與聲波頻率有關(guān),低頻聲波的簡(jiǎn)正波模態(tài)較少,對(duì)孤立子內(nèi)波產(chǎn)生繞射作用,使孤立子內(nèi)波對(duì)聲場(chǎng)的影響很??;聲源位于躍層上、下方時(shí),孤立子內(nèi)波使躍層上、下方的聲壓傳播損失產(chǎn)生了明顯的變化.而在聲源位于躍層時(shí),傳播損失并沒有很大的變化.孤立子內(nèi)波對(duì)聲波有一定的散射作用,可使聲場(chǎng)能量分布更均勻. 對(duì)聲矢量場(chǎng)進(jìn)行的仿真研究發(fā)現(xiàn),孤立子內(nèi)波會(huì)使聲場(chǎng)能量小處水平振速的傳播損失減小,使垂直振速的傳播損失略微增加.在實(shí)際應(yīng)用中,可采用聲壓信號(hào)聯(lián)合處理等方法來減少孤立子內(nèi)波對(duì)聲傳播的影響.

    猜你喜歡
    內(nèi)波聲線聲壓
    孤立內(nèi)波對(duì)過渡海域聲場(chǎng)干涉結(jié)構(gòu)的影響分析
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    水聲中非直達(dá)聲下的聲速修正方法①
    內(nèi)波與死水,連潛艇都怕的海浪
    基于聲線法的特殊體育館模型中聲場(chǎng)均勻性分析
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    糾纏的曲線
    優(yōu)雅(2017年3期)2017-03-09 17:02:52
    三維溫度梯度場(chǎng)中本征聲線軌跡的求取*
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    SAR極化特征圖像與σ0圖像的海洋內(nèi)波可視性對(duì)比
    在线天堂最新版资源| 国产视频首页在线观看| 亚洲av电影在线观看一区二区三区| 91精品伊人久久大香线蕉| 黄色视频在线播放观看不卡| 午夜福利视频在线观看免费| 日本免费在线观看一区| 中文天堂在线官网| 一区在线观看完整版| 18禁裸乳无遮挡动漫免费视频| 97超视频在线观看视频| 国产精品国产三级国产av玫瑰| 国产男人的电影天堂91| 国产高清国产精品国产三级| 永久网站在线| 成年人免费黄色播放视频| 天天影视国产精品| 欧美精品高潮呻吟av久久| 日本欧美视频一区| 成人综合一区亚洲| 久久毛片免费看一区二区三区| 国产精品久久久久久久电影| 中文欧美无线码| √禁漫天堂资源中文www| 大片免费播放器 马上看| 最近的中文字幕免费完整| 天天影视国产精品| av网站免费在线观看视频| 99热网站在线观看| 免费观看的影片在线观看| 99久久人妻综合| 在线播放无遮挡| a级毛片在线看网站| 亚洲综合色网址| 中文字幕精品免费在线观看视频 | 王馨瑶露胸无遮挡在线观看| 91精品国产九色| 少妇的逼水好多| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 一级毛片电影观看| 啦啦啦中文免费视频观看日本| 色哟哟·www| 黄片播放在线免费| 国产成人a∨麻豆精品| 99久国产av精品国产电影| 亚洲精品第二区| av线在线观看网站| 免费av不卡在线播放| 视频区图区小说| 国产成人免费无遮挡视频| 嘟嘟电影网在线观看| 亚洲性久久影院| 亚洲精品久久久久久婷婷小说| √禁漫天堂资源中文www| 国产成人a∨麻豆精品| 欧美国产精品一级二级三级| 精品久久久久久久久亚洲| 欧美bdsm另类| 一级片'在线观看视频| 校园人妻丝袜中文字幕| 美女福利国产在线| 黑人猛操日本美女一级片| 91精品伊人久久大香线蕉| 午夜激情av网站| 国产乱来视频区| 51国产日韩欧美| 美女视频免费永久观看网站| 久久久久久久亚洲中文字幕| 亚洲第一av免费看| a级片在线免费高清观看视频| .国产精品久久| 欧美激情 高清一区二区三区| 91在线精品国自产拍蜜月| 少妇人妻 视频| 国产精品 国内视频| 国产午夜精品久久久久久一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲国产欧美在线一区| av一本久久久久| 18禁在线无遮挡免费观看视频| 亚洲欧美清纯卡通| 午夜视频国产福利| 欧美最新免费一区二区三区| av免费观看日本| 亚洲,一卡二卡三卡| 亚洲精品国产av蜜桃| 日韩熟女老妇一区二区性免费视频| 又大又黄又爽视频免费| 欧美 亚洲 国产 日韩一| 欧美一级a爱片免费观看看| 国产av码专区亚洲av| 免费不卡的大黄色大毛片视频在线观看| 夜夜骑夜夜射夜夜干| 亚洲国产毛片av蜜桃av| 2022亚洲国产成人精品| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| 亚洲精品乱码久久久久久按摩| 国产片内射在线| 一边摸一边做爽爽视频免费| 国产精品99久久久久久久久| 日日撸夜夜添| 99久久精品国产国产毛片| 亚洲三级黄色毛片| 国产永久视频网站| 国产黄色视频一区二区在线观看| 蜜桃在线观看..| 如日韩欧美国产精品一区二区三区 | 国产毛片在线视频| 91久久精品电影网| 啦啦啦啦在线视频资源| 国产精品麻豆人妻色哟哟久久| 97超碰精品成人国产| 亚洲欧美一区二区三区黑人 | 亚洲国产精品一区三区| 日韩欧美精品免费久久| 人妻少妇偷人精品九色| 18禁观看日本| 欧美激情国产日韩精品一区| 男女边摸边吃奶| 秋霞伦理黄片| 精品一品国产午夜福利视频| 成人亚洲精品一区在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品视频女| 五月开心婷婷网| av.在线天堂| 成人无遮挡网站| 亚洲精品av麻豆狂野| 丝袜脚勾引网站| 高清毛片免费看| 在线观看免费视频网站a站| 黑丝袜美女国产一区| 欧美日韩视频精品一区| 99久国产av精品国产电影| 日韩欧美一区视频在线观看| 哪个播放器可以免费观看大片| 午夜福利视频精品| 久久毛片免费看一区二区三区| 伊人亚洲综合成人网| 一本色道久久久久久精品综合| a级毛片免费高清观看在线播放| 天堂8中文在线网| 国模一区二区三区四区视频| 51国产日韩欧美| 青春草视频在线免费观看| 成年人免费黄色播放视频| 在线看a的网站| 少妇人妻久久综合中文| 日韩成人av中文字幕在线观看| av女优亚洲男人天堂| 新久久久久国产一级毛片| 国产精品一区二区三区四区免费观看| 人妻夜夜爽99麻豆av| www.av在线官网国产| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 91久久精品国产一区二区三区| 99热6这里只有精品| 亚洲av成人精品一区久久| 大陆偷拍与自拍| 日本vs欧美在线观看视频| 国产亚洲精品久久久com| 免费少妇av软件| 午夜福利网站1000一区二区三区| 高清av免费在线| 亚洲人与动物交配视频| 久久精品熟女亚洲av麻豆精品| 久久人妻熟女aⅴ| 国产成人91sexporn| 国产综合精华液| 99久久精品国产国产毛片| 自线自在国产av| 成年人免费黄色播放视频| 街头女战士在线观看网站| √禁漫天堂资源中文www| 成年人午夜在线观看视频| 日本91视频免费播放| 国产免费福利视频在线观看| 一区二区三区精品91| 秋霞伦理黄片| 女人久久www免费人成看片| 制服丝袜香蕉在线| 日韩中字成人| 黑人欧美特级aaaaaa片| 国产 精品1| 又黄又爽又刺激的免费视频.| 欧美老熟妇乱子伦牲交| 欧美xxxx性猛交bbbb| 中文字幕亚洲精品专区| 亚洲性久久影院| 蜜臀久久99精品久久宅男| 十分钟在线观看高清视频www| 亚洲色图综合在线观看| 国产精品一二三区在线看| av视频免费观看在线观看| 纵有疾风起免费观看全集完整版| 黑人高潮一二区| 狂野欧美激情性bbbbbb| 成年人免费黄色播放视频| 少妇丰满av| 久久青草综合色| 一边摸一边做爽爽视频免费| 午夜福利,免费看| 日韩av不卡免费在线播放| 多毛熟女@视频| 欧美日韩av久久| 国产日韩欧美在线精品| 色网站视频免费| 精品久久久久久电影网| 亚洲,一卡二卡三卡| 如何舔出高潮| 性高湖久久久久久久久免费观看| 国产男女内射视频| 极品少妇高潮喷水抽搐| 中文欧美无线码| 日本欧美国产在线视频| 亚洲国产日韩一区二区| 男女国产视频网站| 久久久久久久久久久丰满| 18在线观看网站| 18禁在线无遮挡免费观看视频| 热99久久久久精品小说推荐| xxx大片免费视频| 亚洲国产欧美日韩在线播放| 国产黄色免费在线视频| 成年人午夜在线观看视频| 欧美激情极品国产一区二区三区 | 亚洲国产精品一区三区| 久久精品熟女亚洲av麻豆精品| 最新中文字幕久久久久| 欧美xxxx性猛交bbbb| 人人妻人人澡人人爽人人夜夜| 国产精品欧美亚洲77777| 久久久久久人妻| 男女啪啪激烈高潮av片| 免费日韩欧美在线观看| 免费av不卡在线播放| 久久影院123| 又黄又爽又刺激的免费视频.| av免费观看日本| 国产亚洲精品久久久com| 人人妻人人爽人人添夜夜欢视频| 国产精品免费大片| 国产乱人偷精品视频| 国产极品天堂在线| 国产在视频线精品| 伊人久久精品亚洲午夜| 狂野欧美激情性xxxx在线观看| 国产一级毛片在线| 国产极品天堂在线| 国产成人aa在线观看| 国产av精品麻豆| 18禁在线播放成人免费| 日韩一区二区视频免费看| 国产熟女午夜一区二区三区 | 亚洲激情五月婷婷啪啪| 成年人午夜在线观看视频| 国产白丝娇喘喷水9色精品| 一级毛片电影观看| 日韩精品免费视频一区二区三区 | 国产亚洲午夜精品一区二区久久| 99热这里只有是精品在线观看| 天堂中文最新版在线下载| 日日摸夜夜添夜夜爱| 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| 久久免费观看电影| a级毛片黄视频| 国产免费现黄频在线看| 高清不卡的av网站| 精品少妇黑人巨大在线播放| 亚洲av电影在线观看一区二区三区| 国产成人精品一,二区| 成人国产麻豆网| 妹子高潮喷水视频| 这个男人来自地球电影免费观看 | 夫妻性生交免费视频一级片| 日本wwww免费看| 女性被躁到高潮视频| 热99久久久久精品小说推荐| 国产在线免费精品| 少妇高潮的动态图| 亚洲色图 男人天堂 中文字幕 | 国产成人午夜福利电影在线观看| 插逼视频在线观看| 大码成人一级视频| 春色校园在线视频观看| 午夜免费观看性视频| 91午夜精品亚洲一区二区三区| 久久久久久久久大av| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| 91久久精品国产一区二区成人| 亚洲欧美成人精品一区二区| 一级毛片电影观看| av视频免费观看在线观看| 国产一区亚洲一区在线观看| videos熟女内射| 最新中文字幕久久久久| 亚洲精品国产av蜜桃| 只有这里有精品99| 久久99热6这里只有精品| 永久免费av网站大全| 最后的刺客免费高清国语| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 考比视频在线观看| 久久久精品94久久精品| 日韩精品免费视频一区二区三区 | 制服人妻中文乱码| 热99久久久久精品小说推荐| av一本久久久久| 青春草视频在线免费观看| 一级毛片 在线播放| 亚洲精品国产av蜜桃| 制服诱惑二区| 超碰97精品在线观看| freevideosex欧美| 三级国产精品片| 天天躁夜夜躁狠狠久久av| 国产亚洲精品久久久com| 国产高清不卡午夜福利| av国产久精品久网站免费入址| 国产在线免费精品| 国产免费视频播放在线视频| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 不卡视频在线观看欧美| 青青草视频在线视频观看| 国产又色又爽无遮挡免| 成人影院久久| 久久亚洲国产成人精品v| 一级a做视频免费观看| 国产色婷婷99| 爱豆传媒免费全集在线观看| 一级爰片在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 丁香六月天网| 精品久久蜜臀av无| 成人综合一区亚洲| 国产男人的电影天堂91| 国产淫语在线视频| 大片电影免费在线观看免费| 男人添女人高潮全过程视频| 色哟哟·www| 亚洲丝袜综合中文字幕| 大码成人一级视频| 美女中出高潮动态图| 五月玫瑰六月丁香| 青春草视频在线免费观看| 国产精品人妻久久久影院| 精品亚洲成国产av| 桃花免费在线播放| 97精品久久久久久久久久精品| 久久久久久久久久人人人人人人| 亚洲情色 制服丝袜| 99久久综合免费| 免费观看性生交大片5| 久久久久久久国产电影| 国产成人精品福利久久| www.色视频.com| 国产精品一二三区在线看| 在线天堂最新版资源| 啦啦啦视频在线资源免费观看| 午夜视频国产福利| 中国三级夫妇交换| 亚洲av综合色区一区| 美女脱内裤让男人舔精品视频| 日韩强制内射视频| av黄色大香蕉| 91国产中文字幕| 在现免费观看毛片| 国产成人精品福利久久| 夫妻午夜视频| 伊人久久国产一区二区| 免费人妻精品一区二区三区视频| 一本大道久久a久久精品| 免费看光身美女| 国产国拍精品亚洲av在线观看| 欧美日韩亚洲高清精品| 久久免费观看电影| 一区二区三区乱码不卡18| 一区二区日韩欧美中文字幕 | 极品少妇高潮喷水抽搐| 国产永久视频网站| 午夜福利视频精品| 亚洲欧美成人精品一区二区| 爱豆传媒免费全集在线观看| 中文字幕精品免费在线观看视频 | 国产成人午夜福利电影在线观看| 午夜福利网站1000一区二区三区| 国产精品人妻久久久影院| 婷婷色av中文字幕| 成年美女黄网站色视频大全免费 | 伊人久久精品亚洲午夜| 免费久久久久久久精品成人欧美视频 | 午夜激情久久久久久久| 高清毛片免费看| 久久99精品国语久久久| 一本色道久久久久久精品综合| 国产乱来视频区| 狠狠婷婷综合久久久久久88av| 日韩成人av中文字幕在线观看| 亚洲人成网站在线观看播放| 亚洲精品日韩av片在线观看| 国产精品嫩草影院av在线观看| 91精品三级在线观看| 少妇的逼好多水| 亚洲成人一二三区av| av网站免费在线观看视频| 能在线免费看毛片的网站| 久久精品国产亚洲av涩爱| √禁漫天堂资源中文www| 伊人久久国产一区二区| 欧美日韩精品成人综合77777| 另类精品久久| 这个男人来自地球电影免费观看 | av在线app专区| 国产成人freesex在线| 高清黄色对白视频在线免费看| 观看美女的网站| 日日啪夜夜爽| 亚洲精品,欧美精品| 大香蕉97超碰在线| 免费大片黄手机在线观看| 在线观看三级黄色| 女的被弄到高潮叫床怎么办| 国产黄色免费在线视频| 精品一区二区三区视频在线| 精品国产一区二区三区久久久樱花| xxxhd国产人妻xxx| 在线观看免费视频网站a站| 国产成人免费观看mmmm| 91精品国产九色| av在线观看视频网站免费| 欧美日韩av久久| 久久久久久久久久人人人人人人| 视频中文字幕在线观看| 赤兔流量卡办理| 久久人人爽av亚洲精品天堂| 久久人人爽人人爽人人片va| 在线观看一区二区三区激情| 国产一区亚洲一区在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲一级一片aⅴ在线观看| 18禁观看日本| av免费观看日本| 亚洲国产欧美日韩在线播放| 欧美日韩成人在线一区二区| 99热全是精品| 国产黄色免费在线视频| 午夜激情福利司机影院| 国产精品国产三级专区第一集| 久久精品国产亚洲网站| freevideosex欧美| 青春草国产在线视频| 欧美日韩成人在线一区二区| 丝袜美足系列| 人妻制服诱惑在线中文字幕| 校园人妻丝袜中文字幕| 午夜精品国产一区二区电影| 精品国产乱码久久久久久小说| 国产精品蜜桃在线观看| 五月伊人婷婷丁香| 欧美亚洲日本最大视频资源| 大香蕉久久成人网| a级毛色黄片| 久久99精品国语久久久| 啦啦啦啦在线视频资源| 亚洲精品久久久久久婷婷小说| 国产精品一区www在线观看| 成年美女黄网站色视频大全免费 | 丁香六月天网| 18+在线观看网站| 日本免费在线观看一区| 国产黄片视频在线免费观看| 在线观看一区二区三区激情| 超色免费av| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 国产永久视频网站| 精品少妇久久久久久888优播| 国产亚洲精品第一综合不卡 | 免费黄网站久久成人精品| 男女啪啪激烈高潮av片| 久久精品夜色国产| 一级,二级,三级黄色视频| 简卡轻食公司| 久久免费观看电影| 免费黄色在线免费观看| freevideosex欧美| 少妇熟女欧美另类| 婷婷色麻豆天堂久久| 久久人人爽人人片av| 亚洲精品视频女| 午夜激情福利司机影院| 欧美亚洲 丝袜 人妻 在线| 国产老妇伦熟女老妇高清| 日本黄大片高清| 久久久久网色| a级毛片在线看网站| 久久久久国产精品人妻一区二区| 亚洲av国产av综合av卡| 在线观看www视频免费| 黄片无遮挡物在线观看| 嘟嘟电影网在线观看| 欧美人与性动交α欧美精品济南到 | 美女大奶头黄色视频| 亚洲成人一二三区av| 亚洲一级一片aⅴ在线观看| a级片在线免费高清观看视频| 毛片一级片免费看久久久久| 伊人亚洲综合成人网| 在线亚洲精品国产二区图片欧美 | 99re6热这里在线精品视频| 亚洲综合精品二区| 肉色欧美久久久久久久蜜桃| 国内精品宾馆在线| 大话2 男鬼变身卡| 国产片特级美女逼逼视频| 亚洲av欧美aⅴ国产| 丝袜脚勾引网站| 一级,二级,三级黄色视频| 久久久精品区二区三区| 欧美日韩精品成人综合77777| 国产精品国产三级国产专区5o| 日产精品乱码卡一卡2卡三| 91精品国产九色| 最黄视频免费看| 久久精品国产鲁丝片午夜精品| 女人精品久久久久毛片| 国产亚洲精品第一综合不卡 | 男人操女人黄网站| 久热久热在线精品观看| kizo精华| .国产精品久久| 欧美日韩国产mv在线观看视频| 日本黄色日本黄色录像| 免费观看性生交大片5| 黄色配什么色好看| 国精品久久久久久国模美| 国产白丝娇喘喷水9色精品| 能在线免费看毛片的网站| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 国产成人精品无人区| 欧美日韩av久久| 国产极品天堂在线| 超碰97精品在线观看| 国产 一区精品| 亚洲人成网站在线播| 能在线免费看毛片的网站| 最近中文字幕2019免费版| 2021少妇久久久久久久久久久| 美女大奶头黄色视频| 久久精品国产亚洲网站| 国产 一区精品| av黄色大香蕉| 一边亲一边摸免费视频| 欧美老熟妇乱子伦牲交| 国产av国产精品国产| 亚洲人成77777在线视频| 国产精品久久久久成人av| 亚洲精品aⅴ在线观看| 国产乱来视频区| 亚洲色图 男人天堂 中文字幕 | 国产黄片视频在线免费观看| 国产精品一区二区在线观看99| 一区二区三区免费毛片| 一级毛片黄色毛片免费观看视频| 高清欧美精品videossex| 我的老师免费观看完整版| 狂野欧美激情性bbbbbb| 国产精品无大码| 丝瓜视频免费看黄片| av专区在线播放| 色哟哟·www| 蜜桃在线观看..| 国产成人av激情在线播放 | 性高湖久久久久久久久免费观看| videossex国产| 中文天堂在线官网| 一边摸一边做爽爽视频免费| 亚洲精品日韩在线中文字幕| 狂野欧美激情性bbbbbb| 91国产中文字幕| 亚洲美女搞黄在线观看| 久久久午夜欧美精品| 黑人猛操日本美女一级片| 久久久欧美国产精品| 国产日韩一区二区三区精品不卡 | 日韩在线高清观看一区二区三区| 亚洲色图综合在线观看| 国产有黄有色有爽视频| 人妻少妇偷人精品九色| 亚洲精品,欧美精品| 欧美日韩亚洲高清精品| 777米奇影视久久| 亚洲国产精品999| 一本久久精品| 视频在线观看一区二区三区| 最近的中文字幕免费完整| 亚洲天堂av无毛| 99久久人妻综合| 午夜久久久在线观看| 亚洲,欧美,日韩| 欧美日韩国产mv在线观看视频| 伊人亚洲综合成人网| 亚洲av综合色区一区| 美女国产视频在线观看| 日本欧美视频一区| 亚洲av电影在线观看一区二区三区|