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).
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:
concentrated force,
double force without moment (two concentrated forces acting contrariwise along one line and applied at small distance 2 h from each other),
two double forces without moment (double couple without moment) acting at right angles and in a single plane,
three double forces without moment acting in three orthogonally related directions (center of expansion),
double forces with moment (couple of forces),
combination of two couples of forces with sum of moments equal to zero,
uniform pressures in alveole at a section in length d,
tangential stresses applied at a section of alveole contour in length d,
tangential stresses applied to alveole contour in length d.
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)μUj,ii+λ+μUi,ij-ρU¨j=fj,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)μUj,ii+λ+μUi,ij-ρU¨j=μUjcosn,xiδs,i+λ+μUicosn,xiδsj+σjiscosn,xiδs-ρUkt=0×δ˙t+ρUjt=Tδ˙t-T-ρU˙kt=0δt+ρU˙kt=Tδt-T,
where δij is Kronecker symbol, δs is delta function at the area boundary, and δt, δ˙t are delta function and its time derivative. Ujs and σijs are jump of Uj and σij functions upon outside passage through the boundary of area Ω. Since outside of this area these functions are equal to zero, symbols Uj and σij present values of these functions at the area boundary.
Functions Ujt=0, Ujt=T, U˙jt=0, U˙jt=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)Xj=μUjcosn,xiδs,i+λ+μUicosn,xiδs,j+σjiscosn,xiδs-ρUkt=0×δ˙t+ρUjt=Tδ˙t-T-ρUkt=0×δt+ρU˙kt=Tδt-T+Fj.
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.” Fj coincides with fj 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 by expiαkxk+ωτ and integrate with respect to four variables: x1,x2,x3,τ; 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 σij,Uj,Fj let in symbol σ-ij,U-j,F-j. The following properties of the Fourier transform are used upon integration: ∂Ui/∂xj⇒-iαjU~i or in tonsorial notations Ui,j⇒-iαjU~i and U¨i⇒-ω2U~i.
Let us divide the left and right members of system of (2) by μ and introduce the notations: γ2=α12+α22+α32 and β=(λ+2μ)/μ ratio of P-wave velocity c1=(λ+2μ)/ρ to S-wave velocity c2=μ/ρ. Having used the agreed notations and the Fourier transform let us represent the system of (2) in the form
(4)γ2-ω2β2U~k+β2-1αkαlU~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β2X~j-β2-1αjαkX~kρc12γ2-ω2β2γ2-ω2.
For detection of displacements it is necessary to perform an inverse Fourier transform:
(6)Ujx1,x2,x3,τ=14π2·∫Wβ2γ2-ω2β2X~j-β2-1αjαkX~kρc12γ2-ω2β2γ2-ω2=14π2··e-xkαk+ωτdW,
where W means a space of variables α1,α2,α3,ω and dW=dα1,dα2,dα3,dω.
For the infinite space symbols Xj and X~j have only a generalized load of Fj 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 force Ft effect for comparison with famous Love solution [5].
Upon effect of the concentrated force at the origin of coordinates and in x3-direction expressions (3) take the form
(7)Xl=0,X2=0,X3=Ftδx3,
and correspondingly the Fourier transform X~3=F~ω.
Under such conditions the displacement components take the form
(8)U1x1,x2,x3,τ=-β2-14π2ρc12∫Wα1α3X~3γ2-ω2β2γ2-ω2=-β2-14π2ρc12·e-xkαk+ωτdW,U2x1,x2,x3,τ=-β2-14π2ρc12∫Wα2α3X~3γ2-ω2β2γ2-ω2=-β2-14π2ρc12·e-xkαk+ωτdW,U3x1,x2,x3,τ=14π2ρc12·∫WX~3β2γ2-ω2β2-β2-1α32γ2-ω2β2γ2-ω2=-β2-14π2ρc12·e-xkαk+ωτdW.
Let us introduce the notation: r=x12+x22+x32.
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)U1=14πρr∂2r∂x1∂x3·1c12Ft-rc1-1c22Ft-rc2,U2=14πρr∂2r∂x2∂x3·1c12Ft-rc1-1c22Ft-rc2,U3=14πρr∂2r∂x32·1c12Ft-rc1-1c22Ft-rc2.
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 axis x3 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 x3:
(10)Ur=cosφ4πρrc12Ft-rc1,Uθ=0,Uφ=sinφ4πρrc22Ft-rc2.
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,π/2Urr,0=c22c12.
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. Figures 1–10 have been obtained with use of MATLAB program complex [6–11].
Diagrams of waves propagation from concentrated force.
P waves
S waves
Diagrams of waves propagation from combination of two forces.
P waves
S waves
Clarification to diagram of waves propagation from combination of two forces without moment.
Diagrams of waves propagation from combination of two double forces.
P waves
S waves
Diagrams of waves propagation from combination of three couples of forces.
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.
P waves
S waves
Diagrams of propagation of waves created by uniform pressure.
Diagrams of propagation of waves created by tangential stress at alveole contour in length d.
Propagation of waves created by tangential stress applied to alveole contour.
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 as 1/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 as 1/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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
KurbatskiiE. H.1995Moscow, RussiaMIITVerruijtA.BookerJ. R.Surface settlements due to deformation of a tunnel in an elastic half plane19964647537562-s2.0-003046862810.1680/geot.1996.46.4.753LoganathanN.PoulosH. G.Analytical prediction for tunneling-induced ground movements in clays1998124984685610.1061/(asce)1090-0241(1998)124:9(846)2-s2.0-0032170289Wan gJ. G.KongS. L.LeungC. F.Twin tunnels-induced ground settlement in soft soils2003Tsinghua University Press241244MyshkisA. D.Mathematics for technical colleges1996156632644LeeK. M.RoweR. K.LoK. Y.Subsidence owing to tunnelling. I. Estimating the gap parameter199229692994010.1139/t92-1042-s2.0-0027010176ParkK.-H.Analytical solution for tunnelling-induced ground movement in clays200520324926110.1016/j.tust.2004.08.0092-s2.0-11944270443FranziusJ. N.2003London, UKDepartment of Civil and Environmental Engineering, Imperial College of Science, Technology and MedicineJohnO. B.KeuselT. R.KingE. H.1996Chapman & HallLeblaisY.1995Association Française des Travaux en Souterrain (AFTES)HulmeT. W.ShirlawJ. N.HwangR. N.Settlements during the underground constructions of the Singapore MRTProceedings of the 10th Southeast Asian Geo-Technical Conference1990Taipei, Taiwan1620