張典典 雷 浩 吳月勇
(1.河海大學(xué) 港口航道與近海工程學(xué)院,江蘇 南京 210098;2.河海大學(xué) 土木與交通學(xué)院,江蘇 南京 210098)
土力學(xué)中的瑞典條分是邊坡穩(wěn)定分析的一種經(jīng)典算法,目前仍然是工程上經(jīng)常應(yīng)用的一種方法。但是其計算過程較為復(fù)雜,尤其是在進(jìn)行最危險滑弧的搜索時,每給定一個新的滑弧,都要重新進(jìn)行條分和確定參數(shù)。在傳統(tǒng)的土力學(xué)教學(xué)方式下,學(xué)生主要把時間都花費在了繁瑣的數(shù)據(jù)計算上,讓很多學(xué)生失去了學(xué)習(xí)土力學(xué)的興趣。而MATLAB 軟件具有高效的數(shù)值計算能力和完備的圖形功能。將MATLAB 運用到土力學(xué)的學(xué)習(xí)和研究中,不僅可以減輕學(xué)生處理數(shù)據(jù)的壓力,還可以提高學(xué)生應(yīng)用計算機處理數(shù)據(jù)的能力。
條分法就是是先假定若干可能的剪切面、滑裂面。然后將滑裂面以上土體分成若干垂直土條,對作用于各土條上的力進(jìn)行力與力矩的平衡分析,求出在極限平衡狀態(tài)下土體穩(wěn)定的安全系數(shù),并通過一定數(shù)量的試算,找出最危險滑裂面位置及相應(yīng)的(最低的)安全系數(shù)。
而瑞典條分法是最古老的方法之一,除假定滑動面在圓柱面及滑動土體為不變的剛體外,還假定不考慮土條側(cè)面上的作用力。計算公式為:
式中:Fs為邊坡穩(wěn)定系數(shù);c 為邊坡土體內(nèi)聚力;φ 為邊坡土體內(nèi)摩擦角;wi為第i 條土條重;αi為第i 條土條底滑面的傾角;li為第i 條土條底滑面的長度。
具體計算土坡的穩(wěn)定安全系數(shù)時,要先按比例繪出土坡,在可能滑動范圍內(nèi)選取圓心O1,作出相應(yīng)的滑動圓弧。然后將土體分成n 個土條,量出各土條中心高度hi,寬度bi并列表計算sinαi,cosαi以及wi等值,計算與圓心、半徑對應(yīng)的安全系數(shù)。對圓心O1選不同半徑,得到O1對應(yīng)的最小安全系數(shù)。在可能滑動范圍內(nèi),選擇其他圓心,重復(fù)上列計算,從而求出最小的安全系數(shù),即為該土坡的穩(wěn)定安全系數(shù)。
用手算完成上述過程,計算量是相當(dāng)大的。
若已知一粘性土土坡,坡高為h,坡比為b,填土粘聚力為c,內(nèi)摩擦角為φ,重度為γ,試計算土坡的穩(wěn)定安全系數(shù)。
用MATLAB 求解程序如下:
function [Fsmin,xb,yb,Rb]=ruidian(b,h,gama,phi0,c)
%輸入 b:坡度h:坡高gama:土的重度phi0:土的內(nèi)摩擦角c:土的粘聚力
%輸出 Fsmin:最小安全系數(shù)xb,yb:圓心坐標(biāo)Rb:滑弧半徑
我們只需將坡度,坡高,土的重度,內(nèi)摩擦角,土的粘聚力的數(shù)值依次輸入,然后運行程序,將直接得到安全系數(shù)。
已知h 為15m,坡比b 為1:2,填土粘聚力c 為40kPa,內(nèi)摩擦角φ為8 度,重度γ 為19.5kN/m3。將以上數(shù)據(jù)帶入程序,得到結(jié)果:
Fsmin=1.3573;xb=13.0884;yb=23.9249;Rb=27.2710。
與大量手算后的結(jié)果相同,說明該程序是正確的。
傳統(tǒng)手算時我們一般是列表計算,但是當(dāng)n 非常大時,數(shù)據(jù)將非常多,手算很容易出錯,且耗時長。如果運用了以上MATLAB 程序,我們只需輸入已知值,將直接得到結(jié)果,不僅快捷方便而且不會出錯。
本文成功的將MATLAB 應(yīng)用在瑞典條分法中,操作方便,且不易出錯。除了本文中的瑞典條分法,土力學(xué)中的地基豎向附加應(yīng)力計算、地基沉降計算等同樣可以運用MATLAB 解決。
[1]錢家歡,殷宗澤.土工原理與計算[M].北京:中國水利水電出版社,1996.
[2]錢家歡.土力學(xué)[M].南京:河海大學(xué)出版社,1995.
[3]清源工作室.MATLAB6.0 基礎(chǔ)及應(yīng)用[M].北京:機械工業(yè)出版社,2001.