
- •1 Фізичне тлумачення перехідних процесів та їх вплив на
- •Вихідні дані та припущення для запису рівнянь
- •Рівняння за другим законом Кірхгофа
- •Зв’язок параметрів еквівалентної заступної схеми асинхронного двигуна з коефіцієнтами диференційних рівнянь:
- •Вивод диференційних рівнянь машини для нових змінних
- •Можливі варіації системи диференційних рівнянь на практиці.
- •Комплексна система рівнянь узагальненої машини
- •Система диференційних рівнянь, записана через струми
- •Моделювання синхронних двигунів Система диференційних рівнянь реальної машини
- •Рівняння потокозчеплень та індуктивностей обмоток
- •Індуктивність фазних обмоток
- •Взаємоіндуктивності
- •Взаємоіндуктивності з обмотками ротора
- •Коефіцієнти взаємоіндукції статорних обмоток зі статорними
Зв’язок параметрів еквівалентної заступної схеми асинхронного двигуна з коефіцієнтами диференційних рівнянь:
L1,L2,M1,M2,M12
Для аналітичного розрахунку перехідного процесу треба знати числові величини всіх коефіцієнтів диференційного рівняння, індуктивностей та взаємоіндуктивностей. В довідниках по АД можна знайти параметри еквівалентної заступної схеми (Т-подібна або Г-подібна уточнена).
Звичайні параметри Т-подібної схеми є: x1,x2,R1,R2,xm.
x2 – зведений параметр роторної обмотки
R1,R2 – активні опори (R2-зведений) фаз статора і ротора
Для встановлення зв’язку між параметрами Т-схеми та коефіцієнтами систем слід записати диференційне рівняння в усталеному режимі, та порівняти їх з аналітичними виразами, які звичайно використовуються в загальній теорії електричних машин.
Для складання такої системи запишемо функціональні залежності в часі струмів статора і ротора за умови довільного значення початкової фази.
де: Ім1(Ім2) – амплітуди струмів статора (ротора)
() – початкові фази струмів статора (ротора)
Усталене значення напруги статора (без ротора):
UA=UM1sin(0t+)
UB=UM1sin(0t+-120o)
UC=UM1sin(0t++120o)
Роторні напруги дорівнюють нулю.
А=-
Рис.
В зв’язку з тим, що при перетворенні координат реальної машини в модельну, миттєві значення змінних знаходяться як проекції зображуючих векторів на відповідні координатні осі. Кожна проекція є косинусна функція прилеглого кута між зображуючим вектором і відповідною віссю. Тому будемо вважати, що усталені значення струмів та напруг АД теж записані через косинус , при цьому будемо вважати, що їх початкові фази такі самі, як і для синусних величин. Тому підставимо значення струмів статора і ротора в рівняння потокозчеплень, враховуючи функції cos та =t . Будемо мати:
А=Im1 [L1cos(0t+)+М1cos(0t+-120o)+М1cos(0t++120o)]+
+Im2[M12 costcos[(0 -)t+]+M12 cos(t+120o)cos[(0 -)t+
+-120o]+M12 cos(t-120o)cos[(0 -)t++120o]]
Виконаємо перетворення гармонічних функцій, що входять в рівняння згідно з формулами тригонометрії:
Cos(+)=coscos-sinsin
Cos(-)=coscos+sinsin
Для статора:
=(t+), =120o
Для ротора:
=[(0-)t+]
Перетворимо вираз для струму Im1:
L1cos(0t+)+М1cos(0t+-120o)+М1cos(0t++120o)= L1cos+
+М1cos(-)+М1cos(+)=L1cos(0t+)+М1[cos(0t+)cos120o+
+sin(0t+)sin120o]+М1[cos(0t+)cos120o-sin(0t+)sin120o]=
=L1cos(0t+)+М1[cos(0t+)(-0,5)+sin(0t+)
]+М1[cos(0t+
+)(-0,5)-sin(0t+) ]= (L1-М1)cos(0t+)
Перетворимо вираз біля струму Im2:
M12 costcos[(0 -)t+]+M12 cos(t+120o)cos[(0 -)t+-120o]+
+M12 cos(t-120o)cos[(0 -)t++120o]=M12 coscos+
+M12 cos(+)cos( -)+M12 cos(-)cos( +)=M12 coscos+
+M12 [coscos120o-sinsin120o][coscos-sinsin]+
+M12
[coscos-sinsin][coscos-sinsin]=
M12
cos(0
t+)
Для потокозчеплення А маємо вираз:
А=Im1[(L1-M1)cos(0 t+)]+Im2[ M12cos(0 t+)]
Виконаємо аналогічне перетворення для ротора, в результаті будемо мати:
а=Im2[(L1-M1)cos[(0 -)t+]]+Im1[ M12cos[(0-) t+]]
Перейдемо до комплексного виду запису рівняння електричної рівноваги. Для виразимо миттєві значення через комплексні амплітуди та оператори обертання зображуючого вектора.
UA=Ume(t+)
iA=Im1e(t+)
ia=Im2e[(-)t+)
Підставимо миттєві значення в рівняння електричної рівноваги фази. Будемо мати:
UA=Іm1e(t+)RA+
Ume(t+)=Іm1e(t+)RA+
Запишемо потокозчеплення фази А в комплексному вигляді:
А=(L1-M1)Im1e(t+)+ Im2M12e(t+)
Знайдемо похідну
=j0(L1-M1)Im1e(t+)+j0
Im2M12e(t+)
Введемо комплексні амплітуди для струмів статора і ротора:
Im1=Im1ej
Im2=Im2ej
Підставимо комплексні амплітуди струмів:
=j0(L1-M1)Im1et+j0 Im2M12et
Підставимо значення похідної в рівняння електричної рівноваги в комплексному вигляді, будемо мати:
=
RA+[j0(L1-M1)Im1et+j0
Im2M12et]
Перетворимо одержаний вираз таким чином:
=
RA+j0(L1-M1)Im1+j0
Im2M12+j0
m1M12-j0
Im1M12
Згрупуємо подібні члени:
= RA+j[0(L1-M1)-0 M12 ]Im1+j0 M12( m1+ m2)
Введемо позначення:
x1=0(L1-M1)-0 M12
xm= M120
Підставляємо
введені параметри і ділимо праву і ліву
частини на
:
=
1(R1+jx1)+jxm(
1+
2)
В загальній теорії електричних машин відомо, що струм холостого хода :
( 1+ 2)= 0
Тоді маємо:
= 1(R1+jx1)+jxm 0
Для ілюстрації зв’язку між параметрами еквівалентної заступної схеми :
=-
1+j
1x1+
1R1
=j0[(L1-M1) 1- M12 1]- 1R1- j0 M12 10
Розглянемо еквівалентну Т-подібну заступну схему :
Рис.
Rнав=R2
(R2
+Rнав)=
xm=j0 M12
x1=xS-xm
xS=x1+xm
З виразу для х1 видно, що x1=xS-xm або xS=x1+xm
xS – синхронна реактивність фази статора коли ротор рухається зі швидкістю поля.
Для ротора виконаємо аналогічні перетворення як і для статора. Тобто запишемо рівняння електричної рівноваги в комплексному вигляді для усталеного режиму.
Uа=іаR2+
=0
Вважаємо, що:
І2=
ej
Це рівняння буде мати вигляд:
0=j(0 -)(L1-M1) m2ej +j(0 -)( M12 m1ej )+
+ m2R2ej )
Введемо комплекси діючих значень:
0=j(0 -)(L1-M1) 2+j(0 -)( M12 1)+ 2R2)
Для того, щоб позбавитись від струму 1 в рівнянні добавимо і віднімемо вираз: j(0 -) M12 2
0=j(0 -)(L1-M1) 2+j(0 -) M12 2+R2I2+ j(0 -) M12 2-
-j(0 -) M12 2
Остаточно маємо:
j(0 -)(L1-M1) 2- j(0 -) M12 2+j(0 -) M12 10+R2I2
З Т-подібної схеми видно, що
Замінимо
–Е2
на -
=jI10xm
, та будемо мати:
-jI10xm
jI10xm
Для порівняння виразу з еквівалентної схеми з формулою через коефіцієнти L2, M2, M12 введемо ковзання в це рівняння. Для цього помножимо кожну складову на 0 і поділимо на 0-. В результаті будемо мати:
j0(L2-M2)
2-
j0
M12
2+j0
M12
10
З порівняння одержаного виразу з формулою по еквівалентній заступній схемі видно, що:
jx2=j0(L2-M2) - j0 M12
xm=j0 M12
Тепер ясно, що x2=xr-xm, або xr=x2+xm
де:
xr – синхронна реактивність фаз обмотки ротора при живленні її напругою частоти 0 коли ротор розімкнутий
Перетворення диференційних рівнянь машин.
Фізичний смисл перетворень складається в тім, що наявність періодичних коефіцієнтів рівнянь вимушує дослідників шукати шляхів спрощення рівнянь з метою запису їх з постійними коефіцієнтами. Для цього треба здійснити заміну змінних в рівнянні. Будемо вважати, що нова система диференційних рівнянь записана відносно нових змінних відображає ідеалізовану синхронну машину, у якої напруги, струми і потокозчеплення всіх фаз мають аналітичний зв’язок з реальними спеціальними формулами заміни змінних. Оскільки система рівнянь тепер не має періодичних коефіцієнтів, то фізично це означає, що ротор і статор цієї ідеалізованої машини є взаєморухомі тому, що змінні коефіцієнти в системі з’явились як результат руху роторних обмоток відносно статорних. Будемо вважати, що статор і ротор ідеалізованої машини обертаються в просторі з довільною частотою к . При цьому конкретний вибір цієї частоти обертання може змінювати зовнішній вигляд рівнянь зпрощюючі, або ускладнюючі їх.
Припустимо також для спрощення системи, що ідеалізована машина є двофазна реальної трьохфазної машини за величиною ЕРС зазора створеними струмами статора і ротора. Схематично така двофазна машина має вигляд:
Рис.
ЕРС фази, яка враховує наявність обертання ротора відносно статора в реальній машині, а також враховується при цьому задана швидкість обертання к координатних осей UV.
Напруги відповідних фаз мають індексацію відповідної осі обмотки. Всі електричні величини також прив’язуються до назви координатних осей. Цифровими індексами 1(2) позначаються статорні (роторні) величини.
Для перетворення системи диференційних рівнянь трьохфазної реальної машини вирішується 2 задачі:
1) необхідно знайти формули зміни старих змінних на нові,
2) використовуючи ці формули одержати систему диференційних рівнянь в нових змінних.
Для вирішення цих задач розглянемо рис.2, на якому показані координатні системи реальних фаз статора і ротора та координатна система ідеалізованої машини. Крім того, будемо вважати, що зображуючі вектори ЕРС статора і ротора представлені на цьому рисунку в вигляді струму ( вектора еквівалентного відповідній ЕРС ).
Рис.2
Як відомо з загальної теорії електричних машин миттєві значення струмів в реальних фазах машини визначаються як проекції зображуючого вектора струму статора (ротора) на координатні осі реальних фаз статора (А,В,С), ротора (а,в,с). Аналогічний висновок буде і для координатної системи UV.Струм iu1 буде проекція струму Іі на координатну вісь U. А струм іv1 є проекція цього ж вектора на вісь V.
Для ротора миттєві значення теж визначаються як проекції зображуючого вектора І2 на координатні осі для відповідних фаз.
Якщо тепер виразимо аналітично миттєві значення струмів фаз як проекції зображуючого вектора, то ми можемо знайти необхідні формули зв’язку між старими змінними і новими змінними. Для того, щоб можна було записати аналітичні вирази для проекції струмів Іі та І2 треба позначити величини кутів між зображуючим вектором та осями нерухомих та рухомих фаз. За точку відліку кутів приймаємо напрямок фази А статора.
Запишемо миттєві значення струмів фаз через відповідні проекції:
іА=І1cos1
і В=І1cos(1-120)
і С=І1cos(1+120)
де 1- кут між вектором І1 та фазою А
І1-
модуль вектора
Аналогічно знайдемо вирази миттєвих значень струмів фаз еквівалентної машини, які носять назву U,V можна записати також як проекції вектора струму на координатні осі U і V.
іu1= cos(кt-1)
іV1=-
sin(кt-1)
Тепер задача складається в тім, щоб із цих двох систем знайти вирази для струмів іu1 та іV1 в залежності від реальних струмів іА, іВ, іС, тобто , знайти формули перетворення реальної машини в еквівалентну двофазну. Розпишемо в першому рівнянні другої системи:
іu1= cos(кt-1)= cos1 cosкt+ sin1 sinкt
З
цього запису видно, що
cos1=
іА
,
а значення
sin1=
,останній
вираз можна довести наступним чином: з
ТОЕ відомо, що іА+
і
В+
і
С=0
і В- і С= І1cos(1-120)-І1cos(1+120)
Якщо розкладемо cos різниці та суми двох кутів, то вважаючи на те, що іА=-( і В+і С)= І1cos1, то після підстановки cos 120та sin 120 будемо мати вираз:
І1cos1=
- І1(
cos1-
sin1+
cos1+
sin1)
Після перетворень знайдемо, що:
sin1=
Після підстановки значень в перше рівняння будемо мати:
іu1= іА сos кt + sin кt
Останній вираз може бути перетворений :
іu1=
іА
сos кt
+ і
Вcos(кt
-120)+і
Сcos(кt
+120)
Таким чином, для того, щоб записати струм фази U, вісь якої обертається в просторі з довільною швидкістю к , треба знати струми реальних фаз та конкретне значення швидкості к.
Виконаємо аналогічні перетворення для знаходження струму фази U статора (іV1), який має такий вигляд:
іV1= - іА sin кt + і Вsin(кt -120)+і Сsin(кt +120)
Таким чином, ми одержали формули прямого перетворення від реальних струмів машини до еквівалентних струмів двофазної ідеалізованої машини. Якщо машина буде m-фазна перетворюватись в двофазну, то коефіцієнт буде не 2/3, а 2/ m.
Розглянемо
аналогічні перетворення реальних
струмів ротора іа
, ів
, іс
в струми іu2
та іv2.
Як видно з рисунку координатних систем
лінії а, в, с позначають осі фазних
обмоток ротора. Вектор струму
є зображуючий вектор сумарної ЕРС
ротора. Згідно з загальною теорією,
вектор струму
повинен бути в протифазі відносно струму
так, щоб їх сума завжди дорівнювала
намагнічуючому струмові:
+
=
const.
Зв’язок між струмом та реальними струмами фаз запишемо як проекції струму і2 на реальні фазні осі ротора.
В цих виразах:
2 – кут між вектором віссю фази А статора
- кут, миттєве значення положення ротора відносно статора =t.. Відраховується від осі а ротора в формулах перетворення.
Аналогічно як для статора виразимо струми двофазних обмоток ротора еквівалентної машини як проекції вектора на координатні осі U і V.
Тепер з останніх двох систем можна записати залежності іu2(іа,ів,іс) , іv2(іа,ів,іс). Для спрощення процедур приймаємо 2 - =.
Розписуються аналогічно косинуси різниці та суми двох кутів (- 120) та (- 120). Після всіх перетворень, таких же, як і для статора, струми фаз ротора еквівалентної машини мають вигляд:
іu2= іа сos (кt -)+ і вcos( (кt-) -120)+і сcos( (кt-) +120)
cos (кt-t)= cos (к-) t
іV2= - іа sin (кt-) + і вsin((кt-) -120)+і сsin((кt-) +120)
Вирази для струмів статора, а також для струмів ротора в осях U і V, являють собою формули прямого перетворення, тобто формули перехода від реальних струмів двофазної машини.
Виникає питання чи досить двох нових змінних іu1 , іV1 для повної заміни трьох реальних струмів іА , іВ , іС, при переході до еквівалентної двофазної машини. Відповідь на поставлене питання можна одержати шляхом зворотнього перетворення, тобто шляхом визначення струмів фаз А, В, С реальної машини через струми еквівалентної машини. Для цього помножимо струм іu1
іu1cosкt= cosкtіА сos кt + і Вcos(кt -120)+і Сcos(кt +120)
Струм іV1помножимо:
іV1(- sinкt)= sin кt іА sin кt + і Вsin(кt -120)+і Сsin(кt +120)
Виконаємо відомі перетворення:
іА сos2 кt + і Вcosкt (cosкt cos 120+sinкt sin 120)+ і Сcosкt(cosкt cos 120-sinкt sin 120)
Якщо з 1 рівняння визначимо струм фази А за допомогою його перетворень, то будемо мати:
іА=
іu1cosкt-
іu1
sin кt+
(іА
+( іВ
+ іС))
Введемо іV1= (іА +( іВ + іС)) , де цей струм є струм нульового слідування фаз при розкладі несиметричної системи реальних струмів на симетричні складові: пряму, зворотню та нульову.
З останніх виразів видно, що струм і01 не може бути виражений через іu1 та іV1. Це означає, що ні і01, ні струми фаз А,В,С не можуть в загальному випадку бути записані через нові змінні, т.ч. для трьохфазної симетричної системи з нульовим проводом треба враховувати і01, а для системи без нульового провода струм і01=0, але в загальному випадку формули перетворення статорних величин мають 3 рівняння, тобто:
іu1= іА сos кt + і Вcos (кt -120)+і Сcos(кt+120)
іV1= - іА sin кt + і Вsin(кt -120)+і Сsin(кt +120)
і01== (іА + іВ + іС)
Якщо виконати всі аналогічні процедури для ротора, то одержимо систему з трьох рівнянь для формул прямого та зворотного перетворень в загальному випадку( несиметрія ротора) будемо мати наступну систему:
1) іu2= іа сos (кt -)+ і вcos( (кt-) -120)+і сcos( (кt-) +120)
2) іV1= - іа sin (кt-) + і вsin((кt-) -120)+і сsin((кt-) +120)
і02== (іА + іВ + іС)
Висновки:
Фізичний смисл формул перетворення для АМ в тім, що реальна трьохфазна машина зводиться до еквівалентної двофазної машини за величиною МРС, створеної струмами статора і ротора.
При цьому перетворенні завжди статор і ротор взаємно нерухомі і обертаються в просторі з довільною швидкістю к.
Величина швидкості ротора відносно статора реальної машини та величина швидкості к враховується додатковими ЕРС в фазах статора та ротора (еu1, еV1, еu2, еV2).
Оскільки любу трьохфазну систему напруг, струмів, ЕРС, або потокощеплень можна представити відповідними зображуючими векторами, то кожна з цих величин може бути перетворена з трьохфазної системи в двофазну згідно з формулами узагальненого перетворення. Ці формули мають вигляд:
fu1(t)= fA(t) сos кt + fB(t)cos(кt-120)+ fC(t)cos(кt+120)
fV1(t)= - fA(t)sin кt + fB(t) sin (кt-120)+ fC(t)sin (кt+120)
fu2(t)= fа(t) сos кt + fв(t)cos(кt-120)+ fс(t)cos(кt+120)
fV2(t)= - fа(t)sin кt + fв(t)sin (кt-120)+ fс(t)sin (кt+120)
Перетворення диференційних рівнянь асинхронної машини
Запишемо диференційне рівняння для потокозчеплень статора та ротора в нових координатах, тобто в системі осей, що обертаються з довільною швидкістю к. Для цього А множимо на сos кt, рівняння для В множимо на cos(кt-120), та рівняння для С множимо на cos(кt+120)
А сos кt= L1іA сos кt + M1i Вcosкt+M1і Сcosкt+ M12cosкt cos + +M12i вcos кt cos(+120o) + M12i сcos кt cos(-120o)
В сos(кt-120o)= М1іA сos(кt-120o) + L 1i Вcos(кt-120o)+M1і Сcos(кt--120o)+M12cos(кt-120o)cos(-120o)+M12i вcos(кt- 120o)cos+ +M12i сcos (кt-120o)cos(+120o)
Ссos(кt+120o)= М1іAсos(кt+120o)+М1iВcos(кt+120o)+L1іСcos(кt+ +120o)+M12іаcos(кt+120o)cos(-120o)+M12iвcos(кt+120o)cos(-120о)+M12i сcos (кt+120o)cos
Якщо окремо складемо ліві частини рівнянь, то одержимо повне потокозчеплення фази U статора, а саме:
U1= А сos кt+В сos(кt-120o)+ Ссos(кt+120o)
Праві частини будемо складати відповідно до членів з струмами фаз А,В,С статора та відповідними коефіцієнтами при цих струмах, які повинні бути однаковими. Записуємо окремо:
L1іA сos кt+L 1i Вcos(кt-120o)+ L1іСcos(кt+ +120o) = L1 іu1
Виконаємо подібну операцію і для коефіцієнта М1:
М1(iВ+іС) сos кt+(iА+іС) cos(кt-120o) +(iА+іВ) cos(кt+120o) =- М1 іu1
iА +iВ+іС=0
iВ+іС=- iА
Виконаємо відповідну процедуру для коефіцієнта М12:
М12 іаcosкtcos+i вcos кt cos(+120o)+i сcos кt cos(-120o)+
+іаcos(кt-120o)cos(-120o)+iвcos(кt-120o)cos+iсcos(кt-120o)
cos(+120o)+іаcos(кt+120o)cos(+120o)+iвcos(кt+120o)cos(-120о)+
+i сcos (кt+120o)cos = М12 іu2
Згідно з формулами узагальненого перетворення та розкладання косинуса суми та різниці двох кутів та групування подібних членів одержимо М12 іu2
U1= L1 іu1- М1 іu1+ М12 іu2=( L1- М1) іu1+ М12 іu2
Якщо виконати аналогічні перетворення по осі V статора, перемноживши потокозчеплення реальних фаз статора А, В, С
V1= Аsin кt +В sin (кt-120)+ С sin (кt+120)
V1==( L1- М1) і V 1+ М12 іV2
Перетворення роторних рівнянь
Якщо реальне потокозчеплення фаз ротора а, в, с помножимо відповідно на сosкt, сos(кt-120), сos(кt+120) і виконаємо аналогічні як для статора перетворення, то запишемо вирази для потокозчеплень фаз ротора еквівалентної двофазної машини по осям U та V.
U2=( L2- М2) іu2+ М12 іu1
V2=( L2- М2) і V 2+ М12 іV1
Якщо перемножимо праву і ліву частини кожного рівняння для потокозчеплення на 0, то будемо мати:
0 U1=xs іu1 +xm іu2
0 V1= xs іV1 +xm іV2 (а)
0 U2= xp іu2+xm іu1
0 V2= xp і V 2 +xm іV1
xs (xp) – синхронні реактивності статора (ротора) – знаходяться з еквівалентної заступної схеми.
xm – реактивність взаємоіндукції між фазою статора і ротора - знаходяться з еквівалентної заступної схеми.