I tried both of reducing timestep and mesh size, the negative strain energy is still there, as follows (only those areas with negative strain energy is plotted),n

nThis is a blasting simulation with preload (preload is obtained by implicit solution), thus, it seems that it could not be simulated with SPH.nThanksnn