Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 41, 1, pp. 155-167, Warsaw 2003 ACTIVE AND PASSIVE DAMPING OF VIBRATIONS BY TORSIONAL DAMPER DURING STEADY STATE MOTION OF A POWER TRANSMISSION SYSTEM Zbigniew Skup Institute of Machine Design Fundamentals, Warsaw University of Technology e-mail: zskup@ipbm.simr.edu.pl This paper presents a theoretical study of the process of damping of nonlinear vibrations in a two-massmodel of amechanical systemwith a torsion damper. The steady-state motion of the system subject to har- monic excitation is considered on the assumption of a uniform frequency and constant amplitude of the forcing torque. Simultaneous structural frictionphenomena (passivedamping) and thepiezoelectric effect (active damping) are been considered as well. The problem is considered on the assumption of a uniform unit pressure distribution between the contac- ting surfaces of the friction discs and plunger. The aim of the analysis is to asses the influence of geometric parameters, external load, unit pres- sure and electric parameters on the resonance curves of the steady-state vibrations. The equations of motion of the examined system are solved bymeans of the Van der Pol method. Key words: active damping, vibrations, torsional damper, structural friction 1. Introduction Facing the complexity of problems concerned with structural friction in mechanical systems, including transmission systems with clutches, simplifi- cations of the friction model are being assumed. The elastic strain of discs and plunger as well as the influence of the piezoelectric effect on damping of vibrations in such systems are taken into account. The author based his con- siderations on previously derived physical formulations of a frictional torsion damper. In the design process of a power transmission system and selection of 156 Z.Skup existingmodels of dampers, basic standard calculationmethods are applied. It is essential, however, to take into account the natural source of vibration dam- ping bymeans of the structural friction and piezoelectric effect. Development in new materials have made it possible to create the so-called ”intelligent” materials (electrorheological fluids FL,ML; alloys with shapememory, piezo- electric polymers).Thesematerials change their properties under the influence of magnetic and electric field, temperature or mechanical stress. The basic electromechanical properties of the piezoelectric material – lead zirconate titanate (PZT) are presented in the paper in terms of their capacity to damp torsional vibrations in the considered mechanical system. Piezoce- ramics exhibit natural shear effect, which is three times stronger than the longitudinal one. This phenomenon can be successfully utilized in torsional systems for vibration control, as piezoceramics constitute perfect elements for actuator applications. Active control is achieved through closed loopwith pro- portional and velocity feedback (proportional-plus-derivative controller). Meng-Kao Yeh and Chih-Yuan Chin (1994) proved the applicability of piezoelectric sensors for measuring torsional vibration of shafts. Introduc- tory theoretical studies of shafts vibration induced by piezoelements ba- sed on PZT ceramics were presented by Kurnik and Przybyłowicz (1995), Przybyłowicz (1995), Tzou (1991). Experimental investigations dealing with an actively controlled torsional system were carried out by Chia-Chi Sung et al. (1994). 2. Equations of motion We assume a two-mass model of a mechanical system which contains a frictional torsiondamper,as shown inFig.1. Structural frictionoccursbetween the cooperating surfaces of discs 2 and plunger 1. Discs 2 are pressed down to plunger 1 bymeans of springs 3. The shaft is equipped with two piezoelectric elements: actuator 4 and sensor 5. The actuator is posed by a ring-shaped element of considerable thickness, thus, its moment of inertia must be taken into account. The sensor may be made of PZT or with PVDF (piezoelectric foil), yet it must be thin enough to be neglected in the total inertia balance of the system. The actuator and the sensor are electronically coupled with proportional andvelocity feedback ruling theperformanceof the thus arranged control system. Therefore, the equations of motion of the considered system can be written down as follows Active and passive damping of vibrations... 157 Fig. 1. Physical model of the considered friction torsional damper I1ϕ̈1+M(ϕ,A,ϕ̇)= Mm +M(t)+Ma (2.1) I2ϕ̈2−M(ϕ,A,ϕ̇)= 0 where ϕ1,ϕ2 – angulardisplacementsof theactive andpassivedamper ϕ,A – relative angular displacement of the discs and plunger and its amplitude, respectively I1,I2 – reducedmoments of inertia of themovable parts of the engine, actuator and plunger and of the discs in the damper M(t)+Mm – variable engine torque described by a constant ave- rage value Mm and a discrete torque M(t) in the form of harmonic excitation with a uniform frequ- ency and constant amplitude as in Giergiel (1990), Grudziński et al. (1992), Osiński (1986), Osiński and Kosior (1976), Skup (1987, 1991, 2002), Szadkowski andMorford (1992), Zagrodzki (1994) Ma – torque generated by the actuator. 158 Z.Skup M(ϕ,A,ϕ̇) – damper torque in a cycle represented by a structural hysteresis loop (Fig.1) dependent on the relative dis- placement, amplitude and its velocity signum. Therefore M(t)= M0cosωt (2.2) where M0 – excitation amplitude of the forcing torque ω – angular velocity of the excitation torque t – time. According to the studies presented in works by Skup (1976, 2002), the relationship between the damper torque and the relative angular displacement is as follows M(ϕ,A,ϕ̇)= k κ1 sgnϕ̇ [ 2 √ 1+κ1(ϕ−A)sgnϕ̇−1− √ 1−2κ1Asgnϕ̇ ] (2.3) Ma = κ2 sgnϕ̇ ( 1− √ 1−2κ1ϕs sgnϕ̇ ) where k = GI0 l I0 = πd4 32 κ1 = G2I20 2l2kzπ 2pµR5 kz = k1k2 k1+k2 k1 = Gh1 k2 = Gh2 κ2 = πGaGsd a 15d s 15lsd(D 3−d3)kpk 12GlaI0κ1ε0εs (2.4) The following denote: k – stiffness of the elastic shaft of the length l and diameter d, κ1 –nondimensional parameter, k1, k2,h1,h2 – discs andplunger stiffness and their thickness, µ – friction coefficient, p – pressure per unit area, R – external radius of the discs, G – shearmodulus, I0 – cross-sectional moment of inertia of the shaft, Ga – shear modulus of PZT material (for actuator), la –width of the actuator, d a 15 – coupling constant of the actuator, kp – gain factor introduced by the electronic circuit, d, D – inner and outer actuator diameter, ε0, εs – absolute and relative dielectric permittivity of the sensor, Gs – shear modulus of the sensor material, ls – width of the sensor, ds15 – electromechanical coupling constant of the sensor. Active and passive damping of vibrations... 159 3. Solution to the equation of motion After introducing the relative angular displacement ϕ = ϕ1−ϕ2 (3.1) and the reducedmoment of inertia Iz = I1I2 I1+ I2 (3.2) we can transform equations of motion (2.1) into the following form ϕ̈+ M(ϕ,A,ϕ̇) Iz = 1 I1 [M(t)+Mm +Ma] (3.3) Let the solution to equation (3.3) be approximated by ϕ = Acosz z = ωt+ϕ0 (3.4) where: z – forcing phase, ϕ0 – initial forcing phase, A, ϕ0 – are slowly varying functions of time t. Then ϕ̇ = Ȧcosz−Aϕ̇0 sinz −Aω sinz (3.5) By analogy to Lagrange’smethod of variation of a parameter, it is permissible to set Ȧcosz−Aϕ̇0 sinz =0 (3.6) Thus ϕ̈ =−Ȧω sinz −Aω2cosz −Aωϕ̇0 cosz (3.7) Substituting equation (3.7) into the equation of motion (3.3), using formula (3.4)2, gives −Ȧω sinz−Aω2cosz−Aωϕ̇0cosz+ M(ϕ,A,ϕ̇) Iz = Mm +Ma I1 + M0 I1 cos(z−ϕ0) (3.8) Bymultiplying equation (3.6) by ωcosz, equation (3.8) by sinz and substrac- ting the sides while using formula (3.4)2, we obtain −Ȧω−Aω2 sinzcosz + M(ϕ,A,ϕ̇) Iz sinz = (3.9) = Mm +Ma I1 sinz + M0 I1 sinzcos(z −ϕ0) 160 Z.Skup Since A and ϕ0 are slowly varying parameters in equation (3.3), equation (3.9) takes, after integrating over the interval (0,2π), the following form −2πȦω+ 1 Iz 2π ∫ 0 M(ϕ,A,ϕ̇)sinz dz = 2π ∫ 0 Mm +Ma I1 sinz dz + M0π sinϕ0 I1 (3.10) Multiplying equation (3.6) by ω sinz, Eq (3.8) by cosz, adding the sides, using formula (3.4)2 and averaging over one cycle of z, gives −2πAϕ̇0ω−πAω 2+ 1 Iz 2π ∫ 0 M(ϕ,A,ϕ̇)cosz dz = (3.11) = 2π ∫ 0 Mm +Ma I1 cosz dz + M0πcosϕ0 I1 Steady-state equations (3.10) and (3.11) can be obtained when Ȧ = ϕ̇0 = 0, therefore these equations are reduced to the form sinϕ0 = 1 πβM0 2π ∫ 0 M(ϕ,A,ϕ̇)sinz dz − 1 πM0 2π ∫ 0 (Mm +Ma)sinz dz (3.12) Izω 2+ βM0 A cosϕ0 = 1 πA 2π ∫ 0 M(ϕ,A,ϕ̇)cosz dz − β πA 2π ∫ 0 (Mm +Ma)cosz dz In accordance with the Ritz method, the integral with equation (3.3) as the integrand must be equal zero. Therefore 2π ∫ 0 [ ϕ̈+ M(ϕ,A,ϕ̇) Iz − Mm I1 − Mm +Ma I1 − M0cos(z −ϕ0) I1 ] ∂ϕ dt =0 (3.13) where the variation ∂ϕ is equal ∂ϕ = ∂Acosz +∂ϕ0Asinz (3.14) Substituting equation (3.4) and (3.14) into (3.13) and basing on the linear independence of the variations ∂A and ∂ϕ0, we can obtain two independent equations for the displacement amplitude A and phase shift angle ϕ0. We obtain a result which is identical with equations (3.12). Active and passive damping of vibrations... 161 While integrating Eqs (3.12) there is a discontinuity of M(ϕ,A,ϕ̇) and Ma encounterred for ϕ̇ =0. To avoid this, we confine our considerations to a single half-period (motion between two stops). Thus, the integration interval (from 0 to 2π) of the right-hand terms of the above equations is divided into two sub-intervals, from 0 to π for negative dϕ/dt and from π to 2π for positive dϕ/dt. This is, for instance, theprocedure adopted by Giergiel (1990), Osiński (1986), Osiński and Kosior (1976), Skup (1991, 2002). Therefore, substitution of formulas (2.3) into equations (3.12), and subse- quent integration gives after some transformations sinϕ0 = 1 πβM0 ( π ∫ 0 M(ϕ,A,ϕ̇)sinz dz ∣ ∣ ∣ sgn ϕ̇<0 + + 2π ∫ π M(ϕ,A,ϕ̇)sinz dz ∣ ∣ ∣ sgn ϕ̇>0 ) + − 1 πM0 ( π ∫ 0 (Mm +Ma)sinz dz ∣ ∣ ∣ sgn ϕ̇<0 + 2π ∫ π (Mm +Ma)sinz dz ∣ ∣ ∣ sgn ϕ̇>0 ) = = 2k πM0 { 1 βκ1 [√ x+ √ y+2+ 2 3κ1A ( √ y3− √ x3) ] −κ2( √ x+ √ y−2) } x =1+2κ1A y =1−2κ1A (3.15) Izω 2+ βM0 A cosϕ0 = 1 πA ( π ∫ 0 M(ϕ,A,ϕ̇)cosz dz ∣ ∣ ∣ sgn ϕ̇<0 + + 2π ∫ π M(ϕ,A,ϕ̇)cosz dz ∣ ∣ ∣ sgn ϕ̇>0 ) + − β πA ( π ∫ 0 (Mm +Ma)cosz dz ∣ ∣ ∣ sgn ϕ̇<0 + 2π ∫ π (Mm +Ma)cosz dz ∣ ∣ ∣ sgn ϕ̇>0 ) = = k ( 1+ 15 32 κ21A 2 ) Introducing the notations: a – dimensionless vibration amplitude, a = A/ϕst ϕst – static displacement in the form of the relative angular displa- cement of the damper plates, ϕst = M0/k 162 Z.Skup γ – dimensionless frequency, γ = ω/ω0 ω0 – frequency of the free vibration of the system, ω0 = √ k/Iz we find sinϕ0 = 2k πM0 { 1 βκ1 [ η1+2+ 2 3aψ ( √ y31 − √ x31) ] −κ2(η1−2) } (3.16) γ2+ β a cosϕ0 =1+ 15 32 a2ψ2 where ψ = M0κ1 k κ1A = aψ x1 =1+2aψ y1 =1−2aψ η1 = √ x1+ √ y1 (3.17) Finally, inbasic equations (3.16)wehave the relations for the tangent of the phase angle ϕ0 and the dimensionless frequency γ as functions of the external load, geometric and electric parameters, friction forces and the dimensionless amplitude a. Therefore tanϕ0 = 2kβ { 1 βκ1 [ η1+2+ 2 3aψ ( √ y31 − √ x31) ] −κ2(η1−2) } πM0a ( 1+ 15 32 a2ψ2−γ2 ) (3.18) γ = √ 1+ 15 32 a2ψ2∓ β a √ 1− sin2ϕ0 4. Numerical results The following data has been taken for numerical calculations: h1 = 0.004m, h2 = 0.010m, R = 0.070m, D = 0.070m, M0 = 20Nm, µ = 0.25, d = 0.054m, I1 = 0.560kgm 2, I2 = 0.04kgm 2, l = 0.65m, p = 0.8 · 105N/m2, la = 0.04m, lc = 0.04m, G = 8.2 · 10 10N/m2, Ga = 6.3 · 10 9N/m2, Gs = 2 · 10 9N/m2, da15 = 5.6 · 10 −10m/V – for PZT (PIC255), ds15 =0.23·10 −10m/V–forPVDF, εs =12, ε0 =0.088·10 −10F/m, kp =1.0, a =0.25. On the basis of the results of numerical analysis it has been found that all resonance curves start from zero, pass through a resonance and tend again asymptotically to zero in the superresonance range (Fig.2 to Fig.4). They Active and passive damping of vibrations... 163 also becomemore steep in that range. In Fig.2a we can see that the growing amplitude M0 leads to a decrease in the dimensionless amplitudes a. It obvio- usly results from the fact that increment in M0 entails enlargement of the slip zone, thus the amount of the dissipated energy grows. Fig. 2. Resonant curves for various values: (a) of the excitation torque amplitude M0, (b) of the friction coefficient µ Fig. 3. Resonant curves for various values: (a) of the unit pressure p, (b) of the external radius R 164 Z.Skup Fig. 4. Resonant curves for various values: (a) of the discs and plunger equivalent rigidity kz, (b) of the gain factor kp An extreme damping intensity can be observed for a particular value of the friction forces q = pµ (the largest zone of the slip between the discs and plunger surfaces). Figures 2a, 3a disclose a clear effect of the changing of va- lues of µ and p on the nondimensional resonant amplitudes a (all the other parameters remain fixed). The response curves for the two-degree-of-freedom power transmission system reduced to a one-degree-of-freedom non-linear hy- steretic system are reveal a typical ”soft” resonance (Fig.2 to Fig.4). For the excitation frequency ω close to thenatural frequencyof vibrations ω0, thenon- dimensional amplitudes a assume big values. As the presented graphs (Fig.2 to Fig.4) and numerical calculations show, themost dangerous frequency ran- ge for the real values of the parameters of the damper is 0.98 < γ < 1.02. The effect of damping is also best for a suitable value of the external radius R and a reduced discs rigidity kz because of the zone of the biggest slide of the discs in the plunger. The graphs in Fig.3b andFig.4a can serve as an illustration of this conclusion as they show that with the same geometric and electric parameters as well as loading and friction coefficient but varied radius R and reduced rigidity kz, their influenceon the resonance amplitude a is clearly visible. The greater amplifier gain kp is the more visible the effect (Fig.4b). The phase shift angle ϕ0 is a measure of the vibration damping in a me- chanical system. If it increases, then the energy dissipation also increases and, consequently, so does the damping. The dependence of the phase shift an- gle ϕ0 upon the dimensionless frequency γ for different values of the external radius R is represented by diagrams in Fig.5a and, for different values of the Active and passive damping of vibrations... 165 friction coefficient µ, in Fig.5b. A clear influence of the external radius R and friction coefficient µ on the angle ϕ0 can be observed there. The energy dissipation is the largest for big values of R (Fig.5a) and µ (Fig.5b) (big slide zone). Fig. 5. Pphase displacement angle ϕ0 as a function of the dimensionless frequency γ for various values: (a) of the external radius R, (b) of the friction coefficient µ 5. Concluding remarks Structural friction between contacting surfaces of the discs and plunger, causes increased performance of the examined system as for as vibration dam- ping is concerned.Furtherdevelopment of thedamping level is achievedbyma- king use of piezoelectric materials controlled by an electronic circuit. Proper- ties of piezoelectric materials can be utilized in torsional systems for vibration control as piezoceramics constitute perfect elements for actuator applications. Active control by piezoelectric elements and structural friction proves to be a powerful tool in reducing the vibration amplitude of torsional systems. Application of the electronic damping gives excellent results, especially in sys- tems with velocity feedback – even for different points of actuator application and spatial sensor/actuator dislocation. This is highly important since real technical conditions may not always allow arbitrary piezoelements applica- tion. The paper emphasizes some aspects related to the active control under more realistic conditions. Admittelly, it should be said that the vibration damping by friction dam- pers is considerably influenced by the following factors: forcing amplitude, 166 Z.Skup stiffness of the discs and plunger, unit pressure, friction coefficient and ga- in. The examined system has ”soft” frequency characteristic and attenuation diagram. Acknowlegement This paper has been dane within the framework of the KBN Grant No. 5T07C04523 supported by the State Committee for Scientific Research. References 1. Chia-Chi Sung, Vasundara V. Varadan, Xiao-Qi Bao, Vijay K. Va- radan, 1994, Active torsional vibration control experiments using shear-type piezoceramic sensors and actuators, Journal of Intelligent Material Structures and Systems, 5, 3, 436-442 2. Giergiel J., 1990,Damping of Mechanical Vibrations, PWN, Cracow 3. Grudziński K., Kissing W., Zapłata L., 1992, Untersuchung Selbstterreg- ter Reibungsschwingungenmit Hilfe eines Numerischen Simulationsverfahrens, Technische Mechanik, 13, 7-14 4. Kurnik W., Przybyłowicz P.M., 1995, Torsional vibration of a tube with piezoelectric actuators, Zeitschrift für Angewandte Matematik und Mechanik, 75, 1, S55-S56 5. Lee C.K., 1990, Theory of laminated piezoelectric plates for the esign of di- stributed sensors/actuators. Part I: Governing equations and reciprocal rela- tionships, Journal of Acoustical Society of America, 87, 3, 1144-1158 6. Meng-KaoYeh,CHih-YuanChin, 1994,Dynamic response of circular shaft withpiezoelectric sensor,Journal of InteligentMaterial Structures andSystems, 5, 11, 833-840 7. Osiński Z., 1978,Teoria drgań, PWN,Warszawa 8. Osiński Z., 1986,Tłumienie drgań mechanicznych, PWN,Warszawa 9. Osiński Z., 1998, Damping of Vibrations, A.A. BALKEMA/Rotter- dam/Brookfield 10. Osiński Z., Kosior A., 1976, Wpływ tarcia konstrukcyjnego na krzywe re- zonansowe układu przy wymuszeniu harmonicznym,X Konferencja Dynamiki Maszyn, JPBM Politechniki Warszawskiej, 185-195 11. PrzybyłowiczP.M., 1995,Torsional vibration control by active piezoelectric system, Journal of Theoretical and Applied Mechanics, 33, 4, 809-823 Active and passive damping of vibrations... 167 12. Skup Z., 1976, Analiza tarciowego tłumika drgań skrętnych poddanego rów- nomiernemu i nierównomiernemu rozkładowi nacisków z uwzględnieniem tarcia konstrukcyjnego, Rozprawa doktorska, PolitechnikaWarszawska 13. Skup Z., 1987,Drgania swobodne układu ze sprzęgłem ciernymprzy uwzględ- nieniu tarcia konstrukcyjnego,Mechanika, Wyd. AGH, Kraków, 6, 2, 105-122 14. SkupZ., 1991,Forcedvibrationof anonlinear stationary systemwith a friction clutch, structural friction, being taken into consideration, Engineering Trans- actions, 39, 69-87 15. Skup Z., 2002, Active and passive damping of vibrations through torsional damper during starting power transmission system,ArchiwumBudowyMaszyn (in print) 16. Szadkowski A.,MorfordR.B., 1992,Clutch engagement simulation: enga- gement without throttle, SAE Technical Paper Series, 920766, 103-116 17. Tzou H.S., 1991, Distributed piezoelectric neurons andmuscles for shell con- tinua, The 1991 ASME Design Technical Conferences – 13th Biennial Confe- rence on Mechanical Vibration and Noise, Structural Vibration and Acoustics, DE-Vol. 34, 1-6 18. Zagrodzki P., 1994, The simulation of transient processes in hydromecha- nical power transmission systems of vehicles and working machines, (original title in Polish: Symulacja procesówprzejściowychwhydromechanicznychukła- dach napędowych pojazdów i maszyn roboczych), Works of VIII Symposium. Symulation of Dynamical Processes, WarsawDepartment PTETiS, 605-612 Aktywne i pasywne tłumienie drgań poprzez tłłumik drgań skrętnych podczas ruchu ustalonego układu napędowego Streszczenie Praca dotyczy teoretycznego badania tłumienia drgań nieliniowych układu me- chanicznego o dwóch stopniach swobody zawierającego tłumik drgań skrętnych. Roz- ważany jest ruch ustalony układu przywymuszeniu harmonicznymz jednostajną czę- stością o stałej amplitudziemomentuwymuszającego.Uwzględniono tłumienie drgań jednocześnie poprzez zjawisko tarcia konstrukcyjnego (tłumienie pasywne) i elementy piezoelektryczne (tłumienie aktywne). Zagadnienie rozpatrywane jest przy założeniu równomiernego rozkładunacisków jednostkowychwystępującychpomiędzywspółpra- cującymi powierzchniami tarcz ciernych i bezwładnika. Zbadano wpływ parametrów geometrycznych układu, nacisku jednostkowego, obciążenia zewnętrznego oraz pa- rametrów elektrycznych na krzywe rezonansowe drgań ustalonych. Równania ruchu badanego układumechanicznego rozwiązanometodą Van der Pola. Manuscript received June 26, 2002; accepted for print October 5, 2002