Published: 30 September 2016

Vector matching-based disturbance rejection method for load simulator

Xuesong Yang1
Changchun Li2
Hao Yan3
Jing Huang4
1, 2, 3, 4Beijing Jiaotong University, Beijing, China
Corresponding Author:
Xuesong Yang
Views 36
Reads 16
Downloads 1542


Position disturbance is the most troublesome problem for load simulator. It is difficult to work properly for a loading system, when the load-carrying object is moving. Therefore, it is of high significance to compensate for the position disturbance. To achieve this objective, a kind of vector matching-based disturbance rejection (VMBDR) method is proposed. If the input is periodic and the load simulator system satisfies linear superposition, a torque command can be constructed to eliminate the position disturbance of the load-carrying object, according to VMBDR. Fortunately, linear superposition of electric load simulator is verified by experiments, and most of the position disturbance is cyclical. Compared with other methods, the VMBDR does not depend on building a precise mathematical model. Therefore, it has better maneuverability and adaptability. Finally, experimental results show that VMBDR has effective disturbance rejection performance.

1. Introduction

Load simulator is a kind of hardware-in-the-loop simulation equipment [1]. For the most part, they have three components: a loading system, a load-carrying object, a measuring-control system. The purpose of load simulation is to test the performance of a load-carrying object under the condition as real as possible. It helps designers to foresee and detect problems related to mechanics and control of the load-carrying object before it even is created, hence saving a lot of time and cost. There are many successful examples of the application of load simulators relate to aerospace [2, 3], robotics, mechatronics [4], automation [5] and other fields.

In the process of load simulation, a load-carrying object will be moving according to the test signal. For example, in the field of aerospace, the rudder actuator is often implemented as a load-carrying object. It must move, just like when the rocket is in the air. This motion is called as position disturbance. In this condition, the load simulator is not able to accurately control the loading force unless the control system takes the motion into account and makes an appropriate compensation. There are many kinds of position disturbance rejection methods. These methods can be divided into four categories, i.e., the position synchronization information feedforward method, the disturbance force/torque observation method, the composite control method and the disturbance model matching method.

Synchronization is the key to the first method. That is to say, if a loading system can be kept relatively static with a load-carrying object, disturbance torque could be eliminated completely. Under the guidance of this idea, researchers conducted a variety of attempts, including installing the loading motor on a position synchronous motor, as known as the hardware synchronous method [6]. To a certain extent, this method can eliminate the surplus torque caused by speed interference. It is suitable for low frequency and small amplitude disturbance case. But it also makes the system more complex, causes higher cost. In this respect, software synchronization method has obvious advantages. Among the software methods, structure invariance principle [7] is the most widely developed and has the most profound influence [8]. This principle is derived from invariance principle which was proposed by Soviet scholar Г. В. Щипанов in 1939 [9]. The author pointed out that, in a control system with external disturbances, the controller should be divided into two parts. One part is used to compensate the disturbance. The other is used to realize tracking. But all the rejection methods based on the structural invariance principle require a precise mathematical model of the system, which is hard to achieve.

The focus of the second method is concerned with disturbance torque itself. Fang Qiang proposed a disturbance observer design method for electric load simulator [10]. This observer can estimate generalized disturbance information which is a decoupling disturbance signal, through the disturbance torque and the loading command signal. This design makes great progress, considering uncertainties of the mathematical model. However, compared with the real disturbance, a time delay of the estimated disturbance is inevitable, making it difficult to completely eliminate the disturbance.

The characteristic of composite control strategy is that it is necessary to handle position synchronization information and feedback torque at the same time. The composite controller includes two parts. One part is an intelligent controller to eliminate the influence of position disturbance, and the other is a PID controller to ensure stability of the system during learning process of the intelligent controller. Neural networks [11], fuzzy logic [12] and iterative learning theory [13] are employed in the design of intelligent controllers. Although this method can deal with nonlinearity and uncertainty of the system, the convergence is still needed to be careful in.

The matching relation between position disturbance model and loading system model is what the last method looks for. A dynamic fuzzy neural network is applied to design the controller of load simulator [14, 15]. The idea of this method is to construct an inverse disturbance model, based on which a torque command to counteract the effect of position disturbance is generated. This is similar to the deconvolution method. However, it is hard to achieve due to the difficulty of the shock response acquisition, model uncertainty, nonlinearity and time variation. Robust methods and game theory have been used to solve the above problems [16, 17], and make deconvolution problem solved partly.

In summary, each of these methods has both advantages and disadvantages. Classical control theory is simple and powerful in form, but lack of adaptability. While modern control theory and intelligent control theory are enough to deal with complex problems, it is lack of stability and maneuverability. The key to success is to maximize the advantages of the method, according to the actual situation. For a load simulator, sine and square waves are the most widely used test commands of the load-carrying object. The common point is that these signals are cyclical. In this way, the test commands are able to be expanded into Fourier series, and each frequency component can be regarded as a signal vector.

Motivated by the above discussion, a new idea is presented in this paper, which is based on the matching relationship between load-carrying object position signal and torque signal of the loading system. Then, the matching problem of the disturbance model is simplified to a single or finite frequency sequence matching problem, called as vector matching problem. The critical part of the method is to find the relationship between the phase and amplitude of different signal vectors. This method is much easier to operate and needs no accurate mathematical model.

The rest of this paper is organized as follows. In Section 2, the structure of a load simulator and the principle of vector matching are introduced. The algorithm of vector matching-based disturbance rejection (VMBDR) method is presented in Section 3. In Section 4, some experiments are designed to verify the algorithm. Finally, results and further studies are discussed in the conclusion.

2. System description and vector matching principle

The electric load simulator (abbr. ELS) consists of three parts. 1) A servo mechanism under test, which is the actuator of a rudder in this experimental project, is used as the load-carrying object. It is connected to a torque sensor through a coupling. During the experiment, the servo mechanism rotates following the position commands. 2) A loading system, consisting of a PMSM and a motor driver, is connected to the torque sensor as well. 3) A measuring-control system mainly includes a torque sensor and an industrial computer with data acquisition cards, collecting the actual load torque in real time. The outputs of measuring-control system are angular displacement commands and loading system torque commands, respectively.

ELS is a typical double input and single output linear system. For a linear system, as showed in Fig. 1, use yellow on behalf of the position command signal and blue on behalf of the torque command signal, then the feedback loading torque signal should look like green. This is the superposition of linear systems, which is the basis of vector matching.

Fig. 1Construction diagram of an ELS

Construction diagram of an ELS

In order to explain the principle of vector matching, assume that the position command and the torque command are sinusoidal signals of the same frequency. These two commands generate two components of the resultant torque respectively. Fig. 2 shows the relationship between these three torque signals in polar coordinate system.

Fig. 2Schematic diagram of vector matching relation

Schematic diagram of vector matching relation

The symbols in Fig. 2 are defined as follows:

T0 – resultant torque vector measured by the torque sensor,

T1 – torque vector produced by the movement of the load-carrying object,

T2 – torque vector produced by the loading system under the torque command Tc,

T20 – one component of T2,

T21 – another component of T2,

β – phase angle of T0,

φ – phase angle of T1,

θ – angle between T1 and T2.

According to the invariance principle, Tc should consist of two parts, i.e. Tc=Tc0+Tc1. Tc0 is used to compensate for the position disturbance, and Tc1 is used to realize the loading. The purpose of vector matching is to make T20 equal to -T1.When the load target is 0 Nm, T21 should be equal to 0, and T2 should be equal to -T1.

Therefore, the vector relation described in Fig. 2 is the basic of the vector matching method, from which we can see that obtaining the angle θ is the key of vector matching. Fig. 3 demonstrates the working principle of VMBDR for ELS. Some symbols are defined in the figure and the others are defined as follows:

Lm – Motor stator cyclic inductance,

Rm – Motor stator resistance,

Jm – Moment of inertia on the motor shaft,

Bm – Viscous friction coefficient on motor shaft,

Ke – Back EMF coefficient,

Kt – Torque coefficient,

U – Armature voltage,

E – Back- EMF,

I* – Current command,

i – Armature current,

TL – Output load torque,

θout – Elastic deformation of the shaft.

Fig. 3VMBDR of the electric load simulator

VMBDR of the electric load simulator

As shown in Fig. 3, the vector matching process is divided into three basic steps. The first two steps are the initialization of VMBDR. Its purpose is to gather the necessary information for step 3, i.e., the amplitude and phase of T1 and T2.

In the control system, there is a switch, which makes the control work in accordance with steps of 1, 2 and 3. In order to facilitate description, set the loading command Tc1 to 0 Nm at first. Then, when the switch is at position 1, the compensating signal Tc0 is zero, and the original position disturbance torque vector T1 could be obtained. As the switch is at position 2, the compensating signal Tc0 is a filtered position command θc, and the resultant torque vector T0 could be gotten. Usually, the filter is a proportion link. So far, vector matching unit uses T0 and T1 to compute the loading torque vector T2 and θ, then calculates a new compensating signal, which can eliminate T1 completely in step 3.

Changes in the environment will affect many physical parameters of ELS – motor stator resistance Rm, for instance. This effect cannot be ignored when the test time is long enough, although this process is usually slow. Therefore, the control error of T20-T1 should be monitored in real time, and T20 should be fine-tuned by vector matching unit, when the error is accumulated to a certain extent.

Noise is another factor which has a great influence on the accuracy of vector matching, in addition to the parameter variability. The three inputs of the vector matching unit are torque command Tc, position command θc and output load torque TL measured by the torque sensor, as shown in Fig. 3. Tc and θc are generated by ELS controller computer itself. So only TL could be contaminated by the additive noise. A Butterworth filter could be used to process the contaminated TL, when the noise is too large. However, the influence of the filter on the position disturbance torque vector T1 and the loading torque vector T2 are the same. Therefore, the matching relation between T1 and T2 does not change.

3. Vector-matching algorithm

ELS is a typical double input linear system. For a linear system, the output torque should obey the principle of superposition. So transfer function of the TL is:


In the following, WMs and WRs are called non-disturbing load transfer function and position disturbance transfer function, respectively. Their frequency domain expressions are described as:


Assuming there is a sinusoidal position disturbance sinωt. When the loading command Tc1 is equal to 0 Nm, the output torque is only T1, and T1=ARωsinωt+φRω. The total output TL is expected to be 0 Nm. Thus, according to the principle of vector matching and the frequency characteristics of the model, a torque command could be constructed as follows:


So the core of VMBDR is to identify the amplitude and phase of each vector signal. There are two kinds of VMBDR algorithms, corresponding to the sine position disturbance and the other periodic position disturbance respectively.

3.1. Single vector matching algorithm

Single vector matching algorithm (abbr. SVM) is designed to compensate for the position disturbance of sine type. This algorithm is an on-line identification algorithm based on recursive least square method for the amplitude and phase of sinusoid. Moreover, the least square method also has excellent anti-noise ability.

Assuming that there is a sinusoidal signal vt, as Eq. (4). The sampling frequency is fs. A set of M sampling points are obtained.


Eq. (4) can be expanded to:


In Eq. (5), C0=Acosϕ, C1=Asinϕ, C2=W. Let D0t=sinωt, D1t=cosωt, vt=TLt, θ^Tk=C0k C1k C2k, hTk=D0k D1k 1, k=1,2,,M.

The parameters to be identified are C0, C1and C2.The measured sample sequence is v1, v2, …, vM. Then, the residual error δ of v^ is:


In Eq. (6). Sum of squares of residuals is:


According to recursive least square method [18], the recurrence formula for the amplitude and phase is:


In Eq. (8), K and P are respectively the recursive weighting factor vector and K update vector. I is a unit matrix. Finally, we can get the optimal value of C0M, C1M and C2M. The amplitude and phase angle could be calculated as follows:


Take the time domain data of the torque vector T0 and T1 into the recursive formula. The amplitudes A0, A1 and phases φ, β could be calculated. The amplitude A2 and angle θ could be gotten by cosine theorem, as Eq. (10) and Eq. (11). The direction of angle θ depends on β-φ, the value range of φ and β is π, π.


So, Tc0 needed by step 3 of vector matching method is described as follows:


In Eq. (12), Ac0' and ϕc0' are amplitude and phase of Tc0', which is the compensating signal applied by the system controller in step 2.

3.2. Complex vector matching algorithm

While the position disturbance is a non-sinusoidal periodic signal, complex vector matching algorithm (abbr. CVM) could be used. Although the Fourier series of the periodic signal includes infinite harmonic components, in practice, the load simulator is equivalent to a low pass filter. Thus the disturbance signal only retains a finite number of harmonic components. Batch processing of a finite number of vectors is realizable. This algorithm is an off-line identification algorithm based on FFT and IFFT.

The basic idea of CVM is consistent with SVM. The difference is that CVM works in the complex domain. The process of CVM is illustrated in Fig. 4.

Fig. 4Principle diagram of the complex vector matching method

Principle diagram of the complex vector matching method

Fig. 4 introduces a number of vector symbols as follows (n is the length of FFT). Each of them represents a Fourier transform of a signal. Each of the elements in the vectors is written in complex index form.

θc1=θc11, θc12, ,θc1nT – Position command complex vector in step 1,

θc2=θc21, θc22, ,θc2nT – Position command complex vector in step 2,

θc3=θc31, θc32, ,θc3nT – Position command complex vector in step 3,

Tc2=Tc21, Tc22, ,Tc2nT – Compensating command complex vector in step 2,

Tc3=Tc31, Tc32, ,Tc3nT – Compensating command complex vector in step 3,

TL1=TL11, TL12, ,TL1nT – Output load torque complex vector in step 1,

TL2=TL21, TL22, ,TL2nT – Output load torque complex vector in step 2,

H1=H11, H12, ,H1nT – Relative load torque complex vector in step 1,

H2=H21, H22, ,H2nT – Relative load torque complex vector in step 2,

E=E1, E2, ,EnT – Pure load torque complex vector in step 2,

A=A1, A2, ,AnT – Gain of complex vector matching used in step 3.

In step 1 the compensation command is zero, so TL1 contains the position disturbance torque only. However, in step 2, the compensation command is Tc2, then TL2 turns out to be:


The object of complex vector matching is to construct a torque command vector, so that the error torque is close to zero in step 3. Thus, Eq. (14) should be satisfied.


The key to achieve the complex vector matching is to obtain WRjkω0 and WMjkω0, which is essentially a model matching problem. Reference Eq. (13) and Eq. (14), the vector matching relationship can be described as follows:


Eq. (15) shows the calculation process of the matching compensation torque command Tc3. Normally, rudder position command θc3 is set in advance, according to experimental requirements. Once the complex gain vector A is determined, the compensation torque command Tc3 can be matched successfully.

The Fourier series of a periodic signal is a discrete sequence, and the non-zero components only exist on the harmonic frequencies of the periodic signal. Hence the actual calculation procedure of CVM is as follows:

Ak=-H1kθc2kEkTc2k,Ek0 and Tc2k0,0,Ek=0 or Tc2k=0,
i=1,2, k=1,2,,n.

As shown in Eq. (16) to Eq. (18), only the components of harmonic frequencies are involved in the CVM operation. On the other hand, the noise components are distributed on the spectrum weakly and evenly, when the ELS has some measures to resist interference. Therefore, the effect of noise on the CVM operation is limited.

4. Experiment

Experiments are carried out on a servo motor limit capacity test bench, shown in Fig. 5. Similar to the structure described in Fig. 1, the test bench mainly includes a loading motor, a motor under test which is the load-carrying object, and a torque sensor. Meanwhile, the torque sensor is grounded and isolated to reduce the influence of electromagnetic interference on torque measurement.

Fig. 5Photos of servo motor limit capacity test bench

Photos of servo motor limit capacity test bench

The experiments including three sub-experiments could be divided into two categories. One is constant torque loading test, which means that the loading motor should load a constant torque on the motor under test. The loading torque is independent from the motion of the motor under test. The other is called as spring torque loading test. In this experiment, the loading motor is supposed to provide a torque proportional to the rotation angle of the motor under test. So the loading motor operates like a spring for the tested motor. Table 1 shows the arrangement of these experiments.

Table 1Arrangement of these experiments

Test No.
Waveform of the tested motor position command
Test category
Constant torque loading
Constant torque loading
Spring torque loading

4.1. Sinusoid position – constant torque loading

In this experiment, position command of the tested motor is a series of sinusoidal signals with the frequencies of 5 Hz, 10 Hz, 15 Hz. Control object of load torque is 0 Nm. Fig. 6, Fig. 7, and Fig. 8 demonstrate that, by using SVM, the residual torque error is obviously smaller than the original torque error. Detail data are provided in Table 2.

Fig. 65 Hz sinusoid position disturbance constant torque load curve

5 Hz sinusoid position disturbance constant torque load curve

Fig. 710 Hz sinusoid position disturbance constant torque load curve

10 Hz sinusoid position disturbance constant torque load curve

Fig. 815 Hz sinusoid position disturbance constant torque load curve

15 Hz sinusoid position disturbance constant torque load curve

Table 2The detail data of the 1st experiment

Frequency (Hz)
Original torque standard deviation (Nm)
Residual torque standard deviation (Nm)
Elimination %

4.2. Triangle position – constant torque loading

A 5 Hz triangle signal xt, whose magnitude is represented by Atri, is adopted as the position command for the tested motor in this experiment. The Fourier series of xt is shown in Eq. (19). CVM is utilized to improve the loading precision. By comparing Fig. 9 and Fig. 10, we can find that CVM greatly reduces the torque control error.


Fig. 9Torque error caused by position disturbance of triangular wave

Torque error caused by position disturbance of triangular wave

Fig. 10CVM effect of triangle wave position disturbance

CVM effect of triangle wave position disturbance

Fourier transform of the test data illustrates the effect of CVM more intuitively. Fig. 11 shows amplitude spectrum of main signals, including the position command θc3, the compensation torque command Tc3, the original torque error and the residual torque error. The effect of CVM is expressed numerically in Table 3.

Fig. 11Spectrum of triangular disturbance rejection experimental data

Spectrum of triangular disturbance rejection experimental data

a) Position command

Spectrum of triangular disturbance rejection experimental data

b) Original torque error

Spectrum of triangular disturbance rejection experimental data

c) Matching torque command

Spectrum of triangular disturbance rejection experimental data

d) Residual torque error

Table 3Rejection effect of disturbance torque at each frequency component

Frequency (Hz)
Original torque standard deviation (Nm)
Residual torque standard deviation (Nm)
Elimination %

4.3. Sinusoid position – spring torque loading

In this experiment, position command of the tested motor is a 5 Hz sinusoidal signals. And loading motor must provide a load torque, with a load gradient of 0.05 Nm/deg. Fig. 12 shows that when only the feedback controller is working, load torque is obviously lagging behind the rotation of the tested motor. As demonstrated in Fig. 13, load torque and motor rotation are basically synchronous after applying the vector matching algorithm. The original torque error displays in Fig. 14 including two parts. One results from the lag of loading system. The other is caused by the position disturbance coming from the tested motor. Fig. 14 shows that VMBDR is highly effective at suppressing the two kinds of loading errors. The standard deviation of the load torque before and after compensation is: 1.4536 Nm and 0.25432 Nm.

Fig. 125 Hz spring torque loading before vector matching

5 Hz spring torque loading before vector matching

Fig. 135 Hz spring torque loading after vector matching

5 Hz spring torque loading after vector matching

Fig. 14Error contrast before and after vector matching

Error contrast before and after vector matching

5. Conclusion

Load simulator is a passive torque servo system and it is necessary to suppress the position disturbance. In this study, vector matching-based disturbance rejection (VMBDR) method is proposed. The effectiveness of the proposed algorithm is verified on an electric load simulator.

Single vector matching (SVM) method can eliminate sinusoidal position disturbance, and its elimination ratio is related to the precision of identification. Complex vector matching (CVM) method based on FFT can eliminate cyclical position disturbance. The experimental results show that elimination ratio of CVM is higher than SVM. In addition, the vector matching method can compensate the control error caused by the response lag of the loading system as well. This means that VMBDR is not only applicable to constant torque loading, but also to spring torque loading.

Finally, VMBDR works effectively on the electric load simulator which has high linearity. This suggests that, for other linear multi-input systems, vector matching method will be still effective in great probability.


  • Temeltas H., Gokasan M., Bogosyan O. S. A nonlinear load simulator for robot manipulators. The 27th Annual Conference of Industrial Electronics Society, IEEE, Vol. 1, 2001, p. 357-362.
  • Ullah N., Shaoping W. High performance direct torque control of electrical aerodynamics load simulator using adaptive fuzzy backstepping control. Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, Vol. 229, Issue 2, 2014, p. 369-383.
  • Nam Y., Hong S. K. Force control system design for aerodynamic load simulator. Control Engineering Practice, Vol. 10, Issue 5, 2002, p. 549-558.
  • Ryu H. M., Kim S. J., Sul S. K., et al. Dynamic load simulator for high-speed elevator system. Proceedings of the Power Conversion Conference, IEEE, Vol. 2, 2002, p. 885-889.
  • Wang X., Wang S., Yao B. Adaptive robust torque control of electric load simulator with strong position coupling disturbance. International Journal of Control, Automation and Systems, Vol. 11, Issue 2, 2013, p. 325-332.
  • Meng Q., Zhang L., Liu Q., et al. Experimental and theoretic study on eliminating the disturbance torque and widening the frequency bandwidth of the load simulator using position synchro compensation. Journal of Astronautics, Vol. 18, Issue 1, 1997, p. 120-123.
  • Li Y. Development of hybrid control of electrohydraulic torque load simulator. Journal of Dynamic Systems, Measurement, and Control, Vol. 124, Issue 3, 2002, p. 415-419.
  • Wang C., Jiao Z., Wu S., et al. An experimental study of the dual-loop control of electro-hydraulic load simulator (EHLS). Chinese Journal of Aeronautics, Vol. 26, Issue 6, 2013, p. 1586-1595.
  • Shipanob G. V. Theory and methods for the design of automatic controls. Automatica i Telemechanica, Vol. 1, 1939, p. 49-66, (in Russian).
  • Fang Q., Yao Y., Wang X. C. Disturbance observer design for electric aerodynamic load simulator. Proceedings of 2005 International Conference on Machine Learning and Cybernetics, IEEE, Vol. 2, 2005, p. 1316-1321.
  • Wang M., Ben G., Yudong G., et al. Design of electric dynamic load simulator based on recurrent neural networks. Electric Machines and Drives Conference, IEEE, Vol. 1, 2003, p. 207-210.
  • Truong D. Q., Ahn K. K. Force control for press machines using an online smart tuning fuzzy PID based on a robust extended Kalman filter. Expert Systems with Applications, Vol. 38, Issue 5, 2011, p. 5879-5894.
  • Wang M., Gun B. Electric servo load simulator based on iterative learning control. Proceedings of the CSEE, Vol. 23, Issue 12, 2003, p. 123-126.
  • Wu S., Er M. J. Dynamic fuzzy neural networks-a novel approach to function approximation. IEEE Transactions on Systems Man and Cybernetics. Part B Cybernetics a Publication of the IEEE Systems Man and Cybernetics Society, Vol. 30, Issue 2, 2000, p. 358-364.
  • Ni Z., Wang M. A method of restraining surplus torque based on dynamic fuzzy neural network. Journal of Harbin Institute of Technology, Vol. 44, Issue 10, 2012, p. 79-83.
  • Zhang H., Xie L., Soh Y. C. H∞ deconvolution filtering, prediction, and smoothing: a Krein space polynomial approach. Signal Processing IEEE Transactions, Vol. 48, Issue 3, 2000, p. 888-892.
  • Rho H., Hsu C. S. A game theory approach to H∞ deconvolution filter design. Proceedings of the American Control Conference, IEEE, Vol. 4, 1999, p. 2891-2895.
  • Åström K. J., Wittenmark B. Computer-Controlled Systems: Theory and Design. Dover Publications Inc., New York, 1993.

About this article

08 April 2016
18 August 2016
30 September 2016
Oscillations in electrical engineering
position disturbance
load simulator
vector matching