• <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
    av在线观看视频网站免费| 国产精品一区二区三区四区久久| 18禁在线播放成人免费| 一级二级三级毛片免费看| 亚洲婷婷狠狠爱综合网| 免费黄网站久久成人精品| 久久亚洲国产成人精品v| 最近中文字幕2019免费版| 91aial.com中文字幕在线观看| 97在线视频观看| 在现免费观看毛片| 久久久精品大字幕| 日本-黄色视频高清免费观看| 少妇丰满av| 丰满乱子伦码专区| 国产高清国产精品国产三级 | 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 少妇被粗大猛烈的视频| 免费观看人在逋| 国产麻豆成人av免费视频| 一级爰片在线观看| 1000部很黄的大片| 麻豆国产97在线/欧美| 亚洲精品国产av成人精品| 精品99又大又爽又粗少妇毛片| 国产精品99久久久久久久久| 成人三级黄色视频| 亚洲欧美日韩卡通动漫| 欧美高清性xxxxhd video| 搞女人的毛片| www日本黄色视频网| 九九热线精品视视频播放| 亚洲欧美清纯卡通| 久久久精品大字幕| 欧美一区二区精品小视频在线| av又黄又爽大尺度在线免费看 | 两性午夜刺激爽爽歪歪视频在线观看| 日日啪夜夜撸| 免费看美女性在线毛片视频| 国产精品蜜桃在线观看| 少妇熟女aⅴ在线视频| 亚洲真实伦在线观看| 欧美一区二区国产精品久久精品| 久久人人爽人人片av| 麻豆一二三区av精品| 国产免费福利视频在线观看| 国产精品三级大全| 久久人人爽人人片av| 亚洲人成网站在线观看播放| 神马国产精品三级电影在线观看| 超碰av人人做人人爽久久| 久久久久性生活片| 伦精品一区二区三区| 免费看a级黄色片| 久久久精品94久久精品| a级毛色黄片| 18+在线观看网站| 男女啪啪激烈高潮av片| 国产精品精品国产色婷婷| 麻豆乱淫一区二区| 欧美日韩国产亚洲二区| 国产午夜精品久久久久久一区二区三区| 日本猛色少妇xxxxx猛交久久| 久久久国产成人免费| 夜夜爽夜夜爽视频| ponron亚洲| 男女边吃奶边做爰视频| a级毛色黄片| 久久亚洲国产成人精品v| 只有这里有精品99| 一区二区三区四区激情视频| 尤物成人国产欧美一区二区三区| 国产午夜福利久久久久久| 亚洲av免费在线观看| 春色校园在线视频观看| 一级爰片在线观看| 欧美一区二区精品小视频在线| 亚洲国产欧美人成| 美女高潮的动态| 色5月婷婷丁香| 中文欧美无线码| 久久久久久大精品| 国产v大片淫在线免费观看| 内射极品少妇av片p| 国产精品无大码| 搡老妇女老女人老熟妇| 精品久久久久久久久久久久久| 黄色日韩在线| 中文精品一卡2卡3卡4更新| 禁无遮挡网站| 天堂av国产一区二区熟女人妻| 国内少妇人妻偷人精品xxx网站| 九色成人免费人妻av| 国产人妻一区二区三区在| 午夜视频国产福利| 五月伊人婷婷丁香| 性色avwww在线观看| 国产高清三级在线| 免费看a级黄色片| 99九九线精品视频在线观看视频| 午夜精品在线福利| 国产欧美日韩精品一区二区| 永久网站在线| 久久这里只有精品中国| 最后的刺客免费高清国语| 又黄又爽又刺激的免费视频.| 午夜福利在线观看免费完整高清在| 国产亚洲91精品色在线| 1024手机看黄色片| 欧美+日韩+精品| 黄片wwwwww| 精品久久久久久久末码| 成人美女网站在线观看视频| 日产精品乱码卡一卡2卡三| 久久综合国产亚洲精品| 国产av一区在线观看免费| 亚洲中文字幕一区二区三区有码在线看| 天堂√8在线中文| 51国产日韩欧美| 在线播放无遮挡| 国产精品无大码| 一个人观看的视频www高清免费观看| 高清在线视频一区二区三区 | 深爱激情五月婷婷| 嘟嘟电影网在线观看| 国产精品一区二区三区四区免费观看| 国产精品人妻久久久影院| 在线免费十八禁| 99久久精品一区二区三区| 欧美又色又爽又黄视频| 18禁在线无遮挡免费观看视频| 丝袜喷水一区| 亚洲av中文av极速乱| 国产成人免费观看mmmm| 狠狠狠狠99中文字幕| 国产精品永久免费网站| 日韩,欧美,国产一区二区三区 | 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 永久网站在线| 国产黄色小视频在线观看| 精品99又大又爽又粗少妇毛片| 亚洲中文字幕日韩| 男女边吃奶边做爰视频| 色尼玛亚洲综合影院| 国产免费福利视频在线观看| 嫩草影院精品99| 极品教师在线视频| 夫妻性生交免费视频一级片| 男女国产视频网站| 国产av一区在线观看免费| 午夜福利网站1000一区二区三区| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 精品熟女少妇av免费看| 国产亚洲最大av| 中文字幕制服av| 日本免费在线观看一区| 一级毛片久久久久久久久女| 国产精品麻豆人妻色哟哟久久 | 亚洲精品国产成人久久av| 啦啦啦啦在线视频资源| 小说图片视频综合网站| 国产男人的电影天堂91| 秋霞在线观看毛片| 亚洲国产高清在线一区二区三| 精品久久国产蜜桃| 免费黄网站久久成人精品| 国内精品美女久久久久久| 国产视频首页在线观看| 级片在线观看| 国产一区二区三区av在线| 亚洲精品影视一区二区三区av| 国产免费又黄又爽又色| 汤姆久久久久久久影院中文字幕 | 老司机影院成人| 欧美高清成人免费视频www| 日本午夜av视频| 禁无遮挡网站| 日韩成人伦理影院| 国产成年人精品一区二区| 少妇人妻精品综合一区二区| 精品国产一区二区三区久久久樱花 | 国产精品人妻久久久影院| 久久久久性生活片| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 欧美成人午夜免费资源| 蜜臀久久99精品久久宅男| 99久久无色码亚洲精品果冻| 尾随美女入室| 国产亚洲91精品色在线| 直男gayav资源| 精品人妻偷拍中文字幕| av在线播放精品| 久久久a久久爽久久v久久| 国产亚洲5aaaaa淫片| 一边摸一边抽搐一进一小说| 久久精品国产亚洲av涩爱| 亚洲精品色激情综合| 又爽又黄无遮挡网站| 中文亚洲av片在线观看爽| 舔av片在线| 亚洲图色成人| 久久欧美精品欧美久久欧美| 亚洲四区av| 国产精品99久久久久久久久| 变态另类丝袜制服| 国产精品久久电影中文字幕| 国产三级在线视频| 午夜爱爱视频在线播放| 少妇高潮的动态图| 狂野欧美激情性xxxx在线观看| 亚洲精品乱码久久久v下载方式| 成人毛片60女人毛片免费| 日本免费a在线| 最近最新中文字幕大全电影3| 人人妻人人澡欧美一区二区| 久久久久网色| 在线观看66精品国产| 丰满乱子伦码专区| 亚洲在线观看片| 亚洲欧美日韩无卡精品| 在线免费十八禁| 色综合色国产| 亚洲av电影在线观看一区二区三区 | 亚洲欧美日韩卡通动漫| 久久亚洲国产成人精品v| 日韩大片免费观看网站 | 国内少妇人妻偷人精品xxx网站| 91久久精品国产一区二区成人| 中文在线观看免费www的网站| 亚洲精品日韩在线中文字幕| 国产片特级美女逼逼视频| 国产毛片a区久久久久| 亚洲内射少妇av| 久久久久久久国产电影| 久久99精品国语久久久| 久久久国产成人免费| 成年版毛片免费区| 少妇裸体淫交视频免费看高清| 国产亚洲av片在线观看秒播厂 | 日韩制服骚丝袜av| 久久精品国产99精品国产亚洲性色| 亚洲国产高清在线一区二区三| 天美传媒精品一区二区| 在线播放无遮挡| 国产欧美另类精品又又久久亚洲欧美| 如何舔出高潮| 在线观看av片永久免费下载| 男人的好看免费观看在线视频| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 久久婷婷人人爽人人干人人爱| 禁无遮挡网站| 亚洲成人av在线免费| 色综合站精品国产| 国产成人a区在线观看| av国产免费在线观看| 99在线视频只有这里精品首页| 岛国毛片在线播放| 成人漫画全彩无遮挡| 日韩视频在线欧美| 久久亚洲精品不卡| 久久精品国产99精品国产亚洲性色| 国产高潮美女av| 岛国在线免费视频观看| 国产免费男女视频| 国产精品精品国产色婷婷| 欧美一级a爱片免费观看看| 日韩中字成人| 国产精品伦人一区二区| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 国产av码专区亚洲av| 我要搜黄色片| 尾随美女入室| 最近视频中文字幕2019在线8| 亚洲国产精品国产精品| 日韩精品有码人妻一区| 国产精品伦人一区二区| 黄片无遮挡物在线观看| 能在线免费观看的黄片| 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 精品久久久久久久久亚洲| 男女视频在线观看网站免费| 天天一区二区日本电影三级| 免费大片18禁| 天堂av国产一区二区熟女人妻| 亚洲一级一片aⅴ在线观看| 欧美极品一区二区三区四区| 91aial.com中文字幕在线观看| 能在线免费观看的黄片| 纵有疾风起免费观看全集完整版 | 国产精品美女特级片免费视频播放器| 国产一级毛片七仙女欲春2| 亚洲精品影视一区二区三区av| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 91狼人影院| АⅤ资源中文在线天堂| 国产乱人偷精品视频| 精品久久久久久久久亚洲| 桃色一区二区三区在线观看| 久久99精品国语久久久| a级毛色黄片| 中文天堂在线官网| 久久久色成人| 国产免费福利视频在线观看| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 成人特级av手机在线观看| 久久久精品大字幕| 国产亚洲91精品色在线| 精品国内亚洲2022精品成人| 欧美精品国产亚洲| 精品久久久久久久久久久久久| 国产一区二区在线观看日韩| 午夜福利成人在线免费观看| 少妇被粗大猛烈的视频| 晚上一个人看的免费电影| 国产极品天堂在线| 日韩一区二区三区影片| 国产av在哪里看| 热99re8久久精品国产| 日韩欧美 国产精品| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人中文| 亚洲欧美日韩卡通动漫| 国产精品.久久久| 精品免费久久久久久久清纯| 日本熟妇午夜| 一区二区三区四区激情视频| 少妇高潮的动态图| 毛片一级片免费看久久久久| 久久久久久久久久久免费av| 亚洲高清免费不卡视频| 夫妻性生交免费视频一级片| 最后的刺客免费高清国语| 久久久久久久久久久免费av| 18禁在线无遮挡免费观看视频| 伦理电影大哥的女人| 一级毛片久久久久久久久女| 不卡视频在线观看欧美| 狠狠狠狠99中文字幕| 日韩国内少妇激情av| 免费观看人在逋| 欧美变态另类bdsm刘玥| 亚洲伊人久久精品综合 | 亚洲av二区三区四区| 亚洲怡红院男人天堂| 久久国产乱子免费精品| 天天躁日日操中文字幕| 久久久a久久爽久久v久久| 久久久久久久久久久免费av| 97在线视频观看| 超碰97精品在线观看| 日韩一本色道免费dvd| 亚洲国产精品成人综合色| 久久久精品大字幕| 热99在线观看视频| 青春草国产在线视频| 少妇高潮的动态图| 亚洲av电影在线观看一区二区三区 | 午夜免费男女啪啪视频观看| 在线观看av片永久免费下载| 国产高清有码在线观看视频| 国产色婷婷99| 美女大奶头视频| 2021少妇久久久久久久久久久| 亚洲欧美日韩东京热| 中文精品一卡2卡3卡4更新| 久热久热在线精品观看| 精品一区二区三区人妻视频| 最近中文字幕高清免费大全6| 久久这里有精品视频免费| 精品一区二区免费观看| 少妇熟女欧美另类| 日本-黄色视频高清免费观看| 亚洲av熟女| 小蜜桃在线观看免费完整版高清| 国产精品,欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲5aaaaa淫片| 国产真实乱freesex| 亚洲成人av在线免费| 我要搜黄色片| 久久精品熟女亚洲av麻豆精品 | av在线观看视频网站免费| 在线免费十八禁| .国产精品久久| 国产亚洲最大av| 国产在视频线在精品| 日韩一本色道免费dvd| 亚洲精品国产av成人精品| 欧美不卡视频在线免费观看| 激情 狠狠 欧美| 直男gayav资源| 精品不卡国产一区二区三区| 国产精品日韩av在线免费观看| 欧美成人精品欧美一级黄| 成年女人看的毛片在线观看| 亚洲av电影在线观看一区二区三区 | 三级经典国产精品| 国产av不卡久久| 亚洲精品456在线播放app| 亚洲电影在线观看av| 国产精品久久久久久av不卡| a级一级毛片免费在线观看| 国内揄拍国产精品人妻在线| 日韩高清综合在线| 久久这里有精品视频免费| 婷婷六月久久综合丁香| 欧美丝袜亚洲另类| 午夜福利网站1000一区二区三区| 亚洲精品乱久久久久久| 麻豆成人午夜福利视频| 插逼视频在线观看| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 久久久欧美国产精品| 男女边吃奶边做爰视频| 免费观看精品视频网站| 日韩制服骚丝袜av| 国产v大片淫在线免费观看| 免费看a级黄色片| 乱系列少妇在线播放| 国产午夜福利久久久久久| 久久午夜福利片| 人妻制服诱惑在线中文字幕| 一级二级三级毛片免费看| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| 国产一区亚洲一区在线观看| 免费看美女性在线毛片视频| 一级毛片aaaaaa免费看小| 国产私拍福利视频在线观看| 极品教师在线视频| 久久99蜜桃精品久久| 国产精品嫩草影院av在线观看| 一个人看的www免费观看视频| 天美传媒精品一区二区| 日韩在线高清观看一区二区三区| 精品久久久久久久久av| 国产一区有黄有色的免费视频 | 免费观看在线日韩| 一级毛片我不卡| 亚洲经典国产精华液单| 一二三四中文在线观看免费高清| 在线播放无遮挡| 一级毛片电影观看 | 久久久午夜欧美精品| 精品一区二区三区视频在线| 精品国产露脸久久av麻豆 | 国产国拍精品亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产午夜精品久久久久久一区二区三区| 午夜激情福利司机影院| 国产免费福利视频在线观看| 禁无遮挡网站| 在线天堂最新版资源| 一级毛片我不卡| 91久久精品电影网| 日本黄色片子视频| 免费av毛片视频| 大又大粗又爽又黄少妇毛片口| 噜噜噜噜噜久久久久久91| 黄色配什么色好看| 午夜精品在线福利| 2021天堂中文幕一二区在线观| 韩国高清视频一区二区三区| 久久这里只有精品中国| 欧美潮喷喷水| 国产黄色小视频在线观看| av福利片在线观看| 日韩 亚洲 欧美在线| av免费观看日本| 老女人水多毛片| 黑人高潮一二区| 秋霞伦理黄片| 变态另类丝袜制服| 久久鲁丝午夜福利片| 性色avwww在线观看| 国产乱人视频| 亚洲av电影不卡..在线观看| 国产高清视频在线观看网站| 色5月婷婷丁香| 国产成人a∨麻豆精品| 日本黄大片高清| 国产视频内射| 亚洲国产精品专区欧美| 国产精品野战在线观看| 国产黄片视频在线免费观看| 国产午夜福利久久久久久| 中文字幕精品亚洲无线码一区| 男人舔女人下体高潮全视频| 国产淫语在线视频| av在线播放精品| 日本黄大片高清| 国产69精品久久久久777片| 少妇丰满av| 联通29元200g的流量卡| 精品免费久久久久久久清纯| 亚洲中文字幕一区二区三区有码在线看| 国产视频内射| av免费在线看不卡| 午夜免费男女啪啪视频观看| 18禁在线无遮挡免费观看视频| 精品熟女少妇av免费看| 麻豆av噜噜一区二区三区| 久久久久久久久大av| 少妇的逼水好多| 国产精品久久久久久av不卡| 日韩一本色道免费dvd| 精品99又大又爽又粗少妇毛片| 久久精品影院6| 日本熟妇午夜| 久久久国产成人免费| 精华霜和精华液先用哪个| 国产色爽女视频免费观看| 欧美性猛交╳xxx乱大交人| 国产老妇女一区| 精品无人区乱码1区二区| 九九久久精品国产亚洲av麻豆| 99久久精品一区二区三区| 色播亚洲综合网| 永久网站在线| 国产v大片淫在线免费观看| 国产男人的电影天堂91| 搞女人的毛片| 1024手机看黄色片| 3wmmmm亚洲av在线观看| 精品国内亚洲2022精品成人| 日日摸夜夜添夜夜爱| 国产精品综合久久久久久久免费| 成年av动漫网址| 精品久久国产蜜桃| 国产综合懂色| 婷婷六月久久综合丁香| 欧美高清成人免费视频www| 丰满少妇做爰视频| 国产三级在线视频| 99热这里只有是精品50| 久久99热这里只频精品6学生 | 亚洲高清免费不卡视频| 波多野结衣高清无吗| av黄色大香蕉| 偷拍熟女少妇极品色| 日本-黄色视频高清免费观看| 久久精品国产99精品国产亚洲性色| 精品午夜福利在线看| 久久精品国产亚洲av涩爱| 国产老妇伦熟女老妇高清| 3wmmmm亚洲av在线观看| 免费一级毛片在线播放高清视频| 精品久久久久久久末码| 只有这里有精品99| 久久人妻av系列| 观看美女的网站| 久久久久久久久久久免费av| 精品人妻熟女av久视频| 国产av在哪里看| 亚洲不卡免费看| 男女国产视频网站| 九九在线视频观看精品| 免费看日本二区| 麻豆国产97在线/欧美| av福利片在线观看| 国产69精品久久久久777片| 亚洲国产欧洲综合997久久,| 女人久久www免费人成看片 | 九九久久精品国产亚洲av麻豆| 亚洲成av人片在线播放无| 久久国内精品自在自线图片| 亚洲国产精品sss在线观看| 免费人成在线观看视频色| 亚洲五月天丁香| 日本五十路高清| 最近最新中文字幕大全电影3| 嫩草影院新地址| 大香蕉久久网| av在线播放精品| 国产熟女欧美一区二区| 中文乱码字字幕精品一区二区三区 | 欧美一区二区精品小视频在线| 亚洲人与动物交配视频| 3wmmmm亚洲av在线观看| 天堂av国产一区二区熟女人妻| 欧美xxxx黑人xx丫x性爽| 日本与韩国留学比较| 国产淫语在线视频| 欧美成人一区二区免费高清观看| 人人妻人人澡欧美一区二区| 欧美97在线视频| 欧美成人a在线观看| 国产精品福利在线免费观看| 国产成人精品一,二区| 一个人观看的视频www高清免费观看| 午夜福利在线在线| 国产免费又黄又爽又色| 性插视频无遮挡在线免费观看| 汤姆久久久久久久影院中文字幕 | 国产亚洲一区二区精品| 国产精品一区二区三区四区久久| h日本视频在线播放| 欧美3d第一页| 最近最新中文字幕免费大全7| 亚洲国产精品合色在线| 免费在线观看成人毛片| 国产精品人妻久久久影院| 亚洲av二区三区四区| 久久精品国产99精品国产亚洲性色|