Using molecular dynamics (MD) methods combining with two-step radiation heating model, the mechanisms of ablation and the thermodynamic states at Ni surface under femtosecond laser irradiation are investigated. Simulation results show that the main mechanisms of ablation are evaporation and tensile stresses generated inside the target. The velocity of stress wave is predicted to be nearly equal to sound velocity. The rates of ablation at different fluences obtained from simulations are in good agreement with experimental data. Superheating phenomenon is also discovered.