• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.105-115
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.14 no.4
• /
• pp.293-308
• /
• 2012

The shear behavior at the particle/surface interface such as rock joint can determine the mechanical behavior of whole structure. Therefore, a fundamental understanding of the mechanisms governing its behavior and accurately estimation of the interface strength is essential. In this paper, PFC, a numerical analysis program of discrete element method was used to investigate the effects of the surface roughness on interface strength. The surface roughness was characterized by smooth, intermediate, and rough surface, respectively. In order to investigate the effects of particle shape and crushing on particle/surface interface behavior, one ball, clump, and cluster models were created and their results were compared. The shape of particle was characterized by circle, triangle, square, and rectangle, respectively. The results showed that as the surface roughness increases, interface strength and friction angle increase and the void ratio increases. The one ball model with smooth surface shows lower interface strength and friction angle than the clump model with irregular surface. In addition, a cluster model has less interface strength and friction angle than the clump model. The failure envelope of the cluster model shows non-linear characteristic. From these findings, it is verified that the surface roughness and particle shape effect on the particle/surface interface shear behavior.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.14 no.4
• /
• pp.321-336
• /
• 2012

The shear behavior at the particle/surface interface such as rock joint can determine the mechanical behavior of whole structure. Therefore, a fundamental understanding of the mechanisms governing its behavior and accurately estimation of the interface strength is essential. In this paper, PFC, a numerical analysis program of discrete element method was used to investigate the effects of the surface roughness crushing on interface strength. The surface roughness was characterized by smooth, intermediate, and rough surface, respectively. Particle shape was classified into one ball model of circular shape and 3 ball model of triangular shape. The surface shape was modelled by wall model of non-crushing surface and ball model of crushing surface. The results showed that as the bonding strength of ball model decreases, lower interface strength is induced. After the surface roughness crushing was occurred, the interface strength tended to converge and higher bonding strength induced lower surface roughness crushing. Higher friction angle was induced in wall model and higher surface roughness induced the higher friction angle. From these findings, it is verified that the surface roughness and surface roughness crushing effect on the particle/surface interface shear behavior.

• Tunnel and Underground Space
• /
• v.29 no.6
• /
• pp.456-467
• /
• 2019

In this study, the discrete element method was used to simulate the excavation of spoke-type shield TBM. The horizontal stress coefficient was used for the ground to simulate the increase of the horizontal stress according to the depth, and the driving conditions were set based on the torque generated from the cutterhead of the TBM to excavate within the operating range. That is, when the value of the torque generated at the cutterhead exceeds the given operating condition, the speed of excavation is constantly reduced, and conversely, the method of increasing the speed of excavation is considered. The change speed of the excavation was given the minimum change requirement in consideration of the driver's review time, and the change was possible according to the excavation conditions. In order to use these conditions, the user-subroutine was considered separately, and the results show that the DEM model were able to analyze the excavation within the considered operating range.

• Steel and Composite Structures
• /
• v.10 no.2
• /
• pp.129-149
• /
• 2010

The objective of this study is to formulate a general 3D material-structural analysis framework for the thermomechanical behavior of steel-concrete structures in a fire environment. The proposed analysis framework consists of three sequential modeling parts: fire dynamics simulation, heat transfer analysis, and a thermomechanical stress analysis of the structure. The first modeling part consists of applying the NIST (National Institute of Standards and Technology) Fire Dynamics Simulator (FDS) where coupled CFD (Computational Fluid Dynamics) with thermodynamics are combined to realistically model the fire progression within the steel-concrete structure. The goal is to generate the spatial-temporal (ST) solution variables (temperature, heat flux) on the surfaces of the structure. The FDS-ST solutions are generated in a discrete form. Continuous FDS-ST approximations are then developed to represent the temperature or heat-flux at any given time or point within the structure. An extensive numerical study is carried out to examine the best ST approximation functions that strike a balance between accuracy and simplicity. The second modeling part consists of a finite-element (FE) transient heat analysis of the structure using the continuous FDS-ST surface variables as prescribed thermal boundary conditions. The third modeling part is a thermomechanical FE structural analysis using both nonlinear material and geometry. The temperature history from the second modeling part is used at all nodal points. The ABAQUS (2003) FE code is used with external user subroutines for the second and third simulation parts in order to describe the specific heat temperature nonlinear dependency that drastically affects the transient thermal solution especially for concrete materials. User subroutines are also developed to apply the continuous FDS-ST surface nodal boundary conditions in the transient heat FE analysis. The proposed modeling framework is applied to predict the temperature and deflection of the well-documented third Cardington fire test.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.9 no.2
• /
• pp.133-141
• /
• 2007

The excavation damaged zone (EDZ) is an area around an excavation where in situ rock mass properties, stress condition, displacement, groundwater flow conditions have been altered due to the processes induced by the excavation. Various studies have been carried out on EDZ, but most studies have focused on the mechanical bahavior of EDZ by in situ experiment. Even though the EDZ could potentially form a high permeable pathway of groundwater flow, only a few studies were performed on the analysis of groundwater flow in EDZ. In this study, the 'hydraulic EDZ' was defined as the rock zone adjacent to the excavation where the hydraulic aperture has been changed due to the excavation by using H-M coupling analysis. Fundamental principles of distinct element method (DEM) were used in the analysis. In the same groundwater level, the behavior of hydraulic aperture near the cavern was analyzed for different stress ratios, initial apertures, fracture angles and fracture spacings by using a two-dimensional DEM program. We evaluate the excavation induced hydraulic aperture change. Using the results of the study, hydraulic EDZ was defined as an elliptical shape model perpendicular to the joint.

• Transactions of Materials Processing
• /
• v.18 no.7
• /
• pp.556-564
• /
• 2009

Flexible forming process for sheet material using reconfigurable die is introduced based on numerical simulation. In general, this flexible forming process using the reconfigurable die has been utilized for manufacturing of curved thick plates used for hull structures, architectural structures and so on. In this study, numerical simulation of sheet metal forming process is carried out by using flexible dies model instead of conventional matched die set. The numerical simulation and experimental verification for sheet metal forming process using a flexible forming machine that is more suitable for thick plate forming process are carried out to confirm the appropriateness of the simulation process. As an elastic cushion, urethane pads are utilized using hyperelastic material model in the simulation for smoothing the forming surface which is discrete due to characteristics of the flexile die. In the flexible forming process for sheet metal, effect of a blank holder is also investigated according to blank holding methods. Formability in view of occurrence of dimples is compared with regard to the various punch sizes. Consequently, it is confirmed that the flexible forming for sheet material using urethane pad has enough capability and feasibility for manufacturing of smoothly curved surface instead of conventional die forming method.

• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.453-463
• /
• 2018

Three points bending flexural test was modelled numerically to study the crack propagation in the pre-cracked beams. The pre-existing double internal cracks inside the beam models were considered to investigate the crack propagation and coalescence paths within the modelled samples. Notch configuration effects on the failure stress were considered too. This numerical analysis shown that the propagation of wing cracks emanating from the tips of the pre-existing internal cracks caused the final breaking of beams specimens. It was also shown that when two notches were overlapped, they both mobilized in the failure process and the failure stress was decreased when the notches were located in centre line. However, the failure stress was increased by increasing the bridge area angle. Finally, it was shown that in all cases, there were good agreements between the discrete element method results and, the other numerical and experimental results. In this research, it is tried to improve the understanding of the crack propagation and crack coalescence phenomena in brittle materials which is of paramount importance in the stability analyses of rock and concrete structures, such as the underground openings, rock slopes and tunnel construction.

• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.97-108
• /
• 2020

The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.