Frozen soil creep is a key factor in the settlement of cold region subgrades. Clarifying its macro and micromechanical deformation and damage mechanisms is essential for mitigating subgrade distress. To overcome the limitations of the existing creep contact model in PFC for accurately simulating the non-attenuating creep behavior of frozen soil, this study proposes a new discrete element creep contact model that incorporates damage evolution. By introducing a variable-stiffness damage element, the model effectively captures the third-stage creep behavior of frozen soil. Using the central difference method, the study derives the mechanical response equations for each stage of the contact model and compiles a dynamic link library (DLL) in C++ for direct use within the PFC3D software. A comprehensive description of the model development process is provided, and the contact model is validated using a double-sphere model, reproducing the non-attenuating creep curve and revealing the influence of new parameters on creep behavior. As a representative case, a series of triaxial creep tests on frozen subgrade soil under varying temperatures, confining pressures, and deviatoric stresses are conducted, leading to the establishment and calibration of a corresponding discrete element method (DEM) model. The results demonstrate that the model can accurately reproduce the creep mechanical properties of frozen soil across all stages under complex conditions. Furthermore, microscopic analysis of frozen soil reveals the evolution of particle displacement, damage development, internal structure, and pore distribution during different creep stages. The findings extend the application of DEM in studying the mechanical properties of frozen soil and provide a more precise contact model for DEM simulations of creep in frozen soil. It can be further applied to engineering-scale DEM studies of creep in cold region subgrades.
扫码关注我们
求助内容:
应助结果提醒方式:
