2 Программа для расчёта режима токов короткого замыкания
В данном разделе представлен текст программы для расчёта токов короткого замыкания приведённой цепи, в приложении А. Блок схема для исходной программы представлена в приложении В. Схема цепи представлена в приложении С.
Пояснения к тексту программы: для начала возьмём исходные данные из таблицы топологии. Затем зададим количество узлов, ветвей и введём в схему точку короткого замыкания. После этого отключим исходную ветвь и приведём схему к одному уровню напряжения. Далее зададим вектор – столбец ЭДС. Поскольку мы добавили три ветви и один узел, увеличим размерность матриц соединений, проводимостей и переберем в цикле все ветви. Далее заполним матрицы g и A. Добавим в программу вектор – столбец источников тока и выполним расчёт режима токов короткого замыкания.
Заключение
Расчетно-графическая работа включает выполнение теоретического и практического задания. Практическое задание заключается в разработке программы на языке MATLAB, обеспечивающей расчет начального значения периодической составляющей тока трехфазного короткого замыкания в трехфазной электрической сети переменного тока.
При выполнении данной расчетно-графической работы был составлен текст программы. Был рассмотрен теоретический вопрос о методе узлового напряжения для эквивалентных преобразований.
Оформление расчетно-графической работы соответствует стандарту АГТУ СТО 01.04-2005 «Работы студентов. Общие требования и правила оформления» [2].
Список использованных источников
1) Нейман Л.Р. Теоретические основы электротехники: Часть 1/ Л. Р. Нейман, К. С. Демирчян; Под ред. Л. Р. Неймана. – Ленинград: Изд-во «Энергия» Ленинградское отделение, 1967. – 523с.
2) АГТУ СТО 01.04-2005. Работы студентов. Общие требования и правила оформления. – Архангельск.: Изд-во АГТУ, 2006. – 80 с.
Приложения приложение а
(обязательное)
Текст программы
Navvet=97;
lkz=0.34;
rd=0.001;
% исходные данные - таблица топологии
% №вет №уз1 №уз2 R E J
Topol = importdata('C:\Users\fpe-2-1\Desktop\Топология.txt', ' ');
% количество ветвей
n = size (Topol,1);
% количество узлов
m = max(max(Topol(1:n,3:4)));
% Ввод в схему точки короткого замыкания
Nom = find (Topol(:,1)==Navvet);
Topol (n+1,:) = [0 0 Topol(Nom,3) m+1 1 Topol(Nom,6)*lkzTopol(Nom,7)*lkz 0 0 Topol(Nom,10)];
Topol (n+2,:) = [0 0 Topol(Nom,4) m+1 1 Topol(Nom,6)*(1-lkz) Topol(Nom,7)*(1-lkz) 0 0 Topol(Nom,10)];
Topol (n+3,:) = [0 0 0 m+1 1 0 rd 0 0 Topol(Nom,10)];
%Отключим исходную ветвь
Topol(Nom,6) = 1E6;
%Приведение схемы к одному уровню напряжения
z = Topol(:,6)+1i*Topol(:,7);
Zpriv = z.*(Topol(Nom,10)./Topol(:,10)).^2;
% вектор-столбецЭДС
e = Topol(:,8)+1i*Topol(:,9);
Epriv = e.*Topol(Nom,10)./Topol(:,10);
%Поскольку мы добавили три ветви и один узел, увеличим размерность матриц
m = m+1;
n = n+3;
% матрица соединений
A = zeros (m,n);
% матрица проводимостей g
g = zeros (n);
% переберем в цикле все ветви, заполним g и A
fori=1:n
g(i,i) = Zpriv(i)^(-1);
% заполняемматрицу A
ifTopol (i,3)>0
A(Topol(i,3),i) = 1;
end
ifTopol(i,4)>0
A (Topol(i,4),i) = -1;
end
end
% вектор-столбец источников тока
J = zeros(n,1);
%Расчет режима
guz = A*g*A';
Juz = -A*(g*Epriv+J);
Uuz = guz\Juz;
Ivet = g*(A'*Uuz+Epriv)+J;
I = Ivet.*Topol(:,10)/Topol(Nom,10)/sqrt(3);