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

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

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

322

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

ЗАДАЧИ ДЛЯ САМОСТОЯТЕЛЬНОГО РЕШЕНИЯ

Задача 26.5

Вдоль линии с волновым сопротивлением zc1 = 500 Ом распространяется

волна uпад1(t)= U0(1 eαt); U0 = 500 кВ; α = 25·103 Нп/с. Линия через индуктивность L = 5 мГн (последовательное включение) соединена со второй лини-

ей, волновое сопротивление которой zc2 = 300 Ом.

Определить напряжение преломленной uпад2 и отраженной uотр1 волн.

Ответ: uпад2(t) = 375 444е–25000t + 69e–160000t кВ;

uотр1(t) = 125 + 240е–25000t 115e–160000t кВ.

Задача 26.6

К воздушной линии длиной l = 30 км с неразветвленной активноемкостной нагрузкой (r = 100 Ом, С = 0,25 мкФ) подключается генератор по-

стоянного напряжения U0 = 35 кВ. Волновое сопротивление линии zc = 400 Ом. Определить напряжение на сопротивлении r в момент времени, когда вол-

на, отразившись от конца линии, дойдет до ее середины.

Ответ: ur = 9,38 кВ.

27.ПЕРЕХОДНЫЕ ПРОЦЕССЫ

ВНЕЛИНЕЙНЫХ ЦЕПЯХ

Основные сведения

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

 

 

Метод графического интегрирования

 

 

 

Метод применим к цепям с разделяющимися переменными т. е. в основном к нелиней-

ным цепям постоянного тока. Суть метода состоит в следующем:

 

 

 

1) разделяются переменные

 

1

 

 

 

d f ( )

dt 1

 

 

 

 

d ;

f ( )

 

 

 

dt

 

f ( )

 

 

 

 

 

 

 

 

1

и

 

b

 

2) строится функция искомой переменной

 

 

 

 

 

 

f ( )

a

 

 

 

определяется время переходного режима

 

 

 

 

 

 

 

 

 

t1

1

1

 

f

 

 

 

 

 

 

 

 

t1 dt

d S0abcmf m

 

 

c

 

0

0

f ( )

 

0

d

1

 

 

 

 

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

Метод кусочно-линейной аппроксимации

 

3

Нелинейная характеристика заменяется ломаной (чем

 

больше линейных участков ломаной, тем точнее расчет). Для

2

каждого из линейных участков ломаной нелинейное диффе-

 

 

 

 

 

ренциальное уравнение вырождается в линейное. Решения

1

 

линейных дифференциальных уравнений для каждого из

 

 

линейных участков «сшиваются» друг с другом в точках

 

 

излома ломаной исходя из законов коммутации

i

0

324

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

Метод последовательных интервалов

Нелинейное дифференциальное уравнение d f ( ) , где искомая переходная dt

функция, представляется в конечно-разностной форме

k 1 k fk ( k )t t

при допущении, что производная определяет функцию в начале временно´го интервала. По значению искомой функции в начале произвольно выбранного временного интервала t с помощью уравнения

k 1 k t fk ( k )

отыскивается значение переходной величины в конце рассматриваемого интервала, являющееся начальным для последующего временного интервала. Расчет ведется шаг за шагом до получения с заданной погрешностью установившегося значения переходной функции. Точность расчета определяется величиной выбранного временного интервала t

Метод переменных состояния

Состоит в решении тем или иным методом (главным образом численным) системы дифференциальных уравнений 1-го порядка, в левой части которых представлены производные функций, не допускающих скачкообразных изменений (так называемых переменных состояния iL èëè ,uC èëè q ), характеризующих энергетический режим цепи. Число урав-

нений состояния цепи определяется числом элементов, независимо накапливающих энергию. Каждое из уравнений состояния представляется в конечно-разностной форме при допущении, что производная определяет функцию в начале временного интервала t . Система ко- нечно-разностных уравнений решается многократно. При этом каждый раз при анализе очередного временно´го интервала t определяются значения искомых функций в конце рассматриваемого интервала по их значениям в начале интервала. Найденные значения искомых функций рассматриваются как начальные при анализе последующего интервала. Расчет ведется шаг за шагом до получения с заданной погрешностью установившегося значения переходной функции

Библиографический список к разделу 27

1.Зевеке Г.В. Основы теории цепей / Г.В. Зевеке, П.А. Ионкин, Л.В. Не-

тушил, С.В. Страхов. – М.: Энергия, 1989. – § 15.5–15.7.

2.Анализ переходных процессов в нелинейных цепях методом переменных состояния / Л.И. Малинин, Н.А. Юрьева. – Новосибирск, 1989.

3.Бессонов Л.А. Теоретические основы электротехники / Л.А. Бессонов. –

М.: Гардарики, 2002. – § 16.2–16.5.

325

ПРИМЕРЫ

Задача 27.1

Методом графического интегрирования найти ток в цепи (рис. 27.1) после коммутации.

Параметры цепи: L = 10 Гн; r = 200 Ом; Е = 300 В. Вольт-амперная характеристика нелинейного резистивного сопротивления uн(i) представлена на рис. 27.2.

Решение

Методом графического интегрирования решаются лишь дифференциальные уравнения, допускающие разделение переменных.

В соответствии со вторым законом Кирхгофа рассматриваемая цепь описывается следующим дифференциальным уравнением:

di

E L dt uн(i) ri.

Рис. 27.1

Данное уравнение может быть представлено в виде, удобном для интегрирования:

А

i

 

 

)

 

 

 

 

 

 

 

(i

 

 

 

1

 

 

 

uн

 

 

 

iуст

 

 

 

 

 

 

 

0,8

 

 

 

 

 

 

 

0,6

 

 

 

 

E-ri

 

 

 

 

 

 

 

0,4

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

0

60

120

180

240

300

B

 

 

uнуст

 

ur

 

 

 

 

 

 

 

уст

 

 

 

 

 

 

Рис. 27.2

 

 

 

di E uн(i) ri uL (i), dt L L

откуда следует: dt di (i).

После разделения переменных появляется возможность определения времени переходного режима, соответствующего тому или иному значению переходной величины, по соотношению

t

i

di

.

(1)

t dt

 

(i)

0

0

 

 

В соответствии с уравнением (1) строится зависимость 1 (i) f (i) и определяется время, соответствующее конкретному значению переходного тока через площадь, ограниченную кривой f(i) и осью абсцисс.

326

 

 

 

 

 

 

 

 

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

с

1

 

 

 

 

 

 

 

На

рис.

27.2

показано

графическое

 

 

 

 

 

 

построение, из которого находится

 

А

(i)

 

 

 

 

 

 

 

0,4

 

 

1

 

f (i)

 

 

 

 

uL (E ri) uн(i),

 

 

 

 

 

 

 

 

 

 

 

 

 

0,3

 

(i)

 

 

 

на рис. 27.3. построена кривая 1 (i) f (i).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

Все расчеты, необходимые для построения

0,1

 

 

 

b

 

 

 

зависимости f(i) и функции i(t), сведены в табл.

 

 

 

 

 

 

 

27.1.

 

 

 

 

а

S

 

 

с

 

 

i

 

 

 

 

d

0,2

0,4

 

0,6

0,8

1

A

 

При подсчете площади, ограниченной кри-

0

 

 

 

Рис. 27.3

 

 

вой

f (i) 1 (i)

и

осью абсцисс,

пропорцио-

нальной времени t, необходимо учитывать масштаб тока (mi) и масштаб функ-

ции 1 (i) mf (i) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 27.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

uн(i)

E ri

 

uL= E – uн(i)

(i)

uL

 

 

1

 

 

t

 

 

 

(i)

 

 

 

 

 

– ri

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

В

В

 

 

 

 

В

А/с

 

 

с/А

 

с

0

 

0

300

 

 

 

300

30

 

 

0,033

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,1

 

7,5

280

 

 

 

272,5

27,25

 

0,037

 

0,35·10–2

0,2

 

15

260

 

 

 

245

24,5

 

0,041

 

0,74·10–2

0,3

 

22,5

240

 

 

 

217,5

21,75

 

0,046

 

1,175·10–2

0,4

 

33

220

 

 

 

187

18,7

 

0,0535

 

1,673·10–2

0,5

 

47

200

 

 

 

153

15,3

 

0,0654

 

2,27·10–2

0,6

 

60

180

 

 

 

120

12

 

 

0,0833

 

3,01·10–2

0,7

 

77,5

160

 

 

 

82,5

8,25

 

0,121

 

4,03·10–2

0,8

 

100

140

 

 

 

40

4

 

 

0,25

 

 

5,89·10–2

 

 

 

 

ik

1

 

 

 

 

 

 

 

 

 

 

 

В результате tk

 

 

di Skmimf (i) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

(i)

 

 

 

 

 

 

 

 

(В качестве примера на рис. 27.3 показан подсчет площади S, соответствующей времени переходного режима, по истечении которого ток дости-

гает значения 0,6 А.) График изменения искомого тока i(t) представлен на рис. 27.4.

327

A

i

 

 

iуст

 

 

1

 

 

 

 

 

0,8

 

 

 

 

 

 

0,6

 

 

 

i(t)

 

 

0,4

 

 

 

 

 

 

0,2

 

 

 

 

 

t

 

 

 

 

 

 

0

 

 

 

 

 

6 10 2с

0

1

2

3

4

5

 

 

 

Рис. 27.4

 

 

Задача 27.2

Методом кусочно-линейной аппроксимации определить переходный ток i(t) в схеме, представленной на рис. 27.5: Е = 12,1 В; r = 7,5 Ом; L = 0,5 Гн. Вольт-амперная характеристика диода приведена на рис. 27.6.

A

i

 

 

 

 

 

 

0,8

I2

 

 

 

2

 

 

 

 

 

 

 

 

 

0,6

 

u

**(i

**)

 

 

 

 

 

д

д

 

 

 

0,4

 

 

 

 

 

 

 

 

I1

 

1

 

uд(iд)

 

 

0,2

 

 

 

 

 

uд*(iд*)

 

U1

 

U2

u

 

 

 

 

0

1

2

 

3

4

5

B

Рис. 27.5

 

Рис. 27.6

 

 

Решение

Переходный режим в рассматриваемой цепи описывается нелинейным дифференциальным уравнением первого порядка

E L

di

uд(i) ri.

(1)

 

 

dt

 

328

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

Для решения уравнения (1) реальная нелинейная характеристика диода аппроксимируется ломаной линией (рис. 27.6), состоящей из двух линейных участков 0–1 и 1–2. Уравнения для линейных участков:

 

u*

 

U1

i*

 

2,5

i*

12,5i* ;

(2)

 

 

0,2

 

д

 

 

I

д

 

 

д

д

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

uд**

 

 

 

 

U U

iд** I1

 

 

U1

 

 

2

1

 

 

 

 

I

I

 

 

 

 

 

 

 

 

 

2

1

 

 

 

= 2,5

4 2,5

iд** 0,2 2,5 2,5 iд** 0,2 .

(3)

0,8 0,2

 

 

 

 

 

 

 

 

 

 

 

Для каждого из линейных участков аппроксимированной характеристики нелинейное дифференциальное уравнение (1) вырождается в линейное, решение которого можно найти, например, классическим методом.

1. Для участка 0–1 аппроксимированной характеристики дифференциальное уравнение (1) с учетом (2) принимает следующий вид:

Ldi* i*(12,5 r) E , dt

или

0,5

di*

20i* 12,1.

(4)

 

 

dt

 

Уравнение (4) описывает переходный режим в интервале времени (0 t t1), где t1 соответствует моменту времени, когда переходный ток достигает значе-

ния I1 = 0,2 A.

Решение уравнения (4):

i*(t) iпр* А*ер*t .

 

 

 

 

 

 

 

di

*

 

Принужденная составляющая переходного тока при

 

0

 

 

 

 

 

 

dt

 

 

 

 

 

 

i*

 

12,1

0,605А.

 

 

 

 

 

 

 

пр

20

 

 

 

 

329

Характеристическое уравнение 0,5р* + 20 = 0 имеет решениер* = 40 с–1.

Постоянная интегрирования А* определяется с учетом закона коммута-

ции i*(0 ) i*(0 ) 0. В результате

A* iпр* –0,605 А.

Таким образом, переходный ток в интервале времени (0 t t1) определяется соотношением

i*(t) 0,605(1 е 40t ) .

(5)

Момент времени t1, к которому ток i*(t) должен достичь значения I1 = 0,2 A, может быть определен из выражения (5):

i*(t1) 0,605(1 е 40t1) 0,2А,

откуда t1 = 0,01 с.

2. Для участка 1–2 аппроксимированной характеристики (рис. 27.6) дифференциальное уравнение (1) с учетом (3) принимает следующий вид:

Ldi** i**(2,5 r) E 2,5 0,5 dt

или

0,5

di**

10i** 10,1.

(6)

 

 

dt

 

Уравнение (6) описывает переходный режим в интервале времени (t1 t ). Решение уравнения (6):

i**(t) i** А**ер**(t t1) .

 

 

(7)

 

пр

 

 

 

 

 

 

 

di

**

 

Принужденная составляющая переходного тока при

 

0

 

 

 

 

 

 

dt

 

 

 

 

 

 

i**

 

10,1

1,01 А.

 

 

(8)

 

 

 

пр

10

 

 

 

 

Характеристическое уравнение 0,5р** + 10 = 0 имеет корень р** = 20 с–1.

330

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

Постоянная интегрирования А** определяется из условия невозможности скачка переходного тока (первый закон коммутации) в момент времени t = t1, т. е. из условия

i**(t

) i*(t

) 0,2 A.

(9)

1

1

 

 

Из соотношения (7) с учетом (8) и (9) при t1 = 0,01 с следует: 1,01 А** 0,2 А.

A i

 

 

 

 

 

 

0,8

 

 

 

 

 

 

0,6

 

 

 

 

 

 

0,4

 

i**(t)

 

 

 

 

i1

 

 

 

 

 

 

0,2

 

i*(t)

 

 

 

 

 

t1

 

 

 

t

0

0,02

0,04

0,06

0,08

с

В результате

A** 0,2 1,01 0,81A.

Окончательное выражение для переходного тока в интервале времени (t1 t ) определится соотношением

Рис. 27.7

i**

(t) 1,01 0,81е 20(t t1)

А.

 

 

 

3. Рассчитанные зависимости искомого переходного тока представлены на рис. 27.7.

Задача 27.3

Используя метод последовательных интервалов (численный метод), определить переходные функции тока i(t) и потокосцепления (t) в цепи, представленной на рис. 27.8.

r = 200 Ом; E = 100 B.

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

тивности задана выражением i 675 4 , где [i] = A; [ ] = Вб.

(i)

i(t)

Рис. 27.8

Решение

Переходный режим в рассматриваемой цепи описывается нелинейным дифференциальным уравнением, составленным на основании второго закона Кирхгофа:

d

ri E.

(1)

 

dt

 

331

Уравнение (1) может быть приведено к виду

d

E ri.

(2)

 

dt

 

При допущении, что производная функции определяет ее значение в начале каждого из временны´х интервалов t , уравнение (2) может быть представлено в конечно-разностной форме:

 

d

 

k 1 k

E ri .

(3)

 

 

 

 

dt

t

k

 

 

 

 

Из тождества (3) следует, что

 

 

 

k 1

k (E rik ) t ,

(4)

где k 1 значение функции в конце рассматриваемого временнóго интервала;

k и ik – значения функций в начале рассматриваемого временнóго интервала. Уравнение (4) решается в следующем порядке:

задаются достаточно малым приращением t (чем меньше t, тем точ-

нее расчет);

в начале первого интервала (k = 0), как следует из первого закона коммутации и вебер-амперной характеристики, i0 = 0, 0 = 0;

значение потокосцепления в конце первого интервала (k = 0) определяется соотношением (4):

1 0 (E ri0) t E t .

По характеристике (i) определяется значение тока i1, соответствующее найденному значению потокосцепления 1;

значения 1 и i1, являющиеся начальными значениями функций на втором временнóм интервале (k = 1), позволяют определить с помощью уравнения

(4) значение потокосцепления в конце этого интервала:

2 1 (E ri1) t

ивслед за этим по характеристике (i) найти ток i2 в конце рассматриваемого интервала;

далее расчет повторяется до тех пор, пока потокосцепление и ток в конце очередного временнóго интервала не совпадут с требуемой точностью со значениями этих функций в начале рассматриваемого интервала (фактически расчет заканчивается тогда, когда очередное найденное значение тока