A simulation program for vehicle dynamic analysis was developed. The Cartesisn coordinate system was used for translational motion and the Euler angle system was used for rotational motion. A three dimensional multi-wheeled vehicle model and equations of motion were derived. Also static equilibrium analysis was added for initial vehicle condition setting. The program user can describe the exact characteristics of suspension spring force and damping force in the user subroutine. A wheel-ground contact model which represents geometrical effect was developed. Two cases of simulation for 16 D.O.F. vehicle model were conducted to validate the developed program by comparing the simulation results with the experimental data.