Mudstone interlayer is a weak layer in rock engineering. When it is subjected to continuous stress higher than its damage threshold, due to the dislocation of particles in mudstone crystals and the expansion of cracks, mudstone strength is gradually damaged and deteriorated and the strain gradually increases, thus accelerating the phenomenon of creep damage. In order to describe the characteristics of the whole process of mudstone aging deformation, based on the damage evolution of strength parameters (cohesion and internal friction coefficient) with stress and time in mudstone creep tests, a novel damage nonlinear viscoelastoplastic body (D-NVPB) is proposed through improving traditional plastic element. D-NVPB describes the nonlinear characteristics of the accelerated creep stage of mudstone. With the element combination method, D-NVPB is connected with the Burgers model in series to form a new nonlinear damage creep model (D-NVEP model). The analysis results of creep characteristics theoretically verified the rationality of the model in describing the instantaneous elasticity, viscoelasticity, and nonlinear viscoplastic characteristics of the complete creep curve of mudstone. With the data obtained in the uniaxial compression creep test of mudstone under the action of a stress level of 14 MPa, based on the Levenberg-Marquardt nonlinear least squares method, the fitting calculation was performed through piecewise fitting and overall fitting. The correlation coefficient was 0.9909, which verified the applicability of the model. The obtained model parameters by the identification were used to predict the mudstone creep curve under the stress levels of 13 MPa and 15 MPa. The good prediction results further verified the feasibility of the model. Compared with the traditional creep model, the D-NVEP model can better describe the nonlinear characteristics of the accelerated creep stage and quantitatively display the strength damage evolution process of rock in the creep failure process.