• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.314-323
• /
• 2010

Recently, granular soils having a large particle size are frequently used as a filling material in the construction of foundation, harbor, dam, and so on. The shear behavior of this granular soil plays a key role in the stability of structures. For example, soil particle crushing occurring at the interface between structure and soil and/or within soil mass can cause the disturbance of ground characteristics and consequently induce an 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 PFC(Particle Flow Code). Using the crushing model and non-crushing model which were created in this study, numerical analyses of ring shear test were conducted and their results were analyzed and compared. In general, landslide and slope stability are accompanied by a large displacement and consequently not only a peak strength but also a residual strength are very important in the analysis of landslide and slope stability. However the direct shear test which has been commonly used in the determination of shear strength parameters has a limitation on displacement therefore the residual strength parameters can not be obtained. The characteristics of residual shear behavior were investigated through the numerical analyses in this study.

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

• 한국공간정보시스템학회:학술대회논문집
• /
• 2004.05a
• /
• pp.75-82
• /
• 2004

Many researcher's efforts have made a significant advancement of space frame structure with various portion, and it becomes the most outsanding one of space structures. However, with the characteristics of thin and long term of spacing, the unstable behavior of space structure is shown by initial imperfection, erection procedure or joint, especially space frame structure represents more. This kind of unstable problem could not be set up clearly and there is a huge difference between theory and experiment. Moreover, the discrete structure such as space frame has more complex solution, this it is not easy to derive the formulation of design about space structure. In this space frame structure, the character of rise-span ratio or load mode is represented by the instability of space frame structure with initial imperfection, and snap-through or bifurcation might be the main phenomenon. Therefore, in this study, space frame structure which has a lot of aesthetic effect and profitable for large space covering single layer is dealt. And because that the unstable behavior due to variation of inner force resistance in the elastic range is very important collapse mechanism, I would like to investigate unstable character as a nonlinear behavior with a geometric nonlinear. In order to study the instability. I derive tangent stiffness matrix using finite element method and with displacement incremental method perform nonlinear analysis of unit space structure, star dome and 3-ring star dome considering rise-span $ratio(\mu}$ and load $ratio(R_L)$ for analyzing unstable phenomenon.

• Explosives and Blasting
• /
• v.23 no.1
• /
• pp.47-54
• /
• 2005

An explosion modeling technique was developed by using the spherical discrete element code, $PFC^{3D}$, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a $PFC^{3D}$ particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). A test blast was conducted for a RC column, whose dimension was $600\times300\times1800$ in millimeters. The initial velocities of the surface movements were measured to be in the range of $14\~18\;m/s$ with the initiation times of $1.5\~2.0m$. Then the blasting procedure was simulated by using the modeling technique. The particle assembly representing the concrete was made of cement mortar and coarse aggregates, whose mirco-properties were obtained from the calibration processes. As a result, the modeling technique developed in this study made it possible for the burden to move with the velocity of $17\~24\;m/s$, which are slightly higher values compared to those of the test blast.

• Tunnel and Underground Space
• /
• v.8 no.3
• /
• pp.234-248
• /
• 1998

This paper presents a new constitutive model for numerical modeling the shear behaviour of rough rock joints. The model incorporates the dilatancy of joints on the basis of elasto-plastic theory. Barton's empirical shear strength formular are adopted in the formulation process. The mobilized JRC concept is evoked to address the shear strength hardening and sofrening phenomena. The mobilized JRC in the pre- and post-peak range is approximated by assuming that the variation of JRC is a function of tangential plastic work. Discrete finite joint element is used to implement the proposed constitutive model. The model is validated by the numerical direct shear test on a single joint which is subjected to different boundary conditions. The test results are in good agreement with the experimental observations reported by other authors. The numerical tests also exhibit that the proposed model can simulate the salient features envisaged in the behaviour of rough rock joints.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.4
• /
• pp.465-474
• /
• 2012

In this study, a dynamic simulation model that considers many spherical particles and multibody dynamics (MBD) entities is developed. Plenteous spherical particles are solved using the Discrete Element Method (DEM) technique and simulated on a GPU board in a PC. A fast algorithm is used to calculate the Hertzian contact forces between many spherical particles, and NVIDIA CUDA is used to increase the calculation speed. The explicit integration method is applied to solve the many spheres. MBD entities are simulated by recursive formulation. Constraints are reduced by recursive formulation, and the implicit generalized alpha method is applied to solve the dynamic model. A new algorithm is developed to simulate the DEM and MBD models simultaneously. As a numerical example, a truck car model and gear model are developed. The results show that the proposed algorithm using a general-purpose GPU in a PC has many advantages.

• Explosives and Blasting
• /
• v.22 no.4
• /
• pp.7-15
• /
• 2004

An explosion modeling technique was developed by using the spherical discrete element code, PFC3D, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a PFC3D particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). According to this concept, the explosion pressure is applied to the wall particles by the scheme of radius expansion/contraction of inner-hole particles. The output wall force is compared to the input hole pressure in every time step, and a correction routine is activated to control the radius multiplier of the inner-hole particles. A comparative blast simulation far a cement mortar block of $80\times90\times80mm$ was conducted by using the conventional explosion modeling method and the new one. The results of the simulation are presented in a qualitative fashion.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.1
• /
• pp.113-124
• /
• 2008

Distributed operation of overall structural design process, by which product optimization and process parallelization are simultaneously implemented, is presented in this paper. The database-interacted hybrid method, which selectively takes the accustomed procedure of the conventional method in the framework of the optimal design, is utilized here. The staged application of design constraints reduces the computational burden for large complex optimization problems. Two kinds of numeric and graphic processes are simultaneously implemented by concurrent engineering approach in the distributed environment of PC networks. The former is based on finite element optimization method and the latter is represented by AutoCAD using AutoLISP programming language. Numerical computation and database interaction on servers and graphic works on independent clients are communicated through message passing. The numerical experiments for some steel truss models show the validity and usability of the method. This study has sufficient adaptability and expandability, in that it is based on general methodologies and industry standard platforms.

• Journal of the Korean Vacuum Society
• /
• v.19 no.2
• /
• pp.155-160
• /
• 2010

The effect of process temperature of a final annealing step in the fabrication of phase change memory (PCM) devices was investigated. Discrete PCM devices employing $Ge_2Sb_2Te_5$ (GST) films as an active element were made in a pore-style configuration, and they were annealed at various temperatures ranging from 160 to $300^{\circ}C$. The behaviors of cell resistance change from SET resistance to RESET resistance were totally different according to the annealing temperatures. There was a critical annealing temperature for the fabrication of normal PCM devices and abnormal operations were observed in some devices annealed at temperatures lower or higher than the critical temperature. Those influences of annealing temperature seem closely related to the thermal stability of a top electrode/GST/heating layer multilayer structure in the PCM devices.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.12
• /
• pp.59-69
• /
• 2015

This study examined the magnitude and distribution of earth pressure on the support system in a jointed rock mass due to the different joint sets as well as varying the rock type and joint condition (joint shear strength and joint inclination angle). Based on a physical model test and its numerical simulation, a series of numerical parametric analyses were conducted using a discrete element method. The results showed that the induced earth pressure was affected significantly by a joint set depending on the inclusion of the joint inclination angle, which induces a joint sliding condition, but the number of joint sets alone was not important, even though the earth pressure could be increased slightly as the number of joint sets is increased. In addition, the study results were compared with Peck's earth pressure for soil ground, which indicated that the earth pressure in a jointed rock mass could be considerably different from that in soil ground. The study suggests that the effects of joint sets as well as rock type and joint condition are important factors affecting the earth pressure in a jointed rock mass and they should be considered when designing a support system in a jointed rock mass.