CC BY-NC
12
DOI: 10.24937/2542-2324-2019-2-S-I-198-204 УДК 621.313.12
H.A, Доброскок, О.В. Мохова
Филиал «ЦНИИ СЭТ» ФГУП «КПЩ», Санкт-Петербург, Россия
АДАПТИВНОЕ УПРАВЛЕНИЕ ГЕНЕРАТОРОМ ПОСТОЯННОГО ТОКА С НЕЗАВИСИМЫМ ВОЗБУЖДЕНИЕМ, РАБОТАЮЩИМ В ИМПУЛЬСНОМ РЕЖИМЕ
Рассматривается вопрос о возможности повышения качества формирования импульсной последовательности тока на выходе генератора постоянного тока с независимым возбуждением для создания сильного импульсного магнитного поля. Генератор располагается на одном валу с разгонным асинхронным двигателем и приводным синхронным двигателем и входит в состав импульсного источника тока большой мощности. Ввиду наличия нелинейностей в математическом описании объекта и широкого диапазона изменения значений формируемого тока предложен вариант прямой беспоисковой нелинейной адаптивной системы управления с эталонной моделью с известными функциями роста с алгоритмами сигнальной настройки, позволяющей воспроизводить с высокой точностью программно заданную импульсную последовательность.
Ключевые слова: импульсный источник тока большой мощности, генератор постоянного тока с независимым возбуждением, адаптивная система управления.
Авторы заявляют об отсутствии возможных конфликтов интересов.
DOI: 10.24937/2542-2324-2019-2-S-I-198-204 UDC 621.313.12
N. Dobroskok, О. Mokhova
TSNII SET branch of KSRC, St. Petersburg, Russia
ADAPTIVE CONTROL OF SEPARATELY EXCITED DC GENERATOR WORKING IN PULSE MODE
This paper discusses the possibility of improving output current pulse pattern of separately excited DC generator to obtain a strong pulsing magnetic field. The generator has a common shaft with its accelerating synchronous motor and forms a part of heavy-duty pulsing power supply. In view of the non-linearities in the mathematical description of the object, as well as due to wide variation range of obtained current, this paper suggests a variant of direct search-free non-linear adaptive control system with the reference model that will feature known growth functions with signal tuning algorithms, which will provide for high-fidelity reproduction of software-defined pulse pattern.
Keywords: heavy-duty pulsing power supply, separately excited DC generator, adaptive control system. Authors declare lack of the possible conflicts of interests.
Одним из основных требований к качеству импульсных источников питания, в частности импульсных источников тока (ИИТ) большой мощности, является высокая точность воспроизведения импульсного тока в нагрузке, обуславливающего формирование сильного импульсного магнитного поля, которое предна-
значено для электромагнитной обработки судов. Для выполнения данного требования решается ряд задач, связанных с определением структуры ИИТ, способного прецизионно воспроизводить программно заданную импульсную последовательность тока нагрузки, состоящую из разнополярных импульсов
Для цитирования: Доброскок Н.А., Мохова О.В. Адаптивное управление генератором постоянного тока с независимым возбуждением, работающим в импульсном режиме. Труды Крыловского государственного научного центра. 2019; Специальный выпуск 2: 198-204.
For citations'. Dobroskok N. A., Mokhova О. V. Adaptive control of separately excited DC generator working in pulse mode. Transactions of the Krylov State Research Center. 2019; Special Edition 2: 198-204 (in Russian).
трапецеидальной формы с заданными параметрами и законом убывания их амплитуды. При этом выбор структуры должен быть согласован с параметрами первичного источника электроэнергии. В работе рассмотрен случай сети бесконечной мощности. Полагаем, что в случае, когда мощность первичного источника электроэнергии недостаточна для формирования импульсов максимальной амплитуды либо динамические характеристики набора мощности не соответствует параметрам импульсной последовательности, в структуру ПИТ возможно ввести накопитель энергии.
На основании результатов предварительных исследований выбран электромеханический вариант исполнения ПИТ большой мощности, поэтому в данной работе рассматривается трехмашинный агрегат, состоящий из механически связанных разгонного асинхронного двигателя (РАД), приводного синхронного двигателя (СД) и генератора постоянного тока с независимым возбуждением (ГПТ НВ) с двумя обмотками возбуждения (ОВ). Существующие аналоги электромеханических ПИТ большой мощности позволяют создавать импульсную последовательность тока с точностью в несколько процентов от амплитуды текущего импульса, однако для минимизации времени и повышения качества электромагнитной обработки судов требуется точность порядка долей процента от амплитуды формируемого импульса. Стоит отметить, что решение задачи обеспечения такой точности напрямую связано с возможностями средств измерения силы постоянного электрического тока, которые в ограниченном диапазоне условий окружающей среды позволяют получить значение абсолютной статической погрешности, составляющее единицы ампер при измерении тока величиной несколько килоампер. При составлении математического описания исследуемого объекта должны быть учтены особенности его конструкции и обусловленные ими нелинейности различного рода, а также широкий диапазон изменения значения формируемого тока. Поэтому в качестве системы управления предложена прямая беспоисковая нелинейная адаптивная система управления с известными функциями роста с эталонной моделью с алгоритмом сигнальной настройки, которая позволяет обеспечить требуемое качество формирования выходного тока ГПТ НВ.
Работа ПИТ большой мощности, функциональная схема которого представлена на рис. 1, состоит из подготовительного и рабочего этапов.
Подготовительный этап подразумевает под собой разгон ПИТ большой мощности от неподвиж-
ного состояния до синхронной скорости согласно следующему алгоритму: приведение во вращение за счет РАД; подача питания на ОВ СД, работающего в генераторном режиме, при достижении половины от номинальной скорости вращения системы; синхронизация СД, работающего в генераторном режиме, с сетью с использованием синхронизирующего устройства; подача питания на трехфазную обмотку статора СД, работающего в генераторном режиме; переход СД, работающего в генераторном режиме, в двигательный режим; вхождение СД в синхронизм; отключение питания РАД.
В течение рабочего этапа решаются две задачи. Первая - стабилизация скорости вращения системы посредством управления ОВ СД по сигналу обратной связи от датчика положения ротора при питании статорной обмотки от сети бесконечной мощности. Вторая - воспроизведение с высокой точностью программно заданной импульсной последовательности тока нагрузки, состоящей из разнополярных импульсов трапецеидальной формы с заданными значениями длительностей фронта Гф, поддержания вершины Гв, среза Тс импульса, паузы между импульсами в соответствии с экспоненциальным законом убывания амплитуды А импульсов посредством управления ОВ ГПТ НВ по сигналу обратной связи от датчика тока якорной цепи.
В рамках исследования процессов ПИТ большой мощности в импульсном режиме работы допускается, что система на момент начала формирования выходной последовательности тока вращается со стабилизированной синхронной скоростью.
Рис. 1. Функциональная схема ИИТ большой мощности: АВ - автоматический выключатель; ДТ -датчик тока; УС - устройство стабилизации; МСУ -микропроцессорная система управления; Ucд -напряжение питания статора приводного СД; LHf гн - индуктивность и активное сопротивление нагрузки соответственно
Fig. 1. Functional layout of heavy-duty pulsing power supply: AB - automatic switch; ДТ - current gauge; УС - stabilizer; МСУ - microchip-based control system; Ucд - stator feed voltage of synchronous drive; LHr rH - inductance and active resistance of load respectively.
Синтезу системы автоматического управления предшествует составление математического описания ГПТ НВ. Представление экспериментальной характеристики холостого хода (XXX) генератора в виде непрерывной функции может быть произведено путем аппроксимации зависимости электродви-жущсй силы (ЭДС) /-.'я от тока возбуждения /0в полиномом вида
1 = /»
•3 -5 -7 -9
a\lOB + а37ОВ + а57ОВ + allOB + а91ОВ ■
(1)
На основании [1] параметры кривой намагничивания (КН) могут быть рассчитаны через значения ЭДС:
Еа 1
/кн (гов) -
30с,
%B+4pWOBfl)
со
гов
%B+4pWOBfl
_А
j со
(2)
(Xmaх + Хп )
(/<((/- 1)Г + 7ф)) л (/>1); tanh(xx + а0)-й,
(i > (U - W + 7ф)) А (t < и -1 )Т + 0,5Т + 7ф); tanh(xx -а0) + й,
(i > (О' -1)7- + 0,5Т + 7ф)) л (t < и -1 )Т + Т\
(3)
где
v 'ов(0
'ов тах
Т — т
^ _ max ^
Хтах0 Д Хтах Д
•sign(xmax)x
fxxx О'ов (0
F
V ятах
- tanh(xmax0 - sign(xmax )а0)
(4)
Т — Т
sign(*max)x
fxxx Оов (0
F
V ятах
- tanh(xmax - sign(xmax )a0)
где /с,,, - коэффициент трансформации, учитывающий соотношение числа витков основной ОВ м>0в и дополнительной ОВ (ОВД) >совд; Рт ~ число пар полюсов; сэ - конструктивная постоянная ЭДС; со -угловая частота вращения ИИТ большой мощности.
Будем полагать, что при формировании фронта и полки первого из пары разнополярных импульсов ЭДС будет изменяться по основной КН, полученной полиномиальной аппроксимацией экспериментальной XXX (1). Далее вплоть до момента окончания формирования этой пары ЭДС определяется гистере-зисной кривой ГПТ НВ, математическое описание которой получено на основе 7(х)-\юдсли Такакса [2]. Следовательно, процесс изменения Ил для любой пары импульсов можно представить в виде системы уравнений
(О = /г ('ОВ ) = ¿ятах Х
/ххх(хД (;<7ф)лО =1);
(х„ +х„)
/ххх(хтгк'<3в тах)
где /'овтах, ¿ятах - номинальные значения тока воз-буяедения и ЭДС генератора соответственно; х8 -нормированное значение тока возбуждения; хтах0 -значение исходной (предыдущей) вершины петли гистерезиса, которая в зависимости от знака производной тока якоря принимает либо положительное значение амплитуды текущей пары импульсов А (нисходящая кривая петли гистерезиса), либо отрицательное (восходящая кривая петли гистерезиса); Хщях - значение конечной (текущей) вершины петли гистерезиса; а „ /> - коэффициенты, характеризующие коэрцитивную силу и остаточную намагниченность соответственно и принимающие значения от 0 до 1; хп - нормированное значение точки пересечения нисходящей кривой оси абсцисс текущей пары импульсов; ] - номер текущей пары импульсов.
Математическая модель ГПТ НВ с двумя ОВ с учетом явления магнитного гистерезиса (3) может быть представлена в виде системы нелинейных дифференциальных уравнений относительно тока возбуждения /0в и тока якоря /я:
_ _ г и л _
ш
(¡Х2 _ гг л_
,, _С21/ЛХ1-> а22Х2>
Ш
(5)
где
Х1 — 'ов' -*2 — гя' а\ 1 —
ГОВ +ГОВД
(Аэв +4>вд)
Кв+^овд). _ 1 ^11 — > — >
®о (А:)В + Ажд) (А:)В + Азвд)
(6)
«22 =
Гя+Гн-; с21 =]/(4+4).
4+V
При этом функция
Ея W = /г ('ов) = (-^ятахАоВ max ) Х
fl =
ах + 3 в3х2 + 5a5xs + + 9а9х^,
(?<Гф)л(/=1);
■/хХх(Хтах'ов max)
< 1 J
(Хтах + Хп ) +-
max п > Х„
(*тах+*и)^ (7)
(/<((/-1)Г + 7ф))л(/>1);
/и' = sech2 + я0),
(# > ((7 -1 )Т+Гф)) л (/ < (/ - 1)Г + 0,5Г + Гф); /m' = sech2(x,-a0),
(/> § - Х)Т + 0,5 Т + Гф)) л (/<(/- \)Т + Т)
получена путем взятия производной от (3).
Матричная форма записи (5) будет выглядеть следующим образом:
dx
— = А(х, i)x + В(х, mi, dt
(8)
Еяи) = fT(xl) = xlfTl(xl ) =
: .V, X
/ш) -
2 4, 6 , 8
Т (Л'хрСу Т Xj Т ClqX-^ Т CIqX-^ 7
(/<ГФ)А(/ = 1);
-^ятах/ххх(-*тах'0в шах
(/<((/-1)Г + П,))д(7>1);
-vmax + хп )) f
/il(lUHff) -Ефвх
+sign(.vmax)
1
■^lmax
17
max / 2
•^lmax
1 4 —I—o,
3 3
2-og + - —ao 1 +
17
on--On H--On +
15
595
+sigll(„Ymax ) — , ^lmax "
&r\ C(r\ ~
595 ---1
315
(12)
7«n +
2 357 . , ч Щ (17 ^
lmax
17 )
—рЩ +3<73.V[ +5o5.vf + +9 o9.V[ +
1
- x (tanh(.vmax - sign(.imax )o0) -
■^lmax("^max ^ Xn)
~ tanh(-Wo " sigpirW )«o))]
(t > ш-\)т+гф))л(/ < (/-dt+гф).
где
A(\.f) =
all Cll/r2(Xl) c2l/rl(*l)
22.
B(x,/) =
0
(9)
При этом;||;(л-;) и /,';(Х|) выведены за счет вынесения XI из/т (хО (3) и/г(х\) (7) соответственно.
В связи с наличием в/,\(х\) и/т(х\) тригонометрических функций в виде гиперболического тангенса 1ап11(х) и косеканса 5ес1г(х) при
х = хх ± а0 = (Щ/xlmas) -siffli('.vnia4 )я0
(10)
дляи /т2 (х,) в работе допускается использование разложения указанных функций в степенной ряд Тейлора:
tanli(x) = х- —х3 + —х5 --^-х7; w 3 15 315
\ л 2 4 17 л 41 seclr (X) = 1-х" + —х----х +-х"
m
45
320
Если учесть (3), (7) и (11) с (10). тогда ДЩ) и /Т2(х\) могут быть описаны с помощью следующих систем уравнений:
E%V) = /г '(*i) = Jfi /ЙЦ ) =
ЗйГз-V! +5а3х^ +7cijXi +9адх^, (/<Гф)А(/ = 1); 0. (/<((7-1)Г+Гф))л(/>1);
/l(l)'-
/п(1)Л11 (1) ' ~~ -^яшах
sign(.Tmax)
н 3 102 5 н-1 7 |
тах'х1тах^°"зЙГо+'4?ао"40й,° ' +
41
, 17 4 1148 б1
■1 +4о„--о„ +-о„ +
о 3 о 320 о
- / 4 ( 8
+sign(xmax)—- --о0 +
X____^ J
340 , 2296 , , о„ +
45
40
2 17 -—-я'п +"
2870 320
-оп +
х 4 Г102
+sign(.*max)^— —о0 -ЛТ____^ 45
2296 320 '
17 1148 2 45 320
Яп +
■ / ч "Vl6 ( 4i Xl f 41
lmax lmax
(/>((/-1)Г + Гф))л(/<(/-1)Г+Гф).
(13)
При составлении эталонной модели ГПТ НВ допускается, что XXX /ххх^-г1 )• ВВИДУ малости коэффициентов полинома (1) при старших степенях, может быть сведена к равенству
/ххх(*1) = «1*1> (14)
где а\ - коэффициент ./ххх(-г1) при х\ в первой степени.
Матричная форма записи эталонной модели ГПТ НВ на основании (8) и (14) примет вид
а
dt
хэ = Аэхэ +Вэиа(0,
(15)
А =
С 21 Ct^ &'
22
= Вч = С12
' э 0
(16)
lim<2(/) = lim (о,5етРе) = О,
t—Ьсс t—bcc V '
(17)
А, Р + РА, = -G,
(18)
lim
q+l,r\
I f
| Jqt
= q = Q,p, г = 1, и,
(19)
Необходимо отметить, что при нулевой степени роста скалярные функции принимают постоянное значение и равны единице:
|/о9,(*,)|
= 1, г = 1, п.
(20)
После введения понятия функций роста элементы матрицы А(х.[) (9) агу(х,7) могут быть представлены через функции в следующем виде:
Yaqr(X^fqr(Xr)
,i,j=\n, (21)
где матрицы динамики Аэ и входов Вэ выглядят следующим образом:
—а,, 0
где а(11(\.1) - глобально ограниченные функции.
Объединив разложения (21) по всем индексам у, можно получить разложение А(х.[):
Цель управления, как указано в [3], может быть задана через положительно определенную квадратичную форму, характеризующую сближение решений эталонной модели (15) и объекта (8):
ам=Е
д=0
П
_г—1
q = 0,p, г = \п,
где е - вектор рассогласования фактических х и эталонных хэ значений переменных состояния; Р -постоянная положительно определенная матрица, являющаяся единственным решением уравнения Ляпунова
где Aqr (х, 0 = { dJqr (х, /)}, j = 1, п,
при этом А0Дх,/) = А0(х,/), г = \.п.
(22)
(23)
(24)
где С - симметричная положительно определенная матрица.
Таким образом, выполнение предельного соотношения (17) означает обеспечение асимптотической устойчивости тривиального решения е(/)=0 системы (8).
Далее используется понятие функций роста [4], представляющих собой скалярные функции /?ег(х), ранжированные по первому индексу по степеням роста:
В связи с тем, что не все элементы матрицы А(х.[) могут иметь максимальную степень роста, равную р. а также не обязательно содержать в своих разложениях (21) функции роста /?ег(х) со всеми степенями роста, некоторые или многие элементы матриц А(1Г(х.[) (23) могут оказаться равными
нулю; также некоторые или многие матрицы А(/1.(х.[) могут вовсе отсутствовать.
Схема прямого адаптивного управления ГПТ НВ с функциями роста представлена на рис. 2.
В соответствии с (9) и с учетом (24) матрицы (23) для системы (8) будут иметь вид
Aqi - А02 -
где (п/11) - размерность матрицы А(х.[) (9); р -наибольшая степень роста А(х.[) по аргументу х, равная количеству функции по элементу хг
Ац =
Ai2 -
41
"С1 0
0 0" 0 0
-«11 0
0 - «22
0" 0
0 ' A2I =
С21
"0 0"
^-22 = 0 0
(25)
При этом функции роста могут быть заданы следующим образом:
(26)
/0ei(*i) = /02Ы = 1; /n(*i) = /г2<л); fl1 (*i) = /н (*i); /12 й) = /22 ) = о.
Структура адаптивного закона для случая параметрической настройки, согласно [4], определяется уравнением
^а(0 = КАх + Квг/ов, (27)
где т - размерность вектора и (8); КА, Кв - матрицы настраиваемых параметров размерностью (т*п) и (тхт) соответственно.
Вид закона (27) определен условиями инвариантности к эталонной модели, получаемыми при приравнивании правых частей объекта (8) и эталонной модели (15):
(А + ВКА )х + (В + ВКв )иов = Аэх + Вэ^0в- (28)
Требование асимптотической устойчивости системы (8) будет характеризоваться выполнением условий адаптируемости:
гапк{В,Вэ -В} = гапк{В,Аэ - А} = rank В. (29)
Закон адаптивного управления с сигнальной, или релейной, настройкой, в соответствии с [4] может быть представлен алгебраическим уравнением
j (t) = -h sign |^BTPe J,
(30)
где h - постоянный положительный коэффициент усиления, который характеризует амплитуду (ограничение) релейного управления; знаковая функция sign от векторной величины, например f={f\,..-,fn],
вычисляется покомпонентно:
sign/ = [sign sign /J.
(31)
При наличии в разложении матрицы А(х,/) функций роста закон управления (30) приобретет нелинейный характер и преобразуется к системе алгебраических уравнений:
Yjhqr&^fqrM
г=1
lsign[ßTPe]. (32)
"а(0= I
[<7=0
По аналогии с (24) справедливы равенства к0г=к0г=1,п. (33)
Для математической модели генератора постоянного тока (8) выбраны следующие коэффициенты усиления:
//0 =0,214/й1тах; Нп =0,023/а1шах; Й2! =0,039/^^2=0; /^=0.
(34)
Рис. 2. Функциональная схема адаптивного управления ГПТ НВ с эталонной моделью с функциями роста
Fig. 2. Functional layout of adaptive control for separately excited DC generator featuring reference model with growth functions.
E*
P
^almax ^ яЗтах
О
-E
яЗтах ^я2тах
^яЗтах
'OB
Рис. 3. Кривая намагничивания генератора постоянного тока с независимым возбуждением Fig. 3. Magnetization curve of separately excited DC generator
/ /я
1 я1тах
I я2тах I яЗтах
О
I яЗтах 1 я2тах
-I,
[ я1тах
а)
О 20 40 60 80 100 с
^тах ^ Зтах
0
6)
Li -L-. L L_ 1 ■
20
40
60
80
100 t, с
Рис. 4. Ток якоря (а) и модуль ошибки по току якоря (б) генератора постоянного тока с независимым возбуждением при использовании адаптивного закона управления с алгоритмами сигнальной настройки с функциями роста
Fig. 4. Armature current (a) and its error modulus (b) for separately excited DC generator governed by adaptive control law featuring signal tuning algorithms with growth functions.
В качестве доказательства работоспособности исследуемого объекта приведены результаты моделирования системы (8) с прямой беспоисковой адаптивной системой автоматического управления с нелинейным адаптивным законом управления с алгоритмом сигнальной настройки (30) с коэффициентами усиления (34) и с известными функциями роста (26) с эталонной моделью (15) в виде графиков КН. тока якоря и модуля ошибки по току якоря на рис. 3 и 4 соответственно.
Проанализировав результаты моделирования, молено сделать вывод, что прямая беспоисковая адаптивная система управления с функциями роста с эталонной моделью с алгоритмами сигнальной настройки позволяет получить рассогласование заданного и фактического сигналов в установленных пределах (рис. 46) и удовлетворяет требованиям относительно формы и параметров импульса и закона убывания их амплитуды (рис. 4а).
Заключение
Conclusion
В рамках данной работы с целью повышения точности отработки программно заданной импульсной последовательности тока в нагрузке ПИТ большой мощности, предназначенного для электромагнитной обработки судов, синтезирована прямая беспоисковая адаптивная система управления ГИГ НВ с нелинейным адаптивным законом управления с эталонной моделью с известными функциями роста с алгоритмом сигнальной настройки. Результаты проведенных в САПР Matlab/Simulink иссле-дований показали, что указанная система автоматического управления позволяет обеспечить выполнение требований относительно формы, закона убывания амплитуды, значений параметров импульсной последовательности тока, поэтому она может быть использована для высокопрецизионного создания сильного импульсного магнитного поля с учетом особенностей структуры ПИТ большой мощности.
Библиографический список
1. ВолъдекАМ. Электрические машины. СПб.: Энергия,
1978.-832 с.
2. TakácsJ. Mathematics of hysteresis phenomena. Weinhem: WILEY-VCH Verkag GmbH&Co. KGaA, 2003. - 175 p.
3. КимД.П. Теория автоматического управления. Т. 2. Многомерные, нелинейные, оптимальные и адаптивные системы. М: Физматлит, 2004. -464 с.
4. ПутовВ.В., ШелудъкоВМ. Адаптивные и модальные системы управления многомассовыми нелинейными упругими механическими объектами. СПб.: ООО «Тех-номедиа» jИзд-во «Элмор», 2007. -244 с.
References
1. A. Voldek. Electric machines. St. Petersburg, Eilergiya, 1978, 832 pp. (in Russian).
2. Talajes J. Mathematics of hysteresis phenomena. Weinhem: WILEY-VCH Verkag GmbH&Co. KGaA, 2003. - 175 pp.
3. D. Kim. Automatic control theory. Vol. 2: multidimensional, non-linear, optimal and adaptive systems. Moscow, Fizmatlit, 2004, 464 pp. (in Russian).
4. Щ Putov, V. Sheludko. Adaptive and modal control systems for multi-mass non-linear elastic mechanical objects. St. Petersburg, JSC Technomedia / Elmor Publishing House, 2007,244 p. (in Russian).
Сведения об авторах
Доброскок Никита Александрович, к.т.н., начальник сектора - заместитель начальника отдела филиала «ЦНИИ С" )Т» ФГУП «Крыловский государственный научный центр». Адрес: 196128, Россия, Санкт-Петербург, ул. Благодатная, 6. Тед,: +7 906 247-04-26. E-mail: n_dobroskok®!ksrc.ru.
Мохова Ольга Васильевна, инженер 1 категории филиала «ЦНИИ СЭТ» ФГУП «Крыловский государственный научный центр». Адрес: 196128, Россия, Санкт-Петербург, ул. Благодатная, 6. Тел.: +7 921 412-27-85. E-mail: o.kurdyaeva@gmail.com.
About the authors
Dobroskok, Nikita A., Cand. Sci. (Eng), Head of Sector -Deputy Head of Department, SET Branch of KSRC, address: 6, Blagodatnaya si, St. Petersburg, Russia, post code 196128, tel.: +7 906 247-04-26. E-mail: n_dobro skok(o!ksrc .ru.
Uokhova, Olga V., 1st category Engineer, SET Branch of KSRC, address: 6, Blagodatnaya St., St. Petersburg, Russia, post code 196128, tel.: +7 921 412-27-85. E-mail: о .kurdy aeva(®.gmail. com.
Поступила / Received: 20.06.19 Принята в печать t Accepted: .30.08.19 © Доброскок H.A., Мохова О.В., 2019
