Skip to main content

The effect of initial stress and magnetic field on wave propagation in human dry bones

Abstract

The aim of the present paper is to study the influence of initial stress and magnetic field on the propagation of harmonic waves in a human long dry bone as transversely isotropic material, subject to the boundary conditions that the outer and inner surfaces are traction free. The equations of elastodynamics are solved in terms of displacements. The natural frequency of plane vibrations in the case of harmonic vibrations has been obtained. The frequencies and phase velocity are calculated numerically, the effects of initial stress and magnetic field are discussed. Comparisons are made with the result in the absence of initial stress and magnetic field.

MSC:74B05.

1 Introduction

The investigation of wave propagation over a continuous medium has very important application in the fields of engineering, medicine and in bioengineering. Application of the poroelastic materials in medicinal fields such as cardiovascular, dental and orthopedics is well known. The dry bone is piezoelectric in the classical sense [1, 2], i.e., mechanical stress results in electric polarization (the indirect effect); and an applied electric field causes strain (the converse effect). Since that time, many others have confirmed the capacity of bones to produce piezoelectric potentials [3]. Electrical properties of bone are relevant not only as a hypothesized feedback mechanism for bone adaptation and remodeling, but also in the context of external electrical stimulation of bone in order to aid its healing and repair [4]. In orthopedics, the propagation of wave over bone is used in monitoring the rate of fracture healing. There are two types of osseous tissue such as trabecular or cancellous and cortical or compact bone, which are of different materials with respect to their mechanical behavior. In macroscopic terms, the porosity percentage in the compact bone is 3-5%, whereas in the cancellous or trabecular the porosity percentage is up to 90% [1].

Mahmoud [1, 5, 6] investigated the wave propagation under the effects of initial stress, rotation and magnetic field in cylindrical poroelastic bones, a granular medium and a porous medium. Theoretical analyses of bone piezoelectricity may be relevant to the issue of bone remodeling. Recent thorough studies have explored electromechanical effects in wet and dry bone [7, 8]. They suggest that two different mechanisms are responsible for these effects: classical piezoelectricity mainly due to the molecular asymmetry of collagen in dry bone and streaming potentials found in moist or living bone and generated by the flow of a liquid across charged surfaces. The second mechanism was argued by dielectric measurements, and it was suggested that the electromechanical effect in wet (fluid saturated) bone is not due to a piezoelectric effect [9]. Abd-Alla and Mahmoud [10, 11] solved a magneto-thermoelastic problem in a rotating non-homogeneous orthotropic hollow cylinder under the hyperbolic heat conduction model and investigated analytical solution of wave propagation in non-homogeneous orthotropic rotating elastic media. Abd-Alla et al. [12] studied the propagation of S-wave in a non-homogeneous anisotropic incompressible and initially stressed medium under the influence of gravity field. Honarvarla et al. [13], Ding et al. [14] studied the elasticity of transversely isotropic materials. Chen et al. [15, 16] investigated the free vibration and general solution of non-homogeneous transversely isotropic magneto-electroelastic hollow cylinders. Abd-Alla et al. [17, 18] studied the problem of transient coupled thermoelasticity of an annular fin and the problem of radial vibrations in a non-homogeneous isotropic cylinder under the influence of initial stress and magnetic field. Mofakhami et al. [19] studied the finite cylinder vibrations with different end boundary conditions. Abd-Alla et al. [20, 21] studied the effect of rotation, magnetic field and initial stress on peristaltic motion of micropolar fluid and investigated the effect of rotation on a non-homogeneous infinite cylinder of orthotropic material.

In this paper, the equations of elastodynamics for transversely isotropic material under the effect of initial stress and magnetic field are solved in terms of displacement potentials. Also, this paper is concerned with the determination of phase velocity and the eigenvalues of natural frequency of plane vibrations of bones under the effect of initial stress and magnetic field for different boundary conditions in the cases of harmonic vibrations. The numerical results of the frequency equation are discussed in detail for transversely isotropic material and the effect of initial stress and magnetic field for different cases is indicated by figures.

2 Formulation of the problem

Consider a homogeneous and transversely isotropic long bone as a hollow cylinder of inner radius a and outer radius b taking the cylindrical polar coordinates (r,θ,z) such that the z-axis points vertically upward along the bone axis.

The equations of the elastodynamic medium in the presence of magnetic field are as follows:

τ r r r + r 1 τ r θ θ + τ r z z + r 1 ( τ r r τ θ θ )+ ( J × B ) r =ρ 2 u r t 2 ,
(1a)
τ r θ r + r 1 τ θ θ θ + τ θ z z +2 r 1 τ θ r + ( J × B ) θ =ρ 2 u θ t 2 ,
(1b)
τ r z r + r 1 τ θ z θ + τ z z z + r 1 τ θ r + ( J × B ) z =ρ 2 u z t 2 ,
(1c)

where ρ is the density of bones, τ r r , τ θ θ , τ z z and τ r z are the stresses, u r , u θ and u z are the displacement components and B = μ e H , where μ e is the magnetic permeability, H =(0, H 0 ,0) and H 0 is the intensity of the uniform magnetic field parallel to θ-axes, F is Lorentz’s body forces vector where F can take the following form:

F =( f r , f θ , f z )= ( ( J × B ) r , ( J × B ) θ , ( J × B ) z ) .
(2)

The relations of stresses-displacement for homogeneous transversely isotropic bone in two dimensions are in the form:

τ r r = ( c 11 + P ) u r r + ( c 12 + P ) u r r + ( c 13 + P ) u z z ,
(3a)
τ θ θ = ( c 12 + P ) u r r + ( c 11 + P ) u r r + ( c 13 + P ) u z z ,
(3b)
τ z z = c 13 u r r + c 13 u r r + c 33 u z z ,
(3c)
τ r z = c 44 ( u z r + u r z ) ,
(3d)

where c 11 ; c 12 ; c 22 ; c 13 ; c 23 ; c 44 are the elastic constants of bone, c 11 = c 22 , c 13 = c 23 , and P is the initial stress compression. The electromagnetic field is governed by Maxwell’s equations considering that the medium is a perfect electric conductor taking into account the absence of the displacement current SI.

curl h = J , curl E = μ e h t , div h = 0 , div E = 0 , h = curl ( u × H 0 ) , H = H 0 + h , H 0 = ( 0 , H 0 , 0 ) .
(4)

Maxwell’s stresses τ ¯ i j take the following form:

σ ¯ i j = μ e [ H i h j + H j h i ( H h ) δ i j ] ,i,j=1,2,3.
(5)

Two-dimentional equations of elastodynamics are as follows:

σ r r r + τ r z z + 1 r ( τ r r τ θ θ )+ H 0 2 μ e r 1 [ u z z + u r r + r 2 u z r z + r 2 u r r 2 ] =ρ 2 u r t 2 ,
(6a)
τ r z r + σ z z z + 1 r τ r z + H 0 2 μ e ( 2 u z z 2 + 2 u r r z ) =ρ 2 u z t 2 .
(6b)

Substituting equations (3a)-(3d) into equations (6a) and (6b), we obtain

2 u r r 2 + 1 r u r r 1 r 2 u r + c 44 c 11 2 u r z 2 + c 13 + c 44 c 11 2 u z r z + 1 2 P ( u r z u z ) z + H 0 2 μ e r 1 [ u z z + u r r + r 2 u z r z + r 2 u r r 2 ] = ρ c 11 ( 2 u r t 2 ) ,
(7a)
2 u z z 2 + c 44 c 33 ( 2 u z r 2 + 1 r u z r ) + c 13 + c 44 c 33 ( 2 u r r z + 1 r u r z ) + P 2 r ( r ( u r z u z ) ) z + H 0 2 μ e ( 2 u z z 2 + 2 u r r z ) = ρ c 33 ( 2 u z t 2 ) .
(7b)

System (7a)-(7b) can be simplified as follows:

1 2 r 2 [ 2 c 11 u r + r 2 H 0 2 μ e u z z + 2 ( c 11 + H 0 2 μ e ) u r r + r 2 ρ 2 u r t 2 + ( c 44 + P ) 2 u r z 2 + ( 2 c 13 + c 44 P + 2 H 0 2 μ e ) 2 u z r z + 2 ( c 11 + H 0 2 μ e ) 2 u r r 2 ] = 0 ,
(8a)
1 2 r [ ( 2 c 13 + c 44 + P ) u r z + ( c 44 P ) u z r + r ( 2 ρ 2 u z t 2 + 2 ( c 33 + H 0 2 μ e ) 2 u z z 2 + ( 2 c 13 + c 44 + P + 2 H 0 2 μ e ) 2 u r r z + ( c 44 P ) 2 u z r 2 ) ] = 0 .
(8b)

3 Solution of the problem

By Helmohltz’s theorem [22] the displacement vector u can be written as

u = ̲ ϕ 1 + ̲ Ψ ,
(9)

where the two functions ϕ 1 and Ψ =(0, ψ 1 ,0) are known in the theory of elasticity, by Lame’s potentials irrotational and rotational parts of the displacement vector u , respectively. The cylinder being bounded by the curved surface, the stress distribution includes the effect of both ϕ 1 and ψ . It is possible to take only one component of the vector ψ to be nonzero as follows: ψ =(0, ψ 1 ,0). From (9) we obtain

u r = ϕ 1 r ψ 1 z ,
(10a)
u z = ψ 1 r + ϕ z + ψ 1 r .
(10b)

Substituting equations (10a) and (10b) into equations (8a) and (8b), we get two independent equations for ϕ 1 and ψ 1 as follows:

1 2 r 2 [ 2 c 11 ( ψ 1 z ϕ 1 r ) + r ( 2 H 0 2 μ e ( ψ 1 z r 1 + 2 ϕ 1 z 2 + 2 ψ 1 r z ) + 2 ( c 11 + H 0 2 μ e ) ( 2 ψ 1 r z + 2 ϕ 1 r 2 ) + r ( 2 ρ ( 3 ψ 1 z t 2 3 ϕ 1 r t 2 ) + ( c 44 + P ) ( 3 ψ 1 z 3 + 3 ϕ 1 r z 2 ) + ( 2 c 13 + c 44 P + 2 H 0 2 μ e ) × ( r 2 ψ 1 z + r 1 2 ψ 1 r z + 3 ϕ 1 r z 2 + 3 ψ 1 r 2 z ) + 2 ( c 11 + H 0 2 μ e ) ( 3 ψ 1 r 2 z + 3 ϕ 1 r 3 ) ) ) ] ,
(11a)
1 2 r ( 2 c 13 + c 44 + P ) ( 2 ψ 1 z 2 + 2 ϕ 1 r z ) + ( c 44 P ) ( ψ 1 r 2 + ψ 1 r r 2 + 2 ϕ 1 r z + 2 ψ 1 r 2 ) + r ( 2 ρ ( 2 ψ 1 t 2 r 1 + 3 ϕ 1 z t 2 + 3 ψ 1 r t 2 ) + 2 ( c 33 + H 0 2 μ e ) × ( 2 ψ 1 z 2 r 1 + 3 ϕ 1 z 3 + 3 ψ 1 r z 2 ) + ( 2 c 13 + c 44 + P + 2 H 0 2 μ e ) ( 3 ψ 1 r z 2 + 3 ϕ 1 r 2 z ) + ( c 44 P ) ( ( 2 ψ 1 2 r ψ 1 r + r 2 2 ψ 1 r 2 ) r 3 + 3 ϕ 1 r 2 z + 3 ψ 1 r 3 ) ) = 0 .
(11b)

To study the propagation of harmonic waves in the z-direction, we assume a solution in the form

Φ 1 (r,z,t)= Φ 2 (r) e i ( γ z ω 5 t ) ,
(12a)
Ψ 1 (r,z,t)= Ψ 2 (r) e i ( γ z ω 5 t ) ,
(12b)

where γ is the wave number, ω is the angular frequency. Substituting equations (12a) and (12b) into equations (11a) and (11b) and omitting the factor exponential throughout, we have

1 2 r 2 i e i ( z γ t ω 5 ) [ 2 i H 0 2 r γ 2 μ e Φ 2 ( r ) γ ( 2 c 11 2 c 13 + P + c 44 ( 1 + r 2 γ 2 ) + r 2 ( P γ 2 2 ρ ω 5 2 ) ) Ψ 2 ( r ) 2 i c 11 2 i c 13 r 2 γ 2 d Φ 2 d r 2 i c 44 r 2 γ d Φ 2 d r 2 2 i H 0 2 r 2 γ 2 μ e d Φ 2 d r + 2 i r 2 ρ ω 5 2 d Φ 2 d r + 2 c 11 r γ d Ψ 2 d r 2 c 13 r γ d Ψ 2 d r c 14 r γ d Ψ 2 d r + P r γ d Ψ 2 d r 2 H 0 2 r γ μ e d Ψ 2 d r + 2 i c 11 r d 2 Φ 2 d r 2 + 2 i H 0 2 r μ e d 2 Φ 2 d r 2 + 2 c 11 r 2 γ d 2 Ψ 2 d r 2 2 c 13 r 2 γ d 2 Ψ 2 d r 2 c 44 r 2 γ d 2 Ψ 2 d r 2 + P r 2 γ d 2 Ψ 2 d r 2 + 2 i r 2 ( c 11 + H 0 2 μ e ) d 3 Φ 2 d r 3 ] = 0 ,
(13a)
1 2 r 3 e i ( z γ t ω 5 ) [ 2 i r 3 γ ( γ 2 ( c 33 + H 0 2 μ e ) ρ ω 5 2 ) Φ 2 ( r ) + ( c 44 + c 44 r 2 γ 2 + P ( 1 + r 2 γ 2 ) + 2 r 2 ( γ 2 ( c 13 c 33 H 0 2 μ e ) + ρ ω 5 2 ) ) Ψ 2 ( r ) + r ( 2 i ( c 13 + c 44 ) r γ d Φ 2 d r + ( P + P r 2 γ 2 + c 44 ( 1 + r 2 γ 2 ) + 2 r 2 ( ( c 13 c 33 ) γ 2 + ρ ω 5 2 ) ) d Ψ 2 d r + r ( 2 i r γ ( c 13 + c 44 + H 0 2 μ e ) d 2 Φ 2 d r 2 + ( c 44 P ) ( 2 d Ψ 2 d r + r d 3 Ψ 2 d r 3 ) ) ) ] = 0 .
(13b)

Similar results were obtained by Elnagar and Abd-Alla [23], the former deriving the constitutive equation for Rayleigh waves in an elastic medium under initial stress, and the latter deriving the constitutive equation for thermoelastic problems in an infinite cylinder under initial stress. From equations (13a) and (13b), we get equations (14a) and (14b) as follows:

2 ( H 0 2 r γ 2 μ e Φ 2 ( r ) ( c 11 + r 2 ( γ 2 ( c 13 + c 44 + H 0 2 μ e ) ρ ω 5 2 ) ) d Φ 2 d r + r ( c 11 + H 0 2 μ e ) ( d 2 Φ 2 d r 2 + r d 3 Φ 2 d r 3 ) ) = 0 ,
(14a)
( c 44 + c 44 r 2 γ 2 + P ( 1 + r 2 γ 2 ) + 2 r 2 ( γ 2 ( c 13 c 33 H 0 2 μ e ) + ρ ω 5 2 ) ) Ψ 2 ( r ) + r ( ( P + P r 2 γ 2 + c 44 ( 1 + r 2 γ 2 ) + 2 r 2 ( ( c 13 c 33 ) γ 2 + ρ ω 5 2 ) ) d Ψ 2 d r + ( c 44 P ) r ( 2 d 2 Ψ 2 d r 2 + r d 3 Ψ 2 d r 3 ) ) = 0 .
(14b)

Equation (14a) represents the shear wave, and equation (14b) represents the longitudinal wave. The solution of equations (14a) and (14b) can be written in the following form:

Ψ 2 = A 1 r d 1 J n 1 ( d 1 r)+ B 1 r d 1 Y n 1 ( d 2 r),
(15a)
Φ 2 = A 2 J n 2 ( e 1 r)+ B 2 Y n 2 ( e 2 r),
(15b)

where

n 1 = d 0 4 c 11 c 44 + 4 c 13 c 44 + c 44 2 + 4 c 11 P 4 c 13 P 2 c 44 P + P 2 + H 0 4 μ e 2 2 c 11 2 c 13 c 44 + P , d 1 = H 0 2 μ e 2 c 11 2 c 13 c 44 + P , d 2 = i r c 44 γ 2 + P γ 2 2 ρ ω 5 2 2 c 11 2 c 13 c 44 + P , e 1 = H 0 2 μ e 2 ( c 13 + c 44 + H 0 2 μ e ) , d 0 = 4 c 11 2 8 c 11 c 13 + 4 c 13 2 , n 2 = H 0 2 μ e 2 ( c 13 + c 44 + H 0 2 μ e ) , e 2 = c 33 γ 2 + H 0 2 γ 2 μ e ρ ω 5 2 c 13 + c 44 + H 0 2 μ e .

A 1 , A 2 , B 1 and B 2 are arbitrary constants, J is the Bessel function of the first kind and of order zero, Y is the Bessel function of the second kind and of order zero. J n is the Bessel function of the first kind and of order n, Y n is the Bessel function of the second kind and of order n. From equations (12a), (12b) and (15a), (15b) we get

Ψ 1 (r,z,t)= e i ( z γ t ω 5 ) ( A 1 r d 1 J n 1 ( d 1 r ) + B 1 r d 2 Y n 1 ( d 2 r ) ) ,
(16a)
Φ 1 (r,z,t)= e i ( z γ t ω 5 ) ( A 2 J n 2 ( e 1 r ) + B 2 Y n 2 ( e 2 r ) ) .
(16b)

Substituting equations (16a) and (16b) into equations (10a) and (10b), we obtain the final solution of displacement components in the following form:

u r = 1 r e i ( z γ t ω 5 ) [ A 2 e 1 J 1 + n 2 ( e 1 r ) J 1 + n 2 ( e 1 r ) 2 i r d 1 γ A 1 J n 1 ( d 2 r ) + B 1 Y n 1 ( d 2 r ) u r = + B 2 e 2 Y 1 + n 2 ( e 2 r ) Y 1 + n 2 ( e 2 r ) ] ,
(17a)
u z = 1 r e i ( z γ t ω 5 ) [ A 1 d 2 r 1 + d 1 J 1 + n 1 ( d 2 r ) + A 1 ( 1 + d 1 n 1 ) r d 1 J n 1 ( d 2 r ) + i A 2 r γ J n 2 ( e 1 r ) u z = + B 1 d 2 r 1 + d 1 Y 1 + n 1 ( d 2 r ) + B 1 ( 1 + d 1 n 1 ) r d 1 Y n 1 ( d 2 r ) + i B 2 r γ Y n 2 ( e 2 r ) ] .
(17b)

Substituting equations (17a) and (17b) into equations (3a)-(3d), we obtain the final solution of the stress components of solid in the following form:

τ r r = 1 4 r e i ( z γ t ω 5 ) [ A 2 c 11 e 1 2 r J 2 + n 2 ( e 1 r ) + 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 r ) 2 A 2 c 11 e 1 2 r J n 2 ( e 1 r ) 4 A 2 c 13 r γ 2 J n 2 ( e 1 r ) 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 r ) + A 2 c 11 e 1 2 r J 2 + n 2 ( e 1 r ) 4 i r d 1 γ ( A 1 ( c 11 c 13 ) d 2 r J 1 + n 1 ( d 2 r ) + A 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 r ) + B 1 ( c 11 c 13 ) d 2 r Y 1 + n 1 ( d 2 r ) + B 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 r ) ) + B 2 c 11 e 2 2 r Y 2 + n 2 ( e 2 r ) + 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 r ) 2 B 2 c 11 e 2 2 r Y n 2 ( e 2 r ) 4 B 2 c 13 r γ 2 Y n 2 ( e 2 r ) 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 r ) + B 2 c 11 e 2 2 r Y 2 + n 2 ( e 2 r ) ] ,
(18a)
τ θ θ = 1 4 r e i ( z γ t ω 5 ) [ A 2 c 12 e 1 2 r J 2 + n 2 ( e 1 r ) + 2 A 2 c 11 e 1 J 1 + n 2 ( e 1 r ) 2 A 2 c 12 e 1 2 r J n 2 ( e 1 r ) 4 A 2 c 13 r γ 2 J n 2 ( e 1 r ) 2 A 2 c 11 e 1 J 1 + n 2 ( e 1 r ) + A 2 c 12 e 1 2 r J 2 + n 2 ( e 1 r ) 4 i r d 1 γ ( A 1 ( c 12 c 13 ) d 2 r J 1 + n 1 ( d 2 r ) + A 1 ( c 11 + c 12 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 r ) + B 1 ( c 12 c 13 ) d 2 r Y 1 + n 1 ( d 2 r ) + B 1 ( c 11 + c 12 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 r ) ) + B 2 c 12 e 2 2 r Y 2 + n 2 ( e 2 r ) + 2 B 2 c 11 e 2 Y 1 + n 2 ( e 2 r ) 2 B 2 c 12 e 2 2 r Y n 2 ( e 2 r ) 4 B 2 c 13 r γ 2 Y n 2 ( e 2 r ) 2 B 2 c 11 e 2 J 1 + n 2 ( e 2 r ) + B 2 c 12 e 2 2 r Y 2 + n 2 ( e 2 r ) ] ,
(18b)
τ z z = 1 r e i ( z γ t ω 5 ) [ i ( c 13 c 33 ) r d 1 γ ( A 1 d 2 r J 1 + n 1 ( d 2 r ) + A 1 ( 1 + d 1 n 1 ) J n 1 ( d 2 r ) + B 1 d 2 r Y 1 + n 1 ( d 2 r ) + B 1 ( 1 + d 1 n 1 ) Y n 1 ( d 2 r ) ) + r 1 ( A 2 ( c 13 ( n 2 e 1 r ) ( n 2 + e 1 r ) c 33 r 2 γ 2 ) J n 2 ( e 1 r ) + B 2 ( c 13 ( n 2 e 2 r ) ( n 2 + e 2 r ) c 33 r 2 γ 2 ) Y n 2 ( e 2 r ) ) ] ,
(18c)
τ r z = 1 2 r 2 c 44 e i ( z γ t ω 5 ) [ A 1 r d 1 ( 1 + ( d 1 n 1 ) 2 + r 2 ( d 2 2 + γ 2 ) ) J 2 + n 1 ( d 2 r ) + 2 A 1 r 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 r 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + r 2 ( d 2 2 + γ 2 ) ) ) J 1 + n 1 ( d 2 r ) d 2 + 2 i A 2 e 1 r 2 γ J 1 + n 2 ( e 1 r ) 2 i A 2 n 2 r γ J n 2 ( e 1 r ) B 1 r d 1 ( 1 + ( d 1 n 1 ) 2 + r 2 ( d 2 2 + γ 2 ) ) Y 2 + n 1 ( d 2 r ) + 2 B 1 r 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 r 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + r 2 ( d 2 2 + γ 2 ) ) ) Y 1 + n 1 ( d 2 r ) d 2 + 2 i B 2 e 2 r 2 γ Y 1 + n 2 ( e 2 r ) 2 i B 2 n 2 r γ Y n 2 ( e 2 r ) ] ,
(18d)

where J 1 and Y 1 are the Bessel functions of the first order. In the following section, solutions of hollow circular cylinders with three different boundary conditions are performed.

4 Boundary conditions and frequency equation

In this case, we are going to obtain the frequency equation for the boundary conditions. Plane vibrations cylindrical bone free surface traction. In this case, we have

τ r r (r)+ σ r r (r)=0, τ r z (r)=0,at r=a,
(19a)
τ r r (r)+ σ r r (r)=0, τ r z (r)=0,at r=b,
(19b)

which correspond to free inner and outer surfaces, respectively. From equations (18a)-(18d), (5) and (19a)-(19b) we obtain four homogeneous linear equations in A 1 , B 1 , A 2 and  B 2

A 2 c 11 e 1 2 r J 2 + n 2 ( e 1 a ) + 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 a ) 2 A 2 c 11 e 1 2 a J n 2 ( e 1 a ) 4 A 2 c 13 a γ 2 J n 2 ( e 1 a ) 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 a ) + A 2 c 11 e 1 2 a J 2 + n 2 ( e 1 a ) 4 i a d 1 γ ( A 1 ( c 11 c 13 ) d 2 a J 1 + n 1 ( d 2 a ) + A 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 a ) + B 1 ( c 11 c 13 ) d 2 a Y 1 + n 1 ( d 2 a ) + B 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 a ) ) + B 2 c 11 e 2 2 a Y 2 + n 2 ( e 2 a ) + 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 a ) 2 B 2 c 11 e 2 2 a Y n 2 ( e 2 a ) 4 B 2 c 13 a γ 2 Y n 2 ( e 2 a ) 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 a ) + B 2 c 11 e 2 2 a Y 2 + n 2 ( e 2 a ) = 0 ,
(20a)
A 1 a d 1 ( 1 + ( d 1 n 1 ) 2 + a 2 ( d 2 2 + γ 2 ) ) J 2 + n 1 ( d 2 a ) + 2 A 1 a 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 a 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + a 2 ( d 2 2 + γ 2 ) ) ) J 1 + n 1 ( d 2 a ) d 2 + 2 i A 2 e 1 a 2 γ J 1 + n 2 ( e 1 a ) 2 i A 2 n 2 a γ J n 2 ( e 1 a ) B 1 a d 1 ( 1 + ( d 1 n 1 ) 2 + a 2 ( d 2 2 + γ 2 ) ) Y 2 + n 1 ( d 2 a ) + 2 B 1 a 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 a 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + a 2 ( d 2 2 + γ 2 ) ) ) Y 1 + n 1 ( d 2 a ) d 2 + 2 i B 2 e 2 a 2 γ Y 1 + n 2 ( e 2 a ) 2 i B 2 n 2 a γ Y n 2 ( e 2 a ) = 0 ,
(20b)
A 2 c 11 e 1 2 b J 2 + n 2 ( e 1 b ) + 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 b ) 2 A 2 c 11 e 1 2 b J n 2 ( e 1 b ) 4 A 2 c 13 b γ 2 J n 2 ( e 1 b ) 2 A 2 c 12 e 1 J 1 + n 2 ( e 1 b ) + A 2 c 11 e 1 2 b J 2 + n 2 ( e 1 b ) 4 i b d 1 γ ( A 1 ( c 11 c 13 ) d 2 b J 1 + n 1 ( d 2 b ) + A 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 b ) + B 1 ( c 11 c 13 ) d 2 b Y 1 + n 1 ( d 2 b ) + B 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 b ) ) + B 2 c 11 e 2 2 b Y 2 + n 2 ( e 2 b ) + 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 b ) 2 B 2 c 11 e 2 2 b Y n 2 ( e 2 b ) 4 B 2 c 13 b γ 2 Y n 2 ( e 2 b ) 2 B 2 c 12 e 2 Y 1 + n 2 ( e 2 b ) + B 2 c 11 e 2 2 b Y 2 + n 2 ( e 2 b ) = 0 ,
(20c)
A 1 b d 1 ( 1 + ( d 1 n 1 ) 2 + b 2 ( d 2 2 + γ 2 ) ) J 2 + n 1 ( d 2 b ) + 2 A 1 b 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 b 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + b 2 ( d 2 2 + γ 2 ) ) ) J 1 + n 1 ( d 2 b ) d 2 + 2 i A 2 e 1 b 2 γ J 1 + n 2 ( e 1 b ) 2 i A 2 n 2 b γ J n 2 ( e 1 b ) B 1 b d 1 ( 1 + ( d 1 n 1 ) 2 + b 2 ( d 2 2 + γ 2 ) ) Y 2 + n 1 ( d 2 b ) + 2 B 1 b 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 b 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + b 2 ( d 2 2 + γ 2 ) ) ) Y 1 + n 1 ( d 2 b ) d 2 + 2 i B 2 e 2 b 2 γ Y 1 + n 2 ( e 2 b ) 2 i B 2 n 2 b γ Y n 2 ( e 2 b ) = 0 .
(20d)

These are a set of four homogeneous algebraic equations involving four unknown integration constants A 1 , B 1 , A 2 and B 2 . For a nontrivial solution of these equations, the determinant of the coefficient matrix must vanish. The zero determinant of the coefficient matrix will give the frequency equation for the surface waves. Thus, elimination of these unknowns would give us the frequency equation as follows:

Δ=| B 11 B 12 B 13 B 14 B 21 B 22 B 23 B 24 B 31 B 32 B 33 B 34 B 41 B 42 B 43 B 44 |=0.
(21)

The roots of equation (21) give the values of natural frequency for the free surfaces of the cylinder.

B 11 = 4 i a d 1 γ ( ( c 11 c 13 ) d 2 a J 1 + n 1 ( d 2 a ) + ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 a ) ) , B 12 = c 11 e 1 2 r J 2 + n 2 ( e 1 a ) + 2 c 12 e 1 J 1 + n 2 ( e 1 a ) 2 c 11 e 1 2 a J n 2 ( e 1 a ) 4 c 13 a γ 2 J n 2 ( e 1 a ) 2 c 12 e 1 J 1 + n 2 ( e 1 a ) + c 11 e 1 2 a J 2 + n 2 ( e 1 a ) , B 13 = 4 i a d 1 γ ( + B 1 ( c 11 c 13 ) d 2 a Y 1 + n 1 ( d 2 a ) + B 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 a ) ) , B 14 = c 11 e 2 2 a Y 2 + n 2 ( e 2 a ) + 2 c 12 e 2 Y 1 + n 2 ( e 2 a ) 2 c 11 e 2 2 a Y n 2 ( e 2 a ) 4 c 13 a γ 2 Y n 2 ( e 2 a ) 2 c 12 e 2 Y 1 + n 2 ( e 2 a ) + c 11 e 2 2 a Y 2 + n 2 ( e 2 a ) , B 21 = a d 1 ( 1 + ( d 1 n 1 ) 2 + a 2 ( d 2 2 + γ 2 ) ) J 2 + n 1 ( d 2 a ) + 2 a 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 a 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + a 2 ( d 2 2 + γ 2 ) ) ) J 1 + n 1 ( d 2 a ) d 2 , B 22 = 2 i e 1 a 2 γ J 1 + n 2 ( e 1 a ) 2 i n 2 a γ J n 2 ( e 1 a ) , B 23 = a d 1 ( 1 + ( d 1 n 1 ) 2 + a 2 ( d 2 2 + γ 2 ) ) Y 2 + n 1 ( d 2 a ) + 2 a 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 a 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + a 2 ( d 2 2 + γ 2 ) ) ) Y 1 + n 1 ( d 2 a ) d 2 , B 24 = 2 i e 2 a 2 γ Y 1 + n 2 ( e 2 a ) 2 i n 2 a γ Y n 2 ( e 2 a ) , B 31 = 4 i b d 1 γ ( ( c 11 c 13 ) d 2 b J 1 + n 1 ( d 2 b ) + ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) J n 1 ( d 2 b ) ) , B 32 = c 11 e 1 2 r J 2 + n 2 ( e 1 b ) + 2 c 12 e 1 J 1 + n 2 ( e 1 b ) 2 c 11 e 1 2 b J n 2 ( e 1 b ) 4 c 13 b γ 2 J n 2 ( e 1 b ) 2 c 12 e 1 J 1 + n 2 ( e 1 b ) + c 11 e 1 2 b J 2 + n 2 ( e 1 b ) , B 33 = 4 i b d 1 γ ( + B 1 ( c 11 c 13 ) d 2 b Y 1 + n 1 ( d 2 b ) + B 1 ( c 12 + c 11 ( d 1 n 1 ) + c 13 ( 1 d 1 + n 1 ) ) Y n 1 ( d 2 b ) ) , B 34 = c 11 e 2 2 b Y 2 + n 2 ( e 2 b ) + 2 c 12 e 2 Y 1 + n 2 ( e 2 b ) 2 c 11 e 2 2 b Y n 2 ( e 2 b ) 4 c 13 b γ 2 Y n 2 ( e 2 b ) 2 c 12 e 2 Y 1 + n 2 ( e 2 b ) + c 11 e 2 2 b Y 2 + n 2 ( e 2 b ) , B 41 = b d 1 ( 1 + ( d 1 n 1 ) 2 + b 2 ( d 2 2 + γ 2 ) ) J 2 + n 1 ( d 2 b ) + 2 b 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 b 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + b 2 ( d 2 2 + γ 2 ) ) ) J 1 + n 1 ( d 2 b ) d 2 , B 42 = 2 i e 1 b 2 γ J 1 + n 2 ( e 1 b ) 2 i n 2 b γ J n 2 ( e 1 b ) , B 43 = b d 1 ( 1 + ( d 1 n 1 ) 2 + b 2 ( d 2 2 + γ 2 ) ) Y 2 + n 1 ( d 2 b ) + 2 b 1 + d 1 ( d 1 2 ( 1 + n 1 ) + d 1 ( 2 ( 1 + n 1 ) n 1 + d 2 2 b 2 ) + ( 1 + n 1 ) ( 1 + n 1 2 + b 2 ( d 2 2 + γ 2 ) ) ) Y 1 + n 1 ( d 2 b ) d 2 , B 44 = 2 i e 2 b 2 γ Y 1 + n 2 ( e 2 b ) 2 i n 2 b γ Y n 2 ( e 2 b ) .

5 Numerical results and discussion

For the numerical calculation of dimensionless frequency and phase velocity under the effect of initial stress P and magnetic field H 0 , one shall investigate the frequency equations given by (21) numerically for a particular model. Since these equations are an implicit function relation of dimensionless frequency, one can proceed with finding the variation of frequency with ratio h. Once the frequency has been computed, the corresponding effect of initial stress P and magnetic field H 0 on the frequency equations for dimensionless frequency (the eigenvalues) can be studied by taking values of ratio h. As an illustrative purpose of the foregoing solutions, the cylinder has the following geometric and material constants which are used in the computations given by [1, 2, 14]: c 11 =2.33× 10 8 , c 12 =0.85× 10 8 , c 13 =1.3× 10 8 , c 44 =0.62× 10 8 , c 33 =4.981×1 0 8 , ρ=3.986×1 0 3 .

For various values of dimensionless frequency, the phase velocity is obtained from frequency equation flexural modes (n=1). Here, we explain graphically our results of the previous applications. The dimensionless frequency and phase velocity were calculated with the aid of an electronic computer by using the half-interval method. Figure 1 shows the variation of phase velocity with respect to initial stress P for different values of magnetic field H 0 for inner and outer free surfaces. The phase velocity decreases with the increase in initial stress, but the phase velocity increases when the values of magnetic field H 0 increase. Figure 2 shows variation of the dimensionless frequency with respect to initial stress P for different values of magnetic field H 0 for inner and outer free surfaces. The dimensionless frequency increases with the increase in initial stress, and the dimensionless frequency increases when the values of magnetic field H 0 increase. Figure 3 shows variation of phase velocity with respect to the ratio (h) of inner and outer free surfaces for different values of initial stress P when magnetic field H 0 =0.8× 10 3 . The phase velocity decreases with the increase in the ratio (h) of inner and outer free surfaces, and the phase velocity increases when the values of initial stress increase. Figure 4 shows variation of non-dimensional frequency with respect to the ratio (h) of inner and outer free surfaces for different values of initial stress P when magnetic field H 0 =0.8× 10 3 . The dimensionless frequency increases with the increase in the ratio (h) of inner and outer free surfaces, and the dimensionless frequency increases when the values of initial stress P decrease. It should be noticed that the effects of initial stress P and magnetic field H 0 on the dimensionless frequency tend to be the increasing dimensionless frequency. These results are specific for the case considered, but other cases may have different trends because of the dependences of the results on the mechanical and constants of the material [2426]. The results indicate that the effect of initial stress P and magnetic field H 0 is very pronouced.

Figure 1
figure 1

Variation of phase velocity with respect to initial stress P for different values of magnetic field H 0 for inner and outer free surfaces.

Figure 2
figure 2

Variation of the dimensionless frequency with respect to initial stress P for different values of magnetic field H 0 for inner and outer free surfaces.

Figure 3
figure 3

Variation of phase velocity with respect to the ratio (h) of inner and outer free surfaces for different values of initial stress P when magnetic field H 0 =0.8× 10 3 .

Figure 4
figure 4

Variation of non-dimensional frequency with respect to the ratio (h) of inner and outer free surfaces for different values of initial stress P when magnetic field H 0 =0.8× 10 3 .

6 Conclusion

This study has presented the effect of initial stress P and magnetic field H 0 on surface wave dispersion in bone. The phase velocity and the dimensionless frequency for this problem are obtained from the dimensionless frequency equation. A numerical method has been presented for obtaining the estimates of phase velocity and dimensionless frequencies of vibration of transversely isotropic bone using the half-interval method. The eigenvalues are calculated for different cases and compared with those reported in the absence of initial stress P and magnetic field H 0 . The effects of initial stress P and magnetic field H 0 on the dimensionless frequencies and the phase velocity were indicated by figures.

References

  1. Mahmoud SR: Wave propagation in cylindrical poroelastic dry bones. Appl. Math. Inf. Sci. 2010, 4(2):209-226.

    MathSciNet  Google Scholar 

  2. Abd-Alla AM, Mahmoud SR, Abo-Dahab SM: Wave propagation modeling in cylindrical human long wet bones with cavity. Meccanica 2011, 46(6):1413-1428. 10.1007/s11012-010-9398-5

    Article  MathSciNet  Google Scholar 

  3. Ramtani S: Electro-mechanics of bone remodeling. Int. J. Eng. Sci. 2008, 46: 1173-1182. 10.1016/j.ijengsci.2008.06.001

    Article  MathSciNet  Google Scholar 

  4. Ramtani S, Zidi M: A theoretical model of the effect of continuum damage on a bone adaptation model. J. Biomech. 2001, 34: 471-479. 10.1016/S0021-9290(00)00215-3

    Article  Google Scholar 

  5. Mahmoud SR: Influence of rotation and generalized magneto-thermoelastic on Rayleigh waves in a granular medium under effect of initial stress and gravity field. Meccanica 2012, 47(7):1561-1579. 10.1007/s11012-011-9535-9

    Article  MathSciNet  Google Scholar 

  6. Mahmoud SR: Effect of rotation and magnetic field through porous medium on peristaltic transport of a Jeffrey fluid in tube. Math. Probl. Eng. 2011., 2011: Article ID 971456

    Google Scholar 

  7. Nickerson DP, Smith NP, Hunter PJ: A model of cardiac cellular electromechanics. Philos. Trans. R. Soc. A, Math. Phys. Eng. Sci. 2001, 359(1783):1159-1172. 10.1098/rsta.2001.0823

    Article  Google Scholar 

  8. Kohl P, Sachs F: Mechanoelectric feedback in cardiac cells. Philos. Trans. R. Soc. A, Math. Phys. Eng. Sci. 2001, 359(783):1173-1185.

    Article  Google Scholar 

  9. Eringen AC: Electromagnetic theory of microstretch elasticity and bone modeling. Int. J. Eng. Sci. 2004, 42: 231-242. 10.1016/S0020-7225(03)00288-X

    Article  MathSciNet  Google Scholar 

  10. Abd-Alla AM, Mahmoud SR: Magneto-thermoelastic problem in rotating non-homogeneous orthotropic hollow cylindrical under the hyperbolic heat conduction model. Meccanica 2010, 45(4):451-462. 10.1007/s11012-009-9261-8

    Article  MathSciNet  Google Scholar 

  11. Abd-Alla AM, Mahmoud SR: Analytical solution of wave propagation in non-homogeneous orthotropic rotating elastic media. J. Mech. Sci. Technol. 2012, 26(3):917-926. 10.1007/s12206-011-1241-y

    Article  Google Scholar 

  12. Abd-Alla AM, Mahmoud SR, Abo-Dahab SM, Helmi MIR: Propagation of S -wave in a non-homogeneous anisotropic incompressible and initially stressed medium under influence of gravity field. Appl. Math. Comput. 2011, 217(9):4321-4332. 10.1016/j.amc.2010.10.029

    Article  MathSciNet  Google Scholar 

  13. Honarvarla F, Enjilela E, Sinclair A, Mirnezami S: Wave propagation in transversely isotropic cylinders. Int. J. Solids Struct. 2007, 44: 5236-5246. 10.1016/j.ijsolstr.2006.12.029

    Article  Google Scholar 

  14. Ding HJ, Chen WQ, Zhang L: Elasticity of Transversely Isotropic Materials. Springer, Berlin; 2006.

    Google Scholar 

  15. Chen WQ, Lee KY, Ding HJ: On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic hollow cylinders. J. Sound Vib. 2005, 279: 237-251. 10.1016/j.jsv.2003.10.033

    Article  Google Scholar 

  16. Chen WQ, Lee KY, Ding HJ: General solution for transversely isotropic magneto-electro-thermo-elasticity and potential theory method. Int. J. Eng. Sci. 2004, 42: 1361-1379. 10.1016/j.ijengsci.2004.04.002

    Article  Google Scholar 

  17. Abd-Alla AM, Mahmoud SR, Abo-Dahab SM: On problem of transient coupled thermoelasticity of an annular fin. Meccanica 2012, 47(5):1295-1306. 10.1007/s11012-011-9513-2

    Article  MathSciNet  Google Scholar 

  18. Abd-Alla AM, Mahmoud SR: On problem of radial vibrations in non-homogeneity isotropic cylinder under influence of initial stress and magnetic field. J. Vib. Control 2013, 77(12):269-276.

    MathSciNet  Google Scholar 

  19. Mofakhami MR, Toudeshky HH, Hashemi SH: Finite cylinder vibrations with different end boundary conditions. J. Sound Vib. 2006, 297: 293-314. 10.1016/j.jsv.2006.03.041

    Article  Google Scholar 

  20. Abd-Alla AM, Yahya GA, Mahmoud SR, Alosaimi HS: Effect of the rotation, magnetic field and initial stress on peristaltic motion of micropolar fluid. Meccanica 2012, 47(6):1455-1465. 10.1007/s11012-011-9528-8

    Article  MathSciNet  Google Scholar 

  21. Abd-Alla AM, Mahmoud SR, AL-Shehri NA: Effect of the rotation on a non-homogeneous infinite cylinder of orthotropic material. Appl. Math. Comput. 2011, 217(22):8914-8922. 10.1016/j.amc.2011.03.077

    Article  Google Scholar 

  22. Mahmoud SR, Abd-Alla AM, AL-Shehri NA: Effect of the rotation on plane vibrations in a transversely isotropic infinite hollow cylinder. Int. J. Mod. Phys. B 2011, 25(26):3513-3528. 10.1142/S0217979211100928

    Article  Google Scholar 

  23. El-Naggar AM, Abd-Alla AM: On a generalized thermo-elastic problems in an infinite cylinder under initial stress. Earth Moon Planets 1987, 37: 213-223. 10.1007/BF00116637

    Article  Google Scholar 

  24. Marin M, Agarwal RP, Mahmoud SR: Nonsimple material problems addressed by the Lagrange’s identity. Bound. Value Probl. 2013., 2013: Article ID 135

    Google Scholar 

  25. Marin M, Agarwal RP, Mahmoud SR: Modeling a microstretch thermoelastic body with two temperatures. Abstr. Appl. Anal. 2013., 2013: Article ID 583464 10.1155/2013/583464

    Google Scholar 

  26. Lu H, Sun L, Sun J: Existence of positive solutions to a non-positive elastic beam equation with both ends fixed. Bound. Value Probl. 2012., 2012: Article ID 56

    Google Scholar 

Download references

Acknowledgements

This article (project) was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant No. (130-075-D1434). The authors, therefore, acknowledge with thanks DSR technical and financial support.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Samy R Mahmoud.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors, SRM, AT, ATA and KSA contributed to each part of this work equally and read and approved the final version of the manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mahmoud, S.R., Tounsi, A., Ali, A.T. et al. The effect of initial stress and magnetic field on wave propagation in human dry bones. Bound Value Probl 2014, 135 (2014). https://doi.org/10.1186/1687-2770-2014-135

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-135

Keywords