Предложен метод разбиения рабочего пространства шахтной печи на зоны, являющиеся основой для построения математической модели с сосредоточенными параметрами. Составлена схема эквивалентных химических реакций процесса. Выполнено условное разбиение рабочего пространства на идеализированные зоны: зона окисления, зона восстановления, зона нагрева и зоны запаздывания, позволяющие упростить математическое описание. Разработана структура математической модели процесса переработки свинца в шахтных печах. Выбран метод компьютерного моделирования для исследования поведения системы в динамике на его математической модели. Получена итоговая математическая модель работы шахтной печи свинцовой плавки в виде системы уравнений, учитывающей материальный баланс, тепловой баланс восстановительной зоны, тепловой баланс окислительной зоны и опытные данные. Смоделирована на языке визуального программирования математическая модель переработки свинцового агломерата в шахтных печах, полученной аналитическим способом. Приведены результаты исследования поведения системы в динамике на его математической модели.