Next Article in Journal
Developments in the Use of Hinfinity Control and μ-Analysis for Reducing Vibration in Intelligent Structures
Next Article in Special Issue
Obtaining Vortex Formation in Blood Flow by Particle Tracking: Echo-PV Methods and Computer Simulation
Previous Article in Journal
Optimal Siting and Sizing of Electric Vehicle Energy Supplement Infrastructure in Highway Networks
Previous Article in Special Issue
Mathematical Model of the Working Processes of the Gas Cap of a Piston Pump Installed in the Discharge Line
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Support Stiffness Nonlinearity on the Low-Frequency Vibro-Acoustic Characteristics for a Mechanical Equipment—Floating Raft—Underwater Cylindrical Shell Coupled System

1
School of Mechanical Engineering, Shandong University, Jinan 250061, China
2
Key Laboratory of High Efficiency and Clean Mechanical Manufacture, Shandong University, Ministry of Education, Jinan 250061, China
*
Author to whom correspondence should be addressed.
Inventions 2023, 8(5), 118; https://doi.org/10.3390/inventions8050118
Submission received: 11 August 2023 / Revised: 13 September 2023 / Accepted: 18 September 2023 / Published: 21 September 2023

Abstract

:
The low-frequency vibro-acoustic characteristics of a mechanical equipment—floating raft—cylindrical shell—underwater acoustic field coupled system with nonlinear supports are studied in this paper. Firstly, the state space equations were established by a modal superposition theory for the coupled system, and a modal parameter identification method was deduced and verified for the cylindrical shell—underwater acoustic field coupled subsystem. On this basis, the formulas were derived for transmitted power flow in the coupled system, and the nonlinear stiffness constitutive relation of the vibration isolation supports was expressed by softening and hardening characteristics. Finally, dynamic simulations were carried out by the Runge—Kutta method to analyze the effect of nonlinear stiffness characteristic parameters on the low-frequency vibration modes and vibro-acoustic transfer characteristics in the coupled system. The research shows that a superharmonic phenomenon is common in the steady vibration mode of the coupled system with a nonlinear softening (or hardening) stiffness characteristic under harmonic excitation. The stronger the softening (or hardening) stiffness characteristic is, the more complex the vibration form is, and the smaller (or larger) the low-frequency vibro-acoustic transfer level in resonance regions is.

1. Introduction

Nonlinear vibration mode is an important factor to be considered in the study of low-frequency vibration isolation for underwater cylindrical shell systems. Due to relatively high vibration energy at low frequencies, the elastic support, which presents linear stiffness in an original micro-amplitude vibration mode, will show nonlinear characteristics. Strictly speaking, stiffness characteristics of various elastic vibration isolation elements are nonlinear. In a micro-amplitude vibration mode, a vibration isolator has a certain linear stiffness working interval. Stiffness characteristics are nonlinear if amplitude is sufficiently large.
In recent years, many scholars have studied the vibration characteristics of nonlinear vibration isolation systems. Smirnov et al. [1] studied and compared the efficiency of linear and nonlinear vibration isolation systems. Santhosh [2] studied the dynamic performance of an asymmetric nonlinear vibration isolation mechanism under harmonic excitation. Dutta et al. [3] analyzed a nonlinear vibration isolator with magnetorheological fluid dampers and cubic nonlinear springs to determine the effectiveness of isolation during force or displacement excitation. Suman et al. [4] proposed a nonlinear vibration control device based on a negative stiffness mechanism, which can improve the isolation ability of the vibration suspension system in the low-frequency range. Araki et al. [5] proposed a vertical quasi-zero stiffness (QZS) vibration isolator with adjustable restoring force. Zhao et al. [6,7] successively designed QZS vibration isolators with two and three pairs of inclined springs. Ye et al. [8] proposed a QZS system based on a cam—roller mechanism. Chang et al. [9] proposed a QZS dynamic shock absorber. Bouna et al. [10] studied the vibration control of multi-span girder bridges under the excitation of pier foundation vibration by using nonlinear QZS isolators. Drezet et al. [11] proposed a nonlinear vibration energy harvester based on the concept of high static and low dynamic stiffness (HSLDS). Lu et al. [12,13,14] studied an orthogonal six-degree-of-freedom vibration isolation system with HSLDS and an electromagnetic Stewart platform to suppress the low-frequency vibration of high-precision instruments, provided modal coefficients by improved Fourier series and the Rayleigh—Ritz method, and studied the nonlinear energy transfer of flexible plates with arbitrary boundaries coupled to HSLDS isolators. Sun et al. [15] proposed a method for designing HSLDS mounts based on target force curves, which are obtained by placing negative stiffness structures in parallel with positive stiffness linear springs. Yao et al. [16] proposed an HSLDS isolator with a cam—roller—spring mechanism, in which the cam profile can be specially designed to suit different operating requirements. Chong et al. [17] designed a nonlinear X-combined structure with HSLDS in a large displacement range, which avoided negative stiffness and instability factors and improved its vibration isolation characteristics.
Chaos is a unique dynamical behavior of nonlinear systems. Leutcho et al. [18] systematically studied a new type of chaotic system with single parametric nonlinearity. Yan et al. [19] proposed a novel bi-state nonlinear vibration isolator, whose motion can be either periodic or chaotic. Rahman et al. [20] proposed a new fractional-order chaotic system containing multiple nonlinear terms, which can excite hidden chaotic attractors or self-excited chaotic attractors depending on the chosen system parameters or the value of its fractional-order derivative. Li et al. [21] proposed an SRN feature extraction method combining maximum relevance minimum redundancy, which can effectively classify different chaotic signals. Karimov et al. [22] proposed a novel implementation of the chaotic Duffing oscillator as a simple mechanical system with a translational movement. Dreau et al. [23] introduced a stochastic modeling method for nonlinear systems with response surface discontinuity, which improved the efficiency of multi-variate polynomial chaos expansion.
The power flow characteristic is an important index to study nonlinear vibration isolation. Cassidy et al. [24] proposed a general nonlinear control synthesis method for power flow constrained energy harvesters, which is analytically guaranteed to be superior to the optimal static admittance in the steady state random response. Silva et al. [25] proposed a technique for evaluating and optimizing vibrational power flow involving beam-and-plate coupling. Varghese et al. [26] applied the combined transient power flow balancing and acceleration matching techniques to detect and quantify the crack damage of beam structures at different positions. Yang et al. [27,28] successively studied the power flow and force transmissibility characteristics of a two-degree-of-freedom nonlinear vibration isolation system and a coupled nonlinear oscillator. Ren et al. [29] compared the power transmission between rigid, single-stage isolation and double-stage isolation gearbox installation configurations with the vibration power flow of the system as the evaluation index. Zhang et al. [30] proposed a power flow analysis method for quantitatively evaluating the dynamic performance of a nonlinear energy sink in the frequency domain. Xu et al. [31] analyzed the bending vibration and power flow of axially loaded beams with arbitrary boundaries and nonuniform elastic foundations via energy principle in conjunction with the Rayleigh—Ritz procedure. Mahapatra et al. [32] studied the influence of general coupling conditions on the vibration and power flow characteristics of a two-plate composite plate structure.
There has been a lot of literature on the dynamics of nonlinear vibration isolation systems at home and abroad and, especially, the application of quasi-zero stiffness or high static and low dynamic stiffness nonlinear vibration isolators for low-frequency vibration control has become a research hotspot in the last decade. However, compared with the wide application of nonlinear vibration isolators, the mechanism of vibro-acoustic transfer in nonlinear vibration isolation systems, especially the influence of nonlinear factors in low-frequency vibration modes, has not been fully studied. In this paper, the low-frequency vibro-acoustic characteristics of a mechanical equipment—floating raft—underwater cylindrical shell coupled system with nonlinear stiffness supports are studied. It can not only provide theoretical guidance for low-frequency vibration analysis, diagnosis, and design of dynamics of underwater nonlinear coupled systems but also provide reference for the design of nonlinear vibration isolators.
The rest of this paper is organized as follows. In Section 2, based on the parametric modeling of an admittance method of a mechanical equipment—floating raft—cylindrical shell—underwater acoustic field coupled system, the state space equations of the subsystems are established by using modal identification and modal superposition theories; and then the nonlinear stiffness constitutive relation of the vibration isolation supports is expressed by softening and hardening characteristics, and the calculation method of transmitted power flow is derived. In Section 3, the influence of nonlinear stiffness characteristic parameters on the vibro-acoustic characteristics of the coupled system is studied by using MATLAB numerical simulation. Finally, some conclusions are presented in Section 4.

2. Materials and Methods

2.1. Model of a Mechanical Equipment—Floating Raft—Cylindrical Shell—Underwater Acoustic Field Coupled System

Figure 1 describes the model of a mechanical equipment—floating raft—cylindrical shell—underwater acoustic field coupled system.
As shown in Figure 1a,b, vibration source mechanical equipment A1 and A2 (simulating multiple units) are installed on a floating raft R. Each unit is connected to the raft frame through four elastic supports, and the raft frame is connected to a cylindrical shell through four elastic supports. A rectangular coordinate system Oxyz is established, where the Oz axis is the center line of the cylinder and Oy is the vertical direction. To simplify analysis, only the vertical translational vibration and the rotational vibration around the Oz and Ox axes of the system are considered. Therefore, vertical excitation forces F1, F2 and excitation torques T1, T2 around the Oz axis are applied at the center of mass of each unit. The units A1, A2 and the floating raft R are all rigid bodies symmetrical about the yOz plane. In Figure 1a,b, b1, b2, h1z, h2z, h1x, h2x, Hz, Hx, and LR are the installation size parameters of the system, and σ1, σ2, σ3, and σ4 are the connection points of the lower vibration isolation supports with the cylindrical shell.
Figure 1c abstracts the complex coupled system with nonlinear connecting elements. A (the vibration source subsystem) is subjected to l excitations Fen (n = 1, 2, …, l) and has displacement responses wen at the excitations. A is connected to R (the raft frame) by q elastic connecting elements. At the connection points of the elastic elements with A and R, there are excitation forces FAi and Fti (i = 1, 2, …, q), as well as displacement responses wAi and wti. The elastic elements have constitutive relation fi(δi, δ ˙ i ), where δi = wAiwti. Similarly, R is connected to Y (the cylindrical shell–underwater acoustic field coupled subsystem) by r elastic connecting elements. At the connection points of the elastic elements with R and Y, there are excitation forces Fbj and Fyj (j = 1, 2, …, r), as well as displacement responses wbj and wyj. The elastic elements have constitutive relation gj(κj, κ ˙ j ), where κj = wbjwyj. For specific practical problems, the number of elastic coupled elements in series in Figure 1c can be increased or decreased according to different modeling needs.

2.2. State Space Model of the Coupled System

For the unit A1 in Figure 1, its vertical vibration displacement is represented by y1; its angular displacement around the Ox axis is represented by θ1x; and its angular displacement around the Oy axis is represented by θ1y. Define a state variable uA1 = [u1, u2, u3, u4, u5, u6]T = [y1, θ1x, θ1z, y ˙ 1 , θ ˙ 1 x , θ ˙ 1 z ]T, then the state space of the unit A1 equation is
u ˙ A 1 = u ˙ 1 u ˙ 2 u ˙ 3 u ˙ 4 u ˙ 5 u ˙ 6 = 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 u 1 u 2 u 3 u 4 u 5 u 6 + 0 0 0 I 1 × 4 · F 1 / m 1 Z 1 · F 1 / J 1 x X 1 · F 1 / J 1 z + 0 0 0 F 1 / m 1 0 T 1 / J 1 z = Λ 0 · u A 1 + f A 1 F 1 + f in 1 t
where F1 = [FA1, FA2, FA3, FA4]T is an excitation force column vector of the unit A1 to its four supports; m1, J1x, and J1z are the mass of the unit A1 and the moment of inertia that bypasses its centroid and is parallel to the axis of Ox and Oz, respectively; X1 = [0.5h1x, 0.5h1x, −0.5h1x, −0.5h1x]; Z1 = [−0.5h1z, 0.5h1z, −0.5h1z, 0.5h1z]; I1×4 denotes a 4-dimensional row vector with all elements 1; fA1(F1) denotes a vector function with respect to F1; fin1 is a known external input vector function; Λ0 is a constant matrix.
Similarly, for the unit A2, its vertical vibration displacement is represented by y2; its angular displacement around the Ox axis is represented by θ2x; and its angular displacement around the Oy axis is represented by θ2y. Define a state variable uA2 = [u7, u8, u9, u10, u11, u12]T = [y2, θ2x, θ2z, y ˙ 2 , θ ˙ 2 x , θ ˙ 2 z ]T, then the state space equation of the unit A2 is
u ˙ A 2 = Λ 0 · u A 2 + f A 2 F 2 + f in 2 t
where F2 = [FA5, FA6, FA7, FA8]T is an excitation force column vector of the unit A2 to its four supports; fA2(F2) = [0, 0, 0, −I1×4F2/m2, Z2F2/J2x, −X2F2/J2z]T denotes a vector function with respect to F2; fin2 = [0, 0, 0, F2/m2, 0, T2/J2z]T is a known external input vector function; m2, J2x, and J2z are the mass of the unit A2 and the moment of inertia that bypasses its centroid and is parallel to the axis of Ox and Oz, respectively; X2 = [0.5h2x, 0.5h2x, −0.5h2x, −0.5h2x]; Z2 = [−0.5h2z, 0.5h2z, −0.5h2z, 0.5h2z].
Combining the units A1 and A2 into the vibration source equipment subsystem, Equations (1) and (2) are merged into the following state space equation:
u ˙ A = u ˙ A 1 u ˙ A 2 = Λ 0 O 6 × 6 O 6 × 6 Λ 0 u A 1 u A 2 + f A 1 F 1 f A 2 F 2 + f in 1 t f in 2 t = Λ · u A + f A F 1 , F 2 + f in t
For the raft frame R, its vertical vibration displacement is represented by yR; its angular displacement around the Ox axis is represented by θRx; and its angular displacement around the Oy axis is represented by θRy. Define a state variable uR = [uR1, uR2, uR3, uR4, uR5, uR6]T = [yR, θRx, θRz, y ˙ R , θ ˙ R x , θ ˙ R z ]T, then the state space equation of the raft frame R is
u ˙ R = u ˙ R 1 u ˙ R 2 u ˙ R 3 u ˙ R 4 u ˙ R 5 u ˙ R 6 = 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 u R 1 u R 2 u R 3 u R 4 u R 5 u R 6 + 0 0 0 I 1 × 8 · F t F b / m R Z R 2 · F b Z R 1 · F t / J R x X R 1 · F t X R 2 · F b / J R z = Λ 0 · u R + f R F t , F b
where Ft = [Ft1, Ft2, Ft3, Ft4, Ft5, Ft6, Ft7, Ft8]T and Fb = [Fb1, Fb2, Fb3, Fb4]T are excitation force column vectors of the upper supports to the raft frame and the raft frame to the lower supports, respectively; mR, JRx, and JRz are the mass of the raft frame and the moment of inertia that bypasses its centroid and is parallel to the axis of Ox and Oz, respectively; XR1 = [0.5h1x, 0.5h1x, −0.5h1x, −0.5h1x, 0.5h2x, 0.5h2x, −0.5h2x, −0.5h2x]; XR2 = [0.5Hx, 0.5Hx, −0.5Hx, −0.5Hx]; ZR1 = [−b1 − 0.5h1z, −b1 + 0.5h1z, −b1 − 0.5h1z, −b1 + 0.5h1z, b2 − 0.5h2z, b2 + 0.5h2z, b2 − 0.5h2z, b2 + 0.5h2z]; ZR2 = [−0.5Hz, 0.5Hz, −0.5Hz, 0.5Hz]; fR(Ft, Fb) denotes a vector function about Ft, Fb.
For the connection supports between the vibration source equipment and the raft frame, ignoring the effect of their mass, the power/response transfer relation is
F t i = F A i = f i ( δ i ,   δ ˙ i ) ,   i = 1 ,   2 ,   ,   8
where δi = wAiwti, and there are linear relationships between wAi, wti and the state variables uA and uR:
w A 1 = u 1 I 4 × 1 u 2 Z 1 T + u 3 X 1 T , w ˙ A 1 = u 4 I 4 × 1 u 5 Z 1 T + u 6 X 1 T w A 2 = u 7 I 4 × 1 u 8 Z 2 T + u 9 X 2 T , w ˙ A 2 = u 10 I 4 × 1 u 11 Z 2 T + u 12 X 2 T w t = u R 1 I 8 × 1 u R 2 Z R 1 T + u R 3 X R 1 T , w ˙ t = u R 4 I 8 × 1 u R 5 Z R 1 T + u R 6 X R 1 T
where wA1 = [wA1, wA2, wA3, wA4]T; wA2 = [wA5, wA6, wA7, wA8]T; wt = [wt1, wt2, wt3, wt4, wt5, wt6, wt7, wt8]T.
Let wA = [wA1T, wA2T]T, then Equation (5) is expressed as:
F A = F 1 F 2 = F t = f w A u A w t u R , w ˙ A u A w ˙ t u R = f ¯ u A , u R
Similarly, for the connection supports between the raft frame and the cylindrical shell, ignoring the effect of their mass, the power/response transfer relation is
F b j = F y j = g j ( κ j ,   κ ˙ j ) ,   ( j = 1 ,   2 ,   3 ,   4 )
where κj = wbjwyj and there is a linear relationship between wbj and the state variable uR:
w b = u R 1 I 4 × 1 u R 2 Z R 2 T + u R 3 X R 2 T , w ˙ b = u R 4 I 4 × 1 u R 5 Z R 2 T + u R 6 X R 2 T
where wb = [wb1, wb2, wb3, wb4]T.
The cylindrical shell is a continuum, and its modal parameters of all orders can be solved analytically. Further, a state space equation is established by using the modal parameters. The above is a common method for dealing with continuum problems. However, the cylindrical shell coupled with the external underwater acoustic field becomes a dissipative system, so it is difficult to obtain an analytical expression of the modal parameters through theoretical deduction. In view of the above problem, a discrete numerical solution of displacement admittance functions is first calculated for the cylindrical shell—underwater acoustic field coupled subsystem by referring to Equations (11)–(18) in the literature [33]. Then, analytical expressions are fitted according to the admittance functions data, and then system modal parameters are identified from the expressions. Finally, a state space equation is established by using these modal parameters.
For the modal parameter identification method and modeling process mentioned above, please refer to Appendix A. Here, the state space equation after modeling is directly given for the cylindrical shell—underwater acoustic field coupled subsystem:
u ˙ C = S C u C + Φ ¯ C F y
where uC = [uC1, uC2, uC3, …, uCf]T is a state variable (f is first f-order modes intercepted for the low-frequency vibration problem); Fy = [Fy1, Fy2, Fy3, Fy4]T is an excitation force column vector of the connection supports to the cylindrical shell; SC = diag[sC1, sC2, sC3,…, sCf] is a diagonal matrix composed of poles sCi obtained by modal parameter identification (sCi sorted by |sC1| < |sC2| < |sC3| < … < |sCf|); Φ ¯ C is a matrix composed of identified modal parameter ϕ ¯ i j (i = 1, 2, …, f is modal order; j = 1, 2, 3, 4 is the position of the connection point between the lower support and the cylindrical shell), and
Φ ¯ C = ϕ ¯ 11 ϕ ¯ 12 ϕ ¯ 13 ϕ ¯ 14 ϕ ¯ 21 ϕ ¯ 22 ϕ ¯ 23 ϕ ¯ 24 ϕ ¯ 31 ϕ ¯ 32 ϕ ¯ 33 ϕ ¯ 34 ϕ ¯ f 1 ϕ ¯ f 2 ϕ ¯ f 3 ϕ ¯ f 4
The functional relation between wyj and the state space variable uC is
w y = 2 Re Φ ¯ C T u C = 2 Re Φ ¯ C T · Re u C Im Φ ¯ C T · Im u C
where wy = [wy1, wy2, wy3, wy4]T.
By combining the state space equations (Equations (3), (4), and (10)) of the vibration source equipment, raft frame, and cylindrical shell into one, the total state space equation of the system is formed. Then, through the constitutive equation of the elastic elements (Equations (5) and (8)), the subsystems can be connected into a whole.

2.3. Expression of Nonlinear Stiffness for Vibration Isolation Supports

In Equations (5) and (8), the constitutive relation of the elastic supports is only expressed in general terms, which can be a linear or nonlinear function. For the nonlinear stiffness characteristics of vibration isolation supports, they are generally classified qualitatively into two categories: softening and hardening characteristics, as shown in Figure 2a,b. The blue curve in Figure 2 indicates the nonlinear stiffness characteristics expressed by a piecewise linear method, where δ0 is an approximate linear stiffness working interval control parameter, and K1 and K2 are approximate stiffness control parameters inside and outside the linear working interval, respectively. In practical problems, the transition from linear to nonlinear generally has a gradual process, rather than a sudden change, and the piecewise curve is not easy to deal with in numerical calculation because of its unsmoothness. So, the piecewise turning function is modified to a smooth transition function, as shown by the red curve in Figure 2.
The softening and hardening characteristics are expressed in this paper as Equations (13) and (14), respectively:
δ = τ + K 1 c 1 1 + τ / c 2 3 3 1 1 + K 1 2 sgn τ , F I = K 1 τ c 1 1 + τ / c 2 3 3 1 1 + K 1 2 sgn τ c 1 = K 1 K 2 1 + K 1 K 2 c 2 , c 2 = 1 + K 1 2 δ 0 , τ = τ sgn δ K 1 c 1 / c 2 3 + 1 τ 3 3 1 + K 1 2 δ + K 1 c 1 τ 2 + 3 1 + K 1 2 δ + K 1 c 1 2 τ 1 + K 1 2 δ + K 1 c 1 3 + K 1 c 1 3 = 0
where sgn(τ) is the sign function; c1 and c2 are shape control parameters determined by K1, K2, and δ0; τ is a parameter variable related to the positive real root |τ| of the polynomial equation.
δ = τ K 1 c 1 1 + τ / c 2 3 3 1 1 + K 1 2 sgn τ , F I = K 1 τ + c 1 1 + τ / c 2 3 3 1 1 + K 1 2 sgn τ c 1 = K 2 K 1 1 + K 1 K 2 c 2 , c 2 = 1 + K 1 2 δ 0 , τ = τ sgn δ K 1 c 1 / c 2 3 1 τ 3 + 3 1 + K 1 2 δ K 1 c 1 τ 2 3 1 + K 1 2 δ K 1 c 1 2 τ + 1 + K 1 2 δ K 1 c 1 3 + K 1 c 1 3 = 0
It is assumed that the three groups of vibration isolation supports have the same constitutive relation under the unit A1, the unit A2, and the floating raft R. The K1 of the three groups of vibration isolators are denoted as KA11 = Re{k1}, KA21 = Re{k2}, and KR1 = Re{kR}, respectively, where k1, k2, and kR are the linear complex stiffness. The K2 of the three groups of isolators (denoted as KA12, KA22, and KR2, respectively) is set as α (stiffness ratio) times K1, that is, KA12 = αKA11, KA22 = αKA21, and KR2 = αKR1, where α ∈ [0.1, 3.0] (softening characteristic when α < 1.0 and hardening characteristic when α > 1.0). The δ0 (denoted as δ10, δ20, and δR0, respectively) of the three groups of isolators is estimated to be β times their respective static settlement after installation of the units, that is, δ10 = βm1g/Re{4k1}, δ20 = βm2g/Re{4k2}, and δR0 = β(m1 + m2 + mR)g/Re{4kR}, where β ∈ [0.2, 5.0]. In addition, the structural damping coefficients of the vibration isolators are approximately converted to the equivalent viscous damping coefficients by c1 = Im{k1}/ω1, c2 = Im{k2}/ω1, and cR = Im{kR}/ω1 (where ω1 is the fundamental frequency of the system).

2.4. Evaluation Index of Vibro-Acoustic Transfer Characteristics for the Nonlinear Coupled System

The transmitted power flow (which refers to the vibrational power flow transferred to the cylindrical shell through the supports under the floating raft in this paper) can reflect the process of vibration energy transfer and conversion into underwater acoustic radiation of the coupled system [34], which can be regarded as a comprehensive reflection of the vibro-acoustic transfer characteristics.
For the coupled system shown in Figure 1, the instantaneous power of the excitation Fy acting on the cylindrical shell by the vibration isolation support is
I t = F y T · w ˙ y
The vibrational power flow is the average power of the excitation over a period T = 2π/ω:
P T = 1 T 0 T I t d t
Similar to this, the time average power is introduced:
I ¯ t = 1 t 0 t I τ d τ
The initial value I ¯ t can be obtained by the above equation.
I ¯ 0 = lim t 0 1 t 0 t I τ d τ = lim t 0 I t = I 0 = 0
The derivative of time t on both sides of Equation (17) is calculated:
I ¯ ˙ t = 1 t 2 0 t I τ d τ + I t t = I ¯ t + I t t
From the above equation, combined with Equations (3), (4), (10), and (15), the initial value I ¯ ˙ t can be obtained:
I ¯ ˙ 0 = 1 2 lim t 0 F ˙ y T · w ˙ y + F y T · w ¨ y = 1 2 lim t 0 F ˙ y T · 2 Re Φ ¯ C T S C u C + g T w y w b , w ˙ y w ˙ b · w ¨ y = 0
where the initial values of the system state variables uA, uR, and uC have been set to be 0 (the system is in a static state).
On the other hand, assuming that the vibration enters the steady state after a certain time t0 and has a period T, it can be seen from Equation (17) that for any natural number N, the following equation holds:
t 0 + N T I ¯ t 0 + N T t 0 I ¯ t 0 = t 0 t 0 + N T I τ d τ = 0 N T I τ d τ = N T P T
Therefore, the vibration power flow can be calculated by the time average power:
P T = I ¯ t 0 + T + t 0 T I ¯ t 0 + T I ¯ t 0
If t0 is large enough, the system still does not exhibit perfect periodic oscillation. But as long as there is a mathematical expectation in Equation (16), then Equation (22) still holds. At this time, T can be understood as some kind of virtual value related to the excitation period, and PT still represents the average power level of the excitation. Then, Equation (22) is rewritten as:
P T = T + t 0 I ¯ t 0 + T t 0 I ¯ t 0 / T

3. Results and Discussion

Based on the theories in Section 2, the effect of the nonlinear stiffness control parameters of the vibration isolation supports on the low-frequency vibration modes and vibro-acoustic transfer characteristics of the coupled system are analyzed by numerical simulation methods in the following.
To obtain sufficient amplitude for the system to fully demonstrate its nonlinear vibration characteristics, the amplitude of the excitation force applied to the units A1 and A2 is set to 500 N. The values of the variable stiffness isolators and system structural parameters involved in the numerical simulations are set as shown in Table 1 and Table 2, and other control parameters are shown in Appendix B.

3.1. Modal Parameter Identification for the Cylindrical Shell—Underwater Acoustic Field Coupled Subsystem

Taking the displacement derivative functions H11, H12, and H14 as examples, their fitting spectra obtained from the modal parameter identification method are compared with the numerical spectra obtained by the calculation method in the literature [33], as shown in Figure 3. It can be seen that the fitting results can match well with the original data, and the difference between the two is mainly related to the modal truncation (residual admittance) and the computational error of the original data.
The first 10 orders of poles identified for the cylindrical shell—underwater acoustic field coupled subsystem and their corresponding modal frequencies are listed in Table 3. The values of the modal parameters ϕ ¯ i j at the support connection points are listed in Table 4.

3.2. Effect of Nonlinear Stiffness with Softening Characteristic on the Low-Frequency Vibro-Acoustic Characteristics for the Coupled System

3.2.1. Effect of Stiffness Ratio α

Figure 4 shows the bifurcation diagrams of the average vibration displacement and velocity for the coupled system when β is 1.0 and α is 1.0, 0.9, 0.7, 0.5, 0.3, and 0.1, respectively. Taking α = 0.1 as an example, the discontinuity points near 6.75 Hz in the figure are caused by the system resonating at the natural frequency, while the other discontinuity points are caused by the simulation program setting different steps in different frequency bands to improve the calculation efficiency. The discontinuity points in Figure 4 and in the bifurcation diagram below are the same and will not be explained one by one. It can be seen from Figure 4 that the steady vibration of the coupled system with nonlinear softening stiffness supports is single-periodic motion.
Figure 5 shows the motion phase diagrams of the coupled system corresponding to the resonance frequencies of Figure 4. According to Figure 5a–d, when the softening characteristic is weak (α = 0.9), the units A1, A2 and the raft frame all are presented as approximate harmonic vibration modes. As the stiffness ratio decreases, the motion phase diagrams gradually show irregular closed graphics, indicating that the motions of the units A1, A2 and the raft frame are still single-periodic, but there are superharmonic components. According to Figure 5e, the vibration form at the connection points of the lower vibration isolators with the cylindrical shell is much more complex than that of the units A1, A2 and the raft frame, and the smaller the α, the more complex the vibration form. When α is equal to 0.1, the motion phase diagrams have become disordered or even chaotic. It can also be seen from Figure 5 that in a certain range of stiffness ratio α < 1.0, the resonance displacement amplitudes of the units A1, A2 and the raft frame are larger than those of α = 1.0 (linear stiffness), and the smaller the α, the larger the amplitudes. The resonance velocity amplitudes show the opposite trend. In contrast, at the connection points of the lower vibration isolators with the cylindrical shell, the amplitudes of both the resonance displacement and resonance velocity decrease with the decrease in α. But the above law no longer holds when α decreases to a certain extent, as shown in Figure 5e, which is broken when α is equal to 0.1.
In addition, it can be found from Figure 5 that, when the vibration isolation supports have softening stiffness characteristics, the motion laws at the connection points of the upper vibration isolators with units A1 and A2 are basically the same, and the motion laws at the connection points of the upper and lower vibration isolators with the raft frame are also basically similar. In fact, through the simulation analysis, it can be seen that the above laws are also established when the vibration isolation supports have hardening stiffness characteristics. Therefore, in the following, only the motion phase diagrams are shown and analyzed at the connection points of the upper vibration isolators with the unit A1, the lower vibration isolators with the raft frame, and the lower vibration isolators with the cylindrical shell.
Figure 6 shows the transmitted power flow spectrum through the lower vibration isolators input to the cylindrical shell in the coupled system when β is 1.0 and α is 1.0, 0.9, 0.7, 0.5, 0.3, and 0.1, respectively. As can be seen from Figure 6, compared to the linear stiffness, the softening nonlinear characteristic of the support stiffness makes the resonance region of the system shift to lower frequencies, and the resonance peak and the second peak decrease. The smaller α is, the more significant the above effect is.

3.2.2. Effect of Approximate Linear Working Interval Length Control Parameter β

Figure 7 shows the bifurcation diagrams of the average vibration displacement and velocity for the coupled system when α is 0.5 and β is 5.0, 3.0, 2.0, 1.0, and 0.2, respectively. Comparing Figure 7 with Figure 4, it can be seen that the two are quite similar, and the steady state vibration of the coupled system is still single-periodic motion.
Figure 8 shows the motion phase diagrams of the coupled system corresponding to the resonance frequencies of Figure 7. Since the stiffness ratio α = 0.5 here has determined that the vibration isolation supports have softening stiffness characteristics, reducing the approximate linear working interval will further strengthen their nonlinear stiffness characteristics. The results are that the amplitudes of the resonance displacement and resonance velocity of the unit, the raft frame, and the connection points of the lower vibration isolators with the cylindrical shell generally tend to increase with the decrease in β. From the perspective of motion form, the motion phase diagram of the unit is approximately a closed ellipse, that is, it is approximately a harmonic vibration. The motion phase diagram of the raft frame is irregular and closed, indicating that it contains strong superharmonic components. The motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell is relatively complex, indicating that the superharmonic components are complex and strong, and the phase diagram trajectory tends to be messy or even chaotic due to the enhancement of the nonlinear softening characteristic.
Figure 9 shows the transmitted power flow spectrum in the coupled system when α is 0.5 and β is 5.0, 3.0, 2.0, 1.0, and 0.2, respectively. Comparing Figure 9 with Figure 6, it can be seen that the effect of reducing β is similar to that of reducing α. Both of them enhance the nonlinear softening stiffness effect and make the resonance region move to the lower frequency for the coupled system. Reducing α has a more obvious effect on reducing the fundamental frequency of the system and also weakens the resonance sharpness. But reducing β is more effective for reducing the transmitted power flow at higher frequencies after the fundamental frequency. For example, when α is equal to 0.5 and β is equal to 0.2 in Figure 9, the transmitted power flow level is significantly reduced on the band after the fundamental frequency and, especially, the second peak tends to disappear.

3.3. Effect of Nonlinear Stiffness with Hardening Characteristic on the Low-Frequency Vibro-Acoustic Characteristics for the Coupled System

3.3.1. Effect of Stiffness Ratio α

Figure 10 shows the bifurcation diagrams of the average vibration displacement and velocity for the coupled system when β is 1 and α is 1.0, 1.2, 1.5, 2.0, 2.5, and 3.0, respectively. It can be seen that the steady vibration of the coupled system with nonlinear hardening stiffness supports is dominated by single-periodic motion, but when the nonlinear hardening characteristic is enhanced (α = 2.0, 2.5, 3.0), aperiodic motion (quasi-periodic or chaotic) appears in a frequency band after the fundamental frequency resonance region.
Figure 11 shows the motion phase diagrams of the coupled system corresponding to the resonance frequencies of Figure 10. According to Figure 11a, when the hardening characteristic is weak (α = 1.2, 1.5, 2.0), the unit is presented as an approximate harmonic vibration mode. But as the stiffness ratio further increases (α = 2.5, 3.0), the motion phase diagrams gradually become nonelliptical closed graphics, indicating that the motion of the unit is still single-periodic, but there are superharmonic components. According to Figure 5b,c, the motion phase diagrams of the raft and the connection points of the lower vibration isolators with the cylindrical shell are all irregular closed graphics, indicating that their vibration modes are all single-period motion with superharmonics, and the vibration form at the connection points of the lower vibration isolators with the cylindrical shell is particularly complex. Through the calculation and analysis at different frequencies, it can be found that the above vibration modes and laws are universal in the coupled system with hardening characteristic stiffness supports. It can also be seen from Figure 11 that, when stiffness ratio α > 1.0, the resonance displacement amplitudes of the unit and the raft frame are smaller than those of α = 1.0 (linear stiffness), and the larger the α, the smaller the amplitudes. The resonance velocity amplitudes have no consistent change rule. In contrast, at the connection points of the lower vibration isolators with the cylindrical shell, the amplitudes of both the resonance displacement and resonance velocity increase with the increase in α. But they do not conform when α is equal to 3.0, indicating that the above law only holds in a certain range of α > 1.0.
Figure 12 shows the transmitted power flow spectrum in the coupled system when β is 1.0 and α is 1.0, 1.2, 1.5, 2.0, 2.5, and 3.0, respectively. Comparing Figure 12 with Figure 6, it can be seen that the hardening and softening characteristics of the support stiffness have a significantly opposite effect on the transmitted power flow of the coupled system. The hardening nonlinear characteristic of the support stiffness makes the resonance region of the system shift to higher frequencies, and the resonance peak and the second peak increase slightly. The larger α is, the more significant the above effect is.

3.3.2. Effect of Approximate Linear Working Interval Length Control Parameter β

Figure 13 shows the bifurcation diagrams of the average vibration displacement and velocity for the coupled system when α is 1.5 and β is 5.0, 3.0, 2.0, 1.0, and 0.2, respectively. Comparing Figure 13 with Figure 10, it can be seen that the two are quite similar, and the steady state vibration of the coupled system is always single-periodic motion when α is equal to 1.5.
Figure 14 shows the motion phase diagrams of the coupled system corresponding to the resonance frequencies of Figure 13. In the case of a given α, the amplitudes of the resonance displacement and resonance velocity of the unit and the raft frame generally tend to decrease with the decrease in β. From the perspective of motion form, the motion phase diagram of the unit is approximately a closed ellipse, that is, it is approximately a harmonic vibration. The motion phase diagram of the raft frame has obvious irregularity, indicating that it contains strong superharmonic components. The motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell is the most complex, and with the decrease in β, the stability of the phase diagram orbits deteriorates and tends to be messy or even chaotic.
Figure 15 shows the transmitted power flow spectrum in the coupled system when α is 1.5 and β is 5.0, 3.0, 2.0, 1.0, and 0.2, respectively. Comparing Figure 15 with Figure 12, it can be seen that the effect of reducing β is similar to that of increasing α. Both of them enhance the nonlinear hardening stiffness effect and make the resonance region move to the higher frequency for the coupled system. Increasing α leads to a larger range of fundamental frequency variation of the system. But decreasing β extends the effect of nonlinear stiffness characteristics to higher frequencies. For example, when α is equal to 1.5 and β is equal to 0.2 in Figure 15, the effect of the nonlinear hardening characteristics of the vibration isolators extends significantly to higher frequency bands, increasing the transmitted power flow level on the band after the fundamental frequency for the coupled system.

4. Conclusions

The present study was designed to study the effect of support stiffness nonlinearity on the low-frequency vibro-acoustic characteristics for the mechanical equipment—floating raft–cylindrical shell—underwater acoustic field coupled system. Firstly, the state space equations were established for the coupled system containing the vibration isolation supports with nonlinear stiffness characteristics, and in this process, the modal parameter identification method was deduced theoretically and verified by numerical simulation for the cylindrical shell—underwater acoustic field coupled subsystem. Then, the nonlinear stiffness constitutive relation was expressed in general by the softening and hardening characteristics for the vibration isolation supports, and the calculation formulas of the transmitted power flow in the nonlinear vibration isolation system were deduced. Finally, applying the Runge—Kutta method for the dynamics simulation, the effect of the variation of the nonlinear stiffness characteristics parameters on the low-frequency vibration modes and vibro-acoustic transfer characteristics was analyzed and discussed for the coupled system. The main findings are as follows:
  • For the coupled system with a nonlinear softening stiffness characteristic of the vibration isolation supports, the steady state vibration mode under harmonic excitation is single-periodic motion, but the superharmonic phenomenon is common. And the smaller the stiffness ratio is, the more superharmonic components are and the more complex the vibration form is. Compared with the linear stiffness case, the softening nonlinear stiffness can reduce the fundamental frequency and decrease the low-frequency vibro-acoustic transmission level within and after the resonance region for the vibration isolation system. And the stronger the softening characteristic is, the more significant the above effect is;
  • For the coupled system with a nonlinear hardening stiffness characteristic of the vibration isolation supports, the vibration mode under harmonic excitation is single-periodic motion at most frequencies, and the existence of superharmonics is its main vibration feature. And the larger the stiffness ratio is, the more superharmonic components are and the more complex the vibration form is. In a frequency band slightly higher than the fundamental frequency, it may show a nonperiodic motion mode (quasi-periodic or chaotic). Compared with the linear stiffness case, the hardening nonlinear stiffness can increase the fundamental frequency, broaden the resonance band and slightly raise the low-frequency vibro-acoustic transmission level within and after the resonance region for the vibration isolation system. And the stronger the hardening characteristic is, the more significant the above effect is;
  • For the determined nonlinear stiffness ratio, the variation of the length of the approximate linear working interval has no effect on the periodicity of the vibration of the coupled system. For the coupled system with a nonlinear softening stiffness characteristic of the vibration isolation supports, the smaller the length of the approximately linear working interval is, the lower the low-frequency vibro-acoustic transfer level after the fundamental frequency resonance region is; while for the coupled system with a nonlinear hardening stiffness characteristic of the vibration isolation supports, the smaller the length of the approximately linear working interval is, the higher the low-frequency vibro-acoustic transfer level after the fundamental frequency resonance region is.

Author Contributions

Conceptualization, L.W. and R.H.; methodology, L.W. and R.H.; software, L.W. and R.H.; validation, L.W.; formal analysis, L.W. and R.H.; investigation, L.W.; resources, R.H.; data curation, L.W. and R.H.; writing—original draft preparation, L.W.; writing—review and editing, L.W. and R.H.; visualization, L.W.; supervision, R.H.; project administration, R.H.; funding acquisition, R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. State Space Model of the Cylindrical Shell—Underwater Acoustic Field Coupled Subsystem

Appendix A.1. Solutions of the Admittance Functions for the Coupled Subsystem

A cylindrical coordinate system Orϕz is established on the middle surface of the cylindrical shell. Let the coordinate origin be located at the center of the left end face of the cylindrical shell, r and z be the distance from the spatial point σ(r, ϕ, z) to the central axis of the cylindrical shell and the left end face, respectively, and ϕ be the angle of σ(r, ϕ, z) deviating from the vertical direction. The transformation relationship between the cylindrical coordinates and the rectangular coordinates: x = rcosϕ, y = rsinϕ, z = z.
Let the length, thickness, and central radius of the cylindrical shell be L, d, and a, respectively, and the density, elastic modulus, and Poisson’s ratio of the material of the cylindrical shell be ρ, E, and μ, respectively. A harmonic excitation force Fz⋅ejωtalong the Oz axis direction is applied to the inside of the cylindrical shell at σe(φe, ze). Let u(z, ϕ, t) = U(z, ϕ)⋅ejωt, v(z, ϕ, t) = V(z, ϕ)⋅ejωt, and w(z, ϕ, t) = W(z, ϕ)⋅ejωt denote the vibration displacement of the point σ(ϕ, z) on the cylindrical shell in z, ϕ, and r coordinate directions, respectively. Then, u, v, and w satisfy the following differential equations of motion:
S 1 ( u , v , w ) + ρ d · 2 u / t 2 = F z · δ φ φ e · δ z z e · e j ω t S 2 ( u , v , w ) + ρ d · 2 v / t 2 = 0 S 3 ( u , v , w ) + ρ d · 2 w / t 2 = p ( r , φ , z ) r = a
where p is the radiated underwater acoustic pressure generated outside the cylindrical shell due to its vibration; S1, S2, and S3 are differential operators, and
S 1 U , V , W = E d 1 μ 2 2 z 2 + 1 μ 2 a 2 2 φ 2 U + 1 + μ 2 a 2 V z φ + μ a W z S 2 U , V , W = E d 1 μ 2 1 + μ 2 a 2 U z φ + 1 μ 2 2 z 2 + 1 a 2 2 φ 2 V + 1 a 2 φ 1 d 2 12 2 W S 3 U , V , W = E d 1 μ 2 μ a U z + 1 a 2 φ 1 d 2 12 2 V + 1 a 2 + d 2 12 4 W
where 2 = 2 z 2 + 1 a 2 2 φ 2 , 4 = 4 z 4 + 2 a 2 4 z 2 φ 2 + 1 a 4 4 φ 4 .
The pressure p = P(r, ϕ, z)⋅ejωt of the underwater acoustic field outside the cylindrical shell satisfies the Helmholtz equation:
2 P + k 2 P = 0
There are boundary conditions on the interface between the cylindrical shell and the external underwater acoustic field:
P r r = a = ρ 0 ω 2 j = 0 , k = 1 l = 1 3 q j k i H j k l w cos j φ φ e sin k π z L 0 z L 0 z < 0   or   z > L
where ρ0 is the density of the acoustic field medium; Hjk (j, k = 1, 2, 3, 4) denotes a displacement admittance function, which means that when a unit excitation is applied along the direction of the support force at the connection point σk of the support and the cylindrical shell, the displacement response is generated at the connection point σj.
The displacement admittance matrix of the cylindrical shell—underwater acoustic field coupled subsystem is denoted by:
H j ω = H j k 4 × 4
For the calculation of the admittance functions of the cylindrical shell—underwater acoustic field coupled subsystem, see also Equations (11)–(18) in [33]. It can be seen from these equations that, due to the complexity of the calculation formulas, only discrete numerical solutions of the admittance functions can be given. Figure A1 shows examples of calculated results for the admittance functions of the cylindrical shell along the support direction (taking H11, H12, and H14 as examples).
Figure A1. Examples of the calculation results for the admittance along the support direction: (a) H11; (b) H12; (c) H14.
Figure A1. Examples of the calculation results for the admittance along the support direction: (a) H11; (b) H12; (c) H14.
Inventions 08 00118 g0a1

Appendix A.2. Modal Parameter Identification of the Coupled Subsystem

Since the main focus is on low-frequency vibration problems in this paper, only the first f-order modes are intercepted. Also, the admittance function is expressed as a ratio of two polynomial functions:
H j k j ω a m j ω m 1 + a m 1 j ω m 2 + + a 2 j ω + a 1 j ω n 1 + b n 1 j ω n 2 + + b 2 j ω + b 1
where j, k = 1, 2, 3, 4; n = 2f + 1, m = 2f − 1.
Let H ^ j k j ω be the discrete numerical spectrum of Hjk(jω) obtained by theoretical calculation, and for a number of frequencies ω1, ω2,ω3, …, ωN, there are H ^ j k , i = H ^ j k j ω i , (i = 1, 2, …, N). Noting that a = [a1, a2, …, am]T, b = [b1, b2, …,bn−1]T, pi = [1, jωi, …, (jωi)m−1], qi= H ^ j k , i [1, jωi, …, (jωi)n−2], hi = (jωi)n−1  H ^ j k , i , ei = piaqib− hi, P = [p1, p2, …, ps]T, Q = [q1, q2, …, qs]T, and h = [h1, h2, …, hs]T, the matrix expression and the total mean squared error can be obtained:
e = P a Q b h
J = e T e = x T A R x 2 B x + h 2
where the superscript ‘*’ represents conjugation relation; x = a b , A = P T P * P T Q * Q T P * Q T Q * , B = h R T P R + h I T P I h R T Q R + h I T Q I , PR and PI are the real and imaginary parts of P, respectively; QR and QI are the real and imaginary parts of Q, respectively; hR and hI are the real and imaginary parts of h, respectively; AR is the real part of A.
Consider the total mean squared error J as a function of x. To obtain a minimum of J, let the gradient of Equation (A8) to x be 0:
J = A R x B T = 0 x = A R 1 B T
The polynomial coefficients of the admittance function numerator and denominator can be fitted using the above equation.

Appendix A.3. Establishment of State Space Equation for the Coupled Subsystem

For a continuous structure, its vibration displacement response w(σ, t) can be expressed in the form of vibration model superposition:
w σ , t = i = 1 ψ i σ q i t
where ψi(σ) is the ith-order vibration mode function and qi(t) is the ith-order modal coordinate.
A modal damping ratio ξi is introduced, and the modal coordinate equation is obtained by using the orthogonality of the vibration mode function:
q ¨ i + 2 ξ i ω i q ˙ i + ω i 2 q i = j = 1 4 F j t ψ ¯ i σ j ,   ( i =   1 ,   2 ,   )
where ωi is the i-th order modal frequency; Fj(t) is the concentrated excitation force of the cylindrical shell subjected to the support at point σj (j = 1, 2, 3, 4); ψ ¯ i = ψ i / M   i is the regular vibration mode function; M i = Σ m Σ ψ i 2 σ d σ is the i-th order modal mass, mΣ is the areal density.
According to Equations (A10) and (A11), the admittance function can be expressed as:
H j k j ω i = 1 f ψ ¯ i σ j ψ ¯ i σ k j 2 1 ξ i 2 ω i 1 j ω + ξ i ω i j 1 ξ i 2 ω i 1 j ω + ξ i ω i + j 1 ξ i 2 ω i
According to the above relation, if H(jω) is known, each Hjk(jω) has the same poles s i = ξ i ω i + j 1 ξ i 2 ω i , (i = 1, 2, …, f). Since |si| = ωi, the poles are sorted by|s1| < |s2| < |s3| <…, and the modal parameters are determined as follows:
ω i = s i ,   ξ i = Re s i / s i
ψ ¯ i σ 1 = j 2 Im s i · H 11 · j ω s i j ω = s i
ψ ¯ i σ j = j 2 Im s i · H j 1 · j ω s i j ω = s i ψ ¯ i σ 1 = j 2 Im s i H 11 · j ω s i j ω = s i H j 1 · j ω s i j ω = s i
Substituting the above relations into Equations (A11) and (A10), the modal coordinate equation and the modal superposition formula of vibration displacement expressed in terms of the admittance functions are obtained:
q ¨ i 2 Re s i q ˙ i + s i 2 q i = j 2 Im s i j = 1 4 F j t H j 1 · j ω s i j ω = s i H 11 · j ω s i j ω = s i ,   ( i =   1 ,   2 ,   ,   f )
w σ j , t = i = 1 f q i t j 2 Im s i H 11 · j ω s i j ω = s i H j 1 · j ω s i j ω = s i ,   ( j =   1 ,   2 ,   3 ,   4 )
The above theories are based on proportional damping assumption, and the vibration mode function is a real function by default. But applying Equation (A14) in practical analysis mostly yields complex numbers, that is, complex modes appear. If the complex modes are substituted into Equations (A11) or (A16), the modal coordinate equation will be made to have complex terms, and ultimately the displacement response obtained from Equation (A17) will also become a complex number, which is inconsistent with its physical meaning. Considering the universality of complex modes in this kind of problem, Equations (A11), (A12), and (A17) are modified as:
H j k j ω i = 1 f ϕ ¯ i j ϕ ¯ i k j ω s i + i = 1 f ϕ ¯ i j * ϕ ¯ i k * j ω s i *
υ ˙ 1 υ ˙ f υ ˙ 1 * υ ˙ f * s 1 s f s 1 * s f * υ 1 υ f υ 1 * υ f * = j = 1 4 ϕ ¯ 1 j F j t j = 1 4 ϕ ¯ f j F j t j = 1 4 ϕ ¯ 1 j * F j t j = 1 4 ϕ ¯ f j * F j t
w σ j , t = i = 1 f 2 Re υ i t · Re ϕ ¯ i j Im υ i t · Im ϕ ¯ i j ,   ( j =   1 ,   2 ,   3 ,   4 )
where the poles si and si* (i = 1, 2, …, f) are the f pairs of conjugate latent roots obtained in solving the eigenvalue problem for the state space equation; υi is obtained by variable substitution q = Φυ;Φ is a complex modal matrix consisting of Φi and Φi*; Φi and Φi* are complex eigenvectors corresponding to si and si*; ϕ ¯ i j is obtained by using the residue method, ϕ ¯ i j = H j 1 · j ω s i j ω = s i H 11 · j ω s i j ω = s i and ϕ ¯ i j * = H j 1 · j ω s i * j ω = s i * H 11 · j ω s i * j ω = s i * .
The upper part and the lower part of Equation (A19) are repeated (conjugate relation), so only the upper part is abbreviated to a state space equation:
u ˙ C = S C u C + Φ ¯ C F y
where uC = [uC1, uC2, uC3, …, uCf]T is a state variable; Fy = [Fy1, Fy2, Fy3, Fy4]T is an excitation force column vector of the connection supports to the cylindrical shell; SC = diag[sC1, sC2, sC3,…, sCf] is a diagonal matrix composed of poles (sorted by |sC1| < |sC2| < |sC3| < … < |sCf|); Φ ¯ C is a matrix consisting of the modal parameter ϕ ¯ i j , and
Φ ¯ C = ϕ ¯ 11 ϕ ¯ 12 ϕ ¯ 13 ϕ ¯ 14 ϕ ¯ 21 ϕ ¯ 22 ϕ ¯ 23 ϕ ¯ 24 ϕ ¯ 31 ϕ ¯ 32 ϕ ¯ 33 ϕ ¯ 34 ϕ ¯ f 1 ϕ ¯ f 2 ϕ ¯ f 3 ϕ ¯ f 4
The functional relation between wyj and the state space variable uC is
w y = 2 Re Φ ¯ C T u C = 2 Re Φ ¯ C T · Re u C Im Φ ¯ C T · Im u C
where wy = [wy1, wy2, wy3, wy4]T.

Appendix B. Control Parameters in Simulation Calculation

MATLAB is used in the simulation calculation in this paper. The specific theory of mathematical operations has been described in detail in Section 2 and Appendix A. Here, some control parameters related to computational convergence are supplemented.
As mentioned above, the displacement response of the coupled system at any time t1 is wA1 = [wA1, wA2, wA3, wA4]T, wA2 = [wA5, wA6, wA7, wA8]T, wt = [ wt1, wt2, wt3, wt4, wt5, wt6, wt7, wt8]T, wb = [wb1, wb2, wb3, wb4]T, and wy = [wy1, wy2, wy3, wy4]T. The above displacement responses are assembled into a total displacement vector w = w A 1 T ,   w A 2 T ,   w t T ,   w b T ,   w y T T , whose corresponding velocity vector is w ˙ . Suppose there is a time t2 = t1 + T, the displacement response is w , and the velocity response is w ˙ . Considering w and w ˙ as two points in phase space, according to Poincaré mapping theory, if these two points can coincide, then the system has a period T. However, due to the error caused by the precision of numerical calculation, it is difficult to strictly meet this condition. Therefore, a tolerance criterion ε is set so that:
w w 2 w 2 < ε
If the displacement responses w and w satisfy the above condition, the calculation is considered to be convergent, that is, the motion of the system is periodic (T is a period). Otherwise, it is divergent, that is, the motion of the system is nonperiodic (quasi-periodic or chaotic).
Similarly, the velocity responses w ˙ and w ˙ should also satisfy the condition, that is:
w ˙ w ˙ 2 w ˙ 2 < ε
The above tolerance criterion ε is set to 1.0 × 10−5 in the simulation example of this paper.
Figure A2 shows the convergence process of simulation calculation under α = 0.5 and β = 1.0 as an example of the connection points of the upper vibration isolators with unit A1. It can be seen from Figure A2a,b that convergence is obtained when t = 1.7 s. Figure A2c shows a sequence of state as it move to the attractor.
Figure A2. The convergence process of the calculation at the connection points of the upper vibration isolators with the unit A1 (α = 0.5, β = 1.0): (a) time that the displacement response calculation took to obtain convergence; (b) time that the velocity response calculation took to obtain convergence; (c) the sequence of state as it moves to the attractor.
Figure A2. The convergence process of the calculation at the connection points of the upper vibration isolators with the unit A1 (α = 0.5, β = 1.0): (a) time that the displacement response calculation took to obtain convergence; (b) time that the velocity response calculation took to obtain convergence; (c) the sequence of state as it moves to the attractor.
Inventions 08 00118 g0a2
Similarly, at the connection points of the upper vibration isolators with the unit A2, the upper vibration isolators with the raft frame, the lower vibration isolators with the raft frame, and the lower vibration isolators with the cylindrical shell, and under other nonlinear stiffness conditions, simulation calculation of the periodic motion has similar convergence properties, which will not be shown here one by one.

References

  1. Smirnov, V.; Mondrus, V. Comparison of linear and nonlinear vibration isolation system under random excitation. Procedia Eng. 2016, 153, 673–678. [Google Scholar] [CrossRef]
  2. Santhosh, B. Dynamics and performance evaluation of an asymmetric nonlinear vibration isolation mechanism. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 169. [Google Scholar] [CrossRef]
  3. Dutta, S.; Chakraborty, G. Performance analysis of nonlinear vibration isolator with magneto-rheological damper. J. Sound Vib. 2014, 333, 5097–5114. [Google Scholar] [CrossRef]
  4. Suman, S.; Balaji, P.S.; Selvakumar, K.; Kumaraswamidhas, L.A. Nonlinear vibration control device for a vehicle suspension using negative stiffness mechanism. J. Vib. Eng. Technol. 2021, 9, 957–966. [Google Scholar] [CrossRef]
  5. Araki, Y.; Asai, T.; Kimura, K.; Maezawa, K.; Masui, T. Nonlinear vibration isolator with adjustable restoring force. J. Sound Vib. 2013, 332, 6063–6077. [Google Scholar] [CrossRef]
  6. Zhao, F.; Ji, J.C.; Ye, K.; Luo, Q. Increase of quasi-zero stiffness region using two pairs of oblique springs. Mech. Syst. Signal Process. 2020, 144, 106975. [Google Scholar] [CrossRef]
  7. Zhao, F.; Ji, J.C.; Ye, K.; Luo, Q. An innovative quasi-zero stiffness isolator with three pairs of oblique springs. Int. J. Mech. Sci. 2021, 192, 106093. [Google Scholar] [CrossRef]
  8. Ye, K.; Ji, J.C.; Brown, T. Design of a quasi-zero stiffness isolation system for supporting different loads. J. Sound Vib. 2020, 471, 115198. [Google Scholar] [CrossRef]
  9. Chang, Y.; Zhou, J.; Wang, K.; Xu, D. A quasi-zero-stiffness dynamic vibration absorber. J. Sound Vib. 2021, 494, 115859. [Google Scholar] [CrossRef]
  10. Bouna, H.S.; Nbendjo, B.R.N.; Woafo, P. Isolation performance of a quasi-zero stiffness isolator in vibration isolation of a multi-span continuous beam bridge under pier base vibrating excitation. Nonlinear Dyn. 2020, 100, 1125–1141. [Google Scholar] [CrossRef]
  11. Drezet, C.; Kacem, N.; Bouhaddi, N. Design of a nonlinear energy harvester based on high static low dynamic stiffness for low frequency random vibrations. Sens. Actuators A Phys. 2018, 283, 54–64. [Google Scholar] [CrossRef]
  12. Lu, Z.; Wu, D.; Ding, H.; Chen, L. Vibration isolation and energy harvesting integrated in a Stewart platform with high static and low dynamic stiffness. Appl. Math. Model. 2021, 89, 249–267. [Google Scholar] [CrossRef]
  13. Hao, R.; Lu, Z.; Ding, H.; Chen, L. Orthogonal six-DOFs vibration isolation with tunable high-static-low-dynamic stiffness: Experiment and analysis. Int. J. Mech. Sci. 2022, 222, 107237. [Google Scholar] [CrossRef]
  14. Hao, R.; Lu, Z.; Ding, H.; Chen, L. A nonlinear vibration isolator supported on a flexible plate: Analysis and experiment. Nonlinear Dyn. 2022, 108, 941–958. [Google Scholar] [CrossRef]
  15. Sun, Y.; Zhou, J.; Thompson, D.; Yuan, T.; Gong, D.; You, T. Design, analysis and experimental validation of high static and low dynamic stiffness mounts based on target force curves. Int. J. Non-Linear Mech. 2020, 126, 103559. [Google Scholar] [CrossRef]
  16. Yao, Y.; Li, H.; Li, Y.; Wang, X. Analytical and experimental investigation of a high-static-low-dynamic stiffness isolator with cam-roller-spring mechanism. Int. J. Mech. Sci. 2020, 186, 105888. [Google Scholar] [CrossRef]
  17. Chong, X.; Wu, Z.; Li, F. Vibration isolation properties of the nonlinear X-combined structure with a high-static and low-dynamic stiffness: Theory and experiment. Mech. Syst. Signal Process. 2022, 179, 109352. [Google Scholar] [CrossRef]
  18. Leutcho, G.D.; Kengne, J. A unique chaotic snap system with a smoothly adjustable symmetry and nonlinearity: Chaos, offset-boosting, antimonotonicity, and coexisting multiple attractors. Chaos Solitons Fractals 2018, 113, 275–293. [Google Scholar] [CrossRef]
  19. Yan, B.; Ma, H.; Jian, B.; Wang, K.; Wu, C. Nonlinear dynamics analysis of a bi-state nonlinear vibration isolator with symmetric permanent magnets. Nonlinear Dyn. 2019, 97, 2499–2519. [Google Scholar] [CrossRef]
  20. Rahman, Z.A.S.A.; Jasim, B.H. Hidden Dynamics Investigation, Fast Adaptive Synchronization, and Chaos-Based Secure Communication Scheme of a New 3D Fractional-Order Chaotic System. Inventions 2022, 7, 108. [Google Scholar] [CrossRef]
  21. Li, Y.; Yi, Y.; Wu, J.; Gu, Y. A novel feature extraction method for ship-radiated noise based on hierarchical refined composite multi-scale dispersion entropy-based Lempel-Ziv complexity. Deep Sea Res. Part I Oceanogr. Res. Pap. 2023, 199, 104111. [Google Scholar] [CrossRef]
  22. Karimov, A.; Rybin, V.; Dautov, A.; Karimov, T.; Bobrova, Y.; Butusov, D. Mechanical Chaotic Duffing System with Magnetic Springs. Inventions 2023, 8, 19. [Google Scholar] [CrossRef]
  23. Dreau, J.; Magnain, B.; Batailly, A. Multi-element polynomial chaos expansion based on automatic discontinuity detection for nonlinear systems. J. Sound Vib. 2023, 567, 117920. [Google Scholar] [CrossRef]
  24. Cassidy, I.L.; Scruggs, J.T. Nonlinear stochastic controllers for power-flow-constrained vibratory energy harvesters. J. Sound Vib. 2013, 332, 3134–3147. [Google Scholar] [CrossRef]
  25. Silva, O.M.; Neves, M.M.; Jordan, R.; Lenzi, A. An FEM-based method to evaluate and optimize vibration power flow through a beam-to-plate connection. J. Braz. Soc. Mech. Sci. Eng. 2017, 39, 413–426. [Google Scholar] [CrossRef]
  26. Varghese, C.K.; Shankar, K.K. Damage identification using combined transient power flow balance and acceleration matching technique. Struct. Control. Health Monit. 2014, 21, 135–155. [Google Scholar] [CrossRef]
  27. Yang, J.; Xiong, Y.; Xing, J. Vibration power flow and force transmission behaviour of a nonlinear isolator mounted on a nonlinear base. Int. J. Mech. Sci. 2016, 115, 238–252. [Google Scholar] [CrossRef]
  28. Shi, B.; Yang, J. Quantification of vibration force and power flow transmission between coupled nonlinear oscillators. Int. J. Dyn. Control. 2020, 8, 418–435. [Google Scholar] [CrossRef]
  29. Ren, Y.; Chang, S.; Liu, G.; Wu, L.; Wang, H. Vibratory power flow analysis of a gear-housing-foundation coupled system. Shock. Vib. 2018, 2018, 5974759. [Google Scholar] [CrossRef]
  30. Zhang, Y.; Zhou, L.; Wang, S.; Yang, T.; Chen, L. Vibration power flow characteristics of the whole-spacecraft with a nonlinear energy sink. J. Low Freq. Noise Vib. Act. Control 2019, 38, 341–351. [Google Scholar] [CrossRef]
  31. Xu, D.; Du, J.; Zhao, Y. Flexural vibration and power flow analyses of axially loaded beams with general boundary and non-uniform elastic foundations. Adv. Mech. Eng. 2020, 12, 1687814020921719. [Google Scholar] [CrossRef]
  32. Mahapatra, K.; Panigrahi, S.K. Effect of general coupling conditions on the vibration and power flow characteristics of a two-plate built-up plate structure. Mech. Based Des. Struct. Mach. 2021, 49, 841–861. [Google Scholar] [CrossRef]
  33. Liu, S.; Huo, R.; Wang, L. Vibroacoustic Transfer Characteristics of Underwater Cylindrical Shells Containing Complex Internal Elastic Coupled Systems. Appl. Sci. 2023, 13, 3994. [Google Scholar] [CrossRef]
  34. More, S.; Padmanabhan, C. Power flow based analysis of a floating raft vibration isolation system. Am. Phys. Soc. 2014, 72, 691–708. [Google Scholar]
Figure 1. The model of a mechanical equipment—floating raft—cylindrical shell—underwater acoustic field coupled system: (a) −x direction; (b) z direction; (c) complex coupled system with nonlinear connecting elements.
Figure 1. The model of a mechanical equipment—floating raft—cylindrical shell—underwater acoustic field coupled system: (a) −x direction; (b) z direction; (c) complex coupled system with nonlinear connecting elements.
Inventions 08 00118 g001
Figure 2. Nonlinear stiffness characteristics and their control parameters: (a) softening characteristic; (b) hardening characteristic.
Figure 2. Nonlinear stiffness characteristics and their control parameters: (a) softening characteristic; (b) hardening characteristic.
Inventions 08 00118 g002
Figure 3. Comparison of the fitting spectra and the numerical spectra: (a) H11; (b) H12; (c) H14.
Figure 3. Comparison of the fitting spectra and the numerical spectra: (a) H11; (b) H12; (c) H14.
Inventions 08 00118 g003
Figure 4. Effect of the strength variation of the nonlinear softening stiffness characteristic on the periodic bifurcation properties for the coupled system (β = 1.0): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Figure 4. Effect of the strength variation of the nonlinear softening stiffness characteristic on the periodic bifurcation properties for the coupled system (β = 1.0): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Inventions 08 00118 g004
Figure 5. Effect of the strength variation of the nonlinear softening stiffness characteristic on the low-frequency vibration modes for the coupled system (β = 1.0): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the upper vibration isolators with the unit A2; (c) motion phase diagram at the connection points of the upper vibration isolators with the raft frame; (d) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (e) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Figure 5. Effect of the strength variation of the nonlinear softening stiffness characteristic on the low-frequency vibration modes for the coupled system (β = 1.0): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the upper vibration isolators with the unit A2; (c) motion phase diagram at the connection points of the upper vibration isolators with the raft frame; (d) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (e) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Inventions 08 00118 g005
Figure 6. Effect of the strength variation of the nonlinear softening stiffness characteristic on the vibro-acoustic transfer characteristics for the coupled system (β = 1.0).
Figure 6. Effect of the strength variation of the nonlinear softening stiffness characteristic on the vibro-acoustic transfer characteristics for the coupled system (β = 1.0).
Inventions 08 00118 g006
Figure 7. Effect of the length variation of the approximate linear working interval on the periodic bifurcation properties for the coupled system (α = 0.5): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Figure 7. Effect of the length variation of the approximate linear working interval on the periodic bifurcation properties for the coupled system (α = 0.5): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Inventions 08 00118 g007
Figure 8. Effect of the length variation of the approximate linear working interval on the low-frequency vibration modes for the coupled system (α = 0.5): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Figure 8. Effect of the length variation of the approximate linear working interval on the low-frequency vibration modes for the coupled system (α = 0.5): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Inventions 08 00118 g008
Figure 9. Effect of the length variation of the approximate linear working interval on the vibro-acoustic transfer characteristics for the coupled system (α = 0.5).
Figure 9. Effect of the length variation of the approximate linear working interval on the vibro-acoustic transfer characteristics for the coupled system (α = 0.5).
Inventions 08 00118 g009
Figure 10. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the periodic bifurcation properties for the coupled system (β = 1.0): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Figure 10. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the periodic bifurcation properties for the coupled system (β = 1.0): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Inventions 08 00118 g010
Figure 11. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the low-frequency vibration modes for the coupled system (β = 1.0): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Figure 11. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the low-frequency vibration modes for the coupled system (β = 1.0): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Inventions 08 00118 g011
Figure 12. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the vibro-acoustic transfer characteristics for the coupled system (β = 1.0).
Figure 12. Effect of the strength variation of the nonlinear hardening stiffness characteristic on the vibro-acoustic transfer characteristics for the coupled system (β = 1.0).
Inventions 08 00118 g012
Figure 13. Effect of the length variation of the approximate linear working interval on the periodic bifurcation properties for the coupled system (α = 1.5): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Figure 13. Effect of the length variation of the approximate linear working interval on the periodic bifurcation properties for the coupled system (α = 1.5): (a) bifurcation diagram of the average vibration displacement; (b) bifurcation diagram of the average vibration velocity.
Inventions 08 00118 g013
Figure 14. Effect of the length variation of the approximate linear working interval on the low-frequency vibration modes for the coupled system (α = 1.5): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Figure 14. Effect of the length variation of the approximate linear working interval on the low-frequency vibration modes for the coupled system (α = 1.5): (a) motion phase diagram at the connection points of the upper vibration isolators with the unit A1; (b) motion phase diagram at the connection points of the lower vibration isolators with the raft frame; (c) motion phase diagram at the connection points of the lower vibration isolators with the cylindrical shell.
Inventions 08 00118 g014
Figure 15. Effect of the length variation of the approximate linear working interval on the vibro-acoustic transfer characteristics for the coupled system (α = 1.5).
Figure 15. Effect of the length variation of the approximate linear working interval on the vibro-acoustic transfer characteristics for the coupled system (α = 1.5).
Inventions 08 00118 g015
Table 1. Control parameter settings of variable stiffness vibration isolator.
Table 1. Control parameter settings of variable stiffness vibration isolator.
Control ParametersSymbolsUnitsValuesRemarks
Stiffness control parameters in the approximate linear working intervalKA11N/m2.0 × 105
KA21N/m2.7 × 105
KR1N/m8.6 × 105
Stiffness control parameters outside the linear working intervalKA12N/mα⋅2.0 × 105α ∈ [0.1, 3.0]
α < 1, softening characteristic
α > 1, hardening characteristic
KA22N/mα⋅2.7 × 105
KR2N/mα⋅8.6 × 105
Linear working interval length control parametersδ10mβ⋅6.125 × 10−4β ∈ [0.2, 5.0]
δ20mβ⋅6.352 × 10−4
δR0mβ⋅6.267 × 10−4
Damping coefficientsc1kg/s159.155
c2kg/s214.859
cRkg/s684.366
Table 2. System structure parameter settings.
Table 2. System structure parameter settings.
System Structure ParametersSymbolsUnitsValues
Characteristic parameters of the cylindrical shellElastic modulusEN/m22.1 × 1011
Densityρkg/m37.8 × 103
Poisson’s ratioμ 0.28
Dissipation factorξ 0.01
Length, radius, thicknessL, a, dm2, 0.4, 0.02
Structure parameters of the units and the raft frameMass of the unit A1m1kg50
Mass of the unit A2m2kg70
Mass of the raft frame RmRkg100
Moment of inertia of the unit A1J1x, J1zKg·m21, 0.5
Moment of inertia of the unit A2J2x, J2zKg·m21.5, 1
Moment of inertia of the raft frame RJRx, JRzKg·m28.5, 2.5
Linear complex stiffness of the vibration isolation supportsSupport stiffness of the unit A1k1N/m2.0 × 105 × (1 + 0.1 j)
Support stiffness of the unit A2k2N/m2.7 × 105 × (1 + 0.1 j)
Support stiffness of the raft frame RkRN/m8.6 × 105 × (1 + 0.1 j)
Installation position parametersCentroid positioning of the unit A1b1m0.25
Centroid positioning of the unit A2b2m0.25
Centroid positioning of the raft frame RLRm0.8
Support spacing of the unit A1h1x, h1zm0.2, 0.3
Support spacing of the unit A2h2x, h2zm0.2, 0.3
Support spacing of the raft frame RHx, Hzm0.35, 0.7
Table 3. The first 10 orders of poles and modal frequencies of the cylindrical shell—underwater acoustic field coupled subsystem.
Table 3. The first 10 orders of poles and modal frequencies of the cylindrical shell—underwater acoustic field coupled subsystem.
Orders iPoles sCi, sCi* *Modal Frequencies fCi (Hz)
1−0.366 ± 8.277 × 102 j131.73
2−0.574 ± 1.219 × 103 j194.95
3−0.660 ± 1.672 × 103 j266.24
4−0.924 ± 1.906 × 103 j303.33
5−1.246 ± 2.302 × 103 j366.30
6−1.753 ± 2.364 × 103 j376.21
7−1.307 ± 2.621 × 103 j417.11
8−1.586 ± 3.065 × 103 j487.86
9−3.624 ± 3.304 × 103 j525.88
10−2.863 × 103 ± 1.973 × 103 j553.34
* The poles sCi and sCi* are a pair of conjugate complex numbers.
Table 4. The values of the modal parameters ϕ ¯ i j at the support connection points.
Table 4. The values of the modal parameters ϕ ¯ i j at the support connection points.
Orders ij = 1j = 2j = 3j = 4
i = 10.475 − 4.630 j0.686 − 0.703 j0.104 − 0.117 j0.163 − 0.172 j
i = 20.462 − 0.429 j0.601 − 0.647 j−0.317 + 0.312 j−0.448 + 0.469 j
i = 30.436 − 0.208 j0.271 − 0.328 j0.068 − 0.226 j0.214 − 0.335 j
i = 40.510 − 0.529 j−0.274 + 0.240 j−0.398 + 0.384 j0.170 − 0.169 j
i = 50.367 − 0.352 j0.402 − 0.441 j−0.347 + 0.339 j−0.453 + 0.479 j
i = 60.159 − 0.407 j−0.342 + 0.117 j0.412 − 0.031 j0.006 + 0.032 j
i = 70.473 − 0.475 j−0.218 + 0.218 j−0.464 + 0.464 j0.215 − 0.215 j
i = 80.357 − 0.356 j−0.317 + 0.318 j−0.261 + 0.261 j0.233 − 0.234 j
i = 90.470 − 0.434 j−0.254 + 0.271 j−0.361 + 0.357 j0.248 − 0.270 j
i = 100.716 − 0.762 j0.082 − 0.105 j−0.244 + 0.233 j−0.053 + 0.012 j
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, L.; Huo, R. Effect of Support Stiffness Nonlinearity on the Low-Frequency Vibro-Acoustic Characteristics for a Mechanical Equipment—Floating Raft—Underwater Cylindrical Shell Coupled System. Inventions 2023, 8, 118. https://doi.org/10.3390/inventions8050118

AMA Style

Wang L, Huo R. Effect of Support Stiffness Nonlinearity on the Low-Frequency Vibro-Acoustic Characteristics for a Mechanical Equipment—Floating Raft—Underwater Cylindrical Shell Coupled System. Inventions. 2023; 8(5):118. https://doi.org/10.3390/inventions8050118

Chicago/Turabian Style

Wang, Likang, and Rui Huo. 2023. "Effect of Support Stiffness Nonlinearity on the Low-Frequency Vibro-Acoustic Characteristics for a Mechanical Equipment—Floating Raft—Underwater Cylindrical Shell Coupled System" Inventions 8, no. 5: 118. https://doi.org/10.3390/inventions8050118

APA Style

Wang, L., & Huo, R. (2023). Effect of Support Stiffness Nonlinearity on the Low-Frequency Vibro-Acoustic Characteristics for a Mechanical Equipment—Floating Raft—Underwater Cylindrical Shell Coupled System. Inventions, 8(5), 118. https://doi.org/10.3390/inventions8050118

Article Metrics

Back to TopTop