Шарик из стали и тонкостенная стеклянная сфера одинакового радиуса катятся с вершины наклонной плоскости: стальной шар приходит к низу быстрее — какие физические причины (момент инерции, деформация, трение качения, аэродинамика) могут объяснить разницу, как спроектировать эксперимент для количественного разделения вкладов каждого фактора и какие погрешности следует учесть?
Коротко — какие факторы и почему, затем план эксперимента для разделения вкладов и список основных погрешностей. 1) Физические причины (кратко) - Момент инерции: для чистого качения без проскальзывания ускорение зависит только от отношения I/(mr2)I/(mr^2)I/(mr2): a=gsinα1+I/(mr2).
a=\frac{g\sin\alpha}{1+I/(mr^2)}. a=1+I/(mr2)gsinα.
Для твёрдого шара I=25mr2⇒a=57gsinαI=\tfrac{2}{5}mr^2\Rightarrow a=\tfrac{5}{7}g\sin\alphaI=52mr2⇒a=75gsinα, для тонкостенной оболочки I=23mr2⇒a=35gsinαI=\tfrac{2}{3}mr^2\Rightarrow a=\tfrac{3}{5}g\sin\alphaI=32mr2⇒a=53gsinα. Разница 5/73/5=2521≈1,19\dfrac{5/7}{3/5}=\dfrac{25}{21}\approx1{,}193/55/7=2125≈1,19 — тонкостенная сфера теоретически ~19% медленнее при идеальном качении. - Деформация / роллинговое сопротивление (hysteresis): упругие потери в материале и на контактном пятне создают дополнительный момент MrM_rMr (или эквивалентную тангенциальную силу). Это снижает ускорение и сильнее действует на более мягкие/комплайантные материалы (тонкое стекло/подложка). - Трение качения / сцепление: чтобы было чистое качение, требуется достаточный статический коэффициент трения. При частичном проскальзывании кинетические потери замедляют шар сильнее. - Аэродинамика: сопротивление воздуха Fdrag=12ρCdAv2F_{drag}=\tfrac12\rho C_d A v^2Fdrag=21ρCdAv2 (или при малых Re — линейный закон). Для лёгкой пустой стеклянной оболочки аэродинамика влияет сильнее, т.к. тормозящее усилие не уменьшается при малой массе. 2) Простейшая модель с учётом основных потерь Для поступательного ускорения с учётом аэродинамики и дополнительного тормозного момента MrM_rMr: a=gsinα−12mρCdAv2−Mrmr1+I/(mr2).
a=\frac{g\sin\alpha-\dfrac{1}{2m}\rho C_d A v^2-\dfrac{M_r}{mr}}{1+I/(mr^2)}. a=1+I/(mr2)gsinα−2m1ρCdAv2−mrMr.
Условие отсутствия проскальзывания: Ff=I/(mr2)1+I/(mr2) mgsinα≤μsmgcosα⇒μs≥I/(mr2)1+I/(mr2)tanα.
F_f=\frac{I/(mr^2)}{1+I/(mr^2)}\,m g\sin\alpha\le\mu_s m g\cos\alpha \quad\Rightarrow\quad \mu_s\ge\frac{I/(mr^2)}{1+I/(mr^2)}\tan\alpha. Ff=1+I/(mr2)I/(mr2)mgsinα≤μsmgcosα⇒μs≥1+I/(mr2)I/(mr2)tanα. 3) Эксперимент — порядок действий для количественного разделения вкладов - Подготовка: - Использовать одинаковый радиус rrr для обоих шаров и жёсткую ровную направляющую платы (чтобы минимизировать вариации поверхности). - Измерить массу mmm и геометрию; оценить/измерить III (по формуле по геометрии или экспериментально через крутильный маятник). - Основные измерения: - Снимать движение высокоскоростной камерой, отслеживать центр и угловую скорость ω(t)\omega(t)ω(t) (метка на шаре) → вычислять v(t)v(t)v(t) и отношение v/(rω)v/(r\omega)v/(rω) для обнаружения проскальзывания. - Измерять время пробега ttt для набора углов α\alphaα и многократно повторять. - Выделение вкладов: 1. Проверка чистого качения: если v/(rω)≈1v/(r\omega)\approx1v/(rω)≈1 — проскальзывание мало; можно использовать модель с III и дополнительными диссипативными силами. 2. Изолировать вклад момента инерции: сравнить шарики с одинаковой массо‑геометрией, но разным распределением массы (например, солидный шар и равного веса полый шар). Если аэродинамика и деформация минимальны (жёсткая дорожка, невысокие скорости), разница скоростей ≈ расчётной из I/(mr2)I/(mr^2)I/(mr2). 3. Оценить аэродинамику: провести серии в нормальной атмосфере и в вакуумной камере (или при пониженной плотности воздуха). Разность даёт вклад CdC_dCd. Или варьировать длину/угол так, чтобы скорости были сильно различны и подогнать терминологию к уравнению с v2v^2v2 — подгонка параметра CdC_dCd по профилю v(t)v(t)v(t). 4. Оценить роллинговое сопротивление: использовать очень жёсткую дорожку (минимум деформации) и «мягкую» дорожку; разница в ускорениях при прочих равных даёт вклад MrM_rMr или коэффициента CrrC_{rr}Crr. Можно измерить контактное пятно (оптически или давлением) и/или использовать датчик момента на оси (если шар ведётся по валу) для прямого измерения MrM_rMr. 5. Проскальзывание: при высоких углах/низком μs\mu_sμs измерять ω\omegaω и vvv — явное проскальзывание видно по v>rωv>r\omegav>rω (или наоборот). В этом режиме добавить кинетическое трение в модель. - Анализ данных: - Измерить v(t)v(t)v(t) из видео; численно подогнать дифференциальную модель mv˙=mgsinα−12ρCdAv2−Mrr−(в случае проскальзывания) Fk
m\dot v = m g\sin\alpha - \tfrac12\rho C_d A v^2 - \frac{M_r}{r} - \text{(в случае проскальзывания) }F_{k} mv˙=mgsinα−21ρCdAv2−rMr−(вслучаепроскальзывания) Fk
с учётом множителя 1+I/(mr2)1+I/(mr^2)1+I/(mr2) в эффективной инерции. Подгонкой по всем наборам углов извлечь CdC_dCd и MrM_rMr (или CrrC_{rr}Crr). - Сначала фиксируйте/измеряйте III и μs\mu_sμs отдельно, чтобы уменьшить число свободных параметров. 4) Какие погрешности надо учесть и как их уменьшить - Измерение угла α\alphaα: привести к минимальной систематике (лазерный уровень, цифровой нивелир). - Начальное условие: старт без импульса — централизованный держатель с электромагнитным отпуском. - Деление трения поверхности: неоднородности дорожки — полировать/замерить профиль. - Погрешность в r,m,Ir,m,Ir,m,I: измерять с калибровкой; III лучше проверять экспериментально. - Камера/тайминг: частота и разрешение, параллакс; калибровка расстояния по эталонному масштабу. - Внешние потоки воздуха, температура/плотность воздуха (влияет на CdC_dCd). - Статистические: повторять много запусков, оценивать стандартное отклонение. - Модельные допущения: постоянство CdC_dCd и MrM_rMr в зависимости от vvv; при больших скоростях/амплитудах эти параметры могут меняться. 5) Практическое упрощение и порядок приоритета - Сначала проверить чистое качение (видео и ω\omegaω vs vvv). Если соблюдается — основной вклад даст I/(mr2)I/(mr^2)I/(mr2) и роллинговое сопротивление; аэродинамика вторична для короткой низкой рампы. - Для количественной разборки: 1) измерить/зафиксировать III и μs\mu_sμs; 2) провести вакуум/атмосфера для отделения CdC_dCd; 3) менять жёсткость дорожки, чтобы оценить MrM_rMr. Итог: более быстрая стальная сфера, скорее всего, объясняется сочетанием меньшего относительного момента инерции (твердый шар против тонкой оболочки), меньших потерь на деформацию (жёсткость стали vs стекла и контакт с дорожкой) и меньшего влияния аэродинамики из‑за большей массы; чтобы разделить вклады количественно — измеряйте v(t)v(t)v(t) и ω(t)\omega(t)ω(t), сравните в вакууме/атмосфере и на жёсткой/мягкой дорожке и подгоняйте модель с параметрами CdC_dCd и MrM_rMr.
1) Физические причины (кратко)
- Момент инерции: для чистого качения без проскальзывания ускорение зависит только от отношения I/(mr2)I/(mr^2)I/(mr2):
a=gsinα1+I/(mr2). a=\frac{g\sin\alpha}{1+I/(mr^2)}.
a=1+I/(mr2)gsinα . Для твёрдого шара I=25mr2⇒a=57gsinαI=\tfrac{2}{5}mr^2\Rightarrow a=\tfrac{5}{7}g\sin\alphaI=52 mr2⇒a=75 gsinα, для тонкостенной оболочки I=23mr2⇒a=35gsinαI=\tfrac{2}{3}mr^2\Rightarrow a=\tfrac{3}{5}g\sin\alphaI=32 mr2⇒a=53 gsinα. Разница 5/73/5=2521≈1,19\dfrac{5/7}{3/5}=\dfrac{25}{21}\approx1{,}193/55/7 =2125 ≈1,19 — тонкостенная сфера теоретически ~19% медленнее при идеальном качении.
- Деформация / роллинговое сопротивление (hysteresis): упругие потери в материале и на контактном пятне создают дополнительный момент MrM_rMr (или эквивалентную тангенциальную силу). Это снижает ускорение и сильнее действует на более мягкие/комплайантные материалы (тонкое стекло/подложка).
- Трение качения / сцепление: чтобы было чистое качение, требуется достаточный статический коэффициент трения. При частичном проскальзывании кинетические потери замедляют шар сильнее.
- Аэродинамика: сопротивление воздуха Fdrag=12ρCdAv2F_{drag}=\tfrac12\rho C_d A v^2Fdrag =21 ρCd Av2 (или при малых Re — линейный закон). Для лёгкой пустой стеклянной оболочки аэродинамика влияет сильнее, т.к. тормозящее усилие не уменьшается при малой массе.
2) Простейшая модель с учётом основных потерь
Для поступательного ускорения с учётом аэродинамики и дополнительного тормозного момента MrM_rMr :
a=gsinα−12mρCdAv2−Mrmr1+I/(mr2). a=\frac{g\sin\alpha-\dfrac{1}{2m}\rho C_d A v^2-\dfrac{M_r}{mr}}{1+I/(mr^2)}.
a=1+I/(mr2)gsinα−2m1 ρCd Av2−mrMr . Условие отсутствия проскальзывания:
Ff=I/(mr2)1+I/(mr2) mgsinα≤μsmgcosα⇒μs≥I/(mr2)1+I/(mr2)tanα. F_f=\frac{I/(mr^2)}{1+I/(mr^2)}\,m g\sin\alpha\le\mu_s m g\cos\alpha
\quad\Rightarrow\quad
\mu_s\ge\frac{I/(mr^2)}{1+I/(mr^2)}\tan\alpha.
Ff =1+I/(mr2)I/(mr2) mgsinα≤μs mgcosα⇒μs ≥1+I/(mr2)I/(mr2) tanα.
3) Эксперимент — порядок действий для количественного разделения вкладов
- Подготовка:
- Использовать одинаковый радиус rrr для обоих шаров и жёсткую ровную направляющую платы (чтобы минимизировать вариации поверхности).
- Измерить массу mmm и геометрию; оценить/измерить III (по формуле по геометрии или экспериментально через крутильный маятник).
- Основные измерения:
- Снимать движение высокоскоростной камерой, отслеживать центр и угловую скорость ω(t)\omega(t)ω(t) (метка на шаре) → вычислять v(t)v(t)v(t) и отношение v/(rω)v/(r\omega)v/(rω) для обнаружения проскальзывания.
- Измерять время пробега ttt для набора углов α\alphaα и многократно повторять.
- Выделение вкладов:
1. Проверка чистого качения: если v/(rω)≈1v/(r\omega)\approx1v/(rω)≈1 — проскальзывание мало; можно использовать модель с III и дополнительными диссипативными силами.
2. Изолировать вклад момента инерции: сравнить шарики с одинаковой массо‑геометрией, но разным распределением массы (например, солидный шар и равного веса полый шар). Если аэродинамика и деформация минимальны (жёсткая дорожка, невысокие скорости), разница скоростей ≈ расчётной из I/(mr2)I/(mr^2)I/(mr2).
3. Оценить аэродинамику: провести серии в нормальной атмосфере и в вакуумной камере (или при пониженной плотности воздуха). Разность даёт вклад CdC_dCd . Или варьировать длину/угол так, чтобы скорости были сильно различны и подогнать терминологию к уравнению с v2v^2v2 — подгонка параметра CdC_dCd по профилю v(t)v(t)v(t).
4. Оценить роллинговое сопротивление: использовать очень жёсткую дорожку (минимум деформации) и «мягкую» дорожку; разница в ускорениях при прочих равных даёт вклад MrM_rMr или коэффициента CrrC_{rr}Crr . Можно измерить контактное пятно (оптически или давлением) и/или использовать датчик момента на оси (если шар ведётся по валу) для прямого измерения MrM_rMr .
5. Проскальзывание: при высоких углах/низком μs\mu_sμs измерять ω\omegaω и vvv — явное проскальзывание видно по v>rωv>r\omegav>rω (или наоборот). В этом режиме добавить кинетическое трение в модель.
- Анализ данных:
- Измерить v(t)v(t)v(t) из видео; численно подогнать дифференциальную модель
mv˙=mgsinα−12ρCdAv2−Mrr−(в случае проскальзывания) Fk m\dot v = m g\sin\alpha - \tfrac12\rho C_d A v^2 - \frac{M_r}{r} - \text{(в случае проскальзывания) }F_{k}
mv˙=mgsinα−21 ρCd Av2−rMr −(в случае проскальзывания) Fk с учётом множителя 1+I/(mr2)1+I/(mr^2)1+I/(mr2) в эффективной инерции. Подгонкой по всем наборам углов извлечь CdC_dCd и MrM_rMr (или CrrC_{rr}Crr ).
- Сначала фиксируйте/измеряйте III и μs\mu_sμs отдельно, чтобы уменьшить число свободных параметров.
4) Какие погрешности надо учесть и как их уменьшить
- Измерение угла α\alphaα: привести к минимальной систематике (лазерный уровень, цифровой нивелир).
- Начальное условие: старт без импульса — централизованный держатель с электромагнитным отпуском.
- Деление трения поверхности: неоднородности дорожки — полировать/замерить профиль.
- Погрешность в r,m,Ir,m,Ir,m,I: измерять с калибровкой; III лучше проверять экспериментально.
- Камера/тайминг: частота и разрешение, параллакс; калибровка расстояния по эталонному масштабу.
- Внешние потоки воздуха, температура/плотность воздуха (влияет на CdC_dCd ).
- Статистические: повторять много запусков, оценивать стандартное отклонение.
- Модельные допущения: постоянство CdC_dCd и MrM_rMr в зависимости от vvv; при больших скоростях/амплитудах эти параметры могут меняться.
5) Практическое упрощение и порядок приоритета
- Сначала проверить чистое качение (видео и ω\omegaω vs vvv). Если соблюдается — основной вклад даст I/(mr2)I/(mr^2)I/(mr2) и роллинговое сопротивление; аэродинамика вторична для короткой низкой рампы.
- Для количественной разборки: 1) измерить/зафиксировать III и μs\mu_sμs ; 2) провести вакуум/атмосфера для отделения CdC_dCd ; 3) менять жёсткость дорожки, чтобы оценить MrM_rMr .
Итог: более быстрая стальная сфера, скорее всего, объясняется сочетанием меньшего относительного момента инерции (твердый шар против тонкой оболочки), меньших потерь на деформацию (жёсткость стали vs стекла и контакт с дорожкой) и меньшего влияния аэродинамики из‑за большей массы; чтобы разделить вклады количественно — измеряйте v(t)v(t)v(t) и ω(t)\omega(t)ω(t), сравните в вакууме/атмосфере и на жёсткой/мягкой дорожке и подгоняйте модель с параметрами CdC_dCd и MrM_rMr .