Метод преобразования многолучевой звезды в полный многоугольник (Лабораторная работа № 4), страница 5

Подпрограмма ZHORTKA выполняет свёртку схемы до заданного узла и создает трыангуляваную матрицу коэффициентов и матрицу правых частей системы узловых уравнений (9), причем в трыангуляванай матрице коэффициентов запоминаются только ненулевые элементы. Матрица коэффициентов Y в явной форме подпрограммой не создается, а свёртка схемы выполняется путем соответствующих преобразований информации в списках отраслей схемы. При удаленни каждого узла из этих списков избираются ветви, которые присоединены к удаляемому узлу, затем по формулам (4,5) рассчитываются параметры ветвей полного многоугольника и, наконец,ветви многоугольника вносятся назад в список ветвей. В результате последовательного исключения n-1 узлов (где n - количество линейно-независимых узлов) схема оказывается свернутой до последнего n-го узла, в качестве которого может быть выбран узел короткого замыкания (КЗ).

Входными данными для подпрограммы являются списки с информацией о схеме замещения (массива номеров концов N1K (KH), N2K (KH), проводимостей Y (KH) и токов источников тока J (KH) отраслей, массив номеров узлов NW (KW), количество ветвей KH и количество узлов KW схемы. В массиве номеров узлов узел КЗ должен быть записан последним.  . Выходными данными является информация о системе с трыангуляванай матыцай коэффициентов (массивы: ненулевых элментов трыангуляванай матрицы коэффициентов YHMPZWW (KNNE), номеров вершин многолучевых звезд в удаляемых узлах NWMPZWW (KNNE), правых частей трыангуляванай системы JMPZWW (KW), собственных проводимостей удаляемых узлов YSMPZWW (KW), количества ветвей в удаляемых узлах на момент их удаления KHMPZWW (KW) и число ненулевых элементов в трыангуляванай матрице KNNE). После выполнения обнуления переменных и массива (блок 101) начинается выполнение цикла (конец на метке 32) по удалению узлов схемы. Цикл включает блоки 102-110. В этом цикле используются рабочие массива: KHDNW (KW) - количество ветвей, присоединенных к неудаленым узлам; NWMPZ (KHMPZ) - номера узлов многолучевой звезды (номера вершин многоугольника); YHMPZ (KHMPZ) - проводимости отраслей многолучевой звезды; JHMPZ (KHMPZ) - токи источников тока в областях многолучевой звезды. В рабочих массивах находится информация о параметрах отраслей многолучевой звезды в узле, который удаляется в текущем цикле; на новом цикле информация в рабочих массива обновляется.

Алгоритм преобразования многолучевой звезды в полный многоугольник повторяется в каждом цикле подпрограммы ZHORTKA для каждого удаляемого узла. В блоке 102 для каждого узла схемы определяется количество ветвей, присоединенных к каждому из неудаленых узлов, и запоминаются в массиве KHDNW (KW). В блоке 103 определяется наибольшее количество ветвей, присоединенных к одному или нескольким узлам, оставшихся в списке. В блоке 104 выполняется перестановка номров узлов в массиве NW (KW) таким образом, чтобы в начале массива оказались узлы, к которым подключено наименьшее количество отраслей. Удаляемые узлы избираются из массива NW (KW) по порядку, поэтому следующий удаляемый узел будзе иметь наименьшее количество ветвей, и эквивалентный многоугольник при накладыванни его на непеределанну часть схемы даст наименьший прирост в списках с информацией.

В блоке 105 из списков информации о схеме выбирается информация о ветвях многолучевой звезды и записывается в рабочий массив. Информация об остальных ветвях схемы в списках перемещается на свободные места. В блоке 107 выполняется сложение параллельных ветвей многолучевой звезды, если они в ней есть. В блоке 108 выполняется преобразование многолучевой звезды в полный многоугольник. При этом рассматриваются все возможные сочетания из n по 2 и по формулам (4,5) рассчитываются приросты проводимость и токов. Если в списке информации соответствующая приростом ветвь присутствует, то к существующей ветви добавляются рассчитанные приросты проводимости и тока. Если в списке информации ветвь отсутствует, то в списке информации образуется новая ветвь с проводимостью и током, равными рассчитанный приростом проводимости и тока.В блоке 109 в выходные массивы подпрограммы записываются значения коэффициентов и правые части того уравнения трыангуляванай системы, соответствующие удаленному в данном цикле узлу.