“…In this section, we introduce a phase field formulation to model the behaviors of intact glass under triaxial compression . As experiments showed that stress–strain follows a nonlinear relationship under triaxial compression, a nonlinear form of the elastic energy density f el is employed: where K is the material bulk modulus. The strains e i ( i = 1,…, 6) are defined as e 1 = ɛ 11 + ɛ 22 + ɛ 33 , e 2 = ɛ 11 − ɛ 22 , , e 4 = 2ɛ 23 , e 5 = 2ɛ 13 , e 6 = 2ɛ 12 , and , .…”