  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
D.Christodoulou\cite{christodoulou2007}. 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 energy sense). Similar to the undamped
case, the formation of shock 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.
