#### Abstract

This paper is about the nonlinear analysis of a piezoelectric controlled Ziegler column. The piezoelectric controller, here referred to as Tuned Piezoelectric Damper (TPD), possesses evanescent characteristics and, moreover, it is tuned to the first natural frequency of the mechanical system, thus resembling the well-known Tuned Mass Damper. This means that the flow of energy between mechanical and electrical subsystems is driven by the resonance (Den Hartog principle) and magnified by the singularity of the evanescent electrical characteristics. Numerical simulations, showing how the proposed control strategy is effective in increasing the linear stability domain and decreasing the amplitude of the limit-cycles in the postcritical range, are presented.

#### 1. Introduction

Structures subjected to follower forces exhibit rich and interesting dynamic features. Among these latter, the most studied are the loss of stability and the onset of several paradoxes (e.g., related to the damping). Moreover, it is interesting to remark that control strategies that are optimal for reducing vibrations can have a detrimental effect on stability (see, e.g., [1]). Indeed, developing a unified control strategy optimized for damping vibrations in the linear and postcritical regimes, and at the same time increasing the stability zone, is a challenging task.

Since their first introduction in the early 1990s, piezoelectric based controllers have attracted the interest of researchers and their industrial applications have increased over the past years. These controllers are usually divided into two main categories: active and passive controllers. The distinction is based on whether or not the controller is able to pump energy into the mechanical system. As it is well known, active controllers can be really effective but also suffer from stability issues. Examples of active controllers can be found in [2, 3]. This is why our choice is to focus our attention on passive controllers of the type firstly introduced in [3]. As they are usually implemented by tuning the resonance(s) of the electric circuit to one (or more) natural frequency(ies) of the mechanical system, they can represent a valid alternative to the classic Tuned Mass Dampers (TMDs) [4–8]. Their advantage is mostly due to the fact that they are more versatile, simpler to tune, and in most configurations lighter than classic TMD. Moreover, piezoelectric controllers can be used also in microdevices [9].

In our previous studies [10, 11], we focused our attention on several kinds of piezoelectric control strategies, applied to discrete linear structures. In these papers, mechanical and electrical subsystems were studied in their linear and nonlinear dynamics, respectively. In [10], three different controllers were proposed and it was proved that, at various extents, they all have a beneficial effect on stability. In [11], attention was focused on a two-DOF system, namely, the nonlinear Ziegler column, which was taken as case study, and the amplitude of the limit-cycle was investigated with and without the piezoelectric controller: in particular, in that paper we used one of the three controllers defined in [10], namely, the SNRC (Singular Nonresonant Controller), which remembers the Nonlinear Energy Sink (see, e.g., [12, 13]) even if, differently from the NES, the coupling between the mechanical system and the controller is linear and of gyroscopic type.

The objective of this paper is to study the behavior of the nonlinear Ziegler column equipped with another one of the three controllers introduced in [10], which is called the SRC (Singular Resonant Controller), which resembles the working mechanism of the TMD, since it is resonant and possesses evanescent electric properties. This type of controller, therefore, can be also referred to as Tuned Piezoelectric Damper (TPD). It is important to remark that we include into the category of TPD only passive circuits, that is, circuits that cannot pump energy into the mechanical systems. Few examples of controllers having these characteristics can be found in [14, 15] for collocated control and in [16–24] for distributed control; moreover, if the circuit parameters are well optimized, good damping results can be obtained even without resonant elements (see, e.g., [25]).

The paper is structured as follows. In Section 2, the discrete model of a piezoelectric controlled nonlinear Ziegler column is recalled. Section 3 is devoted to the presentation of the numerical results. Finally, in Section 4 some conclusions and a summary of the main results are exposed.

#### 2. The Model

The system under study is a piezoelectric controlled upward double-pendulum, shown in Figure 1. The equations of motion for this system are derived in [11] for finite rotations and will be here briefly recalled. The system consists of a two hinged weightless rigid bars of equal length , carrying two concentrated masses, at the common hinge (point in the figure) and at the tip (point in the figure). The bars are viscoelastically constrained at the hinges by (a) linear springs of stiffness and and (b) linear dashpots having viscosity coefficients and , respectively. The system is loaded at the free end by a follower force of intensity , which keeps its direction parallel to the upper bar. Moreover, the pendulum is equipped with a piezoelectric device of stiffness (here considered negligible with respect to the stiffness of the springs), capacitance , and coupling coefficient , which is placed at the ground hinge, and it is connected to a one-node LR resonant circuit, of inductance and resistance (sketched in Figure 1), and to the ground. When considering infinitesimal rotations of the two bars, the system reduces (except for the masses values adopted here) to that analyzed in [10]; moreover, when the controller is removed, the system degenerates into the well-known Ziegler column [26].

The state of the system is described by taking as Lagrangian coordinates the rotations of the two bars, and , and the flux linkage associated with the piezoelectric actuator. We recall here that the flux linkage is defined as the time integral of the voltage ; that is, .

The equations of motion for the coupled system are obtained by following these steps:(1)The internal dissipation (due to damping) and nonconservative action (due to the follower force) are neglected, so that we can describe the system in terms of its Lagrangian, which is the sum of three components: (a) the mechanical Lagrangian , concerning the structure without any control devices, (b) the electrical Lagrangian , relevant to the electrical circuit, and (c) the piezoelectric Lagrangian , referred to as the piezoelectric device.(2)The internal dissipation and the nonconservative action are accounted for through the Extended Hamilton Principle. For more details about this process, the reader can refer to [10, 11].

The obtained nondimensional system of equations has the following form:where the superposed dot denotes a time derivative and where we introduced the following quantities (accounting for , , and ) for nondimensionalization:in which and are scaling capacitance and flux linkage, respectively. Here, is the load parameter (taken as bifurcation parameter); is the electromechanical coupling parameter; , , and (all positive and reals) are electrical parameters, here referred to as the “electrical mass,” “electrical damping,” and “electrical stiffness,” respectively.

#### 3. Numerical Results

Numerical results, relevant to the uncontrolled and controlled Ziegler column, in both linear and nonlinear regimes are discussed in this section. It should be remarked that all the bifurcation diagrams discussed in the following have been obtained via a numerical continuation procedure, which has been directly applied to system (1), in order to compute the amplitudes of the limit-cycles occurring in the postcritical regime.

##### 3.1. Uncontrolled Column

This section is devoted to briefly resume the critical and postcritical behavior of the uncontrolled Ziegler column with the aim of preparing the discussion on the effects of the controller, given in the next section.

The linear stability domain of the uncontrolled column, that is, when , can be exactly determined by making use of the Routh-Hurwitz criterion on the (complete degree-4) characteristic equation of the algebraic eigenvalue problem associated with the linearized equations (1) (see, e.g., [27]). The criterion leads to the finding of a critical locus in the -space, known in the literature as the “Whitney’s umbrella” surface [28, 29], whose equation reads as follows:where is the critical load of the undamped system at which a circulatory Hopf bifurcation occurs.

The linear bifurcation diagram of the uncontrolled Ziegler column is displayed in Figure 2. In particular, Figure 2(a) shows the critical surface, given by (3), in the 3D domain ; it separates the stable states (the region marked with in the figure) from the unstable ones (the region marked with in the figure). The points on the surface represent the Hopf bifurcation states, that is, except for those on the -axis that, indeed, are regarded as marginally stable (undamped) systems for which a circulatory Hopf bifurcation occurs at (see, e.g., [30]). In Figure 2(b), the contour lines are displayed: it is apparent that the effect of damping is mostly detrimental in the whole -plane, since, except for a small region (filled in gray in the figure) close to an “optimal direction” (dashed line in the figure), the critical load of the damped system is lower with respect to . Moreover, it can be shown that on the left side of the optimal direction the stability is governed by the first mode, while on the right side it is governed by the second one.

**(a)**

**(b)**

The postcritical behavior of the uncontrolled Ziegler column is now discussed. In particular, we will refer to two systems, marked with a black dot and a label I or II in Figure 2(b), respectively, far and close enough to the critical load of the circulatory system, namely,(i)case study I: and , entailing , for which damping has a* strong destabilizing* effect (−62%);(ii)case study II: and , entailing , for which damping has a* moderate destabilizing* effect (−13%). In Figures 3 and 4, the bifurcation diagrams of the uncontrolled Ziegler column are represented for the case studies I and II, respectively. Here, the amplitudes of the motion components, and , are plotted* versus* the bifurcation parameter . It is apparent that also when the bifurcation parameter slightly exceeds its critical value , large amplitudes limit-cycles manifest themselves. This is the consequence of the destabilizing effect of damping yet discussed with regard to the linear behavior of the column, which persists also in the postcritical regime. As a matter of fact, when the increment of the load with respect to the critical value, , is equal, for example, to , we found rad, rad in case I and rad, rad. Therefore, the higher the destabilizing effect of damping on linear stability, the higher the amplitude of the limit-cycle (for the same increment of load with respect the critical one).

**(a)**

**(b)**

**(a)**

**(b)**

##### 3.2. Controlled Column

In this section, we analyze the behavior of the Ziegler column equipped with the SRC (Singular Resonant Controller) issued from the paper [10]. As it was proved there through a perturbation technique, the energy flow between mechanical and electrical subsystems is driven, in this controller, by the resonance (Den Hartog principle [31]) and magnified by the singularity of the evanescent electrical characteristics. Thus, according to this control strategy we need to have , , , and , being a small parameter; moreover, in order to render the electrical oscillator resonant to the mechanical one, the tuning , where is one of the two column’s frequencies at , must be performed. It is important to remark that such a controller resembles the well-known Tuned Mass Damper [4–8] even if the coupling is here linear and is of gyroscopic type; thus, the controller discussed in the present paper can be regarded as a Tuned Piezoelectric Damper (TPD).

According to this strategy, the following electrical parameters have been fixed: and ; moreover, has been chosen in order to tune the electrical frequency to the first mechanical one since, in both case studies, the first mode is the candidate to become unstable, namely, () for case study I and () for case study II. The coupling parameter has been considered variable, with the aim of estimating the sensitivity of the controller’s performances with respect to an increase of the electromechanical coupling and, in particular, it has been fixed to or .

Concerning the linear stability of the controlled system, the effects of the controller on the critical load can be evaluated in Figures 5 and 6 where the bifurcation diagrams of the controlled (black curves) and uncontrolled (gray curves) Ziegler column, for cases I and II, respectively, are represented. It is apparent that, in both case studies, the controller is able to* shift forward* the uncontrolled bifurcation diagrams, thus entailing an increase of the stability region; in particular, we found the following in the two cases: (I) (+54%), when , and (+98%), when ; (II) (+5%), when , and (+9%), when . Therefore, the controller is much more efficient in case study I where, remarkably, the damping has a strong destabilizing effect.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

Figures 5 and 6 illustrate also the influence of control on the postcritical behavior of the column: indeed, the TPD is able to* slightly flatten* the curves on the -axis, at least in a range of load values, thus reducing the limit-cycles amplitude of the uncontrolled system. As a matter of fact, we found, when , for the case study I rad (−18%), rad (−16%), when , and rad (−33%), rad (−21%), when .

It is important to remark that, as a result of the shift, the controller reduces the amplitude of the limit-cycle occurring at the same level of load. This effect is displayed in Figure 7 where, for case study I, when , and for the two considered values of , the time histories, relevant to uncontrolled (light gray curves in the figure) and controlled (dark gray curves in the figure) systems, show the benefits of the control. Finally, the flattening of the curves can be seen in Figure 8 where numerical integrations relevant to the uncontrolled system (first column in the figure) and the controlled system (second column in the figure), obtained for two different values of , at the same distance from the respective critical loads, are displayed, for case study I, when and . In Figures 8(a), 8(b), and 8(c) the projections of the trajectories onto phase planes are plotted, showing the contraction of the amplitude of the limit-cycle.

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

#### 4. Conclusions

In this paper, the effects of a TPD controller on the critical and postcritical behavior of the Ziegler column have been addressed.

The nonlinear equations of motion, for the coupled Piezo-Electro-Mechanical (PEM) system, have been recalled. The adopted control strategy first introduced in [10] is the SRC, namely, the Singular Resonant Controller, and its effects on the nonlinear dynamics of the system have been deeply investigated.

The main results can be summarized in the following points:(i)The analysis of the postcritical scenario shows that the controller is able to shift forward the bifurcation diagrams of the uncontrolled system.(ii)The controller induces a contraction of the limit-cycle amplitudes.(iii)Both previous effects are much more evident when increasing the value of the electromechanical coupling. The obtained results are satisfying and show that the proposed control strategy has a beneficial effect on the stability of the system. Nevertheless, to evaluate the real performances of this controller, a further optimization of the parameters is required. This will be the object of forthcoming papers.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work was granted by the Italian Ministry of University and Research (MIUR), under the PRIN10-11 program, Project no. 2010MBJK5B.