應(yīng)力旋轉(zhuǎn)條件下軟黏土非共軸角計(jì)算方法

嚴(yán)佳佳，林清輝，傅了一

?

應(yīng)力旋轉(zhuǎn)條件下軟黏土非共軸角計(jì)算方法

嚴(yán)佳佳1，林清輝2，傅了一1

(1. 浙江大學(xué)濱海和城市巖土工程研究中心，浙江杭州，310058；2. 臺州職業(yè)技術(shù)學(xué)院建筑工程學(xué)院，浙江臺州，318000)

基于空心圓柱扭剪儀對軟黏土開展主應(yīng)力純旋轉(zhuǎn)、往復(fù)旋轉(zhuǎn)及旋轉(zhuǎn)同時(shí)增加剪應(yīng)力這3類典型應(yīng)力旋轉(zhuǎn)路徑下的不排水試驗(yàn)，重點(diǎn)分析應(yīng)力旋轉(zhuǎn)引起的應(yīng)變增量方向的變化規(guī)律及其影響因素。研究結(jié)果表明：3類應(yīng)力旋轉(zhuǎn)路徑引起的應(yīng)變增量方向與主應(yīng)力方向之間存在顯著的非共軸現(xiàn)象；隨著主應(yīng)力軸的連續(xù)旋轉(zhuǎn)，非共軸角基本以90°為周期波動(dòng)變化，當(dāng)旋轉(zhuǎn)方向反向時(shí)，應(yīng)變增量方向也跟隨突變，應(yīng)力旋轉(zhuǎn)同時(shí)增加剪應(yīng)力，非共軸角有減小的趨勢；剪應(yīng)力水平、應(yīng)力旋轉(zhuǎn)方向及應(yīng)力路徑對非共軸角的變化規(guī)律影響較顯著，中主應(yīng)力系數(shù)對非共軸角的影響則很小。基于Gutierrze提出的應(yīng)力映射方法建立的非共軸角計(jì)算方法能較好地模擬應(yīng)力旋轉(zhuǎn)路徑下軟黏土非共軸角的變化規(guī)律，并能反映剪應(yīng)力水平、應(yīng)力旋轉(zhuǎn)方向及應(yīng)力路徑的影響。

應(yīng)力軸旋轉(zhuǎn)；軟黏土；應(yīng)變；非共軸

應(yīng)力方向旋轉(zhuǎn)是巖土工程中普遍存在的一種典型應(yīng)力路徑。在經(jīng)典塑性理論中基于金屬材料的塑性特性，圣維南提出了材料塑性應(yīng)變增量方向與當(dāng)前應(yīng)力方向一致的假設(shè)，即圣維南假設(shè)。但在很多條件下該假設(shè)對巖土材料并不成立。已有大量的試驗(yàn)結(jié)果都表明土體塑性應(yīng)變增量方向與應(yīng)力方向存在不一致的現(xiàn)象，尤其在主應(yīng)力軸旋轉(zhuǎn)等復(fù)雜非比例加載應(yīng)力路徑條件下主應(yīng)力和主應(yīng)變增量方向存在顯著的偏 差[1?11]，土體塑性應(yīng)變增量方向一般介于應(yīng)力方向和應(yīng)力增量方向之間。土體材料的共軸性假設(shè)一般僅在比例加載條件下近似成立。相對于傳統(tǒng)的共軸性假設(shè)，將土體的這種塑性應(yīng)變增量方向與應(yīng)力方向的不一致性定義為土體的非共軸特性。在應(yīng)力方向旋轉(zhuǎn)路徑下，塑性應(yīng)變增量方向一般介于應(yīng)力和應(yīng)力增量方向之間。合理描述應(yīng)力旋轉(zhuǎn)條件下土體的變形規(guī)律，其關(guān)鍵在于確定應(yīng)力旋轉(zhuǎn)過程中塑性流動(dòng)方向，即土體的非共軸塑性流動(dòng)規(guī)律。TSUTSUMI等[12]的分析表明主應(yīng)力在旋轉(zhuǎn)應(yīng)力路徑下，只有同時(shí)考慮非共軸和各向異性特性的本構(gòu)模型才能較好地模擬試驗(yàn)結(jié)果。此外，非共軸對土體變形分叉點(diǎn)的預(yù)測也有重要影響[13]?；谠囼?yàn)研究結(jié)果很多研究者從塑性流動(dòng)理論出發(fā)建立了相應(yīng)的土體非共軸塑性流動(dòng)模型[13?18]。RUDNICKI等[14]最先在Drucker?Prager模型上附加了非共軸塑性應(yīng)變分量，得到角點(diǎn)屈服非共軸模型，計(jì)算結(jié)果表明考慮非共軸修正后可以大大促進(jìn)土體局部化變形失穩(wěn)。GUTIERREZ等[2]結(jié)合邊界面理論的思想給出了一個(gè)確定砂土塑性應(yīng)變增量方向的方法，該方法通過建立塑性應(yīng)變增量方向與應(yīng)力和應(yīng)力增量方向的聯(lián)系，可以反映土體的非共軸特性，但不能很好地反映非共軸角波動(dòng)變化的規(guī)律。LASHKARI等[15]也是基于邊界面理論，同時(shí)考慮各向異性影響提出了砂土的非共軸模型，將塑性應(yīng)變增量分解為2個(gè)正交的分量并考慮兩者比值隨主應(yīng)力軸旋轉(zhuǎn)的變化，得到的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合得更好，但該方法的模型參數(shù)較多且不易確定。錢建固等[16]從理論上證明了塑性應(yīng)變率的流動(dòng)法則只是Drucker穩(wěn)定性公設(shè)的充分條件，而不是必要條件，即在非共軸條件下也能滿足Drucker公設(shè)的條件，認(rèn)為非共軸特性并不是土體材料所特有的性質(zhì)，在金屬材料中同樣也存在。扈萍等[17]在已有的本構(gòu)模型中引入非共軸塑性流動(dòng)理論建立了砂土非共軸模型，將應(yīng)力增量分解為共軸和非共軸分量，相應(yīng)的塑性應(yīng)變增量則與應(yīng)力增量分量是共軸的。但對于應(yīng)力純旋轉(zhuǎn)路徑，非共軸角始終為45°。董建國等[18]基于三維模型的塑性流動(dòng)特點(diǎn)，提出三維彈塑性模型的非共軸修正項(xiàng)，該模型主要考慮土體的各向異性及非共軸對剪切帶形成的影響。在實(shí)際工程設(shè)計(jì)方面，YANG等[1]的研究結(jié)果表明不考慮土體非共軸特性可能會使設(shè)計(jì)結(jié)果偏于不安全。隨著考慮主應(yīng)力軸旋轉(zhuǎn)影響軟黏土變形特性試驗(yàn)研究的不斷深入，軟黏土非共軸變形特性也引起了廣大研究者的關(guān)注。但從以上的研究結(jié)果可以發(fā)現(xiàn)：一方面合理確定非共軸角對準(zhǔn)確預(yù)測應(yīng)力旋轉(zhuǎn)引起的變形有重要作用，另一方面目前對簡單有效的軟黏土非共軸角的確定方法的研究還較少。因此，本文作者首先總結(jié)3類典型應(yīng)力旋轉(zhuǎn)路徑下軟黏土非共軸特性的試驗(yàn)研究結(jié)果，再基于邊界面理論，根據(jù)應(yīng)力映射法則，對應(yīng)力旋轉(zhuǎn)路徑下軟黏土非共軸角確定方法進(jìn)行研究，以便為考慮主應(yīng)力軸旋轉(zhuǎn)影響土體本構(gòu)建立提供研究基礎(chǔ)。

1 試驗(yàn)研究

1.1 試驗(yàn)儀器和用土

采用浙江大學(xué)ZJU?HCA空心圓柱扭剪儀來實(shí)現(xiàn)相應(yīng)應(yīng)力路徑下的試驗(yàn)。該儀器可以對空心圓柱試樣施加獨(dú)立控制的軸力、扭矩和內(nèi)、外壓力，通過軸力和扭矩的耦合以及內(nèi)、外壓力的變化，實(shí)現(xiàn)主應(yīng)力軸連續(xù)旋轉(zhuǎn)等復(fù)雜試驗(yàn)應(yīng)力路徑[19]。其中軸向應(yīng)力z和扭剪應(yīng)力zθ分別由軸力和扭矩獨(dú)立控制，徑向應(yīng)力r和切向應(yīng)力θ的控制則由內(nèi)、外壓力的變化來實(shí)現(xiàn)。關(guān)于儀器的具體描述及應(yīng)力應(yīng)變分量的計(jì)算以及試驗(yàn)數(shù)據(jù)的可靠性和可重復(fù)性等參見文獻(xiàn)[19?20]。

試驗(yàn)所用土樣為取自杭州某基坑開挖工地的原狀黏土，取土深度為地表以下約4 m處。土樣天然含水率=45.0%，相對密度s=2.70，密度=1.75 g/cm3，塑限L=22.6%，液限P=46.0%，塑性指數(shù)P=23.4。通過專用制樣設(shè)備將原狀黏土試塊制備成高度×外徑×內(nèi)徑為200 mm×100 mm×60 mm的空心試樣，具體取樣及制樣方法參見文獻(xiàn)[19]?？招脑嚇咏?jīng)反壓飽(100 kPa)和后孔壓系數(shù)都能達(dá)到0.97以上，然后對試樣進(jìn)行固結(jié)。

1.2 試驗(yàn)方案和應(yīng)力路徑

對原狀黏土開展3類典型應(yīng)力旋轉(zhuǎn)路徑，分別為主應(yīng)力純旋轉(zhuǎn)、主應(yīng)力旋轉(zhuǎn)的同時(shí)增加剪應(yīng)力以及主應(yīng)力旋轉(zhuǎn)后逆向旋轉(zhuǎn)。試樣飽和完成后先在150 kPa有效固結(jié)壓力下進(jìn)行等向固結(jié)，固結(jié)完成后在排水條件下調(diào)節(jié)各參數(shù)(剪應(yīng)力為，大主應(yīng)力方向角為)到達(dá)目標(biāo)值，并始終保持總應(yīng)力150 kPa不變，然后在保持=150 kPa，=0(起始剪應(yīng)力)條件下調(diào)節(jié)中主應(yīng)力系數(shù)到目標(biāo)值(表1)，最后開展相應(yīng)應(yīng)力路徑下的應(yīng)力旋轉(zhuǎn)試驗(yàn)。試驗(yàn)方案見表1，通過該試驗(yàn)方案分析，剪應(yīng)力水平、中主應(yīng)力系數(shù)及應(yīng)力旋轉(zhuǎn)方向?qū)?yīng)力旋轉(zhuǎn)條件下軟黏土非共軸角變化規(guī)律的影響。平均主應(yīng)力；剪應(yīng)力；中主應(yīng)力系數(shù)；大主應(yīng)力方向角，其中，，和分別為軸向、徑向、環(huán)向和扭剪應(yīng)力。

表1 試驗(yàn)方案設(shè)計(jì)

在應(yīng)力旋轉(zhuǎn)坐標(biāo)平面內(nèi)試樣應(yīng)力路徑如圖1所示，等向固結(jié)完成時(shí)，試樣應(yīng)力狀態(tài)點(diǎn)位于坐標(biāo)原點(diǎn)，沿方向增加剪應(yīng)力至旋轉(zhuǎn)起始點(diǎn)時(shí)的應(yīng)力狀態(tài)點(diǎn)為橫坐標(biāo)上的點(diǎn)如點(diǎn)。應(yīng)力點(diǎn)至原點(diǎn)的距離即表示剪應(yīng)力的大小，應(yīng)力點(diǎn)與橫坐標(biāo)正向的夾角為2倍大主應(yīng)力方向角。純旋轉(zhuǎn)試驗(yàn)PR1~PR5的應(yīng)力路徑為→→→→(=37.5 kPa)，不同的試驗(yàn)在該平面內(nèi)的旋轉(zhuǎn)路徑是相同的。應(yīng)力正逆向旋轉(zhuǎn)試驗(yàn)(CR)的應(yīng)力路徑為→→。主應(yīng)力旋轉(zhuǎn)的同時(shí)增加剪應(yīng)力試驗(yàn)(TR)的路徑為→。

圖1 應(yīng)力路徑示意圖

2 試驗(yàn)研究結(jié)果分析

2.1 塑性的應(yīng)變增量方向

土體的非共軸特性指的是塑性應(yīng)變增量方向與當(dāng)前應(yīng)變增量方向的不一致性，一般采用非共軸角來描述土體變形過程中的非共軸程度。因此，首先需要確定應(yīng)力旋轉(zhuǎn)引起土體塑性應(yīng)變增量的方向。應(yīng)力旋轉(zhuǎn)路徑條件下土體應(yīng)變增量可以表示為

而總應(yīng)變增量又包含了彈性應(yīng)變和塑性應(yīng)變，對于平面主應(yīng)力旋轉(zhuǎn)路徑(2個(gè)主應(yīng)力在垂直第3個(gè)主應(yīng)力的平面內(nèi)旋轉(zhuǎn))，參考ZDRAVKOVIC等[21?22]的研究結(jié)果，主應(yīng)力僅繞某一應(yīng)力主軸旋轉(zhuǎn)時(shí)引起的彈性扭剪應(yīng)變?yōu)?/p>

式中：為彈性剪切模量，而其他各應(yīng)變分量為0，即

因此，彈性應(yīng)變張量為

與砂土相比，軟黏土的彈性區(qū)域很小，從工程實(shí)際角度可以忽略彈性應(yīng)變分量的影響。ZHOU等[9]采用循環(huán)加載的方式確定軟土彈性模量的方法分析了彈性應(yīng)變分量對軟土非共軸角的影響，發(fā)現(xiàn)其影響主要表現(xiàn)在旋轉(zhuǎn)初期應(yīng)變較小的階段，隨著應(yīng)變增加，彈性應(yīng)變分量的影響基本可以忽略。因此，本文參考CAI等提出的[3, 10]的方法，采用總應(yīng)變增量代替塑性應(yīng)變增量來研究軟黏土的非共軸特性。主應(yīng)力旋轉(zhuǎn)引起的主應(yīng)變增量方向角為

圖2 非共軸示意圖

2.2 應(yīng)力旋轉(zhuǎn)路徑下軟黏土非共軸角變化規(guī)律

在3種典型非比例加載應(yīng)力路徑條件下，隨主應(yīng)力旋轉(zhuǎn)原狀軟黏土應(yīng)變增量方向角的變化規(guī)律如圖3所示。

(a) 純旋轉(zhuǎn)應(yīng)力路徑；(b) 應(yīng)力反向旋轉(zhuǎn)路徑；(c) 剪應(yīng)力同步增加應(yīng)力路徑

其中圖3(a)所示為應(yīng)力純旋轉(zhuǎn)路徑下，不同中主應(yīng)力系數(shù)和不同剪應(yīng)力條件下原狀軟黏土非共軸角隨大主應(yīng)力方向旋轉(zhuǎn)的變化規(guī)律。從圖3(a)可以看到：主應(yīng)力旋轉(zhuǎn)1周(180°)的過程中，軟黏土非共軸角基本以90°的周期變化，1個(gè)周期內(nèi)非共軸角表現(xiàn)出先減小后增大的變化規(guī)律。圖中試樣PR1~PR5的非共軸角的變化規(guī)律基本相同，因此，中主應(yīng)力系數(shù)對非共軸角的影響不顯著。而對比試樣PR5和PR6的試驗(yàn)結(jié)果，可以發(fā)現(xiàn)剪應(yīng)力水平對非共軸角的影響則相對較大，=25 kPa(PR6)時(shí)非共軸角比=37.5 kPa(PR5)時(shí)的角大。這些試驗(yàn)規(guī)律與相同應(yīng)力路徑條件下砂土的試驗(yàn)研究結(jié)果[2?4]基本相同，因此，非共軸角變化規(guī)律主要由應(yīng)力路徑?jīng)Q定。

圖3(b)所示為應(yīng)力往復(fù)純旋轉(zhuǎn)路徑條件下軟黏土非共軸角的變化規(guī)律，圖中也給出了純旋轉(zhuǎn)試驗(yàn)PR6 在0~45°階段非共軸角的變化規(guī)律。從圖3(b)可以看到：不同試樣(PR6和CR)在相同應(yīng)力路徑條件下，非共軸角的變化規(guī)律基本一致，這說明試驗(yàn)結(jié)果的可重復(fù)性。試驗(yàn)結(jié)果表明：當(dāng)應(yīng)力旋轉(zhuǎn)方向發(fā)生逆轉(zhuǎn)時(shí)(此時(shí)應(yīng)力增量方向反向(圖2))，應(yīng)變增量的方向也同時(shí)發(fā)生逆轉(zhuǎn)，非共軸角發(fā)生突變。這表明應(yīng)變增量方向受應(yīng)力增量方向的影響顯著，PRADEL等[11]對砂土開展的應(yīng)力探測試驗(yàn)(應(yīng)力沿不同方向施加1個(gè)微小的增量)也發(fā)現(xiàn)了相同的規(guī)律，即應(yīng)變增量的方向隨應(yīng)力增量方向的變化而變化。

圖3(c)所示為應(yīng)力旋轉(zhuǎn)同時(shí)增加剪應(yīng)力路徑條件下原狀軟黏土非共軸角的變化規(guī)律。由圖3(a)可見：應(yīng)力純旋轉(zhuǎn)時(shí)在50°~90°旋轉(zhuǎn)區(qū)間內(nèi)非共軸角隨著應(yīng)力旋轉(zhuǎn)先增加后基本保持不變，而當(dāng)旋轉(zhuǎn)的同時(shí)增加剪應(yīng)力時(shí)，非共軸角則隨著應(yīng)力旋轉(zhuǎn)而逐漸減小。這說明在旋轉(zhuǎn)過程中剪應(yīng)力的同步增加也會使非共軸角趨于減小，這與CAI等[3]對砂土進(jìn)行的試驗(yàn)得到的規(guī)律是一致的。

2.3 試驗(yàn)結(jié)果總結(jié)

由以上不同應(yīng)力旋轉(zhuǎn)路徑下的軟黏土非共軸特性試驗(yàn)研究結(jié)果可知，在應(yīng)力方向旋轉(zhuǎn)過程中，軟黏土?xí)a(chǎn)生相應(yīng)的應(yīng)變，并且應(yīng)變增量方向介于當(dāng)前應(yīng)力和應(yīng)力增量方向之間，即應(yīng)變增量方向與應(yīng)力方向存在顯著的非共軸現(xiàn)象。應(yīng)力旋轉(zhuǎn)條件下軟黏土應(yīng)變非共軸角的變化規(guī)律及影響因素可總結(jié)為：

1) 隨著應(yīng)力方向的旋轉(zhuǎn)非共軸角呈波動(dòng)變化的規(guī)律，且變化周期為90°左右；

2) 隨著剪應(yīng)力水平的增大，非共軸角有減小的趨勢；

3) 中主應(yīng)力系數(shù)對軟黏土非共軸角的影響不 顯著；

4) 應(yīng)變增量方向隨應(yīng)力增量方向的突變而突變，非共軸角變化規(guī)律主要由應(yīng)力路徑的決定；

5) 應(yīng)力方向旋轉(zhuǎn)的同時(shí)增加剪應(yīng)力使非共軸角趨于減小。

3 應(yīng)力旋轉(zhuǎn)下軟黏土非共軸角計(jì)算方法

基于以上試驗(yàn)研究結(jié)果的總結(jié)和分析，參考GUTIERREZ等[2]提出的砂土非共軸角計(jì)算方法，在邊界面理論框架下，對應(yīng)力旋轉(zhuǎn)條件下原狀軟黏土非共軸角計(jì)算方法進(jìn)行研究。

3.1 主應(yīng)力旋轉(zhuǎn)條件下應(yīng)變增量方向角計(jì)算

已有的砂土及本文黏土試驗(yàn)研究的結(jié)果都表明主應(yīng)力旋轉(zhuǎn)引起的土體應(yīng)變增量方向一般介于應(yīng)力方向和應(yīng)力增量方向之間，并同時(shí)由當(dāng)前應(yīng)力方向和應(yīng)力增量方向共同決定。BLANC等[23]的研究結(jié)果也表明主應(yīng)力軸純旋轉(zhuǎn)條件下應(yīng)變增量由共軸和非共軸應(yīng)變機(jī)理共同決定。因此，計(jì)算非共軸角時(shí)需要綜合考慮主應(yīng)力及其增量方向的影響。本節(jié)在應(yīng)力旋轉(zhuǎn)平面內(nèi)研究應(yīng)力旋轉(zhuǎn)引起軟黏土應(yīng)變增量方向角的計(jì)算方法。

也有人說，自焚是有道理的，李駟峋是個(gè)饕餮之徒，大魚大肉吃多了，腸肥腦滿，肚子里板油太厚，整個(gè)人就像一根蠟燭，沾火就著。

具體確定方法如圖4所示，圖中的圓和橢圓曲線分別表示偏應(yīng)力平面內(nèi)純旋轉(zhuǎn)應(yīng)力路徑和破壞邊界面。對于應(yīng)力旋轉(zhuǎn)應(yīng)力路徑上旋轉(zhuǎn)至任意點(diǎn)()處引起的應(yīng)變增量方向，可以由點(diǎn)處應(yīng)力增量方向與破壞邊界面交點(diǎn)到原點(diǎn)的方向來確定。由破壞邊界面上應(yīng)力映射點(diǎn)(，)與坐標(biāo)原點(diǎn)連線的方向角2dε，即可確定該點(diǎn)處應(yīng)變增量的方向角dε。

圖4 應(yīng)變增量方向確定方法

首先需要確定破壞邊界面上應(yīng)力映射點(diǎn)()的坐標(biāo)。對于應(yīng)力旋轉(zhuǎn)路徑上任意點(diǎn)()，沿著應(yīng)力增量方向可以將其映射到破壞邊界面上()。應(yīng)力旋轉(zhuǎn)平面內(nèi)主應(yīng)力軸純旋轉(zhuǎn)的應(yīng)力路徑方程為(任意旋轉(zhuǎn)應(yīng)力路徑都可以表示為和的函數(shù))：

式中：d為應(yīng)力點(diǎn)到原點(diǎn)的距離；為大主應(yīng)力方 向角。

該應(yīng)力路徑上任意點(diǎn)處的應(yīng)力增量方向的斜率為

應(yīng)力增量方向所在直線的方程可以表示為

應(yīng)力映射點(diǎn)(，)同時(shí)也是破壞邊界面上的應(yīng)力點(diǎn)，其坐標(biāo)可以通過聯(lián)立式(7)~(9)，得

由以上的分析可知，對于某一應(yīng)力狀態(tài)下的軟黏土，其破壞邊界面是確定的。旋轉(zhuǎn)過程中剪應(yīng)力的大小決定了應(yīng)力點(diǎn)與破壞邊界面之間的距離。大主應(yīng)力及其增量的方向決定了應(yīng)力映射點(diǎn)的位置。因此，本文試驗(yàn)研究得到的軟黏土非共軸特性的主要影響因素(剪應(yīng)力、應(yīng)力路徑、應(yīng)力旋轉(zhuǎn)方向)都能在該方法中得到體現(xiàn)。

中主應(yīng)力系數(shù)對土體的強(qiáng)度有一定的影響，故中主應(yīng)力系數(shù)的影響體現(xiàn)在破壞邊界面的位置不同。由于目前尚缺少這方面的試驗(yàn)數(shù)據(jù)，還不能在破壞邊界面方程中體現(xiàn)中主應(yīng)力的影響，因此，該方法目前還不能考慮中主應(yīng)力系數(shù)對非共軸角的影響。

3.2 應(yīng)力旋轉(zhuǎn)平面內(nèi)軟黏土破壞邊界面

應(yīng)力旋轉(zhuǎn)路徑一般通過對空心試樣施加耦合的軸力和扭矩實(shí)現(xiàn)，其受力狀態(tài)與三軸試驗(yàn)不同，應(yīng)力在與徑向應(yīng)力垂直的平面內(nèi)旋轉(zhuǎn)。沈揚(yáng)等[19?20]采用空心圓柱扭剪儀對杭州原狀黏土開展了=0.50時(shí)沿不同大主應(yīng)力方向的不排水定向剪切試驗(yàn)，對杭州軟黏土的破壞規(guī)律進(jìn)行研究。本文根據(jù)他們的試驗(yàn)結(jié)果，將廣義剪應(yīng)變q=8%或峰值剪應(yīng)力處的剪應(yīng)力定義為破壞應(yīng)力，將破壞點(diǎn)處的應(yīng)力值繪于偏應(yīng)力坐標(biāo)平面內(nèi)，如圖5所示。對試驗(yàn)點(diǎn)的擬合可以得到應(yīng)力旋轉(zhuǎn)平面內(nèi)破壞邊界面方程為

圖5 原狀黏土各向異性強(qiáng)度破壞邊界面

式(11)表示為1條橢圓曲線。式中：q0為沿沉積方向剪切時(shí)土體的破壞強(qiáng)度；為沿垂直沉積方向剪切時(shí)土體的破壞強(qiáng)度。由于＞，因此，橢圓的圓心偏向原點(diǎn)右側(cè)。雖然試驗(yàn)點(diǎn)較少，但式(11)還是能較好地體現(xiàn)原狀黏土強(qiáng)度各向異性的特性。考慮軟黏土為橫觀各向同性材料，則0°~90°和0°~?90°的破壞邊界面是對稱的，即在應(yīng)力旋轉(zhuǎn)平面內(nèi)原狀軟黏土的破壞邊界面可以表示為

4 計(jì)算結(jié)果與試驗(yàn)結(jié)果的對比

按照本文試驗(yàn)應(yīng)力路徑，計(jì)算得到的不同應(yīng)力旋轉(zhuǎn)路徑條件下軟黏土非共軸角與試驗(yàn)值的對比結(jié)果如圖6所示。從圖6可以看到：采用本文提出的原狀軟黏土各向異性破壞邊界面及考慮應(yīng)力旋轉(zhuǎn)影響應(yīng)變增量方向確定方法計(jì)算得到不同應(yīng)力旋轉(zhuǎn)路徑下非共軸角的變化規(guī)律，基本能較好地模擬非共軸角的變化規(guī)律及趨勢。

(a) 純旋轉(zhuǎn)應(yīng)力路徑；(b) 往復(fù)旋轉(zhuǎn)應(yīng)力路徑；(c) 應(yīng)力旋轉(zhuǎn)同時(shí)增加剪應(yīng)力

由于本文提出的非共軸角確定方法主要為便于實(shí)際應(yīng)用的簡化方法，因此，僅考慮基于試驗(yàn)研究確定的主要影響因素，從計(jì)算結(jié)果也能看到該方法能較好地反映剪應(yīng)力水平、應(yīng)力旋轉(zhuǎn)方向，應(yīng)力旋轉(zhuǎn)與剪切耦合作用等主要因素對非共軸角的影響。

由于破壞邊界面尚未考慮中主應(yīng)力系數(shù)的影響，因此，該方法還不能考慮中主應(yīng)力系數(shù)的影響。但從試驗(yàn)規(guī)律可知，中主應(yīng)力系數(shù)的影響相對不顯著。此外，該方法也未考慮應(yīng)力旋轉(zhuǎn)引起的共軸和非共軸變形的耦合作用。因此，計(jì)算結(jié)果在數(shù)值上與試驗(yàn)結(jié)果有一定的誤差。采用該方法沒有額外的土體參數(shù)，破壞邊界面則可以采用常規(guī)三軸試驗(yàn)(沿不同方向切取試樣)確定。

5 結(jié)論

1) 應(yīng)力旋轉(zhuǎn)引起的應(yīng)變增量方向介于應(yīng)力和應(yīng)力增量方向之間，隨著應(yīng)力方向的純旋轉(zhuǎn)，非共軸角呈波動(dòng)變化的規(guī)律；隨著剪應(yīng)力水平的增大，非共軸角有減小趨勢；中主應(yīng)力系數(shù)對軟黏土非共軸角的影響不顯著，但剪應(yīng)力水平對非共軸角的影響較顯著；應(yīng)變增量方向隨應(yīng)力增量方向的突變而突變。

2) 根據(jù)定向剪切試驗(yàn)結(jié)果，原狀黏土在應(yīng)力旋轉(zhuǎn)平面內(nèi)的破壞邊界面可以表示為1條橢圓曲線。結(jié)合Gutierrez提出的應(yīng)力映射方法，建立的軟黏土非共軸角計(jì)算方法能較好地描述3類典型應(yīng)力旋轉(zhuǎn)路徑下軟黏土的非共軸角變化規(guī)律，并能反映剪應(yīng)力水平和應(yīng)力旋轉(zhuǎn)方向的影響，采用該方法不需要額外的模型 參數(shù)。

[1] YANG Yunming, YU Haisui.Application of a non?coaxial soil model in shallow foundation[J]. Geomechanics and Geoengineering, 2006, 1(2): 139?150.

[2] GUTIERREZ M, ISHIHARA K. Non?coaxiality and energy dissipation in granular materials[J]. Soils and Foundations, 2000, 40(2): 49?59.

[3] CAI Yanyan, YU Haisui, WAMATOWSKI D, et al. Noncoaxial behavior of sand under various stress paths[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(8): 1509?1518.

[4] MIURA K, MIURA S, TOKI S. Deformation behavior of anisotropic dense sand under principal stress axes rotation[J]. Soils and Foundations, 1986, 26(1): 36?52.

[5] ISHIHARA K, TOWHATA I. Sand response to cyclic rotation of principal stress directions as induced by wave loads[J]. Soils and Foundations, 1983, 23(4): 11?26.

[6] SYMES M J, GENS A, HIGHT D W. Undrained anisotropy and principal stress rotation in saturated sand[J]. Geotechnique, 1984, 34(1): 11?27.

[7] SYMES M J, GENS A, HIGHT D W. Drained principal stress rotation in saturated sand[J]. Geotechnique, 1988, 33(4): 355?383.

[8] 嚴(yán)佳佳, 周建, 管林波, 等. 杭州原狀軟黏土非共軸特性及其影響因素研究[J]. 巖土工程學(xué)報(bào), 2013, 35(1): 96?102. YAN Jiajia, ZHOU Jian, GUAN Linbo, et al. Experimental study of non?coaxial property and influence facts of Hangzhou intact soft clay[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(1): 96?102.

[9] ZHOU Jian, YAN Jiajia, LIU Zhengyi, et al. Undrained anisotropy and non?coaxial behavior of clayey soil under principal stress rotation[J]. Journal of Zhejiang University: Science A, 2014, 26(4): 241?254.

[10] TONG Zhaoxia, ZHANG Jianmin, YU Yilin, et al. Drained deformation behavior of anisotropic sands during cyclic rotation of principal stress axes[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 31(6): 1509?1518.

[11] PRADEL D, ISHIHARA K, GUTIERREZ M. Yielding and flow of sand under principal stress axes rotation[J]. Soils and Foundations, 1990, 30(1): 87?99.

[12] TSUTSUMI S, HASHIGUCHI K. General non?proportional loading behavior of soils[J]. International Journal of Plasticity, 2005, 21(10): 1941?1969.

[13] 呂璽琳, 黃茂松, 錢建固. 基于非共軸本構(gòu)模型的砂土真三軸試驗(yàn)分叉分析[J]. 巖土工程學(xué)報(bào), 2008, 30(5): 646?651. Lü Xilin, HUANG Maosong, QIAN Jiangu. Bifurcation analysis in true trixial tests on sands based on non?coaxial elasto?plasticity model[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(5): 646?651.

[14] RUDNICKI J W, RICE J R. Conditions for the localization of the deformation in pressure sensitive dilatant materials[J]. Journal of the Mechanics and Physics of Solids, 1975, 23(6): 371?394.

[15] LASHKARI A, LATIFI M. A simple plasticity model for prediction of non?coaxial flow of sand[J]. Mechanics Research Communications, 2007, 34(1): 191?200.

[16] 錢建固, 黃茂松. 復(fù)雜應(yīng)力狀態(tài)下巖土體的非共軸塑性流動(dòng)理論[J]. 巖石力學(xué)與工程學(xué)報(bào), 2006, 25(6): 1259?1264. QIAN Jiangu, HUANG Maosong. Non?coaxial plastic flow theory in multi?dimensional stress state[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(6): 1259?1264.

[17] 扈萍, 黃茂松, 錢建固, 等. 砂土非共軸特性的本構(gòu)模擬[J]. 巖土工程學(xué)報(bào), 2009, 31(5): 793?798. HU Ping, HUANG Maosong, QIAN Jiangu, et al. Non?coaxial plasticity constitutive modeling of sands[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(5): 793?798.

[18] 董建國, 張啟輝, 李蓓, 等. 上海軟土的非共軸修正的彈塑性模型[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2001, 29(7): 774?778. DONG Jianguo, ZHANG Qihui, LI Bei, et al. Modified non?coaxial elastoplastic model of soft soil in Shanghai[J]. Journal of Tongji University(Natural Science), 2001, 29(7): 774?778.

[19] 沈揚(yáng). 考慮主應(yīng)力方向變化的原狀軟黏土試驗(yàn)研究[D]. 杭州: 浙江大學(xué)建筑工程學(xué)院, 2007: 40?82. SHEN Yang. Experimental study on effect of variation of principal stress orientation on undisturbed soft clay[D]. Hangzhou: Zhejiang University. School of Architectural and Civil Engineering, 2007: 40?82.

[20] 鄭鴻鑌. 主應(yīng)力軸旋轉(zhuǎn)下重塑黏土與原狀黏土特性試驗(yàn)研究[D]. 杭州: 浙江大學(xué)建筑工程學(xué)院, 2011: 52?54. ZHENG Hongbin. Experimental study of reconstituted clay and intact clay under principal stress rotation[D]. Hangzhou: Zhejiang University. School of Architectural and Civil Engineering, 2011: 52?54.

[21] ZDRAVKOVIC L, JARDINE R J. Some anisotropic stiffness characteristics of a silt under general stress conditions[J]. Geotechnique, 1997, 47(3): 407?437

[22] 劉元雪, 鄭穎人. 考慮主應(yīng)力軸旋轉(zhuǎn)對土體應(yīng)力應(yīng)變關(guān)系影響的一種新方法[J]. 巖土工程學(xué)報(bào), 1998, 20(2): 45?47. LIU Yuanxue, ZHENG Yingren. A new method to analyze the influence of principal stress axes rotation on the stress strain relation of soils[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(2): 45?47.

[23] BLANC M, BENEDETTO D, TIOUAJNI S. Deformation characteristics of dry Hostun sand with principal stress axes rotation[J]. Soils and Foundations, 2011, 51(4): 749?760.

Calculation of non-coaxial angle of soft clay under stress rotation

YAN Jiajia1, LIN Qinghui2, FU Liaoyi1

(1. Research Center of Coastal and Urban Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China; 2. Department of Architecture, Taizhou Vocational & Technical College, Taizhou 318000, China)

Three typical series of stress rotation path tests of soft clay with pure principal stress rotation, cyclic principal stress rotation and principal stress rotation combined with shearing were carried out on intact clay using a hollow cylinder apparatus in undrained condition. The direction of the strain increment induced by the principal stress rotation as well as the influence factors were studied in particular. The results show that the obvious deviation between the directions of principal strain increment (non?coaxial) and principal stress is observed in the three series of stress rotation path tests. The non?coaxial angle varies with principal stress rotation with a period of 90°, and the direction of the strain increment changes with the reverse of the principal stress rotation direction. The non?coaxial angle decreases with the increase of shear stress combined with stress rotation. The direction of stress rotation, path and shear stress level have significant impact on the non?coaxial behavior of soft clay, while the influence of the intermediate principal stress parameter is unobvious. The comparison between the calculation and test results indicates that the calculation of non?coaxial angle under stress rotation based on the stress mapping method proposed by Gutierrze can seize the non?coaxial behavior of soft clay subjected to principal stress rotation, and it also can reflect the influence of stress rotation direction, path and shear stress level.

principal stress rotation; soft clay; strain; non-coaxial

10.11817/j.issn.1672-7207.2016.09.028

TU443

A

1672?7207(2016)09?3117?08

2015?09?01；

2015?11?10

國家自然科學(xué)基金資助項(xiàng)目(51608477, 51338009)；廣東省院士工作站資助項(xiàng)目(2013B090400024)；同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室開放基金資助項(xiàng)目(KLE?TJGE?B1501) (Projects(51608477, 51338009) supported by the National Natural Science Foundation of China; Project(2013B090400024) supported by the Academician Workstation of Guangdong Province; Project(KLE?TJGE?B1501) supported by the Open Foundation of Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University)

嚴(yán)佳佳，助理研究員，從事巖土力學(xué)特性試驗(yàn)及理論研究；E-mail: yanjia.01@163.com

(編輯 劉錦偉)