Определяется среднее значение установившегося тока для скольжения

, (1)

где .

Рассчитывается пульсирующая составляющая установившегося тока

, (2)

где ;   - угол между поперечной осью ротора и вектором напряжения сети.

Изменение во времени полного установившегося тока статора с учетом различия параметров ротора по осям магнитной и электрической симметрии

,(3)

Определяется апериодическая составляющая

(4)

и периодический ток частоты близкой к двойной

. (5)

Изменение во времени апериодического тока и тока, обусловленного асимметрией ротора, подчинено следующему закону

, (6)

Собственная частота вращения апериодического тока и постоянная времени его затухания определяются по значению среднего комплекса проводимости  для скольжения .

.     (7)

. (8)

Начальное значение переходного периодического тока  рассчитывается из условия

.

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

, (9)

где  - начальные значения и постоянные времени экспоненциальных составляющих, аппроксимирующих огибающие периодических токов по продольной и поперечной осям соответственно; N - количество контуров на роторе по продольной оси; M – то же по поперечной оси.

Значения параметров указанных периодических составляющих тока  в условиях дискретного задания частотных характеристик целесообразно рассчитывать с использованием эквивалентных схем замещения, исследование которых приведено в [5]. При совпадении вектора тока  с продольной или поперечной осью ротора его затухание будет определяться только параметрами соответствующей схемы замещения:

,

. (10)

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

,

, (11)

где -аргумент вектора тока  (угол между вектором тока и осью действительных чисел комплексной плоскости. Тогда:

,

. (12)

Полный вектор изображающего тока статора представляет собой векторную сумму отдельных составляющих:

, (13)

где U - напряжение на выводах статора в опыте включения невозбужденного генератора в сеть.

Мгновенные значения фазных токов определяются как проекции изображающего вектора на неподвижные оси времени соответствующих фаз:

,

, (14)

где  - аргумент полного изображающего вектора тока статора, зависящий от времени.

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

            Электромагнитный момент в соответствии с общим подходом к его определению по значениям комплексов токов и потокосцеплений будет определяться для ЭМ с симметричным или асимметричным ротором следующим выражением:

(15)

Потокосцепления в (15) имеют следующие значения:

. (16)

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

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

, (17)

где М_мех – внешний механический момент, приложенный к ротору;  T_j – механическая постоянная инерции.

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

- Весь процесс разбивают на ряд равномерных интервалов по времени.

- На каждом интервале в соответствии с (17) находят изменение скольжения

.

- Установившуюся составляющую тока  рассчитывают как при постоянной скорости, подставляя для каждого момента времени действительную скорость вращения ротора в качестве переменного параметра.

- Амплитуду составляющей , затухающей с постоянной времени , рассчитывают для значения скорости в начальный момент времени.

- Собственную частоту вращения  и постоянную  пересчитывают по выражениям (7) и (8) в зависимости от изменения скольжения.

- Амплитуду составляющей тока , обусловленной переходными токами в роторных контурах, определяют для значения скольжения в начальный момент времени.

- Скорость вращения тока , учитываемая в (9), изменяется при изменении скольжения.

- qУчитывается зависимость от скольжения собственной частоты вращения и постоянной времени затухания потокосцепления  в (16).

Необходимым условием для реализации предложенного способа является аналитическое описание исходных частотных характеристик и закона изменения во времени переходной периодической составляющей тока статора. В [5] предложена методика использования для этих целей эквивалентных схем замещения.

Эффективность применения частотных характеристик к анализу динамических режимов исследовалась на примере расчетов следующих переходных процессов:

1. Включение в сеть невозбужденного модельного турбогенератора типа МТ-30. Указанный СГ обладает существенно различными динамическими свойствами по осям магнитной симметрии ротора и имеет следующие параметры: .

При этом анализировалось влияние величины и знака скольжения в момент подключения СМ к сети, а также влияние изменения скорости вращения ротора в переходном режиме.

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

3. Для асинхронного электродвигателя типа ДАЗО-1914-10/12А (Р_ном = 1500 кВт; I_ном = 204 А; U_ном = 6000 В) были выполнены расчеты электромеханических переходных процессов при пуске из ненагруженного состояния.

Оценка степени достоверности результатов расчетов по предлагаемой в работе методике осуществлялась путем сравнения с соответствующими результатами расчетов, выполненных по алгоритмам, реализующим численное интегрирование полных дифференциальных уравнений Парка-Горева. Частотные характеристики для исследуемых ЭМ определялись по методу затухания постоянного тока в обмотке статора.

Программы, реализующие алгоритмы расчета электромеханических переходных процессов, была разработаны  в среде математической системы MathCAD 7.0.

Расчеты зависимости изменения изображающего тока статора и электромагнитного момента вращения при включении невозбужденного модельного генератора МТ-30 в сеть для различных начальных скольжений ротора выполнялись без учета и с учетом изменения скорости вращения ротора после подключения СГ к сети. Постоянство скорости вращения моделировалось путем задания механической постоянной инерции  большой величиной. Как следует из сопоставления результатов расчета, выполненных для скольжений 0; 0,01; 0,03; 0,05 и 0,1о.е., изменение зависимостей изображающего тока статора и электромагнитного момента практически совпали с соответствующими расчетами, выполненными по уравнениям Парка-Горева при постоянном скольжении. При этом отличия в максимальных значениях сопоставляемых параметров не превышает 10,4%. Характер процессов оказался абсолютно идентичным. Это позволяет утверждать, что электромагнитные процессы, отображаются при использовании статических частотных характеристик достаточно точно.

Анализ результатов расчета с учетом изменения скорости вращения ротора после включения СГ в сеть указывает на некоторые отличия в зависимостях изменения параметров динамического режима. На рис.1, в качестве примера, приведены результаты расчета переходных процессов при включении в сеть невозбужденного модельного синхронного генератора типа МТ-30 при короткозамкнутом состоянии обмотки возбуждения.

Рисунок 1 - Переходные процессы при включении в сеть невозбужденного модельного синхронного генератора типа МТ - 30

В расчетах, выполненных по уравнениям Парка-Горева (кривые 2), отражается динамика процесса, связанная с изменением токов и электромагнитных моментов при изменении скорости вращения. Изменение скольжения при этом носит колебательный характер. Моделирование, использующее статические характеристики (кривые 1), приводит к монотонному изменению скольжения ротора в переходном режиме. Однако, общее время приближения скорости к синхронной, оцениваемое по времени второго прохождения скольжения через ноль, практически совпало со временем достижения ротором синхронной скорости, полученное при моделировании по предложенной методике. В целом характер изменения исследуемых параметров режима совпал. Максимальные их значения, возникающие в начальной стадии переходного процесса, практически не отличаются от данных, полученных без учета изменения скорости.

На рис.2 приведены кривые изменения параметров режима при пуске асинхронного двигателя типа ДАЗО–1914–10/12А.

Рисунок 2 - Переходные процессы при пуске ненагруженного асинхронного двигателя типа ДАЗО - 1914 - 10/12А

Параметры трехконтурной схемы замещения с вынесенной ветвью намагничивания, синтезированные по [6], имеют следующие значения (о.е.):

;     ;      ;      ;

               ;    ;   .

Сопоставительный анализ результатов расчета, выполненных различными методами, свидетельствует об их удовлетворительном совпадении (с точки зрения длительности процесса пуска). Время пуска, определенное по предложенной методике (кривые 1), отличается от рассчитанного по полным уравнениям Парка-Горева, на 9,2%.

Переходные процессы в течение времени, определяемым протеканием электромагнитных явлений, оказались близкими по характеру и значениям. В области времен от 0,5 с до 1,8 с. параметры исследуемого режима практически совпадают между собою. Это, очевидно, объясняется тем, изменение скорости вращения ротора происходит с малым ускорением и использование статических характеристик практически не вносит погрешность в динамический процесс. В области времен, когда скорость вращения ротора близка к номинальной, в характере изменения тока статора и электромагнитного момента имеются существенные отличия. Исследования показали, что эти отличия уменьшаются при увеличении  и коэффициента загрузки двигателя, поскольку в этом случае процессы сопровождаются меньшими ускорениями ротора.

Исследования, выполненные при анализе переходных процессов, вызванных включением невозбужденного турбогенератора типа ТГВ-200 в сеть, для пятиконтурных и трехконтурных эквивалентных схем замещения [5], показали, что погрешности, вносимые сокращением количества контуров в изменения параметров режима во временной области, не превышают погрешностей, в частотной области, и составляют 7-11%.

Выводы

1. Предложенные математические модели ЭМ, основанные на экспериментальных статических частотных характеристиках, позволяют без решения дифференциальных уравнений исследовать динамические переходные режимы с учетом изменения скорости вращения.

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

Литература

1. Казовский Е.Я. Переходные процессы в электрических машинах переменного тока. - М.: Изд-во АН СССР, 1962.-624 с.

2. Применение частотных методов в электроэнергетическом исследовании. Труды СибНИИЭ / Под ред Лукашова Э.С. – Новосибирск, 1972, вып. 21. – 294 с 

3. Рогозин Г.Г. Определение электромагнитных параметров машин переменного  тока. – К.: Техніка, 1992. – 168 с.

 4. Ларин А.М., Ларина И.И., Гуедиди Фаузи.Упрощенная модель синхронной машины для расчета трехфазных коротких замыканий с учетом многоконтурности и несимметрии ротора. Научные труды ДонГТУ, серия Электротехника и энергетика, вып.4. Донецк, 1999. – С. 60-64.

5. Ларин А.М., Ламари Абдессалем.Определение параметров схем замещения турбогенераторов, отображающих электромагнитную несимметрию их массивных роторов. Научные труды ДонГТУ, серия: Электротехника и энергетика, вып.4. Донецк, 1999. – С. 65-68.

6. Рогозин Г.Г., Ларин А.М. Расчет параметров эквивалентных роторных контуров синхронных машин по их экспериментальным частотным характеристикам  /  Электричество, 1974,  № 6.  С.10 - 13.