Кравченко. Практикум
.pdf322 |
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.
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.
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 в конце рассматриваемого интервала;
далее расчет повторяется до тех пор, пока потокосцепление и ток в конце очередного временнóго интервала не совпадут с требуемой точностью со значениями этих функций в начале рассматриваемого интервала (фактически расчет заканчивается тогда, когда очередное найденное значение тока

(
(