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

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(уеапре)
(Исходный код программы)
 
(не показано 7 промежуточных версий 2 участников)
Строка 1: Строка 1:
'''''Курсовой проект по [[Механика дискретных сред|Механике дискретных сред]]'''''
+
'''''Курсовой проект по [[Введение в механику дискретных сред|Введению в механику дискретных сред]]'''''
  
'''Исполнитель:''' [[Михеев Евгений]]
+
'''Исполнитель:''' [[Стешова Юлия]]
  
 
'''Группа:'''  5030103/90101
 
'''Группа:'''  5030103/90101
Строка 8: Строка 8:
  
 
==Постановка задачи==
 
==Постановка задачи==
Требутся смоделировать статическое провисание и падение цепочки, состоящей из частиц, соединенных упругими связями под воздействем силы тяжести, а также исследовать зависимость ускорения крайней свободной частицы от времени.
+
Цепочка состоит из частиц, имеющих линейные упругие связи, и находится в статическом равновесии. На первом шаге две крайние частицы не закреплены. Требутся смоделировать падение под действем силы тяжести одного из концов цепочки и исследовать зависимость ускорения крайней свободной частицы от времени и сравнить с ускорением свободно падающей частицы.
  
 
==Математическая модель==
 
==Математическая модель==
Строка 33: Строка 33:
 
</math>,  где <math>c</math> - коэффициент жесткости пружины.
 
</math>,  где <math>c</math> - коэффициент жесткости пружины.
  
Поделив уравнение на mg можно произвести обезразмеривание. Обезразмеренное уравнение будет иметь вид:
+
Обезразмеренное уравнение будет иметь вид:
  
 
<math>
 
<math>
Строка 43: Строка 43:
  
 
==Исходный код программы==
 
==Исходный код программы==
Исходный код программы на языке программирования С# можно найти на GitHub:
+
Исходный код программы на языке программирования Matlab, находится на GitHub:
 
+
https://github.com/jstesh/new_rep
https://github.com/EMcHave/Oscillator
 
 
 
<div align = "center">{{#widget:Iframe |url =  http://tm.spbstu.ru/htmlets/Mikheev/chain.png|width=800|height=800}}</div>
 

Текущая версия на 12:22, 25 января 2023

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

Исполнитель: Стешова Юлия

Группа: 5030103/90101

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

Постановка задачи[править]

Цепочка состоит из частиц, имеющих линейные упругие связи, и находится в статическом равновесии. На первом шаге две крайние частицы не закреплены. Требутся смоделировать падение под действем силы тяжести одного из концов цепочки и исследовать зависимость ускорения крайней свободной частицы от времени и сравнить с ускорением свободно падающей частицы.

Математическая модель[править]

[math] m\underline{\ddot{r}}_i(t)=\underline{F}_{R_1}+\underline{F}_{R_2} + \underline{F}_{g}\\ \underline{r}_i(0)=\underline{r}_i^0,~\underline{v}_i(0)=0~~~i=1,\ldots,n [/math]


где [math] \underline{F}_{R_1}, \underline{F}_{R_2}\\ [/math] - силы упругости действующие на [math]i[/math]-ую частицу со стороны [math]i-1[/math] и [math]i+1[/math] соответственно;

[math] \underline{F}_{g} = -mg\underline{k} \\ [/math] - сила тяжести;

Сила упругости, возникающая в пружине соединяющей частицу 1 и 2, вычисляется по следующей формуле:

[math] \underline{F}_{R}= c(||\underline{r}_2-\underline{r}_1|| - l_0)\frac{(\underline{r}_2-\underline{r}_1)}{||\underline{r}_2-\underline{r}_1||} [/math], где [math]c[/math] - коэффициент жесткости пружины.

Обезразмеренное уравнение будет иметь вид:

[math] \underline{\ddot{r}}_i(t)= \frac{cl_0}{mg}(||\underline{r}_{i+1}-\underline{r}_i|| - 1)\frac{(\underline{r}_{i+1}-\underline{r}_i)}{||\underline{r}_{i+1}-\underline{r}_i||} + \frac{cl_0}{mg}(||\underline{r}_{i-1}-\underline{r}_i|| - 1)\frac{(\underline{r}_{i-1}-\underline{r}_i)}{||\underline{r}_{i-1}-\underline{r}_i||} - \underline{k}\\ [/math]

Интегрирование по времени производится неявной схемой интегрирования.

Исходный код программы[править]

Исходный код программы на языке программирования Matlab, находится на GitHub: https://github.com/jstesh/new_rep