• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.4
• /
• pp.1541-1549
• /
• 2013

This study investigates the response of structures to tunnelling-induced ground movements in sand ground, varying tunnel excavation condition (tunnel depth and diameter), tunnel construction condition (ground loss), ground condition (loose sand and dense sand). Four-story block-bearing structures have been used because the structueres can easily be characterized of the extent of dmages with crack size and distribution. Numerical parametric studies have been used to investigae of the response of structures to varying tunnelling conditions. Numerical analysis has been conducted using Discrete Element Method (DEM) to have real cracks when the shear and tensile stress exceed the maximum shear and tensile strength. The results of structure responses from various parametric studies have been integrated to consider tunnel excavation condition, tunnel construction condition, and ground condition and provided as a relationship chart. Using the chart, the response of structures to tunnelling can easily be evaluated in practice in sand ground.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.11
• /
• pp.39-45
• /
• 2011

This study presents the thermal behaviors of glass beads-rubber mixtures depending on the volumetric fraction of each constituent and relative size between them. The transient plane source method is used to measure the effective thermal conductivity of mixtures. The discrete element method (DEM) and the thermal network model are integrated to investigate the particle-scale mechanism of heat transfer in granular packings. Results show that 1) the effective thermal conductivity decreases as the rubber fraction increases, and 2) the relative size between two solid particles dominates the spatial configuration of inter-particle contact condition that in tum determines the majority of heat propagation path through particle contacts. For the mixtures whose volumetric fraction of rubber is identical, the less conductive materials (e.g., rubber particles) with a large size facilitate heat transfer in granular materials. The experimental results and particle-scale observation highlight that the thermal conduction behavior is dominated not only by the volumetric fraction but also the spatial configuration of each constituent.

• Geomechanics and Engineering
• /
• v.17 no.1
• /
• pp.19-30
• /
• 2019

Recent research on the behavior of silty sand tends to advocate the use of equivalent skeleton void ratio to characterize the density state of this type of soil. This paper presents an investigation to explore the physical meaning of the equivalent skeleton void ratio by means of DEM simulations for assemblies of coarse and fine particles under biaxial shear. The simulations reveal that the distribution pattern of fine particles in the soil skeleton plays a crucial role in the overall macroscopic response: The contractive response observed at the macro scale is mainly caused by the movement of fine particles out of the force chains whereas the dilative response is mainly associated with the migration of fine particles into the force chains. In an assembly of coarse and fine particles, neither all of the fine particles nor all of the coarse ones participate in the force chains to carry the external loads, and therefore a more reasonable definition for equivalent skeleton void ratio is put forward in which a new parameter d is introduced to take into account the fraction of coarse particles absent from the force chains.

• Journal of the Korea Convergence Society
• /
• v.12 no.1
• /
• pp.185-191
• /
• 2021

A powder metallurgy process with granules is the manufacturing technology that can achieve higher-density sintered parts than conventional powder metallurgy processes. However, there is a disadvantage in that the production cost increases significantly due to the additional granulation step. High granule productivity must be guaranteed for affordable material costs in this manufacturing technology. This paper performed a series of scattering, collision, and adhesion simulations of agglomerated powders to investigate the characteristics of granulation process affecting the manufacturing yield rate. The results of this simulation-based convergence study can contribute to improving productivity in the metal powder spray granulation process.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.2
• /
• pp.93-103
• /
• 2023

Warships widely spread numerous chaffs using a blast, which form chaff clouds that create false radar cross-sections to deceive enemy radars. In this study, we established a numerical framework based on a one-way coupling of computational fluid dynamics and discrete element method to simulate the spatiotemporal distribution of chaff clouds for warships in the air. Using the framework, we investigated the effects of wind, initial chaff cartridge angle, and blast pressure on the distribution of chaff clouds. We observed three phases for the chaff cloud diffusion: radial diffusion by the explosion, omnidirectional diffusion by turbulence and collision, and gravity-induced diffusion by the difference in the fall speed. The wind moved the average position of the chaff clouds, and the diffusion due to drag force did not occur. The direction of radial diffusion by the explosion depended on the initial angle of the cartridge, and a more vertical angle led to a wider distribution of the chaffs. As the blast pressure increased, the chaff clouds spread out more widely, but the distribution difference in the direction of gravity was not significant.

• Geomechanics and Engineering
• /
• v.12 no.6
• /
• pp.1003-1020
• /
• 2017

Estimation of groundwater inflow into underground opening is of critical importance for the design and construction of underground structures. Groundwater inflow into a pilot underground storage facility in China was estimated using analytical equations, numerical modeling and field measurement. The applicability of analytical and numerical methods was examined by comparing the estimated and measured results. Field geological investigation indicated that in local scale the high groundwater inflows are associated with the appearance of open joints, fractured zone or dykes induced by shear and/or tensile tectonic stresses. It was found that 8 groundwater inflow spots with high inflow rates account for about 82% of the total rate for the 9 caverns. On the prediction of the magnitude of groundwater inflow rate, it was found that could both (Finite Element Method) FEM and (Discrete Element Method) DEM perform better than analytical equations, due to the fact that in analytical equations simplified assumptions were adopted. However, on the prediction of the spatial distribution estimation of groundwater inflow, both analytical and numerical methods failed to predict at the present state. Nevertheless, numerical simulations would prevail over analytical methods to predict the distribution if more details in the simulations were taken into consideration.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.6
• /
• pp.35-40
• /
• 2020

The ballast layers play a key role in distributing and supporting a trainload. On the other hand, it settles down by dynamic train loading due to large void ratios. Consequently, it requires continuous maintenance. In this paper, ballast layers mixed with three types of specially shaped artificial ballast (AB) (Rectangular, Tetrapod, Hexagonal) were modeled by using a two dimensional DEM (Discrete Element Method). Repeated loading tests were performed to evaluate the settlement behavior of the ballast layers. The smallest settlement was observed in the case of the ballast layer mixed with Tetrapod AB than in other cases, according to an analysis of the force transfer routes. In addition, contact force analysis showed that the Tetrapod AB, which has a concave shape, could easily make small and multi-channel force-transfer routes. This means that the stress in the ballast layer by the train loading transferred through the sleeper uniformly was distributed well by the AB. Therefore, the settlement of the ballast layer mixed with the concave-shaped Tetrapod AB could be reduced effectively under a repeated train loading. The effects of a decrease in settlement of the ballast layer highlight the possibility of a maintenance-free ballasted track.

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

• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.51-62
• /
• 1993

A computational program ［TDET］ of the particle transport equation is developed on radiation shielding problem in two-dimensional cartesian geometry based on the discrete element method. Not like the ordinary discrete ordinates method, the quadrature set of angles is not fixed but steered by the spatially dependent angular fluxes. The angular dependence of the scattering source term in the particle transport equation is described by series expansion in spherical harmonics, and the energy dependence of the particles is considered as well. Three different benchmark tests are made for verification of TDET : For the ray effect analysis on a square absorber with a flat isotropic source, the results of TDET calculation are quite well conformed to those of MORSE-CG calculation while TDET ameliorates the ray effect more effectively than S$_{N}$ calculation. In the analysis of the streaming leakage through a narrow vacuum duct in a shield, TDET shows conspicuous and remarkable results of streaming leakage through the duct as well as MORSE-CG does, and quite better than S$_{N}$ calculation. In a realistic reactor shielding situation which treats in two cases of the isotropic scattering and of linearly anisotropic scattering with two groups of energy, TDET calculations show local ray effect between neighboring meshes compared with S$_{N}$ calculations in which the ray effect extends broadly over several meshes.eshes.