Аннотация. Рассмотрен алгоритм расчета собственных колебаний бурильной колонны в скважине, ось которой задается произвольной пространственной кривой. Бурильная колонна моделируется геометрически нелинейными балочными конечными элементами с использованием модифицированного подхода Крейга-Бемптона. Уравнения движения колонны линеаризуются в окрестности положения равновесия. Для учета влияния жидкости используется полученное авторами S.S. Chen, M.W. Wambsganss, J.A. Jendrzejczyk аналитическое выражение для коэффициентов вязкого сопротивления жидкости и присоединенной массы в случае плоской задачи о гармонических колебаниях кругового цилиндра в трубе, заполненной несжимаемой ньютоновской жидкостью. 135 © Сетевое издание «Нефтегазовое дело». 2018. №5 http://ogbus.ru С учетом модели Чена матрицы масс и демпфирования линеаризованной модели зависят от частоты колебаний. В данной работе приведено решение Чена, анализируются алгоритмы расчета коэффициентов диссипации и присоединенных масс через модифицированные функции Бесселя в условиях плохой обусловленности выражений для этих коэффициентов. Предложен быстрый итерационный алгоритм расчета собственных частот и форм колебаний колонны, а также метод оценки степени демпфирования каждой частоты. Разработанные алгоритмы реализованы в программном комплексе «Универсальный механизм» (ПК УМ).