Tens of monographs devoted to the solution of various tasks of the classical theory of elasticity are published. Authors apply different methods of modern mathematics. However, modern mathematical programs for computers give the real chance to effectively apply one more method: statement and the solution of tasks in the matrix form after bidimensional Fourier’s transform. This work is devoted to it. We will consider the task for the layer on the elastic half-space which has the particular interest for sciences of the Earth.

Bidimensional Fourier’s transform will be used in the look

f ˜ ( ξ , η ) = ∫ − ∞ ∞ ∫ − ∞ ∞ f ( x , y ) e ( ξ x + η y ) i d ξ d η = ∫ 0 ∞ ∫ − π π r f ( x , y ) e i r ρ sin τ d τ d r (2.1)

where i = − 1 , r = x 2 + y 2 , ρ = ξ 2 + η 2 and in the second integral x = r ρ ( ξ sin τ + η cos τ ) , y = r ρ ( − ξ cos τ + η sin τ ) .

Then inverse transformation receives the form

f ( x , y ) = 1 4 π 2 ∫ − ∞ ∞ ∫ − ∞ ∞ f ˜ ( ξ , η ) e − ( ξ x + η y ) i d ξ d η = 1 4 π 2 ∫ 0 ∞ ∫ − π π ρ f ˜ ( ξ , η ) e i r ρ sin τ d τ d ρ (2.2)

wherein the second integral ξ = − ρ r ( x sin τ + y cos τ ) , η = ρ r ( x cos τ − y sin τ ) .

Transformations of derivatives are defined by formulas

F o ( ∂ n f ( x , y ) ∂ x n ) = ( − i ξ ) n ∂ f ˜ ∂ ξ , F o ( ∂ n f ( x , y ) ∂ y n ) = ( − i η ) n ∂ f ˜ ∂ η (2.3)

The following representations of Bessel functions are useful to inverse transformation of Fourier

∫ − π π e i z sin τ d τ = 2 π J 0 ( z ) , ∫ − π π i sin ( τ ) e i z sin τ d τ = − 2 π J 1 ( z ) , ∫ − π π cos ( τ ) e i z sin τ d τ = 0 (2.4)

Homogeneous equations of balance of homogeneous isotropic elastic environment in displacements have the form

drad   d i v   u + ( 1 − 2 ν ) ∇ 2 u = 0 (3.1)

where u = [ u , v , w ] is displacement vector, ∇ is the nabla-operator, ν is Poisson’s coefficient.

If to (3.1) to apply bidimensional Fourier’s transform on x, y, then the general solution of such system receives the form

u ˜ = ( χ A 1 + ξ z G ) e − ρ z + ( χ B 1 + ξ z Q ) e ρ z v ˜ = ( χ A 2 + η z G ) e − ρ z + ( χ B 2 + η z Q ) e ρ z w ˜ = ( χ A 3 − ρ z G i ) e − ρ z + ( χ B 3 + ρ z Q i ) e ρ z (3.2)

where A_i and B_i are arbitrary constants, G = − ξ A 1 + η A 2 − ρ A 3 i ρ , Q = ξ B 1 + η B 2 + ρ B 3 i ρ , χ = 3 − 4 ν .

For the first time such solution was received in work [1]. In the matrix type of the transform of movements and stresses take the form

u ˜ ( z ) = e − ρ z ρ ⋅ U × A + e ρ z ρ ⋅ V × B , s ˜ ( z ) = μ e − ρ z ρ ⋅ S × A + μ e ρ z ρ ⋅ P × B (3.3)

where μ is the shear modulus, × is multiplication of matrixes,

u ˜ ( z ) = ‖ u ˜ v ˜ w ˜ ‖ , s ˜ ( z ) = ‖ σ ˜ z x σ ˜ z y σ ˜ z z ‖ , A = ‖ A 1 A 2 A 3 ‖ , B = ‖ B 1 B 2 B 3 ‖ , α = 1 − 2 ν , β = 1 − ν

U = ‖ χ ρ − ξ 2 z − ξ η z ξ ρ z i − ξ η z χ ρ − η 2 z η ρ z i ξ ρ z i η ρ z i ρ ( χ + ρ z ) ‖ , V = ‖ χ ρ + ξ 2 z ξ η z ξ ρ z i ξ η z χ ρ + η 2 z η ρ z i ξ ρ z i η ρ z i ρ ( χ − ρ z ) ‖

S = ‖ 2 ξ 2 z ρ − χ ρ 2 − ξ 2 ξ η ( 2 ρ z − 1 ) − 2 ξ ρ ( ρ z + α ) i ξ η ( 2 ρ z − 1 ) 2 η 2 z ρ − χ ρ 2 − η 2 − 2 η ρ ( ρ z + α ) i − 2 ξ ρ ( ρ z − α ) i − 2 η ρ ( ρ z − α ) i 2 ρ 2 ( ρ z + 2 β ) ‖

P = ‖ 2 ξ 2 z ρ + χ ρ 2 + ξ 2 ξ η ( 2 ρ z + 1 ) 2 ξ ρ ( ρ z − α ) i ξ η ( 2 ρ z + 1 ) 2 η 2 z ρ + χ ρ 2 + η 2 2 η ρ ( ρ z − α ) i 2 ξ ρ ( ρ z + α ) i 2 η ρ ( ρ z + α ) i − 2 ρ 2 ( ρ z − 2 β ) ‖

We consider it expedient to remind that the algebra of matrixes differs from algebra of numbers. The main difference consists that matrix multiplication is noncommutativity.

Let’s consider the problem about the elastic layer on the surface of the elastic half-space. To simplify calculations, we will consider Poisson’s coefficient identical in the layer and the half-space and ν = 1 / 4 . At calculations it is necessary to carry out the main simplification: ξ 2 + η 2 = ρ 2 .

The layer − h ≤ z ≤ 0 is defined by formulas

u ˜ ( z ) = e − ρ z ρ ⋅ U × A + e ρ z ρ ⋅ V × B , s ˜ ( z ) = μ e − ρ z ρ ⋅ S × A + μ e ρ z ρ ⋅ P × B , − h ≤ z ≤ 0 (4.1)

and the half-space z ≥ 0

w ˜ ( x ) = e − ρ z ρ ⋅ U × C , g ˜ ( z ) = m μ e − ρ z ρ ⋅ S × C (4.2)

where w ˜ ( x ) is matrix of displacements, g ˜ ( z ) is matrix of stresses, mμ is the shear modulus in half-space, m ≠ 1 and C is matrix of arbitrary constants.

On a layer surface at z = − h the following single forces are possible:

σ z z ( x , y , − h ) = − δ ( x ) δ ( y ) , σ z x ( x , y , − h ) = σ z y ( x , y , − h ) = 0 , matrix is G;

σ z x ( x , y , − h ) = − δ ( x ) δ ( y ) , σ z z ( x , y , − h ) = σ z y ( x , y , − h ) = 0 , matrix is Gx;

σ z y ( x , y , − h ) = − δ ( x ) δ ( y ) , σ z z ( x , y , − h ) = σ z x ( x , y , − h ) = 0 , matrix is Gy.

where δ is delta-function.

The mentioned matrixes have the form

G = ‖ 0 0 − 1 ‖ , G x = ‖ − 1 0 0 ‖ , G y = ‖ 0 − 1 0 ‖ (4.3)

Thus, we have two boundary conditions: on demarcation of at z = 0 the continuity of movements and stresses is observed and on the surface of the layer at z = − h single force is applied. As a result, system of equations for definition of matrixes of A, B and C receive the look

u ˜ ( 0 ) = w ˜ ( 0 ) , s ˜ ( 0 ) = g ˜ ( 0 ) , s ˜ ( − h ) = G (4.4)

In the right-hand member of the third Equation (4.4) it is possible to substitute matrixes of Gx or Gy depending on the objective.

The first two equations of system (4.4) in expanded form have the form

A + B = C , S o × A + P o × B = m ⋅ S o × C (4.5)

where So and Po are matrixes S и P at z = 0.

The system (4.4) has the solution

B = ( m − 1 ) ⋅ N b × C , A = N a × C (4.6)

where N b = ( S o − 1 × P o − E ) − 1 , N a = E − ( m − 1 ) ⋅ N b .

Then the equations for the layer receive the form

u ˜ ( z ) = e − ρ z ρ ⋅ U × N a × C + e ρ z ( m − 1 ) ρ ⋅ V × N b × C s ˜ ( z ) = μ e − ρ z ρ ⋅ S × N a × C + μ e ρ z ( m − 1 ) ρ ⋅ P × N b × C (4.7)

and the third equation of system (4.4) receives the form

μ e ρ h ρ ⋅ S h × A + μ e − ρ h ρ ⋅ P h × B = G (4.8)

where Sh and Ph are matrixes S и P at z = − h .

Solution of (4.8) is

C = ρ μ ⋅ ( e ρ h ⋅ S h × N a + e − ρ h ⋅ ( m − 1 ) ⋅ P h × N b ) − 1 × G (4.9)

or

C = ρ μ ⋅ e − ρ h ⋅ C n = ρ μ ⋅ e − ρ h ⋅ ‖ 3 i ξ ( 4 ρ h ( m − 1 ) e − 2 ρ h + 2 ρ h m − 3 m + 4 ρ h ) 8 ρ 3 F 3 i η ( 4 ρ h ( m − 1 ) e − 2 ρ h + 2 ρ h m − 3 m + 4 ρ h ) 8 ρ 3 F − 3 ( 4 ( m − 1 ) ( ρ h − 1 ) e − 2 ρ h − 4 − 4 ρ h − 5 m − 2 ρ h m ) 8 ρ 2 F ‖ (4.10)

where F = 2 ( m − 1 ) 2 e − 4 ρ h + ( m − 1 ) ( 4 ρ 2 h 2 m + 5 m + 4 + 8 ρ 2 h 2 ) e − 2 ρ h             + ( m + 2 ) ( 2 m + 1 ) .

Now solutions for a layer and a half-space take a form

u ˜ ( z ) = e − ρ ( z + h ) μ ⋅ U × N a × C n + e − ρ ( h − z ) ( m − 1 ) μ ⋅ V × N b × C n

s ˜ ( z ) = e − ρ ( z + h ) ⋅ S × N a × C n + e − ρ ( h − z ) ( m − 1 ) ⋅ P × N b × C n (4.11)

w ˜ ( z ) = e − ρ ( z + h ) μ ⋅ U × C n , g ˜ ( z ) = m e − ρ ( z + h ) ⋅ S × C n

It is reasonable to carry out calculations on formulas (4.11) at the numerical values m. In the final solution it makes sense to make substitution z → z − h , and then the layer and the half-space will occupy areas 0 ≤ z ≤ h and z ≥ h .

We will not provide all solutions which can be received on formulas (4.11). Let’s give only the formula for displacements in the half-space. Calculations were made at m = 5 and substitution z → z − h was made. Displacements have the form

w ˜ ( z ) = e − ρ z 2 μ ‖ 3 i ξ ( f 1 ρ z + f 2 ) 2 ρ 2 f 0 3 i η ( f 1 ρ z + f 2 ) 2 ρ 2 f 0 3 ( f 1 ρ z + f 3 ) 2 ρ f 0 ‖ (5.1)

where z ≥ h , f 0 = 32 e − 4 ρ h + ( 112 ρ 2 h 2 + 116 ) e − 2 ρ h + 77 , f 1 = ( 8 − 16 ρ h ) e − 2 ρ h + 22 , f 2 = ( 8 ρ h + 16 ρ 2 h 2 ) e − 2 ρ h − 8 ρ h − 15 , f 3 = ( 16 ρ 2 h 2 − 24 ρ h + 16 ) e − 2 ρ h − 8 ρ h + 29 .

Inverse transformations are carried out by means of section 2 formulas. Let’s give two examples for arbitrary H(ρ) function

F o − 1 ( H ( ρ ) ) = 1 4 π 2 ∫ 0 ∞ ρ H ( ρ ) ∫ − π π e i r ρ sin ( τ ) d τ d ρ = 1 2 π ∫ 0 ∞ ρ H ( ρ ) J 0 ( r ρ ) d ρ (5.2)

and

F o − 1 ( i ξ ρ H ( ρ ) ) = 1 4 π 2 ∫ 0 ∞ ρ H ( ρ ) ∫ − π π i ξ ρ e i r ρ sin ( τ ) d τ d ρ = − 1 4 π 2 r ∫ 0 ∞ ρ H ( ρ ) ∫ − π π i ( x sin τ + y cos τ ) e i r ρ sin ( τ ) d τ d ρ = − x 4 π 2 r ∫ 0 ∞ ρ H ( ρ ) ∫ − π π i sin τ e i r ρ sin ( τ ) d τ d ρ = x 2 π r ∫ 0 ∞ ρ H ( ρ ) J 1 ( r ρ ) d ρ (5.3)

Calculation of the received integrals takes 1 - 2 seconds of machine time and allows to investigate the decision in detail.

By the offered method, it is possible to solve several other problems for the layer. The solution in all cases turns out rather simple and is uniform. It is sometimes simpler to solve again the known problem, than to look for its solution in literature.

The author declares no conflicts of interest regarding the publication of this paper.

Dobrovolsky, I.P. (2021) Elastic Layer on the Elastic Half-Space: The Solution in Matrixes. Open Access Library Journal, 8: e7191. https://doi.org/10.4236/oalib.1107191