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 l, carrying two concentrated masses, m1≔2m at the common hinge (point B in the figure) and m2≔m at the tip (point C in the figure). The bars are viscoelastically constrained at the hinges by (a) linear springs of stiffness k1≔k and k2≔k and (b) linear dashpots having viscosity coefficients c1 and c2, respectively. The system is loaded at the free end by a follower force of intensity F, which keeps its direction parallel to the upper bar. Moreover, the pendulum is equipped with a piezoelectric device of stiffness kp (here considered negligible with respect to the stiffness of the springs), capacitance Cp, and coupling coefficient g, which is placed at the ground hinge, and it is connected to a one-node LR resonant circuit, of inductance L and resistance R (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].
Nonlinear controlled Ziegler column.
The state of the system is described by taking as Lagrangian coordinates the rotations of the two bars, ϑ1t and ϑ2t, and the flux linkage ψ(t) associated with the piezoelectric actuator. We recall here that the flux linkage is defined as the time integral of the voltage V; that is, ψ=∫Vdt.
The equations of motion for the coupled system are obtained by following these steps:
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 Lm, concerning the structure without any control devices, (b) the electrical Lagrangian Le, relevant to the electrical circuit, and (c) the piezoelectric Lagrangian Lp, referred to as the piezoelectric device.
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:(1)3ϑ¨1+ϑ¨2cosϑ1-ϑ2+ϑ˙22sinϑ1-ϑ2+ξ1+ξ2ϑ˙1-ξ2ϑ˙2+2ϑ1-ϑ2-μsinϑ1-ϑ2-γψ˙=0,ϑ¨2+ϑ¨1cosϑ1-ϑ2-ϑ˙12sinϑ1-ϑ2+ξ2ϑ˙2-ϑ˙1+ϑ2-ϑ1=0,νeψ¨+ξeψ˙+κeψ+γϑ˙1=0,where the superposed dot denotes a time derivative and where we introduced the following quantities (accounting for m1=2m, m2=m, and k1=k2=k) for nondimensionalization:(2)τ=ωt,ω2=kml2,μ=Fmlω2,ξ1=c1ml2ω,ξ2=c2ml2ω,ψ~=ψψ0,νe=CpC0,ξe=1RC0ω,κe=1LC0ω2,γ=gωlmC0,ψ0=lmC0,in which C0 and ψ0 are scaling capacitance and flux linkage, respectively. Here, 0<μ∈R is the load parameter (taken as bifurcation parameter); 0<γ∈R is the electromechanical coupling parameter; νe, ξe, and κe (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 γ=0, 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 μ,ξ1,ξ2-space, known in the literature as the “Whitney’s umbrella” surface [28, 29], whose equation reads as follows:(3)μ=μc+ξ1ξ22-3-222ξ1-4+52ξ22ξ1+ξ2ξ1+6ξ2,where μc:=7/2-2≃2.09 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 μ,ξ1,ξ2; it separates the stable states (the region marked with S in the figure) from the unstable ones (the region marked with U 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 μ=μc (see, e.g., [30]). In Figure 2(b), the contour lines μ=const are displayed: it is apparent that the effect of damping is mostly detrimental in the whole ξ1,ξ2-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 μd is lower with respect to μc. 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.
Uncontrolled Ziegler column: (a) critical manifold in the μ,ξ1,ξ2-parameter space: S is stable region and U is unstable region; (b) μ-isolines.
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,
case study I: ξ1=0.05 and ξ2=0.2, entailing μd≃0.80, for which damping has a strong destabilizing effect (−62%);
case study II: ξ1=0.3 and ξ2=0.15, entailing μd≃1.81, 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, maxϑ1 and maxϑ2, are plotted versus the bifurcation parameter μ. It is apparent that also when the bifurcation parameter slightly exceeds its critical value μd, 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, δμ≔μ-μd, is equal, for example, to δμ=0.3, we found maxϑ1≃1.03 rad, maxϑ2≃1.62 rad in case I and maxϑ1≃0.52 rad, maxϑ2≃1.13 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).
Bifurcation diagrams of the uncontrolled Ziegler column for the case study I: (a) maxϑ1 versus μ; (b) maxϑ2 versus μ.
Bifurcation diagrams of the uncontrolled Ziegler column for the case study II: (a) maxϑ1 versus μ; (b) maxϑ2 versus μ.
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 κe=O(ε), νe=O(ε), ξe=Oε3/2, and γ=Oε, ε being a small parameter; moreover, in order to render the electrical oscillator resonant to the mechanical one, the tuning κe/νe=ωd2, where ωd is one of the two column’s frequencies at μ=μd, 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: νe=0.1 and ξe=0.03; moreover, κe 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, κe=0.02 (ωd≃0.447) for case study I and κe=0.0375 (ωd≃0.612) 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 γ=0.05 or γ=0.1.
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) μd≃1.23 (+54%), when γ=0.05, and μd≃1.58 (+98%), when γ=0.1; (II) μd≃1.90 (+5%), when γ=0.05, and μd≃1.98 (+9%), when γ=0.1. Therefore, the controller is much more efficient in case study I where, remarkably, the damping has a strong destabilizing effect.
Bifurcation diagrams of the controlled Ziegler column for the case study I: (a) maxϑ1 versus μ; (b) maxϑ2 versus μ; (c) maxψ versus μ.
Bifurcation diagrams of the controlled Ziegler column for the case study II: (a) maxϑ1 versus μ; (b) maxϑ2 versus μ; (c) maxψ versus μ.
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 δμ=0.3, for the case study I maxϑ1≃0.84 rad (−18%), maxϑ2≃1.36 rad (−16%), when γ=0.05, and maxϑ1≃0.69 rad (−33%), maxϑ2≃1.28 rad (−21%), when γ=0.1.
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 μ=1.7, 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 δμ=0.3 and γ=0.1. 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.
Time histories of ϑ1 and ϑ2 for the case study I, when μ=1.7, for the uncontrolled (light gray curves) and the controlled (dark gray curves) systems when (a) γ=0.05 and (b) γ=0.1.
Numerical phase portraits for the case study I, when γ=0.1 and δμ=0.3, for the uncontrolled (first column, μ≃1.105) and the controlled (second column, μ≃1.88) systems.
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:
The analysis of the postcritical scenario shows that the controller is able to shift forward the bifurcation diagrams of the uncontrolled system.
The controller induces a contraction of the limit-cycle amplitudes.
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.
D’AnnibaleF.RosiG.LuongoA.On the failure of the ‘Similar Piezoelectric Control’ in preventing loss of stability by nonconservative positional forces20156641949196810.1007/s00033-014-0477-72-s2.0-84910047601ElliottS. J.GardonioP.SorsT. C.BrennanM. J.Active vibroacoustic control with multiple local feedback loops2002111290891510.1121/1.14338102-s2.0-0036158569HagoodN. W.von FlotowA.Damping of structural vibrations with piezoelectric materials and passive electrical networks1991146224326810.1016/0022-460X(91)90762-92-s2.0-0026142048BlevinsR. D.1990New York, NY, USAVan Nostrand ReinholdAbeM.FujinoY.Dynamic characterization of multiple tuned mass dampers and some design formulas199423881383510.1002/eqe.42902308022-s2.0-0028668353RanaR.SoongT. T.Parametric study and simplified design of tuned mass dampers199820319320410.1016/S0141-0296(97)00078-32-s2.0-0032019631GattulliV.Di FabioF.LuongoA.Simple and double Hopf bifurcations in aeroelastic oscillators with tuned mass dampers20013382-318720110.1016/s0016-0032(00)00077-62-s2.0-0035276277GattulliV.Di FabioF.LuongoA.One to one resonant double Hopf bifurcation in aeroelastic oscillators with tuned mass dampers2003262220121710.1016/S0022-460X(02)01135-52-s2.0-0037667775AndreausU.PlacidiL.RegaG.Microcantilever dynamics in tapping mode atomic force microscopy via higher eigenmodes analysis2013113221422430210.1063/1.48084462-s2.0-84879197166D'AnnibaleF.RosiG.LuongoA.Linear stability of piezoelectric-controlled discrete mechanical systems under nonconservative positional forces201550382583910.1007/s11012-014-0037-42-s2.0-84906521645D'AnnibaleF.RosiG.LuongoA.Piezoelectric control of Hopf bifurcations: a nonlinear discrete case studyIn pressVakakisA. F.BergmanL. A.GendelmanO. V.GladwellG. M.KerschenG.LeeY. S.McFarlandD. M.2009156Dordrecht, The NetherlandsSpringerSolid Mechanics and Its ApplicationsLuongoA.ZulliD.Dynamic analysis of externally excited NES-controlled systems via a mixed Multiple Scale/Harmonic Balance algorithm20127032049206110.1007/s11071-012-0597-62-s2.0-84870878488GiorgioI.CullaA.Del VescovoD.Multimode vibration control using several piezoelectric transducers shunted with a multiterminal network200979985987910.1007/s00419-008-0258-x2-s2.0-68449089760RosiG.PaccapeliR.OllivierF.PougetJ.Optimization of piezoelectric patch positioning for passive sound radiation control of plates201319565867310.1177/10775463124372362-s2.0-84875919026AlessandroniS.AndreausU.Dell'IsolaF.PorfiriM.Piezo-ElectroMechanical (PEM) Kirchhoff-Love plates200423468970210.1016/j.euromechsol.2004.03.0032-s2.0-2942606291MauriniC.dell'IsolaF.Del VescovoD.Comparison of piezoelectronic networks acting as distributed vibration absorbers20041851243127110.1016/s0888-3270(03)00082-72-s2.0-2142749438dell'IsolaF.MauriniC.PorfiriM.Passive damping of beam vibrations through distributed electric networks and piezoelectric transducers: prototype design and experimental validation200413229930810.1088/0964-1726/13/2/0082-s2.0-2142815266dell'IsolaF.VidoliS.Continuum modelling of piezoelectromechanical truss beams: an application to vibration damping199868111910.1007/s0041900501422-s2.0-0031996955Dell'IsolaF.HennekeE. G.PorfiriM.Synthesis of electrical networks interconnecting PZT actuators to damp mechanical vibrations2002141–44174242-s2.0-0035788814Dell'IsolaF.VestroniF.VidoliS.A class of electro-mechanical systems: linear and nonlinear dynamics20024014771GiorgioI.GalantucciL.Della CorteA.Del VescovoD.Piezo-electromechanical smart materials with distributed arrays of piezoelectric transducers: current and upcoming applications20154741051108410.3233/jae-140148PorfiriM.dell'IsolaF.SantiniE.Modeling and design of passive electric networks interconnecting piezoelectric transducers for distributed vibration control200521269872-s2.0-17744400088VidoliS.dell'IsolaF.Modal coupling in one-dimensional electromechanical structured continua20001411-2375010.1007/bf011768062-s2.0-0033685234RosiG.PougetJ.Dell'IsolaF.Control of sound radiation and transmission by a piezoelectric plate with an optimized resistive electrode201029585987010.1016/j.euromechsol.2010.02.0142-s2.0-77955471091ZieglerH.Die Stabilitätskriterien der Elastomechanik1952201495610.1007/bf005367962-s2.0-34347142065KirillovO. N.A theory of the destabilization paradox in non-conservative systems20051743-414516610.1007/s00707-004-0194-y2-s2.0-15544379723WhitneyH.The general type of singularity of a set of 2n-1 smooth functions of n variables1943101161172KirillovO. N.VerhulstF.Paradoxes of dissipation-induced destabilization or who opened Whitney's umbrella?201090646248810.1002/zamm.2009003152-s2.0-77952638613LuongoA.D'AnnibaleF.On the destabilizing effect of damping on discrete and continuous circulatory systems2014333246723674110.1016/j.jsv.2014.07.0302-s2.0-84907546216Den HartogJ.1956New York, NY, USAMcGraw-Hill