The present paper discusses the accuracy of finite-element solutions for metals possessing dominant plasticity, resulting in an incompressible response in a large deformation field. It is known that poor numerical solutions are obtained in constrained problems due to incompressibility of deformed metals, but can be improved by selecting element shape function and numerical integration technique, as well as applying the mixed method derived from Lagrangian multipilers. Many studies have been made for rigid-plastic finite-element solutions so far, but largely deformed elastoplastic structural analysis is rarely discussed in the literature. In the present study, we discussed the advantages of such techniques in large deformation analysis using Jaumann stress rate and isotropic hardening hypoelasticity model.