Vibrations of beams with a variable cross-section fixed on rotational rigid disks

Vibrations of beams with a variable cross-section fixed on rotational rigid disks

Slawomir Zolkiewski^*

Silesian University of Technology, Institute of Engineering Processes Automation and Integrated Manufacturing Systems, Division of Mechatronics and Designing of Technical Systems, Faculty of Mechanical Engineering, 18a Konarskiego Street, 44-100 Gliwice, Poland

ABSTRACT

The work is focused on the problem of vibrating beams with a variable cross-section fixed on a rotational rigid disk. The beam is loaded by a transversal time varying force orthogonal to an axis of the beam and simultaneously parallel to the disk's plane. There are many ways of usage of the technical moveable systems composed of elements with the variable cross-sections. The main applications are used in numerous types of turbines and pumps. The paper is a kind of introduction to the dynamic analysis of above mentioned beam systems. The equations of motion of rotational beams fixed on the rigid disks were derived. After introducing the Coriolis forces and the centrifugal forces, the transportation effect in the mathematical model was considered. This particular project is the first stage research, where there were proposed certain solutions of problems connected with the linear variable cross-sections systems. The further investigation considering the nonlinear systems has been proceeding. The results, analysis and comparison will be presented in the future works.

Keywords: variable cross-section, beam, rotary motion, dynamic flexibility

1 INTRODUCTION

The fundamental elements of the numerous mechanisms or machines, namely the beams and rods, are most commonly considered as the vibrating stationary systems. In these cases they are treated as the immovable ones without the coupling between the local oscillations and the main operational motion. The mechanical systems are predominantly considered both in the kinematical and dynamical aspects. The main problems resulting from the analysis are connected with the issues of controlling and stabilizing mechanical systems as such. Numerous research can be found in the literature [1, 3, 4, 7, 9, 14-16, 18-22] and so far the solutions have been found by considering the transportation (considered as the main working motion) and the local oscillations separately [2, 3, 5, 6, 8, 10-12, 17]. Mesut Simsek in his paper [14] presents the nonlinear dynamic analysis of a functionally graded beam with pinned-pinned supports. Simsek uses the Timoshenko beam theory with the von Karman's nonlinear strain-displacement relationships. The beam was analyzed in relation to its variable thickness of a cross-section according to the power-law formula. The system of equations of motion is derived by using the Lagrange's equations. Transverse and axial deflections and rotation of the cross-sections of the beam were expressed in polynomial forms. The boundary conditions are taken into account by using the Lagrange multipliers. The Newmark-b method in conjunction with the direct iteration method was used to solve nonlinear equations of motion. In the work, large deflection, velocity of the moving load and excitation frequency on the beam displacements were analyzed. Some detailed cases were compared with each other. Another method of analysis of high-speed rotating beams is presented in [9] by Gunda. The author introduces a new finite element for free vibration analysis of rotating beams. The basis shape functions which apply a linear combination of the solution of the governing static differential equation of a stiff-string and a cubic polynomial are used. The introduced functions depend on the angular velocity and position of the analyzed section, what accounts for the centrifugal stiffening effect. In paper [1] the effects of rotary inertia were presented. The author shows the rotation impact on the extensional tensile force and the eigenvalues of beams. The analyzed beam was rotated uniformly about a transverse axis taking into account longitudinal elasticity. The perturbation technique and Galerkin's method were used. In the example of a typical helicopter rotor blade it was indicated that the extensional tensile force increases up to ten percent, when the rotary inertia contribution is retained in the modelling. Szefer in [15, 16] claimed that the vibrations from flexibility of elements of the mechanical composition are much smaller than a main dislocation of this composition. There is an increased range of velocities and accelerations of such a type systems. All those are caused by using more efficient drives or/and less and less weighting materials. In order to constraint power output of drives for motion of mechanical systems, materials with the lower mass density and stiffness are used. All those things are reasons for creating the new models of designing technical beamlike systems. The variable cross-section beams are applied in many practical applications nowadays. Both the linearly and nonlinearly variable cross-sections are used in technical systems. These systems should be taken under consideration with respect to a realistic distribution of mass and flexibility. Any variations in cross-sections cannot be neglected and should be considered in the mathematical model assumptions. In spite of many works [e.g. 1, 9-12, 17, 22], where their authors analyzes the beams in rotations and beams with variable cross-sections, this work is devoted to general analyses both the linear and non-linear ones of vibrating transversally beams with the variable cross-section. It is believed that new equations and results are presented in this paper.

2 MODEL OF THE ANALYZED BEAMS FIXED ON THE ROTATIONAL RIGID DISK

In this section the way of modelling of beams in rotational transportation is presented. The beam was assumed as the uniform one with the variable cross-section. The beam was fixed on the rotational disk that was assumed as a rigid element (Fig. 1). The disk rotates round its axis and the motion was described in two reference frames, the global stationary reference frame and the local one. The vibrations occurring in the local reference frame are transferred into the global reference frame with taking into consideration so called the transportation effect [3, 16, 19, 20].

where: ρ - mass-density, A(x) - cross-section, l - length of the beam, x - location of the analyzed cross-section, ω - angular velocity, Ω - frequency, Q - rotation matrix, S - position vector, F - harmonic force, E - Young modulus, w - vector of displacement.

There is defined a set of generalized coordinates allowing to make an assumption that generalized coordinates are the orthogonal projections on the coordinate axes of the global reference frame and can be written thus:

and after differentiating (Eq. 1) and (Eq. 2) can be obtained the generalized velocities in the form:

or in the equivalent vector form expressing the displacements:

and the equivalent vector expressing the velocities:

The cross-section is dependent on the x coordinate and for instance making an assumption that the analyzed beam is a cone (a base with radius equal x/l) it can be written as follows:

The angular velocity of the rotational disk treated as the transportation velocity is defined as:

The location of the analyzed cross-section is described by the position vector defined as follows:

The xyz reference frame is the rotational one and the matrix allowing expressing vectors is the rotation matrix. The xyz local reference frame is fixed to the rotating beam and its z-axis coincides with the rotation axis of the rigid disk, thus the rotation matrix is defined as such:

Different applications of such a type beam systems may have different boundary conditions and the cross-section can be considered as the linearly or nonlinearly variable one. For example the beam can be considered as the fixed at one end, the clamped-clamped one or with other boundary conditions etc. For example beams with one end clamped and the second free-free are considered (Fig. 1). The beam is fixed on the rotational rigid disk. An analyzed section of the beam is loaded by a harmonic time varying distributed force with unitary amplitude on a direction perpendicular to the centre line of the beam. The forces on origins of the beams are assumed as equal zero and also the displacements are equal zero, because the system is clamped in one end. Because of all these the specific boundary conditions for this beam should be written as follows:

in every time moment t ≥ 0.

For the beams with the linearly variable cross-section the geometric moment of inertia can be assumed as the function proportional to the constant initial moment I₀ for example as:

After solving the boundary value problem, the eigenfunction for the displacement can be derived in the form:

where [13]:

J_n - the Bessel function of the first kind, Y_n - the Bessel function of the second kind, I_n - the hyperbolic Bessel function of the first kind, K_n - the hyperbolic Bessel function of the second kind, C_1,2,3,4 - the integrations constants, where: n is a mode of vibrations of the beam with a variable cross-section. The k coefficient is a linear approximation of the realistic eigenvalues of the analyzed systems.

3 EQUATIONS OF MOTION OF THE ANALYZED SYSTEMS

The equations of motion of the analyzed beamlike systems are obtained with the help of the classical method. The derived equations of motion can be presented in the matrix form as follows:

in each point of the range the Eq. (3.1) coincides with boundary conditions and initials conditions.

The projections of the equations of motion in the global reference frame are obtained from the Eq. (14) and in the XYZ global reference frame equations for individual axes can be expressed as:

Third governing equation (17) is a well-known one and represents the equation of motion of the immovable beam with the variable cross-section.

As we can see the boundary conditions and the loading direction does not provide any new aspects to the dynamical analysis and the beam (Fig. 2) can be treated as equivalent to the stationary one. Much more interesting case of the loading is the beam loaded by a transversal time varying force orthogonal to an axis of the beam and simultaneously parallel to the disk's plane. This case (Fig. 3) where the beam with a variable cross-section A(x) and a variable moment of inertia I(x) is considered.

Analogically to the previous equations of motion it can be derived by means of the classical method. Interesting case of the loading is the beam loaded by a transversal time varying force orthogonal to an axis of the beam and simultaneously parallel to the disk's plane. The obtained equations of motion are coupled ones with the Coriolis elements. After loading the beam in the direction of the y axis of the local reference frame it can be derived the equations of motion of the beam in the matrix form. This system of equations of motion is the fourth order partial differential equations (PDE) and can be expressed as follows:

in each point of the range the Eq. (18) coincides with boundary conditions and initials conditions.

The projection onto the X and Y axes, can be obtained, where the projection onto the X axis of the global reference frame is as follows:

and for the Y axis of the global reference frame, the projection can be written:

where using the real coordinates the individual displacements onto the axes of the global reference frame are:

It is obvious that the vibrating beams can be described by different boundary conditions for different loading and fixations types and this case can be treated as the general one.

4 MATHEMATICAL MODEL OF THE SYSTEM

The mathematical model of the analyzed beam is provided by the orthogonalization of the equations of motion (19) and (20). The orthogonalization is provided by the Eqs. (19) and (20) which are multiplied using the eigenfunction for the displacements and computing integrals from the equations of motion in the beams limits of integration from the origin of the beam (zero) to the end (a length of the beam). After assuming a constant angular velocity, an angular acceleration equals zero and after introducing Eqs. (21), then the projecting equations of motion Eq. (18) onto the axes of the global reference frame can be written:

and for the second equation, also multiplied by the eigenfunction for displacement and integrated along the beam:

The presented mathematical model is universal for different functions describing the variability of the beam's cross-section. The cross-section of the analyzed beam can be now treated as both linearly and nonlinearly variable. If the beam has the nonlinearly variable cross-section (Fig. 4-5) the appropriate function describing this variation ought to be adopted. The different loading terms can be considered as well. For example the beam is loaded by a transversal force orthogonal to an axis of the beam and orthogonal to the disk's plane (Fig. 4) or the beam is loaded by a transversal force orthogonal to an axis of the beam and in parallel to the disk's plane (Fig. 5).

These cases can be considered separately or combined. According to the scientific literature the case presented in figure 4 is well-known in the mathematical sense, likewise the case in figure 1.

It is assumed that elements of the Eqs. (22) and (23) are the continuously differentiable functions. Using the Eq. (22) and multiple integrating both sides of this equation by parts, the following formula can be obtained:

where the variable cross-section A(x) and the moment of inertia I(x) depend on the geometry of the analyzed beam. There is a beam with the linearly variable cross-section Eqs. (7) and (12) assumed for further calculations. The eigenfunction associated with transversal vibrations X(x) of the beam is obtained by applying undamped conditions. Similarly to the first equation, the integration results for the second equation (23) can be applied as follows:

It is possible to simplify the obtained equations as shown in the Eq. (26). Additionally, after organizing the elements in order, the Eqs. (24) and (25) can be expressed as:

and the coupled second equation:

There are introduced the equivalent notation for the eigenfunctions, cross-sections and the moments of inertia as:

The Eqs. (26) and (27) after substituting (28) are simplified to:

and analogically to the second equation (27) for the Y axis of the global reference frame can be obtained as:

The Eqs. (29) and (30) can be simplified as well and the following two equations yield:

the second equation (30) additionally after simplifying the integrals:

By introducing the acting external force, it can be noticed:

Based on the boundary conditions (11) it can be written:

and

By introducing (34-35) the first equation yields:

Where the searched solutions are defined as displacements in the global reference frame and the ones are separable in time and space domains. The system is symmetric and the motion of the analyzed cross-section is given by the multiplication of the projected harmonic time dependent eigenfunctions with the frequency Ω and the eigenfunctions for displacements with amplitudes XY. Both the variables coincide with the responses of two decoupled harmonic oscillators. In the present case the searched solutions are:

The norm for Eqs. (36) is assumed as:

for the beam with the linearly variable cross-section (Fig. 3) the norm is calculated as integral.

Both I₁ and I₂ coefficients are connected with the moments of inertia and the eigenfunctions for displacements. The I and c (velocity) coefficients are as follows:

It can be also noticed that:

The a coefficients (amplitudes) can be obtained as:

Knowing these amplitudes the total displacements can be calculated and the solution in the global reference frame can be obtained.

5 DYNAMICAL FLEXIBILITIES

According to the author the dynamical flexibility is understood as the amplitude of generalized displacement in the direction of "i" generalized coordinate changed by generalized harmonic force with amplitude equals 1 in the direction of "k" generalized coordinate. This definition can be written in the mathematical form as:

Finally the dynamical flexibility of the rotating beam with the variable cross-section is presented as:

For comparison, the dynamical flexibility of the rotating beam with a constant cross-section is as follows:

On the other side the dynamical flexibility of the immovable beam is presented as:

where for the Eqs. (44) and (45):

and the norm is:

Both the derived dynamical flexibilities, the Eq. (44) and Eq. (46), are presented in figures (Fig. 6-7). For numerical examples the beam of 1m length is assumed. The beam is the clamped-free one and fixed on the rotational disk. Distinctively from the beam with a constant cross-section, the beam with the variable cross-section has different inconstant bifurcations for the individual modes of vibrations (Fig. 7).

Derivations of suitable dynamical flexibilities are applied to the systems in the transportable motion both vibrating longitudinally and transversally. The solutions can be presented on charts of amplitude (dynamical flexibility) in function of frequency. There are many numerical applications that can be used for this purpose, e.g. GENTA'S DynRot or author's program the Modyfit [18].

CONCLUSIONS

The paper concerns a very crucial problem of beams vibrating in transportation with the variable geometric parameters such as: tapered or coned beams, blades with various shape functions, and others. The analyzed systems could be applied in many technical implementations, obviously after proper adjustment of mathematical models and considering the necessary parameters and conditions. The most popular applications could be connected with pumps and rotors; especially the high speed ones, blades of helicopters and arms of robots. The obtained solutions were also presented in the form of equations of motion. The dynamical analysis was provided using the dynamical flexibility method and the Galerkin's method. The dynamical flexibility method applies solutions presentation as the dynamical characteristics of amplitude in the function of frequency. The solutions were determined as the compact mathematical formulae. The dynamical characteristics were presented on charts. The transportation effect and shape functions describing the cross-section variation were taken into consideration. The centrifugal forces and the Coriolis forces were taken into consideration in the obtained mathematical model. The natural frequencies and the minimal amplitudes can be easily read from the dynamical flexibility characteristics. In these characteristics the beam with the constant cross-section and the beam with the variable cross-section are put together. The juxtaposition of the rotating beam and the immovable one is presented. The results presented are very significant and enable controlling of such a type rotating systems. In future works more numerical examples for beamlike systems with the variable cross-section will be presented and the rotating beam with damping forces considered will be analyzed. The rotating beam with some imperfections and defects or cracks can be the potential material for the further analysis and research.

Acknowledgments: the paper was presented during 11th Conference on Dynamical Systems - Theory and Applications.

• Al-Ansary M. D.: Flexural vibrations of rotating beams considering rotary inertia. Computers and Structures, 69, 1998, 321-328.
• Auciello N. M.: Comment on. A note on vibrating tapered beams. Journal of Sound and Vibration, 1995, 187-724.
• Bokhonsky A. I., Zolkiewski S.: Modelling and Analysis of Elastic Systems in Motion. Monograph, Gliwice, Silesian University of Technology Press, 2011, p. 171.
• Chang T.-P., Chang H.-C.: Vibration and buckling analysis of rectangular plates with nonlinear elastic end restraints against rotation. International Journal of Solids and Structures, 34, 18, June 1997, 2291-2301.
• Cheung Y. K., Zhou D.: The free vibrations of tapered rectangular plates using a new set of beam functions with the Rayleigh - Ritz method. Journal of Sound and Vibration, 223, 5, 1999, 703-722.
• Dyniewicz B., Bajer C. I.: New feature of the solution of a Timoshenko beam carrying the moving mass particle. Archives of Mechanics, 62, 5, 2010, 327-341.
• Genta G.: Dynamics of Rotating Systems. New York, Springer, 2005.
• Grossi RO, Bhat RB. A note on vibrating tapered beams. Journal of Sound and Vibration 1991, 147-174.
• Gunda J. B., Gupta R. K., Gangul R.: Hybrid stiff-string - polynomial basis functions for vibration analysis of high speed rotating beams. Computers and Structures, 87, 2009, 254-265.
• Kumar R., Kansal T.: Effect of rotation on Rayleigh waves in an isotropic generalized thermoelastic diffusive half-space. Archives of Mechanics, 60, 5, 2008, 421-443.
• Ozturk H.: In-plane free vibration of a pre-stressed curved beam obtained from a large deflected cantilever beam. Finite Elements in Analysis and Design, 47, 2011, 229-236.
• Rahman M. A., Kowser M. A., Hossain S. M. M.: Large deflection of the cantilever steel beams of uniform strength - experiment and nonlinear analysis. International Journal of Theoretical and Applied Mechanics, 1, 2006, 21-36.
• Solecki R., Szymkiewicz J.: The rod and surface systems. Warszawa, Arkady Budownictwo, Sztuka, Architektura, 1964 (in Polish).
• Simsek M.: Non-linear vibration analysis of a functionally graded Timoshenko beam under action of a moving harmonic load. Composite Structures, 92, 2010, 2532-2546.
• Szefer G.: Dynamics of elastic bodies undergoing large motions and unilateral contact. Journal of Technical Physics. XLI, 4, 2000.
• Szefer G.: Dynamics of elastic bodies in terms of plane frictional motion. Journal of Theoretical and Applied Mechanics, 2, 39, 2001.
• Zhou D.: Vibrations of point-supported rectangular plates with variable thickness using a set of static tapered beam functions. International Journal of Mechanical Sciences, 44, 2002, 149-164.
• Zolkiewski S.: Numerical Application for Dynamical Analysis of Rod and Beam Systems in Transportation. Solid State Phenomena, 164, 2010, 343-348.
• Zolkiewski S.: Dynamic flexibility of the supported-clamped beam in transportation. Journal of Vibroengineering. 13, 4, December 2011, 810-816.
• Zolkiewski S.: Dynamical Flexibility of Complex Damped Systems Vibrating Transversally in Transportation. Solid State Phenomena, 164, 2010 339-342.
• Zolkiewski S.: Damped Vibrations Problem Of Beams Fixed On The Rotational Disk. International Journal of Bifurcation and Chaos, 21, 10, October 2011, 3033-3041.
• Zone-Ching L., Don-Tsun L.: Dynamic deflection analysis of a planar robot Computers and Structures, 53, 4, 17 November 1994, 947-960.

Publication Dates

• Publication in this collection
  25 Feb 2013
• Date of issue
  Jan 2013