- •Редактор Козвонина Е. Г.
- •Текст напечатан с оригинал-макета, предоставленного составителем.
- •1. Программа курса
- •2. Библиографический список
- •4. Задание
- •Рис.1. Схема включения ДПТ
- •Таблица 1
- •5.Методические указания к выполнению контрольной работы
- •Перепишем данную систему относительно производных:
9
5.Методические указания к выполнению контрольной работы
Решить задачу с данными двигателя: Uн=220 В; Рн=6 кВт; Iн=33,5 А;
ωн=78,54 1/с; Rя=0,73 Ом; Lя=26,2 мГн; J=0,162 кг·м2. Построить переходный процесс при iя0=3,35 А; ω0=7,854 1/с.
5.1. Система дифференциальных уравнений
В качестве системы дифференциальных уравнений взяты уравнение равновесия напряжений, составленное для якорной цепи схемы рисунка1, а также уравнение равновесия моментов для вращательного движения ( уравнение движения привода), которые имеют вид
u я |
= R я ×i я |
+ L я |
× |
di я |
+ k ×Ф ×w, |
|
|||||
|
|
|
|
dt |
|
|
M - M с |
= J × |
dw |
, |
||||
|
|
|
|
||||||
|
|
|
|
|
|
|
dt |
||
где uя,iя – напряжение и ток якорной цепи, В и А; |
|||||||||
Rя |
- активное сопротивление якоря, Ом; |
||||||||
Lя |
- |
индуктивность якоря, Гн; |
|||||||
k |
- коэффициент, зависящий от конструктивных данных двигателя; |
||||||||
Ф |
- |
магнитный поток, Вб; |
|||||||
ω |
- |
угловая скорость двигателя, 1/с; |
|||||||
M |
- |
электромагнитный момент двигателя, Н·м; |
|||||||
Мс |
- |
момент сопротивления, Н·м; |
|||||||
J |
- |
момент инерции якоря двигателя, кг·м2. |
|||||||
Перепишем данную систему относительно производных: |
|||||||||
|
|
ì |
di |
я |
|
|
|
|
|
|
|
ïLя × |
|
=u я -R я × i я - k × Ф ×w, |
|||||
|
|
|
|
||||||
|
|
í |
dt |
||||||
|
|
|
|
|
dw |
|
|
||
|
|
ï |
|
J × |
= М - М с . |
||||
|
|
|
|
||||||
|
|
î |
|
|
|
dt |
Обозначим здесь |
M = k ×Ф×iя; |
М с |
= |
|
М н |
|
×w и разделим первое уравне- |
||||||||||||||||
wн |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ние на Lя, а второе – на J, получим окончательное выражение дифференциаль- |
|||||||||||||||||||||||
ных уравнений системы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ì di |
я |
= |
|
u |
я |
- |
R |
я |
× i |
|
- |
k × |
Ф |
× w , |
|||||||||
ï |
|
|
|
|
|
|
я |
|
|
|
|
||||||||||||
|
|
|
|
L я |
L |
|
L я |
||||||||||||||||
ï |
|
dt |
|
|
я |
|
|
|
|
|
|||||||||||||
í |
|
|
d w |
|
|
|
k × Ф |
|
|
|
|
|
М н |
|
|
|
|||||||
ï |
|
|
= |
|
× i я |
- |
|
|
|
×w . |
|||||||||||||
ï |
|
|
dt |
|
|
J |
|
|
|
J × w |
|
||||||||||||
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10
5.2.Система дифференциальных уравнений в переменных состояния
Введём переменные состояния iя=х1, ω=х2 и обозначим
|
|
|
|
|
|
|
|
|
|
|
|
æ |
|
|
|
R |
|
|
|
|
|
|
|
|
k × Ф |
ö |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
æ |
а |
а |
ö |
|
|
|
ç |
- |
|
|
|
|
|
я |
|
|
|
|
- |
|
|
|
|
|
÷ |
|
|
|
|
|
æ 1 |
ö |
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
ç |
|
11 |
|
|
|
12 |
÷ |
|
|
|
ç Lя |
|
|
|
|
|
|
|
|
|
Lя |
÷ |
|
|
|
|
|
ç |
|
|
÷ |
|
|
|
|
u я |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
А = ç |
а |
|
а |
|
|
÷ = |
ç k × Ф |
|
|
|
|
|
|
|
|
М н |
÷; |
|
|
В = |
ç Lя |
÷; |
|
|
|
|
= U . |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
è |
|
21 |
|
|
|
22 |
ø |
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
|
|
|
|
|
|
|
÷ |
|
|
|
|
|
ç 0 |
÷ |
|
|
|
|
L я |
|
|||||||||
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
è |
|
|
|
|
|
|
|
|
|
|
|
J ×wн ø |
|
|
|
|
è |
|
|
ø |
|
|
|
|
|
|
|||||||||||||||||
Тогда получим систему |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ì |
dx1 |
= U + a × x + a × x , |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï dt |
|
|
|
11 |
|
1 |
|
12 |
2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
í |
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
2 |
= a21 |
×x1+a22 × x2. |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
î |
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
При наших данных входящие сюда величины равны |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
a |
|
|
= - |
|
|
Rя |
; a |
|
|
= - |
0,73 |
|
|
|
= -27,863 |
Ом |
; |
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
11 |
|
|
|
|
|
|
|
Lя |
|
|
11 |
|
|
|
|
|
|
0,0262 |
|
|
|
|
|
|
|
|
|
Гн |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
a |
|
|
|
= - |
|
|
k × Ф |
; a |
|
|
|
= - |
|
2,4897 |
= -95,029 |
В × с |
, |
|
|
|
|||||||||||||||||||||||||||||||||
|
|
12 |
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Lя |
|
|
|
|
|
|
|
|
|
|
|
0,0262 |
|
|
|
|
|
|
Гн |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
k ×Ф = |
|
uн - I н × Rя |
, |
|
k × Ф = |
220 - 33,5 × 0,73 |
= 2,4897 В × с; |
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
wн |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
78,54 |
|
В × с |
|
|
|
|
|||||||||||
|
|
|
a 21 |
= |
k × Ф |
; a |
|
21 |
= |
2,4897 |
|
|
= 15 ,369 |
|
|
; |
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
0,1625 |
|
|
кг × |
м 2 |
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
a |
22 |
|
= - |
|
|
М |
н |
|
|
|
|
; a 22 = - |
|
76,394 |
= -6,004 |
|
1 |
, |
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
J ×w н |
12,723 |
|
с |
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
М н |
= |
Рн |
; М н |
= |
6000 |
= 76,394Н × м. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
w н |
|
|
|
|
|
|
|
|
|
|
|
|
78,54 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Витоге получилась система
ìdx1 = U - 27,863× х - 95,029× х ,
ï dt 1 2
í |
dx |
|
|
|
ï |
2 |
= 15,369 |
× х1 |
- 6,004× х2. |
|
||||
î |
dt |
|
|
11
5.3. Определение собственных чисел и векторов
Найдём собственные числа матрицы А. Для этого составим характеристическое уравнение
- 27,863 - l - 95,029 = 0. 15,369 - 6,004 - l
Раскроем определитель: (-27,863-λ)·(-6,004-λ)-(-95,029·15,369)=0.
Решим уравнение относительно λ:
λ2+33,867·λ+1627,79=0,
получим корни уравнения:
λ1,2=-16,993±36,620·i.
Найдём собственные векторы:
λ1=-16,933+36,62·i,
æ - 27,863 |
+ 16,993 - 36,62 × i |
- 95,029 |
ö |
æV (1) |
ö |
æ |
0 |
ö |
|
ç |
|
|
÷ |
× ç |
1 |
÷ |
= ç |
|
÷. |
ç |
15,369 |
|
÷ |
ç |
(1) |
÷ |
ç |
0 |
÷ |
è |
- 6,004 +16,993 - 36,62 ×i ø |
èV2 |
ø |
è |
ø |
Уравнение для нахождения координат вектора имеет вид
|
|
|
(-10,93 - 36,62 × i) ×V (1) |
- 95,029 ×V (1) |
= 0. |
|||||||
|
|
|
|
|
|
|
|
|
1 |
|
2 |
|
Примем |
V (1) |
= 2,595, тогда V (1) |
= -0,298 - i. |
|
||||||||
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
Тогда выражение для х1 |
примет вид |
|
|
|
||||||||
х |
|
= |
æ |
х |
11 |
ö |
æ |
2,595 |
|
ö |
× е ( -16 ,933 + 36 , 62 ×i )×t . |
|
|
ç |
|
÷ = |
ç |
|
|
÷ |
|||||
|
1 |
|
ç |
х 21 |
÷ |
ç |
- 0,298 |
- i |
÷ |
|
|
|
|
|
|
è |
ø |
è |
ø |
|
|
Аналогично найдём решение х2.
Переходя к тригонометрической форме комплексного числа и вводя постоянные интегрирования, получим общее решение:
Х 1 = е-16,933×t × (2,595 ×С1 × cos 36,62 ×t + 2,595 ×C2 × sin 36,62 ×t),
Х 2 = е -16,933×t × ((-0,298 × С1 - С 2 ) × cos 36,62 × t + (C1 - 0,298 × С 2 ) × sin 36,62 × t ).
5.4. Частное решение неоднородной системы
По заданным частным решениям системы находим константы:
при t0=0 x1=3,35 A, x2=7,854 1/c.
Составляем систему
ì |
3,35 = е0 × (2,595 × С1 ×1 + 2,595 × С2 × 0), |
|
|
|||
í7,854 = е0 × ((-0,298 × С |
- С |
2 |
) ×1 + (С - 0,298 × С |
2 |
) × 0). |
|
î |
1 |
|
1 |
|
12
Решая систему, находим С1=1,291 и С2=-8,232. Подставляя в общее решение, находим частное решение:
Х 1 = е -16 , 933 ×t × (3,35 × cos 36 ,62 × t - 21,362 × sin 36 ,62 × t ),
Х 2 = е -16 ,933 ×t × (7,847 × cos 36 ,62 × t + 3,744 × sin 36 ,62 × t ).
Строим графики, используя систему MathCAD.
|
|
|
|
|
|
|
9 |
e |
16.933 |
.t . |
. |
. |
21.362 |
. |
. |
|
( 3.35 cos( 36.62 t ) |
sin ( 36.62 |
t ) ) |
||||
|
16.933 |
.t . |
|
|
|
|
0 |
e |
. |
. |
3.744 |
. |
. |
||
|
( 7.847 cos( 36.62 t ) |
sin ( 36.62 |
t ) ) |
||||
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
13 |
|
|
|
|
|
|
|
0 |
i(t) , W(t)
0.1 |
0.2 |
0 |
t |
0.3 |
|
t |
|
5.5. Фазовые траектории системы
Построение фазовых траекторий выполнено на MathCADе.
Исследование системы: х'=-27,863*х-95,029*у
у'= 15,369*x- 6,004*y
ORIGIN |
|
1 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
A |
|
|
|
|
|
27.863 |
|
95.029 |
|
eigenvals ( A ) = |
|
|
16.934 + |
36.62i |
||||
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|||||||||||||||||
|
|
15.369 |
6.004 |
|
|
|
||||||||||||
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
16.934 |
|
|
36.62i |
||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Собственные значения комплексные -16,934+36,62*i и -16,934-36,62*i
и реальная часть <0, поэтому точка покоя является "устойчивым фокусом"
|
|
|
|
|
|
27.863.x |
|
95.029.x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
f( t , x) |
|
|
1 |
2 |
|
a |
|
0 |
b |
|
20 |
h |
|
0.01 |
||||
|
|
|
|
|
15.369.x |
|
6.004.x |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
3.35 |
|
X |
|
rkfixed x, a , b , |
b |
, f |
|
|
|
7.854 |
|
|
|||||
|
|
|
||||||||
|
|
|
|
|
|
|
h |
13
На рис.2 представлена фазовая траектория системы.
8 |
|
|
|
5 |
|
|
|
X<3> |
|
|
|
0 |
|
|
|
3 |
|
|
|
10 |
5 |
0 |
5 |
12 |
|
X<2 > |
3.35 |
Рис.2. Фазовая траектория системы |
|