Вынужденные колебания цепочки в вязкой среде — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Теоретическая сводка)
(Граничные и начальные условия)
Строка 33: Строка 33:
 
:Периодические граничные условия:
 
:Периодические граничные условия:
 
:<math>V_{N-1}=c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{N-1}+Qcos(ωt)</math>
 
:<math>V_{N-1}=c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{N-1}+Qcos(ωt)</math>
:<math>V_{0}=-c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{0}
+
:<math>V_{0}=-c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{0}</math>
  
 
==Решение==
 
==Решение==
Строка 47: Строка 47:
  
 
[[File:Пуля3.gif|500px|]]
 
[[File:Пуля3.gif|500px|]]
 
 
  
 
=== См. также ===
 
=== См. также ===

Версия 19:26, 22 января 2020

Курсовой проект по Механике дискретных сред

Исполнитель: Шпетный Даниил

Группа: 3630103/60101

Семестр: осень 2019


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

Исследовать вынужденные колебания в вязкой среде. Построить график перемещений частиц от времени.

Построение модели

Процесс моделируется как одномерные колебания цепочки частиц. Уравнение взаимодействия : (строчка с F)

Теоретическая сводка

Парное взаимодействие определяется формулой:

[math]F = c(x_{k+1}-x_{k}-1)[/math]

где [math]F[/math] — сила взаимодействия, [math]c[/math] — жесткость связи, [math]x_k[/math] — перемещение частицы, [math]k[/math] — номер частицы.

Уравнение скорости частиц [math]k[/math] и [math]k+1[/math]:

[math]V_{k}=\frac{F}{m}dt-2μV_{k}[/math]
[math]V_{k+1}=-\frac{F}{m}dt-2μV_{k+1}[/math]

где [math]F[/math] — сила взаимодействия, [math]m[/math] — масса частицы, [math]μ[/math] — коэффициент вязкости, [math]V_{k}[/math] — скорость частицы.

Граничные и начальные условия

Начальные условия нулевые

Периодические граничные условия:
[math]V_{N-1}=c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{N-1}+Qcos(ωt)[/math]
[math]V_{0}=-c(x_{0}+aN-x_{N-1}-1)\frac{dt}{m}-2μV_{0}[/math]

Решение

При построении модели были приняты следующие значения параметров:

[math]a=1[/math]

[math]D=1[/math]

[math]m=1[/math]

[math]dt=0.001[/math] -шаг по времени

Пуля3.gif

См. также