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

    基于DEM的碎屑流運動特性數(shù)值模擬

    2017-03-15 10:54:35付成華
    水利水電科技進(jìn)展 2017年2期
    關(guān)鍵詞:堆積體碎屑滑坡

    趙 川,付成華

    (西華大學(xué)能源與動力工程學(xué)院,四川 成都 610039)

    基于DEM的碎屑流運動特性數(shù)值模擬

    趙 川,付成華

    (西華大學(xué)能源與動力工程學(xué)院,四川 成都 610039)

    采用顆粒離散元法建立了碎屑流三維數(shù)值模型,分別從整體和局部對碎屑流運動特性進(jìn)行深入分析。結(jié)果表明:碎屑流的運動過程可分為啟動加速、高速滑動和減速堆積3個階段;邊坡底面平整時,堆積體中大顆粒堆積在表面,而細(xì)小顆粒位于中下部位;碎屑流運動速度和能量損失受顆粒間的接觸碰撞影響大;碎屑流不同部位運動速度最大值相同,達(dá)到速度峰值的時間與顆粒的初始水平位置有關(guān),與深度無關(guān);堆積體分區(qū)明顯,堆積體中的顆粒按初始水平位置依次排布;坡腳位置堆積深度最大,往兩側(cè)逐漸減小;碎屑流在運動初期和末期橫向運動明顯,擴(kuò)大了災(zāi)害影響范圍。

    碎屑流;三維顆粒離散元法;碰撞接觸;堆積體;數(shù)值試驗

    碎屑流是一種滑坡巖體在初始狀態(tài)或者高速運動過程中轉(zhuǎn)化為碎屑物質(zhì)的運動形式,其某些運動特性與流體相似。碎屑流是一種常見的自然災(zāi)害,有運動速度快、攜帶巖體量大和破壞力強等特點,對人民生命財產(chǎn)安全構(gòu)成巨大威脅。我國地形條件復(fù)雜,高山丘陵廣泛分布,發(fā)生地質(zhì)災(zāi)害的頻率很高,尤其是我國西南地區(qū)在遭受了汶川地震后,邊坡表層巖體破碎嚴(yán)重,更是碎屑流災(zāi)害的高發(fā)地帶[1]。黃潤秋[2]分析了20世紀(jì)以來中國的大型滑坡,將其分為巖質(zhì)滑坡、土層滑坡和松散堆積層滑坡3類。本文研究的碎屑流不考慮顆粒單元間的黏結(jié)作用,屬于松散堆積層滑坡。

    20世紀(jì)以來,國內(nèi)外學(xué)者對顆粒碎屑物質(zhì)進(jìn)行了大量研究,《Science》于2005年把顆粒物質(zhì)列為100個科學(xué)難題之一[3-4];Iverson等[5]基于復(fù)雜的三維地形,對沙粒崩塌(granular avalanche)進(jìn)行了室內(nèi)物理實驗,并采用連續(xù)介質(zhì)模型進(jìn)行模擬,模擬結(jié)果與實驗結(jié)果基本吻合;Hussin等[6]以MassMov2D為工具,通過反分析確定了某滑坡帶的計算參數(shù),對潛在的滑坡進(jìn)行了預(yù)測;王玉峰等[7]對汶川地震觸發(fā)的高速遠(yuǎn)程碎屑流堆積反粒序特征進(jìn)行了分析,揭示了碎屑流堆積體內(nèi)部所特有的堆積規(guī)律;張明等[8]對目前碎屑流的研究成果進(jìn)行了系統(tǒng)的歸納和總結(jié)。總的看來,已有的碎屑流研究對象主要針對室內(nèi)小尺寸模型或已經(jīng)發(fā)生的具體滑坡。具體的滑坡案例受地形條件影響大,室內(nèi)模型難以反映真實高速碎屑流的某些運動特性,其研究結(jié)論難以反映一般規(guī)律[9-12]。

    為了解決上述碎屑流研究存在的問題,本文采用三維顆粒離散元法(discrete element method,DEM)建立了標(biāo)準(zhǔn)的碎屑流數(shù)值模型,選用含有棱角的顆粒單元,并通過漏斗模型試驗標(biāo)定計算參數(shù),對碎屑流運動特性進(jìn)行了深入分析探討。

    1 碎屑流數(shù)值模型建立及參數(shù)選取

    1.1 離散元模型構(gòu)建

    碎屑流數(shù)值模型主要由碎屑堆積槽、滑動面和堆積面組成?;瑒用婧投逊e面構(gòu)成一個邊坡系統(tǒng),坡高為10 m,坡比為1∶2,堆積面長為10 m;碎屑堆積槽用于放置顆粒物質(zhì),槽外形為三棱柱體,尺寸為1.5 m×3 m×5 m,如圖1所示。

    圖1 碎屑流離散元模型尺寸(單位:m)

    由于真實巖體碎屑外形并非標(biāo)準(zhǔn)球體,而是含有一定棱角的球體,因此模型試驗中構(gòu)建了含有棱角的顆粒單元,如圖2所示,其中圖2(a)為真實碎屑流塊體,圖2(b)為構(gòu)建的離散元塊體。顆粒大小按高斯分布隨機(jī)生成,直至充滿堆積槽,顆??倲?shù)為3 707,總填充體積22.5 m3,最終顆粒模型試樣如圖3所示。

    圖2 數(shù)值計算顆粒單元外形

    圖3 碎屑流離散元模型

    滑坡碎屑流運動過程中,巖塊之間存在大量的碰撞和接觸,為準(zhǔn)確描述DEM中顆粒單元運動過程的力學(xué)行為,計算采用Hertz-Mindlin無滑動接觸模型,主要包括法向力Fn、法向阻尼力Fnd、切向力Ft和切向阻尼力Ftd4個參數(shù),計算公式[13-14]如下:

    (1)

    (2)

    Ft=-Stδt

    (3)

    (4)

    式中:R*為等效顆粒半徑,m;Y*為等效楊氏模量,Pa;δn、δt分別為兩個顆粒接觸面法向與切向重疊量,m;m*為等效質(zhì)量,kg;vn、vt分別為兩個顆粒接觸面法向和切向相對速度,m/s;β為阻尼力系數(shù);Sn、St分別為兩個顆粒接觸面法向與切向剛度。

    1.2 數(shù)值計算參數(shù)標(biāo)定

    碎屑流的運動受摩擦系數(shù)的影響很大,為校核相關(guān)摩擦系數(shù),對試驗選用的顆粒采用漏斗模型試驗進(jìn)行標(biāo)定,如圖4所示。由于本試驗中碎屑顆粒間不存在黏結(jié)力,其天然休止角即為松散堆積體的內(nèi)摩擦角φ。數(shù)值試驗中滑動面為平面,為模擬邊坡表面覆蓋的松散堆積層[15],基底摩擦系數(shù)取0.42,顆粒間摩擦系數(shù)取0.4。具體方法是首先設(shè)置好基底摩擦系數(shù)和顆粒間摩擦系數(shù),通過不斷調(diào)整滾動摩擦系數(shù)重復(fù)計算,將堆積體的天然休止角為22°(arctan 0.4=arctanH/R=21.8°,其中H為堆積體中心高度,R為堆積體底半徑)時所采用的參數(shù)作為碎屑流模型的參數(shù),最終得到的計算參數(shù)見表1。

    圖4 離散元顆粒漏斗數(shù)值模型試驗

    表1 DEM模型相關(guān)參數(shù)取值

    圖5 不同時刻碎屑流運動速度分布(單位:m/s)

    2 數(shù)值試驗結(jié)果及分析

    2.1 整體運動特性

    為分析巖體碎屑流不同時刻的運動狀態(tài),圖5給出了碎屑流的運動速度分布。在時間t=0 s時,撤掉試樣前端擋板,在重力作用下,顆粒開始沿著坡面向下運動。由于顆粒并不是標(biāo)準(zhǔn)的球形顆粒,而是含有棱角的顆粒單元,所以整個運動過程中主要以平面滑動為主,僅在顆粒流表面存在滾動現(xiàn)象。運動過程中,前端顆粒由于所受約束小,運動速度較大,最大顆粒速度超過7 m/s。t=4 s時,前端顆粒運動到堆積面,并與底面發(fā)生碰撞,動能損失較大,速度逐漸減小,并堆積在坡腳位置。最終形成的堆積角在10°左右,僅為初始坡角26.6°的38%。

    圖6為整個運動過程顆粒物質(zhì)平均速度隨時間的變化,10 s時碎屑的平均速度僅為0.002 m/s,運動基本停止。根據(jù)曲線變化趨勢,可將碎屑流的運動過程劃分為啟動加速、高速滑動和減速堆積3個運動階段,分別對應(yīng)的時間段為0~4 s、4~6 s和6~10 s,加速和減速階段所用的時間大體相當(dāng),均占總耗時的40%左右,而高速滑動僅占總耗時的20%。

    圖6 碎屑流平均速度的變化過程

    另外,從最后的堆積體顆粒分布可以發(fā)現(xiàn),大顆粒位于堆積體表面,細(xì)小顆粒位于堆積體中下部位。究其原因,是由于在運動過程中,顆粒間相互碰撞摩擦,顆粒間存在一定大小的空隙,細(xì)小顆粒沿著空隙下落到堆積體底部,從而形成大顆粒物質(zhì)在表面、細(xì)小顆粒在中下部的現(xiàn)象,這與實際碎屑流形成的堆積序列相符。

    模擬時不考慮巖體間的黏結(jié)力,僅考慮相互碰撞作用。圖7為碎屑流運動過程中顆粒單元相互接觸數(shù)量隨時間的變化,初始時刻堆積較為緊密,總的接觸數(shù)量大約有11 500;隨著運動速度的增大接觸數(shù)量逐漸減小,在t=5.4 時接觸數(shù)量減小到最小值1 841,說明此時顆粒最為稀疏,顆粒間相互作用最少;此后隨著顆粒流前緣堆積體逐漸形成,接觸數(shù)量又逐漸增加,最后可達(dá)13 500左右,略大于初始時刻的數(shù)量。對比圖6和圖7可見,碎屑流的平均速度隨時間的變化曲線與單元接觸數(shù)量隨時間的變化曲線的變化趨勢正好相反,且速度的最大值正好對應(yīng)接觸數(shù)的最小值,均發(fā)生在t=5.4 s時,說明顆粒越稀疏,運動速度越大,顆粒越密集意味著更多的碰撞,而碰撞將導(dǎo)致碎屑流內(nèi)部能量大量耗散,故運動速度也隨之減小。

    圖7 顆粒單元接觸數(shù)量的變化過程

    圖8為單元的運動跡線和速度,可以發(fā)現(xiàn)大部分單元的速度分布基本是連續(xù)變化的,表現(xiàn)出一定的流體性質(zhì),但少數(shù)表面不連續(xù)的紅色線條顯示出存在顆粒單元在某個時段內(nèi)脫離碎屑流運動,而單獨以較大的速度在表面發(fā)生跳動,這是由于表面部分高速運動的顆粒之間發(fā)生偏心碰撞和回彈造成的。

    圖8 碎屑流DEM模型單元運動跡線及速度(單位:m/s)

    2.2 局部運動特性分析

    為分析碎屑流不同部位在運動過程中的差異,將模型試樣分成5個區(qū)域,每部分顆??倲?shù)占總顆粒數(shù)的20%左右,并分別編號A、B、C、D、E,如圖9所示,其中A區(qū)代表后部,B區(qū)代表中上部,C區(qū)代表中下部,D區(qū)代表前端上部,E區(qū)代表前端下部。圖10為5個區(qū)域內(nèi)顆粒的運動速度隨時間的變化過程,可以發(fā)現(xiàn),位于前端的C和E區(qū)的顆粒最先啟動,在前0.5 s內(nèi),曲線斜率很大,說明速度增長很

    圖9 顆粒分區(qū)

    圖10 不同分區(qū)顆粒的運動速度

    快,這是由于前端顆粒在初始時刻存在掉落現(xiàn)象,隨后曲線斜率保持穩(wěn)定,在t=4.5 s時速度達(dá)到最大值;而位于中部的B和D區(qū)顆粒啟動稍晚,在加速運動過程中,速度增長保持穩(wěn)定,在t=6.5 s左右速度達(dá)到峰值;位于最后端的A區(qū)則最后啟動,且初始時刻速度增加幅度不大,這是因為中部的B區(qū)和D區(qū)顆粒阻擋和延遲了A區(qū)顆粒的運動,隨后速度逐漸增加,并在t=8 s時速度達(dá)到最大值。綜合分析可見,5個區(qū)域的顆粒速度在不同時刻達(dá)到最大值,且最大值基本相同,約為4.5 m/s;B區(qū)顆粒速度稍小,為4.4 m/s。

    圖11為5個區(qū)域內(nèi)顆粒在水平方向x、豎直方向y、橫向z的運動速度分量,可以發(fā)現(xiàn)x方向的速度最大,最大峰值速度約4 m/s;y方向的速度次之,最大速度約2 m/s;z方向的速度最小,最大也僅有0.04 m/s。各部分顆粒的x和y方向速度幾乎同時達(dá)到最大值,且曲線變化趨勢也基本相同,其中C和E區(qū)顆粒在前0.5 s內(nèi)y方向速度快速增大,這是顆粒突然釋放且前緣部分受到阻力很小造成的(圖11(b))。z方向顆粒運動速度一直較小,具體分析曲線變化趨勢可以發(fā)現(xiàn):在碎屑流運動初期,位于前緣的C和E區(qū)顆粒橫向運動較為明顯;運動中期各部分顆粒z方向運動速度均很小,此階段內(nèi)顆粒堆積密度最小,相互之間的碰撞擠壓也最少;在最后的堆積期間,A、C和D區(qū)的顆粒橫向運動速度相對大一些。總體來說,在運動初期和最后堆積期,顆粒沿運動方向速度較小時,z方向的速度相對較大;而當(dāng)顆粒沿運動方向高速運動時,z方向速度相對較小(圖11(c))。

    圖12為碎屑流最終的堆積狀態(tài),A區(qū)顆粒完全堆積在坡面上,B和D區(qū)大部分堆積在坡面上,極少部分堆積在坡腳以下; C和E區(qū)的顆粒完全堆積在底面上,堆積長度較大,越往前堆積厚度越小;坡腳位置的堆積厚度最大,最大堆積厚度可達(dá)1 m。

    圖12 碎屑流最終堆積狀態(tài)

    3 結(jié) 論

    a. 根據(jù)碎屑流速度時間曲線,可將其運動過程劃分為3個運動階段:啟動加速段、高速滑動段和減速堆積段。

    b. 碎屑流在平整坡面上運動時,存在大量接觸和碰撞,小顆粒不斷通過空隙向下運動,最終造成了明顯的堆積分布序列:大顆粒位于堆積體表面,小顆粒位于堆積體底部,在坡腳位置堆積深度最大;碎屑流的速度隨時間的變化曲線與單元接觸數(shù)量隨時間變化曲線的變化趨勢恰好相反,且速度的最大值正好對應(yīng)接觸數(shù)量的最小值。

    c. 不同區(qū)域內(nèi)顆粒速度最大值基本相同,但達(dá)到最大值的時刻有所不同,這與顆粒的初始水平位置相關(guān),而與其深度無關(guān);各區(qū)域內(nèi)顆粒水平和豎直方向運動速度較大,且變化趨勢基本同步,而橫向速度相對較小;堆積體分區(qū)明顯,堆積體中的顆粒按初始水平位置依次排布。

    d. 碎屑流沿運動方向高速運動時,橫向速度很小;在碎屑流運動初期和堆積末期,水平和豎直方向運動速度較小,橫向速度較大,碎屑流橫向運動更為明顯,擴(kuò)大了災(zāi)害橫向影響范圍。

    [1] 鄭穎人,陳祖煜,王恭先,等.邊坡與滑坡工程治理[M].北京:人民交通出版社,2007.

    [2] 黃潤秋.20世紀(jì)以來中國的大型滑坡及其發(fā)生機(jī)制[J].巖石力學(xué)與工程學(xué)報,2007,26(3):433-454.(HUANG Runqiu.Large-scale landslides and their sliding mechanisms in China since the 20th century[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(3): 433-433.(in Chinese))

    [3] HEWITT K.Catastrophic landslide deposits in the Karakoram Himalaya[J].Science,1988,242:64-67.

    [4] 孫其誠,王光謙.顆粒物質(zhì)力學(xué)導(dǎo)論[M].北京:科學(xué)出版社,2009.

    [5] IVERSON R M,LOGAN M,DENLINGER R P.Granular avalanches across irregular three-dimensional terrain: 2.experimental tests[J].Journal of Geophysical Research: Earth Surface,2004,109(1):337-357.

    [6] HUSSIN H Y,CIUREAN R,FRIGERIO S,et al.Assessing the effect of mitigation measures on landslide hazard using 2D numerical runout modeling[J].Landslide Science for a Safer Geoenvironment,2014,2(6):679-684.

    [7] 王玉峰,程謙恭,朱圻.汶川地震觸發(fā)高速遠(yuǎn)程滑坡-碎屑流堆積反粒序特征及機(jī)制分析[J].巖石力學(xué)與工程學(xué)報,2012,31(6): 1089-1106.(WANG Yufeng,CHENG Qiangong,ZHU Qi.Inverse grading analysis of deposit from rock avalanches triggered by Wenchuan Earthquake[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(6): 1089-1106.(in Chinese))

    [8] 張明,殷躍平,吳樹仁,等.高速遠(yuǎn)程碎屑流運動機(jī)理研究發(fā)展現(xiàn)狀與展望[J].工程地質(zhì)學(xué)報,2010,18(6):805-817.(ZHANG Ming,YIN Yueping,WU Shuren,et al.Development status and prospects of studies on kinematics of long runout rock avalanches[J].Journal of Engineering Geology,2010,18(6): 805-817.(in Chinese))

    [9] BANTON J,VILLARD P,JONGMANS D,et al.Two-dimensional discrete element models of debris avalanches: parameterization and the reproducibility of experimental results[J].Journal of Geophysical Research: Earth Surface,2009,114(4):60-71.

    [10] LO Chaiming,LIN Minglang,TANG Chaolung,et al.A kinematic model of the Hsiaolin landslide calibrated to the morphology of the landslide deposit[J].Engineering Geology,2011,123(1):22-39.

    [11] HARALD T,WANG Y,CHIOU M C,et al.Flow-obstacle interaction in rapid granular avalanches: DEM simulation and comparison with experiment[J].Granular Matter,2009,11(4):209-220.

    [12] 畢鈺璋,付躍升,何思明,等.牛眠溝地震滑坡碎屑化全過程離散元模擬[J].中國地質(zhì)災(zāi)害與防治學(xué)報,2015,26(3):17-25.(BI Yuzhang,FU Yuesheng,HE Siming,et al.Simulation of the whole process of Niumiangou creek rock avalanche triggered by the earthquake using a distinct element method[J].The Chinese Journal of Geological Hazard and Control,2015,26(3):17-25.(in Chinese))

    [13] 徐泳,孫其誠,張凌,等.顆粒離散法研究進(jìn)展[J].力學(xué)進(jìn)展,2003,33(2):251-260.(XU Yong,SUN Qicheng,ZHANG Ling,et al.Advances in discrete element methods for particulate materials[J].Advances in Mechanics,2003,33 (2): 251-260.(in Chinese))

    [14] 王國強,郝萬軍,王繼新.離散單元法及其在EDEM上的實踐[M].西安: 西北工業(yè)大學(xué)出版社,2010.

    [15] 錢家歡,殷宗澤.土工原理與計算[M].2版.北京:中國水利水電出版社,1996.

    Numerical simulation of movement characteristics of debris flow based on DEM

    ZHAO Chuan, FU Chenghua

    (SchoolofEnergyandPowerEngineering,XihuaUniversity,Chengdu610039,China)

    A three-dimensional numerical model for debris flow was established with the particle discrete element method (DEM), and the movement characteristics of debris flow were analyzed from the full and local aspects. The results show that the movement process of debris flow can be divided into three stages, including starting acceleration, high-speed sliding, and reduced-speed deposit. When the bottom of the slope is flat, large particles are found on the surface of the deposition body, but small particles are in the middle and low areas. The velocity and energy loss of debris flow is significantly influenced by the contacts and collisions between particles. The maximum speeds of different parts of the debris flow are almost the same, and the time required to reach the speed peak is related to the initially horizontal position of particles, but not related to the deposit depth. Area division is distinct in the deposition body, and the particles of the deposition body are arranged according to their initially horizontal positions. The maximum deposit height is at the slope toe, and the deposit height gradually decreases towards both sides. Significant lateral movement occurs in the debris flow in the early and late stages, which can expand the scope of a disaster.

    debris flow; 3-D particle discrete element; numerical experiment; contact collision; deposition body

    西華大學(xué)流體及動力機(jī)械教育部重點實驗室(重點研究基地)開放課題(szjj2015-039)

    趙川(1989—),男,碩士研究生,主要從事巖土工程穩(wěn)定性研究。E-mail:zhaochuanvip@163.com

    付成華(1978—),女,副教授,博士,主要從事水利水電工程研究。E-mail:fuchh_xhu@163.com

    10.3880/j.issn.1006-7647.2017.02.008

    P642.21

    :A

    :1006-7647(2017)02-0043-05

    2016-02-26 編輯:鄭孝宇)

    猜你喜歡
    堆積體碎屑滑坡
    Sweet Carrots甜甜的胡蘿卜
    一種具有塑料碎屑收集機(jī)構(gòu)的水磨床
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    后退式注漿技術(shù)在隧道突涌堆積體加固處治中的應(yīng)用
    隧道黃土堆積體施工技術(shù)
    清潔攪拌機(jī)的小妙招
    伴侶(2018年2期)2018-02-28 20:38:24
    渤中X油田沙一段生物碎屑白云巖井場識別方法
    錄井工程(2017年4期)2017-03-16 06:10:28
    大型堆積體滑坡治理施工技術(shù)研究
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    基于Fluent的滑坡入水過程數(shù)值模擬
    在线观看免费午夜福利视频| 久久精品91蜜桃| 午夜精品一区二区三区免费看| 色av中文字幕| 级片在线观看| 桃色一区二区三区在线观看| 久久草成人影院| 12—13女人毛片做爰片一| 香蕉久久夜色| 午夜视频精品福利| 免费一级毛片在线播放高清视频| 精品一区二区三区av网在线观看| 一级毛片高清免费大全| 免费观看精品视频网站| 国产亚洲精品综合一区在线观看| 91九色精品人成在线观看| 天堂网av新在线| 美女被艹到高潮喷水动态| 久久热在线av| 亚洲一区高清亚洲精品| 国产麻豆成人av免费视频| 国产成人av教育| 国产成人欧美在线观看| 国产精品 欧美亚洲| 国产高清三级在线| 真人做人爱边吃奶动态| 99热精品在线国产| 99久久精品一区二区三区| 日韩国内少妇激情av| 婷婷精品国产亚洲av| 精品福利观看| 动漫黄色视频在线观看| 国产亚洲av高清不卡| 亚洲av五月六月丁香网| 免费在线观看影片大全网站| 精品福利观看| 全区人妻精品视频| 女同久久另类99精品国产91| 欧美日韩黄片免| 亚洲真实伦在线观看| 精华霜和精华液先用哪个| 国产三级在线视频| 国产精品久久久av美女十八| 1024香蕉在线观看| 91在线精品国自产拍蜜月 | 在线观看午夜福利视频| 女警被强在线播放| 黄色丝袜av网址大全| 中文资源天堂在线| 99久久精品一区二区三区| 亚洲国产精品合色在线| 国内精品美女久久久久久| av国产免费在线观看| 18禁黄网站禁片免费观看直播| 国产一区在线观看成人免费| 日韩三级视频一区二区三区| 色播亚洲综合网| 舔av片在线| 亚洲精品美女久久久久99蜜臀| 深夜精品福利| 久久久久久久久久黄片| 精品久久久久久久毛片微露脸| 国产精华一区二区三区| xxxwww97欧美| 真实男女啪啪啪动态图| 超碰成人久久| 欧美日韩精品网址| 国产99白浆流出| 热99在线观看视频| 91av网站免费观看| 亚洲精品久久国产高清桃花| 久久久国产成人免费| 精品日产1卡2卡| 国产综合懂色| 中文资源天堂在线| 久久精品人妻少妇| 天堂√8在线中文| 色播亚洲综合网| 99热6这里只有精品| 欧美日韩亚洲国产一区二区在线观看| 校园春色视频在线观看| 久久久久精品国产欧美久久久| 日日摸夜夜添夜夜添小说| 久久国产精品影院| 嫩草影视91久久| 国产成人系列免费观看| 12—13女人毛片做爰片一| 人人妻人人澡欧美一区二区| 免费av毛片视频| 变态另类丝袜制服| 国产高清有码在线观看视频| 18禁裸乳无遮挡免费网站照片| 很黄的视频免费| 精品久久久久久成人av| 真实男女啪啪啪动态图| 亚洲国产精品合色在线| 欧美三级亚洲精品| 女人高潮潮喷娇喘18禁视频| 免费电影在线观看免费观看| cao死你这个sao货| 欧美日韩一级在线毛片| 亚洲第一电影网av| 不卡av一区二区三区| 午夜两性在线视频| 久久人人精品亚洲av| 中文字幕人成人乱码亚洲影| 听说在线观看完整版免费高清| 狂野欧美白嫩少妇大欣赏| 午夜福利成人在线免费观看| 久久久久久久久久黄片| 最好的美女福利视频网| 亚洲国产精品sss在线观看| 国产欧美日韩精品亚洲av| 狠狠狠狠99中文字幕| 国产视频内射| 日日夜夜操网爽| 欧美激情久久久久久爽电影| 村上凉子中文字幕在线| 丰满人妻一区二区三区视频av | 午夜福利18| 亚洲中文av在线| 99riav亚洲国产免费| 免费在线观看成人毛片| 欧美成狂野欧美在线观看| 首页视频小说图片口味搜索| 一本精品99久久精品77| 欧美绝顶高潮抽搐喷水| 国产一区二区三区视频了| 18禁黄网站禁片免费观看直播| 在线免费观看的www视频| а√天堂www在线а√下载| 久久久国产成人免费| 色综合欧美亚洲国产小说| 国产精品电影一区二区三区| 欧美日本视频| 国产亚洲精品久久久久久毛片| 久久国产乱子伦精品免费另类| 久久中文看片网| 少妇裸体淫交视频免费看高清| 久久午夜亚洲精品久久| 国产精品永久免费网站| 一级毛片女人18水好多| 激情在线观看视频在线高清| 日韩免费av在线播放| 午夜福利欧美成人| 亚洲国产欧美人成| 99国产精品一区二区蜜桃av| 99热6这里只有精品| 国产欧美日韩精品一区二区| 午夜福利在线在线| 亚洲精品久久国产高清桃花| 一级毛片女人18水好多| 最好的美女福利视频网| 99国产综合亚洲精品| 成人无遮挡网站| 亚洲五月天丁香| 午夜激情欧美在线| 国内精品久久久久精免费| 男人舔奶头视频| 日本a在线网址| 午夜久久久久精精品| 免费在线观看影片大全网站| 国产伦人伦偷精品视频| 伦理电影免费视频| 伊人久久大香线蕉亚洲五| 91在线精品国自产拍蜜月 | 亚洲精品在线观看二区| 欧美在线一区亚洲| 久久久精品大字幕| av欧美777| 香蕉国产在线看| 国产午夜精品论理片| 黄色女人牲交| 女人被狂操c到高潮| 麻豆成人av在线观看| 久久精品影院6| www国产在线视频色| 女同久久另类99精品国产91| 少妇裸体淫交视频免费看高清| 国产人伦9x9x在线观看| 97人妻精品一区二区三区麻豆| 这个男人来自地球电影免费观看| 亚洲自偷自拍图片 自拍| 美女黄网站色视频| 99久久综合精品五月天人人| 午夜福利高清视频| 中国美女看黄片| 国产欧美日韩精品一区二区| 美女午夜性视频免费| 午夜精品在线福利| 精品福利观看| 日韩欧美免费精品| 怎么达到女性高潮| 精华霜和精华液先用哪个| 成年人黄色毛片网站| 国产伦在线观看视频一区| 久久久成人免费电影| 久99久视频精品免费| 岛国在线观看网站| 亚洲18禁久久av| 日本黄色片子视频| 日韩免费av在线播放| 日韩av在线大香蕉| netflix在线观看网站| 九九久久精品国产亚洲av麻豆 | 亚洲人成网站在线播放欧美日韩| 观看美女的网站| 黑人欧美特级aaaaaa片| 村上凉子中文字幕在线| 女警被强在线播放| 999久久久精品免费观看国产| 少妇裸体淫交视频免费看高清| 日本一本二区三区精品| 国产精品香港三级国产av潘金莲| 欧美色欧美亚洲另类二区| 欧美日韩综合久久久久久 | 啦啦啦韩国在线观看视频| 精华霜和精华液先用哪个| 麻豆国产av国片精品| 亚洲国产精品合色在线| 黄色视频,在线免费观看| 制服人妻中文乱码| 日韩 欧美 亚洲 中文字幕| 午夜精品久久久久久毛片777| 国产精品一区二区三区四区免费观看 | 中文字幕人成人乱码亚洲影| x7x7x7水蜜桃| 国产精品综合久久久久久久免费| 久久久久九九精品影院| 激情在线观看视频在线高清| 在线观看美女被高潮喷水网站 | 最近在线观看免费完整版| 在线十欧美十亚洲十日本专区| 国产高清视频在线观看网站| 亚洲国产精品久久男人天堂| 亚洲成人中文字幕在线播放| 99热6这里只有精品| 国产视频一区二区在线看| 亚洲无线在线观看| 高清在线国产一区| 久久久国产欧美日韩av| 亚洲成人久久爱视频| 国产精品99久久99久久久不卡| 淫秽高清视频在线观看| 国产激情欧美一区二区| 深夜精品福利| 亚洲天堂国产精品一区在线| 精品国产三级普通话版| 91麻豆av在线| 无遮挡黄片免费观看| 精品一区二区三区视频在线 | 国产精品久久久人人做人人爽| xxxwww97欧美| 欧美+亚洲+日韩+国产| a在线观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久毛片微露脸| 两个人看的免费小视频| 国产伦精品一区二区三区视频9 | 法律面前人人平等表现在哪些方面| 久久欧美精品欧美久久欧美| 日韩有码中文字幕| 三级国产精品欧美在线观看 | 国产亚洲精品综合一区在线观看| 哪里可以看免费的av片| 国产成人av激情在线播放| 一区福利在线观看| 亚洲第一欧美日韩一区二区三区| 日韩免费av在线播放| 18禁裸乳无遮挡免费网站照片| 日本熟妇午夜| 性色av乱码一区二区三区2| 日日夜夜操网爽| 1024香蕉在线观看| 国内久久婷婷六月综合欲色啪| 久久国产乱子伦精品免费另类| 午夜精品久久久久久毛片777| 一级黄色大片毛片| 69av精品久久久久久| 午夜日韩欧美国产| 国产一区二区在线av高清观看| 精品久久久久久久久久久久久| 一级a爱片免费观看的视频| 无人区码免费观看不卡| 亚洲成人久久性| 在线观看午夜福利视频| av福利片在线观看| 亚洲国产精品成人综合色| 亚洲中文日韩欧美视频| 19禁男女啪啪无遮挡网站| 伦理电影免费视频| 日韩欧美一区二区三区在线观看| 婷婷丁香在线五月| av国产免费在线观看| aaaaa片日本免费| 99久久久亚洲精品蜜臀av| 国产精品久久视频播放| 日韩精品中文字幕看吧| 欧美激情在线99| 久久精品亚洲精品国产色婷小说| 国产精品爽爽va在线观看网站| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 国产成人aa在线观看| 亚洲 欧美 日韩 在线 免费| 最近最新中文字幕大全免费视频| 91九色精品人成在线观看| 久久久久久国产a免费观看| 精品一区二区三区av网在线观看| 亚洲美女视频黄频| 两性夫妻黄色片| 国产精品国产高清国产av| 黄色视频,在线免费观看| 夜夜躁狠狠躁天天躁| 亚洲激情在线av| 国产精品久久电影中文字幕| 亚洲片人在线观看| 成人亚洲精品av一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美激情性xxxx| 91字幕亚洲| 久久午夜亚洲精品久久| 麻豆久久精品国产亚洲av| 日韩高清综合在线| 午夜两性在线视频| 青草久久国产| 一二三四在线观看免费中文在| 观看美女的网站| 美女被艹到高潮喷水动态| 97超视频在线观看视频| 天堂网av新在线| 在线免费观看的www视频| 伦理电影免费视频| 国产精品,欧美在线| 亚洲第一电影网av| 在线观看一区二区三区| 午夜成年电影在线免费观看| 色综合亚洲欧美另类图片| 法律面前人人平等表现在哪些方面| av国产免费在线观看| 亚洲一区高清亚洲精品| 一级a爱片免费观看的视频| 首页视频小说图片口味搜索| 97人妻精品一区二区三区麻豆| 国产97色在线日韩免费| 午夜福利18| 夜夜夜夜夜久久久久| 亚洲av成人不卡在线观看播放网| 麻豆成人午夜福利视频| 国产乱人视频| 在线免费观看不下载黄p国产 | 精品久久蜜臀av无| 啦啦啦免费观看视频1| 久久久久国产一级毛片高清牌| 日韩欧美精品v在线| 欧美色欧美亚洲另类二区| 久久久久亚洲av毛片大全| 99热6这里只有精品| 日韩精品青青久久久久久| 亚洲人与动物交配视频| 国产三级中文精品| 国产成人av教育| www.自偷自拍.com| 日本撒尿小便嘘嘘汇集6| 嫩草影视91久久| 午夜激情欧美在线| 村上凉子中文字幕在线| 色综合欧美亚洲国产小说| 国产一区在线观看成人免费| 国产久久久一区二区三区| 久久这里只有精品中国| 一本精品99久久精品77| 亚洲国产欧洲综合997久久,| 国产欧美日韩精品一区二区| 国产一区二区三区在线臀色熟女| 又黄又粗又硬又大视频| 日日干狠狠操夜夜爽| 亚洲国产精品合色在线| 中亚洲国语对白在线视频| 亚洲专区国产一区二区| 国产亚洲av嫩草精品影院| 欧美日韩黄片免| 久9热在线精品视频| 久久久久久九九精品二区国产| 熟女人妻精品中文字幕| 久久欧美精品欧美久久欧美| 性欧美人与动物交配| 97人妻精品一区二区三区麻豆| 亚洲欧美精品综合久久99| 色噜噜av男人的天堂激情| 法律面前人人平等表现在哪些方面| 亚洲美女视频黄频| 两个人看的免费小视频| 国产午夜精品久久久久久| 天堂av国产一区二区熟女人妻| 久久久国产精品麻豆| 99久久99久久久精品蜜桃| 国产精品av久久久久免费| 国产免费男女视频| a级毛片在线看网站| 免费av不卡在线播放| 18美女黄网站色大片免费观看| 久久香蕉精品热| 日本一二三区视频观看| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看亚洲国产| www日本黄色视频网| 99久久国产精品久久久| 黄色 视频免费看| 极品教师在线免费播放| 亚洲七黄色美女视频| 国产成人av激情在线播放| 99re在线观看精品视频| 后天国语完整版免费观看| 国产毛片a区久久久久| 中文字幕熟女人妻在线| 最近最新中文字幕大全电影3| 天天躁狠狠躁夜夜躁狠狠躁| 别揉我奶头~嗯~啊~动态视频| 国产精品1区2区在线观看.| 亚洲av日韩精品久久久久久密| 午夜福利成人在线免费观看| 五月伊人婷婷丁香| 在线观看免费午夜福利视频| 九九在线视频观看精品| 狠狠狠狠99中文字幕| 一区二区三区国产精品乱码| 国产精品久久久久久精品电影| 国产亚洲欧美98| 美女午夜性视频免费| 成人18禁在线播放| 久久精品影院6| 日韩欧美在线乱码| а√天堂www在线а√下载| 黄色日韩在线| 后天国语完整版免费观看| www.自偷自拍.com| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 精品久久久久久久久久免费视频| 色综合站精品国产| 又黄又爽又免费观看的视频| 禁无遮挡网站| 人人妻,人人澡人人爽秒播| 久久草成人影院| 国产蜜桃级精品一区二区三区| 国产高清三级在线| 噜噜噜噜噜久久久久久91| 国产精品亚洲美女久久久| 亚洲黑人精品在线| 老汉色∧v一级毛片| 亚洲人成伊人成综合网2020| 美女午夜性视频免费| 亚洲天堂国产精品一区在线| a级毛片a级免费在线| 看片在线看免费视频| 欧美性猛交╳xxx乱大交人| 亚洲国产精品合色在线| 国产亚洲精品久久久com| ponron亚洲| 首页视频小说图片口味搜索| 欧美午夜高清在线| 麻豆一二三区av精品| 观看美女的网站| 一二三四在线观看免费中文在| 后天国语完整版免费观看| 午夜免费成人在线视频| 日本三级黄在线观看| 亚洲熟妇中文字幕五十中出| 亚洲成av人片在线播放无| 日日干狠狠操夜夜爽| 中文亚洲av片在线观看爽| 久久午夜亚洲精品久久| 村上凉子中文字幕在线| 午夜久久久久精精品| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 成人无遮挡网站| 国内少妇人妻偷人精品xxx网站 | 日韩免费av在线播放| 精品久久蜜臀av无| 岛国视频午夜一区免费看| 少妇的逼水好多| 国产成人精品久久二区二区91| 日本黄色视频三级网站网址| 日本撒尿小便嘘嘘汇集6| 久久精品aⅴ一区二区三区四区| 日韩大尺度精品在线看网址| 黄频高清免费视频| 精品国产亚洲在线| 欧美丝袜亚洲另类 | 亚洲 欧美一区二区三区| 香蕉丝袜av| av天堂中文字幕网| bbb黄色大片| 一a级毛片在线观看| 亚洲欧美日韩高清在线视频| 观看美女的网站| 国产精品久久久久久久电影 | netflix在线观看网站| 天天躁日日操中文字幕| 黄色 视频免费看| 久久精品国产99精品国产亚洲性色| 亚洲激情在线av| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看 | 99热只有精品国产| 香蕉久久夜色| 美女cb高潮喷水在线观看 | www.熟女人妻精品国产| 九九热线精品视视频播放| 757午夜福利合集在线观看| 香蕉国产在线看| 精品久久久久久,| xxx96com| 欧美日韩乱码在线| 国产一级毛片七仙女欲春2| 欧美丝袜亚洲另类 | 国产亚洲精品av在线| 午夜福利在线观看吧| www.自偷自拍.com| www.熟女人妻精品国产| 噜噜噜噜噜久久久久久91| 在线观看66精品国产| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 国产激情欧美一区二区| 国产精品久久久久久人妻精品电影| 国内精品美女久久久久久| 亚洲精品456在线播放app | 欧美乱色亚洲激情| 黄片大片在线免费观看| 亚洲精品一区av在线观看| 欧美大码av| 亚洲九九香蕉| 成人特级黄色片久久久久久久| 国产视频内射| 岛国视频午夜一区免费看| 国产免费男女视频| 好男人在线观看高清免费视频| 亚洲成人中文字幕在线播放| 中文字幕久久专区| av视频在线观看入口| 国产高潮美女av| 精品99又大又爽又粗少妇毛片 | 亚洲va日本ⅴa欧美va伊人久久| 一本一本综合久久| 97人妻精品一区二区三区麻豆| 怎么达到女性高潮| 欧美在线一区亚洲| 97碰自拍视频| 国产 一区 欧美 日韩| 免费大片18禁| 精品一区二区三区视频在线观看免费| 国产单亲对白刺激| 老汉色av国产亚洲站长工具| 手机成人av网站| 日韩有码中文字幕| 热99在线观看视频| 国产亚洲精品综合一区在线观看| 俺也久久电影网| 色在线成人网| 精品国内亚洲2022精品成人| 久久亚洲真实| 亚洲精品456在线播放app | 国产亚洲精品综合一区在线观看| 国产成人影院久久av| 村上凉子中文字幕在线| 人妻丰满熟妇av一区二区三区| 国产成人精品久久二区二区91| 亚洲专区中文字幕在线| av黄色大香蕉| 午夜亚洲福利在线播放| av在线蜜桃| 欧美中文综合在线视频| 亚洲av电影在线进入| 国产亚洲av高清不卡| www国产在线视频色| 国产黄a三级三级三级人| 黄色日韩在线| 日本黄色片子视频| 88av欧美| 免费无遮挡裸体视频| 在线看三级毛片| 亚洲,欧美精品.| 岛国在线观看网站| 美女 人体艺术 gogo| 在线免费观看的www视频| 欧美日韩瑟瑟在线播放| 热99re8久久精品国产| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 99精品久久久久人妻精品| 男女下面进入的视频免费午夜| 人人妻人人看人人澡| 国产精华一区二区三区| 色噜噜av男人的天堂激情| 亚洲国产欧美网| 给我免费播放毛片高清在线观看| 午夜福利欧美成人| 亚洲午夜理论影院| 精品欧美国产一区二区三| 男女下面进入的视频免费午夜| netflix在线观看网站| 真人一进一出gif抽搐免费| 国产av在哪里看| 亚洲天堂国产精品一区在线| 亚洲av日韩精品久久久久久密| 精品福利观看| 在线观看日韩欧美| 亚洲国产看品久久| 精品国产乱子伦一区二区三区| 免费大片18禁| 国产极品精品免费视频能看的| av片东京热男人的天堂|