Редактирование: Модель растяжения материала

==Цель исследования==
Исследование коэффициента Пуассона модельного материала
==Дано:==
#Модель материала - множество частиц, взаимодействующих друг с другом по определённому закону
#Программа визуализации расчетных данных
==Задача №1:==
#Сымитировать растяжение материала вдоль оси У
#Сымитировать то же растяжение, только под действием фиксированной силы.
<gallery widths=800px heights=341px perrow = 1>
Файл:KP_1_Task.jpg
</gallery>

==Как взаимодействуют частицы== 
Взаимодействие частиц описывается потенциалом Леннарда-Джонса который записывается в следующем виде:<br>
<math>U(r) = D \left[ \left(\frac{a}{r}\right)^{12} - \left(\frac{a}{r}\right)^{6} \right],</math><br>
r — расстояние между центрами частиц<br>
D  — глубина потенциальной ямы<br>
a — равновесное расстояние<br>
Параметры  D и a являются характеристиками вещества. Характерный вид 
потенциала показан на рисунке, его минимум лежит в точке <math>r_{min} = a</math> <br>
<gallery widths=400px heights=283px perrow = 1>
Файл:Lennard_Jones_potential.jpg
</gallery>

==Реализация на компьютере==
===Расчет силы взаимодействия частиц===
<math>F = -U'(r)</math><br>
<math>\b{F} =  f_0\left[ \left(\frac{a}{r}\right)^{14} - \left(\frac{a}{r}\right)^{8}\right]\b{r}</math><br>
<math>f_0 = 12\left(\frac{D}{a}\right)</math><br>
Где D=1 , a=1 для данного материала , а <math>\b{r}</math> – радиус-вектор, соединяющий 2 частицы
===Расчет ускорения частицы (по II закону Ньютона)===
<math>\b{a} = \b{F}/m</math><br>
m – масса частицы, для данного материала m=1<br>
===Расчет скорости частицы===
<math>\underline{\vartheta} = \underline{\vartheta}_{pre} +\underline{a}\cdot dt</math><br>
<math>dt = \frac{1.45 \cdot T_0}{100}</math> - шаг интегрирования<br>
<math>T_0 = 2 \pi \sqrt{\frac{a}{6\cdot f_0}}</math> - период колебаний частиц<br>
<math>f_0</math> – характерная сила межатомного взаимодействия<br>
===Расчет новых координат частицы===
<math>\underline{r} = \underline{r}_{pre} + \underline{\vartheta} \cdot dt</math><br>
==Задача №2==
Измерить коэффициент Пуассона для полученной модели<br>
'''Коэффициент Пуассона''' - одна из физических характеристик материала упругого тела, равная отношению абсолютных значений относительной поперечной деформации элемента тела к его относительной продольной деформации.<br>
Величина коэффициента Пуассона для большинства металлических материалов близка к 0,3.<br>
<gallery widths=556px heights=320px perrow = 1>
Файл:Deformations_KP.jpg
</gallery>
За начальное состояние материала взято его состояние на 1000м шаге. До него сила, растяжения материала = 0. Таким образом система уравновешивается. А уже начиная с этого шага мы нагружаем систему.

==Результат: коэффициент Пуассона для модельного материала без пористости==
В результате проделанной работы было определено, что для данной модели коэффициент Пуассона получился равным 0.33 ,что близко к результату, полученному аналитическим путем <math>\nu = \frac{1}{3}</math>  <br>
[1]Кривцов А. М. Деформирование и разрушение твердых тел с микроструктурой.  М.: ФИЗМАТЛИТ,  (2007)<br>

==Другой способ расчета коэффициента Пуассона==
Другой способ измерения основан на измерении напряжений возникающих в материале в напряженно-деформированном состоянии. И здесь вместо растяжения мы будем сжимать материал. Этим исключается искажение результатов вызванное разрывом материала. Также используются периодичные граничные условия (мы окружаем наш кусочек материала со всех сторон идентичными ему кусочками материала т.е. как только частица улетела вправо за границу рассматриваемой области то она прилетает с левой границы) таким образом создается эффект рассмотрения маленького участка вместо рассмотрения маленького кусочка материала как в предыдущем методе. Здесь мы деформируем материал не динамически прикладывая силу а уже изначально создаем деформированный материал зажатый границами рассматриваемой области.<br>
В этом состоянии мы измеряем напряжения возникшие в рассматриваемом материале по вертикали и по горизонтали.<br>
<math>C_e = \frac {F_e} {S}</math><br>
<math>C_e</math> — механическое напряжение вдоль оси <math>\underline{e}</math><br>
<math>F_e</math> — сила, возникшая в теле вдоль оси <math>\underline{e}</math> при деформации <br>
S — площадь.
<math>\nu = \frac{C_x}{C_y}</math><br>
==Результаты==
===Без пористости===
Используя этот метод удалось с высокой точностью (до <math>10_{-6}</math>) получить значение коэффициента Пуассона для данного материала.
<math>\nu =0.333333 \approx \frac{1}{3}</math>  <br>
Что соответствует аналитическому результату.
===Зависимость коэффициента Пуассона от пористости материала===
В материале было задано случайное распределение пустот(извлечен случайный набор частиц) и для разного процента пористости был измерен коэффициент Пуассона. Для каждого значения процента пористости делалось несколько измерений коэффициента Пуассона для разного распределения пустот. В результате чего был построен '''График Зависимости коэффициента Пуассона от пористости материала'''<br>
<gallery widths=800px heights=417px perrow = 1>
Файл:Porous_KP.jpg
</gallery>
Как видно зависимость возрастающая что само по себе странно т.к. все экспериментальные(для пенобетона) и аналитические данные говорят обратное что зависимость должна быть убывающая. В результате долгих поисков ошибка не обнаружилась и было решено измерить зависимость модуля Юнга от пористости таким же способом и проверить что будет с ней...

===Зависимость модуля Юнга от пористости материала===
<gallery widths=800px heights=417px perrow = 1>
Файл:Porous_mJ.jpg
</gallery>