Second-order work analyses, based on elasto-plastic models, have been frequently carried out leading to the result that failure may occur the limit yield condition is encountered. In this article, second-order work investigations are carried out with barodesy regarding standard element tests and finite element applications. In barodesy, it was shown-like in hypoplasticity and elasto-plasticity-that second-order work may vanish at stress states inside the critical limit surface. For boundary value problems, an end-to-end shear band of vanishing second-order work marks situations, where failure is imminent.

The intergranular strain concept (IGS) and intergranular strain anisotropy formulation (ISA) are state of the art extensions to describe small-strain effects. The main conceptional difference between ISA and IGS is the purely elastic strain range introduced by ISA. In addition, the ISA formulation used in this article includes an additional state variable in order to reduce accumulation effects for cyclic loading with a larger number of repetitive cycles. Barodesy is enhanced here with ISA to improve its small-strain predictions. The performance of this new model is compared with barodesy enhanced with IGS. It turned out that the small-strain extensions do not negatively influence predictions under monotonic loading. Differences between ISA and ISG are only remarkable for very small-strain cycles and even there they are negligible for certain parameter values.

A data-driven framework was used to predict the macroscopic mechanical behavior of dense packings of polydisperse granular materials. The discrete element method, DEM, was used to generate 92,378 sphere packings that covered many different kinds of particle size distributions, PSD, lying within 2 particle sizes. These packings were subjected to triaxial compression and the corresponding stress-strain curves were fitted to Duncan-Chang hyperbolic models. An artificial neural network (NN) scheme was able to anticipate the value of the model parameters for all these PSDs, with an accuracy similar to the precision of the experiment and even when the NN was trained with a few hundred DEM simulations. The estimations were indeed more accurate than those given by multiple linear regressions (MLR) between the model parameters and common geotechnical and statistical descriptors derived from the PSD. This was achieved in spite of the presence of noise in the training data. Although the results of this massive simulation are limited to specific systems, ways of packing and testing conditions, the NN revealed the existence of hidden correlations between PSD of the macroscopic mechanical behavior.

Soil flow is involved in many earth surface processes such as debris flows and landslides. It is a very challenging task to model this large deformational phenomenon because of the extreme change in material configurations and properties when soil flows. Most of the existing models require a two-dimensional (2D) simplification of actual systems, which are however three-dimensional (3D). To overcome this issue, we develop a novel 3D particle finite element method (PFEM) for direct simulation of complex soil flows in 3D space. Our PFEM model implemented in a fully implicit solution framework based on a generalised Hellinger-Reissner variational principle permits the use of a large time step without compromising the numerical stability. A mixed quadratic-linear element is used to avoid volumetric locking issues and ensure computational accuracy. The correctness and robustness of our 3D PFEM formulation for modelling large deformational soil flow problems are demonstrated by a series of benchmarks against analytical or independent numerical solutions. Our model can serve as an effective tool to support the assessment of catastrophic soil slope failures and subsequent runout behaviours.

Sustainable biopolymer additives offer a promising soil stabilisation methodology, with a strong potential to be tuned to soil's specific nature, allowing the tailoring of mechanical properties for a range of geotechnical applications. However, the biopolymer chemical characteristics driving soil mechanical property modifications have yet to be fully established. Within this study we employ a cross-scale approach, utilising the differing galactose:mannose (G:M) ratios of various Galactomannan biopolymers (Guar Gum G:M 1:2, Locust Bean Gum G:M 1:4, Cassia Gum G:M 1:5) to investigate the effect of microscale chemical functionality upon macroscale soil mechanical properties. Molecular weight effects are also investigated, utilising Carboxy Methyl Cellulose (CMC). Soil systems comprising of SiO (100%) (SiO) and a Mine Tailing (MT) exemplar composed of SiO (90%) + FeO (10%) (SiO + Fe) are investigated. The critical importance of biopolymer additive chemical functionality for the resultant soil mechanical properties, is demonstrated.For Galactomannan G:M 1:5 stabilised soils the 'high-affinity, high-strength', mannose-Fe interactions at the microscale (confirmed by mineral binding characterisation) are attributed to the 297% increase in the SiO + Fe systems Unconfined Compressive Strength (UCS), relative to SiO only. Conversely for SiO Galactomannan-stabilised soils, when increasing the G:M ratio from 1:2 to 1:5, a 85% reduction in UCS is observed, attributed to mannose's inability to interact with SiO. UCS variations of up to a factor of 12 were observed across the biopolymer-soil mixes studied, in line with theoretically and experimentally expected values, due to the differences in the G:M ratios. The limited impact of molecular weight upon soil strength properties is also shown in CMC-stabilised soils. When considering a soil's stiffness and energy absorbance, the importance of biopolymer-biopolymer interaction and is discussed, further deciphering biopolymer characteristics driving soil property modifications. This study highlights the importance of biopolymer chemistry for biopolymer stabilisation studies, illustrating the use of simple low-cost, accessible chemistry-based instrumental tools and outlining key design principles for the tailoring of biopolymer-soil composites for specific geotechnical applications.

The peak pressure of a granular debris flow at low Froude conditions can be calculated with knowledge of the stress anisotropy and the bulk density as well as the run-up height at impact. Based on a small-scale physical model, measurements of stress anisotropy and flow density values at impact are presented and applied to existing run-up prediction models, and further compared with back-calculated run-up coefficients from measured maximum impact pressures. For this purpose, we conducted 17 experiments with impact measurements and six experiments without impact measurements at Froude numbers, ranging from 0.84 to 2.41. Our results indicate that run-up heights are best reproduced by predictive models, either based on energy or mass and moment conservation, when anisotropic stress conditions, found in this study to range from 1.2 to 5.0, and bulk density variations due to impact, ranging in this study from 0.8 to 2.3, are considered. The influence of stress anisotropy and density variation on the run-up prediction differs, depending on the modelling approach. For the calculation of run-up heights based on the energy conservation concept, the influence of stress anisotropy becomes more significant with increasing Froude number, whereas for models based on mass and momentum conservation, bulk density variations have a greater influence on the estimation of the potential run-up.

This paper introduces a new hypoplastic model characterized by a simple and elegant formulation. It requires only 7 material parameters to depict salient mechanical behaviors of granular materials. The numerical implementation employs an explicit integration method, enhanced by a best-fit stress correction algorithm in a smoothed particle hydrodynamics code. The performance of this model in capturing soil behavior across a range of scenarios is demonstrated by conducting various numerical tests, including triaxial and simple shear at low strain rates, as well as granular collapse, rigid penetration and landslide process at high strain rates.

We propose a constitutive model for both the solid-like and fluid-like behavior of granular materials by decomposing the stress tensor into quasi-static and collisional components. A hypoplastic model is adopted for the solid-like behavior in the quasi-static regime, while the viscous and dilatant behavior in the fluid-like regime is represented by a modified rheology model. This model effectively captures the transition between solid-like and fluid-like flows. Performance and validation of the proposed model are demonstrated through numerical simulations of element tests.