• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.66-81
• /
• 2021

The effects of directional fracture tensor components and first invariant of fracture tensor on deformation moduli and shear moduli of fractured rock masses is analyzed based on regression analysis performed between 3-D fracture tensor parameters and deformability of DFN blocks. Using one or two deterministic joint sets, a total of 224 3-D discrete fracture network (DFN) cube blocks were generated with various configurations of deterministic density and probabilistic size distribution. The fracture tensor parameters were calculated for each generated DFN systems. Also, deformability moduli with respect to three perpendicular direction of the DFN cube blocks were estimated based on distinct element method. The larger the first invariant of fracture tensor, the smaller the values for the deformability moduli of the DFN blocks. These deformability properties present an asymptotic pattern above the certain threshold. It is found that power-law function describes the relationship between the directional deformability moduli and the corresponding fracture tensor components estimated in same direction.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.4
• /
• pp.123-131
• /
• 2006

The majority of pedestrian fatalities and injuries are caused by vehicle-pedestrian accidents. Recently, it has been recognized as a serious problem. Injuries of occupants in a vehicle have been decreased considerably. However, efforts for protection of pedestrians are still insufficient. These days, many advanced industries are striving for a better protection of pedestrians by using an active hood lift system, rather than reforming the existing structure. In this research, the active hood lift system is designed to enhance the performance for protection. The active hood lift system is analyzed by using the nonlinear finite element method. An optimization problem is formulated by incorporation of the analysis results. Orthogonal arrays are utilized to solve the formulated problem. An iterative optimization algothrithm using orthogonal arrays is utilized for design in the discrete space. It is found that the method can remarkably decrease the number of function evaluations.

• Tunnel and Underground Space
• /
• v.15 no.4 s.57
• /
• pp.250-263
• /
• 2005

Progressive and localized brittle failures around an excavated opening by the overstressed condition can act as a serious obstacle to ensure the stability and the economical efficiency of construction work. In this paper, the characteristics of the brittle failure around an circular opening with stress level was studied by the biaxial compressive test using sealed specimen and by the numerical simulation with $PFC^{2D}$, one of the discrete element codes. The occurring pattern and shape of the brittle failure around a circular opening monitored during the biaxial loading were well coincided with those of the stress induced failures around the excavated openings observed in the brittle rock masses. The crack development stages with stress level were evaluated by the detailed analysis on the acoustic emission event properties. The microcrack development process around a circular opening was successfully visualized by the particle flow analysis. It indicated that the scaled test had a good feasibility in understanding the mechanism of the brittle failure around an opening with a high reliability.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.1
• /
• pp.17-23
• /
• 2014

Graphics processing units (GPUs) are ideal for solving problems involving parallel data computations. In this study, the GPU is used for effectively carrying out a multi-body dynamic simulation with particle dynamics. The Hilber-Hushes-Taylor (HHT) implicit integration algorithm is used to solve the integral equations. For detecting collisions among particles, the spatial subdivision algorithm and discrete-element methods (DEM) are employed. The developed program is verified by comparing its results with those of ADAMS. The numerical efficiencies of the serial program using the CPU and the parallel program using the GPU are compared in terms of the number of particles, and it is observed that when the number of particles is greater, more computing time is saved by using the GPU. In the present example, when the number of particles is 1,300, the computational speed of the parallel analysis program is about 5 times faster than that of the serial analysis program.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.5
• /
• pp.169-177
• /
• 2011

In this study, an impact resistance of concrete segmented composites adopted shell's structures which have the excellent impact resistance was assessed. In order to enhance the performances of concrete segmented composite, the bond strength of mortar between the concrete blocks should be improved. Hence, in this study polymer mortars were applied to increase the bond strength of mortar. From the results of bond tests, the 15% latex mortar was selected and static and low-velocity impact tests were carried out for the specimens applied the plain and latex mortar. The concrete segmented composites, of which the bond strength of mortar was enhanced, showed improved low-velocity impact resistances. A Nonlinear finite element analysis using the discrete crack model showed similar energy dissipating capacities to the impact test's results. Consequently, by improving the analysis models for segmented composites, the impact resistances for manifold variables can be predicted and assessed.

• Proceedings of the KSR Conference
• /
• 2002.10a
• /
• pp.648-653
• /
• 2002

2차원 개별요소법을 이용하여 철도사면에 있어서 낙석이 발생되었을 경우 낙석에 의한 운동에너지, 위치에너지 및 마찰에 의한 에너지변화를 사면과 낙석사이의 마찰계수를 매개변수로 하여 분석하였다. 낙석의 낙하양태(미끄러짐과 회전)에 따라 변화하는 1차지면 충돌시간, 2차지면 충돌시간, 2차충돌까지의 비산거리, 1차지면 충돌 후 도약높이, 사면 운동 동안의 낙석 회전횟수 등을 구하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.1
• /
• pp.143-149
• /
• 2016

This study evaluated the shear behavior of geogrid reinforced ballast material using a large scale direct shear test and discrete element method (DEM) based on PFC 3D program. The direct shear test was conducted on ballast materials that have different particle size distributions. Whereas the test results revealed that the shear strength generally increased with the larger particle size of ballast material without geogrid reinforcement, the shear behavior of ballast material was found to change pertaining to the relationship between particle size distribution and geogrid rib size. Generally, it is deemed the effectiveness of reinforcement can be achieved when the rib size is two times greater than average particle size. A numerical analysis based on DEM was conducted to verify the test results. The geogrid modeling was successfully completed by calibration process along with sensitivity analysis to have actual tensile strength provided by manufacturer. With a given geogrid model, the parametric evaluation was further carried out to examine the interactive behavior between geogrid and ballast material. Consequently, it was found that the effectiveness zone of geogrid reinforcement generated within a specific depth.

• Explosives and Blasting
• /
• v.41 no.3
• /
• pp.62-72
• /
• 2023

The aging of plant structures due to industrialization in the 1970s has increased the demand for blast demolition. While blasting can reduce exposure to environmental pollution by shortening the demolition period, improper blasting design and construction plans pose significant safety risks. Thus, it is vital to consider optimal blasting demolition conditions and other factors through collapse behavior simulation. This study utilizes a 3-D combined finite-discrete element method (FDEM) code-based 3-D DFPA to simulate the collapse of a chimney structure in a thermal power plant in Seocheon, South Korea. The collapse behavior from the numerical simulation is compared to the actual structure collapse, and the numerical simulation result presents good agreement with the actual building demolition. Additionally, various numerical simulations have been conducted on the chimney models to analyze the impact of the duct size in the pre-weakening area. The no-duct, duct, and double-area duct models were compared in terms of crack pattern and history of Z-axis displacement. The findings show that the elapse-time for demolition decreases as the area of the duct increases, causing collapse to occur quickly by increasing the load-bearing area.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.6
• /
• pp.12-19
• /
• 2003

Corrugated bulkheads for a bulk carrier are divided into watertight bulkheads and deep tank bulkheads. Design of the watertight bulkheads is principally determined by the permissible limit of Classification and IACS requirements. But, the verification of strength through finite element analysis is indispensable for design of the deep tank bulkheads. A stage for stress evaluation of corrugated part is required for optimum structural design of the deep tank bulkheads. Since the finite element analysis for real model requires excessive amount of calculation time, in this study one corrugated structure is replaced with beam element and is idealized as 2 dimensional frame structure connected to upper and lower stool Minimum weight design of the deep tank bulkheads is performed through generalized sloped deflection method(GSDM) as direct calculation method. The purpose of this study is the development of design system for the minimization of steel weight of deep tank bulkheads as well as watertight bulkheads. Discrete variables are used as design variables for the practical design. Evolution strategies(ES) is used as an optimization technique.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.