土體熱力學性質及本構關系研究綜述

中圖分類號：TU432 文獻標志碼：A 文章編號：2096-6717（2025）03-0021-15

Review on thermomechanical behaviors and constitutive relationships of soil

LIU Hong ^1，2a ， ， XIAO Yang2a，2b，WU Huanran ^2a，2b （1.School of Architectural Engineering，Neijiang Normal University，Neijiang 64l112，Sichuan，P.R.China; 2a.Key Laboratory of New Technology forConstruction of Cities in Mountain Area，Ministry of Education; 2b. School of Civil Engineering，Chongqing University，Chongqing 4OoO45，P.R. China）

Abstract：Exploitation of the shallow geothermal energy has gradually become a key topic in geotechnical engineering worldwide with the rapid development of renewable energy and geothermal heat pump techniques. However，the theory of the shallow geothermal energy falls far behind its application.In particular，the mechanism of soil under the complex coupling of stressand temperature fields remains unclear.For a single soil such as sand，clay，and silt， some results have been found by carrying out the temperature-controlled experiments and proposing the related thermal constitutive models. To promote more comprehensive and indepth thermal behavior of soil and contribute to more practical application of the thermal constitutive model in energy engineering，this paper firstly summarizes the experimental research on the thermomechanical properties of soil.Then，research progress and status quo of thermomechanical constitutive relations of soil，based on different theoretical frameworks，arereviewed in detail.Thereafter，theapplication of these thermomechanical constitutive relations in engineering practiceis introduced briefly.Finally，suggestions in further research are provided，and the development trendof the thermomechanical constitutive relationshipof soil is discussed and prospected，in view of current existing problems.

Keywords： geothermal energy； geothermal heat pump； stress fields；temperature fields； constitutive relation;temperature-controlled test

近年來，清潔能源的高效率開發(fā)逐漸成為能源開發(fā)進程中的重中之重。作為清潔能源的一種，淺層地熱能的開發(fā)利用自然成為關注的重點課題[2-4]。淺層地熱能是指通過地源熱泵換熱技術，利用蘊藏在地表以下 200m 以內、溫度低于 25^°C 的熱能。其主要蘊藏在地表以下一定深度范圍內的土體、地下水和地表水中。利用地源熱泵技術，在夏季，可實現(xiàn)將建筑物或構筑物中的熱源供給溫度較低的土體，達到夏季制冷的效果；在冬季，可將地表以下土體中的熱源供給地面的建筑物或構筑物中，達到冬季供暖的效果，從而使建筑物或構筑物達到“冬暖夏涼”，功能類似于空調，且該方法具有環(huán)保無污染、耗能低、可循環(huán)利用等優(yōu)點。另外，淺層地熱能還被用于農業(yè)中的溫室大棚，養(yǎng)殖業(yè)中的水產養(yǎng)殖，旅游業(yè)中的溫泉療養(yǎng)，工業(yè)中的產品加工（類似于鍋爐），電力行業(yè)中的地熱發(fā)電等方面。中國淺層地熱能分布廣泛，發(fā)展?jié)摿薮?，為了更有效地開發(fā)淺層地熱能，越來越多的學者開始關注淺層地熱能開發(fā)過程中溫度對周圍土體力學性質的影響。此外，隨著核廢料處置[5]能源樁[6-8]跨越寒區(qū)和酷暑區(qū)的高速公路[9-11]、儲熱結構[12-15]等的不斷發(fā)展和廣泛應用，溫度對土體工程性質的影響逐漸成為巖土工程界的研究熱點。

對于土體熱力學性質的研究，通常包括試驗研究、理論研究和數(shù)值研究三方面。通過開展考慮溫度影響的室內試驗，即溫控試驗，探索溫度場和應力場耦合作用下土體的應力變形特性，并將其與常規(guī)應力場中土體的應力特性進行比較分析，得出土體的熱力學特性，探索溫度對土體變形、強度、剪脹方程、屈服方程、臨界狀態(tài)等方面的影響[16-20]。應用彈塑性理論，基于一定的假定，結合溫控試驗發(fā)現(xiàn)的土體應力變形規(guī)律，開展理論研究，提出合理的土體熱力學本構關系[21-24]。利用有限元軟件的二次開發(fā)平臺，將熱力學本構關系應用于有限元軟件中[25-27]，并對縮尺的模型試驗[28-31]或足尺的現(xiàn)場試驗[32-34]進行建模分析，通過數(shù)值分析結果和試驗監(jiān)測結果的對比分析，驗證數(shù)值分析方法的可靠性和準確性。最終，將數(shù)值分析方法應用于實際能源工程的設計和施工中，預測能源工程開發(fā)過程中可能出現(xiàn)的危險點或薄弱點，提出應對措施，提高能源結構的安全性，從而減少安全事故的發(fā)生，提高能源結構的開發(fā)效率，為淺層地熱能的開發(fā)利用奠定堅實的理論基礎。

基于近年來對土體熱力學本構關系的研究，結合筆者的部分研究成果，綜述土體的溫控試驗、基于不同理論框架建立的熱力學本構關系以及本構關系在數(shù)值分析中的應用等方面的主要進展，分析目前土體熱力學本構關系仍需進一步突破的方向，為淺層地熱能、能源樁等溫度相關實際工程提供重要的理論依據(jù)。

1土體熱力學試驗研究

土體熱力學本構關系的提出需以土體的應力-應變關系為研究對象，在大量溫控試驗的基礎上，得出土體的屈服特性、剪脹特性、固結特性等隨應力和溫度的演變規(guī)律，在彈塑性理論框架下，基于一定的基本假設，建立能夠合理描述土體的基本特性、亞基本特性和關聯(lián)基本特性[35的數(shù)學表達形式。因此，在提出熱力學本構關系之前，開展一系列溫控試驗，總結土體的熱力學性質是不可或缺的一環(huán)。

為明確溫度對土體本構關系的影響，確定本構模型相關參數(shù)，進而建立可靠的土體熱力學本構模型，不少學者在常規(guī)杠桿式固結儀、四聯(lián)直剪儀、靜三軸儀、動三軸儀、空心扭剪儀和真三軸等室內試驗儀器的基礎上，設置螺旋型加熱管、熱電偶、水浴箱、不銹鋼金屬壓力罩和隔熱罩等設備，賦予常規(guī)試驗儀器加熱和溫度監(jiān)測功能，從而能夠進行一維[16]或三維[20]的溫控固結試驗、溫控直剪試驗[36-38]溫控靜三軸試驗[39-44]、溫控真三軸試驗[45]、溫控動三軸試驗[46和溫控空心扭剪試驗[47-49]。

通過一維溫控固結試驗，Abuel-Naga等[6]發(fā)現(xiàn)，飽和黏土的軸向變形隨溫度的增加而線性增加。通過三維溫控固結試驗， Ng 等[20]發(fā)現(xiàn)，在熱固結過程中，飽和砂土發(fā)生的體積變形與砂土的相對密實度和有效應力有關：有效應力一定時，隨著相對密實度的增加，試樣逐漸從熱壓縮向熱膨脹轉化；相對密實度一定時，隨著有效應力的增加，試樣從熱膨脹轉變?yōu)闊釅嚎s。因此，在建立飽和砂土的熱力學本構模型時，應考慮試樣初始相對密實度、有效應力和溫度的耦合作用。

通過溫控直剪試驗，李春紅等[36-38]發(fā)現(xiàn)，樁-土接觸面的剪切強度低于周圍土體的剪切強度，且法向應力越大，該差值越明顯。此外，隨著溫度的升高，剪切強度略有增大。相對于樁基和土體而言，樁基與土體的接觸面是最薄弱的環(huán)節(jié)，在之后的本構理論和工程應用研究中，應給予樁-土接觸面足夠的重視。

通過溫控靜三軸試驗，郭楨等[50]發(fā)現(xiàn)，飽和黏土的剪切強度隨著圍壓的增加而增加，隨著溫度的增加而降低，應力-應變曲線呈典型的雙曲線特征。通過摩爾應力圓得出，其黏聚力和內摩擦角隨著溫度的增加而降低，呈現(xiàn)熱軟化特征，且溫度效應與圍壓大小有關，圍壓越大時，軟黏土的溫度效應越明顯。祁良等[51]發(fā)現(xiàn)，在不排水剪切條件下，孔隙水壓力隨著圍壓的增大而增大，隨著溫度的增大而減小。此外，臨界狀態(tài)下的應力比隨著溫度的增加呈線性增加趨勢。Cekerevac等[17]發(fā)現(xiàn)，飽和黏土的熱體積變形與土體的應力歷史有關，即與超固結比有關，隨著超固結比的增加，飽和黏土的熱體積變形從熱壓縮變形逐漸向熱膨脹變形過渡。此外，還發(fā)現(xiàn)飽和黏土的屈服應力與溫度有關，且隨著溫度的升高呈對數(shù)形式降低，但正常固結線（NCL）和臨界狀態(tài)線（CSL）的斜率與溫度無關。通過溫控真三軸試驗，潘旋[45]發(fā)現(xiàn)，隨著溫度的降低，鈣質黏土的應變硬化現(xiàn)象更明顯，土體的強度增加。

考慮到實際能源工程中經(jīng)常遇到車輛荷載、地震荷載等循環(huán)荷載的作用，不少研究者通過開展飽和與非飽和狀態(tài)下的溫控動三軸試驗，發(fā)現(xiàn)非飽和王體的累積塑性變形與溫度和基質吸力有關，且基質吸力一定時，累積塑性變形隨著溫度的增加而增加，當溫度一定時，累積塑性變形隨著基質吸力的增加而減小[46]?？紤]到尺寸效應的影響，Liu等[47-49]進行了一系列溫控空心扭剪試驗，發(fā)現(xiàn)密砂的熱體積變形僅與溫度有關，與有效應力無關，且隨著溫度的增加呈線性增加趨勢，不排水剪切強度隨著有效應力的增加而增加，隨著溫度的增加而線性減小，呈現(xiàn)出熱軟化現(xiàn)象，臨界狀態(tài)線未發(fā)現(xiàn)明顯變化。

2 土體熱力學本構關系研究

本構關系又稱為土體的力學本構方程，或土體的應力-應變模型，是描述土體力學特性的數(shù)學表達式。通過開展溫控試驗，學者們對土體在不同工況下的熱力學性質有所認識，結合溫控試驗結果，提出力學或熱屈服面方程、剪脹方程、臨界狀態(tài)方程等，隨后，基于不同的彈塑性理論框架，建立考慮溫度效應的土體本構關系。

2.1 Duncan-Chang模型

作為非線性彈性模型的典型代表，鄧肯-張模型包括 E-2 和 E?B 兩種模型，其模型參數(shù)均可由常規(guī)三軸試驗結果確定。

郭楨等[50通過溫控三軸不排水剪切試驗數(shù)據(jù)，分別得出不同溫度作用下的鄧肯-張模型參數(shù)，再將參數(shù)代入應力-應變關系式，則得到考慮溫度效應的鄧肯-張非線性彈性模型，并將該模型預測曲線與試驗數(shù)據(jù)進行比較分析。在Duncan-Chang模型的基礎上，謝云等[52考慮基質吸力的影響，建立常溫下非飽和土的非線性本構方程，通過非飽和土的溫控試驗數(shù)據(jù)擬合，提出考慮溫度影響的非飽和土總黏聚力方程，從而得到考慮溫度影響的非飽和土強度方程。結合試驗數(shù)據(jù)，得出切線體積模量隨溫度和基質吸力的變化規(guī)律，提出考慮溫度和基質吸力影響的非飽和膨脹土的非線性模型。當溫度為常溫時，該模型可以退化為常溫下非飽和土的非線性模型，當基質吸力為O時，該模型可退化為Duncan-Chang雙曲線模型。Liu等[53]從雙曲線模型[54和指數(shù)模型[55-57]的切線模量出發(fā)，引人參數(shù)θ，建立冪函數(shù)模型的切線模量一般表達式，同時引入摩爾-庫倫強度準則，提出適用于黏土的常溫冪函數(shù)型應力-應變關系式。依據(jù)寧波軟黏土[50]、Bourke粉土[58]、瀝青瓦混合物[59]、高嶺土[等不同類型和性質的土體溫控三軸試驗數(shù)據(jù)，分別建立溫度與該冪函數(shù)模型中的黏聚力、內摩擦角、冪指數(shù)和切線體積模量的關系式，并進行歸一化處理，從而得到考慮溫度影響的適用于黏土的冪函數(shù)型非線性彈性模型。

2.2 Cam-Clay模型

英國劍橋大學的Roscoe等[6提出了適用于正常固結飽和黏土的經(jīng)典Cam-Clay模型，因其模型參數(shù)均有明確的物理意義，且可由常規(guī)室內三軸試驗進行確定，從而廣泛應用于巖土工程中。

Hueckel等1通過三軸試驗，探究常溫下碳酸鹽含量對飽和黏土峰值強度或表觀最大預壓各向同性應力以及剪切破壞的影響。發(fā)現(xiàn)飽和黏土的熱彈性變形主要受圍壓影響，與碳酸鹽含量無明顯關系；但其熱塑性變形則需考慮碳酸鹽的影響，碳酸鹽含量越高，土體熱軟化現(xiàn)象越不明顯。此外，通過溫控試驗，還研究了溫度變化過程中飽和黏土的強度、溫度引起的超固結現(xiàn)象、剪切破壞以及不排水剪切過程中孔隙水壓力的變化規(guī)律。依據(jù)試驗數(shù)據(jù)，在劍橋模型的基礎上，提出了一種考慮碳酸鹽含量的飽和黏土熱彈塑性本構模型。劉紅等[62通過不同溫度作用下飽和黏土的壓縮曲線發(fā)現(xiàn)，雖然熱固結和力學固結的作用機理不同，但熱固結階段和常溫下力學固結階段發(fā)生的變形量相同。因此，將土體的熱力學固結過程等效為純力學固結過程，在修正劍橋模型的基礎上，通過推導純力學固結過程中彈塑性體積變形，結合非關聯(lián)流動法則，推導王體的塑性剪切變形，從而得出土體總的體積和剪切變形，最終提出考慮溫度影響的彈塑性本構模型。

在修正劍橋模型的基礎上，在臨界土力學理論框架下，Abuel-Naga等[63]得出常溫下土體的屈服函數(shù)和塑性勢函數(shù)分別為

式中： ?_c 為各向同性先期固結應力； η 為應力比，即η=dq/dp;ψ 為塑性勢函數(shù)； dε_v^p 和 dε_s^p 分別為塑性體積變形增量和塑性剪切變形增量； α 和 β 為材料參數(shù)。通過溫控試驗數(shù)據(jù)，建立溫度與材料參數(shù) α 和 β 的關系式，從而提出一種適用于飽和黏土的熱彈塑性本構模型。該模型中的材料參數(shù) α 隨著溫度的增加而線性增加，當 α 為0時，該屈服面退化為修正劍橋模型的屈服面，該模型退化為不考慮溫度影響的修正劍橋模型，如圖1所示。

Cui等4通過分析已有的溫控試驗數(shù)據(jù)發(fā)現(xiàn)，土體壓縮曲線斜率與溫度無關，先期固結壓力隨著溫度增加呈非線性降低趨勢，從而導致土體彈性區(qū)域減少，而隨著超固結比（OCR）的增加，土體的熱體積變形從熱壓縮向熱膨脹過渡。隨后，在修正劍橋模型的基礎上，分別建立了力學固結和熱固結屈服函數(shù)（如圖2所示），通過三維空間化，利用關聯(lián)流動法則，提出了一種考慮應力歷史和溫度影響的適用于飽和黏土的熱力學本構模型。

在Cui等64的熱本構模型基礎上，利用雙屈服面塑性理論[66]，Hong等[65]提出了一種可以描述天然飽和黏土的有限熱力學彈性區(qū)域以及彈性特性向塑性特性平穩(wěn)轉換的熱力學模型。該模型引入了兩種塑性機制：一種是再現(xiàn)高超固結比下飽和黏土發(fā)生膨脹和收縮的熱塑性特性，另一種是描述溫度對屈服行為的影響。此外，該模型采用與上述兩種塑性機制相關的附加屈服面，即內屈服面，用于考慮現(xiàn)有熱屈服面內部的塑性行為，同時定義了與內屈服面相關的漸進塑性硬化機制，使塑性模量在屈服面內的熱加載過程中平滑變化。Tang等運用Alonso模型描述常溫下非飽和膨脹土的力學特性，利用Cui等4的模型描述飽和土體的熱力學特性，結合這兩種模型，運用雙層結構的概念，將土體劃分為微觀結構水平（即由各向異性團聚體組成的黏土顆粒水平）和宏觀結構水平（即更大尺寸的土體結構水平），提出了一種適用于非飽和膨脹土的熱力學模型。

Graham等[提出了一種考慮溫度影響的修正劍橋模型，該模型能夠預測正常固結和超固結飽和黏王在排水加熱后產生更高的不排水剪切強度和更低的孔隙水壓力，而經(jīng)歷不排水加熱過程的飽和黏土則產生更高的孔隙水壓力和更低的不排水剪切強度，這主要是由于屈服點區(qū)域和超孔隙比的減小所致。Hamidi等僅用一個屈服面，在修正劍橋模型基礎上，提出了一種考慮應力水平和溫度影響的熱力學本構模型。在該模型基礎上，利用加熱后土體的壓縮曲線和常溫時壓縮曲線與臨界狀態(tài)線之間的關系，結合臨界狀態(tài)理論，Hamidi等2建立了一種考慮應力歷史的熱力學本構模型。Wang等[73在修正劍橋模型的基礎上引入兩個參數(shù) α 和 θ 分別用于考慮應力引起的土體各向異性和溫度對壓縮曲線的影響，同時，采用兩個屈服面，即當前屈服面和參考屈服面來描述土體在力學固結和熱固結過程中產生的塑性變形，最終提出了一種適用于飽和黏土的各向異性熱彈塑性本構模型。

2.3 邊界面模型

前述熱力學本構模型常用于模擬土體在大應變情況下的應力變形特性，但這些模型通常假定屈服面形狀與溫度無關，僅尺寸隨溫度變化，導致預測土體剪切特性時產生不可忽略的誤差。因此，Zhou等[74]運用邊界面理論[75-77]，考慮屈服面的形狀和尺寸皆隨溫度變化而發(fā)生相應變化，提出了一種可考慮小應變和大應變的飽和土體熱彈塑性模型。該模型可描述升溫和降溫過程中土體產生的體積變形以及正常固結和超固結土體在排水和不排水條件下的剪切特性，特別是小應變（小于 0.001% ）下土體剪切模量的退化過程。

在該模型基礎上，在邊界面內部區(qū)域考慮塑性應變，Zhou等78推導出溫度-水-循環(huán)荷載耦合作用下非飽和土體的應力變形關系。結合Bishop應力、偏應力、基質吸力、溫度、比熱容和飽和度的變化對非飽和土體的影響，提出了一種可模擬非飽和土體在溫度變化、循環(huán)荷載作用下力學特性的邊界面模型。

此外，Zhou等79采用邊界面模擬先期固結應力，利用記憶面模擬最大主應力隨溫度循環(huán)變化的演化過程，通過溫控試驗發(fā)現(xiàn)，邊界面與記憶面之間的距離和塑性模量都隨著溫度循環(huán)過程中塑性應變的累積而增加。隨后，在雙屈服面理論和單屈服面邊界面熱力學模型的基礎上，提出了一種可考慮多重溫度循環(huán)的飽和土體邊界面模型。Laloui等[23通過分析已有文獻的溫控試驗數(shù)據(jù)，得出溫度變化對土體的體積變形、先期固結應力、正常固結線、剪切強度、孔隙水壓力等特性的影響規(guī)律。在臨界狀態(tài)理論框架下，運用多機制塑性理論和邊界面理論，提出了一種適用于復雜應力路徑、排水和不排水剪切條件的熱力學模型。陳艷妮等[80]基于熱力學理論，建立屈服面，采用非關聯(lián)流動法則，同時，在邊界面模型框架下，提出了一種適用于超固結黏土的考慮溫度效應的邊界面模型。

2.4 UH模型

在劍橋模型的基礎上，結合廣義Mises、SMP、Lade等強度準則，姚仰平等[81-83]提出考慮漸進狀態(tài)、顆粒破碎、各向異性條件下的系列UH模型。

通過已有的溫控試驗數(shù)據(jù)，姚仰平等[84總結土體的壓縮特性、先期固結壓力特性以及剪切特性等，采用超固結土潛在強度的確定方法[85]，建立某一溫度作用下土體的臨界狀態(tài)應力比表達式

其中：

式中： M_T 和 M₀ 分別為當前溫度 T 和初始溫度 T_?0 作用下的臨界狀態(tài)應力比； ξ^_T，γ 和 χ₀ 均為材料參數(shù)。最后，基于SMP準則，將模型進行三維化，提出適用于超固結土體的考慮溫度影響的本構模型。該模型可描述飽和超固結土體的熱變形特性，以及溫度不變時，正常固結土和超固結土體在復雜應力路徑下的應力變形特性。在該模型的基礎上，姚仰平等8將不同溫度作用下的屈服線與溫度不變時不同基質吸力作用下的屈服線進行結合，建立了可同時考慮溫度和基質吸力作用的屈服面，在臨界土力學框架內，提出了適用于正常固結和超固結非飽和土體的熱力學模型。

在熱UH模型84的基礎上，姚仰平等8推導了考慮溫度影響的不排水剪切強度表達式，并提出了一種可考慮熱破壞現(xiàn)象的本構模型，即在不排水條件下，常溫剪切至某狀態(tài)時，應力比保持不變，增加溫度至土體達到臨界狀態(tài)應力比。此外，孔立明等[88]在考慮時間效應的三維彈-黏-塑性本構模型[89]的基礎上，建立溫度與模型參數(shù)的關系式，從而推導考慮溫度和時間效應的應力-應變關系式，并將屈服面與UH模型[9中的屈服方程和流動法則相結合，提出了一種適用于超固結土體的考慮溫度和時間效應的熱-黏-彈-塑性本構模型。

2.5 熱力學模型

熱力學是從宏觀角度研究物質的熱運動性質及其規(guī)律的學科，主要從能量轉化的角度研究物質的熱性質，揭示能量從一種形式轉換為另一種形式時所遵循的宏觀規(guī)律。

Zhang等[91通過定義彈性勢能密度函數(shù)，推導有效應力、彈性應變和干密度之間的超彈性關系式，運用熱力學理論[92-96]，定量描述材料非彈性變形發(fā)展等不可恢復的能量過程；采用顆粒熵的概念，描述由于顆粒不規(guī)則運動引起的顆粒尺度上的動能漲落和彈性勢能漲落，提出了一種能夠模擬溫度循環(huán)作用下飽和黏土的非等溫固結過程的熱力學模型。張志超等[9引入耗散力、耗散流和顆粒熵的概念，推導出一種非彈性變形的多物理場力學模型，該模型可以描述溫度升高時飽和土體內部結合水向自由水轉化而發(fā)生的能量耗散過程（即飽和土體的非等溫固結過程），還可以描述經(jīng)歷不同應力歷史后，飽和土體的不排水剪切過程。在該模型基礎上，Zhang等[98引入適用于單相固體材料的流體動力學理論（GSH），構建了能夠描述黏土在溫度和應力耦合作用（T-M）下應力變形機理的理論模型（TTS）。

Bai等[99]基于多孔顆粒材料在熱-力耦合作用下發(fā)生粒子重排的現(xiàn)象，引入粒子熵、粒子溫度、遷移系數(shù)和勢能密度函數(shù)等概念，建立了適用于巖土材料的熱力學模型。該模型能考慮不同溫度和飽和度條件下巖土材料的能量耗散過程和土-水特征曲線的變化規(guī)律，能較準確地描述正常固結飽和土體在升溫過程中的不可逆固結現(xiàn)象以及多重加熱-降溫過程中土體固體顆粒發(fā)生不可逆重組所導致的巖土材料老化現(xiàn)象；此外，該模型還能較準確地描述圍壓、基質吸力和溫度對非飽和土固結過程的影響，其物理機制與飽和土體熱固結過程的物理機制相同。Yang等[o]將該模型擴展到一種特殊的非飽和土體（細粒氣態(tài)土），該種土體中含有尺寸較正常顆粒尺寸大很多的密封狀態(tài)的氣泡，結合考慮大氣泡的理論模型[101]，建立了描述細粒氣態(tài)土力學特性和溫度效應的熱力學模型。該模型能夠考慮氣體對土體骨架塑性變形的影響以及不同排水條件和超固結比作用下溫度對細粒氣態(tài)土固結和不排水剪切特性的影響。Yang等[102]基于顆粒物質的熱力學機理，考慮溫度對細觀能量耗散的影響，建立了描述飽和砂土熱力學特性的模型。Golchin等[103在能量勢中引入熱彈性，將熱力學原理嵌入到新提出的耗散勢中，提出了一種能夠從微觀尺度定性描述細粒土熱力學性質的模型。

2.6 其他模型

Kurz等[104]在常溫彈-黏-塑性本構模型的基礎上，將飽和黏土的應變分為彈性和塑性兩個部分，且塑性應變部分包含黏性應變，從而提出了一種適用于飽和黏土的半經(jīng)驗彈-熱-黏-塑性本構模型。在該模型中，黏性應變由蠕變率系數(shù)定義，其值與溫度和塑性指數(shù)有關，與圍壓和超固結比無關，Laloui等[105通過先期固結壓力在加熱過程中的變化特征，研究屈服極限隨溫度的演化規(guī)律，從而將等溫條件下的多機構彈塑性模型推廣到非等溫條件，用于考慮黏土因加熱而呈現(xiàn)的非線性和硬化特性。該模型能夠描述黏土在不同溫度和應力作用下的完整行為，即等溫或非等溫、排水或不排水剪切、不同溫度和應力歷史條件下土體的力學行為。在該模型的基礎上，Laloui等[106]將其延伸到循環(huán)荷載條件，提出了一種循環(huán)荷載作用下飽和土體的非等溫塑性模型。在雙硬化模型[108-110]和溫控試驗數(shù)據(jù)的基礎上，Liu等1]引入兩個硬化參數(shù)，分別考慮溫度和超固結比對土體應力變形的影響，從而提出了適用于超固結飽和黏土的熱彈塑性本構模型。

在溫控循環(huán)試驗數(shù)據(jù)的基礎上，Ma等[1]運用安定理論[112-114]，提出一種熱穩(wěn)定線（TSL）。該線可表示土體在溫度循環(huán)作用下所處的穩(wěn)定狀態(tài)，且與溫度循環(huán)階段有關，當土體經(jīng)歷加熱階段時，若土體位于TSL上方，表明土體產生了不可逆的體積壓縮變形（熱塑性變形），若土體位于TSL下方，則產生了熱彈性變形；當土體經(jīng)歷降溫階段時，將產生熱彈性變形;利用亞塑性理論[115]，提出了一種適用于正常固結或輕微超固結飽和土體的本構模型，該模型可以描述飽和土體經(jīng)歷一重或多重溫控循環(huán)后的累積不可逆體積變形。Masin[115運用修正劍橋模型確定正常固結線（NCL）和臨界狀態(tài)線（CSL），運用廣義亞塑性理論表征屈服面和不同超固結比的土體所呈現(xiàn)的非線性，運用晶間應變概念表述土體小應變階段的特性，提出了一種可考慮應力歷史和小應變的亞塑性本構模型。在此基礎上，Masin等[116提出了一種可描述非飽和土體非等溫加熱條件下的熱力學亞塑性本構模型，該模型可用于正常固結和超固結土體。

Zhang等[117將加熱過程中溫度變化使土體產生的熱彈性變形部分等效為土體發(fā)生力學卸載過程。兩種不同作用過程（加熱和力學卸載）中產生的彈性變形保持相等，并將該理論上的力學卸載應力稱為等效應力，用以表示溫度對土體彈性變形的影響，如圖3所示。運用上加載屈服面理論，推導出等效的下加載屈服面，用以表示溫度和應力對土體塑性變形的影響。在修正劍橋模型的基礎上，僅引入一個參數(shù)，提出了一種適用于超固結土體的熱彈塑性本構模型。在此模型基礎上，張升等[118引入一個參量，即線性熱膨脹系數(shù)，提出了一個考慮結構性的軟巖熱彈塑性本構模型。該模型可考慮軟巖熱增強和熱減弱兩種現(xiàn)象，且初始結構性的不同可使軟巖在熱減弱和熱加強之間相互轉化。

在實際能源結構工程中，土體長期經(jīng)歷應力和溫度的影響，時間對土體力學特性的影響顯得格外重要。因此，非常有必要研究溫度和時間作用下軟巖等土體的蠕變破壞特性。Zhang等[119開展一系列軟巖的熱蠕變試驗，并根據(jù)試驗結果，運用已有的熱力學模型[117]，提出了一種適用于超固結飽和土體的考慮溫度和時間效應的本構模型。Coccia等22通過溫控試驗發(fā)現(xiàn)，加熱正常固結土體將使土體加熱前的二次壓縮過程加速完成，且僅發(fā)生壓縮現(xiàn)象，但加熱超固結土體時，將發(fā)生熱膨脹或熱壓縮現(xiàn)象；基于溫控試驗數(shù)據(jù)，提出了一種可描述飽和、非飽和土體在不同應力狀態(tài)和加熱速率條件下熱體積變形的本構模型。Laloui等[120]提出了一種可描述一維固結溫控試驗中考慮溫度和應變率變化、適用于飽和黏土的熱-黏-塑性本構模型。

2.7 模型的應用

通過溫控試驗總結土體的熱力學性質，結合相關理論建立土體的本構模型，其最終目的是將本構模型通過二次開發(fā)引入有限元軟件中，并通過數(shù)值分析方法開展工程應用。

祁良等51通過溫控試驗得出，軟黏土的臨界狀態(tài)應力比 M 隨溫度的增加而線性增加，同時，根據(jù)試驗數(shù)據(jù)進行擬合，提出 M 與溫度 T 的線性方程。利用ABAQUS有限元軟件，采用修正劍橋模型[121]對三軸試樣進行有限元模擬，固結和剪切條件與實際試驗一致，同時在固結過程添加溫度條件。通過模擬值和試驗值的對比分析發(fā)現(xiàn)，常溫（ 25^°C） 下試驗值與模擬值吻合較好，但當溫度為 45，65^°C 時，試驗值與模擬值有較大誤差，該誤差是由于模擬過程中采用的修正劍橋模型并未考慮溫度影響所致。此研究也間接說明，在溫度相關的實際工程中，提出合適的考慮溫度影響的本構模型非常必要。

費康等2在有限元軟件COMSOL5.0內置的修正劍橋模型的基礎上，將屈服函數(shù)更改為Zhou模型[4]中的屈服面，并通過修改其他材料參數(shù)，對Zhou熱邊界面模型[4進行二次開發(fā)；通過對一維固結問題、常溫等向壓縮問題、排水條件下的溫度循環(huán)問題等進行模擬，驗證了二次開發(fā)的可靠性。陸浩杰等[122將ACMET-T本構模型[23]在ABAQUS中進行二次開發(fā)，并對溫控三軸試驗結果進行模擬分析，將模擬結果與試驗結果進行對比，驗證該數(shù)值分析方法的可行性和準確性。馮興等[123將熱UH模型[84]中的非對稱剛度矩陣修正為對稱剛度矩陣，采用Newton-Raphson算法和半隱式回映的應力更新算法，編制熱UH模型的有限元程序；利用有限元非線性求解平臺，實現(xiàn)熱UH模型在有限元中的應用，并與溫控試驗數(shù)據(jù)進行對比分析，驗證了程序的可靠性。王興等[25引入一個附加硬化參數(shù)，在邊界面塑性理論框架下，提出考慮水合物和溫度對土體應力變形影響的邊界面模型；通過ABAQUS的用戶材料子程序接口，開發(fā)該模型的UMAT子程序，并模擬了不同水合物條件下能源土的應力變形特性。費康等2通過ABAQUS子程序對Zhou模型4進行二次開發(fā)，材料的熱膨脹系數(shù)利用ABAQUS內置的材料性質定義，其余力學參數(shù)則通過UMAT接口定義，應力積分過程采用Euler算法，研究了溫度對黏土壓縮特性和剪切特性的影響；隨后，建立能源樁和地基的數(shù)值模型，其中樁-土接觸面采用薄層單元，材料模型及參數(shù)與周圍地基土一致，地基土的應力-應變關系采用二次開發(fā)的熱邊界面模型，并利用該數(shù)值模型，計算分析了樁頂及樁身不同位置的溫度響應、樁頂位移、樁身側摩阻力與樁身軸力等。王浩等[124在TTS模型[98框架下建立了可描述水泥石在熱-力耦合作用下力學特性的黏-彈-塑性模型，并基于MATLAB平臺，開發(fā)瞬態(tài)傳熱平面有限元程序，對地熱井的 1/4 橫斷面進行數(shù)值模擬，研究了升溫過程中套管與水泥環(huán)接觸面以及水泥環(huán)與圍巖界面的應力變形發(fā)展規(guī)律。

3 研究建議

3.1 熱力學試驗研究

1）現(xiàn)有的溫控試驗大多局限于黏土、粉土、砂土等單一土體，是一種理想狀態(tài)下的土體，缺乏對粉砂、黏土-膨潤土等混合土體熱力學行為的研究。實際工程中的土體并不是純砂土、純粉土或純黏土等單一土體，而是一種具有不同細粒含量的砂土、粉土或黏土的混合土體，且粉土或黏土占整個土體的比例不同時，土體呈現(xiàn)的熱力學特性也有所不同。在今后的溫控試驗研究中，建議針對不同混合體或實際工程中的土體開展熱力學試驗，如此，室內的熱力學試驗結果更接近于實際工程中的監(jiān)測值，更具有參考價值。

2）目前已有的溫控試驗主要集中于重塑土，該試樣具有分層沉積特性，可考慮初始各向異性，但對于原狀土的溫控試驗研究較少，即考慮土體的次生各向異性的研究較少。在今后的溫控試驗研究方面，建議多針對原狀土，多考慮土體的結構性進行。

3）已有的溫控靜三軸、動三軸或空心扭剪試驗均考慮土體在溫度作用下的力學性質，僅考慮由溫度引起的土體本身體積或強度變化，未考慮與溫度有關的結構物以及結構物與周圍土體接觸面的熱力學性質。而能源結構物與土體接觸面往往是整個結構的薄弱環(huán)節(jié)。在今后的試驗研究方面，建議多開展能源結構物和土體的模型試驗以及現(xiàn)場的熱響應試驗，這樣既能考慮溫度變化對能源結構物和土體本身的力學性質的影響，又能考慮兩者之間接觸面的影響，從而更全面地考慮能源結構物和土體之間的相互作用。其試驗結果更貼合實際工程，可以更好地應用于實際能源工程的設計與施工中。

4）以溫控靜三軸試驗為例，已有的熱力學試驗過程常常為：土體先在常溫下進行力學固結，隨后，壓力不變，進行加熱或溫度循環(huán)等熱固結過程，然后，溫度不變，在排水或不排水條件下進行剪切過程。在該過程中，溫度和力學是分開加載的，不是真正意義上的熱力耦合作用，而在實際能源工程中，土體中的應力場和溫度場是同時存在的。為了減少室內試驗和實際工程的誤差，開展熱破壞、應力-溫度耦合加載、非比例加載、多重溫控循環(huán)等復雜載荷作用下土體熱力學行為試驗研究尤為重要，同時為相關工程應用提供堅實的理論支持。

5）已有的土體熱力學性質，主要是通過溫控試驗研究土體在溫度變化下的變形和強度等宏觀響應，是土體整體上的力學特征，但對于土體內部，特別是溫度變化下試樣不同部位孔隙水和土體顆粒之間的力學變化特性的研究相對較少。在之后的熱力學特性研究方面，可利用CT、PFC等手段，研究在溫度變化過程中土體內部發(fā)生的熱力學特性，即從微觀維度探索土體的熱力學機理。

3.2 熱力學模型研究

1）已有的熱力學模型主要在劍橋或修正劍橋模型基礎上建立先期固結應力、二次壓縮系數(shù)、超固結比等參數(shù)與溫度之間的關系，再結合臨界狀態(tài)理論，提出考慮時間效應、應力歷史、各向異性等適用于飽和與非飽和粉土（黏土）的熱力學模型。雖然考慮的作用形式較多，但土體種類較為單一，主要集中于黏土，對砂土或實際工程中較為常見的砂王-黏土混合物的熱力學模型的研究相對較少。在今后的熱力學模型方面，建議針對含有一定量粉土或黏土的粉砂類混合土體的力學模型進行研究，提出飽和粉砂在應力場及溫度場耦合作用下的狀態(tài)相關力學模型具有重要的理論意義。

2）已有的熱力學模型建模過程大多較為復雜，參數(shù)不易確定，只適用于某些特定加載條件下土體的應力-應變響應。此外，大多數(shù)模型是基于一維固結或三軸壓縮應力狀態(tài)建立的，關于三維空間中土體強度和變形特性的研究較少。在今后的熱力學模型研究中，建立易于確定參數(shù)的、考慮三維空間的熱力學模型，從而更好地預測土體復雜的熱力學特性是非常必要的。

3）目前對土體熱力學本構關系的研究仍處于探索階段。在溫控試驗的基礎上，不少研究者已經(jīng)建立了一些土體熱力學本構關系，但如何在現(xiàn)有的理論基礎上，面對復雜載荷條件下（如應力場和溫度場耦合加載，非比例加載和升溫、多重溫控循環(huán)、循環(huán)荷載等）土體的力學響應，發(fā)展出更為合理準確的本構關系，還有待進一步深入研究。

3.3 模型應用研究

已有的數(shù)值應用，通過不同的數(shù)值積分算法，將熱力學模型應用于ABAQUS、COMSOL、MATLAB等軟件中，建立考慮溫度效應的有限元模型，并對能源結構進行數(shù)值模擬。由于對數(shù)值積分算法的認識還不夠統(tǒng)一，熱力學模型在實際工程中的應用遠遠滯后于本構理論的發(fā)展。對于不同的有限元軟件，采用的數(shù)值積分算法各不相同。此外，已有的數(shù)值應用常模擬溫控三軸試驗過程和結果，數(shù)值模擬過于簡單，沒有能源結構物的參與，與實際工程差異較大。

在今后的數(shù)值應用中，針對不同框架下建立的熱力學模型，建立可靠且通用的數(shù)字積分算法，并將該算法開發(fā)到有限元中，利用有限元對縮尺的模型試驗或足尺的現(xiàn)場熱響應試驗進行模擬，將模擬值與監(jiān)測值進行對比分析，從而驗證該數(shù)值方法的可靠性，使模擬值更具有普適性，更好地應用于實際能源工程的設計與施工中。

4 結論與展望

淺層地熱能是一種典型的可持續(xù)能源，其開發(fā)技術具有形式多樣、高效節(jié)能、綠色低碳等特點。中國淺層地熱資源豐富、類型齊全、分布廣泛，能夠形成或催生多領域、更廣泛的利用方式，具有廣闊的工程應用前景。地熱資源開發(fā)過程中涉及應力場和溫度場的耦合作用，土體、結構物以及兩者之間接觸面的力學行為的外在表現(xiàn)與內在機理均非常復雜。

微觀維度上的熱力學特性研究尚不多見，宏觀維度上的熱力學試驗研究多為力學荷載和溫度荷載分開加載，是“假\"熱-力耦合作用，而對于力學荷載和溫度荷載同時加載的“真\"熱-力耦合作用的試驗研究較少。因此，可在“假\"熱-力耦合溫控試驗基礎上，通過分析力學固結和熱固結兩個階段土體所發(fā)生的應力變形特性，分析該特性產生的力學機理，從而分析更為復雜的“真\"熱-力耦合溫控試驗中土體的熱力學特性及其機理。室內溫控試驗研究還可以從常規(guī)的應力路徑（力學固結-熱固結-剪切）過渡到主應力旋轉等復雜載荷下土體的熱力學響應。室內溫控試驗主要研究土體的熱力學特性，未考慮結構物的影響，其結果僅適用于能源結構工程的地基變形，即僅考慮能源結構物周圍土體的熱應力變形，而無法考慮能源結構物本身及其與周圍土體接觸面的熱力學特性。因此，通過縮尺模型槽試驗（如能量樁模型試驗）和足尺現(xiàn)場熱響應試驗可以監(jiān)測能源結構物、周圍土體及其接觸面的應力變形特性，結果更接近于實際能源工程。此外，將微觀尺度上的溫控試驗數(shù)值模擬，宏觀尺度上的室內常規(guī)溫控試驗，模型槽試驗與現(xiàn)場熱響應試驗結合起來，通過開展不同尺度的熱力學試驗，全方位考慮尺寸效應對土體及結構物熱力學性質的影響，從而基于相關理論，提出考慮土體熱力學微觀機理與宏觀性質的本構關系，從而加深對土體熱塑性變形機理的認識，這也是今后土體熱力學性質研究的重中之重。

對于土體熱力學本構模型的研究，其目的是進行工程應用，因此，加強土體熱力學本構關系的二次開發(fā)，充分利用有限元計算分析軟件開展相關工程數(shù)值模擬研究，建立考慮溫度效應的土體數(shù)值分析方法非常必要。

總之，盡管目前已經(jīng)取得了一定的研究成果，這些成果也在淺層地熱能開發(fā)過程中發(fā)揮著非常重要的理論指導作用，但關于土體熱力學本構關系的研究還處于起步階段，在熱力學試驗、本構理論和數(shù)值應用方面仍存在諸多問題，亟須克服。這也需要更多的巖土工作者加倍努力，積極創(chuàng)新，不畏困難，從不同角度開展研究，從而建立土體熱力學性質的研究體系，進一步認清土體的熱力學機理，將其更好地應用于能源工程中。

參考文獻

[1」李政，張東杰，潘玲穎，等.\"雙碳\"目標卜我國能源低 碳轉型路徑及建議[J].動力工程學報，2021，41（11）： 905-909，971. LI Z， ZHANG D J， PAN L Y，et al. Low-carbon transition of China’s energy sector and suggestions with the “Carbon-peak and carbon-neutrality” target [J]. Journal of Chinese Society of Power Engineering，2O21，41（11）： 905-909，971.（in Chinese）

［2］馬冰，賈凌霄，于洋，等.世界地熱能開發(fā)利用現(xiàn)狀與 展望[J].中國地質，2021，48（6）：1734-1747. MA B，JIAL X，YU Y，et al.The development and utilization of geothermal energy in the world[J]. Geology in China，2021，48（6）： 1734-1747. （in Chinese）

[3］馬峰，王貴玲，魏帥超，等.2018年地熱勘探開發(fā)熱點 回眸[J].科技導報，2019，37（1）：134-143. MAF，WANGGL，WEISC，etal.Summary of hot research topics in geothermal exploitation in 2018 [J]. Science amp; Technology Review，2019，37（1）：134-143. （in Chinese）

[4］龐忠和，胡圣標，汪集旸.中國地熱能發(fā)展路線圖[J]. 科技導報，2012，30（32）：18-24. PANG Z H，HU SB，WANG JY.A roadmap to geothermal energy development in China [J]. Science amp;. Technology Review，2012，30（32）： 18-24.（in Chinese）

[5]GENS A，SANCHEZ M，GUIMARAES L D，et al. A full-scale in situ heating test for high-level nuclear waste disposal： Observations，analysis and interpretation [J].Géotechnique，2009，59（4）： 377-399.

[6]BOURNE-WEBB P J，AMATYA B， SOGA K， et al. Energy pile test at Lambeth College，London： Geotechnical and thermodynamic aspects of pile response to heat cycles[J].Geotechnique，2009，59（3）： 237-248.

[7］FADEJEV J，SIMSON R，KURNITSKI J， et al.A review on energy piles design，sizing and modelling [J]. Energy，2017，122： 390-407.

[8]KNELLWOLF C，PERON H，LALOUIL.Geotechnical analysis of heat exchanger piles [J]. Journal of Geotechnical and Geoenvironmental Engineering，2O11，137 （10）： 890-902.

[9]CROLL JG A. The role of thermal ratcheting in pavement failures [J].Proceedings of the Institution of Civil Engineers- Transport，2009，162（3）： 127-140.

[10] KERTESZ R， SANSALONE J. Hydrologic transport of thermal energy from pavement [J]. Journal of Environ mental Engineering，2014，140（8）： 04014028.

[11]POTHIRAKSANON C，BERGADO D T，ABUELNAGA H M. Full-scale embankment consolidation test using prefabricated vertical thermal drains [J]. Soils and Foundations，2010，50（5）： 599-608.

[12] BASER T，LU N，MCCARTNEY J S. Operational re sponse of a soil-borehole thermal energy storage system [J]. Journal of Geotechnical and Geoenvironmental Engi neering，2016，142（4）： 04015097.

[13] BOURNE-WEBB P，BURLON S，JAVED S，et al. Analysis and design methods for energy geostructures [J]. Renewable and Sustainable Energy Reviews，2016， 65： 402-419.

[14] BRANDL H. Energy foundations and other thermoactive ground structures [J]. Géotechnique，2OO6，56（2）： 81-122.

[15] LALOUI L，DI DONNA A. Understanding the behav iour of energy geo-structures[J]. Proceedings of the Insti tution of Civil Engineers - Civil Engineering，2Ol1，164 （4）： 184-191.

[16]ABUEL-NAGA H M，BERGADO D T，RAMANA G V，et al. Experimental evaluation of engineering behavior of soft bangkok clay under elevated temperature [J]. Journal of Geotechnical and Geoenvironmental Engineering，2006，132（7）： 902-910.

[17] CEKEREVAC C，LALOUI L. Experimental study of thermal efects on the mechanical behaviour of a clay [J]. International Journal for Numerical and Analytical Methods in Geomechanics，2004，28（3）： 209-228.

[18] KINOSHITA N，YASUHARA H. Thermally induced behavior of the openings in rock mass affected by high temperatures [J]. International Journal of Geomechanics， 2011，11（2）： 124-130.

[19]NG C W W， ZHOU C.Cyclic behaviour of an unsaturated silt at various suctions and temperatures [J]. Geotechnique， 2014， 64（9）： 709-720.

[20] NG C W W，WANG S H， ZHOU C. Volume change behaviour of saturated sand under thermal cycles [J]. Géotechnique Letters， 2016， 6（2）： 124-131.

[21]HUECKEL T，F(xiàn)RANCOIS B，LALOUI L. Explainingthermal failure in saturated clays[J].Geotechnique， 2009，59（3）： 197-212.

[22] COCCIA C JR，MCCARTNEY JS. Thermal volume change of poorly draining soils II ： Model development and experimental validation [J]. Computers and Geotechnics，2016，80： 16-25.

[23] LALOUI L，F(xiàn)RANCOIS B. ACMEG-T： Soil thermo plasticity model [J]. Journal of Engineering Mechanics， 2009，135（9）： 932-944.

[24] YAO Y P， ZHOU A N. Non-isothermal unified hardening model： A thermo-elasto-plastic model for clays [J]. Géotechnique，2013， 63（15）： 1328-1345.

[25］王興，孔亮，袁慶盟，等.能源土邊界面模型的應力更 新算法及算例分析[J].土木工程學報，2019，52（Sup2）： 155-161. WANG X，KONG L，YUAN QM，et al. The stress update algorithm of bounding surface model on gas hydrate-bearing soil and case study [J]. China Civil Engineering Journal，2019，52（Sup2）：155-161.（in Chinese）

[26]費康，錢健，洪偉，等.黏土地基中能量樁力學特性數(shù) 值分析[J].巖土力學，2018，39（7）：2651-2661. FEI K，QIANJ，HONG W，et al. Numerical analysis of mechanical behavior of energy piles in clay [J].Rock and Soil Mechanics，2018，39（7）： 2651-2661.（in Chinese）

[27]費康，劉漢龍，孔綱強，等.熱力耦合邊界面模型在 COMSOL中的開發(fā)應用[J].巖土力學，2017，38（6）： 1819-1826. FEI K，LIU HL，KONG G Q，et al. Implementation of a thermo-bounding surface model in COMSOL [J]. Rock and Soil Mechanics，2017，38（6）：1819-1826. （in Chinese）

[28]劉漢龍，王成龍，孔綱強，等.U型、W型和螺旋型埋管 形式能量樁熱力學特性對比模型試驗[J].巖土力學， 2016，37（Supl）： 441-447. LIU H L，WANG C L， KONG G Q， et al. Comparative model test on thermomechanical characteristics of energy pile with U-shape，W-shape and spiral-shape [J]. Rock and Soil Mechanics，2016，37（Supl）： 441-447. （in Chinese）

[29]黨政，關文，程曉輝，等.CFG能源樁用于混凝土路面 除冰降溫的試驗研究[J].中國公路學報，2019，32（2）： 19-30. DANG Z，GUAN W，CHENG X H，et al. Experimental study on CFG energy pile for concrete pavement deicing and cooling [J]. China Journal of Highway and Transport，2019，32（2）：19-30.（in Chiese）

[30]吳道勇，賴遠明，馬勤國，等.季節(jié)凍土區(qū)水鹽遷移及 土體變形特性模型試驗研究[J].巖土力學，2016，37 （2）： 465-476. WU DY，LAI Y M，MAQG，et al. Model test study of water and salt migration and deformation characteris tics in seasonally frozen soil [J].Rock and Soil Mechan ics，2016，37（2）： 465-476.（in Chinese）

[31]黃旭，孔綱強，劉漢龍，等.循環(huán)溫度場作用下PCC能 量樁熱力學特性模型試驗研究[J].巖土力學，2015，36 （3）： 667-673. HUANG X，KONG G Q，LIU HL，et al. Experimental research on thermomechanical characteristics of PCC energy pile under cyclic temperature field [J]. Rock and Soil Mechanics，2015，36（3）： 667-673. （in Chinese）

[32]桂樹強，程曉輝.能源樁換熱過程中結構響應原位試驗 研究[J].巖土工程學報，2014，36（6）：1087-1094. GUIS Q，CHENG X H.In-situ tests on structural responses of energy piles during heat exchanging process [J]. Chinese Journal of Geotechnical Engineering，2014， 36（6）： 1087-1094. （in Chinese）

[33］陳家威，張國柱，郭易木，等.層狀地層能源管樁傳熱 性能試驗研究[J].巖石力學與工程學報，2020，39 （Sup2）： 3615-3626. CHEN JW， ZHANG G Z， GUO Y M， et al. Investigation on heat transfer characteristics of PHC energy piles in multi-layer strata [J]. Chinese Journal of Rock Mechan ics and Engineering，2020，39（Sup2）：3615-3626.（in Chinese）

[34]CASAGRANDEB，SABOYAFJr，MCCARTNEYJ S，et al. Investigation of a field-scale energy micropile in stratified soil under cyclic temperature changes [J]. Geomechanics for Energy and the Environment，2022， 29： 100263.

[35］姚仰平，張丙印，朱俊高.土的基本特性、本構關系及 數(shù)值模擬研究綜述[J].土木工程學報，2012，45（3）： 127-150. YAO Y P， ZHANG B Y， ZHU JG.Behaviors， constitutive models and numerical simulation of soils [J]. China Civil Engineering Journal，2Ol2，45（3）： 127-150. （in Chinese）

[36]李春紅，孔綱強，車平，等.能量樁樁-土接觸面力學特 性室內試驗研究[J].建筑節(jié)能，2016，44（3）：99- 105，114. LI C H，KONG G Q，CHE P，et al.Laboratory experimental on interface mechanical properties of energy pilesoil [J].Building Energy Efficiency，2Ol6，44（3）：99- 105，114.（in Chinese）

[37]李春紅，孔綱強，劉漢龍，等.樁-紅黏土接觸面溫控測 試及應力-應變關系研究[J].土木工程學報，2019，52 （Sup 2）： 89-94，101. LI C H，KONG G Q，LIU HL，et al.Study of temperature-controlled pile-red clay interface tests and stress-strain relationship [J]. China Civil Engineering Journal，2019，52（Sup 2）：89-94，101.（in Chinese）

[38]李春紅，孔綱強，張鑫蕊，等.溫控樁-土接觸面三軸試 驗系統(tǒng)研制與驗證[J].巖土力學，2019，40（12）：4955- 4962. LIC H，KONG G Q，ZHANG X R，et al. Development and verification of temperature-controlled pile-soil interface triaxial shear test system [J]. Rock and Soil Mechanics，2019，40（12）： 4955-4962. （in Chinese）

[39]BALDIG，HUECKEL T，PELLEGRINI R. Thermal volume changes of the mineral-water system in lowporosity clay soils [J]. Canadian Geotechnical Journal， 1988， 25（4）： 807-825.

[40] BAI B，GUO L J，HAN S. Pore pressure and consolidation of saturated silty clay induced by progressively heating/cooling [J]. Mechanics of Materials，2O14，75： 84-94.

[41]DE BRUYN D，THIMUS J F. The influence of temperature on mechanical characteristics of Boom clay： The results of aninitial laboratoryprogramme [J]. Engineering Geology，1996，41（1-4）： 117-126.

[42]BURGHIGNOLI A，DESIDERIA，MILIZIANO S. Alaboratory study on the thermomechanical behaviour of clayey soils [J]. Canadian Geotechnical Journal，2000， 37（4）： 764-780.

[43]FAVERO V，F(xiàn)ERRARI A，LALOUI L. Thermomechanical volume change behaviour of Opalinus Clay [J]. International Journal of Rock Mechanics and Mining Sciences，2016，90：15-25.

[44] ROMERO E，VILLAR M V，LLORET A. Thermohydro-mechanical behaviour of two heavily overconsolidated clays [J]. Engineering Geology，2OO5，81（3）： 255-268.

[45]潘旋.凍結鈣質黏土應力-應變特性的真三軸試驗研究 [D].安徽 淮南：安徽理工大學，2020. PAN X. True tri-axial test on stress-strain characteristics of frozen calcareous clay [D]. Huainan，Anhui： Anhui Universityof Science amp;. Technology，2020.（in Chinese）

[46] ZHOU C，NG C W W. Effects of temperature and suction on plastic deformation of unsaturated silt under cy clic loads [J]. Journal of Materials in Civil Engineering， 2016，28（12）： 04016170.

[47]LIU H，LIU HL，XIAO Y，et al.Effectsof temperature on the shear strength of saturated sand [J]. Soils and Foundations，2018，58（6）： 1326-1338.

[48] LIU H，LIU H L，XIAO Y，et al. Influence of temperature on the volume change behavior of saturated sand [J]. Geotechnical Testing Journal，2Ol8，41（4）： 20160308.

[49]LIU H，MCCARTNEY J S，XIAO Y. Thermal volume changes of saturated sand during loading-unloadingheating phase [J]. E3S Web of Conferences，2O20，205： 08002.

[50]郭楨，劉干斌，尹鐵鋒，等.考慮溫度的軟黏土鄧肯-張 模型參數(shù)試驗研究[J].建筑結構，2014，44（3）：93- 96，83. GUO Z，LIU G B，YIN T F，et al. Test study on parameters of Duncan-Chang model for soft clay with considerationof temperature influence [J]. Building Structure，2014， 44（3）： 93-96，83.（in Chinese）

[51］祁良，鄭榮躍，陶海冰，等.考慮溫度影響的寧波軟土 臨界狀態(tài)參數(shù)研究[J].水文地質工程地質，2015，42 （5）： 79-83. QI L，ZHENG RY，TAO H B，et al. A study of the critical state parameter of the Ningbo soft clay considering the effect of temperature [J]. Hydrogeology amp; Engineering Geology，2015，42（5）： 79-83.（in Chinese）

[52]謝云，陳正漢，李剛.考慮溫度影響的重塑非飽和膨脹 土非線性本構模型[J].巖土力學，2007，28（9）：1937- 1942. XIE Y，CHEN Z H，LIG. Thermo-nonlinear model for unsaturated expansive soils [J]. Rock and Soil Mechanics，2007，28（9）： 1937-1942.（in Chinese）

[53] LIU H，LIU HL，XIAO Y，et al. Non-linear elastic model incorporating temperature effects [J]. Geotechnical Research，2018，5（1）： 22-30.

[54] XIAO Y，LIU H，ZHANG W G，et al. Testing and modeling of rockfill materials：A review [J]. Journal of Rock Mechanics and Geotechnical Engineering，2016，8 （3）： 415-422.

[55] GITAU A N，GUMBE L O，BIAMAH E K. Influence of soil water on stress-strain behaviour of a compacting soil in semi-arid Kenya [J]. Soil and Tillage Research， 2006，89（2）： 144-154.

[56]王麗琴，鹿忠剛，邵生俊.巖土體復合冪-指數(shù)非線性 模型[J].巖石力學與工程學報，2017，36（5）：1269-1278. WANG L Q，LU Z G， SHAO S J. A composite power exponential nonlinear model of rock and soil[J]. Chinese Journal of Rock Mechanics and Engineering，2O17，36 （5）：1269-1278.（in Chinese）

[57]王偉，宋新江，凌華，等.濱海相軟土應力-應變曲線復 合指數(shù)-雙曲線模型[J].巖土工程學報，2010，32（9）： 1455-1459. WANG W， SONG X J，LING H，et al. Composite exponential-hyperbolic model for stress-strain curve of seashore soft soil [J]. Chinese Journal of Geotechnical Engineering，2010，32（9）：1455-1459.（in Chinese）

[58]UCHAIPICHAT A，KHALILI N. Experimental investigation of thermo-hydro-mechanical behaviour of an unsaturated silt [J]. Géotechnique， 2009，59（4）： 339-353.

[59]SOLEIMANBEIGIA，EDIL TB，BENSONCH. Effect of temperature on geotechnical properties of recycled asphalt shingle mixtures [J]. Journal of Geotechnical and Geoenvironmental Engineering，2015，141（2）： 04014097.

[60]ROSCOE K H，BURLAND JB.On the yielding of soils [J].Geotechnique，1958，8（1）： 22-53.

[61]HUECKEL T，PELLEGRINIR，DELOLMO C.A constitutive study of thermo-elasto-plasticity of deep carbonatic clays [J]. International Journal for Numerical and Analytical Methods in Geomechanics，1998，22（7）： 549-574.

[62］劉紅，陳琴梅，盧黎，等.考慮溫度影響的非關聯(lián)彈塑 性飽和黏土本構模型[J].土木與環(huán)境工程學報（中英 文），2020，42（4）： 53-59. LIU H，CHEN Q M，LU L，et al.Temperaturedependent non-associated elastic-plastic consititutive model forsaturated clay [J]. Journal of Civil and Environmental Engineering，2020，42（4）：53-59.

[63]ABUEL-NAGA H M，BERGADO D T，BOUAZZA A，et al. Thermomechanical model for saturated clays [J].Géotechnique，2009，59（3）： 273-278.

[64]CUIY J， SULTAN N，DELAGE P.A thermomechanical model for saturated clays [J]. Canadian Geotechnical Journal，2000，37（3）： 607-620.

[65]HONGPY，PEREIRAJM，CUIYJ，etal.Atwosurface thermomechanical model for saturated clays [J]. International Journal for Numerical and Analytical Methods in Geomechanics，2016，40（7）： 1059-1080.

[66]殷宗澤.一個土體的雙屈服面應力-應變模型[J].巖土 工程學報，1988，10（4）：64-71. YIN Z Z.A double yield surface stress-strain model of soil [J]. Chinese Journal of Geotechnical Engineering， 1988，10（4）： 64-71. （in Chinese）

[67]TANG A M，CUI Y J.Modelling the thermomechanical volume change behaviour of compacted expansive clays [J].Géotechnique，2009，59（3）： 185-195.

[68] ALONSO E E，VAUNAT J， GENS A. Modelling the mechanical behaviour of expansive clays [J]. Engineering Geology，1999，54（1/2）： 173-183.

[69]KROHN K P. New conceptual models for the resaturation of bentonite [J]. Applied Clay Science，2003，23： 25-33.

[70] GRAHAM J， TANAKA N，CRILLY T，et al. Modified Cam-Clay modelling of temperature effects in clays [J]. Canadian Geotechnical Journal， 2Ool，38（3）： 608-621.

[71]HAMIDI A，KHAZAEI C.A thermo-mechanical con stitutive model for saturated clays [J]. International Jour nal of Geotechnical Engineering，20l0，4（4）： 445-459.

[72]HAMIDI A，TOURCHI S，KHAZAEI C. Thermomechanical constitutive model for saturated clays based on critical state theory [J]. International Journal of Geome chanics，2015，15（1）： 04014038.

[73] WANG L Z， WANG K J， HONG Y. Modeling temper ature-dependent behavior of soft clay [J]. Journal of Engi neering Mechanics，2016，142（8）： 04016054.

[74] ZHOU C，NG C W W.A thermomechanical model for saturated soil at small andlarge strains [J]. Canadian Geotechnical Journal， 2015，52（8）： 1101-1110.

[75]DAFALIAS Y F. Bounding surface plasticity. I： Mathematical foundation and hypoplasticity [J]. Journal of Engineering Mechanics，1986，112（9）： 966-987.

[76] ANANDARAJAH A， DAFALIAS Y F. Bounding surface plasticity. II Application to anisotropic cohesive soils [J]. Journal of Engineering Mechanics，1986，112 （12）： 1292-1318.

[77]DAFALIAS Y F，HERRMANN L R. Bounding surface plasticity. Ⅱ Ψ：Π Application to isotropic cohesive soils [J]. Journal of Engineering Mechanics，1986，112 （12）： 1263-1291.

[78] ZHOU C，NG C W W. Simulating the cyclic behaviour of unsaturated soil at various temperatures using a bounding surface model [J]. Géotechnique，2Ol6，66（4）： 344-350.

[79] ZHOU C，F(xiàn)ONG KY，NG C W W.A new bounding surface model for thermal cyclicbehaviour[J]. International Journal for Numerical and Analytical Methodsin Geomechanics，2017，41（16）： 1656-1666.

[80］陳艷妮，楊仲軒.基于熱力學理論的超固結黏土邊界面 模型[J].巖土工程學報，2017，39（3）： 547-553. CHEN Y N，YANG Z X. Thermodynamics-based bounding surface model for overconsolidated clay [J]. Chinese Journal of Geotechnical Engineering，2017，39 （3）： 547-553. （in Chinese）

[81]姚仰平，田雨，劉林.三維各向異性砂土UH模型[J]. 工程力學，2018，35（3）：49-55. YAO Y P，TIAN Y，LIU L. Thre-dimensional aniso tropic UH model for sands[J]. Engineering Mechanics， 2018，35（3）： 49-55.（in Chinese）

[82]姚仰平，劉林，羅汀.砂土的UH模型[J].巖土工程學 報，2016，38（12）： 2147-2153. YAO Y P，LIU L，LUO T.UH model for sands [J]. Chinese Journal of Geotechnical Engineering，2O16，38 （12）： 2147-2153.（in Chinese）

[83]孔令明，姚仰平.考慮時間效應的 K₀ 各向異性UH模型 [J].巖土工程學報，2015，37（5）：812-820. KONG L M， YAO Y P. K₀ -anisotropic UH model considering time effects [J]. Chinese Journal of Geotechnical Engineering，2015，37（5）： 812-820.（in Chinese）

[84］姚仰平，楊一帆，牛雷.考慮溫度影響的UH模型[J]. 中國科學：技術科學，2011，41（2）：158-169. YAO Y P， YANG YF， NIU L. UH model considering temperature effect [J]. Scientia Sinica （Technologica）， 2011，41（2）： 158-169.（in Chinese）

[85］姚仰平，侯偉，周安楠.基于Hvorslev面的超固結土本 構模型[J].中國科學（E輯：技術科學），2007，37（11）： 1417-1429. YAOYP，HOUW，ZHOUAN.Constitutivemodel of overconsolidated soil based on Hvorslev surface [J]. Science in China （Series E：Technological Sciences）， 2007，37（11）： 1417-1429. （in Chinese）

[86]姚仰平，牛雷，楊一帆，等.考慮溫度影響的非飽和土 本構模型[J].巖土力學，2011，32（10）：2881-2888. YAOYP，NIUL，YANG YF，etal.Constitutive model for unsaturated clays considering temperature effects [J]. Rock and Soil Mechanics，20ll，32（10）： 2881-2888.（in Chinese）

[87]姚仰平，萬征，楊一帆，等.飽和黏土不排水剪切的熱 破壞[J].巖土力學，2011，32（9）：2561-2569. YAOY P，WAN Z，YANG YF，et al. Thermal failure for saturated clay under undrained condition [J]. Rock and Soil Mechanics，2011，32（9）： 2561-2569. （in Chinese）

[88]孔令明，姚仰平.超固結土熱黏彈塑性本構關系[J].巖 土力學，2015，36（Sup 1）： 1-8. KONG L M，YAO Y P. Thermo-visco-elastoplastic constitutive relation for overconsolidated clay [J].Rock and Soil Mechanics，2015，36（Sup 1）： 1-8.（in Chinese）

[89]YAO Y P，KONG L M，ZHOU A N，et al.Timedependent unified hardening model：Three-dimensional elastoviscoplastic constitutive model for clays [J]. Journal of Engineering Mechanics，2015，141（6）： 04014162.

[90]YAOYP，HOUW，ZHOUAN.UH model： Threedimensional unified hardening model for overconsolidated clays [J]. Géotechnique，2009，59（5）： 451-469.

[91]ZHANG Z C，CHENG X H. Simulation of nonisother mal consolidation of saturated soils based on a thermodynamic model [J]. The Scientific World Journal，2013， 2013：192163.

[92]趙成剛，張雪東，郭璇.土的本構方程與熱力學[J].力 學進展，2006，36（4）：611-618. ZHAO C G，ZHANG X D，GUO X. Constitutive equations of soils and thermodynamics [J]. Advances in Mechanics，2006，36（4）：611-618.（in Chinese）

[93]李志敏.基于能量耗散理論凍土本構模型的研究[J]. 低溫建筑技術，2009，31（8）：87-89. LI Z M. Study on elastic-plastic damage constitutive model of frozen soil based on energy dissipation theory [J].Low Temperature Architecture Technology，2009， 31（8）： 87-89. （in Chinese）

[94]魯曉剛，王卓，CUIYuWen，等.計算熱力學、計算動 力學與材料設計[J].科學通報，2013，58（35）：3656- 3664. LUXG，WANGZ，CUIYW，etal.Computational thermodynamics，computational dynamics and material design [J]. Chinese Science Bulletin， 2O13，58（35）： 3656-3664.（in Chinese）

[95]李亮，李彥.基于熱力學原理的混凝土熱-力耦合本構 模型[J].北京工業(yè)大學學報，2016，42（4）：554-560. LI L，LI Y.Thermo-mechanical coupling constitutive model of concrete based on thermodynamics [J]. Journal of Beijing University of Technology，20l6，42（4）： 554- 560.（in Chinese）

[96]何敏，馮孝鵬，李寧，等.飽和正凍土水熱力耦合模型 的改進[J].巖土工程學報，2018，40（7）：1212-1220. HE M，F(xiàn)ENG XP，LI N，et al.Improvement of coupled thermo-hydro-mechanical model for saturated freezing soil [J]. Chinese Journal ofGeotechnical Engineering，2018，40（7）：1212-1220.（in Chinese）

[97]張志超，程曉輝.飽和土非等溫固結和不排水剪切的熱 力學本構模型[J].巖土工程學報，2013，35（7）：1297- 1306. ZHANG Z C， CHENG X H. Thermodynamic constitutive model for non-isothermal consolidation and undrained shear behaviors of saturated soils [J].Chinese Journal of Geotechnical Engineering，2Ol3， 35（7）： 1297- 1306.（in Chinese）

[98] ZHANG Z C，CHENG X H.A thermo-mechanical coupled constitutive model for clay based on extended granularsolid hydrodynamics [J]. Computersand Geotechnics，2016，80：373-382.

[99]BAI B，YANG G C，LI T，et al. A thermodynamic constitutive model with temperature effect based on particle rearrangement for geomaterials [J]. Mechanics of Materials，2019，139：103180.

[100]YANG G C，BAI B.A thermodynamic model to simulate the thermo-mechanical behavior of fine-grained gassy soil[J]. Bulletin of Engineering Geology and the Environment，2020，79（5）： 2325-2339.

[101] WHEELER S J. A conceptual model for soils contain ing large gas bubbles [J]. Géotechnique，1988，38（3）： 389-397.

[102]YANG G C，LIU Y，CHEN P P.Thermodynamic modeling of stress-strain behavior of saturated sand considering temperature effect [J]. AIP Advances， 2021，11（12）： 125206.

[103]GOLCHIN A，VARDON P J，HICKS M A.A thermo-mechanical constitutive model for fine-grained soilsbased on thermodynamics[J]. International Journal of Engineering Science，2022，174： 103579.

[104]KURZ D，SHARMA J，ALFARO M，et al. Semiempirical elastic-thermoviscoplastic model for clay [J]. Canadian Geotechnical Journal，2016，53（1O）：1583- 1599.

[105] LALOUI L，CEKEREVAC C.Numerical simulation of the non-isothermal mechanical behaviour of soils[J]. Computers and Geotechnics，2008，35（5）： 729-745.

[106]LALOUI L， CEKEREVAC C. Non-isothermal plasticity model for cyclic behaviour of soils [J]. International Journal for Numerical and Analytical Methods in Geomechanics，2008，32（5）： 437-460.

[107]LIU E L，XING H L.A double hardening thermomechanical constitutive model for overconsolidated clays[J].Acta Geotechnica，20O9，4（1）： 1-6.

[108]沈珠江.粘土的雙硬化模型[J].巖土力學，1995，16 （1）： 1-8. SHEN Z J. A double hardening model for clays [J]. Rock and Soil Mechanics，1995，16（l）：1-8.（in Chinese）

[109]劉祖典.巖土材料的簡易雙硬化模型[J].水資源與水 工程學報，1993，4（2）：88-96. LIU Z D. Simple double hardening model of geotechni cal materials [J]. Journal of Water Resources and Water Engineering，1993，4（2）： 88-96.（in Chinese）

[110]姚仰平，謝定義，俞茂宏.復雜應力下砂土的廣義雙 剪應力破壞準則及雙硬化本構模型[J].西安冶金建筑 學院學報，1994，26（4）：392-397. YAO Y P，XIED Y，YUM H. Generalized twin shearstress failure criterion and double-hardening constitutivc model of sand under the complex stress state [J]. Journal of Xi'an Institutte of Metallurgy and Constructural Engineering，1994，26（4）： 392-397.（in Chinese）

[11] MA Q J， NG C WW， MASIN D， et al. An approach for modelling volume change of fine-grainedsoil subjected to thermal cycles [J]. Canadian Geotechnical Journal，2017，54（6）： 896-901.

[112] COLLINS I F，BOULBIBANE M. Geomechanical analysis of unbound pavements based on shakedown theory [J]. Journal of Geotechnical and Geoenvironmental Engineering，2000，126（1）：50-59.

[113] HABIBALLAH T，CHAZALLON C. An elastoplastic model based on the shakedown concept for flexible pavements unbound granular materials [J]. International Journal for Numerical and Analytical Methods in Geomechanics，2005，29（6）： 577-596.

[114] NOWAMOOZ H，LI K，CHAZALLON C. Shakedown modeling of unsaturated expansive soils subjected to wetting and drying cycles [J]. E3S Web of Conferences，2016，9： 08007.

[115]MASIN D.A hypoplastic constitutive model for clays [J].International Journal for Numerical and Analytical Methods in Geomechanics，2005，29（4）： 311-336.

[116] MASiN D，KHALILI N. A thermo-mechanical model for variably saturated soils based on hypoplasticity [J]. International Journal for Numerical and Analytical Methods in Geomechanics，2012，36（12）： 1461-1485.

[117] ZHANG S，LENG W M， ZHANG F， et al. A simple thermo-elastoplasticmodelforgeomaterials[J]. International Journal of Plasticity，2Ol2，34： 93-113.

[118］張升，賀佐躍，滕繼東，等.考慮結構性的軟巖熱彈塑 性本構模型研究[J].巖石力學與工程學報，2017，36 （3）： 571-578. ZHANG S，HE ZY，TENGJD，etal.A thermoelasto-plastic model for soft rocks considering structure [J].Chinese Journal of Rock Mechanics and Engineer ing，2017，36（3）： 571-578.（in Chinese）

[119] ZHANG S， XU S，TENG JD，et al. Effect of temperature on the time-dependent behavior of geomaterials [J]. ComptesRendus Mecanique， 20l6，344（8）： 603-611.

[120]LALOUI L，LEROUEIL S，CHALINDARS. Modelling the combined effect of strain rate and temperature on one-dimensional compression of soils [J].Canadian Geotechnical Journal， 2OO8， 45（12）： 1765-1777.

[121]ROSCOEKH，SCHOFIELD AN，THURAIRAJAHA.Yieldingofclaysin stateswetterthan critical [J].Géotechnique，1963，13（3）：211-240.

[122]陸浩杰，孔綱強，劉漢龍，等.黏土熱-力學特性對能 量樁力學特性的影響[J].巖土工程學報，2022，44（1）： 53-61. LUHJ，KONGGQ，LIUHL，etal.Influencesof thermo-mechanical properties of clay on mechanical responses of energypiles[J].Chinese Journalof GeotechnicalEngineering，2022，44（1）：53-6l.（in Chinese）

[123]馮興，姚仰平，李汝寧，等.考慮溫度UH模型的有限 元應用[J].巖土工程學報，2015，37（Sup2）：181-185. FENG X，YAO YP，LIR N，et al.Application of UH model considering temperature to finite element method[J].Chinese Journal of Geotechnical Engineer ing，2015，37（Sup 2）：181-185.（in Chinese）

[124]王浩，張志超，程曉輝.基于熱力學理論的水泥石熱- 力耦合本構模型和有限元分析[J].巖石力學與工程學 報，2018，37（1）：67-76. WANGH，ZHANGZC，CHENGXH.A thermalmechanical constitutivemodel forcement rockbased on thermodynamics and its finite elementapplication [J]. Chinese Journal of Rock Mechanics and Engineering， 2018，37（1）：67-76.（in Chinese）

（編輯王秀玲）