The stability of geotechnical structures after an earthquake is primarily determined by the residual strength of surrounding soils that have not fully liquefied. This research employs the discrete element method (DEM) to study the undrained post-cyclic shear behaviour of sand under triaxial conditions, focusing on the effect of varying degrees of liquefaction (LD) simulated by subjecting the samples to different lengths of cyclic loading. Different types of cyclic loading, i.e. symmetric (fully reversal), partially reversal, and non-reversal ones, as well as the effect of sample density, have been considered. The results indicate that the samples under fully or partially reversal cyclic loading eventually liquefied, displaying a cyclic mobility failure mode. In contrast, samples under non-reversal cyclic loading develop plastic strain accumulation (PSA) failure without liquefaction. The post-cyclic shear stiffness of the samples is affected by both LD and the type of cyclic loading. For samples under reversal cyclic loading, the post-cyclic shear stiffness decreases as LD increases. Notably, the liquefied samples (LD = 1) initially exhibit near-zero stiffness during post-liquefaction shear until highly anisotropic force chains are formed along the loading direction, with their buckling leading to stiffness recovery. The length of the low-stiffness stage is influenced by the static shear stress and the relative density of the sample, which determines the rate of anisotropy accumulation during cyclic loading. The onset and completion of stiffness recovery are marked by a peak in anisotropy and an abrupt increase in effective anisotropy, respectively. For samples under non-reversal cyclic loading, the post-cyclic shear stiffness initially decreases with the increase in LD but increases at higher LDs due to the significant anisotropy developed during the cyclic loading stage.

Although the mechanical response of granular materials strongly depends on the interplay between their anisotropic internal structure (fabric) and loading direction, such coupling is not explicitly considered in existing high-cycle experimental datasets and models. High-cycle experiments on granular specimens specifically prepared with various fabric orientations are presented. It is found that the high-cycle strain accumulation behavior can change remarkably, from shakedown to ratcheting, when the fabric orientation deviates more from the loading direction. Inspired by the experimental observations, a fabric-dependent anisotropic high-cycle model is proposed, by proper recasting of an existing model formulated within Critical State Theory, into the framework of Anisotropic Critical State Theory. The model explicitly accounts for the fabric evolution, which is linked to plastic modulus, dilatancy and kinematic hardening rules. The model can quantitatively reproduce the high-cycle strain accumulation (i.e., shakedown and ratcheting) under drained conditions, as well as pre-liquefaction and post-liquefaction responses granular materials having widely ranged fabric anisotropy, densities and cyclic loading types using a unified set of constants. It exhibits a unique feature of simulating the distinct high-cycle strain accumulation and liquefaction of granular material with various fabric anisotropy, while the existing high-cycle models treat them equally. The successful reproduction of the anisotropic sand element response under high-cycle drained and undrained conditions makes it possible to perform whole life analysis of various foundations on granular soil subjected to high-cycle loading events.