Multisupport, multidimension, and nonuniform excitation seismic experiments have new requirements for shaking table array system in synchronous tracking control. Therefore, this article proposed a novel synchronous tracking strategy, differential movement synchronous tracking control (DMSTC) strategy, for double-shaking table system while taking the interaction between shaking tables and specimen into consideration. DMSTC Simulink model of the double-shaking table with specimen was established and simulations were conducted in various conditions. The results demonstrate the viability of the proposed DMSTC in that the frequency bandwidth of the double-shaking tables is expanded from 3.27 Hz to 64.57 Hz, the maximum value of differential movement synchronous error is decreased from 1.682 mm to 0.482 mm, and the maximum tracking errors of the two shaking tables decrease from 1.138 mm to 0.044 mm and from 1.030 mm to 0.497 mm, respectively.
National Natural Science Foundation of China5173800151508526514210051. Introduction
Electrohydraulic shaking table is a piece of the most important test equipment for replicating actual vibration situations in seismic tests and many other applications [1–6]. Multisupport, multidimension, and nonuniform excitation experiments have new requirements for shaking table array system in synchronous tracking control [2, 7–9]. Seismic tests of such large-span structures as bridges, dams, railways, pipelines, and so on should consider the traveling wave effect, partial coherent effect, and local site effect to accurately simulate their motion in earthquake. Therefore, synchronous tracking control strategy has become increasingly important and is in need of further development.
The common synchronous control strategies are the parallel control, the master/slave control, the virtual-shaft control, the cross-coupling control, the relative coupling control, adjacent cross-coupling control, ring coupling control, and so on [10–18]. Koren initially proposed the cross-coupling control [12], which, yet, is difficult to be applied to more than two systems. To overcome this limitation, Perez-Pinal et al. proposed the relative coupling control which has stronger synchronous performance and better development [13, 14]. To reduce the complexity of multiple motor synchronization control systems, the adjacent cross-coupling control strategy and ring coupling control strategy were proposed by researchers [15, 16].
To improve synchronous control precision, the control algorithms are also further studied by researchers, such as traditional PID control, adaptive feed forward control, H∞ control, iterative learning control, sliding mode control, fuzzy control, and neural network control [10, 11, 17–25]. Sun and Chiu developed a nonlinear control algorithm to address the motion synchronization of a dual-cylinder electrohydraulic lifting system [23]. Chen et al. proposed an integrated fuzzy controller design that consisted of a pair of cylinder controllers and a motion synchronization controller for synchronous motion and positioning in a hydraulic lifting system [24]. The proposed fuzzy control design effectively handled complex, ill-defined, and uncertain systems. In addition, Liu proposed an intelligent proportional-integral-derivative (PID) algorithm based on back propagation (BP) neural network arithmetic for a three-motor system and obtained satisfactory synchronous decoupling control [25].
A typical synchronization control scheme consists of synchronization control parameters, synchronization control strategy to calculate the parameters errors, and a control algorithm to improve synchronization control precision. C. Chen and L. Chen classified it into three types of control mode: position mode, velocity mode, and torque mode [20]. Tsuruta et al. considered motor torsion, displacement, and speed to improve the stability and synchronous control precision of the studied system [26, 27]. Moreover, Zhang et al. focused on the deviation of the array displacement, speed, and acceleration parameters of shaking tables to elevate the synchronous properties of the system [6]. Fang et al. utilized a synchronous control strategy that combined phase-shifting and amplitude-varying methods to achieve the synchronous motion of two parallel servos [28].
Many studies demonstrate that the feedforward compensator control is an effective method to improve a system’s synchronous tracking accuracy. Zhang et al. applied the adaptive inverse control (AIC) algorithm to two shaking tables’ synchronous tracking controls to improve the tracking precision [6]. In addition, a fuzzy tracking controller (FTC) was designed by Chen et al. to reduce the position-tracking error caused by the parameter uncertainties and nonlinearities of the hydraulic components [24].
However, to date, there are few synchronous tracking control researches in shaking table system with the consideration of the interaction between the specimen and the shaking table system. Particularly, the study of differential movement synchronous tracking control is rare. Therefore, a mathematical model of the double-shaking table system with specimen was built in this research; and the shaking table system of Beijing University of Technology (BJUT) was taken as the research object to verify the synchronous tracking strategy of the DMSTC for shaking table’s motion, thereby establishing a Simulink model of the DMSTC. In this paper, the double-shaking table system was described in Section 1, and the hydraulic system was established in Section 2. Subsequently, the mathematical model of the double-shaking table system with specimen was established and the Simulink model of the DMSTC was presented in Section 3, and the simulation analysis was conducted in Section 4. Finally, in Section 5, a brief conclusion of the paper was given.
2. System Description2.1. Introduction of the Double-Shaking Table System
The system of nine subarray shaking tables, composed of nine single shaking tables of 1 m × 1 m, can be used in shaking table’s array experiments with diversified forms and in various positions. Figure 1 shows the structure of the double-shaking table system. Each single shaking table is composed of platform, supporting guide mechanism, servo valves, hydraulic cylinders, connecting rod, and base. As the core control unit, the servo valves use small control current to control the flow of high pressure oil. The hydraulic cylinders, as the force transform unit, can impose huge force on the platform through the supporting guide mechanism. Each system, with hydraulic cylinders and servo valves, has three exciters, one in X direction and two in Y horizontal direction, and three connecting rods in Z vertical direction.
The structure diagram of the double-shaking table system with specimen.
2.2. The Establishment of the Hydraulic Cylinder Model
The hydraulic system is composed of the hydraulic cylinder, the servo valve, pressure sensor, energy accumulator, and so on. The hydraulic system can be presented in a series of difference equations in the establishment of shaking table’s dynamic model. As the core component of shaking table, the servo valve’s property is critical for the motion of shaking table.
2.2.1. The Servo Valve Model
The servo valve model can be established with the basis of the second-order oscillation link:(1)Gqs=Kvs2/wv2+2ξv/wvs+1,where wv is the inherent frequency of the servo valve, ξv is the damping ratio of the servo valve, and Kv is the flow gain of the servo valve.
2.2.2. The Continuity Equations of Hydraulic Cylinder
Assuming that damping and spring effect are not involved in the hydraulic cylinder system [29], the continuity equations of hydraulic cylinder are(2)QL=kqGqE-KCpL,AppL=Md2xpdt2,QL=Apdxpdt+V4βdpLdt+CCpL,where QL is represented as load flow, kq is flow gain of slid valve near static operation point, Gq is servo value transfer function, E is control error signal, KC is flow pressure parameter of slid valve near static operation point, pL is load pressure, Ap is effective bearing area of valve, M is the gross mass of valve and its load deduction, xp is displacement of valve, V is gross capacity of two hydraulic cylinder chamber, β is oil elastic modulus, and CC is gross leakage coefficient of hydraulic cylinder.
Regarding the platform and specimen as an integrity, and load property as inertial load, the continuity equations of hydraulic cylinder can be obtained after the processing of Laplace transform. They are shown as follows:(3)Ms2x=AppL,QL=Apsx+V4βspL+CCpL,QL=kqGqE-KCpL.
3. The DMSTC Model of the Double-Shaking Table System with Specimen3.1. The Deduction of the Interaction between the Double-Shaking Tables and Specimen
Figure 2 shows the simplified mechanical model of the double-shaking table system with specimen. M1 and M2 are the masses of the two platforms, respectively, and M1=M2=MT; x1 and x2 are their displacements; F1 and F2 are exciting forces of their hydraulic systems; MS is the mass of the specimen; xS is the specimen’s displacement; C is the specimen’s damping coefficient; and K is the specimen’s stiffness coefficient.
Simplified mechanical model.
The force equilibrium equations are(4)MSs2xS+Cs+kxs-x1+Cx+kxs-x2=0,M1s2x1+Cs+kx1-xs=F1,M2s2x2+Cs+kx2-xs=F2.Then(5)MSs2+2Cs+kxs-Cx+kx1-Cs+kx2=0,M1s2+Cs+kx1-Cs+kxs=F1,M2s2+Cs+kx2-Cs+kxs=F2.
The force equilibrium equations are(6)AppL1=F1,AppL2=F2.
Presuming that the parameters of the two exciters are the same, the output flow of salve valves is(7)QL1=kqGqE1-KCpL1,QL2=kqGqE2-KCpL2.
The flow of exciters is(8)QL1=Apsx1+V4βspL1+CCpL1,QL2=Apsx2+V4βspL2+CCpL2.
Equalizing (5) and (6), the result is(9)V4βs+KC+CCpL1=kqGqE1-Apsx1,V4βs+KC+CCpL2=kqGqE2-Apsx2.
We substitute (9) into (4) to get(12)G2M1s2+Cs+kx1-G2Cs+kxs=kqGqE1Ap-sx1,G2M2s2+Cs+kx2-G2Cs+kxs=kqGqE2Ap-sx2.
We simplify (4) to get(13)xs=Cs+kMSs2+2Cs+kx1+x2.
The PID control transfer function is(14)GPID=KP1+1Tis+Tds.
We combine (12) and (13) to get(15)G2M1s2+Cs+k+s-G2Cs+k2MSs2+2Cs+kx1-G2Cs+k2MSs2+2Cs+kx2=GqkqE1Ap,G2M2s2+Cs+k+s-G2Cs+k2MSs2+2Cs+kx2-G2Cs+k2MSs2+2Cs+kx1=GqkqE2Ap.
We take the displacement feedback gain, PID control transfer function, and the displacement meter transfer function into consideration to get(16)E1=GPIDu1-GdKdx1,E2=GPIDu2-GdKdx2.
We write (16) into (15) to get(17)G2M1s2+Cs+k+s-G2Cs+k2MSs2+2Cs+k+GqkqApGdKdx1-G2Cs+k2MSs2+2Cs+kx2=GPIDGqkqu1Ap,G2M2s2+Cs+k+s-G2Cs+k2MSs2+2Cs+k+GqkqApGdKdx2-G2Cs+k2MSs2+2Cs+kx1=GPIDGqkqu2Ap.
Given(18)G3=G2MTs2+Cs+k+s-G2Cs+k2MSs2+2Cs+k+GqkqApGdKd,G4=G2Cs+k2MSs2+2Cs+k,then(19)G3x1-G4x2=GPIDGqkqu1Ap,G3x2-G4x1=GPIDGqkqu2Ap.Finally, the system transfer functions can be obtained:(20)x1=GPIDGqkqApG3u1+G4u2G32-G42,x2=GPIDGqkqApG3u2+G4u1G32-G42.
3.2. The Proposition of the DMSC
Figure 3 shows the model of the double-shaking table system in different coordinates of the double-shaking table system, in which the input coordinate system and the output coordinate system were established. In the input coordinate system, u1 and u2 are input signals and ϕ is the desired rotation angle. In the output coordinate system, x1 and x2 are output signals, l is the distance between the shaking table and the central line of the double-shaking table system, and φ is the actual rotation angle of the double-shaking table system.
The model of the double-shaking table system in different coordinates.
We substitute (13) into (12) to get(21)G2M1s2+Cs+k+sx1-G2Cs+k2MSs2+2Cs+kx1+x2=GqkqE1Ap,G2M2s2+Cs+k+sx2-G2Cs+k2MSs2+2Cs+kx1+x2=GqkqE2Ap.
Given(22)G5=G2MTs2+Cs+k+s,G4=G2Cs+k2MSs2+2Cs+k
and writing (22) into (21) to get(23)G5x1-G4x1+x2=GqkqE1Ap,G5x2-G4x1+x2=GqkqE2Ap,then(24)G5-2G4x1+x2=GqkqApE1+E2,G5x1-x2=GqkqApE1-E2.
We deduct from the output signal coordinate system of Figure 3 with geometric method to acquire(25)x1=x+lφ,x2=x-lφ.Then(26)x1+x2=2x,x1-x2=2lφ.
We deduct from the input signal coordinate system in Figure 3 to acquire(27)u1=u+lϕ,u2=u-lϕ.Then(28)u1+u2=2u,u1-u2=2lϕ.
In order to get the control error signal of the double-shaking table system, the controlling block diagram is designed as shown in Figure 4, in which ed is the translation error and ea is the rotation error.
Controlling block diagram.
We deduct from Figure 4 to get(29)E1=ed+ea,E2=ed-ea.Then(30)E1+E2=2ed,E1-E2=2ea.
We take the displacement feedback gain, PID control transfer function, and the displacement meter transfer function into consideration to get(31)ed=12GPIDu1+u2-GdKdx1+x2,ea=12GPIDu1-u2-GdKdx1-x2.
We substitute (26) and (28) into (31) to get(32)ed=GPIDu-GdKdx,ea=GPIDlϕ-GdKdlφ.
We write (32) into (30) to get(33)E1+E2=2GPIDu-GdKdx,E1-E2=2GPIDlϕ-GdKdlφ.
We write (26) and (33) into (24) to get(34)G5-2G4+GqkqApGdKdx=GPIDGqkqApu,G5+GqkqApGdKdφ=GPIDGqkqApϕ.
Equation (34) can be transformed into(35)xu=GPIDGqkqAp1G5-2G4+Gqkq/ApGdKd,φϕ=GPIDGqkqAp1G5+Gqkq/ApGdKd.
We write (28) and (33) into (24) to get(36)G5-2G4+GqkqApGaKdx1+x2=GPIDGqkqApu1+u2,G5+GqkqApGaKdx1-x2=GPIDGqkqApu1-u2.
Finally, the system’s transfer function can be acquired:(39)x1=GPIDGqkq2Apu11G6+1G7+u21G6-1G7,x2=GPIDGqkq2Apu11G6-1G7+u21G6+1G7.
We write (20) and (39) into the form of transfer function matrix as follows:(40)x1x2=H11H12H21H22u1u2.Then the system’s transfer functions of double-shaking table system with specimen are(41)H11=H11=H22=GPIDGqkqApG3G32-G42,H12=H12=H21=GPIDGqkqApG4G32-G42.Then we apply the DMSC strategy to the double-shaking table system with specimen to get the transfer function: (42)H11=G11=G22=GPIDGqkq2Ap1G6+1G7,H12=G12=G21=GPIDGqkq2Ap1G6-1G7.
Equation (35) shows that the displacement and the rotation are independent freedom degrees without coupling. Equation (41) shows the system’s transfer functions of double-shaking table system with specimen. Equation (42) shows the system’s transfer functions of double-shaking table system with specimen utilizing DMSC strategy.
The above-mentioned differential movement synchronous control is divided into two parts: error processing and error coupling. The Simulink model is shown in Figure 5.
Simulink model of the DMSC.
Error processing in detail
Error coupling in detail
3.3. The Design of the AIC
In Section 3.2, the DMSC is developed to improve the synchronization precision of the two shaking tables. In this section, this study combined the AIC based on Normalized Least Mean Square (NLMS) algorithm with the proposed synchronous strategy to realize synchronous tracking control of the double-shaking table system.
As shown in Figure 6, the error signal calculated by the desired signal and the exciter output signal is utilized to adjust the weight coefficient of the controller, thereby realizing the waveform reproduction of the shaking table.
The adaptive inverse control system structure of the shaking table.
The weight updating formula of the NLMS algorithm is shown as follows:(43)wn+1=wn+2μenxn,μ=αβ+xn2,where w(n) is the weight coefficient of the N-rank finite impulse response (FIR) digital filter, μ is the convergence factor constant, e(n) is the error of the filter, x(n) is the input signal of the filter, and x(n)2 is the input signal’s Euclidean norm square.
The designed adaptive inverse controller is shown in Figure 7. In this Simulink model, α=0.25 and β=0.5.
The adaptive inverse controller.
3.4. The Model of the DMSTC3.4.1. The Model of the Double-Shaking Table System with Specimen
The model of the double-shaking table system with specimen is shown in Figure 8. Table 1 shows the parameters of the model of the double-shaking table system and the PID control parameters are shown in Table 2. According to the hydraulic system model established in Section 2.1 and the deduction of the interaction between the shaking table and the specimen in Section 3.1, the model is established in Matlab/Simulink with the parameters in Table 1. In the simulations, the mass of the specimen is 2000 kg, the frequency is 3 Hz, and the damping ratio is 0.02.
The parameters of the double-shaking table model.
Parameters
Values
Mass of the platform
1t
Bulk modulus of oil
7×108N/m2
Flow coefficient of servo valve pressure
1.0×10-11m5/N/s
Spool flow gain
6×10-3m3/s⋅V
Inherent frequency of servo valve
628rad/s
Leakage coefficient of hydraulic cylinder
7.5529×10-15m5/N/s
Damping ratio of servo valve
0.7
Effective bearing area of valve
4.73×10-3m2
Displacement feedback gain
0.21V/m
Equivalent volume of cylinder
3.9023×10-3m3
Damping ratio of sensor
0.7
Inherent frequency of sensor
942rad/s
The PID control parameters of the double-shaking table system.
PID parameters
Kp
Ki
Kd
Shaking table
0.6
0.2
0.025
The model of double-shaking table system with specimen.
3.4.2. The Model of the DMSTC
On the ground of the above deduction, analysis, and modeling, a synchronous tracking control system (including the double-shaking table system and specimen, the DMSC block, and the AIC blocks) is established in Matlab/Simulink, as shown in Figure 9.
The DMSTC model.
4. Results and Discussions4.1. The System’s Frequency Response Performance
The frequency transform functions (FRFs) of the double-shaking table system without DMSTC and with DMSTC are shown in Figure 10. The frequency bandwidth of the system without DMSTC-H11 is only 3.27 Hz, which leads to different dynamic responses and synchronization error, while the frequency performance of the system with DMSTC-G11 is increased from 3.27 Hz to 64.57 Hz, which improves the accuracy of the signal representation and reduces the synchronization error simultaneously.
The FRFs of the double-shaking table system.
4.2. Random Wave Test
The random waves 1 and 2, which are input into shaking tables 1 and 2, respectively, are shown in Table 3. The results without DMSTC and with DMSTC are shown in Figure 11. For comparative study, the details of Figure 11(b) are enlarged.
Different random waves.
Wave parameters
Amplitude
Frequency
Random wave 1
±2.0 mm
0.4 Hz–5 Hz
Random wave 2
±1.0 mm
20 Hz–40 Hz
Random wave test.
Random wave 1
Random wave 2
Table 4 shows that, with DMSTC, the replication accuracy of each shaking table is improved. For shaking table 1, the correlative coefficient rises from 72.7% to 99.8%. For shaking table 2, the correlative coefficient rises from 58.6% to 82.7%.
Random waves’ correlative coefficient of the double-shaking table system.
Shaking table
Without DMSTC
With DMSTC
1
72.7%
99.8%
2
58.6%
82.7%
Figures 12 and 14 show that, with DMSTC, the tracking accuracy of each shaking table is improved. For shaking table 1, the maximum tracking error, max(u1-x1), decreases from 1.138 mm to 0.044 mm and the root mean square tracking error, RMS(u1-x1), decreases from 0.454 mm to 0.017 mm. For shaking table 2, the maximum tracking error, max(u2-x2), decreases from 1.030 mm to 0.497 mm and the root mean square tracking error, RMS(u2-x2), decreases from 0.296 mm to 0.144 mm.
The tracking error of the double-shaking table.
Tracking error of shaking table 1
Tracking error of shaking table 2
Figures 13 and 14 demonstrate that the differential movement synchronous performance is greatly improved. Namely, the error of the desired deviation and the accurate deviation decreases. The maximum value of the differential movement synchronous error, max(E), decreases from 1.682 mm to 0.482 mm, and the root mean square value of the differential movement synchronous error, RMS(E), decreases from 0.541 mm to 0.145 mm.
The differential movement synchronous error.
Data analysis of the tracking errors and synchronous error.
5. Conclusion
In conclusion, this research proposed a novel synchronous tracking strategy, DMSTC strategy, for shaking table’s motion to generate the control error signals while taking the interaction between shaking tables and specimen into consideration and to improve the system’s tracking precision. DMSTC Simulink model of the double-shaking table with specimen was established and simulations were conducted in various conditions. The results of the study are as follows:
(1) The proposed DMSTC strategy extends the double-shaking table system’s frequency bandwidth from 3.27 Hz to 64.57 Hz, which improves the frequency bandwidth of the double-shaking table system.
(2) The replication accuracy of each shaking table is improved. For shaking table 1, the correlative coefficient rises by 27.1%. For shaking table 2, the correlative coefficient rises from 58.6% to 82.7%.
(3) The tracking accuracy of each shaking table is improved. For shaking table 1, the maximum tracking error decreases from 1.138 mm to 0.044 mm. For shaking table 2, the maximum tracking error decreases by 51.74%.
(4) The differential movement synchronous error decreases, the maximum absolute value of that decreases by 71.34%, and the root mean square value of that decreases by 73.20%.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China under Grant nos. 51738001, 51508526, and 51421005.
ShenG.ZhuZ.-C.LiX.TangY.HouD.-D.TengW.-X.Real-time electro-hydraulic hybrid system for structural testing subjected to vibration and force loading20163349702-s2.0-8495950903610.1016/j.mechatronics.2015.10.009GaoC. H.JiJ. B.YanW. M.LiN.Developments of shaking table technology in China2014478919ShenG.LiX.ZhuZ.TangY.ZhuW.LiuS.Acceleration tracking control combining adaptive control and off-line compensators for six-degree-of-freedom electro-hydraulic shaking tables20177032233710.1016/j.isatra.2017.07.0182-s2.0-85025650844ShenG.ZhengS.-T.YeZ.-M.HuangQ.-T.CongD.-C.HanJ.-W.Adaptive inverse control of time waveform replication for electrohydraulic shaking table20111711161116332-s2.0-8005473357910.1177/1077546310380431ShenG.ZhuZ. C.ZhangL.YuT.Chi-fuY.Jin-songZ.Guang-daL.Jun-WeiH.Adaptive feed-forward compensation for hybrid control with acceleration time waveform replication on electro-hydraulic shaking table20132181128114210.1016/j.conengprac.2013.03.0072-s2.0-84882714237ZhangL. P.CongD. C.YangZ. D.Attitude synchronous tracking control of double shaking tables based on hybrid fuzzy logic cross-coupled controller and adaptive inverse controller20152925372546JiJ. B.LiF. F.LiZ. B.Research and advances on the control technology of the multiple shaking table array system201228696101JiJ. B.LiX. L.YanW. M.Research on shaking table array and dynamic model test201127s13136ZongZ. H.ChenL.HuangY.Research and application of testing technology of the earthquake simulation multiple shaking table facilities201127s1614LiL.SunL.ZhangS.Mean deviation coupling synchronous control for multiple motors via second-order adaptive sliding mode control20166222223510.1016/j.isatra.2016.01.0152-s2.0-84968761527LiL.-B.SunL.-L.ZhangS.-Z.YangQ.-Q.Speed tracking and synchronization of multiple motors using ring coupling control and adaptive sliding mode control20155863564910.1016/j.isatra.2015.07.0102-s2.0-84943583292KorenY.Cross-coupled biaxial computer control for manufacturing systems1980102426527210.1115/1.3149612Zbl0471.930312-s2.0-0019247725Perez-PinalF. J.CalderonG.Araujo-VargasI.Relative coupling strategyProceedings of the IEEE International Electric Machines and Drives Conference, IEMDC 2003June 2003Madison, WI, USA1162116610.1109/IEMDC.2003.12103872-s2.0-84946026142Perez-PinalF.NunezC.AlvarezR.CervantesI.Comparison of multi-motor synchronization techniquesProceedings of the 30th Annual Conference of IEEE Industrial Electronics Society, 2004. IECON 2004Busan, South Korea1670167510.1109/IECON.2004.1431832ZhaoD.LiC.RenJ.Speed synchronization of multiple induction motors with adjacent cross coupling controlProceedings of the 48th IEEE Conference on Decision and Control held jointly with 2009 28th Chinese Control Conference, CDC/CCC 2009December 2009chn680568102-s2.0-7795083313910.1109/CDC.2009.5400835SunJ.LiuR.LuoY.SunW.Research on Multi-motor Synchronization Control for Cutter Head of Shield Machine Based on the Ring Coupled Control Strategy20095928Berlin, HeidelbergSpringer Berlin Heidelberg345354Lecture Notes in Computer Science10.1007/978-3-642-10817-4_34LaiX.XieW.Theoretical and Experimental Study on Electromechanical Coupling Properties of Multihammer Synchronous Vibration System2016201611670726410.1155/2016/67072642-s2.0-84962736225ZhangX.-L.ZhaoC.-Y.WenB.-C.Theoretical and experimental study on synchronization of the two homodromy exciters in a non-resonant vibrating system20132023273402-s2.0-8487639065910.3233/SAV-2012-00747OuyangP. R.DamT.PanoV.Cross-coupled PID control in position domain for contour tracking2015336135113742-s2.0-8493038152910.1017/S0263574714000769ChenC.ChenL.Robust cross-coupling synchronous control by shaping position commands in multiaxes system20125912476147732-s2.0-8486376077110.1109/TIE.2011.2182012BartonK. L.AlleyneA. G.A cross-coupled iterative learning control design for precision motion control2008166121812312-s2.0-5534908506210.1109/TCST.2008.919433ZhaoD.LiS.GaoF.ZhuQ.Robust adaptive terminal sliding mode-based synchronised position control for multiple motion axes systems20093113615010.1049/iet-cta:200702722-s2.0-58149522273SunH.ChiuG. T.-C.Motion synchronization for dual-cylinder electrohydraulic lift systems20027217118110.1109/TMECH.2002.10112542-s2.0-0036613152ChenC.-Y.LiuL.-Q.ChengC.-C.ChiuG. T.-C.Fuzzy controller design for synchronous motion in a dual-cylinder electro-hydraulic system20081666586732-s2.0-4084913935210.1016/j.conengprac.2007.08.005LiuX.ZhangX.HuJ.LiuY.LiuG.ZhaoL.Three-motor synchronous decoupling control based on BP neural networkProceedings of the 2009 Chinese Control and Decision Conference, CCDC 2009June 2009Guilin, China5365536810.1109/CCDC.2009.51950722-s2.0-70449449679TsurutaK.OjiroT.HondaH.Synchronization control for twin-axes table drive systemProceedings of the 41st Annual Conference of the IEEE Industrial Electronics Society, IECON 2015November 2015Yokohama, Japan5162516710.1109/IECON.2015.73929102-s2.0-84973138036TsurutaK.OjiroT.UshimiN.ShibataH.Vibration suppression control for multiaxis table drive systemProceedings of the 2013 10th IEEE International Conference on Mechatronics and Automation, IEEE ICMA 2013August 2013Takamatsu, Japan1135114010.1109/ICMA.2013.66180742-s2.0-84887971832FangJ. H.WeiJ. H.KongX. W.Synchronous control strategy for paralleled servo valves201246610541059HuangH. H.2008Beijing, ChinaSeismological Press