Численное моделирование теплового расширения композиционных материалов на основе метода асимптотического осреднения Текст научной статьи по специальности «Физика»

УДК 539.3

Численное моделирование теплового расширения композиционных материалов на основе метода асимптотического осреднения

© Ю.И. Димитриенко, Е.А. Губарева, С.В. Сборщиков МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

Предложен вариант метода асимптотического осреднения композиционных материалов с термоупругими характеристиками, позволяющий вычислять эффективные коэффициенты линейного теплового расширения (КЛТР) композитов. Сформулированы локальные задачи термоупругости на ячейке периодичности (ЯП) композитов. Дана вариационная формулировка задач термоупругости на ЯП, для численного решения которых применен метод конечных элементов. Для программной реализации этого метода использован программный комплекс, разработанный в Научно-образовательном центре «Симплекс» и на кафедре «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Представлены примеры численного решения локальной задачи термоупругости для композиционного материала на основе керамических волокон и полимерной матрицы. Рассчитаны эффективные КЛТР композиционных материалов с пространственным расположением керамических волокон и полимерной матрицей при различных температурах. Показано, что наличие процессов термодеструкции полимерной матрицы при высоких температурах приводит к немонотонной зависимости КЛТР от температуры. Предложенный алгоритм позволяет вычислять КЛТР композитов практически с произвольными структурами армирования волокнами и матрицами, претерпевающими физико-химические превращения при высоких температурах. В отличие от большинства известных приближенных методов расчета КЛТР разработанный метод дает возможность находить точные в математическом смысле значения этих коэффициентов.

Ключевые слова: многослойные тонкие оболочки, метод асимптотического осреднения, асимптотическая теория оболочек.

Введение. Для расчета термонапряжений в конструкциях из композиционных материалов при воздействии повышенных температур кроме сведений об упругих характеристиках необходима информация о коэффициентах линейного теплового расширения (КЛТР). Для расчета КЛТР композитов обычно применяют приближенные методы [1-6], которые обеспечивают приемлемую для инженерных расчетов точность, в основном при простых структурах армирования (однонаправленные и слоистые композиты). Для композитов с более сложными структурами (дисперсно-армированных, тканевых, пространственно-армированных) точность приближенных методов заметно снижается. В настоящее время наиболее точным методом расчета упругих и неупругих характеристик композитов, по-видимому, является ме-

тод гомогенизации (метод асимптотического осреднения) [7-14]. В соответствии с этим методом для расчета эффективных характеристик композитов необходимо решение специальных задач на ячейке периодичности (ЯП) композита. Конечно-элементные методы решения таких локальных задач были разработаны в [10, 11]. С помощью указанного метода были рассчитаны упругие модули пространственно-армированных композитов. Целью настоящей работы является дальнейшее развитие метода гомогенизации и численного алгоритма решения локальных задач термоупругости для расчета КЛТР композитов со сложными структурами армирования.

Метод асимптотического осреднения для задач термоупругости. Рассмотрим пространственно-армированный композиционный материал (рис. 1, 2), которому в пространстве Я3 соответствует область V, имеющая поверхность Е. Композит состоит из N компонентов: компоненты с индексами а = 1, ..., N - 1 представляют собой волокна различных типов, ориентированные по а-му направлению в Я3 , а компонент с индексом а = N — матрица. Обозначим как Га, а = области в пространстве Я3, соответствующие а-му

компоненту композита, а также: Еа — поверхности областей Га; Е^ — поверхности контакта матрицы и волокон (волокна полагаем не контактирующими между собой); Еае — часть поверхности Е композита, занятая а-м компонентом (причем Еа = Е^ и Еае для

N-1

волокон и ЕN = и ЕаМ и Едля матрицы). Волокна и матрицу

а=1

полагаем изотропными термоупругими. Тогда в каждой области Га, а = можно рассмотреть следующую задачу термоупругости:

а,, = 0 в Га;

°а = (гЫ -г<ы) в ГаиЕа;

4=2 (+<1) в Га; (1)

иа = и?; (а,-аЦ = 0 на ЕаN;

„I = и1е на Еае; аап, = на Еа е.

Здесь аа, гО,-, иа — напряжения, деформации и перемещения в а-м

С а *-» а0

уЫ — компоненты тензора модулей упругости; г^ —

тепловые деформации [15], которые, вообще говоря, являются нелинейными функциями от температуры, т. е.

0

s?0 = J (0')d0', (2)

ак/ (б') — компоненты тензора теплового расширения, зависящие от температуры.

Последние два условия в системе (1) — это условия идеального контакта матрицы и волокон, а также условия на границе композита.

Пусть композиционный материал обладает периодической структурой (см. рис. 1, 2), ЯП V которого состоит из N компонентов , а = 1,...,N. Введем малый параметр к = //Ь << 1 как отношение характерного размера ЯП к характерному размеру всего композита, а также глобальные хк и локальные координаты. Будем полагать, что матрица является связной областью. Обозначим также Е^ = Е^ п V поверхности раздела матрицы и волокон в ЯП.

В этом случае для такой структуры может быть применен метод асимптотического осреднения [7-10], согласно которому решение задачи (1) для матрицы и волокон строится в виде асимптотических разложений вида

иа = м(0) (хк ) + киа(1) (хк, I/ ) + к2... ;

еа =еа(0)(хк,^) + кега(1)(хк,^) + к2...; (3)

аа =о;а (0)( хк, ^) + < (1)( хк, ^) + к2...,

причем по аргументу ^ эти функции полагают периодическими. Деформации и напряжения «нулевого уровня» имеют следующий вид:

еа<°>=%+2 («а™+«а1,1'); (4)

_ 1

Sj

аи

sj = 2 («¡3} + «3 >); (5)

^а(0) = j (-s?0), если е е^, а = 1,..., N. (6)

Здесь ,¡ = 8/8x' и и = 8/8 — производные по двум типам координат. При выводе формул (3)-(5) и далее используется правило дифференцирования асимптотических разложений:

df (x, Ъ) 8f (xi, Ъ) | 8f (Xi, Ъ) =8f(XüÍL) + i 8f (X, Ъ)

dxj 8xj 8^ k 9xj 8xj к 8Ъ k

0

Подставляя разложения (2) в систему (1), применяя правило дифференцирования и собирая члены при одинаковых степенях к, получаем так называемую локальную задачу термоупругости на ЯП:

га(0) _

'у/У

_ 0 в V,;

а(0) _ га („а(0)

'у

_ у (гкГв ^^;

еа (0) _„.. +1 2

,,а(1) _ „«(1),

("Й"

+ "У! ) в V, ;

(

(",а (1Л _ 0,

(«> -аN">) _ 0

на Е

,аМ;

У

"

а (1)

_ 0.

Здесь обозначен оператор осреднения:

N N

("а)_Х/„ (.у )

а_1

а_1

(7)

(8)

В(7)условие

_ 0 — это условие периодичности функций на

границе ЯП, а условие 0 вызвано требованием единственно-

сти решения локальной задачи [7]. В силу периодичности функций

а(1)

", имеем следующее соотношение:

*(0)\_ 1/" (0)

„У \„9

_ ^ "9 + ").

(9)

Решение локальных задач термоупругости Ьрч на ЯП. Будем полагать далее, что ЯП является симметричной при преобразованиях

из группы ортотропии относительно осей координат О,5 [16]. Тогда вместо решения локальной задачи (7) на всей ЯП области V, можно перейти к решению задачи на области V,, представляющей 1/8 часть ЯП.

Перемещения «первого уровня», деформации и напряжения «нулевого уровня», а также тепловые деформации представим в виде

следующих сумм:

3

"а(1) _ V "а •

и1 ~ и,(рч)> Р ,Ч _1

а(0) Г .а ; 9 Аи 9 (РЧ)' Р ,Ч _1

_а(0) _ ^ _а ;

аУ _ ^ аУ(рч);

Р ,Ч _1

(10)

а 0

^„а 0 ^У(РЧ)"

РЧ _1

; ,к £ , а _ 1,..., N

где обозначено:

еа0 = Л еа0

ьУ (РЧ) ^УРЯ^рд ,

(11)

а Луп = 1/2 (бгр5уд +5^5ур) — единичный тензор 4-го ранга [16] (суммирования пор и д в этих формулах нет).

Для функций «Орд) для каждой комбинации (рд) может быть выделена линейная часть по локальным координатам:

«аРЧ) = _еРЧ (5'р£д + 5гд£р) + ^рд)(^/), (12)

где 5гр — символ Кронекера; Ц(рд) ) — новые неизвестные

функции, для которых при каждом фиксированном наборе индексов (рд) получаем локальные задачи термоупругости на ЯП:

а У( рд)/У = 0 в ^;

а _ г^а ¡ а а0 \ тт у г у .

аУ(рд) = СУк/ (еу(рд) " еу(рд)) в ^ ;

'у(рд)

2 (и

а тто.

(рд)/У + иУ (рд)/'

■)

в V;

(13)

и' (рд) = и'

■(рд)

.о.

(рд)'

N

(ау(рд) _ аУ(рд)) пУ = 0

на Е

^аЛ.

Здесь обозначено:

ек/ (рд) = 1 е рд(5 кр5 /д +5 А).

(14)

Кроме того, к системе (13) присоединяются условия на координатных плоскостях Е5 = = 0} и на торцевых поверхностях ЯП Е^ = {£, ^ = 1/2}, ^ = 1, 2, 3, которые записываем следующим образом [10]:

при (р = д) и ' Ф у Ф к Ф ■

и1(рд) = 1/ 2ерд5'р; $°(рд) = 0; рд) = 0 на Е'';

при (р Ф д) и и у ={p, д}

пс -

■ (рд) 1>~1р"гр при (р Ф д) и ■ Ф у Ф к Ф'

<(рд) = (1/ 4)егр5гр; $°(рд) = 0; ик(рд) = 0 на Еу

рд) = 0; $°(рд) = 0; ик°( рд) = 0 на Ек

а также

при (р = д) и ■ Ф у Ф к Ф'

^'^Срд) = 0; рд) = 0; рд) = 0 на Е ■

(15)

при (Р Ф Ч) и К 9 _ {Р, Ч}

и

при (р ф ч) и , Ф 9 Ф к Ф,

а(РЧ) _ 0; РЧ) _ 0; ик(РЧ) _ 0 на Е9;

(15)

^рч) _ 0; 89(РЧ) _ 0; РЧ) _ 0 на Е к .

Здесь введены обозначения для векторов усилий: 8,( РЧ) = ^^аа^а})".

Решение задачи (13), (15) разыскиваем в области V,, представляющей 1/8 часть ЯП: V, _ V, п (,, ^ 0), здесь также как Е,ар обозначена поверхность контакта компонентов внутри V, : Е,ар _Е,ар п V,. Представим это решение в виде суммы решений двух задач:

тта _ тта{1} , тта{2} иг(РЧ) - и,(РЧ)+ и,(РЧ).

(16)

Первая — это задача (13), (15) с нулевыми тепловыми деформаци-

ями вОдрч) _ 0 (задача Ьрч):

7а{1} _

9 (РЧ)/9

_ 0 в V,;

9}д) _ 9 еа {РЧ) в ^иЕ, иЕ,; (РЧ) _ 2((Ч)/9 + и9срЧ)/,) в ^;

°9 (РЧ)

тта{1} _ т^{1} (рч) - (РЧ),

(а 9'

а{1} 9 (РЧ)

-а9(Рл) _ 0

^ {1} 9 ( рч) ) "9

на Е

,аN.

(17)

Здесь

при (р = ч) и , Ф 9 Ф к Ф г

иурч) _1/ 2ерч^,р; 89(рч) _ 0; "ч) _ 0 на Е>'

,(РЧ) ""РЧ^У

при (Р Ф Ч) и и 9 _{p, ч}

и

Ч (РЧ)

"РЧ) _ (1/ 4) .гр^ ,р• 89( рч) _ 0; ик( РЧ) _ 0 на Е'

а также

при (р _ ч) и г ф 9 ф к ф г

и?( рч) _ 0; 8я рч) _ 0; 8к(рч) _ 0 на е''

,(РЧ) ' °9'( РЧ)

при (р ф ч) и г ф 9 ф к ф г

(18)

(19)

8К рч) _ 0; рч) _ 0; ика( рч) _ 0 на Ек

при (р Ф д) и у ={p, д}

при (р Ф д) и ■ Ф у Ф к Ф'

*(рд) = 0; рд) = 0; и*(рд) = 0 на Еу;

(19)

рд) = 0; рд) = 0; иО( рд) = 0 на Е к

Вторая задача — с нулевыми средними деформациями ек/(рд) = 0

(задача Ь„д)

а

а{2}

■У(рд

а{2} = с а (еа{2}

аУ ( рч)/У = 0 в ;

урч) = С°к/ ^(м) "еа/(0рд)) в ^ ;

еЙ?д) = 2 ((Крд)/У + идргУ» ) в ;

тт а{2} = ПЛ {2}

■ (рд) ' (рд)'

(а?{2>

-аЛ({22)) = 0

■У (рд) °У( рд )ГУ

на ЕЕ

4аЛ.

Здесь

при (р = д) и ■ Ф у Ф к Ф'

рд) = 0; $у( рд) = 0; рд) = 0 на Е

при (р Ф д) и ^у ={p, д}

и?( рд) = 0; $у( рд) = 0; ик( рд) = 0 на Е';

при (р Ф д) и ■ Ф у Ф к Ф'

рд) = 0; ^Я рд) = 0; ик( рд) = 0 на Е к

а также:

при (р = д) и ■ Ф у Ф к Ф'

рд) = 0; ^Я рд) = 0; рд) = 0 на Е ■

при (р Ф д) и у ={p, д}

рд) = 0; $°( рд) = 0; рд) = 0 на Е j

при (р Ф д) и ■ Ф у Ф к Ф'

(20)

(21)

(22)

рд) = 0; $°(рд) = 0; и%( рд) = 0 на Е к.

Расчет эффективных модулей упругости и КЛТР композита.

После решения двух серий задач Ьрд и Ь°,д для всех комбинаций (рд)

проинтегрируем напряжения по областям, занятым компонентами композита. Тогда эффективные упругие модули композита могут быть вычислены по следующим формулам:

^{1} -

9 (РЧ)

N

_ а

9(рч) /

а_1 К„

(23)

С- _

^11ПП

9( РЧ)

"9РЧ „

РЧ

После расчета тензора модулей упругости Сурч вычисляем тен-

\-1

зор эффективных упругих податливостей П,рч _ (С дрч) , а по его компонентам находим технические константы: Ег _ 1/ Пт — модули Юнга; У9 = —П 9 / Пш — коэффициенты Пуассона; О" _ Сщ — модули сдвига (здесь по г, 9 суммирования нет). Затем с помощью решения второй серии задач 1?рч для тепловых напряжений находим эффективные тепловые деформации композита „к:

а{2} _

а.

а{2} .

Р,Ч _1

9(рч)

(24)

„» _ П а{2}

„к/ _ - Пк/г9аг9 .

Вариационная формулировка локальной задачи Ьп. Вариационная формулировка задачи Ьрч была получена в [10]. Вариационная

формулировка задачи Ь°,ч при фиксированных значениях р и ч для произвольного конечного объема V ^ V, имеет вид

|5.таdV _|5ТтSdЕ-{5.т(0dV. (25)

V Е V

Здесь обозначены координатные столбцы: псевдоперемещений т , напряжений а, деформаций . и поверхностных усилий 8:

и _

о _

е _

8 _

-0

о _

и1 РЧ), и2(РЧ), и3(РЧ)

а _а _а _а / Г^ _а

I / >/2

о11(pч), о22(pч), о33(рч^ о13(рч)■

°23( рч)

/л/2,

012( РЧ)

/72

Г^Ч), „22(рч), е33(рч), 813(РЧ.23(РЧ.а2(РЧ)^2! ;

11( РЧ)

П т

па па па I

_S1(рч), РЧ), S3(рч)! ;

(26)

>а0

■а 0

>а0

>а0

■а 0

■а 0

о11(рч^ о22(pч), о33(pч), о13(pч), о23(pч), о12(рч)

(а 0 _ с а „а 0 а9(РЧ) _ С9к/ .к/(рч).

Определяющие соотношения, следующие из второй группы уравнений в (13), с использованием координатных столбцов записываем таким образом: а = Cs, где С — матрица упругости размером

[6x6], составленная из компонентов тензора С^ы стандартным образом [9]. Соотношения Коши (третья группа уравнений в (13)) в матричном виде записываем так: s = DU, где D — матрица линейных

дифференциальных операторов дифференцирования (9' =9/ ) [10]. Тогда вариационное уравнение (25) можно представить в виде

J (DSU)T CDUdV = J bUTSdE-J(DSU)Tё0dV. (27)

Результаты численного моделирования. Решение вариационного уравнения (27) осуществляли методом конечных элементов с помощью программного обеспечения, разработанного в Научно-образовательном центре «Симплекс» и на кафедре «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана.

В качестве примера был рассмотрен композиционный материал на основе керамических волокон и полимерной матрицы. Волокна в ЯП были ориентированы по четырем главным диагоналям куба и соприкасались в нескольких точках. Термоупругие свойства волокон и матрицы предполагали изменяющимися при нагреве до высоких температур. При расчетах были использованы значения констант термоупругости и КЛТР и волокон, матрицы приведенные в табл. 1 и 2.

Таблица 1

Значения констант термоупругости волокон и матрицы, принятые в расчетах

Температура, Модуль Коэффициент Модуль упругости Коэффициент

оС упругости Пуассона матрицы, Пуассона

волокон, ГПа волокон ГПа матрицы

20 122 0,17 3 0,37

300 120 0,17 0,5 0,37

1000 100 0,17 0,1 0,37

2000 80 0,17 0,1 0,37

2500 40 0,17 0,1 0,32

Таблица 2

Значения КЛТР матрицы и волокон, принятые в расчетах, а также значения КЛТР для композиционного материала

Температура, оС КЛТР матрицы, КЛТР волокон, КЛТР композита,

К-1 • 10-6 К-1 • 10-6 К-1 • 10-6

20 20 5 11

300 1 5 4

600 -5 5 4,5

2000 4 5,5 5,2

2500 4,4 5,6 5,3

В силу линейности задачи термоупругости (20)-(22) задачи ¿°1, 1?22 , ¿3 могут быть решены совместно как одна задача ¿0. Матрица и волокна предполагаются изотропными, поэтому .Ц0рч) _ 0 при

^ 7-0

р Ф ч, в результате решения задач Ьп являются нулевыми, т. е.

аа^Ч) _ 0 и иОрЧ) = 0 при р Ф ч. В связи с этим для расчета КЛТР

композита достаточно решить одну задачу ¿0. Некоторые результаты

решения локальной задачи ¿1 для композиционного материала на основе керамических волокон и полимерной матрицы показаны на рис. 3, 4. Было установлено, что в зонах соприкосновения волокон реализуются наибольшие значения нормальных и касательных напряжений.

При нагреве полимерной матрицы значения КЛТР существенно изменяются. При температурах 300...800 оС в матрице происходят процессы термодеструкции, сопровождающиеся усадкой [17, 18], и значение КЛТР становится отрицательным. Были проведены расчеты КЛТР композиционного материала при различных температурах: 20, 300, 600, 2000 и 2500 оС. Результаты расчетов КЛТР композита в зависимости от температуры представлены на рис. 5 и в табл. 2.

Рис. 1. Компьютерная модель и КЭ-сетка композита с ориентацией волокон 4Б (показаны только волокна)

ВяааааДй

шшшшит

ШШщщШ Щ

Рис. 2. Компьютерная модель и КЭ-сетка матрицы композита с ориентацией волокон 4Б (волокна не показаны)

Рис. 3. Поле компонента напряжения о{1} в ЯП композита

{2}

Рис. 4. Поле компонента напряжения о^ в ЯП композита

Рис. 5. График изменения расчетных значений КЛТР композиционного материала вдоль направления волокон в зависимости от температуры нагрева

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

Выводы. Предложен вариант метода асимптотического осреднения композиционных материалов с термоупругими характеристиками, позволяющий вычислять эффективные КЛТР композитов. Сформулированы локальные задачи термоупругости на ЯП композитов. Дана вариационная формулировка задач термоупругости на ЯП, для численного решения которых применен метод конечных элементов.

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

Показано, что наличие процессов термодеструкции полимерной матрицы при высоких температурах приводит к немонотонной зависимости КЛТР от температуры.

ЛИТЕРАТУРА

[1] Кристенсен Р. Введение в механику композитов. Москва, Мир, 1982, 336 с.

[2] Тарнопольский Ю.М., Жигун И.Г., Поляков В.А. Пространственно-армированные композиционные материалы. Москва, Машиностроение, 1987, 223 с.

[3] Зарубин В.С., Кувыркин Г.Н., Савельева И.Ю. Оценка методом самосогласования температурного коэффициента линейного расширения с дисперсными включениями. Наука и образование, 2015, № 2, с. 197-215.

[4] Ran Z., Yan Y., Li J., Qi Z., Yang L. Determination of thermal expansion coefficients for unidirectional fiber-reinforced composites. Chinese Journal of Aeronautics, October 2014, vol. 27, is. 5, рр. 1180-1187.

[5] Karadeniz Z.H., Kumlutas D. A numerical study on the coefficients of thermal expansion of fiber reinforced composite materials. Composite Structures, March 2007, vol. 78, is. 1, рр. 1-10.

[6] Rupnowskia P., Gentza M., Sutterb J.K., Kumosaa M. An evaluation of the elastic properties and thermal expansion coefficients of medium and high modulus graphite fibers. Composites. Part A: Applied Science and Manufacturing, March 2005, vol. 36, is. 3, рр. 327-338.

[7] Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. Москва, Наука, 1984, 352 с.

[8] Санчес-Паленсия Э. Неоднородные среды и теория колебаний. Москва, Мир, 1984, 472 с.

[9] Победря Б.Е. Механика композиционных материалов. Москва, Изд-во МГУ, 1984, 243 с.

[10] Димитриенко Ю.И., Кашкаров А.И. Конечно-элементный метод для вычисления эффективных характеристик пространственно-армированных композитов. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2002, № 2, с. 95-108.

[11] Димитриенко Ю.И., Сборщиков С.В., Беленовская Ю.В., Анискович В. А., Перевислов С. Н. Моделирование микроструктурного разрушения и прочности керамических композитов на основе реакционно-связанного SiC. Наука и образование, 2013, № 11. doi: 10.7463/1113.0659438 (дата обращения 06.11.2015).

[12] Димитриенко Ю.И., Яковлев Н.О., Ерасов В.С., Федонюк Н.Н., Сборщиков С.В., Губарева Е.А., Крылов В. Д., Григорьев М.М., Прозоровский А.А. Разработка многослойного полимерного композиционного материала с дискретным конструктивно-ортотропным заполнителем. Композиты и наноструктуры, 2014, т. 6, № 1, с. 32-48.

[13] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В., Федонюк Н.Н. Моделирование вязкоупругих характеристик слоисто-волокнистых полимерных композиционных материалов. Наука и образование, 2014, № 11. doi: 10.7463/1114.0734246 (дата обращения 03.11. 2015).

[14] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой. Математическое моделирование и численные методы, 2014, № 1, с. 36-57.

[15] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4: Основы механики твердого тела. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2013, 624 с.

[16] Димитриенко Ю.И. Тензорное исчисление. Москва, Высшая школа, 2001, 576 с.

[17] Dimitrienko Yu.I. A Structural Thermo-Mechanical Model of Textile Composite Materials at High Temperatures. Composites Science and Technology, 1999, vol. 59, no. 7, pp. 1041-1053.

[18] Dimitrienko Yu.I. Thermal Stresses in Ablative Composite Thin-Walled Structures Under Intensive Heat Flows. International Journal of Engineering Science, 1997, vol. 35, no. 1, pp. 15-31.

Статья поступила в редакцию 01.07.2015

Ссылку на эту статью просим оформлять следующим образом: Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Численное моделирование теплового расширения композиционных материалов на основе метода асимптотического осреднения. Инженерный журнал: наука и инновации, 2015, вып. 12. URL: http ://engjournal.ru/catalog/msm/pmcm/1452. html

Димитриенко Юрий Иванович — д-р физ.-мат. наук, заведующий кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана, директор Научно-образовательного центра «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана, действительный член Академии инженерных наук. Автор более 300 научных работ в области механики сплошной среды, вычислительной механики, газодинамики, механики и термомеханики композитов, математического моделирования в науке о материалах. e-mail: dimit.bmstu@gmail.com

Губарева Елена Александровна — канд. физ.-мат. наук, доцент, заместитель заведующего кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 40 научных работ в области механики сплошных сред, механики контактного взаимодействия, математического моделирования, механики композитов. e-mail: gubareva_ea@pochta.ru

Сборщиков Сергей Васильевич — аспирант кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор 20 научных работ в области численного моделирования в механике композитов. e-mail: servasbor@gmail.com

Numerical simulation of composite material thermal expansion by homogenization method

© Yu.I. Dimitrienko, E.A. Gubareva, S.V. Sborschikov Bauman Moscow State Technical University, Moscow, 105005, Russia

The article considers a variant of the asymptotic homogenization method for calculation of effective thermal expansion coefficients of composite materials with thermoelastic properties. We formulate problems of local thermoelasticity over a periodicity cell of composites. A variational formulation of the thermoelasticity problem over a periodicity cell is proposed. A finite element method for computational solving of these problems of thermoelasticity is applied. For software implementation of the finite element method we use the software package developed by the Scientific and Educational Center of the BMSTU. We also give examples of numerical solution of the local problems of thermoelasticity for composites based on ceramic fibers and the polymer matrix. Effective coefficients of thermal expansion for composite materials with spatial arrangement of ceramic fibers and a polymer matrix were calculated for different temperatures. We show that processes of thermal decomposition of polymer matrix result in nonmonotonic dependence of the thermal expansion coefficient on temperature. The proposed algorithm allows to calculate the thermal expansion coefficients for composites with almost any structures of fiber reinforced matrices undergoing physicochemical transformations at high temperatures. Unlike a large number of the well-known approximate methods for calculating thermal expansion coefficients the proposed method allows to obtain the mathematically accurate values for these coefficients.

Ключевые слова: multilayer thin shell, asymptotic homogenization method, asymptotic theory of shells.

REFERENCES

[1] Christensen R.M. Mechanics of Composite Materials. Wiley-Interscience, New York, 1979, 348 p. [In Russian: Christensen R. M. Vvedenie v mekhaniku kompozitov. Moscow, Mir Publ., 1982, 336 p.].

[2] Tarnopolskiy Yu.M., Zhigun I.G., Polyakov V.A. Prostranstvenno armirovannye kompozitsionnye materialy [Spatially Reinforced Composite Materials]. Moscow, Mashinostroenie Publ., 1987, 224 p.

[3] Zarubin V.S., Kuvyrkin G.N., Savelyeva I.Y. Nauka i obrazovanie: electronnyy nauchno-tekhnicheskiy zhurnal - Science and Education: Electronic Scientific and Technical Journal, 2015, no. 2, pp. 197-215.

[4] Zhiguo Ran, Ying Yan, Jianfeng Li, Zhongxing Qi, Lei Yang. Chinese Journal of Aeronautics, October 2014, vol. 27, iss. 5, pp. 1180-1187.

[5] Z. Haktan Karadeniz, Dilek Kumlutas. Composite Structures, March 2007, vol. 78, iss. 1, pp. 1-10.

[6] Rupnowskia P., Gentza M., Sutterb J.K., Kumosaa M. Composites Part A: Applied Science and Manufacturing, March 2005, vol. 36, iss. 3, pp. 327-338.

[7] Bakhvalov N.S., Panasenko G.P. Homogenisation: Averaging Processes in Periodic Media. Mathematical Problems in the Mechanics of Composite Materials. Springer Publ., 1989, 352 p. [In Russian: Bakhvalov N.S., Panasenko G.P. Osrednenie protsessov v periodicheskikh sredakh. Moscow, Nauka Publ., 1984, 352 p.].

ffuMumpueHKO, E.A. ryGapeBa, C.B. C6OP^UKOB

[8] Sanchez-Palencia E. Non-homogeneous media and vibration theory. Lecture Notes in Physics. Berlin, Springer Publ., 1980, vol. 127, 398 p. [In Russian: Neodnorodnye sredy i teoriya kolebaniy. Moscow, Mir Publ., 1984].

[9] Pobedrya B.E. Mekhanika kompozitsionnykh materialov [Mechanics of Composite Materials]. Moscow, MGU Publ., 1984, 352 p.

[10] Dimitrienko Yu.I., Kashkarov A.I. Vestnic MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2002, no. 2, pp. 95-108.

[11] Dimitrienko Yu.I., Sborshchikov S.V., Belenovskaya Yu.V., Aniskovich V.A., Perevislov S.N. Nauka i obrazovanie: electronnyy nauchno-tekhnicheskiy zhurnal - Science and Education: Electronic Scientific and Technical Journal, 2013, no. 11. doi: 10.7463/1113.0659438

[12] Dimitrienko Yu. I, Yakovlev N.O., Erasov V.S., Fedonuk N.N., Sborshchikov S.V., Gubareva E.A., Krylov V.D., Grigoryev M.M., Prozorovskiy A.A. Kompozity i nanostruktury — Composites and Nanostructures, 2014, vol. 6, no. 1, pp. 32-48.

[13] Dimitrienko Yu.I., Gubareva E.A., Sborshchikov S.V., Fedonuk N.N. Modeling of viscoelastic properties of laminated fibrous polymer composite materials. Science and Education: Electronic Scientific and Technical Journal, 2014, no. 11. doi: 10.7463/1113.0659438

[14] Dimitrienko Yu.I., Gubareva E.A., Sborschikov S.V. Matematicheskoe modelirovanie i chislennye menody - Mathematical modeling and Numerical Methods, 2014, no. 1, pp. 36-57.

[15] Dimitrienko Yu.I. Mekhanika sploshnoy sredy. Vol. 4. Osnovy mekhaniki tverdogo tela [Continuum Mechanics. Fundamentals of Solid Mechanics]. Moscow, BMSTU Publ., 2013, vol. 4, 624 p.

[16] Dimitrienko Yu.I. Tenzornoye Ischislenie [Calculus of Tensors]. Moscow, Vysshaya shkola Publ., 2001, 576 p.

[17] Dimitrienko Yu.I. Composites Science and Technology, 1999, vol. 59, no. 7, pp. 1041-1053.

[18] Dimitrienko Yu.I. International Journal of Engineering Science, 1997, vol. 35, no. 1, pp. 15-31.

Dimitrienko Yu. I. (b.1962) graduated from Lomonosov Moscow State University in 1984. Dr. Sci.(Phys. & Math.), Professor, Head of the Department of Computational Mathematics and Mathematical Physics, Director of Scientific-Educational Center of Supercomputer Engineering Modeling and Software Development, Bauman Moscow State Technical University. Member of the Russian Academy of Engineering Science. Author of over 300 publications in the field of computational mechanics, gasdynamics, mechanics and thermomechanics of composite materials, mathematical simulations in material science. e-mail: dimit.bmstu@gmail.com

Gubareva E.A. (b.1982) graduated from Lomonosov Moscow State University in 2004. Cand. Sci. (Phys. & Math.), Associated Professor, the Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University. Author of over 30 publications in the field of continuum mechanics, tribology, mathematical modeling and composite mechanics. e-mail: gubareva_ea@pochta.ru

Sborschikov S.V. (b.1989) graduated from Bauman Moscow State Technical University in 2012. Postgraduate student at the Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University. The author of 20 publications in the field of numerical simulation in mechanics of composites. e-mail: servasbor@gmail.com