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