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

    降雨入滲條件下軟巖邊坡穩(wěn)定性分析

    2012-01-08 07:12:44付宏淵王桂堯蔣中明楊晴霞
    巖土力學(xué) 2012年8期
    關(guān)鍵詞:非飽和軟巖暫態(tài)

    付宏淵,曾 鈴,王桂堯,蔣中明,楊晴霞

    (1. 長(zhǎng)沙理工大學(xué) 土木與建筑學(xué)院,長(zhǎng)沙 410076;2. 長(zhǎng)沙理工大學(xué) 水利工程學(xué)院,長(zhǎng)沙 410076)

    1 引 言

    統(tǒng)計(jì)資料表明[1],邊坡失穩(wěn)的主要誘因在于降雨入滲。在日本,每年由于降雨引起的邊坡失穩(wěn)多達(dá)上萬(wàn)個(gè)。而在我國(guó)[2],天然巖坡和鐵路、公路等開(kāi)挖形成的人工邊坡發(fā)生失穩(wěn)破壞也多在雨季。近些年來(lái),廣大學(xué)者普遍認(rèn)識(shí)到了邊坡失穩(wěn)的發(fā)生與降雨存在著關(guān)聯(lián)性。并認(rèn)為降雨導(dǎo)致邊坡原有非飽和區(qū)負(fù)孔隙水壓力發(fā)生比較強(qiáng)烈的變化。按照非飽和抗剪強(qiáng)度理論分析,認(rèn)為負(fù)孔隙水壓力的降低甚至消失將在很大程度上降低邊坡土體抗剪強(qiáng)度,并且入滲雨水將在一定程度上暫時(shí)地增大降雨入滲區(qū)域邊坡土體自重,這對(duì)邊坡穩(wěn)定性是極其不利的。文獻(xiàn)[3-6]利用有限元軟件Geostudio(seep)對(duì)降雨條件下的邊坡滲流場(chǎng)進(jìn)行了分析,并以此為基礎(chǔ)計(jì)算了邊坡安全系數(shù)在降雨過(guò)程中的變化,文中主要考慮了邊坡非飽和區(qū)負(fù)孔隙水壓力的變化對(duì)邊坡穩(wěn)定性的影響。文獻(xiàn)[7-8]對(duì)降雨的入滲過(guò)程及暫態(tài)飽和區(qū)的形成及消散過(guò)程進(jìn)行了深入的研究,提出了暫態(tài)飽和區(qū)中由于降雨入滲所增加的邊坡巖土體自重將影響邊坡穩(wěn)定性的觀點(diǎn)。由此可見(jiàn),目前研究降雨對(duì)邊坡穩(wěn)定性的影響,主要從降雨入滲將增加邊坡表層暫態(tài)飽和區(qū)的自重及雨水使入滲深度內(nèi)的負(fù)孔隙水壓力減小這兩個(gè)方面來(lái)考慮。但對(duì)于軟巖邊坡而言,考慮這兩點(diǎn)還不能對(duì)邊坡穩(wěn)定性進(jìn)行全面的分析。軟巖邊坡在降雨條件下由于雨水對(duì)巖石的軟化作用,多表現(xiàn)為邊坡表面局部垮塌與整體滑動(dòng)相結(jié)合的破壞模式,所以在分析降雨條件下邊坡穩(wěn)定性時(shí)考慮巖石遇水產(chǎn)生的軟化效應(yīng)是必不可少的。另一方面,在邊坡穩(wěn)定性計(jì)算中,大家普遍采用的商用極限平衡計(jì)算軟件在考慮降雨過(guò)程中的滲透作用時(shí),計(jì)算結(jié)果存在一定誤差[9]。一些學(xué)者在認(rèn)識(shí)這一缺陷的基礎(chǔ)上,采用自行編制計(jì)算程序的方法加以解決,但由于其計(jì)算程序編制的復(fù)雜性以及運(yùn)用的針對(duì)性較強(qiáng),所以很難得到廣泛的推廣。

    本文對(duì)軟巖邊坡在降雨條件下的穩(wěn)定性分析既考慮了軟巖遇水軟化導(dǎo)致的巖石物理力學(xué)參數(shù)的變化,同時(shí)還能如實(shí)地反映孔隙水壓力對(duì)邊坡穩(wěn)定性的影響。本文研究的目的是通過(guò)數(shù)值計(jì)算來(lái)研究降雨條件下邊坡中孔隙水壓力的存在及暫態(tài)飽和區(qū)內(nèi)巖體物理力學(xué)參數(shù)的變化共同作用下的邊坡破壞模式及穩(wěn)定性在空間和時(shí)間上的演化規(guī)律,為軟巖邊坡失穩(wěn)的預(yù)測(cè)和采取合理的治理措施提供參考。

    2 降雨入滲條件下軟巖邊坡穩(wěn)定性分析方法

    根據(jù)軟巖邊坡在降雨條件下存在的實(shí)際情況,采用如下方法來(lái)對(duì)降雨條件下的穩(wěn)定性進(jìn)行分析。

    (1)采用飽和-非飽和滲流計(jì)算理論對(duì)設(shè)計(jì)降雨條件下邊坡滲流場(chǎng)進(jìn)行研究,得出不同降雨時(shí)刻邊坡體內(nèi)的孔隙水壓力、暫態(tài)飽和區(qū)在時(shí)間和空間上的分布。

    (2)以二維飽和-非飽和滲流計(jì)算結(jié)果為基礎(chǔ),將不同時(shí)刻孔隙水壓力wu 導(dǎo)入FLAC3D中,利用強(qiáng)度折減法計(jì)算相應(yīng)時(shí)刻邊坡穩(wěn)定性。對(duì)于Mohr-Coulomb 材料,強(qiáng)度折減安全系數(shù)可以表示為

    其中:

    式中:σ 、 'σ 分別為破壞面上的法向應(yīng)力及法向有效應(yīng)力;τ 、 'τ 分別為任一平面的抗剪強(qiáng)度及有效剪應(yīng)力;c φ、 為土體實(shí)際的抗剪強(qiáng)度指標(biāo);c φ′ ′、為土體達(dá)到極限狀態(tài)時(shí)的抗剪強(qiáng)度指標(biāo);sF 為折減系數(shù)。

    這一方法很好地避免了普通商用極限平衡計(jì)算軟件計(jì)算邊坡穩(wěn)定性時(shí)不能將飽和-非飽和滲流計(jì)算結(jié)果中節(jié)點(diǎn)孔隙水壓力進(jìn)行精確表達(dá)這一缺陷。

    (3)利用室內(nèi)巖石飽水軟化試驗(yàn)所得的軟巖物理力學(xué)參數(shù)變化規(guī)律來(lái)對(duì)降雨條件下邊坡暫態(tài)飽和區(qū)域內(nèi)巖體的物理力學(xué)參數(shù)進(jìn)行取值,再結(jié)合巖體內(nèi)孔隙水壓力的變化(步驟(2))來(lái)研究降雨條件下暫態(tài)飽和區(qū)內(nèi)巖石軟化對(duì)邊坡穩(wěn)定性的影響。

    考慮降雨條件下的邊坡穩(wěn)定性計(jì)算是在二維飽和-非飽和計(jì)算軟件及FLAC3D軟件平臺(tái)上實(shí)現(xiàn)的,通過(guò)在FLAC3D中編寫(xiě)用戶子程序,將某一時(shí)刻飽和-非飽和滲流計(jì)算所得到的孔隙水壓力導(dǎo)入FLAC3D中,同時(shí),將現(xiàn)場(chǎng)取樣的巖石進(jìn)行軟化試驗(yàn)獲得的不同飽水時(shí)刻的巖石力學(xué)強(qiáng)度參數(shù)賦值于對(duì)應(yīng)時(shí)刻的暫態(tài)飽和區(qū)中。最終運(yùn)用強(qiáng)度折減法對(duì)該時(shí)刻的邊坡穩(wěn)定性進(jìn)行計(jì)算。計(jì)算流程如圖1 所示。

    圖1 邊坡穩(wěn)定性計(jì)算流程 Fig.1 Slope stability calculation process

    3 降雨入滲飽和-非飽和計(jì)算

    3.1 計(jì)算模型及分析方案

    湖南省郴州境內(nèi)某高速公路K59+060 處路塹邊坡,邊坡表層有1 m 厚含碎石粉質(zhì)黏土,其下為 粉砂質(zhì)泥巖,土黃色原巖結(jié)構(gòu)清晰可辨,長(zhǎng)度為 5.9 m 巖芯呈土柱狀?yuàn)A塊狀。地下水類型主要有:松散堆積層孔隙水、基巖裂隙水2 大類型。由于該地區(qū)降雨量較大,孔隙水直接受大氣降水補(bǔ)給,大氣降水進(jìn)入表層后,大部分儲(chǔ)積在表層的孔隙中,其含水性、透水性較好,其補(bǔ)給、排泄受地形條件及本身土體結(jié)構(gòu)的控制,總體上表現(xiàn)為順坡向坡腳多向滲流,部分以散流形式向地勢(shì)較低處匯集排泄,相當(dāng)一部分通過(guò)裂隙下滲補(bǔ)給基巖裂隙水,進(jìn)行較深層的水循環(huán)。邊坡剖面及初始水位線如圖2 所示。網(wǎng)格模型見(jiàn)圖3。計(jì)算網(wǎng)格單元數(shù)量為7 036,節(jié)點(diǎn)數(shù)量為7 237。

    圖2 邊坡典型剖面及初始地下水位 Fig.2 Profile of typical slope and initial water level

    圖3 二維有限元計(jì)算模型 Fig.3 2D finite element model

    3.2 飽和-非飽和滲流計(jì)算原理

    在降雨條件下軟巖邊坡滲流場(chǎng)分析中,由于雨水入滲或者地下水位埋藏較淺,飽和區(qū)與非飽和區(qū)的地下水的流動(dòng)是相互聯(lián)系的,需將兩者統(tǒng)一起來(lái)研究,可將此表述為飽和-非飽和滲流問(wèn)題,所以以土-水勢(shì)為因變量來(lái)建立統(tǒng)一的控制方程。飽和-非飽和滲流與飽和滲流一樣,滿足達(dá)西定律土體內(nèi)非恒定滲流,其偏微分方程形式如下[10]:

    式中:h 為巖體裂隙中的總水頭;xk 和yk 為x 和y方向的滲透系數(shù);w 為源匯項(xiàng);wm 為比水容量;wρ為水的密度;g 為重力加速度;t 為時(shí)間。

    巖土體非恒定滲流有限元方程為

    式中:[ ]K 為單元特征矩陣;[ ]M 為單元質(zhì)量矩陣;{ }Q 為節(jié)點(diǎn)流量向量。

    有限元方程求解結(jié)合以下邊界條件進(jìn)行:

    式中:1S 為已知水頭邊界;2S 為已知流量邊界。

    非恒定滲流分析的初始條件為

    在飽和-非飽和滲流數(shù)值分析中,巖土體的滲透系數(shù)與單元的飽和度、孔隙水壓力相關(guān),在缺乏現(xiàn)場(chǎng)試驗(yàn)資料的情況下較為常用的確定滲透系數(shù)的方法為利用經(jīng)典公式Van Genuchten[11]模型進(jìn)行擬合,模型表達(dá)式為

    式中:θ 為體積含水率;rθ 為殘余含水率;sθ 為飽和含水率;H 為負(fù)壓;sK 為飽和滲透系數(shù);α 、m、n 為土-水特征曲線形狀參數(shù)。

    3.3 計(jì)算參數(shù)、初始條件

    全風(fēng)化-中風(fēng)化粉砂質(zhì)泥巖飽和滲透系數(shù)選取現(xiàn)場(chǎng)壓水試驗(yàn)5 個(gè)試驗(yàn)孔的中平均值:1.0×10-8cm/s。粉質(zhì)黏土飽和滲透系數(shù)根據(jù)現(xiàn)場(chǎng)取樣進(jìn)行室內(nèi)滲透試驗(yàn)確定:取1.0×10-7cm/s。飽和含水率則根據(jù)室內(nèi)試驗(yàn)測(cè)得[12],粉質(zhì)黏土取0.2,全風(fēng)化-中風(fēng)化粉砂質(zhì)泥巖取0.1。由于坡體內(nèi)非飽和區(qū)滲透系數(shù)的大小由基質(zhì)吸力及初始含水率共同確定,本文中非飽和區(qū)滲透系數(shù)的大小采用前文所述Van Genuchten模型進(jìn)行描述,滲透系數(shù)及含水率隨基質(zhì)吸力變化如圖4、5 所示,由圖可知,土體體積含水率和滲透系數(shù)都與基質(zhì)吸力存在非線性關(guān)系,基質(zhì)吸力絕對(duì)值越大,土體體積含水率與滲透系數(shù)越小。

    雨水的入滲對(duì)初始滲流場(chǎng)較為敏感,因此,需要通過(guò)求解一個(gè)穩(wěn)態(tài)的滲流場(chǎng)來(lái)獲得瞬態(tài)滲流計(jì)算的初始條件。穩(wěn)態(tài)滲流計(jì)算結(jié)果則作為飽和-非飽和計(jì)算的初始滲流場(chǎng)。飽和-非飽和計(jì)算過(guò)程中,邊坡表面邊界設(shè)置為單位流量邊界q(根據(jù)降雨強(qiáng)度計(jì)算所得),其他邊界設(shè)置為不透水邊界(Q =0)。

    限于篇幅原因,初始狀態(tài)下的滲流場(chǎng)示意圖未列出。

    圖4 體積含水率與負(fù)孔隙水壓力關(guān)系曲線圖(θ -p 曲線) Fig.4 Curve of volatile water content and negative pore water pressure

    圖5 滲透系數(shù)與負(fù)孔隙水壓力關(guān)系曲線圖(k -p 曲線) Fig.5 Curve of conductivity and negative pore water pressure

    3.4 計(jì)算方案

    根據(jù)郴州氣象局提供的氣象資料顯示,郴州地區(qū)近50 a 來(lái)最大21 d 降雨量為1 556.67 mm,換算為單位降雨強(qiáng)度為q=3.23×10-7m/s,降雨歷時(shí)21 d,計(jì)算總時(shí)長(zhǎng)為30 d,單位時(shí)間步長(zhǎng)為1 h。

    3.5 飽和-非飽和滲流計(jì)算結(jié)果分析

    按照前文所述計(jì)算方法,運(yùn)用二維有限元軟件對(duì)邊坡在降雨條件下的滲流場(chǎng)進(jìn)行分析,由于本文主要考慮孔隙水壓力的變化及暫態(tài)飽和區(qū)內(nèi)巖石軟化對(duì)邊坡穩(wěn)定性的影響,故只對(duì)上述兩個(gè)方面進(jìn)行分析。

    (1)降雨入滲對(duì)孔隙水壓力的影響

    邊坡體內(nèi)孔隙水壓力的大小及分布在降雨后發(fā)生了較大的變化,圖6 為邊坡在整個(gè)降雨過(guò)程中及降雨結(jié)束后一段時(shí)間內(nèi)的3 個(gè)時(shí)刻的孔隙水壓力分布情況。由圖可知,邊坡在降雨作用下邊坡體內(nèi)滲流場(chǎng)將發(fā)生明顯的變化,變化主要發(fā)生在邊坡表層。降雨21 d 時(shí)入滲深度達(dá)到5 m。降雨導(dǎo)致邊坡表層非飽和區(qū)負(fù)孔隙水壓力變?yōu)檎?。降雨停止后由于雨水的繼續(xù)入滲,正孔隙水壓力區(qū)逐漸變小,最后完全消失,但過(guò)程緩慢。

    圖6 孔隙水壓力等值線 (單位:kPa) Fig.6 Contours of pore water pressure (unit: kPa)

    (2)降雨入滲作用下暫態(tài)飽和區(qū)的變化規(guī)律

    由圖7 可知,在降雨持續(xù)一段時(shí)間后,邊坡表面出現(xiàn)暫態(tài)飽和區(qū),其面積隨著降雨歷時(shí)的增長(zhǎng)而逐漸擴(kuò)大,在降雨停止一段時(shí)間后其面積開(kāi)始減小。通過(guò)研究暫態(tài)飽和區(qū)的面積大小與形成時(shí)間的關(guān)系,對(duì)于巖石飽和軟化室內(nèi)試驗(yàn)確定巖石飽水時(shí)間具有指導(dǎo)意義。

    圖7 暫態(tài)飽和區(qū)及地下水位的變化 Fig.7 Change processes of transient saturated zone area and the movement of groundwater

    4 軟巖邊坡穩(wěn)定性分析

    4.1 模型的建立與參數(shù)的取值

    數(shù)值分析選擇K59+060 處典型剖面,以圖3 所示二維網(wǎng)格為基礎(chǔ)進(jìn)行準(zhǔn)三維等效模型的建立。為了考慮暫態(tài)飽和區(qū)內(nèi)巖石遇水軟化的時(shí)效問(wèn)題,邊坡表層采用分層建模的形式,以此來(lái)準(zhǔn)確地對(duì)暫態(tài)飽和區(qū)遇水軟化后強(qiáng)度大小進(jìn)行賦值。整個(gè)計(jì)算模型網(wǎng)格單元數(shù)為7 036,節(jié)點(diǎn)數(shù)為14 474。圖8 為數(shù)值計(jì)算網(wǎng)格圖。

    圖8 三維數(shù)值計(jì)算模型網(wǎng)格圖 Fig.8 Mesh of 3D FEM

    邊界條件:鉛直邊界上施加水平位移約束;水平邊界上施加鉛直位移約束。

    初始條件:為正確模擬巖體參數(shù)變化對(duì)邊坡穩(wěn)定性的影響,初始地應(yīng)力場(chǎng)按自重應(yīng)力場(chǎng)考慮,即首先采用彈性模型進(jìn)行分析,然后將計(jì)算模型中的位移置0,保留單元應(yīng)力,再導(dǎo)入滲流計(jì)算網(wǎng)格節(jié)點(diǎn)信息(孔隙水壓力)。圖9 為導(dǎo)入孔隙水壓力后的孔隙水壓力分布圖(降雨21 d)。

    本構(gòu)模型:采用摩爾-庫(kù)侖屈服準(zhǔn)則的彈塑性本構(gòu)模型。

    計(jì)算力學(xué)參數(shù):根據(jù)邊坡的地質(zhì)勘探報(bào)告選取不考慮巖石遇水軟化的邊坡物理力學(xué)參數(shù),取值見(jiàn)表1。根據(jù)試驗(yàn)規(guī)范[13]將現(xiàn)場(chǎng)采取的粉砂質(zhì)泥巖巖樣制備成若干Φ50 mm×100 mm 圓柱體試件進(jìn)行不同飽水時(shí)間下的單軸壓縮試驗(yàn),得到單軸抗壓強(qiáng)度,并在此基礎(chǔ)上利用基于 GSI 法的廣義 Hoek- Brown[14]準(zhǔn)則來(lái)確定本試驗(yàn)試件的抗剪強(qiáng)度及相關(guān)計(jì)算參數(shù)隨飽水時(shí)間的變化。粉砂質(zhì)泥巖飽水軟化試驗(yàn)所得主要物理力學(xué)性質(zhì)參數(shù)及變化規(guī)律如表1和圖10 所示,由圖可知,隨著飽水時(shí)間的持續(xù),其主要物理力學(xué)參數(shù)呈現(xiàn)下降趨勢(shì),飽水試驗(yàn)前期下降幅度大于試驗(yàn)后期。

    計(jì)算方案:選取降雨持續(xù)3、6、12、18、21 d、停止4、7、9 d 為穩(wěn)定性計(jì)算時(shí)間點(diǎn),對(duì)邊坡在降雨過(guò)程中的穩(wěn)定性進(jìn)行模擬,與此同時(shí),考慮邊坡表層暫態(tài)飽和區(qū)軟巖力學(xué)參數(shù)軟化效應(yīng)對(duì)邊坡穩(wěn)定性的影響??偨Y(jié)出邊坡在降雨條件下安全系數(shù)及塑性區(qū)變化規(guī)律。

    圖10 粉砂質(zhì)泥巖飽水試驗(yàn)主要力學(xué)參數(shù)變化圖 Fig.10 Experimental results of saturated silty mudstones

    4.2 邊坡塑性區(qū)分布及發(fā)展規(guī)律分析

    圖11 為考慮降雨及巖石軟化共同作用下的邊坡塑性區(qū)分布圖。(限于篇幅,文中只列出了最具代表性的降雨6 d、21 d、降雨停止9 d 這3 個(gè)時(shí)間點(diǎn)的邊坡塑性區(qū)分布圖)由圖可知,降雨初期(T=6 d)塑性區(qū)主要集中在邊坡表層,并隨著降雨時(shí)間的持續(xù)沿著坡腳向上延伸。這是由于在降雨作用下,入滲雨水由于重力作用迅速向坡腳匯集,在邊坡表層形成暫態(tài)飽和區(qū)。暫態(tài)飽和區(qū)內(nèi)孔隙水負(fù)壓力逐漸變?yōu)檎紫端畨毫?,且自重增大,從而使邊坡表層巖體有效應(yīng)力降低。與此同時(shí),邊坡表層巖體在雨水的軟化作用下也進(jìn)一步加速了表層巖體的塑性區(qū)的形成。這在實(shí)際工程中表現(xiàn)為軟巖邊坡表層局部分層坍塌。當(dāng)持續(xù)降雨到21 d 的過(guò)程中,邊坡表層塑性區(qū)進(jìn)一步沿坡面向上及向內(nèi)延伸,塑性區(qū)面積及連通性進(jìn)一步加大。在降雨入滲深度的進(jìn)一步擴(kuò)大所引起的入滲深度內(nèi),在基質(zhì)吸力降低甚至消失和巖土體重度增大的共同作用下,邊坡內(nèi)部出現(xiàn)了深層的、具有一定貫通性的塑性區(qū)分布,且塑性區(qū)在坡腳區(qū)域分布面積最大,這對(duì)邊坡整體穩(wěn)定性是極為不利的。降雨停止9 d 時(shí),由于邊坡表層暫態(tài)飽和區(qū)內(nèi)雨水向坡內(nèi)繼續(xù)入滲及排出坡外,使邊坡表層塑性區(qū)面積有一定的擴(kuò)大,但并不明顯。由于雨水繼續(xù)入滲的影響,邊坡深層塑性區(qū)分布面積變化擴(kuò)大,形成了具有貫通性塑性區(qū)分布帶,這為邊坡潛在滑動(dòng)面的產(chǎn)生提供了客觀條件。

    圖11 邊坡塑性區(qū)分布圖 Fig.11 Distribution of plastic zone

    4.3 邊坡瞬態(tài)安全系數(shù)變化分析

    在對(duì)邊坡在降雨條件下的塑性區(qū)分布及發(fā)展規(guī)律分析的基礎(chǔ)上,結(jié)合強(qiáng)度折減法計(jì)算邊坡瞬態(tài)安全系數(shù)。安全系數(shù)變化規(guī)律如表2 所示。由表可知,隨著降雨時(shí)間的持續(xù),邊坡整體安全系數(shù)有顯著的降低,由1.65 降為1.15。且由計(jì)算可知,在降雨結(jié)束后仍然有一定的降低,但降低速率明顯小于降雨過(guò)程中安全系數(shù)的降低速率。這一現(xiàn)象說(shuō)明:降雨持續(xù)過(guò)程中,暫態(tài)飽和區(qū)影響面積內(nèi)暫態(tài)孔隙水壓力的升高、邊坡巖體自重的增大以及雨水對(duì)巖石的軟化對(duì)軟巖邊坡穩(wěn)定性是不利的;降雨結(jié)束后,由于暫態(tài)飽和區(qū)雨水的繼續(xù)入滲,使邊坡內(nèi)負(fù)壓區(qū)進(jìn)一步受到影響,入滲區(qū)域內(nèi)非飽和抗剪強(qiáng)度降低。同時(shí),暫態(tài)飽和區(qū)部分雨水由坡面排出,這在一定程度上又減緩了巖石軟化作用以及邊坡自重對(duì)整體穩(wěn)定性的影響,因此,在降雨結(jié)束后邊坡安全系數(shù)仍有小幅降低。

    表2 不同時(shí)刻軟巖邊坡安全系數(shù) Table 2 Different moments soft rock slope stability safety coefficient

    5 結(jié) 論

    (1)本文提出的降雨條件下軟巖邊坡穩(wěn)定性分析方法既能考慮到滲流作用下坡體內(nèi)孔隙水壓力的變化對(duì)邊坡穩(wěn)定性的影響,又能將軟巖遇水所存在的軟化效應(yīng)在邊坡穩(wěn)定性計(jì)算中較好地體現(xiàn)。為分析降雨條件下的軟巖邊坡穩(wěn)定性提供了一種新的方法。

    (2)降雨開(kāi)始后,邊坡表層非飽和區(qū)出現(xiàn)暫態(tài)正孔隙水壓力。隨著降雨歷時(shí)的增長(zhǎng),邊坡表層開(kāi)始出現(xiàn)暫態(tài)飽和區(qū),暫態(tài)飽和區(qū)面積開(kāi)始從坡腳處沿著坡面向上延伸,同時(shí)也朝著與坡面垂直的方向向邊坡內(nèi)部擴(kuò)大。

    (3)邊坡在降雨及巖石軟化效應(yīng)的共同作用 下,在降雨初期,邊坡表層率先出現(xiàn)塑性區(qū),表現(xiàn)為邊坡表面局部分層坍塌;隨著降雨歷時(shí)的增長(zhǎng),邊坡表層塑性區(qū)面積增大,同時(shí)在邊坡內(nèi)部也出現(xiàn)具有貫通性的深層塑性區(qū);降雨停止后,邊坡內(nèi)各位置塑性區(qū)面積還有所增加,但增幅較小。

    (4)邊坡在降雨過(guò)程中,安全系數(shù)有顯著降低,這是邊坡體內(nèi)暫態(tài)孔隙水壓力與巖石軟化效應(yīng)共同影響的結(jié)果。降雨停止后,雨水由坡面緩慢排出及暫態(tài)孔隙水壓力的逐漸降低對(duì)軟巖邊坡安全系數(shù)的繼續(xù)降低具有一定的延緩作用。

    [1] NG CW, SHI Q. A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage[J]. Computers and Geotechanics, 1998, 22(1): 1-28.

    [2] 孫廣忠. 中國(guó)典型滑坡[M]. 北京: 科學(xué)出版社, 1998.

    [3] 劉俊新, 劉育田, 胡啟軍. 非飽和地表徑流-滲流和流固體耦合條件下降雨入滲對(duì)路堤邊坡穩(wěn)定性研究[J]. 巖土力學(xué), 2010, 31(3): 903-910. LIU Jun-xin, LIU Yu-tian, HU Qi-jun. Stability of embankment slope subjected to rainfall infiltration considering both runoff-underground seepage and fluid-solid coupling[J]. Rock and Soil Mechanics, 2010, 31(3): 903-910.

    [4] 周家文, 徐衛(wèi)亞, 鄧俊曄, 等. 降雨入滲條件下邊坡的穩(wěn)定性分析[J]. 水利學(xué)報(bào), 2008, 39(9): 1066-1073. ZHOU Jia-wen, XU Wei-ya, DENG Jun-ye, et al. Analysis of slope stability under intensive rainfall[J]. Journal of Hydraulic Engineering, 2008, 39(9): 1066-1073.

    [5] 吳長(zhǎng)富, 朱向榮, 尹小濤, 等. 強(qiáng)降雨條件下土質(zhì)邊坡瞬態(tài)穩(wěn)定性分析[J]. 巖土力學(xué), 2008, 29(2): 386-391. WU Chang-fu, ZHU Xiang-rong, YIN Xiao-tao, et al. Analysis of soil slope’ transient stability under intensive rainfall[J]. Rock and Soil Mechanics, 2008, 29(2): 386-391.

    [6] 張國(guó), 朱濟(jì)祥, 莊軍, 等. 巖質(zhì)邊坡穩(wěn)定性的滲流耦合分析[J]. 天津大學(xué)學(xué)報(bào), 1999,32(5): 560-564. ZHANG Guo, ZHU Ji-xiang, ZHUANG Jun, et al. Rocky slope hydromechanical coupling analysis[J]. Journal of Tianjin University, 1999,32(5): 560-564.

    [7] 戚國(guó)慶, 黃潤(rùn)秋, 速寶玉, 等. 巖質(zhì)邊坡降雨入滲過(guò)程的數(shù)值模擬[J]. 巖石力學(xué)與工程學(xué)報(bào), 2003, 22(4): 625-629. QI Guo-qing, HUANG Run-qiu, SU Bao-yu, et al. Numeric simulation on rainfall infiltration on rock slope[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(4): 625-629.

    [8] 周桂云, 李同春. 飽和-非飽和非穩(wěn)定滲流作用下巖質(zhì)邊坡穩(wěn)定性分析[J]. 水電能源科學(xué), 2006, 24(5): 79-82. ZHOU Gui-yun, LI Tong-chun. The rock slope stability analysis under saturated-unsaturated seepage[J]. Water Resources and Power, 2006, 24(5): 79-82.

    [9] 毛昶熙, 段祥寶, 李祖貽. 滲流數(shù)值計(jì)算與程序應(yīng)用[M]. 南京: 河海大學(xué)出版社, 1999.

    [10] 顧慰慈. 滲流計(jì)算原理及應(yīng)用[M]. 北京: 中國(guó)建材工業(yè)出版社, 2000.

    [11] VAN GENUCHTEN M T. A closed form equation for prediction the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1979, 44(5): 892-898.

    [12] 錢家歡, 殷宗澤. 土工原理與計(jì)算[M]. 北京: 中國(guó)水利水電出版社, 1996.

    [13] 中華人民共和國(guó)水利部. SL264-2001 水利水電工程巖石試驗(yàn)規(guī)范[S]. 內(nèi)蒙古: 內(nèi)蒙古大學(xué)出版社, 2003.

    [14] HOEK E, BROWN E T. Empirical strength criterion for rock masses[J]. Journal of the Geotechnical Engineering Division, 1980, 106(9): 1013-1035.

    猜你喜歡
    非飽和軟巖暫態(tài)
    300Mvar空冷隱極同步調(diào)相機(jī)暫態(tài)特性仿真分析
    非飽和原狀黃土結(jié)構(gòu)強(qiáng)度的試驗(yàn)研究
    電力系統(tǒng)全網(wǎng)一體化暫態(tài)仿真接口技術(shù)
    電子制作(2018年14期)2018-08-21 01:38:28
    軟巖作為面板堆石壩填筑料的探討
    官帽舟水電站軟巖筑壩技術(shù)的成功探索
    非飽和多孔介質(zhì)應(yīng)力滲流耦合分析研究
    除氧器暫態(tài)計(jì)算研究
    非飽和土基坑剛性擋墻抗傾覆設(shè)計(jì)與參數(shù)分析
    非飽和地基土蠕變特性試驗(yàn)研究
    基于PSD-BPA的暫態(tài)穩(wěn)定控制批處理計(jì)算方法的實(shí)現(xiàn)
    中文在线观看免费www的网站| 又黄又爽又免费观看的视频| 怎么达到女性高潮| 变态另类成人亚洲欧美熟女| 久久这里只有精品中国| 亚洲av第一区精品v没综合| 成人性生交大片免费视频hd| 亚洲在线自拍视频| 欧美区成人在线视频| 此物有八面人人有两片| 婷婷六月久久综合丁香| 亚洲精品影视一区二区三区av| 最新在线观看一区二区三区| 性插视频无遮挡在线免费观看| 久久久久久久午夜电影| 69人妻影院| 国产精品久久久久久人妻精品电影| 国内精品久久久久久久电影| 在线播放无遮挡| 一区二区三区激情视频| 国产伦精品一区二区三区视频9| 国产伦人伦偷精品视频| 极品教师在线视频| 白带黄色成豆腐渣| www.色视频.com| 亚洲黑人精品在线| 成人亚洲精品av一区二区| 久久久久久久久久成人| 国产精品久久久久久亚洲av鲁大| 我要看日韩黄色一级片| 久久草成人影院| 国产熟女xx| 国产精品影院久久| 男女之事视频高清在线观看| 最新在线观看一区二区三区| 国产精品国产高清国产av| 无遮挡黄片免费观看| 嫩草影院精品99| 18禁黄网站禁片午夜丰满| 人妻丰满熟妇av一区二区三区| 全区人妻精品视频| 变态另类丝袜制服| 免费在线观看亚洲国产| 亚洲欧美日韩东京热| 久久国产乱子伦精品免费另类| 中文字幕熟女人妻在线| 久久精品国产99精品国产亚洲性色| 日韩欧美在线乱码| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 一级黄片播放器| 色综合婷婷激情| 久久久久亚洲av毛片大全| 欧美日韩中文字幕国产精品一区二区三区| 亚洲一区高清亚洲精品| 好看av亚洲va欧美ⅴa在| av在线观看视频网站免费| a级毛片a级免费在线| 国产精品美女特级片免费视频播放器| 91在线观看av| eeuss影院久久| 亚洲av第一区精品v没综合| 99久国产av精品| 一个人看视频在线观看www免费| 日韩欧美精品v在线| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 亚洲av二区三区四区| 最新中文字幕久久久久| 国产成人影院久久av| 黄色丝袜av网址大全| 亚洲经典国产精华液单 | 亚洲欧美日韩卡通动漫| www.熟女人妻精品国产| 成年女人看的毛片在线观看| 亚洲熟妇熟女久久| 99热这里只有是精品在线观看 | 一个人观看的视频www高清免费观看| 校园春色视频在线观看| 亚洲精品在线美女| 国产精品爽爽va在线观看网站| 午夜激情福利司机影院| 欧美激情久久久久久爽电影| 久久精品久久久久久噜噜老黄 | 久久久久精品国产欧美久久久| 亚洲精品一卡2卡三卡4卡5卡| 黄色女人牲交| 午夜精品久久久久久毛片777| 99精品在免费线老司机午夜| av中文乱码字幕在线| 欧美性猛交╳xxx乱大交人| 色视频www国产| 欧美激情国产日韩精品一区| 精品一区二区免费观看| 欧美色视频一区免费| 美女被艹到高潮喷水动态| 日本免费一区二区三区高清不卡| 天美传媒精品一区二区| 精品日产1卡2卡| 看十八女毛片水多多多| 国产国拍精品亚洲av在线观看| 国产美女午夜福利| 91狼人影院| 亚洲精品日韩av片在线观看| a在线观看视频网站| 十八禁国产超污无遮挡网站| 99国产极品粉嫩在线观看| 美女xxoo啪啪120秒动态图 | 亚洲无线在线观看| 国产精品久久久久久久久免 | 我要看日韩黄色一级片| 欧美高清性xxxxhd video| 国产乱人伦免费视频| 国产欧美日韩一区二区三| 日本一二三区视频观看| 99热精品在线国产| 一级a爱片免费观看的视频| 色播亚洲综合网| 国产在视频线在精品| 成年女人看的毛片在线观看| 国产精品女同一区二区软件 | 国语自产精品视频在线第100页| 久久久成人免费电影| 亚洲最大成人手机在线| 特级一级黄色大片| 美女黄网站色视频| 国产视频内射| 午夜两性在线视频| av在线天堂中文字幕| 婷婷丁香在线五月| 动漫黄色视频在线观看| 免费大片18禁| 尤物成人国产欧美一区二区三区| 成年女人看的毛片在线观看| av在线天堂中文字幕| 欧美另类亚洲清纯唯美| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 国产成人欧美在线观看| 亚洲 欧美 日韩 在线 免费| 有码 亚洲区| 欧美3d第一页| 国产精品98久久久久久宅男小说| 亚洲精品粉嫩美女一区| 国产成人欧美在线观看| www日本黄色视频网| 日本a在线网址| 又爽又黄a免费视频| 亚洲成人久久性| 久久久久久国产a免费观看| 国产视频内射| 一级黄色大片毛片| 国产不卡一卡二| 国产高清视频在线观看网站| 精品一区二区三区视频在线观看免费| 日韩精品青青久久久久久| 亚洲精品粉嫩美女一区| 亚洲人与动物交配视频| 99精品在免费线老司机午夜| 亚洲精品亚洲一区二区| 久久久久久久午夜电影| 久9热在线精品视频| 内地一区二区视频在线| 少妇熟女aⅴ在线视频| 在线观看一区二区三区| 男女视频在线观看网站免费| 久久人妻av系列| 亚洲中文字幕一区二区三区有码在线看| 99在线人妻在线中文字幕| 精品人妻熟女av久视频| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 丰满的人妻完整版| 久久九九热精品免费| 免费人成在线观看视频色| 成人毛片a级毛片在线播放| 我要看日韩黄色一级片| 亚洲av不卡在线观看| 国语自产精品视频在线第100页| 久久婷婷人人爽人人干人人爱| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品久久男人天堂| 国产精品亚洲美女久久久| 深夜a级毛片| av在线天堂中文字幕| 国产精品伦人一区二区| 18+在线观看网站| 婷婷精品国产亚洲av在线| 噜噜噜噜噜久久久久久91| 禁无遮挡网站| 直男gayav资源| 国产伦在线观看视频一区| 色哟哟哟哟哟哟| 久久这里只有精品中国| 亚洲av熟女| 成年女人毛片免费观看观看9| 亚洲av中文字字幕乱码综合| 国产老妇女一区| 少妇丰满av| 日本a在线网址| 日本黄色片子视频| 波多野结衣巨乳人妻| 精品一区二区三区人妻视频| 亚洲,欧美精品.| 一边摸一边抽搐一进一小说| 精品久久久久久成人av| 久久久久国产精品人妻aⅴ院| 婷婷色综合大香蕉| 天美传媒精品一区二区| 特大巨黑吊av在线直播| 天天一区二区日本电影三级| 12—13女人毛片做爰片一| 99久久99久久久精品蜜桃| 757午夜福利合集在线观看| 18禁黄网站禁片免费观看直播| 免费大片18禁| 成年女人毛片免费观看观看9| 欧美绝顶高潮抽搐喷水| 亚洲精品日韩av片在线观看| 一区二区三区四区激情视频 | 欧美色视频一区免费| 久久精品国产亚洲av天美| 国产亚洲精品综合一区在线观看| 久久热精品热| 日本精品一区二区三区蜜桃| 99riav亚洲国产免费| 在线播放国产精品三级| 欧美最新免费一区二区三区 | 亚洲五月天丁香| 国产精品亚洲av一区麻豆| 内射极品少妇av片p| 午夜免费成人在线视频| 国产精品三级大全| 99在线视频只有这里精品首页| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 亚洲av日韩精品久久久久久密| 国产成人aa在线观看| netflix在线观看网站| 一进一出抽搐gif免费好疼| 亚洲精品色激情综合| 国模一区二区三区四区视频| 丰满人妻熟妇乱又伦精品不卡| 三级男女做爰猛烈吃奶摸视频| 深夜a级毛片| 51午夜福利影视在线观看| bbb黄色大片| netflix在线观看网站| 精品人妻1区二区| 舔av片在线| 成年免费大片在线观看| netflix在线观看网站| 久久精品人妻少妇| 久久久久久国产a免费观看| 美女 人体艺术 gogo| 免费av观看视频| 欧美+亚洲+日韩+国产| a级毛片a级免费在线| 免费在线观看成人毛片| 搞女人的毛片| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 欧美高清性xxxxhd video| 在线观看av片永久免费下载| 欧美在线一区亚洲| 久久人妻av系列| 日韩免费av在线播放| 亚洲狠狠婷婷综合久久图片| 色视频www国产| 欧美xxxx黑人xx丫x性爽| 午夜精品在线福利| 亚洲av中文字字幕乱码综合| 一级a爱片免费观看的视频| 精品一区二区三区av网在线观看| 99久久99久久久精品蜜桃| 亚洲va日本ⅴa欧美va伊人久久| av天堂中文字幕网| 日韩欧美一区二区三区在线观看| 男插女下体视频免费在线播放| 亚洲av不卡在线观看| 午夜福利在线观看免费完整高清在 | 亚洲av二区三区四区| 国产野战对白在线观看| 99热只有精品国产| 国模一区二区三区四区视频| 99精品在免费线老司机午夜| 全区人妻精品视频| 久久人人爽人人爽人人片va | 亚洲av免费高清在线观看| 九色国产91popny在线| av福利片在线观看| 日韩欧美三级三区| 久久久国产成人免费| 午夜免费男女啪啪视频观看 | 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡| 嫩草影院精品99| 如何舔出高潮| 变态另类成人亚洲欧美熟女| 亚洲18禁久久av| 99在线视频只有这里精品首页| 9191精品国产免费久久| 一a级毛片在线观看| 俺也久久电影网| 美女被艹到高潮喷水动态| 国产真实伦视频高清在线观看 | 无人区码免费观看不卡| 日韩中文字幕欧美一区二区| a在线观看视频网站| 国语自产精品视频在线第100页| 一个人免费在线观看电影| 成人精品一区二区免费| АⅤ资源中文在线天堂| 啦啦啦观看免费观看视频高清| 黄色女人牲交| 在线播放国产精品三级| 香蕉av资源在线| 亚洲精品乱码久久久v下载方式| 亚洲成人精品中文字幕电影| 欧美丝袜亚洲另类 | 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区| 国产精品三级大全| 午夜视频国产福利| 亚洲无线在线观看| 国产亚洲精品综合一区在线观看| 舔av片在线| 中文字幕高清在线视频| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添小说| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 麻豆国产97在线/欧美| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 精品人妻一区二区三区麻豆 | 全区人妻精品视频| 给我免费播放毛片高清在线观看| 国产午夜福利久久久久久| 十八禁人妻一区二区| 在线a可以看的网站| 一个人观看的视频www高清免费观看| 人妻制服诱惑在线中文字幕| 欧美日韩乱码在线| 亚洲最大成人中文| 男女视频在线观看网站免费| 麻豆成人av在线观看| 天天躁日日操中文字幕| 在线a可以看的网站| 中亚洲国语对白在线视频| 成年版毛片免费区| 在线观看午夜福利视频| 中文亚洲av片在线观看爽| 成人无遮挡网站| 国产黄片美女视频| 99久久无色码亚洲精品果冻| 欧美最新免费一区二区三区 | 伊人久久精品亚洲午夜| 欧美zozozo另类| 日韩国内少妇激情av| 性色av乱码一区二区三区2| 亚洲无线观看免费| 嫩草影院新地址| 欧美+日韩+精品| 18禁在线播放成人免费| 亚洲欧美日韩东京热| 日本a在线网址| 99国产精品一区二区三区| 亚洲性夜色夜夜综合| 在线观看舔阴道视频| av在线天堂中文字幕| 国产午夜精品论理片| 免费观看人在逋| 国产精品亚洲美女久久久| 免费观看人在逋| 99国产精品一区二区三区| 嫩草影院新地址| 欧美精品啪啪一区二区三区| 日日干狠狠操夜夜爽| 欧美日韩瑟瑟在线播放| 日本 欧美在线| 变态另类成人亚洲欧美熟女| 毛片女人毛片| 国产高清视频在线播放一区| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 在线看三级毛片| 亚洲乱码一区二区免费版| 欧美黄色淫秽网站| 国内揄拍国产精品人妻在线| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 国内精品久久久久精免费| www.999成人在线观看| 亚洲三级黄色毛片| 丰满乱子伦码专区| 欧美日韩国产亚洲二区| 91在线精品国自产拍蜜月| 激情在线观看视频在线高清| 午夜福利在线观看吧| 又紧又爽又黄一区二区| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 一个人看视频在线观看www免费| а√天堂www在线а√下载| 黄色一级大片看看| 俺也久久电影网| 亚洲中文字幕一区二区三区有码在线看| 麻豆国产av国片精品| 国产乱人视频| 欧美中文日本在线观看视频| 国产高清激情床上av| 中文在线观看免费www的网站| 亚洲18禁久久av| 国产伦精品一区二区三区四那| 色综合站精品国产| 国产伦一二天堂av在线观看| 欧美黑人巨大hd| 久久香蕉精品热| 一本一本综合久久| 人人妻,人人澡人人爽秒播| 女人被狂操c到高潮| 久久精品国产自在天天线| 国产精品永久免费网站| 久久人人爽人人爽人人片va | 亚洲国产欧洲综合997久久,| 人人妻人人看人人澡| 国内毛片毛片毛片毛片毛片| 一个人看的www免费观看视频| 18美女黄网站色大片免费观看| 精品国产亚洲在线| 狠狠狠狠99中文字幕| 一本综合久久免费| 丝袜美腿在线中文| 久久精品91蜜桃| 欧美日韩福利视频一区二区| 亚洲乱码一区二区免费版| 欧美日韩亚洲国产一区二区在线观看| 黄片小视频在线播放| av专区在线播放| 亚洲av电影在线进入| 国产精品亚洲美女久久久| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 午夜a级毛片| 日韩 亚洲 欧美在线| 欧美区成人在线视频| 欧美黑人巨大hd| 午夜视频国产福利| 亚洲av.av天堂| 色av中文字幕| 国产男靠女视频免费网站| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 日日干狠狠操夜夜爽| 91狼人影院| 嫩草影院精品99| 丁香六月欧美| 亚洲av免费在线观看| 精品久久久久久久久久久久久| 午夜a级毛片| 国产精品久久久久久久久免 | 搡老妇女老女人老熟妇| 9191精品国产免费久久| 国产一区二区在线av高清观看| 国产成人a区在线观看| 婷婷六月久久综合丁香| 国产精品久久久久久久久免 | 久久亚洲真实| 757午夜福利合集在线观看| 日日干狠狠操夜夜爽| 欧美乱色亚洲激情| 舔av片在线| 韩国av一区二区三区四区| 日本免费a在线| 91狼人影院| 嫩草影院新地址| 欧美日韩瑟瑟在线播放| 搡老妇女老女人老熟妇| 欧美乱色亚洲激情| 美女xxoo啪啪120秒动态图 | 琪琪午夜伦伦电影理论片6080| 久久久久久久亚洲中文字幕 | 在线播放国产精品三级| 久久亚洲精品不卡| 亚洲第一欧美日韩一区二区三区| 村上凉子中文字幕在线| 老司机午夜十八禁免费视频| 91午夜精品亚洲一区二区三区 | 一级av片app| 亚洲成人精品中文字幕电影| 99久久精品国产亚洲精品| 又粗又爽又猛毛片免费看| 亚洲人与动物交配视频| 中文字幕av成人在线电影| 18禁黄网站禁片午夜丰满| 好男人电影高清在线观看| 日韩 亚洲 欧美在线| 国产精品久久久久久亚洲av鲁大| 美女 人体艺术 gogo| 欧美+亚洲+日韩+国产| 99在线视频只有这里精品首页| 国内精品一区二区在线观看| 免费看光身美女| 欧美日韩乱码在线| 午夜视频国产福利| 免费高清视频大片| 在线播放国产精品三级| 国产精品av视频在线免费观看| 国产成人啪精品午夜网站| 一区福利在线观看| 国产精品久久电影中文字幕| 亚洲片人在线观看| 亚洲美女搞黄在线观看 | 久久国产精品人妻蜜桃| 国产成年人精品一区二区| 成人美女网站在线观看视频| 精品午夜福利在线看| 中文字幕久久专区| 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| 欧美性猛交黑人性爽| 国产精品人妻久久久久久| 欧美性感艳星| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 国产精品98久久久久久宅男小说| 久久九九热精品免费| 深夜a级毛片| 成人美女网站在线观看视频| а√天堂www在线а√下载| 国产精品一区二区性色av| 精品午夜福利视频在线观看一区| 99在线人妻在线中文字幕| 少妇高潮的动态图| 能在线免费观看的黄片| 热99在线观看视频| 看片在线看免费视频| 毛片一级片免费看久久久久 | 一本久久中文字幕| 国产精品亚洲美女久久久| 久久久久久大精品| 免费av观看视频| 亚洲性夜色夜夜综合| 久久精品国产亚洲av涩爱 | 特大巨黑吊av在线直播| 热99在线观看视频| 观看免费一级毛片| 亚洲熟妇熟女久久| 在线国产一区二区在线| 欧美性感艳星| 黄色一级大片看看| www.www免费av| 色哟哟哟哟哟哟| 中文字幕精品亚洲无线码一区| 色在线成人网| 99热这里只有精品一区| 亚洲激情在线av| 天堂网av新在线| 90打野战视频偷拍视频| 内射极品少妇av片p| 欧美成人a在线观看| 一区二区三区激情视频| 可以在线观看的亚洲视频| 深夜a级毛片| 99久久成人亚洲精品观看| 精品午夜福利在线看| 国产精品亚洲美女久久久| 欧美色视频一区免费| 亚洲片人在线观看| netflix在线观看网站| 日本一二三区视频观看| av国产免费在线观看| 免费人成在线观看视频色| 最近最新免费中文字幕在线| 欧美日本亚洲视频在线播放| 国产在线精品亚洲第一网站| 99国产精品一区二区蜜桃av| 97碰自拍视频| 国产高清三级在线| 99国产精品一区二区蜜桃av| 国产精品久久久久久亚洲av鲁大| 男女那种视频在线观看| 亚洲欧美日韩卡通动漫| 久久欧美精品欧美久久欧美| 亚洲人成网站高清观看| 69av精品久久久久久| 神马国产精品三级电影在线观看| 日本一本二区三区精品| 女人十人毛片免费观看3o分钟| 亚洲专区中文字幕在线| 一本精品99久久精品77| 亚洲七黄色美女视频| 亚洲精品亚洲一区二区| 国产精品一区二区免费欧美| 久久久久久久精品吃奶| 日本五十路高清| 免费在线观看日本一区| 中亚洲国语对白在线视频| 精品久久久久久,| 国产欧美日韩精品一区二区| 性色avwww在线观看| 亚洲精品成人久久久久久| 赤兔流量卡办理| 尤物成人国产欧美一区二区三区| 免费高清视频大片| 免费av不卡在线播放| 国产成人欧美在线观看| 琪琪午夜伦伦电影理论片6080| 日韩中文字幕欧美一区二区| 欧美另类亚洲清纯唯美| 俄罗斯特黄特色一大片| 12—13女人毛片做爰片一|