Нормальные координаты играют ключевую роль в описании малых колебаний сложных механических систем. Они позволяют разложить движение системы на независимые гармонические моды, каждая из которых колеблется с собственной частотой. Рассмотрение нормальных координат упрощает анализ динамики и существенно облегчает вычисления для систем с множеством степеней свободы.
Рассмотрим механическую систему с n степенями свободы, координаты которой обозначены q1, q2, …, qn. Пусть система находится в устойчивом равновесии при qi = 0 для всех i. Малые отклонения от равновесия обозначим δqi = qi.
Для малых колебаний потенциальная энергия V системы может быть разложена в ряд Тейлора около положения равновесия, оставляя только члены второго порядка:
$$ V \approx V_0 + \frac{1}{2} \sum_{i,j=1}^n \frac{\partial^2 V}{\partial q_i \partial q_j} \bigg|_{q=0} q_i q_j = V_0 + \frac{1}{2} \sum_{i,j=1}^n K_{ij} q_i q_j, $$
где Kij — элементы матрицы жесткости.
Кинетическая энергия системы в линейном приближении записывается как
$$ T = \frac{1}{2} \sum_{i,j=1}^n M_{ij} \dot{q}_i \dot{q}_j, $$
где Mij — элементы матрицы масс, которая в общем случае может быть не диагональной.
Таким образом, уравнения движения для малых колебаний имеют вид
$$ \sum_{j=1}^n M_{ij} \ddot{q}_j + \sum_{j=1}^n K_{ij} q_j = 0, \quad i = 1, \dots, n. $$
Это система линейных дифференциальных уравнений второго порядка с постоянными коэффициентами.
Идея нормальных координат заключается в том, чтобы найти линейное преобразование
$$ q_i = \sum_{k=1}^n a_{ik} Q_k, $$
такое, что новые координаты Qk описывают независимые гармонические осцилляторы, т.е. кинетическая и потенциальная энергии в этих координатах становятся диагональными:
$$ T = \frac{1}{2} \sum_{k=1}^n \dot{Q}_k^2, \quad V = \frac{1}{2} \sum_{k=1}^n \omega_k^2 Q_k^2, $$
где ωk — собственные частоты нормальных мод.
Для нахождения матрицы aik необходимо решить обобщённую задачу о собственных значениях:
det (K − ω2M) = 0.
Решения ωk2 определяют частоты нормальных колебаний, а собственные векторы этой задачи ak = (a1k, a2k, …, ank) задают форму соответствующей нормальной моды.
Каждая нормальная мода характеризуется тем, что все координаты системы колеблются с одной и той же частотой ωk и фиксированными относительными амплитудами, заданными вектором ak.
$$ E = \sum_{k=1}^n \left( \frac{1}{2} \dot{Q}_k^2 + \frac{1}{2} \omega_k^2 Q_k^2 \right). $$
aiTMaj = δij.
$$ \ddot{Q}_k + \omega_k^2 Q_k = 0, \quad Q_k(t) = A_k \cos (\omega_k t + \phi_k). $$
Две массы m1 и m2, соединённые тремя пружинами (концы закреплены), имеют два нормальных колебания:
Kak = ωk2Mak.
aiTMaj = δij.
$$ q_i(t) = \sum_{k=1}^n a_{ik} Q_k(t). $$