In this work the formulation of a finite element axisymetric thin wall shell is extended to the numerical analisys including material nonlinear effects such as elasto-plastic phenomena. The algorithm employs the incremental solution scheme in an explicit integration with a load sub-incrementation in the procedure for equilibrium equation solution. This solution integration scheme allows restraining of stress state field with the equivalent stress withn or on the yielding surface. The modified Newton- Raphson method is used in the solution of resulting system of equilibrium equation. The kinematics of deformation employed in the model uses only translation degrees-of-freedom refered to the generatrix of the element midsurface to represent geometric compatibility conditions of the shell. Boudary conditions such as normal rotation fixity and coontinuity are imposed to the model through the penalty method. Some numerical analyses are presented to evalute the model capabilities in nonlinear elastoplastic examples, which results are compared to other numerical analytic solutions available in the literature. Solution convergebce, within controlling of the error allowed in the numerical solution, in discussed and evaluated.