Theoretical Design of a Novel Vibration Energy Absorbing Mechanism for Cables

: A novel design of a vibration energy absorbing mechanism (VEAM) that is based on multi-physics (magnetic spring, hydraulic system, structural dynamics, etc.) for cable vibration is proposed. The minimum working force of the hydraulic cylinder has been exploited in this design in order to combine a non-linear sti ﬀ ness vibration isolation module that is composed of permanent magnetic springs with hydraulic viscous vibration damping modules. In response to di ﬀ erent environmental vibration impacts, VEAM can automatically switch the vibration control modes without an electronic mechanism. Additionally, the non-contact design e ﬀ ectively reduces the wear that is induced by the reciprocating motion of the small amplitude of the hydraulic viscous dampers. The proposed mechanism is explained and a theoretical model is established. The transmissibility of the two modules at a single degree of freedom is derived using the harmonic balance method. After that, a series of variable control numerical simulations were performed for each important parameter. Empirical rules for designing the system were created by comparing the inﬂuence of each parameter on the vibration isolation performance of the entire system. peak transmissibility value. When the frequency is 13.0 Hz, the maximum transmissibility of the sub-system with a spring structural damping ratio of 0.06 is 65.2% lower than that with a structural damping ratio of 0.02. This result shows that structural damping has a considerable e ﬀ ect on the overall vibration damping characteristic of the system. S.-K.L.; data curation, Z.Q., Y.-T.W. and A.H.; writing—original draft preparation, Z.Q. and Y.-T.W.; writing—review and editing, A.H., S.-K.L. and J.W.S.; visualization, Z.Q.; supervision, Y.-T.W.; project administration, S.-K.L.; funding acquisition, S.-K.L. All authors have the published of


Introduction
Cables are widely used in construction, cable cars, traction, and other fields due to their extremely high tensile strength, long-distance capabilities, and low cost. This article will conduct research that is based on the cable-stayed bridge, which is the most representative and technically demanding application of cables, and illustrate a novel vibration energy absorbing mechanism that effectively responds to common vibration that is caused by environmental impacts with the automatic switching modes of the non-stiffness vibration isolators with the hydraulic viscous dampers.
In recent years, the cable-stayed bridge design has gained popularity as the design of choice for the construction of long-span bridges. This can be attributed to its simple and beautiful design, and its ability to span longer distances at relatively low construction costs [1]. From 2016 to 2019, a total of 20 cable stayed bridges with spans over 500 m were built around the world. Currently, the longest span cable-stayed bridge in the world is the Russky Bridge built over the Eastern Bosphorous in Russia, which has a span of 1104 m. Although the span of this bridge is quite short when compared to suspension bridges, the trend in construction seems to be moving towards an increasing number and require external power supply and can handle small amplitude vibrations with minimal wear while also being able to cater to larger amplitude vibrations.
Consequently, this paper proposes a hybrid damping system design that combines a PM based non-linear stiffness vibration isolator and a hydraulic viscous damper. In the proposed design, when the amplitude of vibration in the cable is less than a preset threshold value, the non-linear stiffness isolator damps the vibration transmitted to the deck and the viscous damper is not actuated. This will help to increase the overall system life, as it protects the viscous damper from thermal deformation and excessive wear of the seal ring due to its high frequency reciprocating motion that is associated with low amplitude vibrations. When the vibration amplitude becomes greater than the threshold value, the PMs travel until there is a minimum relative distance between them and then the viscous dampers are actuated to further dissipate the vibration energy. Thus, the presented system, which we have called a vibration energy absorbing mechanism (VEAM), is developed to meet three objectives: (1) reduce the wear due to inefficient working conditions of the viscous dampers used in the system.
(2) Provide a non-contact link between the viscous damper and the cable so that shear damage to the damper due to multi-directional cable motion can be avoided. (3) Achieve damping of a wide range of vibration amplitudes without using complex electronic components, such as sensors and controllers. The detailed design of the VEAM and validation through simulations are presented in this paper. Figure 1 presents the 3D model of VEAM. The cable passing through the center is supported on six sides by the developed multi-physic dampers, which are in turn connected to the bridge deck. Traditional viscous damper systems usually adopt a three damper arrangement [16]. In contrast, the VEAM adopts an arrangement of three opposing pairs of dampers resulting in a six-directional cable support system. This has been done to better counter the randomness of stay cable vibration. Figure 2a shows the internal configuration of the VEAM. Because the system is symmetrically constructed using six identical modules, we will only use one of those modules to illustrate the design principles. The complete system is fitted to the cable via a central hex piece that is divided into two parts which are assembled using long through bolts. The six damper assemblies are attached to this hex piece, one on each side. The damper assemblies consist mainly of a ball-joint connected to the hex piece, a PM based non-linear stiffness isolator, a viscous damper, and a rotary mount to allow the rotation of the damper assembly about an axis perpendicular to the damper's axis of travel. The vibration of the stay cable moves the magnets, which dampen it due to non-linear stiffness. The vibration force is transmitted to the viscous damper through repulsion between the PMs placed with like poles facing each other. If the force magnitude is greater than the minimum force required to actuate the viscous damper, then the damper starts moving and dissipates the vibration energy through viscous drag.

Mechanism of the VEAM
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 24 the amplitude of vibration in the cable is less than a preset threshold value, the non-linear stiffness isolator damps the vibration transmitted to the deck and the viscous damper is not actuated. This will help to increase the overall system life, as it protects the viscous damper from thermal deformation and excessive wear of the seal ring due to its high frequency reciprocating motion that is associated with low amplitude vibrations. When the vibration amplitude becomes greater than the threshold value, the PMs travel until there is a minimum relative distance between them and then the viscous dampers are actuated to further dissipate the vibration energy. Thus, the presented system, which we have called a vibration energy absorbing mechanism (VEAM), is developed to meet three objectives: (1) reduce the wear due to inefficient working conditions of the viscous dampers used in the system. (2) Provide a non-contact link between the viscous damper and the cable so that shear damage to the damper due to multi-directional cable motion can be avoided. (3) Achieve damping of a wide range of vibration amplitudes without using complex electronic components, such as sensors and controllers. The detailed design of the VEAM and validation through simulations are presented in this paper. Figure 1 presents the 3D model of VEAM. The cable passing through the center is supported on six sides by the developed multi-physic dampers, which are in turn connected to the bridge deck. Traditional viscous damper systems usually adopt a three damper arrangement [16]. In contrast, the VEAM adopts an arrangement of three opposing pairs of dampers resulting in a six-directional cable support system. This has been done to better counter the randomness of stay cable vibration. Figure  2a shows the internal configuration of the VEAM. Because the system is symmetrically constructed using six identical modules, we will only use one of those modules to illustrate the design principles. The complete system is fitted to the cable via a central hex piece that is divided into two parts which are assembled using long through bolts. The six damper assemblies are attached to this hex piece, one on each side. The damper assemblies consist mainly of a ball-joint connected to the hex piece, a PM based non-linear stiffness isolator, a viscous damper, and a rotary mount to allow the rotation of the damper assembly about an axis perpendicular to the damper's axis of travel. The vibration of the stay cable moves the magnets, which dampen it due to non-linear stiffness. The vibration force is transmitted to the viscous damper through repulsion between the PMs placed with like poles facing each other. If the force magnitude is greater than the minimum force required to actuate the viscous damper, then the damper starts moving and dissipates the vibration energy through viscous drag.  and damage to the joint. At the other end of the damper assembly, the cylinder of the viscous damper is connected to the system frame using a ball bearing based rotary joint shown in Figure 2c that allows for rotation of the damper assembly around an axis perpendicular to the cylinder axis. The addition of these two joints ensures that the two PMs are always coaxial and only axial loads are applied to the viscous damper. Furthermore, to provide structural damping while keeping the PMs in parallel [49,50], the design also includes a spring connecting the hex side PM and the viscous damper cylinder.  It can be seen from Figure 3a that when the hex moves due to cable movement, the six damper assemblies rotate to different angles to compensate for this displacement. However, the PMs and viscous damper always remain aligned with each other. The rotation angle of each damper assembly can be derived while using the geometric relationship between the key elements. Figure 3b shows The thrust force applied to the piston should be aligned with the cylinder axis and should have no radial component in order to keep the viscous damper moving smoothly with minimal amount of wear. In the proposed system, the source of this thrust force is the PM based isolator attached to the end of the hydraulic piston. Thus, the necessary conditions for the smooth operation of the hydraulic cylinder can only be guaranteed if the force applied to the PM assembly has no radial component. However, in reality, the direction of vibration is not always aligned with the center axis of the damper assembly. Therefore, a compliant attachment between the damper assembly and the cable is required to ensure that only axial forces are transmitted to the damper assembly. In the proposed design, a ball joint is added between the hex and the PM side of the damper assembly in order to allow for the cable fixing hex to change angles with respect to the damper assembly. This joint, as shown in Figure 2b, consists of a ball at the center with four retaining springs around it to prevent excessive movement and damage to the joint. At the other end of the damper assembly, the cylinder of the viscous damper is connected to the system frame using a ball bearing based rotary joint shown in Figure 2c that allows for rotation of the damper assembly around an axis perpendicular to the cylinder axis. The addition of these two joints ensures that the two PMs are always coaxial and only axial loads are applied to the viscous damper. Furthermore, to provide structural damping while keeping the PMs in parallel [49,50], the design also includes a spring connecting the hex side PM and the viscous damper cylinder.

Mechanism of the VEAM
It can be seen from Figure 3a that when the hex moves due to cable movement, the six damper assemblies rotate to different angles to compensate for this displacement. However, the PMs and viscous damper always remain aligned with each other. The rotation angle of each damper assembly can be derived while using the geometric relationship between the key elements. Figure 3b shows the simplified schematic diagram of this scenario. In this figure, points A, B, C, D, E, and F represent the rotating centers of the bearings between the hydraulic cylinders and the system structure, while d and ∠θ represent the magnitude and angle (acute angle with respect to the horizontal) of displacement of the center point of the cable fixing hex. Circles with black solid hatching represent the joint balls. G is a ball in the initial (default) position and G' is the final position of ball G after the displacement.  It can be seen from Figure 3a that when the hex moves due to cable movement, the six damper assemblies rotate to different angles to compensate for this displacement. However, the PMs and viscous damper always remain aligned with each other. The rotation angle of each damper assembly can be derived while using the geometric relationship between the key elements. Figure 3b shows the simplified schematic diagram of this scenario. In this figure, points A, B, C, D, E, and F represent the rotating centers of the bearings between the hydraulic cylinders and the system structure, while d and ∠θ represent the magnitude and angle (acute angle with respect to the horizontal) of displacement of the center point of the cable fixing hex. Circles with black solid hatching represent the joint balls. G is a ball in the initial (default) position and G' is the final position of ball G after the displacement.
(a) Comparison of system positions before and after displacement  It can be observed that the displacement GG′ is equal to d and the angle ∠G′GA is equal to ∠θ. The distances from the cable fixing hex center to the ball and to the rotating center of the bearing at the initial position are known. Thus, the distance between the ball at position G and the bearing rotating center at the initial position is GA = OA − OG.
According to the law of cosines, G′A can be deduced as: Here, we assume that the radial deformation of the spring is zero, which means that it can only perform telescopic movement along the axial direction. Through the law of sines, the rotation angle of the hydraulic cylinder can be expressed as: From the geometric relationship, it can be seen that the rotation angle of the cable fixing hex side magnet around the ball center and the rotation angle of the hydraulic cylinder are alternate interior angles. Thus, combining Equations (1) and (2), the ∠A can be rewritten as: The rotation angles of the other hydraulic cylinders can also be calculated with Equations (4)-  It can be observed that the displacement GG is equal to d and the angle ∠G GA is equal to ∠θ. The distances from the cable fixing hex center to the ball and to the rotating center of the bearing at the initial position are known. Thus, the distance between the ball at position G and the bearing rotating center at the initial position is GA = OA − OG.
According to the law of cosines, G A can be deduced as: Here, we assume that the radial deformation of the spring is zero, which means that it can only perform telescopic movement along the axial direction. Through the law of sines, the rotation angle of the hydraulic cylinder can be expressed as: From the geometric relationship, it can be seen that the rotation angle of the cable fixing hex side magnet around the ball center and the rotation angle of the hydraulic cylinder are alternate interior angles. Thus, combining Equations (1) and (2), the ∠A can be rewritten as: The rotation angles of the other hydraulic cylinders can also be calculated with Equations (4)-(8): Appl. Sci. 2020, 10, 5309 A conventional viscous damper forms a hydraulic circuit by connecting the two hydraulic chambers cut off by the piston by opening an orifice in the piston itself. In this case, mainly the friction between the hydraulic oil and the orifice provides the damping effect. However, because the orifice is inside the hydraulic cylinder, it is difficult to change the orifice size after assembly. In order to facilitate the adjustment of the damping ratio, this proposed design utilizes an external circuit design, as shown in Figure 4. The double triangle marks in the figure indicate the orifice regulating valves, while LL, LR, RL, and RR represent the hydraulic chambers divided by the pistons. When the cable fixing hex is in the default position, the left and right subsystems of the VEAM are in a balanced state and the pistons on both sides are positioned in the middle of their respective hydraulic cylinders with equal pressures on both sides of the pistons. The LL chamber is connected to the RR chamber and the LR chamber is connected to the RL chamber. The purpose of this inter-connection is to ensure the continuity of the damping process. Assuming that the repulsive force between the PMs is large enough, when the cable fixing hex moves leftwards the magnet pushes the left side piston, reducing the LL chamber volume and increasing its internal pressure. Therefore, the hydraulic oil in the LL chamber flows through the orifice and the piping to the RR chamber, which is at a relatively lower pressure. At the same time, the pressure in the LR chamber is reduced, which sucks in the hydraulic oil from the RL chamber. This creates a pressure difference between the RL and RR chambers that causes the right side piston to move out. This mechanism ensures that when the cable fixing hex moves, the pistons on both sides simultaneously move in the same direction. In this way, when the hex moves to the maximum amplitude in one direction and then moves in the opposite direction, the two PMs on the opposite side will be close enough to generate a repulsive force form the movement onset, thus preventing the occurrence of sudden cable movement or the jumping phenomenon. Because the damper assemblies are designed to move, the hydraulic lines used in the prototype have been coiled to reduce the stress at coupling points. Conventional viscous dampers usually use internal or external springs to restore the piston to its default position. Because of the design of the cable-stayed bridge, the cables are always in tension and so always tend to return to their default position. The VEAM design exploits this characteristic and uses the cable movement to automatically restore the hydraulic dampers to their default positions, thus doing away with the need for an additional spring.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 24 A conventional viscous damper forms a hydraulic circuit by connecting the two hydraulic chambers cut off by the piston by opening an orifice in the piston itself. In this case, mainly the friction between the hydraulic oil and the orifice provides the damping effect. However, because the orifice is inside the hydraulic cylinder, it is difficult to change the orifice size after assembly. In order to facilitate the adjustment of the damping ratio, this proposed design utilizes an external circuit design, as shown in Figure 4. The double triangle marks in the figure indicate the orifice regulating valves, while LL, LR, RL, and RR represent the hydraulic chambers divided by the pistons. When the cable fixing hex is in the default position, the left and right subsystems of the VEAM are in a balanced state and the pistons on both sides are positioned in the middle of their respective hydraulic cylinders with equal pressures on both sides of the pistons. The LL chamber is connected to the RR chamber and the LR chamber is connected to the RL chamber. The purpose of this inter-connection is to ensure the continuity of the damping process. Assuming that the repulsive force between the PMs is large enough, when the cable fixing hex moves leftwards the magnet pushes the left side piston, reducing the LL chamber volume and increasing its internal pressure. Therefore, the hydraulic oil in the LL chamber flows through the orifice and the piping to the RR chamber, which is at a relatively lower pressure. At the same time, the pressure in the LR chamber is reduced, which sucks in the hydraulic oil from the RL chamber. This creates a pressure difference between the RL and RR chambers that causes the right side piston to move out. This mechanism ensures that when the cable fixing hex moves, the pistons on both sides simultaneously move in the same direction. In this way, when the hex moves to the maximum amplitude in one direction and then moves in the opposite direction, the two PMs on the opposite side will be close enough to generate a repulsive force form the movement onset, thus preventing the occurrence of sudden cable movement or the jumping phenomenon. Because the damper assemblies are designed to move, the hydraulic lines used in the prototype have been coiled to reduce the stress at coupling points. Conventional viscous dampers usually use internal or external springs to restore the piston to its default position. Because of the design of the cable-stayed bridge, the cables are always in tension and so always tend to return to their default position. The VEAM design exploits this characteristic and uses the cable movement to automatically restore the hydraulic dampers to their default positions, thus doing away with the need for an additional spring.

Theoretical Modeling of the VEAM
The VEAM is composed of six identical damper assemblies, where each pair of opposing damper assemblies is connected through a hydraulic circuit to form a sub-system. Thus, each sub-system is composed of two magnetic non-linear stiffness modules and two hydraulic viscous damping modules. The isolation and damping performance of the VEAM can be evaluated by calculating its transmissibility [51]. The harmonic balance method (HBM) is an efficient method for modeling periodic problems and it is widely used to derive transmissibility [41]. Therefore, in this paper, HBM is used to establish theoretical models for the non-linear stiffness and the viscous damping modules. In this section, a mathematical model for the mechanism of a subsystem of VEAM in a single degree of freedom is established.

Automatic Switching Between Non-Linear Stiffness and Viscous Damping Modules
Figure 5a presents a simplified single degree of freedom theoretical model of a subsystem of the VEAM along the direction of its central axis. The thick line in the middle indicates the cable, with its two ends fixed between the deck and the pylon under a tensile force (T) acting at each end. The total length of the cable is L. The VEAM is fixed at a distance L1 from the deck. The cable fixing hex is connected to the cable and the damping system structure is fixed to the deck. In the figure, y c represents the excitation of the stay cable and y d represents the absolute displacement of the deck. As we are only considering the single degree of freedom case, the role of the ball joint and the rotary joint of the hydraulic cylinder is not considered. Hence, VEAM is simplified as a model in which the magnetic force that is generated by the PMs is connected in series with the hydraulic cylinder and then in parallel with the spring.
The minimum working pressure of the hydraulic cylinder is mainly caused by the friction of the piston seal, which is mostly dependent on the manufacturing and assembly accuracies, etc. In recent years, in order to improve the efficiency of hydraulic cylinders, researchers have tried to minimize this minimum working pressure through various methods, but static friction that is caused by seals is inevitable. With the development of industrial technology, it is not difficult to maintain the minimum working pressure at a fixed value by controlling the precision of machining and seal manufacturing. This minimum working pressure can be easily determined through experiments. The VEAM realizes automatic switching between the non-linear vibration isolator and the hydraulic viscous damper based on the minimum working pressure of the hydraulic cylinder.

Theoretical Modeling of the VEAM
The VEAM is composed of six identical damper assemblies, where each pair of opposing damper assemblies is connected through a hydraulic circuit to form a sub-system. Thus, each sub-system is composed of two magnetic non-linear stiffness modules and two hydraulic viscous damping modules. The isolation and damping performance of the VEAM can be evaluated by calculating its transmissibility [51]. The harmonic balance method (HBM) is an efficient method for modeling periodic problems and it is widely used to derive transmissibility [41]. Therefore, in this paper, HBM is used to establish theoretical models for the non-linear stiffness and the viscous damping modules. In this section, a mathematical model for the mechanism of a subsystem of VEAM in a single degree of freedom is established.

Automatic Switching Between Non-Linear Stiffness and Viscous Damping Modules
Figure 5a presents a simplified single degree of freedom theoretical model of a subsystem of the VEAM along the direction of its central axis. The thick line in the middle indicates the cable, with its two ends fixed between the deck and the pylon under a tensile force (T) acting at each end. The total length of the cable is L. The VEAM is fixed at a distance L1 from the deck. The cable fixing hex is connected to the cable and the damping system structure is fixed to the deck. In the figure, y represents the excitation of the stay cable and y represents the absolute displacement of the deck.
As we are only considering the single degree of freedom case, the role of the ball joint and the rotary joint of the hydraulic cylinder is not considered. Hence, VEAM is simplified as a model in which the magnetic force that is generated by the PMs is connected in series with the hydraulic cylinder and then in parallel with the spring.
(a) Schematic of the sub-system mechanism  The minimum working pressure of the hydraulic cylinder is mainly caused by the friction of the piston seal, which is mostly dependent on the manufacturing and assembly accuracies, etc. In recent years, in order to improve the efficiency of hydraulic cylinders, researchers have tried to minimize this minimum working pressure through various methods, but static friction that is caused by seals is inevitable. With the development of industrial technology, it is not difficult to maintain the minimum working pressure at a fixed value by controlling the precision of machining and seal manufacturing. This minimum working pressure can be easily determined through experiments. The VEAM realizes automatic switching between the non-linear vibration isolator and the hydraulic viscous damper based on the minimum working pressure of the hydraulic cylinder.
For each subsystem, it can be observed from Figures 2 and 4 that one PM is fixed on the piston rod and the other is connected to the vibration excitation side, such that like the poles of the two PMs are opposite to each other. When a large amplitude vibration occurs at the excitation side, the distance between the PMs becomes small and the repulsive force increases. This repulsive force is transmitted to the piston rod as a thrust force (F ). When this thrust force is less than the minimum working force (F ) of the hydraulic cylinder only the non-linear vibration isolator works. Conversely, only the For each subsystem, it can be observed from Figures 2 and 4 that one PM is fixed on the piston rod and the other is connected to the vibration excitation side, such that like the poles of the two PMs are opposite to each other. When a large amplitude vibration occurs at the excitation side, the distance between the PMs becomes small and the repulsive force increases. This repulsive force is transmitted to the piston rod as a thrust force (F m ). When this thrust force is less than the minimum working force (F f ) of the hydraulic cylinder only the non-linear vibration isolator works. Conversely, only the hydraulic viscous damper works when the thrust force is greater than the minimum working force of the hydraulic cylinder. Yamakawa et al. [52] proposed that the repulsive force between two PMs is inversely proportional to the nth power of the distance between them. The two magnets that are discussed in this article are identical and, therefore, satisfy the conditions for the application of this relationship. Thus, the repulsive force between the two PMs can be taken as, where, a (Unit: Nm −n ) is the equivalence coefficient of the repulsive force, n denotes the characteristic index, l (Unit: m) represents the preset distance between the PMs, and y (Unit: m) denotes the relative displacement between the cable and ground. When only the magnetic non-linear stiffness module works, the inequality F m ≤ F f exists. Substituting Equation (9) into this inequality, we get the relationship between the preset distance of the two PMs and the excitation amplitude:

Nonlinear Stiffness Isolator Model
In the VEAM sub-system, there is a set of mutually exclusive PMs at each end of the cable fixing hex, these two sets of PMs form a magnetic spring. The force exerted on the cable fixing hex by the repulsion of the magnetic spring can be expressed as the following Duffing type cubic polynomial equation [40,41].
F ms = k b y + by 3 (11) where, F ms represents the resultant force provided by the magnetic spring and k b and b represent the equivalent linear and non-linear stiffness coefficients, respectively. When F m ≤ F f , the hydraulic cylinder does not work, and the PM at the end of the hydraulic piston rod can be assumed as grounded. The theoretical model of the non-linear stiffness isolator module of the VEAM can be established as a system composed of; (1) the mass (m), (2) the linear stiffness of the coil spring (k S1 , k S2 ), (3) the structural damping from the two coil springs on each side of the sub-system (which is simplified as Rayleigh damping in this study) (c 1 , c 2 ), and (4) the equivalent linear and non-linear stiffness due to the magnetic spring (k b ), as shown in Figure 5b. The linear stiffness coefficient of the VEAM subsystem can be integrated as k l = k S1 + k S2 + k b .
To simplify the calculations, the vibration excitation from the stay cable is calculated as the following sine function: y c = Y c sin(ωt + ϕ) where, Y c is the vibration amplitude, is the angular frequency, t is the time, and ϕ is the phase shift. Thus, the displacement of the VEAM subsystem (y), which is the relative displacement between the stay cable and the deck, can be expressed as: When the VEAM sub-system receives the vibration transmitted from the stay cable and enters a steady state, the equation of motion of the sub-system according to the HBM can be written as the following Duffing type equation [53,54]: Because the same springs are used on both sides of the sub-system, we assume that their structural damping coefficients are equal, i.e., c = c 1 = c 2 .
Substituting Equations (12) and (13) into Equation (14) and deriving the time, the following equation is obtained.
− mYω 2 sin(ωt) + 2cYω cos(ωt) + k l Ysin(ωt) + bY 3 sin 3 (ωt) = mY c ω 2 [sin(ωt)cos ϕ + cos(ωt)sin ϕ] The cubic term in this equation can be broken down as: sin 3 (ωt) = [3 sin(ωt) − sin(3ωt)]/4. The non-harmonic high frequency term that, in general, does not simply decay with time can be omitted and the constant terms with sin(ωt) and cos(ωt) can be integrated as: and 2cYω = mY c ω 2 sin ϕ The absolute displacement of the deck is the sum of the vibration excitation of the stay cable and the relative displacement between the stay cable and the deck: Therefore, the displacement transmissibility of the VEAM sub-system can be sorted as: Substituting Equation (16) into Equation (19), we get:

Hydraulic Viscous Damper Model
When F f is less than F m , the PM pushes the piston and the hydraulic viscous damper begins to work. The damping due to the hydraulic cylinder is much greater than that due to the structural damping of the coil springs. Therefore, the structural damping is neglected in the theoretical model. The VEAM sub-system can be simplified as an SDOF model, as shown in Figure 5c. The figure graphically shows that the damping of the hydraulic cylinder (c h ) is connected in parallel with the stiffness of the two coil springs (k s = k S1 + k S2 ). Thus, the equation of motion of the VEAM sub-system can be expressed as: m ..
The derivation process to get the displacement transmissibility is similar to that described in Section 3.2. Thus, the equation is obtained, as follows:

Numerical Simulation
VEAM is a multi-physics system composed of a non-linear permanent magnet vibration isolator module and a hydraulic viscous damping module. The vibration isolation or damping performance of the VEAM may be greatly influenced by the different design factors. It is necessary to explore the influence of each design factor on system performance in order to design the VEAM system, so that it meets the highest requirements for cable-stayed bridges. In this section, we use the control variates method to change each of the key parameter values to observe and discuss its impact on the overall system performance. Table 1 lists the default parameters used in the simulation. In each simulation, one of these values is changed, while the rest remain constant. The simulations were coded using Python 3.7 and processed on a workstation with a 3.6 GHz Intel Xeon processor and 16GB RAM. The non-linear equations were solved using the newton_krylov method in the scipyoptimize module [55]. Section 4.1 discusses the relationship between the required distance between PMs and the starting force of the hydraulic cylinder with respect to different magnetic factors. Section 4.2 shows the simulation and discussions addressing the characteristics of the non-linear stiffness magnetic isolation module when the parameters of the coil spring and magnet are varied. In Section 4.3, the characteristics of hydraulic viscous damper under different design parameters are reported.

Item Value
Stiffness of the springs, k s1 = k s2 Nm

Preset Distance of Magnet
In Section 3.1, we mentioned that the automatic switching of the non-linear stiffness vibration isolator and hydraulic viscous damper is mainly determined by the distance between the magnets (l) and the minimum working force (F f ) of the hydraulic cylinder. However, different permanent magnets exhibit different magnetic repulsive force characteristics. We must understand the effect of the hydraulic cylinder working force on the distance between PMs under different PM parameters in order to design a VEAM that can switch modes when the input vibration amplitude is above the specified value. In order to evaluate this, we set the mode switching amplitude y as 0.002 m in this simulation. Figure 6 presents the relationship between the designed distance between the PMs and the minimum working force of the hydraulic cylinder according to Equation (10). We performed simulations that are based on three kinds of PMs with different magnetic repulsive force characteristics. The solid line represents the curve obtained using the PMs with a = 0.4 and n = 1.0, where a and n are consistent with Equation (9). The dotted line and the dash-dot-dash line represent the curves for when a = 0.6, n = 1.1, and a = 0.8, n = 1.2, respectively. Because the magnet calculation formula involves the power function calculation that may get an infinite value, the ordinate is only taken to an effective distance of 0.05 m.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 24 stronger magnets in the future by possibly utilizing some sort of motion limiting mechanism or collision prevention buffers between the magnets. A Simulink model is established using the Simscape module available in Matlab 2019, which is a special numerical simulation software for solving multi-physical problems, in order to verify the automatic switching function between the nonlinear vibration isolator and the hydraulic viscous damper modules. Figure 7 shows the Simulink model of the entire VEAM subsystem. Table 1 shows It can be observed that the three curves show an exponential downward trend as the minimum working force required by the hydraulic cylinder increases, which also reflects the characteristics of the magnet repulsive force. We can also observe that, for the curve with a larger magnetic repulsion factor, the allowable distance between PMs is also greater. Assuming that the minimum working force of the hydraulic cylinder is 200 N, the curve obtained using the magnet repulsive force factor of a = 0.4 and n = 1.0 shows a distance of 0.004 m, which is very close to the preset amplitude value. The reciprocating motion that is caused by vibration is likely to cause the PMs to collide with each other, which may cause damage to the magnet or other components of the VEAM. During the design stage, the allowable distance between PMs should be greater than the preset amplitude to ensure the safety of the system. At the same minimum working force, the allowable distance on the dash-dot-dash curve with a strong magnetic repulsive force is 0.012 m, which is two times higher than the solid curve. Even so, when considering the limitations of geometry and safety factor in actual application, the allowable distance between PMs provided by the simple structural Neodymium magnets (Nd2Fe14B) cannot fully guarantee the safety of the system. However, this work only provides a basis for the design optimization of the VEAM. Based on these methods, we can devise a safer system using stronger magnets in the future by possibly utilizing some sort of motion limiting mechanism or collision prevention buffers between the magnets.
A Simulink model is established using the Simscape module available in Matlab 2019, which is a special numerical simulation software for solving multi-physical problems, in order to verify the automatic switching function between the nonlinear vibration isolator and the hydraulic viscous damper modules. Figure 7 shows the Simulink model of the entire VEAM subsystem. Table 1 shows the default parameter values used in this simulation. The nonlinear stiffness isolator module, hydraulic viscous damper module, and the entire VEAM subsystem are simulated separately in order to discuss whether the developed model is valid. In the simulation of the nonlinear stiffness isolator module, the input amplitude and frequency are set at 0.001 m and 1 Hz. In the simulations of hydraulic viscous damper module and the entire VEAM subsystem, the input amplitude and frequency are set at 0.005 m and 1 Hz, respectively. In order to reduce the calculation time, the symmetrically connected hydraulic cylinders are simplified as general viscous damping and the simulation is performed using ode23t solver. No errors or warnings were reported during the entire simulation processes, which shows that each mechanism can operate normally, and that the nonlinear vibration isolator and hydraulic viscous damper can be switched without any problems. Figure 8 shows the simulation results for the acceleration and the corresponding Fourier spectrum of each mechanism. As shown in Figure 8a, due to the nonlinear preloading force, the acceleration is relatively large at the beginning, and it tends to gradually reduce to a steady state. It can be seen from the corresponding Fourier spectrum (Figure 8b) that the acceleration presents a non-linear decreasing curve over the entire bandwidth. When the nonlinear stiffness isolation module is removed and only the hydraulic damping module is operating, the output becomes a familiar damped simple harmonic motion form, this simulation result is shown in Figure 8c,d. It can be seen from Figure 8e that the high-frequency harmonics appear at the beginning of the bandwidth and their period cannot be clearly distinguished when the entire VEAM is in a state where the nonlinear stiffness isolation module and the hydraulic damping module operate together. This may be attributed to the nonlinear dynamics of the magnetic force that is generated by the permanent magnets. Figure 8f       (e) Simulation results for the acceleration of the mechanism sub-system (f) The corresponding Fourier spectrum of the mechanism sub-system

Non-Linear Stiffness Isolator
When the magnetic repulsive force is less than the minimum working force of the hydraulic cylinder, only the non-linear stiffness vibration isolator composed of the magnetic springs works. In this section, we will individually alter the key parameters of the coil springs and PMs used in the non-linear vibration isolator in order to observe the influence of the parameters k , k , b, and ζ on the vibration isolation performance of the system and summarize the design experience.
In the first case, the linear stiffness of the coil springs is set as 1000 Nm −1 , 2000 Nm −1 , and 3000 Nm −1 , while the other parameters are kept at the default values that are shown in Table 1, and the displacement transmissibility of the non-linear stiffness isolator module is observed. The simulation results are shown in Figure 9. It can be observed that the simulated transmissibility curve moves to the high frequency region with the increase in linear stiffness of the spring. This is accompanied by an increase in the peak value of transmissibility. When the stiffness is 1000 Nm −1 , the maximum displacement transmissibility of 3.9 occurs at a frequency of 16.5 Hz. When the stiffness is 3000 Nm −1 , the maximum displacement transmissibility of 4.5 occurs at a frequency of 21.7 Hz. This shows that increasing the linear stiffness of the spring in the VEAM can effectively reduce the displacement

Non-Linear Stiffness Isolator
When the magnetic repulsive force is less than the minimum working force of the hydraulic cylinder, only the non-linear stiffness vibration isolator composed of the magnetic springs works. In this section, we will individually alter the key parameters of the coil springs and PMs used in the non-linear vibration isolator in order to observe the influence of the parameters k s , k b , b, and ζ on the vibration isolation performance of the system and summarize the design experience.
In the first case, the linear stiffness of the coil springs is set as 1000 Nm −1 , 2000 Nm −1 , and 3000 Nm −1 , while the other parameters are kept at the default values that are shown in Table 1, and the displacement transmissibility of the non-linear stiffness isolator module is observed. The simulation results are shown in Figure 9. It can be observed that the simulated transmissibility curve moves to the high frequency region with the increase in linear stiffness of the spring. This is accompanied by an increase in the peak value of transmissibility. When the stiffness is 1000 Nm −1 , the maximum displacement transmissibility of 3.9 occurs at a frequency of 16.5 Hz. When the stiffness is 3000 Nm −1 , the maximum displacement transmissibility of 4.5 occurs at a frequency of 21.7 Hz. This shows that increasing the linear stiffness of the spring in the VEAM can effectively reduce the displacement transmissibility in the low frequency region. This means that the vibration isolation performance of the VEAM at a tuning frequency can be optimized using springs with higher stiffness values.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 24 transmissibility in the low frequency region. This means that the vibration isolation performance of the VEAM at a tuning frequency can be optimized using springs with higher stiffness values. The results of the second case are presented in Figure 10, which shows the effect of different equivalent linear stiffness obtained using different magnetic springs on the displacement transmissibility of the VEAM. The solid, dotted, and dash-dot-dash curves represent the magnetic spring equivalent linear stiffness values of 2000 Nm −1 , 3500 Nm −1 , and 5000 Nm −1 , respectively. It can be seen that the shifting of the transmissibility curve is similar to that observed while varying the linear stiffness of the spring; i.e., as the equivalent linear stiffness increases the curve shifts to the high frequency region. At frequencies below 18.5 Hz, the VEAM sub-system with a k value of 5000 Nm −1 has a lower transmissibility than that of the sub-system with a k value of 2000 Nm −1 . At 18.5 Hz, the transmissibility in both cases is equal with a value of 3.2. Beyond this frequency, the transmissibility with a k value of 5000 Nm −1 becomes higher than that with a k value of 2000 Nm −1 . Therefore, according to these results the VEAM sub-system with larger k value will have better vibration isolation performance in the low frequency range, but it will not work satisfactorily in the high frequency range  The results of the second case are presented in Figure 10, which shows the effect of different equivalent linear stiffness obtained using different magnetic springs on the displacement transmissibility of the VEAM. The solid, dotted, and dash-dot-dash curves represent the magnetic spring equivalent linear stiffness values of 2000 Nm −1 , 3500 Nm −1 , and 5000 Nm −1 , respectively. It can be seen that the shifting of the transmissibility curve is similar to that observed while varying the linear stiffness of the spring; i.e., as the equivalent linear stiffness increases the curve shifts to the high frequency region. At frequencies below 18.5 Hz, the VEAM sub-system with a k b value of 5000 Nm −1 has a lower transmissibility than that of the sub-system with a k b value of 2000 Nm −1 . At 18.5 Hz, the transmissibility in both cases is equal with a value of 3.2. Beyond this frequency, the transmissibility with a k b value of 5000 Nm −1 becomes higher than that with a k b value of 2000 Nm −1 . Therefore, according to these results the VEAM sub-system with larger k b value will have better vibration isolation performance in the low frequency range, but it will not work satisfactorily in the high frequency range.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 24 transmissibility in the low frequency region. This means that the vibration isolation performance of the VEAM at a tuning frequency can be optimized using springs with higher stiffness values. The results of the second case are presented in Figure 10, which shows the effect of different equivalent linear stiffness obtained using different magnetic springs on the displacement transmissibility of the VEAM. The solid, dotted, and dash-dot-dash curves represent the magnetic spring equivalent linear stiffness values of 2000 Nm −1 , 3500 Nm −1 , and 5000 Nm −1 , respectively. It can be seen that the shifting of the transmissibility curve is similar to that observed while varying the linear stiffness of the spring; i.e., as the equivalent linear stiffness increases the curve shifts to the high frequency region. At frequencies below 18.5 Hz, the VEAM sub-system with a k value of 5000 Nm −1 has a lower transmissibility than that of the sub-system with a k value of 2000 Nm −1 . At 18.5 Hz, the transmissibility in both cases is equal with a value of 3.2. Beyond this frequency, the transmissibility with a k value of 5000 Nm −1 becomes higher than that with a k value of 2000 Nm −1 . Therefore, according to these results the VEAM sub-system with larger k value will have better vibration isolation performance in the low frequency range, but it will not work satisfactorily in the high frequency range  b is the only non-linear stiffness parameter in the VEAM system, so it is necessary to consider its influence on the vibration isolation performance of the whole sub-system. The solid, dotted, and dash-dot-dash curves in Figure 11 indicate the simulated transmissibility of the VEAM sub-system with respect to b values of 0 (no non-linear stiffness component in the system), default value (8.1 × 10 6 ), and 10 × 10 9 (an ideal value chosen in order to clearly see the effect of increasing the value), respectively. It can be observed that, as the value of b increases, the peak value of the transmissibility decreases and shifts towards the high-frequency region. However, the curve appears to bend to the right and, as the value of b gets larger, the amount of bending also becomes greater (as shown in the dash-do-dash curve). A similar sort of transmissibility curve bending towards the high frequency region has also been reported in other works [40,41]. We speculate that the excessive curve bending caused by the large b value will cause multiple transmissibility values to exist in a part of the frequency range, which means that there will be a multiple states phenomenon. This may cause the dynamics of the system to become complicated with the change in amplitude and frequency of excitation and chaotic motion may occur, as shown in Figure 8b.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 19 of 24 b is the only non-linear stiffness parameter in the VEAM system, so it is necessary to consider its influence on the vibration isolation performance of the whole sub-system. The solid, dotted, and dash-dot-dash curves in Figure 11 indicate the simulated transmissibility of the VEAM sub-system with respect to b values of 0 (no non-linear stiffness component in the system), default value (8.1 × 10 6 ), and 10 × 10 9 (an ideal value chosen in order to clearly see the effect of increasing the value), respectively. It can be observed that, as the value of b increases, the peak value of the transmissibility decreases and shifts towards the high-frequency region. However, the curve appears to bend to the right and, as the value of b gets larger, the amount of bending also becomes greater (as shown in the dash-do-dash curve). A similar sort of transmissibility curve bending towards the high frequency region has also been reported in other works [40,41]. We speculate that the excessive curve bending caused by the large b value will cause multiple transmissibility values to exist in a part of the frequency range, which means that there will be a multiple states phenomenon. This may cause the dynamics of the system to become complicated with the change in amplitude and frequency of excitation and chaotic motion may occur, as shown in Figure 8b. In the non-linear stiffness isolator module consisting of a magnetic spring, there is no additional damping, except for the structural damping of the springs. We simulated the transmissibility of the VEAM sub-system with spring structural damping ratios of 0.02, 0.04, and 0.06. Figure 12 shows the results of these simulations. It can be observed that, the larger the structural damping ratio, the smaller the transmissibility of the VEAM sub-system. It is also apparent that the structural damping has almost no effect on the change of the frequency domain of the peak transmissibility value. When the frequency is 13.0 Hz, the maximum transmissibility of the sub-system with a spring structural damping ratio of 0.06 is 65.2% lower than that with a structural damping ratio of 0.02. This result shows that structural damping has a considerable effect on the overall vibration damping characteristic of the system. In the non-linear stiffness isolator module consisting of a magnetic spring, there is no additional damping, except for the structural damping of the springs. We simulated the transmissibility of the VEAM sub-system with spring structural damping ratios of 0.02, 0.04, and 0.06. Figure 12 shows the results of these simulations. It can be observed that, the larger the structural damping ratio, the smaller the transmissibility of the VEAM sub-system. It is also apparent that the structural damping has almost no effect on the change of the frequency domain of the peak transmissibility value. When the frequency is 13.0 Hz, the maximum transmissibility of the sub-system with a spring structural damping ratio of 0.06 is 65.2% lower than that with a structural damping ratio of 0.02. This result shows that structural damping has a considerable effect on the overall vibration damping characteristic of the system. Appl. Sci. 2020, 10, x FOR PEER REVIEW 20 of 24 Figure 12. Transmissibility of the non-linear stiffness isolator module of VEAM with varying values of ζ.

Hydraulic Damper
Under normal conditions, the stay cable vibration has a small amplitude. However, under special environmental impacts, such as earthquakes and typhoons, the amplitude of the stay cable vibration can become large. The VEAM incorporates hydraulic viscous dampers in its design to dampen such large amplitude vibrations by exploiting their high damping performance. This section explores the effects of damping and spring stiffness on the performance of this damper module. In this section, the excitation amplitude Y is set as 0.005 m.
In practical applications, the damping ratio of the hydraulic damper used in VEAM can be adjusted by changing the size of the orifice in the flow path. To make performance comparisons, we set the damping ratios of the hydraulic cylinders at 0.2, 0.6, and 0.8. Figure 13 shows the simulated transmissibility results of the VEAM sub-system. It can be seen that the higher the damping ratio of the hydraulic cylinder the lower the peak value of the transmissibility. However, although the curve with the lowest damping ratio of 0.2 has a higher transmissibility in the low-frequency region (region of amplification) within a characteristic frequency, the frequency region after this frequency, which is called the "isolation region", shows better vibration isolation performance.  Similar to Figure 9, the frequency domain of the peak wave of transmissibility curve becomes wider as the value of k increases, and the maximum transmissibility occurs at a higher frequency.

Hydraulic Damper
Under normal conditions, the stay cable vibration has a small amplitude. However, under special environmental impacts, such as earthquakes and typhoons, the amplitude of the stay cable vibration can become large. The VEAM incorporates hydraulic viscous dampers in its design to dampen such large amplitude vibrations by exploiting their high damping performance. This section explores the effects of damping and spring stiffness on the performance of this damper module. In this section, the excitation amplitude Y c is set as 0.005 m.
In practical applications, the damping ratio of the hydraulic damper used in VEAM can be adjusted by changing the size of the orifice in the flow path. To make performance comparisons, we set the damping ratios of the hydraulic cylinders at 0.2, 0.6, and 0.8. Figure 13 shows the simulated transmissibility results of the VEAM sub-system. It can be seen that the higher the damping ratio of the hydraulic cylinder the lower the peak value of the transmissibility. However, although the curve with the lowest damping ratio of 0.2 has a higher transmissibility in the low-frequency region (region of amplification) within a characteristic frequency, the frequency region after this frequency, which is called the "isolation region", shows better vibration isolation performance.

Hydraulic Damper
Under normal conditions, the stay cable vibration has a small amplitude. However, under special environmental impacts, such as earthquakes and typhoons, the amplitude of the stay cable vibration can become large. The VEAM incorporates hydraulic viscous dampers in its design to dampen such large amplitude vibrations by exploiting their high damping performance. This section explores the effects of damping and spring stiffness on the performance of this damper module. In this section, the excitation amplitude Y is set as 0.005 m.
In practical applications, the damping ratio of the hydraulic damper used in VEAM can be adjusted by changing the size of the orifice in the flow path. To make performance comparisons, we set the damping ratios of the hydraulic cylinders at 0.2, 0.6, and 0.8. Figure 13 shows the simulated transmissibility results of the VEAM sub-system. It can be seen that the higher the damping ratio of the hydraulic cylinder the lower the peak value of the transmissibility. However, although the curve with the lowest damping ratio of 0.2 has a higher transmissibility in the low-frequency region (region of amplification) within a characteristic frequency, the frequency region after this frequency, which is called the "isolation region", shows better vibration isolation performance.  Similar to Figure 9, the frequency domain of the peak wave of transmissibility curve becomes wider as the value of k increases, and the maximum transmissibility occurs at a higher frequency.   Figure 9, the frequency domain of the peak wave of transmissibility curve becomes wider as the value of k s increases, and the maximum transmissibility occurs at a higher frequency. However, the difference between the two figures is that the peak value of the transmissibility is not affected by the change in the value of k s . This is because, when the hydraulic shock absorber is in operation, the magnetic spring stops working, and the equivalent stiffness of the entire system no longer changes due to the nonlinear stiffness. Thus, for the damping module, the spring with a lower k s value shows better isolation performance in the high frequency region.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 24 However, the difference between the two figures is that the peak value of the transmissibility is not affected by the change in the value of k . This is because, when the hydraulic shock absorber is in operation, the magnetic spring stops working, and the equivalent stiffness of the entire system no longer changes due to the nonlinear stiffness. Thus, for the damping module, the spring with a lower k value shows better isolation performance in the high frequency region. Figure 14. Transmissibility of the hydraulic viscous damper module of VEAM with varying k .

Conclusions
This study has proposed a novel damper design for cable stayed bridges that can respond to different vibration amplitudes that are caused by different environmental phenomena. It is a hybrid system that combines non-linear stiffness vibration isolation modules and hydraulic viscous damper modules. The VEAM proposes a novel way to utilize the minimum working force of the hydraulic cylinders, which is an unavoidable characteristic that can be leveraged to combine the two modules together to form a hybrid system that is able to automatically switch the module in operation based on the input vibration amplitude. This paper presented the mechanism design and theoretical model of the VEAM. Multiple numerical case studies were performed in order to observe the influence of various factors on the vibration isolation and damping performance of the entire system. The key findings of this research are as follows: 1. VEAM mechanism, through the combined action of hydraulic cylinder bearings, ball joints, and springs, can maintain the alignment of PMs at all times, and will not cause vibration force transmission instability due to tangential shift. 2. The external hydraulic circuit connecting the opposing damper pairs can maintain motion synchrony to avoid the jumping phenomenon, thus preventing system damage. This method can also easily make adjustments of the orifice regulating valves to achieve a variable damping ratio. Through the Simscape simulation of the entire VEAM subsystem, the design theory of the automatic switching mechanism was tested. The results showed that, when the whole system is working, due to the non-linear stiffness, chaotic output is observed in the high frequency bandwidth. The implications of this behavior need to be verified in future experiments with the real system. 3. We can use strong repulsive magnets to ensure the maintenance of a safe distance between the PM pairs. The use of magnets with stronger linear stiffness will reduce the transmissibility of the nonlinear stiffness module in the low frequency region, but it will increase the transmissibility in the high frequency region. However, magnets with stronger nonlinear stiffness can effectively reduce the peak of transmissibility. In the design stage, according to the needs of the actual situation, the selection of PMs of appropriate magnetic strength can ensure a high performing safe design, while also ensuring a safe distance. However, during prototype manufacturing, we found that, due to space constraints, it is difficult to only use the magnetic

Conclusions
This study has proposed a novel damper design for cable stayed bridges that can respond to different vibration amplitudes that are caused by different environmental phenomena. It is a hybrid system that combines non-linear stiffness vibration isolation modules and hydraulic viscous damper modules. The VEAM proposes a novel way to utilize the minimum working force of the hydraulic cylinders, which is an unavoidable characteristic that can be leveraged to combine the two modules together to form a hybrid system that is able to automatically switch the module in operation based on the input vibration amplitude. This paper presented the mechanism design and theoretical model of the VEAM. Multiple numerical case studies were performed in order to observe the influence of various factors on the vibration isolation and damping performance of the entire system. The key findings of this research are as follows: 1.
VEAM mechanism, through the combined action of hydraulic cylinder bearings, ball joints, and springs, can maintain the alignment of PMs at all times, and will not cause vibration force transmission instability due to tangential shift.

2.
The external hydraulic circuit connecting the opposing damper pairs can maintain motion synchrony to avoid the jumping phenomenon, thus preventing system damage. This method can also easily make adjustments of the orifice regulating valves to achieve a variable damping ratio. Through the Simscape simulation of the entire VEAM subsystem, the design theory of the automatic switching mechanism was tested. The results showed that, when the whole system is working, due to the non-linear stiffness, chaotic output is observed in the high frequency bandwidth. The implications of this behavior need to be verified in future experiments with the real system.

3.
We can use strong repulsive magnets to ensure the maintenance of a safe distance between the PM pairs. The use of magnets with stronger linear stiffness will reduce the transmissibility of the nonlinear stiffness module in the low frequency region, but it will increase the transmissibility in the high frequency region. However, magnets with stronger nonlinear stiffness can effectively reduce the peak of transmissibility. In the design stage, according to the needs of the actual situation, the selection of PMs of appropriate magnetic strength can ensure a high performing safe design, while also ensuring a safe distance. However, during prototype manufacturing, we found that, due to space constraints, it is difficult to only use the magnetic force between the two PMs to actuate the hydraulic cylinder. Accordingly, it may be necessary to use crash pads between the magnets to achieve rigid transmission of force. The magnetic resistance of the surrounding environment and design involving stronger PMs will be considered in future studies. 4.
In the low-frequency region, the greater the spring linear stiffness, the lower will be the transmissibility of both the non-linear stiffness vibration isolation module and the hydraulic damper module, i.e., there will be better vibration isolation performance. However, in the high-frequency region, the stiffness of the spring will affect the vibration isolation performance when only the nonlinear stiffness vibration isolator is in operation. Therefore, the common effect of both modules must be considered when selecting the spring stiffness. 5.
The structural damping of the spring will undoubtedly provide good isolation performance when the non-linear stiffness module is working. However, the viscous damping of the hydraulic cylinder is different. Excessive viscous damping will affect the vibration isolation performance of the VEAM and increase the possibility of structural damage. On the other hand, small viscous damping has a large area of the isolation region, which can improve the poor vibration isolation that is caused by the shift towards the high frequency region of the transmissibility caused by excessive spring stiffness. Funding: This work was supported by the Regional Leading Research Center of NRF and MOCIE (NRF-2019R1A5A808320112).

Conflicts of Interest:
The authors declared that they have no conflict of interest to this work.