吳海清

摘要：將液化后的飽和砂土視為一種流體，其流動(dòng)特性表現(xiàn)出類似剪切稀化非牛頓流體的特性。根據(jù)非牛頓流體力學(xué)理論分析液化后飽和砂土的流動(dòng)特性發(fā)現(xiàn)，當(dāng)處于零有效應(yīng)力狀態(tài)時(shí)，可以用純粘性流動(dòng)本構(gòu)模型較好地描述砂土的剪應(yīng)變率表觀黏度關(guān)系。通過(guò)擬合動(dòng)扭剪試驗(yàn)結(jié)果的流變曲線，對(duì)比分析了幾種常用的純粘性流動(dòng)本構(gòu)模型，發(fā)現(xiàn)Gross模型可以較好且較簡(jiǎn)潔地描述液化后砂土零有效應(yīng)力狀態(tài)時(shí)的流動(dòng)特性。通過(guò)對(duì)擬合參數(shù)的分析，建立了基于Gross模型的液化后飽和砂土零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型，得到時(shí)間量綱下的參數(shù)K以及零剪切表觀黏度η0和極限剪切表觀黏度η∞的函數(shù)關(guān)系?？紤]相對(duì)密實(shí)度、固結(jié)應(yīng)力和應(yīng)力歷史對(duì)模型參數(shù)的影響，闡釋了模型參數(shù)的函數(shù)關(guān)系以及物理意義。

關(guān)鍵詞：飽和砂土;剪切稀化;零有效應(yīng)力;流動(dòng)本構(gòu)模型;剪應(yīng)變率表觀黏度

中圖分類號(hào)：TU447 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：20966717（2020）03003208

Abstract：

When the liquefied watersaturated sand is regarded as a fluid， it can be idealized as a kind of shear thinning nonNewtonian fluid. Based on nonNewtonian fluid mechanics， the theoretical analysis was conducted on the liquefied watersaturated sand. According to this， the relationship between the shear strain rate and the viscosity under zero effective stress state can be properly described by pure viscous flow constitutive model. The Gross model of liquefied watersaturated sand is established by fitting the flow curve of the dynamic torsional shear test results on the basis of predictions obtained from several commonly used pure viscous flow constitutive models. The representative model parameters include the relative density， the consolidation stress and the impact of stress history. The functional relationships among time dimension K， zero shear apparent viscosity η0 and limit shear apparent viscosity η∞ are fitted based on the dynamic torsional shear test results， and the physical meanings of the parameters are introduced in the composed function. The research landmarks a theoretical foundation for further analysis of large deformation of liquefied flows based on fluid mechanics theory.

Keywords：watersaturated sand; shear thinning; zero effective stress; flow constitutive model; shear strain rateapparent viscosity

地震是人類面臨的主要自然災(zāi)害之一。近年來(lái)，世界范圍內(nèi)地震活動(dòng)頻繁，如2008年中國(guó)“5·12”汶川大地震（里氏震級(jí)M=8.0）[1]，2010年智利“2·27”大地震（里氏震級(jí)M=8.8）以及2011年日本“3·11”本州島海域大地震（里氏震級(jí)M=9.0）[2]，2018年印尼“9·28”大地震（里氏震級(jí)M=7.4）等，這些大地震中均存在因飽和砂土液化而引起的震害，因此，地震液化問(wèn)題的研究仍然是土動(dòng)力學(xué)與巖土地震工程領(lǐng)域重要的課題之一。

Seed等[3]把不排水循環(huán)剪切試驗(yàn)中有效應(yīng)力第一次為0的狀態(tài)稱為“初始液化”，從而將液化過(guò)程分為“液化前（初始液化前）”和“液化后（初始液化后）”兩個(gè)階段。已有的震害調(diào)查發(fā)現(xiàn)，砂土層“液化后”引起的地基大變形是導(dǎo)致強(qiáng)震區(qū)各種基礎(chǔ)設(shè)施和生命線工程遭到破壞的主要原因，因此，隨著液化研究的深入，學(xué)者們?cè)絹?lái)越關(guān)注“液化后”砂土層應(yīng)力應(yīng)變的響應(yīng)以及“液化后”地基大變形產(chǎn)生的機(jī)理。

根據(jù)力學(xué)理論的不同，目前對(duì)砂土液化后大變形本構(gòu)模型的研究主要有兩類：一類基于固體力學(xué)理論，如Shamoto等[4]提出描述砂土液化后不排水剪單調(diào)剪切大變形的本構(gòu)模型;另一類將“液化后”的砂土層視為一種流體，基于流體力學(xué)的理論研究大變形的問(wèn)題[5]。

將液化后砂土視為一種流體是一種較新的研究思路。Sasaki等[6]利用振動(dòng)臺(tái)試驗(yàn)研究砂土液化變形問(wèn)題時(shí)發(fā)現(xiàn)，液化砂土具有非常類似于流體的性質(zhì);Miyajima等[7]采用拖球及落球試驗(yàn)，根據(jù)流體力學(xué)中的公式獲得液化后砂土的表觀黏度。Tamate等[8]在拖管試驗(yàn)中，根據(jù)鋼管受到的拖拽力和相似性，研究液化后砂土的流動(dòng)特性。Hadush等[9]總結(jié)了幾組由不同測(cè)量方法獲得的表觀黏度與剪應(yīng)變率關(guān)系后，發(fā)現(xiàn)液化后砂土的表觀黏度隨剪應(yīng)變率的增大而減小的特性。Hwang等[10]等通過(guò)沉球試驗(yàn)及拖管試驗(yàn)研究液化后地表流對(duì)樁基的影響時(shí)，同樣發(fā)現(xiàn)液化砂土的這種特性。陳育民等[5， 1112]通過(guò)動(dòng)扭剪試驗(yàn)和振動(dòng)臺(tái)拖球試驗(yàn)，王志華等[13]通過(guò)自由場(chǎng)地基振動(dòng)臺(tái)試驗(yàn)，其他一些學(xué)者[1420]通過(guò)空心扭剪實(shí)驗(yàn)、離心機(jī)振動(dòng)臺(tái)實(shí)驗(yàn)等方法，研究液化后砂土流動(dòng)特性時(shí)，也獲得了同樣的結(jié)果，即液化后砂土是一種“剪切稀化非牛頓流體”。由于對(duì)于液化后砂土流體動(dòng)力學(xué)特性的研究起步較晚，理論研究相對(duì)缺乏，在流動(dòng)本構(gòu)模型方面的研究尚處于初步探索階段[14]。

筆者采用流體力學(xué)理論進(jìn)一步解釋砂土液化后大變形的機(jī)理，通過(guò)對(duì)動(dòng)扭剪試驗(yàn)結(jié)果的分析，探討可以描述液化后砂土零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型。

1流體力學(xué)方法分析液化后砂土流動(dòng)特性的理論基礎(chǔ)

1687年，牛頓在其所著的《自然哲學(xué)的數(shù)學(xué)原理》中提出牛頓內(nèi)摩擦定律

式中：μ為比例系數(shù)，稱為動(dòng)力黏度（簡(jiǎn)稱黏度），單位是Pa·S，該系數(shù)即反映了流體的粘稠性。dudy為流速在流層法線方向的變化率，稱為流速梯度。

由于液化砂土的流速梯度很難獲得，因此，需要進(jìn)一步說(shuō)明該項(xiàng)的物理意義。如圖1所示，在厚度為dy的上、下兩流層間取矩形流體微團(tuán)。因上、下層的流速相差du，經(jīng)過(guò)時(shí)間dt后，微團(tuán)除位移外，還有剪切應(yīng)變dγ，即

陳育民等[5]通過(guò)不同試驗(yàn)發(fā)現(xiàn)，液化后砂土的流動(dòng)特性表現(xiàn)為“非牛頓流體”的性質(zhì)。而只有牛頓流體才具有一種嚴(yán)格意義的黏度概念，即μ。非牛頓流體則需要兩個(gè)或兩個(gè)以上的參數(shù)來(lái)描述其黏稠特性。但為了方便起見(jiàn)，引入動(dòng)力表觀黏度（簡(jiǎn)稱表觀黏度）η來(lái)近似描述非牛頓流體的黏稠特性，即η=f（γ·）（γ·=dγdt），為方便起見(jiàn)，采用γ·代替dγdt）。于是式（4）可改寫為

由于液化后砂土的表觀黏度η隨剪應(yīng)變率γ·的增大而減小，即表現(xiàn)為一種“剪切稀化”的特性。剪切稀化非牛頓流體材料表現(xiàn)為：在低剪應(yīng)變率γ·0下，表觀黏度為常數(shù)或趨向于一個(gè)常數(shù)，該常數(shù)η0稱為“零剪切表觀黏度”;與此同時(shí)，在高剪應(yīng)變率γ·∞下，表觀黏度再次為常數(shù)或趨向于一個(gè)常數(shù)（且η∞<η0），該常數(shù)η∞稱為“極限剪切表觀黏度”[22]。

可以推斷，砂土液化后處于零有效應(yīng)力狀態(tài)時(shí)，顆粒間失去土骨架結(jié)構(gòu)，處于類似“懸浮”的狀態(tài)，當(dāng)沒(méi)有受到剪應(yīng)力作用時(shí)，顆粒在溶液中處于靜止?fàn)顟B(tài)，由于顆粒本身在液化過(guò)程中受到動(dòng)荷載的作用，處于雜亂無(wú)章的分布狀態(tài)，一旦受到剪應(yīng)力的作用，隨著流動(dòng)的進(jìn)行，它們將沿著流動(dòng)方向開始有序地排列起來(lái)，此時(shí)剪應(yīng)變率開始增大，這種定向排列將隨著流動(dòng)的進(jìn)行，越來(lái)越整齊，流動(dòng)的阻力就越來(lái)越小，即其表觀黏度將越來(lái)越小，這就反映出一種剪切稀化的趨勢(shì)。當(dāng)剪應(yīng)變率較小時(shí)，不足以破壞原有的結(jié)構(gòu)，即不能夠使排列得雜亂無(wú)章的顆粒定向排列，此時(shí)表觀黏度為常數(shù)η0且與剪應(yīng)變率無(wú)關(guān)，而當(dāng)剪應(yīng)變率很大時(shí)，已經(jīng)最大限度地使顆粒定向排列，從而顆粒已經(jīng)處于有序整齊的狀態(tài)，此時(shí)再增大剪應(yīng)變率，表觀黏度也無(wú)法再小了，即趨向于一個(gè)常數(shù)η∞，且η∞<η0。

2非牛頓流體的流動(dòng)本構(gòu)模型

描述流體剪切應(yīng)力和剪應(yīng)變率之間關(guān)系的方程，稱為流動(dòng)本構(gòu)模型。當(dāng)研究非牛頓流體時(shí)，流動(dòng)本構(gòu)模型即式（5），由于表觀黏度η受到非牛頓流體類型以及環(huán)境等眾多因素的影響，表現(xiàn)為關(guān)于剪應(yīng)變率γ·的非線性函數(shù)。

由于目前還沒(méi)有一個(gè)普適的流動(dòng)本構(gòu)模型可以統(tǒng)一描述非牛頓流體的流動(dòng)特性。因此，在研究非牛頓流體剪應(yīng)變率表觀黏度之間的關(guān)系時(shí)，存在著許多流動(dòng)本構(gòu)模型。歸納起來(lái)可以分為兩大類：一類是純粘性流動(dòng)本構(gòu)模型，即在撤除剪切應(yīng)力后，它們?cè)谑芗羟袘?yīng)力作用期間的任何形變都不會(huì)恢復(fù);另一類是粘彈性流動(dòng)本構(gòu)模型，即在撤除剪切應(yīng)力后，它們?cè)谑芗羟袘?yīng)力作用期間所產(chǎn)生的形變會(huì)完全或部分得到恢復(fù)。

張建民等[23]在研究砂土不排水循環(huán)扭剪試驗(yàn)時(shí)發(fā)現(xiàn)：初始液化后零有效應(yīng)力狀態(tài)時(shí)產(chǎn)生的剪應(yīng)變分量γ0的變化特點(diǎn)為：γ0只在液化后有效應(yīng)力過(guò)零點(diǎn)時(shí)產(chǎn)生;γ0的大小與當(dāng)前的剪應(yīng)力大小無(wú)關(guān)，只依賴于應(yīng)力應(yīng)變歷史。由于筆者研究的對(duì)象為液化后處于零有效應(yīng)力狀態(tài)的砂土，只需要考慮剪應(yīng)變分量γ0，即在此狀態(tài)下，砂土不存在土骨架的概念，因此，該研究對(duì)象屬于純粘性流體。

由于液化后砂土的流動(dòng)特性比較接近側(cè)黏流動(dòng)，即在流動(dòng)過(guò)程中只需要考慮平移和轉(zhuǎn)動(dòng)，不需要考慮拉伸等其他因素。因此Π=γ·，即描述液化后砂土流動(dòng)特性的流動(dòng)本構(gòu)模型可以寫成η=η（γ·）的形式，流動(dòng)本構(gòu)模型的書寫將采用這種統(tǒng)一的形式。

3基于動(dòng)扭剪試驗(yàn)結(jié)果分析液化后砂土的流動(dòng)本構(gòu)模型

3.1試驗(yàn)方法

試驗(yàn)采用河海大學(xué)巖土所和日本圓井株式會(huì)社共同研制的靜動(dòng)多功能三軸試驗(yàn)儀。采用的砂樣為重塑試樣（顆粒比重G1為2.56，有效粒徑d10=0.16 mm，限制粒徑d60為0.42 mm，最大孔隙比emax為0.974，最小孔隙比emin為0.568，不均勻系數(shù)Cu為2.63，曲率系數(shù)Cc為1.17，粘粒含量為0%），試樣采用空中砂雨法制備。試驗(yàn)中，采用先通CO2氣體排除試樣中空氣，再通脫氣水，然后再施加反壓的方法進(jìn)行飽和，確保飽和度B≥0.97。試樣完成固結(jié)操作后，在不排水條件下以一定的動(dòng)剪應(yīng)力比給試樣施加動(dòng)扭剪應(yīng)力，當(dāng)試樣發(fā)生初始液化后，繼續(xù)施加動(dòng)扭剪應(yīng)力，以此來(lái)模擬液化后砂土的流動(dòng)特性。試驗(yàn)中考慮了不同初始相對(duì)密實(shí)度和前期固結(jié)應(yīng)力對(duì)流動(dòng)特性的影響，具體方案如表2所示。

3.2分析方法

根據(jù)流體力學(xué)中常用的建立流動(dòng)本構(gòu)模型的方法，即采用實(shí)驗(yàn)數(shù)據(jù)建立流變曲線進(jìn)行參數(shù)擬合。為了確定描述液化后砂土非零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型，將根據(jù)動(dòng)扭剪試驗(yàn)結(jié)果繪制成的流變曲線作為研究對(duì)象。根據(jù)所采用流動(dòng)本構(gòu)模型的形式，流變曲線采用剪應(yīng)變率表觀黏度之間的關(guān)系曲線。剪應(yīng)變率和表觀黏度均由試驗(yàn)數(shù)據(jù)計(jì)算獲得，公式為[11]

液化后動(dòng)加載時(shí)，每一個(gè)荷載周期包括兩次加載和兩次卸載過(guò)程，根據(jù)以往學(xué)者的研究，只有在加載階段超孔壓比FLFL=Δuσ3c′接近1.0時(shí)，液化后砂土才會(huì)在極低的剪應(yīng)力作用下使得剪應(yīng)變迅速增長(zhǎng)，即表現(xiàn)出近似粘性流體的特性。此外，文本選擇的研究對(duì)象為液化后零有效應(yīng)力狀態(tài)下的砂土也滿足了這一特性。因此，將每組試驗(yàn)中超孔壓比接近1.0（這里取大于0.99）時(shí)的數(shù)據(jù)點(diǎn)篩選出來(lái)，并借鑒陳育民等[11]的觀點(diǎn)，將1個(gè)周期的兩次加載過(guò)程分為“第一半周加載”和“第二半周加載”，具體的數(shù)據(jù)組如表3所示。

3.3剪應(yīng)變率表觀黏度曲線擬合

將試驗(yàn)結(jié)果中超孔壓比接近1.0的液化后砂土剪應(yīng)變率表觀黏度的關(guān)系采用不同的流動(dòng)本構(gòu)模型進(jìn)行曲線擬合，繪出流變曲線。采用的數(shù)據(jù)組為No.9、No.10、No.13、No.14、No.25和No.26。

圖2為采用不同流動(dòng)本構(gòu)模型進(jìn)行曲線擬合的結(jié)果，從圖2可以發(fā)現(xiàn)，Carreau方程和Gross方程能較好地?cái)M合出試驗(yàn)結(jié)果。盡管兩個(gè)方程均為4參數(shù)方程，但Carreau方程在形式上比Gross方程更加復(fù)雜，因此，考慮在工程上本構(gòu)模型越簡(jiǎn)單越方便應(yīng)用的原則，選擇Gross方程作為液化后砂土零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型。

3.4Gross方程曲線擬合的結(jié)果分析

采用Gross方程進(jìn)行曲線擬合，擬合的方式為最小二乘法。調(diào)整相關(guān)度后，擬合度均達(dá)到了0.99，說(shuō)明擬合出來(lái)的曲線和實(shí)驗(yàn)數(shù)據(jù)點(diǎn)之間具有良好的相關(guān)度。

圖3為不同相對(duì)密實(shí)度下Gross方程的擬合結(jié)果，從圖中可知，在同一個(gè)周期的兩次加載過(guò)程中，“第一半周加載”和“第二半周加載”為兩個(gè)獨(dú)立的加載過(guò)程，剪應(yīng)變率表觀黏度關(guān)系不一致。主要原因是，兩個(gè)加載過(guò)程的方向相反，應(yīng)力應(yīng)變歷史不一致，因此，有必要把兩個(gè)加載過(guò)程分開分析。

從圖3（a）中擬合曲線可知，當(dāng)剪應(yīng)變率相同時(shí)，相對(duì)密實(shí)度不同的液化后砂土，其表觀黏度不同，其中，相對(duì)密實(shí)度為40%的液化后砂土表觀黏度大于相對(duì)密實(shí)度為30%的液化后砂土，反映了松砂的表觀黏度隨相對(duì)密實(shí)度的提高而增加，即提高相對(duì)密實(shí)度可以增強(qiáng)疏松砂層液化后抵抗變形的能力。產(chǎn)生這一現(xiàn)象的原因是：松砂在振動(dòng)受剪時(shí)，顆?？梢詽L落到平衡位置，使得排列更加緊密，表現(xiàn)出剪縮性，當(dāng)液化后砂土流動(dòng)時(shí)，砂顆粒會(huì)重新排列，趨于整齊，砂顆粒越緊密，趨于整齊的過(guò)程中受到的流動(dòng)阻力就越大。相對(duì)密實(shí)度為50%的液化后砂土，其表觀黏度與相對(duì)密實(shí)度的關(guān)系沒(méi)有反映出松砂的特性，產(chǎn)生這一現(xiàn)象的原因是：緊砂在液化后大變形的過(guò)程中，主要體現(xiàn)了剪脹性，即緊砂在振動(dòng)受剪時(shí)，顆粒只能通過(guò)抬高或者錯(cuò)位來(lái)離開原有位置，表現(xiàn)為剪脹性，當(dāng)液化后砂土流動(dòng)時(shí)，砂顆粒同樣會(huì)重新排列，趨于整齊，趨于整齊的過(guò)程中受到的流動(dòng)阻力隨著砂顆粒的齊整程度而減小。圖3（b）中也存在類似的情況。

由圖3（c）中擬合曲線的走勢(shì)可知，在相對(duì)密實(shí)度和固結(jié)應(yīng)力相同時(shí)，同一剪應(yīng)變率對(duì)應(yīng)的表觀黏度隨著循環(huán)周期的增長(zhǎng)而增長(zhǎng)，圖3（d）中也存在類似的情況。產(chǎn)生這一現(xiàn)象的原因是：液化后零有效應(yīng)力狀態(tài)下砂土的流動(dòng)過(guò)程即是砂顆粒重新排列整齊的過(guò)程，進(jìn)行這一過(guò)程前，就已經(jīng)存在過(guò)一次砂顆粒排列整齊后再被打亂的循環(huán)，本次重新排列過(guò)程勢(shì)必受到上一次的影響，并且這樣的影響是一種累加的過(guò)程，即建立在前一次循環(huán)基礎(chǔ)上的后一次重新排列過(guò)程要更困難，張建民等[17]的研究也發(fā)現(xiàn)，液化后零有效應(yīng)力狀態(tài)時(shí)，砂土經(jīng)歷的前期最大剪應(yīng)變?cè)酱?，此后發(fā)生的剪應(yīng)變也越大。此外，從圖中相同條件下不同循環(huán)周期的擬合曲線分布較密可知，這種應(yīng)力應(yīng)變歷史導(dǎo)致表觀黏度的增長(zhǎng)幅度不大。

圖4為不同固結(jié)應(yīng)力下Gross方程的擬合結(jié)果，從擬合曲線的走勢(shì)可知，在剪應(yīng)變率相同時(shí)，同一相對(duì)密實(shí)度的液化后砂土表觀黏度在固結(jié)應(yīng)力為100 kPa時(shí)比在固結(jié)應(yīng)力為50 kPa時(shí)大，即相同條件下，固結(jié)應(yīng)力越大，液化后砂土流動(dòng)時(shí)受到的阻力越大，產(chǎn)生這一現(xiàn)象的原因是：前期固結(jié)應(yīng)力越大，砂顆粒之間越密實(shí)，當(dāng)液化后砂土開始流動(dòng)，砂顆粒開始重新排列整齊時(shí)，如果砂顆粒之間接觸得越緊密，則這種重新排列的過(guò)程越困難，從而流動(dòng)阻力增大。

3.5液化后砂土基于Gross方程的本構(gòu)模型探討

由式（12）可知，Gross方程中包含4個(gè)參數(shù)：零剪切表觀黏度η0，極限剪切表觀黏度η∞，時(shí)間量綱下的參數(shù)K和無(wú)量綱參數(shù)m。這些參數(shù)在描述液化后砂土零有效應(yīng)力狀態(tài)下的流動(dòng)特性時(shí)，具體表現(xiàn)為：

1）無(wú)量綱參數(shù)m不受初始相對(duì)密實(shí)度、前期固結(jié)應(yīng)力、同一荷載周期內(nèi)加載方向以及循環(huán)荷載周期的影響。

2）時(shí)間量綱下的參數(shù)K，單位為s，受初始相對(duì)密實(shí)度、前期固結(jié)應(yīng)力以及同一荷載周期內(nèi)加載方向的影響，但與循環(huán)荷載周期無(wú)關(guān)。在試驗(yàn)結(jié)果中，當(dāng)處于“第一半周加載”時(shí)，參數(shù)K隨前期固結(jié)應(yīng)力單調(diào)遞減，初始相對(duì)密實(shí)度越高，參數(shù)K越小;當(dāng)處于“第二半周加載”時(shí)，參數(shù)K隨前期固結(jié)應(yīng)力單調(diào)遞減，初始相對(duì)密實(shí)度越高，參數(shù)K越大。參數(shù)K的表達(dá)式可以寫為

3）零剪切表觀黏度η0和極限剪切表觀黏度η∞，單位為kPa·s，受初始相對(duì)密實(shí)度、前期固結(jié)應(yīng)力、同一荷載周期內(nèi)加載方向以及循環(huán)荷載周期的影響，在試驗(yàn)結(jié)果中，當(dāng)處于“第一半周加載”時(shí)，η0和η∞隨循環(huán)荷載周期單調(diào)遞增;η0和η∞隨固結(jié)應(yīng)力單調(diào)遞增;η0和η∞隨相對(duì)密實(shí)度單調(diào)遞增。當(dāng)處于“第二半周加載”，η0和η∞隨循環(huán)荷載周期的增加在低固結(jié)應(yīng)力時(shí)單調(diào)遞增，在高固結(jié)應(yīng)力時(shí)單調(diào)遞減（由于遞減的幅度很小而且數(shù)據(jù)組數(shù)量有限，這里考慮到試驗(yàn)精度以及物理性質(zhì)的統(tǒng)一性，函數(shù)關(guān)系依然采用遞增關(guān)系）;η0和η∞隨固結(jié)應(yīng)力單調(diào)遞增;η0和η∞隨相對(duì)密實(shí)度單調(diào)遞增。η0和η∞的表達(dá)式可以寫為

由于篩選出的試驗(yàn)組數(shù)量有限，式（17）和式（18）具體的形式還有待更多的試驗(yàn)數(shù)據(jù)來(lái)進(jìn)行擬合推導(dǎo)。

4結(jié)論

開展了液化后砂土零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型研究，得出以下主要結(jié)論：

1）液化后砂土在零有效應(yīng)力狀態(tài)時(shí)，表現(xiàn)為一種固液混合流體，其流動(dòng)特性可以歸屬于混合粘性流體的研究范疇，這種混合流體結(jié)構(gòu)導(dǎo)致其剪應(yīng)變率表觀黏度的流動(dòng)曲線上出現(xiàn)了“零剪切表觀黏度”η0和“極限剪切表觀黏度”η∞。

2）對(duì)動(dòng)扭剪試驗(yàn)結(jié)果的分析表明：液化后砂土在零有效應(yīng)力狀態(tài)下，相對(duì)密實(shí)度較低時(shí)，即松砂相對(duì)密實(shí)度越高，液化后砂土流動(dòng)時(shí)表觀黏度越大，相對(duì)密實(shí)度較高時(shí)，則反映了緊砂的剪脹性;液化后砂土流動(dòng)過(guò)程中受到的阻力隨著受剪循環(huán)周期的增加而增加;前期固結(jié)應(yīng)力越高，液化后砂土流動(dòng)過(guò)程中受到的阻力越大。

3）對(duì)動(dòng)扭剪試驗(yàn)結(jié)果的擬合表明，Gross方程可以較好且相對(duì)比較簡(jiǎn)潔地描述液化后砂土零有效應(yīng)力狀態(tài)時(shí)的流動(dòng)特性，其方程中的參數(shù)在加載過(guò)程的“第一半周加載”和“第二半周加載”是相互獨(dú)立的，并且通過(guò)對(duì)擬合參數(shù)的分析，得到了時(shí)間量綱下的參數(shù)K以及零剪切表觀黏度η0和極限剪切表觀黏度η∞的函數(shù)關(guān)系。

動(dòng)扭剪試驗(yàn)獲得的有效數(shù)據(jù)數(shù)量有限，為獲得液化后砂土零有效應(yīng)力狀態(tài)的流動(dòng)本構(gòu)模型的具體形式，還需要開展更多的室內(nèi)試驗(yàn)進(jìn)行驗(yàn)證和校正。參考文獻(xiàn)：

[1] 袁曉銘， 曹振中. 汶川大地震液化的特點(diǎn)及帶來(lái)的新問(wèn)題[J]. 世界地震工程， 2011， 27（1）： 18.

YUAN X M，CAO Z Z. Features and new aspects of liquefaction in the Wenchuan Earthquake [J]. World Earthuqake Engineeering， 2011， 27（1）： 18. （in Chinese）

[2] ?VERDUGO R. Comparing liquefaction phenomena observed during the 2010 maule， Chile earthquake and 2011 great east Japan earthquake [C] // Proceedings of the International Symposium， Tokyo， Japan， 2012.

[3] ?SEED H B， LEE K L. Liquefaction of saturated sands during cyclic loading [J]. Journal of Soil Mechanics and Foundation Engineering Division， ASCE， 1966， 92（6）： 105134.

[4] ?SHAMOTO Y， ZHANG J M， GOTO S. Mechanism of large postliquefaction deformation in saturated sand [J]. Soils and Foundations， 1997， 37（2）： 7180.

[5] 陳育民， 劉漢龍， 周云東. 液化及液化后砂土的流動(dòng)特性分析[J]. 巖土工程學(xué)報(bào)， 2006， 28（9）： 11391143.

CHEN Y M， LIU H L， ZHOU Y D. Analysis on flow characteristics of liquefied and postliquefied sand [J]. Chinese Journal of Geotechnical Engineering， 2006， 28（9）： 11391143. （in Chinese）

[6] ?SASAKI Y， TOWHATA I， TOKIDA K， et al. Mechanism of permanent displacement of ground caused by seismic liquefaction [J]. Soils and Foundations， 1992， 32（3）： 7996.

[7] ?MIYAJIMA M， KITAURA M， KOIKE T， et al. Experimental study on characteristics of liquefied ground flow[C] // The First International Conference on Earthquake Geotechnical Engineering， Balkema， 1995： 969974.

[8] ?TAMATE S， TOWHATA I. Numerical simulation of ground flow causedby seismic liquefaction [J]. Soil Dynamics and Earthquake Engineering， 1999， 18（7）： 473485.

[9] ?HADUSH S， YASHIMA A， UZUOKA R. Importance of viscous fluid characteristics in liquefaction induced lateral spreading analysis [J]. Computers and Geotechnics， 2000， 27（3）： 199224.

[10] ?HWANG J I， KIM C Y， CHUNG C K， et al. Viscous fluid characteristics of liquefied soils and behavior of piles subjected to flow of liquefied soils [J]. Soil Dynamics and Earthquake Engineering， 2006， 26（2/3/4）： 313323.

[11] 劉漢龍， 陳育民. 動(dòng)扭剪試驗(yàn)中砂土液化后流動(dòng)特性分析[J]. 巖土力學(xué)， 2009， 30（6）： 15371541.

LIU H L， CHEN Y M. Analysis of flow characteristics of dynamic torsional tests on post liquefied sand [J]. Rock and Soil Mechanics， 2009， 30（6）： 15371541. （in Chinese）

[12] 陳育民， 劉漢龍， 邵國(guó)建， 等. 砂土液化及液化后流動(dòng)特性試驗(yàn)研究[J]. 巖土工程學(xué)報(bào)， 2009， 31（9）： 14081413.

CHEN Y M， LIU H L， SHAO J G. Laboratory tests on flow characteristics of liquefied and postliquefied sand [J]. Chinese Journal of Geotechnical Engineering， 2009， 31（9）： 14081413. （in Chinese）

[13] 王志華， 周恩全， 陳國(guó)興. 孔壓增長(zhǎng)后的飽和砂土流體特性及其孔壓相關(guān)性[J]. 巖土工程學(xué)報(bào)， 2012， 34（3）： 528533.

WANG Z H， ZHOU E Q， CHEN G X. Fluid characteristics dependent on excess pore water pressure of saturated sand after growth of pore pressure [J]. Chinese Journal of Geotechnical Engineering， 2012， 34（3）： 528533. （in Chinese）

[14] 黃雨， 毛無(wú)衛(wèi). 液化后土體的流態(tài)化特性研究進(jìn)展[J]. 同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2011， 39（4）： 501506.

HUANG Y， MAO W W. State of art of fluidization behavior of postliquefied soils [J]. Journal of Tongji University（Natural Science）， 2011， 39（4）： 501506. （in Chinese）

[15] ?周葛. 緩傾砂土場(chǎng)地液化側(cè)向變形離心模型試驗(yàn)與預(yù)測(cè)方法研究[D]. 南京： 河海大學(xué)， 2016.

ZHOU G. Centrifuge tests and numerical predictions of liquefaction induced lateral spreading on mildly sloping sandy ground [D]. Nanjing： Hohai University， 2016. （in Chinese）

[16] ?李晶， 陳育民， 方志， 等. 減飽和砂土緩傾場(chǎng)地的液化性狀分析[J]. 巖土力學(xué)， 2019， 40（9）： 19.

LI J， CHEN Y M， FANG Z，et al. Liquefaction characteristics analysis on gently tilting desaturated sandy ground [J]. Rock and Soil Mechanics， 2019， 40（9）： 19. （in Chinese）

[17] 周葛， 劉漢龍， 井合進(jìn)， 等. 基于應(yīng)變空間多機(jī)構(gòu)CG模型的地震液化大變形數(shù)值分析[J]. 巖土力學(xué)， 2016， 37（8）： 24172425.

ZHOU G， LIU H L， IAI S， et al. Liquefactioninduced lateral spreading analysis in gently sloping ground from earthquakes based on CG model [J]. Rock and Soil Mechanics， 2016， 37（8）： 24172425. （in Chinese）

[18] ?周正龍， 陳國(guó)興， 吳琪. 四向振動(dòng)空心圓柱扭剪儀模擬主應(yīng)力軸旋轉(zhuǎn)應(yīng)力路徑能力分析[J]. 巖土力學(xué)， 2016， 37（Sup1）： 126132.

ZHOU Z L， CHEN G X， WU Q. Analysis of capabilities of stress paths of HCA to simulate principal stress rotation under fourdirection dynamic loads [J]. Rock and Soil Mechanics， 2016， 37（Sup1）： 126132. （in Chinese）

[19] 周正龍， 陳國(guó)興， 趙凱， 等. 循環(huán)加載方向角對(duì)飽和粉土不排水動(dòng)力特性的影響[J]. 巖土力學(xué)， 2018， 39（1）： 3644.

ZHOU Z L， CHEN G X， ZHAO K， et al. Effect of the direction angle of cyclic loading on undrained cyclic behavior of saturated silt [J]. Rockand Soil Mechanics， 2018， 39（1）： 3644. （in Chinese）

[20] ?方志， 陳育民， 何森凱. 基于單相流的減飽和砂土流固耦合改進(jìn)算法[J]. 巖土力學(xué)， 2018， 39（5）： 18511857.

FANG Z， CHEN Y M， HE S K. A modified method of mechanicalhydraulical coupled simulation based on singlephase fluid for desaturated sand [J]. Rock and Soil Mechanics， 2018， 39（5）： 17. （in Chinese）

[21] ?陳玉璞. 流體動(dòng)力學(xué)[M]. 南京： 河海大學(xué)出版社， 1990： 9.

CHEN Y P. Fluid dynamics [M]. Nanjing： Hohai University Press， 1990： 9. （in Chinese）

[22] ?韓式方. 非牛頓流體本構(gòu)方程和計(jì)算解析理論[M]. 北京： 科學(xué)出版社， 2000： 343.

HAN S F. Costitutive equation and computational analytical theory of NonNewtonian fluid [M]. Beijing： Science Press， 2000： 343. （in Chinese）

[23] 張建民， 王剛. 砂土液化后大變形的機(jī)理[J]. 巖土工程學(xué)報(bào)， 2006， 28（7）： 835840.

ZHANG J M， WANG G. Mechanism of large postliquefaction deformation in saturated sand [J]. Chinese Journal of Geotechnical Engineering， 2006， 28（7）： 835840. （in Chinese）

（編輯王秀玲）