• <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
    丰满饥渴人妻一区二区三| 日本av手机在线免费观看| 五月开心婷婷网| 欧美bdsm另类| 高清不卡的av网站| 晚上一个人看的免费电影| 国产欧美亚洲国产| 最新中文字幕久久久久| 啦啦啦在线观看免费高清www| 一本一本综合久久| 日本欧美视频一区| 2022亚洲国产成人精品| 亚洲,一卡二卡三卡| 国产 精品1| 男的添女的下面高潮视频| 午夜精品国产一区二区电影| 如日韩欧美国产精品一区二区三区 | 免费观看av网站的网址| 午夜久久久在线观看| 在线观看www视频免费| 老司机影院毛片| 久久人人爽人人爽人人片va| 精品一区二区三区视频在线| 性高湖久久久久久久久免费观看| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 女人精品久久久久毛片| 成人综合一区亚洲| 午夜激情久久久久久久| www.av在线官网国产| 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 国产不卡av网站在线观看| 亚洲国产av新网站| 欧美日韩av久久| 看免费成人av毛片| 免费不卡的大黄色大毛片视频在线观看| 国产黄色视频一区二区在线观看| 性高湖久久久久久久久免费观看| av福利片在线| a 毛片基地| 国产免费现黄频在线看| 777米奇影视久久| 国内精品宾馆在线| 亚洲国产日韩一区二区| 亚洲精品日韩av片在线观看| 女性被躁到高潮视频| 十八禁网站网址无遮挡| 777米奇影视久久| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 在线观看国产h片| 97在线人人人人妻| 男的添女的下面高潮视频| 精品少妇久久久久久888优播| 中文字幕制服av| 两个人的视频大全免费| 婷婷成人精品国产| 久久精品夜色国产| 欧美老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 精品人妻熟女av久视频| 黑丝袜美女国产一区| 国产精品免费大片| 97精品久久久久久久久久精品| 国产不卡av网站在线观看| av又黄又爽大尺度在线免费看| 国产精品一二三区在线看| 大陆偷拍与自拍| 亚洲欧美成人综合另类久久久| 精品一区在线观看国产| 久久99热6这里只有精品| 精品国产一区二区三区久久久樱花| 亚洲欧洲日产国产| 黑人巨大精品欧美一区二区蜜桃 | 午夜免费观看性视频| 草草在线视频免费看| 国产黄色视频一区二区在线观看| 久久久久久久久久人人人人人人| 一区二区三区精品91| 2021少妇久久久久久久久久久| kizo精华| 成人国产av品久久久| 丰满乱子伦码专区| 国产欧美日韩一区二区三区在线 | 免费大片黄手机在线观看| √禁漫天堂资源中文www| 国精品久久久久久国模美| av视频免费观看在线观看| 99视频精品全部免费 在线| av福利片在线| 寂寞人妻少妇视频99o| 国产欧美日韩一区二区三区在线 | 欧美激情国产日韩精品一区| 午夜激情av网站| 欧美激情国产日韩精品一区| av一本久久久久| 久久国产精品男人的天堂亚洲 | 久久久久精品性色| 最近中文字幕2019免费版| 久久 成人 亚洲| 人人妻人人爽人人添夜夜欢视频| 亚洲少妇的诱惑av| 少妇猛男粗大的猛烈进出视频| 久久热精品热| 九九爱精品视频在线观看| 大码成人一级视频| 欧美亚洲日本最大视频资源| 亚洲av日韩在线播放| 婷婷色麻豆天堂久久| 国产精品 国内视频| 中文字幕av电影在线播放| 97超碰精品成人国产| 精品久久国产蜜桃| 在线观看www视频免费| 国产免费又黄又爽又色| 18禁裸乳无遮挡动漫免费视频| 成人影院久久| 国产精品久久久久久精品古装| 亚洲激情五月婷婷啪啪| 少妇人妻精品综合一区二区| 亚洲av日韩在线播放| 视频中文字幕在线观看| 日本欧美视频一区| 成人国语在线视频| 乱码一卡2卡4卡精品| 满18在线观看网站| 黄色毛片三级朝国网站| 五月开心婷婷网| 视频中文字幕在线观看| 精品人妻熟女毛片av久久网站| 26uuu在线亚洲综合色| 久久久久久久精品精品| 中文字幕制服av| 美女主播在线视频| 国产精品蜜桃在线观看| 成人黄色视频免费在线看| 久久99热6这里只有精品| 亚洲成色77777| 亚洲精品av麻豆狂野| 国产亚洲最大av| 80岁老熟妇乱子伦牲交| 国产精品 国内视频| 三级国产精品片| 亚洲无线观看免费| 九色成人免费人妻av| 国产欧美日韩一区二区三区在线 | 成年女人在线观看亚洲视频| 色5月婷婷丁香| 成人影院久久| 欧美激情国产日韩精品一区| 久久午夜综合久久蜜桃| 婷婷色av中文字幕| 国产黄频视频在线观看| 女性被躁到高潮视频| 精品人妻在线不人妻| 91成人精品电影| 日韩成人伦理影院| 在线观看免费视频网站a站| 春色校园在线视频观看| 久久热精品热| 超碰97精品在线观看| 精品国产乱码久久久久久小说| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 国产精品99久久99久久久不卡 | 亚洲美女搞黄在线观看| 欧美日韩视频精品一区| 内地一区二区视频在线| 久久久亚洲精品成人影院| 大码成人一级视频| 天美传媒精品一区二区| 成人亚洲精品一区在线观看| 亚洲国产欧美在线一区| 久久久久精品性色| 边亲边吃奶的免费视频| 午夜激情福利司机影院| 黑丝袜美女国产一区| 韩国av在线不卡| 成年美女黄网站色视频大全免费 | 亚洲av福利一区| 黄色视频在线播放观看不卡| 日本色播在线视频| 亚洲婷婷狠狠爱综合网| 黑丝袜美女国产一区| 99九九线精品视频在线观看视频| 97精品久久久久久久久久精品| 中文字幕人妻丝袜制服| 久久久久网色| 午夜日本视频在线| 欧美人与善性xxx| 高清视频免费观看一区二区| 91久久精品国产一区二区三区| 久久国产精品大桥未久av| 中文乱码字字幕精品一区二区三区| 狂野欧美激情性xxxx在线观看| 观看av在线不卡| 午夜福利视频在线观看免费| 欧美性感艳星| 肉色欧美久久久久久久蜜桃| 秋霞在线观看毛片| 日韩,欧美,国产一区二区三区| 卡戴珊不雅视频在线播放| 亚洲国产精品一区二区三区在线| 嘟嘟电影网在线观看| √禁漫天堂资源中文www| 老司机影院成人| 亚洲国产精品999| 精品久久国产蜜桃| 中文字幕av电影在线播放| 高清黄色对白视频在线免费看| 超碰97精品在线观看| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 国产成人91sexporn| 国产探花极品一区二区| 丰满乱子伦码专区| 夜夜看夜夜爽夜夜摸| 国产黄色视频一区二区在线观看| 欧美日韩视频精品一区| 亚洲国产欧美日韩在线播放| 99久久中文字幕三级久久日本| 秋霞伦理黄片| 免费观看无遮挡的男女| 性高湖久久久久久久久免费观看| 亚洲欧美日韩另类电影网站| 99久久中文字幕三级久久日本| 在线观看免费高清a一片| 国产精品 国内视频| 国产深夜福利视频在线观看| av播播在线观看一区| 成人二区视频| 日韩av免费高清视频| 国产在线视频一区二区| 日韩亚洲欧美综合| 久久女婷五月综合色啪小说| 麻豆精品久久久久久蜜桃| 国产成人av激情在线播放 | 永久免费av网站大全| 自线自在国产av| 欧美三级亚洲精品| 丰满饥渴人妻一区二区三| 免费久久久久久久精品成人欧美视频 | 涩涩av久久男人的天堂| 少妇人妻 视频| 丝瓜视频免费看黄片| 综合色丁香网| 国产欧美亚洲国产| 一本大道久久a久久精品| 亚洲成人av在线免费| 美女xxoo啪啪120秒动态图| 日韩av不卡免费在线播放| 丝瓜视频免费看黄片| 久久97久久精品| 婷婷色综合大香蕉| 新久久久久国产一级毛片| 中国国产av一级| 2022亚洲国产成人精品| 午夜视频国产福利| 草草在线视频免费看| 中国美白少妇内射xxxbb| 熟女人妻精品中文字幕| 亚洲精品色激情综合| 亚洲成色77777| 特大巨黑吊av在线直播| 午夜免费鲁丝| 大香蕉久久成人网| 熟女av电影| av国产精品久久久久影院| 少妇猛男粗大的猛烈进出视频| 亚洲av电影在线观看一区二区三区| 日本vs欧美在线观看视频| 肉色欧美久久久久久久蜜桃| 亚洲精品456在线播放app| av免费在线看不卡| 欧美+日韩+精品| 国产老妇伦熟女老妇高清| 亚洲一级一片aⅴ在线观看| 最近中文字幕2019免费版| 亚洲国产精品国产精品| 午夜免费鲁丝| 免费高清在线观看视频在线观看| 18+在线观看网站| 国产 精品1| 夫妻午夜视频| 最黄视频免费看| kizo精华| 国产黄片视频在线免费观看| 26uuu在线亚洲综合色| 成人亚洲精品一区在线观看| 精品人妻一区二区三区麻豆| 成人无遮挡网站| 亚洲精品乱久久久久久| 天天操日日干夜夜撸| 啦啦啦啦在线视频资源| 美女福利国产在线| 99热网站在线观看| 99精国产麻豆久久婷婷| 老女人水多毛片| 一本大道久久a久久精品| 免费久久久久久久精品成人欧美视频 | av不卡在线播放| 多毛熟女@视频| 国产免费现黄频在线看| 中文精品一卡2卡3卡4更新| av专区在线播放| 91精品国产国语对白视频| 国产精品无大码| 欧美激情国产日韩精品一区| 亚洲欧美中文字幕日韩二区| 男女高潮啪啪啪动态图| 伦精品一区二区三区| 一边摸一边做爽爽视频免费| 视频中文字幕在线观看| 成人毛片60女人毛片免费| h视频一区二区三区| 久久午夜综合久久蜜桃| 日韩一区二区视频免费看| 亚洲一区二区三区欧美精品| 免费黄网站久久成人精品| 午夜福利影视在线免费观看| 久久人人爽人人爽人人片va| 亚洲国产精品999| 免费看av在线观看网站| 99九九线精品视频在线观看视频| .国产精品久久| 亚洲图色成人| 国产精品熟女久久久久浪| 国产欧美另类精品又又久久亚洲欧美| av.在线天堂| 国产精品人妻久久久久久| 曰老女人黄片| 9色porny在线观看| 亚洲欧洲国产日韩| 大话2 男鬼变身卡| 欧美精品高潮呻吟av久久| 精品一区在线观看国产| 黑丝袜美女国产一区| 中文字幕久久专区| 久久人人爽av亚洲精品天堂| 亚洲av.av天堂| 一级毛片黄色毛片免费观看视频| 欧美老熟妇乱子伦牲交| 亚洲精品国产av成人精品| 99热网站在线观看| 久久久亚洲精品成人影院| 久久久久久久久久久丰满| 中国国产av一级| 亚洲精品自拍成人| 色5月婷婷丁香| 成年av动漫网址| a级毛片免费高清观看在线播放| 亚洲精品视频女| 国产精品久久久久久久久免| 免费看不卡的av| 亚洲av福利一区| 午夜福利网站1000一区二区三区| 日本-黄色视频高清免费观看| 男女免费视频国产| 亚洲av日韩在线播放| 亚洲av欧美aⅴ国产| 一级毛片黄色毛片免费观看视频| 国产免费福利视频在线观看| 成年美女黄网站色视频大全免费 | 亚洲精品乱码久久久久久按摩| 日韩av免费高清视频| 亚洲精品国产av蜜桃| 成人亚洲精品一区在线观看| 国产一区二区在线观看av| 亚洲欧美一区二区三区国产| 一级毛片aaaaaa免费看小| 国产欧美日韩综合在线一区二区| 美女福利国产在线| 男女无遮挡免费网站观看| 插逼视频在线观看| 久久ye,这里只有精品| 妹子高潮喷水视频| 日日啪夜夜爽| 久久久久久久亚洲中文字幕| 波野结衣二区三区在线| 少妇的逼好多水| 水蜜桃什么品种好| 国产在线免费精品| 日韩制服骚丝袜av| 最近中文字幕高清免费大全6| 国产黄色视频一区二区在线观看| 黄色欧美视频在线观看| 看十八女毛片水多多多| 91午夜精品亚洲一区二区三区| 久久久久久久久大av| 校园人妻丝袜中文字幕| 男人爽女人下面视频在线观看| av在线app专区| 亚洲在久久综合| 日韩一本色道免费dvd| 视频中文字幕在线观看| 亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区 | 黑丝袜美女国产一区| 免费不卡的大黄色大毛片视频在线观看| 亚洲高清免费不卡视频| 久久久久久久大尺度免费视频| 成年女人在线观看亚洲视频| 最近中文字幕2019免费版| 最近中文字幕高清免费大全6| 欧美 日韩 精品 国产| 久久久亚洲精品成人影院| 久久ye,这里只有精品| 久久精品国产鲁丝片午夜精品| 日本黄大片高清| 久久精品久久精品一区二区三区| 亚洲av在线观看美女高潮| 午夜日本视频在线| 久久综合国产亚洲精品| 丝袜脚勾引网站| 99视频精品全部免费 在线| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成色77777| 午夜福利,免费看| 国产精品国产三级专区第一集| 久久国产亚洲av麻豆专区| 精品亚洲成国产av| 国产一区二区在线观看日韩| 一本大道久久a久久精品| 欧美人与善性xxx| 亚州av有码| 日日啪夜夜爽| 欧美日韩一区二区视频在线观看视频在线| 五月伊人婷婷丁香| 久久 成人 亚洲| 丁香六月天网| 日韩一本色道免费dvd| av天堂久久9| 美女国产高潮福利片在线看| 蜜桃久久精品国产亚洲av| 高清视频免费观看一区二区| 婷婷色麻豆天堂久久| 视频区图区小说| 国产精品免费大片| 国产成人a∨麻豆精品| 国产精品国产三级国产专区5o| 精品视频人人做人人爽| 午夜日本视频在线| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 日本黄大片高清| 一本—道久久a久久精品蜜桃钙片| xxxhd国产人妻xxx| 国产精品久久久久久精品古装| 国产精品99久久久久久久久| 国产精品无大码| 国产精品欧美亚洲77777| xxx大片免费视频| 熟女人妻精品中文字幕| 两个人免费观看高清视频| 欧美精品人与动牲交sv欧美| 欧美bdsm另类| 国产高清有码在线观看视频| 中文字幕亚洲精品专区| 日韩精品免费视频一区二区三区 | 国产高清有码在线观看视频| 国产无遮挡羞羞视频在线观看| 大香蕉久久成人网| 亚洲第一av免费看| 你懂的网址亚洲精品在线观看| 国产伦精品一区二区三区视频9| 在线观看美女被高潮喷水网站| 日韩一区二区三区影片| 亚洲国产精品999| 久久精品国产自在天天线| 卡戴珊不雅视频在线播放| 精品亚洲成a人片在线观看| 国产女主播在线喷水免费视频网站| 免费久久久久久久精品成人欧美视频 | 男男h啪啪无遮挡| 丁香六月天网| 亚洲欧美成人精品一区二区| 婷婷色麻豆天堂久久| 视频中文字幕在线观看| 亚洲av不卡在线观看| 黑人巨大精品欧美一区二区蜜桃 | 日韩欧美一区视频在线观看| 国产精品 国内视频| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 熟女电影av网| 国产精品一区二区在线不卡| a级毛色黄片| 精品视频人人做人人爽| 久久精品国产亚洲av涩爱| 亚洲精品,欧美精品| 亚洲色图综合在线观看| 国产成人aa在线观看| 国产高清三级在线| 免费观看在线日韩| 国产成人91sexporn| 日韩成人av中文字幕在线观看| 国产深夜福利视频在线观看| 久久久久久久久久人人人人人人| 大码成人一级视频| 性色av一级| 纵有疾风起免费观看全集完整版| 美女cb高潮喷水在线观看| 日韩熟女老妇一区二区性免费视频| 精品国产国语对白av| 99re6热这里在线精品视频| 国产成人91sexporn| 99热全是精品| 欧美人与善性xxx| 又粗又硬又长又爽又黄的视频| 三级国产精品片| 欧美最新免费一区二区三区| 性色avwww在线观看| 成人毛片60女人毛片免费| 麻豆乱淫一区二区| 99热国产这里只有精品6| 亚洲精品一区蜜桃| 亚洲内射少妇av| 国产又色又爽无遮挡免| 久久久久久久精品精品| 夜夜看夜夜爽夜夜摸| 人妻少妇偷人精品九色| 欧美日韩视频精品一区| 国产日韩欧美亚洲二区| 黑人高潮一二区| 蜜桃国产av成人99| 91国产中文字幕| 男男h啪啪无遮挡| 曰老女人黄片| 交换朋友夫妻互换小说| 日韩亚洲欧美综合| 在线天堂最新版资源| 丰满少妇做爰视频| 岛国毛片在线播放| 国产欧美另类精品又又久久亚洲欧美| 2021少妇久久久久久久久久久| 热99久久久久精品小说推荐| 欧美人与性动交α欧美精品济南到 | 熟女人妻精品中文字幕| 秋霞伦理黄片| 在线观看免费视频网站a站| 日本与韩国留学比较| 久久人人爽av亚洲精品天堂| 亚洲精品国产av蜜桃| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | tube8黄色片| 久久久久久久久久久免费av| 视频区图区小说| 极品人妻少妇av视频| 久久久久网色| 国产日韩欧美亚洲二区| 大话2 男鬼变身卡| 亚洲国产av新网站| 草草在线视频免费看| 午夜福利影视在线免费观看| 久久久亚洲精品成人影院| 亚洲美女视频黄频| 国产精品久久久久成人av| 美女内射精品一级片tv| 国产在视频线精品| 亚洲性久久影院| 黄片播放在线免费| 亚洲国产精品国产精品| 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 久久人妻熟女aⅴ| 亚洲欧美日韩卡通动漫| 亚洲三级黄色毛片| 国产精品无大码| 国产老妇伦熟女老妇高清| 国产精品国产av在线观看| 91aial.com中文字幕在线观看| 伦理电影免费视频| 欧美日韩亚洲高清精品| 全区人妻精品视频| 人人澡人人妻人| av线在线观看网站| 极品人妻少妇av视频| 亚洲欧美日韩另类电影网站| 欧美激情极品国产一区二区三区 | 韩国高清视频一区二区三区| 国产欧美日韩综合在线一区二区| 特大巨黑吊av在线直播| 国产精品人妻久久久影院| a级毛片黄视频| 看非洲黑人一级黄片| 美女视频免费永久观看网站| 亚洲精品乱久久久久久| 又黄又爽又刺激的免费视频.| 亚洲伊人久久精品综合| 久久99热6这里只有精品| 精品久久久久久电影网| 久久久久久久久久久丰满| av播播在线观看一区| 免费黄频网站在线观看国产| 七月丁香在线播放| 亚洲欧美色中文字幕在线| 天美传媒精品一区二区| 日韩熟女老妇一区二区性免费视频| 热re99久久国产66热| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡 | 午夜激情福利司机影院| 3wmmmm亚洲av在线观看| 成年av动漫网址| 国产片内射在线| 黄色欧美视频在线观看| 丝袜美足系列| 免费高清在线观看视频在线观看| 免费黄频网站在线观看国产| √禁漫天堂资源中文www| 亚洲人成77777在线视频| 另类精品久久| 蜜臀久久99精品久久宅男|