On dynamics of a rotating hubflexible rodconcentrated mass system considering rigidflexible coupling effect
Jianshi Fang^{1} , Jianqiang Xiao^{2} , Liang Li^{3} , Dingguo Zhang^{4}
^{1}School of Materials Science and Engineering, Nanjing Institute of Technology, Nanjing, 211167, P. R. China
^{1}Jiangsu Key Laboratory of Advanced Structural Materials and Application Technology, Nanjing, 211167, P. R. China
^{2}School of Architectural Engineering, Nanjing Institute of Technology, Nanjing, 211167, P. R. China
^{3, 4}School of Science, Nanjing University of Science and Technology, Nanjing, 210094, P. R. China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 22, Issue 5, 2020, p. 10061021.
https://doi.org/10.21595/jve.2020.21088
Received 9 October 2019; received in revised form 30 March 2020; accepted 14 April 2020; published 15 August 2020
JVE Conferences
Rotating hubflexible rod system is a typical rigidflexible coupling dynamic mechanism, which has a wide range of industrial applications. In this paper, a comprehensive nonlinear dynamical model of a rotating hubflexible rodconcentrated mass system considering rigidflexible coupling effect is established to study its dynamic properties. By employing the Hamilton principle and classical beam theory, a set of differential equations of motion are derived including the couplings of the elastic deformation of the rod and the rigid rotation of the hub. The additional centrifugal force, tangential force and Coriolis force due to the rigidflexible coupling effect are elaborated. The derived governing partial differential equations are solved by the Galerkin method. The validity of the present model is verified by a comparative study. The tip motion trajectories of the rod for the prescribed rotation, the dynamic responses of the hub and the rod for an external torque acting on the hub and the dimensionless natural frequencies of the system for the steadystate rotation are graphically presented. The influences of parameters such as rotational speed ratio, concentrated mass ratio, concentrated mass location ratio and initial eccentricity ratio on the dynamics are discussed in detail.
 A dynamic model of a rotating hubflexible rodconcentrated mass system is developed.
 The mechanism of rigidflexible coupling effect is expounded.
 The centrifugal, tangential and Coriolis effects on the tip motion trajectories are discussed.
 The influences of related parameters on the dynamics and vibrations are studied.
Keywords: rotating hubflexible rod, concentrated mass, rigidflexible coupling, dynamic response, free vibration.
1. Introduction
The dynamic problems (such as transient response and free vibration) of rotating flexible beams (rods) are the key research topics in a lot of engineering applications such as aerospace, robotics, propellers, turbine blades etc., which are crucial for design purposes, optimization and control. The study on the dynamic behavior and stability of these structures is an essential part of the design process, especially when dealing with the design of control and condition monitoring systems. Rotating beam (rod) structures generally include both the small elastic deflections of the flexible beam and the overall rotating motion of the hub, which are referred to the rigidflexible coupling mechanisms. The rotating hubflexible rod presented herein is a typical example. Therefore, more reliable mechanical models and their solution strategies for these coupling mechanisms have been required for accurate speed control and accurate operation [13].
The dynamic modeling of the rotating hubbeam (rod) structure has roughly gone through three stages. The dynamic model of the first stage was known as the traditional hybrid coordinate model [48], in which the small deformation hypothesis of the traditional structural dynamics was applied. It is assumed that the axial and lateral deformations of the beam were not coupled to each other. This model was generally used in the case that the large overall motion was known. In this case, only the influence of the overall motion of the hub on the elastic deformation of the rod was discussed while ignoring the influence of the flexible deformation of the rod on the rigid overall motions.
The dynamic models of the second stage were called the dynamic stiffening models, in which the dynamic stiffening effect was first proposed by Kane et al. [5] in 1987. It was demonstrated that the traditional hybrid coordinate model failed to handle the dynamics of rotary hubrod system with high angular velocity. Thereafter, many models were established to capture the dynamic stiffening effect, which can be roughly divided into three categories. The dynamic stiffening models of the first category added the additional potential energy induced by centrifugal force to describe the centrifugal stiffening effect [914]. In the second type of dynamic stiffening models [1522], the axial deformation is described by using the stretch deformation (i.e. a nonCartesian variable) and taking into account the secondorder coupling deformation caused by the transverse deformation. The dynamic stiffening models of the third category used a geometrically exact approach to describe the centrifugal stiffening effect [2326]. Although the introduction of dynamic stiffening enables the dynamic model to deal with the dynamics of highspeed rotating beam systems, there are still limitations on understanding of the rigidflexible coupling mechanism of such systems.
The dynamic model of the third stage is the rigidflexible coupling model (RFC model) [27], which is developed through the classical continuum mechanics. The RCF model considers that the essence of the dynamic stiffening effect is the structural dynamics of the noninertial system. The coupling effect of the large overall motion of the rigid body and the elastic deformation of the beam results in the additional stiffness, which leads to the dynamic stiffness of the system. Liu and Hong [28] investigated the dynamic stiffening effect on the dynamic behaviors of a elastic slender beam undergoing free overall motions based on the rigidflexible coupling theory. Cai et al. [29] developed the firstorder approximation coupling (FOAC) model for rotary hubbeamtip mass system and investigated its dynamics of the system when the rotation motion was unknown. Based on the Timoshenko beam theory, You et al. [30] proposed the dynamic model of a rotary flexible hubbeam system and studied the influence of the transverse shear deformation on the dynamic behaviors. On the basis of a new rigidflexible coupling dynamical model, Li and Zhang [31] investigated the dynamic responses and free vibrations of rotary tapered beams made of axially functionally graded materials by utilizing the Bspline approach. Liu et al. [32] presented a rigidflexible coupling dynamic model for a rotary rigid bodyflexible rodtip mass system. Since the mass of the flexible rod was neglected, the model presented by Liu et al. [32] was actually a twodegreeoffreedom model. Recently, by utilizing the slope angle of the beam, Zhang et al. [33] developed a new rigidflexible coupling dynamical model of a rotary hubflexible beamtip mass system under gravity loads, and determined the natural frequencies and nonlinear frequency responses through the incremental harmonic balance method.
In this study, a rigidflexible coupling dynamical model of a rotary hubflexible rodmovable concentrated mass system is developed by employing Hamilton’s principle. Based on the accurate geometric nonlinearity, the secondorder coupling deformations caused by the transverse deformations are employed to capture the dynamic stiffening effect. The derived partial differential governing equations are solved numerically by the Galerkin method. A comparative study is carried out to verify the present model. The dynamic responses of the system for situations of a prescribed rotation and an arbitrary external torque acting on the hub are presented. The influences of rotational speed ratio, concentrated mass location ratio, concentrated mass ratio and initial eccentricity ratio on the responses are studied. Finally, the natural frequencies of the rotating flexible rod under steadystate rotation are calculated.
2. Nonlinear dynamic modeling of the system
Fig. 1 shows a rotating hubflexible rodconcentrated mass system. The flexible rod with length $L$, crosssectional area $A$, mass density $\rho $, Young’s modulus $E$, is attached to the rigid hub, which rotates about the fixed vertical axis (i.e., $X$axis). Parameters ${J}_{H}$ and $\theta $ are the rotary inertia and rotation angle of the hub with respect to $X$axis. The $O$$XYZ$ is the inertial reference frame, whereas the $o$$xyz$ is the floating one that is fixed on the rod. In the present study, a concentrated mass $m$ considered can be moved along the centroid axis of the rod.
Since a slender rod is considered herein, the classical beam theory is used to describe its deformations. The deformation field of a generic point $P$ on the deformed rod with respect to the floating reference frame can be expressed by:
where $u(x,t)$, $v(x,t)$ and $w(x,t)$ are the deformation components at the centroid axis along $x$, $y$ and $z$ directions, respectively.
Fig. 1. A rotating hubflexible rodmass system
Based on the relation of the geometric nonlinearity, the axial deformation $u(x,t)$ could be described by using the arclength $s(x,t)$ and the transverse deformations, as [16, 21]:
in which the terms of $0.5{\int}_{0}^{x}(\partial v(\chi ,t)/\partial \chi {)}^{2}d\chi $ and $0.5{\int}_{0}^{x}(\partial w(\chi ,t)/\partial \chi {)}^{2}d\chi $ are called as the nonlinear secondorder coupling deformations.
Based on the von Kármán strain theory, the axial normal strain ${\epsilon}_{xx}$ is given by:
Utilizing the uniaxial stressstrain relations, the corresponding normal stress is:
Without considering the effect of gravity, the potential energy of the system can be written by:
where ${I}_{z}$ and ${I}_{y}$ are the cross section moments of inertia with respect to $z$ and $y$ axes, respectively.
The position vector of the point $P$ in the inertial reference frame can be written by:
in which, $\mathbf{A}$ is the rotation transformation matrix. ${\mathbf{r}}_{0}$ is the position vector of the point $o$ in the inertial frame. ${\mathbf{r}}_{1}$ is the position vector of the point $P$ on the undeformed rod in the floating frame. The matrix $\mathbf{A}$, the vectors ${\mathbf{r}}_{0}$ and ${\mathbf{r}}_{1}$ are given by:
Accordingly, the absolute velocity of the point $P$ can be derived by:
where $\dot{\theta}$ is the rotational speed. The matrix $\stackrel{~}{\mathbf{I}}$ can be expressed by:
Therefore, the kinetic energy of the system is:
An external rotating torque ${E}_{tor}$ is assumed to act on the hub and thus external virtual work is derived by:
where ${p}_{v}(x,t)$ and ${p}_{w}(x,t)$ are the external distributed forces acted on the rod along the $y$ and $z$ directions.
The rigidflexible coupling dynamic equations of the rotating hubflexible rod system are derived by the following Hamilton’s principle, as:
Substituting Eqs. (5), (10) and (11) into Eq. (12), one can obtain the comprehensive nonlinear differential equations of the system:
$={\int}_{0}^{L}\rho A\left({r}_{0}\ddot{w}+v\ddot{w}+w\ddot{v}\right)dx+{E}_{tor},$
where ${J}_{1}={\int}_{A}\rho ({y}^{2}+{z}^{2})dA$.
The corresponding boundary conditions of the rod are give by:
– Geometric boundary conditions: $u=v=w=\frac{\partial v}{\partial x}=\frac{\partial w}{\partial x}=0$ at $x=0$.
– Force boundary conditions: $\frac{\partial u}{\partial x}=\frac{{\partial}^{2}v}{\partial {x}^{2}}=\frac{{\partial}^{2}w}{\partial {x}^{2}}=0$ at $x=L$.
It is seen from Eq. (13) that rigid rotation of the hub and the elastic deformation of the slender rod are coupled with each other. This is referred to the rigidflexible coupling effect, which leads to the strong nonlinearity of the system. In the rest of the present study, the axial motion of the rod is neglected since it is not coupled with the radial motion $v$ and tangential motion $w$ of the rod and the rotation of the hub.
Then, a Dirac’s delta function is applied to take into account a movable concentrated mass on the rod, as follows:
where $b$ is the distance from the mass to the root of the rod. The values of the Dirac’s delta function are given by:
Using ${\rho}^{*}\left(x\right)$ instead of $\rho A$, one can rewrite the Eqs. 13(a), 13(c) and 13(d), as:
$={\int}_{0}^{L}{\rho}^{*}\left({r}_{0}\ddot{w}+v\ddot{w}+w\ddot{v}\right)dx+{E}_{tor},$
The terms of $2{\rho}^{*}\dot{\theta}\dot{w}$ and $2{\rho}^{*}\dot{\theta}\dot{v}$ in Eqs. (16b) and (16c) represent the Coriolis effect of the rotating rod. The centrifugal stiffening effect of the rotating rod is reflected by the terms of ${\rho}^{*}{\dot{\theta}}^{2}v$ and ${\rho}^{*}{\dot{\theta}}^{2}w$. The terms of ${\rho}^{*}\ddot{\theta}w$ and ${\rho}^{*}\ddot{\theta}v$ represent the tangential effect. The terms of ${\rho}^{*}{\dot{\theta}}^{2}{r}_{0}$ and ${\rho}^{*}\ddot{\theta}{r}_{0}$ in right hand of the Eqs. (16b) and (16c) are, respectively, the centrifugal force and tangential force. Obviously, both centrifugal and tangential forces depend on the rotation of the hub and the initial eccentricity of the rod (i.e. ${r}_{0}$).
The Galerkin method is employed to solve the Eq. (16) numerically, hence the deformations of $v$ and $w$ are approximately given by:
where $N$ is the number of terms of the trial functions. ${c}_{i}\left(t\right)$ and ${e}_{i}\left(t\right)$ denote the timedependent generalized coordinates of the radial and tangential deformations, respectively. ${V}_{i}\left(x\right)$ and ${W}_{i}\left(x\right)$ are the trial functions satisfying both the natural and essential boundary conditions. In the present investigation, the transverse bending mode functions of a stationary cantilever EulerBernoulli beam are chosen as the trial functions, as:
where ${\lambda}_{i}$ is $i$th root of:
Before applying the Galerkin method to discrete the obtained governing equations, the variational equations of motion should be derived by multiplying Eq. (16b) and (16c) by corresponding weighting functions for the radial and tangential deformations, as:
where ${\widehat{c}}_{i}\left(t\right)$ and ${\widehat{e}}_{i}\left(t\right)$ are the arbitrary timedependent functions. Applying integration by parts to Eqs. 16(c) and 16(d) over length $L$, the variational equations can be rewritten by:
$={\int}_{0}^{L}{\rho}^{*}\left({r}_{0}\ddot{w}+v\ddot{w}+w\ddot{v}\right)dx+{E}_{tor}m,$
Within the process of presentation of the numerical results, dimensionless analysis is preferable. To this end, the following dimensionless parameters are defined:
where $\delta $, $\alpha $, $\beta $, $\gamma $, ${\stackrel{}{E}}_{rot}$ and ${\eta}_{H}$ represent the initial eccentricity ratio, concentrated mass ratio, concentrated mass location ratio, rotational speed ratio, external torque ratio and rotary inertia of the hub, respectively. Moreover, $I={I}_{z}={I}_{y}$ is fixed to be used.
By using the above dimensionless parameters, the discretized governing equations of motion can be written in a dimensionless matrix form:
where $\mathbf{M}$, $\mathbf{C}$ and $\mathbf{K}$ represent the generalized mass, Coriolis and stiffness matrices, respectively. $\mathbf{f}$ and $\stackrel{}{\mathbf{q}}$ are the generalized force and coordinate arrays, respectively:
$\mathbf{K}=\left[\begin{array}{lll}0& 0& 0\\ 0& {\mathbf{K}}^{vv}& 0\\ 0& 0& {\mathbf{K}}^{ww}\end{array}\right],\mathbf{f}=\left[\begin{array}{l}{f}^{\theta}\\ {\mathbf{f}}^{v}\\ {\mathbf{f}}^{w}\end{array}\right],\stackrel{}{\mathbf{q}}=\left[\begin{array}{l}\theta \\ \stackrel{}{\mathbf{c}}\\ \stackrel{}{\mathbf{e}}\end{array}\right],$
in which the elements in the submatrices ${\mathbf{M}}^{\varphi \phi}$, ${\mathbf{C}}^{\varphi \phi}$, ${\mathbf{K}}^{\varphi \phi}$ and ${\mathbf{f}}^{\varphi}$ ($\phi ,\varphi =\theta ,v,w$) are given by:
$+\sum _{i=1}^{N}\sum _{j=1}^{N}{M}_{ij}^{vv}{\stackrel{}{c}}_{i}{\stackrel{}{c}}_{j}+\sum _{i=1}^{N}\sum _{j=1}^{N}{M}_{ij}^{ww}{\stackrel{}{e}}_{i}{\stackrel{}{e}}_{j},$
${M}_{i}^{\theta v}={M}_{i}^{v\theta}=\sum _{j=1}^{N}{C}_{ij}^{vw}{\stackrel{}{e}}_{j},{M}_{i}^{\theta w}={M}_{i}^{w\theta}=\sum _{j=1}^{N}{C}_{ij}^{wv}{\stackrel{}{c}}_{j}+\delta \left({\int}_{0}^{1}{W}_{i}\left(\xi \right)d\xi +\alpha {W}_{i}\left(\beta \right)\right),$
${M}_{ij}^{vv}={\int}_{0}^{1}{V}_{i}\left(\xi \right){V}_{j}\left(\xi \right)d\xi +\alpha {V}_{i}\left(\beta \right){V}_{j}\left(\beta \right),{M}_{ij}^{ww}={\int}_{0}^{1}{W}_{i}\left(\xi \right){W}_{j}\left(\xi \right)d\xi +\alpha {W}_{i}\left(\beta \right){W}_{j}\left(\beta \right),$
${C}_{ij}^{vw}={\int}_{0}^{1}{V}_{i}\left(\xi \right){W}_{j}\left(\xi \right)d\xi \alpha {V}_{i}\left(\beta \right){W}_{j}\left(\beta \right),{C}_{ij}^{wv}={\int}_{0}^{1}{W}_{i}\left(\xi \right){V}_{j}\left(\xi \right)d\xi +\alpha {W}_{i}\left(\beta \right){V}_{j}\left(\beta \right),$
${K}_{ij}^{vv}={\int}_{0}^{1}{V}_{i}^{\text{'}\text{'}}\left(\xi \right){V}_{j}^{\text{'}\text{'}}\left(\xi \right)d\xi {\gamma}^{2}{M}_{ij}^{vv},{K}_{ij}^{ww}={\int}_{0}^{1}{W}_{i}^{\text{'}\text{'}}\left(\xi \right){W}_{j}^{\text{'}\text{'}}\left(\xi \right)d\xi {\gamma}^{2}{M}_{ij}^{ww},$
${f}^{\theta}={\stackrel{}{E}}_{rot}2\gamma \left[\sum _{i=1}^{N}\sum _{j=1}^{N}{M}_{ij}^{vv}{\stackrel{}{c}}_{i}{\dot{\stackrel{}{c}}}_{j}+\sum _{i=1}^{N}\sum _{j=1}^{N}{M}_{ij}^{ww}{\stackrel{}{e}}_{i}{\dot{\stackrel{}{e}}}_{j}+\delta \sum _{i=1}^{N}\left({\int}_{0}^{1}{V}_{i}\left(\xi \right)d\xi +\alpha {V}_{i}\left(\beta \right)\right)\right],$
${f}_{i}^{v}={\int}_{0}^{1}{\stackrel{}{p}}_{v}{V}_{i}\left(\xi \right)d\xi +{\gamma}^{2}\delta \left({\int}_{0}^{1}{V}_{i}\left(\xi \right)d\xi +\alpha {V}_{i}\left(\beta \right)\right)+2\gamma \sum _{j=1}^{N}{C}_{ij}^{vw}{\dot{\stackrel{}{e}}}_{j}\sum _{j=1}^{N}\left({\int}_{0}^{1}{{V}^{\text{'}\text{'}}}_{i}\left(\xi \right){{V}^{\text{'}\text{'}}}_{j}\left(\xi \right)d\xi +{\gamma}^{2}{M}_{ij}^{vv}\right){\stackrel{}{c}}_{j},$
If the rotation of the hub is prescribed, the influence of the elastic motion of the rod on the rotation of the hub can be neglected. Therefore, the dynamic model of the system for an arbitrary prescribed rotation can be obtained from Eq. (23), as follows:
In the following study, it is assumed that no external force acts on the rod (i.e., ${p}_{v}={p}_{w}=0$).
3. Dynamic response analysis of system
3.1. Comparative study
Before studying the dynamic response properties of the rotary hubflexible rodconcentrated mass system, a comparative study should be presented between the present study and Liu et al. [32], in which the mass of the rod was neglected and the concentrated mass is located at the end of the rod. In essence, the model developed by Liu et al. [32] was a twodegreeoffreedom model. The physical parameters given by Liu et al. [32] are tabulated in Table 1 alongside the mass density of the rod. The rotation of the hub is prescribed by $\dot{\theta}=\omega =\varpi t/15(\varpi /2\pi \left)\mathrm{s}\mathrm{i}\mathrm{n}\right(2\pi t/15)$ in the time region of 0 s $\le t\le $15 s. Fig. 2 shows the tip motion trajectories of the rod obtained from the present study and Liu et al. [32] with $\varpi =$ 2 rad/s. As shown in Fig. 2, the discrepancy between the results of the present and Liu’s models is very small. This not only proves the correctness of the present model, but also illustrates that the mass of the rod in this simulation is negligible.
Fig. 2. Tip motion trajectories of the rotating rod ($m=$ 0.3 kg, $L=$ 0.3 m and $\varpi =$ 2 rad/s)
Fig. 3. Tip motion trajectories of the rotating rod ($m=$ 0.1 kg, $L=\mathrm{}$0.8 m and $\varpi =$ 1 rad/s)
Fig. 3 shows the tip motion trajectories of the rod obtained from the present study and Liu et al. [32] with $\varpi =$ 1 rad/s. In this calculation, the concentrated mass and the length of the rod are set to 0.1 kg and 0.8 m, respectively, and the other parameters are the same as those in Table 1. It can be seen from Fig. 3 that the results of the present model and Liu’s model [32] have great differences. This reveals that the effect of the mass of the rod on the exact dynamic analysis cannot be neglected.
Table 1. The physical parameters of the rotating hubrodconcentrated mass system for comparative study
Concentrated mass

Rod
length

Rod
diameter

Young’s modulus

Eccentricity
of rod

Mass density
of rod

$m=$ 0.3 kg

$L=$ 0.3 m

$d=$ 0.003 m

$E=$ 100 GPa

${r}_{0}=$ 0.1 m

$\rho =$7.2×10^{3} kg/m^{3}

3.2. Dynamic response analysis for a prescribed rotation
The Eq. (26) is used to determine the elastic deformations of the rod when the rotation of the hub is assumed to be known. A prescribed rotating motion of the hub is defined by:
where ${\gamma}_{0}$ denotes the steadystate rotational speed ratio and $T$ represents a given dimensionless time parameter. Here ${\gamma}_{0}=$ 0.2 and $T=$ 15 are used in this section. According to Eq. (27), the prescribed motion can be divide into four different motion processes. The hub is initially at rest and then a spinup process is performed from 0 to 15. The steadystate rotational speed ratio is maintained when $\text{15}\le \tau \le \text{30}$ and the steadystate rotation follows a spindown process when $\text{30}\le \tau \le \text{45}$, finally the rotating motion stops after $\tau =$ 45. The variations of the rotational speed ratio and the rotational acceleration ratio of the hub are illustrated in Fig. 4. The rotational acceleration ratio undergoes two updown changes in the positive and negative directions during $\text{0}\le \tau \le \text{15}$ and $\text{30}\le \tau \le \text{45}$, respectively.
Fig. 4. Variations of a) the rotational speed ratio and b) the rotational acceleration ratio of the hub
Fig. 5 shows the tip motion trajectory of the rod throughout the whole process of the hub motion prescribed in Eq. (27) with $\alpha =$ 0.1, $\beta =$ 1 and $\delta =$ 0.1. During the spinup process (i.e., $\text{0}\le \tau \le \text{15}$), the tip of the rod moves tangentially backwards ($w/L$ is negative) and radially outwards, as shown in Fig. 6(a). The radial deformation gradually increases and the tangential deformation increases first and then decreases since the centrifugal force depending on the rotational speed ratio increases and the tangential force depending on the rotational acceleration ratio undergoes an updown change.
Fig. 5. Tip motion trajectory of the rod for the whole process of the hub motion ($\gamma =\text{0.2}$, $T=\text{15}$, $\alpha =\text{0.1}$, $\beta =\text{1}$ and $\delta =\text{0.1}$)
Fig. 6. Tip motion trajectory of the rod for the four different motion processes ($\gamma =$ 0.2, $T=$ 15, $\alpha =$ 0.1, $\beta =$ 1 and $\delta =$ 0.1): a) the spinup process ($\text{0}\le \tau \le \text{15}$), b) the steadystate rotation process ($\text{15}\le \tau \le \text{30}$), c) the spindown process ($\text{30}\le \tau \le \text{45}$), d) the static process ($\tau \ge $ 45)
Fig. 6(b) shows the tip motion trajectory of the rod in the steadystate rotation process (i.e., $\text{15}\le \tau \le \text{30}$). An attractive petal pattern is presented since the oval path of the tip trajectory rotates clockwise around a fixed point, which is located in the radial direction. The radial deformation is about two orders of magnitude larger than the tangential deformation. This should be attributed to the fact that the tangential force vanishes due to the steadystate rotation. Therefore, the petal pattern is only affected by Coriolis force.
In the spindown process (i.e., $\text{30}\le \tau \le \text{45}$), the tip motion trajectory of the rod is displayed in Fig. 6(c). The tip of the rod moves radially inwards and tangentially forwards ($w/L$ is positive). The radial deformation gradually decreases and the tangential deformation first increases and then decreases.
In the final process (i.e., $\tau \ge $ 45), the hub is at rest so that the motion of the rod is only influenced by its own elastic restoring force. A steadystate tip elliptical motion trajectory of the rod is depicted in Fig. 6(d). This is a typical twodimensional steadystate oscillation phase diagram.
Fig. 7. Effects of parameters on the tip motion trajectory of the rod: a) effect of concentrated mass ratio, b) effect of concentrated mass location ratio, c) effect of initial eccentricity ratio
Figs. 7(a), 7(b) and 7(c) show the effects of concentrated mass, concentrated mass location and initial eccentricity ratios on the tip motion trajectories of the rod, respectively. As shown in Fig. 7, both the radial and tangential deformations increase as the concentrated mass, concentrated mass location and initial eccentricity ratios increase due to the increase of the centrifugal force and the tangential force.
3.3. Dynamic response analysis for an arbitrary external torque acting on the hub
For the case of an arbitrary external torque acting on the hub, the rotation of the hub is unknown and should be solved. Evidently, the rigid rotation of the hub can lead to the elastic deformation of the rod, and the elastic deformation of the rod will also affect the rigid rotation of the hub. The Eq. (23) is used to determine the elastic deformations of the rod and the rigid rotation of the hub. The following external torque is assumed to act on the hub:
where ${\stackrel{}{E}}_{tor}^{0}$ is a given external torque ratio. Here, ${\stackrel{}{E}}_{tor}^{0}=$ 1 and $T=$ 40 are fixed to be used. According to Eq. (28), the external torque is acted on the hub following a sinusoidal law during $\text{0}\le \tau \le \text{40}$ and the external torque is removed when $\tau \ge \text{40}$. A numerical example for a rotating hubflexible circular crosssection rod with a concentrated mass is presented below when $\alpha =$ 0.5, $\beta =$ 1, $\delta =$ 0.05, $d/L=$ 0.01 and ${\eta}_{H}=$ 20.
Figs. 8(a) and 8(b) show the responses of the rotation angle $\theta $ and the rotational speed ratio $\gamma $ of the hub, respectively. As shown in Fig. 8(a), the rotation angle gradually increases when $\text{0}\le \tau \le \text{40}$ and then fluctuates within a very small range and decreases slightly when $\text{40}\le \tau \le \text{50}$. The rotational speed ratio of the hub increases when $0\le \tau \le 20$, then decreases when $\text{20}\le \tau \le \text{40}$ and finally oscillates around zero when $\text{40}\le \tau \le \text{50}$. As a matter of fact, both the fluctuation of the rotation angle and the oscillation of the rotational speed ratio when $\text{40}\le \tau \le \text{50}$ that directly reflect the influences of elastic deformations of the rod on the motion of the hub.
Fig. 8. Rotation responses of the hub (${\stackrel{}{E}}_{tor}^{0}=$ 1, $T=$ 40, $\alpha =$ 0.5, $\beta =$ 1, $\delta =$ 0.05, $d/L=$0.01 and ${\eta}_{H}=$ 20): a) the rotation angle, b) the rotational speed ratio
Figs. 9(a) and 9(b) show the tip motion trajectories of the rod for the whole process (i.e., $\text{0}\le \tau \le \text{50}$) and the no external torque process (i.e., $\text{40}\le \tau \le \text{50}$), respectively. As shown in Fig. 9(a), when $\text{0}\le \tau \le \text{20}$, the radial deformation gradually increases and the tangential deformation increases first and then decreases to zero. When $\text{20}\le \tau \le \text{40}$, the radial deformation gradually decreases and the tangential deformation increases first and then decreases. When $\text{40}\le \tau \le \text{50}$, a steadystate tip elliptical motion trajectory is depicted in Fig. 9(b). Obviously, the rotation of the hub causes the corresponding elastic deformations of the rod.
Fig. 9. Tip motion trajectory of the rod (${\stackrel{}{E}}_{tor}^{0}=$ 1, $T=$ 40, $\alpha =$ 0.5, $\beta =$ 1, $\delta =$ 0.05, $d/L=$ 0.01 and ${\eta}_{H}=$ 20): a) the whole process ($0\le \tau \le \text{50}$), b) the no external torque process ($\text{40}\le \tau \le \text{50}$)
4. Free vibration analysis for steadystate rotation
In this section, the hub rotates at a constant rotational speed (i.e., $\dot{\gamma}=\text{0}$). The free vibration properties of the rotating rod with a concentrated mass could be explored by transforming Eq. (26) into a state space equation, as follows:
where:
$\mathbf{r}={\left[\begin{array}{l}\dot{\stackrel{}{\mathbf{c}}}\\ \dot{\stackrel{}{\mathbf{e}}}\\ \stackrel{}{\mathbf{c}}\\ \stackrel{}{\mathbf{e}}\end{array}\right]}_{4N\times 1},$
in which $\mathbf{r}$ and $\mathbf{I}$ are the state space vector and an identity matrix, respectively. Substituting $\mathbf{r}={e}^{j\mu \tau}\mathbf{R}$ into Eq. (29) yields:
where $\mu $ and $\mathbf{R}$ are the eigenfrequency and the eigenvector, respectively.
Fig. 10 presents the first four nondimensional natural frequencies of the rotating rodconcentrated mass system versus rotational speed ratio with $\alpha =$ 0.1 and $\beta =$ 1 alongside the results reported by Liu et al. [32]. As shown in Fig. 10, when $\gamma =$ 0 (i.e., the hub is at rest), the first and second dimensionless natural frequencies are equal to 2.9678, whereas the results of Liu’s model are 5.4772. Obviously, the consideration of the mass of the rod reduces the natural frequencies of the system. With the increase of the rotational speed ratio, the first and third natural frequencies decrease while the second and fourth natural frequencies increase. The dimensionless fundamental natural frequency vanishes at $\gamma =$ 2.9678, which is called as the buckling speed ratio. Evidently, the buckling speed ratio of system is equal to the first dimensionless natural frequency of $\gamma =$ 0. Furthermore, an interesting relationship between the adjacent dimensionless natural frequencies can be obtained as:
where ${\mu}_{2i}$ and ${\mu}_{2i1}$ are the 2$i$th and ($2i1$)th nondimensional natural frequencies of the system with the rotational speed ratio $\gamma $, respectively.
Fig. 10. The first four dimensionless natural frequencies of the rotating rodconcentrated mass system versus rotational speed ratio ($\alpha =$ 0.1 and $\beta =$ 1)
The variation of the dimensionless fundamental natural frequency of the rotating rodconcentrated mass system versus concentrated mass location $\beta $ is depicted in Fig. 11 with $\gamma =$ 0.2 for five values of $\alpha =$ 0, 0.1, 0.5, 1.0, 2.0. It is observed from Fig. 11 that the dimensionless fundamental natural frequencies decrease as the concentrated mass location ratio $\beta $ and/or the concentrated mass ratio $\alpha $ increase.
Fig. 11. The variation of the dimensionless fundamental natural frequency of the rotating rodconcentrated mass system versus concentrated mass location ratio ($\gamma =$ 0.2)
5. Conclusions
In this study, a comprehensive rigidflexible coupling dynamic model of a rotating hubflexible rod with a concentrated mass is established to investigate the dynamic characteristics of the system. On the basis of the classical beam theory and von Kármán geometric nonlinearity, the governing equations of the system are derived by using the Hamilton’s principle. The Galerkin method is employed to solve the derived equations numerically. The rigidflexible coupling effect between the rigid rotation of the hub and the flexible deformation of the rod leads to the strong nonlinearity of the system. For the case of the prescribed rigid rotation, the dynamic responses of the rod are numerically calculated, and the tip motion trajectories of the rod are graphically presented. The centrifugal, tangential and Coriolis effects on the tip motion trajectories are discussed. The attractive petal pattern for the steadystate rotation and the conventional elliptical pattern for the case without rotation are shown. For the case of the external torque acting on the hub, a numerical example is presented to illustrate the rigidflexible coupling effect of the system. For the case of steadystate rotation, as the rotational speed ratio increases, the first and third dimensionless natural frequencies decrease and the second and fourth ones increase due to the Coriolis effect. The dimensionless fundamental natural frequency vanishes at the buckling speed ratio. The dimensionless fundamental natural frequency decreases as the concentrated mass location ratio and/or the concentrated mass ratio increase. The modeling method and related conclusions in this paper can provide theoretical basis for the design, optimization and dynamic control of such rotating mechanisms in industrial applications.
Acknowledgements
This research work is supported by the National Natural Science Foundation of China (Grant No. 11872031) and the Natural Science Foundation of Nanjing Institute of Technology (Grant No. ZKJ201701).
References
 Shabana A. Flexible multibody dynamics: review of past and recent developments. Multibody System Dynamics, Vol. 1, 1997, p. 189222. [Publisher]
 Wasfy T., Noor A. Computational strategies for flexible multibody systems. Applied Mechanics Reviews, Vol. 56, Issue 6, 2003, p. 553613. [Publisher]
 Schielen W. Computational dynamics: theory and applications of multibody systems. European Journal of Mechanics ASolids, Vol. 25, Issue 4, 2006, p. 566594. [Publisher]
 Singh R., Vandervoort R., Likins P. Dynamics of flexible bodies in tree topologya computeroriented approach. Journal of Guidance Control and Dynamics, Vol. 8, Issue 5, 1985, p. 584590. [Publisher]
 Kane T. R., Ryan R., Banerjee A. K. Dynamics of a cantilever beam attached to a moving base. Journal of Guidance Control and Dynamics, Vol. 10, Issue 2, 1987, p. 139151. [Publisher]
 Krishnaprasad P. S., Marsden J. E. Hamiltonian structures and stability for rigid bodies with flexible attachments. Archive for Rational Mechanics and Analysis, Vol. 98, Issue 1, 1987, p. 7193. [Publisher]
 Meirovitch L. Hybrid state equations of motion for flexible bodies in terms of quasicoordinates. Journal of Guidance Control and Dynamics, Vol. 14, Issue 5, 1991, p. 10081013. [Publisher]
 Choura S., Jayasuriya S., Medick M. A. On the modeling, and openloop control of a rotating thin flexible beam. Journal of Dynamic Systems, Measurement, and Control, Vol. 113, Issue 1, 1991, p. 2633. [Publisher]
 Yokoyama T. Free vibration characteristics of rotating Timoshenko beams. International Journal of Mechanical Sciences, Vol. 30, Issue 10, 1988, p. 743755. [Publisher]
 Shahba A., Attarnejad R., Hajilar S. Free vibration and stability of axially functionally graded tapered EulerBernoulli beams. Shock and Vibration, Vol. 18, Issue 5, 2011, p. 683696. [Publisher]
 Banerjee J. R., Kennedy D. Dynamic stiffness method for inplane free vibration of rotating beams including Coriolis effects. Journal of Sound and Vibration, Vol. 333, Issue 26, 2014, p. 72997312. [Publisher]
 DehrouyehSemnani A.M. The influence of size effect on flapwise vibration of rotating microbeams. International Journal of Engineering Science, Vol. 94, 2015, p. 150163. [Publisher]
 Ghafarian M., Ariaei A. Free vibration analysis of a system of elastically interconnected rotating tapered Timoshenko beams using differential transform method. International Journal of Mechanical Sciences, Vol. 107, 2016, p. 93109. [Publisher]
 Azimi M., Mirjavadi S. S., Shafiei N., Hamouda A. M. S., Davari E. Vibration of rotating functionally graded Timoshenko nanobeams with nonlinear thermal distribution. Mechanics of Advanced Materials and Structures, Vol. 25, Issue 6, 2018, p. 467480. [Publisher]
 Wright A. D., Smith C. E., Thresher R. W., Wang J. L. C. Vibration modes of centrifugally stiffened beams. Journal of Applied Mechanics, Vol. 49, Issue 1, 1982, p. 197202. [Publisher]
 Yoo H. H., Shin S. H. Vibration analysis of rotating cantilever beams. Journal of Sound and Vibration, Vol. 212, Issue 5, 1998, p. 807828. [Publisher]
 Li L., Zhang D. G., Zhu W. D. Free vibration analysis of a rotating hubfunctionally graded material beam system with the dynamic stiffening effect. Journal of Sound and Vibration, Vol. 333, Issue 5, 2014, p. 15261541. [Publisher]
 Fang J., Zhou D. Free vibration analysis of rotating axially functionally graded tapered Timoshenko beams. International Journal of Structural Stability and Dynamics, Vol. 16, Issue 5, 2016, p. 1550007. [Publisher]
 Oh Y., Yoo H. H. Vibration analysis of rotating pretwisted tapered blades made of functionally graded materials. International Journal of Mechanical Sciences, Vol. 119, 2016, p. 6879. [Publisher]
 Zhao G., Wu Z. Coupling vibration analysis of rotating threedimensional cantilever beam. Computers and Structures, Vol. 179, 2017, p. 6474. [Publisher]
 Fang J., Gu J., Wang H. Sizedependent threedimensional free vibration of rotating functionally graded microbeams based on a modified couple stress theory. International Journal of Mechanical Sciences, Vol. 136, 2018, p. 188199. [Publisher]
 Fang J., Gu J., Wang H., Zhang X. Thermal effect on vibrational behaviors of rotating functionally graded microbeams. European Journal of Mechanics ASolids, Vol. 75, 2019, p. 497515. [Publisher]
 Huang C. L., Lin W. Y., Hsiao K. M. Free vibration analysis of rotating Euler beams at high angular velocity. Computers and Structures, Vol. 88, Issues 1718, 2010, p. 9911001. [Publisher]
 Arvin H., Bakhtiari Nejad F. Nonlinear modal analysis of a rotating beam. International Journal of NonLinear Mechanics, Vol. 46, Issue 6, 2011, p. 877897. [Publisher]
 Lacarbonara W., Arvin H., BakhtiariNejad F. A geometrically exact approach to the overall dynamics of elastic rotating blades – part 1: linear modal properties. Nonlinear Dynamics, Vol. 70, Issue 1, 2012, p. 659675. [Publisher]
 Kim H., Yoo H. H., Chung J. Dynamic model for free vibration and response analysis of rotating beams. Journal of Sound and Vibration, Vol. 332, Issue 22, 2013, p. 59175928. [Publisher]
 Yang H., Hong J., Yu Z. Dynamics modelling of a flexible hub–beam system with a tip mass. Journal of Sound and Vibration, Vol. 266, Issue 4, 2003, p. 759774. [Publisher]
 Liu J. Y., Hong J. Z. Geometric stiffening effect on rigidflexible coupling dynamics of an elastic beam. Journal of Sound and Vibration, Vol. 278, Issues 45, 2004, p. 11471162. [Publisher]
 Cai G. P., Hong J. Z., Yang S. X. Dynamic analysis of a flexible hubbeam system with tip mass. Mechanics Research Communications, Vol. 32, Issue 2, 2005, p. 173190. [Publisher]
 You C., Hong J., Cai G. Modeling study of a flexible hubbeam system with large motion and with considering the effect of shear deformation. Journal of Sound and Vibration, Vol. 295, Issues 12, 2006, p. 282293. [Publisher]
 Li L., Zhang D. Dynamic analysis of rotating axially FG tapered beams based on a new rigidflexible coupled dynamic model using the Bspline method. Composite Structures, Vol. 124, 2015, p. 357367. [Publisher]
 Liu Z., Ye P. X., Guo Y., Guo X. Rigidflexible coupling dynamic analysis on a mass attached to a rotating flexible rod. Applied Mathematical Modelling, Vol. 38, Issues 2122, 2014, p. 49854994. [Publisher]
 Zhang D. W., Liu J. K., Huang J. L., Zhu W. D. Periodic responses of a rotating hubbeam system with a tip mass under gravity loads by the incremental harmonic balance method. Shock and Vibration, Vol. 2018, 2018, p. 8178274. [Search CrossRef]