CC BY
64
Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2015. №2. С. 197-215.
Б01: 10.7463/0215.0757059
Представлена в редакцию: 09.02.2015 Исправлена: 21.02.2015
© МГТУ им. Н.Э. Баумана ХДК 539.3
Оценка методом самосогласования
температурного коэффициента линейного расширения
композита с дисперсными включениями
Зарубин В. С.1'*, Кувыркин Г.Н.1, Савельева И.Ю.1 *fii2@bmstu.ru
1МГТУ им. Н.Э. Баумана, Москва, Россия
С применением метода самосогласования получены расчетные зависимости, устанавливающие связь температурного коэффициента линейного расширения композита, армированного изотропными включениями шаровой формы, с термоупругими характеристиками матрицы и включений и их объемной концентрацией. Проведено сравнение этих зависимостей с известными формулами, полученными с использованием других подходов, и с двусторонними оценками возможных значений этого коэффициента, которые следуют из вариационной формулировки задачи термоупругости. Представленные расчетные зависимости позволяют прогнозировать значение температурного коэффициента линейного расширения композита с дисперсными изотропными включениями, в том числе с наноструктурными элементами.
Ключевые слова: композит; метод самосогласования; двусторонние оценки; изотропные шаровые включения; температурный коэффициент линейного расширения
Наука и Образование
МГТУ им. Н.Э. Баумана
Сетевое научное издание
ISSN 1994-0408
Введение
Композиты в качестве конструкционных и функциональных материалов находят широкое применение во многих отраслях техники. Среди механических свойств, определяющих конкретные области применения композитов при наличии тепловых воздействий, важное место занимают термоупругие характеристики, в том числе температурный коэффициент линейного расширения. Оценку этого коэффициента по известным свойствам матрицы композита и армирующих включений можно получить, используя различные подходы [1, 2, 3, 4, 5, 6, 7]. Один из таких подходов, называемый методом самосогласования [8, 9, 10, 11, 12], состоит в учете взаимодействия отдельно взятого элемента неоднородности материала (включения или частицы матрицы композита) с однородной изотропной средой, характеристики которой подлежат определению [13]. Последующее осреднение по представительному объему материала возмущений распределения параметров в этих элементах позволяет получить расчетные зависимости для оценки искомых характеристик рассматриваемого композита.
В технике наряду с использованием композитов, армированных волокнистыми включениями, широко применяют дисперсно упрочненные композиты, в которых включениями являются частицы из высокопрочных и высокомодульных материалов, в том числе нано-структурные элементы [14, 15, 16, 17]. Обычно дисперсные частицы имеют близкие размеры во всех направлениях, что позволяет форму таких частиц в первом приближении считать шаровой [18, 19, 20]. Достаточно широкое распространение получили композиты, называемые сферопластиками [21], которые содержат дисперсные частицы строго шаровой формы, например полые стеклосферы [22, 23, 24].
Близкими по тематике являются работы [25, 26, 27, 28, 29, 30, 31, 32], посвященные термоупругим свойствам пористых металлов, металлокерамики, полимеров и строительных материалов, так как отдельно взятую пору допустимо рассматривать как абсолютно податливое включение с нулевым температурным коэффициентом линейного расширения. Достоверные оценки теплового расширения композита особенно важны при его использовании в технических устройствах, к которым предъявляются высокие требования по размерной стабильности [33] и в которых используемый материал включает разнородные компоненты [34].
Значительное число публикуемых работ по оценке термоупругих характеристик композитов с дисперсными включениями подтверждает актуальность исследований в этом направлении. Такие исследования требуют комплексного подхода и использования методов компьютерного моделирования, поскольку построение математической модели, описывающей тепловое расширение композита, неразрывно связано с анализом его напряженно-деформированного состояния [35, 36, 37]. В данной работе применение метода самосогласования позволило реализовать такой подход и получить расчетные зависимости, связывающие температурный коэффициент линейного расширения дисперсно упрочненного композита с термоупругими характеристиками матрицы и включений и их объемной концентрацией. Для проверки достоверности этих зависимостей построены двусторонние границы области возможных значений этого коэффициента, основанные на ранее разработанной двойственной вариационной формулировке задачи термоупругости неоднородного твердого тела.
1. Основные соотношения
В случае шаровых включений одинакового размера при их наиболее плотной упаковке, когда в каждом слое центры контактирующих между собой шаров находятся в вершинах правильных треугольников, объемная концентрация Су таких включений принимает наиболь-
п
шее значение Су = —^ ~ 0,7405 [38]. Если не учитывать конкретный размер включений
3\/ 2
и считать, что их диаметр может принимать значения от некоторого конечного до бесконечно малого, т. е. возможно заполнение всех пустот между крупными включениями, то допустимо рассматривать значения Су € [0, 1].
В более общем случае материал включений может быть различным. Примем, что композит армирован включениями из N материалов, причем
N
ЕС = СУ е [0, 1], (1)
где С? — объемная концентрация включений из материала с номером я. Для каждого материала считаем известными объемный модуль упругости , модуль сдвига Оя и температурный коэффициент линейного расширения ая, а также соответствующие параметры Кт, От и ат для материала матрицы композита. Расположение включений примем случайным, что позволяет считать композит макроскопически изотропным с объемным модулем упругости К и модулем сдвига О. Через эти модули можно выразить продольный модуль упругости (модуль Юнга) Е = ^^ и коэффициент Пуассона V 1 3С
3 К + С ^ ^ 2 6 К + 2 С
композита [39]. Искомый температурный коэффициент линейного расширения композита
обозначим а.
Упругие характеристики включений, матрицы и композита в целом можно представить в одинаковой тензорной форме [ 1, 40] (на примере характеристик композита)
С = 3К V + 2ОБ, (2)
где С — тензор четвертого ранга коэффициентов упругости, а V и б — тензоры четвертого ранга, являющиеся соответственно объемной и девиаторной составляющими единичного тензора четвертого ранга I = V + б. Если использовать единичный тензор второго ранга с компонентами 5^, г, ] = 1, 2, 3 (5^ = 1 при г = ] и 5^ = 0 при г = ]), то компоненты этих составляющих, определенные в выбранной прямоугольной декартовой системе координат Ох1х2х3, примут вид
= 3^^ = ^^ + ^^ - т,п = 1 2 3.
С учетом этих формул, правила суммирования слагаемых по повторяющимся в сомножителях латинским индексам [41] и равенства 5ц = 3 можно установить, что У^ы 5к1 = 5^,
Б^ы5к1 = 0, к, I = 1, 2, 3, а также [1]
У^тпУтпк1 ^щк^ Б1]тпБтик1 Dijkl, У^тпБтпк1 ° (3)
Тензоры V и б и представленный через них в формуле (2) тензор С являются изотропными, т. е. их компоненты не зависят от ориентации осей выбранной системы прямоугольных декартовых координат.
2. Метод самосогласования
В процедуре применения метода самосогласования для оценки температурного коэффициента линейного расширения композита можно выделить ряд последовательных этапов. В качестве представительного объема рассматриваемого композита, осредненные параметры
которого отражают свойства этого композита в целом, примем неограниченную область, заполненную линейно упругой средой. В этой среде при исходной однородной по этому объему температуре Т0 первоначально отсутствуют напряжения и связанные с ними деформации.
На первом этапе рассмотрим упругое взаимодействие шарового включения с занимающей неограниченную область окружающей средой из того же линейно упругого изотропного материала с тензором С коэффициентов упругости. Если мысленно извлечь включение из окружающей среды и создать в нем однородную деформацию, определяемую тензором ё", то после возвращения включения на прежнее место деформированное состояние этого включения также будет однородным [13], определяемым тензором ё^. Этот тензор связан с тензором ё" соотношением [1, 13] ё" = Л¥ •• ё% где
* =3^ V 6. (4)
1 + V 2 4 — 5и
При этом напряженное состояние включения будет характеризовать тензор
§• = С •• (ё^ — ё") = С •• ё^ — С •• * •• ё\ (5)
Здесь каждая из точек между сомножителями в произведении тензоров означает свертывание по индексу, одинаковому в обоих сомножителях [41]. В соотношение (5) входит внутреннее произведение изотропных тензоров С и АУ. С учетом формул (2)-(4) и представления коэффициента Пуассона через модули К и С получим
5 3 К + 4С
С •• * = (3К + 4С^ + 5 зк+24с6. (6)
На втором этапе повысим температуру представительного объема на величину ДТ, что приведет к возникновения однородной по этому объему температурной деформации, определяемой тензором ё(т) с компонентами аДТ^-. Поскольку в представительном объеме на данном этапе находится однородный изотропный линейно упругий материал, изменение температуры этого материала не повлияет на напряженное состояние включения, но его полную деформацию теперь будет определять сумма тензоров ё* = ё^ + ё(т).
На третьем этапе при однородном распределении температуры Т0 + ДТ в представительном объеме и условии сохранения полной деформации ё* проведем замену исходного
включения анизотропным шаровым включением из материала с номером с = 1, N. При этом в силу изменения термоупругих свойств нового включения произойдет изменение и его напряженно-деформированного состояния по сравнению с исходным включением. Теперь деформированное состояние нового включения будет определять тензор ё? = ё* — ё(т) = = ё^ + ё(т) — ё(т), где ё(т) — тензор температурной деформации этого включения с компонентами ая ДТ5^, а его однородное напряженное состояние будет характеризовать тензор
в? = С •• ё? = С(ё^ + ё(т) — ё(т)). (7)
Эквивалентность замены включений требует не только совпадения их полной деформации, но и совпадения напряженных состояний, т. е. равенства £• = . Исключая с учетом этого равенства из соотношений (5) и(7) тензор ё^, получаем
8? = (С, - С + С •• Л¥)-1 •• (С - С •• Л¥) •• С? •• (ё(т) - ё(ГГ)). (8)
Форма частиц матрицы, связывающей в композите шаровые включения, при хаотическом расположении включений является неопределенной. Примем в качестве средней статистической форму этих частиц шаровой. Взаимодействие шаровой частицы матрицы с изотропной окружающей средой будут определять соотношения, аналогичные полученным выше для шарового включения, а именно напряженное состояние в частице матрицы будет описывать тензор
8т = (Ст - С + С ■■ Л¥)-1 ■■ (С - С ■■ Л¥) ■■ Ст ■■ (ё^Г} - ё(Т)), (9)
где Ст — тензор коэффициентов упругости материала матрицы, а ёт — тензор ее температурной деформации с компонентами атДТ#у.
Заключительным этапом применения метода самосогласования для получения расчетных зависимостей, определяющих искомый температурный коэффициент линейного расширения а рассматриваемого композита, является осреднение по представительному объему возмущений напряженного состояния во всех включениях и частицах матрицы. Для тензора напряжений как тензора второго ранга такое осреднение эквивалентно вычислению его первого инварианта, т. е. умножению этого тензора на единичный тензор второго ранга 12 [41, 42]. Тогда для включений из материала с номером с, учитывая формулы (2), (3), (6) и (8), получаем
Л ? 12СКЯа - а) 8? "12 = 3К + 4С . Аналогично, используя формулы (2), (3), (6) и (9), находим
Л I = 12СКт(ат - а) 8т " 12 = 3Кт + 4С . Итоговый результат осреднения по всем включениям и частицам матрицы в соответствии с методом самосогласования [10] должен быть равен нулю, т. е. с учетом двух последних равенств и формулы (1)
_ N _
12СК"' зга(1 - > + 12С£ К = 0
Разрешив это равенство относительно искомого температурного коэффициента линейного расширения рассматриваемого композита, запишем
= (К 1 - Су + Л к,с у1/ К 1 - + * аКСЛ (10)
а = ^Кт3Кт + 4С + ^ 3К, + 4С) ^атКт3Кт + 4С + ^ 3К? + 4С). (10)
Следует отметить, что к такому же результату приводит осреднение по представительному объему тензоров деформации во всех включениях и частицах матрицы композита.
3. Двусторонние оценки
Для построения двусторонних оценок температурного коэффициента линейного расширения рассматриваемого композита используем двойственную вариационную форму модели линейной термоупругости применительно к неоднородному линейно упругому твердому телу [36, 38]. Эта форма содержит два альтернативных функционала (минимизируемый и максимизируемый), которые на истинных распределениях перемещений и напряжений достигают равных по значению экстремумов.
Рассмотрим представительный объем V композита, заключенный в абсолютно жесткую оболочку, обеспечивающую при однородном приращении температуры композита на величину ДТ отсутствие перемещений на поверхности, ограничивающей этот объем. В этом случае на допустимом для минимизируемого функционала однородном распределении и(М) = 0 (М Е V) перемещений он примет вид [43]
9 г / ч2 -V / N N
-у K'(M)(a'(M}ДТ) dV(M) = — (AT)2\Кта2т(1 - CV) + £Кa\C4
v ?=1 y
где K'(M) и a'(M) — зависящие от положения точки M G V значения объемного модуля упругости и температурного коэффициента линейного расширения в представительном объеме композита.
Для максимизируемого функционала одним из допустимых распределений напряжений является однородное по объему V, определяемое шаровым тензором s' = aSj, а = const. На таком распределении этот функционал имеет вид [43]
-=-] Ым+aa'(M )AT)dV (M >.
J2 = -
Значение а найдем из необходимого условия существования экстремума этого функционала
■ = - /(-КОМ) + а'<М)ДТ)(М>=0-
Отсюда следует
1 — C N C \-1/ N
а = -9AT( + £ КО (ат(1 - Cv)+ £аC
После подстановки а в функционал J2 получим
оV /1 — С N с )-!/ N ) 2
■ = -г(Д^1-^ + £ - + £ аС) ■ (12)
Для среды, имеющей искомый температурный коэффициент а линейного расширения рассматриваемого композита, но неоднородной по упругим характеристикам, минимизируемый функционал будет равен
QV ( N ) QV
■ = — (аДТ)2 уКт(1 - СУ) + £ КяС,] = —(аДТ)2Ку, (13)
а максимисзируемый функционал на допустимом однородном распределении напряжений
"*=-9аЛТ (+£ К )-1
примет значение
Отметим, что из сравнения а* с а и сопоставления формул (14) и (12) можно получить оценку
N
а* = ат(1 - Cv) + asC, (15)
?=i
которая также следует из теории смесей [44].
Из свойств двойственной вариационной формы модели линейной термоупругости для значений альтернативных функционалов на допустимых распределениях перемещений и напряжений следуют не только неравенства J1 ^ J2, J* ^ J2*, но и J1 ^ J2*, J* ^ J2. Из двух последних неравенств с учетом формул (11)—(15) следуют гарантированные двусторонние оценки
1 ( N ) 1/2 К~
а+ = [Кта2т{1 — Cv) + ^КalC^ > а > а*^-R = а_. (16)
Можно ожидать, что эти оценки будут достаточно грубыми, поскольку при вычислении функционалов были использованы простые и поэтому весьма приближенные допустимые распределения перемещений и напряжений.
4. Результаты расчетов
Количественный анализ полученных расчетных зависимостей проведем для композита, модифицированного включениями лишь из одного материала, т. е. при я = N =1. Введем безразмерные параметры а = а\/ат, а = а/ат, К = Кх/Кт и С = С/Кт. Тогда из формулы (15) получим
Г\!
(17)
а
а* = — = 1 — CV + а
am
из формулы (10) —
а = а* +
(а — 1)(К — 1)Cv (1 — Cv )
Kv + :
4G
где KV = 1 — CV + KCV, а соотношение (16) примет вид
(18)
а+
1 — Cv + K Cx2CV
К
> а > а*
R
Kr
Kv
а_,
(19)
где
К я =_1_..
1 - Су + Су /к
Поскольку Кя/Ку ^ 1 как отношение нижней и верхней оценок объемного модуля упругости композита [39], из соотношения (19) следует, что а- ^ а*. При одинаковых значениях объемных модулей упругости включений и матрицы (К = 1) имеем Кя/Ку = 1 и а- = а*, а верхняя оценка а + = ^1-Су+<-2Су > а* при а = 1 и Су € (0, 1). В этом случае, согласно формуле (18), а = а*.
Если совпадают модули сдвига включений и матрицы (С1 = Ст = С), то для объемного модуля упругости К * такого композита известна точная формула [2]
^ К* 3К + 4С Ку
К = 7Т- = —=-- - .
т 3 К + 4С К я
В этом случае получена и точная формула для температурного коэффициента а* линейного расширения композита [1,2]. Эту формулу можно привести к виду
а* - а - 1 ( 1 1 \
а = — = а +--- —---— .
ат 1 - 1/К \Кя К*/
Поскольку точное значение К в общем случае не известно, последнюю формулу можно применить для построения двусторонних оценок. Например, при использовании вместо значения К * его двусторонних оценок Кя или Ку в первом случае оценка совпадает со значением <5 *, а во втором случае имеет вид
а*0 = 5* + а - - . (20)
0 1 - 1/К- \Кц Ку)
Оценка а0 будет нижней при а > 1 и К < 1 или а < 1 и К > 1 и верхней при а > 1 и К > 1
или а < 1 и К < 1.
На рис. 1 для фиксированного знчения а = 2 представлены построенные по формулам (17) и (19) графики звисимостей а * (пунктирная прямая), а+ (штрихпунктирные кривые) и а- (штриховые кривые) от объемной концентрации Су включений. Кривые для оценок без символов построены при К = 1,9, а с темными кружками — при К = 0,3. По формуле (20) рассчитаны зависимости а0 от Су при К = 1,9 и К = 0,3 (соответственно штрих-пунктирная и штриховая кривые со светлыми кружками). Сплошные кривые соответствуют рассчитанным по формуле (18) зависимостям а от Су при фиксированном значении а = 2 и сочетаниях К = 1,9, (С = 2 (кривые без символов), К=1,9, С = 0,7 (кривые с квадратами), К = С = 0,3 (с ромбами), К = 0,3, С = 0,1 (с треугольниками). Видно, что все сплошные кривые расположены в пределах границ, определяемых соответствующими двусторонними оценками. При уменьшении параметра С сплошные кривые приближаются к пунктирной прямой. В случае К = 1,9 область возможных значений а между пунктирной линией и
о 0,2 0,4 0,6 0,8 Су
Рис. 1. Зависимости безразмерных оценок температурного коэффициента линейного расширения композита от объемной концентрации Су включений при <5 = 2
0 0,2 0,4 0,6 0,8 Су
Рис. 2. Зависимости безразмерных оценок температурного коэффициента линейного расширения композита от объемной концентрации Су включений при а = 0,5
графиком для верхней оценки йд заметно меньше области, ограниченной оценками, определяемыми соотношением (19). Однако при К = 0,3 эти области примерно одинаковы, причем первая из них лишь несколько смещена в сторону меньших возможных значений а.
На рис. 2 при аналогичных с предыдущим рисунком обозначениях приведены графики тех же зависимостей для фиксированного значения а = 0,5. Кривые двусторонних оценок без символов построены при К = 3, а с темными кружками — при К = 0,6. По формуле (20) рассчитаны зависимости йд от Су при К = 3 и К = 0,6 (соответственно штриховая и штрихпунктирная кривые со светлыми кружками). Построенные по формуле (18) сплошные кривые для й соответствуют сочетаниям К = 3, (5 = 2 (кривые без символов), К = 3, (5 = 0,7 (кривые с квадратами), К = 0,6, (5 = 1,2 (с ромбами), К = 0,6, (5 = 0,4 (с треугольниками). Как и на рис. 1, все сплошные кривые не выходят за пределы границ, определяемых соответствующими двусторонними оценками, а при уменьшении параметра (5 приближаются к пунктирной прямой. Соотношение между областями возможных значений а аналогично рис. 1.
Следует отметить, что по мере увеличении отклонения значений параметров а и К от единицы область возможных значений а, определяемая соотношением (19), увеличивает свои размеры быстрее по сравнению с областью, границы которой определены неравенствами (20). При этом несмотря на выполнение равенств а + = а + = 1 при Су = 0 и а + = а + = а при Су = 1 для промежуточных значений объемной концентрации включений эти оценки могут выходит за пределы промежутка, ограниченного значениями 1 и а. Поэтому в случае рассматриваемого композита с включениями только из одного материала целесообразно отдать предпочтение оценкам, определяемым неравенствами (20).
В работе [1] полученную на основе сингулярного приближения теории случайных функций формулу, аналогичную по структуре формуле (18), предложено использовать для построения двусторонних оценок путем подбора параметра О. Один из вариантов состоит в использовании двух значений О у = 1 — Су + О1Су и О и = (1 — Су + Су/01)-1, где 01 = 01/0т, 01 и 0т — модули сдвига включений и материала матрицы соответственно. Первое значение является верхней оценкой безразмерного модуля сдвига композита [34] и его использование в формуле вида (18) должно дать верхнюю оценку для а, а второе значение, будучи нижней оценкой этого модуля, должно привести к нижней оценке для а. Поскольку второе слагаемое в правой части формулы (18) является строго монотонной функцией параметра О, при любом значении О е (О О у) результат вычислений по этой формуле будет заключен между указанными двусторонними оценками.
Заключение
Установлено, что результаты вычисления температурного коэффициента линейного расширения композита с изотропными шаровыми включениями по формуле, полученной в данной работе методом самосогласования, не выходят за пределы двусторонних оценок этого коэффициента. Эта формула определяет зависимость указанного коэффициента от объемных концентраций включений из набора различных материалов и их термоупругих свойств, а также от термоупругих характеристик матрицы композита. Представленные в статье расчетные зависимости могут быть использованы для прогноза температурного коэффициента линейного расширения дисперсно армированных композитов.
Работа выполнена по грантамНШ-1432.2014.8 иМК-6573.2015.8 программ Президента РФ государственной поддержки ведущих научных школ и молодых кандидатов наук, а также в рамках проекта 1712 в сфере научной деятельности в части государственного задания № 2014/104 Минобрнауки РФ и государственного задания по проекту № 1.2640.2014.
Список литературы
1. Шермергор Т.Д. Теория упругости микронеоднородных сред. М.: Наука, 1977. 400 с.
2. Механика композиционных материалов: пер. с англ. / под ред. Дж. Сендецки: М.: Мир, 1978. 564 с.
3. Кристенсен Р. Введение в механику композитов: пер. с англ. М.: Мир, 1982. 336 с.
4. Хорошун Л.П., Солтанов Н.С. Термоупругость двухкомпонентных смесей. Киев: Нау-кова думка, 1985. 109 с.
5. Ванин Г.А. Микромеханика композиционных материалов. Киев: Наукова думка, 1985. 304 с.
6. Anglin B.S., Lebensohn R.A., Rollett A.D. Validation of a numerical method based on Fast Fourier Transforms for heterogeneous thermoelastic materials by comparison with analytical solutions // Computational Materials Science. 2014. Vol.87. P. 209-217. DOI: 10.1016/j.commatsci.2014.02.027
7. Lv J., Yang K., Zhang H., Yang D., Huang Y. A hierarchical multiscale approach for predicting thermo-electro-mechanical behavior of heterogeneous piezoelectric smart materials // Computational Materials Science. 2014. Vol.87. P. 88-99. DOI: 10.1016/j.commatsci.2014. 01.059
8. Hill R. A self-consistent mechanics of composite materials // Journal of the Mechanics and Physics of Solids. 1965. Vol. 13, no. 4. P. 213-222. DOI: 10.1016/0022-5096(65)90010-4
9. Паньков A.A. Методы самосогласования механики композитов. Пермь: Изд-во Перм. гос. техн. ун-та, 2008. 253 с.
10. Зарубин B.C., Кувыркин Г.Н., Савельева И.Ю. Оценка эффективной теплопроводности композита с шаровыми включениями методом самосогласования // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. №9. С. 435-444. DOI: 10.7463/0913.0601512
11. Зарубин B.C., Кувыркин Г.Н., Савельева И.Ю. Оценка эффективной теплопроводности однонаправленного волокнистого композита методом самосогласования // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. №11. С. 519-532. DOI: 10.7463/1113.0622927
12. Зарубин B.C., Кувыркин Г.Н. Оценки упругих свойств композита с анизотропными шаровыми включениями // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 8. С. 237-255. DOI: 10.7463/0814.0720691
13. Эшелби Дж. Континуальная теория дислокаций: пер. с англ. М.: Изд-во иностр. лит., 1963. 248 с.
14. ГоловинН.Н., ЗарубинВ.С., КувыркинГ.Н. Оценка эффективных упругих характеристик материалов, модифицированных фуллереном // Композиты и наноструктуры. 2011. № 4. С. 21-31.
15. Лурье С.А., Миронов Ю.М., Нелюб B.A., Бородулин А.С., Чуднов И.В., Буянов И.А., Соляев Ю.О. Моделирование зависимостей физико-механических характеристик от параметров микро- и наноструктуры полимерных композиционных материалов // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2012. №6. С. 38-60. DOI: 10.7463/0612.0431339
16. Lehmann В., Schlarb A.K., Friedrich K., Zhang M.Q., Rong M.Z. Modelling of mechanical properties of nanoparticle-filled polyethylene // International Journal of Polymeric Materials. 2008. Vol. 57, no. 1.P. 81-100.
17. Ivanisenko Y., Darbandi A., Dasgupta S., Kruk R., Hahn H. Bulk Nanostructured Materials: Non-Mechanical Synthesis // Advanced Engineering Materials. 2010. Vol. 12, no. 8. P. 666676. DOI: 10.1002/adem.201000131
18. Лурье C.A., Соляев Ю.О. Модифицированный метод Эшелби в задаче определения эффективных свойств со сферическими микро- и нановключениями // Вестник ПГТУ Механика. 2010. № 1. С. 80-90.
19. Зайцев А.В., Фукалов A.A. Эффективные модули объемного сжатия дисперсно-упрочненных композитов со сплошными и полыми анизотропными сферическими включениями//Вестник ПГТУ. Механика. 2010. №4. С. 46-54.
20. Prnger S., Mehlhorn L., Muhlich U., Kuna M. Study of Reinforcing Mechanisms in TRIP-Matrix Composites under Compressive Loading by Means of Micromechanical Simulations // Advanced Engineering Materials. 2013. Vol.15, no. 7. P. 542-549. DOI: 10.1002/adem.201200323
21. Соколов И.И., Долматовский М.Г. Сферопластики // Полимерные материалы. 2005. №9. С. 20-21.
22. Соколов И.И., Долматовский М.Г., Деев И.С., Стеценко В.Я. Влияние физико-механических характеристик полых стеклянных микросфер на свойства сферопластиков //Пластические массы. 2005. № 7. С. 16-18.
23. Weise J., SalkN., Jehring U., Baumeister J., Lehmhus D., Bayoumi M.A. Influence of Powder Size on Production Parameters and Properties of Syntactic Invar Foams Produced by Means of Metal Powder Injection Moulding // Advanced Engineering Materials. 2013. Vol. 15, no. 3. P. 118-122. DOI: 10.1002/adem.201200129
24. Cunha S., Aguiar J.B., Ferreira V.M., Tadeu A. Influence of the Type of Phase Change Materials Microcapsules on the Properties of Lime-Gypsum Thermal Mortars // Advanced Engineering Materials. 2014. Vol. 16, no. 4. P. 433-441. DOI: 10.1002/adem.201300278
25. Jin H.-J., Weissmuller J. Bulk Nanoporous Metal for Actuation // Advanced Engineering Materials. 2010. Vol. 12, no. 8. P. 714-723. DOI: 10.1002/adem.200900329
26. Schmidt K., Becker J. Generating Validated 3D Models of Microporous Ceramics // Advanced Engineering Materials. 2013. Vol. 15, no. 1-2. P. 40-45. DOI: 10.1002/adem.201200097
27. Giraud A., Sevostianov I. Micromechanical modeling of the effective elastic properties of oolitic limestone // International Journal of Rock Mechanics and Mining Sciences. 2013. Vol. 62. September. P. 23-27.
28. Cho Y.J., Wook Jin Lee W.J., Park S.K., Park Y.H. Effect of Pore Morphology on Deformation Behaviors in Porous Al by FEM Simulations // Advanced Engineering Materials. 2013. Vol. 15, no. 3. P. 166-169. DOI: 10.1002/adem.201200145
29. Deqing W. Relation of Cell Uniformity and Mechanical Property of a Close Cell Aluminum Foam // Advanced Engineering Materials. 2013. Vol. 15, no. 3. P. 175-179. DOI: 10.1002/adem.201200135
30. Montero-Chacon F., Marin-Montin J., Medina F. Mesomechanical characterization of porosity in cementitious composites by means of a voxel-based finite element model // Computational Materials Science. 2014. Vol. 90. P. 157-170. DOI: 10.1016/j.commatsci.2014.03.066
31. GoehlerH., JehringU., Meinert J., HauserR., QuadbeckP., Kuemmel K., Stephani G., Kieback
B. Functionalized Metallic Hollow Sphere Structures// Advanced Engineering Materials. 2014. Vol. 16, no. 3. P. 335-339. DOI: 10.1002/adem.201300057
32. Schumacher Th.C., Klein T.Y., Treccani L., Rezwan K. Rapid Sintering of Porous Monoliths Assembled from Microbeads with High Specific Surface Area and Multimodal Porosity // Advanced Engineering Materials. 2014. Vol.16, no. 2. P. 151-155. DOI: 10.1002/adem.201300220
33. Шумов A.B., Плахотниченко A.A., Титова М.В. Применение композитных материалов в радиолокационной технике. Отражатель антенного устройства доплеровского метеорологического радиолокатора // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 12. С. 593-606. DOI: 10.7463/1214.0749463
34. Синев Л.С. Оценка механических напряжений в соединенных при повышенной температуре кремнии и стекле // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 12. С. 951-965. DOI: 10.7463/1214.0745310
35. Димитриенко Ю.И., Сборщиков CB., Еголева Е.С., Матвеева A.A. Моделирование термоупругих характеристик композитов на основе алюмо-хромофосфатных связующих // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 11. С. 497-518. DOI: 10.7463/1113.0623564
36. Димитриенко Ю.И., Сборщиков CB., Беленовская Ю^., Aнискович B.A., Перевислов
C.Н. Моделирование микроструктурного разрушения и прочности керамических композитов на основе реакционно-связанного SiC // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 11. С. 475-496. DOI: 10.7463/1113.0659438
37. Димитриенко Ю.И., Федонюк Н.Н., Губарева E.A., Сборщиков CB., Прозоровский A.A. Многомасштабное конечно-элементное моделирование трехслойных сотовых композитных конструкций//Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. №7. С. 243-265. DOI: 10.7463/0714.0717805
38. Шаскольская М.С. Кристаллография. М.: Bысшая школа, 1976. 392 с.
39. Зарубин B.C, Кувыркин Г.Н., Савельева И.Ю. Сравнительный анализ оценок модулей упругости композита. Изотропные шаровые включения // Bестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2014. №5. С. 53-69.
40. Зарубин В.С., Кувыркин Г.Н., Савельева И.Ю. Сравнительный анализ оценок модулей упругости композита. Анизотропные шаровые включения //Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2014. №6. С. 20-31.
41. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 512 с.
42. Введение в микромеханику: пер. с япон. / под ред. М. Онами. М.: Металлургия, 1987. 280 с.
43. Зарубин В.С., Станкевич И.В. Расчет теплонапряженных конструкций. М.: Машиностроение, 2005. 352 с.
44. Головин Н.Н., Зарубин В.С., Кувыркин Г.Н. Смесевые модели механики композитов. Ч. 1. Термомеханика и термоупругость многокомпонентной смеси // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2009. №3. С. 36-49.
Science and Education of the Bauman MSTU, 2015, no. 2, pp. 197-215.
DOI: 10.7463/0215.0757059
Received: 09.02.2015
Revised: 21.02.2015
© Bauman Moscow State Technical University
Self-Consistency Method to Evaluate a Linear Expansion Thermal Coefficient of Composite with Dispersed Inclusions
Zarubin V. S1*, Kuvyrkin G. N.1, Savelieva I. Yu.1 * fn2@bmstu.ru
*Bauman Moscow State Technical University, Russia
Keywords: composite, method of self-consistency, two-sided estimates, isotropic spherical inclusion, temperature coefficient of linear expansion
The rational use of composites as structural materials, while perceiving the thermal and mechanical loads, to a large extent determined by their thermoelastic properties. From the presented review of works devoted to the analysis of thermoelastic characteristics of composites, it follows that the problem of estimating these characteristics is important. Among the thermoelastic properties of composites occupies an important place its temperature coefficient of linear expansion.
Along with fiber composites are widely used in the technique of dispersion hardening composites, in which the role of inclusions carry particles of high-strength and high-modulus materials, including nanostructured elements. Typically, the dispersed particles have similar dimensions in all directions, which allows the shape of the particles in the first approximation the ball.
In an article for the composite with isotropic spherical inclusions of a plurality of different materials by the self-produced design formulas relating the temperature coefficient of linear expansion with volume concentration of inclusions and their thermoelastic characteristics, as well as the thermoelastic properties of the matrix of the composite. Feature of the method is the self-accountability thermomechanical interaction of a single inclusion or matrix particles with a homogeneous isotropic medium having the desired temperature coefficient of linear expansion. Averaging over the volume of the composite arising from such interaction perturbation strain and stress in the inclusions and the matrix particles and makes it possible to obtain such calculation formulas.
For the validation of the results of calculations of the temperature coefficient of linear expansion of the composite of this type used two-sided estimates that are based on the dual variational formulation of linear thermoelasticity problem in an inhomogeneous solid containing two alternative functional (such as Lagrange and Castigliano). These functionals on the true distribution
Science ^Education
of the Bauman MSTU
of strains and stresses in an inhomogeneous body matching meaningfully reach extremes (minimum and maximum respectively). On the convergence of the distribution of the minimized functional application gives an upper estimate of the desired values of the temperature coefficient of linear expansion of the composite, and maximized use - their lower bound. The difference between these estimates allows to predict the maximum possible error in calculations obtained dependences.
References
1. Shermergor T.D. Teoriia uprugosti mikroneodnorodnykh sred [Theory of elasticity of micro-inhomogeneous media]. Moscow, Nauka Publ., 1977. 400 p. (in Russian).
2. Sendetski G., ed. Kompozitsionnye materialy. V 8 t. T.2. Mekhanika kompozitsionnykh mate-rialov [Composite Materials. In 8 vols. Vol. 2. Mechanics of Composite Materials]. Moscow, Mir Publ., 1978. 564 p. (Russian translation).
3. Christensen R.M. Mechanics of composite materials. New York, Wiley-Interscience publication, 1979. 348 p. (Russ. ed.: Christensen R.M. Vvedenie v mekhaniku kompozitov. Moscow, Mir Publ., 1982. 334 p.).
4. Khoroshun L.P., SoltanovN.S. Termouprugost' dvukhkomponentnykh smesei [Thermoelastic of twocomponent mixtures]. Kiev, Naukova dumka Publ., 1985. 109 p. (in Russian).
5. Vanin G.A. Mikromekhanika kompozitsionnykh materialov [Micromechanics of composite materials]. Kiev, Naukova dumka Publ., 1985. 304 p. (in Russian).
6. Anglin B.S., Lebensohn R.A., Rollett A.D. Validation of a numerical method based on Fast Fourier Transforms for heterogeneous thermoelastic materials by comparison with analytical solutions. Computational Materials Science, 2014, vol. 87, pp. 209-217. DOI: 10.1016/j.commatsci.2014.02.027
7. Lv J., Yang K., Zhang H., Yang D., Huang Y. A hierarchical multiscale approach for predicting thermo-electro-mechanical behavior of heterogeneous piezoelectric smart materials. Computational Materials Science, 2014, vol. 87, pp. 88-99. DOI: 10.1016/j.commatsci.2014.01.059
8. Hill R. A self-consistent mechanics of composite materials. Journal of the Mechanics and Physics of Solids, 1965, vol. 13, no. 4, pp. 213-222. DOI: 10.1016/0022-5096(65)90010-4
9. Pan'kov A.A. Metody samosoglasovaniya mekhaniki kompozitov [Self-Consistent Methods for Composites]. Perm', PSTU Publ., 2008. 253 p. (in Russian).
10. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Evaluation of effective thermal conductivity of composites with ball inclusions by the method of self-consistency. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2013, no. 9, pp. 435-444. DOI: 10.7463/0913.0601512 (in Russian).
11. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Evaluation of effective thermal conductivity of unidirectional fiber composites by the method of self-consistency. Nauka i obrazovanie MGTU
im. N.E. Baumana = Science and Education of the Bauman MSTU, 2013, no. 11, pp. 519-532. DOI: 10.7463/1113.0622927 (in Russian).
12. Zarubin V.S., Kuvyrkin G.N. Estimating of the Elastic Properties of the Composite with Anisotropic Ball Inclusions. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 8, pp. 237-255. DOI: 10.7463/0814.0720691 (in Russian).
13. Eshelby J.D. The continuum theory of lattice defects. In: Seitz F., Turnbull D., eds. Progress in Solid State Physics. Vol. 3. New York, Academic Press, 1956, pp. 79-144. (Russ. ed.: Eshelby J.D. Kontinual'naia teoriia dislokatsii. Moscow, Izd-vo inostrannoi literatury, 1963. 248 p.).
14. Golovin N.N., Zarubin V.S., Kuvyrkin G.N. Estimation of effective moduli of materials modified by fullerene. Kompozity i nanostruktury = Composites and Nanostructures, 2011, no. 4, pp. 21-31. (in Russian).
15. Lur'e S.A., Mironov Yu.M., Nelyub V.A., Borodulin A.S., Chudnov I.V., Buyanov I.A., Solyaev Yu.O. Modeling of dependences of physical-mechanical properties on parameters of micro- and nanostructure polymer composite materials. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2012, no. 6, pp. 38-60. DOI: 10.7463/0612.0431339 (in Russian).
16. Lehmann B., Schlarb A.K., Friedrich K., Zhang M.Q., Rong M.Z. Modelling of mechanical properties of nanoparticle-filled polyethylene. International Journal of Polymeric Materials, 2008, vol.57, no. 1, pp. 81-100.
17. Ivanisenko Y., Darbandi A., Dasgupta S., Kruk R., Hahn H. Bulk Nanostructured Materials: Non-Mechanical Synthesis. Advanced Engineering Materials, 2010, vol. 12, no. 8, pp. 666676. DOI: 10.1002/adem.201000131
18. Lur'e S.A., Soliaev Iu.O. Modified Eshelby method in the problem of definition of effective properties with spherical micro- and nanoincludes. Vestnik PGTU. Mekhanika = PNRPU Mechanics Bulletin, 2010, no. 1, pp. 80-90. (in Russian).
19. Zaitsev A.V., Fukalov A.A. Effective Bulk Modules of 3D Matrix-Inclusion Composites with Continuous and Hollow Anisotropic Spherical Particles. Vestnik PGTU. Mekhanika = PNRPU Mechanics Bulletin, 2010, no. 4, pp. 46-54. (in Russian).
20. Prüger S., Mehlhorn L., Mühlich U., Kuna M. Study of Reinforcing Mechanisms in TRIP-Matrix Composites under Compressive Loading by Means of Micromechanical Simulations. Advanced Engineering Materials, 2013, vol.15, no. 7, pp. 542-549. DOI: 10.1002/adem.201200323
21. Sokolov I.I., Dolmatovskii M.G. Spheroplasts. Polimernye materialy, 2005, no. 9, pp. 20-21. (in Russian).
22. Sokolov I.I., Dolmatovskii M.G., Deev I.S., Stetsenko V.Ia. Influence of the hysicomechanical characteristics of hollow glass microspheres on the properties of sphere-filled plastics. Plas-
ticheskie massy, 2005, no. 7, pp. 16-18. (English translation: International Polymer Science and Technology, 2006, vol. 33, no. 3, pp. 31-35).
23. Weise J., SalkN., Jehring U., Baumeister J., Lehmhus D., Bayoumi M.A. Influence of Powder Size on Production Parameters and Properties of Syntactic Invar Foams Produced by Means of Metal Powder Injection Moulding. Advanced Engineering Materials, 2013, vol. 15, no. 3, pp. 118-122. DOI: 10.1002/adem.201200129
24. Cunha S., Aguiar J.B., Ferreira V.M., Tadeu A. Influence of the Type of Phase Change Materials Microcapsules on the Properties of Lime-Gypsum Thermal Mortars. Advanced Engineering Materials, 2014, vol. 16, no. 4, pp. 433-441. DOI: 10.1002/adem.201300278
25. Jin H.-J., Weissmuller J. Bulk Nanoporous Metal for Actuation. Advanced Engineering Materials, 2010, vol. 12, no. 8, pp. 714-723. DOI: 10.1002/adem.200900329
26. Schmidt K., Becker J. Generating Validated 3D Models of Microporous Ceramics. Advanced Engineering Materials, 2013, vol. 15, no. 1-2, pp. 40-45. DOI: 10.1002/adem.201200097
27. Giraud A., Sevostianov I. Micromechanical modeling of the effective elastic properties of oolitic limestone. International Journal of Rock Mechanics and Mining Sciences, 2013, vol. 62, September, pp. 23-27.
28. Cho Y.J., Wook Jin Lee W.J., Park S.K., Park Y.H. Effect of Pore Morphology on Deformation Behaviors in Porous Al by FEM Simulations. Advanced Engineering Materials, 2013, vol. 15, no. 3, pp. 166-169. DOI: 10.1002/adem.201200145
29. Deqing W. Relation of Cell Uniformity and Mechanical Property of a Close Cell Aluminum Foam. Advanced Engineering Materials, 2013, vol.15, no. 3, pp. 175-179. DOI: 10.1002/adem.201200135
30. Montero-Chacon F., Marin-Montin J., Medina F. Mesomechanical characterization of porosity in cementitious composites by means of a voxel-based finite element model. Computational Materials Science, 2014, vol. 90, pp. 157-170. DOI: 10.1016/j.commatsci.2014.03.066
31. GoehlerH., JehringU., Meinert J., HauserR., QuadbeckP., Kuemmel K., Stephani G., Kieback B. Functionalized Metallic Hollow Sphere Structures. Advanced Engineering Materials, 2014, vol. 16, no. 3, pp. 335-339. DOI: 10.1002/adem.201300057
32. Schumacher Th.C., Klein T.Y., Treccani L., Rezwan K. Rapid Sintering of Porous Monoliths Assembled from Microbeads with High Specific Surface Area and Multimodal Porosity. Advanced Engineering Materials, 2014, vol. 16, no. 2, pp. 151-155. DOI: 10.1002/adem.201300220
33. Shumov A.V., Plakhotnichenko A.A., Titova M.V. Composite Technology in Radar Equipment. Dopler Meteo Radar Reflector Device. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 12, pp. 593-606. DOI: 10.7463/1214.0749463 (in Russian).
34. Sinev L.S. Mechanical Stresses Estimation in Silicon and Glass Bonded at Elevated Temperature. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 12, pp. 951-965. DOI: 10.7463/1214.0745310 (in Russian).
35. Dimitrienko Yu.I., Sborshchikov S.V., Egoleva E.S., Matveeva A.A. Modeling of thermo-elastic properties of composites with alumino-chromic phosphate matrices. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2013, no. 11, pp. 497-518. DOI: 10.7463/1113.0623564 (in Russian).
36. Dimitrienko Yu.I., Sborshchikov S.V., Belenovskaya Yu.V., Aniskovich V.A., Perevislov S.N. Modeling microstructural destruction and strength of ceramic composites based on the reaction-bonded SiC. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2013, no. 11, pp. 475-496. DOI: 10.7463/1113.0659438 (in Russian).
37. Dimitrienko Yu.I., Fedonyuk N.N., Gubareva E.A., Sborshchikov S.V., Prozorovskii A.A. Multiscale Finite-Element Modeling of Sandwich Honeycomb Composite Structures. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 7, pp. 243-265. DOI: 10.7463/0714.0717805 (in Russian).
38. Shaskol'skaia M.S. Kristallografiia [Crystallography]. Moscow, Vysshaia shkola Publ., 1976. 392 p. (in Russian).
39. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Comparative Analysis Estimates of Elastic Moduli for Composite. Isotropic Spherical Inclusions. Vestnik MGTU im. N.E. Baumana. Ser. Mashinostroenie = Herald of the Bauman Moscow State Technical University. Ser. Mechanical Engineering, 2014, no. 5, pp. 53-69. (in Russian).
40. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Comparative Analysis Modulus Elasticity Estimates for Composite Anisotropic Spherical Inclusions. Vestnik MGTU im. N.E. Baumana. Ser. Mashinostroenie = Herald of the Bauman Moscow State Technical University. Ser. Mechanical Engineering, 2014, no. 6, pp. 20-31. (in Russian).
41. Zarubin V.S., Kuvyrkin G.N. Matematicheskie modeli mekhaniki i elektrodinamiki sploshnoi sredy [Mathematical models of mechanics and electrodynamics of continuous media]. Moscow, Bauman MSTU Publ., 2008. 512 p. (in Russian).
42. Onami M., ed. Vvedenie v mikromekhaniku [Introduction to micro-mechanics]. Trans. from Japanese. Moscow, Metallurgiia Publ., 1987. 280 p. (in Russian).
43. Zarubin V.S., Stankevich I.V. Raschet teplonapryazhennykh konstruktsii [Calculation of heat-stressed designs]. Moscow, Mashinostroenie Publ., 2005. 352 p. (in Russian).
44. Golovin N.N., Zarubin V.S., Kuvyrkin G.N. Mixture Models of Composite Mechanics. P. 1. Thermal Mechanics and Thermoelasticity of Multicomponent Mixture. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki = Herald ofthe Bauman Moscow State Technical University. Ser. Natural science, 2009, no. 3, pp. 36-49. (in Russian).
