This paper presents a rigorous, semi-analytical solution for the drained cylindrical cavity expansion in transversely isotropic sand. The constitutive model used for the sand is the SANISAND-F model, which is developed within the anisotropic critical state theory framework that can account for the essential fabric anisotropy of soils. By introducing an auxiliary variable, the governing equations of the cylindrical expansion problem are transformed into a system of ten first-order ordinary differential equations. Three of these correspond to the stress components, three are associated with the kinematic hardening tensor, three describe the fabric tensor, and the last one represents the specific volume. The solution is validated through comparison with finite element analysis, using Toyoura sand as the reference material. Parametric analyses and discussion on the impact of initial void ratio, initial mean stress level, at-rest earth pressure coefficient and initial fabric anisotropy intensity are presented. The results demonstrate that the fabric anisotropy of sand significantly influences the distribution of stress components and void ratio around the cavity. When fabric anisotropy is considered, the solution predicts lower values of radial, circumferential and vertical stresses near the cavity wall compared to those obtained without considering fabric anisotropy. The proposed solution is expected to enhance the accuracy of cavity expansion predictions in sand, which will have significant practical applications, including interpreting pressuremeter tests, predicting effects of driven pile installation, and improving the understanding of sand mechanics under complex loading scenarios.

In geotechnical engineering, the development of efficient and accurate constitutive models for granular soils is crucial. The micromechanical models have gained much attention for their capacity to account for particle-scale interactions and fabric anisotropy, while requiring far less computational resources compared to discrete element method. Various micromechanical models have been proposed in the literature, but none of them have been conclusively shown to agree with the critical state theory given theoretical proof, despite the authors described that their models approximately reach the critical state. This paper modifies the previous CHY micromechanical model that is compatible with the critical state theory based on the assumption that the microscopic force-dilatancy relationship should align with the macroscopic stress-dilatancy relationship. Moreover, under the framework of the CHY model, the fabric anisotropy can be easily considered and the anisotropic critical state can be achieved with the introduction of the fabric evolution law. The model is calibrated using drained and undrained triaxial experiments and the results show that the model reliably replicates the mechanical behaviors of granular materials under both drained and undrained conditions. The compatibility of the model with the critical state theory is verified at both macroscopic and microscopic scales.

The application of prefabricated assembly technology in underground structures has increasingly garnered attention due to its potential for urban low-carbon development. However, given the vulnerability of such structures subjected to unexpected seismic events, a resilient prefabricated underground structure is deemed preferable for mitigating seismic responses and facilitating rapid recovery. This study proposes a resilient slip-friction connection-enhanced self-centering column (RSFC-SCC) for prefabricated underground structures to promote the multi-level self-centering benefits against multi-intensity earthquakes. The RSFC-SCC is composed of an SCC with two sub-columns and a series of multi-arranged replaceable RSFCs, intended to substitute the fragile central column. The mechanical model and practical manufacturing approach are elucidated, emphasizing its potential multi-level self-centering benefits and working mechanism. Given the established simulation model of RSFC-SCC-equipped prefabricated underground structures, the seismic response characteristics and mitigation capacity are investigated for a typical underground structure, involving robustness against various earthquakes. A multi-level self-centering capacity-oriented design with suggested parameter selection criteria is proposed for the RSFC-SCC to ensure that prefabricated underground structures achieve the desired vibration mitigation performance. The results show that the SCC with multi-arranged replaceable RSFCs exhibits a significant vibration isolating effect and enhanced self-centering capacity for the entire prefabricated underground structure. Benefiting from the multi-level self-centering process, the RSFC-SCC illustrates a robust capacity that adapts to varying intensities of earthquakes. The multi-level self-centering capacity-oriented design effectively facilitates the target seismic response control for the prefabricated underground structures. The energy dissipation burden and residual deformation of primary structures are mitigated within the target performance framework. Given the replacement ease of RSFCs and SCC, a rapid recovery of the prefabricated underground structure after an earthquake is ensured.

Liquefaction resistance and post-liquefaction shear deformation are key aspects of the liquefaction behavior for granular soil. In this study, 3D discrete element method (DEM) is used to conduct undrained cyclic triaxial numerical tests on specimens with diverse initial fabrics and loading history to associate liquefaction resistance and post-liquefaction shear deformation with the fabric of granular material. The influence of several fabric features on liquefaction resistance is first analyzed, including the void ratio, particle orientation fabric anisotropy, contact normal fabric anisotropy, coordination number, and redundancy index. The results indicate that although the void ratio and anisotropy strongly influence liquefaction resistance, the initial coordination number or redundancy index can uniquely determine liquefaction resistance. Regarding post-liquefaction shear deformation, the above quantities do not dictate the shear strain induced after initial liquefaction. Instead, the mean neighboring particle distance (MNPD), a fabric measure previously introduced in 2D and extended to 3D in this study, is the governing factor for post-liquefaction shear. Most importantly, a unique relationship between the initial MNPD and ultimate saturated post-liquefaction shear strain is identified, providing a measurable state parameter for predicting the post-liquefaction shear of sand.

Anisotropic soils exhibit complex mechanical behaviours under various loadsing conditions, e.g., reversible dilatancy, three-dimensional failure strength, fabric anisotropy, small-strain stiffness, cyclic mobility, making it difficult to comprehensively capture these characteristics within a single constitutive model. Failure to capture anisotropic soil behavious may result in poor predictions in geotechnical engineering. Hence, to provide a unified prediction for the mechanical responses of anisotropic sand and clay under both monotonic and cyclic loading conditions, a fabric-based anisotropic constitutive model, i.e., the CASM-CF, is developed within the framework of the standard Clay and Sand Model (CASM) in this paper. Effects of small-strain stiffness and anisotropic elasticity are incorporated into the stiffness matrix to capture the stiffness variation over a wide strain range and reversible dilation. The fabric tensor defined by particle orientation and its evolution law are integrated into the CASM-CF model through the Anisotropic Transformed Stress (ATS) method. The plastic modulus is modified by considering cyclic loading history and stress reverse to better predict the mechanical responses of soils when subjected to cyclic loadings. The newly proposed model is employed to predict the mechanical behaviours of clay and sand under various strain scales and stress paths, including monotonic, cyclic, proportional, and non-proportional loading conditions, in the literature. Conclusions can be drawn that the model performs satisfactorily under various stress paths, and it has the potential to be used in the analysis of practical geotechnical applications of wide range.

In practical engineering, earthquake-induced liquefaction can occur more than once in sandy soils. The existence of low-permeable soil layers, such as clay and silty layers in situ, may hinder the dissipation of excess pore pressure within sand (or reconsolidation) after the occurrence of liquefaction due to the mainshock and therefore weaken the reliquefaction resistance of sand under an aftershock. To gain more mesomechanical insights into the reduced reliquefaction resistance of the reconsolidated sand under aftershock, a series of discrete element simulations of undrained cyclic simple shear tests were carried out on granular specimens with different degrees of reconsolidation. During both the first (mainshock) and second (aftershock) cyclic shearing processes, the evolution of the load-bearing structure of the granular specimens was quantified through a contact-normal-based fabric tensor. The interplay between mesoscopic structure evolutions and external loadings can well explain the decrease in reliquefaction resistance during an aftershock.

The stress state and density of soil have been considered as the key factors to determine the liquefaction resistance. However, the results of seismic liquefaction case histories, laboratory tests and centrifuge model tests show that the fabric characteristics also influence liquefaction resistance, even more significantly than the contributions of stress state and density. In this study, anisotropic specimens with different consolidation histories were prepared using the 3D Discrete Element Method (DEM) to investigate the influence of fabric characteristics on the mechanical behavior of granular materials and the underlying mechanisms. The simulations revealed that under monotonic shear conditions, horizontally anisotropic specimens exhibited strain hardening and dilatancy characteristics, as well as higher peak strength. Under cyclic shear condition, the normalized liquefaction resistance of the specimens showed a strong linear relationship with the degree of anisotropy, independent of confining pressures and density. Microscopic results indicate that the fabric arrangement aligned with the loading direction leads to the evolution of the mechanical coordination number and average contact force in a manner favorable to resisting loads, which is the underlying mechanism influencing macroscopic mechanical properties. Additionally, the evolution patterns of contact normal magnitude and angle in anisotropic granular materials under cyclic loading conditions were also analyzed. The results of this study provided a new perspective on the macroscopic mechanical properties and the evolution of the microstructure of granular soils under anisotropic conditions.

Macro- and micromechanical interactions between the geogrid and granular aggregates considering particle shape effects are essential for the performance of reinforced soil structures under cyclic normal loading (CNL). Crushed limestone and spherical granular media were mixed to obtain samples with different overall regularities (OR = 0.707, 0.774, 0.841, 0.908, and 0.975). Direct shear tests under CNL were conducted at various overall regularities, normal loading frequencies, and waveforms. Consistent with experiment tests, a discrete-element method (DEM) simulation was performed, incorporating authentic particle shapes obtained through three-dimensional (3D) scanning technology. The results showed that the macroscopic interface shear strength and volume change decreased with an increase in the overall regularity and normal loading frequency. The interface shear strength and deformation under the square waveform are bound to be higher than that under other waveforms. The coordination number, porosity, and fabric anisotropy were used to explain the macroscopic interface shear behavior in relation to the overall regularity. A higher coordination number and stronger contact force were observed with a decrease in the overall regularity. As the overall regularity decreased, the interface integrity and stability became stronger, with the result that the reinforced soil structure can withstand a larger principal stress deflection. Through experimental and DEM analyses, the underlying explanation for the effect of particle shape on the mechanical interaction of reinforced soil was revealed.

A series of finite element analyses, conducted on the basis of modified triaxial tests incorporating radial drainage, were carried out to investigate the lateral deformation and stress state characteristics of prefabricated vertical drain (PVD) unit cells under vacuum preloading. The analyses revealed that the inward horizontal strain of the unit cell increases approximately linearly with the vacuum pressure (Pv) but decreases non-linearly with an increase in the initial vertical effective stress (sigma ' v0). The variations in the effective stress ratio, corresponding to the median excess pore water pressure during vacuum preloading of the PVD unit cell, were elucidated in relation to the Pv and sigma ' v0 using the simulation data. Relationships were established between the normalized horizontal strain and normalized effective stress ratio, as well as between the normalized stress ratio and a composite index parameter that quantitatively captures the effects of vacuum pressure, initial effective stress, and subsoil consolidation characteristics. These relationships facilitate the prediction of lateral deformation in PVD-improved grounds subjected to vacuum preloading, utilizing fundamental preloading conditions and soil properties. Finally, the proposed methodology was applied to analyze two field case histories, and its validity was confirmed by the close correspondence between the predicted and measured lateral deformation.

Laboratory experiments have shown that the proportional shearing of granular materials along arbitrary strain path directions will lead to stress states that converge asymptotically to proportional stress paths with constant stress ratios. The macro- and microscopic characteristics of this asymptotic behaviour, as well as the existence of asymptotic states exhibiting a constant stress ratio and a steady strain-rate direction, have been studied using the discrete element method (DEM). Proportional shearing along a wide range of strain-rate directions and from various initial stress/density states has been conducted. The simulation results suggest that general contractive asymptotic states (except for isotropic states) do exist but may be practically unattainable. Dilative strain path simulations, on the other hand, result in continuously changing stress ratios until static liquefaction occurs, indicating the absence of dilative asymptotic states. Despite this difference, a unique relationship between the stress increments and the current stress ratio gradually emerges from all strain path simulations, regardless of strain path direction and initial stress/density conditions. At the particle scale, the granular assembly sheared along proportional strain paths exhibits a constant partition ratio between strong and weak contacts. Although general proportional strain paths are associated with changing geometric and mechanical anisotropies, the rates of change in these anisotropies for contractive strain paths are synchronised to maintain a constant ratio of their contributions to the mobilised shear strength of the material, with a higher proportion being contributed by geometric anisotropy for more dilative strain paths.