CN108182335A - 一种基于abaqus的岩石力学试验数值仿真方法 - Google Patents

Abstract

本发明公开了一种基于ABAQUS的岩石力学试验数值仿真方法，涉及岩石力学试验技术领域。本发明利用Python语言对ABAQUS软件进行二次开发，对前处理过程(包括模型尺寸、加载参数的确定等)以及后处理过程(模拟结果的输出、数据的采集)进行编程，建立岩石力学试验过程的参数化数值模拟仿真系统；将考虑剪切效应的Drucker‑Prager破坏准则编入ABAQUS软件的vusdfld子程序中，引入失效单元删除算法研究荷载作用下岩石试块变形破坏过程，实现了岩石力学试验过程的数值仿真，并以花岗岩试块三轴压缩试验为例，验证了仿真系统的合理性和有效性；以及本发明的岩石力学试验数值仿真系统可根据用户需要与实验室加载平台进行数据匹配，极大提高了建模和分析效率，同时丰富了力学试验内容。

Description

一种基于ABAQUS的岩石力学试验数值仿真方法

技术领域

本发明涉及岩石力学试验技术领域，更具体的涉及一种基于ABAQUS的岩石力学试验数值仿真方法。

背景技术

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

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

分析以上文献可知，目前针对岩石力学试验的模拟主要侧重于岩石材料本构模型的二次开发或针对特定尺寸、特定加载方式的模型开展分析，而目前针对试验过程尤其是模拟分析的前处理过程(包括模型尺寸、加载参数的确定等)以及后处理过程(模拟结果的输出、数据的采集)等方面的有效研究成果尚未出现。

发明内容

本发明实施例提供一种基于ABAQUS的岩石力学试验数值仿真方法，用以解决现有技术中存的问题。

本发明实施例提供一种基于ABAQUS的岩石力学试验数值仿真方法，包括：

通过Python语言编程对ABAQUS软件进行前处理二次开发和后处理二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统；

将考虑剪切效应的Drucker-Prager破坏准则编入ABAQUS软件中的vusdfld子程序中，并引入失效单元删除算法模拟荷载作用下岩石试块变形破坏过程。

进一步地，所述前处理二次开发，包括：通过编程实现用户自定义交互式模型尺寸、材料属性、边界条件参数的输入，并通过有限元计算，输出不同参数下的计算结果文件。

进一步地，所述后处理二次开发，包括：通过编程实现对计算结果文件的读取、数据分析和用户指定结果的显示。

进一步地，所述通过Python语言编程对ABAQUS软件进行前处理二次开发和后处理二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统，具体包括：

将加载参数轴压P1和围压P2、试样几何参数高H和直径D、以及试块力学参数密度、弹性模量、泊松比、粘聚力和摩擦角设定为变化参数，利用Python语言编程进行参数化建模；

在参数化建模程序开头添加交互式模块程序代码，使交互式模块程序代码中的参数与试验参数相匹配，当程序调用后首先进入用户自定义交互式窗口，即对数值模拟试验参数进行修改；

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

进一步地，所述vusdfld子程序，用于读取有限元模型迭代计算过程的单元状态；其中，单元状态包括：应力、应变、等效塑性应变。

进一步地，所述将考虑剪切效应的Drucker-Prager破坏准则编入ABAQUS软件中的vusdfld子程序中，并引入失效单元删除算法模拟荷载作用下岩石试块变形破坏过程，具体包括：

当岩石试块有限元模型的某单元应力状态满足Drucker-Prager破坏准则进入塑性状态时，认为该单元失去了承载能力发生破坏，将该单元设定为失效单元，在下一加载步迭代计算中将其删除，反复迭代计算，直至所有加载步计算结束。

本发明实施例中，提供一种基于ABAQUS的岩石力学试验数值仿真方法，与现有技术相比，其有益效果如下：

本发明利用Python语言对ABAQUS软件进行二次开发，对前处理过程(包括模型尺寸、加载参数的确定等)以及后处理过程(模拟结果的输出、数据的采集)进行编程，建立岩石力学试验过程的参数化数值模拟仿真系统；将考虑剪切效应的Drucker-Prager破坏准则编入vusdfld子程序中，引入失效单元删除算法研究荷载作用下岩石试块变形破坏过程，实现了岩石力学试验过程的数值仿真，并以花岗岩试块三轴压缩试验为例，验证了仿真系统的合理性和有效性；以及本发明的岩石力学试验数值仿真系统可根据用户需要与实验室加载平台进行数据匹配，极大提高了建模和分析效率，同时丰富了力学试验内容。

附图说明

图1为本发明实施例提供的一种基于ABAQUS的岩石力学试验数值仿真方法流程图；

图2为本发明实施例提供的ABAQUS软件的二次开发流程图；

图3为本发明实施例提供的岩石试件有限元模型；

图4a为本发明实施例提供的用户自定义交互式窗口的加载参数窗口；

图4b为本发明实施例提供的用户自定义交互式窗口的试样几何参数窗口；

图4c为本发明实施例提供的用户自定义交互式窗口的试块力学参数窗口；

图5为本发明实施例提供的三轴压缩试验条件下试样应力-应变曲线图；

图6为本发明实施例提供的花岗岩试样破坏图；

图7为本发明实施例提供的加载后试块位移矢量图；

图8为本发明实施例提供的试样模拟过程应力-应变曲线图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的一种基于ABAQUS的岩石力学试验数值仿真方法流程图。如图1所示，该方法包括：

步骤S 1，通过Python语言编程对ABAQUS软件进行前处理二次开发和后处理二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统。

需要说明的是，ABAQUS大型有限元软件具有强大的工程问题分析能力，拥有丰富的模型材料库和单元库，能够模拟复杂的非线性问题。同时，ABAQUS软件具备很强的通用性和扩展性，允许用户通过自行编程实现对软件前、后处理的二次开发。Python语言作为ABAQUS软件的标准设计语言，为用户自行编程实现对ABAQUS的操控提供了手段。通过基于Python语言的ABAQUS脚本接口，可以直接调用软件丰富的库函数，并绕过用户界面直接操纵内核，实现工程建模、材料属性定义、模型网格划分、边界条件设置、结果数据查看和分析等功能，还可以对用户界面进行修改和二次开发。

在步骤S1中，前处理二次开发，包括：通过编程实现用户自定义交互式模型尺寸、材料属性、边界条件参数的输入，并通过有限元计算，输出不同参数下的计算结果文件。后处理二次开发，包括：通过编程实现对计算结果文件的读取、数据分析和用户指定结果的显示。具体流程见图2。

需要说明的是，本发明适用于模拟不同模型尺寸、材料属性、边界条件等参数下岩石试件试验破坏及过程分析的数值仿真系统，实现了不用重复进行有限元分析问题时繁琐的前后处理操作，仅通过程序调用后在交互窗口中输入不同计算参数就可以自动进行建模、计算，并形成相对应的计算结果文件。岩石试件有限元模型如图3所示。岩石力学试验过程中，通常将试块设计为圆柱形，然后将其放置在加载系统上，通过电脑人工操作对试块进行加载，研究其力学特性。通常，试验参数主要包括试样几何参数(试样的尺寸)和加载参数(围压、轴压的量值)。为与实验室试验过程一致，本发明岩石力学试验数值仿真系统开发的具体实施步骤包括：

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

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

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

表1交互式模块程序代码

需要说明的是，表1中第一部分程序段为加载参数模块代码；表1中第二部分程序段为试样几何参数模块代码；表1中第三部分程序段为试块力学参数模块代码。

步骤S 2，将考虑剪切效应的Drucker-Prager破坏准则编入ABAQUS软件中的vusdfld子程序中，并引入失效单元删除算法模拟荷载作用下岩石试块变形破坏过程。

需要说明的是，岩石试样在压缩作用下，破坏模式主要包括剪切破坏、张拉破坏和拉剪混合破坏等。基于此，设置考虑剪切效应的Drucker-Prager准则(简称D-P准则)作为模型单元破坏的判断准则，以模型加载计算过程单元体是否达到以上状态作为判断其失效的依据。

D-P准则可表示为：

其中，I₁＝σ₁+σ₂+σ₃表示应力张量的第一不变量；

表示偏应力张量的第二不变量；

α和K表示岩石的材料参数，可由以下公式计算：

式中，C和分别表示岩石的粘聚力和内摩擦角。

进一步地，ABAQUS中vusdfld子程序可以读取模型迭代计算过程的单元状态(如应力、应变、等效塑性应变等)，将考虑剪切效应的Drucker-Prager破坏准则编入其中，当试块模型某单元应力状态满足D-P准则进入塑性状态时，认为该单元失去承载能力发生破坏，将该单元设定为失效单元，在下一加载步迭代计算中将其删除。这样反复迭代计算，直至所有加载步计算结束。

实施例1：

花岗岩质地坚硬致密、强度高，是一种较为常见的工程岩体。选取某工程现场取得的岩芯加工成φ50mm×100mm标准试样2块，在RMT多功能岩石试验机上进行三轴压缩试验(围压10MPa)。由三轴压缩试验得出试样偏应力——轴向应变关系曲线，如图5所示。花岗岩试样破坏如图6所示。可以看出，三轴压缩条件下花岗岩试样的破坏形式较为一致，属于剪切破坏，因此采用考虑剪切效应的D-P准则作为岩石破坏的判断准则是合理的。

对两个试块加载过程的试验数据进行整理，得到每个试样的物理力学参数如表2所示。同时根据室内试验分析，得到岩石试块的粘聚力为27MPa，摩擦角为45°。

表2围压10MPa时三轴压缩试验结果

实施例2：

调用岩石力学试验数值仿真系统，利用图4a、图4b和图4c所示用户自定义交互式窗口输入相应系数，对花岗岩试样破坏过程进行数值模拟。用户自定义交互式窗口参数输入如下：

加载参数设置：轴压P1设定为250000000Pa，围岩为10000000Pa。

试样几何参数设置：试样高H设定为100mm，直径D为50mm。

试块力学参数设置：见表2，粘聚力设定为27000000Pa，摩擦角为45°。

通过岩石力学试验数值仿真系统对试块加载过程进行模拟，得到的模拟结果如图7所示。可以看出，试块加载结束后，整体分为两部分，两者间存在较大的相对滑移，说明试块破坏形式为剪切破坏，这与试验得出的结论基本一致(如图6所示)。同时，在试样顶部取一关键点，通过自动读取改点在模拟过程应力应变关系，得到模拟过程中试样应力-应变曲线如图8所示。可以看出，当模拟加载偏应力达到230MPa左右时，试块失去承载能力，发生整体破坏，模拟过程试样关键点应力-应变曲线与试验过程得到的曲线具有较好的一致性。

因此，通过以上分析可知，岩石力学试验数值仿真系统计算结果与试验结果在试样破坏形式和应力-应变关系方面具有较好的一致性，这也验证了本发明开发的岩石力学试验数值仿真系统的合理性和有效性。

综上所述，本发明利用Python语言对ABAQUS软件进行二次开发，对前处理过程(包括模型尺寸、加载参数的确定等)以及后处理过程(模拟结果的输出、数据的采集)进行编程，建立岩石力学试验过程的参数化数值模拟仿真系统；将考虑剪切效应的Drucker-Prager破坏准则编入vusdfld子程序中，引入失效单元删除算法研究荷载作用下岩石试块变形破坏过程，实现了岩石力学试验过程的数值仿真，并以花岗岩试块三轴压缩试验为例，验证了仿真系统的合理性和有效性；以及本发明的岩石力学试验数值仿真系统可根据用户需要与实验室加载平台进行数据匹配，极大提高了建模和分析效率，同时丰富了力学试验内容。

以上公开的仅为本发明的几个具体实施例，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (6)

1.一种基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，包括：

通过Python语言编程对ABAQUS软件进行前处理二次开发和后处理二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统；

将考虑剪切效应的Drucker-Prager破坏准则编入ABAQUS软件中的vusdfld子程序中，并引入失效单元删除算法模拟荷载作用下岩石试块变形破坏过程。

2.如权利要求1所述的基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，所述前处理二次开发，包括：通过编程实现用户自定义交互式模型尺寸、材料属性、边界条件参数的输入，并通过有限元计算，输出不同参数下的计算结果文件。

3.如权利要求1所述的基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，所述后处理二次开发，包括：通过编程实现对计算结果文件的读取、数据分析和用户指定结果的显示。

4.如权利要求1所述的基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，所述通过Python语言编程对ABAQUS软件进行前处理二次开发和后处理二次开发，建立岩石力学试验过程的参数化数值模拟仿真系统，具体包括：

将加载参数轴压P1和围压P2、试样几何参数高H和直径D、以及试块力学参数密度、弹性模量、泊松比、粘聚力和摩擦角设定为变化参数，利用Python语言编程进行参数化建模；

在参数化建模程序开头添加交互式模块程序代码，使交互式模块程序代码中的参数与试验参数相匹配，当程序调用后首先进入用户自定义交互式窗口，即对数值模拟试验参数进行修改；

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

5.如权利要求4所述的基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，所述vusdfld子程序，用于读取有限元模型迭代计算过程的单元状态；其中，单元状态包括：应力、应变、等效塑性应变。

6.如权利要求5所述的基于ABAQUS的岩石力学试验数值仿真方法，其特征在于，所述将考虑剪切效应的Drucker-Prager破坏准则编入ABAQUS软件中的vusdfld子程序中，并引入失效单元删除算法模拟荷载作用下岩石试块变形破坏过程，具体包括：

当岩石试块有限元模型的某单元应力状态满足Drucker-Prager破坏准则进入塑性状态时，认为该单元失去了承载能力发生破坏，将该单元设定为失效单元，在下一加载步迭代计算中将其删除，反复迭代计算，直至所有加载步计算结束。