Взаимодействие гиперупругих тел с жесткими телами

Размещено на http://www.allbest.ru/

Министерство образования и науки Российской Федерации

КАЗАНСКИЙ (ПРИВОЛЖСКИЙ) ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ

ИНСТИТУТ МАТЕМАТИКИ И МЕХАНИКИ ИМ. Н.И.ЛОБАЧЕВСКОГО

КАФЕДРА ТЕОРЕТИЧЕСКОЙ МЕХАНИКИ

Специальность: 010800.62 - механика и математическое моделирование

Выпускная квалификационная работа

(бакалаврская работа)

Взаимодействие гиперупругих тел с жесткими телами

Работа завершена: Т.Б. Финк

Научный руководитель Л.У. Султанов

заведующего кафедрой Ю.Г. Коноплев

Казань -- 2015



Оглавление

Введение

1. Определяющие соотношения гиперупругости

2. Модели гиперупругих материалов

3. Определение констант гиперупругого материала

4. Методы решения нелинейных задач

5. Построение в ANSYS

Заключение



Введение

Изучение поведения гиперупругих материалов в современной жизни имеет большое значение. Гиперупругие материалы являются важнейшими объектами исследования нелинейной теории упругости.

Гиперупругие материалы из-за своих свойств широко используются в различных отраслях. Они не зависят от скорости деформации. К таким материалам относятся резина, пена, биологические материалы, полимерные материалы.

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

Целью данной работы является введение в нелинейную теорию упругости и, создание трехмерной модели с помощью ANSYS для исследования напряженно-деформированного состояния гиперупругих тел на примере деформации кольца.



1. Определяющие соотношения гиперупругости

Рассмотрим закон Гука линейной теории упругости

,

- единичный тензор, и - константы Ламэ.

Перепишем уравнение (1.1) в виде:

где

, , - базовые изомеры.

Запишем (1.2) в потенциальном виде

(1.3)

Для гиперупругого материала нельзя прямо использовать уравнения (1.1), (1.2), (1.3). Необходимо использовать нелинейную теорию упругости. Может быть несколько причин нелинейности, например, нелинейная зависимость деформаций и перемещений, либо большие перемещения, которые можно соизмерить с размерами конструкции, либо непропорциональная связь напряжений и деформаций.

Определяющие соотношения гиперупругости записываются в виде

- удельная потенциальная энергия деформации. Для изотропного гиперупругого материала эта функция запишется в виде:

Сопряженную пару тензоровможно записать через правый тензор напряжений Генки в виде .

2. Модели гиперупругих материалов

Из соотношения (1.4) строятся модели гиперупругих материалов.

Модель Муни-Ривлина ( для несжимаемого материала)

где , - константы;

Общая форма записи Муни-Ривлина:



Главные инварианты можно записать через главные удлинения лi в четных степенях:

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

Из уравнения (10) получаем неогуковскую модель:

Модель Муни-Ривлина можно записать в расширенной форме в пяти- и девяти константном виде:

Модель третьего порядка Джеймса-Грина-Симпсона:



Модель Огдена:

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

Кирхгофа - Сен-Венана:

Генки:

Соотношения (2.3) и (2.4) обобщают закон Гука на случай больших деформаций гиперупругой среды.

Для использования (2.1),(2.3),(2.4) в уравнениях гиперупругости требуется для заданного закона движения определить и который связывает скорости изменения тензоров и в соотношении



3. Определение констант гиперупругого материала

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

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

Напряжения Коши представим в виде:

Инварианты тензора деформации:

I1 = 12 + 22 + 32

I2 = 1222 + 1232 + 2232

I3 = 122232

Для одноосного растяжения главные степени деформации в направлении перпендикулярном оси образца будут одинаковы. Учитывая уравнение несжимаемости, получим:

Тогда инварианты деформаций запишутся как:



Значения главной степени деформации подставляем в уравнение напряжения Коши:

Вычитая уравнение (3.2) из уравнения (3.1), получим главное истинное напряжение для одноосного растяжения:

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

,

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



Частные производные упругого потенциала по главным инвариантам равны: нелинейный гиперупругий тело деформация

Переписав (3.3) с учетом (3.4) , (3.5) и (3.6), получим:

Вычислив частные производные среднеквадратичного отклонения по коэффициентам Муни-Ривлина в соответствии с (3.4), получим:

Решая систему уравнений (3.7) относительно и, приходим к формулам:



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

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

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

4. Методы решения нелинейных задач

Предположим, что в k-ом состоянии известными являются:

- вектор перемещений



- тензор деформаций

Будем считать, что под действием приращений внешних массовых и поверхностных сил,перемещения и напряжения тоже получают приращения, (). При этом эти приращения считаются малыми величинами, что позволяет пренебречь квадратами градиентов приращений перемещений, т.е.

В качестве базового разрешающего уравнения используется уравнение виртуальных работ для (k+1)-го состояния

,

где

Если ввести обозначения

,

,



то после пренебрежениями квадратами градиентов приращений уравнение (4.1) преобразуется к виду

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

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

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

5. Построение в ANSYS

Рассмотрим модель гиперупругого материала при различных степенях сжимаемости. Построим трехмерную конечно-элементарную модель с помощью ANSYS с использованием конечного элемента solid185.

Дано кольцо внешний радиус, которого R1=20 мм, внутренний радиус R2=15мм, толщина h=3мм. Рассмотрим модель Муни -Ривлина для изотропного материала.

Рис.1. Гиперупругая модель

Рассматривалась только четверть кольца, так как имеется симметрия. Задаем запрет перемещений у основания по оси Y и верхнего конца по оси X. Отсутствие перемещений по оси Z соответствует плоской деформации. Подвергаем модель последовательным наложением перемещений в направлении оси Y. Данная задача геометрически нелинейная. В ANSYS включаем большие перемещения и выбираем модель Муни-Ривлина, которая описывает несжимаемые и слабосжимаемые материалы. Этот материал подходить для умеренных растяжений. Для данной задачи можно использовать крупную сетку для простоты вычислений.

Рис.2.Начальное состояние



Возьмем модель Муни - Ривлина с двумя постоянными ,- имеют смысл модуля высокоэластичности, описывают зависимость напряжения и деформации, -почти несжимаемый, то есть объем не меняется, d=0,0042. Где d-это параметр сжимаемости материала.

Для изотропных материалов , то есть для материалов имеющих одинаковые свойства во всех направлениях.

Рис.3. Деформированное состояние четверти кольца

Прикладываемое перемещение, деформировала элемент в сторону растяжения кольца (рис.3.). Отметим, что задача является линейной в самом начале деформирования. Далее начинает проявляться геометрическая нелинейность.



Рис.4. Распределение поля напряжений

На рис.4. видно, что максимальные напряжения появляются у основания тела, особенно на внутреннем изгибе.

Рис.5.Перемещения по x



Рис.6.Перемещения по Y

Рис.7. Суммарные перемещения

Максимальное перемещение на величину 11,956 мм совершает верхний конец модели, при удалении от которого значения перемещений уменьшаются.

Построим эту же модель, при ,и , d=0,323. При малейшем изменении , имеем большую разницу в значениях d. При увеличении или увеличивается жесткость материала. Это значит, что материал увеличивает свою способность сопротивляться деформациям, пытаясь сохранить свою геометрическую форму при внешних воздействиях. Зависит от таких параметров как модуль упругости и модуль сдвига.

Рис.8. Перемещение по x

Рис.9. Перемещение по y



Рис.10. Суммарные перемещения

Рис.11. Распределение поля напряжений



Заключение

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

Задача решалась шаговым методом с учетом геометрической нелинейности, после каждого шага определяли положение координат узлов в пространстве. На каждом шаге переходили к новой расчетной схеме. На каждой ступени прикладывалось вертикальное перемещение.

Размещено на Allbest.ru