The paper proposes a multi-body dynamic simulation to numerically evaluate the generated axial force(G.A.F) and plunging resistant force(P.R.F) practically related to the shudder and idling vibration of an automobile. A numerical analysis of two plunging types of CV joints, tripod joint(TJ) and very low axial tripod joint(VTJ), is conducted using the commercial program DAFUL. User-defined subroutines of a friction model illustrating the contacted parts of the outboard and inboard joint are subsequently developed to overcome the numerical instability and improve the solution performance. The Coulomb friction effect is applied to describe the contact models of the lubricated parts in the rolling and sliding mechanisms. The numerical results, in accordance with the joint articulation angle variation, are validated with experimentation. The offset between spider and tulip housing is demonstrated to be the critical role in producing the 3rd order component of the axial force that potentially causes the noise and vibration in vehicle. The VTJ shows an excellent behavior for the shudder when compared with TJ. In addition, a flexible nonlinear contact analysis coupled with rigid multi-body dynamics is also performed to show the dynamic strength characteristics of the rollers, housing, and spider.