КП: Диск Эйлера — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Общие сведения по теме)
(Решение)
Строка 83: Строка 83:
 
<br> <math> \b{a} \times \b{N} = \frac{a \eta}{ e} N \b{k} \times \b{n} </math>
 
<br> <math> \b{a} \times \b{N} = \frac{a \eta}{ e} N \b{k} \times \b{n} </math>
 
<br> В результате получаем систему, описывающую движение тела:
 
<br> В результате получаем систему, описывающую движение тела:
<br>
+
<br> <math>\left\{   
<math>\left\{   
 
 
           \begin{array}{lcl}   
 
           \begin{array}{lcl}   
 
             \Theta_{12} \b{n} \times \ddot{\b{n}} + \Theta_3 \dot{ \left(\Omega \b{n} \right) } + a \frac{\eta}{\epsilon} N \b{n} \times \b{k} = 0 \\  
 
             \Theta_{12} \b{n} \times \ddot{\b{n}} + \Theta_3 \dot{ \left(\Omega \b{n} \right) } + a \frac{\eta}{\epsilon} N \b{n} \times \b{k} = 0 \\  
Строка 93: Строка 92:
 
           \end{array}   
 
           \end{array}   
 
           \right.  </math>
 
           \right.  </math>
 +
 +
<br>Будем искать решение для стационарного движения:
 +
<br> <math>\left\{ 
 +
          \begin{array}{lcl} 
 +
          \underline{\omega} = \dot{\psi} \b{k} + \dot{\varphi} \b{k} \\
 +
          \underline{\omega} \times \b{a} = 0
 +
\end{array}
 +
\right. </math>
 +
<br> Тогда:
 +
<br> <math> \dot{ \b{n}} = \dot{\psi} \b{k} \times \b{n} </math>
 +
<br> <math> \ddot{ \b{n}} =  \underline{\omega} \times \dot{ \b{n}} =
 +
\left( \dot{\psi} \b{k} + \dot{\varphi} \b{k} \right) \times  \dot{\psi} \b{k} \times \b{n} = (\dot{\psi})^2 \left( \b{k} \cos(\vartheta) - \b{n} \right) + \dot{\psi} \dot{\varphi} \left( \b{k} - \b{n} \cos(\vartheta) \right) </math>
 +
<br> <math> \b{n} \times \ddot{ \b{n}} = \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) \b{n} \times \b{k} </math>
 +
<br> <math> \Omega = \underline{ \omega } \cdot \b{n} = \dot{ \psi } \cos( \vartheta ) + \dot{ \phi } </math>
 +
<br> <math> \dot{\left( \Omega \b{n} \right)} = \Omega \dot{ \b{n}} = \Omega \dot{ \psi } \b{k} \times \b{n} </math>
 +
<br> Подставляем полученные выражения в уравнение (1)
 +
<br> <math> \left[ \left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N \right] \b{n} \times \b {k} = 0 </math>
 +
<br> Получаем уравнение:
 +
<br> <math>\left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N = 0 </math>
 +
<br> Из условия неподвижности точки касания <math> \underline{\omega} \times \b{a} = 0 </math> получаем:
 +
<br> <math> \underline{\omega} \times \b{a} = 0 </math>
 +
<br> <math> \left( \dot{\psi} \b{k} + \dot{\varphi} \b{k}\right) \times \b{a} = 0 </math>
 +
<br> <math> \dot{\psi} \b{k} \times \b{a} + \dot{\varphi} \b{n} \times \b{a} = 0 </math>
 +
<br> <math> a \left( \dot{ \psi } \cot( \vartheta) + \frac{ \dot{ \varphi } }  \sin( \vartheta) \right) = 0 </math>
 +
<br> Получаем систему:
 +
<br> <math>\left\{ 
 +
          \begin{array}{lcl} 
 +
          \left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N = 0 \\
 +
            m \ddot{z} = N - mg; z = - a \sin( \vartheta ) \\
 +
          \dot{ \psi } \sin( \vartheta ) + \dot{ \varphi } = 0
 +
\end{array}
 +
\right. </math>
 +
<br> Имеем: <math> \ddot{z} = 0 </math>
 +
<br> Таким образом <math> N = mg </math>
 +
<br> Имеем
 +
<br> <math> (\dot{\psi})^2 \left( \Theta_{12} - \Theta_3 \right) \left( \cos( \vartheta) - \sin( \vartheta ) \right) + a \cot( \vartheta ) mg = 0 </math>
 +
<br> Получаем следующую зависимость:
 +
<br> <math> \dot{ \psi } = \sqrt{ \frac{a m g \cos( \vartheta)}{ \left( \Theta_{12} - \Theta_3 \right) \left( \sin^2( \vartheta) - \sin( \vartheta) \cos(vartheta) \right) } } </math>
  
 
== Обсуждение результатов и выводы ==
 
== Обсуждение результатов и выводы ==

Версия 19:16, 3 июня 2013

А.М. Кривцов > Теоретическая механика > Курсовые проекты 2013 > Диск Эйлера
Eulerdisc.png

Курсовой проект по Теоретической механике

Исполнитель: Соколов Алексей

Группа: 07 (20510)

Семестр: весна 2013

Аннотация проекта

Диск Эйлера – это твёрдое цилиндрическое тело, имеющее единственную точку контакта с горизонтальной поверхностью, по которой он одновременно и катится, и вращается. Этой простой механической системе присущи две любопытные особенности — быстрое возрастание звуковой частоты в процессе контакта диска и поверхности на финальном этапе движения, а также последующая внезапная остановка движения. Эти эффекты может наблюдать каждый, раскрутив на столе обыкновенную монету: «потеряв равновесие», монета окажется лежащей на столе плашмя, при этом финальная стадия ее движения, когда ее плоскость почти горизонтальна, сопровождается характерным «дрожанием». Однако физическое объяснение этого своеобразного поведения оказалось сложной задачей, решение которой еще требует дальнейших экспериментальных и теоретических исследований.

Постановка задачи

Euler.jpg
Рассмотрим стационарное движение твердого диска с одной точкой касания без диссипаций. Для описания тела будем использовать единичный вектор нормали [math] \b{n} [/math] и вектор [math] \b{a} [/math] , находящегося в плоскости диска с длиной, равной радиусу диска.

Общие сведения по теме

Повышенный интерес к диску Эйлера возник после опубликования работы К. Моффатта, в которой резкое прекращение движения диска объясняется прежде всего воздействием силы вязкого сопротивления воздуха. Эта теория вызвала значительную дискуссию и развитие конкурентных гипотез о природе основного механизма диссипации энергии для диска Эйлера. Ван ден Энг (van den Engh) с соавторами провели некоторые эксперименты по движению дископодобных тел, свидетельствующие не в пользу гипотезы о вязком трении. В ответ на работу Моффатта они опубликовали свое объяснение, ключевым пунктом которого является наличие скольжения при контакте диска и поверхности. Критика выводов и краткое обсуждение возможных источников диссипации содержатся также в неопубликованной заметке Э. Руины. Авторами последовавших экспериментальных и теоретических работ исследовался вопрос о наличии и степени влияния определенных типов трения на различных этапах движения диска. Были численно и аналитически исследованы модели системы с различными ограничениями, проведен анализ полученных экспериментальных данных. Эти исследования в основном указывают на то, что основными диссипативными силами, вызывающими вибрацию и остановку диска, являются силы трения качения и трения скольжения, нежели силы вязкого трения. Из немногих предположений о физике этих эффектов стоит упомянуть, как наиболее естественную, гипотезу Кесслера и О’Рейли, что резкая остановка диска происходит в результате потери контакта между диском и поверхностью в процессе вибраций при малом угле наклона диска. (Для проверки этой гипотезы потребуется рассматривать деформируемую модель для контактирующих тел.)

Решение

Запишем кинетический момент сиcтемы:
[math] \dot{\b{K}} = \dot{ \left( \underline{\underline{\Theta}} \cdot \underline{\omega} \right) } = \b{M} [/math]
[math] \dot{\b{K}} = \Theta_{12} \b{n} \times \ddot{\b{n}} + \Theta_3 \dot{ \left(\Omega \b{n} \right) } [/math]
Положим
[math] \eta = \b{k} \cdot \b{n} = \cos(\vartheta) [/math]
[math] \epsilon = \b{k} \times \b{n} = \sin(\vartheta) [/math]
Имеем далее
[math]\left\{ \begin{array}{lcl} \b{a} \perp \b{k} \times \b{n} \\ \b{a} \perp \b{n} \\ \end{array} \right. [/math]

Из этого следует:
[math] \b{a} \parallel \left[ \b{n} \times ( \b{k} \times \b{n}) \right ] [/math]
Имеем
[math] \b{n} \times ( \b{k} \times \b{n} ) = \b{k} ( \b{n} \cdot \b{n}) - \b{n} (\b{n} \cdot \b{k}) = \b{k} - \eta \b{n} [/math]
Таким образом можем представить вектор [math] \b{a} [/math] следующим образом:
[math] \b{a} = \lambda( \b{k} - \eta \b{n}) [/math]
Имеем далее:
[math] a^2 = \lambda^2(1 - \eta^2) = \lambda^2 \epsilon [/math]
Таким образом
[math] \lambda = \frac{a}{ \sqrt{1 - \eta^2} } [/math]
Недтрудно понять что вектор [math] \b{a} [/math] в таком случае записывается в следующем виде:
[math] \b{a} = \frac{a}{ \sqrt{1 - \eta^2} } ( \b{k} - \eta \b{n}) [/math]
Имеем далее:
[math] \frac{ \b{a} } {a} = \frac{ \b{k} - \eta \b{n} } { \left| \b{k} - \eta \b{n} \right| } [/math]
Получаем таким образом:
[math] \b{a} = \frac{a}{\epsilon} (\b{k} - \eta \b{n}) [/math]
Теперь запишем соотношения для сил и моментов:
[math] m \ddot{\b{r}} = m \b{g} + \b{N} = (N - mg) \b{k} [/math]
[math] \dot{ \left( \b{\b{O}} \cdot \b{w} \right) } = \b{a} \times \b{N} [/math]
[math] \b{N} = N \b{k} [/math]
Таким образом имеем:
[math] \b{a} \times \b{N} = \frac{a \eta}{ e} N \b{k} \times \b{n} [/math]
В результате получаем систему, описывающую движение тела:
[math]\left\{ \begin{array}{lcl} \Theta_{12} \b{n} \times \ddot{\b{n}} + \Theta_3 \dot{ \left(\Omega \b{n} \right) } + a \frac{\eta}{\epsilon} N \b{n} \times \b{k} = 0 \\ m \ddot{z} = N - mg \\ z = - a \epsilon \\ \end{array} \right. [/math]


Будем искать решение для стационарного движения:
[math]\left\{ \begin{array}{lcl} \underline{\omega} = \dot{\psi} \b{k} + \dot{\varphi} \b{k} \\ \underline{\omega} \times \b{a} = 0 \end{array} \right. [/math]
Тогда:
[math] \dot{ \b{n}} = \dot{\psi} \b{k} \times \b{n} [/math]
[math] \ddot{ \b{n}} = \underline{\omega} \times \dot{ \b{n}} = \left( \dot{\psi} \b{k} + \dot{\varphi} \b{k} \right) \times \dot{\psi} \b{k} \times \b{n} = (\dot{\psi})^2 \left( \b{k} \cos(\vartheta) - \b{n} \right) + \dot{\psi} \dot{\varphi} \left( \b{k} - \b{n} \cos(\vartheta) \right) [/math]
[math] \b{n} \times \ddot{ \b{n}} = \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) \b{n} \times \b{k} [/math]
[math] \Omega = \underline{ \omega } \cdot \b{n} = \dot{ \psi } \cos( \vartheta ) + \dot{ \phi } [/math]
[math] \dot{\left( \Omega \b{n} \right)} = \Omega \dot{ \b{n}} = \Omega \dot{ \psi } \b{k} \times \b{n} [/math]
Подставляем полученные выражения в уравнение (1)
[math] \left[ \left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N \right] \b{n} \times \b {k} = 0 [/math]
Получаем уравнение:
[math]\left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N = 0 [/math]
Из условия неподвижности точки касания [math] \underline{\omega} \times \b{a} = 0 [/math] получаем:
[math] \underline{\omega} \times \b{a} = 0 [/math]
[math] \left( \dot{\psi} \b{k} + \dot{\varphi} \b{k}\right) \times \b{a} = 0 [/math]
[math] \dot{\psi} \b{k} \times \b{a} + \dot{\varphi} \b{n} \times \b{a} = 0 [/math]
[math] a \left( \dot{ \psi } \cot( \vartheta) + \frac{ \dot{ \varphi } } \sin( \vartheta) \right) = 0 [/math]
Получаем систему:
[math]\left\{ \begin{array}{lcl} \left( \Theta_{12} - \Theta_3 \right) \left( (\dot{\psi})^2 \cos( \vartheta ) + \dot{\psi} \dot{\varphi} \right) + a \cot( \vartheta )N = 0 \\ m \ddot{z} = N - mg; z = - a \sin( \vartheta ) \\ \dot{ \psi } \sin( \vartheta ) + \dot{ \varphi } = 0 \end{array} \right. [/math]
Имеем: [math] \ddot{z} = 0 [/math]
Таким образом [math] N = mg [/math]
Имеем
[math] (\dot{\psi})^2 \left( \Theta_{12} - \Theta_3 \right) \left( \cos( \vartheta) - \sin( \vartheta ) \right) + a \cot( \vartheta ) mg = 0 [/math]
Получаем следующую зависимость:
[math] \dot{ \psi } = \sqrt{ \frac{a m g \cos( \vartheta)}{ \left( \Theta_{12} - \Theta_3 \right) \left( \sin^2( \vartheta) - \sin( \vartheta) \cos(vartheta) \right) } } [/math]

Обсуждение результатов и выводы

В результате работы найден простой метод описания данной динамической системы. Получены уравнения, описывающие движение диска.

Ссылки по теме

См. также