Open access peer-reviewed chapter

Cylindrical Surface Wave: Revisiting the Classical Biot’s Problem

Written By

Jeremiah Rushchitsky

Submitted: 22 February 2019 Reviewed: 17 May 2019 Published: 27 September 2019

DOI: 10.5772/intechopen.86910

From the Edited Volume

Seismic Waves - Probing Earth System

Edited by Masaki Kanao and Genti Toyokuni

Chapter metrics overview

863 Chapter Downloads

View Full Metrics

Abstract

The problem on a surface harmonic elastic wave propagating along the free surface of cylindrical cavity in the direction of cavity axis is considered. In the case of isotropic medium, this is the classical Biot’s problem of 1952. First, the Biot pioneer work is revisited: the analytical part of Biot’s findings is shown in the main fragments. The features are using two potentials and representation of solution by Macdonald functions of different indexes. Then the new direct generalization of Biot’s problem on the case of transversely isotropic medium within the framework of linear theory of elasticity is proposed. Transition to the transverse isotropy needs some novelty—necessity of using the more complex representations of displacements through two potentials. Finally, a generalization of Biot’s problem on the case of isotropic and transversely isotropic media in the framework of linearized theory of elasticity with allowance for initial stresses is stated. This part repeats briefly the results of A.N. Guz with co-authors of 1974. The main features are using the linearized theory of elasticity and one only potential. All three parts are shown as analytical study up to the level when the numerical methods have to be used.

Keywords

  • surface harmonic cylindrical wave
  • classical Biot’s problem
  • generalization to the case of transversely isotropic medium

1. Introduction

Note first that the seismic waves include mainly the primary and secondary body waves and different kinds of surface waves. This chapter is devoted to one kind of surface waves. The problem is stated as follows: the infinite medium with cylindrical circular cavity having the symmetry axis Oz and constant radius is analyzed. An attenuating in depth of medium surface harmonic wave propagates along the cavity surface in direction Oz . In this case, the problem becomes mathematically the axisymmetric one. This problem is solved by Biot in 1952 [1] with assumption that the medium is isotropic. The context of this chapter includes four parts. The subchapter 1 “Introduction” is the standard one. The subchapter 2 is named: “Main Stages of Solving the Classical Biot’s Problem on Surface Wave along Cylindrical Cavity.” Here, the analytical part of Biot’s findings is shown in the main fragments. The features are using two potentials and representation of solution by Macdonald functions of different indexes. The subchapter 3 “Direct Generalization of Biot’s Problem on the Case of Transversely Isotropic Media within the Framework of Linear Theory of Elasticity” contains the new approach to the classical Biot’s problem and represents the direct generalization of this problem that uses the Biot’s scheme of analysis. Transition to the case of transverse isotropy needs some novelty—necessity of using the more complex representations of displacements through two potentials. The subchapter 4 “Genera-lization of Biot’s Problem on the case of Isotropic and Transversely Isotropic Media within the framework of Linearized Theory of Elasticity with Allowance for Initial Stresses” repeats briefly the results of A.N. Guz with co-authors (1974). They considered a generalization of the Biot’s problem on the case of elastic media with allowance for the initial stresses. The main features are using the linearized theory of elasticity, one only potential, and Macdonald function of one index.

Advertisement

2. Main stages of solving the classical Biot’s problem on surface wave along a cylindrical cavity

2.1 Statement of problem and main equations in potentials

A cylindrical system of coordinates Orϑz is chosen, and a harmonic wave is considered that has the phase variable σ = k z vt , unknown wave number k = ω / v , unknown phase velocity v , and arbitrary (but given) frequency ω and ampitude A . It is supposed that the wave propagates in an infinite medium with cylindrical cavity of constant radius r o in the direction of vertical coordinate z and possibly attenuates in the direction of radial coordinate r . In this linear statement and in assumption that deformations are small, the problem is axisymmetric, and deformations are described by two displacements u r r z t u φ r z t = 0 and u z r z t and two Lame equations of the form

C 11 C 12 2 Δ rz u r 1 r 2 u r + C 11 + C 12 2 u r , r + 1 r u r + u z , z , r = ρ u r , tt , E1
1 2 C 11 C 12 Δ rz u z + 1 2 C 11 + C 12 u r , r + 1 r u r + u z , z , z = ρ u z , tt , E2

or

λ + 2 μ u r , rr + 1 r u r , r 1 r 2 u r + u z , rz + μ u r , zz u z , rz = ρ u r , tt E3
λ + 2 μ u r , rz + 1 r u r , z + u z , zz μ 1 r u r , z u z , r + u r , rz u z , rr = ρ u z , tt . E4

Further the potentials Φ r z t , Ψ r z t are introduced

u r = Φ , r Ψ , z , u z = Φ , z + Ψ , r + 1 / r Ψ . E5

When Eq. (5) is substituted into Eqs. (3) and (4), then two uncoupled linear wave equations are obtained:

Δ rz Φ 1 / v L 2 Φ , tt = 0 , E6
Δ rz Ψ 1 / r 2 Ψ 1 / v T 2 Ψ , tt = 0 . E7

Here the standard notations of Laplace operator Δ rz and velocities of longitudinal and transverse waves in isotropic elastic medium v L = λ + 2 μ / ρ , v T = μ / ρ are used.

2.2 Solving the wave equations in the form of Macdonald functions

The solution of Eqs. (6) and (7) is found in the form of harmonic waves in the direction of vertical coordinate:

Φ r z t = Φ r e i kz ωt , Ψ r z t = Ψ r e i kz ωt , E8

Φ r z t = Φ r cos k z vt , Ψ r z t = Ψ r sin k z vt .

A substitution of representations (8) into the wave Eqs. (6) and (7) gives the equations relative to the unknown amplitudes Φ r , Ψ r Φ , rr + 1 / r Φ , r k 2 k L 2 Φ = 0 .

Φ , rr + 1 / r Φ , r k 2 1 v / v L 2 Φ = 0 , E9
Ψ , rr 1 / r Ψ , r k 2 k T 2 + 1 / r 2 Ψ = 0 Ψ , rr 1 / r Ψ , r k 2 1 v / v T 2 + 1 / r 2 Ψ = 0 E10

These equations correspond to the Bessel equation for Macdonald functions K λ x (modified Bessel functions of the second kind [2, 3, 4])

y + 1 / x y 1 + λ 2 / x 2 y = 0 E11

More exactly, Eqs. (9) and (10) have the solutions in the form of Macdonald functions, if the conditions.

k 2 k L 2 > 0 , k 2 k T 2 > 0 k 2 1 v / v L 2 > 0 k 2 1 v / v T 2 > 0 E12

are fulfilled. According to conditions (12), the wave number of cylindrical wave must be real, and the wave velocity must be less of the velocities of classical longitudinal and transverse plane waves.

Further the wave Eqs. (9) and (10) are considered separately. The first equation is written in the form

Φ , rr + 1 / r Φ , r m L 2 Φ = 0 m L = k 1 v / v L 2 E13

This equation has the solution in the form of Macdonald function:

Φ r = A Φ K 0 m L r E14

of zeroth order and unknown argument x = m L r , which includes the unknown phase velocity of wave.

The second equation can be written in the form

Ψ , rr 1 / r Ψ , r m T 2 + 1 / r 2 Ψ = 0 m T = k 1 v / v T 2 . E15

The corresponding solution under conditions (12) is expressed by the Macdonald function K 1 r k 2 k T 2

Ψ r = A Ψ K 1 m T r E16

of the first order and unknown argument x = m T r , which includes the unknown wave velocity. The amplitude coefficient A Ψ is assumed to be constant and arbitrary.

Note that the Macdonald functions have the property of attenuation with increasing arguments which is shown in Figure 1. Therefore, the propagation along the vertical coordinate z waves (15) and (16) can be considered as the waves with amplitudes Φ r , Ψ r , which attenuate with increasing the radial coordinate r .

Figure 1.

Plots of the first five Macdonald functions.

This means that amplitudes can decrease essentially with increasing the distance from the surface of cylindrical cavity. In this sense, the waves (15) and (16) are the surface ones. This forms also the sense of conditions (12). The same conditions are used in the analysis of classical Rayleigh surface wave which propagates along the plane surface of isotropic elastic medium [5, 6, 7, 8, 9]. But the Rayleigh wave attenuates as an exponential function when being moved from the free surface, whereas the cylindrical surface Biot’s wave attenuates as the Macdonald functions. At that, the arguments in exponential function and Macdonald functions are identical and depend on the wave velocity.

2.3 Boundary conditions: equations for unknown wave number

The boundary conditions correspond to the absence of stresses on surface r = r o

σ rr r = r o z t = 0 , σ rz r = r o z t = 0 . E17

The stresses

σ rr = 2 μ u r , r + λ u r / r + u r , r + u z , z , σ rz = μ u r , z + u z , r E18

are written through the potentials

σ rr = λ + 2 μ Φ , rr Ψ , rz + λ 1 / r Φ , r Ψ , z + Φ , zz + Ψ , rz + 1 / r Ψ , z , E19
σ rz = μ Φ , rz Ψ , zz + Φ , zr + Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ . E20

Then the boundary conditions (17) can be written in the form.

2 μ Φ , rr Ψ , rz + λ ΔΦ r = r o = 0 , μ 2 Φ , rz Ψ , zz + ΔΨ 1 / r 2 Ψ r = r o = 0 E21

In the work [1], Biot has used the expressions.

ΔΦ 1 / v L 2 Φ , tt = 0 , ΔΨ 1 / r 2 Ψ 1 / v T 2 Ψ , tt = 0 and rewrite Eq. (21) in such a way Φ , rr Ψ , rz + λ / 2 μ 1 / v L 2 Φ , tt r = r o = 0 , 2 Φ , rz Ψ , zz + 1 / v T 2 Ψ , tt r = r o = 0 .

Then the substitution of solutions (14) and (16) into the boundary conditions (21) gives two homogeneous algebraic equations relative to the unknown constant amplitude coefficients

1 v / v L 2 λ μ v / v L 2 K 0 m L r o K 0 m L r o + K 2 m L r o A Φ 1 v / v L 2 A Ψ = 0 , E22
2 1 v / v L 2 A Φ + 2 v / v T 2 K 1 m T r o K 1 m L r o A Ψ = 0 . E23

An analysis of these equations that describe the cylindrical surface wave is very similar to the analysis that has been carried out by Rayleigh for the classical wave propagating along the plane surface. Some novelty in analysis of systems (22) and (23) is consideration of the system relative to quantities K 1 m L r o A Φ and K 1 m S r o A Ψ

1 v / v L 2 K 0 m L r o K 1 m L r o + 1 m L r o λ 2 μ v / v L 2 K 0 m L r o K 1 m L r o K 1 m L r o A Φ + 1 v / v S 2 K 0 m T r o K 1 m T r o + 1 m T r o K 1 m T r o A Ψ = 0 , E24
2 1 v / v L 2 K 1 m L r o A Φ + 2 v / v T 2 K 1 m T r o A Ψ = 0 . E25

Solving of systems (24) and (25) gives two results. First, the solution is found accurate within one amplitude factor. Second, an equation for determination of phase velocity of cylindrical surface wave can be obtained in an explicit form.

The work of Biot (1952) has demonstrated some art in handling the Macdonald functions and has written Eq. (24) through only functions of the zeroth and first orders. For that, the known formulas

K 0 x = K 1 x , K 1 x = K 0 x , K 0 x + 1 / x K 0 x = K 0 x , K 0 x = 1 / x K 1 x + K 0 x E26

have been used [3]. As a result, the equation for determination of phase velocity of cylindrical wave has the form

2 v / v T 2 2 v / v T 2 K 0 m L r o K 1 m L r o + 1 v / v L 2 m L r o 4 1 v / v L 2 1 v / v T 2 K 0 m T r o K 1 m T r o + 1 m T r o = 0 . E27

Let us write the corresponding equation for the Rayleigh wave [5, 6, 7, 8, 9] as

4 1 v / v L 2 1 v / v S 2 2 v / v S 2 2 = 0 . E28

Thus, a presence of Macdonald functions in Eq. (27) complicates essentially an analysis of this equation because according to relations m L = k 1 v / v L 2 , m S = k 1 v / v T 2 these functions have the unknown velocity in argument.

If the cavity radius is not small, then the Macdonald functions can be represented by the simple formula K 0 r = K 1 r = e r π / 2 r , and Eq. (27) is reduced to the Rayleigh Eq. (28).

Strictly speaking, the analytical part of analysis is ended by obtaining Eq. (27). Further analysis can be continued with the aim of the numerical methods. Biot in [1] has shown some comments and conclusions based on resources of the 1950s.

A possibility of analytical approach is still saved in the problem on existence of the appropriate wave velocity. First of all, Eq. (27) depends on the elastic constants, and this dependence can be shown in the form of dependence on the ratio of known velocities v L / v T . If the notation v 2 / v T 2 = z is used, then Eq. (27) can be written in the form

2 z v L / v T 2 2 z K 0 r o k 1 z v L / v T 2 K 1 r o k 1 z v L / v T 2 + 1 z v L / v T 2 r o k 1 z v L / v T 2 4 1 z 1 z v L / v T 2 K 0 r o k 1 z K 1 r o k 1 z + 1 r o k 1 z = 0 . E29

It seems appropriate to recall here the most known ways of proving the existence of velocity of the classical Rayleigh wave. An initial equation is always Eq. (28). Two different notations v 2 / v T 2 = z and v = 1 / θ are used, which generate two different representations of Eq. (28)

z z 3 8 z 1 z 2 1 v T 2 / v L 2 = 0 , E30
2 θ 2 1 / v T 2 2 4 θ 2 θ 2 1 / v T 2 θ 2 1 / v L 2 = 0 . E31

Finding the real root of Eq. (30) is the key step in the analysis of the Rayleigh wave [5, 6, 7, 8, 9] . For more than 100 years of analysis of this wave, many methods of proving the existence of real velocity of wave were elaborated.

First of all, the sufficiently useful and exact empirical Viktorov’s formula [5].

v / v T = z 0.87 + 1.12 υ 1 + υ υ the Poisson ratio E32

should be shown.

Let us show further briefly some phenomenological methods. Note that the restriction on the Rayleigh wave velocity is already obtained from a statement of the problem—it is less of the velocity of plane transverse wave. This restriction can be written in the form z < 1 or θ > 1 / c T .

Method 1 (graphical method [10, 11]). Eq. (30) is considered as a sum of two summands Z 1 + Z 2 = 0 . The first summand Z1 = z3 describes a cubic parabola; the lower branch of which lies in the first quadrant of the plane zOZ 1 . The second summand describes a quadratic parabola Z 2 = 8 z 1 z 2 1 c T 2 / c L 2 , which is concave in the direction of coordinate axis OZ 2 . Further the ratio c T 2 / c L 2 = μ / λ + 2 μ can be estimated from below and top 0 c T 2 / c L 2 1 / 2 with allowance for the shear modulus μ that is positive. These parabolas are intersected on the interval 0 1 . More exactly, one of the roots z = z C of Eq. (30) can be estimated 0.764 z = c / c T 2 0.912 . Here, the minimal value corresponds to the case when the parabola is tangent to the abscissa axis, and the maximal value corresponds to the case when the parabola is moved partially into the fourth quadrant. Thus, the velocity of Rayleigh wave is close to the velocity of plane transverse wave, but always less of its 0.874 ≤ (c/cT) ≤ 0.955.

Method 2 (method of finding the interval, on ends of which the equation possesses the different by sign values [2, 11]). This method is based on the analysis of Eq. (30). The value of equation that corresponds to the point c R = c T is positive and equal to 1. The second point is chosen as c R = ε c T , where ε is assumed as the small quantity (this point is close to 0). When this value is substituted into Eq. (30), then expression 2 ε 2 1 c T 2 / c L 2 is always negative. Hence, at least one root of equation lies in the interval ε c T c T .

Method 3 (another method of finding the interval, on ends of which the equation possesses the different by sign values [5]). This method is based on the analysis of Eq. (31). The right point is chosen as θ = 1 / c T (similar to method 2). Then Eq. (31) possesses the positive value. The left point corresponds to θ . Further an expression (31) is expanded into the power series near the point at infinity. This series starts with the term 2 θ 2 1 / c T 2 1 / c L 2 , which is always negative. So this equation possesses in the chosen points the different sign values. Thus, at least one root of the equation lies in the interval 1 / c T .

Method 4 (method based on assumption relative to the Poisson ratio [7]). This assumption consists in the choice of value of Poisson ratio that is often used in the analysis of seismic waves in Earth’s crust ν = λ / 2 λ + μ = 1 / 4 λ = μ . Then cubic Eq. (31) (the zeroth root θ 1 = 0 is ignored from a physical considerations) can be solved exactly, and the roots possess the values θ 2 = 4 , θ 3 = 2 + 2 / 3 , θ 4 = 2 2 / 3 . Since the condition θ < 1 has been fulfilled, then the corresponding root is equal to θ 4 = 0.8453 .

The main conclusion from the shown above methods is that they really allow to establish an existence of real root of Rayleigh equation (the real value of velocity of harmonic Rayleigh wave). They give the positive answer on the question whether the Rayleigh wave exists. In the case of other surface waves including the cylindrical wave under consideration, the experience of the classical Rayleigh wave analysis can be quite useful.

Advertisement

3. Cylindrical wave propagating along the surface of the cylindrical cavity in the direction of vertical axis: The case of transversal isotropy of medium

Let us return to the initial statement of problem and consider an infinite medium with cylindrical circular cavity that has the symmetry axis Oz and radius r o . The medium is assumed to be the transversely isotropic elastic one. It is assumed further that the wave is harmonic in time, and attenuating deep into medium wave propagates in the direction of axis Oz along the cavity surface. Such a problem can be considered as some generalization of Biot’s [1] problem that is solved in the assumption of isotropy of medium on the case of transversal isotropy of medium. Therefore, it seems expedient to recall some facts from the theory of elasticity of transversally isotropic medium.

3.1 Some information on transversally isotropic medium

Let us consider the case when Ox 3 is the axis of symmetry and Ox 1 x 2 is the plane of isotropy. This symmetry corresponds to the hexagonal crystalline system. The matrix of elastic properties is characterized by 5 independent elastic constants C 11 , C 12 , C 13 , C 33 , C 44 and 12 non-zero components [11, 12, 13]:

C IK = C 11 C 12 C 13 0 0 0 C 12 C 11 C 13 0 0 0 C 13 C 13 C 33 0 0 0 0 0 0 C 44 0 0 0 0 0 0 C 44 0 0 0 0 0 0 1 / 2 C 11 C 12 . E33

Then the constitutive relations σ ε have the form [12, 14].

σ 11 = C 11 kl ε kl = C 11 ε 11 + C 12 ε 22 + C 13 ε 33 , σ 22 = C 22 kl ε kl = C 12 ε 11 + C 22 ε 22 + C 13 ε 33 , σ 33 = C 33 kl ε kl = C 13 ε 11 + C 13 ε 22 + C 33 ε 33 , σ 12 = C 11 C 12 ε 12 , σ 13 = 2 C 44 ε 13 , σ 23 = 2 C 44 ε 23 , E34

or in notations σ u [12, 14].

σ 11 = C 11 u 1 , 1 + C 12 u 2 , 2 + C 13 u 3 , 3 , σ 22 = C 12 u 1 , 1 + C 11 u 2 , 2 + C 13 u 3 , 3 , σ 33 = C 13 u 1 , 1 + C 13 u 2 , 2 + C 11 u 3 , 3 , σ 12 = 1 / 2 C 11 C 12 u 1 , 2 + u 2 , 1 , σ 13 = C 44 u 1 , 3 + u 3 , 1 , σ 23 = 1 / 2 C 44 u 2 , 3 + u 3 , 2 . E35

Also, five independent elastic technical constants are often used.

E x = E y , E x = E y , E z , G xy , G xz = G yz , υ xy , υ xz = υ yz , G xy = E x / 1 + 2 υ xy . They are evaluated through C NM by the following formulas:

The longitudinal Young modulus that corresponds to tension along the symmetry axis Oz

E z = C 33 2 C 13 2 / C 11 + C 12 . E36

The transverse Young modulus that corresponds to tension in the isotropy plane Oxy

E x = C 11 C 12 C 11 + C 12 C 33 2 C 13 2 / C 11 C 33 + C 13 2 . E37

The shear modulus that corresponds to the shear along the isotropy plane Oxy

G xy = C 66 = 1 / 2 C 11 C 12 . E38

The shear modulus that corresponds to the shear along the symmetry axis Oz

G xz = C 44 . E39

The Poisson ratio that corresponds to the shear along the symmetry axis Oz under tension in the isotropy plane and characterizes the shortening in this plane

υ xz = C 13 / C 11 + C 12 . E40

Sometimes, the corresponding Lame moduli are used.

λ xy + 2 μ xy = C 11 , λ xy = C 12 , μ xy = 1 / 2 C 11 C 12 λ xz + 2 μ xz = C 33 , λ xz = C 13 , μ xz = C 44 . E41

The Poisson ratio (40) is determined by the known formula of isotropic theory υ xz = λ xz / 2 λ xy + μ xy .

The Poisson ratio υ xy that corresponds to the shear along the symmetry axis Oz under tension along the isotropy plane is determined also by the classical formula υ xy = λ xy / λ xy + μ xy .

The constants C 11 , C 12 , C 13 , C 33 , C 44 are represented through the technical constants E , E , ν , ν , G by the formulas.

C 11 = 1 ν 2 E / E 1 ν 2 + 1 + 2 ν ν 2 E / E E , C 12 = v ν 2 E / E 1 ν 2 + 1 + 2 ν ν 2 E / E E , C 13 = ν 1 ν 1 ν 2 + 1 + 2 ν ν 2 E / E E , C 33 = 1 ν 2 1 ν 2 + 1 + 2 ν ν 2 E / E E , C 44 = G . E42

Let us comment briefly some features of transversally isotropic materials. They can be divided on the natural and artificial ones. An example of the classical natural material is the rock. An example of the modern material is a family of fibers “Kevlar®.” Kevlar® KM2 [15] is characterized by elastic constants E x = 1.34 GPa , E z = 84.62 GPa , G xz = 24.40 GPa , υ xy = 0.24 , υ xz = 0.60 .

An example of composite materials can be four fibrous composites of micro- and nanolevels, which are described in [15]. The corresponding elastic constants for some variants of these materials are as follows [15]:

10% of carbon microfibers E x = 3.59 GPa , E z = 25.22 GPa , G xz = 1.17 GPa , υ xy = 0.39 , υ xz = 0.58 .

10% of graphite microwhiskers E x = 3.69 GPa , E z = 102.4 GPa , G xz = 1.14 GPa , υ xy = 0.39 , υ xz = 0.62 .

10% of zig-zag carbon nanotubes E x = 3.70 GPa , E z = 67.21 GPa , G xz = 1.14 GPa , υ xy = 0.39 , υ xz = 0.62 .

10% of chiral carbon nanotubes E x = 3.67 GPa , E z = 126.4 GPa , G xz = 1.14 GPa , υ xy = 0.39 , υ xz = 0.62 .

The shown above values are typical for the transversally isotropic materials, and therefore they are briefly commented below.

Comment 1. The Young modulus in the direction along the symmetry axis E z exceeds essentially the Young modulus in the isotropy plane E x (from 6 to 34 times in examples above but can in some cases exceed 100 times).

Comment 2. The Lame moduli λ x and λ z repeat the relations between E x and E z .

Comment 3. The Poisson ratio υ xz along the symmetry axis Oz exceeds the classical red line in 0.5 for values of this ratio.

Comment 4. The shear moduli G xy and G xz are differed quite moderately.

3.2 The basic formulas for elastic transversely isotropic medium with axial symmetry

Let us write the basic formulas for the case of symmetry axis Oz . Then displacements are characterized by two components u r r z t , u z r z t . The motion equations in stresses have the form.

σ rr , r + σ rz , z + 1 / r σ rr σ φφ = 0 , σ rz , r + 1 / r σ φz , φ + σ zz , z + 1 / r σ rz = 0 . E43

The substitution of constitutive equations.

σ rr = C 11 u r , r + C 12 1 / r u r + C 13 u z , z , σ zz = C 13 u r , r + C 13 1 / r u r + C 33 u z , z , σ rz = 1 / 2 C 44 u z , r + u r , z , σ φz = σ = 0 E44

into the motion Eqs. (43) gives the motion equations in displacements

C 11 u r , rr + 1 / r u r , r 1 / r 2 u r + C 44 u r , zz + C 13 + C 44 u z , rz = ρ u r , tt , E45
C 44 u z , rr + 1 / r u z , r + C 33 u z , zz + C 13 + C 44 u r , rz + 1 / r u r , z = ρ u z , tt . E46

Note that Eqs. (45) and (46) include only four constants (the constant C 12 is not represented in these equations). This means that displacements and strains are described by only four constants. But the stress state is already described by all five constants.

3.3 Three classical ways of introducing the potentials in transversely isotropic elasticity

The basic equations of the theory of transversely isotropic elasticity are frequently analyzed by the use of potentials. The potentials are introduced in theory of elasticity mainly for static problems. Transition to the dynamic problems is associated with complications that are sometimes impassable. Because the problem on waves is related to the dynamic ones, let us show further the possible complications with introducing the potentials.

Way 1 [12]. It is proposed for the axisymmetric problems of equilibrium (not motion) and is based on introducing one only potential φ r z as the function of stresses. The formulas for stresses include four unknown parameters a , b , c , d , which is characteristic for representations in the transversely isotropic elasticity.

σ rr = φ , rr + b 1 / r φ , r + a φ , zz , z , σ θθ = b φ , rr + 1 / r φ , r + a φ , zz , z , E47
σ zz = c φ , rr + c 1 / r φ , r + d φ , zz , z , σ rz = φ , rr + 1 / r φ , r + a φ , zz , r . E48

The next step consists in substitution of formulas (47) and (48) into the first equation of equilibrium and the equations that are obtained from the Cauchy relations and formulas for the strain tensor. This permits to determine the unknown parameters through the elastic constants represented in the equilibrium equations. Further, the second equation of equilibrium gives the biharmonic equation for finding the potentials

Δ s 1 Δ s 2 φ = 0 , E49

where Δ sN φ = φ , rr + 1 / r φ , r + 1 / s N 2 φ , zz N = 1 ; 2 are some “complicated” copies of classical expressions Δ φ = φ , rr + 1 / r φ , r + φ , zz associated with the Laplace operator. Two constants s N are determined from the algebraic equations

s 4 a + c / d s 2 + 1 / d = 0 , s 1 , 3 = ± a + c + a + c 2 4 d 2 d , s 2 , 4 = ± a + c a + c 2 4 d 2 d . E50

Thus, a transition from the isotropic case to the transversally isotropic one complicates the procedure of solving the static problems. Here a necessity of solving the classical biharmonic equation is changed on necessity of solving some generalization of this equation in the form (49).

Way 2 [12, 16]. This way is also proposed for the static problems. Here, two potentials are introduced which are linked immediately with displacements

u r = ϕ 1 , r + ϕ 2 , r , u z = k 1 ϕ 1 , z + k 2 ϕ 2 , z . E51

A substitution of representations (51) into equations of equilibrium (45), (46)

C 11 u r , rr + 1 / r u r , r 1 / r 2 u r + C 44 u r , zz + C 13 + C 44 u z , rz = 0 ,
C 44 u z , rr + 1 / r u z , r + C 33 u z , zz + C 13 + C 44 u r , rz + 1 / r u r , z = 0

allows to determine the unknown constants k 1 , k 2 . An idea consists in that both equations must be transformed in identical equations relative to the potentials by comparing some coefficients

k 1 2 C 13 + C 44 + C 44 C 11 = kC 33 k 1 2 C 44 + C 13 + C 44 = V .

This expression gives the quadratic equation for k 1 2 and V

V 2 + C 13 2 C 44 + C 33 C 11 C 33 C 11 C 44 V + C 33 C 11 = 0 . E52

Note that the simple link V N = 1 / s N exists between constants V N and s N , which makes the ways 1 and 2 very similar. Then the potentials fulfill the equations

Δ rzN φ N = φ N , rr + 1 / r φ N , r + 1 / V N 2 φ N , zz . E53

The stresses are expressed through new potentials in such a way

σ rr = C 11 C 12 1 / r ϕ 1 , rr + ϕ 2 , rr C 44 1 + k 1 ϕ 1 , zz + 1 + k 2 ϕ 2 , zz , σ θθ = C 11 C 12 1 / r ϕ 1 , rr + ϕ 2 , rr C 13 k 1 C 12 V 1 ϕ 1 , zz + C 13 k 2 C 12 V 2 ϕ 2 , zz , σ zz = C 33 k 1 C 13 V 1 ϕ 1 , zz + C 33 k 2 C 13 V 2 ϕ 2 , zz , σ rz = C 44 1 + k 1 ϕ 1 , rz + 1 + k 2 ϕ 2 , rz . E54

Way 3 [1, 16]. This way is proposed for equations of motion, but only for the isotropic theory of elasticity. It can be used for the static problems of transversely isotropic theory of elasticity. The initial equations here are the equations of motion (43) without inertial summands

C 11 u r , rr + 1 / r u r , r 1 / r 2 u r + C 44 u r , zz + C 13 + C 44 u z , rz = 0 , E55
C 44 u z , rr + 1 / r u z , r + C 33 u z , zz + C 13 + C 44 u r , rz + 1 / r u r , z = 0 . E56

The potentials are introduced like (51), but the representations are complicated by necessity of introducing two new unknown parameters:

u r = Φ , r Ψ , z , u z = n Φ , z + m Ψ , r + m 1 / r Ψ , E57

A substitution of representations (57) into equations of motion (45) and (46) gives equations relative to the potentials. Eq. (45) gives two equations:

Φ , rr + 1 / r Φ , r + C 44 + n C 13 + C 44 C 11 Φ , zz = 0 , E58
Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ + C 44 C 11 m C 13 + C 44 Ψ , zz = 0 , E59

whereas Eq. (46) gives three equations:

Φ , rr + 1 / r Φ , r + nC 33 nC 44 + C 13 + C 44 Φ , zz = 0 , E60
Ψ , rrz + 1 / r Ψ , rz 1 / r 2 Ψ , z + C 33 m C 13 + C 44 C 44 m Ψ , zzz = 0 , E61
Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ + C 33 m C 13 + C 44 C 44 m Ψ , zz = 0 . E62

The last two equations are identical. Eqs. (58) and (60) and (59) and (62) have to be identical. This means that the coefficients in these equations have to be identical. As a result, two equations can be obtained for the determination of unknown constants n , m .

C 44 + n C 13 + C 44 C 11 = nC 33 nC 44 + C 13 + C 44 n 2 n C 11 C 33 C 44 2 C 13 + C 44 2 C 44 C 13 + C 44 + 1 = 0 , E63
n 1 , 2 = C 11 C 33 C 44 2 C 13 + C 44 2 2 C 44 C 13 + C 44 × 1 ± 1 4 C 44 C 13 + C 44 C 11 C 33 C 44 2 C 13 + C 44 2 2 , E64
C 44 C 11 m C 13 + C 44 = C 33 m C 13 + C 44 C 44 m m 2 + m C 44 2 + C 11 C 33 + C 13 + C 44 2 C 33 C 13 + C 44 + C 11 C 33 = 0 . E65
m 1 , 2 = C 44 2 + C 11 C 33 + C 13 + C 44 2 2 C 33 C 13 + C 44 × 1 ± 1 4 C 11 C 33 C 33 C 13 + C 44 C 44 2 + C 11 C 33 + C 13 + C 44 2 2 . E66

The unknown potentials Φ r z and Ψ r z have to be determined from the simple Eqs. (63) and (65) which are the classical Bessel equations of orders 0 and 1 and arguments depending on some rational combination of elastic constants.

Thus, three ways of introduction of potentials in the static problems of transversely isotropic theory of elasticity are shown. The different attempts to transfer these ways into the dynamic problems meet some troubles—the presence of inertial summands generates new additional conditions for the unknown constants in representations of potentials. Introducing the new constants does not help—the number of conditions is still more than the number of all constants.

3.4 Solving the problem on the propagation in the direction of vertical axis surface cylindrical wave for the case of transversal isotropy of medium

Consider now equations of motion (45) and (46) and introduce the potentials by the formula (57). A substitution of formula (57) into equations of motion gives five equations relative to the potentials. Eq. (57) gives two equations:

Φ , rr + 1 / r Φ , r + C 44 + n C 13 + C 44 C 11 Φ , zz = ρ C 11 Φ , tt , E67
Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ + C 44 C 11 m C 13 + C 44 Ψ , zz = ρ C 11 m C 13 + C 44 Ψ , tt . E68

Eq. (46) gives three equations:

Φ , rr + 1 / r Φ , r + nC 33 nC 44 + C 13 + C 44 Φ , zz = nC 44 + C 13 + C 44 Φ , tt , E69
Ψ , rrz + 1 / r Ψ , rz 1 / r 2 Ψ , z + C 33 m C 13 + C 44 C 44 m Ψ , zzz = ρ C 44 Ψ , ztt , E70
Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ + C 33 m C 13 + C 44 C 44 m Ψ , zz = ρ C 44 Ψ , tt . E71

Two last equations are identical. Also the equations for potential Φ must be identical as well as the equations for potential Ψ must be identical. Let us assume additionally that the problem in hand considering the solution in the form of harmonic in time cylindrical wave with unknown wave number k and known frequency ω :

Φ r z t = Φ r e i kz ωt , Ψ r z t = Ψ r e i kz ωt . E72

Note that characterization of an attenuation of wave depth down functions Φ r , Ψ r is unknown. They must be found from equations, which are obtained by substitution of representations (72) into Eqs. (67) and (71):

Φ , rr + 1 / r Φ , r C 44 + n C 13 + C 44 C 11 k 2 k L 11 2 Φ = 0 , E73

k L 11 = ω / v L 11 , v L 11 = C 11 / ρ ,

Φ , rr + 1 / r Φ , r n nC 44 + C 13 + C 44 C 33 k 2 C 11 k L 11 2 Φ = 0 , E74
Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ C 44 C 11 m C 13 + C 44 k 2 k T 44 2 Ψ = 0 , E75

k T 44 = ω / v L 44 , v L 44 = C 44 / ρ ,

Ψ , rr + 1 / r Ψ , r 1 / r 2 Ψ C 33 m C 13 + C 44 C 44 m k 2 k T 44 2 Ψ = 0 . E76

As a result, two equations can be obtained that permit to determine the constants n , m

n 2 2 N 1 n + N 2 = 0 , m 2 + 2 M 1 m + M 2 = 0 , E77
N ± M ± = N 1 M 1 ± N 1 M 1 2 N 2 M 2 , E78
2 N 1 = C 11 C 33 C 13 + C 44 2 k 2 C 11 C 11 C 44 k L 11 2 C 44 2 C 44 C 13 + C 44 k 2 , N 2 = C 44 C 11 k L 11 2 C 44 k 2 = 0 E79
2 M 1 = C 44 2 C 11 C 33 C 13 + C 44 2 k 2 C 44 2 C 11 C 44 k T 44 2 C 13 + C 44 C 33 k 2 C 44 k T 44 2 , M 2 = C 11 C 33 k 2 C 44 k T 44 2 k 2 . E80

Note that restriction on the kind of solution (it has to be a wave) allows to unite two different conditions into one—conditions for equaling coefficients in summands with the second derivative by time t and vertical coordinate z . In this case, the number of unknown constants coincides with the number of conditions which are necessary for the determination of potentials. As a result, the wave attenuation-transformed potentials can be determined from the equations of Bessel type:

Φ , rr + 1 / r Φ , r M L 11 2 Φ = 0 , M L 11 = C 44 + n C 13 + C 44 C 11 k 2 k L 11 2 , E81
Ψ , rr + 1 / r Ψ , r 1 / r 2 + M T 44 2 Ψ = 0 , M T 44 = C 44 C 11 m C 13 + C 44 k 2 k T 44 2 , E82

A success in the determination of transformed potentials is accompanied by a complication of conditions which provide the wave attenuation. They have the form.

C 44 + n C 13 + C 44 C 11 k 2 k L 11 2 > 0 , C 33 m C 13 + C 44 C 44 m k 2 k T 44 2 > 0 . E83

Let us recall that the similar conditions for the case of isotropic medium k 2 k L 2 > 0 , k 2 k T 2 > 0 are slightly simpler and coincide with the corresponding conditions of classical Rayleigh surface wave [5, 6, 7, 8, 9, 17]. A complexity of conditions (83) is increased by the complex form of dependence of constants n , m on the wave number k .

If the conditions (83) are fulfilled, then the solution of wave equations for potentials can be written in the form.

Φ r = A Φ K 0 M L 11 r , Ψ r = A Ψ K 1 M T 44 r . E84

With allowance for formulas (84), the representations of potentials becomes more definite

Φ r z t = A Φ K 0 M L 11 r e i kz ωt , Ψ r z t = A Ψ K 1 M T 44 r e i kz ωt . E85

The formula (85) completes the first analytical part of solving the problem on cylindrical surface wave.

3.5 Boundary conditions: equations for the unknown wave number

This part of analysis can be treated as the second analytical part. The boundary conditions have the form identical for all kinds of symmetry of properties. That is, they have the form (17) or (21). The formulas for stresses depend already on the symmetry of medium. The expressions for stresses through the potential reflect the features of introducing the potentials. In this case, they have the form

σ rr = λ + 2 μ Φ , rr Ψ , rz + λ 1 / r Φ , r Ψ , z + + n Φ , zz + m Ψ , rz + m 1 / r Ψ , z , E86
σ rz = μ Φ , rz Ψ , zz + n Φ , zr + m Ψ , rr + m 1 / r Ψ , r m 1 / r 2 Ψ . E87

Further, the representations (86) and (87) should be substituted into the boundary conditions, and the formulas on differentiation of Macdonald functions [3] should be taken into account:

dK 0 M L 11 rx / dr = M L 11 K 1 M L 11 rx ,
d 2 K 0 M L 11 r / dr 2 = M L 11 1 / r K 1 M L 11 r + M L 11 2 K 0 M L 11 r ,
dK 1 M T 44 r / dr = 1 / r K 1 M T 44 r M T 44 K 0 M T 44 r .

Then the boundary conditions are transformed into the algebraic equations relative to quantities K 1 M L 11 r o A Φ , K 1 M T 44 r o A Ψ

M L 11 1 r o + v L 2 v T 2 M L 11 2 v L 2 v T 2 v L 2 nk 2 K 0 M L 11 r o K 1 M L 11 r o A Φ K 1 M L 11 r o ik v L 2 v T 2 v T 2 × 2 1 m + v T 2 v L 2 v T 2 1 r o + + 1 m + v T 2 v L 2 v T 2 M T 44 K 0 M T 44 r o K 1 M T 44 r o A Ψ K 1 M T 44 r o = 0 , E88
1 + n ik K 0 M L 11 r o K 1 M L 11 r 0 K 1 M L 11 r o A Φ + m M T 44 2 + k 2 K 1 M T 44 r A Ψ = 0 . E89

When the determinant of linear homogeneous system of Eqs. (88) and (89) is equaled to zero, then the equations for the unknown wave number can be obtained:

1 + n k 2 v L 2 v T 2 v T 2 K 0 M L 11 r o K 1 M L 11 r 0 2 1 m + v T 2 v L 2 v T 2 1 / r o + 1 m + v T 2 v L 2 v T 2 M T 44 K 0 M T 44 r o K 1 M T 44 r o m M T 44 2 + k 2 M L 11 1 r o + v L 2 v T 2 M L 11 2 v L 2 v T 2 v L 2 nk 2 K 0 M L 11 r o K 1 M L 11 r o = 0 . E90

Note that the sufficiently complex expression relative to the wave number is hidden coefficients M L 11 , M T 44 of Macdonald’s functions K 0 M L 11 r o K 1 M L 11 r 0 , K 0 M T 44 r o K 1 M T 44 r o . Therefore, the analytical part of analysis is finished on these formulas. Further, the numerical approaches have to be utilized.

Note also that the simple and convenient condition from analysis of classical surface Rayleigh wave [6, 7, 8, 9, 10, 17], when the wave number depends only on ratio v L 2 / v T 2 , does not exist in the analysis of cylindrical surface wave. Here, the parameters M L 11 , M T 44 depend on the complicated form on all elastic constants. Of course, the Macdonald functions can be represented approximately through their arguments. But only the numerical methods can give the final result—the value of wave number or phase velocity.

Advertisement

4. Solving the problem on propagating in the direction of symmetry axis surface wave within the framework of linearized theory of elasticity with allowance for initial stresses

Note that analysis of cylindrical surface wave in isotropic medium was first carried out by Biot [1] in 1952 and the transversally isotropic medium with initial stresses was first carried out by Guz et al. in 1974 [18].

Let us show below an analysis of the problem in hand that is carried out in Subchapter “Longitudinal Waves” of Chapter 4 “Waves in Cylindrical Media” of volume 2 of edition [19]. Here, the cylinder of circular cross-section is considered, and the longitudinal wave is defined as the wave propagating in the direction of cylinder axis Oy 3 . The problem is assumed to be axisymmetric and is described within the framework of linearized theory of elasticity for bodies with initial stresses. The cylindrical coordinates r ' θ y 3 are introduced, and displacements are taken in the form

u r ' = u r ' r ' y 3 t , u θ = 0 , u y 3 = u y 3 r ' y 3 t E91

The medium is assumed isotropic or transversally isotropic. The main relations for transversal isotropy are described by independent constant

ω 1111 , ω 1122 , ω 1133 , ω 1221 , ω 1313 , ω 1331 , ω 3113 , ω 3333 . E92

Note that as shown in (92), eight constants are necessary in the linearized theory, but in the framework of linear theory, they have the form (33), and their number is five.

Further, the general solutions of basic equations in displacements are utilized. These equations have the form (3.174) [19]

ω lmαβ 2 u α / x k x β = ρδ 2 u α / τ 2 E93

where only eight independent constants (92) must be taken into account.

The corresponding equations of linear theory of elasticity for the case of transversally isotropic medium without of initial stresses are written above as Eqs. (45) and (46).

The general solutions for the case of axial symmetry are expressed through one potential in the form (4.13) [19]

u r ' = 2 / r y 3 X , E94
u 3 = ω ' 1111 + ω ' 1313 1 ω ' 1111 Δ ' 1 + ω ' 3113 2 / y 3 2 ρ ' 2 / τ 2 X ' , Δ ' = 2 / r ' 2 + 1 / r ' / r ' .

Note that in Section 3 of this chapter, two potentials Φ , Ψ are introduced by formula (57), which corresponds and generalizes the procedure used in Biot’s analysis [1].

The longitudinal harmonic wave is described analytically through the potential in the form (101) [19]

X ' r ' y 3 τ = X ' 1 r ' e i ky 3 ωτ , E95

where the unknown amplitude X ' 1 r ' has to be determined by substitution of solution (4.13) [19] into the second Eq. (3.362) [19] (for potential X ' ). This gives Eq. (4.16) [19]:

ω ' 1111 ω ' 1331 Δ ' 1 k 2 ξ ' 2 2 Δ ' 1 k 2 ξ ' 3 2 k 2 ρ ' C cp 2 ω ' 1111 + ω ' 1331 Δ ' 1 k 2 ω ' 1111 + ω ' 3113 + ρ ' 2 C cp 2 X ' 1 = 0 , E96
C cp = ω / k , ξ 2 , 3 ' 2 = c ' ± c ' 2 ω ' 3333 ω ' 3113 / ω ' 1111 ω ' 1331 , c ' = 1 / 2 ω ' 3333 / ω ' 1331 + ω ' 3113 / ω ' 1111 ω ' 1111 + ω ' 1331 2 / ω ' 1111 ω ' 1331 ,

which further is written in the form

Δ ' 1 ζ 2 ' 2 Δ ' 1 ζ 3 ' 2 = 0 E97

The unknown quantities ζ ' 2 , 3 must be found from the linear algebraic equation of the fourth degree (4.20) [19].

ω ' 1111 ω ' 1331 ζ ' 4 + k 4 ρ ' C cp 2 ω ' 3333 ρ ' C cp 2 ω ' 3113 + k 2 ω ' 1111 ρ ' C cp 2 ω ' 3333 + ω ' 1331 ρ ' C cp 2 ω ' 3113 + ω ' 1111 + ω ' 3113 2 ζ ' 2 = 0 , E98

The solution (95) describes the surface wave, if amplitude X ' 1 r ' attenuates with increasing the radius. This is provided by the condition that quantities ζ ' 2 , 3 is unequal and pure imaginary. Then the potential gains the form (4.22) [19].

X ' 1 r ' = B 10 J 0 ζ ' 2 r ' + B 20 K 0 ζ ' 2 r ' + B 30 J 0 ζ ' 3 r ' + B 40 K 0 ζ ' 3 r ' , E99

The shown part of analysis from introducing the potential by formula (94) to representation of solution by formula (99) inclusive can be compared with analogous part of analysis from Section 3 of this chapter (from introducing the potentials by formula (57) to the solution in the form of (85)). It is easy to see a difference in representations (99) and (85): formula (99) uses the Bessel functions and in particular the Macdonald function of zero index, whereas formula (85) uses (like the Biot’s solution (14)) the Macdonald functions (16) of the zero and first indexes.

The next part of analysis of cylindrical wave consists in substitution of solution into boundary conditions of the form (99) [19]

Q ' r ' r ' = 0 , Q ' r ' 3 = 0 when r ' = R ' 1 , R ' 2 . E100

The case of oscillatory behavior of wave in the direction of radius is considered with pointing that the case of surface wave is the same type. A substitution of solution (99) into conditions (4.79) [19] gives the dependence of velocity of surface wave or its wave number on frequency—a dispersion equation in the form of determinant of the fourth order in the form (4.26) [19].

det α ij Δ ω k = 0 ; i , j = 1 , 2 , 3 , 4 . E101

This finishes the analytical part of analysis shown in [19]. It corresponds to the part of Section 3.5 of this chapter, where the explicit form of dispersive equations is proposed in the form (90) that includes the Macdonald functions of the zero and first orders which represent some generalization of dispersion Eq. (27) obtained by Biot.

Advertisement

5. Conclusions

This chapter proposes three fragments of analytical analysis of the cylindrical surface wave propagating in the vertical direction of circular cylindrical cavity. The first fragment shows the analytical part of pioneer work of Biot. It represents the classicism of mathematical procedures and physical comments of Biot. Properly speaking, the clear and understandable Rayleigh’s scheme is saved, but it is complemented by some findings reflecting the features of cylindrical waves. Two next fragments show the more late development of the Biot’s problem. They are different by influence of the Biot’s procedure. The approach shown in Section 3 is more close to the Biot’s analytical scheme, whereas Section 4 proposes as an independent scheme that is more close to the Rayleigh scheme. Nevertheless, all fragments testify the mathematical complexity in solving the problem on the cylindrical surface waves. Thus, revisiting the old Biot’s problem shows that it still generates new scientific and practical problems.

References

  1. 1. Biot MA. Propagation of elastic waves in a cylindrical bore containing a fluid. Journal of Applied Physics. 1952;23(9):997-1005. DOI: 10.1063/1.1702365
  2. 2. Achenbach JD. Wave Propagation in Elastic Solids. Amsterdam: North-Holland; 1973. 425 p
  3. 3. FWJ O, Lozier DW, Bousvert RF, Clark CW, editors. NIST (National Institute of Standards and Technology). Handbook of Mathematical Functions. Cambridge: Cambridge University Press; 2010. 968 p
  4. 4. Rushchitsky JJ. Nonlinear Elastic Waves in Materials. Series: Foundations of Engineering Mechanics. Heidelberg: Springer; 2014. 455 p. DOI: 10.1007/978-3-319-00464-8
  5. 5. Sedov LI. A Course in Continuum Mechanics. Vol. I-IV. Amsterdam: Wolters Noordhoff Publishing; 1971
  6. 6. Viktorov IA. Rayleigh and Lamb Waves. NY: Plenum Press; 1967. 168 p
  7. 7. Nowacki W. Theoria sprazystosci. Warszawa: PWN; 1970. 769 p (in Polish)
  8. 8. Rushchitsky JJ, Tsurpal SI. Waves in Materials with Microstructure. Kiev: S.P. Timoshenko Institute of Mechanics; 1998. 377 p (in Ukrainian)
  9. 9. Fedorov FI. Theory of Elastic Waves in Crystals. NY: Academic Press; 1975. 388 p
  10. 10. Dieulesaint E, Royer D. Ondes elastiques dans les solides. Application au traitement du signal. Paris: Masson et Cie; 1974. 424 p (in French)
  11. 11. Royer D, Dieulesaint E. Elastic Waves in Solids. Vols. I and II. Advanced Texts in Physics. Berlin: Springer; 2000
  12. 12. Lekhnitsky SG. Theory of Elasticity of Anisotropic Elastic Body. San Francisco: Golden Day Inc; 1963. 404 p
  13. 13. Kiselev AP. Rayleigh wave with a transverse structure. Proceedings of the Royal Society of London. Series A. 2004;460(2050):3059-3064. DOI: 10.1098/rspa.2004.1353
  14. 14. Guz AN, Rushchitsky JJ. Short Introduction to Mechanics of Nanocomposites. Rosemead, CA: Scientific & Academic Publishing; 2013. 280 p
  15. 15. Cheng M, Chen W. Weerasoorlya Mechanical properties of Kevlar ® KM2. Journal of Engineering Materials and Technology. 2005;127(2):197-203. DOI: 10.1115/1.1857937
  16. 16. Elliot HA. Three-dimensional stress distribution in hexagonal aelotropic crystals. Mathematical Proceedings of the Cambridge Philosophical Society. 1948;44(4):522-533. DOI: 10.1017/50305004100024531
  17. 17. Rushchitsky JJ. Theory of Waves in Materials. Copenhagen: Ventus Publishing ApS; 2011. 270 p. DOI: 10.13140/RG.2.1.3162.8647
  18. 18. Guz AN, Kushnir VP, Makhort FG. On propagation of waves in cylinders with initial stresses. Izvestiya, Academy of Sciences, USSR. Seriya Mekhanika Tverdogo Tiela. 1974;5:67-74 (In Russian)
  19. 19. Guz AN. Elastic Waves in Bodies with Initial Stresses. 2 Vols. Naukova Dumka: Kiev; 1986. 376 and 536 p (in Russian)

Written By

Jeremiah Rushchitsky

Submitted: 22 February 2019 Reviewed: 17 May 2019 Published: 27 September 2019