Femtosecond laser ablation possesses a variety of applications due to its better control, high power density, smaller heat-affected zone, minimal collateral material damage, lower ablation thresholds, and excellent mechanical properties. The non-Fourier effect in heat conduction becomes significant when the heating time becomes extremely small. In order to analyze the femtosecond laser ablation process, a hyperbolic heat conduction model is established based on the dual-phase-lag model. Taken into account in the model are the effect of heat source, laser heating of the target, the evaporation and phase explosion of the target material, the formation and expansion of the plasma plume, and interaction of the plasma plume with the incoming laser. Temperature-dependent optical and thermophysical properties are also considered in the model due to the fact that the properties of the target will change over a wide range in the femtosecond laser ablation process. The effects of the plasma shielding, the ratio of the two delay times, and laser fluence are discussed and the effectiveness of the model is verified by comparing the simulation results with the experimental results. The results show that the plasma shielding has a great influence on the femtosecond laser ablation process, especially when the laser fluence is high. The ratio between the two delay times (the ratio B) has a great influence on the temperature characteristic and ablation characteristic in the femtosecond laser ablation process. The augment of the ratio Bwill increase the degree of thermal diffusion, which will lower down the surface temperature and accelerate the ablation rate after the ablation has begun. The ablation mechanism of femtosecond laser ablation is dominated by phase explosion. The heat affected zone of femtosecond laser ablation is small, and the heat affected zone is less affected by laser fluence. The comparison between the simulation results and the experimental results in the literature shows that the model based on the dual-phase-lag model can effectively simulate the femtosecond laser ablation process.