Application of Genetic Algorithm Elements to Modelling of Rotation Processes in Motion Transmission Including a Long Shaft

: Genetic algorithms are used to parameter identiﬁcation of the model of oscillatory processes in complicated motion transmission of electric drives containing long elastic shafts as systems of distributed mechanical parameters. Shaft equations are generated on the basis of a modiﬁed Hamilton–Ostrogradski principle, which serves as the foundation to analyse the lumped parameter system and distributed parameter system. They serve to compute basic functions of analytical mechanics of velocity continuum and rotational angles of shaft elements. It is demonstrated that the application of the distributed parameter method to multi-mass rotational systems, that contain long elastic elements and complicated control systems, is not always possible. The genetic algorithm is applied to determine the coefﬁcients of approximation the system of Rotational Transmission with Elastic Shaft by equivalent differential equations. The ﬁtness function is determined as least-square error. The obtained results conﬁrm that application of the genetic algorithms allow one to replace the use of a complicated distributed parameter model of mechanical system by a considerably simpler model, and to eliminate sophisticated calculation procedures and identiﬁcation of boundary conditions for wave motion equations of long elastic elements.


Introduction
Electric drives are among the top consumers of electricity, estimated to use 60% of the energy generated at present. A variety of drives are used in almost every industrial plant. A generalised electric drive consists of a drive motor, motion transmission system, loading system, and drive control system [1,2]. Each component has an impact on effective operations of the drive. Their adequate cooperation is crucial for optimal operation of the drive.
Drive motion transmission is emphasized in this paper. There are a number of motion transmission types in practice, depending on: motor power, type of motor loading system, system operating conditions, and the requirements and objectives of an electric drive [3][4][5][6][7]. Motion transmission in systems comprising long shafts is an important practical issue, since such systems are often used, e.g., in industrial conditions. These systems are used where loading mechanisms are far away from drive motors, most commonly in connection with production process requirements. This is particularly true of: high-power vertical pumps, turbine sets, lifts, rolling mills, diggers, etc. As motion transmission is the main factor of correct and efficient drive operation, dynamic processes across shafts must definitely be taken into account.
Mathematical modelling of dynamic processes in long motion transmission shafts of electric drives is well discussed [4,[8][9][10] using Lagrange's theory for lumped parameter systems and mechanics of physical continuum. In simple cases, shafts are divided into several absolutely rigid inertial units connected by means of elastic clutches. Each unit is described with second type Lagrange equations. In another, more complicated case, the shaft is treated as a single-dimensional (one spatial coordinate) continual system and described with fractional equations [4,[8][9][10][11]. This approach poses the problem of computing boundary conditions, which is not always an easy task [4,[11][12][13].
Another approach is applied to describe real physical processes in a drive system, since the long shaft equation might address wave processes in the shaft continuum. Application of Lagrange's theory to lumped parameter systems may be sufficient in a majority of cases, while results are satisfactory and agree with reality to a sufficient extent. The theory of continuum systems must be used in some precise drives, however. If, in line with Lagrange's theory, a shaft model is represented as a dual mass holonomic system (two degrees of freedom), velocity and rotation angle diagrams can be plotted for both the masses, in this case, the left and right shaft ends. If the shaft is presented as a continuum, on the other hand, and the results for velocity and rotation angle are compared to results obtained from an equivalent dual mass holonomic system, there will be significant discrepancies between the amplitudes of rotor and loading machine at the shaft ends, and the frequency of changes of these speeds. [14]. An angle δ is a measure of absolute error of the velocities. It will expand as time increases. The extreme of velocity or rotation angle function calculated by means of Lagrange's equations will in fact be shifted by δ. This is very important in drive control. Supplying signal to a drive motor at a wrong instant shifted by δ, for instance to damp resonance, may cause malfunction of the drive system. A somewhat different approach is proposed here to compute real velocities and rotational angles of shaft ends with the aid of the genetic algorithm.

Mathematical Model of an Elastic-Dissipation Long Shaft of Distributed Parameters
A mathematical model of the long shaft is based on an interdisciplinary method of modelling that involves modifications to the integral variational Hamilton-Ostrogradski principle by expanding the conservative Lagrange's function with two components. The first includes energy of dissipation forces, the other energies of external non-potential forces. This method has been introduced in [4]. The idea of modifying the Hamilton-Ostrogradski principle is similar to that applied in [8,9], where two additional components were introduced to the conservative Lagrangian formally, while the introduction gained mathematical foundations in [4]. Additionally, [8,9] analyse lumped parameter dynamic systems, and both lumped and distributed parameter dynamic systems.
Accordingly, the long shaft equation will be formulated as follows [4]: where a-mechanical wave velocity; ρ-density of shaft material; J p -polar shaft moment of inertia; G-shear modulus; ν e , ζ-external and internal shaft dispersion factors; M 0 -load torque; m 0 -load torque along each shaft point. Direct boundary conditions for Euler-Poisson Equation (1), normally based on d'Alembert's theory [4,8,9], are used here to develop a shaft model.
In general, a long shaft is analysed as the motion transmission part of an electric drive. The shaft connects the drive motor and the loading system ( Figure 1).  The following nomenclature is introduced in Figure 1 MEM-electromagnetic torque of drive motor; MN-drive loading torque; JEM, JN-moments of inertia of motor rotor and loading system; x-spatial coordinate; N-number of discrete nodes in (2); Δx-discretisation step of spatial derivatives; {(0) and (N + 1)}-fictitious discretisation nodes of spatial derivatives [4].
For the drive in Figure 1, (1) is simplified to: Boundary conditions for (2) are found from equality of electromagnetic torque, elasticity and dissipation across elementary cross-sections of shaft ends and equality of loading torques, elasticity and dissipation.
The system of Equations (2)-(4) is solved by discretising the spatial derivatives using the straight line method. This ultimately results in: Velocity and rotation angle of shaft start (rotor)-ω1, ϕ1 and their equivalents for the shaft end-ωN, ϕN are the most interesting functions in (5)- (8). These are starting parameters for determining the instant of supplying the control signal in precise electric drives.
The logical question arises of a possible similarity between a long shaft in drive motion transmission and a dual mass elastic dissipation system [15][16][17][18]. It is not always the case, since similarity of a lumped and dispersed parameter system does not address velocity of mechanical wave motion, which leads to inaccuracies in determination of real system parameters. The angle δ appears among the calculation results, introducing a difference between velocities and rotation angles.
A model of the system in Figure 1 is developed in its dual mass version. The following nomenclature is introduced in Figure 1 M EM -electromagnetic torque of drive motor; M N -drive loading torque; J EM , J N -moments of inertia of motor rotor and loading system; x-spatial coordinate; N-number of discrete nodes in (2); ∆x-discretisation step of spatial derivatives; {(0) and (N + 1)}-fictitious discretisation nodes of spatial derivatives [4].
For the drive in Figure 1, (1) is simplified to: Boundary conditions for (2) are found from equality of electromagnetic torque, elasticity and dissipation across elementary cross-sections of shaft ends and equality of loading torques, elasticity and dissipation.
The system of Equations (2)-(4) is solved by discretising the spatial derivatives using the straight line method. This ultimately results in: Velocity and rotation angle of shaft start (rotor)-ω 1 , ϕ 1 and their equivalents for the shaft end-ω N , ϕ N are the most interesting functions in (5)- (8). These are starting parameters for determining the instant of supplying the control signal in precise electric drives.
The logical question arises of a possible similarity between a long shaft in drive motion transmission and a dual mass elastic dissipation system [15][16][17][18]. It is not always the case, since similarity of a lumped and dispersed parameter system does not address velocity of mechanical wave motion, which leads to inaccuracies in determination of real system parameters. The angle δ appears among the calculation results, introducing a difference between velocities and rotation angles. A model of the system in Figure 1 is developed in its dual mass version.
where c 1,2 -rigidity coefficient of equivalent elastic clutch; ν 1,2 -dispersion coefficient of equivalent elastic clutch; L-shaft length; D-shaft diameter. The mathematical model of electric drive motion transmission including long elastic elements is described with (2) and is the most accurate, which means results of computer simulation of functional dependences sought on its basis can be with great adequacy regarded as an experiment. The functions computed in this way will be approximated by Equation (18).

Principles of Calculating Approximation Coefficients
Initial attempts at solving the problem employed the classic approximation method using the least-squares method [19]. Since the proposed approximation function cannot be reduced by substitution to a linear or polynomial function with their ready-made equations, the definition was considered. This approximation problem can be approached as a non-linear least-squares method [19].
The approximation function is: where the vector β stores m of sought factors; x-function argument.
Parameter values need to be found for this function that 'best' match the data. Matching of a model and a data point is measured with its residual value, defined as the difference between the real dependent variable and its value anticipated by the model: where y i − i th value of the real function; r i − i th difference (residue); n-number of measurement points. The least-squares method finds optimum parameter values by minimising the total residues squared: Extreme of function S SUM can be found with the aid of the necessary condition, that is, by comparing its gradient to zero [19]. Since the model contains 'm' parameters β i , the same number of equations results.
The gradient of (16) will be computed considering (15) Dependences of the approximated functions relating to the long motion transmission shaft are determined using the theory of a dual mass rotation system [14]. We additionally where C 1 , . . . , C 4 -approximation factors.
The approximation function is introduced as a general solution to a homogeneous (at C4 = 0) ordinary second-order equation with constant coefficients which describes oscillation processes in a dual mass system. If the right-hand side of the equation is constant, then C4 = 0.
Application of the non-linear least-squares method to find the parameters C1, C2, C3, C4, yields: The iterative Gauss-Newton method serves to solve this system: where F'(C)-Jacobi's matrix, k-number of iteration. The procedure of numerically solving a similar set of non-linear equations in Mathematica software (Wolfram Research, Champaign, IL, USA) was presented in [20,21].
In previous research it was proven that the system could not be solved numerically using the above method. The method often does not converge in practical calculations, since the fitness function (sum of squared residual values) is strongly non-linear, which is the case of our calculations. Thus in previous work we experienced problems with fitting the model to experimental. An improved method was then applied: Levenberg-Marquardt iterative algorithm that combines features of the greatest (gradient) decline and Gauss-Newton methods. It additionally uses the Hessian matrix of the system (H), that is, information on its second fractional derivatives.
where H(C)-Hessian matrix; ∇-Hamilton operator; λ-a certain positive constant. Levenberg-Marquardt algorithm is a powerful calculation method that includes convergence of Newton's method on the one hand and methods of gradient descent on the other. The Levenberg-Marquardt method can be applied to solve equations with significant non-linearity.
Yet the Levenberg-Marquardt method is difficult to be implemented for this application. Therefore we have turned our attention to the alternative method of tackling the problem, namely we focused on the possibility to use Genetic algorithm (GA) for finding the values of model parameters.

Evolutionary Algorithm Approximating Angular Velocity of a Long Shaft
The issue of determining parameters of an approximation function of shaft motion can be seen as a typical problem of global optimisation. We applied classical version of a genetic algorithm [20] to find parameters of a function describing the shaft's angular velocity. Genetic algorithm (GA) is a metaheuristic algorithm that helps to solve a range of complex problems, where a method of solving is not specified or known in advance, yet ways of evaluating quality of solutions are known.
In work [22], the usefulness of the genetic algorithms was for offline identification of electromagnetic and mechanical parameters of a mathematical model of induction machine, taking into account system dynamics (e.g., Coulomb's friction). Identification was performed using data acquired during a test consisting of a transient from standstill to a certain speed and successive free motion to standstill. The work [23] compares the performance of genetic algorithms and the Levenberg-Marquardt method for identification of parameters in the lumped model of an induction motor using transient stator current waveforms. The authors have found out that both methods are directly applicable to a variety of estimation problems found in the analysis of equivalent diagrams of electric machines. In work [24] it was used genetic algorithms for parameter identification for field control of induction motors. The authors have found out that the solutions of the problem depended heavily on the mathematical model and its measured data. They suggested to use the model equations in relative units in order to improve the accuracy and the efficiency of the identification method. The work [25] uses the real version of genetic algorithms implemented in the free Matlab toolkit (Mathworks Inc., Natick, MA, USA).
In the context of the subject of the present paper it is worth to notice that a somewhat similar identification problem in civil engineering (parameter estimation of internal thermal mass of building dynamic models for prediction of transient heating or cooling in a distributed network model of a building) has been tackled successfully in the paper by [26]. In work [27], the parameters of the hysteresis model for both soft and hard ferromagnets were estimated using the genetic algorithm method. The authors compared the performance of three meta-heuristic approaches (Genetic Algorithms and two more recent optimization techniques, namely Differential Evolution and Artificial Bee Colony).
On the basis of the afore-given literature review it can be stated that despite the genetic algorithms technique is a bit outdated and in some applications it has been successfully superseded with more recent optimization methods, it is still a useful tool for engineering purposes. A state-of-the-art paper, [28], comments on the advances in the field of natureinspired metaheuristics and we allow ourselves to refer the readers more interested in the subject to it for a detailed information on the advantages and the shortcomings of these methods.
Limiting possible values of all approximation parameters to a relatively narrow range is key to the optimisation of a problem. This means certain restrictions must be placed on values subject to genetic operators. These restrictions were initially determined on the basis of shaft motion simulation graphs.
This approach sees the process of optimisation as universal and nearly independent from a specific problem at hand. The objective function is designed in such a way as to meet typical requirements of solutions proposed. Finding a set of approximation function factors for which the total reference value deviations from values derived from approximation (absolute errors) squared is minimum is a universally acceptable criterion of the approximation problem solution. In our case, data from the shaft simulation computed on their basis of (5)-(8) are the point of reference. In earlier studies the objective function referred to the maximum relative error of approximation.
The objective function in the case of C 4 = 0, see (18), becomes: and in the case C 4 = 0 A general diagram of the genetic algorithm employed in this paper is illustrated in Figure 2.
A general diagram of the genetic algorithm employed in this paper is illustrated in Figure 2. In our simulations we have availed of freely available GAOT (Mathworks Inc., Natick, MA, USA) toolbox, which can be applied with practically any version of Matlab (more recent versions of Matlab offer a Global Optimization toolbox, which contains implementations of genetic algorithms and pattern search). We have found that the GAOT toolbox is simple to use and sufficient for our needs.
Evaluation of all population members is based on selection of a function objective minimisation. A range of simulations were undertaken before finally selecting values of the different parameters included in Table 1 and used in our approach based on the genetic algorithm.
The upper and lower boundaries were determined on the basis of simulation testing. C1 depends on function amplitude, C2 on dissipation coefficient (where dissipation is absent, C2 = 0), C3 depends on frequency of process oscillation (i.e., coefficient of shaft rigidity). Root mean square error (Δ) is calculated: In our simulations we have availed of freely available GAOT (Mathworks Inc., Natick, MA, USA) toolbox, which can be applied with practically any version of Matlab (more recent versions of Matlab offer a Global Optimization toolbox, which contains implementations of genetic algorithms and pattern search). We have found that the GAOT toolbox is simple to use and sufficient for our needs.
Evaluation of all population members is based on selection of a function objective minimisation. A range of simulations were undertaken before finally selecting values of the different parameters included in Table 1 and used in our approach based on the genetic algorithm. The upper and lower boundaries were determined on the basis of simulation testing. C1 depends on function amplitude, C2 on dissipation coefficient (where dissipation is absent, C2 = 0), C3 depends on frequency of process oscillation (i.e., coefficient of shaft rigidity). Root mean square error (∆) is calculated: where ω a -function approximated by means of GA. Rotation angles for the shaft ends are important dependences. They are derived from: dt [4]: Integration of (18) by parts and the necessary transformations produce an expression serving to compute the function of shaft rotation angle:

Results of Computer Simulation
The electromechanical system was presented as a motion transmission including a long shaft. The torque generated by the electric motor was applied to the left end and the loading system to the right end. To simplify the calculations, the state with a constant starting torque and loading is taken into consideration. Three experiments were undertaken with the above torques applied to the shaft ends (c.f. Figure 1). In the first experiment, the torques are: M EM = M N = 1000 N·m, and in the second M EM = 1000 N·m M N = 1004 N·m. If moments of inertia are equal, coefficients of approximation will be equal in terms of their absolute values (after the modulus), but with different signs. This is the general approach that is used in the paper for different moments of inertia. The drive will exhibit zero velocity and a constant shaft torsion angle in a steady state as part of experiment one, whereas it will accelerate in the steady state as part of the second experiment. When the third experiment is in progress, given the same torques as in the second experiment, the rotation angle is determined on the basis of the dual mass model (see Equations (9)-(13)) and the genetic algorithm. The parameters of the long shaft and the reference unit are shown in Table 2.

Results of Computer Simulation for the Dual Mass Model
The system of differential Equations (5)-(8) was solved by means of 4th order Runge-Kutta method including the time step ∆t = 5·10 −7 s. The system (9)-(11) was also solved by means of 4th order Runge-Kutta method, though with a different time step ∆t = 5 × 10 -5 s. This shift of the time step helped to avoid latent integration methods [4].   (18). This dependence is virtually identical in both the cases of modelling. For the right shaft end (towards the load), the speed function is practically the same, though in counter-phase to the left end's velocity. The graph shows the system enters a steady state after 5 s. Since both the torque moments are identical in experiment one, the shaft is not turning.   (18). This dependence is virtually identica in both the cases of modelling. For the right shaft end (towards the load), the speed function is practically the same, though in counter-phase to the left end's velocity. Th graph shows the system enters a steady state after 5 s. Since both the torque moments ar identical in experiment one, the shaft is not turning.   8) and (29). Comparative analysis of the graphs demon strates operating principles of the transmission. First, oscillation rotation angles are sim ilar and both the functions are in counter-phase to each other. Second, the left and righ shaft ends in the steady state can be seen to be turned by an angle of approximately 0.00 rad relative to the zero point of the system of coordinates. This means the total angle o shaft torsion in relation to its ends is approximately 0.008 rad.                   Energies 2021, 14, x FOR PEER REVIEW 11 Instant rotation angles of both the shaft in the time interval t ∈ (0; 1) s are prese in Figures 7 and 8. A comparison detects major differences. Both the functions not are in counter-phase but also display other amplitudes in the time range under dis sion.    Energies 2021, 14, x FOR PEER REVIEW 11 Instant rotation angles of both the shaft in the time interval t ∈ (0; 1) s are prese in Figures 7 and 8. A comparison detects major differences. Both the functions not are in counter-phase but also display other amplitudes in the time range under dis sion.       For the sake of improved visualisation of physical processes in the drive motion transmission in Figures 9 and 10, the same dependences are illustrated in the time interval t ∈ (0.85; 1) s.
The time delay angle δ can be seen in both the graphs. The physical interpretation of this process is well known. A mechanical wave moving along a shaft produces the delay. Meanwhile, the transmission model based on a dual mass system assumes no such notion as velocity of wave motion. It is evident the rotation angle function based on the dual mass model does not correspond to the real value of δ. This causes supply of control signal to the drive motor to be inaccurate or occasionally incorrect in time. This is one of the main reasons why other approaches to simplification of the motion transmission model are employed as part of distributed parameter systems. Application of the genetic algorithm has broad prospects in distributed parameter systems.
Results of relative error analysis of approximating formulas for the first and second experiments at given angular velocities of both the shaft ends are included in Table 3.  For the sake of improved visualisation of physical processes in the drive motion transmission in Figures 9 and 10, the same dependences are illustrated in the time interval t ∈ (0.85; 1) s.
The time delay angle δ can be seen in both the graphs. The physical interpretation of this process is well known. A mechanical wave moving along a shaft produces the delay. Meanwhile, the transmission model based on a dual mass system assumes no such notion as velocity of wave motion. It is evident the rotation angle function based on the dual mass model does not correspond to the real value of δ. This causes supply of control signal to the drive motor to be inaccurate or occasionally incorrect in time. This is one of the main reasons why other approaches to simplification of the motion transmission model are employed as part of distributed parameter systems. Application of the genetic algorithm has broad prospects in distributed parameter systems.
Results of relative error analysis of approximating formulas for the first and second experiments at given angular velocities of both the shaft ends are included in Table 3. For the sake of improved visualisation of physical processes in the drive motion transmission in Figures 9 and 10, the same dependences are illustrated in the time interval t ∈ (0.85; 1) s.
The time delay angle δ can be seen in both the graphs. The physical interpretation of this process is well known. A mechanical wave moving along a shaft produces the delay. Meanwhile, the transmission model based on a dual mass system assumes no such notion as velocity of wave motion. It is evident the rotation angle function based on the dual mass model does not correspond to the real value of δ. This causes supply of control signal to the drive motor to be inaccurate or occasionally incorrect in time. This is one of the main reasons why other approaches to simplification of the motion transmission model are employed as part of distributed parameter systems. Application of the genetic algorithm has broad prospects in distributed parameter systems.
Results of relative error analysis of approximating formulas for the first and second experiments at given angular velocities of both the shaft ends are included in Table 3. Duration of the experiments was t ∈ (0; 5) s. The number of approximation points: n = 1000. The table provides simulation results for the time interval t ∈ (0; 25) s. Relative errors ε are computed as follows:  Table 3 indicates values of the angular velocity are virtually identical. Root mean square error for the approximating formulas in both the experiments was also calculated using (28). The results are part of Table 4.

Experimental Testing and Verification of Calculation Results
The experimentation was conducted on a laboratory stand for testing of drive systems including long elastic elements [29]. A flow diagram of the stand is shown in Figure 11.   Table 3 indicates values of the angular velocity are virtually identical. Root mean square error for the approximating formulas in both the experiments was also calculated using (28). The results are part of Table 4.

Experimental Testing and Verification of Calculation Results
The experimentation was conducted on a laboratory stand for testing of drive systems including long elastic elements [29]. A flow diagram of the stand is shown in Figure  11. Figure 11. Flow diagram of a laboratory stand for testing of drive systems including elastic connections.
The drive system in testing consists of a three-phase induction motor connected to loading with an elastic element. The moment of motion transmission system inertia is Figure 11. Flow diagram of a laboratory stand for testing of drive systems including elastic connections.
The drive system in testing consists of a three-phase induction motor connected to loading with an elastic element. The moment of motion transmission system inertia is modified by mounting additional steel discs on the shaft towards the loading system. The laboratory stand is equipped with two pulse angle position encoders (incremental encoders) with 3600 pulses per rotation. The moment is measured with a torque meter affixed to the motor shaft. Signals from measurement sensors are transmitted to a computer having two National Instruments PCle-6351 measurement cards. LabVIEW (National Instruments, Austin, TX, USA) software is employed to gather data.
These    modified by mounting additional steel discs on the shaft towards the loading system. The laboratory stand is equipped with two pulse angle position encoders (incremental encoders) with 3600 pulses per rotation. The moment is measured with a torque meter affixed to the motor shaft. Signals from measurement sensors are transmitted to a computer having two National Instruments PCle-6351 measurement cards. LabVIEW (National Instruments, Austin, TX, USA) software is employed to gather data.
These Numerical calculations were carried out for a shaft length 4.5 m and diameter 0.05 m. Experimentation on a shaft of such large dimensions involves a great risk of torsional moments that may cause a breakdown. Therefore, a long elastic shaft of a length 0.66 m and thickness 0.025 m was used in the testing.
Laboratory results are presented graphically and compared to the numerical calculations. Examples of laboratory test results compared to the corresponding numerical results are illustrated in Figure 12, Figure 13, Figure 14, Figure 15. The drive system was loaded with a default moment ML = 81 N·m.    A simulation for these measurement cases was undertaken using the authors' mathematical models in order to determine degree of matching between the calculation and measurement results. The verification concerned instant values of elasticity moment across the element connecting the motor and loading and the motor's angular velocity.
In order to compare calculation results and measurement results recorded on the test stand, mean values of the first five amplitudes of elasticity moment and angular velocity oscillations were determined. The difference between numerical calculation results and measurement results was 3.1% for elasticity moment and 19.1% for angular velocity, evidence of a high adequacy of our mathematical models based on variational approaches. Such accuracy is sufficient for the purposes of designing new and modernisation of existing drive systems.

Conclusions
Elements of genetic algorithm can be used with success with the aid of parametric identification of approximation coefficients for functions of rotational speed and rotation angles of ends of motion transmission across long elastic elements. This avoids the complicated theory of oscillation process modelling in long shafts at the level of the distributed parameter system.
Experience shows replacement of long elastic elements with distributed mechanical parameters with multiple mass elastic-dissipative units of lumped parameters, for instance elastic clutches, is not always right. This is connected with the fact a multiple (e.g., dual) mass interpretation of an equivalent model fails to include the important phenomenon of mechanical wave motion along a shaft, which will delay the oscillation  A simulation for these measurement cases was undertaken using the authors' mathematical models in order to determine degree of matching between the calculation and measurement results. The verification concerned instant values of elasticity moment across the element connecting the motor and loading and the motor's angular velocity.
In order to compare calculation results and measurement results recorded on the test stand, mean values of the first five amplitudes of elasticity moment and angular velocity oscillations were determined. The difference between numerical calculation results and measurement results was 3.1% for elasticity moment and 19.1% for angular velocity, evidence of a high adequacy of our mathematical models based on variational approaches. Such accuracy is sufficient for the purposes of designing new and modernisation of existing drive systems.

Conclusions
Elements of genetic algorithm can be used with success with the aid of parametric identification of approximation coefficients for functions of rotational speed and rotation angles of ends of motion transmission across long elastic elements. This avoids the complicated theory of oscillation process modelling in long shafts at the level of the distributed parameter system.
Experience shows replacement of long elastic elements with distributed mechanical parameters with multiple mass elastic-dissipative units of lumped parameters, for instance elastic clutches, is not always right. This is connected with the fact a multiple (e.g., dual) mass interpretation of an equivalent model fails to include the important phenomenon of mechanical wave motion along a shaft, which will delay the oscillation A simulation for these measurement cases was undertaken using the authors' mathematical models in order to determine degree of matching between the calculation and measurement results. The verification concerned instant values of elasticity moment across the element connecting the motor and loading and the motor's angular velocity.
In order to compare calculation results and measurement results recorded on the test stand, mean values of the first five amplitudes of elasticity moment and angular velocity oscillations were determined. The difference between numerical calculation results and measurement results was 3.1% for elasticity moment and 19.1% for angular velocity, evidence of a high adequacy of our mathematical models based on variational approaches. Such accuracy is sufficient for the purposes of designing new and modernisation of existing drive systems.

Conclusions
Elements of genetic algorithm can be used with success with the aid of parametric identification of approximation coefficients for functions of rotational speed and rotation angles of ends of motion transmission across long elastic elements. This avoids the complicated theory of oscillation process modelling in long shafts at the level of the distributed parameter system.
Experience shows replacement of long elastic elements with distributed mechanical parameters with multiple mass elastic-dissipative units of lumped parameters, for instance elastic clutches, is not always right. This is connected with the fact a multiple (e.g., dual) mass interpretation of an equivalent model fails to include the important phenomenon of mechanical wave motion along a shaft, which will delay the oscillation process in a clutch being modelled. This physical phenomenon is at the root of the theory of infinite degrees of freedom in lumped parameter systems; contrary to multiple mass models where number of degrees of freedom is pre-determined and equal to number of generalised coordinates for holonomic systems.
The computer simulation results give reasons to conclude the genetic algorithm is correctly used in this paper to identify function parameters of rotation speed and angle of shaft ends in the drive's motion transmission. The error of approximation is to the order of 10 −4 for the first experiment and 10 −6 for the second, which confirms the long shaft model developed by the authors can be used in precise operating conditions of electric drives containing long elastic elements of distributed parameters.