Based on the virtual work increment equation for the three dimensionalcontinum and stress-resultant yield surface plastic rule,the updated Lagrangian formulation is developed,and a new type of geometric stiffness matrix is presented.The numerical examples