杜圣賢,陳軍,于學(xué)峰,賈超,張尚坤,宋香鎖，羅文強(qiáng)

(1.山東省地質(zhì)科學(xué)研究院，國(guó)土資源部金礦成礦過程與資源利用重點(diǎn)實(shí)驗(yàn)室，山東省金屬礦產(chǎn)成礦地質(zhì)過程與資源利用重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250013；2．山東大學(xué)土木工程學(xué)院，山東 濟(jì)南 250013)

0 引言

恐龍化石和其他古生物化石一樣，它是由多種礦物成分所構(gòu)成[1-5]。隨著化石圍巖不同而發(fā)生變化，質(zhì)地也不均勻?；纬珊?，在漫長(zhǎng)的地質(zhì)作用(特別是構(gòu)造運(yùn)動(dòng))過程中，必然會(huì)產(chǎn)生各種裂隙、孔洞和變形等(圖1)。恐龍化石為脆性材料，其變形破壞呈現(xiàn)彈脆性。

為了研究?jī)?nèi)含裂隙角的巖石在受到壓力作用下的破壞規(guī)律，前人做了大量的模擬試驗(yàn)。通過新方法與傳統(tǒng)理論相結(jié)合，研究總結(jié)了裂隙缺陷的擴(kuò)展貫通機(jī)制、巖橋貫穿機(jī)制以及失穩(wěn)模式[6-10]。

圖1 恐龍化石裂隙示意圖

為了探討恐龍化石的風(fēng)化機(jī)理[11-13]，該文以斷裂力學(xué)的巖石強(qiáng)度理論為依據(jù)，對(duì)內(nèi)含裂隙恐龍化石在壓力作用下的破壞情況進(jìn)行數(shù)值模擬，研究?jī)?nèi)含裂隙恐龍化石在各種不同狀態(tài)下壓縮試驗(yàn)的應(yīng)力-應(yīng)變規(guī)律。

1 巖石斷裂力學(xué)介紹

斷裂力學(xué)是近幾十年才發(fā)展起來的一支新興學(xué)科，它從宏觀的連續(xù)介質(zhì)力學(xué)角度出發(fā)，研究含缺陷或裂紋的物體在外界條件作用下宏觀裂紋的擴(kuò)展、失穩(wěn)開裂、傳播和止裂規(guī)律。斷裂力學(xué)成功地研究了含缺陷材料和結(jié)構(gòu)的破壞問題，雖然起步較晚，但是它與材料或結(jié)構(gòu)的安全問題直接相關(guān)，因此實(shí)驗(yàn)與理論發(fā)展迅速，并且在工程上得到了廣泛應(yīng)用[14-17]。

巖石強(qiáng)度理論是一個(gè)復(fù)雜的科學(xué)問題，它的目的在于給出巖石在復(fù)雜應(yīng)力狀態(tài)下失穩(wěn)破壞時(shí)巖石參數(shù)與外部荷載及環(huán)境因素所滿足的條件。建立一種科學(xué)合理的巖石強(qiáng)度理論，對(duì)于工程設(shè)計(jì)、災(zāi)害預(yù)防、資源開發(fā)等領(lǐng)域具有重要意義[18-20]。

早期的巖石強(qiáng)度理論研究屬于傳統(tǒng)固體力學(xué)研究范疇，沿用了經(jīng)典連續(xù)統(tǒng)力學(xué)的研究手段。由于巖石材料特殊，受地質(zhì)構(gòu)造的影響，巖石的組分結(jié)構(gòu)極度不均勻，孔隙、裂隙、夾雜、節(jié)理、斷層等大量缺陷充斥其中，因而均勻連續(xù)假設(shè)與巖石的實(shí)際情況并不相符，建立在連續(xù)統(tǒng)力學(xué)基礎(chǔ)上的巖石強(qiáng)度理論受到了嚴(yán)重挑戰(zhàn)。隨著非線性科學(xué)等相關(guān)學(xué)科的迅猛發(fā)展，巖石力學(xué)的研究融合了經(jīng)典彈塑性力學(xué)、斷裂力學(xué)、損傷力學(xué)、熱力學(xué)、物理學(xué)、化學(xué)、地質(zhì)學(xué)、礦物學(xué)、信息論、控制論、系統(tǒng)論等學(xué)科，使巖石力學(xué)的研究逐漸超越了經(jīng)典固體力學(xué)的框架，極大地豐富了巖石力學(xué)的研究?jī)?nèi)涵。巖石強(qiáng)度理論的研究也逐漸從古典強(qiáng)度理論、廣義強(qiáng)度理論等經(jīng)典強(qiáng)度理論發(fā)展到將斷裂、損傷過程考慮進(jìn)去的強(qiáng)度理論，從宏觀唯象研究發(fā)展到跨尺度多層次的理性研究[18,21-26]。

在古典強(qiáng)度理論的發(fā)展過程中，唯象[27-28]的試驗(yàn)研究是主要手段，古典強(qiáng)度理論僅適用于簡(jiǎn)單應(yīng)力狀態(tài)。隨著實(shí)驗(yàn)測(cè)試手段和數(shù)學(xué)分析方法的發(fā)展，提出了經(jīng)典強(qiáng)度理論[18,29-30]，建立了復(fù)雜應(yīng)力狀態(tài)下固體材料屈服破壞[31-36]的臨界準(zhǔn)則。利用了現(xiàn)代數(shù)學(xué)、力學(xué)的研究成果，克服了古典理論唯象性的缺陷，將理性分析推導(dǎo)與實(shí)驗(yàn)經(jīng)驗(yàn)相結(jié)合，較好地解決了強(qiáng)度理論的計(jì)算問題[18]。

經(jīng)典巖石強(qiáng)度理論采用了連續(xù)介質(zhì)的假定，與巖石材料的實(shí)際不符，經(jīng)典理論未能解決巖石強(qiáng)度的離散性、隨機(jī)性等問題，也沒有回答巖石強(qiáng)度特性與巖石組織結(jié)構(gòu)之間的關(guān)系問題。因此迫切需要從巖石的組織結(jié)構(gòu)出發(fā)，采用新的研究手段，發(fā)展經(jīng)典巖石強(qiáng)度理論,將斷裂力學(xué)引入到巖石強(qiáng)度理論研究中[18,37-38]。

巖石是一種非均質(zhì)的多相復(fù)合結(jié)構(gòu)的材料，天然存在各種大量的缺陷，且這些缺陷的分布具有一定的隨機(jī)性。巖石受到外力作用，巖石內(nèi)部的微缺陷不斷發(fā)生變化，在部分區(qū)域出現(xiàn)貫通，進(jìn)而形成宏觀裂縫導(dǎo)致巖石失穩(wěn)破壞。單純用經(jīng)典彈塑性力學(xué)或斷裂力學(xué)的方法來描述，很難獲得理想結(jié)果，因此在巖石強(qiáng)度理論研究中引入了損傷力學(xué)，建立了宏-細(xì)-微觀多層次耦合的巖石強(qiáng)度理論。巖石細(xì)觀模型和細(xì)觀模型與宏觀整體的映射關(guān)系建立，非線性科學(xué)的應(yīng)用，系統(tǒng)論、控制論、信息論等在巖石力學(xué)研究中的引入，計(jì)算機(jī)技術(shù)在數(shù)值計(jì)算、虛擬現(xiàn)實(shí)等方面的利用，促進(jìn)了巖石強(qiáng)度理論分析及應(yīng)用的新發(fā)展[18,23,25,38-45]。

2 內(nèi)含裂隙恐龍化石數(shù)值壓縮試驗(yàn)

內(nèi)含裂隙恐龍化石數(shù)值模擬試驗(yàn)主要是以斷裂力學(xué)理論等為依據(jù)，在試驗(yàn)的基礎(chǔ)上，根據(jù)已獲得的試驗(yàn)參數(shù)，借助ABAQUS數(shù)值模擬軟件，全面地揭示恐龍化石在加載受力過程中的應(yīng)力、應(yīng)變的變化特征，研究恐龍化石的各項(xiàng)物理力學(xué)性質(zhì)對(duì)不同狀態(tài)下的內(nèi)含裂隙恐龍化石的影響程度，據(jù)此評(píng)價(jià)各種因素對(duì)恐龍化石風(fēng)化的作用效果。

2.1 恐龍化石物理力學(xué)參數(shù)

根據(jù)實(shí)驗(yàn)測(cè)得恐龍化石的物理力學(xué)參數(shù)如表1、表2所示。

表1 化石力學(xué)參數(shù)

表2 化石物理性質(zhì)

2.2 內(nèi)含裂隙恐龍化石軸向受壓數(shù)值試驗(yàn)

本試驗(yàn)是對(duì)恐龍化石在單軸壓力下對(duì)其受力狀態(tài)進(jìn)行測(cè)試，單軸試驗(yàn)數(shù)值采用有限元軟件ABAQUS進(jìn)行模擬。首先建立恐龍化石數(shù)值計(jì)算模型并劃分網(wǎng)格，然后定義力學(xué)參數(shù)、邊界條件，施加荷載，最終獲得恐龍化石的應(yīng)力變化情況。

(1)研究?jī)?nèi)容

數(shù)值試驗(yàn)?zāi)M研究了恐龍化石內(nèi)含裂隙角度分別為15°，30°，45°，60°，75°時(shí)，在軸向受壓情況下恐龍化石的應(yīng)力分布。內(nèi)含裂隙恐龍化石樣品數(shù)值計(jì)算模型如圖2所示，在裂隙兩端處對(duì)網(wǎng)格進(jìn)行了細(xì)分。

(2)結(jié)果分析

裂隙角度變化對(duì)裂紋擴(kuò)展的影響:壓應(yīng)力作用下中心斜裂紋的開裂角和極限載荷的影響如圖3—圖5所示。

圖3 裂隙角度的變化對(duì)開裂角的影響

圖4 裂隙角度的變化對(duì)裂紋擴(kuò)展步長(zhǎng)的影響

圖5 裂隙角度的變化對(duì)極限載荷的影響

由圖中可見，隨著裂隙角度增大，開裂角逐漸減小。裂隙角度的變化對(duì)裂紋擴(kuò)展步長(zhǎng)的影響不大，在裂隙角度大于30°后，裂紋擴(kuò)展步長(zhǎng)減小；裂隙角度大于45°之后，裂紋擴(kuò)展步長(zhǎng)略有增大。裂隙角度在15°～25°時(shí)，極限載荷呈緩降趨勢(shì)；裂隙角度在25°～45°時(shí)，極限載荷逐漸增大；裂隙角度在45°～75°時(shí)，極限載荷逐漸減小。

裂隙角度變化對(duì)應(yīng)力場(chǎng)的影響:軸壓作用下裂隙角度變化對(duì)應(yīng)力分布的影響如圖6所示。由圖中可見，裂隙處特別是兩端應(yīng)力集中，其他部位應(yīng)力較低、分布均勻。

3 內(nèi)含裂隙的恐龍化石裂紋擴(kuò)展研究

3.1 試驗(yàn)理論依據(jù)

恐龍化石組分、結(jié)構(gòu)不均勻，本身又具有很多節(jié)理和裂隙，為非連續(xù)性介質(zhì)。如果采用傳統(tǒng)的材料力學(xué)觀點(diǎn)來計(jì)算化石的承載力，其計(jì)算結(jié)果必將產(chǎn)生很大的誤差。為了更精確了解內(nèi)含裂隙角恐龍化石裂紋擴(kuò)展規(guī)律，本次試驗(yàn)的理論依據(jù)主要為巖石斷裂力學(xué)。直接或間接地采用保角變換法、復(fù)變函數(shù)法、willians展開法、Riemann-Hilbert邊值問題法、Fourier變換法等經(jīng)典方法，對(duì)恐龍化石在壓力作用下裂紋的起裂特性進(jìn)行分析，得出裂紋的起裂規(guī)律、裂紋尖端的位移場(chǎng)、應(yīng)力場(chǎng)分布以及裂尖的應(yīng)力強(qiáng)度因子參數(shù)[46]。

①—裂隙角度15°；②—裂隙角度30°；③—裂隙角度45°；④—裂隙角度60°；⑤—裂隙角度75°圖6 軸壓作用下裂隙角度變化對(duì)應(yīng)力分布的影響

3.2 試驗(yàn)及結(jié)果分析

本試驗(yàn)主要運(yùn)用數(shù)值分析方法對(duì)單軸壓力作用下內(nèi)含裂隙恐龍化石斷裂損傷過程進(jìn)行模擬分析，研究裂紋萌生、擴(kuò)展、貫通規(guī)律和裂隙化石的斷裂損傷機(jī)制。根據(jù)前面的分析結(jié)果可知，裂隙角度為30°時(shí)較為危險(xiǎn)，更易開裂，因此下面以內(nèi)含裂隙角度為30°的恐龍化石作為分析對(duì)象。

裂隙角度為30°時(shí)裂紋擴(kuò)展路徑及應(yīng)力場(chǎng)分布如圖7所示。由圖中可見，裂紋擴(kuò)展從裂紋尖端起裂，預(yù)制裂紋端部產(chǎn)生微小翼型裂紋，當(dāng)繼續(xù)施加荷載時(shí)翼型裂紋變寬延伸，隨后轉(zhuǎn)向上擴(kuò)展，最終裂紋擴(kuò)展到邊緣，產(chǎn)生的次生裂紋明顯呈鋸齒狀，最大應(yīng)力區(qū)分布在裂紋初始尖端及裂紋擴(kuò)展轉(zhuǎn)向區(qū)。

4 結(jié)論

該文以斷裂力學(xué)的巖石強(qiáng)度理論為指導(dǎo)，以數(shù)值模擬試驗(yàn)為主要手段，對(duì)內(nèi)含裂隙恐龍化石在壓力作用下?lián)p傷進(jìn)行了模擬分析，可得到如下結(jié)論：

(1)通過對(duì)內(nèi)含裂隙恐龍化石進(jìn)行數(shù)值壓縮試驗(yàn)得出：隨著裂隙角度增大，開裂角逐漸減小。裂隙角度的變化對(duì)裂紋擴(kuò)展步長(zhǎng)的影響不大，在裂隙角度大于30°后，裂紋擴(kuò)展步長(zhǎng)減小；裂隙角度大于45°之后，裂紋擴(kuò)展步長(zhǎng)略有增大。裂隙角度在15°～25°時(shí)，極限載荷呈緩降趨勢(shì)；裂隙角度在25°～45°時(shí)，極限載荷逐漸增大；裂隙角度在45°～75°時(shí)，極限載荷逐漸減小。

(2)裂隙處特別是兩端應(yīng)力集中，其他部位應(yīng)力較低、分布均勻。

(3)通過對(duì)裂隙角度為30°的恐龍化石進(jìn)行裂紋擴(kuò)展路徑及應(yīng)力場(chǎng)分析表明，裂紋擴(kuò)展是從裂紋尖端起裂，預(yù)制裂紋端部產(chǎn)生微小翼型裂紋，繼續(xù)施加荷載時(shí)翼型裂紋變寬延伸，隨后轉(zhuǎn)向上擴(kuò)展，最終裂紋擴(kuò)展到邊緣，產(chǎn)生的次生裂紋明顯呈鋸齒狀，最大應(yīng)力區(qū)分布在裂紋初始尖端及裂紋擴(kuò)展轉(zhuǎn)向區(qū)。

(4)以上2個(gè)試驗(yàn)表明裂隙角是引起恐龍化石破壞的主要因素之一，該試驗(yàn)為恐龍化石風(fēng)化機(jī)理研究提供了重要基礎(chǔ)數(shù)據(jù)。