3.3.35. Поверхность тел, колеблющихся с ультразвуковой частотой, кажется скользкой на ощупь, а предметы, помещенные на эту поверхность, «плывут» по ней от малейшего приложенного к ним усилия. Объясните это.
Решение
Пусть имеется горизонтальная поверхность с коэффициентом трения скольжения $\mu$, которая колеблется вдоль прямой, проходящей через нее, с большой частотой $\omega$. Указанную прямую назовем $y$, прямую, перпендикулярную ей и лежащую на поверхности, назовем $x$.
Далее, на поверхности в соприкосновении оказывается тело, движущееся в плоскости, параллельной колеблющейся поверхности. Обозначим скорость тела $\vec V$; проекции этой скорости на оси $x,y$ равны $V_x,V_y$ соответственно.
Кроме того, действие на тело усилия моделируется внешней силой $\vec F$, лежащей в плоскости, параллельной колеблющейся поверхности. Ее проекции на оси $x,y$ равны $F_x,F_y$.
Скорость колеблющейся плоскости $\vec v$. Ее проекция на ось $y$ меняется по закону
$$v_y=v_0\sin(\omega t),$$
где $v_0$ есть амплитуда скорости, разумеется ($v_x=0$).
Изобразим качественный рисунок для определения силы трения скольжения $\vec f$, действующей на тело в плоскости его движения. Начнем с момента времени $t=0$.
К задаче $3.3.35$
Сила $\vec f$, как можно показать, направлена вдоль скорости $\vec v-\vec V$ (это скорость плоскости относительно тела!).
Рассмотрим приближение постоянной скорости тела в течение одного периода $T$ колебаний плоскости — то есть будем интересовать законом движения тела в элементах времени, равных $T$. Сам закон движения будем выводить из уравнения изменения импульса тела за время $T$[1,2].
В пределах одного периода начиная с момента $t=0$ сила $\vec f$ будет меняться по направлению. Можно убедиться, что проекция этой силы на ось $y$ будет равна
При условии $V\ll v_0$ зависимость $|f_y|(t)$ достаточно прямоугольная и, конечно, периодическая. Вот пример для $V_x=V_y=0.03v_0$ (Asymptote графика).
К задаче $3.3.35$
Величина $|f_y|$ очень быстро устанавливается на значении $\mu mg$, есть только кратковременные острые снижения до нуля. Случай $V_x\neq V_y$ не отличается принципиально. Красным цветом обозначено, когда проекция $f_y>0$, синим — когда $f_y<0$.
Теперь не будет большой ошибкой записать уравнение изменения импульса по оси $y$ за время $T$:
при рассмотрении $T\equiv dt$ для наблюдателя, не замечающего, что происходит со скоростью тела в течение периода $T$.
Величина $4\mu mg\dfrac{V_y}{v_0\omega}$ эквивалентна силе сопротивления, зависящей от скорости $V_y$. Итак, для указанного наблюдателя по оси $y$ скорость меняется так, как будто вдоль этого направления телу препятствует сила вязкого трения (как по воде).
Движение перпендикулярно оси колебаний
Учитывая предыдущее, проекция силы трения на ось $x$:
Посмотрим на качественный график $|f_x|(t)$, допустим, для $V_x=V_y=0.03v_0$.
К задаче $3.3.35$
Величина $|f_x|$ за период $T$ практически постоянно держится на уровне малом по сравнению с $\mu mg$. Лишь два раза за это время она испытывает кратковременные скачки до значения $\mu mg$. Случай $V_x\neq V_y$ не отличается принципиально.
Уравнение изменения импульса за время $T$ по оси $x$ можно записать так:
$$F_xT+\Delta p_x=m\Delta V_x,\tag{2}$$
где $\Delta p_x$ есть импульс проекции силы $f_x$ за время $T$, равный отрицательной площади под графиком на рисунке зависимости $|f_x|(t)$ на интервале $T$.
Оценивая график указанной зависимости, говорим: практически $\Delta p_x\ll\mu mgT$. Можно представить, что $\Delta p_x$ набирается некоторой эквивалентной постоянной проекцией силы трения на ось $x$, обозначим ее $\tilde f_x$, так что
$$\Delta p_x=\tilde f_xT,$$
и для этой проекции силы ввести эквивалентный коэффициент трения $\tilde \mu$ такой, что
$$\Delta p_x=-\tilde \mu mgT.$$
Тогда из вышесказанного $\tilde\mu\ll\mu$. Для наблюдателя, не замечающего, что происходит со скоростью тела в течении периода $T\equiv dt$, формула $(2)$ переходит во второй закон Ньютона:
$$F_x-\tilde \mu mg=m\frac{dV_x}{dt},$$
где $\tilde\mu mg$ есть эквивалент силы трения скольжения с малым коэффициентом трения. Движение по оси $x$ для указанного наблюдателя напоминает движение по неподвижной поверхности, но с малым коэффициентом трения.
Связанные темы: маятник Капицы, быстрые и медленные движения, вибрация.
[2] Павел Виктор. Решение задачи Савченко 3.3.36. youtube.com.
[3] Блехман И. И. Вибрационная механика. 1994. § 8.1. "Об эффекте вибрационного перемещения, его теории и приложениях".
[4] Мицкевич А.М. «Движение тела по тангенциально колеблющейся поверхности с учетом трения». Акустический журнал, 13, 3, с. 411-416 (1967).
[5] Ю.М. Волченко. Научная графика на языке Asymptote. 2018.
[6] Charles Staats III. Пособие по Asymptote. 2022. (англ.)
Альтернативное понимание условия задачи
Эпиграф: «Вингардиум Левиоса» (Гарри Поттер).
В моём понимании поверхность тела колеблется с ультразвуковой частотой $\omega$ в вертикальном направлении! И если на тело поставить предмет, то это своего рода подбрасывание тела плитой, которая колеблется. В задаче $3.3.30.$ было показано, что предмет будет отрываться от тела на высоте $y=\frac{g}{\omega^{2}}$. Поскольку для ультразвука $\omega$ очень большое, то высота отрыва очень маленькая. Предмет, как бы постоянно парит (левитирует) над телом - между телом и предметом есть воздушная подушка. Через очень быстрое изменение высоты воздушного пространства воздух не успевает выходить из областей сжатия и разрежения, возникает подъёмная сила избыточного давления. Так как паря предмет не касается тела, то нет и силы трения, поэтому тело легко скользит без сопротивления.
Пусть поверхность совершает гармонические колебания:
$$z(t)=A\sin\omega t$$
При малом зазоре $h(t)$ воздух сжимается. Давление можно оценить адиабатически:
где $P_{0}$ - атмосферное давление, $h_{0}$ - средний зазор между предметом и вибрирующей поверхностью, $\gamma$ - показатель адиабаты для воздуха (воздух можно считать двухатомным газом $\gamma=7/5$).
Сила реакции опоры: под$N=mg-F_{под}$, тогда сила трения: трпод$F_{тр}=\mu N=\mu(mg-F_{под})$. Когда подъёмная сила сравняется с силой тяжести, сила трения устремится к 0.
На мой взгляд это и есть объяснение эффекта описанного в задаче. Можно пойти дальше и оценить подъёмную силу.
$$P(t)=P_{0}\left( \frac{h_{0}}{h_{0}-A\sin\omega t} \right)^{\gamma}=P_{0}\left( 1-\frac{A}{h_{0}}\sin\omega t \right)^{-\gamma}$$
Разложим в ряд при $\frac{A}{h_{0}}\ll 1$:
$$P(t)\approx P_{0}\left[ 1+\gamma\frac{A}{h_{0}}\sin\omega t+\frac{\gamma(\gamma+1)}{2}\left( \frac{A}{h_{0}} \right)^{2}\sin^{2}\omega t \right]$$
Усредняя по времени, воспользуемся известными соотношениями: $\left\langle \sin\omega t \right\rangle=0$ и $\left\langle \sin^{2}\omega t \right\rangle=\frac{1}{2}$, тогда:
$$\left\langle P \right\rangle-P_{0}=P_{0}\frac{\gamma(\gamma+1)}{4}\left( \frac{A}{h_{0}} \right)^{2}$$
При $\gamma=7/5$ получим:
$$\left\langle P \right\rangle-P_{0}=\frac{21}{25}P_{0}\left( \frac{A}{h_{0}} \right)^{2}$$