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

Курсовой ЭТФ

.pdf
Скачиваний:
10
Добавлен:
03.03.2016
Размер:
1.55 Mб
Скачать

В программе предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц (согласовать с руководителем курсовой работы) проверить результаты работы программы.

а)

б)

A

 

IA

 

 

 

*

 

 

a

 

EAB

V

 

Z1

*

W1

Iab

 

Ica

UBC

 

V

 

 

xL

 

R

B

1

 

3

Ik1

 

 

 

 

IB

 

b *

xM

 

xC

UAB

 

V2

Z1

Ik2

 

* xL

 

R

ECB

 

 

Ibc

 

UCA

IC

 

 

 

 

 

c

C

 

 

*

 

 

 

 

 

 

W2

 

 

 

 

 

 

Z1

 

 

 

 

 

 

 

 

 

*

Ik3

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 2.8. Электрическая схема к задаче 8

 

Математическая модель

С использованием теоремы косинусов и с помощью построенной векторной диаграммы линейных напряжений (рис. 2.8, б) определим их комплексы, направив вдоль оси вещественных чисел вектор напряжения UAB:

U12

U22 U32

:= acos

 

 

 

 

 

 

2 U1 U2

 

 

 

UAB := U1

UBC := U2 ej ( - ) UCA := – UAB UBC

Припишем систему линейных напряжений двум ЭДС (рис. 2.8,а) EAB = UAB і ECB = –UBC, а расчет токов в этой схеме выполним методом контурных токов. Определим действительные и мнимые составляющие контуров.

Z11 = 2 Z1 + R + j xL

Z22 = 2 Z1 + R + j xL

Z33 = 2 Z1 j xC

Z12 = – Z1 j xM

Z13 = Z1

Z23 = Z13

Матрицы контурных составляющих, ЭДС и токов

Z

 

Z

 

Z

 

 

 

U

 

 

 

11

 

12

 

13

 

 

 

AB

 

Zk = Z12

Z22

Z23

Ek = UBC

 

 

 

Z23

Z

 

 

 

 

 

 

Z13

33

 

UCA

Необходимо решить систему уравнений

Zk Ik = Ek

( 2.14)

21

Токи в ветвях цепи

IA = Ik1 + Ik3,

IB = Ik2 Ik1,

IC = – Ik2 Ik3

Iab = Ik1,

Ibc = Ik2,

Icz = – Ik3

Показания ваттметров

 

 

 

Uab = UAB + Z1 (IB IA)

Ucb = – UBC + Z1 (IB IC)

P1 = Re(Uab

 

)

P2 := Re(Ucb

 

)

IA

IC

Тепловые потери в треугольнике

Pt := Rе (|Iab|2 + |Ibc|2)

Задача 9. Определение характеристик генератора

 

 

 

Постановка задачи

В конце линии (рис. 2.9) фаза А замкнута на землю. Фазная ЭДС

генератора Е = 20 кВ,

контурные сопротивления: генератора ZГ1 = j9 Ом, ZГ2 =j1

Ом, ZГ0 = j0,5 Ом; линии ZЛ1 = ZЛ2 = j1 Ом, ZЛ0 = j2 Ом; нагрузка ZН1 = j10 Ом,

ZН2 = j2 Ом; сопротивление заземления нулевой точки ZN = j0,5 Ом. Определить:

1)значения токов во всех фазах генератора и нагрузку;

2)напряжение по отношению к земле на зажимах генератора и нагрузку.

 

A

а

O

E,ZГ В

ZЛ b

 

С

c

ZN

 

 

ZН

КЗ

Рисунок 2.9. Схема к задаче 9

Необходимо написать программу определения характеристик генератора. Для решения системы линейных алгебраических уравнений (2.15) воспользоваться методом Гаусса.

В программе предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц (согласовать с руководителем курсовой работы) проверить результаты работы программы.

Математическая модель

В задаче имеет место так называемая поперечная несимметрия. Место несимметрии имитируем введением источников с напряжениями Uа, Ub, Uс, через которые протекают токи IA, IB, IC (рис. 2.10). Указаны напряжения и токи образуют несимметричные системы векторов и могут быть разложены на симметричные составляющие.

22

A

 

а

 

 

I

 

 

IаН

 

E,ZГ

ZЛ

b

 

ZН

В

 

I

I

Сc

I

Uа

Ub

Uс

IсН

ZN

IA

IB

IC

 

 

 

Рисунок 2.10. Схема к задаче 9

В итоге в соответствии с принципом наложения один несимметричный трехфазный круг распадается на три симметричные – прямой, обратной и нулевой последовательностей (рис. 2.11, а, б, в). Упростим схемы рис. 2.11, а, б, в, приведя их к виду рис. (рис. 2.12, а, б, в).

a)

Е IГ1 ZЛ1

ZГ1 U1 I1

 

б)

IН2

в)

 

 

 

IН1

IГ2 ZЛ2

 

IГ0

ZН1

ZГ2 U2

ZН2

ZГ0

 

I2

 

 

 

 

3ZN

ZЛ0

U0

I0

а)

I1 Z1

E1 U1

Рисунок 2.11. Схема к задаче 9

б)

в)

I2

Z2

U2

I0 Z0

U0

Рисунок 2.12. Схема к задаче 9

Параметры схем рис. 2.12 определяются по формулам

Z1 = ZН1 (ZГ1 ZЛ1) ;

ZН1 ZГ1 ZЛ1

Z2 = ZН2 (Z Г2 ZЛ2 );

ZН2 Z Г2 ZЛ2

Z0 = ZГ0 + ZЛ0 + 3ZN ;

23

Е1 =

E (Z Г1 ZЛ1 ) 1

 

(в соответствии с методом двух узлов).

(Z Г1 ZЛ1 ) 1 ZН1 1

На основании второго закона Кирхгофа получаем три уравнения:

U1 + I1 Z1 = Е1,

U2 + I2 Z2 = 0,

 

U0 + I0 Z0 = 0.

Следующие три уравнения составляются по условиям конкретной

несимметрии :

 

 

 

 

IВ = a2 I1 + a I2 + I0 = 0,

 

I

С = a I1 + a2 I2 + I0 = 0,

U

а = U1 + U2 + U0 = 0

Получаем систему из шести уравнений

 

 

U1 + I1 Z1 = Е1,

 

 

 

U2 + I2 Z2 = 0,

 

 

 

 

U0 + I0 Z0 = 0,

 

 

 

 

IВ = a2 I1 + a I2 + I0 = 0,

(2.15)

IС = a I1 + a2 I2 + I0 = 0,

 

 

Uа = U1 + U2 + U0 = 0

Решив систему, получаем симметричные составляющие, с помощью которых находим все искомые в задаче величины.

2.2. Задачи численного интегрирования

Задача 10. Расчет напряженности магнитного поля

Постановка задачи

В устройстве (рис. 2.13) постоянный ток I протекает по трём проводам, как указано на этом рисунке. Радиусы проводов r0 = 1см, расстояние между

проводами d = 70 см.

Требуется рассчитать напряжённости магнитного поля в точках А и В, а также магнитное напряжение между ними. Считая, что А и В являются точками сечения длинных сторон прямоугольной рамки длиной l = 5м и с числом витков W = 200, найти магнитный поток рамки и взаимную индуктивность устройства и рамки.

I = 140 А, d= 0,35,

хА = – 0,05 см; yА = – 0,25 см; хВ = 0,15 см; yB = -0,3 см, 0 4 10 7

Рисунок 2.13. Схема к задаче 10

24

Необходимо написать программу определения напряженности магнитного поля. Для вычисления интегралов в формулах (2.18), (2.19) воспользоваться методами согласно своему варианту. Точность вычисления интеграла согласовать

сруководителем курсовой работы . Для каждого метода вывести значение N.

Впрограмме предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц (согласовать с руководителем курсовой работы) проверить результаты работы программы.

Вариант 1. Методы Симпсона и Чебышева Вариант 2. Методы трапеций и Гаусса Вариант 3. Методы Монте-Карло и Ньютона-Котеса

Математическая модель поставленной задачи

На основании рис. 2.14 для любой точки с координатами x, y выражаем проекцию вектора напряженности магнитного поля на ось x от действия, например, левого провода

cos

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d x 2 y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проекция напряженности поля в этой точке на ось х :

 

H

 

(x, y)

2I

 

 

cos

 

 

 

H

 

(x, y)

I

1

 

2y

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

d x

2 y2

 

x

 

 

 

 

d x 2

y2

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично определяются проекции вектора Н и на ось х и на ось y от действия других проводов.

Hy H

α

Hx

 

 

 

2I

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

Рисунок 2.14. Проекции напряженности поля к задаче 10

 

Тогда результат для Нх и Ну будет следующий

 

 

 

 

 

 

 

I

 

2 y

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

d y

 

 

 

 

Hx(x, у):=

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 2.16)

2

(d x)

2

y

2

(d x)

2

y

2

 

 

 

 

 

 

 

 

 

2

 

2

 

 

3 d y)

x

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

Hy(x, у):=

I

 

2 (d x)

 

d x

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 2.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(d x)2 y2

(d x)2 y2

 

( 3 d y )2 x2

 

 

Полные значения напряженностей у точках А и В:

25

HA :=Hx( xA ,yA )2 Hy( xA ,yA )2

HB :=Hx( xB ,yB )2 Hy( xB ,yB )2

Магнитное напряжение между точками определим в соответствии с соотношением:

B

UM = H dl

A

Интегрирование выполним по горизонтали и вертикали. Тогда

xB

yB

 

UM := Hx( x,yA )dx+

Hx( xB ,y)dy

( 2.18)

xA

yA

 

Магнитный поток, замыкающийся через рамку в Вб и взаимная

индуктивность провода и рамки в Гн

 

 

 

yB

xB

 

 

 

 

 

 

Hx (xB , y)dy Hy (x, yA )dx

 

( 2.19)

Ф l 0

 

 

W Ф

yA

xA

 

 

M

 

 

 

 

( 2.20)

I1

 

 

 

 

 

 

 

 

Задача 11. Расчет напряженности электрического поля

Постановка задачи

Между двумя плоскими электродами напряженность поля изменяется по

 

 

2

 

Ey = 0, Ez = 0.

закону E=Ex = E0· 1 x

2d2

,

 

 

 

 

Расстояние между электродами d=5мм намного меньше размеров пластин,

причем a = 25 см,

b = 1 м,

E0 = 12 кB/см, диэлектрическая проницаемость

диэлектрика = 4.

 

 

 

 

Нужно найти напряжение между электродами, объемную плотность свободного заряда и весь свободный заряд, который есть между электродами.

Необходимо написать программу определения заданных характеристик. Для вычисления интегралов (2.21), (2.22) воспользоваться методами согласно своему варианту. Точность вычисления интеграла согласовать с руководителем курсовой работы. Для каждого метода вывести значение N.

В программе предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц (согласовать с руководителем курсовой работы) проверить результаты работы программы.

Вариант 1. Методы трапеций и Чебышева Вариант 2. Методы Симпсона и Ньютона-Котеса Вариант 3. Методы Монте-Карло и Гаусса

26

Математическая модель поставленной задачи

Состояние поля определяется уравнением Пуассона

2 = –

 

0 .

 

 

 

 

 

 

 

 

Так как

d << а и

d << b, то можно пренебрегать краевым эффектом. Тогда

зависит только от х и

2 =

d2

 

 

 

 

 

 

dx2

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

Напряженность поля E(x)

= –

 

= E0· 1 x

 

2d2

 

,

dx

 

 

 

 

 

 

 

 

 

 

Откуда

d2

= –

dE

=

E0x

= –

 

,

dx2

 

dx d2

 

0

Напряженность между электродами

d

0E0 x d2

U = E( x)dx

 

( 2.21)

0

 

 

 

Суммарный заряд между электродами

 

 

 

d

 

 

Q = dV = ( x)abdx

 

( 2.22)

V

0

 

 

Задача 12. Определение энергии в электрическом поле

 

Постановка задачи

 

В изоляции коаксиального кабеля

(рис. 2.15) есть равномерно

распределенный заряд с объемной плотностью =

= 10 -10 Кл/м3. Внутренний и

внешний радиусы диэлектрического слоя – r1 = 1 см і r2 = 10 см, длина l = 100 м, относительная диэлектрическая проницаемость = 4. Определить энергию, которая запасена в электрическом поле.

Необходимо написать программу определения энергии в электрическом поле. Для вычисления интегралов (2.24), (2.25) воспользоваться методами согласно своему варианту. Точность вычисления интеграла согласовать с руководителем курсовой работы. Для каждого метода вывести значение N.

В программе предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц (согласовать с руководителем курсовой работы) проверить результаты работы программы.

Вариант 1. Методы Монте-Карло и Симпсона Вариант 2. Методы Гаусса и Ньютона-Котеса Вариант 3. Методы трапеций и Чебышева

27

r2

r

Рисунок 2.15. Схема изоляции коаксиального кабеля к задаче 12

Математическая модель поставленной задачи

Задачу решим, проинтегрировав уравнение Лапласа для области r1 < r < r2. Запишем:

(r) = –

 

 

r2 + А1ln(r) + А2;

 

 

 

 

 

 

4 0

 

Е(r) =

 

 

r

A1

.

( 2.23)

2 0

 

 

 

 

r

 

Граничные условия: (r=r2) = 0, Е(r=r1) = 0. Из полученной системы уравнений найдем постоянные интегрирования А1 і А2.

Энергию, которая запасена в электрическом поле, можно найти двумя способами.

В первом способе воспользуемся формулою для энергии заряда dq слоя толщиной dr, расположенного на расстоянии r от оси:

dW = 0,5dq· = 0,5 ·2 rl·dr· ;

r2

 

W = dW = rl (r )dr

( 2.24)

Vr1

Вдругом способе воспользуемся формулой для объемной плотности

энергии электростатического поля: w = 0,5 0Е2.

wdV = 0,5 0 Е2·2 r l·dr = 0Е2· rl·dr,

 

r2

 

W = wdV = 0 rl (E(r ))2dr

( 2.25)

V

r1

 

28

2.3. Обработка экспериментальных данных

Задача 13. Использование метода наименьших квадратов при решении электротехнических задач

Постановка задачи

Вариант 1

В результате опыта холостого хода получена зависимость потребляемой из сети мощности (P0, Вт) от входного напряжения (U1, В) для асинхронного кранового двигателя (АД МТН 111-6), которая хранится в текстовом файле mth111.txt.

U1, В

132

140

150

162

170

180

190

200

211

220

232

240

251

P0, Вт

330

350

385

425

450

485

540

600

660

730

920

1020

1350

Построить линию регрессии, вычислить коэффициент корреляции. Методом наименьших квадратов подобрать зависимость вида P0 a1 a2U12 , с помощью которой найти механические потери pмех P0(0) a1 . Вычислить индекс корреляции для подобранной зависимости. Для линии регрессии и зависимости P0 a1 a2U12 вычислить суммарную квадратичную ошибку, среднюю ошибку в точке и относительную ошибку. Все вычисления проводить в относительных единицах.

В программе исходные данные можно считывать с клавиатуры или из файла mth111.txt. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

Вариант 2

В результате расчета характеристик общепромышленного асинхронного двигателя АД АИУМ225М4 по экспериментальным данным получена табличная зависимость величины тока статора (I1, А) от коэффициента нагрузки (β) при входном напряжении U1=0.9U, которая хранится в текстовом файле aium225m4.txt.

 

х.х. ( )

0,2

0,5

0,8

0,9

1

1,1

1,3

1,5

I1, А

18

21

36

54

59

67

73

88

105

Построить линию регрессии, вычислить коэффициент корреляции. Методом

наименьших

квадратов

подобрать

зависимость вида I

1

a

0

a a

2

a

3 .

 

 

 

 

 

1

2

 

3

 

Вычислить

ожидаемое

значение

тока при коэффициенте

нагрузке

1,17.

Вычислить индекс корреляции для подобранной зависимости. Для линии

регрессии и зависимости I

1

a

0

a a

2

a

3

вычислить суммарную

 

 

1

2

 

3

 

 

квадратичную ошибку, среднюю ошибку в точке и относительную ошибку. Все вычисления проводить в относительных единицах.

В программе исходные данные можно считывать с клавиатуры или из файла aium225m4.txt. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

29

Вариант 3

В результате расчета характеристик общепромышленного асинхронного двигателя АД АИУМ225М4 по экспериментальным данным получена табличная зависимость величины тока статора (I1, А) от коэффициента нагрузки (β) при входном напряжении U1=0,9U, которая хранится в текстовом файле aium225m4.txt.

 

 

х.х.

 

0,2

 

0,5

 

0,8

0,9

1

 

 

 

 

1,1

1,3

1,5

I1, А

 

18

 

21

 

36

 

54

59

67

 

 

 

73

88

105

Вычислить коэффициент корреляции. Методом наименьших квадратов

подобрать

зависимости

вида

I

1

a

0

a a

2

a

3 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

3

 

I1 b0 b1 b2 2 I1

b0 b1 b2 2 .

Вычислить ожидаемое значение тока при

коэффициенте

нагрузках

0,65;

0,84.

Для

подобранных

зависимостей

вычислить индекс корреляции, суммарную квадратичную ошибку, среднюю ошибку в точке и относительную ошибку. Все вычисления проводить в относительных единицах.

В программе исходные данные можно считывать с диалогового окна или из файла aium225m4.txt. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

Вариант 4

Известны экспериментальные кривые нагрева для асинхронного кранового двигателя АД МТН 111-6 при номинальном напряжении, которые хранятся в текстовом файле mth111.txt.

t,мин

2

4

6

8

10

18

20

25

30

35

40

50

60

70

90

110

140

170

200

230

, С

11

20

27

32

36

48

50

56

61

66

70

76

82

87

94

101

108

112

113

113

и напряжении U1=0,9U

t,мин

2

4

6

8

10

14

18

20

25

30

40

60

80

100

120

150

180

210

240

, С

16

25

30

35

39

46

52

55

61

66

77

92

102

108

116

123

126

129

130

 

 

Методом

наименьших

квадратов

подобрать

зависимость

вида

(t)

уст

1 a e t /T1 .

Вычислить продолжительность

переходного процесса

 

 

 

1

 

 

T ln

 

a 1 e 60/T1 . Вычислить ожидаемое значение

t

пп

по формуле

t

пп

уст

 

 

 

 

 

1

 

1

 

 

 

температуры в момент времени 57 мин.

В программе исходные данные можно считывать из файла mth111.txt или с клавиатуры . Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

30

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