吳朝義

摘要：為研究巖石在循環(huán)加卸載作用下?lián)p傷信息與巖石破壞之間存在的關(guān)系，開(kāi)展了針對(duì)大理巖試件的室內(nèi)力學(xué)試驗(yàn)，結(jié)合彈性模量法、Felicity比值、加卸載響應(yīng)比（LURR）三種損傷信息來(lái)對(duì)巖石的損傷程度進(jìn)行判定，并在時(shí)間尺度下進(jìn)行分析。結(jié)果表明：三種損傷信息均對(duì)巖石破壞變形過(guò)程中的損傷情況進(jìn)行了良好的反饋；巖石處在穩(wěn)定損傷階段的時(shí)間大約為總時(shí)間的60%-80%，處在臨近失穩(wěn)的絕對(duì)損傷階段的時(shí)間大約為總時(shí)間的20-40%；三種損傷計(jì)算方法各有優(yōu)劣，綜合利用可對(duì)工程施工進(jìn)行良好的指引。

Abstract： In order to study the relationship between rock damage and rock destruction under cyclic loading and unloading， an indoor mechanics test was conducted for marble samples， combining elastic modulus method， Felicity ratio， and loading and unloading response ratio （LURR）. The damage information is used to determine the degree of rock damage and analyze it on the time scale. The results show that the three kinds of damage information provide good feedback on the damage conditions during rock failure and deformation； the time for the rock to be in the stable damage stage is approximately 60% to 80% of the total time. The time of the absolute damage stage near the instability is about 20-40% of the total time； Three kinds of damage calculation methods have their own advantages and disadvantages. Comprehensive utilization can provide good guidance for project construction.

關(guān)鍵詞：循環(huán)加卸載；損傷信息；巖石破壞；時(shí)間尺度

Key words： cyclic loading and unloading；damage information；rock destruction；time scale

中圖分類號(hào)：P581 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1006-4311（2018）28-0205-03

0 引言

巖石破壞機(jī)理研究與巖石損傷有著不可分割的關(guān)系。關(guān)于巖石損傷方面的研究，可為巖石破壞機(jī)理和巖石內(nèi)部微觀變化做出很好的鋪墊。

隨著我國(guó)公路設(shè)施建設(shè)的高速發(fā)展，地表路基及深埋隧道等工程經(jīng)常會(huì)遇到人為的反復(fù)的加卸荷作用，這種作用會(huì)致使巖土工程內(nèi)部產(chǎn)生變形、損傷甚至是失穩(wěn)，因此對(duì)于巖石（體）損傷程度的研究十分具有現(xiàn)實(shí)工程意義。周輝[1]等在考慮未損傷部分的應(yīng)力應(yīng)變?nèi)匀粷M足胡克定律的同時(shí)，將Rabotnov損傷參量和Lemaitre應(yīng)變等效性假設(shè)，創(chuàng)造了簡(jiǎn)化后的關(guān)于硬巖的損傷本構(gòu)方程；李樹(shù)春[2]等運(yùn)用損傷力學(xué)理論，立足于巖石疲勞損傷，建立了一種新的巖石疲勞損傷變量表達(dá)式，并針對(duì)砂巖進(jìn)行室內(nèi)力學(xué)加以驗(yàn)證，證明了所建立的模型可良好的反映巖石疲勞損傷演化規(guī)律；付宏淵[3]等在引入損傷因子的同時(shí)根據(jù)經(jīng)典彈塑性理論創(chuàng)建了一種新的損傷耦合模型，并通過(guò)有限元分析，得出在同等初始圍壓下，卸圍壓過(guò)程中的最大損傷因子全部大于常規(guī)試驗(yàn)；劉飛躍[4]等從聲發(fā)射入手，將力學(xué)試驗(yàn)與數(shù)值模擬相結(jié)合，得出了一種真實(shí)的巖石破裂動(dòng)態(tài)應(yīng)力場(chǎng)搜尋方法，對(duì)巖石的損傷破壞過(guò)程研究具有良好的指導(dǎo)意義。

本文在前人研究的基礎(chǔ)上，利用電液伺服試驗(yàn)機(jī)就大理巖開(kāi)展了不同應(yīng)力路徑下的單軸循環(huán)加卸載試驗(yàn)，實(shí)時(shí)收集應(yīng)力、應(yīng)變、聲發(fā)射等多種基本信息，采用包括彈性模量法、Felicity比值法、加卸載響應(yīng)比法共同探索巖石破壞的前兆損傷信息，并完成了對(duì)其基于時(shí)間尺度下的分析與研究。

1 試驗(yàn)準(zhǔn)備

1.1 試驗(yàn)設(shè)備和試件制備

試驗(yàn)所用儀器為T(mén)AW-2000D微機(jī)控制電液伺服巖石三軸試驗(yàn)機(jī)，該試驗(yàn)機(jī)主要由門(mén)框式剛性主機(jī)、圍壓加載系統(tǒng)、自平衡壓力室、控制柜、計(jì)算機(jī)等幾部分組成，可滿足單軸壓縮、三軸壓縮、剪切等多種試驗(yàn)。同時(shí)，搭載由北京聲華科技有限公司生產(chǎn)的SDAES數(shù)字聲發(fā)射檢測(cè)儀，該設(shè)備可實(shí)現(xiàn)多通道聲波收集且波形穩(wěn)定。

試驗(yàn)試件為云南地區(qū)白色大理巖，按照相關(guān)國(guó)際力學(xué)試驗(yàn)規(guī)范，將試件加工成直徑為50mm，高為100mm的標(biāo)準(zhǔn)圓柱體。同時(shí)，為保證較小的離散型，試驗(yàn)所用試件均取自于同一塊大理巖巖塊，且擇優(yōu)挑選光滑無(wú)節(jié)理的試件作為試驗(yàn)對(duì)象。

1.2 試驗(yàn)方案

為探究單軸循環(huán)加卸載作用下大理巖多種損傷信息的變化情況，在完成了單軸壓縮試驗(yàn)，測(cè)定大理巖的單軸抗壓強(qiáng)度、彈性模量、泊松比等一系列參數(shù)后，設(shè)定以下2種試驗(yàn)方案：

①應(yīng)力上限逐級(jí)增加下限恒定的循環(huán)加卸載試驗(yàn)（恒下限）。首先，施加2kN的初始荷載使試件與試驗(yàn)機(jī)緊密接觸，然后進(jìn)行以10MPa為應(yīng)力增量的循環(huán)加卸載試驗(yàn)，且每次均卸載至2kN，以此往復(fù)直至巖石破壞。加載速率為500N/s，卸載速率為300N/s。

②應(yīng)力上下限恒定增幅的循環(huán)加卸載試驗(yàn)（恒差值）。待施加2kN的初始荷載后，仍以10MPa的應(yīng)力增量進(jìn)行加載，然而在卸載階段，除了第一個(gè)循環(huán)卸載至2kN荷載處，其他所有循環(huán)均卸載至上一循環(huán)的最大荷載處，加卸速率同上。

2 損傷結(jié)果分析

為研究巖石破壞過(guò)程當(dāng)中損傷演化特征的變化情況，引入彈性模量法、Felicity比、加卸載響應(yīng)比三種損傷信息作為研究對(duì)象，并且分階段、分時(shí)間的對(duì)巖石破壞整個(gè)過(guò)程的損傷信息來(lái)進(jìn)行探討。

2.1 彈性模量法結(jié)果分析

謝和平[5]等在完善了前人的研究成果，給出了考慮了不可逆變形的更適合巖石這類彈塑性材料的一維損傷公式：

式中：ε為軸向應(yīng)變；εe為循環(huán)荷載下卸載后的殘余塑性變形；E'、E分別為彈塑性材料的卸載彈性模量和加載彈性模量。

圖1為利用彈性模量法計(jì)算出來(lái)的損傷變化情況圖，可以看出利用彈性模量法所得出的損傷變量數(shù)值具有明顯的階段性特征，隨著循環(huán)數(shù)目的不斷增多，損傷變量也不斷增大。按數(shù)值變化情況，將其損傷過(guò)程分為3個(gè)階段：初始損傷階段、穩(wěn)定損傷階段、絕對(duì)損傷階段。

初始損傷階段：第1、2個(gè)循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的0-25%，此階段巖石破壞階段由巖石壓密期漸入線彈性階段，損傷變量示數(shù)稍有跳動(dòng)或初始值較高但數(shù)值較小，說(shuō)明巖石內(nèi)部出現(xiàn)損傷但并不劇烈；

穩(wěn)定損傷階段：第2至第5、6個(gè)循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的25-85%，此階段損傷變量在數(shù)值上所體現(xiàn)出來(lái)的最大特征便是較為穩(wěn)定，發(fā)展速度一般，此現(xiàn)象B組試件反應(yīng)尤為明顯。此階段雖所經(jīng)歷的循環(huán)個(gè)數(shù)最多，但其累積的損傷卻不及整體損傷的30%；

絕對(duì)損傷階段：第5、6個(gè)循環(huán)至最終循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的70-100%，此階段的損傷變量數(shù)值出現(xiàn)跳躍性的上升且增速極大，這種突變式的增長(zhǎng)代表著巖石內(nèi)部損傷程度的嚴(yán)重加深，也預(yù)示著巖石破裂失穩(wěn)的即將發(fā)生。

2.2 Felicity比值結(jié)果分析

聲發(fā)射是應(yīng)力波的一種，它是伴隨著巖石部分區(qū)域顆粒間的作用而產(chǎn)生的，包含了巖石在破壞過(guò)程中內(nèi)部損傷的豐富信息[6]。利用聲發(fā)射技術(shù)來(lái)衡量巖石損傷的方法較多，F(xiàn)elicity比就是其中之一。

Felicity比值對(duì)于損傷的刻畫(huà)十分直觀，其數(shù)值越小代表著損傷程度越大，具體公式如下：

式中：Fσ代表Felicity比值，σAE代表再次聲發(fā)射的起始應(yīng)力，σmax代表上次加載的最大應(yīng)力。

圖2為兩組試驗(yàn)的Felicity比值結(jié)果示意圖。利用Felicity比所得損傷結(jié)果具有良好的變化特征，隨著循環(huán)數(shù)目的增多，F(xiàn)elicity比逐漸減少，整體上呈現(xiàn)“快速跌落-緩慢減少-繼續(xù)跌落”的特點(diǎn)，同樣分為初始損傷、穩(wěn)定損傷、絕對(duì)損傷3個(gè)階段進(jìn)行分析。

初始損傷階段：第2至第4個(gè)循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的0-30%，兩組試驗(yàn)的Felicity比值變化頗為一致，均以快速跌落為主，數(shù)值變化范圍在0.4-0.6左右，說(shuō)明在低應(yīng)力水平少循環(huán)的情況下聲發(fā)射現(xiàn)象仍然較為活躍，巖石內(nèi)部損傷情況變化較大，微顆粒之間的作用也較為劇烈。

穩(wěn)定損傷階段：第4至第6個(gè)循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的30-60%，此階段巖石進(jìn)入彈性期，兩組數(shù)據(jù)均維持在0.1-0.2左右的波動(dòng)，損傷參量波動(dòng)較小。巖石內(nèi)部微顆粒間的作用雖有發(fā)展但較為溫和。

絕對(duì)損傷階段：第6至最后循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的60-100%，此階段巖石性質(zhì)由彈轉(zhuǎn)塑，各組數(shù)據(jù)再次進(jìn)入急速下降的趨勢(shì)，數(shù)值變化范圍大致在0.1-0.3左右，損傷程度逐漸增大，巖石內(nèi)部活動(dòng)十分劇烈。對(duì)比兩組試驗(yàn)數(shù)據(jù)可發(fā)現(xiàn)，A組試件Felicity比值的最終值接近0.3左右，而B(niǎo)組試件Felicity比值的最終值只接近0.9左右，說(shuō)明A組試驗(yàn)的損傷程度更深。

2.3 加卸載響應(yīng)比結(jié)果分析

加卸載響應(yīng)比（LURR）是尹祥礎(chǔ)[7-8]提出的一種地震預(yù)報(bào)方法，是一種可定量分析非線性系統(tǒng)失穩(wěn)趨勢(shì)的參數(shù)。加卸載響應(yīng)比主要思路為：系統(tǒng)處于穩(wěn)定狀態(tài)時(shí)的加卸兩階段的響應(yīng)量的比值與處于失穩(wěn)狀態(tài)時(shí)的加卸兩階段的響應(yīng)量的比值差別較大。本文，我們選擇彈性模量作為響應(yīng)量來(lái)反映巖石損傷情況。

對(duì)處于循環(huán)荷載下的系統(tǒng)，可將其加卸載響應(yīng)比Y定義為：

加卸載響應(yīng)比具有更加明顯的數(shù)值波動(dòng)情況，呈現(xiàn)出“相對(duì)穩(wěn)定—急速驟升—迅速跌落”的數(shù)值狀態(tài)，按照大理巖LURR數(shù)值變化情況，將其損傷階段分為穩(wěn)定損傷和絕對(duì)損傷2個(gè)階段。

穩(wěn)定損傷階段：第1至第5個(gè)循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的0-80%，此階段巖石經(jīng)歷了壓密期、線彈性期及部分彈性期，所有試件的加卸載響應(yīng)比參數(shù)基本在“1”左右徘徊，說(shuō)明巖石在此階段的損傷程度較小，體現(xiàn)為巖石內(nèi)部顆粒開(kāi)始萌生、發(fā)育、摩擦但較為溫和，微裂隙進(jìn)行小幅度張開(kāi)或閉合；

絕對(duì)損傷階段：第6至最終循環(huán)，對(duì)應(yīng)時(shí)間段大約為巖石破壞總時(shí)間的80-100%，此階段巖石經(jīng)歷了部分彈性期及最后的塑性期，所有試件的加卸載響應(yīng)比參數(shù)均出現(xiàn)跳躍性突增，說(shuō)明巖石在此階段的損傷程度較大，即將面臨失穩(wěn)，體現(xiàn)為巖石內(nèi)部顆粒的衍生十分活躍，微顆粒間的作用十分劇烈。

3 結(jié)論

①本文利用彈性模量法、Felicity比值法、加卸載響應(yīng)比法三種方法共同對(duì)大理巖在循環(huán)加卸荷作用下的損傷情況進(jìn)行討論，并分階段分時(shí)段的對(duì)巖石損傷情況進(jìn)行分析，結(jié)果證明三種方法均對(duì)巖石變形破壞過(guò)程的損傷情況所有反饋，巖石處在穩(wěn)定損傷階段的時(shí)間大約為總時(shí)間的60-80%，處在臨近失穩(wěn)的絕對(duì)損傷階段的時(shí)間大約為總時(shí)間的20-40%；

②對(duì)比兩組試驗(yàn)數(shù)據(jù)可發(fā)現(xiàn)，恒下限試驗(yàn)組巖石的損傷結(jié)果大于恒差值試驗(yàn)組，造成這種情況的原因在于，恒下限試驗(yàn)組的加卸跨度較大，加卸時(shí)間較長(zhǎng)，從而造成巖石內(nèi)部微顆粒的作用更強(qiáng)；

③三種損傷前兆信息各有特色和優(yōu)點(diǎn)，彈性模量法具有明顯的階段性，F(xiàn)elicity比值法利用了具有預(yù)兆能力的聲發(fā)射技術(shù)，加卸載響應(yīng)比法可產(chǎn)生強(qiáng)烈的數(shù)據(jù)跳動(dòng)，在實(shí)際工程當(dāng)中，應(yīng)將多種方法相結(jié)合以便做出更加精準(zhǔn)的損傷評(píng)估。

參考文獻(xiàn)：

[1]周輝，李震，朱國(guó)金，等.基于巖石統(tǒng)一能量屈服準(zhǔn)則的硬巖損傷模型[J].巖土力學(xué)，2016，37（3）：609-615.

[2]李樹(shù)春，許江，陶云奇，等.巖石低周疲勞損傷模型與損傷變量表達(dá)方法[J].巖土力學(xué)，2009，30（6）：1611-1614.

[3]付宏淵，邱祥，史振寧，等.基于彈塑性損傷模型的巖石加卸載計(jì)算研究[J].地下空間與工程學(xué)報(bào)，2018（1）.

[4]劉飛躍，楊天鴻，張鵬海，等.基于聲發(fā)射的巖石破裂應(yīng)力場(chǎng)動(dòng)態(tài)反演[J].巖土力學(xué)，2018（04）：1-7.

[5]謝和平，董毓利.經(jīng)典損傷定義中的“彈性模量法”探討[J].力學(xué)與實(shí)踐，1997，19（2）：1-5.

[6]劉建坡，徐世達(dá)，李元輝，等.預(yù)制孔巖石破壞過(guò)程中的聲發(fā)射時(shí)空演化特征研究[J].巖石力學(xué)與工程學(xué)報(bào)，2012，31（12）：2538-2547.

[7]尹祥礎(chǔ)，陳學(xué)忠，宋治平，等.加卸載響應(yīng)比——一種新的地震預(yù)報(bào)方法[J].地球物理學(xué)報(bào)，1994，37（6）：767-775.

[8]Yin X C， Zhang L P， Zhang H H， et al. LURR's Twenty Years and its Perspective[J]. Pure & Applied Geophysics， 2007，163 （11-12）：2317-2341.