Прикладные задачи теоретической электротехники. Учебное пособие
.pdfИсходное уравнение Дуффинга: |
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
d 2q |
|
|
|
R |
dq |
|
1 |
(aq bq3) |
|
R |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1 |
|||||
|
|
|
|
|
|
|
|
|
dt2 |
L1 dt |
L1 |
|
L1 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
где q – заряд; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Im – амплитуда тока источника тока; |
|
|||||||||||||||||||||||
, |
|
R1, |
L1 – описаны ранее; |
|
|
|
|
|
|
|||||||||||||||
R1 : 100 – сопротивление, Ом; |
|
|
|
|
|
|||||||||||||||||||
L : 100 10 3 |
– индуктивность, Гн; |
|
||||||||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f0 : 50 – частота, Гц; |
|
|
|
|
|
|
|
|
||||||||||||||||
: 2 f0 ; |
– верхний предел по времени; |
|||||||||||||||||||||||
tmax : 0,5 |
||||||||||||||||||||||||
tmin : 0,05 – нижний предел по времени; |
||||||||||||||||||||||||
Im0 : 0, 22 |
– максимальный ток, А; |
|
||||||||||||||||||||||
i1 : 0, 22, |
i2 : 0, 2; |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
q : 6 10 4 , q |
2 |
|
: 6,5 10 4 ; |
|
|
|
|
|
||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
U С1 : 8,9, |
UС2 : 9; |
|
|
|
|
|
|
|
|
|||||||||||||||
С |
|
|
: |
|
|
q1 |
6,742 10 5; |
|
C |
2 |
: |
q2 |
|
|||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
1 |
|
|
UС1 |
|
|
|
|
|
|
|
|
|
|
|
UС2 |
|
||||||||
|
|
|
|
|
|
|
|
) 7, 222 10 5 ; |
|
|
||||||||||||||
C |
0 |
: max(C , C |
2 |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
a := 1 |
|
|
b := 1. |
|
|
|
|
|
|
|
|
|
|
|||||||||||
Given |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
U |
С1 |
aq bq3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
UС2 aq2 |
bq23 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
a |
|
|
|
|
|
|
|
|
|
|
2,052 104 |
|
|
|
|
|||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
: Find (a,b) |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
b1 |
|
|
|
|
|
|
|
|
|
1,579 10 |
|
|
|
|
|
|||||||||
Im sin( t) ,
7,222 10 5;
31
RES : |
1 |
372,104 |
|
L1C0 |
|||
|
|
T : 2
TRES : 2 RES
t : |
min(T ,TRES ) |
|
|
|
|
|||
|
|
100 |
|
|
|
|
|
|
points : |
t |
max |
t |
min |
|
|||
floor |
|
|
|
|||||
|
|
t |
|
|||||
|
|
|
|
|
|
|
||
K : 0,9; |
Im0 0, 22. |
|
||||||
Решение дифференциального уравнения методом Адамса:
f 1(t) : R1 Im0K sin( t) ;
L1
|
|
|
|
|
|
y11 |
|
|
|
|
|
|
|
R |
|
1 |
|
|
|
|
; |
D01( y, y1) : |
|
|
a y1 b ( y1 )3 |
|
|
|
||||
|
|
|
1 |
y1 |
|
f (t) |
|
|
||
|
L |
L |
|
|||||||
|
|
1 |
1 0 1 0 |
|
1 |
|
||||
|
|
|
1 |
|
1 |
|
|
|
|
|
tol : 10 5; |
|
|
|
|
|
|
|
|
|
|
y10 : 0,001; |
|
|
|
y11 : 0,001; |
|
|
|
|
||
Z1: Adams( y1, tmin , tmax , points, D01, tol); |
|
|
|
|||||||
t : Z1 0 , |
yL2 : Z1 2 , |
yL1: Z1 1 ; |
|
|
|
|
||||
K : 0... ROWS(Z1) 1.
Вывод графиков аналогичен описанному в примере 2.3.1.
32
Вольт-кулонная характеристика и ее таблица (рис. 2.7 и таблица)
Вольт-кулонная характеристика
q |
Кл 10 4 |
0 |
0,5 |
1,0 |
1,5 |
2,0 |
2,5 |
3,0 |
3,5 |
4,0 |
4,5 |
5,0 |
5,5 |
6,0 |
6,5 |
uC |
В |
0 |
1,5 |
2,8 |
4,0 |
5,0 |
5,9 |
6,6 |
7,2 |
7,8 |
8,2 |
8,5 |
8,7 |
8,9 |
9,0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
А |
2,0 |
1,7 |
1,44 |
1,2 |
1,0 |
0,82 |
0,68 |
0,56 |
0,44 |
0,36 |
0,3 |
0,26 |
0,22 |
0,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
uс, В
10
9
8
7
6
5
4
3
2
1 q, Кл·10–4
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7
Рис. 2.7
Выбранные значения uс и q в дальнейшем используются для рас-
четаемкостиконденсатора, которая, всвоюочередь, служитдлярасчета коэффициентов аппроксимирующего полинома.
33
Контрольные вопросы и задания
1.Перечислить основные численные методы анализа нелинейных цепей.
2.Перечислить основные функции, предназначенные для численного решения нелинейных уравнений в программе Mathcad.
3.Перечислить основные типы нелинейных дифференциальных уравнений.
3.Вариационный метод анализа линейных
инелинейных электрических цепей
Вариационное исчисление
Вариационное исчисление [1] ставит перед собой задачу отыскания функций x, y, z, … от s, t, …, при которых интеграл
s2 t2
S ∫ ∫ V s, t, , x, y, , xs , xt , , ds dt
s1 t1
x x ,t t
имеет минимальное или максимальное значение (экстремум). Пределы интегрирования при этом могут предполагаться фиксированными или жепеременнымиприопределенныхограничивающихусловиях. Случай переменных пределов интегрирования, однако, приводится к постоянным пределам. Функции из рассматриваемого класса, дающие решение задачи, называются экстремалями.
Вариационные задачи могут решаться прямыми методами (например, с помощью аппроксимации) или косвенным путем – с помощью приведения к дифференциальным уравнениям. В то же время часто может оказаться полезным обратный путь – приведение дифференциального уравнения к вариационной задаче с последующим ее решением прямыми методами вариационного исчисления.
34
Приведение к дифференциальным уравнениям
Уравнения Эйлера – Лагранжа
Далее приводятся наиболее употребимые случаи.
а) V (x, x, t) содержит одну функцию x(t) одной независимой пере-
менной t и ее производную x(t) dxdt . Если будем варьировать функ-
цию x(t), заменяя ее выражением
x(t) (t) x(t) x,
то соответствующая вариация S определится формулой
|
S(x ) |
t2 |
V |
|
V |
|
S |
|
∫ |
|
|
dt, |
|
|
|
x |
||||
|
t |
x |
|
|
||
|
|
1 |
|
|
|
|
откуда, интегрируя по частям член, содержащий , будем иметь:
dV t2 |
t2 |
V |
|
d V |
|||||
S |
|
|
∫ |
|
x |
|
|
|
dt. |
|
|
||||||||
dx |
t |
t |
|
|
dt |
x |
|||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
Если функция x(t) задана для значений t1 и t2 , то в этих точках 0
и, следовательно, первоеслагаемоевыпадает. Изнеобходимогоусловия экстремума S, согласно которому для всех допустимых функций
должнобыть S 0, следуетдифференциальноеуравнениеЭйлера – Ла- гранжа:
V |
— |
d V |
0. |
||
x |
|
|
|
||
|
|||||
|
dt |
x |
|
||
Оно определяет x(t) с точностью до двух произвольных постоянных, значения которых находятся с помощью краевых условий.
35
б) V (x, y, x, y, |
, t) содержит несколько n функций одной незави- |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
симой переменной t. Здесь будем иметь |
S 0 , |
если функции x(t), |
|||||||||
y(t) и другие удовлетворяют уравнениям |
|
|
|
|
|
|
|||||
V |
|
d V |
V |
|
d |
V |
|
|
|||
x |
|
|
|
0, |
y |
|
|
|
|
|
0, |
|
|
||||||||||
|
dt |
x |
|
dt |
y |
|
|
||||
Это n дифференциальных уравнений второго порядка. Они определяют x(t), y(t), … с точностью до 2n постоянных интегрирования, которые находятся с помощью краевых условий.
в) V (x, x, x, x , t) содержит одну функцию и ее производные до третьего порядка включительно. В этом случае получим:
V |
|
|
d V |
|
|
|
d 2 V |
|
|
V |
|
|
d |
V |
|
|
V t2 |
|||||||||||||||||||||||
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
x |
dt x |
|
dt |
|
|
x |
|
|
x |
dt |
x |
|
x |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
t2 |
|
|
V |
|
|
d |
|
|
V |
|
|
|
d 2 |
|
V |
|
|
|
d3 V |
|
|
|
|
||||||||||||
|
|
|
∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
3 |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
t1 |
|
|
|
|
dt x |
|
|
|
|
|
|
x |
|
|
|
|
|
x |
|
|
|
|
|||||||||||||
Таким образом, |
S будет равно нулю, |
|
если при t1 и t2 |
заданы, кроме |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
и удовлетворяется уравнение |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
x(t), также x и x, |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
V |
|
d |
|
V |
|
|
|
d 2 |
|
V |
|
|
|
d3 |
|
V |
0. |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
x |
|
|
|
|
|
|
dt |
2 |
|
|
|
|
dt |
3 |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
dt x |
|
|
|
|
|
|
x |
|
|
|
|
|
|
x |
|
|
|
|
|
|
||||||||||||||
г) V (s, t, |
, |
x, |
y, |
, |
xs , |
xt , |
|
|
, ys |
, yt , ) содержит несколько не- |
||||||||||||||||||||||||||||||
зависимых переменных s, t, …, несколько зависимых переменных x, y, … и их производные xs , ys , Вариация S будет равна нулю, если
x(s, t, …) и другие переменные определяются из уравнений Лагранжа:
V |
|
|
|
V |
|
|
|
V |
0 |
, |
||
x |
|
|
|
|
|
|
|
x |
||||
s |
x |
|
t |
|||||||||
|
|
|
|
|
s |
|
|
|
t |
|
|
|
36
V |
|
|
|
V |
|
|
|
V |
0 |
||
y |
|
|
|
|
|
|
|
y |
|||
s |
y |
|
t |
||||||||
|
|
|
|
|
s |
|
|
|
t |
|
|
и т. д.*
3.1. Линейные уравнения
Рассмотрим последовательную линейную цепь, состоящую из резистора сопротивлением R, катушки индуктивности с индуктивностью L и конденсатора емкостью C.
Найдем свободную составляющую тока i: uL uR uC U 0 ,
где U0 – напряжение, до которого предварительно заряжен конденсатор; uR , uL и uC – напряжения на элементах R, L и C.
iR L dtdi C1 ∫idt U0 .
Продифференцируем это уравнение один раз по времени:
L d 2i R di 1 i 0. dt2 dt C
Приведем это уравнение к следующему виду:
d 2i |
|
R di |
|
1 |
i 0. |
|||
|
|
|
|
|
||||
dt2 |
L dt |
LC |
||||||
|
|
|
||||||
Рассмотрим характеристическое уравнение p2 RL p LC1 0.
* В предположении, что на границе области изменения независимых переменных s, t, … функции x, y, … принимают заданные значения.
37
Рассмотрим такую же цепь
uR uL uC U0 ,
|
|
|
iR L di |
|
1 |
|
∫idt u0 . |
||||||||||||||||||
|
C |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|||||
Зная, что i C |
duC |
, получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
du |
C |
|
|
|
|
|
|
|
d 2u |
|
|
|
|||||||||
|
|
CR |
|
|
|
LC |
|
|
|
|
|
C |
|
u |
0, |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
dt |
|
|
|
|
|
|
|
dt2 |
|
C |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
d 2u |
|
|
|
|
R du |
|
|
|
|
|
1 |
|
|
|
||||||||
|
|
|
|
C |
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
u |
0. |
|||||
|
|
|
dt2 |
|
L |
dt |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
LC C |
|
|||||||||||||||
Окончательно будем иметь |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
p2 |
R |
p |
|
|
1 |
|
|
0. |
|
|||||||||||
|
|
|
|
|
L |
|
LC |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Составим функцию энергий для рассматриваемой выше цепи:
CUC2 |
|
L |
i2 |
R∫i2dt |
CU02 |
. |
|
2 |
2 |
2 |
|||||
|
|
|
|
Эта цепь последовательная, а значит
iL iC C dudtC ,
откуда получим функцию F:
CU 2 |
L |
du |
2 |
|
du |
2 |
CU 2 |
|
||||
C |
|
|
C |
C |
|
R∫ C |
C |
|
dt |
0 |
, |
|
2 |
2 |
dt |
dt |
2 |
||||||||
|
|
|
|
|
|
|
||||||
38
|
|
|
|
|
|
|
|
|
|
|
dF |
CUC , |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
duC |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
dF |
|
|
|
|
|
d |
|
2 |
du |
|
|
|
|
2 |
du |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
R∫ |
|
|
C |
|
|||
dt |
du |
|
|
dt |
LC |
|
dt |
2C |
|
dt |
dt |
|||||||||||||
|
|
d |
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
LC2 |
|
d 2u |
2RC2 |
du |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
C |
|
C |
. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
dt2 |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|||
Сформируем уравнение Эйлера – Лагранжа:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dF |
|
|
|
d |
|
|
|
|
|
dF |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0. |
|
(а) |
|||||||
|
du |
|
|
|
|
|
|
du |
|
|
|
|
|||||||||||
|
|
|
|
dt |
|
|
|
|
С |
|
|
|
|
|
|||||||||
|
|
С |
|
d |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
||||||
Cu LC2 |
d 2u |
2RC2 |
du |
|
|
||||||||||||||||||
|
|
|
|
С |
|
С |
0, |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
C |
|
|
|
|
|
dt |
|
|
|
|
|
|
|
dt |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
преобразуем его к виду |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2u |
|
|
R du |
|
|
|
|
u |
|
|
|
|
|||||||||
|
|
|
С |
|
2 |
|
|
|
|
|
|
С |
|
|
|
С |
|
0 |
|||||
|
dt |
2 |
|
L |
dt |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
LC |
|
|
||||||||||||
и получим характеристическое уравнение
p2 2 RL p LC1 0.
Рассмотрим параллельную цепь с элементами ветвей R, L и C без источника питания. Конденсатор предварительно заряжен до напряжения U0 :
iL iR iC 0,
где iL , iR , iC – точки ветвей с элементами R, L и С.
39
Подробно опишем каждое слагаемое этого выражения:
uС L didtL ;
didtL L1 uС;
diL L1 uСdt;
iL L1 ∫uСdt;
uRС L1 ∫uСdt C dudtС 0.
Продифференцируем приведенное выше уравнение:
|
|
d 2u |
|
|
1 du |
|
|
|
1 |
|
|
|
|||||||||
|
C |
|
С |
|
|
|
|
|
|
С |
|
|
|
|
u |
0, |
|||||
dt2 |
R |
|
|
|
dt |
|
L |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
С |
|
||||||||
упростим его: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2u |
|
1 |
|
|
du |
|
|
|
1 |
|
|
|
||||||||
|
|
С |
|
|
|
|
|
|
|
|
С |
|
|
|
|
|
|
|
|
u |
0 |
|
dt2 |
RC |
|
|
|
dt |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
LC С |
|||||||||||||
и получим характеристическое уравнение |
|
|
|
||||||||||||||||||
|
|
p2 |
|
1 |
|
|
|
p |
|
1 |
|
|
0. |
|
|||||||
|
|
RC |
|
LC |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Составим функцию энергий для параллельной R–L–C-цепи с предварительно заряженной емкостью до напряжения U0 :
CU 2 |
L |
iL2 |
u |
2 |
CU 2 |
|||||
С |
|
|
R∫ |
С |
|
|
0 |
. |
||
2 |
2 |
R |
2 |
|||||||
|
|
|
|
|
|
|||||
40
