
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38589564
58469
10.1038/s41598-024-58469-x
Article
Mechanical characteristics analysis of high dimensional vibration isolation systems based on high-static-low-dynamic stiffness technology
Li Bu-yun 12
Shuai Chang-geng 12
Ma Jian-guo mjianguo0722@163.com

12
1 https://ror.org/056vyez31 grid.472481.c 0000 0004 1759 6293 Institute of Noise & Vibration, Naval University of Engineering, Wuhan, 430030 People’s Republic of China
2 National Key Laboratory on Ship Vibration & Noise, Wuhan, 430030 People’s Republic of China
8 4 2024
8 4 2024
2024
14 81957 8 2023
29 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Large floating raft vibration isolation systems (FRVISs) based on high-static-low-dynamic stiffness (HSLDS) technology offer excellent low frequency vibration isolation performance with broad application prospects. However, the design process for these complex high-dimensional coupled nonlinear systems remains poorly developed, particularly when applied for ocean-going vessels that experience rolling and pitching motions. The present work addresses this issue by establishing a six-degree-of-freedom HSLDS vibration isolation model for FRVISs composed of eight isolators, and the model is applied to fully analyze the swing stability and multidimensional vibration isolation performance of these systems. The influence of nonlinearity on the mechanical properties of the vibration isolators is analyzed more clearly by assuming that each vibration isolator realizes nonlinear HSLDS characteristics in the z direction and linear characteristics in the x and y directions. The results demonstrate that the swing displacement responses of the system are greatly reduced under weak nonlinearity, which reflects the high static stiffness and high static stability characteristics of an HSLDS system. The multidimensional vibration isolation performance of the system is evaluated according to the impacts of nonlinearity, the installation height Hz of the isolators, and the relative position Dr of the two middle isolators. The results of analysis demonstrate that applying a value of Hz = 0 produces the best vibration isolation performance overall under strong nonlinearity by avoiding unnecessary secondary peaks in the force transmission rate under harmonic mechanical excitation and ensuring a maximum high-frequency vibration isolation effect. However, applying a weak nonlinearity is better than a strong nonlinearity if Hz is not zero. In contrast, Dr has little effect on the vibration isolation effect of the raft in the x, y, and z directions. Therefore, an equidistant installation with Dr = 0.5 would be considered ideal from the standpoint of installation stability.

Keywords

High-static-low-dynamic stiffness
Quasi-zero stiffness
High dimensional nonlinear system
Floating raft vibration isolation system
Swing stability
Subject terms

Engineering
Mechanical engineering
the National Key Research and Development Program of ChinaHJ2019C020499 issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Floating raft vibration isolation systems (FRVISs) represent a state-of-the-art technology for minimizing the mechanical vibrations of ocean going vessels with a range of goals, such as for enhancing their acoustic stealth performance. Of particular importance in this regard is to maximize the low-frequency vibration isolation performance of these systems1. From a theoretical perspective, the low-frequency vibration isolation performance of these systems increases as their natural frequency decreases2. The primary means of reducing the natural frequency of an FRVIS is to reduce the stiffness of the vibration isolators applied therein. However, this strategy suffers from at least two problems. First, the lowest natural frequency of an isolator is limited by the physical properties, and cannot be reduced indefinitely. Second, system stability is increasingly compromised under decreasing stiffness because this also increases the deformation experienced by the system under the rolling and pitching conditions of the vessel.

These issues have been addressed in recent years by the development of high-static-low-dynamic stiffness (HSLDS) vibration isolation structures, which can reduce dynamic stiffness while ensuring static stability3,4. Ideally, the dynamic stiffness of this type of vibration isolation structure can approach a value of zero, and therefore represents quasi-zero stiffness (QZS) performance5. At present, a number of HSLDS vibration isolator structures have been designed, such as inclined spring structures6–10, cam-roller structures11–15, pneumatic structures16–18, magnetic structures19–22, and structures inspired by biological organisms23–26. Moreover, Li and Xu27,28 were the first to design an FRVIS using QZS isolators, and the vibration isolation performance of the resulting system was analyzed. However, while the low-frequency vibration isolation performance of the system was demonstrated to be substantially improved via the use of QZS isolators, the model established was a little simple, and some factors affecting the vibration isolation effect of the system were not considered, such as the number of isolators, which can be very many in an FRVIS, and their installation positions. In addition, the influence of different excitation conditions and the nonlinearities of the isolators on the vibration isolation performance of FRVISs were also not considered. For example, the vibration isolation characteristics of the system were evaluated under non-eccentric excitations when the isolators realized QZS condition. However, these characteristics are not at all certain under all possible excitation conditions with non-QZS characteristics. The past research results inspired the research of this paper. This paper further studies the situations that have not been considered before to ensure engineering practicability.

Moreover, in contrast to land-based applications, ocean-going vessels experience rolling and pitching motions during operation, which introduce more rigorous stability requirements for HSLDS-FRVISs. However, the swing stability of these complex vibration isolation systems remains poorly evaluated. Meanwhile, some scholars have analyzed the swing stability of linear systems. For example, He et al.29 analyzed the rolling stability of a vibration isolation system designed for a vessel propulsion system based on a single-layer linear vibration isolation model. However, the stability characteristics of such systems are quite different from those of highly nonlinear HSLDS systems. In addition, most studies focused on analyzing the dynamic stability of nonlinear systems have considered relatively simple HSLDS vibration isolator systems, including their nonlinear dynamic behaviors, such as jumping and bifurcation30–32. Accordingly, the design process for these complex high-dimensional coupled nonlinear systems remains poorly developed in ocean-going vessel applications. As a result, the design of efficient, safe, and stable HSLDS-FRVISs remains beyond the reach of the current state of the art.

The present work addresses this issue by establishing a six degrees-of-freedom (6-DOF) HSLDS-FRVIS model, and the model is applied to fully analyze the swing stability and multidimensional vibration isolation performance of these systems. Therefore, the current work lays a sound theoretical foundation for the subsequent design of HSLDS-FRVISs.

Simplified HSLDS-FRVIS model

For the convenience of analysis, it is assumed that the floating raft is rigid and supported only vertically by vibration isolators. For large floating rafts, vibration isolators will be symmetrically arranged at the four corners and the middle of the raft to ensure stability. According to the actual engineering needs, 6, 8, 10 or even more vibration isolators can be symmetrically installed. For general discussion, it is assumed that the floating raft is supported by eight vibration isolators, which is also a common isolator configuration for ship vibration isolation systems. Therefore, the simplified structure of the HSLDS-FRVIS investigated in the present study is illustrated schematically in Fig. 1. As can be seen, the mechanical equipment is rigidly installed on top of the floating raft, and the floating raft is supported by eight HSLDS vibration isolators. The eight vibration isolators are numbered ①–⑧ counterclockwise from the lower right corner. The global coordinate system OXYZ is located at the center of gravity jointly determined by the mechanical equipment and floating raft. The coordinate system of each vibration isolator coincides with the global coordinate system. The influence of nonlinearity on the mechanical properties of the vibration isolators is analyzed more clearly by assuming that each vibration isolator realizes nonlinear HSLDS characteristics in the z direction and linear characteristics in the x and y directions. The other model parameters include lc, bc, and hc, which are one-half of the length, width, and height of the raft, respectively. In addition, installation positions ax, ay1, and az represent the coordinates of the vibration isolators in the coordinate system OXYZ, where ax=lc and ay1=bc. As can be seen, the positions of vibration isolators installed at the four corners of the floating raft structure are usually determined by the length and width of the floating raft. However, the positions of the middle isolators are arbitrarily adjustable. Therefore, we define a coordinate ay2, which represents the position of the middle vibration isolators (i.e., isolators ②, ③, ⑥, and ⑦) relative to the center of gravity along the y axis. In addition, while the installation heights of the isolators are arbitrarily adjustable, we apply a standard position of az=hc herein unless otherwise specified. az is the distance between the upper end of the isolator with rated load and the center of gravity jointly determined by the mechanical equipment and floating raft. Therefore, a standard position of az=hc means that the vibration isolators are installed on the bottom of the floating raft. Furthermore, the installation height of the isolators and the positions of the middle isolators are defined according to an installation height ratio Hz=aziazihchc (i = 1, 2, …, 8) and a distance ratio Dr=aymaymbcbc (m = 2, 3, 6, 7), respectively. The standard values of Hz and Dr applied are 1 and 0.5 respectively. Unless otherwise specified, the structural parameters applied herein are listed in Table 1. The coordinates of an applied excitation force F are denoted as (sx,sy,sz), where the absolute values of sx, sy, and sz respectively represent the eccentric distances of F in the x, y, and z directions. If sx=sy=sz=0, the excitation force is located at the center of gravity of the system (i.e., at O).Figure 1 Schematic illustrating the simplified structure of a high-static-low-dynamic stiffness (HSLDS) floating raft vibration isolation system (FRVIS).

Table 1 Standard structural parameters applied herein for the HSLDS-FRVIS illustrated in Fig. 1.

Parameter	axaxlclc	ay1ay1lclc	ay2ay2lclc	azazlclc	Hz	Dr	
Value	1.0	2.4	1.2	0.07	1	0.5	

The dynamic equation of the system is given as follows33,34. The nonlinearity includes the linear stiffness term and the cubic nonlinear stiffness term, which is a part of a typical duffing equation.1 Mx¨+∑i=1niGiTCi(Gix˙)+∑i=1niGiTKil(Gix)+∑i=1niGiTKin((Gix)⊗(Gix)⊗(Gix))=F

Here, the operation ⊗ is defined as the multiplications of elements in the same position of different matrices. M is a matrix containing the mass m, the moments of inertia Iii (i = x, y, z), and the products of inertia Iij (i, j = x, y, z, i≠j) of the machinery and floating raft, which is defined as follows.2 M=m000000m000000m000000Ixx-Ixy-Ixz000-IyxIyy-Iyz000-Izx-IzyIzz

x=xcyczcαcβcγcT is the displacement vector of the center of gravity, including translation terms ([xc yc zc]T) and angle terms ([α β γ]T), ni is the number of the isolators, Gi=0azi-ayiE3×3-azi0axiayi-axi0 is the position transformation matrix from the upper end of the i-th isolator to O, where E3×3 is a 3×3 identity matrix, and aυi(υ=x,y,z) are the above-defined vibration isolator coordinates, Ci=TiTdiag(cix,ciy,ciz)Ti is the linear damping matrix, where Ti is an identity matrix because the coordinate system of each vibration isolator coincides with the global coordinate system, Kil=TiTdiag(kilx,kily,kilz)Ti is the linear stiffness matrix, and Kin=TiTdiag(0,0,kinz)Ti is the nonlinear stiffness matrix. Accordingly, we can define the excitation force vector as F=FxFyFzMxMyMzT.

Introducing dimensionless parameters yields the following terms:

x^c=xclc, y^c=yclc, z^c=zclc, α^c=αcθ, β^c=βcθ, γ^c=γcθ, ωn=8kzm, τ=ωnt, Ω=ωωn, ξυ=ciυωnkz, k^ilυ=kilυkz, k^inz=kinzlc2kz, F^υ=Fυ8kzlc, M^υ=Mυ8kzlc2, and a^υ=aυlc (υ=x,y,z).

Here, ω is the excitation frequency, θ is a unit angle for rendering angle terms [α β γ]T dimensionless, τ is dimensionless time, ξ is the dimensionless damping ratio and kz is the static stiffness given by the static force divided by static displacement. It is assumed that the damping of an isolator is the same in all directions. For convenience, dimensionless parameters are still represented by original parameter variables except dimensionless time τ and damping ratio ξ.

These dimensionless parameters transform Eq. (1) into a dimensionless dynamic equation of the system. It is further noted that the linear and nonlinear stiffness elements can be alternatively defined as kilz=1-21-l^l^ and kinz=1-l^2l^334,35. The nonlinearity of the system can be conveniently represented by the dimensionless compression factor l^ which is the ratio of the compressed length to the original length for springs, because kinz increases while kilz decreases with increasing l^. Here, the system attains maximum nonlinearity when l^=0.667, and is accordingly a QZS system. In contrast, a value l^=1 represents an equivalent linear system (ELS).

Static swing stability analysis

Swing motion includes rolling and pitching motions. Under rolling motion, the force exerted on the XOZ plane produces only translational displacements in the x and z directions and rotational displacements β around the y axis, as shown in Fig. 1. Regardless of damping, the dimensionless equations defining the rolling motion of the 3-DOF system based on the 6-DOF system described in “Simplified HSLDS-FRVIS model” section are given as follows.3 x¨c-14azβcklx+14klxxc=Fxz¨c+14klzzc+18knz-axβc+zc3+18knzaxβc+zc3=Fzβ¨c+14βcaz2klx+ax2klz-14azklxxc-18axknz-axβc+zc3+18axknzaxβc+zc3=My

Here, klx, klz are the linear stiffness in the x, z directions respectively. knz is the nonlinear stiffness in the z direction.

Under pitching motion, the force exerted on the YOZ plane produces only translational displacements in the y and z directions and rotational displacements α around the x axis, as shown in Fig. 1. Regardless of damping, the dimensionless equations defining the pitching motion of the corresponding 3-DOF system are defined as follows.4 y¨c+12klyyc+12klyazαc=Fyz¨c+12klzzc+34αc2ay12+ay22knzzc+12knzzc3=Fzα¨c+12az2klyαc+14ay12+ay22klz+3knzzc2αc+14ay14+ay24knzαc3+12azklyyc=Mx

Here, kly is the linear stiffness in the y direction.

The swing stability of an FRVIS is usually determined according to its static response because the swing frequency of a ship is quite small, and swing motion can therefore be regarded as a quasi-static process29. In addition, we assume that all vibration isolators realize QZS characteristics in the z direction (i.e., l^=0.667) because this enables us to obtain an analytical expression of the displacement response of the system that facilitates a clear analysis of the influence of system parameters on swing stability.

Ignoring the differential terms in Eq. (3) yields the following rolling displacement response of the FRVIS.5 xc=4Fx+azβklxklxzc=(-4azFx+4axFzaxknz)1/3+(4azFx+4axFzaxknz)1/32βc=(4azFx+4axFzaxknz)1/3-(-4azFx+4axFzaxknz)1/32ax

In contrast, the greater complexity of the equations of pitching motion in Eq. (4) make it impossible to obtain an explicit analytical pitching displacement response for the FRVIS. Therefore, the pitching displacement response of the system is given as follows. This equation can be solved using the solve function of the MATLAB® software.6 12klyyc+12klyαcaz=Fy34αc2ay12+ay22knzzc+12knzzc3=Fz12az2klyαc+34knzay12+ay22zc2αc+14ay14+ay24knzαc3+12azklyyc=0

Similarly, the dimensionless rolling displacement responses xr-l, zr-l, and βr-l obtained for an ELS (l^=1) with only linear stiffness elements Kil=TiTdiag(kz,kz,kz)Ti and Kin=0 (i = 1, 2, …,8) can be given as follows.7 xr-l=4Fx+4Fxaz2ax2zr-l=4Fzβr-l=4azFxax2

The dimensionless pitching displacement responses yp-l, zp-l, and αp-l are also defined as follows.8 yp-l=2Fy+4Fyaz2ay12+ay22zp-l=2Fzαp-l=-4azFyay12+ay22

To further simplify the analysis, we assume that the load is concentrated at the origin O of the global coordinate system (i.e., at the center of gravity of the system). This enables us to define the balance between the load and the support force of the vibration isolator during the rolling and pitching of the ship as being equivalent to applying the following rolling and pitching forces at O.9 Fx=GsinρFz=G(1-cosρ)

10 Fy=GsinϕFz=G(1-cosϕ)

Here, G=Ga8kzlc is the dimensionless load of the system, where Ga is the load, ρ=ρaθ is the dimensionless rolling angle, where ρa is the rolling angle, and ϕ=ϕaθ is the dimensionless pitching angle, where ϕa is the pitching angle.

The impacts of ρ and G on the dimensionless displacement responses of the HSLDS-FRVIS (l^=0.667) and the ELS (l^=1) under rolling conditions are presented in Fig. 2a and b as functions of ρ for G values of 1 and 10, respectively. Similarly, these impacts on the HSLDS-FRVIS and the ELS under pitching conditions are presented in Fig. 3a and b as functions of φ for G values of 1 and 10, respectively. It is found that the responses of the HSLDS-FRVIS and ELS are nearly equivalent in the x and y directions under the swing conditions. This demonstrates that nonlinearity in the z direction has little influence on the system responses in the x and y directions. Under the rolling condition, the response of the HSLDS-FRVIS in the z direction is generally greater than that of the ELS with increasing ρ when G = 1, and is less than that of the ELS only at relatively large values of ρ. Similar behavior is observed under the pitching condition, except that the response of the HSLDS-FRVIS in the z direction is always greater than that of the ELS. However, this behavior changes when G = 10, where the responses of the HSLDS-FRVIS in the z direction are nearly always less than or equal to that of the ELS under both rolling or pitching conditions. Moreover, the gap between the two responses becomes increasingly obvious with increasing ρ or φ. This demonstrates that the displacement offset of the HSLDS-FRVIS is smaller than that of the ELS under large load and large swing angle conditions, which is beneficial for ensuring the stability of the vibration isolation system. This can be qualitatively explained from an analysis of Eqs. (5) and (7) combined with Eqs. (9) and (10) respectively, where we note that the load and swing angle have less influence on the responses of the HSLDS-FRVIS in the z direction than the ELS. Therefore, the HSLDS-FRVIS provides improved swing stability performance over that of the ELS.Figure 2 Dimensionless displacement responses obtained under rolling motion for the HSLDS-FRVIS (l^=0.667) and the ELS (l^=1) defined in Fig. 1 as functions of the dimensionless rolling angle ρ under different dimensionless loads G: (a) G=1; (b) G=10.

Figure 3 Dimensionless displacement responses obtained under pitching motion for the HSLDS-FRVIS and the equivalent linear system (ELS) as functions of the dimensionless pitching angle φ for different G: (a) G=1; (b) G=10.

We also analyzed the influence of the degree of nonlinearity l^ on the displacement responses of the HSLDS-FRVIS and the ELS in the z direction under swing motion for G=10, and the results obtained under roll and pitch motions are presented as functions of ρ and φ in Fig. 4a and b, respectively. It is found that even slight nonlinearity can greatly reduce the displacement response of the system and improve its swing stability under heavy load. In fact, the displacement response observed for l^=0.9 is little different from that obtained under the minimum displacement response at l^=0.667, and is much less than that of the ELS. For example, the responses observed for l^=0.9 at the maximum rolling and pitching angles considered were about 12% and 30% those of the linear system (l^=1.0), respectively. Accordingly, these results reflect the high static stiffness and high static stability characteristics of HSLDS technology.Figure 4 Dimensionless displacement responses obtained in the z direction under swing motions for the HSLDS-FRVIS and ELS as functions of ρ and φ for different degrees of nonlinearity l^ with G=10: (a) rolling motion; (b) pitching motion.

Taking rolling motion as an example, we further analyzed the influence of the installation height ratio Hz of the vibration isolators on the displacement responses of the vibration isolation systems under rolling motion for G=10, and the results obtained are presented for the HSLDS-FRVIS (l^=0.667) and the ELS as functions of ρ in Fig. 5a and b, respectively. It is found that different installation positions have relatively little effect on the displacement response of the HSLDS-FRVIS compared to that of the ELS. We further note that Hz affects only the response of the ELS in the x direction, but has no effect on the response in the z direction. In addition, a value of Hz = 0 (i.e., az = 0) produces the minimum displacement response for the ELS in the x direction. This is more intuitively conveyed from an analysis of Eqs. (7) and (9). Therefore, nonlinear isolators weaken the influence of their installation position on the displacement response of the system compared with that of the ELS.Figure 5 Dimensionless displacement responses obtained under rolling motion for different systems as functions of ρ with G=10 and different vibration isolator installation positions: (a) HSLDS-FRVIS (l^=0.667); (b) ELS.

Multidimensional dynamic analysis

Under the established condition ξ=ξx=ξy=ξz, the dimensionless equations defining the dynamic motion of the 6-DOF system described in “Simplified HSLDS-FRVIS model” section are given as follows.11 x¨c+ξx˙c-ξazβ˙c+klxx-klxazβc=Fxy¨c+ξy˙c+ξazα˙c+klyy+klyazαc=Fyz¨c+ξz˙c+klzzc+knzzc3+32knzay12zcαc2+32knzay22zcαc2+3knzax2zcβc2=Fzα¨c+azklyyc+azξy˙c+az2klyαc+12ay12klzαc+12ay22klzαc+32ay12knzzc2αc+32ay22knzzc2αc+12ay14knzαc3+12ay24knzαc3+12ξay12α˙c+12ξay22α˙c+ξaz2α˙c+32ax2ay12knzαcβc2+32ax2ay22knzαcβc2=Mxβ¨c-azklxxc-ξazx˙c+az2klxβc+ax2klzβc+3ax2knzzc2βc+32ax2ay12knzαc2βc+32ax2ay22knzαc2βc+ax4knzβc3+ξax2β˙c+ξaz2β˙c=Myγ¨c+12ay12klxγc+12ay22klxγc+ax2klyγc+ξax2γ˙c+12ξay12γ˙c+12ξay22γ˙c=Mz

Past research has demonstrated that little difference is observed between solutions obtained by the harmonic balance analytical method and the numerical method in the higher excitation frequency range27. However, these solutions differ greatly in the low frequency range because the analytical solution includes a truncation error that arises as the solution is assumed to be first order, but the actual solution has infinite order terms. Therefore, the solution obtained by the numerical method is more accurate in the low frequency range. Moreover, it is difficult to obtain explicit analytical solutions for high-dimensional nonlinear systems. Collectively considering both accuracy and simplicity, the present work applies the numerical method to solve Eq. (11) by using the ode45 solver of the MATLAB® software, which also facilitates a clear and accurate analysis of the influences of various parameters on the vibration isolation performance of the system.

In addition, the present work assumes that the external excitation is simple harmonic force excitation. Therefore, the vibration isolation effect of an FRVIS is analyzed in the x, y, and z directions according to the force transmission rate, which is defined as12 Tfυ=20lgRMSFtυRMS(Fυ)(dB)

Here, the function lg(∙) represents the logarithm with base 10 of its argument, the function RMS(∙) represents the root mean square of its argument, Ftυ is the nonlinear force transmitted by the FRVIS to the foundation, and Fυ=fυcosΩt is the external harmonic excitation force, where fυ is the dimensionless amplitude of the excitation force. According to Eq. (11), the forces transmitted to the foundation in the x, y, and z directions can be given as follows.13 Ftx=ξx˙c-ξazβ˙c+klxx-klxazβcFty=ξy˙c+ξazα˙c+klyy+klyazαcFtz=ξz˙c+klzzc+knzzc3+32knzay12zcαc2+32knzay22zcαc2+3knzax2zcβc2

These forces can be analyzed further by considering the following two types of excitations.Unidirectional excitation:14 Fx=Fy=0,Fz=0.01cos(Ωt),Mx=My=Mz=0.

Multidirectional eccentric excitation:15 Fx=Fy=0.005cos(Ωt),Fz=0.01cos(Ωt),Mx=syFz-szFy,My=szFx-sxFz,Mz=sxFy-syFx.

As can be seen, the unidirectional excitation mode defined in Eq. (14) facilitates an analysis of the vibration isolation performance of a vibration isolation system in the z direction when the load is concentrated at the origin O of the global coordinate system (i.e., at the center of gravity of the system). We first analyze the influence of l^ on the z-directional vibration isolation performance by plotting the force transmission rate Tz obtained in the z direction as a function of Ω under various values of l^ with ξ=0.06 in Fig. 6. The values of Ω at which Tz = 0 under the four values of l^ are given as Ω1–Ω4. It is found that the vibration isolation effect of the system in the z direction increases with decreasing l^, and achieves a maximum vibration isolation performance at l^=0.667. We also considered the differences between the Ω1–Ω4 values relative to the vibration isolation frequency Ω4 of the linear system, and the results were (Ω4-Ω1)(Ω4-Ω1)Ω4Ω4=76.8%, (Ω4-Ω2)(Ω4-Ω2)Ω4Ω4=30.3% and (Ω4-Ω3)(Ω4-Ω3)Ω4Ω4=12.7%. Accordingly, an increasing degree of nonlinearity is found to expand the range of vibration isolation frequencies and improve the low-frequency vibration isolation performance of the isolator. This can be explained from an analysis of Eq. (11), where we note that the nonlinear stiffness term is zero at l^=1, such that the system response is mainly composed of the dominant harmonic solution. However, the proportion of the subharmonic solution increases with decreasing l^ because this decreases the linear stiffness term and increases the nonlinear stiffness term, and the frequency corresponding to the peak system transmittance shifts to the left, resulting in an expanded vibration isolation frequency band.Figure 6 Force transmission rate Tz obtained in the z direction for an HSLDS-FRVIS as a function of the dimensionless excitation frequency Ω under various values of l^ with a dimensionless damping ratio ξ=0.06.

In contrast to the above unidirectional excitation analysis, the multidirectional excitation mode defined in Eq. (15), which is common in practical engineering settings, enables the influences of nonlinearity and the installation positions of the vibration isolators on the vibration isolation performance of the HSLDS-FRVIS to be analyzed under eccentric loading at coordinates sx, sy, and sz. The eccentricity was set as sx/lc=0.8, sy/bc=0.8, and sz/hc=-1.

We first analyze the influence of Hz on the vibration isolation performance of an HSLDS-FRVIS with l^ = 0.667 and otherwise standard parameters (i.e., those in Table 1 and ξ=0.06) by plotting the force transmission rates Tx, Ty, and Tz obtained in the x, y, and z directions in Fig. 7a–c, respectively, as functions of Ω under various values of Hz. Corresponding results are presented in Figs. 8 and 9 for l^ = 0.8 and l^ = 0.9, respectively. First of all, it is found that Hz has little effect on the vibration isolation performance of the raft in the z direction, and mainly affects the isolation performance in the x and y directions. The relative impacts of Hz in the x, y, and z directions can be explained from an analysis of Eq. (11), which indicates that Hz (i.e., az) mainly affects the vibration response of the raft in the x, y, α, and β directions. Moreover, we note that the apparent impact of Hz on the isolation performance in the x and y directions increases with increasing l^. In fact, the respective plots of Tx and Ty in Fig. 10a and b as functions of Ω for different values of l^ with Hz=0 demonstrate that the degree of nonlinearity has no effect on the response of the raft in the x and y directions under this isolator installation condition. However, the degree of nonlinearity clearly affects the response of the raft in the x and y directions for Hz values of 0.5 and 1. In fact, secondary peaks are observed in the Tx and Ty spectra under these isolator installation conditions, which weakens the vibration isolation performance of the raft. Moreover, the magnitudes of these secondary peaks increase with increasing Hz, which is an increasing detriment to the vibration isolation performance. In addition, the high-frequency vibration isolation effect of the raft in the x and y directions increases with decreasing Hz. Therefore, an HSLDS-FRVIS design with Hz=0 provides an overall optimal vibration isolation effect by avoiding secondary peaks in the Tx and Ty spectra and ensuring the best high-frequency vibration isolation effect in the x and y directions.Figure 7 Force transmission rates Tx, Ty, and Tz obtained for an HSLDS-FRVIS in the x, y, and z directions, respectively, as functions of Ω under various values of Hz with l^=0.667 and otherwise standard conditions: (a) Tx; (b) Ty; (c) Tz.

Figure 8 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of Hz with l^=0.8 and otherwise standard conditions: (a) Tx; (b) Ty; (c) Tz.

Figure 9 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of Hz with l^=0.9 and otherwise standard conditions: (a) Tx; (b) Ty; (c) Tz.

Figure 10 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of l^ with Hz=0 and otherwise standard conditions: (a) Tx; (b) Ty.

However, the mechanisms by which secondary peaks arise in the Tx and Ty spectra require some analysis. This can be ascertained based on the respective plots of Tx and Ty in Fig. 11a and b as functions of Ω for different values of l^ with Hz=1 as an example. We first consider the ELS (l^=1.0), which yields two strong peaks in the Tx and Ty spectra. We must further note that, according to Eq. (11), the responses in the x and y directions are functions solely of a 2-DOF system. The response in the x direction is defined according to the following typical dynamic equations for a 2-DOF system.16 x¨c+ξx˙c-ξazβ˙c+klxx-klxazβc=Fxβ¨c-azklxxc-ξazx˙c+az2klxβc+ax2klzβc+ξax2β˙c+ξaz2β˙c=My

Figure 11 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of l^ with Hz=1 and otherwise standard conditions: (a) Tx; (b) Ty.

The y-direction responses are the same. Therefore, the system generates two transmission peaks in both the x and y directions. We further note that the strongly nonlinear QZS system (l^=0.667) generates a large number of secondary spectra at frequencies less than the primary peak, which detracts from the low-frequency vibration isolation effect of the raft. These secondary spectra arise because the linear stiffness term in the z direction is zero when l^=0.667, and the strong nonlinearity in the z direction leads to subharmonic solutions with large magnitudes in the x and y directions. As can be observed, these secondary peaks do not appear for weakly nonlinear systems (l^=0.8 and l^=0.9), which either exhibit responses that represent a combination of the responses of the ELS and the QZS system for l^=0.8, or responses that are quite similar to the ELS for l^=0.9. This is because the nonlinear stiffness term of the system decreases with increasing l^, while the linear stiffness term increases, and this reduces the magnitude of the subharmonic responses obtained in the solution.

Finally, the influence of the distance ratio Dr on the vibration isolation performance of the HSLDS-FRVIS was analyzed according to the plots of Tx, Ty, and Tz presented as functions of Ω under various values of Dr with l^=0.667 in Fig. 12a–c, respectively. Corresponding results are presented in Fig. 13 for l^ = 0.8. It can be found that Dr has little effect on the vibration isolation performance of the HSLDS-FRVIS in any direction. Therefore, an equidistant installation with Dr=0.5 would be considered ideal from the standpoint of installation stability.Figure 12 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of Dr with l^=0.667 and otherwise standard conditions: (a) Tx; (b) Ty; (c) Tz.

Figure 13 Force transmission rates obtained for an HSLDS-FRVIS as functions of Ω under various values of Dr with l^=0.8 and otherwise standard conditions: (a) Tx; (b) Ty; (c) Tz.

Conclusion

The present work addressed the poorly developed process for designing the complex high-dimensional HSLDS-FRVISs applied in ocean-going vessels by establishing a 6-DOF HSLDS-FRVIS model, and applying that model for fully analyzing the swing stability and multidimensional vibration isolation performance of these systems. The results of extensive analysis demonstrate the following conclusions.Weak nonlinearity greatly reduces the swing displacement responses of the system and improves the swing stability under heavy loads and large swing angles. This reflects the high static stiffness and high static stability characteristics of HSLDS systems. In addition, applying different installation positions for the vibration isolators has little effect on the displacement responses of the HSLDS systems in comparison with that of the corresponding linear systems. Meanwhile, the influence of the installation position on the swing displacement responses of the system decreases with increasing nonlinearity.

In terms of the multidimensional vibration isolation performance of the systems, the low-frequency vibration isolation effect can increase with increasing nonlinearity. However, a strongly nonlinear QZS system is more sensitive to parameter changes than weakly nonlinear HSLDS systems. The results of analysis demonstrate that applying a value of Hz = 0 produces the best vibration isolation performance overall under strong nonlinearity by avoiding unnecessary secondary peaks in the force transmission rate under harmonic mechanical excitation and ensuring a maximum high-frequency vibration isolation effect. However, applying a weak nonlinearity is better than a strong nonlinearity if Hz is not zero. Therefore, a weakly nonlinear HSLDS system is preferred to a highly nonlinear QZS system in engineering design practice to ensure an optimal vibration isolation performance. Meanwhile, Dr has little impact on the vibration isolation effects of the system in x, y, and z directions. Therefore, applying an equidistant installation with Dr=0.5 would be considered ideal from the standpoint of installation stability.

Accordingly, the current work lays a sound theoretical foundation for the subsequent design of HSLDS-FRVISs.

Acknowledgements

This research was funded by the National Key Research and Development Program of China, Grant Number HJ2019C020499.

Author contributions

BL and CS initiated this study and proposed the model of this work; BL designed the structure, performed research, derived the modeling, and wrote the manuscript; JM assist with analysis. All authors read and approved the final manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, Grant Number HJ2019C020499.

Data availability

Some of the data and models generated during the study are available from the corresponding author by reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Howard, C. Q. Recent developments in submarine vibration isolation and noise control. In Proceedings of 1st Submarine Science Technology and Engineering Conference, Adelaide, SA, Australia (2011).
2. Vane, F. A Guide for the Selection and Application of Resilient Mountings to Shipboard Equipment (Revised). David Taylor Model Basin Report 880 (1958).
3. Carrella A Brennan MJ Waters TP Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic J. Sound Vib. 2007 301 678 689 10.1016/j.jsv.2006.10.011
Carrella, A., Brennan, M. J. & Waters, T. P. Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic. J. Sound Vib. 301, 678–689 (2007).10.1016/j.jsv.2006.10.011
4. Carrella A Brennan MJ Waters TP Shin K On the design of a high-static–low-dynamic stiffness isolator using linear mechanical springs and magnets J. Sound Vib. 2008 315 712 720 10.1016/j.jsv.2008.01.046
Carrella, A., Brennan, M. J., Waters, T. P. & Shin, K. On the design of a high-static–low-dynamic stiffness isolator using linear mechanical springs and magnets. J. Sound Vib. 315, 712–720 (2008).10.1016/j.jsv.2008.01.046
5. Alabuzhev P Rivin E Vibration Protecting and Measuring Systems with Quasi-Zero Stiffness 1989 Hemisphere Publishing Corporation
Alabuzhev, P. & Rivin, E. Vibration Protecting and Measuring Systems with Quasi-Zero Stiffness (Hemisphere Publishing Corporation, 1989).
6. Gatti G Shaw AD Gonçalves PJP Brennan MJ On the detailed design of a quasi-zero stiffness device to assist in the realisation of a translational Lanchester damper Mech. Syst. Signal Process. 2022 164 108258 10.1016/j.ymssp.2021.108258
Gatti, G., Shaw, A. D., Gonçalves, P. J. P. & Brennan, M. J. On the detailed design of a quasi-zero stiffness device to assist in the realisation of a translational Lanchester damper. Mech. Syst. Signal Process. 164, 108258 (2022).10.1016/j.ymssp.2021.108258
7. Shaw AD Gatti G Gonçalves PJP Tang B Brennan MJ Design and test of an adjustable quasi-zero stiffness device and its use to suspend masses on a multi-modal structure Mech. Syst. Signal Process. 2021 152 107354 10.1016/j.ymssp.2020.107354
Shaw, A. D., Gatti, G., Gonçalves, P. J. P., Tang, B. & Brennan, M. J. Design and test of an adjustable quasi-zero stiffness device and its use to suspend masses on a multi-modal structure. Mech. Syst. Signal Process. 152, 107354 (2021).10.1016/j.ymssp.2020.107354
8. Ishida S Uchida H Shimosaka H Hagiwara I Design and numerical analysis of vibration isolators with quasi-zero-stiffness characteristics using bistable foldable structures J. Vib. Acoust. 2017 139 031015 10.1115/1.4036096
Ishida, S., Uchida, H., Shimosaka, H. & Hagiwara, I. Design and numerical analysis of vibration isolators with quasi-zero-stiffness characteristics using bistable foldable structures. J. Vib. Acoust. 139, 031015 (2017).10.1115/1.4036096
9. Ishida S Suzuki K Shimosaka H Design and experimental analysis of origami-inspired vibration isolator with quasi-zero-stiffness characteristic J. Vib. Acoust. 2017 139 051004 51012 10.1115/1.4036465
Ishida, S., Suzuki, K. & Shimosaka, H. Design and experimental analysis of origami-inspired vibration isolator with quasi-zero-stiffness characteristic. J. Vib. Acoust. 139, 051004–051012 (2017).10.1115/1.4036465
10. Xu D Zhang Y Zhou J Lou J On the analytical and experimental assessment of the performance of a quasi-zero-stiffness isolator J. Vib. Control 2014 20 2314 2325 10.1177/1077546313484049
Xu, D., Zhang, Y., Zhou, J. & Lou, J. On the analytical and experimental assessment of the performance of a quasi-zero-stiffness isolator. J. Vib. Control 20, 2314–2325 (2014).10.1177/1077546313484049
11. Zhang Y Wei G Wen H Jin D Hu H Design and analysis of a vibration isolation system with cam–roller–spring–rod mechanism J. Vib. Control 2022 28 1781 1791 10.1177/10775463211000516
Zhang, Y., Wei, G., Wen, H., Jin, D. & Hu, H. Design and analysis of a vibration isolation system with cam–roller–spring–rod mechanism. J. Vib. Control 28, 1781–1791 (2022).10.1177/10775463211000516
12. Yao Y Wang X Li H Design and analysis of a high-static-low-dynamic stiffness isolator using the cam-roller-spring mechanism J. Vib. Acoust. 2020 142 021009 10.1115/1.4045583
Yao, Y., Wang, X. & Li, H. Design and analysis of a high-static-low-dynamic stiffness isolator using the cam-roller-spring mechanism. J. Vib. Acoust. 142, 021009 (2020).10.1115/1.4045583
13. Han J Meng L Sun J Design and characteristics analysis of a nonlinear isolator using a curved-mount-spring-roller mechanism as negative stiffness element Math. Probl. Eng. 2018 2018 1359461
Han, J., Meng, L. & Sun, J. Design and characteristics analysis of a nonlinear isolator using a curved-mount-spring-roller mechanism as negative stiffness element. Math. Probl. Eng. 2018, 1359461 (2018).
14. Wang K Zhou J Xu D Sensitivity analysis of parametric errors on the performance of a torsion quasi-zero-stiffness vibration isolator Int. J. Mech. Sci. 2017 134 336 346 10.1016/j.ijmecsci.2017.10.026
Wang, K., Zhou, J. & Xu, D. Sensitivity analysis of parametric errors on the performance of a torsion quasi-zero-stiffness vibration isolator. Int. J. Mech. Sci. 134, 336–346 (2017).10.1016/j.ijmecsci.2017.10.026
15. Zhou J Xu D Bishop S A torsion quasi-zero stiffness vibration isolator J. Sound Vib. 2015 338 121 133 10.1016/j.jsv.2014.10.027
Zhou, J., Xu, D. & Bishop, S. A torsion quasi-zero stiffness vibration isolator. J. Sound Vib. 338, 121–133 (2015).10.1016/j.jsv.2014.10.027
16. Vo NYP Le TD Dynamic analysis of quasi-zero stiffness pneumatic vibration isolator Appl. Sci. 2022 12 2378 10.3390/app12052378
Vo, N. Y. P. & Le, T. D. Dynamic analysis of quasi-zero stiffness pneumatic vibration isolator. Appl. Sci. 12, 2378 (2022).10.3390/app12052378
17. Vo NYP Nguyen MK Le TD Analytical study of a pneumatic vibration isolation platform featuring adjustable stiffness Commun. Nonlinear Sci. Numer. Simul. 2021 98 105775 10.1016/j.cnsns.2021.105775
Vo, N. Y. P., Nguyen, M. K. & Le, T. D. Analytical study of a pneumatic vibration isolation platform featuring adjustable stiffness. Commun. Nonlinear Sci. Numer. Simul. 98, 105775 (2021).10.1016/j.cnsns.2021.105775
18. Shuai C Li B Ma J A novel multi-directional vibration isolation system with high-static–low-dynamic stiffness Acta Mech. 2022 233 5199 5214 10.1007/s00707-022-03387-0
Shuai, C., Li, B. & Ma, J. A novel multi-directional vibration isolation system with high-static–low-dynamic stiffness. Acta Mech. 233, 5199–5214 (2022).10.1007/s00707-022-03387-0
19. Zhou Z Zhou M Dai Z Liu X Li Z Design and experimental validation of a vibration isolator with high-static low-dynamic stiffness and operating point variable property J. Vib. Control 2022 28 1341 1350 10.1177/1077546321990524
Zhou, Z., Zhou, M., Dai, Z., Liu, X. & Li, Z. Design and experimental validation of a vibration isolator with high-static low-dynamic stiffness and operating point variable property. J. Vib. Control 28, 1341–1350 (2022).10.1177/1077546321990524
20. Jiang Y Song C Ding C Xu B Design of magnetic-air hybrid quasi-zero stiffness vibration isolation system J. Sound Vib. 2020 477 115346 10.1016/j.jsv.2020.115346
Jiang, Y., Song, C., Ding, C. & Xu, B. Design of magnetic-air hybrid quasi-zero stiffness vibration isolation system. J. Sound Vib. 477, 115346 (2020).10.1016/j.jsv.2020.115346
21. Yuan S Tunable negative stiffness spring using Maxwell normal stress Int. J. Mech. Sci. 2021 193 106127 10.1016/j.ijmecsci.2020.106127
Yuan, S. et al. Tunable negative stiffness spring using Maxwell normal stress. Int. J. Mech. Sci. 193, 106127 (2021).10.1016/j.ijmecsci.2020.106127
22. Yuan S A tunable quasi-zero stiffness isolator based on a linear electromagnetic spring J. Sound Vib. 2020 482 115449 10.1016/j.jsv.2020.115449
Yuan, S. et al. A tunable quasi-zero stiffness isolator based on a linear electromagnetic spring. J. Sound Vib. 482, 115449 (2020).10.1016/j.jsv.2020.115449
23. Sun B Jing X A tracked robot with novel bio-inspired passive “legs” Robot. Biomim. 2017 4 18 10.1186/s40638-017-0070-6
Sun, B. & Jing, X. A tracked robot with novel bio-inspired passive “legs”. Robot. Biomim. 4, 18 (2017).10.1186/s40638-017-0070-6
24. Jiang G Jing X Guo Y A novel bio-inspired multi-joint anti-vibration structure and its nonlinear HSLDS properties Mech. Syst. Signal Process. 2020 138 106552 10.1016/j.ymssp.2019.106552
Jiang, G., Jing, X. & Guo, Y. A novel bio-inspired multi-joint anti-vibration structure and its nonlinear HSLDS properties. Mech. Syst. Signal Process. 138, 106552 (2020).10.1016/j.ymssp.2019.106552
25. Jing X Zhang L Feng X Sun B Li Q A novel bio-inspired anti-vibration structure for operating hand-held jackhammers Mech. Syst. Signal Process. 2019 118 317 339 10.1016/j.ymssp.2018.09.004
Jing, X., Zhang, L., Feng, X., Sun, B. & Li, Q. A novel bio-inspired anti-vibration structure for operating hand-held jackhammers. Mech. Syst. Signal Process. 118, 317–339 (2019).10.1016/j.ymssp.2018.09.004
26. Dai H Jing X Wang Y Yue X Yuan J Post-capture vibration suppression of spacecraft via a bio-inspired isolation system Mech. Syst. Signal Process. 2018 105 214 240 10.1016/j.ymssp.2017.12.015
Dai, H., Jing, X., Wang, Y., Yue, X. & Yuan, J. Post-capture vibration suppression of spacecraft via a bio-inspired isolation system. Mech. Syst. Signal Process. 105, 214–240 (2018).10.1016/j.ymssp.2017.12.015
27. Li Y Xu D Vibration attenuation of high dimensional quasi-zero stiffness floating raft system Int. J. Mech. Sci. 2017 126 186 195 10.1016/j.ijmecsci.2017.03.029
Li, Y. & Xu, D. Vibration attenuation of high dimensional quasi-zero stiffness floating raft system. Int. J. Mech. Sci. 126, 186–195 (2017).10.1016/j.ijmecsci.2017.03.029
28. Li Y Xu D Force transmissibility of floating raft systems with quasi-zero-stiffness isolators J. Vib. Control 2018 24 3608 3616 10.1177/1077546317708460
Li, Y. & Xu, D. Force transmissibility of floating raft systems with quasi-zero-stiffness isolators. J. Vib. Control 24, 3608–3616 (2018).10.1177/1077546317708460
29. He L Xu W Bu W Shi L Dynamic analysis and design of air spring mounting system for marine propulsion system J. Sound Vib. 2014 333 4912 4929 10.1016/j.jsv.2014.05.045
He, L., Xu, W., Bu, W. & Shi, L. Dynamic analysis and design of air spring mounting system for marine propulsion system. J. Sound Vib. 333, 4912–4929 (2014).10.1016/j.jsv.2014.05.045
30. Xu X Liu H Jiang X Atindana AV Uncertainty analysis and optimization of quasi-zero stiffness air suspension based on polynomial chaos method Chin. J. Mech. Eng. 2022 35 93 10.1186/s10033-022-00758-5
Xu, X., Liu, H., Jiang, X. & Atindana, A. V. Uncertainty analysis and optimization of quasi-zero stiffness air suspension based on polynomial chaos method. Chin. J. Mech. Eng. 35, 93 (2022).10.1186/s10033-022-00758-5
31. Wang X Liu H Chen Y Gao P Beneficial stiffness design of a high-static-low-dynamic-stiffness vibration isolator based on static and dynamic analysis Int. J. Mech. Sci. 2018 142–143 235 244 10.1016/j.ijmecsci.2018.04.053
Wang, X., Liu, H., Chen, Y. & Gao, P. Beneficial stiffness design of a high-static-low-dynamic-stiffness vibration isolator based on static and dynamic analysis. Int. J. Mech. Sci. 142–143, 235–244 (2018).10.1016/j.ijmecsci.2018.04.053
32. Lu Z Chen L Brennan MJ Li J-M Ding H The characteristics of vibration isolation system with damping and stiffness geometrically nonlinear J. Phys. Conf. Ser. 2016 744 012115 10.1088/1742-6596/744/1/012115
Lu, Z., Chen, L., Brennan, M. J., Li, J.-M. & Ding, H. The characteristics of vibration isolation system with damping and stiffness geometrically nonlinear. J. Phys. Conf. Ser. 744, 012115 (2016).10.1088/1742-6596/744/1/012115
33. Harris’ Shock and Vibration Handbook. (McGraw-Hill, 2010).
34. Tang B Brennan MJ On the shock performance of a nonlinear vibration isolator with high-static-low-dynamic-stiffness Int. J. Mech. Sci. 2014 81 207 214 10.1016/j.ijmecsci.2014.02.019
Tang, B. & Brennan, M. J. On the shock performance of a nonlinear vibration isolator with high-static-low-dynamic-stiffness. Int. J. Mech. Sci. 81, 207–214 (2014).10.1016/j.ijmecsci.2014.02.019
35. Carrella A Brennan MJ Waters TP Lopes V Force and displacement transmissibility of a nonlinear isolator with high-static-low-dynamic-stiffness Int. J. Mech. Sci. 2012 55 22 29 10.1016/j.ijmecsci.2011.11.012
Carrella, A., Brennan, M. J., Waters, T. P. & Lopes, V. Force and displacement transmissibility of a nonlinear isolator with high-static-low-dynamic-stiffness. Int. J. Mech. Sci. 55, 22–29 (2012).10.1016/j.ijmecsci.2011.11.012
