趙文娟

摘要：巖石力學(xué)試驗(yàn)是研究巖石力學(xué)特性，獲取力學(xué)參數(shù)的重要途徑。為研究不同荷載作用下不同尺寸試塊的破壞過(guò)程，利用Python語(yǔ)言對(duì)ABAQUS軟件進(jìn)行二次開(kāi)發(fā)，建立巖石力學(xué)試驗(yàn)過(guò)程的參數(shù)化數(shù)值模擬仿真系統(tǒng)。將考慮剪切效應(yīng)的Drucker-Prager破壞準(zhǔn)則編入vusdfld子程序中，引入失效單元?jiǎng)h除算法研究荷載作用下巖石試塊變形破壞過(guò)程，實(shí)現(xiàn)了巖石力學(xué)試驗(yàn)過(guò)程的數(shù)值仿真。以花崗巖試塊三軸壓縮試驗(yàn)為例，驗(yàn)證了仿真系統(tǒng)的有效性。該系統(tǒng)可根據(jù)用戶需要與實(shí)驗(yàn)室加載平臺(tái)進(jìn)行數(shù)據(jù)匹配，極大提高了建模和分析效率，同時(shí)豐富了力學(xué)試驗(yàn)內(nèi)容。

關(guān)鍵詞：巖石試樣;二次開(kāi)發(fā);數(shù)值仿真;力學(xué)試驗(yàn)

中圖分類號(hào)：G642423;TU455 ? 文獻(xiàn)標(biāo)志碼：A ? 文章編號(hào)：1005-2909（2019）05-0160-06

巖石力學(xué)試驗(yàn)是研究巖石力學(xué)特性，獲取力學(xué)參數(shù)的重要途徑，其應(yīng)用涉及土木建筑、采礦、水利水電、地下工程等眾多工程領(lǐng)域。同時(shí)，隨著現(xiàn)代力學(xué)、數(shù)學(xué)、計(jì)算機(jī)技術(shù)的不斷發(fā)展，數(shù)值仿真以其模擬分析的有效性、計(jì)算結(jié)果的三維可視化等優(yōu)點(diǎn)，為工程問(wèn)題研究提供了有效手段[1-2]。

近年來(lái)，關(guān)于巖石力學(xué)試驗(yàn)的數(shù)值模擬已成為巖石力學(xué)特性分析領(lǐng)域的研究熱點(diǎn)。賈善坡等[3]利用數(shù)值模擬方法建立泥巖試樣二維簡(jiǎn)化模型，并根據(jù)非排水條件下泥巖室內(nèi)三軸試驗(yàn)結(jié)果，反演研究巖石本構(gòu)模型參數(shù)。田洪銘等[4]結(jié)合泥質(zhì)紅砂巖的三軸蠕變?cè)囼?yàn)結(jié)果，建立非線性蠕變損傷模型，并利用有限元軟件建立與室內(nèi)試驗(yàn)標(biāo)準(zhǔn)巖樣尺寸一致的模型，對(duì)試驗(yàn)過(guò)程進(jìn)行模擬。王杰等[5]通過(guò)單軸壓縮試驗(yàn)， 對(duì)巖石的破壞過(guò)程進(jìn)行分析， 并通過(guò)數(shù)值模擬獲得了裂紋萌生、擴(kuò)展、貫通直至形成宏觀裂縫的過(guò)程。王青元等[6]采用所提出的模型對(duì)7個(gè)不同尺寸的綠砂巖巖樣進(jìn)行了單軸壓縮蠕變數(shù)值模擬，并對(duì)巖石長(zhǎng)期強(qiáng)度尺寸效應(yīng)進(jìn)行了分析。趙志濤等[7]建立了砂卵石-砂土復(fù)合試樣有限元數(shù)值模型，通過(guò)不同圍壓、不同砂層厚度的復(fù)合試樣數(shù)值試驗(yàn)，探究砂層對(duì)復(fù)合試樣強(qiáng)度和破壞規(guī)律的影響。劉洪磊等[8]模擬了花崗巖單軸壓縮試驗(yàn)過(guò)程，對(duì)其破裂機(jī)理進(jìn)行研究。

通過(guò)以上文獻(xiàn)可知，目前針對(duì)巖石力學(xué)試驗(yàn)的模擬主要側(cè)重于巖石材料本構(gòu)模型的二次開(kāi)發(fā)或針對(duì)特定尺寸、特定加載方式模型的分析，而針對(duì)試驗(yàn)過(guò)程尤其是模擬分析的前處理過(guò)程（模型尺寸、加載參數(shù)的確定等）和后處理過(guò)程（模擬結(jié)果的輸出、數(shù)據(jù)的采集）等方面的研究成果較少。文章利用Python語(yǔ)言對(duì)ABAQUS軟件進(jìn)行二次開(kāi)發(fā)，建立巖石力學(xué)試驗(yàn)過(guò)程的參數(shù)化數(shù)值模擬仿真系統(tǒng)，將考慮剪切效應(yīng)的Drucker-Prager破壞準(zhǔn)則編入vusdfld子程序中，引入單元生死算法研究荷載作用下巖石試塊變形破壞過(guò)程，實(shí)現(xiàn)了巖石力學(xué)試驗(yàn)過(guò)程的數(shù)值仿真，最后，以花崗巖試塊三軸壓縮試驗(yàn)為例，驗(yàn)證了仿真系統(tǒng)的有效性。

一、巖石力學(xué)試驗(yàn)數(shù)值仿真系統(tǒng)的開(kāi)發(fā)

（一） 基于Python語(yǔ)言的ABAQUS二次開(kāi)發(fā)流程

ABAQUS大型有限元軟件具有強(qiáng)大的工程問(wèn)題分析能力，擁有豐富的模型材料庫(kù)和單元庫(kù)，能模擬復(fù)雜的非線性問(wèn)題。同時(shí)，ABAQUS軟件具備很強(qiáng)的通用性和擴(kuò)展性，允許用戶通過(guò)自行編程實(shí)現(xiàn)對(duì)軟件前、后處理的二次開(kāi)發(fā)。Python語(yǔ)言作為ABAQUS軟件的標(biāo)準(zhǔn)設(shè)計(jì)語(yǔ)言，為用戶自行編程，實(shí)現(xiàn)對(duì)ABAQUS的操控。

通過(guò)基于Python語(yǔ)言的ABAQUS腳本接口，可以直接調(diào)用軟件豐富的庫(kù)函數(shù)，并繞過(guò)用戶界面直接操縱內(nèi)核，實(shí)現(xiàn)工程建模、材料屬性定義、模型網(wǎng)格劃分、邊界條件設(shè)置、結(jié)果數(shù)據(jù)查看和分析等功能，還可對(duì)用戶界面進(jìn)行修改和二次開(kāi)發(fā)[9]。

利用Python語(yǔ)言編程對(duì)ABAQUS進(jìn)行前后處理的二次開(kāi)發(fā)，實(shí)現(xiàn)如下功能：

（1）前處理。通過(guò)編程實(shí)現(xiàn)用戶自定義交互式模型尺寸、材料屬性、邊界條件等參數(shù)的輸入，并通過(guò)有限元計(jì)算，輸出不同參數(shù)下的計(jì)算結(jié)果文件。

（2）后處理。通過(guò)編程實(shí)現(xiàn)對(duì)計(jì)算結(jié)果文件的讀取、數(shù)據(jù)分析和用戶指定結(jié)果的顯示。

具體流程見(jiàn)圖1。

（二） 系統(tǒng)開(kāi)發(fā)

適用于模擬不同模型尺寸、材料屬性、邊界條件的巖石試件試驗(yàn)破壞及過(guò)程分析數(shù)值仿真系統(tǒng)，實(shí)現(xiàn)了不重復(fù)有限元分析時(shí)繁瑣的前后處理操作，僅通過(guò)程序調(diào)用，在交互窗口中輸入不同計(jì)算參數(shù)就可以自動(dòng)進(jìn)行建模、計(jì)算，并形成相對(duì)應(yīng)的計(jì)算結(jié)果文件。巖石試件有限元模型如圖2所示。

巖石力學(xué)試驗(yàn)過(guò)程中，通常將試塊設(shè)計(jì)為圓柱形，然后將其放置在加載系統(tǒng)上，通過(guò)電腦操作對(duì)試塊進(jìn)行加載，研究其力學(xué)特性。通常，試驗(yàn)參數(shù)主要包括試樣幾何參數(shù)（試樣的尺寸）和加載參數(shù)（圍壓、軸壓的量值）。為與實(shí)驗(yàn)室試驗(yàn)過(guò)程一致，巖石力學(xué)試驗(yàn)數(shù)值仿真系統(tǒng)開(kāi)發(fā)的具體實(shí)施步驟包括：

（1）將圖2中加載參數(shù)（軸壓P1、圍壓P2）、試樣幾何參數(shù)（高H、直徑D）以及試塊力學(xué)參數(shù)（密度、彈性模量、泊松比、粘聚力、摩擦角）設(shè)定為可變化的參數(shù)，利用Python語(yǔ)言編程進(jìn)行參數(shù)化建模[10]。

（2）在參數(shù)化建模程序中加上如圖3所示的交互式模塊程序代碼，使程序代碼中的參數(shù)與試驗(yàn)參數(shù)相匹配，當(dāng)程序調(diào)用后首先進(jìn)入用戶自定義交互式窗口，即可對(duì)數(shù)值模擬試驗(yàn)參數(shù)進(jìn)行修改，如圖4所示。

（3）對(duì)不同試驗(yàn)參數(shù)下的模型進(jìn)行有限元計(jì)算。

####Enter the parameter of pressure#####

fields = （（'Axial pressure P1：′，′5000′）， （′Confining pressure P2：′，′2000′））

iP1 = None

while iP1 == None：

[6]王青元，朱萬(wàn)成，劉洪磊，等. 單軸壓縮下綠砂巖長(zhǎng)期強(qiáng)度的尺寸效應(yīng)研究[J]. 巖土力學(xué)，2016，37（4）：981-990.

[7]趙志濤，張子新，朱葉艇. 砂卵石-砂土復(fù)合試樣大三軸數(shù)值試驗(yàn)及驗(yàn)證[J]. 地下空間與工程學(xué)報(bào)，2014，10（6）：1300-1306.

[8]劉洪磊，王培濤，楊天鴻，等. 基于離散元方法的花崗巖單軸壓縮破裂過(guò)程的聲發(fā)射特性[J]. 煤炭學(xué)報(bào)，2015，40（8）：1790-1795.

[9]岳永保，楊合，詹梅，等. 基于Python的ABAQUS數(shù)控彎管數(shù)值模擬后處理[J]. 塑性工程學(xué)報(bào)，2009，16（3）：48-54.

[10]王輝，陳衛(wèi)忠，李廷春，等. 地下工程智能連續(xù)優(yōu)化方法的實(shí)現(xiàn)Ⅰ：參數(shù)化設(shè)計(jì)及優(yōu)化方法[J]. 水利學(xué)報(bào)，2014，45（1）：42-49.

[11]杜明瑞， 靖洪文， 蘇海健，等. 孔洞形狀對(duì)砂巖強(qiáng)度及破壞特征的影響[J]. 工程力學(xué)， 2016，33（7）：190-196.

[12]孫華飛，楊永明，鞠楊，等. 開(kāi)挖卸荷條件下煤巖變形破壞與能量釋放的數(shù)值分析[J]. 煤炭學(xué)報(bào)，2014，39（2）：258-272.

Abstract： Rock mechanics test is an important way to study rock mechanics properties and obtain mechanical parameters. In order to research the damage process of rock specimens with various sizes under different load levels， the Python language is used to make secondary development for ABAQUS and the parameterized numerical simulation system for rock mechanics test is established. The Drucker-Prager failure criterion considering shear effect is programmed into VUSDFLD subroutine and the failure element deletion algorithm is introduced to study the deformation and damage process of rock specimens. Then the triaxial compression test of granite block is taken as an example to verify the validity of the simulation system. The system can match with the laboratory loading platform according to the users need， improve the efficiency of modeling and analysis， and enrich the content of mechanics test.

Key words： rock specimens; secondary development; numerical simulation; mechanics experiment

（責(zé)任編輯 周 沫）