Описание работы программы. Укрупненный алгоритм программы приведен в Приложении. Программа использует стандартный модуль CRT, где сосредоточены процедуры и функции управляющие текстовым режимом работы экрана.
В нашей программе модуль используется для перемещения курсора в произвольную позицию экрана при выводе результатов расчета. Программа расчета имеет процедуру SYSTRD для решения системы алгебраических уравнений методом прогонки. Обратимся к укрупненной блок-схеме процедуры SYSTRD и рассмотрим ее работу. Блок 1. Проверяется условие минимального достаточного числа решаемых уравнений.
Процедура используется при nn 2. В противном случае происходит выход из процедуры, остановка основной программы с выдачей соответствующего сообщения подробнее - в следующем разделе. Блок 2. Вычисляется mnn-1 N-1. Передаются в процедуру значения фактических параметров bb1, aa1, dd1 коэффициенты b1, a1, d1, вычисляются xx1, gg1. Блок 3. Организован цикл для внутренних точек от 2 до N-1. Блок 4. Передается в процедуру значения фактического параметра cc1 и в цикле последовательно вычисляются xxi, qqi Блок 5. Передается в процедуру значения фактического параметра ccnn для точки N т.е.cN. Вычисляются xxm, xxnn Блок 6. Организуется цикл в порядке убывания от точки N-1 к 1. Блок 7. В цикле вычисляются значения параметра-переменной xxi, и передаются в вызывающую программу. Ввод данных соответствует блоку 1 общего алгоритма.
Исходными данными являются шаг разбиения h Объмная мощность тепловыделения qv коэффициент теплоотдачи с боковых поверхностей 0 теплофизические параметры нагретой зоны l1, c1, p1 теплофизические параметры слоя теплоизоляции l2, c2, p2 время расчета tmax шаг времени измерения dt количество узлов n. Блок 1. Организован цикл для выбора стержня от 1 до n. Блок 2. Организован цикл для выбора узла в стержне от 1 до n. Блок 3, 5 Проверка совпадения значения выбранного условия решения с заданным значением варианта решения.
При совпадении вычисляются коэффициенты a, b, c, d для выбранного узла выбранного стержня.
При не совпадении переход к проверке следующего совпадения. Блок 4, 6. Вычисление коэффициентов a, b, c, d. Блок 7. Обращение к подпрограмме SYSTRD для решения системы уравнений при разбивке по X и Y соответственно. Блок 9,10. Получение нового распределения температур. Блок 11. Увеличение временного шага tame. Блок 12. Проверка условия окончания измерения. Блок 13. Вывод результатов.