吳順川, 賈文松, 韓龍強(qiáng), 劉興雷
(1.北京科技大學(xué)金屬礦山高效開采與安全教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083; 2.昆明理工大學(xué)國土資源工程學(xué)院,云南 昆明 650093; 3.自然資源部高原山地地質(zhì)災(zāi)害預(yù)報(bào)預(yù)警與生態(tài)保護(hù)修復(fù)重點(diǎn)實(shí)驗(yàn)室,云南 昆明 650093)
擋土墻土壓力理論計(jì)算普遍采用以極限平衡為理論基礎(chǔ)的朗肯土壓力和庫侖土壓力計(jì)算方法,但在擋土墻后為有限成坡土體條件下,這2 種經(jīng)典理論計(jì)算方法均不適用。
目前大多采用模型試驗(yàn)或理論推導(dǎo)方法研究擋土墻有限土體土壓力問題[1-4],但利用模擬手段建立三維模型分析露天礦工程中某些因素對土壓力分布規(guī)律影響的研究較少。 本文以河北騰龍露天礦止水固坡工程為研究對象,在擋土墻后有限成坡土體簡化模型的基礎(chǔ)上,利用FLAC3D軟件建立擋土墻和有限成坡土體的三維數(shù)字模型,通過有限成坡土體滑裂面傾角變化規(guī)律、有限成坡土體被動土壓力隨墻后土體深度的變化規(guī)律來驗(yàn)證該三維數(shù)字模型的有效性,然后依據(jù)該模型模擬分析平動條件下?lián)跬翂ν馏w參數(shù)、邊坡坡角、平臺寬度等參數(shù)對土壓力的影響。
騰龍露天鐵礦位于河北省遷安市西南部,地處平原,因北接燕山山脈,整體地勢南低北高。 礦區(qū)內(nèi)工程地質(zhì)情況較好,巖石種類繁多,如片麻巖、石英巖等,但巖體風(fēng)化嚴(yán)重,呈破碎狀態(tài),最上層由第四系地層覆蓋。 第四系地層主要為粗砂、卵石、砂礫等顆粒,透水性較好,大氣降水通過第四系地層補(bǔ)給地下水。 采場附近有灤河和沙河兩條河流,受河流以及大氣降水的影響,該采場地下水豐富,嚴(yán)重影響了采場開采工作以及靠幫邊坡穩(wěn)定性,容易導(dǎo)致滑坡,甚至無法形成邊坡[5]。 為了解決地下水對邊坡穩(wěn)定性的重大影響,必須采取合理的堵排水措施。 文獻(xiàn)[6]提出利用大型單一結(jié)構(gòu)地下連續(xù)墻進(jìn)行固坡止水,該措施成功解決了騰龍露天礦邊坡高水壓-低強(qiáng)度復(fù)雜條件下的止水固坡技術(shù)難題。 單一結(jié)構(gòu)地下連續(xù)墻布設(shè)如圖1 所示。本文在該措施的基礎(chǔ)上簡化得到有限成坡土體模型(見圖2)進(jìn)行土壓力分布規(guī)律的模擬研究。
圖1 單一結(jié)構(gòu)地下連續(xù)墻結(jié)構(gòu)布設(shè)圖
圖2 有限成坡土體簡化模型
模型由擋土墻和有限成坡土體兩部分組成,二者以六面塊體的方式建立,通過兩個(gè)六面塊體之間的相互作用來模擬墻土之間的實(shí)際關(guān)系。 擋土墻為矩形,其中墻長設(shè)為b,墻高設(shè)為h。 墻后土體為有限成坡土體,土體截面為梯形,為了減小邊界效應(yīng)對模擬計(jì)算結(jié)果的影響,模型中墻后有限成坡土體Y方向長度設(shè)為15b,Z方向高度設(shè)為3h,X方向上部平臺寬度為a、下部土體寬度為a+Hcotβ,其中β為墻后土體的坡角。建模時(shí),剛性擋土墻的一面與填土豎直面相接觸,坐標(biāo)原點(diǎn)位于墻土接觸面頂部,且模型整體關(guān)于X-Z面對稱,擋土墻平動方向與X軸方向一致。 三維實(shí)體模型如圖3 所示。 為了保證模型更加符合實(shí)際工程,依據(jù)騰龍鐵礦南側(cè)一剖面確定擋土墻和土體尺寸,并在模擬計(jì)算過程中保持不變,最終確定擋土墻墻長b=3 m,墻高h(yuǎn)=12 m,墻后有限成坡土體Y方向長度為15b=45 m,Z方向高度H=36 m。
圖3 三維實(shí)體模型示意圖
圖4 為三維數(shù)字模型的網(wǎng)格劃分情況,以六面體單元對模型進(jìn)行網(wǎng)格劃分,為了保證墻土之間更好地相互作用,墻土接觸面上的節(jié)點(diǎn)需要重合,墻后土體為模擬觀察部分,故該部分劃分的網(wǎng)格間距更小,以便更好地反映土體的變形。墻土間的接觸面假定為剛性的,在擋土墻平動過程中可以更好地將墻土間的摩擦表現(xiàn)出來。 擋土墻也設(shè)置為絕對剛性,保證擋土墻在平動過程中仍然保持穩(wěn)定狀態(tài),土體整體保持均勻狀態(tài)并且各個(gè)方向性質(zhì)一致,破壞過程遵從摩爾-庫倫屈服準(zhǔn)則[7]。
圖4 三維數(shù)字模型網(wǎng)格劃分圖
模型土體上表面邊界為自由狀態(tài),不受約束。 模型X=0 邊界有限元網(wǎng)格節(jié)點(diǎn)處約束其X方向位移,允許其豎向沉降;在Y=±7.5b邊界有限元網(wǎng)格節(jié)點(diǎn)處約束其Y方向位移;模型底部Z=-3h邊界有限元網(wǎng)格節(jié)點(diǎn)處約束其Z方向位移,允許其在X、Y方向有變形。
影響因素和參數(shù)取值如表1 所示。 土體本構(gòu)模型采用FLAC3D中的Mohr-Coulomb 模型,針對所選剖面的實(shí)際工程地質(zhì)情況,為了使模擬結(jié)果更加符合實(shí)際,確定擋土墻后有限土體為無黏性土,即c=0 kPa,有限土體其他參數(shù)列于表2 中。 將有限成坡土體和擋土墻之間的接觸以接觸面的形式呈現(xiàn),該接觸面參數(shù)如表3所示。 表2 和表3 參數(shù)取值均依據(jù)騰龍露天礦實(shí)際勘察資料確定。
表1 影響因素及參數(shù)取值一覽
表2 有限土體參數(shù)
表3 接觸面參數(shù)
在自重以及外界初始條件下,經(jīng)過一定步驟計(jì)算,模型達(dá)到初始平衡。 達(dá)到平衡后,開始對平動條件下的擋土墻和有限成坡土體進(jìn)行整體計(jì)算,得到墻背上的土壓力。 擋土墻沿X方向勻速移動,速度施加在擋土墻有限元網(wǎng)格節(jié)點(diǎn)上,模型計(jì)算再次達(dá)到平衡,墻后土體達(dá)到臨界狀態(tài),此時(shí)作用在墻背上的土壓力為極限值,即為進(jìn)行土壓力分布規(guī)律分析的參考值。 本文以下分析所需要的參考值均沿Z軸負(fù)方向獲取,即沿墻后土體深度方向,并且土體深度用h0表示。
為了保證模擬分析結(jié)果的準(zhǔn)確性,在進(jìn)行墻背上土壓力分布規(guī)律影響分析之前需要對有限元模型的有效性進(jìn)行驗(yàn)證。 下面通過有限成坡土體滑裂面傾角變化規(guī)律以及有限成坡土體被動土壓力解析解和數(shù)值解隨深度的變化規(guī)律這兩方面來驗(yàn)證本文根據(jù)實(shí)際工程建立的有限元模型的有效性。
以平臺寬度a為變量,按照表1 進(jìn)行參數(shù)取值,其他因素為定值,分別為φ=40°,β=30°,來模擬分析不同平臺寬度a對有限成坡土體滑裂面傾角的影響規(guī)律。 選取模型中Y=0 位置處的剖面來觀察滑裂面,不同有限填土上部寬度的最大剪切應(yīng)變速率云圖如圖5所示。
圖5 不同平臺寬度下有限成坡土體剪切應(yīng)變速率云圖
由圖5 可知,不同平臺寬度a的滑裂面與水平面的夾角隨著平臺寬度增大而逐漸減小,逐漸接近水平面,當(dāng)平臺寬度達(dá)到一定值時(shí),滑裂面沿斜向上方向發(fā)展。 滑裂面的剖面基本上呈直線形式,在坡體內(nèi)沿斜向下或斜向上的方向發(fā)展,逐漸接近坡面并實(shí)現(xiàn)貫通,該結(jié)果符合無黏性土情況下的平面滑裂面假定。利用該模型模擬得到的破裂角變化規(guī)律以及滑裂面的變化情況與文獻(xiàn)[6]結(jié)論一致,該模型的有效性得到驗(yàn)證。
根據(jù)式(1)[6],可以確定該模型下沿墻后土體深度h0分布的被動土壓力解析解。 平臺寬度a=9 m 時(shí)有限成坡土體被動土壓力解析解與數(shù)值解的變化規(guī)律對比如圖6 所示。
圖6 平臺寬度a=9 m 時(shí)有限成坡土體被動土壓力解析解與數(shù)值解對比
由圖6 可知,墻背上被動土壓力數(shù)值解和解析解均隨h0增大而增大,與實(shí)際土壓力分布規(guī)律一致。 在同一土體深度處,被動土壓力數(shù)值解和解析解基本一致,并且變化規(guī)律也基本一致,充分驗(yàn)證了該有限元模型的有效性。
4.1.1 平臺寬度
φ=40°,β=30°時(shí),不同平臺寬度a下有限成坡土體被動土壓力隨墻后土體深度h0的變化規(guī)律如圖7所示。
圖7 不同平臺寬度下有限成坡土體被動土壓力變化規(guī)律
由圖7 可知,不同平臺寬度a下有限成坡土體被動土壓力變化趨勢大致相同,均隨h0增大而增大,且為非線性變化。h0<2 m 時(shí),不同平臺寬度a下有限成坡土體被動土壓力大小基本一致;h0=2 ~10 m 時(shí),隨著平臺寬度a增大,有限成坡土體被動土壓力逐漸增大,但增長幅度越來越??;h0>10 m 后,隨著h0增大,不同平臺寬度a下有限成坡土體被動土壓力呈現(xiàn)陡升現(xiàn)象,在擋土墻最底部達(dá)到最大值。
4.1.2 有限土體內(nèi)摩擦角
a=3 m,β=30°時(shí),不同內(nèi)摩擦角φ下有限成坡土體被動土壓力隨墻后土體深度h0的變化規(guī)律如圖8所示。
圖8 不同內(nèi)摩擦角下有限成坡土體被動土壓力變化規(guī)律
由圖8 可知,不同內(nèi)摩擦角φ下有限成坡土體被動土壓力變化趨勢大致相同,均隨h0增大而增大,且為非線性變化;φ<40°時(shí),有限成坡土體被動土壓力呈“S”形增長,隨著h0增大,有限成坡土體被動土壓力增長速度先快后慢;φ=40°時(shí),h0>4 m 后有限成坡土體被動土壓力基本呈線性增長;φ>40°時(shí),有限成坡土體被動土壓力呈折線形式增長,增長速度在h0=4~10 m 時(shí)隨深度增大基本不變,在h0>10 m 后增長速度變大,出現(xiàn)陡升現(xiàn)象。 同一土體深度下,隨著內(nèi)摩擦角φ增大,有限成坡土體被動土壓力逐漸增大。
4.1.3 有限成坡土體邊坡坡角
a=3 m,φ=40°時(shí),不同邊坡坡角β下有限成坡土體被動土壓力隨墻后土體深度h0的變化規(guī)律如圖9 所示。
圖9 不同邊坡坡角下有限成坡土體被動土壓力變化規(guī)律
由圖9 可知,不同邊坡坡角β下有限成坡土體被動土壓力變化趨勢大致相同,均隨h0增大而增大,且為非線性變化。β>30°,h0<4 m 時(shí),不同邊坡坡角下的有限成坡土體被動土壓力大小基本一致,h0>4 m 后,不同邊坡坡角下有限成坡土體被動土壓力基本呈線性增長;β=30°,h0<4 m 時(shí),不同邊坡坡角下有限成坡土體被動土壓力增長速度較慢,h0>4 m 后,隨著h0增大,有限成坡土體被動土壓力基本呈線性增長;β<30°時(shí),隨著h0增大,有限成坡土體被動土壓力呈“S”形增長,增長速度先快后慢。 在同一土體深度處,隨著邊坡坡角逐漸增大,有限成坡土體被動土壓力逐漸減小。
4.1.4 小結(jié)
綜合以上分析可得,針對平動條件下的擋土墻,在不同平臺寬度a、不同內(nèi)摩擦角φ、不同邊坡坡角β影響下,墻后有限成坡土體被動土壓力沿墻后土體深度h0的分布規(guī)律基本不變。 邊坡坡角β對有限成坡土體被動土壓力的影響較大,β較大時(shí),在平動條件下,擋土墻向著土體方向產(chǎn)生較小位移,墻后有限成坡土體可能已經(jīng)被破壞,坡體的穩(wěn)定性難以保證。 為了保證露天礦開挖安全,應(yīng)盡量減小邊坡坡角β[8]。 根據(jù)模擬分析結(jié)果,并綜合考慮經(jīng)濟(jì)效益情況下,建議將邊坡坡角β取為30°,該結(jié)果也得到了實(shí)際工程的驗(yàn)證,說明本文數(shù)值分析結(jié)果的可靠性。
該工程措施在實(shí)際工作過程中,主要承受來自墻后有限成坡土體的被動土壓力,為了保證墻后有限成坡土體的穩(wěn)定性,針對此情況并結(jié)合圖7,墻后平臺寬度a=9 m 左右即可。 該分析結(jié)果與布設(shè)地下連續(xù)墻后邊坡上部平臺的寬度基本一致,進(jìn)一步證明了本文數(shù)值分析結(jié)果的可靠性。
4.2.1 平臺寬度
φ=40°,β=30°時(shí),不同平臺寬度a下有限成坡土體主動土壓力隨墻后土體深度h0的變化規(guī)律如圖10 所示。
圖10 不同平臺寬度下有限成坡土體主動土壓力變化規(guī)律
由圖10 可知,h0<10 m 時(shí),主動土壓力隨h0增大而增大;h0>10 m 后,主動土壓力隨h0增大而減小,變化曲線出現(xiàn)凸起。h0<8 m 時(shí),不同平臺寬度a的有限成坡土體主動土壓力的大小基本上呈線性增長;h0>8 m時(shí)出現(xiàn)陡升現(xiàn)象,同一土體深度處有限成坡土體主動土壓力大小隨著平臺寬度a增大而增大。
4.2.2 有限土體內(nèi)摩擦角
a=3 m,β=30°時(shí),不同內(nèi)摩擦角φ下有限成坡土體中的主動土壓力隨墻后土體深度h0的變化規(guī)律如圖11 所示。
圖11 不同內(nèi)摩擦角下有限成坡土體主動土壓力變化規(guī)律
由圖11 可知,h0<10m時(shí),主動土壓力隨h0增大而增大;h0>10 m 后,主動土壓力隨h0增大而減小,但減小幅度不同,隨著內(nèi)摩擦角φ增大,主動土壓力減小幅度逐漸減小,變化曲線出現(xiàn)凸起。 擋土墻平動條件下,隨著內(nèi)摩擦角φ增大,h0<10 m 時(shí),同一土體深度處的有限成坡土體主動土壓力整體呈減小變化。
4.2.3 有限成坡土體邊坡坡角
a=3 m,φ=40°時(shí),不同邊坡坡角β下有限成坡土體中的主動土壓力隨墻后土體深度h0的變化規(guī)律如圖12 所示。
圖12 不同邊坡坡角下有限成坡土體主動土壓力變化規(guī)律
由圖12 可知,h0<8 m 時(shí),不同邊坡坡角β下有限成坡土體主動土壓力均隨h0增大而增大,且大小基本一致;h0=8~10 m 時(shí),隨著邊坡坡角β增大,有限成坡土體主動土壓力逐漸減??;h0>10 m 后,隨著h0增大,有限成坡土體主動土壓力均減小,減小幅度隨著邊坡坡角β增大而逐漸增大。
4.2.4 小結(jié)
綜合以上分析可得,針對平動條件下的擋土墻,在不同平臺寬度a、不同內(nèi)摩擦角φ、不同邊坡坡角β影響下,墻后有限成坡土體主動土壓力沿墻后土體深度h0先增大再減小,與試驗(yàn)結(jié)果以及現(xiàn)場監(jiān)測結(jié)果基本一致[9-10]。 內(nèi)摩擦角φ對有限成坡土體主動土壓力分布值的影響較大。 在平動條件下,擋土墻向遠(yuǎn)離土體的方向產(chǎn)生位移時(shí),φ較小時(shí),擋土墻上主動土壓力較大,土體可能會發(fā)生破壞,邊坡穩(wěn)定性受到影響;φ較大時(shí),擋土墻上主動土壓力較小,墻后土體能夠較好地保持穩(wěn)定狀態(tài)。 為了使擋土墻和露天礦邊坡更好地保持穩(wěn)定狀態(tài),可以采取有效措施(如注漿[11]或置換[12]等)增大墻后填土內(nèi)摩擦角,從而減小墻背上的主動土壓力,為礦山安全開采提供保障。
1) 有限成坡土體滑裂面傾角隨平臺寬度增大先逐漸減小到零后逐漸增大,滑裂面由斜向下發(fā)展變?yōu)樾毕蛏习l(fā)展,平臺寬度a=9 m 條件下,利用該模型得到的墻背上被動土壓力數(shù)值解和已有推導(dǎo)公式得到的解析解基本一致,且變化規(guī)律基本相同,充分驗(yàn)證了本文數(shù)值模型的有效性。
2) 擋土墻在平動條件下,墻后有限成坡土體達(dá)到極限狀態(tài)時(shí),被動土壓力在不同平臺寬度a、不同內(nèi)摩擦角φ、不同邊坡坡角β影響下均隨h0增大而增大,呈非線性增加分布。 內(nèi)摩擦角φ和邊坡坡角β均會影響有限成坡土體被動土壓力的增長形式,平臺寬度a不會影響有限成坡土體被動土壓力的增長形式。 邊坡坡角β對墻后有限成坡土體被動土壓力的影響較大。 綜合考慮礦山實(shí)際情況,建議將邊坡坡角與平臺寬度分別定為30°和9 m。
3) 在不同平臺寬度a、不同內(nèi)摩擦角φ、不同邊坡坡角β影響下,墻后有限成坡土體主動土壓力沿h0的分布規(guī)律基本一致,在h0=10 m 處發(fā)生轉(zhuǎn)折,h0<10 m時(shí)擋土墻后主動土壓力均隨深度增大而增大,h0>10 m時(shí)均隨深度增大而減小。 有限填土上部寬度a和邊坡坡角β對有限成坡土體主動土壓力分布影響較小,內(nèi)摩擦角φ對有限成坡土體主動土壓力分布影響較大,為了減小擋土墻主動土壓力,可以采取有效措施(如注漿或置換等)增大墻后填土內(nèi)摩擦角,以保證擋土墻和邊坡的穩(wěn)定性,為礦山安全開采提供保障。