To accurately represent crack evolution and hardening-softening in rock elastoplastic constitutive models: (1) Based on conventional triaxial compression and the Mohr-Coulomb strength criterion, the elastoplastic characterization of crack evolution and hardening-softening under plastic strain and confining pressure is examined, along with its relationships with dilatancy angle psi, elastic modulus E, Poisson's ratio mu, cohesion c, and internal friction angle phi; (2) Through experiments on silty mudstone, these mechanical parameters are inverted as bivariate functions of plastic shear strain gamma p and confining pressure sigma 3: psi(gamma p,sigma 3), E(gamma(p),sigma(3)), mu(gamma(p),sigma(3)), c(gamma(p),sigma(3)) and phi(gamma(p),sigma(3)), establishing a model that describes crack evolution and hardening-softening; (3) The predictive capability of the model is validated through secondary development in FLAC(3D) and comparison with experimental results. The study reveals that crack volumetric strain is equal to plastic volumetric strain, and crack evolution can be indirectly quantified through psi(gamma(p),sigma(3)). Disregarding elastoplastic coupling, the elastic phase of the stress-strain curve is represented by E(gamma(p),sigma(3)) and mu(gamma(p),sigma(3)), while the hardening-softening phase is described by c(gamma(p),sigma(3)) and phi(gamma(p),sigma(3)). The simulated results of crack volumetric strain and stress-strain curves closely align with experimental data, indicating that the model adequately characterizes crack evolution and hardening-softening.