A method is introduced for efficient reliability-based design of laser peening (LP) surface treatment to extend fatigue life of metal components. The method includes nonparametric probability density estimation, surrogate modeling using a new finite element (FE or FEA) approach, and reliability analysis with correlated random variables (RVs). Efficient LP simulation is achieved via a new technique termed single explicit analysis using time-dependent damping (SEATD), which reduces simulation times by a factor of 6. The example study of a three-point bend coupon reveals that fatigue life reliability significantly affects optimal LP design, as 52 laser spots are needed for 99% reliability versus 44 spots for 95%.