Contact force of Multi-body dynamics (MBD) system can be classified two parts. First is a normal force and the other is a tangential force called friction force. And the friction force can be represented by two states such as stick and slip. The stick-slip phenomenon is simply described as a simple contact model which is a rigid body contacted on a sloped surface. If the calculated friction coefficient between the body and sloped surface is less than the static friction coefficient, the body should be stuck. If the calculated friction coefficient is greater than the static friction coefficient, the body will be sliding along the surface. The phenomenon is called as stick and slip state of friction, respectively. Usually many researchers and commercial MBD software used a coulomb friction force model which is defined with an only function of relative velocity. This kind of friction force model will be called a conventional friction force model in this paper. A big problem of the conventional model can not describe a stick state of friction phenomenon. In the case of conventional friction force model, the body will be sliding even though friction state is stick. Because, the relative velocity must have a non-zero value in order to generate the friction force. To solve this kind of problem, we propose a stick-slip friction force model including a spring like force. In the case of stick-slip friction force model, the body can be stuck on the sloped surface because the friction force will be a non-zero value, even though the relative velocity approaches zero. We defined a relative displacement variable called stiction deformation. In this paper, the stick-slip friction model is proposed and applied in the contact algorithm of MBD system. And then two friction models are compared with numerical examples. With the proposed stick-slip friction model, more realistic results are achieved.