巖石損傷本構(gòu)方程的研究和應(yīng)用

余秀平，劉保川，吳軍國(guó)，胥克明，高飛，楊旭

（中鐵四局集團(tuán)第二工程有限公司，江蘇 蘇州 215131）

1 引言

巖石本構(gòu)關(guān)系的研究是巖石力學(xué)工作者關(guān)注的重點(diǎn)。實(shí)際工程活動(dòng)中，巖體基本上是經(jīng)過(guò)地質(zhì)活動(dòng)后的損傷巖體，因此，對(duì)巖石損傷本構(gòu)方程的研究顯得尤為重要。由曹文貴教授，徐衛(wèi)亞教授[1,2]引進(jìn)的巖石微元強(qiáng)度服從Weibull分布的概念，建立了巖石損傷統(tǒng)計(jì)本構(gòu)模型，發(fā)展了損傷本構(gòu)方程的研究，但是其模型的建立是基于莫一條試驗(yàn)曲線的，不具有普遍性。本文在前人研究的基礎(chǔ)上，基于D_P準(zhǔn)則和材料存在應(yīng)變閥值[3]的事實(shí)，對(duì)損傷本構(gòu)方程進(jìn)行了詳細(xì)的推導(dǎo)，并對(duì)相關(guān)參數(shù)進(jìn)行了擬合處理，提高了普遍適用性，以期對(duì)實(shí)際巖石工程實(shí)踐有更好的指導(dǎo)作用。

2 損傷本構(gòu)關(guān)系的建立

Lemaitre應(yīng)變等價(jià)性原理[4,5]認(rèn)為在有效應(yīng)力作用下產(chǎn)生的應(yīng)變與同種材料無(wú)損時(shí)發(fā)生的應(yīng)變等價(jià)。傳統(tǒng)的損傷理論認(rèn)為力不能通過(guò)損傷傳播，然而在巖石受壓過(guò)程中，巖石微元破壞后還可以傳遞部分壓應(yīng)力和剪應(yīng)力。根據(jù)文獻(xiàn)[2,6]，假定破壞面積對(duì)壓應(yīng)力和剪應(yīng)力的傳遞效果一樣，設(shè)傳壓系數(shù)為cn。因?yàn)閹r石破壞微元在空間3個(gè)主方向的投影面積和總面積比率都一樣大，所以各個(gè)方向的損傷都用D來(lái)表示，則巖石損傷本構(gòu)關(guān)系為：

式中：{ε}為應(yīng)變矢量，{σ*}為有效應(yīng)力矢量，{σ}為名義應(yīng)力矢量，[C]為彈性柔度矩陣，D為損傷變量。

3 巖石微元強(qiáng)度和損傷變量的計(jì)算

3.1 基于D_P準(zhǔn)則的巖石微元強(qiáng)度的求取

基于D_P準(zhǔn)則表示的巖石微元強(qiáng)度為：

常規(guī)三軸試驗(yàn)中，可以測(cè)得巖石名義應(yīng)力σ1，σ1=σ3以及應(yīng)變?chǔ)?，由（1）式和虎克定律可以求出對(duì)應(yīng)的有效應(yīng)力和應(yīng)變：

把式（3）～（8）代入式（2），得巖石微元強(qiáng)度表達(dá)式為：

3.2 基于Weibull分布的巖石損傷變量的確定

巖石是一種非均質(zhì)材料，內(nèi)含大量隨機(jī)分布的裂隙、空洞、界面等缺陷，因此，巖石微元在壓力作用下的破壞也應(yīng)該是隨機(jī)的。假定巖石微元破壞的概率與應(yīng)力、應(yīng)變有關(guān)，且微元強(qiáng)度服從Weibull分布，則概率密度函數(shù)為：

定義巖石損傷變量D為：

根據(jù)文獻(xiàn)[7]的研究：任何形式的材料內(nèi)部缺陷或微損傷，其演化過(guò)程同時(shí)依賴于應(yīng)變和應(yīng)變率，即D=D(ε,)。對(duì)于恒應(yīng)變率過(guò)程，損傷演化存在某個(gè)應(yīng)變閾值。由剛性壓力機(jī)做巖石力學(xué)試驗(yàn)得到的全應(yīng)力-應(yīng)變曲線如圖1。

圖1 巖石全應(yīng)力-應(yīng)變曲線

可知在OA和AB階段內(nèi)，巖石通常呈彈性性質(zhì)，有時(shí)在OA段也有殘余變形，AB段可能出現(xiàn)細(xì)微的開(kāi)裂，但是巖石的結(jié)構(gòu)和性質(zhì)這兩階段中并無(wú)大的改變，所以我們把屈服強(qiáng)度B點(diǎn)對(duì)應(yīng)的應(yīng)變?chǔ)舃定義為應(yīng)變閥值[8]，則損傷變量可以表示為：

聯(lián)立式（6），（7），（8），（11）可得巖石三軸全應(yīng)力-應(yīng)變曲線為：

再聯(lián)立式（12），（13）可得巖石損傷統(tǒng)計(jì)本構(gòu)方程為：

4 參數(shù)m,F0的計(jì)算與修正

4.1 參數(shù)m,F0的計(jì)算

參數(shù)m，F(xiàn)0的求取一般有兩種方法。

①通過(guò)對(duì)本構(gòu)方程進(jìn)行2次對(duì)數(shù)變換和線性擬合求m，F(xiàn)0。過(guò)程如下：

對(duì)式（13）進(jìn)行對(duì)數(shù)變換得：

再對(duì)式（15）進(jìn)行對(duì)數(shù)變換得：

令b=ln(a)，x=ln(k)，F(xiàn)0則可以表示為：

擬合曲線可以表示為：

通過(guò)線性回歸得到m，b的具體值，再代入到式（17）即可求得F0。

②通過(guò)求應(yīng)力對(duì)應(yīng)變的偏微分，再代入應(yīng)力-應(yīng)變曲線峰值點(diǎn)求m，F(xiàn)0。過(guò)程如下：

式中：Fc為巖石在峰值點(diǎn)的微元強(qiáng)度，通過(guò)把

σc，εc代入式（9）可以計(jì)算出具體值。再把σc，εc帶入式（13）得：

聯(lián)立式（19），（20），整理后可得：

其中：

然后把巖石應(yīng)力-應(yīng)變?cè)囼?yàn)的相關(guān)數(shù)據(jù)代入即可得到具體值。

4.2 參數(shù)m,F0的修正

由式（14）表示的損傷本構(gòu)方程是建立在單條試驗(yàn)曲線上的，因此在反映不同應(yīng)力狀態(tài)下的本構(gòu)關(guān)系時(shí)仍然存在誤差[9]。為了反映多條試驗(yàn)曲線的綜合情況，通過(guò)建立m，F(xiàn)0與圍壓σ3的關(guān)系，對(duì)其數(shù)據(jù)進(jìn)行擬合處理，得到他們之間的關(guān)系式，達(dá)到參數(shù)修正的目的。引用文[10]中的數(shù)據(jù)，見(jiàn)表1。

不同圍壓下?lián)p傷統(tǒng)計(jì)擬合參數(shù) 表1

對(duì)m-σ3，F(xiàn)0-σ3的散點(diǎn)分布分別進(jìn)行二次多項(xiàng)式進(jìn)行擬合，結(jié)果如圖2，圖3所示。

翻轉(zhuǎn)課堂課前學(xué)習(xí)的效果無(wú)法得到保障，語(yǔ)文教師可通過(guò)軟件的管理系統(tǒng)實(shí)施的查看學(xué)生課前的學(xué)習(xí)動(dòng)態(tài)，對(duì)哪些學(xué)生進(jìn)行課前預(yù)習(xí)，哪些學(xué)生沒(méi)有做到課前的統(tǒng)計(jì)，對(duì)自主學(xué)習(xí)的學(xué)生進(jìn)行鼓勵(lì)和表?yè)P(yáng)，并對(duì)平時(shí)成績(jī)進(jìn)行統(tǒng)計(jì)，作為課程改革考試的一項(xiàng)重要的指標(biāo)，由此可以促進(jìn)督促學(xué)生課前的預(yù)習(xí)，在一定程度上保障了教學(xué)的質(zhì)量，保障了學(xué)生學(xué)習(xí)的效果。[4]

圖2 參數(shù)m隨圍壓σ3變化的關(guān)系

圖3 參數(shù)F0隨圍壓σ3變化的關(guān)系

擬合后的關(guān)系式分別為：

相關(guān)擬合系數(shù)分別為0.9145，0.9758。把式（26），（27）代入式（14），再選取合適的傳壓系數(shù)，如令cn=0.80，即可得修正后完整的損傷本構(gòu)方程：

5 損傷本構(gòu)方程的應(yīng)用舉例探討

巖石微元體拉，剪作用是同時(shí)存在的，因此哪一個(gè)先達(dá)到臨界狀態(tài)，便以哪種形態(tài)發(fā)生破壞，所以用單一的破壞準(zhǔn)則來(lái)研究巖石破壞是不全面的。

5.1 張拉破壞

最大拉應(yīng)變理論認(rèn)為巖石破壞是由拉應(yīng)變引起的，而在巖石三軸壓縮試驗(yàn)中，最大拉應(yīng)變發(fā)生在最小主應(yīng)力方向上，因此定義巖石張拉破壞時(shí)的強(qiáng)度指標(biāo)T為：

式中：ε3為巖石最小主應(yīng)變，壓為正，拉為負(fù)。

5.2 剪切破壞

當(dāng)巖石剪切應(yīng)變?chǔ)贸^(guò)臨界值即可認(rèn)為發(fā)生剪切破壞。由八面體剪應(yīng)變理論定義臨界值γ0：

常規(guī)三軸試驗(yàn)中，σ2=σ3，所以，ε2=ε3，可化簡(jiǎn)為：

式中：ε1和ε3是巖石的最大主應(yīng)變和最小主應(yīng)變，壓為正，拉為負(fù)。把式（28）表示的曲線中某一點(diǎn)所對(duì)應(yīng)的相關(guān)數(shù)據(jù)分別代入到式（33），（35），然后再比較T和γ0，即可知道巖石以何種形式破壞。

巖石發(fā)生張拉破壞還是剪切破壞已經(jīng)有了相應(yīng)的判別方式，而通過(guò)試驗(yàn)擬合T和γ0與圍壓的關(guān)系，以期建立比較精確的巖石強(qiáng)度準(zhǔn)則是下一步研究的重點(diǎn)。

6 結(jié)論

①為了反映巖石破壞后能夠繼續(xù)傳遞力的能力，在傳統(tǒng)有效應(yīng)力計(jì)算方程中加入傳遞系數(shù)cn，更加符合實(shí)際情況。

②針對(duì)材料在恒應(yīng)變率過(guò)程中損傷演化存在應(yīng)變閾值，對(duì)巖石損傷本構(gòu)方程進(jìn)行了分段處理，與傳統(tǒng)方法相比提高了對(duì)巖石全應(yīng)力-應(yīng)變曲線的擬合度。

③對(duì)Weibull分布模型參數(shù)m，F(xiàn)0的求取提出了兩種方法，即傳統(tǒng)的二次取對(duì)數(shù)法和應(yīng)力對(duì)應(yīng)變?nèi)∑⒎址?，并?duì)參數(shù)進(jìn)行了擬合修正。

④舉例說(shuō)明了損傷本構(gòu)方程的應(yīng)用，同時(shí)也提出了不足之處，以及討論了當(dāng)前工作的難點(diǎn)所在，指明了下一步研究的方向。