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

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

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

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

Группа: 3630103/60101

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

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

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

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

Процесс моделируется как одномерные колебания цепочки частиц.

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

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

[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}=c(x_{1}+aN-x_{N}-1)\frac{dt}{m}-2μV_{N}+Qcos(ωt)[/math]

[math]V_{1}=-c(x_{1}+aN-x_{N}-1)\frac{dt}{m}-2μV_{1}[/math]

где [math]Q[/math] — амплитуда возмущающей силы, [math]ω[/math] — частота возмущающей силы, [math]N[/math] — количество частиц.

Решение

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

[math]a=1[/math]

[math]D=1[/math]

[math]m=1[/math]

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

См. также

• Метод динамики частиц
• Механика дискретных сред
• Введение в механику дискретных сред
• Виртуальная лаборатория
• Курсовые работы 2012-2013 учебного года
• Курсовые работы 2013-2014 учебного года
• Курсовые работы 2014-2015 учебного года
• Курсовые работы 2016-2017 учебного года