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

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

1 引言

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

2 損傷本構關系的建立

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

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

3 巖石微元強度和損傷變量的計算

3.1 基于D_P準則的巖石微元強度的求取

基于D_P準則表示的巖石微元強度為：

常規(guī)三軸試驗中，可以測得巖石名義應力σ1，σ1=σ3以及應變ε1，由（1）式和虎克定律可以求出對應的有效應力和應變：

把式（3）～（8）代入式（2），得巖石微元強度表達式為：

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

巖石是一種非均質材料，內含大量隨機分布的裂隙、空洞、界面等缺陷，因此，巖石微元在壓力作用下的破壞也應該是隨機的。假定巖石微元破壞的概率與應力、應變有關，且微元強度服從Weibull分布，則概率密度函數為：

定義巖石損傷變量D為：

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

圖1 巖石全應力-應變曲線

可知在OA和AB階段內，巖石通常呈彈性性質，有時在OA段也有殘余變形，AB段可能出現(xiàn)細微的開裂，但是巖石的結構和性質這兩階段中并無大的改變，所以我們把屈服強度B點對應的應變εb定義為應變閥值[8]，則損傷變量可以表示為：

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

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

4 參數m,F0的計算與修正

4.1 參數m,F0的計算

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

①通過對本構方程進行2次對數變換和線性擬合求m，F(xiàn)0。過程如下：

對式（13）進行對數變換得：

再對式（15）進行對數變換得：

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

擬合曲線可以表示為：

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

②通過求應力對應變的偏微分，再代入應力-應變曲線峰值點求m，F(xiàn)0。過程如下：

式中：Fc為巖石在峰值點的微元強度，通過把

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

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

其中：

然后把巖石應力-應變試驗的相關數據代入即可得到具體值。

4.2 參數m,F0的修正

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

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

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

翻轉課堂課前學習的效果無法得到保障，語文教師可通過軟件的管理系統(tǒng)實施的查看學生課前的學習動態(tài)，對哪些學生進行課前預習，哪些學生沒有做到課前的統(tǒng)計，對自主學習的學生進行鼓勵和表揚，并對平時成績進行統(tǒng)計，作為課程改革考試的一項重要的指標，由此可以促進督促學生課前的預習，在一定程度上保障了教學的質量，保障了學生學習的效果。[4]

圖2 參數m隨圍壓σ3變化的關系

圖3 參數F0隨圍壓σ3變化的關系

擬合后的關系式分別為：

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

5 損傷本構方程的應用舉例探討

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

5.1 張拉破壞

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

式中：ε3為巖石最小主應變，壓為正，拉為負。

5.2 剪切破壞

當巖石剪切應變γ超過臨界值即可認為發(fā)生剪切破壞。由八面體剪應變理論定義臨界值γ0：

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

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

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

6 結論

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

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

③對Weibull分布模型參數m，F(xiàn)0的求取提出了兩種方法，即傳統(tǒng)的二次取對數法和應力對應變取偏微分法，并對參數進行了擬合修正。

④舉例說明了損傷本構方程的應用，同時也提出了不足之處，以及討論了當前工作的難點所在，指明了下一步研究的方向。