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

Результаты численно-аналитического решения уравнения движения неоднородного твэла сравниваются с решением по методу конечных элементов (ANSYS). Рассматриваются перспективы использования метода в исследовании динамики реактора НЕПТУН.