2
.pdf2) Метод последовательных приближений. Рассмотренная выше задача может быть решена аналитически методом последовательных приближений. В решении задачи может быть использована как графическая, так и аналитическая форма задания ВАХ нелинейной катушки, например, I = aU + bU5. Решение задачи выполняется в следующей последовательности:
а) жении б)
Задаются напряжением на нелинейном элементе в первом прибли-
U ;
L
Ток в цепи определяется из уравнения аппроксимации
I aU |
bU |
5 |
|
||
L |
L |
|
или графически по диаграмме ВАХ U(I);
в) Далее следуют вычисления:
г) Выполняется сравнение |
E |
|
U |
|
R |
|
E |
RI |
|
E |
U |
2 |
U |
2 |
; |
|
|
|
|
|
|||
|
|
|
R |
|
L |
|
, по результатам которого за-
даются напряжением |
U |
во втором приближении и повторяют расчет до |
|||||
L |
|||||||
|
|
|
|
E E |
0,01 |
|
|
достижения требуемой точности, например, |
E |
. |
|||||
|
|
||||||
|
|
|
|
|
|
Метод последовательных приближений применим к расчету схем любой сложности. Для сложной схемы составляется система уравнений Кирхгофа в комплексной форме, которая решается методом последовательных приближений, при этом все вычисления в отдельном цикле выполняются в комплексной форме. В качестве примера приведем расчет схемы рис. 73.
|
I1 U1(I1) |
|
R |
|
|||||
|
|
|
|
|
I2 |
|
|
|
I3 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
U3 |
|||
|
|
|
|
|
|
||||
|
|
|
U1 |
|
|
||||
|
|
|
|
|
|
|
|||
E |
|
|
|
|
U2 |
XC |
|
U4 |
I3(U4) |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
Рис. 73
Заданы параметры линейных элементов Е, R, XC. ВАХ нелинейных элементов заданы аналитически в виде уравнений аппроксимации:
U |
|
aI bI |
3 |
; |
I |
|
cU dU |
5 |
. |
R |
|
L |
|
||||||
|
|
|
|
|
|
|
Для исследуемой схемы система комплексных уравнений Кирхгофа совместно с уравнениями аппроксимации имеет вид:
101
I |
1 |
I |
2 |
I |
3 |
0 |
|
|
|
|
U |
1 |
U |
2 |
|
E |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
U |
2 |
U |
3 |
U |
4 |
0 |
|||||||
|
|
|
|
|
|
|
|
|
|
||||
U |
|
aI |
|
bI |
3 |
|
|
|
|||||
2 |
|
1 |
|
|
|
||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|||
I |
|
|
cU |
|
dU |
|
5 |
||||||
3 |
4 |
4 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
(1)
(2)
(3)
(4)
(5)
Один из вариантов решения полученной системы уравнений методом последовательных приближений представлен ниже.
а) Задаются в первом приближении комплексным напряжением на не-
линейной катушке, например: |
U |
|
50e |
j 0 |
. |
|
4 |
|
|||||
|
|
|
|
|||
б) Определяется модуль тока |
I |
аналитически из уравнения (5) или |
||||
3 |
графически по диаграмме функции IL(UL). Аргумент этого комплекса принимается равным 90о (в катушке ток отстает от напряжения на угол =
|
|
|
|
|
|
|
|
|
|
I |
I e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
j90 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 90 ). В комплексной форме |
3 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
в) Определяется напряжение на линейном резисторе по закону Ома: |
||||||||||||||||||||||||
U |
I |
R |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
г) |
Из |
|
уравнения |
(3) |
находится |
напряжение |
на конденсаторе: |
|||||||||||||||||
U 2 U 3 U 4 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
д) По закону Ома определяется ток конденсатора: |
I |
|
|
|
U |
2 |
|
. |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
2 |
|
jX |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
I |
|
I |
|
|
|
|
|
|
||||
|
|
е) Из уравнения (1) находится ток источника |
|
I e |
j |
. |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
1 |
|
|
2 |
|
|
3 |
|
|
|
1 |
|
|
||||||||||||
|
|
ж) Определяется модуль напряжения |
U |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
1 аналитически из уравнения |
||||||||||||||||||||||||
(4) |
или графически по диаграмме функции UR(IR). Аргумент этого ком- |
|||||||||||||||||||||||||
плекса принимается равным аргументу комплекса тока |
|
I |
|
(в резисторе |
||||||||||||||||||||||
|
1 |
|||||||||||||||||||||||||
ток совпадает с напряжением). В комплексной форме |
U |
U e |
j |
. |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
1 |
|
|
1 |
|
|
|
|
|
|
||||||||||||||||
|
|
з) |
Из |
|
уравнения |
(2) |
|
находится расчетное |
|
значение |
|
ЭДС: |
||||||||||||||
E U |
U |
|
E e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
j |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и) Сравнивают найденное в первом приближении значение модуля |
||||||||||||||||||||||||
ЭДС |
E |
с заданным значением ЭДС Е и с учетом вида полученного нера- |
||||||||||||||||||||||||
|
||||||||||||||||||||||||||
венства |
E E |
задаются новым значением напряжения |
U |
4 |
|
во втором |
||||||||||||||||||||
|
|
|
|
|
|
|
|
приближении и повторяют расчет по тому же алгоритму. Циклы расчета (итерации) повторяют до достижения желаемой точности. В результатах последнего цикла корректируют аргументы комплексных токов и напряжений путем добавления к ним значения .
3) Расчет нелинейных цепей на ЭВМ. Система нелинейных уравнений, составленная для исследуемой схемы по законам Кирхгофа, может быть решена на ЭВМ методом последовательных приближений как по ин-
102
дивидуальной программе, так и по стандартной программе из пакета MathCAD или MATLAB. При этом система уравнений Кирхгофа дополняется нелинейными уравнениями аппроксимации в комплексной форме. В комплексной форме аппроксимации учитывается как нелинейная связь между модулями величин U и I, так и угол сдвига между ними. Ниже приведены примеры аппроксимации ВАХ U(I) или I(U) для нелинейных элементов разного рода:
Ur |
b |
Ir |
|
Il |
|
d |
Ul |
|
|||
|
|||
Ic |
d |
Uc |
3 ej (arg(Ir))
5 ej (arg(Ul) 90 deg)
5 ej (arg(Uc) 90 deg)
–для резистора,
–для катушки,
–для конденсатора.
Последовательность расчета нелинейной цепи на ЭВМ по стандартной программе из пакета MathCAD приведен ниже.
а) Задана схема нелинейной цепи (рис. 73) и параметры линейных элементов Е, R, XC. ВАХ нелинейных элементов заданы аналитически в виде уравнений аппроксимации в комплексной форме:
|
E 70 |
|
|
R 25 |
U1 |
b I1 |
|
2 |
j arg(I1) |
|
e |
Xc 30 |
|
I3 |
d |
|
b 26 |
|||
|
U4 |
|
4 |
j |
|
e |
d |
1.8 10 |
3 |
|
||
(arg(U4) 90deg) |
б) Составляется система уравнений по законам Кирхгофа, которая дополняется уравнениями аппроксимации в комплексной форме для нелинейных элементов. Полученная таким образом система уравнений решается по стандартной программе:
Given
I1 I2 |
I3 |
0 |
U1 j I2 Xc |
E |
U1 I3 R U4 |
E |
|
I3 |
|
|
I1 |
|
|
I2 |
|
|
||
|
||
|
I3 |
|
|
U1 |
|
|
||
|
||
|
U4 |
d |
|
U4 |
|
|
4 |
j (arg(U4) 90deg) |
|
|
|||||
|
|
|
e |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Find (I1 I2 I3 U1 U4) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
I2 ( j Xc) |
|
|
|
U2 |
||||
|
|
|
|
|
|
U1 |
b |
I1 |
|
2 |
j arg(I1) |
|
e |
U3 I3 R
в) Результаты расчета приводятся к стандартной форме:
I1 |
|
I2 |
|
I3 |
|
U1 |
|
U2 |
|
U3 |
|
U4 |
|
1.315 0.936 1.105
44.984
28.073
27.628
4.978
arg(I1) |
13.036 |
deg |
arg(I2) |
68.812 |
deg |
arg(I3) 31.402 deg arg(U1) 13.036 deg arg(U2) 21.188 deg arg(U3) 31.402 deg
arg(U4) 58.598 deg
103