重塑土和原狀土的應力狀態和屈服準則

王杏杏，尚　軍，夏錦紅，李　琳

(1.天津城建大學土木工程學院， 天津300384； 2.新鄉學院土木工程與建筑系， 河南新鄉453003)

?

重塑土和原狀土的應力狀態和屈服準則

王杏杏1，尚軍1，夏錦紅2，李琳1

(1.天津城建大學土木工程學院， 天津300384； 2.新鄉學院土木工程與建筑系， 河南新鄉453003)

摘要：為尋求適用于結構性土應力狀態的合適的表示方法，分析了土的結構性和對應應力狀態之間的關系，認為，土的應力狀態決定于其結構，結構性記憶了前期應力狀態。在此基礎上，研究了重塑土和原狀土的應力狀態和不同加載方式引起的應力狀態改變，分析了總應力、附加應力在描述巖土力學性質方面的客觀性和合理性以及不足。借用加載比和實質應力的概念，研究了簡單條件和復雜條件下土中實質應力計算方法和特征，并建立了基于實質應力表示方法的結構性土的屈服條件。研究表明，與基于一般應力狀態的屈服準則相比，基于實質應力的屈服準則能反映原狀土的既有結構，因此更為合理。

關鍵詞：結構性土；應力狀態；附加應力；實質應力；屈服條件

0引言

采用連續介質力學理論研究不連續介質(巖土材料)的關鍵在于建立適用于該材料的工程本構模型。即使是最簡單的飽和正常固結黏土也具有系列復雜的物理力學特性，這些特性可以概括為三類，即基本特性、亞基本特性和關聯基本特性[1]。

對于飽和巖土，當其礦物成分、顆粒級配確定時，以上諸方面的特性均取決于其物理狀態，并最終取決于材料本身及材料所處的環境。比如，摩擦性依賴于接觸面兩側的材料和施加于接觸面上的正應力[2-3]，壓硬性依賴于土顆粒組成和前期應力狀態，原生各向異性源于土顆粒組成和礦物成分及其沉積過程中的應力狀態[4]，次生各向異性取決于原生各向異性和后期應力狀態[5-6]。對于理想的各向同性材料，其應力狀態的改變不會導致不同方向上力學特性的差異。

傳統觀點認為，附加應力是地基變形的主要原因，總應力是地基強度和穩定性計算的依據。比如，在巖土工程中，以附加應力作為地基沉降計算的依據。在地基基礎設計中，以總應力來對地基的強度和屈服進行校核[7-9]，而較少考慮初始應力狀態及其各向不等性的影響和作用。實際上，初始應力狀態在很大程度上決定了土的基本特性、亞基本特性和關聯基本特性，而總應力和附加應力作用下的強度、穩定和變形均依賴于土的一系列上述特性。因此，合理研究巖土體在總應力和附加應力作用下的力學行為不能簡單地將初始應力狀態從之中減去，而應該綜合考慮當前應力狀態與初始應力狀態的差別、聯系和繼承性。

本文從土的結構性和應力狀態分析入手，在考慮由初始基本特性和應力狀態決定的結構性基礎上，首先，研究了重塑土和原狀土的應力狀態和不同加載方式條件下應力狀態的改變；隨后，基于初始應力誘導的巖土特性，運用實質應力概念研究了初始應力張量與附加應力張量的關系；最后，建立了基于實質應力的巖土屈服條件。

1土的結構和應力狀態

根據成土條件的不同，土可以分為重塑土和原狀土。原始應力指巖體未經受人工擾動前的天然狀態應力，即在施加所考慮的荷載之前土體中已存在的應力。可見，初始應力是一個相對的概念，也即在人為施加荷載作用之前已經處于穩定狀態的土體應力狀態。因此，初始應力狀態引起的土體變形、壓密或固結過程已經完成，且有效應力等于總應力。

在長期的地質作用過程中，顆粒之間形成了不同的膠結作用和排列方式，即形成了土的結構性。巖體的初始應力主要是由自重和地質構造運動引起的，故巖體的地質構造應力是與巖體的裂隙發育、彈塑粘性等特性密切相關的。因此，巖體中的初始應力隨位置不同而變化。此外，影響巖體初始應力狀態的次要因素還有地形、地震力、水壓力、熱應力等。

目前對強度和本構關系的研究尚難以定量化考慮初始應力狀態各向異性對后期強度和變形特性的影響[10-11]。原狀K0試樣經過固結之后再進行三軸試驗，與天然沉積土K0固結狀態的實際應力狀態及受力后的應力路徑不符。

對于自然沉積的成層土質地基，自重應力是不會引起剪應變的，而水平荷載常常引起剪應力，并促使土顆粒體的運動調整并穩定在一個特有狀態。初始應力不同，則土顆粒體的排列、顆粒體間作用、連接的宏觀表現就不一樣。

巖土體經開挖以后將改變周圍相關區域的應力狀態，從而引起變形甚至破壞，且變形和破壞依賴于前期狀態[12-13]。所以，研究巖土體的初始結構和應力狀態，對正確認識巖土體開挖過程中內部應力變化規律，以及合理設計巖體工程和科學施工具有極其重要的意義。

海拔是指某點相對于平均海平面的垂直高度，是一個相對的概念。某點應力的作用也應該與其基點或零點的選擇有關，其原因就在于巖土體當前的微觀結構取決于這個零點或基點。附加應力的基點就是初始應力水平。

土力學中的基底附加壓力是上部結構和基礎傳到基底的地基反力與基底處原先存在于土中的自重應力之差。附加應力是荷載在地基內引起的應力增量，傳統觀點認為它是地基變形的唯一原因。由于天然土層在自重作用下的變形已經完成，故只有超出基底處原有自重應力的那部分應力才使地基產生新的附加變形。將傳統附加應力在三維條件下進行推廣，可以得到廣義附加應力。廣義附加應力等于總應力與初始應力之差，三者之間的關系為：

(1)

其中：

(2)

(3)

(4)

式中，i、j為1、2、3。對于正常固結的K0應力狀態，初始狀態為：

(5)

可見，在采用分層總和法計算地基沉降過程中涉及到的附加應力只是式(1)中的一項，即：

(6)

圖1　大主應力與沉積面成不同角度時土的內摩擦角Fig.1　Soil internal friction angle ofmaximum principal stress anddeposition surface at different angles

對由撒砂法得到的試樣，Matsuoka等取不同方向為大主應力方向進行了常規三軸試驗。沉積面與大主應力作用面間的夾角θ與土的強度參數存在密切關系，如圖1所示[14]。

由圖1可見，各向異性土的抗剪強度隨θ的不同而變化。由于撒砂法得到的試樣各方向上的壓縮模量不同，故在后期等向壓縮條件下各方向的變形也是不同的。因此，土在各方向壓縮模量的不同也可理解為各方向已壓縮程度的不同。顯然，附加應力的概念無法解釋這種由初始應力狀態引起的各向異性。

如前所述，在K0固結狀態下，土體表現為固有各向異性，或初始狀態各向異性。微觀結構上表現為顆粒體的定向分布；力學性質上則表現為不同方向上施加相同的作用，其后果具有顯著差別。這可以理解為初始狀態各向異性能導致應力誘發各向異性，其原因主要在于K0狀態在幾何上表現為顆粒體的定向分布，并由此引起電化學性和力學性質的定向。

2簡單應力條件下的實質應力

這里所說的簡單應力條件是指式(1)中只存在正應力的過程。K0加載、常規三軸實驗和真三軸實驗均屬于此種情況。在簡單應力條件下，初始應力和當前應力可以表示為：

(7)

(8)

任何土類都會在初始應力作用下產生相應的變形，并在后期荷載作用下表現出記憶性。比如，一維壓縮條件下，后期荷載如果加載在前期的壓縮方向，則前期荷載作用下形成的既有密實度和結構對后期荷載的變形影響甚大，相反，如果后期荷載作用在垂直前期壓縮方向上，則會由于該方向上前期的較小壓力而發生較大變形。所以，實質應力的構造應滿足兩個條件，一是能反映在前期荷載作用下已經形成的結構性，但其影響隨著后期荷載的增大逐漸減弱。借用加載比的概念[15]，即加載比等于真實應力與初始應力的比值，則：

(9)

式中，κ*為某一真實應力，可以是正應力、剪應力，也可以是平均應力、偏應力、八面體應力等；κ0為某一初始應力，可以是正應力、剪應力，也可以是平均應力、偏應力、八面體應力等。

某方向上的實質應力σs定義為該方向上的加卸載度除以各方向上加卸載度絕對值之和所得比值后再乘以各加載分量絕對值之和，即：

(10)

(11)

(12)

(13)

如果各方向均為加載過程或均為卸載過程，則：

(14)

可見，實質應力不但可以表示應力絕對值的大小，還考慮了前期作用和其他應力分量的影響。

圖2　重塑土樣和原狀土樣的應力狀態對比　Fig.2　Comparison remodeling soil samples and　undisturbed soil samples stress state

比如，某重塑三軸試樣在圍壓100 kPa下固結完成(初始應力狀態100 kPa，100 kPa，100 kPa)，若逐漸施加0～200 kPa圍壓，則名義應力曲線為l1，3個方向上的附加應力曲線為l2。若某一原狀試樣(初始應力狀態200 kPa，100 kPa，100 kPa)，同樣逐漸施加圍壓0～200 kPa，則豎直方向的附加應力曲線為l3；水平方向的附加應力曲線為l2。對應重塑試樣，由于初始應力狀態各向相等，所以3個方向的實質應力均與其附加應力一致，即曲線l2；而對于原狀試樣，由于初始應力狀態各向異性，在施加圍壓后，豎向實質應力曲線為l4，水平方向的實質應力曲線為l5，如圖2所示。

根據以上各式，可以獲得與初始應力和當前應力對應的實質應力。可見，實質應力不但與該方向上的當前應力和初始應力有關，還依賴于其他方向上的初始應力和當前應力，即能夠反映其他方向應力在當前方向的影響，這恰恰是巖土區別于其他材料的一個顯著力學特點。

(15)

(16)

(17)

3復雜應力條件下的實質應力

復雜應力條件是指式(6)中的正應力分量和剪應力分量可取任意實數且可以單獨變化的過程。嚴格說來，簡單應力條件是工程的一種理想化和簡化，而復雜應力條件則是普遍存在的。在復雜應力條件下，初始應力和現有應力是一個六維向量，即：

(18)

(19)

復雜應力條件下的實質正應力與式(14)具有相同的形式。同樣，為了考慮初始剪應力條件土的穩定變形對后續剪應力的影響，實質剪應力定義為：

為了確保建設項目水土保持“三同時”制度的有效落實，使水土保持設施與主體工程一并招 （發）標、監測、監理同步實施，要建立發改、經信、住建、規劃、環保、國土等相關職能部門聯合監管機制。相關職能部門在建設項目規劃選址、土地預審、環評、立項、核準、備案、開工許可、財政評審、房產預售、竣工驗收、權屬頒證、檔案管理等環節實行一站式服務、并聯制審批、高效率運轉，并在事前向建設單位履行好書面告之義務，確保建設業主的知曉權，有利于業主在開工前能夠順利完成水土保持方案審批手續和水土保持方案后續設計文件，保障建設工程能如期順利推進，合法開工建設。

圖3　實質應力在p-q平面內的定義Fig.3　The definition of essencestress on p-q plane

(20)

其中：

(21)

4基于實質應力的屈服準則

在巖土工程中，屈服準則多是基于重塑土建立起來的，將其直接應用將會帶來較大的系統誤差。如果考慮前期初始應力及其對土結構的影響，理論上是可以消除這種影響的。

三軸實驗數據在π平面的位置為：

(22)

(23)

在-30°≤θσ≤30°范圍內，Mohr-Coulomb準則在π平面上的直線方程為：

(24)

圖4　結構性土的等向壓縮實驗結果Fig.4　Test results of structural soilisotropic compression

Motohisa由撒砂法制得沉積方向為水平的試樣，對其進行等向壓縮，試驗結果如圖4所示[16]。可見，沉積方向的應變始終小于另一正交方向，且其差值Δε基本保持不變，即：

Δε=εx-εz。

(25)

初始狀態為K0固結狀態，相應的初始應力為{K0σz，K0σz，σz}；等p壓縮對應的當前應力狀態為{p，p，p}；等p壓縮對應的附加應力為{p-K0σz，p-K0σz，p-σz}；根據本文算法施加的實質應力為：

(26)

(27)

(28)

(29)

(30)

所以

(31)

因此，可以得到

(32)

當σz?p，即施加的等壓固結荷載很小時，有：

(33)

當p?σz，即對應于重塑土的等壓固結條件時，有：

(34)

且有K0=1，所以

(35)

所以，直線εx=εz是(εz，εx)的漸近線。這一點與三向等壓固結條件下3個方向的變形相等，即εx=εy=εz是一致的。

5結語

顆粒碎散材料的力學性質和工程性質不但與當前環境有關，還取決于其狀態并最終決定于初始應力狀態。為了建立適用于結構性土的應力狀態表示方法，分析了土的總應力、附加應力和實質應力。研究表明，土的應力狀態決定于其結構，結構性記憶了前期應力狀態。在此基礎上，研究了重塑土和原狀土的應力狀態和不同加載方式引起的應力狀態改變，研究了簡單條件和復雜條件下土中實質應力計算方法和特征，并建立了基于實質應力表示方法的結構性土屈服條件。結果表明，與基于一般應力狀態屈服準則相比，基于實質應力的屈服準則能反映原狀土的既有結構，因此更為合理。

參考文獻：

[1]姚仰平，張丙印，朱俊高.土的基本特性本構關系及數值模擬研究綜述[J]. 土木工程學報，2012，45(3)：127-150.

[2]湯連生，張鵬程，王洋，等.土體內外摩擦及摩擦強度試驗研究[J]. 巖石力學與工程學報，2004，23(6)：974-979.

[3]胡學永.粉體材料動摩擦系數的實驗研究[D]. 北京：北京化工大學，2013.

[4]周建，劉正義，嚴佳佳.原生和次生各向異性對軟黏土強度和變形特性影響研究[J]. 巖土工程學報，2013，35(S2)：666-670.

[5]葉朝良，朱永全，劉堯軍，等.原狀黃土各向異性及卸載變形特征試驗研究[J]. 中國鐵道科學，2014，35(6)：1-6.

[6]張坤勇，殷宗澤，梅國雄.土體兩種各向異性的區別與聯系[J]. 巖石力學與工程學報，2005，24(9)：1599-1604.

[7]俞茂宏，彭一江.強度理論百年總結[J]. 力學進展，2004，34(4)：529-560.

[8]陳宏東，都華，馮林平，等.深挖填土基坑監測及安全性分析[J]. 廣西大學學報(自然科學版)，2010，35(2)：204-209.

[9]MINGFENG L， LIMIN P， CHENGHUA S.Calculation of the surrounding rock pressure on a shallow buried tunnel using linear and nonlinear failure criteria[J]. Automation in Construction，2014，37：191-195.

[10]陳志輝，程曉輝.飽和黏土不排水抗剪強度各向異性的熱力學本構模型研究[J]. 巖土工程學報，2014，36(5)：836-846.

[11]陶帥，王學濱，潘一山，等.基于摩爾-庫侖模型的非線性本構模型的開發及其在應變局部化中的應用[J]. 巖土力學，2011，32(S1)：403-407.

[12]侯公羽，劉宏偉，李晶晶，等.基于開挖卸荷效應的地鐵隧道施工過程數值分析[J]. 巖石力學與工程學報，2013，32(S1)：915-2924.

[13]嚴鵬，盧文波，陳明，等.深部巖體開挖方式對損傷區影響的試驗研究[J]. 巖石力學與工程學報，2011，30(6)：1097-1106.

[14]MATSUOKA H，JUNICHI H，KIYOSHI H.Deformation and failure of anisotropic sand deposit[J]. Soil Mechanics and Foundation Engineering，1984，32(11)：31-36.

[15]李順群，張來棟，夏錦紅，等.土的應力狀態和基于mohr-coulomb準則的結構性模型[J]. 中國公路學報，2014，27(11)：1-10.

[16]MOTOHISA H.Drained deformation-strength characteristics of loose shirasu (volcanic sandy soil) under three dimensional stresses[J]. Soils and Foundations，1985，25(1)：65-76.

(責任編輯唐漢民裴潤梅)

Stress state and yield criterion of remolded soil and undisturbed soil

WANG Xing-xing1, SHANG Jun2, XIA Jin-hong2, LI Lin1

(1. School of Civil Engineering, Tianjin Chengjian University, Tianjin 300384, China;

2. Department of civil Engineering and Architecture, Zinxiang University, Xinxiang 453003, China)

Abstract：In order to find a suitable method that could describe the stress state of structural soil, the relationship between structure and corresponding stress state of soil was analyzed. The research shows that the stress state of soil determines its structure, while structure of soil memorizes former stress state. On the basis, stress state caused by different loadings on remolded soil and undisturbed soil were studied for the objectivity, rationality and disadvantages of applying total stress and subsidiary stress to describe geotechnical mechanical properties. Referencing the concepts of load ratio and essence stress, the calculation method and characteristics of essence stress under simple and complex conditions were researched, and the yield condition of structural soil was established based on the representation method of real stress. The research shows that, compared to the yield criterion accounting on the common stress state, the yield criterion accounting on essence stress is more reasonable, because it could describe the present structure.

Key words：structural soils，stress state，subsidiary stress，essence stress，yield criterion

中圖分類號：TU411.7；TU431

文獻標識碼：A

文章編號：1001-7445(2016)01-0196-07

doi:10.13624/j.cnki.issn.1001-7445.2016.0196

通訊作者：尚軍(1968—)，男，天津人，天津城建大學副教授；E-mail: ttsjll@126.com。

基金項目：國家自然科學基金資助項目(41172233，41472253)

收稿日期：2015-05-15；

修訂日期：2015-12-26

引文格式：王杏杏，尚軍，賈紅晶，等.重塑土和原狀土的應力狀態和屈服準則[J].廣西大學學報(自然科學版)，2016，41(1)：196-202.