Differential evolution algorithm is a new global optimization algorithm. DE does not require an initial value, and it has rapid convergence, strong adaptability to a nonlinear function, the features of parallelcalculation, especially in adoption to the complex problem of multivariable optimization. The constitutive integration algorithm affecting the incremental calculation step, and convergence and accuracy of the results is a key of finite element analysis. It is usually divided into an explicit and implicit integration. Return mapping algorithm is an implicit integration to avoid solving the equivalent plastic strain directly so that we achieve a fast and accurate solution for the constitutive equations. Making use of DE and return mapping algorithm to program, the elasticplastic finite element simulation and parameter inversion, the inversion and simulation results are verificated, the results show that it is closed to the actual situation, indicating usefulness and correctness of the program.