Редактирование: Колебания одномерной цепочки
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 1: | Строка 1: | ||
− | + | '''''Курсовой проект по Механике дискретных сред''''' | |
− | ''''' | + | '''Исполнитель:''' [[Лобанов Илья]] |
− | ''' | + | '''Группа:''' 43604/1 |
− | |||
− | + | == Постановка задачи== | |
+ | Рассматривается цепочка из N материальных частиц P0, P1, ... Pn, ..., PN–1, каждая из которых обладает одинаковой массой m. Частицы соединены линейными пружинками. Рассматриваются продольные колебания образующих цепочку частиц под действием сил взаимодействия между частицами цепочки, а также параллельных направлению цепочки внешних сил. | ||
+ | Уравнение движения имеет вид: | ||
+ | ::<math> | ||
+ | {m}\ddot{\bf U}_{i} = {С}({\bf U}_{i-1}-2{\bf U}_{i} + {\bf U}_{i+1}), | ||
+ | </math> | ||
− | + | где С - жёсткость одной пружинки, m - масса одной частицы, <math> {\bf U}_{i} </math> - перемещение частицы, a - расстояние между двумя соседними частицами в начальный момент времени. | |
− | |||
Период одного колебания: <math> {T}_{o} = 2{\pi}\sqrt\frac {m}{C} </math> | Период одного колебания: <math> {T}_{o} = 2{\pi}\sqrt\frac {m}{C} </math> | ||
− | + | Движение частицы с номером n описывается зависимостью от времени t её смещения Un относительно положения равновесия этой частицы (узла цепочки с номером n).Взаимодействие задаётся с помощью потенциала Леннарда-Джонса | |
− | Рассмотрим модель колебаний одномерной многоатомной цепочки равных масс. Пусть в этой цепочке находится N атомов. Обозначим смещение n-го атома | + | |
+ | ==Решение задачи== | ||
+ | Рассмотрим модель колебаний одномерной многоатомной цепочки равных масс. Пусть в этой цепочке находится N атомов. Обозначим смещение n-го атома un, а атома, отстоящего от него на p узлов, – un+p. Примем в качестве положительных смещения атомов вправо от положения равновесия, а отрицательных – влево. | ||
− | Каждый атом смещается только вдоль цепочки, что следует из требования одномерности модели. Пусть атомы связаны между собой упругой силой | + | Каждый атом смещается только вдоль цепочки, что следует из требования одномерности модели. Пусть атомы связаны между собой упругой силой F с коэффициентом упругости с. |
::<math> | ::<math> | ||
F(r) = \frac{12D}{a}\left[-\left(\frac{a}{r}\right)^{13} + \left(\frac{a}{r}\right)^{7}\right]. | F(r) = \frac{12D}{a}\left[-\left(\frac{a}{r}\right)^{13} + \left(\frac{a}{r}\right)^{7}\right]. | ||
</math> | </math> | ||
− | Найдем уравнение движения n-го и n+1-го атома в цепи. В равновесном положении силы, действующие на атомы, равны нулю. При смещениях на каждый n-й атом будет действовать сила со стороны соседних атомов. | + | Найдем уравнение движения n-го и n+1-го атома в цепи. В равновесном положении силы, действующие на атомы, равны нулю. При произвольных смещениях на каждый n-й атом будет действовать сила со стороны соседних атомов. В соответствии с элементарным законом Гука эту силу можно представить в виде: |
+ | |||
+ | В качестве начальных условий заданы случайные начальные скорости таким образом, что средняя скорость всех частиц равна 0. Перемещения всех частиц в начальный момент времени равны нулю. Также заданы периодические граничные условия на перемещения. | ||
− | + | Для решения данного дифференциального уравнения использовали метод Верле | |
− | |||
− | |||
− | + | [[File:Energy12.jpg|thumb]] | |
− | == | + | ==Анализ результатов== |
− | + | '''График продольных колебаний''' | |
− | + | [[File:rWqGLWc-jcg.jpg]] | |
− | + | [[File:energy12.jpg|thumb|График изменения энергии в системе]] | |
− | |||
− |