Next Article in Journal
Variable Structure PID Controller for Satellite Attitude Control Considering Actuator Failure
Previous Article in Journal
Numerical Analysis of the Calcaneal Nail C-NAIL
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Modeling of Acoustic Emission Due to an Internal Point Source in a Transversely Isotropic Cylinder

Rm #306, Integrity Diagnostics Korea (IDK), IT Venture Town, 35, Techno 9-ro, Yuseong-gu, Daejeon 34027, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(10), 5272; https://doi.org/10.3390/app12105272
Submission received: 5 April 2022 / Revised: 18 May 2022 / Accepted: 18 May 2022 / Published: 23 May 2022
(This article belongs to the Section Acoustics and Vibrations)

Abstract

:
In this paper, the displacement fields responsible for acoustic emission (AE), excited from a point source in a transversely isotropic cylinder, are derived by solving the Navier-Lamé (NL) equation. The point source as an internal defect is represented by a spatiotemporal concentrated force. The introduction of three potentials correlated with the point source to displacement field vector decouples the coupled NL equation in cylindrical coordinates. Under these conditions, we solve the radial, tangential, and axial displacement fields. Analytical simulations of AE were carried out at several point source locations. Our results demonstrate that analytical modeling is a powerful tool for characterizing AE features generated from an internal defect source.

Graphical Abstract

1. Introduction

Acoustic emission (AE) refers to the elastic wave generated by a material as a result of a sudden release of energy (other than heat) from localized sources within a solid. Material degradation related to deformation and fracture development is classified as the primary AE source, and is usually distinguished from secondary AE sources including leak and flow, chemical reactions, and the fabrication process. Crack formation and growth are the most important AE sources in practical non-destructive testing (NDT) investigations. The formation of a new crack face is accompanied by sudden changes in stress and displacement of the material in the vicinity of the crack. Consequently, an elastic wave is generated from the tip of the crack source [1]. AE testing has several unique advantages over other NDT methods, the most significant of which is its real-time monitoring capability; this allows for the detection of damage and degradation in various materials and structures [2,3,4,5,6]. AE data, refined with appropriate algorithms [7,8,9,10,11], provide useful information about the source of the emission and progression of crack growth; however, a quantitative interpretation of the AE signal is required to understand the physical processes underlying AE features. Analytical modeling of the AE signal is very important; however, in cylindrical structures, it is difficult to predict the AE generated from an internal crack.
The point source, defined as a body force f acting at a point, has been adapted as an AE excitation source in seismic displacements, crack formation and fracture, and concentrated vertical step force [12,13,14,15,16]. In the elastic field, the point source has been treated as a concentrated time harmonic source [17]. Although the AE generated by the point source is important for characterizing real signals observed in practical NDT investigations, theoretical modeling has been limited to spherical geometries with an infinite domain, given that a homogeneous solution is obtainable from the three-dimensional wave equation in spherical coordinates [12,13,17]. Most theoretical works on elastic wave propagation in cylindrical coordinates have focused on situations with or without external perturbations [18,19,20,21,22]. In cylindrical structures, the displacement field has been formulated in terms of one compression (P) potential, and two shear (vertically polarized, SV, and horizontally polarized, SH) potentials, by using the models proposed by Morse and Feshbach [23] and Buchwald [24], respectively. The basic difference between the two models is that the compression and shear parts are separated in one model, but not the other. The two models have been examined comprehensively by Honarvar et al. [18,19] and Sakhr et al. [21,22]. In contrast to cases that are unperturbed or perturbed by external forces, the three potentials should be incorporated along with the concentrated force (CF) exerted by the internal defects [17,25].
In this study, a “concentrated force-incorporated potential” (CFIP) is introduced into the displacement field on the basis of the model proposed by Morse and Feshbach [23], and the Navier-Lamé (NL) equation is constructed (which involves the CF in a transversely isotropic cylinder [TIC]). The CFIP allows the NL equation to be separated into three partial differential equations (PDEs), representing the potentials for P, SV, and SH waves. As an internal defect, the CF has both spatial and temporal properties, represented by the delta function and harmonic oscillation of the point source, respectively. For solving the NL equation, determining Green’s function for the delta function in cylindrical coordinates is the first task. Mohammad et al. proposed a Green’s function for a closed cylinder, using “the method of separation of variables” [26]. We reconstructed Green’s function for the delta function by applying the continuity and discontinuity principles to the boundary around the point source. Complete formation of the CF allowed us to solve the three potentials generated by the point source in the closed cylinder. Absolute values of the displacement fields were calculated by applying a fundamental set of boundary conditions in the cylinder to the solutions, in which two modes of CF along the radial and axial directions were considered. To our knowledge, no theoretical work on AE, excited by an internal point source in cylindrical geometries, has been presented in the literature. This paper establishes a mathematical model that provides insight into the overall process of the AE signal from generation and propagation to reception. These studies can be employed for evaluating the AE signal generated from an internal crack in a cylindrical structure.

2. Green’s Function

The body force arising from the point source located at x 0 is mathematically formulated in terms of an oscillating impulse with natural frequencies of the material in a given geometry, as follows:
f = P ( t ) δ ( x x 0 ) e i ω t ,
where P(t) is the CF vector, δ ( x x 0 ) is the delta function, and ω is the predominant angular frequency ( ω = 2 π ν ) of AE. Green’s function g ( x ; x 0 ) , as the solution of the delta function, is defined as
𝛻 2 g ( x ; x 0 ) = δ ( x x 0 ) .
The force can be rewritten as
f = 𝛻 2 P g ( x ; x 0 ) = 𝛻 [ 𝛻 ·   P g ( x ; x 0 ) ] 𝛻 × [ 𝛻 × P g ( x ; x 0 ) ] .
In cylindrical coordinates, Equation (2) is expressed as
𝛻 2 g ( r , θ , z ; r 0 , θ 0 , z 0 ) = δ ( r r 0 ) δ ( θ θ 0 ) δ ( z z 0 ) r .  
If r r 0 0 , θ θ 0 0 , and z z 0 0 ,
𝛻 2 g ( x ; x 0 ) = 0 .
Green’s function is separable in a cylindrical coordinate system as
g ( x ; x 0 ) = g r ( r ; r 0 ) g θ ( θ ; θ 0 ) g z ( z ; z 0 ) .
Substituting the Laplacian in cylindrical coordinates into Equation (5) gives the following:
r 2 g r 2 g r r 2 + r g r g r r + 1 g θ 2 g θ θ 2 + r 2 g z 2 g z z 2 = 0 .
Letting
1 g z 2 g z z 2 = κ z 2 ,
1 g θ 2 g θ θ 2 = v 2 ,
then Equation (7) has the form
r 2 2 g r r 2 + r g r r + ( κ z 2 r 2 v 2 ) g r = 0 .  
The solution of Equation (8) can be formulated in two regions:
g z ( z ; z 0 ) = { A z ( z 0 ) e κ z z                       0 < z < z 0 < l B z ( z 0 ) e κ z z                     0 < z 0 < z < l
The continuity and discontinuity principles are applied to find coefficients A z and B z . The continuity principle states that Green’s function is continuous around the point source, i.e.,
g z ( z ; z 0 ) | z = z 0 + = g z ( z ; z 0 ) | z = z 0 .
Equation (12) leads to
A z ( z 0 ) e κ z z 0 = B z ( z 0 ) e κ z z 0 .
Although Green’s function is discontinuous at the point source, the difference in Green’s function between the adjacent two points, z 0 ± ε , around the point source is unity. When ε becomes 0, the discontinuity principle leads to
g z ( z ; z 0 ) z | z = z 0 + g z ( z ; z 0 ) z | z = z 0 = 1 .
As a result, we obtain
κ z B z ( z 0 ) e κ z z 0 κ z A z ( z 0 ) e κ z z 0 = 1 .
From Equations (13) and (15),
A z ( z 0 ) = 1 2 κ z e κ z z 0       and       B z ( z 0 ) = 1 2 κ z e κ z z 0 .
Equation (11) becomes
g z ( z ; z 0 ) = { 1 2 κ z e κ z ( z 0 z )                       0 < z < z 0 < l 1 2 κ z e κ z ( z z 0 )                     0 < z 0 < z < l .
The solution of Equation (9) can be written as
g θ ( θ ; θ 0 ) = { A θ ( θ 0 ) cos v θ + B θ ( θ 0 )   sin v θ                   0 < θ < θ 0 < 2 π C θ ( θ 0 ) cos v θ + D θ ( θ 0 )   sin v θ                     0 < θ 0 < θ < 2 π .
Applying the continuity principle to Equation (17) results in
A θ ( θ 0 ) = C θ ( θ 0 )       and       B θ ( θ 0 ) = D θ ( θ 0 ) .
The second restriction for g θ ( θ ; θ 0 ) is the symmetry between the point source and observation point, i.e.,
g θ ( θ 0 ; θ ) = g θ ( θ ; θ 0 ) .
Equation (19) leads to
A θ ( θ ) cos v θ 0 + B θ ( θ )   sin v θ 0 = A θ ( θ 0 ) cos v θ + B θ ( θ 0 )   sin v θ .
To satisfy the conservation law,
A θ ( θ 0 ) = cos v θ 0   and   B θ ( θ 0 ) = sin v θ 0 .
Equations (18) and (20) result in
g θ ( θ ; θ 0 ) = [ cos v θ 0 cos v θ + sin v θ 0 sin v θ ] = cos [ v ( θ θ 0 ) ] .
For a complete cylinder, g θ ( θ ; θ 0 ) is periodic, with a period of 2π:
cos [ v ( 2 π + θ θ 0 ) ] = cos [ v ( θ θ 0 ) ] .
Therefore,
v = 0 , ± 1 , ± 2 ,   .
The final form of Equation (17) is
g θ ( θ ; θ 0 ) = cos [ v ( θ θ 0 ) ]             ( v = 0 , ± 1 , ± 2 ,   ) .
The solution of Equation (10) is a typical Bessel function, given by
g r ( r ; r 0 ) = { A v r ( r 0 ) J v ( κ z r )                       0 < r < r 0 < a B v r ( r 0 ) J v ( κ z r )                       0 < r 0 < r < a
where J v is the first kind of Bessel function of v-th order. Note that the second kind of Bessel function is excluded because it has a singularity at the origin, which is included in the domain of the cylinder. Applying the continuity principle to Equation (22), i.e., g r ( r ; r 0 ) | r = r 0 + = g r ( r ; r 0 ) | r = r 0 , results in
A v r ( r 0 ) J v ( κ z r 0 ) = B v r ( r 0 ) J v ( κ z r 0 ) .
We find
A v r ( r 0 ) = B v r ( r 0 ) ,
where r 0 = x 0 i 2 + x 0 j 2 . In addition, applying the symmetry principle to Equation (22) gives g r ( r 0 ; r ) = g r ( r ; r 0 ) , with
A v r ( r ) J v ( κ z r 0 ) = A v r ( r 0 ) J v ( κ z r ) .
we obtain
A v r ( r 0 ) = J v ( κ z r 0 ) .
Substituting Equations (23) and (24) into Equation (22) gives
g r ( r ; r 0 ) = J v ( κ z r 0 ) J v ( κ z r )     ( 0 < r < a ) .
Applying the boundary condition, i.e.,
g r ( r ; r 0 ) | r = a r 0 = 0 ,                 g r ( r ; r 0 ) | r = r 0 0 ,
to Equation (25) gives J v [ κ z ( a r 0 ) ] = 0 . Denoting the n-th root of the first kind of Bessel function of the v-th order as r v n , i.e., J v ( r v n ) = 0 , then r v n = κ z ( a r 0 )   ( n = 1 ,   2 ,   ) . Since the first root of the first kind of Bessel function is most significant, we select n = 1.
κ z = r v 1 a r 0 .
Introducing a parameter, A v 1 , to Equation (25) gives
g r ( r ; r 0 ) = A v 1   J v ( r v 1 a r 0 r 0 ) J v ( r v 1 a r 0 r )         ( 0 < r < a ) .
From Equations (16), (21) and (27), the formula of the Green’s function for the Kronecker delta function becomes
g ( r , θ , z ; r 0 , θ 0 , z 0 ) = A v 1   J v ( r v 1 a r 0 r 0 ) J v ( r v 1 a r 0 r ) cos [ v ( θ θ 0 ) ] × { a r 0 2 r v 1 e r v 1 a r 0 ( z 0 z )           0 < z < z 0 < l a r 0 2 r v 1 e r v 1 a r 0 ( z z 0 )           0 < z 0 < z < l .  
The next task is to find the constant A v 1 in Equation (28). Equation (4) is rewritten as
𝛻 2 g ( r , θ , z ; r 0 , θ 0 , z 0 ) = 𝛻 2 g r ( r ;   r 0 ) g θ ( θ ; θ 0 ) g z ( z ; z 0 ) = g θ g z r r ( r g r r ) + g r g z r 2 2 g θ θ 2 + g r g θ 2 g z z 2 = δ ( r r 0 ) δ ( θ θ 0 ) δ ( z z 0 ) r .
First, solve δ ( z z 0 ) by integrating both sides over (0, l)
z = 0 z = l [ g θ g z r r ( r g r r ) + g r g z r 2 2 g θ θ 2 + g r g θ 2 g z z 2 ] d z = g r g θ z = 0 z = l 2 g z z 2 d z .
Therefore,
g r g θ g z z | 0 z = l = δ ( r r 0 ) δ ( θ θ 0 ) r .
From Equations (16), (26), and (30), we obtain
g r g θ g z z | 0 z = l = g r g θ [ 1 2 e r v 1 a r 0 ( z 0 z ) | 0 z 0 + 1 2 e r v 1 a r 0 ( z z 0 ) | z 0 l ] = g r g θ 2 [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] ,
and
g r g θ 2 [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] = δ ( r r 0 ) δ ( θ θ 0 ) r .  
Substituting Equation (21) into the above equation and integrating both sides over (0, 2π) gives
g r 2 [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] v = 0 2 π cos [ v ( θ θ 0 ) ] d θ = δ ( r r 0 ) r 0 2 π δ ( θ θ 0 ) d θ .
Since
0 2 π cos [ v ( θ θ 0 ) ] d θ = { 2 π                   if   θ = θ 0 0               otherwise ,
the above equation is given by
π g r [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] = δ ( r r 0 ) r .
Substituting Equation (27) into Equation (32) results in
π [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] A v 1   J v ( r v 1 a r 0 r 0 ) J v ( r v 1 a r 0 r ) = δ ( r r 0 ) r .
Multiplying both sides of Equation (33) by r J p ( r v 1 a r 0 r ) and integrating over (0, a r 0 ) gives
π [ e r v 1 a r 0 z 0 + e r v 1 a r 0 ( l z 0 ) 2 ] × A v 1   J v ( r v 1 a r 0 r 0 ) 0 a r 0 J v ( r v 1 a r 0 r ) J p ( r v 1 a r 0 r ) r d r = 0 a δ ( r r 0 ) J p ( r v 1 a r 0 r )   d .
Given [26]
0 a r 0 J μ ( α μ a r 0 r ) J ν ( α ν a r 0 r ) r d r = ( a r 0 ) 2 2 [ J μ + 1 ( α μ ) ] 2 δ μ ν
and
0 a r 0 δ ( r r 0 ) J μ ( r ) d r = J μ ( r 0 ) ,
Equation (34) can be rewritten as
π [ e r v 1 z 0 / ( a r 0 ) + e r v 1 ( l z 0 ) / ( a r 0 ) 2 ] A v 1   J v ( r v 1 a r 0 r 0 ) ( a r 0 ) 2 2 [ J v + 1 ( r v 1 ) ] 2 = J v ( r v 1 a r 0 r 0 ) .
From Equation (35), we obtain the constant
A v 1 = 2 π ( a r 0 ) 2 [ e r v 1 z 0 / ( a r 0 ) + e r v 1 ( l z 0 ) / ( a r 0 ) 2 ] [ J v + 1 ( r v 1 ) ] 2 .
Finally, we obtain the Green’s function for the Kronecker delta function
g ( r , θ , z ; r 0 , θ 0 , z 0 ) = G v 1   J v ( r v 1 a r 0 r 0 ) J v ( r v 1 a r 0 r ) { cos [ v ( θ θ 0 ) ] } × { e r v 1 ( z 0 z ) a r 0           ( 0 < z < z 0 < l ) e r v 1 ( z z 0 ) a r 0           ( 0 < z 0 < z < l ) ,
where
G v 1 = a r 0 2 r v 1 A v 1 = 1 π ( a r 0 ) r v 1 [ e r v 1 z 0 / ( a r 0 ) + e r v 1 ( l z 0 ) / ( a r 0 ) 2 ] [ J v + 1 ( r v 1 ) ] 2 .
For an azimuthally independent force (v = 0), Equation (37) becomes
g ( r , z ; r 0 , z 0 ) = G 1   J 0 ( r 01 a r 0 r ) ,
where
G 1 = J 0 ( r 01 a r 0 r 0 ) π ( a r 0 ) r 01 [ e r v 1 z 0 / ( a r 0 ) + e r v 1 ( l z 0 ) / ( a r 0 ) 2 ] [ J 1 ( r 01 ) ] 2   × { e r 01 ( z 0 z ) a r 0           ( 0 < z < z 0 < l ) e r 01 ( z z 0 ) a r 0           ( 0 < z 0 < z < l ) .
For convenience, we introduce new coordinates, defined as ξ = r r 0 , ϑ = θ θ 0 and η = z z 0 , where the location of the point source is the new origin at ξ 0 = 0 , ϑ 0 = 0 and η 0 = 0 . Thus, Equations (39) and (40) are rewritten as
g ( ξ , ϑ , η ) = G 1 J 0 ( r 01 a r 0 ξ ) ,
G 1 = J 0 ( r 01 a r 0 ξ 0 ) π ( a r 0 ) r 01 [ e r v 1 z 0 / ( a r 0 ) + e r v 1 ( l z 0 ) / ( a r 0 ) 2 ] [ J 1 ( r 01 ) ] 2 × { e r 01 η a r 0 (   z 0 l < η < 0 ) e r 01 η a r 0 ( 0 < η < l z 0 ) .

3. Displacement Fields Generated by a Point Source

The NL equation governing a wave in an elastic and homogeneous medium subject to a local body force f can be written in vector form, as follows [27]
( λ + 2 μ ) 𝛻 ( 𝛻 · u ) μ 𝛻 × ( 𝛻 × u ) + f = ρ 2 u t 2 ,
where u is the displacement vector, λ and μ are the Lamé constants, and ρ is the density of the media. In this paper, the displacement field u in cylindrical coordinates proposed by Morse and Feshbach [23] is given as
u = 𝛻 Φ + 𝛻 × ( Χ e ^ z ) + a 𝛻 × 𝛻 × ( Ψ e ^ z ) ,
where Φ is a scalar potential for the compressional wave (P), Χ e ^ z is a vector potential for the SH wave, Ψ e ^ z is a vector potential for the SV wave, and a is the radius of the cylinder. When the displacement field is generated by an intrinsic point defect, the three potentials are correlated with the force vector P (referred to as CFIP) as
Φ = 𝛻 · P ϕ ,
Χ e ^ z = 𝛻 × P χ ,
Ψ e ^ z = 𝛻 × P ψ ,
where ϕ, χ and ψ are unknown scalar functions, and the negative sign for the SH potential is used for further applications. Factoring the spatial and temporal parts in the potentials leads to
ϕ ( x , t ;   x 0 ) = ϕ ( x ; x 0 ) e i ω t ,
χ ( x , t ; x 0 ) = χ ( x ; x 0 ) e i ω t ,
ψ ( x , t ; x 0 ) = ψ ( x ; x 0 ) e i ω t ,
The displacement vector can be written as:
u ( x , t ; x 0 , t 0 ) = u ( r , θ , z , t ; r 0 , θ 0 , z 0 , t 0 ) = [ 𝛻 ( 𝛻 · P ϕ ) 𝛻 × ( 𝛻 × P χ ) a 𝛻 × 𝛻 × ( 𝛻 × P ψ ) ] e i ω t ,
where ϕ = ϕ ( x ; x 0 ) , χ = χ ( x ; x 0 ) and ψ = ψ ( x ; x 0 ) . By substituting Equation (51) into Equation (43) and applying some manipulation, such as 𝛻 ( 𝛻 ϕ ) = 𝛻 2 ϕ , 𝛻 × ϕ = 0 , 𝛻 ψ = 𝛻 χ = 0 , and 𝛻 × ( 𝛻 𝛻 ψ ) = 𝛻 × ( 𝛻 𝛻 χ ) = 0 , Equation (43) can be rewritten as
𝛻 𝛻 · [ ( λ + 2 μ ) ( 𝛻 2 P ϕ ) + ρ ω 2 ( P ϕ ) + P g ] + a 𝛻 × [ μ 𝛻 4 ( P ψ ) + ρ ω 2 𝛻 2 ( P ψ ) ] + 𝛻 × 𝛻 × [ μ 𝛻 2 ( P χ ) + ρ ω 2 ( P χ ) + P g ] = 0 .
In the above equation, the components in the gradient, the curl, and the curl of curl, can be independently zero
( λ + 2 μ ) 𝛻 2 ϕ + ρ ω 2 ϕ + g = 0 ,
μ 𝛻 4 ψ + ρ ω 2 𝛻 2 ψ = 0 ,
μ 𝛻 2 χ + ρ ω 2 χ + g = 0 .
Introducing the longitudinal wave (P) speed ( c P = λ + 2 μ / ρ ) and transverse wave (S) speed ( c S = μ / ρ ) into the above equations leads to
𝛻 2 ϕ + k p 2 ϕ + k p 2 ρ ω 2 g = 0 ,
𝛻 2 ψ + k s 2 ψ = 0 ,
𝛻 2 χ + k s 2 χ + k s 2 ρ ω 2 g = 0 ,
where k p = ω c P and k s = ω c S , correspond to the angular wavenumbers of the P and the S waves, respectively. Taking into account the Laplacian operator in the ( ξ , ϑ , η ) cylindrical coordinates, the three scalar spatial potentials have the following forms:
ϕ ( ξ , ϑ , η ) = ϕ r ( ξ ) ϕ θ ( ϑ ) ϕ z ( η ) ,
χ ( ξ , ϑ , η ) = χ r ( ξ ) χ θ ( ϑ ) χ z ( η ) ,
ψ ( ξ , ϑ , η ) = ψ r ( ξ ) ψ θ ( ϑ ) ψ z ( η ) .
First, let us solve PDE of Equation (56). Using a similar method to that applied to the case of Green’s function, the axial and the angular parts can be obtained as
ϕ z ( η ) = A z e i k η η ,
ϕ θ ( ϑ ) = A m θ cos m ϑ + B m θ sin m ϑ     ( m = 0 , ± 1 , ± 2 ,   ) ,
respectively. Substituting the Laplacian operator and Equations (41), (59), (62), and (63), and introducing some straightforward algebra to Equation (56) leads to
2 ϕ r ξ 2 + 1 ξ ϕ r ξ + ( α 2 m 2 ξ 2 ) ϕ r = k p 2 ρ ω 2 ( 1 ϕ ϑ ϕ η ) G 1 J 0 ( r 01 a r 0 ξ ) ,
where α 2 = k p 2 k η 2 . Equation (64) is a second-order linear nonhomogeneous PDE, the solution of which is a linear combination of the homogeneous ( ϕ r h ) and particular ( ϕ r p ) solutions. The homogeneous solution of Equation (64) is
ϕ r h = { A m r [ J m ( α ξ ) I m ( α ξ ) ] + B m r [ Y m ( α ξ ) K m ( α ξ ) ] }                 ( α 2 0 )           ( α 2 < 0 ) ,
where J m and I m are the first kind and its modified Bessel functions, respectively, and Y m and K m are the second kind and its modified Bessel functions, respectively. The Bessel functions of the second kind are excluded due to the singularity at the origin included in the domain of the cylinder. For the case of α 2 0 ,
ϕ r h = A m r J m ( α ξ ) .
According to Korenev’s solutions [28], the particular solution of the inhomogeneous equation has the form
ϕ r p = k p 2 ρ ω 2 ( 1 ϕ ϑ ϕ η ) G 1 1 1 [ r 01 α ( a r 0 ) ] 2 J m [ r 01 α ( a r 0 ) ξ ] ,                 r 01 α ( a r 0 ) 1 .
The linear combination of Equations (65) and (66) results in
ϕ ( ξ , ϑ , η ) = ϕ r ( ξ ) ϕ θ ( ϑ ) ϕ z ( η ) = [ A m r J m ( α ξ ) ] [ A m θ cos ( m ϑ ) ] ( A z e i k η η ) k p 2 ρ ω 2 G 1 1 1 [ r 01 α ( a r 0 ) ] 2 J m [ r 01 α ( a r 0 ) ξ ] .
For further applications, the above equation is rewritten as
ϕ ( ξ , ϑ , η ) = A m   J m ( α ξ ) cos ( m ϑ ) e i k η η k p 2 ρ ω 2 G 1 1 1 [ r 01 α ( a r 0 ) ] 2 J m [ r 01 α ( a r 0 ) ξ ] ,
where the coupling constant A m = A m r A m θ A z .
Following a similar procedure, the functions of χ ( ξ , ϑ , η ) and ψ ( ξ , ϑ , η ) can be obtained as
χ ( ξ , ϑ , η ) = B m J m ( β ξ ) sin ( m ϑ ) e i k η η k p 2 ρ ω 2 G 1 1 1 [ r 01 β ( a r 0 ) ] 2 J m [ r 01 β ( a r 0 ) ξ ] ,
ψ ( ξ , ϑ , η ) = C m   J m ( β ξ ) cos ( m ϑ ) e i k z η .
As shown in Figure 1, we introduce force vector P acting in the radial and the axial directions to solve for Φ. Equation (45) can be rewritten as
Φ = 𝛻 · P ϕ = ( P ϕ ) ξ i + ( P ϕ ) ξ j + ( P ϕ ) η ,
where ξ i = x i x 0 i , ξ j = x j x 0 j , and η = x z x 0 z .
For P = P j ,
Φ j = P j Ξ ϕ ξ ,
where
Ξ = ξ i ξ δ i j + ξ j ξ = a   cos φ x 0 i ξ i 2 + ξ j 2 δ i j + a   sin φ x 0 j ξ i 2 + ξ j 2
φ is the angle between an observation point and the x i axis.
For P = P z ,
Φ z = P z ϕ η .
By substituting Equation (67) into Equations (71) and (73), we obtain
Φ j = P j Ξ { A m J m ( α ξ ) ξ cos ( m ϑ ) e i k η η k p 2 ρ ω 2 G 1 1 1 [ r 01 α ( a r 0 ) ] 2 J m [ r 01 α ( a r 0 ) ξ ] ξ } ,
Φ z = P z { ( i k η ) A m   J m ( α ξ ) cos ( m ϑ ) e i k z η k p 2 ρ ω 2 G 1 η 1 1 [ r 01 α ( a r 0 ) ] 2 J m [ r 01 α ( a r 0 ) ξ ] } .
In Equation (75),
G 1 η = J 0 ( r 01 a r 0 ξ 0 ) π ( a r 0 ) r 01 [ e r 01 z 0 a r 0 + e r 01 ( l z 0 ) a r 0 2 ] [ J 1 ( r 01 ) ] 2   × { ( r 01 a r 0 ) e r 01 η / ( a r 0 )                     z 0 l < η < 0 ( r 01 a r 0 ) e r 01 η / ( a r 0 )           0 < η < l z 0 .
From Equations (46) and (68), we obtain the SH potential:
Χ e ^ z = 𝛻 × P χ = [ ( P k χ ξ j P j χ η ) i ^ + ( P i χ η P k χ ξ i ) j ^ + ( P j χ ξ i P i χ ξ j ) e ^ z ] = ( P j χ ξ i P i χ ξ j ) e ^ z
For P = P j ,
Χ j = P j ( χ ξ i δ j i χ ξ j ) = P j ( ξ i ξ ξ j ξ δ j i ) χ ξ = P j Σ χ ξ ,
Χ j = P j Σ { B m J m ( β ξ ) ξ sin ( m ϑ ) e i k z η k p 2 ρ ω 2 G 1 1 1 [ r 01 β ( a r 0 ) ] 2 J m [ r 01 β ( a r 0 ) ξ ] ξ } .
In the above equation,
Σ = ξ i ξ ξ j ξ δ i j = a   cos φ x 0 i ξ i 2 + ξ j 2 a   sin φ x 0 j ξ δ i j .
For P = P z ,
Χ z = P z ( χ ξ i δ i z χ ξ j δ j z ) = 0 .
Similarly, for the SV potential, we find
Ψ j = P j Σ C m J m ( β ξ ) ξ cos ( m ϑ ) e i k z η ,
Ψ z = 0 .
All three potentials for the P, SH, and SV waves have now been completely determined, and we can thus derive the displacement components of u r , u θ and u z in terms of the three potentials. The relations between the displacements and potentials are given in Appendix A. First, let us derive the components generated by Pj. It can be easily shown that by substituting Equations (74), (77), and (80) into Equations (A5)–(A7), and taking the location of the point source as the origin ( ξ = 0 ,   ϑ = 0 , η = 0 ) , the displacement d component is reduced to
u d j = P j ( A m j F d j 1 + B m j F d j 2 + C m j F d j 3 + F d j 4 ) e i ω t ,
where subscript d represents the radial (r), the tangential (θ), or the axial (z) components.
For the radial component u r j ,
F r j 1 = [ Ξ 2 J m ( α ξ ) ξ 2 ] cos ( m ϑ ) e i k z η ,
F r j 2 = m Σ ξ J m ( β ξ ) ξ cos ( m ϑ ) e i k z η ,
F r j 3 = i k z a Σ [ 2 J m ( β ξ ) ξ 2 cos ( m ϑ ) ] e i k z η ,
F r j 4 = Ξ ( k p 2 ρ ω 2 ) G 1 1 1 [ r 01 α ( a r 0 ) ] 2 [ 2 J m [ r 01 α ( a r 0 ) ξ ] ξ 2 ] .
For the tangential component u θ j ,
F θ j 1 = m Ξ ξ J m ( α ξ ) ξ sin ( m ϑ ) e i k z η ,
F θ j 2 = [ Σ 2 J m ( β ξ ) ξ 2 ] sin ( m ϑ ) e i k z η ,
F θ j 3 = i k z a m Σ ξ J m ( β ξ ) ξ sin ( m ϑ ) e i k z η ,
F θ j 4 = Σ ( k s 2 ρ ω 2 ) { G 1 1 1 [ r 01 β ( a r 0 ) ] 2 [ 2 J m [ r 01 β ( a r 0 ) ξ ] ξ 2 ] } .
For the axial component u z j ,
F z j 1 = i k z   Ξ J m ( α ξ ) ξ cos ( m ϑ ) e i k z η ,
F z j 2 = 0 ,
F z j 3 = a Σ { 3 J m ( β ξ ) ξ 3 + 1 ξ 2 J m ( β ξ ) ξ 2 m 2 ξ 2 J m ( β ξ ) ξ } cos ( m ϑ ) e i k z η ,
F z j 4 = Ξ ( k p 2 ρ ω 2 ) { G 1 η 1 1 [ r 01 α ( a r 0 ) ] 2 [ J m [ r 01 α ( a r 0 ) ξ ] ξ ] } .
Following similar procedures, the three displacements generated by Pz can be derived by substituting Equations (75), (79), and (81) into Equations (A6)–(A8):
u d z = P z ( A m z F d z 1 + B m z F d z 2 + C m z F d z 3 + F d z 4 ) e i ω t
For the radial component u r z ,
F r z 1 = i k z J m ( α ξ ) ξ cos ( m θ ) e i k z η ,
F r z 2 = F r z 3 = 0 ,
F r z 4 = ( k p 2 ρ ω 2 ) { G 1 η 1 1 [ r 01 α ( a r 0 ) ] 2 [ J m [ r 01 α ( a r 0 ) ξ ] ξ ] } .
For the tangential component u θ z ,
F θ z 1 = i k z m ξ J m ( α ξ )   sin ( m ϑ ) e i k z η ,
F θ z 2 = F θ z 3 = F θ z 4 = 0 .
For the axial component u z z ,
F z z 1 = k z 2 J m ( α ξ ) cos ( m ϑ ) e i k z η ,
F z z 2 = F z z 3 = 0 ,
F z z 4 = ( k p 2 ρ ω 2 ) { 2 G 1 η 2 1 1 [ r 01 α ( a r 0 ) ] 2 [ J m [ r 01 α ( a r 0 ) ξ ] ξ ] } .
The only remaining task to complete the displacement fields is to determine the coupling constants Am, Bm, and Cm. These constants can be determined directly by applying a fundamental set of linear elastic boundary problems. The outer surface of the cylinder studied in the present paper is stress-free. Thus, the following stress components are zero under these circumstances, i.e., ξ = a r 0 :
σ r r = σ r θ = σ r z = 0 .
Substituting Equations (A17)–(A21) into Equation (104) yields the following algebraic equations
[ a 11 f a 12 f a 13 f a 21 f a 22 f a 23 f a 31 f a 32 f a 33 f ] [ A m f B m f C m f ] = [ b 1 f b 2 f b 3 f ] ,
where f = j for Pj and f = z for Pz. The elements in Equation (105) are given in Appendix A. For Pj, the elements of b1j, b2j and b3j are non-zero; therefore, the coupling constants can be determined by solving Equation (105).
The CF in Equations (1) and (43) is the force impulse defined as
P ( t ) = P 0 t   exp ( b t ) ,
where P 0 and b are parameters that determine the amplitude and the duration of the wave, respectively [29]. The arrival time τ of the impulse at the position ( ξ i , ξ j , η ) must be considered, since the force impulse is acting in the f direction at position ξ i = ξ j = 0 and η = 0 at t = 0. Replacing t with (tτ) yields
P f ( t τ ) = P 0 , f ( t τ ) e b ( t τ ) .
The displacement components generated by the Φ f potential, and the Χ f and Ψ f potentials, correspond to the compressional (P), and the shear (SH and SV) waves, respectively; therefore, the arrival times of the two waves are given as
τ P = ξ i 2 + ξ j 2 + η 2 c P ,   τ S = ξ i 2 + ξ j 2 + η 2 c S ,
where c P and c S are the velocities of the P and S waves, respectively. The displacements generated by Pf can be rewritten as
u r f = [ P f ( t τ P )   ( A m f F r f 1 + F r f 4 ) + P f ( t τ s ) ( B m f F r f 2 + C m f F r f 3 ) ] e i ω t ,
u θ f = [ P f ( t τ P )   ( A m f F θ j 1 ) + P f ( t τ s ) ( B m f F θ f 2 + C m f F θ f 3 + F θ f 4 ) ] e i ω t ,
u z f = [ P f ( t τ P )   ( A m f F z f 1 + F z f 4 ) + P f ( t τ s ) ( B m f F z f 2 ) ] e i ω t .
Furthermore, the P, SH, and SV waves can be obtained from Equations (109)–(111) as
u f P = P f ( t τ P ) ( A m f F d f 1 + F r f 4 ) e i ω t   ( d = r ,   z   and   θ ) ,
u f S V = P f ( t τ S ) ( C m f F d f 3 ) e i ω t   ( d = r ,   z   and   θ ) ,
u f S H = P f ( t τ S ) ( B m f F d f 2 + F θ f 4 ) e i ω t   ( d = r ,   z   and   θ ) .

4. Simulation

In this study, we consider only the azimuthal dependence of the wave propagation arising entirely from the direction of the force vector, m = 0 . A stainless-steel (SS) cylinder (a = 0.50 m, l = 2.0 m, ρ = 7.80 × 103 kg/m3, c P = 5.98 km/s, and c S = 3.30 km/s) was used as the test specimen. Stiffness parameters ( c 11 = c 22 = c 33 = 2.086 × 10 11 ,   c 12 = c 13 = c 23 = 1.465 × 10 11 ,   c 44 = c 55 = c 66 = 1.269 × 10 11   N   m 2 ) can be found for austenitic stainless steel [30]. Previously, the natural frequencies of SS active on AE in the range of 0–300 kHz were determined experimentally by a tensile test (Appendix B). The frequency of ν = 155.4 kHz (the angular frequency ω = 2 π ν ) was found to be predominant.
First, we determined k z by applying the boundary conditions to Equations (A20) and (A22) which were associated with the displacements due to P z . From Equation (105), two solutions of A 0 were obtained, as follows:
A 0 z 1 = b 1 z a 11 z       or       A 0 z 2 = b 3 z a 31 z  
since a 13 z = a 33 z = 0 .
When m = 0,
a 11 z = i k z [ c 13 k z 2 J 0 ( α ξ ) c 12 ξ J 0 ( α ξ ) ξ c 11 2 J 0 ( α ξ ) ξ 2 ] e i k z η = i k z { c 13 k z 2 J 0 ( α ξ ) + c 12 ξ α J 1 ( α ξ ) c 11 2 ( α 2 ) 2 [ J 0 ( α ξ ) + J 2 ( α ξ ) ] } e i k z η ,
b 1 z = ( k p 2 ρ ω 2 ) ( 1 1 [ r 01 α ( a r 0 ) ] 2 ) { c 11 G 1 η 2 J 0 [ r 01 α ( a r 0 ) ξ ] ξ 2 + c 12 1 ξ G 1 η J 0 [ r 01 α ( a r 0 ) ξ ] r + c 13 3 G 1 η 3 J 0 [ r 01 α ( a r 0 ) ξ ] } = ( k p 2 ρ ω 2 ) ( 1 1 [ r 01 α ( a r 0 ) ] 2 ) c 11 G 1 η 2 [ r 01 2 α ( a r 0 ) ] 2 { J 0 [ r 01 α ( a r 0 ) ξ ] + J 2 [ r 01 α ( a r 0 ) ξ ] } c 12 1 ξ G 1 η r 01 α ( a r 0 ) J 1 [ r 01 α ( a r 0 ) ξ ] + c 13 3 G 1 η 3 J 0 [ r 01 α ( a r 0 ) ξ ] ,
a 31 z = c 44 k z 2 α J 1 ( α ξ ) e i k z η ,
b 3 z = c 44 ( k p 2 ρ ω 2 ) { 2 G 1 η 2 1 1 [ r 01 α ( a r 0 ) ] 2 2 J 0 [ r 01 α ( a r 0 ) ξ ] ξ } = 2 c 44 ( k p 2 ρ ω 2 ) 2 G 1 η 2 1 1 [ r 01 α ( a r 0 ) ] 2 r 01 α ( a r 0 ) J 1 [ r 01 α ( a r 0 ) ξ ] .
In Equations (117) and (119), the partial derivatives of G 1 with respect to η can easily be obtained from Equation (42). On the circumference, ξ = a r 0 and η = l z 0 , we solved the roots of the function
f ( k z ) = b 1 z a 11 z b 3 z a 31 z = 0 ,
as a function of the shortest distance from the point source to the end plate of the cylinder (η) at a given ξ . It should be noted that f ( k z ) is independent of Ξ and Σ. Figure 2 shows the η dependence of three roots (n = 1–3) of Equation (120) for r 0 = 0   ( x 0 i = x 0 j = 0   m ) and 0.25 ( x 0 i = 0 , x 0 j = 0.25   m ) . In the simulation, we selected the first root (n = 1) at the given r 0 and z 0 values.
Next, we determined P0 and b in Equation (107). The envelope of a given wave was evaluated by fitting its normalized form to determine b. As shown in Figure 3, when b = 6.0 × 104 s−1, the duration (Δτ) of the wave was approximately 1 ms. By increasing the b value, the duration of the envelope decreased exponentially at b = 4.0 × 105 s−1 with Δτ ≅ 25 μs. In the simulation, the force exerted by the point source was fixed to 1 N with P0 = 1.0 × 1010 N s−1 and b = 1.0 × 105 s−1.
In the simulation two positions of the point source of 1 N were considered, with coordinates of (0 m, 0 m, 1 m) and (0 m, 0.25 m, 1 m). Figure 4a shows the displacements, and their wave properties, at the (0.5 m, 45°, 1 m) position, generated by the P f point source located at the center of the cylinder (0 m, 0 m, 1 m). The Pj excitation produces an axial displacement that is stronger than the radial and tangential displacements. These displacements result in the P wave being the main wave, with a minor SH wave and very weak SV wave (Figure 4b). As shown in Figure 4c, the displacement features generated by the Pz excitation differ significantly from those generated by the Pj excitation. For the Pj excitation, the maximum values of u r j and u z j at the (0.5 m, 45°, 1 m) position were 0.29 and 0.71 nm, respectively, whereas for the Pz excitation, they were 28.1 and 5.0 nm, respectively. The amplitudes of the displacements due to the Pz excitation were much stronger than those due to the Pj excitation. It should be noted that the Pz excitation produces only the P wave (Figure 4d). The angular dependence of the displacement was also calculated (Figure 5). When the point source is located at the center of the circular plane, the angular dependences of u r j , u z j and u t j arise from only Ξ and Σ as defined in Equations (72) and (78), respectively. For Pz excitation, the displacements of u r z and u z z are free from these factors. When the distances from the point source to the circumference are not equivalent, the angular dependences of the radial and axial displacements are highly significant. These effects are due not only to Ξ and Σ but also the superposition of the Bessel functions involved in the displacement equations. At a certain angle, some a i j values in Equation (105) become too small to cause the sudden increase in displacement.

5. Conclusions

In this paper, we present a mathematical model for AE generated by a point source in a TIC. The point source as the CF vector is a model for an initiation state of crack formation in a solid medium. In the cylindrical system, the displacement field vector is very complex: compressional (P), and vertical and horizontal shear (SV and SH) wave potentials are coupled. Introducing CFIPs into the NL equation leads to the solutions of the radial, tangential, and axial displacements in the cylindrical system. The solutions obtained in this study can be used for analyzing experimental data for the NDT of cylinders. In the near future, similar mathematical models will be developed for shells and multilayered cylinders.

Author Contributions

Conceptualization, methodology and investigation, K.B.K. and J.-G.K.; software, validation, and formal analysis, K.B.K.; resources and data curation, S.G.L.; writing—original draft preparation, K.B.K. and J.-G.K.; writing—review and editing, visualization, and supervision, J.-G.K.; project administration and funding acquisition, B.K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) and the Ministry of Trade, Industry & Energy (MOTIE) of the Republic of Korea (NO. 20202910100070).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Some gradient, divergence and curl operators on the potentials are as follows:
𝛻 Φ = Φ r r ^ + 1 r Φ θ θ ^ + Φ z z ^ ,
𝛻 × ( Χ e ^ z ) = [ 1 r ( Χ e ^ z ) z θ ( Χ e ^ z ) θ z ] r ^ + [ ( Χ e ^ z ) r z ( Χ e ^ z ) z r ] θ ^ + 1 r { [ r ( Χ e ^ z ) θ ] r ( Χ e ^ z ) r θ } z ^ = 1 r Χ θ r ^ Χ r θ ^ ,
𝛻 × 𝛻 × ( Ψ e ^ z ) = 𝛻 × ( 1 r Ψ θ r ^ Ψ r θ ^ ) = [ ( Ψ r ) z ] r ^ + [ ( 1 r Ψ θ ) z ] θ ^ + 1 r [ [ r ( Ψ r ) ] r ( 1 r Ψ θ ) θ ] z ^ = ( 2 Ψ r z ) r ^ + ( 1 r 2 Ψ θ z ) θ ^ 1 r ( Ψ r + r 2 Ψ r 2 + 1 r 2 Ψ θ 2 ) z ^
Substituting the above equations into Equation (44) gives
u = ( Φ r + 1 r Χ θ + a 2 Ψ r z ) r ^ + ( 1 r Φ θ + Χ r + a 1 r 2 Ψ θ z ) θ ^ + [ Φ z a ( 2 Ψ r 2 + 1 r Ψ r + 1 r 2 2 Ψ θ 2 ) ] z ^ .
From Equation (A5), we obtain the displacement components as
u r = Φ r + 1 r Χ θ + a 2 Ψ r z ,
u θ = 1 r Φ θ + Χ r + a 1 r 2 Ψ θ z ,
u z = Φ z a ( 2 Ψ r 2 + 1 r Ψ r + 1 r 2 2 Ψ θ 2 ) ,
in which ur, uθ and uz denote the radial, tangential, and axial displacements, respectively. The radial component can be obtained by substituting Equations (71), (77), and (79) into Equation (A6) in the ( ξ , ϑ , η ) coordinates.
u r j = Φ j ξ + 1 ξ Χ j ϑ + a 2 Ψ j ξ η = P j Ξ { A m cos ( m ϑ ) e i k z η 2 J m ( α ξ ) ξ 2 ( k p 2 ρ ω 2 ) G 1 1 1 [ r 01 α ( a r 0 ) ] 2 2 J m [ r 01 α ( a r 0 ) ξ ] ξ 2 } + m Σ ξ [ B m J m ( β ξ ) ξ cos ( m ϑ ) e i k z η ] + a ( i k z ) Σ C m [ 2 J m ( β ξ ) ξ 2 cos ( m ϑ ) ] e i k z η e i ω t ,
where Ξ ξ 0 and Σ ξ 0 . By rearranging the above equation in terms of the Am, Bm and Cm coupling constants, we obtain
u r j = P j [ m = 0 ( A m F r j 1 + B m F r j 2 + C m F r j 3 ) + F r j 4 ] e i ω t ,
where
F r j 1 = [ Ξ 2 J m ( α ξ ) ξ 2 ] cos ( m ϑ ) e i k z η ,
F r j 2 = m ξ J m ( β ξ ) ξ Σ cos ( m ϑ ) e i k z η ,
F r j 3 = i k z a Σ [ 2 J m ( β ξ ) ξ 2 cos ( m ϑ ) ] e i k z η ,
F r j 4 = ( k p 2 ρ ω 2 ) Ξ [ G 1 1 1 [ r 01 α ( a r 0 ) ] 2 2 J m [ r 01 α ( a r 0 ) ξ ] r 2 ] .
For TIC, the stress-strain [18] and the strain-displacement [31] relations are given as
[ σ r r σ θ θ σ z z σ θ z σ r z σ r θ ] = [ 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 ( c 11 c 12 ) / 2 ] [ ϵ r r ϵ θ θ ϵ z z 2 ϵ θ z 2 ϵ r z 2 ϵ r θ ] ,
and
ϵ r r = u r r ,   ϵ θ θ = 1 r u θ θ + u r r ,   ϵ θ z = 1 2 ( u θ z + 1 r u z θ ) ϵ z z = u z z ,   ϵ r θ = 1 2 ( 1 r u r θ + u θ r u θ r ) ,   ϵ r z = 1 2 ( u z r + u r z ) ,
respectively. From these relations, one can obtain three stresses of σ r r , σ r ϕ and σ r z as [18]:
σ r r = c 11 ( u r r ) + c 12 ( u r r + 1 r u θ θ ) + c 13 ( u z z ) ,
σ r θ = ( c 11 c 12 ) 2 ( u θ r u θ r + 1 r u r θ ) ,
σ r z = c 44 ( u z r + u r z )
From u r j , u θ j and u z j
σ r r j = P j ( a 11 j A m j + a 12 j B m j + a 12 j C m j + b 1 j ) e i ω t ,
where
a 11 j = c 11 F r j 1 ξ + c 12 F r j 1 ξ + c 12 1 ξ F θ j 1 ϑ + c 13 ( F z j 1 η ) = Ξ { ( c 12 m 2 ξ 2 + c 13 k z 2 ) J m ( α ξ ) ξ + c 12 1 ξ 2 J m ( α ξ ) ξ 2 + c 11 3 J m ( α ξ ) ξ 3 } cos ( m ϑ ) e i k z η ,
a 12 j = c 11 F r j 2 ξ + c 12 F r j 2 ξ + c 12 1 ξ F θ j 2 ϑ = ( c 11 + c 12 ) m Σ ξ [ 1 ξ J m ( β ξ ) ξ 2 J m ( β ξ ) ξ 2 ] cos ( m ϑ ) e i k z η ,
a 13 j = c 11 F r j 3 r + c 12 F r j 3 r + c 12 1 r F θ j 3 θ + c 13 F z j 3 z = i k z a Σ { ( c 12 + c 13 ) 1 ξ [ m 2 ξ J m ( β ξ ) ξ + 2 J m ( β ξ ) ξ 2 ] + ( c 11 + c 13 ) 3 J m ( β ξ ) ξ 3 } e i k z η ,
b 1 j = c 11 F r j 4 r + c 12 F r j 4 r + c 13 F z j 4 η = ( k p 2 ρ ω 2 ) Ξ 1 1 [ r 01 α ( a r 0 ) ] 2 { c 11 G 1 3 J 0 [ r 01 α ( a r 0 ) ξ ] ξ 3 + c 12 1 ξ G 1 2 J 0 [ r 01 α ( a r 0 ) ξ ] ξ 2 + c 13 2 G n η 2 J 0 [ r 01 α ( a r 0 ) ξ ] ξ } ,
σ r θ j = P j ( a 21 j A m j + a 22 j B m j + a 23 j C m j + b 2 j ) e i ω t ,
where
a 21 j = ( c 11 c 12 ) 2 ( F θ j 1 ξ F θ j 1 ξ + 1 ξ F r j 1 ϑ ) = ( c 11 c 12 ) 2 ( 2 m Ξ ) [ 1 ξ 2 J m ( α ξ ) ξ 1 ξ 2 J m ( α ξ ) ξ 2 ] sin ( m ϑ ) e i k z η ,
a 22 j = ( c 11 c 12 ) 2 ( F θ j 2 r F θ j 2 r + 1 r F r j 2 θ ) = ( c 11 c 12 ) 2 Σ [ m 2 ξ J m ( β ξ ) ξ + ( 2 1 ξ ) 2 J m ( β ξ ) ξ 2 + 3 J m ( β ξ ) ξ 3 ] sin ( m ϑ ) e i k z η ,
a 23 j = ( c 11 c 12 ) 2 ( F θ j 3 r F θ j 3 r + 1 r F r j 3 θ ) = ( c 11 c 12 ) 2 ( i k z a m Σ ξ ) [ 2 ξ 2 J m ( β ξ ) ξ 2 2 J m ( β ξ ) ξ 2 ] sin ( m ϑ ) e i k z η ,
b 2 j = ( c 11 c 12 ) 2 ( F ϑ j 4 ξ F ϑ j 4 ξ + 1 ξ F r j 4 ϑ ) = ( c 11 c 12 ) 2 ( k s 2 ρ ω 2 ) Σ G 1 1 [ r 01 β ( a r 0 ) ] 2 { 1 ξ 2 J m [ r 01 β ( a r 0 ) ξ ] ξ 2 + 3 J m [ r 01 β ( a r 0 ) ξ ] ξ 3 } ,
and
σ r z j = ( a 31 j A m j + a 32 j B m j + a 33 j C m j + b 3 j ) e i ω t ,
where
a 31 j = c 44 ( F z j 1 ξ + F r j 1 η ) = 2 c 44 ( i k z Ξ ) 2 J m ( α ξ ) ξ 2 cos ( m ϑ ) e i k z z ,
a 32 j = c 44 ( F r j 2 η ) = c 44 ( i k z m Σ ξ ) J m ( β ξ ) ξ cos ( m ϑ ) e i k z η ,
a 33 j = c 44 ( F z j 3 ξ + F r j 3 η ) = c 44 a Σ [ ( m 2 ξ 2 ) ( 1 2 ξ i k z ) J m ( β ξ ) ξ + 1 ξ ( 1 + m 2 ξ + i k z ) 2 J m ( β ξ ) ξ 2 ( 1 ξ i k z ) 3 J m ( β ξ ) ξ 3 4 J m ( β ξ ) ξ 4 ] cos ( m ϑ ) e i k z η ,
b 3 j = c 44 ( F z j 4 ξ + F r j 4 η ) = c 44 ( k p 2 ρ ω 2 ) Ξ { 1 1 ( r 0 n α a ) 2 [ G 1 η 2 J m [ r 01 α ( a r 0 ) ξ ] ξ 2 + 2 G n η 2 J m [ r 01 α ( a r 0 ) ξ ] ξ ] } .
From u r z , u θ z and u z z
σ r r z = P z ( a 11 z A m z + a 12 z B m z + a 12 j C m z + b 1 z ) e i ω t ,
where
a 11 z = i k z [ ( c 12 m 2 ξ 2 + c 13 k z 2 ) J m ( α ξ ) c 12 J m ( α ξ ) ξ c 11 2 J m ( α ξ ) ξ 2 ] cos ( m ϑ ) e i k z η ,
a 12 z = 0 ,
a 13 z = 0 ,
b 1 z = ( k p 2 ρ ω 2 ) 1 1 [ r 01 α ( a r 0 ) ] 2 { c 11 G 1 η 2 J m [ r 01 α ( a r 0 ) ξ ] ξ 2 + c 12 1 ξ G 1 η J m [ r 01 α ( a r 0 ) ξ ] r + c 13 3 G 1 η 3 J 0 [ r 01 α ( a r 0 ) ξ ] } ,
σ r θ z = P z ( a 21 z A m + a 22 z B m + a 23 z C m + b 2 z ) e i ω t ,
where
a 21 z = ( c 11 c 12 ) 2 [ i k z m ξ 2 J m ( α ξ ) + i k z m ξ J m ( α ξ ) ξ i k z m ξ 2 J m ( α ξ ) + i k z m ξ J m ( α ξ ) ξ ] cos ( m ϑ ) e i k z η ,
a 22 z = 0 ,
a 23 z = 0 ,
b 2 z = 0 ,
and
σ r z z = ( a 31 z A m z + a 32 z B m z + a 33 z C m j + b 3 z ) e i ω t ,
where
a 31 z = c 44 k z 2 J m ( α ξ ) ξ cos ( m θ ) e i k z z ,
a 32 z = 0 ,
a 33 z = 0 ,
b 3 z = c 44 ( k p 2 ρ ω 2 ) { 2 G 1 η 2 1 1 [ r 01 α ( a r 0 ) ] 2 2 J m [ r 01 α ( a r 0 ) ξ ] ξ } .

Appendix B

Figure A1. Distribution of resolved frequency components and the 5 groups selected using a 10 kHz wide window. The experimental data are obtained from a tensile test using AISI 316 austenitic stainless steel (KS D 3698 standard).
Figure A1. Distribution of resolved frequency components and the 5 groups selected using a 10 kHz wide window. The experimental data are obtained from a tensile test using AISI 316 austenitic stainless steel (KS D 3698 standard).
Applsci 12 05272 g0a1
Table A1. Resolved frequency components (ν’s) and their relative weights of 316 stainless steel.
Table A1. Resolved frequency components (ν’s) and their relative weights of 316 stainless steel.
GroupElements aν/kHz b
120(0.07)103.2(2.62)
258(0.21)127.6(5.59)
3114(0.41)155.4(5.37)
436(0.13)173.6(3.28)
551(0.18)190.6(4.72)
a Value in parenthesis represents the fraction. b Value in parenthesis represents the standard deviation.

References

  1. Seruby, C.B. An Introduction to acoustic emission. J. Phys. E Sci. Instrum. 1987, 20, 946–953. [Google Scholar]
  2. Godin, N.; Reynaud, P.; Fantozzi, G. Challenges and Limitations in the Identification of Acoustic Emission Signature of Damage Mechanisms in Composites Materials. Appl. Sci. 2018, 8, 1267. [Google Scholar] [CrossRef] [Green Version]
  3. Domaneschi, M.; Niccolini, G.; Lacidogna, G.; Cimellaro, G.P. Nondestructive Monitoring Techniques for Crack Detection and Localization in RC Elements. Appl. Sci. 2020, 10, 3248. [Google Scholar] [CrossRef]
  4. Oh, T.-M.; Kim, M.-K.; Lee, J.-W.; Kim, H.; Kim, M.-J. Experimental Investigation on Effective Distances of Acoustic Emission in Concrete Structure. Appl. Sci. 2020, 10, 6051. [Google Scholar] [CrossRef]
  5. Goldaran, R.; Turer, A.; Kouhdaragj, M.; Ozlutas, K. Identification of Corrosion in a Prestressed Concrete Pipe Utilizing Acoustic Emission Technique. Constr. Build. Mater. 2020, 242, 118053. [Google Scholar] [CrossRef]
  6. Dong, L.; Tong, X.; Ma, J. Quantitative Investigation of Tomographic Effects in Abnormal Regions of Complex Structures. Engineering 2021, 7, 1011–1022. [Google Scholar] [CrossRef]
  7. Chai, S.; Zhang, Z.; Duan, Q. A New Qualitative Acoustic Emission Parameter Based on Shannon’s Entropy for Damage Monitoring. Mech. Syst. Signal Process. 2018, 100, 617–629. [Google Scholar] [CrossRef]
  8. Krampikowska, A.; Pała, R.; Dzioba, I.; Świt, G. The Use of the Acoustic Emission Method to Identify Crack Growth in 40CrMO Steel. Materials 2019, 12, 2140. [Google Scholar] [CrossRef] [Green Version]
  9. Kovač, J.; Legat, A.; Zajec, B.; Kosec, T.; Govekar, E. Detection and Characterization of Stainless Steel SCC by the Analysis of Crack Related Acoustic Emission. Ultrasonics 2015, 62, 312–322. [Google Scholar] [CrossRef] [Green Version]
  10. Kek, T.; Kusič, D.; Svečko, R.; Hančič, A.; Grum, J. Acoustic Emission Signal Analysis for the Integrity Evaluation. J. Mech. Eng. 2018, 64, 665–671. [Google Scholar]
  11. Mukherjee, A.; Banerjee, A. Analysis of Acoustic Emission Signal for Crack Detection and Distance Measurement on Steel Structure. Acoust. Aust. 2021, 49, 133–149. [Google Scholar] [CrossRef]
  12. Aki, K.; Richards, P.G. Quantitative Seismology, 2nd ed.; University Science Books: Mill Valley, CA, USA, 2009; Chapter 4. [Google Scholar]
  13. Pujol, J.; Herrmann, B. A Student’s Guide to Point Sources in Homogeneous Media. Seismolog. Res. Lett. 1990, 61, 209–221. [Google Scholar] [CrossRef]
  14. Scruby, C.B.; Wadley, H.N.G.; Hill, J.J. Dynamic Elastic Displacements at the Surface of an Elastic Half-space due to Defect Sources. J. Phys. D Appl. Phys. 1983, 16, 1069–1083. [Google Scholar] [CrossRef]
  15. Ohtsu, M.; Ono, K. The Generalized Theory and Source Representation of Acoustic Emission. J. Acoust. Emiss. 1986, 5, 124–133. [Google Scholar]
  16. Weaver, R.L.; Pao, Y.-H. Axisymmetric Elastic Waves Excited by a Point Source in a Plate. J. Appl. Mech. 1982, 49, 821–836. [Google Scholar] [CrossRef]
  17. Mal, A.K.; Singh, S.J. Deformation of Elastic Solids; Prentice Hall: New York, NY, USA, 1991; pp. 312–317. [Google Scholar]
  18. Honarvar, F.; Sinclair, A.N. Acoustic Wave Scattering from Transversely Isotropic Cylinders. J. Acoust. Soc. Am. 1996, 100, 57–63. [Google Scholar] [CrossRef]
  19. Honarvar, F.; Enjilela, E.; Sinclair, A.N.; Mirnezami, S.A. Wave Propagation in Transversely Isotropic Cylinders. Int. J. Solids Struct. 2007, 44, 5236–5246. [Google Scholar] [CrossRef] [Green Version]
  20. Ponnusamy, P.; Selvamani, R. Wave Propagation in a Homogeneous Isotropic Cylindrical Panel. ATAM 2011, 2, 171–181. [Google Scholar]
  21. Sakhr, J.; Chronik, B.A. Solving the Navier-Lamé Equation in Cylindrical Coordinates using the Buchwald Representation: Some Parametric Solutions with Applications. Adv. Appl. Math. Mech. 2018, 10, 1025–1056. [Google Scholar]
  22. Sakhr, J.; Chronik, B.A. Constructing Separable non-2π-Periodic Solutions to the Navier-Lamé Equation in Cylindrical Coordinates using the Buchwald Representation: Theory and Applications. Adv. Appl. Math. Mech. 2020, 12, 694–728. [Google Scholar]
  23. Morse, P.M.; Feshbach, H. Methods of Theoretical Physics; McGraw-Hill: New York, NY, USA, 1953; pp. 1764–1767. [Google Scholar]
  24. Buchwald, V.T. Rayleigh Waves in Transversely Isotropic Media. Q. J. Mech. Appl. Math. 1961, 14, 293–317. [Google Scholar] [CrossRef]
  25. Kund, T. Mechanics of Elastic Waves and Ultrasonic Nondestructive Evaluation; CRC Press: Boca Raton, FL, USA, 2019; pp. 76–78. [Google Scholar]
  26. Nezhad Mohammad, A.M.; Abdipour, P.; Bababeyg, M.; Noshad, H. Derivation of Green’s Function for the Interior Region of a Closed Cylinder. Amirkabir Int. J. Sci. Res. 2014, 46, 23–33. [Google Scholar]
  27. Eringen, A.C.; Suhubi, E.S. Fundamentals of Linear Elastodynamics. In Elastodynamics; Academic Press: New York, NY, USA, 1975; Volume II, pp. 343–366. [Google Scholar]
  28. Korenev, B.G. Bessel Functions and Application; Taylor & Francis: New York, NY, USA, 2002; p. 55. [Google Scholar]
  29. Hora, P.; Červená, O. Acoustic Emission Source Modeling. Appl. Comput. Mech. 2010, 4, 25–36. [Google Scholar]
  30. Dou, Y.; Luo, H.; Zhang, J. Elastic Properties of FeCr20Ni8Xn (X = Mo, Nb, Ta, Ti, V. W and Zr) Austenitic Stainless Steels: A First Principles Study. Metals 2019, 9, 145. [Google Scholar] [CrossRef] [Green Version]
  31. Boresi, A.P.; Chong, K.P.; Lee, J.D. Elasticity in Engineering Mechanics, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011; pp. 150–151. [Google Scholar]
Figure 1. Two forms of the point source vector along the xj and xz directions used in analytical modeling.
Figure 1. Two forms of the point source vector along the xj and xz directions used in analytical modeling.
Applsci 12 05272 g001
Figure 2. Spectra of three sets of k z values calculated as a function of η(lz0): (a) r0 = 0 ( x 0 i = x 0 j = 0 , z 0 = 1   m ) and (b) r0 = 0.25 m ( x 0 i = 0 , x 0 j = 0.25 ,   z 0 = 1   m ) .
Figure 2. Spectra of three sets of k z values calculated as a function of η(lz0): (a) r0 = 0 ( x 0 i = x 0 j = 0 , z 0 = 1   m ) and (b) r0 = 0.25 m ( x 0 i = 0 , x 0 j = 0.25 ,   z 0 = 1   m ) .
Applsci 12 05272 g002
Figure 3. Effect of the b parameter of the point source on the envelopes of the simulated u r j waveforms: (1) 6 × 104, (2) 8 × 104, (3) 1 × 105, (4) 2 × 105, and (5) 4 × 105.
Figure 3. Effect of the b parameter of the point source on the envelopes of the simulated u r j waveforms: (1) 6 × 104, (2) 8 × 104, (3) 1 × 105, (4) 2 × 105, and (5) 4 × 105.
Applsci 12 05272 g003
Figure 4. Displacements at the (0.5 m, 45°, 1 m) position generated by the point source of 1.0 N located at the center of the cylinder: (a,c) radial (black line), axial (red line), and tangential (green line) components; and (b,d) compression P (black line) wave and vertically polarized (SV, red line) and horizontally polarized (SH, green line) shear waves.
Figure 4. Displacements at the (0.5 m, 45°, 1 m) position generated by the point source of 1.0 N located at the center of the cylinder: (a,c) radial (black line), axial (red line), and tangential (green line) components; and (b,d) compression P (black line) wave and vertically polarized (SV, red line) and horizontally polarized (SH, green line) shear waves.
Applsci 12 05272 g004
Figure 5. Angular dependences of the maximum displacements at the (0.5 m, 45°, 1 m) position: radial (red line), axial (green line), and tangential (blue line) displacements generated by (a,c) Pj and (b,d) Pz of 1.0 N, located at (a,b) (0 m, 0 m, 1 m) and (c,d) (0 m, 0.25 m, 1 m) positions.
Figure 5. Angular dependences of the maximum displacements at the (0.5 m, 45°, 1 m) position: radial (red line), axial (green line), and tangential (blue line) displacements generated by (a,c) Pj and (b,d) Pz of 1.0 N, located at (a,b) (0 m, 0 m, 1 m) and (c,d) (0 m, 0.25 m, 1 m) positions.
Applsci 12 05272 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, K.B.; Kim, B.K.; Lee, S.G.; Kang, J.-G. Analytical Modeling of Acoustic Emission Due to an Internal Point Source in a Transversely Isotropic Cylinder. Appl. Sci. 2022, 12, 5272. https://doi.org/10.3390/app12105272

AMA Style

Kim KB, Kim BK, Lee SG, Kang J-G. Analytical Modeling of Acoustic Emission Due to an Internal Point Source in a Transversely Isotropic Cylinder. Applied Sciences. 2022; 12(10):5272. https://doi.org/10.3390/app12105272

Chicago/Turabian Style

Kim, Kwang Bok, Bong Ki Kim, Sang Guk Lee, and Jun-Gill Kang. 2022. "Analytical Modeling of Acoustic Emission Due to an Internal Point Source in a Transversely Isotropic Cylinder" Applied Sciences 12, no. 10: 5272. https://doi.org/10.3390/app12105272

APA Style

Kim, K. B., Kim, B. K., Lee, S. G., & Kang, J. -G. (2022). Analytical Modeling of Acoustic Emission Due to an Internal Point Source in a Transversely Isotropic Cylinder. Applied Sciences, 12(10), 5272. https://doi.org/10.3390/app12105272

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop