• 터널과지하공간
• /
• 제31권5호
• /
• pp.309-332
• /
• 2021

고준위방사성폐기물 심층처분장 내 천연방벽은 방사성핵종의 누출을 방지 및 지연할 수 있는 수리학적 특성을 갖춰야 한다. 결정질 암반의 경우 불연속면에 의해 수리학적 성능이 결정되기 때문에, 불연속면의 수리-역학적 복합거동에 대한 자세한 모사가 필요하다. 불연속체 기반 해석기법은 불연속면의 생성, 전파, 변형, 미끄러짐과 같은 복잡한 거동을 구현할 수 있어 결정질 암반 모사에 적합하다. 본 연구에서는 불연속면에서의 수리-역학 복합거동에 초점을 맞추어, UDEC, 3DEC, PFC, DDA, FRACOD, TOUGH-UDEC과 같은 상용화된 불연속체 기반 수리-역학 복합거동 해석기법을 조사하였다. 블록 기반 불연속체 해석기법의 경우 주로 불연속면 상에서 진행되는 유체 유동을 바탕으로 수리-역학 복합거동을 해석하였고, 그중 일부는 다른 수리학적 해석기법과의 결합을 통하여 모델 전체에 대한 수리-역학적 복합거동을 제공하였다. 입자 기반 불연속체 해석기법의 경우에는 입자 사이로 흐르는 유체를 반영하여 불연속체 모델 전체에 해당하는 수리-역학적 복합거동 모사가 가능하다. 현재까지 상용화된 불연속체 기반 복합거동 해석기법은 2차원 해석만 제공하거나, 수리학적 해석 성능이 떨어지고, 불연속면에서의 유체 유동만 고려되거나, 자세한 수리학적 해석을 지원하지 않는 등의 한계점이 있어 고준위방사성폐기물 심층처분시스템의 정확한 수리-역학 모델링에는 적합하지 않을 수 있다. 본 기술보고에서 검토한 다양한 해석기법들의 장단점을 참고하여 향후 처분시스템을 정확하고 자세하게 모사할 수 있는 불연속체 기반 수리-역학 복합거동 해석기법의 개발이 필요하다.

• 터널과지하공간
• /
• 제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.