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

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

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

282

24. РАСЧЕТ НЕЛИНЕЙНЫХ ЦЕПЕЙ С ПОМОЩЬЮ АППРОКСИМАЦИИ ХАРАКТЕРИСТИК

 

u

 

 

 

Участок «0–1» (0 t

t1)

 

1

 

 

 

 

 

 

 

 

u (i )

,

 

 

 

u1

 

 

2

u (i ) i r

 

 

 

 

 

 

 

r u1

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

u (i )

i

 

 

i1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

i1

 

 

Участок «1–2»

t1

t T

2

 

 

 

 

 

 

 

 

 

 

 

 

const

 

 

 

 

 

 

 

u (i ) u1

 

 

 

u

 

u (i )

 

Участок «0–1» (0 t

t1)

 

1

 

2

 

 

 

 

 

 

 

 

 

u

(i ) 0 ,

 

 

 

 

u1

 

 

 

(r )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u (i )

 

i

 

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

2

 

0

 

 

 

Участок «1–2»

t T

 

 

 

 

 

 

 

const

 

 

 

 

 

 

 

u (i ) u1

 

 

 

 

 

 

Порядок расчета

 

 

 

 

 

 

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

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

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

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

283

Аппроксимация нелинейной характеристики степенным полиномом

Характер

Порядок определе-

Вид функции

аппроксимирующего

ния коэффициентов

при синусоидальном аргументе

полинома

степенного полинома

 

u

ai bi3

Если ток синусоидален

 

1

1

1

i(t) Im sin t ,

 

 

 

 

ai2

3

напряжение несинусоидально

u2

bi2

 

a ?

 

(содержит 1- и 3-ю гармоники)

 

b ?

 

u(t) U1m sin t U3m sin3 t

u ai bi3

u

 

ai bi3

ci5

Если ток синусоидален

 

1

1

1

 

1

i(t) Im sin t ,

 

 

 

ai

bi3

ci5

u

2

напряжение

несинусоидально

 

 

2

2

 

2

 

 

 

 

 

 

(содержит 1-, 3- и 5-ю гармони-

u

 

ai

bi3

 

ci5

 

3

3

3

 

3

ки)

 

 

 

 

a ?

 

 

u(t) U1m sin t U3m sin3 t

 

 

 

b ?

 

 

U5m sin5 t

 

 

 

c ?

 

 

u ai bi3 ci5

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

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

шил, С.В. Страхов. – М.: Энергия, 1989. – § 22.6, 25.4, 25.5.

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

М.: Гардарики, 2002. – § 15.12, 15.46–15.50.

284 24. РАСЧЕТ НЕЛИНЕЙНЫХ ЦЕПЕЙ С ПОМОЩЬЮ АППРОКСИМАЦИИ ХАРАКТЕРИСТИК

ПРИМЕРЫ

Задача 24.1

В цепи (рис. 24.1) e(t) = 100 sin 1000t B; r1 = r2 = 1 кОм; кулонвольтная ха-

рактеристика qC(uC) нелинейного конденсатора представлена на рис. 24.1; q0 = 10–5 Кл.

Рассчитать и построить q(t), uC(t), i(t).

i r1

ir

iC

e(t)

r2

C

uC

 

q0

q

b

c

q

b

c

q

0

 

 

q0

 

 

 

uC

 

 

 

uC

 

 

 

 

0

 

 

a

 

q0

 

 

d

d

 

a

 

 

 

 

 

 

 

 

Рис. 24.1

Решение

1. Нелинейные дифференциальные уравнения, описывающие цепь

(рис. 24.2):

i = ir + iC;

(1)

e = ir1 +uC;

(2)

uC = irr2;

(3)

iC dq dt CduC dt .

(4)

2. Кулонвольтная характеристика qC(uC) нелинейного

конденсатора

(рис. 24.1) представляет собой петлю гистерезиса с относительно малой площадью. Для упрощения расчета она может быть аппроксимирована ломаной dabc (рис. 24.1). Это дает возможность для каждого из участков da, ab и bc получить решение уравнений (1)–(4) методами линейной электротехники.

 

 

 

 

 

 

 

 

 

 

285

На участке da (рис. 24.1) q = const. Следовательно,

iC dq dt 0. Для

точки a этого участка (см. аппроксимированную кулонвольтную характери-

стику) uC = 0. Из уравнений (1)–(3) для

e(t)

 

 

 

 

 

этой же точки имеем: ir = 0; i = 0;

 

 

 

t1

 

 

t

 

e(t) = 100 sin t = uC = 0.

 

0

 

 

 

 

 

 

Т/4

Т/2

3Т/4

Т

Отсюда вытекает, что точка a соот-

q(t)

 

 

 

 

 

ветствует моменту времени t = 0. Кроме

+q0

 

 

 

 

 

того, в силу периодичности режима рабо-

0

 

 

 

 

t

ты цепи и симметрии характеристики

 

Т/4

Т/2

3Т/4

Т

 

 

q0

 

 

dabc (рис. 24.1) относительно начала ко-

 

 

 

 

 

ординат расчет может быть выполнен за

uC (t)

 

 

 

 

 

половину периода, т. е. для времени 0 t

 

 

 

 

 

t

T/2. Изменение заряда q в этом интер-

0

 

 

 

 

 

Т/4

Т/2

3Т/4

Т

вале

времени

определяется

участком

 

 

 

 

 

 

abcb

аппроксимированной

 

ха-

 

 

 

 

 

 

рактеристики (рис. 24.1), что соот-

i(t)

 

 

 

 

 

ветствует участку acb исходной харак-

 

 

 

 

 

 

теристики (рис. 24.1).

 

 

 

 

 

 

 

t

3. Определение q(t), uC(t), i(t).

 

0 t1

Т/4

Т/2

3Т/4

Т

 

 

 

На участке ab (0 t t1) uC = 0,

 

 

T/2+t1

 

 

 

q0 q +q0. Из уравнений (1)–(4) для

 

 

Рис. 24.2

 

 

этого участка получим

 

 

 

 

 

 

 

 

 

ir uC r2

0;

 

 

 

 

 

 

 

i e r1 100sinωt 1000 0,1sinωt (А);

 

 

 

 

 

iС = i = 0,1 sin t

A;

 

 

 

 

 

 

q iCdt 0,1sinωtdt 10 4 cosωt C1 (Кл).

 

(5)

На участке bc–cb (t1 t T/2)

q = +q0 = const. Из уравнений (1) (4) для

этого участка имеем: iC dq dt 0;

 

 

 

 

 

 

 

 

 

i ir e r1 r2 100sinωt 2000 0,05sinωt ;

 

 

 

 

uC = irr2 = 1000 0,05sin t =

50sin t (В).

 

 

 

28624. РАСЧЕТ НЕЛИНЕЙНЫХ ЦЕПЕЙ С ПОМОЩЬЮ АППРОКСИМАЦИИ ХАРАКТЕРИСТИК

Определение постоянной интегрирования С1 и момента времени t1 перехода рабочей точки режима с участка ab на участок bc характеристики

(рис. 24.1).

В точке а (t = 0)

q = q0 =–10–5 Кл.

Соотношение (5) для этого момента времени принимает вид

–10–5 = –10–4 cos(0) + C1,

из которого следует

С1 = 10–4 – 10–5 = 9 10–5 Кл.

В точке b (t = t1)

q = + q0 = 10–5 Кл.

Из соотношения (5) для момента времени t1 (с учетом найденного значе-

ния постоянной интегрирования С1) вытекает

10–5 = –10–4cos ( t1) + 9 10–5.

Полученное тождество позволяет определить временной момент t1 перехода рабочей точки с участка ab на участок bc характеристики (рис. 24.1):

t

 

1

arccos

8 10 5

 

 

1

arccos(0,8)

0,6435

0,6435 10 3 с.

 

10 4

1000

 

1

 

 

 

1000

 

4. Временны´е зависимости q(t), uC(t), i(t) представлены на рис. 24.2.

Задача 24.2

В цепи (рис. 24.3) ток источника электри-

ческой энергии i(t) = 4 sin t, где = 1000 с–1.

Линейное сопротивление r = 5 Ом. Вебер-

амперная характеристика (i) нелинейной ин-

дуктивности аппроксимирована ломаной dcab

(рис. 24.4),

Рис. 24.3 0 = 10–2 Вб; I0 = 2 А.

Рассчитать и построить (t), uL(t), uab(t).

 

 

 

 

 

287

Решение

 

 

 

 

 

1. Графический расчет рассматриваемой цепи базируется на решении сле-

дующей системы нелинейных дифференциальных уравнений:

 

 

 

 

ir uL uab;

 

 

 

 

(1)

uL d dt .

 

 

 

 

(2)

2. При t = 0 ток источника энергии i(0) = 4 sin(0) = 0.

 

 

 

a

b

Этому моменту времени соответствует точка 0 характе-

 

 

0

 

 

 

 

ристики (i) нелинейного элемента (рис. 24.4). Вебер-

 

I0

 

 

i

амперная характеристика (i) нелинейного элемента

 

0

I0

симметрична относительно начала координат. При пе-

 

 

 

 

 

 

 

 

риодическом изменении тока i(t) источника энергии в

 

 

0

 

 

течение первой половины периода (в пределах которого

d

c

 

 

 

 

 

i 0) рабочая точка режима перемещается из точки 0 в

 

Рис. 24.4

 

точку b, и обратно из точки b в точку 0 вдоль участков

 

 

 

 

 

 

 

0a и ab (рис. 24.4) вебер-амперной характеристики (i).

 

 

 

 

 

В течение второй половины периода (в пределах которого i

0) рабочая

точка режима перемещается из точки 0 в точку d, и обратно из точки d в точку

0 вдоль участков 0c и cd вебер-амперной характеристики (i). Поскольку в

пределах одного полупериода рабочая точка дважды перемещается вдоль уча-

стков характеристики 0a и ab (или 0c и cd), искомые временные зависимости

(t), uL(t), uab(t) могут быть рассчитаны при рассмотрении лишь четверти пе-

риода, в пределах которого рабочая точка перемещается из точки 0 в точку b

вдоль участков 0a и ab.

 

 

 

 

 

3. Определение (t), uL(t), uab(t).

 

 

 

 

 

На участке 0a (0 t t1) характеристика (i) представляется соотноше-

нием

=

0

i

10 2

i 0,5 10 2i.

 

I0

 

 

 

2

 

 

С учетом того, что i(t) = 4 sin 1000t А,

 

(t) = 0,5 10–2i(t) = 0,5 10–2 4 sin t = 0,02 sin t Вб.

(3)

288 24. РАСЧЕТ НЕЛИНЕЙНЫХ ЦЕПЕЙ С ПОМОЩЬЮ АППРОКСИМАЦИИ ХАРАКТЕРИСТИК

Напряжение на зажимах нелинейной индуктивности в рассматриваемом интервале времени (см. уравнение (2))

uL(t) d dt = 0,02 cos t = 20 cos t В,

а напряжение uab на зажимах источника тока (см. уравнение (1))

uab(t) = i(t)r + uL(t)= 5 4 sin t + 20 cos t = 202 sin ( t + 45º) В.

На участке ab (t1 t T/2) характеристика (i) представляется выраже-

нием = 0 = 10–2 Вб. Следовательно, в рассматриваемом диапазоне времени

uL(t) d dt d 0 dt

= 0,

а

напряжение на зажимах источника

тока

(см. уравнение (1)) uab(t) = i(t)r = 5 4 sin t = 20 sin t В.

 

Определение времени t1

перехода рабочей точки с участка 0a на участок

ab вебер-амперной характеристики (i).

 

 

В точке a характеристики (i) (момент времени t1) потокосцепление

должно достигнуть значения (t1) = 0 = 10–2 Вб.

 

С учетом выражения (3),

 

 

 

 

 

 

 

(t1) = 0,02 sin t1 = 0 = 10–2 Вб.

(4)

Из уравнения (4) следует

 

 

 

 

 

 

 

 

 

1

 

 

10

2

 

 

t

 

arcsin

 

0,524 10 3 с.

 

 

 

 

 

 

1

1000

 

 

0,02

 

 

 

 

 

 

 

Характеристики (t), uL(t), uab(t) представлены на рис. 24.5.

Задача 24.3

На рис. 24.6 изображена схема однополупериодного выпрямителя переменного тока. Напряжение источника

u(t) = 1152 sin t В,

где = 314 с–1; Е0 = 42 В; r = 200 Ом.

Вольт-амперная характеристика вентиля uв(i) представлена на рис. 24.7. Рассчитать и построить ur(t) и uв(t).

289

 

(t)

 

 

 

 

 

 

Рис. 24.5

 

 

 

i

 

мА

i

 

 

 

10

I0

a

b

 

 

 

 

0

 

 

 

u(t)

uв

50

 

 

 

 

r

ur

 

 

 

 

E0

c

 

 

uв

 

 

0

 

40

В

 

Рис. 24.6

 

Рис. 24.7

 

290

24. РАСЧЕТ НЕЛИНЕЙНЫХ ЦЕПЕЙ С ПОМОЩЬЮ АППРОКСИМАЦИИ ХАРАКТЕРИСТИК

Решение

 

1.

Нелинейное уравнение, описывающее работу цепи (рис. 24.6),

 

где

ur(t) + uв(t) = u(t) – Е0,

(1)

 

 

 

 

 

ur(t) = ri(t).

 

2.

В момент времени t = 0

 

 

u(t) = 115

 

sin t = 0

 

 

2

 

и, следовательно, i(t) = 0.

Из уравнения (1) для данного момента времени следует, что

uв = –Е0 = – 42 В 0.

Таким образом, первым рабочим участком вольт-амперной характеристики uв(i) нелинейного элемента является участок c0 (рис. 24.7).

На участке c0 (0 t t1) i = 0 и, следовательно, ur = ri = 0. Из (1) вытека-

ет, что

 

uв(t) = u(t) – Е0 = 115 2 sin t – 42 (В).

(2)

Точка 0 вольт-амперной характеристики соответствует моменту времени t1 перехода рабочей точки с участка с0 на участок 0а (рис. 24.7). В этой точке

uв(t1) = 0. Тогда из (2) получим

 

0 = 115 2 sin t1 – 42 (В), откуда

 

t

1

arcsin

 

42

 

 

 

1

arcsin0,2586 8,33 10 4с = 0,833 мс.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

ω

115 2

 

 

314

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в диапазоне времени 0 t t1 = 0,833 мс

 

 

 

 

 

ur(t) = 0; uв(t) = 115

 

sin t – 42 (В).

 

 

 

 

 

2

 

На участке 0а (t1 t t2) вольт-амперная характеристика вентиля пред-

ставляется линейной зависимостью (рис. 24.7)

 

 

 

 

 

 

 

 

i

I0

 

u

 

100 10 3

u

 

 

2,5 10 3u

 

.

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U0

 

в

40

 

 

в

 

в

 

 

Из уравнения (1) с учетом (3) имеем:

uв(t) = u(t) – Е0 ur(t) = u(t) – Е0 ri(t) = u(t) – Е0 2,5 10–3uв(t)r,

откуда в пределах рассматриваемого временнóго диапазона t1 t t2

 

 

 

 

 

 

 

 

 

 

 

291

 

 

u

в

 

u(t) E0

 

115

 

2sinωt 42

108sinωt 28 (В).

(4)

 

 

 

 

 

1 2,5 10 3r

 

 

 

 

 

 

 

 

1 2,5 10 3 200

 

Соотношение (3) с учетом (4) позволяет определить выражения для тока в

цепи:

 

 

 

 

 

 

 

 

 

 

 

i(t) = 2,5 10–3u (t) = 2,5 10–3(108 sin t – 28) = 270 sin t – 70 (мА),

(5)

 

 

 

 

 

в

 

 

 

 

 

 

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

 

t

1

t t

2

u (t) = ri(t) = 200(270 sin t – 70)10–3 = 54 sin t – 14 (В).

 

 

 

 

r

 

 

 

 

 

 

Выражение (5) для тока i(t) во временнóм диапазоне t1 t t2 позволяет определить время t2 перехода рабочей точки с участка 0а на участок ab характеристики uв(i) вентиля (рис. 24.7). В точке а, соответствующей моменту времени t2,

 

 

i(t2) = 270 sin t2 – 70 = 100 мА, откуда

1

170

1

 

 

0,68

2,17 10 3 с.

t2

 

arcsin

 

 

 

 

 

arcsin(0,63)

 

ω

 

314

314

 

 

270

 

 

 

 

На участке ab – ba (t2 t t3) i(t) = I0 = 100 мА.

 

 

Следовательно, для этого диапазона времени

 

 

 

 

 

 

ur(t) = ri(t) = 200 100 10–3 = 20 В.

 

 

Тогда, как следует из уравнения (1),

 

 

 

 

uв(t) = u(t) – Е0 ur(t) = 115

 

 

sin t – 42 – 20 = 115

 

sin t– 62 (В). (6)

 

2

2

В точке а (перехода рабочей точки с участка ab – ba на участок а0), соответствующей моменту времени t3, uв = 40 В. Для этой точки выражение (6) приобретает следующий вид:

uв(t3) = 1152 sin t3 – 62 = 40 (В),

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

четверти периода,

 

 

 

 

 

 

 

 

 

 

1

 

62 40

1

 

T

t3

 

 

 

arcsin

 

 

 

 

 

 

arcsin(0,628)

 

2,16 10 3 с.

 

 

 

 

 

2

115

2

 

 

 

 

 

314