開挖巖體邊坡力學(xué)性狀及塑性區(qū)演變特征研究

陳志峰

巫山至巫溪(巫溪段)二級公路位于重慶市巫溪縣南部地區(qū)。起點位于龍溪金家溝，里程K86+000，終點在花栗路口，里程K109+875.993，全長23.875km。路線路段主要跨越大泉山山脈，地形較復(fù)雜，沿線多高陡邊坡，受巖性影響，有很多高邊坡在開挖過程中出現(xiàn)失穩(wěn)現(xiàn)象。研究路段位于柚子樹境內(nèi)，起止里程樁號K88+840～K88+950，全長110m。該段內(nèi)邊坡為路塹邊坡，切坡最高為32m，最低為21m。場地內(nèi)出露地層主要為第四系全新統(tǒng)殘坡積碎石土（Q4el+dl），厚度0.4~1.5m；以及三疊系中統(tǒng)巴東組（T2b3）泥灰?guī)r，厚度大于30m。

1理論模型

1.1 基本假設(shè)

(1) 處于塑性區(qū)(極限平衡區(qū))內(nèi)任何節(jié)點的的正應(yīng)力與剪應(yīng)力滿足莫爾-庫侖強(qiáng)度準(zhǔn)則：

其中，為斜面上的抗剪強(qiáng)度；，分別為巖體的粘聚力和內(nèi)摩擦角；為斜面上的正應(yīng)力，其值；，為方向余弦。

(2) 分析過程中不考慮體積力的影響。

(3) 巖體是彈塑性材料，且各向同性。

1.2 基本方程

在處于塑性區(qū)(極限平衡狀態(tài))的巖體，應(yīng)力滿足平衡微分方程：

1.3 本構(gòu)模型

本研究采用ansys有限元法求解邊坡穩(wěn)定問題時，采用了理想彈塑性模型，其本構(gòu)模型采Drucker-Prager 準(zhǔn)則：

式中：,分別表示應(yīng)力張量的第一不變量和應(yīng)力偏張量的第二不變量。、為與巖土材料內(nèi)摩擦角和粘聚力有關(guān)的常數(shù)，，。屈服面在平面上為不等角度的六邊形的外接圓。

2高邊坡開挖應(yīng)力場

為了解邊坡開挖的力學(xué)效應(yīng),對K88+680斷面邊坡進(jìn)行了二維彈塑性有限元計算分析。計算模型邊界為：底部為173m,高為112m,分別為開挖寬度和深度的4倍多，基本可以消除邊界對應(yīng)力的影響。底部取x、y方向

位移約束，側(cè)面施加x方向位移約束。模型網(wǎng)格的稀密對二維彈塑性有限元計算有著一定的影響，為了提高計算精度，所以在開挖區(qū)域及周邊敏感部位對網(wǎng)格采取加密措施。模型是由5544個節(jié)點組成的5688個單元。用“殺死”單元法開挖掉3710個單元。

2.1 計算參數(shù)的選取

本次模擬考慮到了地層巖性的差異，計算區(qū)域中所涉及的巖體主要有坡積土(Q4dl+el) 和強(qiáng)風(fēng)化泥灰?guī)r(T2b)，通過巖體物理力學(xué)試驗和工程地質(zhì)類比，最后確定了各巖體的計算參數(shù)(表1)

摘要：邊坡開挖過程的力學(xué)性狀變化是一個復(fù)雜的過程，同時其塑性區(qū)演變趨勢也是一個復(fù)雜的過程。文章通過假設(shè)邊坡完全處于理想彈塑性狀態(tài)，并以Drucker-Prager 準(zhǔn)則為本構(gòu)模型，運(yùn)用ANSYS有限元軟件對巫山至巫溪（巫溪段）公路K88段高邊坡開挖過程應(yīng)力調(diào)整過程及塑性區(qū)變化過程進(jìn)行了模擬，得出了開挖過程中應(yīng)力最大處為坡腳。邊坡塑性區(qū)是一個動態(tài)調(diào)整過程，最終位于強(qiáng)風(fēng)化泥灰?guī)r與弱風(fēng)化泥灰?guī)r交界處。

關(guān)鍵詞：開挖巖體邊坡；力學(xué)性狀；塑性區(qū)演變特征

巫山至巫溪(巫溪段)二級公路位于重慶市巫溪縣南部地區(qū)。起點位于龍溪金家溝，里程K86+000，終點在花栗路口，里程K109+875.993，全長23.875km。路線路段主要跨越大泉山山脈，地形較復(fù)雜，沿線多高陡邊坡，受巖性影響，有很多高邊坡在開挖過程中出現(xiàn)失穩(wěn)現(xiàn)象。研究路段位于柚子樹境內(nèi)，起止里程樁號K88+840～K88+950，全長110m。該段內(nèi)邊坡為路塹邊坡，切坡最高為32m，最低為21m。場地內(nèi)出露地層主要為第四系全新統(tǒng)殘坡積碎石土（Q4el+dl），厚度0.4~1.5m；以及三疊系中統(tǒng)巴東組（T2b3）泥灰?guī)r，厚度大于30m。

1理論模型

1.1 基本假設(shè)

(1) 處于塑性區(qū)(極限平衡區(qū))內(nèi)任何節(jié)點的的正應(yīng)力與剪應(yīng)力滿足莫爾-庫侖強(qiáng)度準(zhǔn)則：

其中，為斜面上的抗剪強(qiáng)度；，分別為巖體的粘聚力和內(nèi)摩擦角；為斜面上的正應(yīng)力，其值；，為方向余弦。

(2) 分析過程中不考慮體積力的影響。

(3) 巖體是彈塑性材料，且各向同性。

1.2 基本方程

在處于塑性區(qū)(極限平衡狀態(tài))的巖體，應(yīng)力滿足平衡微分方程：

式中：,分別表示應(yīng)力張量的第一不變量和應(yīng)力偏張量的第二不變量。、為與巖土材料內(nèi)摩擦角和粘聚力有關(guān)的常數(shù)，，。屈服面在平面上為不等角度的六邊形的外接圓。

2高邊坡開挖應(yīng)力場

為了解邊坡開挖的力學(xué)效應(yīng),對K88+680斷面邊坡進(jìn)行了二維彈塑性有限元計算分析。計算模型邊界為：底部為173m,高為112m,分別為開挖寬度和深度的4倍多，基本可以消除邊界對應(yīng)力的影響。底部取x、y方向

位移約束，側(cè)面施加x方向位移約束。模型網(wǎng)格的稀密對二維彈塑性有限元計算有著一定的影響，為了提高計算精度，所以在開挖區(qū)域及周邊敏感部位對網(wǎng)格采取加密措施。模型是由5544個節(jié)點組成的5688個單元。用“殺死”單元法開挖掉3710個單元。

2.1 計算參數(shù)的選取

本次模擬考慮到了地層巖性的差異，計算區(qū)域中所涉及的巖體主要有坡積土(Q4dl+el) 和強(qiáng)風(fēng)化泥灰?guī)r(T2b)，通過巖體物理力學(xué)試驗和工程地質(zhì)類比，最后確定了各巖體的計算參數(shù)(表1)

巖體 (kg/m3) (MPa) (KPa)

坡積土 2100 150 0.38 25 21

強(qiáng)風(fēng)化泥灰?guī)r 2300 1800 0.32 60 23

弱風(fēng)化泥灰?guī)r 2400 2600 0.26 130 24

灰?guī)r(基巖) 2400 4500 0.18 400 35

2.2 初始應(yīng)力場

巖體的初始應(yīng)力場, 取正應(yīng)力以壓為正（在ansys中數(shù)值上表示為負(fù)）, 其大主應(yīng)力方向在近地表呈不規(guī)則的鋸齒型，深部接近水平，深部應(yīng)力值為2.6MPa左右(圖1)；小主應(yīng)力方向在近地表處基本上與坡面輪廓線平行, 深部接近水平，深部應(yīng)力值為0.57MPa左右(圖2)。

2.3 開挖應(yīng)力場

經(jīng)計算, 在開挖過程中, 初始應(yīng)力場不斷受到擾動與調(diào)整, 開挖區(qū)左右及下部的擾動范圍在1倍開口寬度以內(nèi), 開挖面附近大主應(yīng)力方向接近垂直開挖面方向, 坡面局部地區(qū)由壓應(yīng)力變?yōu)槔瓚?yīng)力,小主應(yīng)力方向接近于平行開挖面方向,路基近表面是處于受拉狀態(tài)。在左右坡腳處存在不同程度的應(yīng)力集中現(xiàn)象,右邊坡腳處最大壓應(yīng)力值為0.17MPa，左邊坡腳處最大壓應(yīng)力0.98MPa。

2.4 受拉區(qū)域

邊坡切削完成后，由于卸荷回彈，整個路基近表面及左切坡的第一臺階中部出現(xiàn)了拉應(yīng)力區(qū)，并受地質(zhì)構(gòu)造等因素的影響而呈現(xiàn)出不同的分布形式。在路基上出現(xiàn)呈矩形狀拉應(yīng)力區(qū)，其大主應(yīng)力值為30～75kPa，左切坡上有一個呈等邊三角形狀的拉應(yīng)力區(qū)，其大主應(yīng)力值為5～27kPa。

3塑性區(qū)演變特征

根據(jù)計算區(qū)內(nèi)的地層結(jié)構(gòu)，在建模的時候分成四層不同巖性的巖體，該邊坡分四個階段開挖來分析其塑性區(qū)的演變情況。隨開挖深度的增加, 塑性區(qū)范圍不斷增大。當(dāng)開采深度達(dá)到一定深度, 邊坡的穩(wěn)定性就會受到很大威脅。

3.1 第一次開挖

第一次開挖主要是把近地表的坡積土挖除，兼挖強(qiáng)風(fēng)化泥灰?guī)r層上部，因為原地表斜坡比較陡峭，且開挖坡比高大(1∶0.5~1∶0.3)。此部分開挖后，觀察有限元計算結(jié)果，可以看出在坡腳除出現(xiàn)了應(yīng)力集中，并有小范圍的塑性變形，其等效塑性應(yīng)變區(qū)的值為：0.395×10-5～0.673×10-4 。

3.2 第二次開挖

第二次開挖是將強(qiáng)風(fēng)化泥灰?guī)r層切掉，坡腳已達(dá)到強(qiáng)風(fēng)化泥灰?guī)r與弱風(fēng)化泥灰?guī)r交界處，坡腳雖然有應(yīng)力集中現(xiàn)象，但是塑性區(qū)是出現(xiàn)在層間，與第一次開挖的塑性區(qū)是相連接的。塑性區(qū)外的彈性區(qū)應(yīng)力有增加的趨勢，這是因為此處巖體發(fā)生塑性變形，將應(yīng)力釋放轉(zhuǎn)移到彈性區(qū)巖體內(nèi)。塑性應(yīng)變出現(xiàn)在坡面臨空面，其等效塑性區(qū)內(nèi)的值為：0.163×10-4～0.277×10-3 。

3.3 第三次開挖

第三次開挖是沿第二次開挖的基礎(chǔ)上往下開挖5m左右，坡角處出現(xiàn)應(yīng)力集中現(xiàn)場，并出現(xiàn)小范圍的塑性區(qū)，在強(qiáng)風(fēng)化泥灰?guī)r與弱風(fēng)化泥灰?guī)r的交界處的塑性區(qū)范圍進(jìn)一步擴(kuò)大，并出現(xiàn)滑移變形。層間塑性區(qū)的等效塑性應(yīng)變?yōu)椋?.124×10-4～0.210×10-3；開挖坡腳塑性區(qū)的等效塑性應(yīng)變?yōu)椋?.124×10-4～0.111×10-3。

3.4 第四次開挖

此次開挖是將邊坡切削到路基設(shè)計標(biāo)高，整個路基是處在弱風(fēng)化泥灰?guī)r層中。強(qiáng)風(fēng)化泥灰?guī)r和弱風(fēng)化泥灰?guī)r交界處出現(xiàn)大面積的塑性區(qū)，坡腳排水溝處出現(xiàn)應(yīng)力集中，有塑性區(qū)分布，并在路肩上了出現(xiàn)塑性變形。第三階段的開挖坡腳塑性區(qū)消失，這是因為隨著開挖的深入，此處的應(yīng)力集中消失。層間塑性區(qū)的等效塑性應(yīng)變?yōu)?.208×10-4～0.187×10-3 ，其值較第三階段小，這是由于開挖卸荷后應(yīng)力調(diào)整，使得部分變形反彈；坡腳塑性區(qū)的等效塑性應(yīng)變?yōu)?.208×10-4～0.229×10-3 。

4結(jié)論

（1）隨著自上而下開挖推進(jìn), 應(yīng)力不斷調(diào)整，位移、塑性區(qū)范圍也不斷增大。開挖結(jié)束后, 左邊坡腳處最大壓應(yīng)力0.98MPa 。雖然坡體總體上處于穩(wěn)定狀態(tài), 但通過對應(yīng)力、位移及塑性區(qū)計算結(jié)果分析知道,坡面出現(xiàn)了拉應(yīng)力區(qū),有局部破壞的危險，應(yīng)力集中區(qū)坡腳處。

（2）隨著開挖的進(jìn)行，在強(qiáng)風(fēng)化泥灰?guī)r與弱風(fēng)化泥灰?guī)r交界處塑性區(qū)范圍逐漸擴(kuò)大，從云圖中可以看出有向下滑移的趨勢；坡腳處塑性區(qū)隨著開挖的深入而位置也發(fā)生移動。