• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.751-758
• /
• 2005

In this study, the Discrete Element Code was applied to the analysis of falling rock blocks. The simulation was performed using the PFC2D computer code. Falling rock blocks should be applied as additional force to each others. The force affect the motion of falling rock blocks. This was used to find out the behavior of each blocks. This study revealed that the DEM can successfully capture the behavior of falling rock blocks.

• Journal of the Korean Society for Railway
• /
• v.9 no.3 s.34
• /
• pp.271-276
• /
• 2006

Insight into behaviour of pedestrians as welt as tools to assess passenger flow condition is important in such instances as planning and geometric design of railway station under regular and safety-critical circumstances. Algorithm for passenger flow analysis based on DEM (Discrete Element Method) is newly developed. There are lots of similarity between particle-laden two phase flow and passenger flow. The velocity component of 1st phase corresponds to the unit vector of calculation cell, each particle to passenger, volume fraction to population density and the particle velocity to the walking velocity, etc. And, the walking velocity of passenger is also represented by the function of population density. Key algorithms are developed to determine the position of passenger, population density and numbering to each passenger. To verify the effectiveness of new algorithm, passenger flow analysis for the basic models of railway station is conducted.

• Computers and Concrete
• /
• v.8 no.6
• /
• pp.677-692
• /
• 2011

Particle packing on meso-level has a significant influence on workability of fresh concrete and also on the mechanical and durability properties of the matured material. It was demonstrated earlier that shape exerts but a marginal influence on the elastic properties of concrete provided being packed to the same density, which is not necessarily the case with different types of aggregate. Hence, elastic properties of concrete can be treated as approximately structure-insensitive parameters. However, fracture behaviour can be expected structure-sensitive. This is supported by the present study based on discrete element method (DEM) simulated three-phase concrete, namely aggregate, matrix and interfacial transition zones (ITZs). Fracture properties are assessed with the aid of a finite element method (FEM) based on the damage materials model. Effects on tensile strength due to grain shape and packing density are investigated. Shape differences are shown to have only modest influence. Significant effects are exerted by packing density and physical-mechanical properties of the phases, whereby the ITZ takes up a major position.

• Structural Engineering and Mechanics
• /
• v.66 no.5
• /
• pp.569-582
• /
• 2018

Effects of crack separation, bridge area, on the tensile behaviour of concrete are studied experimentally and numerically through the Brazilian tensile test. The physical data obtained from the Brazilian tests are used to calibrate the two-dimensional particle flow code based on discrete element method (DEM). Then some specially designed Brazilian disc specimens containing two parallel cracks are used to perform the physical tests in the laboratory and numerically simulated to make the suitable numerical models to be tested. The experimental and numerical results of the Brazilian disc specimens are compared to conclude the validity and applicability of these models used in this research. Validation of the simulated models can be easily checked with the results of Brazilian tests performed on non-persistent cracked physical models. The Brazilian discs used in this work have a diameter of 54 mm and contain two parallel centred cracks ($90^{\circ}$ to the horizontal) loaded indirectly under the compressive line loading. The lengths of cracks are considered as; 10 mm, 20 mm, 30 mm and 40 mm, respectively. The visually observed failure process gained through numerical Brazilian tests are found to be very similar to those obtained through the experimental tests. The fracture patterns demonstrated by DEM simulations are mostly affected by the crack separation but the tensile strength of bridge area is related to the fracture pattern and failure mechanism of the testing samples. It has also been shown that when the crack lengths are less than 30 mm, the tensile cracks may initiate from the cracks tips and propagate parallel to loading direction till coalesce with the other cracks tips while when the cracks lengths are more than 30 mm, these tensile cracks may propagate through the intact concrete itself rather than that of the bridge area.

• Geomechanics and Engineering
• /
• v.23 no.5
• /
• pp.455-466
• /
• 2020

This paper presents analysis of the interaction between tunnel and Qanat with a particular interest for the optimization of Qanat shape using the discrete element code, PFC2D, and the results will be compared with the FEM results of PLAXIS2D. For these concerns, using software PFC2D based on Discrete Element Method (DEM), a model with dimension of 100m ^⁎ 100 m was prepared. A circular tunnel with dimension of 9 m was situated 20 m below the ground surface. Also one Qanat was situated perpendicularly above the tunnel roof. Distance between Qanat center and ground surface was 8 m. Five different shapes for Qanat were selected i.e., square, semi-circular, vertical ellipse, circular and horizontal ellipse. Confining pressure of 5 MPa was applied to the model. The vertical displacement of balls situated in ground surface was picked up to measure the ground subsidence. Also two measuring circles were situated at the tunnel roof and at the Qanat roof to check the vertical displacements. The properties of the alluvial soil of Tehran city are: γ_dry=19 (KN/㎥), E= 750 (kg/㎠), ν=0.35, c=0.3(kg/㎠), φ=34°. In order to validate the DEM results, a comparison between the numerical results (obtained in this study) and analytical and field monitoring have been done. The PFC2D results are compared with the FEM results. The results shows that when Qanat has rectangular shape, the tensile stress concentration at the Qanat corners has maximum value while it has minimum value for vertical ellipse shape. The ground subsidence for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. The vertical displacements at the tunnel roof for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. Historical shape of Qante approved the finding of this research.

• 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 Korean Tunnelling and Underground Space Association
• /
• v.23 no.6
• /
• pp.423-438
• /
• 2021

The blade of motor grader is used for scattering and leveling the aggregates on the foundation of road construction site. The paper performed a design improvement research of the blade part to enhance the working efficiency of motor graders. The scattering works of aggregates by blade driving were simulated by DEM (discrete element method) of a dynamic code. The four design parameters were selected and a specific leveling scenario for the simulation was determined. The nine blade models were numerically experimented, and the sensitivity of each factors was analyzed. Next, the design factors that influence a blade performance have been selected by ANOVA, and these key design factors were applied to the progressive quadratic response surface method (PQRSM). The optimum set of design factors of the blade was finally proposed.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.339-352
• /
• 2023

Prediction of the time-related traits of pebble flow inside pebble-bed HTGRs is of great significance for reactor operation and design. In this work, an image-driven approach with the aid of a convolutional neural network (CNN) is proposed to predict the remaining time of initially loaded pebbles and the time interval of paired flow images of the pebble bed. Two types of strategies are put forward: one is adding FC layers to the classic classification CNN models and using regression training, and the other is CNN-based deep expectation (DEX) by regarding the time prediction as a deep classification task followed by softmax expected value refinements. The current dataset is obtained from the discrete element method (DEM) simulations. Results show that the CNN-aided models generally make satisfactory predictions on the remaining time with the determination coefficient larger than 0.99. Among these models, the VGG19+DEX performs the best and its CumScore (proportion of test set with prediction error within 0.5s) can reach 0.939. Besides, the remaining time of additional test sets and new cases can also be well predicted, indicating good generalization ability of the model. In the task of predicting the time interval of image pairs, the VGG19+DEX model has also generated satisfactory results. Particularly, the trained model, with promising generalization ability, has demonstrated great potential in accurately and instantaneously predicting the traits of interest, without the need for additional computational intensive DEM simulations. Nevertheless, the issues of data diversity and model optimization need to be improved to achieve the full potential of the CNN-aided prediction tool.

• Geomechanics and Engineering
• /
• v.33 no.2
• /
• pp.133-140
• /
• 2023

Soil erosion can cause scouring and failures of underwater structures, therefore, various soil improvement techniques are used to increase the soil erosion resistance. The microbially induced calcium carbonate precipitation (MICP) method is proposed to increase the erosion resistance, however, there are only limited experimental and numerical studies on the use of MICP treatment for improvement of surface erosion resistance. Therefore, this study investigates the improvement in surface erosion resistance of sands by MICP through laboratory experiments and numerical modeling. The surface erosion behaviors of coarse sands with various calcium carbonate contents were first investigated via the erosion function apparatus (EFA). The test results showed that MICP treatment increased the overall erosion resistance, and the contribution of the precipitated calcium carbonate to the erosion resistance and critical shear stress was quantified in relation to the calcium carbonate contents. Further, these surface erosion processes occurring in the EFA test were simulated through the coupled computational fluid dynamics (CFD) and discrete element method (DEM) with the cohesion bonding model to reflect the mineral precipitation effect. The simulation results were compared with the experimental results, and the developed CFD-DEM model with the cohesion bonding model well predicted the critical shear stress of MICP-treated sand. This work demonstrates that the MICP treatment is effective in improving soil erosion resistance, and the coupled CFD-DEM with a bonding model is a useful and promising tool to analyze the soil erosion behavior for MICP-treated sand at a particle scale.