метод расчета переходных процессов

advertisement
Оглавление
ВВЕДЕНИЕ ....................................................................................................................................2
Раздел 1. КЛАССИЧЕСКИЙ МЕТОД РАСЧЕТА ПЕРЕХОДНЫХ ПРОЦЕССОВ................3
Раздел 2. РАСЧЕТ ПЕРЕХОДНЫХ ПРОЦЕССОВ ПРИ ПРОИЗВОЛЬНЫХ ВХОДНЫХ
ВОЗДЕЙСТВИЯХ С ИСПОЛЬЗОВАНИЕМ ИНТЕГРАЛОВ НАЛОЖЕНИЯ ......................19
КОНТРОЛЬНЫЕ ВОПРОСЫ.....................................................................................................27
СПИСОК ЛИТЕРАТУРЫ...........................................................................................................27
ВВЕДЕНИЕ
Процесс перехода от одного энергетического состояния
электрической цепи к другому называется переходным процессом.
Переходный процесс вызывается коммутацией, т.е.мгновенным
изменением параметров цепи, ее схемы или параметров источников
энергии в схеме.
Переход от одного (докоммутационного) состояния к другому
обычно происходит не мгновенно, а в течение некоторого времени –
времени переходного процесса. Это объясняется тем, что каждому
состоянию
цепи
соответствует
определенный
запас
электромагнитной энергии. Изменение же энергии в реактивных
элементах не может происходить мгновенно, так как в этом случае
dw
мощность p =
, развиваемая в цепи, достигала бы бесконечно
dt
больших значений. Следовательно, не могут изменяться мгновенно
и переменные, связанные с энергией. Следствием этих положений
являются законы коммутации:
q(t − ) = q(t + ) , ψ (t − ) = ψ (t + ) .
Для линейных цепей законы коммутации чаще записывают так:
uC (t − ) = uC (t + ) , i L (t − ) = i L (t + ) .
Цель анализа переходных процессов в электрических цепях –
определение временных законов изменения токов или напряжений
на заданных участках цепи в переходном режиме.
Для расчета переходных процессов во временной области
используются два метода: классический метод и метод интегралов
наложения. Классический метод рекомендуется применять для
анализа
цепей,
процессы
в
которых
описываются
дифференциальными уравнениями не выше третьего порядка, при
действии в схеме постоянных или гармонических источников
энергии. Метод интегралов наложения – при действии источников
произвольной формы.
2
Раздел 1. КЛАССИЧЕСКИЙ МЕТОД РАСЧЕТА ПЕРЕХОДНЫХ
ПРОЦЕССОВ
Расчет переходных процессов классическим методом сводится
к решению системы линейных дифференциальных уравнений с
постоянными коэффициентами, составленных на основании законов
Кирхгофа для мгновенных значений токов и напряжений. Эта
система приводится к неоднородному дифференциальному
уравнению n-го порядка, общее решение которого
y (t ) = y частн. (t ) + y общ. (t ),
где y частн. (t ) - частное решение неоднородного дифференциального
уравнения, y общ . (t ) - общее решение однородного уравнения.
y (t ) понимается любой искомый ток или
Здесь под
напряжение.
Частное
решение
неоднородного
уравнения
определяется видом функции, стоящей в правой части уравнения, и
поэтому называется вынужденной составляющей yвын. (t ) . Для цепей
с постоянными или периодическими напряжениями (токами)
источников энергии вынужденное решение совпадает с
установившимися значениями искомых функций. Общее решение
однородного уравнения описывает электромагнитный процесс,
происходящий в схеме без воздействия внешних источников, и
называется свободной составляющей y св . (t )
Из теории дифференциальных уравнений известно, что решение
однородного уравнения ищется в виде
n
yсв (t ) = ∑ As A pst ,
s =1
где As - постоянные интегрирования, определяемые из начальных
условий; ps - корни характеристического уравнения.
y (t )
Основными составляющими расчета
являются:
определение начальных условий, вынужденной составляющей,
корней
характеристического
уравнения,
постоянных
интегрирования. Более подробно остановимся на определении
начальных условий и корней характеристического уравнения.
1.1.ОПРЕДЕЛЕНЕ НАЧАЛЬНЫХ УСЛОВИЙ
Независимые начальные условия определяются
законами
коммутации, которые для линейных цепей можно записать в виде
u C (0 + ) = u C (0 − ) , iL (0 + ) = iL (0 − ).
3
Здесь учитывается, что обычно момент коммутации совмещают
с началом отсчета, то есть полагают t k = 0.
Значения остальных токов и напряжений до и после
коммутации в общем случае не одинаковы:
iC (0 + ) ≠ iC (0 − ), u L (0 + ) ≠ u L (0 − ),
du C
du di
di
(0 + ) ≠ C , L (0 + ) ≠ L (0 − ) и т. д.
dt
dt di
dt
Их значения в момент времени t = 0 + определяются
независимыми начальными условиями, характером (видом)
коммутации и другими факторами. Поэтому они получили название
зависимых начальных условий.
Порядок расчета
1. Определяем независимые начальные условия.
Для схемы до коммутации, находящейся в установившемся
режиме, находят мгновенные значения токов в индуктивностях и
напряжений на емкостях, после чего подставляют в выражения для
u C (t ) и i L (t ) значение времени t = 0 − . В соответствии с законами
коммутации u C (0 − ) = u C (0 + ) , iL (0 − ) = iL (0 + ).
2. Для схемы, полученной после коммутации, составляем
систему уравнений по законам Кирхгофа для мгновенных значений
токов и напряжений. Подставив в эти уравнения время t = 0 + и
uC (0 + ) и iL (0 + ), определяем зависимые
найденные значения
начальные условия. При составлении уравнений контуры
необходимо выбирать таким образом, чтобы можно было
использовать уже найденные независимые начальные условия.
Пример 1.1.
Для схемы рис. 1.1 определить значения всех токов и
напряжений, а также их производных для моментов времени
t = 0 − , t = 0 + , ∞ , если E = 30 B, R = 10 Ом, L = 10 Гн, C = 100 мкФ.
i3
i1
R
i2
R
E
R
C
uC
L
Рис. 1.1
4
Решение
Рассчитаем схему до коммутации, находящуюся в
установившемся режиме, и определим значения токов и
напряжений. Так как в заданной схеме действует постоянный
источник ЭДС, то все токи и напряжения будут постоянными
величинами. На рис. 1.2 представлена расчетная схема для t < 0 .
i3
i1
R
i2
R
E
R
uC
Рис. 1.2
Имеем:
E
E
= 1 A, i2 (0 − ) =
= 1 A, i3 (0 − ) = 0.
3R
3R
E
u L (0 − ) = 0, u C (0 − ) = i2 (0 − ) R =
R = 10 B,
3R
du
i (0 )
di1
u (0 )
(0 − ) = L − = 0, C (0 − ) = 3 − = 0.
dt
L
dt
C
i1 (0 − ) =
В соответствии с законами коммутации
u C (0 − ) = u C (0 + ) = 10 В , i1 (0 − ) = i1 (0 + ) = 1 А.
Для схемы, полученной после коммутации (сопротивление R
закорочено), составим уравнения по законам Кирхгофа для
мгновенных значений токов и напряжений:
E = R i1 + uC + u L , 0 = uC − R i2 , i1 = i2 + i3 .
Откуда при t = 0 + получим:
E = R i1 (0 + ) + uC (0 + ) + u L (0 + ) ,
0 = u C (0 + ) − R i2 (0 + ) ,
i1 (0 + ) = i2 (0 + ) + i3 (0 + ) .
Следовательно,
5
di1
u (0 )
A
( 0 + ) = L + = 1000 ,
dt
L
c
du C
i (0 )
u (0 )
i 2 ( 0 + ) = C + = 1 A, i3 ( 0 + ) = i1 ( 0 + ) − i 2 ( 0 + ) = 0,
(0 + ) = 3 + = 0.
R
dt
C
u L ( 0 + ) = E − R i1 ( 0 + ) − u C ( 0 + ) = 10 B ,
Примечание.
di2
(0 + ), то
dt
уравнения следует составлять таким образом, чтобы можно было
обойтись без определения вторых производных.
Например, если используется левый контур, то для определения
di2
d 2i1
(0 + ) потребуется найти
(0 + ), продифференцировав
dt
dt 2
di1
di1
di
d 2 i1
+L 2 +R 2.
0=
уравнение E = R i1 + uC + L .
dt
dt
dt
dt
Уравнение 0 = uC − R i2 позволяет легко определить зависимые
начальные условия:
Если требуется определить, например, значение
i (0 )
di2
1 duC
(0 + ) =
(0 + ) = C + = 0.
dt
R dt
dt
При расчете тока i3 (t ) его производная по времени в момент
времени t = 0 + равна:
diC
di
di
A
(0 + ) = 1 (0 + ) − 2 (0 + ) = 1000 .
dt
dt
dt
c
Вынужденную составляющую токов определим при t = ∞ для
схемы (рис. 1.3).
i3
i1
i2
R
E
R
uC
Рис. 1.3
E
E
= 1.5 A, i2вын =
= 1.5 A, iCввы = 0,
2R
2R
= i2 R = 15 B, u Lввы = 0.
i1вын =
uCввы
Пример 1.2.
6
Для схемы рис. 1.4 найти независимые и зависимые начальные
условия,
если
дано:
J = 2 A, E = 50 B, R1 = R2 = R3 = 10 Oм, L = 0.01 Гн, C = 10 мкФ.
i2
R2
C
E
uC
J
R3
i3
L
R1
i1
Рис. 1.4
Решение
Рассчитаем токи и напряжения в схеме до коммутации и
определим независимые начальные условия, учитывая, что
i1 (0 − ) = 0, u L (0 − ) = 0, i3 (0 − ) = 0.
Следовательно,
i2 (0 − ) = J = 2 A = i2 (0 + ), uC (0 − ) = i2 (0 − ) R2 = 20 B = uC (0 + ).
Для схемы после коммутации составим систему уравнений по
законам Кирхгофа для момента времени t = 0 + :
J + i3 (0 + ) = i1 (0 + ) + i 2 (0 + ),
E = i3 (0 + ) R3 + u C (0 + ) + i1 (0 + ) R1 ,
E = i 3 ( 0 + ) R3 + u L ( 0 + ) + i 2 ( 0 + ) R 2 .
Из первого уравнения получим: i3 (0 + ) = i1 (0 + ).
Подставив это соотношение во второе уравнение, определим
значение тока i3 (0 + ) : i3 (0 + ) = 1.5 A, i1 (0 + ) = i3 (0 + ) = 1.5 A.
du C
(0 + ) рассчитаем по соотношению
Значение производной
dt
du C
i ( 0 ) i (0 )
B
(0 + ) = C + = 1 + = 15000 .
dt
C
C
c
Значение u L (0 + ) найдем из третьего уравнения:
u L (0 + ) = E − i3 (0 + ) R3 − i2 (0 + ) R2 = 15 B .
7
1.2 ОПРЕДЕЛЕНИЕ КОРНЕЙ ХАРАКТЕРИСТИЧЕСКОГО
УРАВНЕНИЯ
Расчет переходных процессов классическим методом требует
определения
корней
характеристического
уравнения.
Характеристическое уравнение наиболее просто можно получить
методом входного сопротивления из равенства Z вх ( p) = 0.
Для чего:
1) составляем схему для свободных токов. Для этого в схеме,
полученной после коммутации, все источники энергии заменяем их
внутренними сопротивлениями ( внутреннее сопротивление
идеального источника ЭДС равно нулю, внутреннее сопротивление
идеального источника тока равно бесконечности), а элементы R,L,C
– сопротивлениями, равными соответственно R, pL, 1/pC;
2) если в составленной схеме нет короткозамкнутых ветвей, то
размыкаем любую ветвь, определяем входное сопротивление со
стороны разомкнутой ветви и приравниваем его нулю
Z вх ( p) = 0. Для упрощения алгебраических преобразований следует
размыкать ветвь с наибольшим числом элементов, отдавая при этом
предпочтение ветви с сопротивлением 1/pC.
Если в схеме для свободных токов есть короткозамкнутая ветвь,
то размыкаем ту ветвь, в которой рассчитываем переходный ток. В
цепях с магнитосвязанными индуктивностями для определения
входного сопротивления следует в схеме для свободных токов
предварительно устранить магнитную связь.
Число корней характеристического уравнения равно степени
характеристического уравнения и не может превышать числа
накопителей электромагнитной энергии. Число корней (или порядок
уравнения) можно определить без составления этого уравнения по
упрощенной схеме, которая получается после замены идеальных
последовательно или параллельно соединенных индуктивностей или
емкостей соответственно. Тогда порядок характеристического
уравнения равен числу основных независимых начальных условий
i L (0), uC (0) в послекоммутационной схеме после максимального ее
упрощения (пример 1.4).
Пример 1.3.
Для схем рис. 1.5а, 1.6а, 1.7а составить характеристические
уравнения.
8
i3
i1
i2
R1
R1
R2
R2
C
E
1
Cp
uC
RO
L
Lp
а)
б)
Рис. 1.5
Решение
После коммутации составляем схему для свободных токов
(рис.1.5б). Размыкаем любую ветвь (например, R-Lp) и определяем
входное сопротивление схемы
R ( R + Lp)
1
.
Z 3вх ( p ) =
+ 2 1
Cp R2 + R1 + Lp
Приравняв его нулю, получим характеристическое уравнение
R2 LCp 2 + ( R1 R2 C + L) p + R1 + R2 = 0.
Это же характеристическое уравнение можно получить,
разомкнув другие ветви.
R1
R1
R2
R2
i J (t )
Lp
L
а)
б)
Рис. 1.6
Схеме рис. 1.6а соответствует схема для свободных токов рис.
1.6б и характеристическое уравнение Z вх ( p ) = R1 + R2 + Lp = 0.
Схеме рис. 1.7а соответствует схема для свободных токов рис.
1.7б. В этой схеме средняя и правая ветви замкнуты накоротко, т.е.
схема состоит их двух электрически независимых контуров (рис.
1.7в,г).
9
i1
RO
i2
R1
C
e(t)
L
i
R2
а)
1
Cp
R1
Lp
б)
R1
R2
Lp
в)
1
Cp
R2
г)
Рис. 1.7
Для схемы рис. 1.7в Z1вх ( p ) = R1 + Lp = 0, p1 = −
R1
.
L
1
1
= 0, p2 = −
.
Cp
R2 C
Таким образом, степень характеристического уравнения в
каждом контуре равна 1.
При размыкании левой ветви (рис. 1.7б)
1
( R1 + Lp)( R2 +
)
Cp
= 0.
Z вх ( p) =
1
R1 + Lp + R2 +
Cp
R1
1
.
Из этого уравнения получим два корня p1 = − , p 2 = −
L
Cp
Для схемы рис. 1.7г Z 2вх ( p) = R2 +
Это объясняется тем, что iсв = i1св + i2 св .
Пример 1.4.
Для
схемы
рис.
1.8а
определить
число
корней
характеристического уравнения, не составляя самого уравнения.
10
R1
R1
C2
L3
L3
C1
u (t)
M
R3
L1
u (t)
C3
CЭ
LЭ
R3
L2
а)
б)
Рис. 1.8
Схему рис. 1.8а приведем после коммутации к схеме рис. 1.8б,
в которой LЭ = L1 + L2 ± 2M , C Э = C1 +
C 2 C3
.
C 2 + C3
Так как схема рис. 1.8б имеет три основных независимых
i LЭ (0), i L 3 (0), u CЭ (0),
начальных
условия
то
порядок
характеристического уравнения равен трем.
1.3 ПОРЯДОК РАСЧЕТА ПЕРЕХОДНЫХ ПРОЦЕССОВ
КЛАССИЧЕСКИМ МЕТОДОМ
Перед началом расчета необходимо указать на схеме стрелками
положительные направления токов и напряжений в схеме.
1. Рассчитываем схему до коммутации в установившемся
режиме и определяем независимые начальные условия. Это
единственный этап расчета, в котором используется схема до
коммутации. Все остальные этапы расчета проводятся для схемы
после коммутации.
2. Для t ≥ 0 составляем характеристическое уравнение и
определяем его корни.
3. Записываем уравнение для рассчитываемого тока или
напряжения в виде y (t ) = yвын (t ) + yсв (t ).
Рекомендуется проводить расчет для тока в индуктивности или
напряжения на емкости, для которых известны независимые
начальные условия, так как это упрощает нахождение постоянных
интегрирования. Вид корней характеристического уравнения (см.
п.2) позволяет определить вид свободной составляющей y св (t ).
11
4. Для схемы после коммутации записываем систему
дифференциальных уравнений для мгновенных значений токов и
напряжений. Эта система уравнений позволяет определить
вынужденные составляющие токов или напряжений (в общем
случае) и зависимые начальные условия.
5. Пользуясь системой уравнений п.4, при t = ∞ определяем
y вын (t ) известными методами расчета установившихся режимов.
6. Подставляя в систему уравнений п.4, записанную для
t = 0 + , найденные в п.1 независимые начальные условия,
определяем зависимые начальные условия.
7. Пользуясь начальными условиями, находим постоянные
интегрирования.
8. Записываем y (t ) в окончательном виде и строим график
полученной временной функции.
9. Остальные токи и напряжения целесообразно искать,
пользуясь системой уравнений п.4, причем напряжение на
индуктивности и ток в емкости наиболее просто определяются из
du
di L
, iC = C C .
соотношений u L = L
dt
dt
Пример 1.5.
Для схемы рис. 1.9а определить законы изменения напряжения
на
емкости
и
токи
в
ветвях,
если
дано:
E = 120 В, J = 4 А, R1 = 10Ом, R2 = 30Ом, L = 50 мГн, C = 500 / 3 мкФ.
а
C
E
i1
uC
L
i1
J
R2
uC
C
E
J
uab
i2
R2
R1
R1
в
а)
б)
Рис. 1.9
Решение
1.Рассматриваем схему до коммутации (рис. 1.9б) и определяем
ток в ветви с индуктивностью и напряжение на емкости для момента
времени t = 0 − из уравнений:
12
uC (0 − ) = E = 120 В,
u ав (0 − ) =
Eg 2 + J
= 60 В,
g1 + g 2
u ав (0 − )
= 6 А.
R1
В соответствии с законами коммутации получим:
uC (0 − ) = 120 В = uC (0 + ) , i1 (0 − ) = 6 А = i1 (0 + ).
2. Пользуясь схемой после коммутации для свободных токов
(рис. 1.9в), составляем характеристическое уравнение и находим его
корни:
1
Lp + R1 + R2 +
= 0,
Cp
i1 (0 − ) =
LCp 2 + ( R1 + R2 ) Cp + 1 = 0,
p1 = −200
1
1
, p 2 = −600 .
c
c
1
Cp
C
i1
Lp
uC
L
J
R2
R2
R1
i2
R1
в)
г)
Рис. 1.9
3. Выбираем в качестве искомой функции напряжение на
емкости (можно выбрать и ток в индуктивности):
uC = uСвын + uСсв = uСвын + A1A −200 t + A2 A −600 t .
4. Для схемы после коммутации (рис. 1.9г) записываем систему
независимых уравнений для мгновенных значений токов и
напряжений:
J = i1 + i2 ,
L
di1
+ R1i1 − R2 i2 − uC = 0.
dt
5. В установившемся режиме ( t = ∞ ) ток i2вын = 0 , так как
схема питается от источника постоянного тока. Следовательно,
i1вын = J = 4 A, u L вын = L
di1вын
= 0, uСвын = i1вын R1 = 40 B.
dt
13
6. При t = 0 + имеем: J = i1 (0 + ) − i2 (0 + ) . Откуда определяем
du C
i (0 ) i (0 )
B
(0 + ) = C + = 2 + = −1200 .
i2 (0 + ) = J − i1 (0 + ) = −2 A и
dt
C
C
c
Примечание. Если расчет ведем для тока
i1 , то значение
di1
(0 + ) находим из уравнения
dt
di1
(0 + ) + R1i1 (0 + ) − R2 i2 (0 + ) − u C (0 + ) = 0 ,
dt
R i (0 ) + u C (0 + ) − R1i1 (0 + )
di1
(0 + ) = 2 2 +
= 0.
dt
L
L
7. Постоянные
уравнений
интегрирования
определяем
из
системы
uC = uСвын + uСсв = 40 + A1A −200t + A2 A −600t ,
duC duСвын duСсв
=
+
= −200 A1A −200t − 600 A2 A −600t ,
dt
dt
dt
которая при t = 0 + принимает вид:
uC (0 + ) = 120 = 40 + A1 + A2 ,
du C
(0 + ) = −12000 = −200 A1 − 600 A2 .
dt
Откуда
находим
значения
постоянных
интегрирования:
A1 = 90, A2 = −10.
8. Окончательно получим закон изменения напряжения на
емкости
uC (t ) = 40 + 90A −200 t − 10A −600 t .
9. Покажем, как определять остальные переменные:
du C
= −3A −200t + A −600t ,
dt
i1 = J − i2 = 4 + 3A −200t − A −600t ,
i 2 = iC = C
uL = L
di1
= −30A − 200t + 30A −600t .
dt
Пример 1.6.
Для схемы рис. 1.10 определить закон изменения напряжения
на
емкости,
если
дано:
E = 120 В, R1 = R2 = R0 = 10Ом, L = 0.1Гн, C = 100 мкФ.
14
R2
RO
E
i1
i3
C
i2
uC
L
R1
Рис. 1.10
Решение
1. Так как цепь подключена к источнику постоянного
напряжения, то в установившемся режиме до коммутации
индуктивность имеет нулевое сопротивление, а емкость –
бесконечно большое. Поэтому
i3 (0 − ) = 0 , u L (0 − ) = 0 , i1 (0 − ) =
i2 (0 − ) =
E
= 4 A,
R0 + R1 + R2
E
= 4 A , u C (0 − ) = i2 (0 − ) R2 = 40 B .
R0 + R1 + R2
В
соответствии с
законами
коммутации
имеем
i2 (0 − ) = 4 A = i2 (0 + ) , u C (0 − ) = 40 B = u C (0 + ) .
2.Корни
характеристического
уравнения
Z вх ( p ) =
R ( R + Lp )
1
= 0 равны p1, 2 = −δ ± jωсв = −100 ± j100.
+ 1 2
Cp R1 + R2 + Lp
При комплексно-сопряженных корнях характеристического
уравнения решение ищем в виде:
u C (t ) = u Свын + ( A1 cos ω св t + A2 sin ω св t )A −δ t .
3. Используя уравнения для схемы после коммутации,
определим требуемое для нахождения постоянных интегрирования
duC
(0 + ) :
dt
E = i1 (0 + ) R1 + u C (0 + ) → i1 (0 + ) = 8 A ,
значение производной
i1 (0 + ) = i2 (0 + ) + i3 (0 + ) → i3 (0 + ) = 4 A ,
i ( 0 ) i (0 )
du C
A
(0 + ) = C + = 3 + = 4000 .
dt
C
C
c
t=
4. В установившемся режиме при
составляющая напряжения на емкости равна
E
u Свын = i2 вын R2 =
R2 = 60 B .
R1 + R 2
∞ вынужденная
15
5.
Определим
постоянные
интегрирования, используя
duC
(0 + ) . Для чего в систему
найденные начальные условия uC (0 + ) ,
dt
уравнений
u C = 60 + ( A1 cos100t + A2 sin 100t )A −100t ,
du C
= ( A1100 sin 100t + A2 100 cos100t )A −100t +
dt
+ ( A1 cos100t + A2 sin 100t )(−100)A −100t .
для t = 0 + подставим их значения:
u C (0 + ) = 40 = 60 + A1,
du C
(0 + ) = 4000 = A2 ⋅ 100 − A1 ⋅ 100 .
dt
Откуда: A1 = −20 , A2 = 20 .
Ответ:
u C (t ) = 60 + (20 sin 100t − 20 cos 100t )A −100t =
= 60 + 20 2 sin(100t −
π
4
)A −100t .
Пример 1.7.
Рассчитать все токи в схеме рис. 1.11а, если дано:
E = 120 В, R1 = R2 = R3 = 10 Ом, L1 = L2 = 0.2 Гн, M = 0.1 Гн .
R1
R1
R3
R2
E
i2
i1
i3
(L2 + M) p
M
L2
R3
R2
L3
( L2 + M ) p
− Mp
а)
б)
Рис. 1.11
Решение
1. Определим независимые начальные условия из схемы до
коммутации i2 (0 − ) =
E
= 6 A = i2 (0 + ), i3 (0 − ) = 0 = i3 (0 + ) .
R1 + R2
16
2.
Методом
входного
сопротивления
составляем
характеристическое
уравнение,
устранив
предварительно
магнитную связь (рис. 1.11б):
[R + ( L2 + M ) p ][R3 + ( L3 + M ) p ] = 0 ,
Z ( p ) = R1 − Mp + 2
R2 + ( L2 + M ) p + R3 + ( L3 + M ) p
100 1
1
0.03 p 2 + 10 p + 300 = 0 , p1 = −
, p 2 = −300 .
3 c
c
3. Выражение для тока i2 запишем в виде
i2 = i2вын + i2 св = i2 вын + A1A
−
100
t
3
+ A2 A − 300 t .
4. Так как R2 = R3 , то вынужденная составляющая тока равна
E
i2вын =
= 4 A.
R2 R3
2( R1 +
)
R2 + R3
5. Составим систему уравнений по законам Кирхгофа для
послекоммутационной схемы:
di
di
i1 R1 + i2 R2 + L2 2 − M 3 = E ,
dt
dt
di
di
di
di
i2 R2 + L2 2 − M 3 − i3 R3 − L3 3 + M 2 = 0,
dt
dt
dt
dt
i1 = i2 + i3 .
При t = 0 + имеем:
di
di2
(0 + ) − M 3 (0 + ) = E ,
dt
dt
di
di
di
di
i2 (0 + ) R2 + L2 2 (0 + ) − M 3 (0 + ) − i3 (0 + ) R3 − L3 3 (0 + ) + M 2 (0 + ) = 0,
dt
dt
dt
dt
i1 (0 + ) = i2 (0 + ) + i3 (0 + ).
i1 (0 + ) R1 + i2 (0 + ) R2 + L2
Учитывая, что i3 (0 + ) = 0, i2 (0 + ) = 6 A, получим:
i1 (0 + ) = 6 A,
di2
A di3
A
(0 + ) = 200 ,
(0 + ) = 400 .
dt
c dt
c
6. Постоянные интегрирования определяются из следующей
системы уравнений:
i2 = 4 + A1A
−
100
⋅t
3
−
di2
100
=−
A1A
dt
3
+ A2 A −300 t ,
100
t
3
− 300 A2 A −300 t .
17
При t = 0 +
i2 (0 + ) = 6 = 4 + A1 + A2 ,
di2
100
(0 + ) = 200 = −
A1 − 300 A2 .
dt
3
Откуда А1=3, А2=-1.
7. Окончательно получим
i2 = 4 + 3A
8. Ток
100
t
3
− 1A − 300 t .
i3 рассчитываем аналогично:
i3 = 4 − 3A
9. Ток
−
−
100
t
3
− 1A − 300 t .
i1 = i2 + i3 .
Примечание.
В зависимости от вида корней характеристического уравнения
функция y(t) имеет различный вид.
−100⋅t
+ 4 ⋅ A −200⋅t .
Пусть i(t) изменяется по закону: i (t ) = 5 − 3 ⋅ A
Построим его график по составляющим i(t) (рис. 1.12), причем
по расчету i(0-)=2 А.
i
6
4
τ2
2
t
−2
−4
τ1
Рис. 1.12
При построении графика учтено, что: 1) при t=0+ ток i(0+)=6 А
(ток в момент коммутации изменился скачком от 2 А до 6 А); вторая
экспонента затухает в два раза быстрее, чем первая; 3)
di
(0 + ) = −500 A / c. Это значение и его знак определяют поведение
dt
тока в момент времени t=0+.
18
−100 t
.
Пусть ток iL(t) изменяется по закону: iL (t ) = 5 sin 400 t ⋅ A
График тока представлен на рис. 1.13.
i
T
5 ⋅ A −100⋅t
t
− 5 ⋅ A −100⋅t
τ
Рис. 1.13
При построении графика учтено, что: 1) iL(0+)=0; 2) по
di L
A
(0 + ) = +2500 : 3) период колебаний равен
результатам расчета
dt
c
2π
T=
= 0,0157 с − 15,7 мс, постоянная
времени цепи равна
ω св
1
1
с = 0.01мс.
τ= =
δ 100
Раздел 2. РАСЧЕТ ПЕРЕХОДНЫХ ПРОЦЕССОВ ПРИ
ПРОИЗВОЛЬНЫХ ВХОДНЫХ ВОЗДЕЙСТВИЯХ С
ИСПОЛЬЗОВАНИЕМ ИНТЕГРАЛОВ НАЛОЖЕНИЯ
Этот метод основан на сочетании расчета переходных
процессов при постоянных источниках и принципа суперпозиции.
При этом предполагается, что в схеме действует только один
внешний источник сигнала.
2.1. ОПРЕДЕЛЕНИЕ ПЕРЕХОДНОЙ И ИМПУЛЬСНОЙ
ХАРАКТЕРИСТИК ЦЕПИ
Переходная характеристика цепи h(t) – это отношение реакции
этой цепи при нулевых начальных условиях на ступенчатое
воздействие к величине воздействия, при этом в цепи отсутствуют
внутренние источники энергии.
19
Обозначим входное воздействие f1(t), а выходное - f2(t). Тогда
переходная характеристика цепи равна
f (t )
h(t ) = 2 .
f1
Так как обычно входное воздействие является единичным
ступенчатым воздействием, то переходная характеристика цепи h(t)
численно равна реакции цепи на единичное ступенчатое
воздействие при нулевых начальных условиях.
f (t )
h(t ) = 2 .
1
Размерность h(t) равна отношению размерности отклика к
размерности воздействия, поэтому переходная характеристика
может иметь размерность сопротивления, проводимости или быть
безразмерной.
Например, если откликом является ток, а воздействием ЭДС
(1В), то переходная характеристика является переходной
проводимостью и имеет размерность Сименс.
Примечание.
Условимся
нижний
индекс
переходной
характеристики цепи выбирать соответствующим искомому
отклику. Так, если требуется рассчитать, например, напряжение на
емкости, то переходная характеристика обозначается huC (t ) ,а если,
например, ток i3(t), то hi3 (t ) .
Импульсная характеристика k(t) численно равна реакции цепи
на единичное импульсное воздействие.
Обычно при расчете классическим методом импульсную
характеристику определяют как производную от переходной
/
характеристики k (t ) = h (t ).
Пример 2.1.
Для схемы рис. 2.1 рассчитать переходную характеристику для
напряжения на емкости и тока в неразветвленной части схемы i1,
если на входе схемы действует источник напряжения , а параметры
схемы следующие: R1 = 100 Ом, R2 = 400 Ом, C = 125 мкФ.
i1
R1
e(t)
R2
C
uC
Рис. 2.1
20
Решение
Полагая, что схема подключается к источнику напряжения 1В
(рис. 2.2), рассчитываем u C и i1 классическим методом:
uC (t ) = uСвын + uСсв = uСвын + A ⋅ A p t .
t ≥0
i1
R1
R2
E = 1B
C
uC
Рис. 2.2
Для определения значения р составляем характеристическое
уравнение
1
RR
1
Z вх ( p) =
+ 1 2 = 0, откуда p = −100 .
Cp R1 + R2
c
Вынужденная составляющая напряжения
иС
равна
E = 1B
u Свын =
⋅ R2 = 0.8 B .
R1 + R2
Учитывая, что u C (0 − ) = 0 = u C (0 + ) , определим постоянную
интегрирования А:
u C ( 0 + ) = 0 = 0.8 + А → А = − 0.8.
Следовательно,
uC (t ) = 0.8 − 0.8 ⋅ A −100 t .
Ток i1 можно найти аналогично, можно из уравнения
E = i1 R1 + uC → i1 (t ) = 0.02 + 0.08 A −100 t .
Переходные характеристики для напряжения на емкости и тока
i1 равны соответственно
huC (t ) = u uC (t ) ⋅1(t ) = 0.8 − 0.8 ⋅ A −100 t ⋅1(t ) ,
[
[
hi1 (t ) = ii1 (t ) ⋅1(t ) = 0.02 + 0.08 ⋅ A −100 t
]
]⋅1(t ) .
Причем, переходная характеристика для напряжения на
емкости является безразмерной функцией, а переходная
характеристика для тока имеет размерность Сименс.
В выражениях для переходных и импульсных характеристик
единичная функция играет роль ключа.
21
Пример 2.2.
Для схемы (рис. 2.3) определить переходную и импульсную
характеристики цепи по напряжению при R=10 Ом и L=0.1 Гн.
i
R
u (t)
uL
L
Рис. 2.3
Решение
Переходную характеристику рассчитываем, полагая, что схема
подключается к источнику постоянного напряжения 1В, причем
начальные условия нулевые, т.е. i L (0) = 0 .
R
i L (t ) = i L вын + i L св
− t
1
= + A ⋅ A L = 0.1 + A ⋅ A −100 t .
R
При t = 0 + i L (0 + ) = 0 = 0.1 + A. Откуда А=-0.1.
Следовательно,
di
iL (t ) = 0.1 − 0.1 ⋅ A −100 t , u L (t ) = L = 1 ⋅ A −100 t , hu L (t ) = 1 ⋅ A −100 t ⋅ 1(t ) .
dt
Импульсную характеристику просто
определить путем
дифференцирования переходной характеристики.
2.2. ПРИМЕНЕНИЕ ИНТЕГРАЛОВ НАЛОЖЕНИЯ ДЛЯ
РАСЧЕТА ПЕРЕХОДНЫХ ПРОЦЕССОВ
В зависимости от того, какая характеристика цепи (переходная
или импульсная) используется для расчета переходного процесса
при произвольных воздействиях, интегралы наложения имеют
различные формы записи. Одной из форм интегралов наложения с
использованием переходной характеристики, получившей большое
распространение в электротехнике, является интеграл Дюамеля.
Обозначили входное воздействие через f1(t), а выходной сигнал
(отклик) через f2(t). Для случая, когда f1(t) является непрерывной
функцией при всех t, за исключением точки t=0 , где f1(t) может
иметь разрыв первого рода, т.е. f1 (0) ≠ 0 , интеграл Дюамеля имеет
вид:
22
t
t
f 2 (t ) = f1 (0)h(t ) + ∫ f (τ )h(t − τ )dτ = f1 (0)h(t ) + ∫ f1/ (t − τ )h(τ )dτ =
/
1
0
0
t
t
= f1 (t )h(0) + ∫ f1 (τ )h (t − τ )dτ = f 1 (t )h(0) + ∫ f 1 (t − τ )h / (τ )dτ .
/
0
0
С помощью интеграла Дюамеля можно определить реакцию
цепи на входное воздействие и в случае, если входное воздействие
описывается кусочно-непрерывной функцией, имеющей конечное
число разрывов первого рода. В этом случае интервал
интегрирования следует разбить на несколько интервалов в
соответствии с интервалами непрерывности воздействия f1(t) и
учесть при этом реакцию цепи на конечные скачки воздействия f1(t)
в точках разрыва.
f1 (t )
f11
t1
t
f12
Рис. 2.4
Например, если воздействие f1(t) имеет вид рис. 2.4, то решение
записывают для каждого интервала:
0 ≤ t < t1
t
f 2 (t ) = ∫ f11/ (τ )h(t − τ )dτ ,
0
t ≥ t1
t1
t
f 2 (t ) = ∫ f (τ )h(t − τ )dτ + [ f11 (t1− ) − f12 (t1+ )]h(t − t1 ) + ∫ f12/ (τ )h(t − τ )dτ .
0
/
11
t1
Примечание. Так как расчет переходных процессов с
использованием интеграла Дюамеля проводится по интервалам, то в
выражении переходной характеристики цепи множитель 1(t)
отсутствует.
Порядок расчета
1. Определяем переходные и (или) импульсные характеристики
цепи.
2. Проводим расчет, используя интеграл Дюамеля.
23
Пример 2.3.
Для схемы рис. 2.5 определить закон изменения тока в цепи,
если входное напряжение изменяется по закону рис. 2.6. При этом
U = 10 B, t1 = 0.001c, R = 10 Ом, C = 100 мкФ.
i
u (t)
R
u (t)
C
U
uC
t1
Рис. 2.5
t
Рис. 2.6
Решение
1. Определяем переходную характеристику hi (t ) . Для чего
рассчитываем классическим методом напряжение на емкости,
считая, что схема с нулевыми начальными условиями (речь идет о
независимых начальных условиях, т.е. в данном случае о
напряжении на емкости) подключается к источнику постоянного
напряжения 1 В.
В результате расчета получим
u C (t ) = 1 − 1 ⋅ A −1000 t .
−1000 t
i(t ) = 0.1 ⋅ A
.
Ток в схеме равен
Следовательно, ему соответствует переходная характеристика
hi (t ) = 0.1 ⋅ A −1000 t .
2. Так как входное воздействие является разрывной функцией,
то решение проводим по интервалам:
0 ≤ t < t1
t
t
i (t ) = u (0)hi (t ) + ∫ u (τ )hi (t − τ )dτ = ∫ k ⋅ 0.1 ⋅ A −1000 ( t −τ ) dτ =
/
0
=
0
3
10
⋅ A −1000 t ⋅ A1000τ
3
10
t
0
= 1 − A −1000 t .
При интегрировании учтено, что
входное напряжение
изменяется по линейному закону u (t ) = k t , где k = 1000.
24
t ≥ t1
t1
i (t ) = ∫ u / (τ )hi (t − τ )dτ + Δu hi (t − t1 ) = A −1000 t ⋅ A1000τ
t1
0
− 10 ⋅ 0.1 ⋅ A −1000 (t −t1 ) =
0
= A −1000 t ⋅ (A1 − 1) − A −1000 t ⋅ A1 = −A −1000 t .
На рис .2.7 представлен график рассчитанного отклика.
i
0.63
t1
t
− 0.37
Рис. 2.7
Примечание. Следует обратить внимание на то, что напряжение
на емкости не может менять скачком при любом изменении
входного воздействия. В то время как ток в емкости может меняться
скачком, если скачком меняется входное воздействие, причем,
величина скачка тока определяется величиной скачка входного
напряжения так же, как и момент времени, в который этот скачок
происходит. В рассматриваемом примере входное напряжение
уменьшается на 10 В, следовательно, скачок тока равен –1 А. Это же
значение получим, рассчитав i (t1+ ) − i(t1− ) = −0.37 − 0.63 = −1 A. .
Пример 2.4.
Схема рис. 2.3 подключается к напряжению рис. 2.8. Рассчитать
закон изменения тока в схеме и напряжения на индуктивности, если
U = 10 B, t1 = 0.01с, R = 10 Ом, L = 0.1 Гн.
u (t)
U
t1 2t1
t
−U
Рис. 2.8
25
Решение
1. Переходные характеристики были рассчитаны в примере 2.2,
−100 t
−100 t
.
они равны hi (t ) = 0.1 − 0.1 ⋅ A
, huL (t ) = 1 ⋅ A
2. Так как входное воздействие является разрывной функцией,
то решение проводим для каждого интервала:
0 ≤ t < t1
t
u L (t ) = u (0)huL (t ) + ∫ u / (τ )huL (t − τ )dτ =10 ⋅ A −100 t ,
0
i (t ) = 1 − 1 ⋅ A −100 t .
При интегрировании учтено, что непрерывная часть входного
/
напряжения u(t ) = U , u (t ) = 0.
t1 ≤ t < t2
/
Так как и на втором и третьем интервалах u (t ) = 0 , то
опустим в записи интеграла Дюамеля соответствующие слагаемые.
u L (t ) = u (0)huL (t ) + Δu (t1 ) huL (t − t1 ) = 10 ⋅ A −100 t − 20 ⋅ A −100(t −t1 ) =
= −44 ⋅ A −100 t ,
i (t ) = −1 + 4.4 ⋅ A −100 t .
t ≥ t2
u L (t ) = u (0)huL (t ) + Δu (t1 ) huL (t − t1 ) + Δu (t 2 )huL (t − t 2 ) =
= 10 ⋅ A −100 t − 20 ⋅ A −100(t −t1 ) + 10 ⋅ A −100 (t −t2 ) =
= 30 ⋅ A −100 t ,
i (t ) = −3 ⋅ A −100 t .
На рис. 2.9 представлены графики рассчитанных откликов.
i
uL
10
0.63
3.7
t1
4.05
2t1
− 5.95
t
2t1
t1
t
− 0.405
Рис. 2.9
− 16.3
26
Проверка правильности расчета проводится так же, как в
предыдущем примере. Необходимо только учесть, что в этой схеме
скачком не может меняться ток в индуктивности, а напряжение на
индуктивности может.
КОНТРОЛЬНЫЕ ВОПРОСЫ
1. Что такое переходной процесс?
2. В чём суть классического метода расчёта переходных процессов?
3. Как определяются вынужденная и свободная составляющие искомой
функции?
4. От чего зависит вид вынужденной и свободной составляющих искомой
функции?
5. Что такое независимые и зависимые начальные условия и как они
определяются? Для какого момента времени и для какой схемы?
6. Как составляется характеристическое уравнение методом входного
сопротивления? От чего зависит порядок характеристического
уравнения?
7. Как определяется длительность переходного процесса?
8. Что такое переходная и импульсная характеристики цепи?
9. Как рассчитываются переходные и импульсные характеристики? При
каких начальных условиях? Какую размерность они имеют?
10. Как связаны между собой переходная и импульсная характеристики?
11. В чём заключается особенность расчёта переходных процессов при
произвольных входных воздействиях методом интегралов Дюамеля?
12. Как можно проверить правильность расчёта методом интегралов
Дюамеля?
СПИСОК ЛИТЕРАТУРЫ
1. Попов В.П. Основы теории цепей. – М.: Высшая школа,
2000. 576с.
2. Кузовкин В.А. Теоретическая электротехника. – М.: Логос,
2002. 480с.
3. Матханов П.Н. Основы анализа электрических цепей.
Линейные цепи. – М.: Высшая школа, 1990. 364с.
27
Download