姚仰平，萬(wàn) 征，楊一帆，牛 雷

（北京航空航天大學(xué) 交通科學(xué)與工程學(xué)院，北京 100191）

1 引 言

已有數(shù)據(jù)表明：溫度升高會(huì)改變巖土體的應(yīng)力重分布，進(jìn)而導(dǎo)致巖土體的變形，當(dāng)土體應(yīng)力減小到不足以平衡外載荷時(shí)便會(huì)產(chǎn)生破壞現(xiàn)象。這種在外部荷載保持不變而由于溫度變化所誘發(fā)的破壞可稱之為熱破壞現(xiàn)象。受溫度影響的地下巖土結(jié)構(gòu)廣泛地存在于各種地下設(shè)施中，如熱交換樁以及地下核廢料掩埋裝置等等。因此，研究溫度變化對(duì)于巖土的應(yīng)力-應(yīng)變關(guān)系的影響規(guī)律尤其重要。國(guó)外學(xué)者對(duì)熱破壞現(xiàn)象進(jìn)行了大量研究工作。Baldi[1]早在1988年就對(duì)熱破壞現(xiàn)象的細(xì)觀機(jī)制進(jìn)行了解釋說(shuō)明。Bourros[2]從變溫對(duì)于飽和黏土的抗剪強(qiáng)度的影響角度進(jìn)行了相應(yīng)的研究，此后，Jefferson[3]也對(duì)飽和黏土的不排水抗剪強(qiáng)度在變溫條件下進(jìn)行了相關(guān)研究。Hueckel[4]等對(duì)正常固結(jié)黏土以及超固結(jié)黏土在不排水變溫條件下的應(yīng)力-應(yīng)變特性進(jìn)行了比較系統(tǒng)的試驗(yàn)研究。

Hueckel針對(duì)飽和正常固結(jié) Boom黏土開(kāi)展了升溫下的三軸壓縮不排水試驗(yàn)工作。等向固結(jié)壓力值為5.75 MPa，試驗(yàn)結(jié)果如圖1、2所示。為考查溫度變化對(duì)于變形以及強(qiáng)度的影響，采用了對(duì)比試驗(yàn)，將飽和Boom黏土重復(fù)做了兩次試驗(yàn)，第1次是固定常溫下的常規(guī)不排水試驗(yàn)；第2次是常溫加載到一定值，然后做常載升溫試驗(yàn)。圖1、2中圓圈表示在室溫21 °C下所做的關(guān)于飽和正常固結(jié)土的不排水三軸壓縮試驗(yàn)結(jié)果，而正方形表示了加載-升溫的試驗(yàn)結(jié)果，在偏應(yīng)力值達(dá)到2 MPa之前，保持溫度在21 °C下做不排水三軸壓縮加載試驗(yàn)；當(dāng)偏應(yīng)力值達(dá)到2 MPa時(shí)，保持外載不變，此時(shí)將溫度由21°C升高到92 °C。由圖1的有效應(yīng)力路徑可見(jiàn)，在加載-升溫試驗(yàn)中，在偏應(yīng)力值達(dá)到 2 MPa前，應(yīng)力路徑與三軸壓縮不排水應(yīng)力路徑基本相同，圖2中孔壓也基本與常溫21 °C的試驗(yàn)孔壓相同。而當(dāng)做常載升溫的試驗(yàn)時(shí)，偏應(yīng)力值基本保持在2 MPa不變，而有效球應(yīng)力迅速減小，有效應(yīng)力路徑呈水平直線向左延伸，最終應(yīng)力狀態(tài)達(dá)到了臨界狀態(tài)。由圖1所示，升溫不僅改變了有效應(yīng)力路徑，同時(shí)，臨界狀態(tài)應(yīng)力比也有所提高，但在最終的溫度92 °C時(shí)，試驗(yàn)點(diǎn)的有效應(yīng)力狀態(tài)點(diǎn)達(dá)到了臨界狀態(tài)線，已經(jīng)導(dǎo)致了試樣的破壞。由此可見(jiàn)，單純地升溫過(guò)程同樣應(yīng)視為一種加載過(guò)程。由試驗(yàn)結(jié)果可見(jiàn)，升溫過(guò)程導(dǎo)致了排水條件下的體積收縮趨勢(shì)，使黏土顆粒間孔隙減小，但由于不排水邊界條件的限制，使得正孔壓持續(xù)增大，如圖2中正方形點(diǎn)線所示，在升溫過(guò)程中孔壓迅速增大，有效球應(yīng)力迅速降低，當(dāng)應(yīng)力比達(dá)到臨界狀態(tài)應(yīng)力比時(shí)，即使得土體達(dá)到破壞條件。

圖1 升溫對(duì)飽和黏土不排水路徑的影響Fig.1 Influence of heating on stress paths for saturated clay under undrained condition

圖2 升溫對(duì)飽和正常固結(jié)土的孔壓影響Fig.2 Influence of heating on pore pressure for saturated normal consolidated clay

姚仰平等[5－6]基于超固結(jié) UH 模型，提出了考慮溫度影響的UH模型[7]，也即熱UH模型。熱UH模型能否反映上述不排水條件下加載-升溫的應(yīng)力-應(yīng)變關(guān)系以及臨界狀態(tài)強(qiáng)度變化規(guī)律，以及能否基于熱 UH模型建立三軸壓縮不排水熱抗剪強(qiáng)度公式，這些都是本文所要探討的內(nèi)容。

2 熱UH模型

2.1 當(dāng)前屈服面與參考屈服面

由于熱UH模型是建立在UH模型基礎(chǔ)上，因此，在某一溫度下當(dāng)前應(yīng)力點(diǎn)A所在的屈服面為當(dāng)前屈服面，可表示為

式中：p、q分別為球應(yīng)力與剪應(yīng)力；px0為T(mén)0溫度下初始狀態(tài)時(shí)屈服面與 p軸右交點(diǎn)的橫坐標(biāo)；MT表示在某一溫度 T下的臨界狀態(tài)應(yīng)力比。cp=(λ-κ)/(1+e0)，其中λ為e-lnp坐標(biāo)系下正常固結(jié)線(NCL)的斜率；κ為回彈線的斜率；e0為初始孔隙比；H為硬化參數(shù)；ξT是溫度影響先期固結(jié)壓力系數(shù)，表示為溫度的對(duì)數(shù)函數(shù)：

式中：T為當(dāng)前溫度；T0為初始溫度；γ為材料參數(shù)，可以通過(guò)不同溫度下的等向壓縮試驗(yàn)確定。

某一溫度T下的先期固結(jié)壓力pxT與T0下的先期固結(jié)壓力px存在如下關(guān)系：

由式(2)、(3)可以看出，隨著溫度T的增大，其相應(yīng)的先期固結(jié)壓力隨之減小，這表明了溫度升高導(dǎo)致了超固結(jié)土一部分超固結(jié)度的喪失。參考屈服面可表示為

2.2 超固結(jié)參數(shù)RT

由于在任一溫度下，兩屈服面具有相同的MT，因此，屈服面幾何形狀相似，以當(dāng)前屈服面與參考屈服面的相似比來(lái)表示超固結(jié)參數(shù)RT，可表示為

將參考屈服面的表達(dá)式(4)變形得到：

則將式(6)代入到式(5)，得到：

超固結(jié)狀態(tài)下的RT是一個(gè)小于1的數(shù)，隨著加載的進(jìn)行RT值遞增，當(dāng)RT為1時(shí)到達(dá)臨界狀態(tài)。

2.3 臨界狀態(tài)強(qiáng)度

在式(1)中，MT為某一溫度 T下的臨界狀態(tài)應(yīng)力比。由Hueckel等的試驗(yàn)結(jié)果可知，升溫不僅會(huì)使正常固結(jié)線在e-lnp空間下移，導(dǎo)致孔隙比變小，還對(duì)臨界狀態(tài)時(shí)強(qiáng)度特性造成影響。對(duì)于飽和黏土的試驗(yàn)結(jié)果表明，升溫會(huì)提高臨界狀態(tài)應(yīng)力比。考慮到ξT是一種由溫度所導(dǎo)致的溫度超固結(jié)參數(shù)，則類比由應(yīng)力所形成的超固結(jié)參數(shù)R，仿照潛在強(qiáng)度的表達(dá)式形式構(gòu)建了受溫度超固結(jié)影響的臨界狀態(tài)強(qiáng)度表達(dá)式。受溫度影響的臨界狀態(tài)應(yīng)力比可表示為

式中：

式中：M0為初始溫度T0下的臨界狀態(tài)應(yīng)力比。隨著溫度的升高，臨界狀態(tài)強(qiáng)度值逐漸增大。這主要是由于溫度升高導(dǎo)致孔隙水活動(dòng)性增強(qiáng)，從而增加有效應(yīng)力，使土體變密實(shí)，內(nèi)摩擦角變大，提高了臨界狀態(tài)強(qiáng)度值。

2.4 熱彈性模量與泊松比

基于筆者等提出的熱UH模型，在文獻(xiàn)[7]中已經(jīng)仔細(xì)分析了Campanella和Mitchell得到的不同溫度下飽和伊利黏土在e-lnp坐標(biāo)系下的等向壓縮-回彈試驗(yàn)結(jié)果。由試驗(yàn)規(guī)律可知，不同溫度下回彈線基本保持平行狀態(tài)，只是隨著溫度升高，回彈線保持平行狀態(tài)向下平移，則對(duì)應(yīng)的回彈線斜率κ也應(yīng)完全相同，即不受溫度影響。而由κ所決定的彈性模量E也就不受溫度影響，仍與原超固結(jié)UH模型中彈性模量所用公式相同，泊松比ν也仍然作為一個(gè)常數(shù)使用，確定為0.3。

3 考慮溫度影響的三軸壓縮不排水抗剪強(qiáng)度表達(dá)式

基于熱UH模型，推導(dǎo)飽和黏土在不排水條件下受溫度影響的不排水抗剪強(qiáng)度表達(dá)式。其基本思路是根據(jù)屈服面方程，由體變?yōu)?以及臨界狀態(tài)強(qiáng)度特性條件，將不排水抗剪強(qiáng)度式由土性參數(shù)以及初始應(yīng)力條件表示出來(lái)。由于處于不排水條件下，因此體變?yōu)?0?？傮w變可分解為塑性體變與彈性體變之和：

則可利用參考屈服面中塑性體應(yīng)變與當(dāng)前屈服面中塑性體應(yīng)變相同這個(gè)特征，通過(guò)求取參考屈服面中的塑性體應(yīng)變得到其表達(dá)式。則參考屈服面方程經(jīng)過(guò)變形可得到：

又因?yàn)椋?/p>

聯(lián)立式(11)～(14)可得到：

由彈性部分：

考慮初始 K0固結(jié)狀態(tài)，則由當(dāng)前屈服面式(1)，以及K0的定義：

式中：σ110、σ220、σ330分別是K0固結(jié)時(shí)的大、中、小主應(yīng)力。

聯(lián)立式(10)、(15)、(16)、(19)得到：

在不排水條件下，且達(dá)到臨界狀態(tài)，則此時(shí)超固結(jié)應(yīng)力比參數(shù)RT的值為1，且應(yīng)滿足下述條件：

將式(21)、(22)代入到式(20)中，得到：

將式(19)代入到式(23)中，得到：

因?yàn)樵谂R界狀態(tài)，滿足：

不排水抗剪強(qiáng)度su在三軸壓縮下表示為

式(27)為可考慮超固結(jié)度升溫條件下的不排水抗剪強(qiáng)度表達(dá)式。下面分幾種情況進(jìn)行討論：

①當(dāng)溫度為常溫T0時(shí)，則ξT=1，若土體處于等向正常固結(jié)狀態(tài)，則RT0=1，K0=1，則表達(dá)式為

式(28)即退化為Wood所提出的抗剪強(qiáng)度表達(dá)式；

②當(dāng)溫度保持為常溫T0時(shí)，若考慮超固結(jié)度的影響，則式(27)可表示為

式(29)即為只考慮超固結(jié)度影響的不排水抗剪強(qiáng)度表達(dá)式，其中，M0和R0分別為T(mén)0溫度下的臨界狀態(tài)應(yīng)力比以及初始超固結(jié)應(yīng)力比參數(shù)；

③當(dāng)土體為正常固結(jié)狀態(tài)，溫度由常溫T0升高到T時(shí)，則式(27)可寫(xiě)為

式(30)即為只考慮溫度升高時(shí)正常固結(jié)土的不排水抗剪強(qiáng)度表達(dá)式。

采用Kuntiwattanakul等[8]關(guān)于一種黏土的不排水強(qiáng)度試驗(yàn)結(jié)果來(lái)模擬所提出的抗剪強(qiáng)度公式。由于文獻(xiàn)并未提供相關(guān)參數(shù)，因此，可采取表1所示的參數(shù)值來(lái)進(jìn)行計(jì)算，利用所式(27)來(lái)對(duì)公式進(jìn)行功能上的模擬，初始等向固結(jié)球應(yīng)力為196 kPa。

表1 黏土材料參數(shù)Table 1 Material parameters for clay

如圖 3所示，分別采用OCR為1、2.2、4、8的黏土在不同設(shè)定溫度下進(jìn)行不排水剪切測(cè)試，由圖可見(jiàn)，隨著溫度的升高，抗剪強(qiáng)度呈增大趨勢(shì)變化，隨超固結(jié)度的提高，抗剪強(qiáng)度也隨之提高。利用式(3)、(13)分別計(jì)算出ζT和MT，由初始K0條件代入到式(32)中，得到在不同超固結(jié)度下升溫過(guò)程的強(qiáng)度值。圖4為采用抗剪強(qiáng)度式(27)模擬的結(jié)果，不僅能夠考慮超固結(jié)度提高抗剪強(qiáng)度提高的規(guī)律，而且也隨著溫度的提高抗剪強(qiáng)度也相應(yīng)的提高?？梢?jiàn)，所提公式可簡(jiǎn)單、合理地反映溫度以及超固結(jié)度對(duì)于不排水抗剪強(qiáng)度的影響規(guī)律。

圖3 溫度與不同超固結(jié)度黏土抗剪強(qiáng)度的關(guān)系試驗(yàn)結(jié)果Fig.3 Relationships between temperature and shear strength under undrained condition for clay with various OCRs

圖4 溫度與不同超固結(jié)度黏土抗剪強(qiáng)度的關(guān)系模擬Fig.4 Simulation of relationships between temperature and shear strength under undrained condition for clay with various OCRs

4 熱UH模型不排水剪切分析及模擬

由于在某一固定溫度下超固結(jié)土體的變形以及強(qiáng)度特性與超固結(jié)UH模型所描述的相類似，因此，本節(jié)著重討論在固定應(yīng)力狀態(tài)下，溫度的升高對(duì)于土體的變形以及強(qiáng)度等方面的影響。

4.1 三軸不排水常載升溫過(guò)程

如圖5所示，在某一溫度T0下，考慮正常固結(jié)黏土在不排水條件下由A點(diǎn)加載到B點(diǎn)，然后保持外加載荷不變，保持不排水條件，則此時(shí)升高溫度，由T0持續(xù)升溫，當(dāng)升高到T1時(shí)，則達(dá)到C點(diǎn)。在初期不排水加載過(guò)程中，應(yīng)力路徑如AB所示，在升溫過(guò)程中，由于保持外載不變，因此，剪應(yīng)力保持不變，而在點(diǎn)劃線所示的等q平面內(nèi)，應(yīng)力路徑如圖BC所示，升溫過(guò)程造成正塑性體應(yīng)變產(chǎn)生，由于不排水條件的限制，因此，產(chǎn)生正孔壓，造成有效應(yīng)力減小，球應(yīng)力從 pB減小到 pC，此時(shí)，應(yīng)力狀態(tài)達(dá)到臨界狀態(tài)。對(duì)于超固結(jié)黏土，則在升溫過(guò)程中造成超固結(jié)度的損失，使超固結(jié)土逐步還原為正常固結(jié)土，同時(shí)，臨界狀態(tài)強(qiáng)度從 M0增大為MT。對(duì)于某些對(duì)溫度敏感的黏土類型，比如ξT會(huì)隨溫度增大而迅速減小，則使超固結(jié)性很快退化為正常固結(jié)土，使應(yīng)力狀態(tài)很快達(dá)到臨界狀態(tài)。

圖5 三軸不排水常溫加載-常載升溫過(guò)程Fig.5 Process of loading with constant temperatureheating with constant loading under triaxial compression undrained condition

與上述在p-q空間中的加載過(guò)程相對(duì)應(yīng)，如圖6所示，在e-lnp空間中，當(dāng)溫度在某一溫度T0下時(shí)，當(dāng)進(jìn)行不排水三軸壓縮時(shí)，應(yīng)力狀態(tài)從A點(diǎn)達(dá)到B點(diǎn)，此時(shí)，正常固結(jié)線與臨界狀態(tài)線相互平行。達(dá)到B點(diǎn)后，從T0開(kāi)始升溫到達(dá)T1，則對(duì)于正常固結(jié)土而言，若繼續(xù)升溫，則應(yīng)力狀態(tài)點(diǎn)從B點(diǎn)平移到達(dá)C點(diǎn)，在升溫過(guò)程中，臨界狀態(tài)線向下平行移動(dòng)到T1點(diǎn)。若此時(shí)狀態(tài)點(diǎn)與臨界狀態(tài)線相交，則表明土體達(dá)到臨界狀態(tài)。

圖6 常溫加載然后升溫常載過(guò)程的不排水應(yīng)力路徑Fig.6 Stress path in e-lnp space under undrained condition with constant loading-heating process

4.2 三軸不排水模擬

在常溫剪切至某一預(yù)定值后，外載不變、升溫不排水的過(guò)程已有了相關(guān)的測(cè)試結(jié)果，采用Hueckel等[9－10]分別在1987、1990年對(duì)正常固結(jié)黏土的測(cè)試結(jié)果，以及Hueckel等[11]和Del Olmo等[12]關(guān)于超固結(jié)黏土的測(cè)試結(jié)果，試驗(yàn)條件如引言中所述。對(duì)試驗(yàn)過(guò)程進(jìn)行模擬，由于原文獻(xiàn)并未提供所需的模型參數(shù)數(shù)值，因此，為便于熱UH模型對(duì)Boom黏土在不排水常載升溫下的應(yīng)力-應(yīng)變關(guān)系進(jìn)行功能上的演示，根據(jù)黏土的常用取值范圍，選取模型參數(shù)如表2所示。

表2 Boom黏土材料參數(shù)Table 2 Material parameters for Boom clay

利用所提的熱UH模型模擬上述不排水升溫加載過(guò)程，強(qiáng)度結(jié)果如圖7(a)中的M21與M92所示，顯示了升溫過(guò)程中臨界狀態(tài)強(qiáng)度提高的過(guò)程。由應(yīng)力路徑的模擬結(jié)果可見(jiàn)，黑實(shí)線與正常固結(jié)土常溫不排水模擬基本一致。對(duì)于正常固結(jié)土的常溫加載常載升溫的模擬如圖中虛線所示，對(duì)于升溫過(guò)程中球應(yīng)力沿水平直線減小的規(guī)律能夠明顯反映出來(lái)，并最終達(dá)到臨界狀態(tài)，使土體產(chǎn)生了破壞作用。圖7(b)的縱坐標(biāo)表示孔壓用初始平均應(yīng)力來(lái)歸一化，正常固結(jié)土在常溫下不排水剪切與試驗(yàn)結(jié)果基本吻合，常溫加載升溫常載模擬的軸向應(yīng)變偏小，但反映出升溫后孔壓增大的規(guī)律。對(duì)于超固結(jié)土也存在相似的規(guī)律，圖8(a)所示的是超固結(jié)度為3的Boom黏土的不排水試驗(yàn)路徑，在偏應(yīng)力達(dá)到0.9 MPa時(shí)開(kāi)始升溫過(guò)程，圖8(b)為相對(duì)應(yīng)的模擬結(jié)果，由此可見(jiàn)，對(duì)于超固結(jié)土的常載升溫應(yīng)力-應(yīng)變關(guān)系也能得到合理的模擬。

圖7 模擬結(jié)果Fig.7 Simulation results

圖8 超固結(jié)土試驗(yàn)與模擬結(jié)果Fig.8 Test and simulation results for overconsolidated clay

5 不排水條件下熱應(yīng)力-應(yīng)變特性分析

為了便于分析不同應(yīng)力比條件下升溫過(guò)程對(duì)于變形和強(qiáng)度特性的影響以及超固結(jié)度對(duì)升溫不排水過(guò)程的作用，選取一組參數(shù)，對(duì)應(yīng)力比條件、超固結(jié)度、先期固結(jié)壓力等因素對(duì)于升溫不排水的變形與強(qiáng)度特性做一比較。選取參數(shù)如表3所示。

表3 黏土材料參數(shù)Table 3 Material parameters for clay

初始等向平均應(yīng)力為200 kPa，設(shè)置常溫條件為20 °C。在20 °C條件下的臨界狀態(tài)強(qiáng)度記為M20，溫度在90 °C條件下的臨界狀態(tài)強(qiáng)度記為M90。則如圖9所示，在20 °C條件下做不排水三軸剪切至應(yīng)力比分別為0.2、0.4、0.6、0.8時(shí)保持外載不變。此時(shí)再做載荷不變，升溫到90 °C的升溫控制。

模擬結(jié)果如圖 9(a)所示，從下到上依次為升溫過(guò)程的有效應(yīng)力路徑，升溫過(guò)程中，球應(yīng)力迅速減小，而剪應(yīng)力保持不變，同時(shí)臨界狀態(tài)強(qiáng)度從 M20增大到 M90。當(dāng)應(yīng)力比設(shè)定值較大時(shí)，如 0.6、0.8時(shí)都使應(yīng)力比達(dá)到了臨界狀態(tài)，形成了破壞現(xiàn)象。圖9(b)為相應(yīng)的偏應(yīng)力與偏應(yīng)變關(guān)系。從圖9(c)可看出，在應(yīng)力比較小時(shí)若開(kāi)始升溫，則其剪切模量較常溫時(shí)呈增大規(guī)律，偏應(yīng)變較常溫時(shí)呈減小規(guī)律。應(yīng)力比越大或者剪應(yīng)力越大，則升溫對(duì)其剪應(yīng)變有加速增大趨勢(shì)。

圖9 不同應(yīng)力比下升溫不排水應(yīng)力-應(yīng)變特性模擬Fig.9 Simulation of heating under undrained condition with various stress ratios

當(dāng)設(shè)定應(yīng)力比為 0.7時(shí)，則在應(yīng)力比小于 0.7時(shí)保持20 °C條件下的不排水三軸剪切過(guò)程，直到應(yīng)力比達(dá)到0.7。當(dāng)達(dá)到0.7時(shí)，保持外載不變，并實(shí)施升溫過(guò)程。由圖10(a)可見(jiàn)，在溫度依次升高到40 °C、60 °C、80 °C、90 °C 時(shí)，臨界狀態(tài)強(qiáng)度逐漸增大，有效應(yīng)力持續(xù)沿水平方向向左減小，在溫度達(dá)到90 °C時(shí)，達(dá)到了臨界狀態(tài)，構(gòu)成了熱破壞現(xiàn)象。圖10(b)所示為相應(yīng)的偏應(yīng)力隨剪應(yīng)變的發(fā)展過(guò)程，圖10(c)為相應(yīng)的應(yīng)力比偏應(yīng)變關(guān)系。

圖10 應(yīng)力比為0.7下升溫不排水應(yīng)力-應(yīng)變特性模擬Fig.10 Simulation of heating under undrained condition with constant stress ratio 0.7

為便于分析超固結(jié)度對(duì)于升溫不排水的影響規(guī)律，只將參數(shù)γ由0.6設(shè)定為0.15，其他參數(shù)不變。當(dāng)設(shè)定不同先期固結(jié)壓力時(shí)，分別制取先期固結(jié)壓力為200、400、800、1 600 kPa的式樣，然后回彈至平均應(yīng)力為200 kPa。仍然設(shè)定控制應(yīng)力比為0.7，在達(dá)到0.7后保持外載不變并做升溫的模擬。其應(yīng)力路徑如圖11(a)所示，當(dāng)超固結(jié)度較小時(shí)，則有效應(yīng)力減小迅速，應(yīng)力比也偏小，隨著超固結(jié)度的增大，則有效應(yīng)力減小的幅度很小，同時(shí)，剪應(yīng)變也在升溫過(guò)程中變化很小。其偏應(yīng)力、應(yīng)力比與偏應(yīng)變關(guān)系分別如圖11(b)、(c)所示。

圖11 不同先期固結(jié)壓力下升溫不排水應(yīng)力-應(yīng)變特性模擬Fig.11 Simulation of heating under undrained condition with various preconsolidation pressures

在相同先期固結(jié)壓力，不同的當(dāng)前平均應(yīng)力下做不排水剪切的模擬結(jié)果如圖12所示，圖12(a)為初始平均應(yīng)力分別為1 600、800、400、200 kPa下的常溫剪切到應(yīng)力比為0.7，然后外載不變升溫的應(yīng)力路徑模擬結(jié)果。控制應(yīng)力比仍設(shè)定為 0.7，由圖可以看出，影響應(yīng)力路徑的仍然是超固結(jié)度，而先期固結(jié)壓力對(duì)升溫不排水的變形以及強(qiáng)度特性影響很小，圖 12(b)為對(duì)應(yīng)的偏應(yīng)力與偏應(yīng)變關(guān)系，圖12(c)也顯示了同圖 11(c)相同的受超固結(jié)度影響的規(guī)律。

圖12 相同先期固結(jié)壓力下升溫不排水應(yīng)力-應(yīng)變特性模擬Fig.12 Simulation of heating under undrained condition with constant preconsolidation pressure

6 結(jié) 論

（1）根據(jù)已有的試驗(yàn)結(jié)果可知，在不排水條件下，升溫會(huì)造成土體的孔壓持續(xù)增大，使有效球應(yīng)力迅速減小，并最終造成熱破壞現(xiàn)象，因此，不排水條件下的升溫過(guò)程應(yīng)視為一種加載過(guò)程。

（2）根據(jù)所提的熱UH模型，由不排水條件推導(dǎo)得到了能夠反映溫度影響的不排水抗剪強(qiáng)度表達(dá)式，并通過(guò)相關(guān)試驗(yàn)結(jié)果的規(guī)律驗(yàn)證了公式的合理性。

（3）由所提的熱UH模型成功模擬了熱破壞過(guò)程中土體的應(yīng)力-應(yīng)變關(guān)系以及臨界狀態(tài)強(qiáng)度變化規(guī)律。

[1]BALDI G, HUECKEL T, PELLEGRINI R. Thermal volume change of the mineral-water system in low-porosity clay soils[J]. Canadian Geotechnical Journal, 1988, 25(4): 807－825.

[2]BOURROS C M. The effect of temperature change on consolidation and shear strength of saturated cohesive soils[D]. USA: University of Washington, 1973.

[3]JEFFERSON I, ROGERS C D F, SMALLEY I J.Discussion: ‘Temperature effects on undrained shear characteristics of clay’ by Kuntiwattanakul et al[J].Soils and Foundations, 1996, 36(3): 141－143.

[4]HUECKEL T, PEANO A, PELLEGRINI R. A constitutive law for thermo-plastic behaviour of rocks: An analogy with clays[J]. Surveys in Geophysics, 1993,15(5): 643－671.

[5]姚仰平, 李自強(qiáng), 侯偉, 等. 基于改進(jìn)伏斯列夫線的超固結(jié)土本構(gòu)模型[J].水利學(xué)報(bào), 2008, 39(11): 1244－1250.YAO Yang-ping, LI Zi-qiang, HOU Wei, et al.Constitutive model of over-consolidated clay based on improved Hvorslev envelope[J]. Journal of Hydraulic Engineering, 2008, 39(11): 1244－1250.

[6]姚仰平, 馮興, 黃祥, 等. UH模型在有限元分析中的應(yīng)用[J].巖土力學(xué), 2010, 31(1): 237－245.YAO Yang-ping, FENG Xing, HUANG Xiang, et al.Application of UH model to finite element analysis[J].Rock and Soil Mechanics, 2010, 31(1): 237－245.

[7]姚仰平, 楊一帆, 牛雷. 考慮溫度影響的 UH模型[J].中國(guó)科學(xué)(E輯), 2011, 41(2): 158－169.YAO Yang-ping, YANG Yi-fan, NIU Lei. UH model considering temperature effect[J]. Science in China(Series E), 2011, 41(2): 158－169.

[8]PISIT K, IKUO T, KANTA O, et al. Temperature effects on undrained shear characteristics of clay[J]. Soils and Foundations, 1995, 35(1): 147－162.

[9]HUECKEL T, BORSETTO M, PEANO A. Modelling of coupled thermo-elasto-plastic-hydraulic response of clays subjected to nuclear waste heat[M]. Chichester: John Wiley, 1987: 213–235.

[10]HUECKEL T, BALDI G. Thermoplasticity of saturated clays: Experimental constitutive study[J]. Journal of Geotechnical Engineering, 1990, 116(12): 1778－1796.

[11]HUECKEL T. Water-mineral interaction in hydromechanics of clays exposed to environmental loads:A mixture-theory approach[J]. Canadian Geotechnical Journal, 1992, 29(6): 1071－1086.

[12]DEL OLMO C, FIORAVANTE V, GERA F, et al.Thermomechanical properties of deep argillaceous formations[J]. Engineering Geology, 1996, 41(1－4): 87－102.