Abstract
Position disturbance is the most troublesome problem for load simulator. It is difficult to work properly for a loading system, when the loadcarrying object is moving. Therefore, it is of high significance to compensate for the position disturbance. To achieve this objective, a kind of vector matchingbased 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 loadcarrying 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 hardwareintheloop simulation equipment [1]. For the most part, they have three components: a loading system, a loadcarrying object, a measuringcontrol system. The purpose of load simulation is to test the performance of a loadcarrying object under the condition as real as possible. It helps designers to foresee and detect problems related to mechanics and control of the loadcarrying 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 loadcarrying object will be moving according to the test signal. For example, in the field of aerospace, the rudder actuator is often implemented as a loadcarrying 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 loadcarrying 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 loadcarrying 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 loadcarrying 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 matchingbased 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 loadcarrying 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 measuringcontrol 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 measuringcontrol 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
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
The symbols in Fig. 2 are defined as follows:
${\mathbf{T}}_{0}$ – resultant torque vector measured by the torque sensor,
${\mathbf{T}}_{1}$ – torque vector produced by the movement of the loadcarrying object,
${\mathbf{T}}_{2}$ – torque vector produced by the loading system under the torque command ${\mathbf{T}}_{\mathrm{c}}$,
${\mathbf{T}}_{20}$ – one component of ${\mathbf{T}}_{2}$,
${\mathbf{T}}_{21}$ – another component of ${\mathbf{T}}_{2}$,
$\beta $ – phase angle of ${\mathbf{T}}_{0}$,
$\phi $ – phase angle of ${\mathbf{T}}_{1}$,
$\theta $ – angle between ${\mathbf{T}}_{1}$ and ${\mathbf{T}}_{2}$.
According to the invariance principle, ${\mathbf{T}}_{\mathrm{c}}$ should consist of two parts, i.e. ${\mathbf{T}}_{\mathrm{c}}={\mathbf{T}}_{\mathrm{c}0}+{\mathbf{T}}_{\mathrm{c}1}$. ${\mathbf{T}}_{\mathrm{c}0}$ is used to compensate for the position disturbance, and ${\mathbf{T}}_{\mathrm{c}1}$ is used to realize the loading. The purpose of vector matching is to make ${\mathbf{T}}_{20}\mathrm{}$equal to ${\mathbf{T}}_{1}$.When the load target is 0 Nm, ${\mathbf{T}}_{21}$ should be equal to 0, and ${\mathbf{T}}_{2}$ should be equal to ${\mathbf{T}}_{1}$.
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 $\theta $ 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:
${L}_{\mathrm{m}}$ – Motor stator cyclic inductance,
${R}_{\mathrm{m}}$ – Motor stator resistance,$\mathrm{}$
${J}_{\mathrm{m}}$ – Moment of inertia on the motor shaft,
${B}_{\mathrm{m}}$ – Viscous friction coefficient on motor shaft,
${K}_{\mathrm{e}}$ – Back EMF coefficient,$\mathrm{}$
${K}_{\mathrm{t}}$ – Torque coefficient,$\mathrm{}$
$U$ – Armature voltage,
$E$ – Back EMF,
${I}^{\mathrm{*}}$ – Current command,
$i$ – Armature current,
${T}_{\mathrm{L}}$ – Output load torque,
${\theta}_{\mathrm{o}\mathrm{u}\mathrm{t}}$ – Elastic deformation of the shaft.
Fig. 3VMBDR 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 ${\mathbf{T}}_{1}$ and ${\mathbf{T}}_{2}$.
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 ${T}_{\mathrm{c}1}$ to 0 Nm at first. Then, when the switch is at position 1, the compensating signal ${T}_{\mathrm{c}0}$ is zero, and the original position disturbance torque vector ${\mathbf{T}}_{1}$ could be obtained. As the switch is at position 2, the compensating signal ${T}_{\mathrm{c}0}$ is a filtered position command ${\theta}_{\mathrm{c}}$, and the resultant torque vector ${\mathbf{T}}_{0}$ could be gotten. Usually, the filter is a proportion link. So far, vector matching unit uses ${\mathbf{T}}_{0}$ and ${\mathbf{T}}_{1}$ to compute the loading torque vector ${\mathbf{T}}_{2}$ and $\theta $, then calculates a new compensating signal, which can eliminate ${\mathbf{T}}_{1}$ completely in step 3.
Changes in the environment will affect many physical parameters of ELS – motor stator resistance${R}_{\mathrm{m}}$, 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 ${\mathbf{T}}_{20}{\mathbf{T}}_{1}$ should be monitored in real time, and ${\mathbf{T}}_{20}$ should be finetuned 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 ${T}_{\mathrm{c}}$, position command ${\theta}_{\mathrm{c}}$ and output load torque ${T}_{\mathrm{L}}$ measured by the torque sensor, as shown in Fig. 3. ${T}_{\mathrm{c}}$ and ${\theta}_{\mathrm{c}}$ are generated by ELS controller computer itself. So only ${T}_{\mathrm{L}}$ could be contaminated by the additive noise. A Butterworth filter could be used to process the contaminated ${T}_{\mathrm{L}}$, when the noise is too large. However, the influence of the filter on the position disturbance torque vector ${\mathbf{T}}_{1}$ and the loading torque vector ${\mathbf{T}}_{2}$ are the same. Therefore, the matching relation between ${\mathbf{T}}_{1}$ and ${\mathbf{T}}_{2}$ does not change.
3. Vectormatching 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 ${T}_{\mathrm{L}}$ is:
In the following, ${W}_{\mathrm{M}}\left(s\right)$ and ${W}_{\mathrm{R}}\left(s\right)$ are called nondisturbing load transfer function and position disturbance transfer function, respectively. Their frequency domain expressions are described as:
Assuming there is a sinusoidal position disturbance $\mathrm{sin}\left(\omega t\right)$. When the loading command ${T}_{\mathrm{c}1}$ is equal to 0 Nm, the output torque is only ${T}_{1}$, and ${T}_{1}={A}_{\mathrm{R}}\left(\omega \right)\mathrm{sin}\left[\omega t+{\phi}_{\mathrm{R}}\left(\omega \right)\right]$. The total output ${T}_{\mathrm{L}}$ 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 online identification algorithm based on recursive least square method for the amplitude and phase of sinusoid. Moreover, the least square method also has excellent antinoise ability.
Assuming that there is a sinusoidal signal $v\left(t\right)$, as Eq. (4). The sampling frequency is ${f}_{\mathrm{s}}$. A set of $\mathrm{M}$ sampling points are obtained.
Eq. (4) can be expanded to:
In Eq. (5), ${C}_{0}=A\mathrm{cos}\varphi $, ${C}_{1}=A\mathrm{sin}\varphi $, ${C}_{2}=W$. Let ${D}_{0}\left(t\right)=\mathrm{sin}\left(\omega t\right)$, ${D}_{1}\left(t\right)=\mathrm{cos}\left(\omega t\right)$, ${v\left(t\right)=T}_{\mathrm{L}}\left(t\right)$, ${\widehat{\mathbf{\theta}}}^{\mathrm{T}}\left(k\right)=\left[{C}_{0}\left(k\right){C}_{1}\left(k\right){C}_{2}\left(k\right)\right]$, ${\mathbf{h}}^{\mathrm{T}}\left(k\right)=\left[{D}_{0}\left(\mathrm{k}\right)\mathrm{}{\mathrm{D}}_{1}\left(\mathrm{k}\right)1\right]$, $k=\mathrm{1,2},\cdots ,M$.
The parameters to be identified are ${C}_{0}$, ${C}_{1}$and ${C}_{2}$.The measured sample sequence is $v\left(1\right)$, $v\left(2\right)$, …, $v\left(M\right)$. Then, the residual error $\delta $ of $\widehat{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), $\mathbf{K}$ and $\mathbf{P}$ are respectively the recursive weighting factor vector and $\mathbf{K}$ update vector. $\mathbf{I}$ is a unit matrix. Finally, we can get the optimal value of ${C}_{0}\left(\mathrm{M}\right)$, ${C}_{1}\left(\mathrm{M}\right)$ and ${C}_{2}\left(\mathrm{M}\right)$. The amplitude and phase angle could be calculated as follows:
Take the time domain data of the torque vector ${\mathbf{T}}_{0}$ and ${\mathbf{T}}_{1}$ into the recursive formula. The amplitudes ${A}_{0}$, ${A}_{1}$ and phases $\phi $, $\beta $ could be calculated. The amplitude ${A}_{2}$ and angle $\theta $ could be gotten by cosine theorem, as Eq. (10) and Eq. (11). The direction of angle $\theta $ depends on $\left(\beta \phi \right)$, the value range of $\phi $ and $\beta $ is $\left[\u2013\pi ,\pi \right]$.
So, ${T}_{\mathrm{c}0}$ needed by step 3 of vector matching method is described as follows:
In Eq. (12), ${A}_{\mathrm{c}0}^{\text{'}}$ and ${\varphi}_{c0}^{\text{'}}$ are amplitude and phase of ${T}_{\mathrm{c}0}^{\mathrm{\text{'}}}$, 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 nonsinusoidal 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 offline 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
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.
${\mathbf{\theta}}_{\mathrm{c}1}={\left({\theta}_{\mathrm{c}11},{\theta}_{\mathrm{c}12},\cdots ,{\theta}_{\mathrm{c}1n}\right)}^{\mathrm{T}}$ – Position command complex vector in step 1,
${\mathbf{\theta}}_{\mathrm{c}2}={\left({\theta}_{\mathrm{c}21},{\theta}_{\mathrm{c}22},\cdots ,{\theta}_{\mathrm{c}2n}\right)}^{\mathrm{T}}$ – Position command complex vector in step 2,
${\mathbf{\theta}}_{\mathrm{c}3}={\left({\theta}_{\mathrm{c}31},{\theta}_{\mathrm{c}32},\cdots ,{\theta}_{\mathrm{c}3n}\right)}^{\mathrm{T}}$ – Position command complex vector in step 3,
${\mathbf{T}}_{\mathrm{c}2}={\left({T}_{\mathrm{c}21},{T}_{\mathrm{c}22},\cdots ,{T}_{\mathrm{c}2n}\right)}^{\mathrm{T}}$ – Compensating command complex vector in step 2,
${\mathbf{T}}_{\mathrm{c}3}={\left({T}_{\mathrm{c}31},{T}_{\mathrm{c}32},\cdots ,{T}_{\mathrm{c}3n}\right)}^{\mathrm{T}}$ – Compensating command complex vector in step 3,
${\mathbf{T}}_{\mathrm{L}1}={\left({T}_{\mathrm{L}11},{T}_{\mathrm{L}12},\cdots ,{T}_{\mathrm{L}1n}\right)}^{\mathrm{T}}$ – Output load torque complex vector in step 1,
${\mathbf{T}}_{\mathrm{L}2}={\left({T}_{\mathrm{L}21},{T}_{\mathrm{L}22},\cdots ,{T}_{\mathrm{L}2n}\right)}^{\mathrm{T}}$ – Output load torque complex vector in step 2,
${\mathbf{H}}_{1}={\left({H}_{11},{H}_{12},\cdots ,{\mathrm{H}}_{1\mathrm{n}}\right)}^{\mathrm{T}}$ – Relative load torque complex vector in step 1,
${\mathbf{H}}_{2}={\left({H}_{21},{H}_{22},\cdots ,{\mathrm{H}}_{2\mathrm{n}}\right)}^{\mathrm{T}}$ – Relative load torque complex vector in step 2,
$\mathbf{E}={\left({E}_{1},{E}_{2},\cdots ,{\mathrm{E}}_{\mathrm{n}}\right)}^{\mathrm{T}}$ – Pure load torque complex vector in step 2,
$\mathbf{A}={\left({A}_{1},{A}_{2},\cdots ,{A}_{n}\right)}^{\mathrm{T}}$ – Gain of complex vector matching used in step 3.
In step 1 the compensation command is zero, so ${\mathbf{T}}_{\mathrm{L}1}$ contains the position disturbance torque only. However, in step 2, the compensation command is ${\mathbf{T}}_{\mathrm{c}2}$, then ${\mathbf{T}}_{\mathrm{L}2}$ turns out to be:
$\mathbf{}\mathbf{}\mathbf{}\mathbf{}\mathbf{}\mathbf{}\mathbf{}={A}_{\mathrm{R}}\left(k{\omega}_{0}\right){A}_{\mathrm{\theta}2}\left(k{\omega}_{0}\right){e}^{j\left[{\phi}_{\mathrm{R}}\left(k{\omega}_{0}\right)+{\phi}_{\mathrm{\theta}2}\left(k{\omega}_{0}\right)\right]}+{A}_{\mathrm{M}}\left(k{\omega}_{0}\right){A}_{\mathrm{T}2}\left(k{\omega}_{0}\right){e}^{j\left[{\phi}_{\mathrm{M}}\left(k{\omega}_{0}\right)+{\phi}_{\mathrm{T}2}\left(k{\omega}_{0}\right)\right]},$
$k=\mathrm{1,2},\cdots ,n$
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 ${W}_{\mathrm{R}}\left(jk{\omega}_{0}\right)$ and ${W}_{\mathrm{M}}\left(jk{\omega}_{0}\right)$, 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 ${\mathbf{T}}_{\mathrm{c}3}$. Normally, rudder position command ${\mathbf{\theta}}_{\mathrm{c}3}$ is set in advance, according to experimental requirements. Once the complex gain vector $\mathbf{A}$ is determined, the compensation torque command ${\mathbf{T}}_{\mathrm{c}3}$ can be matched successfully.
The Fourier series of a periodic signal is a discrete sequence, and the nonzero components only exist on the harmonic frequencies of the periodic signal. Hence the actual calculation procedure of CVM is as follows:
$i=\mathrm{1,2},k=\mathrm{1,2},\cdots ,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 loadcarrying 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
The experiments including three subexperiments 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 
1  Sinusoid  Constant torque loading 
2  Triangle  Constant torque loading 
3  Sinusoid  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
Fig. 710 Hz sinusoid position disturbance constant torque load curve
Fig. 815 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 % 
5  0.48954  0.14467  70.45 
10  0.48782  0.17098  64.96 
15  0.72031  0.12330  82.88 
4.2. Triangle position – constant torque loading
A 5 Hz triangle signal $x\left(t\right)$, whose magnitude is represented by ${A}_{\mathrm{t}\mathrm{r}\mathrm{i}}$, is adopted as the position command for the tested motor in this experiment. The Fourier series of $x\left(t\right)$ 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.
$n=\mathrm{1,3},\mathrm{5,7},\cdots .$
Fig. 9Torque error caused by position disturbance of triangular wave
Fig. 10CVM 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 ${\mathbf{\theta}}_{\mathrm{c}3}$, the compensation torque command ${\mathbf{T}}_{\mathrm{c}3}$, 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
a) Position command
b) Original torque error
c) Matching torque command
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 % 
5  0.287  0.014  95.135 
15  0.452  0.040  91.230 
25  0.613  0.072  88.217 
30  0.507  0.028  94.407 
45  0.190  0.060  68.163 
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
Fig. 135 Hz spring torque loading after vector matching
Fig. 14Error 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 matchingbased 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 multiinput systems, vector matching method will be still effective in great probability.
References

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. 357362.

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. 369383.

Nam Y., Hong S. K. Force control system design for aerodynamic load simulator. Control Engineering Practice, Vol. 10, Issue 5, 2002, p. 549558.

Ryu H. M., Kim S. J., Sul S. K., et al. Dynamic load simulator for highspeed elevator system. Proceedings of the Power Conversion Conference, IEEE, Vol. 2, 2002, p. 885889.

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. 325332.

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. 120123.

Li Y. Development of hybrid control of electrohydraulic torque load simulator. Journal of Dynamic Systems, Measurement, and Control, Vol. 124, Issue 3, 2002, p. 415419.

Wang C., Jiao Z., Wu S., et al. An experimental study of the dualloop control of electrohydraulic load simulator (EHLS). Chinese Journal of Aeronautics, Vol. 26, Issue 6, 2013, p. 15861595.

Shipanob G. V. Theory and methods for the design of automatic controls. Automatica i Telemechanica, Vol. 1, 1939, p. 4966, (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. 13161321.

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. 207210.

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. 58795894.

Wang M., Gun B. Electric servo load simulator based on iterative learning control. Proceedings of the CSEE, Vol. 23, Issue 12, 2003, p. 123126.

Wu S., Er M. J. Dynamic fuzzy neural networksa 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. 358364.

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. 7983.

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. 888892.

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. 28912895.

Åström K. J., Wittenmark B. ComputerControlled Systems: Theory and Design. Dover Publications Inc., New York, 1993.