The total Lagrangian formulation and implementation of radial point interpolation method (RPIM) was presented for the simulation of 2D bulk forming processes. In the present method, the displacement field function was approximated by using RPIM shape functions possessing Kronecker delta function property, which can facilitates the enforcement of essential boundary conditions. A simple contact searching algorithm was used to compute the contact force between the tools and the workpiece in the metal forming analysis and to solve the problems of the penalty method and Lagrangian multiplier method. In order to remedy the volumetric locking, the internal nodal forces of the workpiece were split into deviatoric part and volumetric part. The geometrical nonlinearities and material nonlinearities for the dynamic solution were treated by using the total Lagrangian description and an elasticplastic constitutive law. In order to remedy the volumetric locking, the workpiece was deformed into two parts, deviatoric part and volumetric part. The numerical examples show that RPIM is more efficient in dealing with large deformation problems such as bulk forming with severe distortion.