Диссипативная функция Рэлея: вывод и применение

Когда маятник раскачивается в воздухе, он постепенно останавливается. Колебания пружины с грузом в вязкой среде затухают. В каждом из этих случаев механическая энергия рассеивается в тепло через силы трения, пропорциональные скорости. Такие силы называют вязкими или силами Стокса: они линейны по скорости и направлены противоположно движению.

Для систематического учёта таких сил в аналитической механике Рэлей ввёл специальную скалярную функцию, которая позволяет добавить диссипацию прямо в уравнения Лагранжа, не выписывая каждую силу трения вручную. Это особенно ценно в задачах с несколькими степенями свободы: достаточно найти одну функцию $\mathcal{F}$ и взять от неё производные - уравнения движения с затуханием готовы. Ниже разберём, что такое функция Рэлея, как её вывести и как с ней работать. Чтобы сразу увидеть, как меняется мощность потерь при разных параметрах системы, покрутите калькулятор ниже.

Что такое диссипативная функция Рэлея

Диссипативная функция Рэлея $\mathcal{F}$ - это скалярная функция обобщённых скоростей, определённая так, что её частные производные по обобщённым скоростям дают обобщённые диссипативные силы:

$$Q_k^{\text{дис}} = -\frac{\partial \mathcal{F}}{\partial \dot{q}_k}.$$

Для системы с вязким трением, где силы сопротивления пропорциональны скоростям, функция Рэлея имеет вид квадратичной формы:

$$\mathcal{F} = \frac{1}{2} \sum_{i,\,j} r_{ij}\,\dot{q}_i\,\dot{q}_j,$$

где $r_{ij}$ - коэффициенты вязкого сопротивления (симметричная матрица). Для простейшего случая одной степени свободы:

$$\mathcal{F} = \frac{1}{2}\,b\,\dot{q}^2,$$

где $b$ - коэффициент вязкого трения (Н·с/м или кг/с). Именно эта формула чаще всего встречается в задачах механики первого и второго курсов.

Рассчитать для своей задачи

Вывод через уравнения Лагранжа

Стандартные уравнения Лагранжа для консервативной системы:

$$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_k} - \frac{\partial L}{\partial q_k} = 0.$$

Если силы трения действуют, то в правой части появляются обобщённые силы $Q_k$:

$$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_k} - \frac{\partial L}{\partial q_k} = Q_k.$$

Для вязкого трения сила $F_{\text{тр}} = -b\,\dot{x}$ в декартовых координатах. Переходя к обобщённым координатам через виртуальную работу $\delta W = F_{\text{тр}}\,\delta x = -b\,\dot{x}\,\delta x$, получаем обобщённую силу:

$$Q_k = -b\,\dot{q}_k.$$

Рэлей заметил, что это в точности $-\partial\mathcal{F}/\partial\dot{q}_k$ при $\mathcal{F} = \frac{1}{2}b\dot{q}^2$. Это ключевое наблюдение: нет нужды вычислять обобщённые силы напрямую через виртуальные перемещения - достаточно построить квадратичную форму $\mathcal{F}$ по скоростям. Таким образом, уравнения Лагранжа с диссипацией записываются компактно:

$$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_k} - \frac{\partial L}{\partial q_k} + \frac{\partial \mathcal{F}}{\partial \dot{q}_k} = 0.$$

Рассчитать для своей задачи

Мощность диссипации и связь с потерями

Ключевое свойство функции Рэлея: скорость рассеяния механической энергии равна удвоенному значению $\mathcal{F}$:

$$P_{\text{дис}} = -\frac{dE}{dt} = 2\,\mathcal{F} = b\,\dot{q}^2.$$

Это следует из теоремы об изменении энергии: полная производная $E = T + U$ по времени равна мощности диссипативных сил:

$$\frac{dE}{dt} = \sum_k Q_k^{\text{дис}}\,\dot{q}_k = -\sum_k \frac{\partial \mathcal{F}}{\partial \dot{q}_k}\,\dot{q}_k = -2\mathcal{F},$$

где последнее равенство использует теорему Эйлера об однородных функциях (для квадратичной $\mathcal{F}$). Это важно: зная $\mathcal{F}$, можно сразу оценить тепловые потери без интегрирования уравнений движения.

Затухающий осциллятор: полный расчёт

Рассмотрим пружинный маятник: груз массой $m$ на пружине с жёсткостью $k$, движущийся в вязкой среде с коэффициентом $b$. Лагранжиан:

$$L = T - U = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2.$$

Функция Рэлея:

$$\mathcal{F} = \frac{1}{2}b\dot{x}^2.$$

Подставляем в уравнение Лагранжа с диссипацией:

$$m\ddot{x} + b\dot{x} + kx = 0.$$

Это и есть стандартное уравнение затухающих колебаний. Вводим параметры: коэффициент затухания $\beta = b/(2m)$ и собственная частота $\omega_0 = \sqrt{k/m}$. При $\beta < \omega_0$ (недодемпфирование) решение:

$$x(t) = A\,e^{-\beta t}\,\cos(\omega t + \varphi_0),$$

где $\omega = \sqrt{\omega_0^2 - \beta^2}$ - частота затухающих колебаний. Средняя мощность потерь за один период:

$$\langle P \rangle = b\,\langle\dot{x}^2\rangle = \frac{1}{2}b\omega^2 A^2 e^{-2\beta t}.$$

Логарифмический декремент и добротность

Из функции Рэлея легко получить важные характеристики затухания. Амплитуда убывает по закону $A(t) = A_0 e^{-\beta t}$, поэтому отношение двух последовательных амплитуд:

$$\frac{A_n}{A_{n+1}} = e^{\beta T},$$

где $T = 2\pi/\omega$ - период. Логарифмический декремент затухания:

$$\vartheta = \ln\frac{A_n}{A_{n+1}} = \beta T = \frac{\pi b}{m\omega}.$$

Добротность системы $Q$ связана с функцией Рэлея через соотношение:

$$Q = \frac{\omega_0}{2\beta} = \frac{m\omega_0}{b} = \frac{\pi}{\vartheta}.$$

Высокая добротность означает малые потери за период: $Q \gg 1$ соответствует $\mathcal{F}_{\text{ср}} \ll T_{\text{ср}}$ (средняя рассеиваемая мощность много меньше средней кинетической энергии, делённой на единицу времени).

Пример расчёта: маятник с вязким трением

Рассмотрим конкретный числовой пример. Дано: груз массой $m = 1$ кг на пружине жёсткостью $k = 16$ Н/м движется в среде с вязким трением $b = 2$ кг/с. Начальная амплитуда $A = 0{,}1$ м.

Шаг 1 - записываем функцию Рэлея:

$$\mathcal{F} = \frac{1}{2} \cdot 2 \cdot \dot{x}^2 = \dot{x}^2 \text{ (Дж)}.$$

Шаг 2 - определяем параметры:

$$\omega_0 = \sqrt{\frac{k}{m}} = \sqrt{\frac{16}{1}} = 4\ \text{рад/с}, \quad \beta = \frac{b}{2m} = \frac{2}{2} = 1\ \text{с}^{-1}.$$

Поскольку $\beta = 1 < \omega_0 = 4$ (недодемпфирование):

$$\omega = \sqrt{\omega_0^2 - \beta^2} = \sqrt{16 - 1} = \sqrt{15} \approx 3{,}873\ \text{рад/с}.$$

Шаг 3 - максимальная скорость равна $v_{\max} = A\,\omega_0 = 0{,}1 \cdot 4 = 0{,}4$ м/с (точнее - через $\sqrt{\beta^2+\omega^2} \approx 4{,}0$ м/с). Максимальная мощность потерь:

$$P_{\max} = b\,v_{\max}^2 = 2 \cdot (0{,}4)^2 = 0{,}32\ \text{Вт}.$$

Шаг 4 - логарифмический декремент и добротность:

$$T = \frac{2\pi}{\omega} \approx \frac{6{,}283}{3{,}873} \approx 1{,}622\ \text{с}, \quad \vartheta = \beta T \approx 1 \cdot 1{,}622 \approx 1{,}62, \quad Q = \frac{\omega_0}{2\beta} = \frac{4}{2} = 2.$$

Малая добротность $Q = 2$ говорит о сильном затухании: уже за несколько периодов амплитуда упадёт до пренебрежимо малого значения. Все эти числа можно проверить и поменять в калькуляторе выше, двигая ползунки.

Обобщение на несколько степеней свободы

Для системы с $n$ степенями свободы и вязким трением в общем случае:

$$\mathcal{F} = \frac{1}{2} \begin{pmatrix} \dot{q}_1 & \cdots & \dot{q}_n \end{pmatrix} \mathbf{R} \begin{pmatrix} \dot{q}_1 \\ \vdots \\ \dot{q}_n \end{pmatrix} = \frac{1}{2}\,\dot{\mathbf{q}}^T \mathbf{R}\,\dot{\mathbf{q}},$$

где $\mathbf{R}$ - симметричная положительно полуопределённая матрица коэффициентов диссипации. В декартовых координатах для системы $N$ частиц:

$$\mathcal{F} = \frac{1}{2}\sum_{i=1}^N b_i\,\dot{\mathbf{r}}_i^2,$$

где $b_i$ - коэффициент трения для $i$-й частицы. Переход к обобщённым координатам осуществляется через матрицу Якоби преобразования координат. Обобщённая диссипативная сила для $k$-й координаты выражается через суперпозицию вкладов всех частиц:

$$Q_k^{\text{дис}} = -\sum_{i=1}^N b_i\,\dot{\mathbf{r}}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_k},$$

и это в точности совпадает с $-\partial\mathcal{F}/\partial\dot{q}_k$, если подставить матричную запись. Такой подход особенно удобен при анализе колебаний стержней, мембран и других распределённых систем, где каждая точка среды вносит вклад в общую диссипацию энергии через локальную вязкость.

Связь функции Рэлея с физическими наблюдаемыми прямая: интегрирование $2\mathcal{F}$ по времени от нуля до бесконечности даёт полную механическую энергию, рассеявшуюся за всё время движения. Для затухающего осциллятора это в точности исходная полная механическая энергия $E_0 = \frac{1}{2}kA^2$, что служит удобной проверкой корректности выбора функции Рэлея в конкретной задаче.

Частые ошибки

• Путаница знака в уравнении Лагранжа. Правильно: $+\partial\mathcal{F}/\partial\dot{q}_k = 0$ в левой части. Если поставить минус, диссипация будет накачивать энергию.
• Подстановка силы трения, а не её модуля. $\mathcal{F}$ всегда положительна; обобщённая диссипативная сила отрицательна ($Q_k^{\text{дис}} = -\partial\mathcal{F}/\partial\dot{q}_k < 0$).
• Применение к кулонову трению. Функция Рэлея применима только к линейному (вязкому) трению. Сухое трение Кулона нелинейно и не допускает такого представления.
• Забыть множитель 2 при вычислении мощности. $P_{\text{дис}} = 2\mathcal{F}$, а не $\mathcal{F}$.
• Игнорирование условия положительной полуопределённости $\mathbf{R}$. Матрица коэффициентов должна быть положительно полуопределённой, иначе $\mathcal{F}$ может принимать отрицательные значения, что физически бессмысленно.

FAQ

Зачем нужна функция Рэлея, если можно просто добавить силу трения? Для одной степени свободы разницы нет. Но при $n > 1$ степенях свободы выписывать все $n$ обобщённых сил трения вручную трудно и легко ошибиться. Функция $\mathcal{F}$ сводит задачу к нахождению одного скаляра, из которого все силы получаются дифференцированием автоматически.

Как связаны функция Рэлея и теорема о флуктуационно-диссипационных соотношениях? Флуктуационно-диссипационная теорема связывает матрицу $\mathbf{R}$ (из функции Рэлея) с матрицей спектральных плотностей тепловых флуктуаций: $S_{ij}(\omega) = 2k_B T\,r_{ij}$. Это позволяет рассчитать тепловой шум в механических системах зная только $\mathcal{F}$.

Можно ли применять функцию Рэлея для нелинейного трения? В общем случае нет. Если сила трения зависит от скорости нелинейно (например, $F \propto \dot{q}^2$), то квадратичная $\mathcal{F}$ не воспроизводит правильную обобщённую силу. Можно формально ввести обобщённую диссипативную функцию, не квадратичную, но она теряет многие удобные свойства функции Рэлея.

Коротко

Диссипативная функция Рэлея $\mathcal{F} = \frac{1}{2}b\dot{q}^2$ - это скалярная функция скоростей, частные производные которой дают обобщённые силы вязкого трения с обратным знаком. Она добавляется в уравнения Лагранжа в виде $+\partial\mathcal{F}/\partial\dot{q}_k$, а мощность диссипации равна $2\mathcal{F}$. Метод позволяет единообразно учитывать затухание в системах с любым числом степеней свободы, переходить от декартовых к обобщённым координатам автоматически и получать логарифмический декремент и добротность из одной формулы.