基于t 準則的各向異性強度準則及變換應(yīng)力法1)

萬征 孟達

(中國建筑科學(xué)研究院地基基礎(chǔ)研究所,北京 100013)

引言

自然界中的很多天然材料具備強烈的各向異性性質(zhì),比如巖石、木材、土壤等材料在宏觀尺度的各個方向上表現(xiàn)出較大的差異性.顯然,從變形以及破壞機理角度解釋,則是由于材料在細觀層次上具有顯著的差異性所導(dǎo)致的.Oda 等[1]的試驗結(jié)果證實:對于某些條件下,巖土材料的各向異性則是由于微觀顆粒在自然沉積作用過程中顆粒的空間排列定向性以及土壤團粒膠結(jié)過程中的復(fù)雜作用而形成的結(jié)果.一般而言,對于層狀水平分布的巖土材料,由于在水平方向內(nèi)顆粒間的隨機分布狀態(tài),顆粒長軸一般會平行于水平沉積面,因而形成了正交各向異性,也可稱之為橫觀各向同性性質(zhì).

Abelev 等[2]針對大主應(yīng)力與沉積面呈不同夾角下的試樣開展了真三軸排水加載測試,結(jié)果表明:大主應(yīng)力方向與沉積面法向一致的試樣表現(xiàn)出了更高的應(yīng)力比強度,而當(dāng)大主應(yīng)力方向與沉積面法向相互垂直時,則表現(xiàn)出了較低的應(yīng)力比強度.Kirkgard等[3]針對舊金山灣區(qū)黏土的試驗結(jié)果也得到了相同的結(jié)論.Duncan 等[4]首先發(fā)現(xiàn)飽和自然沉積黏土在不排水剪切加載測試中隨著大主應(yīng)力與沉積面夾角的變化會產(chǎn)生差異顯著的應(yīng)力應(yīng)變關(guān)系結(jié)果.Yong等[5]針對靈敏性黏土開展的無側(cè)限壓縮強度測試表明,最小強度值僅為最高強度值的60%.Nishimura等[6]發(fā)現(xiàn)自然沉積黏土的強度具有很強的沉積方向依賴性.同樣,Yamada 等[7],Ochiai 等[8],Miura 等[9],Hight 等[10],Tatsuoka 等[11]以及Pradhan 等[12]針對砂土所開展的一系列三軸壓縮、真三軸加載以及空心圓柱扭剪等復(fù)雜加載測試結(jié)果表明:自然沉積砂土強度以及應(yīng)力應(yīng)變關(guān)系特性同樣具有強烈的沉積方向依賴性.為描述橫觀各向同性性質(zhì),從細觀機理上揭示與描述上述各向異性性質(zhì)強弱程度,Oda 等[13]建議采用組構(gòu)張量來描述微觀顆粒形狀信息的長細比以及長軸在空間中分布方向等宏觀統(tǒng)計信息.為了能夠考慮各向異性性質(zhì)對于砂土應(yīng)力應(yīng)變關(guān)系以及強度特性的影響,Li 等[14]和Dafalias 等[15]建議將一般應(yīng)力量σi j與組構(gòu)張量Fij進行并乘以得到一組反映各向異性性質(zhì)的聯(lián)合應(yīng)力不變量,使用上述聯(lián)合應(yīng)力不變量來構(gòu)建彈塑性本構(gòu)模型.Pietruszczak 等[16]和Mroz 等[17]建議了一個隨組構(gòu)量變化的黏聚力和內(nèi)摩擦角的強度指標,以此反映各向異性對于臨界狀態(tài)的影響.Hashiguchi 等[18]為了能夠反映原生各向異性性質(zhì),將組構(gòu)張量中的三個主分量直接弓入到本構(gòu)模型中屈服面的轉(zhuǎn)軸分量中,通過屈服面在應(yīng)力空間中初始位置的不同來反映各向異性程度.

國內(nèi)學(xué)者對于各向異性的研究一直處于熱點狀態(tài).張連衛(wèi)等[19]基于SMP 準則,將構(gòu)成SMP 面的三個平面的摩擦角表達為該平面隨沉積面夾角變化的變量,構(gòu)造了反映沉積面信息的ASMP 準則.曹威等[20]在SMP 準則基礎(chǔ)上,將組構(gòu)張量參量弓入到摩擦法則表達式中,以此能夠反映橫觀各向同性砂土的各向異性性質(zhì).姚仰平等[21]采用SMP 空間滑移面與沉積面之間的夾角作為基本變量,構(gòu)造了反映橫觀各向同性材料強度的應(yīng)力比公式.Kong 等[22]采用考慮了巖土組構(gòu)張量修正的加載應(yīng)力量,并與各向同性SMP 準則相結(jié)合,得到了考慮微觀結(jié)構(gòu)影響的橫觀各向同性強度準則.路德春等[23]采用微觀結(jié)構(gòu)張量在滑動面法線方向上的投影,定義了一個反映三維橫觀各向同性強度參數(shù),并利用此參數(shù)來修正SMP 準則,據(jù)此建立了一個反映橫觀各向同性土體的強度準則.劉洋等[24-26]也采用反映微觀信息的組構(gòu)張量與應(yīng)力相結(jié)合的方法對各向異性性質(zhì)進行了研究.此外,黃茂松等[27]從微觀機理出發(fā),也對各向異性性質(zhì)進行了探討.王國盛等[28]針對混凝土的加載速率效應(yīng)提出了S 強度準則.高江平等[29]考慮菱形十二面單元體主剪面上三個主剪應(yīng)力與三個正應(yīng)力都會對材料破壞產(chǎn)生影響,建議了三剪應(yīng)力統(tǒng)一強度準則.

上述已提出的各類各向異性準則,按照構(gòu)建的思路來分類,可大致分為如下幾類:(1)將大主應(yīng)力與沉積面夾角作為反映各向異性程度的變量來對已有的各向同性強度準則進行修正; (2)根據(jù)組構(gòu)張量來構(gòu)造各向異性狀態(tài)變量,用來修正既有各向同性準則;(3)利用組構(gòu)張量與應(yīng)力量進行并乘得到聯(lián)合應(yīng)力不變量來修正各向同性準則;(4)基于SMP 準則或者某一強度準則,用物理破壞面與沉積面夾角作為狀態(tài)變量來修正各向同性準則;(5)擴展經(jīng)典各向同性準則為各向異性強度準則.

上述各種方法各有利弊,采用大主應(yīng)力與沉積面夾角來作為判斷各向異性程度發(fā)揮的變量,顯然比較直觀,然而在某些情況下,比如平面應(yīng)變條件下,各向異性強度值與上述夾角并不呈現(xiàn)單調(diào)關(guān)系,而是形成“V” 字形曲線,顯然若采用上述加載角與沉積面夾角作為各向異性狀態(tài)變量存在弊端.利用組構(gòu)張量表達的各向異性程度狀態(tài)變量來直接修正既有的強度準則,可能具有一定的實用性,但仍然不具有普適性.利用組構(gòu)張量與普通應(yīng)力量按照某種規(guī)則相結(jié)合為聯(lián)合應(yīng)力不變量的方法能夠在一定程度上考慮了各向異性對于應(yīng)力量的修正作用,但聯(lián)合應(yīng)力不變量的構(gòu)造方法目前只是猜想階段,缺乏嚴謹?shù)睦碚撟C明.采用基于某種準則的破壞面與沉積面的夾角作為狀態(tài)變量,由于破壞面與沉積面都具有明確的物理概念,且夾角與各向異性強度呈現(xiàn)單調(diào)關(guān)系,表明具有較好的歸一性.但在某些路徑下,如三軸壓縮或三軸伸長加載情況下,會出現(xiàn)較多的破壞面,需要進行篩選,選擇其中最小強度值,這可能會違背物質(zhì)客觀性原理.利用經(jīng)典準則擴展為各向異性準則,只是對原有準則在一定程度上的修補,缺乏物理機理以及普適性.

基于已提出的適用于具有凝聚力以及摩擦性材料的t 強度準則[30-34],在有效滑移面的物理概念上,考慮沉積面的位置關(guān)系,利用有效滑移面與材料沉積面之間夾角作為反映各向異性程度的狀態(tài)變量.由于t 準則在偏平面上是能夠反映包括金屬材料的Von-Mises 準則到巖土材料的SMP 準則,在子午面上是冪函數(shù)表達式,因而其適用范圍大,且具有明確的物理涵義.按照上述各向異性狀態(tài)變量所構(gòu)建的各向異性強度準則,能夠反映諸多材料的各向異性性質(zhì),如金屬、巖石、混凝土、黏土、砂土等.通過上述各向異性強度準則公式,在主應(yīng)力空間中按照將各向異性t 準則變換為Von-Mises 準則的變換思路,推導(dǎo)了變換應(yīng)力公式,該變換應(yīng)力公式實質(zhì)上是從各向異性應(yīng)力空間到各向同性應(yīng)力空間的轉(zhuǎn)換方程.

1 t 強度準則

根據(jù)已提出的t 強度準則,如圖1 所示,立方體表示材料單元,當(dāng)材料發(fā)生破壞時,按照應(yīng)力比的思想,則假設(shè)存在空間有效滑移面ABC,在ABC面上等效主剪應(yīng)力τen與主法向力σen之比成為一個材料常數(shù).其中φe12表示由第一、第二主應(yīng)力所構(gòu)成的有效摩擦角,而φe23表示由第二、第三主應(yīng)力所構(gòu)成的有效摩擦角.

圖1 空間坐標系中的滑移面Fig.1 The slip plane in three dimensional space coordinate

令EA=1,根據(jù)三角函數(shù)關(guān)系,可推導(dǎo)得到如下關(guān)系:

在直角三角形AEB中,tan(45°-φe13/2)=EA/EB,根據(jù)三角函數(shù)關(guān)系求解上式可得到如下公式

同理亦可得到

可設(shè)置一個表征摩擦性與凝聚性權(quán)重分配的參數(shù)t,且0 ＜t＜1,由此可得

弓入?yún)?shù)t用來反映摩擦力與凝聚力的比例權(quán)重,兩者所對應(yīng)的有效摩擦角可分別由τ–σ 空間內(nèi)的摩爾圓外切直線的反正切值來表示,相應(yīng)的直線截距則可分別表示為

由圖可見,由于確定了∠ABE與∠ACE,則相應(yīng)的空間滑移面為平面的假設(shè)前提下則已完全確定下來,另外的σ1,σ2之間構(gòu)成的夾角則由三角函數(shù)的正切值定義可求出

可解出

根據(jù)式(1)、式(2)的關(guān)系,將其代入可得

式(10)即為空間有效滑移面中用以確定有效滑移面(ESMP)空間位置的三個角的關(guān)系式.

如圖2 所示,若需確定有效滑移面ESMP,則需要先確定該面的法線方向,法線方向可由余弦來表示.

圖2 空間坐標系中的滑移面與沉積面Fig.2 The slip plane and depositional plane in three dimensional space coordinate

空間滑移面的方向余弦可表示為ω1(l1,m1,n1),而巖土材料在空間坐標系中存在一個空間沉積面(DP),此空間沉積面可以用該面的方向向量表示,則可令空間沉積面的方向余弦表示為:ω2(l2,m2,n2),則兩空間平面的夾角可由兩方向向量點積反余弦值表示為

強度可由應(yīng)力比表示,由于強度值隨上述夾角α變化,顯然,當(dāng)α=αmin=0°時,則空間滑移面與空間沉積面相重合,此時由于沉積面之間聯(lián)結(jié)最為薄弱,因而強度最低,當(dāng)α=αmax時,則空間滑移面與空間沉積面呈現(xiàn)最大夾角狀態(tài),此時最難破壞,因而強度最高,而當(dāng)0°＜α ＜αmax時,則應(yīng)力比介于上述兩者之間,由此可見,可利用上述兩者極端應(yīng)力比強度值,選用合理的插值函數(shù)進行表達,根據(jù)上述思路,可利用比較簡單的內(nèi)插函數(shù)來表示各向異性強度表達式

空間沉積面法向單位向量為ω2,法向量在3 個空間坐標方向的投影分別如圖3 所示,ω2與xy平面夾角為α1,其在xy平面的投影向量與x軸夾角為α2.則顯然,空間沉積面法向方向可表示為

圖3 空間坐標系中沉積面方向投影Fig.3 The project vector of depositional plane in three dimensional space coordinate

首先先求出有效空間滑移面與沉積面夾角為90°時的強度值.

令

則有效空間滑移面法向向量分量可表示為

根據(jù)三角函數(shù)關(guān)系,可知

根據(jù)對正四面體AEBC 的力平衡關(guān)系,并利用上述公式,可得到等效正應(yīng)力公式.

則等效正應(yīng)力可表示為

經(jīng)推導(dǎo)可得到如下公式

其中,φmo表示空間有效滑移面的內(nèi)摩擦角.

當(dāng)處于三軸壓縮時,則式(27)可表達為

此時,大小主應(yīng)力分別可表示為

將式(28)代入式(26)中,可得到關(guān)于p,qc的函數(shù)為

采用破壞時應(yīng)力比的表示方法,M=qc/p,則可表示為

其中,M表示三軸壓縮時所對應(yīng)的破壞應(yīng)力比.若三軸伸長下的破壞應(yīng)力比為Me,且設(shè)三軸伸長破壞應(yīng)力比與三軸壓縮破壞應(yīng)力比之比值為λ,則Me=λM

由于在三軸壓縮路徑下,方程式(26)、式(29)完全相等,因此得到

式(35)即為偏平面上廣義偏應(yīng)力強度公式.

子午面上可采用考慮靜水壓力效應(yīng)的關(guān)于平均應(yīng)力的雙曲線函數(shù)來作為強度表達式

2 三維正交各向異性的特征沉積面

由上節(jié)可知,對于橫觀各向同性而言,其三維立方體只存在一個沉積面,對于一般情況而言,存在三維正交各向異性固體材料,其3 個方向都存在物理上的沉積面.如圖4 所示.考慮沿著xyz三方向都存在著沉積面,3 個沉積面兩兩正交,且沿著3 個方向的特征沉積面的沉積厚度分別為hx,hy,hz.則顯然由圖可見

圖4 三維空間坐標系中的三維特征沉積面Fig.4 Three dimension characteristic depositional plane in three dimensional space coordinate

上述3 個方向的沉積面厚度為沉積面的特征尺寸,表征沿該長度方向的性質(zhì)均勻化的最小單位.而同理,在y以及z方向都存在類似的沉積特征厚度.由上述3 個特征沉積厚度,假定3 個正交沉積面同時發(fā)揮抵抗破壞的貢獻影響力,可由此得到空間的三維特征沉積面,空間四面體為LMOP.空間特征沉積面的法向向量可由MO 與MP 的叉積得到.

空間特征沉積面法向單位向量為

當(dāng)巖土材料為二維正交各向異性情況時,此時其中一個方向的特征尺寸顯然為無窮大,此時存在3種情形:

(1)當(dāng)hx=+∞時,則顯然式(42)可退化為如下公式

顯然由式(43)可知,該二維特征沉積面平行于x軸.

(2)當(dāng)hy=+∞時,則顯然式(42)可退化為如下公式

上述沉積面平行于y軸.

(3)當(dāng)hz=+∞時,則顯然式(42)可退化為如下公式

當(dāng)巖土材料退化為一維各向異性,即為橫觀各向同性情形時,此時其沉積面仍分為3 種情形:

(1)當(dāng)hx=+∞,hy=+∞同時滿足,此時式(42)退化為顯然此時法向方向為沿著z軸方向.

(2)當(dāng)hy=+∞,hz=+∞同時滿足,此時式(42)退化為

顯然此時法向方向為沿著x軸方向.

(3)當(dāng)hz=+∞,hx=+∞同時滿足,此時式(42)退化為

顯然此時法向方向為沿著y軸方向.

3 橫觀各向同性t 強度準則

對于各向異性的考慮,需要將式(36) 中應(yīng)力比強度Mf表示為有效滑移面與三維特征沉積面夾角的函數(shù)即可.則顯然需要構(gòu)造一個表示夾角的關(guān)系式.首先需滿足以下兩個條件:

(1) 首先滿足夾角越大,則強度值越大的單調(diào)規(guī)律;

(2)需要構(gòu)造一個各向同性函數(shù),以此滿足物質(zhì)客觀性原理.

其中,右端的φmo表示有效滑移面的內(nèi)摩擦角,在各向異性情形下,通常tanφmo并非恒定值,可由函數(shù)tanφmo=F(α,M)來表示.α 表示滑移面與三維特征沉積面的夾角,各向異性程度的度量參數(shù),因此上述函數(shù)包含了描述各向異性方向以及各向異性程度雙重信息.

同理由有效滑移面與三維特征沉積面夾角公式由兩個面法向向量的點積反余弦來表示,可表示為如下公式

由于在三軸壓縮條件下,可得到式(49) 的簡化表達式為

采用大主應(yīng)力與沉積面之間夾角δ 作為衡量各向異性對強度應(yīng)力比的主要影響變量是慣常的一種做法,上述方法可以部分的考慮各向異性沉積面對于強度的貢獻.然而,由Matsuoka 等[35]關(guān)于砂土在平面應(yīng)變的強度測試結(jié)果表明,平面應(yīng)變強度與δ 是非單調(diào)關(guān)系,兩者關(guān)系出現(xiàn)先減小后增大的規(guī)律特點.但若是采用空間滑動面與沉積面夾角ζ 與強度值建立關(guān)系,則會發(fā)現(xiàn)兩者具有單調(diào)關(guān)系.

通過類比橫觀各向同性沉積面與空間滑移面之間夾角與強度應(yīng)力比的單調(diào)試驗關(guān)系結(jié)果,進一步可將其推廣為三維正交各向異性的情形,此時,三維特征沉積面綜合考慮了3 個主方向的沉積面的影響,在各自3 個主方向上由沉積面特征厚度來確定該方向上的特征尺度,綜合反映在三維特征沉積面上.該三維特征沉積面與空間滑移面之間的夾角與最終強度應(yīng)力比表現(xiàn)為單調(diào)關(guān)系.可將其簡單地由如下拋物線型公式來表達

其中,x 表示max 的簡稱,表示最大值,而n 為min 的縮寫,表示為最小值.分析式(51) 可知,強度應(yīng)力比隨夾角β 呈現(xiàn)二次函數(shù)的單調(diào)遞增關(guān)系,完全符合試驗點的單調(diào)遞增規(guī)律特點.

當(dāng)為橫觀各向同性情形時,大主應(yīng)力方向垂直于沉積面時,此時存在最大的角度βx.由類比可知,當(dāng)大主應(yīng)力方向與三維特征沉積面垂直時,此時存在最大角度βx.而當(dāng)ESMP 面與三維特征沉積面相重合時,則顯然存在最小的夾角βn=0.對于Mx的確定,可由圖5 所示加載工況,根據(jù)三維特征沉積面與大主應(yīng)力軸相垂直條件下的三軸壓縮試驗來確定最終的內(nèi)摩擦角,由摩擦角可通過式(52)計算得到.而Mn則針對的是當(dāng)β=0 時候的強度應(yīng)力比,顯然直接確定較為困難,可間接通過由順沉積面的三軸壓縮試驗得到的強度值M0來求出,由圖6 可見,采用的是大主應(yīng)力方向與沉積面法向相互垂直時候的加載工況,M0為此種加載條件下的強度應(yīng)力比.其中表達式可表示為

強度應(yīng)力比最小值Mn可通過式(52)變形導(dǎo)出,可由通過常規(guī)試驗測試得到的強度值M0以及最大值Mx的關(guān)系表達式表達.由于三軸壓縮條件,因而EC0=1,聯(lián)合式(51)可得到

圖5 大主應(yīng)力與三維特征沉積面垂直時的三軸壓縮示意圖Fig.5 Triaxial compression condition with major principal stress vertical to the depositional plane

圖6 大主應(yīng)力與三維特征沉積面平行時的三軸壓縮示意圖Fig.6 Triaxial compression condition with major principal stress parallel to the depositional plane

其中

考察各向異性性質(zhì)對于三軸壓縮強度的影響規(guī)律,由于微觀顆粒在重力場及外部自然作用下會形成一定的排序分布,在空間中微觀顆粒會產(chǎn)生定向性,若用空間橢球體作為對長方體顆粒的近似,則長軸會平行于沉積方向,而垂直于沉積面的法向方向是沉積面的空間對稱軸.由于長軸方向平行于沉積面的分布形態(tài)是形成巖土體的一種穩(wěn)態(tài)結(jié)構(gòu),因而,在這種作用下形成的巖土體在自然界非常普遍.下面關(guān)于空間沉積面在空間中的排布方式,按照研究方便考慮3 種特殊情況,分別對應(yīng)空間有效滑移面與沉積面的3 個特殊位置.

當(dāng)為橫觀各向同性情況時,此時只沿著一個方向分布沉積面.

(1) 考察式(49),若當(dāng)沉積面法向方向與z軸相一致時,則可知,hx=∞,hy=∞,則滑移面與沉積面夾角表示為

(2)當(dāng)沉積面法向方向與x軸相一致時,則可知,hz=∞,hy=∞,則滑移面與沉積面夾角表示為

(3)當(dāng)沉積面法向方向與y軸相一致時,則可知,hz=∞,hx=∞,則滑移面與沉積面夾角表示為

而當(dāng)處于二維各向異性情況時,此時也存在3 種情況.

(1)當(dāng)沉積面沿x與y軸均存在分布時,則此時可知hz=∞,由此可知滑移面與特征沉積面夾角表示為

(2)當(dāng)沉積面沿y與z軸均存在分布時,則此時可知hx=∞,由此可知滑移面與特征沉積面夾角表示為

(3)當(dāng)沉積面沿z與x軸均存在分布時,則此時可知hy=∞,由此可知滑移面與特征沉積面夾角表示為

綜合考慮偏平面與子午面上強度表達式,聯(lián)立式(26)、式(35)、式(36)、式(51),可得到最終關(guān)于各向異性的巖土非線性強度準則

式(65)即為表示各向異性強度準則表達式.

下面具體討論當(dāng)參量取值不同時強度準則的表現(xiàn)形式.

(1) 對于各向異性強度準則偏平面表示式,式(51) 中當(dāng)Mβ與Mx,M0各不相等時,則表示的是一般的橫觀各向異性強度應(yīng)力比,此時,Mβ為一隨β 角變化的強度值,此時強度準則為三維正交各向異性a準則.

(2) 對于各向異性強度準則偏平面表示式,式(51)中當(dāng)Mβ=Mx=M0時,則由于其兩相互垂直方向的應(yīng)力比強度相同,則介于兩者之間的任意應(yīng)力比強度值都相同,則此時,Mβ為一恒定值.當(dāng)a13=a23時,則此時強度準則退化為a準則.

(3)進一步,當(dāng)a13=a23=0 時,則強度準則退化為偏平面上為SMP 準則的非線性強度準則.

(4)當(dāng)Mx與M0都不為零,而a13=a23=0 時,則強度準則退化為偏平面上為各向異性SMP 準則的非線性強度準則.

4 主要參數(shù)意義及確定方法

4.1 參數(shù)t 對于偏平面上強度準則曲線的影響

由于該強度準則的核心是由比例因子參數(shù)t作為控制偏平面上形狀的主要因素,因此命名該準則為t 準則.由圖7 可見,反映摩擦性與凝聚性權(quán)重的比例因子參數(shù)t對于偏平面上不同應(yīng)力羅德角下材料的強度特性影響顯著.當(dāng)t=0 時,此時,由于有效滑移角為零,因而材料破壞只受到偏差應(yīng)力強度控制,因而退化為廣義Mises 強度準則,反映的是金屬材料的宏觀破壞性質(zhì).當(dāng)t=1 時,則t 準則退化為SMP 準則,反映的是純摩擦性的材料破壞特性,而當(dāng)0 ＜t＜1 時,則反映的是具有部分摩擦部分凝聚性材料的破壞特性.而當(dāng)t＞1 時,則反映的是受應(yīng)力羅德角影響更為顯著的材料破壞特性,如圖中t=2 時的破壞曲線,反映了偏平面上破壞曲線逐漸趨向于等三角形曲線的特點.

圖7 各向同性不同t 值下偏平面上的破壞曲面Fig.7 Failure surfaces with different a value at deviatoric plane in isotropic stress space

由圖7 可見,顯然對于任意t值下,在三軸壓縮所對應(yīng)的偏應(yīng)力強度值恒定,但對于其他應(yīng)力羅德角下的廣義偏應(yīng)力強度值,卻均不同,因而,顯然可將t值與描述偏平面上曲線形狀特性相聯(lián)系.觀察圖7 可知,由于對應(yīng)三軸伸長的偏應(yīng)力強度值各個不同,因而,可將三軸伸長偏應(yīng)力強度值與對應(yīng)的三軸壓縮偏應(yīng)力強度值的比值λ=qc/qe,作為確定t值參數(shù)的依據(jù).根據(jù)近似,建議可取如下公式作為確定參數(shù)t的表達式

分析上式,顯然當(dāng)λ=1 時,則此時t=0,對應(yīng)的是廣義Mises 準則,而λ=3/(3+M)時,則此時t=1,對應(yīng)的是SMP 準則.

圖8 所示分別為不同t值下的破壞面形態(tài),為了便于觀察,取Mf不同值時,即對Mf由小到大分別為0.7,1.0,1.4,1.8,2.4 下的破壞面分別對應(yīng)著圖8 中由內(nèi)到外的空間破壞面,其所對應(yīng)的t值分別為0,0.3,0.6,0.8,1.0.由圖可見,當(dāng)t值由小到大變化時,對應(yīng)的破壞面的形態(tài)在偏平面上由圓形逐漸過渡為較為尖銳的曲邊三角形形態(tài).

圖8 主應(yīng)力空間中參數(shù)t 及Mf 不同值下的三維空間破壞面Fig.8 Failure surfaces in 3d space under different values of parameters t and Mf in the principal stress space

對于各向同性a準則,當(dāng)參數(shù)a變化時,由圖9可見,隨著a值由0～1 逐漸增大,偏平面上強度曲線由表示SMP 準則的曲邊三角形逐步過渡到表示Von-Mises 準則的圓形曲線.由此可見,參數(shù)a關(guān)于中主應(yīng)力對于強度應(yīng)力比的影響非常明確直接,可以直接決定中主應(yīng)力對于強度的貢獻程度.

4.2 參數(shù)n 意義及確定

圖9 大主應(yīng)力與沉積面垂直時的三軸壓縮及三軸伸長示意圖Fig.9 Triaxial compression and triaxial extension conditions with major principal stress vertical to the depositional plane

圖10 不同冪參數(shù)n 值下破壞曲線形態(tài)Fig.10 Failure curves with different n value

從圖10 可看出,在子午面上,由下到上分別對應(yīng)不同n值下的破壞曲線,分別對應(yīng)n=0.2,0.4,0.6,0.8,1.0 時的破壞曲線.由此可看出,參數(shù)n對于破壞形態(tài)的影響主要體現(xiàn)在兩個方面.第一,在子午面上,隨著n值增大,則子午面上破壞線逐漸趨近于一條斜直線,隨著n值減小逐漸趨近于零,則破壞線逐漸趨近于平行于球應(yīng)力軸的水平直線,當(dāng)為水平直線時,則破壞時對應(yīng)的廣義偏應(yīng)力恒定不變,即不受靜水壓力影響,起決定作用的是廣義偏應(yīng)力.而當(dāng)n值介于0 與1 之間時,則子午面上對應(yīng)的是一條冪函數(shù)曲線,隨著球應(yīng)力的增大,則廣義偏應(yīng)力也逐漸增大,但增大率逐漸減小,體現(xiàn)了非線性特點,表明靜水壓力對于三軸壓縮下的剪切強度有直接影響.第二,在偏平面上,隨著n值的增大,則三軸壓縮對應(yīng)的廣義偏應(yīng)力強度相應(yīng)增大,在相同的球應(yīng)力下,隨著三軸壓縮剪切強度的增大,則不同羅德角作用下的強度值差異逐漸顯著,也就是應(yīng)力誘導(dǎo)各向異性程度逐漸加深.

通過式(36)變形可得到如下的關(guān)系式

以等式左右兩端的對數(shù)分別為變量,則上述公式顯然是關(guān)于參數(shù)n的一次函數(shù),利用三軸壓縮試驗結(jié)果,可將對應(yīng)不同靜水壓力下的剪切強度整理在對數(shù)坐標系內(nèi),擬合出的直線型曲線,斜率為n,而截距值為lnMf.

4.3 參數(shù)σ 涵義及確定

σ0為強度曲線與靜水壓力軸的左交點值,其物理意義表示材料在拉伸條件下的強度值,可以反映出材料的凝聚力特性.在實際狀態(tài)中三向拉伸作用下的強度值,一般很難實現(xiàn),對于無黏性土,則取為0.對于具有拉伸強度的材料,如混凝土等,根據(jù)過鎮(zhèn)海的建議,可將其取為單軸拉伸強度值的0.9 倍.

4.4 參數(shù)pr 涵義及確定

參數(shù)pr反映在一定靜水壓力下,將剪切強度q歸一化的特征壓力.另外,也起到將靜水壓力無量綱化的作用,通常在砂土等散粒體條件下,可取為一個大氣壓力值.參數(shù)pr的確定可根據(jù)式(36),由對試驗結(jié)果整理在對數(shù)坐標系內(nèi),根據(jù)擬合出來的直線來確定參數(shù)pr的值.

4.5 參數(shù)Mx 與M0 涵義及確定

作為反映各向異性的參數(shù),通常只需確定沉積面與空間有效滑移面的呈現(xiàn)最大夾角和最小夾角狀態(tài)下的三軸壓縮強度值,即可得到相應(yīng)的Mx與Mn值.Mx根據(jù)大主應(yīng)力方向垂直于沉積面加載條件下的三軸壓縮強度值來測試得到,而Mn可間接通過M0用式(53)求取得到,而M0則是利用大主應(yīng)力與沉積面順層方向的三軸壓縮強度測試得到.

5 各向異性對強度曲線影響

由圖3 可知,當(dāng)空間沉積面的法線不與三維空間坐標軸的任一軸重合時,即空間沉積面處于三維空間的一般位置時,則沉積面與有效滑移面的夾角可由式(49)來表述,通過式(51)可知,應(yīng)力比強度值是夾角β 的單增函數(shù).當(dāng)沉積面空間位置一旦確定后則是固定不變的,而隨著主應(yīng)力大小以及方向的調(diào)整,顯然有效滑移面始終是動態(tài)變化的,在非三軸壓縮條件下,則其強度值始終是小于Mx,隨著夾角β在零度附近變化,則強度值也會相應(yīng)的升高或降低,若按照大主應(yīng)力與沉積面法向夾角為加載角的定義,則強度值會隨著加載角產(chǎn)生先減小后增大的規(guī)律.

如圖11 所示,為便于描述,按照主應(yīng)力的變化在偏平面上將全部區(qū)域分為3 個象限I,II,III.

圖11 描述橫觀各向同性狀態(tài)的偏平面空間區(qū)劃圖Fig.11 Zones describing transverse isotropic state in deviatoric plane

(1)當(dāng)沉積面法線方向與x軸,即σ1方向相一致時,則顯然,由圖12 可見,當(dāng)為各向同性強度準則時,如圖12 中的虛線與點劃線,即當(dāng)Mf為恒定值時,則破壞曲線在偏平面上為一個幾何對稱圖形,對于實線,及Mf為一個變量時,則此時的破壞曲線僅僅關(guān)于材料主軸對稱.

由圖13 可見,當(dāng)沉積面垂直于y軸時,則此時最大應(yīng)力比為位于與σ2軸相重合的縱軸上,而兩側(cè)的幾何破壞線關(guān)于該軸幾何對稱分布.

由圖14 可見,當(dāng)沉積面垂直于z軸時,則此時最大應(yīng)力比為位于與σ3軸相重合的軸線上,而兩側(cè)的幾何破壞線關(guān)于該軸幾何對稱分布.

圖12 各向異性強度準則(沉積面垂直于x 軸)Fig.12 Failure surface of anisotropic criterion(depositional plane vertical to the x axis)

圖13 各向異性強度準則(沉積面垂直于y 軸)Fig.13 Failure surface of anisotropic criterion(depositional plane vertical to the y axis)

圖14 各向異性強度準則(沉積面垂直于z 軸)Fig.14 Failure surface of anisotropic criterion(depositional plane vertical to the z axis)

綜合圖12～圖14 可見,采用考慮有效滑移面與沉積面的夾角作為反映各向異性程度的狀態(tài)量,能夠有效地反映出材料本身沉積面對于強度應(yīng)力比的影響規(guī)律.

6 基于橫觀各向同性t 準則的變換應(yīng)力公式

為了可簡單合理的將已有的基于p-q空間的二維彈塑性模型推廣到能應(yīng)用于描述原生各向異性特性對本構(gòu)關(guān)系的影響規(guī)律,基于所提的橫觀各向同性a準則,將其應(yīng)用于構(gòu)造變換應(yīng)力方程[36-38],由此可將橫觀各向同性應(yīng)力空間轉(zhuǎn)變?yōu)楦飨蛲詰?yīng)力空間.其基本思路是,先將一般的考慮各向異性的廣義偏應(yīng)力q進行歸一化處理,即將對于某一特定加載條件下的q轉(zhuǎn)換為qf,即第一步實現(xiàn)了將各向異性轉(zhuǎn)換為各向同性,qf即為對應(yīng)大主應(yīng)力方向與沉積面法向方向一致的三軸壓縮的廣義偏應(yīng)力強度值.然后再通過將各個偏應(yīng)力分量放大到對應(yīng)von-Mises圓的對應(yīng)點上,第二步是將各向同性應(yīng)力空間中主應(yīng)力對應(yīng)力比強度的影響歸一化處理.通過上述兩步驟,即可實現(xiàn)由各向異性應(yīng)力空間與各向同性應(yīng)力空間的轉(zhuǎn)換.

由Satake 等針對純摩擦特性所提出的基于SMP準則的形狀函數(shù),基于摩擦角的概念而提出,因而本文中的有效摩擦角可將其替換以得到基于有效摩擦角的形狀函數(shù)

EB,EC是σ1,σ2,σ3的函數(shù),具體見式(6)、式(7)、式(26).

其中,θ 為應(yīng)力羅德角,可表示為

對應(yīng)a準則的偏平面上的形狀函數(shù)可表示為

由于已知有基于a準則的形狀函數(shù),因此,可得到在對應(yīng)任意一個球應(yīng)力p下的三軸壓縮路徑下的廣義偏應(yīng)力qc,表示為

考慮到偏平面上各向異性的修正,最終得到偏平面上考慮各向異性的a準則

基于各向異性強度準則的轉(zhuǎn)換應(yīng)力公式可表達為

將其推廣為一般應(yīng)力表示的轉(zhuǎn)換應(yīng)力公式

式(77)中?qac/?σi可表示為

式(55)～式(63)為將普通應(yīng)力轉(zhuǎn)換為變換應(yīng)力空間的變換應(yīng)力公式,而式(64)～式(88)為將變換應(yīng)力空間中變換應(yīng)力應(yīng)用到具體本構(gòu)模型中時微分的導(dǎo)函數(shù)公式.

7 平面應(yīng)變條件下的t 準則公式

平面應(yīng)變條件是自然界以及工程實踐中常常遇到的某種約束條件,指的是三維條件下材料某一方向的尺寸遠大于另兩個相近的方向上的尺寸,比如堤壩、路基或者擋土墻等工程中的土體約束作用.將所提t(yī) 準則應(yīng)用于平面應(yīng)變約束下的工況,一方面可將其具體應(yīng)用化,另一方面可對于所提準則進行相應(yīng)的檢驗.

假設(shè)屈服準則表達式為

且屈服準則為理想彈塑性屈服面,則對中主應(yīng)變的塑性增量,可由一致性條件表達為

由于理想彈塑性的流動法則可知,當(dāng)進入塑性流動狀態(tài)后,應(yīng)變增量完全為塑性應(yīng)變增量,此時,彈性應(yīng)變增量為零,而由平面應(yīng)變條件可知,中主應(yīng)變增量為零,因此可得

由于塑性因子λ 在加載階段為大于零的數(shù)值因此顯然可知,平面應(yīng)變條件的等價方程為

由鏈式法則,可依次求取得到

其余變量可參見式(85)～式(88).

8 準則及變換應(yīng)力法的試驗驗證

為了便于對所提各向異性t 準則及其變換應(yīng)力公式進行驗證,分別采用如下的巖土材料,在真三軸以及平面應(yīng)變加載條件下的破壞以及應(yīng)力應(yīng)變關(guān)系試驗結(jié)果進行對比分析.表1 為所有巖土材料的參數(shù).

8.1 真三軸強度準則預(yù)測

圖15 中圓圈表示Lade 等[39]關(guān)于Grundite 黏土的真三軸加載測試結(jié)果,對該試驗結(jié)果在圖15 所示的偏平面進行預(yù)測對比.黏土除了具備應(yīng)力誘導(dǎo)各向異性,也就是隨著應(yīng)力羅德角的不同而表現(xiàn)出不同的偏應(yīng)力強度以外,還具有一定程度的原生各向異性,也就是由于沉積過程的特點所形成的近似為橫觀各向同性的材料.采用所提t(yī) 準則可以較好地描述該黏土在真三軸下的強度特性.

表1 巖土材料參數(shù)Table 1 Geomaterial parameters

圖15 偏平面Grundite 黏土的真三軸測試結(jié)果與預(yù)測對比Fig.15 Comparison results with test data and prediction for Grundite clay in deviatoric plane

圖16 中測試點為Lam 等[40]關(guān)于Toyoura 砂土的真三軸測試結(jié)果.由圖可見,測試點在偏平面上形成了較為尖銳的曲邊三角形形態(tài),且三軸伸長較三軸壓縮強度值具有顯著的減小特點.采用所提準則可以較好地模擬砂土的這種特性.

圖17 中離散點為Mogi[41]關(guān)于粗面巖巖石的真三軸測試結(jié)果.巖石具有更為顯著的原生各向異性性質(zhì),且應(yīng)力誘導(dǎo)各向異性對于強度值的影響也很明顯.圖中采用所提準則可以很好地預(yù)測巖石的真三軸強度特性.

圖16 偏平面Toyoura 砂土的真三軸測試結(jié)果與預(yù)測對比Fig.16 Comparison results with test data and prediction for Toyoura sand in deviatoric plane

圖17 偏平面粗面巖石的真三軸測試結(jié)果與預(yù)測對比Fig.17 Comparison results with test data and prediction for Trachyte rock in deviatoric plane

8.2 平面應(yīng)變條件強度準則預(yù)測

圖18 中離散點為Park 等[42]關(guān)于3 種砂土在平面應(yīng)變條件下的內(nèi)摩擦角特性得到的試驗研究結(jié)果.其中角度δ 表示大主應(yīng)力與沉積面法線的夾角,而φ 表示內(nèi)摩擦角.由試驗曲線可見,隨著大主應(yīng)力與沉積面法線夾角δ 的單調(diào)增大,從0°到90°區(qū)間,內(nèi)摩擦角呈現(xiàn)先減小后增大的規(guī)律特點,而對于Monetary 砂土,則呈現(xiàn)了單調(diào)減小的規(guī)律,與Hostun 砂土和Silica 砂土的規(guī)律不一致,表明Monetary 砂土在δ趨近90°時的測試值有待商榷.采用所提的t 準則在平面應(yīng)變條件下,3 條曲線也呈現(xiàn)了先減小后增大的規(guī)律特點,且δ=0°對應(yīng)的值要大于δ=90°對應(yīng)的值.

圖18 平面應(yīng)變下3 種砂土內(nèi)摩擦角測試結(jié)果與預(yù)測對比Fig.18 Comparison results with test data and prediction for three kinds of sand in plane strain condition

圖19 平面應(yīng)變下兩種砂土內(nèi)摩擦角測試結(jié)果與預(yù)測對比Fig.19 Comparison results with test data and prediction for two kinds of sand in plane strain condition

圖19 中測試點為Oda 等[43]關(guān)于兩種Toyoura砂土在平面應(yīng)變條件下的內(nèi)摩擦角值.δ 角仍與前述相同,由兩種測試結(jié)果可見,內(nèi)摩擦角曲線隨δ 角的曲線形態(tài)表現(xiàn)為明顯的“湯匙”形狀.一方面仍然是δ=0°對應(yīng)的值要大于δ=90°對應(yīng)的內(nèi)摩擦角值,另一方面,最小的內(nèi)摩擦角所對應(yīng)的δ 角向右側(cè)偏移,集中在60°～80°之間.說明當(dāng)δ 角在最低曲線時刻,此時破壞面與沉積面相互平行,此時達到最小的內(nèi)摩擦角.采用所提準則可以較好地描述上述特性.

8.3 變換應(yīng)力法應(yīng)用預(yù)測

圖20 中測試點為Ward 等[44]和Indraratna 等[45]關(guān)于倫敦裂隙黏土所做的常規(guī)不排水三軸壓縮測試結(jié)果.顯然裂隙可以視為某種沉積面的弱化物態(tài)狀態(tài),則圖中廣義偏應(yīng)力強度分為大、中、小3 條曲線,分別對應(yīng)著大主應(yīng)力方向與裂隙面法線夾角為0°以及某種銳角,和90°三種情況下的測試結(jié)果.采用基于t 準則的變換應(yīng)力法修正UH 模型[46],用于反映超固結(jié)度以及裂隙等原生各向異性對應(yīng)力應(yīng)變關(guān)系的影響.由圖可見,采用基于t 準則的變換應(yīng)力法能夠有效地反映原生各向異性所帶來的對應(yīng)力應(yīng)變關(guān)系的影響.

圖20 倫敦裂隙黏土應(yīng)力應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.20 Comparison results with test data and prediction of stress-strain relationship for fissured London clay

圖21～圖28 中離散點為Yang 等[47]關(guān)于Leighton Buzzard 砂土在原生各向異性以及應(yīng)力羅德角影響下的應(yīng)力應(yīng)變關(guān)系測試結(jié)果.圖中連續(xù)曲線為采用所提的變換應(yīng)力法修正的UH 模型做出的預(yù)測結(jié)果.其中,對于原生各向異性的影響,利用大主應(yīng)力與沉積面法線夾角δ 呈現(xiàn)為0°,30°,60°,90°四種情況作為初始加載條件,在上述4 種工況下再分別考察b=0,0.2,0.5,1 下的應(yīng)力應(yīng)變關(guān)系曲線.

圖21 為大主應(yīng)力垂直于沉積面加載時的測試與預(yù)測對比結(jié)果,由圖21 可見,應(yīng)力比與偏應(yīng)變的關(guān)系,采用修正后模型可以較好地反映b值增大,而應(yīng)力比強度減小的規(guī)律,以及應(yīng)變軟化現(xiàn)象.圖22 為所對應(yīng)的體應(yīng)變隨偏應(yīng)變的關(guān)系對比結(jié)果,隨著b值的增大剪脹體應(yīng)變逐漸減小,所用模型可較好地反映上述特點.

圖23 與圖24 為對應(yīng)δ=30°下的應(yīng)力應(yīng)變關(guān)系對比結(jié)果.除了反映隨b值增大應(yīng)力比峰值減小的規(guī)律以外,由于δ 的增大,對應(yīng)的各條曲線的峰值應(yīng)力比較δ=0°工況下有所減小,采用所提變換應(yīng)力法修正后模型可較好地反映上述規(guī)律.圖24 中體變特性也能用模型較好的描述.

圖21 δ=0° 條件下Leighton Buzzard 砂土應(yīng)力比與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.21 Comparison results with test data and prediction of stress ratio and deviatoric strain relationship for Leighton Buzzard sand under δ=0°

圖22 δ=0° 條件下Leighton Buzzard 砂土體應(yīng)變與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.22 Comparison results with test data and prediction of volume strain and deviatoric strain relationship for Leighton Buzzard sand under δ=0°

圖25 與圖26 為對應(yīng)δ=60°下的應(yīng)力應(yīng)變關(guān)系對比結(jié)果.由圖可見,隨著δ 的增大,各條曲線的應(yīng)力比強度值出現(xiàn)小幅度的減小,用修正后模型能夠反映上述規(guī)律.圖26 中的體應(yīng)變的剪脹量較δ=30°工況有所減小,用修正模型可以反映上述特點.

圖27 與圖28 為對應(yīng)δ=90°下的應(yīng)力應(yīng)變關(guān)系對比結(jié)果.由圖可見,測試結(jié)果表明,各條曲線峰值應(yīng)力比仍出現(xiàn)了小幅度的減小,但所用修正模型過高地預(yù)估了應(yīng)力比強度值.圖28 中的預(yù)測的剪脹體變值也稍大于實測值.

圖23 δ=30° 條件下Leighton Buzzard 砂土應(yīng)力比與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.23 Comparison results with test data and prediction of stress ratio and deviatoric strain relationship for Leighton Buzzard sand under δ=30°

圖24 δ=30° 條件下Leighton Buzzard 砂土體應(yīng)變與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.24 Comparison results with test data and prediction of volume strain and deviatoric strain relationship for Leighton Buzzard sand under δ=30°

9 結(jié)論

基于已提出的t 強度準則,考察橫觀各向同性性質(zhì),利用由主應(yīng)力構(gòu)成的空間滑移面與材料宏觀上的物理沉積面之間的夾角作為度量各向異性程度的狀態(tài)量.同理,將其推廣于三維正交各向異性的描述,提出了用于描述三維各向異性的狀態(tài)量.并在此基礎(chǔ)上,給出了基于各向異性t 準則的變換應(yīng)力方程.主要完成了如下幾點工作:

圖25 δ=60° 條件下Leighton Buzzard 砂土應(yīng)力比與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.25 Comparison results with test data and prediction of stress ratio and deviatoric strain relationship for Leighton Buzzard sand under δ=60°

圖26 δ=60° 條件下Leighton Buzzard 砂土體應(yīng)變與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.26 Comparison results with test data and prediction of volume strain and deviatoric strain relationship for Leighton Buzzard sand under δ=60°

(1)基于已有的t 準則公式,通過弓入有效滑移面與沉積面之間的夾角作為各向異性程度狀態(tài)量,建立了三維各向異性t 準則.

(2) 基于上述所提出的各向異性t 準則,得到了變換應(yīng)力方程,可方便的將既有的以p-q為變量的二維模型推廣為可考慮各向異性影響的三維應(yīng)力應(yīng)變關(guān)系模型.

(3)通過真三軸、平面應(yīng)變條件的強度測試對比,以及真三軸下的應(yīng)力應(yīng)變關(guān)系預(yù)測對比,進一步驗證了所提準則及變換應(yīng)力法的合理性與適用性.

圖27 δ=90° 條件下Leighton Buzzard 砂土應(yīng)力比與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.27 Comparison results with test data and prediction of stress ratio and deviatoric strain relationship for Leighton Buzzard sand under δ=90°

圖28 δ=90° 條件下Leighton Buzzard 砂土體應(yīng)變與偏應(yīng)變關(guān)系測試結(jié)果與預(yù)測對比Fig.28 Comparison results with test data and prediction of volume strain and deviatoric strain relationship for Leighton Buzzard sand under δ=90°