# Match of Negative Stiffness and Viscous Damping in a Passive Damper for Cable Vibration Control

Authors: Peng Zhou; Qinghe Fang
Journal: Shock and Vibration
 (2019)

---

## Abstract

Match of negative stiffness and viscous damping in a passive negative stiffness damper (NSD) is studied for the vibration control of stay cables in this paper. At first, a discrete model of the stay cable with an NSD attached perpendicularly near the support is established. Under sinusoidal excitations, forced responses of the system are derived theoretically, which results in an asymptotic form for the additional modal damping ratios. Then, experimental results are presented to verify the discrete model and the corresponding theoretical derivations. Subsequently, numerical analysis is performed further to show the optimal match of negative stiffness and viscous damping, which is a function of the attachment location. The energy dissipated by the NSD and the cable energy are analyzed, thereby demonstrating the change trend of the additional modal damping ratios. Moreover, the energy distribution along the cable is investigated to reveal the effect of the negative stiffness and viscous damping. This study demonstrates the control mechanism of negative stiffness and viscous damping in the passive damper and is of practical significance for designing the optimal match of the damper parameters for cable vibration control.


---

## Body

## 1. Introduction

As critical components of cable-stayed bridges, stay cables are usually extremely flexible with low inherent damping, which makes them prone to various unwanted vibrations induced by kinds of external disturbances. For example, the deck motion might interact with the cables to result in nonlinear resonance [1, 2]. These vibrations would significantly affect the safety of stay cables by accelerating the corrosion or destroying the protection. Sometimes, even a pedestrian’s nonconfidence could be caused to interrupt the normal operation of the bridges. In this regard, it turns out to be of importance to mitigate these undesired vibrations of stay cables.To achieve such an aim, several countermeasures have been proposed, including the aerodynamic methods [3, 4], cross-ties [5–9], and mechanical dampers attached near the support [10, 11], among which negative stiffness dampers (NSDs) have been widely studied in recent years because of their superior performance. At the very beginning, the semiactive magnetorheological (MR) damper is found to yield higher effects in controlling cable vibrations [12, 13], which is a result of the negative stiffness behavior in the damping force as observed in in situ experiments [14]. Furthermore, Li et al. [15] analyzed the negative stiffness behavior in active or semiactive dampers and systematically demonstrated its ability to amplify the damper displacement response for dissipating more energy to provide a higher modal damping ratio. Moreover, Weber and Boston [16] investigated clipped viscous damping with negative stiffness numerically and experimentally. When tuned properly, additional modal damping was twice as much as the case of only optimal viscous damping. After that, an approximate applicable collocated method was developed for efficient mitigation of stay cable multimode vibrations based on the control force under clipped linear quadratic regulators [17].The negative stiffness provided by semiactive dampers depends on the real-time feedback and power supply, which makes the entire system complex and unreliable. To draw these drawbacks, passive NSDs were suggested to replace the semiactive ones. For example, prepressed springs were incorporated to produce negative stiffness behavior [18–23]. As shown by the theoretical analysis, more damping was added to stay cables by a self-contained highly compressed spring in parallel with a viscous damper [21]. Furthermore, numerical and experimental results indicated that a negative stiffness damper, which consisted of two pressed springs and an oil damper, provided a larger modal damping ratio to the cable regardless of single-mode and multimode vibrations [22]. Even for the cable network, this kind of NSD exhibited superior damping enhancement for two stay cables simultaneously through interconnecting them together [23]. Magnetic NSDs using the magnetic interaction to produce negative stiffness behavior are also developed and investigated for the cable vibration mitigation [24, 25]. Kalathur and Lakes [26] proposed a postbuckled column NSD, which could achieve high structural damping combined with high initial stiffness at small amplitudes. Lu et al. [27] and Wang et al. [28] utilized the inertial mass to realize negative stiffness for cable vibration control.In the above studies, much attention was paid to the design and realization of passive NSDs. The control mechanism and the optimal match of negative stiffness and viscous damping in the passive damper are seldom mentioned. This paper investigates the optimal match of the negative stiffness and viscous damping, followed by revelation of their control mechanism for stay cables from the perspective of energy. Initially, a discrete model with infinite degrees is established to describe the stay cable with an NSD located close to the support. Steady-state forced responses of the system under harmonic loads are derived, and additional modal damping ratios are presented approximately. Then, experimental tests are presented to compare with the above derivations. After verifying the model, numerical simulation is conducted further to study the relationship between negative stiffness and viscous damping in a passive damper for the optimal control of stay cables. The energy dissipated by the passive NSD and the cable energy, which show the changing law of additional modal damping, are analyzed numerically. Furthermore, the spatial distribution of vibrational energy along the cable is formulated to reveal the effect of negative stiffness and viscous damping. Finally, conclusions from the current study are addressed.


## 2. Finite Degrees of Discrete Model for the Stay Cable with a Passive NSD

### 2.1. Equation of Motion for the System

A discrete model withn+1 masses connected by n springs, as shown in Figure 1, is explored to investigate the dynamic behavior of a stay cable. The mass on two supports is m/2, while on the other place is m. Each spring has the identical length l and tension T in the horizontal condition. The passive NSD with negative stiffness kd and viscous damping cd is located at the node nd close to the left support, while the external excitation Fet acts at the node ne near the right support. Neglecting the gravity and flexural effects, only the transverse vibration yi for the ith node is considered.Figure 1
Discrete model of a stay cable with a passive NSD attached in the transverse direction.For small deflections, three adjacent nodesi−1, i, and i+1 are described in Figure 2. Based on d’Alembert’s principle, the motion equation for the node i 1≤i≤n−1,≠nd,ne in the transverse direction yields(1)my¨i+Tsinθi−1−Tsinθi=0,where sinθi−1 and sinθi are approximated as(2a)sinθi−1≈yi−yi−1l,(2b)sinθi≈yi+1−yil,when yi−yi−1/l≪1 and yi+1−yi/l≪1 are satisfied. Further simplification leads to(3)y¨i+Tlm−yi−1+2yi−yi+1=0.Figure 2
Force and motion diagram for theith node 1<i<n,≠nd,ne.With the damping force, the equation of motion for the nodend changes into(4)y¨nd+Tlm−ynd−1+2ynd−ynd+1=−cdmy˙nd−kdmynd.Similarly, at the nodene, it can be obtained that(5)y¨ne+Tlm−yne−1+2yne−yne+1=Fetm.The nodes0 and n are clamped to the left and the right support, respectively. Hence, the boundary conditions(6a)y0=0,(6b)yn=0,are always satisfied, which results in [29](7)ωr=2Tlmsinrπ2n,(8)uir=sinirπn,where ωr is the rth modal angular frequency of the undamped discrete cable model and uir is the corresponding amplitude of the rth mode at node i.Rewriting in the vector form, the modal shape for therth mode becomes(9)ur=u1ru2r⋯un−1r,which has the orthogonal relation(10)ur⋅us=n−12,r=s,0,r≠s.


### 2.2. Steady-State Responses under Sinusoidal Excitation

Sinusoidal excitation is assumed to act on the nodene as(11)Fet=Fsinωjt,where ωj is the jth modal angular frequency depicted in equation (7) and F is the amplitude of the external excitation.Then, the motion equation for the system is rewritten in the matrix form as(12)y¨+Ky=LdFdt+LeFet,where y=y1y2⋯yn−1T is the displacement vector of n−1 nodes; y¨ is the corresponding acceleration; Ld=0⋯1⋯0T is the damper location vector in which only the nd element is unity; Le=0⋯1⋯0T is the excitation position vector in which the ne element is unity, while others are zero; and K is the stiffness matrix:(13)K=Tlm2−1−12−1−1⋱⋱⋱⋱−1−12,which satisfies the orthogonality relation(14)ur⋅K⋅us=2n−1sin2rπ2n,r=s,0,r≠s.After reaching steady-state vibration, the system response has the following form:(15)y=u1u2⋯un−1a1sinωjt+b1cosωjta2sinωjt+b2cosωjt⋮an−1sinωjt+bn−1cosωjt.Substitution of equation (15) into (12) leads to a linear algebraic equation giving the solution of a1,a2,…,an−1 and b1,b2,…,bn−1 as follows:(16a)ar=2n−11ωr2−ωj21undjundjuner−unejundrFm,r=1,2,…,n−1,r≠j,(16b)aj=2n−11undj2∑r=1,r≠jn−1undrωr2−ωj2unejundr−undjunerFm+unejundj2kdkd2+cd2ωj2F,(16c)br=0,(16d)bj=−unejundj2cdωjkd2+cd2ωj2F.Hence, the response of each node in the discrete model can be obtained. In particular, the displacement of the nodend with the NSD attached is(17)ynd=unejundjFkd2+cd2ωj2sinωjt−φ,in which φ is the phase difference depicted as(18)φ=arctancdωjkd.


### 2.3. Additional Modal Damping Ratios for the Cable

Using the displacement on the nodend, the energy dissipated by the NSD in one period of the jth mode is obtained as(19)ΔEj=∫02π/ωjcdy˙nd+kdynd⋅y˙nddt=unejundj2πcdωjkd2+cd2ωj2F2.Meanwhile, the potential energy of the stay cable in thejth mode is(20)Ej=max12yTKy≈n−14mωj2∑2m2+1kd2+cd2ωj22unejundj2∑mkd+unej2undj4F2,in which ∑=2/n−11/undj2∑r=1,r≠jn−1undr/ωr2−ωj2unejundr−undjuner.Then, thejth modal damping ratio provided by the NSD is(21)ξj=ΔEj4πEj=unejundj2mn−1ωj∑2cdkd+munej/∑undj22+ωj2cd2.Setting the derivative ofξj with respect to cd equal to zero, the modal damping ratio gains its peak value when(22)cdopt=1ωjkd+m∑unejundj2,in which cdopt is the optimal damper coefficient maximizing the modal damping ratio, and the corresponding optimal damper stiffness is(23)kdopt=−m∑unejundj2.


## 2.1. Equation of Motion for the System

A discrete model withn+1 masses connected by n springs, as shown in Figure 1, is explored to investigate the dynamic behavior of a stay cable. The mass on two supports is m/2, while on the other place is m. Each spring has the identical length l and tension T in the horizontal condition. The passive NSD with negative stiffness kd and viscous damping cd is located at the node nd close to the left support, while the external excitation Fet acts at the node ne near the right support. Neglecting the gravity and flexural effects, only the transverse vibration yi for the ith node is considered.Figure 1
Discrete model of a stay cable with a passive NSD attached in the transverse direction.For small deflections, three adjacent nodesi−1, i, and i+1 are described in Figure 2. Based on d’Alembert’s principle, the motion equation for the node i 1≤i≤n−1,≠nd,ne in the transverse direction yields(1)my¨i+Tsinθi−1−Tsinθi=0,where sinθi−1 and sinθi are approximated as(2a)sinθi−1≈yi−yi−1l,(2b)sinθi≈yi+1−yil,when yi−yi−1/l≪1 and yi+1−yi/l≪1 are satisfied. Further simplification leads to(3)y¨i+Tlm−yi−1+2yi−yi+1=0.Figure 2
Force and motion diagram for theith node 1<i<n,≠nd,ne.With the damping force, the equation of motion for the nodend changes into(4)y¨nd+Tlm−ynd−1+2ynd−ynd+1=−cdmy˙nd−kdmynd.Similarly, at the nodene, it can be obtained that(5)y¨ne+Tlm−yne−1+2yne−yne+1=Fetm.The nodes0 and n are clamped to the left and the right support, respectively. Hence, the boundary conditions(6a)y0=0,(6b)yn=0,are always satisfied, which results in [29](7)ωr=2Tlmsinrπ2n,(8)uir=sinirπn,where ωr is the rth modal angular frequency of the undamped discrete cable model and uir is the corresponding amplitude of the rth mode at node i.Rewriting in the vector form, the modal shape for therth mode becomes(9)ur=u1ru2r⋯un−1r,which has the orthogonal relation(10)ur⋅us=n−12,r=s,0,r≠s.


## 2.2. Steady-State Responses under Sinusoidal Excitation

Sinusoidal excitation is assumed to act on the nodene as(11)Fet=Fsinωjt,where ωj is the jth modal angular frequency depicted in equation (7) and F is the amplitude of the external excitation.Then, the motion equation for the system is rewritten in the matrix form as(12)y¨+Ky=LdFdt+LeFet,where y=y1y2⋯yn−1T is the displacement vector of n−1 nodes; y¨ is the corresponding acceleration; Ld=0⋯1⋯0T is the damper location vector in which only the nd element is unity; Le=0⋯1⋯0T is the excitation position vector in which the ne element is unity, while others are zero; and K is the stiffness matrix:(13)K=Tlm2−1−12−1−1⋱⋱⋱⋱−1−12,which satisfies the orthogonality relation(14)ur⋅K⋅us=2n−1sin2rπ2n,r=s,0,r≠s.After reaching steady-state vibration, the system response has the following form:(15)y=u1u2⋯un−1a1sinωjt+b1cosωjta2sinωjt+b2cosωjt⋮an−1sinωjt+bn−1cosωjt.Substitution of equation (15) into (12) leads to a linear algebraic equation giving the solution of a1,a2,…,an−1 and b1,b2,…,bn−1 as follows:(16a)ar=2n−11ωr2−ωj21undjundjuner−unejundrFm,r=1,2,…,n−1,r≠j,(16b)aj=2n−11undj2∑r=1,r≠jn−1undrωr2−ωj2unejundr−undjunerFm+unejundj2kdkd2+cd2ωj2F,(16c)br=0,(16d)bj=−unejundj2cdωjkd2+cd2ωj2F.Hence, the response of each node in the discrete model can be obtained. In particular, the displacement of the nodend with the NSD attached is(17)ynd=unejundjFkd2+cd2ωj2sinωjt−φ,in which φ is the phase difference depicted as(18)φ=arctancdωjkd.


## 2.3. Additional Modal Damping Ratios for the Cable

Using the displacement on the nodend, the energy dissipated by the NSD in one period of the jth mode is obtained as(19)ΔEj=∫02π/ωjcdy˙nd+kdynd⋅y˙nddt=unejundj2πcdωjkd2+cd2ωj2F2.Meanwhile, the potential energy of the stay cable in thejth mode is(20)Ej=max12yTKy≈n−14mωj2∑2m2+1kd2+cd2ωj22unejundj2∑mkd+unej2undj4F2,in which ∑=2/n−11/undj2∑r=1,r≠jn−1undr/ωr2−ωj2unejundr−undjuner.Then, thejth modal damping ratio provided by the NSD is(21)ξj=ΔEj4πEj=unejundj2mn−1ωj∑2cdkd+munej/∑undj22+ωj2cd2.Setting the derivative ofξj with respect to cd equal to zero, the modal damping ratio gains its peak value when(22)cdopt=1ωjkd+m∑unejundj2,in which cdopt is the optimal damper coefficient maximizing the modal damping ratio, and the corresponding optimal damper stiffness is(23)kdopt=−m∑unejundj2.


## 3. Experimental Validation

### 3.1. Equivalent Parameters for the NSD

To suppress the cable vibration, Zhou and Li [22] developed a passive NSD, which is composed of two prepressed springs and an oil damper. The damping force consists of the negative stiffness and the Coulomb friction as(24)fNSD=kdx+Frsgnx˙,where fNSD is the damping force, kd=−2.73N/mm is the negative stiffness approximately determined by the NSD’s geometrical configuration, Fr=17N is the Coulomb friction level, and x and x˙ are the displacement and velocity of the damper, respectively.Making the dissipation energy per cycle the same, an equivalent viscous damping is adopted to replace the Coulomb friction part of the NSD when the cable-damper system vibrates in thejth mode as follows:(25)ceq=4FrπωjAd,in which ceq is the viscous damping coefficient and Ad is the amplitude of displacement at the damper position.Figure3 shows the equivalent hysteretic loops for one experimental case. The areas surrounded by the two curves are the same. Equivalent viscous damping coefficient ceq is inversely proportional to the product of the jth modal angular frequency ωj and the damper’s displacement amplitude Ad when the Coulomb friction force is kept unchanged. In the following verification, analytical results using the equivalent viscous damping coefficients will be compared with the experimental results.Figure 3
Equivalent hysteretic loops for the experimental NSD.


### 3.2. Discretization of the Continuous Cable

The length of the scaled stay cable used in the experiment [22] was 14m with the mass per unit length of 1.33kg/m. The tension force was adjusted to 10kN. Two damper locations were considered in the tests, which are 0.535m and 1.06m to the left anchorage, respectively. An electric shaker was employed to produce sinusoidal excitation forces acting at 0.345m from the right anchorage. For each damper location, the first and the second modal vibration was excited for the cable. After reaching the steady state, the electric shaker was cut off, and then the cable began its free vibration with exponential decaying amplitudes, which are explored to identify the additional modal damping ratios provided by the passive NSD.For the discrete model, the parameters, such as the node massm and the length l, are equal to the total mass and total length divided by the node number n. Then, the natural frequencies are determined by equation (7). Figure 4(a) depicts the frequency ratio of the discrete model to the continuous cable for the first three modes. With the node number n increase, the ratio converges to unity quickly and the relative error decreases dramatically. When the node number n is larger than 100, the relative error becomes less than 0.1%. Hence, n=200 is utilized in the discrete model to guarantee the precision. As shown in Figure 4(b), the first three modal frequencies of the discrete model are almost the same as the experimental results.Figure 4
(a) Frequency ratios for the first three modes; (b) comparison of the natural frequencies withn=200 to the experimental results.
(a)
(b)


### 3.3. Verification for the Displacement Response and Additional Modal Damping Ratio

The excitation force in the experiment is used to calculate the theoretical responses of the discrete model. The excitation position isne=195. Figure 5 compares the first modal displacement responses of the discrete model with experimental results when the damper is located at nd=8. The calculated displacement response at the damper location coincides well with the experimental result. Moreover, the calculated displacement response at the midspan of the cable agrees well with the experimental result.Figure 5
Displacements of the cable under the first modal vibration with the NSD: (a) damper location; (b) midspan.
(a)
(b)Specially, the relationship of displacement amplitude to damping at the damper’s location is shown in Figure6. The vertical axis is the dimensionless displacement amplitude according to equation (17). Under different levels of damping, experimental results are very close to the theoretical curve. The value of dimensionless displacement decreases with increasing damping cdω, as high damping will clamp the cable to make the damper’s location response small.Figure 6
Relationship of displacement amplitude to damping at the damper’s location.The additional modal damping ratios provided by the NSD are demonstrated in Figure7. In general, the experimental damping ratios at two locations agree with the theoretical curves according to the above derivation, although some points deviate a little largely. Comparing with damper stiffness kd=0, more damping ratios can be achieved when kd=−2.73N/mm, which means the negative stiffness is beneficial to enhance damping ratios. However, the beneficial effect varies with the level of the damping cdω.Figure 7
Additional damping ratios provided by the damper: (a)nd=8; (b) nd=15.
(a)
(b)The absolute increment of the additional damping ratio reaches its maximum near the region of optimal damping values. To be different, the curves of the relative increment of the additional damping ratio and damping are plotted in Figure8. The vertical axis represents the relative ratio increased by the negative stiffness kd=−2.73N/mm to the case kd=0. The relative increment decreases with the increasing damping. For small values of damping, the negative stiffness remarkably increases the damping ratio. Also, a larger damper position enhances the damping ratio much more. For large values of damping, the negative stiffness has little effect on the improvement of the damping ratio wherever the damper is located.Figure 8
Relative increment of the additional damping ratio with negative stiffnesskd=−2.73N/mm.


## 3.1. Equivalent Parameters for the NSD

To suppress the cable vibration, Zhou and Li [22] developed a passive NSD, which is composed of two prepressed springs and an oil damper. The damping force consists of the negative stiffness and the Coulomb friction as(24)fNSD=kdx+Frsgnx˙,where fNSD is the damping force, kd=−2.73N/mm is the negative stiffness approximately determined by the NSD’s geometrical configuration, Fr=17N is the Coulomb friction level, and x and x˙ are the displacement and velocity of the damper, respectively.Making the dissipation energy per cycle the same, an equivalent viscous damping is adopted to replace the Coulomb friction part of the NSD when the cable-damper system vibrates in thejth mode as follows:(25)ceq=4FrπωjAd,in which ceq is the viscous damping coefficient and Ad is the amplitude of displacement at the damper position.Figure3 shows the equivalent hysteretic loops for one experimental case. The areas surrounded by the two curves are the same. Equivalent viscous damping coefficient ceq is inversely proportional to the product of the jth modal angular frequency ωj and the damper’s displacement amplitude Ad when the Coulomb friction force is kept unchanged. In the following verification, analytical results using the equivalent viscous damping coefficients will be compared with the experimental results.Figure 3
Equivalent hysteretic loops for the experimental NSD.


## 3.2. Discretization of the Continuous Cable

The length of the scaled stay cable used in the experiment [22] was 14m with the mass per unit length of 1.33kg/m. The tension force was adjusted to 10kN. Two damper locations were considered in the tests, which are 0.535m and 1.06m to the left anchorage, respectively. An electric shaker was employed to produce sinusoidal excitation forces acting at 0.345m from the right anchorage. For each damper location, the first and the second modal vibration was excited for the cable. After reaching the steady state, the electric shaker was cut off, and then the cable began its free vibration with exponential decaying amplitudes, which are explored to identify the additional modal damping ratios provided by the passive NSD.For the discrete model, the parameters, such as the node massm and the length l, are equal to the total mass and total length divided by the node number n. Then, the natural frequencies are determined by equation (7). Figure 4(a) depicts the frequency ratio of the discrete model to the continuous cable for the first three modes. With the node number n increase, the ratio converges to unity quickly and the relative error decreases dramatically. When the node number n is larger than 100, the relative error becomes less than 0.1%. Hence, n=200 is utilized in the discrete model to guarantee the precision. As shown in Figure 4(b), the first three modal frequencies of the discrete model are almost the same as the experimental results.Figure 4
(a) Frequency ratios for the first three modes; (b) comparison of the natural frequencies withn=200 to the experimental results.
(a)
(b)


## 3.3. Verification for the Displacement Response and Additional Modal Damping Ratio

The excitation force in the experiment is used to calculate the theoretical responses of the discrete model. The excitation position isne=195. Figure 5 compares the first modal displacement responses of the discrete model with experimental results when the damper is located at nd=8. The calculated displacement response at the damper location coincides well with the experimental result. Moreover, the calculated displacement response at the midspan of the cable agrees well with the experimental result.Figure 5
Displacements of the cable under the first modal vibration with the NSD: (a) damper location; (b) midspan.
(a)
(b)Specially, the relationship of displacement amplitude to damping at the damper’s location is shown in Figure6. The vertical axis is the dimensionless displacement amplitude according to equation (17). Under different levels of damping, experimental results are very close to the theoretical curve. The value of dimensionless displacement decreases with increasing damping cdω, as high damping will clamp the cable to make the damper’s location response small.Figure 6
Relationship of displacement amplitude to damping at the damper’s location.The additional modal damping ratios provided by the NSD are demonstrated in Figure7. In general, the experimental damping ratios at two locations agree with the theoretical curves according to the above derivation, although some points deviate a little largely. Comparing with damper stiffness kd=0, more damping ratios can be achieved when kd=−2.73N/mm, which means the negative stiffness is beneficial to enhance damping ratios. However, the beneficial effect varies with the level of the damping cdω.Figure 7
Additional damping ratios provided by the damper: (a)nd=8; (b) nd=15.
(a)
(b)The absolute increment of the additional damping ratio reaches its maximum near the region of optimal damping values. To be different, the curves of the relative increment of the additional damping ratio and damping are plotted in Figure8. The vertical axis represents the relative ratio increased by the negative stiffness kd=−2.73N/mm to the case kd=0. The relative increment decreases with the increasing damping. For small values of damping, the negative stiffness remarkably increases the damping ratio. Also, a larger damper position enhances the damping ratio much more. For large values of damping, the negative stiffness has little effect on the improvement of the damping ratio wherever the damper is located.Figure 8
Relative increment of the additional damping ratio with negative stiffnesskd=−2.73N/mm.


## 4. Numerical Analysis

The discrete model of the stay cable with an NSD has been deduced theoretically and validated experimentally in the former sections. In this section, numerical analysis will be further performed to reveal the optimal match of the negative stiffness and viscous damping based on the aforementioned derivations of the discrete model. Moreover, the influence and mechanism of the negative stiffness and viscous damping will be demonstrated in the perspective of energy.


### 4.1. Optimal Negative Stiffness and Viscous Damping

The effect of the excitation position on the optimal negative stiffness for the cable is analyzed firstly. As shown in Figure9, the optimal negative stiffness for the first and second modes is kept unchanged when the excitation position ne/n varies from 0.8 to about 1. So, the excitation position has no effect on the optimal negative stiffness.Figure 9
Change trend of the optimal negative stiffness to the excitation position.Being different from the excitation position, the damper location mainly determines the value of the optimal negative stiffness. Figure10 plots the curves of the optimal negative stiffness, which is inversely proportional to the damper location, agreeing with the following:(26)kdopt=−Tndl.Figure 10
Change trend of the optimal negative stiffness to the damper location.The relationship between the optimal damping and the negative stiffness is further described in Figure11. When the negative stiffness reaches its optimum, the corresponding optimal damping is zero for mode 1 and mode 2. For other values of negative stiffness, the optimal damping cdopt differs with the mode number. And less damping is needed for the higher mode.Figure 11
Optimal damping corresponding to the negative stiffness.The phase difference, as depicted in equation (18), caused by the optimal negative stiffness and damping at the damper location is depicted in Figure 12. The phase difference is zero when the optimal damping is zero, which means that the response of the damper location and the excitation are kept unchanged simultaneously. In contrast, when the damping tends to infinity, this phase difference approaches −π/4, regardless of the damper location.Figure 12
Phase difference at the damper location.


### 4.2. Energy Dissipated by the NSD

From the initial state, the excitation force maintains inputting energy until the cable-NSD system reaches its steady-state vibration. Under the steady state, the cable’s mechanical energy stays constant, while the excitation energy is always dissipated by the damper. Obviously, the cable mechanical energy or the energy dissipated by the NSD depends on the parameters of the damping and stiffness. The change tendencies of the first modal vibration are shown in Figure13, in which the horizontal axis is the damping and the vertical axis is the cable mechanical energy or the energy dissipated by the NSD itself. Four different levels of the damper stiffness, kd=2kdopt,kdopt,0, and kdopt, are employed to compare with each other. When the damper stiffness is negative, such as kd=2kdopt or kdopt, the cable mechanical energy increases with the damping. The energy dissipated by the NSD rises until reaching a peak value and then decreases with the damping. Especially, for the optimal negative stiffness and zero damping, the cable mechanical energy becomes zero in Figure 13(b). In contrast, the cable mechanical energy tends to infinity when the damping is zero and then decreases with the damping, as shown in both Figures 13(c) and 13(d), under the nonnegative conditions. When the damping tends to zero, the dissipated energy tends to infinity for kd=0 while tends to zero for kd=kdopt. For all the values of stiffness, the cable mechanical energy moves toward a constant value at last, while the energy dissipated reduces to zero as the damping increases to infinity.Figure 13
Change trend of the cable mechanical energy or the energy dissipated by the NSD to the damping with different levels of stiffness: (a)kd=2kdopt; (b) kd=kdopt; (c) kd=0; (d) kd=kdopt.
(a)
(b)
(c)
(d)Figure14 shows the change trend of the cable mechanical energy or dissipated energy with the normalized stiffness under four levels of viscous damping. It can be seen clearly that the minimum cable mechanical energy is zero when the normalized stiffness is negative unity and the damping is zero in Figure 14(a). In general, the cable mechanical energy decreases with stiffness when kd/kdopt≤−1 is satisfied. For kd/kdopt>−1, the cable mechanical energy firstly increases to the peak value and then declines with the increasing stiffness. The viscous damping can remarkably reduce the maximum cable energy, as shown in Figures 14(a)–14(d). With the same damping, the cable energy is generally much less when the stiffness is negative. The energy dissipated by the NSD achieves the maximum value around the kd/kdopt=0 region.Figure 14
Change trend of the cable mechanical energy or the energy dissipated by the NSD to the stiffness with different levels of damping: (a)cdω=0; (b) cdω=10N/mm; (c) cdω=20N/mm; (d) cdω=30N/mm.
(a)
(b)
(c)
(d)The results of other modes are similar to those of the first mode, so the curves of other modes are not drawn here. Under the steady-state forced vibration, the cable mechanical energy means the vibration intensity of the cable. Only the energy dissipated by the NSD is not enough to evaluate the control performance. The cable mechanical energy has a minimum value of zero under the condition of optimal negative stiffness and zero damping. This zero mechanical energy means that the cable has the weakest vibration. In other words, the parameters for best control performance are the combination of the optimal negative stiffness and zero damping.Figure15 plots the additional modal damping ratio for the stay cable with various levels of stiffness, which could be achieved by altering the stiffness coefficient of the springs and their compression degree in the proposed damper in [22], while the viscous damping depends on the oil damper. So, the proposed apparatus could be effectively employed for a real cable by choosing proper parameters. For kd/kdopt=−1, the additional modal damping ratio reaches infinity under zero damping. When the damping cdω≥100N/mm, the additional modal damping ratio is basically unchanged regardless of the level of stiffness.Figure 15
Additional modal damping ratio for the stay cable with various levels of stiffness.Furthermore, Figure16 compares the additional modal damping ratio with various levels of damping. For the same level of damping, the additional modal damping ratio achieves its peak value at the point kd/kdopt=−1. That is to say, the optimal negative stiffness enhances the effect of damping remarkably.Figure 16
Additional modal damping ratio for the stay cable with various levels of damping.


### 4.3. Energy Distribution along the Cable

This section will analyze the energy distribution along the cable when the NSD is attached for vibration control.The energy ratio of each node in the discrete model is demonstrated in Figure17. The horizontal axis is the dimensionless position, in which zero represents the left support and unity corresponds to the right support. The vertical axis is the energy ratio of each node to the total mechanical energy of the cable. The negative stiffness and viscous damping significantly affect the energy distribution along the cable and increase the energy at the damper location. In particular, the energy ratio at the damper location is only 0.3% when kd/kdopt=−1.2 and cd=1N⋅s/m. For the case kd/kdopt=−1, this ratio increases to 1%, while the energy ratio at the cable midspan declines to zero, which is obviously different from the original first mode. When the damping cd=108N⋅s/m, the energy distribution is similar to that of the first mode and the peak value of the energy ratio appears at the midspan of the cable.Figure 17
Distribution of the cable mechanical energy under the first modal vibration: (a)kd/kdopt=−1.2; (b) kd/kdopt=−1; (c) kd/kdopt=−0.8; (d) kd/kdopt=0.
(a)
(b)
(c)
(d)Under the second modal vibration, the energy ratio at the damper location is increased by the negative stiffness, as shown in Figure18. In particular, the energy ratio at the damper location reaches 1%, which is its peak value, when kd/kdopt=−1 and the damping cd approaches zero. For the high damping, energy distribution along the cable is unchanged and even the negative stiffness is optimal.Figure 18
Distribution of the cable mechanical energy under the second modal vibration: (a)kd/kdopt=−1.2; (b) kd/kdopt=−1; (c) kd/kdopt=−0.8; (d) kd/kdopt=0.
(a)
(b)
(c)
(d)


## 4.1. Optimal Negative Stiffness and Viscous Damping

The effect of the excitation position on the optimal negative stiffness for the cable is analyzed firstly. As shown in Figure9, the optimal negative stiffness for the first and second modes is kept unchanged when the excitation position ne/n varies from 0.8 to about 1. So, the excitation position has no effect on the optimal negative stiffness.Figure 9
Change trend of the optimal negative stiffness to the excitation position.Being different from the excitation position, the damper location mainly determines the value of the optimal negative stiffness. Figure10 plots the curves of the optimal negative stiffness, which is inversely proportional to the damper location, agreeing with the following:(26)kdopt=−Tndl.Figure 10
Change trend of the optimal negative stiffness to the damper location.The relationship between the optimal damping and the negative stiffness is further described in Figure11. When the negative stiffness reaches its optimum, the corresponding optimal damping is zero for mode 1 and mode 2. For other values of negative stiffness, the optimal damping cdopt differs with the mode number. And less damping is needed for the higher mode.Figure 11
Optimal damping corresponding to the negative stiffness.The phase difference, as depicted in equation (18), caused by the optimal negative stiffness and damping at the damper location is depicted in Figure 12. The phase difference is zero when the optimal damping is zero, which means that the response of the damper location and the excitation are kept unchanged simultaneously. In contrast, when the damping tends to infinity, this phase difference approaches −π/4, regardless of the damper location.Figure 12
Phase difference at the damper location.


## 4.2. Energy Dissipated by the NSD

From the initial state, the excitation force maintains inputting energy until the cable-NSD system reaches its steady-state vibration. Under the steady state, the cable’s mechanical energy stays constant, while the excitation energy is always dissipated by the damper. Obviously, the cable mechanical energy or the energy dissipated by the NSD depends on the parameters of the damping and stiffness. The change tendencies of the first modal vibration are shown in Figure13, in which the horizontal axis is the damping and the vertical axis is the cable mechanical energy or the energy dissipated by the NSD itself. Four different levels of the damper stiffness, kd=2kdopt,kdopt,0, and kdopt, are employed to compare with each other. When the damper stiffness is negative, such as kd=2kdopt or kdopt, the cable mechanical energy increases with the damping. The energy dissipated by the NSD rises until reaching a peak value and then decreases with the damping. Especially, for the optimal negative stiffness and zero damping, the cable mechanical energy becomes zero in Figure 13(b). In contrast, the cable mechanical energy tends to infinity when the damping is zero and then decreases with the damping, as shown in both Figures 13(c) and 13(d), under the nonnegative conditions. When the damping tends to zero, the dissipated energy tends to infinity for kd=0 while tends to zero for kd=kdopt. For all the values of stiffness, the cable mechanical energy moves toward a constant value at last, while the energy dissipated reduces to zero as the damping increases to infinity.Figure 13
Change trend of the cable mechanical energy or the energy dissipated by the NSD to the damping with different levels of stiffness: (a)kd=2kdopt; (b) kd=kdopt; (c) kd=0; (d) kd=kdopt.
(a)
(b)
(c)
(d)Figure14 shows the change trend of the cable mechanical energy or dissipated energy with the normalized stiffness under four levels of viscous damping. It can be seen clearly that the minimum cable mechanical energy is zero when the normalized stiffness is negative unity and the damping is zero in Figure 14(a). In general, the cable mechanical energy decreases with stiffness when kd/kdopt≤−1 is satisfied. For kd/kdopt>−1, the cable mechanical energy firstly increases to the peak value and then declines with the increasing stiffness. The viscous damping can remarkably reduce the maximum cable energy, as shown in Figures 14(a)–14(d). With the same damping, the cable energy is generally much less when the stiffness is negative. The energy dissipated by the NSD achieves the maximum value around the kd/kdopt=0 region.Figure 14
Change trend of the cable mechanical energy or the energy dissipated by the NSD to the stiffness with different levels of damping: (a)cdω=0; (b) cdω=10N/mm; (c) cdω=20N/mm; (d) cdω=30N/mm.
(a)
(b)
(c)
(d)The results of other modes are similar to those of the first mode, so the curves of other modes are not drawn here. Under the steady-state forced vibration, the cable mechanical energy means the vibration intensity of the cable. Only the energy dissipated by the NSD is not enough to evaluate the control performance. The cable mechanical energy has a minimum value of zero under the condition of optimal negative stiffness and zero damping. This zero mechanical energy means that the cable has the weakest vibration. In other words, the parameters for best control performance are the combination of the optimal negative stiffness and zero damping.Figure15 plots the additional modal damping ratio for the stay cable with various levels of stiffness, which could be achieved by altering the stiffness coefficient of the springs and their compression degree in the proposed damper in [22], while the viscous damping depends on the oil damper. So, the proposed apparatus could be effectively employed for a real cable by choosing proper parameters. For kd/kdopt=−1, the additional modal damping ratio reaches infinity under zero damping. When the damping cdω≥100N/mm, the additional modal damping ratio is basically unchanged regardless of the level of stiffness.Figure 15
Additional modal damping ratio for the stay cable with various levels of stiffness.Furthermore, Figure16 compares the additional modal damping ratio with various levels of damping. For the same level of damping, the additional modal damping ratio achieves its peak value at the point kd/kdopt=−1. That is to say, the optimal negative stiffness enhances the effect of damping remarkably.Figure 16
Additional modal damping ratio for the stay cable with various levels of damping.


## 4.3. Energy Distribution along the Cable

This section will analyze the energy distribution along the cable when the NSD is attached for vibration control.The energy ratio of each node in the discrete model is demonstrated in Figure17. The horizontal axis is the dimensionless position, in which zero represents the left support and unity corresponds to the right support. The vertical axis is the energy ratio of each node to the total mechanical energy of the cable. The negative stiffness and viscous damping significantly affect the energy distribution along the cable and increase the energy at the damper location. In particular, the energy ratio at the damper location is only 0.3% when kd/kdopt=−1.2 and cd=1N⋅s/m. For the case kd/kdopt=−1, this ratio increases to 1%, while the energy ratio at the cable midspan declines to zero, which is obviously different from the original first mode. When the damping cd=108N⋅s/m, the energy distribution is similar to that of the first mode and the peak value of the energy ratio appears at the midspan of the cable.Figure 17
Distribution of the cable mechanical energy under the first modal vibration: (a)kd/kdopt=−1.2; (b) kd/kdopt=−1; (c) kd/kdopt=−0.8; (d) kd/kdopt=0.
(a)
(b)
(c)
(d)Under the second modal vibration, the energy ratio at the damper location is increased by the negative stiffness, as shown in Figure18. In particular, the energy ratio at the damper location reaches 1%, which is its peak value, when kd/kdopt=−1 and the damping cd approaches zero. For the high damping, energy distribution along the cable is unchanged and even the negative stiffness is optimal.Figure 18
Distribution of the cable mechanical energy under the second modal vibration: (a)kd/kdopt=−1.2; (b) kd/kdopt=−1; (c) kd/kdopt=−0.8; (d) kd/kdopt=0.
(a)
(b)
(c)
(d)


## 5. Conclusions

The match of negative stiffness and viscous damping for cable vibration mitigation is investigated in this paper. A discrete model for the stay cable with an NSD was firstly established, and the system dynamic characteristics, including the forced steady-state vibration and additional modal damping ratio, were performed theoretically. After that, experimental results were utilized to verify the discrete model and the corresponding derivations. Then, numerical analysis was carried out further to reveal the effect of negative stiffness and viscous damping in the cable vibration mitigation. Moreover, the energy dissipated by the NSD and the energy distribution along the cable were evaluated. Based on the results of this study, the following conclusions could be drawn.A discrete model, which can effectively describe the system dynamic behavior, was established for the stay cable with an NSD attached near the support. The forced steady-state vibration response and the additional modal damping ratio based on this model agree with the experimental results well.The cable energy and the energy dissipated by the NSD depend on the negative stiffness and viscous damping. With the optimal negative stiffness and zero damping, the cable mechanical energy declines to the weakest vibration condition.In particular, the match of negative stiffness and viscous damping also has a significant impact on the energy distribution along the cable. More precisely, the negative stiffness concentrates the energy to the damper location. For the optimal negative stiffness and zero damping, the energy ratio at the damper location approaches its peak value.



---

*Source: 3208321-2019-08-04.xml*