Search results

Simple discrete element models using PFC^2D models with bonded assemblies of particles were used to numerically simulate direct tension and block punching tests on thin spray‐on tunnel liner materials to gain insight about the liner support mechanisms. PFC^2D input parameters were calibrated such that the rupture load and elongation at rupture were similar to the laboratory test data. The calibrated model of the liner material was then used to simulate a liner around a highly stressed tunnel in rock where stresses caused extensive fracturing near the top of the tunnel. The effect of the liner was analysed by modelling the tunnel with and without the liner and showed that the liner had minimal impact on fracture propagation in the rock because of the liner's highly deformable nature. However, the liner was able to retain the fractured rock in place.

This paper aims to simulate the effect of counterface roughness on the friction transfer and wear of the polymer material sliding against steel.

The dynamic process of friction transfer and wear of polytetrafluoroethylene (PTFE) sliding against steel 45 was simulated by the software of particle flow code in two dimensions and a discrete element method. The effect of the counterface roughness was considered in the simulation. The definitions of the transferred particle and worn particle were given.

The simulation results showed that a transferred particle layer was formed on the surface of steel 45 during friction. The wear rate of PTFE can be effectively reduced by the formation of the transferred particle layer. The formation and stability of this particle layer is certainly affected by the counterface roughness (Rz). In this paper, the transferred particle numbers increased with Rz increase. And so did the worn particle numbers. However, there was little effect of Rz on the wear rate of PTFE.

The dynamic process of the friction transfer and wear of the PTFE/ steel 45 friction pair was reproduced at the micro-level. Then, the transfer and wear were quantitatively exhibited. The relations between the transfer or wear and counterface roughness was simulated and discussed. It will be meaningful for the optimization and effective control of friction and wear of polymer/metal sliding system.

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.

The interaction between rock mass structural planes and dynamic stress levels is important to determine the stability of rock mass structures in underground geotechnical engineering. In this work, the authors aim to focus on the degradation effects of fracture geometric parameters and unloading stress paths on rock mechanical properties.

A three-dimensional Particle Flow Code (PFC3D) was used for a systematic numerical simulation of the strength failure and cracking behavior of granite specimens containing prefabricated cracks under conventional triaxial compression and triaxial unilateral unloading. The authors demonstrated the unique mechanical response of prefabricated fractured rock under two conditions. The crack initiation, propagation, and coalescence process of pre-fissured specimens were analyzed in detail.

The authors show that the prefabricated cracks and unilateral unloading conditions not only deteriorate the mechanical strength but also have significant differences in failure modes. The degrading effect of cracks on model strength increases linearly with the decrease of the dip angle. Under the condition of true triaxial unilateral unloading, the deterioration effect of peak strength of rock is very significant, and unloading plays a role in promoting the instability failure of rock after peak, making the rock earlier instability failure. Associating with the particle vector diagram and crack coalescence process, the authors find that model failure mode under unilateral loading conditions is obviously distinct from that in triaxial loading. The peak strain in the unloading direction increases sharply, resulting in a new shear slip.

This study is expected to improve the understanding of the strength failure and cracking behavior of fractured rock under unilateral unloading.

Rock-fill dams are embankments of compacted free-draining granular earth containing an impervious zone. Earth utilized in such dams often contains a high percentage of large particles – hence the term rock-fill. Mass stability of these dams results from friction and particle interactions rather than through a cementing agent binding the particles together. However, high-stress conditions and prolonged exposure to the elements can severely damage rock-fill. Therefore, understanding and modeling rock-fill breakage is important for dam engineering. The purpose of this paper is to improve discontinuous deformation analysis (DDA) techniques for modeling rock-fill breakage, proving the new method using simulations of spherical particle crushing.

This work models rock-fill as bonded ellipsoid particles, and develops an improved DDA method to model the breakage of particle assemblies. The paper starts by describing the principles of three-dimensional DDA for spherical particles, and then derives the submatrices for normal contact, shear contact, and frictional force. The new algorithm incorporates a bond model with a revised open-close iteration algorithm into the DDA method to simulate particle crushing. To validate the improved DDA method, calculated particle contacts and movements are validated against theoretical results. Finally, this work performs a series of point-loading experimental tests for cement ellipsoid particles of both high and low compression strengths, with the test results compared against the results from corresponding DDA simulations.

In particle crushing tests, the force and displacement show an approximately linear relationship until the crushing point, at which point low compression ellipsoid particles split into several large pieces while the high-compression particles break into many small fragments. The DDA simulation results are in good agreement with the crushing tests, demonstrating the validity of the DDA method for solving particle crushing problems. Although the improved DDA model is applicable to rock-fill particle crushing studies, some issues remain, particularly in increasing calculation efficiency and performing large-scale computations and long real-time simulations. Future research should address these issues.

A bond model with a revised open-close iteration algorithm is incorporated into the DDA method. The simulated results shed insight into rock-fill crushing mechanisms, an element of concern in engineering practices.

To use distinct element simulation (PFC2D) to investigate the relationships between microparameters and macroproperties of the specimens that are modeled by bonded particles. To determine quantitative relationships between particle level parameters and mechanical properties of the specimens.

A combined theoretical and numerical approach is used to achieve the objectives. First, theoretical formulations are proposed for the relationships between microparameters and macroproperties. Then numerical simulations are conducted to quantify the relationships.

The Young's modulus is mainly determined by particle contact modulus and affected by particle stiffness ratio and slightly affected by particle size. The Poisson's ratio is mainly determined by particle stiffness ratio and slightly affected by particle size. The compressive strength can be scaled by either the bond shear strength or the bond normal strength depending on the ratio of the two quantities.

The quantitative relationships between microparameters and macroproperties for parallel‐bonded PFC2D specimens are empirical in nature. Some modifications may be needed to model a specific material. The effects of the particle distribution and bond strength distribution of a PFC2D specimen are very important aspects that deserve further investigation.

The results will provide guidance for people who use distinct element method, especially the PFC2D, to model brittle materials such as rocks and ceramics.

This paper offers some new quantitative relationships between microparameters and macroproperties of a synthetic specimen created using bonded particle model.

Granular materials possess multiscale structures, i.e. micro-scales involving atoms and molecules in a solid particle, meso-scales involving individual particles and their correlated structure, and macroscopic assembly. Strong and abundant dissipations are exhibited due to mesoscopic unsteady motion of individual grains, and evolution of underlying structures (e.g. force chains, vortex, etc.), which defines the key differences between granular materials and ordinary objects. The purpose of this paper is to introduce the major studies have been conducted in recent two decades.

The main properties at individual scale are introduced, including the coordination number, pair-correlation function, force and mean stress distribution functions, and the dynamic correlation function. The relationship between meso- and macro-scales is analyzed, such as between contact force and stress, the elastic modulus, and bulk friction in granular flows. At macroscales, conventional engineering models (i.e. elasto-plastic and hypo-plastic ones) are introduced. In particular, the so-called granular hydrodynamics theory, derived from thermodynamics principles, is explained.

On the basis of recent study the authors conducted, the multiscales (both spatial and temporal) in granular materials are first explained, and a multiscale framework is presented for the mechanics of granular materials.

It would provide a paramount view on the multiscale studies of granular materials.

The main purpose of this paper is to present and use the particle simulation method to explicitly simulate the spontaneous crack initiation phenomenon in brittle materials, and to compare the particle simulation results with experimental ones on the laboratory scale.

Using the particle simulation method, the brittle material is simulated as an assembly of particles so that the microscopic mechanism of inter‐ and intra‐particle crack initiation can be straightforwardly considered on the microscopic scale. A laboratory test has been conducted using a gypsum sample model to validate the particle simulation method for explicitly simulating the spontaneous crack initiation phenomenon.

The paper finds that in terms of simulating the macroscopic sliding surface along or around the contact plane between a block and its foundation, both the laboratory test and the particle simulation have produced consistent results. This indicated that the particle simulation method is capable of simulating macroscopic cracks through simulating conglomerations and accumulations of microscopic crack initiation within the brittle material. Compared with other numerical methods, the particle simulation method is more suitable for explicitly and effectively simulating spontaneous crack initiation problems on the microscopic scale in brittle materials.

The particle simulation method can be used to explicitly and effectively simulate the spontaneous crack initiation on the microscopic scale in brittle materials. It can be also used to simulate the macroscopic sliding surface along or around the contact plane between a block and its foundation. The experimental results of simulating the spontaneous crack initiation on the laboratory scale in brittle materials are very valuable for validating the numerical simulation results obtained not only from the particle simulation method, but also from other numerical simulation methods.

Accurate 3D experimental particle trajectory data, acquired from a laboratory tumbling mill using bi‐planar X‐ray filming, are used to validate the discrete element method (DEM). Novel numerical characterisation techniques are presented that provide a basis for comparing the experimental and simulated charge behaviour. These techniques are based on fundamental conservation principles, and provide robust, new interpretations of charge behaviour that are free of operator bias. Two‐ and three‐dimensional DEM simulations of the experimental tumbling mill are performed, and the relative merits of each discussed. The results indicate that in its current form DEM can simulate some of the salient features of the tumbling mill charge, however, comparison with the experiment indicate that the technique requires refinement to adequately simulate all aspects of the system.

The purpose of this paper is to present an investigation of the effect of different gravity conditions on the penetration mechanism using the two-dimensional Distinct Element Method (DEM), which ranges from high gravity used in centrifuge model tests to low gravity incurred by serial parabolic flight, with the aim of efficiently analyzing cone penetration tests on the lunar surface.

Seven penetration tests were numerically simulated on loose granular ground under different gravity conditions, i.e. one-sixth, one-half, one, five, ten, 15 and 20 terrestrial gravities. The effect of gravity on the mechanisms is examined with aspect to the tip resistance, deformation pattern, displacement paths, stress fields, stress paths, strain and rotation paths, and velocity fields during the penetration process.

First, under both low and high gravities, the penetration leads to high gradients of the value and direction of stresses in addition to high gradients in the velocity field near the penetrometer. In addition, the soil near the penetrometer undergoes large rotations of the principal stresses. Second, high gravity leads to a larger rotation of principal stresses and more downward particle motions than low gravity. Third, the tip resistance increases with penetration depth and gravity. Both the maximum (steady) normalized cone tip resistance and the maximum normalized mean (deviatoric) stress can be uniquely expressed by a linear equation in terms of the reciprocal of gravity.

This study investigates the effect of different gravity conditions on penetration mechanisms by using DEM.