Grain-Based Distinct Element Modelling of the Mechanical Behavior of a Single Fracture Embedded in Rock: DECOVALEX-2023 Task G (Benchmark Simulation)

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

입자기반 개별요소모델을 통한 결정질 암석 내 균열의 역학적 거동 모델링: 국제공동연구 DECOVALEX-2023 Task G(Benchmark Simulation)

본 논문에서는 국제공동연구인 DECOVALEX-2023 프로젝트 Task G의 연구 현황과 현재까지 수행된 benchmark 해석 결과를 소개하였다. Task G의 명칭은 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks(SAFENET)'로, 결정질 암반 내 균열의 생성과 성장 메커니즘 및 균열에서 발생하는 열-수리-역학적 복합거동을 해석하기 위한 수치해석기법을 개발하는 데에 목표가 있다. Task G의 첫 번째 연구 테마는 결정질 암석 내 단일 균열의 역학적 거동에 대한 해석해(analytical solution)를 바탕으로 각 연구팀의 수치모델링기법을 개발 및 검증하는 Benchmark 해석이다. 본 연구에서는 3차원 입자기반 개별요소모델을 이용하여 단일 균열을 포함한 암석의 역학적 거동 특성을 모델링하고자 하였다. 이 모델에서는 상호독립적으로 거동하는 개별입자의 집합체를 통해 암석의 구조적 특징을 모사하고, 입자와 입자간 접촉에서 발생하는 역학적 거동을 개별요소해석모델인 3DEC을 통해 계산하게 된다. 해석 결과, 도메인의 경계응력으로 인해 균열에 유도되는 수직응력과 전단응력 수준은 변위 구속과 응력 재배치로 인해 이론적인 수치보다 낮게 나타났다. 그러나 수치모델에서 계산된 수직변위와 전단변위는 실제 균열의 유도 응력을 통해 추정된 해석해와 비교할 때 상당히 유사한 결과를 보였으며 균열의 응력-변위 관계를 합리적으로 재현할 수 있음을 확인하였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 다양한 조건의 실내시험에 적용하여 타당성을 검증할 예정이다.