ПОИСК Статьи Чертежи Таблицы Программная реализация расчета стационарных средних температур из "Применение ЭВМ для решения задач теплообмена " Рассмотрим процедуру формирования матрицы А и столбца Т. Сначала двумерный массив А и одномерный массив Т обнуляются, а затем производится расчет их ненулевых элементов путем последовательного суммирования отдельных членов, входящих в формулы 1.26)—(1.28). Организация этой процедуры суммирования зависит от используемого способа описания теплового взаимодействия между элементами системы. [c.23] Для примера поясним способ описания теплового взаимодействия между телами. В случае, когда каждое тело взаимодействует со всеми остальными, все отличны от нуля, и для описания тепловых связей можно использовать двумерный массив, элементами которого являются а . Однако обычно на практике число взаимодействующих пар тел значительно меньше максимально возможного значения, н поэтому большинство тепловых проводимостей равно нулю. С целью экономии машинной памяти и сокращения объема исходных данных целесообразно приводить информацию только об отличных от нуля ст . В приведенной программе эта информация задается с помощью двух одномерных массивов IJ и S1.I. В массиве IJ длиной 2 N1J, где NIJ — число связей между телами, парами записаны номера взаимодействующих тел элементы IJ (2 М — 1), IJ (2 М) указывают номера тел, участвующих в М-м взаимодействии. Массив SIJ длиной NIJ содержит значения тепловых проводимостей о , ooi-ветствующих этим тепловым связям. [c.23] Совершенно аналогичным образом представлена информация о тепловых связях 0J между телами и теплоносителями с помощью массивов IL и SIL, а также о связях oJt между телами и средами с.-помощью массивов 1 К и SIK. [c.23] Для формирования матрицы А и столбца свободных членов Т организуются пять циклов по телам по связям между телами (o j/) по связям между телами и теплоносителями (ст] ) по связям между телами и средами (aj ) и по потокам теплоносителя между объемами Gmi, Gft(). В каждом цикле последовательно выбирается связь, определяются номера взаимодействующих тел и объемов, находятся номера строк и столбцов, определяющие коэффициенты матрицы и столбца свободных членов, в которые вносит вклад данная тепловая связь, и, наконец, к ранее вычисленному значению коэффициента матрицы прибавляется вклад от рассматриваемого взаимодействия. [c.26] Далее в коэффициенты матрицы А с указанными номерами добавляются слагаемые о]Тfi, rJ (1 —ft) согласно уравнениям (1.26). (1.27). [c.27] После формирования матрицы А и столбца Т выполняется обращение к стандартной подпрограмме GELG, выходным параметром которой является массив Т. В этом массиве первые Л т элементов соответствуют температурам тел Т,-, а затем парами расположены температуры теплоносителей i/ . Эти температуры выводятся на печать. Если в подпрограмме GELG параметр ошибки IER принимает значение, отличное от нуля, то печатается соответствующее сообщение об ошибке. [c.27] Вернуться к основной статье