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

Прикладные задачи теоретической электротехники. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Исходное уравнение Дуффинга:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Ri2dt

CU02

.

2

2

2

 

 

 

 

Эта цепь последовательная, а значит

iL iC C dudtC ,

откуда получим функцию F:

CU 2

L

du

2

 

du

2

CU 2

 

C

 

 

C

C

 

RC

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]