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

    海塘越浪過程SPH模擬

    2014-11-14 11:41楊洲顏君來曾甄
    珠江水運(yùn) 2014年18期
    關(guān)鍵詞:越浪數(shù)值模擬

    楊洲+顏君來+曾甄

    摘 要:隨著計(jì)算機(jī)軟硬件水平的不斷發(fā)展,數(shù)值模擬技術(shù)已經(jīng)越來越廣泛地應(yīng)用到海岸工程領(lǐng)域中?,F(xiàn)今波浪爬高越浪的數(shù)值模擬技術(shù)主要有網(wǎng)格方法和無網(wǎng)格方法兩種。和網(wǎng)格方法相比,無網(wǎng)格方法在計(jì)算海塘越浪的連續(xù)性和大變形問題上體現(xiàn)出了優(yōu)勢。

    關(guān)鍵詞:越浪 SPH方法 數(shù)值模擬

    浙江省地理位置特殊,海岸線綿長,臺風(fēng)災(zāi)害頻繁,一旦海塘受損,將會造成不可估量的損失。根據(jù)以往潰壩的經(jīng)驗(yàn),越浪是造成海塘潰損的主要原因之一。隨著計(jì)算機(jī)軟硬件水平的不斷發(fā)展,海塘越浪的研究當(dāng)中不斷地引入數(shù)值模擬技術(shù)。數(shù)值模擬海塘的越浪以及計(jì)算海塘的越浪量,通過更加符合實(shí)際的邊界條件等計(jì)算的越浪量,比通過規(guī)范公式計(jì)算的越浪量更加符合工程安全要求,同時模擬的越浪過程可以直觀形象地了解海塘越浪的情況;數(shù)值模擬不需要專門的試驗(yàn)室、儀器和相關(guān)的工作人員,計(jì)算周期要短得多,方便進(jìn)行多次計(jì)算,相比傳統(tǒng)的物理模型試驗(yàn),數(shù)值模擬節(jié)約了大量的人力物力,同時在試驗(yàn)結(jié)果的精確度上又有一定的保證。

    如今海塘越浪的數(shù)值模擬計(jì)算方法主要有兩種:網(wǎng)格方法和無網(wǎng)格方法。 1977年Lucy、Gingold、Monaghan]等分別提出了光滑質(zhì)點(diǎn)流體動力學(xué)(SPH)方法,近幾年來,光滑粒子流體動力學(xué)(SPH)方法因?yàn)樗睦窭嗜仗匦?,在模擬自由表面流體的應(yīng)用正變得熱門,但目前而言在海塘(堤)越浪的定量問題應(yīng)用研究還并不廣泛。

    本文基于DualSPHysics開源平臺,建立二維數(shù)值水槽,模擬規(guī)則波在海塘上的爬坡和越浪,并與物理模型試驗(yàn)結(jié)果進(jìn)行了對比,驗(yàn)證了數(shù)值模型的可靠性。再通過數(shù)值模型計(jì)算越浪量,通過和規(guī)范公式計(jì)算得到的以及物理模型試驗(yàn)得到的越浪量進(jìn)行對比,實(shí)現(xiàn)SPH方法在工程實(shí)踐中的應(yīng)用。

    1.SPH方法介紹

    SPH方法有兩個主要的步驟:第一步是積分表達(dá)區(qū)域函數(shù)的函數(shù)近似;第二個步驟就是粒子近似。通過函數(shù)近似,可以將任意場函數(shù)(質(zhì)量場、溫度場等)進(jìn)行積分表達(dá)形式,并進(jìn)行粒子近似,將在粒子i處的場函數(shù)積分表達(dá)式進(jìn)一步離散為迭加形式。由連續(xù)性方程、動量方程和能量方程組成流體的基本控制方程即N-S方程,并分別通過流體的質(zhì)量守恒、動量守恒和能量守恒原理推導(dǎo)得出。

    2.物理模型驗(yàn)證

    選取某工程一個斷面進(jìn)行物理模型試驗(yàn)和SPH數(shù)值模擬試驗(yàn),然后對比兩個試驗(yàn)?zāi)P椭胁ɡ说倪\(yùn)動過程。當(dāng)試驗(yàn)開始波浪沿著海塘護(hù)面向上傳播,在經(jīng)過海塘平臺處時,由于水深變淺,波浪開始變形翻卷并破碎。水體沖擊防浪墻,并沿防浪墻向上爬,由于防浪墻擋浪面的弧線形式,水體沖上一定高度開始翻卷并下落。接著擋浪墻前面的水體退去,高度降低,又開始重復(fù)上面的過程。通過對比物理模型和數(shù)值模型波浪的運(yùn)動過程,可以發(fā)現(xiàn)兩個過程基本一致。

    在防浪墻前面選取測量點(diǎn)(如圖1所示)測量水體高度,由于SPH模擬的越浪過程剛開始由于水體內(nèi)部還沒有產(chǎn)生持續(xù)且穩(wěn)定的紊流使得前面幾個波浪一般偏大,到第五個開始水體周期內(nèi)的運(yùn)動過程趨于一致,故這次測量從第5個波浪開始。測量18s的結(jié)果如圖2所示。

    通過觀察上圖,在所設(shè)測量點(diǎn)的水體高度變化中SPH模擬的值和物模試驗(yàn)中的值非常吻合,最大誤差范圍在16%以內(nèi),主要是由于SPH方法在模擬海塘越浪過程DualSPHysics對固液邊界處理不夠完善以及物理模型試驗(yàn)中存在一個人為誤差。驗(yàn)證了SPH方法在海塘越浪過程模擬中的可靠性。

    3.工程應(yīng)用

    3.1越浪過程模擬

    通過dualSPHysics開源平臺計(jì)算的海塘越浪過程大致可以分為3個階段,分別為:爬坡階段、飛濺階段、越浪流過程。首先當(dāng)水流上爬至海塘中段的平臺,水流及其速度都比較穩(wěn)定,再往上爬,波浪翻卷變形,水舌產(chǎn)生較大的速度。水流遇到防浪墻并沖擊防浪墻,水流改變方向沿著防浪墻往上爬,水流在防浪墻的頂端處產(chǎn)生局部的速度峰值,但防浪墻未能阻擋全部水流;接著進(jìn)入飛濺階段,躍入空中的水流質(zhì)點(diǎn)速度方向均朝著斜上方,并有水滴飛濺,之后雖然水流還在不斷爬坡,但是水流的速度開始變小,這是仍在不斷爬坡,但水流仍保持一個較小的指向斜上方的塘頂速度,在到達(dá)海塘上空的水流下降的同時,前坡上的水流匯合,這時在前堤上的大部分速度接近于零。之后前坡上的水流開始回落,此時坡前平臺附近的水流速度剛剛降落至零點(diǎn),回落水流與之相沖,形成一個明顯的回落波形。最后是越浪流過程,越浪流過程和前面的飛濺過程時間上有部分重合,當(dāng)波高較大時,有些水質(zhì)點(diǎn)會直接落在后坡上,而下落的水流砸落在防浪墻頂上后逐漸形成越浪流這時水流將分成兩部分,靠近防浪墻的水流將沖回原水流當(dāng)中,并卷入氣泡,靠近后坡的水流年則順著塘頂向后坡流去,飛濺在空中的水花也會有部分回落這部分水流當(dāng)中。

    在消浪平臺的左端、折角處、右端以及防浪墻前端分別設(shè)置水體高度測量點(diǎn),記為1#、2#、3#和4#(如圖3所示)。測量點(diǎn)水體高度測量值見表1,測量點(diǎn)水體高度變化圖見圖4。

    1#測量點(diǎn)初段水體高度變化曲線較平時由于前一個波浪的回落力度和下一個波浪的雍高力度相當(dāng),是此處的水體高度能夠保持一段時間穩(wěn)定。接著隨著波浪的進(jìn)一步傳播由于1#測量點(diǎn)為一個消浪平臺,波浪的傳播形態(tài)能夠基本保持一個完整的波形,同時由于消浪平臺的高程太高,會使水體有雍高現(xiàn)象,故水體高度的變化曲線的峰值比谷值大;2#測量點(diǎn)之后就是有坡度的消浪坡面,高程逐漸抬高,波浪到2#測量點(diǎn)會受到阻礙,水體高程較平穩(wěn)。當(dāng)水體回落,由于前后高程差較大,產(chǎn)生較大的谷值;3#測量點(diǎn)由于靠近擋浪墻,水體遇到擋浪墻會沿著防浪墻往上爬,高度升高比較快,同時也是因?yàn)閾趵藟Φ淖饔?,雍高了水位,是水體高度下降較緩慢;4#測量點(diǎn)由于在擋浪墻前,水體高度變化圖跟3#測量點(diǎn)相差不大,但是水體高度的變化不管是升高還是降低都平緩一點(diǎn)。

    3.2越浪量計(jì)算

    和物理模型一樣,越浪過程中的前幾個波浪往往偏大,本文考察第5個和第6個波浪作用過程,以此計(jì)算越浪量。無風(fēng)條件下三種方式得到的越浪量比較見表2。

    通過對比發(fā)現(xiàn)通過規(guī)范公式計(jì)算的越浪量值比物理模型和SPH模型得到的值都大,這是因?yàn)槲锢砟P秃蚐PH模擬試驗(yàn)采用的是更加符合實(shí)際情況的邊界條件,而計(jì)算公式中沒有考慮反弧挑浪式防浪墻的挑浪作用和消浪平臺的消浪作用,因此認(rèn)為物理模型試驗(yàn)和SPH模擬試驗(yàn)的結(jié)果更符合真實(shí)情況。

    通過SPH方法得到的越浪量幾乎都略大于物理試驗(yàn)值,比如本文SPH方法所模擬的結(jié)果值比物理模型試驗(yàn)值要大大約29%,造成這樣誤差的原因是建立數(shù)值模型的時候簡化了消浪平臺,使其消浪作用減弱,另外本文采用設(shè)置一個水槽來接水最后按照所接水體體積來計(jì)算越浪量,這個方法存在人為誤差,另一個計(jì)算誤差是Dual SPHysics對于固液邊界處理不夠完善,首先,液體粒子和固壁粒子之間并未添加摩阻機(jī)制,使得波浪在斜坡堤的爬坡過程中波高偏大。其次,在使用連續(xù)密度法求解粒子密度導(dǎo)數(shù)時,將粒子相對速度引入了算式。由于固壁粒子和液體粒子的相對速度要明顯大于液體粒子之間的相對速度,直接導(dǎo)致邊界粒子的密度和壓強(qiáng)大于正常值,液體粒子收排斥而在固壁粒子之間逐漸形成無粒子層。該無粒子層的厚度約在2個光滑長度左右,起到一種類似“氣墊”的效果,將越浪水流托起,使其能更輕松地越過堤頂,如圖4所示。

    4.結(jié)論

    本文基于DualSPHysics開源平臺,采用光滑質(zhì)點(diǎn)流體動力學(xué)(SPH)進(jìn)行嘗試,通過選取某工程一個斷面的物理模型試驗(yàn)的結(jié)果和SPH方式模擬的海塘擋浪過程以及預(yù)設(shè)在防浪墻前的水體高度測量點(diǎn)位置的水體高度變化過程的對比來驗(yàn)證SPH方法在水動力模型試驗(yàn)中的可靠性,最后把SPH方法應(yīng)用在漩門三期工程中,模擬海塘的越浪過程并驗(yàn)證越浪量的準(zhǔn)確性。

    本文SPH方法在海塘越浪中作了初步的研究,為保證該方法的穩(wěn)定性和可靠性,仍需進(jìn)行大量的測試和改進(jìn)工作。另外,SPH方法的基礎(chǔ)理論尚存很多待解決的問題,特別是在固壁粒子與液體粒子相互作用方面,亟需引入一種更為有效的處理方法,以便實(shí)現(xiàn)準(zhǔn)確而快速的邊界壓強(qiáng)計(jì)算。

    參考文獻(xiàn):

    [1] Lucy L B. A Approach to the Testing of the Fission Hypothesis. The Astron. J. , 1977, 8 (12): 1013- 1024p.

    [2]Gingold RA, Monaghan JJ. Smoothed Partide Hydrodynamics: Theory and Applications to Non-spherical Stars. Mon. Not. Roy. Astrou. Soc. , 1977, 18: 375 -389P.

    [3] Monaghan J J. Smoothed particle hydrodynamics[j]. Reports On Progress In Physics. 2005, 68(8): 1703-1759.

    通過對比發(fā)現(xiàn)通過規(guī)范公式計(jì)算的越浪量值比物理模型和SPH模型得到的值都大,這是因?yàn)槲锢砟P秃蚐PH模擬試驗(yàn)采用的是更加符合實(shí)際情況的邊界條件,而計(jì)算公式中沒有考慮反弧挑浪式防浪墻的挑浪作用和消浪平臺的消浪作用,因此認(rèn)為物理模型試驗(yàn)和SPH模擬試驗(yàn)的結(jié)果更符合真實(shí)情況。

    通過SPH方法得到的越浪量幾乎都略大于物理試驗(yàn)值,比如本文SPH方法所模擬的結(jié)果值比物理模型試驗(yàn)值要大大約29%,造成這樣誤差的原因是建立數(shù)值模型的時候簡化了消浪平臺,使其消浪作用減弱,另外本文采用設(shè)置一個水槽來接水最后按照所接水體體積來計(jì)算越浪量,這個方法存在人為誤差,另一個計(jì)算誤差是Dual SPHysics對于固液邊界處理不夠完善,首先,液體粒子和固壁粒子之間并未添加摩阻機(jī)制,使得波浪在斜坡堤的爬坡過程中波高偏大。其次,在使用連續(xù)密度法求解粒子密度導(dǎo)數(shù)時,將粒子相對速度引入了算式。由于固壁粒子和液體粒子的相對速度要明顯大于液體粒子之間的相對速度,直接導(dǎo)致邊界粒子的密度和壓強(qiáng)大于正常值,液體粒子收排斥而在固壁粒子之間逐漸形成無粒子層。該無粒子層的厚度約在2個光滑長度左右,起到一種類似“氣墊”的效果,將越浪水流托起,使其能更輕松地越過堤頂,如圖4所示。

    4.結(jié)論

    本文基于DualSPHysics開源平臺,采用光滑質(zhì)點(diǎn)流體動力學(xué)(SPH)進(jìn)行嘗試,通過選取某工程一個斷面的物理模型試驗(yàn)的結(jié)果和SPH方式模擬的海塘擋浪過程以及預(yù)設(shè)在防浪墻前的水體高度測量點(diǎn)位置的水體高度變化過程的對比來驗(yàn)證SPH方法在水動力模型試驗(yàn)中的可靠性,最后把SPH方法應(yīng)用在漩門三期工程中,模擬海塘的越浪過程并驗(yàn)證越浪量的準(zhǔn)確性。

    本文SPH方法在海塘越浪中作了初步的研究,為保證該方法的穩(wěn)定性和可靠性,仍需進(jìn)行大量的測試和改進(jìn)工作。另外,SPH方法的基礎(chǔ)理論尚存很多待解決的問題,特別是在固壁粒子與液體粒子相互作用方面,亟需引入一種更為有效的處理方法,以便實(shí)現(xiàn)準(zhǔn)確而快速的邊界壓強(qiáng)計(jì)算。

    參考文獻(xiàn):

    [1] Lucy L B. A Approach to the Testing of the Fission Hypothesis. The Astron. J. , 1977, 8 (12): 1013- 1024p.

    [2]Gingold RA, Monaghan JJ. Smoothed Partide Hydrodynamics: Theory and Applications to Non-spherical Stars. Mon. Not. Roy. Astrou. Soc. , 1977, 18: 375 -389P.

    [3] Monaghan J J. Smoothed particle hydrodynamics[j]. Reports On Progress In Physics. 2005, 68(8): 1703-1759.

    通過對比發(fā)現(xiàn)通過規(guī)范公式計(jì)算的越浪量值比物理模型和SPH模型得到的值都大,這是因?yàn)槲锢砟P秃蚐PH模擬試驗(yàn)采用的是更加符合實(shí)際情況的邊界條件,而計(jì)算公式中沒有考慮反弧挑浪式防浪墻的挑浪作用和消浪平臺的消浪作用,因此認(rèn)為物理模型試驗(yàn)和SPH模擬試驗(yàn)的結(jié)果更符合真實(shí)情況。

    通過SPH方法得到的越浪量幾乎都略大于物理試驗(yàn)值,比如本文SPH方法所模擬的結(jié)果值比物理模型試驗(yàn)值要大大約29%,造成這樣誤差的原因是建立數(shù)值模型的時候簡化了消浪平臺,使其消浪作用減弱,另外本文采用設(shè)置一個水槽來接水最后按照所接水體體積來計(jì)算越浪量,這個方法存在人為誤差,另一個計(jì)算誤差是Dual SPHysics對于固液邊界處理不夠完善,首先,液體粒子和固壁粒子之間并未添加摩阻機(jī)制,使得波浪在斜坡堤的爬坡過程中波高偏大。其次,在使用連續(xù)密度法求解粒子密度導(dǎo)數(shù)時,將粒子相對速度引入了算式。由于固壁粒子和液體粒子的相對速度要明顯大于液體粒子之間的相對速度,直接導(dǎo)致邊界粒子的密度和壓強(qiáng)大于正常值,液體粒子收排斥而在固壁粒子之間逐漸形成無粒子層。該無粒子層的厚度約在2個光滑長度左右,起到一種類似“氣墊”的效果,將越浪水流托起,使其能更輕松地越過堤頂,如圖4所示。

    4.結(jié)論

    本文基于DualSPHysics開源平臺,采用光滑質(zhì)點(diǎn)流體動力學(xué)(SPH)進(jìn)行嘗試,通過選取某工程一個斷面的物理模型試驗(yàn)的結(jié)果和SPH方式模擬的海塘擋浪過程以及預(yù)設(shè)在防浪墻前的水體高度測量點(diǎn)位置的水體高度變化過程的對比來驗(yàn)證SPH方法在水動力模型試驗(yàn)中的可靠性,最后把SPH方法應(yīng)用在漩門三期工程中,模擬海塘的越浪過程并驗(yàn)證越浪量的準(zhǔn)確性。

    本文SPH方法在海塘越浪中作了初步的研究,為保證該方法的穩(wěn)定性和可靠性,仍需進(jìn)行大量的測試和改進(jìn)工作。另外,SPH方法的基礎(chǔ)理論尚存很多待解決的問題,特別是在固壁粒子與液體粒子相互作用方面,亟需引入一種更為有效的處理方法,以便實(shí)現(xiàn)準(zhǔn)確而快速的邊界壓強(qiáng)計(jì)算。

    參考文獻(xiàn):

    [1] Lucy L B. A Approach to the Testing of the Fission Hypothesis. The Astron. J. , 1977, 8 (12): 1013- 1024p.

    [2]Gingold RA, Monaghan JJ. Smoothed Partide Hydrodynamics: Theory and Applications to Non-spherical Stars. Mon. Not. Roy. Astrou. Soc. , 1977, 18: 375 -389P.

    [3] Monaghan J J. Smoothed particle hydrodynamics[j]. Reports On Progress In Physics. 2005, 68(8): 1703-1759.

    猜你喜歡
    越浪數(shù)值模擬
    復(fù)式海堤越浪水體數(shù)值模擬研究
    臺風(fēng)過程下復(fù)式海堤越浪量計(jì)算方法研究
    中外港口規(guī)范對比研究Ⅱ:海堤越浪量標(biāo)準(zhǔn)
    拋石斜坡堤后坡面穩(wěn)定重量的模型試驗(yàn)研究
    海堤背水坡加筋草皮抗沖蝕能力試驗(yàn)研究
    允許部分越浪海堤設(shè)計(jì)在工程中的應(yīng)用
    夜夜爽夜夜爽视频| 亚洲精品一二三| 狂野欧美白嫩少妇大欣赏| 国产高清国产精品国产三级 | 一级a做视频免费观看| 美女主播在线视频| 五月伊人婷婷丁香| 最新中文字幕久久久久| 哪个播放器可以免费观看大片| 亚洲精品日韩在线中文字幕| 成年美女黄网站色视频大全免费 | av一本久久久久| 亚洲一级一片aⅴ在线观看| 黑人猛操日本美女一级片| 一本色道久久久久久精品综合| 日本色播在线视频| 22中文网久久字幕| 我的老师免费观看完整版| av天堂中文字幕网| 性高湖久久久久久久久免费观看| 涩涩av久久男人的天堂| 国产一区有黄有色的免费视频| 51国产日韩欧美| 777米奇影视久久| 十分钟在线观看高清视频www | 中文字幕久久专区| 晚上一个人看的免费电影| 国产爽快片一区二区三区| 国产乱来视频区| 人妻系列 视频| 久久人人爽人人片av| 麻豆成人午夜福利视频| 国产成人a∨麻豆精品| 国产成人免费无遮挡视频| 国产欧美另类精品又又久久亚洲欧美| 一区在线观看完整版| 国产亚洲5aaaaa淫片| 国产免费福利视频在线观看| 夫妻性生交免费视频一级片| 少妇高潮的动态图| 国产精品国产av在线观看| 在线观看一区二区三区激情| 狠狠精品人妻久久久久久综合| 我要看黄色一级片免费的| 亚洲高清免费不卡视频| 久久97久久精品| 男女国产视频网站| 国产成人精品婷婷| 99热网站在线观看| 亚洲精品乱码久久久久久按摩| 国产精品一区二区性色av| 国产精品一区二区性色av| 天天躁夜夜躁狠狠久久av| 菩萨蛮人人尽说江南好唐韦庄| 91在线精品国自产拍蜜月| 国产中年淑女户外野战色| 日韩视频在线欧美| 黑丝袜美女国产一区| 日本色播在线视频| 18禁在线无遮挡免费观看视频| 亚洲国产成人一精品久久久| 国产白丝娇喘喷水9色精品| 在线观看免费日韩欧美大片 | 亚洲av二区三区四区| 观看av在线不卡| 一级毛片久久久久久久久女| 亚洲精品456在线播放app| 人人妻人人爽人人添夜夜欢视频 | 啦啦啦视频在线资源免费观看| 久久久久国产精品人妻一区二区| 国产午夜精品一二区理论片| 22中文网久久字幕| 国产一区亚洲一区在线观看| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 国产亚洲5aaaaa淫片| 久久ye,这里只有精品| 水蜜桃什么品种好| 国产精品国产av在线观看| 永久免费av网站大全| 国产精品av视频在线免费观看| 亚洲av电影在线观看一区二区三区| 深爱激情五月婷婷| 亚洲综合色惰| 热99国产精品久久久久久7| 大香蕉97超碰在线| 亚洲精品456在线播放app| a级一级毛片免费在线观看| 人体艺术视频欧美日本| 99久久精品热视频| 99久久综合免费| 在线免费观看不下载黄p国产| 免费少妇av软件| 伊人久久精品亚洲午夜| 精华霜和精华液先用哪个| 有码 亚洲区| 99久久精品国产国产毛片| 成年免费大片在线观看| 伊人久久精品亚洲午夜| 观看av在线不卡| 一二三四中文在线观看免费高清| 99热这里只有是精品在线观看| 80岁老熟妇乱子伦牲交| 制服丝袜香蕉在线| 国产精品国产三级国产av玫瑰| 在线观看免费高清a一片| 国产精品一区www在线观看| 一级爰片在线观看| 亚洲国产日韩一区二区| 亚洲美女视频黄频| 久久韩国三级中文字幕| 亚洲国产精品国产精品| 欧美+日韩+精品| 亚洲欧美精品自产自拍| 日韩大片免费观看网站| 欧美成人一区二区免费高清观看| 又爽又黄a免费视频| 黄色日韩在线| av国产免费在线观看| 啦啦啦中文免费视频观看日本| 欧美精品亚洲一区二区| 国产av一区二区精品久久 | 国产成人91sexporn| 99九九线精品视频在线观看视频| 一个人看的www免费观看视频| 国产精品福利在线免费观看| 欧美一区二区亚洲| 又爽又黄a免费视频| 国产精品不卡视频一区二区| 黄片wwwwww| 一级毛片久久久久久久久女| 妹子高潮喷水视频| 黄色日韩在线| 爱豆传媒免费全集在线观看| 少妇熟女欧美另类| 精品少妇黑人巨大在线播放| 观看av在线不卡| 亚洲国产欧美人成| 少妇人妻精品综合一区二区| 在线精品无人区一区二区三 | 国产中年淑女户外野战色| 国产深夜福利视频在线观看| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| 如何舔出高潮| 九草在线视频观看| 丰满人妻一区二区三区视频av| 卡戴珊不雅视频在线播放| 久久精品国产a三级三级三级| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 狂野欧美激情性xxxx在线观看| 午夜老司机福利剧场| av线在线观看网站| 成人亚洲欧美一区二区av| 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂 | 制服丝袜香蕉在线| 国产亚洲91精品色在线| 人妻少妇偷人精品九色| 秋霞伦理黄片| 亚洲欧美一区二区三区黑人 | 亚洲精品视频女| 丰满乱子伦码专区| 国产精品熟女久久久久浪| 老师上课跳d突然被开到最大视频| 波野结衣二区三区在线| 日本vs欧美在线观看视频 | 国产v大片淫在线免费观看| 久热这里只有精品99| 国产乱人视频| 校园人妻丝袜中文字幕| 国产精品蜜桃在线观看| 精品熟女少妇av免费看| 18禁裸乳无遮挡动漫免费视频| 国产片特级美女逼逼视频| 一个人免费看片子| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 国产91av在线免费观看| 国产男女超爽视频在线观看| 精品久久久久久电影网| 欧美精品亚洲一区二区| 久久亚洲国产成人精品v| 黄片wwwwww| 久久韩国三级中文字幕| 亚洲国产毛片av蜜桃av| 亚洲国产高清在线一区二区三| 免费久久久久久久精品成人欧美视频 | 中文欧美无线码| 国产精品爽爽va在线观看网站| 国产精品一区www在线观看| 中文字幕久久专区| 亚洲av中文av极速乱| 日韩av免费高清视频| 中文乱码字字幕精品一区二区三区| 成人综合一区亚洲| 亚洲精品日韩在线中文字幕| 丰满少妇做爰视频| 国产成人a区在线观看| 亚洲av日韩在线播放| 少妇的逼好多水| 亚洲欧美一区二区三区国产| 人体艺术视频欧美日本| 少妇的逼好多水| 国产淫片久久久久久久久| 成人一区二区视频在线观看| 欧美日韩在线观看h| 久久午夜福利片| 老师上课跳d突然被开到最大视频| 午夜福利在线观看免费完整高清在| 亚洲综合色惰| 麻豆成人午夜福利视频| 久久婷婷青草| 国内揄拍国产精品人妻在线| 成人无遮挡网站| 中文字幕制服av| 在线精品无人区一区二区三 | 中文精品一卡2卡3卡4更新| 好男人视频免费观看在线| 99热这里只有精品一区| 九九在线视频观看精品| 精品久久久久久久久av| 成年人午夜在线观看视频| 老师上课跳d突然被开到最大视频| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的 | 国产伦精品一区二区三区四那| 中文天堂在线官网| 天美传媒精品一区二区| 久久这里有精品视频免费| 22中文网久久字幕| 日本av手机在线免费观看| 香蕉精品网在线| 岛国毛片在线播放| 日韩中文字幕视频在线看片 | 国产欧美亚洲国产| 亚洲欧洲国产日韩| 简卡轻食公司| 伦理电影大哥的女人| 成人国产麻豆网| 精品久久久精品久久久| 免费人成在线观看视频色| 最近中文字幕2019免费版| 亚州av有码| 少妇熟女欧美另类| 国产黄片美女视频| 午夜激情久久久久久久| 女性被躁到高潮视频| 十八禁网站网址无遮挡 | 亚洲成色77777| 国产国拍精品亚洲av在线观看| 国产视频首页在线观看| 亚洲不卡免费看| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级 | 乱系列少妇在线播放| 国产探花极品一区二区| 丰满人妻一区二区三区视频av| 青春草国产在线视频| 亚洲久久久国产精品| 高清日韩中文字幕在线| 国产中年淑女户外野战色| 国产白丝娇喘喷水9色精品| 最近最新中文字幕免费大全7| 丰满少妇做爰视频| 人妻制服诱惑在线中文字幕| 日本wwww免费看| 一级毛片我不卡| 超碰97精品在线观看| 男人添女人高潮全过程视频| 成人高潮视频无遮挡免费网站| 久久久午夜欧美精品| 国产成人精品婷婷| 久久精品夜色国产| 久久精品熟女亚洲av麻豆精品| 波野结衣二区三区在线| 国内揄拍国产精品人妻在线| 大片免费播放器 马上看| 国产无遮挡羞羞视频在线观看| h视频一区二区三区| 全区人妻精品视频| 久久99蜜桃精品久久| 亚洲国产精品专区欧美| 久久久久人妻精品一区果冻| 九色成人免费人妻av| 一级毛片久久久久久久久女| 午夜福利影视在线免费观看| 久久99热6这里只有精品| 国产成人aa在线观看| 亚洲精品乱码久久久v下载方式| 成人综合一区亚洲| 人妻制服诱惑在线中文字幕| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| a级毛色黄片| 少妇高潮的动态图| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 人妻制服诱惑在线中文字幕| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 亚洲美女视频黄频| 国产v大片淫在线免费观看| 久久99精品国语久久久| 校园人妻丝袜中文字幕| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看| 成人二区视频| 免费观看a级毛片全部| 麻豆乱淫一区二区| 午夜福利在线在线| 青春草国产在线视频| 国产乱人视频| 亚洲精品第二区| 在线看a的网站| 能在线免费看毛片的网站| 看十八女毛片水多多多| 久久 成人 亚洲| 不卡视频在线观看欧美| 国产精品精品国产色婷婷| 精品一区二区三卡| 亚洲精品色激情综合| 久久久久久人妻| freevideosex欧美| 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看| 亚洲,一卡二卡三卡| 久久久精品免费免费高清| av又黄又爽大尺度在线免费看| 精品99又大又爽又粗少妇毛片| av播播在线观看一区| 欧美精品亚洲一区二区| 精品人妻偷拍中文字幕| 欧美精品亚洲一区二区| 成人黄色视频免费在线看| 建设人人有责人人尽责人人享有的 | 22中文网久久字幕| 水蜜桃什么品种好| 国产永久视频网站| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 免费大片18禁| 国产欧美日韩精品一区二区| 亚洲av福利一区| 亚洲国产精品国产精品| 日本vs欧美在线观看视频 | 三级经典国产精品| 高清不卡的av网站| 亚洲最大成人中文| 久久97久久精品| 99re6热这里在线精品视频| 高清av免费在线| 日本欧美视频一区| 国产男人的电影天堂91| 久久国产精品男人的天堂亚洲 | 性色avwww在线观看| 麻豆成人av视频| 国产综合精华液| 人妻一区二区av| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 久久99热6这里只有精品| 高清黄色对白视频在线免费看 | 中文字幕精品免费在线观看视频 | 一区二区av电影网| 国产久久久一区二区三区| 国产极品天堂在线| 99久国产av精品国产电影| 欧美+日韩+精品| 午夜免费鲁丝| 免费看不卡的av| 亚洲国产av新网站| 人妻系列 视频| 亚洲精品,欧美精品| 中文乱码字字幕精品一区二区三区| 精品酒店卫生间| 欧美另类一区| 老女人水多毛片| 黑人高潮一二区| 亚洲美女搞黄在线观看| 国产精品成人在线| 97精品久久久久久久久久精品| 国产女主播在线喷水免费视频网站| 好男人视频免费观看在线| 91狼人影院| 色哟哟·www| 国产av一区二区精品久久 | 老司机影院成人| 免费大片黄手机在线观看| 高清不卡的av网站| 十八禁网站网址无遮挡 | 国产精品秋霞免费鲁丝片| 亚洲精品自拍成人| 久久久久久久久大av| 国产欧美日韩一区二区三区在线 | 亚洲欧美一区二区三区国产| 久久精品人妻少妇| 91精品国产国语对白视频| 岛国毛片在线播放| 综合色丁香网| 人妻制服诱惑在线中文字幕| 精品一区二区免费观看| 超碰97精品在线观看| 日本黄色日本黄色录像| 精品一区在线观看国产| 在线观看人妻少妇| 国产精品免费大片| 高清毛片免费看| 一区二区av电影网| 纯流量卡能插随身wifi吗| 欧美bdsm另类| 91精品国产九色| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产| 蜜臀久久99精品久久宅男| av不卡在线播放| 日本黄大片高清| av国产久精品久网站免费入址| 秋霞伦理黄片| 狂野欧美激情性xxxx在线观看| 成年女人在线观看亚洲视频| 高清av免费在线| 免费黄频网站在线观看国产| 2018国产大陆天天弄谢| 免费少妇av软件| 亚洲欧洲日产国产| 日韩视频在线欧美| 国内揄拍国产精品人妻在线| 极品少妇高潮喷水抽搐| 免费观看av网站的网址| 国产精品久久久久久精品电影小说 | 国产成人精品婷婷| 高清毛片免费看| 最近最新中文字幕免费大全7| 日韩视频在线欧美| 亚洲国产av新网站| 又爽又黄a免费视频| 精品午夜福利在线看| 国产欧美日韩一区二区三区在线 | 亚洲第一av免费看| 国内揄拍国产精品人妻在线| 涩涩av久久男人的天堂| 七月丁香在线播放| 国产欧美日韩精品一区二区| 亚洲国产精品一区三区| 中国美白少妇内射xxxbb| 国产av国产精品国产| 久久久午夜欧美精品| 18禁裸乳无遮挡免费网站照片| 天堂中文最新版在线下载| 2021少妇久久久久久久久久久| 蜜桃在线观看..| 国产av国产精品国产| 国产探花极品一区二区| 色视频在线一区二区三区| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 18+在线观看网站| 亚洲色图综合在线观看| 九草在线视频观看| 久久国内精品自在自线图片| 大香蕉久久网| 精品熟女少妇av免费看| 一级毛片我不卡| 亚洲成人av在线免费| 亚洲婷婷狠狠爱综合网| 日韩欧美 国产精品| 菩萨蛮人人尽说江南好唐韦庄| 丝瓜视频免费看黄片| 久久国产乱子免费精品| 97精品久久久久久久久久精品| 日本与韩国留学比较| 在线免费十八禁| 久久精品国产亚洲av涩爱| 久久久久人妻精品一区果冻| 尾随美女入室| 18禁在线播放成人免费| 成人亚洲精品一区在线观看 | 日本色播在线视频| 亚洲精品日本国产第一区| 日韩人妻高清精品专区| 女性被躁到高潮视频| 久久99蜜桃精品久久| 天天躁日日操中文字幕| 成人一区二区视频在线观看| 国产精品爽爽va在线观看网站| 久久这里有精品视频免费| 成人影院久久| 亚洲不卡免费看| 亚洲天堂av无毛| 久久久a久久爽久久v久久| 久久久久久伊人网av| 尤物成人国产欧美一区二区三区| 欧美bdsm另类| 久久精品熟女亚洲av麻豆精品| 少妇的逼好多水| 美女视频免费永久观看网站| 色吧在线观看| kizo精华| 青青草视频在线视频观看| 欧美区成人在线视频| 在线免费观看不下载黄p国产| 国产精品国产三级专区第一集| 女性被躁到高潮视频| 国产真实伦视频高清在线观看| 国产精品伦人一区二区| 亚洲欧美成人综合另类久久久| 能在线免费看毛片的网站| 亚洲精品亚洲一区二区| 日本黄色片子视频| 免费高清在线观看视频在线观看| 久久99热6这里只有精品| 99热网站在线观看| 国产成人精品一,二区| 在线观看国产h片| 国产淫语在线视频| 女人十人毛片免费观看3o分钟| 日韩国内少妇激情av| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 久久久久久久久大av| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看| 亚洲成人中文字幕在线播放| 国内精品宾馆在线| 黄片无遮挡物在线观看| 乱系列少妇在线播放| 日韩av在线免费看完整版不卡| 青春草国产在线视频| www.色视频.com| 少妇精品久久久久久久| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久影院| 少妇精品久久久久久久| 中文精品一卡2卡3卡4更新| 建设人人有责人人尽责人人享有的 | 午夜视频国产福利| 日韩不卡一区二区三区视频在线| 成人午夜精彩视频在线观看| 免费av中文字幕在线| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 肉色欧美久久久久久久蜜桃| 在线 av 中文字幕| 久久久午夜欧美精品| 青春草国产在线视频| 高清黄色对白视频在线免费看 | 各种免费的搞黄视频| 在线观看人妻少妇| 校园人妻丝袜中文字幕| 三级国产精品片| 少妇精品久久久久久久| 亚洲一区二区三区欧美精品| 男人狂女人下面高潮的视频| 熟女电影av网| 欧美人与善性xxx| 国产精品一及| 人人妻人人爽人人添夜夜欢视频 | 免费观看在线日韩| 日韩免费高清中文字幕av| 制服丝袜香蕉在线| av女优亚洲男人天堂| 搡女人真爽免费视频火全软件| 99久久人妻综合| 99九九线精品视频在线观看视频| 国产一区二区三区综合在线观看 | 97在线人人人人妻| 日本免费在线观看一区| 国产老妇伦熟女老妇高清| 久久久久精品久久久久真实原创| 国产高清三级在线| 国产黄片视频在线免费观看| 国产欧美日韩精品一区二区| 国产精品无大码| 欧美日韩亚洲高清精品| 免费不卡的大黄色大毛片视频在线观看| 99久久精品一区二区三区| 搡老乐熟女国产| 国产精品一区二区三区四区免费观看| 一二三四中文在线观看免费高清| 夜夜看夜夜爽夜夜摸| 欧美日韩亚洲高清精品| 午夜免费男女啪啪视频观看| 欧美日韩视频高清一区二区三区二| 网址你懂的国产日韩在线| 嫩草影院入口| 国产爱豆传媒在线观看| 亚洲精品第二区| 一区二区三区精品91| 联通29元200g的流量卡| 美女福利国产在线 | 成人毛片60女人毛片免费| 亚洲人与动物交配视频| 欧美日韩精品成人综合77777| 99热这里只有精品一区| 国产精品精品国产色婷婷| 国产成人午夜福利电影在线观看| 色婷婷av一区二区三区视频| 人妻 亚洲 视频| 91精品国产国语对白视频| 蜜臀久久99精品久久宅男| 精品久久久久久电影网| 秋霞伦理黄片| 综合色丁香网| 亚洲激情五月婷婷啪啪| 亚洲国产毛片av蜜桃av| 中文精品一卡2卡3卡4更新| www.av在线官网国产| 一区在线观看完整版|