Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 51, 1, pp. 225-234, Warsaw 2013 CONTROL OF CHARACTERISTICS OF MECHATRONIC SYSTEMS USING PIEZOELECTRIC MATERIALS Andrzej Buchacz, Marek Płaczek, Andrzej Wróbel Silesian University of Technology, Faculty of Mechanical Engineering, Gliwice, Poland e-mail: andrzej.buchacz @polsl.pl; marek.placzek@polsl.pl; andrzej.wrobel@polsl.pl This paper presents two examples of systems with piezoelectric transducers used in order to control vibration. First of them is a combination of a single piezoelectric plate glued on a mechanical subsystem surface. The piezoelectric plate is applied to suppress vibration of the system. The secondmodel consists of a few piezoelectric plates vibrating in a thickness mode. Such kind of systems is called the piezostack and is used to enhance the transducers sensitivity or to increase the displacement. In both cases, characteristics of the considered systems were generated and presented on charts. Key words: dynamics, vibration control, piezoelectricity, approximatemethods 1. Control of vibrations with piezoelectric transducers The piezoelectric effect discovered in 1880 by Pierre and Jacques Curie has found a lot of ap- plications in modern technical devices. Nowadays, there are a lot of applications of both – the direct and reverse piezoelectric effects. The phenomenon of production of an electric potential when stress is applied to the piezoelectric material is called the direct piezoelectric effect.While the phenomenon of production of strain when an electric field is applied is the reverse piezoelec- tric effect (Moheimani and Fleming, 2006). The direct piezoelectric effect application in passive vibration damping is being considered in this section. 1.1. Idea of passive damping of vibration In the nineteen eighties, a new idea for piezoelectric transducers application – passive dam- ping of vibrations – was developed. It is the possibility of dissipating mechanical energy with piezoelectric transducers shunted with passive electrical circuits. A piezoelectric transducer is bonded to some mechanical subsystem (for example a beam or a shaft) and external electric circuit is adjoined to the transducer clamps. This idea was described and experimentally in- vestigated by Forward (1981), and the first analytical models for such a kind of systems were developed byHagood and von Flotow (1991, 1994). There aremany commercial applications of this idea, such as products of ACX (Active Control eXperts inc.) – skis, snowboards or baseball bats with piezoelectric dampers (Yoshikawa et al., 1998). Two cases of passive damping of vibrations with piezoelectric transducers were described. In the first case, only a resistor is used as a shunting circuit, and in the second case, it is a circuit composed of a resistor and an inductor. Since this research, many authors have worked to improve this idea. For example Fleming et al. (2002) describedmultimode piezoelectric shunt damping systems. This work presents analysis of a mechatronic system with piezoelectric shunt damping of vibration using a singlePZTpiezoelectric transducerwith a resistor used as the shunting circuit. The piezoelectric transducer is used as a passive vibration damper. The resistor is adjoined to the transducer clamps. The dynamic characteristic of this system was assigned on the basis of 226 A. Buchacz et al. the approximate Galerkin method. Dynamic flexibilities of the considered mechatronic system and mechanical subsystem without piezoelectric transducer were juxtaposed. The influence of geometrical and material properties of the considered system on its dynamic flexibility was determined. The correct description of a designed system by its mathematical model during the design phase is a fundamental condition for its proper operation. A very precise mathematical model and method of the analysis of the system is very important. It is indispensable to take into consideration all geometrical andmaterial parameters of system components (including the glue layer between these elements) because the omission of the influence of one of them results in inaccuracy in the analysis of the system (Pietrzakowski, 2001; Buchacz and Płaczek, 2010a). Therefore, the processes of modelling and testing of one-dimensional vibrating mechatronic systems are presented. It is a part of work that was done at Gliwice research centre, where tasks of analysis and synthesis of mechatronic systems with piezoelectric transducers were considered. Tasks of ana- lysis and synthesis were also considered and presented in previous publications, see Dzitkowski andDymarek (2005, 2011), Zolkiewski (2010, 2011). For example, passive and activemechanical andmechatronic systems were analyzed byBiałas (2008, 2009) and investigations supported by computer-aided methods by Buchacz andWróbel (2010) and Dzitkowski (2004). 1.2. Considered system and its mathematical model The considered system (presented in Fig. 1) consists of a cantilever beam which has a rec- tangular constant cross-section, length l and Young’s modulus E. The piezoelectric transducer of length lp and thickness hp is bonded to the top surface of the beam by a glue layer of thickness hk andKirchhoff’s modulus G. The glue layer has homogeneous properties in overall length. The external shunting circuit is adjoined to the transducer clamps. The system is loaded with a harmonic force F(t), described by the equation F(t)=F0cos(ωt) (1.1) Fig. 1. The considered systemwith the external shunting circuit In accordance with the idea of passive damping of vibrations with piezoelectric materials, the mechanical system affects the piezoelectric transducer which generates electric charge and produces additional stiffness of electromechanical nature, dependent on the capacitance of the transducer and the adjoined external circuit. This shunting circuit introduces electric dissipation resulting in electronic damping of vibration (Kurnik, 2004). The piezoelectric transducer with the external shunting circuit was treated as a linear RC circuit as it is shown in Fig. 2, where RZ is the external electric circuit resistance and Cp is the capacitance of the piezoelectric transducer (Buchacz and Płaczek, 2010b, 2011a,b). Control of characteristics of mechatronic systems... 227 Fig. 2. Electrical representation of the piezoelectric transducer with the external shunting circuit The equation of flexural vibration of the beam was assigned in accordance with Bernoulli- Euler’s model of the beam, and the following assumptions were made: • the material of which the system is made is subjected to Hooke’s law, • the system has a continuous, linear mass distribution, • the system vibration is harmonic, • planes of sections that are perpendicular to the axis of the beam remain flat during defor- mation of the beam, • displacements are small compared with the dimensions of the system. Displacements of the beam cross sections were described by the coordinate y(x,t). The dynamic equation of motion of the beamwas assigned on the basis of elementary beam and transducer section dynamic equilibrium in agreement with d’Alembert’s principle with the assumption on the uniaxial, homogeneous strain of the transducer and pure shear of the glue layer. Heaviside’s function H(x) was introduced to curb the working space of the transducer within the range from x1 to x2 (Buchacz andPłaczek, 2011a,b). The obtained equation of beam motion can be described as ∂2y(x,t) ∂t2 =−a4 ∂4y(x,t) ∂x4 + c2 ∂ ∂x [εb(x,t)H−εk(x,t)H+λ1(t)H]+ δ(x− l) ρbAb F(t) (1.2) where a= 4 √ EbJb ρbAb c2 = Glp 2ρbhk H =H(x−x1)−H(x−x2) λ1(t)= d31 UC(t) hp (1.3) The symbols εb and εk are strains of the beam and glue layer surfaces, d31 is the piezoelectric constant and UC(t) is the electric voltage on the capacitance Cp of the piezoelectric transducer. The distribution of the external force was determined using Dirac’s delta function δ(x− l). Taking intoaccount constitutive equationswhich includethe relationshipbetweenmechanical and electrical properties of transducers (Preumont, 2006; Moheimani and Fleming, 2006), the equation of the piezoelectric transducer with the external electric circuit can be described as RZCp ∂UC(t) ∂t +UC(t)= lpbd31 Cps E 11 S1(x,t)+ lpbε T 33 UC(t) Cphp (1−k231) (1.4) where k231 = d231 sE11ε T 33 (1.5) 228 A. Buchacz et al. is the electromechanical coupling constant of the transducer and S1(x,t) is the strain of the transducer resulting from deflection of the beam. The symbols Es11 and ε T 33 are respectively Young’s modulus of the transducers under zero/constant electric field and dielectric constant under zero/constant stress. The electric voltage generated by the transducer as a result of its strain was assumed as UC(t)= |Up| ωCp|Z| sin(ωt+ϕ) (1.6) where |Z| and ϕ are absolute value and argument of the serial circuit impedance. Up is the electric voltage generated as a result of the strain of the transducer. Equations (1.2) and (1.4) create the mathematical model of the system. Dynamic flexibilities of the consideredmechatronic systemand itsmechanical subsystemwe- re assigned on the basis of approximateGalerkin’smethod.The solution to the beamdifferential motion equationwas defined as a product of time anddisplacement eigenfunctions (Buchacz and Płaczek, 2010c) y(x,t)=A ∞ ∑ n=1 sin [ (2n−1) π 2l x ] cos(ωt) (1.7) This equationmeets only two boundary conditions – displacements of the clamped and free ends of the beam, but the approximate method was verified – dynamic flexibility of the mechanical subsystem (the beam without the piezoelectric transducer) was assigned using exact Fourier’s and approximate methods and the obtained results were juxtaposed. Inexactness of the appro- ximate method was reduced by the correction coefficients introduced for the first three values of periodicity of the system (Buchacz andPłaczek, 2011a,b). It was proved that after correction this method can be used to analyse such a kind of systems, and the obtained results are very precise.Derivatives of the assumed solutionwere substituted into themathematicalmodel of the system. After transformations, the dynamic flexibility was calculated, and taking into account parameters of the system it was presented on charts. 1.3. Parameters of the system and obtained results Parameters of the systemarepresented inTable 1.Thedynamicflexibility of themechatronic system, assigned on the free end of the beam (x = l), separately for the first three natural frequencies (n= 1,2,3) is presented in Fig. 3. It is juxtaposed with the dynamic flexibility of the mechanical subsystem. Table 1.Parameters of the mechatronic system l=0.24m Eb =210000MPa b=0.04m ρb =7850kgm −3 hb =0.002m d31 =−240 ·10 −12mV−1 x1 =0.01m ε T 33 =2.65 ·10 −9Fm−1 x2 =0.06m s E 11 =63000MPa hp =0.0005m RZ =6000Ω hk =0.0001m G=1000MPa Using the developedmathematicalmodel of the considered system and the proposedmethod of analysis of the system, it is possible to analyse the influence of the system parameters on its characteristic. Effects of some selected parameters on the dynamic flexibility are presented in Figs. 4 and 5. Control of characteristics of mechatronic systems... 229 Fig. 3. Dynamic flexibility, (a) for the first natural frequency, (b) for the second natural frequency, (c) for the third natural frequency Fig. 4. Influence of modulus of shear elasticity of the glue layer, for n=1 Fig. 5. Influence of the piezoelectric constant d31 of the transducer material, for n=1 The presented mathematical model of the considered system is one of a series of models developed and presented in other publications. The developed models have different degrees of precision of the real system representation because different simplificationswere introduced.The aim of this study was to develop a series of mathematical models of the considered system, to verify them and to indicate an adequatemodel to accurately describe the phenomena occurring 230 A. Buchacz et al. in the system and effectively simplifymathematical calculations with reducing the required time for them (Buchacz and Płaczek, 2010a). 2. Vibration control of the piezoelectric stack In the construction and operation of mamymachatronic subassemblies, the piezoelectric pheno- menon has been utilised for a long time. For testing of machine elements, piezoelectric sensors are glued on their surfaces, whereas for health monitoring, transducers made of piezoelectric films are used. Piezoelectrics are often used in machine building in form of assemblies, subas- semblies or actuators. In somecases, theparameters of simplepiezoelectric plates are insufficient; sometimes we want to have a greater sensitivity for transducers, or expect more displacement for a piezoelectric actuator. In such cases, single piezoelectric plates are replaced by piezoelec- tric stacks. The piezoelectric stack is a composition of few thin piezoelectric layers, which are connected in parallel. This kind of monolithic ceramic construction is characterized by: • good conversion of energy, • low voltage control, • high force, • rapid response of the system. The view of piezoelectric stacks with electrical connectors was shown in Fig. 6. The stacks may have different shapes and dimensions – it depends on the kind of characteristics we want to get. Piezoelectric stacks are used if bigger displacement or bigger force on the piezoelectric surface is required. Fig. 6. A stack of piezoelectric plates Analyzing the reverse piezoelectric phenomenon, in other words, the generation of charge on the piezoelectric linings, it is worth noting that in the construction of piezoelectric transducers the longitudinal phenomenon ismost commonlyused.Plates used for buildingof suchaconverter are usually in the form of plates of rolls and connected in parallel. Through such connection, an increase in the obtained charge is obtained, and thereby an increase in the sensitivity of the transducer. Proper implementation of the piezoelectric systemwhich acts as a sensor or actuator depends on selection of geometric dimensions of the plate, or a group of plates, their numbers and their basic physical parameters. Figures 7, 8 and 9 show the influence of various parameters such as thickness of the plate, density of thematerial of which the plate is made on the characteristics of the system (Buchacz andWróbel, 2007; Wróbel, 2012). Control of characteristics of mechatronic systems... 231 Fig. 7. Characteristic of a piezoelectric plate in the frequency domain depending on thickness of the plate Fig. 8. Characteristic of a piezoelectric plate in the frequency domain depending on piezoelectric density Fig. 9. 3D characteristic of the piezoelectric plate in the frequency domain in function of material density of the plate Upto now, the proper selection has been realised by checkingwhether the resultsmeasure up the expectations through applying another, electrically connected, piezoelectric layer in order to increase the transducer sensitivity until achieving a satisfactory result. This section presents the possibility of controlling the characteristics of the piezoelectric stack consisting of PZT plates. The control in this case is done by changing the supply frequency, geometrical and material parameters. 232 A. Buchacz et al. Table 2.Parameters used in numerical simulations Number Symbol Value Unit 1 A 3 [cm2] 2 d 3 [mm] 3 cE33 5.8 [GPa] 4 ρ 7500 [g/cm3] All characteristics inFigs. 7 to9 showndisplacements in the thicknessmode. Infiniteamplitu- des at the resonances are due to disregarding damping in thematerial. In the current researches by the authors, this problem is under consideration and will be published in the future work. Analysing a systemconsisting of only twopiezoelectric plates and assuming that the two sub- systems are characterized by identical physicochemical and geometric parameters, and assuming the same thickness of individual plates d (i) 1 = 0.5d and d (i+1) 2 = 0.5d, the same characteristics as in the case of a single piezoelectric plate with thickness d1+d2 were obtained. 3. Conclusions The discussed subjet is important due to the increasing number of applications of both the simple and reverse piezoelectric phenomena in variousmodern technical devices. The process of modelling of technical deviceswithpiezoelectricmaterials is complex and requires large amounts of time because of complexity of the phenomena occurring in such systems. The correct descrip- tion of a given system by its mathematical model during the design phase is the fundamental condition for proper operation of the designed system. Amathematical model that provides the most accurate analysis of the system together with maximum simplification of used mathema- tical tools and minimum required amount of time should be incorporated. The identification of the optimal mathematical model that meets the assumed criteria is one of the purposes of research works realized at Gliwice research centre. It is an introduction to the task of synthesis of one-dimensional vibrating complex mechatronic systems. Tasks of analysis and synthesis of mechanical and mechatronic systems were considered at Gliwice centre in previous works (Bu- chacz andWróbel, 2008, 2010). Passive and activemechanical systems andmechatronic systems with piezoelectric transducers were analyzed byBuchacz andPłaczek (2011). The examinations were also supported by computer-aided methods (Buchacz andWróbel, 2010). Acknowledgements This work was partially supported by Polish Ministry of Science and Higher Education – National Science Centre as a part of the research project No. N501064440 (2011-2013). References 1. Białas K., 2009, Reverse task of passive and active mechanical system in torsional vibrations, Journal of Achievements in Materials and Manufacturing Engineering, 35, 2, 129-137 2. BiałasK., 2008,Polargraphsandstructuralnumbers in synthesis of activeandpassivemechanical systems, Journal of Achievements in Materials and Manufacturing Engineering, 30, 1, 43-50 3. Buchacz A., PłaczekM., 2010a,Development ofmathematicalmodel of amechatronic system, Solid State Phenomena, 164, 319-322 4. Buchacz A., Płaczek M., 2010b, Selection of parameters of external electric circuit for control of dynamic flexibility of a mechatronic system, Solid State Phenomena, 164, 323-326 Control of characteristics of mechatronic systems... 233 5. BuchaczA.,PłaczekM., 2010c,The approximateGalerkin’smethod in the vibratingmechatro- nic system’s investigation,Proceedings of The 14th International ConferenceModern Technologies, Quality and Innovation ModTech 2010, 147-150 6. Buchacz A., Płaczek M., 2011a,Modeling and investigation of one-dimensional flexural vibra- ting mechatronic systems with piezoelectric transducers, Advances in Piezoelectric Transducers, Farzad Ebrahimi (Edit.), InTech, Rijeka 7. Buchacz A., Płaczek M., 2011b,Modelling and Testing of One-Dimensional Vibrating Mecha- tronic Systems, Silesian University of Technology Publishing House, Gliwice [in Polish] 8. Buchacz A., Wróbel A., 2007, The investigation of the piezoelectric effect of bimorph system and his practical application in level sensors,Proceedings of the XII The Scientific Conference, of Vibroacoustics and Vibrotechnique, Jachranka, 63-66 9. Buchacz A., Wróbel A., 2008, Homogenization of mechatronical systems in modeling piezo- electric layer,Machine Dynamics Problems, 32, 3, 7-11 10. Buchacz A., Wróbel A., 2010, Computer-aided analysis of piezoelectric plates, Solid State Phenomena, 164, 339-342 11. Dymarek A., Dzitkowski T., 2005,Modelling and synthesis of discrete-continuous subsystems of machines with damping, Journal of Materials Processing Technology, 164, 1317-1326 12. Dzitkowski T., 2004, Computer-aided synthesis of discrete-continuous subsystems of machines with the assumed frequency spectrum representedbygraphs,Journal ofMaterials ProcessingTech- nology, 157, 144-149 13. Dzitkowski T., Dymarek A., 2011, Active synthesis of multiaxial drive systems using a com- parative method, Journal of Achievements in Materials and Manufacturing Engineering, 49, 2, 275-284 14. FlemingA.J., Behrens S., RezaMoheimani S.O., 2002,Optimization and implementation of multimodepiezoelectric shuntdamping systems,Proc. IEEE/ASMETransactions onMechatronics, 7, 87-94 15. Forward R.L., 1981, Electronic damping of orthogonal bending modes in a cylindrical mast – experiment, Journal of Spacecraft and Rockets, 18, 11-17 16. Hagood N.W., von Flotow A., 1991, Damping of structural vibrations with piezoelectric ma- terials and passive electric networks, Journal of Sound and Vibration, 146, 2, 243-268 17. HagoodN.W., vonFlotowA., 1994,Design of Passive Piezoelectric Damping for Space Struc- tures,NationalAeronauticsandSpaceAdministrationLangleyResearchCenter,Hampton,Virginia 18. Kurnik W., 2004, Damping of mechanical vibrations utilising shunted piezoelements, Machine Dynamics Problems, 28, 15-26 19. Moheimani S.O.R., Fleming A.J., 2006, Piezoelectric Transducers for Vibration Control and Damping, Springer, Londyn 20. Pietrzakowski M., 2001, Active damping of beams by piezoelectric system: effects of bonding layer properties, International Journal of Solids and Structures, 38, 7885-7897 21. Preumont A., 2006, Mechatronics: Dynamics of Electromechanical and Piezoelectric Systems, Springer 22. Yoshikawa S., Bogue A., Degon B., 1998, Commercial application of passive and active pie- zoelectric vibration control,Proc. of the Eleventh IEEE International Symposium on Applications of Ferroelectrics, 293-294 23. Wróbel A., 2012,Model of piezoelectric including material dumping, Proceedings of 16th Inter- national Conference ModTech, ISSN 2069-6736, 1061-1064 24. Zolkiewski S., 2010, Attenuation-frequency characteristics of beam systems in spatial motion, Solid State Phenomena, 164, 349-354 25. Zolkiewski S., 2011, Dynamic flexibility of the supported-clamped beam in transportation, Journal of Vibroengineering, 13, 810-816 234 A. Buchacz et al. Kontrola charakterystyk układów mechatronicznych przy użyciu materiałów piezoelektrycznych Streszczenie W artykule przedstawiono dwa przykłady układów mechatronicznych zawierających przetworniki piezoelektryczne stosowane w celu kontroli i stabilizacji drgań.W pierwszym układzie pojedyncza płyt- ka piezoelektryczna naklejana jest na powierzchni podukładu mechanicznego. Do przetwornika piezo- elektrycznego dołączony jest zewnętrzny, pasywny obwód elektryczny, przez co pełni on rolę szeroko- pasmowego, pasywnego tłumika drgań. Drugi z rozpatrywanych układów złożony jest ze stosu płytek piezoelektrycznych drgających grubościowo. Tego typu układy stosowane są w celu zwiększenia czułości przetworników stosowanych jako sensory lub zwiększenia odkształceń układuwprzypadku stosowania go w roliwzbudnika.Wprzypadkuobuomawianychukładówdrgających zaproponowanometodę ich analizy oraz wyznaczono charakterystyki, których przebiegi przedstawiono w formie graficznej. Manuscript received June 6, 2011; accepted for print June 23, 2012