Abstract:Based on updated Lagrangian formulation of elasto-plastic large deformation theory,the numerical simulation method of Biot′s consolidation theory on elasto-plastic large deformation is developed using the idea of parametric variation principle. The advantages of this method are to convert elasto-plastic analysis into finding out the minimum of elasto-plastic potential energy with constrains of yield condition and avoid elasto-plastic iteration at every load step. Furthermore,this method can solve the problems of non-associate flow and stress softening of elasto-plasticity.