Теплопроводность текстурированного композита с анизотропными включениями в виде эллипсоидов вращения Текст научной статьи по специальности «Математика»

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 48211. Государственная регистрация №0421200025. ИБН 1994-0408

электронный научно-технический журнал

Теплопроводность текстурированного композита

с анизотропными включениями

в виде эллипсоидов вращения

# 06, июнь 2013

Б01:10.7463/0613.0569312

Зарубин В. С., Кувыркни Г. Н.

УДК 541.124

Россия, МГТУ им. Н.Э. Баумана fn2@bmstu.ru

Введение

Текстура, возникающая при модификации композита армирующими элементами, может существенным образом повлиять на величину эффективных коэффициентов теплопроводности такого композита. В идеализированном случае одинаковой ориентации анизотропных включений эллипсоидальной формы, главные оси тензора теплопроводности которых соосны с осями симметрии эллипсоидов, значения ла (а =1, 2, 3) эффективных коэффициентов теплопроводности композита зависят лишь от соответствующих главных значений Ла этого тензора, коэффициента теплопроводности Лт матрицы и объемной концентрации Су таких включений. Но в общем случае произвольной текстуры главные значения Ла тензора эффективной теплопроводности текстурированного композита будут отличны от значений Ла.

1. Математическая модель

Рассмотрим вариант текстурированного композита с анизотропными включениями в виде эллипсоидов вращения, являющимися трансверсально изотропными [1] относительно оси вращения, которую обозначим 3. Эта ось будет одной из главных осей тензора теплопроводности включения, которой будет соответствовать главное значение Л3 этого тензора. Расположение двух остальных ортогональных главных осей О^1 и О^2 этого тензор в плоскости, перпендикулярной оси вращения, произвольно и им соответствуют равные главные значения Л1 = Л2. Таким образом, тензор теплопроводности включения будет определен лишь двумя независимыми параметрами Л1 и Л3.

Построенная в [2] математическая модель (ММ) переноса тепловой энергии в композите в предположении, что одинаково ориентированные анизотропные эллипсоидальные включения не контактируют между собой, т.е. отделены друг от друга слоем изотропного материала матрицы, позволила получить оценки для значений Ла эффективных коэффициентов теплопроводности композита в виде безразмерных соотношений

Л = Ла = 1 + (Ла - 1)(Да + (1 - Р°а)СУ) (1)

а = Лт = 1 + (Ла - 1)^а(1 - Су) ,

где Ла — Ла/Лm,

= ь^ъ^ Г (2)

а 2 У (ьа + и)/(и), ()

Ьа — полуоси эллипсоида и / (и) = ^/(Ь? + и)(Ь2 + и)(Ь3 + и). Отметим, что + Б2 + Б3 = = 1 (в частности, для шарового включения Ба = 1/3). Интегралы в формуле (2) можно выразить через эллиптические интегралы [3, 4].

Для включений в виде эллипсоидов вращения интегралы в формуле (2) представимы в элементарных функциях [5]. В случае удлиненных эллипсоидов вращения (Ь1/Ь3 = Ь2/Ь3 =

= Ь < 1)

= 1 - Ь2 АгсМуЬ), ^3 = 1 - 2Б1, (3)

1 2 2(1 - Ь2) 4(1 - Ь2)3/2’ 3 1

а в случае сплющенного эллипсоида (сфероида) при Ь > 1

Ь2 аrccos(1/b) 1

2(Ь2 - 1)3/2 - Ь2- 1’

= Б2 = —^---------^= 1 - 2Б1,

При наличии текстуры, определенной в прямоугольной декартовой системе координат с осями Ох^ («макроосями») и ортами е^ г = 1, 2, 3, необходимо для каждого эллипсоидального включения задать ориентацию его осей симметрии О'£к, к = 1, 2, 3 («микроосей»), относительно «макроосей». Эту ориентацию можно задать ортогональной матрицей А направляющих косинусов. Элементы этой матрицы связаны дополнительным соотношением (с учетом используемого здесь и далее правила суммирования по повторяющимся латинским индексам)

а1к ajk 8ij, ] 1 2 3 (4)

где 8^- — символ Кронекера (единичный тензор второго ранга), 8^- = 1 при г = ] и 8^- = 0 при г = ]. В этом случае каждая из строк матрицы А задает компоненты орта ei в ортогональных «микроосях», а при г = ] равенство (4) является условием ортогональности двух ортов ei и ej. Таким образом, из девяти элементов Щк матрицы А с учетом равенства (4) независимыми будут только три. Это связано, в частности, с тем, что любую ориентацию «микроосей» относительно «макроосей» можно задать тремя угловыми координатами Эйлера: ^ —угол собственного вращения, ф — угол прецессии и § — угол нутации [6]. Углы ф и § задают направление одной из «микроосей» (на рис. 1 оси О'£3) относительно «макроосей», а угол <^,

Рис. 1. Ориентация «микроосей» относительно «макроосей»

отсчитываемый от положения меридиана с долготой ф на сфере единичного радиуса, задает поворот репера «микроосей» вокруг этой «микрооси». Элементы матрицы A выражают через углы Эйлера следующим образом [7]:

йц = cos cos ф cos $ — sin sin ф,

a12 = — cos sin ф cos $ — sin ^ cos ф,

ai3 = cos ^ sin $, a21 = sin ^ cos ф cos $ + cos ^ sin ф,

a22 = cos ^ cos ф — sin ^ sin ф cos $, a23 = sin ^ sin$,

a31 = — cos ф sin $, a32 = sin ф sin$, a33 = cos $.

Применительно к текстурированному композиту модифицируем построенную в [2] ММ композита следующим образом. Композит представим состоящим из различно ориентированных эллипсоидальных составных частиц, в каждой из которых анизотропное включение в виде эллипсоида вращения окружено материалом матрицы, причем составная частица и включение геометрически подобны с коэффициентом подобия K = СУ3. Поскольку каждая такая составная частица является представительным элементов композита с одинаковой ориентацией эллипсоидальных включений, главные оси тензора эффективной теплопроводности такой частицы соосны с ее осями симметрии, а главные значения Л^ этого тензора определены формулой (1). Составная частица, являясь по форме эллипсоидом вращения, аналогично включению обладает трансверсальной изотропией по отношению к свойству теплопроводности, т.е. в системе координат с «микроосями» O'£h имеем Л1 = Л2.

Рассмотрим произвольное расположение составной частицы относительно координат с «макроосями» O'xi, определяемое углами Эйлера ф и $, так как поворот частицы на угол ^ относительно ее оси вращения O'£h не влияет на ее ориентацию относительно «макроосей». По компонентам Л11 = Л22 = Л1 и Л33 = Л3 тензора эффективной теплопроводности этой частицы в «микроосях» найдем компоненты Л*,- (j = 1, 2, 3) этого тензора в «макроосях», воспользовавшись элементами транспонированной матрицы A и соотношениями [7]

Л*, = ЛЫahialj , 1 = 1 2 3.

Ли = Al(all + a^) + Л3«31 = Л1(sin2 ф + cos2 ф cos2 §) + Л3 cos2 ф sin2 §, (5)

Л22 = Al(al2 + a^) + Л3«32 = Л1(cos2 ф + sin2 ф cos2 §) + Л3 sin2 ф sin2 §, (б)

Л33 = Лl (al3 + a23) + Л3а33 = Л1 sin § + Л3 cos § (7)

Л12 = л!(anal2 + «21^22) + Л3азі«з2 = (Лі - Л3) sinф cos ф sin2 §, (В)

Л13 = Л! (апаіз + а2і«2з) + Л3азі«зз = (Л! - Л3) cos ф sin § cos §, (9)

Л23 = Л1(аі2аіз + а22Й2з) + Л3аз2«зз = (Л3 - Л!) sin ф sin § cos §. (10)

2. Процедура осреднения

Для перехода к оценке эффективных коэффициентов теплопроводности композита в целом с конкретной заданной текстурой необходимо предварительно определить в общем виде процедуру осреднения компонент Aj (ф, $) тензора эффективной теплопроводности отдельный составных частиц, зависящих, согласно формулам (5)-(10) от углов ф и $, определяющих ориентацию каждой из этих частиц относительно «макроосей» Oxj. Сначала рассмотрим вариант случайной (хаотической) ориентации составных частиц, при которой любое расположение «микроосей» частиц относительно «макроосей» равновероятно. Выделим в композите объем с достаточно большим числом геометрически подобных частиц с одинаковым тензором эффективной теплопроводности. Тогда при ^ то можно считать, что относительная плотность по распределения микрообъемов композита по возможным ориентациям «микроосей» будет постоянна, а в силу независимости компонент Ajj- от угла ^ плотность распределения точек пересечения единичной сферы на рис. 1 «микроосью» будет одинакова на всей ее поверхности. Это означает, что каждому элементарному телесному углу = sin $ d$ ^ф соответствует доля п0 d^ рассматриваемого объема композита, для которой ориентация «микрооси» O'£3 определена текущими значениями углов ф и $. После интегрирования по всем возможным ориентациям получим

2п П

J п0 d^ = J n0d^y sin $d$ = 4пп0 = 1.

П 0 0

Отсюда следует, что п0 = 1/(4п), где 4п — площадь поверхности сферы единичного радиуса.

Если составные частицы обладают каким-либо свойством, заданным скалярной величиной F(ф, $), не зависящей в пределах рассматриваемого объема композита от координат Xj, то среднее значение этой величины для такого композита со случайной ориентацией частиц будет равно

2п П

F = п0 J d-ф J F(ф, $) sin$d$ = n0(F).

00

10.7463/0613.0569312

ЗбВ

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

В данном случае интерес представляет характеристика составных частиц, заданная тензором второго ранга с компонентами Л*,- (ф, $), операцию осреднения которой можно представить в виде

Тогда, используя формулы (5)-(10), находим = Л22 = Л23 = 2Л1/3 + Л3/3 и Л£2 = Л^3 = = Л23 = 0, т.е. тензор эффективных коэффициентов теплопроводности рассматриваемого объема композита с хаотической ориентацией составных частиц в виде эллипсоидов вращения является шаровым, а этот объем будет изотропным по отношению к свойству теплопроводности. В этом случае совпадают значения А** и Л°кк первых инвариантов тензоров эффективных коэффициентов теплопроводности композита и составных частиц, что может служить косвенным подтверждением корректности процедуры осреднения. Отметим, что полученный результат справедлив при такой ориентации для частиц любой формы с произвольной анизотропией [8].

Если ориентация составных частиц в рассматриваемой модели композита упорядочена, т.е. композит текстурирован, то относительная плотность п распределения по ориентациям «микроосей» частиц в виде эллипсоидов вращения зависит от угловых координат ф и $. В общем случае она может быть неоднородной в рассматриваемом объеме композита, т.е. зависеть и от координат я*. В этом случае необходимо проводить осреднение не только по ориентациям частиц, но и по объему композита. Если же функция п(ф, $) не зависит от координат я* и удовлетворяет условию нормирования (п) = 1/(4п), то при осреднении характеристики частиц, определяемой компонентами Л-(<£, ф, $) тензора теплопроводности, вместо формулы (11) следует записать

Функцию f (ф, $) = п(Ф, $)/по называют текстурной [8, 9]. В данном случае ее текущее значение равно доле объема композита, занятого составными частицами, ориентация которых определена текущими значениям аргументов этой функции. Вероятностная трактовка этой функции связана с понятием совместной плотности распределения случайного вектора [10], координатами которого являются углы Эйлера ф и $.

При хаотической ориентации составных частиц, т.е. при отсутствии у композита текстуры, зависимость от углов Эйлера исчезает и текстурная функция тождественно равна

2п

п

о

о

(11)

о

о

2п

п

/

^ф / п(Ф, ^)Л*?-(ф, $) йіп= (пЛ*).

/

(12)

о

о

единице, что соответствует переходу формулы (12) в равенство (11). Эту функцию также можно принять равной единице в случае так называемой идеальной текстуры, когда ориентация всех анизотропных частиц одинакова [2], а «микрооси» совпадают с соответствующими «макроосями», т.е. ^ = ф = $ = 0 (см. рис. 1) и матрица А станет единичной, что приведет к формулам (1).

В общем случае текстура композита может быть комбинированной, включающей набор из п дискретных идеальных текстур и непрерывно распределенную, характеризуемую текстурной функцией /1 (ф, $). Пусть для каждой из дискретных идеальных структур ори-

ентация «микроосей» определена значениями фв и '^в (в = 1, n ) углов Эйлера. Тогда вместо формулы (12) следует использовать соотношение

2п п n

Kj = По / <# / Л(ф, d)Kj sintfdtf + Е veA*j(фв, '^в),

о о

в=1

где ve — доля объема композита, в которой идеальная текстура определена значениями фв и

При этом текстурную функцию /1(ф, tf) следует нормировать согласно условию

2п п n

[ f /1(ф, tf) sin tfdtf = 1 — Е vв.

о о в=1

3. Коническая текстура композита

Среди непрерывно распределенных текстур композита с эллипсоидальными включениями можно выделить так называемую коническую текстуру [9], когда одноименные «микрооси» всех частиц являются образующими соосных круговых конических поверхностей и равномерно распределены по этим поверхностям, а направления двух остальных ортогональных «микроосей» случайны. Если ось этих поверхностей совместить с «макроосью» Ox3, то текстурная функция F(tf) в этом случае будет зависеть лишь от одной угловой координаты tf. После осреднения по углу ф из формул (5) и (6) получим

A° А°

A1i = A22 = у (1 + cos2 tf) + у sin2 tf, (13)

формула (7) для A33 останется без изменения, а из формул (8)-(10) (как и при хаотической ориентации «микроосей») следует A12 = A13 = A23 = 0, т.е. «макроось» Ox3 композита будет главной осью тензора эффективной теплопроводности композита с конической текстурой.

В рассматриваемом случае составных частиц в виде эллипсоидов вращения при осреднении допустимо ограниться интегрированием по углу tf в в интервале (0; п/2). Тогда текстурную функцию F (tf) следует нормировать из условия

п/2

[ F (tf) sin tfdtf = 1, (14)

причем равномерному распределению «микроосей» O'£3 по углу $ соответствует значение F($) = 1. Таким образом, для главных значений тензора эффективной теплопроводности композита с конической текстурой с учетом формул (7) и (13) получим

7г/2

Л- = Л2 = 1/^1 (1 + cos2 $) + Л° sin2 $) F($)sin $d$, (15)

0

7п/2

Л3 = J(Л- sin2 $ + Л° cos2 $) F($) sin$d$. (16)

0

При идеальной конической текстуре микрооси O'£3 всех составных частиц будут образующими одной конической поверхности с заданным значением y G (0; п/2] полуугла раствора конуса. Тогда из условия (14) следует F($) = 1/ sin y и формулы (15) переходят в равенства

А° ло

А- = л2 = у (1 + cos2 y) + у sin2 y, Л3 = Л- sin2 y + Л3 cos2 y. (17)

Отметим, что первая формула (17) идентична формуле (13), а вторая — формуле (7).

При y = п/2 текстуру называют кольцевой [9]. В этом случае А- = л2 = (Л- + Л°)/2 и Л3 = Л°. В случае y = 0 имеем частный случай идеальной конической текстуры, называемой аксиальной [9], для которой главные значения тензора эффективной теплопроводности композита совпадают со значениями эффективных коэффициентов теплопроводности составных частиц, т.е. А- = л2 = Л- и Л3 = Л3.

В отличие от идеальной конической текстуры для реальной текстуры возможно ее некоторое рассеяние, вызванное тем, что не все одноименные «микрооси» частиц (в данном случае оси O'£3) строго направлены по образующим одной конической поверхности. При упрощенном описании слабого рассеяния идеальной конической текстуры можно принять, что «микрооси» O'£3 составных частиц равномерно заполняют зазор между двумя соосными круговыми коническими поверхностями, образующие которых составляют с осью углы y + 8 и y — 8, причем 8 ^ y и y + 8 ^ п/2. Совместим ось этих поверхностей с «макроосью» Ox3, а зазор между ними равномерно заполним «микроосями» O£3 составных частиц. Тогда условие (14) примет вид

7+5

J F5 ($) sin $d$ = 1,

7 —5

но в пределах интервала интегрирования F5($) = F(y, 8) = const. После вычисления интеграла находим F(y, 8) = 1/(2 sin y sin 8) и вместо формул (15) и (16) получаем соответственно

А- = л2 = Л- (1 + cos2 y cos2 8 + - sin2 y sin2 8^ + sin2 y + sin2 8 — - sin2 y sin2 8^ , (18)

2 \ 3 /2\ 3 /

Л3 = Л-(sin2 y + sin2 8 — - sin2 y sin2 8^ + Л3 ^cos2 y cos2 8 + 1 sin2 y sin2 8^. (19)

В качестве примера на рис. 2 и 3 представлены построенные с использованием формул (1), (3) и (15)—(17) графики зависимостей Ад = \\/Л и Л3 = Ад/Лт от Су для некоторых конических текстур композита при различных сочетаниях значений Ь, А і и А3. Ромбы на кривых соответствуют удлиненным эллипсоидам вращения при Ь = 0,2, а квадраты — сфероидам при Ь = 5. На рис. 2 утолщенные кривые отвечают хаотической ориентации эллипсоидальных включений при значениях Аі = 8, Аз = 2 (позиция 1) и А1 = 2, Аз = 8 (позиция 3), а на рис. 3 — при значениях А1 = 0, 9, А3 = 0, 3 (позиция 1) и А1 = 0, 2, А3 = 0, 8 (позиция 3).

Т*

6

5

4

3

2

1

0 0,2 0,4 0,6 0,8 Су

Рис. 2. Графики зависимостей безразмерных эффективных коэффициентов теплопроводности композита от объемной концентрации включений при различных сочетаниях значений А1 > 1 и А3 > 1

Штриховые линии со светлыми ромбами и квадратами построены на рис. 2 и 3 при указанных сочетаниях значений А1 и А3 по формулам [11]

Л1

Л2

Л3

2 + (01(1 - р;)(1 + 0) + вз(1 - Р3°)(1 - в))Су

2 — (в1^3(1 + 0) + в3^3°(1 — 0))СУ

1 + (в1(1 - Р3)(1 - в) + в3(1 - Р)3)0))Су 1 - (в1^з(1 - в) + в3^з)0))Су ,

(20)

(21)

1,0

0,6

0,4

0,2

\

ш \ \\ \ ■ і А Л 4 а.

V •н* 3^

\ ■ чЧ V у V ’ К ч<

■ 1ч ' ч *1 Ч ' 1 < ч< ч

ч ^ *4 'м уч Ч

Ъ \ чи* * \ 'Й.ч 1

0

0,2

0,4

0,6

О,

С

V

Рис. 3. Графики зависимостей безразмерных эффективных коэффициентов теплопроводности композита от объемной концентрации включений при различных сочетаниях значений А1 < 1 и А3 < 1

где

ві =

Аі - 1

1 + (Аі - 1)£>і

вз =

Аз — 1

1 + (А3 - 1)£>3

0

^($) сое2 $ від

причем в случае хаотической ориентации включений 0 = 1/3. Эти формулы получены с использованием сингулярного приближения теории случайных функций [8] и результаты расчета по ним при 0 = 1/3 практически совпадают с результатами вычислений по формулам (15) и (16) для сфероидов при А1 > А3 > 1 (позиция 1 на рис. 2) и А1 < А3 < 1 (позиция

3 на рис. 3), а также для удлиненных эллипсоидов вращения при А3 > А1 > 1 (позиция 3 на рис. 2) и А3 < А1 < 1 (позиция 1 на рис. 3). Но для сфероидов при А3 > А1 > 1 (позиция 3' на рис. 2) и А3 < А1 < 1 (позиция 1' на рис. 3), а также для удлиненных эллипсоидов вращения при А1 > А3 > 1 (позиция 1' на рис. 2) и А1 < А3 < 1 (позиция 3' на рис. 3) отличие в результатах значительно. Корректность формул (20) и (21) вызывает сомнение вследствие несовпадения результатов расчета по ним при Су = 1 со значением (2А1 + А3)/3 для эффективного коэффициента теплопроводности материала, состоящего из хаотически ориентированных анизотропных фрагментов любой формы [8, 9].

П

П

Кольцевой текстуре (7 = п/2) отвечают построенные по первой формуле (17) сплошные кривые для зависимости А от Су при А1 = 2, Аз = 8 (позиция 2 на рис. 2) и при А1 = 0, 8, Аз = 0, 2 (позиция 2 на рис. 3). Близкие результаты при 0 = 0 дает формула (20) в первом случае для удлиненных эллипсоидов вращения, а во втором — для сфероидов. Однако для сфероидов в первом случае и удлиненных эллипсоидов вращения во втором (позиции 2' на рис. 2 и 3) расчеты по формуле (20) приводят к заниженным результатам.

Для идеальной конической текстуры при 7 = п/4, А1 = 1 и А3 = 3 на рис. 2 тонкие линии (позиция 4) соответствуют зависимостям А^ от Су, а штрихпунктирные (позиция 5) — зависимостям А3 от Су. На рис. 3 (позиции 4 и 5) представлены аналогичные графики при 7 = п/4, А1 = 0, 3 и А3 = 0,1.

Как и выше, ромбы на кривых на рис. 4 и 5 соответствуют удлиненным эллипсоидам вращения при Ь = 0, 2, а квадраты — сфероидам при Ь = 5. На рис. 4 сплошные и штрихпунктирные линии построены по формулам (18) и (19) и отвечают зависимостям соответственно А1 и А3 от угла 5 (в градусах) рассеяния конической текстуры (7 = п/4) при Су = 0, 8 (позиции 1—4) и Су = 0, 5 (позиции 5—8). При этом позициям 1,2 и 5,6 соответствуют значения А1 = 8 и А3 = 2, а остальным позициям — значения А1 = 2 и А3 = 8. Графики аналогичных зависимостей представлены на рис. 5, причем позициям 1,2 и 5, 6 отвечают значения А1 = 0, 5 и А3 = 0,1, а остальным позициям — значения А1 = 0,1 и А3 = 0, 5.

Л ^ Л 'Е'

А1, А?

4.8

4.0

3.2

2,4

1.6

Л ^ Лз_, г ж''

1 Кн л^.и4

!Ь*г.

► - *-1 р *■

► £ —

1—■ ► а —в-. -Л * 1

♦- ♦ ♦ | » ■ ♦ * > V —♦—

► ф, 4...... ♦ ■ -•■и

8

0 10 20 30 40 б°

Рис. 4. Графики зависимостей А! (сплошные линии) и А3 (штрихпунктирные линии) от угла 8 (в градусах) при различных сочетаниях значений Су, А1 > 1 и А3 > 1

0,7

0,6

0,5

0,4

0.3

Л ^

1--И- «Л

►—+ 1- . И с

■4

1--И- >—+ |-в Г' Л.л шГ ■*“

1—■ ♦—

10

20

30

40 8°

Рис. 5. Графики зависимостей Л1 (сплошные линии) и Л3 (штрихпунктирные линии) от угла 8 (в градусах) при различных сочетаниях значений Су, Л1 > 1 и А3 > 1

Отметим, что при 6 = п/4 = 45° данная коническая текстура с рассеянием эквивалентна хаотической ориентации эллипсоидальных составных частиц, причем Л1 = Л3 = Л*(Су) для каждого значения Су и каждого сочетания значений А1 и А3. Непосредственная проверка показывает, что Л*(Су) = 2А1(Су)/3 + Л3(Су)/3 для фиксированного значения Су, что является косвенным подтверждением корректности использованной выше процедуры учета рассеяния конической текстуры.

Заключение

На основе построенной в работе [2] математической модели переноса тепловой энергии в композите с анизотропными включениями эллипсоидальной формы разработана модель текстурированного композита с включениями в виде эллипсоидов вращения, главные оси тензора теплопроводности которых совпадают с осями симметрии таких включений. Эта модель в сочетании с известной процедурой осреднения характеристик текстурированной среды использована для оценки эффективных коэффициентов теплопроводности рассматриваемого композита. Приведены расчетные зависимости для этих коэффициентов в случае различных вариантов конической текстуры композита с учетом возможного рассеяния такой текстуры. Представленные результаты могут быть использованы для прогноза эффектив-

ных коэффициентов теплопроводности текстурированных композитов, модифицированных эллипсоидальными включениями (в том числе наноструктурными элементами).

Работа выполнена по гранту НШ-255.2012.8 программы Президента РФ поддержки ведущих научных школ.

Список литературы

1. Димитриенко Ю.И. Тензорное исчисление. М.: Высшая школа, 2001. 575 с.

2. ЗарубинВ.С., КувыркинГ.Н. Эффективные коэффициенты теплопроводности композита с анизотропными эллипсоидальными включениями // Наука и образование. МГТУ им.

Н.Э. Баумана. Электрон. журн. 2013. №4. DOI: 10.7463/0413.0541050.

3. Эшелби Дж. Континуальная теория дислокаций: пер. с англ. М.: Изд-во иностр. литературы, 1963. 248 с.

4. ЗарубинВ.С., КувыркинГ.Н. Эффективные коэффициенты теплопроводности композита с эллипсоидальными включениями // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2012. №3. С. 76-85.

5. Карслоу Г., Егер Д. Теплопроводность твердых тел: пер. с англ. М.: Наука, 1964. 488 с.

6. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 512 с.

7. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров: пер. с англ. М.: Наука, 1968. 720 с.

8. Шермергор Т.Д. Теория упругости микронеоднородных сред. М.: Наука, 1977. 400 с.

9. Адамеску Р.А., Гельд П.В., Митюшов Е.А. Анизотропия физических свойств металлов. М.: Металлургия, 1985. 136 с.

10. Теория вероятностей / А.В. Печинкин, О.И. Тескин, Г.М. Цветкова и др.; Под ред. В.С. Зарубина, А.П. Крищенко. М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. 456 с. (Сер. Математика в техническом университете; Вып. XVI).

11. Nan C.-W., Birringer R., Clarke D.R., Gleiter H. Effective thermal conductivity of particulate composites with interfasial thermal resistance // Journal of Applied Physics. 1997. Vol. 81, no. 10. P. 6692-6699.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Thermal conductivity of the textured composite

with anisotropic inclusions

in the form of ellipsoids of rotation

# 06, June 2013

DOI: 10.7463/0613.0569312

Zarubin V. S., Kuvyrkin G. N.

Bauman Moscow State Technical University 105005, Moscow, Russian Federation

fn2@bmstu.ru

On the basis of a developed mathematical model of heat transfer in composite with anisotropic inclusions of ellipsoidal shape, a procedure for calculating components of an effective heat-conduction tensor of the textured composite was proposed. A variant of inclusions in the form of ellipsoids of rotation was considered. The main axes of the heat-conduction tensor of such inclusions coincide with symmetric axes of an ellipsoid. In case of conic texture of a composite its possible dispersion was taken into account. The obtained results could be used for estimating effective thermal conductivity of composites modified with nanostructured elements (including carbon nanotubes). Due to electrothermal analogy these results are applicable when consider-ingcharacteristics of electrical conductivity and dielectric capacitance of textured composites.

References

1. Dimitrienko Iu.I. Tenzornoe ischislenie [Tensor calculus]. Moscow, Vysshaia shkola, 2001. 575 p.

2. Zarubin V.S., Kuvyrkin G.N. Effektivnye koeffitsienty teploprovodnosti kompozita s anizotropnymi ellipsoidal’nymi vkliucheniiami [Effective thermal conductivity coefficients of composite with anisotropic ellipsoidal inclusions]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2013, no. 4. DOI: 10.7463/0413.0541050.

3. 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.: Eshelbi Dzh. Kontinual’naia teoriia dislokatsii. Moscow, Izd-vo inostrannoi literatury, 1963. 248 p.).

4. Zarubin V.S., Kuvyrkin G.N. Effektivnye koeffitsienty teploprovodnosti kompozita s ellip-soidal’nymi vkliucheniiami [Effective coefficients of thermal conductivity of a composite with ellipsoidal inclusions]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural Sciences], 2012, no. 3, pp. 76-85.

5. CarslawH.S., Jaeger J.C. Conduction of Heat in Solids. 2nd ed. Oxford University Press, 1959. (Russ. ed.: Karslou G., Eger D. Teploprovodnost’ tverdykh tel. Moscow, Nauka, 1964. 488 p.).

6. 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.

7. Korn G., Korn T. Mathematical Handbook for scientists and engineers. Definitions, Theorems and Formulas for Reference and Review. McGraw-Hill Book Company, Inc., New-York, Toronto, London, 1961. 943 p. (Russ. ed.: Korn G., Korn T. Spravochnikpo matematike dlia nauchnykh rabotnikov i inzhenerov. Moscow, Nauka, 1968. 720 p.).

8. Shermergor T.D. Teoriia uprugosti mikroneodnorodnykh sred [Theory of elasticity of microin-homogeneous media]. Moscow, Nauka, 1977. 400 p.

9. Adamesku R.A., Gel’d P.V., Mitiushov E.A. Anizotropiia fizicheskikh svoistv metallov [Anisotropy of physical properties of metals]. Moscow, Metallurgiia, 1985. 136 p.

10. Pechinkin A.V., Teskin O.I., Tsvetkova G.M. Teoriia veroiatnostei [Probability theory]. Moscow, Bauman MSTU Publ., 2004. 456 p. (Ser. Matematika v tekhnicheskom universitete [Mathematics in Technical University], vol. 16).

11. Nan C.-W., Birringer R., Clarke D.R., Gleiter H. Effective thermal conductivity of particulate composites with interfasial thermal resistance. Journal of Applied Physics, 1997, vol. 81, no. 10, pp. 6692-6699.