RusAtom
youtube
РУС
Home Articles Development of a model of DC electric drive double-loop system based on REPEAT
08.06.2023
Development of a model of DC electric drive double-loop system based on REPEAT

Electric drive as a device is widely used in various applied and scientific fields. The most common ones are electrical machinery, radio electronics, automobile industry, automation engineering, and computer engineering. Multiloop circuits are used to simulate the operation of electric drives. The best and effective way is mathematical modeling of double-loop electric drive system, and REPEAT software being a model-oriented design and mathematical modeling environment successfully copes with these tasks of digital model creation. This paper describes the process of developing a DC double-loop electric drive system. Based on the given parameters, a suitable DC motor is selected. Using REPEAT software, a dynamic model of the electric drive is built based on the transfer characteristics of the electric motor, power supply (PS), current sensors (CT) and tachogenerator (TG) and a current and speed controller is synthesized. The simulation is examined using transient analysis and comparison with design values.


Input information for creating a simulation

Load moment of inertia Jl

215 kg ∙ m2

Static moment of load resistance Мl

145 N ∙ m

Angular velocity of load Ωl

44 deg/s

Required angular acceleration of load εl

9 deg/s2

Power supply transfer factor image001.png

24

Number of rectified voltage ripples per period m

3

Filter time constant

0.006 s

Input voltage of the current loop summing amplifierimage002.png

8 V

Current sensor time constantimage003.png

0.008 s

Input voltage of the speed loop summing amplifierimage004.png

6 V

Tachogenerator time constantimage005.png

0.007 s

Frequency of inverter supply voltageimage006.png

400 Hz

Gearbox efficiencyimage007.png

0.85


Functional diagram of electric drive

The electric drive consists of a potentiometer, an amplifier, a DC motor (DC EM), a gearbox and a tachogenerator (see Figure 1).

The potentiometer is a control element. The control signals from the potentiometer output are transmitted to the input of the amplifier, with the armature winding of the DC EM with separate excitation being its load.

The EM turn the mechanism with angular speed proportional to the reference stimulus through the reducer.

image015.gif

Figure 1. Structural diagram of servo drive

The tachogenerator forms a rigid negative feedback (NFB) on angular velocity and ensures the signal generation at the amplifier input to deviate the load angular velocity from the design values.

This paper presents a closed-loop electric drive (ED) system based on the principle of REPEAT-based subordinate coordinate control. The system is multi-circuit, it consists of a speed loop and a current loop. Each loop of this system is configured separately.


Calculating of power and selecting of a DC motor

Based on the required parameters listed below, calculations are made for the selection of a DC EM and its dynamic model:

- required angular velocity of the load image008.png

- required angular acceleration of the load image009.png

- load inertia moment image010.png

- load static resistance torque image011.png;

- gearbox efficiency image012.png.

1. First, one should convert the angular velocity of the load rotation from “s” to “rad/s” and the angular acceleration of the load rotation from “deg/s2” to “rad/s2”:

image013.png

image014.png

The required power is calculated as follows:

image015.png

2. EM with shaft rated power greater than the required power (image016.png ˃ image017.png) is selected.

Relative to the calculated parameters, electric motor MI-22F is selected. This motor is a DC reversible parallel-excitation actuator motor. It is designed for operation in automatic control circuits. Specifications are given in Table 1.

Table 1. Specifications of EM MI-22F

Motor type

Shaft power Prat, kW

Rotation frequency nrat,

min-1

Supply voltage Urat, V

Armature current

 Ia, А

Armature winding resistance

Ra, Ohm

Rated torque Мrat,

N ∙ m

Inertia moment Jmtr∙10–4

kg ∙ m2

MI-22

0,37

3000

110

4,4

0,546

1,2

40,8

3. The best gear ratio  image037.gif is calculated using the formula:

image019.png

4. The selected EM is checked to ensure that it meets the speed requirements.

The following formula determines the rated angular velocity image020.png:

image021.png

and the load angular velocity reduced to the EM shaft:

image022.png

Since image023.png the speed requirements are not met.

     A new gear ratio is calculated using the formula:

image024.png

5. The required motor torque image025.png is calculated:

image026.png

image027.png

6. The selected EM is checked for compliance with the torque (moment) requirements:

image028.png

image029.png

The torque requirements are met.

7. The parameters of the EM dynamic model are calculated.

   7.1. Counter EMF factor  image030.pngis determined:

image031.png

7.2. Torque factor image032.png is determined:

image033.png

7.3. The EM electromechanical time constant image034.pngis determined:

image035.png

7.4. To determine the EM electromagnetic time constant the armature inductance  image036.pngis calculated:

image037.png

La = 0.0022 H is adopted to calculate the electromagnetic time constant image038.png used in EM transfer function:

image039.png

Electric motor parameters and factors are calculated:

Electromechanical factor:

image040.png

Motor torque factor:

image041.png

Electromagnetic factor:

image042.png

Taking into account the obtained numerical values of the dynamic model structural diagram, the EM DSSD appears as follows (see Figure 2):

image043.png

Figure 2. Structural diagram of dynamic model of electric motor with numerical values


Transient characteristics

Below the results of simulation with the following input parameters of EM Mi-22 are given (see Figure 3 and Figure 4).

- supply voltage Urat = 110 V;

- rated torque Mrat = 1.2 image044.png

image089.jpg

Figure 3. Reference stimulus transient characteristic of EM MI-22

image091.jpg

Figure 4. Disturbance transient characteristic of EM MI-22


Analysis of transients

Relative error image093.gif is determine by the formula:

image048.png

The data sheet rated value of angular rotation speed of EM MI-22 is image049.pngand slightly differs from the simulation results (323.5 rad/s). It allows to conclude that the performed calculations are correct.


Current loop setting

The dynamic properties of PS are described highly accurately by a lag with a transfer function image050.png:

image051.png

The dynamic properties of CT are also described by a lag image052.png

image053.png

We can now present a dynamic simulation of the current loop based on REPEAT software (see Figure 5).

image054.png

Figure 5. Structural diagram of the dynamic current loop simulation

Based on the structural diagram of the current loop simulation, the transfer function of the open current loop is found image055.png:

image056.png

The current loop (CL) shall be adjusted to optimum modulo (OM). CL transfer function adjusted to OM:

image057.png

where image058.png- is a total CL short time constant. To find it, a PS time constant image059.png is determined::

image060.png

Then

image061.png

To find the transfer function of the current controller (CC) it is necessary to equate the right parts of equations (3.3) and (3.4):

image062.png

In terms of structure, the obtained equation is a transfer function of PI-controller:

image063.png

Based on comparison of the last two formulas, we can obtain formulas to calculate a transfer factor Kcc and a time constant Тcc of СС:

image064.png

image065.png

To calculate the current transducer transfer ratio, use the formula:

image066.png

To build a current loop according to (3.6), the current controller transfer function is calculated:

image067.png


Transient characteristic

Preamble with indication (see Figure 6).

image135.jpg

Figure 6. Transient characteristic of the current loop according to the reference stimulus (a CL reference stimulus is 8 V)


Transient analysis

Overshoot σcl and buildup time image069.pngare determined..

The overshoot was calculated using the formula:

image070.png

The expected overshoot image141.gif when adjusting to OM image143.gif

Based on Fig. 7 a maximum deviation of the armature current image145.gifis determined:

image074.png

and the steady state current value image075.png:

image076.png

Overshoot is determined:

image077.png

The overshoot corresponds to the correct setting.

The buildup time image069.png is determined in the first point of intersection of the transient function graph and the armature steady state current value:

image078.png.

Based on the graph (Fig. 6) a buildup time image069.pngis determined:

image079.png

The buildup time can be calculated and should meet the requirement:

image080.png.

Deviation of measured and calculated value of image069.png:

image081.png

The deviation is small and acceptable. One can conclude that the adjustment of the current loop to the optimum modulo was successful.


Adjustment of speed loop to symmetrical optimum

The speed loop (SL) consists of the following elements:

-       speed controller (SC) F1;

-       CL adjusted to OM;

-       motor electromechanical part;

-       tachogenerator (TG).

Before proceeding to the building of CL diagram, it is necessary to represent all loop elements in as links.

Dynamic properties of the TG (inertial link):

image082.png

Adjustment of speed loop to symmetrical optimum (SO).

When adjusting SL to SO, due to short time constants Tct andimage083.png transfer function of the current loop image084.png


image085.png

Transfer function of open SL image086.png:

image087.png

where image088.png– is a SL total short time constant. image089.png

The transfer function of SL image086.png adjusted to SO is as follows:

image090.png

To determine the SC structure it is necessary to equate the right parts of equations (4.3) and (4.4):

image091.png

Let’s introduce designations to obtain the PI controller formula:

SC transfer ratio image092.png:

image093.png

SC time constant:

image094.png

Then the transfer function of the speed loop controller will be as follows:

image095.png

The required TG transfer ratio is calculated using the formula:

image096.png

The transfer function image097.png of the speed controller is calculated using the formula (4.8):

image098.png

A gearbox gear ratio image099.png is calculated:

image101.png

image201.jpg

Figure 7. Structural diagram of the dynamic speed loop simulation

Transient characteristics

Let’s build the output characteristics and analyze the results.

1. Building of SL transient characteristics based on reference stimulus. The value of SL reference stimulus is image103.png= 6 В.

image205.jpg

Figure 8. Reference stimulus transient characteristic of speed loop

2. Building of disturbance transient characteristic of SL (static moment of load resistance). Load Ml = 145 N/m.

image207.jpg

Figure 9. Disturbance transient characteristic of speed loop

The steady-state error component image209.gif is reduced to zero.

Analysis of transients

Maximum deviation of the angular speed of motor shaft rotation image107.png and a steady-state value image108.png

Expected overshoot image141.gif when adjusting to SO image215.gif

Overshoot σsl:

image110.png

Deviation from the expected value is image219.gif and is adopted as acceptable.

and is adopted as acceptable. image112.png:

image113.png

Calculated value image112.png:

image114.png

Deviation of measured and calculated value image112.png:

image115.png

Deviation of image112.png is acceptable. The speed loop adjustment to the optimum modulo is successful.


Simulation results

This paper describes the development of the DC servo drive dynamic model and the synthesis of current and speed controller (double-loop system) using standard tuning methods such as modulo optimum and symmetric optimum. In the course of work, a DC electric motor MI-22F was selected.

REPEAT software was successfully used for creating a digital model. The simulation results and calculated values showed that the speed and current loops were set correctly with acceptable deviation values.