朱德勝, 朱佳能, 李志剛, 高 越
(揚(yáng)州大學(xué)建筑科學(xué)與工程學(xué)院, 江蘇 揚(yáng)州 225127)
邊坡穩(wěn)定問題對(duì)城市地下基坑開挖、水利水電、公路路基等工程建設(shè)具有十分重要的意義.眾多專家學(xué)者提出了各種邊坡穩(wěn)定性的評(píng)價(jià)方法,并采用安全系數(shù)作為評(píng)價(jià)邊坡穩(wěn)定性的指標(biāo);然而,由于邊坡工程具有不確定性的特點(diǎn), 采用單一的安全系數(shù)不能對(duì)邊坡的可靠度進(jìn)行合理評(píng)價(jià).邊坡工程的可靠度研究最早出現(xiàn)在20世紀(jì)70年代(如Alonso[1]), 目前邊坡穩(wěn)定的可靠度分析方法主要有: 可靠度估算法[2-3]、響應(yīng)面法[4-5]、蒙特卡羅法[6-7].在蒙特卡羅法的框架下, Griffiths等[8]提出的邊坡可靠度隨機(jī)有限元方法, 結(jié)合了強(qiáng)度折減有限元分析和隨機(jī)場理論,在邊坡穩(wěn)定可靠度研究中被廣泛應(yīng)用[9-12].現(xiàn)有的邊坡可靠度隨機(jī)有限元分析主要還是基于各向同性隨機(jī)場模型,很少考慮土體參數(shù)的各向異性.由于土體在水平方向和豎直方向上經(jīng)歷不同的地質(zhì)作用,使得土體參數(shù)往往表現(xiàn)出各向異性[13].本文基于各向異性隨機(jī)場,對(duì)黏土邊坡的可靠度進(jìn)行分析,探究各向異性對(duì)黏土邊坡可靠度及破壞機(jī)制的影響規(guī)律.
圖1 黏土邊坡幾何模型Fig.1 Geometry model of a clay slope
圖1為黏土邊坡幾何模型.圖中β為邊坡坡度,H為邊坡高度, 內(nèi)摩擦角φu=0°,γ為土體重度.本文取邊坡高度H=10 m, 土體重度γ=20 kN·m-3.
本文采用隨機(jī)有限單元法計(jì)算給定參數(shù)下的黏土邊坡可靠度.在蒙特卡羅模擬的框架下,采用該方法模擬2 000次,每次模擬首先生成不排水強(qiáng)度的各向異性隨機(jī)場, 隨后對(duì)黏土邊坡加重力荷載.如果隨機(jī)有限元程序500次以內(nèi)不收斂,則定義邊坡失穩(wěn)破壞, 2 000次模擬中邊坡失穩(wěn)破壞次數(shù)的比例定義為黏土邊坡的失效概率pf.
由于土體經(jīng)歷復(fù)雜的地質(zhì)演化過程,其參數(shù)往往具有空間變異性,可采用隨機(jī)場理論對(duì)土體參數(shù)空間變異性進(jìn)行描述.黏土邊坡穩(wěn)定的主要抗力是土體不排水強(qiáng)度參數(shù),其空間變異性表現(xiàn)為不排水強(qiáng)度取值隨空間位置而改變,即土體空間內(nèi)兩點(diǎn)之間的不排水強(qiáng)度值有所不同;同時(shí),黏土邊坡區(qū)域內(nèi)的不排水強(qiáng)度也具有相關(guān)性,距離越近,不排水強(qiáng)度取值越接近.相關(guān)距離θ是衡量土體不排水強(qiáng)度空間變異性的重要特征: 相關(guān)距離越小,土體不排水強(qiáng)度的空間變異性越大.在相關(guān)距離之內(nèi),黏土邊坡的土體不排水強(qiáng)度之間相關(guān)性較強(qiáng);而在相關(guān)距離之外,黏土邊坡的土體不排水強(qiáng)度之間則相關(guān)性較弱.此外,采用方差折減系數(shù)描述土體不排水強(qiáng)度的“點(diǎn)”變異性和“空間”變異性之間的關(guān)系,即土體不排水強(qiáng)度在某一區(qū)域內(nèi)的“空間”平均值的方差等于土體不排水強(qiáng)度“點(diǎn)”的方差與方差折減系數(shù)的乘積.
為了避免各向異性隨機(jī)場模擬生成土體不排水強(qiáng)度為負(fù)值,本文假設(shè)黏土邊坡的不排水強(qiáng)度cu服從對(duì)數(shù)正態(tài)分布, 其變異系數(shù)vcu通常為0.1~0.5[13].考慮各向異性的影響, 本文參數(shù)分析中取θy=2 m,ξ為1~10.
圖2 不同安全系數(shù)均值下邊坡失效 概率與各向異性系數(shù)的關(guān)系Fig.2 Slope failure versus degree of anisotropy with different mean factors of safety
圖3 不同vcu下邊坡失效概率與各向異性系數(shù)的關(guān)系Fig.3 Slope failure versus degree of anisotropy with different coefficients of variation
圖4 不同坡度下邊坡失效概率與各向異性系數(shù)的關(guān)系Fig.4 Slope failure versus degree of anisotropy with different slope angles
圖5 各向同性隨機(jī)場下算例邊坡破壞機(jī)制Fig.5 Failure mechanism for the example slope in isotropic random field
圖6 各向異性隨機(jī)場下算例邊坡破壞機(jī)制Fig.6 Failure mechanism for the example slope in anisotropic random field
1) 當(dāng)黏土邊坡安全系數(shù)均值大于等于1.3時(shí), 忽略各向異性的影響會(huì)低估實(shí)際邊坡的失效概率,得到非保守的結(jié)果.黏土邊坡坡度越緩,越不能忽略各向異性的影響.
2) 考慮各向異性時(shí), 黏土邊坡的破壞機(jī)制會(huì)表現(xiàn)出部分的“線性”形狀.