Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Reliability Optimization Design of Constrained Metamorphic Mechanism Based on the Augmented Assur Groups
Previous Article in Journal
Quantitative Research on Roof Deformation and Temporary Support Stiffness in Deep-Mine Gob-Side Entry Retaining by Roof Cutting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stress Wave Propagation in a Semi-Infinite Rayleigh–Love Rod under the Collinear Impact of a Striker Rod with Different General Impedances

by
Nguyen Ngoc Thang
* and
Chung-Yue Wang
Department of Civil Engineering, National Central University, No. 300, Chungda Road, Chung-Li District, Taoyuan 32001, Taiwan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(15), 6523; https://doi.org/10.3390/app14156523
Submission received: 30 June 2024 / Revised: 21 July 2024 / Accepted: 23 July 2024 / Published: 26 July 2024

Abstract

:
This paper studies elastic stress wave propagation generated by an impact in a system consisting of a moving striker rod and an initially stationary semi-infinite rod. This research emphasizes the role of different general impedances in affecting the response during the wave propagation process. The Rayleigh–Love rod theory is used in this research to consider lateral inertia and Poisson’s effects on longitudinal waves in rods, as these factors lead to greater stress results compared to the traditional wave equation. To solve the complex wave equations in a Rayleigh–Love rod with different general impedances, the numerical inversion of Laplace transformation is applied and verified using the results of previous research. This study demonstrates that variations in general impedances cause different wave reflection and transmission behaviors at the interface. As the wave interacts with this discontinuity of impedance, it may be amplified or attenuated, and changes in impedance can significantly affect wave propagation behaviors.

1. Introduction

The study of stress wave propagation is a fundamental aspect of materials science and engineering, with implications for the design, analysis, and application of mechanical systems subjected to dynamic loading conditions [1,2,3,4,5,6,7,8,9]. In particular, the interaction of stress waves at the interfaces of materials with differing impedances is critical for understanding the behavior of waves under impact loading. The traditional models often consider idealized conditions, neglecting factors such as lateral inertia and the Poisson effect, which can significantly influence the accuracy of the predictions in real-world situations [10,11]. The Rayleigh–Love rod theory [12,13,14], an advanced theoretical framework, incorporates these aspects to provide a more comprehensive understanding of one-dimensional stress wave propagation behavior.
This paper delves into the complex dynamics of stress wave propagation in a system composed of a moving striker rod and an initially stationary semi-infinite Rayleigh–Love rod. The focus is primarily on the role of different general impedances in modifying the wave’s response during its propagation through the system. The general impedance Z = Aρc, where A is the cross-sectional area, ρ is the mass density, and c is the wave velocity of a rod, is a parameter that embodies the mechanical resistance against wave transmission at an interface, and it is a crucial factor in determining the wave’s behavior upon encountering a discontinuity of material properties. This research utilizes the Rayleigh–Love rod theory, which accounts for lateral inertia and the Poisson effect, to provide a more detailed insight into stress wave propagation in rods.
The novelty of this paper lies in addressing the challenges associated with the inversion of Laplace transforms for analyzing stress wave propagation, taking into account lateral inertia and Poisson’s effects. This represents a significant advancement over previous research [13,14], which only discussed impacts between a striker rod and a semi-infinite rod with the same impedance. To solve complex wave equations derived from the Rayleigh–Love rod theory with varying general impedances, we employ the approximate formula for numerical inversion of Laplace transforms introduced by Valsa et al. [15]. The precision achieved by this method can be very high, even in cases in which the transformations are irrational or transcendental functions of the complex variable. The algorithms can be easily and efficiently implemented using programming packages such as MATLAB or Mathematica. Based on the results obtained from the proposed method, this study demonstrates that alterations in general impedances at the interface between the striker rod and the semi-infinite rod lead to distinct reflection and transmission behaviors of the stress waves. These interactions can result in the amplification or attenuation of the wave, underscoring the significant impact that impedance variations can have on the propagation of stress waves.
Through a thorough analysis and comparison with the results of previous research, this paper confirms the accuracy of the findings and contributes to the body of knowledge by providing detailed insights into the effects of different general impedances on stress wave propagation in rods. The implications of this research are vast, offering valuable guidance for the design and evaluation of mechanical systems and materials subjected to dynamic impacts, ensuring their integrity and performance in practical applications.

2. Materials and Methods

The research model in this paper is presented in Figure 1. The governing equation of the Rayleigh–Love theory [12] for wave propagation in rods is as follows:
2 u 1 t 2 = c 01 2 2 u 1 x 2 + υ 1 2 κ 1 2 4 u 1 x 2 t 2 , L x 0 ,
2 u 2 t 2 = c 02 2 2 u 2 x 2 + υ 2 2 κ 2 2 4 u 2 x 2 t 2 , 0 x < ,
where c 01 = E 1 / ρ 1 and c 02 = E 2 / ρ 2 are the constant wave velocities without considering Poisson’s effect in the striker rod and the the semi-infinite rod, respectively. E 1 , E 2 are the Young’s moduli values of the striker rod and the semi-infinite rod, respectively. ρ1 and ρ2 are the mass densities; υ1 and υ2 are the Poisson ratios; and κ 1 , κ 2 are the radii of gyration of the cross-section of the striker rod and the semi-infinite rod, respectively. The impedances are Z 1 = A 1 ρ 1 c 1 and Z 2 = A 2 ρ 2 c 2 for the striker rod and the semi-infinite rod, respectively. c1 and c2 are the constant wave velocities considering Poisson’s effect [16] in the striker rod and the semi-infinite rod, respectively, and are determined as follows:
c 1 = 1 υ 1 E 1 1 + υ 1 1 2 υ 1 ρ 1 ,
c 2 = 1 υ 2 E 2 1 + υ 2 1 2 υ 2 ρ 2 .
The initial conditions in terms of displacement u can be written in a piecewise manner:
u 1 ( x , 0 ) = 0 , u 1 ( x , 0 ) t = 2 v 0 , L x 0 ,
u 2 ( x , 0 ) = 0 , u 2 ( x , 0 ) t = 0 , 0 x < ,
The boundary conditions and continuity conditions at x = L , x = 0 , and x = are given by:
E 1 A 1 u 1 x ( L , t ) = 0 , x = L ,
u 1 ( 0 , t ) = u 2 ( 0 + , t ) , x = 0 ,
E 1 A 1 u 1 ( 0 , t ) x = E 2 A 2 u 2 ( 0 + , t ) x , x = 0 ,
lim x u 2 ( x , t ) = 0 , x = , t 0 .
To make the Laplace transform process more convenient, we can convert the quantities to dimensionless quantities for substitution into Equations (1) and (2). This is achieved by introducing the following dimensionless variables and parameters:
u ¯ 1 = u 1 D 2 ,   u ¯ 2 = u 2 D 2 , x ¯ = x D 2 ,   t ¯ = t ¯ 2 = c 2 t D 2 ,   t ¯ 1 = c 1 t D 2 = c 1 c 2 t ¯ , L ¯ = L D 2 ,   κ ¯ 1 = κ 1 D 2 ,   κ ¯ 2 = κ 2 D 2 ,   b 1 = υ 1 κ ¯ 1 ,   b 2 = υ 2 κ ¯ 2 ,   v ¯ 0 = v 0 c 1 ,
where D2 is the diameter of the semi-infinite rod.
Substituting Equation (11) into Equations (1) and (2), we obtain:
2 u ¯ 1 t ¯ 2 = m 2 2 u ¯ 1 x ¯ 2 + b 1 2 4 u ¯ 1 x ¯ 2 t ¯ 2 , L ¯ x ¯ 0
2 u ¯ 2 t ¯ 2 = 2 u ¯ 2 x ¯ 2 + b 2 2 4 u ¯ 2 x ¯ 2 t ¯ 2 , 0 x ¯ <
where m = c 1 c 2 .
Apply Laplace transformation to find the displacement U 1 ( x ¯ , s ) and U 2 ( x ¯ , s ) in s domain (see in Appendix A) as follows:
U 1 ( x ¯ , s ) = 2 v ¯ 0 E 2 A 2 m 2 + b 1 2 s 2 e x ¯ s m 2 + b 1 2 s 2 + e ( x ¯ + 2 L ¯ ) s m 2 + b 1 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 e 2 L ¯ s m 2 + b 1 2 s 2 + E 2 A 2 m 2 + b 1 2 s 2 + E 1 A 1 1 + b 2 2 s 2 + 2 v ¯ 0 s 2 , L ¯ x ¯ 0
U 2 ( x ¯ , s ) = 2 v ¯ 0 1 e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 e 2 L ¯ s m 2 + b 1 2 s 2 + E 2 A 2 m 2 + b 1 2 s 2 + E 1 A 1 1 + b 2 2 s 2 e x ¯ s 1 + b 2 2 s 2 , 0 x ¯ <
When E1 = E2, A1 = A2, ρ1 = ρ2, υ1 = υ2, c1 = c2, m = 1, and b1 = b2 = b, Equations (14) and (15) are respectively reduced to:
U 1 ( x ¯ , s ) = v ¯ 0 s 2 2 e x ¯ s 1 + b 2 s 2 e ( x ¯ + 2 L ¯ ) s 1 + b 2 s 2 , L ¯ x ¯ 0
U 2 ( x ¯ , s ) = v ¯ 0 s 2 e x ¯ s 1 + b 2 s 2 e ( x + 2 L ¯ ) s 1 + b 2 s 2 , 0 x ¯ <
Equations (16) and (17) are the transformed functions of a semi-infinite rod under the impact of a striker rod of the same material and cross-sectional area, as provided by previous research [14]. This verifies the correctness of Equations (14) and (15).

3. Inverse Laplace Transform

Inverse Laplace transformation of Equations (14) and (15) is challenging, so a numerical method was utilized to perform this task. By using the inverse Laplace transform method provided by Valsa et al. [15], we can compute the inverse Laplace transforms directly and effectively. The derivation of the inverse Laplace transform by Valsa et al. was carried out from the basic transform, as follows:
f ( t ) = L 1 { F ( s ) } = 1 2 π i γ i γ + i F ( s ) e s t d s ,
where the integration variable is s = γ + .
The transform F(s) is supposed to fulfil the following assumptions:
  • F(s) is regular for Re{s} > 0.
  • lim | s | F ( s ) = 0 .
  • F*(s) = F(s*), where the asterisk denotes the complex conjugate value.
To simplify the evaluation of the integral, we approximate the exponential function as follows:
e s t E c ( s t , a ) = e a 2 cosh ( a s t ) = e a e a e s t + e a e s t = e s t 1 + e 2 a e 2 s t ,
where E c ( s t , a ) is an approximate expression of e st using the hyperbolic cosine function cosh(x), with x = a − st, and a is a real constant.
When a > γt, then e 2 a e 2 s t < < 1 , and the fraction 1 1 + e 2 a e 2 s t can be expanded as a convergent Maclaurin series as follows:
E c ( s t , a ) = e s t + n = 1 ( 1 ) n e 2 n a e ( 2 n + 1 ) s t .
The sum in the expression represents the error of approximation of the exponential function that can be arbitrarily suppressed by a choice of the parameter a. Then, Equation (18) becomes:
f ( t ) f c ( t , a ) = 1 2 π i γ i γ + i F ( s ) E c ( s t , a ) d s = 1 2 π i γ i γ + i F ( s ) e s t d s + 1 2 π i γ i γ + i F ( s ) n = 1 ( 1 ) n e 2 n a e ( 2 n + 1 ) s t d s = f ( t ) + n = 1 ( 1 ) n e 2 n a f ( 2 n + 1 ) t ) = f ( t ) e 2 a f ( 3 t ) + e 4 a f ( 5 t ) = f ( t ) + ε c ( t , a ) ,
where f c ( t , a ) is the approximate function of f(t).
Moreover, the error ε c ( t , a ) in Equation (21) of the inverse transform can be kept at low values by appropriate choice of the parameter a. The function E c ( s t , a ) can be expressed as an infinite sum of rational functions in st, since (see [17]):
1 cosh z = 2 π n = 0 ( 1 ) n ( n + 1 2 ) ( n + 1 2 ) 2 π 2 + z 2 .
Combining Equation (19), Equation (21), and Equation (22), we have:
E c ( s t , a ) = π e a n = 0 ( 1 ) n ( n + 1 2 ) ( n + 1 2 ) 2 π 2 + ( a s t ) 2 .
Substituting Equation (23) into Equation (21), we obtain:
f c ( t , a ) = e a 2 i γ i γ + i F ( s ) n = 0 ( 1 ) n ( n + 1 2 ) ( n + 1 2 ) 2 π 2 + ( a s t ) 2 d s .
By interchanging the sequence of summation and integration, we obtain:
f c ( t , a ) = e a 2 i n = 0 ( 1 ) n ( n + 1 2 ) I n ,
where   I n = γ i γ + i F ( s ) ( n + 1 2 ) 2 π 2 + ( a s t ) 2 d s = γ i γ + i F ( s ) G n ( s ) d s = γ i γ + i H n ( s ) d s .
To evaluate the integral, we completed the integration path using an arc contour with an infinite radius. Due to Assumption 2 of Laplace transform, the function F(s) equals zero along this semicircle. The integral In is then equal to the sum of the residuals of Hn(s) in the poles lying inside the integration contour (see Figure 2), multiplied by 2πi. The poles of Hn(s), in turn, are the poles of F(s) (these are concentrated in R e { s } 0 plus the roots of the equation Gn(s) = 0, as follows:
s 1 , 2 = a ± i ( n + 1 2 ) π t ,
which lie in the right half-plane, to the right of the original path of integration Re{s} = γ. It is therefore advantageous to close the path of integration here and to calculate I n as follows:
I n = 2 π i res H n ( s ) | s = s 1 + res H n ( s ) | s = s 2 = 2 π i F ( s 1 ) G n ( s 1 ) + F ( s 2 ) G n ( s 2 ) ,
where G ( s ) = d G ( s ) d s .
The negative sign results from the negative sense of integration. Obviously, G n s = 2 t ( a s t ) ; therefore,
I n = 2 π i F ( s 1 ) 2 π i ( n + 1 2 ) F ( s 2 ) 2 π i ( n + 1 2 ) = F ( s 1 * ) F ( s 1 ) t ( n + 1 2 ) = F * ( s 1 ) F ( s 1 ) t ( n + 1 2 ) = 2 l Im F ( s 1 ) t ( n + 1 2 ) = 2 i Im F a + i ( n + 1 2 ) π / t ) t ( n + 1 2 ) .
Finally, we obtain approximate formulae for the numerical inversion of Laplace transformation, as follows:
f c ( t , a ) = e a t n = 0 1 n Im F a t + i ( n + 1 2 ) π t .
In the programming process, as the value of n cannot start from 0 but must begin with 1, Equation (30) is adjusted as follows:
f c ( t , a ) = e a t n = 1 n e n d ( 1 ) n Im F a t + i ( n 1 ) π t ,
where f c ( t , a ) denotes the calculated inverse Laplace transformation of the function F(s), dependent on time t and the convergence parameter a. F a / t + i n 1 π / t is Laplace transform function F(s), evaluated at specific points in the complex plane. The argument of F is a complex number, with a / t being the real part and i n 1 π / t being the imaginary part; the term i is the imaginary unit; nend is the final step of the loop.
The process of numerically inverting Laplace transforms is essential in dissecting complex systems, and MATLAB offers a key function, denoted as f ( t ) = INVLAP F ( s ) [18]. To accomplish this task, we utilized the inverse Laplace transform technique offered by Valsa et al. [15], which enabled us to carry out the calculation of inverse Laplace transforms directly and efficiently. The flexibility of the F(s) function within INVLAP is notable; it adeptly handles a diverse array of complex variable functions where s is raised to an exponent a. This exponent can be any real value, including both integers and fractions, providing INVLAP with the ability to solve fractional problems and invert functions that consist of both rational and irrational components. One of the primary benefits of utilizing INVLAP is its ability to bypass the often-complex step of calculating the poles or zeros of F(s). The function instead evaluates F(s) at chosen complex values of s. Based on Equation (31), we wrote a code for the inverse Laplace transformation of Equations (14) and (15) to find the displacement functions u ¯ 1 x ¯ , t ¯ and u ¯ 2 x ¯ , t ¯ .
Then, the velocities in the striker rod and the semi-infinite rod are determined:
v ¯ 1 ( x ¯ , t ¯ ) = u ¯ 1 t ¯ , L ¯ x ¯ 0 ,
v ¯ 2 ( x ¯ , t ¯ ) = u ¯ 2 t ¯ , 0 x ¯ < .
The stress wave propagation in the striker rod and the semi-infinite rod is determined as follows:
σ ¯ 1 ( x ¯ , t ¯ ) = u ¯ 1 x ¯ + b 1 2 m 2 3 u ¯ 1 x ¯ t ¯ 2 , L ¯ x ¯ 0 ,
σ ¯ 2 ( x ¯ , t ¯ ) = u ¯ 2 x ¯ + b 2 2 3 u ¯ 2 x ¯ t ¯ 2 , 0 x ¯ < .
The algorithm schematic for the numerical inversion Laplace transform process is illustrated in Figure 3. Applying the numerical inversion Laplace transform to Equations (14) and (15), we obtain the displacement functions u ¯ 1 x ¯ , t ¯ and u ¯ 2 x ¯ , t ¯ , respectively. From the obtained results, Equations (32) and (33) are used to calculate the corresponding velocities v ¯ 1 x ¯ , t ¯ and v ¯ 2 x ¯ , t ¯ , respectively. Furthermore, stress functions σ ¯ 1 x ¯ , t ¯ and σ ¯ 2 x ¯ , t ¯ are derived by applying Equations (34) and (35), respectively. Finally, these non-dimensional quantities are converted to dimensional forms using Equation (11), resulting in the physical measures of displacements u 1 x , t and u 2 x , t , velocities v 1 x , t and v 2 x , t , and stresses σ 1 x , t and σ 2 x , t for the striker and the semi-infinite rod, respectively. The obtained results should be compared to previous research to confirm the accuracy of the proposed method.

4. Validation of the Proposed Numerical Algorithm for the Inversion of Laplace Transformation

4.1. Exampe 1

To confirm the accuracy of the proposed method and the computational analysis code, consider the initial value problem as follows
f ( t ) + 9 f ( t ) = 2 sin 3 t ,
given the initial conditions are:
f ( 0 ) = 1 , f ( 0 ) = 0 .
Using Laplace transform to find f(t), take Laplace transform of both sides of Equation (36):
L f ( t ) + L 9 f ( t ) = L 2 sin 3 t ,
s 2 F ( s ) s f ( 0 ) f ( 0 ) + 9 F ( s ) = 6 s 2 + 9 .
Substituting Equation (37) into Equation (39), we obtain F(s) as follows:
F ( s ) = 6 s 2 + 9 2 + s s 2 + 9 .
Take the inverse Laplace transform Equation (40) to find f(t) as follows:
f ( t ) = L 1 F ( s ) = L 1 6 s 2 + 9 2 + L 1 s s 2 + 9 .
The inverse Laplace transform of the second term is:
L 1 s s 2 + 9 = cos 3 t .
Since the inverse Laplace transform of the first term is more complicated, we can use the convolution theorem [19] to find the Laplace transform of the first term:
L 1 6 s 2 + 9 2 = L 1 2 3 3 s 2 + 9 3 s 2 + 9 = 2 3 f g ( t ) ,
where f(t) = g(t) = sin3t, the convolution of Equation (43) is determined as:
L 1 6 s 2 + 9 2 = 2 3 0 t sin 3 τ sin 3 ( t τ ) d τ .
Evaluating the integral of Equation (44), we obtain:
2 3 0 t sin 3 τ sin 3 ( t τ ) d τ = 1 9 sin 3 t 1 3 t cos 3 t .
Thus, the exact solution of Equation (36) is:
f ( t ) = 1 9 sin 3 t 1 3 t cos 3 t + cos 3 t .
Instead of using the convolution theorem [19], we can apply the numerical inversion of the Laplace transform for F(s) as presented in Section 3. The results between the two methods are shown in Figure 4, which demonstrates a high degree of agreement between the exact and numerical solutions, thereby validating the effectiveness and accuracy of the numerical approach for the numerical inversion of the Laplace transform.

4.2. Exampe 2

Example 2 demonstrates the application of the proposed method for solving the Rayleigh–Love wave equation, as modeled in Figure 5. A comparison was performed between the results obtained from approximate formulae for numerical inversion of Laplace transforms and Yang’s solution as presented in his research paper [14]. The parameters used in the model are shown in Table 1. The wave velocity (c) considering Poisson’s effect [16] in the striker rod and the semi-infinite rod is determined as follows
c = 1 υ E 1 + υ 1 2 υ ρ ,
A detailed comparison of the stress propagation over time at varying positions, with the results obtained from both methods, is illustrated in Figure 6, Figure 7 and Figure 8. These figures illustrate a comparison between the results from the proposed numerical solution and Yang’s analytical solution [14]. These results validate the reliability and accuracy of the numerical approach for analyzing stress wave propagation in the semi-infinite Rayleigh–Love rod. Figure 6, Figure 7 and Figure 8 illustrate the dynamic behavior of displacement, velocity, and stress at different positions in the striker rod and the semi-infinite Rayleigh–Love rod due to the impact of an initial velocity 2v0. This demonstrates that the results from the proposed numerical method adequately match those of Yang’s analytical solution. This agreement highlights the validity of the results obtained from the numerical analysis algorithm and confirms its correctness in predicting the behavior of stress waves in Rayleigh–Love rods. As shown in Figure 8, the number of oscillations of the stress decreases as x increases, due to the effects of lateral inertia and Poisson’s ratio. This results in the dispersion of the signal during the propagation in the x-axis because of the effect of the term υ 2 κ 2 4 u x 2 t 2 in the governing equation of Rayleigh–Love rod theory.

5. Investigation of Wave Propagation in a Rayleigh–Love Rod under the Collinear Impact of a Striker Rod with Different Impedances

This section verifies stress wave propagation in a rod under the collinear impact of a striker rod with different impedances. The proposed numerical solution is compared to the results provided by Shin et al. [11], which do not include Poisson’s effect, the stress of the striker rod in the kth cycle is:
σ 1 , k = ρ 2 A 2 c 2 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 ρ 1 c 1 v p , k = ρ 2 A 2 c 2 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 ρ 1 A 1 c 1 ρ 2 A 2 c 2 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 k 1 ρ 1 c 1 v i , k = 1 , 2 , 3 , 4 ,
The stress of the semi-infinite rod in the kth cycle is:
σ 2 , k = ρ 1 A 1 c 1 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 ρ 2 c 2 v p , k = ρ 1 A 1 c 1 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 ρ 1 A 1 c 1 ρ 2 A 2 c 2 ρ 1 A 1 c 1 + ρ 2 A 2 c 2 k 1 ρ 2 c 2 v i , k = 1 , 2 , 3 , 4 .
where v 1 , k is the particle velocity of the striker rod in the cycle kth and v 2 , k is the particle velocity of the semi-infinite rod in the cycle kth.

5.1. Investigation of Stress Wave Propagation in a Semi-Infinite Rod and a Striker Rod at the Interface x = 0 m

Based on the stress equation with cyclically varying values without considering the effects of lateral inertia and Poisson’s ratio on the striker rod (Equation (48)) and the semi-infinite rod (Equation (49)), we will now compare the results with those of the proposed numerical solution. This comparison aims to evaluate the results obtained from modifications in the cross-sectional area and material properties of the rods, as indicated in Table 2. In the first case, we investigate the position at x = 0, where the interaction between the two rods occurs upon impact. The results obtained are illustrated in Figure 9 and Table 3.
In Table 3 and Figure 9, the results of comparing stress wave propagation when varying the cross-sectional area of the striker rod for three cases (Z1/Z2 = 0.5, Z1/Z2 = 1, and Z1/Z2 = 2) are presented. The results obtained from the proposed numerical solution and the solutions with no Poisson’s effect reported in Shin and Kim’s paper [11] are in good agreement. When considering the influence of Poisson’s effect, the stress wave shapes oscillate around the without Poisson’s effect line until reaching a steady state, at which point the two results are the same. When comparing the three cases of varying cross-sectional areas, a common observation is that compressive stress is generated in the first cycle for all cases, with a significant variation observed in the second cycle. For Z1/Z2 = 1, the stresses in both the striker and the semi-infinite rod are equal and amount to −263.286 MPa. The stress returns to zero in the second cycle and continues to oscillate around this value. In the case of Z1/Z2 = 0.5, the stress in the striker is amplified to −351.05 MPa, while the stress in the semi-infinite rod is attenuated to −175.52 MPa. Tensile stress appears in the second cycle, followed by compressive stress, alternating in this manner and gradually approaching zero. This indicates that the two rods will separate after impact when Z1 < Z2. However, the constraint condition expressed by Equation (8) allows the existence of tensile stress at the interface x = 0. When Z1/Z2 = 2, the stress in the striker rod is attenuated to −175.52 MPa, and the stress in the semi-infinite rod is amplified to −351.05 MPa. The stress in each cycle is compressive, decreasing gradually and approaching zero. The results of this research are similar to those reported in Shin and Kim’s paper [11]. Notably, the duration of each cycle in all three cases is identical, equal to 2L/c1 = 69.17 μs. The oscillation observed in wave propagation is not due to the numerical algorithm effect but is the actual effect of lateral inertia and Poisson’s ratio.

5.2. Investigation of Stress Wave Propagation in a Striker Rod at x = −0.1 m

This case examines the propagation of stress waves at the position x = −0.1 m on the striker rod for three scenarios related to changes in cross-sectional area, Young’s modulus, and density, with ratios of 0.5, 1, and 2, respectively, as shown in Figure 10 and Table 4. A remarkable aspect of this case is that a half-cycle is generated with either compressive or tensile stress, and there is no stress during the half-cycle when the rod is released. After the impact in the first half-cycle, compressive stress is generated, and then, depending on the contact surface characteristics between the two rods (cross-sectional area and material properties), the subsequent cycle will produce either compressive or tensile stress. This indicates that the duration of the stress pulse Ti in the striker rod is linearly decaying in each cycle as T j = 2 L c 1 1 + x L , where L x 0 . The wave velocity (c1) is calculated as Equation (3).

5.3. Investigation of Stress Wave Propagation in a Semi-Infinite Rod at x = 0.5 m

This case investigates stress wave propagation at the position x = 0.5 m on the semi-infinite rod with different impedances, as shown in Figure 11 and Table 5. All cases initially generate compressive stress, with notable variations in the second cycle. Specifically, the Z1 = 0.5Z2 case alternates between tensile stress and compressive stress, approaching zero; the Z1 = 0.5Z2 case sees stress oscillating around zero; and Z1 = 2Z2 consistently shows decreasing compressive stress towards zero. Particularly, when Z1 = Z2, the stress wave appears in only one stress compression cycle, after which the stress returns to oscillate around zero. Each cycle’s duration remains constant at T i = 2 L / c 1 = 69.17 μ s in all cases.
To further validate the accuracy of the method we propose, the results obtained are compared to finite element method (FEM) results. This 3D finite element analysis is a kind of numerical experiment to simulate the real performance of the system. The 3D FEM model is constructed as shown in Figure 12. The parameters for material properties and the geometries of the striker rod and the finite rod used in this model are detailed in Table 1; the length of the striker is 0.2 m, and that of the finite rod is 2.5 m. The approximate element size is 0.01 m, with a maximum deviation factor of 0.05 and a minimum element size of 0.001 m. The boundary conditions for the finite rod are set with one end free at x = 0 and fixed at x = L. In the Abaqus finite element analysis code, the solid element type C3D8R was employed and had 11,193 nodes and 8672 elements in the model, as shown in Figure 12b. In the simulation conducted using the Abaqus/Explicit [20], the system of equations was solved through an explicit time integration method. The time increment selected for the simulation is Δt = 0.5 × 10−7 (s). This setting guarantees that the condition ΔtLmin/c = 1.729 × 10−7 (s), as recommended for stable explicit dynamic analysis in Abaqus [20], where Lmin is the smallest element dimension in the mesh and c is the wave speed in the material of the rod, ensuring that the simulation remains stable and accurate. Since the rods are axisymmetric, resulting in no friction force in the interaction of this model, the contact interaction between bodies is defined as frictionless in the Abaqus contact model. Figure 13 demonstrates the accuracy of the proposed numerical solution in simulating displacement, velocity, and stress in a rod, as shown by its comparison to FEM results. Particularly, in axial displacement, the proposed solution aligns closely with the FEM, indicating high accuracy in predicting displacement behavior over time. Despite minor discrepancies in velocity oscillations, the overall trend closely matches that of the FEM, underscoring the method’s reliability. The stress response also shows a strong correlation with the FEM standard, solidifying the proposed method’s potential as a valuable tool for dynamic stress analysis. In addition to its comparison with the FEM, the proposed numerical solution demonstrates considerable consistency with Shin and Kim’s solution [11], which—for Poisson’s ratio—is equal to zero. This agreement also confirms the accuracy of the proposed numerical solution, reinforcing its potential as a reliable method for understanding the impact dynamics of rods with varying cross-sectional areas.
In Figure 14, when considering the cases with the same impedance but varying cross-sectional areas A1 and A2, as well as differing ρ1c1 and ρ2c2, we can observe that the generated stress wave undergoes only one cycle, yet the duration of the stress signal varies. It is evident that as the modulus of elasticity E1 increases, the duration of the stress signal decreases. When both rods have the same impedance, the generated incident waveforms are simpler compared to those produced by rods with differing impedances, as discussed in Section 5.1. Thus, for practical purposes, to simplify the wave propagation analysis and to obtain a clearer output signal, it is advantageous to use rods with identical impedance. The results obtained in this study are compared with those of Shin et al. [11] to evaluate the accuracy of the proposed method. Shin et al. [11] present the impact of a striker on a bar with different general impedance, leading to various wave behaviors that can cause the striker and bar to separate or bond together. Shin et al. [11] do not consider the effects of lateral inertia and Poisson’s ratio, resulting in the transmission of rectangular pulses without oscillation in wave propagation as shown in Figure 14. Our study is more comprehensive; it includes the oscillation in wave propagation. This oscillation leads to higher stress magnitudes, which can potentially cause different structural performances. Therefore, it is crucial to consider the effects of lateral inertia and Poisson’s ratio during wave propagation. We verify our results through comparisons with the results in Shin et al. [11] and FEM results. The result of Shin et al. [11] is just a special case of our general analysis model.

6. Travel Time History

In Figure 15, Figure 16 and Figure 17, points P1 to P5L are positions located in the striker rod and P5R to P10 are positions located in the semi-infinite rod. These illustrate the stress wave travel time histories at positions P1 to P10 of two rods with various cross-sectional areas (A1 = A2, A1 = 2A2, and A1 = 0.5A2) at the initial impact velocity 2v0. These cases facilitate a fine understanding of the interaction between stress waves with varying cross-sectional areas.
From the results, it is clear that the amplitude of the stress waves differs in each case, with the wave amplitude increasing as the cross-sectional area A1 increases. This demonstrates the sensitivity of wave transmission to the cross-sectional area; specifically, when A1 is larger than A2, the amplitude of the stress wave across the rods increases, implying that a larger cross-sectional area facilitates the transmission of more impact energy. In contrast, when A1 is smaller than A2, the amplitude decreases, which may indicate increased reflection rather than transmission of wave energy at the interface between the rods. Positions P1 to P5L in the striker rod have pre-existing wave energy at the interface between the rods. In the initial state, the striker rod moves with a velocity of 2v0, creating a stress in the semi-infinite rod, which will cause a reflection back to the striker rod. The waveform shape at points P5R to P10 is the same in form, but the timing of the wave initiation differs among the positions, and the amplitude of wave propagation is different in all cases, reflecting the impact of the varying cross-sectional areas on the amplitude of wave propagation at these positions. This phenomenon underscores the critical role of the cross-sectional area in the modulation of wave energy transmission, where a larger A1 facilitates the conveyance of a greater amount of wave energy across the rods.
Due to the traction-free boundary condition at the left end, the striker stress is zero at any instance at this point. The stress at each point reaches its maximum value when a stress wave travels through and then gradually decreases to zero. The duration of the stress pulse Tj in the striker rod decays linearly in each cycle as T j = 2 L c 1 1 + x L , where L x 0 , and the duration of the stress pulse Ti in the semi-infinite rod is T i = 2 L c 1 , where x 0 . The wave velocity c1 in all of the cases presented here is the same, because only the cross-sectional area of the striker rod is changed in the three cases, while all of the material properties are the same. In the case where Z1 = Z2, as shown in Figure 15, the stress wave propagation occurs in a single cycle, resulting in a clear signal after cycle 1. In the case where Z1Z2, as shown in Figure 16 and Figure 17, the wave propagation occurs over multiple consecutive cycles, creating a complex wave signal. This indicates that in order to easily identify the stress wave signal, using the condition Z1 = Z2 at the material interface is the most appropriate. Specifically, Figure 16 shows that when Z1 > Z2, a series of compression waves are displayed over consecutive cycles, causing the two rods to remain joined. Conversely, Figure 17 demonstrates that when Z1 < Z2, a tensile wave appears after the first cycle, leading to the separation of the two rods after the impact.

7. Conclusions

This study presents an in-depth analysis of stress wave propagation in a semi-infinite Rayleigh–Love rod system when impacted by a striker rod with varying impedances. The different impedances between the striker and the semi-infinite rod give rise to distinct wave signals, occurring in one cycle or multiple impact cycles. These observations offer significant insights into the effect of impedance contrast on the characteristics of the waves. The influence of lateral inertia and Poisson’s ratio on longitudinal waves in rods leads to higher stress results compared to scenarios where these factors are not considered.
This research applied approximate formulae for numerically inverting Laplace transformations to handle complex equations. This is an effective method for solving the problem of wave propagation in rods with different impedances; it provides a reliable and accurate tool to analyze stress wave behaviors in rods involving interactions with different impedances. The comparison analyses with Yang’s analytical solution with the same impedance [14], Shin and Kim’s solution [11] without effects of lateral inertia and Poisson’s ratio, and numerical experiment by finite element analysis affirmed the high accuracy of this numerical inversion Laplace transform approach.
This study provides examples showing that changing the material parameters of the striker rod results in different waveforms. These findings enable the development of customized experimental models that can precisely achieve goals in materials testing and analysis.
An interesting finding from this study is that when Z1 < Z2, a positive reflection wave occurs after the first cycle, causing the two rods to separate; thus, there is no need to continue the analysis after Cycle 1. Conversely, when Z1 > Z2, multiple cycles of compressive waves occur, leading to the rods adhering together. Notably, with the same impedance, the clarity of the signal is optimal, making it easier to analyze and identify the geometric and material properties of the sample.

Author Contributions

C.-Y.W.: advisor, conceptualization, methodology, writing—review, and editing. N.N.T.: conceptualization, methodology, formal analysis, numerical simulation, validation, data curation, visualization, and writing—original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Solution of Laplace Transform of Equations (11) and (12)

This section derives the displacement U 1 ( x ¯ , s ) and U 2 ( x ¯ , s ) in the s domain.
We use Laplace transformation properties for the derivative as follows
L d n f d t n = s n F ( s ) s n 1 f ( 0 ) s n 2 f ( 0 ) s f n 2 ( 0 ) f n 1 ( 0 ) .
We can convert the initial and boundary condition terms in Laplace transformation. First, the initial conditions in terms of displacement and velocity can be written as follows:
u ¯ 1 ( x ¯ , 0 ) = 0 , u ¯ 1 ( x ¯ , 0 ) t ¯ = 2 v ¯ 0 , L ¯ x ¯ < 0 ,
u ¯ 2 ( x ¯ , 0 ) = 0 , u ¯ 2 ( x ¯ , 0 ) t ¯ = 0 , 0 x ¯ < .
The boundary conditions at x ¯ = L ¯ , x ¯ = 0 , and x ¯ = are given by:
u ¯ 1 x ¯ ( L ¯ , t ¯ ) = 0 ,
u ¯ 1 ( 0 , t ¯ ) = u ¯ 2 ( 0 + , t ¯ ) , u ¯ 1 ( 0 , t ¯ ) x ¯ = E 2 A 2 E 1 A 1 u ¯ 2 ( 0 + , t ¯ ) x ¯ , 0 t ¯ < ,
lim x ¯ u ¯ 2 ( x ¯ , t ¯ ) = 0 , t ¯ 0 .
If we apply Laplace transformation to both sides of Equations (7) and (8), we have:
L 2 u ¯ 1 t ¯ 2 = m 2 L 2 u ¯ 1 x ¯ 2 + b 1 2 L 4 u ¯ 1 x ¯ 2 t ¯ 2 , L ¯ x ¯ 0 ,
L 2 u ¯ 2 t ¯ 2 = L 2 u ¯ 2 x ¯ 2 + b 2 2 L 4 u ¯ 2 x ¯ 2 t ¯ 2 , 0 x ¯ < ,
where L is the operator of Laplace transformation.
L u ¯ 1 ( x ¯ , t ¯ ) = U 1 ( x ¯ , s ) ,
L u ¯ 2 ( x ¯ , t ¯ ) = U 2 ( x ¯ , s ) .
From Equations (A2) and (A3), we have:
L 2 u ¯ 1 ( x ¯ , t ¯ ) t ¯ 2 = s 2 U 1 ( x ¯ , s ) s u ¯ 1 ( x ¯ , 0 ) u ¯ 1 ( x ¯ , 0 ) t ¯ = s 2 U 1 ( x ¯ , s ) 2 v ¯ 0 ,
L 2 u ¯ 1 x ¯ 2 = 2 U 1 ( x ¯ , s ) x ¯ 2 ,
L 4 u ¯ 1 x ¯ 2 t ¯ 2 = L 2 2 u ¯ 1 / t ¯ 2 x ¯ 2 = L 2 s 2 U 1 ( x ¯ , s ) s u ¯ 1 ( x ¯ , 0 ) u ¯ 1 ( x ¯ , 0 ) t ¯ x ¯ 2 = s 2 2 U 1 ( x ¯ , s ) x ¯ 2 .
Then, Equation (A7) can be rewritten as follows:
d 2 U 1 d x ¯ 2 s 2 m 2 + b 1 2 s 2 U 1 = 2 v ¯ 0 m 2 + b 1 2 s 2 , L ¯ x ¯ 0 .
Similarly, Equation (A8) can be rewritten as follows:
d 2 U 2 d x ¯ 2 s 2 1 + b 2 2 s 2 U 2 = 0 , 0 x ¯ < .
We now present step-by-step methods to solve Equations (A14) and (A15). First, find the homogeneous solution of Equation (A14). We have the following homogeneous differential equation:
d 2 U 1 d x ¯ 2 s 2 m 2 + b 1 2 s 2 U 1 = 0 .
Assuming U 1 = e r x ¯ , we have:
U 1 = r e r x ¯ , U 1 = r 2 e r x ¯ .
Then, substituting Equation (A17) into Equation (A16), we obtain:
r 2 e r x ¯ s 2 m 2 + b 1 2 s 2 e r x ¯ = r 2 s 2 m 2 + b 1 2 s 2 e r x ¯ = 0 .
Solving the characteristic equation of Equation (A18), we obtain two roots r1 and r2 as follows:
r 1 = s m 2 + b 1 2 s 2 , r 2 = s m 2 + b 1 2 s 2 .
Therefore, the homogeneous solution is:
U 1 h ( x ¯ , s ) = B 1 e r 1 x ¯ + B 2 e r 2 x ¯ = B 1 e x ¯ s m 2 + b 1 2 s 2 + B 2 e x ¯ s m 2 + b 1 2 s 2 ,
where B1 and B2 are constants.
Now, we need to find a particular solution for Equation (A14). Let us assume that the particular solution is of the following form:
U 1 p ( x ¯ , s ) = A ( s ) .
where A(s) is a function of s that does not depend on x ¯ . Using the method of undetermined coefficients, we can find:
A s = 2 v ¯ 0 s 2 .
Thus, we can obtain the particular solution as follows:
U 1 p ( x ¯ , s ) = 2 v ¯ 0 s 2 .
The general solution for Equation (A14) is as follows:
U 1 ( x ¯ , s ) = U 1 h ( x ¯ , s ) + U 1 p ( x ¯ , s ) = B 1 e x ¯ s m 2 + b 1 2 s 2 + B 2 e x ¯ s m 2 + b 1 2 s 2 + 2 v ¯ 0 s 2 , L ¯ x ¯ 0 .
For Equation (A15), d 2 U 2 d x ¯ 2 s 2 1 + b 2 2 s 2 U 2 = 0 , 0 x ¯ < .
Similarly, we can find a solution to this homogeneous differential Equation (A15) as follows:
U 2 ( x ¯ , s ) = B 3 e x ¯ s 1 + b 2 2 s 2 + B 4 e x ¯ s 1 + b 2 2 s 2 , 0 x ¯ < ,
where B1, B2, B3, and B4 in Equations (A20) and (A25) are constants to be determined by the transformed boundary conditions.
With the boundary condition U 1 ( L ¯ , s ) = 0 , from Equation (A20), we have:
U 1 ( x ¯ , s ) = B 1 e x ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 + B 2 e x ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 ,
U 1 ( L ¯ , s ) = B 1 e L ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 + B 2 e L ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 = 0 ,
B 1 = B 2 e 2 L ¯ s m 2 + b 1 2 s 2 .
With the boundary condition Equation (10), U 2 ( , s ) = 0 , from Equation (A25), we have:
U 2 ( + , s ) = B 3 e s 1 + b 2 2 s 2 + B 4 e + s 1 + b 2 2 s 2 = 0 .
So, B4 = 0 for a nontrivial solution.
Now, U 2 ( x ¯ , s ) becomes:
U 2 ( x ¯ , s ) = B 3 e s x ¯ 1 + b 2 2 s 2 .
With the boundary condition U 1 ( 0 , s ) = U 2 ( 0 + , s ) , from Equations (A24) and (A25), we have:
B 1 e 0 s m 2 + b 1 2 s 2 + B 2 e 0 s m 2 + b 1 2 s 2 + 2 v ¯ 0 s 2 = B 3 e 0 s 1 + b 2 2 s 2 + B 4 e 0 s 1 + b 2 2 s 2 ,
B 1 + B 2 + 2 v ¯ 0 s 2 = B 3 , ( Note   B 4 = 0 ) .
With the boundary condition U 1 ( 0 , t ) x ¯ = E 2 A 2 E 1 A 1 U 2 ( 0 + , t ) x ¯ , from Equations (A24) and (A25), we have:
B 1 e x ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 + B 2 e x ¯ s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 = E 2 A 2 E 1 A 1 B 3 e x ¯ s 1 + b 2 2 s 2 s 1 + b 2 2 s 2 + B 4 e x ¯ s 1 + b 2 2 s 2 s 1 + b 2 2 s 2 ,
B 1 e 0 s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 + B 2 e 0 s m 2 + b 1 2 s 2 s m 2 + b 1 2 s 2 = E 2 A 2 E 1 A 1 B 3 e 0 s 1 + b 2 2 s 2 s 1 + b 2 2 s 2 + B 4 e 0 s 1 + b 2 2 s 2 s 1 + b 2 2 s 2 ,
B 1 B 2 s m 2 + b 1 2 s 2 = E 2 A 2 E 1 A 1 B 3 s 1 + b 2 2 s 2 ,
B 1 B 2 = B 3 E 2 A 2 E 1 A 1 m 2 + b 1 2 s 2 1 + b 2 2 s 2
Solving Equations (A28), (A32) and (A36), we obtain B1, B2, and B3 as follows:
B 1 = 2 v ¯ 0 E 2 A 2 m 2 + b 1 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 ,
B 2 = 2 v ¯ 0 E 2 A 2 m 2 + b 1 2 s 2 e 2 L ¯ s m 2 + b 1 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 ,
B 3 = 2 v ¯ 0 1 e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 + e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 .
Thus, we have the transformed solutions of Equations (A24) and (A25) for the striker rod and the semi-infinite rod, as follows:
U 1 ( x ¯ , s ) = 2 v ¯ 0 E 2 A 2 m 2 + b 1 2 s 2 e x ¯ s m 2 + b 1 2 s 2 + e ( x ¯ + 2 L ¯ ) s m 2 + b 1 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 e 2 L ¯ s m 2 + b 1 2 s 2 + E 2 A 2 m 2 + b 1 2 s 2 + E 1 A 1 1 + b 2 2 s 2 + 2 v ¯ 0 s 2 , L ¯ x ¯ 0 ,
U 2 ( x ¯ , s ) = 2 v ¯ 0 1 e 2 L ¯ s m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 s 2 E 2 A 2 m 2 + b 1 2 s 2 E 1 A 1 1 + b 2 2 s 2 e 2 L ¯ s m 2 + b 1 2 s 2 + E 2 A 2 m 2 + b 1 2 s 2 + E 1 A 1 1 + b 2 2 s 2 e x ¯ s 1 + b 2 2 s 2 , 0 x ¯ <

References

  1. Tesfamariam, S.; Martín-Pérez, B. Stress wave propagation for evaluation of reinforced concrete structures. In Non-Destructive Evaluation of Reinforced Concrete Structures; Woodhead Publishing: Sawston, UK, 2010; pp. 417–440. [Google Scholar]
  2. Sawangsuriya, A. Wave propagation methods for determining stiffness of geomaterials. In Wave Processes in Classical and New Solids; IntechOpen: London, UK, 2012. [Google Scholar]
  3. Ostachowicz, W.; Radzieński, M. Structural health monitoring by means of elastic wave propagation. J. Phys. Conf. Ser. 2012, 382, 012003. [Google Scholar] [CrossRef]
  4. Jamaludin, N.; Mba, D.; Bannister, R.H. Condition monitoring of slow-speed rolling element bearings using stress waves. Proc. Inst. Mech. Eng. Part E J. Process Mech. Eng. 2001, 215, 245–271. [Google Scholar] [CrossRef]
  5. Lee, B.C.; Staszewski, W.J. Modelling of Lamb waves for damage detection in metallic structures: Part II. Wave interactions with damage. Smart Mater. Struct. 2003, 12, 815. [Google Scholar] [CrossRef]
  6. Palacz, M.; Krawczuk, M. Analysis of longitudinal wave propagation in a cracked rod by the spectral element method. Comput. Struct. 2002, 80, 1809–1816. [Google Scholar] [CrossRef]
  7. Gaul, L.; Bischoff, S.; Sprenger, H.; Haag, T. Numerical and experimental investigation of wave propagation in rod-systems with cracks. Eng. Fract. Mech. 2010, 77, 3532–3540. [Google Scholar] [CrossRef]
  8. Wang, X.X.; Yu, T.; Yan, H.P.; Ding, J.F.; Li, Z.; Qin, Z.Y.; Chu, F.L. Application of stress wave theory for pyroshock isolation at spacecraft-rocket interface. Chin. J. Aeronaut. 2021, 34, 75–86. [Google Scholar] [CrossRef]
  9. Fang, X. A one-dimensional stress wave model for analytical design and optimization of oscillation-free force measurement in high-speed tensile test specimens. Int. J. Impact Eng. 2021, 149, 103770. [Google Scholar] [CrossRef]
  10. D’Alembert, J. Researches on the curve that a tense cord forms when set into vibration. Hist. Acad. R. Des Sci. BL Berl. 1747, 3, 214–249. [Google Scholar]
  11. Shin, H.; Kim, D. One-dimensional analyses of striker impact on bar with different general impedance. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2020, 234, 589–608. [Google Scholar] [CrossRef]
  12. Love, A.E.H. A Treatise on the Mathematical Theory of Elasticity; Dover Publications: New York, NY, USA, 1944. [Google Scholar]
  13. Yang, H.; Li, Y.; Zhou, F. Propagation of stress pulses in a Rayleigh-Love elastic rod. Int. J. Impact Eng. 2021, 153, 103854. [Google Scholar] [CrossRef]
  14. Yang, H.; Li, Y.; Zhou, F. Stress waves generated in a Rayleigh–Love rod due to impacts. Int. J. Impact Eng. 2022, 159, 104027. [Google Scholar] [CrossRef]
  15. Valsa, J.; Brančik, L. Approximate formulae for numerical inversion of Laplace transforms. Int. J. Numer. Model. Electron. Netw. Devices Fields 1998, 11, 153–166. [Google Scholar] [CrossRef]
  16. Wang, C.Y.; Thang, N.N.; Wang, H. Stress Wave Propagation in a Rayleigh–Love Rod with Sudden Cross-Sectional Area Variations Impacted by a Striker Rod. Sensors 2024, 24, 4230. [Google Scholar] [CrossRef] [PubMed]
  17. Carrier, G.F.; Krook, M.; Pearson, C.E. Functions of a Complex Variable: Theory and Technique; McGraw Hill: New York, NY, USA, 1966. [Google Scholar]
  18. Juraj. Numerical Inversion of Laplace Transforms in Matlab. MATLAB Central File Exchange. 2024. Available online: https://www.mathworks.com/matlabcentral/fileexchange/32824-numerical-inversion-of-laplace-transforms-in-matlab (accessed on 24 January 2024).
  19. Schiff, J.L. The Laplace Transform: Theory and Applications; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  20. Dassault Systèmes. Abaqus Analysis User’s Guide, Volume II: Analysis; Dassault Systèmes: Providence, RI, USA, 2016. [Google Scholar]
Figure 1. Stress wave propagation in a semi-infinite rod under the collinear impact of the striker rod with different impedances.
Figure 1. Stress wave propagation in a semi-infinite rod under the collinear impact of the striker rod with different impedances.
Applsci 14 06523 g001
Figure 2. Contour integration path with residues in the complex plane.
Figure 2. Contour integration path with residues in the complex plane.
Applsci 14 06523 g002
Figure 3. Flowchart of numerical inversion of Laplace transformation.
Figure 3. Flowchart of numerical inversion of Laplace transformation.
Applsci 14 06523 g003
Figure 4. Comparison of the exact solution and the proposed numerical solution of Equation (36).
Figure 4. Comparison of the exact solution and the proposed numerical solution of Equation (36).
Applsci 14 06523 g004
Figure 5. Stress wave propagation in the rods under the impact of a striker rod with the same impedances.
Figure 5. Stress wave propagation in the rods under the impact of a striker rod with the same impedances.
Applsci 14 06523 g005
Figure 6. Comparison of axial displacement history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Figure 6. Comparison of axial displacement history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Applsci 14 06523 g006
Figure 7. Comparison of velocity history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Figure 7. Comparison of velocity history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Applsci 14 06523 g007
Figure 8. Comparison of stress wave history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Figure 8. Comparison of stress wave history between the proposed numerical solution and Yang’s analytical solution at various locations in the rods.
Applsci 14 06523 g008
Figure 9. Stress wave propagation was observed at the contact surface (x = 0 m) with the different impedances of the striker rod (A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s).
Figure 9. Stress wave propagation was observed at the contact surface (x = 0 m) with the different impedances of the striker rod (A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s).
Applsci 14 06523 g009
Figure 10. Comparison of stress wave propagation in the striker rod observed at x = −0.1 m under the impact of the striker rod with different impedances by varying the cross-sectional area (Z2 = 32,087.204 kg/s, A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m).
Figure 10. Comparison of stress wave propagation in the striker rod observed at x = −0.1 m under the impact of the striker rod with different impedances by varying the cross-sectional area (Z2 = 32,087.204 kg/s, A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m).
Applsci 14 06523 g010
Figure 11. Comparison of stress wave propagation in a semi-infinite rod observed at x = 0.5 m under the impact of different impedances of the striker rod (Z2 = 32,087.204 kg/s, A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m).
Figure 11. Comparison of stress wave propagation in a semi-infinite rod observed at x = 0.5 m under the impact of different impedances of the striker rod (Z2 = 32,087.204 kg/s, A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m).
Applsci 14 06523 g011
Figure 12. The 3D finite element model under the initial impact velocity 2v0 using solid elements in Abaqus.
Figure 12. The 3D finite element model under the initial impact velocity 2v0 using solid elements in Abaqus.
Applsci 14 06523 g012
Figure 13. Comparison of stress wave propagation in the semi-infinite rod observed at x = 0.5 m under the impact of the striker rod with different impedances by varying the cross-sectional area (Z2 = 32,087.204 kg/s, A1 = 2A2 = 2 × 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m). The results of no Poisson’s effect are determined based on the research by Shin and Kim [11].
Figure 13. Comparison of stress wave propagation in the semi-infinite rod observed at x = 0.5 m under the impact of the striker rod with different impedances by varying the cross-sectional area (Z2 = 32,087.204 kg/s, A1 = 2A2 = 2 × 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m). The results of no Poisson’s effect are determined based on the research by Shin and Kim [11].
Applsci 14 06523 g013
Figure 14. Comparison of stress wave propagation in the semi-infinite rod observed at x = 0.5 m under the impact of a striker rod with the same impedance by varying the cross-sectional area and material properties (Z1 = Z2 = 32,087.204 kg/s, A1 = 0.5A2 = 0.5 × 7.069 × 10−4 m2, E2 = 195 GPa, υ1 = υ2 = 0.3, ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m). The results of no Poisson’s effect are determined based on the research by Shin and Kim [11].
Figure 14. Comparison of stress wave propagation in the semi-infinite rod observed at x = 0.5 m under the impact of a striker rod with the same impedance by varying the cross-sectional area and material properties (Z1 = Z2 = 32,087.204 kg/s, A1 = 0.5A2 = 0.5 × 7.069 × 10−4 m2, E2 = 195 GPa, υ1 = υ2 = 0.3, ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m). The results of no Poisson’s effect are determined based on the research by Shin and Kim [11].
Applsci 14 06523 g014
Figure 15. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = Z2 = 32,087.204 kg/s, A1 = A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, impact energy KE = 74.67 J).
Figure 15. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = Z2 = 32,087.204 kg/s, A1 = A2 = 7.069 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, impact energy KE = 74.67 J).
Applsci 14 06523 g015
Figure 16. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = 2Z2 = 64,174.408 kg/s, A1 = 2A2 = 14.138 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, impact energy KE = 149.34 J).
Figure 16. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = 2Z2 = 64,174.408 kg/s, A1 = 2A2 = 14.138 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, impact energy KE = 149.34 J).
Applsci 14 06523 g016
Figure 17. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = 0.5Z2 =16,043.602 kg/s, A1 = 0.5A2 = 3.5345 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, input energy = 37.34 J).
Figure 17. Stress wave propagation histories in the striker rod and the semi-infinite rod at various locations (Z1 = 0.5Z2 =16,043.602 kg/s, A1 = 0.5A2 = 3.5345 × 10−4 m2, E1 = E2 = 195 GPa, υ1 = υ2 = 0.3, ρ1 = ρ2 = 7850 kg/m3, 2v0 = 11.6 m/s, L = 0.2 m, input energy = 37.34 J).
Applsci 14 06523 g017
Table 1. Parameters of material properties and geometries of the striker rod and the semi-infinite rod under an initial impact velocity of 2v0.
Table 1. Parameters of material properties and geometries of the striker rod and the semi-infinite rod under an initial impact velocity of 2v0.
ParametersStriker RodSemi-Infinite Rod
Diameter, D (mm)30 (A1/A2 = 1)30
Young’s modulus, E (GPa)195 (E1/E2 = 1)195
Mass density, ρ (kg/m3)7850 (ρ1/ρ2 = 1)7850
Poisson’s ratio, υ0.30.3
Rod length, L (m)0.2
Initial velocity, 2v0 (m/s)11.60
Wave speed, c (m/s)5782.695782.69
Impedance, Z = Aρc (kg/s)32,087.20432,087.204
Table 2. Parameters of the striker rod and the semi-infinite rod used in the model under the impact of different impedances.
Table 2. Parameters of the striker rod and the semi-infinite rod used in the model under the impact of different impedances.
ParametersStriker RodSemi-Infinite Rod
Diameter, D (mm)21.213 (A1/A2 = 0.5)
30 (A1/A2 = 1)
42.426 (A1/A2 = 2)
30
Young’s modulus, E (GPa)195195
Mass density, ρ (kg/m3)78507850
Poisson’s ratio, υ0.30.3
Yield stress, σ0 (MPa)12751275
Table 3. Stress variation in rods without Poisson’s effect (Equations (48) and (49)) over the cycles at x = 0 m with different impedances in the striker rod and the semi-infinite rod under the initial impact velocity 2v0.
Table 3. Stress variation in rods without Poisson’s effect (Equations (48) and (49)) over the cycles at x = 0 m with different impedances in the striker rod and the semi-infinite rod under the initial impact velocity 2v0.
CycleZ1 = 0.5Z2Z1 = Z2Z1 = 2Z2
Time
(μs)
σ1
(MPa)
σ2
(MPa)
Time
(μs)
σ1
(MPa)
σ2
(MPa)
Time
(μs)
σ1
(MPa)
σ2
(MPa)
10~69.17 −351.05−175.520~69.17 −263.286−263.2860~69.17−175.52−351.05
269.17~138.34117.0258.50869.17~138.340069.17~138.34−58.51−117.02
3138.34~207.52−39.01−19.50138.34~207.5200138.34~207.52−19.50−39.01
4207.52~276.6913.006.50207.52~276.6900207.52~276.69−6.501−13.00
5276.69~345.86−4.33−2.17276.69~345.8600276.69~345.86−2.17−4.33
Table 4. Stress variation in rods with υ = 0.3 and the solution without Poisson’s effect (Equation (48)) over the cycles at x = −0.1 m with different impedances, under the initial impact velocity 2v0 = 11.6 m/s, Z2 = 32,087.204 kg/s.
Table 4. Stress variation in rods with υ = 0.3 and the solution without Poisson’s effect (Equation (48)) over the cycles at x = −0.1 m with different impedances, under the initial impact velocity 2v0 = 11.6 m/s, Z2 = 32,087.204 kg/s.
CycleZ1 = 0.5Z2Z1 = Z2Z1 = 2Z2
Time
(μs)
σ1
(MPa)
Time
(μs)
σ1
(MPa)
Time
(μs)
σ1
(MPa)
11/217.29~51.88−351.0517.29~51.88−263.28617.29~51.88−175.52
1/251.88~86.46051.88~86.46051.88~86.460
21/286.46~121.05117.0286.46~121.05086.46~121.05−58.51
1/2121.05~155.640121.05~155.640121.05~155.640
31/2155.64~190.22−39.01155.64~190.220155.64~190.22−19.5
1/2190.22~224.810190.22~224.810190.22~224.810
41/2224.81~259.3913.00224.81~259.390224.81~259.39−6.50
1/2259.39~293.980259.39~293.980259.39~293.980
51/2293.98~328.57−4.33293.98~328.570293.98~328.57−2.17
1/2328.57~363.150328.57~363.150328.57~363.150
Table 5. Stress variation in rods without Poisson’s effect (Equation (49)) over the cycles at x = 0.5 m with different impedances in the semi-infinite rod under the initial impact velocity 2v0.
Table 5. Stress variation in rods without Poisson’s effect (Equation (49)) over the cycles at x = 0.5 m with different impedances in the semi-infinite rod under the initial impact velocity 2v0.
CycleZ1 = 0.5Z2Z1 = Z2Z1 = 2Z2
Time
(μs)
σ2
(MPa)
Time
(μs)
σ2
(MPa)
Time
(μs)
σ2
(MPa)
186.46~155.64−175.5286.46~155.64−263.28686.46~155.64−351.05
2155.64~224.8158.51155.64~224.810155.64~224.81−117.02
3224.81~293.98−19.50224.81~293.980224.81~293.98−39.01
4293.98~363.156.50293.98~363.150293.98~363.15−13.00
5363.15~432.32−2.167363.15~432.320363.15~432.32−4.33
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Thang, N.N.; Wang, C.-Y. Stress Wave Propagation in a Semi-Infinite Rayleigh–Love Rod under the Collinear Impact of a Striker Rod with Different General Impedances. Appl. Sci. 2024, 14, 6523. https://doi.org/10.3390/app14156523

AMA Style

Thang NN, Wang C-Y. Stress Wave Propagation in a Semi-Infinite Rayleigh–Love Rod under the Collinear Impact of a Striker Rod with Different General Impedances. Applied Sciences. 2024; 14(15):6523. https://doi.org/10.3390/app14156523

Chicago/Turabian Style

Thang, Nguyen Ngoc, and Chung-Yue Wang. 2024. "Stress Wave Propagation in a Semi-Infinite Rayleigh–Love Rod under the Collinear Impact of a Striker Rod with Different General Impedances" Applied Sciences 14, no. 15: 6523. https://doi.org/10.3390/app14156523

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