Research Article  Open Access
Periodic Responses of a Rotating HubBeam System with a Tip Mass under Gravity Loads by the Incremental Harmonic Balance Method
Abstract
Dynamic characteristics of a flexible hubbeam system with a tip mass under gravity loads are investigated. The slope angle of the centroid line of the beam is utilized to describe its motion. Hamiltonâ€™s principle is used to derive the equations of motion and their boundary conditions. By using Lagrangeâ€™s equations, spatially discretized equations based on assumed mode method are derived, and the equations of motion are expressed in nondimensional matrix form. The incremental harmonic balance (IHB) method is used to solve for periodic responses of a highdimensional model of the rotating hubbeam system with a tip mass for which convergence is reached. A frequency equation is derived giving the relationship between the nondimensional natural frequencies and three nondimensional parameters, that is, the rotating angular velocity, the tip mass, and the hub radius ratio. A comparative study is performed for nonlinear frequency responses of the system with a tip mass under different values of tip masses and damping ratios.
1. Introduction
A number of systems in the fields of navigation and mechanical engineering can be modeled as rotating hubbeam system with a tip mass. For example, flexible manipulators, spacecraft structures, and cranes carrying moving loads can be studied in this way. In these systems, the beams carrying a tip mass are rotating in a horizontal plane with the whole systems mounted on a rotating hub. In order to study their dynamic characteristics, the rotating beam systems are simplified as a rotating hubbeam system with a tip mass model. In recent years much attention has been placed on linear and nonlinear dynamic characteristics of rotating beam with a tip mass, and earlier work has been done with the linear analysis. For instance, Conrad and MorgÃ¼l [1] used linearized feedback law to study the stabilization of a flexible beam with a tip mass. Rao [2] derived the timedependent equations of motion that governs the vibration of an EulerBernoulli beam; it was used to obtain the linear dynamic response of the beam under moving load mass. Demetriou [3] presented a method for construction of observer for linear secondorder lumped and distributed parameter systems using parameterdependent Lyapunov functions. Furta [4] proved that the attached point mass on a thin elastic beam plays a destabilizing role for any values of the problem parameters and studied the dynamical stability of the rectilinear shape of the beam by means of the direct Lyapunov method. In [5], the extended Hamilton principle was employed to derive the equations of motion of a rotating beam with a tip mass undergoing coupled torsionalbending vibrations and analyzing the exact frequencies leading to a better control of the system.
Recently, the nonlinear vibration of a rotating beam with a tip mass has been studied by numerous researchers. Yang et al. [6] presented a finite element model for a flexible hubbeam system with a tip mass where viscous damping of the hub and the air drag force were considered. They showed that the traditional linear model cannot account for the dynamic stiffening and it may lead to erroneous results in highspeed systems. Simulation results on the vibration measurements of the rotating flexible beam carrying a tip mass for different combinations of shaft and root flexibilities and arm properties were reported by Ismail et al. [7]. Gregory et al. [8] have pointed out the shortcomings of the sequential singleaxis vibration method. They have also reported evidence of differences in failure modes and fatigue life for multiaxial loadings versus singleaxis inputs by utilizing multiaxial electrodynamic shakers. Sunar and AlBedoor [9] tested the suitability of a piezoelectric (PZT) sensor in measuring vibrations of the rotating beams. A comparison between the results of finite element and experimental was reported and indicated that the rootembedded PZT sensor can be effectively used in measurements of blade vibration. In [10], Sinha et al. investigated a simultaneously precessing and nutating beam with a tip mass using a variant of Hillâ€™s method and found that the stability of an only precessing beam depends on the inclination of the beamcenterline with the axis of precession. Patil and Gandhi [11] demonstrated chaotic behavior for an inverted flexible pendulum with a tip mass on a cart system; the dynamic model was validated with experimental data for a couple of cases of beam excitation. Li et al. [12] established the dynamic equations of the rigidflexible coupling system in a noninertial coordinate system; the numerical simulation results showed that dynamic stiffening is produced by the coupling effect of the centrifugal inertial load distributed on the flexible beam and the transverse vibration deformation of the beam. Yong et al. [13] studied the rigidflexible coupling system with a hub and concentrated mass; they utilized the second Lagrange equation and the assumed mode method to establish the dynamic equations of the system. The results showed that the concentrated mass mainly suppresses the bending vibration of the beam and exhibits damping characteristics.
In a rotating hubbeam system with a tip mass, the effect of the gravity load of beam was not accounted for in the previous literature [14]. As the size and height of rotating machineries such as manipulators and cranes increase, it is necessary to consider the effect of the gravity load of the beam during the rotation. Park and Kim [15] derived the equations of motion of the rotating beam which include all dynamic effects and mainly studied the effects of curvature and tip mass that can change the dynamic response of the beam. Flatness based controller design techniques for a rotating hubbeam system with a payload attached to the tip of the beam in a gravitational field were utilized in [16]. Cai et al. [17] made a comparison between the firstorder approximation coupling model and the zerothorder approximation coupling model when prescribed torque drives the beam with a tip mass, which takes into account the influence of the gravitational force filed of the beam in the analysis of the dynamic characteristics of the system. As for a flexible hub geometrical nonlinearity beam with a tip mass, Emam [18] employed a flexural model to study the dynamic responses of a flexible hub geometrical beam with a tip mass of which the hub is restrained by a translational and a rotational spring; such a model accounts for the geometrical coupling between the axial and lateral deformations. If gravity of the beam is much less than the gravity of tip mass, it is not necessary to consider the effect of the gravity load of beam [19].
The main objective of this work is to use the IHB method to analyze nonlinear dynamic responses of a rotating hubbeam system with a tip mass for which convergence is reached. The vibration analysis is performed under the slope angle model with five included trial functions, and the model considers the effects of gravity loads of the beam and tip mass during the rotation. To deal with the rotating hubbeam system with a tip mass, the IHB method is applied, which was first introduced by Lau and Cheung [20, 21]. The IHB method was successfully applied to the analysis of nonlinear structural vibrations problems. Xu and Zhu [22] used the IHB method to determine parametric instability boundaries of a parametrically excited system. In order to improve calculation efficiency of the IHB method, the fast Fourier transform (FFT) was used to efficiently calculate coefficients of a discrete Fourier transform in practice [23]. Wang and Zhu [24] developed a modified IHB method that incorporated FFT and Broydenâ€™s method to obtain the dynamic responses of beltdrive system with a noncircular sprocket. The configuration of the rotating beam is described by the slope angle of its centroid line; the slope angle model has been widely used in studying the nonlinear dynamics [25]. The nonlinear dynamic equation of the rotating hubbeam system with a tip mass is derived in Section 2 using Hamiltonâ€™s principle. Spatially discretized equations of the slope angle model neglecting the fourth and higherorder nonlinear terms are derived by Lagrangeâ€™s equations in Section 3. The IHB method for calculating the periodic responses of the rotating hubbeam system with a tip mass is formulated in Section 4. The critical buckling load of the rotating beam under an axial tip mass is calculated in Section 5. Natural frequencies of the rotating beam with a constant rotating angular velocity are calculated in Section 6.1 using the slope angle model, and compared with those from the dynamical model in [14]. Nonlinear frequency responses of the system with a tip mass for different numbers of trial functions are calculated in Section 6.2. The influences of tip mass and damping ratio on nonlinear frequency responses are studied in Sections 6.3 and 6.4, respectively.
2. Equations of Motion
A schematic diagram of a rotating beam with a tip mass is shown in Figure 1, a rotating planar beam is attached to a rigid hub of radius , and the beam rotates in a vertical plane around the central axis in the inertial coordinates . The effects of gravity loads of the beam and tip mass are considered. The beam has a slender shape that the shear deformation effect and rotary inertia of the beam are ignored. The position of a cross section of the beam is described by the arclength coordinate along its centroid line. For the whole beam, the bending stiffness , crosssection area , and linear mass density can be approximated as constant. is the length of the beam, and is the tip mass at the end of the beam. The hub rotates with a given rotating angular velocity , where is time; the gravitational acceleration is . represents the space attitude motion of hub, andThe slope angle of the beam at the position is expressed as , then .
The Cartesian coordinates of any one particle of the beam with the arclength coordinate areThe Cartesian coordinates are functions of ; hence, one can use to describe the configuration of the beam. The kinetic energy of the rotating hubbeam system with a tip mass iswhere a dot over and denotes differentiation of them with respect to and is the rotary inertia of the hub.
The potential energy of the system with a tip mass can be expressed aswhere and represent the elastic potential energy and the gravitational potential energy, respectively. For brevity in expression, the variables and are written as and . Substituting (2) into (3) and (4), one obtains the variations of the kinetic energy and the potential energywhere denotes partial differentiation of with respect to . The boundary term , which gives the following boundary conditions:The relationship between and can be expressed as follows:where is the slope angle of the centroid line of the beam measured from a rotating axis attached to the hub. The boundary conditions in (6) can then becomewhere and denote partial differentiation of and with respect to , respectively. Using Hamiltonâ€™s principle and (7), the governing integropartial differential equation is obtained as follows:
Equation (9) presents a new rigidflexible dynamic model of a rotating hubbeam system with a tip mass described by rotating angle of the hub and slope angle .
3. Spatial Discretization
Assume the solution of (9) can be expressed bywhere the superscript denotes transpose of a vector or matrix, are trial functions used, which are eigenvalue functions of a fixedfree rod, is the number of included trial functions, are the generalized coordinates, , and .
The position vector of the particle and tip mass in the inertial coordinate system can be expressed by , , and , and their Cartesian coordinates corresponding to arclength coordinate are as follows, respectively:where is a vector tangent to the centroid line of the beam at the base. By differentiating (11) and (12) with respect to , one obtains the velocity of the particle and tip masswherein which
The trigonometric functions and in (11) and (12) are approximated using Taylor expansion around ; one obtainswhere
The tip mass was considered to be a particle so that the moment of inertia of tip mass is set to be zero. The kinetic energy of the beam and tip mass can be approximated by neglecting the nonlinear terms over fourth orderThe elastic potential energy of the beam can be written asGravity loads on the beam and tip mass areThe virtual displacements of the particle on the beam and tip mass areThe virtual work done by on the beam and the virtual work done by on the tip mass areand generalized forces associated with gravity loads on the beam and tip mass areAs a matter of convenience, for a rotating beam with a constant rotating angular velocity, introduce the following dimensionless variables:Substitute (19), (20), (24), and (25) into Lagrangeâ€™s equationsand dropping the over bar of each variable for brevity, one can obtain the following dimensionless equations of motion:where , , , and are constant coefficients for the linear terms; and are constant coefficients for the quadratic geometrical nonlinear terms; , , and are constant coefficients for the cubic geometrical nonlinear terms; are amplitudes of the harmonic forcing terms; and are constant coefficient of a damping term that has been introduced via a structural damping ratio . All the coefficients in (27) can be seen in Appendix A. The trial functions in terms of the dimensionless variable are
4. IHB Method for Determining Periodic Responses
Among the solution methods of ordinary differential equations, the IHB method is relatively efficient and reliable for treating the vibrations of nonlinear systems. The method has been developed and successfully applied to analyze periodic nonlinear structural vibrations. With the new dimensionless time variable defined aswhere is the unknown fundamental frequency of the steadystate response, and this paper assumes that the rotating angular velocity can be , in which is the positive integer introduced to analyze possible subharmonic responses of the nonlinear system. Equation (27) becomeswhere and respect the first and second differentiation of with respect to , respectively.
The IHB method consists of two steps; the first step of the IHB method is an incremental process, which is a NewtonRaphson procedure. Let and denote the solutions of the equation, where ; a neighboring state can be expressed by adding increments to and
Substituting (31) into (30) and neglecting higherorder terms of the increments yield the following linearized incremental equation with and as the unknown quantity:whereis a correction term which goes to zero when the solution is exact. Equations (32) and (33) can be written in the matrix formwhereis the correction vector, , , , , , , , and have entries , , , , , , , and , respectively; and
The second step of the IHB method is a harmonic balance procedure; one assumes and can be expressed by expanding them in truncated Fourier serieswhere is the transformation matrix, it consists of trigonometric functions, is the harmonic term coefficient matrix. , , and can be written asin which and in (37) are Fourier coefficients and and are numbers of cosine and sine harmonic terms reserved, respectively.
Hence, vectors of generalized coordinates and their increments can be written aswhere , , and ; diag represents extracting only diagonal items in matrix .
Differentiating (40) yieldsUsing Galerkin procedure to balance harmonics in (34) yieldsSubstituting (40) and (41) into (42) yields a series linear equations with and as the unknown quantitieswherein whichThe vectors and entries of matrices in (45) are given in detail in Appendix B.
Note that the number of incremental unknowns ( and ) in (43) is one more than the number of equations; the solution process starts from a guessed solution. The nonlinear amplitudefrequency response is solved point by point by adding the small increments and to the current solution. The NewtonRaphson iterative procedure is repeated until the convergent solution is reached. One can also choose the arclength increment method combined with a cubic extrapolation technique in [26] to trace the solutions around some peaks in the nonlinear frequency response curve.
5. Analysis of Buckling of the Beam
Buckling and large deformation are classic problems in flexible link systems, especially those in vertical planes. For exploring the critical buckling load, the rotating beam with a tip mass is setting in vertical position and rotating angular velocity , as shown in Figure 2. It is well known that the critical buckling load for a fixedfree cantilever beam can be calculated as follows:
The properties of the rotating hubbeam system with a tip mass are given with the length ; the hub radius ; the mass density ; the modulus of elasticity ; the crosssection area and the area rotary inertia ; and the dimensionless parameter . Substituting into (25) yields the dimensionless tip massThis indicates that the buckling of the beam will not happen to the rotating system with the dimensionless tip mass .
6. Results and Discussion
6.1. Natural Frequencies of the Rotating HubBeam System with a Tip Mass
To investigate the bending vibration natural frequencies of the rotating beam at various rotational velocities under different values of tip masses and hub radius ratios, natural frequencies of the rotating beam can be approximated by neglecting the nonlinear terms, the damping terms, and gravity loads in (27)where and are the static stiffness term and the dynamic stiffness term relating to rotation of the hub, respectively. To derive the eigenvalue problem associated with (48), one assumeswhere is an imaginary number, stands for a bending natural frequency, and is a constant column matrix characterizing the mode shape. Substituting (49) into (48) yieldswhere and have entries and , respectively. The dimensionless natural frequency is assessed by calculating the eigenvalues of (50).
Results of the first three dimensionless natural frequencies of the rotating hubbeam system with a tip mass from the proposed model are compared with those from [14] in Table 1. In the case of low velocity rotation, there is no need to consider the stretching of beam. And five modes are considered to reach convergence. As can be noted from Table 1, there are only small differences between the results from the two models. This indicates the use of slope angle model with five modes is sufficient to obtain right results.

It can be seen from Figure 3 that when hub radius ratio , an increase in rotating angular velocity causes an increase in the first four natural frequencies for the five modes of vibration under different values of tip masses. Figure 3(a) shows that an increase in the tip mass always causes a decrease in the first natural frequency for all rotating angular velocities. In Figures 3(b)â€“3(d), the tip mass tends to depress the frequencies at low velocities of rotation; however, it tends to increase the frequencies at high velocities of rotation. Similar consequence was also obtained in [27].
(a)
(b)
(c)
(d)
Figure 4 shows the first four dimensionless natural frequencies veering behavior of the rotating beam associated with different hub radius ratios. It can be seen that the natural frequencies increase with the hub radius ratio and the rotating angular velocity.
(a)
(b)
(c)
(d)
6.2. Nonlinear Frequency Response of the Rotating HubBeam System with a Tip Mass under Gravity Loads
Consider the damped vibration of the rotating hubbeam system with a tip mass with gravity loads. The properties of the beam are chosen the same as in Section 5. With and , (37) can be developed aswhere
The nonlinear frequency responses of the system with a tip mass and when using different numbers of included trial functions are shown in Figures 5(a) and 5(b), respectively, where and are the amplitudes of the first harmonic terms of the first two generalized coordinates and , respectively, and is the dimensionless excitation frequency. Subharmonic responses corresponding to are not found in the frequency range shown with the numbers of included trial functions. The higher harmonics are small and not shown here. Comparison of various included trial functions for nonlinear frequency responses indicates that convergence is reached when . Therefore, it verifies that five modes used in this work are sufficient to obtain convergent results. Due to the quadratic nonlinearities, the response of the system with a tip mass exhibits a typical softeningspring nonlinear characteristic.
(a)
(b)
6.3. FiveMode Model for Bending Vibration Analysis with Consideration of Variation in Tip Mass
Consider the damped vibration of the rotating hubbeam system with a tip mass with respect to the variation of tip mass: , , and . In Section 6.2, it is indicated that a fivemode model is efficient and accurate to calculate nonlinear frequency response. In all cases below, , , and in (37) are used. The tip deformation of the rotating beam is given by
Figure 6 shows the effect of tip mass on tip amplitude of the rotating beam , where the tip amplitude of the beam is defined as . It can be seen that the tip amplitude of the rotating flexible beam gets larger with increasing the tip mass. Moreover, it can be seen that the region which resonant behavior appears tends to be located in the lower frequency range as the tip mass increases. Resonant behavior can be produced when the system is excited at a frequency close to the natural frequency. The systemâ€™s natural frequency strongly depends on the tip mass.
The comparison of bending vibration analysis between the gravity of beam considered and ignored model is also presented in Figure 6. It can be seen that discrepancies of results from the two models are remarkable, and tip amplitude calculated from the gravity of beam ignored model is much smaller than that from the gravity of beam considered model. This indicates that consideration of taking into account the gravity of beam in the model of rotating hubbeam system with a tip mass is necessary.
6.4. FiveMode Model for Bending Vibration Analysis with Consideration of Variation in Damping Ratio
Figure 7 shows nonlinear frequency response curves of the rotating flexible beam for different values of damping ratios with , , and . As shown in Figure 7, the tip amplitude of the rotating flexible beam decreases as damping ratio increases. It indicates that damping ratio plays an important role to decrease the tip amplitude and enhance the stabilizing effect on the system. In addition, the tip amplitudes of the rotating flexible beam have the same region which resonant behavior appears for all of damping ratios, and the tip amplitudes are small and almost identical when excitation frequency is far separated from the natural frequency.
7. Conclusions
The dynamic characteristics of a flexible hubbeam system with a tip mass have been investigated in this work. The motion of the beam with a tip mass is described based on the slope angle model, and gravity loads of beam and tip mass introduce combined parametric and forced harmonic excitations in the resulting equations of motion. The IHB method is used to solve for periodic responses of a fivedimensional model of the rotating hubbeam system with a tip mass for which convergence is reached.
The buckling will not happen to the rotating beam if the value of tip mass within the critical buckling load range. Natural frequencies of the rotating hubbeam system with a tip mass calculated using the slope angle model with five assumed modes show good agreement with those in [14]. Comparisons of the nonlinear frequency response curves under different values of tip masses and damping ratios are obtained by using the IHB method. The nonlinear frequency response curves show that the tip amplitude of the rotating beam increases with the tip mass, while the tip amplitude of that decreases with damping ratio. On the other hand, the frequency response curves show different structures of resonance regions according to the values of the tip mass and damping ratio. In addition, large discrepancy exists between the results obtained from the gravity of beam considered and those from the gravity of beam ignored.
Appendix
A. Coefficients in the Spatially Discretized Equations in (27)
The coefficients in (27) about the linear, quadratic nonlinear, cubic nonlinear, damping, and excitation force terms are written as