Soil liquefaction poses a significant risk to both human lives and property security. Recent in-situ cases have shown that clayey sand can experience multiple liquefaction events during mainshock-aftershock sequences, known as repeated liquefaction. While existing studies have focused on the cyclic behavior of initial liquefaction events, there is a lack of research on the mechanisms and cyclic response of repeated liquefaction in clayey sand. The factors that control repeated liquefaction in clayey sand are still not fully understood. In this study, a series of cyclic triaxial tests were conducted on sand with varying clay content (0 %, 5 %, 10 %, 15 %, and 20 %) under earthquake sequences. The test results showed that the liquefaction resistance initially decreased significantly and then increased with the number of liquefaction events. Sands with higher clay content exhibited earlier recovery of resistance during continuous liquefaction events. The analysis of the test results revealed that the repeated liquefaction resistance of clayey sand was quite intricate. Sands with a relative density (after reconsolidation) below 80 % were primarily influenced by the degree of stress-induced anisotropy, while sands with a relative density above 80 % were mainly affected by relative density.

Triaxial tests are performed for remolded, artificially isotropic, and anisotropic structured samples under undrained conditions at confining pressures of 25, 100, and 200 kPa. Based on these test results, a binary-medium constitutive model is formulated based on homogenization theory and a breakage mechanism to describe the behaviors of structured soils. In this model, the binary-medium material is idealized as a representative volume element (RVE) composed of bonded elements, whose mechanical behaviors are expressed by the linearly elastic model, and frictional elements, whose mechanical behaviors are described by the double-yield surfaces constitutive model. The parameters of the bonded and frictional elements are determined from the test results of structured and remolded samples, respectively. The expressions for the breakage ratio and local stress coefficient matrix are introduced, and their parameters are provided. The computed results are compared with the test results, demonstrating that the model can reflect the main deformation features of structured soil relatively well, including the influence of anisotropy, gradual damage to particle bonding, and pore development.