劉振賢 丁瑜 夏振堯 彭博識(shí)

摘要：合理描述根系和巖面粗糙度影響下土-巖接觸面的剪切行為是開展土巖混合邊坡穩(wěn)定評(píng)價(jià)的關(guān)鍵?；赪eibull分布并考慮含根率、粗糙度的影響建立了土-巖接觸面剪切統(tǒng)計(jì)損傷模型。根據(jù)試驗(yàn)結(jié)果，確定了含根率、粗糙度對(duì)接觸面剪切剛度的影響函數(shù)及模型參數(shù)，并驗(yàn)證了模型的可行性；進(jìn)一步，探討了含根率、粗糙度和模型參數(shù)m、w0對(duì)接觸面剪切變形和損傷演化的影響。研究表明：根系與粗糙度均能抑制接觸面損傷發(fā)展從而提高接觸面強(qiáng)度；提出的損傷模型能較好地描述根系和粗糙度影響下土-巖接觸面剪切損傷演化和剪切力學(xué)特性。研究成果可為土-巖接觸面剪切性能分析及土巖混合邊坡穩(wěn)定評(píng)價(jià)提供理論參考。

摘要：土-巖接觸面； 含根率； Weibull分布； 粗糙度； 剪切損傷模型； 生態(tài)防護(hù)； 土巖邊坡

中圖法分類號(hào)： TU44

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.01.030

0 引 言

巖土工程中普遍存在土與各類結(jié)構(gòu)物的相互作用，形成土-巖、墻-土、樁-土、筋-土等各類土與結(jié)構(gòu)物的接觸面問(wèn)題［1-4］。土與結(jié)構(gòu)物的接觸面剪切特性被認(rèn)為是控制相關(guān)工程整體安全與長(zhǎng)期穩(wěn)定的關(guān)鍵因素，受到學(xué)者和工程技術(shù)人員的密切關(guān)注。對(duì)于土與結(jié)構(gòu)物接觸面的剪切變形行為，國(guó)內(nèi)外學(xué)者圍繞土體級(jí)配與類型［5-6］、接觸面粗糙度［7］、環(huán)境條件（如侵蝕［8］、生物作用［9］）等因素開展了大量試驗(yàn)研究。

為了描述土與結(jié)構(gòu)物接觸面的剪切應(yīng)力-位移特性，眾多學(xué)者開展了本構(gòu)模型研究?；谠囼?yàn)，早期提出的模型主要有雙曲線剪切模型［10］、雙折線模型［11］、彈塑性模型［12］。20世紀(jì)90年代后，損傷理論逐步被引入接觸面剪切本構(gòu)研究［13］。胡黎明［14］、張嘎［15］等相繼對(duì)土與結(jié)構(gòu)物接觸面的損傷本構(gòu)進(jìn)行了理論分析和驗(yàn)證。近年來(lái)，在各類土與結(jié)構(gòu)物接觸面剪切變形研究中，統(tǒng)計(jì)損傷理論應(yīng)用日益廣泛，其在硫酸腐蝕作用下的樁-土接觸面剪切變形［16］、筋-土接觸面剪切軟化［17］、土與結(jié)構(gòu)物接觸面厚度對(duì)強(qiáng)度及變形影響［18］等研究中均具有良好的適用性。

在土巖混合邊坡與巖質(zhì)邊坡開展生態(tài)防護(hù)工程后，常有護(hù)坡基材脫落現(xiàn)象發(fā)生，脫落程度主要受土-巖接觸面力學(xué)特性影響。試驗(yàn)研究表明，土-巖接觸面剪切力學(xué)特性受土體性質(zhì)、巖面粗糙度、植物根系等影響。目前，對(duì)根系作用、粗糙度影響下的土-巖接觸面剪切特性以試驗(yàn)研究為主，相關(guān)理論模型研究較為缺乏。為此，本文基于統(tǒng)計(jì)損傷理論，考慮根系與粗糙度影響，建立了土-巖接觸面剪切統(tǒng)計(jì)損傷模型；根據(jù)原位剪切試驗(yàn)結(jié)果確定參數(shù)，對(duì)模型進(jìn)行驗(yàn)證，分析并探討含根率、粗糙度對(duì)土-巖接觸面剪切力學(xué)特性的影響。

1 土-巖接觸面原位剪切試驗(yàn)

采用混凝土試塊代替巖塊（尺寸30 cm×30 cm），以不同鋸齒狀鋼模在表面制造4種粗糙度，采用相對(duì)凸起高度RS表示，分別為3.85，4.12，4.84，5.00 mm。各試塊豎向預(yù)留9個(gè)孔隙用于種植灌木。在樣地鋪設(shè)種植土并種植多花木藍(lán)，以無(wú)植物為對(duì)照進(jìn)行試驗(yàn)（見圖1）。設(shè)置4種植株梯度，1 a后進(jìn)行原位剪切試驗(yàn)，按實(shí)際植株根系計(jì)算含根率，含根率以RAR表示。用于試驗(yàn)的原位直剪儀由框架、直剪盒、千斤頂、液壓泵、位移傳感器、壓力傳感器、顯示器等組成。剪切盒尺寸為30 cm×30 cm，最大剪切位移為100 mm。試驗(yàn)測(cè)定了不同粗糙度、含根率下界面剪切應(yīng)力-位移之間的關(guān)系，具體試驗(yàn)過(guò)程與結(jié)果詳見文獻(xiàn)［19］。

2 土-巖接觸面剪切損傷模型

2.1 基本假定

為建立土-巖接觸面的統(tǒng)計(jì)損傷模型，本文做出如下假定：土-巖接觸面微元強(qiáng)度服從Weibull分布；土-巖接觸面由理想原狀微元和損傷微元組成，接觸面荷載由損傷和未損傷部分共同承擔(dān)；接觸面微元強(qiáng)度服從Weibull分布函數(shù)；無(wú)損原狀微元具有線彈性性質(zhì)，所受剪切應(yīng)力遵循線彈性法則；損傷部分所承受荷載等于接觸面剪切殘余應(yīng)力。

2.2 考慮根系與粗糙度影響的統(tǒng)計(jì)損傷模型

接觸面微元強(qiáng)度F服從Weibull分布時(shí)，其概率密度函數(shù)表示為

f（F）=mF0FF0m-1exp［-FF0m］（1）

式中：m、F0為Weibull分布的分布參數(shù)。

土-巖接觸面由N個(gè)微元組成，剪切加載至某應(yīng)力水平［F，F(xiàn)+dF］時(shí)，損傷微元數(shù)目為Nd（F）。根據(jù)統(tǒng)計(jì)損傷理論［20-21］，損傷微元數(shù)目可由失效概率確定，即

Nd（F）=N∫F0f（x）dx（2）

相應(yīng)地，定義土-巖接觸面損傷變量D為剪切作用下?lián)p傷微元數(shù)量Nd（F）與微元總數(shù)N的比值，即

D=Nd（F）N=N∫F0f（x）dxN（3）

由式（1）和式（3），可得土-巖接觸面損傷變量：

D=1-exp-FF0m（4）

式中：無(wú)損原狀接觸面微元強(qiáng)度F為理想線彈性介質(zhì)。根據(jù)相關(guān)研究［20-21］，F(xiàn)可由接觸面的剪切剛度k和剪切位移w來(lái)表示：

F=kw（5）

研究認(rèn)為，式（4）中參數(shù)F0代表微元體平均強(qiáng)度。參考Xie等［22］研究，將F0表示為

F0=kw0（6）

式中：w0為剪切作用下的特征位移。

將式（5）、（6）代入式（4）中，損傷變量D可進(jìn)一步表示為

D=1-exp-ww0m（7）

由于剪切荷載由土-巖接觸面的損傷及未損傷部分共同承擔(dān)，可記為

τ=kw（1-D）+τrD（8）

式中：τr為接觸面的殘余剪切應(yīng)力。

由式（7）、（8），得到土-巖接觸面的統(tǒng)計(jì)損傷剪切模型為

τ=（kw-τr）exp-ww0m+τr（9）

根系與粗糙度均對(duì)土-巖接觸面的剪切變形具有重要影響［19］，為在式（9）基礎(chǔ)上進(jìn)一步考慮二者對(duì)其作用，本文采用剪切剛度表征二者的影響。對(duì)于不含根系的光滑接觸面，其剪切剛度為k0，而對(duì)于含根系的粗糙接觸面，其剪切剛度k可通過(guò)根系和粗糙度的影響函數(shù)表示，即

k=k0·f（RAR）·g（RS）（10）

式中：f（RAR）為含根率對(duì)剪切剛度的影響函數(shù)；g（RS）為粗糙度對(duì)剪切剛度的影響函數(shù)。

將式（10）代入式（9）可得到受根系與粗糙度共同影響的土-巖接觸面統(tǒng)計(jì)損傷剪切模型如下：

τ=k0·f（RAR）·g（RS）w-τrexp-ww0m+τr（11）

3 參數(shù)確定

3.1 損傷參數(shù)m、w0及影響函數(shù)確定

通過(guò)原位剪切試驗(yàn)獲取的位移-應(yīng)力曲線，首先確定式（11）中的k0、τr。參數(shù)m、w0及影響函數(shù)f（RAR）、g（RS）可通過(guò)計(jì)算進(jìn)一步確定。根據(jù)峰值點(diǎn)處應(yīng)力對(duì)位移的一階導(dǎo)數(shù)為零且峰值點(diǎn)滿足損傷模型式（11），可確定m、w0，即

m=k0·f（RAR）·g（RS）wpk0·f（RAR）·g（RS）wp-τrlnk0·f（RAR）·g（RS）wp-τrτp-τr（12）

w0=wplnk0·f（RAR）·g（RS）wp-τrτp-τr1m（13）

式中：τp、wp可由試驗(yàn)曲線峰值點(diǎn)獲取。

對(duì)含根率和粗糙度影響函數(shù)，采用單因素分析法確定。對(duì)于f（RAR），根據(jù)試驗(yàn)結(jié)果，在每一粗糙度下，將不同含根率（0.51‰，0.97‰，1.30‰，1.99‰）時(shí)的剛度與含根率為零時(shí)做歸一化處理，如表1所列。分析發(fā)現(xiàn)各組粗糙度下，歸一化剛度值與含根率線性正相關(guān)，得到含根率影響函數(shù)：f（RAR）=0.093 5（RAR）+0.994（R2=0.996）。

獲取f（RAR）之后，則可確定粗糙度對(duì)剛度的影響函數(shù)g（RS）。由式（10）可得：

g（RS）=kk0·f（RAR）（14）

由試驗(yàn)測(cè)知k0=0.995，則根據(jù)式（14）可求解不同粗糙度情況下的g（RS）值如表2所列。經(jīng)分析，剛度與粗糙度呈二次函數(shù)關(guān)系，經(jīng)擬合得粗糙度影響函數(shù)：g（RS）=-0.29（RS）2+2.46（RS）-4.44（R2=0.99）。

3.2 參數(shù)m、w0的計(jì)算結(jié)果

由式（12）、（13）計(jì)算得到各組粗糙度下不同含根率對(duì)應(yīng)的模型參數(shù)，如表3所列。

4 模型驗(yàn)證及分析

4.1 RAR=0時(shí)的粗糙土-巖接觸面

含根率為0時(shí)，模型計(jì)算得到的各組粗糙度下的接觸面剪切應(yīng)力-位移關(guān)系曲線如圖2所示。試驗(yàn)結(jié)果表明［19］，無(wú)根系作用時(shí)，土-巖接觸面宏觀剪切變形經(jīng)歷線彈性、彈塑性、剪切軟化、殘余變形4個(gè)階段，剪切應(yīng)力-位移曲線具有典型的剪切軟化特性。在剪切至一定加載水平時(shí)，土-巖接觸面微元發(fā)生破損；隨著剪切變形發(fā)展，損傷會(huì)逐步累積演化。

從唯象角度分析，隨著剪切變形不斷增長(zhǎng)，接觸面剪切剛度會(huì)不斷弱化。由圖2不難發(fā)現(xiàn)，不同粗糙度下，試驗(yàn)曲線與模型曲線吻合較好。這表明無(wú)根系作用時(shí)，通過(guò)考慮粗糙度對(duì)剛度的影響，本文的損傷模型能較好地描述接觸面剪切應(yīng)力-位移發(fā)展演化過(guò)程。

4.2 不同含根率與粗糙度的土-巖接觸面

在4種粗糙度下，考慮不同含根率，模型計(jì)算和試驗(yàn)得到的剪應(yīng)力-位移曲線如圖3所示。

根系通過(guò)部分分擔(dān)荷載發(fā)揮錨固作用，從而顯著提升接觸面的抗剪性能［21］。根據(jù)試驗(yàn)結(jié)果，對(duì)于同一粗糙度接觸面，根系增強(qiáng)作用表現(xiàn)為接觸面的峰值應(yīng)力、殘余應(yīng)力均隨含根率增加而增大。從損傷的角度而言，根系因深入裂隙而對(duì)周圍土體具有一定錨固作用，通過(guò)自身變形或拉伸破壞抵抗一部分剪切力（見圖4），相對(duì)于無(wú)根系的接觸面，根系可明顯降低接觸面損傷，在增強(qiáng)接觸面強(qiáng)度的同時(shí)延遲接觸面的剪切破壞；粗糙度的提高可增加接觸面摩擦力，增加土-巖實(shí)際接觸面積，提升咬合力，兩者均能抑制接觸面微元的破損。

對(duì)于根系作用下的粗糙接觸面，本文提出的同時(shí)考慮含根率和粗糙度影響的損傷模型可以較好描述其剪切變形過(guò)程以及強(qiáng)度特征。由圖3不難發(fā)現(xiàn)，模型在彈性與彈塑性變形階段的計(jì)算結(jié)果與試驗(yàn)結(jié)果高度吻合，兩者之間誤差在9.6%以內(nèi)，且主要差異存在于軟化階段。對(duì)比發(fā)現(xiàn)，軟化階段接觸面實(shí)際承載能力稍高于模型計(jì)算結(jié)果。造成差異的原因可能在于，根系作用下接觸面破損時(shí)仍有較好的承載能力，含根率越高時(shí)越明顯。

由此可見，本文提出的損傷模型能較好地描述含根率、粗糙度對(duì)土-巖接觸面剪切變形的影響，可為分析土-巖接觸面剪切性能、評(píng)價(jià)土巖混合邊坡穩(wěn)定性提供一定的理論指導(dǎo)。

5 討 論

5.1 含根率與粗糙度對(duì)損傷的影響

試驗(yàn)結(jié)果和模型計(jì)算顯示，含根率、粗糙度均會(huì)影響土-巖接觸面的破損和損傷演化，進(jìn)而改變接觸面剪切變形過(guò)程和抗剪性能。含根率與粗糙度的變化直接影響模型參數(shù)w0與m，可由式（7）分析含根率與粗糙度對(duì)損傷的影響。

某一粗糙度下，以RS=4.12 mm為例（見圖5），不同含根率下?lián)p傷變量D隨位移呈負(fù)指數(shù)演化，其演化曲線呈“S”形。在演化初期，含根率越大時(shí)損傷越低，但損傷總體較微弱，不同含根率之間損傷差異較小。在損傷累積發(fā)展階段，無(wú)根系作用時(shí)，損傷發(fā)展最快；含根率越大，損傷延滯越明顯。相應(yīng)地，隨著含根率增大，達(dá)到完全損傷（損傷累積趨近于1）的位移增大，從含根率為0時(shí)的約40 mm變?yōu)楹?.98‰時(shí)的60 mm左右。

分析表明，根系在剪切變形初期尚未發(fā)揮作用，但土-巖接觸面微元發(fā)生一定數(shù)量破損后，損傷不斷累積發(fā)展，根系對(duì)接觸面損傷的抑制作用逐步突顯。

含根率為零時(shí)，4種粗糙度對(duì)損傷演化的影響如圖6所示。可以看出，粗糙度對(duì)初期損傷、損傷發(fā)展、最終損傷各階段的影響總體與根系作用類似。

4種粗糙度下的損傷演化曲線差異明顯，粗糙度增大可顯著延緩損傷發(fā)展速度。分析認(rèn)為，粗糙度的增大能有效抑制損傷、增強(qiáng)接觸面強(qiáng)度，其根本原因在于粗糙度增大增加了土巖之間的接觸面積。以相對(duì)凸起高度表征粗糙度增大，意味著土巖接觸面積明顯增大。因此，相同大小的試樣粗糙度越大，土巖相互作用面積越大，土巖嚙合作用增強(qiáng)，更能抑制損傷發(fā)展，提高接觸面強(qiáng)度。

5.2 參數(shù)w0的影響

模型參數(shù)w0為表征接觸面微元平均強(qiáng)度的特征位移。為探究其對(duì)接觸面損傷和剪切變形的影響，在表3給出的m值（2.77），取不同w0值進(jìn)行分析。

由圖7可以看出，相同m時(shí)，隨著w0增大，損傷演化速率減緩，達(dá)到完全損傷時(shí)的剪切位移增大；此外，w0增大不僅引起峰值剪應(yīng)力增大，對(duì)應(yīng)的峰值剪切位移逐漸向后推移。由此，可推測(cè)參數(shù)w0實(shí)質(zhì)上反映了接觸面峰值剪應(yīng)力與峰值位移的平均水平。

5.3 參數(shù)m的影響

w0一定時(shí)（按表3取w0為30.37 mm），分析不同m值對(duì)接觸面損傷和剪切性能的影響。如圖8（a）所示，相同w0時(shí)，損傷變量D隨m增大在某一剪切位移時(shí)相交于A點(diǎn)（23.97，0.64）：A點(diǎn)之前，m越小，D越大；A點(diǎn)之后，完全損傷之前，m越小，D越小，達(dá)到完全損傷對(duì)應(yīng)的剪切位移逐漸增大。

已有研究表明［23］，m是反映土-巖接觸面微元強(qiáng)度分布的集中程度和表征接觸面脆性程度的參數(shù)。相應(yīng)地，當(dāng)w0為定值時(shí)，m越大，峰值剪應(yīng)力越大，于圖8（b）中B點(diǎn)（23.97，16.15）后可明顯看出，峰后曲線越陡，土-巖接觸面剪切損傷軟化越嚴(yán)重，微元強(qiáng)度分布越集中，接觸面的脆性越明顯。分析表明，參數(shù)m是影響接觸面損傷演化和應(yīng)力-位移曲線形態(tài)的敏感指標(biāo)，模型計(jì)算時(shí)應(yīng)綜合考慮含根率、粗糙度對(duì)接觸面強(qiáng)度的影響從而合理確定。

6 結(jié) 論

本文通過(guò)考慮含根率、粗糙度的影響，分析和探討了土-巖接觸面剪切損傷變形特征，得到以下結(jié)論：

（1） 本文提出的損傷本構(gòu)模型的計(jì)算結(jié)果與試驗(yàn)吻合較好，兩者之間最大誤差小于9.6%，表明該模型可用于描述無(wú)根系和有根系的粗糙土-巖接觸面的剪切變形過(guò)程和剪切力學(xué)性能。

（2） 根系含量、粗糙度的增大會(huì)不同程度地抑制損傷發(fā)展累積。通過(guò)確定具體土-巖接觸面含根率、相對(duì)凸起高度，可分析和揭示根系作用、巖面粗糙度對(duì)接觸面損傷演變過(guò)程的影響，結(jié)果表明粗糙度較含根率對(duì)損傷的影響更為明顯。

（3） 模型參數(shù)m主要影響剪切過(guò)程中的軟化階段，m越大，接觸面脆性越大，因此軟化階段的損傷速率越大；參數(shù)w0主要影響峰值應(yīng)力對(duì)應(yīng)的剪切位移，w0越大，損傷速率越小，峰值應(yīng)力對(duì)應(yīng)的剪切位移越大。

參考文獻(xiàn)：

［1］ 張朝暉，婁宗科.土巖接觸復(fù)雜地層中接觸面力學(xué)性質(zhì)研究［J］.人民長(zhǎng)江，2010，41（17）：16-18，36.

［2］ MUNOZ J J，LYAMIN A V，HUERTA A.Stability of anchored sheet wall in cohesive-frictional soils by FE limit analysis［J］.International Journal for Numerical and Analytical Methods in Geomechanics，2013，37（9）：1213-1230.

［3］ WANG Y，SANG S，ZHANG M，et al.Field test of earth pressure at pile-soil interface by single pile penetration in silty soil and silty clay［J］.Soil Dynamics and Earthquake Engineering，2021，145（1）：106666.

［4］ OURIA A，KARAMZADEGAN S，EMAMI S.Interface properties of a cement coated geocomposite［J］.Construction and Building Materials，2021，266（1）：121014.

［5］ WANG H L，ZHOU W H，YIN Z Y，et al.Effect of grain size distribution of sandy soil on shearing behaviors at soil-structure interface［J］.Journal of Materials in Civil Engineering，2019，31（10）：04019238.

［6］ SHEN W，ZHAO T，DAI F.Influence of particle size on the buffering efficiency of soil cushion layer against rockfall impact［J］.Natural Hazards，2021，108（2）：1469-1488.

［7］ 陳美婷，趙光思，雪青華，等.粗糙度對(duì)接觸面處土石混合體剪切特性的影響［J］.地下空間與工程學(xué)報(bào)，2020，16（5）：1328-1337.

［8］ CHEN C，MEI S，CHEN S，et al.Laboratory investigation of erosion behavior at the soil-structure interface affected by various structural factors［J］.Natural Hazards，2021，111（1）：1-20.

［9］ BAK H M，KARIMINIA T，SHAHBODAGH B，et al.Application of bio-cementation to enhance shear strength parameters of soil-steel interface［J］.Construction and Building Materials，2021，294（5）：123470.

［10］ CLOUGH G W，DUNCAN J M.Finite element analyses of retaining wall behavior［J］.Journal of Soil Mechanics & Foundation Engineering，1971，97（12）：1657-1673.

［11］ BRANDT J R T.Behavior of soil-concrete interface［D］.Edmonton：University of Alberta，1985.

［12］ 陳慧遠(yuǎn).摩擦接觸單元及其分析方法［J］.水利學(xué)報(bào)，1985（4）：44-50.

［13］ DESAI C S，MA Y.Modeling of joints and interfaces using the disturbed state concept［J］.International Journal for Numerical and Analytical Methods in Geomechanics，1992，16（9）：623-653.

［14］ 胡黎明，濮家騮.土與結(jié)構(gòu)物接觸面損傷本構(gòu)模型［J］.巖土力學(xué)，2002，23（1）：6-11.

［15］ 張嘎，張建民.粗粒土與結(jié)構(gòu)接觸面統(tǒng)一本構(gòu)模型及試驗(yàn)驗(yàn)證［J］.巖土工程學(xué)報(bào)，2005，27（10）：1175-1179.

［16］ CAO X，LI J.A shear damage model of the interface between soil and sulfate-corroded concrete［J］.Construction and Building Materials，2021，293（6）：123450.

［17］ 成浩，王晅，張家生，等.加筋粗粒土筋土界面剪切特性與統(tǒng)計(jì)損傷軟化模型研究［J］.鐵道科學(xué)與工程學(xué)報(bào)，2018，15（11）：2780-2787.

［18］ 楊林德，劉齊建.土-結(jié)構(gòu)物接觸面統(tǒng)計(jì)損傷本構(gòu)模型［J］.地下空間與工程學(xué)報(bào)，2006（1）：79-86.

［19］ 丁瑜，夏振堯，許文年，等.根系作用下的基材土-巖接觸面原位剪切試驗(yàn)研究［J］.巖土工程學(xué)報(bào)，2016，38（11）：2107-2113.

［20］ 魏檣.含砂高液限黏土統(tǒng)計(jì)損傷硬化模型研究［J］.人民長(zhǎng)江，2015，46（10）：52-55.

［21］ 石泉彬，楊平.凍結(jié)粉細(xì)砂與鋼板接觸面剪切統(tǒng)計(jì)損傷模型構(gòu)建［J］.鐵道科學(xué)與工程學(xué)報(bào)，2021，18（10）：2591-2599.

［22］ XIE S，LIN H，WAN Y，et al.A statistical damage constitutive model considering whole joint shear deformation［J］.International Journal of Damage Mechanics，2020，29（6）：988-1008.

［23］ 丁瑜，魏偉兵，潘波，等.纖維加筋植被混凝土基材統(tǒng)計(jì)損傷模型研究［J］.巖土工程學(xué)報(bào)，2022，44（4）：652-659.

（編輯：鄭 毅）

Research on shear damage model of rough soil-rock interface influenced by root

LIU Zhenxian1，2，DING Yu1，2，XIA Zhenyao1，2，PENG Boshi1，2

（1.Key Laboratory of Geological Hazards on Three Gorges Reservoir Area of Ministry of Education（China Three Gorges University），Yichang 443002，China； 2.Hubei Key Laboratory of Disaster Prevention and Mitigation（China Three Gorges University），Yichang 443002，China）

Abstract：

To evaluate the safety and stability of soil-rock mixed slopes，it is vital to reasonably describe the shear behavior of soil-rock interface influenced by roots and roughness.In this study，based on Weibull distribution，a statistical shear damage model of soil-rock interface was established by taking the influences of root content and roughness into account.According to the experimental results，the influence functions of root content and roughness on the shear stiffness of the contact interface and model parameters were determined，and the model feasibility was verified.Then，the influences of root content，roughness，and model parameters（m and w0），on the shear deformation of contact interface and damage evolution were further discussed.The results showed that root and roughness could inhibit the damage development of interface and improve the strength of interface.The damage model proposed in this study could describe the shear damage evolution and shear mechanics characteristics of soil-rock contact surface under the influence of roots and roughness，and provide a theoretical basis for accurate analysis and evaluation of shear behavior of soil-rock contact surface and stability of soil-rock mixed slope.

Key words：

soil-rock interface；root content；Weibull distribution；roughness；shear damage model；ecological slope protection；soil-rock slope