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

    滑坡運(yùn)動(dòng)二維黏壺分析模型

    2023-10-28 09:17:12雷先順
    電力勘測(cè)設(shè)計(jì) 2023年10期
    關(guān)鍵詞:滑體分析模型模型試驗(yàn)

    雷先順,李 漪

    (1.湖北省電力勘測(cè)設(shè)計(jì)院有限公司,湖北 武漢 430040;2.合肥工業(yè)大學(xué),安徽 合肥 230009)

    0 引言

    在我國,每年發(fā)生的各種地質(zhì)災(zāi)害中,滑坡災(zāi)害發(fā)生的頻率最高,產(chǎn)生的危害最大[1]。為了對(duì)滑坡災(zāi)害進(jìn)行有效的預(yù)測(cè)和防治,需要對(duì)其發(fā)生機(jī)理、運(yùn)動(dòng)過程和災(zāi)害評(píng)估等問題進(jìn)行深入研究。

    關(guān)于滑坡的靜力學(xué)分析和災(zāi)害評(píng)估,已經(jīng)有大量的研究成果[2-3]。而滑坡運(yùn)動(dòng)過程分析是滑坡靜力學(xué)分析的后續(xù)步驟,并且關(guān)系到最終的滑坡災(zāi)害評(píng)估結(jié)果,具有重要的理論研究意義和實(shí)際用途。在滑坡運(yùn)動(dòng)過程理論研究方面,主要有質(zhì)點(diǎn)能量法、流體理論、條分法和離散元法等理論,但是各有其局限性。質(zhì)點(diǎn)能量法[4]將整個(gè)滑體當(dāng)做一個(gè)整體,運(yùn)用能量轉(zhuǎn)化關(guān)系求得最終結(jié)果,但是沒有考慮滑道的性質(zhì)與滑體各位置的變形和速度不同。流體理論[5]運(yùn)用流體力學(xué)的觀點(diǎn),將滑體看作流動(dòng)性的連續(xù)介質(zhì),但是不能反映滑坡過程中的破碎、相互作用等現(xiàn)象,且不適用于一些整體性較好的滑坡。條分法[6-9]是一種運(yùn)用較廣泛的方法,它把滑體分成若干個(gè)豎向條塊,可以和滑坡靜力學(xué)分析相結(jié)合,但是沒有考慮滑體在運(yùn)動(dòng)過程的變形和內(nèi)部耗能。離散元理論[10]基于微觀顆粒間的相互作用原理,用一個(gè)個(gè)顆粒組成的堆積體來模擬實(shí)際的滑體,可以十分方便地分析滑體內(nèi)部顆粒的力學(xué)行為,是一種較準(zhǔn)確的方法。

    此外,基于災(zāi)變鏈?zhǔn)嚼碚?,肖盛燮[11]等提出一種考慮孔隙水和空氣阻尼耦合作用的滑坡力學(xué)模型。江巍[12]等對(duì)DDA分析方法進(jìn)行了改進(jìn),可以很好地模擬滑坡運(yùn)動(dòng)過程。

    滑體在下滑過程中,是其重力勢(shì)能轉(zhuǎn)化成滑體動(dòng)能和熱能,并最終全部轉(zhuǎn)化成熱能、聲能等其他能量的過程。其中熱能包括:滑體和滑道摩擦產(chǎn)生的熱能和滑體自身變形、內(nèi)部相互摩擦碰撞產(chǎn)生的熱能。特別是對(duì)于高速巖石滑坡,在其運(yùn)動(dòng)過程中,內(nèi)部巖塊往往會(huì)發(fā)生崩解,相互碰撞作用[13]十分頻繁。所以,在對(duì)滑體運(yùn)動(dòng)過程進(jìn)行分析時(shí),有必要考慮這一部分能量損耗。本文基于條分法的特點(diǎn),考慮滑體的內(nèi)部耗能[14],采用黏壺連接代替彈簧連接,再根據(jù)條塊運(yùn)動(dòng)方程,得到滑坡最終的各項(xiàng)評(píng)估指標(biāo)。

    1 滑坡運(yùn)動(dòng)分析模型

    1.1 分析模型的提出

    在條分法中,將滑坡體離散成垂直條,且假定條塊為剛體,或者在滑動(dòng)過程中滑塊仍然保持垂直狀態(tài)。但事實(shí)上,這些垂直塊體在運(yùn)動(dòng)過程中,會(huì)變成曲線狀,并且顆粒組成越松散,滑速越大,滑體變形越普遍,甚至?xí)霈F(xiàn)碎屑流現(xiàn)象,這與傳統(tǒng)條分法的假定顯然有很大的不同。

    根據(jù)滑體在運(yùn)動(dòng)過程中能量轉(zhuǎn)化的規(guī)律,對(duì)于消耗的熱能,其中滑體和滑道摩擦消耗的能量可以通過滑塊底部摩擦力做功求得。而滑體由于自身變形,內(nèi)部顆粒間互相作用消耗的能量則無法定量求出,吳越[14]等通過試驗(yàn)發(fā)現(xiàn),在滑坡過程中,滑體內(nèi)部耗能所占比例較大,是不能忽略的,且內(nèi)部變形越大,內(nèi)部耗能所占比例越大。這里,基于滑坡過程中內(nèi)部耗能的特點(diǎn),并結(jié)合傳統(tǒng)條分法的一些優(yōu)點(diǎn),提出一種滑坡運(yùn)動(dòng)過程黏壺分析模型。

    將滑坡運(yùn)動(dòng)考慮為平面問題,即不考慮滑體的橫向擴(kuò)展問題,將滑體分成等質(zhì)量的若干個(gè)滑塊,滑塊間用黏壺連接。相鄰前后滑塊間速度不同時(shí),會(huì)導(dǎo)致黏壺變形,利用黏壺變形做功來模擬滑坡過程中的內(nèi)部能耗。

    另外,在滑體運(yùn)動(dòng)軌跡方面,也做出了相應(yīng)的簡(jiǎn)化處理,將運(yùn)動(dòng)軌跡簡(jiǎn)化成滑坡啟動(dòng)區(qū)、滑坡加速區(qū)和滑坡減速堆積區(qū)?;逻\(yùn)動(dòng)模型如圖1所示。

    圖1 滑坡運(yùn)動(dòng)模型

    由滑體剪切試驗(yàn)可知,滑體剪切段的剪應(yīng)力[15]為:

    式中:G為滑體剪切模量;γ為滑體某處剪應(yīng)變,μ為滑體內(nèi)某處剪切位移;μm為滑體剪切破壞時(shí)的臨界位移;τm為滑體的殘余剪切強(qiáng)度。

    只考慮滑體發(fā)生剪切破壞后的過程,即滑體內(nèi)部的剪切力達(dá)到其殘余剪切強(qiáng)度。另外,當(dāng)相鄰滑塊速度不同時(shí),兩者會(huì)出現(xiàn)一定的相對(duì)剪切滑動(dòng),產(chǎn)生剪切力;反之,則沒有剪切力。從而,可以得到滑體內(nèi)崩解剪切面上的剪切力:

    式中:τ為崩解面上的剪切力;V1和V2分別為相鄰兩滑塊的速度。

    而對(duì)崩解剪切面進(jìn)行分析,由庫倫—摩爾準(zhǔn)則:

    式中:σ為剪切面上的正壓力;θ為每個(gè)剪切面上的殘余內(nèi)摩擦角??梢杂涩F(xiàn)場(chǎng)勘察資料得到。

    大量滑坡實(shí)例和室內(nèi)試驗(yàn)發(fā)現(xiàn)滑體運(yùn)動(dòng)過程中,其自身發(fā)生了很大的變形,必然存在顆粒間的相互摩擦,滑體內(nèi)部會(huì)損失部分能量,這里將這部分內(nèi)部耗能轉(zhuǎn)化成滑塊間黏壺軸向變形做功耗能。

    則黏壺的黏滯力T等于崩解面上的總剪力Q,黏壺軸向變形量等于上下相鄰滑塊的相對(duì)運(yùn)動(dòng)距離。可以得到:

    將式(3)帶入式(4)中,可得:

    這里,將滑體分成n個(gè)相互平行的等質(zhì)量塊體,從下到上依次編號(hào)為1至n號(hào),假設(shè)滑體在運(yùn)動(dòng)過程中相鄰塊體間均發(fā)生了剪切破壞,即滑體內(nèi)出現(xiàn)n-1個(gè)內(nèi)部崩解剪切面,則有n-1個(gè)黏壺。并且,假設(shè)滑體中位于上層的滑塊先滑出,滑塊由上到下依次運(yùn)動(dòng),則崩解剪切面上的正壓力σ即為塊體的自身重力在垂直作用在剪切面上的分量,然后推導(dǎo)得到黏壺模型中的黏滯力為:

    巖土體的殘余強(qiáng)度大小與剪切速率快慢關(guān)系密切,式(6)中k為與剪切速率有關(guān)的修正系數(shù),針對(duì)不同的滑坡形式,可以選取不同的值。

    從而可以得到黏壺黏滯力的表達(dá)式:

    1.2 建立運(yùn)動(dòng)方程

    假設(shè)滑塊所受到的黏壺黏滯力均是平行于滑面方向,且取沿滑面向下方向?yàn)檎较?。?duì)于相鄰的兩個(gè)滑塊,若V1>V2則滑塊1所受的黏壺黏滯力為抗滑力,滑塊2所受的黏壺黏滯力則為下滑力;反之,若V1<V2,則結(jié)果相反。從而,可以建立每個(gè)滑塊的運(yùn)動(dòng)學(xué)方程,沿滑面方向,滑塊所受的合力為:

    則有:

    在具體計(jì)算時(shí),可以假設(shè)所有滑塊的初速度均為零,而且滑坡開始運(yùn)動(dòng)時(shí),所有黏壺均處于拉伸狀態(tài)。那么,可以得出初始時(shí),滑塊的加速度。然后,再取一個(gè)很小的時(shí)間間隔Δt,可認(rèn)為在(t,t+Δt)時(shí)段內(nèi)滑塊的加速度不變,即滑塊做勻加速直線運(yùn)動(dòng),則變形塊在t+Δt時(shí)刻的速度為:

    在t+Δt時(shí)刻的位置為:

    當(dāng)1號(hào)滑塊的位移達(dá)到最大值時(shí),停止計(jì)算,即可以得到最終滑坡的沖程和運(yùn)動(dòng)時(shí)間t總。相應(yīng)地,也可以得到滑坡過程中的內(nèi)部耗能和滑體自身的變形量。其中,內(nèi)部耗能可以通過黏壺軸向力做功求出,第i個(gè)黏壺耗能:

    則整個(gè)滑體的內(nèi)部耗能為:

    滑體自身的變形量即為1號(hào)滑塊和n號(hào)滑塊運(yùn)動(dòng)距離的差值:

    1.3 計(jì)算步驟

    通過C語言編程,可以非常方便地進(jìn)行計(jì)算。且不需要往復(fù)迭代,最終可以得到滑坡的沖程、任意時(shí)刻的滑體變形量、最大速度、耗能情況等參數(shù)。這里,將第一塊滑塊速度為0時(shí)作為程序終止的判定條件。具體計(jì)算步驟如圖2所示。

    圖2 計(jì)算步驟流程圖

    2 計(jì)算與分析

    根據(jù)前文,該滑坡分析模型的具體推導(dǎo)過程,現(xiàn)運(yùn)用該黏壺分析模型對(duì)前期的室內(nèi)模型試驗(yàn)進(jìn)行計(jì)算分析,并和離散元軟件(PFC 2D)的計(jì)算結(jié)果進(jìn)行對(duì)比。

    2.1 模型計(jì)算結(jié)果

    模型試驗(yàn)裝置如圖3所示。其中,滑坡啟動(dòng)區(qū)、滑坡加速區(qū)和減速堆積區(qū)的坡度分別為30°、30°和0°,三段長(zhǎng)度分別為0.4 m、3.6 m和10 m。試驗(yàn)所用顆粒材料是碎石子,其粒徑范圍為0~20 mm,堆積體重度為13.5 kN/m3。室內(nèi)試驗(yàn)的結(jié)果滑體的最終沖程為0.92 m,堆積長(zhǎng)度為1.1 m,總共耗時(shí)4.5 s。在其運(yùn)動(dòng)過程中,出現(xiàn)了明顯的顆粒跳動(dòng)、相互碰撞等現(xiàn)象,滑體自身發(fā)生了很大的變形。如果用常規(guī)的條分法,顯然是不符合實(shí)際的。

    圖3 室內(nèi)模型試驗(yàn)圖

    首先,運(yùn)用本文提出的分析模型進(jìn)行計(jì)算。因?yàn)槟P驮囼?yàn)中滑體的運(yùn)動(dòng)路徑均為直線段,在分析模型中,也只考慮3個(gè)直線段。在該計(jì)算模型中,參數(shù)的取值(包括各位置殘余內(nèi)摩擦角θ、修正系數(shù)k和滑面的摩擦系數(shù))對(duì)最終結(jié)果的影響較大?;w和滑面的參數(shù)取值需要經(jīng)過試驗(yàn)和反分析相結(jié)合的方法得到?;w的體積為40 L,滑體的質(zhì)量為55.1 kg。這里,需要確定的系數(shù)包括滑體各位置殘余內(nèi)摩擦角、修正系數(shù)k和滑面摩擦系數(shù)。由于滑體的組成比較均勻,所以可以認(rèn)為滑體各位置處顆粒的殘余內(nèi)摩擦角θ相同,即在確定黏壺黏滯力時(shí),可以將每個(gè)崩解面上的tanθ取相同的數(shù)值。三段運(yùn)動(dòng)路徑中,啟動(dòng)區(qū)和加速區(qū)的滑面材料相同,均為光滑板,而減速區(qū)滑面采用粗糙板。三段路徑的長(zhǎng)度和坡度均和模型試驗(yàn)保持一致,這里三段路徑的長(zhǎng)度S1=0.4 m,S2=3.6 m,S3=10 m;三段路徑坡度分別為α1=α2=30°,α3=0°;三段路徑的摩擦系數(shù)分別為0.51、0.36和0.62。碎石子的直剪摩擦角為40.1°。由于該模型試驗(yàn)中滑體全部由無黏性顆粒組成,其運(yùn)動(dòng)形式類似于碎屑流,內(nèi)部崩解速率相對(duì)較快。綜合參考文獻(xiàn)[14]中的取值,這里將修正系數(shù)K設(shè)定為0.15。程序計(jì)算中,將計(jì)算時(shí)步取為0.05 s。

    這里先探究滑塊不同分塊數(shù)n的影響,如圖4所示。從圖4中可以看出,當(dāng)滑塊數(shù)n取不同的數(shù)值時(shí),對(duì)1號(hào)滑塊最終的運(yùn)動(dòng)距離S有一定的影響。當(dāng)n=2時(shí),S為5.23 m;當(dāng)n=7時(shí),S為5.33 m;當(dāng)n=12時(shí),S為5.27 m;當(dāng)n繼續(xù)增大時(shí),S最終結(jié)果保持在5.22左右。即當(dāng)n取不同值時(shí),對(duì)最終結(jié)果的影響在10%之內(nèi)。實(shí)際計(jì)算時(shí),可以將n取一系列的數(shù)值,然后再找出計(jì)算結(jié)果的范圍。另外,可以看出滑體的運(yùn)動(dòng)距離都是隨時(shí)間的變化呈先增加后逐漸趨于穩(wěn)定。當(dāng)運(yùn)動(dòng)到2.7 s時(shí),1號(hào)滑塊的運(yùn)動(dòng)距離幾乎保持不變,而后面的滑塊繼續(xù)運(yùn)動(dòng),和模型試驗(yàn)的試驗(yàn)現(xiàn)象比較接近。

    圖4 滑塊分塊數(shù)對(duì)1號(hào)滑塊運(yùn)動(dòng)距離的影響規(guī)律

    因此,以n=7為例,則模型中黏壺的黏滯力為:

    T=σ×tanθ×k= 9.7 N

    然后,將各參數(shù)帶入編制好的程序中,經(jīng)過計(jì)算可以得到,7個(gè)滑塊的最終運(yùn)動(dòng)距離依次為:5.33 m、5.05 m、4.90 m、4.67 m、4.42 m、4.14 m和3.84 m。即滑體最終的沖程是1.33 m,相對(duì)于模型試驗(yàn)的結(jié)果,誤差為44.5%?;w的堆積長(zhǎng)度為1.49 m,相對(duì)于模型試驗(yàn)的結(jié)果,誤差為35.4%;滑體的內(nèi)部變形為1.09 m,相對(duì)于初始堆積長(zhǎng)度,變形比例達(dá)到172%。表1為計(jì)算結(jié)果和模型試驗(yàn)結(jié)果對(duì)比情況。另外,分析滑體堆積長(zhǎng)度變化情況,如圖5所示。

    表1 理論計(jì)算與模型試驗(yàn)結(jié)果對(duì)比 m

    圖5 滑體長(zhǎng)度變化歷時(shí)圖

    可以看出,開始階段滑體前端的滑塊先運(yùn)動(dòng),而后面滑塊運(yùn)動(dòng)較慢,導(dǎo)致前后滑塊的距離越來越遠(yuǎn)。當(dāng)運(yùn)動(dòng)到1.8 s左右時(shí),后方滑塊速度大于前端滑塊,前后滑塊的相對(duì)距離逐漸減小,直至到達(dá)穩(wěn)定值。這和試驗(yàn)中觀察到的現(xiàn)象相同,即滑體先逐漸在斜面上展開,然后再逐漸在減速區(qū)上堆積,該模型可以較好地模擬滑體的內(nèi)部變形特性。

    2.2 計(jì)算結(jié)果對(duì)比分析

    離散元是一種比較好的分析方法,可以模擬物體發(fā)生的變形,并且可以分析其運(yùn)動(dòng)過程。PFC2D的數(shù)值模型和室內(nèi)試驗(yàn)的模型大致相同,如圖6所示。主要由裝料槽、滑動(dòng)加速面和減速堆積面組成,且尺寸和模型試驗(yàn)尺寸相同,滑道強(qiáng)度摩擦角為31.8°。堆積體顆??偣舶?00個(gè)球,半徑在0~2 cm范圍內(nèi),并按照實(shí)際砂子的顆粒級(jí)配分布。球間摩擦系數(shù)為0.15,球間法向和切向剛度均為1e8,球的密度為2.49×103kg/m3。通過計(jì)算得到滑體的最終沖程為1.10 m,和試驗(yàn)結(jié)果0.92 m的差距較小。

    圖6 PFC模型

    圖7為兩種計(jì)算方法下,模型試驗(yàn)中滑坡體的最終堆積狀態(tài)。本文計(jì)算的結(jié)果和試驗(yàn)現(xiàn)象較為吻合,且計(jì)算模型方法較簡(jiǎn)單,易于編程,但是由于沒有考慮滑體的形狀特征,不能模擬出滑體運(yùn)動(dòng)過程中滑體各位置處高度的變化,從而無法得到滑體最終堆積形態(tài),這是需要改進(jìn)的地方。而運(yùn)用PFC分析軟件可以較好地模擬室內(nèi)模型試驗(yàn)的結(jié)果,最終的堆積形態(tài)十分相似,但由于在離散元方法中,各試驗(yàn)參數(shù)取值較困難,需要經(jīng)過試算,方法不易掌握。

    圖7 計(jì)算結(jié)果對(duì)比圖

    而后,分析滑體運(yùn)動(dòng)過程中各參數(shù)值的變化情況。在滑坡評(píng)估中,比較重要的是滑坡最終的沖程,圖8為兩種計(jì)算方法下,1號(hào)滑塊運(yùn)動(dòng)距離隨時(shí)間變化圖??芍?,2種方法中滑體運(yùn)動(dòng)距離隨時(shí)間的變化規(guī)律大致相同,滑體的運(yùn)動(dòng)距離都是隨時(shí)間的變化呈先增加后逐漸趨于穩(wěn)定。本文計(jì)算模型得出的滑坡運(yùn)動(dòng)距離相對(duì)于PFC 2D計(jì)算結(jié)果較大,相對(duì)誤差為20%。另一方面,2種方法計(jì)算的滑坡運(yùn)動(dòng)時(shí)間不同,本文計(jì)算的運(yùn)動(dòng)時(shí)間為2.7 s,而運(yùn)用PFC 2D計(jì)算的運(yùn)動(dòng)時(shí)間為3.9 s,實(shí)際運(yùn)動(dòng)時(shí)間為4.5 s??偟膩碚f,本文計(jì)算模型的計(jì)算結(jié)果均使最終的風(fēng)險(xiǎn)偏大,且誤差相對(duì)較小。

    圖8 1號(hào)滑塊運(yùn)動(dòng)距離歷時(shí)圖

    3 結(jié)論

    本文基于滑坡運(yùn)動(dòng)過程中能量轉(zhuǎn)化規(guī)律,考慮滑體運(yùn)動(dòng)過程中自身崩解、碰撞等內(nèi)部耗能,提出一種滑坡運(yùn)動(dòng)分析模型,并編制了相關(guān)計(jì)算程序。

    1)采用本文分析模型,可以簡(jiǎn)單、準(zhǔn)確地計(jì)算滑坡過程中的最大滑速、滑坡沖程、運(yùn)動(dòng)時(shí)間等滑坡致災(zāi)指標(biāo)。

    2)基于室內(nèi)模型試驗(yàn),分別利用該分析模型和離散元分析進(jìn)行計(jì)算,并對(duì)比計(jì)算結(jié)果,證明了該分析模型的適用性。

    3)本文計(jì)算模型可以較好地模擬出現(xiàn)變形的滑坡現(xiàn)象,相比傳統(tǒng)剛體運(yùn)動(dòng)模型,更加符合實(shí)際情況。

    4)由于沒有考慮滑體的結(jié)構(gòu)組成,滑道定義過于簡(jiǎn)化等原因,本文提出的滑坡分析模型存在一定的不足。如不能模擬出滑坡過程中的滑體高度變化情況、不能分析滑體起始破壞情況等。

    猜你喜歡
    滑體分析模型模型試驗(yàn)
    基于BERT-VGG16的多模態(tài)情感分析模型
    反推力裝置模型試驗(yàn)臺(tái)的研制及驗(yàn)證
    滑坡碎屑流顆粒分選效應(yīng)的數(shù)值模擬
    立式旋壓機(jī)縱向進(jìn)給機(jī)構(gòu)液壓配重設(shè)計(jì)
    萬梁高速某滑坡降雨入滲穩(wěn)定性及處治技術(shù)研究*
    露天礦反鏟挖掘機(jī)處理滑體的方式
    層次分析模型在結(jié)核疾病預(yù)防控制系統(tǒng)中的應(yīng)用
    臺(tái)階式短加筋土擋墻行為特征的離心模型試驗(yàn)
    全啟發(fā)式語言分析模型
    巨厚堅(jiān)硬巖漿巖不同配比的模型試驗(yàn)研究
    搞女人的毛片| 网址你懂的国产日韩在线| 亚洲第一区二区三区不卡| 国产 一区精品| 久久人人爽人人爽人人片va| 午夜免费激情av| 国产单亲对白刺激| 欧美人与善性xxx| 久久久久久九九精品二区国产| 日韩欧美三级三区| 国产精华一区二区三区| 日本在线视频免费播放| 国内精品一区二区在线观看| 免费看a级黄色片| 九色成人免费人妻av| 97热精品久久久久久| 免费看日本二区| 赤兔流量卡办理| 午夜福利在线观看免费完整高清在 | 国内精品宾馆在线| 天堂网av新在线| 久久久久国产精品人妻aⅴ院| 久久久午夜欧美精品| 十八禁网站免费在线| 免费搜索国产男女视频| 久久精品国产亚洲av天美| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 久久久久久大精品| 51国产日韩欧美| 国产色婷婷99| 麻豆久久精品国产亚洲av| 他把我摸到了高潮在线观看| 在线观看午夜福利视频| 亚洲天堂国产精品一区在线| 永久网站在线| 亚洲熟妇熟女久久| 麻豆国产av国片精品| 九色国产91popny在线| 欧美日韩综合久久久久久 | 51国产日韩欧美| 亚洲在线观看片| 欧美区成人在线视频| 久久99热这里只有精品18| 久久6这里有精品| 国产亚洲欧美98| 午夜日韩欧美国产| 人人妻人人澡欧美一区二区| 噜噜噜噜噜久久久久久91| 成人性生交大片免费视频hd| 看片在线看免费视频| 国产一级毛片七仙女欲春2| 国产伦精品一区二区三区视频9| 最近中文字幕高清免费大全6 | 日本黄大片高清| 久久久久久久久大av| 日韩av在线大香蕉| 免费黄网站久久成人精品| 看免费成人av毛片| 国产主播在线观看一区二区| 久久久久国内视频| 国产成人a区在线观看| 日韩中文字幕欧美一区二区| 国产男人的电影天堂91| 乱系列少妇在线播放| 国产色婷婷99| 波多野结衣高清无吗| 色综合站精品国产| 最新在线观看一区二区三区| 看免费成人av毛片| 国产91精品成人一区二区三区| 国产精品久久视频播放| 亚洲av成人精品一区久久| 九九在线视频观看精品| 精品乱码久久久久久99久播| 毛片女人毛片| 国产在线精品亚洲第一网站| 嫩草影视91久久| 成人特级av手机在线观看| 最后的刺客免费高清国语| 狠狠狠狠99中文字幕| 国产真实伦视频高清在线观看 | 在线免费观看不下载黄p国产 | 很黄的视频免费| 三级男女做爰猛烈吃奶摸视频| 久久久久国产精品人妻aⅴ院| 亚洲av中文字字幕乱码综合| 亚洲av不卡在线观看| 久久久国产成人免费| 干丝袜人妻中文字幕| 高清在线国产一区| 老师上课跳d突然被开到最大视频| 伦精品一区二区三区| 亚洲国产精品合色在线| ponron亚洲| 日本 av在线| 亚洲内射少妇av| 国模一区二区三区四区视频| 熟女人妻精品中文字幕| 欧美色欧美亚洲另类二区| 欧美成人免费av一区二区三区| 国产在线精品亚洲第一网站| 一级黄色大片毛片| 国产精品福利在线免费观看| 99久久久亚洲精品蜜臀av| 亚洲精品成人久久久久久| 欧美日韩精品成人综合77777| 亚洲性夜色夜夜综合| 中亚洲国语对白在线视频| 久久精品国产99精品国产亚洲性色| 国产精品电影一区二区三区| 国语自产精品视频在线第100页| 日本五十路高清| a级毛片免费高清观看在线播放| 18禁裸乳无遮挡免费网站照片| 一进一出抽搐动态| 日本三级黄在线观看| 亚洲精华国产精华液的使用体验 | 久9热在线精品视频| 舔av片在线| 国产亚洲av嫩草精品影院| 欧美色欧美亚洲另类二区| 国产成人a区在线观看| 免费看光身美女| 国产 一区精品| 最近中文字幕高清免费大全6 | 黄色配什么色好看| 悠悠久久av| 国内少妇人妻偷人精品xxx网站| 久久人人爽人人爽人人片va| 麻豆国产av国片精品| 两人在一起打扑克的视频| 91麻豆av在线| 日日摸夜夜添夜夜添小说| 日韩国内少妇激情av| 白带黄色成豆腐渣| 亚洲精品成人久久久久久| 午夜福利高清视频| 国产精品无大码| 亚洲国产欧美人成| 国产精品免费一区二区三区在线| 男女做爰动态图高潮gif福利片| 亚洲一级一片aⅴ在线观看| 精品久久久久久久末码| 九九爱精品视频在线观看| 真实男女啪啪啪动态图| 网址你懂的国产日韩在线| 亚洲自拍偷在线| 美女xxoo啪啪120秒动态图| 在线免费观看不下载黄p国产 | 精品一区二区免费观看| 成人午夜高清在线视频| 亚洲成a人片在线一区二区| 欧美一级a爱片免费观看看| 国产高清不卡午夜福利| 一区福利在线观看| 久久久久久国产a免费观看| 亚洲五月天丁香| 亚洲图色成人| 亚洲国产精品合色在线| ponron亚洲| 色吧在线观看| 国产精品久久视频播放| 亚洲av不卡在线观看| 在线国产一区二区在线| 黄色配什么色好看| 国产免费男女视频| 国产伦在线观看视频一区| 亚洲一区二区三区色噜噜| 18禁黄网站禁片免费观看直播| 国产精品免费一区二区三区在线| 校园人妻丝袜中文字幕| 熟女人妻精品中文字幕| 最近最新中文字幕大全电影3| 日本熟妇午夜| 2021天堂中文幕一二区在线观| 他把我摸到了高潮在线观看| 亚洲性夜色夜夜综合| 美女cb高潮喷水在线观看| 日本三级黄在线观看| 国产免费男女视频| 久久久久久久久久黄片| 久久草成人影院| 亚洲人成网站在线播放欧美日韩| 九九热线精品视视频播放| 美女高潮的动态| 亚洲经典国产精华液单| 丰满乱子伦码专区| 色尼玛亚洲综合影院| 国产高清视频在线播放一区| 亚洲五月天丁香| 成人一区二区视频在线观看| 在线观看一区二区三区| 日韩 亚洲 欧美在线| 久久精品国产亚洲av涩爱 | 亚洲专区中文字幕在线| 日本爱情动作片www.在线观看 | 在线观看免费视频日本深夜| 我要搜黄色片| 免费人成在线观看视频色| 国产真实伦视频高清在线观看 | 国产精品一区二区性色av| 亚洲最大成人av| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 亚洲无线在线观看| 又爽又黄a免费视频| 欧美一区二区精品小视频在线| 欧美日韩乱码在线| 国产色爽女视频免费观看| 久久久午夜欧美精品| 男女视频在线观看网站免费| 在线播放无遮挡| 一区福利在线观看| 国产精品美女特级片免费视频播放器| 国产白丝娇喘喷水9色精品| 国产精品一区二区三区四区久久| 亚洲久久久久久中文字幕| 日本黄大片高清| 国产黄色小视频在线观看| 久久国产精品人妻蜜桃| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 一个人看视频在线观看www免费| 国产精品免费一区二区三区在线| 亚洲欧美精品综合久久99| 我的老师免费观看完整版| 亚洲av一区综合| 精品人妻一区二区三区麻豆 | 久久国产乱子免费精品| 亚洲国产色片| 黄色日韩在线| 美女xxoo啪啪120秒动态图| 日韩欧美国产在线观看| 午夜a级毛片| 可以在线观看的亚洲视频| 国产精品1区2区在线观看.| 哪里可以看免费的av片| 韩国av在线不卡| 精品国内亚洲2022精品成人| 色综合色国产| 国产精品一及| 午夜精品在线福利| 中文字幕久久专区| 变态另类成人亚洲欧美熟女| 日韩精品有码人妻一区| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 国产精品电影一区二区三区| 日韩一本色道免费dvd| 天堂动漫精品| 一本一本综合久久| 国内精品久久久久精免费| 99久久久亚洲精品蜜臀av| 少妇的逼好多水| 午夜爱爱视频在线播放| 啦啦啦韩国在线观看视频| 日韩人妻高清精品专区| 女生性感内裤真人,穿戴方法视频| 国产真实乱freesex| av专区在线播放| 神马国产精品三级电影在线观看| 99riav亚洲国产免费| 韩国av一区二区三区四区| 97超视频在线观看视频| 中文字幕免费在线视频6| 一区二区三区高清视频在线| 久久久久国内视频| 亚洲va日本ⅴa欧美va伊人久久| 九色成人免费人妻av| 亚洲美女搞黄在线观看 | 欧美精品国产亚洲| 久久精品国产亚洲av天美| 人人妻人人澡欧美一区二区| 欧美xxxx性猛交bbbb| 欧美最新免费一区二区三区| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 欧美黑人欧美精品刺激| 别揉我奶头 嗯啊视频| 国产三级在线视频| 热99re8久久精品国产| 成年女人看的毛片在线观看| 在线观看午夜福利视频| 白带黄色成豆腐渣| 日日干狠狠操夜夜爽| 看片在线看免费视频| 久久久久久久久大av| 成年女人毛片免费观看观看9| 两个人的视频大全免费| 18禁裸乳无遮挡免费网站照片| 看免费成人av毛片| 制服丝袜大香蕉在线| 99热精品在线国产| 日本精品一区二区三区蜜桃| 国产精品永久免费网站| 五月伊人婷婷丁香| 97超级碰碰碰精品色视频在线观看| 国产极品精品免费视频能看的| 亚洲av熟女| 18+在线观看网站| 国产成年人精品一区二区| 精品99又大又爽又粗少妇毛片 | 亚洲久久久久久中文字幕| 91午夜精品亚洲一区二区三区 | 国产成人一区二区在线| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站在线播| av在线观看视频网站免费| 午夜福利在线观看吧| 乱码一卡2卡4卡精品| 精品久久久久久,| 成人无遮挡网站| 国产综合懂色| 国产真实伦视频高清在线观看 | 在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| av天堂在线播放| 真人一进一出gif抽搐免费| av福利片在线观看| 女人十人毛片免费观看3o分钟| 久久久久久久久中文| 五月玫瑰六月丁香| 又黄又爽又刺激的免费视频.| 国产 一区精品| 午夜激情福利司机影院| 国产极品精品免费视频能看的| 亚洲精品国产成人久久av| 99精品久久久久人妻精品| 嫩草影院入口| 一本久久中文字幕| 国产日本99.免费观看| 特大巨黑吊av在线直播| 变态另类成人亚洲欧美熟女| 亚洲av成人av| 免费在线观看成人毛片| 搞女人的毛片| 亚洲av.av天堂| 可以在线观看的亚洲视频| 久久精品国产亚洲av涩爱 | 中文字幕av在线有码专区| 黄片wwwwww| 性色avwww在线观看| 可以在线观看毛片的网站| 午夜影院日韩av| 日韩欧美精品v在线| 国国产精品蜜臀av免费| 老女人水多毛片| 搡女人真爽免费视频火全软件 | 精品久久久久久久久久久久久| 国产 一区精品| 亚洲欧美激情综合另类| 99国产极品粉嫩在线观看| 成人av在线播放网站| 欧美三级亚洲精品| 久久久久国内视频| 日本成人三级电影网站| 日韩高清综合在线| 网址你懂的国产日韩在线| 日日干狠狠操夜夜爽| 国产伦在线观看视频一区| or卡值多少钱| 午夜日韩欧美国产| 亚洲精品国产成人久久av| 国产欧美日韩精品一区二区| 亚洲av免费在线观看| 色尼玛亚洲综合影院| 少妇丰满av| 婷婷丁香在线五月| 亚洲自拍偷在线| 日本撒尿小便嘘嘘汇集6| 又粗又爽又猛毛片免费看| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线观看免费完整高清在 | 亚洲美女黄片视频| 99热这里只有是精品在线观看| 亚洲精品亚洲一区二区| 99热6这里只有精品| 久久精品人妻少妇| 三级毛片av免费| 国产高潮美女av| 真人一进一出gif抽搐免费| 日本-黄色视频高清免费观看| 亚洲va在线va天堂va国产| 国产老妇女一区| 色精品久久人妻99蜜桃| 欧美日韩综合久久久久久 | 国内精品一区二区在线观看| 丝袜美腿在线中文| 国产精品1区2区在线观看.| 少妇猛男粗大的猛烈进出视频 | 久久久久免费精品人妻一区二区| 免费观看人在逋| 亚洲性久久影院| 舔av片在线| 欧美日本视频| 国内毛片毛片毛片毛片毛片| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩卡通动漫| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6 | 免费高清视频大片| 国产伦精品一区二区三区四那| 女生性感内裤真人,穿戴方法视频| 国产精品电影一区二区三区| 精品久久久久久久久av| 日韩一区二区视频免费看| 女生性感内裤真人,穿戴方法视频| av专区在线播放| 男女边吃奶边做爰视频| 成人精品一区二区免费| 亚洲性夜色夜夜综合| 国产免费男女视频| 亚洲精品影视一区二区三区av| 欧美不卡视频在线免费观看| 成人无遮挡网站| 日本-黄色视频高清免费观看| 亚洲va在线va天堂va国产| 午夜福利在线观看免费完整高清在 | 亚洲国产欧洲综合997久久,| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 亚洲色图av天堂| 最后的刺客免费高清国语| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 亚洲专区中文字幕在线| 成人二区视频| 精品一区二区免费观看| 免费在线观看日本一区| 国产一区二区亚洲精品在线观看| 在线国产一区二区在线| 两个人视频免费观看高清| 黄色日韩在线| 欧美色视频一区免费| 男女啪啪激烈高潮av片| 老司机午夜福利在线观看视频| 国内揄拍国产精品人妻在线| 亚洲在线自拍视频| 国产一区二区三区在线臀色熟女| 人妻制服诱惑在线中文字幕| 午夜免费激情av| 国产高潮美女av| 在线看三级毛片| 精品久久久久久久久亚洲 | a级毛片免费高清观看在线播放| 国产伦精品一区二区三区四那| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 一区二区三区高清视频在线| 欧美3d第一页| 日韩欧美一区二区三区在线观看| 欧美最黄视频在线播放免费| 欧美成人免费av一区二区三区| 国产午夜精品久久久久久一区二区三区 | 给我免费播放毛片高清在线观看| 1000部很黄的大片| 亚洲一区高清亚洲精品| 亚洲综合色惰| 啦啦啦韩国在线观看视频| 我要看日韩黄色一级片| 伦理电影大哥的女人| 午夜老司机福利剧场| 国内精品美女久久久久久| 免费在线观看成人毛片| 免费看a级黄色片| 少妇人妻精品综合一区二区 | 熟妇人妻久久中文字幕3abv| 精品久久久久久久末码| 久久热精品热| 成熟少妇高潮喷水视频| 亚洲av美国av| 日日摸夜夜添夜夜添av毛片 | 色哟哟哟哟哟哟| 中文字幕精品亚洲无线码一区| 亚洲自拍偷在线| 亚洲第一电影网av| 国产一区二区亚洲精品在线观看| 国产探花极品一区二区| 天天躁日日操中文字幕| 国产爱豆传媒在线观看| 成熟少妇高潮喷水视频| 国产欧美日韩精品亚洲av| 日韩一区二区视频免费看| 亚洲最大成人av| 久久久久国产精品人妻aⅴ院| 精品免费久久久久久久清纯| 国产精品乱码一区二三区的特点| 高清毛片免费观看视频网站| 国产极品精品免费视频能看的| 婷婷六月久久综合丁香| 亚洲av不卡在线观看| 亚洲av美国av| 亚洲性久久影院| 亚洲精品久久国产高清桃花| 中文字幕av在线有码专区| 欧美在线一区亚洲| 色哟哟哟哟哟哟| 干丝袜人妻中文字幕| av在线观看视频网站免费| 直男gayav资源| 国产一区二区激情短视频| 麻豆成人午夜福利视频| 22中文网久久字幕| 精品午夜福利视频在线观看一区| 日本精品一区二区三区蜜桃| 真实男女啪啪啪动态图| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 好男人在线观看高清免费视频| 欧美最新免费一区二区三区| 18禁黄网站禁片午夜丰满| 一a级毛片在线观看| 无人区码免费观看不卡| 亚洲av中文av极速乱 | 久久婷婷人人爽人人干人人爱| 能在线免费观看的黄片| av国产免费在线观看| 亚洲久久久久久中文字幕| 久久人人爽人人爽人人片va| 亚洲va日本ⅴa欧美va伊人久久| 99久久精品国产国产毛片| 亚洲精品在线观看二区| 久久久国产成人精品二区| 97超视频在线观看视频| 欧美激情国产日韩精品一区| 给我免费播放毛片高清在线观看| 在线观看免费视频日本深夜| 中文字幕精品亚洲无线码一区| 男女做爰动态图高潮gif福利片| 国产日本99.免费观看| 99国产极品粉嫩在线观看| 高清在线国产一区| 亚洲狠狠婷婷综合久久图片| 欧美日本视频| 国产精华一区二区三区| 亚洲欧美日韩无卡精品| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 天天一区二区日本电影三级| 99在线视频只有这里精品首页| 免费av不卡在线播放| 在线免费观看不下载黄p国产 | 国产精品伦人一区二区| 国产一区二区三区av在线 | 国产高清有码在线观看视频| 国产一区二区三区视频了| 女人十人毛片免费观看3o分钟| 国产精品综合久久久久久久免费| 欧美成人免费av一区二区三区| 国产av不卡久久| 免费一级毛片在线播放高清视频| 性插视频无遮挡在线免费观看| 嫩草影视91久久| 亚洲av免费高清在线观看| 精华霜和精华液先用哪个| 欧美潮喷喷水| 色播亚洲综合网| 搞女人的毛片| 一个人观看的视频www高清免费观看| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 最近视频中文字幕2019在线8| 国产精品野战在线观看| 男女视频在线观看网站免费| 亚洲精品在线观看二区| 国产高清三级在线| 丰满的人妻完整版| 天堂影院成人在线观看| 欧美日韩乱码在线| 91久久精品电影网| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品久久男人天堂| 丰满乱子伦码专区| 国产亚洲精品久久久久久毛片| 丝袜美腿在线中文| 婷婷精品国产亚洲av| 午夜激情福利司机影院| 在线免费观看的www视频| 伊人久久精品亚洲午夜| 男女视频在线观看网站免费| 国产av在哪里看| 精品久久久久久久久久免费视频| 久久6这里有精品| 国产免费一级a男人的天堂| 免费在线观看成人毛片| 五月玫瑰六月丁香| 国产白丝娇喘喷水9色精品| 久久久久性生活片| 搡女人真爽免费视频火全软件 | 中文字幕高清在线视频| 久久久久国内视频| 久久精品91蜜桃| 校园春色视频在线观看| 丰满乱子伦码专区| 亚洲电影在线观看av| 很黄的视频免费| eeuss影院久久| 国产真实乱freesex| 久久亚洲精品不卡| 久久精品国产鲁丝片午夜精品 | 亚洲中文日韩欧美视频| 日韩欧美精品免费久久| 美女免费视频网站| 国产日本99.免费观看| 免费高清视频大片| 亚洲专区国产一区二区| 久久精品国产99精品国产亚洲性色| 精品久久国产蜜桃| 成人精品一区二区免费| 国产亚洲av嫩草精品影院| 成人国产麻豆网| 国产精品亚洲美女久久久| 精品人妻熟女av久视频|