In this paper a rigid plastic mesh free method for slightly compressible materials is presented for the simulation of bulk metal forming processes. In this model MLS shape functions are used for the approximation of velocity field. These shape functions are defined in initial coordinates of particles. By using this method the support of shape functions covers the same set of particles during material deformation. Transformation method is used for imposition of essential boundary conditions. Transformation matrix is formed only once at the initial stage, because shape functions are defined in initial coordinates system. For evaluation of the presented model, disk upsetting and ring upsetting processes are simulated. Results obtained from this model are compared with experimental data and FEM results and good agreement is observed between them.