本發(fā)明涉及工程巖體細(xì)觀時(shí)效破裂分析技術(shù)領(lǐng)域，具體地指巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型的構(gòu)建方法。

背景技術(shù)：

深部巖體工程開挖后的失穩(wěn)與破壞往往不是在開挖后立刻發(fā)生的，一般都存在著明顯的變形破裂時(shí)效性和災(zāi)變(巖爆、大變形等)的滯后性，嚴(yán)重危害工程的施工安全與長期運(yùn)營。目前，在細(xì)觀方面的時(shí)效力學(xué)研究成果相對(duì)較少?！渡盥翊罄韼r破裂擴(kuò)展時(shí)間效應(yīng)的顆粒流模擬》一文對(duì)錦屏大理巖破裂的時(shí)間效應(yīng)進(jìn)行了試驗(yàn)和二維數(shù)值分析(巖石力學(xué)與工程學(xué)報(bào)，2011，Vol.30 No.10：1989-1996)；《錦屏大理巖蠕變損傷演化細(xì)觀力學(xué)特征的數(shù)值模擬研究》一文應(yīng)用二維蠕變細(xì)觀力學(xué)模型對(duì)錦屏大理巖短期和長期強(qiáng)度特征進(jìn)行了數(shù)值研究(巖土力學(xué)，2013，Vol.34 No.12：3601-3608)。這類模型是以指數(shù)型構(gòu)建驅(qū)動(dòng)應(yīng)力和裂紋擴(kuò)展速度間的關(guān)系，用來描述巖石細(xì)觀層面上的時(shí)效破裂，適用于應(yīng)力和裂紋擴(kuò)展速度之間符合指數(shù)表達(dá)方式的這類巖體。另外，這類模型還存在如下不足之處：(1)顆粒間的剪切破裂準(zhǔn)則是一條與平行粘結(jié)正應(yīng)力平行的水平直線，也即這種剪切破裂準(zhǔn)則與平行粘結(jié)正應(yīng)力狀態(tài)無關(guān)，只要平行粘結(jié)剪切應(yīng)力大于或等于固定平行粘結(jié)剪切破裂強(qiáng)度，顆粒間即可發(fā)生剪切破裂，無法體現(xiàn)巖體中不同平行粘結(jié)正應(yīng)力具有不同平行粘結(jié)剪切破裂強(qiáng)度的客觀事實(shí)；(2)沒有考慮粘結(jié)力矩的差異作用對(duì)接觸破壞的影響，將粘結(jié)力矩的貢獻(xiàn)度對(duì)不同巖性的影響均視為一致；(3)對(duì)于應(yīng)力和裂紋擴(kuò)展速度之間不合符指數(shù)表達(dá)方式的巖體，這類模型缺乏適應(yīng)性。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的在于針對(duì)上述缺陷，提出了一種巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型的構(gòu)建方法，該模型結(jié)構(gòu)包括考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒粘結(jié)應(yīng)力二維模式、考慮彎矩貢獻(xiàn)因子的細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式、考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則和考慮阻尼效應(yīng)的細(xì)觀顆粒線性接觸二維模型。本發(fā)明適應(yīng)于應(yīng)力和裂紋擴(kuò)展速度之間的關(guān)系符合冪函數(shù)型的這類巖體，對(duì)于平面狀態(tài)下這類深部巖體工程的圍巖長期穩(wěn)定性預(yù)測、評(píng)價(jià)以及優(yōu)化設(shè)計(jì)提供技術(shù)支持。

本發(fā)明的目的是通過如下措施來達(dá)到的：巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型的構(gòu)建方法，其特殊之處在于，包括如下步驟：

步驟1：設(shè)定巖體內(nèi)部細(xì)觀顆粒粘結(jié)接觸的幾何參數(shù)量，包括顆粒粘結(jié)面積和顆粒粘結(jié)慣性矩，R^a、R^b分別為二維粘結(jié)接觸a端的顆粒半徑、b端的顆粒半徑，為巖體細(xì)觀顆粒粘結(jié)直徑乘數(shù)或半徑乘數(shù)，在二維情況下，粘結(jié)單位厚度為1時(shí)的顆粒粘結(jié)面積A和粘結(jié)慣性矩I分別通過公式(2)、公式(3)來確定：

其中：為巖體細(xì)觀顆粒二維粘結(jié)半徑，為巖體細(xì)觀顆粒二維粘結(jié)直徑乘數(shù)或半徑乘數(shù)，A為巖體細(xì)觀顆粒二維粘結(jié)面積，I為巖體細(xì)觀顆粒二維粘結(jié)慣性矩；

步驟201：利用巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的初始時(shí)間步長增量Δt，通過冪函數(shù)形式計(jì)算巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的直徑由公式(4)來確定；

其中：為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)法向應(yīng)力，為判斷巖體細(xì)觀顆粒二維粘結(jié)開始時(shí)效劣化衰減時(shí)的應(yīng)力閥值，為巖體細(xì)觀顆粒二維粘結(jié)的拉伸強(qiáng)度，為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)應(yīng)力比，β₁、β₂分別為巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化的第一控制參數(shù)、第二控制參數(shù)，為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的直徑，為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的直徑，Δt為巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的時(shí)間增量；

步驟202：根據(jù)步驟201中的公式(4)，設(shè)定細(xì)觀顆粒二維粘結(jié)直徑的時(shí)效衰減因子β，見公式(5)：

其中：A′、I′、分別表示為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)，A、I、為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)；

步驟203：根據(jù)步驟1中的公式(1)～公式(3)以及步驟202中的公式(5)，設(shè)定巖體細(xì)觀顆粒二維粘結(jié)幾何參數(shù)時(shí)效劣化衰減模式，在平面二維情況下，巖體細(xì)觀顆粒二維粘結(jié)直徑隨著時(shí)間增加而不斷劣化衰減，粘結(jié)單位厚度為1時(shí)的顆粒粘結(jié)的面積和慣性矩也隨著時(shí)間增加而不斷劣化衰減，分別見公式(6)、公式(7)；

其中：A′、I′分別表示為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩，A、I為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的粘結(jié)面積、粘結(jié)慣性矩；

步驟204：依次計(jì)算第j個(gè)至第k個(gè)的巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的二維粘結(jié)法向彎矩增量，在二維情況下，由粘結(jié)兩端顆粒的速度、角速度和給定的循環(huán)計(jì)算時(shí)間步長增量Δt，通過公式(8)、公式(9)、公式(10)確定第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)相對(duì)轉(zhuǎn)角巖體細(xì)觀顆粒二維粘結(jié)法向增量位移以及巖體細(xì)觀顆粒二維粘結(jié)切向增量位移再結(jié)合步驟203中的公式(7)和步驟202中的公式(5)，確定第i個(gè)巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的二維粘結(jié)彎矩增量，具體見公式(11)；

其中：ff、j、k、i是自然數(shù)，且2≤j≤ff≤k，j為每次循環(huán)計(jì)算中，包含時(shí)間效應(yīng)的巖體細(xì)觀顆粒二維粘結(jié)衰減后未破裂的初始標(biāo)號(hào)值，ff為中間標(biāo)號(hào)值，k為每次循環(huán)計(jì)算中，包含時(shí)間效應(yīng)的巖體細(xì)觀顆粒二維粘結(jié)衰減后未破裂的最末標(biāo)號(hào)值，i為第一個(gè)至最后一個(gè)二維粘結(jié)顆粒標(biāo)號(hào)值，分別為第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)接觸的a端和b端的絕對(duì)運(yùn)動(dòng)速度和角速度，n_n、n_s分別為巖體細(xì)觀顆粒二維粘結(jié)接觸面的法向單位向量和切向單位向量，分別為巖體細(xì)觀顆粒二維粘結(jié)法向位移增量和切向位移增量，為巖體細(xì)觀顆粒二維粘結(jié)法向剛度，為巖體細(xì)觀顆粒二維粘結(jié)彎矩增量。

步驟205：根據(jù)步驟203中的公式(6)和公式(7)、步驟204中的公式(8)、公式(9)和公式(11)以及步驟202中的公式(5)，依次更新計(jì)算第j個(gè)至第k個(gè)巖體細(xì)觀顆粒粘結(jié)未破裂的粘結(jié)接觸并包含時(shí)間效應(yīng)的二維粘結(jié)法向力、切向力和切向彎矩，通過公式(12)、公式(13)、公式(14)計(jì)算第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)接觸的粘結(jié)法向力、切向力和切向彎矩，在二維情況下，通過公式(15)來確定巖體細(xì)觀顆粒粘結(jié)法向彎矩，

法向力：

切向力：

切向彎矩：

法向彎矩：

其中：分別為第i個(gè)巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的粘結(jié)法向力、粘結(jié)切向力、包含時(shí)間效應(yīng)的粘結(jié)法向彎矩、粘結(jié)切向彎矩、粘結(jié)法向位移增量和粘結(jié)切向位移增量，為巖體細(xì)觀顆粒二維粘結(jié)切向剛度，+＝為加法自反運(yùn)算符，-＝為減法自反運(yùn)算符。

步驟206：設(shè)置彎矩貢獻(xiàn)因子考慮彎矩對(duì)巖體細(xì)觀顆粒粘結(jié)法向應(yīng)力的貢獻(xiàn)程度，根據(jù)巖體細(xì)觀顆粒粘結(jié)二維正應(yīng)力計(jì)算公式和二維粘結(jié)剪應(yīng)力計(jì)算公式同時(shí)將這兩個(gè)公式中A、I以及用A′、I′及替換，然后將步驟203中的公式(6)和公式(7)以及步驟202中的公式(5)代入，獲得包含冪函數(shù)型時(shí)間效應(yīng)和細(xì)觀顆粒粘結(jié)考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力計(jì)算公式和二維粘結(jié)剪應(yīng)力計(jì)算公式，分別見公式(16)和公式(17)，

步驟207：將步驟206中包含時(shí)間效應(yīng)的代入公式(18)，確定考慮彎矩貢獻(xiàn)因子且?guī)Ю旖刂瓜薜哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則，并且依次計(jì)算第j個(gè)至第k個(gè)的巖體細(xì)觀顆粒二維粘結(jié)應(yīng)力，用于判斷巖體細(xì)觀顆粒粘結(jié)是否破裂以及破裂模式，在該準(zhǔn)則的巖體細(xì)觀顆粒粘結(jié)應(yīng)力中包含了冪函數(shù)型時(shí)間效應(yīng)和考慮彎矩貢獻(xiàn)因子，

其中：f^s、fⁿ分別為巖體細(xì)觀顆粒二維粘結(jié)的時(shí)效剪切破裂準(zhǔn)則、時(shí)效拉伸破裂準(zhǔn)則，為第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)的二維粘結(jié)剪應(yīng)力，為第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)且考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力，分別為巖體細(xì)觀顆粒二維粘結(jié)的拉伸強(qiáng)度、抗剪強(qiáng)度，為巖體細(xì)觀顆粒二維粘結(jié)的粘聚力，為巖體細(xì)觀顆粒二維粘結(jié)的內(nèi)摩擦角；f^s大于等于0表示巖體細(xì)觀顆粒粘結(jié)剪切破裂，小于0表示巖體細(xì)觀顆粒粘結(jié)未發(fā)生剪切破裂；fⁿ大于等于0表示巖體細(xì)觀顆粒粘結(jié)拉伸破裂，小于0表示巖體細(xì)觀顆粒粘結(jié)未發(fā)生拉伸破裂；

步驟208：當(dāng)步驟207中的公式(18)中的f^s或fⁿ大于等于0時(shí)，表明巖體細(xì)觀顆粒粘結(jié)發(fā)生破裂，此時(shí)巖體細(xì)觀顆粒的運(yùn)動(dòng)模式采用考慮阻尼效應(yīng)的二維線性接觸模型來表達(dá)；當(dāng)步驟207中的公式(18)中的f^s和fⁿ都小于0時(shí)，表明巖體細(xì)觀顆粒粘結(jié)未破裂，繼續(xù)循環(huán)步驟201至207，計(jì)算、更新、判斷巖體細(xì)觀顆粒接觸的粘結(jié)狀態(tài)，直至巖體不產(chǎn)生新的細(xì)觀顆粒粘結(jié)破裂或者巖體細(xì)觀顆粒粘結(jié)破裂加速發(fā)展而形成宏觀破壞，循環(huán)終止。

優(yōu)選地，所述巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的初始時(shí)間步長增量Δt的確定，是采用考慮彎矩貢獻(xiàn)因子的細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式，由每個(gè)時(shí)間步內(nèi)的細(xì)觀顆粒二維粘結(jié)首次衰減破裂所損耗的時(shí)間來確定，即通過第一個(gè)細(xì)觀顆粒粘結(jié)按冪函數(shù)型模式進(jìn)行衰減破裂所歷時(shí)的時(shí)間除以直至第一個(gè)細(xì)觀顆粒粘結(jié)破裂所需要的計(jì)算循環(huán)次數(shù)來估算初始時(shí)間步長增量Δt，見公式其中，為第i個(gè)接觸的巖體細(xì)觀顆粒二維粘結(jié)直徑乘數(shù)，n_c為第一個(gè)巖體細(xì)觀顆粒二維粘結(jié)破裂所需的循環(huán)計(jì)算的次數(shù)，β_σ、β_τ分別為巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子、二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子，i依次為第一個(gè)至最后一個(gè)巖體細(xì)觀顆粒粘結(jié)數(shù)，∞為無窮大。

優(yōu)選地，所述巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_σ和巖體細(xì)觀顆粒二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_τ的確定包括如下步驟，其中這些步驟中包含的公式下標(biāo)1代表第一個(gè)按冪函數(shù)型模式進(jìn)行時(shí)效衰減劣化的細(xì)觀顆粒二維粘結(jié)破裂標(biāo)號(hào)：

步驟211：在二維情況下，由巖體細(xì)觀顆粒粘結(jié)兩端顆粒的速度、角速度和給定的循環(huán)計(jì)算時(shí)間步長增量Δt_c，通過公式確定細(xì)觀顆粒粘結(jié)接觸的相對(duì)轉(zhuǎn)角通過公式確定細(xì)觀顆粒粘結(jié)法向增量位移通過公式確定顆粒粘結(jié)切向增量位移通過公式確定顆粒粘結(jié)接觸的彎矩增量；

步驟212：根據(jù)步驟211中的公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)法向力；根據(jù)步驟211中的公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)切向力；根據(jù)步驟211中的公式和公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)切向彎矩；通過公式確定細(xì)觀顆粒接觸的粘結(jié)法向彎矩，其中，+＝為加法自反運(yùn)算符，-＝為減法自反運(yùn)算符；

步驟213：在二維情況下，通過公式確定細(xì)觀顆粒接觸的粘結(jié)正應(yīng)力，通過公式確定細(xì)觀顆粒接觸的粘結(jié)剪應(yīng)力，將這兩個(gè)公式中A、I以及用A′、I′及替換，然后將步驟203中的公式(6)和公式(7)以及步驟202中的公式(5)代入，獲得細(xì)觀顆粒粘結(jié)包含冪函數(shù)型時(shí)間效應(yīng)和考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力計(jì)算公式和包含冪函數(shù)型時(shí)間效應(yīng)的二維粘結(jié)剪應(yīng)力計(jì)算公式

步驟214：將代入公式并令β＝β_σ；將代入公式并令β＝β_τ，據(jù)此，分別得到所述巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子以及巖體細(xì)觀顆粒二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子

優(yōu)選地，所述巖體細(xì)觀顆粒粘結(jié)發(fā)生破裂后，巖體細(xì)觀顆粒的運(yùn)動(dòng)模式采用考慮阻尼效應(yīng)的二維線性接觸模型來表達(dá)，用于描述巖體細(xì)觀顆粒粘結(jié)時(shí)效破裂后細(xì)觀顆粒的應(yīng)力、變形及運(yùn)行規(guī)律，考慮阻尼效應(yīng)的二維線性接觸模型的構(gòu)建包括如下步驟：

步驟301：通過Monte Carlo搜索算法，遍歷尋找?guī)r體細(xì)觀顆粒每個(gè)二維線性接觸端a、二維線性接觸端b(顆粒與顆粒、顆粒與墻)的中心坐標(biāo)，在二維情況下，通過公式(19)計(jì)算接觸點(diǎn)a端，接觸點(diǎn)b端的中心距離：

其中：d為二維線性接觸兩端顆粒與顆?；蛘哳w粒與墻之間的中心距離，為二維線性接觸端a的坐標(biāo)，為二維線性接觸端b的坐標(biāo)。

步驟302：二維平面狀態(tài)下巖體細(xì)觀顆粒間每個(gè)接觸點(diǎn)的單位向量通過公式(20)計(jì)算，如果是顆粒與顆粒之間的接觸，利用步驟301中得到二維線性接觸兩端的中心點(diǎn)坐標(biāo)及距離計(jì)算；如果是顆粒與墻接觸，直接利用墻體的法向量等效替換來計(jì)算；確定每個(gè)接觸端的單位向量：

其中：n_i為接觸的單位矢量，為接觸端b的方向向量，為接觸端a的方向向量，n_wall為約束墻的方向向量；

步驟303：巖體細(xì)觀顆粒粘結(jié)破裂后，每一個(gè)二維線性接觸點(diǎn)的接觸重疊量U，通過步驟301計(jì)算顆粒間距d以及二維線性接觸兩端的顆粒半徑R^a、R^b，再利用公式(21)來確定；通過設(shè)定顆粒二維線性接觸參考距離g_r，并結(jié)合公式(22)，確定顆粒二維線性接觸的距離g_s：

g_s＝|U|-g_r (22)

步驟304：確定巖體細(xì)觀顆粒接觸點(diǎn)法向、切向等效剛度，利用接觸兩端顆粒實(shí)體或是墻體的剛度k^a，k^b等效代替為接觸點(diǎn)的剛度，按公式(23)計(jì)算：

其中：K_n、K_s為巖體細(xì)觀顆粒接觸點(diǎn)等效的法向剛度和切向剛度，為顆粒與顆粒或者顆粒與墻的接觸a端的法向剛度和切向剛度，為顆粒與顆?；蛘哳w粒與墻的接觸b端的法向剛度和切向剛度。

步驟305：確定巖體中接觸兩端顆粒間的相對(duì)運(yùn)動(dòng)速度，利用公式(24)、公式(25)來計(jì)算。其中e_pqz為Ricci指標(biāo)置換符號(hào)，按照公式(26)來計(jì)算：

其中：V_p與V_q等價(jià)，V_p與V_q為巖體中接觸兩端顆粒間的相對(duì)運(yùn)動(dòng)速度，p、q為指標(biāo)等價(jià)符號(hào)，p＝1，q＝1表示顆粒與顆粒接觸，p＝2，q＝2時(shí)表示顆粒與墻接觸，為顆粒與顆?；蚴穷w粒與墻的接觸b端單元的速度，為顆粒與顆粒或是顆粒與墻的接觸a端單元的速度，為顆粒與顆粒或是顆粒與墻的接觸a端單元的角速度，為顆粒與顆?；蚴穷w粒與墻的接觸b端單元的角速度，為顆粒與顆粒或是顆粒與墻的接觸a端的位移，為顆粒與顆?；蚴穷w粒與墻的接觸b端的位移，為位移指標(biāo)變換的中間過渡符號(hào)，表示指標(biāo)符號(hào)為p時(shí)顆粒-顆?；蛘哳w粒-墻的接觸a端的速度，表示指標(biāo)符號(hào)為q時(shí)顆粒-顆?；蛘哳w粒-墻的接觸a端的速度，表示指標(biāo)符號(hào)為p時(shí)顆粒-顆?；蛘哳w粒-墻的接觸b端的速度，表示指標(biāo)符號(hào)為q時(shí)顆粒-顆?；蛘哳w粒-墻的接觸b端的速度(只有a端和b端兩個(gè)接觸端)。

步驟306：對(duì)于巖體細(xì)觀顆粒線性接觸模型的初始時(shí)間步長增量Δt的取值，通過公式(29)估計(jì)最小的時(shí)間步長Δt，確保所構(gòu)建模型的計(jì)算時(shí)間步長小于該值，即可保證系統(tǒng)積分計(jì)算趨于穩(wěn)定；通過公式(30)、公式(31)、公式(32)確定每個(gè)線性接觸的總位移增量、法向位移增量和切向位移增量：

R＝min(R_a，R_b) (27)

ΔU_p1＝V_p1Δt (30)

其中：R為巖體細(xì)觀顆粒的等效半徑，m為巖體細(xì)觀顆粒質(zhì)量，J1為巖體細(xì)觀顆粒的轉(zhuǎn)動(dòng)慣量；k^平為巖體細(xì)觀顆粒系統(tǒng)平移剛度，k^轉(zhuǎn)為巖體細(xì)觀顆粒系統(tǒng)旋轉(zhuǎn)剛度；ΔU_p1為巖體細(xì)觀顆粒二維線性接觸的總位移增量，Δδ_n、物理意義相同，均表示巖體細(xì)觀顆粒二維線性接觸的法向位移增量，Δδ_s、物理意義相同，均表示巖體細(xì)觀顆粒二維線性接觸的切向位移增量，V_p1與V_q1為巖體細(xì)觀顆粒接觸兩端的相對(duì)運(yùn)動(dòng)速度，n為單位法向量，p1，q1為張量指標(biāo)變換符號(hào)。

步驟307：由公式(22)判定巖體細(xì)觀顆粒表面接觸允許存在的最大距離，計(jì)算法向和切向位移更新因子，另外，巖體細(xì)觀顆粒二維線性接觸法向位移增量的更新是采用前一步的法向位移增量與更新因子α的乘積獲得，巖體細(xì)觀顆粒二維線性接觸切向位移增量的更新是采用前一步的切向位移增量與更新因子α的乘積獲得。

其中：(g_s)₀為模型計(jì)算初始時(shí)刻的表面接觸距離，g_s為巖體細(xì)觀顆粒接觸的距離，α為位移更新因子。

步驟308：巖體細(xì)觀顆粒二維線性接觸的法向線性力采取相對(duì)矢量累加(M_l＝1)和絕對(duì)矢量累加(M_l＝0)模式，通過公式(33)、(34)計(jì)算獲得；巖體細(xì)觀顆粒二維線性接觸的切向線性力采用巖體細(xì)觀顆粒接觸滑動(dòng)來表示，通過公式(35)計(jì)算獲得：

其中：分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向線性力、切向線性力，k_n、k_s分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向、切向線性剛度，Δδ_n、Δδ_s分別為巖體細(xì)觀顆粒二維線性接觸的法向位移增量、切向位移增量，分別為巖體細(xì)觀顆粒二維線性接觸的初始法向力增量值、切向力增量值，為巖體細(xì)觀顆粒未滑動(dòng)時(shí)的靜摩擦力，為巖體細(xì)觀顆?；瑒?dòng)摩擦力，其值可通過摩擦系數(shù)μ與乘積得到。

步驟309：巖體細(xì)觀顆粒線性接觸的法向阻尼采用全法向模式M_d＝{0，2}和無拉伸模式M_d＝{1，3}兩種，通過公式(36)計(jì)算，其中m_c為等效顆粒質(zhì)量，按公式(37)計(jì)算，巖體細(xì)觀顆粒線性接觸的切向阻尼采用全剪切模式M_d＝{0，1}和滑-剪模式M_d＝{2，3}，按照公式(38)來計(jì)算，

其中：分別為巖體細(xì)觀顆粒線性接觸的法向線性阻尼力、切向線性阻尼力，β_n為巖體細(xì)觀顆粒線性接觸的法向阻尼系數(shù)，β_s為巖體細(xì)觀顆粒線性接觸的切向阻尼系數(shù)，k_n為巖體細(xì)觀顆粒線性接觸的法向線性剛度，k_s為巖體細(xì)觀顆粒線性接觸的切向線性剛度，分別為巖體細(xì)觀顆粒線性接觸的法向速率和切向速率，F(xiàn)^*為巖體細(xì)觀顆粒線性接觸的全法向阻尼力，表達(dá)式為m_c為巖體細(xì)觀等效顆粒質(zhì)量，m⁽¹⁾為巖體細(xì)觀顆粒接觸端1的細(xì)觀顆粒質(zhì)量，m⁽²⁾為巖體細(xì)觀顆粒接觸端2的細(xì)觀顆粒質(zhì)量，F(xiàn)^d為巖體細(xì)觀顆粒線性接觸的總阻尼力。

本發(fā)明所提出的一種巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型及構(gòu)建方法，其有益效果和優(yōu)勢主要體現(xiàn)在：

(1)本發(fā)明中巖體細(xì)觀顆粒粘結(jié)正應(yīng)力二維計(jì)算公式中設(shè)置了彎矩貢獻(xiàn)因子，不僅考慮了彎矩對(duì)細(xì)觀顆粒粘結(jié)正應(yīng)力的貢獻(xiàn)程度，而且還考慮了彎矩對(duì)巖體長期強(qiáng)度的影響，適合描述平面應(yīng)力或平面應(yīng)變條件下巖體的細(xì)觀力學(xué)破裂行為。

(2)本發(fā)明中構(gòu)建了考慮彎矩貢獻(xiàn)因子的細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式，包括在巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減時(shí)，設(shè)置了冪函數(shù)型與考慮彎矩貢獻(xiàn)因子的粘結(jié)應(yīng)力相關(guān)的細(xì)觀顆粒粘結(jié)二維劣化衰減模式，設(shè)置了細(xì)觀顆粒粘結(jié)隨時(shí)間逐步劣化衰減的二維模式，設(shè)置了細(xì)觀顆粒粘結(jié)的面積和截面慣性矩時(shí)效劣化衰減二維模式；同時(shí)按照這種時(shí)效劣化衰減模式估算巖體細(xì)觀顆粒粘結(jié)破裂的初始時(shí)間步長。這種冪函數(shù)型構(gòu)建模式適合描述平面狀態(tài)下一類深部巖體的細(xì)觀力學(xué)時(shí)效破裂機(jī)制和響應(yīng)規(guī)律。

(3)本發(fā)明中在所構(gòu)建的二維細(xì)觀時(shí)效破裂冪函數(shù)型模型中，嵌入考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜牡哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則。在巖體細(xì)觀顆粒粘結(jié)時(shí)效破裂時(shí)，采用所嵌入的考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜牡哪?庫倫時(shí)效破裂準(zhǔn)則來判斷；在該準(zhǔn)則顆粒粘結(jié)應(yīng)力中包含冪函數(shù)型時(shí)間效應(yīng)且增加了彎矩貢獻(xiàn)因子，不僅可以描述與顆粒粘結(jié)正應(yīng)力相關(guān)時(shí)效剪切破裂強(qiáng)度的差異，還可以對(duì)時(shí)效拉伸破裂進(jìn)行合理的表達(dá)，且考慮了彎矩對(duì)粘結(jié)時(shí)效破裂的影響，符合平面條件下一類巖體細(xì)觀時(shí)效破裂模式。

(4)本發(fā)明中在所構(gòu)建的二維細(xì)觀時(shí)效破裂冪函數(shù)型模型中，嵌入考慮阻尼效應(yīng)的細(xì)觀顆粒線性接觸二維模型結(jié)構(gòu)，在巖體時(shí)效破裂后，通過指定二維接觸參考距離設(shè)定巖體細(xì)觀顆粒間接觸距離，設(shè)置考慮巖體細(xì)觀顆粒間受力變形的二維接觸模式以及在巖體細(xì)觀顆粒之間設(shè)置考慮二維滑動(dòng)摩擦的作用模式，同時(shí)設(shè)置二維接觸的阻尼模式，可合理描述平面條件下一類深部工程巖體在時(shí)效破裂后的顆粒運(yùn)動(dòng)和受力特征。

附圖說明

圖1是本發(fā)明模型中細(xì)觀顆粒與顆粒接觸示意圖。

圖2是本發(fā)明模型中細(xì)觀顆粒與剛性墻接觸示意圖。

圖3是本發(fā)明模型中細(xì)觀顆粒重疊狀態(tài)示意圖。

圖4是本發(fā)明模型中細(xì)觀顆粒剛度計(jì)算示意圖。

圖5是本發(fā)明模型中細(xì)觀顆粒粘結(jié)線性切向力和切向位移示意圖。

圖6是本發(fā)明模型中細(xì)觀顆粒粘接觸本構(gòu)物理模型示意圖。

圖7是本發(fā)明模型中細(xì)觀顆粒線性粘結(jié)結(jié)構(gòu)示意圖。

圖8是本發(fā)明模型中考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜牡哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則示意圖。

圖9是本發(fā)明模型中細(xì)觀顆粒粘結(jié)直徑(或半徑)時(shí)效劣化衰減示意圖。

圖10是本發(fā)明模型中細(xì)觀顆粒二維接觸面的法向和切向單位向量示意圖。

圖11是本發(fā)明模型的構(gòu)建方法的流程示意圖。

圖12是本發(fā)明模型基礎(chǔ)模型試樣

圖13是本發(fā)明模型蠕變位移與時(shí)間關(guān)系曲線圖

圖中：1代表兩顆粒的中心距離d，2代表巖體細(xì)觀顆粒間的半接觸距離，3代表巖體細(xì)觀顆粒間的半?yún)⒖季嚯xg_r，4代表巖體細(xì)觀顆粒a的坐標(biāo)，5代表巖體細(xì)觀顆粒b的坐標(biāo)，6巖體細(xì)觀顆粒表面接觸距離的中心坐標(biāo)，7代表巖體細(xì)觀顆粒表面接觸距離g_s，8代表巖體細(xì)觀顆粒間的單位接觸法向量，9代表巖體細(xì)觀顆粒a的半徑R^a，10代表巖體細(xì)觀顆粒b的半徑R^b，11代表巖體細(xì)觀顆粒接觸點(diǎn)的接觸重疊量U，12代表b(巖體顆粒或是邊界墻體)的剛度(法向、切向剛度統(tǒng)稱)k^b，13代表a(巖體顆粒或是邊界墻體)的剛度(法向、切向剛度統(tǒng)稱)k^a，14代表巖體細(xì)觀顆粒接觸點(diǎn)的等效剛度，15代表總位移增量ΔU_i，16代表初始法向力增量值，17代表初始接觸力矢量和，18代表初始切向力增量值，19代表所構(gòu)建的二維時(shí)效破裂模型法向位移增量Δδ_n或20代表所構(gòu)建的二維時(shí)效破裂模型切向位移增量Δδ_s或21代表巖體細(xì)觀顆粒粘結(jié)的拉伸強(qiáng)度22代表巖體細(xì)觀顆粒粘結(jié)法向剛度23代表巖體細(xì)觀顆粒接觸點(diǎn)的法向剛度K_n，24代表巖體細(xì)觀顆粒粘結(jié)切向剛度25代表巖體細(xì)觀顆粒粘結(jié)剪切強(qiáng)度，25.1代表為巖體細(xì)觀顆粒粘結(jié)的粘聚力，25.2代表巖體細(xì)觀顆粒粘結(jié)內(nèi)摩擦角26代表巖體細(xì)觀顆粒接觸點(diǎn)的切向剛度K_s，27代表巖體細(xì)觀顆粒線性接觸滑動(dòng)摩擦系數(shù)，28代表巖體細(xì)觀顆粒線性接觸法向阻尼系數(shù)β_n，29代表巖體細(xì)觀顆粒線性接觸切向阻尼系數(shù)β_s，30代表為巖體細(xì)觀顆粒粘結(jié)半徑乘數(shù)31代表巖體細(xì)觀顆粒粘結(jié)直徑32代表考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜哪?庫倫時(shí)效破裂準(zhǔn)則，33代表第i個(gè)接觸的包含時(shí)間效應(yīng)的粘結(jié)剪應(yīng)力34代表第i個(gè)接觸的包含時(shí)間效應(yīng)且考慮彎矩貢獻(xiàn)因子的粘結(jié)正應(yīng)力35代表巖體細(xì)觀顆粒粘結(jié)時(shí)效衰減的半徑36代表巖體細(xì)觀顆粒粘結(jié)時(shí)效衰減的直徑37代表巖體細(xì)觀顆粒粘結(jié)未衰減時(shí)的直徑38代表巖體細(xì)觀顆粒粘結(jié)未衰減時(shí)的半徑39代表巖體細(xì)觀顆粒接觸面的法向向量n_n，40代表巖體細(xì)觀顆粒接觸面切向單位向量n_s。

具體實(shí)施方式

下面結(jié)合附圖和具體構(gòu)建步驟及實(shí)施實(shí)例，對(duì)本發(fā)明模型進(jìn)行詳細(xì)的闡述。實(shí)例說明僅是輔助用于本發(fā)明的理解，而不限定本發(fā)明的實(shí)際應(yīng)用范圍。在閱讀了本發(fā)明以后，本領(lǐng)域的技術(shù)人員對(duì)本發(fā)明的各種等價(jià)形式的修改均屬于本發(fā)明所申請(qǐng)的權(quán)利要求所限定的范圍。

注：說明書中的所有的標(biāo)號(hào)前面寫明了公式，如公式(1)，均為公式標(biāo)號(hào)。

如圖1～圖10所示，本發(fā)明所述的巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型，適應(yīng)于二維顆粒離散元、二維顆粒不連續(xù)變形分析方法、二維顆粒流形元；冪函數(shù)型模型包括考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒粘結(jié)應(yīng)力二維模式、考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式、考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則以及考慮阻尼效應(yīng)的細(xì)觀顆粒線性接觸二維模型。

考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒粘結(jié)應(yīng)力二維模式是指巖體細(xì)觀顆粒粘結(jié)正應(yīng)力二維計(jì)算公式中設(shè)置了彎矩貢獻(xiàn)因子考慮彎矩對(duì)巖體細(xì)觀顆粒二維粘結(jié)正應(yīng)力的貢獻(xiàn)程度；在上述公式中，為第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)正應(yīng)力，分別為第i個(gè)接觸的巖體細(xì)觀顆粒二維粘結(jié)法向力、切向彎矩，為巖體細(xì)觀顆粒二維粘結(jié)半徑，為彎矩貢獻(xiàn)因子，I為巖體細(xì)觀顆粒二維粘結(jié)的慣性矩，A為巖體細(xì)觀顆粒二維粘結(jié)面積，i依次為第一個(gè)至最后一個(gè)巖體細(xì)觀顆粒粘結(jié)數(shù)。第i個(gè)接觸的巖體細(xì)觀顆粒二維粘結(jié)法向力粘結(jié)切向彎矩的計(jì)算方法為：法向力式中，為巖體細(xì)觀顆粒二維粘結(jié)法向位移增量，為巖體細(xì)觀顆粒二維粘結(jié)法向剛度，切向彎矩式中，為巖體細(xì)觀顆粒二維粘結(jié)切向相對(duì)轉(zhuǎn)角增量，+＝為加法自反運(yùn)算符，-＝為減法自反運(yùn)算符，法向彎矩

考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式包括在巖體細(xì)觀顆粒二維粘結(jié)時(shí)效劣化衰減時(shí)，設(shè)置了與考慮彎矩貢獻(xiàn)因子的粘結(jié)應(yīng)力相關(guān)的冪函數(shù)型模式，這種冪函數(shù)型模式中的細(xì)觀顆粒二維粘結(jié)直徑隨時(shí)間逐步劣化衰減，見粘結(jié)直徑公式，

式中，為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的直徑，為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的直徑，為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)法向應(yīng)力，為判斷巖體細(xì)觀顆粒二維粘結(jié)開始時(shí)效劣化衰減時(shí)的應(yīng)力閥值，為巖體細(xì)觀顆粒二維粘結(jié)的拉伸強(qiáng)度，為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)應(yīng)力比，β₁、β₂為控制巖體細(xì)觀顆粒二維粘結(jié)時(shí)效劣化的兩個(gè)指定指數(shù)，Δt為巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的時(shí)間增量；設(shè)置了巖體細(xì)觀顆粒粘結(jié)面積和面慣性矩時(shí)效劣化衰減二維模式，分別見粘結(jié)單位厚度為1時(shí)的巖體細(xì)觀顆粒隨時(shí)間劣化衰減的粘結(jié)面積計(jì)算公式粘結(jié)單位厚度為1時(shí)的慣性矩I′計(jì)算公式其中，β為巖體細(xì)觀顆粒二維粘結(jié)直徑的時(shí)效衰減因子，其計(jì)算公式見

其中，A′、I′、分別表示為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)，A、I、為巖體細(xì)觀顆粒粘結(jié)未衰減時(shí)的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)；同時(shí)按照這種冪函數(shù)型時(shí)效劣化衰減模式估算巖體細(xì)觀顆粒粘結(jié)破裂的初始時(shí)間步長Δt，見公式

其中，為第i個(gè)接觸的巖體細(xì)觀顆粒二維粘結(jié)直徑乘數(shù)，n_c為第一個(gè)巖體細(xì)觀顆粒二維粘結(jié)破裂所需的循環(huán)計(jì)算的次數(shù)，β_σ、β_τ分別為巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子、二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子，i依次為第一個(gè)至最后一個(gè)巖體細(xì)觀顆粒粘結(jié)數(shù)，∞為無窮大。

巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_σ和二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_τ的計(jì)算公式分別為

其中，分別為第i個(gè)顆粒接觸的粘結(jié)法向力、切向力和粘結(jié)切向彎矩，為巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度，為巖體細(xì)觀顆粒二維粘結(jié)的粘聚力，為巖體細(xì)觀顆粒二維粘結(jié)的內(nèi)摩擦角。

考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則是指在巖體細(xì)觀顆粒二維粘結(jié)時(shí)效破裂時(shí)，采用所嵌入的考慮彎矩貢獻(xiàn)效應(yīng)且?guī)Ю旖刂瓜薜哪?庫倫時(shí)效破裂準(zhǔn)則來判斷，見公式

其中，f^s、fⁿ分別為巖體細(xì)觀顆粒二維粘結(jié)的時(shí)效剪切破裂準(zhǔn)則、時(shí)效拉伸破裂準(zhǔn)則，分別為巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度、抗剪強(qiáng)度，分別為第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)且考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒二維粘結(jié)正應(yīng)力、剪應(yīng)力。

第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)且考慮彎矩貢獻(xiàn)因子的巖體細(xì)觀顆粒二維粘結(jié)正應(yīng)力的計(jì)算公式為

第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)的巖體細(xì)觀顆粒二維粘結(jié)剪應(yīng)力的計(jì)算公式為

在該準(zhǔn)則的二維粘結(jié)應(yīng)力中包含了冪函數(shù)型時(shí)間效應(yīng)，見巖體細(xì)觀顆粒二維粘結(jié)直徑的時(shí)效衰減因子β計(jì)算公式

β₁、β₂分別為控制巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化的第一控制參數(shù)、第二控制參數(shù)；f^s大于等于0表示巖體細(xì)觀顆粒二維粘結(jié)剪切破裂，小于0表示巖體細(xì)觀顆粒二維粘結(jié)未發(fā)生剪切破裂；fⁿ大于等于0表示巖體細(xì)觀顆粒二維粘結(jié)拉伸破裂，小于0表示巖體細(xì)觀顆粒二維粘結(jié)未發(fā)生拉伸破裂。

考慮阻尼效應(yīng)的細(xì)觀顆粒線性接觸二維模型是指在巖體細(xì)觀顆粒粘結(jié)時(shí)效破裂后，通過給定的二維線性接觸參考距離g_r設(shè)置了細(xì)觀顆粒二維線性接觸距離g_s，見巖體細(xì)觀顆粒二維線性接觸距計(jì)算公式

其中，為巖體內(nèi)部顆粒與顆粒二維線性接觸端a的坐標(biāo)，為巖體內(nèi)部顆粒與顆粒二維線性接觸端b的坐標(biāo)，R^a、R^b分別為巖體細(xì)觀二維線性接觸端a的顆粒半徑和二維線性接觸端b的顆粒半徑；設(shè)置了考慮巖體細(xì)觀顆粒間變形的二維線性接觸模式，在巖體細(xì)觀顆粒之間設(shè)置了考慮二維滑動(dòng)摩擦線力的作用模式，巖體細(xì)觀顆粒間受力變形的二維線性接觸法向線性力計(jì)算公式取M_l＝1為相對(duì)矢量累加模式，取M_l＝0為絕對(duì)矢量累加模式，巖體細(xì)觀顆粒間受力變形的二維線性接觸切向線性力計(jì)算公式為其中，分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向線性力、切向線性力，k_n、k_s分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向、切向線性剛度，Δδ_n、Δδ_s分別為法向位移增量、切向位移增量，分別為初始法向力增量值和切向力增量值，為顆粒未滑動(dòng)時(shí)的靜摩擦力，為巖體細(xì)觀顆?；瑒?dòng)摩擦力，通過摩擦系數(shù)μ與乘積得到；

同時(shí)設(shè)置二維線性接觸的阻尼模式，其中法向阻尼采用全法向模式M_d＝{0，2}和無拉伸模式M_d＝{1，3}兩種，通過公式計(jì)算，其中m_c為等效顆粒質(zhì)量，按公式計(jì)算，切向阻尼采用全剪切模式M_d＝{0，1}和滑-剪模式M_d＝{2，3}，按照公式來計(jì)算，其中：分別為法向阻尼力、切向阻尼力，β_n為法向阻尼系數(shù)、β_s為切向阻尼系數(shù)，k_n為法向線性剛度、k_s為切向線性剛度，為法向速率、切向速率，m_c為等效顆粒質(zhì)量。F^*為巖體細(xì)觀顆粒線性接觸的全法向阻尼力，表達(dá)式為m⁽¹⁾為二維線性接觸端1的顆粒質(zhì)量，m⁽²⁾為二維線性接觸端2的顆粒質(zhì)量。

如圖11所示，本發(fā)明所述的巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型的構(gòu)建方法，包括如下步驟：

步驟1：設(shè)定巖體內(nèi)部細(xì)觀顆粒粘結(jié)接觸的幾何參數(shù)量，包括顆粒粘結(jié)面積和顆粒粘結(jié)慣性矩，R^a、R^b分別為二維粘結(jié)接觸a端的顆粒半徑、b端的顆粒半徑，為巖體細(xì)觀顆粒粘結(jié)直徑乘數(shù)或半徑乘數(shù)，在二維情況下，粘結(jié)單位厚度為1時(shí)的顆粒粘結(jié)面積A和粘結(jié)慣性矩I分別通過公式(2)、公式(3)來確定：

其中：為巖體細(xì)觀顆粒二維粘結(jié)半徑，為巖體細(xì)觀顆粒二維粘結(jié)直徑乘數(shù)或半徑乘數(shù)，A為巖體細(xì)觀顆粒二維粘結(jié)面積，I為巖體細(xì)觀顆粒二維粘結(jié)慣性矩；

步驟201：利用巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的初始時(shí)間步長增量Δt，通過冪函數(shù)形式計(jì)算巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的直徑由公式(4)來確定；

其中：為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)法向應(yīng)力，為判斷巖體細(xì)觀顆粒二維粘結(jié)開始時(shí)效劣化衰減時(shí)的應(yīng)力閥值，為巖體細(xì)觀顆粒二維粘結(jié)的拉伸強(qiáng)度，為考慮彎矩貢獻(xiàn)因子的二維粘結(jié)應(yīng)力比，β₁、β₂分別為巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化的第一控制參數(shù)、第二控制參數(shù)，為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的直徑，為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的直徑，Δt為巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的時(shí)間增量；

步驟202：根據(jù)步驟201中的公式(4)，設(shè)定細(xì)觀顆粒二維粘結(jié)直徑的時(shí)效衰減因子β，見公式(5)：

其中：A′、I′、分別表示為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)，A、I、為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的粘結(jié)直徑、粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩、粘結(jié)直徑乘數(shù)；

步驟203：根據(jù)步驟1中的公式(1)～公式(3)以及步驟202中的公式(5)，設(shè)定巖體細(xì)觀顆粒二維粘結(jié)幾何參數(shù)時(shí)效劣化衰減模式，在平面二維情況下，巖體細(xì)觀顆粒二維粘結(jié)直徑隨著時(shí)間增加而不斷劣化衰減，粘結(jié)單位厚度為1時(shí)的顆粒粘結(jié)的面積和慣性矩也隨著時(shí)間增加而不斷劣化衰減，分別見公式(6)、公式(7)；

其中：A′、I′分別表示為巖體細(xì)觀顆粒二維粘結(jié)隨時(shí)間劣化衰減的粘結(jié)半徑、粘結(jié)面積、粘結(jié)慣性矩，A、I為巖體細(xì)觀顆粒二維粘結(jié)未衰減時(shí)的粘結(jié)面積、粘結(jié)慣性矩；

步驟204：依次計(jì)算第j個(gè)至第k個(gè)的巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的二維粘結(jié)法向彎矩增量，在二維情況下，由粘結(jié)兩端顆粒的速度、角速度和給定的循環(huán)計(jì)算時(shí)間步長增量Δt，通過公式(8)、公式(9)、公式(10)確定第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)相對(duì)轉(zhuǎn)角巖體細(xì)觀顆粒二維粘結(jié)法向增量位移以及巖體細(xì)觀顆粒二維粘結(jié)切向增量位移再結(jié)合步驟203中的公式(7)和步驟202中的公式(5)，確定第i個(gè)巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的二維粘結(jié)彎矩增量，具體見公式(11)；

其中：ff、j、k、i是自然數(shù)，且2≤j≤ff≤k，j為每次循環(huán)計(jì)算中，包含時(shí)間效應(yīng)的巖體細(xì)觀顆粒二維粘結(jié)衰減后未破裂的初始標(biāo)號(hào)值，ff為中間標(biāo)號(hào)值，k為每次循環(huán)計(jì)算中，包含時(shí)間效應(yīng)的巖體細(xì)觀顆粒二維粘結(jié)衰減后未破裂的最末標(biāo)號(hào)值，i為第一個(gè)至最后一個(gè)二維粘結(jié)顆粒標(biāo)號(hào)值，分別為第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)接觸的a端和b端的絕對(duì)運(yùn)動(dòng)速度和角速度，n_n、n_s分別為巖體細(xì)觀顆粒二維粘結(jié)接觸面的法向單位向量和切向單位向量，分別為巖體細(xì)觀顆粒二維粘結(jié)法向位移增量和切向位移增量，為巖體細(xì)觀顆粒二維粘結(jié)法向剛度，為巖體細(xì)觀顆粒二維粘結(jié)彎矩增量。

其中，巖體細(xì)觀顆粒二維粘結(jié)時(shí)效衰減劣化的初始時(shí)間步長增量Δt的確定，是采用考慮彎矩貢獻(xiàn)因子的細(xì)觀顆粒粘結(jié)時(shí)效劣化衰減的二維冪函數(shù)型模式，由每個(gè)時(shí)間步內(nèi)的細(xì)觀顆粒二維粘結(jié)首次衰減破裂所損耗的時(shí)間來確定，即通過第一個(gè)細(xì)觀顆粒粘結(jié)按冪函數(shù)型模式進(jìn)行衰減破裂所歷時(shí)的時(shí)間除以直至第一個(gè)細(xì)觀顆粒粘結(jié)破裂所需要的計(jì)算循環(huán)次數(shù)來估算初始時(shí)間步長增量Δt，見公式

其中，為第i個(gè)接觸的巖體細(xì)觀顆粒二維粘結(jié)直徑乘數(shù)，n_c為第一個(gè)巖體細(xì)觀顆粒二維粘結(jié)破裂所需的循環(huán)計(jì)算的次數(shù)，β_σ、β_τ分別為巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子、二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子，i依次為第一個(gè)至最后一個(gè)巖體細(xì)觀顆粒粘結(jié)數(shù)，∞為無窮大。

巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_σ和巖體細(xì)觀顆粒二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子β_τ的確定包括如下步驟，其中這些步驟中包含的公式下標(biāo)1代表第一個(gè)按冪函數(shù)型模式進(jìn)行時(shí)效衰減劣化的細(xì)觀顆粒二維粘結(jié)破裂標(biāo)號(hào)：

步驟211：在二維情況下，由巖體細(xì)觀顆粒粘結(jié)兩端顆粒的速度、角速度和給定的循環(huán)計(jì)算時(shí)間步長增量Δt_c，通過公式確定細(xì)觀顆粒粘結(jié)接觸的相對(duì)轉(zhuǎn)角通過公式確定細(xì)觀顆粒粘結(jié)法向增量位移通過公式確定顆粒粘結(jié)切向增量位移通過公式確定顆粒粘結(jié)接觸的彎矩增量；

步驟212：根據(jù)步驟211中的公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)法向力；根據(jù)步驟211中的公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)切向力；根據(jù)步驟211中的公式和公式通過公式確定細(xì)觀顆粒接觸的粘結(jié)切向彎矩；通過公式確定細(xì)觀顆粒接觸的粘結(jié)法向彎矩，其中，+＝為加法自反運(yùn)算符，-＝為減法自反運(yùn)算符；

步驟213：在二維情況下，通過公式確定細(xì)觀顆粒接觸的粘結(jié)正應(yīng)力，通過公式確定細(xì)觀顆粒接觸的粘結(jié)剪應(yīng)力，將這兩個(gè)公式中A、I以及用A′、I′及替換，然后將步驟203中的公式(6)和公式(7)以及步驟202中的公式(5)代入，獲得細(xì)觀顆粒粘結(jié)包含冪函數(shù)型時(shí)間效應(yīng)和考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力計(jì)算公式和包含冪函數(shù)型時(shí)間效應(yīng)的二維粘結(jié)剪應(yīng)力計(jì)算公式

步驟214：將代入公式并令β＝β_σ；將代入公式并令β＝β_σ，據(jù)此，分別得到所述巖體細(xì)觀顆粒二維粘結(jié)拉伸強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子以及巖體細(xì)觀顆粒二維粘結(jié)剪切強(qiáng)度對(duì)應(yīng)的時(shí)效劣化因子

步驟205：根據(jù)步驟203中的公式(6)和公式(7)、步驟204中的公式(8)、公式(9)和公式(11)以及步驟202中的公式(5)，依次更新計(jì)算第j個(gè)至第k個(gè)巖體細(xì)觀顆粒粘結(jié)未破裂的粘結(jié)接觸并包含時(shí)間效應(yīng)的二維粘結(jié)法向力、切向力和切向彎矩，通過公式(12)、公式(13)、公式(14)計(jì)算第i個(gè)巖體細(xì)觀顆粒二維粘結(jié)接觸的粘結(jié)法向力、切向力和切向彎矩，在二維情況下，通過公式(15)來確定巖體細(xì)觀顆粒粘結(jié)法向彎矩，

法向力：

切向力：

切向彎矩：

法向彎矩：

其中：分別為第i個(gè)巖體細(xì)觀顆粒包含時(shí)間效應(yīng)的粘結(jié)法向力、粘結(jié)切向力、包含時(shí)間效應(yīng)的粘結(jié)法向彎矩、粘結(jié)切向彎矩、粘結(jié)法向位移增量和粘結(jié)切向位移增量，為巖體細(xì)觀顆粒二維粘結(jié)切向剛度，+＝為加法自反運(yùn)算符，-＝為減法自反運(yùn)算符。

步驟206：設(shè)置彎矩貢獻(xiàn)因子考慮彎矩對(duì)巖體細(xì)觀顆粒粘結(jié)法向應(yīng)力的貢獻(xiàn)程度，根據(jù)巖體細(xì)觀顆粒粘結(jié)二維正應(yīng)力計(jì)算公式和二維粘結(jié)剪應(yīng)力計(jì)算公式同時(shí)將這兩個(gè)公式中A、I以及用A′、I′及替換，然后將步驟203中的公式(6)和公式(7)以及步驟202中的公式(5)代入，獲得包含冪函數(shù)型時(shí)間效應(yīng)和細(xì)觀顆粒粘結(jié)考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力計(jì)算公式和二維粘結(jié)剪應(yīng)力計(jì)算公式，分別見公式(16)和公式(17)，

步驟207：將步驟206中包含時(shí)間效應(yīng)的代入公式(18)，確定考慮彎矩貢獻(xiàn)因子且?guī)Ю旖刂瓜薜哪?庫倫細(xì)觀顆粒粘結(jié)時(shí)效破裂準(zhǔn)則，并且依次計(jì)算第j個(gè)至第k個(gè)的巖體細(xì)觀顆粒二維粘結(jié)應(yīng)力，用于判斷巖體細(xì)觀顆粒粘結(jié)是否破裂以及破裂模式，在該準(zhǔn)則的巖體細(xì)觀顆粒粘結(jié)應(yīng)力中包含了冪函數(shù)型時(shí)間效應(yīng)和考慮彎矩貢獻(xiàn)因子，

其中：f^s、fⁿ分別為巖體細(xì)觀顆粒二維粘結(jié)的時(shí)效剪切破裂準(zhǔn)則、時(shí)效拉伸破裂準(zhǔn)則，為第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)的二維粘結(jié)剪應(yīng)力，為第i個(gè)接觸的含冪函數(shù)型時(shí)間效應(yīng)且考慮彎矩貢獻(xiàn)因子的二維粘結(jié)正應(yīng)力，分別為巖體細(xì)觀顆粒二維粘結(jié)的拉伸強(qiáng)度、抗剪強(qiáng)度，為巖體細(xì)觀顆粒二維粘結(jié)的粘聚力，為巖體細(xì)觀顆粒二維粘結(jié)的內(nèi)摩擦角；f^s大于等于0表示巖體細(xì)觀顆粒粘結(jié)剪切破裂，小于0表示巖體細(xì)觀顆粒粘結(jié)未發(fā)生剪切破裂；fⁿ大于等于0表示巖體細(xì)觀顆粒粘結(jié)拉伸破裂，小于0表示巖體細(xì)觀顆粒粘結(jié)未發(fā)生拉伸破裂；

步驟208：當(dāng)步驟207中的公式(18)中的f^s或fⁿ大于等于0時(shí)，表明巖體細(xì)觀顆粒粘結(jié)發(fā)生破裂，此時(shí)巖體細(xì)觀顆粒的運(yùn)動(dòng)模式采用考慮阻尼效應(yīng)的二維線性接觸模型來表達(dá)；當(dāng)步驟207中的公式(18)中的f^s和fⁿ都小于0時(shí)，表明巖體細(xì)觀顆粒粘結(jié)未破裂，繼續(xù)循環(huán)步驟201至207，計(jì)算、更新、判斷巖體細(xì)觀顆粒接觸的粘結(jié)狀態(tài)，直至巖體不產(chǎn)生新的細(xì)觀顆粒粘結(jié)破裂或者巖體細(xì)觀顆粒粘結(jié)破裂加速發(fā)展而形成宏觀破壞，循環(huán)終止。

所述巖體細(xì)觀顆粒粘結(jié)發(fā)生破裂后，巖體細(xì)觀顆粒的運(yùn)動(dòng)模式采用考慮阻尼效應(yīng)的二維線性接觸模型來表達(dá)，用于描述巖體細(xì)觀顆粒粘結(jié)時(shí)效破裂后細(xì)觀顆粒的應(yīng)力、變形及運(yùn)行規(guī)律，考慮阻尼效應(yīng)的二維線性接觸模型的構(gòu)建包括如下步驟：

步驟301：通過Monte Carlo搜索算法，遍歷尋找?guī)r體細(xì)觀顆粒每個(gè)二維線性接觸端a、二維線性接觸端b(顆粒與顆粒、顆粒與墻)的中心坐標(biāo)，在二維情況下，通過公式(19)計(jì)算接觸點(diǎn)a端，接觸點(diǎn)b端的中心距離：

其中：d為二維線性接觸兩端顆粒與顆?；蛘哳w粒與墻之間的中心距離，為二維線性接觸端a的坐標(biāo)，為二維線性接觸端b的坐標(biāo)。

步驟302：二維平面狀態(tài)下巖體細(xì)觀顆粒間每個(gè)接觸點(diǎn)的單位向量通過公式(20)計(jì)算，如果是顆粒與顆粒之間的接觸，利用步驟301中得到二維線性接觸兩端的中心點(diǎn)坐標(biāo)及距離計(jì)算；如果是顆粒與墻接觸，直接利用墻體的法向量等效替換來計(jì)算；確定每個(gè)接觸端的單位向量：

其中：n_i為接觸的單位矢量，為接觸端b的方向向量，為接觸端a的方向向量，n_wall為約束墻的方向向量；

步驟303：巖體細(xì)觀顆粒粘結(jié)破裂后，每一個(gè)二維線性接觸點(diǎn)的接觸重疊量U，通過步驟301計(jì)算顆粒間距d以及二維線性接觸兩端的顆粒半徑R^a、R^b，再利用公式(21)來確定；通過設(shè)定顆粒二維線性接觸參考距離g_r，并結(jié)合公式(22)，確定顆粒二維線性接觸的距離g_s：

g_s＝|U|-g_r (22)

步驟304：確定巖體細(xì)觀顆粒接觸點(diǎn)法向、切向等效剛度，利用接觸兩端顆粒實(shí)體或是墻體的剛度k^a，k^b等效代替為接觸點(diǎn)的剛度，按公式(23)計(jì)算：

其中：K_n、K_s為巖體細(xì)觀顆粒接觸點(diǎn)等效的法向剛度和切向剛度，為顆粒與顆?；蛘哳w粒與墻的接觸a端的法向剛度和切向剛度，為顆粒與顆?；蛘哳w粒與墻的接觸b端的法向剛度和切向剛度。

步驟305：確定巖體中接觸兩端顆粒間的相對(duì)運(yùn)動(dòng)速度，利用公式(24)、公式(25)來計(jì)算。其中e_pqz為Ricci指標(biāo)置換符號(hào)，按照公式(26)來計(jì)算：

其中：V_p與V_q等價(jià)，V_p與V_q為巖體中接觸兩端顆粒間的相對(duì)運(yùn)動(dòng)速度，p、q為指標(biāo)等價(jià)符號(hào)，p＝1，q＝1表示顆粒與顆粒接觸，p＝2，q＝2時(shí)表示顆粒與墻接觸，為顆粒與顆?；蚴穷w粒與墻的接觸b端單元的速度，為顆粒與顆?；蚴穷w粒與墻的接觸a端單元的速度，為顆粒與顆?；蚴穷w粒與墻的接觸a端單元的角速度，為顆粒與顆?；蚴穷w粒與墻的接觸b端單元的角速度，為顆粒與顆?；蚴穷w粒與墻的接觸a端的位移，為顆粒與顆粒或是顆粒與墻的接觸b端的位移，為位移指標(biāo)變換的中間過渡符號(hào)，表示指標(biāo)符號(hào)為p時(shí)顆粒-顆?；蛘哳w粒-墻的接觸a端的速度，表示指標(biāo)符號(hào)為q時(shí)顆粒-顆?；蛘哳w粒-墻的接觸a端的速度，表示指標(biāo)符號(hào)為p時(shí)顆粒-顆?；蛘哳w粒-墻的接觸b端的速度，表示指標(biāo)符號(hào)為q時(shí)顆粒-顆?；蛘哳w粒-墻的接觸b端的速度(只有a端和b端兩個(gè)接觸端)。

步驟306：對(duì)于巖體細(xì)觀顆粒線性接觸模型的初始時(shí)間步長增量Δt的取值，通過公式(29)估計(jì)最小的時(shí)間步長Δt，確保所構(gòu)建模型的計(jì)算時(shí)間步長小于該值，即可保證系統(tǒng)積分計(jì)算趨于穩(wěn)定；通過公式(30)、公式(31)、公式(32)確定每個(gè)線性接觸的總位移增量、法向位移增量和切向位移增量：

R＝min(R_a，R_b) (27)

ΔU_p1＝V_p1Δt (30)

其中：R為巖體細(xì)觀顆粒的等效半徑，m為巖體細(xì)觀顆粒質(zhì)量，J1為巖體細(xì)觀顆粒的轉(zhuǎn)動(dòng)慣量；k^平為巖體細(xì)觀顆粒系統(tǒng)平移剛度，k^轉(zhuǎn)為巖體細(xì)觀顆粒系統(tǒng)旋轉(zhuǎn)剛度；ΔU_p1為巖體細(xì)觀顆粒二維線性接觸的總位移增量，Δδ_n、物理意義相同，均表示巖體細(xì)觀顆粒二維線性接觸的法向位移增量，Δδ_s、物理意義相同，均表示巖體細(xì)觀顆粒二維線性接觸的切向位移增量，V_p1與V_q1為巖體細(xì)觀顆粒接觸兩端的相對(duì)運(yùn)動(dòng)速度，n為單位法向量，p1，q1為張量指標(biāo)變換符號(hào)。

步驟307：由公式(22)判定巖體細(xì)觀顆粒表面接觸允許存在的最大距離，計(jì)算法向和切向位移更新因子，另外，巖體細(xì)觀顆粒二維線性接觸法向位移增量的更新是采用前一步的法向位移增量與更新因子α的乘積獲得，巖體細(xì)觀顆粒二維線性接觸切向位移增量的更新是采用前一步的切向位移增量與更新因子α的乘積獲得。

其中：(g_s)₀為模型計(jì)算初始時(shí)刻的表面接觸距離，g_s為巖體細(xì)觀顆粒接觸的距離，α為位移更新因子。

步驟308：巖體細(xì)觀顆粒二維線性接觸的法向線性力采取相對(duì)矢量累加(M_l＝1)和絕對(duì)矢量累加(M_l＝0)模式，通過公式(33)、(34)計(jì)算獲得；巖體細(xì)觀顆粒二維線性接觸的切向線性力采用巖體細(xì)觀顆粒接觸滑動(dòng)來表示，通過公式(35)計(jì)算獲得：

其中：分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向線性力、切向線性力，k_n、k_s分別為巖體細(xì)觀顆粒間受力變形的二維線性接觸法向、切向線性剛度，Δδ_n、Δδ_s分別為巖體細(xì)觀顆粒二維線性接觸的法向位移增量、切向位移增量，分別為巖體細(xì)觀顆粒二維線性接觸的初始法向力增量值、切向力增量值，為巖體細(xì)觀顆粒未滑動(dòng)時(shí)的靜摩擦力，為巖體細(xì)觀顆?；瑒?dòng)摩擦力，其值可通過摩擦系數(shù)μ與乘積得到。

步驟309：巖體細(xì)觀顆粒線性接觸的法向阻尼采用全法向模式M_d＝{0，2}和無拉伸模式M_d＝{1，3}兩種，通過公式(36)計(jì)算，其中m_c為等效顆粒質(zhì)量，按公式(37)計(jì)算，巖體細(xì)觀顆粒線性接觸的切向阻尼采用全剪切模式M_d＝{0，1}和滑-剪模式M_d＝{2，3}，按照公式(38)來計(jì)算，

其中：分別為巖體細(xì)觀顆粒線性接觸的法向線性阻尼力、切向線性阻尼力，β_n為巖體細(xì)觀顆粒線性接觸的法向阻尼系數(shù)，β_s為巖體細(xì)觀顆粒線性接觸的切向阻尼系數(shù)，k_n為巖體細(xì)觀顆粒線性接觸的法向線性剛度，k_s為巖體細(xì)觀顆粒線性接觸的切向線性剛度，分別為巖體細(xì)觀顆粒線性接觸的法向速率和切向速率，F(xiàn)^*為巖體細(xì)觀顆粒線性接觸的全法向阻尼力，表達(dá)式為m_c為巖體細(xì)觀等效顆粒質(zhì)量，m⁽¹⁾為巖體細(xì)觀顆粒接觸端1的細(xì)觀顆粒質(zhì)量，m⁽²⁾為巖體細(xì)觀顆粒接觸端2的細(xì)觀顆粒質(zhì)量，F(xiàn)^d為巖體細(xì)觀顆粒線性接觸的總阻尼力。

下面以深部巖體為實(shí)例，結(jié)合附圖詳述本發(fā)明模型的數(shù)值實(shí)現(xiàn)的具體過程，請(qǐng)參閱實(shí)例附圖說明中的圖12至圖13以及模型附圖說明中的圖1至圖10，來理解本發(fā)明模型的數(shù)值實(shí)現(xiàn)步驟及效果：

步驟1：采用C++編程語言，并結(jié)合fish語言，根據(jù)本發(fā)明的模型結(jié)構(gòu)構(gòu)建流程圖(圖11)，在數(shù)值平臺(tái)上實(shí)現(xiàn)了巖體二維細(xì)觀時(shí)效破裂冪函數(shù)型模型。

步驟2：初步確定巖體時(shí)效破裂模型的細(xì)觀參數(shù)

粒徑比Rratio、線性接觸法向剛度kn(圖6)、線性接觸切向剛度ks(圖6)、顆粒密度ba_rho、顆粒接觸模量b_Ec、粘結(jié)法向剛度pb_kn(圖6)、粘結(jié)切向剛度pb_ks(圖6)、粘結(jié)模型pb_Ec、顆粒的摩擦系數(shù)ba_fric、粘結(jié)拉伸強(qiáng)度的平均值pb_sn_mean、粘結(jié)拉伸強(qiáng)度的標(biāo)準(zhǔn)差pb_sn_sdev、粘聚力平均值pb_coh_mean、粘聚力標(biāo)準(zhǔn)差pb_coh_sdev、粘結(jié)半徑乘數(shù)gamma(圖7)、粘結(jié)彎矩奉獻(xiàn)因子beta、法向阻尼系數(shù)Apfan(圖6)、切向阻尼系數(shù)Apfas(圖6)以及內(nèi)摩擦角pb_phi(圖8)等19個(gè)參數(shù)，參數(shù)具體值見表一。

步驟3：生成巖體模型

根據(jù)高斯分布或weibull分布確定模型的粘結(jié)拉伸強(qiáng)度和粘聚力分布，通過均勻隨機(jī)函數(shù)分布法確定顆粒的粒徑分布；通過各向同性應(yīng)力調(diào)整法及自適應(yīng)動(dòng)態(tài)膨脹法，調(diào)整顆粒的位置，減少顆粒重疊量；通過懸浮顆粒刪除法，刪除孤立顆粒，提高模型樣本的整體性，減少缺陷模型的生成。最后賦予模型材料粘結(jié)強(qiáng)度參數(shù)，生成具有真實(shí)巖體結(jié)構(gòu)的巖體模型。巖體模型的直徑為50mm、高度為100mm(圖12)。

步驟4：精確標(biāo)定本發(fā)明中模型的細(xì)觀物理力學(xué)參數(shù)

通過室內(nèi)單軸和三軸壓縮試驗(yàn)得到的應(yīng)力-應(yīng)變曲線，確定巖體的宏觀彈性模量峰值強(qiáng)度σ_p，以及泊松比通過優(yōu)化方法，使巖體單、三軸壓縮模型的應(yīng)力-應(yīng)變曲線與室內(nèi)試驗(yàn)的應(yīng)力-應(yīng)變以及宏觀變形參數(shù)和強(qiáng)度參數(shù)吻合，獲得本發(fā)明所構(gòu)建模型的細(xì)觀物理力學(xué)參數(shù)。

步驟5：巖體時(shí)效力學(xué)參數(shù)標(biāo)定

對(duì)巖體進(jìn)行一系列不同應(yīng)力強(qiáng)度比條件下的時(shí)效力學(xué)試驗(yàn)，得到不同應(yīng)力強(qiáng)度比條件下巖體變形隨時(shí)間演化的曲線。通過參數(shù)擬合法，匹配實(shí)際巖體的時(shí)效變形過程，確定巖體細(xì)觀顆粒粘結(jié)時(shí)效劣化的第一控制參數(shù)β₁、第二控制參數(shù)β₂。

步驟6：巖體時(shí)效力學(xué)數(shù)值試驗(yàn)

在荷載一定的條件下，分別設(shè)定不同的彎矩貢獻(xiàn)因子，獲得了巖體時(shí)效變形破壞的演化規(guī)律(圖13)。

本發(fā)明模型的參數(shù)名稱及值如表一所示。

表一：本發(fā)明模型的參數(shù)名稱及值

上述實(shí)施例中，公式的符號(hào)與圖1～圖10以及附圖說明中的符號(hào)相互對(duì)應(yīng)。

其它未詳細(xì)說明的部分均為現(xiàn)有技術(shù)，以上所有參數(shù)均可通過查閱手冊(cè)或計(jì)算得到。本發(fā)明并不嚴(yán)格地局限于上述實(shí)施例。以上所述僅為本發(fā)明的特定實(shí)施例，并不用于限制本發(fā)明。凡在本發(fā)明的精神和原則之內(nèi)所做的任何修改、等同替換及改進(jìn)等，都在本發(fā)明的保護(hù)范圍之內(nèi)。