Перераспределение энергии между степенями свободы в нелинейной двухатомной цепочке — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Решение)
Строка 21: Строка 21:
 
В задаче на каждом шаге по времени находилась кинетическая энергия всей системы, после чего был построен график зависимости энергии системы от времени.  
 
В задаче на каждом шаге по времени находилась кинетическая энергия всей системы, после чего был построен график зависимости энергии системы от времени.  
 
Для частиц одинаковой массы был получен следующий график:
 
Для частиц одинаковой массы был получен следующий график:
 +
 +
Если рассматривать частицы разной массы, то график зависимости будет иметь следующий вид (m1=1, m2=1.3):
 +
 +
Рассмотрим нелинейную постановку задачи. В выражении для силы будет присутствовать слагаемое третьего порядка:
 +
[[File:Формула8.jpg|center]]
 +
Решается данная задача в нелинейной постановке аналогичным образом. если мы возьмем нелинейный коэффициент b равным нулю, получим решение, соответствующее предыдущей задаче. Начнем постепенно увеличивать коэффициент b и заметим, что система начнет затухать быстрее.
 +
При одинаковой массе частиц:
 +
b=0.01
 +
b=0.1
 +
b=0.5
 +
b=1
 +
Разные массы частиц (m1=1, m2=1.3):
 +
b=0.01
 +
b=0.1
 +
b=0.5
 +
b=1

Версия 01:22, 24 декабря 2018

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

Рассмотреть перераспределение энергии между степенями свободы в нелинейной двухатомной цепочке, построить графики зависимости энергии частиц от времени.

Решение

Рассмотрим модель колебаний одномерной двухатомной цепочки массами m1 и m2. Пусть в этой цепочке находится N атомов. Обозначим смещение n-го атома un, а атома, отстоящего от него на p узлов, – un+p. Примем в качестве положительных смещения атомов вправо от положения равновесия, а отрицательных – влево (рис.1).

Ри1111.png

Каждый атом смещается только вдоль цепочки, что следует из требования одномерности модели. Пусть атомы связаны между собой упругой силой F с коэффициентом упругости с. Найдем уравнение движения n-го и n+1-го атома в цепи. В равновесном положении силы, действующие на атомы, равны нулю. При произвольных смещениях на каждый n-й атом будет действовать сила со стороны соседних атомов. В соответствии с элементарным законом Гука эту силу можно представить в виде:

Формула1.jpg

Тогда суммарная сила, действующая на n-й атом со стороны соседних атомов, будет:

Формула2.jpg

Уравнение движения n-го атома под действием силы F_n выглядит следующим образом:

Формула3.jpg

Аналогичное уравнение записывается для частиц с массой m1. Таким образом получим систему уравнений:

Формула4.jpg

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

Формула5.jpg

От системы (4) с начальными и граничными условиями (5) мы перешли к системе (6, 7):

Формула6.jpg
Формула7.jpg

Система (6) решалась в Matlab методом конечных элементов Эйлера. В задаче на каждом шаге по времени находилась кинетическая энергия всей системы, после чего был построен график зависимости энергии системы от времени. Для частиц одинаковой массы был получен следующий график:

Если рассматривать частицы разной массы, то график зависимости будет иметь следующий вид (m1=1, m2=1.3):

Рассмотрим нелинейную постановку задачи. В выражении для силы будет присутствовать слагаемое третьего порядка:

Формула8.jpg

Решается данная задача в нелинейной постановке аналогичным образом. если мы возьмем нелинейный коэффициент b равным нулю, получим решение, соответствующее предыдущей задаче. Начнем постепенно увеличивать коэффициент b и заметим, что система начнет затухать быстрее. При одинаковой массе частиц: b=0.01 b=0.1 b=0.5 b=1 Разные массы частиц (m1=1, m2=1.3): b=0.01 b=0.1 b=0.5 b=1