基于ABAQUS二次开发的岩石力学试验数值仿真系统

赵文娟

摘要：岩石力学试验是研究岩石力学特性，获取力学参数的重要途径。为研究不同荷载作用下不同尺寸试块的破坏过程，利用Python语言对ABAQUS软件进行二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统。将考虑剪切效应的Drucker-Prager破坏准则编入vusdfld子程序中，引入失效单元删除算法研究荷载作用下岩石试块变形破坏过程，实现了岩石力学试验过程的数值仿真。以花岗岩试块三轴压缩试验为例，验证了仿真系统的有效性。该系统可根据用户需要与实验室加载平台进行数据匹配，极大提高了建模和分析效率，同时丰富了力学试验内容。

关键词：岩石试样;二次开发;数值仿真;力学试验

中图分类号：G642423;TU455   文献标志码：A   文章编号：1005-2909（2019）05-0160-06

岩石力学试验是研究岩石力学特性，获取力学参数的重要途径，其应用涉及土木建筑、采矿、水利水电、地下工程等众多工程领域。同时，随着现代力学、数学、计算机技术的不断发展，数值仿真以其模拟分析的有效性、计算结果的三维可视化等优点，为工程问题研究提供了有效手段[1-2]。

近年来，关于岩石力学试验的数值模拟已成为岩石力学特性分析领域的研究热点。贾善坡等[3]利用数值模拟方法建立泥岩试样二维简化模型，并根据非排水条件下泥岩室内三轴试验结果，反演研究岩石本构模型参数。田洪铭等[4]结合泥质红砂岩的三轴蠕变试验结果，建立非线性蠕变损伤模型，并利用有限元软件建立与室内试验标准岩样尺寸一致的模型，对试验过程进行模拟。王杰等[5]通过单轴压缩试验， 对岩石的破坏过程进行分析， 并通过数值模拟获得了裂纹萌生、扩展、贯通直至形成宏观裂缝的过程。王青元等[6]采用所提出的模型对7个不同尺寸的绿砂岩岩样进行了单轴压缩蠕变数值模拟，并对岩石长期强度尺寸效应进行了分析。赵志涛等[7]建立了砂卵石-砂土复合试样有限元数值模型，通过不同围压、不同砂层厚度的复合试样数值试验，探究砂层对复合试样强度和破坏规律的影响。刘洪磊等[8]模拟了花岗岩单轴压缩试验过程，对其破裂机理进行研究。

通过以上文献可知，目前针对岩石力学试验的模拟主要侧重于岩石材料本构模型的二次开发或针对特定尺寸、特定加载方式模型的分析，而针对试验过程尤其是模拟分析的前处理过程（模型尺寸、加载参数的确定等）和后处理过程（模拟结果的输出、数据的采集）等方面的研究成果较少。文章利用Python语言对ABAQUS软件进行二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统，将考虑剪切效应的Drucker-Prager破坏准则编入vusdfld子程序中，引入单元生死算法研究荷载作用下岩石试块变形破坏过程，实现了岩石力学试验过程的数值仿真，最后，以花岗岩试块三轴压缩试验为例，验证了仿真系统的有效性。

一、岩石力学试验数值仿真系统的开发

（一） 基于Python语言的ABAQUS二次开发流程

ABAQUS大型有限元软件具有强大的工程问题分析能力，拥有丰富的模型材料库和单元库，能模拟复杂的非线性问题。同时，ABAQUS软件具备很强的通用性和扩展性，允许用户通过自行编程实现对软件前、后處理的二次开发。Python语言作为ABAQUS软件的标准设计语言，为用户自行编程，实现对ABAQUS的操控。

通过基于Python语言的ABAQUS脚本接口，可以直接调用软件丰富的库函数，并绕过用户界面直接操纵内核，实现工程建模、材料属性定义、模型网格划分、边界条件设置、结果数据查看和分析等功能，还可对用户界面进行修改和二次开发[9]。

利用Python语言编程对ABAQUS进行前后处理的二次开发，实现如下功能：

（1）前处理。通过编程实现用户自定义交互式模型尺寸、材料属性、边界条件等参数的输入，并通过有限元计算，输出不同参数下的计算结果文件。

（2）后处理。通过编程实现对计算结果文件的读取、数据分析和用户指定结果的显示。

具体流程见图1。

（二） 系统开发

适用于模拟不同模型尺寸、材料属性、边界条件的岩石试件试验破坏及过程分析数值仿真系统，实现了不重复有限元分析时繁琐的前后处理操作，仅通过程序调用，在交互窗口中输入不同计算参数就可以自动进行建模、计算，并形成相对应的计算结果文件。岩石试件有限元模型如图2所示。

岩石力学试验过程中，通常将试块设计为圆柱形，然后将其放置在加载系统上，通过电脑操作对试块进行加载，研究其力学特性。通常，试验参数主要包括试样几何参数（试样的尺寸）和加载参数（围压、轴压的量值）。为与实验室试验过程一致，岩石力学试验数值仿真系统开发的具体实施步骤包括：

（1）将图2中加载参数（轴压P1、围压P2）、试样几何参数（高H、直径D）以及试块力学参数（密度、弹性模量、泊松比、粘聚力、摩擦角）设定为可变化的参数，利用Python语言编程进行参数化建模[10]。

（2）在参数化建模程序中加上如图3所示的交互式模块程序代码，使程序代码中的参数与试验参数相匹配，当程序调用后首先进入用户自定义交互式窗口，即可对数值模拟试验参数进行修改，如图4所示。

（3）对不同试验参数下的模型进行有限元计算。

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

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

iP1 = None

while iP1 == None：

[6]王青元，朱万成，刘洪磊，等. 单轴压缩下绿砂岩长期强度的尺寸效应研究[J]. 岩土力学，2016，37（4）：981-990.

[7]赵志涛，张子新，朱叶艇. 砂卵石-砂土复合试样大三轴数值试验及验证[J]. 地下空间与工程学报，2014，10（6）：1300-1306.

[8]刘洪磊，王培涛，杨天鸿，等. 基于离散元方法的花岗岩单轴压缩破裂过程的声发射特性[J]. 煤炭学报，2015，40（8）：1790-1795.

[9]岳永保，杨合，詹梅，等. 基于Python的ABAQUS数控弯管数值模拟后处理[J]. 塑性工程学报，2009，16（3）：48-54.

[10]王輝，陈卫忠，李廷春，等. 地下工程智能连续优化方法的实现Ⅰ：参数化设计及优化方法[J]. 水利学报，2014，45（1）：42-49.

[11]杜明瑞， 靖洪文， 苏海健，等. 孔洞形状对砂岩强度及破坏特征的影响[J]. 工程力学， 2016，33（7）：190-196.

[12]孙华飞，杨永明，鞠杨，等. 开挖卸荷条件下煤岩变形破坏与能量释放的数值分析[J]. 煤炭学报，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

（责任编辑 周 沫）