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 thermodynamics-based three-scale constitutive model for partially saturated granular materials

A three-scale constitutive model for unsaturated granular materials based on thermodynamic theory is presented.The three-scale yield locus,derived from the explicit yield criterion for solid matrix,is developed from a series of discrete interparticle contact planes.The three-scale yield locus is sensitive to porosity changes;therefore,it is reinterpreted as a corresponding constitutive model without phenomenological parameters.Furthermore,a water retention curve is proposed based on special pore morphology and experimental observations.The features of the partially saturated granular materials are well captured by the model.Under wetting and isotropic compression,volumetric compaction occurs,and the degree of saturation increases.Moreover,the higher the matric suction,the greater the strength,and the smaller the volumetric compaction.Compared with the phenomenological Barcelona basic model,the proposed three-scale constitutive model has fewer parameters;virtually all parameters have clear physical meanings.

Evaluation of the coefficient of lateral stress at rest of granular materials under repetitive loading conditions

Although the internal stress state of soils can be affected by repetitive loading,there are few studies evaluating the lateral stress(or K_(0))of soils under repetitive loading.This study investigates the changes in K_(0) and directional shear wave velocity(V_(s))in samples of two granular materials with different particle shapes during repetitive loading.A modified oedometer cell equipped with bender elements and a diaphragm transducer was developed to measure the variations in the lateral stress and the shear wave velocity,under repetitive loading on the loading and unloading paths.The study produced the following results:(1)Repetitive loading on the loading path resulted in an increase in the K_(0) of test samples as a function of cyclic loading number(i),and(2)Repetitive loading on the unloading path resulted in a decrease in K_(0) according to i.The shear wave velocity ratio(i.e.V_(s)(HH)/V_(s)(VH),where the first and second letters in parentheses corresponds to the directions of wave propagation and particle motion,respectively,and V and H corresponds to the vertical and horizontal directions,respectively)according to i supports the experimental observations of this study.However,when the tested material was in lightly over-consolidated state,there was an increase in K_(0) during repetitive loading,indicating that it was the initial K_(0),rather than the loading path,which is responsible for the change in K_(0).The power model can capture the variation in the K_(0) of samples according to i.Notably,the K_(0)=1 line acts as the boundary between the increase and decrease in K_(0) under repetitive loading.

Finite Element Simulations on Failure Behaviors of Granular Materials with Microstructures Using a Micromechanics-Based Cosserat Elastoplastic Model

This paper presents a micromechanics-based Cosserat continuum model for microstructured granular materials.By utilizing this model,the macroscopic constitutive parameters of granular materials with different microstructures are expressed as sums of microstructural information.The microstructures under consideration can be classified into three categories:a medium-dense microstructure,a dense microstructure consisting of one-sized particles,and a dense microstructure consisting of two-sized particles.Subsequently,the Cosserat elastoplastic model,along with its finite element formulation,is derived using the extended Drucker-Prager yield criteria.To investigate failure behaviors,numerical simulations of granular materials with different microstructures are conducted using the ABAQUS User Element(UEL)interface.It demonstrates the capacity of the proposed model to simulate the phenomena of strain-softening and strain localization.The study investigates the influence of microscopic parameters,including contact stiffness parameters and characteristic length,on the failure behaviors of granularmaterials withmicrostructures.Additionally,the study examines themesh independence of the presented model and establishes its relationship with the characteristic length.A comparison is made between finite element simulations and discrete element simulations for a medium-dense microstructure,revealing a good agreement in results during the elastic stage.Somemacroscopic parameters describing plasticity are shown to be partially related to microscopic factors such as confining pressure and size of the representative volume element.

Modelling of the variation of granular base materials resilient modulus with material characteristics and humidity conditions

This study aims to quantify the susceptibility of granular materials used in pavements to changes in moisture content and propose a correlation model to incorporate this susceptibility into seasonal analyses.The fines content and the percentage of fractured coarse aggregates were identified as direct indicators of the resilient modulus susceptibility to changes in water content.The results showed that the percentage of fractured coarse aggregates particles(FR)has a more significant impact on the resilient modulus(Er)of crushed granular materials used in pavement construction than the combined indicator of the fines content and sample volumetrics(nf).Crushed granular materials with a higher percentage of fractured coarse aggregates are relatively insensitive to changes in the degree of saturation,but become more sensitive as the fine fraction porosity decreases.An adjusted model was proposed based on the existing formulation,but considers a complex parameter to describe and adjust the sensitivity of base granular materials to variations in moisture content with respect to fabrication charac-teristics,fines content and volumetric properties.The model shows that the variation of Er values is below10%for fully crushed granular materials.However,it reaches approximately±12%for materials with 75%of crushed coarse aggregates andþ40%and-25%for materials with FR=50%.This model could help select good ag-gregates characteristics and adjust grain-size distribution for environments where significant moisture content variations can occur in the pavement system,such as in the Province of Quebec(Canada).As it is based on pa-rameters that can be easily determined or estimated,it also represents a valuable tool for detailed design and analysis that can consider material characteristics.

Constitutive modelling of idealised granular materials using machine learning method

Predicting the constitutive response of granular soils is a fundamental goal in geomechanics.This paper presents a machine learning(ML)framework for the prediction of the stress-strain behaviour and shearinduced contact fabric evolution of an idealised granular material subject to triaxial shearing.The MLbased framework is comprised of a set of mini-triaxial tests which provide a benchmark for the setup and validation of the discrete element method(DEM)model of the granular materials,a parametric DEM simulation programme of virtual triaxial tests which provides datasets of micro-and macro-mechanical information,as well as a multi-layer perceptron(MLP)neural network which is trained and tested using the DEM-based datasets.The ML model only requires the initial void ratio of the granular sample as the input for predicting its constitutive response.The excellent agreement between the ML model prediction and experimental test and DEM simulation results indicates that the MLebased modelling approach is capable of capturing accurately the effects of initial void ratio on the constitutive behaviour of idealised granular materials,bypassing the need to incorporate the complex micromechanics underlying the macroscopic mechanical behaviour of granular materials.Lastly,a detailed comparison between the used MLP model and long short-term memory(LSTM)model was made from the perspective of technical algorithm,prediction accuracy,and computational efficiency.

The Effects of the Particle Size Ratio on the Behaviors of Binary Granular Materials

The particle size ratio(PSR)is an important parameter for binary granular materials,which may aect the microstructure and macro behaviors of granular materials.However,the eect of particle ratio on granular assemblies with dierent arrangements is still unclear.To explore and further clarify the eect of PSR in dierent packing structures,three types of numerical samples with regular,layered,and random packing are designed.Numerical results show that PSR has signi􀀀cant eects on binary granular samples with regular packing.The larger the PSR,the stronger the strength,the larger the modulus,and the smaller the angle between the shear band and the load direction.And a theoretical solution of the peak stress ratio vs.PSR is obtained for regular packing,and the results by DEM are in good agreement with the theoretical solution.Under layered packing,PSR has little eect on peak stress ratio due to similar microstructure obtained with the changing of PSR.The modulus slightly increased with the increase of PSR.Under random packing with small grain content of 50%,PSR has little eect in the range of 0.5–0.9,but in a larger range,larger PSR leads to greater modulus.

Granular behaviour under bi-directional shear with constant vertical stress and constant volume

This paper aims to investigate the role of bi-directional shear in the mechanical behaviour of granular materials and macro-micro relations by conducting experiments and discrete element method(DEM)modelling.The bi-directional shear consists of a static shear consolidation and subsequent shear under constant vertical stress and constant volume conditions.A side wall node loading method is used to exert bi-directional shear of various angles.The results show that bi-directional shear can significantly influence the mechanical behaviour of granular materials.However,the relationship between bidirectional shear and mechanical responses relies on loading conditions,i.e.constant vertical stress or constant volume conditions.The stress states induced by static shear consolidation are affected by loading angles,which are enlarged by subsequent shear,consistent with the relationship between bidirectional shear and principal stresses.It provides evidence for the dissipation of stresses accompanying static liquefaction of granular materials.The presence of bi-directional principal stress rotation(PSR)is demonstrated,which evidences why the bi-directional shear of loading angles with components in two directions results in faster dissipations of stresses with static liquefaction.Contant volume shearing leads to cross-anisotropic stress and fabric at micro-contacts,but constant vertical stress shearing leads to complete anisotropic stress and fabric at micro-contacts.It explains the differentiating relationship between stress-strain responses and fabric anisotropy under these two conditions.Micromechanical signatures such as the slip state of micro-contacts and coordination number are also examined,providing further insights into understanding granular behaviour under bi-directional shear.

Fluidization and Transport of Vibrated Granular Matter:A Review of Landmark and Recent Contributions

We present a short retrospective review of the existing literature about the dynamics of(dry)granular matter under the effect of vibrations.The main objective is the development of an integrated resource where vital information about past findings and recent discoveries is provided in a single treatment.Special attention is paid to those works where successful synthetic routes to as-yet unknown phenomena were identified.Such landmark results are analyzed,while smoothly blending them with a history of the field and introducing possible categorizations of the prevalent dynamics.Although no classification is perfect,and it is hard to distillate general properties out of specific observations or realizations,two possible ways to interpret the existing results are defined according to the type of forcing or the emerging(ensuing)regime of motion.In particular,first results concerning the case where vibrations and gravity are concurrent(vertical shaking)are examined,then the companion situation with vibrations perpendicular to gravity(horizontal shaking)is described.Universality classes are introduced as follows:(1)Regimes where sand self-organizes leading to highly regular geometrical“pulsating”patterns(thin layer case);(2)Regimes where the material undergoes“fluidization”and develops an internal multicellular convective state(tick layers case);(3)Regimes where the free interface separating the sand from the overlying gas changes inclination or develops a kind a patterned configuration consisting of stable valleys and mountains or travelling waves;(4)Regimes where segregation is produced,i.e.,particles of a given size tend to be separated from the other grains(deep containers).Where possible,an analogy or parallelism is drawn with respect to the companion field of fluid-dynamics for which the assumption of“continuum”can be applied.

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.

Effect of granular shape on radial segregation in a two-dimensional drum

Granular segregation is widely observed in nature and industry.Most research has focused on segregation caused by differences in the size and density of spherical grains.However,due to the fact that grains typically have different shapes,the focus is shifting towards shape segregation.In this study,experiments are conducted by mixing cubic and spherical grains.The results indicate that spherical grains gather at the center and cubic grains are distributed around them,and the degree of segregation is low.Through experiments,a structured analysis of local regions is conducted to explain the inability to form stable segregation patterns with obviously different geometric shapes.Further,through simulations,the reasons for the central and peripheral distributions are explained by comparing velocities and the number of collisions of the grains in the flow layer.

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.

New double yield surface model for coarse granular materials incorporating nonlinear unified failure criterion

A new double-yield-sarface （DYS） model was developed to characterize the strength and deformation behaviors of coarse granular materials （CGMs）. Two kinds of deformation mechanisms, including the shear and compressive plastic deformation, were taken into account in this model, These two deformation mechanisms were described by the shear and compressive yield functions, respectively. The Lode angle dependent formulations of proposed model were deduced by incorporating a 3D nonlinear unified failure criterion. Some comparisons were presented between the numerical predictions of proposed model and test data of true triaxial tests on the modeled rockfills. The model predictions are in good agreement with the test data and capture the strain hardening and plastic volumetric dilation of CGMs. These findings verify the reasonability of current DYS model, and indicate that this model is well suited to reproduce the stress-strain-volume change behavior of CGMs in general.

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.

Investigations of the effects of particle morphology on granular material behaviors using a multi-sphere approach

This article studies the influences of particle morphology on the behaviors of granular materials at both macroscopic and microscopic levels based on the discrete element method(DEM).A set of numerical tests under drained triaxial compression was performed by controlling two morphological descriptors,i.e.ratio of the smallest to the largest pebble diameter,x,and the maximum pebbleepebble intersection angle,b.These descriptors are vital in generating particle geometry and surface textures.It was found that the stress responses of all assemblies exhibited similar behavior and showed post-peak strainsoftening.The normalized stress ratio and volumetric strains flatten off and tended to reach a steady value after an axial strain of 40%.While the friction angles at peak state varied with different morphological descriptors,the friction angles at critical state showed no significant variation.Moreover,evolution of the average coordination numbers showed a dramatic exponential decay until an axial strain of about 15%after which it stabilized and was unaffected by further increase of axial strain.In addition,stress ratio q/p and strong fabric parameter fs d=fs m were found to follow an approximately linear relationship for each assembly.These findings emphasized the significance of the influences of particle morphology on the macroscopic and microscopic responses of granular materials.

Probabilistic analysis of random contact force between geomembrane and granular material

A probabilistic method based on principle of maximum entropy was employed to analyze the randomness of contact force between geomembrane and granular material.The contact force distribution is exponential according to the proposed method and the grain size is the most important factor that affects the distribution of contact force.The proposed method is then verified by a series of laboratory experiments using glass beads and cobbles as granular material and a very thin pressure,indicating that film is firstly used in these experiments which give a reliable method to measure the contact force at each contact point.

Experimental and numerical study of the blast wave decrease using sandwich panel by granular materials core

Among the intrinsic properties of some materials,e.g.,foams,porous materials,and granular materials,are their ability to mitigate shock waves.This paper investigated shock wave mitigation by a sandwich panel with a granular core.Numerical simulations and experimental tests were performed using Autodyn hydro-code software and a shock tube,respectively.The smoothed particle hydrodynamics(SPH)method was used to model granular materials.Sawdust and pumice,whose properties were determined by several compression tests,were used as granular materials in the sandwich panel core.These granular materials possess many mechanisms,including compacting(e.g.,sawdust)and crushing(e.g.,pumice)that mitigate shock/blast wave.The results indicated the ineffectiveness of using a core with low thickness,yet it was demonstrated to be effective with high thickness.Low-thickness pumice yielded better results for wave mitigation.The use of these materials with a core with appropriate core reduces up to 88%of the shock wave.The results of the experiments and numerical simulations were compared,suggesting a good agreement between the two.This indicates the accuracy of simulation and the ability of the SPH method to modeling granular material under shock loading.The effects of grain size and the coefficient of friction between grains have also been investigated using simulation,implying that increasing the grain size and coefficient of friction between grains both reduce overpressure.

The Material Deformation and Internal Structure Development of Granular Materials under Different Cyclic Loadings

Common structures in engineering such as slopes,roadbeds,ballasts,etc.,are closely related to granular materials.They are usually subjected to long-term cyclic loads.This study mainly focused on the mechanical behaviors of randomly arranged granular materials before they reach a stable state under different cyclic loads.The variation of the maximum axial strain and the influence of CSR(cyclic stress ratio)were analyzed.The energy consumed in each cycle under constant confining stress loading condition is significantly greater than that of the fixed wall loading condition.The internal deformation evolution of granular materials is studied in detail.The deformation mode of granular material under cyclic loading at different positions inside the material is different according to the strain variation.In addition,the strain,force chain structure and contact force magnitude are combined to explore their effects on local deformation of granular materials under cyclic loading.From the perspective of the deformation form,the material sample can be divided into several regions,and the ability to adjust particle positions determines the deformation mode of different regions.The changes of local strain with the cyclic loading also reflect the contribution of particle displacements to the evolution ofmicrostructure.This research will provide insights into the understanding of granular materials behaviors under cyclic loading.

Characterization of Elastic Modulus of Granular Materials in a New Designed Uniaxial Oedometric System

A simple uniaxial oedometric system is developed to test the elastic modulus of granular materials. The stress- strain relationship is first measured under conditions of uniaxial compression with additional lateral stress and strain, then the elastic modulus of the material is determined by the linear fitting method. It is found that the modulus is positively correlated to the grain size and negatively correlated to the container size. Arching and dragging are revealed to be the mechanism of such correlations by using the digital image correlation method and the pressure film technology based on the statistical method.

Numerical implementation of suction-dependent resilient modulus constitutive model for subgrade granular material

In order to investigate the suction-dependent properties of subgrade granular material and its effect on pavement responses,coupled hydro-mechanical simulations were conducted in Abaqus.A suction-dependent resilient modulus model was integrated into the commercial finite element(FE)code Abaqus by developing a user-defined material(UMAT)subroutine.The developed model was validated by triaxial test results under different suction conditions and good agreement was achieved.A three-dimensional(3D)FE pavement model was established and the suction-dependent properties of subgrade granular material was characterized by the developed constitutive model.Hydro-mechanical pavement responses subjected to three moisture states and the falling weight deflectometer(FWD)load were calculated.Simulation results reveal that the resilient modulus of subgrade granular material is sensitive to suction and stress states;high groundwater table decreases the overall resilient moduli of subgrade structure due to suction reduction,leading to the increase of the maximum surface deflection,the tensile strain at bottom of the surface layer,compressive strain on top of subgrade,and consequently,deterioration in pavement performance.