李 博, 吳潤(rùn)江, 高為超

(西北核技術(shù)研究所,陜西 西安 710024)

?

巖石損傷軟化的修正本構(gòu)模型

李 博, 吳潤(rùn)江, 高為超

(西北核技術(shù)研究所,陜西 西安 710024)

利用服從Weibull分布的巖石微元強(qiáng)度表示方法,基于Drucker-Prager破壞準(zhǔn)則,通過(guò)引入損傷變量的修正系數(shù),利用應(yīng)力-應(yīng)變關(guān)系曲線峰值點(diǎn)處的應(yīng)力、應(yīng)變確定Weibull分布參數(shù)的關(guān)系式,建立不同圍壓下的損傷軟化本構(gòu)模型。該模型參數(shù)較少且易于確定,其參數(shù)的確定方法揭示了模型參數(shù)的物理意義。與前人的研究成果進(jìn)行對(duì)比分析,結(jié)果表明該模型比未修正前有著明顯的優(yōu)越性,且與實(shí)測(cè)結(jié)果吻合較好,分析結(jié)果顯示出該模型的合理性。

巖石; 損傷; 軟化; 修正; 本構(gòu)模型

0 引言

巖石損傷是指在加載條件下,由于巖石內(nèi)部存在原始微觀缺陷,其承載后又產(chǎn)生新的裂紋,并隨荷載增加呈逐漸擴(kuò)展趨勢(shì),巖石材料體積單元發(fā)生破壞、力學(xué)性能劣化的現(xiàn)象。巖石損傷主要指巖石內(nèi)部微裂紋的形成、擴(kuò)展和聚合。巖土材料的強(qiáng)度與應(yīng)變關(guān)系可以通過(guò)損傷力學(xué)來(lái)解決,而巖石材料損傷本構(gòu)關(guān)系可通過(guò)兩種方式研究:一是以實(shí)驗(yàn)為基礎(chǔ),假設(shè)損傷變量與應(yīng)力、應(yīng)變服從某種關(guān)系,再以假設(shè)模型模擬試驗(yàn)所得應(yīng)力-應(yīng)變關(guān)系,建立損傷本構(gòu)模型;二是根據(jù)巖石微元強(qiáng)度符合概率分布,先建立損傷變量與應(yīng)力、應(yīng)變之間的關(guān)系,進(jìn)而建立起基于概率分布的本構(gòu)模型[1-2]。

1 損傷變量與損傷本構(gòu)關(guān)系

三軸試驗(yàn)通常為假三軸,取σ2=σ3,利用J.Lemaitre應(yīng)變等效性假說(shuō)[3]建立巖石損傷的本構(gòu)關(guān)系[4]:

[σ*]=[σ]/(1-D)=[C][ε]/(1-D)

(1)

式中:[σ*]為有效應(yīng)力矩陣;[C]為巖石材料彈性矩陣;[σ]為名義應(yīng)力矩陣;[ε]為應(yīng)變矩陣;D為巖石損傷變量。

假定巖石微元強(qiáng)度F服從Weibull分布[4-6],其概率密度函數(shù)可以表示為:

(2)

式中:F為微元破壞Weibull分布的分布變量;m、F0為Weibull分布參數(shù)。

巖石損傷變量D可定義為損傷面積與無(wú)損時(shí)材料全面積之比,考慮巖石材料損傷的隨機(jī)性,巖石損傷變量D可表示為:

(3)

曹文貴等[2-7]率先提出基于巖石破壞準(zhǔn)則的新型巖石微元強(qiáng)度度量方法,此方法較為合理,能反映應(yīng)力狀態(tài)對(duì)巖石微元強(qiáng)度的影響。由于Drucker-Prager破壞準(zhǔn)則參數(shù)形式簡(jiǎn)單且適用于巖石介質(zhì)而得到廣泛應(yīng)用,因此本文巖石微元強(qiáng)度采用該準(zhǔn)則,即

(4)

(5)

巖石微元體破壞前符合虎克定律,即:

(6)

利用式(1),并考慮式(6),由式(5)可得:

(7)

將式(7)代入到式(4)可得:

(8)

2 巖石損傷本構(gòu)模型的建立及其參數(shù)的確定

2.1 巖石損傷本構(gòu)模型的建立

謝和平[8]認(rèn)為三維條件下?lián)p傷變量為一個(gè)代表性體積元素內(nèi)的損傷等效面積與該截面總面積的比值,如果巖石材料各向同性,比值與截面方位無(wú)關(guān),各應(yīng)力分量的損傷程度相同,為同性損傷。但這只有理想巖石材料才滿足,因此假設(shè)損傷情況滿足Weibull分布是不合適的,本文嘗試引入一個(gè)修正系數(shù)對(duì)1-D進(jìn)行修正,使其能夠反映巖石破裂后的殘余強(qiáng)度特征。

利用式(1)及式(6)可得:

σ1=qEε(1-D)+2μσ3=

(9)式中:q為修正系數(shù)。

2.2 巖石損傷本構(gòu)模型參數(shù)的確定方法

分別對(duì)式(8)、式(9)求導(dǎo)可得:

(10)

(11)

(12)

將式(12)代入到式(11)可得:

(13)

由式(13),并取峰值強(qiáng)度處F=Fc可得:

(14)

在峰值點(diǎn)處,將式(14)代入式(9)可得:

(15)

由式(15)可得:

(16)

(17)

3 修正的巖石損傷本構(gòu)模型的驗(yàn)證

為驗(yàn)證修正模型,引用文獻(xiàn)[9]中的數(shù)據(jù)資料,巖石為大冶鐵礦大理巖,可得其彈性模量E=25.34 GPa,泊松比μ=0.25,內(nèi)摩擦角φ=25.73°,黏聚力c=23.9 MPa,單軸抗壓強(qiáng)度σc=70.0 MPa。

將文獻(xiàn)[8]中實(shí)驗(yàn)數(shù)據(jù)匯總,利用Fc、F0、m等的公式,可求出Fc和m值,F0用m表示。

在不同圍壓條件下,利用式(9)可求出σ1-σ3隨應(yīng)變?chǔ)?變化的關(guān)系式。試驗(yàn)曲線與理論曲線的比較圖見(jiàn)圖1。

從圖1中可以看出不同的q值擬合所得的曲線圖不同,但本文通過(guò)修正得出的擬合曲線與試驗(yàn)曲線比較接近,通過(guò)取不同q值的方法可以得到與試驗(yàn)曲線最為接近的擬合曲線。不同圍壓下的對(duì)比分析結(jié)果表明該修正方法相比q=1時(shí)的方法具有一定的優(yōu)越性,但與實(shí)際試驗(yàn)曲線仍有一定程度上的出入,需要進(jìn)一步地研究改進(jìn)。

圖1 試驗(yàn)曲線與理論曲線的比較Fig.1 Comparison between experimental and theoretical curves

4 結(jié)論

基于巖石破壞準(zhǔn)則的巖石微元強(qiáng)度度量方法,從巖石微元強(qiáng)度服從Weibull分布的角度出發(fā),通過(guò)對(duì)損傷變量進(jìn)行修正,建立了不同圍壓下巖石損傷軟化本構(gòu)模型,并得出如下結(jié)論:

(1) 通過(guò)對(duì)損傷變量進(jìn)行修正,建立起損傷軟化本構(gòu)模型,該模型能夠反映巖石在復(fù)雜應(yīng)力狀態(tài)下變性破壞的全過(guò)程,尤其可反映巖石的軟化特性;

(2) 該模型可通過(guò)峰值點(diǎn)應(yīng)力、應(yīng)變確定模型參數(shù),使得模型參數(shù)的物理意義更加明確,增強(qiáng)參數(shù)計(jì)算的準(zhǔn)確性,應(yīng)用比較簡(jiǎn)便;

(3) 與未修正前模型的對(duì)比分析結(jié)果表明:該模型與試驗(yàn)曲線吻合較好,具有一定的優(yōu)越性。

References)

[1] 楊友卿.巖石強(qiáng)度的損傷力學(xué)分析[J].巖石力學(xué)與工程學(xué)報(bào),1999,18(1):23-27.

YANG You-qing.Continuum Damage Mechanics Analysis on Strength of Rock[J].Chinese Journal of Rock Mechanics and Engineering,1999,18(1):23-27.(in Chinese)

[2] 曹文貴,方祖烈,唐學(xué)軍.巖石損傷軟化統(tǒng)計(jì)本構(gòu)模型之研究[J].巖石力學(xué)與工程學(xué)報(bào),1998,17(6):628-633.

CAO Wei-gui,FANG Zu-lie,TANG Xue-jun.A Study on Statistical Constitutive Model for Soft and Damage Rocks[J].Chinese Journal of Rock Mechanics and Engineering,1998,17(6):628-633.(in Chinese)

[3] Lemaitre J.A Continuous Damage Mechanics Model for Ductile Fracture[J].J Eng Mater Tech,1985,107:83-89.

[4] 唐春安.巖石破裂過(guò)程中的災(zāi)變[M].北京:煤炭工業(yè)出版社,1993.

TANG Chun-an.Catastrophe in Rock Unstable Failure[M].Beijing:China Coal Industry Publishing House,1993.(in Chinese)

[5] Krajcinovic D,Silva M A G .Statistical Aspects of the Continuous Damage Theory[J].Int J Solids Structures,1982,18(7):551-562.

[6] Tang C A,Chen Z H,Xu X H,et al.A Theoretical Model for Kaiser Effect in Rock[J].Pure Appl Geophys,1997,150:203-215.

[7] 徐衛(wèi)亞,韋立德.演示損傷統(tǒng)計(jì)本構(gòu)模型的研究[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(6):787-791.

XU Wei-ya,WEI Li-de.Study on Statistical Damage Constitutive Model of Rock[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(6):787-791.(in Chinese)

[8] 謝和平.巖石混凝土損傷力學(xué)[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,1990.

XIE He-ping.Damage Mechanics of Rocks and Concrete[M].Xuzhou:China University of Mining and Technology Press,1990.(in Chinese)

[9] 林卓英,吳玉山,關(guān)玲莉,等.巖石在三軸壓縮下脆—延性轉(zhuǎn)化的研究[J].巖土力學(xué),1992,13(2/3):45-53.

LIN Zhuo-ying,WU Yu-shan,GUAN Ling-li,et al.Research on the Brittle-ductile Transition Property of Rocks under Tri-axial Compression[J].Rock and Soil Mechanics,1992,13(2/3):45-53.(in Chinese)

A Modified Damage Softening Constitutive Model for Rock

LI Bo, WU Run-jiang, GAO Wei-chao

(NorthwestInstituteofNuclearTechnology,Xi'an710024,Shaanxi,China)

Using a rock micro-unit strength,which satisfies the Weibull random distribution and is based on the Drucker-Prager failure criterion,a correction coefficient for damage variables is introduced,which reflects the full stress-strain process and residual strength of rock.The relationship between Weibull distribution parameters under different confining pressures is established using the stress and strain at peak points on the stress-strain curve.A damage softening constitutive model of rock is then established,accounting for the rock triaxial stress-strain test curve.The model has few parameters,and the way of determining the parameters indicates the physical meaning of the model parameters.Compared with previous studies,this model,after modification,has obvious advantages.The rationality of the model proposed here is demonstrated using comparative analysis between the theoretical and experimental results.

rock; damage; softening; modified; constitutive model

2015-05-31

李 博(1986-)，男，助理工程師，主要從事巖土工程方面的研究。E-mail：libocare@126.com。

吳潤(rùn)江，主要從事構(gòu)造地質(zhì)學(xué)方面的研究。E-mail：future_wrj@yahoo.com.cn。

TU452

A

1000-0844(2016)05-0783-04

10.3969/j.issn.1000-0844.2016.05.0783