РАЗРАБОТКА МЕТОДА ЭКСПЕРИМЕНТАЛЬНОГО ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ СИНХРОННЫХ МАШИН

     Желанова Е.В., Ткаченко А.А.
     Донецкий национальный технический университет

     Введение. Широкое внедрение средств вычислительной техники и разработка различных математических моделей переходных процессов в электрических машинах переменного тока (ЭМ) создают предпосылки для повышения точности методов экспериментального определения их электромагнитных параметров. Для этой цели наиболее часто используют опыты внезапного трехфазного короткого замыкания (ВКЗ).

     Процедура нахождения значений параметров синхронных машин (СМ) из опытов ВКЗ регламентирована рекомендациями международной электротехнической комиссии [1] и действующим в Украине и странах СНГ ГОСТ 10169-77 [2]. Однако найденные по стандартной методике из опытов ВКЗ значения электромагнитных параметров имеют существенный разброс даже при обработке исходных данных с помощью ЭВМ [3]. Вопросам точности экспериментального определения параметров из опытов ВКЗ посвящено достаточное количество работ [4-7]. В [4,5] сделана попытка выяснить причины погрешностей в допущениях при упрощенном решении дифференциальных уравнений СМ, на которых основана стандартная методика и при получении самих уравнений. В качестве причин разброса в [6,7] выявлено влияние допущений, связанных с пренебрежением медленным вращением апериодической составляющей и неучетом второй гармоники тока ВКЗ. Сформулированы предложения по уточнению стандартной методики путем контроля точности проведения экспериментов и определения параметров отдельно по значениям токов в каждой фазе. Во всех рассматриваемых случаях предполагается, что СМ кроме обмотки возбуждения (ОВ) имеет только по одному демпферному контуру по каждой оси несимметрии ротора.

     При этом не оценивается влияние параметров СМ по поперечной оси на токи в фазах при ВКЗ. Как показано в [8], такое влияние всегда имеет место даже при КЗ из режима холостого хода (ХХ). Это обусловлено влиянием активных сопротивлений роторных контуров. Действующий в Украине стандарт [2] предусматривает также определение электромагнитных параметров и частотных характеристик с учетом многоконтурности ротора.

     Вопросы исследования математической модели, положенной в основу стандартной методики определения совокупности электромагнитных параметров СМ с учетом многоконтурности ротора по данным опытов ВКЗ в литературе не рассматривались.

     Целью настоящей работы является уточнение математической модели ВКЗ СМ, позволяющей повысить точность определения электромагнитных параметров и частотных характеристик проводимости со стороны обмотки статора с учетом многоконтурности ротора.

     Анализ стандартного алгоритма. Математическая модель, принятая в [6, 7], представляет аналитическое выражение для фазных токов. При КЗ из режима холостого хода, например, для фазы А обмотки статора оно имеет вид:

(1)

     где I_d, I'_d, I"_d, I₂, I_am — установившаяся, начальные значения переходной, сверхпереходной, тока двойной частоты и максимально возможной апериодической составляющих фазного тока; T_d, T'_d, T_a — постоянные времени затухания свободных токов; γ₀ — угол между осью фазы А статора и осью d ротора в момент КЗ; ω - частота вращения ротора; ω_a — частота вращения апериодической составляющей потока обмотки статора.

     Уравнение (1) представляет собой приближенное решение уравнений Парка-Горева при неучете активных сопротивлений обмоток статора и ротора на амплитуды и фазы соответствующих составляющих токов. Из анализа (1) следует, что токи ВКЗ отражают влияние только параметров СМ по продольной оси.

     Выражение (1) при отказе от учета тока двойной частоты и вращения апериодического потока статора является математической моделью, которой соответствует алгоритм нахождения электромагнитных параметров, описанный пунктами 17 и 24 ГОСТ 10169-77 [2].

     Выбор математической модели для определения электромагнитных параметров и частотных характеристик (ЧХ) с учетом многоконтурности ротора дан в [2]. В соответствии с пунктом 25.1 ГОСТ -10169-77 параметры машины по продольной оси следует определять по переходной функции, представляющей изменение огибающей периодического тока якоря из опыта ВКЗ. Эта составляющая определяется как полуразность ординат верхней и нижней огибающих фазного тока. Представляя периодическую составляющую тока суммой экспонент закон ее изменения во времени может быть описан функцией вида:

(2)

     где I_d — установившееся значение тока КЗ; I_d^k, T_d^k — начальное значение и постоянная времени затухания k-ой экспоненциальной составляющей периодического тока; n-количество экспонент, равное числу учитываемых контуров на роторе по продольной оси.

     Переходной функции (2) соответствует следующее выражение для определения ЧХ проводимости со стороны обмотки статора по продольной оси ротора y_d(js):

(3)

     Параметры входящие в (3) определяются следующим образом:

, .

     Соотношению (3) соответствует эквивалентная схема замещения с вынесенной на зажимы источника ветвью намагничивания, приведенная на рис. 1.

     Таким образом с учетом многоконтурности ротора закон изменения фазного тока статора ВКЗ имеет вид, аналогичный (1), с той лишь разницей, что периодическая составляющая тока представляется не двумя, а тремя-четырьмя (по рекомендации [2]) экспонентами. Тогда можно предположить, что все рекомендации по уточнению стандартной методики нахождения электромагнитных параметров, предложенные в [6, 7], могут быть распространены и для определения частотных характеристик с учетом многоконтурности. Однако, необходимость представления переходного периодического тока тремя-четырьмя экспоненциальными составляющими, требует дополнительной оценки погрешности математической модели, соответствующей выражению (1).

     Исследования проводились на примере расчета токов при трехфазном КЗ на выводах модельного синхронного генератора типа МТ-30 (S_ном = 30 кВА, U_ном = 230 В, I_ном = 47 А), имеющего демпферную систему и обладающего существенно различными динамическими свойствами по осям магнитной симметрии ротора. Параметры трехконтурных эквивалентных схем замещения по каждой оси ротора указанного генератора определялись методом затухания постоянного в обмотке статора неподвижной машины.

     На рис. 2 приведены результаты расчета тока статора при переходном процессе вызванном ВКЗ на выводах генератора МТ-30 при использовании различных математических моделей.

     Сопоставление результатов моделирования, выполненного путем решения полных уравнений Парка-Горева (кривая 1) с расчетом по выражению (1) указывает на их несовпадение (кривая 2). Если в области максимальных мгновенных значений тока погрешность не превышает 6,2%, то в области минимальных — она достигает 44%. Использование выражения (1) при неучете составляющей двойной частоты дает качественно неверные результаты (кривая 3).

     Для повышение точности необходимо получение математической модели, основанной на описании закона изменения изображающего вектора тока статора. В этом случае регистрация переходных функций возможна с любой дискретизацией.

     Уточненное определение переходного периодического тока статора СМ с помощью частотных характеристик. Получим уточненную модель на основе установленных в [8] и уточненных для АМ в [9] связях между частотными характеристиками, эквивалентными схемами замещения и составляющими переходного тока при коротких замыканиях.

     Исследование переходных процессов с помощью частотных характеристик в СМ имеет свои особенности вследствие наличия асимметрии ротора. Это приводит к необходимости характеризовать машину двумя частотными характеристиками по продольной и поперечной осям магнитной симметрии ротора. Раcсмотрим определение составляющих токов при КЗ на выводах СМ. Будем предполагать, что на роторе имеется несколько демпферных контуров и ротор имеет магнитную и электрическую несимметрию.

     При КЗ в обмотке статора будут протекать: установившийся — i_-s0, апериодический — i_-s1 и переходный периодический — i_-s2 токи (рис.3).

     Математическая модель СМ представляет аналитическое выражение:

(4)

     В (4) отдельные составляющие переходного тока статора определяются по следующим соотношениям:

(5)

(6)

(7)

     где δ₀ — угол, характеризующий положение поперечной оси ротора по отношению к вектору напряжения, подводимого к обмотке статора.

     При КЗ из режима холостого хода δ₀ = 0.

     На рис. 2 приведены результаты расчета по полученной математической модели (кривая 4).

     Сопоставление с результатами расчета, выполненного путем решения полных уравнений Парка-Горева (кривая 1) свидетельствует об их практическом совпадении. Отличие в течение всего переходного процесса не превышает 1%. Аналогичные исследования, проведенные для синхронной машины СГ-18,75 показали практическое совпадение результатов расчета по уравнениям Парка-Горева, соотношению (1) и предложенной методике. Отличия не превышают 3%.

     Это обусловлено тем, что вследствие малого значения активного сопротивления ОВ вектор переходного периодического тока i_-s2 практически не имеет проекции на поперечную ось q.

     Вывод

     Полученная в работе уточненная модель синхронной машины позволит повысить точность определения электромагнитных параметров с учетом многоконтурности ротора по данным обработки переходных функций токов статора при внезапном трехфазном коротком замыкании.

Список литературы