## Abstract

The Lame system describing the dynamics of an isotropic elastic medium affected by a steady transport load moving at subsonic, transonic, and supersonic speed is considered. Its fundamental and generalized solutions in a moving frame of reference tied to the transport load are analyzed. Shock waves arising in the medium at supersonic speeds are studied. Conditions on the jump in the stress, displacement rate, and energy across the shock front are obtained using distribution theory. Transport boundary value problem for an elastic medium bounded by a cylindrical surface of arbitrary cross section and subjected to transport loads is considered in the subsonic and supersonic case with regard to shock waves. To solve problems, the generalized functions method is developed. In the space of generalized functions, generalized solutions are constructed and their regular integral presentations are obtained. Singular boundary equations solving the boundary value problems are presented.

### Keywords

- elastic medium
- transport load
- subsonic
- transonic
- supersonic speed
- shock waves
- boundary value problem
- generalized functions method
- generalized solutions
- singular boundary equations

## 1. Introduction

A widespread source of wave generation in continuous media is transport loading, i.e., moving loads whose form does not change over time. The velocity of a transport load has a large effect on the type of differential equations describing the dynamics of the medium. The equations depend parametrically on the Mach numbers, i.e., on the ratio of the speed of motion to the propagation speeds of perturbations in the medium (sound speeds). It is well known [1] that, in an isotropic elastic medium, there are two sound speeds (*c*_{1}*,c*_{2}), which determine the velocities of dilatational and shear waves propagation. This has a large effect on the type of equations and leads to systems of elliptic, hyperbolic, or mixed equations. For transport problems, typical factors are shock effects generated by supersonic loading. At shock fronts, the stresses, displacement rates, and energy density are discontinuous. A convenient research method for such problems is provided by the theory of generalized functions (distributions), which makes it possible to significantly expand the class of processes amenable to study by using singular generalized functions in the simulation of observed phenomena. In this chapter, methods of this theory are used to solve boundary value problems using motion equations of the theory of elasticity in cylindrical domains under the action of transport loads, moving at supersonic and supersonic speeds.

## 2. Motion equation of elastic medium

We consider an isotropic elastic medium with Lame’s parameters *λ,μ,* and a density *ρ*. Let us denote *x* = *xjej, ej* as the unit vectors of Cartesian coordinate system in the space *R*^{3}; displacements vector *u*(x*,t*) = *ujej*; stress tensors *σij* deformation tensor *εij*. These tensors are connected by Hook’s law [1]:

The elastic constant tensor has the symmetry properties.

In the case of an isotropic medium, it is equal to

and Hook’s law has the form

Here

Motion equations for material continuum

for elastic medium by using Eqs. (1) and (2) have the form:

Here *L* is the matrix Lame’s operator:

*c*_{1} > *c*_{2}), *G*(x*,t*) is the mass force, ∆ is the Laplace operator.

The system shown in Eq. (4) was fairly well studied by Petrashen [2]. Since the elastic potential of the medium is positive definite, this system is strictly hyperbolic. Such systems can have solutions with discontinuous derivatives. The discontinuity surface *F* in *R4* = *R*^{3} × *t*(−∞ *< t <* ∞) coincides with a characteristic surface of the system. It corresponds to a wave front *Ft* moving in *R*^{3} at the velocity *V*:

We note that *ν*(*x,t*) = (*ν*_{1}*,ν*_{2}*,ν*_{3}*,νt*) is a normal vector to *F* in *R*^{4}, satisfying the characteristic equation

This equation has the roots:

From Eqs. (5) and (7), we get that *R*^{3} at the sound velocity *V* = *c*_{1} or *V* = *c*_{2}.

We introduce a wave vector *m* = (*m*_{1}*, m*_{2}*, m*_{3}). It is a unit normal vector to *Ft* in *R*^{3} for fixed *t* in the direction of wave propagation. By virtue of Eq. (7),

Let *νt* = *ν*_{4}. The requirement that the displacements be continuous across the wave front, i.e.,

which is associated with the preservation of the continuity of the medium, leads to kinematic consistency conditions for solutions at the wave front:

(the continuity of the tangent derivatives on *Ft*). Additionally, Eq. (4) implies dynamical consistency conditions for solutions at the wave front, which are equivalent to the momentum conservation law in its neighborhood:

*Definition*. A wave is called a *shock wave* if the jump in the stresses across the wave front is finite: *eimj*[*σij*]*Ft* = 0̸. If *mj*[*σij*]*Ft* = 0, then this is a *weak shock* wave. If *mj*[*σij*]*Ft* = ∞, then this is a *strong shock* wave.

Velocity suffers a jump discontinuity across a shock front. At fronts of weak shock waves, the velocities are continuous, but the second derivatives of solutions are not. Strong shock waves (in the sense of the aforementioned definition) do not occur in actual media, since, at large stress jumps, the medium is destroyed and ceases to be elastic. However, strong shock waves in elastic media play an important theoretical role in the construction of solutions, specifically, fundamental solutions of Eq. (4).

## 3. Lame transport equations and Mach numbers

Suppose that the force affecting the medium moves at a constant velocity *c* along the *X*_{3} axis (for convenience, in its negative direction) and, in a moving coordinate system *t*:

Transport solutions are solutions of Eq. (4) with the same structure:

The speed of transport loads is called *subsonic* if *c < c*_{2},*transonic* if *c*_{2} *< c < c*_{1}, and *supersonic* if *c > c*_{1}. A speed is called the *first* or *second sound* speed if *c* = *cj, j* = 1, 2, respectively.

In the new variables, the equations of motion are brought to the form

Here *M*_{1} *< M*_{2}).

As *Mj <* 1(*j* = 1, 2) the load is subsonic and the system of equations is elliptic. If the load is supersonic, i.e., *Mj > l, j* = 1, 2, then the system becomes hyperbolic. In the case of transonic speeds, i.e., *M1 <* 1 and *M*_{2} *>* 1, the equations are hyperbolic-elliptic. In the case of sound speeds, the equations are parabolic-elliptic if *M*_{2} = 1 and parabolic-hyperbolic if *M*_{1} = 1. We will show this later when considering fundamental solutions of Eq. (14).

Since the original system is hyperbolic, Eq. (14) can also have discontinuous solutions. Let *F* be a discontinuity surface in the space of variables *x*^{′} such that it is stationary in this space and moves at one of the sound velocities *V* = *c*_{1}*,c*_{2} in the space of (*x*_{1}*,x*_{2}*,x*_{3}). It follows from Eq. (7) that *V* = *cn*_{3}, where *n* = (*n*_{1}*,n*_{2}*,n*_{3}) is the unit normal to *F* in *R*^{3}. Therefore, since *c* = *cj/n*_{3} and |*n*_{3}| ≤ 1, such surfaces can arise only at supersonic speeds: *c* ≥ *cj*.

It follows from Eqs. (9) to (11) and Eq. (13) that the kinematic and dynamical consistency conditions for solutions at discontinuities in the mobile coordinate system have the form:

*wave vector*, *k* = 1 for shock dilatational waves, *k* = 2 for shock shear waves. Here and hereafter, the derivative with respect to *xj* is denoted by the index *j* after a comma in the function notation or by the variable itself.

*Definition.* If *c > c*_{2,} the solution of the system in Eq. (14) is called *classical* if it is continuous and twice differentiable everywhere, except for, possibly, wave fronts. The number of fronts is finite at any fixed *t* and the conditions on the gaps, Eqs. (15) and (16), are satisfied on the wave fronts.

At first, we construct the solutions of the transport Lame equation using methods of generalized functions theory.

## 4. Shock waves as generalized solutions of transport Lame equations: conditions on wave front

Consider Eq. (14) and its solutions on the space of generalized vector functions *D*’_{3} (*R*^{3}) with components being generalized functions from *D*^{′}(*R*^{3}) (see [3]). Obviously, if *u* is a solution of Eq. (14) that is twice differentiable, then it is also a generalized solution of Eq. (14). If a vector function *u* satisfies Eq. (14) in the classical sense almost everywhere, except for some surfaces, on which its derivatives are discontinuous, then, generally speaking, *u* is not a generalized solution of Eq. (14).

Let *u*(*x,z*) be a shock wave (*x* = (*x*_{1}*,x*_{2})), i.e., a classical solution of the Lame transport equations, Eq. (14), that satisfies conditions Eqs. (15) and (16) at the front *F*. Let

*The shock wave* *is a generalized solution of the Lame equation in D*’_{3} (*R*^{3}).

**Proof.** Using the rules for differentiating generalized functions with derivatives having jump discontinuities across some surfaces (see [3]), for the equations of motion in *D*’_{3} (*R*^{3}), we obtain

Here, the right-hand side involves singular generalized functions, namely, single layers *δF* (*x,z*) and double layers on *F*. By virtue of conditions Eqs. (15) and (16), the densities of these layers are equal to zero, so the right-hand side of Eq. (17) vanishes; i.e., the shock wave satisfies the same equations, Eq. (14), but in the generalized sense.

As a result, we obtain a simple formal method for deriving conditions at jumps in solutions and their derivatives across the shock fronts in hyperbolic equations. Namely, these equations are written in the space of generalized functions and the densities of the singular functions corresponding to single, double, etc., layers are set to zero.

Define as follows the *kinetic energy density*

and *elastic potential*

Consider the following functions: the *energy density E* = *K* + *W* of elastic deformations and the *Lagrangian* Λ = *K* − *W*.

*If G is continuous, then the Lagrangian Λ is continuous at the shock waves fronts.*

(*and the jump in the energy density satisfies the relation*

*First formula is equivalent to the equality:*

where *F*, *F*.

The last formula may be easy to get if we write the equation for *E* in

as

## 5. Fundamental Green’s tensors and generalized solutions of transport Lame equations

The matrix of fundamental solutions

This matrix is called *Green’s tensor* for the transport Lame equations if it satisfies the decay conditions at infinity

For a fixed *k*, its components describe the displacements of the elastic medium under a concentrated force moving at the velocity *c* along the axis *Z* = *X*_{3} and acting in the *Xk* direction.

Green’s tensor can be obtained by taking the Fourier transform of Eq. (17) and solving the corresponding system of linear algebraic equations for the Fourier transforms *ξ*_{1}*,ξ*_{2}*,ξ*_{3}). It is reduced to the form (see [4]).

It can be seen that

are of crucial importance in the construction of the original Green’s tensor. It is easy to see that *f*^{¯}_{0m} is the Fourier transform of the fundamental solution to the equation

This equation is similar to the elliptic Laplace equation at subsonic speeds if *z* disappears from the equation and the equation becomes parabolic, since the space dimension is higher by one, which determines the type of Eq. (14), as noted earlier, since the solutions contain waves of two types. Green’s tensor for the Lame transport equation was constructed by Alekseyeva [4] by applying fundamental solutions of the Laplace and wave equations and regularization functions *f*^{¯}*km*, which depends on the speed of transport load. Green’s tensor has the regular form:

where the type of basic function depends on velocity *c.*

*In subsonic case* (*Mk <* 1):

*In sonic case* (*Mk* = 1):

*In supersonic case* (*Mk >* 1):

Here and hereafter, we use the following notation:

The dilatational and shear components of

In the supersonic case, the support of the functions is the cone *z > mk*∥*x*∥. This determines a radiation condition as physical considerations imply that there are no displacements of the elastic medium outside this cone since the perturbations have a finite propagation velocity, which cannot be higher than the corresponding sound velocity for a particular type of deformation. At the fronts of shock waves (*z* = *mk*∥*x*∥), Green’s tensor grows to infinity.

If the following convolution exists,

it is easy to prove that it is the generalized solution of the transport Lame equations, Eq. (14).

If mass forces are regular, then Eq. (28) has an integral presentation:

If mass forces are concentrated on surface *D* and described by singular generalized functions of the type of single layers

Moreover, by the Du Bois-Reymond lemma [3], these solutions are classical. For other types of singular mass forces, to calculate Eq. (28), we use the definition of convolution of a generalized function [3].

It is easy to see from Eqs. (23) to (25) that the solution is represented as a composition of fundamental solutions distributed over the support of the function *f*(*x,z*); their intensities are determined by its value.

In Alexeyeva and Kayshibayeva’s paper [5], there are some numerical examples of calculation of the dynamic of elastic medium at subsonic, transonic, and supersonic speed of transport loads moving along the strip in an elastic medium.

## 6. Subsonic Green’s tensor, fundamental stress tensors, and their properties

In the subsonic case from Eq. (25), we obtain the components of Green’s tensor in the form:

They are regular functions. Since by *x*^{′} → 0 [6]:

these components are bounded for *R*^{−2} is the order of the tensor derivatives asymptotic and the behavior of at

Tensor

Then the elastic constant tensor is presented in the form

Tensor *n* in a point *x*^{′} = (*x,z*). Tensor

*Fundamental stress tensor* *is the generalized solution of the transport Lame equation with singular mass forces of the multipole type:*

where

*For any closed Lyapunov’s surface D, bounding a domain*

*where* *is the characteristic function of* *, which is equal to 0.5 at D;* *is a unit normal vector to D. The integrals are regular for* *and are taken in value principle sense* for

These formulas have been proved by Alexeyeva [6]. The formula in Eq. (35) can be referred to as a dynamic analog of the well-known Gauss formula for a double-layer potential of the fundamental solution of Laplace’s equation ([3]: 403). It plays a fundamental role in the solution of transport boundary value problems (BVP).

## 7. Statement of subsonic transport boundary value problems. Uniqueness of solution

Let *D*^{−} be an elastic medium bounded by a cylindrical surface *D* with generator parallel to the axis *X*_{3}; let *S*^{−} be the cross-section of the cylindrical domain; let *S* be its boundary, and let *n* be the unit outward normal of *D*. Obviously, *n* = *n*(*x*) and *n*_{3} = 0. We assume that *G* is an integrable vector function and ∃*ε >* 0 such that

There is the subsonic transport load *P*(*x,z*) moving along the boundary *D* (*c < c*_{2}):

We assume that ∃*εi >* 0:

A vector function *u*(*x,z*) satisfying the aforementioned conditions is referred to as a classical solution of the BVP. Let *Cab*^{−} = {(*x,z*): *x* ∈*D*^{−}*, a < z < b*}. The two useful energetic equalities have been proved by Alexeyeva [6].

*Classic solution of transport BVP satisfying to the equalities:*

The following assertion is its corollary.

*The solution of the subsonic transport boundary value problem is unique.*

**Proof.** Since the problem is linear, it suffices to prove the uniqueness of the zero solution. Let *u*(*x,z*) satisfy the zero boundary conditions *P(x,z)* = 0 on *D* and be a solution of the homogeneous Lame equations ((Eq. (14)) by *G*(*x,z*) = 0.

Then for ∀*z*

It follows from the formula (Eq. (40)) of Theorem 7.1. The integrand is a positive quadratic form in *ui,j*, since the elastic potential satisfies the relation *W* ≥ 0 ([1]: 589, 591); moreover, *W* = 0 only for displacements of the medium treated as an absolutely rigid body. Therefore, Eq. (40) is true only if *ui,j* = 0 for all *i, j*. This, together with the decay of solutions at infinity and the arbitrary choice of *z*, implies that *u* = 0.

The proof of the theorem is complete. It is valid both for the internal and external boundary value problem. The asymptotic requirements on *G* and the boundary functions may be weakened.

## 8. General functions method: statement of subsonic transport BVP in D 3 ' R 3

Our aim is to construct the solution of BVP by using boundary integral equations (BIE) for *u*(*x,z*). The construction of an analog of Green’s formula for solutions of elliptic equations ([3]: 366), which permits one to determine the values of the desired function inside the domain on the basis of the boundary values of the function and its normal derivative, is the key point in the construction of BIE of boundary value problems. An analog of this formula for equations of the static theory of elasticity is referred to as the *Somigliana formula* [1]. It determines the function *u*(*x,z*) in the domain *D*^{−}, if the boundary values of displacements *uD*(*x,z*) and stresses *p*(*x,z*) are given. We construct a dynamic analog of that formula in the case of transport solutions. To this end, we use the method of generalized functions (GFM).

We introduce the regular generalized solution of BVP

which defines it as a regular vector function on all space *R*^{3}. Here *D*: 1(*z*) ≡ 1, *S*^{−}, which is equal to 0.5 at *S*: *nj*(*x*)*δS*(*x*) is a simple layer at *S*.

By using the properties of the differentiation of regular generalized functions with jumps on *D*, we obtain the equation for

*D*. Since *n*_{3} = 0 on *D*, it follows from the properties of the Green tensor that an analog of the Somigliana formula holds in the space of generalized functions:

which we write in a form more suitable for transformation as:

If we write out this convolution in integral form with regard to the notation introduced here and Eqs. (1) and (2), then we obtain a formula, whose form coincides with the Somigliana formula for problems of elastostatics ([1]: 605):

where we introduce the shift tensors:

This formula permits one to determine displacements in the medium on the basis of known boundary values of displacements and stresses. But the integrals are regular only *for*

## 9. Singular boundary integral equations of subsonic transport BVP

The following assertion provides a solution for the aforementioned boundary value problems.

*If the solution u(x;z) of subsonic transport BVP satisfies the Holder condition on D; namely*,

then u(x;z) satisfies the singular boundary integral equation

**Proof.** Let consider Eq. (45) for *x*^{∗}*,z*^{∗}) ∈ *D*, *x*^{′} → (*x*^{∗}*,z*^{∗}). Then, using Theorem 6.1, we have

In the last relation, we have used the obvious properties: integrals with *u* on *D* and weak singularity *D.* Then if the surface integral exists, its value coincides with the principal value; the principal value of the integral containing the difference of integrated functions is equal to the difference of the principal values of integrals corresponding to each of these functions if they exist.

By transposing the last two terms to the left-hand side of the relation, we obtain the formula of the theorem for the boundary points. The proof of the theorem is complete.

This theorem gives us resolving system of integral equations for defining unknown values of boundary displacements.

Note also that the subsonic analog of the Somigliana formula was obtained for generalized functions. But since they are regular, from the Dubois-Reymond lemma ([3]: 97), the solution is classical. However, if the acting loads are described by singular generalized functions, which often takes place in physical problems, then one should use a representation of a generalized solution in the convolution form (Eq. (43)) with the evaluation of convolutions by the definition (see [3]: 133).

## 10. Supersonic green’s tensor and its antiderivative with respect to *z*

From Eq. (25), we get the regular representation of

Here

One can readily see that its components are zero outside the sonic cones:

On the surfaces of the cones, the components *z*^{2} − *m*^{2}*jr*^{2})^{−1/2}.

For solution of supersonic problems, we introduce the tensor *z*:

They are also fundamental solutions of Eq. (14) for the mass forces of the corresponding

Tensor *Kj*

it continues on fronts *Kj*. *x’ =* 0 and weak logarithmic singularity on Z with respect to ∥*x*∥ by *x* = 0. To single out these singularities, we decompose it into the terms:

The tensors *z* inside the sonic cones *l* = 1, 2) and have a logarithmic singularity with respect to ∥*x*∥ on the *Z*-axis. Unlike the generating tensor *Kl*. One can readily see that the tensor shifts

have the following symmetry properties around the Z-axis:

But for the components with indices (*i,j*) = (1,3), (3,1), (2,3), (3,2)

## 11. Fundamental supersonic antiderivative stress tensor H ̂ and its properties

We introduce antiderivative stress tensor

This tensor can be obtained in a different way, by analogy with *T,* using Hooke’s law, except that the Green tensor should be replaced with its antiderivative *W*. By using the presentation of the basic functions of Green’s tensor construction (Eq. (25)) in the supersonic case, it can be presented in the following form:

Obviously, for *z < τ,* all the introduced shifted tensors are zero. It has the following symmetry properties around the Z-axis:

except for (*i,j*) = (1,3), (2,3), (3,1):

Components *x*∥^{−1} on the axis Z. If we put Eq. (51) in Hook’s law, then we can again single out two terms in

Since the tensors *z* inside the sonic cones *Kl* (*l* = 1,2), we conventionally say that they are *stationary*. Accordingly, the tensors *dynamic*, because they depend essentially on *z*, although they are regular functions. The aforementioned symmetry properties hold for both stationary and dynamic terms in the tensors.

For this type of tensors, the next theorem was proved (see [7]).

*The fundamental stress tensor H satisfies the relation*

*For x∉D all integrals are regular; for x ∈ D the first integral is singular, calculated in value principle sense.*

This theorem enables us to obtain solvable singular boundary integral equations for a supersonic transport boundary value problem.

## 12. Statement of supersonic transport BVP: uniqueness of solutions

We suppose here that supersonic transport loads, moving at supersonic velocity *c > c*_{1}, are known on the boundary *D*:

Functions *pj*(*x,z*) are integrable on *D*_{+}. We assume here *G* = 0 and

For ∥(*x,z*)∥ → ∞

The jump conditions, Eqs. (15) and (16) are satisfied on the shock wave fronts.

*The solution of the supersonic transport boundary value problem is unique*.

**Proof.** Suppose that there exist two solutions. Since the problem is linear, it follows that their difference *u*(*x,z*) satisfies the zero boundary conditions, i.e., *P*(*x*;*z*) = 0, and is a solution of the homogeneous equations of motion (*G* = 0). We note, that Lemma 8.1 is also true in the supersonic case for shock waves as there is Theorem 3.2 for the gaps of energy on their fronts (see full proof). Then together with conditions given in Eq. (59) of decay of the solutions at infinity and the zero conditions for *z* = 0,

The energy density *E* is a positive definite quadratic form of *ui,j* by construction. Therefore, by virtue of the decay of the solution at infinity, the relation only holds if *ui,j* = 0 for all *i* and*j.**u* = 0; i.e., the solutions coincide. The proof of the theorem is complete.

Theorem 12.1 holds for both exterior and interior boundary value problems.

## 13. Statement of supersonic BVP in D 3 ' R 3 and its generalized solution

To solve the problem, we also use the method of generalized functions. We introduce here the regular generalized function with support on

Also using the properties of differentiation of regular generalized functions with gaps at *D*, and taking into account the boundary conditions and the conditions on the fronts, we obtain the transport Lame equations (Eq. (14)) on the space of distributions with singular mass forces:

By using the properties of convolutions with the Green tensor and the boundary conditions, we obtain the generalized solution of BVP in the form:

By analog with the subsonic case, if we use fundamental stress tensor, then the right-hand side of Eq. (61) may be represented in the form of a surface integral over the boundary of the domain. In our notation, on the boundary, it acquires the form

This formula is similar to the Somigliana formula in the static theory of elasticity ([1]: 146), but it is impossible to use this formula to determine the solution of the boundary value problem in the case of supersonic loads, because the second term contains strong non-integrable singularities of the tensor *T* on the shock wave fronts of fundamental solutions; therefore, the integrals are divergent. To construct a regular integral representation of the formula, we must regularize it. For this, we use the tensor *H.*

## 14. Dynamic analog of the Somigliana formula in supersonic case

For regularization of Eq. (61), we put *W,z* instead of *U* in the second term and use the property of differentiation of convolution:

From here on, we use Eq. (57) we get the formula which can be written in integral form.

*The generalized solution of supersonic transport BVP can be presented in the form:*

*which for* *has the next integral presentation*

**Proof.** Formula (65) follows from Eq. (64) in virtue of Eqs. (61) and (32). Its integral form is

If we use Eq. (55) for

Note that

In virtue of this equity, we get from Eq. (66) the last formula of the theorem.

All integrals exist; indeed, the integrands are integrable everywhere, including the fronts of fundamental solutions, because the kernels of the integrands have weak singularities on the fronts of the form *U* and *H*. The proof is completed.

This formula is a dynamic analog of Somigliana formula for supersonic loads. It defines the displacement in elastic medium by using boundary values of stresses and velocity of displacements of boundary surface.

This formula also preserves its form for *D.*

## 15. Singular boundary integral equations of supersonic transport BVP

**Тheorem 15.1**. *If the classical solution of BVP satisfies the Holder’s conditions at D+, i.e., ∃ C > 0, β > 0 that*

*then it satisfies the singular boundary integral equation at D+*

**Proof.** The desired assertion follows from Theorem 14.1 and Theorem 11.1 for tensor *H* by analogy of the proof of Theorem 12.1 about singular boundary integral equations in the subsonic case. Full proofs of these theorems can be found in [7].

This theorem gives us a resolving system of integral equations for definition of unknown values of boundary displacements in the supersonic case.

Moreover, the Somigliana formula for displacements was obtained for generalized functions. But since they are regular, from the Dubois-Reymond lemma ([3]: 97), this solution is classical. However, if the acting loads are described by singular generalized functions, which often takes place in physical problems, then one should use a representation of a generalized solution in the convolution form (Eq. (65)) with the evaluation of convolutions by the definition (see [3]: 133).

## 16. Conclusion

The constructed singular boundary integral equations in the supersonic case are not classical equations because the solution inside a domain is determined by the boundary values of stresses and displacement rates rather than displacements themselves, unlike the Somigliana formula. In addition, the domain of integration over a boundary surface substantially depends on *z*, which is specific for hyperbolic equations. This complicates finding solutions of such problems by the successive approximation method. However, for the numerical discretization of singular boundary integral equations, the method of boundary elements makes it possible to use standard methods of computational mathematics for a computer implementation of the solution of such problems. The aforementioned boundary value problems model the dynamics of underground structures like transport tunnels and extended excavations subjected to the dynamic influence of moving vehicles and seismic loads. They permit one to study the dynamics of a rock mass in a neighborhood of underground structures depending on its physical-mechanical properties, the velocity of moving transport, specific features of the transport load, and the geometric properties of structures in technical computations of displacements and the stress-strain state of the mass away from the tunnel.