Shock and Vibration

Volume 2014 (2014), Article ID 530131, 7 pages

http://dx.doi.org/10.1155/2014/530131

Research Article

Analytical Solution for the Sound Radiation Field of a Viscoelastically Supported Beam Traversed by a Moving Load

Rezgar Shakeri and Davood Younesian

Center of Excellence in Railway Transportation, School of Railway Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran

Received 20 October 2012; Accepted 19 November 2012; Published 23 June 2014

Academic Editor: Abdolreza Ohadi

Copyright © 2014 Rezgar Shakeri and Davood Younesian. 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.

Abstract.

Sound radiation from a beam resting on a viscoelastic foundation is analytically studied when it is subjected to a moving load. The topic could cover a range of applications such as submerged floating tunnels, buried pipelines, and railway tracks. Galerkin’s method is employed to obtain the transverse vibration of the beam. Based on the Rayleigh integral approach, acoustic pressure distribution around the beam is obtained in the time domain. In the second part of this paper, corresponding displacement and acoustic pressure are obtained by the use of the Rayleigh-Ritz approach in conjunction with the Laplace transform method and by the use of the Fourier transform, respectively. Durbin’s numerical Laplace transform inversion scheme is eventually employed to obtain dynamic responses. A parametric study is then carried out and influences of the design parameters as well as the loading conditions on the acoustic pressure field are investigated.

1. Introduction

The investigation of the dynamic response of beams resting on viscoelastic foundation is a very interesting problem, with applications in various fields of engineering and technology. In particular, acoustic radiation can be found in a different array of practical applications such as submerged floating tunnels, railway tracks, and civil structures. Numerical, analytical, and experimental studies have been performed in the past decades to study the dynamic analysis of soil-structure interaction problem. Beam on elastic foundation has been studied by several authors [1–4]. Kenney [5] studied the steady-state deformation when the point load is moving with constant speed. Kargarnovin et al. [6–8] and Muscolino and Palmeri [9] investigated the response of beam on viscoelastic foundation under moving loads. Andersen et al. [10], Weitsman [11], Choros and Adams [12], and Lin and Adams [13] investigated the response of Euler-Bernoulli beams on tensionless Winkler foundations. The response of embedded railway track using the Winkler foundation model was studied by Shamalta and Metrikine [14]. Steele [15] and Chen and Huang [16, 17] investigated the response of Timoshenko beam on Winkler foundation for a variety of beam, foundation, and loading conditions. Yang and Ge [18] and Senalp et al. [19] investigated the dynamic behavior of Euler-Bernoulli beam resting on viscoelastic foundation subjected to moving load by using the mode decomposition method together with the precise time integration method (MDPIM). Vlasov and Leont’ev [20] showed that the mechanical behavior of an elastic continuum can be quite accurately simulated using springs with shear interactions between them. The tensile resistance is generated in the ground due to shear strain of the soil. Basu and Kameswara Rao [21] studied the steady-state response of an infinite beam resting on a viscoelastic foundation, modeled using springs with shear interactions and subjected to a concentrated moving load with a constant velocity. In this paper, Galerkin’s method is proposed to investigate the dynamic response of Euler-Bernoulli beam resting on viscoelastic foundation and shear stiffness of soil is considered. Acoustic radiation is examined by using appropriate 2D images of the sound field for selected foundation parameters. This paper is an invited article selected from the ISAV2012 conference proceedings.

2. Formulation

2.1. Analytical Method

Consider a simply supported uniform beam of length , mass per unit length , and bending rigidity , resting on viscoelastic foundation, indicated by a bending layer coefficient , a shear layer coefficient , and a viscous coefficient subjected to a moving load with constant speed along the beam as shown in Figure 1. The equation of motion for the transverse vibration of the beam is given by [21]

where is transverse deflection of the beam and is applied moving concentrated load. The displacement of the beam is expanded using the eigenfunction of the uniform beam as [22]

where is the natural frequency of the th mode and is the eigenfunction. The orthogonality condition is

where is Kronecker’s delta. The transverse displacement can be expanded as follows:

where is the modal coordinate that must be determined by solving an ordinary differential equation. By direct substitution of expansion (4) into (1) and making use of the general orthogonality relation (3) between the distinct normal modes and the fact that the normal modes, , must satisfy the classical relation , (1) can be reduced to

where . The solution of (5) with the zero initial condition can be expressed as

Figure 1: Beam on foundation under moving load.

2.2. Rayleigh-Ritz Method

The Rayleigh-Ritz method is based on an energy or variational principle, such as those provided by the principles of virtual work or their variants, which account for the natural boundary conditions as a part of the principle. The total strain and kinetic energy and the work done by the nonconservative forces are given by [22]

Substituting for , , and from (7) in Hamilton’s principle, we obtain

In the Ritz method, a dependent unknown (e.g., the displacement) is approximated by a finite linear combination as the form [23]

where are time-dependent parameters to be determined for all times and are called the approximation functions. Since the specified geometric boundary condition is homogeneous, we have . The approximation functions must satisfy the simple support’s condition. Direct substituting of (9) into (8) and equating coefficients of to zero may advantageously be put in matrix form as

where

By taking Laplace transform of (10), one obtains

where is the unknown coefficients vector and is the load vector which are given as

And the coefficient matrix is rewritten as

The unknown coefficient can readily be determined from

In the current work, Durbin’s approach [24] for numerical inversion of the Laplace transforms which involves the discretized form of the complex Laplace inversion formula , where and is an arbitrary real number greater than all the real parts of the singularities of and can be or in the interval , shall be adopted. Accordingly, one could readily employ expansion [25]

where is the truncation constant and the suggested value of “” is between 5 and 10 for sufficient accuracy [24].

2.3. Acoustic Radiation

The acoustic pressure radiated from a vibration beam in an infinite baffle can be obtained by evaluating the Rayleigh surface integral where each elemental area on the beam surface is regarded as a simple point source of an outgoing wave and its contribution is added with an appropriate time delay. Referring to Figure 1, the acoustic pressure at the observation point with Cartesian coordinates at time caused by the vibration of the beam is given by [26]

where and are the mass density and wave velocity of the acoustic medium, respectively, is the acceleration time history of the beam obtained previously, and is the distance between the observation point and the point at on the beam surface.

In addition, the pressure distribution in the acoustic fluid domain can be obtained from the solution of the wave equation as follows:

In frequency domain, (18) reduces to , and referring to Figure 1, the Rayleigh integral in the frequency domain [27] is , where and is the acceleration of the beam in the frequency domain. Taking Laplace transform of (18) leads to , and then, the Rayleigh integral in the Laplace domain is

where and is the displacement of the beam in the Laplace domain where is given by (2) and is obtained from (5). Finally, the acoustic pressure can be calculated by (16). In the following section, the comparison between these two methods is considered.

3. Numerical Results

A parametric study is directed in this section to investigate the effects of different parameters on the acoustic pressure field radiated from the vibrating beam series of numerical simulations that have been carried out and some of the numerical results are presented as a few samples in this section. The physical and geometrical properties of the beam and acoustic medium are listed in Table 1 [8].

Table 1: Properties of the UIC60 rail, foundation, load, and acoustic medium [8].

Item
Notation
Value

Young’s modulus (steel)
E
210 GPa

Mass density of beam
ρ
7850 kg/m³

Cross-sectional area
A
7.69 × 10⁻³ m²

Second moment of area
I
30.55 × 10⁻⁶ m⁴

Bending stiffness

202 MN/m²

Viscous damping
c
141 kN s/m²

Load

65 kN

Mass density of acoustic medium (air)

1.2 kg/m³

Wave velocity

343 m/s

Figure 2 exhibits the midspan displacement of the beam with respect to the time. Because of the foundation under the beam, the vibration around the maximum deflection is zero. The maximum deflection will occur at the location of the load. Furthermore, excellent agreement of the present result and [19] for load velocity  m/s is obtained. Figure 3 displays the comparison between the results of the analytical and Rayleigh-Ritz methods for selected load velocity  m/s and for a beam extended 15 m in longitudinal direction. The results show a good agreement between these two methods. Taking shear layer coefficient into consideration, the maximum deflection of the beam will decrease when the shear layer coefficient is increased due to increment of natural frequencies of the beam.

Figure 2: Comparison of the displacement response with [19].

Figure 3: Comparison of analytical and Rayleigh-Ritz methods.

Figure 4 shows the acoustic pressure response at two selected points and for , when the load passes through the beam. The acoustic pressure decreases as the viscoelastic and bending stiffness and height of selected points increase. Also for , an interesting reduction of the first radiated pressure wave from the vibrating beam is seen at .

Figure 4: Acoustic pressure at different heights (midspan).

Figure 5 displays the acoustic pressure response and comparison between two methods for . The results show good agreement between two equations. Equation (17) exhibits instabilities for time interval of ; therefore, it must be calculated by numerical techniques such as standard trapezoidal rule [26]. On the other hand, robustness of (19) was seen in both of the mentioned time spans.

Figure 5: Comparison between two methods for calculating acoustic radiation.

Figure 6 shows the series of time snapshots of acoustic pressure distribution while the moving load is passing on the beam for shear layer coefficient , length of the beam , bending layer coefficient  MN/m², viscous coefficient  kN·sec/m², and typical vehicle load (65 kN). The acoustic pressure fields around the beam show that maximum distribution pressure appears around the location of the moving load which is a result of maximum deflection of the beam.

Figure 6: Acoustic pressure field (Pa) along the beam (15 m) at different positions of moving load.

4. Conclusion

In this paper, transient response of simply supported, Euler-Bernoulli beam resting on viscoelastic foundation subjected to a moving load is studied. Based on the Rayleigh integral and dynamic response of the beam, the acoustic pressure distribution around the beam was obtained for the first time. The time responses of the beam with and without shear stiffness effects of the soil are presented for various foundations of moving load. The effects of and can be easily observed from Figure 4. The most important observations are summarized as follows:

(i)the beam structural vibration is observed to radiate sound into acoustic medium. In particular, increasing the magnitude of foundation is seen to decrease the acoustic pressure;(ii)increasing the height of selected points decreases the acoustic pressure radiation, while the time response keeps the same trend;(iii)the delay of the first acoustic pressure wave radiated from the beam is easily seen;(iv)the contours of acoustic pressure fields around the beam show that maximum distribution pressure appears around the location of the moving load which is a result of maximum deflection of the beam.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

M. Hetenyi, in Beams on Elastic Foundations, vol. 16 of Scientific Series, The University of Michigan Press, University of Michigan Studies, Ann Arbor, Mich, USA, 1946.

L. Andersen, S. R. K. Nielsen, and P. H. Kirkegaard, “Finite element modelling of infinite Euler beams on Kelvin foundations exposed to moving loads in convected co-ordinates,” Journal of Sound and Vibration, vol. 241, no. 4, pp. 587–604, 2001.

D. Thambiratnam and Y. Zhuge, “Dynamic analysis of beams on an elastic foundation subjected to moving loads,” Journal of Sound and Vibration, vol. 198, no. 2, pp. 149–169, 1996.

L. Fryba, Vibration of Solids and Structures Under Moving Loads, Thomas Telford, London, UK, 1999.

J. T. Kenney, “Steady-state vibrations of beam on elastic foundation for moving load,” Journal of Applied Mechanics, vol. 21, pp. 359–364, 1954.

M. H. Kargarnovin, D. Younesian, D. J. Thompson, and C. J. C. Jones, “Response of beams on nonlinear viscoelastic foundations to harmonic moving loads,” Computers and Structures, vol. 83, no. 23-24, pp. 1865–1877, 2005.

M. H. Kargarnovin and D. Younesian, “Dynamics of Timoshenko beams on Pasternak foundation under moving load,” Mechanics Research Communications, vol. 31, no. 6, pp. 713–723, 2004.

D. Younesian and M. H. Kargarnovin, “Response of the beams on random Pasternak foundations subjected to harmonic moving loads,” Journal of Mechanical Science and Technology, vol. 23, no. 11, pp. 3013–3023, 2010.

G. Muscolino and A. Palmeri, “Response of beams resting on viscoelastically damped foundation to moving oscillators,” International Journal of Solids and Structures, vol. 44, no. 5, pp. 1317–1336, 2007.

L. Andersen, S. R. K. Nielsen, and R. Iwankiewicz, “Vehicle moving along an infinite beam with random surface irregularities on a Kelvin foundation,” Journal of Applied Mechanics, Transactions ASME, vol. 69, no. 1, pp. 69–75, 2002.

Y. Weitsman, “Onset of separation between a beam and a tensionless elastic foundation under a moving load,” International Journal of Mechanical Sciences, vol. 13, no. 8, pp. 707–711, 1971.

J. Choros and G. G. Adams, “Steadily moving load on an elastic beam resting on a tensionless winkler foundation,” Journal of Applied Mechanics, Transactions ASME, vol. 46, no. 1, pp. 175–180, 1979.

L. Lin and G. G. Adams, “Beam on tensionless elastic foundation,” Journal of Engineering Mechanics, vol. 113, no. 4, pp. 542–553, 1987.

M. Shamalta and A. V. Metrikine, “Analytical study of the dynamic response of an embedded railway track to a moving load,” Archive of Applied Mechanics, vol. 73, no. 1-2, pp. 131–146, 2003.

C. R. Steele, “The Timoshenko beam with a moving Load,” Journal of Applied Mechanics, vol. 35, no. 3, pp. 481–488, 1968.

Y.-H. Chen and Y.-H. Huang, “Dynamic stiffness of infinite Timoshenko beam on viscoelastic foundation in moving co-ordinate,” International Journal for Numerical Methods in Engineering, vol. 48, no. 1, pp. 1–18, 2000.

Y.-H. Chen and Y.-H. Huang, “Dynamic characteristics of infinite and finite railways to moving loads,” Journal of Engineering Mechanics, vol. 129, no. 9, pp. 987–995, 2003.

Y. Yang and X. Ge, “A new method for vibration response of beam on foundation under moving load,” International Journal of Modern Physics B, vol. 22, no. 31-32, pp. 5615–5620, 2008.

A. D. Senalp, A. Arikoglu, I. Ozkol, and V. Z. Dogan, “Dynamic response of a finite length euler-bernoulli beam on linear and nonlinear viscoelastic foundations to a concentrated moving force,” Journal of Mechanical Science and Technology, vol. 24, no. 10, pp. 1957–1961, 2010.

V. Z. Vlasov and U. N. Leont’ev, Beams, Plates and Shells on Elastic Foundations, Israel Program for Scientific Translations, Jerusalem, Israel, 1966.

D. Basu and N. S. V. Kameswara Rao, “Analytical solutions for Euler-Bernoulli beam on visco-elastic foundation subjected to moving load,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 37, no. 8, pp. 945–960, 2013.

S. S. Rao, Vibration of Continuous Systems, Wiley, New York, NY, USA, 2007.

J. N. Reddy, Theory and Analysis of Elatic Plates and Shells, CRC Press, New York, NY, USA, 2007.

F. Durbin, “Numerical inversion of laplace transforms: an efficient improvement to dubner and abate's method,” Computer Journal, vol. 17, no. 4, pp. 371–376, 1974.

S. M. Hasheminejad, R. Shakeri, and S. Rezaei, “Vibro-acoustic response of an elliptical plate-cavity coupled system to external shock loads,” Applied Acoustics, vol. 73, no. 8, pp. 757–769, 2012.

F. T. K. Au and M. F. Wang, “Sound radiation from forced vibration of rectangular orthotropic plates under moving loads,” Journal of Sound and Vibration, vol. 281, no. 3–5, pp. 1057–1075, 2005.

M. C. Junger and D. Feit, Sound, Structures and Interaction, MIT Press, Cambridge, Mass, USA, 1986.