• 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.

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.568-585
• /
• 2022

In the present study, we proposed a numerical method for simulating thermally induced fracture slip using a grain-based distinct element model (GBDEM). As a part of DECOVALEX-2023, the thermo-mechanical loading test on a saw-cut rock fracture conducted at the Korea Institute of Civil Engineering and Building Technology was simulated. In the numerical model, the rock sample including a saw-cut fracture was represented as a group of random Voronoi polyhedra. Then, the coupled thermo-mechanical behavior of grains and their interfaces was calculated using 3DEC. The key concerns focused on the temperature evolution, thermally induced principal stress increment, and fracture normal and shear displacements under thermo-mechanical loading. The comparisons between laboratory experimental results and the numerical results revealed that the numerical model reasonably captured the heat transfer and heat loss characteristics of the rock specimen, the horizontal stress increment due to constrained displacement, and the progressive shear failure of the fracture. However, the onset of the fracture slip and the magnitudes of stress increment and fracture displacement showed discrepancies between the numerical and experimental results. 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.

• Journal of Korean Society of societal Security
• /
• v.3 no.1
• /
• pp.43-50
• /
• 2010

Sliding accidents on the road have a high percentage by road freezing, especially, they often have appeared at bridges and Tunnel of freezing areas. Thus, the stability of road operations is enhanced by preventing partial freezing phenomenon. According to the geothermal snow melting system analysis, a pattern of thermal conductivity is found out; pavement materials of concrete and asphalt where the system is buried. The heat transfer simulation is essential when the geothermal snow melting system is applied according to heating exchanger pipe placed in the lower pavements. The model tests are conducted on low temperature in freezer using the manufactured test model which is equal to pavement materials. Many variables are discovered from numerical analyses under the same conditions with model test.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.12
• /
• pp.679-684
• /
• 2017

Friction Stir Welding is a metal welding technique, in which friction heat between a welding tool and a welding material is used to weld parts at temperatures below the melting point of a material. In this study, the temperature and velocity changes in a magnesium alloy (AZ31) during the welding process were analyzed by computational flow dynamics technique while welding the material using a friction stir welding technique. For the analysis, the modeling and analysis were carried out using Fluent as a fluid analysis tool. First, the welding material was assumed to be a temperature-dependent Newtonian fluid with high viscosity, and the rotation region and the stationary region were simulated separately to consider the rotational flow generated by the rotation of the welding tool having a helical groove. The interface between the welding tool and welding material was given the friction and slip boundary conditions and the heat transfer effect to the welding tool was considered. Overall, the velocity and temperature characteristics of the welded material according to time can be understood from the results of transient analysis through the above flow analysis modeling.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.22 no.1
• /
• pp.190-198
• /
• 2018

In recent years, non-welded building hardware has been installed by bolt assembly is used. The non-welded building hardware method can reduce accidents caused by welding, and can be constructed by bolt assembly, which can reduce labor costs and shorten the construction period. However, there is a need for a method to compensate for the occurrence of buckling at the time of construction. The purpose of this study is to evaluate the behavior of joints between steel pipe and fastener and to evaluate the behavior of joints of non-welded and welded hardware frame. As a result, it was found that the foundation steel structure without welded joints was deformed to a rotation angle of member much larger than the allowable interlayer displacement angle 0.01 to 0.02 required according to the seismic load rating in the seismic load resistance system.