Conventional finite element analysis of metal forming processes often breaks down due to severe mesh distortion. Since 1993, meshless methods have been considerably developed for structural applications. The main feature of these methods is that the domain of the problem is represented by a set of nodes, and a finite element mesh is unnecessary. This new generation of computational methods reduces time-consuming model generation and refinement effort, and it provides a higher rate of convergence than that of the conventional finite element methods. A meshless method based on the reproducing kernel particle method (RKPM) is applied to metal forming analysis. The displacement shape functions are developed from a reproducing kernel approximation that satisfies consistency conditions. The use of smooth shape functions with large support size are particularly effective in dealing with large material distortion in metal forming analysis. In this work, a collocation formulation is used in the boundary integral of the contact constraint equations formulated by a penalty method. Metal forming examples, such as ring compression test and upsetting, are analyzed to demonstrate the performance of the method.