2025
.pdf41
Z C |
Z1 1 j |
cos |
U |
|
j sin |
U |
|
1 j |
k |
|
. |
|
|
|
|
|
|
(3.21) |
||||||
|
|
cos |
U |
j |
|
sin |
U |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Так как Z C |
Z C |
|
Z |
|
, |
г д е Z |
1 |
|
r |
j |
|
X . |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
2 |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Учитывая последние выражения, получаем |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Z C |
Z1 1 |
j |
cos |
U |
|
j |
sin |
U |
1 |
|
j k |
|
|
|
1 |
|
r |
j |
X . |
(3.22) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
cos |
U |
|
j |
sin |
U |
|
|
|
|
2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Анализ выражения (3.22) проводится в следующем порядке: |
|
||||||||||||||||||||||
|
1) задаемся углом нагрузки , причем в случае кругового поля, со- |
|||||||||||||||||||||||
гласно векторной диаграмме, показанной на рис.2.3 угол |
|
1= 0 и |
U = |
; |
||||||||||||||||||||
|
2) находим угол |
из выражения (2.41) с учетом того, что в случае кру- |
||||||||||||||||||||||
гового поля U1=U, в относительных единицах U1= |
; |
|
|
|
|
|
|
|
|
3)рассчитываем сопротивление прямой последовательности Z1 по формуле (2.37);
4)подставляем значения , Z1 и параметры в выражения (3.22) и определяем Z1C ;
5)используя соотношения (2.6), находим емкостное сопротивление конденсатора, соответствующее круговому полю, при данных значениях па-
раметров двигателя и угла нагрузки.
4.СТАТИЧЕСКАЯ УСТОЙЧИВОСТЬ
4.1.Линеаризация системы дифференциальных уравнений
Анализ статической устойчивости синхронных двигателей проводится на основе предположения, что ротор совершает малые и гармонические коле-
бания около установившейся скорости вращения, вызванные возмущением установившегося режима работы. Учитывая данное предположение, опреде-
42
ление границ устойчивой работы конденсаторных СДПМ можно выполнить с помощью рассмотренной выше эквивалентной модели, работа которой опи-
сывается системой нелинейных дифференциальных уравнений первого по-
рядка, при этом систему уравнений необходимо линеаризовать [8]. Линеари-
зация уравнений (2.49), (2.61) и (2.74) заключается в представлении перемен-
ных, входящих в эти уравнения, в виде суммы двух составляющих. Одна со-
ставляющая определяет исходный установившийся режим работы, вторая
учитывает изменение переменных под действием малых возмущений.
Таким образом можно записать
1d= |
1d0+ |
1d ; |
1q= 1q0+ |
1q |
; |
|
1yd= |
1yd0+ |
|
1yd ; |
1yq= 1yq0+ |
1yq ; |
|
i1d=i1d0+ i1d |
; |
i1q=i1q0+ i1q; |
|
(4.1) |
||
i1yd=i1yd0+ i1yd |
; i1yq=i1yq0+ |
i1yq ; |
|
= 0+ ; = 0+ .
В данных выражениях величины с индексом "ноль" характеризуют ус-
тановившийся режим, а величины с индексом " " соответствуют приращени-
ям переменных при малых возмущениях.
Линеаризацию будем проводить, полагая, что двигатель работает от сети бесконечно большой мощности, поэтому величина приложенного на-
пряжения и его частота остается постоянной и не зависит от возмущающих воздействий. Запишем уравнения (2.74) с учетом (4.1)
d |
( |
1d0+ |
1d)=( |
1q0+ |
1q)·( |
0+ |
)-r·(i1d0+ |
i1d)-U1·sin( |
0+ |
) ; |
|||
d |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
d |
( |
1q0+ |
1q)=-( |
1d0+ |
1d)·( |
0+ |
)-r·(i1q0+ |
i1q)-U1·cos( |
0+ |
) ; |
|||
d |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
d |
|
( |
1yd0+ |
1yd)= -ryd·(i1yd0+ |
i1yd) ; |
|
|
|
|
||||
d |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
d |
|
( |
1yq0+ |
1yq)= -ryq·(i1yq0+ |
i1yq) ; |
|
|
|
(4.2) |
||||
d |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
43 |
|
|
|
|
d |
( |
0+ )= |
|
1 |
·(M10+ M1-MC0-M20-km·( + - )) ; |
|
|
|
d |
|
|
|
||||
|
|
|
H |
|
|
|||
|
d |
( |
0+ )= -( 0+ ) . |
|
|
|||
|
d |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Подставив (4.1) в уравнение (2.61), получим выражения, определяю- |
|||||
щие приращение электромагнитного момента |
|
|||||||
|
|
|
M1= |
|
1d·iq0+ 1d0· i1q- 1q0· i1d- |
1q·i1d0 ; |
(4.3) |
|
и значение электромагнитного момента в установившемся режиме |
|
|||||||
|
|
|
M10= |
|
1d0·iq0- 1q0·i1d0 . |
|
(4.4) |
|
|
|
|
Уравнения потокосцеплений с учетом (4.1) запишутся в следующем |
|||||
виде |
|
|
|
|
|
|
||
|
|
|
1d0+ |
|
|
1d=(Xad+X )·(i1d0+ i1d)+Xad·(i1yd0+ i1yd) + ; |
|
|
|
|
|
1q0+ |
|
|
1q=(Xaq+X )·(i1q0+ i1q)+Xaq·(i1yq0+ i1yq) ; |
|
|
|
|
|
1yd0+ |
|
1yd=Xad·(i1d0+ i1d)+(Xad+X |
yd)· i1yd + ; |
(4.5) |
|
|
|
|
1yq0+ |
|
1yq=Xaq·(i1q0+ i1q)+(Xaq+X yq)· i1yq . |
|
||
|
|
|
Рассматривая линеаризацию уравнений конденсаторного СДПМ в слу- |
чае малых возмущений, полагаем, что частота приложенного напряжения в
каждый момент времени постоянна ( =const), возмущающее воздействие бес-
конечно мало, поэтому угол |
, на который отклоняется ротор электродвига- |
|||||||||
теля, также мал, и можно считать, что cos |
= 1, |
sin |
= |
|
. Учитывая |
|||||
это, в уравнениях (4.2) можно записать, что |
|
|
|
|
|
|||||
sin( |
0+ |
)=sin |
0·cos |
+sin |
·cos |
0=sin |
0+ |
·cos |
0 ; |
|
cos( |
0+ |
)=cos |
0·cos |
-sin |
·sin |
0=cos |
0- |
·sin |
0 . |
(4.6) |
Преобразуем уравнения (4.2) и (4.5), исключая при этом члены, соот-
ветствующие режиму, пренебрегая произведением приращений, имеющие второй порядок малости и учитывая выражения (4.6), запишем
d |
1d= · |
1q+ 1q0 · -r· i1d-U1 ·cos 0· |
; |
|
|
||||
d |
||||
|
|
|
|
|
44 |
|
d |
1q=- · |
1d- 1d0 · -r· i1q-U1 ·sin 0 · ; |
|
|
|||
d |
|||
|
|
d |
|
|
1yd= -ryd · i1yd ; |
|||
|
|
|
||||
d |
||||||
|
|
|
|
|
||
d |
|
|
1yq= -ryq · i1yq ; |
|||
|
|
|
||||
d |
|
|||||
|
|
|
|
|
||
d |
|
|
= |
1 |
·( M1-km· ) ; |
|
d |
|
|
|
|||
|
|
|
H |
|||
d |
|
=- . |
||||
|
|
|||||
d |
|
|||||
|
|
|
|
|
||
|
|
|
А уравнения (4.5) примут вид |
|||
|
|
|
|
|
1d=(Xad+X )· i1d+Xad i1yd ; |
|
|
|
|
|
|
1q=(Xaq+X )· i1q+Xaq i1yq ; |
|
|
|
|
|
|
1yd=Xad· i1d+(Xad+X yd) i1yd ; |
|
|
|
|
|
|
1yq=Xaq· i1q+(Xaq+X yq) i1yq . |
(4.7)
(4.8)
Решая совместно уравнения (4.8), найдем выражения для токов i1d,
i1q, i1yd, i1yq
i1d |
|
1d |
X yd |
|
1yd |
X ad |
|
; |
|
|
|
|
|
|
X ad |
X |
X yd |
X ad2 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
i1q |
|
1q |
X yq |
|
1yq |
X aq |
|
; |
|
|
|
|
|
|
X aq |
X |
X yq |
X aq2 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
i1yd |
|
1yd |
X ad |
X |
|
|
1d |
X ad |
; |
(4.9) |
|||
|
X ad |
X |
X yd |
X ad2 |
|
||||||||
|
|
|
|
|
|
||||||||
i1yq |
|
1yq |
X aq |
X |
|
|
1q |
X aq |
|
. |
|
||
|
X aq |
X |
X yq |
X aq2 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
45 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Запишем (4.7) с учетом (4.9) |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
d |
|
|
1d= |
· |
|
|
1q+ |
1q0· |
|
-r· |
|
|
1yd |
X yd |
1yd |
X ad |
-U1·cos |
|||||||||||
|
d |
|
|
|
|
|
|
X ad |
X |
X yd |
|
|
X ad2 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
d |
|
|
1q=- |
· |
|
|
1d- |
1d0· |
|
-r· |
|
|
1yq |
X yq |
1yq |
X aq |
|
-U1·sin |
||||||||||
|
d |
|
|
|
|
|
|
X aq |
X |
X yq |
|
|
X aq2 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
d |
|
|
1yd= -ryd· |
|
1yd |
X ad |
|
X |
|
1d |
X ad |
|
; |
|
|
|
|
|||||||||||
|
d |
|
|
|
X ad |
X |
|
X yd |
X ad2 |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
d |
|
|
1yq= -ryq· |
|
|
1yq |
X aq |
|
X |
|
1q |
X aq |
; |
|
|
|
|
|
||||||||||
|
d |
|
|
|
X aq |
X |
|
X yq |
X aq2 |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
d |
|
|
= |
|
1 |
|
·( |
|
|
1d·i1q0+ |
|
|
|
|
|
1q |
X yq |
1q |
X aq |
- |
|
|
||||||
|
d |
|
|
|
H |
|
|
1d0 |
|
|
X aq |
X X yq |
X aq2 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
- |
|
|
|
|
|
1d |
X yd |
1d |
|
X ad |
|
- |
1q·i1d0 - km· |
) ; |
|
|
|
||||||||||||
1q0 |
|
|
X ad |
X |
X yd |
|
X ad2 |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
=- . |
|
|
||
d |
||
|
0· ;
0· ;
(4.10)
Система уравнений (4.10) представляет собой линейную систему диф-
ференциальных уравнений первого порядка с постоянными коэффициентами.
Данные уравнения описывают режим малых гармонических колебаний рото-
ра конденсаторного СДПМ в полном объеме, и могут быть использованы для исследования статической устойчивости .
Определим величины, входящие в данную систему уравнений и соот-
ветствующие установившемуся режиму работы, для этого приравняем нулю производные в выражениях (2.49), (2.74). После промежуточных преобразо-
ваний получим
|
U1 cos 0 |
Xaq |
X |
U1 sin 0 |
r |
|
i1d0 |
|
|
|
|
|
; |
r 2 |
Xad X |
Xaq X |
2 |
|
||
|
|
|
46
|
U1 cos |
0 |
r |
U1 sin |
0 Xad |
|
X |
(4.11) |
|
i1q0 |
|
|
|
|
|
|
|
. |
|
r 2 |
Xad |
X |
Xaq |
X |
2 |
|
|||
|
|
|
|
||||||
|
1d0=(Xad+X )·i1d0+ ; |
|
|
|
|
|
|
||
|
1q0=(Xaq+X )·i1q0 . |
|
|
|
|
|
|
(4.12) |
|
Напряжение прямой последовательности U1 определяется по методике, |
|||||||||
изложенной в пункте 3.2. Угол нагрузки |
0, соответствующий установивше- |
||||||||
муся режиму работы, является заданной величиной. |
|
|
|
|
Электромагнитный момент от напряжений прямой последовательно-
сти в установившемся режиме работы определяется по уравнению (4.4).
Суммарный электромагнитный момент установившегося режима рабо-
ты модели конденсаторного СДПМ, учитывая (4.4), (2,4), а также генератор-
ный тормозной момент, создаваемый постоянным магнитом, равен |
|
M0=M10-M20-M30 , |
(4.13) |
где M20- момент установившегося режима, соответствующий обратной последовательности, определяется по методике, изложенной в [8];
M30 - генераторный тормозной момент от постоянного магнита рассчи-
тывается по формулам (2.66), (2.67).
Конденсаторные СДПМ являются более общим типом синхронного двигателя малой мощности. Поэтому полученные уравнения СДПМ могут быть использованы для анализа статической устойчивости синхронных реак-
тивных двигателей (СРД) при работе их от однофазной сети . Для исследова-
ния конденсаторных СРД необходимо в уравнениях (2.32), (2.33), (2.34), (2,41), (4.11), (4.12) приравнять нулю возбужденность постоянных магнитов ,
а в формуле (2.70) коэффициент km3. В однофазных СРД кроме выше пере-
численных допущений, необходимо принять равной нулю емкость конденса-
тора в выражении (3.1).
4.2. Методика и алгоритм расчета границ статической устойчивости
47
Решение вопроса устойчивости, на основе линеаризованных уравне-
ний доказано А.М.Ляпуновым .
Устойчивость системы линейных дифференциалных уравнений опре-
деляется по корням характеристического уравнения, составленного для дан-
ной системы [8]. Характеристическое уравнение для системы (4.10) будет иметь вид
a0·p6+a1·p5+a2·p4+a3·p3+a4·p2+a5·p+a6=0 . (4.14)
Система уравнений (4.10) будет устойчивой, а следовательно устойчи-
вым будет и режим работы данного электродвигателя, если все вещественные корни и вещественные части комплексных корней характеристического урав-
нения (4.14) будут отрицательными. Нарушение устойчивости происходит,
когда становится положительным хотя бы один из вещественных корней или вещественная часть комплексного корня.
Определение корней характеристического уравнения вызывает опре-
деленные затруднения, связанные с большой трудоемкостью. Провести ана-
лиз устойчивости, не прибегая к вычислению корней характеристического уравнения, позволяют специальные критерии [10]. Наиболее удобным и про-
стым в применении при практическом решении задач устойчивости на вы-
числительной технике является алгебраический критерий Рауса. Границы ус-
тойчивости конденсаторных СДПМ с помощью данного критерия определя-
ются следующим образом. Составляется таблица Рауса, коэффициенты в ней вычисляются через коэффициенты характеристического уравнения (4.14) и
проводится анализ коэффициентов первого столбца таблицы [8].
Таблица 4.1.
a0 |
a2 |
a4 |
a6 |
|
|
|
|
a1 |
a3 |
a5 |
0 |
|
|
|
|
C13=a2-a0/a1·a3 |
C23=a4-a0/a1·a5 |
a6 |
0 |
|
|
|
|
48
C14=a3-a1/C13·C23 |
C24=a5-a1/a13·a6 |
0 |
0 |
|
|
|
|
C15=C23-C13/C14·C24 |
a6 |
0 |
0 |
|
|
|
|
C16=C24-C14/C15·C16 |
0 |
0 |
0 |
|
|
|
|
a6 |
0 |
0 |
0 |
|
|
|
|
Работа конденсаторного СДПМ при данном сочетании параметров бу-
дет устойчивой, если все коэффициенты первого столбца таблицы Рауса бу-
дут положительными. В случае, когда хотя бы один из коэффициентов стано-
вится отрицательным, нарушается устойчивая работа конденсаторного СДПМ.
Практическое исследование статической устойчивости предполагает определение пределов изменения параметров двигателя, сети и нагрузки, при которых работа синхронной машины будет устойчивой. Поэтому границы ус-
тойчивости удобно находить в плоскости двух параметров, от которых зави-
сят коэффициенты характеристического уравнения [8].
Построение найденных областей статической устойчивости будем проводить в прямоугольной системе координат. Одна ось показывает измене-
ние параметра , характеризующего исходный установившийся режим работы,
а вторая варьируемый параметр. В качестве параметра установившегося ре-
жима работы наиболее удобно использовать результирующий электромагнит-
ный момент M0, так как по его знаку можно определить режим работы элек-
трической машины. Варьируемыми параметрами выберем активное сопро-
тивление обмотки статора r, а в случае работы двигателя от регулируемого источника - относителное значение частоты . Следовательно, будем исполь-
зовать системы координат M0(r) и M0 ( ).
Как уже отмечалось, нарушение статической устойчивости может при-
вести к двум различным режимам работы СДПМ. В трехфазных СДПМ апе-
риодическое нарушение устойчивости (сползание) наступает, когда послед-
49
ний коэффициент первого столбца таблицы Рауса принимает отрицательные значения. В этом режиме скорость вращения ротора двигателя становится меньше синхронной. Двигатель выпадает из синхронизма. Таким образом,
сползание определяет предел статической перегружаемости, соответствую-
щий максимуму электромагнитного момента.
В конденсаторном СДПМ анализируемые уравнения записываются для двигателя прямой последовательности эквивалентной модели, и отрицатель-
ный знак коэффициента a6 соответствует максимальному значению электро-
магнитного момента только прямой последовательности M10. В этом случае сползание необходимо определять по суммарному максимальному электро-
магнитному моменту M0, учитывающему действие момента от поля обратной последовательности M20 и от постоянного магнита M30.
Второй режим работы, возникающий в результате нарушения статиче-
ской устойчивости в трехфазных СДПМ, характеризуется появлением само-
возбуждающихся периодических колебаний скорости вращения ротора (са-
мораскачивание). Возникновение самораскачивания характеризует отрица-
тельный знак предпоследнего коэффициента (С16) первого столбца таблицы Рауса.
Для конденсаторных СДПМ граница самораскачивания определяется,
так же, как и в трехфазных двигателях, по знаку коэффициента С16, таблицы Рауса.
Методика расчета границ областей статической устойчивости конден-
саторных СДПМ реализована в виде программы для ЭВМ типа IВМ . Про-
грамма написана на языке ФОРТРАН и приведена в приложении 4. Алго-
ритм работы данной программы можно рассмотреть на основе структурной схемы, представленной на рис.4.1.
Как видно из представленной структурной схемы, программа работает в следующем порядке. В начале задаются исходные данные для расчета,
включающие в себя параметры математической модели r, .Xad, Xaq, X , ryd, ryq,
50
X yd, X yq, H, km1, C, . Задаются также подводимое к исследуемому двигателю напряжение или закон изменения = f( ) . Вводятся начальные значения варьируемых параметров, интервал и шаг их изменения, а также количество циклов расчета. Ввод перечисленных данных осуществляется через текстовый файл PARAM.DAT с помощью любого текстового редактора (блок 1).
В блоке 2 производится расчет напряжений прямой и обратной после-
довательности и угла нагрузки с помощью подпрограммы U1U2N.
Далее в блоке 3 по подпрограмме AMOB(S) вычисляется величина асинхронного тормозного момента M20 и коэффициента электромагнитного демпфирования km2.
Затем в блоке 4, используя подпрограмму АМ3(W0), определяется мо-
мент, обусловленный наличием постоянных магнитов M30 и коэффициент демпфирования от этого момента km3.
На следующем этапе расчета в блоке 5 вычисляются коэффициенты матрицы для системы уравнений.
В блоке 6 производится расчет коэффициентов характеристического уравнения методом Леверье по подпрограмме LEVER.
Затем в блоке 7 вычисляются коэффициенты таблицы Рауса по под-
программе RAUS.
Пуск
1. Ввод исходных D данных из файла
sdpm.dat