Abstract: The large deformation elashc-plastic shell finite element code is developed. In order to improve the convergence of incremental calculations, the C¹-conhninty geometric nonlinear-material nonlinear triangle shallow shell finite elements are applied in the code. Also, a mixed send-inverse finite element contact algorithm is created for effective calculations of contact problems widely existing in sheet metal forming processes. To verify the computational code based on present shell element and contact algorithm, some numerical results are worked out. These results are compared with other numerical or analytical solutions, which turns out to be good in agreement.