media are given by [@DAVISON1966249; @ZIV2019mom]

\[integralCurves\] $$\begin{aligned}\textrm{quasi-pressure:\ \ } & \derivative{\shearstrain}{\pressurestrain}=\eta_{+},\end{aligned}
\label{eq:integralCurvesPressure}$$ \[etea\] $$\begin{aligned}\textrm{\;\;\;\;\ \;\;\;\,quasi-shear:\ \ } & \derivative{\shearstrain}{\pressurestrain}=-\eta_{-}^{-1},\end{aligned}
\label{eq:integralCurvesShear}$$

with the compatibility condition $$\shearstrain\left(\pressurestrain=\initial{\pressurestrain}\right)=\initial{\shearstrain},\label{eq:comp0}$$ where $\initial{\pressurestrain}$ and $\initial{\shearstrain}$ can be associated with axial and shear pre-strains, respectively. Eq.$\ $ is intimately related to the Poynting effect [@Poynting1909prsa], namely, the axial deformation that accompanies finite shear, hence occurs only at nonlinear deformations [@ogden97book]. Thus, this effect is prominent in soft materials, which are capable of undergoing large strains [@Cioroianu2013PRE; @Horgan2017prsa]. The exact manner in which this effect takes place depends on the constitutive response of the material [@Mihai2017nr]. Accordingly, the solution of Eq.$\ $ also depends on the constitutive response of the material, which in our study is modeled using the compressible model with the strain energy function [@olp2007] $$\hyper\left(\defgrad\right)=-\frac{\mu\jm}{2}\ln\left(1-\frac{\tr{\defgradT\defgrad}-3}{\jm}\right)-\mu\ln\det{\defgrad}+\left(\frac{\bulk}{2}-\frac{\mu}{3}-\frac{\mu}{\jm}\right)\left(\det{\defgrad}-1\right)^{2};\label{eq:Gent}$$ here, $\mu$ and $\kappa$ coincide with the shear and bulk moduli in the limit of small strains, respectively, and the dimensionless parameter $\jm$ models a stiffening of strain. The Gent model was originally developed to capture the stiffening of rubber materials . Recently, the model was also shown suitable for capturing nonlinear wave phenomena such as shear shocks and tensile-induced shocks, which were experimentally observed by @Catheline2003PRL and @espindola2017shear, 