УП Нелинейные цепи 2011
.pdfа |
в |
б |
г |
д е
Рис. 3.9, лист 1. Схемы магнитных цепей к самостоятельной работе № 2
40
ж |
з |
|
|
и |
|
|
|
|
|
к |
|
|
|
Рис. 3.9, лист 2. Схемы магнитных цепей к самостоятельной работе № 2 |
|||||||||||
|
|
|
|
|
|
|
|
|
Таблиц а 3.6 |
||
Зависимость для основной кривой намагничивания применяемой стали |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
H, А/м |
20 |
40 |
60 |
80 |
120 |
200 |
400 |
600 |
|
800 |
1200 |
|
|
|
|
|
|
|
|
|
|
|
|
B, Тл |
0,22 |
0,75 |
0,93 |
1,02 |
1,14 |
1,28 |
1,47 |
1,53 |
|
1,57 |
1,6 |
|
|
|
|
|
|
|
|
|
|
|
|
3.3.2Пример расчета
Вразветвленной магнитной цепи (рис. 3.10) магнитопровод выполнен из электротехнической стали марки 1211, кривая намагничивания которой приведена в табл. 3.7 с параметрами обмоток и геометрическими размерами: I1 = 5 A,
w1 = 200, I2 = 3 A, w2 = 100, a = b = 10 мм, a1 = 15 мм, a2 = 25 мм, c = 40 км, l02 = l03 = 1 мм.
41
Требуется выполнить следующее:
1.Составить систему уравнений по законам Кирхгофа для магнитной цепи.
2.Определить магнитные потоки и индукции в стержнях.
Рис. 3.10. Расчетная схема разветвленной магнитной цепи
Таблиц а 3.7
Зависимость для основной кривой намагничивания стали 1211
H, А/м |
140 |
211 |
318 |
502 |
843 |
1580 |
4370 |
7780 |
B, Тл |
0,4 |
0,6 |
0,8 |
1,0 |
1,2 |
1,4 |
1,6 |
1,7 |
Решение
1. Разбиваем магнитную цепь на участки одинакового материала и равно-
го сечения. Длины участков определяем по средней линии:
|
|
|
l1 c 2a1 3a; |
l1 0,1м; |
|
||||
|
|
|
l2 c a l02; |
l2 0,049 м; |
|
||||
|
l3 c 2a2 3a l03; |
|
l3 0,119 м; |
|
|||||
S S |
2 |
S a b; |
S S |
2 |
S |
3 |
10 4 |
м2; |
|
1 |
|
3 |
1 |
|
|
|
|||
|
|
|
S02 S2; |
S03 S3. |
|
|
|
42
|
Определяем по правилу правой руки положительное направление МДС: |
|||||
|
|
m |
охватываем правой рукой витки с |
|||
|
|
током в направлении тока; отогну- |
||||
I1w1 |
R02 |
R03 |
тый большой палец руки указывает |
|||
|
направление |
МДС. |
Направления |
|||
Ф1 |
I2w2 |
Ф3 |
магнитных |
потоков |
выбираются |
|
|
произвольно. |
|
|
|
||
|
Ф2 |
UMmn |
|
|
|
|
RM |
|
Эквивалентная |
схема |
заме- |
||
|
RM3 |
|||||
1 |
RM2 |
щения магнитной цепи представле- |
||||
|
||||||
|
|
n |
на на рис. 3.11. |
|
|
|
|
|
Составляем уравнения по за- |
||||
Рис. 3.11. Эквивалентная схема |
конам Кирхгофа, выбирая направ- |
|||||
|
замещения магнитной цепи |
ления и пути обхода контуров (по |
||||
второму закону Кирхгофа), как показано на рисунке: |
|
|
|
Ф1 Ф2 Ф3 0; |
|
||
|
|
I1w1; |
|
H1l1 UMmn |
(3.24) |
||
|
H2l2 |
|
|
UMmn |
H02l02 I2w2; |
||
U |
H l |
H l |
0. |
Mmn |
3 3 |
03 03 |
|
Из полученных уравнений выразим напряжение UMmn
U |
I w |
H l ; |
|
Mmn |
1 1 |
1 1 |
(3.25) |
UMmn H2l2 |
H02l02 I2w2; |
||
|
H3l3 |
H03l03. |
|
UMmn |
|
2. Необходимо построить характеристики для каждой ветви магнитной
цепи:
Ф1 = f1(UMmn) f1(I1w1 H1l1); |
(3.26) |
|
Ф2 = f2(UMmn) f2 |
(H2l2 H02l02 I2w2); |
(3.27) |
Ф3 = f3(UMmn) |
f3(H3l3 H03l03). |
(3.28) |
С этой целью задаемся значениями магнитной индукции на участках в пределах B 0 1,7 Тл,определяем магнитные потоки Ф= B S, напряженности магнитного поля H – по кривой намагничивания (см. табл. 2) и магнитные на-
пряжения UMmn. Результаты вычислений сводим в табл. 3.8 – 3.10.
43
Таблиц а 3.8
Результаты вычислений для первого участка
|
B1, Тл |
|
Ф , 10 4 Вб |
|
H1, А/м |
|
H1l1, А |
|
|
UMmn I1w1 H1l1, A |
|||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
1000 |
|
|
|||
0,4 |
|
|
0,4 |
|
|
|
|
140 |
|
|
14,0 |
|
|
|
|
|
|
986,0 |
|
|
|||||
0,6 |
|
|
0,6 |
|
|
|
|
211 |
|
|
21,1 |
|
|
|
|
|
|
978,9 |
|
|
|||||
0,8 |
|
|
0,8 |
|
|
|
|
318 |
|
|
31,8 |
|
|
|
|
|
|
968,2 |
|
|
|||||
1,0 |
|
|
1,0 |
|
|
|
|
502 |
|
|
50,2 |
|
|
|
|
|
|
949,8 |
|
|
|||||
1,2 |
|
|
1,2 |
|
|
|
|
843 |
|
|
84,3 |
|
|
|
|
|
|
915,7 |
|
|
|||||
1,4 |
|
|
1,4 |
|
|
|
1580 |
|
|
158,0 |
|
|
|
|
|
842,0 |
|
|
|||||||
1,6 |
|
|
1,6 |
|
|
|
4370 |
|
|
437,0 |
|
|
|
|
|
563,0 |
|
|
|||||||
1,7 |
|
|
1,7 |
|
|
|
7780 |
|
|
778,0 |
|
|
|
|
|
222,0 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблиц а 3.9 |
||
|
|
|
|
|
Результаты вычислений для второго участка |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
B , |
|
|
Ф2, |
|
H |
2 |
, |
|
H l , |
|
|
|
H02, |
|
H |
|
l , |
|
UMmn H2l2 H02l02 |
|||||
|
2 |
|
|
|
|
|
|
2 2 |
|
|
|
|
|
|
|
02 02 |
|
|
|
I2w2, A |
|
|
|||
|
Тл |
10 4Вб |
А/м |
|
|
А |
|
10 4 А/м |
|
|
А |
|
|
|
|
|
|||||||||
|
0 |
|
0 |
|
0 |
|
|
0 |
|
|
|
0 |
|
|
0 |
|
|
|
-300 |
|
|
||||
|
0,4 |
|
0,4 |
|
140 |
|
6,9 |
|
|
|
32 |
|
320 |
|
|
|
26,9 |
|
|
||||||
|
0,6 |
|
0,6 |
|
211 |
|
10,3 |
|
|
|
48 |
|
480 |
|
|
|
190,3 |
|
|
||||||
|
0,8 |
|
0,8 |
|
318 |
|
15,6 |
|
|
|
64 |
|
640 |
|
|
|
355,6 |
|
|
||||||
|
1,0 |
|
1,0 |
|
502 |
|
24,6 |
|
|
|
80 |
|
800 |
|
|
|
524,6 |
|
|
||||||
|
1,2 |
|
1,2 |
|
843 |
|
41,3 |
|
|
|
96 |
|
960 |
|
|
|
701,3 |
|
|
||||||
|
1,4 |
|
1,4 |
|
1580 |
77,4 |
|
|
|
112 |
|
1120 |
|
|
|
897,4 |
|
|
|||||||
|
1,6 |
|
1,6 |
|
4370 |
214,1 |
|
|
128 |
|
1280 |
|
|
|
1194,1 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблиц а 3.10 |
||
|
|
|
|
|
Результаты вычислений для третьего участка |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
B , |
|
|
Ф3, |
|
H |
3 |
, |
|
H l , |
|
|
|
H03, |
|
H03l03, |
|
U |
Mmn |
H l H l , |
|
||||
|
3 |
|
|
|
|
|
|
3 3 |
|
|
|
|
|
|
|
|
|
|
|
3 3 |
03 03 |
|
|||
|
Тл |
10 4Вб |
А/м |
|
|
А |
|
10 4 А/м |
|
|
А |
|
|
|
A |
|
|
||||||||
|
0,4 |
|
0,4 |
|
140 |
|
16,7 |
|
|
|
32 |
|
320 |
|
|
|
336,7 |
|
|
||||||
|
0,6 |
|
0,6 |
|
211 |
|
25,1 |
|
|
|
48 |
|
480 |
|
|
|
505,1 |
|
|
||||||
|
0,8 |
|
0,8 |
|
318 |
|
37,8 |
|
|
|
64 |
|
640 |
|
|
|
677,8 |
|
|
||||||
|
1,0 |
|
1,0 |
|
502 |
|
59,7 |
|
|
|
80 |
|
800 |
|
|
|
859,7 |
|
|
||||||
|
1,2 |
|
1,2 |
|
843 |
|
100,3 |
|
|
96 |
|
960 |
|
|
|
1060,3 |
|
|
|||||||
|
1,4 |
|
1,4 |
|
1580 |
188,0 |
|
|
112 |
|
1120 |
|
|
|
1308,0 |
|
|
||||||||
|
1,6 |
|
1,6 |
|
4370 |
520,0 |
|
|
128 |
|
1280 |
|
|
|
1800,0 |
|
|
||||||||
|
1,7 |
|
1,7 |
|
7780 |
925,8 |
|
135,4 |
|
1350 |
|
|
|
2275,8 |
|
|
44
На основании результатов вычислений строим зависимости Ф1(UMmn), |
|||||
Ф2(UMmn), Ф3(UMmn) (рис. 3.12). |
|
|
|
|
|
Вб |
|
|
|
|
|
10 4 |
|
|
c |
|
|
1,6 |
|
|
|
|
|
|
|
|
|
|
|
1,4 |
|
|
|
|
Ф2 |
|
|
|
|
|
|
1,2 |
Ф2 Ф3 |
|
|
|
|
|
|
|
|
Ф3 |
|
|
|
|
|
|
|
1,0 |
|
|
d |
|
|
|
|
|
|
|
|
0,8 |
|
|
|
|
Ф1 |
|
|
|
e |
|
|
Ф 0,6 |
|
|
|
|
|
|
|
|
|
|
|
0,4 |
|
|
|
|
|
0,2 |
|
|
|
|
|
|
|
|
f |
|
|
0 |
2 |
4 |
6 |
8 |
10 100 A |
|
UMmn |
|
|
|
|
Рис. 3.12. Процесс и результаты расчета магнитных потоков |
Решение для магнитных потоков, удовлетворяющее уравнению
Ф1 Ф2 Ф3, находим графически. Для этого строим вспомогательную кривую
Ф2 Ф3 f (UMmn) (см. рис. 3.12), суммируя ординаты кривых Ф2(UMmn) и
Ф3(UMmn), и находим точку пересечения ее с кривой Ф1(UMmn).
Абсцисса точки с определяет значение UMmn,ординаты точек c, d, e – со-
ответственно значения магнитных потоков Ф1, Ф2, Ф3.
Магнитные индукции на участках цепи, Тл: B1 Ф1 1,61; В2 = В02 = 0,99;
S1
В3 = В03 = 0,62.
45
4. НЕЛИНЕЙНЫЕ ЭЛЕКТРИЧЕСКИЕ ЦЕПИ ПЕРЕМЕННОГО ТОКА
Нелинейные цепи описываются нелинейными дифференциальными урав-
нениями, к которым неприменим принцип наложения, и решение не может быть представлено как сумма принужденной и свободной составляющих. Если все элементы цепи сосредоточенные, то уравнения относятся к классу обыкно-
венных дифференциальных уравнений. Тип дифференциального уравнения за-
висит от способа аппроксимации нелинейных характеристик и вида НЭ в схеме.
Существуют различные способы расчета нелинейных электрических цепей.
4.1. Аппроксимация характеристик нелинейных элементов
Характеристики НЭ, как уже указывалось, определяются эксперимен-
тально. Их непосредственное использование возможно при графических спосо-
бах расчета нелинейных цепей. Если исследование производится аналитиче-
ским или численным методом, то требуется аппроксимация (представление)
нелинейных характеристик математическими выражениями. Сложные характе-
ристики трудно или невозможно аппроксимировать единым выражением, по-
этому используется аналитическая аппроксимация отдельных участков или так называемая кусочно-линейная аппроксимация.
Основными критериями при выборе аппроксимирующего выражения яв-
ляются следующие:
а) точность представления заданной нелинейной характеристики;
б) степень сложности, вносимая в решение задачи выбранным типом ап-
проксимации.
Аналитическое выражение, аппроксимирующее характеристику НЭ, оп-
ределяет тип, а следовательно, и особенности решения нелинейных уравнений.
4.1.1. Аппроксимация степенным полиномом
Аппроксимация характеристик НЭ может осуществляться с помощью степенного ряда, показательной, тригонометрических, гиперболических и дру-
гих функций. Наиболее часто применяется аппроксимация степенными полино-
мами вида
i(u) a |
a (u) a u2 a u3 ... |
(4.1) |
|||
0 |
1 |
2 |
3 |
|
|
или |
|
|
|
|
|
u(i) b |
b (i) b i2 |
bi3 |
... |
(4.2) |
|
0 |
1 |
2 |
3 |
|
|
|
46 |
|
|
|
|
Эти полиномы при большом количестве членов достаточно точно ап-
проксимируют нелинейную характеристику. Но при этом возрастает сложность расчетов. Поэтому, как правило, в (4.1) и (4.2) берут ограниченное количество членов.
Для нахождения коэффициентов ak или bk может быть использовано не-
сколько способов. Простейший из них состоит в следующем.
На исходной кривой выбирается столько точек, сколько коэффициентов входит в аппроксимирующее выражение (4.1) или (4.2). Пусть
i(u) au a |
u2 a u3. |
(4.3) |
|
1 |
2 |
3 |
|
В этом случае на характеристике НЭ выбираются три точки с координа-
тами i1,u1;i2,u2;i3,u3. Для каждой из них записывается (4.3) и формируется система уравнений:
i au a |
u2 |
a u3; |
|
|||||||||
|
1 |
1 |
1 |
2 |
|
|
1 |
3 |
1 |
; |
(4.4) |
|
i |
au |
2 |
a |
|
|
u2 |
a u3 |
|||||
|
2 |
1 |
|
2 |
2 |
3 |
2 |
|
|
|||
i au |
3 |
a |
|
u2 |
a u3. |
|
||||||
|
3 |
1 |
|
2 |
|
3 |
3 |
3 |
|
|
Решение (4.4) дает коэффициенты a1, a2, a3 .
При данном способе совпадение исходной характеристики и аппроксими-
рующей функции гарантируется только в выбранных точках. Чем ближе распо-
ложены точки, тем точнее представляется участок между ними. Поэтому наи-
лучшие результаты могут быть получены, если, например, (4.4) описывает не всю кривую, а только рабочий ее участок, когда таковой можно знать априорно,
то есть до решения задачи.
Кроме рассмотренного способа, применяется еще приближение в среднем и в квадратичном смысле. В последнем случае применяется метод наименьших
квадратов.
Примером использования полиномиальной аппроксимации характеристи-
ки НЭ может служить анализ кривой тока цепи при синусоидальном приложен-
ном напряжении.
Если к НЭ, характеристика которого аппроксимирована полиномом
третьей степени, подвести синусоидальное напряжение u Um sin t, то выра-
жение для тока будет иметь вид:
i a |
aU |
m |
sin t a U2 |
sin2 t aU3 |
sin3 t. |
(4.5) |
0 |
1 |
2 m |
3 m |
|
|
47
Используя известные тождества: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
sin2 t |
1 |
|
1 |
cos2 t; |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
3 |
2 |
2 |
|
1 |
|
|
|
|
|
|
|
(4.6) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
sin3 t |
sin t |
sin3 t, |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
4 |
|
|
4 |
|
|
|
|
|
|
|
|
||||
получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a U2 |
|
|
|
3 |
|
3 |
|
1 |
|
2 |
|
1 |
|
3 |
|
|||||||
i a |
|
|
2 m |
|
aU |
m |
|
|
aU |
m |
sin t |
|
a U |
m |
cos2 t |
|
aU |
m |
sin3 t. (4.7) |
||||||
|
|
|
|
|
|||||||||||||||||||||
|
0 |
|
2 |
1 |
4 |
3 |
|
2 |
2 |
|
4 3 |
|
Следовательно, ток содержит постоянную составляющую и гармоники час-
тот, кратных частоте приложенного синусоидального напряжения, причем наи-
высший коэффициент кратности равен степени аппроксимирующего полинома.
Постоянная составляющая тока I |
|
a |
|
|
a U2 |
|
|
|
|
|
|||||||||||
0 |
0 |
|
2 |
m |
, амплитуда первой гармо- |
||||||||||||||||
|
2 |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ники I |
aU |
|
3 |
aU3 |
, амплитуда второй гармоники I |
|
|
1 |
a U2 |
, амплитуда |
|||||||||||
|
|
|
|||||||||||||||||||
1m |
1 m |
4 3 |
m |
|
|
|
|
|
|
|
|
|
|
2m |
|
2 2 m |
|
||||
третьей гармоники I |
|
|
1 |
aU3 . |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
3m |
|
|
4 3 m |
|
|
|
|
|
|
|
|
|
|
|
|
Приведенные выражения показывают, что постоянная составляющая и амплитуды четных гармоник тока зависят только от коэффициентов при четных степенях.
4.1.2. Метод кусочно-линейной аппроксимации
Суть этого метода состоит в том, что характеристика НЭ заменяется ло-
маной линией, состоящей из отрезков прямых. В пределах линейного участка
характеристики нелинейные дифференциальные |
u |
|
уравнения заменяются линейными с постоян- |
b |
|
ными коэффициентами. Решение нелинейного |
|
|
|
a |
|
уравнения (системы) заменяется последователь- |
|
|
|
|
|
ным решением нескольких линейных уравнений |
|
|
(систем) с припасовыванием в точках стыкова- |
|
|
ния линейных участков характеристики НЭ. |
0 |
i |
Если заменить реальные характеристики |
Рис. 4.1. Кусочно-линейная |
|
НЭ кусочно-линейными характеристиками (на- |
||
|
|
аппроксимация |
пример, рис. 4.1), то для расчета процессов в
характеристики НЭ
48
цепи можно воспользоваться нижеследующим методом.
В отдельные интервалы времени, пока во всех элементах цепи процессы находятся на определенных прямолинейных отрезках их характеристик, про-
цесс во всей цепи описывается совокупностью линейных дифференциальных уравнений, коэффициенты которых определяются параметрами этих линейных отрезков характеристик.
При переходе процесса в любом НЭ через точку излома характеристики
(точки а и b на рис. 4.1), изменяются параметры уравнений. Назовем момент ка-
ждого такого перехода моментом коммутации. Процесс за весь рассматривае-
мый промежуток времени, разбивается на интервалы, заключенные между дву-
мя любыми соседними моментами коммутации. Решения совокупности уравнений внутри каждого интервала содержат некоторое количество своих произвольных постоянных. Эти произвольные постоянные определяются из физических усло-
вий неизменности токов в индуктивных катушках и напряжений на конденса-
торах в моменты коммутации, т. е. путем сопряжения решений, полученных для двух смежных интервалов. Соответственно этот метод можно назвать методом с о п р я ж е н ия интервалов. Подлежат определению также моменты ком-
мутации из условий, что ток или напряжение достигает значения, соответст-
вующего точке излома характеристики.
Периодические процессы повторяются через период Т, и поэтому доста-
точно произвести расчет процессов в течение одного периода, используя усло-
вия, что значения токов в катушках и напряжений на конденсаторах одинаковы в начале и в конце периода. В случае симметричных многофазных цепей про-
цесс может повторяться за промежутки, составляющие целую долю периода.
Такой промежуток можно назвать и н т е р ва л ом повторяемо сти процесса. Очевидно, при этом достаточно произвести расчет в пределах интер-
вала повторяемости.
Метод сопряжения интервалов с успехом может быть применен, когда ха-
рактеристики нелинейных элементов состоят из отрезков близких к прямоли-
нейным, например, в случае элементов с ферритами, обладающими прямо-
угольной кривой намагничивания. Он широко используется для расчета цепей с полупроводниковыми вентилями.
Для примера рассмотрим расчет установившегося режима в вентильной
цепи.
49