CC BY
90
УДК 621.313.001.63
Анализ динамических режимов работы трехфазных трехстержневых трансформаторов в пакете MATLAB
В.А. Мартынов, А.Н. Голубев, А.Е. Евдаков ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»
г. Иваново, Российская Федерация E-mail: eae_rus@mail.ru
Аннотация
Состояние вопроса: В настоящее время расчет динамических режимов различных нелинейных электромагнитных устройств относительно мгновенных токов и напряжений производят численными методами. При этом широкое распространение получил программный продукт MATLAB, являющийся высокоуровневой средой программирования. Однако возможностей встроенных блоков MATLAB недостаточно для адекватного описания цепей, содержащих нелинейные электромагнитные устройства, в частности трехстержневые трансформаторы. Это проявляется при исследовании несимметричных переходных процессов, которые возникают при изменении параметров как самого трансформатора, например при обрыве или коротком замыкании обмоток, так и присоединенной нагрузки. Анализ несимметричных переходных режимов трансформаторов, как правило, проводится с использованием метода наложения принужденного и свободного режимов. При этом расчет принужденного режима производится методом симметричных составляющих, а учет свободных составляющих достигается введением различных коэффициентов, выявленных на основе сопоставления накопленных экспериментальных и расчетных данных. Такой подход базируется на довольно грубых допущениях. В связи с этим актуальной является разработка методов, алгоритмов и программ, позволяющих использовать пакет Simulink для расчета цепей с нелинейными электромагнитными устройствами.
Материалы и методы: Исследование динамических режимов произведено методом переменных состояний с помощью среды MATLAB, в которую интегрирована разработанная программа расчета нелинейной магнитной цепи трехстержневого трансформатора.
Результаты: Разработаны модель и алгоритмы, расширяющие возможности Simulink для расчета цепей с нелинейными электромагнитными устройствами. Произведены расчеты несимметричных режимов работы трехстерж-невого трансформатора с соединением обмоток по схеме Y/Yn.
Выводы: Предложенная модель является универсальной и позволяет анализировать переходные, установившиеся и аварийные режимы при несимметричных нагрузках в трансформаторах с различной схемой соединений. Оценка адекватности произведена на основе верификации результатов численных экспериментов с паспортными данными трансформатора.
Ключевые слова: трехфазный трехстержневой трансформатор, динамические и аварийные режимы, несимметричная нагрузка, программный продукт MATLAB, нелинейное электромагнитное устройство; алгоритмы расчета; схемы соединения обмоток.
MATLAB analysis of dynamic modes of three-phase three-core transformers
V.A. Martynov, A.N. Golubev, A.E. Yevdakov Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: eae_rus@mail.ru
Abstract
Background: Dynamic modes of various nonlinear electromagnetic devices regarding instantaneous currents and voltages are now calculated by numerical methods. That is why the software package MATLAB, a high-level programming environment, has become widely spread. However, its built-in features are not always sufficient to describe circuits containing non-linear electromagnetic devices, in particular three-core transformers. It becomes clear while studying asymmetrical transients which are caused by changing parameters of the transformer itself, such as breakage or short circuit of windings, or by using a connecting load. Asymmetric transients of transformers are normally analysed by the method of superposition of forced and free modes. And the forced mode is calculated by the method of symmetrical components, while the available free components are determined by introducing different coefficients identified by comparing the accumulated experimental and calculated data. This approach is based on rather rough assumptions. The goal of the present research is to develop methods, algorithms and programs which allow us to use Simulink for calculating nonlinear electromagnetic devices.
Materials and methods: The study of dynamical modes in the transformer was performed by the alternating direction method in MATLAB environment that has a program for calculation of nonlinear magnetic circuit of the three-core transformer.
Results: A model and algorithms have been developed extending the capabilities of Simulink for the calculation of circuits with non-linear electromagnetic devices. The asymmetrical modes for the three-core transformer with a Y/Yn winding connection scheme have been calculated.
Conclusions: The proposed model is universal and allows analysing transient, steady-state and emergency modes under asymmetric loads in transformers with different connection diagrams. The evaluation was made by verifying the results of numerical experiments against the passport data of the transformer.
Key words: three-phase three-core transformer, dynamic and emergency modes, assymmetric load, MATLAB software, nonlinear electromagnetic device, calculation algorithms; winding connection schemes.
DOI: 10.17588/2072-2672.2016.4.011-018
В настоящее время расчет динамических режимов сложных цепей, содержащих линейные элементы, не представляет сложности и производится на ЭВМ численными методами с использованием специализированных программных продуктов.
Одним из современных пакетов, реализующих численные методы, является МАНАБ [1], который расшифровывается как «матричная лаборатория». По сути это высокоуровневая среда программирования, содержащая множество функций и встроенных пакетов. Одним из таких пакетов является БтиПпк.
Б1ти!тк имеет множество встроенных блоков для расчета как линейных, так и нелинейных систем. Однако при моделировании цепей, содержащих нелинейные трехфазные трансформаторы [2-4], возникает ряд сложностей, связанных с моделированием электромагнитных полей этих устройств в динамических режимах, особенно несимметричных. Поэтому важной задачей является разработка методов, алгоритмов и программ, легко интегрируемых в Б1ти!тк и позволяющих использовать этот пакет для расчета цепей с нелинейными электромагнитными устройствами.
Ниже предлагается модель, позволяющая использовать Б1ти!тк для расчета цепей с трехфазными трехстержневыми трансформаторами, и программа, позволяющая анализировать динамические и установившиеся режимы работы таких цепей как в симметричных, так и несимметричных режимах.
Активная зона исследуемого трансформатора в разрезе представлена на рис. 1. Вид блока, представляющий этот трансформатор в Б1ти!тк, и схема замещения электрической цепи трехфазного трансформатора с соединением обмоток по схеме У/У0 в этом блоке показаны на рис. 2, где обмотки трансформатора представлены последовательным соединением линейных резисторов, сопротивления которых равны активным сопротивлениям соответствующих обмоток трансформатора, и нелинейными индуктивными катушками, напряжения на зажимах которых [иЕ] обусловлены ЭДС само- и взаимоиндукции.
Для работы блока в Б1ти!тк необходимо на каждом шаге численного интегрирования уравнений переменных состояния определять
уравнения связи между мгновенными напряжениями [иЕ ] и производными токов ветвей
трансформатора [5], что можно сделать на основе анализа его магнитной цепи.
Стержм, Одщткп НН Одгтта ВН
Рис. 1. Активная зона трехфазного трехстержневого трансформатора
Как известно [5], магнитная цепь трех-стержневого трансформатора при любых токах и напряжениях его ветвей может быть представлена аналоговой схемой замещения с сосредоточенными параметрами, которая будет содержать нелинейные и линейные магнитные проводимости и источники МДС.
а) б)
Рис. 2. Разработанный в в^иПпк блок и его схема замещения: а - принципиальный вид блока; б - реальный вид блока в Б1ти!1пк; силовые входы 1, 2 предназначены для подключения к цепи; 3 - измерительные выходы
На рис. 3 представлена магнитная цепь исследуемого трансформатора, на рис. 4 - аналоговая схема замещения его магнитной цепи.
На рис. 3 не показаны цепи магнитных потоков рассеяния, так как они практически не влияют на распределение магнитных потоков в стальном сердечнике. Токи в первичных и вторичных обмотках создают магнитодвижущие силы, представленные источниками МДС: Р1 - для фазы А; Р2, Р3 - для фаз В и С соответственно. МДС фаз А, В и С создают потоки в стержнях сердечника, а при наличии потоков нулевой последовательности могут также создавать поток между ярмами трансформатора Ф4(Ф^ [6]. Стержни и ярма магнитопровода трансформатора представляются в схеме замещения нелинейными элементами, а ветвь, соответствующая потоку между ярмами, представлена в виде линейной магнитной проводимости, которая определяется методами теории поля [7].
Рис. 3. Магнитная цепь трехфазного трехстержневого трансформатора
Рис. 4. Схема замещения магнитной цепи трехфазного трех-стержневого трансформатора
Как показано в [5], каждую ветвь схемы замещения магнитной цепи трансформатора в общем случае можно представить в виде обобщенной ветви, образованной последовательным соединением нелинейного элемента и источника МДС (рис. 4). Вебер-амперные характеристики нелинейных элементов обобщенных ветвей могут быть определены предварительно [7] по известным маркам стали и размерам стержней и ярм магнитопровода трансформатора. При кусочно-линейной аппроксимации характеристик этих нелинейных элементов магнитный поток фк к-й ветви магнитной цепи может быть определен через магнитное напряжение ик этой ветви по формуле
фк =лкик +КкРк + ик, (1)
где лк и Jk - коэффициенты кусочно-линейной аппроксимации вебер-амперных характеристик нелинейных элементов, представляющих соответственно дифференциальную проводимость к-й ветви и управляемый напряжением источник магнитного потока; ик - магнитное напряжение обобщенной ветви.
Пронумеруем слева направо все отрезки аппроксимированной к-й вебер-амперной характеристики. Обозначим через стк один из этих отрезков. Тогда значения Лк и Л в (1) будут определяться номером отрезка стк, который, в свою очередь, зависит от значений фк.
Одним из достоинств кусочно-линейной аппроксимации является то, что если номера отрезков стк всех нелинейных ветвей магнитной цепи трансформатора будут известны, то при изменении значений фк в пределах отрезка стк магнитную цепь можно считать линейной. Естественно, что значения стк, лк и ик каждой ветви будут зависеть от токов ветвей трансформатора [/в ] и должны
определяться на каждом шаге численного интегрирования уравнений переменных состояния на основе анализа магнитной цепи.
Введем следующие обозначения: [Ф] -матрица-столбец магнитных потоков всех ветвей; [и ] - матрица-столбец магнитных напряжений на обобщенных ветвях; [Г] - матрица-
столбец источников МДС всех ветвей.
С учетом принятых обозначений, уравнение (1) можно записать для всех ветвей магнитной цепи в матричной форме:
[фНл][и]+[л][гн4 (2)
где [л] - диагональная матрица дифференциальных проводимостей ветвей; [Л] - матрица-
столбец управляемых напряжением источников магнитного потока.
Как показано в [5], величины МДС ветвей схемы замещения магнитной цепи ЭМУ [Г] могут быть выражены через значения токов ветвей ЭМУ [/в] с помощью сингулярной матрицы преобразования [С]:
[г ]=№ ]. (3)
Матрица [С] имеет число столбцов, равное числу ветвей электрической цепи трансформатора, и число строк, равное числу ветвей п схемы замещения его магнитной цепи. Элементом этой матрицы, расположенным в к-й строке и 7-м столбце, является число wk, равное количеству витков с током 7-й ветви электрической цепи /Ву, охватывающих к-ю ветвь схемы замещения
магнитной цепи. При этом если число витков с током /Ву создают в к-й ветви МДС, совпадающие с выбранным положительным направлени-
ем магнитного потока, то берется со знаком «плюс». В противном случае число wk берется со знаком «минус». Так, для схемы замещения магнитной цепи трансформатора, приведенной на рис. 3, матрицы [Ф], [С ] и [ ¡В ] будут иметь следующую структуру:
[Ф] =
"Ф А" 0 0 0 0]
Фв Фс ; [с ] = 0 0 мв 0 0 Мс 0 0 -ШЬ 0 -ШС ; [и ] =
_ФМ _ 0 0 0 0 0 0 с [и]
[¡в ] = ['
А ¡В 'С
и
1 -
С учетом (3) уравнение (2) можно записать в следующем виде:
[ф] = [л]["Мл][с][/вМ4 (4)
Потокосцепления всех фаз трансформатора определятся выражением [5]
Ы=Сг [Фма/в]. (5)
где элементами матрицы [¿а ] являются линейные само- и взаимоиндуктивности рассеяние катушек трансформатора, не входящие в схему замещения магнитной цепи.
С учетом (2)-(4) выражение (5) можно записать в виде
[*в ] = [С]Г [Л][и ] + [с ]т [л][с ][/в ] + [с]Т [Л ] + + М'в ] = ^м ]['в ] + [С ]т ([л][и ] + [Л ]),
где [¿м ] = [С]Г [л][С] + [^] - (7)
Напряжения ветвей трансформатора, обусловленные ЭДС само- и взаимоиндукции как главного поля, так и полей рассеяния, могут быть выражены через производные потокосцеп-лений ветвей [5]:
(6)
.Ш^г.
[ив ] = ^ = ]
[С] [Л]-
(8)
dt L тл сН 1 ' 1 ' С Выразим магнитные напряжения [и] на обобщенных ветвях схемы замещения магнитной цепи трансформатора через токи [/в] электрической цепи. Для этого составим уравнения по первому закону Кирхгофа для схемы замещения магнитной цепи:
[а][ф] = 0, (9)
где [ А] - редуцированная матрица инциденций [6], которая в нашем случае имеет одну строку:
[А] = [-1 -1 -1 1].
Схема замещения магнитной цепи трансформатора (рис. 4) имеет только один узел, магнитный потенциал которого ф отличен от нуля. Напряжения ветвей [и] связаны с этим потенциалом соотношением [6] [и] = [ А]Т ф. (10)
Подставив в (9) выражения (4) и (10), получим
[ А][Л][ А]Т ф + [ А][Л][С][ ¡в ] + [ А][Л ] = 0. (11) Считая известными значения матриц [л] и [Л], введем обозначение
[Ли ] = [А ][Л][ А ]Т. (12)
С учетом (10), (11) и (12) получим:
^" = -[ А]Т [ли ]-1[А][л][С ] ^. Обозначим
(13)
(14)
[¿-В ] = М-[С]Г [Л][А]Г Ы-1[А][Л][С]. (15) Тогда напряжения ветвей трансформатора, обусловленные ЭДС само- и взаимоиндукции [иЕ] (8), можно выразить через производные токов ветвей:
,с [¡в ]
[ие ] = [^ ]-
сН
(16)
Уравнение (16) устанавливает связь между мгновенными напряжениями [и ] и производными токов ветвей трансформатора, что необходимо для использования данной математической модели в БтиПпк [1]. Однако элементы матрицы [¿В], имеющие размерности индуктив-
ностей, в соответствии с (4), (7), (12) и (15) зависят от значений токов ветвей обмоток трансформатора [/в] и должны определяться на каждом шаге численного интегрирования.
Рис. 5. Алгоритм работы блока
В [5] показано, что при известных токах в обмотках трансформатора [ ] элементы матрицы [¿В ] легко могут быть определены через итерационный расчет магнитной цепи кусочно-линейным методом Ньютона-Рафсона.
На основе изложенного алгоритма была разработана программа, укрупненная блок-схема которой представлена на рис. 5. Для нелинейных ветвей схемы замещения магнитной цепи трансформатора была написана отдельная функция, которая вычисляет вебер-амперные характеристики ветвей и производит их кусочно-линейную аппроксимацию с заданной относительной погрешностью.
В качестве примера по разработанной программе были проведены расчеты воздушного трехфазного трехстержневого трансформатора, паспортные данные которого приведены в таблице.
Паспортные данные трансформатора
Sном, кВА Схема обмоток Напряжение обмоток, кВ Потери кВт /х%
ВН НН Рх Рк
250 У/У0 10 0,4 0,92 3,5 5,5 2,2
Две рассчитанные вебер-амперные характеристики для трех нелинейных ветвей схемы замещения его магнитной цепи приведены на рис. 6.
Рис. 6. Вебер-амперные характеристики нелинейных ветвей магнитной цепи и их кусочно-линейная аппроксимация
На рис. 7 приведена расчетная схема электрической цепи с исследуемым трансформатором. К первичным зажимам подключена система симметричных синусоидальных номинальных линейных напряжений, которая моделируется двумя идеальными источниками ЭДС. В каждую фазу вторичной обмотки подключена резистивная однофазная нагрузка. Расчетная схема в 81ти1тк представлена на рис. 8
Рис. 7. Схема замещения электрической цепи
Рис. 8. Собранная для расчета модель в Simulink
Проверка модели была произведена на основе паспортных данных. Первым опытом был пуск трансформатора на холостом ходу (ХХ). На рис. 9 представлены пусковые токи, здесь мы можем наблюдать броски тока, которые значительно превосходят номинальное значение. Мгновенное значение номинального тока ¡н = 20,4 А, а максимальное значение пускового тока ¡тах = 125 А. Токи трансформатора в принужденном режиме представлены на рис. 10. Действующие значения расчетных токов незначительно отличаются от паспортного.Ток холостого хода по паспорту /ХХ = 2,2 % /н, расчетные токи 1А = 2,8 % /н, 1В= 2,1 % /н, /с= 2,8 % /н.
Рис. 9. Пусковые токи: 1 - /д; 2 - /в; 3 - /с
¡.А 0.0
0,2
/ ! / ; ! ! !
: 1 Л : * / Л /
: / ...1...../ \ ! /......./ \ : / А
: ! 1 ! 1 - /
\ ! / / ■ * 1 : V В
Л" Г С
...../¡д.. ШЧ.
/ : 1 /' '
V ; : 1 : ! ^ 1 " ! Чт- т.;.. ... т / ! ' у .....7.; о град
:
\ ^ и г : Г
! 1: / ! 1: /
; ; ;
Рис. 10. Токи холостого хода: 1 - ¡С; 2 - ¡А; 3 - ¡В
Вторым опытом был опыт короткого замыкания (КЗ). На первичые обмотки подавалось напряжение ик = 5,5 % ин. При данном напряжении в первичных и вторичных обмотках протекали токи, близкие к номинальным. По паспорту, /1кз = 14,4 А, расчетные токи - /А = 14,2 А, /В = 14,2 А, /С = 14,2 А.
На вторичных обмотках по паспорту /2КЗ = 361 А, соответственно, расчетные токи - /а = 360 А, /ь = 360 А, 1С = 361 А (рис. 11).
а) б)
Рис. 11. Токи КЗ: а - на первичных обмотках: 1 - ¡А; 2 - /В; 3 -
¡с, б - 1- ¡а; 2 - ¡ь; 3 - ¡с
При симметричной номинальной резистивной нагрузке на приемнике наблюдается симметрия напряжений (рис. 12), а, например, при обрыве провода фазы С эта симметрия нарушается, что вызвано протеканием токов нулевой последовательности (рис. 13, где иа = 0,96ин; 1УЬ=1,031УН; ис = ин).
Рис. 12 Мгновенные значения фазных напряжений приемника при симметричной номинальной нагрузке (соединение У/Уп): 1 - иа; 2 - иь; 3 - и
Рис. 13. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке - обрыве фазы С (соединение У/Уп): 1 - иа; 2 - иь; 3 - ис
Ряд следующих опытов демонстрирует, как изменяются фазные напряжения приемника в зависимости от магнитной проводимости между ярмами трансформатора при несимметричной нагрузке. Во всех опытах приняты сопротивления: Яа = 1,5 Ян; Яь = 0,5 Ян; Яс = Ян.
Опыт 1. Проводимость между ярмами номинальная (рис. 14, где иа = 0,97ин, иь = 0,96ин, ис = 1,02ин).
Рис. 14. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке Л0 = Л (соединение У/Уп): 1 - иа; 2 - иь; 3 - ис
Опыт 2. Проводимость между ярмами в два раза больше номинальной (рис. 15, где иа = 0,96ин; 1/ь = 0,931/н; ис=1,07ин).
Рис. 15. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке Л0 = 2Л (соединение У/Уп): 1 - иа ; 2 - иь ; 3 - ис
Опыт 3. Проводимость между ярмами в четыре раза больше номинальной (рис. 16, где иа=0,99ин] иь = 0,871/н; ис=1,Ш1н).
Рис. 16. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке Л0 = 4Л (соединение У/Уп): 1 - иа; 2 - иь; 3 - ис
Опыт 4. Проводимость между ярмами в восемь раз больше номинальной (проводимость по воздуху практически равна проводимостям стальных стержней в области насыщения -близко к групповому трансформатору) (рис. 17, где иа = 1,041/н; иь = 0,81 (Ун; ис = 1,17ин).
Рис. 17. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке Л0 = 8Л (соединение У/Уп): 1 - иа; 2 - иь; 3 - ис
При соединении трансформатора по схеме А/Уп такой проблемы не наблюдается, о чем говорит следующий опыт (рис. 18, где иа = 0,99ин; иь = 0,97ин; ис = 0,99ин). Поэтому такое соединение на практике предпочтительнее (количество витков в первичных обмотках увеличено в V3 раз, сопротивления проводов первичных обмоток увеличены в 3 раза, сопротивления полей рассеивания первичных обмоток также увеличе-нывЗ раза).
Рис. 18. Мгновенные значения фазных напряжений приемника при увеличенной магнитной проводимости между торцами ярма (соединение А/Уп): 1 - иа; 2 - иь; 3 - ис
Данная модель может быть подключена и к более сложным цепям (рис. 19) (расчетная модель представлена на рис. 20). Здесь трансформатор подключен к внешней цепи через линию электропередач, а нагрузка носит индуктивный характер. Для повышения КПД включена батарея конденсаторов (длина линии 20 км, Ял = 14,8 Ом, Х1л = 6,2 Ом, Я = 0,7 Ом, 1н = 1,91 10-3 Гн, С = 1,810-3 Ф).
R2A L2A
R2B ¿2B СаВГ
R2C ¿2C
Rha ¿HA
^¡13 LHB
Рис. 19. Схема замещения электрической цепи
Рис. 20. Собранная для расчета модель в Simulink
В качестве несимметричного режима была рассмотрена нагрузка 1а = 1,51н, 1Ь = 0,52н, 1с = 1н (рис. 21, где иа = 1,26ин, иь = 0,68ин, ис = 1,13ин). При отсутствии нейтрального провода из-за смещения нейтрали фазные напряжения могут сильно различаться. Поэтому данный режим работы сети считается очень опасным.
Рис. 21. Мгновенные значения фазных напряжений приемника при несимметричной нагрузке (соединение А/У): 1 - Ua; 2 - Ub; 3 - Uc
Модель позволяет рассчитывать и переходные процессы в различных цепях. Для примера была рассмотрена коммутация в фазе А. Соответствующая схема представлена на рис. 22.
Рис. 22. Схема коммутации
В опыте коммутация происходит в момент времени t = 0,04 с. Сопротивление фазы А изменяется с I = 0,92 Ом (Я-1-нагрузка), до I = 0,01 Ом (Я-нагрузка). Расчетные значения мгновенных значений токов и напряжений приведены на рис. 23-25.
Рис .23. Мгновенные значения фазных напряжений приемника при изменении сопротивления фазы А (соединение А/У): 1 - Uc; 2 - Ua; 3 - Ub. Коммутация в t = 0,04 c
с
CA
C
Рис. 24. Мгновенные значения фазных токов на вторичных обмотках при изменении сопротивления фазы А (соединение
А/У): 1 - /0; 2 - /а; 3 - /ь. Коммутация в t = 0,04 с
¡.А
J1
if
L 2 3 1 Ii 2 1 з/
/ \ г
д \/
'V >2 U и iL 15 >6 , 17
V; \
t.c
-40
Рис. 25. Мгновенные значения фазных токов в лини электропередач при изменении сопротивления фазы А (соединение А/У): 1 - А; 2 - /В; 3 - /С. Коммутация в t = 0,04 c
Список литературы
1. Черных И.В. Моделирование электротехнических устройств в MATLAB, SimPowerSystems и Simulink. - М.: ДМК Пресс; СПб.: Питер, 2008. - 288 с.
2. Алтунин Б.Ю., Кралин А.А., Гуляев В.В. Исследование режимов работы нелинейного трехфазного трансформатора в пакете Simulink // Вестник волжской государственной академии водного транспорта. - 2012. - № 32. -С. 195-198.
3. Бошняга В.А. Моделирование трехфазных трансформаторов для расчета несимметричных режимов работы электрических сетей // Известия вузов и энергетических объединений СНГ. Сер. Энергетика. - 1997. - №1-2. - С. 14-19.
4. Закарюкин В.П., Крюков А.В. Сложно несимметричные режимы электрических систем. - Иркутск: Изд-во Иркутского гос. ун-та, 2005. - 274 с.
5. Мартынов В.А., Голубев А.Н., Алейников А.В. Применение метода переменных состояния к анализу стационарных и динамических режимов нелинейных электромагнитных устройств // Вестник ИГЭУ. - 2015. - Вып. 4. -С. 27-32. doi: 10.17588/2072-2672.2015.4.027-032.
6. Основы теории цепей: учеб. для вузов / Г.В. Зеве-ке, П.А. Ионкин, А.В. Нетушил, С.В. Страхов. - 5-е изд., пе-рераб. - М.: Энергоатомиздат, 1989. - 528 с.
7. Иванов-Смоленский А.В. Электрические машины: учеб. для вузов в 2 т. Т. 1. - 2-е изд., перераб. и доп. - М.: Изд-во МЭИ, 2004. - 656 с.
References
1. Chernykh, I.V. Modelirovanie elektrotekhnicheskikh ustroystv vMATLAB, SimPowerSystems i Simulink [Simulation of electrical devices in MATLAB, SimPowerSystems and Simulink]. Moscow, DMK Press; Saint-Petersburg, Peter, 2008. 288 p.
2. Altunin, B.Yu., Kralin, A.A., Gulyaev, V.V. Issledovanie rezhimov raboty nelineynogo trekhfaznogo transformatora v pakete Simulink [A study of nonlinear three-phase transformer modes in Simulink]. Vestnik volzhskoy gosudarstvennoy akade-mii vodnogo transporta, 2012, no. 32, pp.195-198.
3. Boshnyaga, V.A. Modelirovanie trekhfaznykh trans-formatorov dlya rascheta nesimmetrichnykh rezhimov raboty elektricheskikh setey [Modeling of three-phase transformers for the calculation of asymmetrical modes of operation of electrical networks]. Izvestiya vuzov i energeticheskikh ob"edineniy SNG. Energetika, 1997, no. 1-2, pp. 14-19.
4. Zakaryukin, V.P., Kryukov, A.V. Slozhno nesimme-trichnye rezhimy elektricheskikh system [Complicated asymmetric modes of electrical systems]. Irkutsk, Izdatel'stvo Irkutskogo gosudarstvennogo universiteta, 2005. 274 p.
5. Martynov, V.A., Golubev, A.N., Aleinikov, A.V. Prime-nenie metoda peremennykh sostoyaniya k analizu statsionarnykh i dinamicheskikh rezhimov nelineynykh elektromagnitnykh us-troystv [Application of the state variables method to the analysis of stationary and dynamic regimes of nonlinear electromagnetic devices]. Vestnik IGEU, 2015, issue 4, pp. 27-32. doi: 10.17588/2072-2672.2015.4.027-032.
6. Zeveke, G.V., Yonkin, P.A., Netushil, A.V., Strakhov, S.V. Osnovy teorii tsepey [Fundamentals of circuit theory]. Moscow, Energoatomizdat, 1989. 528 p.
7. Ivanov-Smolensky, A.V. Elektricheskie mashiny: v 2 t, t. 1 [Electrical machines: in 2 vol., vol. 1]. Moscow, Izdatel'stvo MEI, 2004. 656 p.
Мартынов Владимир Александрович,
ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор, зав. кафедрой теоретических основ электротехники и электротехнологий телефон (4932) 26-99-08, e-mail: zav@toe.spu.ru
Голубев Александр Николаевич,
ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,
доктор технических наук, профессор, профессор кафедры теоретических основ электротехники и электротехнологий телефон (4932) 26-99-03, e-mail: alenikgo@yandex.ru
Евдаков Алексей Евгеньевич,
ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», студент электроэнергетического факультета, e-mail: eae_rus@mail.ru
