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

    紫花苜蓿對黃土邊坡淺層破壞防護(hù)時間效應(yīng)的數(shù)值分析

    2022-11-13 08:03:26毛正君張瑾鴿畢銀麗孫偉博
    關(guān)鍵詞:植物

    毛正君,張瑾鴿,畢銀麗,孫偉博,安 寧

    紫花苜蓿對黃土邊坡淺層破壞防護(hù)時間效應(yīng)的數(shù)值分析

    毛正君1,2,3,張瑾鴿1,畢銀麗1,孫偉博4,安 寧1

    (1. 西安科技大學(xué)地質(zhì)與環(huán)境學(xué)院,西安 710054;2. 西安科技大學(xué)煤炭綠色開采地質(zhì)研究院,西安 710054;3. 陜西省煤炭綠色開發(fā)地質(zhì)保障重點(diǎn)實(shí)驗(yàn)室,西安 710054;4. 西安科技大學(xué)能源學(xué)院,西安 710054)

    為探究紫花苜蓿對黃土邊坡淺層破壞防護(hù)的時間效應(yīng),該研究考慮降雨入滲條件并以含不同生長期紫花苜蓿的黃土邊坡為例,基于室內(nèi)土工試驗(yàn)獲取參數(shù),分別采用含植物根系的無限邊坡模型和數(shù)值模擬兩種方法進(jìn)行穩(wěn)定性分析。結(jié)果表明:降雨作用使得邊坡潛在滑動面位置由坡體內(nèi)部轉(zhuǎn)移到降雨最大入滲深度處;隨著坡面植被的生長,邊坡淺層土體的強(qiáng)度逐漸提高,潛在滑動面又逐漸轉(zhuǎn)移到邊坡內(nèi)部。根據(jù)數(shù)值模擬強(qiáng)度折減法,邊坡穩(wěn)定系數(shù)隨著植被生長不斷增大。對比無限邊坡模型與數(shù)值模擬方法計(jì)算的邊坡穩(wěn)定系數(shù),植物根系生長到1.0 m且不論在天然狀態(tài)或降雨條件下黃土邊坡均處于穩(wěn)定狀態(tài),植物根系在0~0.8 m且在降雨條件下時處于不穩(wěn)定狀態(tài),隨著邊坡植物生長期的延長,邊坡穩(wěn)定系數(shù)總體變化趨勢一致,即先降低再逐漸增大。草本植物隨著生長期的延長,對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用,從而草本植物根系生長能夠提高黃土邊坡穩(wěn)定性。當(dāng)草本植物生長時間達(dá)到150 d時,能夠有效防治黃土坡面侵蝕,并提高黃土邊坡穩(wěn)定性,因此草本植物生長的前5個月為關(guān)鍵期。研究成果對于黃土邊坡淺層破壞防護(hù)以及實(shí)現(xiàn)黃土地區(qū)可持續(xù)發(fā)展具有重要的理論意義和應(yīng)用價值。

    邊坡穩(wěn)定性;數(shù)值模擬;黃土;紫花苜蓿;邊坡淺層破壞;時間效應(yīng);無限邊坡模型

    0 引 言

    黃土是一種碎散的多孔多相顆粒介質(zhì),由于其獨(dú)特的物質(zhì)來源、特有的生成環(huán)境及應(yīng)力-應(yīng)變歷史而具有特殊的微結(jié)構(gòu)[1]。當(dāng)黃土遇水浸濕或遭受強(qiáng)震作用時,其結(jié)構(gòu)會出現(xiàn)崩潰性破壞,導(dǎo)致強(qiáng)度突發(fā)性完全喪失,從而誘發(fā)嚴(yán)重的地質(zhì)災(zāi)害或者工程病害[2]。黃土邊坡在降雨作用下往往先出現(xiàn)邊坡坡面侵蝕,隨著降雨次數(shù)的增加將會導(dǎo)致邊坡局部破壞,甚至影響邊坡的整體穩(wěn)定性[3]。黃土梁峁地形切割強(qiáng)烈,不利于降水入滲及地下水賦存,因此降雨直接入滲的影響深度是有限的[4];表層黃土軟化、質(zhì)量增加,造成數(shù)量多、規(guī)模小、多發(fā)于30°~50°黃土斜坡上的邊坡淺層破壞[5]。中國黃土高原作為世界上黃土沉積厚度最大地區(qū)之一,黃土邊坡淺層破壞數(shù)目每年高達(dá)數(shù)百個[6-7]。黃土邊坡淺層破壞一般發(fā)生在坡體表層或坡面下2 m的范圍內(nèi),其破壞類型有坡面侵蝕、風(fēng)化剝落、溜坍、錯落、淺層滑坡等[8],可造成黃土地區(qū)土地可持續(xù)生產(chǎn)力的損失和破壞,是當(dāng)前黃土地區(qū)研究的重點(diǎn)。

    自20世紀(jì)60年代末以來,計(jì)算土力學(xué)一直是熱點(diǎn)研究領(lǐng)域[9]。近年來各種數(shù)值模擬計(jì)算技術(shù)的迅猛發(fā)展,為土體破壞現(xiàn)象的模擬提供了有力的手段[10]。通過合理的對根系加固作用進(jìn)行數(shù)值模擬,進(jìn)而對根系加固土體的機(jī)理進(jìn)行分析,能為植物防護(hù)技術(shù)的應(yīng)用提供強(qiáng)有力的理論依據(jù)[11]。周群華等[12]采用有限元數(shù)值模擬方法,對植物根系錨固邊坡的作用及其力學(xué)機(jī)理進(jìn)行分析,發(fā)現(xiàn)根系固土是因?yàn)楦翟黾恿送馏w的抗剪強(qiáng)度,從而增強(qiáng)了邊坡土體的抗滑力。Temgoua等[13]采用SIMULIA軟件建立三維數(shù)值模型,發(fā)現(xiàn)整個邊坡的穩(wěn)定性主要取決于根系的深度和根系提供的附加黏聚力。Li等[14]采用數(shù)值模擬的方法,將植物根系結(jié)構(gòu)概化,分析了植物根系空間布局對邊坡穩(wěn)定性的影響。黃建坤等[15]基于有限元軟件,模擬了含草邊坡和無草邊坡的應(yīng)力和應(yīng)變規(guī)律,通過邊坡安全系數(shù)反映了根系的固土效應(yīng)。由于高密度的細(xì)根能夠和土體顆粒緊密地纏繞在一起,草本植物能有效地提高根土復(fù)合體的穩(wěn)定性[16],是防治黃土邊坡淺層破壞土木工程解決方案的一種生態(tài)替代方案[17]。紫花苜蓿屬多年生草本植物,具有繁殖能力強(qiáng)、適應(yīng)性強(qiáng)、經(jīng)濟(jì)效益高等特點(diǎn)[18],能夠有效提高土體的抗沖刷能力,并且根系的固土能力強(qiáng),具有顯著的水土保持效果,是中國黃土高原地區(qū)水土保持的先鋒植物之一[19]。

    為了探究紫花苜蓿對黃土邊坡淺層破壞防護(hù)的時間效應(yīng),本文采用套管人工種植紫花苜蓿的方式,獲取含紫花苜蓿根系的原狀土樣進(jìn)行試驗(yàn),并分別利用數(shù)值模擬和無限邊坡模型,對不同生長期紫花苜蓿天然狀態(tài)和降雨條件下黃土邊坡的不利情況進(jìn)行穩(wěn)定性分析,探究不同生長期紫花苜蓿對黃土邊坡淺層破壞防護(hù)作用。研究成果對于黃土邊坡淺層破壞防護(hù)以及實(shí)現(xiàn)黃土地區(qū)可持續(xù)發(fā)展具有重要的理論意義和應(yīng)用價值。

    1 材料與方法

    1.1 試驗(yàn)材料與試樣制備

    試驗(yàn)用土取自陜西省西安市臨潼區(qū),地理坐標(biāo)為109°20′25″E,34°25′46″N。試驗(yàn)用土為離石黃土(Q2黃土),顏色為褐黃色,天然密度為1.67g/cm3,天然含水率為16.3%,塑限為17.4%,塑性指數(shù)為19.8%,液限為37.2%,液限指數(shù)為?0.06,不均勻系數(shù)(C)和曲率系數(shù)(C)分別為8.36和0.78,為不均勻級配土。

    試驗(yàn)所用紫花苜蓿種子來自內(nèi)蒙古赤峰市林西縣草原站,采用人工種植培養(yǎng)的方式,獲得含不同生長期紫花苜蓿根系的黃土試樣。方法如下:首先,將試驗(yàn)用土取回后經(jīng)自然風(fēng)干,挑出其中的礫石與垃圾,把土塊砸碎后過2 mm篩,之后置于陰涼、干燥處保存;其次,在種植箱內(nèi)鋪一層5 cm左右的土,將切割好的內(nèi)徑為67.8 mm,外徑為75 mm的PVC種植管插入種植箱內(nèi),并在PVC種植管中裝滿黃土,在PVC種植管上端澆水浸潤土體后靜置一段時間,繼續(xù)填滿黃土并澆水浸潤,直至黃土裝滿PVC種植管澆水后不再下沉;再次,在每根PVC種植管中撒上數(shù)量大致相同的紫花苜蓿種子,生長期為60與90 d的紫花苜蓿(根長約15~20 cm)種植于20 cm長的PVC管中,生長期為120與150 d的紫花苜蓿(根長約30~40 cm)種植于40 cm長的PVC管中;最后,為了保證PVC種植管中黃土的含水率大致相同,在養(yǎng)護(hù)過程中定時給PVC管內(nèi)澆入相同質(zhì)量的水。

    為了滿足三軸試驗(yàn)中直徑為39.1 mm、高為80 mm的標(biāo)準(zhǔn)三軸樣的要求,且便于對比,對于兩組根土復(fù)合體土樣,取樣位置分別選取20 cm長PVC管的中部和40 cm長PVC管上半部分20 cm管的中部。

    1.2 根土復(fù)合體力學(xué)指標(biāo)求解

    本文考慮邊坡表層1 m內(nèi)完全飽和,當(dāng)在強(qiáng)度相對較高且埋深較淺的根系分布范圍內(nèi)已然發(fā)生了破壞,那么在強(qiáng)度較低且埋深較大的根系分布范圍內(nèi)必然更容易發(fā)生破壞。此外,所取三軸樣即根土復(fù)合體土樣被看做一個整體進(jìn)行試驗(yàn)。因此,考慮植物不同生長期為變量,根系分布對于本文試驗(yàn)抗剪指標(biāo)測定的影響忽略。采用中國南京土壤儀器廠生產(chǎn)的ZJ型應(yīng)變控制式直剪儀,進(jìn)行直剪快剪試驗(yàn)測定天然素黃土的抗剪強(qiáng)度指標(biāo)黏聚力c和內(nèi)摩擦角。采用中國南京土壤儀器廠生產(chǎn)的TSZ-1型三軸儀,進(jìn)行不固結(jié)不排水(UU)三軸試驗(yàn)測定天然素黃土的彈性模量和泊松比,固結(jié)不排水(CU)三軸試驗(yàn)測定飽和素黃土和含不同生長期紫花苜蓿根系黃土試樣的抗剪強(qiáng)度指標(biāo)黏聚力c、內(nèi)摩擦角、彈性模量和泊松比。

    根據(jù)直剪試驗(yàn)剪切應(yīng)力與剪切位移關(guān)系曲線,選取峰值強(qiáng)度作出不同垂直壓力下所對應(yīng)的抗剪強(qiáng)度坐標(biāo)點(diǎn),并通過擬合求解抗剪強(qiáng)度指標(biāo)。本文根據(jù)三軸試驗(yàn)應(yīng)力應(yīng)變曲線畫出了應(yīng)力莫爾圓,但各個應(yīng)力圓之間無明顯規(guī)律,難以繪制強(qiáng)度包絡(luò)線,因此選用最小二乘法求解抗剪強(qiáng)度指標(biāo)。

    根據(jù)三軸試驗(yàn)得出的應(yīng)力-應(yīng)變關(guān)系曲線可以求得土的變形模量和泊松比,它們不是常量,而是隨應(yīng)力水平變化的參數(shù)。在常規(guī)三軸試驗(yàn)中,當(dāng)試樣只承受豎向應(yīng)力增量Δ1時,除了產(chǎn)生豎向應(yīng)變Δ1之外,同時還產(chǎn)生側(cè)向應(yīng)變Δ3(側(cè)向膨脹)。側(cè)向應(yīng)變增量Δ3與豎向應(yīng)變增量Δ1的比值可用彈性力學(xué)中的泊松比表示。

    在不允許有側(cè)向變形的情況下,土樣受到軸向壓力增量Δ1將會引起側(cè)向壓力的相應(yīng)增量Δ3,二者比值(Δ3/Δ1)稱為土的側(cè)壓力系數(shù)或靜止土壓力系數(shù)0:

    側(cè)壓力系數(shù)與泊松比的關(guān)系為

    1.3 含植物根系的無限邊坡模型

    現(xiàn)實(shí)中的天然邊坡和人工邊坡都是有限長的,但對于一些具有平行于坡面的淺層滑動面的長坡,滑動面的上部和坡腳處的滑動面彎曲部分、兩側(cè)邊界的約束都可以忽略,可以視為二維的無限邊坡[20]。無限邊坡穩(wěn)定性分析模型常常應(yīng)用于淺層滑坡的穩(wěn)定性分析[21]和坡面穩(wěn)定性分析[22]。已有研究認(rèn)為植物根系擁有防止淺層滑坡和地表侵蝕的潛能[23],因此無限邊坡穩(wěn)定性分析模型適用于含植物根系邊坡穩(wěn)定性分析[24-26]。此外,根土復(fù)合體與素土的內(nèi)摩擦角基本一致,僅黏聚力發(fā)生了變化[27-28]。黃土地區(qū)地下水位埋深較大[29-30],降雨入滲深度較淺[5],本文假定不存在地下水滲流作用,且降水入滲影響范圍內(nèi)的土體完全飽和,邊坡穩(wěn)定系數(shù)可以表示為

    式中F是邊坡穩(wěn)定系數(shù);是邊坡土體的重度,kN/m3;是滑動面深度,m;是坡面與水平地面的夾角,(°);是內(nèi)摩擦角,(°);c是根土復(fù)合體的黏聚力,kPa;h是植物根系分布深度,m;c是素土的黏聚力,kPa。

    紫花苜蓿根系分布深度較淺[31],本文假定含紫花苜蓿根系的邊坡土體厚度為1 m。已有研究[5,32]表明黃土地區(qū)降水入滲深度約為1~3 m,本文僅考慮一次降雨事件造成的影響,因此將降水入滲影響深度假定為1 m。本文僅考慮含紫花苜蓿根系的淺層邊坡土體穩(wěn)定性,由于邊坡穩(wěn)定系數(shù)與滑動面深度成反比,因此邊坡穩(wěn)定系數(shù)最小值對應(yīng)的滑動面深度=1.0 m。將黏聚力、內(nèi)摩擦角、重度以及滑動面深度=1.0 m代入式(1),發(fā)現(xiàn)坡度在49°~51°范圍內(nèi)邊坡穩(wěn)定系數(shù)F最小,因此本文假定坡度=50°。

    1.4 含植物根系的邊坡穩(wěn)定性數(shù)值模擬分析

    與現(xiàn)場試驗(yàn)和模型試驗(yàn)方法相比,數(shù)值模擬需要的人力和財(cái)力更少,研究周期更短,靈活性和便利性更好[33]。數(shù)值模擬是用數(shù)值方法定量地描述一個物理系統(tǒng)的演化過程,允許分析復(fù)雜的因素,如非線性行為,包含復(fù)雜的不連續(xù)模式、動態(tài)載荷、非均勻性、各向異性、流體流動和力學(xué)行為的耦合、材料軟化、時間依賴的行為等,還可以不受限制地設(shè)置任意時間尺度或長度的數(shù)值,可以快速、有效地改變測試條件和組件[34]。強(qiáng)度折減法是指采用降低巖土強(qiáng)度的方法來計(jì)算巖土工程的穩(wěn)定系數(shù),通過不斷降低邊坡巖土體的抗剪強(qiáng)度參數(shù)來逐步逼近極限破壞狀態(tài),利用數(shù)值模擬計(jì)算程序自動計(jì)算得到滑動破壞面和強(qiáng)度折減穩(wěn)定系數(shù);運(yùn)用強(qiáng)度折減法分析巖土體穩(wěn)定性具有不需要假定滑動面的形狀和位置、不需進(jìn)行條分和離散的優(yōu)勢[8]。

    1.4.1 工況設(shè)定

    含不同生長期紫花苜蓿根系邊坡穩(wěn)定性數(shù)值模擬分析,共設(shè)置6種工況,見表1。植被的適宜播種時間是秋季和春季,黃土地區(qū)的雨季一般在7-9月,為了使植被護(hù)坡能夠在同年中起到預(yù)期效果,春季是進(jìn)行播種的最佳時機(jī)。假定春季3月進(jìn)行播種,則到雨季來臨時植物的生長周期可能是150 d,因此本次試驗(yàn)中采樣周期選擇60、90、120和150 d共4個采樣周期。本文所選草本植物紫花苜蓿,屬多年生深根系草本植物,主根入土深,在播種當(dāng)年根系深度可達(dá)2 m左右[35],為考慮野外邊坡紫花苜蓿的實(shí)際生長情況,反映植物根系隨生長期的變化特征,假定植物生長期為60、90、120、150 d時,根系深度分別為0.4、0.6、0.8、1.0 m。

    表1 工況設(shè)定

    1.4.2 模型建立

    本文基于數(shù)值模擬軟件FLAC3D,進(jìn)行含不同生長期紫花苜蓿根系黃土邊坡淺層破壞防護(hù)效應(yīng)分析。根據(jù)基本假設(shè),含紫花苜蓿根系的邊坡土體厚度隨生長時間分別為0.4、0.6、0.8和1.0 m,降水入滲影響深度為1 m,坡度=50°。為了達(dá)到理想計(jì)算精度,坡腳到左端邊界的距離設(shè)置為坡高12 m的1.5倍即18 m,坡肩到右端邊界的距離為坡高12 m的2.5倍即30 m,且上下邊界總高等于坡高12 m的2倍即24 m[36],建立邊坡幾何模型,如圖1所示。模型網(wǎng)格劃分后共包含87 720個六面體單元和96 938個節(jié)點(diǎn)。圖1中①分區(qū)為受到降雨作用和植物根系共同影響的土體,①分區(qū)隨著植物根系的發(fā)育而增大;①分區(qū)與②分區(qū)的分界線為植物根系的發(fā)育深度,②分區(qū)隨著植物根系的發(fā)育而減小,為僅受到降雨作用影響的土體;③分區(qū)為僅受到降雨影響的土體;④分區(qū)為天然狀態(tài)土體。屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則,并將素土和根-土復(fù)合體都視為均質(zhì)材料進(jìn)行模擬分析。模型的側(cè)面和底面設(shè)置為固定邊界,模型的上表面設(shè)置為自由邊界,模型在計(jì)算過程中僅受重力荷載的作用。在模型中設(shè)置了34個監(jiān)測點(diǎn)(各監(jiān)測點(diǎn)位由上而下依次為邊坡表面、①分區(qū)中點(diǎn)、①②分區(qū)界面、②分區(qū)中點(diǎn)、②④分區(qū)界面和坡面豎直向下2 m處;斷面6在邊坡內(nèi)部,故放置4個點(diǎn)位),用于監(jiān)測計(jì)算過程中的水平位移和豎向位移的變化,如圖2所示。

    注:①~④為4個分區(qū)。  Note: ①-④are 4 partitions.

    注:圓點(diǎn)表示每個斷面下各監(jiān)測點(diǎn)(層位),從上到下依次為1~6號點(diǎn),記為1#~6#。

    1.4.3 參數(shù)取值

    重度通過環(huán)刀法測得,抗拉強(qiáng)度σ參考文獻(xiàn)[37]獲得,剪脹角參考文獻(xiàn)[38]取值為0°。不同狀態(tài)下黃土物理力學(xué)指標(biāo)參數(shù)取值如表2所示,不同工況下各分區(qū)參數(shù)選取如表3所示。

    表2 不同狀態(tài)下黃土物理力學(xué)指標(biāo)參數(shù)取值

    注:素黃土與含紫花苜蓿根系黃土的重度值略有波動,此處將其看作定值。:重度;:內(nèi)摩擦角;c:根土復(fù)合體的黏聚力;c:素土的黏聚力;:彈性模量;σ:抗拉強(qiáng)度。

    Note: The gravity values of plain loess and root loess containing alfalfa fluctuated slightly, which was regarded as a fixed value here.: unit weight of soil;: internal friction angle;c: cohesion of root soil complex;c: cohesion of plain soil;: elastic modulus;σ:tensile strength.

    表3 不同工況下各分區(qū)參數(shù)選取

    注:A、B、C、D、E、F指代的為表2中六種狀態(tài);工況6時植物根系深度為1 m,此時沒有②分區(qū)。

    Note: A, B, C, D, E, F refer to the six states in table 2; under condition 6, the root depth was 1 m, and there was no ②partition.

    2 結(jié)果與分析

    2.1 數(shù)值模擬分析結(jié)果

    2.1.1 自重作用下邊坡監(jiān)測點(diǎn)水平和豎向位移

    本文工況2~工況5的計(jì)算結(jié)果不收斂,工況2~工況5在計(jì)算至450時步時塑性區(qū)都發(fā)生貫通,本文以此作為計(jì)算結(jié)束的標(biāo)準(zhǔn),此時邊坡正處于極限平衡狀態(tài)。限于篇幅,圖3和圖4僅展示了模型的監(jiān)測斷面1和坡面層位1的位移情況,能夠反映模型的整體位移變化規(guī)律。從工況1的監(jiān)測點(diǎn)水平位移(圖3a)和豎向位移(圖4a)可以看出:監(jiān)測點(diǎn)水平位移隨時步增加,位移量逐漸增大,且都為負(fù)值,其中邊坡坡肩處(監(jiān)測斷面5)位移值最大;豎向位移各點(diǎn)位移量出現(xiàn)波動,但量值均很??;邊坡坡面與素黃土邊坡整體的位移無明顯差異。工況2~工況5監(jiān)測點(diǎn)的水平位移均為負(fù)值,在邊坡坡腳處監(jiān)測斷面1-1#點(diǎn)的水平位移量最大,從邊坡坡腳到邊坡坡肩各監(jiān)測斷面水平位移量總體上越來越??;邊坡坡腳處由于受到擠壓隆起飽水層或根土復(fù)合體層與無根系飽水層豎向位移量大且為正值,其余監(jiān)測點(diǎn)豎向位移量較??;邊坡根土復(fù)合體層與無根系飽水層位移量無明顯差異,但兩者與天然狀態(tài)土體位移量具有明顯差異,表明邊坡破壞面可能處于邊坡飽水層與天然土體的交界處。但在工況5時,邊坡根土復(fù)合體層與天然狀態(tài)土體位移量具有明顯差異,無根系飽水層位移量與天然狀態(tài)土體位移量差異較小,表明邊坡破壞面可能處于邊坡飽水層。工況6監(jiān)測點(diǎn)的水平位移均為負(fù)值,在靠近邊坡坡腳處監(jiān)測斷面2的水平位移量最大,從邊坡靠近坡腳到邊坡坡肩各監(jiān)測斷面水平位移量總體上越來越小,邊坡坡腳處監(jiān)測斷面1的水平位移量最?。贿吰缕履_處監(jiān)測斷面1的豎向位移為正值,其余監(jiān)測斷面豎向位移均為負(fù)值,且越靠近邊坡坡肩其絕對值越大;坡面表層與素黃土邊坡整體的位移無明顯差異。

    2.1.2 自重作用下邊坡水平位移與塑性區(qū)

    由圖5和圖6可知,工況2~工況5在自重作用下邊坡的最大位移出現(xiàn)在邊坡坡面;隨著計(jì)算進(jìn)行,首先在邊坡坡面以下1 m處即降水最大入滲深度處形成了平行于坡面的剪切屈服塑性區(qū),之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服,并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動。工況6在自重作用下邊坡的最大位移發(fā)生在邊坡坡底面靠近邊坡坡腳處,邊坡整體有向左滑動的趨勢,邊坡坡底面靠近邊坡坡腳處由于擠壓出現(xiàn)隆起,此時邊坡整體處于穩(wěn)定狀態(tài);從塑性區(qū)分布來看,此時僅在邊坡坡底面靠近邊坡坡腳處有少部分網(wǎng)格發(fā)生剪切屈服。在圖5中隨著紫花苜蓿生長期的延長,邊坡的最大水平位移依次為0.193、0.194、0.203、0.261、0.394 cm,都出現(xiàn)在坡面位置,當(dāng)紫花苜蓿生長至150 d即工況6時,最大水平位移僅出現(xiàn)在坡腳位置,坡面水平位移明顯減小,且相對于整體邊坡而言該水平位移量相對較小。

    2.1.3 強(qiáng)度折減法計(jì)算下的最大剪應(yīng)變增量和塑性區(qū)

    工況1采用強(qiáng)度折減法計(jì)算的邊坡穩(wěn)定系數(shù)為1.457,此時邊坡處于穩(wěn)定狀態(tài)。由圖7a可以看出,邊坡的潛在滑動面在坡體內(nèi)部較深的位置,沿著潛在滑動面出現(xiàn)大量的剪切屈服區(qū)域,即圖8a中紅色區(qū)域,邊坡發(fā)生失穩(wěn)破壞時將從坡腳位置剪出。工況2~工況5的邊坡穩(wěn)定系數(shù)分別為0.600、0.637、0.660、0.695,此時邊坡可能會發(fā)生失穩(wěn)破壞。從圖7和圖8可以看出,在邊坡坡腳和邊坡坡肩處均出現(xiàn)剪切屈服區(qū)域,邊坡潛在滑動面發(fā)育在降雨最大入滲深度處,邊坡發(fā)生失穩(wěn)破壞時將從坡腳位置剪出。工況6的邊坡穩(wěn)定系數(shù)為1.367,此時邊坡處于穩(wěn)定狀態(tài)。從圖7f和圖8f可以看出,此時邊坡失穩(wěn)破壞可能有兩個潛在滑動面,一個在降雨入滲能夠影響的最大深度,即坡面以下1 m處,另一個在坡體內(nèi)部較深位置。沿著坡體深部潛在滑動面出現(xiàn)大量的剪切屈服區(qū)域,同時塑性區(qū)已經(jīng)貫通;降雨最大入滲深度的淺部潛在滑動面還未貫通;對于這兩個潛在滑動面,邊坡發(fā)生失穩(wěn)破壞時都將從坡腳位置剪出。隨著紫花苜蓿生長期的延長,邊坡的最大剪應(yīng)變增量依次為0.033、0.156、0.209、0.362、3.829、0.082,當(dāng)紫花苜蓿的生長期達(dá)到150 d時,邊坡的最大剪應(yīng)變增量驟然下降,此時紫花苜蓿對降雨條件下黃土邊坡坡面侵蝕具有明顯的抑制作用。

    注:圖例中每條曲線的第一個數(shù)字代表斷面的編號,第二個數(shù)字代表監(jiān)測點(diǎn)(層位)序號。下同。

    圖4 自重作用下各工況監(jiān)測點(diǎn)豎向位移變化曲線

    圖5 自重作用下各工況450時步的水平位移云圖

    注:圖中紅色區(qū)域?yàn)榧羟星^(qū)。

    圖7 基于強(qiáng)度折減法的各工況最大剪應(yīng)變增量云圖

    圖8 基于強(qiáng)度折減法的各工況塑性區(qū)分布圖

    由數(shù)值模擬結(jié)果發(fā)現(xiàn):工況1時黃土邊坡的潛在滑動面出現(xiàn)在坡體內(nèi)部較深位置處;當(dāng)處于工況2~工況5時,黃土邊坡潛在滑動面出現(xiàn)在降雨入滲能夠影響到的最大深度處,即坡面以下1 m處,表明降雨作用使得黃土邊坡可能發(fā)生失穩(wěn)破壞的潛在滑動面位置由邊坡坡體內(nèi)部轉(zhuǎn)移到邊坡坡面;隨著植物生長至工況6時,黃土邊坡內(nèi)部又開始發(fā)育潛在滑動面。在工況2~工況5時,黃土邊坡的破壞過程都是首先在邊坡坡面以下1 m處,即降水最大入滲深度處,形成了平行于坡面的剪切屈服塑性區(qū),之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服,并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動,且邊坡的最大位移都出現(xiàn)在邊坡坡腳處。因此,黃土邊坡在天然狀態(tài)下是穩(wěn)定的,而在降雨條件下邊坡坡面土體強(qiáng)度迅速降低,重力作用下坡面土體可能沿著降雨入滲的最大影響深度處發(fā)生失穩(wěn)破壞。隨著植被的生長,植物根系的發(fā)育,植物根系對邊坡土體的加固作用顯著增強(qiáng),增大了邊坡的抗剪強(qiáng)度,降低了邊坡的坡面侵蝕作用,阻止了邊坡的淺層破壞,最終提高了邊坡的整體穩(wěn)定性。

    2.2 無限邊坡模型結(jié)果

    將已測得的黏聚力、內(nèi)摩擦角、重度(表2)、坡面與水平地面的夾角=50°以及滑面深度=1.0 m代入無限邊坡模型,得出無限邊坡模型計(jì)算的6種工況下的邊坡穩(wěn)定系數(shù)F值,如表4所示。當(dāng)邊坡處于工況1時邊坡穩(wěn)定系數(shù)大于1,表明邊坡處于穩(wěn)定狀態(tài)。當(dāng)邊坡處于工況2~工況5即受到降雨影響時,邊坡穩(wěn)定系數(shù)迅速減小,處于不穩(wěn)定狀態(tài)。但隨著邊坡坡面紫花苜蓿生長期的延長,邊坡穩(wěn)定系數(shù)逐漸增大。當(dāng)達(dá)到工況6即紫花苜蓿生長150 d時,邊坡處于穩(wěn)定狀態(tài)。

    2.3 數(shù)值模擬與無限邊坡模型結(jié)果對比

    將含不同生長期紫花苜蓿根系的土體的物理力學(xué)參數(shù)值帶入FLAC3D數(shù)值模擬軟件與無限邊坡模型中,探究植物根系對黃土邊坡淺層破壞的防護(hù)效應(yīng)。將無限邊坡模型計(jì)算的邊坡穩(wěn)定系數(shù)與數(shù)值模擬強(qiáng)度折減法的計(jì)算結(jié)果進(jìn)行對比(表4),可以看出:黃土邊坡在工況1與工況6時處于穩(wěn)定狀態(tài),在工況2~工況5時均處于不穩(wěn)定狀態(tài);采用數(shù)值模擬與分析模型兩種方法計(jì)算結(jié)果總體一致。但是,在工況1條件下兩種方法的計(jì)算結(jié)果相差較大,這是由于兩種方法計(jì)算過程中的潛在滑動面位置不同所引起的,數(shù)值模擬強(qiáng)度折減法計(jì)算的潛在滑動面位置在坡體內(nèi)部較深處(圖7),而無限邊坡模型計(jì)算的潛在滑動面假定為坡面以下1 m處。在工況2~工況5條件下,兩種方法的計(jì)算結(jié)果基本一致,這是由于降雨工況下強(qiáng)度折減法計(jì)算的潛在滑動面位置與無限邊坡模型假定的坡面以下1 m處潛在滑動面位置基本一致。對于工況6,數(shù)值模擬計(jì)算結(jié)果表明坡體的潛在滑動面開始從淺部向深部轉(zhuǎn)移,強(qiáng)度折減法計(jì)算的邊坡穩(wěn)定系數(shù)為深部潛在滑動面對應(yīng)的邊坡穩(wěn)定系數(shù),無限邊坡模型仍然以1 m作為潛在滑動面進(jìn)行計(jì)算對應(yīng)的邊坡穩(wěn)定系數(shù),此時兩種計(jì)算方法結(jié)果數(shù)值上很接近,說明當(dāng)植物根系發(fā)育至1 m時,即便是降雨入滲深度也達(dá)到1 m,含植物根系的黃土邊坡整體處于穩(wěn)定狀態(tài),且黃土邊坡淺層破壞也不會發(fā)生。

    數(shù)值模擬與無限邊坡模型兩種方法的結(jié)果表明,黃土邊坡在降雨條件下,其邊坡坡面土體強(qiáng)度迅速降低,邊坡發(fā)生失穩(wěn)破壞,但隨著植物根系的生長發(fā)育,其加固作用顯著增強(qiáng),并增大了黃土邊坡的抗剪強(qiáng)度,降低了黃土邊坡坡面侵蝕作用,阻止了邊坡的淺層局部破壞,最終提高了黃土邊坡的整體穩(wěn)定性。

    表4 數(shù)值模擬與分析模型邊坡穩(wěn)定系數(shù)對比

    3 討 論

    3.1 根系固土?xí)r間效應(yīng)及其數(shù)值分析

    本文與?wita?a等[39-42]的研究結(jié)果一致,都表明了植物根系對邊坡穩(wěn)定性具有加強(qiáng)作用。張興玲等[43]發(fā)現(xiàn)植物根系的生長指標(biāo)與生長時間之間呈現(xiàn)顯著相關(guān)關(guān)系,并通過對生長期為1~4 a的4種不同灌木進(jìn)行現(xiàn)場拉拔試驗(yàn),指出根系平均抗拔力隨生長時間的增加呈顯著增加趨勢。Zhang等[44]通過進(jìn)行原位不排水抗剪強(qiáng)度測試,指出根系數(shù)量、單位體積根截面積與土體抗剪強(qiáng)度增強(qiáng)值均極顯著相關(guān),并且根土復(fù)合體抗剪強(qiáng)度隨樹齡的增大而增大。Valizade等[45]通過直剪試驗(yàn),發(fā)現(xiàn)隨著植物齡期的增長,植物根系與周圍土壤顆粒之間的結(jié)合更好,從而增強(qiáng)了土體的抗剪強(qiáng)度。本文在不同生長期植物根系對邊坡穩(wěn)定性的影響方面,結(jié)論與以上研究一致。

    本文采用數(shù)值模擬方法,研究了不同生長期紫花苜蓿根系對黃土邊坡的護(hù)坡效果,驗(yàn)證了紫花苜蓿隨著生長期的延長,坡面土體的強(qiáng)度逐漸提高,邊坡穩(wěn)定系數(shù)逐漸增大,這與已有文獻(xiàn)[46-47]結(jié)論一致。將根土復(fù)合體層看做一個整體,通過測定根土復(fù)合體試樣的抗剪強(qiáng)度指標(biāo),這與將根系與土分別考慮進(jìn)行數(shù)值模擬計(jì)算有所不同,但所的結(jié)論一致。通過分別測定多花木蘭根系與土體的物理力學(xué)參數(shù),嵇曉雷等[46]采用點(diǎn)對點(diǎn)錨桿單元,模擬不同生長期多花木蘭根系;鄭明新等[47]采用桁架單元模擬多花木蘭根系。本文發(fā)現(xiàn)了紫花苜蓿生長的前5個月為關(guān)鍵期,當(dāng)紫花苜蓿生長時間達(dá)到150 d時,紫花苜蓿對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用。

    總之,植物根系對邊坡穩(wěn)定性具有加強(qiáng)作用;隨著植物生長期的延長,根系固坡的作用也在逐步加強(qiáng);數(shù)值模擬方法能夠有效分析不同生長期植物根系對邊坡淺層破壞的防護(hù)效果。今后將進(jìn)一步開展模型試驗(yàn)或現(xiàn)場監(jiān)測分析不同生長期植物根系對邊坡淺層破壞的防護(hù)效果。

    3.2 數(shù)值模型中的參數(shù)取值合理性分析

    本文通過三軸壓縮試驗(yàn)測定了含紫花苜蓿根系黃土試樣的黏聚力、內(nèi)摩擦角、彈性模量和泊松比4個參數(shù)。由于植物根系的生長不可避免會存在上下部根密度、根直徑的差異,三軸試驗(yàn)試樣取樣處的根土復(fù)合體強(qiáng)度對于埋深較大的根系分布范圍內(nèi)土樣確實(shí)相對較高。本文考慮降雨強(qiáng)度的最不利因素,即邊坡表層1 m內(nèi)完全飽和,同時在相對較高的土體強(qiáng)度下邊坡坡面已然都發(fā)生了破壞,在強(qiáng)度較低的埋深較大的根系分布范圍內(nèi)發(fā)生破壞更容易。另外,本文模型計(jì)算是針對邊坡整體的穩(wěn)定性,文章研究內(nèi)容重點(diǎn)關(guān)注的是邊坡表層1 m,相對整體邊坡的幾何尺寸,邊坡表層僅1 m幾何尺寸相差較大,因此對于整體邊坡的穩(wěn)定性分析并沒有太大影響。若再將邊坡表層1 m進(jìn)行細(xì)分,不僅會增加數(shù)值模擬模型建立的復(fù)雜度,數(shù)值模擬模型的建立并非越精細(xì)越好,適當(dāng)簡化模型能夠更加清楚地反映計(jì)算結(jié)果的規(guī)律[48-49],并且節(jié)約了時間成本。

    泊松比是側(cè)向應(yīng)變增量與豎向應(yīng)變增量的比值。曾超峰等[50]通過開展基坑抽水模型試驗(yàn)與三維數(shù)值分析,引用其他文獻(xiàn)的材料參數(shù),將圍護(hù)結(jié)構(gòu)、內(nèi)隔墻、降水井的泊松比統(tǒng)一設(shè)定為0.2,探究內(nèi)隔墻對于限制開挖前抽水引發(fā)基坑變形的有效性。常成等[51]開發(fā)了一個概念性的熱-孔隙流體-應(yīng)力應(yīng)變耦合的二維平面應(yīng)變有限元模型,將圍巖、侵入巖體和花崗巖基底的泊松比統(tǒng)一定為0.25,并模擬了斑巖成礦系統(tǒng)的溫度、應(yīng)力、應(yīng)變、孔隙流體壓力的演化及巖漿侵入體周圍的破裂形成過程。楊衛(wèi)波等[52]為了獲得熱力耦合作用下相變混凝土能量樁的熱-力學(xué)特性,建立了其三維數(shù)值模型,引用其他文獻(xiàn)的材料參數(shù),將泊松比設(shè)定為定值,比較了傳統(tǒng)和相變混凝土能量樁熱-力學(xué)特性的差異。本文泊松比參數(shù)為三軸試驗(yàn)實(shí)測獲得,隨著植物生長期的延長,植物根系逐漸生長發(fā)育,根系固土作用逐漸增強(qiáng),試樣的側(cè)向應(yīng)變增量逐漸減小,導(dǎo)致泊松比的逐漸減小。左學(xué)龍[53]采用GDS非飽和土試驗(yàn)儀,針對四類不同含根量的根-土復(fù)合體試樣,分別進(jìn)行了常規(guī)三軸壓縮排水試驗(yàn)和等向固結(jié)試驗(yàn)等,基于鄧肯-張模型和K-G模型理論,探索了含根量對根-土復(fù)合體力學(xué)特性的影響規(guī)律,發(fā)現(xiàn)隨著含根量的增大,根土復(fù)合體的泊松比減小。而其他論文多將泊松比設(shè)置為一定值,如鄭明新等[47]將不同生長期多花木蘭根系的泊松比參數(shù)統(tǒng)一設(shè)置為0.2,采用二維有限元法模擬試驗(yàn)區(qū)不同生長期多花木蘭根系對邊坡穩(wěn)定性的增強(qiáng)作用。

    3.3 黃土邊坡淺層植物防護(hù)的時間效應(yīng)

    中國西北黃土高原地處干旱-半干旱地區(qū),地質(zhì)生態(tài)環(huán)境十分脆弱,水土流失嚴(yán)重,滑坡、崩塌災(zāi)害頻發(fā),給當(dāng)?shù)氐墓まr(nóng)業(yè)生產(chǎn)以及人民的生命財(cái)產(chǎn)安全帶來巨大的威脅[54]。植物護(hù)坡是恢復(fù)黃土邊坡生態(tài)、減少土壤侵蝕和防治黃土邊坡地質(zhì)災(zāi)害的有效手段[55]。

    在開展植被護(hù)坡措施的同時,研究人員已經(jīng)關(guān)注到植物護(hù)坡的時間效應(yīng)問題,如張興玲等[43,56]為研究植物護(hù)坡的時間效應(yīng),對生長期為1~4 a的四種不同灌木進(jìn)行現(xiàn)場拉拔試驗(yàn),指出根系平均抗拔力隨生長時間的增加呈顯著增加趨勢;陳潮[57]從時間尺度,對護(hù)坡灌木根系與土體之間相互作用力學(xué)機(jī)理展開研究,發(fā)現(xiàn)邊坡淺層穩(wěn)定性隨著時間的推移穩(wěn)步提升,且植被防護(hù)效果越來越好;劉治興等[58]通過對素土及植物生長期61、95、125 d的根土復(fù)合體抗剪強(qiáng)度進(jìn)行現(xiàn)場測定和物理模型預(yù)測,發(fā)現(xiàn)抗剪強(qiáng)度隨植物生長階段呈逐漸增大趨勢。

    當(dāng)紫花苜蓿從120 d生長至150 d時,邊坡穩(wěn)定系數(shù)從0.695增大至1.367。究其原因,首先工況5在坡面處分布有3層,分別為飽和根土復(fù)合體層、飽和素黃土層和天然素黃土層,工況6僅有飽和根土復(fù)合體層和天然素黃土層。當(dāng)邊坡處于工況5時,由于飽和素黃土的強(qiáng)度較低,在整個黃土邊坡中成為了軟弱層帶,軟弱層帶通常是邊(斜)坡失穩(wěn)的控滑帶[59],因此其邊坡穩(wěn)定系數(shù)較低。當(dāng)處于工況6時,植物根系生長至1 m,根系的加固作用大大增強(qiáng)了邊坡軟弱層帶的強(qiáng)度,提高了黃土邊坡的穩(wěn)定系數(shù)。

    黃土地區(qū)的雨季一般在7月–9月,雨季坡面侵蝕現(xiàn)象多發(fā)。根據(jù)本文的研究結(jié)果表明:當(dāng)生長時間達(dá)到150 d時,草本植物根系能夠有效防治黃土坡面侵蝕,并保持黃土邊坡穩(wěn)定性。并且隨著植物生長期的延長,植物根系對邊坡土體的固土效果更加顯著[46]。因此,草本植物生長的前5個月為關(guān)鍵期,假定對裸露黃土邊坡于春季(3月或4月)進(jìn)行播種,則同年雨季來臨之時,坡面植被能夠起到防治降雨條件下的坡面侵蝕的作用,若為其他時間播種,在前5個月時需要進(jìn)行一定的人為保護(hù),如鋪麥草、塑料布等,對裸露坡面進(jìn)行保護(hù)。坡面侵蝕普遍發(fā)育面積大,采用傳統(tǒng)工程措施進(jìn)行防治需要龐大的資金支持且容易造成資源浪費(fèi),相比而言造價低廉的植被種植技術(shù)具有廣闊的推廣前景,這一技術(shù)在黃土邊坡淺層破壞防治的同時,還兼顧了景觀效應(yīng)和生態(tài)可持續(xù)發(fā)展。

    4 結(jié) 論

    本文考慮降雨入滲條件并以含不同生長期紫花苜蓿的黃土邊坡為例,分別采用無限邊坡模型和數(shù)值模擬兩種方法進(jìn)行了穩(wěn)定性分析,探究了紫花苜蓿對黃土邊坡淺層破壞防護(hù)的時間效應(yīng)。本文得出以下結(jié)論:

    1)不同工況下黃土邊坡淺層破壞過程都是首先在邊坡坡面以下1 m處,即降水最大入滲深度處,形成了平行于坡面的剪切屈服塑性區(qū),之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服,并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動,且邊坡的最大位移都出現(xiàn)在邊坡坡腳處。

    2)黃土邊坡在天然條件下處于穩(wěn)定狀態(tài),潛在滑動面在坡體內(nèi)部較深位置,降雨作用使得邊坡潛在滑動面位置由坡體內(nèi)部轉(zhuǎn)移到坡體降雨最大入滲深度處,在植物根系的加固作用下,坡面位置的穩(wěn)定性提高,此時邊坡潛在滑動面又轉(zhuǎn)為內(nèi)部。

    3)草本植物隨著生長期的延長,坡面土體的強(qiáng)度逐漸提高。數(shù)值模擬強(qiáng)度折減法計(jì)算的邊坡穩(wěn)定系數(shù)隨著植被生長不斷增大,表明植被對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用,在坡面上種植植物能夠有效的防止黃土邊坡在雨季降雨事件時發(fā)生黃土邊坡淺層破壞。

    4)通過對比無限邊坡模型與數(shù)值模擬方法計(jì)算的邊坡穩(wěn)定系數(shù),發(fā)現(xiàn)隨著植物生長期的延長,采用數(shù)值模擬與無限邊坡模型兩種方法計(jì)算的邊坡穩(wěn)定系數(shù),總體變化趨勢一致,表現(xiàn)出草本植物對黃土邊坡淺層破壞防護(hù)具有時間效應(yīng);黃土邊坡在天然狀態(tài)和降雨條件下植物生長期為150 d時穩(wěn)定系數(shù)大于1,處于穩(wěn)定狀態(tài),在降雨條件下植物生長期為0~120 d時穩(wěn)定系數(shù)小于1,都處于欠穩(wěn)定狀態(tài)。

    5)當(dāng)草本植物生長時間達(dá)到150 d時,植物根系生長能夠有效防治黃土坡面侵蝕,并提高黃土邊坡穩(wěn)定性,因此草本植物生長的前5個月為關(guān)鍵期。

    [1] 王蘭民,蒲小武,吳志堅(jiān),等. 地震和降雨耦合作用下黃土邊坡失穩(wěn)滑移的振動臺試驗(yàn)研究[J]. 巖石力學(xué)與工程學(xué)報(bào),2017,36(增刊2):3873-3883.

    Wang Lanmin, Pu Xiaowu, Wu Zhijian, et al. The shaking table test of the instability sliding of loess slope under the coupling effects of earthquake and rainfall[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(Supp.2): 3873-3883. (in Chinese with English abstract)

    [2] 孫軍杰,王蘭民,龍鵬偉,等. 地震與降雨耦合作用下區(qū)域滑坡災(zāi)害評價方法[J]. 巖石力學(xué)與工程學(xué)報(bào),2011,30(4):752-760.

    Sun Junjie, Wang Lanmin, Long Pengwei, et al. An assessment method for regional susceptibility of landslides under coupling condition of earthquake and rainfall[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(4): 752-760. (in Chinese with English abstract)

    [3] 王鮮,倪萬魁,劉海松,等. 黃土邊坡坡面降雨沖刷破壞機(jī)制模擬試驗(yàn)[J]. 中國地質(zhì)災(zāi)害與防治學(xué)報(bào),2015,26(3):41-45.

    Wang Xian, Ni Wankui, Liu Haisong, et al. Simulation experiment on loess slope rainfall erosion mechanism[J]. The Chinese Journal of Geological Hazard and Control, 2015, 26(3): 41-45. (in Chinese with English abstract)

    [4] Tu X B, Kwong A K L, Dai F C, et al. Field monitoring of rainfall infiltration in a loess slope and analysis of failure mechanism of rainfall-induced landslides[J]. Engineering Geology, 2009, 105(1/2): 134-150.

    [5] 張茂省,李同錄. 黃土滑坡誘發(fā)因素及其形成機(jī)理研究[J]. 工程地質(zhì)學(xué)報(bào),2011,19(4):530-540.

    Zhang Maosheng, Li Tonglu. Triggering factors and forming mechanism of loess landslides[J]. Journal of Engineering Geology, 2011, 19(4): 530-540. (in Chinese with English abstract)

    [6] 董凌勃,海旭瑩,汪曉珍,等. 黃土高原退耕還草地植物群落動態(tài)對生態(tài)系統(tǒng)碳儲量的影響[J]. 生態(tài)學(xué)報(bào),2020,40(23):8559-8569.

    Dong Linbo, Hai Xuying, Wang Xiaozhen, et al. Effects of plant community dynamics on ecosystem carbon stocks since returning farmlands to grasslands on the Loess Plateau[J]. Acta Ecologica Sinica, 2020, 40(23): 8559-8569. (in Chinese with English abstract)

    [7] 邱海軍,崔鵬,曹明明,等. 基于最大熵原理的黃土丘陵區(qū)地質(zhì)災(zāi)害規(guī)模頻率分布研究[J]. 巖土力學(xué),2014,35(12):3541-3549,3555.

    Qiu Haijun, Cui Peng, Cao Mingming, et al. Distribution of frequency-size of geological disaster base on principle of maximum entropy in loess plateau[J]. Rock and Soil Mechanics, 2014, 35(12): 3541-3549, 3555. (in Chinese with English abstract)

    [8] 殷躍平,胡時友,石勝偉,等. 滑坡防治技術(shù)指南[M]. 北京:地質(zhì)出版社,2018:267-268.

    [9] 姚仰平,張丙印,朱俊高. 土的基本特性、本構(gòu)關(guān)系及數(shù)值模擬研究綜述[J]. 土木工程學(xué)報(bào),2012,45(3):127-150.

    Yao Yangping, Zhang Bingyin, Zhu Jungao. Behaviors, constitutive models and numerical simulation of soils[J]. China Civil Engineering Journal, 2012, 45(3): 127-150. (in Chinese with English abstract)

    [10] 張丙印,溫彥鋒,朱本珍,等. 土工構(gòu)筑物和邊坡工程發(fā)展綜述:作用機(jī)理與數(shù)值模擬方法[J]. 土木工程學(xué)報(bào),2016,49(8):1-15,35.

    Zhang Bingyin, Wen Yanfeng, Zhu Benzhen, et al. Action mechanism and numerical simulation methods of soil structures and slope engineering: the state-of-the-art[J]. China Civil Engineering Journal, 2016, 49(8): 1-15, 35. (in Chinese with English abstract)

    [11] 姚成成. 土質(zhì)淺層邊坡的草本植物防護(hù)加固效應(yīng)分析[D]. 成都:成都理工大學(xué),2017.

    Yao Chengcheng. Analysis of Herb Protection and Reinforcement Effect of Soil Shallow Slope[D]. Chengdu: Chengdu University of Technology, 2017. (in Chinese with English abstract)

    [12] 周群華,鄧衛(wèi)東. 植物根系固坡的有限元數(shù)值模擬分析[J]. 公路,2007(12):132-136.

    Zhou Qunhua, Deng Weidong. Finite element numerical simulation analysis of slope reinforcement by plant roots[J]. Highway, 2007(12): 132-136. (in Chinese with English abstract)

    [13] Temgoua A G T, Kokutse N K, Kavazovic Z, et al. 3D model applied to analyze the mechanical stability of real-world forested hillslopes prone to landslides[J]. Ecological Engineering, 2017, 106: 609-619.

    [14] Li Y P, Wang Y Q, Ma C, et al. Influence of the spatial layout of plant roots on slope stability[J]. Ecological Engineering, 2016, 91: 477-486.

    [15] 黃建坤,王學(xué)林,及金楠,等. 基于漸近均勻化理論的黃土高原草本植物固土效果模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(9):168-176.

    Huang Jiankun, Wang Xuelin, Ji Jinnan, et al. Numerical simulation of root reinforcement for herbs in Loess Plateau based on asymptotic homogenization theory[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(9): 168-176. (in Chinese with English abstract)

    [16] Fattet M, Fu Y, Ghestem M, et al. Effects of vegetation type on soil resistance to erosion: Relationship between aggregate stability and shear strength[J]. Catena, 2011, 87(1): 60-69.

    [17] Stokes A, Douglas G B, Fourcaud T, et al. Ecological mitigation of hillslope instability: Ten key issues facing researchers and practitioners[J]. Plant Soil, 2014, 377(1): 1-23.

    [18] 楊吉華,張光燦,劉霞,等. 紫花苜蓿保持水土效益的研究[J]. 土壤侵蝕與水土保持學(xué)報(bào),1997,3(2):91-96.

    Yang Jihua, Zhang Guangcan, Liu Xia, et al. Study on benefits of alfalfa conservation soil and water[J]. Journal of Soil Erosion and Soil and Water Conservation, 1997, 3(2): 91-96. (in Chinese with English abstract)

    [19] 趙祥,朱慧森. 牧草科學(xué)種植與利用[M]. 北京:中國社會出版社,2006.

    [20] 李廣信. 高等土力學(xué)[M]. 北京:清華大學(xué)出版社,2004.

    [21] 雷小芹,劉恩龍,何思明,等. 非飽和堆積土邊坡降雨-滲流潛蝕耦合過程模擬[J]. 工程科學(xué)與技術(shù),2021,53(2):28-37.

    Lei Xiaoqin, Liu Enlong, He Siming, et al. Modelling of coupled process of rainfall infiltration and internal erosion within unsaturated deposited soil slopes[J]. Advanced Engineering Sciences, 2021, 53(2): 28-37. (in Chinese with English abstract)

    [22] 金保明,林鵬,李光敦. 基于運(yùn)動波壤中流理論與無限邊坡穩(wěn)定分析理論的滑坡預(yù)測模型[J]. 工程科學(xué)與技術(shù),2021,53(4):140-148.

    Jin Baoming, Lin Peng, Lee Kwantun. Integrating kinematic subsurface-wave approximation with infinite-slope stability analysis for landslide predictions[J]. Advanced Engineering Sciences, 2021, 53(4): 140-148. (in Chinese with English abstract)

    [23] 吳宏偉. 大氣-植被-土體相互作用:理論與機(jī)理[J]. 巖土工程學(xué)報(bào),2017,39(1):1-47.

    Wu Hongwei. Atmosphere-plant-soil interactions: Theories and mechanisms[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(1): 1-47. (in Chinese with English abstract)

    [24] Arnone E, Caracciolo D, Noto L V, et al. Modeling the hydrological and mechanical effect of roots on shallow landslides[J]. Water Resources Research, 2016, 52(11): 8590-8612.

    [25] Hayati E, Abdi E, Saravi M M, et al. Soil water dynamics under different forest vegetation cover: Implications for hillslope stability[J]. Earth Surface Processes and Landforms, 2018, 43(10): 2106-2120.

    [26] Leung F T Y , Yan W M, Hau B C H, et al. Root systems of native shrubs and trees in Hong Kong and their effects on enhancing slope stability[J]. Catena, 2015, 125: 102-110.

    [27] Waldron L J. The shear resistance of root-permeated homogeneous and stratified soil[J]. Soil Science Society of America Journal, 1977, 41(5): 843-849

    [28] Wu T H, McKinnell W P, Swanston D N. Strength of tree roots and landslides on Price of Wales Island, Alaska[J]. Canadian Geotechnical Journal, 1979, 16(1): 19-33.

    [29] 李國良,邵生俊,靳寶成,等. 黃土隧道地基的濕陷性問題研究[J]. 鐵道工程學(xué)報(bào),2015,32(12):12-16,50.

    Li Guoliang, Shao Shengjun, Jin Baocheng, et al. Research on the Problems of collapsibility of loess tunnel foundation[J]. Journal of Railway Engineering Society, 2015, 32(12): 12-16, 50. (in Chinese with English abstract)

    [30] 邵生俊,楊春鳴,焦陽陽,等. 濕陷性黃土隧道的工程性質(zhì)分析[J]. 巖土工程學(xué)報(bào),2013,35(9):1580-1590.

    Shao Shengjun, Yang Chunming, Jiao Yangyang, et al. Engineering properties of collapsible loess tunnel[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(9): 1580-1590. (in Chinese with English abstract)

    [31] 白文明,左強(qiáng),李保國. 烏蘭布和沙區(qū)紫花苜蓿根系吸水模型[J]. 植物生態(tài)學(xué)報(bào),2001,25(4):431-437.

    Bai Wenming, Zuo Qiang, Li Baoguo. A simulation model for water uptake by alfalfa roots in the Wulanbuhe sandy desert[J]. Acta Phytoecologica Sinica, 2001, 25(4): 431-437. (in Chinese with English abstract)

    [32] Tu X B, Kwong A K L, Dai F C, et al. Field monitoring of rainfall infiltration in a loess slope and analysis of failure mechanism of rainfall-induced landslides[J]. Engineering Geology, 2009, 105(1/2): 134-150.

    [33] Mao Z J, Wang X K, An N, et al. Water leakage susceptible areas in loess multi-arch tunnel operation under the lateral recharge conditions[J]. Environmental Earth Sciences, 2020, 79(15): 1-32.

    [34] Li Z Q, Oyediran I A, Tang C, et al. FEM application to loess slope excavation and support: Case study of Dong Loutian coal bunker, Shuozhou, China[J]. Bulletin of Engineering Geology and the Environment, 2014, 73(4): 1013-1023.

    [35] 趙秋. 北方低產(chǎn)土壤實(shí)用綠肥作物栽培與利用[M]. 天津:天津科技翻譯出版公司,2010.

    [36] 張魯渝,鄭穎人,趙尚毅,等. 有限元強(qiáng)度折減系數(shù)法計(jì)算土坡穩(wěn)定安全系數(shù)的精度研究[J]. 水利學(xué)報(bào),2003(1):21-27.

    Zhang Luyu, Zheng Yingren, Zhao Shangyi, et al. The feasibility study of strength-reduction method with FEM for calculating safety factors of soil slope stability[J]. Journal of Hydraulic Engineering, 2003(1): 21-27. (in Chinese with English abstract)

    [37] 趙勇,李國良,喻渝. 黃土隧道工程[M]. 北京:中國鐵道出版社,2011:32.

    [38] Yin Y P, Li B, Wang W P. Dynamic analysis of the stabilized Wangjiayan landslide in the Wenchuan Ms 8.0 earthquake and aftershocks[J]. Landslides, 2015, 12(3): 537-547.

    [39] ?wita?a B M, Wu W. Numerical modelling of rainfall- induced instability of vegetated slopes[J]. Géotechnique, 2018, 68(6): 481-491.

    [40] ?wita?a B M. Numerical simulations of triaxial tests on soil-root composites and extension to practical problem: rainfall-induced landslide[J]. International Journal of Geomechanics, 2020, 20(11): 04020206.

    [41] Park J, Kim I, Kang J K. Root Reinforcement effect on cover slopes of solid waste landfill in soil bioengineering[J]. Sustainability, 2021, 13(7): 3991.

    [42] Eab K H, Likitlersuang S, Takahashi A. Laboratory and modelling investigation of root-reinforced system for slope stabilization[J]. Soils and Foundations, 2015, 55(5): 1270-1281.

    [43] 張興玲,胡夏嵩,李國榮,等. 青藏高原東北部黃土區(qū)灌木幼林根系護(hù)坡的時間效應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(4):136-141.

    Zhang Xingling, Hu Xiasong, Li Guorong, et al. Time effect of young shrub roots on slope protection of loess area in Northeast Qinghai-Tibetan plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(4): 136-141. (in Chinese with English abstract)

    [44] Zhang D, Cheng J H, Liu Y, et al. Spatio-temporal dynamic architecture of living brush mattress:root system and soil shear strength in riverbanks[J]. Forests, 2018, 9(8): 493.

    [45] Valizade N, Tabarsa A. Laboratory investigation of plant root reinforcement on the mechanical behaviour and collapse potential of loess soil[J]. European Journal of Environmental and Civil Engineering, 2020, 26(4): 1475-1491.

    [46] 嵇曉雷,尹鵬. 不同生長時間根系對邊坡穩(wěn)定性的影響[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào),2020,40(12):144-150,178.

    Ji Xiaolei, Yin Peng. Influence of different root growth stages on the slope anti-erosion performance[J]. Journal of Central South University of Forestry & Technology, 2020, 40(12): 144-150, 178. (in Chinese with English abstract)

    [47] 鄭明新,黃鋼,彭晶. 不同生長期多花木蘭根系抗拉拔特性及其根系邊坡的穩(wěn)定性[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(20):175-182.

    Zheng Mingxin, Huang Gang, Peng Jing. Tensile-pullout properties of roots of Magnolia multiflora in different growth stages and stability of slope with its root[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(20): 175-182. (in Chinese with English abstract)

    [48] 朱家奇,朱晟. 單元網(wǎng)格密度及畸變度對面板壩有限元計(jì)算精度的影響[J]. 水電能源科學(xué),2022,40(6):88-92.

    Zhu Jiaqi, Zhu Sheng. Influence of element mesh density and distortion degree on finite element calculation accuracy of face rockfill dams[J]. Water Resources and Power, 2022, 40(6): 88-92. (in Chinese with English abstract)

    [49] 汪必升,李毅波,廖雅詩,等. 基于擴(kuò)展有限元模型的動態(tài)應(yīng)力強(qiáng)度因子計(jì)算[J]. 中國機(jī)械工程,2019,30(11):1294-1301.

    Wang Bisheng, Li Yibo, Liao Yashi, et al. Calculation of dynamic stress intensity factors based on XFFM model[J]. China Mechanical Engineering, 2019, 30(11): 1294-1301. (in Chinese with English abstract)

    [50] 曾超峰,王碩,宋偉煒,等. 內(nèi)隔墻對開挖前抽水引發(fā)軟土區(qū)地鐵深基坑變形的控制效果[J]. 巖石力學(xué)與工程學(xué)報(bào),2021,40(6):1277-1286.

    Zeng Chaofeng, Wang Shuo, Song Weiwei, et al. Control effect of cross wall on metro foundation pit deformation induced by pre-excavation dewatering in soft soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1277-1286. (in Chinese with English abstract)

    [51] 常成,羅綱. 斑巖礦床侵入體頂部破裂系統(tǒng)形成的力學(xué)機(jī)制:多場耦合數(shù)值模擬的啟示[J]. 地球物理學(xué)報(bào),2022,65(8):3006-3024.

    Chang Cheng, Luo Gang. The fracture forming mechanism at the top of the intrusion body in porphyry deposit: The insights from a multi-physical field coupled numerical model[J]. Chinese Journal of Geophysics, 2022, 65(8): 3006-3024. (in Chinese with English abstract)

    [52] 楊衛(wèi)波,楊彬彬,汪峰. 相變混凝土能量樁熱-力學(xué)特性的數(shù)值模擬與試驗(yàn)驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(2):268-277.

    Yang Weibo, Yang Binbin, Wang Feng. Numerical simulation and experimental validation of the thermo-mechanical characteristics of phase change concrete energy pile[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(2): 268-277. (in Chinese with English abstract)

    [53] 左學(xué)龍. 基于彈性本構(gòu)關(guān)系的根-土復(fù)合體力學(xué)特性試驗(yàn)研究[D]. 長沙:中南林業(yè)科技大學(xué),2016.

    Zuo Xuelong. Experimental Study on Mechanical Properties of Root Soil Composite Based on Elastic Constitutive Relationship[D]. Changsha: Central South University of Forestry Science and Technology, 2016. (in Chinese with English abstract)

    [54] 彭建兵,王啟耀,門玉明,等. 黃土高原滑坡災(zāi)害[M]. 北京:科學(xué)出版社,2019.

    [55] Li Q, Liu G B, Zhang Z, et al. Effect of root architecture on structural stability and erodibility of topsoils during concentrated flow in hilly Loess Plateau[J]. Chinese Geographical Science, 2015, 25(6): 757-764.

    [56] 張興玲,胡夏嵩,李國榮,等. 寒旱環(huán)境草本植物根系護(hù)坡的時間尺度效應(yīng)[J]. 水文地質(zhì)工程地質(zhì),2009,36(4):117-120.

    Zhang Xingling, Hu Xiasong, Li Guorong, et al. Time scale effect of herb roots on slope protection in a cold, arid and semi-arid environment[J]. Hydrogeology and Engineering Geology, 2009, 36(4): 117-120. (in Chinese with English abstract)

    [57] 陳潮. 護(hù)坡灌木根—土力學(xué)作用時間尺度效應(yīng)研究[D]. 成都:西南交通大學(xué),2017.

    Chen Chao. Study on Time Scale Effect of Shrub Root Soil Mechanical Action on Slope Protection[D]. Chengdu: Southwest Jiaotong University, 2017. (in Chinese with English abstract)

    [58] 劉治興,楊建英,楊陽,等. 高速公路不同植物防護(hù)邊坡根土復(fù)合體抗剪能力研究[J]. 生態(tài)環(huán)境學(xué)報(bào),2015,24(4):631-637.

    Liu Zhixing, Yang Jianying, Yang Yang, et al. Research on the root-soil composite anti-shearing strength on expressway slope in different plant protection measures[J]. Ecology and Environmental Sciences, 2015, 24(4): 631-637. (in Chinese with English abstract)

    [59] 崔圣華,裴向軍,黃潤秋. 大光包滑坡啟動機(jī)制:強(qiáng)震過程滑帶非協(xié)調(diào)變形與巖體動力致?lián)p[J]. 巖石力學(xué)與工程學(xué)報(bào),2019,38(2):237-253.

    Cui Shenghua, Pei Xiangjun, Huang Runqiu. An initiation model of DGB landslide: Non-coordinated deformation inducing rock damage in sliding zone during strong seismic shaking[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38 (2): 237-253.

    Numerical analysis of protection time effect on planting alfalfa in loess slope with shallow failure

    Mao Zhengjun1,2,3, Zhang Jinge1, Bi Yinli1, Sun Weibo4, An Ning1

    (1.,,710054,; 2.,,710054,; 3.,710054,; 4.,,710054,)

    The special microstructure of loess has often posed a great risk on the slope in the loess areas. The loess structure can suddenly collapse, due to the complete loss of strength, particularly when the loess is wetted by water or subjected to a strong earthquake. The loess slope can usually be triggered by the slope surface erosion under the action of rainfall. The local failure of the slope body can often occur under heavy rainfall. There is a serious threat to the overall stability of the slope. For instance, hundreds of slope shallow failures have occurred in the Chinese Loess Plateau every year, as one of the regions with the thickest loess deposition in the world. The loess slope shallow failure can cause the loss and destruction of land sustainable productivity in the loess region. Fortunately, the slope protection with vegetation can be expected to alleviate the current slope shallow failure. A reasonable numerical simulation has been used to analyze the mechanism of root reinforcement soil, providing a strong theoretical basis for the application of slope protection with vegetation. This study aims to explore the protection time effect on planting alfalfa on the loess slope with shallow failure. The example was set as the loess slope with the alfalfa in different growth periods (0, 60, 90, 120, and 150 days). The rainfall infiltration was also considered in this case. A series of tests were carried out using the indoor planting of PVC pipe soil samples, including the quick direct shear tests, Consolidated Undrained (CU), and Unconsolidated Undrained (UU) triaxial tests. The physical and mechanical parameters were measured for the root-soil composite under different conditions. The slope stability was analyzed using the infinite slope model with the plant roots and the numerical simulation. The results show that three stages were divided for the loess slope under different conditions in the process of shallow failure. Firstly, a shear yield plastic zone parallel to the slope surface was formed at the maximum infiltration depth of precipitation. Secondly, there was a small amount of soil in the middle and upper position of the slope surface, due to the tension action. Thirdly, the shear yield zone gradually moved along the slope surface toward the slope shoulder. The loess slope was in a stable state under natural conditions. The potential sliding surface was at a deep position inside the slope. Therefore, the rainfall infiltration was attributed to the potential sliding surface of the slope transferring from the inside of the slope to the maximum infiltration depth of rainfall. The soil strength gradually increased in the shallow portion of the slope with the growth of vegetation on the slope. The potential sliding surface of the slope gradually shifted to the inside of the slope. The slope stability coefficient increased with the growth of vegetation, according to the strength reduction of numerical simulation. A comparison was made to determine the slope stability coefficient from the infinite slope model and the numerical simulation. The loess slope was in a stable state, when the plant roots grew to 1.0 m, either under natural or rainfall conditions, whereas, in an unstable state, when the plant roots grew to 0-0.8 m under rainfall conditions. The slope stability coefficient went down and then gradually increased with the growth of vegetation. The overall trend was all the same. The herbaceous plants were used to significantly inhibit the loess slope surface erosion under rainfall conditions, particularly with the extension of the growing period. Consequently, the herbaceous roots can be expected to improve the stability of the loess slope. Once the growth time of herbaceous plants reached 150 days, the loess slope surface erosion was effectively prevented for the better stability of the loess slope. Therefore, the first five months of herbaceous plant growth can be the key period to improving the stability of the loess slope in this case. The finding can provide the theoretical reference and application for the shallow failure protection of the loess slope in the sustainable development of the loess areas.

    slope stability; numerical simulation; loess;alfalfa; shallow slope failure; time effect; infinite slope model

    10.11975/j.issn.1002-6819.2022.15.008

    S157.2

    A

    1002-6819(2022)-15-0072-12

    毛正君,張瑾鴿,畢銀麗,等. 紫花苜蓿對黃土邊坡淺層破壞防護(hù)時間效應(yīng)的數(shù)值分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(15):72-83.doi:10.11975/j.issn.1002-6819.2022.15.008 http://www.tcsae.org

    Mao Zhengjun, Zhang Jinge, Bi Yinli, et al. Numerical analysis of protection time effect on planting alfalfa in loess slope with shallow failure[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(15): 72-83. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.15.008 http://www.tcsae.org

    2022-06-15

    2022-07-30

    陜西省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2020SF-379);寧夏回族自治區(qū)重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2022BEG03059,2022BEG03052)

    毛正君,博士,副教授,研究方向?yàn)榈刭|(zhì)環(huán)境保護(hù)與國土空間生態(tài)修復(fù)。Email:zhengjun_mao@163.com。

    猜你喜歡
    植物
    誰是最好的植物?
    為什么植物也要睡覺
    長得最快的植物
    各種有趣的植物
    植物也會感到痛苦
    會喝水的植物
    植物的防身術(shù)
    把植物做成藥
    哦,不怕,不怕
    將植物穿身上
    首页视频小说图片口味搜索| 国产高清videossex| 亚洲av成人精品一区久久| 亚洲av五月六月丁香网| 国产主播在线观看一区二区| 欧美激情在线99| 成年人黄色毛片网站| 日本在线视频免费播放| 变态另类成人亚洲欧美熟女| 国内揄拍国产精品人妻在线| 亚洲aⅴ乱码一区二区在线播放| 两人在一起打扑克的视频| 成人欧美大片| 日韩精品青青久久久久久| 亚洲18禁久久av| 免费看十八禁软件| 九色国产91popny在线| 国产乱人伦免费视频| 亚洲成人精品中文字幕电影| 精品免费久久久久久久清纯| 亚洲avbb在线观看| 婷婷亚洲欧美| 色综合亚洲欧美另类图片| 女人十人毛片免费观看3o分钟| 欧美丝袜亚洲另类 | 午夜免费激情av| 老司机福利观看| 人人妻人人看人人澡| 久久久久久九九精品二区国产| 男女床上黄色一级片免费看| 男插女下体视频免费在线播放| 国内久久婷婷六月综合欲色啪| 亚洲国产精品合色在线| 久久99热这里只有精品18| 亚洲国产精品成人综合色| 免费看日本二区| 亚洲中文日韩欧美视频| 日韩精品青青久久久久久| 成人特级av手机在线观看| 亚洲人成电影免费在线| 亚洲午夜理论影院| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 可以在线观看的亚洲视频| 99热精品在线国产| 国产精品久久久久久亚洲av鲁大| 国产一区二区三区在线臀色熟女| 亚洲欧美激情综合另类| 天天一区二区日本电影三级| 九九久久精品国产亚洲av麻豆| 啪啪无遮挡十八禁网站| 日韩欧美在线二视频| 亚洲av五月六月丁香网| 午夜免费观看网址| 99在线视频只有这里精品首页| 99热这里只有是精品50| 久久久久久人人人人人| 亚洲无线在线观看| 亚洲精品456在线播放app | 国产伦精品一区二区三区四那| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久| 久99久视频精品免费| 麻豆一二三区av精品| 中文字幕av成人在线电影| 琪琪午夜伦伦电影理论片6080| 黄色片一级片一级黄色片| 99久久精品一区二区三区| 国产精品综合久久久久久久免费| 午夜福利免费观看在线| 中出人妻视频一区二区| 久久国产精品影院| 色吧在线观看| 18禁黄网站禁片午夜丰满| 亚洲av中文字字幕乱码综合| 久久久久久大精品| 国产在线精品亚洲第一网站| 国产精品av视频在线免费观看| 男女午夜视频在线观看| 噜噜噜噜噜久久久久久91| 一级黄片播放器| 在线国产一区二区在线| 午夜福利在线在线| 在线观看美女被高潮喷水网站 | 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www| 欧美一区二区精品小视频在线| 欧美日韩综合久久久久久 | 日本熟妇午夜| 一a级毛片在线观看| 欧美在线一区亚洲| 欧美大码av| 国产亚洲精品久久久久久毛片| 午夜免费观看网址| 久久性视频一级片| 中文字幕av在线有码专区| 亚洲精品亚洲一区二区| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜一区二区| 亚洲欧美日韩高清专用| 国产色婷婷99| 亚洲美女黄片视频| 精品久久久久久久末码| 成人一区二区视频在线观看| 午夜影院日韩av| а√天堂www在线а√下载| 亚洲人成网站在线播| 又紧又爽又黄一区二区| 岛国视频午夜一区免费看| 精品一区二区三区视频在线 | 制服丝袜大香蕉在线| 特大巨黑吊av在线直播| 国产精品久久久久久久电影 | 91麻豆av在线| 精品人妻一区二区三区麻豆 | 亚洲国产精品久久男人天堂| a级一级毛片免费在线观看| 婷婷亚洲欧美| av中文乱码字幕在线| 黄片小视频在线播放| 18禁美女被吸乳视频| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看 | 真实男女啪啪啪动态图| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 法律面前人人平等表现在哪些方面| 在线观看免费视频日本深夜| 露出奶头的视频| 首页视频小说图片口味搜索| 动漫黄色视频在线观看| 欧美一区二区国产精品久久精品| 国产97色在线日韩免费| 精品福利观看| 琪琪午夜伦伦电影理论片6080| 欧美3d第一页| 手机成人av网站| 亚洲,欧美精品.| 亚洲成人久久爱视频| 国产成+人综合+亚洲专区| 少妇丰满av| www.www免费av| 免费观看精品视频网站| av天堂在线播放| 丝袜美腿在线中文| 成人精品一区二区免费| 国产蜜桃级精品一区二区三区| 日韩亚洲欧美综合| 91在线观看av| 国产精品免费一区二区三区在线| 亚洲欧美激情综合另类| 看黄色毛片网站| 在线天堂最新版资源| www日本黄色视频网| 综合色av麻豆| 精品欧美国产一区二区三| 欧美大码av| 亚洲av二区三区四区| 日韩欧美国产一区二区入口| 88av欧美| 久久久久性生活片| 黄色成人免费大全| 舔av片在线| 亚洲av不卡在线观看| 欧美三级亚洲精品| 51午夜福利影视在线观看| 国产高清videossex| 亚洲专区中文字幕在线| 日本 欧美在线| 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 亚洲最大成人手机在线| 中文资源天堂在线| 亚洲精品美女久久久久99蜜臀| 很黄的视频免费| 少妇高潮的动态图| a级毛片a级免费在线| 精品99又大又爽又粗少妇毛片 | 亚洲精品一卡2卡三卡4卡5卡| 天堂√8在线中文| 狂野欧美激情性xxxx| 熟女电影av网| 欧美性猛交╳xxx乱大交人| 午夜福利成人在线免费观看| 日韩精品青青久久久久久| 亚洲黑人精品在线| 亚洲 欧美 日韩 在线 免费| ponron亚洲| 日本免费a在线| www.熟女人妻精品国产| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 脱女人内裤的视频| av国产免费在线观看| 久久九九热精品免费| 无限看片的www在线观看| 午夜免费观看网址| 婷婷丁香在线五月| 中文字幕av成人在线电影| 日韩av在线大香蕉| 99久久精品一区二区三区| 国产精品香港三级国产av潘金莲| 国内精品久久久久久久电影| 成年女人看的毛片在线观看| 69av精品久久久久久| 国内精品久久久久精免费| 精品久久久久久,| 久久性视频一级片| 国产乱人伦免费视频| 久久香蕉国产精品| 最新美女视频免费是黄的| 欧美+亚洲+日韩+国产| 国产精品 国内视频| 国产欧美日韩精品一区二区| www.色视频.com| 亚洲人成伊人成综合网2020| 不卡一级毛片| 亚洲精品影视一区二区三区av| 一级毛片高清免费大全| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 精品人妻偷拍中文字幕| 成人无遮挡网站| 国产亚洲精品av在线| 宅男免费午夜| 亚洲国产中文字幕在线视频| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 99久久九九国产精品国产免费| 九色成人免费人妻av| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 色播亚洲综合网| 人人妻人人看人人澡| 国产私拍福利视频在线观看| 免费看日本二区| 免费在线观看影片大全网站| 俺也久久电影网| 国产免费男女视频| 欧美一区二区国产精品久久精品| 亚洲欧美日韩高清在线视频| 色综合亚洲欧美另类图片| 国产视频内射| 国产乱人视频| 精品欧美国产一区二区三| 婷婷亚洲欧美| 男女之事视频高清在线观看| 窝窝影院91人妻| 一区福利在线观看| 精品无人区乱码1区二区| 久久香蕉精品热| 男女之事视频高清在线观看| 在线播放无遮挡| 欧美黄色片欧美黄色片| svipshipincom国产片| 国产av在哪里看| 久久九九热精品免费| 激情在线观看视频在线高清| 综合色av麻豆| 美女被艹到高潮喷水动态| 草草在线视频免费看| 99热6这里只有精品| 成人特级黄色片久久久久久久| 性色av乱码一区二区三区2| 窝窝影院91人妻| 午夜免费观看网址| 色尼玛亚洲综合影院| 男人舔奶头视频| 男女之事视频高清在线观看| 亚洲av电影不卡..在线观看| 欧美绝顶高潮抽搐喷水| 国产真实伦视频高清在线观看 | 每晚都被弄得嗷嗷叫到高潮| 国产亚洲欧美在线一区二区| 九色国产91popny在线| 午夜福利欧美成人| 久久草成人影院| 精品电影一区二区在线| 欧美色视频一区免费| av福利片在线观看| 欧美成人性av电影在线观看| 欧美乱色亚洲激情| 久久久精品欧美日韩精品| 高潮久久久久久久久久久不卡| 久久婷婷人人爽人人干人人爱| 国产成人欧美在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲成人久久爱视频| 欧美三级亚洲精品| 日本一二三区视频观看| 国产老妇女一区| 欧美zozozo另类| 亚洲午夜理论影院| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 他把我摸到了高潮在线观看| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区 | 国产真实乱freesex| 99视频精品全部免费 在线| 18禁美女被吸乳视频| 日韩欧美在线乱码| 久久久国产精品麻豆| 怎么达到女性高潮| 亚洲av美国av| 久久久久精品国产欧美久久久| 18禁黄网站禁片免费观看直播| 在线免费观看不下载黄p国产 | 国产老妇女一区| 99国产精品一区二区蜜桃av| 91在线观看av| 欧美性感艳星| 色综合婷婷激情| 久久香蕉国产精品| a级一级毛片免费在线观看| 深夜精品福利| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 真实男女啪啪啪动态图| 岛国视频午夜一区免费看| 高清在线国产一区| 国产精品久久久人人做人人爽| 久久草成人影院| 国产三级在线视频| 国产视频内射| 亚洲成人精品中文字幕电影| 欧美中文日本在线观看视频| 少妇人妻一区二区三区视频| 又爽又黄无遮挡网站| 亚洲精品一区av在线观看| 亚洲成人久久爱视频| 两个人的视频大全免费| 欧美日韩一级在线毛片| 大型黄色视频在线免费观看| 亚洲成人久久爱视频| 国产高潮美女av| 久久久色成人| 人妻夜夜爽99麻豆av| 俺也久久电影网| 日本成人三级电影网站| 久久久久国内视频| 亚洲 欧美 日韩 在线 免费| 给我免费播放毛片高清在线观看| 亚洲精品在线美女| 亚洲国产色片| 国内精品美女久久久久久| 51午夜福利影视在线观看| 国产探花在线观看一区二区| 国产真人三级小视频在线观看| 亚洲最大成人手机在线| 动漫黄色视频在线观看| 1024手机看黄色片| 91久久精品电影网| 美女被艹到高潮喷水动态| 美女黄网站色视频| 99热这里只有是精品50| 亚洲精品色激情综合| 日韩欧美在线二视频| avwww免费| 国产成人欧美在线观看| 精品一区二区三区视频在线观看免费| 噜噜噜噜噜久久久久久91| 国产精品 欧美亚洲| av天堂在线播放| 久久久色成人| 国内少妇人妻偷人精品xxx网站| 好男人在线观看高清免费视频| 在线观看免费视频日本深夜| 禁无遮挡网站| 日韩欧美 国产精品| 精品熟女少妇八av免费久了| 欧美+亚洲+日韩+国产| 午夜激情欧美在线| 女人高潮潮喷娇喘18禁视频| 老司机午夜福利在线观看视频| 丁香欧美五月| 蜜桃久久精品国产亚洲av| 男女下面进入的视频免费午夜| 亚洲avbb在线观看| 欧美午夜高清在线| 国产色婷婷99| 90打野战视频偷拍视频| 欧美乱色亚洲激情| 色播亚洲综合网| 国产精品98久久久久久宅男小说| 日韩高清综合在线| 国产一区二区三区视频了| 99在线视频只有这里精品首页| 12—13女人毛片做爰片一| 久99久视频精品免费| 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| 人妻夜夜爽99麻豆av| 久久午夜亚洲精品久久| 两个人的视频大全免费| 在线观看免费午夜福利视频| 伊人久久大香线蕉亚洲五| 亚洲人成电影免费在线| 日韩精品青青久久久久久| 熟女人妻精品中文字幕| 神马国产精品三级电影在线观看| 色视频www国产| 最好的美女福利视频网| 欧美3d第一页| 国产不卡一卡二| 亚洲男人的天堂狠狠| 国产真人三级小视频在线观看| 精品午夜福利视频在线观看一区| 国产一区二区在线观看日韩 | 人人妻,人人澡人人爽秒播| 欧美色欧美亚洲另类二区| 最好的美女福利视频网| 亚洲欧美精品综合久久99| 美女大奶头视频| 亚洲欧美精品综合久久99| 搡老岳熟女国产| 中出人妻视频一区二区| 国产高清激情床上av| 色综合欧美亚洲国产小说| 男人舔女人下体高潮全视频| 国产高清videossex| 高清日韩中文字幕在线| 国产成人a区在线观看| 长腿黑丝高跟| 日韩欧美在线乱码| 日本免费a在线| 欧美绝顶高潮抽搐喷水| 99国产精品一区二区蜜桃av| 欧美日韩综合久久久久久 | 亚洲精品在线观看二区| 欧美黄色淫秽网站| 色噜噜av男人的天堂激情| 婷婷丁香在线五月| 国产主播在线观看一区二区| 最新美女视频免费是黄的| 男女床上黄色一级片免费看| 高清在线国产一区| 在线视频色国产色| 亚洲人成伊人成综合网2020| 国产午夜福利久久久久久| 亚洲av第一区精品v没综合| 少妇的逼水好多| 美女高潮喷水抽搐中文字幕| 亚洲第一电影网av| 搡老岳熟女国产| 老司机午夜十八禁免费视频| 国产精品久久久久久亚洲av鲁大| 色av中文字幕| 深夜精品福利| 在线观看舔阴道视频| 欧美成人性av电影在线观看| 夜夜看夜夜爽夜夜摸| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看 | 成年女人毛片免费观看观看9| 欧美bdsm另类| 亚洲av成人精品一区久久| 成人特级av手机在线观看| 伊人久久精品亚洲午夜| 午夜免费激情av| 大型黄色视频在线免费观看| 色综合欧美亚洲国产小说| or卡值多少钱| 搞女人的毛片| 免费搜索国产男女视频| 亚洲一区二区三区色噜噜| 非洲黑人性xxxx精品又粗又长| www日本黄色视频网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品爽爽va在线观看网站| 亚洲欧美日韩东京热| 国产成人影院久久av| 久久久久久九九精品二区国产| 欧美中文日本在线观看视频| 高清日韩中文字幕在线| 一区二区三区国产精品乱码| 久久6这里有精品| 九九在线视频观看精品| 操出白浆在线播放| av国产免费在线观看| 国产色爽女视频免费观看| 日韩欧美三级三区| 国产精品一区二区免费欧美| 美女高潮喷水抽搐中文字幕| 国产高清videossex| 国产免费av片在线观看野外av| 亚洲人成网站在线播| 亚洲片人在线观看| 色av中文字幕| 亚洲欧美日韩无卡精品| 亚洲成a人片在线一区二区| 国产三级中文精品| 亚洲人成电影免费在线| 国产精品综合久久久久久久免费| 久9热在线精品视频| 三级国产精品欧美在线观看| 毛片女人毛片| 国内精品美女久久久久久| 精品熟女少妇八av免费久了| 亚洲真实伦在线观看| 美女高潮的动态| 国产成人系列免费观看| 色综合欧美亚洲国产小说| 男女做爰动态图高潮gif福利片| 日本在线视频免费播放| 丰满人妻一区二区三区视频av | 91字幕亚洲| 在线免费观看不下载黄p国产 | 亚洲国产欧美网| 色综合站精品国产| 国产av麻豆久久久久久久| 琪琪午夜伦伦电影理论片6080| 国产欧美日韩精品一区二区| 国产一区二区亚洲精品在线观看| 麻豆一二三区av精品| or卡值多少钱| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 国产精品综合久久久久久久免费| 欧美另类亚洲清纯唯美| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩东京热| 日韩亚洲欧美综合| av国产免费在线观看| 夜夜躁狠狠躁天天躁| 一个人免费在线观看电影| 成年女人看的毛片在线观看| 欧美中文日本在线观看视频| 中文字幕人妻丝袜一区二区| 免费看光身美女| 国产精品嫩草影院av在线观看 | 欧美激情久久久久久爽电影| 久久久国产成人免费| 搡老岳熟女国产| 99热这里只有是精品50| 97碰自拍视频| 毛片女人毛片| 成人av一区二区三区在线看| 亚洲人与动物交配视频| 成人国产一区最新在线观看| 久久久久性生活片| 久久国产乱子伦精品免费另类| 成年免费大片在线观看| 亚洲美女黄片视频| 国产单亲对白刺激| 亚洲精品日韩av片在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美日韩高清专用| www日本黄色视频网| 国产精品电影一区二区三区| 免费搜索国产男女视频| 国产三级中文精品| 制服人妻中文乱码| 天美传媒精品一区二区| 12—13女人毛片做爰片一| 欧美极品一区二区三区四区| 变态另类丝袜制服| 婷婷精品国产亚洲av| 国内精品久久久久久久电影| 美女黄网站色视频| 天天一区二区日本电影三级| 在线观看免费午夜福利视频| 久久99热这里只有精品18| 国产亚洲精品av在线| 禁无遮挡网站| 中文字幕人成人乱码亚洲影| 国产精品野战在线观看| 欧美成人性av电影在线观看| 国产乱人视频| 亚洲天堂国产精品一区在线| 午夜亚洲福利在线播放| 12—13女人毛片做爰片一| 欧洲精品卡2卡3卡4卡5卡区| 国产精品自产拍在线观看55亚洲| 国产亚洲精品久久久久久毛片| 又爽又黄无遮挡网站| 少妇的逼好多水| 亚洲av五月六月丁香网| 免费在线观看成人毛片| 亚洲,欧美精品.| 成人高潮视频无遮挡免费网站| 国产成人aa在线观看| 国产欧美日韩一区二区三| 欧美日韩国产亚洲二区| 国产精品嫩草影院av在线观看 | 九九在线视频观看精品| 成年女人毛片免费观看观看9| 亚洲内射少妇av| av福利片在线观看| 麻豆成人午夜福利视频| 男女午夜视频在线观看| 亚洲国产精品合色在线| 一个人免费在线观看的高清视频| 麻豆一二三区av精品| 天天添夜夜摸| 国产精品乱码一区二三区的特点| 亚洲精华国产精华精| 亚洲一区二区三区色噜噜| 免费在线观看成人毛片| 亚洲片人在线观看| 怎么达到女性高潮| 欧美一区二区亚洲| 一级黄色大片毛片| 国产美女午夜福利| 五月伊人婷婷丁香| 两个人看的免费小视频| 乱人视频在线观看| 亚洲专区国产一区二区| 1024手机看黄色片| 搡女人真爽免费视频火全软件 | 亚洲欧美日韩无卡精品| 日日摸夜夜添夜夜添小说| 欧美xxxx黑人xx丫x性爽| 国产成人系列免费观看| 亚洲中文日韩欧美视频|