In this paper, numerical simulation of nanoparticle integrated laser shock peening of aluminum alloys was carried out. A “tied constraint” was used to connect the matrix and nanoparticle assembly in ABAQUS package. Different particle size and particle volumes fraction (PVF) were studied. It was found that there is significant stress concentration around the nanoparticles. The existence of nanoparticle will influence the stress wave propagation and thus the final stress and strain state of the material after LSP. In addition, particle size, PVF and particle orientation all influence the strain rate, static residual stress, static plastic strain and energy absorption during the LSP process.