Search results

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.

Cemented conglomerate accumulation is a weak and heterogeneous medium that occurs in western China. It consists mainly of argillaceous cement that loses strength rapidly upon contact with water, leading to collapse instability failure. Its deformation failure mechanism is complex and poorly understood. In this paper, the erosion failure mechanism of cemented conglomerate accumulation is investigated.

The collapse failure process after erosion of the slope foot for typical cemented conglomerate accumulation is studied based on field investigation using the particle discrete element method. And how the medium composition, slope angle and cementation degree influence the failure mode and process of the cemented conglomerate accumulation is examined.

The foot erosion of slope induces a tensile failure that typically manifests as “erosion at the foot of slope – tensile cracking at the back edge of slope top – integral collapse.” The collapse failure is more likely to occur when the cemented conglomerate accumulation has a higher rock content, a steeper slope angle or a weaker cementation degree.

A model based on rigid blocks and disk particles to simulate the cemented conglomerate accumulation is developed. It shows that the hydraulic erosion at the foot of the slope resulted in a different failure mechanism than that of general slopes. The results can inform the stability management, disaster prevention and mitigation of similar slopes.

The objective of this paper is to quantitatively assess shear band evolution by using two-dimensional discrete element method (DEM).

The DEM model was first calibrated by retrospectively modelling existing triaxial tests. A series of DEM analyses was then conducted with the focus on the particle rotation during loading. An approach based on particle rotation was developed to precisely identify the shear band region from the surrounding. In this approach, a threshold rotation angle ω₀ was defined to distinguish the potential particles inside and outside the shear band and an index g(ω₀) was introduced to assess the discrepancy between the rotation response inside and outside shear band. The most distinct shear band region can be determined by the ω₀ corresponding to the peak g(ω₀). By using the proposed approach, the shear band development of two computational cases with different typical localised failure patterns were successfully examined by quantitatively measuring the inclination angle and thickness of shear band, as well as the microscopic quantities.

The results show that the shear band formation is stress-dependent, transiting from conjugated double shear bands to single shear band with confining stress increasing. The shear band evolution of two typical localised failure modes exhibits opposite trends with increasing strain level, both in inclination angle and thickness. Shear band featured a larger volumetric dilatancy and a lower coordination number than the surrounding. The shear band also significantly disturbs the induced anisotropy of soil.

This paper proposed an approach to quantitatively assess shear band evolution based on the result of two-dimensional DEM modelling.

The study aims to investigate the effects of pre-loading histories (pre-shearing and pre-consolidation) on the liquefaction behaviour of saturated loose sand via discrete element method (DEM) simulations.

The pre-shearing history is mimicked under drained conditions (triaxial compression) with different pre-shearing strain levels ranging from 0% to 2%. The pre-consolidation history is mimicked by increasing the isotropic compression to different levels ranging from 100 kPa to 300 kPa. The macroscopic and microscopic behaviours are analysed and compared.

Temporary liquefaction, or quasi-steady state (QSS), is observed in most samples. A higher pre-shearing or pre-consolidation level can provide higher liquefaction resistance. The ultimate state line is found to be unique and independent of the pre-loading histories in stress space. The Lade instability line prematurely predicts the onset of liquefaction for all samples, both with and without pre-loading histories. The redundancy index is an effective microscopic indicator to monitor liquefaction, and the onset of the liquefaction corresponds to the phase transition state where the value of redundancy index is one, which is true for all cases irrespective of the proportions of sliding contacts.

The liquefaction behaviour of granular materials still remains elusive, especially concerning the effects of pre-loading histories on soils. Furthermore, the investigation of the effects of pre-consolidation histories on undrained behaviour and its comparison to pre-sheared samples is rarely reported in the DEM literature.

The paper aims to find the correlation between the microparameters and the macroparameters of the soil. The study aims to calibrate the macroscale and microscale parameters of rolling resistance contact models to successfully apply the discrete element method to do some research of the geotechnical problem.

The paper opted for an exploratory study using the PFC3D to simulate the triaxial tests that include more than 50 cases and the coupling analysis method, which considering several effect of various factors.

The paper provides a quantitative relationship between the macroparameters and microparameters of the rolling resistance linear model and a method for fast calibration of macroscopic parameters is proposed and verified by a triaxial test example.

This paper provides the quantitative relationship of micro and macroparameters in the rolling resistance linear model by studying a single factor and considering the coupling effect of various factors and a fast method for the calibration of microparameters based on the rolling resistance linear model is proposed.

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.

The research in this paper helps to understand the difference between the Eulerian method and the Lagrangian method in describing the performance of Pelton turbine buckets, so as to improve the design level and design efficiency of the runner.

This paper used DualSPHysics to calculate the unsteady flow of the Pelton turbine runner bucket and compared it with the mesh-based method to explore the difference between mesh-based and particle-based methods in torque curves, jet flow patterns and pressure characteristics.

It is noted that the particle-based method is challenging to compare with the mesh-based method concerning accuracy. In addition to better describing the free water film, the particle method also captures many droplets near the water film, but it cannot well describe the negative pressure region on the bucket back and the resulting jet interference after cutting off the jet. Compared with the mesh-based method, the pressure measurement points obtained by the particle-based method generally have shorter periods and violent fluctuations, and the pressure value of some points is underestimated.

This paper helped to calculate the unsteady characteristics of the Pelton turbine by Fluent, CFX and DualSPHysics; exploration jet flow pattern differences between the mesh and meshfree methods; prediction of the flow interference between the bucket back and the jet and the pressure curve of SPH usually has a shorter period and violent fluctuations.

This study aims to propose optimised function-based evolutionary algorithms in this research to effectively replace the traditional electronic circuitry used in linearising constant temperature anemometer (CTA) and Microbridge mass flow sensor AWM 5000.

The proposed linearisation technique effectively uses the ratiometric function for the linearisation of CTA and Microbridge mass flow sensor AWM 5000. In addition, the well-known transfer relation, namely, the King’s Law is used for the linearisation of CTA and successfully implemented using LabVIEW 7.1.

Investigational results unveil that the proposed evolutionary optimised linearisation technique performs better in linearisation of both CTA and Mass flow sensors, and hence finds applications for computer-based flow measurement/control systems.

The evolutionary optimisation algorithms such as the real-coded genetic algorithm, particle swarm optimisation algorithm, differential evolution algorithm and covariance matrix adopted evolutionary strategy algorithm are used to determine the optimal values of the parameters present in the proposed ratiometric function. The performance measures, namely, the full-scale error and mean square error are used to analyse the overall performance of the proposed approach is compared to a state of art techniques available in the literature.

This paper aims to present a novel approach for computing particle temperatures in simulations coupling computational fluid dynamics (CFD) and discrete element method (DEM) to predict flow and heat transfer in fluidized beds of thermally thick spherical particles.

An improved lumped formulation based on Hermite-type approximations for integrals to relate surface temperature to average temperature and surface heat flux is used to overcome the limitations of classical lumped models. The model is validated through comparisons with analytical solutions for a convectively cooled sphere and experimental data for a fixed particle bed. The coupled CFD-DEM model is then applied to simulate a Geldart D bubbling fluidized bed, comparing the results to those obtained using the classical lumped model.

The validation cases demonstrate that ignoring internal thermal resistance can significantly impact the temperature in cases where the Biot number is greater than 0.1. The results for the fixed bed case clearly demonstrate that the proposed method yields significantly improved outcomes compared to the classical model. The fluidized bed results show that surface temperature can deviate considerably from the average temperature, underscoring the importance of accurately accounting for surface temperature in convective heat transfer predictions and surface processes.

The proposed approach offers a physically more consistent simulation without imposing a significant increase in computational cost. The improved lumped formulation can be easily and inexpensively integrated into a typical DEM solver workflow to predict heat transfer for spherical particles, with important implications for various industrial applications.

Particles bed binding by selective cement activation (SCA) method is a computer-aided manufacturing (CAM) technique used to produce cementitious elements. A computer-aided design file is sliced to generate G-codes before printing. This paper aims to study the effect of key input parameters for slicer software on the final properties of printed products.

The one factor at a time (OFAT) methodology is used to investigate the impact of selected parameters on the final properties of printed specimens, and the causes for the variations in outcomes of each variable are discussed.

Finer aggregates can generate a more compact layer, resulting in a denser product with higher strength. Fluid pressure is directly determined by voxel rate (r_V); however, high pressures enable better fluid penetration control for fortified products; for extreme r_Vs, residual voids in the interfaces between successive layers and single-line primitives impair mechanical strength. It was understood that printhead movement along the orientation of the parts in the powder bed improved the mechanical properties.

The design of experiment (DOE) method assesses the influence of process parameters on various input printing variables at the same time. As the resources are limited, a fractional factorial plan is carried out on a subset of a full factorial design; hence, providing physical interpretation behind changes in each factor is difficult. OFAT aids in analyzing the effect of a change in one factor on output while all other parameters are kept constant. The results assist engineers in properly considering the influence of variable variations for future DOE designs.