基于裂紋擴(kuò)展脆性巖石動(dòng)力壓縮力學(xué)特性研究*

李曉照，戚承志

（1. 北京建筑大學(xué)土木與交通工程學(xué)院，北京 100044；2. 北京未來(lái)城市設(shè)計(jì)高精尖創(chuàng)新中心，北京 100044）

動(dòng)力壓縮荷載作用脆性巖石力學(xué)特性的研究，對(duì)于深部地下工程圍巖受到地震、沖擊及爆炸荷載作用下的穩(wěn)定性分析有著重要的意義。關(guān)于應(yīng)變率對(duì)巖石力學(xué)特性影響的研究，對(duì)于巖石動(dòng)態(tài)力學(xué)性能分析有著重要價(jià)值。眾多學(xué)者，通過(guò)開(kāi)展霍布金森桿動(dòng)力試驗(yàn)，研究了不同應(yīng)變率下巖石的應(yīng)力-應(yīng)變本構(gòu)曲線、巖石動(dòng)態(tài)強(qiáng)度、動(dòng)態(tài)破壞形態(tài)等特性[1-4]。也有學(xué)者研究了在恒定軸壓與圍壓作用下，巖石承受循環(huán)重復(fù)沖擊荷載作用下的力學(xué)性能[5-7]。

在這些試驗(yàn)研究的基礎(chǔ)上，許多理論模型被提出。曹文貴等[8]通過(guò)將非線性動(dòng)態(tài)強(qiáng)度準(zhǔn)則與統(tǒng)計(jì)損傷理論相結(jié)合，建立巖石動(dòng)態(tài)統(tǒng)計(jì)損傷本構(gòu)模型，并根據(jù)試驗(yàn)結(jié)果確定模型參數(shù)，計(jì)算分析了高應(yīng)變對(duì)巖石動(dòng)態(tài)本構(gòu)關(guān)系的影響。李夕兵等[9-10]將統(tǒng)計(jì)損傷模型和黏彈性模型相結(jié)合，提出了中應(yīng)變率作用下巖石動(dòng)載本構(gòu)模型，并分析了重復(fù)沖擊荷載作用下巖石的動(dòng)力特性。謝理想等[11]基于巖石動(dòng)態(tài)強(qiáng)度理論模型、連續(xù)損傷理論及Bingham 模型，蔡燦等[12]基于Maxwell 模型、Bingham 模型及孔隙度損傷模型分別建立了不同組合模型，研究了巖石應(yīng)變率相關(guān)的應(yīng)力應(yīng)變本構(gòu)關(guān)系。戚承志等[13-14]基于熱活化機(jī)制與宏觀黏性機(jī)制，研究了應(yīng)變率與巖石強(qiáng)度之間關(guān)系，并分析了動(dòng)態(tài)荷載作用下巖石破壞塊度分布物理機(jī)理。趙堅(jiān)等[15]、馬林建等[16]基于動(dòng)態(tài)試驗(yàn)結(jié)果并結(jié)合摩爾庫(kù)倫及霍克布朗強(qiáng)度準(zhǔn)則，提出了應(yīng)變率相關(guān)的動(dòng)態(tài)強(qiáng)度準(zhǔn)則。宮鳳強(qiáng)等[17]基于巖石從低到高應(yīng)變率范圍內(nèi)的率效應(yīng)試驗(yàn)結(jié)果，建立了一種基于率效應(yīng)動(dòng)態(tài)增強(qiáng)因子統(tǒng)一模型，探討了應(yīng)變率對(duì)壓縮強(qiáng)度、切線模量和峰值應(yīng)變的影響。以上均是通過(guò)宏觀現(xiàn)象研究了巖石動(dòng)態(tài)變化特性，無(wú)法解釋動(dòng)態(tài)荷載作用下，巖石內(nèi)部細(xì)觀裂紋機(jī)理變化特性。

岑奪豐等[18]通過(guò)細(xì)觀顆粒平行黏結(jié)模型的模擬，研究了高應(yīng)變率壓縮作用，單裂隙巖樣的損傷演化及細(xì)觀位移場(chǎng)，結(jié)果表明不同高應(yīng)變率對(duì)巖樣最終破裂形態(tài)影響不大，但隨應(yīng)變率的增大，細(xì)觀裂紋越多且局部化程度越強(qiáng)。潘紅宇等[19]通過(guò)用數(shù)字圖像處理技術(shù)，分析了不同應(yīng)變率下單軸壓縮下巖石表面裂隙擴(kuò)展下的應(yīng)變場(chǎng)變化特征，然而該研究只是分析準(zhǔn)靜態(tài)壓縮荷載作用下的裂紋特性，沒(méi)有研究動(dòng)態(tài)荷載下巖石裂紋擴(kuò)展特性。劉紅巖等[20]將基于巖石動(dòng)態(tài)損傷TCK 模型與基于斷裂損傷理論推出的節(jié)理巖體損傷本構(gòu)模型相結(jié)合，提出了一個(gè)考慮裂紋節(jié)理幾何特性的巖體動(dòng)態(tài)損傷本構(gòu)模型。

然而對(duì)能夠建立脆性巖石內(nèi)部細(xì)觀動(dòng)態(tài)裂紋擴(kuò)展與宏觀動(dòng)態(tài)應(yīng)力-應(yīng)變本構(gòu)關(guān)系的理論模型研究較少。本文中，將準(zhǔn)靜態(tài)壓縮荷載作用下裂紋擴(kuò)展誘發(fā)的應(yīng)力應(yīng)變本構(gòu)關(guān)系模型、基于裂紋與應(yīng)變關(guān)系得到的裂紋速率與應(yīng)變率關(guān)系及裂紋速率與斷裂韌度關(guān)系相結(jié)合，提出了一個(gè)基于細(xì)觀裂紋擴(kuò)展作用下的動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系模型，并討論巖石應(yīng)變率、初始損傷、圍壓及模型參數(shù)對(duì)應(yīng)力應(yīng)變關(guān)系曲線的影響。

1 理論模型

1.1 準(zhǔn)靜態(tài)裂紋擴(kuò)展應(yīng)力-應(yīng)變關(guān)系

本文中提出了一個(gè)動(dòng)態(tài)壓縮荷載作用下，基于巖石內(nèi)部裂紋擴(kuò)展誘發(fā)的應(yīng)力-應(yīng)變本構(gòu)關(guān)系模型（見(jiàn)圖1），模型中各項(xiàng)參數(shù)詳見(jiàn)下面關(guān)于模型理論公式的敘述。該模型是由準(zhǔn)靜態(tài)作用下裂紋擴(kuò)展模型，與動(dòng)態(tài)裂紋擴(kuò)展與準(zhǔn)靜態(tài)裂紋擴(kuò)展關(guān)系，及新建立的裂紋擴(kuò)展速率與應(yīng)變率關(guān)系相結(jié)合而提出的。在該模型中，當(dāng)不考慮動(dòng)態(tài)裂紋擴(kuò)展作用時(shí)，模型內(nèi)部準(zhǔn)靜態(tài)裂紋擴(kuò)展與承受應(yīng)力狀態(tài)（軸向應(yīng)力σ1、圍壓σ3）之間的關(guān)系可以表示為[21]：

圖 1 動(dòng)態(tài)壓縮荷載作用下的裂紋擴(kuò)展模型Fig. 1 Crack growth model under dynamic compressive loadings

此外，低圍壓脆性巖石宏觀軸向應(yīng)變?chǔ)?與細(xì)觀裂紋擴(kuò)展長(zhǎng)度l 間的關(guān)系可表示為[22]：

式中：m、ε0分別為材料常數(shù)。

將式(7)代入式(1)，裂紋擴(kuò)展下（即l＞0），巖石準(zhǔn)靜態(tài)應(yīng)力-應(yīng)變本構(gòu)關(guān)系可表示為[22]：

式中：

此外，當(dāng)裂紋長(zhǎng)度l=0 時(shí)，基于式（1）和（7），巖石軸向起裂應(yīng)力σ1c和裂紋起裂應(yīng)變?chǔ)?c可以分別表示為：

假設(shè)巖石內(nèi)部裂紋擴(kuò)展前（即0＜ ε1＜ε1c），巖石服從線彈性變形，則應(yīng)力-應(yīng)變關(guān)系可表示為：

因此，綜上所述，準(zhǔn)靜態(tài)荷載作用下，巖石完整的應(yīng)力-應(yīng)變本構(gòu)關(guān)系可以表示為：

將裂紋極限長(zhǎng)度llim代入式（7），巖石極限軸向應(yīng)變?chǔ)?f可以表示為：

該極限應(yīng)變對(duì)巖石失效應(yīng)變提供一定參考。

1.2 準(zhǔn)靜態(tài)裂紋擴(kuò)展與動(dòng)態(tài)裂紋擴(kuò)展的關(guān)系

式中：vl為裂紋速度（即dl/dt），vm為裂紋分叉速度，vm=RcR，該速度隨著材料柔度增加而減小，瑞利波速cR= [E/(2ρ(1+ν))]1/2(0.87+1.12ν)/(1+ν)，cp為P 波速度，cp=[(1-ν)E/((1+ν)(1-2ν)ρ)]1/2，式中ν 為泊松比，ρ 為材料密度，E 為彈性模量，R 是材料常數(shù)（0＜R＜1）。

基于式（7）可以解出裂紋擴(kuò)展長(zhǎng)度表達(dá)式為：

對(duì)式（19）裂紋長(zhǎng)度求時(shí)間導(dǎo)數(shù)，則可以推出巖石內(nèi)部細(xì)觀裂紋速率vl為：

將式（20）代入式（18），則可以建立動(dòng)態(tài)斷裂韌度與宏觀軸向應(yīng)變率之間關(guān)系：

式中：

1.3 動(dòng)態(tài)裂紋擴(kuò)展下的應(yīng)力-應(yīng)變關(guān)系

式（24）中E(d)為動(dòng)彈性模量，下面將對(duì)模型參數(shù)對(duì)動(dòng)彈性模量影響進(jìn)行詳細(xì)討論。此外，值得注意的是在動(dòng)彈性模量E(d)的方程計(jì)算中，參數(shù)Q 中的軸向應(yīng)變?chǔ)?應(yīng)取值為ε1c。

2 結(jié)果與討論

2.1 動(dòng)態(tài)應(yīng)力-應(yīng)變關(guān)系理論結(jié)果合理性驗(yàn)證

圖 2 大理巖與與花崗巖動(dòng)態(tài)應(yīng)力-應(yīng)變曲線的理論與試驗(yàn)結(jié)果[3, 23]對(duì)比Fig. 2 Comparisons between theoretical and experimental [3, 23] dynamic stress-strain curves in marble and granite

圖2 中分別給出基于大理巖與花崗巖材料屬性，得到的動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系的理論與試驗(yàn)結(jié)果[3,23]曲線對(duì)比結(jié)果。可以看到，理論結(jié)果與試驗(yàn)結(jié)果在變化趨勢(shì)上具有較好的相似性，在數(shù)值上具有一定的可比性，但是仍然存在一定的差異性。其中導(dǎo)致這些差異的原因可以分析為以下幾點(diǎn)：（1）本文中提出的理論模型是一個(gè)基于各向同性的均勻彈性體得到的；（2）巖石內(nèi)部隨機(jī)分布裂紋無(wú)法真實(shí)呈現(xiàn)，而是通過(guò)假設(shè)裂紋尺寸、角度相同的均布化思路進(jìn)行計(jì)算分析的。

2.2 應(yīng)變率對(duì)應(yīng)力-應(yīng)變關(guān)系、動(dòng)態(tài)壓縮強(qiáng)度和動(dòng)態(tài)彈性模量的影響

圖 3 應(yīng)變率影響下的應(yīng)力應(yīng)變本構(gòu)關(guān)系Fig. 3 Effect of strain rate on stress-strain curve

以上通過(guò)試驗(yàn)結(jié)果驗(yàn)證了本文基于細(xì)觀力學(xué)的動(dòng)態(tài)本構(gòu)模型合理性。下面將詳細(xì)討論應(yīng)變率對(duì)巖石力學(xué)特性（動(dòng)態(tài)壓縮強(qiáng)度和動(dòng)態(tài)彈性模量）的影響。其中模型參數(shù)選取參考花崗巖參數(shù)。圖3中給出了不同應(yīng)變率下，應(yīng)力-應(yīng)變本構(gòu)關(guān)系的變化趨勢(shì)。可以看到，隨著應(yīng)變率的增加，應(yīng)力應(yīng)變曲線峰值應(yīng)力，即壓縮強(qiáng)度，也在不斷增大，而且是加速增加，該現(xiàn)象也可以通過(guò)圖4 給出的軸向動(dòng)態(tài)壓縮強(qiáng)度與軸向應(yīng)變率關(guān)系清楚地看到，該結(jié)果已經(jīng)被眾多研究動(dòng)態(tài)強(qiáng)度理論的學(xué)者驗(yàn)證[1,3, 13, 15,23]。還可看到，隨著應(yīng)變率的增大，巖石動(dòng)態(tài)強(qiáng)度對(duì)應(yīng)峰值軸向應(yīng)變不斷減小。圖4 中，研究了不同巖石內(nèi)部初始損傷、模型參數(shù)m、ε0、R 下，巖石軸向動(dòng)態(tài)壓縮強(qiáng)度與軸向應(yīng)變率的聯(lián)系。隨著初始損傷D0或參數(shù)ε0、R 的減小、參數(shù)m 的增大，在一定軸向應(yīng)變率下軸向動(dòng)壓縮強(qiáng)度增大。

圖 4 不同模型參數(shù)的應(yīng)變率與巖石動(dòng)態(tài)壓縮強(qiáng)度關(guān)系Fig. 4 Relations between strain rate and dynamic compressive strength under different model parameters

此外，還可以看到，巖石裂紋初始應(yīng)變值ε1c隨著應(yīng)變率增加為一個(gè)恒定值，而裂紋初始應(yīng)變對(duì)應(yīng)的裂紋初始軸向應(yīng)力則是隨著應(yīng)變率的增加也在增大。本文中假設(shè)裂紋初始應(yīng)變前的變形為線彈性變形。隨著應(yīng)變率增加，巖石彈性變形階段曲線的切線模量，即動(dòng)態(tài)彈性模量，在不斷加速增大，該現(xiàn)象也可以在圖5 中動(dòng)態(tài)彈性模量與應(yīng)變率之間關(guān)系看到，并且與眾多學(xué)者的試驗(yàn)結(jié)果[1-3,23]相一致。圖5 研究了不同圍壓、巖石內(nèi)部初始損傷及模型參數(shù)m、ε0、R 下，巖石動(dòng)態(tài)彈性模量與軸向應(yīng)變率之間的聯(lián)系。隨著初始損傷D0或參數(shù)m、ε0、R 減小、圍壓增大，在一定軸向應(yīng)變率下動(dòng)態(tài)彈性模量不斷增大。

圖 5 不同模型參數(shù)的應(yīng)變率與動(dòng)態(tài)彈性模量關(guān)系Fig. 5 Relations between strain rate and dynamic elastic module under different model parameters

2.3 模型參數(shù)敏感性分析

巖石內(nèi)部初始損傷對(duì)巖石力學(xué)性能有著重要的影響。圖6 中給出了不同初始損傷對(duì)于高應(yīng)變率300 s?1下動(dòng)態(tài)應(yīng)力-應(yīng)變本構(gòu)關(guān)系曲線的影響。可以看到，隨著損傷的增大，巖石的動(dòng)態(tài)壓縮強(qiáng)度不斷減小，而且隨著損傷的減小，峰值強(qiáng)度對(duì)應(yīng)的峰值應(yīng)變?cè)诓粩鄿p小。圖7 給出了不同圍壓狀態(tài)下，巖石動(dòng)態(tài)應(yīng)力-應(yīng)變關(guān)系曲線結(jié)果。可以看到，隨著圍壓增大，動(dòng)態(tài)壓縮強(qiáng)度也在不斷增大。此外，圖8、9 分別給出了模型中參數(shù)m 和ε0對(duì)動(dòng)態(tài)應(yīng)力-應(yīng)變關(guān)系曲線的影響，可以看到，隨著參數(shù)m 和ε0減小，峰值強(qiáng)度對(duì)應(yīng)的峰值應(yīng)變值在不斷減小。圖10 中給出了參數(shù)R 對(duì)應(yīng)力-應(yīng)變曲線的影響，隨著R 的減小，峰值應(yīng)力對(duì)應(yīng)的峰值應(yīng)變值都在不斷減小。

圖 6 初始損傷D0 對(duì)動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系影響Fig. 6 Effect of initial damage on dynamic stress-strain constitutive relation

圖 7 圍壓對(duì)動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系影響Fig. 7 Effect of confining pressure on dynamic stress-strain constitutive relation

圖 9 參數(shù)ε0 對(duì)動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系影響Fig. 9 Effect of parameter ε0 on dynamic stress-strain constitutive relation

圖 10 參數(shù)R 對(duì)動(dòng)態(tài)應(yīng)力應(yīng)變本構(gòu)關(guān)系影響Fig. 10 Effect of parameter R on dynamic stress-strain constitutive relation

圖11 中，給出了不同應(yīng)變率下，動(dòng)態(tài)彈性模量與圍壓、初始損傷和參數(shù)m、ε0、R 的關(guān)系。可以看到，隨著圍壓的增大，動(dòng)態(tài)彈性模量線性增大；隨著初始損傷和模型參數(shù)m、ε0、R 的增大，動(dòng)態(tài)彈性模量在減速降低。以上關(guān)于模型參數(shù)影響研究，更加清楚詳細(xì)地分析了本文提出的宏細(xì)觀力學(xué)模型中參數(shù)敏感性。進(jìn)而為更加精確地建立真實(shí)巖石與模型關(guān)系之間聯(lián)系提供了一定的參考。

圖 11 應(yīng)變率對(duì)不同模型參數(shù)與動(dòng)態(tài)彈性模量關(guān)系影響Fig. 11 Effects of strain rate on the relationships between model parameters and dynamic elastic modules

3 結(jié) 論

基于脆性巖石準(zhǔn)靜態(tài)裂紋擴(kuò)展作用下的應(yīng)力-應(yīng)變本構(gòu)關(guān)系、裂紋速率與應(yīng)變率的變化關(guān)系及裂紋速率與斷裂韌度的關(guān)系，建立了基于細(xì)觀力學(xué)的動(dòng)態(tài)應(yīng)力-應(yīng)變本構(gòu)關(guān)系模型。通過(guò)試驗(yàn)結(jié)果驗(yàn)證了模型的合理性，并討論了模型參數(shù)對(duì)巖石動(dòng)態(tài)壓縮力學(xué)特性的影響。具體結(jié)果如下：

（1）隨著應(yīng)變率的增大，巖石動(dòng)態(tài)壓縮強(qiáng)度不斷增大，巖石裂紋初始應(yīng)力不斷加速增大，峰值軸向應(yīng)力對(duì)應(yīng)的峰值軸向應(yīng)變值不斷減小，巖石動(dòng)態(tài)彈性模量不斷加速增大。

（2）隨著圍壓，初始損傷D0，參數(shù)m、R 和ε0的減小，應(yīng)力-應(yīng)變曲線中峰值強(qiáng)度對(duì)應(yīng)的峰值應(yīng)變不斷減小。

（3）隨著初始損傷D0或參數(shù)ε0、R 的減小，圍壓或參數(shù)m 的增大，在一定的軸向應(yīng)變率下，軸向動(dòng)態(tài)壓縮強(qiáng)度在增大。

（4）隨著初始損傷D0或參數(shù)m、ε0、R 的減小，圍壓的增大，在一定軸向應(yīng)變率下，動(dòng)彈性模量不斷增大。