不同建模假定下巖石損傷黏彈模型的參數(shù)轉(zhuǎn)換關(guān)系探討

方 燾，劉新榮，黃 明

（1.重慶大學(xué) 土木工程學(xué)院，重慶 400045；2.華東交通大學(xué) 土木建筑學(xué)院，南昌 330013；3.重慶大學(xué) 山地城鎮(zhèn)建設(shè)與新技術(shù)教育部重點實驗室，重慶 400045；4.福州大學(xué) 土木工程學(xué)院，福州 350108）

1 引 言

目前，針對巖石蠕變本構(gòu)模型的研究已經(jīng)取得了較大的進展，并獲得了較豐富的研究成果[1－3]。對于元件組合模型的應(yīng)用，單軸蠕變試驗時模型的彈性參數(shù)常選用拉壓模量，而在FLAC3D軟件中對三維工程進行模擬計算時一般采用剪切模量，因此，針對工程應(yīng)用或三維數(shù)值計算中的需要，有必要進行一維狀態(tài)下巖石蠕變模型的三維狀態(tài)推廣研究。通常情況下，蠕變方程由一維轉(zhuǎn)化到三維時存在不同的建模思路，不同的建模思路對應(yīng)的模型參數(shù)并不相同，需要建立轉(zhuǎn)換關(guān)系。因此，本文將從巖石的一維損傷黏彈性模型分析入手，將其推廣到復(fù)雜應(yīng)力狀態(tài)，建立三維狀態(tài)下巖石損傷黏彈性模型，并探討其參數(shù)的轉(zhuǎn)換關(guān)系。

2 巖石的一維損傷黏彈性模型

文獻[4]將 Kelvin模型和 Burgers模型進行綜合，提出可描述巖石單軸蠕變特性的nBurgers模型，如圖1所示，可較好地概括這兩種情況，公式如下：

圖1 nBurgers模型Fig.1 nBurgers model

在受到諸如含水率、溫度變化以及輻射等外界條件影響時，巖石材料的變形特性將受到影響，巖石的蠕變特性將有所改變，筆者認為，在考慮巖石內(nèi)部的這些外界因素對變形造成的影響時，可將其分為瞬間彈性損傷和長期蠕變損傷這一特性反映到巖石的蠕變模型中時，瞬間彈性損傷階段將根據(jù)加載瞬間彈性模量的劣化定義瞬間彈性損傷變量，而長期蠕變損傷階段則根據(jù)蠕變模量的劣化定義蠕變損傷變量。

（1）瞬間彈性損傷

對于損傷變量的定義，一種方法是從能量損傷中按彈性模量變化來定義的損傷變量，這種方法認為，材料性質(zhì)的“劣化”是造成材料損傷的主要因素，借鑒這一思路，在nBurgers模型中，可采用描述瞬間變形特性的彈性元件的變形參數(shù)EM的劣化來定義彈性損傷D(w)，有

式中：EM(0)為初始狀態(tài)的瞬間變形模量；EM(w)為受外界因素影響后的變形模量；w為含水率。

（2）長期蠕變損傷

與定義瞬間彈性損傷的思路一致，在定義長期蠕變損傷時將根據(jù)蠕變模量的變化來進行，認為初始狀態(tài)下巖石處于無損狀態(tài)，當(dāng)外部因素影響后，巖石蠕變模量發(fā)生衰減，由此可定義蠕變損傷變量D′(w) 為

因此，對于考慮巖石的這種損傷蠕變特性研究時，在nBurgers模型基礎(chǔ)上建立kBurgers模型，如圖2所示，一維情況時，kBurgers模型蠕變方程為

圖2 一維狀態(tài)下的kBurgers模型Fig.2 kBurgers model in the case of one-dimensional state

若將模型參數(shù)的損傷統(tǒng)一成如下表達式：

由式（6）可得

將式（6）、（7）代入式（5）得 kBurgers模型一維形式為

3 不同假定下巖石損傷蠕變模型的三維推廣

關(guān)于三維狀態(tài)下的損傷表達式，需要由一維狀態(tài)下的損傷表達式拓展而得。因此，在拓展到三維狀態(tài)下時，可做出以下設(shè)定：（1）巖石損傷的主軸與主應(yīng)力軸、應(yīng)變主軸重合；（2）初始狀態(tài)損傷為0。把蠕變本構(gòu)方程擴展到三維狀態(tài)，對材料的損傷不考慮材料各向異性的影響，D(w)和D′(w)拓展到三維狀態(tài)分別表示為 Dijkl(w)和，各損傷主軸上的損傷度仍可采用一維狀態(tài)下?lián)p傷度的求法。

3.1 常泊松比假定建模方法

朱伯芳[5]在研究復(fù)雜應(yīng)力和邊界條件下黏彈性介質(zhì)中的應(yīng)力和位移與其邊界條件的聯(lián)系時，先后提出3個與彈性體相對應(yīng)的定理，其中第1、第3定理中都是假定蠕變泊松比保持常數(shù)且等于彈性泊松比，即 μ(t′)=μ(t,t ′)=μ=常數(shù)。

由一維到三維過程，考慮常泊松比條件也存在兩種方法[6]，第1種就是從拉壓狀態(tài)的蠕變方程直接推廣到三維狀態(tài)，此時模型參數(shù)為拉壓蠕變參數(shù)，與無任何假定時一維拉壓蠕變參數(shù)是一致的；還有一種是通過黏彈性力學(xué)中的對應(yīng)性原理，通過流變本構(gòu)方程的拉普拉斯變換和逆變換得到對應(yīng)的三維蠕變方程，此時模型中的參數(shù)為剪切蠕變參數(shù)，與無任何假定時一維剪切蠕變參數(shù)是一致的[7]。

3.1.1 組合模型彈性參數(shù)選用拉壓模量

三維應(yīng)力情況下，考慮材料各向同性的標(biāo)準(zhǔn)線性固體彈性應(yīng)力-應(yīng)變關(guān)系可表示為[8]

式中：[U]為泊松比的常系數(shù)矩陣，表達式為

考慮到黏彈性模型的流變特性，由流變本構(gòu)方程利用微分算子代換[9]并求解微分方程后得到蠕變方程的基本形式：

對于nBurgers模型，則有

式中：EM、E0、ηk和ηM都是一維拉壓情況下對應(yīng)的模型參數(shù)；開關(guān)函數(shù)在一維情況下可用H(σ-σk)來表示，而在三維情況下則可用表示：

把式（6）、（7）代入式（12）可得

將式（14）改寫為張量形式：

3.1.2 組合模型彈性參數(shù)選用剪切模量

若一維情況下模型參數(shù)通過蠕變剪切試驗得到，則模型參數(shù)將采用一維剪切蠕變試驗參數(shù)，因此，可將模型參數(shù)的損傷統(tǒng)一成如下表達式：

故由式（16）可得

結(jié)合文[8]在常泊松比假定下黏彈性體的三維建模思路，若模型參數(shù)采用剪切蠕變試驗獲得（如圖3所示），則推導(dǎo)可得kBurgers模型的三維張量形式：

式中：μ*為考慮含水影響的泊松比。

因此，基于兩種假定條件下的單軸拉壓蠕變方程可分別通過式（15）和（18）進行變換得到：

圖3 常泊松比假定下的kBurgers模型Fig.3 kBurgers model in the case of static Poisson′s ratio

3.2 常體積模量假定建模方法

物體內(nèi)一點的應(yīng)力狀態(tài)，可用應(yīng)力張量表示，即等于體應(yīng)力張量和偏應(yīng)力張量之和[11]：

將模型的總應(yīng)變 εij(t)分解為由偏應(yīng)力引起的偏應(yīng)變 eij(t)和體應(yīng)力產(chǎn)生的體應(yīng)變εm(t)，即

常體積模量假定，即假設(shè)黏彈性變形只發(fā)生在剪切變形方面，體積變形是彈性的，且假定蠕變過程中體積模量保持常數(shù)且等于彈性變形時的體積模量，即 K(t′)=K(t,t′)=K 為常數(shù)。

彈性狀態(tài)下：

則根據(jù)黏彈性問題的對應(yīng)性原理，黏彈性狀態(tài)下有：

式中：J′(t)為蠕變?nèi)崃?；K為體積模量。

將式（24）代入式（22），得到三維狀態(tài)下nBurgers模型蠕變方程的張量形式：

同理，若考慮損傷效應(yīng)，參數(shù)的損傷同樣可以統(tǒng)一表示為

υ(w)和υ′(w)分別為常體積模量建模假定下的彈性損傷變量和蠕變損傷變量。

由式（27）可得：

如圖4所示，將式（27）、（28）代入式（26）得kBurgers模型的三維張量表達式：

同樣給出單軸拉壓蠕變情況下的方程，則由式（29）可得：

圖4 常體積模量假定下的kBurgers模型Fig.4 kBurgers model in the case of constant bulk modulus

3.3 各種處理方法模型參數(shù)的關(guān)系及損傷規(guī)律

比較式（8）、（19）、（20）和式（30）可知，不同假定情況下蠕變方程的形式有所不同。常規(guī)建模方式與常泊松比假定下模型的一維蠕變方程完全相同，式（8）和式（19）的方程形式及參數(shù)一致，說明若模型參數(shù)在一維和三維狀態(tài)下始終選用拉壓模量，則不考慮任何假定的常規(guī)建模方式實質(zhì)上就是常泊松比假定。

若對同一單軸壓縮條件下的蠕變試驗結(jié)果，分別用常泊松比和常體積模量假定下的模型對其進行描述時[6]，則有：

式中：ε1(t)、ε2(t)、ε3(t)分別對應(yīng)常規(guī)建模、常泊松比和常體積模量時的蠕變應(yīng)變。

于是，常體積模量和常泊松比條件下模型的參數(shù)也存在以下對應(yīng)關(guān)系：

式（32）說明，雖然選定的都是一維狀態(tài)，但在不同假定條件下組合模型的參數(shù)是不相同的。在常泊松比假定下，模型參數(shù)需在常規(guī)建模的基礎(chǔ)上乘以因子 1/[2(1+μ*)]；常體積模量假定下，不同參數(shù)的乘積因子各不相同。

若一維和三維蠕變方程使用的模型參數(shù)完全相同，則黏性參數(shù)的換算關(guān)系也就是常規(guī)方式與三維狀態(tài)中添加不同假定時黏性參數(shù)所需的換算關(guān)系，這與習(xí)慣上將一維模型的黏性參數(shù)直接運用于三維狀態(tài)，保持其值不變的做法是不同的。然而，模型參數(shù)的這些轉(zhuǎn)換做法有時也會造成一定的混亂。比較合理的解決方法，一是直接使用包含常體積模量假定的蠕變方程來處理試驗數(shù)據(jù)，從而保證一維和三維應(yīng)力狀態(tài)的一致性，避開模型參數(shù)轉(zhuǎn)換；二是仍采用常規(guī)建模方式來處理蠕變試驗數(shù)據(jù)，然后利用不同假定下模型參數(shù)的變換關(guān)系，實現(xiàn)常規(guī)建模方式的假設(shè)化。

值得思考的是，對于一維拉壓參數(shù)轉(zhuǎn)換至剪切參數(shù)時，需要獲得不同狀態(tài)下泊松比μ的大小，從而在拉壓蠕變參數(shù)轉(zhuǎn)換成剪切蠕變參數(shù)時，求解剪切蠕變參數(shù)的損傷方程還應(yīng)考慮泊松比的影響。因此，在模型參數(shù)轉(zhuǎn)換關(guān)系中涉及到泊松比μ時，模型參數(shù)的損傷表達式也會變化。例如，反映瞬間彈性變形的參數(shù)有如下計算式：

又有，

則由式（33）、（34）可得到以下關(guān)系：

式中：μ為初始狀態(tài)的泊松比；μ*受外部因素影響后的泊松比。同理可得到模型中其他參數(shù)的損傷變量關(guān)系式。

4 算例分析

文獻[4]針對 T2b2泥質(zhì)粉砂巖進行了不同加載等級對應(yīng)不同含水狀態(tài)下單軸蠕變試驗研究，探討了泥質(zhì)粉砂巖在考慮含水率w變化對巖石蠕變特性的影響，得到了瞬間彈性損傷演化方程：

以及長期蠕變損傷演化方程：

并分析得到了模型三維格式中開關(guān)函數(shù)的統(tǒng)一表達式：

式中：σ1、σ3分別為最大、最小主應(yīng)力。

因此，按照以上3種建模方式，不同含水條件下泥質(zhì)粉砂巖蠕變加載荷載為32 kN時對應(yīng)的模型參數(shù)如表1所示。

表1 3種不同建模方式下不同含水狀態(tài)對應(yīng)kBurgers模型參數(shù)反演結(jié)果Table1 Parameters of kBurgers model with different moisture contents using three different modeling approaches

常泊松比情況下，參數(shù)的損傷變量表達式與常規(guī)建模情況相比有下式：

常體積模量假定下有

因此，將試驗得到的泥質(zhì)粉砂巖不同含水率情況下泊松比μ及不同含水情況下的損傷變量D(w)和D′(w)數(shù)值代入式（39）和式（40）得到的計算結(jié)果見表2所示。由計算結(jié)果可知，蠕變方程的參數(shù)隨含水率變化的變化規(guī)律有所不同，因而損傷演化規(guī)律也不同。對常體積模量假定下模型參數(shù)隨含水損傷的演化規(guī)律進行分析，結(jié)果如圖5、6所示。

得到常體積模量假定下彈性損傷變量方程為

長期蠕變損傷變量方程為

表2 不同含水條件對應(yīng)kBurgers模型損傷變量Table2 Damage variables of kBurgers model with different moisture contents

圖5 υ(w)與w關(guān)系回歸分析結(jié)果Fig.5 Regression analysis of υ(w)and w

圖6 υ′(w)與w關(guān)系回歸分析結(jié)果Fig.6 Regression analysis of υ′(w) and w

5 結(jié) 論

（1）將考慮巖石內(nèi)部含水率變化、溫度變化以及輻射等外界因素對變形造成的影響定義為一種損傷，且包括瞬間彈性損傷和長期蠕變損傷。并根據(jù)加載瞬間彈性模量的劣化定義瞬間彈性損傷變量，根據(jù)蠕變模量的劣化定義蠕變損傷變量，由此建立考慮損傷的kBurgers模型。

（2）假定泊松比μ為常量，從一維拉壓狀態(tài)的蠕變方程入手，轉(zhuǎn)換得到了常泊松比假定的三維表達式；考慮體積模量K為常數(shù)，將黏彈性體的總應(yīng)變 εij(t )分解為由偏應(yīng)力引起的偏應(yīng)變 eij(t)和體應(yīng)力產(chǎn)生的體應(yīng)變εii，推導(dǎo)得到了常體積模量的三維表達式。通過對不同假定建模方法得到的模型參數(shù)與常規(guī)建模模型參數(shù)之間的換算關(guān)系進行研究，表明不同參數(shù)間的換算關(guān)系有所不同。

（3）結(jié)合算例，對常體積模量假定下模型參數(shù)隨含水率變化的損傷演化規(guī)律進行了分析，計算得到了常體積模量假定下模型的瞬間彈性損傷和長期蠕變損傷演化方程，且檢驗了不同假定下模型參數(shù)間轉(zhuǎn)換關(guān)系的有效性。

[1]夏才初,王曉東,許崇幫,等. 用統(tǒng)一流變力學(xué)模型理論辨識流變模型的方法和實例[J]. 巖石力學(xué)與工程學(xué)報,2008,27(8): 1594－1600.XIA Cai-chu,WNAG Xiao-dong,XU Chong-bang,et al.Method to identify rheological models by unified rheological model theory and case study[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(8): 1594－1600.

[2]徐衛(wèi)亞,楊圣奇,褚衛(wèi)江. 巖石非線性黏彈塑性流變模型（河海模型）及其應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報,2006,25(3): 433－447.XU Wei-ya,YANG Sheng-qi,CHU Wei-jiang. Nonlinear viscoelasto-plastic rheological model (Hohai model) of rock and its engineering application[J]. Chinese Journal of Rock Mechanics and Engineering,2006,25(3): 433－447.

[3]丁志坤,呂愛鐘. 巖石黏彈性非定常蠕變方程的參數(shù)辨識[J].巖土力學(xué),2004,25(增刊1): 37－40.DING Zhi-kun,Lü Ai-zhong. Parameter identification of non-stationary creep equation of rocks[J]. Rock and Soil Mechanics,2004,25(Supp.1): 37－40.

[4]黃明. 含水泥質(zhì)粉砂巖蠕變特性及其在軟巖隧道穩(wěn)定性分析中的應(yīng)用研究[D]. 重慶: 重慶大學(xué),2010.

[5]朱伯芳. 在混合邊界條件下非均質(zhì)黏彈性體應(yīng)力與位移[J]. 力學(xué)學(xué)報,1964,7(2): 162－167.ZHU Bo-fang. The stresses and deformations in the nonhomogeneous visco-elastic media undermixed boundary conditions[J]. Acta Mechanica Sinica,1964,7(2): 162－167.

[6]黃小華,馮夏庭,陳炳瑞,等. 蠕變試驗中黏彈組合模型參數(shù)確定方法的探討[J]. 巖石力學(xué)與工程學(xué)報,2007,26(6): 1226－1231.HUANG Xiao-hua,FENG Xia-ting,CHEN Bing-rui,et al.Discussion on parameters determination of viscoelastic model in creep test[J]. Chinese Journal of Rock Mechanics and Engineering,2007,26(6): 1226－1231.

[7]徐平,李云鵬,丁秀麗,等. FLAC3D黏彈性模型的二次開發(fā)及其應(yīng)用[J]. 長江科學(xué)院院報,2004,21(2): 10－13.XU Ping,LI Yun-peng,DING Xiu-li,et al. Secondary development and application of visco-elastic constitutive model in FLAC3Dsoftware[J]. Journal of Yangtze River Scientific Research Institute,2004,21(2): 10－13.

[8]黃小華,馮夏庭. 常泊松比下黏彈性體的算子代換法與黏彈對應(yīng)原理的關(guān)系[J]. 巖石力學(xué)與工程學(xué)報,2006,25(12): 2509－2514.HUANG Xiao-hua,FENG Xia-ting. Relation between operator substitution of viscoelastic body and elastoviscoelastic correspondence principle with constant Poisson ratio[J]. Chinese Journal of Rock Mechanics and Engineering,2006,25(12): 2509－2514.

[9]孫鈞. 巖土材料流變及其工程應(yīng)用[M]. 北京: 中國建筑工業(yè)出版社,1999.

[10]宋飛. 石膏角礫巖非線性流變模型研究及有限元分析[D].西安: 長安大學(xué),2006.

[11]楊桂通. 彈塑性力學(xué)引論[M]. 北京: 清華大學(xué)出版社,2004.