Search results

YADE‐OPEN DEM is an open‐source software based on the discrete element method, (DEM) which uses object oriented programming techniques. The purpose of this paper is to describe the software architecture.

The DEM chosen uses position, orientation, velocity and angular velocity as independent variables of simulated particles which are subject to explicit leapfrog time‐integration scheme (Lagrangian method). The three‐dimensional dynamics equations based on the classical Newtonian approach for the second law of motion are used. The track of forces and moments acting on each particle is kept at every time step. Contact forces depend on the particle geometry overlap and material properties. The normal, tangential and moment components of interaction force are included.

An effort is undertaken to extract the underlying object oriented abstractions in the DEM. These abstractions are implemented in C⁺⁺, conform to object oriented design principles and use design patterns. Based on that, a software framework is developed in which the abstractions provide the interface where the modelling methods can be plugged‐in.

The resulting YADE‐OPEN DEM framework is designed in a generic way which provides great flexibility when adding new scientific simulation code. Some of the advantages are that numerous simulation methods can be coupled within the same framework while plug‐ins can import data from other software. In addition, this promotes code improvement through open‐source development and allows feedback from the community. However implementing such models requires that one adheres to the framework design and the YADE framework is a new emerging software. To download the software see http : //yade.wikia.com webpage.

The purpose of this paper is to describe how a discrete element model is used to predict the penetration depth and the perforation caused by a non‐deformable missile against a thin reinforced concrete slab.

Initial calibration of the model was done with a series of flat‐nose missile tests. Additional simulations were performed with varying the percentage of reinforcement. The present numerical model is compared to experimental test data provided by the French Atomic Energy Agency (CEA) and the French Electrical Power Company (EDF).

For thin concrete slabs, the evolution of the penetration depth in terms of percentage of reinforcement was compared with experimental results: quantitatively the results are very coherent.

The modeling scale is higher than the heterogeneity scale, so the model may be used to simulate real structures, which means that the discrete element method is mainly used here for its ability to account for discontinuities; an identification process based on quasi‐static tests is used, so the quasi‐static behavior of concrete is reproduced. This identification process is the key point, to allow a complete predictive computation for complex impact configurations, especially when the missile diameter and the thickness of the concrete slab are on the same order in size.

The purpose of this paper is to describe a novel implementation of a multispatial method, multitime-scheme subdomain differential algebraic equation (DAE) framework allowing a mix of different space discretization methods and different time schemes by a robust generalized single step single solve (GS4) family of linear multistep (LMS) algorithms on a single body analysis for the first-order nonlinear transient systems.

This proposed method allows the coupling of different numerical methods, such as the finite element method and particle methods, and different implicit and/or explicit algorithms in each subdomain into a single analysis with the GS4 framework. The DAE, which constrains both space and time in multi-subdomain analysis, combined with the GS4 framework ensures the second-order time accuracy in all primary variables and Lagrange multiplier. With the appropriate GS4 parameters, the algorithmic temperature rate variable shift can be matched for all time steps using the DAE. The proposed method is used to solve various combinations of spatial methods and time schemes between subdomains in a single analysis of nonlinear first-order system problems.

The proposed method is capable of coupling different spatial methods for multiple subdomains and different implicit/explicit time integration schemes in the GS4 framework while sustaining second-order time accuracy.

Traditional approaches do not permit such robust and flexible coupling features. The proposed framework encompasses most of the LMS methods that are second-order time accurate and unconditionally stable.

The purpose of this study is the reduction of computational time demand of discrete element based modeling.

The methodology is the systematic changing of particle size and micromechanical parameters to reduce computational time requirements.

In some cases, the computational demand of discrete simulations can be reduced to about 95 per cent.

Based on the results and demonstrated methodology, the enormous computational time demand of discrete element-based modeling can be reduced significantly.

The purpose of this paper is to determine the representative volume element (RVE) size for quasi‐brittle materials using a discrete approach, namely a lattice‐particle model.

Different material samples are generated and subjected to study regarding its size, maximum aggregate size and boundary conditions. In order to determine the mechanical properties such as the elastic modulus, Poisson's ratio or tensile strength, several tension tests are carried out. For this purpose, a lattice‐particle approach is used to model concrete's fracturing behavior. The information provided by the previous simulations is implemented in a statistical analysis to determine the size of the RVE.

The determination of the RVE size for quasi‐brittle materials is successfully achieved by means of a lattice‐particle model. Computed results show a good agreement with other results reported in the bibliography.

Within a general multiscale framework, the determination of the RVE size is of great interest and some studies have been performed for random heterogeneous materials. However, these analyses are mainly continuum‐based. The estimation of the RVE size is important for correctly predicting the mechanical properties and can be used in different multiscale schemes.

To study fracture characteristics of jointed rock masses under blasting load, the RFPA2D analysis software for dynamic fracture of rocks based on the finite element method and statistical damage theory was used.

On this basis, this research simulated the fracture process of rock masses in blasting with different joint geometrical characteristics and mainly analysed the influences of distance from joints to blasting holes, the length of joints, the number of joints and joint angle on fracture of rock masses.

The calculation results show that with the constant increase of the distance from joints to blasting holes, the influences of joints on blasting effects of rock masses gradually reduced. Rock masses with long joints experienced more serious damages than those with short joints. Damages obviously increased with the changing from rock masses without joints to rock masses with joints, and when there were three joints, the further increase of the number of joints had unobvious changes on blasting effects of rock masses. Joints showed significant guidance effect on the propagation of cracks in blasting: promoting propagation of main vertical cracks deflecting to the ends of joints.

The research results are expected to provide some theoretical bases in practical application of engineering blasting.

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 simulate a punch shear test of partly consolidated ice ridge keel by using a three-dimensional discrete element method. The authors model the contact forces between discrete ice blocks with Hertz–Mindlin contact model. For freeze bonds between the ice blocks, the authors apply classical linear cohesion model with few modifications. Based on punch shear test simulations, the authors are able to determine the main characteristics of an ice ridge from the material parameters of the ice and freeze bonds.

The authors introduced a discrete model for ice that can be used for modelling of ice ridges. The authors started with short introduction to current status with ice ridge modelling. Then they introduced the model, which comprises Hertz–Mindlin contact model and freeze bond model with linear cohesion and softening. Finally, the authors presented the numerical results obtained using EDEM is commercial Discrete Element Modeling software (EDEM) and analysed the results.

The Hertz–Mindlin model with cohesive freeze bonds and linear softening is a reasonable model for ice rubble. It is trivial that the ice blocks within the ice ridge are not spherical particles, but according to results, the representation of ice blocks as spheres gave promising results. The simulation results provide information on how the properties of freeze bond affect the results of punch shear test. Thus, the simulation results can be used to approximate the freeze bonds properties within an ice ridge when experimental data are available.

As the exact properties of ice rubble are unknown, more research is required both in experimental and theoretical fields of ice rubble mechanics.

Based on this numerical study, the authors are able to determine the main characteristics of an ice ridge from material parameters of ice and freeze bonds. Furthermore, the authors conclude that the model creates a promising basis for further development in other applications within ice mechanics.

The purpose of this paper is to use the discrete element method (DEM) to model the fracture behaviour of brittle materials in 2D.

The material consists of a set of particles in contact with a close‐packed structure. It allows the derivation of an expression for the stress intensity factor as a function of the contact forces near the crack tip. A classical failure criterion, based on the material's toughness, is then adopted for the analysis of crack propagation, represented by the contact loss between particles.

The DEM approach is compared to two tensile cases (mode I); both presenting a monotonous convergence towards classical solutions for more precise discretization.

The paper proposes a DEM approach in fracture mechanics of isotropic brittle materials entirely compatible with continuous classical theory. Hence the toughness value is directly introduced as a parameter of the material without any previous calibration of the DEM.

Porous concrete is a mixture of open‐graded coarse aggregate, water and cement. It is also occasionally referred to as no‐fines concrete or pervious concrete. Due to its high infiltration capacity, it is viewed as an environmentally sustainable paving material for use in urban drainage systems since it can lead to reduced flooding and to the possibilities of stormwater harvesting and reuse. However, the high porosity is due in the main part to the lack of fine aggregate particles used in the manufacture of porous concrete. The purpose of this paper is to present a numerical method to understand more fully the structural properties of porous concrete. This method will provide a useful tool for engineers to design with confidence higher strength porous concrete systems.

In the method, porous concrete is modelled using a discrete element method (DEM). The mechanical behaviour of a porous concrete sample subjected to compressive and tensile forces is estimated using two‐dimensional Particle Flow Code (PFC2D).

Three numerical examples are given to verify the model. A comprehensive set of micro‐parameters particularly suitable for porous concrete is proposed. The accuracy and effectiveness of simulation are confirmed by comparison with experimental results and empirical equations.

The experimental investigations for porous concrete described in this paper have been designed and conducted by the authors. In addition, the type of two dimensional PFC analysis presented has rarely been used to model porous concrete strength characteristics and from the results presented in this paper, this analysis technique has good potential for predicting its mechanical properties.