Блок массой m скользит по горизонтальной поверхности с переменным коэффициентом трения μ(x); как построить модель движения, оценить потерю механической энергии и выделить условия, при которых блок остановится или войдёт в режим периодического движения при наличии внешней периодической силы?
Уравнение движения (модель) - Для кинетического трения (Coulomb) с пространственно‑зависимым коэффициентом μ(x)\mu(x)μ(x) и внешней силой Fext(t)F_{\mathrm{ext}}(t)Fext(t): mx¨=−μ(x) mg sgn(x˙)+Fext(t).
m\ddot x = -\mu(x)\,m g\,\operatorname{sgn}(\dot x)+F_{\mathrm{ext}}(t). mx¨=−μ(x)mgsgn(x˙)+Fext(t).
Здесь sgn(x˙)=±1\operatorname{sgn}(\dot x)=\pm1sgn(x˙)=±1 при x˙≠0\dot x\neq0x˙=0. Для численной реализации удобно регуляризовать sgn\operatorname{sgn}sgn через tanh(x˙/ε)\tanh(\dot x/\varepsilon)tanh(x˙/ε). Потеря механической энергии - Закон изменения кинетической энергии: ddt(12mx˙2)=−μ(x) mg ∣x˙∣+Fext(t)x˙.
\frac{d}{dt}\Big(\tfrac12 m\dot x^2\Big) = -\mu(x)\,m g\,|\dot x| + F_{\mathrm{ext}}(t)\dot x. dtd(21mx˙2)=−μ(x)mg∣x˙∣+Fext(t)x˙.
- Для случая без внешней силы (Fext=0F_{\mathrm{ext}}=0Fext=0) интегрирование по координате даёт зависимость энергии от пройденного пути: 12mx˙2(x)=12mx˙02−mg∫x0xμ(s) ds.
\tfrac12 m\dot x^2(x)=\tfrac12 m\dot x_0^2 - m g\int_{x_0}^{x}\mu(s)\,ds. 21mx˙2(x)=21mx˙02−mg∫x0xμ(s)ds.
Отсюда потеря кинетической энергии при фазовом перемещении [x0,x][x_0,x][x0,x] равна mg∫x0xμ(s) dsm g\int_{x_0}^{x}\mu(s)\,dsmg∫x0xμ(s)ds. Условие остановки (без внешней силы) - Точкой остановки xstopx_{\mathrm{stop}}xstop является решение 12mx˙02=mg∫x0xstopμ(s) ds.
\tfrac12 m\dot x_0^2 = m g\int_{x_0}^{x_{\mathrm{stop}}}\mu(s)\,ds. 21mx˙02=mg∫x0xstopμ(s)ds.
- Если интеграл по полуоси бесконечно большой, стоп произойдёт обязательно. Если же mg∫x0∞μ(s) ds<12mx˙02,
m g\int_{x_0}^{\infty}\mu(s)\,ds < \tfrac12 m\dot x_0^2, mg∫x0∞μ(s)ds<21mx˙02,
то остаточная кинетическая энергия 12mx˙∞2=12mx˙02−mg∫x0∞μ(s) ds>0\tfrac12 m\dot x_\infty^2=\tfrac12 m\dot x_0^2 - m g\int_{x_0}^{\infty}\mu(s)\,ds >021mx˙∞2=21mx˙02−mg∫x0∞μ(s)ds>0 и тело не остановится (уйдёт на бесконечность). - В реальности нужно учитывать статическое трение μs(x)\mu_s(x)μs(x): если при достижении x˙=0\dot x=0x˙=0 внешние силы меньше μsmg\mu_s m gμsmg, тело закрепляется (прилипание). Условия периодического движения при внешней периодической силе - Для внешней силы периода TTT баланс энергии за период в установившемся режиме: 1T∫0TFext(t)x˙(t) dt=1T∫0Tμ(x(t)) mg ∣x˙(t)∣ dt.
\frac{1}{T}\int_0^T F_{\mathrm{ext}}(t)\dot x(t)\,dt = \frac{1}{T}\int_0^T \mu(x(t))\,m g\,|\dot x(t)|\,dt. T1∫0TFext(t)x˙(t)dt=T1∫0Tμ(x(t))mg∣x˙(t)∣dt.
Периодическое движение (установившийся периодический режим или limit cycle) возможно если средняя подводимая мощность равна средней диссипации. - Амплитуда и форма Fext(t)F_{\mathrm{ext}}(t)Fext(t) определяют режимы: - Если maxt∣Fext(t)∣<minxμs(x) mg\max_t |F_{\mathrm{ext}}(t)| < \min_x \mu_s(x)\,m gmaxt∣Fext(t)∣<minxμs(x)mg, то тело может оставаться в состоянии «прилипания» (нет движения). - Если амплитуда превышает локальную статическую границу, возникают фазы скольжения/прилипания (stick–slip). Такое движение может быть периодическим (с периодом, равным TTT или кратным ему) при совпадении входной работы и средних потерь. - При отсутствии возвратной силы (пружины) возможен также средний дрейф (ненулевая средняя скорость) — если за период суммарная работа FextF_{\mathrm{ext}}Fext больше суммарных потерь и тело не ограничено по пространству. Для получения замкнутого периодического движения часто требуется наличие пространственно‑ограничивающей силы (пружины) или геометрических границ. Практические замечания для анализа - Интегрируйте численно уравнение с учётом сдвига sgn\operatorname{sgn}sgn; учтите μs(x)\mu_s(x)μs(x) для режимов прилипания. - Для поиска периодических решений используйте Poincaré‑отображение (сечения в моменты t=nTt=nTt=nT) и проверку условия средней мощности. - Для аналитических оценок полезно рассматривать энерго‑баланс за период и локальные условия скольжения: ∣Fext(t)∣>μs(x)mg|F_{\mathrm{ext}}(t)|>\mu_s(x)m g∣Fext(t)∣>μs(x)mg. Кратко: модель — уравнение mx¨=−μ(x)mgsgn(x˙)+Fext(t)\;m\ddot x=-\mu(x)m g\operatorname{sgn}(\dot x)+F_{\mathrm{ext}}(t)mx¨=−μ(x)mgsgn(x˙)+Fext(t);\; энергия убывает по закону ddt(12mx˙2)=−μmg∣x˙∣+Fextx˙\frac{d}{dt}(\tfrac12 m\dot x^2)=-\mu m g|\dot x|+F_{\mathrm{ext}}\dot xdtd(21mx˙2)=−μmg∣x˙∣+Fextx˙;\; остановка определяется условием 12mx˙02≤mg∫x0∞μ(s) ds\tfrac12 m\dot x_0^2 \le m g\int_{x_0}^{\infty}\mu(s)\,ds21mx˙02≤mg∫x0∞μ(s)ds (с учётом μs\mu_sμs для прилипания); периодический режим возможен при равенстве средней подводимой и рассеиваемой мощности за период.
- Для кинетического трения (Coulomb) с пространственно‑зависимым коэффициентом μ(x)\mu(x)μ(x) и внешней силой Fext(t)F_{\mathrm{ext}}(t)Fext (t):
mx¨=−μ(x) mg sgn(x˙)+Fext(t). m\ddot x = -\mu(x)\,m g\,\operatorname{sgn}(\dot x)+F_{\mathrm{ext}}(t).
mx¨=−μ(x)mgsgn(x˙)+Fext (t). Здесь sgn(x˙)=±1\operatorname{sgn}(\dot x)=\pm1sgn(x˙)=±1 при x˙≠0\dot x\neq0x˙=0. Для численной реализации удобно регуляризовать sgn\operatorname{sgn}sgn через tanh(x˙/ε)\tanh(\dot x/\varepsilon)tanh(x˙/ε).
Потеря механической энергии
- Закон изменения кинетической энергии:
ddt(12mx˙2)=−μ(x) mg ∣x˙∣+Fext(t)x˙. \frac{d}{dt}\Big(\tfrac12 m\dot x^2\Big) = -\mu(x)\,m g\,|\dot x| + F_{\mathrm{ext}}(t)\dot x.
dtd (21 mx˙2)=−μ(x)mg∣x˙∣+Fext (t)x˙. - Для случая без внешней силы (Fext=0F_{\mathrm{ext}}=0Fext =0) интегрирование по координате даёт зависимость энергии от пройденного пути:
12mx˙2(x)=12mx˙02−mg∫x0xμ(s) ds. \tfrac12 m\dot x^2(x)=\tfrac12 m\dot x_0^2 - m g\int_{x_0}^{x}\mu(s)\,ds.
21 mx˙2(x)=21 mx˙02 −mg∫x0 x μ(s)ds. Отсюда потеря кинетической энергии при фазовом перемещении [x0,x][x_0,x][x0 ,x] равна mg∫x0xμ(s) dsm g\int_{x_0}^{x}\mu(s)\,dsmg∫x0 x μ(s)ds.
Условие остановки (без внешней силы)
- Точкой остановки xstopx_{\mathrm{stop}}xstop является решение
12mx˙02=mg∫x0xstopμ(s) ds. \tfrac12 m\dot x_0^2 = m g\int_{x_0}^{x_{\mathrm{stop}}}\mu(s)\,ds.
21 mx˙02 =mg∫x0 xstop μ(s)ds. - Если интеграл по полуоси бесконечно большой, стоп произойдёт обязательно. Если же
mg∫x0∞μ(s) ds<12mx˙02, m g\int_{x_0}^{\infty}\mu(s)\,ds < \tfrac12 m\dot x_0^2,
mg∫x0 ∞ μ(s)ds<21 mx˙02 , то остаточная кинетическая энергия 12mx˙∞2=12mx˙02−mg∫x0∞μ(s) ds>0\tfrac12 m\dot x_\infty^2=\tfrac12 m\dot x_0^2 - m g\int_{x_0}^{\infty}\mu(s)\,ds >021 mx˙∞2 =21 mx˙02 −mg∫x0 ∞ μ(s)ds>0 и тело не остановится (уйдёт на бесконечность).
- В реальности нужно учитывать статическое трение μs(x)\mu_s(x)μs (x): если при достижении x˙=0\dot x=0x˙=0 внешние силы меньше μsmg\mu_s m gμs mg, тело закрепляется (прилипание).
Условия периодического движения при внешней периодической силе
- Для внешней силы периода TTT баланс энергии за период в установившемся режиме:
1T∫0TFext(t)x˙(t) dt=1T∫0Tμ(x(t)) mg ∣x˙(t)∣ dt. \frac{1}{T}\int_0^T F_{\mathrm{ext}}(t)\dot x(t)\,dt = \frac{1}{T}\int_0^T \mu(x(t))\,m g\,|\dot x(t)|\,dt.
T1 ∫0T Fext (t)x˙(t)dt=T1 ∫0T μ(x(t))mg∣x˙(t)∣dt. Периодическое движение (установившийся периодический режим или limit cycle) возможно если средняя подводимая мощность равна средней диссипации.
- Амплитуда и форма Fext(t)F_{\mathrm{ext}}(t)Fext (t) определяют режимы:
- Если maxt∣Fext(t)∣<minxμs(x) mg\max_t |F_{\mathrm{ext}}(t)| < \min_x \mu_s(x)\,m gmaxt ∣Fext (t)∣<minx μs (x)mg, то тело может оставаться в состоянии «прилипания» (нет движения).
- Если амплитуда превышает локальную статическую границу, возникают фазы скольжения/прилипания (stick–slip). Такое движение может быть периодическим (с периодом, равным TTT или кратным ему) при совпадении входной работы и средних потерь.
- При отсутствии возвратной силы (пружины) возможен также средний дрейф (ненулевая средняя скорость) — если за период суммарная работа FextF_{\mathrm{ext}}Fext больше суммарных потерь и тело не ограничено по пространству. Для получения замкнутого периодического движения часто требуется наличие пространственно‑ограничивающей силы (пружины) или геометрических границ.
Практические замечания для анализа
- Интегрируйте численно уравнение с учётом сдвига sgn\operatorname{sgn}sgn; учтите μs(x)\mu_s(x)μs (x) для режимов прилипания.
- Для поиска периодических решений используйте Poincaré‑отображение (сечения в моменты t=nTt=nTt=nT) и проверку условия средней мощности.
- Для аналитических оценок полезно рассматривать энерго‑баланс за период и локальные условия скольжения: ∣Fext(t)∣>μs(x)mg|F_{\mathrm{ext}}(t)|>\mu_s(x)m g∣Fext (t)∣>μs (x)mg.
Кратко: модель — уравнение mx¨=−μ(x)mgsgn(x˙)+Fext(t)\;m\ddot x=-\mu(x)m g\operatorname{sgn}(\dot x)+F_{\mathrm{ext}}(t)mx¨=−μ(x)mgsgn(x˙)+Fext (t);\; энергия убывает по закону ddt(12mx˙2)=−μmg∣x˙∣+Fextx˙\frac{d}{dt}(\tfrac12 m\dot x^2)=-\mu m g|\dot x|+F_{\mathrm{ext}}\dot xdtd (21 mx˙2)=−μmg∣x˙∣+Fext x˙;\; остановка определяется условием 12mx˙02≤mg∫x0∞μ(s) ds\tfrac12 m\dot x_0^2 \le m g\int_{x_0}^{\infty}\mu(s)\,ds21 mx˙02 ≤mg∫x0 ∞ μ(s)ds (с учётом μs\mu_sμs для прилипания); периодический режим возможен при равенстве средней подводимой и рассеиваемой мощности за период.