Определение упругих модулей материала — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Введение)
(Компьютерный эксперимент с конкретным материалом)
 
(не показано 27 промежуточных версий 4 участников)
Строка 1: Строка 1:
 +
 +
'''''Курсовой проект по [[Механика дискретных сред|Механике дискретных сред]]'''''
 +
 +
'''Исполнитель:''' [[Фомичева Мария]]
 +
 +
'''Группа:''' [[Группа 10|10]] (43604/1)
 +
 +
'''Семестр:''' осень 2017
 +
 
== Введение ==
 
== Введение ==
 
В настоящее время большое внимание уделяется исследованию упругих свойств различных материалов.
 
  
 
[[Файл:Материал для определения упругих модулей.png|thumb|Рис.1. Исследуемый материал]]
 
[[Файл:Материал для определения упругих модулей.png|thumb|Рис.1. Исследуемый материал]]
 +
В настоящее время большое внимание уделяется исследованию упругих свойств разных материалов. Эти исследования позволяют определить поведение материала при различных деформациях и напряжениях.
  
В данной работе проводится исследование двух упругих модулей - коэффициента Пуассона и модуля Юнга. Вычисление модулей ведется с помощью компьютерного эксперимента. Компьютерный эксперимент был поставлен на материале, изображенном на Рис.1. При вычислении модулей используется метод молекулярной динамики (ММД). Кроме того, в задаче ставятся фиксированные граничные условия.
+
В данной работе проводится исследование материала на его упругие характеристики - коэффициента Пуассона и модуля Юнга. Вычисление модулей ведется с помощью компьютерного эксперимента. Пример материала, для которого производились расчеты, показан на Рис.1. При вычислении упругих коэффициентов используется метод молекулярной динамики (ММД). Кроме того, в задаче ставятся фиксированные граничные условия: левая грань материала фиксируется, правая грань растягивается вдоль горизонтально оси, все остальные грани - свободные.
  
 
== Алгоритм компьютерного эксперимента ==
 
== Алгоритм компьютерного эксперимента ==
Строка 11: Строка 19:
 
Весь компьютерный эксперимент можно условно разделить на три этапа.
 
Весь компьютерный эксперимент можно условно разделить на три этапа.
  
''На первом этапе'' вычисления находится положение равновесия материала в растянутом состоянии.
+
''На первом этапе'' находится положение равновесия материала в растянутом состоянии.
При этом задается растяжение вдоль одной из оси симметрии решетки (оси X). На этом этапе решается динамическая задача
+
При этом задается растяжение вдоль одной из осей симметрии материала (оси X). Компьютерный эксперимент производится посредством нахождения
достижения положения равновесия. Компьютерный эксперимент производится посредством вычисления
 
 
радиус векторов и векторов скорости частиц в зависимости от времени. Интегрирование
 
радиус векторов и векторов скорости частиц в зависимости от времени. Интегрирование
ведется методом центральных разностей. Метод состоит в том, что координаты и силы вычисляются
+
ведется методом центральных разностей. Данный метод состоит в том, что координаты и силы вычисляются
во временных точках, разделенных интервалами, равными шагу интегрирования, а скорости
+
во временных точках, разделенных интервалами равными шагу интегрирования, а скорости
 
вычисляются во временных точках, находящихся в серединах вышеупомянутых интервалов:
 
вычисляются во временных точках, находящихся в серединах вышеупомянутых интервалов:
  
Строка 25: Строка 32:
 
где <math>\tau</math> – шаг интегрирования. Ускорение <math>\underline{w}(t)</math>
 
где <math>\tau</math> – шаг интегрирования. Ускорение <math>\underline{w}(t)</math>
 
вычисляется через приложенную к частице силу.
 
вычисляется через приложенную к частице силу.
 +
Кроме того, на первом этапе вычисляется средняя деформация материала после его растяжения.
  
''
+
''Второй этап'' представляет собой определение слагаемых сил, действующих на один атом системы и на
Второй этап'' представляет собой определение слагаемых сил, действующих на один атом системы и на
+
соседние с ним атомы. Зная силы, механические напряжения в решетке можно вычислить по формулам:
соседние с ним атомы. Эта часть содержит вычисление значения производных от одного
 
слагаемого потенциальной энергии системы, т.е. потенциальной энергии, приходящейся на
 
один атом.
 
 
 
Механические напряжения в решетке вычисляются по формулам:
 
  
 
<math>
 
<math>
Строка 38: Строка 41:
 
     \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i \underline{A}_{\alpha}^i =
 
     \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i \underline{A}_{\alpha}^i =
 
     \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i (\underline{r}_{\alpha}^i - \underline{r}_i),
 
     \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i (\underline{r}_{\alpha}^i - \underline{r}_i),
</math>
 
 
иначе
 
 
<math>
 
    {{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i =
 
    \frac{1}{V} \sum_{\alpha} \frac{\partial \Pi^i}{\partial \underline{A}_{\alpha}^i}
 
    (\underline{r}_{\alpha}^i - \underline{r}_i), \quad \Pi^i = \frac{1}{2} \sum_{j (\neq i)} V_{ij}
 
 
</math>
 
</math>
  
 
Здесь <math>{{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i</math> – тензор механических напряжений для частицы с номером <math>i</math>. При однородном поле деформации находится средний тензор напряжений
 
Здесь <math>{{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i</math> – тензор механических напряжений для частицы с номером <math>i</math>. При однородном поле деформации находится средний тензор напряжений
<math>({{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i)</math> по всем частицам. <math>V</math> – объем ячейки периодичности.\frac{\partial \Pi^i}{\partial \underline{A}_{\alpha}^i}</math>, где <math>\alpha</math> – номер
+
<math>({{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i)</math> по всем частицам. <math>V</math> – объем ячейки периодичности. <math>\underline{A}_{\alpha}^i</math> – вектор относительного
соседней частицы к частице с номером <math>i</math>. <math>\underline{A}_{\alpha}^i</math> – вектор относительного
 
 
положения соседней частицы: <math>\underline{A}_{\alpha}^i = \underline{r}_{\alpha}^i - \underline{r}_i</math>,
 
положения соседней частицы: <math>\underline{A}_{\alpha}^i = \underline{r}_{\alpha}^i - \underline{r}_i</math>,
 
где <math>\underline{r}_i</math>
 
где <math>\underline{r}_i</math>
 
– радиус-вектор частицы с номером <math>i</math>, <math>\underline{r}_{\alpha}^i</math> – радиус-вектор соседней
 
– радиус-вектор частицы с номером <math>i</math>, <math>\underline{r}_{\alpha}^i</math> – радиус-вектор соседней
частицы (<math>\alpha</math>). <math>V_{ij}</math> – энергия, приходящаяся на одну связь.
+
частицы (<math>\alpha</math>).
  
  
''Третий этап'' представляет собой нахождение упругих модулей через коэффициенты упругости. Для нахождения коэффициентов упругости воспользуемся формулами для их выражения через компоненты тензоров напряжения и деформации.
+
''Третий этап'' представляет собой нахождение упругих модулей через коэффициенты упругости. Для нахождения коэффициентов упругости воспользуемся формулами их выражения через компоненты тензоров напряжения и деформации.
  
 
В трехмерном материале коэффициенты упругости
 
В трехмерном материале коэффициенты упругости
Строка 73: Строка 67:
  
  
Модули упругости выражаются по формулам:
+
Модули упругости выражаются формулами:
  
 
<math>
 
<math>
 
     \nu = \frac{C_{12}}{C_{11} + C_{12}},\quad
 
     \nu = \frac{C_{12}}{C_{11} + C_{12}},\quad
 
     E = \frac{(C_{11} - C_{12}) (C_{11} + 2 C_{12})}{(C_{11} + C_{12})},
 
     E = \frac{(C_{11} - C_{12}) (C_{11} + 2 C_{12})}{(C_{11} + C_{12})},
где  
+
</math> где  
E  - модуль Юнга,  
+
<math> </math> - модуль Юнга,  
\nu - коэффициент Пуассона
+
<math>\nu </math> - коэффициент Пуассона
</math>
+
 
 +
== Компьютерный эксперимент с конкретным материалом ==
 +
 
 +
При компьютерном эксперименте был рассмотрен изотропный материал со следующими параметрами:
 +
 
 +
Расстояние между частицами - <math>d = 0.33,</math>
 +
 
 +
Количество частиц - <math> N = 8000,</math>
 +
 
 +
Радиус обрезания - <math> A_c = 1.3</math>
 +
 
 +
Масса частиц, жесткость частиц и длина ребра материала приняты равными единице.
 +
 +
При растяжении такого материала упругие модули получаются следующими:
 +
<math>E/E* = 0.808,</math> где <math> E/E*</math> - обезразмеренный модуль Юнга
 +
 
 +
<math>\nu = 0.396</math>
  
При выборе конкретного материалана основе ГЦК с расстоянием между частицами <math>d = 0.33,
+
==Ссылки==
</math> упругие модули получились следующими
+
*Автор проекта: [[ Фомичева Мария]]
<math>E = 0.926682, \nu = 0.2274
+
*[[Виртуальная лаборатория]]
</math>
 

Текущая версия на 14:10, 23 января 2018

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

Исполнитель: Фомичева Мария

Группа: 10 (43604/1)

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

Введение[править]

Рис.1. Исследуемый материал

В настоящее время большое внимание уделяется исследованию упругих свойств разных материалов. Эти исследования позволяют определить поведение материала при различных деформациях и напряжениях.

В данной работе проводится исследование материала на его упругие характеристики - коэффициента Пуассона и модуля Юнга. Вычисление модулей ведется с помощью компьютерного эксперимента. Пример материала, для которого производились расчеты, показан на Рис.1. При вычислении упругих коэффициентов используется метод молекулярной динамики (ММД). Кроме того, в задаче ставятся фиксированные граничные условия: левая грань материала фиксируется, правая грань растягивается вдоль горизонтально оси, все остальные грани - свободные.

Алгоритм компьютерного эксперимента[править]

Весь компьютерный эксперимент можно условно разделить на три этапа.

На первом этапе находится положение равновесия материала в растянутом состоянии. При этом задается растяжение вдоль одной из осей симметрии материала (оси X). Компьютерный эксперимент производится посредством нахождения радиус векторов и векторов скорости частиц в зависимости от времени. Интегрирование ведется методом центральных разностей. Данный метод состоит в том, что координаты и силы вычисляются во временных точках, разделенных интервалами равными шагу интегрирования, а скорости вычисляются во временных точках, находящихся в серединах вышеупомянутых интервалов:

[math]\underline{v} (t + \tau / 2) = \underline{v} (t - \tau / 2) + \underline{w} (t) \tau[/math]

[math]\underline{r} (t + \tau) = \underline{r} (t) + \underline{v} (t + \tau / 2) \tau,[/math]

где [math]\tau[/math] – шаг интегрирования. Ускорение [math]\underline{w}(t)[/math] вычисляется через приложенную к частице силу. Кроме того, на первом этапе вычисляется средняя деформация материала после его растяжения.

Второй этап представляет собой определение слагаемых сил, действующих на один атом системы и на соседние с ним атомы. Зная силы, механические напряжения в решетке можно вычислить по формулам:

[math] {{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i = \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i \underline{A}_{\alpha}^i = \frac{1}{2V} \sum_{\alpha} \underline{F}_{\alpha}^i (\underline{r}_{\alpha}^i - \underline{r}_i), [/math]

Здесь [math]{{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i[/math] – тензор механических напряжений для частицы с номером [math]i[/math]. При однородном поле деформации находится средний тензор напряжений [math]({{\underline{\underline{\tau\hspace{-0.5mm}}}\hspace{0.5mm}}}_i)[/math] по всем частицам. [math]V[/math] – объем ячейки периодичности. [math]\underline{A}_{\alpha}^i[/math] – вектор относительного положения соседней частицы: [math]\underline{A}_{\alpha}^i = \underline{r}_{\alpha}^i - \underline{r}_i[/math], где [math]\underline{r}_i[/math] – радиус-вектор частицы с номером [math]i[/math], [math]\underline{r}_{\alpha}^i[/math] – радиус-вектор соседней частицы ([math]\alpha[/math]).


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

В трехмерном материале коэффициенты упругости определяются через следующие выражения:

[math] \begin{array}{l} \sigma_1 = C_{11} \varepsilon_{11} + C_{12} \varepsilon_{22} + C_{12} \varepsilon_{33},\quad \sigma_2 = C_{12} \varepsilon_{11} + C_{11} \varepsilon_{22} + C_{12} \varepsilon_{33},\\ \sigma_3 = C_{12} \varepsilon_{11} + C_{12} \varepsilon_{22} + C_{11} \varepsilon_{33},\\ \tau_{12} = 2 C_{44} \varepsilon_{12},\quad \tau_{23} = 2 C_{44} \varepsilon_{23},\quad \tau_{31} = 2 C_{44} \varepsilon_{31}. \end{array} [/math]


Модули упругости выражаются формулами:

[math] \nu = \frac{C_{12}}{C_{11} + C_{12}},\quad E = \frac{(C_{11} - C_{12}) (C_{11} + 2 C_{12})}{(C_{11} + C_{12})}, [/math] где [math] E [/math] - модуль Юнга, [math]\nu [/math] - коэффициент Пуассона

Компьютерный эксперимент с конкретным материалом[править]

При компьютерном эксперименте был рассмотрен изотропный материал со следующими параметрами:

Расстояние между частицами - [math]d = 0.33,[/math]

Количество частиц - [math] N = 8000,[/math]

Радиус обрезания - [math] A_c = 1.3[/math]

Масса частиц, жесткость частиц и длина ребра материала приняты равными единице.

При растяжении такого материала упругие модули получаются следующими: [math]E/E* = 0.808,[/math] где [math] E/E*[/math] - обезразмеренный модуль Юнга

[math]\nu = 0.396[/math]

Ссылки[править]