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

    降雨循環(huán)條件下高切坡穩(wěn)定性演化過(guò)程及預(yù)測(cè)方法

    2021-07-22 10:52:18靳紅華王林峰任青陽(yáng)張星星
    土木建筑與環(huán)境工程 2021年4期
    關(guān)鍵詞:預(yù)測(cè)方法演化過(guò)程數(shù)值模擬

    靳紅華 王林峰 任青陽(yáng) 張星星

    摘 要:為進(jìn)一步揭示降雨循環(huán)條件下基巖型臺(tái)階狀高切坡的降雨入滲及穩(wěn)定性演化過(guò)程,以樂(lè)西高速馬邊至昭覺(jué)段某粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡為研究對(duì)象,通過(guò)室內(nèi)干濕循環(huán)試驗(yàn)建立土體抗剪強(qiáng)度參數(shù)劣化數(shù)學(xué)模型;利用Geo-studio數(shù)值模擬軟件研究多工況降雨循環(huán)下高切坡降雨入滲過(guò)程,揭示不同降雨循環(huán)工況下及雨后高切坡內(nèi)部滲流場(chǎng)及穩(wěn)定性變化規(guī)律,建立降雨型高切坡穩(wěn)定系數(shù)逐年劣化方程;結(jié)合室內(nèi)試驗(yàn)及數(shù)值模擬結(jié)果,建立該類(lèi)高切坡穩(wěn)定系數(shù)預(yù)測(cè)方法。研究結(jié)果表明:降雨入滲過(guò)程中高切坡潛在滑移面形狀及位置并不發(fā)生明顯變化,且表現(xiàn)為高切坡深處圓弧面及基巖積水面的組合型滑面;相同降雨時(shí)間內(nèi),高切坡降雨入滲深度及穩(wěn)定性劣化幅度與降雨強(qiáng)度成正相關(guān);單次降雨循環(huán)周期內(nèi),高切坡穩(wěn)定性劣化幅度與降雨循環(huán)次數(shù)成負(fù)相關(guān);高切坡降雨入滲深度或入滲總量越大,雨后高切坡穩(wěn)定系數(shù)回升越小;高切坡穩(wěn)定系數(shù)劣化系數(shù)采用土體抗剪強(qiáng)度參數(shù)劣化系數(shù)平均值較為安全合理,所建降雨型高切坡穩(wěn)定系數(shù)預(yù)測(cè)方法具有較高的精度。

    關(guān)鍵詞:降雨型高切坡;干濕循環(huán)試驗(yàn);數(shù)值模擬;穩(wěn)定性;演化過(guò)程;預(yù)測(cè)方法

    中圖分類(lèi)號(hào):TU457? ?文獻(xiàn)標(biāo)志碼:A ??文章編號(hào):2096-6717(2020)04-0012-12

    收稿日期:2020-07-06

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(51678097、41472262);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0802206-3);交通運(yùn)輸行業(yè)高層次技術(shù)人才培養(yǎng)項(xiàng)目(2018-26);重慶市自然科學(xué)基金(cstc2020jcyj-msxmX0218);重慶市教委科學(xué)技術(shù)研究項(xiàng)目(KJQN201800706)

    作者簡(jiǎn)介:靳紅華(1995- ),男,博士生,主要從事道路工程重大地質(zhì)災(zāi)害防治技術(shù)研究,E-mail:1589864695@qq.com。

    王林峰(通信作者),男,教授,博士,E-mail:wanglinfeng@cqjtu.edu.cn。

    Received:2020-07-06

    Foundation items:National Natural Science Foundation of China (No. 51678097, 41472262); National Key R & D Program of China (No. 2016YFC0802206-3); High Level Technical Personnel Training Project in Transportation Industry (No. 2018-26); Chongqing Natural Science Foundation (cstc2020jcyj-msxmX0218); Scientific and Technological Research Program of Chongqing Municipal Education Commission (No. KJQN201800706)

    Author brief:JIN Honghua (1995- ), PhD candidate, main research interest: prevention and treatment technology of major geological disasters in road engineering, E-mail: 1589864695@qq.com.

    WANG Linfeng (corresponding author), professor, PhD, E-mail: wanglinfeng@cqjtu.edu.cn.

    Analysis of the stability evolution process and prediction method of high cutting slope under rainfall cycle condition

    JIN Honghua1a, WANG Linfeng1b, REN Qingyang1a, ZHANG Xingxing2

    (1a. State Key Laboratory of Mountain Bridge and Tunnel Engineering; 1b. Mountain Road Water Transport Chongqing Municipal Key Laboratory of Geological Disaster Reduction, Chongqing Jiaotong University, Chongqing 400074, P. R. China; 2. School of Civil Engineering, Chongqing University of Arts and Science, Chongqing 402160, P. R. China)

    Abstract: In order to further reveal the process of rainfall infiltration and stability evolution of bedrock-type stepped high cutting slopes under rainfall cycling conditions, taking the bedrock-type stepped high cutting slope of a silty clay cover layer from Mabian to Zhaojue section of Lexi Expressway as the research object, a mathematical model for the deterioration of soil shear strength parameters was established through indoor dry-wet cycle tests; Geo-studio numerical simulation software was used to study the rainfall infiltration process of the high cutting slope under multi-case rainfall cycles, which reveals the internal seepage field and stability change rules of high cutting slopes under different rainfall cycle conditions and after rain. And the annual degradation equation of the stability coefficient of the rainfall-type high cutting slope has been established; Combining the results of indoor experiment and numerical simulation, a method for predicting the stability coefficient of this type of high cutting slope has been established. The research results show that the shape and position of the potential sliding surface of the high cutting slope does not change significantly during the process of rainfall infiltration, and it appears as a combined sliding surface of the circular surface and the bedrock water surface at the depth of the high cutting slope; In the same rainfall period, the depth of rainfall infiltration and the extent of stability degradation of high cutting slopes are positively correlated with rainfall intensity; The degree of stability degradation of high cutting slopes in a single rainfall cycle is negatively correlated with the number of rainfall cycles. The greater the depth of rainfall infiltration or the total rainfall infiltration of high cutting slopes, the smaller the recovery of the stability coefficient of high cutting slopes after rain; It is relatively safe and reasonable to use the average value of the degradation coefficient of the shear strength parameter of the soil for the degradation coefficient of the high cutting slope stability coefficient. And the prediction method for the stability coefficient of rainfall-type high cutting slope has high accuracy.

    Keywords:rainfall-type high cutting slope; dry-wet cycle test; numerical simulation; stability; evolution process; prediction method

    近年來(lái),中國(guó)山區(qū)大型交通工程基礎(chǔ)設(shè)施大規(guī)模建設(shè),隨之而來(lái)的巖土邊坡穩(wěn)定性問(wèn)題也與日俱增。以云貴川渝為例,當(dāng)?shù)厝荷江h(huán)伺、山高谷深、地勢(shì)陡峻、地質(zhì)構(gòu)造及地形地貌條件復(fù)雜特殊,公路工程勢(shì)必依山而建,常在山腰或山腳處形成長(zhǎng)達(dá)幾公里的工程高切坡,受當(dāng)?shù)亟?jīng)濟(jì)發(fā)展、交通量需求、道路等級(jí)要求、財(cái)政投資等因素限制,絕大多數(shù)工程切坡缺乏有效的支護(hù)措施,大量無(wú)支護(hù)裸露的工程切坡為降雨型滑坡的產(chǎn)生提供了一定的物質(zhì)來(lái)源[1]。據(jù)統(tǒng)計(jì),中國(guó)每年平均發(fā)生新老滑坡次數(shù)約3×105起,其中,降雨型滑坡次數(shù)占滑坡總次數(shù)的90%以上,當(dāng)降雨條件(降雨類(lèi)型、降雨強(qiáng)度、降雨歷時(shí))達(dá)到一定閾值時(shí),該類(lèi)邊坡極易發(fā)生失穩(wěn)破壞,進(jìn)而造成巨大財(cái)產(chǎn)、經(jīng)濟(jì)損失。

    目前針對(duì)降雨型邊坡穩(wěn)定性的計(jì)算主要以降雨入滲模型為基礎(chǔ),建立穩(wěn)定性計(jì)算與評(píng)價(jià)方法,如常金源等[2]、覃小華等[3]、汪丁建等[4]、蘇永華等[5]基于Chen-Young改進(jìn)模型建立了適用于降雨入滲條件下無(wú)限長(zhǎng)邊坡的穩(wěn)定性理論計(jì)算與分析方法;Bordoni等[6]、Tsuchida等[7]從誘發(fā)滑坡的水文因素出發(fā),利用現(xiàn)場(chǎng)觀測(cè)調(diào)查數(shù)據(jù)建立了簡(jiǎn)化的降雨型滑坡穩(wěn)定性分析方法。但以上研究?jī)H建立了降雨型邊坡瞬時(shí)穩(wěn)定性計(jì)算方法,無(wú)法獲得降雨型邊坡降雨入滲過(guò)程及穩(wěn)定性演化過(guò)程,因此,近來(lái)年,降雨型邊坡穩(wěn)定性演化過(guò)程分析逐漸成為該領(lǐng)域的研究熱點(diǎn)及技術(shù)難題之一,其旨在揭示邊坡降雨入滲過(guò)程及其失穩(wěn)破壞過(guò)程,常見(jiàn)手段如采用ABSQUS、SEEP/W、FLAC3D等數(shù)值模擬軟件進(jìn)行降雨入滲模擬,如任德斌等[8]、李濤[9]、韓同春等[10]、劉杰等[11]通過(guò)有限元數(shù)值模擬軟件建立邊坡數(shù)值模型,模擬了不同工況下的邊坡降雨入滲過(guò)程,獲得了諸如孔隙水壓力、位移、穩(wěn)定系數(shù)的分布及變化規(guī)律,揭示了降雨作用下邊坡失穩(wěn)機(jī)理;部分學(xué)者從降雨入滲過(guò)程中坡體抗剪強(qiáng)度參數(shù)劣化這一角度,針對(duì)邊坡穩(wěn)定性演化過(guò)程展開(kāi)分析,如董金玉等[12]、李文等[13]、涂義亮等[14]、龍安發(fā)等[15]、朱澤勇等[16]以粉質(zhì)黏土為研究對(duì)象,通過(guò)開(kāi)展室內(nèi)人工降雨模擬干濕循環(huán)試驗(yàn),分析了干濕循環(huán)下粉質(zhì)黏土試樣力學(xué)參數(shù)演變規(guī)律及邊坡破壞機(jī)理;Griffiths等[17]、Conte等[18]、Godt等[19]通過(guò)建立理想的邊坡滑動(dòng)模型,結(jié)合實(shí)地調(diào)查建立了無(wú)限邊坡在降雨入滲過(guò)程中的入滲深度及穩(wěn)定系數(shù)評(píng)價(jià)方法。

    以上研究多立足無(wú)限長(zhǎng)傾斜式均土質(zhì)邊坡,降雨過(guò)程中,邊坡穩(wěn)定性理論計(jì)算多采用瞬時(shí)計(jì)算且基于滑面平行于坡表假定,但忽略了有限長(zhǎng)人為工程切坡具有臺(tái)階,即其坡面為非直線這一事實(shí),導(dǎo)致針對(duì)上覆土體下伏基巖的基巖型臺(tái)階狀土質(zhì)高切坡降雨入滲過(guò)程及穩(wěn)定性分析認(rèn)識(shí)不足;部分學(xué)者雖通過(guò)開(kāi)展室內(nèi)土樣的干濕循環(huán)試驗(yàn)獲得了土體抗剪強(qiáng)度參數(shù)劣化規(guī)律,并以此揭示降雨型邊坡穩(wěn)定系數(shù)演化規(guī)律,但其在邊坡穩(wěn)定性演化模擬分析中多采用單次降雨入滲工況,且受限于現(xiàn)有數(shù)值模擬手段無(wú)法輸入時(shí)變坡體抗剪強(qiáng)度參數(shù),導(dǎo)致其針對(duì)降雨循環(huán)條件下邊坡入滲過(guò)程穩(wěn)定性演化過(guò)程適用性較窄,進(jìn)而缺乏降雨型邊坡穩(wěn)定性狀態(tài)預(yù)測(cè)方法相關(guān)研究。基于此,筆者以樂(lè)西高速馬邊至昭覺(jué)段A1標(biāo)段某降雨型粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡為研究對(duì)象,通過(guò)開(kāi)展室內(nèi)干濕循環(huán)試驗(yàn),建立了土體抗剪強(qiáng)度參數(shù)劣化數(shù)學(xué)模型;利用Geo-studio數(shù)值模擬軟件,建立了高切坡二維滲流數(shù)值模型,通過(guò)輸入2017年實(shí)時(shí)降雨數(shù)據(jù)及自擬多工況降雨循環(huán)函數(shù),開(kāi)展了高切坡降雨入滲過(guò)程研究,揭示不同降雨循環(huán)工況下及雨后高切坡內(nèi)部滲流場(chǎng)及穩(wěn)定性變化規(guī)律,建立降雨型高切坡穩(wěn)定系數(shù)逐年劣化方程;基于室內(nèi)試驗(yàn)及數(shù)值模擬結(jié)果,建立該類(lèi)降雨型高切坡穩(wěn)定性預(yù)測(cè)方法,對(duì)降雨循環(huán)條件下此類(lèi)高切坡穩(wěn)定性評(píng)價(jià)、預(yù)測(cè)預(yù)報(bào)具有重要的科學(xué)指導(dǎo)意義及實(shí)用價(jià)值。

    1 研究工點(diǎn)選取

    公路工程高切坡坡體長(zhǎng)期暴露于自然環(huán)境中,在氣候周期性交替變化作用下,會(huì)經(jīng)歷“降雨”“蒸發(fā)”“降雨”的循環(huán)過(guò)程,長(zhǎng)期處在干燥和加濕的干濕循環(huán)條件下,坡體會(huì)經(jīng)歷“濕”“干”“濕”的循環(huán)過(guò)程,“干”狀態(tài)指坡體受自身排水過(guò)程及日照蒸發(fā),土體含水量逐漸減少,直至天然含水量的過(guò)程;“濕”狀態(tài)指坡體受降雨作用,土體含水量增加,直至趨向飽和狀態(tài)的過(guò)程。周期性降雨潤(rùn)濕及排水蒸發(fā)即干濕交替過(guò)程對(duì)坡體產(chǎn)生了干濕循環(huán)作用,導(dǎo)致土體力學(xué)性質(zhì)的劣化,周期性降雨作用后,該劣化效應(yīng)累積性發(fā)展,高切坡穩(wěn)定性隨之劣化,直至失穩(wěn)。作為常見(jiàn)邊坡工程病害之一,公路工程高切坡淺層失穩(wěn)溜坍土體呈厚度薄、規(guī)模小及頻率高的特點(diǎn),其中,中國(guó)南部粉質(zhì)黏土公路邊坡及中西部粉土公路邊坡中分布較為廣泛,

    筆者以樂(lè)西高速馬邊至昭覺(jué)段A1標(biāo)段K41+538~K41+615右側(cè)粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡為研究對(duì)象,見(jiàn)圖1,研究工點(diǎn)位于樂(lè)山市馬邊縣羅山溪附近,地處四川盆地中部,總體地貌類(lèi)型屬構(gòu)造侵蝕剝蝕型丘陵地貌。邊坡坡體主要為第四系全新統(tǒng)殘坡積(Q4dl+el)粉質(zhì)黏土,土層厚度6.0~13.0 m,呈軟塑~可塑狀;下部不透水基巖為侏羅系上統(tǒng)遂寧組(J2S)泥巖,多成強(qiáng)~中風(fēng)化,初勘顯示,地下水水位線遠(yuǎn)低于基巖底部。工點(diǎn)擬采用二級(jí)平臺(tái)開(kāi)挖方式,高切坡單級(jí)坡高8 m,坡率1∶1,二級(jí)邊坡交接處平臺(tái)寬2 m,坡角為45°,且不透水基巖面傾角30°。工程區(qū)屬亞熱帶濕潤(rùn)季風(fēng)氣候區(qū),年內(nèi)四季分明,雨量充沛,秋季多綿雨,多年平均降雨量1 025.8 mm,降雨多集中在5月-9月,約占全年降雨量的70%,多年平均蒸發(fā)量1 168.3 mm,蒸發(fā)量與降雨量基本持平,可認(rèn)定坡體處于飽和非飽和狀態(tài)。土體力學(xué)參數(shù)見(jiàn)表1。

    2 土體強(qiáng)度參數(shù)劣化數(shù)學(xué)模型建立

    2.1 模型建立

    為描述土體抗剪強(qiáng)度τ在周期性降雨循環(huán)過(guò)程中的變化規(guī)律,定義土體黏聚力cn與內(nèi)摩擦角φn表示土體在周期性降雨循環(huán)過(guò)程中的實(shí)時(shí)抗剪強(qiáng)度參數(shù),假定干濕循環(huán)作用過(guò)程對(duì)土體強(qiáng)度參數(shù)劣化效應(yīng)過(guò)程連續(xù)且不可逆,引入時(shí)間概念,則土體實(shí)時(shí)抗剪強(qiáng)度為

    τn=cn+σn·tan φn(1)

    式中:τn為n次干濕循環(huán)后的土體抗剪強(qiáng)度,kPa;σn為作用于土體單元的法向正應(yīng)力,kPa。

    定義土體黏聚力劣化系數(shù)為ζ,則土體實(shí)時(shí)黏聚力cn與干濕循環(huán)次數(shù)n的變化關(guān)系式為

    cn=c0·ζ? 或ζ=cn/c0(2)

    式中:c0為土體初始黏聚力,kPa。

    定義土體內(nèi)摩擦角劣化系數(shù)為η,則土體實(shí)時(shí)內(nèi)摩擦角φn與干濕循環(huán)次數(shù)n的變化關(guān)系式為[12, 14]

    φn=φ0·η或η=φn/φ0(3)

    式中:φ0為土體初始內(nèi)摩擦角,(°)。

    鄧華鋒等[20]定義干濕循環(huán)過(guò)程中的土體抗剪強(qiáng)度參數(shù)降低程度為劣化度,因此,定義n次干濕循環(huán)后土體抗剪強(qiáng)度參數(shù)劣化程度為劣化度Sn,則土體黏聚力劣化度Scn

    Scn=c0-cnc0×100%或Scn=1-ζ(4)

    內(nèi)摩擦角劣化度Sφ n

    Sφ n=φ0-φnφ0×100%或Sφ n=1-η(5)

    2.2 室內(nèi)干濕循環(huán)試驗(yàn)

    為進(jìn)一步探尋粉質(zhì)黏土坡體周期性降雨過(guò)程與室內(nèi)干濕循環(huán)過(guò)程之間的聯(lián)系,給出高切坡穩(wěn)定性預(yù)測(cè)方法,通過(guò)現(xiàn)場(chǎng)制取粉質(zhì)黏土試樣展開(kāi)室內(nèi)干濕循環(huán)試驗(yàn)及直剪試驗(yàn),以獲得土體抗剪強(qiáng)度參數(shù)干濕循環(huán)劣化規(guī)律。干濕循環(huán)試驗(yàn)共制取土體試樣數(shù)量48個(gè),擬定直剪試驗(yàn)荷載共分為0、50、100、200 kPa共4種,每種荷載下12個(gè)試樣共4組,每組平均3個(gè)試樣,進(jìn)行相應(yīng)荷載條件下的1、2、3、4次干濕循環(huán)試驗(yàn),其中,直剪試驗(yàn)土樣固結(jié)穩(wěn)定時(shí)的垂直變形值標(biāo)準(zhǔn)為每1 h不大于0.005 mm,

    土樣實(shí)際固結(jié)穩(wěn)定時(shí)間實(shí)測(cè)約24 h,剪切速率為0.8 mm/min,并在3~5 min內(nèi)剪損。干濕循環(huán)試驗(yàn)首先將天然土樣放入飽和缸中進(jìn)行抽氣,接著向缸內(nèi)注滿水,浸泡24 h增濕至土樣飽和,然后將飽和試件放入烘箱烘干,減濕至土樣天然含水率,至此,完成一次干濕循環(huán)。重復(fù)以上過(guò)程,進(jìn)行第2~4次循環(huán)試驗(yàn),試驗(yàn)過(guò)程見(jiàn)圖2~圖4。由圖3可見(jiàn),干濕循環(huán)1次時(shí),粉質(zhì)黏土試樣周邊開(kāi)始出現(xiàn)微裂隙,隨著循環(huán)次數(shù)的增加,微裂隙不斷由試樣四周向中間擴(kuò)展,發(fā)展成為貫通裂隙,干濕循環(huán)4次時(shí),試樣貫通裂隙增寬加深,土體結(jié)構(gòu)破壞嚴(yán)重,并且表現(xiàn)有一定的體積收縮,究其原因在于,干濕循環(huán)過(guò)程中,試樣中黏土礦物吸水膨脹和失水收縮,進(jìn)而土體出現(xiàn)干縮濕脹循環(huán),導(dǎo)致土體中原生裂隙的開(kāi)展和次生裂隙的發(fā)育,而產(chǎn)生的裂隙降低了黏土顆粒間的膠結(jié)作用,造成了不可逆的損傷累計(jì),最終破壞了原狀土體的結(jié)構(gòu),土體抗剪強(qiáng)度逐漸降低。所獲粉質(zhì)黏土試樣干濕循環(huán)試驗(yàn)結(jié)果見(jiàn)表2。

    粉質(zhì)黏土試樣黏聚力及內(nèi)摩擦角劣化曲線見(jiàn)圖5。由圖5可知,干濕循環(huán)周期內(nèi)粉質(zhì)黏土試樣黏聚力及內(nèi)摩擦角與干濕循環(huán)次數(shù)呈負(fù)相關(guān),試樣抗剪強(qiáng)度不斷降低,但試樣黏聚力劣化幅度較內(nèi)摩擦角劇烈,整體上看,單次干濕循環(huán)周期內(nèi)試樣黏聚力及內(nèi)摩擦角劣化幅度隨干濕循環(huán)次數(shù)增加逐漸減小,表明后續(xù)的干濕循環(huán)試驗(yàn)對(duì)粉質(zhì)黏土試樣黏聚力及內(nèi)摩擦角的劣化效應(yīng)逐漸減小。

    3 Geo-studio降雨入滲數(shù)值模擬

    3.1 模型建立

    以樂(lè)西高速馬邊至昭覺(jué)段A1標(biāo)段K41+538~K41+615右側(cè)粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡為研究對(duì)象,利用Geo-studio數(shù)值模擬軟件建立該處高切坡二維滲流數(shù)值模型(圖6)。模型邊坡坡頂長(zhǎng)度設(shè)置為24 m(1.5 H′,H′為2個(gè)單級(jí)坡高總和),左邊界至坡腳的距離設(shè)置為24 m(1.5 H′),坡頂部到底部邊界的距離設(shè)置為32 m(2 H′的方式進(jìn)行擴(kuò)大),且不透水基巖面傾角為30°,最終有限元模型尺寸設(shè)計(jì)如圖6所示。模型底部邊界設(shè)定為豎直約束和水平約束邊界;左右邊界為水平約束邊界;高切坡坡表及坡腳左側(cè)路面設(shè)置為降雨入滲及潛在的排水邊界,為考慮雨后高切坡穩(wěn)定性回升,同時(shí)還將坡表設(shè)置為蒸發(fā)輻射邊界;考慮計(jì)算精度與計(jì)算時(shí)步,將模型劃分為6 492個(gè)節(jié)點(diǎn)和6 321個(gè)單元。

    模型邊坡坡體為均質(zhì)非飽和粉質(zhì)黏土,根據(jù)GENUCHTEN研究發(fā)現(xiàn),其滲透系數(shù)為基質(zhì)力的函數(shù),而非常數(shù),即基質(zhì)力隨含水量的變化而變化。土體非飽和特性參數(shù)如表1所示,利用Geo-studio數(shù)值模擬軟件內(nèi)置的非飽和粉質(zhì)黏土樣本函數(shù)及實(shí)驗(yàn)測(cè)定的飽和含水率(θs=0.46)和殘余含水率(θr=0.08)數(shù)據(jù),生成覆蓋層粉質(zhì)黏土的土水特征曲線及滲透系數(shù)函數(shù),見(jiàn)圖7。

    3.2 模擬工況選定

    3.2.1 實(shí)時(shí)降雨工況

    1)降雨工況

    收集了中國(guó)氣象數(shù)據(jù)網(wǎng)關(guān)于樂(lè)山市馬邊縣2017年的降雨資料,將高切坡降雨入滲模擬過(guò)程周期定為1個(gè)單位年,由于研究工點(diǎn)全年降雨量主要集中在6月—8月,僅給出研究工點(diǎn)處6月—8月的降雨數(shù)據(jù),見(jiàn)圖8,其余月份降雨量反映在圖9中。

    其中,研究工點(diǎn)處近10年來(lái)全年水汽蒸發(fā)及輻射強(qiáng)度約為3 340~4 190 MJ/m2·a,考慮蒸發(fā)量與降雨量基本持平,6月—8月模型輻射強(qiáng)度占全年輻射強(qiáng)度的70%,其余月份占30%,故將6月—8月模型邊界條件每日輻射強(qiáng)度設(shè)為28 339 kJ/m2·d,其余月份設(shè)為4 183 kJ/m2·d;空氣濕度對(duì)邊坡穩(wěn)定性分析影響不大,6月—8月設(shè)為30%,其余月份設(shè)為20%;考慮6月—8月溫度較高,分別設(shè)為26、29、28 ℃;各類(lèi)型降雨均取國(guó)家氣象局規(guī)定的降雨量標(biāo)準(zhǔn)表中相應(yīng)降雨類(lèi)型的平均值[21],即小雨、中雨、大雨、暴雨分別取5、17.5、37.5、75 mm/d。

    2)穩(wěn)定性分析

    利用Geo-studio數(shù)值模擬軟件瞬時(shí)分析計(jì)算邊坡在2017年降雨周期內(nèi)各時(shí)刻邊坡的降雨入滲過(guò)程,并獲得該高切坡穩(wěn)定系數(shù)時(shí)程曲線,見(jiàn)圖10。

    由圖10可見(jiàn),該高切坡2017年全年穩(wěn)定狀態(tài)整體處于波動(dòng)狀態(tài),全年穩(wěn)定系數(shù)劣化0.053,高切坡最有利狀態(tài)位于6月中旬,原因在于該階段內(nèi)研究工點(diǎn)處降雨較少,高切坡經(jīng)歷長(zhǎng)時(shí)間太陽(yáng)輻射及排水,坡體內(nèi)部孔隙水壓力逐漸消散,邊坡穩(wěn)定性較好;邊坡最不利狀態(tài)位于8月上旬,原因在于該階段內(nèi)研究工點(diǎn)降雨較多,高切坡經(jīng)歷長(zhǎng)時(shí)間降雨入滲,坡體軟化,坡體自重加大,孔隙水壓力逐漸增大,邊坡穩(wěn)定性較差,兩者相差0.256,可見(jiàn),雨季高切坡穩(wěn)定狀態(tài)劣化較大,應(yīng)及時(shí)加強(qiáng)邊坡監(jiān)測(cè),做好疏排水措施。

    考慮研究工點(diǎn)處若干年內(nèi)降雨條件無(wú)較大差別,則該高切坡穩(wěn)定系數(shù)逐年劣化方程為

    FS=2.451-0.053a(6)

    式中:a為考察年限,按年計(jì)。該式可推廣為

    FS=F0-ΔF·a(7)

    式中:F0為考察年限內(nèi)高切坡穩(wěn)定系數(shù)初始值;ΔF為單年內(nèi)高切坡穩(wěn)定系數(shù)劣化值。

    參照《建筑邊坡工程技術(shù)規(guī)范》(GB 50330—2013)[22]對(duì)永久邊坡安全系數(shù)的規(guī)定,當(dāng)該邊坡工程安全等級(jí)分別為一級(jí)、二級(jí)及三級(jí)時(shí),其穩(wěn)定系數(shù)應(yīng)大于1.35、1.30及1.25,在不加任何支護(hù)措施的前提下,由式(6)可知,研究對(duì)象在常年周期性降雨下其穩(wěn)定系數(shù)將分別在第20.77年、21.72年、22.66年劣化至安全系數(shù)臨界值,在第27.38年時(shí)劣化至失穩(wěn)破壞臨界值。

    3.2.2 典型降雨工況

    1)降雨工況擬定

    為進(jìn)一步揭示粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡在不同降雨循環(huán)作用下其穩(wěn)定性變化規(guī)律,考慮粉質(zhì)黏土坡體低滲透性,初擬3種典型雨型進(jìn)行坡體降雨入滲特征分析及穩(wěn)定性分析,其中降雨循環(huán)結(jié)束后觀察期設(shè)為4 d,期間模型邊界每日輻射強(qiáng)度設(shè)為28 339 kJ/m2·d。滲流計(jì)算過(guò)程中,根據(jù)所擬降雨工況設(shè)置90個(gè)時(shí)間步,每個(gè)時(shí)間步長(zhǎng)時(shí)長(zhǎng)為8 h,計(jì)算過(guò)程共計(jì)30 d,每個(gè)時(shí)間步最大迭代次數(shù)為50步,迭代誤差小于1%視為收斂,以揭示降雨循環(huán)入滲過(guò)程及降雨循環(huán)后高切坡內(nèi)部滲流場(chǎng)及穩(wěn)定性變化規(guī)律,見(jiàn)圖11。

    2)穩(wěn)定性計(jì)算與分析

    各工況降雨周期內(nèi)高切坡內(nèi)部濕潤(rùn)峰動(dòng)態(tài)變化過(guò)程如圖12~圖14所示,可見(jiàn),整個(gè)降雨周期內(nèi),Geo-studio軟件自動(dòng)搜索高切坡潛在滑移面位置及形狀并未發(fā)生明顯變化,邊坡潛在滑移面形態(tài)無(wú)明顯的尺寸效應(yīng),均表現(xiàn)為邊坡深處圓弧面和基巖積水面的組合型滑面;相同降雨時(shí)間內(nèi)高切坡降雨入滲深度zf與降雨強(qiáng)度成正相關(guān);降雨結(jié)束后,觀察期內(nèi)高切坡入滲深度逐漸減少,究其原因在于該階段內(nèi)模型邊界施加了蒸發(fā)與輻射條件,導(dǎo)致坡體孔隙水壓消散加速,土體有效應(yīng)力提高,高切坡穩(wěn)定性逐漸提高。

    各降雨工況高切坡穩(wěn)定性演化過(guò)程見(jiàn)圖15。由圖15可見(jiàn),降雨強(qiáng)度一定時(shí),隨著降雨時(shí)間的增大,邊坡穩(wěn)定性逐漸下降;工況1、工況2及工況3降雨周期內(nèi)導(dǎo)致高切坡穩(wěn)定系數(shù)累計(jì)分別下降0.502、0.602及0.710,表明相同降雨時(shí)間內(nèi),高切坡穩(wěn)定性劣化幅度與降雨強(qiáng)度成正相關(guān)。究其原因,相同時(shí)間內(nèi)降雨強(qiáng)度越大,降雨總量越大,覆蓋層坡體飽和度、孔隙水壓力、飽和區(qū)域越大,降雨入滲越嚴(yán)重,高切坡穩(wěn)定性下降越劇烈;降雨循環(huán)次數(shù)越多,相同降雨時(shí)間內(nèi)高切坡穩(wěn)定系數(shù)下降幅度越小,見(jiàn)圖16,該現(xiàn)象與上述室內(nèi)粉質(zhì)黏土抗剪強(qiáng)度參數(shù)劣化規(guī)律基本吻合,表明單次降雨循環(huán)周期內(nèi)高切坡穩(wěn)定性劣化幅度與降雨循環(huán)次數(shù)成負(fù)相關(guān),高切坡最終會(huì)趨于穩(wěn)定。各降雨工況結(jié)束后,相同觀察期內(nèi),高切坡穩(wěn)定系數(shù)均有回升,工況1、工況2及工況3穩(wěn)定系數(shù)分別回升0.027、0.024及0.019,回升幅度分別為1.5%、1.3%及1%,表明工況1回升幅度

    最大,工況2次之,工況3最小,究其原因在于工況1降雨入滲總量最小,雨后孔壓水壓消散最快,邊坡排水過(guò)程較快,非飽和區(qū)域增大較快,高切坡穩(wěn)定性恢復(fù)最快,相同觀察期內(nèi)高切坡穩(wěn)定系數(shù)回升幅度最大。

    4 高切坡穩(wěn)定系數(shù)預(yù)測(cè)方法

    4.1 室內(nèi)試驗(yàn)法

    結(jié)合室內(nèi)干濕循環(huán)試驗(yàn),將單次干濕循環(huán)后的粉質(zhì)黏土試樣抗剪強(qiáng)度參數(shù)視為單次降雨循環(huán)后高切坡坡體抗剪強(qiáng)度參數(shù)初始值,并計(jì)算高切坡穩(wěn)定系數(shù),見(jiàn)圖17。

    由圖17可見(jiàn),由室內(nèi)干濕循環(huán)試驗(yàn)所獲抗剪強(qiáng)度參數(shù)計(jì)算得到的高切坡穩(wěn)定系數(shù)與Geo-studio軟件模擬工況2循環(huán)降雨下的高切坡穩(wěn)定系數(shù)基本吻合,當(dāng)循環(huán)次數(shù)小于2次時(shí),工況1所示降雨循環(huán)試驗(yàn)高切坡降雨強(qiáng)度較小,坡體內(nèi)濕潤(rùn)鋒下滲較淺,無(wú)法較真實(shí)地模擬室內(nèi)干濕循環(huán)試驗(yàn)“濕”狀態(tài),導(dǎo)致高切坡穩(wěn)定系數(shù)顯著大于室內(nèi)干濕循環(huán)試驗(yàn);當(dāng)循環(huán)次數(shù)小于2次時(shí),工況3降雨循環(huán)周期內(nèi)由于高切坡初始降雨強(qiáng)度較大,其較好地模擬了坡體飽和天然狀態(tài),故兩者吻合較好,當(dāng)循環(huán)次數(shù)大于2次時(shí),

    由于坡體排水較慢,孔隙水未完全消散,坡體未達(dá)到“干”狀態(tài),而室內(nèi)干濕循環(huán)試驗(yàn)可保證試樣處于天然狀態(tài),由室內(nèi)干濕循環(huán)試驗(yàn)所得高切坡穩(wěn)定系數(shù)高于降雨循環(huán)試驗(yàn),故工況2較好地代表了室內(nèi)干濕循環(huán)試驗(yàn),其中,當(dāng)循環(huán)次數(shù)分別為1、2、3及4時(shí),兩者分別相差0.009、0.011、0.022及0.031,表明隨著循環(huán)次數(shù)的增加,兩者差值逐漸增大,究其原因,當(dāng)坡體趨于完全飽和時(shí),在為時(shí)4 d的非降雨時(shí)段內(nèi),輻射強(qiáng)度28 339 kJ/m2·d坡體排水未完全,孔隙水未完全消散,見(jiàn)圖14(d)、(e)、(f),但該誤差在工程范圍內(nèi)可以接受。

    定義高切坡穩(wěn)定系數(shù)劣化系數(shù)ω表示在周期性降雨循環(huán)過(guò)程中其穩(wěn)定性劣化規(guī)律,進(jìn)而引入時(shí)間概念表征隨降雨循環(huán)次數(shù)變化的高切坡穩(wěn)定系數(shù)為FSn

    FSn=F0·ω(8)

    由圖18可發(fā)現(xiàn),所研究高切坡穩(wěn)定系數(shù)劣化規(guī)律與粉質(zhì)黏土土樣cn及φn平均值劣化規(guī)律一致,且劣化程度雖有誤差但誤差范圍內(nèi)高切坡穩(wěn)定系數(shù)整體基本吻合,高切坡穩(wěn)定系數(shù)劣化規(guī)律描述采用坡體cn及φn劣化系數(shù)平均值量化大于數(shù)值模擬劣化結(jié)果,用該值偏于安全,在工程范圍內(nèi)可以接受。

    因此,針對(duì)所研究粉質(zhì)黏土覆蓋層高切坡,式(8)可轉(zhuǎn)換為

    FSn=F0·ζ+η2(9)

    4.2 數(shù)值模擬擬合法

    對(duì)于某一特定高切坡,降雨入滲過(guò)程中,決定高切坡穩(wěn)定性的主要因素為降雨時(shí)間T及降雨強(qiáng)度q,基于此,可獲得高切坡在不同降雨時(shí)間及降雨強(qiáng)度作用后的穩(wěn)定性情況,利用MATLAB 2016a數(shù)據(jù)分析軟件對(duì)模擬數(shù)據(jù)進(jìn)行擬合,其中,為更準(zhǔn)確地描述降雨強(qiáng)度q隨時(shí)間的變化趨勢(shì),仍采用典型雨型函數(shù)并用降雨總量Q刻畫(huà)降雨時(shí)間T內(nèi)的降雨強(qiáng)度,可以發(fā)現(xiàn),針對(duì)所研究粉質(zhì)黏土覆蓋層高切坡,高切坡降雨實(shí)時(shí)穩(wěn)定系數(shù)FS與降雨總量Q、降雨時(shí)間T滿足如下二元函數(shù)關(guān)系,如圖19所示。

    FS=2.389-0.015T-0.001Q(10)

    為進(jìn)一步探究降雨總量及降雨時(shí)間對(duì)高切坡穩(wěn)定性的影響,利用多元回歸理論對(duì)兩變量進(jìn)行SPSS擬合,得到回歸模型參數(shù),見(jiàn)表5。結(jié)果顯示,回歸模型系數(shù)同MATLAB 2016a所得模型系數(shù)一致,降雨總量Q、降雨歷時(shí)T變量均通過(guò)了t檢驗(yàn),且兩變量方差膨脹系數(shù)VIF均小于5,表明兩變量不具有共線性,故兩變量均能引入線性回歸模型,且模型相關(guān)系數(shù)R2=0.976 2,證明模型擬合度較高。

    圖20為模型標(biāo)準(zhǔn)化殘差P-P圖,由圖20可見(jiàn),各點(diǎn)近似一條直線且均分布于直線附近;圖21為模型回歸殘差分布直方圖,由圖21可見(jiàn),直方圖形態(tài)與標(biāo)準(zhǔn)正態(tài)分布近似,即表明模型的殘差符合正態(tài)分布假設(shè)。

    綜上,對(duì)于所建粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡,認(rèn)定其降雨實(shí)時(shí)穩(wěn)定系數(shù)FS與降雨總量Q、降雨時(shí)間T滿足式(10)的函數(shù)關(guān)系是合理的。

    5 結(jié)論

    以降雨型粉質(zhì)黏土覆蓋層基巖型臺(tái)階狀高切坡為研究對(duì)象,通過(guò)室內(nèi)干濕循環(huán)試驗(yàn)建立了土體抗剪強(qiáng)度參數(shù)劣化數(shù)學(xué)模型;利用Geo-studio軟件揭示了其在不同降雨工況下的降雨入滲過(guò)程及穩(wěn)定性演化過(guò)程,并給出了該類(lèi)降雨型高切坡的穩(wěn)定性預(yù)測(cè)方法,主要結(jié)論如下:

    1)通過(guò)室內(nèi)干濕循環(huán)試驗(yàn)及直剪試驗(yàn)獲得了土體抗剪強(qiáng)度參數(shù)干濕循環(huán)劣化規(guī)律,建立了土體抗剪強(qiáng)度參數(shù)劣化數(shù)學(xué)模型。

    2)以高切坡2017年實(shí)時(shí)降雨數(shù)據(jù)為例,利用Geo-studio數(shù)值模擬軟件獲得了高切坡在該年內(nèi)的穩(wěn)定性演化過(guò)程,建立了高切坡穩(wěn)定系數(shù)逐年劣化方程;通過(guò)自擬3種典型雨型函數(shù),利用Geo-studio數(shù)值模擬軟件開(kāi)展了高切坡在不同降雨強(qiáng)度、不同降雨歷時(shí)下的降雨入滲特征及穩(wěn)定性分析,結(jié)果表明:整個(gè)降雨周期內(nèi),高切坡潛在滑移面為高切坡深處圓弧面和基巖積水面的組合型滑面;相同降雨時(shí)間內(nèi),高切坡降雨入滲深度及穩(wěn)定性劣化幅度與降雨強(qiáng)度成正相關(guān);單次降雨循環(huán)周期內(nèi),高切坡穩(wěn)定性劣化幅度與降雨循環(huán)次數(shù)成負(fù)相關(guān);高切坡降雨入滲深度或入滲總量越大,雨后高切坡穩(wěn)定系數(shù)回升越小。

    3)針對(duì)所研究粉質(zhì)黏土覆蓋層高切坡,所擬工況2下高切坡穩(wěn)定系數(shù)與以室內(nèi)干濕循環(huán)試驗(yàn)所獲抗剪強(qiáng)度參數(shù)計(jì)算得到的高切坡穩(wěn)定系數(shù)劣化規(guī)律基本吻合,降雨循環(huán)作用下,高切坡穩(wěn)定系數(shù)劣化系數(shù)采用土樣cn及φn平均值是偏安全合理的;基于Geo-studio數(shù)值模擬結(jié)果給出了該類(lèi)高切坡在不同降雨量及降雨時(shí)間下的穩(wěn)定性預(yù)測(cè)方法,所獲結(jié)果精度較高。參考文獻(xiàn):

    [1] 陳洪凱, 楊世勝, 唐紅梅. 高切坡研究現(xiàn)狀及趨勢(shì)分析[C]//公路邊坡及其環(huán)境工程技術(shù)交流會(huì), 2005.

    CHEN H K, YANG S S, TANG H M. Research status and trend of high-cut slope [C]//Highway Slope and Environmental Engineering Technology Exchange Conference, 2005. (in Chinese)

    [2] 常金源, 包含, 伍法權(quán), 等. 降雨條件下淺層滑坡穩(wěn)定性探討[J]. 巖土力學(xué), 2015, 36(4): 995-1001.

    CHANG J Y, BAO H, WU F Q, et al. Discussion on stability of shallow landslide under rainfall [J]. Rock and Soil Mechanics, 2015, 36(4): 995-1001. (in Chinese)

    [3] 覃小華, 劉東升, 宋強(qiáng)輝, 等. 強(qiáng)降雨條件下基巖型層狀邊坡入滲模型及穩(wěn)定性研究[J]. 巖土力學(xué), 2016, 37(11): 3156-3164.

    QIN X H, LIU D S, SONG Q H, et al. Infiltration model of bedrock laminar slope under heavy rainfall and its stability analysis [J]. Rock and Soil Mechanics, 2016, 37(11): 3156-3164. (in Chinese)

    [4] 汪丁建, 唐輝明, 李長(zhǎng)冬, 等. 強(qiáng)降雨作用下堆積層滑坡穩(wěn)定性分析[J]. 巖土力學(xué), 2016, 37(2): 439-445.

    WANG D J, TANG H M, LI C D, et al. Stability analysis of colluvial landslide due to heavy rainfall [J]. Rock and Soil Mechanics, 2016, 37(2): 439-445. (in Chinese)

    [5] 蘇永華, 李誠(chéng)誠(chéng). 強(qiáng)降雨下基于Green-Ampt模型的邊坡穩(wěn)定性分析[J]. 巖土力學(xué), 2020, 41(2): 389-398.

    SU Y H, LI C C. Stability analysis of slope based on Green-Ampt model under heavy rainfall [J]. Rock and Soil Mechanics, 2020, 41(2): 389-398.(in Chinese)

    [6] ?BORDONI M, MEISINA C, VALENTINO R, et al. Hydrological factors affecting rainfall-induced shallow landslides: From the field monitoring to a simplified slope stability analysis [J]. Engineering Geology, 2015, 193: 19-37.

    [7] ?TSUCHIDA T, ATHAPATHTHU A M R G, HANAOKA T, et al. Investigation of landslide calamity due to torrential rainfall in Shobara City, Japan [J]. Soils and Foundations, 2015, 55(5): 1305-1317.

    [8] 任德斌, 汪瑩, 于世海. 邊坡在降雨及地震作用下的穩(wěn)定性分析[J]. 沈陽(yáng)建筑大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 33(3): 439-446.

    REN D B, WANG Y, YU S H. Analysis on stability of gravel slope under rainfall and earthquake [J]. Journal of Shenyang Jianzhu University (Natural Science), 2017, 33(3): 439-446.(in Chinese)

    [9] 李濤. 考慮降雨及開(kāi)挖影響下的厚覆蓋層邊坡滲流特征及穩(wěn)定性[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 47(5): 1708-1714.

    LI T. Seepage characteristics and stability of overburden slope considering rainfall and excavation [J]. Journal of Central South University (Science and Technology), 2016, 47(5): 1708-1714.(in Chinese)

    [10] 韓同春, 黃福明. 雙層結(jié)構(gòu)土質(zhì)邊坡降雨入滲過(guò)程及穩(wěn)定性分析[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2012, 46(1): 39-45.

    HAN T C, HUANG F M. Rainfall infiltration process and stability analysis of two-layered slope [J]. Journal of Zhejiang University (Engineering Science), 2012, 46(1): 39-45.(in Chinese)

    [11] 劉杰, 曾鈴, 付宏淵, 等. 土質(zhì)邊坡降雨入滲深度及飽和區(qū)變化規(guī)律[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 50(2): 452-459.

    LIU J, ZENG L, FU H Y, et al. Variation law of rainfall infiltration depth and saturation zone of soil slope [J]. Journal of Central South University (Science and Technology), 2019, 50(2): 452-459.(in Chinese)

    [12] 董金玉, 趙志強(qiáng), 楊繼紅, 等. 干濕循環(huán)作用下滑帶土的變形演化和強(qiáng)度參數(shù)弱化試驗(yàn)研究[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2016, 48(Sup2): 1-7.

    DONG J Y, ZHAO Z Q, YANG J H, et al. Research on the deformation evolution and the strength weaking of slip soil under wetting-drying cycle [J]. Journal of Sichuan University (Engineering Science Edition), 2016, 48(Sup2): 1-7.(in Chinese)

    [13] 李文, 曾勝, 趙健, 等. 干濕循環(huán)作用下長(zhǎng)沙繞城高速公路粉質(zhì)黏土的損傷特性[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 48(5): 1360-1366.

    LI W, ZENG S, ZHAO J, et al. Damage characteristics of Changsha ring expressway silty clay during dry-wet cycle [J]. Journal of Central South University (Science and Technology), 2017, 48(5): 1360-1366.(in Chinese)

    [14] 涂義亮, 劉新榮, 鐘祖良, 等. 干濕循環(huán)下粉質(zhì)黏土強(qiáng)度及變形特性試驗(yàn)研究[J]. 巖土力學(xué), 2017, 38(12): 3581-3589.

    TU Y L, LIU X R, ZHONG Z L, et al. Experimental study on strength and deformation characteristics of silty clay during wetting-drying cycles [J]. Rock and Soil Mechanics, 2017, 38(12): 3581-3589.(in Chinese)

    [15] 龍安發(fā), 陳開(kāi)圣, 季永新. 不同降雨強(qiáng)度下紅黏土邊坡干濕循環(huán)試驗(yàn)研究[J]. 巖土工程學(xué)報(bào), 2019, 41(Sup2): 193-196.

    LONG A F, CHEN K S, JI Y X. Experimental study on wetting-drying cycles of red clay slopes under different rainfall intensities [J]. Chinese Journal of Geotechnical Engineering, 2019, 41(Sup2): 193-196. (in Chinese)

    [16] 朱澤勇, 賀桂成, 李豐雄, 等. 干濕交替條件下紅土邊坡破壞機(jī)理試驗(yàn)研究[J]. 長(zhǎng)江科學(xué)院院報(bào), 2018, 35(2): 73-77.

    ZHU Z Y, HE G C, LI F X, et al. Experimental investigation on the failure mechanism of red soil slope under cyclic wetting and drying [J]. Journal of Yangtze River Scientific Research Institute, 2018, 35(2): 73-77. (in Chinese)

    [17] ?GRIFFITHS D V, HUANG J S, DEWOLFE G F. Numerical and analytical observations on long and infinite slopes [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2011, 35(5): 569-585.

    [18] ?CONTE E, TRONCONE A. Simplified approach for the analysis of rainfall-induced shallow landslides [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(3): 398-406.

    [19] ?GODT J W,ENER-KAYA B, LU N, et al. Stability of infinite slopes under transient partially saturated seepage conditions [J]. Water Resources Research, 2012, 48(5): 5505-5518.

    [20] ?鄧華鋒, 肖瑤, 方景成, 等. 干濕循環(huán)作用下岸坡消落帶土體抗剪強(qiáng)度劣化規(guī)律及其對(duì)岸坡穩(wěn)定性影響研究[J]. 巖土力學(xué), 2017, 38(9): 2629-2638.

    DENG H F, XIAO Y, FANG J C, et al. Shear strength degradation and slope stability of soils at hydro-fluctuation belt of river bank slope during drying-wetting cycle [J]. Rock and Soil Mechanics, 2017, 38(9): 2629-2638.(in Chinese)

    [21] ?孔鋒, 王一飛, 方建, 等. 中等排放情景下中國(guó)未來(lái)不同強(qiáng)度降雨變化及其對(duì)總降雨貢獻(xiàn)的預(yù)估(2006—2100年)[J]. 水利水電技術(shù), 2017, 48(12): 14-21, 40.

    KONG F, WANG Y F, FANG J, et al. Pre-estimation on variations of rainfalls with different intensities and their contributions to total rainfall under emission scenarios in China from 2006 to 2100 [J]. Water Resources and Hydropower Engineering, 2017, 48(12): 14-21, 40.(in Chinese)

    [22] ?建筑邊坡工程技術(shù)規(guī)范: GB 50330—2013[S]. 北京: 中國(guó)建筑工業(yè)出版社, 2014.

    Technical code for building slope engineering: GB 50330-2013 [S]. Beijing: China Architecture & Building Press, 2014. (in Chinese)

    (編輯 王秀玲)

    猜你喜歡
    預(yù)測(cè)方法演化過(guò)程數(shù)值模擬
    生命的演化過(guò)程
    模因論視角下韓語(yǔ)“??”表“喝”動(dòng)作演化過(guò)程研究
    時(shí)間非齊次二態(tài)量子游蕩的演化過(guò)程分析
    重慶萬(wàn)盛石林的形成時(shí)代及發(fā)育演化過(guò)程
    風(fēng)電短期發(fā)電功率預(yù)測(cè)方法探討
    臨汾區(qū)塊煤系地層天然氣水合物危害與預(yù)防措施
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    水驅(qū)油效率影響因素研究進(jìn)展
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    精品久久国产蜜桃| 99精国产麻豆久久婷婷| 一区二区三区精品91| 女人精品久久久久毛片| 久久 成人 亚洲| 久久精品夜色国产| 国产淫语在线视频| 精品国产露脸久久av麻豆| 亚洲熟女精品中文字幕| 欧美精品av麻豆av| 狂野欧美激情性bbbbbb| 日韩精品免费视频一区二区三区 | 最新的欧美精品一区二区| www.色视频.com| 日本爱情动作片www.在线观看| 亚洲国产色片| 国产亚洲av片在线观看秒播厂| 人人妻人人爽人人添夜夜欢视频| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 99久久综合免费| 国产精品久久久久久精品电影小说| 在线观看免费高清a一片| 波多野结衣一区麻豆| 亚洲精品aⅴ在线观看| 亚洲av国产av综合av卡| 日韩中字成人| 99香蕉大伊视频| 精品一区二区三区四区五区乱码 | 秋霞在线观看毛片| 99热这里只有是精品在线观看| 51国产日韩欧美| 美国免费a级毛片| 国产男女内射视频| 只有这里有精品99| 日韩人妻精品一区2区三区| 久久免费观看电影| 欧美激情 高清一区二区三区| 免费少妇av软件| 97精品久久久久久久久久精品| 色婷婷久久久亚洲欧美| 国产片内射在线| 亚洲精品,欧美精品| 午夜福利影视在线免费观看| 我要看黄色一级片免费的| 人妻人人澡人人爽人人| 亚洲综合精品二区| 老司机影院成人| 黄色毛片三级朝国网站| 日本爱情动作片www.在线观看| 18禁在线无遮挡免费观看视频| 国产精品嫩草影院av在线观看| 飞空精品影院首页| 夜夜爽夜夜爽视频| 丁香六月天网| 中文字幕制服av| 国产免费视频播放在线视频| 亚洲成国产人片在线观看| 午夜福利在线观看免费完整高清在| 久久久欧美国产精品| 久久久久网色| 91久久精品国产一区二区三区| 91aial.com中文字幕在线观看| 亚洲国产日韩一区二区| 久久婷婷青草| 国产午夜精品一二区理论片| 秋霞在线观看毛片| 9191精品国产免费久久| 校园人妻丝袜中文字幕| 一区二区av电影网| 国产av码专区亚洲av| 乱人伦中国视频| 如何舔出高潮| 最后的刺客免费高清国语| 亚洲美女视频黄频| 欧美xxⅹ黑人| 免费av不卡在线播放| 午夜福利乱码中文字幕| 一二三四在线观看免费中文在 | 男人爽女人下面视频在线观看| 少妇猛男粗大的猛烈进出视频| 在现免费观看毛片| 成人漫画全彩无遮挡| 1024视频免费在线观看| 国产淫语在线视频| 91精品三级在线观看| 熟女电影av网| 视频在线观看一区二区三区| 亚洲国产av影院在线观看| 少妇人妻 视频| 亚洲第一区二区三区不卡| 老女人水多毛片| 国国产精品蜜臀av免费| 在线观看免费日韩欧美大片| 制服诱惑二区| 欧美人与性动交α欧美精品济南到 | 91午夜精品亚洲一区二区三区| freevideosex欧美| 欧美bdsm另类| 国产亚洲欧美精品永久| 国产综合精华液| 亚洲精品国产av成人精品| 51国产日韩欧美| 久久精品国产鲁丝片午夜精品| 亚洲精品久久成人aⅴ小说| 国产在线免费精品| 国产女主播在线喷水免费视频网站| 亚洲天堂av无毛| av天堂久久9| 麻豆精品久久久久久蜜桃| 免费在线观看完整版高清| 国产精品无大码| 青青草视频在线视频观看| 亚洲av综合色区一区| 在线天堂中文资源库| 国产精品无大码| 亚洲 欧美一区二区三区| 如何舔出高潮| 肉色欧美久久久久久久蜜桃| √禁漫天堂资源中文www| 日韩 亚洲 欧美在线| 两个人免费观看高清视频| 1024视频免费在线观看| 欧美xxⅹ黑人| 国产视频首页在线观看| 久久精品久久久久久久性| 亚洲性久久影院| 美女主播在线视频| 亚洲综合精品二区| 18禁国产床啪视频网站| 草草在线视频免费看| 咕卡用的链子| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕| 亚洲情色 制服丝袜| 国产日韩欧美视频二区| 99re6热这里在线精品视频| 欧美日韩视频精品一区| 国产精品久久久久久av不卡| 男女无遮挡免费网站观看| 国产女主播在线喷水免费视频网站| 97在线视频观看| 一边摸一边做爽爽视频免费| 欧美日本中文国产一区发布| 女人精品久久久久毛片| 国产日韩欧美视频二区| 欧美激情极品国产一区二区三区 | 日日啪夜夜爽| 视频在线观看一区二区三区| 亚洲高清免费不卡视频| 亚洲第一av免费看| 国产极品天堂在线| 毛片一级片免费看久久久久| 男女无遮挡免费网站观看| freevideosex欧美| 成人无遮挡网站| 久久久久久久亚洲中文字幕| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| 十八禁高潮呻吟视频| 久久人妻熟女aⅴ| 日韩中文字幕视频在线看片| 久久99热6这里只有精品| 久久久久精品人妻al黑| 久久久久视频综合| 成人二区视频| 免费av不卡在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 狂野欧美激情性bbbbbb| 久久久久久久久久久免费av| 久久久国产一区二区| 在线免费观看不下载黄p国产| 99国产精品免费福利视频| 日韩精品有码人妻一区| 国内精品宾馆在线| 热99久久久久精品小说推荐| 亚洲国产av新网站| 日韩一区二区三区影片| 中文字幕另类日韩欧美亚洲嫩草| 国国产精品蜜臀av免费| 亚洲国产精品国产精品| 婷婷色麻豆天堂久久| 最近中文字幕高清免费大全6| 韩国精品一区二区三区 | 建设人人有责人人尽责人人享有的| 不卡视频在线观看欧美| 亚洲综合色网址| 99热全是精品| 精品国产一区二区三区久久久樱花| 日韩av免费高清视频| 性色av一级| 在线天堂中文资源库| 99久久中文字幕三级久久日本| 国产 精品1| 91成人精品电影| 成年人免费黄色播放视频| 亚洲精品日韩在线中文字幕| 精品午夜福利在线看| 久久97久久精品| 一二三四在线观看免费中文在 | 欧美xxⅹ黑人| 黑人高潮一二区| 国产1区2区3区精品| 国产乱来视频区| 国产xxxxx性猛交| 亚洲成人av在线免费| av免费观看日本| av一本久久久久| 性色avwww在线观看| 天天躁夜夜躁狠狠躁躁| 国产日韩欧美亚洲二区| 中文字幕av电影在线播放| 看十八女毛片水多多多| 岛国毛片在线播放| 日本av免费视频播放| 夜夜骑夜夜射夜夜干| 日日爽夜夜爽网站| 欧美3d第一页| 国产色婷婷99| 日韩熟女老妇一区二区性免费视频| 永久网站在线| 精品午夜福利在线看| 亚洲国产最新在线播放| 亚洲精品日本国产第一区| 中文字幕免费在线视频6| 久久精品国产亚洲av天美| 国产探花极品一区二区| 一级,二级,三级黄色视频| 中文精品一卡2卡3卡4更新| 熟女人妻精品中文字幕| tube8黄色片| 女人久久www免费人成看片| 观看av在线不卡| 校园人妻丝袜中文字幕| 99久久人妻综合| 大陆偷拍与自拍| 大片电影免费在线观看免费| 亚洲人与动物交配视频| 中国美白少妇内射xxxbb| 亚洲美女黄色视频免费看| 黄色毛片三级朝国网站| 精品亚洲成国产av| 午夜福利,免费看| 22中文网久久字幕| 一级片'在线观看视频| 人妻人人澡人人爽人人| 乱人伦中国视频| 22中文网久久字幕| 91久久精品国产一区二区三区| 亚洲欧美日韩另类电影网站| 国产午夜精品一二区理论片| 男的添女的下面高潮视频| 草草在线视频免费看| 久久精品久久久久久噜噜老黄| 国产亚洲午夜精品一区二区久久| 一级毛片我不卡| 一本大道久久a久久精品| a 毛片基地| 日本欧美视频一区| 美女中出高潮动态图| 亚洲久久久国产精品| 在线看a的网站| 国产成人免费观看mmmm| 成人无遮挡网站| 黑人巨大精品欧美一区二区蜜桃 | 亚洲av男天堂| 曰老女人黄片| 成人午夜精彩视频在线观看| 亚洲av欧美aⅴ国产| 国产女主播在线喷水免费视频网站| 女人精品久久久久毛片| 精品一区二区三区视频在线| 精品少妇黑人巨大在线播放| 国产麻豆69| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| av在线老鸭窝| 汤姆久久久久久久影院中文字幕| 久久精品久久精品一区二区三区| 国产深夜福利视频在线观看| 男女下面插进去视频免费观看 | 两个人免费观看高清视频| 91国产中文字幕| 日本爱情动作片www.在线观看| 国产高清三级在线| 高清黄色对白视频在线免费看| 激情五月婷婷亚洲| 欧美亚洲日本最大视频资源| 寂寞人妻少妇视频99o| 中文字幕最新亚洲高清| 飞空精品影院首页| 日本免费在线观看一区| 亚洲av男天堂| 97人妻天天添夜夜摸| 少妇的逼水好多| 蜜臀久久99精品久久宅男| 搡老乐熟女国产| 国产乱人偷精品视频| 一区二区日韩欧美中文字幕 | 亚洲综合色惰| 国产成人aa在线观看| 黄色一级大片看看| 啦啦啦在线观看免费高清www| 天堂俺去俺来也www色官网| 久久久久久久久久成人| 观看av在线不卡| 一级毛片我不卡| 建设人人有责人人尽责人人享有的| 另类亚洲欧美激情| 在线天堂最新版资源| 有码 亚洲区| 国产精品国产三级国产专区5o| 国产探花极品一区二区| www.色视频.com| 五月天丁香电影| 蜜桃在线观看..| 亚洲天堂av无毛| 18禁在线无遮挡免费观看视频| 国产免费又黄又爽又色| 亚洲av国产av综合av卡| 啦啦啦视频在线资源免费观看| 涩涩av久久男人的天堂| 精品人妻偷拍中文字幕| 校园人妻丝袜中文字幕| 免费黄频网站在线观看国产| av网站免费在线观看视频| 国产激情久久老熟女| 国产精品嫩草影院av在线观看| 亚洲国产毛片av蜜桃av| 精品熟女少妇av免费看| 国产片内射在线| 亚洲人与动物交配视频| 在线精品无人区一区二区三| 久久这里只有精品19| 黄色视频在线播放观看不卡| 在线观看人妻少妇| 免费高清在线观看视频在线观看| 国精品久久久久久国模美| 亚洲精品成人av观看孕妇| 亚洲国产精品999| 18在线观看网站| 精品少妇久久久久久888优播| 久热久热在线精品观看| 夫妻性生交免费视频一级片| 五月玫瑰六月丁香| 欧美成人午夜精品| 波野结衣二区三区在线| 高清av免费在线| 999精品在线视频| 国产有黄有色有爽视频| 国国产精品蜜臀av免费| 久久久久久久国产电影| 亚洲一码二码三码区别大吗| 久久久久人妻精品一区果冻| 内地一区二区视频在线| 成人毛片60女人毛片免费| av一本久久久久| 精品视频人人做人人爽| xxx大片免费视频| 老女人水多毛片| 亚洲综合色惰| 乱码一卡2卡4卡精品| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 免费大片黄手机在线观看| av又黄又爽大尺度在线免费看| 免费黄色在线免费观看| 色网站视频免费| www.熟女人妻精品国产 | 天天操日日干夜夜撸| 欧美日韩视频高清一区二区三区二| 亚洲av在线观看美女高潮| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 80岁老熟妇乱子伦牲交| av免费在线看不卡| 亚洲精品456在线播放app| 国产有黄有色有爽视频| av国产久精品久网站免费入址| 一个人免费看片子| 哪个播放器可以免费观看大片| 99热6这里只有精品| 这个男人来自地球电影免费观看 | 国产免费视频播放在线视频| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 欧美人与善性xxx| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 亚洲精品美女久久av网站| 久久久久久久精品精品| 看十八女毛片水多多多| 亚洲综合精品二区| 成人毛片60女人毛片免费| av女优亚洲男人天堂| 国产免费福利视频在线观看| 水蜜桃什么品种好| freevideosex欧美| 黄色视频在线播放观看不卡| 亚洲成人av在线免费| 国产精品一区二区在线不卡| 国产色爽女视频免费观看| 母亲3免费完整高清在线观看 | 美女国产视频在线观看| 久久久久国产精品人妻一区二区| videosex国产| av一本久久久久| 一个人免费看片子| 99香蕉大伊视频| 久久久久国产精品人妻一区二区| 精品国产一区二区三区久久久樱花| 国产精品女同一区二区软件| 欧美精品高潮呻吟av久久| 成人手机av| 亚洲国产精品成人久久小说| 国产一区二区激情短视频 | 日本与韩国留学比较| 男人爽女人下面视频在线观看| 国产精品.久久久| 人人澡人人妻人| 97超碰精品成人国产| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 99热国产这里只有精品6| av线在线观看网站| 五月玫瑰六月丁香| 久久久欧美国产精品| 免费大片18禁| 色94色欧美一区二区| 成人无遮挡网站| 亚洲欧美精品自产自拍| 九色亚洲精品在线播放| 大香蕉久久成人网| 天天躁夜夜躁狠狠躁躁| 亚洲成国产人片在线观看| 日韩电影二区| 9色porny在线观看| 成人免费观看视频高清| 高清欧美精品videossex| 国产女主播在线喷水免费视频网站| 97精品久久久久久久久久精品| 亚洲精品,欧美精品| 欧美成人精品欧美一级黄| 国产麻豆69| 在线观看一区二区三区激情| 久久久精品94久久精品| 欧美激情国产日韩精品一区| 精品卡一卡二卡四卡免费| 精品国产露脸久久av麻豆| 精品一区二区三区视频在线| 在线观看免费日韩欧美大片| 中文字幕制服av| 亚洲av.av天堂| 亚洲第一av免费看| 成人18禁高潮啪啪吃奶动态图| 久久久久国产网址| 最近手机中文字幕大全| 亚洲av欧美aⅴ国产| 青春草国产在线视频| 精品国产乱码久久久久久小说| 国产欧美日韩综合在线一区二区| 欧美日韩国产mv在线观看视频| 寂寞人妻少妇视频99o| 侵犯人妻中文字幕一二三四区| 久久99精品国语久久久| 午夜精品国产一区二区电影| 少妇的丰满在线观看| 秋霞伦理黄片| 看免费av毛片| 人体艺术视频欧美日本| 曰老女人黄片| 男人操女人黄网站| 久久久国产欧美日韩av| 欧美另类一区| 精品酒店卫生间| 性色avwww在线观看| 美女大奶头黄色视频| 久久99热6这里只有精品| 国产国语露脸激情在线看| 国产精品久久久av美女十八| 国产高清国产精品国产三级| 国产白丝娇喘喷水9色精品| 日韩免费高清中文字幕av| 亚洲精品aⅴ在线观看| 啦啦啦啦在线视频资源| 久久久久久伊人网av| 日韩精品免费视频一区二区三区 | 国产精品久久久av美女十八| 中国国产av一级| 久久精品人人爽人人爽视色| 亚洲欧美中文字幕日韩二区| 国产熟女午夜一区二区三区| 久久久久久久亚洲中文字幕| 两个人看的免费小视频| 久久女婷五月综合色啪小说| 九色亚洲精品在线播放| 26uuu在线亚洲综合色| 日本欧美视频一区| 视频中文字幕在线观看| 在线天堂最新版资源| 亚洲精品第二区| 免费播放大片免费观看视频在线观看| 午夜影院在线不卡| 久久综合国产亚洲精品| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 亚洲欧洲国产日韩| 男人添女人高潮全过程视频| 黄片播放在线免费| 国产av精品麻豆| 少妇 在线观看| 视频区图区小说| 国产日韩欧美在线精品| 色视频在线一区二区三区| 精品国产一区二区三区久久久樱花| 熟女人妻精品中文字幕| 久久久久精品久久久久真实原创| 免费少妇av软件| 观看美女的网站| 中文字幕人妻丝袜制服| 国产69精品久久久久777片| 亚洲av男天堂| 欧美激情 高清一区二区三区| 99九九在线精品视频| 久久免费观看电影| 亚洲精品国产av蜜桃| 国产国拍精品亚洲av在线观看| 免费看光身美女| 在现免费观看毛片| 国产熟女午夜一区二区三区| 日韩视频在线欧美| xxxhd国产人妻xxx| 国产一区二区三区综合在线观看 | 亚洲经典国产精华液单| 久久国产亚洲av麻豆专区| kizo精华| 美女国产视频在线观看| 丰满饥渴人妻一区二区三| 永久免费av网站大全| 亚洲国产精品成人久久小说| 免费女性裸体啪啪无遮挡网站| 精品少妇内射三级| av.在线天堂| 精品一区二区三区视频在线| 国产极品天堂在线| 免费看av在线观看网站| 这个男人来自地球电影免费观看 | 91久久精品国产一区二区三区| 深夜精品福利| 人人澡人人妻人| 日本av手机在线免费观看| 久久亚洲国产成人精品v| 欧美性感艳星| 看十八女毛片水多多多| 欧美人与性动交α欧美软件 | 国产亚洲精品久久久com| 日韩人妻精品一区2区三区| 秋霞在线观看毛片| 美女主播在线视频| 久久精品国产鲁丝片午夜精品| 久久久久视频综合| 久久婷婷青草| 久久热在线av| 久久久久久久久久久免费av| 久久精品国产综合久久久 | 啦啦啦视频在线资源免费观看| 精品国产一区二区三区久久久樱花| 99久国产av精品国产电影| 久久免费观看电影| 丝袜在线中文字幕| 国产一区二区在线观看日韩| 高清在线视频一区二区三区| 久久久国产一区二区| 国产成人午夜福利电影在线观看| 在线观看美女被高潮喷水网站| 十八禁网站网址无遮挡| 啦啦啦中文免费视频观看日本| 亚洲成色77777| av不卡在线播放| 边亲边吃奶的免费视频| 老司机亚洲免费影院| 国产免费一级a男人的天堂| 少妇人妻精品综合一区二区| 汤姆久久久久久久影院中文字幕| 欧美3d第一页| 菩萨蛮人人尽说江南好唐韦庄| 国产xxxxx性猛交| www.熟女人妻精品国产 | 亚洲成国产人片在线观看| 美女中出高潮动态图| 国产精品女同一区二区软件| 久久久国产精品麻豆| 欧美人与性动交α欧美软件 | 免费黄色在线免费观看| 国产一区二区在线观看日韩| 涩涩av久久男人的天堂| 免费久久久久久久精品成人欧美视频 | 王馨瑶露胸无遮挡在线观看| 精品酒店卫生间| 国内精品宾馆在线| 日韩av在线免费看完整版不卡| 美女国产高潮福利片在线看| 日韩av不卡免费在线播放| 日本欧美国产在线视频| 22中文网久久字幕| 亚洲成人一二三区av| 日本午夜av视频| 久久久a久久爽久久v久久| 2018国产大陆天天弄谢| 中国美白少妇内射xxxbb| 亚洲精品av麻豆狂野| 亚洲人成77777在线视频|