Курсовой ЭТФ
.pdfВ программе предусмотреть считывание исходных данных с клавиатуры или из текстового файла. Результаты вычислений сохранить в текстовом файле 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 |
|
а |
|
|
IAГ |
|
|
IаН |
|
E,ZГ |
ZЛ |
b |
|
ZН |
В |
|
IBГ |
IbН |
Сc
ICГ |
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.9U1н, которая хранится в текстовом файле 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,9U1н, которая хранится в текстовом файле 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,9U1н
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