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