# Nonlinear Dynamics of Marine Rotor System Coupled with Air Bag-Floating Raft Subjected to the Basement Excitations in Lateral Directions

Authors: Wen Zhao; Ming Li; Yuanbo Liu
Journal: Shock and Vibration
 (2020)

---

## Abstract

This paper presents an investigation into the nonlinear dynamic behaviors of the mechanical isolation system coupled with air-bag and floating-raft subject to basement excitation in lateral directions. First, the coupling effects between the excitation source and isolation system are considered. Also, the mechanical isolation model under basic excitation and its motion equation are deduced, and then the dynamic responses are mainly investigated by using the techniques of displacement response, frequency spectrum, rotor orbit, Poincaré maps, and the bifurcation diagram. Last, the bifurcations of the mechanical isolation system with different parameters are analyzed through numerical methods, especially the effect of excitation frequency and amplitude. The result predicts that period-5 is mainly performed, with the increase of rotor speed, and the system moves into quasi-bifurcation. However, the system stays in chaos state at high rotor speed, and the vibration amplitude rises rapidly until against bearing bush. Furthermore, the effects of basement excitation on the mechanical isolation system are mainly concentrated on the stage of lower rotor speed, but with the increasing speed, the effects become weak and at the same time the vibration amplitude reduces significantly. The points projected on the Poincaré section are five, three, or two solitary attractors, in which the system stays in periodic motion. Above all, the dynamic characteristics can provide the theoretic supporting for the dynamic, vibration control and its parametric optimization of the marine mechanical isolation system coupled with air-bag and floating-raft.


---

## Body

## 1. Introduction

The external environment is relatively complex because of work requirements on the rotating machinery. Usually, base excitation of rotating machinery could occur when its surroundings are equipped with a plurality of mechanical power equipment such as vessels and submarines on ground that is prone to seismic motion or other similar environments. Foundation excitations are the most common source of harmonic excitation of rotating machinery because any form of excitation can be converted into harmonic excitation by fast Fourier transformation (FFT). A rotating machine under motion base could be excited at its base by another harmonic excitation that in turn could arise from, for example, vibration from basement motion on the ground or ship. The reduction of undesired high vibration amplitudes is of paramount importance. It is essential for a designer to have at their disposal dynamic models to predict the vibrational behavior of rotating machines operating under such multifrequency excitations at the design stage to avoid expensive design modifications or in-field catastrophic failure.There are many researches on the basic excitation. Wilson and Brebbia investigated the dynamic behavior of steel foundations for turbo-alternators by the finite element displacement method [1]; Gasch studied the vibration of large turbo-rotors in fluid-film bearings on an elastic foundation by the finite element analysis [2]. Hori reviewed 19 works mainly on the foundation excitation effect on the dynamics/stability of rotor-fluid film bearings systems under basement excitation [3]. Rajan and Davies considered the random primary response of a duffing oscillator subject to narrowband excitation by using the method of multiple scales and stochastic averaging [4]. Nayfeh and Serhan analyzed the stationary mean and mean square responses and its local stability of a Duffing–Rayleigh oscillator excited by the sum of a deterministic harmonic component and a random component by using a second-order closure method [5]; Chen studied the nonlinear dynamical behavior of a single-axis rate gyro with harmonic excitation mounted on a space vehicle [6]; American scientists researched a method for determining the relationship between the structural type of a manipulator and its susceptibility to motion-induced vibrational excitation [7]; Feng and Hu developed a set of nonlinear differential equations by using Kane’s method for the planar oscillation of slender beams subject to a parametric excitation of the base movement [8, 9]. Fawzi developed a formulation for the dynamic analysis of the rigid rotor subject to base excitations plus mass imbalance. Han and Chu acquired the parametric instability regions induced by the periodic base angular motions based on the discrete state transition matrix method [10, 11]. Dakel et al. investigated the dynamics of two different rotor configurations under various support motion through the numerical method [12]. Saeed and Kamel investigated the vibration control of a horizontally suspended Jeffcott-rotor system. A nonlinear restoring force and the rotor weight are considered in the system model [13]. Lei et al. focused on the nonlinear responses of a cracked rotor-ball bearing system caused by aircraft flight maneuvers [14]. Saeed and El-Ganaini deal with the vibration analysis of a Jeffcott-rotor system under various supports. Finally, a comparison with the papers that published previously is included [15–18]. Chen and Wang developed a general model for the flexible rotor system subjected to time-varying base excitation and studied the direct effects of angular base motion on the dynamic behaviors of a simple rotor [19].However, it is worth mentioned that researchers’ works mainly concentrated on the effect of base motions while the mechanical isolation system coupled with air bag-floating raft under base motions is not referred. In addition, the excitation resource from their researches is always in single direction, and the definition of base excitation was not accurate enough to describe the real movement of the base.In this paper, an investigation into the nonlinear dynamic behaviors of the mechanical isolation system coupled with air-bag and floating-raft under basement excitation in lateral directions is presented. First, the coupling effects between the excitation source and isolation system are considered, and the excitation on the longitudinal and transverse directions is also taken into account. Based on this, the mechanical isolation model subjected to basic excitation and its motion equation are developed, and then its dynamic responses are mainly investigated by using the techniques of displacement response, frequency spectrum, rotor orbit, Poincaré maps, and the bifurcation diagram. Last, the bifurcations of the mechanical isolation system with different parameters are analyzed through numerical methods, especially the effect of excitation frequency and amplitude. In particular, this work aims to present design guidelines for predicting dynamics of the mechanical isolation system under various basement excitation.


## 2. Dynamic Model under Basement Excitation

Due to the low frequency vibration of ship hull and other power machinery equipment, mainly motor or diesel engine and other power equipment, the form of its vibrations is mainly harmonic vibrations [20].


### 2.1. Dynamic Model

Figure1 depicts the rotating mechanical isolation system coupled with the air bag-floating raft device subjected to basement excitation, in which the floating-raft and the rotor are considered as a rigid body and the air-bag is regarded as a spring with vertical and horizontal deformations. In general, the air-bag and floating-raft isolation device is a kind of a mixed complex dynamic isolation system which includes air-bags, floating-raft, and rotor-bearings. Particularly, all the system should be fixed on the ground such as the ship and vessel that can be prone to seismic motion.Figure 1
Dynamic motion model.It is well known that any form of excitation can be converted into harmonic excitation by the FFT technique; thus, basement motions in this paper are considered as harmonic motionS(t) driven by the angular frequency Ω1 and Ω2 in the longitudinal and transverse directions, respectively, and by considering the coupling effects between the excitation source and isolation system, the mechanical model of rotating machinery under base excitation is established.Also, based on the dynamics theory, the differential equations of motion [21] under base excitation motion are derived as follows:(1)m1x¨1=Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2+m1eΩ2cosΩt+m1g,m1y¨1=Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2+m1eΩ2sinΩt,m2x¨2=−Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2−k1x2−s1−d1x˙2−s˙1+m2g,m2y¨2=−Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2−k2y2−s2−d2y˙2−s˙2,where m1, m2, e, k, d, Ω, Fx, and Fy are the mass of disc, the mass of floating-raft, eccentricity, stiffness, damping, rotor speed, nonlinear oil film in x and y directions, respectively; s1(t) and s2(t) are the longitudinal and transverse basement excitation/driven-base excitation/foundation excitation, respectively, and its motion parametric settings can be expressed as s1t=ΦasinΩ1t and s2t=ΦbsinΩ2t;the parameters Ω1 and Ω2 are the angular frequency in the longitudinal and transverse directions, respectively; and Φa and Φb are the excitation amplitudes of longitudinal and transverse directions, respectively.


### 2.2. Oil Film Force

Figure2 illustrates the schematic diagram of journal bearing, in which Fr and Fτ are the nonlinear oil forces in the radial and tangential directions, respectively. In this article, the long bearing theory is introduced; accordingly, the Reynolds equation for the oil film pressure p yields [22, 23](2)1Rb2∂∂φh3u∂p∂φ=6Ω−2θ˙∂h∂φ+12e˙cosφ,where h=c+ecosθ=c1+εcosθ.Figure 2
Schematic diagram of journal bearing.After considering the semi-Sommerfeld conditionθ=0, p=0; θ=π, p=0, the oil film forces in radial and tangent directions are obtained as follows:(3)Fr=6uBRRc2Ω−2dθdtE1+2dεdtE2,Fτ=6uBRRc2Ω−2dθdtE3+2dεdtE4,where(4)E1=2α21−α22+α2,E2=11−α23/2π2−8π2+α2,E3=πα1−α21/22+α2,E4=2α1−α22+α2.The oil film forces in equation (3) can be expressed as(5)Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2=−Frcosφ−Ftsinφ,Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2=−Frsinφ+Ftcosφ,where Fx and Fy are the nonlinear oil film forces components of journal bearing in the x and y direction, respectively.


### 2.3. Nondimensional Motion Equation

For the simplification and wide applicability of the following analysis, the nondimensional variables are obtained on the basis of the bearing clearance, weight of disc, and time as shown in Table1.Table 1
The dimensionless expressions of the parameters.
Parameters
Expression
α
α = e/c
τ
τ
=
Ω
t
σ
¯
σ
¯
=
6
u
B
R
3
/
c
2
n
n = m2/m1
ω
ω
=
Ω
c
/
g
D
D
=
d
c
Ω
/
m
1
g
x
¯
i
x
¯
i
=
x
i
/
c
y
¯
i
y
¯
i
=
y
i
/
c
f
j
f
j
=
F
j
/
m
1
g
σ
σ
=
σ
¯
/
m
1
g
c
λ
λ
=
B
/
2
R
K
K
=
kc
/
m
1
g
η
1
η
1
=
Ω
1
/
Ω
η
2
η
2
=
Ω
2
/
Ω
ϕ
a
ϕ
a
=
Φ
a
/
c
ϕ
b
ϕ
b
=
Φ
b
/
c
f
x
f
x
=
−
f
r
cos
φ
−
f
t
sin
φ
f
y
f
x
=
−
f
r
cos
φ
−
f
t
sin
φBy the way, denotedx/dt=x˙anddx/dτ=x¯′. The meaning of the parameters in the above equations is shown in nomenclature, and equation (3) can be put into a convenient dimensionless form as follows:(6)fr=σω1−2φ′E1+2e′E2,fτ=σω1−2φ′E3+2e′E4.Thereby, we obtain the following set of first-order and second-order differential relation:(7)dxdt=dxdτdτdt=Ωdxdτ,d2xdt2=ddtΩdxdτ=ddτdτdtΩdxdτ=Ω2d2xdτ2.We rewritedx/dtandd2x/dt2 as follows:(8)dxdt=Ωx¯′,d2xdt=Ωx¯″.In the same way, we get(9)x¨1=Ω2x¯1″,y¨1=Ω2y¯1″,x¨2=Ω2x¯2″,y¨2=Ω2y¯2″.Substituting equations (3) to (9) into equation (1), we get the nondimensional formulation of motion:(10)x¯1″=1ω2+fxω2+αcosτ,y¯1″=fyω2+αsinτ,x¯2″=1ω2−fxnω2−K1nω2x¯2−D1nω2x¯2′+K1nω2+D1nω2ϕaη1cosη1τ,y¯2″=−fynω2−K2nω2y¯2−D2nω2y¯2′+K2nω2ϕbsinη2τ+D2nω2ϕbη2cosη2τ.Equation (10) is a two order nonlinear differential equation group controlled by 11 varied parameters such as α, ω, n, η1, η2, ϕa, ϕb, D1, D2, K1, and K2. The motion styles of system under the unbalanced excitation and basement excitation are complex, where η1 and η2 are the dimensionless vibration frequency in the vertical and horizontal direction and ϕa and ϕb are the dimensionless excitation amplitude, respectively. Hence, we introduce s1τ and s2τ as the unbalanced excitation and basement excitation, and we have(11)s1τ=ϕasinη1τ,s˙1τ=ϕaη1cosη1τ,s2τ=ϕbsinη2τ,s˙2τ=ϕbη2cosη2τ.Plugging equation (11) into equation (10), we get(12)Z˙1τ=x¯˙1τ=Z2τ,Z˙2τ=x¯¨1τ=1ω2+fxω2+αcosτ,Z˙3τ=y¯˙1τ=Z4τ,Z˙4τ=y¯¨1τ=fyω2+αsinτ,Z˙5τ=x¯˙2τ=Z6τ,Z˙6τ=x¯¨2τ=1ω2−fxnω2−K1nω2Z5τ−D1nω2Z6τ+K1nω2s1τ+D1nω2s˙1τ,Z˙7τ=y¯˙2τ=Z8τ,Z˙8τ=y¯¨2τ=−fynω2−K2nω2Z7τ−D2nω2Z8τ+K2nω2s2τ+D2nω2s˙2τ.Equation (10) belongs to high-dimensional, strong constrained, and nonlinear nonautonomous problems with various parameters; generally, in order to facilitate calculation and solution for convenience. Equation (10) could be discussed on the state space, and hence the four second-order equations can be easily converted into eight first-order ones and the equation of motion governing the behavior of the mechanical isolation system can be written as equation (12).


## 2.1. Dynamic Model

Figure1 depicts the rotating mechanical isolation system coupled with the air bag-floating raft device subjected to basement excitation, in which the floating-raft and the rotor are considered as a rigid body and the air-bag is regarded as a spring with vertical and horizontal deformations. In general, the air-bag and floating-raft isolation device is a kind of a mixed complex dynamic isolation system which includes air-bags, floating-raft, and rotor-bearings. Particularly, all the system should be fixed on the ground such as the ship and vessel that can be prone to seismic motion.Figure 1
Dynamic motion model.It is well known that any form of excitation can be converted into harmonic excitation by the FFT technique; thus, basement motions in this paper are considered as harmonic motionS(t) driven by the angular frequency Ω1 and Ω2 in the longitudinal and transverse directions, respectively, and by considering the coupling effects between the excitation source and isolation system, the mechanical model of rotating machinery under base excitation is established.Also, based on the dynamics theory, the differential equations of motion [21] under base excitation motion are derived as follows:(1)m1x¨1=Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2+m1eΩ2cosΩt+m1g,m1y¨1=Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2+m1eΩ2sinΩt,m2x¨2=−Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2−k1x2−s1−d1x˙2−s˙1+m2g,m2y¨2=−Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2−k2y2−s2−d2y˙2−s˙2,where m1, m2, e, k, d, Ω, Fx, and Fy are the mass of disc, the mass of floating-raft, eccentricity, stiffness, damping, rotor speed, nonlinear oil film in x and y directions, respectively; s1(t) and s2(t) are the longitudinal and transverse basement excitation/driven-base excitation/foundation excitation, respectively, and its motion parametric settings can be expressed as s1t=ΦasinΩ1t and s2t=ΦbsinΩ2t;the parameters Ω1 and Ω2 are the angular frequency in the longitudinal and transverse directions, respectively; and Φa and Φb are the excitation amplitudes of longitudinal and transverse directions, respectively.


## 2.2. Oil Film Force

Figure2 illustrates the schematic diagram of journal bearing, in which Fr and Fτ are the nonlinear oil forces in the radial and tangential directions, respectively. In this article, the long bearing theory is introduced; accordingly, the Reynolds equation for the oil film pressure p yields [22, 23](2)1Rb2∂∂φh3u∂p∂φ=6Ω−2θ˙∂h∂φ+12e˙cosφ,where h=c+ecosθ=c1+εcosθ.Figure 2
Schematic diagram of journal bearing.After considering the semi-Sommerfeld conditionθ=0, p=0; θ=π, p=0, the oil film forces in radial and tangent directions are obtained as follows:(3)Fr=6uBRRc2Ω−2dθdtE1+2dεdtE2,Fτ=6uBRRc2Ω−2dθdtE3+2dεdtE4,where(4)E1=2α21−α22+α2,E2=11−α23/2π2−8π2+α2,E3=πα1−α21/22+α2,E4=2α1−α22+α2.The oil film forces in equation (3) can be expressed as(5)Fxx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2=−Frcosφ−Ftsinφ,Fyx1−x2,y1−y2,x˙1−x˙2,y˙1−y˙2=−Frsinφ+Ftcosφ,where Fx and Fy are the nonlinear oil film forces components of journal bearing in the x and y direction, respectively.


## 2.3. Nondimensional Motion Equation

For the simplification and wide applicability of the following analysis, the nondimensional variables are obtained on the basis of the bearing clearance, weight of disc, and time as shown in Table1.Table 1
The dimensionless expressions of the parameters.
Parameters
Expression
α
α = e/c
τ
τ
=
Ω
t
σ
¯
σ
¯
=
6
u
B
R
3
/
c
2
n
n = m2/m1
ω
ω
=
Ω
c
/
g
D
D
=
d
c
Ω
/
m
1
g
x
¯
i
x
¯
i
=
x
i
/
c
y
¯
i
y
¯
i
=
y
i
/
c
f
j
f
j
=
F
j
/
m
1
g
σ
σ
=
σ
¯
/
m
1
g
c
λ
λ
=
B
/
2
R
K
K
=
kc
/
m
1
g
η
1
η
1
=
Ω
1
/
Ω
η
2
η
2
=
Ω
2
/
Ω
ϕ
a
ϕ
a
=
Φ
a
/
c
ϕ
b
ϕ
b
=
Φ
b
/
c
f
x
f
x
=
−
f
r
cos
φ
−
f
t
sin
φ
f
y
f
x
=
−
f
r
cos
φ
−
f
t
sin
φBy the way, denotedx/dt=x˙anddx/dτ=x¯′. The meaning of the parameters in the above equations is shown in nomenclature, and equation (3) can be put into a convenient dimensionless form as follows:(6)fr=σω1−2φ′E1+2e′E2,fτ=σω1−2φ′E3+2e′E4.Thereby, we obtain the following set of first-order and second-order differential relation:(7)dxdt=dxdτdτdt=Ωdxdτ,d2xdt2=ddtΩdxdτ=ddτdτdtΩdxdτ=Ω2d2xdτ2.We rewritedx/dtandd2x/dt2 as follows:(8)dxdt=Ωx¯′,d2xdt=Ωx¯″.In the same way, we get(9)x¨1=Ω2x¯1″,y¨1=Ω2y¯1″,x¨2=Ω2x¯2″,y¨2=Ω2y¯2″.Substituting equations (3) to (9) into equation (1), we get the nondimensional formulation of motion:(10)x¯1″=1ω2+fxω2+αcosτ,y¯1″=fyω2+αsinτ,x¯2″=1ω2−fxnω2−K1nω2x¯2−D1nω2x¯2′+K1nω2+D1nω2ϕaη1cosη1τ,y¯2″=−fynω2−K2nω2y¯2−D2nω2y¯2′+K2nω2ϕbsinη2τ+D2nω2ϕbη2cosη2τ.Equation (10) is a two order nonlinear differential equation group controlled by 11 varied parameters such as α, ω, n, η1, η2, ϕa, ϕb, D1, D2, K1, and K2. The motion styles of system under the unbalanced excitation and basement excitation are complex, where η1 and η2 are the dimensionless vibration frequency in the vertical and horizontal direction and ϕa and ϕb are the dimensionless excitation amplitude, respectively. Hence, we introduce s1τ and s2τ as the unbalanced excitation and basement excitation, and we have(11)s1τ=ϕasinη1τ,s˙1τ=ϕaη1cosη1τ,s2τ=ϕbsinη2τ,s˙2τ=ϕbη2cosη2τ.Plugging equation (11) into equation (10), we get(12)Z˙1τ=x¯˙1τ=Z2τ,Z˙2τ=x¯¨1τ=1ω2+fxω2+αcosτ,Z˙3τ=y¯˙1τ=Z4τ,Z˙4τ=y¯¨1τ=fyω2+αsinτ,Z˙5τ=x¯˙2τ=Z6τ,Z˙6τ=x¯¨2τ=1ω2−fxnω2−K1nω2Z5τ−D1nω2Z6τ+K1nω2s1τ+D1nω2s˙1τ,Z˙7τ=y¯˙2τ=Z8τ,Z˙8τ=y¯¨2τ=−fynω2−K2nω2Z7τ−D2nω2Z8τ+K2nω2s2τ+D2nω2s˙2τ.Equation (10) belongs to high-dimensional, strong constrained, and nonlinear nonautonomous problems with various parameters; generally, in order to facilitate calculation and solution for convenience. Equation (10) could be discussed on the state space, and hence the four second-order equations can be easily converted into eight first-order ones and the equation of motion governing the behavior of the mechanical isolation system can be written as equation (12).


## 3. Nonlinear Dynamic Analysis

For the strongly nonlinear differential equation (12), no analytical method exists, which is solved by applying the Runge–Kutta integration routine with a variable step in the present work. The calculation parameters of the system are as follows [24, 25]:n = 10, α = 0.1, σ = 3, K1 = K2 = 12, D1 = D2 = 0.1, x1=0.2, x˙1=0.0, y1=0.1, y˙1=0.0, x2=0.1, x˙2=0.0, y2=0.1, and y˙2=0.0.Equation (12) is eight one-order nonlinear differential equations controlled by 11 varied parameters. Solving the analytical solution or the approximate analytical solution of the system is much more difficult to address than we imagine. Thus, numerical methods are widely used to address the nonlinear problems. In order to simulate easily, the initial conditions are given as following: η1 = η2 = η = 0.125; more importantly, the dynamic characteristic of the basement excitation isolation system under given conditions can be analyzed in detail.


### 3.1. Bifurcation Diagram

Because equation (12) is a nonlinear nonautonomous system with multidimensional rotor-bearing system, which contains the nonlinear oil film, nonlinear functions, and excitation function about sinϕ and cosϕ, consequently, the analytic solution or the approximate one cannot be obtained easily. A theoretical analysis has been conducted for the mechanical isolation system. Toward this, (12) is integrated over a long period of time for bifurcation of initial condition by using the Rung–Kutta technique to be introduced to calculate dynamic response of the rotating mechanical isolation system coupled with air bag-floating raft under basement excitation.In the non-linear analysis, for the purpose, the bifurcation diagram is constructed via numerical calculations. Construction of a bifurcation diagram is a standard approach used to analyze various nonlinear systems. In the bifurcation diagram, the dynamical behaviors may be viewed globally over a range of parameter values and allow us to compare simultaneously different types of motion. Therefore, the bifurcation diagram provides a summary of essential dynamics and is a useful tool for acquiring an overview.Figure3 illustrates the bifurcation diagram based on the above initial condition and shows the theoretically calculated bifurcation diagram of the mechanical isolation system with respect to the key motion parameter ω, in which ω varied from 0.60 to 3.50. It can be summarized from the bifurcation diagram (a) that the main influences by basic excitation on the mechanical isolation system are mainly concentrated in the low rotating speed stage, where the system amplitude is large in comparation with the prior research [26] and the response shows a periodic motion. With the rapid increments of rotor speed, the system vibration decreases gradually in the vertical direction. Figure 3(b) shows bifurcation of the same system for the displacement response y1 with respect to the bifurcation parameter ω. The results indicate it has the same dynamics at lower speed, but with the increase of rotor speed, the vibration starts to rise and the dynamic performance demonstrates period-5, and when the system stays at higher speed, strong vibration and chaos will appear.Figure 3
Dynamic bifurcation diagram and its LLE: (a)x; (b) y.
(a)
(b)


### 3.2. Dynamic Steady-State Response

Figures4 and 5 illustrate the steady response and its frequency spectrum, rotor obit, and its Poincaré map when the rotor speed is ω = 1.20 and 2.25, respectively. The numerical results show that the dynamic response of the isolation system contains the rotor frequency and the basement excitation frequency. The basement excitation affects the external force of the motion differential equation; therefore, the trajectory of the system is more complex and the nonlinear factor on the mechanical isolation system has weak effects; the repeatability of the rotor orbit is better expressed clearly, and the same time it is worth mentioning here that points on the Poincaré section reflect five isolated points, which means the dynamic characteristics of the system shows period-5 motion.Figure 4
ω = 1.210: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 5
ω = 2.25: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure6 shows the oscillations at ω = 3.412, in which the nonlinear dynamic behaviors are the quasi-period motion and the displacement response (Figure 6(a)) highlights steady vibration. The 1/2 frequency components mainly dominate on the frequency spectrum (Figure 6(a)), the steady-state response of this system is confined to an annular region, and the rotor makes complex movements as shown in (Figure 6(c)). Points on the Poincaré maps are an approximate closed curve, the dynamic characteristic of the quasi-period motion at this rotor speed.Figure 6
ω = 3.412: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figures7–9 illustrate the displacement response and its frequency spectrum, rotor obit, and its Poincaré map when ω = 3.35, 3.50, and 3.58, respectively. The dynamic behaviors of the system mainly show chaotic motion, in which the displacement response shows asymmetry and a continuous frequency is not expressed in the frequency spectrum clearly, but trajectory is a mess network structure and its Poincaré section shows some break, scattered points that are not closed.Figure 7
ω = 3.35: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 8
ω = 3.50: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 9
ω = 3.58: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)


## 3.1. Bifurcation Diagram

Because equation (12) is a nonlinear nonautonomous system with multidimensional rotor-bearing system, which contains the nonlinear oil film, nonlinear functions, and excitation function about sinϕ and cosϕ, consequently, the analytic solution or the approximate one cannot be obtained easily. A theoretical analysis has been conducted for the mechanical isolation system. Toward this, (12) is integrated over a long period of time for bifurcation of initial condition by using the Rung–Kutta technique to be introduced to calculate dynamic response of the rotating mechanical isolation system coupled with air bag-floating raft under basement excitation.In the non-linear analysis, for the purpose, the bifurcation diagram is constructed via numerical calculations. Construction of a bifurcation diagram is a standard approach used to analyze various nonlinear systems. In the bifurcation diagram, the dynamical behaviors may be viewed globally over a range of parameter values and allow us to compare simultaneously different types of motion. Therefore, the bifurcation diagram provides a summary of essential dynamics and is a useful tool for acquiring an overview.Figure3 illustrates the bifurcation diagram based on the above initial condition and shows the theoretically calculated bifurcation diagram of the mechanical isolation system with respect to the key motion parameter ω, in which ω varied from 0.60 to 3.50. It can be summarized from the bifurcation diagram (a) that the main influences by basic excitation on the mechanical isolation system are mainly concentrated in the low rotating speed stage, where the system amplitude is large in comparation with the prior research [26] and the response shows a periodic motion. With the rapid increments of rotor speed, the system vibration decreases gradually in the vertical direction. Figure 3(b) shows bifurcation of the same system for the displacement response y1 with respect to the bifurcation parameter ω. The results indicate it has the same dynamics at lower speed, but with the increase of rotor speed, the vibration starts to rise and the dynamic performance demonstrates period-5, and when the system stays at higher speed, strong vibration and chaos will appear.Figure 3
Dynamic bifurcation diagram and its LLE: (a)x; (b) y.
(a)
(b)


## 3.2. Dynamic Steady-State Response

Figures4 and 5 illustrate the steady response and its frequency spectrum, rotor obit, and its Poincaré map when the rotor speed is ω = 1.20 and 2.25, respectively. The numerical results show that the dynamic response of the isolation system contains the rotor frequency and the basement excitation frequency. The basement excitation affects the external force of the motion differential equation; therefore, the trajectory of the system is more complex and the nonlinear factor on the mechanical isolation system has weak effects; the repeatability of the rotor orbit is better expressed clearly, and the same time it is worth mentioning here that points on the Poincaré section reflect five isolated points, which means the dynamic characteristics of the system shows period-5 motion.Figure 4
ω = 1.210: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 5
ω = 2.25: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure6 shows the oscillations at ω = 3.412, in which the nonlinear dynamic behaviors are the quasi-period motion and the displacement response (Figure 6(a)) highlights steady vibration. The 1/2 frequency components mainly dominate on the frequency spectrum (Figure 6(a)), the steady-state response of this system is confined to an annular region, and the rotor makes complex movements as shown in (Figure 6(c)). Points on the Poincaré maps are an approximate closed curve, the dynamic characteristic of the quasi-period motion at this rotor speed.Figure 6
ω = 3.412: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figures7–9 illustrate the displacement response and its frequency spectrum, rotor obit, and its Poincaré map when ω = 3.35, 3.50, and 3.58, respectively. The dynamic behaviors of the system mainly show chaotic motion, in which the displacement response shows asymmetry and a continuous frequency is not expressed in the frequency spectrum clearly, but trajectory is a mess network structure and its Poincaré section shows some break, scattered points that are not closed.Figure 7
ω = 3.35: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 8
ω = 3.50: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 9
ω = 3.58: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)


## 4. Parametric Analysis

### 4.1. The Effect of Excitation Frequency

Figure10 illustrates the bifurcation diagram of floating-raft forced by the displacement response x1 and y1 with respect to the ratio of excitation frequency η, varying from 0.1 to 0.5 in certain parameters for ω = 1.20 based on the initial condition as σ=3, α=0.01, n=30, K1=100, K2=40, D1=D2=0.1, λ=0.2, and ϕ=0.0175. It is worth mentioning here that the external force of the motion differential equation is influenced by the basement excitation; thus, the trajectory is no longer a circle or ellipse, but a complex shape with the η increasing. The vertical displacement x1 of the mechanical isolation system at low frequency shows the high amplitude due to the effect of basic excitation frequency, but the amplitude starts to drop gradually with the increase of frequency (however the y2 had a opposite trend); throughout the whole period, the dynamic responses are dominated by period-5, quasi-period, period-2, period-3, and chaos.Figure 10
Dynamic bifurcation diagram: (a)x1; (b) y1.
(a)
(b)Figures11–15 illustrate the displacement response and its frequency spectrum, rotor obit, and its Poincaré map at ω = 1.20 when η = 0.125, 0.167, 0.25, 0.5, and 0.65, respectively.Figure 11
ω = 1.20 and η = 0.125: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 12
ω = 1.20 and η = 0.167: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 13
ω = 1.20 and η = 0.25: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 14
ω = 1.20 and η = 0.27: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 15
ω = 1.20 and η = 0.43: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)From the spectrum, it can be found thatf/2 in the frequency response component is mainly dominated as shown. At this point, the frequency whirl vibration of the system is obvious; meanwhile, the unbalance magnitude of the system plays a dominant role during the whole process of operation.It can be seen from the Poincaré maps that five isolated points project on the Poincaré section whenη = 0.125, which presents period-5 motion. When η = 0.167 and 0.25, the dynamic characters of the system demonstrate period-3 motion because of the three unrelated and isolated attractors shown and the steady response have a steady beat vibration and trajectories are simple. However, with the increase of η, the trajectories start to become complex; furthermore, an unclosed curve is shown in Figure 16, and the Poincaré section also proves a chaos motion.Figure 16
ω = 1.20 and η = 0.50: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Generally, after periodic motion at single frequency disturbing, the system will move into steady state, and its trajectories will acquire a disjoint closed form. However, in the system under the basement external stimuliη = 0.5 (Figure 17), the phase plane will extend, unbalanced force plays a dominant role, the closure of “8” trajectory appears as shown in Figure 17, and its motion trajectory is periodic; two isolated points projected on the Poincaré maps create attractors, that is to say that the system move in period-two motion.Figure 17
ω = 1.20 and η = 0.650: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)


### 4.2. The Effect of Excitation Amplitude

Figure18 illustrates the bifurcation diagram at ω = 1.20 with respect to the ratio of excitation amplitude ϕ, in which ϕ varied from 0.01 to 0.10. The basement excitation affects the external force of the motion differential equation; therefore, the trajectory of the system is more complex. The mechanical isolation system is weakly affected by the excitation amplitude at minor excitation amplitude. However, its vibration amplitude will increase with the excitation amplitude. The dynamic behavior of the system stayed at state of period-5 throughout the whole period.Figure 18
Dynamic bifurcation diagram: (a)x1; (b) y1.
(a)
(b)Figures19–22 illustrate the dynamic steady response at ϕ = 0.020, 0.055, 0.075, and 0.100, respectively. The numerical simulation results indicate that the region value of system vibration, with the increase of excitation amplitude, grows rapidly and vibration amplitude also rises. The beat vibration and 1f are mainly shown on displacement response and frequency, which is the interaction of both system’s unbalance magnitude and basement excitation; thus, oil film oscillation appears. Furthermore, the other frequency components cannot be easily actuated under these specific values compared with the first mode frequency which indicates that the unbalance force and nonlinear oil force of the system play a dominant role during the whole process of operation.Figure 19
ω = 2.20 and ϕ = 0.02: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 20
ω = 2.20 and ϕ = 0.055: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 21
ω = 2.20 and ϕ = 0.075: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 22
ω = 2.20 and ϕ = 0.100: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Due to the basement excitation, the external force of equations is affected, as a result the rotor orbit has changed and clearly shows several single-closed curves. Besides, we note that the trajectory will expand outward iny direction as the excitation amplitude increases gradually when we compare the trajectories under various excitation values. More importantly, the points on the Poincaré map could be numbered and five isolated simple attractors appear apparently on its section, which demonstrates incontestably that the system moves into period-5 motion.


## 4.1. The Effect of Excitation Frequency

Figure10 illustrates the bifurcation diagram of floating-raft forced by the displacement response x1 and y1 with respect to the ratio of excitation frequency η, varying from 0.1 to 0.5 in certain parameters for ω = 1.20 based on the initial condition as σ=3, α=0.01, n=30, K1=100, K2=40, D1=D2=0.1, λ=0.2, and ϕ=0.0175. It is worth mentioning here that the external force of the motion differential equation is influenced by the basement excitation; thus, the trajectory is no longer a circle or ellipse, but a complex shape with the η increasing. The vertical displacement x1 of the mechanical isolation system at low frequency shows the high amplitude due to the effect of basic excitation frequency, but the amplitude starts to drop gradually with the increase of frequency (however the y2 had a opposite trend); throughout the whole period, the dynamic responses are dominated by period-5, quasi-period, period-2, period-3, and chaos.Figure 10
Dynamic bifurcation diagram: (a)x1; (b) y1.
(a)
(b)Figures11–15 illustrate the displacement response and its frequency spectrum, rotor obit, and its Poincaré map at ω = 1.20 when η = 0.125, 0.167, 0.25, 0.5, and 0.65, respectively.Figure 11
ω = 1.20 and η = 0.125: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 12
ω = 1.20 and η = 0.167: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 13
ω = 1.20 and η = 0.25: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 14
ω = 1.20 and η = 0.27: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 15
ω = 1.20 and η = 0.43: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)From the spectrum, it can be found thatf/2 in the frequency response component is mainly dominated as shown. At this point, the frequency whirl vibration of the system is obvious; meanwhile, the unbalance magnitude of the system plays a dominant role during the whole process of operation.It can be seen from the Poincaré maps that five isolated points project on the Poincaré section whenη = 0.125, which presents period-5 motion. When η = 0.167 and 0.25, the dynamic characters of the system demonstrate period-3 motion because of the three unrelated and isolated attractors shown and the steady response have a steady beat vibration and trajectories are simple. However, with the increase of η, the trajectories start to become complex; furthermore, an unclosed curve is shown in Figure 16, and the Poincaré section also proves a chaos motion.Figure 16
ω = 1.20 and η = 0.50: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Generally, after periodic motion at single frequency disturbing, the system will move into steady state, and its trajectories will acquire a disjoint closed form. However, in the system under the basement external stimuliη = 0.5 (Figure 17), the phase plane will extend, unbalanced force plays a dominant role, the closure of “8” trajectory appears as shown in Figure 17, and its motion trajectory is periodic; two isolated points projected on the Poincaré maps create attractors, that is to say that the system move in period-two motion.Figure 17
ω = 1.20 and η = 0.650: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)


## 4.2. The Effect of Excitation Amplitude

Figure18 illustrates the bifurcation diagram at ω = 1.20 with respect to the ratio of excitation amplitude ϕ, in which ϕ varied from 0.01 to 0.10. The basement excitation affects the external force of the motion differential equation; therefore, the trajectory of the system is more complex. The mechanical isolation system is weakly affected by the excitation amplitude at minor excitation amplitude. However, its vibration amplitude will increase with the excitation amplitude. The dynamic behavior of the system stayed at state of period-5 throughout the whole period.Figure 18
Dynamic bifurcation diagram: (a)x1; (b) y1.
(a)
(b)Figures19–22 illustrate the dynamic steady response at ϕ = 0.020, 0.055, 0.075, and 0.100, respectively. The numerical simulation results indicate that the region value of system vibration, with the increase of excitation amplitude, grows rapidly and vibration amplitude also rises. The beat vibration and 1f are mainly shown on displacement response and frequency, which is the interaction of both system’s unbalance magnitude and basement excitation; thus, oil film oscillation appears. Furthermore, the other frequency components cannot be easily actuated under these specific values compared with the first mode frequency which indicates that the unbalance force and nonlinear oil force of the system play a dominant role during the whole process of operation.Figure 19
ω = 2.20 and ϕ = 0.02: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 20
ω = 2.20 and ϕ = 0.055: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 21
ω = 2.20 and ϕ = 0.075: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Figure 22
ω = 2.20 and ϕ = 0.100: (a) displacement response; (b) frequency spectrum; (c) rotor orbit; (d) Poincaré maps.
(a)
(b)
(c)
(d)Due to the basement excitation, the external force of equations is affected, as a result the rotor orbit has changed and clearly shows several single-closed curves. Besides, we note that the trajectory will expand outward iny direction as the excitation amplitude increases gradually when we compare the trajectories under various excitation values. More importantly, the points on the Poincaré map could be numbered and five isolated simple attractors appear apparently on its section, which demonstrates incontestably that the system moves into period-5 motion.


## 5. Conclusions

In this paper, an investigation into the nonlinear dynamic behavior of the mechanical isolation system coupled with air-bag and floating-raft under lateral basement excitation is presented. First, the coupling effects between the excitation source and isolation system are considered, the mechanical isolation model subjected to the basic excitation is built, and its motion equation is developed and dynamic responses are mainly investigated by using the displacement response, frequency spectrum, rotor orbit, Poincaré maps, and the bifurcation diagram. The result predicts that period-5 is mainly performed, and with the increase of rotor speed, the system moves into quasi-bifurcation. However, the system stays in chaos state at high rotor speed, and the vibration amplitude rises rapidly until against bearing bush. Furthermore, the effects of basement excitation on the mechanical isolation system are mainly concentrated on the stage of lower rotor speed, but with the increasing speed, the effects become weak and at the same time, the vibration amplitude reduces significantly. The points projected on the Poincaré section are five, three, or two solitary attractors, in which the system stays in periodic motion. Above all, the dynamic characteristics can provide theoretic supporting for the dynamic, vibration control and its parametric optimization of the marine mechanical isolation system coupled with air-bag and floating-raft on the motion base. Besides, it is of importance for a designer to have at their disposal dynamic models to predict the vibrational behavior of rotating machines operating under base excitations at the design stage to avoid expensive design modifications or in-field catastrophic failure.



---

*Source: 8572178-2020-01-09.xml*