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

Выполнение расчетно-графических работ по ТОЭ в среде MathCAD

.pdf
Скачиваний:
294
Добавлен:
26.03.2016
Размер:
1.09 Mб
Скачать

выбирают условное направление тока в одну точку, например 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.526110i

 

 

210iπ

 

Ec := U e

 

180

Ec = −190.526110i

Так как нейтральный провод обладает сопротивлением, необходимо найти напряжение смещения нейтрали. Для этого достаточно записать только формулу и, если все значения сопротивлений и ЭДС были определены ранее, то 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.224251.253i

Uc := Ic Z6

Uc = −85.878209.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.514145.145i

ψc := Ec Ic Zl

ψc = −109.588103.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.0011.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