Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
МОТС Вариант кнтр. и решения.pdf
Скачиваний:
5
Добавлен:
02.06.2015
Размер:
241.42 Кб
Скачать

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. Фазовая траектория системы