# Study on the Resonant Suppression of Electrohydraulic Servo Shaking Table of a 2-Mass Dynamic System

Authors: Tao Wang; Jinchun Song
Journal: Shock and Vibration
 (2019)

---

## Abstract

Since the load of electrohydraulic servo shaking table is an elastic load, there is a mutual coupling effect that exists between the experimental object and the shaking table, forming the resonance to weaken the dynamic characteristics, and producing resonant peak and resonant valley in the bandwidth required by the system, in which the amplitude is often larger than the stability range of acceleration’s amplitude. In this paper, the mathematical modeling of hydraulic power mechanism in a two-mass dynamic system is established based on electrohydraulic servo shaking table, yielding the frequency characteristic curve in accordance with the transfer function of the model. A multifrequency adaptive notch filter based on the least mean square algorithm is proposed to suppress the resonance, and the suppression effect of the resonance is simulated and verified in line with various values of load stiffness. Finally, the power spectrum is used to demonstrate the effectiveness of the multifrequency adaptive notch filter in the resonant suppression.


---

## Body

## 1. Introduction

In recent years, the shaking table has become a critical test tool in engineering research and plays a crucial role in environmental simulation, e.g., aerospace, automobile, building, and earthquake simulation, especially in modern times, which requires high reliability of the system. The hydraulic vibration test also has become an essential way of modern industrial reliability test. With the improvement of product reliability and the standard of vibration damping, the servo shaking table directly impacts the accuracy of test. Thus, the system’s reliability is becoming increasingly demanding [1, 2].When the elastic load test is performed on the shaking table, the oil liquid in the closed chamber is compressed, and the leakage is equivalent to a “hydraulic spring.” In the meantime, the shaking table test system has a coupling relation with the elastic load, reducing the dynamic characteristics of the system, and a corresponding resonant peak and the resonant valley will exist in the bandwidth required by the system [3, 4]. Due to the existence of this resonant phenomenon, the stability of the system becomes worse, which cannot meet the requirements of system stability. For the vibration test system, how to eliminate the coupling effect on the frequency characteristic of the shaking table’s acceleration, to restrain the resonance and make the frequency characteristic of the shaking table meet the standard of vibration test, has become a critical subject of the present research.The existence of the resonance will not only impact the working performance and the service life of the components but also produce noise and other types of pollution and sometimes even damage the components and sealing devices. The harmonic oscillation will occur if the natural frequency of the hydraulic system happens to the frequency of the pressure oscillation simultaneously. This type of oscillation suddenly strikes in the short time which will not only be reduced by the influence of resonance, whereas it will in turn increase. Under this short-time impact, the pipe, component, and sealing device of the hydraulic system will be higher than the working pressure and make them to be seriously damaged to affect normal work after several times of pressure shock. Besides, the impact of pressure shock and oscillation also lead to noise pollution in the workspace [5, 6].In comparison with the mechanical system and the power system, the resonant suppression method based on the servo system is quite limited by its high requirement and precision. Thus, the new suppression methods should adapt to meet the performance requirements of electrohydraulic servo system. Resonance of the brake pressure control system is a fluid-solid coupling vibration caused by the joint action of the brake control valve and the brake system pipe. Given this, Lu et al. put forward the method of changing the parameters of the pipeline, the back pressure of the system, and the brake control valve in the brake system, capable of effectively restraining the resonant phenomenon of the aircraft brake system [7]. Fu et al. had analyzed the low-frequency resonance of an aero hydraulic engine and designed a feedforward controller composed of low pass filter, PD compensation, and feedforward correction, which successfully suppressed the resonance and improved the overall amplitude and phase margin of the system [8]. Tomoki et al. proposed a parameter identification method to model electrohydrostatic actuators as two-mass resonant systems, and the self-resonant cancellation technique was implemented to suppress the effect of resonance [9]. The notch filters were applied by Younkin in fixed frequency resonances to stabilize the servo drives with unacceptable machine dynamics, and this frequency selective feedback control technique also stabilized the servo drives with unacceptable machine dynamics [10]. Cao et al. put forward a method to suppress the resonance of the excited cylinder by setting an external damper between the two working chambers of the excited cylinder to form a bypass leakage to improve the hydraulic damping ratio [11].In the servo system, change of natural frequency, load characteristic, and additional characteristic can change the frequency and amplitude of the resonance, so the adaptive algorithm is put forward to suppress the resonance. A developed method for resonant suppression based on adaptive filteredx-least mean square algorithm was presented by Wang et al., which obtained the resonance iteratively and reduced the significant resonant oscillations using a finite impulse response filter, which was adapted by the least mean square error principle [12]. Kang and Kim proposed the method of an adaptive digital notch filter that could identify the resonant frequency of the actuator quickly and adjust automatically its center frequency [13]. The resonant online suppressing algorithm based on adaptive IIR notch filter was proposed by Li et al. in accordance with the application requirements; the resonant frequency of the system was online estimated and then the frequency parameters of notch filter was adjusted online by using a recursive least-squares algorithm [14]. Bahn et al. brought a new algorithm for estimating the resonant frequency of adaptive notch filters in servo systems, which was proportional to the difference in the estimated frequency and actual frequency and allowed to select the estimation parameters independent from the value of resonant frequencies [15]. To improve the performance of the servo system to overcome the influence of the previously neglected elastic components, the adaptive notch filter was employed by Yang et al. for online resonant suppression; thus, the characteristics of the mechanical resonance could be identified quickly, and the parameters of the notch filter could be automatically adjusted by the identification results [16].In brief, though there are numerous methods to suppress the resonance, the suppression of the experimental platform based on shaking table is not much, which should keep the real-time suppression effect when this method is applied in the servo system. The rest of the paper is arranged as follows: in Section2, the frequency characteristics of the shaking table are provided based on the power mechanism. In Section 3, the theory of adaptive notch filter methods for resonance suppression is proposed, and the multifrequency adaptive notch filter is analyzed based on LMS. In Section 4, the simulation results of the resonance suppression in line with different values of load stiffness and the power spectrum of the electrohydraulic servo system are presented. The effectiveness and superiority are verified by the comparison of the suppression resonant effect between the multifrequency adaptive notch filter and the digital filter. At last, Section 5 is the conclusion.


## 2. Electrohydraulic Servo Shaking Table

### 2.1. The Equipment

The block diagram of the shaking table system is shown in Figure1. It is a uniaxial shaking table with supply pressure of 8 MPa. A hydraulic cylinder is controlled by a two-stage servo valve. The stroke of the cylinder is ±25 mm, and the diameters of its rod and piston are Ф36 mm and Ф50 mm, respectively. The system is under acceleration control with a frequency range from 0 Hz to 125 Hz to reproduce the desired shaking response.Figure 1
Control system of shaking table.The three A/D channels are mainly used to receive displacement, acceleration signals, and speed signals obtained by conditioning, and all of them are converted into digital signals and transmitted to the upper computer for analysis and processing in the system. A D/A channel is used to convert the control signals of the system into analog signals and input them into the shaking table system, so as to realize the conversion between the signals of the system [17].


### 2.2. Analysis of the Hydraulic Power Mechanism

In Figure2, the load part is abstracted into an elastic load, and the effect of friction on the load is ignored. Under the action of elastic load, the electrohydraulic servo system is a two-degree-of-freedom model [18].Figure 2
Schematic diagram of hydraulic power mechanism.Before modeling the servo system, it is necessary to make some following assumptions: without considering leakage and compressibility of oil, the zero-opening four-side slide valves should be used as the type, in which throttles of the hydraulic can be matched with each other, and they are completely in symmetrical distribution. Because the constant pressure variable oil source is adopted, the oil pressure is kept constant. So, the oil return pressure is at zero [19]. The equation is shown as follows:(1)qL=Kqxv−KcpL.The following assumptions are shown: connection among the components of the servo system remains symmetrical; it does not consider the effect of friction in the hydraulic pipeline [20]. Under this condition, the equation is established as follows:(2)qL=Adypdt+CtppL+Vt4βe⋅dpLdt.The right-hand side of Equation (2) is made up of three parts. One is the flow used to act on the piston rod which plays a leading role in providing power to the hydraulic cylinder for the action of the piston rod. The second is the total leakage flow of the system that consists mainly of two components of internal leakage and external leakage, which form a hydraulic servo system. The third is the total compression flow, which is manifested in the compressibility of hydraulic oil to composing the compression loss of the system, but it can be ignored.The equation of force balance of two degrees of freedom is established, and its characteristics are analyzed by linearization method. Considering the influence of load force, viscous damping force, and inertia force, the force balance equation can be obtained as follows:(3)ApL=mpd2ypdt2+BLdypdt+KLyP−BLdyLdt−KLyL,(4)mLd2yLdt2=BLdypdt+KLyP−BLdyLdt−KLyL.Equations (1)–(4) can be organized by Laplace transformation:(5)YLXV=KqABLs+KLDs,(6)Ds=mpmLKhs5+BLmp+mLKh+KcempmLA2s4+mL+KceBLmp+mLA2+KLmp+mLKhs3+BL+KceKLmp+mLA2s2+KLs.Using equations (3) and (4), we obtain the following equation:(7)YLYp=BLs+KLmLs2+BLs+KL.Using equations (5) and (7), equation (8) can be summarized as follows:(8)YpXV=KqAmLs2+BLs+KLDs.From the above, it can be seen that when the vibration test of the elastic load is carried out, the coupling action of the hydraulic system constitutes a resonant system. The comprehensive stiffness of the system is made up of two parts in series, one is the stiffness of the hydraulic springKh, which is caused by the compressibility of the fluid. The second is the load stiffness KL, which is determined by the load condition of the load. After being connected in series, the composite stiffness is less than Kh and KL.Comparison of equation (8) shows that the elastic system increases two links than the inertial system, one is the oscillation link which can cause the emergence of the resonant peak, and the other is the two-order differential link which can cause the appearance of the resonant valley. Therefore, equation (8) causes the resonance induced in the system bandwidth, and the amplitude of the resonance exceeds the required range due to the smaller damping. The conventional resonant suppression strategy, such as improving the comprehensive resonant frequency and the comprehensive damping ratio, is difficult to obtain accurate coupling parameters because of the uncertainty and nonlinearity of the coupling characteristics. The suppression effect of the resonance is often very limited, which cannot meet the demand of tolerance limit. Therefore, it is necessary to study the nonlinear control strategy with adaptive function to suppress the coupling effect of elastic load on the acceleration frequency characteristics of shaking table.The effective working area of the piston rodA is calculated as follows:(9)A=πD2−d24=945.6×10−6m2,where d represents the diameter of the piston rod and D represents the inner diameter of the hydraulic cylinder. L value is ±200 mm which represents the hydraulic cylinder in stroke. The total volume of the oil tank and the pipeline oil can be obtained as(10)Vt=AL=3.78×10−4m3.Compressed liquid of hydraulic cylinder can be approximately a linear equivalent total hydraulic spring stiffness, soKh is(11)Kh=βeA21Vi+1Vo.And when the piston rod is at the middle of the position, the value is(12)Kh=4βeA2Vt=6.53×106N/m,where the damping ratio is the smallest value. So the parameters based on the hydraulic system power mechanism are shown as Table 1.Table 1
Parameters based on the hydraulic system power mechanism.
Piston diameterD (mm)
50
Piston rod diameterd (mm)
36
Piston effective areaA (m2)
945.6 × 10−6
Stroke of hydraulic cylinderL (mm)
±200
Equivalent volume modulus of elasticityβe (MPa)
690
Pipe diameterd1 (m)
0.01
Equivalent mass of load inertiamp (kg)
100
Elastic load qualitymL (kg)
100
Hydraulic spring stiffnessKh (N/m)
6.53 × 106
Servo valve gainKq (m3/s/A)
2 × 10−2
Flow-pressure coefficientKce (m3/s/Pa)
3 × 10−12
Viscous damping coefficient of the piston and loadBL (N/(m/s))
80


### 2.3. The Frequency Analysis

From the above section, it can be seen that the elastic system increases the oscillating link and the two-order differential link compared with the inertial system, which leads to the emergence of resonance. In order to determine the relationship between the frequency value of the resonant peak and valley and the natural frequency of the dynamic mechanism, the two-degree-of-freedom equivalent vibration model is established that ignores the damping effect, and the dynamic differential equations are established, respectively:(13)mpd2ypdt2+KL+KhyP−KLyL=0,(14)mLd2yLdt2−KLyP+KLyL=0.Equations (13) and (14) can be expressed in matrix form as(15)mp00mlypyL+Kh+KL−KL−KLKLypyL=00.The natural frequency of the two-degree-of-freedom system can be obtained by solving the following characteristic equation.(16)Kh+KL−KL−KLKL−ω2mp00ml=0.The expansions of equation (16) can be obtained as follows:(17)mLmpω4−KhmL+KLmL+KLmpω2+KhKL=0.Thus, two natural frequencies can be solved, respectively.(18)ω02=KhmL+KLmL+KLmp−KhmL+KLmL+KLmp2−4KhKLmLmp2mLmp,ωp2=KhmL+KLmL+KLmp+KhmL+KLmL+KLmp2−4KhKLmLmp2mLmp,where ω0 is the composite frequency and ωp is the frequency of resonant peak. Equation (6) can be converted into the following form:(19)Ds=1ωL2ωh2s5+2ςLωLωx+ςxωxωL2ωh2s4+1ωL2+2ςxςLωxωL+1ωx2s3+2ςLωL+ςxωxs2+s,where ωL is the frequency of the resonant valley, ωL=KL/mL, ςL=BL/2KLmL, ωx is the coupling frequency, ωx=Kh/mp+mL, ςx=2Kce/Aβe/Vtmp+mL, and ωh is the natural frequency of the power mechanism, ωh=Kh/mp.But the actual value is not exactly equal to the calculated value. The reason is that the damping of the hydraulic system is small, so the damping has little influence on the comprehensive natural frequency, which can be used to calculate the integrated fixed frequency.


## 2.1. The Equipment

The block diagram of the shaking table system is shown in Figure1. It is a uniaxial shaking table with supply pressure of 8 MPa. A hydraulic cylinder is controlled by a two-stage servo valve. The stroke of the cylinder is ±25 mm, and the diameters of its rod and piston are Ф36 mm and Ф50 mm, respectively. The system is under acceleration control with a frequency range from 0 Hz to 125 Hz to reproduce the desired shaking response.Figure 1
Control system of shaking table.The three A/D channels are mainly used to receive displacement, acceleration signals, and speed signals obtained by conditioning, and all of them are converted into digital signals and transmitted to the upper computer for analysis and processing in the system. A D/A channel is used to convert the control signals of the system into analog signals and input them into the shaking table system, so as to realize the conversion between the signals of the system [17].


## 2.2. Analysis of the Hydraulic Power Mechanism

In Figure2, the load part is abstracted into an elastic load, and the effect of friction on the load is ignored. Under the action of elastic load, the electrohydraulic servo system is a two-degree-of-freedom model [18].Figure 2
Schematic diagram of hydraulic power mechanism.Before modeling the servo system, it is necessary to make some following assumptions: without considering leakage and compressibility of oil, the zero-opening four-side slide valves should be used as the type, in which throttles of the hydraulic can be matched with each other, and they are completely in symmetrical distribution. Because the constant pressure variable oil source is adopted, the oil pressure is kept constant. So, the oil return pressure is at zero [19]. The equation is shown as follows:(1)qL=Kqxv−KcpL.The following assumptions are shown: connection among the components of the servo system remains symmetrical; it does not consider the effect of friction in the hydraulic pipeline [20]. Under this condition, the equation is established as follows:(2)qL=Adypdt+CtppL+Vt4βe⋅dpLdt.The right-hand side of Equation (2) is made up of three parts. One is the flow used to act on the piston rod which plays a leading role in providing power to the hydraulic cylinder for the action of the piston rod. The second is the total leakage flow of the system that consists mainly of two components of internal leakage and external leakage, which form a hydraulic servo system. The third is the total compression flow, which is manifested in the compressibility of hydraulic oil to composing the compression loss of the system, but it can be ignored.The equation of force balance of two degrees of freedom is established, and its characteristics are analyzed by linearization method. Considering the influence of load force, viscous damping force, and inertia force, the force balance equation can be obtained as follows:(3)ApL=mpd2ypdt2+BLdypdt+KLyP−BLdyLdt−KLyL,(4)mLd2yLdt2=BLdypdt+KLyP−BLdyLdt−KLyL.Equations (1)–(4) can be organized by Laplace transformation:(5)YLXV=KqABLs+KLDs,(6)Ds=mpmLKhs5+BLmp+mLKh+KcempmLA2s4+mL+KceBLmp+mLA2+KLmp+mLKhs3+BL+KceKLmp+mLA2s2+KLs.Using equations (3) and (4), we obtain the following equation:(7)YLYp=BLs+KLmLs2+BLs+KL.Using equations (5) and (7), equation (8) can be summarized as follows:(8)YpXV=KqAmLs2+BLs+KLDs.From the above, it can be seen that when the vibration test of the elastic load is carried out, the coupling action of the hydraulic system constitutes a resonant system. The comprehensive stiffness of the system is made up of two parts in series, one is the stiffness of the hydraulic springKh, which is caused by the compressibility of the fluid. The second is the load stiffness KL, which is determined by the load condition of the load. After being connected in series, the composite stiffness is less than Kh and KL.Comparison of equation (8) shows that the elastic system increases two links than the inertial system, one is the oscillation link which can cause the emergence of the resonant peak, and the other is the two-order differential link which can cause the appearance of the resonant valley. Therefore, equation (8) causes the resonance induced in the system bandwidth, and the amplitude of the resonance exceeds the required range due to the smaller damping. The conventional resonant suppression strategy, such as improving the comprehensive resonant frequency and the comprehensive damping ratio, is difficult to obtain accurate coupling parameters because of the uncertainty and nonlinearity of the coupling characteristics. The suppression effect of the resonance is often very limited, which cannot meet the demand of tolerance limit. Therefore, it is necessary to study the nonlinear control strategy with adaptive function to suppress the coupling effect of elastic load on the acceleration frequency characteristics of shaking table.The effective working area of the piston rodA is calculated as follows:(9)A=πD2−d24=945.6×10−6m2,where d represents the diameter of the piston rod and D represents the inner diameter of the hydraulic cylinder. L value is ±200 mm which represents the hydraulic cylinder in stroke. The total volume of the oil tank and the pipeline oil can be obtained as(10)Vt=AL=3.78×10−4m3.Compressed liquid of hydraulic cylinder can be approximately a linear equivalent total hydraulic spring stiffness, soKh is(11)Kh=βeA21Vi+1Vo.And when the piston rod is at the middle of the position, the value is(12)Kh=4βeA2Vt=6.53×106N/m,where the damping ratio is the smallest value. So the parameters based on the hydraulic system power mechanism are shown as Table 1.Table 1
Parameters based on the hydraulic system power mechanism.
Piston diameterD (mm)
50
Piston rod diameterd (mm)
36
Piston effective areaA (m2)
945.6 × 10−6
Stroke of hydraulic cylinderL (mm)
±200
Equivalent volume modulus of elasticityβe (MPa)
690
Pipe diameterd1 (m)
0.01
Equivalent mass of load inertiamp (kg)
100
Elastic load qualitymL (kg)
100
Hydraulic spring stiffnessKh (N/m)
6.53 × 106
Servo valve gainKq (m3/s/A)
2 × 10−2
Flow-pressure coefficientKce (m3/s/Pa)
3 × 10−12
Viscous damping coefficient of the piston and loadBL (N/(m/s))
80


## 2.3. The Frequency Analysis

From the above section, it can be seen that the elastic system increases the oscillating link and the two-order differential link compared with the inertial system, which leads to the emergence of resonance. In order to determine the relationship between the frequency value of the resonant peak and valley and the natural frequency of the dynamic mechanism, the two-degree-of-freedom equivalent vibration model is established that ignores the damping effect, and the dynamic differential equations are established, respectively:(13)mpd2ypdt2+KL+KhyP−KLyL=0,(14)mLd2yLdt2−KLyP+KLyL=0.Equations (13) and (14) can be expressed in matrix form as(15)mp00mlypyL+Kh+KL−KL−KLKLypyL=00.The natural frequency of the two-degree-of-freedom system can be obtained by solving the following characteristic equation.(16)Kh+KL−KL−KLKL−ω2mp00ml=0.The expansions of equation (16) can be obtained as follows:(17)mLmpω4−KhmL+KLmL+KLmpω2+KhKL=0.Thus, two natural frequencies can be solved, respectively.(18)ω02=KhmL+KLmL+KLmp−KhmL+KLmL+KLmp2−4KhKLmLmp2mLmp,ωp2=KhmL+KLmL+KLmp+KhmL+KLmL+KLmp2−4KhKLmLmp2mLmp,where ω0 is the composite frequency and ωp is the frequency of resonant peak. Equation (6) can be converted into the following form:(19)Ds=1ωL2ωh2s5+2ςLωLωx+ςxωxωL2ωh2s4+1ωL2+2ςxςLωxωL+1ωx2s3+2ςLωL+ςxωxs2+s,where ωL is the frequency of the resonant valley, ωL=KL/mL, ςL=BL/2KLmL, ωx is the coupling frequency, ωx=Kh/mp+mL, ςx=2Kce/Aβe/Vtmp+mL, and ωh is the natural frequency of the power mechanism, ωh=Kh/mp.But the actual value is not exactly equal to the calculated value. The reason is that the damping of the hydraulic system is small, so the damping has little influence on the comprehensive natural frequency, which can be used to calculate the integrated fixed frequency.


## 3. Multifrequency Adaptive Notch Filter

Adaptive notch technology is derived from adaptive noise cancellation, which sets a sinusoidal signal as a reference signal to counteract every component in narrowband noise. When the sine wave is used as the reference input, it can eliminate the main spectrum component at the center of the narrow band near the reference frequency. Because the adaptive notch filter has only one parameter to estimate, the algorithm is simple and easy to control the bandwidth and the depth of zero point, which can accurately track the frequency of interference, and its the stability is guaranteed by limiting the pole to the unit circle [21].The block diagram of the single frequency adaptive notch filter has 2 adaptive weights as shown in Figure3. The original input is a mixed waveform s+n of the real signal s and the monochromatic interference n. Besides, the noise and signal are irrelevant. The reference input is a signal x which is not related to the signal, but it is related to noise. The original input is transferred to the dk end by sampling; the reference input x is sampled to the x1k end and the x2k end. The x2k is a sample value after 90-phase shift of the reference input x, and the purpose is to obtain the two weights of w1k and w2k. Accordingly, the amplitude of the sine wave after the combination can be equal to the amplitude of the interference components in the original input. The weighted output of reference input is yk, i.e., the estimation of noise nk, containing the sampling signals dk and yk interfering with the frequency. After subtraction, the optimal estimate ek of useful signals is obtained.Figure 3
Single frequency adaptive notch filter.Assuming that the input signal is arbitrary, the reference input is a single frequency sine wave which is expressed asCcosΩ0t+φ, where Ω0 in the formula is the angular frequency of the analog signal and C is the amplitude of the reference signal. The original input and the reference input are periodically sampled by the time interval T, while ω0=Ω0T is the sampled signal digital angular frequency. Using the LMS algorithm, the single frequency adaptive notch filter transfer function can be obtained as follows:(20)Hz=z2−2zcosω0+1z2−21−μC2zcosω0+1−2μC2.There is zero at frequencyf0, and it is exactly located at z=e±jω0 of unit flat of Z plane, while the pole is located as follows:(21)z=1−μC2cosω0±j1−μC2−1−μC22cos2ω01/2.The poles in the unit circle, the radial distance of which from the origin is1−μC21/2, approximately equals to 1−μC2 so the system is stable. The multifrequency adaptive notch filter is implemented under the parallel structure so that the original input remains to be disturbed signal, and the reference input is expanded into multiple with the amplitude adjusted by 2 weights to realize multiple frequency notch. The output of the adaptive notch filter is the sum of the output of the two-stage notch filter, and the output of the noise canceller is the difference between the original signal input and the output of the adaptive filter. Figure 4 shows that the multifrequency adaptive notch filter is the first notch filter output as the input of the second notch filter, forming a cascade [22].Figure 4
Multifrequency adaptive notch filter.(22)
H
z
=
H
0
z
H
1
z
,
H
0
z
=
z
2
−
2
z
cos
ω
0
+
1
z
2
−
2
1
−
μ
0
C
0
2
z
cos
ω
0
+
1
−
2
μ
0
C
0
2
,
H
1
z
=
z
2
−
2
z
cos
ω
1
+
1
z
2
−
2
1
−
μ
1
C
1
2
z
cos
ω
1
+
1
−
2
μ
1
C
1
2
.The zero poles of the system are the zero poles of each notch filter, which can directly change the adjustment parameters of one level without affecting the zero pole characteristics of the other notch filter. Therefore, the cascade structure not only has the advantages of the single frequency notch filter but also makes the performance of the notch filter simple and effective and is easier to realize.


## 4. Simulation Results

In this paper, the frequency response curve of electrohydraulic servo shaking table is measured by using the random wave signal generated by the signal generator. The diagram of resonant suppression is given in Figure5. The input signal generates a sinusoidal signal with a frequency of resonant frequency ω, and the error signal is obtained under the difference between acceleration signal and acceleration response signal of the electrohydraulic servo shaking table. The LMS adaptive filtering algorithm adjusts the weight of the adaptive notch filter online using the error signal, and the input signal of the adaptive notch filter is weighted by its weight to obtain the output signal. The output signal of the adaptive notch filter serves as the restraining quantity of the resonance. Moreover, the signal added with the input signal serves as new input signal of the shaking table system. After several iterations, the weight of the adaptive notch filter can significantly suppress the resonance, and its amplitude can also be within the required range.Figure 5
Principle of suppression of the resonance.Since the resonance includes resonant peak and resonant valley, the adaptive notch filter must be employed to suppress two targets simultaneously. In this case, it is necessary to input two reference sinusoidal signals, in which the frequency of sinusoidal signals is the frequency of resonant peak and resonant valley. Subsequently, the resonance can be suppressed using the schematic diagram of the resonant suppression, as shown in Figure5. And the models could be obtained in accordance with previous equation to analyze the performance. The Bode diagram of the open-loop transfer function can be simulated by the noted parameters as shown in Figure 6.Figure 6
The Bode diagram of the open-loop transfer function under the elastic load.As shown in Figure7, the low-frequency segment is integral link and two-order oscillation link. As ωL > ω0, the two-order oscillating link is followed by the two-order differential link. The two-order differential link makes the transfer function appear an resonant valley near the natural frequency ωL of the flexible foundation. And the natural frequency has a cancellation effect on the comprehensive frequency ω0, which leads to the reduction at ω0. As ωL < ωp, two-order differential link is followed by a two-order oscillation link. In brief, due to the impact of the elastic load, the frequency of the resonant peak and the resonant valley of the dynamic mechanism are larger than the comprehensive frequency ω0, which affects the stability of the system. Now given the known Kh, resonant patterns of the system could be used under the different KL, and the result which is influenced by the multifrequency adaptive notch filter is verified according to the frequency characteristic diagram.


### 4.1.KL Is Less than Kh

The load stiffness is relatively small, and the electrohydraulic servo system is mainly based on rigid load. Under the condition of knownKh, the Bode diagram before resonant suppression is shown in Figure 7(a). It is known that resonance occurs at frequencies of 15.9 Hz and 43.8 Hz, respectively. Given the condition of known resonant frequency, the Bode diagram after resonant suppression can be obtained as shown in Figure 7(b).Figure 7
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant valley suppression.
(a)
(b)


### 4.2.KL Is Similar to Kh

As shown in Figure8(a) before resonant suppression, resonance is produced at the frequency of 37.4 Hz and 62.1 Hz, and the stability of the system begins to become worse because of the existence of resonance. For the test system, it is necessary to eliminate the coupling effect of the elastic load on the frequency characteristic of the shaking table’s acceleration and suppress the resonance, so that the frequency characteristic of the shaking table can be in accordance with the performance requirements. The frequency characteristic diagram after resonant suppression is shown in Figure 8(b).Figure 8
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant suppression.
(a)
(b)


### 4.3.KL Is Larger than Kh

NowKL becomes dominant, and the resonant frequency is larger than before. Figure 9(a) shows the frequency characteristic diagram before resonant suppression for which the frequency points are 50.3 Hz and 77.6 Hz, and Figure 9(b) shows the frequency characteristic diagram after resonant suppression.Figure 9
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant suppression.
(a)
(b)Finally, the influence of viscous damping coefficient based on the load of the open-loop transfer function is analyzed, and various damping values are applied for simulation analysis. The values ofBL are 100, 500, and 1000, and other parameters are selected in line with Table 1. As shown in Figure 10 which is the open-loop transfer function, the increase of damping does not impact the comprehensive natural frequency of ω0, and the frequency value does not change although the value of resonant peak and resonant valley decreases. Accordingly, increasing the viscous damping coefficient of the load can effectively increase the comprehensive damping ratio of the vibration system, whereas it cannot restrain the impact of the resonance.Figure 10
The open-loop Bode diagram of different viscous damping coefficients of load.In this paper, a variety of random vibration power spectrum control strategies are implemented by MATLAB software. In order to compare with others, the same reference power spectrum is used to assess the actual effects of various control techniques including stability, rapidity, accuracy, and other statistical properties. By using the power spectrum, the frequency characteristics should be seen more intuitively.According to reference power spectrum, digital notch filter and multifrequency adaptive notch filter are used, respectively, to realize the antiresonance. Figure11 shows that the frequency of resonance is 16.6 Hz and 44.3 Hz. Because of the resonant effect, the frequency response is not satisfactory. Its result is suppressed using digital filter depicted in Figure 12. Figure 13 is the resonant suppression using multifrequency adaptive notch filter. The comparison of the two methods shows that the traditional digital notch filter can ensure precision of control process, but small feedback gain leads to a relatively slow convergence process, and the excessive feedback gain easily leads to the divergence of the control process. The accuracy of the traditional digital notch filter can only reach the tolerance requirement of ±3 dB, and the maximum deviation occurs at the resonance point, but the multifrequency adaptive notch filter achieves an accuracy of ±3 dB after several iterations. The main reason for improving the precision is that the adaptive control can allow the increase of the feedback gain no matter what causes the instability of the system. That is to say, the adaptive control algorithm can keep the stability of the iterative process and also has control precision.Figure 11
Before resonant suppression.Figure 12
After resonant suppression by the digital filter.Figure 13
After resonant suppression by the multifrequency adaptive notch filter.


## 4.1.KL Is Less than Kh

The load stiffness is relatively small, and the electrohydraulic servo system is mainly based on rigid load. Under the condition of knownKh, the Bode diagram before resonant suppression is shown in Figure 7(a). It is known that resonance occurs at frequencies of 15.9 Hz and 43.8 Hz, respectively. Given the condition of known resonant frequency, the Bode diagram after resonant suppression can be obtained as shown in Figure 7(b).Figure 7
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant valley suppression.
(a)
(b)


## 4.2.KL Is Similar to Kh

As shown in Figure8(a) before resonant suppression, resonance is produced at the frequency of 37.4 Hz and 62.1 Hz, and the stability of the system begins to become worse because of the existence of resonance. For the test system, it is necessary to eliminate the coupling effect of the elastic load on the frequency characteristic of the shaking table’s acceleration and suppress the resonance, so that the frequency characteristic of the shaking table can be in accordance with the performance requirements. The frequency characteristic diagram after resonant suppression is shown in Figure 8(b).Figure 8
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant suppression.
(a)
(b)


## 4.3.KL Is Larger than Kh

NowKL becomes dominant, and the resonant frequency is larger than before. Figure 9(a) shows the frequency characteristic diagram before resonant suppression for which the frequency points are 50.3 Hz and 77.6 Hz, and Figure 9(b) shows the frequency characteristic diagram after resonant suppression.Figure 9
Frequency characteristic diagrams of the system: (a) before resonant suppression; (b) after resonant suppression.
(a)
(b)Finally, the influence of viscous damping coefficient based on the load of the open-loop transfer function is analyzed, and various damping values are applied for simulation analysis. The values ofBL are 100, 500, and 1000, and other parameters are selected in line with Table 1. As shown in Figure 10 which is the open-loop transfer function, the increase of damping does not impact the comprehensive natural frequency of ω0, and the frequency value does not change although the value of resonant peak and resonant valley decreases. Accordingly, increasing the viscous damping coefficient of the load can effectively increase the comprehensive damping ratio of the vibration system, whereas it cannot restrain the impact of the resonance.Figure 10
The open-loop Bode diagram of different viscous damping coefficients of load.In this paper, a variety of random vibration power spectrum control strategies are implemented by MATLAB software. In order to compare with others, the same reference power spectrum is used to assess the actual effects of various control techniques including stability, rapidity, accuracy, and other statistical properties. By using the power spectrum, the frequency characteristics should be seen more intuitively.According to reference power spectrum, digital notch filter and multifrequency adaptive notch filter are used, respectively, to realize the antiresonance. Figure11 shows that the frequency of resonance is 16.6 Hz and 44.3 Hz. Because of the resonant effect, the frequency response is not satisfactory. Its result is suppressed using digital filter depicted in Figure 12. Figure 13 is the resonant suppression using multifrequency adaptive notch filter. The comparison of the two methods shows that the traditional digital notch filter can ensure precision of control process, but small feedback gain leads to a relatively slow convergence process, and the excessive feedback gain easily leads to the divergence of the control process. The accuracy of the traditional digital notch filter can only reach the tolerance requirement of ±3 dB, and the maximum deviation occurs at the resonance point, but the multifrequency adaptive notch filter achieves an accuracy of ±3 dB after several iterations. The main reason for improving the precision is that the adaptive control can allow the increase of the feedback gain no matter what causes the instability of the system. That is to say, the adaptive control algorithm can keep the stability of the iterative process and also has control precision.Figure 11
Before resonant suppression.Figure 12
After resonant suppression by the digital filter.Figure 13
After resonant suppression by the multifrequency adaptive notch filter.


## 5. Conclusions

In this paper, when the shaking table is elastic load, the coupling action between elastic load and the hydraulic system constitutes an elastic load-resonant system. In this case, the servo system is analyzed by establishing the mathematical model, and the bandwidth can be improved by adding three-variable controller. To suppress the resonance of the shaking table at the frequency of the trap, we propose a multifrequency adaptive notch filter suppression strategy based on the LMS algorithm by studying the resonant frequency and the frequency characteristic of the adaptive trap at the zero point in the unit circle. Considering the different load stiffness, MATLAB simulates the resonant suppression and verifies the effectiveness of multifrequency adaptive notch filter.In this study, the random vibration control is conducted, and the results suggest that the multifrequency adaptive notch filter has better stability, quicker convergence speed, and higher precision than traditional notch filter. Also, the multifrequency adaptive notch filter meets requirements of the stability of the random process. Besides, it is suggested from the results that this control mode has high control accuracy and real-time performance for different frequency bands and different system characteristics and has high application value.



---

*Source: 8592578-2019-01-02.xml*