Выполнение расчетно-графических работ по ТОЭ в среде MathCAD
.pdfвыбирают условное направление тока в одну точку, например b, которую заземляют;
амперметры в ветвях схемы используют в качестве шунтов для съема сигналов тока;
один из выводов осциллографа, например А, присоединяют к точке а для
измерения узлового напряжения Uab, на второй вход осциллографа поочередно подают токовые сигналы с точек 1, 2, 3.
В результате получают углы сдвига фаз токов I1, I2, I3 с выбранным условным направлением векторов относительно узлового напряжения
Uab, вектор которого принимают с нулевой начальной фазой, как базовый.
Такой эксперимент служит неплохой проверкой правильности расчетов, а, главное, позволяет с помощью осциллографа увидеть синусоидально изменяющиеся величины во времени.
3. Расчетно-графическая работа №3 Расчет и анализ трехфазной цепи синусоидального тока при
симметричном источнике питания
Задание Для заданной электрической цепи выполнить следующее:
1.Рассчитать фазные и линейные токи и напряжения приёмника.
2.Определить мощность потребления энергии и коэффициент мощности приемника.
21
3. Построить топографическую диаграмму, совмещенную с векторной диаграммой токов.
3.1 Выполнение РГР№3 в среде MathCAD
Схема электрической цепи приведена на рис. 3.1, ее исходные данные: Uф=220,
Z4=15+25j –
сопротивление приемника в фазе а, Z5=15-25j–
сопротивление приемника в фазе b, Z6=25j– сопротивление приемника в фазе с, Zl=8+4j – сопротивление линии,
Zn=10+15j –
сопротивление нулевого провода, f=50Гц.
На рабочем листе эти значения следует записать следующим образом:
U := 220 |
Z5 |
:= |
15 − 25i |
Zl := 8 + 4i |
Z4 := 15 + 25i |
Z6 |
:= |
25i |
Zn := 10 + 15i |
Определим значения ЭДС в комплексном виде, учтя, что они смещены относительно друг друга на 120 градусов.
90iπ
Ea := U e |
180 |
Ea = 220i |
|
|
|
330iπ |
|
Eb := U e |
180 |
Eb = 190.526− 110i |
|
|
|
210iπ |
|
Ec := U e |
|
180 |
Ec = −190.526− 110i |
Так как нейтральный провод обладает сопротивлением, необходимо найти напряжение смещения нейтрали. Для этого достаточно записать только формулу и, если все значения сопротивлений и ЭДС были определены ранее, то MathCAD сам подставит цифры, и решение будет выглядеть следующим образом:
22
|
|
|
|
Ea |
|
+ |
|
Eb |
|
+ |
|
Ec |
|
|
|
|
|||
UnN:= |
|
|
Zl + Z4 |
Zl + Z5 |
Zl + Z6 |
|
|
UnN = −23.71+ 106.108i |
|||||||||||
|
|
|
|
|
|
||||||||||||||
|
1 |
|
+ |
|
1 |
|
+ |
|
1 |
|
+ |
|
1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Zl + |
|
|
Zl + |
|
|
Zl + |
|
|
Zn |
||||||||
|
|
Z4 |
Z5 |
Z6 |
Найдя напряжение смещения, определим токи цепи, для этого просто запишем закон Ома в общем виде, а цифры в формулу подставит программа. Тут же найдем модули токов, нажав на панели “калькулятор” кнопку со значком модуля.
Ia := |
Ea − UnN |
Ia = 2.809 + 1.41i |
|
Ia |
|
|
|
|
= 3.143 |
|
|
|
|||||||||
|
Z4 + Zl |
|
|
|
|
|
|
|
|
|
Ib := |
Eb − UnN |
Ib = 9.758 − 0.486i |
|
Ib |
|
|
= 9.771 |
|||
|
|
|||||||||
|
Z5 + Zl |
|
|
|
|
|
|
|
|
|
Ic := |
Ec − UnN |
Ic = −8.4 + 3.435i |
|
Ic |
|
|
= 9.075 |
|||
|
|
|||||||||
|
Z6 + Zl |
|
|
|
|
|
|
|
|
|
IN := |
UnN |
IN = 4.168 + 4.359i |
|
IN |
|
= 6.031 |
||||
|
|
|||||||||
|
Zn |
|
|
|
|
|
|
|
= 6.031 |
|
In := Ia + Ib + Ic |
In = 4.168 + 4.359i |
|
In |
|
|
|||||
|
|
Здесь как и ранее выражения комплексных величин в MathCADe записаны без черточек.
Расчет тока в нейтральном проводе для проверки выполнен двумя способами.
Зная линейные токи (при соединении звездой они же являются и фазными), можно определить фазные напряжения:
Ua := Ia Z4 |
Ua = 6.88 + 91.375i |
Ub := Ib Z5 |
Ub = 134.224− 251.253i |
Uc := Ic Z6 |
Uc = −85.878− 209.99i |
После того как токи и напряжения приёмников определены, можно приступить к определению выделившейся мощности и нахождению показаний ваттметров.
23
определение мощности потребления энергии
и коэффициента мощности
активная и реактивная мощности, выделившиеся в приемниках:
P1:=( |
|
|
Ia |
|
)2 Re(Z4) + ( |
|
|
Ib |
|
)2 Re(Z5) + ( |
|
|
Ic |
|
)2 Re(Z6) |
P1 = 1.58× 103 |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
Q1:=( |
|
|
Ia |
|
)2 Im(Z4) + ( |
|
Ib |
|
|
)2 Im(Z5) + ( |
|
|
Ic |
|
)2 Im(Z6) |
Q1= −80.784 |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
определим активную и реактивную мощности, |
|
||||||||||||||||||||||||||||||||||||||||||
выделившиеся на проводах |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
P2:=( |
|
|
Ia |
|
)2 Re(Zl) + ( |
|
|
Ib |
|
)2 Re(Zl) + ( |
|
|
Ic |
|
)2 Re(Zl) |
P2 = 1.502× 103 |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
Q2:=( |
|
|
|
Ia |
|
)2 Im(Zl) + ( |
|
Ib |
|
|
)2 Im(Zl) + ( |
|
|
Ic |
|
)2 Im(Zl) |
Q2= 750.782 |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
Pn :=( |
|
|
|
|
In |
|
)2 Re(Zn) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pn = 363.724 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
Qn:=( |
|
In |
|
|
)2 Im(Zn) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Qn = 545.586 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
определим полную мощность |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
S := |
(P1 + P2 + Pn)2 + (Q1+ Q2+ Qn)2 |
|
|
|
|
S = 3.654× 103 |
|||||||||||||||||||||||||||||||||||||
определим коэффициент мощности cos(α) |
|
||||||||||||||||||||||||||||||||||||||||||
cos α := |
|
P1 + P2 + Pn |
|
|
|
|
|
|
|
|
|
cos α = 0.943 |
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
найдем показания ваттметров |
(Ua Ia) |
|
|||||||||||||||||||||||||||||||||||||||||
W1 := Ua (2.812 − 1.417i) |
|
|
|
|
|
|
|
|
|
W1 := Re |
Re(W1) = 148.824 |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
W2 := Ub (9.738 + 0.495i) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Re(W2) = 1.431× 103 |
||||||||||||||||||||||||
W3 := Uc (−8.389 − 3.424i) |
|
|
|
|
|
|
|
|
Re(W3) = 1.424 |
||||||||||||||||||||||||||||||||||
сравним показания ваттметров с расчетами |
|
P1 = 1.58× 103
Re(W1) + Re(W3) + Re(W2) = 1.582× 103
Равенство мощностей, посчитанных разными способами, подтверждает правильность расчетов. Теперь можно приступить к расчету потенциалов различных точек схемы, т.е. построению топографической диаграммы.
24
Для этого заземляем любую точку схемы, например точку N, т.е. считаем, что ее потенциал равен нулю, а все остальные потенциалы схемы находим относительно заземленной точки. Например:
ψa := Ea − Ia Zl |
ψa = −16.831+ 197.483i |
ψb := Eb − Ib Zl |
ψb = 110.514− 145.145i |
ψc := Ec − Ic Zl |
ψc = −109.588− 103.883i |
ψn := In Zn |
ψn = −23.71 + 106.108i |
3.2Выполнение РГР 3 в среде EWB
После запуска программы EWB на рабочем поле собирают заданную схему. При указании параметров элементов учитывают следующие обстоятельства:
Начальные фазы источников напряжения могут быть только положительными, поэтому если фаза А принята с начальной фазой 90°, соответственно у фазы В будет 330°, а у фазы С 210°.
Для реактивных элементов по заданным сопротивлениям подсчитывают соответствующее значение емкости и индуктивности:
L = |
X |
L |
(Гн);C = |
10 |
6 |
(мкФ) . |
2πf |
|
|
||||
|
|
2πfX C |
С учетом принятых условий рабочая схема имеет следующий вид:
Приняв потенциал точки n за ноль, т.е. заземлив ее, находим с помощью осциллографа положение векторов фазных токов Ia, Ib, Ic и тока нейтрали In относительно напряжения смещения UNn . Здесь UNn, а в
расчетах приведённых ранее, было UnN. Помнить UNn= -UnN !
С помощью амперметров и вольтметра, указанных на схеме, определяют действующие значения фазных токов и напряжение
25
смещения нейтрали. Для этого опцию измеряемых величин в используемых приборах надо переустановить с DC на AC.
4. Расчетно-графическая работа №4 Расчет переходных процессов в линейных электрических
цепях
Расчет переходных процессов в линейных цепях постоянного тока
При расчете переходных процессов в линейных цепях в основном используют три метода:
•Классический с нахождением постоянных интегрирования.
•Численный метод решения, выполняемый на компьютере.
•Операторный метод.
4.1 Выполнение РГР№4 в среде MathCAD
Последовательность операций проиллюстрируем на конкретном примере. В электрической цепи постоянного тока замыкается ключ, требуется найти переходные функции uc(t), il(t), ic(t), ul(t)
Схема цепи приведена на рис 4.1, ее исходные данные:
Рис 4.1
Способ 1 Классический метод.
1.Составление характеристического уравнения. Составляем выражение входного сопротивления Z(p) записанное относительно источника после коммутации, приняв катушку за элемент p L C
сопротивления pL , а |
Z(p) := R1 + R3 + |
|
|
|
|
|
2 |
|
|
||
конденсатор 1/pC. |
|
p L + R2 + p |
R2 |
C L |
|
|
|
|
26
Z(p) := p2 C L |
R2 |
(R1 + R3) + p L (R1 + R3 + C) + R2 (R1 + R3) |
|||||||||||||||
Z(p) → |
|
9 |
|
p2 + |
|
900001 p + 900 |
|
|
|
|
|
||||||
40000 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
4000000 |
|
|
|
|
|
|
|
|
|||||
|
|
900001 |
|
|
900001 |
2 |
|
9 |
|
|
|
|
|||||
|
− 4000000 |
+ |
|
4000000 |
− 4 900 |
|
|
|
|
|
|||||||
40000 |
|
3 |
|||||||||||||||
p1 := |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p1 |
= −500.001+ 1.936i× 10 |
|
|
|
|
|
|
|
2 |
9 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
40000 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
900001 |
|
|
900001 |
2 |
|
9 |
|
|
|
|
|||||
|
− 4000000 |
− |
|
4000000 |
− 4 900 |
|
|
|
|
|
|
||||||
|
|
40000 |
|
|
3 |
||||||||||||
p2 := |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p2 |
= −500.001− 1.936i× 10 |
|
|
|
|
|
|
|
2 |
9 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
40000 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2.определение времени переходного процесса Т
3 |
|
5 |
− 3 |
|
δ := 500.001 ω := 1.936i 10 |
T := |
|
|
T = 10 × 10 |
δ |
3.по виду корней можно определить характер свободной составляющей
тока и напряжения, а именно:
•если корни p1 и p2 комплексно сопряженные, то свободная составляющая запишется: A e− δt sin(ω t + β)
•в случае отрицательных корней p1 и p2 свободная составляющая запишется: A1 ep1 t + A2 ep2 t
•если корни p1 и p2 отрицательные и равны (для схем с двумя реактивными элементами), то свободная составляющая будет представлена: (A1 + A2 t)ep t .
4.Определение начальных условий. Ток в катушке и напряжение на конденсаторе до коммутации, согласно основным законам коммутации скачком измениться не могут, т.е. il(0-)=il(0+) и uc(0-)=uc(0+). В нашем случае источник до коммутации был отключен, поэтому ток в катушке и напряжение на конденсаторе до коммутации были равны нулю.
Il := 0 |
Uc := 0 |
5.Расчет ul(0+) и ic(0+) произведем на основе законов Кирхгофа в схеме после коммутации
27
Ul := 0 |
|
I := 0 Ic := 0 |
Ir := 0 |
||||
Given |
|
|
|
||||
I |
|
|
|
Il + Ic + Ir |
|
|
|
|
|
|
|
||||
E |
|
|
I (R1 |
+ R3) |
+ Ir |
R2 |
|
|
|
||||||
E |
|
|
I (R1 |
+ R3) |
+ Uc |
|
|
|
|
|
|||||
E |
|
|
I (R1 |
+ R3) |
+ Ul |
|
|
|
|
|
5.556
Find(I,Ir,Ic,Ul) = 05.556
0
6.Расчет принужденных значений токов и напряжения выполняют,
исходя из того, что на постоянном токе емкость – обрыв, а катушка
обладает нулевым сопротивлением.
Ilnp := |
E |
Ilnp = 5.556 |
Ucnp := 0 |
R1 + R3 |
7.Определение постоянных интегрирования производят с учетом корней характеристического уравнения:
Ic := Icnp + A e− δt sin(ω t + γ)
U |
:= U |
+ B e− δt sin(ω t + ν) |
; |
c |
cnp |
|
в программе MathCAD решение этой системы будет выглядеть следующим образом:
A := 0 |
γ := 0 |
|
|
|
|||
Given |
|
|
|
|
|
|
|
5.556 + A sin(γ) |
|
0 |
|
|
|
||
|
|
|
|
||||
|
|
|
|
||||
L A (−δ sin(γ) + ω cos (γ)) |
|
0 |
|||||
|
|||||||
|
|||||||
Find(A ,γ) = |
−5.752 |
|
|
|
|||
1.309 |
|
. |
|||||
|
28
Аналогичным образом найдем постоянные интегрирования для напряжения uc:
B := 0 |
|
|
ν := 0 |
|
|
|
|
Given |
|
|
|
|
|
|
|
B sin(ν) |
|
0 |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|||
C B (−δ sin(ν) + ω cos (ν)) |
|
5.556 |
|||||
|
|||||||
|
|||||||
Find(B,ν) = |
59.582 |
|
|
||||
0 |
|
. |
|||||
|
|
|
8.Представление искомых переходных функций.
Отыскав постоянные интегрирования, и определив вид свободной составляющей, запишем результат и построим графики, зависимостей il(t) и uc(t).
Il(t) := 5.556 − 5.752e− δt sin(ω t + 1.309) |
|
|
Uc(t) := 59.582 e− δt sin(ω t) |
|||||||||||||||
10 |
|
|
|
|
|
|
|
|
|
50 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
5 |
|
|
|
|
|
|
|
Uc(t) |
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
Il(t) |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
50 |
|
|
|
|
|||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0.005 |
|
|||
|
0 |
0.005 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
t |
|
|
|
|
|
|
t |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Недостатки классического метода:
•большой объем и громоздкость вычислений, в результате чего есть большая вероятность ошибиться, где-либо.
Способ 2 Численный метод решения системы дифференциальных уравнений
Система MathCAD позволяет непосредственно решать системы дифференциальных уравнений численным методом, с помощью встроенной функции rkfixed.
29
rkfixed(y, x1, x2, npoints, D) – общий вид функции, возвращающей матрицу, в который первый столбец содержит значение аргумента функции, а оставшийся столбец содержит соответствующие значения искомой функции.
•Y должен быть матрицей nx1 первоначальных значений искомых функций, в нашем случае это значения тока на катушке и напряжения на емкости до коммутации.
•X1, X2 - границы интервала, на котором будет решаться система дифференциальных уравнений, в нашем случае это время переходного процесса (от t=0 до t=T=0.01).
•Npoints - число точек на интервале интегрирования, в которых искомые функции будут аппроксимированы.
•D - это n-элементная вектор-функция (матрица), содержащая первые производные неизвестных функций.
Наиболее важным шагом в решении дифференциального уравнения в программе MathCAD является корректное представление вектора D, в правой части которого кроме напряжения источника могут быть записаны только ток индуктивности и напряжение на емкости.
Составим и решим систему уравнений в дифференциальной форме:
i |
|
|
|
ir + il + C d uc ; |
|
|
|
|
|
dt |
|
E |
|
|
i (R1 + R3) |
+ ir R2 |
|
|
|
||||
E |
|
|
i (R1 + R3) |
+ L d il |
|
|
|
|
|
|
dt |
E |
|
|
i (R1 + R3) |
+ uc . |
|
|
|
;
;
Выразим из полученной системы уравнений производные тока на катушке и напряжения на конденсаторе:
d |
|
|
|
|
|
|
E − uc |
|
|
|
|
|
uc |
|
ic |
; |
||||
dtuc |
|
|
|
|
(R1 + R3) C − |
R2 C − |
C |
|
||||||||||||
d il |
|
|
uc |
. |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Примем : |
|
|
|
|
|
|
|
|
|
|
||||||||||
a |
|
|
|
|
R1 + R3 ; |
y |
1 |
|
|
uc ; |
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
il . |
|
|
||||
b |
|
|
R2 |
|
y |
2 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда в канонической форме матрицы Y и D запишутся следующим образом:
30