# The Analysis of Stress Waves at a Junction of Beam and String.

1. Introduction

Complex structures, such as arch bridge, suspension bridge, and cable-stayed bridge, are used increasingly in the world. Different from the traditional theory, the vehicle load is regarded as the vibration load in the bridge structure. The theory about stress waves is introduced to study the influence of the vibration load. The propagation speed of longitudinal waves and transverse waves in the concrete and steel is in the range of 1000 m/s to 5000 m/s, and in those large structures, the time scale of vehicle load and stress wave in balance are in the same order of magnitude. So, different from the traditional theory of structural dynamics that ignores stress wave propagation and interaction before the static equilibrium, the spread of the stress process is what we concern. Investigation of stress waves at a junction composed of beams and strings is the most important and foundational work.

Lee and Kolsky [1] investigate the generation of stress pulses at the junction of two noncollinear rods. They used the simple one-dimensional theory to describe the propagation of longitudinal waves and used the Timoshenko equation for flexural propagation. The analytic relations between the incident pulse and the four generated pulses are derived both for an incident longitudinal pulse and an incident flexural pulse. Atkins and Hunter [2] study the propagation of longitudinal elastic waves around right-angled corners in rods of square cross section. In their work, the rod is assumed to behave as a rigid body, and it was demonstrated that longitudinal elastic waves in rods of square cross section are transmitted in part around right-angled corners. Doyle [3] used the fast Fourier transform algorithm that allows the force history to be determined to establish a strain relation for the transverse impact of beams.

Mace [4] studied wave reflection and transmission in beams. In his paper, it is seen that incident near fields can give rise to substantial propagation components. Doyle and Kamle [5, 6] had taken an experimental study of the reflection and transmission of flexural waves at an arbitrary T-joint. It could be obtained from their study that the dynamic response is not sensitive to particular values of the joint model, and it depends primarily on the member arrangement and only secondarily on the particular shape and mass of the joint. Horner and White [7] predicted vibrational power transmission through bends and joints in beam-like structures. They built some models which determine the wave type which carries most power in each section of the system. By establishing the wave type, it is then possible to apply the most suitable vibration control technique. Sack [8] studied transverse oscillations in travelling strings. He found that the frequency response to a simple harmonic disturbance and the envelope of the oscillations are similar to those of a clamped string, but the phase of oscillation varies from point to point so that the instantaneous configuration of the string is not a sine curve. Little and Heywoodb [9] had described the axial stress waves in an elastic-perfectly-plastic bar. Different from the classical elastic solution, the load is greater than half the yield load in their study. The paper describes and verifies the model and presents the predictions that it gives for various levels of load. Chen [10] had thoroughly reviewed the research of transverse vibrations of axially moving strings and their control. Chen presents the governing equations with large amplitude and reviews progress on discretized or direct approximate analytical analyses and numerical approaches based on the Galerkin method or the finite difference method. Guo et al. [11] investigated the elastic wave and energy propagation in angled beams. Experimental and numerical methods were taken. It could demonstrate that, after several reflections from and transmission across the bends, energy is progressively smeared throughout the entire beam and does not concentrate at any particular segment. Li et al. [12] developed a fractal damage joint model based on the fractal damage theory to investigate the transmission and reflection of stress waves across joints. The fractal geometrical characteristics of joint surfaces are investigated by using a laser profilometer to scan the joint surfaces in their study. Lancioni et al. [13] studied the scattering of seismic waves induced by the presence of underground cavities in homogeneous soils by using the numerical method. Singh and Goya [14] established the governing equations of a transversely isotropic microstretch material, which are specialized in x-z plane. Plane wave solutions of these governing equations result in a biquadratic velocity equation. They found that, for a specific material, numerical simulation in the presence as well as in the absence of microstretch shows that the coupled longitudinal displacement (CLD) wave is the fastest wave and the coupled transverse microrotational (CTM) wave is the observed slowest wave. Piccolroaza et al. [15] presented a novel analysis of Floquet-Bloch flexural waves in a periodic lattice-like structure consisting of flexural beam ligaments. A special feature of this structure is in the presence of the rotational inertia, which is commonly neglected in conventional models of the Euler-Bernoulli type. Numerical examples were taken to illustrate directional localization, negative refraction, localization at an interface, and neutrality for propagating plane waves across a structured interface for a frequency range corresponding to a Dirac cone.

The previous studies focused on the distribution of the stress wave in the junction of the beams. According to the need of practical engineering, the distribution of the stress wave that was produced by transverse impact in the junction of the beam and string is the focus of the study in this paper.

2. Theoretical Analysis

2.1. Wave Propagation through the Junction. We will here consider the reflection of a stress pulse at the boundary of beams and strings of different diameters and materials (Figure 1(a)). There is an angle [beta] between two beams and [alpha] between beam 1 and string 2. We call the flexural displacement due to the incident wave [v.sub.0]. We denote the displacement of the reflected longitudinal wave by [u.sub.1], the lateral displacement produced by the flexural wave by [v.sub.1] in beam 1, the longitudinal displacement in string 2 by [u.sub.2], the lateral displacement in string 2 by [v.sub.2], the longitudinal displacement in beam 3 by [u.sub.3], and the lateral displacement in beam 3 by [v.sub.3]. The joint is modeled as a rigid body, and the internal forces are shown in Figure 1(b).

Continuity of the displacement of beam 1 and string 2 leads to

[u.sub.1] = [u.sub.2]cos[alpha] - [v.sub.2]sin[alpha], (1)

[v.sub.0] - [v.sub.i] = [u.sub.2]sin[alpha] + [v.sub.2]cos[alpha]. (2)

Continuity of the displacement of beams 1 and 3 leads to

[u.sub.1] = [u.sub.3]cos[beta] - [v.sub.3]sin[beta] + [L/2]sin[beta][[partial derivative][v.sub.3]/[partial derivative][x.sub.3]], (3)

[v.sub.0] - [v.sub.1] = [u.sub.3]sin[beta] + [v.sub.3]cos[beta] - [L/2](1 + cos[beta])[[partial derivative][v.sub.3]/[partial derivative][x.sub.3]], (4)

[[partial derivative][v.sub.0]/[partial derivative][x.sub.1]] - [[partial derivative][v.sub.1]/[partial derivative][x.sub.1]] = -[[partial derivative][v.sub.3]/ [partial derivative][x.sub.3]]. (5)

The balance of the forces at the junction leads to

[mathematical expression not reproducible], (6)

[mathematical expression not reproducible], (7)

[mathematical expression not reproducible], (8)

where [E.sub.1][I.sub.1] and [E.sub.3][I.sub.3] are bending stiffness of beams 1 and 3. [E.sub.1][A.sub.1], [E.sub.2][A.sub.2], and [E.sub.3][A.sub.3] are axial stiffness of beams 1 and 3 and string 2. L is a parameter characterizing the size of the joint. [I.sub.j] and [m.sub.j] are the moment of inertia and quality of the joint.

2.2. Characterization of the Wave. In the paper by Doyle and Kamle [5, 6], the fast Fourier transform (FFT) numerical algorithm was used to express an incident pulse of arbitrary shape. The FFT is defined as

[mathematical expression not reproducible], (9)

where x(n) is the value of the signal at time (n - 1)[DELTA]t and X(k) is the FFT at frequency k/N[DELTA]t. Different from the previous study, the Timoshenko beams are used in this paper. The displacement can be written as

[mathematical expression not reproducible], (10)

where

[mathematical expression not reproducible]. (11)

The displacements and strains can be written as

[mathematical expression not reproducible], (12)

where [A.sub.n], [B.sub.n], and [k.sub.1], [k.sub.2] are understood to be frequency dependent. If the FFT is taken of the strain history, then it can be approximated with a finite set of frequency components, that is,

[mathematical expression not reproducible], (13)

where [??](x, [w.sub.n]) is the FFT component at frequency [w.sub.n] of the strain history at x. If the strain is sampled at two locations, then [A.sub.n] and [B.sub.n] at discrete frequency can be solved as follows:

[mathematical expression not reproducible]. (14)

This is different from that reported by Doyle and Kamle [5, 6]. They assumed that the accuracy of the simpler Bernoulli-Euler model was up to the requirements of the study. This is generally true when the wavelength is larger than the size of the beam section. However, if the wavelength and the size of the beam section are approximately equal, then the Timoshenko model can get more accurate results.

Finally, the incident and transmitted displacements are expanded in the forms:

[mathematical expression not reproducible], (15)

[mathematical expression not reproducible], (16)

[mathematical expression not reproducible], (17)

[mathematical expression not reproducible], (18)

[mathematical expression not reproducible], (19)

[mathematical expression not reproducible], (20)

where [c.sub.2] and [c'.sub.2] are velocities of the longitudinal and flexural waves in string 2. The hypotheses that [v.sub.2] is small and perpendicular to the x-axis at any time are taken in this paper. So, the tension T of the string does not change over time. [c.sub.2] and [c'.sub.2] can be expressed in the forms:

[c.sub.2] = [square root of [E.sub.2]/[[rho].sub.2]],

[c'.sub.2] = [square root of T/[[rho].sub.2]]. (21)

[c.sub.3] in (19) is the velocity of the longitudinal wave and it can be found using the following equation:

[c.sub.3] = [square root of [E.sub.3]/[[rho].sub.3]]. (22)

For convenience, the joint is located at x = 0. From (1)-(8), the eight unknowns [D.sub.1n], [D.sub.2n], ..., [D.sub.8n] can be determined in terms of the coefficients [A.sub.n] and [B.sub.n] which are given by the initial conditions. Due to the limitation of the space, the expression of the coefficients will not be given in this article.

3. Numerical Validation of the Cosine Wave

According to the theory of Fourier transform, the arbitrary wave produced by vehicle vibration load can be explained as a superposition of the different frequency and amplitude cosine waves. In order to simplify calculation, consider a beam-string system subjected to a cosine wave. The length of the string and beam is 1 m. A square whose side length is 25 mm is used as a beam cross section, and a circle whose radius is 1 mm is used as a string cross section. The top of the string and the right of the beam are fixed. The load is applied in the left side of the beam. In order to simplify the process of calculation, a cycle of cosine wave is calculated in this article, as shown in Figure 2(a). The circular frequency and the amplitude of the wave are 2000[pi] and 3 mm. For convenience of the calculation, two beams adopt the same material and cross section. The beam is made of the steel with Young's modulus E = 207,800 N/[mm.sup.2], Poisson's ratio = 0.3, and density = 7.8[e.sup.-9]t/[mm.sup.3]. The string is made of the steel with Young's modulus E= 195,000 N/[mm.sup.2], Poisson's ratio = 0.3, and density = 7.8[e.sup.-9]t/[mm.sup.3].

The angle [alpha] changed with the different mechanical system of the bridge. To investigate the angle [alpha] between the beam and string on the distribution of the force in the joint, three angles corresponding to [pi]/2, [pi]/3, and [pi]/6 are investigated, respectively. The angle [beta] is kept as zero. Figure 2(b) shows the finite element of the model when [alpha] = [pi]/2. The right end of the beam and the left end of the cable are consolidated. The load is applied to the left of the beam. Because the section size of the string is much smaller and the string is always tensioned under the cosine load, a two-node linear Timoshenko beam element (B31) was used and each bent beam and the string were discretized into 100 elements of uniform size.

Figures 3(a)-3(c) show the difference between the two methods (theoretical analysis and finite element analysis by Abaqus programs) in three locations (the left and right of the joint and the string). We choose the half cycle when the string bears the pulling force to discuss. The phase is ignored because it is not our concern, and the time domain is from 0 to 0.5 ms. It can be seen from Figure 3 that the amplitudes of [Q.sub.1], [N.sub.2], and [Q.sub.3] are 23940 N, 405 N, and 23220 N by theoretical analysis and 20836 N, 328 N, and 20128 N by finite element analysis. It is noticed that the amplitude of the internal force by finite element analysis is slightly smaller than that by theoretical analysis. This is not surprising because of the loss in the process of the wave propagation and the error of the algorithm. It can be seen from the results of the comparison that the analytical solutions are in good agreement with finite element analysis prediction.

Figure 4 illustrates the effect of a on the internal force of the joint. It is shown in Figures 4(a), 4(b), 4(e), and 4(f) that the effect of [alpha] on the shearing forces on both sides of the joint is minor compared to that on the pulling force of the string whatever the analytical method is. It means that the pulling force [N.sub.2] exhibits a great relevance to the angle [alpha]. When [alpha] = [pi]/2, [pi]/3, and [pi]/6, [N.sub.2] = 405, 352, and 212 N by theoretical analysis and 328,268, and 148 N by finite element analysis. It can be concluded that the pulling force has a nonlinear relationship with the angle a. It can also be observed that the rate of change of the pulling force [N.sub.2] by theoretical analysis is close to that by finite element analysis. This can also prove the accuracy of the theoretical method on the other side.

The energy density is another important research object in this paper. [alpha] = [pi]/2, [beta] = 0, and the cosine wave are taken to study the influence of the radius of the string on the energy density. The energy density can be calculated from mechanical energy as follows:

[mathematical expression not reproducible], (23)

[mathematical expression not reproducible], (24)

[mathematical expression not reproducible], (25)

[mathematical expression not reproducible], (26)

[mathematical expression not reproducible], (27)

[mathematical expression not reproducible], (28)

[mathematical expression not reproducible], (29)

[mathematical expression not reproducible], (30)

[mathematical expression not reproducible], (31)

where [W.sub.k], [W.sub.p], and w are the kinetic energy, potential energy, and energy density. [xi] is the quality of the unit length. [DELTA]V is the unit volume. The subscript numbers 1, 2, and 3 represent the locations on the left side of the joint, the string, and the right of the joint. From (25), (28), and (31), the energy densities [w.sub.1], [w.sub.2], and [w.sub.3] can be obtained by the coefficients A, [D.sub.2], [D.sub.3], [D.sub.4], [D.sub.7], and [D.sub.8].

Figure 5 illustrates the influence of the radius of the string on the energy density. It is shown in Figure 5 that the energy density [w.sub.3] is much larger than that in locations 1 and 2. It can be explained that the reflected wave [v.sub.1] has an unloading effect on the incident wave [v.sub.0]. It also can be observed that the energy density [w.sub.3] varies little with the radius of the string from Figure 5(c), and that the energy densities [w.sub.1] and [w.sub.2] decrease with the increase of the radius of the string from Figures 5(a) and 5(b). Such a discrepancy is attributed to the sensitivity to the tensile strength of the string. As a matter of fact, with the increase of the radius, the string should take on more energy. This can also explain the phenomenon that the energy density [w.sub.1] decreases. But the growth rate of the tensile strength of the string is much greater than that of the energy they sustain. This indicates that increasing the rigidity of the string properly can make the structure more reasonable. In practical engineering, the choice of the string will also be limited by other factors. Therefore, meeting the various conditions of the diameter of the string is the best choice.

Referring to (25), the expressions of the energy density [w.sub.4] caused by the longitudinal wave [u.sub.1], [w.sub.5] caused by the flexural wave [v.sub.2], and [w.sub.6] caused by the longitudinal wave [u.sub.3] are as follows:

[mathematical expression not reproducible]. (32)

Figure 6 plots the influence of [alpha] on the energy density when the radius of the string r = 1 mm. It can be seen from Figure 6(a) that the energy density [w.sub.4] by the flexural wave [v.sub.0] - [v.sub.1] decreases with the increase of the angle [alpha]. The energy density [w.sub.2] by the longitudinal wave [u.sub.2] is opposite to the change in energy density [w.sub.1]. The energy density [w.sub.3] by the flexural wave [v.sub.3] is not influenced by angular change at all. It can be observed from Figures 6(d) and 6(f) that the energy densities [w.sub.4] and [w.sub.5] by the longitudinal wave reach maximum when the angle [alpha] is [pi]/4. The reason is that the mechanical energy reaches maximum in that angle. The effect of the longitudinal wave on the whole energy is very small because the energy density caused by the longitudinal wave is two orders of magnitude smaller than that caused by the flexural wave. From Figure 6(e), we can obtain that the energy density [w.sub.5] caused by the flexural wave [v.sub.2] decreases with the increase of the angle. On the contrary, the energy density [w.sub.1] increases with the increase of the angle, as shown in Figure 6(b). So, we can come to a conclusion that the flexural wave is the main wave pattern in the string when the angle is relatively small; otherwise, it is the longitudinal wave.

4. Conclusions

The stress waves at a joint of a structure composed of the string and beam are investigated in this paper. The governing equations are established by displacement balance and internal force balance. The fast Fourier transform (FFT) numerical algorithm is used to deal with the incident wave of arbitrary wave type. The unknown coefficients of the reflected wave and transmitted wave are obtained by governing equations. The analytical solutions are compared with the finite element analysis by Abaqus program, and some parameters are analyzed in this paper. We can yield the following conclusions: first, the results of theoretical analysis are slightly larger than those of the finite element analysis, and the error is about 10 percent. Because the wave is dissipative in the finite element analysis, the result is reasonable, and the theoretical analysis is correct. Second, through the analysis of the influence of the angle [alpha] on the internal force and energy density, we can come to a conclusion that the flexural wave is the main wave pattern in the string when the angle is relatively small; otherwise, it is the longitudinal wave. The last conclusion is that the tensile strength of the string has a great influence on the stress of the structure. If the tensile strength is relatively small, then the energy density of the string and beam will be relatively large, which is not conducive to the longterm use of the structure. This is important for practical engineering and further study.

https://doi.org/10.1155/2018/4863503

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

All the authors gratefully acknowledge the support of the National Natural Science Foundation of China (nos. 51408228 and 51378220).

References

[1] J. P. Lee and H. Kolsky, "The generation of stress pulses at the junction of two noncollinear rods," Journal of Applied Mechanics, vol. 39, no. 3, pp. 809-813, 1972.

[2] K. J. Atkins and S. C. Hunter, "The propagation of longitudinal elastic waves around right-angle corners in rods of square cross-section," Quarterly Journal of Mechanics and Applied Mathematics, vol. 28, no. 2, pp. 245-260, 1974.

[3] J. F. Doyle, "Further developments in determining the dynamic contact," Experimental Mechanics, vol. 24, no. 4, pp. 65-270, 1984.

[4] B. R. Mace, "Wave reflection and transmission in beams," Journal of Sound and Vibration, vol. 97, no. 2, pp. 237-246, 1984.

[5] J. F. Doyle and S. Kamle, "An experimental study of the reflection and transmission of flexural waves at discontinuities," Journal of Applied Mechanics, vol. 52, no. 3, pp. 669-673, 1985.

[6] J. F. Doyle and S. Kamle, "An experimental study of the reflection and transmission of flexural waves at an arbitrary T-joint," Journal of Applied Mechanics, vol. 54, no. 1, pp. 136-140, 1987.

[7] J. L. Horner and R. G. White, "Prediction of vibrational power transmission through bends and joints in beam-like structures," Journal of Sound and Vibration, vol. 147, no. 1, pp. 87-103, 1991.

[8] R. A. Sack, "Transverse oscillations in travelling strings," British Journal of Applied and Physics, vol. 5, no. 5, pp. 224-226, 1991.

[9] G. H. Little and M. D. Heywoodb, "Axial stress waves in an elastic-perfectly-plastic bar," International Journal of Mechanical Sciences, vol. 43, no. 12, pp. 2791-2813, 2001.

[10] L. Q. Chen, "Analysis and control of transverse vibrations of axially moving strings," Applied Mechanics Reviews, vol. 58, no. 2, pp. 91-116, 2005.

[11] Y. B. Guo, V. P. W. Shim, and A. Y. L. Yeo, "Elastic wave and energy propagation in angled beams," Acta Mechanics, vol. 214, no. 1-2, pp. 79-94, 2010.

[12] Y. X. Li, Z. M. Zhu, B. X. Li, J. H. Deng, and H. P. Xie, "Study on the transmission and reflection of stress waves across joints," International Journal of Rock Mechanics & Mining Sciences, vol. 48, no. 3, pp. 364-371, 2011.

[13] G. Lancioni, R. Bernetti, E. Quagliarini, and L. Tonti, "Effects of underground cavities on the frequency spectrum of seismic shear waves," Advances in Civil Engineering, vol. 2014, Article ID 934284, 31 pages, 2014.

[14] B. Singh and M. Goya, "Wave propagation in a transversely isotropic microstretch elastic solid," Mechanics of Advanced Materials and Modern Processes, vol. 3, no. 1, p. 10, 2017.

[15] A. Piccolroaza, A. B. Movchanb, and L. Cabrasa, "Dispersion degeneracies and standing modes in flexural waves supported by Rayleigh beam structures," International Journal of Solids and Structures, vol. 109, pp. 152-165, 2017.

Mu Chen [ID], (1) Ronghui Wang, (1) Niujing Ma [ID], (1) and Zhuojie Zhang (2)

(1) School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, China

(2) School of Civil Engineering, Shijiazhuang Tiedao University, Shjiazhuang 050043, China

Correspondence should be addressed to Niujing Ma; maniujing@yahoo.com

Received 3 September 2017; Revised 14 November 2017; Accepted 19 December 2017; Published 7 March 2018

Caption: Figure 1: The displacement and stress of the three-member intersection. (a) The displacement of the components. (b) The internal forces of the junction.

Caption: Figure 2: A cycle of cosine load (a) and the finite element of the model (b).

Caption: Figure 3: The internal force by two methods ([alpha] = [pi]/2). (a) Reflected wave. (b) The string. (c) The transmission wave.

Caption: Figure 4: Effect of [alpha] on the internal force of the joint ([beta] = 0). (a) [Q.sub.1] by theoretical analysis. (b) [Q.sub.1] by finite element analysis. (c) [N.sub.2] by theoretical analysis. (d) [N.sub.2] by finite element analysis. (e) [Q.sub.3] by theoretical analysis. (f) [Q.sub.3] by finite element analysis.

Caption: Figure 5: The influence of the radius of the string on the energy density ([alpha] = [pi]/2, [beta] = 0). (a) The energy density caused by the flexural wave [v.sub.0] - [v.sub.1]. (b) The energy density caused by the longitudinal wave [u.sub.2]. (c) The energy density caused by the flexural wave [v.sub.3].

Caption: Figure 6: The influence of a on the energy density (r = 1 mm). (a) The energy density by the flexural wave [v.sub.0] - [v.sub.1]. (b) The energy density by the longitudinal wave [u.sub.2]. (c) The energy density by the longitudinal wave [v.sub.3]. (d) The energy density by the longitudinal wave [u.sub.1]. (e) The energy density by the flexural wave [v.sub.2]. (f) The energy density by the longitudinal wave [u.sub.3].