• Structural Engineering and Mechanics
• /
• v.80 no.5
• /
• pp.523-538
• /
• 2021

The mechanical behaviour and stability of coal mining engineering underground is significantly affected by ground water. In this study, nuclear magnetic resonance imaging (NMRI) technique was employed to determine the water distribution characteristics in coal specimens during saturation process, based on which the functional rule for water distribution was proposed. Then, using discrete element method (DEM), an innovative numerical modelling method was developed to simulate water-weakening effect on coal behaviour considering moisture content and water distribution. Three water distribution numerical models, namely surface-wetting model, core-wetting model and uniform-wetting model, were established to explore the water distribution influences. The feasibility and validity of the surface-wetting model were further demonstrated by comparing the simulation results with laboratory results. The investigation reveals that coal mechanical properties are affected by both water saturation coefficient and water distribution condition. For all water distribution models, micro-cracks always initiate and nucleate in the water-rich area and thus lead to distinct macro fracture characteristics. With the increase of water saturation coefficient, the failure of coal tends to be less violent with less cracks and ejected fragments. In addition, the core-wetting specimen is more sensitive to water than specimens with other water distribution models.

• Journal of the Korean Geosynthetics Society
• /
• v.10 no.2
• /
• pp.1-12
• /
• 2011

In order to analyze the influence of particle bonding and crushing on the characteristics of shear behavior, especially residual shear behavior of granular soil, ring shear test was simulated by using DEM(Discrete Element Method)-based software program PFC(Particle Flow Code). Total four models including two non-crushing models and two crushing models were created in this study by using clump or cluster model built in PFC. The applicability of Lobo-crushing model proposed by Lobo-Guerrero and Vallejo(2005) was investigated. In addition, the results of ring shear test were analyzed and compared with those of direct shear test. The results showed that the modelling of ring shear test should be conducted to investigate the residual shear behavior. The Lobo-crushing model cannot be applied to investigate the residual shear strength. Finally, it can be concluded that the numerical models excluding Lobo-crushing model suggested in this study can be used extensively for other studies concerning the residual shear behavior of granular soil including soil crushing.

• Journal of the Korean Geotechnical Society
• /
• v.25 no.10
• /
• pp.41-53
• /
• 2009

In order to analyze the influence of particle shape on evolution of particle crushing and characteristic of shear behavior of granular soil, direct shear test was simulated by using DEM (Discrete Element Method). Six particle shapes were generated by clump and cluster model built in PFC (Particle Flow Code). The results of direct shear test for six particle shapes were compared and analyzed with those for circular particle shape. The results of numerical tests showed a good agreement with those of experimental tests, thus the appropriateness of numerical modelling set in this study was proved. As for particle shape, more angular and rougher particle induced larger internal friction angle and more particle crushing than relatively round and smooth particle. When particles were crushed, crushing was concentrated on the shear band adjacent to the shear plane. Finally, it can be concluded that the numerical models suggested in this study can be used extensively for other studies concerning the shear behavior of granular soil including soil crushing.

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

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

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

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.306-319
• /
• 2020

Characterizations of Excavation Damage Zone (EDZ), which is hydro-mechanical degrading the host rock, are the important issues on the geological repository for the spent nuclear fuel. In the DECOVALEX 2019 project, Task G aimed to model the fractured rock numerically, describe the hydro-mechanical behavior of EDZ, and predict the change of the hydraulic factor during the lifetime of the geological repository. Task G prepared two-dimensional fractured rock model to compare the characteristics of each simulation tools in Work Package 1, validated the extended three-dimensional model using the TAS04 in-situ interference tests from Äspö Hard Rock Laboratory in Work Package 2, and applied the thermal and glacial loads to monitor the long-term hydro-mechanical response on the fractured rock in Work Package 3. Each modelling team adopted both Finite Element Method (FEM) and Discrete Element Method (DEM) to simulate the hydro-mechanical behavior of the fracture rock, and added the various approaches to describe the EDZ and fracture geometry which are appropriate to each simulation method. Therefore, this research can introduce a variety of numerical approaches and considerations to model the geological repository for the spent nuclear fuel in the crystalline fractured rock.