Active vibration control of plate by active mass damper and negative acceleration feedback control algorithms

Abstract

In this study, an active mass damper that could track desired displacement was considered as a controller to suppress vibrations of a plate. Multi-input multi-output modal-space negative acceleration feedback control algorithm and decentralized multi-input multi-output negative acceleration feedback control algorithm were proposed and their stability conditions were checked theoretically. Stability conditions for both controllers are static, and they can effectively increase the damping of target mode without causing instability in low-frequency region. It was proved theoretically that both the multi-input multi-output modal-space negative acceleration feedback control and the decentralized multi-input multi-output negative acceleration feedback control can suppress vibrations of a continuous system such as a plate.

Keywords

Active mass damper vibration control of plate multi-input multi-output modal-space negative acceleration feedback decentralized multi-input multi-output negative acceleration control

Introduction

Reduction of vibrations in structures is a main engineering objective for humans, machinery, sensitive instruments, automobiles, airplanes and buildings because excessive vibrations may harm human body and cause structural failure. In some cases, structural modification may be impossible. For such cases, an auxiliary system that can absorb vibrations of the main structure is preferred. A dynamic absorber (DA) or a tuned mass damper (TMD) that can be easily attached to the structure of interest has been developed to suppress vibration of the primary structure. DA or TMD is a passive system that suppresses vibration of the primary structure by tuning its natural frequency to dominant frequency component. Dominant frequency can be either excitation frequency or natural frequency of the primary structure.¹ DA has been used to absorb vibration caused by harmonic disturbance with fixed frequency. If damping of the DA is zero, then vibration can be suppressed theoretically. However, a small amount of damping is more effective when disturbance has varying excitation frequency. Hence, we will use the term of TMD instead of DA in this case. Passive TMDs have been effectively used to suppress vibrations of structures. However, TMD has limited performance because of its fixed structural parameters, a narrow suppression frequency range, ineffective reduction of non-stationary vibration, and a sensitivity problem because of detuning. To improve the performance of TMD, a new light-weight passive control device called Tuned Mass Damper Inerter (TMDI) has been proposed to achieve mass amplification effect.^2,3 Cao and Li⁴ have proposed Tuned Tandem Mass Damper-Inerters by integrating tuned tandem mass dampers with two inerters, resulting in high effectiveness and wide frequency spacing.

To suppress vibrations and improve the performance of system of interest, active methods have been sought. Semi-active mass damper (SMD) and active mass damper (AMD) have been proposed in order to overcome limitations of passive TMD. Both SMD and AMD are equipped with sensors, actuators, and a microprocessor that can operate according to a control algorithm. SMD differs from AMD in that the damping coefficient of the damper is changed by control command, while AMD produces a force according to command signal. One benefit of using SMD or AMD is that it can cope with changes in the dominant frequency such as changes in natural frequency of the structure or excitation frequency of external disturbances by sensing vibrations. However, AMD may destabilize a main structure if it is ill-designed. Hence, reliability of the control system should be guaranteed before implementation.

Control algorithm is very important when applying AMD technology to real structures because instability may occur. Many control algorithms have been proposed and validated both theoretically and experimentally. Kobori et al.^5,6 have proposed a simple control algorithm based on optimal control theory. Li et al.⁷ have applied multiple AMDs to torsional vibration of building structure. Li et al.^8,9 have also used multiple AMD to attenuate vibrations of structures due to ground excitation. They demonstrated that multiple ATMD could significantly reduce vibration of structures. Chung et al.¹⁰ have developed active tendon system driven by hydraulic servo system with optimal control law. Chang and Yang¹¹ have proposed a feedback control using displacement, velocity, and acceleration measurements. Ankireddi and Yang¹² have also investigated feedback control for AMD. Cao et al.¹³ have obtained linear quadratic regulator (LQR) and a nonlinear feedback control algorithm. Cao and Li¹⁴ have developed a robust control algorithm. Samali and Al-Dawod¹⁵ have applied a Fuzzy Logic controller to AMD. Ricciardelli et al.¹⁶ have considered AMD based on linear quadratic regulator theory. Cao and Lin¹⁷ have proposed new control strategies. Wang and Lin¹⁸ have developed variable structure control and fuzzy sliding-mode control.

Research on control algorithms for AMD mentioned above assumed that displacement, velocity, and acceleration were measurable. However, a popular sensor for vibration measurement in real world is an accelerometer. One may think that velocity and displacement can be obtained by integrating accelerometer signal in time. However, this is not possible because of drift and bias in the accelerometer signal. Observer can be constructed for velocity and displacement using measured acceleration. This would require additional computation. Hence, a feedback control algorithm that can use acceleration signal directly is preferred to avoid computational complexities for real-world application.

Recently, negative acceleration feedback (NAF) control has been proposed by Yang et al.¹⁹ to control moving mass based on servo mechanism. The advantage of a servo control system is that movement of the mass can be accurately controlled by the AC servo motor. Yang et al.¹⁹ have proposed an NAF control that uses an accelerometer signal and produces desired displacement of the moving mass. The proposed NAF control is different from other control algorithms that produce control force. The NAF control theory that directly uses an accelerometer signal as a sensor input is different from other control theories that require displacement or velocity. Yang et al.¹⁹ have demonstrated the effectiveness of NAF control experimentally by using a building-like structure. Talib et al.²⁰ have designed multi-input multi-output (MIMO) modal-space NAF control and verified the control law experimentally by using a building-like structure. They used linear servo motors to move masses.

Although AMD with NAF control has been successfully applied to building-like structure, whether AMD can be applied to a continuous structure such as plate remains unclear. Active vibration control of a plate-like structure has been studied by many researchers. However, most of these studies were focused on vibration control by means of piezoelectric wafers. Lazarus et al.²¹ have carried out theoretical modelling and MIMO LQG control for a rectangular plate by using an assumed mode method. Li et al.²² have investigated the use of a piezoelectric wafer attached to a rectangular plate using $μ$ control theory. Kwak and Yang²³ have developed a dynamic model for a submerged rectangular plate equipped with piezoelectric sensors and actuators and designed the active vibration control. The actuating mechanism of AMD is different from the one of a piezoelectric actuator. AMD is an inertial actuator because the force is generated by a moving mass. On the other hand, the piezoelectric actuator is a strain actuator because it deforms when it is subjected to applied voltage. Current application of a piezoelectric actuator is limited to light structures such as thin plate and thin shell because its wafer dimension is small and expensive. In this regard, AMD is more practical than piezoelectric wafer in controlling vibrations of a large plate-like structure.

In this study, active vibration control of a rectangular plate equipped with multiple AMDs was investigated theoretically. In order to increase damping of natural modes, MIMO modal-space NAF control and decentralized MIMO NAF control were proposed for vibration suppression of a plate, and their stability conditions were investigated mathematically. It was proven that stability conditions of the MIMO modal-space NAF control and the decentralized MIMO NAF control are static, meaning that the stability of control algorithms will not depend on frequency. Control and observer spillover problems that might occur when applying MIMO modal-space NAF control were also investigated. The proposed MIMO modal-space NAF control and the decentralized MIMO NAF control were tested numerically. Numerical results showed that the proposed control methods can be successfully applied to a plate structure if the gain matrix is sufficiently small.

Dynamic modelling of rectangular plate with AMDs

We will briefly summarize the work of Kwak and Han²⁴ on free vibration analysis of a rectangular plate. Let us consider a rectangular plate with side lengths $a$ in the $x$ direction and $b$ in the $y$ direction and that it is equipped with AMDs as shown in Figure 1.

Figure 1.

Rectangular plate with AMD.

The assumed mode method is used to express the deflection of the plate as

w (x, y, t) = W (x, y) q (t)

(1)

where

W (x, y) = [W_{1} (x, y) W_{2} (x, y) \dots W_{m} (x, y)]

is a

1 \times m

matrix consisting of admissible functions, and

q (t) = {[q_{1} q_{2} \dots q_{m}]}^{T}

is an

m \times 1

vector consisting of generalized coordinates, for which

m

is the number of admissible functions used to approximate the deflection. In order to easily calculate mass and stiffness matrices, the admissible function matrix of equation (1) can be expressed in terms of admissible function matrices in each direction.²⁵

W (x, y) = Φ (ξ) \otimes Ψ (ζ)

(2)

where

\otimes

indicates element by element multiplication and non-dimensional variables,

ξ = x / a

and

ζ = y / b

are introduced for each coordinate, and

Φ (ξ) = [φ_{1} φ_{2} \dots φ_{m}]

and

Ψ (η) = [ψ_{1} ψ_{2} \dots ψ_{m}]

are

1 \times m

matrices.

Kinetic and potential energies of the rectangular plate can then be expressed as follows²⁵

T_{P} = \frac{ρ_{p} h ab}{2} {\dot{q}}^{T} {\bar{M}}_{p} \dot{q}

(3a)

V_{P} = \frac{Db}{2 a^{3}} q^{T} {\bar{K}}_{p} q

(3b)

where

ρ_{p}

is the mass density,

h

is the thickness,

D = E h^{3} / 12 (1 - v^{2})

E

is the Young’s modulus, and

ν

is Poisson’s ratio.

{\bar{M}}_{p}

and

{\bar{K}}_{p}

are non-dimensionalized mass and stiffness matrices, respectively. Subscript

p

represents the plate. The non-dimensionalized mass and stiffness matrices can be further expressed as follows²⁴

{\bar{M}}_{p} = X \otimes Y

(4a)

{\bar{K}}_{p} = \hat{X} \otimes Y + α^{4} X \otimes \hat{Y} + α^{2} ν ({\tilde{X}}^{T} \otimes \tilde{Y} + \tilde{X} \otimes {\tilde{Y}}^{T}) + α^{2} (1 - ν) \bar{X} \otimes \bar{Y}

(4b)

where

α = a / b

represents the aspect ratio of the plate and

X = \int_{0}^{1} Φ^{T} Φ d ξ

(5a)

\bar{X} = \int_{0}^{1} {Φ^{'}}^{T} Φ^{'} d ξ

(5b)

\hat{X} = \int_{0}^{1} {Φ^{″}}^{T} Φ^{″} d ξ

(5c)

\tilde{X} = \int_{0}^{1} Φ^{T} Φ^{″} d ξ

(5d)

Y = \int_{0}^{1} Ψ^{T} Ψ d η

(5e)

\bar{Y} = \int_{0}^{1} {Ψ^{'}}^{T} Ψ^{'} d η

(5f)

\hat{Y} = \int_{0}^{1} {Ψ^{″}}^{T} Ψ^{″} d η

(5g)

\tilde{Y} = \int_{0}^{1} Ψ^{T} Ψ^{″} d η

(5h)

where

' = d / d ξ, ″ = d^{2} / d ξ^{2} .

If n admissible functions are used in the x and y directions and the combination of these admissible functions is used, a total of

n^{2}

admissible functions can be obtained, yielding

m = n^{2}

. If each type of admissible function is considered to be

χ_{i} (i = 1, 2, \dots, n)

and

γ_{i} (i = 1, 2, \dots, n)

, then the relationship between the sequence of admissible functions introduced in equation (2) and that of separate admissible functions can be expressed as below²⁴

φ_{k} = {\begin{array}{l} χ_{1} 1 \leq k \leq n \\ χ_{2} n + 1 \leq k \leq 2 n \\ χ_{3} 2 n + 1 \leq k \leq 3 n \\ ⋮ \\ χ_{n} (n - 1) n + 1 \leq k \leq n^{2} \end{array}

(6a)

ψ_{k} = {\begin{array}{l} γ_{k} 1 \leq k \leq n \\ γ_{k - n} n + 1 \leq k \leq 2 n \\ γ_{k - 2 n} 2 n + 1 \leq k \leq 3 n \\ ⋮ \\ γ_{k - (n - 1) n} (n - 1) n + 1 \leq k \leq n^{2} \end{array}

(6b)

Therefore, instead of integrating $m^{2} = n^{4}$ elements in equation (5), $n^{2}$ integration and matrix rearrangement will suffice.²⁴

For an all simply-supported plate, the eigenfunction of a simply-supported beam²⁵ can be used.

χ_{i} = γ_{i} = \sin i π ξ i = 1, 2, \dots, n

(7)

Detailed expressions of matrices given in equation (5) have been described by Kwak and Yang.²⁶

Kinetic energy of moving and supporting masses of the AMD located at $(x_{ak}, y_{ak})$ can be expressed as follows

T_{ak} = \frac{1}{2} [m_{s} {\dot{w}}^{2} (x_{ak}, y_{ak}) + m_{a} {(\dot{w} (x_{ak}, y_{ak}) + {\dot{u}}_{ak})}^{2}], k = 1, 2, \dots, n_{a}

(8)

where

m_{s}

and

m_{a}

are supporting and moving masses of the AMD,

u_{ak}

is the relative displacement of the kth AMD, and

n_{a}

is the number of AMDs. It is assumed that the same AMDs are used for control. Using equation (1), equation (8) can be expressed as follows

T_{ak} = \frac{1}{2} {\dot{q}}^{T} (m_{s} + m_{a}) W^{T} (x_{ak}, y_{ak}) W (x_{ak}, y_{ak}) \dot{q} + {\dot{q}}^{T} m_{a} W^{T} (x_{ak}, y_{ak}) {\dot{u}}_{ak} + \frac{1}{2} m_{a} {\dot{u}}_{ak}^{2}, k = 1, 2, \dots, n_{a}

(9)

Hence, the kinetic energy of all AMDs can be expressed as

T_{a} = \frac{1}{2} (m_{s} + m_{a}) {\dot{q}}^{T} {\bar{M}}_{a} \dot{q} + m_{a} {\dot{q}}^{T} B_{a}^{T} \dot{u} + \frac{1}{2} m_{a} {\dot{u}}^{T} \dot{u}

(10)

where

u = {[u_{a 1} u_{a 2} \dots u_{a n_{a}}]}^{T}

and

{\bar{M}}_{a} = \sum_{k = 1}^{n_{a}} W^{T} (x_{ak}, y_{ak}) W (x_{ak}, y_{ak})

(11a)

B_{a} = [\begin{matrix} W (x_{a 1}, y_{a 1}) \\ W (x_{a 2}, y_{a 2}) \\ ⋮ \\ W (x_{a n_{a}}, y_{a n_{a}}) \end{matrix}]

(11b)

Using kinetic and potential energies given by equations (3) and (10), equations of motion for this system including disturbance can be written as

M \ddot{q} + Kq = - m_{a} B_{a}^{T} {\ddot{u}}_{a} + d

(12)

where

M = ρ_{p} h ab {\bar{M}}_{p} + (m_{s} + m_{a}) {\bar{M}}_{a}

(13a)

K = \frac{Db}{a^{3}} {\bar{K}}_{p}

(13b)

are mass and stiffness matrices, and

d

is the disturbance vector.

Free vibration problem results in the following eigenvalue problem.

(K - ω^{2} M) A = 0

(14)

By solving the eigenvalue problem, we can obtain the eigenvector matrix $U$ , which satisfies the orthonormality condition.

U^{T} M U = I

(15a)

U^{T} K U = Λ

(15b)

where

I

is an

n \times n

identity matrix

Λ = [\begin{matrix} ω_{1}^{2} & 0 & \dots & 0 \\ 0 & ω_{2}^{2} & \dots & 0 \\ ⋮ & ⋮ & ⋱ & ⋮ \\ 0 & 0 & ⋮ & ω_{n}^{2} \end{matrix}]

(16)

in which

ω_{i} (i = 1, 2, \dots, n)

is the ith natural frequency. If we apply the modal transformation,

q = Ur

, into equation (12) and use the orthonormality condition given by equation (15), then we can obtain modal equations of motion.

\ddot{r} + 2 Z Ω \dot{r} + Λ r = - m_{a} {\bar{B}}_{a}^{T} {\ddot{u}}_{a} + \bar{d}

(17)

where

{\bar{B}}_{a} = B_{a} U, \bar{d} = U^{T} d

and

Ω = [\begin{matrix} ω_{1} & 0 & \dots & 0 \\ 0 & ω_{2} & \dots & 0 \\ ⋮ & ⋮ & ⋱ & ⋮ \\ 0 & 0 & ⋮ & ω_{n} \end{matrix}]

(18a)

Z = [\begin{matrix} ζ_{1} & 0 & \dots & 0 \\ 0 & ζ_{2} & \dots & 0 \\ ⋮ & ⋮ & ⋱ & ⋮ \\ 0 & 0 & ⋮ & ζ_{n} \end{matrix}]

(18b)

In deriving equation (17), modal damping was assumed.

Displacement at the mounting point can be written as

w_{ak} = w (x_{ak}, y_{ak}, t) = W (x_{ak}, y_{al}) q (t), k = 1, 2, \dots, n_{a}

(19)

Hence, vector equation for displacement measurement can be expressed as

w_{a} = B_{a} q = B_{a} U r = {\bar{B}}_{a} r

(20)

where

w_{a} = {[w_{a 1} w_{a 2} \dots w_{a n_{a}}]}^{T}

. In reality, displacements cannot be measured, although accelerations can be measured by accelerometers. The AMD proposed in this study uses acceleration at the mounting point. Thus, measurement equation in the real world can be written as

{\ddot{w}}_{a} = {\bar{B}}_{a} \ddot{r}

(21)

MIMO modal-space NAF control

There are always more natural modes than available sensors. It is not desirable to control all modes with only a few AMDs. Let us consider a case in which the same number of AMDs as the number of natural modes to be controlled is used. Hence, we may write ${\bar{B}}_{a} = [{\bar{B}}_{ac} {\bar{B}}_{au}]$ . This results in a reduced-order system for control design. If we divide the generalized displacement vector into generalized displacements to be controlled by the AMD and the remaining displacements, such that $r = {[r_{c}^{T} r_{u}^{T}]}^{T}$ , then controlled and uncontrolled equations of motion and the sensor equation can be written as

{\ddot{r}}_{c} + 2 Z_{c} Ω_{c} {\dot{r}}_{c} + Λ_{c} r_{c} = - m_{a} {\bar{B}}_{ac}^{T} {\ddot{u}}_{a} + {\bar{d}}_{c}

(22)

{\ddot{r}}_{u} + 2 Z_{u} Ω_{u} {\dot{r}}_{u} + Λ_{u} r_{u} = - m_{a} {\bar{B}}_{au}^{T} {\ddot{u}}_{a} + {\bar{d}}_{u}

(23)

{\ddot{w}}_{a} = {\bar{B}}_{ac} {\ddot{r}}_{c} + {\bar{B}}_{au} {\ddot{r}}_{u}

(24)

The MIMO modal-space NAF control for equation (22) can then be designed as

\ddot{Q} + 2 Z_{f} Ω_{c} \dot{Q} + Λ_{c} Q = - G_{m} {\ddot{r}}_{c}

(25)

where

Z_{f}

is the damping matrix for NAF control. The MIMO modal-space NAF control can be obtained by the following relation

m_{a} {\bar{B}}_{ac}^{T} u_{a} = G_{m} Q

(26)

where

Q

is an

n_{a} \times 1

modal control force vector, and

G_{m}

is an

n_{a} \times n_{a}

diagonal matrix whose diagonal is the gain for each mode. Hence, the control displacement of AMD can be obtained by the following

u_{a} = \frac{1}{m_{a}} {\bar{B}}_{ac}^{- T} G_{m} Q

(27)

Let us neglect residual modes for the control design, then equation (28) can be derived from equation (21).

{\ddot{r}}_{c} = {\bar{B}}_{ac}^{- 1} {\ddot{w}}_{a}

(28)

Of course, ${\bar{B}}_{ac}$ should be invertible, meaning that AMDs should not be placed at the node of natural modes of interest. Inserting equation (26) into equation (22), we can obtain the following

{\ddot{r}}_{c} + 2 Z_{c} Ω_{c} {\dot{r}}_{c} + Λ_{c} r_{c} = - G_{m} \ddot{Q} + {\bar{d}}_{c}

(29)

Combining equations (25) and (29), we can obtain the following

[\begin{matrix} I_{n_{a} \times n_{a}} & G_{m} \\ G_{m} & I_{n_{a} \times n_{a}} \end{matrix}] {\begin{matrix} {\ddot{r}}_{c} \\ \ddot{Q} \end{matrix}} + [\begin{matrix} 2 Z_{c} Ω_{c} & 0_{n_{a} \times n_{a}} \\ 0_{n_{a} \times n_{a}} & 2 Z_{f} Ω_{c} \end{matrix}] {\begin{matrix} {\dot{r}}_{c} \\ \dot{Q} \end{matrix}} + [\begin{matrix} Λ_{c} & 0_{n_{a} \times n_{a}} \\ 0_{n_{a} \times n_{a}} & Λ_{c} \end{matrix}] {\begin{matrix} r_{c} \\ Q \end{matrix}} = {\begin{matrix} {\bar{d}}_{c} \\ 0_{n_{a} \times 1} \end{matrix}}

(30)

Hence, the stability condition can be obtained as

stable if G_{m}^{2} < I

(31)

The stability condition, equation (31), is static, meaning that the stability does not depend on frequency. Equation (31) also implies that the stability is guaranteed if the gain matrix is small enough.

By using equations (25), (27), and (28), the MIMO modal-space NAF controller using the acceleration measurement can be written as

U_{a} (s) = - \frac{1}{m_{a}} {\bar{B}}_{ac}^{- T} G_{m}^{2} {[s^{2} I + 2 s Z_{f} Ω_{c} + Λ_{c}]}^{- 1} {\bar{B}}_{ac}^{- 1} {\ddot{W}}_{a} (s)

(32)

However, control and observer spill over problems can occur because we take only parts of ${\bar{B}}_{a}$ . Let us take a look at uncontrolled modes. The control force is spilled into the remaining modes and the resulting equation can be written as

{\ddot{r}}_{u} + 2 Z_{u} Ω_{u} {\dot{r}}_{u} + Λ_{u} r_{u} = - {\bar{B}}_{au}^{T} {\bar{B}}_{ac}^{- T} G_{m} Q + {\bar{d}}_{u}

(33)

Hence, higher uncontrolled modes are disturbed because of spill over problems. It should be noted here that positive position feedback (PPF)²⁷ control does not cause this kind of spill over problem.

In addition, modal displacements to be controlled can be calculated based on measured relative displacements of the plate at positions of AMDs. Hence, we may write

{\hat{r}}_{c} = {\bar{B}}_{ac}^{- 1} w_{a}

(34)

where

{\hat{r}}_{c}

represent the measured modal coordinate. Inserting equation (24) into equation (34), it can be seen that computed modal displacements of natural modes to be controlled consist of actual modal displacements to be controlled and uncontrolled actual modal displacements.

{\hat{r}}_{c} = r_{c} + {\bar{B}}_{ac}^{- 1} B_{au} r_{u} = Rr

(35)

where

R = [I {\bar{B}}_{ac}^{- 1} {\bar{B}}_{au}]

. It can be readily seen from equation (35) that uncontrolled modal displacements spill into controlled modal displacements to be used for MIMO modal-space NAF control.

\ddot{r} + 2 Z Ω \dot{r} + Λ r = - m_{a}^{2} {\bar{B}}_{a}^{T} {\bar{B}}_{ac}^{- T} G_{m} Q + \bar{d}

(36)

\ddot{Q} + 2 Z_{f} Ω_{c} \dot{Q} + Λ_{c} Q = - G_{m} {\hat{\ddot{r}}}_{c} = - G_{m} R \ddot{r}

(37)

Equations (36) and (37) represent fully coupled equations of motion with MIMO modal-space NAF control that include both controlled and the uncontrolled natural modes. Combined equations can be written in matrix form as

[\begin{matrix} I_{n \times n} & m_{a}^{2} {\bar{B}}_{a}^{T} {\bar{B}}_{ac}^{- T} G_{m} \\ G_{m} R & I_{n_{a} \times n_{a}} \end{matrix}] {\begin{matrix} \ddot{r} \\ \ddot{Q} \end{matrix}} + [\begin{matrix} 2 Z Ω & 0_{n \times n_{a}} \\ 0_{n_{a} \times n} & 2 Z_{f} Ω_{c} \end{matrix}] {\begin{matrix} \dot{r} \\ \dot{Q} \end{matrix}} + [\begin{matrix} Λ & 0_{n \times n_{a}} \\ 0_{n_{a} \times n} & Λ_{c} \end{matrix}] {\begin{matrix} r \\ Q \end{matrix}} = {\begin{matrix} \bar{d} \\ 0_{n_{a} \times 1} \end{matrix}}

(38)

Hence, the stability condition becomes

stable if m_{a}^{2} {\bar{B}}_{a}^{T} {\bar{B}}_{ac}^{- T} G_{m}^{2} R < I

(39)

The stability condition given by equation (39) when considering control and observer spill overs are still static. However, it tells us that the gain should be small enough to satisfy this condition.

Decentralized MIMO NAF control

The control law given by equation (32) represents a fully coupled MIMO modal-space NAF controller. This implies that modal accelerations are computed from accelerations at each point, and the modal control is then converted into the real command using matrix computations. For simplicity, one can design each AMD to cope with each natural mode in a decentralized way. Hence, we may write

{\ddot{u}}_{a} + 2 Z_{f} Ω_{c} {\dot{u}}_{a} + Λ_{c} u_{a} = - G_{s} {\ddot{w}}_{a} = - G_{s} {\bar{B}}_{a} \ddot{r}

(40)

Equation (40) implies that the desired displacement of the AMD is calculated by the acceleration measured at that point, and each AMD is assigned to tackle a certain mode. Combining equations (17) and (40), we can obtain

[\begin{matrix} I_{n \times n} & m_{a} {\bar{B}}_{a}^{T} \\ G_{s} {\bar{B}}_{a} & I_{n_{a} \times n_{a}} \end{matrix}] {\begin{matrix} \ddot{r} \\ {\ddot{u}}_{a} \end{matrix}} + [\begin{matrix} 2 Z Ω & 0_{n \times n_{a}} \\ 0_{n_{a} \times n} & 2 Z_{f} Ω_{c} \end{matrix}] {\begin{matrix} \dot{r} \\ {\dot{u}}_{a} \end{matrix}} + [\begin{matrix} Λ & 0_{n \times n_{a}} \\ 0_{n_{a} \times n} & Λ_{c} \end{matrix}] {\begin{matrix} r \\ u_{a} \end{matrix}} = {\begin{matrix} \bar{d} \\ 0_{n_{a} \times 1} \end{matrix}}

(41)

Hence, the stability condition can be obtained as

stable if m_{a} {\bar{B}}_{a}^{T} G_{s} {\bar{B}}_{a} < I

(42)

The stability condition, equation (42), is also static. Equation (42) implies that the stability is guaranteed if the gain matrix is small enough. In addition, control or observer spill over problems would not occur in this case.

Using equation (40), the decentralized MIMO NAF controller using the acceleration measurement can be written as

U_{a} (s) = - G_{s} {[s^{2} I + 2 s Z_{f} Ω_{c} + Λ_{c}]}^{- 1} {\ddot{W}}_{a} (s)

(43)

Unlike the fully coupled MIMO modal-space NAF controller given by equation (32), the decentralized MIMO NAF controller does not require matrix computation. Thus, it is simple to use.

Numerical results

As a numerical example, a rectangular aluminum plate having $a = 4 m$ , $b = 3 m$ , $h = 3 mm$ , $ρ = 2770 {kg/m}^{3}$ , and $E = 70 GPa$ was considered. Total mass of the plate was 99.72 kg. Three AMDs having support mass of 0.5 kg and moving mass of 1 kg were mounted at $(2 m, 1.5 m)$ , $(1 m, 1.5 m)$ , $(2 m, 0.75 m)$ , respectively. These positions were selected because they are the anti-nodal points of three natural modes, resulting in maximum efficiency for each mode. Six admissible functions were considered in each direction. Thus, the total degree-of-freedom was 36, which was high enough for computing natural frequencies. Natural frequencies of the plate without AMDs were 1.2446, 2.5887, 3.6341, 4.8289, 4.9782, 7.2185, 7.6167, 7.9652, 8.9608, and 10.3548 Hz. Natural frequencies of the plate with AMDs were 1.1752, 2.5145, 3.5245, 4.5713, 4.9782, 6.9414, 7.2806, 7.9652, 8.7104, and 10.3548 Hz, respectively. There were slight decreases in natural frequencies because of added masses of AMDs. The three lowest natural frequencies were used as filter frequencies of NAF controllers. Figures 2 and 3 show natural mode shapes of the plate with and without AMDs. As shown in Figures 2 and 3, AMDs affected only higher modes, particularly the sixth and seventh modes. The modal damping factor was assumed to be 0.6%, and filter damping factor of 0.3 was chosen. The MIMO modal-space NAF controller was designed with $G_{m} = 0.2 I_{3 \times 3}$ and the decentralized MIMO NAF controller was designed with $G_{s} = 1.5 I_{3 \times 3}$ to suppress the three lowest natural modes. In order to obtain transfer function, the excitation point, $(2.4242, 0.9091) m$ , was chosen. Figure 4 shows frequency response function at the excitation point when the MIMO modal-space NAF controller is used. As shown in Figure 4, the introduction of AMD increased the damping of the three lowest modes. Figure 5 shows impulse response at the excitation point. As can be seen from Figure 5, vibrations were quickly suppressed by the fully coupled MIMO modal-space NAF controller. Figure 6 shows frequency response function at the excitation point when the decentralized MIMO NAF controller is used. Figure 7 shows impulse response. It can be said that the decentralized MIMO NAF controller is equally effective in suppressing vibrations.

Figure 2.

Natural modes of a plate without AMDs: (a) first mode, (b) second mode, (c) third mode, (d) fourth mode, (e) fifth mode, (f) sixth mode, (g) seventh mode and (h) eighth mode.

Figure 3.

Natural modes of a plate with AMDs: (a) first mode, (b) second mode, (c) third mode, (d) fourth mode, (e) fifth mode, (f) sixth mode, (g) seventh mode and (h) eighth mode.

Figure 4.

Frequency response function by MIMO modal-space NAF controller: (a) magnitude and (b) phase.

Figure 5.

Impulse response by MIMO modal-space NAF controller.

Figure 6.

Frequency response function by decentralized MIMO NAF controller: (a) magnitude and (b) phase.

Figure 7.

Impulse response by decentralized MIMO NAF controller.

Conclusions

In this study, vibrations of a rectangular plate are suppressed by using AMDs and NAF control. It was assumed that the AMD could track the desired displacement. The NAF controller uses accelerometer signal directly. Thus, it does not require an additional algorithm to compute velocity or displacement from the acceleration signal.

Two types of NAF controllers were proposed. One was the MIMO modal-space NAF controller and another was the decentralized MIMO NAF controller. The MIMO modal-space NAF controller was designed in modal space. The control command was then converted into a real command. It is a fully coupled controller. It can be used to selectively control damping of a target natural mode. On the other hand, the decentralized MIMO NAF controller was designed without matrix computation. However, the AMD operated by a decentralized MIMO NAF controller was responsible for suppressing only the assigned natural mode.

The MIMO modal-space NAF controller inherently had control and observer spill-over problems. Because modal coordinates were estimated from sensor measurements, control that suppressed a certain mode might have undesirable effects on higher modes. In this study, control and observer spill over problems that could occur when using MIMO modal-space NAF control were investigated theoretically. The stability was also checked. It was found that control and observer spill-over problems could be overcome as long as the gain matrix was kept small.

The decentralized MIMO NAF controller has a practical advantage over the MIMO modal-space NAF controller because it does not require matrix computation. In this study, the stability of the decentralized MIMO NAF controller was checked and found to be effective as long as the gain matrix was small, as in the case of the MIMO modal-space NAF control. Numerical results showed that both NAF controllers were effective in suppressing vibrations.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the research program of Dongguk University-Seoul, Korea.

References

McNamara

RJ.

Tuned mass dampers for buildings. J Struct Div 1977; 103: 1785–1798.

Marian

Giaralis

Optimal design of a novel tuned mass-damper-inerter (TMDI) passive vibration control configuration for stochastically support-excited structural systems. Prob Eng Mech 2014; 38: 156–164.

Pietrosanti

Angelis

Basili

Optimal design and performance evaluation of systems with tuned mass damper inerter (TMDI). Earthquake Eng Struct Dyn 2017; 46: 1367–1388.

Cao

Tuned tandem mass dampers-inerters with broadband high effectiveness for structures under white noise base excitations. Struct Control Health Monit 2019; 26: e2319.

Kobori

Koshika

Yamada

, et al. Seismic-response-controlled structure with active mass driver system Part 1: design. Earthquake Eng Struct Dyn 1991; 20: 133–149.

Kobori

Koshika

Yamada

, et al. Seismic-response-controlled structure with active mass driver system Part 2: verification. Earthquake Eng Struct Dyn 1991; 20: 151–166.

An optimum design methodology of active tuned mass damper for asymmetric structures.

Mech Syst Signal Process 2010; 24: 746–765.

Liu

Wang

Active multiple tune mass dampers: a new control strategy.

J Struct Eng. 2003; 129: 972.

Liu

Active multiple tune mass dampers for structures under the ground acceleration.

Earthquake Eng Struct Dyn 2002; 31: 1041–1052.

10.

Chung

Lin

Soong

, et al. Experimental study of active control for MDOF seismic structures. J Eng Mech 1989; 115: 1609–1627.

11.

Chang

Yang

Control of buildings using active tuned mass damper. J Eng Mech 1995; 121: 355–366.

12.

Ankireddi

Yang

TT.

Simple ATMD control methodology for tall building subject to wind loads. J Struct Eng 1996; 122: 83–91.

13.

Cao

Reinhorn

Soong

TT.

Design of active mass damper for a tall TV tower in Nanjing China. Eng Struct 1998; 21: 134–143.

14.

Cao

and Design of active tuned mass damper based on robust control. IEEE International Conference on Computer Science 2012; 2: 760–764.

15.

Samali

Al-Dawod

Performance of a five-storey benchmark model using active tuned mass damper and a fuzzy controller. Eng Struct 2003; 25: 1597–1610.

16.

Ricciardelli

Pizzimenti

Mattei

Passive and active mass damper control of the response of tall buildings to wind gustiness. Eng Struct 2003; 25: 1199–1209.

17.

Cao

Lin

QS.

New control strategies for active tuned mass damper systems. Comput Struct 2004; 82: 2341–2350.

18.

Wang

Lin

YH.

Vibration control of a tall building subjected to earthquake excitation. J Sound Vibr 2007; 299: 757–773.

19.

Yang

Shin

Lee

, et al. Active vibration control of structure by Active Mass Damper and Multi-Modal Negative Acceleration Feedback control algorithm. J Sound Vibr 2017; 392: 18–30.

20.

Talib

Shin

Kwak

MK.

Designing multi-input multi-output modal-space negative acceleration feedback control for vibration suppression of structures using active mass dampers. J Sound Vibr 2019; 439: 77–98.

21.

Lazarus

Crawley

Lin

CY.

Multivariable high-authority control of plate-like active structures. J Guid Control Dyn 1996; 19: 1357–1363.

22.

Cheng

Modeling and vibration control of a plate coupled with piezoelectric material. Compos Struct 2003; 62: 155–162.

23.

Kwak

Yang

DH.

Dynamic modelling and active vibration control of a submerged rectangular plate equipped with piezoelectric sensors and actuators. J Fluids Struct 2015; 54: 848–867.

24.

Kwak

Han

SB.

Free vibration analysis of rectangular plate with a hole by means of independent coordinate coupling method. J Sound Vibr 2007; 306: 12–30.

25.

Blevins

RD.

Formulas for natural frequency and mode shape. FL: Robert E. Krieger Publishing Co., Inc., 1987.

26.

Kwak

Yang

DH.

Free vibration analysis of cantilever plate partially submerged into a fluid. J Fluids Struct 2013; 40: 25–41.

27.

Fanson

Caughy

TK.

Positive position feedback control for large space structures. AIAA J 1987; 28: 717–724.