Modeling a Mechanical Molecular Spring Isolator with High-Static-Low-Dynamic-Stiffness Properties

A mechanical molecular spring isolator (MMSI) is proposed for the purpose of isolating the low-frequency vibration of a heavy payload. *e MMSI is a passive vibration isolation technique mimicking molecular spring isolator characteristics of high-staticlow-dynamic stiffness (HSLDS). AnMMSI consists of a piston-cylinder container filled with the liquid and some hydraulic spring accumulators. *e piston would support a lump of mass and be subjected to a specific external vibration excitation force. *ose accumulators can get intercommunication by the liquid to produce the transformation from high static stiffness to low dynamic stiffness. *e stiffness model of the MMSI with several identical accumulators is established based on the hydrostatic law. After that, some parameters that significantly influence the stiffness characteristics are studied. Results show that the stiffness property of this kind of MMSI demonstrates a piecewise linearity of three segments. It applies the averaging method to acquire amplitudefrequency and phase-frequency relationships of the piecewise linear vibration isolation system. An inevitable jump phenomenon may occur when the exciting force reaches the critical value. *e vibration isolation performance is evaluated by energy transmissibility. Finally, an experimental prototype was designed to carry out quasi-static and dynamic experiments to verify the stiffness model and the dynamic properties as an HSLDS vibration isolator.


Introduction
Vibration isolation has been increasingly demanded for high-speed vehicles, precision machinery, and megastructures as a result of the higher fatigue endurance for structures and higher human comfort expectations in both static and dynamic environments. A variety of vibration protection means are continuously improving to accommodate these new standards. Passive vibration control technology would always capture engineers' enthusiasm for its low costs and reliable vibration isolation performance [1][2][3]. Concerning an ideal case where a lump of mass m is supported by a linear spring with stiffness k on a rigid foundation, effective vibration isolation of this linear system can be achieved when the excitation frequency ω is greater than ���� � 2k/m √ . A lower stiffness of the supporting linear spring is required for isolating a lower frequency vibration. As a result, the linear spring will inevitably produce a higher static deformation, which may bring about more vertical installation space of the system and worse lateral stability of the linear spring. It causes the dilemma that low natural frequency and small static deformation cannot be achieved at the same time for a traditional linear isolator. Particularly, the dilemma is more remarkable on the heavy payload lowfrequency vibration isolation field. e nonlinear passive vibration isolators with the HSLDS property can easily break through this dilemma [4].
is kind of isolators only requires a smaller deformation while offering a lower stiffness when the payload is at the equilibrium position compared to the traditional isolator. In the recent decade, the study of this type of passive vibration isolators has been progressed. On the one hand, the most prevalent solution is to develop a negative stiffness device to combine with a positive stiffness element (such as the traditional linear spring). Due to the diversity of the negative stiffness device mechanisms, there exist plentiful nonlinear passive vibration isolators with the HSLDS property. Brennan's research group [5][6][7][8] developed two different types of nonlinear isolators. e first one consists of a vertical linear spring and two oblique springs; the second one is configured with coil springs and attracting magnets. Huang et al. [9] applied Euler-bulked beams as the negative stiffness element. Zhou et al. [10] utilized cam-roller-spring mechanisms as the negative stiffness element. Sun and Jing [11] found that scissor-like structures may also bring in the negative stiffness. Besides, Mofidian and Bardaweel [12] employed the magnetic spring as the positive stiffness element connecting with oblique springs to form a different isolator with the HSLDS property. ey all are able to achieve a quasi-zero stiffness vibration isolator by properly choosing system parameters. e aforementioned nonlinear isolators have been proved to be superior on isolation performance than traditional linear springs. However, on account of those mechanisms' geometrical and physical nonlinearities, some harmful nonlinear dynamic phenomena (the most common one is the sudden jump) would inevitably occur around the equilibrium position, and they get more severe as the magnitude of the excitation force increases. Besides, because of the connection of positive stiffness and negative stiffness elements, those isolators are normally with big size and complicated structures.
On the other hand, Chen's research group developed one type of solid and liquid mixture (SALiM) vibration isolators for isolating the low-frequency vibration of massive machinery.
e SALiM consists of the almost incompressible liquid and a number of compressible elastic solid elements. When under shock or vibration, the incompressible liquid can instantly pass the pressure on to all elastic solid elements, which causes all elastic solid elements to be compressed and deformed simultaneously. As a result, it could greatly absorb and dissipate the energy of vibrations or shock. e SALiM vibration isolator possesses a high loading capacity and a high stability of performance; it is especially aimed for the heavy payload low-frequency vibration isolation purpose. Previous work mainly focuses on the solid elastic element and so does the research of this article Teng and Chen [13] applied hollow rubber spheres as the elastic element and estimated the vibration isolation performance of the SALiM isolator. Gao et al. [14] took advantage of the U-shaped bellows as the elastic element to develop the B-SALiM isolator, and they denoted that the isolator can be modified as a nonlinear isolator with the HSLDS property. Inspired by the aforementioned work, Yu et al. [15,16] found that a kind of nanoporous hydrophobic nanofunctional materials, which can absorb and release some water until the hydraulic pressure reaches a certain value, can be applied to the solid elastic element. erefore, the molecular spring isolator based on nanoporous hydrophobic nanofunctional materials is developed, and it naturally presents the HSLDS property. Moreover, its stiffness exhibits nearly a long linearity around the static equilibrium position. Because of the characteristics of high density and high water absorption of the nanoporous hydrophobic nanofunctional materials, the molecular spring isolator owns a very small structure size.
Inspired by the distinct mechanism for the HSLDS property of the molecular spring isolator, we noted that the nanoporous hydrophobic nanofunctional materials only absorb and release water until the hydraulic pressure reaches a certain value. It seems that naturally there exists a switch in this kind of material related to the hydraulic pressure. A conventional hydraulic spring accumulator may acquire that function artificially. So we substitute the particular material of the molecular spring isolator with some hydraulic spring accumulators to form the mechanical molecular spring isolator (MMSI). Each hydraulic spring accumulator consists of a piston, a container, and a linear spring. e piston is sealed and is subjected to a force offered by the precompressed linear spring. As a result, the liquid cannot intrude into the accumulator until the hydraulic pressure reaches a certain value. e MMSI also possesses the HSLDS property. Besides, on account of the fact that the configuration of the accumulator is artificial, there is no necessity to utilize the specific water to get all accumulators intercommunicated. Any liquid is available; and the critical hydraulic pressure is adjustable, so the MMSI may not have to work at that high hydraulic pressure.
It is noticed that the configurations of all accumulators in the MMSI can be designed to be quite different from each other, which may result in a much more complicated stiffness property. In order to understand some simple and essential features of the MMSI, we assume that all accumulators are of the same configuration in this paper. e MMSI stiffness model is studied with the hydrostatic law based on that assumption, and it presents a piecewise linearity of three segments. e equation of motion of the MMSI isolation system is established at the center of the second segment. e averaging method is applied to investigate the primary resonance response of this kind of nonlinear vibration isolation system. e energy transmissibility is applied to evaluate the vibration isolation performance. An experimental prototype has been designed to carry out quasi-static and dynamic experiments to estimate the static and dynamic properties of the MMSI isolation system.

MMSI Stiffness Model.
In analogy with the function of the molecular spring, we mimic the molecular spring with hydraulic spring accumulators and the liquid. e combination of hydraulic spring accumulators and the liquid is referred to as mechanical molecular spring (MMS) correspondingly. e MMS is sealed in the piston-cylinder container.
e schematic diagram of an MMSI is presented in Figure 1(a). Define the MMSI piston diameter as D, the initial volume of the liquid as V 0 , and the initial hydraulic pressure as P 0 . Define n as the number of those accumulators, and there is at least one accumulator; that is, n ≥ 1. e structure of the accumulator is shown in Figure 1(b). e accumulator piston diameter is d, and it is sealed and matched with the container. Meanwhile, the piston is subjected to a force offered by the precompressed linear spring whose stiffness is k, and the precompressed force is F p . e movement of the piston is confined by two stoppers, which results in a maximum displacement δ that the piston can reach. Hence, there exist two critical hydraulic pressures of the accumulator. e first critical hydraulic pressure is formed by the precompressed force, and the definition of the first critical hydraulic pressure is Meanwhile the second critical hydraulic pressure is produced by the finite displacement of the accumulator piston, so the second critical hydraulic pressure is defined as Applying an external force F onto the MMSI container piston would cause the MMSI piston to have a vertical displacement x.
en the hydraulic pressure would become P, and the volume of the liquid would become V correspondingly. Among them, the external force is related to the hydraulic pressure as According to the hydrostatic law, the volume compressibility of a particular liquid is a constant in a specific environment. Its magnitude is decided by the hydraulic pressure and the liquid volume as follows: Since this paper focuses on studying the MMSI with all identical accumulators, there are only two critical hydraulic pressures. erefore, the hydraulic pressure starting from zero to infinity is divided into three segments. We define the segment P < P c1 as the bearing segment.
ere the hydraulic pressure is not able to overcome the precompressed force. So the displacement of the MMSI piston is only produced by the volume change of the liquid. en the vertical displacement x of the MMSI piston is given by Combining equation (5) with equations (3) and (4) derives the force-displacement relationship of the bearing segment as Define the segment P c1 ≤ P ≤ P c2 as the working segment. ere the hydraulic pressure can overcome the precompressed force, and each accumulator piston will produce an extra displacement Δl (0 ≤ Δl ≤ δ).
e displacement of the MMSI piston is produced by both the volume change of the liquid and the displacement of accumulator pistons. en the vertical displacement x of the MMSI piston is given by Combining equations (7) and (8) with equations (3) and (4) derives the force-displacement relationship of the working segment as Define the segment P c2 < P as the limit segment. ere the hydraulic pressure already pushes all accumulator pistons reaching the maximum displacement δ. e displacement of the MMSI piston is produced by both the volume change of the liquid and the maximum space all accumulator pistons can offer.
en the vertical displacement x of the MMSI piston is given by Combining equation (10) with equations (3) and (4) derives the force-displacement relationship of the limit segment as us, the force-displacement relationship of the MMSI can be written as Define the bearing stiffness or the limit stiffness as Define a geometrical parameter piston diameter ratio as Define an equivalent stiffness as Define the working stiffness as Define the length of the working segment as erefore, the force-displacement relationship of the MMSI can be rewritten as e MMSI presents a piecewise linear force-displacement relationship according to equation (18). e bearing stiffness is equal to the limit stiffness. e bearing stiffness or the limit stiffness refers to the stiffness of the liquid spring, which is formed by only the piston-cylinder container filled with the liquid and without any accumulator. An ordinarily high bearing stiffness or the limit stiffness is because of the low compressibility of the liquid, which ensures a relatively high load capacity. As a result, the payload only needs a small static displacement to reach the working segment. e working stiffness can be regarded as the series connection of an equivalent spring with the liquid spring according to equation (16). In particular, the equivalent spring stiffness is only in connection with the configuration of accumulators according to equation (15). Moreover, concerning (17), the length of the working segment is also quite related to the configuration of accumulators. e configuration of accumulators plays a quite critical fundamental function on the MMSI HSLDS property.

MMSI Dynamical Model.
It supposes that the payload locates exactly at the center of the working segment, and a linear viscous damper is added in parallel with the MMSI. Without loss of generality, the cosine wave is applied to express the harmonic excitation, so that the equation of motion of the MMSI vibration isolation system under harmonic excitation is given by where It is possible to obtain an exact quantitative solution by utilizing piecewise analytic integral as for a piecewise linear system. Meanwhile applying the averaging method is a simpler and convenient way to achieve the quantitative solution when the accuracy requirement is limited to the first-order term of the nonlinear factor [17].
Define the nonlinear factor as ε � (k 0 /k w ) − 1; and let , where σ is a detuning parameter. Define a piecewise function as follows: en the equation of motion can be rewritten in the dimensionless form as where In particular, when ε � 0, the solution is given by where a and θ are constants. However, when ε ≠ 0, the solution can still be represented by that solution, but a and θ are functions with respect to t rather than constants. By introducing ψ � ωt − θ and substituting ψ into the solution and then differentiating the solution with respect to t, we obtain

Combining equations (24) and (25) with equation (22) derives
According to the equation set of (26), we can obtain two differential equations of a and θ as Equation (27) denotes that _ a and _ θ have the same order as the parameter ε, and they are slowly varying functions compared to the function _ ψ. It is assumed that the right term of equation (27) can be approximately alternated by its average over a period of ψ, and a and θ are constants in that period. As a result, the averaging equation can be achieved as where α � (2 · a)/l and In order to achieve the steady primary resonance response, letting _ a � 0 and _ θ � 0, we can obtain the amplitudefrequency equation and phase-frequency equation as where B � (2 · f 0 )/(ω 2 0 · l) and s � ω/ω 0 .

Energy Transmissibility (ET).
In the conventional linear vibration isolation theory, force transmissibility is often applied to the vibration isolation performance index. Concerning a nonlinear system, the nonlinear steady response may have some nonexcitation frequency including subharmonics or superharmonics. e force transmissibility may not be applicable to that index. Some literatures [13,14,18] mentioned that the energy transmissibility may be suitable for the vibration isolation performance index. e energy transmissibility is introduced and defined as where P 1 and P 0 are power of transmission force and excitation force, respectively, which can be evaluated by Equation (32) shows that the power of the transmission force can be evaluated in two ways. e first approach needs the displacement and the velocity, which can be obtained with the averaging method theoretically. Meanwhile, as for an experimental test, the second approach is more convenient for the acceleration is easier to be measured. e ET essentially signifies the ratio of root-mean-square of transmitted force over that of external excitation force. It is a generalised form of the force transmissibility. Similarly, if the ET is less than zero, the isolator has positive performance of reducing the energy transmitted from the excitation source to the base.

Experimental Apparatus and Results
An experimental apparatus was designed to validate the MMSI stiffness model; moreover, dynamic experiments were carried out to estimate its vibration isolation Shock and Vibration 5 performance. e section view of the experimental apparatus 3D model is plotted in Figure 2(a). e apparatus consists of four identical accumulators that are embedded within the piston-cylinder container rather than being independent. e container is filled with water so as to intercommunicate all accumulators. e pressure on water may pass to the piston of all accumulators instantly and simultaneously. Each accumulator includes a stair-step piston, a coil spring, a plate, a nut, and a stopper. e piston is confined by the container, the coil spring, and the stopper. It is noticed that the container piston and accumulator pistons are sealed by rubber diaphragms, which may not introduce any friction compared to the conventional seal ring. e container piston is oriented by a ball guide sleeve, while each accumulator piston is oriented by the precompressed force and the rubber diaphragm automatically. e static stiffness property of the MMSI was tested with quasi-static experiments on a fatigue testing machine. e dynamic experiments were carried out with a vibration exciter. e dynamic experimental setup is presented in Figure 2(b); the vibration exciter was fastened on the upper clamping head of the fatigue testing machine, and the MMSI was mounted on the ground. Meanwhile, the payload was fixed on the piston by a threaded rod; the impedance head was fixed between the vibration exciter and the threaded rod by two connecting rods; the acceleration transducer was pasted on the payload. e vibration exciter was controlled by a power amplifier, which was connecting with the dynamic signal analyzer HP35670A. e force and acceleration signals gathered by the impedance head and the acceleration transducer were delivered to the dynamic signal analyzer HP35670A. It employed the sweep sine mode to proceeding dynamic experiments.
Basic parameters of the experimental apparatus applied to the theoretical analysis are listed in Table 1. It is noted that the definition of the effective diameter of a piston is extracted from the theory of rubber diaphragm [19]. As a result, according to the theory above, the theoretical stiffnesses of the bearing segment and the working segment can be derived, and they are listed in the last column of Table 2. e theoretical force-displacement figures are plotted in Figure 3(a). ose curves marked "0," "1," "2," "3," and "4" mean that there are none, one, two, three, and four accumulators being engaged in the isolator. In particular, the curve marked "0" acts as the traditional linear counterpart of those MMSIs with different accumulators. On account of the precompressed force on all identical accumulators, Figure 3(a) shows that the MMSI presents a piecewise linear property of three segments. All curves of MMSIs have the same turning point, which denotes the transformation from the bearing segment to the working segment. However, concerning the transformation from the working segment to the limit segment, they turn at the same force but not at the same displacement. What is more, apparently, the working segment's stiffness gets lower significantly as the number of accumulators increases; however, the working segment's length is widened. In the initial design, the payload is expected to be balanced at the center of the working segment of all MMSIs.
A series of quasi-static experiments on isolators with different number of accumulators were conducted. e loading process and the unloading process are presented in Figure 3(b). Compared to theoretical curves on Figure 3(a), the most obvious difference lies in two points. First, the experimental results show a distinct noncoincidence of the loading process and the unloading process. Second, the experimental slope of the bearing segment is much less than the theoretical one. Looking back to the experimental apparatus as shown in Figure 2(a), we found that only rubber diaphragms may produce both effects on  Shock and Vibration the MMSI stiffness property. First, the rubber diaphragm is inevitable to get some elastic deformation to satisfy the inner water pressure. Second, the rubber in the rubber diaphragm may get some plastic deformation when the loading is much higher so as to produce some hysteresis [19]. Even so, the experimental working stiffnesses of MMSIs fit well with the theoretical one according to Table 2, in which the loading working stiffnesses and the unloading working stiffnesses were extracted from Figure 3(a) and the average stiffnesses of both are listed. After quasi-static experiments, the dynamic experiments were in procedure. In the first place, the theoretical ETs of the isolation system with different isolators were worked out. Based on the experimental apparatus, it neglects the damping in each accumulator while only considering the damping on the piston that supports the payload to conduct the theory calculation. It supposes that there is a linear damping on the piston with the damping coefficient of 200 Ns/m. As shown in Figure 4(a), when the isolation system is under the same relatively low excitation force of 10 N, all dynamic responses appear like a linear system. Meanwhile, compared to the conventional counterpart, the system applying the MMSI with different number of accumulators shows a distinct benefit in the initial frequency of the isolation bandwidth. On the contrary, when the excitation force reaches a relatively higher value of 20 N, an inevitable jump phenomenon may occur according to Figure 4(b). is is a quite common phenomenon for an HSLDS vibration isolator. In the second place, the experimental ETs were also obtained correspondingly. In the test process, we adjusted the magnitude of the power amplifier to change the excitation force. A relatively low signal of 0.5 V produced from power amplifier was applied to conduct the dynamic experiments at first; and the results are shown in Figure 4(c). In comparison with the theoretical results of Figure 4(a), firstly, they both demonstrated that the MMSI is superior to the conventional counterpart in isolating a lower frequency vibration. Secondly, the MMSI possessing more accumulators is better for isolating a much lower frequency vibration.
irdly, the curves marked "0," "1," "2," "3," or "4" dramatically do not fit well. is is because of the application of rubber diaphragm, which may produce much more elastic deformation than the liquid and some plastic deformation. A relatively high signal of 1.0 V was applied to the test to verify the jump phenomenon. As shown in Figure 4(d), the results of both sequential and backward sweeps on the system with the MMSI of one accumulator are presented. It is clear that both results share the low-and high-frequency response, while they do not have the same jump frequency around the resonance frequency.
ese are the typical experimental results of a system with jump phenomenon.  However, there are also some weird properties in contrast to the theoretical results in Figure 4(d), where the experimental results of a higher excitation force get shifted to the left compared to the lower one. We also attribute the mismatching results to the application of the rubber diaphragm. After all, the MMSI is proved to be a typical HSLDS vibration isolator.

Conclusions
is article focused on the stiffness model of the MMSI with several identical accumulators. It found that the stiffness model of MMSI presents a piecewise linear HSLDS property, while the working stiffness was mainly determined by the structure of the accumulator. e working stiffness is biquadratic to the diameter ratio of pistons, directly proportional to the linear spring stiffness, and inversely proportional to the number of accumulators. e length of the working segment is directly proportional to the number of accumulators and the maximum displacement and inversely proportional to the square of the diameter ratio of pistons. A proper design of the length of the segment may allow the payload to always vibrate within the working segment without introducing any nonlinear dynamic phenomenon. Besides, it applies the averaging method to achieve the approximate frequency responses and utilizes the energy transmissibility to evaluate the MMSI vibration isolation system performance.