Design and Simulation of a Damper with Negative Stiffness for Vibration Mitigation from Drilling Equipment to a Semi-Submersible Platform

A negative-stiffness damper (NSD) that incorporates eccentric columns and high-stiffness springs is proposed for vibration mitigation from drilling equipment to a semi-submersible platform. Eccentric columns are considered as negative-stiffness elements, which can provide high static stiffness and achieve negative stiffness with hysteresis energy dissipation caused by buckling mode transition under dynamic conditions. With the assistance of high-stiffness springs, the damper possesses a hysteresis loop that can dissipate vibration energy and mitigate responses. Firstly, the hysteresis characteristics and negative stiffness of eccentric columns are analyzed. ,en, the design principle of the NSD is derived, and several characteristics required for applications are investigated. In addition, an NSD system is designed for a target drilling mud pump based on measured data, and the mitigation performance is analyzed by comparing the vibration responses between cases of the proposed NSD system and the steel pedestal. ,e results demonstrate that the proposed system can effectively reduce the width of affected coverage and the magnitudes of vibration responses. ,e investigations performed in this paper not only provide the design principle of the new damper, but also prove the great potential of the proposed damper in vibration mitigation from drilling equipment to the hull of the platform.


Introduction
Vibrations generated by large-scale equipment on a semisubmersible platform are attracting more and more attentions because they can easily propagate through cabins and may cause harm to the health of the personnel and even cause structural damage [1]. Steel pedestals are commonly used to support large-scale equipment due to the high bearing capacity. However, the low damping of steel pedestals may lead to poor effects on vibration mitigation [2,3]. erefore, a damper with both high bearing capacity and energy dissipation is required to mitigate the responses generated by large-scale equipment.
In recent years, nonlinear dampers have drawn increasing attentions since they possess several advantages compared with linear ones. Ibrahim [4] summarized the characteristics of many different types of nonlinear dampers, in which a group of dampers using buckled columns are being analyzed. When buckled columns are constrained laterally, they exhibit negative stiffness by providing a reaction force in the same direction as the displacement under compressions. Initially, buckled columns were used as auxiliary components in damper design. Liu et al. [5] and Huang et al. [6] proposed nonlinear dampers containing Euler buckled columns, which can be regarded as correctors with negative stiffness. A new damper using buckled columns with fixed ends was proposed and analyzed in [7][8][9]. Brennan et al. [10] and Fulcher et al. [11] proposed a tuned damper consisting of two buckled columns. e mitigating frequency range can be tuned by adjusting the buckling shape of the columns.
Afterwards, a new axial damper of flat-ended column was proposed. e flat-ended column possesses high static stiffness. When axial force of a column reaches critical load, the negative stiffness is generated through the buckling mode transition from buckling with fixed ends to pinned ends. During a loading-unloading process, the axial force-deformation curve of the flat-ended column shows a flag-shaped hysteresis loop, which makes the structure suitable for negative stiffness elements in the design of dampers. Kalathur et al. [12] analyzed the buckling behavior of compressed flat-end columns theoretically. Dong et al. [13][14][15] proposed dampers of columns. Both high damping and high stiffness are achieved by these dampers. Zhang et al. [16] proposed a convex eccentric column brace with wider applicability. Wu and Phillips [17] designed a shape memory alloy (SMA) capped column with elastic buckling mode jump (BMJ) mechanism, obtaining flag-shaped hysteretic damping without material failure. Wang et al. [18] designed an all-metal brace with hysteresis dissipation for impact protection of jacket platforms.
Although the high stiffness and high damping of eccentric columns are attractive for the vibration mitigation of large-scale equipment, a method of harmonizing the damper parameters with the specific equipment and dynamic loads is required before applications. In addition, the mitigation efficiency of the damper utilizing these elements in the cabin case has not been analyzed.
In this paper, a damper incorporating both eccentric columns and high-stiffness springs is proposed and investigated. e hysteresis behavior and multistage stiffness of eccentric columns are analyzed, and the design principle and dynamic analysis of the damper are clarified. Based on the design principle, a negative-stiffness damper (NSD) system is designed for a target drilling mud pump. e buckling behavior of dampers is well captured. Vibration responses of the structure within the cases of the system and the steel pedestal are compared. e results demonstrate that the negative-stiffness damper system not only greatly decreases the amplitude of vibration responses, but also effectively reduces the width of the affected coverage.

Numerical Simulation of the Hysteresis Behavior of an Eccentric
Column. e configuration of an eccentric column is shown in Figure 1. e proposed eccentric column can be regarded as a rectangular column connected to two thick ends. A quarter of circle arc is carved for avoiding stress concentration at the junction between each end and the main column. e thickness of junction equals the difference value between b e and b. e eccentricity of an eccentric column is obtained by such difference value. e buckling behavior of a sample column is analyzed by numerical simulation. A three-dimensional (3D) model consisting of two identical base plates and the sample column is established by utilizing the finite-element program, Abaqus 6.14. Geometric properties of the sample column are listed in Table 1. Both sizes of two base plates are 30 mm × 30 mm × 5 mm. All parts are assigned as the material of steel (60SiMnA), and linear elastic property is applied in the simulation. e material properties are listed in Table 2.
All parts of the model are meshed with the element type of 8-node linear brick with reduced integration (C3D8R), and the mesh size of the column is set as 0.4 mm according to the convergence study. For capturing the buckling phases of the sample column, the boundary conditions of two base plates are assigned as shown in Figure 2. All contact surfaces between ends and plates are modeled as roughly contacted and two surfaces are allowed to separate. e damping of the applied material, 0.004, is transformed to the form of Rayleigh damping and applied.
Instead of a harmonic force, a harmonic displacement control is loaded on the upper base plate to avoid structural failure, which may be caused by negative stiffness. e amplitude of the harmonic displacement is 0.1 mm and the frequency is 1 Hz. e implicit dynamic analysis is used in the simulation because of its advantage in dealing with nonlinear problems. e buckling mode transition process of an eccentric column in a loading-unloading cycle is summarized in Figure 3.
Buckling with fixed ends is generated when axial force reaches the first critical load, as shown in Figure 3. e boundary condition of both ends transforms from fixed to pinned mode when axial force reaches the second critical load. en, the eccentric column goes backward to upright state directly when unloading. Time histories of axial force and deformation and the relationship between them and maximum stresses under two modes are shown in Figures 4 and 5. e time histories of the axial force and deformation of the sample column in a cycle are merged and shown in Figure 4(a), and the relationship between axial force and deformation is shown in Figure 4   Shock and Vibration axial force at points 3 and 6 represent the moments when the buckling mode changes from ends fixed to pinned and then go back to upright state, respectively. It indicates that point 3 (d 3 , P 3 ) and point 6 (d 1 , P 0 ) can be considered as the trigger point and the attenuation point of hysteresis behavior, respectively. A D-shaped loop is formed as shown in Figure 4(b), whose area represents the energy dissipation of the column during a cycle. Figure 5(a) represents the buckling mode with fixed ends triggered when axial force reaches the first critical load, P 2 (in Figure 4(b)). Figure 5(b) shows the buckling shape of the eccentric column when subjected to the peak displacement. e zoom-in windows clearly show the contact modes between ends and plate surfaces. Figure 5(b) shows the maximum stress in a cycle, which is far below the material yielding stress of the material.
Simulation results demonstrate that the hysteresis behavior of the sample column can be triggered in a small displacement range. e excellent bearing capacity and the large margin between the maximum stress and yielding stress proves the high applicability of the eccentric column.

Multistage Stiffness of Eccentric Columns.
e multistage stiffness of eccentric columns is an external expression of the hysteresis loop, which is directly related to the instantaneous bearing capacity. erefore, when designing an NSD, the key points associated with the loop (such as points 1, 2, 3, 4, and 6) are significant factors in evaluating the stiffness of a column at different buckling phases. e buckling phases and the multistage stiffness of the sample column are summarized and shown in Figure 6. e hysteresis loop is divided into four phases according to the buckling mode transition-a: 0-2 (upright), b: 2-3 (buckling mode with fixed ends), c: 3-4-5-6-1 (buckling mode with pinned ends), and d: 1-0 (upright). Firstly, the stiffness value of column in phase a is 21463 N/mm and then gradually decreases to 0 N/mm in phase b. Negative stiffness is triggered when the axial force reaches P 3 (in Figure 4). When the boundary condition of ends completely transforms to pinned, the axial stiffness maintains at 0 N/mm as the axial force in this phase remains the same. Negative stiffness is also triggered at the moment of backward transition, and then the stiffness goes back to positive when the axial force returns to P 1 (in Figure 4) again.
Two buckling modes can be regarded as surface-to-surface and edge-to-surface mode, respectively. e corresponding drastic reduction in the contact area leads to a decrease in the bending stiffness of the column during loading. e restoring of bending stiffness is conducted with the increase of contact area during unloading. Two transitions in phases b and c represent the moments when the bending stiffness is too weak to support the compression due to the increase in deflection and when the increase of the bending stiffness is strong enough to outstrip the axial force reduction caused by the decrease in deflection, respectively. e abrupt snap effect converts the low-frequency input to a much higher-frequency one. Energy loss mechanisms are more effective at a higher frequency [13]. erefore, with extra stiffness provided by auxiliary components, eccentric columns are able to support and mitigate a large-scale equipment. e theoretical model summarized in [17], as listed in Table 3, can be applied to verify the simulation results and predict the stiffness of eccentric columns in damper design. Table 3 lists the axial forces and stiffness of eccentric columns in different phases.
In Table 3, E is the elastic modulus, and A is the area of the column section. e deflection function from point 2 to 3 is v 2 � (v max /2)(1 − cos(πz/L b )), and v max � (b/2). z is the axial coordinate. e stiffness of an eccentric column between point 0 and point 2, K 1 , represents the initial bearing capacity. But when slenderness ratio (L b /b) exceeds 100, the value of K 2 can be regarded as 0 N/mm because of the little difference between P 3 and P 2 . Figure 6(a) shows the theoretical and simulation results of the sample column. Although the discrepancies at point 2 and point 6 are due to the ignorance of stiffness varying in buckling mode transition processes, good agreements between two loops can still verify the simulation result.        Shock and Vibration to the relationship among geometric parameters of columns, positive-stiffness elements, and exciting force. en, the effects of several factors on hysteresis energy of the NSD and responses of the bottom structure are investigated for providing more information for applications.

e Design Principle of the NSD.
e NSD incorporates both eccentric columns and high-stiffness springs, as shown in Figure 7. Eccentric columns can be considered as negative-stiffness elements, while high-stiffness springs can be considered as positive-stiffness elements. Several columns and springs are all arranged in the cavity of nested sleeves. In addition, ends of each column are assembled into the buckling-guiding grooves of base plates for restricting sliding, as shown in the zoom-in window of Figure 7. In this way, buckling in the specified direction is triggered when columns are compressed.
When supporting a piece of equipment, eccentric columns remain upright in equilibrium, and compression forces of both positive-and negative-stiffness parts hold the mass of equipment. When the supported equipment is running, the upper base plate starts to fluctuate up and down. Negative stiffness of each eccentric column would be triggered along with the buckling mode transition, resulting in a decrease in the axial stiffness of the NSD. e dynamic load originally carried by eccentric columns would mostly transfer to high-stiffness springs to maintain the stability of the whole damper. In the unloading process, eccentric columns return to the upright state directly, and the axial stiffness of the whole damper returns to the initial stage.
An NSD consists of several identical eccentric columns and high-stiffness springs to support large-scale equipment. e parameters and numbers of each element should be determined according to the exciting force of the target equipment. It is assumed that the mass of the supported equipment is M, and it generates an exciting force with the amplitude of P A and the rotated speed of n (the unit of n is r/s). e exciting force curve of the equipment can be written as a sinusoid: (1) e first task is to determine the number and the configuration of eccentric columns. e relationship between the geometric parameters and the exciting force can be analyzed based on the recommended locations of significant points of columns as shown in Figure 8(a). e trigger point and attenuation point of eccentric columns should be located at the marked regions of the exciting force curve, which means the negative stiffness should be triggered before the trough and ceased before the crest. e positive-and negative-stiffness parts of the damper are composed of a number of parallel highstiffness springs (with total stiffness of K S ) and several identical eccentric columns, respectively. Besides, the high-stiffness springs can be precompressed by part of equipment weight (M P ) to shorten the axial displacement of the supported equipment. To ensure the triggering of hysteresis loop, the weight of equipment should surpass the maximum bearing capacity of all columns (that is nP 3 < Mg − M P g), and the minimum of dynamic load should surpass the residual bearing capacity of buckled columns (that is nP 4 erefore, the relationship between geometric properties of eccentric columns and the weight of equipment can be written as follows: (2) e second task is to determine the total stiffness of all high-stiffness springs, which can be inferred from the relationship between the static equilibrium position and hysteresis loop of columns, as shown in Figure 8 As high-stiffness springs are important elements to support the running equipment when negative stiffness of columns is triggered, K s should be determined based on K 1 and dynamic loads. Initially, the stiffness of negative stiffness part is nK 1 (n is the number of eccentric columns). erefore, the total stiffness of the damper at the static equilibrium position, K SEP , is calculated by the sum of two parts: nK 1 + K S .

Phase
Axial force Stiffness Eccentric columns Sleeves High-stiffness spring Shock and Vibration e axial force of a column at static equilibrium position, P SEP , is On the one hand, the static equilibrium position should be located between point 0 and point 3, so the value of P SEP should be less than P 3 . On the other hand, the range of the hysteresis loop should include the acting range of the exciting force. erefore, total spring stiffness, axial forces, equipment weight, and exciting force amplitude should meet the following requirements: e relationship between the value of K S and the column factors can be summarized as follows:

Hysteresis Loop of the NSD.
An NSD model utilizing the sample column (described in Section 2.2) is designed according to a specific exciting force. e hysteresis loop of the whole damper and the effects of several factors on the shape of the loop are analyzed. When the 3D column model shown in Figure 2 is subjected to an exciting force, which is written as F E � 3200 − 1800 sin(πt), four springs with stiffness of 4500 N/mm (K S �18000 N/mm) are assigned between the upper and lower base plates according to (5). erefore, the sample column and four springs form an NSD model. Boundary conditions remain the same as described in Section 2.2. e exciting force is loaded to the upper base plate, and implicit analysis is also adopted to capture the axial force and deformation of the NSD model. e axial force-deformation relationship of the NSD model is shown in Figure 9(a), and the time histories of axial force and deformation of each spring are integrated in Figure 9(b) due to the constant stiffness. e clear loop demonstrates a considerable energy dissipation of the NSD model. Compared to the hysteresis loop of the single eccentric column, the loop of the NSD model can be considered as obliquely stretched due to the enlarged axial static stiffness. Initially, eccentric columns remain as the major bearing parts before the first critical load is reached. When the negative stiffness of columns is triggered at point 3, the stiffness drop of the model causes an enlargement in axial displacement of the equipment. en, the major bearing part converts to the springs, but eccentric columns also participate in the bearing with the residual axial force, P 0 . A slight resonance process at point 4, which cannot be analyzed by theoretical model, is motivated by buckling mode transition as shown in the dotted boxes of Figure 9(b). e results demonstrate that changing the stiffness of springs can adjust both of the deformation at static equilibrium position and maximum deformation of the damper. e geometric parameters of the columns can affect the stiffness and the positions of key points of the NSD model [17]. erefore, the configuration of the columns and K S can be adjusted to change the characteristics of the NSD model. To clarify the effects of the important parameters on the hysteresis loop of the NSD model, the hysteresis loops of NSD models with different K S , column slenderness ratio, and depth ratio are simulated by altering one of them at a time, respectively. In addition, to figure out the rate of change in each parameter which can increase the energy dissipation under the same load, the value of K S , column slenderness ratio, and depth ratio are changed by 1600 N/ mm, 4, and 0.8 twice, respectively. e hysteresis loops are shown in Figure 10. K S has a little effect on the initial stiffness and energy dissipation because they are not directly involved in buckling  hysteresis dissipation, as shown in Figure 10(a). e axial force of the trigger point is proportional to K S , and changing K S can effectively adjust the maximum deformation. e slenderness ratio can significantly increase the energy dissipation, but shows little effect on the stiffness of the damper, as shown in Figure 10(b). e axial force of trigger point can be significantly increased by decreasing slenderness ratio, as P 3 of a column is proportional to the cube of b. In terms of depth ratio, only a great increase in depth ratio can bring about an obvious enlargement both in stiffness and in energy dissipation, as shown in Figure 10(c). Both of the stiffness and energy dissipation of an eccentric column can be intensified by decreasing slenderness ratio or increasing depth ratio. Figure 10 also demonstrates that all factors are associated with the stiffness between point 0 and point 4, but only K S can affect the damper stiffness from point 5 to point 6. Effects of the three parameters on design factors are summarized as shown in Figure 11. e energy dissipation of the column can be enlarged by increasing the slenderness and depth ratio, as shown in Figure 11. Axial force of the trigger point can be magnified by increasing all three parameters. Axial deformation of the trigger point can only be changed by adjusting slenderness ratio. e increase in K S leads to a negative effect on energy dissipation, which is due to more potential energy absorbed by springs, resulting in smaller loops [18]. In general, slenderness ratio is the most significant parameter due to its strong effect on all three factors. Depth ratio follows due to its appreciable effects on energy dissipation and axial force of the trigger point. K S can be considered when adjusting maximum axial deformation.

e Effect of Vibration Frequency on NSD.
e resonance process mentioned in Section 3.2 can be affected by vibration frequency directly, and the intensity of the resonance process may cause the distortion of hysteresis loop, leading to a loss in energy dissipation. To clarify the influence of vibration frequency on the resonance process, the NSD model described in Section 3.2 is assumed to be loaded by harmonic forces with the same amplitude of 2500 N and four exciting frequencies, 2.5 Hz, 5 Hz, 7.5 Hz, and 10 Hz. e hysteresis loops are illustrated in Figure 12.
Hysteresis behavior of the NSD model is triggered under exciting forces at all four frequencies, which basically verifies the effectiveness of this model under loads at 2.5∼10 Hz. e resonance process gets more intense and the maximum deformation develops as the vibration frequency increases. e influence of exciting frequency on displacement range of the resonance process and the Shock and Vibration 7 maximum deformation of the NSD model are summarized as shown in Figure 13. e displacement ranges of the resonance process and the maximum deformation of the NSD model are all increased along with the growth of exciting frequency. e histogram in Figure 13 demonstrates that the displacement range is in proportion to the vibration frequency, which means that the duration of resonance process of the column does not change along with the exciting frequency. Line chart in Figure 13 indicates that although the maximum deformation of the NSD model is enlarged, the increment gradually gets smaller. e enlargement in maximum deformation means more potential energy is absorbed. erefore, a smaller increment in axial deformation may result in a poor performance in vibration mitigation at a higher frequency. When designing an NSD, a reasonable margin should be retained between the deformation of trigger point and maximum deformation to avoid excessive instability in practical applications.

e Force Transmissibility of the NSD.
e force transmissibility is the basic index of mitigation performance of a damper. However, the force transmissibility of the NSD can be affected by the resonance process, which cannot be predicted by theoretical model. erefore, numerical simulations are performed to analyze the force transmissibility of the NSD model.
To obtain the transmitted force of the NSD model, a thick bottom plate with a size of 600 mm × 600 mm × 9 mm is modeled under the lower base plate. In this case, the lower base plate is tied to the thick bottom plate, and the boundary condition of edges is pinned. e force transmissibility (T) is defined as the ratio of the amplitude of the transmitted force on bottom plate to the exciting force amplitude. e force transmissibility can be written in decibels as follows: where f t is the amplitude of the transmitted force, and f e is the amplitude of exciting force. e target mitigation frequency range of the NSD is concentrated in the range of 1 Hz to10 Hz, so only the desired frequency range (1 Hz-10 Hz) is discussed. Exciting force inputs at 1 Hz to 10 Hz are assigned to the center of the upper base plate, and implicit dynamic analysis is applied. e force transmissibility curves of this NSD model within force amplitudes of 2000 N, 2500 N, and 3000 N and at the exciting frequency of 1 Hz to 10 Hz are obtained as shown in Figure 14.
All three curves keep horizontal at 1 to 5 Hz and then slightly increase at 5 to 10 Hz, which demonstrates the stability of the NSD model in the target frequency range. e mitigation performance of the model depends on both the potential energy absorbed by column deformation and the hysteresis dissipation. However, the resonance process can affect the hysteresis loop when subjected to an exciting force at a higher frequency, and that is the reason why the transmissibility remains the same at 1 to 5 Hz. e transmissibility is decreased as more energy is absorbed through larger deformation of columns. erefore, the mitigation efficiency can be improved when subjected to a relative larger exciting force, but the increase rate is decreasing.
e turning points of transmissibility curves, which represent the critical points of mitigation efficiency, move to a higher frequency when subjected to a harmonic force with larger amplitude, due to a larger margin obtained between the trigger point and the maximum deformation.
In terms of the effect of geometric parameter on force transmissibility, the slenderness ratio, which is the most significant parameter on the hysteresis behavior, is analyzed. e transmissibility curves of another two cases with slenderness ratios of 134 and 137.5, which are increased by 0.69% and 3.4%, respectively, are analyzed and compared to those of the original case as shown in Figure 15.
e slightly enlarged slenderness ratio shows positive effect on the force transmissibility. e more margin retained between the trigger point and maximum deformation, the larger transmissibility and the higher frequency the turning point located. Although the energy dissipation is decreased, the augment in potential energy caused by the decreased stiffness when increasing slenderness ratio brings an improved transmissibility. However, the trend is weakening due to the nonlinearity between the slenderness ratio and P 3 of the column. It is demonstrated that pursuing energy dissipation maximization may cause a secondary problem of poor transmissibility when subjected to dynamic load at a higher frequency. erefore, when designing an NSD, a reasonable margin between trigger point and maximum deformation should be retained according to the target frequency.

Numerical Simulation of the NSD System Designed for a Drilling Mud Pump
An NSD system is designed to mitigate the vibration from a drilling mud pump to platform structure, and its mitigation efficiency is investigated by comparing the responses of the NSD system case and the steel pedestal case.

Excitation of the Drilling Mud Pump.
Drilling mud pump is an important oil-mining equipment, which can deliver mud from the surface to drilling bit for cooling it and transfer the cuttings back to the surface. However, drilling mud pumps for deep-sea oil exploitation operate at low rotated speeds and they are easily to generate low-frequency exciting force with large amplitude. e target drilling mud pump is a triplex piston pump with the size of 5543 mm × 3290 mm × 4528 mm and a weight of 39000 kg. e drilling mud pump possesses the rotated speed of 1.75 r/s. When the pump starts, a motor motivates the rotation of the pinion and then triggers the rolling of crankshaft assembly. e rolling motion transforms to the unidirectional reciprocating motion of the crossheads through the connection rods. erefore, the exciting force of the cylinder is mainly composed of inertia forces of rolling crankshaft assembly, which consists of three components: crankshaft, pinion, and crosshead [19]. e target pump contains three cylinders, and the phase difference between each cylinder during operation is (2π/3). To obtain the exciting force of the target drilling mud pump, the acceleration of the steel pedestal is measured on a semisubmersible platform as shown in Figure 16(a), and the frequency spectrum of the measured data is shown in Figure 16(b). e measured acceleration contains several-order components in the frequency range of 0∼15 Hz, where the components at 1.75 Hz and 5.25 Hz represent the excitation loads generated by the single cylinder and the whole pump, respectively. e measured acceleration can be transformed to velocity by numerical simulation, and the exciting force of the pump can be estimated according to [20]. e exciting force under rotated speed of 1.75 r/s can be obtained as shown in Figure 17.

e Design of the NSD System.
In the traditional arrangement, the drilling mud pump is placed on steel pedestals welded on the double bottom in pump cabin. e length, width, and height of the cabin are 15 m, 15 m, and 1.8 m, respectively. e steel pedestal consists of several webs and panels. e new proposed system consists of nine dampers in the arrangement of 3 × 3, as shown in Figure 18. All dampers are fixed on the double bottom. It is assumed that all dampers are of the same configuration as they are evenly loaded. e hysteresis behavior of dampers is assigned to be triggered and attenuated at two dotted lines as shown Figure 17, for providing the desired energy dissipation, as well as sufficient bearing capacity and margin between the trigger point and the peak. e details of the single damper are shown in Figure 19.
e NSD consists of two sleeves, eight eccentric columns, and four high-stiffness springs. To display all internal elements, the sleeves, which are utilized to restrict lateral displacement, are set as transparent.
e geometric parameters of eccentric columns are shown in Table 4. e total stiffness of springs is 380,000 N/mm, and K 1 of the column is 125,000 N/mm. For this damper, K s is set to be 380,000 N/mm, and the whole damper can provide the static stiffness of 1,380,000 N/mm. e hysteresis behavior of the single NSD model is simulated on the double bottom model, and the force--deformation relationship is simulated through implicit dynamic analysis to verify the energy dissipation and the effect of resonance process of the proposed damper in a cycle. e double bottom of the pump cabin is established as shown in Figure 18. e double bottom model consists of 23307 quadrilateral elements and 9417 line elements, whose sizes are all 300 mm. Boundary conditions of the edges are assigned as simply supported. In terms of the exciting force inputs, 33,000 N is preloaded to springs of the damper to shrink the axial displacement of trigger point, and the exciting force shown in Figure 17 is divided and loaded on the upper base plate of damper in the form of pressure. e force of an eccentric column and the force-deformation relationship of the damper model under a cycle are shown in Figure 20. e hysteresis behavior of the NSD model is quite stable. e negative stiffness was not triggered at 2 s because the exciting force at this moment did not exceed the trigger line. Force-deformation relationship of the NSD model is shown in Figure 20(b); the maximum deformation of the upper base plate only reaches 0.327 mm, and the maximum stress of columns is 100.9 MPa in a cycle. e energy dissipation of an eccentric column is 418.6 N·mm, and the whole dissipation of an NSD is 3348.8 N·mm.

Mitigation Efficiency of the NSD System.
Steel pedestals and the NSD system are modeled at the center of the cabin. e mesh size of steel pedestals is 10 mm, and that of eccentric columns is 0.25 mm. Exciting force is evenly assigned to the panels of two pedestals in the form of pressure. For the NSD system, each damper is assigned to the exciting force described in the last section. e total analysis time is 60 seconds, and the steady-state responses of the double bottom of the two cases are shown in Figure 21.
e results demonstrate that the NSD system can effectively decrease the maximum displacement by 75% and shrink the width of affected coverage by 46%. Steel pedestals not only possess little energy dissipation, but also centralize vibration transmission, which makes it easy to induce high response over the structure. e maximum response region is divided into two parts, which are located on the front and rear of the pump, respectively. Responses of these two regions are superimposed by displacements transmitted through surrounding dampers. On the one hand, the weight supported by the single damper is reduced. On the other hand, the deformation of columns greatly alleviates the excitation transmitted to the hull.
Responses of four test points are captured and shown in Figure 22. Time histories of four test points in the last 5 seconds are shown in Figure 23. e displacements of static equilibrium position and amplitudes of all responses are decreased significantly, and the displacement decreases as the distance from the test point to the center of the cabin increases. Time history of steel pedestal case in Figure 23(a) represents the maximum displacement responses in the traditional arrangement. Because of the potential energy absorbed by column buckling and the hysteresis dissipation generated by buckling mode transition, the maximum displacement is decreased by 75%, and the displacement amplitude is shortened by 70%.
Time histories of test points 2 and 3 as shown in Figures 23(b) and 23(c) represent the maximum displacement responses in the NSD case. e equilibrium position displacement is decreased by 0.2 mm, and the amplitude is shortened by 0.1 mm, as shown in Figure 23(b). However, the mitigation efficiency is limited for test points 2 and 3. e mitigation efficiencies of peak displacement and response amplitude are 66% and 60%, respectively.
Compared to the responses of the other three test points, response 4 in both cases is all relatively slight. e displacement of equilibrium position in the NSD case slightly decreases. e mitigation efficiencies of peak displacement and response amplitude at test point 4 are 60% and 50%, respectively.
Fourier transform is conducted on four responses to analyze the frequency spectrums. Since the frequency components of all four responses in each case are basically identical, the frequency spectrum of response 1 is taken as an example, as shown in Figure 24. e frequency spectrum in Figure 24 demonstrates that response 1 in steel pedestal case contains a tiny component of 1.75 Hz and a large component of 5.25 Hz. Compared to the frequency spectrum of excitation shown in Figure 16, it is clear that the two components at 1.75 Hz and 5.25 Hz are due to the external excitation.
In terms of the response in the NSD system case, it contains a large component of 6.7 Hz and a smaller component of 7.3 Hz. e 6.7 Hz component is induced from the NSDs, which can be proved by the frequency spectrum of the transmitted force from an eccentric column to the double bottom, as shown in Figure 25. e 7.3 Hz component of response 1 is superimposed by the 6.7 Hz components generated by the surrounding NSDs.

Steel pedestal
The drilling mud pump NSD system Figure 18: e schematic diagrams of the steel pedestal and the proposed NSD system.

Shock and Vibration 13
When the exciting force shown in Figure 17 is loaded on the NSDs, the exciting force can be converted to an output excitation due to the multistage stiffness of eccentric columns. And responses are induced on the double bottom of the cabin by such output excitation. e output excitation contains a large component of 6.7 Hz, as well as three small ones at 11.4 Hz, 15.1 Hz, and 22.5 Hz, respectively. e component of lowest frequency of 1.75 Hz in the exciting force is eliminated in the output excitation due to the buckling mode transition. e 6.7 Hz component is converted from the 5.25 Hz component in exciting force. e original cycle is shortened by removing the duration of resonance process in buckling mode transition. e component of 11.4 Hz is generated by the corresponding component of the exciting force during phase a and phase c of eccentric columns, and the component of 15.1 Hz is generated by the buckling mode transition between phases c and d. e component of 22.5 Hz represents the maximum amplitude of the resonance process. e results demonstrate that the response frequency range of the NSD model is shifted higher than that of the excitation. And the mitigation efficiency is improved because energy loss mechanisms is more effective at higher frequencies.
e response of three high-frequency components in frequency spectrum of the output excitation (as shown in Figure 25) cannot be found evidently in the frequency spectrum of response 1 (as shown in Figure 24). Because the response amplitudes transmitted to the bottom can be alleviated significantly, the mitigation efficiency of the vibration energy is improved to a higher-frequency range. e investigation can preliminarily verify the potential utilization of the eccentric column in cabin structure, and the results also demonstrate that this NSD system can effectively mitigate the vibration responses transmitted from drilling equipment to platform structures.

Conclusion
is paper describes a kind of damper that incorporates both eccentric columns and high-stiffness springs, where the two parts are working as the positive-and negative-stiffness parts, respectively. By allowing the ends of eccentric columns to tilt, the boundary conditions of the two ends can change from fixed to pinned, which generates hysteresis damping. e buckling of an eccentric column under harmonic displacement control is simulated. Based on the four phases distinguished by buckling modes, its hysteresis loop and multistage stiffness are theoretically analyzed. e results not only verify the accuracy of simulation, but also confirm the applicability of eccentric columns in damper design. e design principle of an NSD is derived according to the relationship between parameters under the condition of ensuring the trigger of hysteresis behavior. e hysteresis    loop of an NSD model under harmonic force is captured, and the influences of key factors, such as column configuration and stiffness of springs, are analyzed to provide more information for damper design.
Furthermore, an NSD system is designed for a drilling mud pump. e mitigation efficiency of the system is compared to the steel pedestal. Because of the abrupt snap effect and the hysteresis dissipation generated by the buckling mode transition, the NSD system can effectively decrease the maximum displacement by 75% and reduce the affected coverage on the structure by 46%. e investigations performed in this paper not only clarify hysteresis characteristics of the NSD, but also prove the great potential of the proposed damper in vibration mitigation of drilling equipment to platform structures.

Data Availability
e data used to support the finding of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.