CC BY
26
ISSN 2521-6244 (Online) Po3gi. «E.eKTpoTexmKa»
UDC 621.314
MATRIX-TOPOLOGICAL MODEL OF ELECTROMAGNETIC CIRCUITS
SIDDELEV N.I. Ph.D, Docent, Associate Professor of the faculty of automatics and the computer-integrated technologies of the Black Sea National University of a name of Peter Mogila, Nikolaev, Ukraine, e-mail: [email protected].
Purpose. To develop a digital model of electromagnetic devices for research and optimization of powerful secondary electric power sources and electromagnetic converters.
Methodology. Nodal potential method, Contour current method, Topologically isomorphic transformations.
Findings. The purpose of this work is to create a mathematical apparatus that allows solving problems of modeling and researching electromagnetic devices in parts (by types of accumulated energy). This will simplify the research and optimization of technical characteristics such as efficiency, weight and size indicators, etc. The proposed mathematical model of electromagnetic circuits has the greatest degree of detail of the electric and magnetic circuit. The magnetic circuit is represented in the same detail as the electric circuit, and is described by a contour matrix. A mathematical description of electromagnetic devices is obtained in which inductive parameters are determined by the geometric dimensions and characteristics of magnetic circuits. The topology of the electrical circuit is represented by matrix blocks, which allowed obtaining a mathematical description, which simultaneously takes into account the distribution of currents and charges in the elements of the circuit. The system of equations reduces to the Cauchy form and is composed with respect to increments of magnetic fluxes and potentials on capacitors, which simplifies its solution by numerical methods on a computer. Thus, it is convenient to monitor the energy processes in the reactive power-consuming elements of the circuit. A stable and adaptive digital model of electromagnetic circuits has been developed that makes it possible to combine several methods of integrating a system of differential equations. Feedback is provided through a special parameter. This makes it possible to maximize the correctness of the computations for the energy components in the simulation of the electromagnetic circuit. The originality of the mathematical description lies in the fact that the topology of the electromagnetic circuit is represented in the form of separate matrices that are connected by a matrix of coil connections. The practical value of the digital model of the electromagnetic circuit is that the parameters of the magnetic circuits are introduced in the form of geometric dimensions of the magnetic circuits. This eliminates the need for equivalent transformations to produce data for a specific model. This simplifies the study of secondary power supplies and other powerful electric power consumers by efficiency criteria, weight and size parameters.
Originality. The topology of the electromagnetic circuit is represented in the form of separate matrices that are connected by a matrix of coil connections.
Practical value. The parameters of the magnetic circuits are introduced in the form of geometric dimensions of the magnetic circuits.
Keywords: electric circuit; magnetic circuit; static electromagnetic devices; topological isomorphic modeling; matrix topological description; topological matrices; incidence matrices; block structure of the topological matrix; matrix of helical links; secondary power supplies; pulse current generators; topological isomorphic model.
I. INTRODUTION
Modern computer-aided design systems require the development of special mathematical software, which will provide the most simulation for the devices being developed. The main requirements for the development of models can be the greatest degree of detail, the permissible quality of modeling, the simplicity of obtaining model parameters.
A special role is played by design automation in the electric power industry [6], [11]. In autonomous electric power systems, all devices can be divided into three groups: power supplies, converters (or secondary power supplies), and electricity consumers.
Among them, devices that belong to the second group, by mass and size, are sometimes commensurate with power supplies and often exceed by these parameters consumers of electricity. In addition, electricity converters
© Siddelev N.I., 2018
DOI 10.15588/1607-6761-2018-1-1
are also a kind of energy consumers, which is used to control switching elements and is released as heat. The amount of energy consumed by the converter affects the economics of the autonomous system for the worse, and is therefore often the criterion for research.
Today, computer modeling of electrical and electronic devices continues to play a special role in developing economically and technically efficient devices for various purposes. This allows not only to speed up the calculation work and improve the design quality, but also to create favorable conditions for the development of the theory of modeling energy-intensive parts of electric power converters.
The complexity of studying nonlinear electric and magnetic circuits is that there are no general analytical solutions of systems of nonlinear algebraic and differential equations. It is generally believed that the analysis of essentially non-linear chains, with the exception of a
ISSN 2521-6244 (Online) Po3gi^ «E.roKTpoTexmKa»
small number of simplest cases, can be performed only numerically.
This material was compiled from the results of [12], [20] - [21].
II. ANALYSIS OF LAST RESEARCHES
When modeling devices, it is very important to choose a method that allows you to accurately estimate the accumulation and loss of energy in the elements of the circuit. This may affect the final result of the study, for example, the efficiency, the power factor, and the like.
In the study of devices using computers, the principles of structure representation and methods for describing electromagnetic circuits (EMC) play an important role. The principles of topological-isomorphic modeling (TIM) are well studied and successfully used in the automation of EMC design [12]. However, modern electronic devices with electromagnetic components in their composition usually have a large number of capacitive elements, which requires the further development of TIM methods and methods of EMC description. One way to solve this problem is to divide the EMC into subcircuits.
Choosing the appropriate description for the study of physical processes, you should take into account the purpose, direction of development, the degree of detail of the device nodes, are modeled. One of the most important problems of automation of scientific research of electromagnetic devices (EMD) is the study of energy processes and characteristics (efficiency, power factor, etc.). In this case, it is necessary to take into account the energy losses at each element of the circuit. A number of devices (electromagnetic converters, electric and magnetic energy storage devices, power supplies of electronic equipment, etc.) include powerful magnetic systems where there are losses to magnetic hysteresis and eddy currents. To take into account the losses in steel, it is necessary to represent the magnetic circuit in the same detail as the electric circuit.
At present, methods for representing EMC have been developed, where the separation of the circuit into electric and magnetic circuits is used. To represent the graph of a magnetic circuit, use an incidence matrix or, more often, a contour matrix of magnetic bonds. The connection between the magnetic and electric graphs is provided by the matrix W , the dimension of which depends on the number of turns located on the magnetic rods and the method of connections [12]. There are two ways of forming W . In one case, the connection of the electric branch to the magnetic flux circuit is realized, in the other case the connection of the contour stream to the electric branch is realized. The choice of the matrix is determined by the concrete implementation of the mathematical models of the EMC.
The structure of electrical circuits is described by matrices of incidents, contour currents, cross sections, etc., makes it possible to widely use the theory of matrices in the compilation of EMC descriptions. In turn, the the-
ory of matrices allows solving the problems of modeling electric and magnetic circuits by parts [5], [7], [13], [16] -[17], [23]. If all known types of elements (resistors, capacitors, inductors, etc.) are present in the circuits, then the matrices that retain the electrical circuit structure do not agree well with the mathematical description. In such cases, the electrical circuit is divided into subcircuits in which elements of the same type are localized. For example, it is possible to isolate the linear resistive part of the circuit, and the nonlinear elements and components accumulating energy are replaced by current sources and electromotive forces. Thus, the designer of EMD must deal with the transformation of the original and the construction of a new (equivalent) scheme.
III. FORMULATION OF THE WORK PURPOSE
In view of the above, it is possible to formulate the goal of the work as follows: to develop a digital model of electromagnetic circuits that is convenient for researching and optimizing powerful secondary power supplies and electromagnetic converters. In this model there should be the following:
- separation of the elements of the scheme according to the types of accumulated energy of the magnetic and electric fields and the consumed energy radiated in the form of heat;
- representation of the magnetic components of the circuit in the form of their design parameters and the magnetic properties of the material used;
- eliminating the need for an equivalent transformation scheme;
- adaptability to the use of application packages, which by their organization use a matrix representation of data (for example, MATLAB).
IV. DEVELOPMENT OF THE MATRIX-TOPOLOGICAL MODEL OF ELECTROMAGNETIC CIRCUITS
The review of the diagrams of the power part of the EMU proposed in the theoretical and reference literature [1] - [2], [8], [10], [19], [24], etc. makes it possible to use the following fragment of an electrical circuit, on the basis of which a mathematical model will be developed (see Fig1).
Limitations on the number of branches and the complexity of magnetic structures are not superimposed. Exceptions are schemes with multiple multiplication of voltage (for example, the Latour scheme [12]), unified circuits for high-voltage transformer-rectifier modules [10] and some others, the complexity of which lies in the fact that in them capacitive couplings can not be represented as a "star" (see Fig.2).
In these cases, it is necessary to solve the problem of charge distribution between capacitors or to transform the circuit by introducing an active resistance into some capacitive branches.
ISSN 2521-6244 (Online) Po3gi. «E.eKTpoTexmKa»
Figure 1. Fragment of the electrical circuit of a static electromagnetic device
«Pel 9e3
<Ps2 <Ps4
a
9e 1 •—
"el
b
<Pe4
<Pe5
Uf;5
Figure 2. Capacitive subcircuits
ISSN 2521-6244 (Online) Po3gi. «E.eKTpoTexmKa»
Simulation of electrical circuits is greatly simplified, if it can be solved with respect to electrical voltages on capacitors. This follows from the fact that the number of capacitive elements, as a rule, does not exceed 40-50% of the number of components in the circuit. Therefore, the problem of numerical modeling is simplified in the sense of reducing the number of unknown variables.
Let the electrical circuits for which a mathematical model is develop are represented with two-pole components.
Then, depending on the nature of the components, the branches of the graph of electrical circuits are divided into the following types:
1) branches of passive two-ports (capacitances);
2) branches of dependent and independent voltage sources with linear and nonlinear active resistances.
Considering branches by types, the graph of the electrical circuit can be conditionally divided into subgraphs: "capacitive" and "resistive". In the future, the concept of "subgraph" will be replaced by a "graph" when we consider them as separate structures. If you first number the branches and nodes of the arcs of the capacitive graph, and then the branches and nodes of the resistive graph, then the incidence matrix, which reflects the structure of the electrical circuit, is automatically divided into blocks:
The parameters of the magnetic branch can be calculated as follows (see Fig. 3).
A =
As Aer Ao Ar
where matrix blocks are:
As - incidence matrix of incoming branches of the graph of the capacitive part of the electronic circuit;
Ao - incidence matrix of outgoing branches of the graph of the capacitive part of the electronic circuit;
Acr - incident matrix, where the resistive branches are incident with the capacitive;
Ar - incident matrix, where the resistive branches are not incident with capacitive branches.
Thus, the "capacitive" graph is described by a topological submatrix Ac. The basic nodes of it are placed in a separate block Ao . This block provides a connection between the base nodes of the "capacitive" graph and the "resistive" graph, which is represented by the submatrix Ar . Other nodes of the capacitive graph can be and are most often associated with a "resistive" graph, and this relationship is described by a submatrix Aer.
Description of the structure of the magnetic circuit and its connection with the electrical circuit
The magnetic circuit is represented in as much detail as the electric circuit and is described by a contour matrix rm . In order to describe scattering fluxes, when forming a graph, each magnetic connected chain can be supplemented by an "air" branch [12].
Figure 3. Constructive parameters of the magnetic branch
The rod made of magnetic material (length - I, cross-section - 5, gap - 3) is wound with one or several windings with a number of turns w. The magnetic resistance of the branch is calculated by the formula
l-5 5
■ + -
where ^ - magnetic permeability of a rod, - magnetic permeability of a gap.
The cross-section of the gap and the rod is assumed to be the same under the condition that 5<<l. The inductance of the branch can be determined by the expression:
L = W..
The connection between the magnetic and electric graphs is provided by the matrix W. Its dimension in rows corresponds to the contour matrix rm, and by the col-umns-the conductivity matrix Y. Elements of the matrix are filled in the following way. If the beginning and end of the electrical winding coincide with the chosen direction of the magnetic flux, then the number of turns with a plus sign is inscribed in the corresponding cell, if the directions do not coincide, then the number of turns with a minus sign.
Topological-isomorphic description of the electro-magnetic circuit
Consider the electrical circuit of the static electromagnetic device in Fig.1. The numbering of branches and knots is performed according to the above rules for forming a matrix A. From the general scheme, a capacitive subscheme can be distinguished (see Fig. 2, a).
In what follows we confine ourselves to the consideration of capacitive subgraphs of the "star" type. For graphs such as "triangle" you need a special algorithm for the redistribution of charges. After dividing the circuit into resistive and capacitor nodes, we turn to the problem
r
m
ISSN 2521-6244 (Online) Po3gi. «E.eKTpoTexmKa»
of the distribution of currents between the branches of both subgraphs.
The instantaneous stresses of the branches of the electric graph are determined by the potentials of the nodes and the voltages on the capacitors. Let the node 9m+2 (see Fig. 1) is not an incident capacitor and it is necessary to find the voltage drop on the branch vk+5. The sought voltage is determined
rk+5 - e, + w
k + 5
by the potential difference 9m+2 and 9m+3. The algorithmic complexity consists in determining the node's potential, if the capacitor is incident with it. Define it as:
9m + 2 = 9m + 1 + u
cl '
where uc1 = -tyc1 - voltage on the capacitor (see Fig. 2, a).
Then
9m + 2 = 9m + 1 - 9c1 ■
Using the block structure of the matrix A, determine the distribution of currents in the circuit as follows:
I = Y[aX + Asr (vs- AsA0Vr)+ E - W'O'j, (1)
where I, E, O - respectively, currents, electromotive
forces and magnetic fluxes, and the component W'O' takes into account the voltage drop across the magnetic coils.
Expression (1) implements the method of nodal potentials using separate matrix blocks A: As, Ar> Ao,
Ar.
The relationship between the magnetic flux and the current of the electrical circuit is determined as follows:
WI = r R F®,
m m m '
where Rm - diagonal matrix of magnetic resistances.
(2)
The derivative of the voltage on the capacitor is determined by its charge current:
AECAEVe = ASrI.
(3)
We draw a section covering both nodes to which the capacitor is connected. Then the sum of the currents flowing in the cut branches according to the first Kirgoth law will be zero. Using the topological matrix A, this law is realized as follows:
ArI A0AsAsrI = 0 •
(4)
The above relations relate all the main variables that determine the state of the electric and magnetic circuit. Thus, the matrix equations (1) - (4) form a complete system of differential equations of the EMC.
The matrix form of the representation of the electromagnetic system is ideally suited to the realization of a
model on electronic computers using application software packages, for example, MATLAB, in which all objects of computation are represented by matrices [19].
The mathematical description (1) - (4) provides for the use of resistors, capacitances, inductive elements and EMF sources. However, in the converter devices there are elements with nonlinear characteristics - controlled and uncontrolled semiconductor diodes, transistors, cores of magnetic material and the like. To model the converter devices with these elements, without going beyond the scope of the description (1) - (4), requires the use of special mathematical models of the above elements, which must be added to this model in digital implementation on a computer. In this connection, only numerical methods can be used to solve system (1) - (4). The abundance of state variables in it: I, O, O', VE, VE', Vr, - complicates the solution, in addition, requires the use of special iterative methods.
Iterative methods for solving systems of differential equations are well studied in [20], [22] and are successfully used in numerical study of nonlinear systems [20]. Attempts have been made to develop original simple fast-acting methods [22], but the problems of the dimension of problems, algorithmic difficulties in the iteration step and their number remain.
In our case, the use of iterative methods of investigation would not be effective, since the task posed - energy research, in particular taking into account dispersion and accumulation of electrical energy on each element of the scheme - will require additional computational costs of the computer time [20].
In addition, modeling of nonlinear magnetic characteristics of materials requires the development of special models that should be included in the general algorithm for the study of chains.
On the basis of the foregoing, machine time must be saved, and this is achieved by using explicit methods of numerical integration, the main advantages of which are described in [9].
Reduction of the topologically isomorphic description of the EMC to a form convenient for solving by numerical methods
When developing a mathematical model on a computer using explicit methods of numerical integration, it is necessary to solve the system (1) - (4) with respect to the derivatives of the magnetic flux vectors O' and the capacitor voltages VE'.
To do this, we perform matrix transformations as follows. We re-group the terms in (1) as follows:
I = Y[(Ar - AgrAs A0 X + 4rVs + E - W'O'j (5)
We introduce a new topological matrix:
AY
= A1 - A1 A A1
- Ar sr sA0-
then
ISSN 2521-6244 (Online) Po3gi. «E.eKTpoTexmKa»
AY - Ar A0 Acr ■
If, we substitute expression (5) in (4), then we get: A/A + AJY(AtrVs + E - W*O')- 0 ■
From here you can get Vr , to exclude it from equation (5):
F =
AyYAJ
-1
AyY{W'O' - A'rFs- E
Eliminate the vector Vr from equation (5) and ob-
tain:
I = YAY AyA
-1
AjY\W O' - AerVe - Ej-
- Y(('O' - A'srVs- E)
(6)
We introduce the matrix of active parameters:
Yr = Y
Ay AiYAi AyY
where e - unit matrix.
Then equation (6) takes the form:
I - Yr AV + E - W*O
(7)
Thus, we obtained an equation in the form of nodal potentials, where all the potentials of the electric circuit graph are determined through potentials on capacitances in a capacitive subgraph.
Eliminate the vector I in the description (1) - (4) by substituting (7) into (2) and (3), which gives:
WYr
U'F + E - W'O')= r R rf O
From this we determine the derivative of the magnetic flux:
O' =
WYW-
I"1 [W (
ASrFs+ E - F)], (8)
where
F - r R r O ■
m m m
Vector F can be defined and thus: F-rmLH,
where L - matrix of lengths of magnetic cores; H - vector-column of magnetic tension in rods. We substitute equation (7) in (3):
AeCA<Ve' - -AsrVr (( + E - W*O') ■
Let us single out the derivative of the electric potential vector on capacitors:
VE' - [[CA]~JASrK (W* O'-A^rVs- E )■ (9)
Thus, the system (1) - (4) is solved with respect to
r
the derived vectors O' and Vs . System (8) - (9) is a description of the electromagnetic circuit, resolved with respect to magnetic fluxes and electrical potentials on the capacitors. A feature of this description is that it allows you to model ideal electromagnetic circuits, since it does not have redundancy. Other special cases of calculating such chains are known, for example, in [14]. In them, ideal transformers are considered as additional structural equations. However, this leads to a complication of the mathematical model. In our case, the dimension of the problem is significantly reduced. The system (8) - (9) is written in the Cauchy form and is suitable for investigation by explicit methods of numerical integration.
The problem of rigidity
As the simulation results show, the mathematical model of the EMC (8) - (9) behaved well when calculating the transformer circuit [20]. However, in the general case there can be problems associated with the rigidity of the equations and, in connection with this, with the control of the computational process.
In the general case, the components of the parameters Y, C, Rm of electrotechnical devices can differ greatly from each other. This leads to the rigidity of systems of differential equations. The concept of rigid equations was first introduced in 1952 by K. Curtis and D. Hirschfelder [21]. A feature of such systems is the need to reduce the interation step when using the classical methods of Euler, Runge-Kutta, Adams, and others:
h
(10)
where h - integration step; rbj < b - a - boundary layer [15]; [a,b] - a segment of observation (may correspond to the period of the fundamental harmonic of the process under study).
The fulfillment of condition (10) allows to display high-frequency components, which are a consequence of physical processes in charge circuits with small time constants and in high-frequency oscillatory circuits. However, the
N =
T
where T - the period of the fundamental harmonic component of the process under investigation may turn out to be such that the use of standard methods will be unsuitable because of the long calculation time.
In this case, inverse methods of numerical integration are useful: the inverse Euler method [21], the trapezium method, the differentiation formulas back [21], etc. These methods are good in that, depending on the step of numerical integration, we get the corresponding "detail" of the physical process, while the stability of the computational process is not violated.
It is known that the errors of the explicit and implicit
e
X
n.c
ISSN 2521-6244 (Online) Po3gi^ «E.roKTpoTexmKa»
methods are close in modulus and are proportional h , but are different in sign [21]. Therefore, combining these methods, you can get a good approximation. This is achieved by calculating the derivatives of state vectors with respect to different points within the integration interval [21].
If the linear approximation is performed at the integration step, then increments ao and AVS can be determined by solving a system of algebraic equations. We will show this with an example. Let
X =
o
Vc
then
AX' = BX + CE, (11)
where A, B, C - square matrices; X, E - vector columns. We write the system (11) in the following form:
X ' = A'1 (BX + CE )■
(12)
The increment at the integration step with respect to different points within the interval
aX = h
A^B\X +AX j + A~1CE
(13)
We solve system (6) with respect to AX :
aX =
1 -
hA_1B ^
z
hA'1 (BX + CE). (14)
We represent the derivatives (8) - (9) inside the interval:
O- = Fofo+Oi+^.E+A | (15)
z z
ao
av„
WYpAv + E)- F (O)
- A£PYp(A^pV£ + E)
(17)
change the methods of numerical integration, for example:
a) at z = w - explicit Eulerian method;
b) at z = 1 - inverse Euler method;
c) at z = 2 - trapezium method and others.
Parameter z can also take any non-integer value.
The square matrix in system (17) can be poorly conditioned in one case, when there are capacitors with "hanging" nodes in the circuit. This may occur in the study of chains in parts. In this case, it is recommended to shunt the "hanging" capacitors with resistors in one of the above ways (see Fig. 4). The variant (Fig. 4, b) is preferable, since it requires fewer additional elements.
Thus, a controlled digital EMC model (17) has been developed that satisfies the requirements of stability under stiffness conditions for systems of differential equations. Obviously, it will be necessary to define a criterion for optimizing a numerical solution that allows us to control the computational process by means of z .
v£ = Fv\o',Vs+A^,E+AEI (16)
If you perform operations similar to (12) - (14) over system (8) - (9), you can get:
h ~1WYpWt + z-rmRmrm - h~1ASpYpWt
- z ~1WYpAtSp
h ASCAS +z ASpYpAp
Relatively ao and aVS system (17) can be solved by any numerical method. Varying the value z, we can
Figure 4. Methods for improving the conditionality of the matrix (17) by shunting with resistors of "hanging" capacitances
IV. CONCLUSION
1. The topology of the electrical circuit is represented by matrix blocks As, A,r, Ao , Ar, which made it possible to obtain a mathematical description that simultaneously takes into account the distribution of currents and charges in the elements of the circuit.
2. A mathematical description of electromagnetic devices (1) - (4) is obtained, in which inductive parameters are determined by the geometric dimensions and characteristics of the magnetic cores.
3. The matrix form of the representation of the electromagnetic system ideally approaches the implementation of the model on electronic computers using application software packages, for example, MATLAB, in which all the calculation objects are represented by matrices.
z
X
ISSN 2521-6244 (Online) Po3gi^ «E.roKTpoTexmKa»
4. Proceeding from the given structure of the mathematical description of the EMC, it is expedient to use the proposed model for the study of devices that have magnetic cores in their composition that are powerful enough to account for them. This model was used to investigate secondary power supplies and powerful current generators.
5. In this paper, special attention is paid to the simplicity of presenting the initial information, excluding equivalent transformations of the electrical circuit.
6. A stable mathematical model of electromagnetic circuits in the matrix form (17) is developed, which is convenient for realization on digital computers. The model is composed relative to the increments of magnetic fluxes and potentials on capacitors.
7. Computational methods are controlled by the parameter z.
The developed software is well-algorithmized on digital computers and can be used to create specialized software systems.
REFERENCES
[1] Anisimov, Y. F. (1979). Marine power semiconductor technics [Sudovaya silovaya poluprovodnik-ovaya tehnika], Shipping industry, Leningrad, 192.
[2] Boyko, V. I., Zory A.A., Kornev V.D. (2009). Theoretic bases for the analysis of electronic circuits [Os-novi analIzu elektronnih shem], DonNTU, Donetsk, 306.
[3] Collatz, L. (1966). Functional analysis and computational mathematics, Academic Press, 494.
[4] Curtiss, C. F., Hirschfelder, J. O. (1952). Integration of stiff equations, Proc. Nat. Acad. Sci. USA, 38.
[5] Deskur Jan (1999). Models of magnetic circuits and their equivalent electrical diagrams, COMPEL: The International Journal for Computation and Mathematics in Electrical and Electronic Engineering, 18 4, 600 - 610.
[6] Kabyshev, A. V., Obuhov S.G. (2005). Calculation and designing of power supply systems: Reference materials for electrical equipment [Raschet i proekti-rovanie sistem elektrosnabzheniya: Spravochnyie materialyi po elektrooborudovaniyu], TPU, Tomsk, 168.
[7] Khepp, Kh. (1974). Diakoptika and electrical circuits" [Diakoptika i elektricheskie tsepi],World, Moskow, 343.
[8] Knysh, V. A. (1981). Semiconductor converters in systems capacitor charging [Poluprovodnikovyie preobrazovateli v sistemah zaryada nakopitelnyih kondensatorov], Energoizdat, Leningrad, 160.
[9] Kogan, V.L. (1982). The algorithm for calculating schemes in the basis of nodal potentials using the explicit methods of integrating [Algoritm rascheta shem v bazise uzlovyih potentsialov s primeneniem yav-nyih metodov integrirovaniya], News of high schools: Electronics, 6, 9-12.
[10]Kostikov, V. G., Nikitin I. Y. (1986 ). Power supply
high voltage REA [Istochniki elektropitaniya vy-isokogo napryazheniya REA], Radio and communication, Moskow, 200.
[11]Krasnov, V. V. (1972). Methods of mathematical modeling of ship electric power systems: Study guide" [Metodyi matematicheskogo modelirovaniya sudovyih elektroenergeticheskih sistem], NSI, Ni-kolaev, 90.
[12]Krasnov, V. V. and Siddelev, N. I. (2013). Matrix-topological description of electromagnetic circuits [Matrychno-topologichnyj opys elektromagnitnyh kil], Electrical and Computer Systems, 87, 66-73.
[13]Kron, G. (1972). The study of complex systems in parts - diakoptika [Issledovanie slozhnyih sistem po chastyam - diakoptika],Science, Moskow, 544.
[14]Mustafa G.M., Fedotov Y.B. (1983). The program for the calculation of circuits with ideal transformers [Programma dlya rascheta tsepey s idealnyimi transformatorami], Electrical industry, Conversion equipment, 5, 3-5.
[15]Ortega J., Reynolds V. (2000). Iterative solution of nonlinear equations in several variables, Society for Industrial and Applied Mathematics Philadelphia, 560.
[16]Ortner M. G., Christian Magele, and Klaus Krischan (2010). Solver for a Magnetic Equivalent Circuit and Modeling the Inrush Current of a 3-Phase Transformer, World Academy of Science, Engineering and Technology, 40.
[17]Parmantier J. P. (2007). EM topology: From theory to application", Ultra-wideband, Short-pulse Electromagnetics, 7, 3-12.
[18]Rakitskij, J. V., Ustinov, S. M., Chernoruckij, I. G. (1979). Numerical methods for solving rigid systems" [Chislennye metody reshenija zhestkih sistem], Science, USSR, Moskow, 208.
[19]Riabenkey, V. M., Dragan S. V., Solobuto L. V. (2008). Fundamentals of modeling of systems and processes in electrical engineering (Use the package of applied programs of MATLAB/Simulink) [Osnovi modelyuvannya sistem I protsesIv v elektrotehnItsI (Vikoristannya paketa prikladnih program MAT-LAB/Simulink)], New World, Lvov, 385.
[20]Siddelev, N. I. (2015). Matrix-topological description of electromagnetic circuits in the form Cauchy [Ma-trichno-topologicheskoe opisanie jelektromagnitnyh cepej v forme Koshi], Electrical and Computer Systems, 96, 63-73.
[21] Siddelev, N. I. (2017). Manage digital model based on matrix-topological description electromagnetic circuits [Upravlyaemaya tsifrovaya model na osnove matrichno-topologicheskogo opisaniya elektromag-nitnyih tsepey], Electrical and Computer Systems, 26 102, 32-40.
[22] Siddelev, N.I. (1982). On the control of iterative calculations for modeling physical processes [Ob uprav-lenii iteratsionnyimi vyichisleniyami pri modeliro-vanii fizicheskih protsessov], NSI, Nikolaev, 190, 3639.
[23]Tesche, F. M., M. V. Ianoz, and T. Karlsson. (1997).
ISSN 2521-6244 (Online) Роздш «Електротехшка»
EMC Analysis Methods and Computational Models, John Wiley & Sons Inc, New York, 656 [24] Volkov, K. V., Gubarevich V.N., Isakov V.N., Kaban V.P. (1981). Principles of construction and optimization of schemes for inductive-capacitive transducers [Printsipyi postroeniya i optimizatsii shem induk-
tivno-emkostnyih preobrazovateley], Scientific Thought , 176.
The article was received 10.05.2018
МАТРИЧНО-ТОПОЛОГ1ЧНА МОДЕЛЬ ЕЛЕКТРОМАГН1ТНИХ К1Л
С1ДСЛСВ M.I. канд. техн. наук, доцент, доцент кафедри автоматики та комп'ютерно-штегрованих
технологш Чорноморського нацiонального унiверситету iменi Петра Могили, Микола!в, Украша, e-mail: [email protected];
Мета роботи. Розробити цифрову модель електромагнтних пристрогв для дослгдження та оптимгзаци потужних вторинних джерел живлення та електромагнтних перетворювачгв.
Методи дослгдження. Метод вузлових потенцгалгв, метод контурних струмгв, топологгчночзоморфнг пе-ретворення.
Отримаш результати. Сучаснг системи автоматизованого проектування вимагають розробки спецгаль-ного математичного забезпечення. Основними вимогами до розробки моделей можуть бути найбшьша сту-пгнь деталгзацИ, допустима яюсть моделювання, простота отримання параметргв модел1. В автономних еле-ктроенергетичних системах ва пристрог можна подшити на три групи: джерела живлення, перетворювачi i споживачi електроенергИ. Серед них пристрог, що вiдносяться до другог групи, за масою i габаритами iнодi сумiрнi з джерелами живлення i часто перевищують за цими параметрами споживачi електроенергИ. Шль-юсть спожитог енергИ перетворювачем впливае в гiршу сторону на економiчнiсть автономног системи i тому часто е критерiем до^дження. Метою даноi роботи е створення математичного апарату, що дозволяе виршувати задачi моделювання та до^дження електромагнтних пристрогв по частинах (за видами накопи-чуваног енергп). Це дозволить спростити до^дження i оптимiзацiю таких технiчних характеристик, як ко-ефщент корисног дп, масогабаритш показники тощо. Запропонована математична модель електромагнтних юл мае найбшьший ступiнь деталгзаци електричного i магнтного кола. Магнтне коло представлено так само докладно, що i електричне коло, i описуеться контурног матрицею. Отримано математичний опис електромагнтних пристрогв, в которому тдуктивш параметри визначаються геометричними розмiрами i характеристиками магттного кола. Топологiя електричного кола представлена матричними блоками, що дозволило отримати математичний опис, котрий одночасно враховуе розподш струмiв i зарядiв в елементах схеми. Система рiвнянь зводиться до форми Кошi i складена вiдносно приросту магттних потоюв i потенцiалiв на конденсаторах, що спрощуе и ршення чисельними методами на комп'ютерi. Таким чином, зручно стежити за енергетичними процесами в реактивних енергоемних елементах схеми. Розроблена сттка i адаптивна цифро-ва модель електромагнтних юл, що дозволяе об'еднати юлька методiв ттегрування системи диференщальних рiвнянь. Звороттй зв'язок надаеться через спецiальний параметр. Це дозволяе домогтися максимальног коре-ктностi обчислень для енергетичних компонентiв при моделюваннi електромагнтного кола. Оригiнальнiсть математичного опису полягае в тому, що топологiя електромагнiтного кола представлена у виглядi окремих матриць, котрi пов'язанi мiж собою матрицею виткових зачеплень. Практична цттсть цифровог моделi електромагнтного кола полягае в тому, що параметри магнiтних юл вводяться в виглядi геометричних розмiрiв магнiтопроводiв. Це виключае необхiднiсть проводити еквiвалентнi перетворення для пiдготовки даних конк-ретног моделi. Це спрощуе вивчення вторинних джерел живлення та тших потужних споживачiв електроене-ргп за критерiями ефективностi, вагових i розмiрних параметрiв.
Наукова новизна. Топологiя елекmромагнimного кола представлена у виглядi окремих матриць, з'еднаних матрицею виткових зачеплень.
Практична цттсть. Параметри магнтних юл вводяться у виглядi геометричних розмiрiв магттних юл.
Ключовг слова: статичш електромагнтш пристрой топологгчно-гзоморфне моделювання; матричний топопоачний опис; топологгчнг матриц!; матрицi тциденцШ; блокова структура топологгчно1 матрицу, мат-риця виткових зачеплень; вторинн джерела живлення; генератори тпульсних сmрумiв.
МАТРИЧНО-ТОПОЛОГИЧЕСКАЯ МОДЕЛЬ ЭЛЕКТРОМАГНИТНЫХ
ЦЕПЕЙ
СИДЕЛЕВ НИ. канд. техн. наук, доцент, доцент кафедры автоматики и компьютерно-интегрированных технологий Черноморского национального университета имени Петра Могилы, Николаев, Украина, e-mail: [email protected];
ISSN1607-6761 (Print) ISSN 2521-6244 (Online)
«ЕЛЕКТРОТЕХШКА ТА ЕЛЕКТРОЕНЕРГЕТИКА» №1 (2018) Роздш «Електротехшка»
Цель работы. Разработать цифровую модель электромагнитных устройств для исследования и оптимизации мощных вторичных источников питания и электромагнитных преобразователей.
Методы исследования. Метод узловых потенциалов, метод контурных токов, топологически-изоморфные преобразования.
Полученные результаты. Современные системы автоматизированного проектирования требуют разработки специального математического обеспечения. Основными требованиями к разработке моделей могут быть наибольшая степень детализации, допустимое качество моделирования, простота получения параметров модели. В автономных электроэнергетических системах все устройства можно подразделить на три группы: источники питания, преобразователи и потребители электроэнергии. Среди них устройства, относящиеся ко второй группе, по массе и габаритам иногда соизмеримы с источниками питания и часто превышают по этим параметрам потребители электроэнергии. К тому же преобразователи электроэнергии также является своего рода потребителями энергии, которая используется для управления коммутационными элементами и выделяется в виде тепла. Количество потребленной энергии преобразователем влияет в худшую сторону на экономичность автономной системы и поэтому часто является критерием исследования. Целью данной работы является создание математического аппарата, позволяющего решать задачи моделирования и исследования электромагнитных устройств по частям (по видам накапливаемой энергии). Это позволит упростить исследование и оптимизацию таких технических характеристик, как коэффициент полезного действия, массогабаритные показатели и т.д. Предложенная математическая модель электромагнитных цепей имеет наибольшую степень детализации электрической и магнитной цепи. Магнитная цепь представлена так же подробно, что и электрическая цепь, и описывается контурной матрицей. Получено математическое описание электромагнитных устройств, в которых индуктивные параметры определяются геометрическими размерами и характеристиками магнитных цепей. Топология электрической цепи представлена матричными блоками, что позволило получить математическое описание, которое одновременно учитывает распределение токов и зарядов в элементах схемы. Система уравнений сводится к форме Коши и составлена относительно приращений магнитных потоков и потенциалов на конденсаторах, что упрощает ее решение численными методами на компьютере. Таким образом, удобно следить за энергетическими процессами в реактивных энергоемких элементах схемы. Разработана устойчивая и адаптивная цифровая модель электромагнитных схем, позволяющая объединить несколько методов интегрирования системы дифференциальных уравнений. Обратная связь предоставляется через специальный параметр. Это позволяет добиться максимальной корректности вычислений для энергетических компонентов при моделировании электромагнитной цепи. Оригинальность математического описания заключается в том, что топология электромагнитной цепи представлена в виде отдельных матриц, которые связаны между собой матрицей витковых зацеплений. Практическая ценность цифровой модели электромагнитной цепи заключается в том, параметры магнитных цепей вводятся в виде геометрических размеров магнитопроводов. Это исключает необходимость проводить эквивалентные преобразования для подготовки данных конкретной модели. Это упрощает изучение вторичных источников питания и других мощных потребителей электроэнергии по критериям эффективности, весовых и размерных параметров.
Научна новизна. Топология электромагнитной цепи представлена в виде отдельных матриц, соединенных матрицей витковых зацеплений.
Практическая ценность. Параметры магнитных цепей вводятся в виде геометрических размеров магнитных цепей.
Ключевые слова: статические электромагнитные устройства; топологически-изоморфное моделирование; матрично-топологическое описание; топологические матрицы; матрицы инциденций; блочная структура топологической матрицы; матрица витковых зацеплений; вторичные источники питания; генераторы импульсных токов.
