УДК 621.3.072:519.673(075)
Забудский Е.И.
Московский государственный агроинженерный университет
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЭЛЕКТРОМАГНИТНЫХ
РЕЖИМОВ СОВМЕЩЕННЫХ УПРАВЛЯЕМЫХ УСТРОЙСТВ
Введение
К совмещенным управляемым устройствам относятся насыщающиеся и управляемые реакторы (НР и УР). Являясь средством автоматического регулирования реактивной мощности, они необходимы для управления режимами электроэнергетических систем с целью решения следующих задач: компенсация избыточной зарядной мощности ЛЭП и повышение их пропускной способности, ограничение коммутационных перенапряжений, ограничение токов короткого замыкания, уменьшение колебаний напряжения, рациональное распределение напряжения и тока и др. [1, 2, 3].
Реактор это статическое силовое нелинейное устройство, работа которого основана на явлении электромагнитной индукции. Реактивная мощность, потребляемая реактором, плавно регулируется за счет изменения насыщения (магнитного сопротивления) его магнитопровода.
Активная часть реактора содержит одну или несколько обмоток и магнитопровод, выполненный из электротехнической стали. В реакторе могут быть следующие отдельные обмотки: 1) рабочая обмотка (РО). Она предназначена для включения в электрическую цепь, в которой используется индуктивность реактора; 2) обмотка управления (ОУ). Она предназначена для создания управляющего магнитного поля. В большинстве случаев это постоянное магнитное поле или поле, имеющее значительную постоянную составляющую; 3) компенсационная обмотка; 4) фазосдвигающая обмотка. Функции отдельных обмоток может выполнять одна совмещенная обмотка. Магнитопровод изготовляется или стержневым, подобно магнитопроводу силового трансформатора (реактор трансформаторного типа) или кольцевым, по типу магнитопровода электрической машины переменного тока с неявновыраженными полюсами (реактор электромашинного типа), но с неподвижным ротором.
Достоинство реакторов трансформаторного типа состоит и в том, что они выполняются на трансформаторные напряжения и поэтому могут быть включены в энергосистему без промежуточного трансформатора. Реакторы электромашинного типа отличают такие положительные качества, как компактность и простота конструкции при многофазном исполнении, отсутствие взаимоиндуктивных связей между обмотками и, следовательно, повышенное быстродействие, а также симметричность и синусоидальность рабочего тока.
Насыщающиеся и управляемые реакторы предназначены для использования в воздушных линиях электропередачи традиционного конструктивного исполнения и повышенной натуральной мощности [2], в распределительных электросетях и в системах электроснабжения промышленных предприятий, поэтому важное значение имеет анализ электромагнитных режимов их работы в этих энергосистемах.
Исследованию режимов, характеристик и процессов устройств должна предшествовать реализация триады “модель - алгоритм - программа”. В результате осуществляется замена устройства его моделью, которая затем анализируется посредством экспериментирования на ПЭВМ при помощи вычислительно-логических алгоритмов.
Математическая модель глубже вскрывает внутренние связи устройства, дает его точные количественные характеристики. Вычислительный эксперимент частично или полностью заменяет натурное экспериментирование, позволяя уменьшить сроки и стоимость разработок. Универсальность математических моделей, алгоритмов и программ дает возможность оперативно и без дополнительных затрат переходить от решения одной проблемы к другой.
На основе разработанной обобщенной математической модели электромагнитных режимов и явлений, возникающих при намагничивании магнитопровода, исследуется поведение (“физиология”) в электроэнергетических системах двух типичных реакторов: НР трансформаторного типа и УР электромашинного типа.
Описание изготовленных и установленных на эксплуатацию устройств и разработанных программных продуктов, возможность свободного копирования некоторых из них и др. реализованы на Web-странице автора http://zabudsky.ru
1. Насыщающийся реактор трансформаторного типа
Разработан трехфазный НР активная часть которого состоит из стержневого магнитопровода и одной совмещенной обмотки (рис. 1). Реактор предназначен для работы в схеме параметрического стабилизатора напряжения с целью ограничения колебаний напряжения в системе электроснабжения промышленных предприятий, которые обусловлены нагрузкой. Он также может быть использован для компенсации избыточной реактивной мощности энергосистемы в схеме статического компенсатора реактивной мощности. Устройство и теория НР рассмотрены в [3, 4].
1.1. Математическая модель электромагнитных режимов НР
Математическая модель составляется на основе эквивалентирования реального устройства электрической (рис. 2) и магнитной (рис. 3) схемами замещения с нелинейными сосредоточенными параметрами [4].
На рис. 2 приняты следующие обозначения: 1) uAB, uBC, uCA - мгновенные значения линейных напряжений, подведенных к совмещенной трехфазной обмотке; 2) i1A, i1B, i1C - мгновенные значения линейных токов (искомые величины); 3) RA, RB, RC - активные сопротивления фаз обмотки; 4) LsA, LsB, LsC - индуктивности, обусловленные магнитными потоками рассеивания фаз обмотки; 5) RAМ, RBМ, RCМ - нелинейные сопротивления, их значения приближенно учитывают потери в стали магнитопровода; 6) p = d/dt - символ дифференцирования по времени t; 7) pYA = - eA, pYB, pYC - мгновенные значения ЭДС индуктированных в фазах обмотки основными магнитными потокосцеплениями YA, YB, YC этих фаз.
На рис. 3 приняты обозначения: 1) i1i WA, i1iWB, i1iWC, i = A, B или C - мгновенные значения МДС катушек фаз совмещенной трехфазной обмотки; 2) Ф1, Ф2
,..., Ф9 - магнитные потоки в стержнях (искомые величины); 3) RC1, RC2,..., RC9 - нелинейные дифференциальные магнитные сопротивления стержней. Эти сопротивления определяются на основе аппроксимации кривой намагничивания электротехнической стали магнитопровода кубическим сплайном. Экспериментальная кривая намагничивания сглаживалась методом наименьших квадратов; 4) Rd1, Rd2,..., Rd9 - линейные магнитные сопротивления эквивалентных зазоров.
После ряда преобразований [4], обобщенная математическая модель представляющая собой систему нелинейных обыкновенных дифференциальных уравнений (ОДУ), записывается в виде системы (1).
Разработаны алгоритм и программная реализация модели на ПЭВМ [5]. Результатом последней являются искомые значения: токов обмоток, магнитных потоков, индукций и др. величин, получаемых как функции времени.
Система нелинейных ОДУ (1), относительно производных от потоков pФi, i = 1, 2,.., 9 и от токов pi1A, pi1B, pi1C, является системой линейных алгебраических уравнений (ЛАУ). Перед ее решением необходимо задать все величины, входящие в матрицу коэффициентов и в вектор-столбец свободных членов системы ЛАУ, а нелинейные дифференциальные магнитные сопротивления, и нелинейные сопротивления, приближенно учитывающие потери в стали, рассчитать. После решения системы ЛАУ и нахождения значений производных от потоков и токов, они интегрируются численным методом и определяются искомые мгновенные значения потоков и токов для некоторого момента времени t лежащего на интервале интегрирования. Результаты интегрирования используются при формировании матрицы коэффициентов и вектора-столбца свободных членов системы (1) с целью ее решения для последующего момента времени, который определяется переменным шагом интегрирования. Значение шага интегрирования вычисляется программно для достижения заданной точности. Процесс продолжается до тех пор, пока не будет пройден весь интервал интегрирования.
Математическая модель электромагнитных режимов и явлений, возникающих при намагничивании других статических устройств различного функционального назначения, активная часть которых содержит магнитопровод с обмотками, будет структурно совпадать с (1). Это следует и из сопоставления (1) и (3). Обобщенная математическая модель в структурированной (блочной) матричной форме записывается следующим образом:
Для различных устройств (2) будет отличаться лишь числом уравнений и порядком заполнения выделенных пунктирными линиями блоков в матрице коэффициентов и в векторах-столбцах, которые определяются, в основном, конструктивными и схемотехническими особенностями устройств. В связи с этим, программная реализация математической модели различных устройств с точки зрения алгоритма не претерпевает принципиальных изменений.
Из анализа математической модели (1) следует, что при представлении ее в Международной системе единиц SI в вектор-столбец искомых, а также в матрицу коэффициентов входят величины, которые имеют существенно отличающиеся порядки. Последнее усложняет решение системы ЛАУ методом Гаусса и последующее интегрирование значений производных методом Рунге-Кутта. Для выравнивания порядков величин программная реализация математической модели представляется в преобразованной системе единиц (ПСЕ), отличающейся от SI тем, что единица длины в ПСЕ принята в 102 раз большей единицы длины в SI. После выполнения расчетов на ПЭВМ все величины, выводимые на печать, представляются в системе SI.
Математическая модель (1) позволяет анализировать различные электромагнитные режимы реакторов, в частности, при соответствующем изменении значений активных сопротивлений обмоток, индуктивностей рассеивания обмоток или чисел их витков, реализуется несимметричный режим работы устройств, обусловленный несимметрией обмоток. Если варьировать значения геометрических размеров участков магнитопровода, реализуется несимметричный режим работы реакторов, обусловленный несимметрией магнитопровода. Если же изменять соответствующим образом значения напряжений подводимых к обмоткам устройства, то моделируется работа его при несимметрии подводимых напряжений. Задавая равными нулю число витков катушек той или иной обмотки, можно исследовать роль и воздействие этой обмотки на работу устройства. Изменяя величины сопротивлений нагрузки, можно переходить от режима холостого хода к режиму нагрузки, в том числе, и к несимметричному, и далее - к короткому замыканию.
Кроме того, на основе модели анализируются явления, возникающие при намагничивании магнитопровода устройств. Эти явления обусловлены насыщением магнитной цепи. Они сопровождаются появлением гармоник в намагничивающем токе реактора, то есть в напряженности магнитного поля, и (или) появлением гармоник в ЭДС, то есть в индукции поля. Наличие гармоник в магнитном поле отрицательно сказывается как на работе самого реактора, так и на работе энергосистемы, в которую включен реактор, а также на работе близлежащих линий связи. Насыщение магнитопровода реакторов может быть ~ 2 Тл и несколько больше, поэтому амплитуды 2, 3,...,12, 13-й, ... гармоник насыщения могут достичь заметных величин.
Анализ явлений намагничивания иллюстрирует качество устройств с точки зрения проявления высших гармоник, подтверждает достаточность мер, выработанных для их подавления и позволяет принять рекомендации, необходимые для проектирования устройств с учетом соответствующих гармоник. Результаты исследования могут быть использованы и в теории силовых серийных трансформаторов.
1.2. Анализ электромагнитных режимов НР
Ниже приводятся расчетные зависимости, полученные в результате реализации математической модели (1) на ПЭВМ.
На рис. 4,а показана зависимость линейных токов в фазах трехфазной совмещенной обмотки НР при следующих переходных процессах: 1) t0 = 0.00 с - включение НР на 3-фазное напряжение; 2) t1 =0.06 с - наброс напряжения; 3) t2 = 0.14 с - 10-процентный сброс его. Как видно из рисунка, 1-й процесс длится в течении ~ 0.02 с, 2-й и 3-й процессы заканчиваются за ~ 0.01 с. Как следует из приведенной зависимости, форма кривой тока является практически синусоидальной. В токе отсутствуют 5, 7, 11 и 13-я гармоники. Векторные диаграммы ЭДС по гармоникам этого порядка представляют собой замкнутые многоугольники, по пяти векторов в каждом, соответственно пяти катушкам, из которых состоит фаза обмотки. Отсутствуют в токе и гармоники, кратные трем. В токе проявляются 17 и 19-я гармоники, амплитуды которых невелики.
Рис. 4,б иллюстрирует переходный процесс при обрыве фазы A-X питающего напряжения (начало процесса соответствует моменту времени t1 = 0.06 с). При t > 0.06 с форма кривых токов в фазах B-Y и C-Z искажается за счет третьей гармоники, и нормальная работа НР нарушается.
На рис. 5,а приведено изменение магнитной индукции во времени в средних стержнях 1, 4 и 7 соответственно 1, 2 и 3-го модулей реактора. Уплощенная форма кривых свидетельствует о наличии в магнитном поле 3-й гармоники индукции.
На рис. 5,б представлена зависимость 3 гармоники индукции B3a, B3b, B3c = f(t), составляющей симметричную и синусоидальную трехфазную систему, т.к. в каждых трех стержнях 1, 4, 7; 2, 5, 8 и 3, 6, 9 первая и высшие гармоники индукции порядков 6·k ± 1, k = 1, 2, 3,.. взаимно уравновешиваются. Трехфазную систему ЭДС утроенной частоты можно снять с соответствующим образом выполненной вторичной обмотки. На этом принципе основан утроитель частоты, обладающий приемлемыми технико-экономическими показателями [6].
Благодаря быстрому параметрическому изменению реактивной мощности реактора в функции приложенного напряжения и отсутствию параллельных ветвей в обмотке, НР обладает практически безынерционным принципом действия, что делает особенно эффективным его использование для стабилизации напряжения в сети с резкопеременной нагрузкой. Вследствие пониженного содержания высших гармоник в кривой тока, а также других положительных качеств, НР является наиболее часто используемым типом реакторов.
2. Управляемый реактор электромашинного типа
Управляемые реакторы этого типа изготовляются на основе магнитопровода электрических машин переменного тока общепромышленного назначения. Активная часть реактора содержит магнитопровод и три обмотки: трехфазная пространственно распределенная рабочая обмотка (РО) и в общем случае две кольцевые обмотки управления (ОУ1 и ОУ2) (рис. 6) (фотография). Устройство и теория УР рассмотрены в [4, 9].
Трехфазная РО реактора может быть выполнена на напряжение не большее стандартного напряжения синхронного генератора, поэтому предпочтительнее использовать реактор в распределительных электросетях и в системах электроснабжения промышленных предприятий, т.к. в этом случае он подключается непосредственно, без промежуточного трансформатора. Реактор используется в качестве регулирующего элемента в схеме статического компенсатора реактивной мощности, в схеме параметрического стабилизатора напряжения (в этом случае обмотки управления не нужны, и устройство работает в режиме насыщающегося реактора) и др.
В реакторе органично реализуется (без усложнения устройства) режим симметричного намагничивания (СН), характеризующийся отсутствием четных гармоник пространственно-временного спектра, и в индукции, и в напряженности магнитного поля при подмагничивании магнитопровода постоянным током. Реактор в режиме СН выгодно отличается от реактора в режиме вынужденного намагничивания (ВН) рядом качеств. Для реализации режима СН необходимо, чтобы при равенстве размеров поперечных сечений внешнего и внутреннего ярм, равные МДС F0а и F0i обмоток управления, действовали встречно относительно контура вращающегося магнитного поля.
2.1. Математическая модель электромагнитных режимов УР
Математическая модель составляется на основе эквивалентирования реального устройства электрической (рис. 2 и рис. 7) и магнитной (рис. 8) схемами замещения с нелинейными сосредоточенными параметрами [4, 9]. Схема замещения трехфазной обмотки УР идентична схеме представленной на рис. 2. На рис. 7,а,б приведены схемы замещения кольцевых обмоток управления. Одна из них (ОУ1) уложена вокруг поперечного сечения ярма статора, другая (ОУ2) - вокруг сечения ярма ротора.
На рис. 7,а,б приняты обозначения: 1) U0а, U0i - значения напряжений постоянного тока, подведенных соответственно к обмоткам ОУ1 и ОУ2; 2) i0а , i0i - мгновенные значения токов подмагничивания (искомые величины); 3) R0а , R0i - активные сопротивления обмоток ОУ1 и ОУ2; 4) Ls0а, Ls0i - индуктивности, обусловленные магнитными потоками рассеивания обмоток ОУ1 и ОУ2 в переходном процессе; 5) pY0а = -e0а, pY0i - мгновенные значения э.д.с. индуктированных в обмотках ОУ1 и ОУ2 основными магнитными потокосцеплениями Y0а, Y0i этих обмоток в переходном процессе.
На рис. 8 показана схема замещения магнитной цепи УР. Приняты обозначения: 1) i1А, i1В, i1С - мгновенные значения токов фаз двухслойной трехфазной обмотки; 2) WА, WВ, WС - числа витков пазовых катушек фаз двухслойной трехфазной обмотки; 3) i0а, i0i - мгновенные значения токов обмоток ОУ1 и ОУ2); 4) W0а, W0i - числа витков пазовых катушек обмоток ОУ1 и ОУ2; 5) Ф1, Ф2 ,...,Ф13 - контурные магнитные потоки (искомые величины); 6) Rа1, Rа2,..., Rа12 - нелинейные дифференциальные магнитные сопротивления участков ярма статора; 7) Ri1, Ri2,..., Ri12 - нелинейные дифференциальные магнитные сопротивления участков ярма ротора; 8) Rz1, Rz2,..., Rz12 - нелинейные дифференциальные магнитные сопротивления зубцов; 9) Rd - линейное магнитное сопротивление радиального зазора между статором и ротором.
Математическая модель электромагнитных режимов реактора и явлений, возникающих при намагничивании магнитопровода, составлена аналогично модели для НР (раздел 1.1). Математическая модель - это система нелинейных ОДУ (3), включающая 18 уравнений. Неизвестными в ней являются тринадцать контурных магнитных потоков, три линейных тока двухслойной трехфазной обмотки и два тока кольцевых обмоток управления. Алгоритм и программная реализация математической модели принципиально не отличаются от таковых для НР.
2.2. Анализ электромагнитных режимов УР
Небольшое предварительное подмагничивание УР позволяет сократить время переходного процесса [4, 9]. Приводимые ниже зависимости, полученные в результате реализации математической модели (4) на ПЭВМ, соответствуют предварительно подмагниченному магнитопроводу.
На рис. 9,а,б приведены значения магнитной индукции Bа в поперечном сечении ярма статора и Bz в поперечном сечении зубца изменяющиеся в функции времени. Рис. 9,а соответствует режиму вынужденного намагничивания (ВН), а рис. 9,б - режиму симметричного намагничивания (СН). Как следует из анализа кривых, в режиме ВН в магнитной индукции содержатся четные гармоники, а в режиме СН они отсутствуют, что полностью согласуется с экспериментальными данными (рис. 10,а).
На рис. 11,а,б показана зависимость токов, протекающих в фазах трехфазной обмотки и в обмотках управления, для следующих переходных процессов: 1) t0 = 0.0 с - включение УР только на напряжения U0а, и U0i постоянного тока; 2) t1 = 0.1 с - включение УР на трехфазное напряжение; 3) t2 = 0.17 с - обрыв фазы A питающего напряжения. Рис. 11,а соответствует режиму ВН, а рис. 11,б - режиму СН при прочих одинаковых условиях. Как следует из сопоставления кривых реактор в режиме СН отличается следующими достоинствами: 1) в режиме СН токи подмагничивания и ток трехфазной рабочей обмотки принимают установившееся значение ~ через 0.015 с после подачи трехфазного напряжения, а в режиме ВН это время существенно больше; 2) диапазон регулирования тока в режиме СН больше, чем в режиме ВН примерно на 15...20 %; 3) отсутствием “трясущих” вибраций магнитопровода, что подтверждается данными эксперимента (рис. 10,б). При обрыве фазы А питающего напряжения в обоих режимах намагничивания в токах подмагничивания появляются четные гармоники, а в линейных токах - нечетные третьи гармоники, кратные трем, и работа реактора становится невозможной.
Изготовлены на основе пространственного витого магнитопровода производства Минского трансформаторного завода два совмещенных управляемых реактора-трансформатора на мощность 250 кВ·А, напряжение 10/0,4 кВ (фотография) [10]. На стороне высокого напряжения включены батареи конденсаторов 3x100 квар, 6,3 кВ. Устройства установлены на подстанциях в распределительных электросетях с целью выполнения следующих функций: трансформация напряжения, компенсация реактивной мощности, улучшение качества электроэнергии.
Выводы
1. Разработаны обобщенная математическая модель, ее алгоритм и программная реализация, электромагнитных режимов устройств с шихтованным магнитопроводом произвольной конструкции и уложенными на нем обмотками.
2. Проанализированы важные для практики нормальные установившиеся электромагнитные режимы работы реакторов; переходные процессы, возникающие при коммутациях в цепях обмоток; несимметричные режимы, обусловленные несимметрией питающей сети, самих устройств и нагрузки; явления, возникающие при намагничивании магнитопроводов; даны рекомендации по воздействию на качество электромагнитных процессов осуществлением различных режимов намагничивания и пр.
3. Установлено, что качественно электромагнитные процессы в реакторах трансформаторного и электромашинного типов протекают практически одинаково. Поэтому исследования реакторов этих типов осуществляются параллельно, дополняя друг друга, что способствует более полному изучению устройств, развитию их теории и улучшению технико-экономических показателей. Вместе с тем, выявлены количественные различия в характерах электромагнитных процессов, обусловленные схемотехническими и конструктивными особенностями реакторов указанных типов.
4. Адекватность математической модели подтверждается экспериментальными данными и определяется в основном принятыми допущениями, в соответствии с которыми устройства трактуются как электрические и магнитные цепи с нелинейными сосредоточенными параметрами. При принятых допущениях быстродействие программной реализации модели является вполне приемлемым.
5. При вводе в обобщенную математическую модель дифференциального уравнения движения, возможно исследовать, при соответствующей коррекции ее, не только электромагнитные, но и электромеханические режимы работы электрических машин, контакторов, реле и др.
Список литературы
1. Статические компенсаторы для регулирования реактивной мощности: Пер. с англ. / Под ред. Р.М. Матура. - М.: Энергоатомиздат, 1987. - 160 с.
2. Александров Г.Н. Передача электрической энергии переменным током. - Л.: Энергоатомиздат, 1990. - 176 с.
3. Забудский Е.И. Совмещенные электромагнитные устройства для систем управления в электроэнергетике // Всероссийский Электротехнический конгресс “На рубеже веков: итоги и перспективы” – Москва: МЭИ. – 1999. – C. 43–45.
4. Забудский Е.И. Математическое моделирование управляемых электроэнергетических устройств: Учебное пособие для вузов. - Ульяновск: УлГТУ, 1998. - 120 с.
5. Забудский Е.И. Расчет электромагнитных режимов работы ферромагнитных устройств электроэнергетического назначения: Пакет компьютерных программ / Зарегистрирован в ОФАП НИИВО Минобразования РФ. Регистр. №946999, 16.02.94.
6. А.с. 1663721 СССР, 5H02 M 5/16. Трехфазный статический ферромагнитный утроитель частоты / Е.И. Забудский, Ю.В. Ермураки, С.Ф. Козырин и др. (СССР).- № 4630114/07; Опубл. 15.07.91, Бюл. № 26.
7. Забудский Е.И. Стабилизация напряжения электросети на основе однокристальной микроЭВМ // Повышение эффективности функционирования систем электроснабжения электроэнергетики сельского хозяйства: Сборник научных трудов. – Москва: МГАУ, 2001. – С.50–69.
8. Мишин В.И., Забудский Е.И., Собор И.В. Трехфазные управляемые реакторы. - Кишинев: Штиинца, 1977. - 133 с.
9. Мишин В.И., Собор И.В., Забудский Е.И. Математическое моделирование статических ферромагнитных устройств. - Кишинев: Штиинца, 1981. - 104 с.
10. Забудский Е.И. Электрические машины. Ч. 1. Трансформаторы. Учебное пособие для вузов. – Москва: МГАУ, 2002. – 166 с.
Доктор техн. наук, профессор Забудский Евгений Иванович. Московский государственный агроинженерный университет, профессор кафедры Электроснабжение и Электрические машины. Россия, 127550, Москва, ул. Тимирязевская, 58.
Phone: (095) 407 3740, fax: 976 4396, e-mail: zei@inbox.ru. Web page: http://zabudsky.ru.
Вернуться назад
Возврат на начальную страницу