Next Article in Journal
CBLSTM-AE: A Hybrid Deep Learning Framework for Predicting Energy Consumption
Next Article in Special Issue
Mathematical Modelling of Transient Processes in a Three Phase Electric Power System for a Single Phase Short-Circuit
Previous Article in Journal
Layout Optimization of a Modular Floating Wind Farm Based on the Full-Field Wake Model
Previous Article in Special Issue
The Comparison of Microwave Reflectance of Graphite and Reduced Graphene Oxide Used for Electronic Devices Protection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Voltage Stabilisation of a Drive System Including a Power Transformer and Asynchronous and Synchronous Motors of Susceptible Motion Transmission

1
Faculty of Transport, Electrical Engineering and Computer Science, University of Technology and Humanities in Radom, 26-600 Radom, Poland
2
Faculty of Electrical Engineering, Czestochowa University of Technology, Al. Armii Krajowej 17, 42-201 Czestochowa, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(3), 811; https://doi.org/10.3390/en15030811
Submission received: 25 November 2021 / Revised: 16 January 2022 / Accepted: 20 January 2022 / Published: 23 January 2022

Abstract

:
A mathematical model is developed of the master circuit of an electric driver system including a power transformer and susceptible motion transmission of asynchronous and synchronous drives. Electric motors drive water pumps by means of motion transmission that comprises two elastic couplings of lumped mechanical parameters and a long shaft of distributed mechanical parameters. Differential equations for oscillatory processes for the long shaft and the elastic couplings are different. The shaft is described with partial derivative Euler–Poisson equations, which, combined with the boundary conditions, form mixed problems from the mathematical point of view. The elastic couplings, on the other hand, are described with the ordinary second type Lagrange equations. Based on the theory of electromagnetic field, the partial differential equations describe the skin effects across the rotor age bars. Vertical pumps are presented by means of a loading torque waveform as a function of the input shaft’s angular velocity. The complex mathematical model serves to analyse electromechanical transient processes across the integrated drive system. Starting from there, conditions of stabilisation of the drive system voltage are determined. Electromechanical state equations are presented in the normal Cauchy form and integrated using the fourth-order Runge–Kutta method. Results of computer simulations are shown with graphics, which are interpreted and described.

1. Introduction

The process of water supply plays an important part in contemporary society. There is a global trend towards high-quality centralised water supply for increasing numbers of consumers. High-power pumping systems operate in the metallurgic, steelmaking, chemical, and electric power industries, for instance. Virtually every large plant has dedicated water and sewage divisions tasked with providing water for production processes [1,2,3,4,5]. Asynchronous drives are normally used to power-pumping systems that use susceptible motion transmissions to drive pump rotors. Asynchronous motors have a range of advantages over other electric motors, although they suffer from the fundamental drawback of consuming high magnetising currents that can reach values of several times the rated currents in transient processes [6,7,8]. Motors in pumping systems are supplied from power transformer secondary windings. The way the power supply grid will respond to drive system processes and quality of electricity needs to be analysed and described, in particular, to what value supply voltage of pump drives will reduce, especially in transient states [9,10,11].
A variety of methods are employed to stabilise voltage. One consists in reducing the power transformer turns ratio, which will in turn increase the supply voltage of motors by way of raising voltage across a power transformer’s secondary winding. In other words, depending on a dynamic state of a drive system, voltage supplied to motors can be adjusted by modifying the turn’s ratio. This approach can be applied to rigid grids. Where grids are not rigid, particularly given many medium and high-power drives, currents across power transformers will grow a lot and, consequently, a plant’s electric power grid will sustain a greater loading.
Another approach is based on a completely different idea—not of changing energetic properties of a system but introducing additional stabilising elements directly to the power supply of a drive system. From the perspective of applied physics, currents flow across these extra elements which, in the steady state, are in counterphase to currents of the power supply system. According to the first Kirchhoff law, total current across an electric unit will fall while voltage will rise. This is the so-called passive power compensation in a drive system in the steady state. Static capacitor batteries are commonly used as compensation devices. The method gets more complicated in the case of high-power drive systems supplied from medium-voltage grids, which requires application of rather costly compensation equipment [12,13,14,15].
There is also a third approach to stabilising voltage in a drive system, which is essentially different to the two described above. If practicable, part of asynchronous drive motors should be replaced with synchronous motors as early as the stage of a drive system design [9,16,17,18,19,20]. In static states, the angle between the vector of a drive system supply voltage and the vector of synchronous machine phase current can be controlled with the well-known V-characteristics (Mordey curves). In transient states, varying excitation currents of synchronous machines lead to regulation effects and consequently to stabilisation of a drive system voltage.
Each of these three methods has its advantages and disadvantages. It is not our intention to analyse them in detail here. A final decision to select a method of, and a solution to, voltage stabilisation of a drive system should evidently address both economic and energetic analyses. We are presenting an energetic analysis of a drive system’s voltage stabilisation by replacing some asynchronous motors with their equivalent synchronous motors. The fact that a variety of oscillation processes arise in electric drives including susceptible motion transmission which impair operations of drive systems in resonant and near-resonant (vibration rumble) states is an important argument in the general analysis [21,22,23,24,25].
The aim of this work is mathematical modelling of transient processes in an electric load assembly. The elements of the assembly are a power transformer and asynchronous and synchronous pump drives with flexible transmission of motion with distributed parameters. Additional goals are:
  • Analysis of the power system depending on the value of the excitation current of synchronous motors, taking into account the oscillating mechanical processes in the transmission of motion;
  • Simulation test of load assembly states depending on electrical and mechanical parameters.
We analyse pump systems of susceptible motion transmission and lumped and distributed mechanical parameters in drive systems. Thus, the purpose of this study is to analyse transient processes in a drive system including a power transformer, asynchronous and synchronous pump systems of susceptible motion transmission, and long torsional elements of lumped and distributed mechanical parameters.
In order to achieve the above-mentioned goal, mathematical models of all components of the electric load assembly were developed. The main approach of formulating the model is to use the concept of electromechanical energy conversion: electrical components of devices are described by electromagnetic equations taking into account the phenomenon of current displacement in deep-bar cage induction motors. The mechanical components of the electromechanical system are described by the theory of analytical mechanics as components with distributed parameters. On this basis, one comprehensive model of the electric load assembly was developed, with the help of which simulation tests were carried out.
The paper presents the results of the computer simulation, which are analysed and interpreted. The accuracy of the above-mentioned results was verified on the basis of a comparison with the real waveforms occurring in the components of the electric load assembly and included in our articles [11,23,24,26,27].

2. Mathematical Model of an Electromechanical System

The symbols in Figure 1 stand for: e(t)—columnar vector of phase electromotive forces of a rigid electric power grid, T—power transformer, u10—columnar vector of phase voltages supplied to the power transformer’s primary winding, u20—columnar vector of phase voltages supplied to the power transformer’s secondary winding, V = u20—columnar vector of phase drive system voltages, AM—asynchronous drives, SM—synchronous drives, N—number of asynchronous drives in the drive system, and M—number of synchronous drives in the drive system.
Motion transmission of pumping systems of both the drive types is depicted in Figure 2.
Figure 2 contains the following nomenclature: MEM—electromagnetic torque of the drive motor, MP—pump’s hydraulic torque, SP1—flexible coupling between the rotor and motion transmission long shaft (the first coupling), SP2—flexible coupling between the motion transmission long shaft and the vertical pump’s input rotor (the second coupling), C1,2, ν1,2—rigidity and leakage coefficients of the first coupling, CN−1,N, νN−1,N—rigidity and leakage coefficients of the second coupling, kT—value of the reducer’s gear.
Mathematical models of the drive system elements will be presented on the basis of Figure 1 and Figure 2 and of [9,10,21,28,29,30,31,32]. A general interdisciplinary theory modifies the Hamilton–Ostrogradsky principle by expanding the Lagrangian with two components: energy of internal and external dissipation and energy of external non-potential forces is applied to mathematical modelling of the drive system’s electric elements. The method is described in [9]. As distinct from the well-known methods, discussed, e.g., in [22,23,29], it enables analysis of both lumped and distributed parameter systems addressing non-linear dependencies between the vectoral functions of the scalar argument. Thus, drive system elements are modelled using the modified interdisciplinary Hamilton–Ostrogradsky principle.
Use of the first Kirchhoff law allows to ignore three-phase equation elements in the model. We assume, therefore, the electromechanical system in Figure 1 works in a symmetrical state and uncompensated associated fluxes in the non-magnetic zone of the power transformer can be ignored [9,11].
When the power transformer is described with equations, parameters of the primary transformer winding are reduced to those of the secondary winding including the transformer voltage ratio [11].
d i 10 d t = A 11 , 0 ( u 10 r 10 i 10 ) + A 12 , 0 ( V r 20 i 20 )
d i 20 d t = A 21 , 0 ( u 10 r 10 i 10 ) + A 22 , 0 ( V r 20 i 20 ) , V u 20
A 11 = α σ 1 ( 1 α σ 1 G ) , A 12 = A 21 = α σ 1 α σ 2 G , A 22 = α σ 2 ( 1 α σ 2 G )
G ( α σ 1 A + α σ 2 A + ρ A ) 1 ( α σ 1 B + α σ 2 B + ρ B ) 1
i m A + i m B + i m C = 0 , ψ m A + ψ m B + ψ m C = 0
i m i m A i m B T = i 1 A + i 2 A i 1 B + i 2 B T
τ τ A τ B = diag ( τ A ; τ B ) , τ j = ( ψ m , j i m , j ) 1
ρ τ ρ A ρ B = diag ( ρ A ; ρ B ) , ρ j τ j = ( ψ m , j i m , j ) 1 , j = A , B
where i10, i20—columnar vectors of phase currents of the transformer’s primary and secondary windings, respectively, r0—matrix of winding resistances, ασ1, ασ2—matrices of reverse transformer leakage inductances, A0—matrix of reverse transformer inductances, im—columnar vector of transformer magnetising currents, ψm,j—working associated transformer fluxes, to be determined by means of the power transformer magnetising curve, τ—matrix of reverse magnetising inductances, ρ—matrix of reverse differential magnetising inductances, and j = A,B—indices of three-phase system phases.
When modelling the asynchronous motor, parameters of rotor winding (virtual equivalent three-phase system) are reduced to the stator winding parameters considering the current and voltage gear ki, ku. Deep bar cage asynchronous motors are adopted for the purposes of analysis, where current displacement across the rotor cage bars is addressed [9,30] based on the electromagnetic field theory, which helps to calculate the voltage across the rotor cage bars with ordinary and partial differential equations:
d i S , k d t = A S , k ( V r S , k i S , k ) + A S R , k ( u R , k Ω k Ψ R , k r R L , k i R , k )
d i R , k d t = A R S , k ( V r S , k i S , k ) + A R , k ( u R , k Ω k Ψ R , k r R L , k i R , k ) + Ω k i R , k , V u 20
A S = α σ S ( 1 α σ S G ) , A S R = α σ S α σ R L G Π 1 , A R S = Π 1 1 α σ S α σ R L G , A R = Π 1 1 α σ R L ( 1 α σ R L G ) Π 1
G S = α σ S G , G R = α σ R L G , G = B + b A i A b B i A b A i B B + b B i B
Π 1 = 2 3 sin ( k T p 0 γ + 120 ° ) sin ( k T p 0 γ ) sin ( k T p 0 γ ) sin ( k T p 0 γ 120 ° ) , Π 1 1 = 2 3 sin ( k T p 0 γ 120 ° ) sin ( k T p 0 γ ) sin ( k T p 0 γ ) sin ( k T p 0 γ + 120 ° )
B = 1 / ( α σ S + α σ R L + τ ) , R = 1 / ( α σ S + α σ R L + ρ ) , b = 2 3 R B i m 2 , b A = b ( 2 i A + i B ) , b B = b ( 2 i B + i A )
i m = 2 ( i A 2 + i A i B + i B 2 ) / 3 , τ = 1 L m = ( ψ m i m ) 1 , ρ = ( ψ m i m ) 1 , Ω 1 = k T p 0 ω 3 1 2 1 2
i i A i B T = i S A i S B T + Π 1 i R A i R B T = i S A i S B T + i R A Π 1 i R B Π 1 T , Ψ R = 1 α σ R L i R + 1 τ i
where ΨR—columnar vector of main associated fluxes of the rotor, Π1—matrix of angle transformations, Ω1—matrix of angular rotational speeds of asynchronous motor in angle coordinates, uR—columnar vector of voltages across the rotor cage bars, p0—number of machine pole pairs, γ—angle of rotor rotation, ω—angular velocity of the rotor, im—spatial vector of motor magnetisation, ψm (im)—working flux of machine magnetisation, Lm—inductance of machine magnetisation, and k = N—number of asynchronous motors in the drive system.
Voltages across the rotor cage bars will be determined according to the theorem of voltage on the deep bar wire surface [9]:
u R = E ( 0 ) l
where E(0)—electric intensity on the deep bar surface, and l—length of a deep bar.
E(0) will be derived from Maxwell’s theory [9]:
H t = ν γ × × H , E = 1 γ × H , ν = 1 μ
where H, E—vectors of the magnetic and electric field intensity, respectively, γ—specific electric conductance of a wire in a rotor deep bar; μ—magnetic permeability of conductor material in a rotor deep bar.
In a cylindrical system of coordinates, (18) will be expressed as [9]:
H j t = ν γ 2 H j z 2 , E n = k u k i γ H j z , j = A , B
where ku, ki—voltage and current reduces coefficients of machine transformation, and z—coordinate oriented into a cylindrical system deep bar.
Boundary conditions for the first expression in Equation (19) will be determined according to the law of current surge [9]:
H j | z = 0 = i R , n a , H j | z = h = 0 , j = A , B
where a—deep bar width, and h—bar depth in relation to the coordinate z.
Discretisation of the first expression in Equation (19) by means of the straight-line method will produce:
d H j , n d t = ν γ ( Δ z ) 2 ( H j , n 1 2 H j , n + H j , n + 1 ) , j = A , B , n = 2 , 3 , , T 1
where Δz = h/(T − 1)—discretisation step, and T—number of discretisation nodes relative to the coordinate z. Normally, T ≥ 10.
Approximating the spatial derivative in the second expression of Equation (19) using Lagrange’s formula [9] produces the needed function of voltage across the rotor cage bars, see (17):
u R , j = l k u k i 2 γ Δ z ( 3 H j , 1 + 4 H j , 2 H j , 3 ) , j = A , B
Electromagnetic torque of the asynchronous motor will be determined in the usual way [9]:
M E M = 3 k T p 0 ( i S B i R A Π i S A i R B Π ) / τ
When modelling the synchronous motor, stator (armature) winding parameters are reduced to those of the rotor winding considering the current and voltage gear ki, ku. Synchronous motors with salient rotor poles are analysed. We note the synchronous motor model is different to that of the asynchronous motor. There are three currents in the synchronous motor: two in the longitudinal axis d-d and one in the transverse axis q-q, while the phase C current across the stator is excluded. Then, in order to meet the law of matrix multiplication, a topological matrix B is introduced:
d i S , m d t = A S , m ( V r S , m i S , m Ω m Ψ S , m ) + A S R , m B T ( u R , m r R , m i R , m ) + Ω m i S , m
d i R , m d t = B A R S , m ( V r S , m i S , m Ω m Ψ S , m ) + A R , m ( u R , m r R , m i R , m )
A S = Π 1 α σ S ( 1 α σ S G ) Π , A S R = α σ S B T α σ R B G Π , A R S = Π 1 α σ S B T α σ R B G , A R = α σ R ( 1 α σ R B G )
Ψ S ( Ψ S A , Ψ S B ) T , Ψ R ( Ψ D , Ψ Q , Ψ f ) T , i S ( i S A , i S B ) T , i R ( i D , i Q , i f ) T
u S ( u S A , u S B ) T , u R ( 0 , 0 , u f ) T , r S diag ( r S , r S ) , r R diag ( r D , r Q , r f )
L σ S diag ( L σ S , L σ S ) , L σ R diag ( L σ D , L σ Q , L σ f ) , L m diag ( L d , L q ) , B T [ 1 0 1 0 1 0 ]
α m = L m 1 = diag ( α d , α q ) = diag ( 1 L d , 1 L q ) , G = 1 α d + α σ S + α σ f + α σ D 1 α q + α σ S + α σ Q , Ω = k T p 0 ω 1 1
Π 2 3 sin ( k T p 0 γ 120 ° ) sin ( k T p 0 γ ) cos ( k T p 0 γ 120 ° ) cos ( k T p 0 γ ) , Π 1 cos ( k T p 0 γ ) sin ( k T p 0 γ ) cos ( k T p 0 γ 120 ° ) sin ( k T p 0 γ 120 ° )
ψ S = Π 1 ψ m = L m i m = L m ( i S + Π 1 B T i R )
where Ω—matrix of angular rotational speeds of the synchronous machine in rectangular coordinates, Π—Park matrix of rectangular transformations, d—index related to parameters in the machine’s longitudinal axis, q—index related to parameters in the machine’s transverse axis, f—index related to parameters of the excitation winding, D—index related to parameters of the damping winding in the longitudinal axis d, Q—index related to parameters of the damping winding in the transverse axis q, and m = M—number of synchronous drives in the drive system.
Electromagnetic torque of the synchronous motor will be computed in the ordinary manner [9]:
M E M = 3 k T p 0 ( ψ S B i S A ψ S A i S B )

3. Mathematical Model of a Pumping System’s Motion Transmission

Motion transmissions in asynchronous and synchronous drives are presented as torsional elements of lumped and distributed parameters. We will then also present a mathematical model of an assembly on the basis of the generalised interdisciplinary method [9] that serves to develop the models of a power transformer and drive motors. For an improved visualisation of the modelling process of the mechanical subsystem, a calculation scheme of the motion transmission is shown in Figure 3.
Figure 3 includes the following symbols: JEM—moment of inertia of the motor’s drive rotor, JP—moment of inertia of the hydraulic pump rotor recalculated considering the reducer gear, y, ω1, φ1—angular velocity and rotation angle of the drive motor rotor, ωN, φN—angular velocity and rotation angle of the vertical pump rotor, Δx—discretisation step of the long shaft equation, *—fictitious discretization node, and N—number of discretisation nodes.
The equation of oscillations of the drive motion transmission long shaft is as follows [24,28]:
d ω 1 d t = 1 k T 2 J E M ( c 1 , 2 ( φ 2 φ 1 ) + ν 1 , 2 ( ω 2 ω 1 ) + k T M E M )
d ω N d t = 1 J N ( M P c N 1 , N ( φ N φ N 1 ) ν N 1 , N ( ω N ω N 1 ) )
2 φ ( x , t ) t 2 = G ρ 2 φ ( x , t ) x 2 + ξ ρ J p 2 φ ( x , t ) x 2
where G—transverse torsion module, ξ—coefficient of shaft internal leakage; x—current coordinate along the shaft; ρ—density of shaft material; Jp—polar moment of shaft inertia, and φ(x,t)—angle of shaft rotation.
Poincaré boundary conditions of the third type [25] are formulated on the basis of d’Alembert law [9], cf. Figure 2 and Figure 3 will be attached to the shaft Equation (36):
G J p φ ( x , t ) x | x = 0 + ξ ω ( x , t ) x | x = 0 + c 1 , 2 ( φ 1 φ 2 ) + ν 1 , 2 ( ω 1 ω 2 ) = 0
G J p φ ( x , t ) x | x = l + ξ ω ( x , t ) x | x = l c N 1 , N ( φ N 1 φ N ) ν N 1 , N ( ω N 1 ω N ) = 0
Discretisation of (36)–(38) with the straight-line method will produce:
d ω 2 d t = 2 ( c 1 , 2 ( φ 1 φ 2 ) + ν 1 , 2 ( ω 1 ω 2 ) ) ρ J p Δ x 2 ( J p G ( φ 2 φ 3 ) + ξ ( ω 2 ω 3 ) ) ρ J p ( Δ x ) 2
d ω N 1 d t = 2 ( J p G ( φ N 2 φ N 1 ) + ξ ( ω N 2 ω N 1 ) ) ρ J p ( Δ x ) 2 2 ( c N 1 , N ( φ N 1 φ N ) + ν N 1 , N ( ω N 1 ω N ) ) ρ J p Δ x
d ω j d t = ξ ρ J p ( Δ x ) 2 ( ω j 1 2 ω j + ω j + 1 ) + G ρ ( Δ x ) 2 ( φ j 1 2 φ j + φ j + 1 ) , j = 3 , , N 2
Rotation angle equations of discrete shaft units, rotors, and vertical pumps flywheels will be attached to (39)–(41):
d φ i d t = ω i , i = 1 , 2 , , N k
Analysis of Equations (1), (2), (9), (10), (24), and (25) implies presence of an unknown function: V—drive system voltage. It is not easy to find that function. We shall proceed as follows. For the drive system in Figure 1, the first Kirchhoff equation will be written as:
i 20 + n = 1 N i S , n + m = 1 M i S , m = 0 , n = 1 , , N , m = 1 , , M
Then, the initial conditions for the state Equations (1), (2), (9), (10), (24), and (25) will be considered and (43) will be differentiated over time to produce:
d i 20 d t + d d t n = 1 N i S , n + d d t m = 1 M i S , m = 0
By substituting (2), (9), and (24) to (44), the columnar vector of V assembly voltage will be found:
V = { A 22 , 0 + n = 1 N A S , n + m = 1 N A S , m } 1 { A 22 , 0 r 20 i 20 A 21 , 0 ( u 10 r 10 i 10 ) + + n = 1 N [ A S , n r S , n i S , n A S R , n ( u R , n Ω 1 , n Ψ R , n r R , n i R , n ) ] + + m = 1 M [ A S , m ( r S , m i S , m Ω m Ψ S , m ) A S R , m ( u R , m r R , m i R , m ) Ω m i S , m ] }

4. Calculation of the Loading Torques of Electric Pumping Drive

Loading torques in the drive system in Figure 2 apply to vertical pumps. The torques are determined using head H waveforms as a function of efficiency Q, H = f(Q). We analyse two vertical pumps here: one, OB 16-87, loads asynchronous motors, and the other, OB 2-87, loads synchronous motors.
The well-known method employing classic waveforms H = f(Q) for varying rotational speeds of the pump is used to determine the loading torque [30,31].
The following assumptions are made for the purpose of determining a pump’s loading torque. A pumping system, see Figure 2, pumps water from one to the other tank, both situated at the same levels, yet the distance between them is important. The dependence between efficiency and the mains and pump heads as a function of capacity is illustrated in Figure 4.
The waveforms 1–7 are plotted for the speeds: 585, 550, 500, 450, 400, 350, and 300 rpm, respectively. Waveforms 8, 9, and 10 represent the pipeline heads: Hg = 0 m; Hg = 2 m; and Hg = 4 m, and the pump’s capacity is represented with 11.
In the other case under analysis, the level difference between the tanks is taken into account. Where the geometrical head Hg > 0, the problem is solved with numerical methods. Waveforms of the rotational torque across the shaft of pump OB 16-87 working with a steel pipeline of diameter D = 1200 mm and length l = 500 m for Hg = 0 m, 2 m, and 4 m are shown in Figure 5.
An equivalent thickness of the pipeline wall irregularities is assumed as Δe = 1 mm, and total drag coefficients of the pipeline are Σζ = 1.75 [31]. Since flows in such systems can be expressed with Reynolds numbers to the order of 106, a parabolic dependence of the hydraulic resistance within the range of pipeline capacity is assumed.
In the event, the pipeline Hp equation is written as:
H p = H g + S p Q 2
where Sp—hydraulic resistance of the pipeline.
In the simplest case of a pipeline head Hg = 0, dependences of the rotational torque MP across the pump shaft are determined analytically as a function of rotational speed ω:
M P = ρ g S p Q o 3 η o ω o 3 ω 2
where ρ—water density; and Qo, ηo, ωo—capacity, efficiency, and rotational speed of a pump.
Loading torques of asynchronous drives will be based on Figure 5 as follows:
M P 1 = 18.449 ω + 0.691 ω 2
M P 2 = 23.817 ω + 0.730 ω 2
The method of calculating dependences for head of water and loading torques for synchronous drives is the same as for asynchronous motors. In the case of synchronous motors working with an OB 2-87 vertical pump, the waveforms of efficiency, and pipeline’s and pump’s heads as a function of efficiency, are depicted in Figure 6.
The waveforms 1–6 are plotted for speeds of: 585, 550, 500, 450, 400, and 350 rpm, respectively. Waveforms 7 and 8 illustrate the pipeline’s head for Hg = 0 m and Hg = 4 m, respectively, while the pump’s efficiency is represented with 9.
Figure 7 illustrates the dependence of the loading torque across the pump shaft on rotational speed ω at start-up of a synchronous pumping system.
The difference between the tank levels is considered in this case. Figure 7 presents waveforms of rotational torque across OB 2-87 pump shaft working with a steel pipeline of diameter D = 1200 mm and length l = 500 m for heads Hg = 0 m and 4 m. The synchronous machine’s loading torque will be determined using Figure 7 as follows:
M P = 3.958   ω 2 1.856
The following Equations: (1), (2), (9), (10), (21), (24), (25), (34), (35), and (39)–(42) are jointly integrated considering the algebraic dependences: (3)–(8), (11)–(16), (22), (23), (26)–(33), (45), and (48)–(50).

5. Results of Computer Simulation

The computer simulation relates to the drive system whose diagram is shown in Figure 1, Figure 2 and Figure 3. Operation of two deep bar cage asynchronous motors and of one synchronous motor is shown. These are the parameters of the drive system.
  • Power transformer. Type: TM-1600 (1.6 MVA; 35/6.3 kV).
  • Synchronous motor: SDZB 13-52-8, PN = 630 kW, nN = 750 rpm; UN = 6 kV; uf = 42 V, p0 = 4, J1 = 1700 kg·m2, which drives the left end of a distributed mechanical parameter shaft via a reducer gear kT = 750/585 and the first rigid coupling SP1, cf. Figure 2. The rotational torque is then transmitted by a long shaft with the following parameters: L = 4.5 m, D = 0.2 m, and JP = 8.1 1010 N m and the second rigid coupling SP2 to the input shaft of an OB 2-87 vertical pump. These are the parameters of the hydraulic system: nN = 585 rpm, Hmax = 13.6 m. QN = 10,700 m3/h. The pipeline details: pipe diameter D = 1.2 m, length l = 1500 m, and Hg = 4 m.
  • Asynchronous drives with identical motors and motion transmissions but different heads of water lifted by vertical pumps. Motor type: 12-52-8A, PN = 320 kW, UN = 6 kV, IN = 39. Parameters of the long shafts: L = 4.5 m, D = 0.15б, and JP = 8.1 1010 N m. Vertical pump type: OB 16-87, D = 1200 mm L = 500 m, Hg = 0 m, and 4 m, Δe = 1 mm, and Σζ = 1.75.
Electromechanical analysis of transient processes proceeds as follows. The power transformer is connected to an electricity source with a rated voltage supply of 35 kV. The asynchronous drives are disconnected from the drive system; only the synchronous drive is in operation. The synchronous motor is started as the pump valve on the pipeline is closed. The excitation winding is short circuited by a low resistance. Once a sub-synchronous speed is reached, the rotor’s excitation winding is switched to a DC source and the pump valve is opened at t = 30 s. After the motor is in the synchronous state, two asynchronous drives are switched on at t = 40 s. The testing continues until the pumping drive reaches its steady state.
Three experiments are undertaken for the following excitation voltages of the synchronous motor: uf = 12 V; 26 V; and 42 V. This corresponds to the synchronous motor operation in the following three states: asynchronous state with a resistance-inductance loading R-L, synchronous operation with a resistance loading R, cos φ = 1, and synchronous operation with a resistance-capacitance loading R-C relative to the mains, that is, an overcompensation state.
All the mechanical parameters of the three drives are converted by reducer gears to the input vertical pump shaft considering the gear kT = 750/585. The system of ordinary differential equations is integrated by means of the fourth-order Runge–Kutta method.
The subsequent figures show instantaneous voltage, current, and torque waveforms for the synchronous motor’s excitation voltage uf = 12 V—experiment one.
Figure 8 depicts instantaneous phase A voltage of the drive system and phase A current of the power transformer secondary winding for uf = 12 V.
In the steady state, very high inductive currents flow across the transformer winding. The synchronous motor is fed with passive power from the mains due to the reduced excitation voltage of the rotor, which leads to the loss of synchronization and asynchronous operation of the drive. The drive system is in unstable operation, which causes high currents to flow through the secondary winding of the transformer.
Physically it is understandable. During the start-up of asynchronous drives, after the time t = 40 s, there are very high values of the magnetizing current of asynchronous motors. In addition, the synchronous motor consumes reactive power in steady state. It is obvious that this operating state is not normal for the load assembly.
Instantaneous currents across phase A of the armature and across the excitation winding of the synchronous motor are presented in Figure 9.
In the steady state, non-periodic phase currents of enhanced amplitude flow across the armature winding, nearly three times greater than the rated current, which affirms the emergency state of the motor’s operation. The waveform of the motor’s excitation current is similar. At start-up, the excitation current is oscillatory. It should become constant in the motor’s steady state, on the other hand, which can be noted in Figure 9b. The reason for this state is that the synchronous motor is operating in an asynchronous state. Then, a slip occurs, which in turn leads to the induction of an additional alternating current component in the excitation winding. Since high-amplitude oscillatory currents flow across the excitation wiring in the steady state, the machine’s operation is certainly asynchronous. Analysis of Figure 8 and Figure 9 shows that the operation of the electric part of the drive system is incorrect.
We analyse the mechanical part of the three drives. Figure 10 shows a transitional moment of torsion in the middle shaft section and electromagnetic torque of the synchronous drive.
Considering the principles of electromechanical energy conversion in a drive system, a very close connection between electromagnetic and mechanical motions is obvious. Electric drives based on susceptible transmission of motion are known to experience such dangerous phenomena as resonance or resonance-like processes (rumble, vibrations) in which the natural frequencies of the system are the same or close to the forced oscillation of the input signal. In our case of an electric drive containing a long shaft with distributed mechanical parameters, there is a problem of calculating the determined natural frequency of the drive motion transmission. Therefore, we will continue to use the term “range of natural frequencies.” By analysing the electromagnetic torque versus time it is possible to state the occurrence of oscillatory phenomena in the drive transmission of motion.
Electromagnetic torque of the motor in Figure 10b is the input signal transmitted to the mechanical section of the synchronous drive, that is, its motion transmission. For the assumed parameters of the drive system elements, the resonant frequency of the synchronous drive motion transmission is approximately 11 Hz. For these oscillations of the electromagnetic torque, resonant effects will arise in the synchronous drive, which can be noted in Figure 10a. In the time intervals t ∈ (20; 25) s and t ∈ (45; 50), synchronous drive is in resonance.
Instantaneous angular velocity of the synchronous motor rotor and the waveforms of the moment of torsion in the middle shaft section and of the electromagnetic torque of the synchronous drive are presented in Figure 11.
Instantaneous angular velocity of the synchronous motor is depicted in Figure 11a. This shows that, in the steady state, the drive operates at an asynchronous angular velocity and evidence the motor is under excited. In Figure 11b, the waveforms of moment of torsion in the middle shaft section are blue and of the electromagnetic torque of the synchronous drive are red. The resonance epicenter is within the time range t ∈ (46.8; 48.3) s. At t = 47.55 s, amplitude frequencies of the moment and the torque are virtually identical. This means the resonance reaches its maximum at that time.
The waveforms of the synchronous drive’s moment of torsion in three-time intervals: (a)—t ∈ (14; 15) s, (b)—t ∈ (26; 46) s, and (c)—t ∈ (64; 65) s are illustrated in Figure 12.
Figure 12a demonstrates the moment’s waveform is nearly harmonic and its frequency is ca. 38 Hz, which is substantially more than the synchronous drive’s own resonant frequency of about 11 Hz. Resonance is obviously absent from the system. Figure 12b shows an essentially different situation. In the time interval t ∈ (26; 46) s, the moment of torsion’s waveform is highly non-linear. Fourier series can be used to determine the moment’s frequency. Frequencies will then vary for diverse harmonics of the moment. The first harmonic is the most important in this case. Its frequencies are of course several times lower than resonant frequency of the drive’s motion transmission. There is no resonance here, either. Steady low-amplitude resonant processes can be observed in Figure 12c. The amplitude of resonant oscillation is now around eight lower than the amplitude in the range t ∈ (46.8; 48.3) s, see Figure 10a. In this interval, oscillation frequencies cohere at higher harmonics of the moment of torsion distribution over the Fourier series.
At the time of the asynchronous start-up of synchronous drive, the frequency of electromagnetic torque varies over a rather wide range. Oscillations occur in these states too, occasionally very non-linear, which give rise to resonance at the higher harmonics of electromagnetic torque. The only advantage is that the oscillation amplitudes of electromagnetic torque are significantly lower than the amplitude of the first harmonic.
The waveforms of phase A current and voltage of the armature winding in the interval t ∈ (69.88; 70) s and of phase A current of the asynchronous motor’s stator winding in the interval t ∈ (40; 46) s are presented in Figure 13.
Figure 13a shows the operating state of the synchronous drive from the perspective of the network supplying power to the drive system. The waveform of phase A voltage is marked blue and of phase A current is red. The phase of the voltage leads the current by a significant angle of ca. φ = 50°. This corroborates the conclusion that the synchronous drive operates in an asynchronous state of under excitation.
In such a state, quite permanent transient processes are visible, which are related to the resistive–inductive load of the network.
The moment of torsion across the motion transmission shaft and electromagnetic moment in the time interval t ∈ (40; 46) s are illustrated in Figure 14.
Figure 14 demonstrates electromechanical processes in the synchronous drive, the key part of the drive system, which affects the processes in both the asynchronous drives. The waveforms of current, torque, and moment contain additional oscillations. The electromagnetic torque and the moment of torsion in the steady state are not constant, either, but oscillatory in nature.
The successive figures include instantaneous waveforms of voltages, currents, torques, and moments for the synchronous motor’s excitation voltage of uf = 26 V—the second experiment.
The instantaneous phase A voltage of the drive system and phase A current of the secondary power transformer wiring are presented in Figure 15.
Instantaneous waveforms of phase A current of the synchronous motor’s armature wiring and excitation current are depicted in Figure 16.
A comparative analysis with the results of the first experiment (cf. Figure 8 and Figure 9) suggests a totally different operating condition of the drive system in experiment two. The synchronous drive is in a normal operating state now. The excitation voltage is selected in such a way that cos φ = 1, which means the synchronous motor only consumes the active power. The stator’s current in its steady state is sinusoid, see Figure 9a and Figure 16a. The excitation current in the steady state is constant, see Figure 9b and Figure 16b. The drive in the steady state operates at a synchronous speed, cf. Figure 11a and Figure 17a.
Figure 17a shows the instantaneous angular velocity of the synchronous drive rotor. Compared with Figure 11a, the electric drive in the steady process rotates at a synchronous speed.
The instantaneous moment of torsion in the middle shaft section and electromagnetic torque of the synchronous drive are illustrated in Figure 17b,c, which need to be compared with Figure 10. The moment of torsion across the drive shaft is exclusive of the second resonance once the excitation wiring is connected to a DC source, as the motor has become synchronous and the electromagnetic torque does not include oscillatory motions in the steady state. In addition, the moment of torsion in the steady state is constant, since the stator current is exclusive of higher harmonics and resonant processes are absent from the synchronous drive. The drive becoming synchronous substantially reduces the amplitudes of transitional shock torques and currents. It is also dependent on the motor’s excitation current.
Phase A voltage and current of the armature winding and electromagnetic torque of the first asynchronous motor are shown in Figure 18.
The voltage is blue and the current is red in the time interval t ∈ (54.88; 55) s in Figure 18a. A synchronous motor is a purely resistive load of the network, i.e., φ = 0°. We note again the stator current is sinusoid, as distinct from the first experiment, cf. Figure 13a. From the point of view of applied electrical engineering, it is related to the type of the assembly load. In the case shown in Figure 13a, the under excited synchronous motor is a resistive–inductive load on the network. According to the characteristics of the synchronous motor, the resistive–inductive load reduces the voltage of the assembly the most. However, in the case presented in Figure 18a, the situation is different—the motor is a purely resistive network load. In this case, the time constant is minimal and processes reach a steady state faster.
The instantaneous electromagnetic torque of the first asynchronous motor in the time interval t ∈ (40; 46) s is shown in Figure 18b. Compared with Figure 14b, adverse oscillatory processes are absent from asynchronous drives associated with the asynchronous operation of the synchronous drive in experiment one. Thus, both the asynchronous drives work in the normal state. Note the synchronous drives play the main role in the system under discussion. Therefore, we compare operation of the drive system to diverse operating conditions of the synchronous drive in this study. For instance, [21,24,28] address different operating conditions of susceptible motion transmission asynchronous drives that include long shafts.
The following graphs illustrate instantaneous waveforms of voltages, currents, moments, and torques for the synchronous motor’s rated excitation voltage of uf = 42 V—experiment three.
The instantaneous phase A current across the secondary power transformer winding is shown in Figure 19a.
At uf = 42 V, the synchronous drive works at rated excitation. The supply network is loaded with a volumetric-resistance receiver. Compared with Figure 15b, the current across the synchronous motor armature’s winding is larger in steady states. This is an example of passive power compensation by increasing the excitation current of a synchronous machine. In the final analysis, raising motor excitation current increases the voltage supplied to the drive system.
The instantaneous excitation current of the synchronous motor is shown in Figure 19b. Its waveform is similar to the one in Figure 16b. In the steady state of the second experiment (for uf = 26 V), on the other hand, the excitation current is around if = 100 A, and for experiment three (uf = 42 V), approximately if =180 A.
The instantaneous waveforms of phase A voltage and current across the motor armature winding are presented in Figure 20a.
The voltage is blue and the current red in the time interval t ∈ (54.88; 55) s in Figure 20a. From the perspective of the power network, the synchronous drive acts as a resistant-volumetric load with a phase shift angle of circa ϕ = −40°. The armature current is ahead of the voltage supplied to the driving motors. Instantaneous angular velocities of the first asynchronous motor’s rotor are depicted in Figure 20b, marked blue for the second experiment and red for the third.
Analysis of both of the waveforms proves operation of the drive system in the second and third experiments is correct. Owing to stabilisation, or increasing of the drive system’s voltage, the start-up time of the asynchronous drive is of course around 0.3 s shorter in the third than in the second experiment. This time difference is negligible for the practical purposes of both of the pumps’ operation. As far as the power network is concerned, the synchronous motor in the steady state generates passive power.

6. Conclusions

It is in some cases reasonable to introduce a synchronous drive to stabilise voltage in drive systems containing asynchronous drives. It helps to arrive at an absolutely rigid mechanical characteristic of a drive in the static sate, on the one hand, and to avoid adding compensation devices, for instance, static capacitor batteries, to a drive system.
Operation of synchronous drives in asynchronous states is a quite dangerous phenomenon in drive systems. In steady states, this increases currents several times compared with the rated current in the motor’s armature winding. Furthermore, it causes induction of oscillatory current waveforms in the excitation winding of a motor. Therefore, the excitation current of a motor should be calculated and adapted to all operating states of a synchronous drive as part of an electric power system.
Long elastic elements should ultimately be analysed as systems of distributed mechanical parameters for the purposes of modelling transient processes in pumping electric drives that contain complex motion transmissions including long shafts. This enables to address some latent motions that cannot always be determined, e.g., as part of the classic dual-mass interpretation of long motion transmission shafts.
During asynchronous start-up of synchronous drives with large masses of inertia, the electromagnetic torque of a motor varies its oscillation frequencies over a quite broad range. In the event, there is a risk of the motion transmission entering resonant or nearly resonant states not only for the first but also for the higher harmonics of the torque.
The results of computer simulation of transient electromechanical processes show:
  • In resonant states, the moment of torsion in a drive’s motion transmission occasionally rises as much as several dozen times relative to the nominal load torque of a drive. If the susceptibility of long elastic elements of electric drives’ motion transmission are neglected and resonant motions are not addressed, a motion transmission of a drive may be wrongly selected at the stage of design.
  • In asynchronous drives, the total inertia of drive motion transmission is substantially lower than in equivalent synchronous drives. In asynchronous start-up of a synchronous drive, the system is then likely to enter resonant or nearly resonant states over a quite wide time interval.
  • In asynchronous operating states of synchronous drives, given low excitation currents of a motor, a system is very likely to enter a second wave of resonant oscillations whose amplitude is greater than the amplitude of oscillations in the first resonant wave.
  • In asynchronous operating states of synchronous drives, contrary to synchronous states in steady states, the electromagnetic torque includes complicated oscillatory waveforms. This causes a drive to enter resonant states at the higher harmonics of the torque. The amplitude of the moment of torsion’s oscillations in the motion transmission shaft is substantially greater than the oscillation amplitude for the first harmonic of the Fourier series, which is a good development.
  • Increasing the excitation current of a synchronous drive motor stabilises the drive system’s voltage and changes the nature of the synchronous drive’s load in the drive system to more resistance–volumetric loading.
  • Depending on the values of excitation current, start-up times of all the drives in an electric power drive system can be regulated.
The proposed method of modelling transient processes in drive systems that contain electric load assemblies enables the analysis of transient processes not only in pumping systems, but practically in all electric drives. This makes it possible to extend the scope of our research in other branches of the economy. In our next plans, we will put into practice the simulation results that we obtained in this paper.

Author Contributions

Conceptualization, A.C.; methodology, A.C., M.L. and A.S.; software, A.C. and M.L.; validation, A.C., M.L. and A.S.; formal analysis, A.C., M.L. and A.S.; investigation, A.C. and M.L.; writing—original draft preparation, A.C., M.L. and A.S.; writing—review and editing, A.C., M.L. and A.S.; visualization, M.L. and A.S.; supervision, A.C., M.L. and A.S.; project administration, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Parameters and variables
ασ1, ασ2Matrices of reverse transformer leakage inductances
ρMatrix of reverse differential magnetising inductances
τMatrix of reverse magnetising inductances
γAngle of rotor rotation
ηoEfficiency pump
ξCoefficient of shaft internal leakage
ΠPark matrix of rectangular transformations
Π1Matrix of angle transformations
ΩMatrix of angular rotational speeds of the synchronous machine in rectangular coordinates
Ω1Matrix of angular rotational speeds of asynchronous motor in angle coordinates
μMagnetic permeability of conductor material in a rotor deep bar
φ(x,t)Angle of shaft rotation
ωAngular velocity of the rotor
ωNAngular velocity of the vertical pump rotor
ωoRotational speed of a pump
ω1Angular velocity angle of the drive motor rotor
φNRotation angle of the vertical pump rotor
φ1Rotation angle of the drive motor rotor
ψm(im)Working flux of machine magnetisation
ψm,jWorking associated transformer fluxes, to be determined by means of the power transformer magnetising curve
ΨRColumnar vector of main associated fluxes of the rotor
ΔxDiscretisation step of the long shaft equation
ΔzDiscretisation step
ν1,2Leakage coefficients of the first coupling
νN−1,NLeakage coefficients of the second coupling
aDeep bar width
A0Matrix of reverse transformer inductances
AMAsynchronous drives
C1,2Rigidity coefficients of the first coupling
CN−1,NRigidity coefficients of the second coupling
dIndex related to parameters in the machine’s longitudinal axis
DIndex related to parameters of the damping winding in the longitudinal axis d
EVector of the electric field intensity
E(0)Electric intensity on the deep bar surface
e(t)Columnar vector of phase electromotive forces of a rigid electric power grid
fIndex related to parameters of the excitation winding
GTransverse torsion module
hBar depth in relation to the coordinate z
HVector of the magnetic field intensity
i10Columnar vectors of phase currents of the transformer’s primary windings respectively
i20Columnar vectors of phase currents of the transformer’s secondary windings, respectively
imColumnar vector of transformer magnetising currents
jIndices of three-phase system phases
JEMMoment of inertia of the motor’s drive rotor
JPMoment of inertia of the hydraulic pump rotor recalculated considering the reducer gear
JpPolar moment of shaft inertia
kNumber of asynchronous motor in the drive system
kTValue of the reducer’s gear
ku, kiVoltage and current reduces coefficients of machine transformation
lLength of a deep bar
LmInductance of machine magnetisation
MNumber of synchronous drives in the drive system
mNumber of synchronous drives in the drive system
MEMElectromagnetic torque of the drive motor
MPPump’s hydraulic torque
NNumber of asynchronous drives in the drive system
NNumber of discretisation nodes
p0Number of machine pole pairs
qIndex related to parameters in the machine’s transverse axis
QIndex related to parameters of the damping winding in the transverse axis q
QoCapacity pump
r0Matrix of winding resistances
SMSynchronous drives
SpHydraulic resistance of the pipeline
SP1Flexible coupling between the rotor and motion transmission long shaft (the first coupling)
SP2Flexible coupling between the motion transmission long shaft and the vertical pump’s Input rotor (the second coupling)
TNumber of discretisation nodes relative to the coordinate z
u10Columnar vector of phase voltages supplied to the power transformer’s primary winding
u20Columnar vector of phase voltages supplied to the power transformer’s secondary winding
uRColumnar vector of voltages across the rotor cage bars
VColumnar vector of phase drive system voltages
VDrive system voltage
xCurrent coordinate along the shaft
zCoordinate oriented into a cylindrical system deep bar

References

  1. Martin-Candilejo, A.; Martin-Carrasco, F.J.; Santillán, D. How to Select the Number of Active Pumps during the Operation of a Pumping Station: The Convex Hyperbola Charts. Water 2021, 13, 1474. [Google Scholar] [CrossRef]
  2. Shokoohi, M.; Tabesh, M.; Nazif, S.; Dini, M. Water quality based multi-objective optimal design of water distribution systems. Water Resour. Manag. 2017, 31, 93–108. [Google Scholar] [CrossRef]
  3. Kurek, W.; Ostfeld, A. Multi-objective optimization of water quality, pumps operation, and storage sizing of water distribution systems. J. Environ. Manag. 2013, 115, 189–197. [Google Scholar] [CrossRef] [PubMed]
  4. Salomons, E.; Housh, M.; Asce, M. A Practical Optimization Scheme for Real-Time Operation of Water Distribution Systems. J. Water Resour. Plan. Manag. 2020, 146, 04020016. [Google Scholar] [CrossRef]
  5. Palkin, G.; Suvorov, I. Simulation Modeling of First Rise Section of Water Supply System with Installed Complex of Automatic Pump Performance Control. Machines 2021, 9, 63. [Google Scholar] [CrossRef]
  6. Zhou, X.; Wang, C.; Ma, Y. Vector Speed Regulation of an Asynchronous Motor Based on Improved First-Order Linear Active Disturbance Rejection Technology. Energies 2020, 13, 2168. [Google Scholar] [CrossRef]
  7. Hannan, M.A.; Ali, J.A.; Mohamed, A.; Hussain, A. Optimization techniques to enhance the performance of induction motor drives: A review. Renew. Sustain. Energy Rev. 2018, 81, 1611–1626. [Google Scholar] [CrossRef]
  8. Lei, G.; Zhu, J.; Guo, Y.; Liu, C.; Ma, B. A Review of Design Optimization Methods for Electrical Machines. Energies 2017, 10, 1962. [Google Scholar] [CrossRef] [Green Version]
  9. Chaban, A. Hamilton-Ostrogradski Principle in Electromechanical Systems; Soroki: Lviv, Ukraine, 2015; p. 488. [Google Scholar]
  10. Lis, M.; Chaban, A.; Szafraniec, A.; Figura, R.; Levoniuk, V. Mathematical model of a part of an opened extra-high voltage electrical grid. In E3S Web of Conferences; EDP Sciences: Podlesice, Poland, 2019; Volume 84, p. 02005. [Google Scholar]
  11. Czaban, A.; Rusek, A.; Lis, M. The approach based on variation principles for mathematical modeling of asymmetrical states in a power transformer. Przegląd Elektrotechniczny 2012, 88, 240–242. [Google Scholar]
  12. Lee, Y.; Song, H. A Reactive Power Compensation Strategy for Voltage Stability Challenges in the Korean Power System with Dynamic Loads. Sustainability 2019, 11, 326. [Google Scholar] [CrossRef] [Green Version]
  13. Buitrago-Velandia, A.F.; Montoya, O.D.; Gil-González, W. Dynamic Reactive Power Compensation in Power Systems through the Optimal Siting and Sizing of Photovoltaic Sources. Resources 2021, 10, 47. [Google Scholar] [CrossRef]
  14. Vlahinić, S.; Franković, D.; Komen, V.; Antonić, A. Reactive Power Compensation with PV Inverters for System Loss Reduction. Energies 2019, 12, 4062. [Google Scholar] [CrossRef] [Green Version]
  15. Liu, H.; Wang, Z.; Yang, J.; Li, B.; Ren, A. Circuit Breaker Rate-of-Rise Recovery Voltage in Ultra-High Voltage Lines with Hybrid Reactive Power Compensation. Energies 2018, 11, 100. [Google Scholar] [CrossRef] [Green Version]
  16. Doğru, U.; Özçelik, N.; Gedik, H.; İmeryüz, M.; Ergene, L. Numerical and Experimental Comparison of TLA Synchronous Reluctance Motor and Induction Motor. In Proceedings of the 2016 IEEE International Power Electronics and Motion Control Conference (PEMC 2016), Varna, Bulgaria, 25–28 September 2016; pp. 619–624. [Google Scholar]
  17. Lang, K.; Muetze, A.; Bauer, R.; Pircher, S. Comparison of induction and synchronous reluctance machine based actuators for elevated temperature environments. IEEE Trans. Energy Convers. 2016, 31, 1012–1022. [Google Scholar] [CrossRef]
  18. Kazakbaev, V.; Prakht, V.; Dmitrievskii, V.; Ibrahim, M.N.; Oshurbekov, S.; Sarapulov, S. Efficiency Analysis of Low Electric Power Drives Employing Induction and Synchronous Reluctance Motors in Pump Applications. Energies 2019, 12, 1144. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, W.; Liang, J.; Shi, T. Speed Synchronous Control of Multiple Permanent Magnet Synchronous Motors Based on an Improved Cross-Coupling Structure. Energies 2018, 11, 282. [Google Scholar] [CrossRef] [Green Version]
  20. Hussein, I.; Al-Hamouz, Z.; Abido, M.A.; Milhem, A. On the Mathematical Modeling of Line-Start Permanent Magnet Synchronous Motors under Static Eccentricity. Energies 2018, 11, 197. [Google Scholar] [CrossRef]
  21. Czaban, A.; Lis, M.; Szafraniec, A.; Chrzan, M.; Levoniuk, V. Interdisciplinary modelling of transient processes in local electric power systems including long supply lines of distributed parameters. In Proceedings of the 2018 Applications of Electromagnetics in Modern Techniques and Medicine (PTZE), Raclawice, Poland, 9–12 September 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 105–108. [Google Scholar]
  22. Ekemb, G.; Slaoui-Hasnaoui, F.; Song-Manguelle, J.; Lingom, P.M.; Fofana, I. Instantaneous Electromagnetic Torque Components in Synchronous Motors Fed by Load-Commutated Inverters. Energies 2021, 14, 3223. [Google Scholar] [CrossRef]
  23. Popenda, A.; Lis, M.; Nowak, M.; Blecharz, K. Mathematical Modelling of Drive System with an Elastic Coupling Based on Formal Analogy between the Transmission Shaft and the Electric Transmission Line. Energies 2020, 13, 1181. [Google Scholar] [CrossRef] [Green Version]
  24. Chaban, A.; Lis, M.; Szafraniec, A.; Jedynak, R. Application of Genetic Algorithm Elements to Modelling of Rotation Processes in Motion Transmission Including a Long Shaft. Energies 2020, 14, 115. [Google Scholar] [CrossRef]
  25. Pukach, P.Y.; Kuzio, I.V.; Nytrebych, Z.M.; Ilkiv, V.S. Analytical methods for determining the effect of the dynamic process on the nonlinear flexural vibrations and the strength of compressed shaft. Naukovyi Visnyk Natsionalnoho Hirnychoho Universytetu 2017, 5, 69–76. [Google Scholar]
  26. Popenda, A.; Szafraniec, A.; Chaban, A. Dynamics of Electromechanical Systems Containing Long Elastic Couplings and Safety of Their Operation. Energies 2021, 14, 7882. [Google Scholar] [CrossRef]
  27. Lis, M. Mathematical Modeling of Transient Processes in Electric Drive Systems with Complex Motion Transmission; Publishing House of the Częstochowa University of Technology: Czestochowa, Poland, 2013; p. 258. [Google Scholar]
  28. Lozynskyy, A.; Chaban, A.; Perzyński, T.; Szafraniec, A.; Kasha, L. Application of Fractional-Order Calculus to Improve the Mathematical Model of a Two-Mass System with a Long Shaft. Energies 2021, 14, 1854. [Google Scholar] [CrossRef]
  29. Popenda, A. Mathematical modelling of transmission shafts based on electrical and mechanical similarities. Przegląd Elektrotechniczny 2019, 12, 196–199. [Google Scholar] [CrossRef]
  30. Łukasik, Z.; Czaban, A.; Szafraniec, A.; Żuk, W. The mathematical model of the drive system with asynchronous motor and vertical pump. Przegląd Elektrotechniczny 2018, 1, 133–138. [Google Scholar] [CrossRef] [Green Version]
  31. Szafraniec, A. Mathematical modeling of oscillating processes in the electro-hydraulic drives with susceptible motion transmission. Przegląd Elektrotechniczny 2017, 12, 167–170. [Google Scholar]
  32. Popenda, A. Model-simulation investigations of induction motor with the consideration of skin effect in rotor bars. Przegląd Elektrotechniczny 2012, 88, 29–31. [Google Scholar]
Figure 1. Circuit diagram of the drive system.
Figure 1. Circuit diagram of the drive system.
Energies 15 00811 g001
Figure 2. Basic scheme of the pumping system.
Figure 2. Basic scheme of the pumping system.
Energies 15 00811 g002
Figure 3. Calculation scheme of the shaft continuum.
Figure 3. Calculation scheme of the shaft continuum.
Energies 15 00811 g003
Figure 4. Efficiency and head waveforms as a function of the vertical pump OB 16-87 capacity.
Figure 4. Efficiency and head waveforms as a function of the vertical pump OB 16-87 capacity.
Energies 15 00811 g004
Figure 5. Waveforms of the rotational torque across the pump shaft as a function of angular velocity for: 1—Hg = 0 m; 2—Hg = 2 m; and 3—Hg = 4 m.
Figure 5. Waveforms of the rotational torque across the pump shaft as a function of angular velocity for: 1—Hg = 0 m; 2—Hg = 2 m; and 3—Hg = 4 m.
Energies 15 00811 g005
Figure 6. Waveforms of efficiency and head of water as a function of OB 2-87 vertical pump’s capacity.
Figure 6. Waveforms of efficiency and head of water as a function of OB 2-87 vertical pump’s capacity.
Energies 15 00811 g006
Figure 7. Waveforms of rotational torque across the pump shaft as a function of angular velocity for: 1 − Hg = 0 m; and 2 − Hg = 4 m.
Figure 7. Waveforms of rotational torque across the pump shaft as a function of angular velocity for: 1 − Hg = 0 m; and 2 − Hg = 4 m.
Energies 15 00811 g007
Figure 8. Instantaneous waveforms of: (a) phase A voltage of the drive system; and (b) phase A current of the power transformer secondary winding.
Figure 8. Instantaneous waveforms of: (a) phase A voltage of the drive system; and (b) phase A current of the power transformer secondary winding.
Energies 15 00811 g008
Figure 9. Instantaneous current waveforms of: (a) phase A of the synchronous motor armature; and (b) synchronous motor excitation.
Figure 9. Instantaneous current waveforms of: (a) phase A of the synchronous motor armature; and (b) synchronous motor excitation.
Energies 15 00811 g009
Figure 10. Instantaneous waveforms of: (a) moment of torsion in the middle section of the synchronous drive motion transmission long shaft; and (b) electromagnetic torque of the synchronous motor.
Figure 10. Instantaneous waveforms of: (a) moment of torsion in the middle section of the synchronous drive motion transmission long shaft; and (b) electromagnetic torque of the synchronous motor.
Energies 15 00811 g010
Figure 11. Instantaneous: (a) angular velocity of the synchronous drive rotor; and (b) moment of torsion in the middle shaft section (blue) and electromagnetic torque of the synchronous drive (red).
Figure 11. Instantaneous: (a) angular velocity of the synchronous drive rotor; and (b) moment of torsion in the middle shaft section (blue) and electromagnetic torque of the synchronous drive (red).
Energies 15 00811 g011
Figure 12. Instantaneous waveforms of the synchronous drive’s moment of torsion for the following time intervals: (a) t ∈ (14; 15) s; (b) t ∈ (26; 46) s; and (c) t ∈ (64; 65) s.
Figure 12. Instantaneous waveforms of the synchronous drive’s moment of torsion for the following time intervals: (a) t ∈ (14; 15) s; (b) t ∈ (26; 46) s; and (c) t ∈ (64; 65) s.
Energies 15 00811 g012
Figure 13. Instantaneous waveforms of: (a) phase A current (red) and voltage (blue) of the armature winding; and (b) phase A current of the asynchronous motor’s stator winding.
Figure 13. Instantaneous waveforms of: (a) phase A current (red) and voltage (blue) of the armature winding; and (b) phase A current of the asynchronous motor’s stator winding.
Energies 15 00811 g013
Figure 14. Instantaneous waveforms of the moment and torque in the interval t ∈ (40; 46) s: (a) torsion in the middle section of the motion transmission long shaft in the first drive; and (b) electromagnetic moment of the first asynchronous motor.
Figure 14. Instantaneous waveforms of the moment and torque in the interval t ∈ (40; 46) s: (a) torsion in the middle section of the motion transmission long shaft in the first drive; and (b) electromagnetic moment of the first asynchronous motor.
Energies 15 00811 g014
Figure 15. Instantaneous waveforms of: (a) phase A voltage of the drive system; and (b) phase A current in the secondary power transformer wiring.
Figure 15. Instantaneous waveforms of: (a) phase A voltage of the drive system; and (b) phase A current in the secondary power transformer wiring.
Energies 15 00811 g015
Figure 16. Instantaneous waveforms of: (a) A current of the synchronous motor’s armature wiring; and (b) excitation current of the synchronous motor.
Figure 16. Instantaneous waveforms of: (a) A current of the synchronous motor’s armature wiring; and (b) excitation current of the synchronous motor.
Energies 15 00811 g016
Figure 17. Instantaneous waveforms of: (a) the angular velocity of the synchronous drive rotor; (b) moment of torsion in the middle section of the synchronous drive’s long motion transmission shaft; and (c) electromagnetic torque of the synchronous motor.
Figure 17. Instantaneous waveforms of: (a) the angular velocity of the synchronous drive rotor; (b) moment of torsion in the middle section of the synchronous drive’s long motion transmission shaft; and (c) electromagnetic torque of the synchronous motor.
Energies 15 00811 g017
Figure 18. Instantaneous waveforms of: (a) phase A voltage (blue) and current (red) of the armature winding; and (b) electromagnetic torque of the first asynchronous motor.
Figure 18. Instantaneous waveforms of: (a) phase A voltage (blue) and current (red) of the armature winding; and (b) electromagnetic torque of the first asynchronous motor.
Energies 15 00811 g018
Figure 19. Instantaneous waveforms of: (a) phase A current across the secondary power transformer winding; and (b) excitation current of the synchronous motor.
Figure 19. Instantaneous waveforms of: (a) phase A current across the secondary power transformer winding; and (b) excitation current of the synchronous motor.
Energies 15 00811 g019
Figure 20. Instantaneous waveforms of: (a) phase A voltage (blue) and current (red) across the motor armature winding; and (b) angular velocities of the first asynchronous rotor (blue—the second experiment, red—the third experiment).
Figure 20. Instantaneous waveforms of: (a) phase A voltage (blue) and current (red) across the motor armature winding; and (b) angular velocities of the first asynchronous rotor (blue—the second experiment, red—the third experiment).
Energies 15 00811 g020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chaban, A.; Lis, M.; Szafraniec, A. Voltage Stabilisation of a Drive System Including a Power Transformer and Asynchronous and Synchronous Motors of Susceptible Motion Transmission. Energies 2022, 15, 811. https://doi.org/10.3390/en15030811

AMA Style

Chaban A, Lis M, Szafraniec A. Voltage Stabilisation of a Drive System Including a Power Transformer and Asynchronous and Synchronous Motors of Susceptible Motion Transmission. Energies. 2022; 15(3):811. https://doi.org/10.3390/en15030811

Chicago/Turabian Style

Chaban, Andriy, Marek Lis, and Andrzej Szafraniec. 2022. "Voltage Stabilisation of a Drive System Including a Power Transformer and Asynchronous and Synchronous Motors of Susceptible Motion Transmission" Energies 15, no. 3: 811. https://doi.org/10.3390/en15030811

APA Style

Chaban, A., Lis, M., & Szafraniec, A. (2022). Voltage Stabilisation of a Drive System Including a Power Transformer and Asynchronous and Synchronous Motors of Susceptible Motion Transmission. Energies, 15(3), 811. https://doi.org/10.3390/en15030811

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