Разработка алгоритма и программы расчета электрической цепи
Методика расчета электрической цепи и особенности использования библиотеки стандартных подпрограмм. Описание и структура алгоритма подпрограммы вычисления определенного интеграла методом прямоугольников, угла сдвига фазы между током и напряжением.
Рубрика | Программирование, компьютеры и кибернетика |
Вид | курсовая работа |
Язык | русский |
Дата добавления | 17.05.2017 |
Размер файла | 240,0 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Размещено на http://www.allbest.ru/
Курсовая работа
Разработка алгоритма и программы расчета электрической цепи
Введение
подпрограмма алгоритм электрический цепь
В настоящее время скорость и качество выполнение технологических процессов тесно связано с применением электронной техники, преимущественно электронных вычислительных машин. Поэтому при подготовке высококвалифицированных специалистов отдельное время должно отводиться на изучение ЭВМ. Инженер должен не только уметь работать на ЭВМ, но и знать основы программирования на алгоритмических языках, уметь применять эти навыки при решении различных задач, моделировании энергетических процессов и т.д.
Задачей курсового проекта является углубление знаний студента в области Фортрана, а также приобретения навыков алгоритмизации и моделирования различных процессов. Это позволит студентам применять средства вычислительной техники в дальнейшем, а также, с помощью навыков алгоритмизации, позволит упростить свою работу.
1. Постановка задачи
Требуется составить программу, позволяющую:
а) сопоставить изменение токов в любых частях схемы, вычисленные мгновенные значения использовать для вычисления действующих значений соответствующих токов;
б) сопоставить изменение напряжения на зажимах схемы UАД и напряжений на любых конденсаторе (UCx) и индуктивности (ULx) исследуемой цепи, вычисленные мгновенные значения использовать для вычисления действующих значений соответствующих напряжений;
в) исследовать влияние изменения произвольного параметра Q от значения Q1 до Q2 на действующие значения токов и напряжений на элементах из пункта б.
Закон изменения напряжения на зажимах схемы AD:
.
Данные для расчёта варианта №46 представлены в таблице 1.1. Исходная расчётная схема имеет вид, представленный на рисунке 1.1.
Таблица 1.1. Исходные данные
1 |
R1 |
R2 |
R3 |
R4 |
R5 |
C1 |
C2 |
C3 |
C4 |
|
2 |
Ом |
Ом |
Ом |
Ом |
Ом |
мкФ |
мкФ |
мкФ |
мкФ |
|
3 |
20 |
40 |
500 |
300 |
50 |
4500 |
900 |
1500 |
900 |
1 |
L1 |
L2 |
L3 |
U |
F |
Q |
Q1 |
Q2 |
|
2 |
мГн |
мГн |
мГн |
В |
Гц |
||||
3 |
15 |
44 |
130 |
220 |
50 |
R5 |
50 |
1500 |
Рисунок 1.1. Расчётная схема электрической цепи (вариант 46)
Рисунок 1.2. Расчётная схема электрической цепи (вариант 33)
Данные для расчёта варианта №33 представлены в таблице 1.2. Исходная расчётная схема имеет вид, представленный на рисунке 1.2.
Таблица 1.2. Исходные данные
1 |
R1 |
R2 |
R3 |
R4 |
R5 |
C1 |
C2 |
C3 |
C4 |
|
2 |
Ом |
Ом |
Ом |
Ом |
Ом |
мкФ |
мкФ |
мкФ |
МкФ |
|
3 |
60 |
0 |
70 |
0 |
0 |
200 |
0 |
2500 |
0 |
1 |
L1 |
L2 |
L3 |
U |
F |
Q |
Q1 |
Q2 |
|
2 |
мГн |
мГн |
мГн |
В |
Гц |
||||
3 |
0 |
41 |
0 |
110 |
200 |
U |
110 |
10 |
2. Методика расчета электрической цепи
Индуктивное сопротивление катушек индуктивности L1:
Индуктивное сопротивление катушек индуктивности L2:
Индуктивное сопротивление катушек индуктивности L3:
Емкостное сопротивление конденсаторов С1:
Емкостное сопротивление конденсаторов С2:
Емкостное сопротивление конденсаторов С3:
Емкостное сопротивление конденсаторов С4:
Рассчитываем сопротивление отдельных участков схемы.
Сопротивления 1-7 ветвиУ
Сопротивление участка EС:
Сопротивление участка BE:
Полное сопротивление цепи AD:
Максимальное значение тока всей схемы:
Максимальное значение тока в седьмой ветви схемы:
Так как через последовательно подключенные элементы протекает один и тот же ток, то амплитудное значение напряжения на участке AC:
Значения амплитудных токов в ветвях 2-7 рассчитываются с помощь закона Ома в виде:
Амплитудное значение напряжения на конденсаторе C1:
Амплитудное значение напряжения на катушке L1:
Находим углы сдвига в неразветвленной части цепи и ее ветвях. Для нахождения этих углов определяется вещественная (Re) и мнимая (Im) части комплексного сопротивления искомого участка. В зависимости от этих величин углы сдвига фаз между током и напряжением в ветвях равны:
Закон изменения тока в электрической цепи и ее ветвях имеет вид:
Законы изменения напряжения на конденсаторе C1 и катушке L1:
Действующие значения токов в цепи и ее ветвях находим по формуле:
Действующие значения напряжений в цепи и ее ветвях находим по формуле:
2.1 Использование библиотеки стандартных подпрограмм
Для реализации вывода всех графиков в программе была применена подпрограмма GRAFIK
Стандартная подпрограмма GRAFIK имеет девять параметров:
V, W, X, Y, Z, R - массивы значений исследуемых функций;
N - размерность массивов;
K - количество позиций для вывода графиков (если вывод осуществляется на широкую печать (АЦПУ), то K128, если на принтер, то K80);
NU - номер логического устройства для вывода графиков.
Вызов этой подпрограммы осуществляется следующим оператором: CALL GRAFIK (V, W, X, Y, Z, R, N).
3. Описание алгоритма
3.1 Алгоритм главной программы
Описание переменных.
Ввод исходных данных (численных значений элементов электрической цепи и др.) из файла ZADANIE2.DAT для 46 варианта и ZADANIE1.DAT для 33 варианта.
Описание массивов и переменных.
Присвоение и расчет констант (р, фазовая частота, период).
Расчет комплексных сопротивлений элементов, ветвей и цепи.
Расчет токов и напряжений.
Нахождение углов сдвига фаз токов и напряжений с помощью подпрограммы FAZ.
Вычисление мгновенных значений токов и напряжений в ветвях 1-7 и на зажимах исследуемых элементов, формирование массивов мгновенных значений.
Вывод мгновенных значений токов и напряжений в виде таблицы в файл REZULTATY.DAT.
Расчет действующих значений токов и напряжений с помощью подпрограммы T1.
Вывод действующих токов и напряжений в файл REZULTATY.DAT.
Вывод графиков мгновенных значений токов и напряжений с помощью подпрограмм grafik и grafik1 соответственно.
Расчет переменных значений, для определения влияния изменения параметра Q на элементы C1 и L1 для 46 варианта и на элементы C3 и L2 для 33 варианта.
Расчет мгновенных токов и напряжений на исследуемых элементах.
Вычисление действующих значений токов и напряжений на элементах C1 и L1 для 46 варианта и C3 и L2 для 33 варианта.
Вывод графиков действующих значений токов и напряжений в зависимости от параметра Q на элементах C1 и L1 для 46 варианта и на элементах C3 и L2 для 33 варианта, а также на зажимах АД, используя подпрограмму grafik1. Вывод таблицы действующих значений тока и напряжения на элементах в зависимости от влияния параметра Q.
3.2 Алгоритм подпрограммы вычисления определённого интеграла методом прямоугольников
Описание переменных.
Вычисляется H = T/N, где H - шаг интегрирования; T - период;
N - число разбиений периода синусоидальной функции.
Очистка ячеек, где будут накапливаться значения.
Вычисление суммы членов массива M начиная с первого и заканчивая N элементом.
Вычисление значения определённого интеграла по известной формуле.
Конец (возврат в точку вызова).
3.3 Алгоритм подпрограммы вычисления угла сдвига фазы между током и напряжением
Описываем переменные.
Проверяем значения переменных.
Вычисляем угол сдвига фаз исходя из условий.
Конец (возврат в точку вызова).
3.4 Алгоритм подпрограммы вычисления минимального и максимального значения функции
Минимальному и максимальному значению присваиваем первые значения первого массива.
Сравниваем значения массивов с минимум и максимум, и заменяем в случае выполнения условия.
Конец (возврат в точку вызова).
4. Таблица идентификаторов
Обозначения параметров кода программы (приложение А) сведем в таблицу идентификаторов (таблица 2).
Таблица 2. Таблица идентификаторов
Обозначение |
Определение |
|
R1R2 R3 R4 R5 C1 C2 C3 C4 L1 L2 L3 U Q1 Q2 F W T XL1 XL2 XL3 XC1 XC2 XC3 XC4 ZADZBCZECZBEZ1 Z2 Z3 Z4 Z5 Z6 Z7 I1am, I6am I2am I3am I4am I5am UCXam ULXam ID1, ID6 ID2 ID3 ID4 ID5 UDC UDL I1, I2, I3, I4, I5, I6UC1, UC2, UC3, UC4 UL1, UL2, UL3 I1D, I2D, I3D, I4D, I5D, I6D IM1, IM2, IM3, IM4, IM5UC1D, UC1M-UC4D, UC4M UL2D, UL2MFII1-FII6FIUC1-FIUC4FIUL1-FIUL3POG-POG7 MICXQ MILXQ MUQ MUCX MULX ICXQD ILXQD UCXQD ULXQD UadQD |
Сопротивление на резисторе R1Сопротивление на резисторе R2Сопротивление на резисторе R3Сопротивление на резисторе R4Сопротивление на резисторе R5Емкость на конденсаторе C1Емкость на конденсаторе C2Емкость на конденсаторе C2Емкость на конденсаторе C4Индуктивность на катушке L1Индуктивность на катушке L2Индуктивность на катушке L3Напряжение цепи АDНачальное значение U, R5Конечное значение U, R5Частота Циклическая частота Период Индуктивное сопротивление катушки L1Индуктивное сопротивление катушки L2Индуктивное сопротивление катушки L3Емкостное сопротивление конденсатора C1Емкостное сопротивление конденсатора C2Емкостное сопротивление конденсатора C3Емкостное сопротивление конденсатора C4Сопротивление участка цепи АDСопротивление участка BCСопротивление участка ECСопротивление участка BEСопротивление ветви 1 Сопротивление ветви 2 Сопротивление ветви 3 Сопротивление ветви 4 Сопротивление ветви 5 Сопротивление ветви 6 Сопротивление ветви 7 Амплитудное значение тока на участке АД Амплитудное значение тока на участке 2 Амплитудное значение тока на участке 3 Амплитудное значение тока на участке 4 Амплитудное значение тока на участке 5 Амплитудное значение напряжения на СXАмплитудное значение напряжения на LXДействующее значение тока на участке АД Действующее значение тока в ветви 2 Действующее значение тока в ветви 3 Действующее значение тока в ветви 4 Действующее значение тока в ветви 5 Действующее напряжение на СXДействующие напряжение на LXТоки в соответствующих ветвяхНапряжения на соответствующих конденсаторахНапряжения на соответствующих индуктивностяхДействительные части токов Мнимые части токов Действительные и мнимые части напряжений на соответствующих конденсаторахДействительные и мнимые части напряжений на соответствующих индуктивностяхУглы сдвига фаз соответствующих токов I1-I6Углы сдвига фаз соответствующих напряжений UC1-UC4Углы сдвига фаз соответствующих напряжений UL1-UL3Погрешности вычислений Мгновенный ток через СХ при изменении параметра QМгновенный ток через LХ при изменении параметра QМгновенное напряжение на ADМгновенное напряжение на СХ при изменении параметра QМгновенное напряжение на LХ при изменении параметра QДействующее значение тока через СХ при изменении параметра QДействующее значение тока через LХ при изменении параметра QДействующее значение напряжения на СХ при изменении параметра QДействующее значение напряжения на LХ при изменении параметра QДействующее значение напряжения на AD |
5. Анализ результатов
При выполнении отладки программы (приложение А) были выявлены некоторые синтаксические и логические ошибки. Синтаксические ошибки находились на этапе трансляции программы, их исправление заключалось во внесении в программу необходимых изменении и дополнений. Логические ошибки определялись при выполнении проверочных расчетов. Для этого вводились в программу дополнительные операторы, позволяющих находить значения некоторых параметров электрической цепи, которые легко проверяемых по законам Ома и Кирхгофа и другим известным методам расчёта электрических цепей. Для получения значений этих параметров вводилась отладочная печать.
Результаты расчётов представлены в приложении Б.
Для проверки программы был выполнен параллельны расчет поставленной задачи с помощью MathCAD. Небольшое расхождение результата объясняется различным округлением результатов, однако можно утверждать о правильности решения поставленной задачи.
Заключение
При выполнении данной курсовой работы была составлена программа расчета электрической цепи на языке FORTRAN, которая отлажена на ЭВМ, получены результаты расчёта, представленные в виде таблиц и графиков.
В процессе программирования разработана подпрограмма для расчета действующих значений синусоидальных величин токов и напряжений по их мгновенным значениям.
Составленная программа может быть использована для исследования параметров идентичной электрической цепи. Такое ограниченное использование является недостатком данной программы.
В целом выполнение курсовой работы позволило углубить знания по программированию в целом, а в частности на языке FORTRAN.
Список использованных источников
1 Фурунжиев Р.И. Вычислительная техника ? практикум. - М.: Высшая школа, 1985.
2 Методические указания по выполнению курсовой работы для студентов специальностей 1-43 01 01 «Электрические станции», 1-53 01 04 «Автоматизация и управление энергетическими процессами» / сост. В.А. Булат, Е.А. Дерюгина, П.И. Климкович. - Минск: БНТУ, 2007.
3 Бессонов Л.А. Теоретические основы электротехники. Электрические цепи. - М.: Высшая школа, 1978.
4 Павловец В.В. Информатика: Программирование на Фортране. - Минск: Асконто, 2006.
5 [Электронный ресурс] / Франц Ульрих Теодор Эпинус /. - Электрон. дан. - Режим доступа http://www.persons-info.com/persons/KHEFNER-ALTENEK_Fridrikh_fon/. - Загл. с экрана. - Яз. рус.
Размещено на Allbest.ru
Подобные документы
Вычисление физических параметров реальной электрической цепи посредством преобразования её к эквивалентной. Схема алгоритма программы и ее разработка на языках программирования СИ и С++, результаты расчета зависимостей эквивалентных сопротивлений.
курсовая работа [19,9 K], добавлен 15.10.2010Разработка программы, которая по заданной самостоятельно функции будет выполнять интегрирование методом прямоугольников. Блок-схема алгоритма вычисления интеграла (функция rectangle_integrate). Экспериментальная проверка программы, ее текст на языке C.
курсовая работа [232,0 K], добавлен 27.05.2013Метод хорд решения нелинейных уравнений. Вычисление интеграла методом Симпсона. Процесс численного решения уравнения. Окно программы расчета корней уравнения методом хорд. Алгоритм вычисления интеграла в виде блок-схемы. Выбор алгоритма для вычислений.
курсовая работа [832,6 K], добавлен 24.07.2012Средства Delphi для разработки Windows приложений. Математическая формулировка задачи, описание программы вычисления определенного интеграла по формуле левых прямоугольников. Руководство пользователя, методика испытаний продукта. Листинг программы.
курсовая работа [178,1 K], добавлен 14.11.2010Описание математической модели определения тока в электрической цепи с помощью решения системы алгебраических уравнений методом Гаусса. Описание и разработка блок-схемы программы. Ввод данных задачи, составление программы и анализ результатов решения.
контрольная работа [231,8 K], добавлен 15.08.2012Разработка программы на языке С++ по определению величин и направлений токов в ветвях электрической цепи с использованием метода Гаусса. Блок-схема алгоритма. Контрольный расчет с помощью электронных таблиц Excel, используя метод обратной матрицы.
курсовая работа [30,3 K], добавлен 10.11.2010Защита электрооборудования и сетей от тока короткого замыкания и перегрузок. Разработка программного обеспечения для расчета номинала предохранителя в электрической цепи. Выбор языка программирования. Требования к составу и параметрам технических средств.
курсовая работа [233,2 K], добавлен 31.03.2018Законы электрических цепей, порядок и методы их расчета. Разработка программы на языке программирования Borland C++ Builder 5.0 для анализа разветвленных электрических цепей с использованием матричного метода. Алгоритм решения задачи и описание его работы
курсовая работа [211,5 K], добавлен 08.10.2012Особенности метода численного интегрирования обыкновенных дифференциальных уравнений. Расчет переходного процесса в нелинейной электрической цепи, вызванного ее включением или отключением. Метод численного интегрирования Рунге-Кутта с переменным шагом.
отчет по практике [740,1 K], добавлен 10.10.2011Написание программы вычисления сопротивления электрической цепи, состоящей из двух параллельно и двух последовательно соединенных сопротивлений. Схема машинного алгоритма по условию задачи. Применение операций при написании программ на языке C/C++.
контрольная работа [17,3 K], добавлен 09.11.2010