Dynamic behaviour under moving concentrated masses of simply supported rectangular plates resting on variable Winkler elastic foundation

ARTICLES

Dynamic behaviour under moving concentrated masses of simply supported rectangular plates resting on variable Winkler elastic foundation

S. T. Oni; T. O. Awodola^*

Department of Mathematical Sciences, Federal University of Technology, Akure – Nigeria

ABSTRACT

The response of simply supported rectangular plates carrying moving masses and resting on variable Winkler elastic foundations is investigated in this work. The governing equation of the problem is a fourth order partial differential equation. In order to solve this problem, a technique based on separation of variables is used to reduce the governing fourth order partial differential equations with variable and singular coefficients to a sequence of second order ordinary differential equations. For the solutions of these equations, a modification of the Struble's technique and method of integral transformations are employed. Numerical results in plotted curves are then presented. The results show that response amplitudes of the plate decrease as the value of the rotatory inertia correction factor R₀ increases. Furthermore, for fixed value of R₀, the displacements of the simply supported rectangular plates resting on variable elastic foundations decrease as the foundation modulus F₀ increases. The results further show that, for fixed R₀ and F₀, the transverse deflections of the rectangular plates under the actions of moving masses are higher than those when only the force effects of the moving load are considered. Therefore, the moving force solution is not a safe approximation to the moving mass problem. Hence, safety is not guaranteed for a design based on the moving force solution. Also, the analyses show that the response amplitudes of both moving force and moving mass problems decrease both with increasing Foundation modulus and with increasing rotatory inertia correction factor. The results again show that the critical speed for the moving mass problem is reached prior to that of the moving force for the simply supported rectangular plates on variable Winkler elastic foundation.

Keywords: rectangular plates, winkler foundation, foundation modulus, rotatory inertia, resonance, moving force, moving mass.

1 INTRODUCTION

The analyses of elastic structures (beams, plates and shells), resting on a subgrade, such as railway tracks, highway pavements, navigation locks and structural foundations, constitute an important part of the Civil Engineering and applied Mathematics literatures. In general, such analyses are mathematically complex due to the difficulty in modeling the mechanical response of the subgrade which is governed by many factors. When these structures are acted upon by moving loads, the dynamic analyses of the system become more cumbersome [7]. The crudest approximation known to the literatures to this problem is the so called "moving force" problem, in which the vehicle-track interaction is completely neglected, and the action of the vehicle is described as a concentrated force moving along the beam [12].

Several researchers have considered vehicle-track interaction in their analyses. These they commonly termed moving mass problems. These researchers include Stanisic et al [22], Milornir et al [11], Clastornic et al [3], Sadiku and Leipholz [19] and Gbadeyan and Oni [8]. Douglas et al [5] solved the problem of plate strip of varying thickness and the center of shear. In their work, they considered a free-vibrating strip with classical boundary conditions, precisely, they assumed the plate strip clamped at one end and free at the other end. Pesterev et al [18] came up with a series expansion method for calculating bending moment and shear force in the problem of vibration of a damped beam subject to an arbitrary number of moving loads. This kind of solution, though could be accurate, cannot account for vital information such as the phenomenon of resonance in the dynamical system.

Recently, several other researchers have made tremendous efforts in the study of dynamics of structures under moving loads, these include Oni [13], Oni and Omolofe [17], Oni and Awodola [14], Omer and Aitung [2], Adams [1], Savin [20], Jia-Jang [23]. In all of these, considerations have been limited to cases of one-dimensional (beam) problems. Where two-dimensional (plate) problems have been considered, the foundation moduli are taken to be constants. No considerations have been given to the class of dynamical problems in which the foundation is the type with stiffness variation. In an attempt to solve such two-dimensional problem, all the methods used in the above works break down due to the variation of the foundation model.

It is generally known that the dynamical problems of structures under moving load and resting on a foundation is generally complex, the complexity increases if the foundation stiffness varies along the structure. Aside the problem of singularity brought in by the inclusion of the inertia effects of the moving load, the coefficients of the governing fourth order partial differential equation are no longer constant but variable. Earlier researchers into beam member on variable elastic foundation include Franklin and Scott [6] who presented a closed-form solution to a linear variation of the foundation modulus using contour-integrals. Closely following this, Lentini [10] presented a finite difference method to solve the problem where the foundation stiffness varies along x as a power of x. Much later, Clastornik et al [4] presented a solution for the finite beams resting on a Winkler elastic foundation with stiffness variation that can be presented as a general polynomial of x. Though works in [4, 6, 10] are useful, the loads acting on the beams are not moving loads. In a recent development, Oni and Awodola [15] extended the works of these previous authors to investigate the dynamic response to moving concentrated masses of uniform Rayleigh beams resting on variable Winkler elastic foundation. Their work shows that for all variants of classical boundary conditions, the displacements of a uniform Rayleigh beam resting on variable elastic foundation and traversed by moving masses decrease both with increase in the foundation moduli and the rotatory inertial correction factor.

More recently, Oni and Awodola [16] considered the dynamic response under a moving load of an elastically supported non-prismatic Bernoulli-Euler beam on variable elastic foundation. The technique was based on the generalized Galerkin's method and integral transformations.

In all these previous investigations, extension of the theory to cover two-dimensional (plate) problem has not been effected, when the plate is on variable foundation. Therefore, this study concerns the response to moving concentrated masses of simply supported rectangular plate resting on Winkler elastic foundation with stiffness variation.

2 GOVERNING EQUATION

The equation governing the dynamic transverse displacement W(x, y, t) of a rectangular plate when it is resting on a variable Winkler foundation and traversed by several moving concentrated masses is the fourth order partial differential equation given by

where

is the bending rigidity of the plate, ∇² is the two-dimensional Laplacian operator, W(x, y, t) is the transverse displacement, h is the plate's thickness, E is the Young's Modulus, v is the Poisson's ratio (v < 1), µ is the mass per unit area of the plate, R₀ is the Rotatory inertia correction factor, F(x) is the variable foundation's stiffness, P(x, y, t) is the Moving load, x and y are respectively the spatial coordinates in x and y directions and t is the time coordinate.

When the effect of the mass of the moving load on the response of the plate is taken into consideration, the external moving surface load takes on the form

where P_f(x, y, t) is the continuous moving force, Δ^* is the substantive acceleration operator and g is the acceleration due to gravity.

The structure under consideration is assumed to be carrying an arbitrary number (say N) of concentrated masses M_i moving with constant velocities c_i, i = 1,2,3,..., N along a straight line parallel to the x-axis (no difficulty arises by assuming that masses travel in an arbitrary path) issuing from point y = s on the y-axis. Thus, the moving force acting on the plate is defined as

where δ(.) is the Dirac-Delta function.

The operator Δ^* used in equation (3) for masses traveling in an arbitrary path in the x-y plane is defined as

On the assumptions of the paragraph above, this operator takes the form

As an example in this problem, a variable elastic foundation of the form

where F₀ is the foundation constant, is considered.

Thus, substituting (3), (4), (5) and (7) into (1), one obtains

The initial conditions, without any loss of generality, is taken as

3 ANALYTICAL APPROXIMATE SOLUTION

The method of analysis involves expressing the Dirac-Delta function as a Fourier cosine series. Because of the variable foundation term, the elegant method of the generalized integral transform breaks down while the generalized Galerkin's method used in one-dimensional structural problems (Beam problems) could not handle the two-dimensional structural problem (Plate problems). Hence, the technique based on separation of variables is used to reduce the fourth order partial differential equation governing the motion of the plate to a set of coupled second order ordinary differential equations. Then, the modified asymptotic method of Struble is used to simplify these resulting equations. The method of integral transformation and convolution theory are then employed to obtain the closed form solution of the two-dimensional dynamical problems.

In order to solve equation (8), in the first instance, the deflection is written in the form [21]

where ϕ_n are the known eigenfunctions of the plate with the same boundary conditions. The ϕ_n have the form of

where

Ω_n, n = 1, 2, 3, ... , are the natural frequencies of the dynamical system and T_n(t) are amplitude functions which have to be calculated.

In order to solve the equation (8), it is rewritten as

At this juncture, the right hand side of equation (13) is written in the form of a series, we have

Substituting equation (10) into equation (14) we have

where

Multiplying both sides of equation (15) by ϕ_p(x, y) and integrating on area A of the plate, we have

Considering the orthogonality of ϕ_n(x, y)

where

Using (18), equation (13), taking into account (10) and (11), can be written as

Equation (19) must be satisfied for arbitrary x, y (that is, each point of the plate) and this is possible only when

The system in equation (20) is a set of coupled ordinary differential equations.

Considering the property of the Dirac-Delta function and expressing it in the Fourier cosine series as

and

equation (20) becomes

where

The second order coupled differential equation (23) is the transformed equation governing the problem of a rectangular plate on a variable Winkler elastic foundation. This differential equation holds for all variants of the classical boundary conditions.

In what follows, ϕ_n(x, y) are assumed to be the products of the functions ψ_ni(x) and ψ_nj(y) which are the beam functions in the directions of x and y axes respectively [9]. That is

Since each of these beam functions satisfies all the boundary conditions in its direction, the kernel (the product of these beam functions) in the above integrals satisfies all boundary conditions for any plate problem of practical interest. In particular, these beam functions can be defined respectively, as

and

where A_ni, A_nj, B_ni, B_nj, C_ni and C_nj are constants determined by the boundary conditions. Ω_ni and Ω_nj are called the mode frequencies.

In order to solve equation (23) we shall consider a mass M traveling with uniform velocity c along the line y = s. The solution for any arbitrary number of moving masses can be obtained by superposition of the individual solution since the governing differential equation is linear. Thus for the single mass M₁ equation (23) reduces to

where

Equation (27) is now the fundamental equation of our problem when the rectangular plate has arbitrary end support conditions.

In what follows, we shall solve the equation (27) when the plate has simple supports at all its edges.

An elastic rectangular plate resting on a variable Winkler elastic foundation and having simple supports at all its edges has the boundary conditions given by

Hence for the normal modes

When use is made of equations (31a), (31b), (32a), (32b) and the initial conditions given by equation (9), it can be shown that

Similarly,

Thus, we substitute equations (33), (34), (35) and (36) into the transformed equation (27) to obtain the transformed equation for a rectangular plate, resting on a variable Winkler elastic foundation and having simple supports at all its edges. That is

where

Further simplification and rearrangement of (37), taking into account (33), (34), (35) and (36), yields

where

and

Equation (39) is now the fundamental equation of our problem when the rectangular plate resting on variable Winkler foundation has simple support at all its edges. In what follows, we shall discuss two cases of the equation.

Case I: Simply supported plate traversed by moving force

An approximate model of the system, when the inertia effect of the moving mass M is neglected, that is, when Γ = 0 in equation (39), is the moving force problem associated with the system. Thus the differential equation is given by

To solve equation (42) using the Struble's asymptotic technique [16], first, we neglect the rotatory inertial term and rearrange the equation to take the form

where

Consider a parameter λ < 1 for any arbitrary mass ratio Γ^*, defined as

It can be shown that

Thus, the homogeneous part of equation (43) can be replaced with

where

represents the modified frequency due to the effect of foundation.

Thus using (47), equation (42) can be written as

where

The homogeneous part of equation (49) can be written as

Consider the parameter є₀ < 1 for any arbitrary mass ratio λ₀ defined as

which implies

Following the same argument, (51) can be replaced with

where

Therefore, the moving force problem (42) for the simply supported rectangular plate is reduced to the non-homogeneous ordinary differential equation given as

When equation (56) is solved in conjunction with the initial conditions, one obtains expression for T_n(t). Thus in view of equation (10), one obtains

as the transverse-displacement response to a moving force of a simply supported rectangular plate on a variable Winkler elastic foundation.

Case II: Simply supported rectangular plate resting on variable foundation and traversed by a moving mass

In this section we seek the solution to the entire equation (39) when no term of the equation is neglected. To solve this problem, we use the modified asymptotic method of Struble already alluded to [16]. To this end, we rearrange equation (39) to take the form

where Γ has been written as a function of the mass ratio η₀.

Thus, considering the homogeneous part of the equation (58) and going through the same arguments and analysis as the previous case, the modified frequency corresponding to the frequency of the free system due to the presence of the moving mass is

retaining terms to o(η₀) only.

Therefore, to solve the non-homogeneous equation (58), the differential operator which acts on T_n(t) and T_q(t) is replaced by the equivalent free system operator defined by the modified frequency β_f. That is

where

Clearly, equation (60) is directly analogous to equation (56). Hence when equation (60) is solved in conjunction with the initial conditions, one obtains expression for T_n(t). Thus in view of equation (10), we have

Equation (62) represents the transverse-displacement response to a moving mass of a simply supported rectangular plate on a variable Winkler elastic foundation.

4 DISCUSSION OF THE ANALYTICAL SOLUTIONS

In studying undamped system such as this, it is desirable to examine the phenomenon of resonance. Equation (57) clearly shows that the simply supported rectangular plate on a variable Winkler elastic foundation and traversed by a moving force reaches a state of resonance whenever

while equation (62) shows that the same plate under the action of a moving mass experiences resonance when

where

Equations (64) and (65) imply that

Since < 1 for all n_i, it can be deduced from equation (66) that, for the same natural frequency, the critical speed (and the natural frequency) for the system of a simply supported rectangular plate traversed by a moving mass is smaller than that of the same system traversed by a moving force. Thus, for the same natural frequency of the plate, resonance is reached earlier when we consider the moving mass system than when we consider the moving force system.

5 NUMERICAL CALCULATIONS AND DISCUSSION OF RESULTS

In order to carry out the calculations of practical interests in dynamics of structures and engineering design for the simply supported plate resting on variable Winkler elastic foundation, a rectangular plate of length L_Y = 0.914m and breadth L_X = 0.457m is considered. It is assumed that the mass travels at the constant velocity 0.8123m/s. Furthermore, values for E, S and Γ are chosen to be 2.109×10⁹kg/m², 0.4m and 0.2 respectively. For various values of the foundation modulus F₀ and the rotatory inertia correction factor R₀, the deflections of the simply supported plate are calculated and plotted against time t.

Figures 1 and 2 display the effect of foundation modulus (F₀) on the transverse deflection of the simply supported rectangular plate in both cases of moving force and moving mass respectively. The graphs show that the response amplitude decreases as the value of the foundation modulus increases. Values of F₀ between 0 N/m³ and 3000 N/m³ are used.

The effect of rotatory inertia correction factor (R₀) on the transverse deflection in both cases of moving force and moving mass displayed in figures 3 and 4 respectively show that an increase in the value of the rotatory inertia correction factor decreases the deflection of the simply supported rectangular plate resting on variable Winkler elastic foundation. Here, values of R₀ between 0m and 0.2m are used.

Figure 5 compares the displacement curves of the moving force and moving mass for a simply supported rectangular plate with F₀ = 1000 N/m³ and R₀ = 1m. Clearly, the response amplitude of a moving mass is greater than that of a moving force problem. However, this result holds for other choices of the values of F₀ and R₀.

6 CONCLUSION

The problem of the dynamic behaviour under moving concentrated masses of rectangular plates resting on variable elastic foundation is considered in this work. The governing fourth order partial differential equation is a non-homogenous equation with variable and singular coefficients. The objective of the work has been to study the problem of the dynamic response to moving concentrated masses of rectangular plates on variable Winkler elastic foundations. In particular, the closed form solutions of the fourth order partial differential equations with variable and singular coefficients of the rectangular plate is obtained for both cases of moving force and moving mass. The method is based on (i) Separation of variables (ii) The modified Struble's technique and (iii) The method of integral transformations.

These solutions are analyzed and resonance conditions are obtained for the problem. The numerical analysis for both moving force and moving mass problems carried out show that the moving force solution is not an upper bound for the accurate solution of the moving mass solution and that as the rotatory inertia correction factor increases, the response amplitudes of the plates decrease for both cases of moving force and moving mass problem. When the rotatory inertia correction factor is fixed, the displacements of the simply supported rectangular plates resting on variable Winkler elastic foundations decrease as the foundation modulus increases.

Furthermore, for fixed values of rotatory inertia correction factor and foundation modulus, the response amplitude for the moving mass problem is greater than that of the moving force problem implying that resonance is reached earlier in moving mass problem than in moving force problem of the simply supported rectangular plate resting on variable elastic foundation. Hence, it is dangerous to rely on the moving force solutions.

Finally, for the simply supported rectangular plate resting on Winkler elastic foundation with stiffness variation, for the same natural frequency, the critical speed for moving mass problem is smaller than that of the moving force problem, and as rotatory inertia correction factor and the foundation modulus increase, the critical speeds increase showing that risk is reduced.

Received 10 Dec 2010;

In revised form 15 Sep 2011

• [1] G. G. Adams. Critical speeds and the response of a tensioned beam on an elastic foundation to repetitive moving loads. Int. Jour. Mech. Sci., 7:773781, 1995.
• [2] Omer Civalek and Aitung Yauas. Large deflection static analysis of rectangular plates on two parameter elastic foundations. International Journal of Science and Technology, 1(1):4350, 2006.
• [3] J. Clastornic, M. Eisenberger, D. Z. Yankelevsky, and M. A. Adin. Beams on variable elastic foundation. Journal of Applied Mechanics, 53:925928, 1986.
• [4] J. Clastornik, M. Eisenberger, D. Z. Yankelevsky, and M. A. Adin. Beams on variable elastic foundation. Journal of Applied Mechanics, 53:925928, 1986.
• [5] R. Douglas Gregory, Charles C. Gu, and Frederic Y. M. Wan. The cantilever strip plate of varying thickness and the centre of shear. Q. JI Mech. Appl. Math., 55(1):2948, 2002.
• [6] J. N. Franklin and R. F. Scott. Beam equation with variable foundation coefficient. ASCE. Eng. Mech. Div., 105:811827, 1979.
• [7] L. Fryba. Vibration of solids and structures under moving loads Groningen, Noordhoff, 1972.
• [8] J. A. Gbadeyan and S. T. Oni. Dynamic behaviour of beams and rectangular plates under moving loads. Journal of Sound and Vibration, 182(5):677693, 1995.
• [9] H. P. Lee and T. Y. Ng. Transverse vibration of a plate moving over multiple point supports. Applied Acoustics, 47(4):291301, 1996.
• [10] M. Lentini. Numerical solution of the beam equation with non-uniform foundation coefficient. ASME. Journal of Applied Mechanics, 46:901902, 1979.
• [11] M. Milornir, M. Stanisic, and J. C. Hardin. On the response of beams to an arbitrary number of concentrated moving masses. Journal of the Franklin Institute, 287(2), 1969.
• [12] G. Muscolino and A. Palmeri. Response of beams resting on viscoelastically damped foundation to moving oscillators. International Journal of Solids and Structures, 44(5):13171336, 2007.
• [13] S. T. Oni. Flexural motions of a uniform beam under the actions of a concentrated mass traveling with variable velocity. Abacus, Journal of Mathematical Association of Nigeria, 31(2a):7993, 2004.
• [14] S. T. Oni and T. O. Awodola. Vibrations under a moving load of a non-uniform rayleigh beam on variable elastic foundation. Journal of Nigerian Association of Mathematical Physics, 7:191206, 2003.
• [15] S. T. Oni and T. O. Awodola. Dynamic response to moving concentrated masses of uniform rayleigh beams resting on variable winkler elastic foundation. Journal of the Nigerian Association of Mathematical Physics, 9:151162, 2005.
• [16] S. T. Oni and T. O. Awodola. Dynamic response under a moving load of an elastically supported non-prismatic bernoulli-euler beam on variable elastic foundation. Latin American Journal of Solids and Structures, 7:320, 2010.
• [17] S. T. Oni and B. Omolofe. Dynamic analysis of a prestressed elastic beam with general boundary conditions under moving loads at varying velocities. Journal of Engineering and Engineering Technology, FUTA, 4(1):5572, 2005.
• [18] A. V. Pesterev, C. A. Tan, and L. A. Bergman. A new method for calculating bending moment and shear force in moving load problems. Journal of Applied mechanics, 68:252259, 2001.
• [19] S. Sadiku and H. H. E. Leipholz. On the dynamics of elastic systems with moving concentrated masses. Ing. Archiv., 57:223242, 1981.
• [20] E. Savin. Dynamics amplification factor and response spectrum for the evaluation of vibrations of beams under successive moving loads. Journal of Sound and Vibrations, 248(2):267288, 2001.
• [21] M. R. Shadnam, M. Mofid, and J. E. Akin. On the dynamic response of rectangular plate, with moving mass. Thin-Walled Structures, 39:797806, 2001.
• [22] M.M. Stanisic, J.A. Euler, and S.T. Montgomery. On a theory concerning the dynamical behaviour of structures carrying moving masses. Ing. Archiv., 43:295303, 1974.
• [23] Jia-Jang Wu. Vibration analysis of a portal frame under the action of a moving distributed mass using moving mass element. Int. Jour. for Numerical Methods in Engineering, 62:20282052, 2006.

Publication Dates

• Publication in this collection
  17 Apr 2012
• Date of issue
  2011

History

• Accepted
  15 Sept 2011
• Received
  10 Dec 2010