Research Article | Open Access
Dynamic Analysis of Stockbridge Damper
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 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.
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 , 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  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.
López and Venegas  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 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 continuity, and the degrees of freedom (d.o.f) in each node are the transversal displacement and the rotation, . The dynamic equation for this element can be written as: where is the linear density of the cable, is the finite element length, is the cable flexural (bending) stiffness and is the base acceleration (the acceleration imposed by the cam machine).
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) and another parcel is in the vector force due to base acceleration where is the mass of the Stockbridge damper weight, is the center of mass coordinate and is the inertia moment with the reference fixed in node .
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 (the node of the finite element mesh connected to the messenger wire) as a reference. This expression is where is the velocity of node , is the damper weight angular velocity and is the particle position of mass with the origin fixed in node .
Taking into consideration the hypotheses of plane motion of rigid body, , and symmetry in , (), the expression for the damper weight kinetic energy can be rewritten after integration as
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 where  and  are the mass and stiffness matrices, respectively, and is the force vector. The vector components are the finite-element node displacements and rotations, and and is the acceleration in node 1 (base shaker acceleration).
Considering the base excitation as harmonic, , the solution becomes . Substituting this supposed solution into the motion equation (7) yields
The amplitude of the displacement vector is calculated solving Equation (8) for each frequency and the amplitude of the acceleration vector is easily calculated with the product of .
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): where encloses all the restoring forces and supposing it is a nonlinear function of the velocity and displacement of the mass .
It is possible to find a linearized coefficient which provides the best average of the true restoring force. This coefficient acts on the fundamental harmonic of the nonlinear response for a single load cycle, in such way that 
In order to find the nonlinear coefficient , the restoring force is expanded around via a Fourier series, neglecting all the higher order terms: where the functions are given by
The nonlinear coefficient is defined by (using only the first two terms of (11)):
The mathematical model of a cubic stiffness element can be expressed as : where the coefficient represents the linear component of the spring and the coefficient accounts for the nonlinear effects due to the term . 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 : 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
Expanding the idea of the simple oscillator introduced in (9) to a MDOF system, we have where is the mass matrix; , and are the acceleration, velocity and displacement vectors, respectively, and is the harmonic excitation vector operating at frequency .
For additive nonlinearities, it is possible to expand the nonlinear vector into individual nonlinear restoring forces, as follows: where is the size of the system (in DOFs). Each nonlinear function represents a restoring force acting between DOFs and , while terms with repeated indexes represent a restoring force between DOF and the ground.
Introducing the newly redefined nonlinear coefficients into (20), a matrix of nonlinear coefficients is formed:
The motion equation of a general nonlinear system subjected to harmonic excitation can be described by the following nonlinear ordinary differential equation: where , , , and are the mass, viscous damping, hysteretic damping, and stiffness matrices, respectively, of the underlying linear system; , and are the acceleration, velocity and displacement vectors, respectively, and is the harmonic excitation vector operating at frequency .
The nonlinear component of the system is represented by the nonlinear vector, which is a function of all displacements and velocities in the general case.
Considering a harmonic response, where is a nonlinear complex vector allowing it to accommodate phase, the motion equation is further reduced to
The linear receptance can be defined as and its inverse, , as Inserting (16) and (20) into (18), we have leading to the final compact representation of the nonlinear function: where is a composite matrix, enclosing linear and nonlinear coefficients and it is formulated for the current state . It can be considered the system matrix. Obtaining the determinant roots of yields the frequency-dependent nonlinear natural frequencies and damping ratios. The system response 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 . Satisfactory results were obtained using the genetic algorithms method  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 where is the acceleration ratio (acceleration of the reference sensor/acceleration of the accelerometer on the cable); is the same acceleration ratio estimated using GA; is the number of points (normally ).
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 ≈100 mv/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 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. López and Venegas  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 and the loss factor according to (2).
To simulate the hysteretic damping of the system, it was considered (frequency dependent parameter). For the nonlinear system, these two parameters and the other nonlinear parameters of stiffness ( and ) and damping ( and ), described in (16) and (18) using the procedure described in  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 (bending stiffness) and (loss factor) as a function of the base 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 (region of the first vibration mode). As the frequency testing is lower than 17 Hz, it is not possible to extrapolate the results for .
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.
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 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.
The authors gratefully acknowledge the support for this work provided by the Brazilian Science Foundation CNPq.
- 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.
- H. Verma, The Stockbridge damper as a continuous hysteric system in single overhead transmission lines [M.S. thesis], Indian Institute of Technology Bombay, 2002.
- 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.
- 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.
- 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.
- A. Lara-López and J. C. Venegas, “Endurance of dampers for electric conductors,” International Journal of Fatigue, vol. 23, no. 1, pp. 21–28, 2001.
- 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.
- H. R. E. Siller, Non-linear modal analysis methods for engineering structures [Ph.D. thesis], Imperial College London/University of London, 2004.
- G. Baldo, A. Bonelli, O. Bursi, and S. Erlicher, “The accuracy of the generalized-α method in the time integration of non-linear single- and two-DOF forced systems,” Computational Mechanics, vol. 38, no. 1, pp. 15–31, 2006.
- D. Akgün and I. Çankaya, “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.
- 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.
- H. Grafe, Model updating of large structural dynamics models using measured response functions [Ph.D. thesis], Imperial College London/ University of London, 1998.
- 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 © 2012 Nilson Barbieri and Renato 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.