基于甘肅玉門花崗巖的巖石統計損傷模型研究

Research on Rock Statistical Damage Model Based on Gansu Yumen Granite

LIAO Wenjun，LIU Junxin，DIAO Hongwei，REN Tianbin，WANG Di （School of Civil Engineering and Architecture，Southwest University of Science and Technology， Mianyang 621010，Sichuan，China）

Abstract： To study the evolution behavior of rock damage under diffrent confining pressures，this study conducted triaxial tests on Yumen granite to obtain stress- strain curves and analyzed the variation patterns of its mechanical properties with confining pressure. By assuming that the micro-strength of rocks followed a Weibull probability distribution based onan existing rock damage constitutive model and considering damage thresholds，a modified strength coeffcient accounting for residual stress effects was introduced.The exponential（EXP）strength criterion for micro-element failure characterization was implemented and the statistical damage model representing the complete stress -strain relationship during triaxial compresson was established.The theoretical curves of the established constitutive model were compared with the experimental curves and D-P model curves based on granite tests.The results indicate that the EXP model can reflect the characteristics offull stress -strain curves during rock triaxial compression and performs better in post-peak residual sections than the D-P model.

Keywords：Granite；Triaxial test；Rock damage model;Exponential strength criterion；Weibull distribution

巖石作為一種天然的高強度脆性材料，在國防和民用建筑中具有重要作用。同混凝土材料相比，巖石材料的內部結構更加復雜，存在隨機分布的裂紋、孔隙等原始損傷。一般情況下巖石材料埋藏在較高地應力區，當進行地下開發或開采時產生應力場的二次分布1，其內部孔隙壓密、裂紋擴展導致力學性能不斷發生變化，而這種變化因其埋深不同而受到不同圍壓影響。因此，研究圍壓對巖石損傷變化的影響，建立和完善巖石本構模型，對災害防控和工程建設具有重要意義。

Krajcinovic等2]假定巖石微元強度服從Weibull概率分布，并首次將統計損傷理論應用于構建巖石統計損傷本構模型。基于這種關系，曹文貴等[3-4]、劉樹新等[5]、李樹春等[6采用 Drucker -Prager及Mohr-Coulomb強度準則來描述巖石微元強度的變化，建立了反映巖石破壞的統計損傷本構模型。徐衛亞等建立了能夠反映軟化特性的損傷模型，但描述殘余階段特征的效果不好；薛云亮等[8-9]在此基礎上建立了能夠反映巖石殘余強度特性的統計本構損傷模型。眾多學者深入研究發現：有較明顯的保守性的Drucker-Prager準則使得在一定水平上所確定的微元強度存在偏差，而以Mohr-Coulomb強度準則構建的本構模型則因其未能充分考慮巖石受中主應力影響導致在預測高、低應力區破壞強度時亦有一定局限[10]。付義勝[11]通過比較5種常規三軸強度準則發現， You^[12] 提出的包含3個材料參數的指數強度準則（即EXP強度準則）綜合考慮了多種因素的影響，能夠較好描述巖石加載變化的全過程。

本研究基于甘肅玉門花崗巖三軸試驗結果，在巖石微元強度服從Weibull分布假定的基礎上，引入殘余強度修正系數，并結合EXP強度準則，建立能夠預測圍壓作用下巖石損傷演化全過程的統計本構模型，并通過試驗曲線與理論曲線的對比驗證了模型的合理性和適用性，為地下與巖體工程建設提供參考。

1 中商 試驗與分析

1.1 試驗方案

本次試驗所用材料為甘肅玉門組粗粒花崗巖，主要礦物成分為半自粒形鉀長石（ 47% ）、石英（ 27% ）斜長石（ 23% ）等。根據《工程巖體試驗方法標準》[13]，將所取得的樣品巖芯加工為高度100mm 、直徑 50mm 的標準試驗尺寸試樣，如圖1所示。試驗設備為YSY-2500型三軸試驗系統，該系統的最大豎向荷載為 2500kN ，最高圍壓為120MPa ，試驗控制方式為位移控制，設置試驗的加載速率為 0.06mm/min 。

為了獲取能夠反映圍壓對花崗巖損傷變化影響的應力-應變曲線，本次試驗設計0，15，30，45，60，75MPa6 個圍壓，每個圍壓下設置5組平行試驗，試件編號為 GSYM-x（GSYM 指甘肅玉門， Ψ_，x 為1，230）。

1.2 試驗結果與分析

1.2.1 應力-應變曲線特征

基于試驗結果，繪制各圍壓下的應力-應變曲線并分析其特征變化。

圖2展示了不同圍壓下的應力-應變曲線。從圖2可以看出曲線均呈現出脆性巖石所具有的階段性，即在加載前期試樣內部孔隙逐漸壓密，試樣進入線彈性變形階段；當加載應力持續增長至超過材料的屈服極限，進入塑性變形階段，曲線呈現非線性，內部出現損傷并逐漸積累至失穩破壞。殘余強度與圍壓在試驗范圍呈正相關變化，根據巖石臨界狀態理論14分析，當圍壓由低到高逐漸增加，巖石由脆性向延性狀態轉化，導致殘余強度總體呈非線性增加。

1.2.2 花崗巖物理性質分析

為進一步探討圍壓水平對花崗巖應力-應變曲線特征影響，描述偏應力峰值強度隨圍壓增長的變化規律，通過圖2曲線計算并統計三軸試驗中花崗巖的物理參量，結果如表1所示。基于其結果擬合峰值強度與圍壓關系曲線，計算統計花崗巖黏聚力和摩擦角，如圖3所示。從圖3（a）可以看出，在低圍壓區域花崗巖峰值離散性較小，而隨著圍壓的增加各平行試樣峰值差異開始增加。結合表1結果分析，發現在圍壓 15～45MPa 區間，峰值應力與彈性模量隨圍壓增長的幅值穩定， 45MPa 之后峰值增長幅值與泊松比開始下降，即表明圍壓 45MPa 時巖樣由脆性轉為延性（抗剪強度降低）。圖4給出了不同圍壓下玉門花崗巖的破壞形式，可以清楚看到巖樣的宏觀破壞形式隨圍壓增長而發生改變，由低圍壓至高圍壓，從剪切型破壞逐步向剪切劈裂混合型破壞轉變，而其破壞形態也由低圍壓時的單一斷面轉變為高圍壓狀態下的多斷面破壞。分析認為當處于高圍壓狀態時花崗巖進入延性，抗剪強度與圍壓變化關系轉為非線性，抗剪強度增加不足以支撐相應軸向應力增加，泊松比降低，巖樣斷面最終出現混合型破壞。

2 統計損傷演化

2.1 損傷演化

在圍壓與軸壓作用下，巖體因內部微裂隙、裂紋的發展導致破壞[15]，假定其微元強度變化遵循分布統計規律，采用Weibull分布[1來描述巖石微元體強度分布概率密度。考慮到巖石材料存在損傷閾值的情況[17]，將材料抽象為未受損傷和受損傷兩部分，可得到巖石的損傷函數表達式：

式中： m，F₀ 為Weibull函數參數; F 為微元體強度；

D 為統計損傷變量值，且 0?D?1 。

經過塑性段，處于加載晚期的巖石受到圍壓和摩擦力的影響，仍然存有一定殘余強度，巖石微元破壞后繼續傳遞少量壓力和剪應力[18]，因此本文引入修正系數 γ 來描述這一階段變化。根據應變等價理論和有效應力概念[19]，得到引人殘余修正參數的損傷本構關系：

σ_i^*（1-γD）=σ_i，i=1，2，3

式中： σ_i^* 為有效應力，即物體受到的實際應力值； σ_i 為名義應力，即忽略物體不連續性時有效橫截面的應力； γ 為殘余修正系數。

在巖石破壞進入殘余階段時，可將損傷值 D 近似為0，此時式（2）可轉化為：

σ₁=Eε_r（1-γ）+2μσ₃

變換式（3）即可得到修正系數：

γ=（Eε_r-σ_r+2μσ₃）/（Eε_r）

式中： σ_r 為殘余應力； ε_r 為殘余應變。

在常規三軸試驗中主應力與圍壓關系為 σ₁gt; σ₂，σ₂=σ₃ ，那么假定巖石微元體破壞前服從廣義

的虎克定律，即可得到有效應力-應變關系式：

式中： ε₁ 為豎向主應變； E 為彈性模量； ?_：μμ 為泊松比。

聯立式（1）式（2）式（5）得到基于Weibull分布的巖石損傷軟化本構關系式：

2.2 巖石微元體強度

微元體是指細觀統計損傷的最小單位，該單位既可作為損傷力學中的質點，又能很好體現應力-應變的相關作用特性。巖石材料是一種具有黏結和摩擦特性的脆性材料，因此本文引入能夠較好描述這種特性的文獻[11]提出的EXP強度準則來構建統計損傷模型：

即：

式中： σ_c 為峰值應力； σ₀ 為單軸壓縮強度； K 為圍壓等于 0MPa 時強度隨圍壓的增長速率。

這里將EXP強度準則引入并采用有效應力來描述巖石微元體強度 F 的演化：

F=f（σ）=σ₁^*-σ₃^*-σ_c^*+

將式（5）代入式（9）得到可通過常規三軸試驗計算的基于EXP強度準則的微元強度函數，筆者采用該函數來模擬巖石材料已經發生損傷階段的微元體強度演化：

F=Eε₁+2μσ₃-σ₃-σ_c+

2.3 巖石損傷本構模型

在加載過程中，巖石由內部裂紋產生并發展貫通導致失效破壞。當應力值處于較低區域時，巖石內部裂隙未發生擴展損傷現象，因此，巖石材料的損傷存在一個閾值點，當加載應力還未增長至損傷閾值時，可認為材料沒有產生損傷。在此階段，巖石的微元損傷值 D 近似為零，因此可以采用彈性模量和應變來直接描述這一階段的強度和變形演化，如式（11）所示：

σ₁=Eε₁+2μσ₃，Flt;0

當加載應力持續增加，巖石內部產生裂紋，開始產生損傷。對于發生損傷的階段，這里將式（10）代入式（1），得到巖石統計損傷本構方程（損傷閾值點之后）：

結合未發生損傷段本構方程（11）和已發生損傷段本構方程（12），得到基于EXP強度準則的巖石統計損傷軟化本構方程：

本構方程中的參數 m，F₀ 均可由三軸試驗中所獲得的不同圍壓下巖石的峰值應力和應變確定，在相應狀態下是具有普遍意義的模型參數。結合文獻[20]，根據巖石全應力-應變曲線幾何條件求解得到參數的確定方法：

式中 F_c 為峰值狀態對應巖石微元強度。

3模型適用性驗證

基于花崗巖三軸試驗結果計算不同圍壓下的EXP模型參數，并將試驗曲線與理論模型曲線進行對比，對本文建立的損傷模型的準確性進行驗證。為了更好評價該模型的適用性，本文引用文獻[21]的Drucker-Prager損傷模型（后文簡稱D-P模型）作對比。試驗的相關數據如表2所示，EXP模型中參數 K 取值如圖3所示（經擬合， K=6.84 ）。

根據表2中整理的試驗數據及2.3節參數的計算方法得到EXP模型參數，同時計算確定相應的D-P損傷模型參數，結果如表3所示。在此基礎上依據式（1）即可以獲得花崗巖整個破壞過程中的損傷演化過程，如圖5所示。圖5中損傷變量 D 的演化曲線均呈現“S”型，且表現出隨圍壓增加損傷延緩的規律。特別地，圍壓 60MPa 時曲線與 45MPa 曲線交叉，認為在 45MPa 時玉門花崗巖進入脆延轉換臨界狀態，圍壓延緩損傷作用降低。

根據表3中的本構參數，繪制EXP模型理論曲線，并將試驗曲線同理論曲線進行對比，結果如圖6所示。通過對比分析，可以發現EXP模型理論曲線和試驗曲線的吻合水平較高，且在反映峰后軟化階段有良好效果，反映了圍壓對花崗巖三軸試驗應力－應變關系的影響。與文獻[21]中D-P模型曲線相比，在峰前二者重合度較高，而峰后殘余軟化段

EXP模型更加穩定，與試驗曲線偏差更小，驗證了該模型的適用性及準確性。需指出的是，EXP模型在結果上也存在一定誤差，可能是由于試驗樣本相對較少及參數計算對精度產生了影響。

4結論

基于玉門花崗巖常規三軸試驗結果，在前人巖石損傷理論研究的基礎上，基于巖石微元強度服從Weibull分布的假定，考慮到巖石材料存在損傷閾值及加載后期殘余強度的影響，引入指數強度準則，建立統計損傷本構模型，得到以下結論：（1）玉門花崗巖在圍壓超過 45MPa 后進入脆延轉換的臨界狀態，圍壓效應降低（損傷延緩效果下降），最終呈現出多斷面的破壞模式。（2）假定巖石微元強度服從

Weibull分布，考慮到材料存在損傷閾值及加載后期殘余強度，引入殘余強度修正系數，結合擁有3個材料參數的指數強度準則作為巖石微元強度演化的依據，建立了可反映巖石加載過程全應力曲線的統計損傷本構模型（EXP模型）。EXP模型曲線和試驗曲線具有較高水平的一致性，能夠較好地反映不同圍壓花崗巖三軸試驗加載全過程應力-應變關系，且在峰后軟化階段優于現有的D-P模型，說明EXP巖石損傷模型具有較高的準確性和較好的適用性。

參考文獻

[1]SHOUYD，ZHOUXP，QIANQH.A critical conditionofthezonal disintegration in deeprock masses：Strainen-ergy density approach[J]. Theoretical and Applied Frac-tureMechanics，2018，97：322-332.

[2]KRAJCINOVIC D，SILVA M A G. Statistical aspects ofthe continuous damage theory[J]. International Journal ofSolids and Structures，1982，18（7） ： 551-562.

［3］曹文貴，方祖烈，唐學軍．巖石損傷軟化統計本構模型之研究[J]．巖石力學與工程學報，1998，17（6）：628-633.

[4] 曹文貴，趙明華，唐學軍．巖石破裂過程的統計損傷模擬研究[J]．巖土工程學報，2003，25（2）：184-187.

[5] 劉樹新，劉長武，韓小剛，等．基于損傷多重分形特征的巖石強度Weibull參數研究［J]．巖土工程學報，2011，33（11）：1786-1791.

［6］李樹春，許江，李克鋼，等．基于Weibull分布的巖石損傷本構模型研究［J]．湖南科技大學學報（自然科學版），2007，22（4）：65-68.

[7] 徐衛亞，韋立德．巖石損傷統計本構模型的研究［J].巖石力學與工程學報，2002，21（6）：787-791.

［8］薛云亮，李庶林，徐宏斌，等．考慮殘余強度的巖石損傷統計本構模型[C]/第二屆中國水利水電巖土力學與工程學術討論會論文集（二），2008：3-8.

[9]曹文貴，趙衡，李翔，等．基于殘余強度變形階段特征的巖石變形全過程統計損傷模擬方法[J]．土木工程學報，2012，45（6）：139-145.

［10］王軍保，劉新榮，李鵬.巖石損傷軟化統計本構模型[J].蘭州大學學報（自然科學版），2011，47（3）：24-28

[11］付義勝．五種常規三軸強度準則的比較［J]．采礦技術，2012，12（3）：74-78.

[12]YOUMQ.True-triaxial strengthcriteriaforrock[J].In-ternational Journal of Rock Mechanicsand Mining Sci-ences，2009，46（1）：115-127.

[13］中華人民共和國住房和城鄉建設部.工程巖體試驗方法標準：GB/T50266—2013［S]．北京：中國計劃出版社，2013.

［14］張晶，任鳳玉，何榮興，等．包含巖石臨界狀態的Mo-hr-Wedge準則適用性研究［J]．東北大學學報（自然科學版），2021，42（4）：567-575.

［15］王東，張婧，陳強，等．基于3種破壞類型的巖石損傷軟化統計模型［J]．巖石力學與工程學報，2015，34（S2）：3759-3765.

［16］楊明輝，趙明華，曹文貴.巖石損傷軟化統計本構模型參數的確定方法[J].水利學報，2005，36（3）：345-349.

［17］張超，曹文貴，王江營．考慮損傷閾值影響的巖石脆-延性轉化統計損傷本構模型研究［J]．水文地質工程地質，2013，40（5）：45-50，57.

[18］楊圣奇，徐衛亞，韋立德，等．單軸壓縮下巖石損傷統計本構模型與試驗研究［J]．河海大學學報（自然科學版），2004，32（2）：200-203.

[19] JEAN L. A continuous damage mechanics model for duc-tilefracture[J].Journal of Engineering Materials andTechnology，1985，107（1）：83-89.

[20］曹瑞瑯，賀少輝，韋京，等．基于殘余強度修正的巖石損傷軟化統計本構模型研究［J]．巖土力學，2013，34（6）： 1652-1660，1667.

[21］肖南，劉啟明，邱祥，等．考慮殘余強度階段裂隙演化影響的巖石統計損傷本構模型［J/OL]．西安建筑科技大學學報（自然科學版），2024：1-9（2024-06-20）.https：//kns.cnki.net/kcms/detail/61.1295. TU.20240619.1311.004.html.