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

1Corresponding author

Journal of Vibroengineering, Vol. 18, Issue 6, 2016, p. 4116-4129. https://doi.org/10.21595/jve.2016.17053
Received 8 April 2016; received in revised form 20 July 2016; accepted 18 August 2016; published 30 September 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 23
Abstract.

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.

Keywords: position disturbance, load simulator, vector matching.

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. 1. 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. 2. Schematic 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 load-carrying 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 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,

${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,

${K}_{\mathrm{t}}$ – Torque coefficient,

$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. 3. 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 ${\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, 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 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 ${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. 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 ${T}_{\mathrm{L}}$ is:

(1)
${T}_{\mathrm{L}}={W}_{\mathrm{M}}\left(s\right){T}_{\mathrm{c}}+{W}_{\mathrm{R}}\left(s\right){\theta }_{\mathrm{c}}.$

In the following, ${W}_{\mathrm{M}}\left(s\right)$ and ${W}_{\mathrm{R}}\left(s\right)$ are called non-disturbing load transfer function and position disturbance transfer function, respectively. Their frequency domain expressions are described as:

(2)
$\left\{\begin{array}{l}{W}_{\mathrm{M}}\left(j\omega \right)={A}_{\mathrm{M}}\left(\omega \right){e}^{{j\phi }_{\mathrm{M}}\left(\omega \right)},\\ {W}_{\mathrm{R}}\left(j\omega \right)={A}_{\mathrm{R}}\left(\omega \right){e}^{{j\phi }_{\mathrm{R}}\left(\omega \right)}.\end{array}\right\$

Assuming there is a sinusoidal position disturbance $\mathrm{sin}\left(\omega t\right)$. When the loading command is equal to 0 Nm, the output torque is only , and . 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:

(3)

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 $v\left(t\right)$, as Eq. (4). The sampling frequency is ${f}_{\mathrm{s}}$. A set of $\mathrm{M}$ sampling points are obtained.

(4)
$\stackrel{^}{v}\left(t\right)=A\mathrm{sin}\left(\omega t+\varphi \right)+W.$

Eq. (4) can be expanded to:

(5)
$\stackrel{^}{v}\left(t\right)={C}_{0}\mathrm{sin}\left(\omega t\right)+{C}_{1}\mathrm{cos}\left(\omega t\right)+{C}_{2}.$

In Eq. (5), ${C}_{0}=A\mathrm{cos}\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)$, , , $k=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 $\stackrel{^}{v}$ is:

(6)
$\delta \left(k\right)=v\left(k\right)-{C}_{0}\left(k\right){D}_{0}\left(k\right)-{C}_{1}\left(k\right){D}_{1}\left(k\right)-{C}_{2}\left(k\right).$

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

(7)
$\left[{\delta }^{2}\right]=\sum _{k=1}^{M}{\left[v\left(\mathrm{k}\right)-\stackrel{^}{\mathrm{v}}\left(\mathrm{k}\right)\right]}^{2}.$

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

(8)
$\left\{\begin{array}{l}{\stackrel{^}{\mathbf{\theta }}}^{\mathrm{T}}\left(k+1\right)={\stackrel{^}{\mathbf{\theta }}}^{\mathrm{T}}\left(k\right)+\mathbf{K}\left(k\right)\left[v\left(k+1\right)-{\mathbf{h}}^{\mathrm{T}}\left(k+1\right){\stackrel{^}{\mathbf{\theta }}}^{\mathrm{T}}\left(k\right)\right],\\ \mathbf{K}\left(k\right)=\mathbf{P}\left(k\right)\mathbf{h}\left(k+1\right){\left[{\mathbf{h}}^{\mathrm{T}}\left(k+1\right)\mathbf{P}\left(k\right)\mathbf{h}\left(k+1\right)+1\right]}^{-1},\\ \mathbf{P}\left(k+1\right)=\left[\mathbf{I}-\mathbf{K}\left(k\right){\mathbf{h}}^{\mathrm{T}}\left(k+1\right)\right]\mathbf{P}\left(k\right).\end{array}\right\$

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:

(9)
$\left\{\begin{array}{l}A=\sqrt{{{C}_{0}\left(\mathrm{M}\right)}^{2}+{{C}_{1}\left(\mathrm{M}\right)}^{2}},\\ \varphi =\mathrm{arctan}\left({C}_{1}\left(\mathrm{M}\right)/{C}_{0}\left(\mathrm{M}\right)\right).\end{array}\right\$

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 .

(10)
${A}_{2}=\sqrt{{A}_{0}^{2}+{A}_{1}^{2}-2{A}_{0}{A}_{1}\mathrm{cos}\left(\beta -\phi \right)}.$
(11)
$\theta =\mathrm{s}\mathrm{g}\mathrm{n}\left(\beta -\phi \right)×\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\left(\frac{{A}_{0}^{2}+{A}_{1}^{2}-{A}_{2}^{2}}{2{A}_{0}{A}_{1}}\right).$

So, ${T}_{\mathrm{c}0}$ needed by step 3 of vector matching method is described as follows:

(12)
${T}_{\mathrm{c}0}=\frac{{A}_{\mathrm{c}0}^{\text{'}}{A}_{0}}{{A}_{1}}\mathrm{sin}\left(\omega t+{\varphi }_{\mathrm{c}0}^{\text{'}}+\theta \right).$

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

– Position command complex vector in step 1,

– Position command complex vector in step 2,

– Position command complex vector in step 3,

– Compensating command complex vector in step 2,

– Compensating command complex vector in step 3,

– Output load torque complex vector in step 1,

– Output load torque complex vector in step 2,

– Relative load torque complex vector in step 1,

– Relative load torque complex vector in step 2,

– Pure load torque complex vector in step 2,

– 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:

(13)
${\mathbf{T}}_{\mathrm{L}2}={W}_{\mathrm{R}}\left(jk{\omega }_{0}\right){\mathbf{\theta }}_{\mathrm{c}2}+{W}_{\mathrm{M}}\left(jk{\omega }_{0}\right){\mathbf{T}}_{\mathrm{c}2}$

$k=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.

(14)
$0={W}_{\mathrm{R}}\left(jk{\omega }_{0}\right){\mathbf{\theta }}_{\mathrm{c}3}+{W}_{\mathrm{M}}\left(jk{\omega }_{0}\right){\mathbf{T}}_{\mathrm{c}3}.$

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:

(15)
${\mathbf{T}}_{\mathrm{c}3}=-\frac{{W}_{\mathrm{R}}\left(jk{\omega }_{0}\right)}{{W}_{\mathrm{M}}\left(jk{\omega }_{0}\right)}{\mathbf{\theta }}_{\mathrm{c}3}=-\frac{\frac{{\mathbf{T}}_{\mathrm{L}1}}{{\mathbf{\theta }}_{\mathrm{c}1}}{\mathbf{T}}_{\mathrm{c}2}}{\left(\frac{{\mathbf{T}}_{\mathrm{L}2}}{{\mathbf{\theta }}_{\mathrm{c}2}}-\frac{{\mathbf{T}}_{\mathrm{L}1}}{{\mathbf{\theta }}_{\mathrm{c}1}}\right)}{\mathbf{\theta }}_{\mathrm{c}3}=-\frac{{\mathbf{H}}_{1}{\mathbf{T}}_{\mathrm{c}2}}{\left({\mathbf{H}}_{2}-{\mathbf{H}}_{1}\right){\mathbf{\theta }}_{\mathrm{c}2}}{\mathbf{\theta }}_{\mathrm{c}3}=-\frac{{\mathbf{H}}_{1}{\mathbf{T}}_{\mathrm{c}2}}{\mathbf{E}{\mathbf{\theta }}_{\mathrm{c}2}}{\mathbf{\theta }}_{\mathrm{c}3}=\mathbf{A}{\mathbf{\theta }}_{\mathrm{c}3}.$

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 non-zero components only exist on the harmonic frequencies of the periodic signal. Hence the actual calculation procedure of CVM is as follows:

(16)
${H}_{ik}=\left\{\begin{array}{ll}\frac{{T}_{\mathrm{L}ik}}{{\theta }_{\mathrm{c}ik}},& {\theta }_{\mathrm{c}ik}\ne 0,\\ 0,& {\theta }_{\mathrm{c}ik}=0,\end{array}\right\$
(17)
${E}_{k}={H}_{2k}-{H}_{1k},$
(18)

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. 5. Photos of servo motor limit capacity test bench

Table 1. Arrangement of these experiments

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. 6. 5 Hz sinusoid position disturbance constant torque load curve

Fig. 7. 10 Hz sinusoid position disturbance constant torque load curve

Fig. 8. 15 Hz sinusoid position disturbance constant torque load curve

Table 2. The 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

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.

(19)
$x\left(t\right)=\frac{8{A}_{\mathrm{t}\mathrm{r}\mathrm{i}}}{{\pi }^{2}}\left[\mathrm{cos}\left({\omega }_{0}t\right)+\frac{1}{{3}^{2}}\mathrm{cos}\left(3{\omega }_{0}t\right)+\frac{1}{{5}^{2}}\mathrm{cos}\left(5{\omega }_{0}t\right)+\cdots \right]=\frac{8{A}_{\mathrm{t}\mathrm{r}\mathrm{i}}}{{\pi }^{2}}\sum _{n=1}^{\infty }\left[\frac{1}{{n}^{2}}\mathrm{cos}\left(n{\omega }_{0}t\right)\right]$
$n=1,3,5,7,\cdots .$

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

Fig. 10. 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 ${\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. 11. Spectrum of triangular disturbance rejection experimental data

a) Position command

b) Original torque error

c) Matching torque command

d) Residual torque error

Table 3. Rejection 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

Fig. 14. 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.

References

1. 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. [Publisher]
2. 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. [Search CrossRef]
3. Nam Y., Hong S. K. Force control system design for aerodynamic load simulator. Control Engineering Practice, Vol. 10, Issue 5, 2002, p. 549-558. [Publisher]
4. 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. [Search CrossRef]
5. 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. [Publisher]
6. 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. [Search CrossRef]
7. 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. [Search CrossRef]
8. 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. [Publisher]
9. Shipanob G. V. Theory and methods for the design of automatic controls. Automatica i Telemechanica, Vol. 1, 1939, p. 49-66, (in Russian). [Search CrossRef]
10. 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. [Publisher]
11. 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. [Search CrossRef]
12. 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. [Publisher]
13. 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. [Search CrossRef]
14. 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. [Search CrossRef]
15. 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. [Search CrossRef]
16. 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. [Publisher]
17. 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. [Search CrossRef]
18. Åström K. J., Wittenmark B. Computer-Controlled Systems: Theory and Design. Dover Publications Inc., New York, 1993. [Search CrossRef]

Cited By

 Measurement Zurong Qiu, Jie Xue 2021