1. Introduction
It is well known that excessive vibrations represent a destructive dynamic condition. Repetitive operation or external forces cause simultaneous movement that can resonate through the machine, building or bridge to a dangerous magnitude. One common method to address the vibration control issue on flexible mechanical structures is through linear or nonlinear passive devices, taking advantage of the physical properties of the system itself. In order to prevent undesirable consequences of vibrations, this method modifies, mostly, mass, damping and stiffness properties with regard to initial configuration of the principal structure. Passive control techniques are characterized by implementation of devices in structures that do not require any external energy source to reduce mechanical vibrations [
1,
2].
Within the approach of linear passive vibration control, one of the widely used devices is the tuned mass damper (TMD), which consists of a mass, a spring and a viscous damper. This device is frequently implemented because of its properties such as effectiveness, reliability and low costs, with applications such as machinery and civil structures [
3].
TMD was initially used at the beginning of the past century since its conceptualization was applied for the first time by Frahm to reduce movement of ships as well as vibrations of the ships hull [
4]. After a major development on its dynamic behavior, TMD was designed to control the structural dynamic response on different topics. Ormondroyd and Den Hartog [
5] came upon that a TMD, with a damping element, can suppress the amplitude of the primary system in a wider frequency range, followed by a detailed discussion of the optimization that adjusts the damping parameters. Application of a passive vibration control scheme in flexible structures, using a TMD, is not just for controlling the dynamic response on lateral loads but also to mitigate torsional displacements in buildings with significant torsional coupling [
6].
Simplicity of tuned mass dampers makes them the most used devices for vibration control in buildings with great height. Guo and Chen [
7] proposed an innovative technique for using multiple TMDs to control partial loads on the ground in a limited number of floors. They indicated using numerical results that the use of multiple TMDs can effectively alter the distribution of natural frequencies as well as reduce the frequency/transient responses of the structure. Nowadays, research related to the study and implementation of a TMD remains a current topic, for example, for vibration control of adjacent twin buildings or using it in combination with the tapering method in order to control the dynamic response of super-tall buildings [
8,
9].
However, the use of nonlinear devices for passive vibration control is a relevant issue due to dynamic behaviors that may occur and do not happen in linear vibration systems [
10]. Usually, a nonlinear vibration absorber is implemented in order to overcome possible drawbacks due to the use of a TMD [
11]. There is a classification of nonlinear vibration absorbers called autoparametric absorbers. This type of nonlinear systems differ from the traditional TMD, mainly because these have nonlinear coupling between at least two vibration modes, satisfying the so-called autoparametric condition (external and internal resonance condition), which are certainly related with parametric excitation. Autoparametric absorbers are specifically used where a primary system is being excited close to one of its principal parametric resonances—that is, the worst case situation in a physical structure. When the autoparametric interaction occurs between two subsystems, there is a great energy transfer to the autoparametric absorber.
Autoparametric absorbers have been designed to mitigate resonant oscillations due to the advantages that this type of systems present in their frequency response function in comparison with the classic vibration absorber (TMD). Ibrahim and Heo [
12] and Dahlberg [
13] described how a continuous cantilever beam absorber with tip mass, oriented in the same direction with the motion of the primary system, can be implemented with better attenuation properties than those obtained with classical TMD. Cuvalci et al. [
14] defined an absorption region for an autoparametric vibration absorber for a single degree of freedom primary system under sinusoidal and random excitations. They experimentally determined the parameters that influence the effectiveness of a nonlinear vibration absorber. Hui and Ng [
15] presented the implementation of autoparametric phenomena to reduce symmetrical vibration of a curved beam/panel under external harmonic excitation showing that internal energy transfer of a first symmetric mode into first anti-symmetric mode in a curved panel is one example of autoparametric vibration absorber effect. Abundis-Fong et al. [
16] developed an optimum design of an autoparametric absorber (cantilever beam configuration) coupled to a resonant oscillator where the implementation of the nonlinear absorber was obtained by using an approximation of the nonlinear frequency response function, computed via a perturbation method. Recently, Ting Tan et al. [
17] used the nonlinear saturation principle and 1:2 internal resonance in the design of the piezoelectric autoparametric vibration absorber for vibration suppression and energy harvesting. Moreover, active vibration absorbers can be implemented to suppress undesirable vibrations and simultaneous tracking of reference trajectories by implementing on-line algebraic parameter identification methods [
18]. In this paper, the theoretical framework for algebraic parametrical identification of linear dynamical systems introduced in [
19] is extended to the on-line modal parameter estimation problem of a class of harmonically perturbed flexible structures. It has been theoretically proved that algebraic identification is robust against noise and polynomial disturbances. On-line algebraic identification has been also applied for synthesis of model-free control strategies [
20] and numerical differentiation techniques of noisy measurement signals [
21].
In this article, we are interested in implementing an on-line modal parameters identification technique to two different passive vibration control schemes for a flexible structure, justifying its tuning by means of a modal decomposition, in order to make an experimental comparison from a dynamic, frequency and energy approach. The work is structured as follows. In
Section 2, the tuning conditions of the passive control schemes implemented are required. The dynamic representation and modal decomposition of a flexible structure with
n degrees of freedom are presented in
Section 3. In
Section 4 an on-line algebraic identification scheme for estimating the amplitude and frequency of a harmonic excitation is detailed. Experiments using linear (TMD) and nonlinear (autoparametric system) vibration absorber for a flexible structure discretized in a finite number of degrees of freedom are presented in
Section 5. Finally, conclusions are given in
Section 6.
3. A Flexible Structure with n Degrees of Freedom
Real flexible structures are continuous elastic systems which have an infinite number of degrees of freedom. Therefore, their dynamical analysis commonly entails an approximation. This consists of describing their dynamical behaviour through the use of a finite number of degrees of freedom, as many as necessary to ensure enough accuracy. Flexible structures are usually described as discretized multiple degree-of-freedom systems.
In this work, the primary or main system is composed of a flexible structure which simulates a building-like structure with
n degrees of freedom as depicted in
Figure 2. Dynamical system behavior can be described by [
27]
where
represents the relative displacement vector with respect to a main frame reference.
, and
K are
matrices of mass, damping and stiffness of the primary system.
is the input force vector such that
. It is possible to transform the effect of the ground acceleration,
, as a force acting on each mass by using the vector
. In order to express the forces acting on the system in a compact way, let us define:
Substitution of (
9) into (
8) leads to the linear model of a vibrating mechanical system with
n degrees of freedom with base excitation
A linear or nonlinear dynamic vibration absorber can be then connected over the nth floor of the primary system for passive vibration control purposes. The main objective is to decrease harmonic oscillations disturbing a specific vibration mode, near some natural frequency of the main system.
Primary system (
10) can be then expressed in terms of modal or principal coordinates
as follows [
28,
29]
with
where
stands for a vector of modal or principal coordinates.
represents an external force vector. Parameters
and
are the equivalent modal damping and the natural frequency.
is the
modal matrix given by
The mathematical model of the flexible structure (
10) can be hence described as
Using notation of operational calculus of Mikusiński [
28,
30], modal model (
14) can be also written as
with
,
. Modal coordinates
can be then expressed as
where constants
and
depend on unknown initial conditions of the system at
. From Equations (
11) and (
12), displacements
can be also expressed in notation of operational calculus as
This representation leads to the expression
with
where
is an unique polynomial independent of the output variable
available for parametrical estimation. Furthermore,
is known as the characteristic polynomial of the dynamical system. Here,
are constants that depend on the roots of
. Roots of the characteristic polynomial (
20) provide the natural frequencies and damping ratios of the flexible structure.
Equation (
11) describes the dynamic characteristics of the
ith vibration mode and it constitutes what is known as the modal model, i.e., this describes the system through its modal properties (mode shapes, natural frequencies, and damping ratios), as opposed to the spatial model (
10), where the system is described by its spatial properties (
M,
C, and
K) [
3,
29].
Once the original system (
10) has been decoupled by (
11) and (
12), it is possible to focus on the vibration mode to be attenuated or damped. Therefore, passive vibration control on the flexible mechanical structure converges to the primary system configuration of a single degree of freedom with a secondary system (vibration absorber). The vibration absorption device can be tuned to the original system either using (
3), for the linear case, or (
6) and (
7) for the nonlinear case (autoparametric absorber). It is important to mention that although the vibration absorber is tuned to attenuate a single vibration mode, in fact, it is capable of attenuating the dynamic response of all degrees of freedom of the primary flexible structure. This situation is shown in the experimental results section of the present study.
4. Time-Domain and On-Line Algebraic Identification of the Harmonic Excitation
Consider the mechanical system shown in
Figure 2 with a dynamic behavior described by the mathematical model (
10). It is possible to use measurements of a single position variable
in the synthesis of a modal model based on coefficients
of the characteristic polynomial (
20) [
30,
31]. When the ground acceleration produces the excitation force
as described in (
8), individual component forces
are harmonic. We can then express those forces in notation of operational calculus as follows
where the frequency
and the amplitude
A is the same for each component of the input force vector
. Hence, by substituting (
21) in (
18) we have:
where
is the polynomial defined by (
20). Multiplying (
22) by
we obtain:
with
. Constants
,
, depend on the modal matrix entries
, the mass
of each floor and the initial conditions. Then, we derive
times Equation (
24) with respect to the variable
s in order to annihilate the polynomial disturbance
. Next, the result is multiplied by
and transformed back to the time domain to yield
where
with
Here,
,
, and notation
is employed to represent iterated integrals of the form
where
m is a positive integer. By considering that
is a positive number, we have that
Then, we get an expression for on-line and time-domain estimation of the excitation frequency
. It is possible to avoid singularities and, at the same time, to get a considerably smoother estimation by integrating Equation (
25) two times respect to time as follows
with
where
is used like a low pass filter with cut frequency defined by the parameter
. In this way, the numerical estimation process implementation of the frequency of the excitation force is smoothed.
The proposed identification scheme is capable to identify the frequency of the excitation force, by using measurements of some of the outputs
with a similar performance. In addition, estimation of the parameter
does not depend, neither on the amplitude
A of the ground acceleration nor the unknown system initial conditions. Moreover, since measurement noise could be considered as fast external fluctuations corrupting output signals
, iterated integrals operate like low pass filters [
19,
32]. In this sense, on-line algebraic identification is, at some extent, robust against noise [
19,
20,
21]. However, measurement signals should be suitably conditioned and pre-filtered to reduce harmful noise levels. Otherwise, for operational scenarios where highly large noise is manifested, a considerable degradation on the estimation performance could occur.
5. Experimental Results
A series of experiments were performed considering two different case studies. A six-story building-like structure discretized in six degrees of freedom with a linear vibration absorber is first implemented. The other case is a different structure discretized in three degrees of freedom on which a nonlinear vibration absorber is implemented. Both of the structures are made of aluminum alloy and supported by columns of the same material.
The first primary system to be considered is the six-story building-like structure shown in
Figure 3. The sensing elements are IEPE accelerometers attached to each story of the structure and a high resolution rotary encoder used to take angle measurements of the linear pendulum absorber. An electromechanical shaker acts as a source of excitation in conjunction with the frictionless slider for performing experimental modal analysis and performance tests of the vibrations absorber. The impedance head is used to measure the force input, these measurements are needed for the application of experimental modal analysis techniques.
Modal parameters of the structure are obtained using experimental modal analysis, by averaging two of the most popular excitation techniques: impact hammer and sine sweep [
28,
29]. We then process the corresponding acceleration measurements by applying a FRF-based multiple degrees-of-freedom method (MDoF) [
29]. In particular, we apply the well founded RFP curve fitting method described in [
34], that is, we express the experimental FRF in terms of rational fraction polynomials so the coefficients
of the characteristic polynomial (
20) can be determined. The transfer function or estimated FRF for one measurable degree of freedom of the structure can be expressed as a function of the complex variable
s, so that:
where
. Then, the curve fitting method RFP leads to an estimated or fitted FRF (
37) such that:
where
is the difference or error between the measured FRF and the estimated FRF, which is constructed by using the coefficients
and
,
is defined by Equations (
17) and (
18). The description of the technical details of the curve fitting method are explained in [
35]. The experimental and fitted frequency response functions for the primary system is described in
Figure 4, where the peaks shown correspond only to transverse vibration modes. We use only measurements of the sixth floor of the structure to apply the curve fitting method. Notice the close relation between the experimental FRF (in solid blue line) and the estimated FRF (in dotted black line) especially near the peaks or resonances. The modal parameters, frequency and damping ratios, obtained by solving the roots of the estimated characteristic polynomial or the denominator of (
37) are reported in
Table 1.
The coefficients of the characteristic polynomial are reported in
Table 2. Those coefficients are used for synthesis of an on-line algebraic identifier (
30) for the harmonic excitation force.
Then, the characteristic polynomial of the system is
where the constants
are those reported in
Table 2.
5.1. Non-Linearity Analysis
In order to determine the presence of nonlinearities in the dynamic behavior of the six-story building-like structure, we apply an analysis of the frequency response shown in
Figure 4. Thus, the non-linearity index defined by (
35) is determined. The Nyquist diagram in
Figure 5 shows the comparison between the original FRF in solid blue line and the corresponding Hilbert transformation in dotted lines. The numerical value of the calculated no-linearity index is
Therefore, it is assumed that the dynamic behavior of the six-story building structure is dominantly linear and the experimental determination of the coefficients of the characteristic polynomial are valid and reliable for the synthesis of the algebraic identifier for the excitation frequency. Moreover, this linear behavior adds a foundation to the reliability of the estimated modal parameters.
5.2. Application of a Linear Absorber
The primary system mentioned above is coupled to a pendulum vibration absorber (configured to work as a TMD) on the top floor. In order to excite the complete system, an electromechanical shaker is used as shown in
Figure 3.
Equations of motions related to the system depicted in
Figure 3 which is discretized in seven degrees of freedom, considering small angular displacements in the secondary system and submitted to an harmonic forced excitation, are given by
where
is a vector containing the lateral displacements of each floor,
is the angular displacement of the pendulum absorber. The mass, damping and stiffness matrices are represented by
M,
C and
K, respectively. The input vector is
. According to (
9), the input force
is given by:
where
is the second derivative of ground motion with respect to time. The functions
and
are:
Parameters associated with the pendulum absorber are its mass (
), length (
L) as well as viscous damping (
). In the experiments, the base of the structure is directly affected by a ground motion
with amplitude
and excitation frequency
. Due to we are interested in attenuating the first transverse vibration mode, the value of
is close to the first resonance frequency
. According to (
3), the linear pendulum vibration absorber must be tuned in such a way that
Usually, under stable operating conditions, the linear pendulum vibration absorber can be appropriately designed to passively control any vibration mode associated to () of the main system, depending on the narrow frequency bandwidth where is acting.
The dynamic behavior of the six-story building-like structure when the TMD-pendulum absorber is not tuned is described in
Figure 6, where the amplitudes of vibration in stable state are
mm,
mm,
mm,
mm,
mm and
mm.
Figure 7 shows the dynamic response of the flexible structure when experimentally (
45) is satisfied. Now, the amplitudes of vibration in stable state are
mm,
mm,
mm,
mm,
mm and
mm, resulting in an average absorption percentage close to
.
The linear absorber time history response, when the secondary system is tuned, is described in
Figure 8. The stable state amplitude is
rad.
Figure 9 shows the performance of the on-line estimator (
31) for the excitation force frequency. Acceptable estimations are obtained after
s. The numerical value of the estimated frequency is
Hz.
5.3. Application of a Nonlinear Absorber
With the intention of showing that it is also possible to use the modal decomposition for the tuning of a nonlinear vibration absorber (independently of the number of degrees of freedom in the primary system), certain experiments were carried out with a pendulum type absorber, implemented in autoparametric form.
The experiment consisted in a flexible structure discretized in three degrees of freedom where a nonlinear vibration absorber (pendulum type) is coupled as shown in
Figure 10. In this second case, we also use a high resolution rotational encoder to take angle measurements of the autoparametric pendulum absorber and a IEPE accelerometer, attached to the third story of the structure to take vibrations measurements. Once again, the electromechanical shaker acts as a source of excitation in conjunction with the frictionless slider. In this configuration, it is necessary to add a spring element in the secondary system in order to provide the pendulum with potential energy since its rotational dynamics occurs in a horizontal plane, so there are not gravity effects to be taken into account (see
Figure 10). It results in the following nonlinear dynamic model
where
M,
C and
K are the
matrices of mass, damping and stiffness.
is an input vector. Functions
and
are defined as:
The parameter associated with the pendulum absorber is its equivalent mass , where L is the pendulum length.
It is evident that the nonlinear coupling between both subsystems is by means of
and
which are defined by (
48). These nonlinear functions allow the implementation of the autoparametric pendulum absorber in one specific vibration mode of the primary system.
The frequency response function of the primary system with nonlinear pendulum absorber is described in
Figure 11. In this type of nonlinear vibration systems, the frequency response function remains unchanged, i.e., no additional peaks are added in the aforementioned graph (which is a typical dynamic situation using a TMD absorber). Only the natural frequencies of the flexible structure change slightly in value due to the added mass (mass absorber).
5.4. Nonlinearity Analysis
Both primary systems analyzed in this work are assumed to be linear due to the nominal operating conditions considered, in addition to the characteristics of its construction material. Nevertheless, the vibration absorption scheme used to mitigate undesired effects of the first vibration mode of the structure is inherently non linear. In order to evaluate the influence of the nonlinear autoparametric pendulum absorber, we perform a nonlinearity analysis, using the nonlinearity index defined by (
36), on the complete system (included the nonlinear vibrations absorber). The Nyquist diagram of the corresponding FRF of the system is shown in
Figure 12 with its corresponding Hilbert transformation. The dashed black line, representing the Hilbert transformation of the FRF, does not show an important distortion.
The low distortion in the Nyquist diagram produced by the application of the Hilbert transformation suggests a linear behavior of the system in its nominal operational conditions. The nonlinearity index (
35) of this particular building like structure including the nonlinear vibrations absorber is:
Thus, we can use the coefficients of the characteristic polynomial reported in
Table 3 for the synthesis of an on-line algebraic identifier (
31) for the harmonic excitation force.
The modal parameters obtained from experimental modal analysis by applying the curve fitting method for the case when a pendulum absorber (nonlinear-type) is applied are given in
Table 4.
Coefficients of the characteristic polynomial are reported in
Table 3. Those coefficients are used for the synthesis of an on-line algebraic identifier (
30) for the harmonic excitation force.
Figure 13 shows the performance of the on-line estimator (
31) for the excitation frequency. Estimations are stable after only
s. The numerical value for the on-line estimation parameter is
Hz after
s.
On the other hand, the dynamic performance in energy terms of the system shown in
Figure 10, is described in
Figure 14 and
Figure 15, respectively. Once the transient response has disappeared the autoparametric absorber shows good dynamic performance dissipating a large percentage of the external energy supplied to the main structure.
6. Conclusions
In the present work, a comparison of a passive vibration control scheme applied to a flexible structure is studied experimentally. In order to mitigate resonant excitations, associated with the first mode of vibration, a pendulum absorber was implemented in two different configurations. It was observed that when the passive absorber was configured to work as an autoparametric system (nonlinear case), there were no significant changes in its frequency response function, this is because additional resonances are not introduced when this kind of dynamic vibration absorber is used, in contrast with the TMD configuration where this is a well known disadvantage.
Basically, the main advantage of an autoparametric vibration absorber is its property of high energy absorption exactly in resonant excitation. That is, when the external and internal resonance (autoresonance) conditions are tuned, the external energy affecting directly the main structure is transferred as kinetic energy (motion) to the autoparametric absorber and this situation results reasonable, because then the primary system can be protected from worst case dynamic conditions (resonance). It is important to mention that in both passive vibration absorbers, the effective damping should be small, otherwise the absorption capability would be significantly reduced. It is evident that there is a compromise about the amount of damping existing into the primary system and absorber, because stability and performance depend on this type of criteria. In general, small damping results in better vibration attenuation properties but at the same time makes the complete system more sensitive to endogenous perturbations. In addition, the performance of the on-line algebraic estimator for frequency of the excitation force showed a good efficiency using measurements of only one degree of freedom of the building-like structure, having the advantage of its estimation speed (it is achieved in less than a cycle of the original signal) in contrast with the fast-Fourier-transform (FFT)-based methods where at least one cycle of the signal is needed. Finally, the nonlinearity index tested here is easy to program and compute. In particular, we have assumed a value of to establish that a given system is dominantly linear. Subsequent investigation on passive/active dynamic vibration absorption based on on-line algebraic identification of modal parameters and excitation forces on nonlinear multi-delay flexible structures will be considered in future work. Future research work will consider the implementation of multiple-frequency vibration absorption devices as well.