Соколов Алексей. "Динамика несферических частиц" — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Задание характеристик частиц)
 
(не показано 6 промежуточных версий 3 участников)
Строка 9: Строка 9:
 
Идея метода состоит в том, чтобы переходить в систему отсчета одной из частиц, и проверять, находятся ли углы внутри частицы
 
Идея метода состоит в том, чтобы переходить в систему отсчета одной из частиц, и проверять, находятся ли углы внутри частицы
  
[[Файл:2323.png |300 px|Детектирование столкновений]]
+
[[Файл:Detec.PNG |300 px|Детектирование столкновений]]
  
 
Т.о. если выполняется условие
 
Т.о. если выполняется условие
Строка 23: Строка 23:
 
<math>\left\{   
 
<math>\left\{   
 
           \begin{array}{rcl}   
 
           \begin{array}{rcl}   
             m \frac{d^2r}{dt^2} & = & \sum^{n}_{j=1} {F}_{\j} \\  
+
             m \frac{d^2 \vec{r}}{dt^2} & = & \sum^{n}_{j=1} {\vec{F}}_{\j} \\  
 
               &\\
 
               &\\
             J \frac{d^2w}{dt^2} & = & \sum^{m}_{i=1} {M}_{\i} \\   
+
             J \frac{d^2\vec{w}}{dt^2} & = & \sum^{m}_{i=1} {\vec{M}}_{\i} \\   
 
              
 
              
 
           \end{array}   
 
           \end{array}   
Строка 32: Строка 32:
 
==== Leapfrog интегрирование ====  
 
==== Leapfrog интегрирование ====  
  
[[Файл:leapfrog.png |180 px|Leapfrog интегрирование]]
+
<math>\left\{ 
 +
          \begin{array}{rcl} 
 +
            \vec{v}_{i+1} & = & \vec{v}_i + \vec{a}_i\Delta t \\
 +
              &\\
 +
            \vec{r}_{i+1} & = & \vec{r}_i + \vec{v}_{i+1}\Delta t \\
 +
           
 +
          \end{array} 
 +
          \right. </math>
  
 
==Диссипативная модель==
 
==Диссипативная модель==
Строка 39: Строка 46:
  
 
<br>, где <math> k_{1}, k_{2}</math> - коэффициенты упругости  <br> <math> \beta </math> - коэффициент вязкого трения  <br> <math> \mu </math> - коэффициент сухого трения
 
<br>, где <math> k_{1}, k_{2}</math> - коэффициенты упругости  <br> <math> \beta </math> - коэффициент вязкого трения  <br> <math> \mu </math> - коэффициент сухого трения
 +
 +
== См. также ==
 +
 +
*[[Соколов Алексей]]
 +
  
 
[[Category: Студенческие проекты]]
 
[[Category: Студенческие проекты]]
 +
[[Category: Механика дискретных сред]]
 +
 +
[[Category: Личные страницы]]

Текущая версия на 17:39, 29 мая 2014

Модель взаимодействия квадратных частиц в 2D

Вектора и углы

Задание характеристик частиц[править]

Каждая частица имеет радиус вектор и пару ортогональных векторов. Таким образом определяем положение углов.

Детектирование столкновений[править]

Идея метода состоит в том, чтобы переходить в систему отсчета одной из частиц, и проверять, находятся ли углы внутри частицы

Детектирование столкновений

Т.о. если выполняется условие

Вектора и углы

то частицы находятся в контакте.

Упругие силы и моменты[править]

[math]\vec{F} = kl\vec{n}[/math], где [math]\vec{n}[/math] - нормаль к поверхности

Динамические уравнения[править]

[math]\left\{ \begin{array}{rcl} m \frac{d^2 \vec{r}}{dt^2} & = & \sum^{n}_{j=1} {\vec{F}}_{\j} \\ &\\ J \frac{d^2\vec{w}}{dt^2} & = & \sum^{m}_{i=1} {\vec{M}}_{\i} \\ \end{array} \right. [/math]

Leapfrog интегрирование[править]

[math]\left\{ \begin{array}{rcl} \vec{v}_{i+1} & = & \vec{v}_i + \vec{a}_i\Delta t \\ &\\ \vec{r}_{i+1} & = & \vec{r}_i + \vec{v}_{i+1}\Delta t \\ \end{array} \right. [/math]

Диссипативная модель[править]

диссипативная модель


, где [math] k_{1}, k_{2}[/math] - коэффициенты упругости
[math] \beta [/math] - коэффициент вязкого трения
[math] \mu [/math] - коэффициент сухого трения

См. также[править]