Search results

The purpose of this paper is to propose a novel combined finite-discrete element method (FDEM), based on the cohesive zone model, for simulating rockslide problems at the laboratory scale.

The combined FDEM is realized using ABAQUS/Explicit. The rock mass is represented as a collection of elastic bulk elements glued by cohesive elements with zero thickness. To reproduce the tensile and shear micro-fractures in rock material, the Mohr-Coulomb model with tension cut-off is employed as the damage initiation criterion of cohesive elements. Three simulated laboratory tests are considered to verify the capability of combined FDEM in reproducing the mechanical behavior of rock masses. Three slope models with different joint inclinations are taken to illustrate the application of the combined FDEM to rockslide simulation.

The results show that the joint inclination is an important factor for inducing the progressive failure behavior. With a low joint inclination, the slope failure process is observed to be a collapse mode. As the joint inclination becomes higher, the failure mode changes to sliding and the steady time of rock blocks is shortened. Moreover, the runout distance and post-failure slope angle decrease as the joint inclination increases.

These studies indicate that the combined FDEM performed within ABAQUS can simulate slope stability problems for research purposes and is useful for studying the slope failure mechanism comprehensively.

This paper aims to investigate the responses of laminated glass under soft body impact, including elastic impact and fracture/fragmentation consideration.

The simulation uses the combined finite-discrete element method (FDEM) which combines finite element mesh into discrete elements, enabling the accurate prediction of contact force and deformation. Material rupture is modelled with a cohesive fracture criterion, evaluating the process from continua to discontinua.

Responses of laminated glass under soft impact (both elastic and fracture) agree well with known data. Crack initiation time in laminated glass increases with the increase of the outside glass thickness. With the increase of Eprojectile, failure mode is changing from flexural to shear, and damage tends to propagate longitudinally when the contact surface increases. Results show that the FDEM is capable of modelling soft impact behaviour of laminated glass successfully.

The work is done in 2D, and it will not represent fully the 3D mechanisms.

Elastic and fracture behaviour of laminated glass under soft impact is simulated using the 2D FDEM. Limited work has been done on soft impact of laminated glass with FDEM, and special research endeavours are warranted. Benchmark examples and discussions are provided for future research.

Continua and discontinua coexist in natural rock materials. This paper aims to present an improved approach for addressing the mechanical response of rock masses based on the combined finite-discrete element method (FDEM) proposed by Munjiza.

Several algorithms have been programmed in the new approach. The algorithms include (1) a simpler and more efficient algorithm to calculate the contact force; (2) An algorithm for tangential contact force closer to the actual physical process; (3) a plastic yielding criterion (e.g. Mohr-Coulomb) to modify the elastic stress for fitting the mechanical behavior of elastoplastic materials; and (4) a complete code for the mechanical calculation to be implemented in Matrix Laboratory (MATLAB).

Three case studies, including two standard laboratory experiments (uniaxial compression and Brazilian split test) and one engineering-scale anti-dip slop model, are presented to illustrate the feasibility of the Y-Mat code and its ability to deal with multi-scale rock mechanics problems. The results, including the progressive failure process, failure mode and trajectory of each case, are acceptable compared to other corresponding studies. It is shown that, the code is capable of modeling geotechnical and geological engineering problems.

This article gives an improved FDEM-based numerical calculation code. And, feasibility of the code is verified through three cases. It can effectively solve the geotechnical and geological engineering problems.

This paper aims to present a new numerical model for geometric nonlinear analysis of thin-shell structures based on a combined finite-discrete element method (FDEM).

The model uses rotation-free, three-node triangular finite elements with exact formulation for large rotations, large displacements in conjunction with small strains.

The presented numerical results related to behaviour of arbitrary shaped thin shell structures under large rotations and large displacement are in a good agreement with reference solutions.

This paper presents new computationally efficient numerical model for geometric nonlinear analysis and prediction of the behaviour of thin-shell structures based on combined FDEM. The model is implemented into the open source FDEM package “Yfdem”, and is tested on simple benchmark problems.

In order to avoid the problem of volumetric locking often encountered when using constant strain tetrahedral finite elements, the purpose of this paper is to present a new composite tetrahedron element which is especially designed for the combined finite-discrete element method (FDEM).

A ten-noded composite tetrahedral (COMPTet) finite element, composed of eight four-noded low order tetrahedrons, has been implemented based on Munjiza’s multiplicative decomposition approach. This approach naturally decomposes deformation into translation, rotation, plastic stretches, elastic stretches, volumetric stretches, shear stretches, etc. The problem of volumetric locking is avoided via a selective integration approach that allows for different constitutive components to be evaluated at different integration points.

A number of validation cases considering different loading and boundary conditions and different materials for the proposed element are presented. A practical application of the use of the COMPTet finite element is presented by quantitative comparison of numerical model results against simple theoretical estimates and results from acrylic fracturing experiments. All of these examples clearly show the capability of the composite element in eliminating volumetric locking.

For this tetrahedral element, the combination of “composite” and “low order sub-element” properties are good choices for FDEM dynamic fracture propagation simulations: in order to eliminate the volumetric locking, only the information from the sub-elements of the composite element are needed which is especially convenient for cases where re-meshing is necessary, and the low order sub-elements will enable robust contact interaction algorithms, which maintains both relatively high computational efficiency and accuracy.

The purpose of this paper is to discretely model rockfill materials considering the irregular shape of the particles and their crushability. The scientific goal was to investigate the influence of particle crushability and shape on the mechanical behavior of rockfill materials.

The method of generating irregular-shaped particles was based on the observation that most rockfill grains can be approximately circumscribed by an ellipsoid. Two shape descriptors were used to make the virtual particles closely replicate the geometric features of natural rockfill grains. The combined finite-discrete element method (FDEM) was used to numerically simulate a drained, tri-axial compression test. The particle assemblies were subjected to tri-axial compression under strain controlled conditions while a constant confining pressure was maintained.

The non-breakable particles showed a remarkable ability to dilate as a result of a higher inter-particle locking effect. Dilation forces the particles to move from a lower potential energy state to a higher potential energy state, which causes the micro-structure to become less stable, resulting in a dramatic decline in the angle of friction from the peak state to the residual state. In addition, the elongated particles enhance the interlocking effect, but breakage is also more likely to occur. The net effect of those two mechanisms controls the overall shearing resistance of rockfill materials.

After calibration using a few micro-parameters, the combined FDEM was able to reproduce the typical behavior of rockfill materials without requiring a description of the complex relationship that exists between constituents; this relationship must be described in continuum mechanics. The simulation results showed that this approach is predictive. The combined FDEM also provides an opportunity for a quantitative study of the micro-structure of granular materials, and this study will help us to better understand the mechanical characteristics of rockfill materials.

Contact interaction and contact detection (CD) remain key components of any discontinua simulations. The methods of discontinua include combined finite-discrete element method (FDEM), discrete element method, molecular dynamics, etc. In recent years, a number of CD algorithms have been developed, such as Munjiza–Rougier (MR), Munjiza–Rougier–Schiava (MR-S), Munjiza-No Binary Search (NBS), Balanced Binary Tree Schiava (BBTS), 3D Discontinuous Deformation Analysis and many others. This work aims to conduct a numerical comparison of certain algorithms often used in FDEM for bodies of the same size. These include MR, MR-S, NBS and BBTS algorithms.

Computational simulations were used in this work.

In discrete element simulations where particles are introduced randomly or in which the relative position between particles is constantly changing, the MR and MR-S algorithms present an advantage in terms of CD times.

This paper presents a detailed comparison between CD algorithms. The comparisons are performed for problem cases with different lattices and distributions of particles in discrete element simulations. The comparison includes algorithms that have not been evaluated between them. Also, two new algorithms are presented in the paper, MR-S and BBTS.

The purpose of this paper is to investigate the evolution of the permeability of spherical packing during cold compaction by pore-scale modeling.

The discrete element method (DEM) is used to generate spherical packing structure under different compressive pressures and the Lattice Boltzmann method (LBM) is adopted to calculate the permeability of each spherical assembly.

It is found that the decrease of the porosity is the main reason of the reduction in permeability in the initial compression stage, but its influence becomes insufficient in the late compression stages. Besides, two empirical formulas are obtained, which describe the relation between the permeability and the equivalent mean diameter and the variation of normalized permeability with compressive pressure, respectively.

In this study, the authors study the spherical particles and ignore the non-spherical effects. Besides, the classical contact model, the linear-spring-damping model, is used in DEM, so the plastic deformation cannot be considered.

The DEM and the LBM are well combined to study the compaction effects on permeability of spherical packing. Two simple expressions of the spherical packing structure with uniform diameter distribution are given for the first time.

The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect.

In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method.

The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised.

A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field.

Accurate presentation of the rock microstructure is critical to the grain-scale analysis of rock deformation and failure in numerical modelling. 3D granite microstructure modelling has only been used in limited studies with the mineral pattern often remaining poorly constructed. In this study, the authors developed a new approach for generating 2D and 3D granite microstructure models from a 2D image by combining a heterogeneous material reconstruction method (simulated annealing method) with Voronoi tessellation.

More specifically, the stochastic information in the 2D image is first extracted using the two-point correlation function (TPCF). Then an initial 2D or 3D Voronoi diagram with a random distribution of the minerals is generated and optimised using a simulated annealing method until the corresponding TPCF is consistent with that in the 2D image. The generated microstructure model accurately inherits the stochastic information (e.g. volume fraction and mineral pattern) from the 2D image. Lastly, the authors compared the topological characteristics and mechanical properties of the 2D and 3D reconstructed microstructure models with the model obtained by direct mapping from the 2D image of a real rock sample.

The good agreements between the mapped and reconstructed models indicate the accuracy of the reconstructed microstructure models on topological characteristics and mechanical properties.

The newly developed reconstruction method successfully transfers the mineral pattern from a granite sample into the 2D and 3D Voronoi-based microstructure models ready for use in grain-scale modelling.