In this paper, we show the shock formation of the solutions to the 3-dimensional (3D) compressible isentropic and irrotational Euler equations with damping for the initial short pulse data, which was first introduced by Christodoulou [9]. Due to the damping effect, the largeness of the initial data is necessary for the shock formation, and we will work on the class of large data (in the energy sense). Similar to the undamped case, the formation of shocks is characterized by the collapse of the characteristic hypersurfaces and the vanishing of the inverse foliation density function \(\mu \) , at which the first derivatives of the velocity and the density blow up. However, the damping effect changes the asymptotic behavior of the inverse foliation density function \(\mu \) and then shifts the time of shock formation compared with the undamped case. The methods in the paper can also be extended to a class of 3D quasilinear wave equations for the short pulse initial data.
