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

    波致海底緩傾角無限坡滑動(dòng)穩(wěn)定性計(jì)算分析探討

    2015-06-01 12:30:12劉小麗竇錦鐘英姿霍沿東張旭
    海洋學(xué)報(bào) 2015年3期
    關(guān)鍵詞:坡角坡體波浪

    劉小麗,竇錦鐘,英姿,霍沿東,張旭

    (1.中國海洋大學(xué)環(huán)境科學(xué)與工程學(xué)院,山東青島 266100;2.海洋環(huán)境與生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,山東青島 266100)

    波致海底緩傾角無限坡滑動(dòng)穩(wěn)定性計(jì)算分析探討

    劉小麗1,2,竇錦鐘1,英姿1,霍沿東1,張旭1

    (1.中國海洋大學(xué)環(huán)境科學(xué)與工程學(xué)院,山東青島 266100;2.海洋環(huán)境與生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,山東青島 266100)

    波浪作用下海底無限坡滑動(dòng)穩(wěn)定性計(jì)算的極限平衡法中,忽略了坡體水平向應(yīng)力狀態(tài)的影響,為此,針對波浪作用下海底緩傾角無限邊坡的特點(diǎn),提出直接基于滑動(dòng)面處土體應(yīng)力狀態(tài)的滑動(dòng)穩(wěn)定性計(jì)算方法(簡稱應(yīng)力狀態(tài)法),并分析了其適用范圍。對具體算例的分析表明,應(yīng)力狀態(tài)法計(jì)算得出的安全系數(shù)大于極限平衡法的安全系數(shù),且隨著滑動(dòng)面深度的增加、土體泊松比以及邊坡坡角的增大,兩種計(jì)算方法得出的安全系數(shù)的差異會(huì)逐漸增大;對于波浪作用下的海底緩傾角無限邊坡,在失穩(wěn)時(shí)極可能沿具有一定厚度的滑動(dòng)帶而不是單一的滑動(dòng)面而滑動(dòng),且波致最大剪應(yīng)力所在的深度,常常不是斜坡體最易失穩(wěn)滑移的深度。

    海底無限邊坡;波浪;穩(wěn)定性分析;緩傾角

    1 引言

    與陸地上的滑坡相比,海底滑坡的最顯著特征之一,是大部分的海底滑坡發(fā)生在非常平緩的斜坡上[1]。如位于挪威海域的Storegga滑坡[2]是迄今發(fā)現(xiàn)的最大的海底復(fù)合滑坡體,其滑坡區(qū)西南部坡角1.17°~1.32°,東部坡角0.55°~1.14°;南海北部陸坡區(qū)有較多滑坡發(fā)育的跡象,其東部的臺(tái)灣淺灘陸坡段,總體坡度約為3.4°,珠江海谷段陸坡的下陸坡相對陡峭處的坡度為3°左右[3];楊作升等[4]對黃河口水下滑坡體系的特征及其相互作用進(jìn)行了分析,指出廢棄的1964-1976年的黃河口沙嘴水下底坡的中上部坡度相對較陡,但其平均坡度不超過0.5°。

    除重力作用外,沉積物快速沉積導(dǎo)致的超高孔隙水壓力、淺層氣的存在、地震以及波浪作用等因素都有可能誘發(fā)海底緩傾角滑坡,如甘華陽等[5]介紹了海底天然氣水合物分解導(dǎo)致滑坡的機(jī)制及該種海底滑坡體的特征;張亮和欒錫武[3]對南海北部陸坡的穩(wěn)定性進(jìn)行了詳細(xì)定量分析研究,認(rèn)為地震是引發(fā)南海北部滑坡最主要的觸發(fā)機(jī)制之一。波浪是近岸及淺水區(qū)海底斜坡穩(wěn)定性分析中需要考慮的一個(gè)重要因素,Henkel[6]通過計(jì)算分析認(rèn)為在密西西比河口波浪能引起120 m深度處海底軟土斜坡的破壞。波浪一般所影響的水深相當(dāng)于波長的一半,當(dāng)水深小于1/2波長時(shí),波長變短、波高變大,并最終出現(xiàn)翻卷形成破波[7]。

    Henkel[6]采用圓弧形滑動(dòng)面通過力矩平衡分析了波浪對海底斜坡穩(wěn)定性的影響,Wright和Dunham[8]利用有限元分析了無限坡滑動(dòng)模式下波浪對海底斜坡的作用。對于海底緩傾角斜坡,大部分文獻(xiàn)將其視為無限坡模式利用極限平衡法進(jìn)行了穩(wěn)定性分析[2—3,9—12]。無限坡滑動(dòng)模式是一種較簡單的邊坡穩(wěn)定性分析模型,通常當(dāng)斜坡滑動(dòng)面為平行于坡面的平面(沿滑動(dòng)方向的傾角為常數(shù)),且滑坡體的厚度遠(yuǎn)小于其滑動(dòng)方向的長度(一般滑坡體的長度與厚度之比大于10)時(shí),可將其作為無限坡進(jìn)行分析[12—13]。

    在波浪導(dǎo)致海底無限坡滑動(dòng)穩(wěn)定性計(jì)算中,常采用極限平衡法對單位寬度的條塊進(jìn)行受力分析,將波浪產(chǎn)生的最大剪切應(yīng)力作為下滑應(yīng)力作用于條塊底部[1,13],條塊的重力則分解為平行于滑動(dòng)面的下滑力和垂直于滑動(dòng)面的壓力,與傳統(tǒng)的條分法處理方式相同,這種方法只考慮了坡體重力產(chǎn)生的豎向應(yīng)力作用,而忽略了土體水平向應(yīng)力的影響[14],這是目前海底斜坡穩(wěn)定性極限平衡法分析中普遍存在的問題[2,9—12]。

    對于海底緩傾角邊坡,因其坡度較小,重力作用下土體的水平向應(yīng)力場可近似通過靜止側(cè)壓力系數(shù)進(jìn)行計(jì)算[14]。在波浪導(dǎo)致海底緩傾角無限邊坡的穩(wěn)定性分析中,除了需要考慮波浪產(chǎn)生的最大剪切應(yīng)力外,是否需要考慮坡體中的水平向應(yīng)力,其對穩(wěn)定性分析結(jié)果的影響如何,目前的文獻(xiàn)中尚未見分析。本文在波浪作用下海底緩傾角無限坡穩(wěn)定性分析的極限平衡法基礎(chǔ)上,提出基于坡體應(yīng)力狀態(tài)進(jìn)行穩(wěn)定性分析的方法(應(yīng)力狀態(tài)法),并對該方法的適用性進(jìn)行了分析,通過具體算例對極限平衡法和應(yīng)力狀態(tài)法的計(jì)算結(jié)果進(jìn)行了對比,對波致海底緩傾角無限坡的滑動(dòng)失穩(wěn)特征進(jìn)行了探討。

    2 波致海底無限坡穩(wěn)定分析的極限平衡法

    將海床視為均質(zhì)多孔彈性介質(zhì),海床內(nèi)的滲流符合Darcy定律,波浪作用下海床的應(yīng)力變形滿足Biot固結(jié)方程,波浪為線性波條件下,根據(jù)Yamamoto[15]的理論分析,可得無限深海床應(yīng)力場的簡化計(jì)算公式為:式中,σ′z、σ′x、τ′xz分別為海床面以下某深度處土體的波致豎向正應(yīng)力、水平向正應(yīng)力及波致剪切應(yīng)力,p0是海床表面波浪壓力(波底壓力)幅值,k=2π/L為波數(shù),L為波長,γw為水的容重,z為海床面以下的深度坐標(biāo),x為水平向坐標(biāo),ω為波浪角頻率,H為波高,d為水深。從式(2)可知,海床面以下某深度z處的波致最大剪應(yīng)力為τ′xz-m=p0kze-kz。

    在波浪循環(huán)荷載作用下,不考慮孔隙水壓力的累積效應(yīng),在波峰和波谷處分別出現(xiàn)最大豎向正應(yīng)力和最大水平向正應(yīng)力;在波節(jié)點(diǎn)處剪切應(yīng)力最大,而水平正應(yīng)力和豎向正應(yīng)力為零。

    對于海底無限坡的穩(wěn)定性計(jì)算,除了利用有限元進(jìn)行應(yīng)力變形分析外[8],一般采用簡單的極限平衡法[2,9—12],取單位寬度的土體條塊進(jìn)行受力平衡分析。如圖1所示,波浪作用下,取位于波節(jié)點(diǎn)位置的土條進(jìn)行受力分析,土條底部的下滑力由兩部分組成[1,13],一部分為波浪在滑面處產(chǎn)生的最大剪切應(yīng)力形成的下滑力Tw,另一部分是重力在滑動(dòng)方向的下滑力分量Gt,則總的下滑力為:

    土條的抗滑力為R,在總應(yīng)力法分析中,R=Cul,Cu為土體的不排水抗剪強(qiáng)度;在有效應(yīng)力法的分析中,由于重力在垂直滑面方向的壓力分量為Gn=γ′zlcos2β,則抗滑力根據(jù)Mohr-Coulomb抗剪強(qiáng)度公式進(jìn)行計(jì)算,其值為R=(c′+γ′zcos2βtanφ′)l。

    圖1 波致海底無限坡極限平衡法計(jì)算示意圖Fig.1 Sketch of wave-induced stability analysis on submarine infinite slopes by limit equilibrium method

    因此,極限平衡法計(jì)算中,總應(yīng)力分析時(shí)安全系數(shù)的表達(dá)式為:

    有效應(yīng)力分析時(shí)安全系數(shù)的表達(dá)式為:

    式中,γ′為斜坡土體的有效容重,c′和φ′分別為斜坡土體的有效黏聚力和有效內(nèi)摩擦角。

    波浪作用下海底無限坡穩(wěn)定性計(jì)算的極限平衡法中,考慮了坡體相應(yīng)位置的豎向應(yīng)力(重力產(chǎn)生的豎向應(yīng)力場),而沒有考慮與之對應(yīng)的水平向應(yīng)力[14]。波浪作用下海底無限坡穩(wěn)定性計(jì)算中由于是對波節(jié)點(diǎn)位置的土條進(jìn)行受力平衡分析,而波節(jié)點(diǎn)處波浪產(chǎn)生的剪切應(yīng)力最大,水平向正應(yīng)力和豎向正應(yīng)力均為零,因此,波節(jié)點(diǎn)處土體水平向正應(yīng)力主要是坡體重力作用引起的。

    3 基于坡體應(yīng)力狀態(tài)的穩(wěn)定性計(jì)算方法

    3.1 穩(wěn)定性分析思路與計(jì)算公式

    如上所述,極限平衡法中未考慮土體水平向應(yīng)力對穩(wěn)定性分析的影響,對于海底緩傾角斜坡,重力作用下坡體水平向應(yīng)力場可通過土體的靜止側(cè)壓力系數(shù)K0進(jìn)行近似計(jì)算[14]。

    緩傾角斜坡自重作用下某深度處的水平向有效應(yīng)力可近似計(jì)算為:

    式中,K0為土體靜止側(cè)壓力系數(shù),與土體泊松比v有關(guān),根據(jù)彈性理論

    土體靜止側(cè)壓力系數(shù)K0的范圍一般為0.4~1.0[15]。

    在海底無限坡穩(wěn)定性計(jì)算的極限平衡法基礎(chǔ)上,基于坡體的應(yīng)力狀態(tài),即同時(shí)考慮坡體水平向應(yīng)力和豎向應(yīng)力的影響,并結(jié)合波浪導(dǎo)致的最大剪切應(yīng)力,根據(jù)滑面處土體的應(yīng)力平衡條件計(jì)算沿滑面的切向下滑應(yīng)力及垂直滑面的正應(yīng)力,并在此基礎(chǔ)上計(jì)算坡體的安全系數(shù),這種基于坡體應(yīng)力狀態(tài)進(jìn)行穩(wěn)定性分析的方法,此處稱為“應(yīng)力狀態(tài)法”,其具體計(jì)算公式如下。

    設(shè)海底緩傾角無限坡的坡角為β,坡體重力作用下的水平向有效正應(yīng)力通過K0計(jì)算,在波節(jié)點(diǎn)位置,根據(jù)滑面處微元體的應(yīng)力平衡條件,可得平行于滑面的下滑應(yīng)力τ及垂直于滑面的正應(yīng)力σ′分別為:

    故應(yīng)力狀態(tài)法計(jì)算中,總應(yīng)力分析時(shí)安全系數(shù)的表達(dá)式為:

    τ和σ′的計(jì)算分別如式(9)和式(10)中所示。

    對比波浪作用下海底無限坡滑動(dòng)穩(wěn)定性分析的極限平衡法安全系數(shù)計(jì)算公式,可發(fā)現(xiàn)兩種穩(wěn)定性分析方法中,其安全系數(shù)最終計(jì)算表達(dá)式物理意義的實(shí)質(zhì)是相同的,即二者都是滑動(dòng)面上波致最大剪應(yīng)力位置處坡體抗滑應(yīng)力與下滑應(yīng)力的比值,只是坡體應(yīng)力的計(jì)算方法有所不同,極限平衡法是基于不變形體的受力平衡進(jìn)行下滑力和抗滑力的計(jì)算,忽略了坡體水平向外力的作用,由此得到的安全系數(shù)表達(dá)式中只包含了坡體豎向應(yīng)力分量和波致剪應(yīng)力的影響;而應(yīng)力狀態(tài)法則是直接基于外力作用下變形體的應(yīng)力計(jì)算結(jié)果,根據(jù)坡體的應(yīng)力平衡條件進(jìn)行下滑力和抗滑力的計(jì)算,其安全系數(shù)表達(dá)式中不僅包含坡體豎向應(yīng)力分量和波致剪應(yīng)力,還包含水平向應(yīng)力分量的影響。

    綜上,無論是極限平衡法還是應(yīng)力狀態(tài)法,二者在分析波致海底緩傾角無限邊坡滑動(dòng)穩(wěn)定性時(shí),其安全系數(shù)計(jì)算式均體現(xiàn)為應(yīng)力表達(dá)的抗滑力與下滑力的比值,兩種方法安全系數(shù)所表達(dá)的物理意義是相同的;且應(yīng)力狀態(tài)法的計(jì)算是基于變形體的應(yīng)力分析,更符合坡體的實(shí)際物理狀態(tài),因此,將該方法用于波致海底緩傾角無限坡的滑動(dòng)穩(wěn)定性計(jì)算是可行的。

    3.2 應(yīng)力狀態(tài)法的適用范圍

    由應(yīng)力狀態(tài)法計(jì)算公式可知,坡體水平向正應(yīng)力是通過土體的靜止側(cè)壓力系數(shù)K0計(jì)算得到的,當(dāng)坡角為零時(shí),根據(jù)彈性理論該計(jì)算是正確的,但隨著邊坡坡度的增加(坡角的增大),通過K0計(jì)算得到的坡體水平向應(yīng)力場與其實(shí)際的應(yīng)力狀態(tài)相差會(huì)逐漸增大,故應(yīng)力狀態(tài)法的可靠度在較大程度上依賴于對坡體水平向應(yīng)力場計(jì)算的準(zhǔn)確度,亦即其適用范圍與坡角密切相關(guān)。

    為了分析應(yīng)力狀態(tài)法適用的坡角范圍,將坡體視為彈性介質(zhì),利用Sigma/W有限元軟件[16]計(jì)算重力作用下無限坡的應(yīng)力場(可視為實(shí)際應(yīng)力場),并將其與利用K0計(jì)算得到的坡體水平向應(yīng)力場進(jìn)行比較,分析計(jì)算誤差。

    如圖2所示為有限元計(jì)算的幾何模型,邊坡水平向長400 m,坡腳處土體厚度20 m,坡頂處的厚度依

    有效應(yīng)力分析時(shí)安全系數(shù)的表達(dá)式為:據(jù)坡度的不同而變化。計(jì)算分析了坡角β為3°、4°和6°時(shí)的邊坡應(yīng)力場,與傾角為0°時(shí)的邊坡應(yīng)力場相比,坡體的豎向應(yīng)力場變化很小,水平向應(yīng)力場依據(jù)坡角的不同,在不同的深度范圍內(nèi)有一定變化。總體而言,與K0計(jì)算的坡體水平應(yīng)力相比,靠近坡腳附近的水平應(yīng)力場偏大,而靠近坡頂附近的水平應(yīng)力場偏小,考慮到此處分析的無限坡的性質(zhì),取水平方向坡體中間位置處的截面(如圖2中的A-A截面)進(jìn)行水平應(yīng)力的計(jì)算誤差分析,結(jié)果如表1所示。

    圖2 有限元計(jì)算的幾何模型Fig.2 Geometric model for finite element analysis

    表1 水平應(yīng)力計(jì)算誤差(單位:%)Tab.1 Error analysis of horizontal stress(unit:%)

    因坡體中部A-A截面利用K0計(jì)算出的水平應(yīng)力數(shù)值小于實(shí)際的水平應(yīng)力,因此表1中的誤差百分?jǐn)?shù)為負(fù)值。從表1中數(shù)值可知,應(yīng)力狀態(tài)法的適用坡角與潛在滑動(dòng)面深度有關(guān),如果以水平應(yīng)力計(jì)算誤差10%為判斷標(biāo)準(zhǔn)(由此引起的安全系數(shù)的計(jì)算誤差分析見后述4.3節(jié)內(nèi)容),則當(dāng)滑動(dòng)面深度大于5 m時(shí),其適用的坡角范圍為0°~6°;滑動(dòng)面深度3~5 m時(shí),其適用坡角范圍為0°~4°;滑動(dòng)面深度2~3 m時(shí),其適用坡角為0°~3°。

    4 算例分析與討論

    分別利用極限平衡法與本文提出的應(yīng)力狀態(tài)法,通過某具體算例,從海底緩傾角無限坡的穩(wěn)定性隨滑動(dòng)面深度的變化特征、土體泊松比和斜坡坡度對穩(wěn)定性計(jì)算結(jié)果的影響等方面,對波致海底緩傾角無限坡滑動(dòng)穩(wěn)定性的計(jì)算方法進(jìn)行對比分析,同時(shí)對波浪作用下海底緩傾角無限坡的失穩(wěn)滑動(dòng)特征進(jìn)行探討。

    某海底緩傾角無限邊坡,坡角β=1.8°,土體浮容重γ′=9.4 kN/m3,剪切模量G=3.5×106Pa,泊松比ν=0.35,有效黏聚力c′=1.8 kPa,有效內(nèi)摩擦角φ′=6°;波浪參數(shù)為波高5.2 m、周期8.6 s、水深8 m,對應(yīng)的波長70.6 m。

    4.1 波致海底緩傾角無限坡穩(wěn)定性隨深度的變化特征

    基于前述波致剪應(yīng)力的表達(dá)式,在0.159L深度處剪應(yīng)力達(dá)最大值,其值為0.368P0。此處算例中波長L為70.6 m,最大波致剪應(yīng)力發(fā)生在11 m深度處,為此,對距坡面0.5~13 m深度范圍內(nèi)的坡體穩(wěn)定性及其隨深度的變化特征進(jìn)行分析。

    如圖3所示為坡體安全系數(shù)隨潛在滑動(dòng)面深度的變化曲線,可以看出,應(yīng)力狀態(tài)法計(jì)算的安全系數(shù)大于極限平衡法的計(jì)算結(jié)果,且隨著潛在滑動(dòng)面深度的增加,兩種計(jì)算方法所得安全系數(shù)的差值也在逐漸增加,在深度4 m處二者相差11%,在13 m位置處達(dá)22%。

    圖3 安全系數(shù)隨滑動(dòng)面深度的變化Fig.3 Variations of safety factor with different depths of slip surface

    根據(jù)兩種安全系數(shù)計(jì)算方法的表達(dá)式可知,應(yīng)力狀態(tài)法考慮了坡體水平向應(yīng)力場的影響,由于該應(yīng)力場是重力引起的水平向壓應(yīng)力場,而水平壓力場對坡體的橫向滑動(dòng)變形具有一定的約束作用,因此在應(yīng)力狀態(tài)法安全系數(shù)計(jì)算表達(dá)式中,水平向應(yīng)力的出現(xiàn)明顯減小了沿滑面的下滑力,其計(jì)算所得的安全系數(shù)較極限平衡法偏大,且隨著滑面計(jì)算深度的增加,水平向應(yīng)力場逐漸增大,其對坡體滑動(dòng)變形的約束效應(yīng)也逐漸增加,導(dǎo)致隨著深度增加兩種方法計(jì)算結(jié)果的差異逐漸增大。

    兩種計(jì)算方法所得安全系數(shù)均隨滑面深度的增加出現(xiàn)先減小后增大的變化趨勢,在距坡面4 m深度處安全系數(shù)達(dá)最小值,此時(shí)極限平衡法計(jì)算的安全系數(shù)為0.92,應(yīng)力狀態(tài)法的安全系數(shù)為1.02,表明斜坡最易在該深度處發(fā)生失穩(wěn)滑動(dòng)。對安全系數(shù)隨深度的變化數(shù)值進(jìn)行分析可知,在距坡面3~6 m的深度范圍內(nèi),該斜坡滑動(dòng)穩(wěn)定系數(shù)的數(shù)值非常接近且都相對較小,基本都處于最危險(xiǎn)滑動(dòng)區(qū)域,因此斜坡體極有可能沿具有一定厚度的滑動(dòng)帶失穩(wěn)滑動(dòng),而不只是沿單一的滑動(dòng)面發(fā)生滑動(dòng)。

    本例中最大波致剪切應(yīng)力發(fā)生在距坡面11 m的深度處,而穩(wěn)定性計(jì)算結(jié)果表明,最易滑移的位置發(fā)生在距坡面3~6 m的深度范圍,遠(yuǎn)小于最大波致剪應(yīng)力的深度位置,這是因?yàn)殡S著深度的增加,斜坡土體的抗剪能力也在逐漸增大,在有效應(yīng)力分析中體現(xiàn)在斜坡滑動(dòng)面上的有效正應(yīng)力隨著深度增加而增大,進(jìn)而增大了抗剪強(qiáng)度;而在總應(yīng)力分析中,斜坡土體的不排水抗剪強(qiáng)度一般會(huì)隨著深度的增加而增大[6,11—12],故波致剪應(yīng)力最大的深度位置往往不是斜坡體最易失穩(wěn)滑移的位置。

    4.2 泊松比對波致海底緩傾角無限坡穩(wěn)定性計(jì)算的影響

    以滑動(dòng)面的深度4 m為例(該位置處的坡體安全系數(shù)最?。?,繪制海底緩傾角無限坡安全系數(shù)隨土體泊松比的變化曲線,如圖4所示。

    圖4 安全系數(shù)隨土體泊松比的變化Fig.4 Variations of safety factor with different Poisson ratios of soil

    從圖4中可以看出,極限平衡法中因沒有考慮重力引起的水平向應(yīng)力作用,其計(jì)算出的安全系數(shù)是一個(gè)常值0.92,不隨泊松比變化。應(yīng)力狀態(tài)法中考慮了水平向應(yīng)力作用,且該水平應(yīng)力的數(shù)值直接受泊松比的影響,因此,應(yīng)力狀態(tài)法計(jì)算的坡體安全系數(shù)隨著泊松比的變化而有較明顯的改變,當(dāng)泊松比為0.3時(shí),其安全系數(shù)為1.00,與極限平衡法所得安全系數(shù)相差9%,當(dāng)泊松比為0.5時(shí),其安全系數(shù)為1.13,與極限平衡法的安全系數(shù)相差達(dá)23%。

    從應(yīng)力狀態(tài)法的安全系數(shù)計(jì)算公式可以看出,當(dāng)泊松比較小時(shí),坡體水平向壓應(yīng)力的影響相對減弱,由此導(dǎo)致水平向壓應(yīng)力場對坡體滑動(dòng)變形的約束作用隨之降低,其安全系數(shù)計(jì)算結(jié)果與極限平衡法的計(jì)算結(jié)果就越相近,故泊松比越小,兩種方法的計(jì)算結(jié)果相差越??;反之則相差越大。

    4.3 坡度對波致海底緩傾角無限坡穩(wěn)定性計(jì)算的影響

    改變海底斜坡的坡度,分析坡度的變化對兩種穩(wěn)定性分析方法計(jì)算結(jié)果的影響。以滑動(dòng)面深度4 m為例進(jìn)行分析,則根據(jù)前述應(yīng)力狀態(tài)法適用范圍的分析結(jié)果,選取坡角變化范圍0.3°~4.0°,計(jì)算結(jié)果如圖5所示。

    圖5 安全系數(shù)隨坡角的變化Fig.5 Variations of safety factor with different slope angles

    根據(jù)圖5可知,隨著坡角的增大,兩種計(jì)算方法所得安全系數(shù)的差異逐漸增大,當(dāng)坡角為1.8°時(shí),相差11%,而當(dāng)坡角為4.0°時(shí)相差22%。隨著坡角的增大,重力對坡體穩(wěn)定性計(jì)算結(jié)果的影響逐漸增加,進(jìn)一步由重力作用導(dǎo)致的水平向壓應(yīng)力場對坡體滑動(dòng)變形的約束效應(yīng)相對隨之增強(qiáng),應(yīng)力狀態(tài)法計(jì)算中考慮了水平向應(yīng)力場的這種作用,而極限平衡法沒有考慮,故導(dǎo)致隨著坡角的增加,兩種穩(wěn)定性分析方法計(jì)算結(jié)果之間的差異隨之增大。

    如前所述,應(yīng)力狀態(tài)法計(jì)算穩(wěn)定性時(shí)的坡角適用范圍,是以K0計(jì)算的坡體水平應(yīng)力與實(shí)際水平應(yīng)力相差不大于10%為判斷標(biāo)準(zhǔn)而得出的,則當(dāng)水平應(yīng)力計(jì)算誤差為10%時(shí),其對應(yīng)的安全系數(shù)的誤差是否可以接受,此處將通過算例進(jìn)行分析。

    以海底無限坡坡角4°為例,當(dāng)潛在滑動(dòng)面深度為4 m時(shí),利用K0計(jì)算坡體水平應(yīng)力,此時(shí)應(yīng)力狀態(tài)法計(jì)算得到的安全系數(shù)為0.91;當(dāng)坡體的實(shí)際水平應(yīng)力與利用K0計(jì)算得到的水平應(yīng)力相差10%時(shí),利用坡體實(shí)際水平應(yīng)力計(jì)算得到的安全系數(shù)為0.93,與直接利用K0計(jì)算得到的安全系數(shù)相差2%,誤差較小。上述分析表明,當(dāng)K0計(jì)算的坡體水平應(yīng)力與實(shí)際水平應(yīng)力相差不超過10%時(shí),對坡體安全系數(shù)計(jì)算的影響較小,其誤差是可以接受的。

    5 結(jié)論

    對目前波致海底無限坡滑動(dòng)穩(wěn)定性計(jì)算的極限平衡法中存在的問題進(jìn)行了分析,并在此基礎(chǔ)上,提出了基于坡體應(yīng)力狀態(tài)的穩(wěn)定性計(jì)算方法(應(yīng)力狀態(tài)法)及其適用范圍,通過具體算例,對極限平衡法和應(yīng)力狀態(tài)法的計(jì)算結(jié)果進(jìn)行了分析對比,對波致海底緩傾角無限坡的滑動(dòng)失穩(wěn)特征進(jìn)行了探討,主要得到了以下結(jié)論。

    (1)提出了波浪作用下海底緩傾角無限坡滑動(dòng)穩(wěn)定性分析的應(yīng)力狀態(tài)法,并分析了其適用范圍;

    (2)應(yīng)力狀態(tài)法計(jì)算的安全系數(shù)大于極限平衡法的安全系數(shù),且隨著滑動(dòng)面深度的增加、土體泊松比以及邊坡坡角的增大,兩種計(jì)算方法得出的安全系數(shù)的差異逐漸增大;

    (3)對于海底緩傾角無限邊坡,在波浪作用下失穩(wěn)時(shí),極有可能沿具有一定厚度的滑動(dòng)帶失穩(wěn)滑動(dòng),而不只是沿著單一的滑動(dòng)面滑動(dòng);

    (4)對于波浪作用下的海底緩傾角無限邊坡,由于坡體抗剪強(qiáng)度的影響,波致最大剪應(yīng)力所在的深度位置,常常不是斜坡體最易失穩(wěn)滑移的位置。

    [1] 葉銀燦.中國海洋災(zāi)害地質(zhì)學(xué)[M].北京:海洋出版社,2012.

    Ye Yingcan.Marine hazard geology[M].Beijing:China Ocean Press,2012.

    [2] Leynaud D,Mienert J,Nadim F.Slope stability assessment of the Helland Hansen area offshore the mid-Norwegian margin[J].Marine Geology,2004,213(1/4):457-480.

    [3] 張亮,欒錫武.南海北部陸坡穩(wěn)定性定量分析[J].地球物理學(xué)進(jìn)展,2012,27(4):1443-1453.

    Zhang Liang,Luan Xiwu.Quantitative analysis of submarine slope stability on the northern slope of the South China Sea[J].Progress in Geophsics,2012,27(4):1443-1453.

    [4] 楊作升,陳衛(wèi)民,陳彰榕,等.黃河口水下滑坡體系[J].海洋與湖沼,1994,25(6):573-581.

    Yang Zuosheng,Chen Weimin,Chen Zhangrong,et al.Subaqueous landslide system in the Huanghe River(Yellow River)delta[J].Oceanologia et Limnologia Sinica,1994,25(6):573-581.

    [5] 甘華陽,王家生,胡高韋.海洋沉積物中的天然氣水合物與海底滑坡[J].防災(zāi)減災(zāi)工程學(xué)報(bào),2004,24(2):177-181.

    Gan Huayang,Wang Jiasheng,Hu Gaowei.Submarine landslide related to natural gas hydrate within benthal deposit[J].Journal of Disaster Prevention and Mitigation Engineering,2004,24(2):177-181.

    [6] Henkel D J.The role of waves in causing submarine landslides[J].Geotechnique,1970,20(1):75-80.

    [7] 王樹青,梁丙臣.海洋工程波浪力學(xué)[M].青島:中國海洋大學(xué)出版社,2013.

    Wang Shuqing,Liang Bingchen.Wave mechanics for ocean engineering[M].Qingdao:China Ocean University Press,2013.

    [8] Wright S G,Dunham R S.Bottom stability under wave induced loading[C]//Proceeding 4th Annual Offshore Technology Conference,Houston,1972:853-862.

    [9] L'Heureux JS,Vanneste M,Rise L,et al.Stability,mobility and failure mechanism for landslides at the upper continental slopeoff Vesteralen,Norway[J].Marine Geology,2013,346:192-207.

    [10] Ikari M J,Strasser M,Saffer D M,et al.Submarine landslide potential near the megasplay fault at the Nankai subduction zone[J].Earth and Planetary Science Letters,2011,312:453-462.

    [11] 顧小蕓.粘質(zhì)海底穩(wěn)定性實(shí)例分析[J].工程地質(zhì)學(xué)報(bào),1996,4(1):32-38.

    Gu Xiaoyun.Case study of clayey seabottom stability[J].Journal of Engineering Geology,1996,4(1):32-38.

    [12] 葉銀燦,陳錫土,宋連清,等.浙江北部島嶼海域土體穩(wěn)定性研究[J].東海海洋,1996,14(1):1-18.

    Ye Yincan,Chen Xitu,Song Lianqing,et al.A study for the seafloor sediment stability around the island areas of northern Zhejiang Province[J].Donghai Marine Science,1996,14(1):1-18.

    [13] Hampton M A,Lee H J,Locat J.Submarine landslides[J].Reviews of Geophysics,1996,34(1):33-59.

    [14] Lade P V.Static instability and liquefaction of loose fine sandy slopes[J].Journal of Geotechnical Engineering,1992,118(1):51-71.

    [15] Yamamoto T.Sea bed instability from waves[C]//Proceeding 10th Annual Offshore Technology Conference,Houston,1978:1819-1828.

    [16] GEO-SLOPE International Ltd.Stress-deformation modeling with Sigma/W 2007[M].Canada:GEO-SLOPE International Ltd,2008.

    Wave-induced stability analysis on submarine infinite slope with gentle dip angle

    Liu Xiaoli1,2,Dou Jinzhong1,Ying Zi1,Huo Yandong1,Zhang Xu1
    (1.College of Environmental Science and Engineering,Ocean University of China,Qingdao 266100,China;2.Key Laboratory of Marine Environment and Ecology,Ministry of Education,Qingdao 266100,China)

    For wave-induced stability analysis on submarine infinite slopes,the limit equilibrium method generally ignores horizontal stresses acting in the soil deposit.On basis of describing the characteristics of the submarine infinite slope under wave action and considering the soil stress state at the slip surface,this study proposed a stress analysis method to analyze sliding stability of the submarine infinite slope,and discussed the applicable prospect of the method.A case study of a submarine infinite slope has been conducted by both of the limit equilibrium method and the proposed stress analysis method.The results have shown that the factor of safety by the stress analysis method is greater than that of the limit equilibrium method,and the differences of factors obtained by the two methods become larger with increasing depth of the slip surface,Poisson ratio of the soil deposit and dip angle of the infinite slope.For a submarine infinite slope with gentle dip angle under wave action,it may slide along a slip zone rather than a single slip surface under unstable state,and the depth for the maximum shear stress does not generally indicate the most instability location of the infinite slope.

    submarine infinite slope;wave action;stability analysis;gentle dip angle

    P642.22

    A

    0253-4193(2015)03-0099-07

    劉小麗,竇錦鐘,英姿,等.波致海底緩傾角無限坡滑動(dòng)穩(wěn)定性計(jì)算分析探討[J].海洋學(xué)報(bào),2015,37(3):99—105,

    10.3969/j.issn.0253-4193.2015.03.010

    Liu Xiaoli,Dou Jinzhong,Ying Zi,et al.Wave-induced stability analysis on submarine infinite slopes with gentle dip angle[J].Haiyang Xuebao,2015,37(3):99—105,doi:10.3969/j.issn.0253-4193.2015.03.000

    2014-05-05;

    2014-08-30。

    國家自然科學(xué)基金項(xiàng)目(41272316)。

    劉小麗(1974—),女,山東省青島市人,副教授,主要從事海洋地質(zhì)災(zāi)害相關(guān)方面的研究。E-mail:LXL4791@163.com

    猜你喜歡
    坡角坡體波浪
    坡角對雙坡屋蓋風(fēng)荷載特性影響分析
    降雨對庫區(qū)邊坡入滲規(guī)律的影響研究
    波浪谷和波浪巖
    采動(dòng)-裂隙水耦合下含深大裂隙巖溶山體失穩(wěn)破壞機(jī)理
    烏弄龍水電站庫區(qū)拉金神谷坡體變形成因機(jī)制分析
    波浪谷隨想
    不同開采位置對邊坡穩(wěn)定性影響的數(shù)值模擬分析
    山西煤炭(2019年2期)2019-08-29 05:35:40
    去看神奇波浪谷
    抗滑樁與邊坡夾角對抗滑樁抗滑效果研究
    四川水泥(2017年8期)2017-08-30 17:04:43
    山地光伏電站間距計(jì)算與合理選址
    太陽能(2015年9期)2015-01-01 03:00:20
    757午夜福利合集在线观看| 丰满饥渴人妻一区二区三| 免费人妻精品一区二区三区视频| 午夜福利一区二区在线看| 九色亚洲精品在线播放| 亚洲av成人一区二区三| 夫妻午夜视频| bbb黄色大片| 国产福利在线免费观看视频| 精品人妻在线不人妻| 国产亚洲一区二区精品| 成人精品一区二区免费| 国产亚洲欧美精品永久| 久久精品91无色码中文字幕| 成年人免费黄色播放视频| 搡老岳熟女国产| 亚洲精品在线美女| 男男h啪啪无遮挡| 多毛熟女@视频| 操出白浆在线播放| 久久国产精品大桥未久av| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区三区在线| 咕卡用的链子| 国产日韩一区二区三区精品不卡| 欧美人与性动交α欧美精品济南到| 极品教师在线免费播放| 欧美黑人精品巨大| 国产无遮挡羞羞视频在线观看| 一本综合久久免费| 两个人看的免费小视频| 国产区一区二久久| 亚洲av国产av综合av卡| 亚洲av电影在线进入| 久久毛片免费看一区二区三区| 伊人久久大香线蕉亚洲五| 亚洲专区国产一区二区| 成人精品一区二区免费| 正在播放国产对白刺激| 国产成人免费观看mmmm| 18在线观看网站| 老司机影院毛片| 午夜福利一区二区在线看| 嫩草影视91久久| 国产亚洲av高清不卡| 国产免费av片在线观看野外av| 久久青草综合色| 日韩中文字幕视频在线看片| 麻豆国产av国片精品| 国产免费福利视频在线观看| 亚洲成av片中文字幕在线观看| 免费黄频网站在线观看国产| 成人av一区二区三区在线看| 久热爱精品视频在线9| 女同久久另类99精品国产91| 久热爱精品视频在线9| 精品福利观看| 五月天丁香电影| 丰满人妻熟妇乱又伦精品不卡| 正在播放国产对白刺激| 老司机午夜福利在线观看视频 | 国产精品av久久久久免费| 99国产综合亚洲精品| 精品一品国产午夜福利视频| 999久久久国产精品视频| 日韩免费av在线播放| 老熟妇仑乱视频hdxx| 欧美日韩黄片免| www日本在线高清视频| 久久国产精品人妻蜜桃| 国产成人欧美在线观看 | 色精品久久人妻99蜜桃| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美一区二区综合| 国产成人精品久久二区二区91| 男女免费视频国产| 久久精品人人爽人人爽视色| 99re6热这里在线精品视频| 国产又色又爽无遮挡免费看| 涩涩av久久男人的天堂| 亚洲专区中文字幕在线| 久久久久精品人妻al黑| 黄片小视频在线播放| 日本五十路高清| 久久久久久人人人人人| 日韩制服丝袜自拍偷拍| 精品亚洲乱码少妇综合久久| 韩国精品一区二区三区| av线在线观看网站| videos熟女内射| 美女午夜性视频免费| 12—13女人毛片做爰片一| 久久国产精品影院| 久久久久久免费高清国产稀缺| 国产老妇伦熟女老妇高清| 高清欧美精品videossex| av电影中文网址| 我要看黄色一级片免费的| 丝袜美足系列| 久久影院123| 2018国产大陆天天弄谢| 国产麻豆69| 国产精品成人在线| 大香蕉久久成人网| 欧美午夜高清在线| 国产欧美日韩一区二区三区在线| 精品人妻1区二区| 97在线人人人人妻| 99久久人妻综合| 亚洲国产看品久久| 在线永久观看黄色视频| 久久精品国产亚洲av香蕉五月 | 黄网站色视频无遮挡免费观看| 国产无遮挡羞羞视频在线观看| 婷婷成人精品国产| 纯流量卡能插随身wifi吗| 欧美日韩亚洲高清精品| 天天躁夜夜躁狠狠躁躁| 亚洲 国产 在线| 精品亚洲乱码少妇综合久久| 在线观看免费视频日本深夜| 丰满少妇做爰视频| 窝窝影院91人妻| 我的亚洲天堂| 久久午夜综合久久蜜桃| 777米奇影视久久| 精品一区二区三区四区五区乱码| 亚洲人成伊人成综合网2020| 亚洲免费av在线视频| 亚洲综合色网址| 精品国产一区二区久久| 亚洲五月婷婷丁香| 亚洲av欧美aⅴ国产| 丝袜美腿诱惑在线| 国产精品九九99| 亚洲成人国产一区在线观看| 少妇被粗大的猛进出69影院| 国产精品影院久久| 久久久久久久大尺度免费视频| 国产高清视频在线播放一区| 美女视频免费永久观看网站| 亚洲国产成人一精品久久久| 女同久久另类99精品国产91| 少妇 在线观看| 在线观看免费视频日本深夜| 亚洲专区国产一区二区| 亚洲成人免费av在线播放| 天堂俺去俺来也www色官网| 国产精品98久久久久久宅男小说| 不卡av一区二区三区| 欧美日韩国产mv在线观看视频| 一级毛片女人18水好多| 91大片在线观看| 高清av免费在线| 国产97色在线日韩免费| 18禁黄网站禁片午夜丰满| 欧美一级毛片孕妇| 看免费av毛片| 9热在线视频观看99| 丝袜美足系列| 9191精品国产免费久久| netflix在线观看网站| 我要看黄色一级片免费的| 国产精品久久久久成人av| 一本综合久久免费| 国产亚洲精品一区二区www | 久久久久视频综合| 高清在线国产一区| 国产精品一区二区在线不卡| 我的亚洲天堂| 丰满少妇做爰视频| 91九色精品人成在线观看| 怎么达到女性高潮| 中文字幕精品免费在线观看视频| 国产精品久久久久成人av| 高清视频免费观看一区二区| av有码第一页| 男女之事视频高清在线观看| 久久中文字幕一级| 国产福利在线免费观看视频| 久久久久久久精品吃奶| av网站免费在线观看视频| av超薄肉色丝袜交足视频| 中国美女看黄片| 日本av免费视频播放| 久久久精品国产亚洲av高清涩受| 久久ye,这里只有精品| 国精品久久久久久国模美| 国产1区2区3区精品| 午夜福利影视在线免费观看| 又紧又爽又黄一区二区| 飞空精品影院首页| 婷婷成人精品国产| 精品人妻1区二区| 午夜福利一区二区在线看| 91av网站免费观看| 成年人黄色毛片网站| 日韩大片免费观看网站| 精品一区二区三区av网在线观看 | 中文字幕av电影在线播放| 啦啦啦视频在线资源免费观看| 亚洲 欧美一区二区三区| 99国产精品99久久久久| av电影中文网址| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品美女久久久久99蜜臀| 亚洲久久久国产精品| 中文字幕人妻丝袜一区二区| 美女高潮到喷水免费观看| 国产成人精品久久二区二区91| 亚洲国产欧美一区二区综合| 一边摸一边抽搐一进一小说 | 丝袜美腿诱惑在线| 国产福利在线免费观看视频| 精品福利永久在线观看| 又大又爽又粗| 国产精品美女特级片免费视频播放器 | 久久中文看片网| 狠狠婷婷综合久久久久久88av| 两性夫妻黄色片| 精品久久蜜臀av无| 日本vs欧美在线观看视频| 国产精品亚洲av一区麻豆| 老司机福利观看| 国产一区二区 视频在线| 国产精品电影一区二区三区 | 精品卡一卡二卡四卡免费| 色综合婷婷激情| 国产成+人综合+亚洲专区| 黑人巨大精品欧美一区二区mp4| 一边摸一边做爽爽视频免费| 天天躁夜夜躁狠狠躁躁| 悠悠久久av| 国产xxxxx性猛交| 亚洲国产毛片av蜜桃av| 精品国产一区二区三区四区第35| 在线观看免费视频网站a站| 99re6热这里在线精品视频| 天天添夜夜摸| a级毛片在线看网站| 欧美日韩福利视频一区二区| 夫妻午夜视频| 99国产精品99久久久久| 9热在线视频观看99| 伊人久久大香线蕉亚洲五| 少妇被粗大的猛进出69影院| 欧美日韩中文字幕国产精品一区二区三区 | 日日夜夜操网爽| 成人永久免费在线观看视频 | 人成视频在线观看免费观看| 久久香蕉激情| 成年女人毛片免费观看观看9 | 麻豆av在线久日| 国产成人免费观看mmmm| 中文字幕最新亚洲高清| 黑人操中国人逼视频| 在线观看免费视频日本深夜| 精品亚洲成国产av| 久久av网站| 精品免费久久久久久久清纯 | 高清毛片免费观看视频网站 | av不卡在线播放| 亚洲三区欧美一区| 老熟女久久久| 国产精品香港三级国产av潘金莲| 免费在线观看完整版高清| 老鸭窝网址在线观看| 三级毛片av免费| 国产高清国产精品国产三级| 欧美人与性动交α欧美软件| 国产一区二区三区视频了| 国产一区有黄有色的免费视频| 亚洲第一青青草原| 久久久久久久精品吃奶| 成人黄色视频免费在线看| 精品国产一区二区久久| 成人影院久久| 亚洲自偷自拍图片 自拍| 老汉色∧v一级毛片| 老熟女久久久| 免费在线观看完整版高清| 精品福利永久在线观看| 亚洲免费av在线视频| 啦啦啦视频在线资源免费观看| 手机成人av网站| 久久久国产欧美日韩av| 午夜福利免费观看在线| 最黄视频免费看| 亚洲欧洲日产国产| 免费在线观看视频国产中文字幕亚洲| 三上悠亚av全集在线观看| 久久天堂一区二区三区四区| 国产av精品麻豆| 视频区图区小说| 亚洲欧美日韩高清在线视频 | 国产伦理片在线播放av一区| 成人亚洲精品一区在线观看| 亚洲天堂av无毛| 青青草视频在线视频观看| 国产精品98久久久久久宅男小说| 欧美精品高潮呻吟av久久| 亚洲精品中文字幕一二三四区 | 国产精品亚洲一级av第二区| 午夜免费成人在线视频| 19禁男女啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 中文字幕最新亚洲高清| 一边摸一边抽搐一进一出视频| 无遮挡黄片免费观看| 99在线人妻在线中文字幕 | 国产福利在线免费观看视频| 最黄视频免费看| 国产日韩一区二区三区精品不卡| av电影中文网址| 女人高潮潮喷娇喘18禁视频| 国产成人精品在线电影| 91国产中文字幕| 多毛熟女@视频| 午夜福利影视在线免费观看| 韩国精品一区二区三区| 国产精品免费一区二区三区在线 | 乱人伦中国视频| 女人久久www免费人成看片| 大香蕉久久成人网| 国产99久久九九免费精品| 少妇粗大呻吟视频| 欧美在线黄色| 肉色欧美久久久久久久蜜桃| 丝袜人妻中文字幕| 久久久久视频综合| 精品高清国产在线一区| 亚洲国产成人一精品久久久| 80岁老熟妇乱子伦牲交| 少妇裸体淫交视频免费看高清 | 啦啦啦免费观看视频1| 黄频高清免费视频| 亚洲欧洲日产国产| 黄色视频在线播放观看不卡| 国产成人欧美| 日韩欧美一区二区三区在线观看 | 国产日韩欧美在线精品| 视频区欧美日本亚洲| 国产麻豆69| 亚洲中文日韩欧美视频| 日本五十路高清| 久久久水蜜桃国产精品网| bbb黄色大片| 美国免费a级毛片| 国产精品影院久久| 精品久久久久久久毛片微露脸| 成人手机av| 久久热在线av| 亚洲欧美日韩另类电影网站| 久久毛片免费看一区二区三区| 狠狠婷婷综合久久久久久88av| 十分钟在线观看高清视频www| 高清欧美精品videossex| 一个人免费看片子| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性xxxx| 色婷婷久久久亚洲欧美| 日韩 欧美 亚洲 中文字幕| 天天操日日干夜夜撸| 国产片内射在线| 国产日韩欧美视频二区| 80岁老熟妇乱子伦牲交| a在线观看视频网站| 91精品国产国语对白视频| tube8黄色片| 黄色视频,在线免费观看| 国产麻豆69| 女人爽到高潮嗷嗷叫在线视频| 久久久国产欧美日韩av| 69av精品久久久久久 | 蜜桃国产av成人99| 国产亚洲欧美在线一区二区| 久久久精品国产亚洲av高清涩受| 天天躁日日躁夜夜躁夜夜| videos熟女内射| 日韩欧美一区视频在线观看| 国产成人精品无人区| 国产亚洲欧美精品永久| 国产精品1区2区在线观看. | 国产精品免费视频内射| 80岁老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| 国产一区二区 视频在线| 日韩免费av在线播放| 视频区欧美日本亚洲| 日韩中文字幕视频在线看片| 国产精品98久久久久久宅男小说| 18禁国产床啪视频网站| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| av欧美777| 两性夫妻黄色片| 一级毛片女人18水好多| 每晚都被弄得嗷嗷叫到高潮| 岛国毛片在线播放| 日本一区二区免费在线视频| 精品国产乱子伦一区二区三区| 伊人久久大香线蕉亚洲五| 精品免费久久久久久久清纯 | av福利片在线| 中文字幕高清在线视频| 午夜精品国产一区二区电影| 亚洲欧美一区二区三区黑人| 成人影院久久| 亚洲精品在线美女| 老司机影院毛片| 夫妻午夜视频| 老司机靠b影院| 亚洲av日韩精品久久久久久密| 免费人妻精品一区二区三区视频| 国产成人av激情在线播放| 日韩免费av在线播放| 99精国产麻豆久久婷婷| 成人国产av品久久久| 黄色丝袜av网址大全| 国产亚洲一区二区精品| 色94色欧美一区二区| 午夜老司机福利片| 在线观看舔阴道视频| 久久毛片免费看一区二区三区| 青草久久国产| 少妇猛男粗大的猛烈进出视频| 亚洲少妇的诱惑av| 亚洲欧美一区二区三区黑人| 国产男女内射视频| 麻豆成人av在线观看| 久久毛片免费看一区二区三区| 免费人妻精品一区二区三区视频| 日日爽夜夜爽网站| 999久久久国产精品视频| 俄罗斯特黄特色一大片| 欧美成人免费av一区二区三区 | 99热国产这里只有精品6| 亚洲成人手机| 丁香欧美五月| 99国产精品一区二区蜜桃av | 性少妇av在线| 757午夜福利合集在线观看| 久久中文看片网| 热99re8久久精品国产| 欧美日韩国产mv在线观看视频| 国产淫语在线视频| 亚洲九九香蕉| 99香蕉大伊视频| 久久精品亚洲精品国产色婷小说| 国产精品欧美亚洲77777| 亚洲自偷自拍图片 自拍| 久久天堂一区二区三区四区| 国产成人一区二区三区免费视频网站| 老司机在亚洲福利影院| 19禁男女啪啪无遮挡网站| 国产亚洲欧美精品永久| 伊人久久大香线蕉亚洲五| 免费看十八禁软件| 精品国产国语对白av| 黄片小视频在线播放| 18禁裸乳无遮挡动漫免费视频| 妹子高潮喷水视频| 男女下面插进去视频免费观看| 亚洲精品久久午夜乱码| 国产又爽黄色视频| 一区二区av电影网| 国产精品1区2区在线观看. | 天天躁狠狠躁夜夜躁狠狠躁| 成年人黄色毛片网站| 日韩免费高清中文字幕av| 国产精品免费一区二区三区在线 | 久久国产精品大桥未久av| 国产成人av激情在线播放| 99re在线观看精品视频| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 一个人免费在线观看的高清视频| 欧美国产精品一级二级三级| 91字幕亚洲| 中文欧美无线码| 不卡av一区二区三区| 成人亚洲精品一区在线观看| 两个人看的免费小视频| 国产不卡一卡二| 极品人妻少妇av视频| 999久久久精品免费观看国产| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| 黑丝袜美女国产一区| 国产aⅴ精品一区二区三区波| 午夜福利视频在线观看免费| 桃红色精品国产亚洲av| 国产精品一区二区精品视频观看| 免费高清在线观看日韩| 亚洲色图综合在线观看| 久久天堂一区二区三区四区| 精品久久蜜臀av无| 成人国产一区最新在线观看| 电影成人av| 一区二区三区乱码不卡18| 亚洲精品一卡2卡三卡4卡5卡| kizo精华| 老司机午夜十八禁免费视频| 精品国产亚洲在线| 久久久久久久久久久久大奶| 亚洲av成人一区二区三| 18禁国产床啪视频网站| 精品熟女少妇八av免费久了| 美国免费a级毛片| 日韩精品免费视频一区二区三区| 又大又爽又粗| 久久久久精品国产欧美久久久| 欧美国产精品va在线观看不卡| 99国产精品免费福利视频| av网站在线播放免费| 色老头精品视频在线观看| 少妇精品久久久久久久| tocl精华| 精品少妇内射三级| 国产91精品成人一区二区三区 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成电影免费在线| 夜夜骑夜夜射夜夜干| 亚洲精品美女久久av网站| 一区二区av电影网| 男男h啪啪无遮挡| 99精品欧美一区二区三区四区| 美女主播在线视频| 日本黄色视频三级网站网址 | 国产精品亚洲一级av第二区| 精品免费久久久久久久清纯 | 久久毛片免费看一区二区三区| 大香蕉久久成人网| 老司机亚洲免费影院| 精品国产一区二区三区四区第35| 国产成人一区二区三区免费视频网站| 欧美乱码精品一区二区三区| 国产午夜精品久久久久久| 中文字幕色久视频| 欧美成人午夜精品| 国产精品偷伦视频观看了| 亚洲综合色网址| av线在线观看网站| 亚洲精华国产精华精| 91九色精品人成在线观看| 久久亚洲精品不卡| 一区在线观看完整版| 久久国产精品人妻蜜桃| 在线观看免费视频日本深夜| 人妻一区二区av| 18禁黄网站禁片午夜丰满| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 老司机福利观看| 国产精品久久久久久精品古装| 电影成人av| 午夜福利视频精品| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| 欧美日韩一级在线毛片| 99精品在免费线老司机午夜| 美女视频免费永久观看网站| 久久av网站| 嫁个100分男人电影在线观看| 精品高清国产在线一区| 在线观看人妻少妇| 一区在线观看完整版| 热re99久久国产66热| 少妇精品久久久久久久| 日韩大片免费观看网站| 日韩中文字幕欧美一区二区| 精品免费久久久久久久清纯 | 国产一区二区三区在线臀色熟女 | 亚洲午夜精品一区,二区,三区| 男女无遮挡免费网站观看| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美亚洲二区| 国产男女超爽视频在线观看| 真人做人爱边吃奶动态| 啦啦啦免费观看视频1| 亚洲国产av新网站| 国产精品一区二区精品视频观看| 成在线人永久免费视频| svipshipincom国产片| 欧美精品高潮呻吟av久久| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 亚洲色图av天堂| 精品免费久久久久久久清纯 | 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| av免费在线观看网站| 国产三级黄色录像| 亚洲国产中文字幕在线视频| 无遮挡黄片免费观看| 亚洲五月婷婷丁香| 咕卡用的链子| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 国产亚洲精品第一综合不卡| 国产精品秋霞免费鲁丝片| 97人妻天天添夜夜摸| 人人妻人人澡人人爽人人夜夜| 日本av手机在线免费观看| 国产一区有黄有色的免费视频| 日本av免费视频播放| 高清欧美精品videossex| 国产在线视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 欧美激情 高清一区二区三区| 法律面前人人平等表现在哪些方面| 在线观看免费视频网站a站| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 精品国产一区二区三区久久久樱花|