ПОИСК Статьи Чертежи Таблицы Программная реализация расчета нестационарных средних температур из "Применение ЭВМ для решения задач теплообмена " Отличие неявных схем для одного линейного уравнения и для системы уравнений состоит в том, что разностная схема (1.35) разрешалась в явном виде относительно ы + , а в данном случае мы имеем систему Nj линейных алгебраических уравнений для определения N,r значений сеточной функции == + , u f,. .., и + . [c.43] Таким образом расчет по неявной схеме Эйлера сводится к решению на каждом шаге по времени системы линейных уравнении (1.66), которое может быть выполнено с помощью какой-либо стандартной подпрограммы. В рассматриваемой задаче матрица А является симметричной, так как согласно (1.67) a,j- = aji = —оТТ, и поэтому используется подпрограмма GELS (см. 1.3). [c.44] Однако если или i зависят от температур, то на каждом шаге / приходится вычислять их новые значения на основе температур u . и, следовательно, снова формировать матрицу А и решать систему уравнений. Поэтому выше рассмотрена структура программы, которая пригодна и для таких нелинейных задач, решаемых путем пересчета матрицы на основе температур предыдущего шага по времени. [c.48] Если начальный шаг в подпрограмме R KGS уменьшался вдвое И раз (IHLF = 11), то управление передается в головную программу, расчеты прекращаются и печатается соответствующее сообщение-В этом случае следует в исходных данных изменить начальное значение тага (TAU) или погрешность (PRMT (4)) и повторить расчет. [c.50] Вернуться к основной статье