An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.

Multi-layer perceptron-based data-driven multiscale modelling of granular materials with a novel Frobenius norm-based internal variable

One objective of developing machine learning(ML)-based material models is to integrate them with well-established numerical methods to solve boundary value problems(BVPs).In the family of ML models,recurrent neural networks(RNNs)have been extensively applied to capture history-dependent constitutive responses of granular materials,but these multiple-step-based neural networks are neither sufficiently efficient nor aligned with the standard finite element method(FEM).Single-step-based neural networks like the multi-layer perceptron(MLP)are an alternative to bypass the above issues but have to introduce some internal variables to encode complex loading histories.In this work,one novel Frobenius norm-based internal variable,together with the Fourier layer and residual architectureenhanced MLP model,is crafted to replicate the history-dependent constitutive features of representative volume element(RVE)for granular materials.The obtained ML models are then seamlessly embedded into the FEM to solve the BVP of a biaxial compression case and a rigid strip footing case.The obtained solutions are comparable to results from the FEM-DEM multiscale modelling but achieve significantly improved efficiency.The results demonstrate the applicability of the proposed internal variable in enabling MLP to capture highly nonlinear constitutive responses of granular materials.

A refined deviatoric hardening plastic model for sand

The classical deviatoric hardening models are capable of characterizing the mechanical response of granular materials for a broad range of degrees of compaction.This work finds that it has limitations in accurately predicting the volumetric deformation characteristics under a wide range of confining/consolidation pressures.The issue stems from the pressure independent hardening law in the classical deviatoric hardening model.To overcome this problem,we propose a refined deviatoric hardening model in which a pressure-dependent hardening law is developed based on experimental observations.Comparisons between numerical results and laboratory triaxial tests indicate that the improved model succeeds in capturing the volumetric deformation behavior under various confining/consolidation pressure conditions for both dense and loose sands.Furthermore,to examine the importance of the improved deviatoric hardening model,it is combined with the bounding surface plasticity theory to investigate the mechanical response of loose sand under complex cyclic loadings and different initial consolidation pressures.It is proved that the proposed pressure-dependent deviatoric hardening law is capable of predicting the volumetric deformation characteristics to a satisfactory degree and plays an important role in the simulation of complex deformations for granular geomaterials.

Experimental study on workability and permeability of sandy soils conditioned with thickened foam

Water spewing and muck plugging often occur during earth pressure balance(EPB)shield machines tunnelling in water-rich sandy strata,even though the conventional foam has been employed to condition sandy soils.In this study,a novel thickened foaming agent suitable for EPB shield tunnelling in water-rich sandy strata is developed.In contrast to conventional foam-conditioned sands,the thickened foam-conditioned sand has a low permeability due to the consistent filling of soil pores with the thickened foam,and the initial permeability coefficient decreases by approximately two orders of magnitude.It also exhibits a suitable workability,which is attributed to the enhanced capability of the thickened foam to condition sandy soils.In addition,the effect of concentration on the stability of the foam is explained by the Gibbs-Marangoni effect,and conditioning mechanisms for the thickened foam on sands are discussed from the evolution of foam bubbles.

The role of polyurethane foam compressible layer in the mechanical behaviour of multi-layer yielding supports for deep soft rock tunnels

The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.

Deep Learning Predicts Stress–Strain Relations of Granular Materials Based on Triaxial Testing Data

This study presents an AI-based constitutive modelling framework wherein the prediction model directly learns from triaxial testing data by combining discrete element modelling(DEM)and deep learning.A constitutive learning strategy is proposed based on the generally accepted frame-indifference assumption in constructing material constitutive models.The low-dimensional principal stress-strain sequence pairs,measured from discrete element modelling of triaxial testing,are used to train recurrent neural networks,and then the predicted principal stress sequence is augmented to other high-dimensional or general stress tensor via coordinate transformation.Through detailed hyperparameter investigations,it is found that long short-term memory(LSTM)and gated recurrent unit(GRU)networks have similar prediction performance in constitutive modelling problems,and both satisfactorily predict the stress responses of granular materials subjected to a given unseen strain path.Furthermore,the unique merits and ongoing challenges of data-driven constitutive models for granular materials are discussed.

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method(DEM)has been intensively used to study the constitutive behaviour of granular materials.However,to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear.This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials.Three groups of models with spherical and clumped particles are investigated from four perspectives:(i)deviatoric stress and volumetric behaviour;(ii)critical state behaviour;(iii)stress-dilatancy relationship;and(iv)the evolution of principal stress ratio against axial strain.The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses,volumetric behaviour,and critical state behaviour observed in experiments for granular materials.On the other hand,some qualitative deviations between experiments and the investigated DEM simulations are also observed,in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain,which are proven to be critical for constitutive modelling.The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective.It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

Undrained vane shear strength of sand-foam mixtures subjected to different shear rates

The shear strength of sand-foam mixtures plays a crucial role in ensuring successful earth pressure balance(EPB)shield tunneling.Since the sand-foam mixtures are constantly sheared by the cutterhead and the screw conveyor with varied rotation speeds during tunneling,it is non-trivial to investigate the effect of shear rates on the undrained shear strength of sand-foam mixtures under chamber pressures to extend the understanding on the tunneling process.This study conducted a series of pressurized vane shear tests to investigate the role of shear rates on the peak and residual strengths of sand-foam mixtures at different pore states.Different from the shear-rate characteristics of natural sands or clay,the results showed that the peak strength of sand-foam mixtures under high vertical total stress(σ_(v)≥200 kPa)and low foam injection ratio(FIR30%)decreased with the increase in shear rate.Otherwise,the peak strength was not measurably affected by shear rates.The sand-foam mixtures in the residual state resembled low-viscous fluid with yield stress and the residual strength increased slightly with shear rates.In addition,the peak and residual strengths were approximately linear with vertical effective stress regardless of the total stress and FIR.The peak effective internal friction angle remained almost invariant in a low shear rate(γ′<0.25 s1)but decreased when the shear rate continued increasing.The residual effective internal friction angle was lower than the peak counterpart and insensitive to shear rates.This study unveiled the role of shear rates in the undrained shear strength of sand-foam mixtures with various FIRs and vertical total stresses.The findings can extend the understanding of the rate-dependent shear characteristics of conditioned soils and guide the decision-making of soil conditioning schemes in the EPB shield tunneling practice.

Influence factors and calculation model of the adhesion strength of clayey soil for EPB shield tunnelling

When earth pressure balance (EPB) shield tunnels are constructed through clayey ground, the soil adheres to the cutter, cutterheadand chamber bulkhead due to the high adhesion strength between the steel and the clay, thus clogging the shield. To investigatethe influence of different factors on the adhesion strength, this study used montmorillonite, kaolin and mixtures of the two as testsoils. The adhesion strength between the steel and clay is determined with a customized rotary shear apparatus. The results showthat when the consistency index of the soil specimen is less than 1, the adhesion strength between the steel and clay increaseswith the consistency index. As the consistency index decreases, the effect of the normal pressure on the adhesion strength graduallyweakened. As the contact angle of the shear plate increases, thus reducing the hydrophilicity, the adhesion strength decreases. Whenthe soil specimens with different plasticity index values have the same consistency index value, the adhesion strengths are similar.The adhesion strength increases gradually with increasing surface roughness. Based on grey incidence analysis, the order of thefactors affecting the adhesion strength is as follows: normal pressure > consistency index > contact angle > plasticity index ≈ surfaceroughness. The normal pressure, consistency index and contact angle all have important effects on the adhesion strength betweensteel and clay. However, because there are no large differences in the contact angle among the metal materials, the approximateadhesion strength can be calculated by considering the effects of only the normal pressure and consistency index. Themeasures thatare effective for preventing EPB shield clogging are increasing the soil softness and decreasing the hydrophilicity of the cutterheadby applying new materials for the cutters and cutterhead.