The strong dependent behavior of semi-crystalline polymers can lead to the use of simplified material laws in Finite Element structural calculations for reasons of robustness to the detriment of the quantitative response of the models. This work focuses on numerical integration methods as a solution to overcome the possible convergence and robustness limitations of mean-stress dependent elastoviscoplastic material laws, typical of the semi-crystalline polymers mechanical behavior.What is proposed here is a rational application of three explicit integration methods (fourth and second order Rung-Kutta method, a hybrid schema between Runge-Kutta and Euler method) in engineering structural calculations, which provide a reliable solution for constitutive models of semi-crystalline polymer. These methods are examined for structure creep test and tensile test, in comparison with experimental data. The investigations have been done in terms of the stability toward convergence, the accuracy of results, the plastic consistency, and CPU time efficiency. This work, proposes an easy implementation of integration methods in any computational Finite Element code. It also provides a flexible modular implementation which is applicable to any different constitutive equations. An integration step sub-division technique is recommended. It is a powerful technique to improve the convergence of solution and accuracy of result by damping oscillation around stress Gauss point integration solution. The results obtained illustrate the effect of numerical integration schemas on structural analysis and provide an insight to select suitable method.