# Resonance Suppression for Hydraulic Servo Shaking Table Based on Adaptive Notch Filter

Authors: Jianjun Yao; Zhenshuai Wan; Yue Zhao; Jie Yu; Chen Qian; Yu Fu
Journal: Shock and Vibration
 (2019)

---

## Abstract

In dynamic structure test, the specimen of hydraulic servo shaking table contains not only inertia load but also elastic load. The specimen herein is simplified as a spring-mass-damping system, and the mathematical model of the hydraulic servo shaking table is established by theoretical analysis. The coupling between specimen’s elastic load and shaking table itself produces resonance phenomenon in the required bandwidth when the elastic load is not negligible, which deteriorates the system’s dynamic performance and even leads to the instability of the control system. Also, the time-varying resonance frequency further aggravates the control performance of the system in the shaking test. In this paper, an adaptive notch filter based on least mean square (LMS) error principle is employed to identify the resonance frequencies online and real-time adjust the parameters of the notch filter. Simulation and experiment results show the effectiveness of frequency identification and resonance mode suppression. Compared with the existing resonance suppression scheme, the proposed method can suppress the appeared resonance mode adaptively.


---

## Body

## 1. Introduction

Considering the repeatability, safety, and cost, the mechanical properties of structure and components are usually tested using dynamic laboratory test rather than field test [1–4]. As one of the most important dynamic structure test devices, shaking table is widely used in industrial fields including seismic simulation, aerospace, and infrastructure [5–8]. Compared with mechanical shaker and electrodynamic shaker, hydraulic servo shaking table has such merits as high power to weight ratio, load stiffness, response speed, and control precision, and these features make it invaluable in the dynamic structure test [9, 10]. In previous studies, the specimen of hydraulic servo shaking table is often considered as inertial load, in other words, without consideration of the effect of elastic load. However, the elastic load of the specimen in the dynamic structure test is too big to be ignored [11]. The coupling between the specimen’s elastic load and shaking table itself produces resonance phenomenon in the required bandwidth when the elastic load is not negligible, which reduces the system’s dynamic performance and even leads to the instability of the control system. Increased accuracy in dynamic structure test requires a high bandwidth. However, the bandwidth cannot be increased arbitrarily due to the appeared resonance modes. Moreover, the resonance frequency can vary during the operation of hydraulic servo shaking table because of various factors such as fabrication errors, hydraulic cylinder leakage, and hydraulic fluid compressibility [12]. This leads to that the conventional control method is difficult to achieve exact tracking response and may result in damage to the specimen.Several methods have been studied to solve the resonance problem [13]. Schmidt and Rehm adopted fast Fourier transform (FFT) to measure resonance frequency and calculate antiresonant frequency of a dual inertia spring system [14]. Wang et al. used an adaptive notch filter with a FFT analyzer based on a sorting algorithm and nonlinear model predictive controller to suppress the low-order torsional vibration and compensate for the dynamic control performance of a helicopter system, respectively [15]. Paolo Mercorelli and Nils Werner utilized adaptive control strategy to deal with the servo valve resonance problem. Experimental results showed that adaptive resonance regulator can adaptively change in accordance with any change in the velocity of the revolution of the engine [16]. Wang et al. adopted the modified resonance frequency detection and reduction method based on an adaptive notch filter to extend the middle frequency range of an industrial servo system [17]. Wang et al. presented an improved adaptive filtered x-LMS algorithm to suppress the resonance in the elastic drive system, which proves the better resonance suppression effect as well as convergence speed than the conventional method [18]. Rahman put forward a discrete time adaptive compensator based on an autotuning algorithm to suppress the time-varying resonance characteristics of a hard drive servo system [19]. Wang et al. proposed a model predictive control method to suppress mechanical resonance of a two-mass servo system; in addition, amplitude limit of the shaft torque dynamic performance was improved [20]. Yan et al. proposed an add-on multirate adaptive control scheme, which was based on a polynomial transformation technique and recursive least-squares algorithm, to compensate the uncertain resonance modes beyond the Nyquist frequency in high-performance mechatronic systems, and the vibration attenuation for uncertain resonances was effectively improved [21]. Besides, fuzzy control, robust control, and neural network were also presented to suppress the resonance modes of system whose resonance frequency varies frequently [22].From the above literature, it is easy to find that the resonance suppression methods are mainly focused on mechanical system. As an integrated industrial equipment, the hydraulic servo shaking table needs to deal with electronic signal, hydraulic signal, and mechanical signal simultaneously, so accuracy and real time are highly required. Therefore, an adaptive notch filter based on LMS error principle is employed to identify the resonance frequencies online and real-time adjust the parameters of the notch filter. The output of resonance frequency detector which estimated all the resonance frequency components is injected into the adaptive notch filter to eliminate the appeared resonance modes. Simulation and experiment results validate the effectiveness of frequency identification and resonance suppression.The rest of the paper is organized as follows. Section2 first presents the hydraulic servo shaking table; afterwards, its control strategy and mathematical model are discussed. Section 3 introduces the adaptive notch filter design procedure and its theoretical analysis. The resonance suppression scheme is discussed in Section 4. Then, in section 5, simulation and experimental results are presented to validate the proposed method. Finally, the main points are concluded in section 6.


## 2. Hydraulic Servo Shaking Table System

Figure1 shows the hydraulic servo shaking table system, which is a real-time operation system based on online rapid prototyping control technology xPC-Target. The host computer serves as the user interface and downloads the compiled program to the target computer by Ethernet adapter. The target computer converts voltage analog signal into current drive signal and performs real-time execution of compiled program. The hydraulic pump is used to provide hydraulic energy for shaking table and controls the pressure and flow of servo valve. The signal conditioner regulates the command signal and feedback signal to accurately reproduce the reference signal [23].Figure 1
The hydraulic servo shaking table system.Figure2 shows the control strategy of a hydraulic servo shaking table that mainly includes an input filter, servo valve, hydraulic cylinder, shaking table, and various transducers [24]. As the most common controller for the shaking table, the three variable controller (TVC) includes TVC feedforward controller and TVC feedback controller. The TVC feedback controller comprises of displacement, velocity, and acceleration feedback signals, where the displacement feedback signal is used to improve the closed-loop system stability, the velocity feedback signal is utilized to increase natural frequency of system, and the acceleration feedback signal is adopted to improve system’s damping ratio. The displacement and acceleration feedback signal are measured by displacement transducer and acceleration transducer, respectively. The velocity feedback signal is obtained from displacement feedback signal in low frequency range and acceleration feedback signal in high frequency range by signal conditioner because it is difficult to measure the three feedback signals simultaneously. The TVC feedforward controller including displacement, velocity, and acceleration feedforward signals generated by the input filter is mainly used to extend frequency bandwidth and reduce the tracking error. The error signal, which is the difference between reference signal and feedback signal, drives the hydraulic cylinder generating desired movement by servo valve [25].Figure 2
Control strategy of the hydraulic servo shaking table.The TVC control strategy is based on closed-loop position design, but the shaking table requires acceleration control. So, the reference displacement, velocity, and acceleration signals are obtained by the input filter [26]. The transfer function from input acceleration signal to displacement output signal can be deduced as(1)Gre=drr=Kus2+v0s+d0=Kud0s2/ω02+2ζ0/ω0s+1,where d0=ω02, v0=2ζ0ω0, r is the reference acceleration, Ku is the acceleration gain, ω0=0.4 is the stating frequency of acceleration control, and ζ0=0.6 is the system’s damping ratio.The transfer function of the TVC feedforward controller can be expressed as(2)Gr=Kdr+Kvrs+Kars2,where Kar, Kvr, and Kdr are corresponding acceleration, velocity, and displacement gain of feedforward controller. ar, vr, and dr represent the reference acceleration, velocity, and displacement, respectively.The transfer function of the TVC feedback controller can be described as(3)Gf=Kdf+Kvfs+Kafs2,where Kaf, Kvf, and Kdf are feedback part’s acceleration, velocity, and displacement. af, vf, and df represent the measured acceleration, velocity, and displacement, respectively.The servo valve regulates oil flow into the actuator chambers by changing the position of valve spool. The response of hydraulic oil flowQ to current control signal I can be approximated by the following second order transfer function(4)Gsvs=QI=Ksvs2/ωsv2+2ζsv/ωsv+1,where Ksv, ωsv, and ζsv are gain, natural angular frequency, and damping coefficient of the servo valve.The schematic diagram of hydraulic actuator controlled by servo valve is shown in Figure3, where the specimen is simplified into a spring-mass-damping system. For convenience, ps and po are supply and return pressure, q1 and q2 are input and output flow, p1 and p2 are input and output pressure, V1 and V2 are input and output oil volume, Ci and Ce are internal leakage and external leakage coefficient, m is the mass of the actuator, ms is the mass of the specimen, xv is the valve position, B is the viscous damping coefficient of the specimen, K is the spring stiffness of the specimen, A is the effective area of the piston, xp is the displacement of the actuator, xs is the displacement of the specimen.Figure 3
Schematic diagram of hydraulic actuator controlled by servo valve.The flow rate equation of the servo valve can be expressed as(5)qL=Kqxv−KcpL,where qL=q1+q2/2 is load flow, pL=p1−p2 is load pressure, and Kq and Kc are flow gain and flow-pressure coefficient, which can be defined as(6)Kq=∂qL∂xv=CdW1ρps−pL,Kc=−∂qL∂pL=CdWxv2ps−pL1ρps−pL,where Cd is flow coefficient, W is the gradient of valve orifices, and ρ is the density of hydraulic oil.In consideration of oil elasticity, leakage, and chamber volume variation, the flow equation for each chamber can be formulated as(7)q1=Adxpdt+V1βe⋅dp1dt+Cip1−p2+Cep1,(8)q2=Adxpdt+Cip1−p2−V2βe⋅dp2dt−Cep2,where βe is the effective bulk modulus of hydraulic oil, p1=ps+pL/2, p2=ps−pL/2.Combining equations (7) and (8) gives(9)qL=Adxpdt+CtpL+Vt4βe⋅dpLdt,where Ct=Ci+Ce/2 is the total leakage coefficient and Vt=V1+V2 is the total chamber volume.Ignoring the oil mass and coulomb friction, the force balance equation of the shaking table and specimen can be expressed as(10)ApL=md2xpdt2+Bdxpdt−dxsdt+Kxp−xs,(11)Kxp−xs+Bdxpdt−dxsdt=msd2xsdt2.Combining equations (5) and (9)–(11) and performing Laplace transform, the transfer function from the spool displacement to the displacement of hydraulic actuator can be derived as(12)Gs=Kq/Amss2+Bs+Kmms/Khs5+Bm+ms/Kh+Kcemms/A2s4+Km+ms/Kh+KceBm+ms/A2+mss3+KKcem+ms/A2+Bs2+s,where Kce=Kc+Ct is the total pressure flow coefficient considering leakage. Kh=4βeA2/Vt is the hydraulic spring stiffness.The specimen and shaking table itself form into resonance system due to the specimen’s elastic load. The coupling term occurring in transfer function has a great influence on the specimen in the dynamic structure test of the hydraulic servo shaking table. Therefore, the adaptive notch filter must be designed to accommodate all the frequency variation.


## 3. Adaptive Notch Filter

The notch filter attenuates frequency characteristic to zero at center frequency while keeping the original value unchanged at other frequencies. However, the resonance frequency can vary due to various factors such as friction, dead zone, and manufacturing tolerances. The variation of resonance frequency may deteriorate the suppression performance of the notch filter. The notch filter with wider notch can attenuate to some extent the varying resonance mode. As the notch filter becomes wider, it also induces greater magnitude and phase lag at lower frequencies resulting in a lower bandwidth system. In order to solve the above problem, the adaptive notch filter whose center frequency varies online to track the resonance frequency is proposed to online real-time suppress the appeared resonance model.The LMS algorithm is widely used in the adaptive filter due to its computational simplicity, unbiased convergence, and stable behavior. It is an iterative gradient descent algorithm that changes each iteration based on imperfect gradient estimate to seek the optimum value on the performance surface. The LMS algorithm can be written as follows [27, 28]:(13)yk=wkTxk,εk=dk−yk,wk+1=wk+αεkxk,where yk is output signal, xk is the reference input signal, wk is the updating filter weight vector, dk is the desired signal and εk is the error signal, and α denotes convergence factor that controls the stability and speed of adaptation. The signal propagation process from the reference input to the system output based on the LMS algorithm is shown in Figure 4 in detail.Figure 4
Flow diagram of the LMS adaptive algorithm.The sampled reference inputs are(14)x1k=CsinωrkT+θ=CejωrkTejθ−e−jωrkTe−jθ2j,x2k=CcosωrkT+θ=CejωrkTejθ+e−jωrkTe−jθ2.Following the path through the weight, we obtain(15)Zεkx1k=CejθEze−jωrT−e−jθEzejωrT2j,Zεkx2k=CejθEze−jωrT+e−jθEzejωrT2,where Eze−jωrT, EzejωrT are Ez rotated counterclockwise and clockwise around the unit circle through an angle.The weights are obtained as follows:(16)W1z=αCUzejθEze−jωrT−e−jθEzejωrT2j,W2z=αCUzejθEze−jωrT+e−jθEzejωrT2.The contribution to the output atYz is(17)Y1z=Cejθw1ze−jωrT−e−jθw1zejωrT2j,Y2z=Cejθw2ze−jωrT+e−jθw1zejωrT2.Combining equation (16) into equation (17), we obtain(18)Yz=αCEzUze−jωrT+UzejωrT2.The term in equation (18) represents the time invariant part of the response from Ez to Yz. The time invariant transfer function is(19)Gz=YzEz=αC2zcosωrT−1z2−2zcosωrT+1,where Gz has poles on the unit circle at z=e±jωrT, and a zero at z=1/cosωrT.The transform function fromdk to εk can be expressed as(20)Gnotchz=EzDz=z2−2zcosωrT+1z2−21−αC2/2zcosωrT+1−αC2.For ease of exposition, equation (20) can be rewritten as follows:(21)Gnotchz=α2z2+α1z+α0z2+β1z+β0.The coefficients of the above adaptive notch filter are defined as [29](22)α2=1+2ξaΩc+Ωc21+2ξΩc+Ωc2,α1=2Ωc2−11+2ξΩc+Ωc2,α0=1−2ξaΩc+Ωc21+2ξΩc+Ωc2,β1=2Ωc2−11+2ξΩc+Ωc2,β0=1−2ξΩc+Ωc21+2ξΩc+Ωc2,where ξ and a determine the width and depth of the notch filter, respectively and Ωc=tanωcT/2 is related to the center frequency ωc.Figure5 illustrates the amplitude-phase characteristic of the notch filter for different parameters. It is easy to see that the adaptive notch filter has larger notch width but larger phase delay as the same time. The larger phase delay may affect the system stability. Thus, the notch width should be kept as narrow as possible for the servo control system. In order to effectively suppress the variation of resonance frequency, it is necessary to online adjust the parameters of the notch filter with the varying resonance modes.Figure 5
Bode diagram of the adaptive notch filter with different values. (a)a with constant and ξ with different values. (b)ξ with constant and a with different values.
(a)
(b)


## 4. Resonance Suppression Scheme

Inspired by adaptive noise reduction, which adds the sinusoidal harmonics with just the right amplitude and phase to the primary signal including a noise signal to cancel the noise component, a new resonance suppression scheme based on the adaptive notch filter is shown in Figure6. It is generally known that the peak of the frequency response in magnitude reaches the maximum value at resonance frequency. The signal with various resonance frequency is used as excitation signal to extract the resonance frequency component. Then, the resonance frequency is estimated by resonance frequency detector. With being placed in cascade between control system and hydraulic servo shaking table, the proposed adaptive notch filter can suppress the resonance modes with the variation of resonance frequency in real time.Figure 6
Schematic diagram of the resonance suppression.The resonance frequency detector consists of an excitation signal generator and FIR filter. The excitation signal can be represented as(23)xkexc=∑i=1NaicosωiTk+θi,where ωi denotes a group of resonance frequency candidates, which are spread around the actual resonance frequency. ai and θi are the amplitude and phase of frequency components, respectively.The transfer function of the FIR filter is given by(24)Hz=1−2λz−1+z−2.The frequency response of the FIR filter can be expressed by(25)HejωT=1−2λe−jωT+e−2jωT=e−jωTejωT−2λ+e−jωT=2cosωT−λ.It is easy to see that the frequency is zero whenλ=cosωT for λ≤1. The amplitude of FIR filter output reaches the minimum value when the null frequency coincides with the frequency of the component having the largest amplitude among all the frequency components.According to the linear system theory, the estimation output of the FIR filter is derived as(26)yk=∑i=1NaiHejωTcosωiTk+θi+∠HejωT.The average power of the filter output signalyk is defined as the objective function:(27)J=1L∑k=1Lyk2=12∑i=1Nai2HejωiT2.Substituting equation (25) into equation (27), the objective can be rewritten as(28)J=2∑i=1Nai2cosωiT−λ2.It can be found that the objective function is a quadratic function of the filter coefficientλ. Taking the gradient of the objective function with respect to the filter coefficient and setting it equal to zero:(29)∂J∂λ=4∑i=1Nai2λ−∑i=1Nai2cosωiT=0.Thus, the optimal filter coefficientλ∗ is found when the objective function reaches the minimization:(30)λ=λ∗=∑i=1Nai2cosωiT∑i=1Nai2.Suppose thatωr is the resonance frequency. Because the resonance frequency component has much larger amplitude than the other candidate components, the optimal coefficient λ∗ in equation can be approximated as follows:(31)λ∗=a12cosω1T+a22cosω2T+⋯+aN2cosωNTa12+a22+⋯+aN2=a1/ar2cosω1T+a2/ar2cosω2T+⋯+a2/ar2cosωNTa1/ar2+a2/ar2+⋯+a2/ar2≈cosωreT.The recursive algorithm for updating the filter coefficient can be described as(32)λk+1=λk+12μykxk−1,where μ is the adaptive gain.(33)Ωc2=tanωcT22=tanωreT22=1−cosωreT1+cosωreT=1−λ∗1+λ∗.The optimal coefficients of the adaptive notch filter can be determined by(34)α2∗=1+ξa1−λ∗21+ξ1−λ∗2,α1∗=−2λ∗1+ξ1−λ∗2,α0∗=1−ξa1−λ∗21+ξ1−λ∗2,β1∗=−2λ∗1+ξ1−λ∗2,β0∗=1−ξ1−λ∗21+ξ1−λ∗2.


## 5. Simulation and Experiment Results

### 5.1. Simulation Results

In this subsection, simulation is carried out to verify the effectiveness of the proposed adaptive notch filter for resonance suppression. The transfer function for simulation model consists of three resonance modes (39 Hz, 60 Hz, and 65 Hz), which is given by(35)Gss=s2+2×0.01×2π×60s+2π×602s2+2×0.035×2π×39s+2π×392s2+2×0.05×2π×65s+2π×652.Figure7 shows convergence of estimated frequency parameters. It can be seen that large fluctuation only occurred in initial stage, and frequency parameters quickly converge to reference value in a short time. This indicates that resonance frequency detector is able to estimate the resonance efficiently. The amplitude-phase characteristic of simulation with and without the adaptive notch filter is shown in Figure 8. It is clear that the system has three resonance frequencies, which severally restricts the dynamic performance and stability performance of system. However, the three resonance frequencies are significantly attenuated by the adaptive filter notch, compared with the case of the system without the adaptive notch filter. Although the bandwidth of shaking table is sacrificed, the stability of system is well ensured, which is important for shaking table test. The detailed parameters of the adaptive notch filter are listed in Table 1. This shows that the proposed adaptive notch filter is very effective in suppressing the resonance modes.Figure 7
Convergence of estimated frequency parameters.Figure 8
Amplitude-phase characteristic of simulation.Table 1
The parameters of adaptive notch filter.
α
2
∗
α
1
∗
α
0
∗
β
1
∗
β
0
∗
39 Hz
0.86
1.92
0.97
1.92
0.60
60 Hz
0.91
1.42
0.56
1.42
0.99
65 Hz
0.98
1.84
0.94
1.84
0.94Figure9 shows the sine sweep frequency response with and without the adaptive notch filter. As can be seen, magnitude at the resonance frequencies obviously decreases using the proposed method. Figure 10 shows the acceleration response of 3sin80πtm/s2 excitation signal. It can be seen that the suppressed results have a good agreement with reference signal compared with unsuppressed results, which also proves that the proposed method has better performance on resonance suppression. Combined with the proposed adaptive notch filter, the final maximum tracking error is reduced from 1.8 m/s2 to less than 0.1 m/s2.Figure 9
The sine sweep frequency response of simulation.Figure 10
Sine acceleration response. (a) Tracking performance. (b) Tracking error.
(a)
(b)


### 5.2. Experimental Results

The electrohydraulic servo shake table test used to validate the proposed resonance suppression scheme is shown in Figure1, whose main parameters are listed in Table 2. The adaptive algorithm for the resonance frequency parameters estimation is depicted in Figure 11. Although there is an error at the beginning of the estimation process, the estimation of frequency is rapidly decreased and converges to steady state. Amplitude-phase characteristic of the shaking table system without the adaptive notch filter and with the adaptive notch filter is shown in Figure 12. The parameters calculation of the adaptive notch filter is similar with the simulation, not tried in words here. As can be seen, the generated three resonance modes are significantly attenuated using the proposed scheme. The gain margin and phase margin with the adaptive notch filter are greatly improved compared to without the adaptive notch filter, which is conductive to precise control of the hydraulic shaking table. The power spectrum of the response signal with and without the adaptive notch filter is shown in Figure 13. It is clear that amplitudes of the response signal are significantly attenuated at the resonance frequency by using the proposed resonance suppression scheme.Table 2
Main parameters of the hydraulic shaking table.
Parameters
Values
Flow gainKq
2
×
10
−
3
m
3
/
s
/
A
Flow-pressure coefficientKc
3
×
10
−
12
m
3
/
s
/
Pa
Piston’s effective areaA
1.9
×
10
−
3
m
2
Total leakage coefficientCt
9
×
10
−
17
m
3
/
s
/
Pa
Total chamber volumeVt
3.78
×
10
−
4
m
3
Bulk modulusβe
6
.
9
×
10
8
Pa
Viscous damping coefficientB
3.2
×
10
3
N
/
m
/
s
Spring stiffnessK
1.7
×
10
6
N
/
mFigure 11
Convergence of estimated frequency parameters.Figure 12
Amplitude-phase characteristic of shaking table system.Figure 13
Experiment spectrum of response signal.Figure14 shows the sine acceleration response of tracking performance and the tracking error with and without adaptive notch filter, where the excitation signal is 3sin6πtm/s2. It is obvious that sine acceleration response with the adaptive notch filter matches well with the reference signal compared with the situation without the adaptive notch filter. The tracking error of sine acceleration response reveals that the error has been greatly reduced by the adaptive notch filter from the initial maximum 1.5 m/s2 to 0.3 m/s2.Figure 14
Sine acceleration response. (a) Tracking performance. (b) Tracking error.
(a)
(b)Apart from sine acceleration excitation signal, a recorded earthquake wave, which happened in El Centro in southern California, is also used to validate the proposed resonance suppression scheme. The experiment results of tracking performance and the tracking error are presented in Figure15. As can be inferred from Figure 15, the proposed resonance suppression scheme with the adaptive notch filter has better tracking accuracy than without the adaptive notch filter.Figure 15
Experiment results of El Centro earthquake wave. (a) Tracking performance. (b) Tracking error.
(a)
(b)


## 5.1. Simulation Results

In this subsection, simulation is carried out to verify the effectiveness of the proposed adaptive notch filter for resonance suppression. The transfer function for simulation model consists of three resonance modes (39 Hz, 60 Hz, and 65 Hz), which is given by(35)Gss=s2+2×0.01×2π×60s+2π×602s2+2×0.035×2π×39s+2π×392s2+2×0.05×2π×65s+2π×652.Figure7 shows convergence of estimated frequency parameters. It can be seen that large fluctuation only occurred in initial stage, and frequency parameters quickly converge to reference value in a short time. This indicates that resonance frequency detector is able to estimate the resonance efficiently. The amplitude-phase characteristic of simulation with and without the adaptive notch filter is shown in Figure 8. It is clear that the system has three resonance frequencies, which severally restricts the dynamic performance and stability performance of system. However, the three resonance frequencies are significantly attenuated by the adaptive filter notch, compared with the case of the system without the adaptive notch filter. Although the bandwidth of shaking table is sacrificed, the stability of system is well ensured, which is important for shaking table test. The detailed parameters of the adaptive notch filter are listed in Table 1. This shows that the proposed adaptive notch filter is very effective in suppressing the resonance modes.Figure 7
Convergence of estimated frequency parameters.Figure 8
Amplitude-phase characteristic of simulation.Table 1
The parameters of adaptive notch filter.
α
2
∗
α
1
∗
α
0
∗
β
1
∗
β
0
∗
39 Hz
0.86
1.92
0.97
1.92
0.60
60 Hz
0.91
1.42
0.56
1.42
0.99
65 Hz
0.98
1.84
0.94
1.84
0.94Figure9 shows the sine sweep frequency response with and without the adaptive notch filter. As can be seen, magnitude at the resonance frequencies obviously decreases using the proposed method. Figure 10 shows the acceleration response of 3sin80πtm/s2 excitation signal. It can be seen that the suppressed results have a good agreement with reference signal compared with unsuppressed results, which also proves that the proposed method has better performance on resonance suppression. Combined with the proposed adaptive notch filter, the final maximum tracking error is reduced from 1.8 m/s2 to less than 0.1 m/s2.Figure 9
The sine sweep frequency response of simulation.Figure 10
Sine acceleration response. (a) Tracking performance. (b) Tracking error.
(a)
(b)


## 5.2. Experimental Results

The electrohydraulic servo shake table test used to validate the proposed resonance suppression scheme is shown in Figure1, whose main parameters are listed in Table 2. The adaptive algorithm for the resonance frequency parameters estimation is depicted in Figure 11. Although there is an error at the beginning of the estimation process, the estimation of frequency is rapidly decreased and converges to steady state. Amplitude-phase characteristic of the shaking table system without the adaptive notch filter and with the adaptive notch filter is shown in Figure 12. The parameters calculation of the adaptive notch filter is similar with the simulation, not tried in words here. As can be seen, the generated three resonance modes are significantly attenuated using the proposed scheme. The gain margin and phase margin with the adaptive notch filter are greatly improved compared to without the adaptive notch filter, which is conductive to precise control of the hydraulic shaking table. The power spectrum of the response signal with and without the adaptive notch filter is shown in Figure 13. It is clear that amplitudes of the response signal are significantly attenuated at the resonance frequency by using the proposed resonance suppression scheme.Table 2
Main parameters of the hydraulic shaking table.
Parameters
Values
Flow gainKq
2
×
10
−
3
m
3
/
s
/
A
Flow-pressure coefficientKc
3
×
10
−
12
m
3
/
s
/
Pa
Piston’s effective areaA
1.9
×
10
−
3
m
2
Total leakage coefficientCt
9
×
10
−
17
m
3
/
s
/
Pa
Total chamber volumeVt
3.78
×
10
−
4
m
3
Bulk modulusβe
6
.
9
×
10
8
Pa
Viscous damping coefficientB
3.2
×
10
3
N
/
m
/
s
Spring stiffnessK
1.7
×
10
6
N
/
mFigure 11
Convergence of estimated frequency parameters.Figure 12
Amplitude-phase characteristic of shaking table system.Figure 13
Experiment spectrum of response signal.Figure14 shows the sine acceleration response of tracking performance and the tracking error with and without adaptive notch filter, where the excitation signal is 3sin6πtm/s2. It is obvious that sine acceleration response with the adaptive notch filter matches well with the reference signal compared with the situation without the adaptive notch filter. The tracking error of sine acceleration response reveals that the error has been greatly reduced by the adaptive notch filter from the initial maximum 1.5 m/s2 to 0.3 m/s2.Figure 14
Sine acceleration response. (a) Tracking performance. (b) Tracking error.
(a)
(b)Apart from sine acceleration excitation signal, a recorded earthquake wave, which happened in El Centro in southern California, is also used to validate the proposed resonance suppression scheme. The experiment results of tracking performance and the tracking error are presented in Figure15. As can be inferred from Figure 15, the proposed resonance suppression scheme with the adaptive notch filter has better tracking accuracy than without the adaptive notch filter.Figure 15
Experiment results of El Centro earthquake wave. (a) Tracking performance. (b) Tracking error.
(a)
(b)


## 6. Conclusions

In this paper, the operational principle and hydraulic actuator model are formulated to illustrate the resonance problem existing in the electrohydraulic servo shake table. The coupling between shaking table itself and the specimen’s elastic load deteriorates the system’s dynamic performance and even leads to the instability of the control system. Furthermore, the resonance frequency is a time-varying value during the shaking table test. In order to meet the requirements of resonance online suppressing, the adaptive notch filter based on the LMS adaptive algorithm is designed to online identify the resonance frequency and automatically adjust the filter coefficients. Simulation and experiment results demonstrate that the proposed method can be effectively used in suppressing the varying resonance mode.



---

*Source: 9407520-2019-05-21.xml*