With the gradual warming of the global climate

the possibility of large-scale sliding and collapse disasters of glaciers

large ice sheets or thicker ice layers are increasing. Research on the deformation characteristics of polycrystalline ice and establishing its constitutive model are of great significance for the disaster forecasting of polar regions. Based on the existing triaxial compression test data of artificial polycrystalline ice

it is found that the artificial polycrystalline ice has prominent characteristics of strain softening. The damage variable with equivalent plastic strain as the variable is proposed. By using the Mohr-Coulomb criterion as the initial yield function

and adopting the associated flow rule

an elastoplastic damage constitutive model that can reflect the strain softening characteristics of polycrystalline ice is established. The proposed constitutive model is verified by triaxial compression test data under different confining pressures and different temperatures

and parameter sensitivity analysis of some model parameters is carried out. The established elastoplastic damage constitutive model was written in the FLAC software

and the embedded constitutive model was verified by the triaxial compression of a single unit. Then

the embedded elastoplastic damage constitutive model is used to perform triaxial compression on the cylindrical specimen. The simulation results are analyzed in detail. The results provide a theoretical and numerical basis for the stability analysis of glaciers and the multi field coupling study of glaciers.