Усреднить тангенс наклона касательной на шаге интегрирования

можно следующим образом (рис. 3-2). С помощью метода Эйлера находится точка (t _k + h, x_k + hx^′_k), лежащая на прямой L₁. В этой точке вновь вычисляется тангенс угла наклона касательной (пря­мая L₂). Усреднение тангенсов дает прямую L. Через точку (t_k, x_k) проводится прямая L, параллельная L, и находится искомая точка (t_k+1, x_k+1 ).

В этом случае, очевидно, функция Ф [см. (3-7)] запишется как

Ф = 0,5[ f(t_k, x_k) + f(t_k + h, x_k + hx^′_k)], (3-8)

где x'k==f(tk. Xk)-

Можно поступить иначе, определив тангенс наклона касатель­ной в средней точке интервала экстраполяции (шага интегрирова­ния), как это показано на рис. 7-3. В этом случае

Ф = f(t_k + 0,5 h, x_k + 0,5hx^′_k). (3-9)

Без вывода приведем формулы, описывающие метод четвертого порядка, один из самых употребительных методов интегрирования дифференциальных уравнений. Этот метод применяется настолько широко, что в литературе по вычислениям на ЦВМ он называется «метод Рунге—Кутта» без всяких указаний на тип или порядок. Этот классический метод Рунге—Кутта описывается системой сле­дующих соотношений:

x_k+1 = x_k + h(k₁ + 2k₂ +2k₃ +k₄)/ 6,

где

k₁ = f (t _k, x _k); k₂ = f (t_k + 0,5 h, x_k + 0,5 h k₁); k₃ = f(t_k + 0,5hk₂);

k₄ = f(t_k + h, x_k + hk₃).

Данный метод самоначинающий, но один из его недостатков—

сложность оценки локальной ошибки интегрирования -является об­щим недостатком всех методов Рунге —Кутта. Поэтому при реализации метода с автоматическим выбором шага интегрирования поступают следующим образом. Вычисляют значение искомой функции в каж­дой точке с полным и с половинным шагом и сравнивают два по­лученных значения. Если модуль разности между ними окажется меньше некоторого заданного положительного числа, то считается, что интегрировать с данным шагом можно; если же это не так, то шаг делится пополам. Этот метод выбора шага интегрирования при­водит к увеличению времени счета как минимум в три раза по сравнению с методом постоянного шага.

Некоторые особенности анализа переходных процессов в элект­рической системе. Рассмотрим особенности анализа переходных про­цессов в электрической системе с точки зрения метода решения за­дачи:

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

2) порядок системы уравнений, описывающих поведение слож­ной электрической системы, высок, что, с одной стороны, приводит к естественному увеличению времени счета, а с другой,—к столь же естественному увеличению объема оперативной памяти ЦВМ, нани­маемого программой;

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

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

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

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

T_{J I} d²δ_i/dt² = P_{T i} -- P_Эi, (5—10) где

Рэ_i = y_ii E²_i sin α_ii -- ∑ y_imE_i E_m sin (δ_i –δ_m –α_im), i=l, 2,..., η.

От системы n уравнений второго порядка (7-10) можно перейти к системе 2п уравнении первого порядка

dδ / dt = ω_i; T_J dω/ dt = P_Ti -- P_Эi. (5-11)

Предположим, что в системе происходит короткое замыкание. Тогда, если в исходном режиме δ_i = δ_{i 0} , ω = 0, ΔP = 0 для каж­дого генератора, то после короткого замыкания ΔР не равно О. При этом изменяются и значения проводимостей схемы, а следовательно, и значения Рэ и ΔР==Рг—,Рэ. Обозначая ускорение ротора i-го генератора через α=ΔP/Τ_J и считая его постоянным, запишем приращение угла генератора ηа первом шаге интегрирования h

Δδ_i.= 0,5α₀h².

Применяя для каждого генератора формулу δι=δο+Δδι, найдем значения углов δ к началу второго шага интегрирования. Очевидно, что ко второму шагу интегрирования ω не равна нулю. Прираще­ние ω на первом шаге интегрирования определим по среднему ускорению генератора, т. е. по формуле

Δω₁ == 0,5h (α₁+ α₀),

где α = ΔP₁ / T_J..

k-ι ••Σ s=l

Для второго интервала

Δδ = Δω₁h + 0,5 α₁h².

Таким образом, интегрирование для всех последующих интервалов производится по формулам:

ω_k = ω_k—1 + 0,5h (α_k + α_k—1); (3—12)

δ_k+1 = δ_k + ω_k h + 0,5α_k h², (3—13)

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

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

ω_k = h/2∑(α_S + α_S—1) = 0,5α_kh + h∑ α_s + 0,5 α₀h;

δ_k+1 = δ_k + h² ∑α_s + 0,5α₀ h². (5—14)

Очевидно, что

δ_k = δ_k---1 + Δδ_k = δ_k—1 + h²∑ α_s + 0,5α₀h²,

и Δδ_k = h² ∑α_s + 0,5α₀ h²

Следовательно, (5-14) можно переписать:

δ_{k +1}= δ_k + Δδ_k + α_kh².

Δδ_{k+ 1} = Δδ_k+ α_kh^2;

δ_k+1= δ_k + Δδ_k+1;

α_k = 18000ΔP_k/T_J,

где углы δ измеряются в электрических градусах, время интегри­рования и постоянная инерции T_J—в секундах, небаланс мощно­сти Δρ—в относительных единицах. Алгоритм метода следующий;

1/Установить m = 1 2/ Определить константу k = (18000/Tj} h².

3.ВычислитьΔΡ₀= Рт-Рэ (Пункты 1 – 5 выполняются однократно)

вы-одно-

4. Определить Δδ₁==0,5kΔΡο.

  кратно

5. Найти δ₁=δ₀ + Δδ₁

6. Вычислить ΔР_m=Р_T—Рэ.

7. Определить Δδ_{m+ 1}=Δδ_m +kΔP_m.

8. Найти δ_m+1 = δ_m+ δ_m+ Δδ_m+1.

9. Увеличить m на единицу {m — >m+ 1}.

10. Если mh не больше заданного времени интегрирования, то перейти к п. 6, если mh больше, то счет оканчивается.

Данный алгоритм справедлив как во время короткого замыкания, так и после его ликвидации. Однако непосредственно в момент отключения ускорение α претерпевает разрыв первого рода, по­этому п. 6 алгоритма для моментов коммутаций следует изменить и определить небаланс мощности по формуле

ΔΡ_m=0,5(ΔΡ_{m +} +ΔΡ_M-),

где индекс «+» свидетельствует о том, что ΔΡ необходимо вычислить при послеаварийных значениях проводимостей у, а «_» о том, что значения проводимостей отвечают аварийному состоянию в системе (короткое замыкание).

Данный метод можно использовать и в тех случаях, когда ис­следуется устойчивость регулируемой системы с учетом изменения ЭДС генераторов. Однако в этом случае он уже не является мето­дом второго порядка, так как «новые» значения ЭДС определяются по формуле Эйлера.

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

Для прогноза в данном методе используется формула Эйлера

x_k+₁ = x_k + hf (t_k,x_k),

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

x^(m+1)_k+1 = x_k+ 0,5h[ f(t_k,x_k ) + x′_{k +1}], (3-15)

где x′_k+1 = f(t _k + h, x^m_k+1).

Данный метод имеет сходство с методом Рунге—Кутта второго порядка. Однако при использовании метода Рунге—Кутта формула (5-15) применяется во всех случаях один раз, а при использовании данного метода корректировка проводится до тех пор, пока не вы­полнится условие

[ x^m+1_k+1 -- x^m_k+1 ] ‹ έ (3-16)

где ε—заданная положительная величина.

В данном методе в зависимости от того, сколько итераций не­обходимо сделать для выполнения условия (5-16), можно выбирать шаг интегрирования: если этих итераций много—шаг следует умень­шить, если итерация одна, то шаг можно увеличить.

3-2. ПОНЯТИЯ КАЧЕСТВЕННЫХ ХАРАКТЕРИСТИК ДЛЯ ВОЗМОЖНЫХ ВИДОВ ДВИЖЕНИЯ НЕЛИНЕЙНОЙ СИСТЕМЫ

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

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

Рассмотрим понятие фазовой плоскости для системы, движение которой описывается уравнением второго порядка

d^2.x./ dt² + ω₀² x = 0. (3-17)

Состояние системы в данный момент времени t ₀ определяется двумя координатами:

(x)_{t = t0} = x₀; (dx/dt) _{t =t0} = (dx/dt)_0.

Примем величины x и y = dx/dt в качестве декартовых коорди­нат на плоскости (рис. 3-4). Точка.М на плоскости, имеющая координаты x₀ и (dx/dt)₀ будет соответствовать некоторому состоя­нию системы, движение которой описывается уравнением (5-17). Точка М, положение которой полностью характеризуется состоя­нием системы, называется изображающей точкой, а плоскость (.x, dx/dt} — фазовой плоскостью. При движении системы ее состоя­ние изменяется и изображающая точка М, находящаяся в началь­ный момент времени t==t₀ в положении М₀, описывает некоторую кривую (см. рис. 5-4), которая называется траекторией изобра­жающей точки на фазовой плоскости. Запишем решение уравнения (3-17)

x = A sin (ω₀ t + φ). (3-18)

Здесь А и φ—постоянные, определяемые из начальных условий.

Дифференцируя выражение (3-18), находим

у = dx/dt = Αώ_o cos (ω_o t + φ). (3-19)

Исключая величину ω_ot +φ из (3-18) и (3-19), получаем урав­нение

x²/A² + y²/ (A²ω_o² ) = 1,

представляющее собой уравнение эллипса с полуосями A и ω_o Α.

Итак, если система находится в движении, описываемом урав­нением (5-17), то изображающая точка движется по эллипсу (рис. 5-5). При этом время нигде не отображено на фазовой плоскости и можно лишь представить себе, что с течением времени изображающая точка Μ будет двигаться по траектории в направлении, указанном стрелкой. При других начальных условиях частное решение урав­нения (5-17) дает другую траекторию, например траекторию, на которой находится точка N.

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

Рис. 3-6

н ачальных условий имеет лишь одно определенное движение, а не 'два или более «возможных» движений. Исключение составляют особые точки на фазовой плоскости, соответствующие состояниям равновесия системы. В нелинейной системе может быть несколько состояний равновесия.

Так, в нелинейной системе, описываемой уравнением второго порядка

d²x/dt² + 2βdx/dt + ω_o² sin x = ω_o² sin x_o,

может быть два состояния равновесия, определяемые как sinx= sin x_o,

где х==x_o и х=л – x_o.

Первое состояние равновесия (2β)² < ω_о ² называется устойчивым

фокусом (рис. 3-6, а).

Второе неустойчивое состояние равновесия называется седлом (рис. 3-6. д). На рис. 3-6, б, в, г приведены и другие возможные

a

Рис. 3-7

особые точки- неустойчивый фокус, устойчивый узел и неустойчи­вый узел соответственно.

Исследуя характер движения нелинейной системы на фазовой плоскости для распознания всей картины траекторий, т. е. так на­зываемый фазовый портрет системы, необходимо уделить особое внимание особым траекториям. Заметим, что точки равновесия являются частными случаями особых траекторий.

Периодические движения, при которых через каждый период времени Т исходное состояние системы повторяется, представляются замкнутыми фазовыми траекториями. В консервативных системах кривые, соседние с периодической кривой, обычно также пери­одические. Характерной картиной при существовании периоди­ческих движений в консервативной системе является совокупность замкнутых траекторий, «вложенных» одна в другую (см. рис. 3-5). Нет никаких оснований назвать какую-либо из этих траекторий особой. Однако в неконсервативных системах замкнутые траектории имеют принципиально другое значсние. Здесь существуют только изолированные замкнутые траектории, представляющие собой осо­бые траектории, называемые предельным циклом.

На рис. 3-7, а, б приведены примеры устойчивого и неустой­чивого предельных циклов. Начальная точка М, лежащая в области,

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

Малое отклонение начальной точки, лежащей вокруг или на самом неустойчивом предельном цикле (рис. 5-7. б) приводит к уда­лению движения системы от периодического. Поэтому движение но неустойчивому предельному циклу физически нереализуемо. Устой­чивые и неустойчивые предельные циклы служат «водоразделами»,

Рис. 3.8

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

При достаточно большом удалении от окрестности седла «усы» также служат «водоразделами» между движениями различного ха­рактера и называются сепаратрисами (кривая С =V_гр на рис. 3-8). Особые траектории— точки равновесия, сепаратрисы, предельные циклы разбивают всю фазовую плоскость на ряд областей, которые определяют полную качественную характеристику всех возможных видов движения системы. Для электрической системы с устройст­вами автоматического регулирования характерны все описанные виды особых траекторий. Сепаратриса разделяет затухающие пери­одические движения и расходящиеся апериодические движения. Она имеет важное техническое значение, так как область начальных условий, находящихся внутри сепаратрисы, определяет условия устойчивости после больших возмущений (динамической устойчивости электрической системы).

3-3. НАХОЖДЕНИЕ ДОСТАТОЧНЫХ УСЛОВИЙ УСТОЙЧИВОСТИ ПРИ БОЛЬШИХ ВОЗМУЩЕНИЯХ С ПОМОЩЬЮ ПРЯМОГО (ВТОРОГО) МЕТОДА ЛЯПУНОВА

Для определения устойчивости после больших возмущений (дина­мической устойчивости электрической системы) необходимо знать совокупность начальных отклонений параметров режима и их производных, при которых система вернется в положение равновесия (область устой­чивости). Эта совокупность ограничена в фазовом пространстве сепаратрисной поверхностью. Для отыскания поверхности, лежащей внутри области устойчивости и прибли­жающейся к сепаратрисной поверхности, можно применить прямой (второй) ме­тод Ляпунова, заключающийся в следующем. Обобщив известное положение физики (равновесное положение устойчи­во, если в нем потенциальная энергия имеет минимум), Ляпунов предложил на­ходить при выводе условий устойчивости вспомогательную функцию координат х_i изображающей точки в фазовом простран­стве V{x_i)- Эта функция должна быть одно­значна, дифференцируема, определенно по­ложительна вне положения равновесия и обращаться в нуль в по­ложении равновесия.

Если удается подобрать такую функцию V, что при любом дви­жении системы она уменьшается, т. е. dV/dt < 0, то система оказы­вается «устойчивой в большом».

Примером функции Ляпунова может служить функция

v = x₁² + x₂² +x₃².

В этом случае поверхность V =С в фазовом пространстве (рис. 5-9) есть сфера. Поверхности V=C, построенные для различных значе­ний С(С₁, C₂, С₃), не пересекаются, так как — V --

Рис. 3-9

однозначная функ­ция координат. Поэтому семейство представляет собой вложенные одна в другую замкнутые поверхности, причем поверхность с мень­шим С вложена в поверхность с большим С. Поверхность V=С при С— >0 стремится к точке равновесия О.

Пусть в какой-либо момент времени изображающая точка М₁ (рис. 5-9) находится на поверхности V › C₃; если при движении системы производная dV/di < 0, то с течением времени изображаю­щая точка перемещается к поверхности V =С₂ , где С₂< С₃. Следо­вательно, фазовая траектория пронизывает вторую поверхность, лежащую внутри первой. Если изображающая точка при движении пронизывает поверхности семейства в направлении снаружи внутрь и с течением времени неограниченно приближается к точке равно­весия, то выполняются достаточные условия устойчи­вости в большом. Однако условие dV/dt < 0 не является необ­ходимым признаком устойчивости, так как могут существовать

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

Недостатки прямого метода Ляпунова—отсутствие общих приемов

отыскания функции Ляпунова и невозможность оценки, насколько достаточные условия устойчивости в большом (С_i < V_гр, рис, 5-8) уже, чем необходимые условия устойчивости, определяемые сепара­трисой С=V_гр· Одно из достоинств метода Ляпунова—сокращение времени численного счета переходных процессов после приложения большого возмущения для определения динамической устойчивости.

Численными методами рассчитывается переходный процесс от момента приложения до момента снятия возмущающей силы (на­пример, такая кратковременная аварийная ситуация, как короткое замыкание в электрической системе, длится 0,1—0,4с).

Если значения переменных послеаварийного режима, определен­ные в конце аварийной ситуации, находятся внутри фазовой пло­скости, ограниченной предварительно построенной функцией Ляпу­нова, то система устойчива в большом (или динамически устойчива).

Пример 3-1. Найти прямым методом Ляпунова условия устойчивости в боль­шом электрической системы, представленном схемой «станция - шины».

Запишем нелинейное дифференциальное уравнение движения синхронной ма­шины, соединенной с системой, в простейшей (консервативной) идеализации:

dδ/dt=ω; dω/dt=P_o/T_J—[EU/{xT_J)] sinδ. (3-20) Введем обозначения:

A=P_o /T_J; b=eu/(x t_J).

Найдем функцию Ляпунова. Для электрической системы в консервативной идеализации при выборе функции Ляпунова используем энергетический подход, Функцию Ляпунова запишем в виде первого интеграла движения синхронной ма­шины—интеграла энергии. Разделив второе уравнение (5-20) на первое,получим

(dω/dt)(dδ/df) = (A—Bsinδ)/ ω (5-21)

Разделяя переменные и интегрируя уравнение (5-21), найдем

0,5ω² – Аδ – В cos δ = C. (5-22)

Выражение (5-22) характеризует полную энергию системы (5-20). состоящую из кинетической Ο.5ω² и потенциальной Aδ + B cos δ. Траектории рассмотренной электрической системы на фазовой плоскости [δ, ω] являются контурами (кри­вые C_i на рис. 5-8) постоянной полной энергии. Изменение полной энергии по времени запишем как

dC/dt = (дС/дω) (dω/dt) + (дС/дδ) (dδ/dt).

Используя (3-22) и (3-20), найдем

dC/dt=ω(A—B sinδ)—Aω)+Bωsinδ=0. (3-23)

Следующий из (3-23) вывод о равенстве нулю изменения полной энергии в контуре С физически означает, что система, переведенная в результате какого-либо возмущения из точки 0 в точку М, лежащую на траектории со значением Сi (рис. 3-8), будет совершать движение по этой траектории. Если эта траектория замкнутая, то движение будет устойчивым. Так как функция (3-22) знакоопределенная, то ее можно принять за функцию Ляпунова, поскольку она удовлетворяет требованиям, предъявляемым к функциям этого класса, а именно:

1) V (δ, ω) > О в некоторой области вокруг начала координат, совмещенного с устойчивым положением равновесия;

2) dV/dt > 0 во всем пространстве.

Областью устойчивости системы (3-20) является часть окружающего начало координат пространства,ограниченного самой большой замкнутой поверхностью C ==Vrp. соответствующей граничному значению функции Ляпунова. Эта гранич­ная—сепаратрисная—поверхность проходит через седловую точку (точка О₂ на рис. 5-8), которую необходимо определить.

Как известно, особые точки находятся приравниванием нулю правых частей исходных дифференциальных уравнений переходного процесса (5-20):

ω=0; A -- Bsinδ=0; δ_о=arcsin {А/В).

При А < В имеется два положения равновесия—центр Ο₁ [δ_о, 0] и седло О₂ [π—δ_о, О) (рис. 5-8).

Замкнутая часть сепаратрисной поверхности Vrp, проходящей через седло, определяет область устойчивости системы в большом Критерий устойчивости математически запишется как V = V_гр,

3-4. ИССЛЕДОВАНИЯ ПЕРИОДИЧЕСКИХ РЕШЕНИЙ С ПОМОЩЬЮ ГАРМОНИЧЕСКОЙ ЛИНЕАРИЗАЦИИ НЕЛИНЕЙНОСТЕЙ

Построение длительных колебательных процессов численными методами решения нелинейных дифференциальных уравнений может дать большую ошибку; кроме того, при этом невозможно получить какие-либо общие закономерности. Эти обстоятельства заставляют искать особые приемы для приближенного определения периодиче­ских решений нелинейных дифференциальных уравнений, характер­ных для синхронных машин в электрической системе. В этих при­емах используется разложение нелинейности в ряд Фурье при отыскании периодических решений и совокупного решения уравне­ний для каждой гармоники (метод гармонического баланса).

В случаях, когда амплитуды колебаний высших гармоник по от­ношению к первой малы и когда в ряду Фурье можно оставить только постоянные составляющие и первые гармоники, метод гар­монического баланса превращается в метод гармонической линеари­зации. Заметим, что условиями малости составляющих высших гармоник выполняются, если система обладает свойством фильтра. Объект исследования электрической системы—синхронная машина—обла­дает свойством фильтра тем сильнее, чем меньше ее демпфирование (амплитудно-частотная характеристика колебательной системы имеет резко выраженный резонанс, см, рис. 3-23). Цель метода гармони­ческой линеаризации состоит в том, чтобы вычислять параметры колебательного процесса—координаты средней точки колебаний (постоянные составляющие), максимальные отклонения от нее (ам­плитуды первой гармоники), время одного полного колебания, пе­риод или частоту. Однако предварительно нужно определить, будет ли процесс при данных параметрах системы, данном внешнем воз­действии и данных начальных условиях колебательным или моно­тонным. Иными словами, надо найти границы существования пери­одических решений. Для синхронной машины, работающей в элек­трической системе, синхронные колебания ограничиваются условиями

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

Пример 5-2. Определить максимально допустимую амплитуду колебаний не­линейной колебательной системы, превышение которой приводит к неустойчивости в большом. Применительно к синхронной машине, работающей параллельно с электрический системой, появление колебания ротора (либо свободного—вслед­ствие начального возмущения, либо вынужденного—вследствие внешней гармо­нической силы) с амплитудой выше критической oзначает нарушение синхронной динамической устойчивости, т, е. неустойчивости в большом. Определение макси­мально допустимых амплитуд колебаний, таким образом, является практически важной задачей для электрической системы.

Нелинейное дифференциальное уравнение, описывающее в некоторой идеали­зации электромеханические переходные процессы простейшей электрической системы, имеет вид:

d² δ/dt² + β dδ/dt + ω_o² sin δ = F(t).

Дата добавления: 2014-01-07