- •«Череповецкий государственный университет» Инженерно-технический институт
- •Технология математического моделирования металлургических процессов курс лекций
- •1. Методология разработки математических моделей металлургических процессов
- •1.1. Роль моделей при разработке и освоении новых металлургических технологий и агрегатов
- •1.2. Методология разработки математических детерминированных моделей металлургических процессов
- •Элементы системного подхода
- •Р Рис. 1.2. Структурная схема процесса разработки модели и моделирования объекта ис. 1.3. Условное представление объекта управления
- •2.2. Примеры проблемных ситуаций
- •2. От объекта до формулировки задач исследования
- •2.1. Основные понятия
- •2.3. Пример выполнения этапа разработки «От объекта до формулировки задач исследования»
- •3.Физическое описание объекта
- •Источники информации и пример физического описания
- •Применение системного подхода к формализации физического описания
- •Р ис.3.2. Структурная схема системы процесса вакуумной обработки
- •4. Математическое описание процесса тепломассопереноса в металлургических агрегатах
- •4.1. Виды балансовых уравнений
- •Некоторые упрощения уравнения теплопереноса
- •Выбор начального и граничных условий
- •4.3.1. Выбор начального условия
- •4.3.2. Выбор граничных условий
- •4.4. Учет выделения тепла кристаллизации при затвердевании двойных сплавов и чистых металлов
- •5. Математическое описание гидродинамических процессов в объеме жидкой стали металлургических агрегатов.
- •Вынужденная конвекция жидкости в ограниченном объеме
- •5.2 Свободная конвекция жидкости в ограниченном объеме
- •6. Разработка численной модели
- •6.1. Метод конечных разностей и пример применения
- •Второй производной по координате
- •6.2. Алгоритмы решения конечно-разностной задачи теплопроводности в цилиндре
- •7. Тестирование компьютерной модели тепловых процессов
- •7.1. Введение в тестирование
- •Тестирование алгоритма решения задачи теплопроводности при граничных условиях третьего рода
- •Вывод критериев подобия для граничных условий третьего рода
- •Тестирование алгоритма решения
- •Тестирование алгоритма решения задачи теплопроводности при теплообмене излучением с окружающей средой
- •Вывод критериев подобия для граничных условий излучением
- •Тестирование алгоритмов решения
- •Дополнительные методы настройки алгоритмов численного решения
- •8. Проверка адекватности и адаптация модели
- •Р ис. 8.1. Кривые охлаждения оси (1) и поверхностного слоя (2):
- •9. Выбор исходных данных для моделирования
- •10. Исследование объекта с помощью модели
- •Р ис. 10.3. Схема к расчету количества вариантов моделирования:
- •11. Решение задач исследования
- •Уравнение теплопереноса
- •Уравнение массопереноса
- •Уравнения гидродинамики
- •Уравнение сплошности
- •Точное решение задачи теплопроводности в охлаждаемой пластине при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в пластине, нагреваемой излучением от окружающей среды
- •Точное решение задачи теплопроводности в охлаждаемом цилиндре при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в цилиндре, нагреваемой излучением от окружающей среды
- •Литература
Р ис. 8.1. Кривые охлаждения оси (1) и поверхностного слоя (2):
а – схема расположения термопар, s- расстояние между термопарами
щиты
термопар. Методическую ошибку можно
оценить экспериментальным или
теоретическим путем и учесть в результатах
эксперимента. Тогда погрешность
определяется суммированием отдельных
составляющих и составляет по литературным
данным менее 1 %. В этом случае абсолютная
погрешность измерений, например, при
Тmax=1000°C
составляет
.
В результате с учетом абсолютной
погрешности можно на рис. 8.1 нарисовать
полосу (заштрихована) значений T1±15
,
внутри которой находится кривая (1).
Затем проводят моделирование процесса
затвердевания и охлаждения плоской
симметрично охлажденной заготовки
толщиной 2S
при граничном условии
рода.
При этом в качестве температуры
поверхности принимают значения
температуры в т.1 (рис. 8.1) Т2
=f(t).
Выбор исходных данных для моделирования
производят, используя результаты
собственных исследований, либо обзора
литературных данных (раздел 9).
В результате моделирования получают кривую охлаждения оси заготовки и сравнивают ее с экспериментальной кривой 1 на рис. 8.1. Если расчетная кривая попадает в заштрихованную полосу на рис. 8.1, то модель считается адекватной. Если максимальное отличие результатов моделирования от экспериментальных данных составляет величину δ %:
,
то говорят, что погрешность моделирования составляет δ %.
Если δ>25 %, то вместо модели целесообразно использовать инженерную методику.
При 1≤δ≤25 модель можно использовать, если результаты моделирования устраивают исследователя на данном этапе изучения объекта.
В противном случае модель необходимо адаптировать к объекту. В зависимости от сроков разработки модели возможны три способа адаптации:
параметрический;
параметрический с дополнительной экспериментальной проверкой;
развитие модели.
Если времени недостаточно для развития модели, то выбирают первый способ. Среди исходных данных выбирают параметр, изменением которого добиваются адекватности модели объекта. Следует отметить, что модель будет выдавать недостоверные данные при отклонении от условий эксперимента, по результатам которого выполнена адаптация модели.
При представлении модели необходимо указать, что адаптация модели производилась по такому-то параметру и привести значение параметра, при котором достигнута адекватность.
Во втором способе проводят адаптацию модели как в первом. Затем для других условий проводят дополнительные эксперименты, которые используют для проверки адекватности модели при значении параметра, установленного ранее. Если адекватность будет доказана на основе дополнительных экспериментов, то достоверность модели существенно возрастет.
В последнем способе возвращаются к этапу моделированию – разработке математического описания; находят самое «грубое» допущение; разрабатывают математическое описание, устраняющее это допущение; повторяют все последующие этапы моделирования, включая проверку адекватности. Если адекватность опять не достигнута, то повторяют все изложенное в данном способе.
