Next Article in Journal
An Adaptive Low Computational Cost Alternating Direction Method of Multiplier for RELM Large-Scale Distributed Optimization
Next Article in Special Issue
A Source Identification Problem in Magnetics Solved by Means of Deep Learning Methods
Previous Article in Journal
A Dynamic CGE Model for Optimization in Business Analytics: Simulating the Impact of Investment Shocks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Innovative Approach for the Determination of a DC Motor’s and Drive’s Parameters Using Evolutionary Methods and Different Measured Current and Angular Speed Responses

Faculty of Electrical Engineering and Computer Science, University of Maribor, Koroška Cesta 46, 2000 Maribor, Slovenia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(1), 42; https://doi.org/10.3390/math12010042
Submission received: 21 November 2023 / Revised: 12 December 2023 / Accepted: 21 December 2023 / Published: 22 December 2023
(This article belongs to the Special Issue Numerical Optimization for Electromagnetic Problems)

Abstract

:
The determination is presented of seven parameters of a DC motor’s drive. The determination was based on a comparison between the measured and simulated current and speed responses. For the parameters’ determination, different evolutionary methods were used and compared to each other. The mathematical model presenting the DC drives model was written using two coupled differential equations, which were solved using the Runge–Kutta first-, second-, third- and fourth-order methods. The approach allows determining the parameters of controlled drives in such a way that the controller is taken into account with the measured voltage. Between the tested evolutionary methods, which were Differential Evolution with three strategies, Teaching-Learning Based Optimization and Artificial Bee Colony, the Differential Evolution (DE/rand/1/exp) can be suggested as the most appropriate for the presented problem. Measurements with different sampling times were used, and it was found out that at least some measuring points should be at the speed-up interval. Different lengths of the measured signal were tested, and it is sufficient to use a signal consisting of the drive’s acceleration and a short part of the stationary operation. The analysis showed that the procedure has good repeatability. The biggest deviation of calculated parameters considering 10 repeated measurements was 6% in case of the La calculation. The deviations of all the other parameters’ calculations were less than 2%.

1. Introduction

Direct current motors (DC motors) are still used in many industrial applications. Their good features are a high starting torque and easy control. In many industrial applications, they are connected to power supply units that allow speed or torque control. Correct adjustment of the controllers requires knowledge of the electrical motor parameters. Usually, the electrical motor parameters are specified by the manufacturer. It may happen that these are not accurate enough or have large tolerances, especially in the case of cheaper electrical motors. It also happens that the parameters of older electrical motors are not known, but we need them due to the modernization of the drive in which the electrical motor is preserved, and the older power systems are replaced by more modern ones.
Knowledge of the inertia of the drive, friction, ventilation and load is required in order to determine the transition times and energy balance of the drive correctly. Even if the inertia of the motor is known, the inertia of the load and the friction of the load are often unknown. These data are important for correct dimensioning and optimization of the drive.
We can use a stand-alone program to identify the parameters, or it can be part of the software in the power supply. The stand-alone program offers the possibility of parameters’ determination for any drive with or without a power supply and in the case when we do not have access to the software in the power supply unit. The disadvantage of the stand-alone program is that it is necessary to transfer the measured data from the measuring devices, which is not the case if the parameters identification program is integrated into the power supply unit. If it is in the power supply, we want the identification to be as fast and robust as possible, so that it is not necessary to repeat the runs of the optimization method. According to what has been written, the choice of an appropriate optimization method is important. In this paper, we focus on evolutionary optimization methods. We focus on evolutionary methods, mainly because of their good properties. In the case of measured signals, measurement errors may occur, which could lead to a local minimum in the determination of parameters and thus to the determination of incorrect parameters. Evolutionary methods are capable of avoiding local minima, which is not the case for standard optimization methods. They are also robust and easy to use, because it is not necessary to determine the derivatives of the objective function, which can be difficult. or impossible, for many problems.
Other authors also tried to identify the parameters of the DC motor. They used different methods and approaches. In the literature, the following approaches can be found: the use of conjugate gradient and regularization methods [1], the use of curve fitting [2,3], the use of constraint optimization methods [4], the use of the regression method [5], the use of least-squares-based approaches [6,7,8,9,10] and the use of the moment method [11]. Approaches using evolutionary methods can also be found: the use of the Genetic algorithm [12] and the Multiobjective elitist genetic algorithm [13], the use of Differential Evolution [14], the use of the Flower pollination algorithm [15] and the use of the Particle Swarm Optimization [16].
This work presents the determination of seven parameters of a DC motor drive, which are the motor’s resistance Ra, the motor’s inductance La, the motor’s constant ce, the drive’s inertia J and the load, expressed with the parameters Tla, Tlb and Tlc. The parameters were determined based on the measured current and speed response of the DC motor. Parameters’ determination is an inverse problem. We used the direct approach. The mathematical model of the DC drive was constructed, and simulated current and speed responses were compared to measured values. Described is an optimization problem, where we are searching for the parameters for the mathematical model in such a way that the simulated responses fit with the measured values. The evolutionary optimization method was chosen for the parameter’s determination. Different evolutionary methods have different performance for different problems. Different ones have been tested, with the aim of determining a fast and robust method. The first one was Differential Evolution (DE) [14,17,18,19,20,21,22,23,24,25,26], and three different strategies were tested. Also tested were the state-of-the-art Artificial Bee Colony (ABC) [27,28,29,30,31,32,33,34] and Teaching Learning Based Optimization (TLBO) [35,36,37,38,39,40,41,42,43]. All three methods were compared with each other. In the case of the mathematical model, different numerical integrations were tested, and different sampling times and different measuring signals were tested in the case of the measured signals.
Our contributions in the presented work are:
  • The supply/control unit was considered in the mathematical model with a measured voltage. Not only the current and speed responses but also the armature supply voltage must be measured in our approach. In this way, no block model is needed of the supply/control unit. Such an approach was not found in the literature to the best of our knowledge.
  • Step response [2,3,4,5,13,14,15] or stepped response [1,9] was used in the literature. We made tests for speed-up (SU) and additionally also for speed-up + speed-down (SUD) responses.
  • We compared three evolutionary methods, DE (using three strategies), ABC and TLBO, with the aim to find the most suitable one. The methods were compared considering the SU and also SUD responses.
  • Different numerical integrations (first, second, third and fourth order) were tested considering precision and calculation time.
  • The results for different sampling times of the measured responses were tested to determine the appropriate sampling time according to the speed-up time.
  • The calculations were made in such a way that, for each sampling time, all four integration methods are used, and with that the relation was analyzed between the sampling time and the integration method.
  • The repeatability of the measurement was analyzed. For each of the 10 repeated measurements, the parameters were determined and compared with each other. In this way, it was determined if the measurement was sufficiently reliable that it was not necessary to repeat it.
  • Determination was carried out of the parameters for different length measured signals. These were full signal, signal without a pre-trigger, start only, start + same time operation, and start + operation for twice the start time.
In our work, the supply/control unit voltage was used as the input for the mathematical model; block diagrams of the drive were used in the related works [7,8,11,15,16]. An improved innovative approach is presented in this way. Other authors used a certain evolutionary method to determine the parameters, such as GA [12,13], DE [14] and PSO [16]. We used and compared the three methods with each other.
The paper consists of eight sections. The mathematical model and measurements are described in Section 2. The evolutionary methods and algorithm used are presented in Section 3. The test measurements, together with the parameters’ determination results, are presented in Section 4. The convergence of the evolutionary methods is also presented. Section 5 presents an analysis of the interaction between the order of the integration method and the measurements’ sampling time. The repeatability of the parameters’ determination is shown in Section 6. The parameters’ determination results for different lengths of the measured signals are presented in Section 7. Validation based on the simulated input data is shown in Section 8. The conclusions are given in the last section, Section 9.

2. Mathematical Model and Measurements

The parameters’ determination was based on the comparison between the measured responses and the simulated responses using evolutionary methods. An appropriate mathematical model of the DC motor’s drive was constructed to obtain simulated responses.

2.1. Mathematical Model

The model of the DC motor’s drive is presented schematically in Figure 1.
The DC motor’s drive is described using three coupled differential equations. The field circuit of the motor is described using (1), the armature circuit is described using (2) and the mechanical behavior of the drive is described using (3).
u f = i f R f + L f d i f d t
u a ( t ) = i a R a + L a d i a d t + e
T m T load = J d ω d t
J is the total inertia of the drive, which is the sum of the inertia of the motor and the inertia of the working machine. The induced voltage e depends on the speed and magnetic field generated in the excitation coil. It is expressed using (4).
e = c m ω   ;     c m = f ( i f )
The torque of the motor is expressed using (5), and load is expressed using (6).
T m = c m i a
T load = T la + T lb ω + T lc ω 2
where Tla presents friction and load independent of speed (e.g., Coulomb friction); Tlb presents friction and load linearly, depending on the speed (e.g., viscous friction); and Tlc presents the friction and load quadratically, depending on the speed. In the case that the excitation is switched on before the armature circuit, the equation of the excitation circuit (1) does not need to be taken into account, because the transient phenomenon at switching on ends before the armature is switched on. Also, a constant magnetic field is established and the motor constant cm no longer changes. The AC/DC Supply/Control unit presented in Figure 1 responds to the speed and current changes. Its response results in a ua voltage. By considering the measured voltage ua(t) in Equation (2), the AC/DC power supply is covered and its model is not required. In this way, an improved innovative approach is presented.
By entering (4)–(6) into (2) and (3), the system of two coupled deferential equations is obtained, which presents the DC motor’s drive model. It is written in (7) and (8).
u a ( t ) = i a R a + L a d i a d t + c m ω
c m i a T la + T lb ω + T lc ω 2 = J d ω d t
Three different cases can be considered using (7) and (8):
  • Only the DC motor is considered. Inertia J is the inertia of the motor and Tload is the friction of the motor.
  • The DC motor and working machine with no load were considered. Inertia J is the inertia of the motor + inertia of the working machine, and Tload is the friction of the motor + friction of the working machine.
  • The DC motor and working machine with load were considered. Inertia J is the inertia of the motor + the inertia of the working machine. Tload is the friction of the motor + the friction of the working machine + load.

2.2. Numerical Solution of a System of Differential Equations

The system of differential Equations (7) and (8) was solved using the Runge-Kutta method for numerical integration. In order to obtain the relationship between the accuracy and the calculation time depending on the time step of the measurement, we used the Runge–Kutta method of the first order, second order, third order and fourth order.
The first-order method is presented in (9).
K 1 = f x ( k ) , y ( k ) ;   y ( k + 1 ) = y ( k ) + h K 1
y is the searched value (in our case, ia and ω were searched). h is the calculated step, k is a counter, x is an independent variable (in our case, t) and f is the derivative (in our case, dia/dt and /dt).
The second-order method is presented in (10).
K 1 = f x ( k ) , y ( k ) ;   K 2 = f x ( k ) + h , y ( k ) + h K 1 ;   y ( k + 1 ) = y ( k ) + h 1 2 K 1 + 1 2 K 2
The third-order method is presented in (11).
K 1 = f x ( k ) , y ( k ) ;   K 2 = f x ( k ) + 1 2 h , y ( k ) + 1 2 h K 1 ;   K 3 = f x ( k ) + h , y ( k ) h K 1 + 2 h K 2 ;   y ( k + 1 ) = y ( k ) + h 1 6 K 1 + 2 3 K 2 + 1 6 K 3
The fourth-order method is presented in (12).
K 1 = f x ( k ) , y ( k ) ;   K 2 = f x ( k ) + 1 2 h , y ( k ) + 1 2 h K 1 ;   K 3 = f x ( k ) + 1 2 h , y ( k ) + 1 2 h K 2 ;   K 4 = f x ( k ) + h , y ( k ) + h K 3 y ( k + 1 ) = y ( k ) + h 1 6 K 1 + 1 3 K 2 + 1 3 K 3 + 1 6 K 4  
From (7) and (8) the derivatives of current and speed are expressed, and (13) and (14) are written.
d i a d t = 1 L a u a i a R a c m ω = f t , i a , ω
d ω d t = 1 J c m i a T la + T lb ω + T lc ω 2 = g t , i a , ω
Based on (9) and using (13) and (14), the equations for simulation of current and speed using the first-order method are written in (15) and (16).
K 1 = f t ( k ) , i a ( k ) , ω ( k )   = 1 L a u a ( k 1 ) i a ( k 1 ) R a c m ω ( k 1 ) ;     i a ( k + 1 ) = i a ( k ) + h K 1
L 1 = g t ( k ) , i a ( k ) , ω ( k ) = 1 J c m i a ( k ) T la + T lb ω ( k ) + T lc ω ( k ) 2 ;     ω ( k + 1 ) = ω ( k ) + h L 1
The calculation with the first-order method is presented using (15) and (16). Similar expressions can be written for the second-, third- and fourth-order methods. In the case of the second-order method K1, L1, K2, L2, ia and ω were calculated for each time step; in the case of the third-order method K1, L1, K2, L2, K3, L3, ia and ω were calculated for each time step; and. K1, L1, K2, L2, K3, L3, K4, L4, ia and ω were calculated in the case of the fourth-order method. Since we were interested in calculation times, we would like to mention that, for the considered problem, in the case of the first-order method we solved 4 equations in each time step; in the case of the second-order method we solved 6 equations in each time step; in the case of the third-order method we solved 8 equations in each time step; and in the case of the fourth-order method we solved 10 equations in each time step.

2.3. Objective Function

We are dealing with an optimization problem. The simulated data obtained by solving differential equations must be as similar as possible to the measured data. The Objective Function (OF) is defined as the sum of squares of differences between the speed and current responses. It is written in (17).
OF = 1 N i = 1 N i a _ s i m u l a t e d _ i i a _ m e a s u r e d _ i i a _ m e a s u r e d _ max 2 + ω s i m u l a t e d _ i ω m e a s u r e d _ i ω m e a s u r e d _ max 2
N is the number of measured points. The speed and current are normalized and divided with the maximum measured value to obtain the same influence of speed and current on the OF.
The first part of the OF is the comparison of the simulated current with the measured one. The electrical subsystem that describes the flow of the current is presented by (7) (see Section 2.1). It contains the searched parameters Ra, La and cm, and their selection influences the flow of the current. The second part of the OF is the comparison of the simulated speed with the measured one. The mechanical subsystem that describes the flow of speed is described by (8) (see Section 2.1). It contains the searched parameters J, cm, Tla, Tlb and Tlc, and their selection influences the flow of the speed. Based on the described inclusion of the searched parameters in (7) and (8), it can be seen that the course of the current and speed is influenced by changing Ra, La, cm, J, Tla, Tlb and Tlc. The objective function is designed in such a way that the simulated courses of current and speed are as close as possible to the measured ones.

2.4. Measurements

Measurements are an important part of the parameters’ determination process. Well-executed measurements are a condition for successful parameter determination. For the tests made in the scope of the presented work, a laboratory drive was used, presented in Figure 2.
The drive shown in Figure 2 consists of the following parts:
  • Supply/Control unit: SIEMENS AG, (Munich, Germany), SIMOREG DC-Master 6RA7013-6DV62-0-Z.
  • DC motor with separate cooling: SIEMENS AG, (Munich, Germany), 1GG5104-0ED40-6VV1.
  • Pulse encoder: HUBNER (Berlin, Germany), P0G 9D 1024.
  • Motor used for the load simulation: SIEMENS AG, (Munich, Germany), 1LA7139-4AA10-Z FDB0.
  • The “Trace” function was used for the measurements, which is a part of the SIEMENS AG, (Munich, Germany) “Drive Monitor” software, version V05.05.02.00_00.00.02.99, release V05.05.02.00_27.00.00.00.
An example of the measurement is shown in Figure 3.
As was mentioned, only a part of the measurement after the transient phenomena in the field circuit can be used, because, with that, (1) can be eliminated for the parameters’ determination procedure, cm becomes a constant and only (2) and (3) are being solved.
In order to determine the most appropriate measurement, we performed various tests:
  • Different sampling times were used, which were 3.3 ms, 6.6 ms, 16.5 ms and 33 ms.
  • Not only speed up of the drive was considered but also speed down.
  • Different parts of the measured signal were used: full signal, signal without a pre-trigger, speed up only, speed up + same time stationary operation, speed up + stationary operation for twice the start time.
Different sampling times were used to make the article useful to a wider readership. The acceleration times are different for different drives, and the sampling times are limited by the measuring devices used. The results for a different number of measured points in the drive acceleration area are shown by changing the sampling time.

3. Used Evolutionary Methods and Algorithm

As described, we were dealing with an inverse problem, which was solved by a direct approach, for which we needed an optimization method. Some authors used standard optimization [1,3,4,5,6,7,8,9,10,11], and some of them used evolutionary methods, mainly only one of the methods [12,13,14,15,16].
Due to the good properties of evolutionary methods, such as their capacity to avoid local minima, their independence from explicit function derivations and their establishment of stable calculation procedures, we decided to rely solely on these methods. It is well known that different methods are suitable for different problems. Therefore, we tested several different methods on the presented problem, with the aim of determining the most appropriate one for the given and similar problems.
Although we have selected evolutionary methods, it should be noted that their use also has disadvantages. Due to the large number of objective function calculations, especially in the case of a large number of parameters and, thus, a large population, it is a problem if the calculation of the objective function is time-consuming. The quality of the evolutionary method depends on the setting of the parameters of the method, which are not universal but are set according to the problem being solved. The speed of the method depends on the area in which it searches, and only this is defined by the limits of the search parameters. Too large an area is often set, because the exact boundaries of the problem are not known. Finally, the quality of the solution is related directly to a correctly defined objective function, which can be a problem in the case of challenging real-world problems.
The first selected method was DE [14,17,18,19,20,21,22,23,24,25,26]. Some of DE’s advantages are ease of use, simple structure, speed and robustness. In the case of DE, different strategies can be used, generally written as DE/x/y/z. DE indicates Differential Evolution, x represents a string labeling the vector to be perturbed, y is the number of difference vectors considered for perturbation and z stands for the type of crossover being used. Different strategies perform differently for different problems; therefore, we tested three different strategies, which are DE/best/1/exp, DE/rand/1/exp and DE/best/1/bin. For all calculations we set the differential weight F to 0.6 and the crossover probability CR to 0.8.
The second selected algorithm was the more recent ABC [27,28,29,30,31,32,33,34]. It was proposed by D. Karaboga et al. [27]. The ABC algorithm consists of three steps: sending the employed bees onto the food sources and then measuring their nectar; selection of the food sources by the onlookers after sharing the information of the employed bees and determining the nectar amount of the food; determining the scout bees and then sending them onto possible food sources. Since a scout bee might not be employed for every iteration, the number of fitness evaluations (FEs) cannot be determined statistically. The limit value used as the control parameter for the bee population was set to 100, and only one scout bee was used.
The third selected algorithm was TLBO [35,36,37,38,39,40,41,42,43], which is a recently proposed population-based algorithm, which simulates the teaching-learning process in a classroom. TLBO is divided into a Teacher Phase and a Learner Phase. In our work the TLBO implementation was used, where the duplicate elimination phase was omitted. The number of FEs consumed can be determined statistically as: FEs = 2 × population size × iterations. TLBO has no additional control parameters beside the population size.
This work presents the determination of seven parameters of the DC motor drive, which are the motor’s resistance Ra, the motor’s inductance La, the motor’s constant cm, the drive’s inertia J and load expressed with parameters Tla for the constant part of friction and load, Tlb for the linear part of friction and load and Tlc for the square part of friction and load. The population size was set to 10 times the number of parameters, which was 70. The limits of the parameters were set according to the Laboratory drive used, presented in Section 2.3. The parameters’ limits used are presented in Table 1.
In the case of a larger or smaller drive, the parameter limits can be increased or decreased accordingly.
To make a fair comparison of the evolutionary methods, the stopping criteria were the same for all methods, which was 140,000 FEs. Because the population size was 70, it means the number of iterations for each of the methods was the following:
  • DE: 2000 iterations, determined by (FEs/population size)
  • TLBO: 1000 iterations, determined by (FEs/2 × population size) due to the two phases of the TLBO algorithm.
  • ABC: ≤2000 iterations, determined by (FEs/population size + scout bees) due to the fact that the scout bee was not employed for every iteration and dynamic counting must be used.
The algorithm for the parameters’ determination procedure is presented in Figure 4.
The text in Figure 4 “using one of the methods” means that only one of the methods was used in the scope of one calculation procedure.

4. Parameters’ Determination

4.1. Test Measurements

We conducted tests and analyses based on four sets of measured data. Two of them were only for speed up, and these are marked as M1U and M2U. In the additional two test cases, we added speed down after the speed, marked with M1UD and M2UD, which was not found in the literature. With this, we wanted to find out if we could achieve better results by adding the speed-down response. The data used for M1U, M2U, M1UD and M2UD are presented in Table 2.
The meaning of tspeed_up and tspeed_down is presented in Figure 3. The measured signals will be visible when presented together with the results in Figure 5 and Figure 6.
The steady state part of the M1U is also used for the motor constant cm determination. Considering that, dia/dt equals 0, the (18) is obtained at steady-state operation.
u a ( t ) = i a R a + c m ω                     c m = u a ( t ) i a R a ω  
The value cm = 1.356 Vs is used as a known value in the continuation. The values presented in Table 3 can be seen in Figure 5a.

4.2. Results for M1U, M2U, M1UD and M2UD

To determine the parameters, the transient phenomenon of current and speed is necessary, since in (7) (see Section 2.1) La is multiplied by dia/dt, and in Equation (8) (see Section 2.1) J is multiplied by /dt. Acceleration is used by many authors, so we added stopping, in order to determine the possibility of improving the approach.
The sampling time of 6.6 ms was used for all measurements (M1U, M2U, M1UD and M2UD). The Runge–Kutta fourth-order method was used for the simulation used in the scope of the OF determination. For the parameters’ determination DE/best/1/exp, DE/rand/1/exp, DE/best/1/bin, TLBO and ABC were used (presented in Section 3). For each evolutionary method, 50 independent runs were conducted. In Table 4, Table 5, Table 6 and Table 7 and also in the continuation of the paper, B is used for the best value, W is used for the worst value, M is used for the mean value and SD is used for the Standard Deviation.
The values of OF and parameters for M1U are presented in Table 4, for M2U in Table 5, for M1UD in Table 6 and for M2UD in Table 7.
To see the difference between the measured and simulated ia and ω, the measured transient phenomena and simulated transient phenomena obtained using the mean values of the calculated parameters using DE/rand/1/exp are presented in Figure 5 for M1U and M2U and in Figure 6 for M1UD and M2UD.
The calculation time was also an important factor, which should be as short as possible. The average calculation times for 50 independent runs are presented in Table 8.

4.3. Convergence of the Methods

When comparing evolutionary methods, the convergence of the method is important information. It is known that, due to the use of a random generator, each calculation flow is different. Therefore, we looked for the average convergence speed of each method, which we showed with the average value of FEs needed to achieve a certain value of OF. The results are shown in Figure 7, where the average FEs required to achieve OF < 1, OF < 10−1, OF < 10−2 and OF < 10−3 values of OF are shown. For better image transparency, the y scale in Figure 7 is a logarithmic scale.
Analysis of the results in terms of accuracy and in terms of the evolutionary method used are the following:
  • The accuracy of the results was comparable for all test cases and used methods. The deviations in the Ra calculation compared to the known values are as follows: M1U 25%, M2UD 19%, M1UD 19% and M2U 19%. For the La calculations, the deviations stood at M1U 23%, M2U 22%, M1UD 10% and M2UD 15%. Slightly smaller deviations were noticeable in the case of M1UD and M2UD, but we could not characterize them as essential.
  • When comparing the methods, it is noticeable that DE/rand/1/exp and TLBO were the most robust, as they gave exactly the same result for each of the 50 calculations, which was not the case for DE/best/1/exp, DE/best/1/bin and ABC.
  • Comparing the convergence of the methods, it can be observed that DE/best/1/exp and DE/best/1/bin were the fastest. DE/rand/1/exp and TLBO were slower and ABC was the slowest.
  • Based on Table 8, it can be seen that the DE methods for all three strategies were three to four times faster than TLBO and ABC.
Since DE/rand/1/exp and TLBO were the most robust methods and DE/rand/1/exp was significantly faster, the DE/rand/1/exp method will be used for further analysis shown in the article. Also, since the results for all the test cases were comparable, we will use only the M1U test case for further analysis (sampling time, order of integration method, repeatability of measurement and length of measurement).

5. Analysis of the Interaction of the Order of Integration Method and Sampling Time

The power supply used enabled measurements with different sampling times. We used sampling times of 3.3 ms, 6.6 ms, 16.5 ms and 33 ms. To simulate the motor response, which is necessary for the Objective Function’s evaluation, we used the first-order, second-order, third-order and fourth-order Runge–Kutta methods.
The results using the shortest sampling time of 3.3 ms for all for the Runge–Kutta (RK) methods for M1U and DE/rand/1/exp are presented in Table 9.
To see the difference between the measured and simulated ia and ω, the measured and simulated values obtained using the mean values of the calculated parameters using RK first order and fourth order in the case of 3.3 ms sampling time are presented in Figure 8.
The results obtained using a sampling time of 6.6 ms for all the Runge–Kutta (RK) methods for M1U and DE/rand/1/exp are presented in Table 10.
To see the difference between the measured and simulated ia and ω, the measured and simulated values obtained using the mean values of the calculated parameters using the RK first order and fourth order in the case of 6.6 ms sampling time are presented in Figure 9.
The results using the sampling time of 16.5 ms for all four Runge–Kutta (RK) methods for M1U and DE/rand/1/exp are presented in Table 11.
To see the difference between the measured and simulated ia and ω, the measured and simulated values obtained using the mean values of the calculated parameters using the RK first order and fourth order in the case of a 16.5 ms sampling time are presented in Figure 10.
The results using the sampling time of 33 ms for all four Runge–Kutta (RK) methods for M1U and DE/rand/1/exp are presented in Table 12.
To see the difference between the measured and simulated ia and ω, the measured and simulated values obtained using the mean values of the calculated parameters using the RK first order and fourth order in the case of a 33 ms sampling time are presented in Figure 11.
The analyses of the results presented in Table 9, Table 10, Table 11 and Table 12 and Figure 8, Figure 9, Figure 10 and Figure 11 are the following:
  • For the sampling time 3.3 ms, the OF for the first-order method was 5.7 ∙ 10−4, and all the other methods were in the range from 4.8 to 5 ∙ 10−4. A small difference of the simulated values can be seen in Figure 8, where the simulated values for the first- and fourth-order methods are presented. All the calculated parameters were correct.
  • For the sampling time 6.6 ms, the OF for the first-order method was 5.8 ∙ 10−4, and all the other methods were in the range from 3.4 to 3.8 ∙ 10−4. The parameters calculated with the first-order method showed a larger deviation than with the other methods. A difference of the simulated values can be seen in Figure 9, where the simulated values for the first- and fourth-order methods are presented.
  • For the sampling time 16.5 ms, the OF for the first-order method was 1.5 ∙ 10−3, and all the other methods were in the range from 5.4 to 5.8 ∙ 10−4. In the case of the first-order method, the calculated parameters were not correct. In the case of the second-order method, Ra and La were not accurate. In the case of the third- and fourth-order methods, La was not accurate. In Figure 10 it can be seen that only two measuring points were present at the speed up of the motor. Nevertheless, we can use the second-, third- and fourth-order methods to estimate the parameters. A big difference of the simulated values can be seen in Figure 10, where the simulated values are presented for the first-and fourth-order methods.
  • For the sampling time 33 ms, the OF for the first-order method was 4.5 ∙ 10−3, and all the other methods were in the range from 1.1 to 1.9 ∙ 10−3. The calculated values of the parameters were not calculated well enough with respect to the known values.
Based on the given analysis, we can conclude that the first-order method was not appropriate, except in the case of a small sampling time. The second-, third- and fourth-order methods were comparable. It is advantageous if there are at least a few measurement points during the acceleration, such as in Figure 8 and Figure 9.

6. Repeatability of Parameters’ Determination

It is important that the determination of parameters is a stable and repeatable process. In Section 4.2 we defined a robust evolutionary method, which is DE/rand/1/exp. The question was whether we would obtain comparable results if we repeated the measurement under the same conditions. Therefore, we repeated the measurement under the same conditions 10 times and looked at how the calculated parameters were comparable to each other.
The test was made using the M1U data, a sampling time of 6.6 ms, DE/rand/1/exp and the fourth-order method. The results for the 1st to 5th measurements are presented in Table 13 and for the 6th to 10th measurements in Table 14.
The overall mean values of the results for all 10 measurements were calculated based on the results presented in Table 13 and Table 14.
The deviations of OFs and parameters from the overall mean values for all 10 measurements are presented in Table 15.
From Table 15 it can be seen that the biggest deviation was 6% in case of the La calculation. The deviations of all the other parameters’ calculations were less than 2%. Based on this, we can conclude that the procedure ensures good repeatability of parameter determination.
The statistical distribution of the results is presented in Figure 12, where histograms are presented for OF, Ra, La, cm, J and Tlb.
In the histograms presented in Figure 12, blocks of 50 values can be seen clearly. They appear due to the fact that DE/rand/1/exp is very robust, so for each run with unchanged input data, the results are the same. For each set of measured input data, 50 independent runs were conducted. In Figure 12, each histogram/chart represents 500 values divided into 10 blocks of 50 values, since 10 sets of measured data were used in the calculation.

7. Analysis of the Influence of Measurement Length

The different lengths of the measurement influenced the calculation time. The tests were made for the following measured signals with a sampling time of 6.6 ms using the M1U data:
  • The whole measurement, which contained a pre-trigger, speed up and continuous operation. The whole signal consisted of 400 measurement points.
  • A signal with no pre-trigger that contained speed up and continuous operation. The length of the signal with no pre-trigger was 362 measured points.
  • Only speed up, which contained 56 measured points with a time duration of 0.35 s. The length of the signal was 56 measured points.
  • Speed up with a time duration of 0.35 s and continuous operation with the same time duration of 0.35 s. The length of the signal was 109 measured points.
  • Speed up with a time duration of 0.35 s and continuous operation with a time duration of 0.70 s. The length of the signal was 162 measured points.
The results are presented in Table 16.
If we compare the results for all five measured signals by calculating the overall mean values in the same way as presented in Figure 12, we obtain the results presented in Table 17.
Based on Table 17 it can be seen that the deviations of Ra, La, cm and J were less than 2%. A bigger deviation can be observed in the case of Tlb, which appeared for the START ONLY measured signal, and it was up to 32%. Based on this, it is suggested to use speed up and at least the same time of the continuous operation. Using signal 0.35 start + 0.35 continuous operation, the calculation times decreased from approximately 9 s to approximately 7 s. The measurements and simulations for the whole signal, signal with no pre-trigger, start only signal and signal containing start + 0.35 s of continuous operation are shown in Figure 13.

8. Validation Based on the Simulated Input Data

To validate the presented approach, the DC motor model was used to generate ua(t), ia(t) and ω. Generation was made using the RK fourth-order method and a time step of 6.6 ms. The parameters used for the generation of the simulated input data were the known data from the presented drive, shown in Table 18 in the column “Parameters used for the generation of the responses.” For the input voltage, a linear increase up to 255 V was used, which can be seen in Figure 14.
The generated responses were then used as input data for the parameters’ determination. For the parameters’ determination the RK fourth-order method was used (the same as was used for the simulated data generation). The results of the 50 independent runs are presented in Table 18 in the column “Calculated parameters”.
From Table 18 it can be seen that the calculated parameters are the same as the parameters used for simulated data determination. In this way, the used model was validated against the used data set.
The input data and simulated data obtained as a result of parameters’ determination are presented in Figure 14.

9. Conclusions

This paper introduced an approach that incorporates the measurement of voltage and not only current and speed response, as is common in the literature. It was proven that such an approach is appropriate to be used. Not only the parameters of drives with small DC motors but also drives with big DC motors powered by control/supply units can be determined with such an approach. The good side of this approach is simplicity, because the control/supply unit does not need to be modeled; it is considered with measured voltage.
The presented use of measured values are not only for speed up but also for speed down. It was found that the use of speed down did not improve the parameter determination process significantly.
An important part of the presented work is the determination of the appropriate evolutionary optimization method. The algorithms DE/best/1/exp, DE/rand/1/exp, DE/best/1/bin, TLBO and ABC were used and analyzed. The most robust were DE/rand/1/exp and TLBO, which led to the same calculated parameters for each of the 50 calculations in the case of all 4 test data (M1U, M2U, M1UD, M2UD). The convergence of TLBO was slightly better (Figure 7), but its calculation times were 3–6 times longer that the calculation times of DE/rand/1/exp (Table 8). DE/rand/1/exp can be presented as the most appropriate method for the presented problem. The selected method is the best only between those tested. We did not test other methods, since DE/rand/1/exp satisfied our expectations considering robustness and speed. Due to the NFL theorem, there are a plethora of different evolutionary algorithms and their derivatives, making it impossible to test all of them. Also, most of the algorithms were tested on benchmark problems, but their performance is usually very different on real-world problems.
The comparison of results using four sampling times (3.3 ms, 6.6 ms, 16.5 ms and 33 ms) and four different methods used for solving of coupled differential equations were analyzed (first-order, second-order, third-order and fourth-order methods). It can be concluded that at least some measured points should be at the speed-up interval, such as in the case of 3.3 ms (Figure 8) and 6.6 ms (Figure 9). The sampling times of 16.5 ms and 33 ms were not appropriate in the case of the used drive (Figure 10 and Figure 11). The results obtained using the first-order method were the worst. The results obtained using the second-, third- and fourth-order methods were better and comparable.
The parameters’ determination procedure showed very good repeatability, although three measured signals (speed, current and voltage) were used as the input data. In Section 6 it can be seen that, for 10 repeated measurements, the results were almost identical.
If we want to speed up the parameter determination process, we can use only part of the signal. The analysis in Section 7 showed that it is sufficient to use a signal consisting of a part of the signal when the motor/drive is accelerating and a part of the signal when the motor/drive is in constant operation, which is approximately the same length as the acceleration time.
Future research will involve exploring evolutionary methods and utilizing parameter determination techniques, not only for this specific motor type but also for various other types of motors. As a future work, also the findings of this research could be analyzed using artificial intelligence.

Author Contributions

Conceptualization, M.J.; methodology, M.J. and M.T.; software, M.J. and M.R.; validation, M.T., M.R. and M.J.; formal analysis, M.R. and M.J.; investigation, M.T. and M.J.; writing—original draft preparation, M.J.; writing—review and editing, M.T. and M.R.; supervision, M.J.; project administration, M.J.; funding acquisition, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research and Innovation Agency, grant number P2-0114.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the use of the laboratory equipment.

Conflicts of Interest

The authors declare no conflicts.

References

  1. Hadef, M.; Mekideche, M.R. Parameter identification of a separately excited dc motor via inverse problem methodology. Turk. J. Electr. Eng. Comp. Sci. 2009, 17, 99–106. [Google Scholar] [CrossRef]
  2. Wu, W. DC Motor Parameters Identification Using Speed Step Response. Model. Simul. Eng. 2012, 2012, 189757. [Google Scholar] [CrossRef]
  3. Adewusi, S. Modeling and Parameters Identification of a DC Motor Using Constraint Optimization Technique. IOSR J. Mech. Civ. Eng. 2016, 13, 46–56. [Google Scholar]
  4. Gao, D.; Wu, S.; Yu, J.; Wang, M.; Wang, Y. Parameter identification of DC motor model based on improved dynamic forgetting factor recursive least squares method. In Proceedings of the 2022 IEEE 8th International Conference on Smart Instrumentation, Measurement and Applications (ICSIMA), Melaka, Malaysia, 26–28 September 2022; pp. 282–286. [Google Scholar]
  5. Ivanov, D.V.; Sandler, I.L.; Chertykovtseva, N.V.; Mitroshin, D.I.; Ivanova, O.S.; Kormakov, A.A. Identification of Parameters of DC Motor of Independent Excitation by Noisy Data. In Proceedings of the 2021 3rd International Conference on Control Systems, Mathematical Modeling, Automation and Energy Efficiency (SUMMA), Lipetsk, Russia, 10–12 November 2021; pp. 194–198. [Google Scholar]
  6. Li, M.; Ma, Y. Parameter Identification of DC Motor based on Compound Least Square Method. In Proceedings of the 2020 IEEE 5th Information Technology and Mechatronics Engineering Conference (ITOEC), Chongqing, China, 12–14 June 2020; pp. 1107–1111. [Google Scholar]
  7. Awoda, M.L.; Ramzy, S.A. Parameter Estimation of a Permanent Magnets DC motor. Iraqi J. Electr. Electron. Eng. 2019, 15, 28–36. [Google Scholar] [CrossRef]
  8. Arshad, S.; Qamar, S.; Jabbar, T.; Malik, A. Parameter Estimation of a DC Motor Using Ordinary Least Squares and Recursive Least Squares Algorithms. In Proceedings of the 8th International Conference on Frontiers of Information Technology, Islamabad, Pakistan, 12–13 December 2010. [Google Scholar]
  9. Hadef, M.; Bourouina, A.; Mekideche, M. Parameter identification of a DC motor via moments method. Iran. J. Electr. Comput. Eng. 2008, 7, 159–163. [Google Scholar]
  10. Amiri, M.S.; Ibrahim, M.F.; Ramli, R. Optimal parameter estimation for a DC motor using genetic algorithm. Int. J. Power Electron. Drive Syst. 2020, 11, 1047–1054. [Google Scholar] [CrossRef]
  11. Dupuis, A.; Ghribi, M.; Kaddouri, A. Multiobjective genetic estimation of DC motor parameters and load torque. In Proceedings of the 2004 IEEE International Conference on Industrial Technology, 2004. IEEE ICIT ’04., Hammamet, Tunisia, 8–10 December 2004; Volume 3, pp. 1511–1514. [Google Scholar]
  12. Rodríguez-Molina, A.; Villarreal-Cervantes, M.G.; Aldape-Pérez, M. Optimal Adaptive Control of a DC Motor Using Differential Evolution Variants. In Proceedings of the 2017 IEEE 7th Annual International Conference on CYBER Technology in Automation, Control, and Intelligent Systems (CYBER), Honolulu, HI, USA, 31 July–4 August 2017; pp. 283–288. [Google Scholar]
  13. Puangdownreong, D.; Hlungnamtip, S.; Thammarat, C.; Nawikavatan, A. Application of flower pollination algorithm to parameter identification of DC motor model. In Proceedings of the IEEE International Electrical Engineering Congress (IEECON), Pattaya, Thailand, 8–10 March 2017; pp. 1–4. [Google Scholar]
  14. Hafez, I.; Dhaouadi, R. Parameter Identification of DC Motor Drive Systems using Particle Swarm Optimization. In Proceedings of the 2021 International Conference on Engineering and Emerging Technologies (ICEET), Istanbul, Turkey, 27–28 October 2021; pp. 1–6. [Google Scholar]
  15. Veček, N.; Mernik, M.; Črepinšek, M. A chess rating system for evolutionary algorithms: A new method for the comparison and ranking of evolutionary algorithms. Inf. Sci. 2014, 277, 656–679. [Google Scholar] [CrossRef]
  16. Das, S.; Suganthan, P.N. Differential Evolution: A Survay of the State-of-the-Art. IEEE Trans. Evol. Comput. 2011, 15, 4–31. [Google Scholar] [CrossRef]
  17. Mokan, M.; Sharma, K.; Sharma, H.; Verma, C. Gbest guided differential evolution. In Proceedings of the 9th International Conference on Industrial and Information Systems, Gwalior, India, 15–17 December 2014; pp. 1–6. [Google Scholar]
  18. He, R.J.; Yang, Z.Y. Differential evolution with adaptive mutation and parameter control using Levy probability distribution. J. Comput. Sci. Technol. 2012, 27, 1035–1055. [Google Scholar] [CrossRef]
  19. Mohamed, A.W.; Sabry, H.Z.; Elaziz, T.A. Real parameter optimization by an effective differential evolution algorithm. Egipt. Inform. J. 2013, 14, 27–53. [Google Scholar] [CrossRef]
  20. Reed, H.M.; Nichols, J.M.; Earls, C.J. A modified differential evolution algorithm for damage identification in submerged shell structures. Mech. Syst. Signal Process. 2013, 39, 396–408. [Google Scholar] [CrossRef]
  21. Chattopadhyay, S.; Sanyal, S.K. Optimization of Control Parameters of Differential Evolution Technique for the Design of FIR Pulse-shaping Filter in QPSK Modulated System. J. Common. 2011, 6, 558–570. [Google Scholar] [CrossRef]
  22. Saruhan, H. Differential evolution and simulated annealing algorithms for mechanical systems design. Eng. Sci. Technol. Int. J. 2014, 17, 131–136. [Google Scholar] [CrossRef]
  23. Mallipeddi, R.; Suganthan, P.N.; Pan, Q.K.; Tasgeriren, M.F. Differential evolution algorithm with ensemble of parameters and mutation strategies. Appl. Soft Comput. 2011, 11, 1679–1696. [Google Scholar] [CrossRef]
  24. Rocca, P.; Oliveri, G.; Massa, A. Differential Evolution as Applied to Electromagnetics. IEEE Trans. Antennas Propag. 2011, 50, 38–49. [Google Scholar] [CrossRef]
  25. Karaboga, D.; Basturk, B. On the performance of artificial bee colony (ABC) algorithm. Appl. Soft Comput. 2008, 8, 687–697. [Google Scholar] [CrossRef]
  26. Mernik, M.; Liu, S.H.; Karaboga, D.; Črepinšek, M. On clarifying misconceptions when comparing variants of the Artificial Bee Colony Algorithm by offering a new implementation. Inf. Sci. 2015, 291, 115–127. [Google Scholar] [CrossRef]
  27. Kiran, M.S.; Gündüz, M. The Analysis of Peculiar Control Parameters of Artificial Bee Colony Algorithm on the Numerical Optimization Problems. Int. J. Comput. Commn. 2014, 2, 127–136. [Google Scholar] [CrossRef]
  28. Özyon, S.; Aydin, D. Incremental artificial bee colony with local search to economic dispatch problems with ramp rate limits and prohibited operating zones. Energ. Convers. Manag. 2013, 65, 397–407. [Google Scholar] [CrossRef]
  29. Dwivedl, A.K.; Ghosh, S.; Londhe, N.D. Modified artificial bee colony optimisation based FIR filter design with experimental validation using field-programmable gate array. IET Signal Process. 2016, 10, 955–964. [Google Scholar] [CrossRef]
  30. Jing, B.; Hong, L. Improved Artificial Bee Colony Algorithm and Application in Path Planning of Crowd Animation. Int. J. Control Autom. 2015, 8, 53–66. [Google Scholar] [CrossRef]
  31. Yan, G.; Li, C. An Effective Refinement Artificial Bee Colony Optimization Algorithm Based on Chaotic Search and Application for PID Control Tuning. J. Comput. Inf. Syst. 2011, 7, 3309–3316. [Google Scholar]
  32. Ozturk, C.; Karaboga, D. A novel clustering approach: Artificial Bee Colony (ABC) algorithm. Appl. Soft Comput. 2011, 11, 652–657. [Google Scholar]
  33. Karaboga, B.; Akay, B. A Comparative Study of Artificial Bee Colony Algorithm. Appl. Math. Comput. 2009, 214, 108–132. [Google Scholar] [CrossRef]
  34. Karaboga, D.; Basturk, B. A powerful and efficient algorithm for numerical function optimization: Artificial bee colony (ABC) algorithm. J. Glob. Optim. 2007, 39, 459–471. [Google Scholar] [CrossRef]
  35. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching-Learning-Based Optimization: An optimization method for continuous non-linear large scale problems. Inf. Sci. 2012, 183, 1–15. [Google Scholar] [CrossRef]
  36. Črepinšek, M.; Liu, S.H.; Mernik, L. A note on leaching-learning-based optimization algorithm. Int. Sci. 2012, 212, 79–93. [Google Scholar]
  37. Črepinšek, M.; Liu, S.H.; Mernik, L.; Mernik, M. Is a comparison of results meaningful from the inexact replications of computational experiments. Soft Comp. 2016, 20, 223–235. [Google Scholar] [CrossRef]
  38. Sahu, B.K.; Pati, S.; Mohanty, P.K.; Panda, S. Teaching-learning based optimization algorithm based fuzzy-PID controller for automatic generation control of multi-area power system. Appl. Soft Compt. 2015, 27, 240–249. [Google Scholar] [CrossRef]
  39. Pickard, J.K.; Carreter, J.A.; Bhavsar, V.C. On the convergence and original bias of the Teaching-Learning-Based-Optimization algorithm. Appl. Soft Comp. 2016, 46, 115–127. [Google Scholar] [CrossRef]
  40. Baghlani, A.; Makiabadi, M.H. Teaching-learning based optimization algorithm for shape and size optimization of truss structures with dynamic frequency constraints. IJST Trans. Civ. Eng. 2013, 37, 409–421. [Google Scholar]
  41. Waghmare, G. Comments on A note on teachnig-learning-based optimization algorithm. Inf. Sci. 2013, 229, 159–169. [Google Scholar] [CrossRef]
  42. Rao, R.V.; Patel, V. An elitist teaching-learning-based optimization algorithm for solving complex constrained optimization problems. Int. J. Ind. Eng. Comput. 2012, 3, 535–560. [Google Scholar] [CrossRef]
  43. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching-learning-based optimization: A novel method for constrained mechanical design optimization problems. Comput. Aided Des. 2011, 43, 303–315. [Google Scholar] [CrossRef]
Figure 1. The DC motor’s drive.
Figure 1. The DC motor’s drive.
Mathematics 12 00042 g001
Figure 2. Laboratory drive used for the tests.
Figure 2. Laboratory drive used for the tests.
Mathematics 12 00042 g002
Figure 3. An example of the measurement.
Figure 3. An example of the measurement.
Mathematics 12 00042 g003
Figure 4. Algorithm of the parameters’ determination procedure: (a) algorithm; (b) presentation of the input and calculated data used for the OF calculation.
Figure 4. Algorithm of the parameters’ determination procedure: (a) algorithm; (b) presentation of the input and calculated data used for the OF calculation.
Mathematics 12 00042 g004
Figure 5. Measured and calculated transient phenomena for (a) M1U and (b) M2U.
Figure 5. Measured and calculated transient phenomena for (a) M1U and (b) M2U.
Mathematics 12 00042 g005
Figure 6. Measured and calculated transient phenomena for (a) M1UD and (b) M2UD.
Figure 6. Measured and calculated transient phenomena for (a) M1UD and (b) M2UD.
Mathematics 12 00042 g006
Figure 7. Convergence of the methods in case of (a) M1U, (b) M2U, (c) M1UD and (d) M2UD.
Figure 7. Convergence of the methods in case of (a) M1U, (b) M2U, (c) M1UD and (d) M2UD.
Mathematics 12 00042 g007
Figure 8. Measured and calculated responses for M1U and sampling time 3.3 ms used for the simulation (a) RK first-order and (b) RK fourth-order methods.
Figure 8. Measured and calculated responses for M1U and sampling time 3.3 ms used for the simulation (a) RK first-order and (b) RK fourth-order methods.
Mathematics 12 00042 g008
Figure 9. Measured and calculated responses for M1U and sampling time of 6.6 ms using for the simulation (a) RK first-order and (b) RK fourth-order methods.
Figure 9. Measured and calculated responses for M1U and sampling time of 6.6 ms using for the simulation (a) RK first-order and (b) RK fourth-order methods.
Mathematics 12 00042 g009
Figure 10. Measured and calculated responses for M1U and sampling time 16.5 ms using for the simulation (a) RK first-order and (b) RK fourth-order methods.
Figure 10. Measured and calculated responses for M1U and sampling time 16.5 ms using for the simulation (a) RK first-order and (b) RK fourth-order methods.
Mathematics 12 00042 g010
Figure 11. Measured and calculated responses for M1U and sampling time 33 ms used for the simulation (a) RK first-order and (b) RK fourth-order methods.
Figure 11. Measured and calculated responses for M1U and sampling time 33 ms used for the simulation (a) RK first-order and (b) RK fourth-order methods.
Mathematics 12 00042 g011
Figure 12. Histograms presenting the statistical distribution of the results.
Figure 12. Histograms presenting the statistical distribution of the results.
Mathematics 12 00042 g012
Figure 13. Measured and simulated values for (a) the whole signal, (b) a signal with no pre-trigger, (c) a start only signal, (d) a start + 0.35 s of continuous operation signal.
Figure 13. Measured and simulated values for (a) the whole signal, (b) a signal with no pre-trigger, (c) a start only signal, (d) a start + 0.35 s of continuous operation signal.
Mathematics 12 00042 g013
Figure 14. Input and simulated signals (input and simulated data of ia being covered and also input and simulated data of ω being covered).
Figure 14. Input and simulated signals (input and simulated data of ia being covered and also input and simulated data of ω being covered).
Mathematics 12 00042 g014
Table 1. Parameters’ limits.
Table 1. Parameters’ limits.
ParameterLower LimitUpper Limit
Motor’s resistance Ra (Ω)0100
Motor’s inductance La (H)0100
Motor’s constant cm (Vs)05
Drive’s inertia J (kg∙m2)01
Constant part of friction and load Tla (Nm)020
Linear part of friction and load Tlb (Nm∙s)09.55 ∙ 10−2 (20 Nm torque at speed 209 s−1; 182 s−1 is motor-rated speed)
Square part of friction and load Tlc (Nm∙s2)04.56 ∙ 10−4 (20 Nm torque at speed 209 s−1; 182 s−1 is motor-rated speed)
Table 2. The values used to obtain the measured data.
Table 2. The values used to obtain the measured data.
Measured Input Data
ValueM1UM2UM1UDM2UD
tspeed_up (s)00.500.5
tspeed_down (s)//00.5
ωfinal (s−1)182182182182
ia_limit (A)12.48
(120% of Ia_rated)
12.48
(120% of Ia_rated)
12.48
(120% of Ia_rated)
12.48
(120% of Ia_rated)
Load (Nm)no loadno loadno loadno load
Number of measured points400400400400
Table 3. Determination of the motor constant cm based on the measured data MU1.
Table 3. Determination of the motor constant cm based on the measured data MU1.
t (s)Measured
ua (V)
Measured
ia (A)
Measured
ω (s−1)
Known
Ra (Ω)
Calculated
cm (Vs)
0.5254.760.605181.665.661.384
1255.080.780182.815.661.371
1.5251.640.636181.625.661.366
2246.260.674181.205.661.338
Average cm determined, based on 286 measured points for steady state operation from 0.5 s1.356
Table 4. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M1U with a 6.6 ms sampling time.
Table 4. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M1U with a 6.6 ms sampling time.
OF and
Parameters
Method
Known ValueDE/best/1/expDE/rand/1/expDE/best/1/binTLBOABC
OFB-3.8094 · 10−43.8094 · 10−43.8094 · 10−43.8094 · 10−43.8127 · 10−4
W-3.9247 · 10−43.8094 · 10−43.9247 · 10−43.8094 · 10−43.8094 · 10−4
M-3.8164 · 10−43.8094 · 10−43.8141 · 10−43.8094 · 10−43.8097 · 10−4
SD-2.7376 · 10−65.4210 · 10−202.2589 · 10−65.4210 · 10−206.7087 · 10−8
Ra (Ω)M5.664.2144.2134.2144.2134.213
La (H)M0.04720.05830.05830.05830.05830.0583
cm (Vs)M1.3561.3601.3601.3601.3601.360
J (kgm2)M≈3.725 · 10−22.578 · 10−22.576 · 10−22.577 · 10−22.576 · 10−22.576 · 10−2
Tla (Nm)M≈00.03.764 · 10−160.05.684 · 10−160.0
Tlb (Nms)M≈4.8 · 10−34.530 · 10−34.819 · 10−34.626 · 10−34.819 · 10−34.811 · 10−3
Tlc (Nms2)M≈06.244 · 10−66.871 · 10−201.052 · 10−51.984 · 10−194.0320 · 10−8
Table 5. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M2U with a 6.6 ms sampling time.
Table 5. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M2U with a 6.6 ms sampling time.
OF and
Parameters
Method
Known ValueDE/best/1/expDE/rand/1/expDE/best/1/binTLBOABC
OFB-3.3639 · 10−43.3639 · 10−43.3639 · 10−43.3639 · 10−43.4054 · 10−4
W-3.5599 · 10−43.3639 · 10−43.5600 · 10−43.3639 · 10−43.3640 · 10−4
M-3.3757 · 10−43.3639 · 10−43.3757 · 10−43.3639 · 10−43.3685 · 10−4
SD-4.6540 · 10−65.4210 · 10−204.6540 · 10−65.4210 · 10−209.0725 · 10−7
Ra (Ω)M5.664.5304.5304.5304.5304.533
La (H)M0.04720.05770.05770.05770.05770.0578
cm (Vs)M1.3561.3581.3581.3581.3581.358
J (kgm2)M≈3.725 · 10−22.534 · 10−22.531 · 10−22.534 · 10−22.531 · 10−22.531 · 10−2
Tla (Nm)M≈01.466 · 10−45.309 · 10−141.466 · 10−41.385 · 10−136.868 · 10−3
Tlb (Nms)M≈4.8 · 10−34.432 · 10−34.715 · 10−34.432 · 10−34.715 · 10−34.6025 · 10−3
Tlc (Nms2)M≈01.539 · 10−66.710 · 10−201.539 · 10−61.911 · 10−193.9619 · 10−7
Table 6. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M1UD with a 6.6 ms sampling time.
Table 6. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M1UD with a 6.6 ms sampling time.
OF and
Parameters
Method
Known ValueDE/best/1/expDE/rand/1/expDE/best/1/binTLBOABC
OFB-8.0983 · 10−48.0983 · 10−48.0983 · 10−48.0983 · 10−48.0983 · 10−4
W-8.6327 · 10−48.0983 · 10−48.6327 · 10−48.0983 · 10−48.0983 · 10−4
M-8.1090 · 10−48.0983 · 10−48.1304 · 10−48.0983 · 10−48.0983 · 10−4
SD-7.4820 · 10−63.2526 · 10−191.2691 · 10−53.2526 · 10−193.2526 · 10−19
Ra (Ω)M5.664.5514.5514.5514.5514.551
La (H)M0.04720.05220.05220.05220.05220.0522
cm (Vs)M1.3561.3541.3541.3541.3541.354
J (kgm2)M≈3.725 · 10−22.561 · 10−22.561 · 10−22.561 · 10−22.561 · 10−22.561 · 10−2
Tla (Nm)M≈02.093 · 10−262.383 · 10−162.771 · 10−173.553 · 10−161.214 · 10−15
Tlb (Nms)M≈4.8 · 10−39.509 · 10−53.320 · 10−182.853 · 10−42.392 · 10−170.0
Tlc (Nms2)M≈02.775 · 10−52.831 · 10−52.661 · 10−52.831 · 10−52.831 · 10−5
Table 7. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M2UD with a 6.6 ms sampling time.
Table 7. OF and Mean values of the calculated parameters for 50 independent runs using the Runge–Kutta fourth-order method for M2UD with a 6.6 ms sampling time.
OF and
Parameters
Method
Known ValueDE/best/1/expDE/rand/1/expDE/best/1/binTLBOABC
OFB-4.9342 · 10−44.9432 · 10−44.9432 · 10−44.9432 · 10−44.9432 · 10−4
W-6.0877 · 10−44.9432 · 10−46.0877 · 10−44.9432 · 10−44.9432 · 10−4
M-4.9890 · 10−44.9432 · 10−45.0119 · 10−44.9432 · 10−44.9432 · 10−4
SD-2.2426 · 10−51.0842 · 10−192.7178 · 10−51.0842 · 10−196.8147 · 10−13
Ra (Ω)M5.664.5364.5354.5364.5354.535
La (H)M0.04720.05410.05410.05410.05410.0541
cm (Vs)M1.3561.3591.3591.3591.3591.359
J (kgm2)M≈3.725 · 10−22.490 · 10−22.490 · 10−22.490 · 10−22.490 · 10−22.490 · 10−2
Tla (Nm)M≈02.856 · 10−32.911 · 10−154.284 · 10−39.344 · 10−150.0
Tlb (Nms)M≈4.8 · 10−34.979 · 10−35.187 · 10−34.876 · 10−35.187 · 10−35.187 · 10−3
Tlc (Nms2)M≈01.097 · 10−65.645 · 10−201.645 · 10−61.612 · 10−190.0
Table 8. Mean values of calculation times for 50 independent runs using the RK fourth-order method for M1U, M2U, M1UD and M2UD input data.
Table 8. Mean values of calculation times for 50 independent runs using the RK fourth-order method for M1U, M2U, M1UD and M2UD input data.
Measured Data Method
DE/best/1/expDE/rand/1/expDE/best/1/binTLBOABC
M1UM t(s)9.6539.2128.73933.11136.121
M2UM t(s)9.1669.0959.30133.36735.734
M1UDM t(s)9.0999.0369.32233.67736.738
M2UDM t(s)9.0359.6638.67933.91736.603
Table 9. OF and Mean value of the calculated parameters for 50 independent runs: data M1U, sampling time 3.3 ms, simulation RK first, second, third and fourth order, solving method DE/rand/1/exp.
Table 9. OF and Mean value of the calculated parameters for 50 independent runs: data M1U, sampling time 3.3 ms, simulation RK first, second, third and fourth order, solving method DE/rand/1/exp.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueRK First OrderRK Secons OrderRK Third OrderRK Fourth Order
OFB-5.7374 · 10−44.8270 · 10−44.9176 · 10−45.0743 · 10−3
W-5.7374 · 10−44.8270 · 10−44.9176 · 10−45.0743 · 10−4
M-5.7374 · 10−44.8270 · 10−44.9176 · 10−45.0743 · 10−4
SD-0.05.4210 · 10−201.0842 · 10−191.0842 · 10−19
Ra (Ω)M5.664.3184.1084.1444.160
La (H)M0.04720.05640.05150.05360.0549
cm (Vs)M1.3561.3621.3641.3641.364
J (kgm2)M≈3.725 · 10−22.542 · 10−22.560 · 10−22.557 · 10−22.555 · 10−2
Tla (Nm)M≈04.114 · 10−164.571 · 10−164.374 · 10−163.113 · 10−16
Tlb (Nms)M≈4.8 · 10−35.107 · 10−35.087 · 10−35.099 · 10−35.109 · 10−3
Tlc (Nms2)M≈05.226 · 10−207.082 · 10−206.399 · 10−59.638 · 10−20
t (s)M-6.4917.1768.8768.785
Table 10. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 6.6 ms, simulation RK first, second, third and fourth order, solving method DE/rand/1/exp.
Table 10. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 6.6 ms, simulation RK first, second, third and fourth order, solving method DE/rand/1/exp.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueRK First OrderRK Second OrderRK Third OrderRK Fourth Order
OFB-5.8057 · 10−43.4525 · 10−43.5215 · 10−43.8094 · 10−4
W-5.8057 · 10−43.4525 · 10−43.5215 · 10−43.8094 · 10−4
M-5.8057 · 10−43.4525 · 10−43.5215 · 10−43.8094 · 10−4
SD-1.0842 · 10−195.4210 · 10−205.4210 · 10−205.4210 · 10−20
Ra (Ω)M5.664.5464.0674.1764.213
La (H)M0.04720.06430.05040.05480.0583
cm (Vs)M1.3561.3571.3601.3601.360
J (kgm2)M≈3.725 · 10−22.548 · 10−22.604 · 10−22.583 · 10−22.576 · 10−2
Tla (Nm)M≈06.767 · 10−163.563 · 10−53.939 · 10−163.764 · 10−16
Tlb (Nms)M≈4.8 · 10−34.824 · 10−34.798 · 10−34.807 · 10−34.819 · 10−3
Tlc (Nms2)M≈04.580 · 10−204.712 · 10−204.271 · 10−196.871 · 10−20
t (s)M-7.4537.0858.1399.212
Table 11. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 16.5 ms, simulation RK first, second, third and fourth orders, solving method DE/rand/1/exp.
Table 11. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 16.5 ms, simulation RK first, second, third and fourth orders, solving method DE/rand/1/exp.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueRK First OrderRK Second OrderRK Third OrderRK Fourth Order
OFB-1.5423 · 10−35.4764 · 10−45.4974 · 10−45.7259 · 10−4
W-1.5423 · 10−35.4764 · 10−45.4974 · 10−45.7259 · 10−4
M-1.5423 · 10−35.4764 · 10−45.4974 · 10−45.7259 · 10−4
SD-1.5423 · 10−190.01.0842 · 10−191.0842 · 10−19
Ra (Ω)M5.665.4523.7604.1274.237
La (H)M0.04720.12930.06670.07590.0883
cm (Vs)M1.3561.3681.3761.3751.374
J (kgm2)M≈3.725 · 10−22.239 · 10−22.697 · 10−22.551 · 10−22.522 · 10−2
Tla (Nm)M≈01.897 · 10−11.883 · 10−22.216 · 10−25.566 · 10−15
Tlb (Nms)M≈4.8 · 10−37.837 · 10−174.927 · 10−34.924 · 10−35.064 · 10−3
Tlc (Nms2)M≈02.212 · 10−59.463 · 10−199.258 · 10−196.246 · 10−18
t (s)M-8.7037.3137.7999.196
Table 12. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 33 ms, simulation RK first, second, third and fourth orders, solving method DE/rand/1/exp.
Table 12. OF and Mean values of the calculated parameters for 50 independent runs: data M1U, sampling time 33 ms, simulation RK first, second, third and fourth orders, solving method DE/rand/1/exp.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueRK First OrderRK Second OrderRK Third OrderRK Fourth Order
OFB-4.5163 · 10−31.1598 · 10−31.9605 · 10−31.1697 · 10−3
W-4.5163 · 10−31.1598 · 10−31.2945 · 10−31.1697 · 10−3
M-4.5163 · 10−31.1598 · 10−31.9160 · 10−31.1697 · 10−3
SD-8.6736 · 10−190.01.2690 · 10−42.1684 · 10−19
Ra (Ω)M5.668.0493.3204.0765.011
La (H)M0.04720.29360.06640.1090.0764
cm (Vs)M1.3561.3491.3671.3651.361
J (kgm2)M≈3.725 · 10−21.969 · 10−23.479 · 10−22.431 · 10−22.486 · 10−2
Tla (Nm)M≈03.625 · 10−11.269 · 10−12.486 · 10−11.115 · 10−1
Tlb (Nms)M≈4.8 · 10−34.443 · 10−176.333 · 10−173.339 · 10−35.411 · 10−17
Tlc (Nms2)M≈01.554 · 10−52.223 · 10−58.088 · 10−72.262 · 10−5
t (s)M-6.9457.5788.0389.603
Table 13. OF and Mean values of the calculated parameters for 50 independent runs of DE/rand/1/exp using the RK fourth-order method for M1U with a sampling time of 6.6 ms for the first to fifth measurements.
Table 13. OF and Mean values of the calculated parameters for 50 independent runs of DE/rand/1/exp using the RK fourth-order method for M1U with a sampling time of 6.6 ms for the first to fifth measurements.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueFirst Meas.Second Meas.Third Meas.Fourth Meas.Fifth Meas.
OFB-3.8094 · 10−42.9718 · 10−43.2691 · 10−43.5838 · 10−42.9635 · 10−4
W-3.8094 · 10−42.9718 · 10−43.2691 · 10−43.5838 · 10−42.9635 · 10−4
M-3.8094 · 10−42.9718 · 10−43.2691 · 10−43.5838 · 10−42.9635 · 10−4
SD-5.4210 · 10−205.4210 · 10−205.4210 · 10−200.05.4210 · 10−20
Ra (Ω)M5.664.2134.2494.1894.3194.230
La (H)M0.04720.05830.05430.05990.05330.0541
cm (Vs)M1.3561.3601.3591.3601.3591.359
J (kgm2)M≈3.725 · 10−22.576 · 10−22.582 · 10−22.585 · 10−22.580 · 10−22.588 · 10−2
Tla (Nm)M≈03.764 · 10−162.340 · 10−163.327 · 10−163.102 · 10−161.540 · 10−16
Tlb (Nms)M≈4.8 · 10−34.819 · 10−34.834 · 10−34.838 · 10−34.723 · 10−34.765 · 10−3
Tlc (Nms2)M≈06.871 · 10−208.843 · 10−203.989 · 10−205.079 · 10−204.787 · 10−20
t (s)M-9.2128.7958.6689.1358.919
Table 14. OF and Mean values of the calculated parameters for 50 independent runs of DE/rand/1/exp using the RK 4th order for M1U with a sampling time of 6.6 ms for the 6th to 10th measurements.
Table 14. OF and Mean values of the calculated parameters for 50 independent runs of DE/rand/1/exp using the RK 4th order for M1U with a sampling time of 6.6 ms for the 6th to 10th measurements.
OF and
Parameters
MethodDE/rand/1/exp
Known Value6th Meas.7th Meas.8th Meas.9th Meas.10th Meas.
OFB-3.1850 · 10−43.9034 · 10−43.1248 · 10−43.2844 · 10−43.2252 · 10−4
W-3.1850 · 10−43.9034 · 10−43.1248 · 10−43.2844 · 10−43.2252 · 10−4
M-3.1850 · 10−43.9034 · 10−43.1248 · 10−43.2844 · 10−43.2252 · 10−4
SD-0.01.0842 · 10−195.4210 · 10−201.0842 · 10−195.4210 · 10−20
Ra (Ω)M5.664.2514.2674.2064.2404.258
La (H)M0.04720.05460.05740.05630.05890.0596
cm (Vs)M1.3561.3591.3581.3591.3591.359
J (kgm2)M≈3.725 · 10−22.580 · 10−22.578 · 10−22.584 · 10−22.585 · 10−22.583 · 10−2
Tla (Nm)M≈02.023 · 10−162.606 · 10−162.883 · 10−162.284 · 10−163.962 · 10−16
Tlb (Nms)M≈4.8 · 10−34.740 · 10−34.729 · 10−34.744 · 10−34.752 · 10−34.721 · 10−3
Tlc (Nms2)M≈08.943 · 10−208.483 · 10−204.873 · 10−207.036 · 10−206.796 · 10−20
t (s)M-8.6579.2338.8759.2028.919
Table 15. Deviation of OFs and parameters from the overall mean values for 50 × 10 independent runs using the RK fourth-order method for M1U and sampling 6.6 ms.
Table 15. Deviation of OFs and parameters from the overall mean values for 50 × 10 independent runs using the RK fourth-order method for M1U and sampling 6.6 ms.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueLowestHighestMeanSD% dev. to Mean
OF-2.9635 · 10−43.9034 · 10−43.3320 · 10−43.1073 · 10−5−11.06% to +17.15%
Ra (Ω)5.664.1894.3194.2423.4588 · 10−2−1.25% to +1.82%
La (H)0.04720.05330.06000.05672.3667 · 10−3−6.00% to +5.82%
cm (Vs)1.3561.3581.3601.3594.4725 · 10−4−0.07% to +0.07%
J (kgm2)≈3.725 · 10−22.576 · 10−22.588 · 10−22.582 · 10−23.4700 · 10−5−0.25% to +0.23%
Tla (Nm)≈00.03.569 · 10−152.783 · 10−165.3603 · 10−16-
Tlb (Nms)≈4.8 · 10−34.721 · 10−34.838 · 10−34.766 · 10−34.3791 · 10−5−0.94% to +1.51%
Tlc (Nms2)≈00.08.511 · 10−196.570 · 10−201.1562 · 10−19-
t (s)-7.56311.0638.9614.5968 · 10−1−15.60% to 23.46%
Table 16. OF and Mean values of the calculated parameters for 50 independent runs using the RK fourth-order method for M1U, sampling time of 6.6 ms and different lengths of the measured values.
Table 16. OF and Mean values of the calculated parameters for 50 independent runs using the RK fourth-order method for M1U, sampling time of 6.6 ms and different lengths of the measured values.
OF and
Parameters
Method
Known ValueWhole Meas.No Pre-TriggerStart Only0.35 s Start + 0.35 s Operation0.35 s Start + 0.70 s Operation
OFB-3.8094 · 10−44.2093 · 10−41.8684 · 10−31.0974 · 10−37.9138 · 10−4
W-3.8094 · 10−44.2093 · 10−41.8684 · 10−31.0974 · 10−37.9138 · 10−4
M-3.8094 · 10−44.2093 · 10−41.8684 · 10−31.0974 · 10−37.9138 · 10−4
SD-5.4210 · 10−200.02.1684 · 10−192.1684 · 10−192.1584 · 10−19
Ra (Ω)M5.664.2134.2134.2764.2064.208
La (H)M0.04720.05830.05830.05780.05790.0580
cm (Vs)M1.3561.3601.3601.3401.3561.358
J (kgm2)M≈3.725 · 10−22.576 · 10−22.576 · 10−22.534 · 10−22.559 · 10−22.569 · 10−2
Tla (Nm)M≈03.764 · 10−161.544 · 10−155.264 · 10−165.692 · 10−161.083 · 10−15
Tlb (Nms)M≈4.8 · 10−34.819 · 10−34.819 · 10−33.087 · 10−35.060 · 10−34.927 · 10−3
Tlc (Nms2)M≈06.871 · 10−207.404 · 10−203.234 · 10−205.122 · 10−207.601 · 10−20
t (s)M-9.2129.5576.3096.8947.043
Data 40036256109162
Table 17. Differences of OFs and parameters from the overall mean values of the calculated parameters for 50 × 5 independent runs using the RK fourth-order method for M1U with a sampling time of 6.6 ms.
Table 17. Differences of OFs and parameters from the overall mean values of the calculated parameters for 50 × 5 independent runs using the RK fourth-order method for M1U with a sampling time of 6.6 ms.
OF and
Parameters
MethodDE/rand/1/exp
Known ValueLowestHighestMeanSD% dev. to Mean
OF-3.8094 · 10−41.8684 · 10−39.1182 · 10−45.4532 · 10−4−58.22% to +104.91%
Ra (Ω)5.664.2064.2764.2232.6351 · 10−2−0.40 to +1.24%
La (H)0.04720.05780.05830.05812.1877 · 10−4−0.52% to +0.34%
cm (Vs)1.3561.3401.3601.3557.4982 · 10−3−1.11% to +0.37%
J (kgm2)≈3.725 · 10−22.534 · 10−22.576 · 10−22.563 · 10−21.5774 · 10−4−1.13% to +0.51%
Tla (Nm)≈00.01.323 · 10−148.197 · 10−161.6478 · 10−15-
Tlb (Nms)≈4.8 · 10−33.087 · 10−35.060 · 10−34.542 · 10−37.3313 · 10−4−32.03% to +11.40%
Tlc (Nms2)≈00.06.757 · 10−196.046 · 10−201.0705 · 10−19-
t (s)-5.71911.5637.8031.4173−26.71% to +48.19
Table 18. Parameters used for the generated responses and OF and Mean values of the calculated parameters for 50 independent runs using the RK fourth-order method for simulated responses, with a sampling time of 6.6 ms.
Table 18. Parameters used for the generated responses and OF and Mean values of the calculated parameters for 50 independent runs using the RK fourth-order method for simulated responses, with a sampling time of 6.6 ms.
OF and
Parameters
DE/rand/1/exp
Parameters Used for Generation of the ResponsesCalculated Parameters
OFB-2.9207 · 10−18
W-2.9207 · 10−18
M-2.9207 · 10−18
SD-5.4553 · 10−27
Ra (Ω)M5.665.66
La (H)M0.04720.0472
cm (Vs)M1.3561.356
J (kgm2)M3.725 · 10−23.725 · 10−2
Tla (Nm)M04.03 · 10−13
Tlb (Nms)M4.8 · 10−34.80 · 10−3
Tlc (Nms2)M01.22 · 10−12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jesenik, M.; Ravber, M.; Trbušić, M. Innovative Approach for the Determination of a DC Motor’s and Drive’s Parameters Using Evolutionary Methods and Different Measured Current and Angular Speed Responses. Mathematics 2024, 12, 42. https://doi.org/10.3390/math12010042

AMA Style

Jesenik M, Ravber M, Trbušić M. Innovative Approach for the Determination of a DC Motor’s and Drive’s Parameters Using Evolutionary Methods and Different Measured Current and Angular Speed Responses. Mathematics. 2024; 12(1):42. https://doi.org/10.3390/math12010042

Chicago/Turabian Style

Jesenik, Marko, Miha Ravber, and Mislav Trbušić. 2024. "Innovative Approach for the Determination of a DC Motor’s and Drive’s Parameters Using Evolutionary Methods and Different Measured Current and Angular Speed Responses" Mathematics 12, no. 1: 42. https://doi.org/10.3390/math12010042

APA Style

Jesenik, M., Ravber, M., & Trbušić, M. (2024). Innovative Approach for the Determination of a DC Motor’s and Drive’s Parameters Using Evolutionary Methods and Different Measured Current and Angular Speed Responses. Mathematics, 12(1), 42. https://doi.org/10.3390/math12010042

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop