CC BY
43
CALCULATION OF THREE-PHASE CIRCUIT WITH ARBITRARY
NUMBER OF LOADS
Berzan V.
Institute of Power Engineering of the Academy of Sciences of Moldova
Abstract. An algorithm for calculating the normal and emergency modes of three-phase circuits with distributed loads is elaborated. We propose to divide the required unknown currents into two groups: load currents are considered as principal, and currents in the other circuit branches are considered as complementary. This approach allows calculating at the first phase the load currents, and then other currents using the first Kirchhoffs law. The form of the matrix equation and rules of matrix coefficients composition for each phase depends on the number of loads connected to the phases of the circuit. Matrix equation includes 3n unknown currents. The proposed method allows calculating both the symmetric and asymmetric modes in the three-phase unbalanced network, using the same calculation algorithm. The fault modes, such as a single-phase, two-phase, three-phase and phase-to-phase short-circuits as well as the open circuits with some phase loads (simulation of the no-load mode in any node and in any phase of the circuit), could be also calculated.
Keywords: Matrix method, three-phase circuit, symmetrical and non-symmetrical distributed load.
CALCULAREA CIRCUITULUI TRIFAZAT CU UN NUMÄR ARBITRAR DE SARCINI
Berzan V.
Institutul de Energeticä al Academiei de §timte a Moldovei Rezumat. În lucrare se propune un algoritm generalizat de calcul a regimurilor stationare normale çi de avarie în circuitele cu trei faze. Se propune separarea marimilor necunoscute în doua grupe: curentii în sarcinile conectate la fazele circuitului care sunt denumite marimi necunoscute principale çi ceilalti curent din ramurile circuitului pe care le denumim grupa marimilor necunoscute secundare. Urmare a acestei separari procedura de calcul a marimilor necunoscute consta din doua etape. La prima etapa se calculeaza curentii în sarcini, utilizând a doua teorema a lui Kirchhoff, iar apoi curentii clasifica^ ca marimi secundare dupa prima teorema a lui Kirchhoff. S-a propus forma matriciala generalizata a sistemului ecuatiilor de calcul çi recomandarile de întocmire a matricelor coeficientilor de pe lânga marimile necunoscute pentru fiecare faza. Ecuatia matriciala generalizata include 3n ecuatii. Metoda propusa permite calcularea regimurilor simetrice çi nesimetrice în circuitele trifazate pentru un numar nelimitat de sarcini conectate la fazele lui. Regimurile de limita cum sunt scurtcircuitele monofazate, bifazate, trifazate çi intre faze, precum çi de mers în gol pentru unele din sarcini se efectueaza în baza unei descrieri unice a circuitului çi algoritmului generalizat de calcul a marimilor necunoscute.
Cuvinte-cheie: Metoda matriciala, circuit trifazat, sarcini distribuite simetrice çi nesimetrice.
РАСЧЕТ ЦЕПИ ТРЕХФАЗНОГО ТОКА ПРИ ПРОИЗВОЛЬНОМ ЧИСЛЕ ПОДКЛЮЧЕННЫХ
НАГРУЗОК Берзан В.П.
Институт энергетики Академии Наук Молдовы Аннотация. В работе предложен алгоритм расчета нормальных и аварийных режимов трехфазных цепей при распределенном характере подключения нагрузок. Предложено разделить искомые неизвестные (токи) на две группы. Токи, которые относятся к нагрузкам считать главными искомыми величинами, а токи в остальных ветвях цепи - дополняющими. Такая постановка задачи позволяет рассчитать на первом этапе токи в нагрузках, а затем и другие токи цепи на основе первого законна Кирхгофа. Предложены записи матричного уравнения и правила составления матриц коэффициентов для каждой фазы в зависимости от числа подключенных нагрузок к фазам цепи. Матричное уравнение включает 3 n неизвестных токов. Предложенный метод позволяет рассчитать режим как в симметричной, так и в несимметричной трехфазной сети на основе одного и того же алгоритма расчета. Предельные режимы такие как: однофазные, двухфазные, трехфазные и междуфазные короткие замыкания, а также, и режимы при отключении некоторых фазных нагрузок (моделирование режима ХХ в любом узле фазы или цепи), также рассчитываются по установленному общему алгоритму. Ключевые слова: Матричный метод, трехфазная цепь, распределенные симметричные и несимметричные нагрузки.
1. Введение
В учебно-методической литературе присутствует утверждение, что при расчете трехфазной цепи применимы все методы, используемые для расчета линейных однофазных цепей[1-4]. При выполнении расчетов принимаются некоторые допущения: сопротивления фазных проводов, нейтрали и внутреннее сопротивление генератора меньше сопротивлений приемников и при упрощенных расчетах режима цепи они не учитываются. Нагрузки могут быть однофазные, трехфазные, симметричные и несимметричные, сосредоточенные и распределенные по длине питающей линии, что имеет в качестве следствия усложнение методики расчета. Также, нагрузки трехфазной цепи могут быть активные, реактивные или иметь комплексный характер, например, иметь активно-индуктивный или активно-емкостной характер.
В схемах замещения эти нагрузки представляются комплексными сопротивлениями, которые могут быть подключены по схеме «звезда» или треугольник. Отметим, что схему соединения нагрузок в «треугольник» можно преобразовать в «звезду» и наоборот [1-3]. Это свойство позволяет привести схемы замещения трехфазных цепей к однотипной структуре, что является несомненным преимуществом при выполнении расчетов режимов таких цепей.
Учет реальных параметров цепи и распределенного характера нагрузок в цепи приводит к определенным сложностям и необходимости применения специальных методов, что усложняет решение поставленной расчетной задачи. Вследствие этого, приходится каждый раз проводить анализ цепи, выбирать метод расчета и формально, каждый раз решать как бы новую задачу, особенно при переменном характере нагрузки. Все это представляется достаточно сложным и трудоемким процессом при расчете трехфазных цепей с произвольной структурой и несимметричными нагрузками.
Поэтому, представляется рациональным и обоснованным рассмотреть задачу разработки подхода, который бы позволил, не изменяя общую структуру расчетного метода обеспечить расчет трехфазной цепи при переменной ее структуре в нормальных и в аварийных режимах.
2. Рабочая гипотеза
При решении задачи применим подход, изложенный в работе [5]. Основная идея этого подхода состоит в разделении токов цепи на две группы: на токи, которые текут в нагрузках и токи, протекающие по проводам фаз и нейтрали. Будем полагать, что токи нагрузок являются главными (базовыми), а оставшиеся токи цепи составляют вспомогательную (второстепенную) группу искомых величин. Нагрузки образуют поперечные ветви (их могут назвать и ветвями связи) в схеме замещения трехфазной четырехпроводной цепи, а участки проводов фаз и нейтрального провода, расположенные между двумя близлежащими точками i и i +1, к которым подключены фазные нагрузки ZAS(i} и ZAS(i+1) образуют продольные
ветви цепи (их можно назвать ветвями дерева). Отличительной особенностью такого разделения цепи и ее параметров состоит в том, что продольные параметры для данной цепи являются постоянными и заданы при проектировании. В то же время, значения параметров поперечных ветвей изменяются в широких пределах, например, от нуля (режим КЗ для данной поперечной ветви) до режима холостого хода, например, при отключении данной нагрузки (поперечной ветви) от цепи. В общем случае к фазам подключены n однофазных и/или трехфазных нагрузок, которые в схеме замещения цепи представляются комплексными сопротивлениями.
Расчет линейных электрических цепей имеет в основу первый и второй законы Кирхгофа. Поэтому, для расчета цепи мы должны иметь число уравнений равных числу искомых неизвестных. В рамках сформулированной задачи рассмотрим схему замещения, представленную на рис. 1, заранее не ограничивая число нагрузок подключенных к цепи.
Рис. 1. Схема замещения трехфазной четырехпроводной цепи с п нагрузками в фазах
Для сохранения системности и однообразности в обозначениях величин, предлагается величины которые относятся к источнику питания обозначить следующим образом:
1 ЕЛ = /А01;1ВА = 1В 01; 1СА = 1С 01;1ЕЫ = /Ы 01; ^ЕА = ^Л0\; ^ЕВ = ^ЕС = ^С 01. В представленной
схеме каждая фаза имеет п нагрузок, а цепь фазы содержит 2п ветвей, в том числе щ поперечных и п продольных ветвей. Однофазные нагрузки, подключенные к той или иной фазе в узле I предлагается дополнить добавлением недостающих поперечных ветвей составленных из импедансов 25чтобы в схеме замещения они приобрели структуру трехфазных нагрузок. Значения импедансов добавленных ветвей должны намного превышать значение самого большого из наблюдаемых комплексных сопротивлений ^тах
цепи, т.е. должно быть выполнено условие 25 } » 25^. Таким способом расчетная схема
преобразуется к виду, чтобы она имела симметричную структуру. Необходимо также преобразовать нагрузки подключенные по схеме «треугольник» в «звезду», т.е Л5 ^ ^ . При выполнении этих условий, число продольных ветвей в нулевом проводе будет равно числу поперечных ветвей в фазах цепи. Выполнение этих условий приводит к тому, что для случая четырехпроводной трехфазной цепи число неизвестных величин равно 7п, т.е. выражается через число нагрузок подключенных к фазе цепи. Количество неизвестных для главной группы токов, т.е. токов ^ нагрузках равно 3п, а 4п токов относятся к
вспомогательной группе токов г, протекающие в продольных ветвях с импедансами г и нулевого провода (г._мгде / = 1,2,3,...,п .
При учете нулевого провода мы можем записать для представленной схемы (т +1) п уравнений по первому закону Кирхгофа, где т - число фаз цепи; п - число узлов
подключения нагрузок к фазе после преобразования и симметрирования структуры цепи. Для того, чтобы замкнуть систему нам требуется составить еще 3п уравнений по второму
закону Кирхгофа. В эту группу уравнений должны войти уравнения, которые включают источник и неизвестные переменные, которые мы включили в главную группу искомых величин-токи в нагрузках.
Нашу задачу мы сформулируем следующим образом. Необходимо предложить обобщенный метод расчета токов в нагрузках трехфазной цепи с нулевым проводом с учетом влияния реальных параметров цепи при подключении к цепи n трехфазных нагрузок.
3. Метод расчета
Для эквивалентной схемы (рис.1) можно записать различные комбинаций уравнений для токов в узлах и падений напряжения в контурах.
В общем виде соотношения для токов в узлах ветвления в каждой фазе и нулевом проводе записываются следующим образом:
IA(i-1,i) = IAS(i) + IA(i,i +1); ^В(i-1,i) =IBS(i) + 1B(i,i+1); 1N(i-1,i) = ^AS(i) + 1BS(i) + 1CS(i), (1)
где i = 1,2,3,..., n.
При включении в цепи n трехфазных нагрузок в четырехпроводной цепи число уравнении для токов будет равно 4n . Для падений напряжения в замкнутых контурах можем записать 3n уравнений. Эти уравнения разделим на 2 группы. Первая группа, которая состоит из 3-х уравнений, представляет сумму падений напряжений в главном контуре фазы и включает фазные ЭДС источника питания. Главный контур, кроме продольных ветвей фаз и нулевого провода источник питания включает только одну поперечную ветвь (ветвь связи),которая также является наиболее удаленной от источника питания. Другая группа включает к подгрупп из трех уравнений падений напряжения в замкнутых контурах включающих две близлежащие поперечные ветви (ветви связи) и двух продольных ветвей (ветви дерева фазы и нулевого провода) фазного провода и нулевого провода соединяющие концы поперечных ветвей с образованием замкнутого локального контура цепи между узлами i и i +1. Через к обозначим порядковый номер рассматриваемого замкнутого контура, причем к = 2,3,4,...., n.
Для главного контура в случае подключения n нагрузок на фазу можем записать следующую подгруппу уравнений для фаз А, В и С:
n n
2 IA(i-1,i)ZA(i-1,i) + 1 ASnZASn +2 1N(i-1,i)ZN(i-1,i) EA ; i=1 i=1
2 IB(i-1,i)ZB(i-1,i) + 1 BSnZBSn +2 1N(i-1,i)ZN(i-1,i) EB ; (2)
i=1 i=1
2 -1,г)^(¿-1,0 + lCSnZCSn + 2 1Ы(г-1,г)^Ы(г-1,г) ЕС 1=1 г=1
где /' =1,2,3,....., п.
В общем виде уравнения для контуров смежных поперечных ветвей можно записать для случая подключения п нагрузок следующим образом:
/А(г,г+1) ^А(г,г+1) + /А8 (г +1) (г+1) + /Ы (г,г+1) ^N (г,г+1) /А8 (г) (г) 0;
/В (г,г+1) ^В (г,г +1) + (г +1) (г +1) + (г,г+1) (г ,г +1) /В8 (г) (г) 0; (3)
/С (г,г+1) ^С (г,г+1) + /С8 (г +1) (г +1) + (г ,г +1) (г ,г +1) - /С8 (г) (г) = 0,
где /' = 1,2,3,....., п -1.
Последовательной подстановкой токов из уравнения (1) в уравнения подгрупп (2) и (3) преобразуем их к виду, чтобы уравнения (2) и (3) содержали в качестве независимых переменных только токи нагрузок / , т.е. токи которые мы определили как главные:
п-1 п п
^ /А8г (^А01 + (г-1)г ) + /А8п (^А0п + ^А8п + (п-1)п ) + ^ /В8г^Ы(г-1)г + ^ /С81^Ы(г-1)г = ЕА ;
г =1 г=1 г
п п-1 п
^ /А81^Ы(г-1)г + ^ (^В0г + (г-1)г ) + /В8п (^В0п + ^В8п + (п-1)п ) + ^ /С81^Ы(г-1)г = ЕВ ; г =1 г =1 г
п п п-1
^ /А81^Ы(г-1)г /В8г^Ы(г-1)г /С5г (^С0г + (г-1)г ) + /С8п (^С0п + ^С5п + (п-1)п ) = ЕС ;
г =1 г =1 г=1
к-2 п
^ -^(г )0 - /А5(к-1)^А5(к-1) + /А8(к)(^А(к-1,к) + ^А5(к ) + (к-1,к )) + ^ (г )^А(к-1,к ) + г =1 г=к+1
к-2 п к-2 п
/В8 (г)0 + /В8 (к) (к-1,к ) + ^ /В8 (г )0 + ^ (г )0 + (к) (к-1,к) + ^ (г )0 = 0; (4)
г=1 г=к+1 г =1 г=к+1
к-2 п к-2
^ (г )0 + /А8(к)(к-1,к ) + ^ /А8(г )0 + ^ (г )0 - (к-1)^В8(к-1) + г =1 г = к +1 г =1
п к-2 п
+/В8(к)(^В5(к-1,к) + (к) + (к-1,к)) + ^ (г)^В (к-1,к) ^(г )0 + (к)(к-1,к) + ^ ^ (г )0 = 0;
г=к+1 г=1 г =к+1
к-2 п к-2 п
^ (г )0 + /А8(к)(к-1,к ) + ^ (г )0 + ^ (г )0 + /В8(к)(к-1,к) + ^ (г )0 +
г =1 г = к +1 г =1 г = к+1
к-2 п
/С8(г)0 - /С8(к-1)^С8(к-1) + /С8(к) (к-1,к) + ^С8(к) + (к-1,к)) + ^ (г)^ВС(к-1,к) = 0,
г=1 г=к+1
где г = 1,2,3,.....,п; к = 2,3,4,....,п.
В зависимости от значения показателя к система уравнений (4), представленная в общем виде, после ее развертывания содержать 3п уравнении с 3п независимыми переменными в каждой строке. Эти независимые переменные можно представить следующим вектором:
Х = (/А51, /А82 ,...., /А8(п-1) ,1 А8п, /В51, /В82 ,...., /В8(п-1), /В5п , /С51, /С52 ,"", /С5(п-1), /С8п ) , (5)
а систему уравнений (4) в матричном виде
ХВ = 7 (6)
где В -квадратная матрица пял коэффициентов при независимыми переменными; 7-вектор известных величин. Вектор 7 имеет п компонент, т.е. 7 = (Ел, Ев, Ес ,0,0,........,0).
Квадратную матрицу В удобно представить в следующем виде, т.е. выделив через отдельные матрицы коэффициентов независимых переменных для каждой фазы:
34
Б]
ВЕ
БЕ
^ _ -^^ASi BSi CSi
BAS(k-1,k) BBS(k-1,k) BCS(k-1,k)
Матрицы BESi, BESi, BESi представляют коэффициенты при неизвестного переменных в уравнениях главного контура (соответственно для фазы А, Б и С) имеют размерность тхп, где т - число фаз цепи, т.е. для трехфазной цепи т =3. Матрицы BAS(к-и), BBS^-и), BCS(к-и)
подразделяются на группы по индексу к = 2,3,4,.....п, т.е. на (п -1) подгрупп размерностью
тхп.
Для главного контура матрицы Bí^S,i, BBSi, BCSi представляются следующими коэффициентами:
^ =
7 + 7
7A01 ^ 7Ы01
N12
N12
7 + 7
7А02 ^ 7Ы12
N12
N12
7 + 7
7А0(1-1) ^ 7N(i-2,i-1)
JN (i■ -2,i-1)
7А0п + 7АЯп + 7N(п-1,п)
JN ( п-1,п)
JN ( п-1,п)
BBSi =
7
7N12
7 + 7
^01 ^ 7N01 7
7
7 + 7
^02 ^ 7N12 7
ЛГ1 о
N (i-2,i-1)
7 + 7
7B0(i-1) ^ 7N(i—2,i—1)
^'N (i-2,i-1)
N ( п-1,п) ZB0n + + ZN(п-1,п)
^JN ( п-1,п)
(8)
^ =
N12
N12
7 + 7
7C01 ^ 7N01
N12
N12
7 + 7
7C 02 ^ 7N12
^ N О^,!-!)
'JN (i—2,i—1)
7 + 7
7C0(i—1) ^ 7N(г—2,i—1)
JN ( п-1,п)
^ ( п-1,п)
0п + ZCSn + (п-1,п)
Для контуров включающие две смежные продольные ветви 7 и 7 структура
матриц зависит от порядкового номера контура к = 2,3,....., п. Первая строка матрицы
-^£(¿-1 к) имеет ненулевые элементы начиная с индекса k — 1. Элементам расположенных
слева от элемента с индексом 7
А£ (k-1)
присваиваются нулевые значения, т.е. 7
AS ^-2)
= 0.
Вторая и третья строка содержат только один ненулевой элемент 7Н^) в столбце i = k,
является комплексным сопротивлением продольной ветви (участка) нулевого провода соединяющие узлы i — 1 и i при i = k. По таким же правилам составляются остальные матрицы для фаз Б и С, которые имеют подобную структуру, что и матрица BAS (к—1Л). Учтем
только следующую особенность, что заполненной строкой с ненулевыми элементами в матрице BBS^-и) являются вторая строка, также начиная с элемента (ь—1). В матрице
BCS(k-1,k) ненулевые элементы третьей строки начинаются с 7^ (к— 1), а элементы 7Н(k_1,k)
расположены также в столбце с номером i = k.
При построении матриц BAS^^),BBS(k-1,k),BCS(k-1,k) коэффициентам при искомых
неизвестных, имеющие индексы , i = k — 2 < 0 присваивают нулевые значения для получения матрицы с размерностью тхп. Матрицы коэффициентов при искомых неизвестных для подгрупп уравнений контуров включающих смежные комплексные сопротивления двух нагрузок 7 и 7 можно с учетом изложенных правил представить в общем виде:
B
AS (k - 1,k)
0 0 0
- Z
AS(k-1) 0
0
(ZA( k-1,k) + ZAS (k) + ZN (k-1,k)
) Z
JN (k - 1,k)
JN (k-1,k)
A(k-1,k) 0
0
'A(k-1,k) 0
0
B,
BS (k-1,k)
0 0 0
-Z
0
BS (k-1) 0
( Z
N (k-1,k) B (k-1,k) + ZBS (k) + ZN (k-1,k)
N(k-1,k)
) Z
0
B( k -1,k) 0
0
7
B(k-1,k) 0
B
CS (k - 1,k)
0 0 0
-Z
0 0
CS(k-1)
N(k-1,k)
JN (k - 1,k)
(ZC
C (k-1,k) + ZCS (k) + ZN (k-1,k)
) Z,
0 0
C (k-1,k)
0 0
7
C(k-1,k)
Выполним проверку соответствия числа уравнений для определения искомых неизвестных, т.е. токов в каждой из п нагрузок. Пусть число нагрузок на фазу равно п = 7. Число искомых неизвестных токов в нагрузках будет т • п = 21 (см. урав. (5)). Количество дополняющих контуров &тах = п -1 = 7 -1 = 6. Поскольку имеем т • £тах = 3 • 6 = 18 б следовательно можем составить 18 уравнений для падений напряжения в локальных контурах и три уравнения системы (2) главных контуров, что дает нам в итоге для 21-ой неизвестных величин и систему из 21-го уравнения, что и обеспечивает возможность определения этих искомых неизвестных.
Искомые значения неизвестных токов 1СЯ при г = 1,2,3,...,п вычисляются на
основе матричного уравнения (6) с применением, например, правила Крамера:
1ASi
В,
Br
Br
B
1 BSi
B
1CSi
B
(10)
где i = 1,2,3,...,n .
Матрицы B7 , B7 B7 получаем заменой столбца из коэффициентов при величинах
^ASi, I est, I est на вектор известных величин Y = (EA, EB, Ec ,0,0,........,0). При известных
значениях комплексных сопротивлений схемы замещения (рис.1) и применение стандартных программ, например MATLAB, легко можно вычислить искомые токи в нагрузках, а зная их вычислить и другие токи цепи, используя соотношения, составленные на основе первого закона Кирхгофа. Если нам известны токи, можно вычислить, также, и другие энергетические характеристики цепи, например, потери в различных участках (ветвях) цепи, активную и реактивную мощность, КПД, cos р итп.
Предложенный подход и алгоритм расчета позволяет вычислять и установившиеся аварийные режимы, например, однофазные, двухфазные и трехфазные короткие замыкания в трехфазной цепи. Преимуществом подхода является то, что при расчете нет необходимости формулировать каждый раз задачу по новому в зависимости от структуры трехфазной цепи и соотношений значений фазных нагрузок.
Выводы
ASi
1. Предложенный подход разделения токов цепи на главные и вспомогательные позволяет алгоритмизировать процедуру расчета режима трехфазной цепи с распределенными нагрузками.
2. Использование матричного метода расчета позволяет формализовать процесс составления уравнений электромагнитного баланса цепи, упорядочить подготовку исходных данных, что особенно существенно при расчете трехфазных цепей со сложной структурой.
3. Изменение структуры осуществляется за счет задания детерминированных значений параметров фазных нагрузок цепи, например, для режима КЗ, конкретной нагрузки (комплексному сопротивлению) присваивают в исходных данных нулевое значение, а в случае отключения от фазы присваивается очень большое значение.
4. Предложен простой алгоритм формирования исходных матриц на основе структуры схемы замещения цепи, который позволяет выполнить и проверку правильности их составления методом сопоставления с аналогичными матрицами для других контуров.
Литература.
[1] Основы теории цепей: Учеб. для вузов /Г.В.Зевеке, П.А.Ионкин, А.В.Нетушил, С.В.Страхов. -5-е изд., перераб. -М.: Энергоатомиздат, 1989. -528с.
[2] Бессонов Л.А. Теоретические основы электротехники: Электрические цепи. Учеб. для студентов электротехнических, энергетических и приборостроительных специальностей вузов. -7-е изд., перераб. и доп. -М.: Высш. шк., 1978. -528с.
[3] Теоретические основы электротехники. Учеб. для вузов. В трех т. Под общ. ред. К.М.Поливанова. Т.1. К.М.Поливанов. Линейные электрические цепи с сосредоточенными постоянными. -М.: Энергия, 1972. -240 с.
[4] http://model.exponenta.ru/electro/0042.htm
[5] Berzan V. Electrofizica §i energetica. Universitatea Academiei de §tiin{e a Moldovei; IE A§M. - Chi§inau, 2014. - 262p. ISBN 978-99-75-62
[6] Berzan V., Tirsu M., Iliescu P. Calculation of Electric Circuit rules of Sources and Distributed Loads. ICHQP 2014. 16th International Conference on Harmonics and Quality of Power (ICHQP). Bucharest, Romania, 25-28 May 2014. SESSION 6C - Renewable Generation/Distributed Generation and Power Quality. IEEE (Journal, Magazine, Conference, Book). ICHQP_187. 978-1-4673-6487-4/14/$31.00 ©2014 IEEE
Сведения об авторе.
Берзан В.П. Доктор-хабилитат технических наук. Область научных интересов: диагностика энергетического оборудования, нестационарные процессы в неоднородных электрических цепях, математическое моделирование, транспорт электроэнергии на большие расстояния, возобновляемые источники энергии. Автор более 220 научных публикаций, в том числе 25 патентов и 12 монографий. E-mail: [email protected]
