Chaotification can be employed to weaken or eliminate the feature of line spectra of waterborne noise. The efficiency of this method lies on the use of small control. The analysis reveals that the critical control gain depends on the stiffness of vibration isolation systems. Thus, an isolation raft system based on quasi-zero-stiffness (QZS) property is proposed for line spectrum chaotification. A nonlinear time-delay controller is derived accordingly. Comparative analysis shows that the new approach allows much smaller control, and the intensity of line spectra is further reduced. Numerical simulations also indicate other advantages with the introduction of QZS system into chaotification.
National Natural Science Foundation of China11602084Natural Science Foundation of Hunan Province2017JJ30962016JJ4036Education Department of Hunan Province17A0881. Introduction
The spectra of the radiated waterborne noises of underwater vehicles are usually divided into two categories. One is continuous broadband spectra corresponding to noise induced by water current and bubble burst. The other is discrete narrowband line spectra mainly caused by propeller and machinery vibrations. Line spectra consist of closely spaced spectral lines that characterize the intensity of noise components. Line spectra of radiated noise from machinery vibration signify the features of operating underwater vehicles, such as the speed and distance of moving objects. Thus, line spectra have been generally utilized for identifying underwater vehicles and thus regarded as harmful signals for the vehicle safety.
Extensive efforts have been made to attenuate machinery vibrations for the reduction of noise spectra through vibration isolation techniques [1–3]. It has been known that the linear vibration isolation techniques can reduce the intensity of line spectra in a certain frequency bandwidth but cannot change the structure of line spectra of vibration noises due to the frequency fidelity. It means that a sinusoidal input of vibration source to an isolation system will simply turn out a sinusoidal output unable to eliminate the line spike at the frequency. To cope with the problem, an innovative idea [4] was proposed to change the spectra configuration by chaotification [4–7] rather than only rely on vibration attenuation. With a special mechanism, the response of the simple harmonic excitation from the vibration source to the base can present a chaotic state. The process can be called line spectrum chaotification, and the line spectrum structure without control and the line spectrum structure with control are shown in Figure 1. Chaotification is a new technique of making use of broadband nature of chaos to blur and change the feature of line spectra. The mechanism is that a nonlinear isolation system, installed in between the vibrating machine and supporting base, is designed with control which enables to convert a periodic excitation (narrowband input of machinery vibration) into a chaotic one (broadband output on a supporting base) to improve the concealment capability of underwater vehicles.
Comparison of line spectrum structure before (a) and after (b) chaotification.
Lou et al. [4] first reported that a Duffing-type vibration isolator excited by harmonic forces might generate a broadband continuous spectrum at special parameter settings. Through a chaotification process, the intensity of line spectra transmitted to the base could be reduced. The feasibility of this approach was verified thoroughly in an experiment [5]. To enable chaotification for wider range of parameter settings, a parametric perturbation scheme [6] was proposed by using generalized chaos synchronization [8]. However, the persistence of chaotification may not be guaranteed since this scheme could deteriorate the system installation stability. This problem could be avoided by using a state perturbation approach [9]. Inspired by the work [8], chaotification [10] could be realized by employing the projective synchronization method [11], where a nonlinear vibration isolation system (VIS) is driven into a chaotic state by coupling a Duffing system. However, this approach requires large control energy and seems impractical for applications. We will elaborate later that control energy is very much concerned because not only there is a limitation of energy on board but also larger control energy may result in poor quality when reconstructing line spectra. A discrete impact method [12] based on Lyapunov exponents was proposed for chaotifying a Duffing type of VIS. This method also suffers the limitation for chaotification persistence and small controls.
It is known that the existence of time-delay feedback can extend a simple dynamic system into high dimensional one. This nature could make chaotification of a time-delay system readily applicable. In this regard, a new strategy by using time-delay feedback control was introduced for line spectrum chaotification [13–18]. Stability analysis [14] of a double-layer vibration isolation floating raft system (VIFRS) with a linear time-delay feedback control was carried out, and a set of critical criteria for stability switches are derived. These criteria provide a theoretical guidance for the setting of the system parameters and control parameters to achieve chaotification. Strictly speaking, these criteria are not necessary conditions but chaotification of VIFRS seems more likely with the criteria. For a better understanding, Li et al. [16] studied the dynamic behavior of the double-layer nonlinear isolation raft and revealed a variety of response solutions under time-delay control. Furthermore, chaotification of a two-dimensional VIFRS with dual time-delay feedback control was investigated [18]. A more sophisticated methodology originally developed from antichaos control [19, 20] could be employed for line spectrum reconstruction [17] of double-layer VIFRS. The significance of this work [17] is to provide a standard procedure for time-delay controller design which promises the occurrence of chaotification in the sense of Li and Yoke criteria [19]. This method based on nonlinear time-delay control notably outperforms the linear time-delay approaches [14, 18], especially in the requirement of small control. Different from the previous ideas, Zhou et al. [15] proposed a chaotification method totally based on spectrum optimization. This method can realize chaotification of VIFRS without exactly knowing system parameter settings or operational conditions, which is effective, easy to use, and handy in real applications.
The existing research works mainly focus on the generation of chaotic response in the floating raft isolation systems. In the purpose of altering the line spectrum configuration and eliminating spectrum line spikes, abundant results from numerical simulations [13, 15, 17] indicate that small control energy can lead to a better quality of chaotification in terms of the line spectrum intensity and broadness of spectrum bandwidth. Pursuing small control for chaotification is one of the key issues in the line spectrum reconfiguration. In this paper, we will propose time-delay control on a quasi-zero-stiffness (QZS) raft system for the improvement of both chaotification and isolation efficiency. The QZS system is a kind of nonlinear VIS with the characteristics of high static and low dynamic stiffness [21–25]. In the following study, we will first elaborate what factors affect the control energy. Then, we will introduce a newly developed QZS system [26] for line spectra control. Based on the new system, we will derive the nonlinear time-delay controller for chaotification. Numerical simulation will show the advance of the approach. By combining time-delay scheme with the QZS system, we can significantly decrease the critical control gain and meanwhile can effectively change the line spectrum configuration.
2. Critical Control Gain
In this section, we will state the key factors that dominate the minimum required control gain for chaotification, which is referred to as critical control gain. Firstly, we look into a single degree-of-freedom (DOF) VIS and then extend the discussion to a double-layer VIS. Consider a nonlinear isolation system with a time-delay control:(1)x˙=fx+gxKtφτd,where x is a state vector, fx and gx∈Rn are vector field, Kt is a control gain, and φτd is a time-delay feedback control function. The form of φτd is not restricted which may be either linear or nonlinear function.
As shown in Figure 2(a), the mass m is supported by a nonlinear spring with quadratic and cubic nonlinearity, a linear damper, and an actuator which is utilized to implement time-delay feedback control. The governing equation of the single DOF mass-spring system with time-delay feedback control can be written as(2)Mx¨+Cx˙+K1x−Δx+K3x−Δx3=F0cosω0t+Ktφτd,where Kt is the feedback gain, φτd is the time-delay feedback control function, C is the damping coefficient, K1 and K3 are the stiffness coefficients, F0 and ω0 are the amplitude and frequency of the harmonic excitation, respectively, and Δx denotes the static deformation of the spring subjected to the system static weight.
A single DOF nonlinear VIS: (a) system model and (b) load-displacement curve.
We define the following nondimensional parameters, and setting X=x−Δx,(3)ωn=K1M,ξ=C2Mωn,η=K3M,f0=F0M,ωt=KtM,the governing equation can be rewritten as(4)X¨+2ξωnX˙+ωn2X+ηX3=f0cosω0t+ωt2φτd.
By defining x1=X and x2=X˙, the equation of motion without external excitation force can be written in the autonomous form:(5)x˙1=x2,x˙2=−2ξωnx2−ωn2x1−ηx13+ωt2φτd.
The system possess the equilibrium point A0,0. Giving a perturbed motion Δx1,Δx2 at the equilibrium point A0,0, it may yield two linearized time-delay differential equations:(6)Δx˙1=Δx2,Δx˙2=−2ξωnΔx2−ωn2Δx1+ωt2Δx1t−τd.
Note that the form of time-delay feedback control function φτd may be arbitrary. If it is linear, the linearized part of φτd is Δx1t−τd. If it is nonlinear, φτd can be expressed in Taylor series of ∑anΔx1nt−τd. Because of the vibration amplitude of a real VIS is often very small, the high order terms can be ignored and the dominant term is still Δx1t−τd. A very useful method for determining system stability is bifurcation analysis. According to the theorem proposed [27], the characteristic equation can be written as(7)P1λ+P2λe−λτ=0,where(8)P1=λ2+2ξωnλ+ωn2,P2=−ωt2.
When the real part of a certain eigenvalue changes from negative to zero or even to positive, the time-delay increased may cause the occurrence of Hopf bifurcation. That is to say, there exists a critical time-delayτdc, at which the system with time-delay feedback control will loose stability and then the following characteristic equation has a purely imaginary root [27]. Suppose that λ=iν is the purely imaginary root of the characteristic equation:(9)R1ν+iQ1ν+R2ν+iQ2νcosντ−isinντ=0,where(10)R1ν=−ν2+ωn2,R2ν=−ωt2,Q1ν=2ξωnν,Q2ν=0.
Separating the real and imaginary parts of the abovementioned equation may give(11)R1ν+R2νcosντ+Q2νsinντ=0,Q1ν−R2νsinντ+Q2νcosντ=0.
By squaring the above equations and then summing the results, we may obtain a quadratic equation after defining a new variable μ=ν2:(12)μ2+4ξ2−2ωn2μ+ωn4−ωt4=0.
The roots of equation (12) can be obtained as(13)μ1,2=2−4ξ2ωn2±16ξ2ξ2−1ωn4+4ωt42=1−2ξ2ωn2±Δ.
Due to the fact that 1−2ξ2>0 is always satisfied for engineering applications and if Δ≥0 equation (12) has one positive root at least, then the critical time-delayτdc can also be acquired. So, the system goes into the state of chaos through Hopf bifurcation only when(14)Ktc≥2ξ1−ξ2K1,where ξ is the damping ratio. The relationship in equation (14) reveals that the critical control gain Ktc is only related to ξ and K1 of the isolation system.
Figure 3 plots the variation of the critical control gain versus the equivalent linear stiffness for different damping ratios. Obviously, Ktc is linearly dependent on system stiffness K1 and nonlinearly dependent on the damping ratio ξ. It delivers an important clue that the system can be much easily disturbed into the chaotic state when the system’s stiffness and damping are small. This finding motivates us to look for small stiffness isolation system, which is beneficial not only for chaotification but also for vibration attenuation. We will show later that small control gain leads to small control energy and results in low intensity of chaotified line spectra.
The critical control gain Ktc of a single DOF VIS with time-delay control for chaotification versus different equivalent linear stiffnesses K1 and different damping ratios ξ.
The nonlinear VIFRS [14] can be regarded as a double DOF mass-spring system, as shown in Figure 4. M1 and M2 denote the isolated equipment and the floating raft, respectively. M1 is supported by a linear damper and a nonlinear spring which possesses quadric and cubic nonlinearity. M2 is supported by a linear damper and a linear spring which are connected with a fixed ground base. There is an actuator installed between M1 and M2, which is utilized to implement time-delay feedback control for chaotification.
The structure diagram of a double DOF VIFRS.
The governing equations with time-delay feedback control are given by(15)M1x¨1=−C1x˙1−x˙2−K1−2U1H+3U2H2x1−x2+U1−3U2Hx1−x22−U2x1−x23+F0cosω0t+Ktφτd,M2x¨2=−C2x˙2−K2x2+C1x˙1−x˙2+K1−2U1H+3U2H2x1−x2−U1−3U2Hx1−x22+U2x1−x23−Ktφτd,where C1 is the damping coefficient of the nonlinear vibration isolator; K1, U1, and U2 are the linear, quadric, and cubic stiffness coefficients of the nonlinear vibration isolator, respectively. The equivalent linear stiffness at the static equilibrium point is K0=K1−2U1H+3U2H2, where H=h1−h2, C2 is the damping coefficient of the damper between the floating raft and fixed ground, K2 is the stiffness coefficient of the linear spring, F0 and ω0 are the amplitude and frequency of the harmonic excitation, respectively, Kt is the feedback control gain, and φτd is the time-delay feedback control function.
For the double-layer VIS with time-delay control [14], the derivation of critical control gain is much complicated, and it can only be expressed in an implicit form. We are interested in the critical boundary for the control gain where chaotification can be implemented effectively. To find the critical condition for chaotification in the time-delay system (15), we firstly impose a linearization on the double-layer nonlinear VIFRS at an equilibrium point. Giving a perturbation to the system, it yields a set of linearized time-delay differential equations. By taking the Laplace transform, we can get the characteristic equation from which the eigenvalues are investigated for the stability switches on the parameter domain of time-delay settings and control gains.
Figure 5 illustrates the skeleton region for the stability switches on the parameter plane of the time-delay and control gain. The white region indicates the feasible area where we can apply chaotification, while the region marked in gray color denotes the infeasible area for chaotification. The critical boundaries consist of a number of solid curves. In Figure 5(a), there is a region bracketed by two dashed lines in vertical direction. As the time-delay increases to infinity τ⟶∞, these dashed lines define the boundary for chaotification. The corresponding value Ktc of the control gain is the required minimum control gain, regarded as the critical control gain. For the case where the system stiffness is set at K0=9600N/m, as shown in Figure 5(a), the critical value of the control gain is Ktc=2112N/m. When the system stiffness is decreased to K0=2400N/m, as shown Figure 5(b), the skeleton structure is similar but the critical control gain becomes much smaller, as Ktc=1224N/m. It reveals that when the system stiffness decreases the critical control gain decreases as well, very similar to the case of the single DOF system.
The skeleton for stability switches on the parameter plane of the time delay and control gain, indicating the feasible area (white) and infeasible area (gray) for chaotification when the system stiffness is (a) K0=9600N/m and (b) K0=2400N/m.
To understand the relationship between the critical control gain and the system’s settings, Figure 6 shows the effects of the variation of the equivalent linear stiffness and damping on the critical control gains. With the increase of the equivalent linear stiffness K0, the critical control gain Ktc increases nonlinearly unlike the linear relationship of a single DOF VIS. The curve trend also shows the significant effect for the variation of the damping coefficients C1 and C2 on the critical control gain. It strongly implies that small control energy required for chaotification greatly depends on small stiffness K0 and small damping of the system. Especially, when K0 is small, Ktc drops quickly with the decrease of K0.
The critical control gain for chaotification versus different equivalent linear stiffnesses for different damping (C1 and C2 are the damping coefficients of the double-layer system, respectively).
We know that the feedback control gain Kt represents the level of control energy to be inputted. The increase of the control gain leads to the increase of vibration amplitude generally, and accordingly the intensity of the line spectra increases as well, surely harmful to the concealment capability of underwater vehicles. In order to explain this fact quantitatively, we define the intensity of line spectra as the root mean square values of the normal difference between the spectrum peak values with control and the spectrum value at excitation frequency without control. We can obtain the relationship between the intensity of the line spectra and the feedback control gain, as shown in Figure 7.
The intensity of the line spectra versus feedback control gain σ=50,τ=20.
Figure 7 plots the intensity of the line spectra under the variation of feedback control gain Kt. Obviously, the intensity of the line spectra is linearly dependent on the Kt approximately. As the absolute value of control gain increases, the intensity of line spectra increases proportionally. This diagram indicates that the smaller control leads to the lower intensity of line spectra.
3. Controller for a QZS Vibration Isolation System
In this section, the QZS system [26] is introduced to a double-layer VIS for chaotification due to its favorable feature of quasi-zero-stiffness. Based on this system, the analytical function of time-delay feedback control is derived by using the methodology [17].
In the double-layer isolation system, as shown in Figure 8, we consider a QZS system placed on the upper layer and a linear isolation system placed on the bottom layer. It is known that the QZS system [26] has the characteristics of high static and low dynamic stiffness. At the equilibrium state, the QZS system theoretically has zero stiffness, which is a perfect model for chaotification in terms of small energy control.
The structure diagram of a double-layer QZS system.
The QZS isolation system uses the negative stiffness of magnetic springs to offset the positive stiffness at the equilibrium state yielding the property of zero stiffness, which offers high static stability but very low dynamic stiffness around the equilibrium state. The theoretical analysis and experimental results [26] show the excellent attenuation performance in terms of the force transmissibility in comparison with conventional vibration isolation technology. Note that this mechanical property provides a combination of two advantages. Number one, it can most effectively isolate vibration due to its inherent low natural frequency, especially outperform in the low-frequency band; number two, it allows us to use small control for chaotification.
The stiffness of the QZS system can be modeled by a unique cubic term of displacement [26], while there is no linear stiffness. The equations for the double-layer QZS system with time-delay control can be formulated as follows:(16)M1x¨1=−C1x˙1−x˙2−Kx1−x23+F0cosω0t+Ktφτd,M2x¨2=−C2x˙2−K2x2+C1x˙1−x˙2+Kx1−x23−Ktφτd,where C1 and C2 are the damping coefficient of the QZS isolator and the bottom isolator, respectively, K and K2 are the stiffness coefficient of the QZS isolator and linear spring, respectively, F0 and ω0 are the amplitude and frequency of the harmonic excitation, respectively, Kt is the feedback control gain, and φτd is the time-delay feedback control function to be derived.
The governing equations (16) can be transformed to a standard form of the first-order governing equations, given by(17)x˙1=y1,y˙1=−C1M1y1−y2−KM1x1−x23+F0M1cosω0t+KtM1φτd,x˙2=y2,y˙2=−C2M2y2−K2M2x2+C1M2y1−y2+KM2x1−x23−KtM2φτd.
We use vector x=x1y1x2y2T to denote the system state, where x1 and x2 are the displacements and y1 and y2 are the velocities of M1 and M2, respectively. The controlled system (17) can be expressed in a general form of the single-input and single-output system given by(18)x˙=fx+gxδxt,y=hx,where δxt=Ktφτd is the input of the feedback control, hx is the output function of the system, and(19)fx=y1−C1y1−y2M1−Kx1−x23M1+F0cosω0tM1y2−C2y2M2−K2x2M2+C1y1−y2M2+Kx1−x23M2,gx=01M10−1M2.
A nonlinear time-delay feedback controller for chaotification can be derived based on the differential geometry theory and the definition of Li and Yoke Chaos [19, 20]. It follows the procedure. Firstly, the abovementioned nonlinear isolation system (18) is transformed into a standard linear system by a set of nonlinear transformation functions. The set of nonlinear transformation functions are defined by(20)z=Φx=φx1,x2,···,xn1φx1,x2,···,xn2⋮φx1,x2,···,xnn=hxLfhx⋮Lfn−1hx,where Φx is a partial diffeomorphism and Lfhx,…,Lfn−1hx are Lie derivatives. If the relative degree in a neighborhood of equilibrium point is exactly equal to the degree of the system, the nonlinear isolation system can be exactly transformed into a standard linear system as(21)z˙=Az+Bv,where z is the state vector of the linear system, v is the control function, A is the state coefficient matrix, and B is the control coefficient matrix. A nonlinear time-delay controller for chaotification in this linearized system can be derived through the technique [19, 20].
Based on the lemma below, the output function of hx can be obtained and the control function of δxt for chaotification can be designed accordingly. The Lemma from nonlinear control theory is listed as follows.
Lemma 1.
(see [28]) A nonlinear control system is feedback linearizable on a neighborhood D of an equilibrium point if and only if
rankgxadfgx⋯adfn−1gx=n,x∈D
spangx,adfgx,···,adfn−2gx is involutive on D
Then, a solution of output y=hx from the following partial differential equations can be determined:(22)∂hx∂xgxadfgx…adfn−2gx=0.
Point A0,0,0,0 is an equilibrium point of system (17). Based on this, we can acquire an analytical time-delay control function with Lemma 1. According to the definition, we can obtain adfgx, adf2gx, and adf3gx, respectively. Using the parameter a=x1−x2, we have(23)adfgx=∂g∂xf−∂f∂xg=−1M1C1M12+C1M1M21M2−C1M1M2−C1M22−C2M22,(24)adf2gx=∂adfgx∂xf−∂f∂xadfgx=−C1M1+M2M12M2C1C2M12−3a2KM1M2M1+M2+C12M1+M22M13M22C2M1+C1M1+M2M1M22−C12M1+M22−C1C2M12M1+M2+M1−C22M1+M2K2M1+3a2KM1+M2M12M23,(25)adf3gx=∂adf2gx∂xf−∂f∂xadf2gx=Λ1Λ2Λ3Λ4,where(26)Λ1=−C1C2M12−3a2KM1M2M1+M2+C12M1+M22M13M22,Λ2=2C12C2M12M1+M2+C13M1+M23+C1M1C22M12−M2K2M12+6a2KM1+M22−3aKM12M2aC2M1+2M2M1+M2y1−y2M14M23,Λ3=C12M1+M22+C1C2M12M1+M2+M1C22M1−M2K2M1+3a2KM1+M2M12M23,Λ4=−C13M1+M23−C12C2M13M12+4M1M2+M12+C1M1−C22M13M1+M2+M26a2KM1+M22+K2M12M1+M2+M12−C23M1+C2M22K2M1+3a2K2M1+M2+6aKM22M1+M2y1−y2/M13M24.
After simplification of matrix transformation and inserting the particular values of the system parameters, we can determine the rank of matrix gadfgadf2gadf3g.
Next, we will determine the rank of matrix gadfgadf2gadfg,adf2g in order to estimate the second condition of Lemma 1:(27)g,adf2g=∂adf2gx∂xg−∂g∂xadf2gx=0000,(28)adfg,adf2g=∂adf2gx∂xadfg−∂adfg∂xadf2gx=06aKM1+M22M13M220−6aKM1+M22M12M23.
For system (16), we assign a set of parameters:(29)M1=85 kg,M2=42.5 kg,C1=C2=260Ns/m,K=2.80681×107N/m3,K2=3.84×104N/m,and we can evaluate the relative degree:(30)rankgadfgadf2gadf3g=4.
Based on equations (23)–(25) and (28) and the parameter values, we can obtain(31)rankgadfgadf2gadfg,adf2g=3.
Thus, we conclude that the linear subspace Δ=spangadfgadf2g is also involutive on a neighborhood D of the equilibrium point A0,0,0,0. The controlled system (18) has a relative degree of 4 at A0,0,0,0.
From condition (22) in Lemma 1, we can derive the output function of y=hx from the partial differential equations:(32)∂hx∂xgx=1M1∂hx∂y1−1M2∂hx∂y2=0,(33)∂hx∂xadfgx=−1M1∂hx∂x1+C1M12+C1M1M2∂hx∂y1+1M2∂hx∂x2−C1M1M2+C1M22+C2M22∂hx∂y2=0,(34)∂hx∂xadf2gx=−C1M12+C1M1M2∂hx∂x1+C1C2M12+C1M1+M22M13M22∂hx∂y1+C1M1M2+C1M22+C2M22∂hx∂x2+−C12M1+M22−C1C2M12M1+M2+M12−C22+K2M2M12M23∂hx∂y2=0.
There are multiple solutions for the set of equations (32)–(34). One of the solutions for the abovementioned equations could be derived as(35)y=hx=−K2M1C2M2x1−K2M2−C22C2M2x2+M1M2y1+y2.
According to Wang et al. [19, 20], if the map associated with δxt is a bounded chaotic map and the time-delay is sufficiently large, then the output yt of the system equation with time-delay (17) could be chaotic. The solution of y=hx is not unique and accordingly there are multiple solutions for δxt. We may choose various bounded functional forms. Here, the sinusoidal form of δxt=K˜tsinσyt−τd is utilized for the controller, and we have(36)δxt=K˜tsinσ−K2M1C2M2x1t−τd−K2M2−C22C2M2x2t−τd+M1M2x˙1t−τd+x˙2t−τd.
By inserting the values of the system parameters, we obtain the control function:(37)δxt=K˜tsinσ−295.4x1t−τd−141.6x2t−τd+2x˙1t−τd+x˙2t−τd,where the control gain K˜t, the feedback frequency σ, and the time-delay τd are the control parameters for the nonlinear time-delay feedback controller (37).
4. Chaotification on the QZS Vibration Isolation System
We will illustrate the benefits by proposing the 2-DOF QZS isolation system for chaotification, in comparison with the case of nonlinear isolation system used in [17]. The discussion will be focused on the three aspects including the effective reduction of critical control gain, the suppression of line spectra, and how the feature of line spectra of the system would be weakened or even be eliminated by using the approach.
4.1. Effective Reduction of Critical Control Gain
We are interested in the actual reduction of the control gain when applying chaotification on the 2-DOF QZS VIS. By setting the control parameters at σ=50 and τd=20 s, we will examine the persistence of chaotification across the variation of control gain K˜t, with particular concerns on the critical control gain for the onset of chaos. We will compare the difference between the nonlinear model in [14, 17] and the QZS model proposed in this paper under the same loading conditions.
Figure 9(a) shows the global bifurcation of the nonlinear model [17], and Figure 9(b) shows the global bifurcation of the QZS model (17). Line dots in the vicinity of the origin represent periodic motions, and the cloudy dots represent the onset of chaos. From Figure 9(a), we can see that chaotification occurs when K˜t≥0.452 N. For the corresponding QZS system (17), as shown in Figure 9(b), chaotification starts when K˜t≥0.008 N and maintains across the whole parametric domain. Note that there is a significant difference in the required control gain between the two systems; the minimum control for the QZS system is only about 1.77% of that, for the nonlinear system. It indicates that the QZS system outperforms the nonlinear system [17] in terms of using small control in chaotification.
The system response versus feedback control gain K˜t; the cloudy dots represent chaotic motion and line dots denote periodic motion. (a) Response of nonlinear system [17] and (b) response of the QZS system (17).
Figure 9 also indicates that the small control gain gives rise to relatively low amplitude of system responses, which means small control for chaotification can lead to low intensity of line spectra as well. Apart from the benefit of using small control, the QZS system inherently enables to attenuate vibration in low-frequency bandwidth. Thus, the utilization of the QZS system allows for the synergistic benefits of both chaotification and vibration isolation and is greatly favorable and attractive to applications.
4.2. Suppression of Line Spectra
We may benefit from the introduction of the QZS system (17) into chaotification. A comparison between the system [17] and the QZS system (17) will be carried out to indicate the performance in the reduction of the intensity of line spectra.
We show a numerical example in Figure 10, where the power spectral density (PSD) of chaotified system responses is plotted for system [17] and the QZS system (17), respectively, under the same loading conditions and the external excitation. In general, the spike line of the power spectra induced by the external excitation is covered by the continuous power spectra of chaotified responses. Thus, there is no significant spike lines protruded from the base of power spectra. The profile of the power spectra also indicates the broadband nature of system responses, implying the effectiveness of chaotification for the both systems under control.
Power spectral densities of chaotified responses for (a) the nonlinear system [17] and (b) the QZS system, under the excitation at the frequency of ω/2π=15.9155.
Figure 10(a) shows that the peak value for system [17] is about −71.68 dB, while Figure 10(b) shows that the peak value for the QZS system is about −81.47 dB. Obviously, in terms of line spectrum suppression, the QZS system outperforms system [17] by about 10 dB. For the interval from 0 to 20 in the frequency domain, we can see that the QZS system is more efficient in the reduction of the intensity of line spectra. This advantage lies on two factors. Number one, the control gain required for system [17] is much larger than that of the QZS system. Small control results in low intensity of line spectra. Number two, the QZS system has better attenuation ability, especially in the low-frequency bandwidth. This nature reduces the intensity of power spectra of the response.
4.3. Reconstruction of Line Spectra
In the course of numerical analysis of line spectrum reconstruction, we can observe that chaotification on the QZS vibration isolation system constantly persists across a large range of settings for the time-delayτd and the feedback frequency σ. The profile of power spectra could be altered by choosing different settings of the control parameter pair τd,σ. We are interested in the results where the reconstructed pattern of power spectra is smooth and broadened without line spikes protruded from the base of the power spectra when applying time-delay control. The whole purpose of chaotification is to eliminate the signature of line spectra induced by vibration excitations. In this regard, optimal design [15] of the control parameter pair τd,σ will be useful for obtaining a favorable configuration of line spectra.
In what follows, we will examine the effect of chaotification on the reconstruction of power spectra. There is a significant difference in the patterns of power spectra before and after the implementation of control. Since the formation of line spectra much depends on the excitation frequency, we consider three typical frequencies at the low frequency at 1.6711 which is smaller than the first natural frequency at 1.7311 of system (17), the intermediate frequency at 4.7746 which is in between the first and second natural frequency, and the high frequency at 22.2817 which is far away from the second natural frequency at 5.7234. We will see how the line spikes of power spectra induced by excitations, as shown in the 1st row of Figure 10, can be masked by chaotification at the three frequencies, as shown in the 2nd row of Figure 11.
Profiles of power spectral densities of the QZS system response transmitted to the base without control (the 1st row) and with control (the 2nd row) at different excitation frequencies.
Case 1.
Chaotification at low-excitation frequency at ω/2π=1.6711.
Figure 11(a) illustrates the patterns of power spectra without control. There are two line spikes protruded from the base of power spectra, where the first peak represents a periodic motion at the excitation frequency and the second peak represents a subharmonic motion near the first resonance induced by the excitation. Surely, the system experiences periodic oscillations and the line feature of the power spectra is obvious. On the contrary, Figure 11(b) shows the reconstructed pattern of the power spectra after applying control, where the parameter pair of optimal control [15] is set at τd,σopt=0.8266,1.8026 and the control gain isK˜tc=50N. It can be observed that the new shape of the power spectra completely masks the second line spike but is unable to cover up the first spike. Based on extensive numerical tests, we found that it is relatively difficult to eliminate the line spikes in the low-frequency band in comparison with high-frequency band. Nevertheless, chaotification can weaken the signature of line spectra in general, and the feature of the second spike is completely eliminated.
Case 2.
Chaotification at intermediate frequency at ω/2π=4.7746.
Similarly, Figure 11(c) shows the power spectra without control when applying an excitation at ω/2π=4.7746. The line spike protruded from the base indicates a periodic motion of the system at excitation frequency, and the intensity of the line spike is −68.98 dB. This line spectrum signifies the features of the machinery vibration of vehicles, for example, the operating speed according to the frequency of the line spike and the distance between the noise source and signal detector by the intensity of the line spectrum. To eliminate the features, we implement chaotification by setting the parameter pair of optimal control at τd,σopt=0.2789,1.1964, and the control gain K˜tc=45N. Figure 11(d) shows the reconstructed pattern of the power spectra of the system with the control. It is clear to be seen that the feature of the line spike is completely eliminated and replaced by the chaotic broadband spectra. Moreover, the intensity of the spectra around the excitation frequency is reduced to −79.03 dB. There is about 10 dB reduction. It means that, for certain frequency, the optimal chaotification not only covers the line spectra but also is possible to decrease its intensity.
Case 3.
Chaotification at high frequency at ω/2π=22.2817.
Figure 11(e) shows the line spectrum without control when excited at ω/2π=22.2817, while Figure 11(f) shows the reformed pattern of the power spectra by applying the control at τd,σopt=0.3970,65.0355 and the control gain K˜tc=8.7N. The feature of the line spectrum at the excitation frequency is covered by the chaotic broadband spectra. However, the intensity of power spectra corresponding to the low-frequency bandwidth increases in this case. Nevertheless, the power spectra after chaotification do not carry the information of the signal feature of the vibration source.
The required control gain for effective chaotification in general decreases as the external excitation frequency increases. This characteristic can be observed from Cases 1–3, where the control gain for chaotification reduces corresponding to the frequencies. Figure 12 plots the trend to describe the relationship between the required control gains and the excitation frequency for chaotification. In general, the higher excitation frequency, the less control energy required for chaotification. Meanwhile, the line spectra at high excitation frequency can be eliminated more easily.
The required control gains for chaotification versus the excitation frequency.
5. Conclusion
The most valuable part of this paper is to propose a QZS system for chaotification that can be used to reduce the feature of line spectra induced by machinery vibration of underwater vehicles. This work combines the characteristic of small stiffness of the QZS system with the need of small control in chaotification, leading to the advantages of the reduction of the intensity of line spectra and the improvement of efficiency in the reconstruction of line spectra.
A standard procedure has been presented for the chaotification on a double-layer QZS system. We have shown the derivation of the nonlinear time-delay controller for such a system. Numerical simulations have illustrated the superior performance of the QZS system in terms of the required critical control gain and the suppression of line spectra in comparison with a nonlinear model previous studied in [17], especially outperformed in low-frequency band. In the reconstruction of line spectrum patterns, chaotification can effectively change a narrowband line spectrum induced by a harmonic excitation into a spectrum pattern with broad bandwidth. We have examined the performance of the new system in low-, intermediate-, and high-frequency bands divided according to the natural frequencies of the system. It showed that the feature of line spectra could be completely eliminated for the intermediate- and high-frequency bandwidth. It is relatively difficult to remove the line spikes in the low-frequency band, but the signal feature can be greatly weakened through a chaotification process.
Data Availability
Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to appreciate Professor Daolin Xu and Chunlai Li for their help and support. This research work was supported by the National Natural Science Foundation of China (11602084), Natural Science Foundation of Hunan Province, China (2017JJ3096, 2016JJ4036), and Scientific Research Foundation of the Education Department of Hunan Province, China (17A088).
HillerM. W.BryantM. D.UmegakiJ.Attenuation and transformation of vibration through active control of magnetostrictive terfenol1989134350751910.1016/0022-460x(89)90571-32-s2.0-0024965188JenkinsM. D.NelsonP. A.PinningtonR. J.ElliottS. J.Active isolation of periodic machinery vibrations1993166111714010.1006/jsvi.1993.12872-s2.0-0027912220SpeltaC.PrevidiF.SavaresiS. M.FraternaleG.GaudianoN.Control of magnetorheological dampers for vibration reduction in a washing machine200919341042110.1016/j.mechatronics.2008.09.0062-s2.0-62249170533LouJ.-J.ZhuS.-J.HeL.YuX.Application of chaos method to line spectra reduction2005286364565210.1016/j.jsv.2004.12.0182-s2.0-21444433574LouJ.-J.ZhuS.-J.HeL.HeQ.-W.Experimental chaos in nonlinear vibration isolation system20094031367137510.1016/j.chaos.2007.09.0532-s2.0-65349166654YuX.ZhuS. J.LiuS. Y.A new method for line spectra reduction similar to generalized synchronization of chaos20073063–583584810.1016/j.jsv.2007.06.0342-s2.0-34547916857LiuS. Y.YuX.ZhuS. J.Study on the chaos anti-control technology in nonlinear vibration isolation system20083104-585586410.1016/j.jsv.2007.08.0062-s2.0-37849054092RulkovN. F.SushchikM. M.TsimringL. S.AbarbanelH. D. I.Generalized synchronization of chaos in directionally coupled chaotic systems199551298099410.1103/physreve.51.9802-s2.0-33744833457LiY.XuD.FuY.ZhouJ.Chaotification and optimization design of a nonlinear vibration isolation system201218142129213910.1177/10775463114290542-s2.0-84869045233WenG. L.LuY. Z.ZhangZ. Y.MaC. S.YinH. F.CuiZ.Line spectra reduction and vibration isolation via modified projective synchronization for acoustic stealth of submarines20093243–595496110.1016/j.jsv.2009.02.0492-s2.0-67349236734XuD. L.Control of projective synchronization in chaotic systems200163202720110.1103/physreve.63.0272012-s2.0-85035246715ZhangZ. H.ZhuS. J.LiuS. Y.A new method of chaotic anti-control for line spectra reduction of vibration isolation system3Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering ConferenceAugust 2010New York, NY, USAThe American Society of Mechanical Engineers313610.1115/detc2009-872512-s2.0-77953708145ZhouJ. X.XuD. L.LiY. L.Chaotifing duffing-type system with large parameter range based on optimal time-delay feedback controlProceedings of 2010 International Workshop on Chaos-Fractal Theories and ApplicationsOctober 2010Kunming, China12112610.1109/iwcfta.2010.382-s2.0-78751480710LiY. L.XuD. L.FuY. M.ZhouJ. X.Stability and chaotification of vibration isolation floating raft systems with time-delayed feedback control201121303311510.1063/1.36157102-s2.0-80053421133ZhouJ.XuD.ZhangJ.LiuC.Spectrum optimization-based chaotification using time-delay feedback control201245681582410.1016/j.chaos.2012.02.0152-s2.0-84860196415LiY.XuD.FuY.ZhouJ.Nonlinear dynamic analysis of 2-DOF nonlinear vibration isolation floating raft systems with feedback control2012459-101092109910.1016/j.chaos.2012.06.0102-s2.0-84864071963ZhangJ.XuD.ZhouJ.LiY.Chaotification of vibration isolation floating raft system via nonlinear time-delay feedback control2012459-101255126510.1016/j.chaos.2012.05.0122-s2.0-84864770841LiY.XuD.FuY.ZhouJ.Chaotification of a nonlinear vibration isolation system by dual time-delayed feedback control2013236135009610.1142/s021812741350096x2-s2.0-84880249124WangX. F.ChenG.YuX.Anticontrol of chaos in continuous-time systems via time-delay feedback200010477177910.1063/1.13223582-s2.0-0034366942WangX. F.ChenG. R.YuX.Generating chaos in continuous-time systems via feedback control2003292Berlin, GermanySpringer17920410.1007/978-3-540-44986-7_9CarrellaA.BrennanM. J.WatersT. P.Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic20073013–567868910.1016/j.jsv.2006.10.0112-s2.0-33846430520KovacicI.BrennanM. J.WatersT. P.A study of a nonlinear vibration isolator with a quasi-zero stiffness characteristic2008315370071110.1016/j.jsv.2007.12.0192-s2.0-44649156272CarrellaA.BrennanM. J.KovacicI.WatersT. P.On the force transmissibility of a vibration isolator with quasi-zero-stiffness20093224-570771710.1016/j.jsv.2008.11.0342-s2.0-63749096013ZhouJ.DouL.WangK.XuD.OuyangH.A nonlinear resonator with inertial amplification for very low-frequency flexural wave attenuations in beams201996164766510.1007/s11071-019-04812-12-s2.0-85061326993WangK.ZhouJ.CaiC.XuD.OuyangH.Mathematical modeling and analysis of a meta-plate for very low-frequency band gap20197358159710.1016/j.apm.2019.04.0332-s2.0-85064719082XuD.YuQ.ZhouJ.BishopS. R.Theoretical and experimental analyses of a nonlinear magnetic vibration isolator with quasi-zero-stiffness characteristic2013332143377338910.1016/j.jsv.2013.01.0342-s2.0-84876283039FordeJ.NelsonP.Applications of Sturm sequences to bifurcation analysis of delay differential equation models2004300227328410.1016/j.jmaa.2004.02.0632-s2.0-8644249819IsidoriA.19992Berlin, GermanySpringer