ACE Advances in Civil Engineering 1687-8094 1687-8086 Hindawi 10.1155/2018/4863503 4863503 Research Article The Analysis of Stress Waves at a Junction of Beam and String http://orcid.org/0000-0002-7350-6198 Chen Mu 1 Wang Ronghui 1 http://orcid.org/0000-0002-1855-2773 Ma Niujing 1 Zhang Zhuojie 2 Dede Tayfun 1 School of Civil Engineering and Transportation South China University of Technology Guangzhou 510640 China scut.edu.cn 2 School of Civil Engineering Shijiazhuang Tiedao University Shjiazhuang 050043 China stdu.edu.cn 2018 732018 2018 03 09 2017 14 11 2017 19 12 2017 732018 2018 Copyright © 2018 Mu Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In the bridge engineering, there are some problems about the dynamics that traditional theory cannot solve. So, the theory about stress waves is introduced to solve the related problems. This is a new attempt that the mechanic theory is applied to practical engineering. The stress wave at a junction of the structure composed of beams and strings is investigated in this paper. The structure is studied because the existence of a soft rope makes the transmission of the force in the bridge structure different from the traditional theory, and it is the basis for further research. The equilibrium equations of the displacement and the internal force are built based on the hypothesis. The fast Fourier transform (FFT) numerical algorithm is used to express an incident pulse of arbitrary shape. The analytical solutions are substantiated by comparing with the finite element programs. The conclusion that if the cross section of the string is relatively small, then the energy density of the structure is relatively large, which is disadvantageous to the structure, can be obtained from this paper.

National Natural Science Foundation of China 51408228 51378220
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  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  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  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  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  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  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  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  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.  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.  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.  studied the scattering of seismic waves induced by the presence of underground cavities in homogeneous soils by using the numerical method. Singh and Goya  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.  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 β between two beams and α between beam 1 and string 2. We call the flexural displacement due to the incident wave v0. We denote the displacement of the reflected longitudinal wave by u1, the lateral displacement produced by the flexural wave by v1 in beam 1, the longitudinal displacement in string 2 by u2, the lateral displacement in string 2 by v2, the longitudinal displacement in beam 3 by u3, and the lateral displacement in beam 3 by v3. The joint is modeled as a rigid body, and the internal forces are shown in Figure 1(b).

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

Continuity of the displacement of beam 1 and string 2 leads to(1)u1=u2cosαv2sinα,(2)v0v1=u2sinα+v2cosα.

Continuity of the displacement of beams 1 and 3 leads to(3)u1=u3cosβv3sinβ+L2sinβv3x3,(4)v0v1=u3sinβ+v3cosβL21+cosβv3x3,(5)v0x1v1x1=v3x3.

The balance of the forces at the junction leads to(6)E1I12v0v1x12L2E1I13v1x13=E3I33v3x33+L2E3I33v3x33Ijv1x1,(7)E1A1u1x1=E2A2u2x2cosα+E3A3u3x3cosβ+E3I33v3x33sinβmj2u1t2,(8)E1I13v0v1x13=E2A2u2x2sinαE3A3u3x3sinβE3I33v3x33cosβmju2t,where E1I1 and E3I3 are bending stiffness of beams 1 and 3. E1A1, E2A2, and E3A3 are axial stiffness of beams 1 and 3 and string 2. L is a parameter characterizing the size of the joint. Ij and mj 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(9)Xk=n=0N1xne2πikn/N,n,k=0,1,2,,N1,where x(n) is the value of the signal at time (n − 1)Δt and X(k) is the FFT at frequency k/NΔt. Different from the previous study, the Timoshenko beams are used in this paper. The displacement can be written as(10)v=Aeik1x+Bek2xeiwt,where(11)k1=k1+14wkcs4+12wkcs21/2,k=ρAEI1/4w,k2=k1+14wkcs412wkics21/2,cs=Gρ.

The displacements and strains can be written as(12)vx,t=n=0NAneik1x+Bnek2xeiwt,ex,t=h22vx2=h2n=0NAnk12eik1xBnk22ek2xeiwt,where An, Bn, and k1, k2 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,(13)ex,t=n=0Ne^x,wneiwnt,where e^x,wn is the FFT component at frequency wn of the strain history at x. If the strain is sampled at two locations, then An and Bn at discrete frequency can be solved as follows:(14)An=1k2n2e^x1,wnek2nx1e^x2,wnek2nx2Q,Q=k1n2k2n2eik1nx1k2nx1eik1nx2k2nx2,Bn=1k1n2e^x1,wneik1nx1e^x2,wneik1nx2L,L=k2n2k1n2ek2nx2+ik1nx1ek2nx1+ik1nx1.

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:(15)u1=n=0ND1neiwtx1/c0,(16)v1=n=0ND2neik11x1+D3nek21x1eiwt,(17)u2=n=0ND4neiwtx2/c2,(18)v2=n=0ND5neiwtx2/c2,(19)u3=n=0ND6neiwtx3/c3,(20)v3=n=0ND7neik13x3+D8nek23x3eiwt,where c2 and c2 are velocities of the longitudinal and flexural waves in string 2. The hypotheses that v2 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. c2 and c2 can be expressed in the forms:(21)c2=E2ρ2,c2=Tρ2.

c 3 in (19) is the velocity of the longitudinal wave and it can be found using the following equation:(22)c3=E3ρ3.

For convenience, the joint is located at x = 0. From (1)–(8), the eight unknowns D1n,D2n,,D8n can be determined in terms of the coefficients An and Bn 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π 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/mm2, Poisson’s ratio = 0.3, and density = 7.8e−9 t/mm3. The string is made of the steel with Young’s modulus E = 195,000 N/mm2, Poisson’s ratio = 0.3, and density = 7.8e−9 t/mm3.

A cycle of cosine load (a) and the finite element of the model (b).

The angle α changed with the different mechanical system of the bridge. To investigate the angle α between the beam and string on the distribution of the force in the joint, three angles corresponding to π/2, π/3, and π/6 are investigated, respectively. The angle β is kept as zero. Figure 2(b) shows the finite element of the model when α = π/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 Q1, N2, and Q3 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.

The internal force by two methods (α = π/2). (a) Reflected wave. (b) The string. (c) The transmission wave.

Figure 4 illustrates the effect of α 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 α 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 N2 exhibits a great relevance to the angle α. When α = π/2, π/3, and π/6, N2 = 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 α. It can also be observed that the rate of change of the pulling force N2 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.

Effect of α on the internal force of the joint (β = 0). (a) Q1 by theoretical analysis. (b) Q1 by finite element analysis. (c) N2 by theoretical analysis. (d) N2 by finite element analysis. (e) Q3 by theoretical analysis. (f) Q3 by finite element analysis.

The energy density is another important research object in this paper. α = π/2, β = 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:(23)Wk1=12ξΔxv0v1t2=12ξΔxw2AD2D3eiwt2,(24)Wp1=12c0ξΔxv0v1x12=12c0ξΔxAiwc0+D2ik11+D3k212eiwt2,(25)w1=Wp1+Wk1ΔV=12ρ1AD2D3w2c0Aiwc0+D2ik11+D3k212eiwt2,(26)Wk2=12ξΔxπu2t2=12ρΔxD42w2eiwt2,(27)Wp2=12c22ξΔxu2x22=12ρ2ΔxD42w2eiwt2,(28)w2=Wk1+Wp1ΔV=ρ2D42w2eiwt2,(29)Wk3=12ξΔxv3t2=12ξΔxw2D7+D82eiwt2,(30)Wp3=12c32ξΔxv3x32=12c32ξΔxD7ik13+D8k232eiwt2,(31)w3=Wk3+Wp3ΔV=12ρ3D7+D82w2c32D7ik31+D8k322eiwt2,where Wk, Wp, and w are the kinetic energy, potential energy, and energy density. ξ is the quality of the unit length. Δ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 w1,w2, and w3 can be obtained by the coefficients A, D2, D3, D4, D7, and D8.

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 w3 is much larger than that in locations 1 and 2. It can be explained that the reflected wave v1 has an unloading effect on the incident wave v0. It also can be observed that the energy density w3 varies little with the radius of the string from Figure 5(c), and that the energy densities w1 and w2 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 w1 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.

The influence of the radius of the string on the energy density (α = π/2, β = 0). (a) The energy density caused by the flexural wave v0 − v1. (b) The energy density caused by the longitudinal wave u2. (c) The energy density caused by the flexural wave v3.

Referring to (25), the expressions of the energy density w4 caused by the longitudinal wave u1, w5 caused by the flexural wave v2, and w6 caused by the longitudinal wave u3 are as follows:(32)w4=ρ1D12w2eiwt2,w5=ρ2D52w2eiwt2,w6=ρ1D62w2eiwt2.

Figure 6 plots the influence of α 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 w4 by the flexural wave v0 − v1 decreases with the increase of the angle α. The energy density w2 by the longitudinal wave u2 is opposite to the change in energy density w1. The energy density w3 by the flexural wave v3 is not influenced by angular change at all. It can be observed from Figures 6(d) and 6(f) that the energy densities w4 and w5 by the longitudinal wave reach maximum when the angle α is π/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 w5 caused by the flexural wave v2 decreases with the increase of the angle. On the contrary, the energy density w1 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.

The influence of α on the energy density (r = 1 mm). (a) The energy density by the flexural wave v0 − v1. (b) The energy density by the longitudinal wave u2. (c) The energy density by the longitudinal wave v3. (d) The energy density by the longitudinal wave u1. (e) The energy density by the flexural wave v2. (f) The energy density by the longitudinal wave u3.

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 α 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 long-term use of the structure. This is important for practical engineering and further study.

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).

Lee J. P. Kolsky H. The generation of stress pulses at the junction of two noncollinear rods Journal of Applied Mechanics 1972 39 3 809 813 10.1115/1.34227932-s2.0-0015397565 Atkins K. J. Hunter S. C. The propagation of longitudinal elastic waves around right-angle corners in rods of square cross-section Quarterly Journal of Mechanics and Applied Mathematics 1974 28 2 245 260 10.1093/qjmam/28.2.2452-s2.0-0016507052 Doyle J. F. Further developments in determining the dynamic contact Experimental Mechanics 1984 24 4 65 270 10.1007/bf023239862-s2.0-0021595541 Mace B. R. Wave reflection and transmission in beams Journal of Sound and Vibration 1984 97 2 237 246 10.1016/0022-460x(84)90320-12-s2.0-0021521981 Doyle J. F. Kamle S. An experimental study of the reflection and transmission of flexural waves at discontinuities Journal of Applied Mechanics 1985 52 3 669 673 10.1115/1.31691192-s2.0-0022132022 Doyle J. F. Kamle S. An experimental study of the reflection and transmission of flexural waves at an arbitrary T-joint Journal of Applied Mechanics 1987 54 1 136 140 10.1115/1.31729472-s2.0-0023306587 Horner J. L. White R. G. Prediction of vibrational power transmission through bends and joints in beam-like structures Journal of Sound and Vibration 1991 147 1 87 103 10.1016/0022-460x(91)90685-d2-s2.0-0026152477 Sack R. A. Transverse oscillations in travelling strings British Journal of Applied and Physics 1991 5 5 224 226 10.1088/0508-3443/5/6/3072-s2.0-0001186943 Little G. H. Heywoodb M. D. Axial stress waves in an elastic–perfectly-plastic bar International Journal of Mechanical Sciences 2001 43 12 2791 2813 10.1016/s0020-7403(01)00069-82-s2.0-0035546016 Chen L. Q. Analysis and control of transverse vibrations of axially moving strings Applied Mechanics Reviews 2005 58 2 91 116 10.1115/1.18491692-s2.0-25844516751 Guo Y. B. Shim V. P. W. Yeo A. Y. L. Elastic wave and energy propagation in angled beams Acta Mechanics 2010 214 1-2 79 94 10.1007/s00707-010-0317-62-s2.0-77957929985 Li Y. X. Zhu Z. M. Li B. X. Deng J. H. Xie H. P. Study on the transmission and reflection of stress waves across joints International Journal of Rock Mechanics & Mining Sciences 2011 48 3 364 371 10.1016/j.ijrmms.2011.01.0022-s2.0-79952625348 Lancioni G. Bernetti R. Quagliarini E. Tonti L. Effects of underground cavities on the frequency spectrum of seismic shear waves Advances in Civil Engineering 2014 2014 31 934284 10.1155/2014/9342842-s2.0-84919740448 Singh B. Goya M. Wave propagation in a transversely isotropic microstretch elastic solid Mechanics of Advanced Materials and Modern Processes 2017 3 1 10 10.1186/s40759-017-0023-3 Piccolroaza A. Movchanb A. B. Cabrasa L. Dispersion degeneracies and standing modes in flexural waves supported by Rayleigh beam structures International Journal of Solids and Structures 2017 109 152 165 10.1016/j.ijsolstr.2017.01.0172-s2.0-85009754306