# Analysis of Disturbing Influence of Traffic Load on Soil Body

Authors: Janat Musayev; Algazy Zhauyt
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

Stress waves propagate in soil in case of earthquake and man-made effects (traffic flow, buried explosions, shield-driven pipes and tunnels, etc.). The wave point-sources are those located at the distances equal to more than two waves lengths, which significantly simplifies solving of a problem of these waves’ strength evaluation. Distribution of stress and displacement by the stress waves propagation in elastic medium is a complex pattern. The stress distribution in propagating waves depends on the type and form of source, conditions of the source contact with medium, and properties of mediums in the vicinity of the source. The point-sources and their combinations are selected in such a way to model an influence of machines and processes on soil body in case of shield-driven pipes (tunnels).


---

## Body

## 1. Introduction

In confined environment the wave pattern becomes more difficult due to the boundary reflection of waves. For evaluation of different factors influence on the nature of waves propagation it is convenient to divide the waves propagation problem into several stages.At the first stage tasks of waves propagation in infinite medium from the point-sources of the following different types are reviewed:(i)
concentrated force,
(ii)
double force without moment (two concentrated forces acting contrariwise along one line and applied at small distance 2 h from each other),
(iii)
two double forces without moment (double couple without moment) acting at right angles and in a single plane,
(iv)
three double forces without moment acting in three orthogonally related directions (center of expansion),
(v)
double forces with moment (couple of forces),
(vi)
combination of two couples of forces with sum of moments equal to zero,
(vii)
uniform pressures in alveole at a section in lengthd,
(viii)
tangential stresses applied at a section of alveole contour in lengthd,
(ix)
tangential stresses applied to alveole contour in lengthd.Solution of tasks enables evaluation of types of wave sources’ influence on the type and parameters of waves propagating over a distance.At the second stage tasks of wave propagation from the point-sources in semi-infinite elastic medium are reviewed. At this stage tasks of vibration of semi-infinite elastic medium surface from different point-sources of waves at various depths with account of waves reflected from free surface are solved. The solutions make it possible to evaluate a dynamic impact on the environment in case of earthquake, as well as the dynamic impact of different devices and machines used for underground works. This paper deals with a task of waves propagation from the point-sources of different types in the infinite space. In the performance of the tasks the Fourier integral transform and generalized functions are used [1–3].


## 2. Materials and Methods

For dynamic tasks both infinite and semi-infinite spaces and for bounded areas we will use differential equations of motion(1)
μ
U
j
,
i
i
+
λ
+
μ
U
i
,
i
j
-
ρ
U
¨
j
=
f
j
,
i
,
j
=
1,2
,
3
.Suppose that in the generalized functions outside areaΩ occupied by the elastic area the displacement and stress functions are equal to zero [1]. Then (1) may be in the form
(2)
μ
U
j
,
i
i
+
λ
+
μ
U
i
,
i
j
-
ρ
U
¨
j
=
μ
U
j
cos
⁡
n
,
x
i
δ
s
,
i
+
λ
+
μ
U
i
cos
⁡
n
,
x
i
δ
s
j
+
σ
j
i
s
cos
⁡
n
,
x
i
δ
s
-
ρ
U
k
t
=
0
×
δ
˙
t
+
ρ
U
j
t
=
T
δ
˙
t
-
T
-
ρ
U
˙
k
t
=
0
δ
t
+
ρ
U
˙
k
t
=
T
δ
t
-
T
,
where δ
i
j is Kronecker symbol, δ
s is delta function at the area boundary, and δ
t, δ
˙
t are delta function and its time derivative. U
j
s and σ
i
j
s are jump of U
j and σ
i
j functions upon outside passage through the boundary of area Ω. Since outside of this area these functions are equal to zero, symbols U
j and σ
i
j present values of these functions at the area boundary.FunctionsU
j
t
=
0, U
j
t
=
T, U
˙
j
t
=
0, U
˙
j
t
=
T represent initial and final conditions, that is, displacements and velocities of medium points at t
=
0, t
=
T.In the expressions below and above the summation over repeated indices is performed.Let us introduce the notations:(3)
X
j
=
μ
U
j
cos
⁡
n
,
x
i
δ
s
,
i
+
λ
+
μ
U
i
cos
⁡
n
,
x
i
δ
s
,
j
+
σ
j
i
s
cos
⁡
n
,
x
i
δ
s
-
ρ
U
k
t
=
0
×
δ
˙
t
+
ρ
U
j
t
=
T
δ
˙
t
-
T
-
ρ
U
k
t
=
0
×
δ
t
+
ρ
U
˙
k
t
=
T
δ
t
-
T
+
F
j
.Expression (3) contains information on mass load and impacts on the medium boundary and initial conditions. If the medium motion is considered at the finite interval (0, T) of time, the final conditions of displacements and velocities functions of the medium points at time T are included, which represent unknown quantities; in the used method the solutions fulfill a role of “integration constants.” F
j coincides with f
j in area Ω and in interval (0, T) and is equal to zero outside of this area and interval.Since all functions represented in the finite functions are equal to zero outside the area, consequently jumps of functions at the area boundary are boundary conditions. Depending on the assigned tasks a part of these functions is set; the other is defined in the course of solution [4].Let us apply the Fourier transform for the differential equation system solution. Let us multiply the left and right member of the equation byexp
⁡
i
α
k
x
k
+
ω
τ and integrate with respect to four variables: x
1
,
x
2
,
x
3
,
τ; in other words let us apply the Fourier transform to system of (2). Here α
1
,
α
2
,
α
3 mean parameters of the Fourier transform spacewise, ω being frequency.For description of the Fourier functionsσ
i
j
,
U
j
,
F
j let in symbol σ
-
i
j
,
U
-
j
,
F
-
j. The following properties of the Fourier transform are used upon integration: ∂
U
i
/
∂
x
j
⇒
-
i
α
j
U
~
i or in tonsorial notations U
i
,
j
⇒
-
i
α
j
U
~
i and U
¨
i
⇒
-
ω
2
U
~
i.Let us divide the left and right members of system of (2) by μ and introduce the notations: γ
2
=
α
1
2
+
α
2
2
+
α
3
2 and β
=
(
λ
+
2
μ
)
/
μ ratio of P-wave velocity c
1
=
(
λ
+
2
μ
)
/
ρ to S-wave velocity c
2
=
μ
/
ρ. Having used the agreed notations and the Fourier transform let us represent the system of (2) in the form
(4)
γ
2
-
ω
2
β
2
U
~
k
+
β
2
-
1
α
k
α
l
U
~
l
=
X
~
k
μ
.
In the right member of (4)  X
~
k includes not only the force impacts on medium but also kinematic ones.Solving of algebraic equations system (4) may be written down in the following form:
(5)
U
~
k
=
β
2
γ
2
-
ω
2
β
2
X
~
j
-
β
2
-
1
α
j
α
k
X
~
k
ρ
c
1
2
γ
2
-
ω
2
β
2
γ
2
-
ω
2
.
For detection of displacements it is necessary to perform an inverse Fourier transform:
(6)
U
j
x
1
,
x
2
,
x
3
,
τ
=
1
4
π
2
·
∫
W
β
2
γ
2
-
ω
2
β
2
X
~
j
-
β
2
-
1
α
j
α
k
X
~
k
ρ
c
1
2
γ
2
-
ω
2
β
2
γ
2
-
ω
2
=
1
4
π
2
·
·
e
-
x
k
α
k
+
ω
τ
d
W
,
where W means a space of variables α
1
,
α
2
,
α
3
,
ω and d
W
=
d
α
1
,
d
α
2
,
d
α
3
,
d
ω.For the infinite space symbolsX
j and X
~
j have only a generalized load of F
j and F
~
j.


## 3. Results and Discussion

Let us consider the method of solution of the tasks of waves propagation in the infinite elastic medium from the point-sources with application of the Fourier transform and generalized functions from the point-sources using the example of a task of the concentrated forceF
t effect for comparison with famous Love solution [5].Upon effect of the concentrated force at the origin of coordinates and inx
3-direction expressions (3) take the form
(7)
X
l
=
0
,
X
2
=
0
,
X
3
=
F
t
δ
x
3
,
and correspondingly the Fourier transform X
~
3
=
F
~
ω.Under such conditions the displacement components take the form(8)
U
1
x
1
,
x
2
,
x
3
,
τ
=
-
β
2
-
1
4
π
2
ρ
c
1
2
∫
W
α
1
α
3
X
~
3
γ
2
-
ω
2
β
2
γ
2
-
ω
2
=
-
β
2
-
1
4
π
2
ρ
c
1
2
·
e
-
x
k
α
k
+
ω
τ
d
W
,
U
2
x
1
,
x
2
,
x
3
,
τ
=
-
β
2
-
1
4
π
2
ρ
c
1
2
∫
W
α
2
α
3
X
~
3
γ
2
-
ω
2
β
2
γ
2
-
ω
2
=
-
β
2
-
1
4
π
2
ρ
c
1
2
·
e
-
x
k
α
k
+
ω
τ
d
W
,
U
3
x
1
,
x
2
,
x
3
,
τ
=
1
4
π
2
ρ
c
1
2
·
∫
W
X
~
3
β
2
γ
2
-
ω
2
β
2
-
β
2
-
1
α
3
2
γ
2
-
ω
2
β
2
γ
2
-
ω
2
=
-
β
2
-
1
4
π
2
ρ
c
1
2
·
e
-
x
k
α
k
+
ω
τ
d
W
.Let us introduce the notation:r
=
x
1
2
+
x
2
2
+
x
3
2.Considering radiation conditions and using asymptotic development [2] of the Fourier integrals (8) and neglecting components representing fluctuation of near-field we will get
(9)
U
1
=
1
4
π
ρ
r
∂
2
r
∂
x
1
∂
x
3
·
1
c
1
2
F
t
-
r
c
1
-
1
c
2
2
F
t
-
r
c
2
,
U
2
=
1
4
π
ρ
r
∂
2
r
∂
x
2
∂
x
3
·
1
c
1
2
F
t
-
r
c
1
-
1
c
2
2
F
t
-
r
c
2
,
U
3
=
1
4
π
ρ
r
∂
2
r
∂
x
3
2
·
1
c
1
2
F
t
-
r
c
1
-
1
c
2
2
F
t
-
r
c
2
.We point out that it is possible to get more complete transforms of integrals (8). In such case there are solutions fully coinciding with Love solutions [5] for the concentrated force in the infinite space. This paper deals with propagation of stress waves from the point-sources of different types in the infinite space based on the asymptotic development of Fourier integrals.In consequence of symmetry of displacements and stresses about axisx
3 it is possible to get rather convenient expressions if a spherical coordinate system going through the coordinates center is applied and φ will be defined as an angle between the radial coordinate and positive axis x
3:
(10)
U
r
=
cos
⁡
φ
4
π
ρ
r
c
1
2
F
t
-
r
c
1
,
U
θ
=
0
,
U
φ
=
sin
φ
4
π
ρ
r
c
2
2
F
t
-
r
c
2
.An interesting fact will be pointed out: over all distances from the point of force cross motions perpendicular to the force lineφ
=
π
/
2 exceed the longitudinal motion directed along the force line φ
=
0:
(11)
U
φ
r
,
π
/
2
U
r
r
,
0
=
c
2
2
c
1
2
.The achieved method is used upon getting of functions describing wave propagation from different sources. Wave radiation patterns from such sources are provided below with the required clarifications. Propagation of stress waves from the point-sources of different types acting in the infinite elastic medium. Figures1–10 have been obtained with use of MATLAB program complex [6–11].Diagrams of waves propagation from concentrated force.
(a)
P waves
(b)
S wavesDiagrams of waves propagation from combination of two forces.
(a)
P waves
(b)
S wavesFigure 3
Clarification to diagram of waves propagation from combination of two forces without moment.Diagrams of waves propagation from combination of two double forces.
(a)
P waves
(b)
S wavesFigure 5
Diagrams of waves propagation from combination of three couples of forces.Figure 6
Diagrams of waves propagation from two couples of forces.Diagrams of waves propagation from combination of two couples of forces with sum of moments equal to zero.
(a)
P waves
(b)
S wavesFigure 8
Diagrams of propagation of waves created by uniform pressure.Figure 9
Diagrams of propagation of waves created by tangential stress at alveole contour in lengthd.Figure 10
Propagation of waves created by tangential stress applied to alveole contour.


### 3.1. Waves Propagation from Concentrated Force

Consider(12)
u
r
=
cos
⁡
ϕ
4
π
ρ
c
1
2
r
F
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
sin
ϕ
4
π
ρ
c
2
2
r
F
t
-
r
c
2
.


### 3.2. Propagation of Waves from Combination of Two Forces (Double Force without Moment)

Consider(13)
u
r
=
2
h
cos
⁡
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
2
h
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


### 3.3. Propagation of Waves from Combination of Two Double Forces (Two Double Forces without Moment)

Consider(14)
u
r
=
2
h
sin
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
2
h
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


### 3.4. Propagation of Waves from Combination of Three Couples of Forces without Moment Directed in Parallel with Three Orthogonal Axes

Consider(15)
u
r
=
2
h
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
0
.


### 3.5. Propagation of Waves from Two Couples of Forces

Consider(16)
u
r
=
0
,
u
θ
=
-
2
h
sin
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
,
u
ϕ
=
0
.


### 3.6. Propagation of Waves from Combination of Two Couples of Forces with Sum of Moments Equal to Zero

Consider(17)
u
r
=
2
h
sin
θ
cos
⁡
θ
si
n
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
2
h
sin
ϕ
cos
⁡
2
θ
-
sin
2
θ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
,
u
ϕ
=
2
h
sin
θ
cos
⁡
θ
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


### 3.7. Propagation of Waves Created by Uniform Pressure in Alveoleat a Sections in Lengthd

Consider(18)
u
r
=
π
a
2
d
4
π
μ
c
1
r
1
-
2
c
2
2
c
1
2
cos
⁡
2
ϕ
p
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
π
a
2
d
θ
sin
ϕ
cos
⁡
ϕ
2
π
μ
c
2
r
p
′
t
-
r
c
2
.


### 3.8.  Propagation of Waves Created by Tangential Stress at Alveole Contour in Lengthd

Consider(19)
u
r
=
0
,
u
θ
=
π
a
2
d
4
π
μ
c
2
r
τ
′
t
-
r
c
1
,
u
ϕ
=
0
.


### 3.9. Propagation of Waves Created by Tangential Stress Applied to Alveole Contour at Section in Lengthd

Consider(20)
u
r
=
2
π
a
d
cos
⁡
ϕ
4
π
ρ
c
1
2
r
p
′
t
-
r
c
1
,
u
θ
=
0
,
u
θ
=
0
u
ϕ
=
2
π
a
d
sin
ϕ
4
π
ρ
c
2
2
r
τ
t
-
r
c
2
.Remark 1.
Derivative with respect to all expressions above between the brackets is marked with a prime.


## 3.1. Waves Propagation from Concentrated Force

Consider(12)
u
r
=
cos
⁡
ϕ
4
π
ρ
c
1
2
r
F
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
sin
ϕ
4
π
ρ
c
2
2
r
F
t
-
r
c
2
.


## 3.2. Propagation of Waves from Combination of Two Forces (Double Force without Moment)

Consider(13)
u
r
=
2
h
cos
⁡
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
2
h
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


## 3.3. Propagation of Waves from Combination of Two Double Forces (Two Double Forces without Moment)

Consider(14)
u
r
=
2
h
sin
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
2
h
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


## 3.4. Propagation of Waves from Combination of Three Couples of Forces without Moment Directed in Parallel with Three Orthogonal Axes

Consider(15)
u
r
=
2
h
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
0
.


## 3.5. Propagation of Waves from Two Couples of Forces

Consider(16)
u
r
=
0
,
u
θ
=
-
2
h
sin
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
,
u
ϕ
=
0
.


## 3.6. Propagation of Waves from Combination of Two Couples of Forces with Sum of Moments Equal to Zero

Consider(17)
u
r
=
2
h
sin
θ
cos
⁡
θ
si
n
2
ϕ
4
π
ρ
c
1
3
r
F
′
t
-
r
c
1
,
u
θ
=
2
h
sin
ϕ
cos
⁡
2
θ
-
sin
2
θ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
,
u
ϕ
=
2
h
sin
θ
cos
⁡
θ
sin
ϕ
cos
⁡
ϕ
4
π
ρ
c
2
3
r
F
′
t
-
r
c
2
.


## 3.7. Propagation of Waves Created by Uniform Pressure in Alveoleat a Sections in Lengthd

Consider(18)
u
r
=
π
a
2
d
4
π
μ
c
1
r
1
-
2
c
2
2
c
1
2
cos
⁡
2
ϕ
p
′
t
-
r
c
1
,
u
θ
=
0
,
u
ϕ
=
π
a
2
d
θ
sin
ϕ
cos
⁡
ϕ
2
π
μ
c
2
r
p
′
t
-
r
c
2
.


## 3.8.  Propagation of Waves Created by Tangential Stress at Alveole Contour in Lengthd

Consider(19)
u
r
=
0
,
u
θ
=
π
a
2
d
4
π
μ
c
2
r
τ
′
t
-
r
c
1
,
u
ϕ
=
0
.


## 3.9. Propagation of Waves Created by Tangential Stress Applied to Alveole Contour at Section in Lengthd

Consider(20)
u
r
=
2
π
a
d
cos
⁡
ϕ
4
π
ρ
c
1
2
r
p
′
t
-
r
c
1
,
u
θ
=
0
,
u
θ
=
0
u
ϕ
=
2
π
a
d
sin
ϕ
4
π
ρ
c
2
2
r
τ
t
-
r
c
2
.Remark 1.
Derivative with respect to all expressions above between the brackets is marked with a prime.


## 4. Conclusions

Solutions of stress waves propagation in elastic medium from different sources have been provided.The point-sources and their combinations have been selected in such a way to model the influence of machines and processes on soil body in case of shield-driven pipes (tunnels).P waves and S waves propagate from each source of waves.The waves propagate from concentrated force; the forms of these waves comply with the law of force variation. Wave amplitude decreases as1
/
r due to radiation in space.Upon propagation of waves from double force without moment the wave forms represent a derivative of functions describing the law of force variation. The wave amplitude also decreases as1
/
r due to radiation in space.The forms of waves propagating from different sources comply with the law of force variation and are derivatives of functions describing variation of forces and stresses.



---

*Source: 318289-2015-01-26.xml*