Published: 31 December 2017

Study on active vibration isolation system using neural network sliding mode control

Lihua Yang1
Kai Su2
Shuyong Liu3
Haiping Wu4
Haifeng Li5
1, 4, 5Power Control Department, Navy Submarine Academy, Qingdao, 266042, P. R. China
2Department of Equipment Economics and Management, Naval University of Engineering, Wuhan 430033, P. R. China
3College of Power Engineering, Naval University of Engineering, Wuhan, 430033, P. R. China
Corresponding Author:
Kai Su
Views 303
Reads 156
Downloads 1720

Abstract

In this paper, an active vibration isolation system based on hybrid algorithm is presented in a wide frequency band. Initially, a nonlinear magnetostrictive actuator model is used to establish the appropriate parameters by experiments, which make the actuator using in vibration isolation system work in a better linear dynamic performance, then the sliding mode algorithm modified by neural network, a hybrid algorithm is proposed as the active control controller, and its stability is also analyzed by Lyapunov theory. Furthermore, a dynamic virtual prototype model of active vibration isolation is established to carry out the co-simulation with Adams and Matlab/Simulink, and the results show that under the difference excitations, the neural network sliding mode controller exhibits a good control performance, and the active vibration isolation can effectively improve the vibration isolation effect, reduce the force transmitted to the base and broaden the vibration isolation bandwidth.

1. Introduction

Since the traditional passive vibration isolation technology, such as rubber isolation and spring isolation, was investigated, it has been applied widely to control the vibration in practical engineering for a long time, but its effect is still a little deficient, especially the vibration isolation band is additionally narrow in a low-frequency excitation, which severely affects the device life and reliability [1, 2]. With the development of intelligent actuators and control theory, the passive vibration isolation can’t meet the requirements of vibration control. It gets much more attention to the active vibration isolation [3, 4], which also can be known as the active control.

Magnetostrictive actuator (MA), an intelligent actuator with characteristics of a faster response, much more accuracy positioning, wider bandwidth and larger displacement, has played an important role in micro-precision positioning and active vibration control [5-8]. Francesco [9] calculated the MA amplitude-frequency curve of transfer function, and studied the active control in a single degree of freedom isolation system, the results demonstrate that the force transmitted to the foundation can be effectively reduced by using MA. Based on the LMS algorithm, Zhang Xuhui [10] also applied a magnetostrictive actuator to the active control, and the simulations show that the MA is very suitable for active control. However, the linear model above applications cannot describe the complex dynamic characteristics of the actuator, which aims at maximizing the output force in a wide frequency range and minimizing the using time at the same times. In other words, the linear model, its working frequency band is too narrower to control the vibration efficiently. Therefore, it is necessary to utilize a nonlinear actuator model which can illustrate the transformation between magnetic field, electric field and mechanical field, and further describe its dynamic characteristics in a wide frequency range.

Currently, active control algorithms mainly consist of feedback or feedforward control [11, 12], H2/LQG control [13], LMS control [14], etc. While there are many uncertainties causes in vibration isolation system, such as dynamic stiffness and external time-varying excitation, all of which can affect the isolation effect. It is well known that the control algorithm has a significant influence on the isolation effect. Therefore, these algorithms must be designed with strong robust and excellent control effect against external interferences. Such a solution introduced in the present paper study is based on the sliding mode control, which is independent with system parameters and external disturbances when the system is in sliding mode. However, a strong robustness coexists with certain chattering in control algorithm [15, 16]. In order to reduce the chattering, a hybrid algorithm is proposed by modifying and improving the sliding mode control using neural network theory. It can enhance the dynamic response and achieve a more stable control effect.

In addition, many scholars are mainly concerned about the simulation and verification based on the accurate dynamics analysis, but the theoretical derivation, simulation and calculation are more intricate. The Adams software can carry out co-simulation with MATLAB/Simulink to complete movement and control analysis, which can be used for dynamic characteristic simulation and calculation to solve mechanical problems. This provides a new method for studying the vibration isolation system [17, 18]. Based on the aforementioned analysis, an Adams virtual prototype model is established for vibration isolation system, and then it is utilized to execute the co-simulation using Adams and Matlab/Simulink, which can conveniently evaluate the active control effect.

2. The magnetostrictive actuator model

The MA structure and schematic diagram are indicated by Fig. 1 and Fig. 2. It is mainly made of the push rod, the preload mechanism, the magnetostrictive rod, the driving and bias coil, where the preload mechanism is comprised of preload spring, push rod and preload screw. The polarization magnetic field of the bias coil can not only make magnetostrictive rod deform in a linear region but also prevent the multiplier phenomenon affecting control precision. Then, when an alternating current is passed through the driving coil, the magnetostrictive rod can generate an alternating magnetic field as well, which causes a dynamic redistribution of magnetic field domains in the magnetostrictive rod. Furthermore, the telescopic length can be obtained microscopically to promote a displacement, which can thrust the output push rod. Thus, electromagnetic energy can be converted into mechanical energy.

In Fig. 1 – (1) the top cover screw, (2) the output push rod, (3) the top cover, (4) the preload spring, (5) the jacket, (6) the driving coil, (7) the bias coil, (8) the magnetostrictive rod, (9) the preload screw, (10) the bottom cover screw, (11) the bottom cover.

Fig. 1Structure diagram of the magnetostrictive actuator

Structure diagram  of the magnetostrictive actuator

Fig. 2Schematic diagram of the magnetostrictive actuator

Schematic diagram  of the magnetostrictive actuator

Based on the Jiles-Atherton theory, which reflects nonlinear hysteresis characteristics of the magnetostrictive actuator, it can accurately match the theoretical simulation and experiment by the structure optimization and parameter identification. Then, the partial differential equations of the magnetization process are established by magnetic domain and piezomagnetic theory [19, 20]. Moreover, its dynamics equations are also set up with the consideration of the mass and damping effect of magnetostrictive rod. The specific equations are as given below:

1
He=H+α~M=H+α+9λsσ02μ0Ms2M,
2
Man=MscothHea-aHe,
3
dMirrdH=Man-Mirrδk-α~Man-Mirr,
4
Mrev=bMan-Mirr,
5
M=Mirr+Mrev,
6
λ=32λsMs2M2,
7
x=1(m1+m2)s2+(c1+c2)s+(k1+k2)ASHλ,
8
f=m2s2+c2s+k2(m1+m2)s2+(c1+c2)s+(k1+k2)ASHλ,

where He is the effective magnetic field strength, H is the external driving magnetic field strength. it’s value is the sum of the driving magnetic field Hd=nI (n is the number of turns in per unit coil length, I is a driving current) and bias magnetic field Hb, α is the molecular field parameter of interacting magnetic torque, λs is the saturation magnetostrictive coefficient, σ0 is the constant preload stress, μ0 is the vacuum permeability, Ms is the saturation magnetization, λ is the axial magnetostriction strain, a is the shape factor of non-hysteresis magnetization, k is the irreversible loss coefficient, α~ is the integrated magnetic domain coefficient, b is the reversible coefficient, the sign parameter δ +1 when dH> 0 and δ –1 when dH< 0, Man is the non-hysteresis magnetization, Mirr is the irreversible magnetization, Mrev is the reversible magnetization, M is the total magnetization.c, ρ, L, A, SH are respectively the damping, density, length, cross-sectional area and axial compliance coefficient of the magnetostrictive rod, m1=ρLA/3, c1=cA/L, k1=A/SHL are respectively the equivalent mass, damping and stiffness coefficient of the magnetostrictive rod. m2, c2, k2 are respectively the equivalent mass, damping and stiffness coefficient of the end load. Based on Eqs. (1-8) and experiments, the displacement responses to the preload and the bias magnetic field are shown in Fig. 3 and Fig. 4, then we can optimize and choose proper parameters of the preload and the bias magnetic field, which can make the actuator work perfectly within an approximately linear range and a larger displacement.

Fig. 3Displacement changes with the preload

Displacement changes  with the preload

Fig. 4Displacement changes with the bias magnetic field

Displacement changes  with the bias magnetic field

3. The active control system model

In the present paper, the active vibration isolation system model is shown in Fig. 5, M2, C2, K2 are respectively the mass, damp and stiffness of the upper isolated device.

Fig. 5The active vibration isolation system model based on the magnetostrictive actuator

The active vibration isolation system model based on the magnetostrictive actuator

Here, all of them are respectively equal to m2, c2 and k2 in Eq. (7) or Eq. (8) as well. M1, C1 and K1 are the mass, damp and stiffness of the middle isolation platform. f is the output force of the magnetostrictive actuator, z2, z˙2, z¨2, z1, z˙1 and z¨1 are respectively the displacement, velocity and acceleration of the upper isolated device (M2) and the middle platform (M1). Consequently, the isolation system dynamics equations can be obtained as Eq. (9):

9
M1z¨1+C1z˙1+K1z1+C2z˙1-z˙2+K2z1-z2=f,M2z¨2+C2z˙2-z˙1+K2z2-z1=p-f.

Take Z=[z˙1,z1,z˙2,z2]T as the state variable, and the force transmitted to the foundation as the state space output F=K1z1+C1z˙1. Then, the state equation of Eq. (9) can be written as:

10
Z˙=AZ+B1U+B2W,Y=CZ+D1U+D2W,

where:

A=-C1+C2M1-K1+K2M1C2M1K2M11000C2M2K2M2-C2M2-K2M20010, B1=1/M10-1/M20,
B2=001/M20, C=C1K100T, D1=D2=0000, U=f, W=p.

By Eq. (10), we can get that the system is completely passive vibration isolation if U= 0.

4. Design and analysis of the controller

4.1. System descriptions

Let the system state equation as Eq. (11):

11
x˙(t)=Ax(t)+Bu, xRn, uRn, tR,

where x is the state variable, and u is the control input, when the system reaches the ideal sliding mode S˙= 0, which can be described as Eq. (12):

12
S˙=sxx˙=0.

The solution of Eq. (12) is known as an equivalent control in the sliding mode region, the equivalent control is designed for the absence of external perturbations in certainty system. Its robustness will be significantly affected by disturbances. Therefore, it is important to study on the controller design with disturbances. Then, the system state space is expressed as follows:

13
x˙=Ax+Bu+ξ, xRn, uRn, tR, ξξmax,

where x, u are the state variable and the control input, ξ is the external disturbances. So, the active vibration isolation system is combined specifically with Eq. (10):

14
Z˙=AZ+B1U+B2W,

where:

x=Z, B=B1, u=U, ξ=B2W.

4.2. The design of sliding mode controller using neural network

The sliding mode control is an effective control way, its advantages are that the system is totally unrelated with external disturbance and parameter perturbations in sliding mode [21, 22]. Based on the state space Eq. (14), the switching function can be designed as Eq. (15):

15
St=GZ,

where St is the switching surface, Z=[z˙1,z1,z˙2,z2]T is the state variable, G is the switching matrix. In this paper, the approaching exponential law of the saturated function is adopted as a control method:

16
S˙t=-εsatSt-kStε>0, k>0,

where ε is the speed of the moving point approaching the switching surface, sats is the saturation function, k is the approaching exponential coefficient. As a switch control, there is chattering phenomenon in traditional sliding mode. Therefore, the control parameter ε and k should be modified to improve the controller stability. Here, a linear function is established between ε and S as Eq. (17), which can meet the adaptive requirements for the sliding mode controller:

17
ε=S(t)β,

where β is the modification coefficient, k also has a significant influence on system chattering. Therefore, a hybrid algorithm, the RBF neural network sliding mode controller is designed to eliminate chattering phenomenon, the RBF neural network is a kind of three-layer feedforward neural network and its training speed is 103-104 times as fast as traditional BP neural network [23,24]. With a good generalized ability, it can approximate any continuous function with a higher accuracy. In this paper, it is used to improve the sliding mode controller. Take the Δ=[S(t),S˙(t)] as the network input, and the exponential approaching coefficient k~ as the network output. The details are as followings:

18
k~=wTh(Δ),

where wT is the weight of RBF neural network, and h(Δ) is the Gaussian function:

19
hiΔ=exp-Δ-ci2σi2, i=1,2,3, j=1,2,

where ci=[ci1ci2ci3ci4cij] is the network center vector, σi is the neurons width. Take the displacement of the middle platform as the adjustment indicator:

20
E=12z12.

Then the adjustment algorithm of network weights can be expressed as:

21
Δw=-λEw=-λz1z1w=-λz1z1UUk~k~w-λz1sgnz1UUk~k~w,

where z1/U depends on the sign, its value can be compensated by the weight values. z1 is proportional with U in active control, so it can be written as sgn(z1/U)= 1, U/k~=-sgn(S), k~/w=h(Δ)sgn(wTh(Δ)) and then the weight adjustment algorithm can be modified as:

22
Δwλz1sgnShΔsgnwThΔ.

Moreover, the learning algorithm of network weights is expressed as:

23
wt=wt-1+Δwt+αwt-wt-1,

where λ is the neural networks learning rate, λ(0,1), α is inertial coefficient, α(0,1). Substituting Eq. (14) and Eqs. (16-18) into Eq. (15), we get the RBF neural network sliding mode control law as Eq. (24):

24
U=-GB1-1GAZ+GB2W+StsatStβ+k~St.

4.3. Controller stability analysis

In this paper, the Lyapunov theory is used to prove the stability of the RBF neural network sliding mode controller [25], take the Lyapunov function as Eq. (25):

25
Vz=12S2(t)V˙z<0.

In continuous conditions, the first derivative of the Lyapunov function can be written as:

26
V˙(z)=StS˙t.

By substituting Eqs. (15-17) and Eq. (18) into Eq. (26), we can get:

27
V˙z=StS˙t=StGAZ+GB1U+GB1B2W=S(t)-StsatStβ-k~St=-StStsatStβ+k~S2t0,

where S(t)= 0 is the only equilibrium point when V˙z 0, so the RBF neural network sliding mode controller can satisfy the conditions of Lyapunov asymptotic and overall stability, which shows the controller designed above is stable.

5. Co-simulation

First of all, we need to build a virtual prototype model of the active vibration isolation system in Adams Software (the parameters will be shown in numerical simulation). And the active control force of the magnetostrictive actuator will be simplified to an interaction force between the upper isolated device and the middle platform, external excitation and system output are defined as well. In addition, the RBF neural network sliding mode control algorithm will be compiled in MATLAB/Simulink environment. It can connect input and output interfaces defined by Adams and MATLAB/Simulink. Thus, the closed-loop co-simulation model can be acquired. The Adams virtual prototype model, a double-layer vibration isolation platform consisting of a water pump, is illustrated in Fig. 6, and Adams and MATLAB/Simulink co-simulation model is shown in Fig. 7.

Fig. 6The Adams virtual prototype model of the vibration isolation system

The Adams virtual prototype model of the vibration isolation system

Fig. 7The co-simulation model of Adams and MATLAB/Simulink

The co-simulation model of Adams and MATLAB/Simulink

6. Numerical simulations

The system modal frequency is f1= 3.2 Hz, f2= 16.7 Hz by Adams analysis, in order to verify the model correctness and control reasonableness. p1=sin2π*4*t,p2=sin2π*30*t,p3=sin2π*5*t+0.2sin2π*30*t+0.1sin2π*60*tand random excitation are applied to calculate the force transmitted to the foundation, which can obtain the isolation system effect, and some simulation parameters are: n= 1200, Hb= 20 kA/m, a= 7102 A/m, α~= –0.001, Ms= 7.65×105 A/m, λs= 1.005×106,k= 3283 A/m, b= 0.18, SH= 1.3×1011, d= 1.0×108 m2/N, c= 3000 kNs/m2, ρ= 9250 kg/m3, L= 8.6×104m, A= 78.5×106 m2, M1= 22.5 kg, C1= 290.544 Ns/m, K1= 190000 N/m, M2= 105.42 kg, C2= 302.7 Ns/m, K2= 57016.2 N/m, the specific simulation results are shown as Figs. 8-12.

The control-off and control-on respectively denote the passive vibration isolation and the active vibration isolation. Fig. 8-11 and Table 1 show that, compared with the passive vibration isolation under single-frequency, multi-frequency and random signal excitation, the RMS values of the force transmitted to the base are greatly reduced by active control, and the external excitation is also effectively attenuated. While the former is respectively 51.3769 %, 4.7628 %, 93.0617 % and 32.9105 % in passive vibration isolation, and 4.3784 %, 0.5231 %, 2.0644 % and 7.6101 % in active vibration isolation. Therefore, the active vibration isolation system, based on the neural network sliding mode control, has an obvious control effect.

Fig. 8The force transmitted to the foundation with p1=sin(2π*4*t) excitation

The force transmitted to the foundation  with p1=sin(2π*4*t) excitation

Fig. 9The force transmitted to the foundation with p2=sin(2π*30*t) excitation

The force transmitted to the foundation  with p2=sin(2π*30*t) excitation

Fig. 10The force transmitted to the foundation with p3=sin(2π*5*t)+0.2sin(2π*30*t)+0.1sin(2π*60*t) excitation

The force transmitted to the foundation with p3=sin(2π*5*t)+0.2sin(2π*30*t)+0.1sin(2π*60*t) excitation

Fig. 11The force transmitted to the foundation with white noise excitation

The force transmitted to the foundation with white noise excitation

Table 1Vibration isolation performance comparison

4 Hz
15 Hz
5 + 30 + 60 Hz
Random excitation
RMS of excitation
2.8869
2.8869
2.8869
3.1537
RMS of passive isolation (control-off)
1.4832
0.1375
2.6866
1.0379
RMS of active isolation (control-on)
0.1264
0.0151
0.0596
0.2407

The force transmissibility curve is shown in Fig. 12, which demonstrates the natural frequencies corresponding to the two peaks of the passive isolation system, are perfectly matched with theoretical calculation, and the active isolation in first and second-order modal obtain a perfect isolation effect with 20 dB and 15 dB. The cross point frequency is 1.23 Hz, and there is a poor effect on the active vibration isolation when the excitation frequency is lower than the cross point, but the active vibration isolation has a better isolation effect when the excitation frequency is higher than the one, which will not only make up for the lack of passive vibration isolation in low-frequency, but also effectively broaden the system isolation frequency bandwidth.

Fig. 12The force transmission with control-on and control-off

The force transmission with control-on and control-off

7. Conclusions

Based on Jiles-Atherton theory and experiments, a nonlinear model of the magnetostrictive actuator is used for the active vibration isolation in a wide frequency range. And the sliding mode controller is modified by neural network to eliminate chattering problems, its stability is analyzed by Lyapunov theory as well. At last, the Adams dynamics model of the active vibration isolation is also established to carry out the co-simulation using Adams and MATLAB/Simulink. The results exhibit that, compared with the traditional passive vibration isolation, the active control with RBF neural network sliding mode control is more effective under the conditions of single-frequency, multi-frequency and random excitation. And the neural network sliding mode controller has a good control performance, widening the isolation band, improving system response speed, effectively reducing the force transmitted to the foundation and eliminating the first and second-order system resonance peak, all of which can improve the isolation performance. Moreover, the algorithm proposed in this paper has some practical meaning for the design and engineering applications on the isolation system, and the co-simulation modelling method also provides a new research method for the isolation system.

References

  • Zuo L., Jacques J., Slotine E. Robust vibration isolation via frequency-shaped sliding control and modal decomposition. Journal of Sound and Vibration, Vol. 285, Issue 2005, 112, p. 3-1149.
  • Hoque M. E., Mizuno T., Ishino Y., Takeshi M. A six-axis hybrid vibration isolation system using active zero-power control supported by the passive weight support mechanism. Journal of Sound and Vibration, Vol. 329, Issue 17, 2010, p. 3417-3430.
  • Kyihwan P., Sangyoo K., Dongyoub C. An active vibration isolation system using a loop shaping control technique. IEEE/ASME International Conference on Mechatronics and Embedded Systems and Applications, Beijing, China, 2008, p. 586-590.
  • Kim S. J., Dean R., Flowers G., Dean R. Active vibration control and isolation for micromachined devices. Journal of Mechanical Design, Vol. 131, Issue 9, 2009, p. 1-6.
  • Vimalajuliet Sathishkumar R. A., Prasath J. S. Terfenol-D: a high power giant magnetostrictive material for submarine mapping. International Journal of Engineering Science and Technology, Vol. 2, Issue 12, 2010, p. 7165-7170.
  • Sato E., Fujita T. Semi-active seismic isolation system with controllable friction dampers using hydraulic system driven by giant magnetostrictive actuators. Transactions of the Japan Society of Mechanical Engineers, Part C, Vol. 73, Issue 6, 2007, p. 1723-1730.
  • Hu S. F., Zhu S. J., Lou J. J. High-precision control of giant magnetostrictive actuator based on cmac neural network. Journal of Vibration and Shock, Vol. 28, Issue 3, 2009, p. 68-72.
  • Lau H. Y., Liu K. P. Feasibility of using GMM based actuators in active control of journal bearing system. Proceedings of the World Congress on Engineering, London, U.K, 2009, p. 1493-1498.
  • Braghin F., Cinquemani S., Resta F. A model of magnetostrictive actuators for active vibration control. Sensors and Actuators A: Physical, Vol. 165, Issue 2, 2011, p. 342-350.
  • Zhang X. H., Liu Y. G., Fu Y. L. Novel model of active vibration isolation based on giant magnetostrictive actuator. Journal of Beijing University of Aeronautics and Astronautics, Vol. 33, Issue 11, 2007, p. 1317-1320.
  • Li Y. L., Xu D. L. Dynamic effects of delayed feedback control on nonlinear vibration isolation floating raft systems. Journal of Sound and Vibration, Vol. 333, Issue 13, 2014, p. 2665-2676.
  • Chang C. M., Wang Z., Spencer B. F. Application of active base isolation control. Sensors and Smart Structures Technologies for Civil, Mechanical, and Aerospace Systems, Vol. 7292, Issue 39, 2009, p. 1-12.
  • Wang Y., Dyke S. Modal based LQG for smart base isolation system design in seismic response control. Structural Control and Health Monitoring, Vol. 20, Issue 5, 2013, p. 753-768.
  • Karthikeyan C., Prakash S. FxLMS algorithm with feedback neutralization for active vibration control. Institute of Smart Structures and Systems, Vol. 1, Issue 1, 2012, p. 23-33.
  • Lee T. Y., Chen P. C. Sliding mode control for nonlinear isolated bridges. Journal of Earthquake Engineering, Vol. 15, Issue 1, 2011, p. 582-600.
  • Yan X. G., Spurgeon S. K., Edwards C. Robust sliding mode observer with parameters estimation for a class of nonlinear time-delay systems. Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference, Atlanta, Georgia, USA, 2005, p. 987-992.
  • Zhu D. L., Qin J. Y., Zhang Y. Research on co-simulation using ADAMS and MATLAB for active vibration isolation system. International Conference on Intelligent Computation Technology and Automation, Changsha, China, Vol. 2, 2010, p. 1126-1129.
  • Zhou J. P., Zhang Z. Y., Feng G. P., Hua H. X. Modeling and simulation of active-passive vibration isolation of machinery with flexible base. Journal of Vibration and Shock, Vol. 27, Issue 11, 2008, p. 97-100.
  • Pop N. C., Caltun O. F. Jiles-atherton magnetic hysteresis parameters identification. Acta Physica Polonica A, Vol. 120, Issue 3, 2011, p. 491-496.
  • Knypiński Ł., Nowak L., Sujka P., Radziuk K. Application of a PSO algorithm for identification of the parameters of Jiles-Atherton hysteresis model. Archives of Electrical Engineering, Vol. 61, Issue 2, 2012, p. 139-148.
  • Liu J. K. MATLAB Simulation for Sliding Mode Variable Structure Control. Tsinghua University Press, Beijing, 2005.
  • David Y. K., Vadim I. U. A control engineer’s guide to sliding mode control. IEEE Transactions on Control System Technology, Vol. 3, 2004, p. 156-164.
  • Yu S., Alici G., Shirinzadeh B., Smith, J. Sliding mode control of a piezoelectric actuator with neural network compensating rate-dependent hysteresis. Proceedings of the 2005 IEEE International Conference on Robotics and Automation Barcelona, Spain, 2004, p. 3641-3645.
  • Ameziane M., Sefriti B., Boumhidi J., Slaoui K. Neural network sliding mode control for a photovoltaic pumping system. Journal of Electrical Systems, Vol. 9, Issue 3, 2013, p. 380-391.
  • Oliveira R. C., de Oliveira M. C., Peres P. L. Convergent LMI relaxations for robust analysis of uncertain linear systems using lifted polynomial parameter-dependent Lyapunov functions. Systems Control Letters, Vol. 57, Issue 8, 2008, p. 680-689.

About this article

Received
09 May 2017
Accepted
21 July 2017
Published
31 December 2017
SUBJECTS
Vibration generation and control
Keywords
passive vibration isolation
magnetostrictive actuator
sliding mode algorithm
neural network
co-simulation
Acknowledgements

This work was supported by the National Natural Science Funds of China (51579242, 51509253) and the Key Scientific Research Foundation Project of Naval Submarine Academy (42551B1408). And the authors also wish to thank the anonymous reviewers for their valuable comments and suggestions.