A physics-based model has been developed for laser shock peening (LSP) with femtosecond (fs) laser pulses (fs-LSP), which has never been reported in literature to the authors' best knowledge. The model is tested by comparing simulations with measured plume/shock wave front transient propagations and the LSP-induced hardness enhancement layer thickness. Reasonably good agreements have been obtained. The model shows that fs-LSP can produce much higher pressure than LSP with nanosecond (ns) laser pulses (ns-LSP), and it can also generate very large compressive residual stress in the workpiece near-surface layer with a thickness up to similar to 100 mu m. The developed model provides a powerful guiding tool for the fundamental study and the practical applications of fs-LSP. This study, together with the recently reported work by Nakano et al. [Journal of Laser Micro/Nanoengineering 4(1) (2009) 35-38], has confirmed the feasibility of fs-LSP on both theoretical and experimental sides. (C) 2010 Elsevier B. V. All rights reserved.