Моделирование включает такую последовательность операций (рис. 12.6):
1) шаг решения для всех границ ячеек (все тепловые потоки в течение шага принимаются постоянными и рассчитываются по состоянию на начало шага):
2) расчет мощности тепловых потоков через границу двух ячеек по формуле (12.16) и количества теплоты, проходящей через границу за шаг;
3) удаление найденного количества теплоты из ячейки по одну сторону границы и добавление в ячейку за границей;
4) расчет для всех ячеек температуры в конце шага:
5) расчет приращений количества теплоты в каждой ячейке за время шага с учетом ее поступлений через все границы ячейки по формуле (12.17); расчет приращения теплосодержания в каждой ячейке за шаг по формуле ; расчет температуры ячейки в конце шага по формуле
. (12.22)
Рис. 12.6. Моделирование тепловых процессов по явной схеме
Расчет по формуле (12.22) - наиболее сложная операция, поскольку требуется найти приращение температуры по известному приращению теплосодержания, часть которого может быть израсходована на фазовые превращения. Приращение температуры в шве и околошовной зоне даже за малый шаг по времени может быть достаточно большим, в этом случае возможно изменение за шаг входящей в формулу теплоемкости материала. В программе используется процедура, позволяющая решить нелинейное уравнение (12.22) без итераций. Процедура такова (рис. 12.7):
1) имеем исходную температуру после предыдущего шага и зависимость теплоемкости от температуры, заданную в табличной форме с достаточно мелким шагом ;
2) если приращение теплосодержания за шаг , начинаем увеличивать температуру от с шагом , учитывая затраты тепла на нагрев и на фазовые превращения, пока не израсходуем все приращение теплосодержания (назовем эту операцию интегрированием диаграммы вверх);
3) если , начинаем уменьшать температуру от с шагом , учитывая поступление тепла от охлаждения и фазовых превращений (назовем эту операцию интегрированием диаграммы вниз);
4) Последнее значение температуры в конце интегрирования и будет решением – температурой в конце шага .
Рис. 12.7. Интегрирование вверх по диаграмме теплоемкости