CC BY
6
Eastern-European Journal of Enterprise Technologies ISSN 1729-3774
H I-
-□ □-
Для аналтичного анализу нелтшних електричних кш електротехнiчних систем був використаний метод малого параметра, реалгзований в частотны область. Для мож-ливостъ реалгзацп розрахунтв в частотны област1 вико-ристовувався автоматизований метод формування орто-гональних гармоншних складових електричних величин на базi алгоритму дискретног згортки. Характеристика нелтшного елемента була представлена полтомгальною функщею третього ступеня. Показано, що застосування методу малого параметра з його реалгзащею в частот-ног област1 дозволяв спростити процес аналгзу електричних кш з нелтшними елементами в аналтичному за раху-нок автоматизацп розрахункв у математичному пакетг. Анаттичт i чисельтрозрахунки кола з активно-тдуктив-ним навантаженням продемонстрували достатню точ-тсть запропонованого методу, вiдносна похибка за основною гармошкою струму не перевищила 6 %. Проведений порiвняльний аналгз запропонованого методу малого параметра з класичним методом малого параметра на приклад розрахунку електричного кола з ЯЬ навантаженням показав, щорозроблений метод забезпечув бЫьшу адекват-тсть результатов та вищу точтсть розрахункв у порiв-няннг Ь кнуючим. Вiдносна похибка за амплтудою пер-шог та третьог гармонж струму не перевищуе 2,5 % та за фазою - 1,04210-3 %. Метод чисельного структурного моделювання застосовувався для визначення еталонних значень струму дослъджуваного кола. Результаты роботи можуть бути використан при розрахунках електротех-нхчних пристрогв, що мктять напiвпровiдниковi елемен-ти i електричн апарати з нелтшними характеристиками. Також отримат результати дозволять удосконалити процеси активног компенсаци вищих гармонж струму в електричних мережах з нелтшним навантаженням, що мгстить, а також розробляти засоби пасивног компенсаци
Ключовi слова: нетнтна система, електричне коло, аналiз, метод малого параметра, частотна область,
автоматизований алгоритм -□ □-
UDC 621.311:621.314
|DOI: 10.15587/1729-4061.2018.136398|
IMPROVEMENT OF THE METHOD FOR ANALYSIS OF NONLINEAR ELECTROTECHNICAL SYSTEMS BASED ON THE SMALL PARAMETER METHOD
M. Maliakova
PhD, Senior Lecturer Department of Electric Machines and Devices Kremenchuk Mykhailo Ostrohradskyi National University Pershotravneva str., 20, Kremenchuk, Ukraine, 39600 E-mail: mariia.maliakova@gmail.com
1. Introduction
At present, analysis of nonlinear electrical circuits (NEC) is one of the most important tasks of modern theoretical electrical engineering [1, 2]. NEC are an integral part of most electrotechnical and electronic devices. The need for complicating the power semiconductor converters and accounting of new nonlinear effects is becoming good stimuli for the development of the theory of calculation and analysis of NEC. The latter is used increasingly often, both in traditional areas of analysis and synthesis [3, 4], and in new directions of modeling, identification, etc. [2, 5]. Analysis of NEC often comes down to determining the currents by known parameters of a circuit [3]. However, the process of analysis of the researched nonlinear systems is accompanied by specific difficulties, among which we can separate the following:
- impossibility of using well-developed methods for analysis of linear circuits to solve nonlinear ones, due to emerging of phenomena not inherent in linear circuits [4, 5];
- necessity of an accurate approximation of a nonlinear characteristic to maximize the accuracy of calculations [6].
That is why researchers have to resort in practice to a combination of several methods, such as the wavelet transformations and iterative methods [7] of modeling the NEC operation, using a large number of methods of numerical calculation with a complicated set of initial conditions [8]. This in turn leads to a complication of calculation dependences and provokes an increase in calculation errors [3, 5, 6-8].
Today, the task of studying the processes occurring in circuits with semiconductor converters, which are the elements of a clearly pronounced nonlinearity, remains relevant. These are nonlinear circuits of power systems [9, 10], nonlinear circuits of power active filters [11], control systems, which are based on the p-q theory [12] and cross-vector [13] theories of instantaneous power. The main issues that are tackled by researchers in the analysis of such systems include: determining the patterns and an increase in effectiveness of
©
energy conversion [11-14, 16-18], determining the nature of occurrence of reactive power [12, 15, 20], and development of effective methods for calculation of currents of researched circuits [21].
It should also be noted that the relevant problem today is selection of tools [19] of active and passive filtration [22, 23]. But the correct choice of the tools for passive filtration requires the analysis of indicators of operation of a connected nonlinear load, which is possible using the analytical methods of calculation of electrical circuits [21].
Given the above, development of the analytical method for calculation that is effective and simple in practical terms is reasonable. The developed method should provide an opportunity to conduct a qualitative analysis of a nonlinear electrotechnical system and will also make it possible to conduct an analysis of the mechanisms of the current components formation.
2. Literature review and problem statement
As it is known, the methods for analysis of nonlinear electrical circuits can be divided into analytical and numerical.
The use of the numerical methods greatly facilitates the task for researchers because they make it possible with a relative ease of implementation to obtain information about almost any processes in the NEC. Due to the rapid development of computer technology and existence of a large number of specialized software products, such as SAPWIN [24, 25], SPICE [26], SEQUEL [27], MATLAB [28-30], which allow implementation of calculations and numerical modeling, numerical methods the most frequently used. However, it should be noted that these methods have the following shortcomings: numerical solution of equations with significantly pronounced nonlinearity requires a large number of iterations, accompanied by considerable time consumption. Another drawback is the fact the accuracy of the calculation significantly depends on correctness of the chosen initial approximations and configuration of a scheme, as well as on a method and a pitch of integration.
In contrast to numerical methods, analytical methods have a significant advantage. They allow us to conduct a qualitative analysis and give an adequate evaluation to the processes in the NEC and estimate the influence of all factors in the course of operation of a system, their contribution and importance [1, 3, 5, 9, 10, 21].
As it is known, the work of circuits is described by nonlinear differential equations, the solution of which allows determining the currents that flow in a circle.
Based on the analysis of the literary sources, to solve such problems, researchers often choose the small parameter method (SPM) [3, 5, 9, 10, 21], which belongs to the perturbation methods, as the most efficient, universal, and simple to implement. SPM is one of the most effective tools of modern applied mathematics, which is widely used in mechanics, physics, and other sciences that operate with differential equations. SPM allows obtaining analytical solutions of complex linear and nonlinear boundary value problems both for ordinary differential equations and for equations in partial ones. It can be used separately or in a combination with other methods [3, 5].
The main advantage of SPM is its effectiveness in solving complicated differential equations that describe the pro-
cesses in NEC. That is why in this paper it is proposed to use SPM to determine a harmonic composition of current of nonlinear load, connected to the electrical supply system. This will make it possible to select correctly the parameters of a linear reactor to reduce the impact of higher harmonics in the current signal.
According to SPM [3, 9, 10], a non-linear differential equation is solved in the time domain by successive approximations, representing the sought parameter x as a series by powers of some coefficient which is called the small parameter:
x = x0 +11X1 + 1 x2 + ... + l"Xn, (1)
where x0 is the solution to an equation of zero approximation (the latter is obtained from the original, assuming all nonlinear terms are absent in them); x1 is the equation of first correction, which takes into consideration the influence of non-linear terms in first approximation; x2 is the solution of the equation of second correction, etc.
The author of paper [3] proved effectiveness of using SPM for analyzing NEC with minimal errors. According to SPM, an approximated solution of nonlinear equations, which describe processes in nonlinear circuits, is found in the form of the functional series by increasing powers of small parameter. Substituting such a series into the original equation and representing it with the system of equations by powers of small parameter, all the functions, included in the solution, are found.
Taking into consideration the above shortcomings of the described numerical methods of analysis of NEC, as well as the fact that SPM is usually implemented in the time domain, which is accompanied by complex trigonometric transformations and decreases its effectiveness, it is necessary to search for a new method or to improve the existing one. That is why in order to enhance efficiency of analytical calculations, it is proposed to implement SPM in the frequency domain using the automated method for forming orthogonal harmonic components of electrical magnitudes [20, 31] based on the algorithm of discrete convolution [32].
3. The aim and objectives of the study
The aim of the research is to develop an effective and simple in practical terms analytical method for analysis of nonlinear of electrotechnical system due to the implementation of calculations in frequency domain with a possibility of their automation, which will give the opportunity to perform a qualitative analysis of the researched system, specifically, will make it possible to assess the impact of the parameters of a circuit on the harmonic composition of current.
To accomplish the set goal, the following tasks were stated:
- to develop SPM through its implementation in frequency domain due to the use of the automated method for forming electrical magnitudes based on the algorithm of discrete convolution;
- to carry out a comparative analysis of SPM with implementation of calculations in frequency domain that is developed with the classic SPM in time domain and with the results of calculation of a mathematical model of the researched electrical circuit.
4. Analysis of calculation of nonlinear electrical circuits of electrotechinical systems with the use of SPM in the frequency domain
4. 1. Algorithm of the analytical method for NEC analysis in the frequency domain based on SPM
SPM in time domain is the prototype of SPM, implemented in frequency domain. It should be noted that the implementation in frequency or time domain does not affect accuracy. In this regard, only accuracy of approximation of VAC of nonlinear function, as well as the number of harmonics, which are analyzed are important.
But it will cause complications of analytical expressions. As a result, calculation in time domain will be time-consuming and accompanied by difficulty associated with multiple multiplication of the current and voltage signals that are assigned in the form of the Fourier series, trigonometric transformations of cosine sine functions, double angle functions, etc.
The ultimate result of trigonometric transformations should be expressions that contain functions cos(krat) and sin(krat), where k is the number of harmonics. In addition, trigonometric transformations are difficult to automate. All this can lead to obtaining an incorrect and unpredictable result. That is why it is proposed to develop the SPM through implementation in frequency domain. The main advantage of SPM in frequency domain is that trigonometric transformations, associated with multiplication of current and voltage signals, are implemented with the use of the operation of discrete convolution, which will ensure adaptation to software automation of analytical calculations. This allows us to significantly decrease time consumption and to simplify the process of analytical calculation. It also makes it possible to obtain the predicted result regardless of degree of an approximating polynomial, as well as the number of analyzed harmonics [9, 10, 13, 23].
The algorithm of analysis of NEC in an analytical form based on SPM in the frequency domain, shown in Fig. 1, includes the following stages:
1. We write down the expression that determines the form of nonlinearity of the researched electrical circuit:
^ = f(i),
(2)
where is the nonlinear parameter. According to the classic SPM, nonlinearity is represented in the form of instantaneous resistance, which nonlinearly depends on current.
2. Using the automated method of forming components of instantaneous variables of electrical circuit in frequency domain based on the operation of discrete convolution, we form the equation with the use of the second law of Kirchhoff. The expression for voltage on a nonlinear element is written down in the left part, while the differences of voltages on the power source and all consumers, except for the nonlinear one, are written down in the right part.
3. To solve the resulting nonlinear equation, current is represented in the form of a power series by nonlinear parameter
It should be noted that for a reliable analysis of the processes occurring in electrical circuits, as well as for taking into account the features, introduced by nonlinear elements, it is necessary to select a higher power of the series that describes the current than the order of the power series that describes the nonlinear element.
4. To represent (3) in the frequency domain, additional designations are introduced. According to SPM, it was accepted that orthogonal cosine aa0 and sine aj0 elements of coefficient a0, which is a generating solution and describes the work of the linear part of the circuit, is determined by the first harmonic:
(4)
where Ao, Bo are the coefficients that depend on parameters of a circuit, cosine and sine, respectively
According to SPM, which is implemented in frequency domain, harmonic composition of coefficients ao, ai, power series (2) is formed by the following principle. Coefficient ai must include harmonic composition, which corresponds to harmonic composition of coefficient ajJ+i raised to power m in frequency domain, where h is the maximum power of current in the expression that describes current on the nonlinear element.
In other words, the cosine aai and sine abi components of ai in the general case in frequency domain will be made up of
' oN ' oN
= Ao = Bo
Vo J Vo J
the sum of coefficients
I Ai„
and
IBim , where Aib
Bik are the components of coefficient ai, which depends on the parameters of a circuit; k is the number of harmonic component that forms the coefficient; K is the minimal number of the harmonic component of current, which form the coefficient, m is the number of harmonic of current.
In its turn, harmonic composition of coefficient a2 will be determined as convolution a0h ■ ai in frequency domain.
Thus, an will be formed according to the maximum degree of current in the expression that describes voltage on the nonlinear element and its harmonic structure:
( k \
I An,
= |I Bn,
(5)
5. Arrays of harmonic components of current in the general form are formed:
Iam = Ao +I^ll(An,
n=i V k=i
hm = Bo +I ^ (I (Bn,
(6)
m
m
m
m
I = a0 + aji.+ a2|i2... + ajan, (3)
where a0, ai, a2, an are the corresponding coefficients of the power series; n is the order of the power series, and of the corresponding number of the coefficient.
where N is the maximum order of the power series.
6. (6) are substituted in the equation of balance of voltages of the circuit, coefficients at equal powers of small parameter are grouped, and the system of equations, from which we determine coefficients A0, B0, Anm, Bnm, is formed.
The expression that determines the form of nonlinearity of the studied electric circuit is assigned: |x=f(i).
The equation using the second law of Kirchhoff is written down.
T
Current is represented in the form of a power series for nonlinear parameter ji:
I = + a\\i +... +an\i" ■ (a)
To represent the power series in frequency region, additional designations are introduced and the arrays of real and imaginary components of correspondent coefficient are written down:
aao = A> ; ato = -So;
«.= IX ; «HZ*.
(b)
The arrays of harmonic components of current in general form are formed:
n= 1 \k=l J,„
^,=5„+lM"fi(S,„„)l ■ (c)
»=1 \k=l J„,
Substituting (c) in the equation of balance of voltage of circuit, coefficients at the same powers of small parameter are grouped and the system of equations, from which coefficients Ao, B0, Anm, Bnm are determined, is formed.
As we know, NEC are described by differential equations, the solution of which allows analysis of the system. For the researched system, shown in Fig. 2, the differential equation of the balance of voltages with dependence of parameters of a circuit on time that describes its operation of the system:
uS (t ) = i (t )R +
d y (t )
dt
(7)
where y is the current linkage.
Nonlinear dependence of current on current linkage is approximated by the polynomial of third degree:
i (y) = ay + by3,
(8)
where a, b are the corresponding coefficients of a polynomial. During the studies that were conducted with the aim of selecting the optimal function to describe a nonlinear characteristic in the software package CurveExpert, it was discovered that function (8), that is, a polynomial of third degree in the absence of quadratic term has the highest coefficient of coincidence with the curve of dependence of current on current linkage.
As a result of numerical construction (8), the magnetization curve of nonlinear inductance was obtained (Fig. 3), in this case, numerical values of coefficients were accepted as a = 0.1, and b = 40.
Fig. 1. The algorithm implementation of SPM in the frequency domain
4. 2. Calculation of the electrical circuit with nonlinear inductance using SPM in the frequency domain
Consider the example of calculation of the electrical circuit, composed of sequentially connected nonlinear inductance and linear active resistance (Fig. 2).
Fig. 3. Magnetization curve of nonlinear inductance
i(t)
u(f)<Q
L(i)
R
To have the possibility of realization of analysis of the researched electrical circuit using SPM in the frequency domain, differential equation (7) is written down in the frequency domain:
USa = URa + ULa USb = URb + ULb,
(9)
Fig. 2. Electrical circuit with nonlinear inductance
Power voltage uS(t) is assigned by sine component of the first harmonic: uS(t) = Ubisin(rat). For analysis, current is represented in the form of the first and the third harmonic components:
i (t) = Ia1 cos (cot) + Ib1 sin (cot) +
+ Ia3 cos (3rat) + Ib3sin (3c).
For an analysis of equivalent scheme of the researched circuit (Fig. 2), the following numerical values of parameters of the circuit were accepted: Ub1 = -12.5 W, R = 50 Ohm.
where USa, URa, ULa are the cosine component of power voltage, voltage on active resistance R, voltage on nonlinear inductance, respectively; USb, URb, ULb are the sine component of power voltage, voltage on active resistance R, voltage on nonlinear inductance, respectively.
Cosine and sine components of power voltage (10), voltage on active resistance (10) in frequency domain were represented according to SPM, implemented in the frequency domain.
Usa =
' 0N ' 0 N
0 0 U = Usbl 0 . (10)
, 0, , 0 ,
Ua =
-R
R
- y ala - 3y 3alb - 3 y aly ^ - 6y aly 2a3b
-6ValV23b - 3V2lVa3b - 6ValV43Vbl4 +3V 2lV a3b 0
Va3a + V3alb - 3 ValVMb + 6V3l V23b +
+6V a 3 V hb + 3V a 3' 0
3 b + +3V23Va3b
-3bR
-Va3Val + 2ValV43V4l + Va3VM -
R
-R
Ur4 =
-3Rb
Ka3 -V43VM -V3alV43" A „, i „,3
-ValV23 - 2V«3V43V41 + VLVal . 0
-3bR [-V«1V23 + ValV23 + 2 Va3 V43 V41 ] 0
4RVa3 [V23 - 3V43 ] 0
V 4ia + 3 V 2lV 414 + 3 V h4 + 6V 4lV234 +
+6V 41V234 - 6V alV a 3 V 414 - 3 V 43 V 24 + +3V 2lV 434 0
-V43a - 3VIlV414 + 3Vh4 - 6V4lV434 -
-6V43 V«i4 - 3Va3V434 - 3VM4 ] (11) 0
-2V 4lV alV,
-2ValV43V«3 - V23V41 + V^V41 0
■¿rn [V 41V 23- 2 V 23 V 41+2 V «lV 43 V «3 ] 0
4RV43 [3V«3 - V23 ]
According to (7), voltage on nonlinear inductance is determined as dy(i)/di, that is why for further representation in frequency domain, cosine and sine components of current linkage, assigned by the first and third harmonics are written down first:
(l2)
where yal, ya3 are the cosine components of the first and third harmonics of current linkage, respectively; y^, V43 are the sine component of the first and third harmonics of current linkage, respectively.
Then, a derivative is determined from the obtained arrays.
' 0 N ' 0 N
V al V 4l
Va = 0 V a = 0
\ Va3 / ,V 43 ,
' 0 N ' 0 N
-V 4lC 0 • d V 4 = V alC 0 . (l3)
v-3V 43Cy , 3V a3™
1am
14m
V ala + 3 V 3alb + 3 V alV Mb + 6V alV 23b +
+6ValV23b + 3V2lVa3b + 6ValV43V414 - 3V2lVa3 0
V a3a + V 3alb - 3 V alV M4 + 6V 2lV a 34 +
+6V a3 V h4 + 3Va3b + 3V 23V a 3b 0
■Va3V2l + 2ValV43V4l + Va3VM -V alV 23 - 2V a3 V 43 V 4l +V LV al _ 0
-3b [-V alV a 3 + V alV 23 + 2V a3 V 43 V 4l ] 0
4Va3 [V23 - 3V23 ] 0
-3b
V 4la + 3V 2lV 4lb + 3V 4lb
"6V 4lV 23b"
(l4)
+6V 4lV 4,3b - 6V alV a 3 V 4lb -
-3V 43 V 24 + 3V 2lV 43b 0
V43a + 3ValV4lb - 3V4lb + 6V2lV43b
+6V 4 3 V 2lb + 3V 33 V 43b + 3 V 4,3b 0
3b i2V 4lV alV a 3 - V 4 3 V 2l + V alV 43 + [+2ValV43 Va3 + VmV4l - V33V4l. 0
34 [V4lVa3 - 2V23V4l + 2ValV43Va3 ] 0
4V43 [-3Va3 + V43 ]
According to the algorithm of SPM (Fig. 1) and taking into account expression (7), we will represent current linkage in the form of a polynomial of the first degree by nonlinear parameter 4:
V = a0 + alb,
(l5)
where a0, al are the corresponding coefficients of the polynomial.
According to the harmonic composition of current linkage y(t) and algorithm of implementation of the SPM in frequency domain, coefficient a0 will be determined from the following expression:
a0 = ^cos (cot ) + B0sin (cot),
(l6)
d V a =
According to (8), using an automated algorithm of formation of components of electrical magnitudes in frequency domain, the cosine and sine components of current linkage in third power are determined. Then, the arrays of cosine and sine components of current, flowing in the researched circuit, are formed:
since they are the generating solutions.
In accordance with the rule of formation of harmonic composition al for a given electrical circuit, after analyzing a03, we obtain that al is formed by third harmonic in addition to first harmonic:
al = ^llcos (cot ) + ^3cos (3c )-+5llsin (cot ) + Bl3 sin (3rat ),
(l7)
where A0, B0, All, Bll, Al3, Bl3 are the coefficients that depend on parameters of a circuit.
Then the polynomial that describes the current linkage in the researched circuit will take the form:
y = A cos (coi ) + B0 sin (cot ) + A11cos (cot ) + Aj3cos (3rai ) + +B41sm (cot ) + B13 sin (3rni )
(18)
' 0 N ' 0 N
A0+Ab B0+Bnb
V a = 0 ; v b = 0
, A13b , , B13b ,
Coefficients in (l8) are grouped by the appropriate trigonometric functions and frequencies, and the arrays of cosine and sine harmonic components of current linkage in frequency domain are written:
(19)
I. e., Val = A0+Anb, Vbl = Bci+Bllb, Va3=A^b, Vb3=Bl3b.
The above expressions are substituted in the equation of balances of voltages in frequency domain, and coefficients at the same powers of the nonlinear parameter are equaled. It should be noted that we will take into consideration only coefficients of zero and first degrees b. And then, the system of equations is obtained:
Ub1 = cA0 + aRB0 ; 0 = aRAj, - raB0 ;
0 = 3RB03 + raA11 + 3RB0 A,2 + aRB11;
0 = aRA13 + RAJ - 3raB13 - 3RA,B02 ;
0 = aRB13 + 3RB0 A|2 + 3raA13 - RB03 ;
Ub1 = A0ra + aB0R + (3RB03 + 3RB0 A + raA11 + aRB11 ) b;
0 = -B0ra + aA^R + (3RA„3 - raB11 + aRA11 + 3RA„B02 ) b;
0 = (aRA13 + RA03 - 3raB13 - 3RA0B02 ) b;
0 = (aRB13 - RB03 + 3mAl3 + 3RB0A„2 ) b,
from which coefficients A0, B0, All, Bll, Al3, Bl3 are determined:
A0 =
B=
A11 =
Ub1C
B=
A13 =
B=
ra2 + R2a2 ' Ub1Ra (Ra -ra2 ) -ra4 + R4a4 ' -6Ub13R2 a ; ra3 + rnR2a2 )2 ' 3Ub13R .
(ra2 + R 2a2 )
(3rn4 - 12R2 a 2rn2 + R 4a4 )Ub13aR2
(9rn2 + R2 a2 )(ra2 + R 2a2 )3 (3rn4 - 12R2 a 2rn2 + R 4a4 )Ub13R ( + R2 a2 )( + R2 a2 )3
The obtained coefficients have the following numerical values: A0 = -0.04, B0 = -6.33740-4, All = 9.59340-7,
b11=-3.0l2^l0-5, a13=-l.776^l0-7, b13=-a.a4a•lo-6.
Substituting the obtained analytical values of the coefficients in arrays (l8), we determine analytical dependences and numerical values of harmonic components of current linkage of the researched yal, Vbl, V«3, Vb3 in the function of coefficients A0, B0, All, Bll, Al3, Bl3. Substituting yal, Vbl, Va3, Vb3 in expressions (l4), analytical dependences and numerical values of harmonic components of current of the researched circuit are determined.
4. 3. A comparative analysis of the developed and classic SPM
The calculated values of harmonic components of current with the use of classic SPM and the proposed SPM using instantaneous resistance frequency domain were compared with the results of numerical modeling of the researched circuit in the environment of MATLAB package (Fig. 4).
Fig. 4. Curves of current of the researched circuit:
(......) is the one, obtained using the proposed SPM,
(— - ■) is the one, obtained using the classic SPM,
(-) is the one, obtained as a result of calculation
of the mathematical model
As evaluation criteria of accuracy assessment for the two compared methods, we selected the relative error by magnitude of harmonic components of the sought current and by amplitude and phase of first and third harmonics (Table 1). In addition, we assessed the degree of coincidence of calculated curves of the current signal with the curve, obtained as a result of the numerical experiment (Fig. 4), by coefficient of determination R2.
Table 1
Results of comparative analysis of the developed and classic SPM
Magnitude Relative error Developed SPM, implemented in frequency domain Classic SPM, implemented in time domain
Harmonic components of current S(/a1>, % 6.3l3 8.5
S(/b1>, % 0.39 4.7
S(/a3>, % 2.3l6 8.l
S(/b3>, % 0.46 24.2
Amplitude of first and third harmonics of current S(/1>, % l.786 3.988
S(/3>, % 2.273 7.577
Phase of first and third harmonics of current S(j/1>, % l.04M0-4 2. 4ll0-4
S(j/3>, % l.042-l0-3 4.5840-3
Coefficient of determination R2 0.993 0.962
The conducted comparative analysis of the proposed and the classic method of the small parameter with the results of numerical calculation of the mathematical model of the electrical circuit with RL-load showed that the developed method provides greater adequacy of results and a higher accuracy of calculations.
5. Discussion of results of using the improved SPM with the implementation in the frequency domain for analysis of nonlinear electrical circuits
It should be noted that calculation accuracy depends on the accuracy of approximation of the characteristic of a nonlinear element. When implementing the analysis of the researched circuit, in order to avoid the awkwardness of analytical expressions, the polynomial dependence of the third degree was selected for calculation. This choice was determined by the conducted experimental research in the mathematical package with the aim of searching for a function that at the minimal number of terms will ensure acceptable accuracy of the description of a nonlinear characteristic.
In this case, the above numerical values of relative errors for values of the harmonic components of current do not exceed 6 % even at a relatively low quality of approximation.
Thus, the advantage of the SPM in frequency domain can be the fact that it is quite simple to implement, ensures high accuracy and effectiveness of obtaining numerical values of amplitude components of current and allows conducting an analysis in the analytical form of the processes of the researched electrical circuit with a nonlinear load.
The disadvantage of the method is the need to ignore higher harmonic components when solving the formed system of equations, which leads to a decrease in accuracy. In addition, to achieve the accurate approximation dependence of a nonlinear characteristic, it is necessary to use a polynomial function of the high degree, which will make the analytical expressions cumbersome.
Presented research shows the possibility of implementation of the SPM in frequency domain with its practical
use, specifically: the presented method enables qualitative analysis of indicators of operation of nonlinear load and the current spectrum, which is a good basis for the development of the method for selection of the tools for active and passive filtering that was already partially presented in paper [21].
In subsequent studies, the obtained analytical expressions will make it possible to explore the processes of compensation of higher harmonics of current in electrical networks with nonlinear load, containing semiconductor converters for more precise selection of the tools for both passive and active filtering.
6. Conclusions
1. The method for analytical analysis using the method of the small parameter, implemented in frequency domain, was developed. It was shown that the method makes it possible to avoid complex trigonometric transformations and to simplify analytical calculations and maximally adapt them to automation through the use of the automated method of formation of components of instantaneous values in frequency domain, which will enable a decrease in time consumption.
2. The algorithm of calculation of nonlinear electrical circuits using the small parameter method in frequency domain with the help of orthogonal components of nonlinear resistance was developed, and effectiveness of its use was proved. A special feature of the presented algorithm is automation of analytical calculation in the environment of the mathematical package. Conducted comparative analysis of the proposed small parameter method with the classic small parameter method on the example of calculation of an electrical circuit with the RL loading showed that the developed method provides greater adequacy of results and higher precision calculations in comparison with the existing one. Relative error for harmonics of current does not exceed 6 %. Relative error by the amplitude of first and third harmonics of current does not exceed 2.5 % and relative error for phase does not exceed 1.04240-3 %.
References
1. Zaezdniy A. M. Bases for nonlinear and parametric radio circuits. Moscow: Svyaz, 1973. 448 p.
2. Shydlovska N. Analysis of nonlinear electric circuits by the small parameter method. Kyiv: Evroindeks, 1999. 192 p.
3. Bessonov L. Theoretical foundations of electrical engineering. Electric circuits. Moscow: Gardariki, 2001. 638 p.
4. Basics of circuit theory / Zeveke G., Ionkin P., Netushil A. et. al. Moscow: Energiya, 1975. 752 p.
5. Energy balance in the electric circuits / Tonkal V., Novoseltsev A., Denisuk S. et. al. Kyiv: Naukova dumka, 1992. 312 p.
6. Steady-state analysis of nonlinear circuits using discrete singular convolution method / Zhou X., Zhou D., Liu J., Li R., Zeng X., Chiang C. // Proceedings Design, Automation and Test in Europe Conference and Exhibition. 2004. doi: https://doi.org/10.1109/ date.2004.1269078
7. A wavelet balance approach for steady-state analysis of nonlinear circuits / Li X., Hu B., Ling X., Zeng X. // ISCAS 2001. The 2001 IEEE International Symposium on Circuits and Systems (Cat. No.01CH37196). 2001. doi: https://doi.org/10.1109/iscas.2001.921249
8. Zhu L. (Lana), Christoffersen C. E. Transient and Steady-State Analysis of Nonlinear RF and Microwave Circuits // EURASIP Journal on Wireless Communications and Networking. 2006. Vol. 2006. P. 1-11. doi: https://doi.org/10.1155/wcn/2006/32097
9. Zagirnyak M., Maliakova M., Kalinov A. Analysis of electric circuits with semiconductor converters with the use of a small parameter method in frequency domain // COMPEL - The international journal for computation and mathematics in electrical and electronic engineering. 2015. Vol. 34, Issue 3. P. 808-823. doi: https://doi.org/10.1108/compel-10-2014-0260
10. Zagirnyak M., Kalinov A., Maliakova M. Analysis of instantaneous power components of electric circuit with a semiconductor element // Archives of Electrical Engineering. 2013. Vol. 62, Issue 3. doi: https://doi.org/10.2478/aee-2013-0038
11. Czarnecki L. S. Effect of Supply Voltage Harmonics on IRP-Based Switching Compensator Control // IEEE Transactions on Power Electronics. 2009. Vol. 24, Issue 2. P. 483-488. doi: https://doi.org/10.1109/tpel.2008.2009175
12. Czarnecki L. S. On some misinterpretations of the instantaneous reactive power p-q theory // IEEE Transactions on Power Electronics. 2004. Vol. 19, Issue 3. P. 828-836. doi: https://doi.org/10.1109/tpel.2004.826500
13. Zagirnyak M., Maliakova M., Kalinov A. Analysis of operation of power components compensation systems at harmonic distortions of mains supply voltage // 2015 Intl Aegean Conference on Electrical Machines & Power Electronics (ACEMP), 2015 Intl Conference on Optimization of Electrical & Electronic Equipment (OPTIM) & 2015 Intl Symposium on Advanced Electromechanical Motion Systems (ELECTROMOTION). 2015. doi: https://doi.org/10.1109/optim.2015.7426958
14. Zagirnyak M., Maliakova M., Kalinov A. Compensation of higher current harmonics at harmonic distortions of mains supply voltage // 2015 16th International Conference on Computational Problems of Electrical Engineering (CPEE). 2015. doi: https:// doi.org/10.1109/cpee.2015.7333388
15. Improvement of compensation method for non-active current components at mains supply voltage unbalance / Al-Mashakbeh A. S., Zagirnyak M., Maliakova M., Kalinov A. // Eastern-European Journal of Enterprise Technologies. 2017. Vol. 1, Issue 8 (85). P. 41-49. doi: https://doi.org/10.15587/1729-4061.2017.87316
16. Zagirnyak M., Kalinov A., Chumachova A. Correction of operating condition of a variable-frequency electric drive with a non-linear and asymmetric induction motor // Eurocon 2013. 2013. doi: https://doi.org/10.1109/eurocon.2013.6625108
17. Research of rnergy processes in circuits containing iron in saturation condition / Prus V., Nikitina A., Zagirnyak M., Miljavec D. // Przegl^d Elektrotechniczny (Electrical Review). 2011. Vol. 87, Issue 3. P. 149-152. URL: http://pe.org.pl/articles/201V3/39.pdf
18. Zagirnyak M. V., Rodkin D. I., Korenkova T. V. Estimation of energy conversion processes in an electromechanical complex with the use of instantaneous power method // 2014 16th International Power Electronics and Motion Control Conference and Exposition. 2014. doi: https://doi.org/10.1109/epepemc.2014.6980719
19. Al-Mashakbeh A. A diagnostic of induction motors supplied using frequency converter basing on current and power signal analysis // Przegl^d Elektrotechniczny. 2016. Vol. 1, Issue 12. P. 7-10. doi: https://doi.org/10.15199/48.2016.12.02
20. Zagirnyak M., Mamchur D., Kalinov A. A comparison of informative value of motor current and power spectra for the tasks of induction motor diagnostics // 2014 16th International Power Electronics and Motion Control Conference and Exposition. 2014. doi: https://doi.org/10.1109/epepemc.2014.6980549
21. Zagirnyak M. An analytical method for calculation of passive filter parameters with the assuring of the set factor of the voltage supply total harmonic distortion // Przegl^d Elektrotechniczny. 2017. Vol. 1, Issue 12. P. 197-200. doi: https:// doi.org/10.15199/48.2017.12.49
22. Salmeron P., Litran S. P. Improvement of the Electric Power Quality Using Series Active and Shunt Passive Filters // IEEE Transactions on Power Delivery. 2010. Vol. 25, Issue 2. P. 1058-1067. doi: https://doi.org/10.1109/tpwrd.2009.2034902
23. Ginn H. L., Czarnecki L. S. An Optimization Based Method for Selection of Resonant Harmonic Filter Branch Parameters // IEEE Transactions on Power Delivery. 2006. Vol. 21, Issue 3. P. 1445-1451. doi: https://doi.org/10.1109/tpwrd.2008.2009899
24. Luchetta A., Manetti S., Reatti A. SAPWIN-a symbolic simulator as a support in electrical engineering education // IEEE Transactions on Education. 2001. Vol. 44, Issue 2. doi: https://doi.org/10.1109/13.925868
25. Huelsman L. P. SAPWIN, Symbolic analysis program for Windows - PC programs for engineers // IEEE Circuits and Devices Magazine. 1996. Vol. 6.
26. Moura L., Darwazeh I. Introduction to linear circuit analysis and modelling from DC to RF, MatLAB and SPICE. Burlington, 2005. 405 p.
27. Raju A. B., Karnik S. R. SEQUEL: A Free Circuit Simulation Software as an Aid in Teaching the Principles of Power Electronics to Undergraduate Students // 2009 Second International Conference on Emerging Trends in Engineering & Technology. 2009. doi: https://doi.org/10.1109/icetet.2009.200
28. Pires V. F., Silva J. F. A. Teaching nonlinear modeling, simulation, and control of electronic power converters using MATLAB/ SIMULINK // IEEE Transactions on Education. 2002. Vol. 45, Issue 3. P. 253-261. doi: https://doi.org/10.1109/te.2002.1024618
29. Chen W.-K. Feedback, Nonlinear, and Distributed Circuits. 3rd ed. CRC Press, New York, NY, 2009. 466 p.
30. Zagirnyak M., Kalinov A., Maliakova M. An algorithm for electric circuits calculation based on instantaneous power component balance // Przegl^d Elektrotechniczny (Electrical Review). 2011. Vol. 87, Issue 12. P. 212-215. URL: http://pe.org.pl/articles/ 2011/12b/59.pdf
31. Nayfeh A. H., Chen C.-Y. Perturbation methods with mathematica // Nonlinear dynamic. Wiley, 1999. 437 p.
