CC BY
92
УДК 539.3
НЕУПРУГОЕ ДЕФОРМИРОВАНИЕ И ПРОЧНОСТЬ ОДНОНАПРАВЛЕННЫХ КОМПОЗИТОВ ПРИ ПРОДОЛЬНОМ СДВИГЕ
Аношкин А.Н.(Пермь)
Abstract
The mathematical model for investigation of the nonelastic and fracture processes in unidirectional fiber-reinforced composites under longitudinal shear loading is presented. The stress and strain fields in fibers and matrix are calculated by method of local approximation from the solution of antiplane boundary-value problems for structural fragments with tetragonal fibers arrangement. The matrix constitutive equations include formulas of deformation theory of plastisity and shear modulas reduction relations for zones where a strength criterion is violated. Fibers were elastic and nonfructured. The boundary forces for structural fragment are given by steps in accord with simply longitudinal shear loading conditions. Epoxy composites with fibers from glass, carbon, boron and kevlar were considered. Volume fraction of fibers were 0.58 and
0.72. Stress and strain fields in matrix and fibers were calculated. Macrostress and macrostrain for the composites were obtained from the averaging ones.Nonlinear macrostress-macrostrain diagrams were constructed. Composite strength was estimated from the analysis of fracture processes in matrix. The results were correlated with experimental and calculated date presented in the literature.
Исследования, проведенные рядом авторов [1,2,3], показали, что момент, когда сдвиговые напряжения в критических зонах армированного пластика достигают значений прочности связующего, не всегда является началом лавинообразного разрушения материала. Микротрещины, возникающие в полимерном связующем перпендикулярно направлению действия максимальных растягивающих напряжений, имеют возможность равновесно развиваться вследствие перераспределения нагрузки на соседние области матрицы. В работе [1] для моделирования указанного механизма использовалось представление об условном течении полимерного связующего при достижении напряжениями сдвига критического значения. Напряженно-деформированное состояние в компонентах пластика рассчитывалось по методу сечений, пренебрегая взаимодействием соседних волокон. По мере нагружения зона условной текучести распространялась по всей матрице, обуславливая значительную нелинейность диаграммы деформирования композита при продольном сдвиге.
В настоящей работе для моделирования процессов неупругого деформирования и разрушения однонаправленных композитов при продольном сдвиге используется метод локального приближения , учитывающий взаимодействие ближайших соседних волокон [4]. Напряженно-деформированное состояние в волокнах и связующем композита рассчитывается на структурном фрагменте с регулярной тетрагональной упаковкой девяти цилиндрических волокон. Математическая постановка антиплоской физически нелинейной задачи для области V, соответствующей структурному фрагменту, заключается в отыскании минимума функционала Лагранжа
81 = 21 (Т вб,:,»! |5 + Т 2:,02з «Т 2:,^ ~1Р ^ • «
V Я
где 8w - вариация перемещений в направлении оси OX3• совпадающей с направлением волокон, у , у , 8 у , 8 у - деформации сдвига и их вариации, являющиеся
ненулевыми в поставленной задаче, F - внешние силы на границе S структурного фрагмента, действующие в направлении оси OX. Многосвязная область V состоит из подобластей V(K)• соответствующих матрице (к=1) и волокнам (к=2), на границе подобластей задаются условия идеального контакта. Дополнительными условиями для функционала (|) являются геометрические уравнения
У,: = дх|"’^2: “ах^. ()
Физические соотношения для матрицы включали уравнения деформационной теории пластичности на участке активного нагружения и редуцирования жесткости, в случае невыполнения критерия прочности (Ф>0),
Т„(г) = С>-ф) У„,(гХ а =1,2 (:)
с <■> -®(®,)х ф—
Ga: I 0, Ф> 0
8 ■ = Т 2: + У 2:).
Волокна считались упругими и неразрушаемыми.В качестве критерия прочности использовался критерий Мизеса, принимающий в данной постановке вид
ф = 2Ти + г2з -1 - 0. (4)
Величину внешних сил F на границе S подбирали таким образом, чтобы напряжения Т,:, т2:, осредненные по центральной ячейке периодичности фрагмента, равнялись
* *
макроскопическим т1:, Т2: [4]. Макроскопические напряжения в данной работе
задавались по шагам в соответствии с условиями простого нагружения :
Т* =Т0 (| + ак)
Т*: Т2: , k=1,2,...,N (5)
Ъ: = 0
с шагом =1 до начала разрушения матрицы и далее с шагом =0.25 по формулам
Т2: = Т2: (I + а^ + вк)
2: : , k=1,2,...,M (6)
т = 0.
І13
Начальное значение макронапряжений т^ принималось равным 10 МПа.
Поставленная антиплоская задача решалась МКЭ с использованием конечноэлементной сетки, содержащей 1400 узлов, 2700 элементов, характер триангуляции структурного фрагмента с объемной долей волокон vf=0.58 показан на рис., .
Рис. I. Триангуляция структурного фрагмента.
Исследовалось несколько типов полимерных пластиков со следующими характеристиками компонентов :
• матрица - эпоксидная смола ЭДТ-10,0=107: МПа ,р=0.0158 МПа (что соответствует
пределу прочности на сдвиг порядка 45 МПа), функция ю(8;) рассчитывалась на
основе диаграммы т-у (рис 2,а), построенной по данным работы [5];
• матрица - эпоксифенольная смола, б =1500 МПа , р=0.0158 , линейная модель [I];
• высокомодульное стекловолокно, б=41:22 МПа [6];
• стекловолокно на основе Е-стекла, б=:0578 МПа [6];
• бороволокно, б=,60 000 МПа [6];
• углеволокно, б=47000 МПа [6];
• органоволокно, б=2400 МПа [7].
Тестирование программного комплекса проводилось путем сравнения эффективного модуля сдвига и полей напряжений в компонентах высокомодульного стеклопластика (б(2) / б(,) =:8.5 , vf=0.58 ), полученных в результате решения упругой
задачи аналитическим методом в рядах [I] и численным по разработанной программе. Сравнение показало, что в области максимальных концентраций отличие напряжений на межфазной поверхности составляло не более 7% , отличие эффективных модулей - 5%.
На рис.2 приведены диаграммы деформирования т2*: -у*2: высокомодульного
стеклопластика с объемной долей волокон 0.58 , построенные по осредненным значениям полей напряжений и деформаций в центральной ячейке периодичности структурного фрагмента. Анализ результатов показывает, что возникновение зон разрушения полимерной матрицы на начальном этапе не влияет на характер диаграмм.
Рис. 2. Диаграммы деформирования эпоксидного связующего (а) и
эпоксидного (1-3) и эпоксифенольного (4,5) стеклопластика (б) vf=0.58: 1-нелинейная модель с разрушением; 2-нелинейная модель без разрушения; 3,4- линейная модель с разрушением; 5-эксперимент [I]
При достижении размера разрушенной зоны Sr порядка !6% от площади матрицы и выше наблюдается более заметное отклонение от линейной зависимости. Сравнение различных вариантов физических соотношений (рис. 2,б) позволяет сказать, что для корректного моделирования диаграммы т*3 - ^23 композита учет возможности разрушения матрицы
посредством редуцирования жесткости в уравнениях (3) является более важным фактором, чем учет нелинейного деформирования. Однако при использовании для
матрицы нелинейной модели процесс развития поврежденных зон является более медленным (кривая I на рис.2,б лежит выше кривой 3). Сравнение результатов расчета с экспериментальной диаграммой для эпоксифенольного стеклопластика показывает, что физические соотношения для матрицы, моделирующие разрушение занулением модуля сдвига в соответствующих зонах, позволяют качественно описать характер деформирования однонаправленного композита при продольном сдвиге.
На рис.3 представлены рассчитанные диаграммы деформирования для нескольких типов эпоксидных пластиков с объемными долями волокон 0.58 и 0.72, а также экспериментальная диаграмма для углепластика из работы [8]. Видно, что с увеличением объемной доли волокон нелинейность диаграмм уменьшается. При этом увеличивается отличие диаграмм для различных типов пластиков, что согласуется с аналитическими оценками эффективных модулей продольного сдвига для однонаправленных композитов [9]. Уровень макронагрузки, соответствующей началу структурного разрушения, для боро-, угле- и стеклопластиков примерно одинаков: 20-2! МПа при vf=0.72 и 27-28 МПа при vf=0.58.
Рис.3. Диаграммы деформирования эпоксидных пластиков с объемной долей волокон 0.58 (Ь4) и 0.72 (5-9) : ^5- боропластик; 2,6- углепластик; 3,8-стеклопластик на основе Е-стекла; 4,9- органопластик; 7-высокомодульный стеклопластик; точками показана экспериментальная диаграмма для эпоксидного
углепластика [5]
Интересно отметить, что хотя разрушение в матрице композитов с объемной долей волокон 0.72 начинается при меньших нагрузках , дальнейшее увеличение разрушенной области происходит более медленными темпами, например,для стеклопластика из Е-стекла vf=0.72 при т**3 =49 МПа размер разрушенной зоны Sr=6% от площади матрицы,
а для аналогичного композита с vf=0.58 ,при Т*3=49 МПа - Sr =35%. Это обстоятельство
приводит к выводу, что с увеличением коэффициента армирования прочность при продольном сдвиге повышается. Данное заключение противоречит работам [^2], в которых отмечается, что прочность композита при продольном сдвиге не зависит от коэффициента армирования , либо уменьшается с его увеличением. Отсутствие данных по коэффициентам армирования и свойствам матриц у большинства экспериментальных результатов не позволяет провести проверку полученного вывода.
В целом диаграммы деформирования исследованных композитов при соответствующей обьемной доле являются подобными, как и характер процессов разрушения, протекающих в матрице. Исключение составляет эпоксидный органопластик. В связи с малым модулем сдвига волокон поле напряжений в структуре данного композита более однородно, поэтому области разрушения существенно увеличиваются при небольших изменениях внешней нагрузки, что приводит к значительной нелинейности диаграмм деформирования. Увеличение коэффициента армирования, в
^ * * связи с вышеуказанным обстоятельством мало влияет на характер диаграмм т23 - у23 и на
уровень нагрузок, соответствующих началу разрушения матрицы: Т*3 =33.7 МПа, vf =0.72 и т*3=34.0 МПа , vf=0.58.
Прочность композита предлагается оценивать анализируя характер развития областей разрушения : можно считать, что композит утратил несущую способность, если разрушенная область расположена вдоль всей стороны ячейки периодичности либо по всей межфазной поверхности. Для высокомодульного стеклопластика с vf=0.58 такая ситуация возникает при т2*3=65 МПа, при этом размер разрушенной зоны составляет
порядка 45% от площади матрицы. Полученная оценка попадает в диапазон экспериментальных значений сдвиговой прочности однонаправленных стеклопластиков 57-76 МПа с коэффициентами армирования 0.5-0.62, приведенный по данным различных источников в работе [1]. Оценка прочности композитов предложенным способом дает для боро-, угле- и стеклопластика приблизительно одинаковый уровень предельных напряжений порядка 63 МПа при vf=0.58 и порядка 120 МПа при vf=0.72. Для органопластика в качестве оценки прочности можно взять уровни нагрузки 40 МПа (vf=0.58) и 48 МПа (vf=0.72)• предшествующие скачкообразному увеличению разрушенной области с 16 до 40% при vf=0.58 и с 7% до 38% при vf=0.72. Приведённые оценки согласуются с диапазонами экспериментальных данных : органопластик - 15-37 МПа [10], 43-46 МПа [3], 20-44 МПа [6], 60-62 МПа [11]; углепластик - 43-53 МПа [11] и 67-127 МПа[3]; боропластик - 46 МПа [12], 60 МПа [6].
Заключение.
Моделирование структурных процессов неупругого деформирования и разрушения на основе уточненных расчетов полей напряжений и микродеформаций по разработанной модели показало возможность описания нелинейных диаграмм деформирования полимерных композитов при продольном сдвиге. Процесс развития зон разрушения в матрице является устойчивым. Сдвиговая прочность композита, прогнозируемая на основе анализа численных фрактограмм, соответствует экспериментальным результатам.
Литература
1. 1.Скудра А.М., Булавс Ф.Я. Структурная теория армированных пластиков. -Рига.:Зинатне ,1978. - 192с.
2. Скудра А.М., Булавс Ф.Я. Прочность армированных пластиков. - М.:Химия,1982. -213с.
3. Арнаутов А.К. Перспективные методы испытаний пространственно армирован-ных композитов на сдвиг// Механика композит. материалов.-1990.- N5.- С.891-898.
4. Соколкин Ю.В., Ташкинов А.А. Механика деформирования и разрушения структурно-неоднородных тел. - М.:Наука, 1984. - 116 с.
5. Удрис АО., Упитис З.Т. Экспериментальное исследование упругих и прочностных свойств эпоксидного связующего ЭДТ-10 в условиях сложного напряженного состояния //Механика композит. материалов.-1988.-№.- С.972-978.
6. Композиционные материалы :Справочник / В.В.Васильев, В.Д.Протасов, В .В. Болотин и др. - М.: Машиностроение, 1990. - 512 с.
7. Максимов Р.Д., Плуме Э.З., Соколов Е.Д. Упругость высокопрочного органичес-кого волокна и органопластика // Механика композит. материалов.- 1988.- N2.- С.211-220.
8. Чамис К. Микромеханические теории прочности // Композиционные материалы. -Т.5. -1978.- С.106...165
9. Победря Б.Е. Механика композиционных материалов. - М.: Изд-во Моск. ун-та, 1984. - 336с.
10. 10.Мешков Е.В., Кулик В.И., Нилов А.С.,Упитис З.Т. Влияние технологических факторов на механические характеристики однонаправленного органопластика //Механика композит. материалов. 1990.- N3.- С.528-536.
11. 11.Хитров В.В., Лапоткин В. А., Суханов А.В. Несущая способность многослойных составных трубчатых стержней из композитов // Механика композит. материалов -1990.- N1.- С.85-92.
12. 12.Тарнопольский Ю.М., Кинцис Т.Я. Сравнительная оценка методов испытаний композитов на сдвиг // Прочность и разрушение композитных материалов.-Рига:Зинатне,1983.- С.75-85.
