• Tunnel and Underground Space
• /
• v.21 no.2
• /
• pp.128-137
• /
• 2011

The shut-in pressure calculated in common hydrofracturing test for vertical borehole equals generally to the minimum horizontal principal stress, so it should be considered as an essential parameter for determining the in-situ stress regime around the rock mass. It shows usually an ambiguous value in pressure-time history curves, however, because of the relationship between the behavior of hydraulic fractures and the condition of remote stress regime. In this study, a series of numerical analyses have been carried out to compare several methods for determining the shut-in pressure during hydrofracturing. The hydraulic-mechanical coupling has been applied to numerical analysis for simulating the fracture propagation by hydraulic pressure, and the different discontinuity geometry has been considered in numerical models to examine the effect of numerical element shape on fracture propagation pattern. From the numerical simulations with the four different discontinuity geometries, it was revealed that the shut-in pressure obtained from graphical methods rather than statistical method was relatively small. Consequently a care should be taken in selecting a method for determining the shut-in pressure when a stress anomaly around borehole and a fracture propagation with complicate mechanism are considered.

• Tunnel and Underground Space
• /
• v.21 no.1
• /
• pp.66-81
• /
• 2011

In this study, it was aimed to develop a thermal-hydraulic-mechanical coupled fracture mechanics code that models a fracture initiation, propagation and failure of underground rock mass due to thermal and hydraulic loadings. The development was based on a 2D FRACOD (Shen & Stephasson, 1993), and newly developed T-M and H-M coupled analysis modules were implemented into it. T-M coupling in FRACOD employed a fictitious heat source and time-marching method, and explicit iteration method was used in H-M coupling. The validity of developed coupled modules was verified by the comparison with the analytical result, and its applicability to the fracture initiation and propagation behavior due to temperature changes and hydraulic fracturing was confirmed by test simulations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.31-37
• /
• 2017

In this paper, computation results of reducible modeling, stress recovery and energy release rate were compared with the results of VABS, Virtual Crack Closure Technique. The result of stress recovery analysis for 1-D model including the stiffness matrix is compared with stress results of three-dimensional 3-D FEM. Energy release rate of composite beam with longitudinal cracks is calculated and compare verifications of numerical analysis results of 3-D FEM and VABS. The procedure of calculating energy release rate through dimensional reduction and stress recovery is intended to be efficient and be utilized in the life-cycle of high-altitude uav's wing, wind blades and tilt rotor blade.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.245-247
• /
• 1999

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

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.

• Journal of the Korean Recycled Construction Resources Institute
• /
• v.4 no.1
• /
• pp.10-18
• /
• 2016

RC(Reinfored Concrete) structures exposed to carbonation in urban city have durability degradation with extended service life and cracks in concrete causes a local accelerated carbonation. In the present work, crack effect on carbonation depth is investigated and the service life of RC structure is evaluated considering cracks from early age and time-dependent cracks based on the previous field investigation. DFP(Durability Failure Probability), safety index, and the related service life are calculated considering the time to crack width reaches to maximum crack width(0.3mm). The results with time effect on crack width show lower DFP and longer service life, which seems to be reasonable compared with conservative results from crack effect from initial stage. Furthermore, crack effect is evaluated to be insignificant on DFP and service life. The technique with time-dependent crack effect on carbonation can be effectively used for RC structure containing cracking in use.

• Computational Structural Engineering
• /
• v.7 no.2
• /
• pp.111-120
• /
• 1994

Cracks in rubber bonded to a rigid material such as steel are analyzed with the aid of a mixed finite element technique. Firstly the weak form is derived for finite element analysis of an incompressible material, and the Mooney-Rivlin form is assumed for the constitutive modeling of rubber. The numerical results from finite element analysis is examined to confirm the accuracy and convergence of solution by way of comparison to other numerical results. The interpretation of the J-integral for large elastic deformation as the energy release rate is confirmed, and the J-integral is calculated for varing crack length. The crack growth stability is discussed using the result of finite element analysis.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.4 no.4
• /
• pp.37-51
• /
• 2000

이 연구는 지진 시 철근콘크리트 교각의 비탄성 거동 및 연성능력을 해석적으로 파악하는데 그 목적이 있다. 재료적 비선형성에 대해서는 균열 콘크리트에 대한 인장, 압축, 전단모델과 콘크리트 속에 있는 철근 모델을 조합하여 고려하였다. 이에 대한 콘크리트의 균열 모델로서의 분산균열모델을 사용하였다. 두께가 서로 다른 부재간의 접합부에 단면강성이 급변하기 때문에 생기는 국소적인 불연속변형을 고려하기 위한 경계면 요소를 도입하였다. 또한, 축방향철근의 유무 및 그 양 등에 따른 구속효과를 적절히 표현할 수 있는 해석 모델을 개발하였다. 본 연구에서는 철근콘크리트 교각의 비탄성 거동 및 연성 능력의 파악을 위해 제안한 해석기법을 신뢰성 있는 연구자의 실험결과와 비교하여 그 타당성을 검증하였다.