Кравченко. Практикум
.pdf282 |
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.
|
|
|
|
|
|
|
|
|
|
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 |
|
|