CC BY
95
НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408
электронный научно-технический журнал
Моделирование поверхностей прочности композитов на основе микроструктурного конечно-элементного анализа
# 11, ноябрь 2012 Б01: 10.7463/1112.0496336
Димитриенко Ю. И., Соколов А. П., Шпакова Ю. В., Юрин Ю. В.
УДК 539.3
Россия, МГТУ им. Н.Э. Баумана dimit.bmstu@gmail.com alsokolo@bmstu.ru shpakovayuliya@bmstu.ru ven8469@ya.ru
Введение. Современные композиционные материалы, как правило, обладают явно выраженной многомасштабной иерархической структурой, в которой каждый последующий структурный уровень, содержит предыдущий уровень. Особенно четко такая структура реализуется в композитах на основе армирующих волокон различного переплетения: тканевого, намоточного, пространственного, у которых сами волокна представляют собой пучки большого числа моноволокон [1]. Традиционные методы расчета [2-3] прочностных характеристик таких композитов основаны на приближенноаналитических подходах, которые, как правило, приводят к приемлемым результатам при расчете прочности в направлении армирования, но обладают значительной погрешностью при расчете прочностных характеристик в поперечных направлениях и при сдвиге.
Значительный успех в проблеме расчета характеристик композитов был достигнут после создания метода гомогенизации или осреднения и разработки эффективных вычислительных алгоритмов решения локальных задач, возникающих в этом методе [4-7]. В работах [8-10] разработаны вычислительные технологии решения локальных задач на ячайках периодичности, которые возникают в этом методе. Решение данных задач методом конечных элементов является достаточно ресурсоемким и потребовало создания высокопроизводительной распределенной программной системы, которая бы обеспечила информационно-вычислительную поддержку процессов совместного проектирования новых композиционных материалов и конструкций на их основе. Теоретические основы метода гомогенизации применительно к различным микроструктурам упругих композитов изложены в [11]. Целью настоящей работы являлось применение метода микроструктурного конечно-элементного анализа для расчета поверхностей прочности композитов, экспериментальное определение которых, как правило, представляет собой сложную техническую проблему.
Метод расчета поверхностей прочности композитов. Для расчета поверхностей прочности композитов (КМ) предлагается следующий метод, состоящий из 4-х основных этапов.
1. Согласно методу гомогенизации [7, 11] осуществляется решение серии так
называемых локальных задач теории упругости на 1/8 ячейке периодичности
композита V, ,
= °’ в ^ в V' и' и'
,>=2К,» + иАП)!< > , в ^ и г; и г. 0)
ии=ии, е', "„>]к =0 ? е'
где Х£ар - поверхности контакта компонентов композита, X, = {,: , = 0} -
координатные плоскости, X; = {1:,= 1/2}, ^ = 1,2,3 - торцовые поверхности.
ТТа (р\ !-Га са
(рЧ >\Ьг/, у у (рЧ), (рч) - перемещения, напряжения и деформации в компонентах
композита Ка,, а = 1,...,N; с;ы - тензоры модулей упругости компонентов, , -локальные координаты, и.(рд>п - производные по локальным координатам. Функции
и.1р, >(£) удовлетворяют следующим граничным условиям на X :
1. при Р = , ,
и(Р,> = "2 £Р,^Р , рч>/* = 0 , (р,>/I = 0, ^ еХ1-, . Ф к (2)
2. при р ^ Ч,
Щр, >= 4 ^, .,>,.. = 0, и“(р,> = 0, , /,. = {р, ч }; (3)
и>(рч> 4 £рч^'р, и .(р,>/к = 0 , ^(р,> “ 0 , , е Хк , * ^ 7 ^ к ^ * .
где £рч - средние значения деформации композита (заданы). Решение задач (1)-(3) ищется с помощью метода конечного элемента [12].
2. После решения серии задач ^рЧ для всех р, Ч вычисляются эффективные упругие модули композита по следующим формулам:
С т = уг1, ЫР> = Ы р,>) = ЕIУ(р,^, у. = ЕУ(р,> (4)
РЧ к=1 V. Р, Ч=1
3. Далее вычисляется тензор эффективных упругих податливостей ~~1]РЧ \~уРЧ,
и по его компонентам - технические константы: Е = 1/ пйй - модули Юнга,
У'у = —Пг® / Пгш - коэффициенты Пуассона, = с - модули сдвига (здесь по 1, 7 суммирования нет). Затем вычисляем компоненты тензора концентрации напряжений:
Вк (£) = С.РЧ .)прф,,
связывающие напряжения в компонентах композита с
а __ г) —
осредненными напряжениями: = ВуыЫы . Это соотношение удобно переписать в
матричном виде, введя следующие обозначения для групп индексов:
11 ^ 1,22 ^ 2,33 ^ 3,12 ^ 4,13 ^ 5,23 ^ 6:
(.) = В(. )ы7, и 7 = 1.6, (5)
4. Одним из наиболее широко применяемых на практике способов моделирования поверхностей прочности анизотропных материалов является использование квадратичного критерия прочности Малмейстера-Ву. Положим, что разрушение в точке
компонента композита наступает, если в этой точке выполняется критерий Малмейстера -Ву:
/<к ) (£ >=у ) (,}+2 ^ у ) (£ у) . >=1, *, 7=1...6, (6)
где 4к), к : - матрицы и векторы констант прочности -й компоненты композита (заданы),
они имеют по 9 и 3 независимых компонент соответственно.
Подставляя (5) в (6), получим
Г ) (., г • >=$ к ) В . (£ у (/•)+1 ¿у ) Вк (. > в. (. Ык (г ■)ы, (г •)=1
(7)
* ^ *
где I - момент времени разрушения компоненты композита в точке ,
5. Зададим теперь серию из М законов изменения компонент средних напряжений
ык(г) в виде линейной функции от времени нагружения г: (г) = , р = 1...М .
Для каждого в- го элемента серии осуществляется нагружение композита с заданным вектором <Ув} при изменяющемся t от 0 до г*в. Осуществляется итерационное приращение аргумента t и на каждом шаге итерации проверяется выполнение условия /( к) (, ^ >. 1.
Как только при некотором {* впервые выполняется условие /. ) (,*, 1*р > = 1, происходит
разрушение компонента композита. Полагается, что это разрушение и есть разрушение композита в целом. В результате в численном эксперименте находятся значения ,*,^и
ы';т(('в)
5) Будем далее аппроксимировать поверхность прочности композита с помощью того же квадратичного критерия прочности Малмейстера-Ву, но с заранее неизвестными
константами
Осуществляя далее всю серию из М нагружений с различными значениями <с{в} в каждом нагружении, на основе (8) можно составить систему уравнений
+ 2 Я, к™ = 1, в = 1..М (9)
к {в} = с{в} (/*) к {в} = с{в} (/* )с{в} (/*) „
где — с»г в)-> у ~ г V*в' у У в - найденные к данному этапу
величины. Полагая М=12, систему (9) можно рассматривать как систему линейных алгебраических уравнений относительно 12 констант , ¿у . Решая эту СЛАУ, находим константы прочности композита в целом.
Применялся также другой способ нахождения констант прочности , ¿у - из условия
наилучшей аппроксимации экспериментально получаемой поверхности прочности композита в 6мерном пространстве осредненных компонент тензора напряжений ы..
Результаты численного моделирования. Численные расчеты проводились для дисперсно-армированного КМ, структура которого представлена на рисунке 1. Исходные
данные: матрица изотропная (алюминиевый сплав)- Е = 70ГПа, и = 0.35, пределы прочности при растяжении, сжатии и сдвиге сГг2) = 0.2ГПа, <ы(С'> = 0.3ГПа, = 0.1ГПа;
армирующий наполнитель - в виде сферических частиц с константами упругости (карбид кремния) Е = 400ГПа, и = 0.2 и пределами прочности сГгУ) = 0.4ГПа, с(С1) = 0.5ГПа, 0^ = 0.1ГПа. По пределам прочности вычислялись значения констант
Ыа ) Ыа )
прочности ¿у , 3 г .
Были получены следующие значения эффективных упругих констант композита:
Е = 121,61511 ГПа, У= 0,298931, в = 39,85084 ГПа. Число конечных
элементов, на которые разбивалась ЯП, равно 4773. На рис. 1-2 показаны распределения
а т
напряжений су у (рч) в ЯП как результат решения локальных задач ърч .
гг
Рис. 1 - Сдвиговые напряжения ^13(13) , Мпа
Рис. 2 - Нормальные напряжения °33(33) в МПа С помощью изложенного выше алгоритма были получены прочностные
также следующая матрица прочности композита:
Л 790.12346 -1577.29722
¿г.
-1577.29722 -1687.99102 867.20867 -1694.99877
1687.99102
0
0
0
1694.99876
0
0
0
1531.10048
0
0
0
0
0
0
42.10526
0
0
= 15,3ГПа-1
0 0
0 0
0 0
0 0
10526 0
0 40
На рис. 3 показаны сечения поверхности прочности матрицы, армирующего наполнителя и композита. Точками показаны «точные» численные значения прочности, а линиями - аппроксимация с помощью критерия Малмейтера-Ву. Получено достаточно хорошее качество аппроксимации.
с2
\
1—" \
/ \
/ і і
/ / /
і /
-л V /
\
Модули Юнга Е, [МПа]
Коэффициенты Пуассона і'
Е, =3883.327 |/,2 = (>.:« — ГЛОТ,. 28
Ь, = 7 |/ц = 0.33 6'и = 1Ж.Щ
/■..( = .'5882.928 /*23 - 0. -ІЗ с; а = 1306. 32
Модули сдвига С, [МПа]
-0.7-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
Рис. 3 - Сечение поверхностей прочности компонент композита и самого КМ, ГПа
Рис. 4 - Распределение компонент тензоров концентрации напряжений В1212 и В1111 для дисперсно-армированного КМ с учетом поверхности раздела
Были также проведены численные расчеты для случая дисперсно-армированного композита, для которого учитывалась поверхность раздела матрицы и наполнителя. Для этой цели вводился дополнительный 3-й элемент в ЯП. Геометрическая структура такого композита показана на рис. 4, а некоторые результаты расчетов задачи Lpq - на рис. 5.
Рис. 5 - Реальная микроструктура дисперсно-армированного композита на основе эластомерной матрицы и стеклянных микросфер, компьютерная модель микроструктуры и конечно-элементная сетка композита,
использованная в расчетах
Были также проведены расчеты для тканевых композиционных материалов с сатиновым типом переплетения (рис. 6 а, б).
расчетной области для тканевого КМ тензора напряжений для ЯП тканевого КМ с
сатиновым типом переплетения
На рис. 7 показаны результаты численного моделирования задачи о напряженно-деформированном состоянии для ячейки периодичности тканевого КМ с сатиновым типом переплетения.
Рис.7 - Поля функции повреждаемости в при решении задач: сдвига в плоскости XY (слева) и растяжения
по оси Y (справа) для тканевого КМ
Выводы. Разработан математический метод и вычислительная технология расчета прочностных характеристик КМ, в основе которого лежит использование квадратичного критерия прочности Малмейстера-Ву. Показано, что предложенный метод является достаточно эффективным с вычислительной точки: он позволяет прогнозировать поверхность прочности композита, используя только одну итерацию решения локальных задач Lpq. В то же время, следует отметить, что описанный метод нуждается в модифицировании при расчете не начальной поверхности прочности композита, которая фактически и находится этим методом, а предельной поверхности прочности, когда нагружение осуществляется до полного разрушения композита. В этом случае необходимо многократное решение задач Lpq на каждом шаге цикла нагружения, учитывая при этом изменение зоны разрушения компонентов композита.
Список литературы
1. Композиционные материалы: Справочник / Под ред. В.В. Васильева,
Ю.М. Тарнопольского. М.: Машиностроение, 1989. 510 с.
2. Малмейстер А.К., Тамуж В.П., Тетерс Г.А. Сопротивление полимерных и композитных материалов. Рига: Зинатне, 1980. 572 с.
3. Тарнопольский Ю.М., Жигун И.Г., Поляков В.А. Пространственно-армированные композиционные материалы. М.: Машиностроение, 1987. 225 с.
4. Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах.
М.: Наука, 1984. 352 с.
5. Победря Б.Е. Механика композиционных материалов. М.: МГУ, 1984. 336 с.
6. Санчес-Паленсия Э. Неоднородные среды и теория колебаний : пер. с англ. М.:
Мир, 1984. 471 с.
7. Димитриенко Ю.И. Механика композиционных материалов при высоких температурах. М.: Машиностроение, 1997. 356 с.
8. Димитриенко Ю.И., Соколов А.П. Разработка автоматизированной технологии вычисления эффективных упругих характеристик композитов методом асимптотического осреднения // Вестник МГТУ им. Н.Э.Баумана. Сер. Естественные науки. 2008. № 2. С. 57-67.
9. Димитриенко Ю.И., Соколов А.П. Автоматизация прогнозирования свойств композиционных материалов на основе метода асимптотического осреднения // Информационные технологии. 2008. № 8. С. 31-38.
10. Соколов А.П., Шпакова Ю.В., Першин А.Ю. Проектирование и разработка распределенной вычислительной системы инженерного анализа свойств композиционных материалов // Материалы IX Международной конференции по неравновесным процессам в соплах и струях (ЫР^Г2012). Алушта, Украина: Московский Авиационный Институт, 2012. С. 518-521.
11. Димитриенко Ю.И., Соколов А.П. Об упругих свойствах композиционных материалов // Математическое моделирование. 2009. Т. 21, № 4. С. 96-110.
12. Димитриенко Ю.И., Кашкаров А.И. Конечно-элементный метод для вычисления эффективных характеристик пространственно-армированных композитов // Вестник МГТУ им. Н.Э.Баумана. Сер. Естественные науки. 2002. № 2. С. 95-108.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Modeling of surfaces of composite strength on the basis of the microstructural finite element analysis
# 11, November 2012 DOI: 10.7463/1112.0496336
Dimitrienko Yu.I., Sokolov A.P., Shpakova Yu.V., Yurin Yu.V.
Russia, Bauman Moscow State Technical University
dimit.bmstu@gmail.com
alsokolo@bmstu.ru
shpakovayuliya@bmstu.ru
ven8469@ya.ru
Publications with keywords: FEM, asymptotic averaging method, multiscale structure, effective elastic-strength properties, surface of strength, fabric composites Publications with words: FEM, asymptotic averaging method, multiscale structure, effective elastic-strength properties, surface of strength, fabric composites
In this work a method of microstructural finite element analysis for calculation of strength surface of composite materials with periodic fabric structure of reinforcement on the basis of a developed software was proposed. This method is a method of asymptotic averaging adapted for multi-scale periodic structures developed by Prof. Y. Dimitrienko at the "Computational Mathematics and Mathematical Physics" department. The method is based on solving a special class of elasticity problems on 1/8th of the periodicity cell of the composite in order to determine effective elastic properties of the material. The results of calculation of ultimate strength for disperse reinforced and textile composites are presented.
References
1. Vasil'ev V.V., Tarnopol'skii Iu.M., eds. Kompozitsionnye materialy: Spravochnik [Composite materials: Handbook]. Moscow, Mashinostroenie, 1989. 510 p.
2. Malmeister A.K., Tamuzh V.P., Teters G.A. Soprotivlenie polimernykh i kompozitnykh materialov [Resistance of polymeric and composite materials]. Riga, Zinatne, 1980. 572 p.
3. Tarnopol'skii Iu.M., Zhigun I.G., Poliakov V.A. Prostranstvenno-armirovannye kompozitsionnye materialy [Space-reinforced composite materials]. Moscow, Mashinostroenie, 1987. 225 p.
4. Bakhvalov N.S., Panasenko G.P. Osrednenie protsessov v periodicheskikh sredakh [Averaging processes in periodic media]. Moscow, Nauka, 1984. 352 p.
5. Pobedria B.E. Mekhanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, MSU Publ., 1984. 336 p.
6. Sanchez-Palencia E. Non-Homogeneous Media and Vibration Theory. Springer Verlag, Berlin, 1980. (Lecture Notes in Physics, vol. 127.). (Russ. ed.: Sanches-Palensiia E. Neodnorodnye sredy i teoriia kolebanii. Moscow, Mir, 1984. 471 p.).
7. Dimitrienko Iu.I. Mekhanika kompozitsionnykh materialov pri vysokikh temperaturakh [Mechanics of composite materials at high temperatures]. Moscow, Mashinostroenie, 1997. 356 p.
8. Dimitrienko Iu.I., Sokolov A.P. Razrabotka avtomatizirovannoi tekhnologii vychisleniia effektivnykh uprugikh kharakteristik kompozitov metodom asimptoticheskogo osredneniia [Development of automated technology of calculation of effective elastic characteristics of composites by method of asymptotic averaging]. Vestnik MGTU im. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science],
2008, no. 2, pp. 57-67.
9. Dimitrienko Iu.I., Sokolov A.P. Avtomatizatsiia prognozirovaniia svoistv kompozitsionnykh materialov na osnove metoda asimptoticheskogo osredneniia [Automation of predicting the properties of composite materials on basis of the method of asymptotic averaging]. Informatsionnye tekhnologii, 2008, no. 8, pp. 31-38.
10. Sokolov A.P., Shpakova Iu.V., Pershin A.Iu. Proektirovanie i razrabotka raspredelennoi vychislitel'noi sistemy inzhenernogo analiza svoistv kompozitsionnykh materialov [Design and development of distributed computing systems of engineering analysis of the properties of composite materials]. Materialy 9 Mezhdunarodnoi konferentsii po neravnovesnym protsessam v soplakh i struiakh (NPNJ’2012) [Proc. of the 9th International Conference on Nonequilibrium Processes in Nozzles and Jets (NPNJ'2012)]. Alushta, Ukraina, Moskovskii Aviatsionnyi Institut Publ., 2012, pp. 518-521.
11. Dimitrienko Iu.I., Sokolov A.P. Ob uprugikh svoistvakh kompozitsionnykh materialov [On elastic properties of composite materials]. Matematicheskoe modelirovanie,
2009, vol. 21, no. 4, pp. 96-110.
12. Dimitrienko Iu.I., Kashkarov A.I. Konechno-elementnyi metod dlia vychisleniia effektivnykh kharakteristik prostranstvenno-armirovannykh kompozitov [Finite-element method to calculate the effective characteristics of spatially reinforced composites]. Vestnik MGTU im. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2002, no. 2, pp. 95-108.
