УДК 539.3
Р. А. Каюмов, А. Р. Мангушева
ПРЕДЕЛЬНЫЙ АНАЛИЗ ДЛЯ ПЛЕНОЧНО-ТКАНЕВОГО КОМПОЗИЦИОННОГО МАТЕРИАЛА
Ключевые слова: пластичность; упрочнение; большие перемещения; предельное равновесие.
Рассмотрен пленочно-тканевый композиционный материал, для которого построена трехмерная конечно-элементная модель представительной ячейки. Приведены методика решения задачи об ее двухосном растяжении и задачи построения для нее предельной поверхности. Результаты представлены в виде графиков.
Key words: plasticity; hardening; large displacements; limiting balance.
The fabric composite material is considered. For its representative cell the 3-D finite-element model is developed. The technique of the task decision about its two axes a lengthening and task to construct a limiting surface is given. The results are shown as the diagrams.
Введение
Рассматривается ячейка пленочно-тканевого композиционного материала (ПТКМ), которая представляет собой ткань из высокопрочных синтетических нитей и пленочного полимерного покрытия. Последние служат обычно для защиты армирующей основы от воздействия атмосферных факторов и придания воздухонепроницаемости материалу. Для анализа ПТКМ можно рассмотреть его представительный элемент, изображенный на рисунке 1.
Рис. 1 - Элементарная ячейка пленочно-тканевого композита
Линии АB, CD и часть эллипса описываются функциями:
7Х
/1 = aCos(—) + сх + а + ёШ / 2 а
Х
7Х
/2 = аСоз(—) + 02 + а - ёпН / 2 а
Х
г3 = + (2 - С0)2 ,
с2 .dn.it.2
у (-)
2
где c0 =
f 1 + f 2 . 2
Эти функции позволяют выделить область, занимаемую нитью ткани, которую называем основой. Поперечную нить назовем утком. Нити на других сторонах представительного элемента описываются аналогично.
Для дискретизации области по пространственным координатам используем метод конечных элементов. Для вычисления параметра предельной нагрузки применим метод вариации упругих характеристик [1], который основан на законе пластического течения.
Основная часть
Рассмотрим задачу построения предельных кривых в пространстве усилий растяжения Px и Ру. Подобные задачи рассматривались в [2] для более простых случаев ячеек, а именно для слоистых элементов оболочек. Ввиду простоты строения ячейки в этом случае авторам удалось получить систему параметрических уравнений предельной поверхности. Однако при этом в большинстве частных случаев между этими параметрами имеется нелинейная связь. Сравнения экспериментальных данных и результатов, полученных по теории предельного равновесия для слоистых образцов, приведены, в частности, в работах [3,4].
Рассмотрим частный случай условия текучести, являющийся обобщением критерия Мизеса на случай ортропного материала:
f (a) = aTAa
(4)
Матрицу А для него можно представить в
виде:
(1) An A12 A13 0 0 0
A21 A22 A23 0 0 0
(2) A = A31 A32 A33 0 0 0
0 0 0 A44 0 0
(3) 0 0 0 0 A55 0
0 0 0 0 A66
При рассмотрении соотношений пластического течения в форме:
закона
а -
А-Ч
£ = е
(5)
можно увидеть аналогию между нелинеино-упругои и жесткопластическоИ задачами. Для этого под 4 нужно понимать приращение деформации, а под матрицей упругих характеристик - матрицу
Е -
А-1
(6)
Этот подход, предложенный в [5] для условия типа (4), в дальнейшем был развит на более общие условия текучести. Для изотропных тел и условия Мизе-са аналогичная методика была рассмотрена в [6].
На основе вышеизложенного разработана методика решения полученной задачи методом итерационного уточнения матрицы Е на каждом шаге по перемещениям. Были использованы следующие исходные данные (все данные в миллиметрах): а=0.4, 4^=0.2, С2=0.3, еу=0.2, Сх=0.2, ах=1, ау=1, а2=1 (см. рис. 1). Предел текучести матрицы
- Цу[3кн /
мм
основы
,2
утка
Тоси - тутка - 10/^/зКН/мм2. Матрица А имела вид:
А-
1.34615 0.576923 0.576923 0 0 0
0.576923 1.34615 0.576923 0 0 0
0.576923 0.576923 1.34615 0 0 0
0 0 0 0.333 0 0
0 0 0 0 0.333 0
0 0 0 0 0 0.333
Здесь Тт - предел текучести при сдвиге.
На рисунке 2 приведены предельные кривые в пространстве усилий растяжения в зависимости от геометрических характеристик ПТКМ.
><*■-*----
Л
\
I1
л )
-<1пй=0,2 -<1гй=0:25 -«=0,35 -<1гй=а:3
Рх
Рис. 2 - Зависимость предельных нагрузок от толщины нити и утка
Для определения реальной несущей способности конструкции необходимо учитывать ее поведение после начала пластического деформирования, поскольку в одних случаях происходит уменьшение
предельной нагрузки, необходимой для дальнейшей деформации конструкции (что приводит к катастрофическому разрушению), а в других - ее увеличение. Такой анализ можно провести при учете конечности перемещений. При решении задач по теории предельного равновесия, как было отмечено выше, для этого можно использовать метод пересчета геометрии конструкции по известному полю перемещений. Результаты данной методики приведены на рисунке 3.
Рис. 3 - Зависимость предельной нагрузки от деформации ячейки 8у в направлении утка (оси у)
Из рисунка 3 видно, что при некотором значении удлинения ячейки, предельная нагрузка достигает своего максимума (ввиду выпрямления нитей), после чего наступает исчерпание несущей способности рассматриваемого элемента.
Поскольку для большинства материалов при больших деформациях появляется упрочнение, то необходимо учитывать изменения предела текучести. В данной работе используется теория изотропного упрочнения, тогда соотношение примем в
виде ат - ат (1 + С * ), где - интенсивность
деформации, ат0 - предел текучести. Результаты приведены на рис.4.
Рис. 4 - Зависимость предельной нагрузки от деформации ячейки 8у при разных параметрах упрочнения С*
и
3
Выводы
Рассмотрена задача оценки несущей способности ячейки ПТКМ при двухосном растяжении. Для вычисления параметра предельной нагрузки применяется метод вариации упругих характеристик, а для дискретизации области - метод конечных элементов. Получены предельные кривые в пространстве усилий растяжения Рх и Ру и зависимости предельной нагрузки от общего удлинении стержня по оси х и у.
Литература
1. Каюмов, Р.А. Об одном методе двусторонней оценки предельной нагрузки / Р.А. Каюмов // Проблемы прочности. - 1992. - N1. - С.51-55.
2. Терегулов, И.Г. Предельное состояние многослойных композитных пластин и оболочек / И.Г. Терегулов, Э.С. Сибгатуллин, О.А. Маркин // Механика композитных материалов. - 1988. - N4. - С.715-720.
3. Терегулов, И.Г Расчет конструкций по теории предельного равновесия / И.Г. Терегулов, Р.А. Каюмов, Сибгатуллин Э.С. Казань: Фэн, 2003. - 180 с.
4. Каюмов, Р.А. Прямые и обратные задачи расчета слоистых оболочечных конструкций /Р.А. Каюмов, С.В. Гусев, Р.О. Нежданов// Казань: Каз.гос.энерг.ун-т. -2004.- 180 с.
5. Каюмов Р.А. Метод вариации упругих характеристик в задаче о предельной нагрузке // Жур. ПМТФ. - 1990. - № 3. - С.134-139.
6. Рутман Ю.Л. Метод псевдожесткостей для решения задач о предельном равновесии жесткопластических конструкций. - С.Петербург. - 1998. - 53 с.
7. Каюмов, Р. А. Физически нелинейная модель тканевых материалов, обработанных в потоке высокочастотной плазма/ Р.А. Каюмов, И.Ш. Абдуллин, М.Н. Серазутдинов, В.В. Хамматова, И.З. Мухамедова, Е.В. Кумпан, А.М, Шакирова // Вестник Казанского технологического университета. №6; Федеральное агентство по образованию, Казан. гос. технол. ун-т. -Казань: КГТУ, 2010. - С. 42 - 53.
© Р. А. Каюмов - проф., д.ф.-м.н., зав. каф. сопротивления материалов и основ теории упругости КГАСУ, [email protected]; А. Р. Мангушева - к.ф.-м. н., доц. каф. интеллектуальных систем и управления информационными ресурсами КНИТУ, [email protected].