CC BY
19
СТРУКТУРА ВЕЩЕСТВА И ТЕОРИЯ ХИМИЧЕСКИХ ПРОЦЕССОВ
УДК 539.3
Е. В. Артюх, Р. А. Каюмов, И. З. Мухамедова ПРИБЛИЖЕННОЕ АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ОДНОЙ ПЛОСКОЙ ЗАДАЧИ НАЛОЖЕНИЯ БОЛЬШИХ ДЕФОРМАЦИЙ ДЛЯ НЕКОТОРЫХ МОДЕЛЕЙ НЕЛИНЕЙНО-УПРУГИХ МАТЕРИАЛОВ
Ключевые слова: тензор напряжений, деформация, нелинейно-упругий материал, полость, контур.
Рассматривается решение плоской задачи об образовании круговой полости в предварительно нагруженном бесконечно протяженном теле из нелинейно-упругого материала средствами компьютерной алгебры на основе пакета «Mathematica». Учитывается, что образование полости приводит к перераспределению в теле конечных деформаций. Форма полости считается заданной в момент ее образования. Свойства материала описываются апробированными соотношениями, предложенными различными исследователями. Постановка задачи осуществляется на основе теории многократного наложения больших деформаций. Приводится приближенное аналитическое решение задачи и его сравнение с известным точным решением для частного случая.
Keywords: stress tensor, deformation, nonlinear elastic material, cavity, circuit.
The decision of using computer algebra package based on «Mathematica» plane problem of educa-tion circular cavity in a preloaded infinitely extended body of the non-linear elastic material. Taken into account that the formation of the cavity leads to a redistribution in the solid of finite deformations. The shape of the cavity defined at the moment of its formation. The material properties are described untested relations proposed by various researchers. Statement of the problem is based on the theory of multiple overlay large deformations. An approximate analytical solution and its comparison with known exact solution for a particular
Введение
Механическая постановка задачи следующая. Пусть в нелинейно-упругом бесконечно протяженном теле, находящемся в начальном (ненапряженном) состоянии, под воздействием внешних нагрузок возникли большие плоские статические деформации. Тело перешло в первое (промежуточное) состояние (рис. 1). Далее в этом теле намечается замкнутый круговой контур (будущая граница полости). Затем область, ограниченная данным контуром, удаляется. Под удалением, например, можно понимать «откол» одной части от другой или изменение свойств «удаляемой» части тела таким образом, что она не взаимодействует с оставшейся частью тела. Действие удаленной части тела на оставшуюся заменяется по принципу освобождаемости от связей силами, распределенными по данному контуру. Это не приводит к изменению напряженно -деформированного состояния в оставшейся части тела. Далее предполагается, что эти силы, перешедшие в разряд внешних, квазистатически значительно изменяются, что вызывает появление в оставшейся части тела дополнительных больших (по крайней мере, в окрестности вновь образованной граничной поверхности) деформаций, которые накладываются на уже имеющиеся большие начальные. Естественно меняется и форма граничного контура полости. Тело переходит в конечное состояние.
Математическая постановка данной краевой задачи записывается в координатах первого (промежуточного) состояния, т.к. в этом состоянии известна форма отверстия.
Уравнение равновесия
i
V-
(1 + Д0,1 ) ¿0,2 -*1,2
= 0
(1)
о
m
444-
о
+++
Рис. 1 - Постановка задачи (схема нагружения). 0 - начальное состояние (в теле отсутствуют напряжения и деформации); 1 - промежуточное состояние; 2 - конечное состояние, с11, сг22, ст12 -начальные напряжения на бесконечности
Граничные условия
i i
N2 -¿0,2
t =-p (i+д0д ))-4v1
(2)
(3)
Условие несжимаемости (для несжимаемых материалов)
(4)
(5)
Д0,2 = 0.
Геометрические соотношения:
1 + Д 0,2 = det W 0,2,
(7
22
21
а
11
12
^0,2 =^од (6) 1
^1,2 = I+ У п2. (7)
Поставленная задача (1)-(7) решена для различных нелинейно-упругих материалов, определяющие соотношения для которых описаны известными исследователями.
Например, материал Мурнагана определяется следующим соотношением:
о ( о Л о (о Л2
Ео,„ = А| Ео,„• -III + 2ОЕо,„ + 3С31 Ео,„• -I I I +
+С4 ^Ео,п •1| I + 2С4 ^Ео,п - I^ Ео,п + 3С5 ^£о,п ^ , (8)
Определяющее соотношение для материала Муни имеет вид:
М
{(1 + Р)Ъп + (1 - Р)[(Ро,п • -I)^ - Ро2п]} - Ро,
I
Материал Трелоара определяется как:
С,п =МРо,п - Ро,nI (9)
Материал Черныха имеет вид:
С = М {{1 + Р)С + (1 - Р)РоЛ - Ро,п! , (1о)
Материал Бартенева-Хазановича определяется как:
С = 2МР,1П2 - Р0nnI , (11)
Определяющее соотношение для материала Валани-са-Ландела имеет вид:
См = 2С 1п С - РоJ , (12)
Материал Муни-Ривлина имеет вид:
См = 2С1Ро,п - 2^2Р-П + Сз<П2 - РоJ . (13)
Решение поставленной нелинейной задачи строилось на основе метода Синьорини. Для этого был разработан алгоритм и написана его программная реализация на внутреннем языке системы компьютерной алгебры «МаШешаИса», что позволило получить приближенное аналитическое решение задачи (расчеты выполнены для нулевого и первого приближений) [2]. Решение задачи - характеристики напряженно-деформированного состояния тела для различных типов материала - представлены элементарными функциями координат, параметров материала и нагружения [3].
Далее приводятся некоторые полученные результаты для материала Мурнагана. Определяются выражения для главных значений тензора полных истинных напряжений как функции полярных координатах при всестороннем растяжении с учетом давления на границе полости.
С =-4 + Ч(Р\-1) -(Р + ч)2(Р-1) X(О + А-2С4 -3С5), (14)
1 рг рг 4О (2О + Х)р4 45
Р+ Ч(Р +1) +_ (Р + Ч)2
-х[р2(О + А-2С4 -3С5) +
- р р 4О (2О + А)р4
+13О + 5А + 6С4 + 9С5 ], (15) Здесь р - давление на контуре; ч - нагрузка по осям х и у ; А, О, С3, С4 и С5 - параметры материала Мурнагана.
Получены соотношения для компонент сги и с22 тензора полных истинных напряжений как функции в декартовых координатах при растяжении вдоль оси без давления на границе полости. Построены соотношения для компонент тензора полных деформаций для случая плоской деформации и одноосного растяжения по оси х нагрузкой ч на бесконечности как функции декартовых координат.
При растяжении тела из материала Мурна-гана вдоль оси х нагрузкой ч = о.3О на бесконечности в точке контура полости с максимальной концентрацией напряжений посчитаны главные значения тензора полных деформаций. Они равны: Е1 = о.25 , Е2 = -о.13 .
В точке контура полости тела из материала Мурнагана с максимальной концентрацией напряжений интенсивность деформаций для случая одноосного растяжения по оси х нагрузкой ч = о.3О достигает 22%.
Аналогичные выражения и результаты получены и для других материалов.
На рис. 2 приведен контур полости до и после деформирования тела из материала Мурнагана при одноосном нагружении с давлением на границе полости и без него.
0.5
О . 5
Рис. 2 - Материал Мурнагана. А = 2.1О, С3 = -о.о7О, С4 = -о.38О, С5 = о.34О. Растяжение вдоль оси х. Контур полости до и после деформирования. а) ч = о.3О, р = о ; б) ч = о.2О, р = о.1О
Из приведенных выше графиков виден качественный результат учета нелинейных эффектов в точках контура [4].
На рис. 3 приводится сравнение полученных результатов с известным точным решением задачи для случая предварительного всестороннего
нагружения при плоской деформации материала Трелоара [5]. Представлены зависимости контурных напряжений от нагрузки на бесконечности. Отметим, что при | q | / / < 0.6 погрешность решения задачи методом Синьорини не превышает 9% для напряжений, при этом для перемещений получаются более точные результаты, чем для напряжений.
Рассматривая напряжения как функции начальных нагрузок, можно непосредственно до начала численных или численно-аналитических расчетов указать диапазон нагрузок, при которых такой подход невозможен.
Полученные результаты, в частности, позволяют в аналитической форме использовать нелокальные критерии прочности, давать предварительно оценку диапазона констант, при которых можно предполагать практическую сходимость численных методов в аналогичных задачах.
с //
Рис. 3 - Материал Трелоара. Зависимость контурных напряжений от нагрузки на бесконечности
Авторы благодарят профессора В.А.Левина (Московский государственный университет им. М.В.Ломоносова) за внимание к поставленной задачи и консультации при ее решении.
Работа выполнена при поддержке Министерства образования и науки РФ в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 20092013 годы (соглашение 14.B37.21.2018).
Литература
1. Левин В.А. Многократное наложение больших деформаций в упругих и вязкоупругих телах. - М.: Наука, 1999. - 224 с.
2. Левин В.А. Расчет напряженно деформированного состояния вблизи повреждения, возникающего в предварительно нагруженном теле из несжимаемого материала, с помощью программного комплекса для аналитических вычислений на ЭВМ «Наложение» /В.А.Левин, Е.В.Рыбалка // Пятнадцатый симпозиум «Проблемы шин и резинокордных композитов». - Т.2. - ООО «НТЦ «НИИШП», Москва, 2004. - С. 40-45.
3. Левин В.А. О приближенном аналитическом решении задачи об образовании кругового концентратора напряжений / В.А.Левин, Е.В.Рыбалка // Известия Тульского государственного университета. Серия «Математика. Механика. Информатика». - 2004. - Т.10. - Вып.2. - Тула: ТулГУ. - С. 131-137.
4. Рыбалка Е.В. Анализ напряженно-деформированного состояния вблизи кругового концентратора напряжений, образованного в нагруженном теле из нелинейно-упругого материала. Конечные деформации // Е.В.Рыбалка // Известия Тульского государственного университета. Серия «Дифференциальные уравнения и прикладные задачи». - 2005. - Вып.1. - Тула: ТулГУ. -С. 233-249.
5. Левин В.А., Зингерман К.М. Плоские задачи многократного наложения больших деформаций. Методы решения. - М.: ФИЗМАТЛИТ, 2002. - 272с.
6. Каюмов Р.А. Моделирование поведения пленочно-тканевого материала при воздействии эксплуатационных факторов/ Р.А. Каюмов, А.М.Сулейманов, И.З.Мухамедова // Механика композиционных материалов и конструкций. 2005. - т.11. - №4. - С.519-530.
7. Каюмов, Р.А. Физически нелинейная модель тканевых материалов, обработанных в потоке высокочастотной плазма/ Р. А. Каюмов, И. Ш. Абдуллин, М.Н.Серазутдинов, В.В.Хамматова, И.З.Мухамедова, Е.В.Кумпан, А.М,Шакирова // Вестник Казанского технологического университета. № 6. 2010. - С. 42 - 53.
8. Каюмов, Р.А. Математическое моделирование деформирования текстильных материалов с содержанием полимерных волокон /Р.А.Каюмов, В.В.Хамматова// Вестник Казанского технологического университета. № 14; 2012. - С. 154 - 156.
© Е. В. Артюх- к.ф.-м.н, доцент, ФГБОУ ВПО «Тульский государственный университет» [email protected], И. З. Мухамедова - к.ф.-м.н., доц. каф. сопротивления материалов и основ теории упругости КГАСУ, тиИате(1(эта[email protected]; Р. А. Каюмов - проф., д.ф.-м.н., зав. каф. сопротивления материалов и основ теории упругости КГАСУ, проф. каф. дизайна КНИТУ, [email protected].
