Compression Deformation Characteristics of Frozen Soil Containing Ice Lenses Under an Asymmetric Temperature Field

Frozen soil on alpine slopes is influenced by inclination and aspect-induced differential solar radiation effects, resulting in non-uniform temperature fields and inclined layered ice lenses that enhance anisotropy and degrade mechanical properties. This study investigates the deformation and strength responses of frozen soil under different temperature-gradient magnitudes and orientations and ice-lens conditions. A stress–strain constitutive model incorporating the magnitude and orientation of the temperature gradient is established. In addition, an equal-scale discrete element model based on the parallel-bond contact model is developed and calibrated against the laboratory results. The numerical specimen is divided into 13 layers, and temperature-dependent interparticle bond properties are assigned layer by layer to reproduce the prescribed magnitude and orientation of the temperature gradient. Results show that the orientation of the temperature gradient significantly alters the mechanical response and failure mode. As the inclination angle increases, the failure mode transitions from compressive dilatancy to combined dilatancy–shear failure and ultimately to shear-dominated failure. At −10 °C, increasing the inclination angle from 0° to 30° reduces the compressive strength by 44.48%. The elastic modulus also decreases with increasing inclination, with a maximum inclination-induced difference of 111.98 kPa. Moreover, the presence of an ice lens further reduces specimen stiffness, and the elastic-modulus difference between ice-lens-bearing and ice-lens-free specimens increases from 5.57 kPa at −1 °C to 75.72 kPa at −10 °C. The DEM results show that particles at the top and bottom of the specimen primarily undergo vertical displacement, whereas particles in the middle region exhibit dominant horizontal displacement, forming an X-shaped shear band. The inclined temperature gradient produces a heterogeneous distribution of interparticle bond strength within each horizontal layer. As inclination increases, the shear band evolves from symmetric to asymmetric; particle displacements on the side toward which the temperature gradient points are larger than those on the opposite side, revealing the microscopic origins of macroscopic mechanical behavior.