Hindawi Publishing Corporation Advances in Acoustics and Vibration Volume 2012, Article ID 659398, 8 pages doi:10.1155/2012/659398

Research Article

Dynamic Analysis of Stockbridge Damper

Nilson Barbieri1,2 and Renato Barbieri1

1 Mechanial Engineering Department, Pontificia Universidade Católica do Parana (PUCPR), Rua Imaculada Conceigao, 1155, 80215-901 Curitiba, PR, Brazil

2 Mechanial Engineering Department, Universidade Tecnológica Federal do Parana (UTFPR), Rua Sete de Setembro, 3165, 80230-901 Curitiba, PR, Brazil

Correspondence should be addressed to Nilson Barbieri, nilson.barbieri@pucpr.br

Received 9 April 2012; Accepted 18 May 2012

Academic Editor: Joseph C. S. Lai

Copyright © 2012 N. Barbieri and R. Barbieri. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In this work a machine cam with five different profiles was used to investigate the linear and nonlinear dynamical behavior of asymmetric Stockbridge damper with excitation frequencies in the range of 5 < f < 17 Hz. The experimental vibration signals were acquired through accelerometers placed along the sample. The loss factor and the Young's modulus were estimated through approximation of the experimental and numerical results using Genetic Algorithms (GAs). Linear and nonlinear mathematical models were used to adjust the data. The two parameters are dependent on the excitation frequency and the amplitude of the base displacement. The results are validated comparing typical impedance curves obtained in conventional testing using an electromechanical shaker.

1. Introduction

The Stockbridge damper is presently the most common type of transmission line damper. In general, the absorber consists of two weights attached to the end of stranded cables, which are known as messenger wires.

A detailed mathematical description of conductor motion is difficult due to the stranded construction of a conductor. An example of this problem is the study carried out by Nawrocki and Labrosse [1], where the cable is modeled using each individual wire model and all possible contacts are investigated. However, in order to get good results for the static analysis, this model was not applied for dynamic problems and the dynamic friction between the individual wires of the cable was also not studied.

In Stockbridge dampers of transmission line, mechanical energy is dissipated in wire cables "damper or messenger cables". The damping mechanism is due to statical hysteresis resulting from Coulomb (dry) friction between the individual wires of the cable undergoing bending deformation. In order to test this dynamical model of Stockbridge damper, the typical experimental impedance curves are compared

with numeric results. This procedure was used by [2] who uses the masing model for modeling the nonlinear damping behavior of the damper cable of the Stockbridge damper.

Zhu and Meguid [3-5] analyzed the behavior of cables considering curved beam models. Due to difficulties in modeling the hysteretic damping, the standard Rayleigh damping was used. The numerical and experimental results presented good approximation.

Lopez and Venegas [6] analyzed the dynamic behavior of Stockbridge dampers through dynamic and fatigue tests and showed that the damping ratio decreases linearly with the excitation motion amplitude.

In this work, the physical parameters are adjusted through the comparison between numerical and experimental results. The numerical models are obtained using the Finite Element Method. The experimental results are obtained using a cam machine with five different eccentricities. This way, the motion amplitude of the excitation is maintained constant. The experimental results are compared with linear and nonlinear numerical results. The nonlinear system contains nonlinear stiffness and damping elements. The dynamic responses of the linear and nonlinear systems are

Figure 1: References and admissible displacements.

also compared with the experimental results (impedance curves) obtained through conventional testing using an electromechanical shaker to excite the system.

2. Mathematical Models

The mathematical models of the Stockbridge damper system are shown in this section.

2.1. Messenger Wire Model. The messenger wire is modeled by the Euler-Bernoulli beam finite element. In this element, the transversal displacement is interpolated using the well-known Hermitian interpolation polynomials with C1 continuity, and the degrees of freedom (d.o.f) in each node are the transversal displacement and the rotation, {v, 9]. The dynamic equation for this element can be written as:

" 156 22L 54 -13L" '>1"

pAL 4L2 13L -3L2 01

420 sym. 156 -22L 4L2 >2 A

12 6L -12 6L >1 6

EI 4L2 -6L 2L2 1 pAL L

L3 12 -6L >2 12~" 6

sym. 4L2 2 -L

where pA is the linear density of the cable, L is the finite element length, EI is the cable flexural (bending) stiffness and yo is the base acceleration (the acceleration imposed by the cam machine).

In order to take into consideration the cable hysteretic damping in (1), it is enough to consider the flexural stiffness as [6];

EI = EIo (1 + nO, (2)

where n is the hysteretic damping constant and i =

2.2. Damper Mass Model. The suspended masses of the Stockbridge damper are modeled with a rigid body plane motion hypothesis and the admissible displacements are shown in Figure 1.

After assembling all the elements of the messenger wire, each weight of the Stockbridge damper contributes to two terms of the dynamical equilibrium. The first contribution is in the mass matrix (inertia force)

[Ms]{ qn} =

m mx mx In

and another parcel is in the vector force due to base acceleration

I f\ = - \m-

us' I mx

,x\yo>

where m is the mass of the Stockbridge damper weight, x is the center of mass coordinate and In is the inertia moment with the reference fixed in node n.

These two terms are obtained using the first variation of the kinetic energy (Hamilton Principles) and the rigid body plane motion hypothesis for modeling the suspended Stockbridge damper weight. With the convention defined in Figure 1, the kinetic energy of each weight of the Stockbridge damper can be written by taking node n (the node of the finite element mesh connected to the messenger wire) as a reference. This expression is

1 - 1 f

T = -mrn ■ rn + rn ■ [w X mx] + -w ■ pn X [w X pn] dm,

2 2 J m

where rn is the velocity of node n, w is the damper weight angular velocity and pn is the particle position of dm mass with the origin fixed in node n.

Taking into consideration the hypotheses of plane motion of rigid body, w = 9nk, and symmetry in y, (y = 0), the expression for the damper weight kinetic energy can be rewritten after integration as

T = - m [ùl + 1>2+ 2 yovn + >0]

[dn>n + Önyo] + 2 InOl.

The variation of T can be written as ST = dT/dün8ün + dT/dvnSvn + dT/d9nS9n. Integrating by parts the terms in Svn and S9n and neglecting the boundary terms, it is found in (3) and (4).

2.3. Finite-Element System of Equations. The Stockbridge damper dynamical equilibrium equation is obtained after assembling all finite elements and it can be conventionally written as

[M]{ q} + [K H q} = {/o} m>

where [M] and [K] are the mass and stiffness matrices, respectively, and {fo} is the force vector. The vector components {q} are the finite-element node displacements and rotations, v and 9 and yo is the acceleration in node 1 (base shaker acceleration).

Considering the base excitation as harmonic, yo(t) = \yo\eiwt, the solution q(t) becomes q(t) = qoeiwt. Substituting this supposed solution into the motion equation (7) yields

[-œ2[M] + [K]] {q0} = |y0|{/0}.

The amplitude of the displacement vector is calculated solving Equation (8) for each frequency w and the amplitude of the acceleration vector is easily calculated with the product of w2{qo}.

2.4. The Nonlinear Formulation. The parameter estimation of the nonlinear system can be made by approximating the numeric and experimental Frequency Response Function (FRF) curves. The analysis can be exemplified considering the motion equation of a simple oscillator subjected to a harmonic excitation [7-11] (Harmonic Balance Method):

my + g( y, y) = / sin œt,

where g(y, y) encloses all the restoring forces and supposing it is a nonlinear function of the velocity and displacement (y, y) of the mass m.

It is possible to find a linearized coefficient v(y, y) which provides the best average of the true restoring force. This coefficient acts on the fundamental harmonic of the nonlinear response (Y1st) for a single load cycle, in such way that [8]

g(y,y) ~ V(y,y) ■ y, for y « Y1st sin(wt + 9) = Y 1stsin t.

In order to find the nonlinear coefficient v(y, y), the restoring force g(y, y) is expanded around y via a Fourier series, neglecting all the higher order terms:

g{y,y) ~ v(y,y) ■ y = O^y + olsty

+ Oclsty + O^y + ■■■ where the o functions are given by

1 ct J-

nY1st Jo

j g^y1st sin t, œ71st cos ^ sin t dr, (12)

oist =

nY1st Jo

g( Y 1stsin r, œY 1stcos ^ cos r dr. (13)

The nonlinear coefficient v(y, y) is defined by (using only the first two terms of (11)):

y) = oist +

The mathematical model of a cubic stiffness element can be expressed as [8]:

&(y> y) = ky + ßy3, (15)

where the coefficient k represents the linear component of the spring and the coefficient ß accounts for the nonlinear effects due to the term y3. Introducing (15) into (12) and (13) and dropping the superscript first and after many manipulations, it can be found that the first-order representation of a cubic stiffness element is [8]:

Vk{y, y) = k +

where the second term on the right side of (16) represents the nonlinear part of the coefficient.

The nonlinear friction damping can be obtained using a similar approach to the cubic stiffness development. The nonlinear restoring force becomes

ïd(y, y) = cy + yjyj ■ (17)

The linearized coefficient vc(y, y) is found to be

vc(y, y) = iwc + i—^. (18)

The g(y, y) restoring forces presents in (9) involves the sum of (15) and (17).

Expanding the idea of the simple oscillator introduced in (9) to a MDOF system, we have

[M]{y} + |g( y, y)j = {F}eiœt,

where [M] is the mass matrix; {y}, {y} and {y} are the acceleration, velocity and displacement vectors, respectively, and {F} is the harmonic excitation vector operating at frequency w.

For additive nonlinearities, it is possible to expand the nonlinear vector into individual nonlinear restoring forces, as follows:

{G( y, y) =

g11 + g12 + g13 + ■ ■ ■ + g1N g21 + g22 + g23 + ■■■ + g2N g31 + g32 + g33 + ■ ■ ■ + g3N

gN1 + gN2 + gN 3 + ■■■ + gNN

where N is the size of the system (in DOFs). Each nonlinear function gij represents a restoring force acting between DOFs i and j, while terms with repeated indexes gu represent a restoring force between DOF i and the ground.

Introducing the newly redefined nonlinear coefficients into (20), a matrix of nonlinear coefficients is formed:

{GO, y)j = [V]{ Yj.

Accelerometer

Signal analyzer

Oscillation with controlled amplitude

-1 4/A2 ■r- ■ I A1 \ 1-

c-1 T j i-' ;

Cam machine Connecting rod

Figure 2: Schematic cam machine.

The motion equation of a general nonlinear system subjected to harmonic excitation can be described by the following nonlinear ordinary differential equation:

[M]{ y} + [C]{ y} + i[D]{ y}

+ [K ]{y} + {G( y, y)} = {F je^,

where [M], [C], [D], and [K] are the mass, viscous damping, hysteretic damping, and stiffness matrices, respectively, ofthe underlying linear system; {y], {y], and {y] are the acceleration, velocity and displacement vectors, respectively, and {F] is the harmonic excitation vector operating at frequency w.

The nonlinear component of the system is represented by the nonlinear vector{G] , which is a function of all displacements and velocities in the general case.

Considering a harmonic response {y(i)] = {7 ]eiwt, where {7 ] = {17 \ei 9 ] is a nonlinear complex vector allowing it to accommodate phase, the motion equation is further reduced to

(-w2[M] + iw[C] + i[D] + [K]){7} + {G(w, 7)} = {F].

The linear receptance can be defined as

[a] = (-w2[M]+ iw[C] + i[D] + [K ])-1 (24) and its inverse, [A] = [a]-1, as

[A] = -w2[M] + iw[C] + i[D] + [K ]. (25) Inserting (16) and (20) into (18), we have

([A]+ [?]){ 7j = {F ] (26)

leading to the final compact representation of the nonlinear function:

[A] {7} = {F], where: [A] = [A] + [v], (27)

where [A] is a composite matrix, enclosing linear and nonlinear coefficients and it is formulated for the current state [y,y, w]. It can be considered the system matrix. Obtaining

the determinant roots of [A] yields the frequency-dependent nonlinear natural frequencies and damping ratios. The system response {7] can be obtained solving (27).

The solution of (27) can be used to adjust and/or update the parameters in numerical models obtained through the finite-Element method by comparing numerical and experimental results [12]. Satisfactory results were obtained using the genetic algorithms method [13] to make the approximation of the numeric and experimental values. The GA's parameters used in this application are mutation rate =0.02, population size =50, and number of generations =5000. The objective function is defined by

f = ^ Paa,

where Paaexp is the acceleration ratio (acceleration of the reference sensor/acceleration of the accelerometer on the cable); PaaFEM is the same acceleration ratio estimated using GA; np is the number of points (normally np = 47).

3. Results

Figure 2 shows the schematic representation of the machine cam used in the experimental testing of the Stockbridge damper with controlled oscillation (displacement). The experimental data are obtained through three accelerometers placed in the half sample. One accelerometer was placed in the center (A2) and another at the end of the messenger wire (A3). A third accelerometer (reference) was placed on the connecting rod (A1).

Figure 3 shows the cam machine with Stockbridge damper coupled and two accelerometers positioned on the weight. The data acquisition system is composed of a signal analyzer (HP 3567) and modal accelerometers (PCB 333AX, sensitivity «100mv/g). The cam machine with radial flat-faced follower is manually operated.

Five different disk cams with eccentricities of 0.25, 0.5, 0.75, 1.25, and 1.5 mm were used. The tests were carried out varying the excitation frequency between 5 and 17 Hz with increments of 0.25 Hz. This lower frequency range was used due to the mechanical limitation of the excitation system. Figures 4 and 5 show the experimental curves of the acceleration ratio obtained for the accelerometers placed in

Figure 3: Cam machine with Stockbridge damper coupled.

4 3.5 3 2.5 2 1.5 1

10 12 Frequency (Hz)

- 1.25

---0.25

...... 0.5

...... 0.75

10 12 14 Frequency (Hz)

- 1.25

Figure 4: Acceleration ratio curves of accelerometer 1.

Figure 5: Acceleration ratio curves of accelerometer 2.

the center (accelerometer 1) and at the end (accelerometer 2) of the messenger wire. The natural frequency changes with the motion amplitude can be noticed. The frequency variation is inversely proportional to the motion amplitude. Figures 6(a) and 6(b) show the variation of the natural frequency and the damping ratio (adjusted with the linear system) versus the amplitude of the base excitation. It can be noted that parameters, the natural frequency, and the damping ratio vary inversely with the increase of the base motion amplitude. Lopez and Venegas [6] found slightly different results in experimental tests.

The linear and nonlinear parameters used in the numerical models were adjusted considering five different cam profiles. The results are shown in Table 1. The linear parameters adjusted are the Young's modulus E and the loss factor n according to (2).

To simulate the hysteretic damping of the system, it was considered n = nt/w (frequency dependent parameter). For

the nonlinear system, these two parameters and the other nonlinear parameters of stiffness (^ and ¡32) and damping (y1 and y2), described in (16) and (18) using the procedure described in [8] and only considering the main terms of the nonlinear stiffness and damping matrixes (main diagonal), were adjusted. The error column (Table 1) represents the sum of the absolute errors (difference between the linear and nonlinear FRF) obtained for a limited range of frequencies (approximately 2 Hz in the resonance region) divided by the number of points.

Figure 7(a) shows the reference base displacement curve obtained in a conventional testing using an electromechanical shaker for excitation. The base displacement variation with the excitation frequency can be noted. The real acceleration measured at Stockbridge damper was considered as input for the mathematical system (Figure 7(b)).

Figures 8(a) and 8(b) show the fitted curves of EI (bending stiffness) and n (loss factor) as a function of the base

0.6 0.8 1 1.2 Displacement (mm)

Accelerometer 1 Accelerometer 2

& 0.09 m

0.6 0.8 1 1.2 Displacement (mm)

Accelerometer 1 Accelerometer 2

Figure 6: Curves of natural frequencies and damping ratio versus base displacement.

30 40 50 Frequency (Hz) (a)

50 45 40

iT 35 (m/s 30

tio 25 tar

lere 20 c

ccA15 10 5 0

30 40 50 Frequency (Hz)

Figure 7: Base displacement and acceleration. Table 1: Parameters adjusted for the linear and nonlinear systems

Eccentricity (mm) Absolute error EI (Nm2) n> ß1 (N/m3) ß2 (N/m3) Y1 (N) Y2(N)

0.25 linear nonlinear 0.589 0.600 1.90 1.98 9.82 11.46 -7.9e9 -0.0036 -0.0264 -0.0502

0.5 linear 0.454 1.74 6.21

nonlinear 0.480 1.69 7.68 -9.59e8 -0.0308 -0.0416 -0.0964

0.75 linear nonlinear 0.496 0.464 1.65 1.67 4.40 6.57 -4.59e8 -0.0077 -0.0445 -0.0084

1.25 linear nonlinear 0.568 0.569 1.52 1.54 2.94 3.72 -1.67e8 -0.0182 -0.0490 -0.0537

1.5 linear 0.828 1.43 2.56

nonlinear 0.797 1.08 11.24 -5.20e8 -0.1601 -0.0147 -0.0996

(a) (b)

Figure 8: Curve of the flexural stiffness and loss factor versus base displacement.

45 40 35 2 30

^ 15 10 5 0

30 40 50 Frequency (Hz)

—•— Numeric - Experimental

Figure 9: Cable acceleration.

8 101 n

0 10 20 30 40 50 60 70 Frequency (Hz)

- Experimental

° Numeric

Figure 10: Impedance curves.

motion amplitude. These curves are used as a reference for adjusting the Stockbridge damper parameters tested in an electromechanical shaker.

Figure 9 shows the experimental and adjusted cable acceleration curves using a linear numerical model. The curves present good agreement for f <= 10 Hz (region of the first vibration mode). As the frequency testing is lower than 17 Hz, it is not possible to extrapolate the results for f > 17 Hz.

Figure 10 shows a typical experimentally measured impedance curve and the corresponding numeric curve fitted (linear model) with the parameters shown in Figure 8. It is possible to notice the good agreement between the numeric and experimental curves.

4. Conclusions

The results showed that the natural frequency of the system varies with the amplitude of the excitation motion, that is, the higher the motion amplitude is the lower the natural frequency gets. The same behavior was verified for the damping ratio of the system, that is, the higher the motion amplitude is the lower the damping ratio gets.

It was verified that the real and imaginary parts of the complex Young modulus decrease with the motion amplitude increase when the parameters are fitted for a linear system.

For the frequencies range 5 < f < 17 Hz and the motion displacement smaller than 3 mm, the system response with

nonlinear parameters did not show improvements in relation to the linear system response.

As the test was done with a full Stockbridge damper, it is possible that the smallest cable and the smallest mass have influenced the results. The mathematical model was based only on the Stockbridge damper half model containing the largest weight.

When the tests are performed with constant base motion amplitude, it is possible to fit a single value of the parameters according to the motion amplitude.

When the tests are carried out using an electromechanical shaker, the motion amplitude is varied according to the excitation frequency. In this case, the fitted parameters are frequency dependent.

The typical numerical and experimental impedance curves present good agreement.

Acknowledgment

The authors gratefully acknowledge the support for this work provided by the Brazilian Science Foundation CNPq.

References

[1] A. Nawrocki and M. Labrosse, "A finite element model for simple straight wire rope strands," Computers and Structures, vol. 77, no. 4, pp. 345-359, 2000.

[2] H. Verma, The Stockbridge damper as a continuous hysteric system in single overhead transmission lines [M.S. thesis], Indian Institute of Technology Bombay, 2002.

[3] Z. H. Zhu and S. A. Meguid, "Elastodynamic analysis of low tension cables using a new curved beam element," International Journal of Solids and Structures, vol. 43, no. 6, pp. 1490-1504, 2006.

[4] Z. H. Zhu and S. A. Meguid, "Modeling and simulation of aerial refueling by finite element method," International Journal of Solids and Structures, vol. 44, no. 24, pp. 8057-8073, 2007.

[5] Z. H. Zhu and S. A. Meguid, "Nonlinear FE-based investigation of flexural damping of slacking wire cables," International Journal of Solids and Structures, vol. 44, no. 16, pp. 5122-5132, 2007.

[6] A. Lara-Lopez and J. C. Venegas, "Endurance of dampers for electric conductors," International Journal of Fatigue, vol. 23, no. 1, pp. 21-28, 2001.

[7] G. Kerschen, K. Worden, A. F. Vakakis, and J. C. Golinval, "Past, present and future of nonlinear system identification in structural dynamics," Mechanical Systems and Signal Processing, vol. 20, no. 3, pp. 505-592, 2006.

[8] H. R. E. Siller, Non-linear modal analysis methods for engineering structures [Ph.D. thesis], Imperial College London/University ofLondon, 2004.

[9] G. Baldo, A. Bonelli, O. Bursi, and S. Erlicher, "The accuracy of the generalized-a method in the time integration of nonlinear single- and two-DOF forced systems," Computational Mechanics, vol. 38, no. 1, pp. 15-31, 2006.

[10] D. Akgun and I. Cankaya, "Frequency response investigations of multi-input multi-output nonlinear systems using automated symbolic harmonic balance method," Nonlinear Dynamics, vol. 61, no. 4, pp. 803-818, 2010.

[11] R. K. Narisetti, M. Ruzzene, and M. J. Leamy, "Study of wave propagation in strongly nonlinear periodic lattices using

a harmonic balance approach," Wave Motion, vol. 49, no. 2, pp. 394-210, 2012.

[12] H. Grafe, Model updating of large structural dynamics models using measured response functions [Ph.D. thesis], Imperial College London/ University of London, 1998.

[13] W. D. Chang, "An improved real-coded genetic algorithm for parameters estimation of nonlinear systems," Mechanical Systems and Signal Processing, vol. 20, no. 1, pp. 236-246, 2006.

Copyright of Advances in Acoustics & Vibration is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.