Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кравченко. Практикум

.pdf
Скачиваний:
88
Добавлен:
11.03.2016
Размер:
4 Мб
Скачать
☆

332

 

 

 

 

 

 

27. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ ЦЕПЯХ

 

совпадет

с заданной погрешностью с

установившимся значением тока

 

Iуст

E

 

100

0,5А).

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

Результаты расчета методом конечных элементов при выбранном значе-

 

нии t = 0,2 мс сведены в табл. 27.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 27.2

 

 

 

 

 

 

 

 

 

 

 

k

 

t

ik

k

E – ikr

k+1

 

ik+1

 

–

 

мс

А

Вб

В

Вб

 

А

 

0

 

 

 

0,2

0

0

100

0,02

 

1,08·10–4

 

1

 

 

 

0,4

1,08·10–4

0,02

99,98

0,0399

 

1,727·10–3

 

2

 

 

 

0,6

1,727·10–3

0,0399

99,7

0,0599

 

8,7·10–3

 

3

 

 

 

0,8

8,7·10–3

0,0599

98,27

0,0796

 

2,7·10–2

 

4

 

 

 

1

2,7·10–2

0,0796

94,62

0,0985

 

6,35·10–2

 

5

 

 

 

1,2

6,35·10–2

0,0985

87,34

0,116

 

1,22·10–1

 

6

 

 

 

1,4

1,22·10–1

0,116

75,6

0,131

 

1,99·10–1

 

7

 

 

 

1,6

1,99·10–1

0,131

59,5

0,143

 

2,83·10–1

 

8

 

 

 

1,8

2,83·10–1

0,143

43,6

0,152

 

3,58·10–1

 

9

 

 

 

2

3,58·10–1

0,152

28,46

0,157

 

4,15·10–1

 

10

 

 

 

2,2

4,15·10–1

0,157

17,2

0,16

 

4,52·10–1

 

11

 

 

 

2,4

4,52·10–1

0,16

9,6

0,1628

 

4,74·10–1

 

12

 

 

 

2,6

4,74·10–1

0,1628

5,4

0,1638

 

4,86·10–1

 

По данным расчета построены временны´е зависимости (t) и i(t), представленные на рис. 27.9. и 27.10.

Рис. 27.9

Рис. 27.10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

333

Задача 27.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью метода переменных состояния

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рассчитать и построить uC(t) в схеме (рис. 27.11).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

r

 

 

 

 

 

 

Параметры цепи:

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

 

 

r2

 

 

 

 

 

uC(q)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = 120 B; r = 30 Ом; r

= 20 Ом.

 

 

 

 

 

 

 

 

 

 

 

i2

 

iC

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кулонвольтная

характеристика нелиней-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 27.11

 

 

 

 

 

ного конденсатора задана в табл. 27.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 27.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

мкК

0

 

10

20

30

40

 

44

 

46

 

47

 

 

 

 

 

 

 

 

 

uс

В

0

 

13

27

43

63

 

80

 

100

 

120

 

 

 

 

 

 

 

 

Решение

Рассматриваемая цепь имеет первый порядок, поэтому для расчета переходного процесса методом переменных состояния требуется лишь одно уравнение состояния цепи. В качестве переменной состояния цепи целесообразно рассматривать заряд q на обкладках нелинейного конденсатора. С учетом это-

го уравнение состояния цепи должно иметь следующий вид: q f (q,E) .

t

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

 

 

 

 

 

 

 

 

i

 

i ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

C

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

E i1r i2r2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i r

 

u q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2

 

 

 

C

 

 

 

 

 

 

 

 

и разрешить ее относительно тока i

 

 

q

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

t

 

 

 

 

 

 

 

 

 

 

i

uC (q)

,

i

1

(E i r )

E uC(q)

;

 

 

 

 

 

2

r

1

r

 

 

2 2

 

 

 

r

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

dq

i

i i

E uC (q)

 

uC (q)

 

E

 

r1 r2

u

 

(q).

 

 

 

 

r

 

 

 

dt C

1

2

 

 

r

 

 

 

 

 

 

r

 

 

rr

C

 

 

 

 

 

 

 

 

1

 

 

 

2

 

1

1 2

 

 

 

 

Рис. 27.12

334

27. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ ЦЕПЯХ

После подстановки числовых значений

 

 

 

 

dq

4

uC(q)

.

(2)

 

 

12

 

dt

 

 

Полученное уравнение является уравнением состояния рассматриваемой цепи

(рис. 27.11).

 

 

 

 

 

 

 

 

 

2. Определение

области

изменения перемен-

 

i1

 

 

 

 

 

 

ных состояния.

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

uC(q)

Как следует

из докоммутационной схемы,

i2

 

iC

 

 

 

 

 

 

 

 

представленной на рис. 27.12, uC(0 ) = Е = 120 В.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании

второго

закона коммутации

uC(0+) = uC(0 ).

Следовательно, начальное значение заряда на обкладках конденсатора,

соответствующее uC(0+) = 120 В (как следует из кулонвольтной характеристики, представленной в табл. 27.4):

q(0+) = 47 мкК.

Установившееся значение напряжения на конденсаторе uCóñò определя-

ется по уравнению состояния (2) при условии dqуст 0. В результате dt

uCóñò = 48 В. Найденное значение установившегося напряжения предопределя-

ет (в соответствии с кулонвольтной характеристикой, представленной в табл. 27.4) установившееся значение заряда

qуст = 33 мкК.

Таким образом, пределы изменения переменной состояния q и напряжения на конденсаторе uC:

q [47; 33] мкК,

uC [120; 48] В.

3. Решение уравнения состояния (в качестве примера рассматривается метод численного интегрирования).

Согласно алгоритму Эйлера

 

 

u (t

k

)

 

 

q(tk t) q(tk ) t

4

C

 

 

q(tk ) qk ,

(3)

12

 

 

 

 

 

 

 

335

где k – номер временнóго интервала; t – шаг интегрирования (может быть увеличен, если q изменяется медленно); tk – время, соответствующее началу

k-го интервала; q

t

 

4

uC(tk )

 

– приращение переменной состояния на

 

 

k

 

12

 

 

 

 

 

 

 

k-м интервале; uC(tk) – напряжение на емкости в начале временнóго интервала, соответствующее величине заряда q(tk) и определяемое по кулонвольтной характеристике (табл. 27.4).

Результаты расчета представлены в табл. 27.4.

Таблица 27.4

Величина

Размерность

 

 

 

Номер интервала k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

3

4

5

6

 

7

8

9

 

 

 

 

 

 

 

 

 

 

 

 

t

мкс

0,5

1

1,5

1,5

1,5

2

 

2

2,5

2,5

tk

мкс

0

0,5

1,5

3

4,5

6

 

8

10

12, 5

qk

мкК

0

–3

–2,67

–2,5

–1,5

–0,875

 

–1

–0,67

–0,416

q(tk)

мкК

47

44

41,33

38,83

37,33

36,455

 

35,455

34,875

34,369

uC(tk)

В

120

80

68,7

60

55

54

 

52

50

49

В процессе расчета шаг интегрирования t изменялся ввиду медленного увеличения q. Расчет прекращен при достижении uС = 49 В, отличающегося от

uCóñò 48 В на 2 %.

Зависимость uC(t) представлена на рис. 27.13.

В

uC

 

 

 

 

 

 

 

140

 

 

 

 

 

 

 

 

120

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

60

 

 

 

 

 

 

 

 

40

 

 

 

 

uC

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

ï ð

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

4

6

8

10

12

мкс

 

 

 

 

Рис. 27.13

 

 

 

336 27. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ ЦЕПЯХ

Задача 27.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В схеме, изображенной на рис. 27.14,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E 180 В,

r1

50

Ом,

r2 r3 200 Ом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iL

( )

 

r4 1200 Ом.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вебер-амперная характеристика нели-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нейной индуктивности задана табл. 27.5.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя метод переменных состояния,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 27.14

 

 

 

 

 

 

рассчитать и построить iL(t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 27.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iL

 

А

 

 

0

0,1

 

0,2

 

0,3

0,4

0,5

 

0,6

0,8

 

1,0

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мВб

 

 

0

10

 

13,5

 

16,5

18

19

 

20

21

 

21,8

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение

1. Составление уравнения состояния цепи (рис. 27.14).

В качестве переменной состояния целесообразно выбрать потокосцепление . Для режима после коммутации составим систему уравнений по законам Кирхгофа и решим ее относительно uL d dt :

i1 i2 iL,

 

 

 

d

 

 

 

 

 

i

r

 

 

i r ,

dt

 

L 3

 

2 2

E i r i r .

 

 

1 1

2 2

E (i2 iL)r1 i2r2 iLr1 i2(r1 r2);

i

E iLr1

;

 

i

 

r

d

 

E iLr1

r 0;

 

 

 

 

 

2

 

r r

 

 

L 3

dt

 

 

 

r r

2

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

d

i

r

r1r2

 

 

E

r2

 

 

0;

 

 

r r

 

 

r r

 

dt

L

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

1

2

 

d 240iL 144 0; dt

 

 

 

337

 

d

144 240iL .

(1)

 

 

 

dt

 

Соотношение (1) является уравнением состояния цепи (рис. 27.14). 2. Определение области изменения переменной состояния .

Начальное значение потокосцепления определяется с помощью первого

закона коммутации.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из расчета цепи до коммутации (рис. 27.15)

 

 

 

 

iL(0 )

iL(0 )

 

 

E

 

 

 

 

r2

 

= 0,1 А;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

(r3 r4)

 

r

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

4

 

 

 

 

 

 

 

 

 

 

 

 

r1 r

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iL( 0)

iL( 0) 0,1 А.

 

Рис. 27.15

По вебер-амперной

характеристике определяем

(0 ) 10 мВб.

Установившееся значение потокосцепления находится из уравнения состояния при условии d dt 0 (также может быть найдено из расчета цепи в

установившемся послекоммутационном режиме).

0 144 240i

, i

L

 

144

0,6 А,

 

уст

20 мВб.

 

Lуст

 

240

 

 

 

 

 

уст

 

 

 

 

Таким образом, в переходном процессе:

iL 0,1; 0,6 А и 10; 20 мВб.

3. Решение уравнения состояния (в качестве примера рассматриваются два метода: метод численного интегрирования и метод кусочно-линейной аппроксимации характеристики нелинейного элемента):

а) метод численного интегрирования уравнения состояния в соответствии с алгоритмом Эйлера:

(tk t) = (tk )+ t 144 240iL tk = (tk ) (tk ),

где k – номер интервала; t – шаг интегрирования (он может быть увеличен, если будет изменяться медленно); tk – время, соответствующее началу k-го интервала; (tk ) t 144 240iL(tk ) – приращение переменной состояния;

338

27. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ ЦЕПЯХ

iL(tk ) – ток в индуктивности, соответствующий величине (tk ) и определяе-

мый по характеристике нелинейной индуктивности (табл. 27.5). Результаты расчета сведены в табл. 27.6.

Таблица 27.6

Значения

 

 

 

 

Номер интервала k

 

 

 

 

функций

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в начале

1

2

3

4

5

6

7

8

9

10

11

временного

интервала

 

 

 

 

 

 

 

 

 

 

 

t

мкс

–

10

10

10

10

10

20

20

30

40

40

 

 

 

 

 

 

 

 

 

 

 

 

tk

мкс

0

10

20

30

40

50

60

70

100

140

1890

(tk )

мВб

10

11,2

12,3

13,3

14,3

15,2

16

17,5

18,7

19,9

20

iL(tk )

A

0,1

0,134

0,166

0,196

0,227

0,257

0,285

0,37

0,466

0,59

0,6

(tk )

мВб

1,2

1,1

1,04

0,97

0,89

0,82

1,5

1,1

0,94

0,05

–

мВб

 

 

 

В процессе расчета шаг интегрирования t из-

 

 

менялся ввиду медленного увеличения . Расчет

24

 

 

 

16,5

2

3

прекращен при достижении (tk )

= 19,9978 мВб,

20

 

 

 

 

 

 

 

16

 

 

 

отличающегося от

уст 20 мВб на величину

12

 

 

 

10 1

 

менее 0,01 %;

 

 

8

 

 

 

4

0,1

iL

б) метод

кусочно-линейной

аппроксимации

 

0,3

характеристики нелинейного элемента.

 

0

0,2

0,4 0,6 0,8 1,0 1,2 A

 

 

 

Рис. 27.16

В пределах изменения переменной состояния

 

 

 

10; 20

мВб

характеристика нелинейного

 

 

 

 

элемента может быть аппроксимирована ломаной

сдвумя линейными участками (рис. 27.16):

1)участок 1–2:

10; 16,5 мВб;

iL 0,1; 0,3 А;

2) участок 2–3:

16,5; 20 мВб;

iL 0,3; 0,6 А.

339

Составим описание аппроксимированной характеристики iL( ) на каж-

дом из участков, используя уравнение прямой, проходящей через две точки.

Участок 1–2 t 0,t1 :

 

0,0165

=

0,01 0,0165

.

(2)

 

 

i 0,3

 

 

 

 

 

 

0,1 0,3

 

Решая уравнение (2) относительно i, получим

 

 

 

i 30,77 0,2077.

(3)

Участок 2–3 t t1, :

 

 

 

 

 

 

 

 

0,02

=

0,0165 0,02

,

 

 

 

 

 

 

 

 

i 0,6

 

 

0,3 0,6

 

откуда

 

 

 

 

 

 

 

 

i 85,7 1,11.

(4)

Расчет переходного режима в интервале (0 t t1).

Подставив уравнение (3) в уравнение состояния (1) для участка характеристики 1–2 , получим

d 7385 194 0. dt

При d dt 0 установившееся значения потокосцепления

уст 0,0263 Вб.

Характеристическое уравнение имеет вид p 7,385 103 0.

Корень характеристического уравнения p 7,385 103 (с–1).

Общее решение для искомого потокосцепления:

(t) 0,0263 Ae1 7,385103t .

340

27. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ ЦЕПЯХ

Определение постоянной

интегрирования (с учетом (0 ) (0 )

10 мВб):

 

(0) 0,01 0,0263 A1,

следовательно, A1 0,0163 Вб, и полное решение для потокосцепления в рассматриваемом временном интервале будет следующим:

(t) 0,0263 0,0163e 7,385103t Вб.

(5)

Из соотношения (3) с учетом (5) следует:

 

i (t) 0,6 0,5e 7,385103t А.

(6)

Момент времени t1 перехода рабочей точки с участка 1–2 на участок 2–3 определяется из (6) с учетом того, что при t t1

1 0,0165 Вб, i1 0,3(точка 2): 0,3 0,6 0,5e 7,385103t1 ,

откуда t1 69 мкс.

 

Расчет переходного режима во временном интервале (t1

t ).

Подставив уравнение (4) в уравнение состояния (1) для участка характе-

ристики 2–3, получим

 

 

d

20,6 103 410 0.

(7)

 

 

 

dt

 

Решение уравнения (7) проводится аналогично расчету, проводимому на первом интервале:

3

уст A2e 20,510 (t t1);

уст

0,02 Вб;

p 20,6 103с–1.

341

В силу непрерывности потокосцепления начальным условием для второго интервала является конечное значение на первом интервале, т.е.

 

(t )

(t ) = 0,0165 Вб.

 

 

1

1

 

При t t1

 

 

 

 

(t ) 0,0165 0,02 A ,

 

 

1

2

 

следовательно,

 

 

 

 

А2 = 0,0035

 

и

(t) 0,02 0,0035e 20,5103(t t1) .

(8)

После подстановки (8) в (4) получим

 

 

 

iL (t)=0,6 0,3e 20,5103(t t1) А.

(9)

Таким образом, искомое решение рассматриваемой задачи имеет следующий вид:

 

i

(t) 0,6 0,5e 7,385103t

А

при

t 0;t

,

 

L

 

 

 

1

 

i

(t) 0,6 0,3e 20,5103(t t1)

А

при

t t ; .

L

 

 

 

 

1

 

График зависимости iL(t) представлен на рис. 27.17.

i**L (t)

i*L(t)

10 5c

Рис. 27.17