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

    考慮抗剪強度參數(shù)空間變異性的邊坡穩(wěn)定性分析

    2016-10-13 19:40:39袁葳常曉林段寅程勇剛馬剛
    中南大學學報(自然科學版) 2016年11期
    關鍵詞:滑面黏聚力摩擦角

    袁葳,常曉林,段寅, 3,程勇剛,馬剛

    ?

    考慮抗剪強度參數(shù)空間變異性的邊坡穩(wěn)定性分析

    袁葳1, 2,常曉林1,段寅1, 3,程勇剛1,馬剛1

    (1. 武漢大學水資源與水電工程科學國家重點實驗室,湖北武漢,430072;2. 湖北省水利水電規(guī)劃勘測設計院,湖北武漢,430064;3.長江勘測規(guī)劃設計研究院,湖北武漢,430010)

    提出基于商業(yè)軟件ABAQUS的邊坡穩(wěn)定分析的隨機有限元法。采用Karhunen-Loeve級數(shù)展開方法建立土體抗剪強度參數(shù)的隨機場模型,在強度折減有限元法的基礎上模擬邊坡失穩(wěn)過程,并采用塑性屈服區(qū)貫通以及計算不收斂判據(jù)作為邊坡失穩(wěn)判據(jù),利用python腳本文件分析最危險滑面分布規(guī)律。運用該方法研究土體內(nèi)摩擦角和黏聚力的空間變異性對邊坡穩(wěn)定性的影響。研究結果表明:所提的邊坡穩(wěn)定分析的隨機有限元法操作簡單易于實現(xiàn),能夠從安全系數(shù)以及最危險滑面2個方面分析土體抗剪強度參數(shù)的空間變異性對邊坡穩(wěn)定性的影響;內(nèi)摩擦角隨機場中相關長度的改變對邊坡的穩(wěn)定性影響不明顯;而黏聚力的空間變異性對邊坡穩(wěn)定性影響則相對更大,黏聚力隨機場中水平與豎直相關長度的比值不同,邊坡穩(wěn)定性的變化規(guī)律不同。此外,內(nèi)摩擦角與黏聚力的相關系數(shù)越大,失效概率越高,最危險滑面分布越廣。

    空間變異性;邊坡穩(wěn)定性;強度折減有限元法;隨機有限元方法

    在巖體工程領域,傳統(tǒng)的考慮邊坡穩(wěn)定性的方法是使用安全系數(shù)評判其風險等級。但是,安全系數(shù)并非評判邊坡風險等級的唯一標準,即使邊坡具有相同的安全系數(shù),由于土體參數(shù)的不均勻性,其風險等級可能是不同的[1]。因此,近年來眾多學者系統(tǒng)地研究了邊坡的可靠度分析方法。GRIFFITH等[2?3]提出了基于強度折減法和隨機場理論的隨機有限元方法,對比了采用隨機變量模型和隨機場模型分析邊坡可靠度的成果。SHEN等[4]結合彈塑性有限差分法和隨機場理論分析非均質(zhì)邊坡的可靠度,采用局部平均法模擬黏聚力參數(shù)的空間變異性。CHO[5]運用隨機極限平衡法研究土體抗剪強度的空間變異性對邊坡穩(wěn)定性的影響,采用Karhunen-Loeve級數(shù)展開方法模擬土體參數(shù)的空間變異性。譚曉慧等[6]基于一階可靠度方法,對邊坡的整體可靠度指標進行了計算。吳振君等[7]采用驗算點法和一維隨機場模擬實現(xiàn)了考慮地質(zhì)成因的土坡可靠度分析。李典慶等[8?11]提出了非侵入式隨機有限元方法,研究了該方法在邊坡穩(wěn)定、地下洞室可靠度分析中的應用。國內(nèi)外學者進行邊坡可靠度分析時多采用傳統(tǒng)極限平衡法(LEM)[1, 5?6],LEM方法計算時間短,適用于蒙特卡洛模擬,但是,采用LEM計算時需要假定臨界滑面,通?;谔囟ǖ幕嬗嬎氵吰率Ц怕省EcLEM相比,彈塑性有限元分析方法無需假定臨界滑面的位置,能夠求得任意形狀的滑面及相應的折減系數(shù),更為真實地反映邊坡的穩(wěn)定性。GRIFFITH等[2?3]針對隨機有限元方法進行了大量的研究,采用局部平均法模擬參數(shù)空間變異性,依據(jù)迭代不收斂判據(jù)得到邊坡的安全系數(shù)。然而,采用局部平均法時,隨機變量的數(shù)目與有限元單位的數(shù)目相同,計算效率不高。肖特等[12]提出基于有限元強度折減法的非侵入式邊坡可靠度分析方法,采用位移突變作為邊坡失穩(wěn)的判據(jù),應用于多層土質(zhì)邊坡可靠度分析。上述基于彈塑性有限元法的邊坡可靠度研究往往局限于安全系數(shù)、失效概率等指標的規(guī)律,針對滑面的分布規(guī)律的研究偏少,然而,潛在的滑坡的規(guī)模與邊坡最危險滑面的分布直接相關。因此,有必要從邊坡安全系數(shù)與最危險滑面分布2方面研究巖土參數(shù)空間變異性對邊坡穩(wěn)定性的影響。本文作者提出了基于有限元軟件ABAQUS的邊坡穩(wěn)定分析的隨機有限元法。該方法采用Karhunen-Loeve級數(shù)展開方法建立土體參數(shù)的隨機場模型;基于強度折減有限元法模擬邊坡失穩(wěn)過程;主要以塑性屈服區(qū)貫通作為邊坡失穩(wěn)判據(jù);利用python腳本文件對計算結果進行后處理以得出最危險滑面分布規(guī)律。然后,在蒙特卡洛隨機有限元法的基礎上,從安全系數(shù)和最危險滑面分布2個方面,研究了土體內(nèi)摩擦角和黏聚力的空間變異性對邊坡穩(wěn)定性的影響。

    1 土體空間變異性的隨機場模擬

    1.1 Karhunen-Loeve級數(shù)展開方法

    采用Karhunen-Loeve級數(shù)展開方法離散隨機場實質(zhì)上將土體參數(shù)隨機場的離散轉(zhuǎn)化為求解第2類Fredhom積分方程的特征值問題,即

    式中:1和2為隨機場區(qū)域內(nèi)任意兩點的坐標;為隨機場區(qū)域內(nèi)任意兩點的特性值之間的相關函數(shù)值;和分別為與相關函數(shù)對應的特征值以及特征函數(shù)。

    根據(jù)求解Fredhom積分方程的方法的不同,K-L級數(shù)展開方法可以分為解析法和數(shù)值法。解析方法在具體工程應用中具有局限性,部分隨機場的相關函數(shù)的特征值以及特征函數(shù)用解析方法求解困難,甚至不可解。LI等[1]研究了不同相關函數(shù)對邊坡可靠度指標的影響,研究表明影響不大。因此,本文采用的相關函數(shù)為指數(shù)型式,并采用解析方法求解Fredhom積分方程,以減小程序的復雜性。指數(shù)型相關函數(shù)如下:

    式中:和為隨機場區(qū)域內(nèi)任意點的坐標;h和v分別為水平向和豎直向相關長度。

    當相關函數(shù)是可分離時,二維隨機場相關函數(shù)的特征值可簡化為一維相關函數(shù)特征值的乘積。以一維標準高斯隨機場為例,詳細說明K-L級數(shù)展開方法的具體過程。假定,特征值和特征函數(shù)滿足第2類Fredhom積分公式:

    其中相關函數(shù)為

    (4)

    將式(4)代入式(3)得

    (6)

    由式(6)解得

    當為偶數(shù)時,由式(5)可得

    (8)

    由式(8)解得

    基于以上計算過程得到的相關函數(shù)的特征解,可將標準高斯隨機場展開為

    (10)

    本文假定土體的內(nèi)摩擦角和黏聚力服從Weibull分布,其概率密度函數(shù)、累積分布函數(shù)、均值、標準差分別為:

    (12)

    (13)

    式中:為均質(zhì)度,越大,參數(shù)分布越均勻。

    將土體參數(shù)在計算區(qū)域內(nèi)模擬為二維Weibull分布的統(tǒng)計均質(zhì)隨機場,通過等概率變換,Weibull隨機場可以表示為

    在計算精度滿足要求的前提下,通常截取前項K-L級數(shù)展開以提高計算效率,式(15)截取前項后可表示為

    1.2 抗剪強度參數(shù)隨機場的模擬

    本文模擬土體參數(shù)的隨機分布分為3種:1) 考慮內(nèi)摩擦角的空間變異性;2) 考慮黏聚力的空間變異性;3)同時考慮兩者的空間變異性。以上3種情況均采用Weibull隨機場表征參數(shù)的空間變異性。隨機場的模擬采用1.1節(jié)介紹的K-L級數(shù)展開方法。邊坡的幾何尺寸以及有限元網(wǎng)格見3.1節(jié)。

    圖1所示為土體內(nèi)摩擦角隨機場在2種均質(zhì)度下的模擬結果,內(nèi)摩擦角的均值均為15°。由圖1可以看出:當=2時,內(nèi)摩擦角取值分散;而當=5時,內(nèi)摩擦角取值集中。

    (a) mφ=2, lφ,h=20 m, lφ,v=2 m; (b) mφ=5, lφ,h=20 m, lφ,v=2 m

    圖2所示為土體黏聚力隨機場在4種不同水平向和豎直向相關長度組合下的模擬結果,黏聚力的均值均為20 kPa。在均質(zhì)度相同的情況下,當時,黏聚力隨機場中高黏聚力帶和低黏聚力帶呈水平向?qū)訝罘植?;當時,呈塊狀分布;相關長度越長,高黏聚力帶、低黏聚力帶越連續(xù)。

    模擬具有一定相關性的內(nèi)摩擦角和黏聚力隨機場時,通過控制隨機向量和之間的相關性進而得到具有相同相關性的內(nèi)摩擦角和黏聚力隨機場。圖3給出了土體內(nèi)摩擦角、黏聚力隨機場在2種不同相關系數(shù)下的模擬結果。其中內(nèi)摩擦角和黏聚力均服從=3的Weibull分布,均值分別為15°和20 kPa,水平相關長度取20 m,豎直相關長度取2 m。當時,內(nèi)摩擦角隨機場與黏聚力隨機場分布規(guī)律大致相反,即內(nèi)摩擦角大的區(qū)域其黏聚力較低;當時,內(nèi)摩擦角隨機場與黏聚力隨機場分布規(guī)律相似。

    (a) mc=3, lc,h=2 m, lc,v=2 m; (b) mc=3, lc,h=20 m, lc,v=20 m; (c) mc=3, lc,h=10 m, lc,v=1 m; (d) mc=3, lc,h=40 m, lc,v=4 m

    (a), (b) r(c,φ)=?0.5; (c), (d) r(c,φ)=0.5

    2 邊坡強度折減有限元法

    ZIENKIEWICS等[13]于1975年提出強度折減有限元法,經(jīng)過幾十年的發(fā)展,強度折減有限元法在邊坡穩(wěn)定分析方面取得了眾多成果,GRIFFITH等[14]運用強度折減非線性有限元法對不同類型的邊坡進行了研究,驗證了該方法的可靠性;DAWSON等[15]討論了網(wǎng)格密度、單元類型、邊界范圍、屈服準則等對計算精度的影響;欒茂田等[16]、劉金龍等[17]、呂慶等[18]對邊坡失穩(wěn)判據(jù)進行了研究。邊坡失穩(wěn)判據(jù)通常分為3類:計算不收斂、塑性屈服區(qū)貫通以及特征點位移突變??紤]到邊坡中不同特征點的位移曲線存在一定的差異,且曲線的突變點不易確定,本文主要采用計算不收斂判據(jù)以及塑性屈服區(qū)貫通判據(jù)。

    本文依托有限元軟件ABAQUS,利用其用戶子程序接口USDFLD實現(xiàn)場變量與時間步長、材料參數(shù)的對應關系,在每個分析步中,材料強度參數(shù)在給定的范圍內(nèi)逐漸折減,利用ABAQUS/Standard模塊模擬邊坡破壞過程。

    圖4給出了編寫的用戶子程序USDFLD的主要計算流程,其中表示計算時間,NOEL表示單元編號,K表示時刻對應的折減系數(shù)。單元的初始強度參數(shù)符合Weibull分布,實現(xiàn)Weibull隨機場的程序嵌入USDFLD用戶子程序中,僅在總時間為0時被調(diào)用。本文通過2個(0,1)內(nèi)均勻分布的獨立隨機數(shù)轉(zhuǎn)換得到標準正態(tài)隨機數(shù),其中均勻分布隨機數(shù)是由fortran自帶的函數(shù)random_number生成。

    圖4 用戶子程序USDFLD計算流程

    3 邊坡穩(wěn)定的確定性分析

    3.1 計算模型

    本文采用的算例為一均質(zhì)土坡,尺寸如圖5所示,坡高為10 m,坡比為1:2。計算參數(shù)如下:黏聚力20 kPa,內(nèi)摩擦角15°,重度20 kN/m3,彈性模量8 MPa,泊松比0.25。

    在ABAQUS中按照平面應變問題建立二維邊坡有限元分析模型,如圖6所示,共剖分1 200個單元和1 281個節(jié)點。土體采用理想彈塑性本構模型和Mohr-Coulomb屈服準則。計算時模型底部采取固定約束,兩側采取法向約束;模型施加自重荷載。

    圖5 邊坡幾何尺寸

    圖6 邊坡有限元模型

    3.2 確定性分析

    本節(jié)計算未考慮參數(shù)的空間變異性,每個單元賦予相同的和。計算時設定的強度折減系數(shù)的變化范圍為0至5.0,每一個分析步內(nèi)S增加0.1。圖7所示為邊坡S=1.54和1.55時的等效塑性應變云圖,從圖7可以看出:當S=1.54時,坡腳開始出現(xiàn)塑性屈服區(qū),隨著的增加,塑性屈服區(qū)逐漸向坡頂擴展;當S=1.55時,塑性屈服區(qū)完全貫通,形成了從坡頂至坡腳的連續(xù)滑動面;當折減系數(shù)為1.77時,計算不收斂。

    (a) FS=1.540; (b) FS=1.550

    為了對比分析強度折減有限元法和傳統(tǒng)極限平衡法的計算結果,本文采用加拿大邊坡分析軟件SLOPE/W按照Morgenstern-Price法計算得到的該邊坡的安全系數(shù)為1.533。確定性分析時,根據(jù)計算不收斂得到的邊坡安全系數(shù)大于依據(jù)塑性區(qū)貫通判據(jù)得到的安全系數(shù),且后者得到的安全系數(shù)與極限平衡法計算得到的安全系數(shù)十分接近,相對誤差僅為1.11%,二者計算得到的滑面位置基本一致。因此,本文確定邊坡安全系數(shù)時以塑性屈服區(qū)貫通判據(jù)為主,以計算不收斂判據(jù)為輔。

    4 邊坡穩(wěn)定的隨機分析

    進行隨機分析時,計算結果數(shù)量大,重復性工作多,為了高效地判斷邊坡的安全系數(shù)以及確定最危險滑面的位置,編寫python腳本對結果文件進行后處理,以得出最危險滑面分布規(guī)律。本文確定最危險滑面的方法與文獻[19]中所描述的方法類似,即將水平坐標相同的單元分為一組,提取每組中等效塑性應變最大的單元的形心坐標,當計算到某一步時,提取的所有單元形心能夠連成完整的滑面,此時對應的折減系數(shù)即為邊坡的安全系數(shù),所連的滑面即為最危險滑面。以第3節(jié)確定性計算結果為例驗證所提后處理方法的可行性,運行python腳本文件,可得塑性屈服區(qū)貫通時對應的安全系數(shù)為1.55,計算不收斂時對應的安全系數(shù)為1.77,最危險滑面如圖8所示,與圖7(b)中貫通的等效塑性應變帶位置一致。

    圖8 最危險滑面

    為了更全面地分析抗剪強度參數(shù)的空間變異性對邊坡穩(wěn)定的影響,定義失效概率為

    式中:f為安全系數(shù)小于1.0的次數(shù);m為總共的隨機次數(shù)。

    4.1 內(nèi)摩擦角的空間變異性的影響

    本節(jié)計算只考慮內(nèi)摩擦角的空間變異性,計算模型與確定性分析中采用的模型一致,計算參數(shù)見表1,黏聚力均取20 kPa。2組試驗分別考慮均質(zhì)度和相關長度對邊坡穩(wěn)定性的影響。通常水平方向相關長度大于豎直方向相關長度[20],故考慮均質(zhì)度對邊坡穩(wěn)定性的影響時,取,h=20 m,,v=2 m。

    表1 考慮內(nèi)摩擦角空間變異性時的參數(shù)

    圖9所示為內(nèi)摩擦角=3,,h=20 m,,v=2 m對應的隨機分析結果。當隨機次數(shù)達到200以上時,平均安全系數(shù)呈降低的趨勢;邊坡失效概率持續(xù)為0,考慮到計算時間和計算精度,每組計算時隨機次數(shù)為200次。

    圖9 蒙特卡洛隨機分析結果(mφ=3,lφ,h=20 m,lφ,v=2 m)

    表2所示為內(nèi)摩擦角隨機場在不同均質(zhì)度下邊坡安全系數(shù)的統(tǒng)計參數(shù)。隨著均質(zhì)度的增加,平均安全系數(shù)由1.39提高到1.51,變異系數(shù)、失效概率降低。圖10所示為不同均質(zhì)度下最危險滑面分布范圍。對比=2和=5時滑面分布可以看出,均質(zhì)度越高,最危險滑面分布越集中。

    表2 考慮內(nèi)摩擦角,lφ,h=20 m,lφ,v=2 m時不同均質(zhì)度對應的安全系數(shù)分布

    圖10 考慮內(nèi)摩擦角,lφ,h=20 m,lφ,v=2 m時不同均質(zhì)度對應的最危險滑面分布

    表3所示為內(nèi)摩擦角隨著機場在不同相關長度下邊坡安全系數(shù)的統(tǒng)計參數(shù)。圖11所示為不同相關長度下最危險滑面分布范圍。當,h/,v=1時,隨著相關長度的增加,平均安全系數(shù)僅提高2%,變異系數(shù)提高,失效概率不超過0.01;最危險滑面分布更加集中。當,h/,v=10時,變異系數(shù)和失效概率的變化規(guī)律與,h/,v=1時一致,平均安全系數(shù)稍微降低;最危險滑面分布范圍擴大。

    表3 考慮內(nèi)摩擦角,mφ=3時不同相關長度對應的安全系數(shù)分布

    (a) lφ,h/lφ,v=1; (b) lφ,h/lφ,v=10

    對比內(nèi)摩擦角隨機場在不同均質(zhì)度與不同相關長度下的計算結果可知,邊坡的整體穩(wěn)定性對于內(nèi)摩擦角均質(zhì)度的敏感性高于相關長度,均質(zhì)度越高,平均安全系數(shù)與確定性計算結果越接近;相關長度的改變對邊坡的平均安全系數(shù)和失效概率影響不大。

    4.2 黏聚力的空間變異性的影響

    本節(jié)計算只考慮黏聚力的空間變異性,計算模型與確定性分析中采用的模型一致,計算參數(shù)見表4,內(nèi)摩擦角均取15°。

    表4 考慮黏聚力空間變異性時的參數(shù)取值

    圖12所示為黏聚力=3,,h=20 m,,v=2 m對應的隨機分析結果。由圖12可知:當隨機次數(shù)超過200次時,曲線趨于平緩,故本節(jié)每組計算時隨機次數(shù)定為200次。

    (a) 平均安全系數(shù);(b) 失效概率

    表5所示為黏聚力隨機場在不同均質(zhì)度下邊坡安全系數(shù)的統(tǒng)計參數(shù),平均安全系數(shù)、變異系數(shù)、失效概率的變化規(guī)律與考慮內(nèi)摩擦角的均質(zhì)度變化時得出的規(guī)律相同,但是變化幅度明顯提高,當=2時,平均安全系數(shù)為1.10,當=5時,平均安全系數(shù)增加到1.49,失效概率由0.40下降到0.01;在相同的均質(zhì)度下,與內(nèi)摩擦角的空間變異性相比,黏聚力的空間變異性對邊坡穩(wěn)定性的影響更大,平均安全系數(shù)更低,變異系數(shù)、失效概率更高。圖13所示為不同均質(zhì)度下最危險滑面分布范圍。隨著均質(zhì)度提高,最危險滑面分布逐漸集中。

    表5 考慮黏聚力,lc,h=20 m,lc,v=2 m時不同均質(zhì)度對應的安全系數(shù)分布

    圖13 考慮黏聚力,lc,h=20 m,lc,v=2 m時不同均質(zhì)度對應的最危險滑面分布

    表6所示為黏聚力隨機場在不同相關長度下邊坡安全系數(shù)的統(tǒng)計參數(shù),圖14所示為不同相關長度下最危險滑面分布范圍。當,h/,v=1時,隨著,h由2 m增加到20 m,平均安全系數(shù)相應提高,變異系數(shù)、失效概率降低,最危險滑面分布范圍縮?。欢?i>,h/,v=10時,隨著相關長度的增加,平均安全系數(shù)等的變化規(guī)律與,h/,v=1時相反。

    表6 考慮黏聚力,mc=3時不同相關長度對應的安全系數(shù)分布

    (a) lc,h/lc,v=1; (b) lc,h/lc,v=10

    與考慮內(nèi)摩擦角的空間變異性的計算結果相比,當黏聚力隨機場中均質(zhì)度與相關長度發(fā)生變化時,邊坡的平均安全系數(shù)、變異系數(shù)以及失效概率的變化幅度相對較高,可見,黏聚力的空間變異性對邊坡穩(wěn)定性影響大于內(nèi)摩擦角對其的影響。

    黏聚力隨機場中水平與豎直相關長度的比值不同,邊坡穩(wěn)定性隨著相關長度的改變而變化的規(guī)律不同。當,h/,v=1時,水平相關長度趨于無窮大時得到的隨機場近似于均質(zhì)分布的情況,隨著,h的增加,平均安全系數(shù)趨于確定性計算結果,變異系數(shù)、失效概率趨于0,最危險滑面逐漸集中于確定性計算得到的滑面。當,h/,v=10時,,h越長,越易形成連續(xù)的低黏聚力區(qū),發(fā)生邊坡滑動的可能性越高,可能出現(xiàn)的滑面分布范圍越廣。

    4.3 同時考慮兩者的空間變異性的影響

    本節(jié)計算同時考慮了內(nèi)摩擦角與黏聚力的空間變異性,假定和之間的相關系數(shù)分別為?0.5、0以及0.5,計算共分為3組,和均服從=3的Weibull分布,均值分別為20 kPa和15°,水平相關長度取20 m,豎直相關長度取2 m。

    表7所示為不同相關系數(shù)下邊坡安全系數(shù)的統(tǒng)計參數(shù),當和由負相關變?yōu)檎嚓P時,平均安全系數(shù)降低,變異系數(shù)、失效概率增加,與文獻[5,9]中的規(guī)律一致。圖16所示為不同相關系數(shù)下最危險滑面分布范圍,隨著相關系數(shù)的增加,最危險滑面分布越分散。

    (a) 平均安全系數(shù);(b) 失效概率

    表7 不同相關系數(shù)對應的安全系數(shù)分布

    圖16 不同相關系數(shù)對應的最危險滑面分布

    5 結論

    1) 所提的隨機有限元方法利用ABAQUS提供的用戶子程序接口USDFLD編寫程序模擬考慮土體抗剪強度參數(shù)空間變異性的邊坡失穩(wěn)過程,運用python腳本進行后處理,便于從邊坡的安全系數(shù)和最危險滑面分布2方面進行邊坡穩(wěn)定性分析,操作簡單易于實現(xiàn),可為復雜邊坡的穩(wěn)定分析提供一條有效的途徑。

    2) 內(nèi)摩擦角隨機場中均質(zhì)度越高,邊坡越穩(wěn)定,而相關長度的改變對邊坡的穩(wěn)定性影響不明顯;內(nèi)摩擦角的空間變異性對邊坡穩(wěn)定性的影響小于黏聚力的影響。

    3) 黏聚力隨機場中水平與豎直相關長度的比值不同,邊坡的穩(wěn)定性隨著相關長度的改變而變化的規(guī)律不同。當比值為1時,水平相關長度越大,平均安全系數(shù)越高,變異系數(shù)、失效概率越低,最危險滑面分布越集中,邊坡越穩(wěn)定;當比值為10時,水平相關長度越大,低黏聚力區(qū)越連續(xù),越易出現(xiàn)邊坡失穩(wěn)破壞現(xiàn)象。

    4) 內(nèi)摩擦角與黏聚力的相關性直接影響邊坡的穩(wěn)定性,相關系數(shù)越大,平均安全系數(shù)越低,變異系數(shù)、失效概率越高,最危險滑面分布越分散,邊坡的穩(wěn)定性越差。

    參考文獻:

    [1] LI K S, LUMB P. Probabilistic design of slopes[J]. Canadian Geotechnical Journal, 1987, 24(4): 520?535.

    [2] GRIFFITHS D V, FENTON G A. Probabilistic slope stability analysis by finite elements[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2004, 130(5): 507?518.

    [3] GRIFFITHS D V, HUANG J S, FENTON G A. Influence of spatial variability on slope reliability using 2-D random fields[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2009, 135(10): 1367?1378.

    [4] SHEN Hong, FU Helin. Spatial variability and slope reliability analysis[J]. Electronic Journal of Geotechnical Engineering, 2011, 16(L): 1261?1276.

    [5] CHO S E. Probabilistic assessment of slope stability that considers the spatial variability of soil properties[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2010, 136(7): 975?984.

    [6] 譚曉慧, 王建國, 劉新榮, 等. 邊坡穩(wěn)定的有限元可靠度計算及敏感性分析[J]. 巖石力學與工程學報, 2007, 26(1): 115?122. TAN Xiaohui, WANG Jianguo, LIU Xinrong, el al. Finite element reliability computation and sensitivity analysis of slope stability[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(1): 115?122.

    [7] 吳振君, 葛修潤, 王水林. 考慮地質(zhì)成因的土坡可靠度分析[J]. 巖石力學與工程學報, 2011, 30(9): 1904?1911. WU Zhenjun, GE Xiurun, WANG Shuilin. Reliability analysis of soil slope stability considering geologic origin[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(9): 1904?1911.

    [8] 李典慶, 蔣水華, 周創(chuàng)兵. 基于非侵入式隨機有限元法的地下洞室可靠度分析[J]. 巖土工程學報, 2012, 34(1): 123?129. LI Dianqing, JIANG Shuihua, ZHOU Chuangbing. Reliability analysis of underground rock caverns using non-intrusive stochastic finite element method[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(1): 123?129.

    [9] 李典慶, 蔣水華, 周創(chuàng)兵, 等. 考慮參數(shù)空間變異性的邊坡可靠度分析非侵入式隨機有限元法[J]. 巖土工程學報, 2013, 35(8): 1413?1422. LI Dianqing, JIANG Shuihua, ZHOU Chuangbing, el al. Reliability analysis of slopes considering spatial variability of soil parameters using non-intrusive stochastic finite element method[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1413?1422.

    [10] 蔣水華, 馮曉波, 李典慶, 等. 邊坡可靠度分析的非侵入式隨機有限元法[J]. 巖土力學, 2013, 34(8): 2347?2354. JIANG Shuihua, FENG Xiaobo, LI Dianqing, et al. Reliability analysis of slope using non-instrusive stochastic finite element method[J]. Rock and Soil Mechanics, 2013, 34(8): 2347?2354.

    [11] 祁小輝, 李典慶, 周創(chuàng)兵, 等. 考慮土體空間變異性的邊坡最危險滑動面隨機分析方法[J]. 巖土工程學報, 2013, 35(4): 745?753. QI Xiaohui, LI Dianqing, ZHOU Chuangbing, et al. Stochastic analysis method of critical slip surfaces in soil slopes considering spatial variability[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(4): 745?753.

    [12] 肖特, 李典慶, 周創(chuàng)兵, 等. 基于有限元強度折減法的多層邊坡非侵入式可靠度分析[J]. 應用基礎與工程科學學報, 2014, 22(4): 718?732. XIAO Te, LI Dianqing, ZHOU Chuangbing, et al. Non-intrusive reliability analysis of multi-layered slopes using strength reduction FEM[J]. Journal of Basic Science and Engineering, 2014, 22(4): 718?732.

    [13] ZIENKIEWICZ O C, HUMPHESON C, LEWIS R W. Associated and non-associated visco-plasticity and plasticity in soil mechanics [J]. Geotechnique, 1975, 25(4): 671?689.

    [14] GRIFFITH D V, LANE P A. Slope stability analysis by finite elements[J]. Geotechnique, 1999, 49(3): 387?403.

    [15] DAWSON E M, ROTH W H, DRESCHER A. Slope stability analysis by strength reduction[J]. Geotechnique, 1999, 49(6): 835?840.

    [16] 欒茂田, 武亞軍, 年延凱. 強度折減有限元法中邊坡失穩(wěn)的塑形區(qū)判據(jù)及其應用[J]. 防災減災工程學報, 2003, 23(3): 1?7. LUAN Maotian, WU Yajun, NIAN Yankai. A criterion for evaluating slope stability based on development of plastic zone by shear strength reduction FEM[J]. Journal of Disaster Prevention and Mitigation Engineering, 2003, 23(3): 1?7.

    [17] 劉金龍, 欒茂田, 趙少飛, 等. 關于強度折減有限元方法中邊坡失穩(wěn)判據(jù)的討論[J]. 巖土力學, 2005, 26(8): 1345?1348. LIU Jinlong, LUAN Maotian, ZHAO Shaofei, el al. Discussion on criteria for evaluating stability of slope in elastoplastic FEM based on shear strength reduction technique[J]. Rock and Soil Mechanics, 2005, 26(8): 1345?1348.

    [18] 呂慶, 孫紅月, 尚岳全. 強度折減有限元法中邊坡失穩(wěn)判據(jù)的研究[J]. 浙江大學學報(工學版), 2008, 42(1): 83?87. Lü Qing, SUN Hongyue, SHANG Yuequan. Slope failure criteria of shear strength reduction finite element method[J]. Journal of Zhejiang University (Engineering Science), 2008, 42(1): 83?87.

    [19] 孫冠華, 鄭宏, 李春光. 基于等效塑性應變的三維邊坡滑面搜索[J]. 巖土力學, 2010, 31(2): 627?632. SUN Guanhua, ZHENG Hong, LI Chunguang. Searching critical surface of three-dimensional slopes based on equivalent plastic strain[J]. Rock and Soil Mechanics, 2010, 31(2): 627?632.

    [20] EL-RAMLY H, MORGENSTERN N R, CRUDEN D M. Probabilistic stability analysis of a tailings dyke on presheared clay-shale[J]. Canadian Geotechnical Journal, 2003, 40(1): 192?208.

    (編輯 陳愛華)

    Stability analysis of slope considering spatial variation of shear strength parameters

    YUAN Wei1, 2, CHANG Xiaolin1, DUAN Yin1, 3, CHENG Yonggang1, MA Gang1

    (1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China;2. Hubei Provincial Water Resources and Hydropower Planning Survey and Design Institute, Wuhan 430064, China;3. Changjiang Institute of Survey, Planning Design and Research, Wuhan 430010, China)

    A stochastic finite element method was proposed to analyze slope stability based on ABAQUS. Firstly, the Karhunen-Loeve (K-L) expansion method was used to build the random field model of shear strength parameters. Thereafter, the strength reduction finite element method was adopted to simulate the process of slope failure. The criteria for evaluating stability of slope were connectivity of plastic zone and convergence of numerical computation. Finally, the distribution of critical slip surfaces was determined by python script. The proposed method was used to study the effect of spatial variation of internal friction angle and spatial cohesion on the slope stability. The results indicate that the proposed stochastic finite element method is easy to operate and implement. It can analyze the effect of spatial variation of soil parameter on the slope stability from two aspects: safety factor and the critical slip surface. The stability of a slope with different correlation length of friction is similar. The variation trend of slope stability is contrary with two different ratio of horizontal correlation length to vertical correlation length while considering spatial variation of cohesion. Besides, as the correlation coefficient of soil parameters increases, the failure probability is higher, and the distribution of critical slip surface is wider.

    spatial variation; slope stability; strength reduction finite element method; stochastic finite element method

    10.11817/j.issn.1672-7207.2016.11.037

    TU43

    A

    1672?7207(2016)11?3899?10

    2016?01?08;

    2016?03?27

    長江科學院開發(fā)研究基金資助項目(CKWV2015224/KY);國家自然科學基金資助項目(51379161) (Project (CKWV2015224/KY) supported by the CRSRI Open Research Program; Project(51379161) supported by the National Natural Science Foundation of China)

    袁葳,博士研究生,從事高壩結構設計理論與數(shù)值仿真方面研究;E-mail: yw_0201@126.com

    猜你喜歡
    滑面黏聚力摩擦角
    季凍黃土地區(qū)框錨邊坡預應力錨桿內(nèi)力計算方法
    應用摩擦角,巧解動力學問題
    基于滑面正應力修正的錨固邊坡穩(wěn)定性分析方法
    河南科學(2020年3期)2020-06-02 08:30:24
    含軟弱土層邊坡穩(wěn)定性分析的全局滑面自動搜索技術
    公路工程(2020年2期)2020-05-15 04:10:38
    基于矢量和法的三維邊坡滑面搜索研究
    借助摩擦角 快解勻速運動問題
    黏聚力強度對滑面作用的差異分析
    科技視界(2016年7期)2016-04-01 20:15:33
    摩擦角在平衡問題中的應用
    用摩擦角巧解靜力學問題
    基于準黏聚力原理的加筋黏土強度理論淺析
    四川建筑(2014年6期)2014-09-03 11:00:08
    欧美 日韩 精品 国产| 狠狠狠狠99中文字幕| 久久久水蜜桃国产精品网| 人妻一区二区av| 91精品三级在线观看| www.999成人在线观看| 男女床上黄色一级片免费看| 国产蜜桃级精品一区二区三区 | 满18在线观看网站| 69av精品久久久久久| 国产精品久久视频播放| 欧美精品人与动牲交sv欧美| 丝袜人妻中文字幕| 久久久久国内视频| 亚洲av欧美aⅴ国产| 日韩人妻精品一区2区三区| 高清黄色对白视频在线免费看| 少妇裸体淫交视频免费看高清 | 免费一级毛片在线播放高清视频 | 日韩视频一区二区在线观看| 十八禁网站免费在线| 欧美黑人欧美精品刺激| 国产成人精品在线电影| 久久久久久久午夜电影 | 黑人巨大精品欧美一区二区mp4| 真人做人爱边吃奶动态| 伦理电影免费视频| 欧美丝袜亚洲另类 | 两人在一起打扑克的视频| 国产精品美女特级片免费视频播放器 | 男人的好看免费观看在线视频 | 看片在线看免费视频| 狂野欧美激情性xxxx| 精品熟女少妇八av免费久了| 最近最新免费中文字幕在线| 亚洲欧美一区二区三区黑人| 制服诱惑二区| 亚洲一区高清亚洲精品| 精品午夜福利视频在线观看一区| 91老司机精品| 亚洲一区二区三区不卡视频| 亚洲专区国产一区二区| 两个人看的免费小视频| 欧美中文综合在线视频| 国产人伦9x9x在线观看| 国产成人av激情在线播放| 国产精品久久电影中文字幕 | 国产熟女午夜一区二区三区| 免费观看人在逋| 国产精品久久久人人做人人爽| 国产成人欧美| 精品卡一卡二卡四卡免费| √禁漫天堂资源中文www| 国产精华一区二区三区| 俄罗斯特黄特色一大片| 露出奶头的视频| 日韩三级视频一区二区三区| 久久亚洲精品不卡| 国产真人三级小视频在线观看| 亚洲国产精品sss在线观看 | 国产亚洲精品一区二区www | av有码第一页| av天堂在线播放| 国产精品98久久久久久宅男小说| 脱女人内裤的视频| 大片电影免费在线观看免费| 熟女少妇亚洲综合色aaa.| 亚洲免费av在线视频| 久久精品国产99精品国产亚洲性色 | 啦啦啦 在线观看视频| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 国产高清激情床上av| 天天影视国产精品| 国产亚洲一区二区精品| 丁香六月欧美| 丰满迷人的少妇在线观看| 亚洲欧洲精品一区二区精品久久久| 免费久久久久久久精品成人欧美视频| 免费观看精品视频网站| 国产精品电影一区二区三区 | 国产亚洲一区二区精品| 脱女人内裤的视频| 黄片播放在线免费| 欧美大码av| 国产精品一区二区免费欧美| av电影中文网址| 美女视频免费永久观看网站| 久久天堂一区二区三区四区| 免费女性裸体啪啪无遮挡网站| 日本a在线网址| 交换朋友夫妻互换小说| a在线观看视频网站| 精品熟女少妇八av免费久了| 一级片'在线观看视频| 成人三级做爰电影| 国产极品粉嫩免费观看在线| 亚洲人成电影免费在线| av片东京热男人的天堂| 中文字幕人妻熟女乱码| 嫁个100分男人电影在线观看| 在线观看日韩欧美| 狂野欧美激情性xxxx| svipshipincom国产片| 欧美久久黑人一区二区| 搡老乐熟女国产| 成年人黄色毛片网站| 亚洲av美国av| 久久精品人人爽人人爽视色| 最近最新中文字幕大全免费视频| 亚洲精华国产精华精| 国产不卡一卡二| 欧美日韩黄片免| 我的亚洲天堂| 国产精华一区二区三区| 天天操日日干夜夜撸| 又大又爽又粗| 久久人妻福利社区极品人妻图片| 国产91精品成人一区二区三区| 日本五十路高清| 国产成人欧美| 成人特级黄色片久久久久久久| 欧美日韩av久久| 国产免费男女视频| 亚洲av成人一区二区三| 日日夜夜操网爽| 人人澡人人妻人| 国产高清videossex| 精品午夜福利视频在线观看一区| 免费在线观看亚洲国产| 免费观看a级毛片全部| 国产黄色免费在线视频| а√天堂www在线а√下载 | 久久人妻av系列| 欧美日韩瑟瑟在线播放| 国产在线一区二区三区精| 国产一区有黄有色的免费视频| 国产真人三级小视频在线观看| 女人高潮潮喷娇喘18禁视频| www.自偷自拍.com| 新久久久久国产一级毛片| 精品一区二区三区视频在线观看免费 | 一进一出好大好爽视频| 免费在线观看完整版高清| 久久国产精品男人的天堂亚洲| 日韩中文字幕欧美一区二区| 国产精品久久久久久精品古装| 天天躁夜夜躁狠狠躁躁| 精品一品国产午夜福利视频| 亚洲男人天堂网一区| 国产亚洲精品久久久久5区| 欧美精品人与动牲交sv欧美| 欧美在线一区亚洲| 桃红色精品国产亚洲av| 精品欧美一区二区三区在线| 精品高清国产在线一区| 黄网站色视频无遮挡免费观看| 搡老岳熟女国产| 另类亚洲欧美激情| 久久精品aⅴ一区二区三区四区| 午夜影院日韩av| av线在线观看网站| av欧美777| 日本黄色日本黄色录像| 亚洲av熟女| 最新在线观看一区二区三区| 制服人妻中文乱码| 精品久久蜜臀av无| 女人精品久久久久毛片| 成人精品一区二区免费| 亚洲欧美一区二区三区久久| 在线观看免费视频日本深夜| 午夜福利视频在线观看免费| 国产欧美亚洲国产| 亚洲精品中文字幕一二三四区| 国产伦人伦偷精品视频| 国产成人影院久久av| 日韩有码中文字幕| 亚洲免费av在线视频| 国产亚洲欧美精品永久| 脱女人内裤的视频| 国产精品98久久久久久宅男小说| 韩国精品一区二区三区| 高清视频免费观看一区二区| 久久精品人人爽人人爽视色| 99精品欧美一区二区三区四区| 国产又色又爽无遮挡免费看| 亚洲第一青青草原| 男女下面插进去视频免费观看| 国产熟女午夜一区二区三区| 欧美一级毛片孕妇| tocl精华| 亚洲色图 男人天堂 中文字幕| 午夜老司机福利片| 大码成人一级视频| 亚洲成国产人片在线观看| 免费av中文字幕在线| www.精华液| 精品少妇久久久久久888优播| 国产在视频线精品| 一本大道久久a久久精品| 久久国产亚洲av麻豆专区| 在线观看免费高清a一片| 91在线观看av| 美国免费a级毛片| 99国产精品一区二区三区| 欧美精品亚洲一区二区| 久久国产精品影院| 国产99白浆流出| 好男人电影高清在线观看| 女人高潮潮喷娇喘18禁视频| 丁香欧美五月| 80岁老熟妇乱子伦牲交| 80岁老熟妇乱子伦牲交| 国产蜜桃级精品一区二区三区 | 色94色欧美一区二区| 美女高潮到喷水免费观看| 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| av一本久久久久| 免费日韩欧美在线观看| 人妻丰满熟妇av一区二区三区 | 亚洲国产欧美网| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区三| 久久久久久久精品吃奶| 婷婷丁香在线五月| 午夜福利影视在线免费观看| 国产亚洲精品一区二区www | 久久中文字幕一级| 在线观看日韩欧美| 亚洲色图综合在线观看| 亚洲全国av大片| 91国产中文字幕| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av高清一级| 久久人妻福利社区极品人妻图片| 中亚洲国语对白在线视频| 久久午夜亚洲精品久久| 久久午夜综合久久蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 久久久精品区二区三区| 亚洲成人手机| 国产在视频线精品| 黑人巨大精品欧美一区二区mp4| 午夜福利视频在线观看免费| 十八禁网站免费在线| 青草久久国产| 丰满迷人的少妇在线观看| 国产精品一区二区精品视频观看| 久久国产精品男人的天堂亚洲| 中文字幕另类日韩欧美亚洲嫩草| 欧美亚洲 丝袜 人妻 在线| 日本一区二区免费在线视频| 精品第一国产精品| 妹子高潮喷水视频| 久热这里只有精品99| 欧美黄色片欧美黄色片| 91av网站免费观看| 黄色片一级片一级黄色片| 两个人看的免费小视频| 一级a爱视频在线免费观看| 国产亚洲精品第一综合不卡| 精品视频人人做人人爽| 亚洲成人国产一区在线观看| 大陆偷拍与自拍| 黄色片一级片一级黄色片| 脱女人内裤的视频| 中文字幕制服av| 国产一卡二卡三卡精品| 精品人妻在线不人妻| 91精品三级在线观看| 免费人成视频x8x8入口观看| 午夜老司机福利片| 久久久久精品国产欧美久久久| 香蕉久久夜色| 悠悠久久av| 亚洲片人在线观看| 人人妻人人爽人人添夜夜欢视频| 中国美女看黄片| 大片电影免费在线观看免费| 免费在线观看亚洲国产| 久久精品国产99精品国产亚洲性色 | 亚洲七黄色美女视频| 男人操女人黄网站| 国产激情久久老熟女| 亚洲中文日韩欧美视频| 国产亚洲av高清不卡| 天天躁狠狠躁夜夜躁狠狠躁| 捣出白浆h1v1| 亚洲五月天丁香| 色婷婷av一区二区三区视频| 人人澡人人妻人| 下体分泌物呈黄色| 涩涩av久久男人的天堂| 欧美黄色片欧美黄色片| 大型av网站在线播放| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲国产一区二区在线观看 | 热99re8久久精品国产| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲av高清不卡| 免费少妇av软件| 99re6热这里在线精品视频| 久久人妻熟女aⅴ| 久久精品亚洲av国产电影网| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品影院| 国产精品免费一区二区三区在线 | 18禁裸乳无遮挡动漫免费视频| 黄色丝袜av网址大全| 欧美精品av麻豆av| 国产成+人综合+亚洲专区| 国产在线观看jvid| 人成视频在线观看免费观看| 亚洲欧洲精品一区二区精品久久久| 少妇 在线观看| 亚洲美女黄片视频| 两个人免费观看高清视频| 亚洲精品美女久久av网站| 成年人免费黄色播放视频| av超薄肉色丝袜交足视频| 搡老乐熟女国产| 日本黄色视频三级网站网址 | 首页视频小说图片口味搜索| 成年女人毛片免费观看观看9 | 无限看片的www在线观看| 一进一出抽搐gif免费好疼 | 欧美黑人精品巨大| 亚洲七黄色美女视频| 国产激情久久老熟女| 中文字幕精品免费在线观看视频| 国产主播在线观看一区二区| 青草久久国产| 成人三级做爰电影| 热re99久久国产66热| 新久久久久国产一级毛片| 中国美女看黄片| 伦理电影免费视频| 天天操日日干夜夜撸| 高清毛片免费观看视频网站 | 欧美日韩中文字幕国产精品一区二区三区 | 久久人妻熟女aⅴ| 国产精品一区二区在线不卡| 999精品在线视频| 久久亚洲精品不卡| 国产精品av久久久久免费| 男女之事视频高清在线观看| 成熟少妇高潮喷水视频| 午夜两性在线视频| 成人亚洲精品一区在线观看| 中出人妻视频一区二区| 久热爱精品视频在线9| 一a级毛片在线观看| 欧美老熟妇乱子伦牲交| 男女之事视频高清在线观看| 亚洲精品乱久久久久久| 国产免费男女视频| 波多野结衣av一区二区av| 精品免费久久久久久久清纯 | 飞空精品影院首页| 黄色视频不卡| 一边摸一边抽搐一进一小说 | av欧美777| 两性夫妻黄色片| 日日夜夜操网爽| 老司机深夜福利视频在线观看| 午夜视频精品福利| 欧美黑人欧美精品刺激| 日韩大码丰满熟妇| 丰满人妻熟妇乱又伦精品不卡| 亚洲色图av天堂| 人妻 亚洲 视频| 国产精品久久久久久人妻精品电影| 国产高清激情床上av| 日韩三级视频一区二区三区| 久久香蕉激情| 久久国产亚洲av麻豆专区| 好男人电影高清在线观看| videosex国产| 精品无人区乱码1区二区| 看黄色毛片网站| 十八禁高潮呻吟视频| 叶爱在线成人免费视频播放| 国产97色在线日韩免费| 在线观看舔阴道视频| 久久国产乱子伦精品免费另类| 一夜夜www| 国产男女内射视频| 国产亚洲av高清不卡| 天天影视国产精品| 精品国产一区二区三区久久久樱花| 日韩有码中文字幕| 欧美国产精品va在线观看不卡| 日韩精品免费视频一区二区三区| 人妻一区二区av| 1024香蕉在线观看| 精品免费久久久久久久清纯 | 手机成人av网站| 国产一区二区三区综合在线观看| 18禁裸乳无遮挡动漫免费视频| 国产亚洲av高清不卡| 国产av一区二区精品久久| 欧美在线一区亚洲| 大码成人一级视频| 国产精品久久久久久精品古装| 极品教师在线免费播放| 欧美国产精品一级二级三级| 成人国产一区最新在线观看| 丝袜美腿诱惑在线| 高清毛片免费观看视频网站 | 成年人免费黄色播放视频| 国产精品美女特级片免费视频播放器 | 亚洲精品国产一区二区精华液| 男女高潮啪啪啪动态图| 曰老女人黄片| 亚洲第一青青草原| 一进一出抽搐动态| 岛国在线观看网站| 在线永久观看黄色视频| 欧美亚洲 丝袜 人妻 在线| 美国免费a级毛片| 日韩欧美三级三区| 久久久久国产一级毛片高清牌| 少妇粗大呻吟视频| 欧美精品人与动牲交sv欧美| 在线国产一区二区在线| 久久精品aⅴ一区二区三区四区| 亚洲中文av在线| 一级作爱视频免费观看| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 丰满人妻熟妇乱又伦精品不卡| 色婷婷久久久亚洲欧美| 精品亚洲成国产av| 久久精品国产亚洲av高清一级| 男女午夜视频在线观看| 女性生殖器流出的白浆| 国产人伦9x9x在线观看| 国产男靠女视频免费网站| 欧美日韩成人在线一区二区| 国产成人精品久久二区二区免费| 丁香欧美五月| 欧美精品一区二区免费开放| 咕卡用的链子| 午夜福利一区二区在线看| 免费观看a级毛片全部| 欧美乱色亚洲激情| 建设人人有责人人尽责人人享有的| 久久久水蜜桃国产精品网| 老司机午夜福利在线观看视频| 又紧又爽又黄一区二区| 视频在线观看一区二区三区| xxxhd国产人妻xxx| 久久精品熟女亚洲av麻豆精品| 午夜老司机福利片| 亚洲熟女毛片儿| 俄罗斯特黄特色一大片| 成人18禁在线播放| ponron亚洲| 亚洲国产精品一区二区三区在线| 制服人妻中文乱码| 亚洲av成人不卡在线观看播放网| 夜夜爽天天搞| 丝袜美足系列| 黄色女人牲交| 午夜老司机福利片| 国产99白浆流出| 电影成人av| 好看av亚洲va欧美ⅴa在| 99热只有精品国产| 婷婷丁香在线五月| 91av网站免费观看| 首页视频小说图片口味搜索| 日本黄色视频三级网站网址 | 久久国产精品人妻蜜桃| 在线观看一区二区三区激情| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲av高清不卡| 日韩欧美免费精品| xxxhd国产人妻xxx| av在线播放免费不卡| 极品人妻少妇av视频| 亚洲精品国产区一区二| 国产精品久久视频播放| 中文字幕精品免费在线观看视频| 国产男女超爽视频在线观看| 脱女人内裤的视频| 黄片小视频在线播放| 亚洲第一欧美日韩一区二区三区| 飞空精品影院首页| 国产有黄有色有爽视频| 欧美久久黑人一区二区| 国产野战对白在线观看| 嫁个100分男人电影在线观看| 亚洲av成人一区二区三| 在线视频色国产色| 国产国语露脸激情在线看| 久久久久久免费高清国产稀缺| 夜夜爽天天搞| 亚洲一区中文字幕在线| 俄罗斯特黄特色一大片| 黄色怎么调成土黄色| 91成人精品电影| 免费在线观看完整版高清| 国产成人欧美在线观看 | 啦啦啦 在线观看视频| 91大片在线观看| 日韩欧美一区视频在线观看| 99久久精品国产亚洲精品| 亚洲中文字幕日韩| 亚洲七黄色美女视频| 美女福利国产在线| 亚洲第一青青草原| 大片电影免费在线观看免费| 国产精品影院久久| 国产精品久久视频播放| 精品久久久久久久久久免费视频 | 操美女的视频在线观看| 日本黄色视频三级网站网址 | 美国免费a级毛片| av有码第一页| 我的亚洲天堂| 少妇被粗大的猛进出69影院| 18禁裸乳无遮挡动漫免费视频| 国产成人欧美| 一级片免费观看大全| 一本综合久久免费| 黄色视频不卡| www.熟女人妻精品国产| 亚洲国产精品sss在线观看 | av欧美777| 夜夜夜夜夜久久久久| 少妇粗大呻吟视频| 最新的欧美精品一区二区| www日本在线高清视频| 成人亚洲精品一区在线观看| xxx96com| 久热这里只有精品99| 久久久久视频综合| 婷婷精品国产亚洲av在线 | 一二三四在线观看免费中文在| 天天操日日干夜夜撸| 国产精品99久久99久久久不卡| 亚洲精品自拍成人| 黑人巨大精品欧美一区二区mp4| 18在线观看网站| 欧美成人午夜精品| 久久精品亚洲av国产电影网| 亚洲国产欧美日韩在线播放| 在线观看免费午夜福利视频| 欧美 日韩 精品 国产| 两性夫妻黄色片| 麻豆成人av在线观看| 国产成人精品无人区| 精品久久久精品久久久| 国产免费现黄频在线看| 在线观看免费高清a一片| 国产精品.久久久| 成人精品一区二区免费| 俄罗斯特黄特色一大片| 99国产精品一区二区三区| 高清黄色对白视频在线免费看| 岛国在线观看网站| 国产精品久久视频播放| 女同久久另类99精品国产91| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品中文字幕在线视频| 欧美日韩福利视频一区二区| 日韩免费av在线播放| 亚洲五月婷婷丁香| 不卡一级毛片| 欧美日韩福利视频一区二区| 日韩制服丝袜自拍偷拍| 国产成人av教育| 婷婷精品国产亚洲av在线 | 伊人久久大香线蕉亚洲五| 三级毛片av免费| 99国产极品粉嫩在线观看| 久久久精品国产亚洲av高清涩受| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 自线自在国产av| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼 | 少妇粗大呻吟视频| 亚洲人成电影免费在线| 中文字幕制服av| 免费在线观看亚洲国产| 国产日韩欧美亚洲二区| 十分钟在线观看高清视频www| 黑人巨大精品欧美一区二区蜜桃| 久久久精品国产亚洲av高清涩受| 国产在线一区二区三区精| 亚洲成a人片在线一区二区| videos熟女内射| 中出人妻视频一区二区| 成年动漫av网址| 久久午夜综合久久蜜桃| 国产欧美日韩精品亚洲av| 免费av中文字幕在线| 一本一本久久a久久精品综合妖精| 精品国产乱子伦一区二区三区| 女人久久www免费人成看片| 亚洲午夜理论影院| 中文字幕人妻丝袜制服| 国产精品国产高清国产av | 99久久99久久久精品蜜桃| 俄罗斯特黄特色一大片| 免费在线观看视频国产中文字幕亚洲| 男女之事视频高清在线观看| 岛国在线观看网站|