The propagation of shock waves in an isotropic viscoelastic homogeneous and inhomogeneous medium was considered in    . The nature of the damping of waves in such media is determined by the geometry of the discontinuity surface and the rheological model described by the corresponding relaxation nuclei.
It is of interest to consider shock waves in a fluid-saturated inhomogeneous two-component medium when one of the components is viscoelastic. Porosity    is understood to be effective porosity, taking into account only interconnected pores. Isolated pores are considered as elements of the viscoelastic part of the porous skeleton. It is assumed that the pore sizes are small in comparison with the distance at which the kinematic and geometric characteristics of the motion change. This allows us to assume that the viscoelastic and liquid components of the medium are continuous and at each point of space in this case there will be two displacement vectors, the displacement vector of the viscoelastic component (the skeleton of the porous medium) and the displacement vector of the liquid. We shall assume that the fluid is compressible.
Differential equations are obtained for determining the amplitude of longitudinal and transverse shock waves in a fluid-saturated inhomogeneous two-component viscoelastic medium using the mathematical theory of discontinuities   .
2. Formulation Problem
The interpenetrating motion of a viscoelastic component and a liquid is considered as a fluid motion in a deformable porous medium.
Rheological relations in such a medium through deformations will be written as
Here , are linear integral operators whose relaxation kernels depend continuously on the spatial coordinates:
here , , , ―continuous functions of coordinates, porosity, ―modulus of fluid compressibility, ―of the stress tensor of a porous medium, N―force acting on the liquid per unit area of the cross section of the porous medium, , ―respectively, the components of the strain tensors of components, δ―the Kronecker symbol. Index 1 above after the letter refers to the viscoelastic component, 2 to the liquid.
Relations (1), (2) together with the equations of motion
and Cauchy formulas
represent a closed system for describing the dynamic behavior of a fluid-saturated inhomogeneous two-component medium.
In formulas (3) and (4), ―the effective density of the viscoelastic component, the ―effective density of the fluid, ―the coefficient of dynamic coupling of the viscoelastic component and the liquid, , ―and the density of the viscoelastic component and liquid, respectively.
The repeated Latin indices assume a summation of one to three, in Greek―from one to two. The dot above the letter indicates the time derivative, and the comma below the letter is the derivative of the coordinate.
3. Velocities of Shock Waves
A shock wave in an inhomogeneous two-component medium saturated with a liquid is determined by an isolated surface on which the displacements are continuous in time and coordinates, and the stresses for ces N, and velocities of phase discontinuities undergo. The parameters of the porous medium and their gradients are continuous.
Dynamic relationships must be satisfied on the wave surface of a two-component medium 
where ―the normal velocity of propagation of the surface, the ―components of the unit normal to the surface, directed into the unperturbed region, which depend on the spatial coordinates.
The sign [ ] denotes the difference in the values of a certain quantity on different sides of the rupture surface.
From the rheological relations (1) and (2) written in the discontinuities, dynamic (5) and kinematic conditions of first-order compatibility for the phases, we obtain the system of equations
Assuming that , on the surface where the values of the discontinuities of the
displacements of the components, multiply (6) by and sum over the repeating index i, we obtain a homogeneous system of equations for , . For system (6) to have a nonzero solution, its determinant must be zero. This condition leads to an equation for determining the propagation velocity of a longitudinal wave
Equation (7), generally speaking, has two distinct roots, each of which determines the propagation velocity of the longitudinal wave and
Thus, in an inhomogeneous viscoelastic two-component medium, there are two types of longitudinal waves, which are determined by Formula (I).
On the other hand, if on the surface , , then from (6) it follows that the propagation velocity of the transverse wave is given by
Thus, it is shown that in a two-component inhomogeneous viscoelastic medium, there are three types of shock waves: two types of longitudinal waves and one transverse, for which, , respectively, and the velocities of these waves are determined by Formulas (7) and (8).
If the connection between the viscoelastic component and the liquid is weak , then from (7) and (8)
, , (9)
It follows from (9) that the propagation velocities of waves in a viscoelastic two-component medium are equal to the velocities of waves propagating separately in a continuous viscoelastic medium and liquid.
4. Determination of Shock Wave Amplitudes
Let us determine the change in the amplitude of the waves. To do this, we differentiate (1) with respect to u and take their difference on different sides of the discontinuity surface, we obtain
Taking into account the compatibility conditions of the first order  and the equations of motion written in the discontinuities, (10) take the following form
The quantities, , , χ and are determined on the surface and characterize the jumps of the first derivative of the stresses, the forces acting on the fluid and the velocity of displacement of the components, respectively, ―the coefficients of the first fundamental quadratic form of the surface, ―the derivatives of the Cartesian coordinates along the curvilinear coordinates and ―the time derivative.
By the standard method, we exclude the quantities , χ from (11). As a result, we get
From the rheological relations (1) written in the discontinuities, the dynamic (5) and kinematic conditions of first-order compatibility, we find, ,
, , , and substitute in (12)-(13). After the transformations,
we obtain a system of equations determining the change in the components of
the wave amplitude vector
where Ω―the average curvature of the wave surface.
From the relations (14) and (15) we obtain differential equations for determining the change in the amplitude of the longitudinal and transverse waves ( , , , ).
For longitudinal waves, we multiply (14) and (15) by and sum over the repeated indices. As a result, we get
We eliminate the quantity from (16) and (17). To do this, we multiply (16) by, and (17) by and subtract, we obtain the following equation of the first component of the two-component medium
When Equation (18) was obtained, the first Equation (6) and Formula (7) were used.
For the transverse wave from (14) and (15), taking into account the second Equation (6) and Formula (8), we similarly obtain
Passing to the variable denoting the distance along the normals to the surface in (18) and (19), we obtain equations for changing the amplitude of the longitudinal and transverse waves during their propagation in the first component
, (20) (III)
where ―for longitudinal waves; for a transverse wave,
―for longitudinal waves;
―for a transverse wave.
The change in the amplitude of the waves in the second component is obtained from (6)
, , (21) (IV)
Then the change in the amplitude of the waves during their propagation in a porous medium will be written in the form
Equations (20) contain mean curvature , which is a function of σ.
The average curvature is related to the first and second quadratic forms and the Gaussian curvature of the wave surface by the equations  , where its method of finding is given.
The amplitude level satisfying Equation (20) is found by the method of successive approximations under the initial conditions, , ―approximation. A homogeneous medium corresponds to the zeroth approximation. In the first approximation, the rate of change in the inhomogeneity of the viscoelastic two-component medium along the beam is taken into account, and in the second, across the beam.
A more detailed solution of equations of the form (20) by the method of successive approximations was considered in   .
We consider an inhomogeneous viscoelastic two-component medium saturated with a liquid, characterized by parameters and relaxation nuclei (a special case)
At the moment of time in the plane xy along the axis OX, the wave front propagates with velocity and , which is determined from (7) (8) with .
Since by the assumption of the problem the wave surfaces form a family of parallel planes, then the average curvature at all times.
Then from (20) for, we obtain the dependence of the level of the wave amplitude on the velocity, the physico-mechanical characteristics and the relaxation kernels of the viscoelastic two-component medium for the first phase
where the value of the function for , , и and are from (18).
By the formula (22) we obtain
Expression (25) determines the change in the level of the wave amplitude in a fluid-saturated inhomogeneous viscoelastic two-component medium.
If we specify in (25) a specific form of the physico-mechanical characteristics of the medium and the relaxation nuclei, we obtain an expression for determining the change in the level of the wave velositi in a fluid-saturated inhomogeneous viscoelastic two-component medium.
Applying the mathematical theory of discontinuities to the basic equations of a two-component viscoelastic porous medium, formulas were first obtained for determining the velocities and changing the level of the amplitude of longitudinal and transverse waves.
In the absence of dynamic coupling between the propagation velocity components of the first and second types of longitudinal waves, the wave velocities propagating separately in a continuous viscoelastic component and a continuous liquid component are equal.
In the medium under consideration, in the absence of a connection between the viscoelastic component and the liquid, the transverse wave propagates only in a continuous viscoelastic component.
Using geometric and dynamic conditions of second-order compatibility of the theory of discontinuities, equations are obtained for determining shock wave amplitudes in a saturated inhomogeneous viscoelastic two-component medium.
The amplitude of the waves in an inhomogeneous viscoelastic two-component medium depends on the porosity of the medium, the viscosity, the initial mean curvature, and the Gaussian curvature of the wave surface.
A List of Symbols
: stress tensor of a porous medium
N: force acting on the liquid per unit area of the cross section of the porous medium
: components of the strain tensors of the elastic component
: components of fluid deformations
: displacement of the elastic component
: fluid movement
: linear elastic integral operators
: wave amplitudes
: average curvature of the wave surface
: the coefficients of the first fundamental quadratic form of the surface
 Biot, M.A. (1956) Theory Propagation of Elastic Waves in a Fluid-Saturated Porous Solid I. Low-Frequency Range. The Journal of the Acoustical Society of America, 28, 168-178. https://doi.org/10.1121/1.1908239
 Biot, M.A. (1956) Theory Propagation of Elastic Waves in a Fluid-Saturated Porous Solid. II. Higher Frequency Range. The Journal of the Acoustical Society of America, 28. 179-191. https://doi.org/10.1121/1.1908241