"Численные методы интегрирования уравнений движения для одномерной линейной цепочки и частицы в потенциальной яме Леннарда-Джонса" — различия между версиями
Материал из Department of Theoretical and Applied Mechanics
(Новая страница: «'''''Курсовой проект по Механике дискретных сред''''' '''Исполни…») |
Totamonik (обсуждение | вклад) |
||
| Строка 8: | Строка 8: | ||
==Постановка задачи:== | ==Постановка задачи:== | ||
| − | # Сравнить различные методы интегрирования уравнений движения (Верле, Рунге-Кутта 4 порядка). Реализовать фиксированные, свободные и периодические условия | + | # Сравнить различные методы интегрирования уравнений движения (Верле, Рунге-Кутта 4 порядка). Реализовать фиксированные, свободные и периодические условия |
# Численно определить скорость диссоциации частицы в потенциальной яме Леннарда-Джонса | # Численно определить скорость диссоциации частицы в потенциальной яме Леннарда-Джонса | ||
==Теоретическая сводка:== | ==Теоретическая сводка:== | ||
| − | + | ==== 1. Одномерная линейная цепочка ==== | |
Рассмотрим модель колебаний одинаковых атомов массой m, находящихся в одномерной цепочке. Пусть в этой цепочке находится N атомов, связанных между собой квазиупругой силой с коэффициентом упругости k. | Рассмотрим модель колебаний одинаковых атомов массой m, находящихся в одномерной цепочке. Пусть в этой цепочке находится N атомов, связанных между собой квазиупругой силой с коэффициентом упругости k. | ||
| Строка 18: | Строка 18: | ||
[[Image:Изображение1.png|top]] | [[Image:Изображение1.png|top]] | ||
| − | + | [[Image:Изображение22331.png|top]] | |
| − | + | Если учитывать взаимодействие только соседних атомов, уравнение движения можно записать в следующем виде: | |
| − | + | [[Image:Изображение22333.png|top]] | |
| − | + | Для решения уравнения движения воспользуемся численными методами интегрирования: | |
| − | + | а. Метод Верле | |
| − | |||
б. Метод Рунге-Кутта 4 порядка | б. Метод Рунге-Кутта 4 порядка | ||
| − | + | [[Image:Изображение22334.png|top]] | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| + | [[Image:Изображение22335.png|top]] | ||
| + | [[Image:Изображение22336.png|top]] | ||
| − | + | [[Image:Изображение22337.png|top]] | |
| − | + | [[Image:Изображение22338.png|top]] | |
| − | + | [[Image:Изображение22339.png|top]] | |
| − | + | [[Image:Изображение223310.png|top]] | |
| − | + | [[Image:Изображение223311.png|top]] | |
| + | [[Image:Изображение223312.png|top]] | ||
| − | + | [[Image:Изображение223313.png|top]] | |
| + | [[Image:Изображение223314.png|top]] | ||
| + | [[Image:Изображение223315.png|top]] | ||
| − | |||
| − | + | Для каждого из методов реализуются 3 вида граничных условий: | |
| − | + | * Фиксированные граничные условия | |
| + | [[Image:Изображение223316.png|top]] | ||
| − | + | [[Image:Изображение223317.png|top]] | |
| + | [[Image:Изображение223318.png|top]] | ||
| + | [[Image:Изображение223319.png|top]] | ||
| − | + | * Свободные граничные условия | |
| − | + | [[Image:Изображение223320.png|top]] | |
| − | + | [[Image:Изображение223321.png|top]] | |
| + | * Периодические граничные условия | ||
| − | + | [[Image:Изображение223322.png|top]] | |
| + | [[Image:Изображение223323.png|top]] | ||
| + | |||
| + | ====2. Частица в потенциальной яме Леннарда-Джонса ==== | ||
Уравнение движения частицы в потенциальной яме Леннарда-Джонса: | Уравнение движения частицы в потенциальной яме Леннарда-Джонса: | ||
| − | + | [[Image:Изображение223324.png|top]] | |
| − | + | [[Image:Изображение223325.png|top]] | |
| − | [[Image: | + | [[Image:Изображение223326.png|top]] |
Скоростью диссоциации будем называть скорость, которую необходимо сообщить частице, чтобы она улетела на бесконечность. | Скоростью диссоциации будем называть скорость, которую необходимо сообщить частице, чтобы она улетела на бесконечность. | ||
| − | [[Image: | + | [[Image:Изображение223327.png|top]] |
| − | Решение: | + | ==Решение:== |
| − | Вывод: | + | ==Вывод:== |
| − | |||
Версия 22:58, 23 января 2020
Курсовой проект по Механике дискретных сред
Исполнитель: Троцкая Валерия
Группа: 3630103/60101
Семестр: осенний семестр 2019-2020 учебного года
Содержание
Постановка задачи:
- Сравнить различные методы интегрирования уравнений движения (Верле, Рунге-Кутта 4 порядка). Реализовать фиксированные, свободные и периодические условия
- Численно определить скорость диссоциации частицы в потенциальной яме Леннарда-Джонса
Теоретическая сводка:
1. Одномерная линейная цепочка
Рассмотрим модель колебаний одинаковых атомов массой m, находящихся в одномерной цепочке. Пусть в этой цепочке находится N атомов, связанных между собой квазиупругой силой с коэффициентом упругости k.
Если учитывать взаимодействие только соседних атомов, уравнение движения можно записать в следующем виде:
Для решения уравнения движения воспользуемся численными методами интегрирования:
а. Метод Верле
б. Метод Рунге-Кутта 4 порядка
Для каждого из методов реализуются 3 вида граничных условий:
- Фиксированные граничные условия
- Свободные граничные условия
- Периодические граничные условия
2. Частица в потенциальной яме Леннарда-Джонса
Уравнение движения частицы в потенциальной яме Леннарда-Джонса:
Скоростью диссоциации будем называть скорость, которую необходимо сообщить частице, чтобы она улетела на бесконечность.
