Finite element modeling and active vibration control of high-speed spinning flexible beam

Zhou, Lanwei; Chen, Guoping; Yang, Jingyu

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 September 2015

Check for updates

Finite element modeling and active vibration control of high-speed spinning flexible beam

Lanwei Zhou¹

Guoping Chen²

Jingyu Yang³

^{1, 2}State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing, China

³Department of Aeronautics and Astronautics, Shenyang, China

Corresponding Author:

Lanwei Zhou

Cite the article Download PDF

Downloads 1340

Abstract

Finite element modeling and active vibration control of a high-speed spinning flexible coupled electromechanical beam is investigated using a first-order approximation coupling (FOAC) model. Due to centrifugal forces caused by eccentricity in a spinning flexible beam, there exists coupling between axial and transverse vibration modes. The partial differential equations of motion of the beam governing this coupling are derived using Hamilton’s principle based on an FOAC model, and a finite element method for discretization is given. It is observed that the zero-order approximate coupling (ZOAC) model is valid for dynamic description of the flexible beam spinning at low speeds, but no longer valid at high speeds. However, the validity of FOAC model is confirmed at different speeds. Piezoelectric elements for active vibration control of the spinning flexible beam are analyzed and a velocity feedback controller is proposed. Simulation results demonstrate good performance of the proposed velocity feedback controller.

1. Introduction

Spinning beams play an important role in many branches of engineering, for instance, shafts of combustion engines, DC and AC motors, turbine blades, propellers, spindles of gyro-scopes for control of satellites and satellite booms. Accurate prediction of the dynamics of such structures is critical for accurate design and analysis of the modern engineering systems, manufacturing accuracy and machine efficiency. Free as well as forced vibrations of beams with an orientation in or perpendicular to the axis of rotation have been investigated since the earliest days of artificial satellites. For example, the spin orientation of Explorer I was changed due to the energy dissipation of the flexible antennas [1]. However, an accurate analysis of a flexible beam spinning at high speeds is not available as yet.

An appropriate beam theory and an accurate model are important for precise prediction of motion of such spinning beams [2]. Early researches on the spinning beam were based on the classical beam, or Euler-Bernoulli beam description. Dimentberg [3] was among the first to investigate the dynamics of pinned-pinned spinning Euler-Bernoulli beams. Bauchau and Eidel [4] studied the dynamic response and stability analysis of spinning beams undergoing large amplitude and rotations using the finite-element method in time and discussed the classical nonlinear beam problems.

It is generally known that, for a non-spinning isotropic beam, however, rotary inertia and transverse shear deformation have about the same order of influence on the dynamic predictions [5, 6]. Due to the introduction of gyroscopic, Coriolis and centrifugal effects in a spinning beam, inertia becomes more important. Therefore, it is essential to use the Rayleigh beam theory in study of the dynamics of spinning beams [7, 8].

Sheu and Yang [9] studied the critical speed, whirl speed and mode shape of a spinning beam using a Rayleigh beam model with rotatory inertia and gyroscopic effects. Qian et al. [10] investigated the non-contact dynamic testing of a highly flexible spinning vertical shaft and time frequency analysis through a camera-based 3D motion analysis system and the Hilbert-Huang Transform (HHT), respectively. Frank et al. [2] derived linear and nonlinear models of spinning Rayleigh beams and investigated dynamic characteristics of downward vertical spinning Rayleigh beams, respectively.

Besides, some researchers adopted Timoshenko beam theory to describe spinning beams. Such as Boukhalfa [11], who analyzed the free vibrations of the spinning composite shafts using the hp-version of the finite element method. Lee et al. [12] studied a spectral element model for the spinning Timoshenko shaft. Mustapha and Zhong [13] analyzed the free vibration of a circular doubly-symmetric spinning micro beam embedded in an elastic medium using a mathematical model and premised the formulation based on the higher-order modified couple stress theory. Chan and Wang [14] investigated the revolving superposed standing waves in a spinning Timoshenko beam, and showed that the revolving waves should be represented by wave functions in a form of four-component column matrix vectors.

In general, a spinning flexible beam usually results in large nonlinearity and large deflections, especially at high speeds which becomes necessary to investigate. Khadem and Shahgholi [15] studied the two-mode combination resonances of a simply supported rotating shaft undergoing large deflections, rotary inertia and gyroscopic effects. Shahgholi et al. [16] studied simultaneous primary resonances of a simply supported nonlinear spinning beam due to unbalance periodic force and parametric resonances arising from excitations due to unequal mass moments of inertia and flexural rigidities in the direction of the principal axes. Shahgholi et al. [17] investigated parametric and main resonances of an asymmetrical spinning beam with in-extensional nonlinearity and large deflections.

The traditional hybrid coordinate model of flexible system is considered as a zero-order approximate dynamic model as it neglects high order coupling deformation kinematics. All of the existing studies ignored the effects of rigid-flexible coupling of the spinning beam. While this may not be an important influencing factor at low spinning speeds, it may lead to erroneous results in some high speed cases. Therefore, it is important to include the effects of rigid-flexible coupling.

In this paper, the effects of rigid-flexible coupling are considered by using a FOAC model to describe the dynamic behavior of the spinning Rayleigh beam. A new formulation of Hamilton’s principle and a finite element method for the case of a spinning flexible beam are proposed. In addition, active vibration control of a high-speed spinning flexible mechanical-electric coupled beam is investigated using velocity feedback control strategy. Numerical simulations are also carried out to verify the validity of FOAC model and the effectiveness of the velocity feedback controller. The comparison of the displacement responses derived from the ZOAC model and the proposed model not only illustrates the deficiency of the ZOAC model, but also validates the proposed model.

This paper is organized as follows: Section 2 first presents the expression of FOAC model for a high-speed spinning flexible beam and the coupled governing partial differential equations of motion of the beam which is subjected to a three-directional force using Hamilton’s principle. Finite element modeling of the above-mentioned beam is given in Section 3, followed by the active vibration control implementation in Section 4. Section 5 provides numerical simulation and comparison studies using FOAC model and ZOAC model, and control simulations by the proposed controller. Finally, concluding remarks are given in Section 6.

2. Mathematical modeling

2.1. First-order approximation model

Consider a spinning square Rayleigh beam with rotatory inertia and gyroscopic effects, as shown in Fig. 1. It is spinning with constant angular velocity $Ω$ about its longitudinal axis, the length is $L$ , and the density is $ρ$ with a prescribed compressive axial load $P$ at both ends. Assuming the displacement of a point on the beam $D = u i_{x} + v i_{y} + w i_{z}$ , the actual displacement along its longitudinal axis based on continuum medium mechanics [18] can be written as:

1

s_{, x} = \frac{\partial u}{\partial x} + \frac{1}{2} [{(\frac{\partial v}{\partial x})}^{2} + {(\frac{\partial w}{\partial x})}^{2}],

D = \{s + w_{c}\} i_{x} + v i_{y} + w i_{z},

where:

2

w_{c} = - \frac{1}{2} \int_{0}^{x} [{(\frac{\partial v}{\partial x})}^{2} + {(\frac{\partial v}{\partial x})}^{2}] d σ .

Fig. 1A uniform pinned-pinned axially loaded spinning square Rayleigh beam

The integral term in Eq. (2) reflects axial shortening resulting from the transverse deformation. The method neglecting the integral term is considered as a zero-order approximate method, the corresponding model is called a ZOAC model. Otherwise, it is considered as a first-order approximate method, and the corresponding model is called a FOAC model. A ZOAC model applies the small deformation theory and is valid at low spinning speeds, but may lead to erroneous results in some high speed cases.

The longitudinal strain based on the Euler-Bernoulli assumptions is as follows:

3

ε_{x x} = \frac{\partial u}{\partial x} - y \frac{\partial^{2} v}{\partial x^{2}} + \frac{1}{2} {(\frac{\partial v}{\partial x})}^{2} - z \frac{\partial^{2} w}{\partial x^{2}} + \frac{1}{2} {(\frac{\partial w}{\partial x})}^{2} .

Substituting Eq. (1) into Eq. (3), the longitudinal strain can be determined as:

4

ε_{x x} = s_{, x} - y \frac{\partial^{2} v}{\partial x^{2}} - z \frac{\partial^{2} w}{\partial x^{2}} .

Eq. (3) contains a nonlinear term of the stress-strain relation $(v_{, x}^{2} + w_{, x}^{2}) / 2$ , which is eliminated in Eq. (4), where $s_{, x}$ reflects axial strain, $y v_{, x x}$ and $z w_{, x x}$ reflect transverse bending strains. The substitution of $s$ with $u$ eliminates the nonlinear term. However, Eq. (4) still represents a nonlinear strain-displacement relation.

2.2. Formulation

The total kinetic energy and strain energy of a spinning flexible beam and the external virtual work acting on it [19] can be referred to Appendix A1. By applying Hamilton’s principle $\int_{0}^{t} δ T - δ E + δ W = 0$ , the nonlinear governing coupled PDEs of motion of the spinning Rayleigh beam are given as follows:

Force relation in $x$ direction:

5a

ρ A (\ddot{s} + {\ddot{w}}_{c}) - E A s_{, x x} = F_{x} δ (x - l_{1}) + e (x) ρ A Ω^{2} c o s θ_{1} s i n (Ω t + ϕ_{1} (x)) .

Force relation in $y$ direction:

5b

ρ A \ddot{v} - ρ I_{d} {\ddot{v}}_{, x x} - ρ I {}_{p} Ω {\dot{w}}_{, x x} + E I_{d} v_{, x x x x} + P v_{, x x} + ρ I_{d} Ω^{2} v_{, x x}

- ρ A \frac{\partial}{\partial x} (v_{, x} \int_{x}^{l} (\ddot{s} + {\ddot{w}}_{c}) d ζ) = F_{y} δ (x - l_{2}) + e (x) ρ A Ω^{2} c o s θ_{2} s i n (Ω t + ϕ_{2} (x)) .

Force relation in $z$ direction:

5c

ρ A \ddot{w} - ρ I_{d} {\ddot{w}}_{, x x} + ρ I {}_{p} Ω {\dot{v}}_{, x x} + E I_{d} w_{, x x x x} + P w_{, x x} + ρ I_{d} Ω^{2} w_{, x x}

- ρ A \frac{\partial}{\partial x} (w_{, x} \int_{x}^{l} (\ddot{s} + {\ddot{w}}_{c}) d ζ) = F_{z} δ (x - l_{3}) + e (x) ρ A Ω^{2} c o s θ_{3} s i n (Ω t + ϕ_{3} (x)) .

The coupling term of deformation $w_{c}$ is a second-order small quantity of transverse deformations $v$ and $w$ , and could be neglected. $e (x)$ is the mass eccentricity distribution of the beam. Therefore, if the second-order coupling terms of axial and transverse deformations are neglected, Eq. (5a)-(5c) may degenerate into the kinetic equations of a spinning Rayleigh beam using a ZOAC model [9].

It is almost impossible to derive an analytic solution to Eq. (5a)-(5c), which are non-linear and time variable with intensive coupling. To solve the equations, the discretization technique of continuous system should be used.

3. Finite element model

Spatial discretization is an important stage of dynamic modeling. Assumed mode method (AMM) and finite element method (FEM) are the most common methods in spatial discretization of flexible systems. Most researchers directly use the analysis results of the vibration frequencies and modals in structural dynamics, i.e. directly use the modes of structures with no global motion as the assumed modes for parts with global motion. While it is possible to get solutions to the problem using AMM, the accuracy of using static modes herein is questionable since large overall motion has an effect on modal frequencies. It is hard to have a quantitative estimate of the solution’s accuracy, even if it can get an accurate solution at small motion. Valembois et al. [20] showed that it appeared to be a very difficult task for AMM to deal with ﬂexible beams that are subject to a complex dynamic loading. However, FEM could get an accurate solution and had a convergence guarantee by the comparison among several discretization methods. Thus, this paper proposes a finite element method for a spinning Rayleigh beam, and compares with the solutions obtained from AMM.

3.1. The discretization of deformation field

A flexible beam is divided into $n$ elements, as shown in Fig. 2. Choose actual axial deformation $s$ and transverse displacements $v$ , $w$ in the former part as generalized coordinates. Therefore, the displacements of point $M$ on the ith element $s$ , $v$ and $w$ are described herein as the linear interpolation of the two nodes coordinates using element shape functions given by:

6

s (\bar{x}, t) = N_{1, i} (\bar{x}) η_{i} (t), v (\bar{x}, t) = N_{2, i} (\bar{x}) η_{i} (t), w (\bar{x}, t) = N_{3, i} (\bar{x}) η_{i} (t),

where $N_{1, i} (\bar{x})$ is the axial displacement shape function, $N_{2, i} (\bar{x})$ and $N_{3, i} (\bar{x})$ are the transverse displacement shape functions. $\bar{x}$ represents the longitudinal coordinate at element local coordinate, while $x$ represents the longitudinal coordinate at the structure coordinate.

Assume $ς = \bar{x} / l_{i}$ :

7

s (\bar{x}) = s_{i} + \frac{s_{i + 1} - s_{i}}{l_{i}} \bar{x} = (1 - ς) s_{i} + ς s_{i + 1},

8

v (\bar{x}) = f_{1} (\bar{x}) v_{i} + f_{2} (\bar{x}) Γ_{i} + f_{3} (\bar{x}) v_{i + 1} + f_{4} (\bar{x}) Γ_{i + 1},

where:

f_{1} (\bar{x}) = 1 - \frac{3 {\bar{x}}^{2}}{l_{i}^{2}} + \frac{2 {\bar{x}}^{3}}{l_{i}^{3}} = 1 - 3 ς^{2} + 2 ς^{3}, f_{2} (\bar{x}) = \bar{x} - \frac{2 {\bar{x}}^{2}}{l_{i}} + \frac{{\bar{x}}^{3}}{l_{i}^{2}} = l_{i} (1 - 2 ς^{2} + ς^{3}),

f_{3} (\bar{x}) = \frac{3 {\bar{x}}^{2}}{l_{i}^{2}} - \frac{2 {\bar{x}}^{3}}{l_{i}^{3}} = 3 ς^{2} - 2 ς^{3}, f_{4} (\bar{x}) = - \frac{{\bar{x}}^{2}}{l_{i}} + \frac{{\bar{x}}^{3}}{l_{i}^{2}} = l_{i} (- ς^{2} + ς^{3}),

and:

9

w (\bar{x}) = h_{1} (\bar{x}) w_{i} + h_{2} (\bar{x}) Β_{i} + h_{3} (\bar{x}) w_{i + 1} + h_{4} (\bar{x}) Β_{i + 1},

where:

h_{1} (\bar{x}) = 1 - \frac{3 {\bar{x}}^{2}}{l_{i}^{2}} + \frac{2 {\bar{x}}^{3}}{l_{i}^{3}} = 1 - 3 ς^{2} + 2 ς^{3}, h_{2} (\bar{x}) = \bar{x} - \frac{2 {\bar{x}}^{2}}{l_{i}} + \frac{{\bar{x}}^{3}}{l_{i}^{2}} = l_{i} (- 1 + 2 ς^{2} - ς^{3}),

h_{3} (\bar{x}) = \frac{3 {\bar{x}}^{2}}{l_{i}^{2}} - \frac{2 {\bar{x}}^{3}}{l_{i}^{3}} = 3 ς^{2} - 2 ς^{3}, h_{4} (\bar{x}) = - \frac{{\bar{x}}^{2}}{l_{i}} + \frac{{\bar{x}}^{3}}{l_{i}^{2}} = l_{i} (ς^{2} - ς^{3}) .

$η_{i} (t)$ is the displacement matrix of the element nodes:

10

η_{i} (t) = [\begin{array}{l} s_{i} & v_{i} & B_{i} & w_{i} & Γ_{i} \end{array} \begin{array}{l} s_{i + 1} & v_{i + 1} & Β_{i + 1} & w_{i + 1} & Γ_{i + 1} \end{array}]^{T},

where $Β_{i} = - \frac{\partial w}{\partial x} (0)$ , $Β_{i + 1} = - \frac{\partial w}{\partial x} (l_{i})$ , $Γ_{i} = \frac{\partial v}{\partial x} (0)$ , $Γ_{i + 1} = \frac{\partial v}{\partial x} (l_{i})$ .

Fig. 2The element division of the beam

The shape functions can be expressed as:

11

N_{1, i} (\bar{x}) = [\begin{array}{l} 1 - ς & 0 & 0 & 0 & 0 & ς & 0 & 0 & 0 & 0 \end{array}],

12

N_{2, i} (\bar{x}) = [\begin{matrix} 0 & 1 - 3 ς^{2} + 2 ς^{3} & 0 & 0 & l_{i} (ς - 2 ς^{2} + ς^{3}) & 0 & 3 ς^{2} - 2 ς^{3} & 0 & 0 & l_{i} (- ς^{2} + ς^{3}) \end{matrix}],

13

N_{2, i} (\bar{x}) = [\begin{matrix} 0 & 0 & l_{i} (- ς + 2 ς^{2} - ς^{3}) & 1 - 3 ς^{2} + 2 ς^{3} & 0 & 0 & 0 & l_{i} (ς^{2} - ς^{3}) & 3 ς^{2} - 2 ς^{3} & 0 \end{matrix}] .

Assume $η (t)$ as the displacement matrix of all nodes, therefore:

14

η_{i} (t) = Α_{i} η (t),

where $Α_{i}$ is Boole indicating matrix determined by element numbering. And:

15

η (t) = [\begin{array}{l} s_{1} & v_{1} & Β_{1} & w_{1} & Γ_{1} \end{array} \dots \begin{array}{l} s_{n + 1} & v_{n + 1} & Β_{n + 1} & w_{n + 1} & Γ_{n + 1} \end{array}]^{T},

16

Α_{i} = {[\begin{array}{l} 0 & 0 & \dots & I_{5 \times 5} & \dots & 0 & \dots & 0 \\ 0 & 0 & \dots & 0 & \dots & I_{5 \times 5} & \dots & 0 \end{array}]}_{10 \times 5 (n + 1)} .

Substituting Eq. (14) into Eq. (6), gives:

17

s (\bar{x}, t) = N_{1} (\bar{x}) η (t), v (\bar{x}, t) = N_{2} (\bar{x}) η (t), w (\bar{x}, t) = N_{3} (\bar{x}) η (t),

where $N_{1} (\bar{x}) = N_{1 . i} (\bar{x}) Α_{i}$ , $N_{2} (\bar{x}) = N_{2 . i} (\bar{x}) Α_{i}$ , $N_{3} (\bar{x}) = N_{3 . i} (\bar{x}) Α_{i}$ .

The displacements of an arbitrary point on the beam can be expressed as:

18

u_{p} = [\begin{array}{l} u \\ v \\ w \end{array}] = [\begin{matrix} N_{1} (\bar{x}) η (t) - \frac{1}{2} η^{T} (t) S (i, \bar{x}) η (t) \\ N_{2} (\bar{x}) η (t) \\ N_{3} (\bar{x}) η (t) \end{matrix}],

where coupling shape function:

19

S (i, \bar{x}) = (Α_{i})^{T} \int_{0}^{\tilde{x}} \frac{\partial {N_{2 . i}}^{T}}{\partial \bar{x}} \cdot \frac{\partial N_{2 . i}}{\partial \bar{x}} d \bar{x} Α_{i} + \sum_{j}^{i - 1} (Α_{j})^{T} \int_{0}^{l_{i}} \frac{\partial {N_{2 . j}}^{T}}{\partial \bar{x}} \cdot \frac{\partial N_{2, j}}{\partial \bar{x}} d \bar{x} Α_{i}

+ (Α_{i})^{T} \int_{0}^{\tilde{x}} \frac{\partial {N_{3 . i}}^{T}}{\partial \bar{x}} \cdot \frac{\partial N_{3 . i}}{\partial \bar{x}} d \bar{x} Α_{i} + \sum_{j}^{i - 1} (Α_{j})^{T} \int_{0}^{l_{i}} \frac{\partial {N_{3 . j}}^{T}}{\partial \bar{x}} \cdot \frac{\partial N_{3, j}}{\partial \bar{x}} d \bar{x} Α_{i} .

Eq. (19) describes a complete and second-order strain displacement field, and $S (i, x)$ is a symmetric and nonnegative definite matrix. Neglecting $S (i, x)$ , it becomes a linear displacement field in structure dynamics.

3.2. Element analysis

3.2.1. The dynamic equations of the element

Substituting Eq. (5), (16) and (17) into the Lagrange equation of the spinning flexible beam, the dynamic equations of spinning flexible beam element are deduced using Hamilton principle.

The total kinetic energy and strain energy of a spinning flexible beam element and the external virtual work acting on it are shown in Appendix A2. Using Hamilton’s principle, the dynamic equations of a spinning flexible beam element are as follows:

20

M_{i} \ddot{η} (t) + C_{i} \dot{η} (t) + K_{i} η (t) = {F_{i}}^{T},

where:

21

M_{i} = \int_{0}^{l_{i}} \{ρ A [N_{1}^{T} (\bar{x}) N_{1} (\bar{x}) + N_{2}^{T} (\bar{x}) N_{2} (\bar{x}) + N_{3}^{T} (\bar{x}) N_{3} (\bar{x})]

+ ρ I_{d} [{(\frac{\partial N_{3} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3} (\bar{x})}{\partial x}) + {(\frac{\partial N_{2} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2} (\bar{x})}{\partial x})]

- ρ A S (i, x) η (t) N_{1} (\bar{x})\} d x,

22

C_{i} = \int_{0}^{l_{i}} \{- 2 ρ I_{p} Ω [{(\frac{\partial N_{2} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3} (\bar{x})}{\partial x}) - {(\frac{\partial N_{3} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2} (\bar{x})}{\partial x})] - ρ A S (i, x) \dot{η} (t) N_{1} (\bar{x})} d x,

23

K_{i} = \int_{0}^{l_{i}} \{E A {(\frac{\partial N_{1} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{1} (\bar{x})}{\partial x}) + (E I_{d} - P) [{(\frac{\partial^{2} N_{2} (\bar{x})}{\partial x^{2}})}^{T} (\frac{\partial^{2} N_{2} (\bar{x})}{\partial x^{2}})

+ {(\frac{\partial^{2} N_{3} (\bar{x})}{\partial x^{2}})}^{T} (\frac{\partial^{2} N_{3} (\bar{x})}{\partial x^{2}})] + ρ I_{p} Ω^{2} \int_{0}^{l_{i}} [{(\frac{\partial N_{3} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3} (\bar{x})}{\partial x})

+ {(\frac{\partial N_{2} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2} (\bar{x})}{\partial x})] - ρ A S (i, x) \ddot{η} (t) N_{1} (\bar{x})\} d x,

24

F_{i} = \int_{0}^{l_{i}} {f_{x i} [N_{1} (\bar{x}) - η^{T} (t) S (i, x)] + f_{y i} N_{2} (\bar{x}) + f_{z i} N_{3} (\bar{x})} d x + ρ A Ω^{2} (K_{1} + K_{2} + K_{3}),

where:

Κ_{1} = - \int_{0}^{l_{i}} e (x) N_{1} (\bar{x}) s i n θ_{1} \sin (Ω t + ϕ_{1} (x)) d x,

Κ_{2} = - \int_{0}^{l_{i}} e (x) N_{2} (\bar{x}) s i n θ_{2} s i n (Ω t + ϕ_{2} (x)) d x

Κ_{3} = - \int_{0}^{l_{i}} e (x) N_{3} (\bar{x}) s i n θ_{3} s i n (Ω t + ϕ_{3} (x)) d x

3.2.2. Element assembling

The finite element equations of the spinning flexible beam are derived by assembling the former mass matrix, damping matrix, stiffness matrix and force matrix of the beam element according to the element division. Substituting Eq. (15), (17) and (18) into Eq. (21)-(24), the mass matrix, damping matrix, stiffness matrix and force matrix of the beam are:

25

M = \sum_{i = 1}^{n} \int_{0}^{l_{i}} \{ρ A (Α_{i})^{T} [N_{1, i}^{T} (\bar{x}) N_{1, i} (\bar{x}) + N_{2, i}^{T} (\bar{x}) N_{2, i} (\bar{x}) + N_{3, i}^{T} (\bar{x}) N_{3, i} (\bar{x})] Α_{i}

+ ρ I_{d} (Α_{i})^{T} [{(\frac{\partial N_{3, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3, i} (\bar{x})}{\partial x}) + {(\frac{\partial N_{2, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2, i} (\bar{x})}{\partial x})] Α_{i}

- ρ A S (i, x) η (t) N_{1, i} (\bar{x}) Α_{i}} d x,

26

C = \sum_{i = 1}^{n} \int_{0}^{l_{i}} \{- 2 ρ I_{p} Ω (Α_{i})^{T} [{(\frac{\partial N_{2, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3, i} (\bar{x})}{\partial x}) - {(\frac{\partial N_{3, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2, i} (\bar{x})}{\partial x})] Α_{i}

- ρ A S (i, x) \dot{η} (t) N_{1, i} (\bar{x}) Α_{i}} d x,

27

K = \sum_{i = 1}^{n} \int_{0}^{l_{i}} \{E A (Α_{i})^{T} {(\frac{\partial N_{1, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{1, i} (\bar{x})}{\partial x}) Α_{i}

+ (E I_{d} - P) (Α_{i})^{T} [{(\frac{\partial^{2} N_{2, i} (\bar{x})}{\partial x^{2}})}^{T} (\frac{\partial^{2} N_{2, i} (\bar{x})}{\partial x^{2}}) + {(\frac{\partial^{2} N_{3, i} (\bar{x})}{\partial x^{2}})}^{T} (\frac{\partial^{2} N_{3, i} (\bar{x})}{\partial x^{2}})] Α_{i}

+ ρ I_{p} Ω^{2} (Α_{i})^{T} [{(\frac{\partial N_{3, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{3, i} (\bar{x})}{\partial x}) + {(\frac{\partial N_{2, i} (\bar{x})}{\partial x})}^{T} (\frac{\partial N_{2, i} (\bar{x})}{\partial x})] Α_{i}

- ρ A S (i, x) \ddot{η} (t) N_{1, i} (\bar{x}) Α_{i}} d x,

28

F = \sum_{i = 1}^{n} \int_{0}^{l_{i}} {f_{x i} [N_{1, i} (\bar{x}) Α_{i} - S (i, x) η (t)] + f_{y i} N_{2, i} (\bar{x}) Α_{i} + f_{z i} N_{3, i} (\bar{x}) Α_{i}} d x .

If the force is independent on spatial coordinates, then:

29

F = \sum_{i = 1}^{n} \{f_{x i} \int_{0}^{l_{i}} [N_{1, i} (\bar{x}) Α_{i} - S (i, x) η (t)] + f_{y i} \int_{0}^{l_{i}} N_{2, i} (\bar{x}) d x Α_{i} + f_{z i} \int_{0}^{l_{i}} N_{3, i} (\bar{x}) d x Α_{i}

+ ρ A Ω^{2} (Κ_{1}^{'} + Κ_{2}^{'} + Κ_{3}^{'})} d x,

where:

Κ_{1}^{'} = - \int_{0}^{l_{i}} e (x) N_{1} (\bar{x}) s i n θ_{1} \sin (Ω t + ϕ_{1} (x)) d x A_{i},

Κ_{2}^{'} = - \int_{0}^{l_{i}} e (x) N_{2} (\bar{x}) s i n θ_{2} \sin (Ω t + ϕ_{2} (x)) d x A_{i},

Κ_{3}^{'} = - \int_{0}^{l_{i}} e (x) N_{3} (\bar{x}) s i n θ_{3} \sin (Ω t + ϕ_{3} (x)) d x A_{i} .

It is important to note that the mass, damping, stiffness and force matrices of the beam based on an FOAC model are different from the original dynamic matrices in which they are not only related to the properties of the beam but also to the real-time dynamic responses.

4. Piezoelectric controller implementation

Since the large displacements caused by the nonlinear vibration are dangerous to the stability and lifespan of structures, suppressing such vibrations are essential. Piezoelectric actuators are widely used in active vibration control because of their low-weight, large bandwidth and fast response speed. Thus, piezoelectric controllers are implemented in this paper to suppress structural vibration.

4.1. Dynamic equation of piezoelectric element

The piezoelectric beam and piezoelectric element are shown as Fig. 3. The width and thickness of the piezoelectric element are $b$ and $h_{p}$ , respectively, and the height of the base is $h_{1}$ .

Fig. 3Piezoelectric beam and piezoelectric element

The strain energy of the beam includes the electrical coenergy, and Eq. (A7) can be expressed as:

30

E_{j} = \frac{1}{2} \int_{0}^{l} [{E A s}_{j, x}^{2} + E I_{d} (v_{j, x x}^{2} + w_{i, x x}^{2}) - P (v_{j, x}^{2} + w_{j, x}^{2})] d x - \frac{1}{2} \int_{0}^{l_{j}} \int_{A} E^{T} D d A d x,

where $E_{u}$ and $D$ reflect electric field and electric displacement, respectively. Both the electric field and the electric displacement vectors only contain a single non-zero component, respectively $E_{3}$ and $D {}_{3}$ . The linear constitutive relationships describing the electrical and mechanical interaction for piezoelectric materials are therefore given as:

31

\{\begin{array}{l} σ_{1} = E_{P} ε_{1} - e_{31} E_{3}, \\ D = e_{31} ε_{1} + ε_{33}^{ε} E_{3}, \end{array}

where $E_{P}$ is Young modulus of the piezoelectric material along its longitudinal axis.

The transverse displacements of the beam element is:

32

v (\bar{x}, t) = N_{2} (\bar{x}) η (t) = N_{2 . j} (\bar{x}) Α_{j} η (t) .

The longitudinal strain of the piezoelectric element based on the plane section assumption can be expressed as:

33

ε_{1} = \frac{h_{p}}{2} \frac{\partial^{2}}{\partial x^{2}} v (x, t) = \frac{h_{p}}{2} V_{j} (x) Α_{j} η (t),

where $V_{j} (x) = \frac{\partial^{2}}{\partial x^{2}} N_{2 . j} (\bar{x})$ .

The electric field through the piezoelectric material is given as:

34

E_{3} = \frac{U}{(h_{p} - h_{1}) / 2},

where $U$ is the voltage applied on the surfaces of the piezoelectric material.

Define the generalized strain of the beam element as the curvature of the beam, and the generalized stress as the section moment:

35

κ_{x} = \frac{\partial^{2}}{\partial x^{2}} v (x, t) = V_{j} (x) A_{j} η (t),

36

M_{p} = E I_{d} \frac{\partial^{2}}{\partial x^{2}} v (x, t) = E I_{d} κ_{x} = E I_{d} V_{j} (x) A_{j} η (t),

where $E I_{d}$ is the elastic matrix corresponding to the generalized strain and the generalized stress.

The dynamic equations of the piezoelectric element are derived:

37

δ \int_{t_{1}}^{t_{2}} L d t = δ \int_{t_{1}}^{t_{2}} (T - U_{e} + W) d t = 0 .

The kinetic energy of the piezoelectric element is given as:

38

T_{j} = \frac{ρ'}{2} \int_{0}^{l} \{A ({\dot{u}}_{j}^{2} + {\dot{v}}_{j}^{2} + {\dot{w}}_{j}^{2}) + I_{d} (Β_{j, t}^{2} + Γ_{j, t}^{2})

+ I_{p} [Ω (Γ_{j} B_{j, t} - B_{j} Γ_{j, t}) + \frac{1}{2} Ω^{2} (Β_{j}^{2} + Γ_{j}^{2}) + Ω^{2}]\} d x,

where $ρ'$ is the density of the piezoelectric element.

The strain energy and the external virtual work are, respectively:

39

{U_{e}}_{j} = \frac{1}{2} \int_{0}^{l} [{E A s_{j}}_{, x}^{2} + E I_{d} (v_{j, x x}^{2} + w_{j, x x}^{2}) - P (v_{j, x}^{2} + w_{j, x}^{2})] d x

- \frac{e_{31} h_{p} b U}{4} {\bar{V}}_{j} (x) A_{j} η (t) - \frac{ε_{33}^{ε} b l_{j}}{h_{p} - h_{1}} U^{2},

40

W_{j} = f_{x j} u_{j} + f_{y j} v_{j} + f_{z j} w_{j} - U Q,

where ${\bar{V}}_{j} (x)$ is the integration of $V_{j} (x)$ with respect to spatial coordinate $x$ , and $Q$ is the free charge on surfaces of single piezoelectric actuator:

41

{\bar{V}}_{j} (x) = \int_{0}^{l_{j}} V_{j} (x) d x,

42

Q = \int_{s} D_{3} d s,

where $s$ means the section of the piezoelectric material.

The dynamic equations of the piezoelectric element are:

43

\{\begin{array}{l} M_{j} \ddot{η} + C_{j} \dot{η} + K_{j} η = F_{j} + K_{d} U, \\ \frac{e_{31} h_{p} b}{4} {\bar{W}}_{j} (x) A_{j} η (t) + 2 \frac{ε_{33}^{ε} b l_{j}}{h_{p} - h_{1}} U = Q, \end{array}

where $K_{d} = (e_{31} h_{p} b U / 4) ({\bar{V}}_{j} (x) A_{j})^{T}$ .

4.2. Piezoelectric sensor analysis

The piezoelectric sensors generate electric charges due to piezoelectric effect when the beam experiences structural deformations. The electric field intensity can be expressed as:

44

E_{3} = \frac{Q}{ε_{33}^{ε} s} .

Let $U$ in Eq. (43) be zero and substitute Eq. (43) into Eq. (44). Then the electric filed intensity of the piezoelectric sensors under only external force can be given as:

45

E_{3} = \frac{e_{31} h_{p} b}{4 ε_{33}^{ε} s} {\bar{W}}_{j} (x) A_{j} η (t) .

The current and the voltage in the amplifying circuit satisfy that:

46

i = \frac{d Q}{d t}, V_{s} = R i,

where $R$ is the resistance, and $V_{s}$ is the output voltage of amplifying circuit.

The output voltage $V_{s}$ can also be derived from Eq. (43):

47

V_{s} = R \frac{e_{31} h_{p} b}{8} {\bar{V}}_{j} (x) A_{j} \dot{η} (t) = K_{s} \dot{η} (t) .

Eq. (47) illustrates the relation between the output voltage of the sensor and the velocity of the node, and $K_{s}$ is the coefficient matrix of the output voltage:

48

K_{s} = R \frac{e_{31} h_{p} b}{8} {\bar{V}}_{j} (x) A_{j} .

4.3. Velocity feedback control

A velocity feedback controller is implemented if sensors and actuators are collocated and a negative proportional feedback controller is used:

49

K_{s} = \frac{R}{2} K_{d}^{T},

50

U = - G V_{s} = - G K_{s} \dot{η} (t) .

Substituting the former equation into Eq. (43), the finite element equations of the closed loop system can be expressed as:

51

M_{j} \ddot{η} + (C_{j} + K_{d} G K_{s}) \dot{η} + K_{j} η = F_{j} .

The use of collocated actuator/sensor pairs and the introduction of negative proportion velocity feedback method result in an additional proportion damping to the original system. This is called active damping and requires relatively little control effort. Provided that the actuator and sensor are collocated, active damping has perfect dynamic characteristics and can achieve guaranteed stability without a model of the structure.

5. Numerical example

Consider a spinning square beam whose physical parameters are as follows: the density of the beam $ρ$ is 2000 Kg/m³, the area of section $A$ is 4×10^-6 m², the length $l$ is 0.6 m, the moment of inertia $I_{d}$ is $1$ .33×10^-12 m⁴, the polar moment of inertia $I_{p}$ is 2.66×10^-12 m⁴, elastic modulus $E$ is 1×10¹⁰ N/m², the axial force $P$ is 0 N. Assume the three-directional centrifugal forces by the eccentricity are:

F_{e x} = ρ A Ω^{2} K_{1} = 0.06 ρ A Ω^{2} s i n (Ω t + 0.5 π),

F_{e y} = ρ A Ω^{2} K_{2} = 0.001 ρ A Ω^{2} s i n (Ω t + 0.1 π),

F_{e z} = ρ A Ω^{2} K_{3} = 0.001 ρ A Ω^{2} s i n (Ω t + 0.1 π) .

The external forces are as shown in Fig. 4, where $F_{y} = 0.05 s i n (6 π t)$ , $F_{x} = F_{z} = 0$ .

Fig. 4The external force acting on the structure

Two pairs of piezoelectric actuators and sensors are collocated at the 7th and 15th element in $y$ direction, respectively. The feedback gain $G$ is set to 6×10⁵, and the typical properties of piezoelectric material are as follows, $e_{31} =$ –7.5, $h_{p} = 0.0005,$ $b = 0.002,$ the voltage $U$ is 100 V, $R$ is 2×10⁴ Ω.

5.1. Dynamic responses of ZOAC and FOAC model

Fig. 5 shows the transverse responses of a beam spinning at 150 rad/s derived by AMM and FEM, respectively. The difference between the displacements derived by two different methods is obvious. The main reason is that AMM directly uses the analysis results of the vibration frequencies and modals in structural dynamics, which neglects the effect of a large overall motion. On the other hand, FEM allows a rigid multibody formalism to study flexible slender bodies directly [20]. Therefore, FEM could obtain more accurate dynamic responses of spinning flexible beams.

Fig. 6 illustrates the transverse responses of a spinning flexible beam at a speed of 50 rad/s based on ZOAC and FOAC model. The differences are so small since the effect of the rigid-flexible coupling terms is small. As a result, the ZOAC model is a valid dynamic model of a spinning flexible beam at low spinning speeds.

At 400 rad/s the differences in transverse responses of a spinning flexible beam predicted by the FOAC and ZOAC models are more apparent (Fig. 7). The differences are large due to the effect of the rigid-flexible coupling terms increases as spinning speed and cannot be neglected.

Fig. 5Dynamic responses of a spinning beam using assumed mode method (AMM) and finite element method (FEM), respectively, at a spinning speed 150 rad/s

Fig. 6Dynamic responses of a spinning beam using finite element method based on FOAC model and ZOAC model, respectively, at a spinning speed 50 rad/s

Dynamic responses of a spinning beam using finite element method based on FOAC model and ZOAC model, respectively, at a spinning speed 50 rad/s

Fig. 7Dynamic responses of a spinning beam using finite element method based on FOAC model and ZOAC model, respectively, at a spinning speed 400 rad/s

Dynamic responses of a spinning beam using finite element method based on FOAC model and ZOAC model, respectively, at a spinning speed 400 rad/s

The differences of transverse responses derived by a ZOAC and an FOAC model are investigated. Define values as following that reflect the similarity of the responses:

52

\{\begin{array}{l} R V = \frac{{|v_{0}^{T} v_{1}|}^{2}}{(v_{0}^{T} v_{0}) (v_{1}^{T} v_{1})}, \\ R W = \frac{{|w_{0}^{T} w_{1}|}^{2}}{(w_{0}^{T} w_{0}) (w_{1}^{T} w_{1})}, \end{array}

where $v_{0}$ , $w_{0}$ represent the displacements of an arbitrary point derived by a ZOAC model, and $v_{1}$ , $w_{1}$ represent the displacements of the same point derived by an FOAC model. If $R V$ and $R W$ are 1, it implies that the displacements derived by a ZOAC model and an FOAC model are the same. Corresponding values of 0 mean that the displacements are totally different.

Fig. 8 shows that the difference between the responses derived by the two models gets larger as the spinning speeds increase. The similarity decreases sharply at a speed over 280 rad/s, to 0.9 at about 320 rad/s, meaning that the difference is so large that the ZOAC model is not accurate any more. It can be seen that at higher spinning speeds, accounting for the effect of the rigid-flexible coupling is a key factor to obtain a more accurate dynamic response.

Fig. 8Similarity of dynamic responses of a spinning beam using finite element method based on FOAC model at different spinning speeds

From the above discussions, the dynamic matrices are influenced by the effect of the rigid-flexible coupling terms, which increases as the spinning speeds increase. The effect is small and has little influence on the dynamic matrices at low speeds. Thus the dynamic responses derived by two different models are almost the same. In these low speed cases, the ZOAC model is valid. However, in high speed cases, the dynamic matrices change a lot since the coupling effect gets lager. The difference of the dynamic responses is obvious and the ZOAC model is no longer accurate.

Fig. 9Dynamic responses of a spinning beam without and with control at a spinning speed 50 rad/s

5.2. Vibration control

Numerical simulation of active vibration control of the spinning flexible beam in Figs. 9-10. It shows that the piezoelectric controllers increase system damping and result in an obvious vibration suppression of the spinning beam. The displacement and the velocity in $y$ direction are suppressed markedly, and the peak of dynamic response could be reduced about 78 %. There is no effect on response in $z$ direction due to the fact that the piezoelectric actuators only suppress the vibration perpendicular to the planes of piezoelectric films. Besides, since the beam is forced by centrifugal forces caused by eccentricity all the time, the transverse responses don’t decay asymptotically to zero.

Fig. 10Displacement frequency response a) and velocity frequency response b) without and with control

a)

b)

6. Conclusions

In this paper, finite element modeling and velocity feedback control strategy for the active vibration control of a high-speed spinning flexible mechanical-electric coupled beam are developed. The dynamic equations of a pinned-pinned spinning flexible Rayleigh beam accounting for the effect of the rigid-flexible coupling are derived in detail on the basis of a first order approximate coupling (FOAC) model, and a finite element method is used to obtain the dynamic responses. The numerical simulation of a spinning beam at different speeds concludes that a zero order approximate coupling (ZOAC) model is valid to be used as the dynamic model of a spinning flexible beam at low speeds, but invalid at high speeds. This problem is largely overcome by use of the FOAC model, which accounts for the effect of the rigid-flexible coupling. In addition, the vibration of the spinning beam can be suppressed effectively by the proposed velocity feedback controller.

References

Bauer Helmut F. Vibration of a rotating uniform beam, part I: orientation in the axis of rotation. Journal of Sound and Vibration, Vol. 72, Issue 2, 1980, p. 177-189.

Search CrossRef
Pai P. Frank, Qian Xin, Du Xingwen Modeling and dynamic characteristics of spinning Rayleigh beams. International Journal of Mechanical Sciences, Vol. 68, 2013, p. 291-303.

Search CrossRef
Dimentberg Fedor M. Flexural Vibrations of Rotating Shafts. Butterworths, London, 1961.

Search CrossRef
Bauer H. F., Eidel W. Vibration of a rotating uniform beam, part II: orientation perpendicular to the axis of rotation. Journal of Sound and Vibration, Vol. 122, Issue 2, 1988, p. 357-375.

Search CrossRef
Pai P. Frank Highly Flexible Structures: Modeling, Computation, and Experimentation. AIAA Education Series, 2007.

Search CrossRef
Irving H. Shames Energy and Finite Element Methods in Structural Mechanics. CRC Press, 1985.

Search CrossRef
Muszynska Agnieszka Rotordynamics. CRC Press, 2010.

Search CrossRef
Yamamoto Toshio, Yukio Ishida Linear and Nonlinear Rotordynamics: a Modern Treatment with Applications. John Wiley & Sons, Vol. 11, 2001.

Search CrossRef
Sheu G. J., Shih-Ming Yang Dynamic analysis of a spinning Rayleigh beam. International Journal of Mechanical Sciences, Vol. 47, Issue 2, 2005, p. 157-169.

Search CrossRef
Qian Xin, Guo Chang Lin, Xing Wen Du Vibration measurement and data analysis of a spinning shaft using a camera-based motion analysis system. Applied Mechanics and Materials, Vol. 29, 2010, p. 203-208.

Search CrossRef
Boukhalfa A. Dynamic analysis of a spinning laminated composite-material shaft using the hp-version of the finite element method. Advances in Vibration Analysis Research, 2011, p. 161-186.

Search CrossRef
Lee U., Jang I. J., Park I. W. Frequency-domain spectral element model of a uniform spinning shaft. Applied Mechanics and Materials, Vol. 224, 2012, p. 264-267.

Search CrossRef
Mustapha K. B., Zhong Z. W. Spectral element analysis of a non-classical model of a spinning micro beam embedded in an elastic medium. Mechanism and Machine Theory, Vol. 53, 2012, p. 66-85.

Search CrossRef
Chan K. T., et al. Revolving superposed standing waves in a spinning Timoshenko beam. Journal of Sound and Vibration, Vol. 331, Issue 4, 2012, p. 815-832.

Search CrossRef
Khadem S. E., Shahgholi M., Hosseini S. A. A. Two-mode combination resonances of an in-extensional rotating shaft with large amplitude. Nonlinear Dynamics, Vol. 65, Issue 3, 2011, p. 217-233.

Search CrossRef
Shahgholi Majid, Khadem Siamak Esmaeilzadeh Resonance analysis of gyroscopic nonlinear spinning shafts with parametric excitations and speed fluctuations. International Journal of Mechanical Sciences, Vol. 64, Issue 1, 2012, p. 94-109.

Search CrossRef
Shahgholi Majid, Khadem Siamak Esmaeilzadeh Resonances of an in-extensional asymmetrical spinning shaft with speed fluctuations. Meccanica, Vol. 48, Issue 1, 2013, p. 103-120.

Search CrossRef
Yoo H. H., Ryan R. R., Scott R. A. Dynamics of flexible beams undergoing overall motions. Journal of Sound and Vibration, Vol. 181, Issue 2, 1995, p. 261-278.

Search CrossRef
Mamandi Ahmad, Kargarnovin Mohammad H. Nonlinear dynamic analysis of an axially loaded rotating Timoshenko beam with extensional condition included subjected to general type of force moving along the beam length. Journal of Vibration and Control, Vol. 19, Issue 16, 2013, p. 2448-2458.

Search CrossRef
Valembois R. E., Fisette Paul, Samin Jean-Claude Comparison of various techniques for modelling flexible beams in multibody dynamics. Nonlinear Dynamics, Vol. 12, Issue 4, 1997, p. 367-397.

Search CrossRef

About this article

Received

06 January 2015

Accepted

05 May 2015

Published

30 September 2015

SUBJECTS

Vibration generation and control

Keywords

high-speed

spinning flexible beam

first-order approximation model

finite element method

active control

Acknowledgements

This work was supported by: (1) Funding of Jiangsu Innovation Program for Graduate Education (Grant No. CXZZ13_0148) and the Fundamental Research Funds for the Central Universities, (2) a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.