Stress Wave Propagation in a Rayleigh–Love Rod with Sudden Cross-Sectional Area Variations Impacted by a Striker Rod

This paper presents an in-depth study of the stress wave behavior propagating in a Rayleigh–Love rod with sudden cross-sectional area variations. The analytical solutions of stress waves are derived for the reflection and transmission propagation behavior at the interface of the cross-sectional area change in the rod, considering inertia and Poisson’s effects on the rod material. Examples solved using the finite element method are provided to verify the correctness of the analytical results. Based on the forward analysis of Rayleigh–Love wave propagation in a rod impacted by a striker rod, an impact-echo-type nondestructive testing (NDT) method is proposed to conduct defect assessment in rod-type structural components with sudden cross-sectional area changes within a cover medium. This proposed NDT method can identify the location, extension, and cross-sectional area drop ratios of an irregular zone in the rod to be inspected.


Introduction
Stress wave propagation, particularly concerning the nondestructive testing (NDT) of structures, is an important field in engineering.The NDT method is crucial for assessing the safety, reliability, and efficiency of structures [1][2][3][4][5] because it allows for the evaluation of structures and component integrity without causing damage.This approach is particularly vital for detecting defects or corrosive damage in engineering structures [6][7][8][9][10][11][12].Corrosive damage, often manifested as reductions in the cross-sectional area of rod-type elements, is a significant concern in structural engineering.This type of deterioration typically results from environmental exposures that accelerate the corrosion process, leading to the weakening of the material and potential structural failure.Detecting and quantifying such damage is critical for maintaining structural integrity and safety.The detection of the damage condition in rod-type structural elements embedded inside a medium is an important topic in the NDT field.Among the various theories introduced to understand the longitudinal wave propagation for the impact-echo NDT method, Rayleigh-Love rod theory [13] provides a very good model of the observed real behavior.This theory combines the effects of lateral inertia and Poisson's ratio on the rod material, which is often neglected in the traditional model of D'Alembert [14].Incorporating lateral inertia and Poisson's ratio is crucial for analyzing longitudinal stress wave propagation in rods because it considerably affects wave behavior, as evidenced by the research provided by Yang et al. [15,16].
The primary goal of this study was to develop an NDT method aimed at identifying the positions of cross-sectional area changes based on the forward analysis of stress wave propagation in a rod.Exhibiting abrupt changes in the cross-sectional area of a rod often indicates damage or defects in the rod-type structural member, like corrosion and bulge.To distinguish damage or defects from the effects of wave propagation modifications, reliance on changes in the amplitude of the wave is crucial; specifically, the amplitude values of reflected waves are indicative.The larger the reflected wave, the greater the damage or defect, as it indicates a discontinuity in the rod's cross-sectional area.This change in the waveform reveals the condition of the rod and helps identify the reasons for the discontinuities affecting the waves.The theoretical foundation of this research is grounded in Rayleigh-Love rod theory [13], which provides a comprehensive understanding of stress wave propagation in rods, particularly in the context of sudden cross-sectional changes.This analytical approach enables a detailed analysis of the complex dynamics of wave propagation, including the critical effects of lateral inertia and Poisson's ratio.
As shown in Figure 1, the model for this research was developed for using the impactecho method to determine the changing cross-sections A 2 and A 3 and length L 1 and l 2 based on stress wave propagation signals, where σ I is the incident stress.σ R1 and σ T1 are the reflected and transmitted stresses due to σ I at L 1 .σ R2 and σ T2 are the reflected and transmitted stresses due to σ T1 at L 2 .σ R3 and σ T3 are the reflected and transmitted stresses due to σ R2 at L 1 .The stress propagation process begins with the impact of the striker on the semi-infinite rod, generating an incident stress (σ I ) that travels along the semi-infinite rod.This incident stress (σ I ) propagates to location L 1 , where a change in cross-sectional area causes reflected stress (σ R1 ) and transmitted stress (σ T1 ).After the stress wave passes through the cross-section change at L 1 , the stress wave continues to move to the changing cross-section at L 2 .At this location, stress superposition occurs as the transmitted stress (σ T1 ) moves toward the cross-section at L 2 and the reflected stress (σ R1 ) travels back to the start of the rod.The transmitted stress (σ T1 ) will generate reflected stress (σ R2 ) and transmitted stress (σ T2 ) due to the change in cross-sectional area.The reflected wave (σ R2 ) will propagate back to position L 1 , generating a transmitted wave (σ R3 ) and a reflected wave (σ R3 ).A 1 , A 2 , and A 3 are the cross-sectional areas of segments 1, 2, and 3 of the semi-infinite rod, respectively.E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the rod, respectively.The striker rod is moved with an impact velocity of 2v 0 .Then, at the observed position, an incident wave is generated, and a reflection wave travels back due to a change in the cross-sectional area.A signal is measured at a sensing site for back-calculation to determine the position and severity of a cross-sectional area change.The analysis delves into the reflection and transmission in rods with double sudden changes in cross-sectional area.Such modifications lead to alterations in the waveform during propagation.These abrupt changes are often indicative of major defects, including cracks, necking, expansion, and corrosion within the rod.This aspect is crucial because the shape and intensity of the stress waves are significantly affected at the discontinuity interface within a structure [17,18].Previous studies [17,18] have documented the substantial influence of boundaries on wave behavior, particularly noting how sudden changes in cross-sectional area can modify wave propagation.However, these studies did not fully consider the specific effects of lateral inertia and Poisson's ratio on stress wave propagation.Results obtained from considering lateral inertia and Poisson's ratio show that the magnitude of stress increases, ranging from 23.84% to 30.94%, compared to results obtained without considering lateral inertia and Poisson's ratio.This higher stress value could cause inaccurate structure analysis results if not considered.In this study, we introduce a method to identify changes in the cross-sectional area and their locations along the rod, considering the effects of lateral inertia and Poisson's ratio.The stress values obtained are consistently higher than those calculated without considering the effects of lateral inertia and Poisson's ratio.This is an important consideration, as exceeding allowable stress values can lead to structural damage.

Stress Wave Propagation in a Rod
The behavior of stress waves in a material is typically described by the wave tion.It captures the relationship between the displacement of particles within the m and the rate at which these displacements change in the spatial-temporal domains dynamic phenomenon is often attributed to the transmission of energy and mome through the medium.The traditional wave equation for one-dimensional wave pro tion in a rod is given by [14]   The wave equation for one-dimensional wave propagation in a rod, which take account the effects of lateral inertia and Poisson's ratio as provided by Rayleigh-Lov is expressed as follows: , where u is the axial displacement in the axial direction of stress wave propagatio Poisson's ratio, and κ is the radius of gyration of the cross-section.The influence of κ to oscillations with larger amplitudes compared to scenarios without consideratio method for determining κ, related to the wave propagation equation, is detai Appendix B. In this theory, the stress equation can be expressed as Various methods, such as the finite difference method [19][20][21][22][23][24], the finite el method (FEM) [25,26], or Laplace transform [27][28][29], can be employed to solve Equ (2).The Laplace transform, a mathematical tool that transitions functions from th domain to the s domain, offers an elegant method for simplifying and solving differ equations [30].When it is applied in solving wave propagation, this transform co the original differential equation problem into an algebraic equation, making man and solving easy.The solution is then converted back to the time domain using the in Laplace transform, which provides an accurate solution to the original wave propa

Analysis of Stress Wave Propagation in a Rayleigh-Love Rod with Sudden Cross-Sectional Area Variations 2.1. Stress Wave Propagation in a Rod
The behavior of stress waves in a material is typically described by the wave equation.It captures the relationship between the displacement of particles within the material and the rate at which these displacements change in the spatial-temporal domains.This dynamic phenomenon is often attributed to the transmission of energy and momentum through the medium.The traditional wave equation for one-dimensional wave propagation in a rod is given by [14] as where u (x, t) represents the axial (longitudinal) displacement of a point located at position x in the rod at time t, and c 0 is the wave velocity without considering Poisson's ratio effect.Herein, wave velocity is related to the mass density ρ and Young's modulus E of the rod by using the equation c 0 = E/ρ.But considering the effect of Poisson's ratio, wave velocity (c) used in the stress wave propagation analysis is calculated as shown in Appendix A.
The wave equation for one-dimensional wave propagation in a rod, which takes into account the effects of lateral inertia and Poisson's ratio as provided by Rayleigh-Love [13], is expressed as follows: where u is the axial displacement in the axial direction of stress wave propagation, υ is Poisson's ratio, and κ is the radius of gyration of the cross-section.The influence of κ leads to oscillations with larger amplitudes compared to scenarios without consideration.
The method for determining κ, related to the wave propagation equation, is detailed in Appendix B. In this theory, the stress equation can be expressed as Various methods, such as the finite difference method [19][20][21][22][23][24], the finite element method (FEM) [25,26], or Laplace transform [27][28][29], can be employed to solve Equation (2).The Laplace transform, a mathematical tool that transitions functions from the time domain to the s domain, offers an elegant method for simplifying and solving differential equations [30].When it is applied in solving wave propagation, this transform converts the original differential equation problem into an algebraic equation, making managing and solving easy.The solution is then converted back to the time domain using the inverse Laplace transform, which provides an accurate solution to the original wave propagation problem.Therefore, the Rayleigh-Love rod model, which considers the effects of lateral inertia, can be effectively solved with the Laplace transform.This analytical approach contributes to a detailed interpretation of the wave propagation dynamics and the con-sequential stress distribution, thereby providing substantial insights into the dynamic response of materials under stress wave propagation.

Analysis of Rayleigh-Love Rod Impacted by a Striker Rod
Based on the results of the Laplace transform method derived from Yang et al. [16], the exact solution for the transmitted and reflected stress waves in a Rayleigh-Love rod with a sudden cross-sectional variation is derived first.Subsequently, a backward analysis method can be developed to determine the position of the changing cross-section and the ratio of the reduced area to the original area, as illustrated in Figure 2.
Sensors 2024, 24, 4230 4 problem.Therefore, the Rayleigh-Love rod model, which considers the effects of la inertia, can be effectively solved with the Laplace transform.This analytical approach tributes to a detailed interpretation of the wave propagation dynamics and the conseq tial stress distribution, thereby providing substantial insights into the dynamic resp of materials under stress wave propagation.

Analysis of Rayleigh-Love Rod Impacted by a Striker Rod
Based on the results of the Laplace transform method derived from Yang et al. the exact solution for the transmitted and reflected stress waves in a Rayleigh-Love with a sudden cross-sectional variation is derived first.Subsequently, a backward ana method can be developed to determine the position of the changing cross-section an ratio of the reduced area to the original area, as illustrated in Figure 2. In Figure 2, σI, σR, and σT represent the incident, reflected, and transmitted stre respectively.D1 is the diameter of the striker rod and rod segment 1. D2 is the diamet rod segment 2. A1 is the cross-sectional area of the striker rod and rod segment 1. A2 i cross-sectional area of rod segment 2. E1, ρ1, and υ1 are the Young's modulus, mass den and Poisson's ratio of the rod, respectively.
Yang et al. [16] focused on building an impact model based on split Hopkinson The process begins with a semi-infinite rod, known as the incident bar, being struck gitudinally by a striker bar.The traditional one-dimensional theory suggests that th sulting impact generates a rectangular pulse if both rods are composed of the same m rial and have identical cross-sectional areas.However, the distortions of the pulse expected on theoretical grounds and observed during experiments.Yang et al. [16 rived the analytical solutions for the longitudinal impact problems in dispersive Their discussion focused on a situation where a semi-infinite Rayleigh-Love rod wa pacted by a striker rod of length L with the same material and cross-sectional area a Rayleigh-Love rod, as shown in Figure 3, where σI is the incident stress, and D1 and A the diameter and the cross-sectional area of the striker rod and the semi-infinite rod spectively.E1, 1, and υ1 are the Young's modulus, mass density, and Poisson's ratio o two rods, respectively.2v0 is the impact velocity.In Figure 2, σ I , σ R , and σ T represent the incident, reflected, and transmitted stresses, respectively.D 1 is the diameter of the striker rod and rod segment 1. D 2 is the diameter of rod segment 2. A 1 is the cross-sectional area of the striker rod and rod segment 1.A 2 is the cross-sectional area of rod segment 2. E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the rod, respectively.
Yang et al. [16] focused on building an impact model based on split Hopkinson tests.The process begins with a semi-infinite rod, known as the incident bar, being struck longitudinally by a striker bar.The traditional one-dimensional theory suggests that the resulting impact generates a rectangular pulse if both rods are composed of the same material and have identical cross-sectional areas.However, the distortions of the pulses are expected on theoretical grounds and observed during experiments.Yang et al. [16] derived the analytical solutions for the longitudinal impact problems in dispersive rods.Their discussion focused on a situation where a semi-infinite Rayleigh-Love rod was impacted by a striker rod of length L with the same material and cross-sectional area as the Rayleigh-Love rod, as shown in Figure 3, where σ I is the incident stress, and D 1 and A 1 are the diameter and the cross-sectional area of the striker rod and the semi-infinite rod, respectively.E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the two rods, respectively.2v 0 is the impact velocity.
Sensors 2024, 24, 4230 4 problem.Therefore, the Rayleigh-Love rod model, which considers the effects of lat inertia, can be effectively solved with the Laplace transform.This analytical approach tributes to a detailed interpretation of the wave propagation dynamics and the consequ tial stress distribution, thereby providing substantial insights into the dynamic respo of materials under stress wave propagation.

Analysis of Rayleigh-Love Rod Impacted by a Striker Rod
Based on the results of the Laplace transform method derived from Yang et al. the exact solution for the transmitted and reflected stress waves in a Rayleigh-Love with a sudden cross-sectional variation is derived first.Subsequently, a backward anal method can be developed to determine the position of the changing cross-section and ratio of the reduced area to the original area, as illustrated in Figure 2. In Figure 2, σI, σR, and σT represent the incident, reflected, and transmitted stres respectively.D1 is the diameter of the striker rod and rod segment 1. D2 is the diamete rod segment 2. A1 is the cross-sectional area of the striker rod and rod segment 1. A2 is cross-sectional area of rod segment 2. E1, ρ1, and υ1 are the Young's modulus, mass den and Poisson's ratio of the rod, respectively.
Yang et al. [16] focused on building an impact model based on split Hopkinson t The process begins with a semi-infinite rod, known as the incident bar, being struck gitudinally by a striker bar.The traditional one-dimensional theory suggests that the sulting impact generates a rectangular pulse if both rods are composed of the same m rial and have identical cross-sectional areas.However, the distortions of the pulses expected on theoretical grounds and observed during experiments.Yang et al. [16] rived the analytical solutions for the longitudinal impact problems in dispersive r Their discussion focused on a situation where a semi-infinite Rayleigh-Love rod was pacted by a striker rod of length L with the same material and cross-sectional area as Rayleigh-Love rod, as shown in Figure 3, where σI is the incident stress, and D1 and A the diameter and the cross-sectional area of the striker rod and the semi-infinite rod spectively.E1, 1, and υ1 are the Young's modulus, mass density, and Poisson's ratio o two rods, respectively.2v0 is the impact velocity.At the initial time t = 0, a finite striker rod moving coaxially at an impact velocity of 2v 0 hits a stationary semi-infinite rod, with both rods assumed to be unstressed.Consequently, the initial conditions relating to the axial displacement u can be expressed as follows: Sensors 2024, 24, 4230 The boundary conditions at the free end, where x = −L, and at infinity are determined by If the two rods are in tight contact in a compressed state following the collision, then the coupled continuity conditions at the contact surface x = 0 can be established as Yang et al. [16] converted the values to dimensionless quantities for easy transformation by setting , and σ = σ E .The Laplace transform and inverse techniques are employed to solve Equation ( 2).The result of the displacement is shown as Equation (15a) in Yang et al.'s study [16].The stress state in a semi-infinite rod is then determined as follows: By performing the calculation and converting the dimensionless quantities from Equation ( 8) into dimensional quantities, the result of stress wave propagation in the semi-infinite rod can be obtained as follows: , η is a variable in the contour integration of the Laplace inverse transform, ρ is the mass density, c is the longitudinal wave velocity of the elastic material as expressed by Equation (A16) in Appendix A, and v 0 is the impact velocity.
In this study, the incident stress (σ I ) result of Equation ( 9) is considered to determine the reflected stress (σ R ) and transmitted stress (σ T ) in a Rayleigh-Love rod with a sudden change in cross-sectional area under the effects of Poisson's ratio and lateral inertia.Based on the balance of forces at the interface between two segments, the total force applied by segment 1 (from the incident and reflected waves) must be equal to the total force applied by segment 2 (from the transmitted wave), as shown in Figure 4, where v I , σ I , and u I are the incident velocity, incident stress, and incident displacement of segment 1, respectively.v R , σ R , and u R are the reflected velocity, reflected stress, and reflected displacement of segment 1, respectively.v T , σ T , and u T are the transmitted velocity, transmitted stress, and transmitted displacement of segment 2, respectively.A 1 and A 2 are the cross-sectional areas of segment 1 and segment 2 of the semi-infinite rod, respectively.E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the two rods, respectively.x I , x R , and x T are the positions determined for the incident, reflected, and transmitted stresses in the semi-infinite rod, respectively.x b is the position-sensing stress in the striker rod.Thus, the following equation is obtained: where A 1 and A 2 are the cross-sectional areas of segment 1 and segment 2, respectively, and σ I , σ R , and σ T are the incident, reflected, and transmitted stresses, respectively.Based on the continuity of velocity in Figure 4, the velocity at the interface of the two segments must be the same, which is expressed as transmitted stresses in the semi-infinite rod, respectively.xb is the position-sensing stress in the striker rod.Thus, the following equation is obtained: where A1 and A2 are the cross-sectional areas of segment 1 and segment 2, respectively, and σI, σR, and σT are the incident, reflected, and transmitted stresses, respectively.Based on the continuity of velocity in Figure 4, the velocity at the interface of the two segments must be the same, which is expressed as The relationship between stress and velocity in the rod is where v is the particle velocity of the material point under stress σ, and c is the longitudinal wave velocity of the elastic material (see Appendix A) and can be determined as , where Based on Equation ( 12), the velocity for incident, reflected, and transmitted waves can be, respectively, presented as Substituting Equations ( 14)-( 16) into Equation ( 11) can obtain Solving Equations ( 10) and ( 17) can obtain the reflected stress R The relationship between stress and velocity in the rod is where v is the particle velocity of the material point under stress σ, and c is the longitudinal wave velocity of the elastic material (see Appendix A) and can be determined as Based on Equation ( 12), the velocity for incident, reflected, and transmitted waves can be, respectively, presented as Substituting Equations ( 14)-( 16) into Equation ( 11) can obtain Solving Equations ( 10) and ( 17) can obtain the reflected stress σ R and transmitted stress σ T as where the reflected ratio at L 1 due to incident stress (σ I ) is where the transmitted ratio at L 1 due to incident stress (σ I ) is According to Figure 2, the stress wave propagation σ 1 traveling along segment 1 is determined as where Here, x ∈ [0, L 1 ], where x is the position determined for the incident and reflected stresses (see Figure 4).
The propagation of the transmitted stress propagation σ 2 traveling along segment 2 is determined as Here, x ∈ [L 1 , ∞], where x is the position determined for the transmitted stress.The stress wave propagation in the striker rod is determined as follows: Here, x ∈ [−L, 0], where x is the position for the determined stress wave in the striker rod.
When Poisson's ratio is set as zero, the incident stress, reflection stress, and transmitted stress are, respectively, expressed as follows: The numerical integration provided by Yang et al. [31] was utilized to obtain the stress wave propagation in a Rayleigh-Love rod to solve Equations ( 21)- (25).When the stress wave propagation does not consider Poisson's ratio and inertia effects in Equations ( 21)-( 25), they will degenerate to Equations ( 26)- (28), and the analytical solution can be obtained.21)- (25) with the FEM Solution This section presents the stress wave propagation in the rods consisting of two segments with different cross-sectional areas, where cross-sectional area A 2 is greater than A 1 , as shown in Figure 2. First, the accuracy of Equations ( 21)-( 25) will be confirmed by comparing the results derived from these equations and the results obtained from the FEM.The parameters of material properties and geometries of the striker rod and the segment rods used in this study are shown in Table 1.The geometrical shape and finite element mesh of the Rayleigh-Love rod model are plotted in Figure 5a,b, and the relevant parameters used in the Rayleigh-Love rod model are shown in Table 1.The elastic material wave velocity c of the rod and striker materials was computed as 5782.69m/s.In the Abaqus finite element analysis code, the solid elementtype C3D8R was employed, with 10,947 nodes and 8448 elements in the model, as shown in Figure 5b.The detailed stress variations are demonstrated in Figure 5c-h at the time from the initial state to 450 µs.
For the cross-sectional area ratio A 2 /A 1 = 2, i.e., A 2 = 1.41 × 10 −3 mm 2 , the comparisons of the stress responses over time at varying positions at the striker rod, segment 1, and segment 2 between both methods are illustrated in detail in Figures 6-8, respectively, through the whole wave propagation duration.A high consistency presents positive evidence for a high degree of match between the analytical solutions expressed by Equations ( 21)-( 25) and the FEM simulation results.This verifies the correctness of the numerical algorithm [31] used to solve Equations ( 21)-( 25) in this study.

Investigation of Stress Wave Propagation in a Rod at Various Positions
Figure 9 presents the two types of stress responses, incident and reflected stresses, at x = 0.5 m, i.e., point A in Figure 5b, in segment 1, within the rod based on Equation (21).The shape of the incident stresses remains unchanged in the period of 90-140 µs, as expected.By contrast, the reflected stress in the period of 260-310 µs significantly varies with cross-sectional area A 2 .In particular, an increase in cross-sectional area A 2 results in a corresponding increase in the magnitude of the reflected stresses.When the crosssectional areas are equal, i.e., A 1 = A 2 , no reflected stresses are observed, signifying the absence of reflected stresses.Moreover, the incident and reflected stresses exhibit the same compression stress responses.For the cross-sectional area ratio A2/A1 = 2, i.e., A2 = 1.41 × 10 mm , the compariso of the stress responses over time at varying positions at the striker rod, segment 1, an segment 2 between both methods are illustrated in detail in Figures 6, 7 and 8, respe tively, through the whole wave propagation duration.A high consistency presents po tive evidence for a high degree of match between the analytical solutions expressed b Equations ( 21)-( 25) and the FEM simulation results.This verifies the correctness of t numerical algorithm [31] used to solve Equations ( 21)-( 25) in this study.For the cross-sectional area ratio A2/A1 = 2, i.e., A2 = 1.41 × 10 −3 mm 2 , the comparison of the stress responses over time at varying positions at the striker rod, segment 1, and segment 2 between both methods are illustrated in detail in Figures 6, 7 and 8, respec tively, through the whole wave propagation duration.A high consistency presents posi tive evidence for a high degree of match between the analytical solutions expressed by Equations ( 21)-( 25) and the FEM simulation results.This verifies the correctness of th numerical algorithm [31] used to solve Equations ( 21)-( 25) in this study.

Reflected stress
Incident + reflected stresses   2 shows that the reflected stress values range from 23.84% to 24.55% for the cross-sectional area ratio, A 2 /A 1 , varying from two to five.This finding demonstrates that Poisson's effect has similar stress responses for different cross-sectional area ratios.2The shape of the incident stresses remains unchanged in the period of 90-140 μs, as e pected.By contrast, the reflected stress in the period of 260-310 μs significantly varies wi cross-sectional area A2.In particular, an increase in cross-sectional area A2 results in a co responding increase in the magnitude of the reflected stresses.When the cross-section areas are equal, i.e., A1 = A2, no reflected stresses are observed, signifying the absence reflected stresses.Moreover, the incident and reflected stresses exhibit the same compre sion stress responses.21The shape of the incident stresses remains unchanged in the period of 90-140 μs, as e pected.By contrast, the reflected stress in the period of 260-310 μs significantly varies wi cross-sectional area A2.In particular, an increase in cross-sectional area A2 results in a co responding increase in the magnitude of the reflected stresses.When the cross-section areas are equal, i.e., A1 = A2, no reflected stresses are observed, signifying the absence reflected stresses.Moreover, the incident and reflected stresses exhibit the same compre sion stress responses.

Transmitted stress
Incident stress Reflected stress   stress values when the cross-sectional area A 1 is equal to A 2 .In particular, the resulting transmitted stress responses become increasingly intensive as the cross-sectional area A 2 increases.Table 3 shows that the stress values have similar stress responses, with and without considering Poisson's effect for different cross-sectional area ratios, A 2 /A 1 .The waveforms of nonzero Poisson's ratio oscillate about the waveform of zero Poisson's value.Figure 10 illustrates the stress responses at x = 1.5 m, i.e., point B in Figure 5b, segment 2 within the rod based on Equation (24).The first peak values of this transmitte stress vary with cross-sectional area A2.The transmitted stresses are equal to the incide stress values when the cross-sectional area A1 is equal to A2.In particular, the resultin transmitted stress responses become increasingly intensive as the cross-sectional area A increases.Table 3 shows that the stress values have similar stress responses, with an without considering Poisson's effect for different cross-sectional area ratios, A2/A1.Th waveforms of nonzero Poisson's ratio oscillate about the waveform of zero Poisson value.Figure 11 illustrates the stress wave travel time history at various positions in a Rayleigh-Love rod.The stress wave propagation in the striker rod is presented at points I 1 , I 2 , I 3 , and P 0 (left).Points P 0 (right) to P 3 represent the presence of the incident and reflected waves.Points P 5 to P 7 describe the transmitted wave in the semi-infinite rod.The dotted lines represent stress in the striker rod, the solid lines represent the incident and reflected waves, and the dashed lines correspond to the transmitted wave.
The interface of the cross-sectional area jumps off the rod, and all three types of waves (incident, reflected, and transmitted waves) are observed at P 4 .The incident and reflected waves are detected on the left side of point P 4 , whereas the transmitted wave is discernible on the right side.The oscillations of the wave, influenced by Poisson's effect, are consistently found at each position, leading to stress values surpassing those from traditional wave equations that disregard Poisson's effect.This finding emphasizes the pertinence of the theoretical approach in aligning closely with the actual wave propagation observed in real tests and verifies its applicability to structural analysis works.The interface of the cross-sectional area jumps off the rod, and all three types of waves (incident, reflected, and transmitted waves) are observed at P4.The incident and reflected waves are detected on the left side of point P4, whereas the transmitted wave is discernible on the right side.The oscillations of the wave, influenced by Poisson's effect, are consistently found at each position, leading to stress values surpassing those from traditional wave equations that disregard Poisson's effect.This finding emphasizes the pertinence of the theoretical approach in aligning closely with the actual wave propagation observed in real tests and verifies its applicability to structural analysis works.

Example 2: Analysis of the Signal Obtained for the Rod Where A 1 ≥ A 2
This section presents stress wave propagation where cross-sectional area A 2 is less than cross-sectional area A 1 , as shown in Figure 12, where σ I , σ R , and σ T represent the incident, reflected, and transmitted stresses, respectively.D 1 is the diameter of the striker rod and rod segment 1. D 2 is the diameter of rod segment 2. A 1 is the cross-sectional areas of the striker rod and rod segment 1.A 2 is the cross-sectional area of rod segment 2. E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the rod, respectively.Except for the diameters of segment 2, the relevant parameters used in the Rayleigh-Love rod model are listed in Table 1.The diameters of segment 2 are set as 30, 21, 17, 15, and 13 mm instead.

Example 2: Analysis of the Signal Obtained for the Rod Where A1 ≥ A2
This section presents stress wave propagation where cross-sectional area A2 than cross-sectional area A1, as shown in Figure 12, where σI, σR, and σT represent cident, reflected, and transmitted stresses, respectively.D1 is the diameter of the rod and rod segment 1. D2 is the diameter of rod segment 2. A1 is the cross-sectiona of the striker rod and rod segment 1. A2 is the cross-sectional area of rod segment 2. and υ1 are the Young's modulus, mass density, and Poisson's ratio of the rod, respec Except for the diameters of segment 2, the relevant parameters used in the Rayleigh rod model are listed in Table 1.The diameters of segment 2 are set as 30, 21, 17, 15, mm instead.Figure 13 depicts the stress responses and incident and reflected stresses, at x = 0.5 m, i.e., point A in Figure 5b, in segment 1, within the rod based on Equation ( 21) when crosssectional area A 1 is greater than cross-sectional area A 2 .The shape of the incident stresses remains invariant in the 90-140 µs period regardless of the changes in cross-sectional area A 2 , as expected.Conversely, the reflected stress in the 260-310 µs period significantly varies with alterations in cross-sectional area A 2 .If A 2 is smaller than A 1 , this leads to intensive reflected stresses.When the cross-sectional areas are equal, i.e., A 1 = A 2 , no reflected stress responses are observed, signifying the absence of reflected stresses.In this scenario, the incident and reflected stresses exhibit opposite stress states, indicating that the reflected wave phase is strongly affected by the cross-sectional area ratio at the interface.Table 4 shows that the stress values have similar stress responses with and without Poisson's effect for different cross-sectional area ratios.Figure 14 illustrates the stress responses at x = 1.5 m, i.e., point B in Figure 5b, in se ment 2 within the rod based on Equation (24).When cross-sectional area A1 is greater tha cross-sectional area A2, it leads to more intensive transmitted stresses than the incide stress.When cross-sectional area A1 is equal to A2, the transmitted stresses are equal to th incident stress values.In particular, the resulting transmitted stress responses become increa ingly intensive as cross-sectional area A2 increases.Table 5 shows that the stress values hav similar stress responses, A2/A1, with and without considering Poisson's effect for differe cross-sectional area ratios.Figure 14 illustrates the stress responses at x = 1.5 m, i.e., point B in Figure 5b, in segment 2 within the rod based on Equation (24).When cross-sectional area A 1 is greater than cross-sectional area A 2 , it leads to more intensive transmitted stresses than the incident stress.When cross-sectional area A 1 is equal to A 2 , the transmitted stresses are equal to the incident stress values.In particular, the resulting transmitted stress responses become increasingly intensive as cross-sectional area A 2 increases.Table 5 shows that the stress values have similar stress responses, A 2 /A 1 , with and without considering Poisson's effect for different cross-sectional area ratios.In some engineering diagnosis problems, the damaged position and the condition the rod inside a cover medium must be identified.The impact-echo method is usual applied for this purpose.The variations in the reflected and transmitted stresses can presented in terms of cross-sectional area ratio A2/A1 to illustrate the influence of the cros sectional area change in the rod.The reflected stress  R and transmitted stress  T c be expressed as follows by letting α = A2/A1 and substituting it into Equation (18) an Equation ( 19):

Table 4. Reflected stresses in case of
where α = A2/A1 is the cross-sectional area ratio.In some engineering diagnosis problems, the damaged position and the condition of the rod inside a cover medium must be identified.The impact-echo method is usually applied for this purpose.The variations in the reflected and transmitted stresses can be presented in terms of cross-sectional area ratio A 2 /A 1 to illustrate the influence of the cross-sectional area change in the rod.The reflected stress σ R and transmitted stress σ T can be expressed as follows by letting α = A 2 /A 1 and substituting it into Equations ( 18) and ( 19): where α = A 2 /A 1 is the cross-sectional area ratio.Equations ( 29) and (30) show the relationship between the cross-sectional area ratio and the reflected and transmitted stress.When the cross-sectional area ratio A 2 /A 1 approaches 0, it follows that σ R approaches −σ I , and σ T approaches 2σ I ; as A 2 /A 1 approaches ∞, it follows that σ R approaches σ I , and σ T approaches 0.
Suppose that one knows the first peak stress values, labeled σ p I as the peak incident stress and σ p R as the peak reflected stress.In this case, one can determine the change in the cross-sectional area A 2 from the acquired signal.From Equation ( 29), the cross-sectional area ratio α can be rewritten as Assume that one measures the wave traveling time t, which is the interval of the first peak of the wave traveling from the sensing point and the first peak of the reflected wave received at the same sensing point.In this case, one can determine the length of L m between the sensing point and the reflected interface as where t is the time duration from the first peak stress to the second peak stress, and c 1 is the longitudinal wave speed of the rod material.The sensing point is located at a distance L s from the rod end that is smaller than the embedded length L e , as shown in Figure 15.
Suppose that one knows the first peak stress values, labeled p I  as the peak incide stress and p R  as the peak reflected stress.In this case, one can determine the change the cross-sectional area A2 from the acquired signal.From Equation ( 29), the cross-section area ratio α can be rewritten as ) Assume that one measures the wave traveling time t, which is the interval of the fir peak of the wave traveling from the sensing point and the first peak of the reflected wa received at the same sensing point.In this case, one can determine the length of Lm b tween the sensing point and the reflected interface as where t is the time duration from the first peak stress to the second peak stress, and c1 the longitudinal wave speed of the rod material.The sensing point is located at a distan Ls from the rod end that is smaller than the embedded length Le, as shown in Figure We assume that We assume that As shown in Figure 16a, Case A indicates the stress responses at the position of x = L s = 0.2 m by using Equation (21).The relevant parameters used in the Rayleigh-Love rod model are identical to those used in Example 1, i.e., Table 1.The diameter of segment 2 is chosen as 67 mm.The first peak incident stress is σ As shown in Figure 16a, Case A indicates the stress responses at the position of x = Ls = 0.2 m by using Equation (21).The relevant parameters used in the Rayleigh-Love rod model are identical to those used in Example 1, i.e., Table 1.The diameter of segment 2 is chosen as 67 mm.The first peak incident stress is 327.257MPa , and the first peak reflected stress is 220.398MPa , leading to a reflected ratio of 0.67

RI
. Figure 16 shows the corresponding cross-sectional area ratio

Case B: Cross-sectional area ratio
As shown in Figure 16b, Case B indicates the stress responses at the position of x = L s = 0.2 m by using Equation (21).The relevant parameters of materials are identical to those used in Case A except for the diameter of segment 2. The diameter of segment 2 is set as 13 mm.The first peak incident stress is σ   6, the length of the striker rod varies between 0.15 m and 0.075 m, respectively.Velocities were used for the cases as 2.9 m/s, 5.8 m/s, and 11.6 m/s, respectively.The impact energy of each case that the striker applied on the rod is also listed in Table 6.Stress responses for the different cases of impact velocity and striker length combination are shown in Figure 17.This figure demonstrates that the first peak stress values varied with the impact energy of the striker input to the rod.It is noted that the time interval t between the first peak of the incident and reflected wave signal is fixed as t = 282.3µs, and the stress ratio obtained between the cases of reflected and incident waves is the same for all cases, as calculated in Equation (33).
Sensors 2024, 24, 4230 20 of 31 sensing point.The measurements of the impact velocity and the length of the striker are not required.The location and cross-sectional area ratio of the rod with cross-sectional area jump can be identified through the curves in Figure 15.15.This case study shows that the impact-echo test method proposed in this paper provides suitable impact energy to the inspected rod and measures the response signal at a sensing point.The measurements of the impact velocity and the length of the striker are not required.The location and cross-sectional area ratio of the rod with cross-sectional area jump can be identified through the curves in Figure 15.

Reflection and Transmission of Stress Wave in a Rod with Double Sudden Cross-Sectional Area Variations
This section discusses a rod with double cross-sectional area variations.Initially, the cross-section is A 1 , which then changes to A 2 and A 3 at positions L 1 and L 2 , respectively.Using the model depicted in Figure 18, a formula is developed to determine the signal received at the sensor location.Based on this formula, changes in the cross-sectional area and the positions of these changes can be identified.A stress wave is generated by a striker, and upon encountering a change in the cro section, it causes a reflected wave and a transmitted wave, as described in Section There are two locations of sudden cross-section change, as illustrated in Figure 18, wh σI is the incident stresses.σR1 and σT1 are the reflected and transmitted stresses due to σ x = L1.σR2 and σT2 are the reflected and transmitted stresses due to σT1 at x = L2.σR3 and are the reflected and transmitted stresses due to σR2 at L1. A1, A2, and A3 are the cro sectional areas of segments 1, 2, and 3 of the semi-infinite rod, respectively.E1, 1, and are the Young's modulus, mass density, and Poisson's ratio of the rod, respectively.The fore, the stress wave at the sensor location is determined as follows: is the reflected stresses obtained at the sensor position due to interface at L2 and determined as follows: ( ) where R is the reflected ratio due to interface at L1 and L2, and the method for determin A stress wave is generated by a striker, and upon a change in the crosssection, it causes a reflected wave and a transmitted wave, as described in Section 2.2.There are two locations of sudden cross-section change, as illustrated in Figure 18, where σ I is the incident stresses.σ R1 and σ T1 are the reflected and transmitted stresses due to σ I at x = L 1 .σ R2 and σ T2 are the reflected and transmitted stresses due to σ T1 at x = L 2 .σ R3 and σ T3 are the reflected and transmitted stresses due to σ R2 at L 1 .A 1 , A 2 , and A 3 are the cross-sectional areas of segments 1, 2, and 3 of the semi-infinite rod, respectively.E 1 , ρ 1 , and υ 1 are the Young's modulus, mass density, and Poisson's ratio of the rod, respectively.Therefore, the stress wave at the sensor location is determined as follows: where σ I (x, t) is the incident stress as shown in Equation ( 22); σ R1 (2L 1 − x, t) is the reflected stresses due to interface at L 1 as shown in Equation ( 23); and σ T3 (2L 2 − x, t) is the reflected stresses obtained at the sensor position due to interface at L 2 and is determined as follows: where R T3 is the reflected ratio due to interface at L 1 and L 2 , and the method for determining the reflection ratio is similar to that described in Section 2.2, calculated as follows: A 2 be the cross-sectional area ratio; substituting them into Equation (36) and combining with Equation (35) obtains β, as follows: To verify the signal received at the sensor location, determined by Equation ( 38), a comparison with the FEM is illustrated in Figure 19.The parameters of the material properties and geometries used are listed in Table 1.Additionally, there is a modification where A 1 = A 3 = 2A 2 and the length L of the striker rod is 0.1 m.The results obtained indicate that there is good agreement between the analytical solution and the FEM results.The general steps of the processed NDT method to identify defect sizes in the rod (see Figure 20) are as follows: 1. First, use a striker rod with a length of L to impact the semi-infinite rod with an initial velocity of 2v0 to generate a stress wave propagating in the rod.2. Place a strain sensor at position x = Ls on the semi-infinite rod to receive the deformation signal generated by the striker rod.The general steps of the processed NDT method to identify defect sizes in the rod (see Figure 20) are as follows: 1.
First, use a striker rod with a length of L to impact the semi-infinite rod with an initial velocity of 2v 0 to generate a stress wave propagating in the rod.

2.
Place a strain sensor at position x = L s on the semi-infinite rod to receive the deformation signal generated by the striker rod.

3.
The first part of the signal received at position x = L s is the incident stress, with a peak stress value of σ p I and a duration of T i .The wave then oscillates around zero. 4.
The second part of the signal obtained is the reflected stress σ p R1 and a time interval from the first stress peak to the second stress peak of t 1 .From this signal, A 2 and L 1 can be determined as follows: the first peak incident stress is σ p I , and the first peak reflected stress is σ p R1 , leading to a reflected ratio of σ p R1 /σ p I .Figure 15 shows the corresponding cross-sectional area ratio A 2 /A 1 = α (or A 2 = αA 1 ).The received signal of the reflected stress at t 1 corresponds to L m in Figure 15.The location of the first cross-sectional area variation, L 1 , is calculated as L m + L s . 5.
The third part of the signal obtained is the reflected stress σ p T3 at a time interval t 2 from the second stress peak to the third stress peak.From this signal, A 3 and L 2 can be determined as follows: the first peak reflected stress σ  As shown in Figure 22 and Table 7, the signals are received at the sensor position with various length of l 2 (area reduction zone).As the length of l 2 increases, the time t 2 to receive the reflected signal also increases.The main wavelength λ of the stress impulse of duration T i is equal to where c is the wave speed of the semi-infinite rod material and c s is the wave speed of the striker rod material.From Equation (38), we know that the wavelength λ can be less than 2L when the striker rod material has a higher Young's modulus than the semi-infinite rod material.When c = c s , the wavelength λ is equal to 2L.It is known from the NDT theory that the defect of size l 2 can be identified only by wavelength λ if the impact wave is smaller than l 2 , i.e., λ = 2L < l 2 .The damaged length l 2 is determined as When l 2 is less than or equal to 2L (2 × 0.1 m = 0.2 m), the received signal is noisy due to the effects of lateral inertia and Poisson's ratio; this is indicated in Table 7 with a large error (>5%), making it difficult to accurately determine l 2 .This result highlights the significant influence of inertia and Poisson's ratio on stress wave propagation in the rod, which causes noise in the signal when the damaged length l 2 is small, thus complicating the determination of changes in length and cross-sectional area.From Equation (38), it is known that the wavelength λ of the detecting incident wave is affected by the striker rod length L and the wave speed ratio of the rod and striker materials.The minimum defect length l 2 that can be identified is equal to the wavelength λ of the incident wave.The result is confirmed through the comparison in Table 7 of the calculated lengths (L 1 , l 2 ) and cross-sectional areas (A 2 , A 3 ) using the Rayleigh-Love theory with the real lengths and real cross-sectional areas.
where c is the wave speed of the semi-infinite rod material and cs is the wave speed of the striker rod material.The examples presented above aim to determine the change in cross-sectional area and its position in a rod based on the Rayleigh-Love rod theory.The model remains theoretical; it assumes the material to be linearly elastic, homogeneous, isotropic, and uniformly varying in cross-section, which may differ from real-world conditions.This theoretical work, however, establishes a foundation for understanding the behavior of stress wave propagation in Rayleigh-Love rods with abrupt variations in the cross-sectional area, considering lateral inertia and Poisson's ratio.This understanding is crucial for developing accurate nondestructive testing methods applicable in real-world scenarios.Based on this research, further developments will be pursued with more complex material models and intricate geometric changes.This approach will enhance the applicability of our method for nondestructive testing (NDT) of rod-type structural elements under more realistic conditions.

Conclusions
This paper presents a detailed study of stress wave propagation in a Rayleigh-Love rod characterized by a sudden change in the cross-sectional area after a distance from the impacted end of the rod.Based on the analytical solution by Yang et al. [16], the analytical solutions of the transmitted and reflected stresses in a Rayleigh-Love rod with sudden cross-section variation were obtained.This study highlights the influence of a sudden cross-sectional variation on wave behavior, triggering reflections and transmissions at the interface of discontinuity.The role of Poisson's effect is emphasized, demonstrating how it modulates wave behavior in the rod with discontinuities.Examples solved using the FEM verify the correctness of the modeling and numerical algorithm in terms of the analytical results of this study.
In addition to the forward analysis of Rayleigh-Love wave propagation in rods impacted by a striker rod with the same impedance as the rod to be inspected, an impact-echotype NDT method is proposed to assess the condition of rod-type structural components with sudden cross-sectional area changes within a cover medium based on the measured signal at the measurable zone of the rod to be inspected.
The investigation results show that when determining the cross-sectional area using the σ R /σ I and σ T3 /σ I ratios, the decision to consider or not consider Poisson's ratio has little influence on the final ratio (Figures 15 and 21).Including Poisson's ratio leads to a similar rate of increase in amplitudes in the incident and reflected stress waves.As a result, the reflected ratio σ p R /σ p I remains nearly the same, regardless of whether Poisson's ratio is accounted for or not.However, in practical applications, the first peak σ values can be easily determined through measured signals.The technique proposed in this study can be easily applied in the field to conduct nondestructive evaluation of the location L 1 , extension l 2 , and ratios of cross-sectional area changes (A 2 /A 1 , A 3 /A 2 ) of a defect in a rod within a cover medium.

Figure 1 .
Figure 1.Model for determining the changing cross-sectional areas A2 and A3 and the length L2 (area reduction zone) in a cover medium (x > Le) based on stress wave propagation theory as x, t) represents the axial (longitudinal) displacement of a point located at pos in the rod at time t, and 0 c is the wave velocity without considering Poisson's ratio Herein, wave velocity is related to the mass density  and Young's modulus E of the r using the equation 0 / cE  = .But considering the effect of Poisson's ratio, wave v (c) used in the stress wave propagation analysis is calculated as shown in Appendix

Figure 1 .
Figure 1.Model for determining the changing cross-sectional areas A 2 and A 3 and the length L 1 and L 2 (area reduction zone) in a cover medium (x > L e ) based on stress wave propagation theory.

Figure 2 .
Figure 2. Scheme of a Rayleigh-Love rod with a sudden cross-sectional area change within a medium of x > Le, impacted by a striker rod.

Figure 3 .
Figure 3. Scheme of a semi-infinite rod under the impact of a striker rod of the same materia cross-sectional area.

Figure 2 .
Figure 2. Scheme of a Rayleigh-Love rod with a sudden cross-sectional area change within a cover medium of x > L e , impacted by a striker rod.

Figure 2 .
Figure 2. Scheme of a Rayleigh-Love rod with a sudden cross-sectional area change within a c medium of x > Le, impacted by a striker rod.

Figure 3 .Figure 3 .
Figure 3. Scheme of a semi-infinite rod under the impact of a striker rod of the same material cross-sectional area.

Figure 4 .
Figure 4.The interface between the two segments of a Rayleigh-Love rod with the same material but different cross-sectional areas.

Figure 4 .
Figure 4.The interface between the two segments of a Rayleigh-Love rod with the same material but different cross-sectional areas.

1 .
Example 1: Study of the Stress Wave Propagation in a Rod with A 1 ≤ A 2 3.1.1.Verification of the Analytical Solution from Equations (

Figure 5 .Figure 5 .
Figure 5. Stress wave propagating in a Rayleigh-Love rod model using solid elements for crosssectional area ratio A2/A1 = 2 analyzed using the finite element analysis code in Abaqus.

Figure 6 .
Figure 6.Comparison of the stress wave propagation in the striker rod between the analytical so tion and FEM results for the cross-sectional area ratio A2/A1 = 2.

Figure 7 .
Figure 7.Comparison of the stress wave propagation in segment 1 between the analytical soluti and FEM results for the cross-sectional area ratio A2/A1 = 2.

Figure 6 .
Figure 6.Comparison of the stress wave propagation in the striker rod between the analytical solution and FEM results for the cross-sectional area ratio A 2 /A 1 = 2.

Figure 6 .
Figure 6.Comparison of the stress wave propagation in the striker rod between the analytical solu tion and FEM results for the cross-sectional area ratio A2/A1 = 2.

Figure 7 .
Figure 7.Comparison of the stress wave propagation in segment 1 between the analytical solutio and FEM results for the cross-sectional area ratio A2/A1 = 2.

Figure 7 .
Figure 7.Comparison of the stress wave propagation in segment 1 between the analytical solution and FEM results for the cross-sectional area ratio A 2 /A 1 = 2.

Figure 9
Figure 9 also shows the stress responses without considering Poisson's effect.The peak stress values of the incident stress wave are −333.24and −263.286MPa, respectively.This finding indicates that the stress under Poisson's effect increases by 25.57% compared with that without Poisson's effect.Table2shows that the reflected stress values range from 23.84% to 24.55% for the cross-sectional area ratio, A 2 /A 1 , varying from two to five.This finding demonstrates that Poisson's effect has similar stress responses for different cross-sectional area ratios.

Figure 8 .
Figure 8.Comparison of the stress wave propagation in segment 2 between the analytical soluti and FEM results for the cross-sectional area ratio A2/A1 = 2.3.1.2.Investigation of Stress Wave Propagation in a Rod at Various PositionsFigure9presents the two types of stress responses, incident and reflected stresses, x = 0.5 m, i.e., point A in Figure5b, in segment 1, within the rod based on Equation (2The shape of the incident stresses remains unchanged in the period of 90-140 μs, as e pected.By contrast, the reflected stress in the period of 260-310 μs significantly varies wi cross-sectional area A2.In particular, an increase in cross-sectional area A2 results in a co responding increase in the magnitude of the reflected stresses.When the cross-section areas are equal, i.e., A1 = A2, no reflected stresses are observed, signifying the absence reflected stresses.Moreover, the incident and reflected stresses exhibit the same compre sion stress responses.

Figure 9 .
Figure 9. Stress wave propagation for the incident and reflected stresses in some cases of A2 ≥ observed at x = 0.5 m (impact velocity 2v0 = 11.6 m/s).

Figure 8 .
Figure 8.Comparison of the stress wave propagation in segment 2 between the analytical solution and FEM results for the cross-sectional area ratio A 2 /A 1 = 2.

Figure 8 .
Figure 8.Comparison of the stress wave propagation in segment 2 between the analytical solutio and FEM results for the cross-sectional area ratio A2/A1 = 2.3.1.2.Investigation of Stress Wave Propagation in a Rod at Various PositionsFigure9presents the two types of stress responses, incident and reflected stresses, x = 0.5 m, i.e., point A in Figure5b, in segment 1, within the rod based on Equation (21The shape of the incident stresses remains unchanged in the period of 90-140 μs, as e pected.By contrast, the reflected stress in the period of 260-310 μs significantly varies wi cross-sectional area A2.In particular, an increase in cross-sectional area A2 results in a co responding increase in the magnitude of the reflected stresses.When the cross-section areas are equal, i.e., A1 = A2, no reflected stresses are observed, signifying the absence reflected stresses.Moreover, the incident and reflected stresses exhibit the same compre sion stress responses.

Figure 9 .
Figure 9. Stress wave propagation for the incident and reflected stresses in some cases of A2 ≥ observed at x = 0.5 m (impact velocity 2v0 = 11.6 m/s).

Figure 9 .
Figure 9. Stress wave propagation for the incident and reflected stresses in some cases of A 2 ≥ A 1 observed at x = 0.5 m (impact velocity 2v 0 = 11.6 m/s).

Figure 10
Figure 10 illustrates the stress responses at x = 1.5 m, i.e., point B in Figure 5b, in segment 2 within the rod based on Equation (24).The first peak values of this transmitted stress vary with cross-sectional area A 2 .The transmitted stresses are equal to the incident

Figure 11
Figure11illustrates the stress wave travel time history at various positions in a Rayleigh-Love rod.The stress wave propagation in the striker rod is presented at points I1, I2, I3, and P0 (left).Points P0 (right) to P3 represent the presence of the incident and reflected waves.Points P5 to P7 describe the transmitted wave in the semi-infinite rod.The dotted lines represent stress in the striker rod, the solid lines represent the incident and reflected waves, and the dashed lines correspond to the transmitted wave.

Figure 11 .
Figure 11.Stress wave travel time history at various positions in the striker rod and Rayleigh-Love rod with a sudden cross-sectional area variation for the incident, reflected, and transmitted waves depending on position and time.

Figure 11 .
Figure 11.Stress wave travel time history at various positions in the striker rod and Rayleigh-Love rod with a sudden cross-sectional area variation for the incident, reflected, and transmitted waves depending on position and time.

Figure 12 .
Figure 12.Scheme of a Rayleigh-Love rod with a sudden cross-sectional area change from l small embedded in a cover medium of x > Le.

Figure 13 Figure 12 .
Figure13depicts the stress responses and incident and reflected stresses, at m, i.e., point A in Figure5b, in segment 1, within the rod based on Equation (21) cross-sectional area A1 is greater than cross-sectional area A2.The shape of the in stresses remains invariant in the 90-140 μs period regardless of the changes in cro tional area A2, as expected.Conversely, the reflected stress in the 260-310 μs perio

Figure 13 .
Figure 13.Stress wave propagation for the incident and reflected waves in some cases of A2 ≤ observed at x = 0.5 m (impact velocity of 2v0 = 11.6 m/s).

Figure 13 .
Figure 13.Stress wave propagation for the incident and reflected waves in some cases of A 2 ≤ A 1 observed at x = 0.5 m (impact velocity of 2v 0 = 11.6 m/s).

Figure 14 .
Figure 14.Stress wave propagation for transmitted waves in some cases of A2 ≤ A1 observed at x 1.5 m (impact velocity of 2v0 = 11.6 m/s).

Figure 14 .
Figure 14.Stress wave propagation for transmitted waves in some cases of A 2 ≤ A 1 observed at x = 1.5 m (impact velocity of 2v 0 = 11.6 m/s).

Figure 15 .
Figure 15.Diagram for determining the length Lm and cross-sectional area A2.

1 Figure 15 .
Figure 15.Diagram for determining the length L m and cross-sectional area A 2 .

Figure 15
shows the relation curve (labeled as a red line) between the cross-sectional area ratio A 2 /A 1 and the reflected ratio σ p R /σ p I and the relation curve (labeled as a blue straight line) between the length L m between the sensing point and the reflected interface and the wave traveling time t.σ p R /σ p I = −1 implies that A 2 = 0, indicating that the reflected wave is equal in magnitude but opposite in sign to the incident wave.It also represents a free surface at position L 1 .σ p R /σ p I = 1 implies that A 2 → ∞ , indicating that the reflected wave is equal to the incident wave in magnitude.When σ p R /σ p I = 0, this scenario implies that A 2 /A 1 = 1.If no reflected wave is received, the cross-sectional area in segment 2 remains unchanged, i.e., A 2 = A 1 .For −1 < σ p R /σ p I < 0, the cross-sectional area A 2 is smaller than A 1 .When 0 < σ p R /σ p I < 1, the cross-sectional Sensors 2024, 24, 4230 17 of 29 area A 2 is greater than A 1 .The length of L 1 = L m + L s can be determined from Figure 15.The application in Figure 15 is demonstrated in the following subsection.4.1.2.Examples: Determination of Cross-Sectional Area A 2 and Length L 1 Based on Signals Obtained with Changing Cross-Sectional Area Ratios and Impact Velocities 1. Case A: Cross-sectional area ratio A 2 /A 1 = 5

pI 1 RI = im- plies that 2 A
= −327.257MPa, and the first peak reflected stress is σ p R = −220.398MPa, leading to a reflected ratio of σ p R /σ p I = 0.67. Figure 16 shows the corresponding cross-sectional area ratio α = A 2 /A 1 = 5, which implies A 2 = 5A 1 .The signal received by the reflected stress at time t is 458.9 µs, indicating that L m = 1.3 m, from the corresponding curve in Figure 15.Given that L s = 0.2 m is chosen, L 1 is calculated as L m + L s , equal to 1.5 m. sign to the incident wave.It also represents a free surface at position L1.→ , indicating that the reflected wave is equal to the incident wave in mag- nitude.When 0 pp RI = , this scenario implies that A2/A1 = 1.If no reflected wave is re- ceived, the cross-sectional area in segment 2 remains unchanged, i.e., A2 = A1.cross-sectional area A2 is smaller than A1.When 01 pp RI   , the cross-sectional area A2 is greater than A1.The length of L1 = Lm + Ls can be determined from Figure 15.The application in Figure 15 is demonstrated in the following subsection.4.1.2.Examples: Determination of Cross-Sectional Area A2 and Length L1 Based on Signals Obtained with Changing Cross-Sectional Area Ratios and Impact Velocities 1. Case A: Cross-sectional area ratio A2/A1 = 5 , which implies A2 = 5A1.The signal received by the reflected stress at time t is 458.9 μs, indicating that Lm = 1.3 m, from the corresponding curve in Figure15.Given that Ls = 0.2 m is chosen, L1 is calculated as Lm + Ls, equal to 1.5 m.

Figure 16 .Figure 16 .
Figure 16.The cross-sectional area A2 and length L1 determined based on the stress signal measured at x = Ls: (a) A2 > A1; (b) A1 > A2. 2. Case B: Cross-sectional area ratio A2/A1 = 0.2 As shown in Figure 16b, Case B indicates the stress responses at the position of x = Ls = 0.2 m by using Equation (21).The relevant parameters of materials are identical to those used in Case A except for the diameter of segment 2. The diameter of segment 2 is set as

pI
= −327.257MPa, and the first peak reflected stress is σ p R = 220.411MPa, indicating the reflected ratio of σ p R /σ p I = −0.67.
Figure 16   shows the corresponding cross-sectional area ratio A 2 /A 1 = 0.2; thus, A 2 = A 1 /5.The received signal of the reflected stress at t = 282.8µs corresponds to L m = 0.8 m in Figure15.Given that L s = 0.2 m is chosen, L 1 is calculated as L m + L s , equal to 1.0 m.3.Case C: Influences of different impact velocities and striker lengths on the characteristics of response signalsAll parameters of Case C are identical to those of Case B, except for the different impact velocities and striker lengths.As shown in Table

Figure 17 .
Figure 17.Determination of the cross-sectional area A2 and length L1 based on the stress signal measured with different velocities and striker lengths at the position of x = Ls, A1 > A2, α = 0.2.

12mvFigure 17 .
Figure 17.Determination of the cross-sectional area A 2 and length L 1 based on the stress signal measured with different velocities and striker lengths at the position of x = L s , A 1 > A 2 , α = 0.2.

Figure 18 .
Figure 18.Scheme of a Rayleigh-Love rod impacted by a striker rod with double sudden cr sectional area variations in a cover medium of x > Le.

Figure 18 .
Figure 18.Scheme of a Rayleigh-Love rod impacted by a striker rod with double sudden crosssectional area variations in a cover medium of x > L e .

Figure 19 .
Figure 19.Comparison of stress wave propagation in segment 1 between the analytical solution and FEM results, observed at x = 0.5 m and x = 0.8 m.

3 .Figure 19 .
Figure 19.Comparison of stress wave propagation in segment 1 between the analytical solution and FEM results, observed at x = 0.5 m and x = 0.8 m.

p
T3 leads to a reflected ratio of σ p T3 /σ p I .Based on Figure 21, the corresponding cross-sectional area ratio is A 3 /A 2 = β with A 2 = αA 1 and A 3 = αβA 1 .The received signal of the reflected stress at t 2 corresponds to l 2 in Figure 21.Sensors 2024, 24, 4230 23 of 31

Figure 22 .Figure 22 .
Figure 22.The signal received at the sensor location when the length l2 varies.

Table 1 .
The parameters of material properties and geometries of the striker rod and the segment rods used in example 1.

. Identification of the Locations of Sudden Cross-Sectional Area Change in the Rod 4
.1.Reflection and Transmission of Stress Wave in a Rod with Single Sudden Cross-Sectional Area Variation 4.1.1.Developing a Formula to Determine the Damaged Zone

. Identification of the Locations of Sudden Cross-Sectional Area Change in the Rod 4
.1.Reflection and Transmission of Stress Wave in a Rod with Single Sudden Cross-Sectional Area Variation 4.1.1.Developing a Formula to Determine the Damaged Zone

Table 6 .
The initial impact kinetic energy ( )

Table 6 .
The initial impact kinetic energy 1/2mv 2 0 of the striker in each case in Figure17., they all have the same length L m and cross-sectional area A 2 as Figure )

Table 7 .
Comparison of calculated length (L 1 , l 2 ) and cross-sectional areas (A 2 , A 3 ) using Rayleigh-Love theory with real lengths and real cross-sectional areas.No.

Table 7 .
Comparison of calculated length (L1, l2) and cross-sectional areas (A2, A3) using Rayleigh-Love theory with real lengths and real cross-sectional areas.