# Development of a Negative Stiffness Bistable Damper for Structural Vibration Control

Authors: Liming Fan; Chen Huang; Linsheng Huo
Journal: Shock and Vibration
 (2022)

---

## Abstract

Linear dampers have been widely applied for suppressing the dynamic responses of structures to mitigate their damage. However, the primary disadvantage of the classical linear damper is that it is vulnerable to detuning, which has become an issue of great importance recently due to a great reduction in vibration control performance. To overcome the shortcoming, this study develops a negative stiffness bistable damper (NSBD) composed of a simple assembly consisting of a bistable buckling beam with a mass. Energy is dissipated through the transformation between the bistable states. The constitutive equation of the NSBD is derived to analyze the effects of the stiffness ratio, the arch-span ratio, and the damping ratio on its restoring capabilities. The vibration reduction effect of the NSBD is experimentally evaluated under different sinusoidal and seismic excitations in shaking table tests. The obtained results reveal that the NSBD can effectively restrain structural displacements.


---

## Body

## 1. Introduction

Vibration control is well-established since first proposed by Yao [1] in 1972, yet still developed to enhance the functionality and safety of structures. Up to now, many structures have been installed with various types of vibration control devices to restrain their responses to wind loads, seismic action, or other dynamic excitations [2–5]. Linear dampers (e.g., tuned mass damper) have been widely applied in engineering practice, such as the vibration control systems of the Akashi Kaikyo Bridge [6], the Taipei 101 tower [7], and the Shanghai Center Tower [8].Generally, the dynamic response of a structure is greatly influenced by the first frequency and vibration mode [9–11]. The conventional tuned mass damper (TMD) must keep consistent with the first frequency of the primary structure to effectively restrain the structural vibration. Therefore, as a type of frequency-sensitive device, the conventional TMD operates within a relatively narrow frequency band for vibration reduction. The TMD is detuned and its effectiveness can be evidently reduced when the natural frequency of the host structure alters.To overcome this deficiency, many negative stiffness vibration control devices have been developed. The concept of negative stiffness was first introduced by Molyneaux [12]. When the directions of restoring force and corresponding deformation are opposite, their ratio is negative. A structure in such a state has negative stiffness. For structural vibration control, the negative stiffness can reduce the stiffness of the control system and thus lower its natural frequency, allowing it to provide vibration reduction over a wide frequency range. Negative stiffness can be achieved via various methods. Traditionally, negative stiffness is achieved by special mechanical devices for energy storage elements (e.g., spring and prebuckled beam) [13–17], magnetic elements [18–20], geometric nonlinearities [21–23], or composite structures and metamaterials [23, 24]. Virgin et al. [25] proposed negative vibration isolation using a highly deformed, slender beam. The isolator had a high adjustable stiffness and a high rate of displacement transmission over a wide range. Shi and Zhu [26] developed and experimentally verified two configurations of magnetic negative stiffness dampers (MNSDs), in which multiple permanent magnets are placed in a conductive tube. Zhou et al. [27]. introduced a negative stiffness vibration isolation damper (NSVID) using a nonlinear structure and utilized the NSVID to control torsional vibration along a shaft. This proposed device can realize quasi-zero stiffness (QZS) by changing the stiffness of rubber and control parameters. Xiang et al. [28] designed a semiactive control torsional vibration damper, which is composed of positive and negative stiffness in parallel. They provided a derivation of the theory behind the damper mechanism and analyzed its elastic properties. Zhou et al. [29] conducted theoretical analysis and experimental studies on a passive negative stiffness damper (PNSD) in series with flexible support. The results indicated that their proposed device can achieve a range of equivalent negative stiffnesses and Coulomb values. Haghpanah et al. [24] invented a negative stiffness element composed of one convex elastic element together with elastic elements that have nonconvex strain energy. This element can greatly enhance damping performance under dynamic excitation. Recently, numerous other negative stiffness dampers have been introduced for vibration control of various types of excitations [30–34] and different types of structures [35–37].Bistable structures are structures that not only have negative stiffness properties, but can also provide cubic stiffness to enhance vibration control. Most previous studies focused on mechanical properties analysis [38, 39] and vibration isolation applications of bistable structures. Johnson et al. [40] studied the vibration-attenuating effect of a bistable oscillator and investigated its snap-through dynamics using the harmonic balance method and experimentation. Experimental results show that the bistable oscillator has a slightly superior vibration attenuation performance compared to the fundamental harmonic snap-through action. Farhangdoust et al. [41]. introduced a bistable tuned mass damper (BTMD) mechanism for suppressing the vortex-induced vibration (VIV) in suspension bridge decks, which are exceptionally sensitive to broadband input of vortex shedding velocity. Xia et al. [42] investigated the enhancement in vibration attenuation bandwidth of a locally resonant metamaterial beam with bistable attachments, the interwell oscillations of the attachments produced an attenuation frequency range of 350% wider than linear locally resonant bandgap for moderate-to-high amplitude excitation levels. Zhao et al. [43] explored multiple snap-through pathways and force-displacement curves of a discontinuous point and verified the results through theoretical analysis and numerical simulation. Liu et al. [44] used equivalent, analytical, numerical, and experimental methods to investigate 1/2 sub-harmonic resonance in a bistable structure and its vibration isolation characteristics. Zhang et al. [45] designed a recoverable inclined beam energy absorption structure that leverages bistable characteristics to improve the safety and efficiency of passive energy absorption systems in automobiles. Yan et al. [46, 47] invented a bistable vibration isolation (BVI) device composed of several ring permanent magnets (PMS). The PMS improved the vibration isolation performance of the BVI through nonlinear electromagnetic shunt damping. Zhang et al. [48] utilized an optimized and varying sectional profile to improve the performance of a quasi-zero stiffness (QZS) isolation system using the nonlinear and negative stiffness generated by the snap-through effect of bistable structures. Yang et al. [49] proposed a feedback control law to modulate the current input to the actuator of an active vibration isolation system. The proposed control law greatly attenuated the transmissibility of the bistable nonlinear electromagnetic actuator from the base excitation.Inspired by the advantages of negative stiffness and bistable structures, this paper proposes a negative stiffness bistable damper (NSBD) composed of a bistable buckling beam assembled with a mass. The NSBD can enhance energy dissipation through the transformation between the bistable states. The key component of the damper is a precompressed spring steel set in a clamped-clamped configuration in the frame structure. The constitutive equation of the NSBD is derived to analyze the effects of the stiffness ratio, the arch-span ratio, and the damping ratio on the restoring force of the NSBD. To validate the vibration reduction effect of the NSBD, a frame structure was subjected to shake table tests with and without the NSBD. Finally, a following discussion of the experimental results is presented to explore the damping effect of different parameters on the NSBD.


## 2. Constitutive Equation of Negative Stiffness Bistable Damper

### 2.1. Mechanical Characterization of Negative Stiffness Bistable Damper

A schematic of a linear single-degree-of-freedom (SDOF) structure equipped with an NSBD is depicted in Figure1. The mass, stiffness, and damping coefficient of the structure are m1, k1, and c1, respectively. In the schematic, kns is the stiffness of the NSBD and consists of the negative stiffness k21 and the cubic stiffness k22. The key component of the NSBD is the bistable buckling beam, which is clamped within the frame structure and precompressed. An attached mass is fixed at the midspan of the bistable buckling beam. When an external excitation acts on the primary structure, the NSBD will be continuously transformed between two equilibrium positions.Figure 1
Schematic diagram of NSBD.Under small transverse displacement, the prestress imposed along the axial direction (y direction) of the straight beam to move the right end of the beam from position b to position a results in displacement wy in the y direction. The prestress also causes the force P to be applied in the x-direction at the midpoint of the beam, resulting in displacement wx in the x-direction, as shown in Figure 2. The deflection curve of the deformed beam is the first-order mode shape of the clamped-clamped beam.(1)wy=wx21−cos2πylc.Figure 2
The buckled beam with prestress.Assuming that the beam section thickness isδ and the beam section width is b, the moment of inertia and cross-sectional area of the beam section are, respectively, I=bδ3/12 and A=bδ. The axial stiffness of the beam is k0=AE/lc. According to Jin et al. [50], the relationship between the external force on the buckling beam and the displacement is as follows:(2)F0=3π4Q22ΔΔ−32−14−43Q2×Δ−32+14−43Q2,where F0=plc3/EIδ, Δ=σ/δ, and Q=wy/δ are normalized parameters. lc is the span of the buckling beam and δ is the thickness. σ is the intermediate deformation of the buckling beam under external force, and wy is the arching height of the buckling beam.The model parameters of the bistable buckling beam are solved by equation (2) and transformed into the ordinary form:(3)Plc3EIwy=3π42wyδ2σwy3−9π42wyδ2σwy2+3π4wyδ2+2π4σwy,(4)P=3π4EI2lc3δ2σ3−9π4EIwy2lc3δ2σ2+3π4wy2EIδ2lc3+2π4EIlc3σ.The formula (5) can be obtained by placing the zero position at the origin.(5)P=3π4EI2lc3δ2x+σ3−9π4EIwy2lc3δ2x+σ2+3π4wy2EIδ2lc3+2π4EIlc3x+σ=3π4EI2lc3δ2x3−3π4wy2EI2δ2lc3−2π4EIlc3x+2π4EIwylc3.Thus, the cubic stiffness coefficient of the buckling beam can be expressed as3π4EI/2lc3δ2, The negative stiffness coefficient of the buckling beam is −3π4EI/2δ2lc3−2π4EI/lc3. When E = 198.6 GPa, lc = 0.36 m, δ = 0.0002 m, wy = 0.025 m, and b = 0.02 m, the relationship between force and displacement of the buckling beam is shown in Figure 3.Figure 3
Force and displacement curve of buckling beam.


### 2.2. Stiffness Analysis

From the above analysis, it can be seen that the stiffness expression of the buckling beam is as follows:(6)k=∂P∂x=3π4EI2lc3δ2x2−3π4wy2EI2δ2lc3−2π4EIlc3=3π4EI2lc3x2δ2−wy2δ2+43.The axial stiffness of the buckling beam is(7)k0=AElc.Then, the equivalent stiffness of the buckling beam is(8)k¯=kk0=π48x¯2−ϕ2+43η2,where k¯=k/k0, η=δ/lc, ϕ=wy/lc, x¯=x/lc.The variation of equivalent stiffness of buckling beam with different camber heights is shown in Figure4.Figure 4
Relationship between equivalent stiffness of buckling beam and camber height.


## 2.1. Mechanical Characterization of Negative Stiffness Bistable Damper

A schematic of a linear single-degree-of-freedom (SDOF) structure equipped with an NSBD is depicted in Figure1. The mass, stiffness, and damping coefficient of the structure are m1, k1, and c1, respectively. In the schematic, kns is the stiffness of the NSBD and consists of the negative stiffness k21 and the cubic stiffness k22. The key component of the NSBD is the bistable buckling beam, which is clamped within the frame structure and precompressed. An attached mass is fixed at the midspan of the bistable buckling beam. When an external excitation acts on the primary structure, the NSBD will be continuously transformed between two equilibrium positions.Figure 1
Schematic diagram of NSBD.Under small transverse displacement, the prestress imposed along the axial direction (y direction) of the straight beam to move the right end of the beam from position b to position a results in displacement wy in the y direction. The prestress also causes the force P to be applied in the x-direction at the midpoint of the beam, resulting in displacement wx in the x-direction, as shown in Figure 2. The deflection curve of the deformed beam is the first-order mode shape of the clamped-clamped beam.(1)wy=wx21−cos2πylc.Figure 2
The buckled beam with prestress.Assuming that the beam section thickness isδ and the beam section width is b, the moment of inertia and cross-sectional area of the beam section are, respectively, I=bδ3/12 and A=bδ. The axial stiffness of the beam is k0=AE/lc. According to Jin et al. [50], the relationship between the external force on the buckling beam and the displacement is as follows:(2)F0=3π4Q22ΔΔ−32−14−43Q2×Δ−32+14−43Q2,where F0=plc3/EIδ, Δ=σ/δ, and Q=wy/δ are normalized parameters. lc is the span of the buckling beam and δ is the thickness. σ is the intermediate deformation of the buckling beam under external force, and wy is the arching height of the buckling beam.The model parameters of the bistable buckling beam are solved by equation (2) and transformed into the ordinary form:(3)Plc3EIwy=3π42wyδ2σwy3−9π42wyδ2σwy2+3π4wyδ2+2π4σwy,(4)P=3π4EI2lc3δ2σ3−9π4EIwy2lc3δ2σ2+3π4wy2EIδ2lc3+2π4EIlc3σ.The formula (5) can be obtained by placing the zero position at the origin.(5)P=3π4EI2lc3δ2x+σ3−9π4EIwy2lc3δ2x+σ2+3π4wy2EIδ2lc3+2π4EIlc3x+σ=3π4EI2lc3δ2x3−3π4wy2EI2δ2lc3−2π4EIlc3x+2π4EIwylc3.Thus, the cubic stiffness coefficient of the buckling beam can be expressed as3π4EI/2lc3δ2, The negative stiffness coefficient of the buckling beam is −3π4EI/2δ2lc3−2π4EI/lc3. When E = 198.6 GPa, lc = 0.36 m, δ = 0.0002 m, wy = 0.025 m, and b = 0.02 m, the relationship between force and displacement of the buckling beam is shown in Figure 3.Figure 3
Force and displacement curve of buckling beam.


## 2.2. Stiffness Analysis

From the above analysis, it can be seen that the stiffness expression of the buckling beam is as follows:(6)k=∂P∂x=3π4EI2lc3δ2x2−3π4wy2EI2δ2lc3−2π4EIlc3=3π4EI2lc3x2δ2−wy2δ2+43.The axial stiffness of the buckling beam is(7)k0=AElc.Then, the equivalent stiffness of the buckling beam is(8)k¯=kk0=π48x¯2−ϕ2+43η2,where k¯=k/k0, η=δ/lc, ϕ=wy/lc, x¯=x/lc.The variation of equivalent stiffness of buckling beam with different camber heights is shown in Figure4.Figure 4
Relationship between equivalent stiffness of buckling beam and camber height.


## 3. Parametric Analysis

### 3.1. Mechanical Constitutive Equation of Buckling Beam

From the above derivation, it can be seen that the mechanical constitutive equation of the buckling beam can be composed of a negative stiffness term and a cubic nonlinear term as follows:(9)F=−k21x+k22x3,where k21 and k22 are, respectively, the negative stiffness coefficient and cubic stiffness coefficient. Therefore, the potential energy function of the buckling beam is as follows:(10)Ux=−12k21x2+14k22x4.LettingF=0, the equilibrium position of the buckling beam can be obtained as shown in equation (11).(11)x0=±k21k22.Although the stiffness of the buckling beam changes with a change in position during vibration, the response characteristics are approximately linear when the buckling beam vibrates at the equilibrium point with a very small amplitude. Therefore, the stiffness of the buckling beam is simplified with the linear stiffness at the equilibrium position in the parametric analysis. In the analysis, the equivalent stiffness is used for simplification. The equivalent stiffness is as follows:(12)ke=d2Uxdx2=−k21+3k22x2|x=x0=2k21.The equivalent frequency of the buckled beam is(13)ωe=kem.From the potential energy function of the buckling beam, the bistable limit position can be solved as(14a)Ux=−12k21x2+14k22x4=0,(14b)xcr=±2k21k22=x02.It can be seen that when the applied external load meets the conditionA0>xcr, the buckling beam is in bistable vibration, and when A0≤xcr, the buckling beam experiencing monostable vibration.


### 3.2. Analytical Derivation

When the buckling beam with fixed ends is subjected to sinusoidal excitation with amplitudeF and frequency ω, the beam can be simplified into an SDOF system. The dynamic equation of the buckled beam is as follows:(15)mx¨+cx˙−k21x+k22x3=Fsinωt,where m, c,and k21, k22 refer to the mass, damping, and negative stiffness and cubic stiffness of the buckling beam, respectively; the equation (16) can be simplified to(16)y¨+2ξy˙−y+αf2y3=sinθτ,where ξ=c/2mω, α=k22/k21, ω1=k21/m, θ=ω/ω1, y=x/f, f=F/k22, and τ=θt.In order to obtain an approximate solution, the harmonic balance method is adopted in the paper, as used to solve equation (16); the first-order approximate solution is assumed to be as follows:(17)y=Ysinθτ+φ.In the analysis, the first-order approximate solution is considered, by substituting equation (17) into equation (16), and approximating sin3θ=3/4sinθ−1/4sin3θ≈3/4sinθ, the following relationships are obtained:(18a)−θ2Y−Y+αf2×34Y3=cosφ,(18b)2ξθY=−sinφ.By combined equations (18a) and (18b) and eliminating φ, equation (19) can be obtained as follows:(19)916α2f4Y6−32αf2Y41+θ2+4ξ2θ2+θ4+2θ2+1Y2−1=0.Rearranging equation (19) to be expressed in terms of θ is(20)Y2θ4+4ξ2Y2−32αf2Y4+2Y2θ2+916α2f4Y6−32αf2Y4+Y2−1=0.The solution of the equation is thus(21)θ=34αf2Y2−2ξ2−1±4ξ4+4ξ2−3αξ2f2Y2+1Y2.The amplitude-frequency response function of the buckled beam can be obtained by solving equation (21), as shown in Figure 5.Figure 5
The amplitude-frequency response curve of the buckled beam.
(a)(b)(c)The amplitude-frequency response curve of a linear system is constant; however, the amplitude-frequency response curve of the nonlinear system depends on the initial value conditions. Figure5 reveals the relation between the amplitude of the response of the nonlinear system and the excitation frequency at the following initial conditions: (a) ξ=0.03, α=0.3, (b) ξ=0.03, f=0.5, and (c) f=0.5, α=0.3. It can be noted that the negative stiffness term and the nonlinear term make the curve bend towards the right, resulting in one frequency value corresponding to multiple amplitudes in the frequency range.Figure5(a) shows the relationship between the amplitude of the system (ξ=0.03, α=0.3) response at different f and the excitation frequency. As the nonlinear coefficient increases, the maximum amplitude of the response decreases, but the degree of curvature does not change. The upward jump point is also getting smaller.Figure5(b) shows the relationship between the amplitude response of the nonlinear system (ξ=0.03, f=0.5) at different ∂ and the excitation frequency. Note that as the nonlinear coefficient increases, the bend of the curve does not change appreciably, but the upward jump point also decreases.Figure5(c) shows the relationship between the amplitude of the system response (f=0.5, α=0.3) with different damping ratio ξ and the excitation frequency. With the increase of the damping ratio, the bending degree of the curve does not change, and the upward jump point is also unchanged.


## 3.1. Mechanical Constitutive Equation of Buckling Beam

From the above derivation, it can be seen that the mechanical constitutive equation of the buckling beam can be composed of a negative stiffness term and a cubic nonlinear term as follows:(9)F=−k21x+k22x3,where k21 and k22 are, respectively, the negative stiffness coefficient and cubic stiffness coefficient. Therefore, the potential energy function of the buckling beam is as follows:(10)Ux=−12k21x2+14k22x4.LettingF=0, the equilibrium position of the buckling beam can be obtained as shown in equation (11).(11)x0=±k21k22.Although the stiffness of the buckling beam changes with a change in position during vibration, the response characteristics are approximately linear when the buckling beam vibrates at the equilibrium point with a very small amplitude. Therefore, the stiffness of the buckling beam is simplified with the linear stiffness at the equilibrium position in the parametric analysis. In the analysis, the equivalent stiffness is used for simplification. The equivalent stiffness is as follows:(12)ke=d2Uxdx2=−k21+3k22x2|x=x0=2k21.The equivalent frequency of the buckled beam is(13)ωe=kem.From the potential energy function of the buckling beam, the bistable limit position can be solved as(14a)Ux=−12k21x2+14k22x4=0,(14b)xcr=±2k21k22=x02.It can be seen that when the applied external load meets the conditionA0>xcr, the buckling beam is in bistable vibration, and when A0≤xcr, the buckling beam experiencing monostable vibration.


## 3.2. Analytical Derivation

When the buckling beam with fixed ends is subjected to sinusoidal excitation with amplitudeF and frequency ω, the beam can be simplified into an SDOF system. The dynamic equation of the buckled beam is as follows:(15)mx¨+cx˙−k21x+k22x3=Fsinωt,where m, c,and k21, k22 refer to the mass, damping, and negative stiffness and cubic stiffness of the buckling beam, respectively; the equation (16) can be simplified to(16)y¨+2ξy˙−y+αf2y3=sinθτ,where ξ=c/2mω, α=k22/k21, ω1=k21/m, θ=ω/ω1, y=x/f, f=F/k22, and τ=θt.In order to obtain an approximate solution, the harmonic balance method is adopted in the paper, as used to solve equation (16); the first-order approximate solution is assumed to be as follows:(17)y=Ysinθτ+φ.In the analysis, the first-order approximate solution is considered, by substituting equation (17) into equation (16), and approximating sin3θ=3/4sinθ−1/4sin3θ≈3/4sinθ, the following relationships are obtained:(18a)−θ2Y−Y+αf2×34Y3=cosφ,(18b)2ξθY=−sinφ.By combined equations (18a) and (18b) and eliminating φ, equation (19) can be obtained as follows:(19)916α2f4Y6−32αf2Y41+θ2+4ξ2θ2+θ4+2θ2+1Y2−1=0.Rearranging equation (19) to be expressed in terms of θ is(20)Y2θ4+4ξ2Y2−32αf2Y4+2Y2θ2+916α2f4Y6−32αf2Y4+Y2−1=0.The solution of the equation is thus(21)θ=34αf2Y2−2ξ2−1±4ξ4+4ξ2−3αξ2f2Y2+1Y2.The amplitude-frequency response function of the buckled beam can be obtained by solving equation (21), as shown in Figure 5.Figure 5
The amplitude-frequency response curve of the buckled beam.
(a)(b)(c)The amplitude-frequency response curve of a linear system is constant; however, the amplitude-frequency response curve of the nonlinear system depends on the initial value conditions. Figure5 reveals the relation between the amplitude of the response of the nonlinear system and the excitation frequency at the following initial conditions: (a) ξ=0.03, α=0.3, (b) ξ=0.03, f=0.5, and (c) f=0.5, α=0.3. It can be noted that the negative stiffness term and the nonlinear term make the curve bend towards the right, resulting in one frequency value corresponding to multiple amplitudes in the frequency range.Figure5(a) shows the relationship between the amplitude of the system (ξ=0.03, α=0.3) response at different f and the excitation frequency. As the nonlinear coefficient increases, the maximum amplitude of the response decreases, but the degree of curvature does not change. The upward jump point is also getting smaller.Figure5(b) shows the relationship between the amplitude response of the nonlinear system (ξ=0.03, f=0.5) at different ∂ and the excitation frequency. Note that as the nonlinear coefficient increases, the bend of the curve does not change appreciably, but the upward jump point also decreases.Figure5(c) shows the relationship between the amplitude of the system response (f=0.5, α=0.3) with different damping ratio ξ and the excitation frequency. With the increase of the damping ratio, the bending degree of the curve does not change, and the upward jump point is also unchanged.


## 4. Experiment Setup and Conditions

### 4.1. Experiment Setup

Based on the above discussion, experiments are designed to verify the attenuating effect of NSBD on the vibration response of frame structure and study the influence of the relevant parameters on the vibration response of frame structure. The main component of the damper is a buckling beam made of spring steel, which can reach different camber heights depending on an external force. The clamped-clamped buckling beam is fixed within a steel frame supported by four circular steel columns. To better enable the bistable effect, the attached mass is added to the midspan of the buckling beam. The whole experimental device is fastened on the shaking table using bolts. Under the transverse acceleration provided by the shaking table, the mass block is used to provide the driving force according toF=ma to cause the device to switch between the two equilibrium positions and thus provide vibration reduction. The experimental device is shown in Figure 6. The experiments are conducted with and without the NSBD coupled with the frame structure.Figure 6
Experimental device.A data acquisition system is (DH5922D dynamic signal test and analysis system) connected to a computer recorded data from a laser displacement sensor that measures the movements of the frame. The setup of the experiment and instrumentation is shown in Figure6. As shown in Figure 6, before installing the NSBD system, holes are reserved in the center of the beam and the square frame, respectively. In the experiment, the buckling beam is guided freely through the clamping parts of the square frame at both ends. Then, prepressure is applied to both ends of the buckling beam to make the buckling beam reach the specified height, and the corresponding mass blocks are bolted in the middle of the buckling beam according to the experimental requirements after the two ends of the buckling beam are fixed; then the NSBD system is completed.The buckling beam in the experiment is arched using spring steel with a prepress at the end, and a mass composed of a square iron block is fixed in the middle of the span. The mass of the block can be changed to adjust the characteristics of the NSBD. In the experiment, the spring steel is restrained by a square steel tube to control the height of the buckling beam arch. The support system of the frame structure consisted of four round steel columns that are fixed on the shaking table with bolts. The experimental parameters are listed in Table1.Table 1
Experimental parameters.
Experimental materialsParameter valueABFrame size20 mm × 40 mm × 150 mm20 mm × 40 mm × 400 mmBuckling beam thickness0.2 mm0.25 mmMass of mass block69.26 g111.6 gBuckling beam width20 mmBuckling beam density7.85 g/cm3Buckling beam elastic modulus198.6 N/mm3Support systemØ4mm × 600 mm


### 4.2. Experimental Conditions

In the experiment, the effects of the arch-span ratio (i.e., the ratio of arch height to the buckling beam span) and the mass ratio (i.e., the ratio of additional mass block mass to the frame mass) on the vibration reduction performance of NSBD are considered. The arch-span ratio and mass used in the experiment are listed in Tables2 and 3.Table 2
Arch-span ratio.
Arch height (mm)Beam span (mm)Ratio (%)1536000.42250.69350.97Table 3
Mass ratio.
Attachment mass (g)Mass of primary structure (g)Ratio (%)69.2650501.37111.62.21The experiment validated the attenuating effect of NSBD on the frame structure under sinusoidal and earthquake excitation. Firstly, the shaking table is used to perform a frequency sweep analysis of the experimental frame structure. The analysis determined that the natural frequency of the experimental frame structure is 2.89 Hz. Based on the switch threshold force of the buckling beams for different thicknesses, the amplitude of sinusoidal excitation applied to the frame structure is 0.05 cm. Secondly, the experiments used sinusoidal excitations to examine the damping provided by the NSBD for frame structures of different thicknesses. Finally, a total of 10 different earthquake excitations selected from two types of sites are used to investigate the vibration control effectiveness of the NSBD at the best damping. According to the code (GB5011-2010) [51], 10 different earthquake excitations are selected, as listed in Table 4, which represent different site conditions.Table 4
Information of the selected seismic records.
Site classificationEarthquake nameYearStationMagnitudeRrup (km)Vs30 (m/s)IImperial Valley-031951EI Centro Array #95.625.24213.44Southern Calif1952San Luis Obispo673.41493.5Imperial Valley-051955El Centro Array #95.414.88213.44Central Calif-021960Hollister City Hall59.02198.77Hollister-011961Hollister City Hall5.619.56198.77San Fernando1971Pasadena - Old seismo lab6.6121.5969.07IIGuy_2010-10-15-HHE2010X102 Temporary Station_ Greenbriar_ AR3.8612.26495.3Guy_2010-10-15-HHN2010X102 Temporary Station_ greenbriar_ AR3.8612.26495.3Mineral_2011-08-232011Keystone College_ La Plume_ PA5.74442.73689.9Sparks_2011-11-052011Smith Ranch_ Marlow_ OK4.73148.55548.4


## 4.1. Experiment Setup

Based on the above discussion, experiments are designed to verify the attenuating effect of NSBD on the vibration response of frame structure and study the influence of the relevant parameters on the vibration response of frame structure. The main component of the damper is a buckling beam made of spring steel, which can reach different camber heights depending on an external force. The clamped-clamped buckling beam is fixed within a steel frame supported by four circular steel columns. To better enable the bistable effect, the attached mass is added to the midspan of the buckling beam. The whole experimental device is fastened on the shaking table using bolts. Under the transverse acceleration provided by the shaking table, the mass block is used to provide the driving force according toF=ma to cause the device to switch between the two equilibrium positions and thus provide vibration reduction. The experimental device is shown in Figure 6. The experiments are conducted with and without the NSBD coupled with the frame structure.Figure 6
Experimental device.A data acquisition system is (DH5922D dynamic signal test and analysis system) connected to a computer recorded data from a laser displacement sensor that measures the movements of the frame. The setup of the experiment and instrumentation is shown in Figure6. As shown in Figure 6, before installing the NSBD system, holes are reserved in the center of the beam and the square frame, respectively. In the experiment, the buckling beam is guided freely through the clamping parts of the square frame at both ends. Then, prepressure is applied to both ends of the buckling beam to make the buckling beam reach the specified height, and the corresponding mass blocks are bolted in the middle of the buckling beam according to the experimental requirements after the two ends of the buckling beam are fixed; then the NSBD system is completed.The buckling beam in the experiment is arched using spring steel with a prepress at the end, and a mass composed of a square iron block is fixed in the middle of the span. The mass of the block can be changed to adjust the characteristics of the NSBD. In the experiment, the spring steel is restrained by a square steel tube to control the height of the buckling beam arch. The support system of the frame structure consisted of four round steel columns that are fixed on the shaking table with bolts. The experimental parameters are listed in Table1.Table 1
Experimental parameters.
Experimental materialsParameter valueABFrame size20 mm × 40 mm × 150 mm20 mm × 40 mm × 400 mmBuckling beam thickness0.2 mm0.25 mmMass of mass block69.26 g111.6 gBuckling beam width20 mmBuckling beam density7.85 g/cm3Buckling beam elastic modulus198.6 N/mm3Support systemØ4mm × 600 mm


## 4.2. Experimental Conditions

In the experiment, the effects of the arch-span ratio (i.e., the ratio of arch height to the buckling beam span) and the mass ratio (i.e., the ratio of additional mass block mass to the frame mass) on the vibration reduction performance of NSBD are considered. The arch-span ratio and mass used in the experiment are listed in Tables2 and 3.Table 2
Arch-span ratio.
Arch height (mm)Beam span (mm)Ratio (%)1536000.42250.69350.97Table 3
Mass ratio.
Attachment mass (g)Mass of primary structure (g)Ratio (%)69.2650501.37111.62.21The experiment validated the attenuating effect of NSBD on the frame structure under sinusoidal and earthquake excitation. Firstly, the shaking table is used to perform a frequency sweep analysis of the experimental frame structure. The analysis determined that the natural frequency of the experimental frame structure is 2.89 Hz. Based on the switch threshold force of the buckling beams for different thicknesses, the amplitude of sinusoidal excitation applied to the frame structure is 0.05 cm. Secondly, the experiments used sinusoidal excitations to examine the damping provided by the NSBD for frame structures of different thicknesses. Finally, a total of 10 different earthquake excitations selected from two types of sites are used to investigate the vibration control effectiveness of the NSBD at the best damping. According to the code (GB5011-2010) [51], 10 different earthquake excitations are selected, as listed in Table 4, which represent different site conditions.Table 4
Information of the selected seismic records.
Site classificationEarthquake nameYearStationMagnitudeRrup (km)Vs30 (m/s)IImperial Valley-031951EI Centro Array #95.625.24213.44Southern Calif1952San Luis Obispo673.41493.5Imperial Valley-051955El Centro Array #95.414.88213.44Central Calif-021960Hollister City Hall59.02198.77Hollister-011961Hollister City Hall5.619.56198.77San Fernando1971Pasadena - Old seismo lab6.6121.5969.07IIGuy_2010-10-15-HHE2010X102 Temporary Station_ Greenbriar_ AR3.8612.26495.3Guy_2010-10-15-HHN2010X102 Temporary Station_ greenbriar_ AR3.8612.26495.3Mineral_2011-08-232011Keystone College_ La Plume_ PA5.74442.73689.9Sparks_2011-11-052011Smith Ranch_ Marlow_ OK4.73148.55548.4


## 5. Experiment Results

### 5.1. Vibration Control Effectiveness under Sinusoidal Excitation

Sweep sine excitations (0.05 displacement) ranged between 2.78 Hz and 3.0 Hz at 0.1 Hz intervals are used to test the effectiveness of the NSBD at and near the natural frequency. The frequency-displacement responses of the frame structure for the different configurations of the NSBD are shown in Figure7 and the corresponding time history responses are shown in Figure 8.Figure 7
Frequency-displacement response of frame structure with different parameters: (a)δ = 0.2 mm, μ = 1.37%, (b) δ = 0.2 mm, μ = 2.21%, (c) δ = 0.25 mm, μ = 1.37%, and (d) δ = 0.25 mm, μ = 2.21%.
(a)(b)(c)(d)Figure 8
Time-history response of frame structure under different parameters: (a)δ = 0.2 mm, μ = 1.37%, (b) δ = 0.2 mm, μ = 2.21%, (c) δ = 0.25 mm, μ = 1.37%, and (d) δ = 0.25 mm, μ = 2.21%.
(a)(b)(c)(d)As shown in Figure7, the inclusion of the NSBD reduced the vibrations of the frame structure especially for those under δ = 0.2 mm, μ = 2.21%, and h = 35 mm. The installation of the NSBD shifted the natural frequency far away from the dominant excitation frequency, thus providing a great damping effect between 2.86 Hz and 3.0 Hz. The vibration of frame structure can therefore be suppressed effectively by utilizing the NSBD, which not only mitigates vibrations an overall wide frequency range but also shifted the natural frequency far from the dominant frequency. Figure 8 illustrates the time-history response of frame structure with different parameters, and all detailed results are compared in Table 5.Table 5
Reduction of the peak displacement.
Arch heightδ = 0.2 mmδ = 0.25 mmμ = 1.37%μ = 2.21%μ = 1.37%μ = 2.21%w/o ctrl (mm)w ctrl (mm)Reduction ratiox (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)h = 15 mm46.8719.7157.9546.8713.7070.7746.8722.2152.6146.8730.5534.82h = 25 mm46.8717.4362.8146.8714.0570.0246.8725.3745.8746.8729.2137.68h = 35 mm46.8715.2867.4046.8715.5866.7646.8724.3847.8046.8731.4732.86Table5 shows that the maximum attenuation of the peak displacement reached up to 70.77%. Among the different tested configurations, the damping effect is best when the thickness of NSBD is 0.2 mm. An explanation is that when the thickness of the NSBD is large, a larger driving force is needed to initiate the snap-through motion between the two equilibrium positions of the NSBD. Therefore, under the same external excitation, NSBD with a lower thickness is more prone to snap-through motion, which improves vibration attenuation. This improvement can be seen in Figure 7(a) vs Figure 7(c) and Figure 7(b) vs Figure 7(d). For the same beam thickness, the difference in arching height leads to significant differences in kinetic energy consumed when snap-through occurs. In these experiments, NSBD with a large arch height has a larger displacement time history, which indicates that the vibration mitigation effect is also more obvious, as shown in Figures 7(a) and 7(b), the same in Figures 7(c) and 7(d). The attenuating effect of NSBD on the frame structure is more intuitively derived from its displacement time history, as shown in Figure 8, and its experimental regulation is similar to that in Figure 7.


### 5.2. Control Effectiveness under Earthquake Excitation

The values of the parameters for the NSBD are selected by the vibration suppression effect of the NSBD under sinusoidal excitation in the shake table tests. It is noted from Figures7(a) and 7(b) that the vibration suppression effect of the NSBD is better when its mass ratio is 2.21%. As shown in Figures 7(a) and 7(c), the vibration suppression effect of the NSBD is better when its thickness and height are 0.2 mm and 35 mm, respectively. The first set of tests uses six earthquake records selected according to Site classification I. To investigate the NSBD effectiveness with other site classes, four other earthquake records are selected according to Site classification II. The details of the selected earthquake excitations are listed in Table 4.Figure9 shows the displacement time histories of the frame structure (with and without optimized NSBD) during different earthquake excitations. The displacements of the frame structure are generally reduced when the NSBD is installed. The peak vibration reduction ratio (P) and RMS vibration reduction ratio are both dependent on seismic motions, demonstrating that NSBD can improve the structure’s resilience against hazardous earthquake ground motion. The peak vibration reduction ratio and RMS vibration reduction ratio can be obtained by equations (22) and (23), compared in Table 6. Table 6 shows that the peak damping ratio and average damping ratio can be attenuated respectively by up to 44.63% and 45.03%, demonstrating that the NSBD effectively suppresses the action of earthquake ground motion on the frame structure.(22)P=amax1amax2,where P refers to the peak vibration reduction ratio (%), amax1 refers to the maximum displacement amplitude of the structure during the earthquake excitation when the NSBD is installed (mm), amax2 refers to the maximum displacement amplitude of the structure during the earthquake excitation (mm).(23)Srms=1/N∑i=1na1i2=a112+a122+⋯+a1n2/N1/N∑i=1na2i2=a212+a222+⋯+a2n2/N,where Srms refers to RMS vibration reduction ratio (%), N refers to total time intervals, a1i refers to the displacement of the structure during the earthquake excitation at per time interval when the NSBD is installed (mm), and a2i refers to the displacement of the structure during the earthquake excitation at per time interval (mm).Figure 9
Time-history response of frame structure under different earthquake excitations: (a) Imperial Valley-03; (b) Imperial Valley-05; (c) Central Calif-02; (d) Mineral_2011-08-23; (e) Sparks_2011-11-05.Table 6
Vibration reduction ratio of NSBD to frame structure under different earthquake excitations.
Earthquake excitationsP (%)Srms (%)Imperial Valley-0325.1231.39Southern Calif26.8326.79Imperial Valley-0526.8925.89Central Calif-0225.1445.03Hollister-0132.8036.34San Fernando26.7729.82Guy_2010-10-15-HHE24.5840.50Guy_2010-10-15-HHN44.6324.58Mineral_2011-08-2328.8825.17Sparks_2011-11-0539.2111.25In summary, the experimental results show that the proposed NSBD has an excellent vibration mitigation effect on the frame structure under both sinusoidal excitation and earthquake excitations.


## 5.1. Vibration Control Effectiveness under Sinusoidal Excitation

Sweep sine excitations (0.05 displacement) ranged between 2.78 Hz and 3.0 Hz at 0.1 Hz intervals are used to test the effectiveness of the NSBD at and near the natural frequency. The frequency-displacement responses of the frame structure for the different configurations of the NSBD are shown in Figure7 and the corresponding time history responses are shown in Figure 8.Figure 7
Frequency-displacement response of frame structure with different parameters: (a)δ = 0.2 mm, μ = 1.37%, (b) δ = 0.2 mm, μ = 2.21%, (c) δ = 0.25 mm, μ = 1.37%, and (d) δ = 0.25 mm, μ = 2.21%.
(a)(b)(c)(d)Figure 8
Time-history response of frame structure under different parameters: (a)δ = 0.2 mm, μ = 1.37%, (b) δ = 0.2 mm, μ = 2.21%, (c) δ = 0.25 mm, μ = 1.37%, and (d) δ = 0.25 mm, μ = 2.21%.
(a)(b)(c)(d)As shown in Figure7, the inclusion of the NSBD reduced the vibrations of the frame structure especially for those under δ = 0.2 mm, μ = 2.21%, and h = 35 mm. The installation of the NSBD shifted the natural frequency far away from the dominant excitation frequency, thus providing a great damping effect between 2.86 Hz and 3.0 Hz. The vibration of frame structure can therefore be suppressed effectively by utilizing the NSBD, which not only mitigates vibrations an overall wide frequency range but also shifted the natural frequency far from the dominant frequency. Figure 8 illustrates the time-history response of frame structure with different parameters, and all detailed results are compared in Table 5.Table 5
Reduction of the peak displacement.
Arch heightδ = 0.2 mmδ = 0.25 mmμ = 1.37%μ = 2.21%μ = 1.37%μ = 2.21%w/o ctrl (mm)w ctrl (mm)Reduction ratiox (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)w/o ctrl (mm)w ctrl (mm)Reduction ratio (%)h = 15 mm46.8719.7157.9546.8713.7070.7746.8722.2152.6146.8730.5534.82h = 25 mm46.8717.4362.8146.8714.0570.0246.8725.3745.8746.8729.2137.68h = 35 mm46.8715.2867.4046.8715.5866.7646.8724.3847.8046.8731.4732.86Table5 shows that the maximum attenuation of the peak displacement reached up to 70.77%. Among the different tested configurations, the damping effect is best when the thickness of NSBD is 0.2 mm. An explanation is that when the thickness of the NSBD is large, a larger driving force is needed to initiate the snap-through motion between the two equilibrium positions of the NSBD. Therefore, under the same external excitation, NSBD with a lower thickness is more prone to snap-through motion, which improves vibration attenuation. This improvement can be seen in Figure 7(a) vs Figure 7(c) and Figure 7(b) vs Figure 7(d). For the same beam thickness, the difference in arching height leads to significant differences in kinetic energy consumed when snap-through occurs. In these experiments, NSBD with a large arch height has a larger displacement time history, which indicates that the vibration mitigation effect is also more obvious, as shown in Figures 7(a) and 7(b), the same in Figures 7(c) and 7(d). The attenuating effect of NSBD on the frame structure is more intuitively derived from its displacement time history, as shown in Figure 8, and its experimental regulation is similar to that in Figure 7.


## 5.2. Control Effectiveness under Earthquake Excitation

The values of the parameters for the NSBD are selected by the vibration suppression effect of the NSBD under sinusoidal excitation in the shake table tests. It is noted from Figures7(a) and 7(b) that the vibration suppression effect of the NSBD is better when its mass ratio is 2.21%. As shown in Figures 7(a) and 7(c), the vibration suppression effect of the NSBD is better when its thickness and height are 0.2 mm and 35 mm, respectively. The first set of tests uses six earthquake records selected according to Site classification I. To investigate the NSBD effectiveness with other site classes, four other earthquake records are selected according to Site classification II. The details of the selected earthquake excitations are listed in Table 4.Figure9 shows the displacement time histories of the frame structure (with and without optimized NSBD) during different earthquake excitations. The displacements of the frame structure are generally reduced when the NSBD is installed. The peak vibration reduction ratio (P) and RMS vibration reduction ratio are both dependent on seismic motions, demonstrating that NSBD can improve the structure’s resilience against hazardous earthquake ground motion. The peak vibration reduction ratio and RMS vibration reduction ratio can be obtained by equations (22) and (23), compared in Table 6. Table 6 shows that the peak damping ratio and average damping ratio can be attenuated respectively by up to 44.63% and 45.03%, demonstrating that the NSBD effectively suppresses the action of earthquake ground motion on the frame structure.(22)P=amax1amax2,where P refers to the peak vibration reduction ratio (%), amax1 refers to the maximum displacement amplitude of the structure during the earthquake excitation when the NSBD is installed (mm), amax2 refers to the maximum displacement amplitude of the structure during the earthquake excitation (mm).(23)Srms=1/N∑i=1na1i2=a112+a122+⋯+a1n2/N1/N∑i=1na2i2=a212+a222+⋯+a2n2/N,where Srms refers to RMS vibration reduction ratio (%), N refers to total time intervals, a1i refers to the displacement of the structure during the earthquake excitation at per time interval when the NSBD is installed (mm), and a2i refers to the displacement of the structure during the earthquake excitation at per time interval (mm).Figure 9
Time-history response of frame structure under different earthquake excitations: (a) Imperial Valley-03; (b) Imperial Valley-05; (c) Central Calif-02; (d) Mineral_2011-08-23; (e) Sparks_2011-11-05.Table 6
Vibration reduction ratio of NSBD to frame structure under different earthquake excitations.
Earthquake excitationsP (%)Srms (%)Imperial Valley-0325.1231.39Southern Calif26.8326.79Imperial Valley-0526.8925.89Central Calif-0225.1445.03Hollister-0132.8036.34San Fernando26.7729.82Guy_2010-10-15-HHE24.5840.50Guy_2010-10-15-HHN44.6324.58Mineral_2011-08-2328.8825.17Sparks_2011-11-0539.2111.25In summary, the experimental results show that the proposed NSBD has an excellent vibration mitigation effect on the frame structure under both sinusoidal excitation and earthquake excitations.


## 6. Conclusions

The negative stiffness bistable damper (NSBD) is currently proposed to dampen the vibration response of the frame structure. The constitutive equation of the NSBD is derived to analyze the effects of stiffness ratio, arch-span ratio, and damping ratio on its restoring capabilities. To verify the vibration suppression performance of the NSBD, shaking table tests of a frame structure are conducted under different sinusoidal and seismic excitations with and without an NSBD. The results reveal that the NSBD can effectively restrain the displacements of the structure. In addition, the buckling beam cannot have significant deflection to allow it to vibrate between two equilibrium positions, meaning that the beam must remain in its elastic range for the NSBD to be effective. The main findings in the present study are summarized as follows:The proposed damper can greatly mitigate the peak displacement of the frame structure at the dominant frequency and has a wide damping bandwidth. Furthermore, the control effectiveness of the NSBD is robust against a range of sinusoidal and seismic excitations, demonstrating that the NSBD can strengthen structural security against hazardous seismic events.This paper only explored the attenuating effect of NSBD on a single-story frame structure experimentally. Future works will be undertaken to perform related analytical and numerical research. Furthermore, the vibration control of structures using the NSBD in actual engineering will be further studied and improved in future work.



---

*Source: 6397602-2022-11-01.xml*