Если текущий угол отклонения груза $\varphi$, скорость груза $v = l \dot{\varphi}$, сила натяжения нити равна
$$
T = mg \cos \varphi + \frac{m v^2}{l} \approx mg\left(1 - \frac{\varphi^2}{2} \right) + ml \dot{\varphi}^2.
$$
Средний квадрат угла и его производной при гармонических колебаниях выражаются через амплитуду $\alpha$ как
$$
\langle \varphi^2\rangle = \frac{\alpha^2}{2}, \quad
\langle \dot{\varphi}^2\rangle = \frac{\alpha^2 \omega^2}{2},
$$
где $\omega^2 = g/l$ — циклическая частота колебаний.
Отсюда получаем
$$
\langle T \rangle = mg \left(1 - \frac{\langle \varphi^2 \rangle}{2} \right) + ml \langle \dot{\varphi}^2 \rangle = mg \left( 1 + \frac{\alpha^2}{4}\right).
$$
Работа силы натяжения нити равна произведению средней силы натяжения нити на перемещение. Учтем, что при увеличении длины нити перемещение направлено противоположно действию силы, поэтому работа отрицательна. Также можно учесть работу силы тяжести $mg \Delta l$, действующей на груз, тогда полная работа будет равна $A = - mg \Delta l \alpha^2/4$.
Изменение энергии груза равно работе, совершенной над ним внешними силами при изменении длины нити. При этом первый член в работе $- mg \Delta l$ отвечает изменению потенциальной энергии груза в положении равновесия. Мы будем отсчитывать потенциальную энергию от положения равновесия, поэтому в дальнейшем этот член учитывать не нужно. Энергия колебаний выражается через их амплитуду
$$
W = m g l \frac{\alpha^2}{2}.
$$
Скорость изменения энергии
$$
\frac{dW}{dt} = - mg \frac{\alpha^2}{4} \frac{dl}{dt}.$$
Выразим здесь амплитуду через энергию:
$$
\frac{dW}{dt} = -\frac{W}{2l}\frac{dl}{dt}.
$$
Отсюда получаем
$$
\frac{dW}{W} = - \frac{1}{2} \frac{dl}{l}, \quad \ln W = -\frac{1}{2} \ln l + const,
$$
а значит
$$
W \sim l^{-1/2}
$$
Из этого соотношения следует, что при медленном изменении длины нити маятников отношение $W/\omega \approx const$ остается постоянной. Она называется адиабатическим инвариантом.
Амплитуда колебаний определяется энергией
$$
W \sim l \alpha^2,
$$
поэтому
$$
\alpha \sim W^{1/2}l^{-1/2} \sim l^{-3/4}.
$$
Тогда конечная амплитуда
Если электрическое поле постоянно и равно $E$, то движение маятника происходит так же, как и в эквивалентном гравитационном поле с ускорением свободного падения $g_1 = g + qE/m$. При изменении поля частота колебаний $\omega =\sqrt{g_1/l}$ также будет меняться. Поскольку поле переменное его работа за период не будет равна нулю, несмотря на то, что груз маятник вернется в исходную точку. Пусть за период колебаний $T$ электрическое поле изменилось на величину $\Delta E$. Найдем соответствующую работу, она будет равна изменению кинетической энергии груза. Важно, что мы будем рассматривать промежуток, который начинается с прохождения грузом положения равновесия. Тогда вся энергия груза в начале и в конце равна кинетической энергии, а значит работа равна увеличению энергии колебаний.
Зависимость угла отклонения маятника от времени имеет вид $\varphi = \alpha \sin \omega t$, а вертикальная координата маятника $z = l(1- \cos \varphi) \approx l \varphi^2/2$ (ось $z$ направлена вверх).
Тогда работа электрического поля
$$
dA = -q E dz = -q El \varphi d \varphi = - q E l \alpha \sin\omega t \omega \alpha \cos\omega t dt = -q El \omega \alpha ^2 \sin \omega t \cos \omega t dt.
$$
Электрическое поле $E = E_0 + \Delta E\frac{t}{T}$, причем вклад в работу за период дает только переменная часть поля, поэтому полная работа
$$
A = -q l \omega \alpha ^2\frac{\Delta E}{T} \int_0^T dt\, t \sin \omega t \cos \omega t
$$
Вычислим стоящий здесь интеграл:
$$
\int_0^T dt\, t \sin \omega t \cos \omega t = \frac{1}{2}\int_0^T tdt\, \sin 2 \omega t = \frac{1}{2}\int_0^T td \left(- \frac{\cos 2 \omega t}{2 \omega} \right) =
$$
$$
\left.- \frac{t \cos 2\omega t}{ 4 \omega }\right|_0^T + \frac{1}{4 \omega}\int_0^T dt \cos 2 \omega t = - \frac{T}{4 \omega}.
$$
$$
A = \frac{1}{4}q l \alpha ^2\Delta E
$$
Тогда скорость изменения энергии
$$
\frac{dW}{dt} = \frac{A}{T} = \frac{1}{4}q l \alpha ^2 \frac{dE}{dt}.
$$
Энергия колебаний выражается через амплитуду
$$
W = (mg + qE)l \alpha^2/2
$$
$$
\frac{dW}{dt} = \frac{qW}{2(mg +qE)} \frac{dE}{dt}
$$
$$
\frac{dW}{W} = \frac{1}{2} \frac{q dE}{q E + mg}, \quad \ln W = \frac{1}{2} \ln(qE + mg) + const
$$
Таким образом получаем
$$
W \sim \sqrt{g+ qE/m}
$$
то есть опять сохраняется адиабатический инвариант $W/\omega$.
Значит амплитуда
$$
\alpha \sim \left( g + \frac{qE}{m}\right)^{-1/4}
$$
Cредняя сила натяжения нити равна $T = m_2 g (1 + \alpha_0^2/4) - m_2 a_0$, здесь первое слагаемое — это среднее натяжение нити из предыдущей части, а последний вклад возникает из-за ускорения второго груза вдоль нити.
Второй закон Ньютона для первого груза (после усреднения по движению второго):
$$
m_1 a_0 = - m_1 g + m_2 g (1 + \alpha_0^2/4) - m_2 a_0.
$$
Подставляя значения масс, получаем
$$
(m_1+m_2) a_0 = - (m + \delta m)g + mg + mg \alpha_0^2/4 = - \delta m g + m g \alpha_0^2/4.
$$
Все члены в правой части малы, поэтому слева можно пренебречь отличием $m_1$ от $m$, тогда ускорение
Среднее ускорение обращается в 0 при
Запишем уравнение движения груза $m_1$, его координату обозначим $z$, начальное значение координаты равно нулю, ось $z$ направлена вверх.
$$
m_1 \ddot{z} = - m_1 g + T.
$$
В первом приближении левый груз не влияет на движение правого. тогда закон движения правого груза
$$
\varphi(t) = \alpha_0 \cos \omega t, \omega = \sqrt{\frac{g}{l_0}},
$$
тогда сила натяжения нити
$$
T = m_2( g \cos \varphi + l_0 \dot{\varphi}^2) - m_2 \ddot{z}= m\left( g - g\frac{\varphi^2}{2}+ l _0\dot{\varphi}^2 \right) - m_2\ddot{z};$$
$$
T = mg\left(1 - \frac{\alpha_0^2}{2} \cos^2 \omega t + \alpha_0^2 \sin^2 \omega t \right) - m \ddot{z}.
$$
Уравнение движения
$$
2m \ddot{z} = - (m + \delta m)g + mg\left(1 - \frac{\alpha_0^2}{2} \cos^2 \omega t + \alpha_0^2 \sin^2 \omega t \right).
$$
Подставим значение $\delta m = \delta m_0$:
$$
2m \ddot{z}= mg \alpha_0^2\left(- \frac{1}{4} -\frac{1}{2} \cos^2 \omega t + \sin^2 \omega t\right) = \frac{3 m g \alpha_0 ^2}{4} \left( 1 - 2 \cos^2 \omega t\right)=- \frac{3 m g \alpha_0 ^2}{4}\cos 2 \omega t.
$$
Решение этого уравнения при заданных начальных условиях
$$
z = \frac{3 g \alpha_0^2}{32 \omega^2}\left( \cos2 \omega t -1\right).
$$
Смещение будет максимальным, когда $\cos \omega t = -1$, при этом
$$
z = - \frac{3 g \alpha_0^2}{16 g/l_0} = - \frac{3}{16} l_0 \alpha_0^2,
$$
знак - говорит о том, что груз сдвинется вниз. Видим, что смещение груза действительно мало, поэтому движение правого груза не меняется.
Снова запишем уравнения движения для левого груза, усредненные по периоду. Пусть в текущий момент времени амплитуда колебаний правого груза равна $\alpha(t)$, тогда
$$
\ddot{z} = - \frac{\delta m}{2m} g + \frac{1}{8} g \alpha^2(t).
$$
В процессе движения длина нити правого груза меняется, при этом согласно A4 меняется и амплитуда колебаний
$$
\alpha(t) = \alpha_0\left( \frac{l_0}{l(t)}\right)^{3/4}.
$$
Текущая длина нити правого груза выражается через смещение левого груза $l(t) = l_0 + z$. Тогда уравнение движение принимает вид
$$
\ddot{z} = - \frac{\delta m}{2m} g + \frac{1}{8} g \alpha_0^2 \left( \frac{l_0}{l_0 + z}\right)^{3/2}
$$
Запишем $\delta m= \delta m_0 + \Delta m$ и разложим правую часть уравнения до первого порядка по $z$:
$$
\ddot{z} = -\frac{\delta m_0}{2m}g - \frac{\Delta m}{2m} g + \frac{1}{8} g \alpha_0^2 \left( 1 - \frac{3}{2} \frac{z}{l_0}\right) = - \frac{\Delta m}{2m} g - \frac{3g \alpha_0^2}{16 l_0} z.
$$
Это уравнение гармонических колебаний с частотой
$$
\Omega = \alpha_0\sqrt{\frac{3 g}{16 l_0}}
$$
Из-за постоянного слагаемого слагаемого в правой части решение уравнения — сумма постоянного смещения и гармонического колебания. С учетом того, что в начальный момент времени смещение равно нулю, закон движения имеет вид
$$
z = - \frac{\Delta m}{m} \frac{8 l_0}{3 \alpha_0^2} (1 - \cos \Omega t) = -B(1 - \cos \Omega t).
$$
Коэффициент перед косинусом — амплитуда колебаний.
Мы получили уравнение гармонических колебаний, раскладывая уравнение движения по $z/l_0$. Приближение работает, если амплитуда колебаний мала $B \ll l_0$, то есть
$$
\frac{\Delta m}{m} \ll \alpha_0^2.
$$
Начальная энергия колебаний правого груза
$$
W_0 = m g l_0 \frac{\alpha_0^2}{2}.
$$
Если в некоторый момент времени длина нити правого груза равна $l$, энергия колебаний в этот момент времени
$$
W = W_0 \frac{l_0^{1/2}}{l^{1/2}}.
$$
При этом изменение потенциальной энергии в гравитационном поле
$$
\Delta U = \delta m g (l - l_0).
$$
Запишем закон сохранения энергии, с учетом того, что в крайних точках кинетическая энергия вертикального движения грузов равна нулю:
$$
W_0 = W + \delta m g (l - l_0) = W_0 \frac{l_0^{1/2}}{l^{1/2}} + \delta m g (l - l_0).
$$
Введем обозначения
$$
x^2 = \frac{l}{l_0}, \quad \xi = \frac{W_0}{\delta m g l_0} = \frac{m \alpha_0^2 }{2 \delta m},
$$
тогда получим уравнение на $x$:
$$
\xi \left( 1 - \frac{1}{x}\right) = x^2 -1.
$$
У этого уравнения всегда есть решение $x=1$, отвечающее начальному положению грузов. Преобразуем:
$$
\xi = \frac{x^2-1}{1-1/x} = x(x+1).
$$
Корни этого уравнения
$$
x = \pm \sqrt{\xi + 1/4} - 1/2.
$$
Имеют смысл только положительные решения, поэтому дальше рассматриваем только положительный знак перед корнем.
Ему отвечает длина нити
$$
l_1 = l_0 \left(\sqrt{\frac{m \alpha_0^2 }{2 \delta m} + \frac{1}{4}}- \frac{1}{2}\right)^2.
$$
При $\xi >2$ $l_1 >l_0$ (то есть $ \alpha_0^2 > 4 \delta m/m$) и значение $l_1$ — максимальное значение длины нити, тогда минимальное значение $l_0$. Если же $\xi <2$, $l_1$ — минимальное значение, а максимальное значение $l_0$.
Заметим, что граничное значение $\xi = 2$ как раз отвечает значению $\delta m/m = \alpha_0^2/4$, при котором левый груз в целом не движется. Если масса левого груза больше критической, он движется вниз, а длина правой нити уменьшается. Если же масса меньше критической, левый груз будет двигаться вверх, а длина нити будет увеличиваться.
Пусть амплитуда колебаний заряда на конденсаторе в некоторый момент времени равна $Q$, тогда энергия колебаний в контуре равна
$$
W = \frac{Q^2}{2 C},
$$
где емкость конденсатора с расстоянием между пластинами $d$ равна
$$C = \frac{\varepsilon_0 S}{d}.$$
Если в данный момент времени заряд пластин конденсатора $q$, электрическое поле в нем равно
$$
E = \frac{q}{\varepsilon_0 S},
$$
они притягиваются друг к другу с силой
$$
F = \frac{q^2}{2 \varepsilon_0 S},
$$
а средняя сила притяжения равна
$$
\langle F \rangle = \frac{Q^2}{4 \varepsilon_0 S}.
$$
Тогда если за время, много большее периода колебаний, расстояние между пластинами изменили на величину $\Delta d \ll d$, при этом совершили работу
$$
A = \langle F \rangle \Delta d = \frac{Q^2}{4 \varepsilon_0 S} \Delta d = \frac{C}{2 \varepsilon_0 S} W \Delta d = \frac{1}{2} W \frac{\Delta d}{d}.
$$
Таким образом, скорость изменения энергии колебаний
$$
\frac{dW}{dt} = \frac{1}{2} \frac{d}{dt} d,
$$
откуда получаем
$$
\ln W = \frac{1}{2}\ln d + const, \quad W \sim d^{1/2} \sim C^{-1/2}.
$$
Такая зависимость энергии напряжения от емкости отвечает адиабатическому инварианту $W/\omega = const$, так как $\omega \sim C^{-1/2}$. Инвариант имеет такой же вид, как и в механике.
Энергия выражается через амплитуду тока:
$$
W = \frac{LI^2}{2},
$$
поэтому $I \sim W^{1/2}$, а значит