第23卷第3期巖土工程學(xué)報(bào)Vol.23 No.32001年.5月Chinese Joumal of Geotechnical EngineeringMay，2001粘土的流變特性分析Research on rheology character of clay陳鐵林陳生水周成沈珠江(南京水利科學(xué)研究院江蘇南京210024)摘要將聚合物網(wǎng)絡(luò)理論和速率過(guò)程理論應(yīng)用到粘土的流變特性分析中得出了--些有益的結(jié)論將雙電層理論應(yīng)用到實(shí)驗(yàn)的設(shè)計(jì)中得出了和理論分析一致的結(jié)果并從理論上研討了砌塊模型的正確性。關(guān)鍵詞聚合物網(wǎng)絡(luò)理論速率過(guò)程理論流變觸變中圖分類號(hào):TU 44文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào):1000- 4548( 2001 )3 - 0279-05作者簡(jiǎn)介陳鐵林男,1970 年生在讀博土。CHEN Tie-lin CHEN Sheng shui ZHOU Cheng SHEN Zhu- jiang( Nanjing Hydraulic Research Institute ,Nanjing 210024 ,China )Abstract :A study is made on rheology character of clay which is assumed to obey the worknet theory of polymer and rate process theory. Thenexperiments designed by using the theory of double layer have been performed. The theoretical predictions is found to agree with experimental re-sults . The study theretically support the masonry model for structural clays.Key words :worknet theory of polymer irate process theory ;theology athixotrpy土體中結(jié)構(gòu)的變化,從而產(chǎn)生了觸變現(xiàn)象。觸變-詞1引言'由A.F. Petefi在1927年提出41。H. Freundlich用觸變粘土的流變性質(zhì)在很大程度上取決于土體內(nèi)部缺來(lái)描述膠體懸液中凝膠-溶膠相互轉(zhuǎn)變的可逆過(guò)陷(孔隙)的存在以及土顆粒邊界性質(zhì)(形狀和表面特程5。J.M. Burgers and G. W. Sctt-Blair 把觸變"定義性)》可以認(rèn)為粘性變形是由顆粒邊界滑移造成的。為重塑造成軟化之后,和時(shí)間相關(guān)的硬化過(guò)程"[6]。結(jié)構(gòu)性粘土流變的特殊性在于伴隨著顆粒的破損接0. Morettd7] A. W . Skemptom和R. D. NortheF 8]以及H.觸點(diǎn)的增加(網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)的增加和與之相伴的連通性B. Seed和C. K. Charf9]對(duì)壓縮粘土的觸變性進(jìn)行了研的變化)。為了把位錯(cuò)運(yùn)動(dòng)(滑移)和缺陷擴(kuò)散(結(jié)構(gòu)破究。特別是以下兩個(gè)結(jié)論為本文的研究提供了依據(jù)。損)以及損傷的概念應(yīng)用到對(duì)結(jié)構(gòu)性粘土的研究中需Kruyt通過(guò)實(shí)驗(yàn)指出觸變材料中顆粒之間相互接觸形要從理論和實(shí)驗(yàn).上同時(shí)進(jìn)行研究。本文進(jìn)行了初步的成網(wǎng)狀結(jié)構(gòu)10。H. Van Olpherf11] ,R. Fahn ,A Weiss和嘗試將聚合物網(wǎng)絡(luò)理論和速率過(guò)程理論應(yīng)用到對(duì)粘Hoffmart 12也都得出了類似的結(jié)論。James K. Mitchell土流變特性的研究中,得出了-些有益的結(jié)論并在實(shí)認(rèn)為在發(fā)生觸變的粘土中,當(dāng)粘土顆?？梢宰杂蛇x擇驗(yàn)中得到了驗(yàn)證。位置時(shí)，這些顆?？偸潜憩F(xiàn)為形成凝絮狀結(jié)構(gòu)的傾向3。Freundliclf 5 ]指出觸變行為無(wú)疑與顆粒之間作用2綜述力的平衡有關(guān);電解質(zhì)的濃度有顯著的影響;任何顆粒土的結(jié)構(gòu)是指土顆粒或集合體的大小、形狀、表面排列發(fā)生的變化都和雙電層有關(guān)。特征、排列形式以及它們之間的連接特征,它對(duì)土的力土是缺陷體材料在受到外力作用時(shí)總會(huì)有許多學(xué)特性有著重要的影響1。土的結(jié)構(gòu)大體上可分為三膠結(jié)的斷裂和顆粒的破碎必然伴隨著內(nèi)應(yīng)力的調(diào)整。類:單粒結(jié)構(gòu)、凝絮狀結(jié)構(gòu)和凝散狀結(jié)構(gòu)。粘土的工程在-定的條件下這些內(nèi)應(yīng)力的調(diào)整與時(shí)間相關(guān)就會(huì)特性由現(xiàn)時(shí)結(jié)構(gòu)決定,而現(xiàn)時(shí)結(jié)構(gòu)依賴于初始結(jié)構(gòu)和產(chǎn)生流變現(xiàn)象和某些流變現(xiàn)象的耦合如蠕變-應(yīng)力它的變化。換句話說(shuō)粘土特性的變化遵循以下順?biāo)沙谕瑫r(shí)存在于十體中。把膠體化學(xué)與速率過(guò)程理論序2]①壓實(shí)或沉積條件決定初始結(jié)構(gòu)②初始結(jié)構(gòu)發(fā)引入對(duì)中國(guó)煤化工勺研究已經(jīng)取得了成生變化形成現(xiàn)時(shí)結(jié)構(gòu)如圖1。功13-16MHCNMHG文把聚合物網(wǎng)絡(luò)理.在本文實(shí)驗(yàn)的設(shè)計(jì)中考慮了觸變因素,因此先對(duì)論引入對(duì)粘土流變的研究17-21]。觸變的研究做一些回顧。 觸變現(xiàn)象-般發(fā) 生在具有凝絮狀結(jié)構(gòu)的粘土中3]。Jamnes K. Mitchell認(rèn)為擾動(dòng)造成基金項(xiàng)目國(guó)家自然科學(xué)基金資助項(xiàng)目( 19772019 )顆粒之間作用為的不平衡擾動(dòng)停止后不平衡力導(dǎo)致收稿日期2000-09- 28280巖土工程學(xué)報(bào)2001年形成因素變化因素(壓 力同時(shí)Po=勢(shì)(5)時(shí)間由于P= Po+△P(6)(壓實(shí)粘土環(huán)境變化壓縮方式和大小壓力孔隙水的特性含水量可得「” (y)xiy= 5+ 5[' rn(γ)xiy(7)2。t溫度電解性濃度其它初始結(jié)構(gòu)+| 陽(yáng)離子價(jià)=現(xiàn)時(shí)結(jié)構(gòu)求導(dǎo)可得d_三__ Krn(r)(8)沉積粘土陰離子大小運(yùn)積因素雙電層特性沉積過(guò)程陰離子通解為(r)= γ[ C- An( ry)](9).沉積環(huán)境pH值則r以y)== C- An(γ)( 10)\擾動(dòng)γ圖1土結(jié)構(gòu)的變化邊界條件為Fig. 1 Development of soil structure(1)y=0時(shí),(0)= ηo=C ;3對(duì) 觸變現(xiàn)象的研究(2)y> Ym時(shí),η=7m= ηo- :Knmax物體的完善流變曲線如圖2所示。曲線rτ)與其中no為結(jié)構(gòu)未 破壞時(shí)的粘度; ηm為結(jié)構(gòu)完全破γ軸之間包括的面積等于產(chǎn)生定常流動(dòng)需要的功率,壞時(shí)的粘度;nmaxk 為η = ηm時(shí)破壞膠結(jié)點(diǎn)總數(shù)。式( 10 )又可記為包括兩部分:-部分是保持具有Newton粘度流動(dòng)所需功率Po ;另一部分是在流動(dòng)過(guò)程中單位時(shí)間內(nèi)土體結(jié)爪r)= ηo-( ηo- 7)4Qr)( 11)構(gòu)破損所需的功率△P:其中ψ y)=KY當(dāng)ψ y)=1時(shí)結(jié)構(gòu)完全破壞。nmaxP=I°t(rXiγ(1)令4( γ)=1-ψ i)代入式11)得i↑以r)= 7m+( ηo- 7m)4( Y )( 12)或ηr)= 7m+Oη( 13)R%式中4( r稱為觸變恢復(fù)函數(shù);Aη是結(jié)構(gòu)粘度增量。由以上分析還可以得出Kn maxη0-7m=( 14)式( 14表明結(jié)構(gòu)強(qiáng)度越高( K值越大)達(dá)到結(jié)構(gòu)完全--7破損時(shí)破壞膠結(jié)點(diǎn)總數(shù)越多所需時(shí)間越短(剪切應(yīng)變速率越大)則結(jié)構(gòu)性表現(xiàn)得越明顯。圖2物體的流變曲線根據(jù)速率過(guò)程理論[ 16] ,如圖3。Fig.2 Rheological curve能以n表示單位土體中膠結(jié)點(diǎn)的數(shù)量則dn( γ )表示應(yīng)變率增量dγ引起膠結(jié)點(diǎn)破壞的數(shù)量,由因次分析可知一個(gè)膠結(jié)點(diǎn)破壞所需的平均功率為(AP)=張(2)21K與γ無(wú)關(guān)與土的結(jié)構(gòu)強(qiáng)度、干密度、顆粒大小等因素有關(guān)結(jié)構(gòu)強(qiáng)度越大,K值越大其因次為質(zhì)量長(zhǎng)度。中國(guó)煤化工產(chǎn)生增量dy所需平均功率為fYH.CNMHG還十過(guò)任理論dOP)= LK. dndy= r"Kn(7)dy (3)Fig.3 Rate process theoryly剪切速率V所以’=氣[ rn(y)ir(4)γ=2 hex(-的)inl(τ( 15)λh202h0第3期陳鐵林等.粘土的流 變特性分析281式中k為Bolzmann常數(shù);0為絕對(duì)溫度;h為Planck兩個(gè)不同的范疇分為兩類結(jié)構(gòu)重整①土顆?；圃斐?shù);R為氣體常數(shù);iF為克服能柵耗費(fèi)(包括結(jié)構(gòu)成接觸點(diǎn)的變化屬拓?fù)鋭?dòng)力學(xué)②伴有接觸點(diǎn)數(shù)變化破損)的功;為兩個(gè)平衡點(diǎn)之間的距離;入為流動(dòng)點(diǎn)的顆粒破損屬擴(kuò)散動(dòng)力學(xué)。之間,垂直于剪切面方向的距離,λ/λ與顆粒形狀有定義二次矩張量M為關(guān);V;為流動(dòng)單元的體積。M: = |m(r t )xndir( 20)V{hT^2k60式中m( r ,t 為材料接觸點(diǎn)分布函數(shù)ix; 表示r的一Aη==v:-ex(-能)xinK( 20 )AFτsVry般分量。定義參數(shù)Mo為τs/BMo= mf(rX 2x)XPr(21 )( 16)sir( τs/B)2h0式中m( r表示材料的靜態(tài)各向同性。式中τs為作用于結(jié)構(gòu)粘度增量的應(yīng)力;B=T% i再定義J=M表示歸一化矩張量。應(yīng)力則可表Mo=y.Hex(的)示為代入式13 )得T= 2bk0M,J + pI .(22)式中b為材料常數(shù)ip為球應(yīng)力;I為單位矩陣。”= 7m+Asinτs/B)將二次矩張量J直接寫成運(yùn)動(dòng)方程τs/B _τ/B _因?yàn)閟inh(τs/B)≈sinh( τ/B)同時(shí)下列邊界條件得d.J= H(J)+ LJ + JLT(23 )到滿足①當(dāng)τ=τs= r=0時(shí),η= η= 7m+A②式中L為應(yīng)變速率;H(J)為重整化張量表示結(jié)構(gòu)當(dāng)τ→∞r(nóng)τs→0,Y→∞Aη→0時(shí),η= 7mo調(diào)整后兩項(xiàng)表示土體的仿射變形滑移2所以由式12)可得二次矩張量M的運(yùn)動(dòng)方程的一般形式是η= 7m+(ηo- 7m)sir( τ/B )IB( 17)dM dMJM(24)lt F dt|,+ Jt對(duì)比式17)式(12和式15可得JM |4》)=tV2X λ/λ )( 18))t| ,是由結(jié)構(gòu)重整引起的,即t I,= H( m)iat式18表明對(duì)于不同的土達(dá)到同樣剪切應(yīng)變速率時(shí)是由土體仿射變形引起的即M = LM+ ML'。t|需要的剪切力越大觸變恢復(fù)函數(shù)越大,在同樣的剪切以上分析與沈珠江提出的砌塊模型22]是-致的。力作用下獲得的剪切應(yīng)變速率越大,觸變恢復(fù)函數(shù)越原狀粘土中第二類結(jié)構(gòu)重整，必然會(huì)增加理論分析的小;λ/λ表明觸變恢復(fù)函數(shù)還與顆粒形狀、排列有關(guān)。復(fù)雜性。以下對(duì)重塑土在不排水條件下流變特性進(jìn)行如對(duì)于橢圓形顆粒當(dāng)長(zhǎng)軸垂直于大主應(yīng)力時(shí)入/λ1達(dá)分析。到最大值ry( γ )也達(dá)到最大值結(jié)構(gòu)最穩(wěn)定;肖長(zhǎng)軸平首先認(rèn)為重塑土是均勻各向同性的,再假定重塑行于大主應(yīng)力時(shí)入/入1達(dá)到最小值結(jié)構(gòu)最不穩(wěn)定, V;土在變形過(guò)程中不存在顆粒的破碎,不排水條件下顆表示觸變恢復(fù)函數(shù)與材料的缺陷(或強(qiáng)度分布不均勻)粒體積密度不變。則以下條件需得到滿足:①等容變有關(guān)。缺陷造成應(yīng)力集中使V;較小則4('γ)也較小。形條件t( γ )=0 ;②統(tǒng)計(jì)不變性條件de( J)=1，這個(gè).對(duì)比式17 )和式16 )還可以得到不變性不涉及顆粒排列的變化,而僅涉及結(jié)構(gòu)參數(shù)的平均值。1入1er AF70-7m=v.→ex(RE).(19)下面對(duì)兩種情況進(jìn)行分析。較大的ηo- nm 表示材料結(jié)構(gòu)性表現(xiàn)明顯。式( 19表取(25)明材料結(jié)構(gòu)性隨著顆粒長(zhǎng)短軸比的增大,以及長(zhǎng)軸與中國(guó)煤化工g^大主應(yīng)力方向的接近而變得顯著。另外材料的缺陷和其中MCNMHG數(shù)[17]。.式25 )兩定山、④網(wǎng)個(gè)余什,歸一化矩張量J的能柵的提高都會(huì)使結(jié)構(gòu)性表現(xiàn)得更加明顯。非零張量分別表示為J小J22、J3和J12 ,以y表示剪切4對(duì)定常剪切流變強(qiáng)度的分析應(yīng)變速率。再引入歸一化剪切速率a= γto和歸-化在原狀和整押發(fā)生的重整效應(yīng)(結(jié)構(gòu)調(diào)整)應(yīng)屬于時(shí)間tn= l/t0o282巖土工程學(xué)報(bào)2001年由式23和( 25 )以及條件②得到根據(jù)式23 )定常剪切方程可寫為dJ1=-(Jn- 3)+ 2aJ12y J + JyT+ H(J)=0(31)dtn8此式的數(shù)值解如圖5。由圖中可以看出剪切強(qiáng)度不再dJ22=_ 2+3g隨剪切速率提高而單調(diào)增加a≈9.6時(shí)強(qiáng)度達(dá)到最大值在此以后強(qiáng)度隨剪切速率的增大而減小。這是因dJ33 .=-J33 +3(26)為觸變效應(yīng)減少造成的。以上分析了不存在顆粒破損的情況。對(duì)于第二類結(jié)構(gòu)調(diào)整的研究必須引入和砌dJ12dt- J12塊模型類似的損傷概念。1.0葉g = J1J2-后+ J:6( J1+ J2)取初始條件JI<0)=0 ,J1(0)= J260)= J36(0)=1。計(jì)算結(jié)果如圖4所示。0.5-1. 5廠1.0|2406-a=2.5a=1.3圖5 Jn-a 關(guān)系曲線Fig.5 The curveof J2-a05實(shí)驗(yàn)結(jié)果圖4不同速率下Jn2-tn 曲線觸變現(xiàn)象-般發(fā)生在凝絮狀結(jié)構(gòu)的土中，這種結(jié)Fig.4 The curve of Jr2- In at diferent strain rate構(gòu)形成時(shí)顆粒之間以吸力為主顆粒呈隨機(jī)排列。根從圖中可以看出隨著歸-化剪切應(yīng)變速率的增據(jù)雙電層理論雙電層越厚越不易形成凝絮狀結(jié)構(gòu)加強(qiáng)度提高而軟化現(xiàn)象明顯。也就是說(shuō)在高剪切速而離子濃度的提高,可以減少雙電層的厚度因此為形率下重塑土也可能會(huì)出現(xiàn)軟化現(xiàn)象。這一趨勢(shì)在下面成凝絮狀結(jié)構(gòu)必須保證-定的離子濃度和足夠多的的實(shí)驗(yàn)中得到了證實(shí)。片狀顆粒。另外如果土的滲透系數(shù)太小在較高的剪從式26 )還可得出定常剪切方程切速率下會(huì)產(chǎn)生不均勻孔壓而限制了剪切速率的選J1 -二= 2aJ12J22--=0]擇范圍只能在較低的速率下實(shí)驗(yàn)使得- -些性狀得不(27)到表現(xiàn)。以上幾點(diǎn)是實(shí)驗(yàn)過(guò)程中重點(diǎn)考慮的問題。J33-2=(J12=aJm.試樣的制備采用人工制備結(jié)構(gòu)性粘土的方法231。式27的解為試樣中的水泥成份,保證了-定的離子濃度大孔隙和高嶺土的存在提高了土樣的滲透系數(shù);高嶺土的存在J1=(a2+ 1)3(2a2+ 1)還提高了片狀顆粒的比例。電鏡掃描照片也證明了凝J22 = J33 =(a2+ 1)言(28)絮狀結(jié)構(gòu)的形成[23。土樣的滲透系數(shù)為8.112x 10-5J12 = a(a2+ 1)3en/s。根據(jù)規(guī)范24],本次實(shí)驗(yàn)最大應(yīng)變率為0.5%/mino分別對(duì)結(jié)構(gòu)性粘土和重塑粘土進(jìn)行了3種圍壓將式28 )代入式( 22 )得到τ12 = 2bk0M{[d(a2+ 1)3](29)下4種剪切應(yīng)變速率共24組CU實(shí)驗(yàn)。以下圖中圖例順序依次為土樣類型、前期固結(jié)壓力、剪切應(yīng)變速率,考慮到r以a)=三,由式(29)得以a )xa-3 ,可見粘度其中J表示結(jié)構(gòu)土,C表示重塑土。圖6是應(yīng)變速率η只隨a的增大而降低,而沒有觸變項(xiàng)。為了分析具-剪切中國(guó)煤化工果(圖5)是一致的。，有觸變性的剪切流動(dòng)過(guò)程考慮到能量耗散必須是彈從圖7qCNMH G中應(yīng)變速率越大軟性能本身的嚴(yán)格增函數(shù)取新的模式17]化現(xiàn)象越明顯在本次實(shí)驗(yàn)選定的速率范圍內(nèi)重塑土雖然沒有出現(xiàn)軟化現(xiàn)象但這種趨勢(shì)也很明顯(如圖7HKJ)=-][6r+(8+2)J-引(30) .(b)屈服以后,速率越大,曲線越平緩)這和圖4的其中l(wèi)互方數(shù)據(jù)理論分析結(jié)果是一致的。第3期陳鐵林等.粘土的流變特性分析283[ 3 ] Mitchell James K. Fundamental aspects of thixotropy in soil[ J].5r士{ 25kPeJoumal of the Soil Mechanics and Foundations ,1960 ( SM3 ):19二{ 100kPe~ 52..... J 200kPa: 200kPa[4 ] Peterfi Arch F. Entwichlungsmech[ J]. D Organism ,1927 ,112 :689[ 5 ]Fruendlich H , Hermann et cie. Thixotropy[ M] Paris , 1935.[6] Burgers J M Scott Blair G W. Report on the principles of rtheo-0.1。2.30.4 0.5logical nomenclature [ A ]. Proceedings Internatl Rheologici/(103n)congress Joint Committee on Rheology of the Intematl Council of圖6強(qiáng)度-應(yīng)變速率關(guān)系曲線scientific Union[ C]. Amsterdam , 1948.Fig.' 6 Strength - strain rate curves[7] Moretto O. Effect of natural hardening on the unconfined com-120pJ 100kPa 0. 075%/minpression strength of remoulded clay{ A ]. Proceedings Second In-100tematl Conference on Soil Mechanics :Vol I[ C ].1948.J 100kPa 0. 25%/min炙80. J 100kPa 0. 5%/min[8] Skempton A W ,Northey R D. The sensitivity of clays[ J ]量60-/1Geotechnique ,1952 3 1).J 100kPe 0. 045%/min[9]Seed H B ,Chan C K. Thixotropic Characteristics of compactedclay{ J] Proceedings , ASCE ,1957 83( SM4 ).[ 10 ]Kruyt H R. Colloid Science ,I , Ireversible System([ M] New0 24681012141618202York Elsevier Pub Co , 1952.應(yīng)變/%[ 11 ]Van Olphen H. Forces between suspended bentonite particles&)結(jié)構(gòu)土(圍壓100kPa)90p[ A] Clays and Clay Minerals , Pocedings , Fourth Natl Con-80-C 25kPa 0. 075%/minference on Clays and Clay Mineral[ C ]. Publication 456 , Natl.Academy of Sciences , Natl. Research council ，1956.生6050-[ 12 ]Fahn R ,Weiss A ,Hofmann H. Thixotropy in clay[ J ] BerichteC 25kPa 0. 25%/min創(chuàng)40-der Deutsche Keramischen Gesellschaft ,1953 30 21.寒30/C 25kPa 0.5%/min[ 13 ]Feda J. Interpretation of creep of soils by rate process theory[ J] .20C 25kPa 0. 045%/minGeotechnique , 1989 (4):667 ~ 677.10[ 14 Mitchell1J K. Shearing resistance of soils as a rate proces[ J].246.810121416Jourmal of Soil Mechanics and Foundations ,1964 ( SM1 ):29 ~(b)重塑土(圈壓25kPa)61[ 15 ]Anclersland O B ,Douglas A G. Soil deformation rates and acti-圖7不同速率下應(yīng)力-應(yīng)變曲線Fig.7 Deviator stress - strain curves at diferent strain ratesvation energie[ J]. Geotechnique ,1970( 1)1~ 16.[ 16 ]Christensen Richard W ,Wu Tien Hsing. Analysis of clay defor-6結(jié)論mation as a rate process[ J ]. Journal of Soil Mechanics andFoundations ,1964 ( SM6 ):125~ 157.本文把聚合物網(wǎng)絡(luò)理論和速率過(guò)程理論,應(yīng)用到[ 17袁龍蔚流變力學(xué)M ].北京科學(xué)出版社,1986.粘土流變分析中實(shí)驗(yàn)證明了這種方法的可行性。通[ 18 ]Tobolsky A V , Mark H F. Polymer Science and Material[ M]常認(rèn)為剪切速率越快,強(qiáng)度越高。本文從實(shí)驗(yàn)中發(fā)現(xiàn)中國(guó)科學(xué)院吉林應(yīng)用化學(xué)研究所《聚合物科學(xué)與材料》翻強(qiáng)度隨剪切速率的增大先增后減現(xiàn)象并從理論上進(jìn)譯組譯.北京科學(xué)出版社,1977.行了解釋。也從理論方面為砌塊模型提供了依據(jù)。理. [ 19岡小天近久芳昭高分子物性M]朝倉(cāng)書店,1974.論的積累是前人辛勤工作的成果,本文只是把前人理[ 20 ]Flory P J. Principles of Polymer Chemistry[ M ] New York :Itha-論推廣的初步嘗試。ca ,1953.[ 21 Lodge A s. Elastic I iauidd M 1 NewYork ,1964.[22沈珠中國(guó)煤化工]巖士力學(xué)200 21參考文獻(xiàn):[1]錢家歡殷宗澤.土工原理與計(jì)算(第二版IM]北京:中MHCNMHG[ 23蔣明鏡.結(jié)構(gòu)性粘土的本構(gòu)模型和土體逐漸破損分析[ D ].國(guó)水利水電出版社,1996. 26.南京南京水利科學(xué)研究院,1996.[2]Lambe T Milliam. The structure of compacted clayC[ J] Jourmal[24 ]SDS01- -79 土工試驗(yàn)規(guī)程S]of the Soil Mechanics and Foundations ,1958 { SM2 ):1 ~ 33.