CC BY
66
UDC 681.513 + 621.365, 62-83-52
doi: 10.20998/2074-272X.2019.4.04
O.Y. Lozynskyi, A.O. Lozynskyi, Y.S. Paranchuk, R.Y. Paranchuk
SYNTHESIS AND ANALYSIS OF ARC FURNACE ELECTRICAL MODE CONTROL SYSTEM ON THE BASIS OF THREE-DIMENSIONAL PHASE CURRENTS VECTOR DISTRIBUTION
Goal. The purpose of the article is to create the method for the operative synthesis of an arc steel-melting furnace (ASF) electric mode (EM) control signal on the basis of a three-dimensional arc currents vector, which takes into account the stochastic nature of the processes in the melting space and power circuit and has low sensitivity to the control object parameters changes, as well as development of the control system structure for its implementation. Method. The basis of the created control method is formed on the statistical theory of dynamical systems, as well as the provisions of the statistical theory of optimal control based on the Fokker-Planck-Kolmogorov equation, which enables to synthesize operational control by the criterion of approaching the regulated coordinate distribution density to the Sfunction, that is to minimize the dispersion of the three-dimensional furnace phases arc currents vector. Results. The system of equations for operational real-time calculation of control influences of the thyristor switch of phase inductors, included in the power supply circuits of three-phase arcs, and the structural scheme of the adaptive contour for the formation of three-dimensional phase currents vector dispersion for the implementation of adaptive optimal control were obtained. Scientific novelty. For the first time, based on the Fokker-Planck-Kolmogorov equation, we obtain a system of equations representing a mathematical model of a stochastic adaptive optimal control of the arc furnace electric mode by the criterion of a minimum dispersion of three-dimensional phase (arcs) currents vector, which enables, in comparison with known methods, to increase dynamic precision of the arc currents stabilization at the level set by the criteria of energy efficiency and electromagnetic compatibility values. Practical value. The use of the proposed adaptive optimal control model and structural system scheme for its implementation allows, in comparison with the serial arc power regulators, to improve the dynamic accuracy of the arc current currents stabilization at the level of given optimal settings and, based on this, to improve the energy efficiency and electromagnetic compatibility indices of the arc furnace and power supply network. References 18, figures 5. Key words: arc furnace, electric mode, three-dimensional vector of phase currents, stochastic control, dispersion, optimization, adaptation, arc current control circuit.
Мета. Метою cmammi e створення методу оперативного синтезу сигналу керування електричним режимом (ЕР) дуговоТ сталеплавильно'1 печ1 (ДСП) на основ1 тривим1рного вектора струм1в фаз, що враховуе стохастичну природу процес1в у плавильному просторi, силовому електричному кол печ1, мае низьку чутлив1сть до змши параметры об'екта керування та розроблення структури системи керування для його реатзаци. Методика. В основу створеного методу керування покладено положення статистичноТ теори динaмiчних систем, а також положення статистичноТ теорп оптимального керування на основi рыняння Фоккера-Планка-Колмогорова, що дае змогу синтезувати оперативне керування за критерieм наближенням густини розподщу регульованоТ координати до S-функцИ, тобто мiнiмiзувaти диспераю тривимiрного вектора струмв фаз (дуг) дуговоТ печi Результати Отримано систему рiвнянь для оперативного в режимi on-line розрахунку керуючих aii.iitaia тиристорного комутатора фазних дросетв, що включен у силове коло живлення трифазних дуг, та структурну схему адаптивного контура формування дисперсИ' тривимiрного вектора струмы дуг дуговоТ печi для реатзаци адаптивного оптимального керування. Наукова новизна. Вперше на основi рiвняння Фоккера-Планка-Колмогорова отримано систему рiвнянь, що подають математичну модель стохастичного адаптивного оптимального керування електричним режимом дуговоТ сталеплавильноТ печi за критерieм мншуму дисперсИ' тривимiрного вектора струмв дуг, що дае змогу у порiвняннi з вЫомими методами пЫвищити динaмiчну точтсть стабШзащТ струмв дуг на рiвнi заданих за критерiями енергоефективностi та електромагштно'Т сумкно<т значень. Практична щншсть. Реал^защя запропонованоТмодел адаптивного оптимального керування та структурноТ схеми системи для ТТ реатзаци дасть змогу у порiвняннi з сершними регуляторами потужностi дуг полтшити динaмiчну точнкть стабШзаци струмы дуг на рiвнi заданих оптимальних уставок i на основi цього комплексно полпшити показники енергоефективностi та електромагштно'Т сумiсностiрежимiв дуговоТпч та електромережь Бiбл. 18, рис. 5.
Ключовi слова: дугова сталеплавильна ni4, електричний режим, тривимiрний вектор crpyMiB фаз, стохастичне керування, дисперая, оптимiзацiя, адаптащя, контур регулювання стрyмiв дуг.
Цель. Целью статьи является создание метода оперативного синтеза сигнала управления электрическим режимом (ЕР) дуговой сталеплавильной печи (ДСП) на основе трехмерного вектора токов фаз, который учитывает стохастическую природу процессов в плавильном пространстве, силовой электрической цепи печи, имеет низкую чувствительность к изменению параметров объекта управления и разработка структуры системы управления для его реализации. Методика. В основе созданного метода управления использованы положения статистической теории динамических систем, а также положение статистической теории оптимального управления на основе уравнения Фоккера-Планка-Колмогорова, что позволяет синтезировать оперативное управление по критерию приближения плотности распределения регулируемой координаты к 8-функции, то есть минимизировать дисперсию трехмерного вектора токов дуг дуговой печи. Результаты. Получена система уравнений для оперативного в режиме on-line расчета управляющих воздействий тиристорного коммутатора фазных дросселей, включенных в силовую цепь питания трехфазных дуг, и структурную схему адаптивного контура формирования дисперсии трехмерного вектора токов фаз дуговой печи для реализации адаптивного оптимального управления. Научная новизна. Впервые на основе уравнения Фоккера-Планка-Колмогорова получена система уравнений, представляющих математическую модель стохастического адаптивного оптимального управления электрическим режимом дуговой сталеплавильной печи по
© O.Y. Lozynskyi, A.O. Lozynskyi, Y.S. Paranchuk, R.Y. Paranchuk
критерию минимума дисперсии трехмерного вектора токов дуг, что позволяет по сравнению с известными методами повысить динамическую точность стабилизации токов дуг на уровне заданных по критериям энергоэффективности и электромагнитной совместимости значений. Практическая ценность. Реализация предложенной модели адаптивного оптимального управления и структурной схемы системы для ее реализации позволяет по сравнению с серийными регуляторами мощности дуг улучшить динамическую точность стабилизации токов дуг на уровне заданных оптимальных уставок и на основе этого комплексно улучшить показатели энергоэффективности и электромагнитной совместимости режимов дуговой печи и электросети. Библ. 18, рис. 5. Ключевые слова: дуговая сталеплавильная печь, электрический режим, трехмерный вектор токов фаз, стохастическое управление, дисперсия, оптимизация, адаптация, контур регулирования токов дуг.
Introduction. Arc furnaces are powerful electrical technological installations that belong to a class of complex systems and are characterized by extremely incident, dynamic, nonlinear, phase by phase non-symmetric nature of loading and continuous action of intense coordinate and parametric disturbances in arc gaps and power supply circuit of three-phase arcs. The specified loading characteristics complicate the control process of such objects and impose appropriate limitations on the system engineering - models, methods and approaches for improving existing systems of control modes and regulation of electric coordinates.
The problem of integrated improvement of energy efficiency and electromagnetic compatibility indices of arc steel-melting furnaces (ASFs) is dictated by the necessity of increasing the competitiveness of electric steels and high alloys on the domestic and foreign markets of metal products. Its state is largely determined by the level of excellence of automatic control systems (ACS) of electrical modes (EM) of arc furnaces, characterized by the speed and dynamic accuracy of coordinate (first of all, arc currents) regulation, the efficiency of electric modes control models, state identification and prediction of the process of electric steelmaking.
It is obvious that for such electrical technological nonlinear stochastic control objects with a range of installed power of power electric equipment of 1...175 MVA, it is most appropriate to use models based on probabilistic characteristics of control processes and that most fully correspond to the nature of the processes which occur in them. Efficiency and completeness of solution of control problems in general and adaptive optimal, in particular, of electrical modes of electric steelmaking are determined, first of all, by the speed and dynamic accuracy of the regulation of the coordinates of the electric mode and, the most important, by currents of arcs.
Problem definition. The complexity of the tasks of the modern theory of adaptive optimal control of stochastic dynamic objects and processes requires the improvement of the mathematical apparatus for their description, identification and control synthesis models, and also requires significant computing power of digital means (microcontrollers, microprocessor devices) for implementation of identification and real-time control.
Unfortunately, the control of modes in the vast majority of arc furnaces is realized on the basis of classical deterministic models of identification of states, parameters, phase by phase regulation of coordinates and control of modes that do not correspond to the stochastic phase by phase interconnected nature of the processes underlying them functioning.
Therefore, in our opinion, the most appropriate approach to solving the above problem is to improve the existing and create new effective methods and approaches for the tasks of controlling modes and coordinate regulation, in particular, the operational formation of control effects based on the three-dimensional vector of phase currents and its probabilistic characteristics.
Review of recent publications. For the first time, the theoretical principles of stochastic models for the formation of control effects for the electromechanical system of displacement of electrodes of three-phase ASFs have been published in [1]. The coefficients of the relationship between the average rectified currents of arcs included in the created model of the ASF are proposed to be determined on the basis of probabilistic analysis and taking into account the correlation interconnections of phase processes in the furnace space of the three-phase arc furnace. The model of electric mode control of the ASF obtained in this work allows to adjust the control signals of typical regulators of the electric mode by current and, thus, to avoid false operations of the regulator of the /'-th phase (to eliminate false electrodes displacement) by perturbations in adjacent phases. Under this control model, the phase by phase autonomy of electric mode control is improved, which, in turn, increases the efficiency of the electrometallurgical plant. But with the help of such a model, it is impossible to implement the above-stated approach of control of three-phase ASF modes.
Similarly, the problems of the phase by phase autonomy of electric mode control by arc furnace phases on the basis of the account of stochastic parameters of perturbations of adjacent phases are considered in the work [2]. According to this work, the signal formed by the displacement of the electrodes in each phase is formed additively from the signals of the dissonance of all three phases, which are normalized by the weight factors. To find the values of these coefficients, the authors obtained a mathematical model that describes the reactions of such a complex object as the ASF on the pre-synthesized control influences and perturbation processes that convert the electric mode of the furnace to a particular state. In the final case, these coefficients are also some averaged over a certain time interval the coefficients of the weight of the signals of discrepancies of phase modes. Nevertheless, it should be noted here that in this paper for the first time it was noted that optimization of the system of controlling the electric mode of the ASF should be carried out with an orientation to such an integral characteristic of the mode as the variance of the three-dimensional vector of arc currents, but specific solutions in this direction in this work are not offered.
In the paper [3] for the analysis and synthesis of control systems of electric drives, which are under the influence of random perturbations, the expediency of using the probabilistic approach is substantiated. This approach ensures an adequate response of the ACS to the processes and perturbations that take place in the control object. In this paper, the method of formation of the variance of the regulated coordinate, that is, the method of stochastic dynamic stabilization for application in the problems of stabilization of coordinates of electromechanical systems with random perturbations is worked out. In its implementation, it is possible to control the dispersion of the output coordinate of the dynamic system in accordance with the current requirements and conditions of operation of the control object, in particular such as an arc furnace.
In work [4] a mathematical description of the electromechanical control system of electrodes position was worked out, on the basis of which «on-line» synthesis of the desired dynamic characteristics of the state change process on the basis of integral quality criteria is carried out. The proposed approach is based on the account of the phase by phase interconnection of electric modes caused by the peculiarities of the parametric non-symmetric power supply circuit of the three-phase arcs and the impedance or differential law of the formation of a control signal on the displacement of the electrodes. But this work, like the previous one, is far from the idea of applying the control of processes in arc furnaces and the idea of forming a three-dimensional vector of arc currents and its stochastic characteristics.
The mathematical and computer models of control of electrical modes of the three-phase arc furnaces proposed in the papers [5-9] have certain advantages and disadvantages among themselves in terms of the completeness and accuracy of the description of modes, identification of states, conveniences in use and readjustment, but they do not meet the above requirements for energy efficiency and electromagnetic compatibility modes in full, primarily due to the inadequacy of the control models to the nature of the real processes of changing the coordinates of the EM, also because of the high sensitivity of the obtained dynamics to changing parameters of the control object, which is highly undesirable in terms of continuous action of stochastic parametric disturbances in the power circuit of the furnace and in phase arc gaps.
The goal of the work is the development of system and structural solutions for the operative synthesis of the control signal of an arc furnace electric mode on the basis of a three-dimensional vector of arc currents, taking into account the stochastic nature of the processes in the melting space and the power supply circuit of the three-phase arcs and has a low sensitivity to the change of the control object parameters, the use of which in comparison with known solutions allows to increase the dynamic accuracy of the stabilization of arc currents at the level specified by the criteria of energy efficiency and electromagnetic compatibility.
The scientific task is to create the method of adaptive optimal control of the electric mode of an arc furnace according to the criterion of the minimum
distribution density of a three-dimensional vector of phase currents and the structural scheme of the control system for its realization, which, in comparison with known control methods, enables to increase the dynamic accuracy of the stabilization of arc currents at the given (in particular, optimal) values and thereby improve the energy efficiency and electromagnetic compatibility indicators of arc furnace and power supply network.
Identification of previously unsolved parts of the general problem. The dynamics of the regulation of the EM coordinates of the ASF, in particular the arc currents, in the overwhelming majority of the existing (serial, typical) control systems of the process of electrical steel molding (by using arc power) does not fully meet the high modern requirements for energy efficiency and electromagnetic compatibility indicators of arc furnace modes and power supply network. Therefore, the problem of creating effective system, circuit and algorithmic solutions, which are based on increasing the speed of processes of regulation of phase (arcs) currents and complex improvement on the basis of indicators of energy efficiency and electromagnetic compatibility of modes for today for the electro-metallurgical industry is important and relevant.
Content of research material. In this paper, the term «optimality» is used in the narrow sense, in which the system of automatic control is evaluated only by the indicators of the quality of dynamic processes, and one of the criteria of this quality is the integral quality index -the generalized dispersion of the regulated coordinate. Such a description of the quality criteria makes it possible to use modern well-developed mathematical optimization apparatus to find optimal control.
Under the action of the flows of random perturbations and control influences, the electric mode (which is estimated by the stochastic characteristics of arc currents) of a three-phase arc furnace can be in different states.
Define these states by a three-dimensional vector of phase currents, which characterizes the electric mode of the arc furnace during the melting company as:
1 - a state characterized by the given (required) values of the coordinates of the electric mode, in particular by the set values of arc currents;
2 - a state characterized by deviations of arc currents in the domain of admissible directive deviations;
3 - a state characterized by deviations of the arc currents into the region of large, in particular extreme, deviations - short circuits, breaks of arcs and close to them.
It is clear that the concept of being in one or another state is associated with some definite time interval T0. The graph of the states of the three-dimensional vector of states of the electric mode of an arc furnace is shown in Fig. 1.
Each of these states is proposed to be identified by means of the integral parameter (indicator) value, which is the generalized variance of the three-dimensional vector of the ASF arc currents. In turn, this integral indicator is defined as the determinant of the matrix of second order moments, or the correlation matrix [10], and it characterizes the value of the deviation of the vector of
arc currents of the three-phase arc furnace from the desired state.
Fig. 1. A graph of states of the electric mode of an arc furnace -kiM M*)> ku(i), hM ki,3(t), k2)3(f), k3,3(i), k3,i(i), k3,2(i) -the intensities of the transition from state to state
So, having adopted the designation of a three-dimensional vector of regulated coordinates (for the considered problem of arc currents) as
Ia = y = |>1 >"2 Уз\,
we obtain an expression for the dispersion of the three-dimensional vector in the form:
DI = D> = det A> ,
where Ay is the correlation matrix, or matrix of moments of the second order of the form:
¿11 ¿12
A> =
¿13 ¿23 ¿33
¿21 ¿22 ¿31 ¿32
The criterion for the functioning of such a dynamic system is the desired level of probability of being the system in state 1 during the time of the melting of charge. In view of the above, we write the expression for the density of the three-dimensional vector of the regulated coordinate as:
PO'U >2, >3) =- 1
(2/^)3/2^detA~>
x exp
1 3 I
-,A>UI \yi j=1
• >j- yj\
(1)
where A>. =
1
det A,
- • Ai j are the elements of the inverse
state 1. Taking into account that in a real object we have a flow of disturbances and a flow of control influences that change the state of the system (the state of the electric mode), we write the model of state dynamics in the form of the Kolmogorov-Chapman equation [11]: 3 3
■ = Z Pji • kji (t) -X pi • h j (t),
j=1 j=1
dp (t) dt
matrix Ay ; Aj are the corresponding elements of the
attached matrix, and the probability of being of the electric mode of an arc furnace in a given state is found as an integral of density (1):
>1* >* >3
P(У1, >2, >3) = jj jp^ >2, >3)d>1 d>2 d>3 , 0 0 0
where >*, >*, >* are the maximum values of the change of the coordinates of the electric mode, for example, the value of the circuits of the short circuit in the ASF phases.
It is clear that the smaller the variance of the three-dimensional vector of the regulated coordinate (the arc currents of the arc furnace), the greater the probability of a state that is identified by such a variance value, i.e.,
and the system of equations for determining the probabilities of individual states is written as:
dp1(t) = -P (t) • [[ (t) + ¿13 (t)] + k21(t) • P2 (t) +
dt
+ k31(t) • P,(t);
dp2(t) = - P2 (t) • [[1 (t) + ¿23 (t)] + ¿12 (t) • P1 (t); (2)
dt
dp3(t) = -P3 (t) • [[31 (t) + ¿32 (t)] + ¿13 (t) • P1 (t) . dt
From the presented system of equations (2) and graph in Fig. 1 it can be concluded that the intensities of transitions ¿21(t); ¿31(t); £32(t) are formed by controlling effects of the control system, which bring the electric mode of the arc furnace from the unwanted states 2 and 3, and especially from the state 3 to state 1. At the same time, the intensities of transitions £12(t); ¿13(t); ¿23(t) are determined by the disturbances that act in the melting space and power circuit of the ASF.
The method of formation of effective control influences was developed in works [12-15], namely, the method of forming a vector of control effects of the so-called second (high-speed electric) circuit of regulating arc currents included in the structure of existing (serial) one-contour electric control systems of electric mode of an arc furnace. The main feature of such a two-contour structure of the electric mode control system of the arc furnace is the high speed of the arc currents regulation (the current regulation time is 0.03-0.04 s), which allows to significantly improve the control dynamics, that is, to obtain a high dynamic accuracy of the arc currents stabilization, and thus, with high accuracy to control the dynamics of the state graph, in particular, to transfer the electrical mode of the arc furnace to states 1 or 2.
The functional scheme of such a two-circuit control system of the EM of the arc furnace is shown in Fig. 2. This system contains the traditional electromechanical (or electrohydraulic) contour for regulating the lengths of the arcs EMCRLA (it is traditionally called the regulator of the arcs power), which has a relatively high inertia, which results in the operation of such a regulator of arcs power is accompanied by a significant dispersion of the coordinates of the EM - lengths, voltages, currents and powers of arcs.
Each phase channel of this circuit includes the arcs voltage sensor (VS) and current sensors (CS), the control signal formation unit (CSFU), the electrode displacement electric drive (EDED) and the electrode displacement mechanism (EDM), and also contains a voltage level switch (VLS) of the furnace transformer FT. Most often,
in such EMCRLAs, the control signal on the displacement of the electrode in each phase is formed according to the differential law.
controllability and dynamic precision of the stabilization of currents (lengths, voltages, powers) of arcs under the conditions of the action of continuous random transient parametric and coordinate disturbances.
Ia, kA laB, kA
IaC, kA
Ia, kA I.B, kA IaC, kA
Fig. 2. Functional scheme of the two-circuit control system of the electric mode of the arc furnace
An example of the dynamics of the regulation of arc currents by the electromechanical circuit, performed in the structure of the ACS of the EM of the arc furnace ^Cn-200 by the serial regulator of the arcs power of the type APflM-T-12, is illustrated in Fig. 3,a by the fragments of the temporal dependencies of the arcs currents Iaj(t) (j=A, B, C) (computer experiments were carried out at time intervals of stationary Tc=180-300 s of random processes of perturbations by the lengths of arcs of the furnace ^Cn-200 for various technological melting stages).
In the composition of each phase channel of the electric high-speed contour for the regulation of arc current (HSCRAC) there are the arc current sensor (CS), adjuster (CA) and the regulator (CR) of the arc current, the thyristor switch (TS), the control effect of which is the time of the shunting of the throttle (Th) on a certain regulated part of the half-period of the voltage given by the angle a of the control of the parallel thyristor-reactor group. At the joint functioning of these two regulating contour, their advantages are combined: reliable ignition of arcs in the treatment of extreme disturbances -operational short circuits and breakdowns of the arcs of the EMCRLA and high speed of regulation of arc current of the HSCRAC. As a result of such a combination of contours in the two-contour structure of the ACS of the EM of the ASF it is possible to achieve high
Fig. 3. Temporal dependencies of arc currents of the flCn-200 furnace at the operation of the regulator AP^M-T-12 (a) and AP^M-T-12 with a high-speed contour at control by the criterion DQ min (b)
To compare dynamics, as an example, Fig. 3,b shows the temporal dependencies of arcs currents in the operation of a two-contour system under the influence of random disturbances with the same parameters of stochastic perturbation characteristics using the proportional integral current regulator CR and EM optimization based on the scalar criterion for the minimum of the dispersion of the reactive power of the furnace DQ ^ min, which is largely correlated with the criterion for the minimum dispersion of arcs currents. The dispersion of arcs currents in the two-contour structure of the ACS of the EM of the ^Cn-200 at the functioning of the proportional-integral regulator was: dIa = 2.95-107 A2;
Dj = 3.15-107 A2; Dj = 2.35-107 A2. The average on
aB 1ac °
the phases the dispersion of the currents of arcs was Dia =2.82-107 A2. The presented temporal dependencies illustrate high speed of regulation and high-quality dynamic stabilization of arcs currents in the structure of the two-contour ACS of the EM of the arc furnace. Comparison of dispersions shown in Fig. 3 for processes of the change of arcs currents Ijt) with the values of dispersions under other parameters of stochastic disturbances by the arcs lengths shows a considerable (almost an order of magnitude) decrease in the dispersion of arcs currents in the joint operation of electromechanical and high-speed electric contours with a proportional integral regulator of arcs currents in the two-contour structure of the ACS of the EM of the arc furnace in comparison with the operation of a one-contour ACS of the EM (regulator APflM-T-12) [10].
For such a two-contour structure of the ACS of the EM of the arc furnace, it is important to develop effective
system-engineering solutions - models of the formation of phase controlling effects of the high-speed contour of regulating the arcs currents of the HSCRAC, which would improve the indicators of the electrical efficiency and electromagnetic compatibility of the arc furnace modes and power grid.
Therefore, it is expedient in the context of the task defined in the work to develop the theoretical bases of a stochastic model of operative synthesis of control influences by a three-dimensional vector of arcs currents in the structure of such a two-contour ACS of the electric mode of the ASF.
For the synthesis of the vector of control effects on the regulation of the arcs currents of the two-contour ACS of the EM of the arc furnace, we apply the principle according to which the optimal change in the distribution density of the three-dimensional arcs currents vector is realized, in particular, in the direction of approximation of this distribution to the form of the ¿-function.
So, in the ideal case, the system will better perform its intended purpose, namely, to ensure the state 1 of the three-dimensional vector of phases currents, the sooner it will convert the initial distribution density p(>A, >B, >c, t0) of the three-dimensional vector of currents of arcs of the furnaces in ¿-function, or in the ¿-distribution, which is combined with the point >a = >A.set, >b = >B.set, >c = >c.set, where >a, >b, >c, >A.set, >B.set, >c.set are the current and set values of the regulated coordinates - arcs currents of the arc furnace.
This conclusion corresponds to the principle of the statistical theory of transients, according to which, on the basis of the Fokker-Planck-Kolmogorov equation [16], optimal control is obtained in relation to the ctransient of changing the density of the «-dimensional probability distribution to the density of the form of the ¿¡-function.
For our case, we write the Fokker-Planck-Kolmogorov equation in the form
dp ^ It
i=1
d (P • Fi) d>i
(3)
Such a control on one of the coordinates, for example, on the first, for our case, is written as:
«1 = d1 • sign(8ln(P(>1, >2, >3))).
8>1
To find the appropriate control we write down the distribution density (1) of the three-dimensional vector of the regulated coordinate of the electric mode of the ASF in the form convenient for differentiation:
P^b >2, >3) =
1
(2/^)3/2^det A>
r • exp
11
2 det A
>
x{A11(>1 ->1) + 2 • A12(>i ->1) • (>2 ->2) + + 2 • A13(>1 ->1) • (>3 ->3) + A22 (>2 ->2)2 + + 2 • A32 (>2 ->2) • (>3 ->3) + A33 (>3 ->3)2}
where >1 , > 2 , > 3 are the mathematical expectations of
phase regulated coordinates, which for the considered object are given by the values of the currents of arcs of individual phases of the ASF.
If we make the necessary mathematical transformations and substitute currents of arcs IaA, IaB, laC instead of generalized regulated coordinates > (i = A, B, C), we obtain the expressions for the mentioned control influences uA, uB, uc for each of the phases of the arc furnace in form:
«A = d1 • sign{ '-• [An(IaA - IaA ) +
>/detAI
+ A12(IaB - IaB ) + An(I aC - IaC )]};
UB = d2 • sign{
1
r—T [A21(IaA - IaA ) + VdetAI (4)
uc = d3 • sign{'
+ A22(IaB - IaB ) + A23(IaC - IaC )]}; 1
yjdetA
•[A31 (IaA - IaA ) +
I
as an equation for a dynamic system in the absence of noise.
In this equation (3) it is indicated: p is the probability distribution density of a three-dimensional vector p = P(>1, >2, >3); >i + Fi(>1, >2, >3) = 0 is the equation describing the system coordinates movement.
For the automatic control system, the functions F can be represented as a set of two functions:
Fi = fi (>b >2, >3) -u(Уl, >2, >3) , where the function f belongs to the control object, and the function u -to the control system.
Synthesized in [16] controls for dynamical systems that provide an optimal change in the density of the distribution of the «-dimensional vector of regulated coordinates >t are calculated by the expression:
. 8 ln( p) u, = di • sign(—-),
8>i
where dt is the maximum permissible or limit value of control influence; >t is the component of the «-dimensional vector of regulated coordinates.
+ A32(IaB - IaB ) + A33(Ic - IaC )]}.
As we see, for the operative formation of the control influences uA(t), uB(t), and uc(t) it is necessary to know the matrix Ai of other moments of the three-dimensional vector of the currents of arcs (phases) of the ASF. The operational calculation of control effects by the received model (4) for modern microprocessor devices is a simple (to some extent trivial) technical task.
We note here that the control signals obtained by the equations (4) are formulated as boundary controls and can only be considered as conventionally optimal ones. However, studies [16] show that these controls can be quite close, and even coincide with strictly optimal controls, which translate the «-dimensional distribution of the regulated coordinate into ¿-function.
The obtained model (4) of operative synthesis of the vector of control signals uA, uB, uc is realized in the proposed structure of the two-contour adaptive control system of arc furnace electric mode, the functional block diagram of which is shown in Fig. 4.
In the presented scheme of the two-contour adaptive ACS by the values of implementations of the average
x
rectified currents of the three phases jaA, jaB, jac of the arc furnace in the block of calculations BC operatively at each interval of stationary (Tc = 3-5 min) of the processes of phases currents changes, the calculation of the matrix Aj of other moments of the three-dimensional vector of
the currents of arcs of arc furnace, its determinant deA and the values of the elements Ati of the adjoint matrix A of the inverse matrix Al _1 is performed.
r"
A„. Ai.-Aa.detA
An, An, Aa,detA
4s. Aa- As. det A
bc
VV
da-
L/l^ V
gcsa 11*
gcsb "B
gcsc Uc
Ifl
edm it) eded csfu
asf
Fig. 4. Functional scheme of the two-contour adaptive control system of arc furnace electric mode for minimization of dispersion of three-dimensional vector of currents
The specified signals from the BC outputs are fed to the third vector input of each phase generator of the control signal for current Iai, namely: Au, A12, A13, detAj - to the third GCSA input; A21, A22, A23, detAj - to the third GCSb input; A31, A32, A33, detAj - to the third GCSC input, and to the first and second vector inputs of GCSA, GCSb, GCSc the three-dimensional vector of currentaverage rectified values of arcs currents Iai from the outputs of currents censors CS and a three-dimensional vector of the averaged on the stationary intervals Tc currents I a/ are input, respectively. The fourth input of the GCS; is fed by the scalar signal d of the maximum value of the control influence. On the outputs of the GCS' continuously online phase control signals u are generated based on the obtained control model (4),
which are fed to the corresponding phase thyristor switches TS''. Output signals at of pulse-phase control systems (PPCSs) determine the moments of shunting/ switching on in the power circuit of the corresponding phase throttles Th'.
The proposed procedure for operative synthesis of control effects uA(t), uB(t), and uc{t) is divided into two parallel processes in time: the first process implements the adaptation of the coefficients of the synthesis model (4) to change the parameters of stochastic characteristics of coordinate and parametric perturbations on the time the interval of melting, and the second one performs online operative with discreteness At calculation of control effects uA(t), uB(t), and uc(t) in GCSA, GCSB, and GCSC units by this model. The first process involves at each interval of stationary Tc, the duration of which is correlated with the technological stages of melting (the physical and chemical state of the melt) and depends on the types of arc furnaces ((Tc = 180-300 s), the calculation in the BC of the coefficients of the model (4): det AI,
I
a, j, Anm, and the second one implements the process of calculating the control effects uA, uB, uc and is performed with the interval At = 0.02 s in the function of averaging the values (rms) of phases currents Ij (tk)
(tk = tk-l +At, m = 1, 2, 3, n = 1, 2, 3, j = A, B, C ) in this interval.
To verify the effectiveness of the proposed stochastic model for synthesizing the control signals uA(t), uB(t), and uc(t), the corresponding mathematical experiments were performed on the three-phase in instantaneous coordinates Simulink model [17, 18] of the two-contour ACS of the EM of the arc furnace ^Cn-200 (Fig. 2). The modelling studies were performed at the functioning of the proportional integral arc voltage regulator and using the proposed stochastic model (4) for synthesizing the control signals uA(t), uB(t), and uc(t) which is implemented in the proposed structure of the adaptive ACS of the EM of the arc furnace (Fig. 4), which implements a strategy for minimizing the dispersion of a three-dimensional vector of phases currents.
To do this, a computational block is introduced into the structure of the Simulink model, in which, according to the well-known model [10], on the intervals of stationary Tc/, the coefficients of the model (4) det AI,
I
An
are calculated. In the next interval of
a,j' ^n,m
stationary Tc/+1, these coefficients are used in GCSA,
GCSB, and GCSc for synthesis of control signals uA, uB, uc by model (4) in the function of current averaged over a period of supply voltage At = 0.02 s phase currents I a , IB , Ic. Obtained control uA, uB, uc are applied to the inputs of thyristor switches. At the same time, the new values of the coefficients of the model (4) det Al, la, j, An,m , are used in the next Tc, +2 interval of online control signal uA, uB, uc, etc. synthesis, are calculated on the current stationary interval Tc, +2 .
Figure 5 shows the initial fragments of temporal dependencies of phase currents obtained on the Simulink
L„, R
E
j
Ps
l
emcrla
model at the control by the above model (4), implementing the strategy of minimizing the dispersion of a three-dimensional vector of phases currents.
laA, kA
Ib, kA
mNM hffi
If II 1 - It a jc ;>; i JP
Fig. 5. Temporal dependencies of the arcs currents of the furnace flCn-200 with the regulation by the obtained model (4) of minimizing the dispersion of the three-dimensional vector of arcs currents in the two-contour structure of the ACS of the EM of the arc furnace ^Cn-200
As a result of working out the time dependenceis of the currents of the arcs I(t) in Fig. 5, the following values of their dispersions are obtained: D1a = 0.95-107 A2;
dIb = 1.05-107 A2; Dlc = 1.03-107 A2. The average by
the phases the dispersion of the currents of arcs is
di = 1.01-107 A2.
Comparative analysis of temporal dependencies of arcs currents in Fig. 3,b and Fig. 5 shows that the average by phases dispersion of arcs currents at the regulation by the obtained model (4) of minimizing the dispersion of the three-dimensional vector of arcs current in the two-contour structure of the ACS of the EM decreased by 2.73 times. The computer researches for perturbation processes by the lengths of arcs of other technological melting stages, which differ in frequency spectrum and amplitude of oscillations, have shown that the average dispersion of arcs currents at the regulation by the obtained model (4) of the minimization the dispersion of the three-dimensional arcs current vector in the two-contour structure of the ACS of the EM compared with the use of known proportionalintegral model of the formation of control signals uA(t), uB(t), and uc(t) decreased in 1.6-3 times.
Due to the cyclical updating of the matrix AI of other moments of the three-dimensional vector of arcs currents of the acr furnace in the block BC and elements of the vector Iai on the outputs of the averaging units AU, the adaptation of the vector of control influences uA(t), uB(t), and uc(t) to the change of parameters of stochastic characteristics of coordinate and parametric perturbations in the power circuit and arc gaps of the arc furnaces in the full range of melting is realized.
Thus, due to the obtained model (4) of the operative formation and adaptation of the vector of control effects and their realization through the adaptive contour of the formation of the dispersion of the three-dimensional arcs current vector ACFDTCV of the arc furnace one can obtain a significant approximation of the distribution of the three-dimensional vector of arcs phase currents of the arc furnace to the type of the ¿-function. And the fact that the reduction of the dispersion of arcs currents
significantly influences the improvement of energy efficiency indicators (electrical efficiency coefficient of the furnace, specific energy consumption, specific ASF productivity, the price of the ton of smelted steel, etc.) of an electro-technological installation, which is an AC arc furnace power is a well-known fact.
It should also be noted that the proposed stochastic model of the regulation of the coordinates of the EM of the arc furnace and the model of adaptive synthesis of control signals based on the three-dimensional vector of phase currents implements a relay, in contrast to the existing control law. Such a law of control of compliance with the conditions of stability, as is known, provides maximum regulation speed and, as a result, high dynamic accuracy of the arc current stabilization at the level of optimal for the selected criterion of values, and also is characterized by a much lower sensitivity to the change of parameters of the control object - the parameters of dynamic volt-ampere characteristics of three-phase arcs and parameters of the elements of the power circuit (short circuit) of an arc furnace. The last feature of the proposed method is particularly important for implementing adaptive optimal control strategies under continuous conditions of intense stochastic parametric perturbations in the power circuit and arc gaps of the furnace during melting.
The stochastic model of operative synthesis and adaptation of the vector of controlling for minimizing the dispersion of the three-dimensional arcs currents vector obtained in the paper is expedient for practical use in the two-contour structures of the ACS of the EM of arc furnaces with a high-speed contour for regulating arc currents.
Conclusions. The developed in the paper the theoretical bases of the method of electric arc furnace control on the basis of the formation of the distribution density of the three-dimensional vector of arc currents make it possible to implement the adaptive optimal control strategy of the electric mode by the criterion of the minimum dispersion of arc currents. As the results of the performed model researches show, the synthesis of the vector of control of arc currents by the obtained stochastic model (4) compared with the use of the proportional-integral regulator of arc current for controlling the EM reduces the dispersion of arc currents by 1.6-3 times. The minimization of the arc current dispersion positively affects the reduction of the power of electrical losses in the elements of the short arc furnace network and, as a result, reduces the specific energy consumption, the corresponding increase in the furnace productivity and the electric efficiency of the arc furnace. In addition, the transformation of the three-dimensional vector of arc currents into the ¿-function in the control of the electric mode greatly reduces the consumption of reactive power of the furnace and, accordingly, increases the power factor, reduces the volatility and deviation of the network voltage on the power supply buses of the furnace and accordingly reduces the dose of the flicker.
REFERENCES
1. Lozinskyi O.Yu., Maruschak Y.Yu. Three-dimensional stochastic model of the electric mode regulation of an arc
furnace. Industrial Process Automation in Engineering and Instrumentation, 1993, no.31, pp. 7-11. (Ukr).
2. Lozinsky O.Yu., Parancuk Ya.S., Lozinsky A.O. Optimization of dynamic regimes of interconnected electromechanical systems. Bulletin of Lviv Polytechnic National University, «Electric Power and Electromechanical Systems» series, 2001, no.421, pp. 98-103. (Ukr).
3. Lozynskyi O.Yu., Parancuk Ya.S., Moroz V.I. Synthesis of the control process for electromechanical systems that are under the influence of random perturbations. Bulletin of NTU «KhPI». Series: Problems of automated electric drive. Theory and practice, 1994, pp. 104-106. (Ukr).
4. Lozinsky O.Yu., Parancuk Ya.S., Tsyapa V.B. Mathematical description of the dynamics of the regulation of the position of the electrodes of the chipboard model in the space of states. Bulletin of Lviv Polytechnic National University, «Electric Power and Electromechanical Systems» series, 2017, no.840, pp. 54-60. (Ukr).
5. Nikolaev A., Povelitsa E., Kornilov G., Anufriev A. Research and Development of Automatic Control System for Electric Arc Furnace Electrode Positioning. Applied Mechanics and Materials, 2015, vol.785, pp. 707-713. doi: 10.4028/www.scientific.net/amm.785.707.
6. Ghiormez L., Panoiu M. Current control of a 3-phase electric arc furnace using fuzzy logic. ANNALS of Faculty Engineering Hunedoara - International Journal of Engineering, 2015, Fascicule 4 Tome XIII, pp. 237-242.
7. Nikolaev A.A., Tulupov P.G. Method of setting optimum asymmetric mode of operation of electric arc furnace. 2016 11th France-Japan & 9th Europe-Asia congress on Mechatronics (MEcATRONicS) / 17th International conference on Research and Education in Mechatronics (REM), Jun. 2016. doi: 10.1109/mecatronics.2016.7547111.
8. Zheng T., Makram E.B. An adaptive arc furnace model. IEEE Transactions on Power Delivery, 2000, vol.15, no.3, pp. 931-939. doi: 10.1109/61.871355.
9. Hooshmand R., Banejad M., Torabian Esfahani M. A New Time Domain Model for Electric Arc Furnace. Journal of Electrical Engineering, 2008, vol.59, no.4, pp. 195-202.
10. Pugachov V.S. Teoriia sluchainykh funktsii i ee primenenie k zadacham avtomaticheskogo upravleniia [The theory of random functions and its application to problems of automatic control]. Moscow, Fizmatizdat Publ., 1960. 883 p. (Rus).
11. Kazakov V.A. Vvedenie v teoriiu markovskikh protsessov i nekotorye radiotekhnicheskie zadachi [Introduction to the theory of Markov's processes and some radio engineering problems]. Moscow, Soviet radio Publ., 1973. 232 p. (Rus).
12. Lozynskyy O., Lozynskyy A., Paranchuk Y., Paranchuk R., Marushchak Y., Malyar A. Analysis and Synthesis of Intelligent System for Electric Mode Control in Electric Arc Furnace. Part of the Lecture Notes in Electrical Engineering book series (vol.452). Analysis and Simulation of Electrical and Computer Systems, 2017, pp. 111-130. doi: 10.1007/978-3-319-63949-9_7.
13. Lozynskyi O., Lozynskyi A., Paranchuk Y., Paranchuk R., Holovach I., Tsyapa V. Fuzzy extreme control and electric mode coordinates stabilization of arc steel-melting furnace. 2016 Xlth International Scientific and Technical Conference Computer Sciences and Information Technologies (CSIT), Sep. 2016. doi: 10.1109/stc-csit.2016.7589866.
14. Lozynskyy O., Paranchuk Y., Paranchuk R. Fuzzy control law of electrode travel in arc steelmaking furnace. 16th International Conference on Computational Problems of Electrical Engineering (CPEE), Sep. 2015. doi: 10.1109/cpee.2015.7333349.
15. Lozinskyy O.Yu., Paranchuk Y.S. Optimization of the modes of the process control system for electric steel-melting in arc steel-smelting furnaces. Electrical engineering, 2004, no.6, pp. 50-54. (Rus).
16. Krasovskyi A.A. Statisticheskaia teoriia perekhodnykh protsessov v sistemakh upravleniia [Statistical theory of transient processes in control systems]. Moscow, Nauka Publ., 1968. 240 p. (Rus).
17. Lozynskyi O.Y., Paranchuk Y.S., Paranchuk R.Y., Matico F.D. Development of methods and means of computer simulation for studying arc furnace electric modes. Electrical engineering & electromechanics, 2018, no.3, pp. 28-36. doi: 10.20998/2074-272X.2018.3.04.
18. Lozynskyy O., Paranchuk Y., Stakhiv P. The Study of Dynamics of the Two-Loop Arc Furnace Electric Mode ACS on a Simulink-model. Przeglqd Elektrotechniczny, 2018, vol.1, no.12, pp. 24-27. doi: 10.15199/48.2018.12.06.
Received 30.04.2019
O.Y. Lozynskyi1, Doctor of Technical Science, Professor, A.O. Lozynskyi1, Doctor of Technical Science, Professor, Y.S. Paranchuk1, Doctor of Technical Science, Professor, R.Y. Paranchuk1, Candidate of Technical Science, 1 Lviv Polytechnic National University, 12, S. Bandera Str., Lviv, 79013, Ukraine, e-mail: yparanchuk@yachoo.com
How to cite this article:
Lozynskyi O.Y., Lozynskyi A.O., Paranchuk Y.S., Paranchuk R.Y. Synthesis and analysis of arc furnace electrical mode control system on the basis of three-dimensional phase currents vector distribution. Electrical engineering & electromechanics, 2019, no.4, pp. 26-34. doi: 10.20998/2074-272X.2019.4.04.
