• 터널과지하공간
• /
• 제24권3호
• /
• pp.243-254
• /
• 2014

본 연구에서는 GBM-UDEC(grain based model combined with Universal Distinct Element Code) 모델을 통해 암석의 역학적 거동을 모사하기 위한 수치해석기법을 소개하였다. 이를 적용하여 암석의 광물학적 구조를 다각형 입자의 집합체로 표현하고, 압축하중 하에서 암석의 파괴 특성 및 인장균열의 전파 양상를 살펴보았다. 제시된 수치해석모델은 단축압축강도시험 및 간접인장강도시험을 통해 관찰되는 암석의 역학적 거동 특성을 합리적으로 모사할 수 있는 것으로 나타났다.

• 터널과지하공간
• /
• 제31권6호
• /
• pp.593-609
• /
• 2021

본 연구에서는 입자기반 개별요소모델(grain-based distinct element model, GBDEM)을 이용하여 결정질 암석 내 포함된 균열의 열-역학적 거동을 평가할 수 있는 수치해석기법을 제시하고 열에 의한 균열의 미끄러짐 거동을 해석하였다. 이는 DECOVALEX-2023 프로젝트 Task G의 일환으로 수행된 벤치마크 모델링 연구로, Task G는 결정질 암반 내 균열의 열-수리-역학적 복합거동을 해석하기 위한 수치해석기법을 개발하는 데에 목표가 있다. 여기에서는 Voronoi diagram을 이용하여 다면체 개별입자의 집합체로서 해석모델을 생성하고, 입자 및 입자간 접촉에서 발생하는 열-역학적 거동을 개별요소프로그램인 3DEC을 통해 해석하였다. 암석 시험편의 탄성거동을 재현하기 위하여 등가연속체 개념을 적용하여 입자와 접촉의 미시물성을 산정하였으며, 균열에 상응하는 접촉에는 Coulomb slip model을 부여하여 인장강도와 전단강도를 갖는 불연속면을 모사하였다. 경계응력과 열응력에 의한 균열의 거동을 수치적으로 모델링하였으며, 경계조건에 따라 균열의 미끄러짐이 발생하는 열-역학적 메커니즘을 정량적으로 분석하였다. 해석 결과, 본 연구에서 제시한 해석모델이 암석 내 열팽창과 열응력의 증가, 균열 응력과 변위, 경계조건의 영향 등을 합리적으로 재현하고 있음을 확인하였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 실내실험에 적용하여 타당성을 검증할 예정이다.

• 터널과지하공간
• /
• 제31권4호
• /
• pp.270-288
• /
• 2021

본 연구에서는 입자기반 개별요소모델(grain-based distinct element model, GBDEM)을 이용하여 암석 균열의 역학적, 수리적 거동을 평가할 수 있는 수치해석기법을 제시하고 해석해와의 비교를 통해 검증하였다. 이는 DECOVALEX-2023 프로젝트 Task G의 일환으로 수행된 벤치마크 모델링 연구로, Task G는 결정질 암반 내 균열의 열-수리-역학적 복합거동을 해석하기 위한 수치해석기법을 개발하는 데에 목표가 있다. 본 연구에서는 사면체 개별 입자들을 이용하여 해석모델을 생성하고 3DEC을 이용하여 입자와 접촉에서의 거동을 해석하였다. 이 과정에서 등가연속체 개념을 적용해 입자기반모델의 미시물성을 산정할 수 있는 새로운 기법을 제시하였다. 한편, 균열 경사각과 거칠기, 경계응력조건 및 압력 조건에 따른 해석을 실시하여 각 해석조건이 균열의 수직, 전단방향 거동에 미치는 영향을 살펴보았다. 해석 결과, 제안된 수치모델은 경계응력에 따른 균열의 미끄러짐(fracture slip)과 유체 압력에 따른 균열의 개방(fracture opening), 균열 경사에 따른 응력 분포, 거칠기로 인한 전단변위의 구속 등을 합리적으로 재현하고 있음을 확인하였다. 수치해석을 통해 계산된 균열의 수직방향, 전단방향 변위는 모두 해석해를 통해 계산된 값과 거의 일치하는 결과를 보였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 다양한 조건의 실내시험에 적용하여 타당성을 검증할 예정이다.

• 터널과지하공간
• /
• 제32권6호
• /
• pp.568-585
• /
• 2022

본 연구에서는 3차원 입자기반 개별요소모델(grain-based distinct element model, GBDEM)을 이용하여 암석 균열의 열에 의한 미끄러짐 거동을 해석하였다. 이는 DECOVALEX-2023 프로젝트 Task G에 참여하여 수행한 연구로, 해석대상은 한국건설기술연구원에서 수행된 saw-cut 균열 시료에 대한 열-역학적 하중 재하 실험 결과이다. 여기에서는 암석 시험편을 Voronoi 다면체의 집합체로 모델링하고, 개별요소법 코드인 3DEC을 통해 입자와 입자 간 경계면, 내부에 포함된 균열에서의 열-역학적 연계거동을 해석하였다. 주요 해석내용은 가열로 인한 암석 표면의 온도 분포, 열응력의 증가에 따른 주응력 변화, 균열의 전단변위와 수직변위이다. 해석 결과, 상기 수치모델은 실내실험에서 관찰된 열전달과 열손실 특성, 열에 의한 균열의 점진적 전단파괴 프로세스, 변위의 제한으로 인한 열응력의 증가 등을 합리적 수준에서 재현하고 있는 것으로 나타났다. 그러나 가열에 의한 전단파괴 시점, 열응력 증분과 변위 크기 등에서는 다소 차이를 보였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류 및 협력을 통해 지속적으로 개선, 검증할 예정이다.

• 터널과지하공간
• /
• 제26권3호
• /
• pp.235-252
• /
• 2016

본 연구에서는 다각형 입자 기반 개별요소모델을 이용하여 실험실 스케일에서 등방성, 횡등방성 암석의 거동과 점진적 파괴 과정을 모델링할 수 있는 시뮬레이션 기법을 구축하였다. 가압에 따른 미세균열의 개시와 성장 과정을 모니터링할 수 있는 기법을 제안하였으며, 이를 통해 전단균열과 인장균열의 개시와 성장이 암석의 점진적 파괴 과정에 미치는 영향을 살펴보았다. 다각형 입자기반 개별요소모델의 거동 및 미세균열의 발생 양상은 실험실에서 관찰되는 암석의 일반적인 특징과 상당 부분 일치하는 것으로 나타났으며, 이를 통해 상기 모델이 암석의 역학적 거동을 합리적인 수준에서 재현할 수 있음을 확인하였다. 다각형 입자기반 개별요소모델에 대한 기초연구로서 접촉면의 미시변수와 시료의 거시물성 간의 상관관계를 살펴보았으며, 미시변수를 조정함으로써 다양한 암석의 강도와 변형 특성을 재현하였다. 한편, 상기 모델을 횡등방성 암석을 모사하기 위한 방법론을 제시하였으며, 이를 국내 횡등방성 암석인 아산편마암에 적용하여 근소한 오차 범위 내에서 실내시험 결과를 재현하였다.

• 터널과지하공간
• /
• 제30권6호
• /
• pp.573-590
• /
• 2020

본 논문에서는 국제공동연구인 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에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 다양한 조건의 실내시험에 적용하여 타당성을 검증할 예정이다.

• 터널과지하공간
• /
• 제33권3호
• /
• pp.189-207
• /
• 2023

In the present study, the thermoshearing experiment on a rough rock fracture were modeled using a three-dimensional grain-based distinct element model (GBDEM). The experiment was conducted by the Korea Institute of Construction Technology to investigate the progressive shear failure of fracture under the influence of thermal stress in a critical stress state. The numerical model employs an assembly of multiple polyhedral grains and their interfaces to represent the rock sample, and calculates the coupled thermo-mechanical behavior of the grains (blocks) and the interfaces (contacts) using 3DEC, a DEM code. The primary focus was on simulating the temperature evolution, generation of thermal stress, and shear and normal displacements of the fracture. Two fracture models, namely the mated fracture model and the unmated fracture model, were constructed based on the degree of surface matedness, and their respective behaviors were compared and analyzed. By leveraging the advantage of the DEM, the contact area between the fracture surfaces was continuously monitored during the simulation, enabling an examination of its influence on shear behavior. The numerical results demonstrated distinct differences depending on the degree of the surface matedness at the initial stage. In the mated fracture model, where the surfaces were in almost full contact, the characteristic stages of peak stress and residual stress commonly observed in shear behavior of natural rock joints were reasonably replicated, despite exhibiting discrepancies with the experimental results. The analysis of contact area variation over time confirmed that our numerical model effectively simulated the abrupt normal dilation and shear slip, stress softening phenomenon, and transition to the residual state that occur during the peak stress stage. The unmated fracture model, which closely resembled the experimental specimen, showed qualitative agreement with the experimental observations, including heat transfer characteristics, the progressive shear failure process induced by heating, and the increase in thermal stress. However, there were some mismatches between the numerical and experimental results regarding the onset of fracture slip and the magnitudes of fracture stress and displacement. This research was conducted as part of DECOVALEX-2023 Task G, and we expect the numerical model to be enhanced through continued collaboration with other research teams and validated in further studies.