• Journal of the Korean Geotechnical Society
• /
• v.33 no.11
• /
• pp.83-95
• /
• 2017

This study examined the magnitude and distribution of earth pressure on the excavation wall in jointed rock mass by considering different groundwater conditions under various rock types, joint inclination angles, and earth pressure coefficients. Based on a physical model test (Son and Park, 2014), extended studies were conducted considering rock-structure interactions based on the discrete element method, which can consider the joints characteristics of rock mass. The results showed that the earth pressure was highly influenced by the groundwater condition as well as the rock type, joint inclination angle, and earth pressure coefficient. The results were also compared with Peck's earth pressure for soil ground, and clearly showed that the earth pressure in jointed rock mass can be greatly different from that in soil ground.

• Geomechanics and Engineering
• /
• v.22 no.5
• /
• pp.385-396
• /
• 2020

A uniaxial compression test was performed to analyse the mechanical properties and macroscale and mesoscale failure mechanisms of sandstone with pyrite concretions. The effect of the pyrite concretions on the evolution of macroscale cracks in the sandstone was further investigated through numerical simulations with Particle Flow Code in 2D (PFC2D). The results revealed that pyrite concretions substantially influence the mechanical properties and macroscale and mesoscale failure characteristics of sandstone. During the initial loading stage, significant stress concentrations occurred around the edges of the pyrite concretion accompanied by the preferential generation of cracks. Meanwhile, the events and cumulative energy counts of the acoustic emission (AE) signal increased rapidly because of friction sliding between the concretion and sandstone matrix. As the axial stress increased, the degree of the stress concentration remained relatively unchanged around the edges of the concretions. The cracks continued growing rapidly around the edges of the concretions and gradually expanded toward the centre of the sample. During this stage, the AE events and cumulative energy counts increased quite slowly. As the axial stress approached the peak strength of the sandstone, the cracks that developed around the edges of the concretion started to merge with cracks that propagated at the top-left and bottom-right corners of the sample. This crack evolution ultimately resulted in the shear failure of the sandstone sample around the edges of the pyrite concretions.

• Korean Journal of Materials Research
• /
• v.31 no.5
• /
• pp.286-295
• /
• 2021

This study attempts to find optimal conditions of the friction coefficient using a discrete element method (DEM) simulation with various friction coefficient conditions and three different grinding media with various ball sizes in a traditional ball mill (TBM). Using ball motion of the DEM simulation are obtained using the optimal friction coefficient compared with actual motion; photographs are taken by the digital camera and the snapshot images are analyzed. In the simulation, the rotation speed of the mill, the materials and velocity of the grinding media, and the friction coefficient between the balls and the wall of the pot are fixed as the actual experimental conditions. We observe the velocity according to the friction coefficient from the DEM simulation. The friction coefficient is found to increase with the velocity. Milling experiments using a traditional ball mill with the same experimental conditions as those of the DEM simulation are conducted to verify the simulated results. In addition, particle morphology change of copper powder is investigated and analyzed using scanning electron microscopy (SEM) for the milling experiment.

• Structural Engineering and Mechanics
• /
• v.69 no.1
• /
• pp.11-20
• /
• 2019

The biaxial failure mechanism of transversally bedding concrete layers was numerically simulated using a sophisticated two-dimensional discrete element method (DEM) implemented in the particle flow code (PFC2D). This numerical modelling code was first calibrated by uniaxial compression and Brazilian testing results to ensure the conformity of the simulated numerical model's response. Secondly, 21 rectangular models with dimension of $54mm{\times}108mm$ were built. Each model contains two transversely bedding layers. The first bedding layer has low mechanical properties, less than mechanical properties of intact material, and second bedding layer has high mechanical properties, more than mechanical properties of intact material. The angle of first bedding layer, with weak mechanical properties, related to loading direction was $0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, $75^{\circ}$ and $90^{\circ}$ while the angle of second layer, with high mechanical properties, related to loading direction was $90^{\circ}$, $105^{\circ}$, $120^{\circ}$, $135^{\circ}$, $150^{\circ}$, $160^{\circ}$ and $180^{\circ}$. Is to be note that the angle between bedding layer was $90^{\circ}$ in all bedding configurations. Also, three different pairs of the thickness were chosen in models, i.e., 5 mm/10 mm, 10 mm/10 mm and 20 mm/10 mm. The result shows that in all configurations, shear cracks develop between the weaker bedding layers. Shear cracks angel related to normal load change from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Numbers of shear cracks are constant by increasing the bedding thickness. It's to be noted that in some configuration, tensile cracks develop through the intact area of material model. There is not any failure in direction of bedding plane interface with higher strength.

• Wind and Structures
• /
• v.29 no.4
• /
• pp.247-270
• /
• 2019

Wind-resistant design of existing cooling tower structures overlooks the impacts of rainfall. However, rainstorm will influence aerodynamic force on the tower surface directly. Under this circumstance, the structural response of the super-large cooling tower (SLCT) will become more complicated, and then the stability and safety of SLCT will receive significant impact. In this paper, surrounding wind fields of the world highest (210 m) cooling tower in Northwest China underthree typical wind velocities were simulated based on the wind-rain two-way coupling algorithm. Next, wind-rain coupling synchronous iteration calculations were conducted under 9 different wind speed-rainfall intensity combinations by adding the discrete phase model (DPM). On this basis, the influencing laws of different wind speed-rainfall intensity combinations on wind-driving rain, adhesive force of rain drops and rain pressure coefficients were discussed. The acting mechanisms of speed line, turbulence energy strength as well as running speed and trajectory of rain drops on structural surface in the wind-rain coupling field were disclosed. Moreover, the fitting formula of wind-rain coupling equivalent pressure coefficient of the cooling tower was proposed. A systematic contrast analysis on its 3D distribution pattern was carried out. Finally, coupling model of SLCT under different working conditions was constructed by combining the finite element method. Structural response, buckling stability and local stability of SLCT under different wind velocities and wind speed-rainfall intensity combinations were compared and analyzed. Major research conclusions can provide references to determine loads of similar SLCT accurately under extremely complicated working conditions.

• Journal of the Korean Society of Mineral and Energy Resources Engineers
• /
• v.55 no.6
• /
• pp.614-623
• /
• 2018

The recent occurrence of consecutive large earthquakes in the southeastern part of the Korean peninsula has brought significant attention to the prevention of earthquake damage in Korea. This article aims to explore a technology-based approach for earthquake drills using state-of-the-art visual intelligence and virtual reality technologies. The technical process consists of several stages, including acquisition of image information in living spaces using a camera, recognition of objects from the acquired image information, extraction of three dimensional geometric information, simulation of virtual earthquakes using dynamic modelling techniques such as the discrete element method, and realization of the simulated earthquake in a virtual reality environment. This article provides a comprehensive analysis of the individual processes at each stage of the technical process, a survey on the current status of related technologies, and discussion of the technical challenges in its execution.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.1
• /
• pp.15-23
• /
• 2008

In this study, modeling shear characteristics of a coarse material mainly containing crushed stones were implemented using PFC2D, a commercially available code based on DEM(Discrete Element Method). Using the DEM code, this study provides the methodology considering the shear characteristics due to a irregular grain shape, GSD(Grain Size Distribution) and porosity of coarse material which are not effectively incorporated in conventional continuum numerical codes. Direct shear test was simulated for the GSD and porosity generated sample using the code and the simulated results showed very good agreement with the laboratory test results. The current modeling approach can be applied to other coarse materials having various GSD and porosities. Using such application, prediction of the strength characteristics of coarse material in field scale would be possible, which is limited in laboratory scale so far.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Journal of the Korean Geosynthetics Society
• /
• v.9 no.1
• /
• pp.47-57
• /
• 2010

Granular soils having a large particle size have been used as a filling material in the construction of foundation, harbor, dam, and so on. Consequently, the shear behavior of this granular soil plays a key role in respect of stability of structures. For example, soil particle crushing occurring at the interface between structure and soil and/or within soil mass can cause a disturbance of ground characteristics and consequently induce issues in respect of stability of structures. In order to investigate the shear behavior according to an existence and nonexistence of particle crushing, numerical analyses were conducted by using the DEM (Discrete Element Method)-based software program PFC2D (Particle Flow Code). By dividing soil particle bonding model into crushing model and noncrushing model, total four particle bonding models were simulated and their results were compared. Noncrushing model included one ball model and clump model, and crushing model included cluster model and Lobo-crushing model. The combinations of soil particle followed the research results of Lobo-Guerrero and Vallejo (2005) which were composed of eight circles. The results showed that the friction angle was in order of clump model > cluster model > one ball model. The particle bonding model compared to one ball model and noncrushing model compared to crushing model showed higher shear strength. It was also concluded that the model suggested by Lobo-Guerrero and Vallejo (2005) is not appropriate to simulate the soil particle crushing.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.7
• /
• pp.265-274
• /
• 2018

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has prompted studies directed toward the development of an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have been obtained traditionally using standard finite difference or finite element methods. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under the conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method (DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as the excessive use of storage due to the conditions of complex geometry and loading. The in-plane buckling of curved beams considering the extensibility of the arch axis was analyzed under uniformly distributed radial loads using the DQM. The critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with the precise results by other methods for cases, in which they were available. The DQM, using only a limited number of grid points, provided results that agreed very well (less than 0.3%) with the exact ones. New results according to diverse variations were obtained, showing the important roles in the buckling behavior of curved beams, and can be used in comparisons with other numerical solutions or with experimental test data.