AbstractA computational method, based on an advanced elasto‐plastic large strain formulation, well suited for the analysis of the cone penetration problem, is presented.A new approach of finite strain elasto‐plastic analysis is employed.13The basic (non‐rate) constitutive relations are developed in a spatial reference space to preserve their physical significance. They are subsequently transformed in Lagrangian co‐ordinates, and through simple time differentiation, their rate equations are introduced. The method is computationally implemented with the finite element method and special provisions are taken to allow for the moving boundary conditions of the