Recent development of the shock analysis on the HDD is briefly reviewed. A flexible multi-body dynamics formulation is developed to simulate the shock response of the HDD. If one component in the HDD is changed, only mode shapes and frequencies of that component should be re-calculated and then used to obtain the system’s response. Steady state Reynolds equation is solved to obtain the air pressure on the slider and disk for various slider positions. An air pressure table is formed and used to model the non-linear air bearing during the simulation. Responses of flying height for different direction and shock duration time are analyzed. Results show that the flying state of the slider is more sensitive to the shock with shorter duration time.