Study of the motion of an electromechanical system in the presence of an elastic link and delay in the negative feedback loop of a servo drive
Sergey Jatsun^{1} , Peter Bezmen^{2} , Andrey Yatsun^{3}
^{1, 2, 3}Southwest State University, Kursk, Russian Federatio
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 8, 2016, p. 403-408.
Received 22 July 2016; accepted 3 September 2016; published 7 October 2016
JVE Conferences
This article presents the results of the study of the motion of an electromechanical system in the presence of an elastic link and delay in the negative feedback loop of a servo drive. We propose a model of the drive with a digital angular position control system of the working body. The servo drive model takes into account the dynamics of the processes occurring in the collector of the DC motor, stiffness of the elastic link that transfers torque from the electric motor to the working body, the influence of the digital sensor of the angular coordinate of the working body and delay in the digital control system on the operation of the drive. Results of mathematical modeling of the operation of the studied model are presented.
Keywords: electromechanical system, mechatronics, servo drive, digital control system, dynamics.
1. Introduction
The use of digital electronic components (microcontrollers, digital sensors of physical quantities) in drives transforms these systems into the digital servo systems. The digital servo system is an example of a nonlinear impulse automatic control system [1-3]. Digital systems operate in the converter mode “digital code to physical quantity” (e.g., position or speed of the current state of the controlled system) and the converter mode “current status of control object to digital code”. An important criterion for evaluating the performance of digital servo systems is the value of the period of delay time in the feedback loop, caused by the operation of digital electronic components in the following processes: analogue-to-digital and digital-to-analogue processing of signal, arithmetic, logic and tabular conversions of a digital code (for example, conversion of Gray code into straight binary code), generation of a pulse-width modulation (PWM) signal for the control of the servo drive actuators and implementation of service procedures of the control system (for example, output of data about the current state of the system). In other words, the delay period is the time spent by the digital components to process signals. The smaller the delay period, the greater the speed of the system.
Sensitivity, range of measurement and bandwidth of the sensors affect the accuracy of the operation of servo drives. When digital sensors are used in servo systems the discreteness of the output digital signal and sampling rate also affect the operation of the systems. The discreteness of the output signal of the sensor is determined by resolution expressed in units and tens of bits. The use of a digital sensor with low resolution in the servo system reduces its sensitivity and can cause oscillatory processes [4-6].
In this work we present a model of an electric servo drive with a digital control system of the working body’s angular coordinate that takes into account the stiffness of the shaft, which transfers torque from the electric motor to the working body and delay in the negative feedback loop of the electric servo drive.
2. Description of the object of study
The object of study is the influence of the actions of delay in the negative feedback loop, stiffness of the elastic link that transfers torque from the motor to the working body and the size of dead zones of digital sensors due to their resolution capabilities on the operation of the servo system of the electromechanical drive (Fig. 1). Thus, we consider a servo system with a DC collector electric motor (EM). The reference signal former, digital signal processing block, digital comparator and digital regulator block (Fig. 1) are programmatically implemented on the microcontroller [7].
The reference signal former (Fig. 1) generates a certain value of the required angular displacement of the working body, ${\stackrel{-}{\phi}}_{21}=f\left(t\right)$. The task of the sensor – digital encoder (DE) is to measure the angular displacement of the working body of the actuating mechanism, ${\phi}_{21}^{*}$.
The difference:
is the position (angular coordinate) error, $\epsilon $. The task of the servo system is to minimize this error i.e. $\epsilon \to $ 0.
Fig. 1. Diagram of the servo system of the electric drive during control of the angular coordinate of the working body
We assume that the preliminary digital signal processing block with a time delay of ${\tau}_{1}$, the digital comparator with a time delay of ${\tau}_{2}$, the digital regulator block with a time delay of ${\tau}_{3}$ and the reference signal former with a time delay of ${\tau}_{4}$ (Figs. 1, 2) are units of pure time delay.
We take ${\tau}_{1}=$ 10·10^{-3} s (when we use an absolute 8-bit encoder Bourns ACE-128 is connected to the ${I}^{2}C$ interface of a microcontroller – Atmel ATmega32u4 via an input-output expander microchip Microchip MCP23008 and using a data rate of 100 kb/s) for the unit of the preliminary digital signal processing block (including the time taken to receive the sensor signal).
For the digital comparator we take ${\tau}_{2}=$ 10·10^{-6} s, for the digital regulator block – the value ${\tau}_{3}=$ 10·10^{-6} s and for the reference signal former – the value ${\tau}_{4}=$ 10·10^{-6} s.
Since the maximal value of delay is period ${\tau}_{1}$ of the unit of the preliminary digital signal processing block, we should take period ${\tau}_{1}$ as the overall delay period in the negative feedback loop of the servo drive. All units operate asynchronous and simultaneously.
The period of delay for digital devices amounts to tens – hundreds of milliseconds.
The functional diagram of the operation of the servo drive with the angular coordinate control of the working body is shown in Fig. 2.
Encoders are one of the digital devices for measuring the angular coordinate. This type of sensor is a coded position sensor. The majority of the commercially available code encoders use the Gray code or similar codes, in which the reading error does not exceed the value of the least significant bit of the angle code. The disadvantage of the sensors using the scale with Gray code, or similar, is the need for decryption of the special code into straight binary code that is proportional to the measured angular coordinate.
Fig. 2. Functional diagram of the operation of the servo drive with the angular coordinate control of the working body
In order to convert a special code into straight binary code, it is necessary to use table lookup, as a result of which each value of the special code, ${\phi}_{21rd}^{*}$, is assigned to a binary code $N$. The shaft of the encoder, Bourns ACE-128, can occupy 128 different positions that correspond to the special code. Thus, the length of the binary code, Bourns ACE-128, obtained after table lookup, is 7 bits.
The expression for the conversion of encoder binary code $N$ into the encoder shaft rotation angle:
where: ${\phi}_{21}^{**}$ is the angle of rotation of the encoder shaft in degrees, $N$ – the value of binary code (number of the shaft’s position).
One of the methods of varying the supply voltage, ${U}_{21}$ is the application of pulse-width modulation (PWM).
We define the pulse ratio of the PWM signal as the ratio of the duration of the electric pulses, μ to the pulse repetition interval, $T$. The pulse ratio, $S$ determines the relation between the peak and average values of voltage and current pulses.
The value of the pulse ratio is less or equal to unity:
A PWM signal, ${u}_{21}$ with pulse ratio, $S$ enters the input of the power amplifier (Figs. 1, 2) that converts it into the supply voltage, ${U}_{21}$ of the motor. The supply voltage, ${U}_{21}$ is maximum when $S=$ 1 and consequently the angular velocity of the shaft’s rotation, ${\dot{\phi}}_{21}$ is also maximum. The value of ${\dot{\phi}}_{21}$ is directly proportional to the value of pulse ratio of the PWM signal. Further, under the PWM, ${u}_{21}$ we imply its pulse ratio, $S$.
3. Mathematical model
The system of equations describing the operation of the electromechanical system (Fig. 1) taking into account the presence of an elastic link, time delay in the negative feedback loop of the servo drive and the digital encoder with a dead zone can be presented in the form of system of Eq. (4):
where: ${\phi}_{21}$ – angular coordinate of rotation of the motor shaft of the electromechanical drive, rad; ${\phi}_{21}^{*}$ – angular coordinate of rotation of the mechanism’s working body, rad; ${J}_{EM}$ – total moment of inertia on the axis of the motor, kg·m^{2}; ${m}_{1}$ – mass of the working body of the actuating mechanism, kg; ${l}_{1}$ – length of the working body of the actuating mechanism, m; $J$ – moment of inertia of the actuating mechanism, kg·m^{2}, $J=\left(1/3\right){m}_{1}{l}_{1}^{2}$; $L$ and ${r}_{r}$ – inductance, $H$ and resistance, Ohm of the electromechanical drive motor respectively; ${M}_{EM}$ – torque of the electromechanical drive motor, N·m; ${i}_{r}$ – current of the armature of the electromechanical drive motor, А; ${C}_{E}$ and ${C}_{M}$ – coefficients of proportionality electric motor of the electromechanical drive; ${U}_{21}$ – control voltage of the motor, V; $\epsilon $ – error in the angular coordinate of rotation of the working body, rad; ${u}_{21}$ – pulse ratio of the PWM signal; ${\stackrel{-}{\phi}}_{21}$ – value of the required angular displacement of the working body, rad; $t$ – time, s; ${\tau}_{1}$ – value of the time delay period of the digital signal processing block, s; ${\tau}_{2}$ – value of the time delay period of the of the digital comparator, s; ${\tau}_{3}$ – value of the time delay period of the digital regulator block, s; ${\tau}_{4}$ – value of the time delay period of the reference signal former, s; ${\phi}_{21prev.}^{**}$ – angular coordinate of the rotation of the working body of the actuating mechanism measured by the encoder and calculated at $\mathrm{m}\mathrm{o}\mathrm{d}\left(t,{\tau}_{1}\right)=$ 0, rad; ${\epsilon}_{prev}$_{.} – error in the angular coordinate of rotation of the working body evaluated at $\mathrm{m}\mathrm{o}\mathrm{d}\left(t,{\tau}_{2}\right)=$ 0, rad; ${u}_{21prev}$_{.} – pulse ratio of the PWM signal evaluated at $\mathrm{m}\mathrm{o}\mathrm{d}\left(t,{\tau}_{3}\right)=$ 0; ${\stackrel{-}{\phi}}_{21prev.}$ – value of the required angular displacement of the working body, evaluated at $\mathrm{m}\mathrm{o}\mathrm{d}\left(t,{\tau}_{4}\right)=$ 0, rad; $c$ – spring constant of the elastic link between the electric motor and the working body, N/m; $P$ – proportional regulator coefficient; ${\phi}_{S}$ – size of the dead zone of the digital encoder, ${\phi}_{S}=$ 0,049 rad; $\delta $ – number of the different positions for the encoder, Bourns ACE-128, $\delta =$ 128; $N$ – number of current positions of the encoder shaft, $N\in {N}_{0}$, $N=$ 0…$\delta -$1.
4. Results
To obtain the data describing the operation of the electromechanical system the mathematical modelling was implemented in MathWorks MATLAB using the system of Eq. (4). The mathematical modelling process was carried out at following values of the model parameters: ${r}_{r}=$ 1 Ohm, $L=$ 1·10^{-6} H, ${J}_{EM}=$ 1·10^{-6} kg·m^{2}, ${U}_{max}=$ 12 V, ${\omega}_{idle}=$ 209 rad/s, ${M}_{nom}=$ 0.6025 N·m, ${i}_{r.nom}=$ 1 A, $c=$ 6.025 N/m, $P=$ 0.1, ${\tau}_{1}=$ 10·10^{-3} s, ${\tau}_{2}=$ 10·10^{-6} s, ${\tau}_{3}=$ 10·10^{-6} s, ${\tau}_{4}=$ 10·10^{-6} s, $\delta =$ 128, ${\phi}_{S}=$ 0,049 rad, ${m}_{1}=$ 0.1 kg, ${l}_{1}=$ 0.2 m, the total period of mathematical modelling is 2.5 s and the fixed time step – 1·10^{-6} s.
The results of mathematical modelling are presented in Fig. 3-7.
Fig. 3. The diagrams ${\stackrel{-}{\phi}}_{21}\left(t\right)$, ${\phi}_{21}^{**}\left(t\right)$ and $\epsilon \left(t\right)$
Fig. 4. The diagrams ${\phi}_{21}\left(t\right)$, ${\dot{\phi}}_{21}\left(t\right)$ and ${\ddot{\phi}}_{21}\left(t\right)$
Fig. 5. The diagrams ${\phi}_{21}^{*}\left(t\right)$, ${\dot{\phi}}_{21}^{*}\left(t\right)$ and ${\ddot{\phi}}_{21}^{*}\left(t\right)$
Fig. 6. The diagrams of ${\epsilon}_{max}\left({\tau}_{1}\right)$ at: 1) ${\tau}_{2}={\tau}_{3}={\tau}_{4}=$ 1·10^{-6} s, 2) ${\tau}_{2}={\tau}_{3}={\tau}_{4}=$ 1·10^{-3} s, 3) ${\tau}_{2}={\tau}_{3}={\tau}_{4}=$ 10·10^{-3} s
Fig. 7. The diagrams of ${\epsilon}_{max}\left(c\right)$ at: 1) ${\tau}_{1}={\tau}_{2}={\tau}_{3}={\tau}_{4}=$ 1·10^{-6} s, 2) ${\tau}_{1}={\tau}_{2}={\tau}_{3}={\tau}_{4}=$1·10^{-3} s, 3) ${\tau}_{1}={\tau}_{2}={\tau}_{3}={\tau}_{4}=$ 10·10^{-3} s
5. Conclusions
The results of mathematical modelling show that an important criterion for evaluating the performance of the digital servo system is the value of the period ${\tau}_{1}$ of delay in the feedback loop. The smaller delay period ${\tau}_{1}$ present, the greater the operating precision of the digital servo system. The repetition rate of the error ε peaks (Fig. 3) depends on the period ${\tau}_{1}$ of delay (${\tau}_{1}=$ 10·10^{-3} s) in the feedback loop. The absolute maximal error $\epsilon $ value is determined by the size of the dead zone of the digital encoder (${\phi}_{S}$), and the angular velocity of rotation of the mechanism’s working body. If the delay periods ${\tau}_{1}$, ${\tau}_{2}$, ${\tau}_{3}$, ${\tau}_{4}$ by orders of magnitude are smaller than the rotation period of the mechanism’s working body, it means that the servo drive can function stably. Figs. 4, 5 demonstrate the time dependences of the mechanical characteristics of the rotation of the motor shaft of the electromechanical drive and rotation of the mechanism working body. These mechanical characteristics are significantly dependent on the delay in the feedback loop. The angular velocity and angular acceleration characteristics are susceptible to the delay in the feedback loop, ${\tau}_{1}$.
Fig. 6 presents the dependence of the absolute maximum value of ε on the value of ${\tau}_{1}$ at the accepted model parameters. The diagram ${\epsilon}_{max}\left({\tau}_{1}\right)$ disposition is explained by various combinations of the delay periods ${\tau}_{2}$, ${\tau}_{3}$ and ${\tau}_{4}$. The smaller the values of these delay periods as compared to ${\tau}_{1}$, the bigger the absolute maximal value of $\epsilon $ at identical values of ${\tau}_{1}$. If the delay period, ${\tau}_{1}$ is smaller than the delay periods, ${\tau}_{2}$, ${\tau}_{3}$ and ${\tau}_{4}$, the absolute maximal value of $\epsilon $ is invariable.
Fig. 7 shows the dependence of the absolute maximal ε value on the $c$ value at the accepted model parameters. The diagram ${\epsilon}_{max}\left(c\right)$ disposition is explained by various combinations of the delay periods ${\tau}_{1}$, ${\tau}_{2}$, ${\tau}_{3}$ and ${\tau}_{4}$. The smaller these delay periods are, the smaller the absolute maximal value of $\epsilon $ at identical values of $c$. Since the value of $c$ increases, it is clear that the absolute maximum value of $\epsilon $ tends to its minimum value that is provided by certain delay periods ${\tau}_{1}$, ${\tau}_{2}$, ${\tau}_{3}$ and ${\tau}_{4}$.
Acknowledgements
This work is financially supported by the Russian Science Foundation, Project No. 14-39-00008 “The Establishment of the Research Laboratory of Modern Methods and Robotic Systems to Improve the Human Environment”.
References
- Jatsun S. F. Locomotion control method for patients verticalization with regard to their safety and comfort. Proceedings of the 26th DAAAM International Symposium on Intelligent Manufacturing and Automation, 2015, p. 107-112. [CrossRef]
- Astrom K. J. Advanced PID Control. ISA – The Instrumentation, Systems, and Automation Society, 2006, p. 460. [CrossRef]
- Jatsun S. F. The modelling of the standing-up process of the anthropomorphic mechanism. Assistive Robotics: Proceedings of the 18th International Conference on CLAWAR, 2015, p. 175-182. [CrossRef]
- Panovko G. Y. Simulation of exoskeleton sit-to-stand movement. Journal of Machinery Manufacture and Reliability, Vol. 45, Issue 3, 2016, p. 206-210. [CrossRef]
- Jatsun S. F. Modelling of exoskeleton movement in verticalization process. New Developments in Pure and Applied Mathematics: Proceedings of the International Conference on Pure Mathematics – Applied Mathematics, Vienna, Austria, 2015, p. 83-87. [CrossRef]
- Jatsun S. Study of Controlled Motion of Exoskeleton Moving from Sitting to Standing Position. Advances in Robot Design and Intelligent Control, Springer International Publishing, 2016, p. 165-172. [CrossRef]
- Jatsun S. Parameter optimization for exoskeleton control system using Sobol sequences. Proceedings of the 21st CISM-IFToMM Symposium on Robot Design, Dynamics and Control, Springer International Publishing, 2016, p. 361-368. [CrossRef]
- Levin Mindy F. Progress in Motor Control – Skill Learning, Performance, Health, and Injury. Springer, 2014. [CrossRef]
- Phillips Charles L. Digital Control System Analysis and Design. 4 Edition, Pearson Publishing, 2014. [CrossRef]