CC BY
52
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Эффективные коэффициенты теплопроводности
композита с анизотропными
эллипсоидальными включениями
# 04, апрель 2013
Б01:10.7463/0413.0541050
Зарубин В. С., Кувыркин Г. Н.
УДК 536.2
Россия, МГТУ им. Н.Э. Баумана fn2@bmstu.ru
Введение
Для повышения механических характеристик таких конструкционных материалов, какими являются композиты, представляется перспективным их модификация высокопрочными и высокомодульными включениями, в том числе наноструктурными элементами (например, фуллеренами и углеродными нанотрубками [1]). При использовании модифицированных композитов в теплонапряженных элементах конструкций, работающих в условиях одновременного воздействия как механических, так и тепловых нагрузок, наряду с механическими характеристиками композита важную роль играют и его теплофизические свойства (в частности, коэффициент теплопроводности). Эффективное значение коэффициента теплопроводности композита, модифицированного армирующими элементами, зависит от их объемной концентрации Су и от соотношения между коэффициентами теплопроводности этих элементов и матрицы. В данной работе рассмотрен композит, армированный анизотропными включениями в виде эллипсоидов. Такие включения могут иметь различную природу (например, образующие новые фазы в поликристаллических материалах при их термической обработке [2, 3] или упомянутые выше наноструктурные элементы).
1. Математическая модель
Математическую модель переноса тепловой энергии в композите построим в предположении, что эллипсоидальные включения в общем случае не контактируют между собой, т.е. отделены друг от друга слоем изотропного материала матрицы. Композит полагаем состоящим из множества одинаково ориентированных анизотропных эллипсоидальных частиц с
тензором теплопроводности, главные осп которого совпадают с осями симметрии эллипсоидов. Главные значения Xk (к = 1, 2, 3) этого тензора и коэффициент теплопроводности Лт материала матрицы считаем известными.
Рассмотрим тепловое взаимодействие отдельно взятого эллипсоидального включения с неограниченным объемом окружающей его матрицы. Начало прямоугольной декартовой системы координат O£iвыберем в центре эллипсоида, причем направления координатных осей совпадают с осями симметрии эллипсоида, уравнение поверхности которого имеет вид
222 м , I 23 = 1
Ь2 + Ь2 + Ь2 = ,
где bk — полуоси эллипсоида.
Примем, что на весьма большом расстоянии от центра включения по сравнению с длиной наибольшей полуоси эллипсоида, составляющие градиента установившегося распределения температуры равны T°k (запятая с последующим нижним индексом к у обозначения T температуры означает производную по направлению оси O2k). Установившееся распределение температуры в анизотропной среде описывает дифференциальное уравнение ЛaTkk = 0 (а = к) с частными производными второго порядка (здесь и далее использовано правило суммирования по повторяющемуся латинскому индексу). Этому уравнению удовлетворяет линейное распределение температуры во включении
т (21,22,23 ) = Ak 2k (1)
с пока неизвестными составляющими Ak градиента температуры.
Установившееся распределение температуры T*(21, 22 , 23) в матрице должно удовлетворять уравнению Лапласа T*kk = 0. Известно [4], что этому уравнению удовлетворяет позволяющая выполнить также и граничные условия идеального теплового контакта на поверхности эллипсоидального включения функция
т * = T°k 2k + Ba Dk 2k, (2)
где Ba = const,
СЮ
D = M2&3 / du
Da = 2 J (ba + u)f (u), (3)
в
в — характеризующий положение точки M с координатами 2k положительный корень уравнения
222 21 I 22 I 2з = 1 (4)
bi + в + Ь2 + в + ьз + в ,
а f (u) = V(b2 + u)(b2 + u)(b3 + u).
Постоянные Ak и Bk в соотношениях (1) и (2) найдем в предположении идеального теплового контакта на поверхности эллипсоидального включения из условий непрерывности
температуры Т = Т* и плотности теплового потока АаТкпк = АтТкпк, где пк — направляющие косинусы вектора внешней нормали к этой поверхности. Из этих условий для произвольной точки поверхности эллипсоида соответственно следует
Ak £k = T°k £k + Bk £k D°°, XaAk nk = T°k nk + Bfc nfc (Da£a);Q; , (5)
в=о
где A« = A«/Am и
сю
D° = 6i66s f (6)
a 2 J (62 + u)f (u) ' ()
0
так как для точек этой поверхности, согласно уравнению (4), в = 0. Отметим, что D° + D° + D° = 1 (в частности, для шарового включения D^ = 1/3). Интегралы в формуле (6) можно выразить через эллиптические интегралы [5, 6].
Равенства (5) должны быть выполнены в любой точке поверхности включения. В точке с координатами = 6i, £2 = £3 = 0 имеем ni = 1, n2 = n3 = 0 и из этих равенств получим
Ai = T°i + BiD°, A«Ai = T°i + Bi(Di£i)
. (7)
Çi=bi, Î2=Î3=0
С учетом правил дифференцирования интеграла с переменным пределом и неявно заданной уравнением (4) функции вычислим в этой точке производную в правой части второго равенства (7):
(Di£i), i
Тогда из равенств (7) получим
= D° + 6iDi ,i
в=о
= D°a - 1.
îi=bi, Î2=Î3=0
T° T°(1 - Ai)
Ai = /т Д , Bi ,i( i)
1 + (А1 - 1)Б?' 1 1 + (А1 - ад
Аналогичным путем, выполнив равенства (5) в точках поверхности эллипсоида с координатами = 62, = = 0 и = 63, = = 0, найдем соответственно коэффициенты А2, В2 и А3, В3. Полученные результаты можно представить в виде
, Т,к д Тк(1 - Аа)
Ак = --—-., „ , в к
1 + (а„ -ад' 1 + (а„ - ад
Таким образом, наличие анизотропного эллипсоидального включения создает в матрице возмущение температурного поля относительно линейного распределения Т^к£к на большом удалении от этого включения, описываемое соотношением
ДТ* = вкБа£к = Тк£к (1 - Аа)Да . (8)
к а^к ,к1 + (А„ - 1)Ба
i
2. Случай одинаковой ориентации включений
Сначала рассмотрим случай, когда T °2 = T3 = 0. Тогда для возмущения температурного поля из формулы (8) получим
AT * = Bk Da 2k = T ° 21—(i-—. (9)
k a2k ,1U1 + (Л1 - 1)D1 W
Предположим, что все эллипсоидальные включения в композите подобны по форме и одинаково ориентированы относительно выбранной системы координат, но размеры таких включений могут быть различны. Это приведет к различию эффективных коэффициентов теплопроводности композита в направлениях различных координатных осей, т.е. к анизотропии свойств композита по отношению к теплопроводности. Пусть N таких включений находятся в объеме VN, ограниченном поверхностью эллипсоида с уравнением 2i/B2 + + 22/B22 + 2I/B3 = 1 и равном 4nB1B2 B3/3, где Ba/bOn) = Cn = const > 1 и bOn) — полуоси
эллипсоидального включения с номером п = 1, N. Поскольку объем п-го включения равен б2га) ЬзП) / 3, объемную концентрацию включений в объеме Ум можно определить
величинои
м 1
^ = £ ^ < 1. (10)
п=1 Сп
Отметим, что для подобных по форме эллипсоидальных включений значения не зависят от размера включения. Например, в соответствии с формулой (6) можно представить как функцию одинаковых для таких включений аргументов Ь2 = Ь2/Ь1 и Ь3 = Ь3/Ь1 в виде
1 ['
1
2 J (1 + u)V(1 + u)(b2 + u)(b2 + u)'
Для точки, удаленной от каждого из включений на весьма большое расстояние по сравнению со значением Вь примем для всх включений 21 ^ В1. Тогда, согласно формуле (9), N весьма удаленных включений, расположенных в объеме Ум в виде эллипсоида с полуосями Вк, вызовут в этой точке возмущение температуры, равное
дтм = V ДТП = т 12^ (1 -Л1)Б . (11)
Если считать эллипсоид объемом Ум представительным элементом композита с рассматриваемыми включениями, то этот элемент с искомым значением Л1 эффективного коэффициента теплопроводности в направлении оси 021 создаст в той же весьма удаленной точке с
координатами 2к с учетом формулы (8) такое же возмущение температуры
ДТм=^ , (12)
где Ai = A1/Am и
= BiBBi i_du__(13)
D = 2 J (B2 + u)F(u)' (13)
ß*
причем в1 — положительный корень уравнения
e2 , й , e2
в2 + в B2 + в B2 + i
а
F (u) = ^(ß? + u)(B? + u)(B2 + u). Приравняв правые части формул (11) и (12), запишем
+ + "тиТ^Т = 1 (14)
(1 - - (1 - A1> ± ^ (.5)
1 + (Ai - 1)D1 1 + (Ai - 1)D1 D1
Для весьма удаленной точки |£11 ^ то, что, согласно уравнениям (4) и (14), равносильно в ^ то, в1 ^ то и стремлению к нулю интегралов D(ra) и D1 в формуле (15). Для раскрытия неопределенности типа 0/0 используем правило Лопиталя, продифференцировав каждый из этих интегралов по нижнему переменному пределу. В итоге с учетом формул (3) и (13) получим
lim = ± lim (В? + ei)F (в) =±
ß,ß*D1 cn ß,ß*— ((б1га))2 + в)f(в) cn
Таким образом, формула (15) с учетом равенства (10) примет вид
A = 1 + (Ai - 1)(D + (1 - d;)CV) (16)
1 1 + (Ai - 1)D(1 - CV) ' ( )
Аналогичным путем можно найти формулы для A2 = A?/Am и A3 = A3/Am, где A? и Ai — эффективные коэффициенты теплопроводности композита в направлении осей O£2 и O£3 соответственно, являющиеся вместе с A 1 главными значениями тензора эффективной теплопроводности композита. В итоге при а = 1, 2, 3 имеем
A = 1 + (Äg - 1)(D + (1 - D )CV)
Aa 1 + (A« - 1)D(1 - Cv) , (17)
что в частном случае изотропных эллипсоидальных включений (Aa = Ao) совпадает с соотношением, полученным в [6].
Если включения абсолютно нетеплопроводны (Aa = 0), то из формулы (17) следует
a = (1 - D;)(1 - cv) (18)
Aa = 1 - (1 - Cv) (18)
Эта формула применима к материалу с коэффициентом теплопроводности Лт, содержащему поры с объемной концентрацией Су. При абсолютно теплопроводных включениях (Л, ^ то) формула (17) примет вид
а = + (1 - юсу
Ла = ад - с,) . (19)
Отметим, что формулы (18) и (19) также совпадают с полученными в [6].
3. Построение двусторонних оценок
Применение двойственной вариационной формулировки задачи стационарной теплопроводности в неоднородном теле [7, 8] при использовании достаточно простых допустимых для альтернативных функционалов распределений температуры и вектора плотности теплового потока позволяет получить двусторонние оценки эффективных коэффициентов теплопроводности рассматриваемого композита в виде
= г\г к < Ла < 1 - Су + Л, Су = Л+. (20)
1 - Су + Су/Л,
Нижняя (Л-) и верхняя (А+) оценки совпадают между собой и со значениями, определяемыми формулой (17), при значениях Су = 0 и Су = 1 объемной концентрации включений. Но при промежуточных значениях Су разность ДЛ, = - А- возрастает от нуля по мере отклонения параметра Ла от единицы, что дает основание предполагать при этом и возрастание возможной погрешности формулы (17). Отметим, что аналогичная тенденция выявлена и в случае композита с изотропными эллипсоидальными включениями [6].
С использованием неравенств (20) оценим наибольшую возможную относительную погрешность среднего значения = (А+ + Л-)/2 при промежуточных значениях Су, которую представим в виде
е = ДЛ = 1__1__(21)
а 2А, 1 + Су (1 - Су)(Ла + 1 /Ла - 2)/2. ( )
На рис. 1 приведены зависимости е, от Су при различных значениях Л, > 1 (кривые для пар значений Л, и 1/Л, совпадают). Наибольших значений е, = 1 -1/(1 + (Л, + 1/Л, - 2)/8) относительная погрешность достигает при С° = 0, 5.
Использованные для получения двусторонних оценок достаточно простые допустимые распределения температуры и плотности теплового потока учитывают лишь объемное содержание эллипсоидальных включений. Поэтому для всех трех направлений координатных осей представленные в формулах (20) и (21) оценки эффективных коэффициентов теплопроводности композита и возможной относительной погрешности среднего значения этих оценок идентичны.
О 0,2 0,4 0,6 0,8 Су
Рис. 1. Графики зависимости от Су при различных значениях Ла > 1
Заключение
Построенная математическая модель переноса тепловой энергии в композите с одинаково ориентированными анизотропными включениями эллипсоидальной формы позволила найти эффективные коэффициенты теплопроводности такого композита. Для оценки возможной погрешности полученных результатов применена двойственная формулировка вариационной задачи стационарной теплопроводности. Эти результаты могут быть использованы для прогноза эффективных коэффициентов теплопроводности композитов, модифицированных наноструктурными элементами (в том числе углеродными нанотрубками).
Работа выполнена по гранту НШ-255.2012.8 программы Президента РФ поддержки ведущих научных школ.
Список литературы
1. Кац Е.А. Фуллерены, углеродные нанотрубки и нанокластеры. Родословная форм и идей. М.: Изд-во ЛКИ, 2008. 296 с.
2. Научные основы материаловедения / Б.Н. Арзамасов, А.И. Крашенинников, Ж.П. Пастухова, А.Г. Рахштадт. М.: Изд-во МГТУ им. Н.Э. Баумана, 1994. 366 с.
3. Ван Флек Л. Теоретическое и прикладное материаловедение: Пер.с англ. М.: Атомиздат, 1975. 472 с.
4. Карслоу Г., Егер Д. Теплопроводность твердых тел: Пер.с англ. М.: Наука, 1964. 488 с.
5. Эшелби Дж. Континуальная теория дислокаций: Пер. с англ. М.: Изд-во иностр. лит., 1963. 248 с.
6. Зарубин В.С., Кувыркин Г.Н. Эффективные коэффициенты теплопроводности композита с эллипсоидальными включениями // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2012. №3. С. 76-85.
7. Зарубин В.С. Инженерные методы решения задач теплопроводности. М.: Энергоатомиз-дат, 1983. 328 с.
8. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 512 с.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Effective thermal conductivity coefficient
of composite with anisotropic
ellipsoidal inclusions
# 04, April 2013
DOI: 10.7463/0413.0541050
Zarubin V. S., Kuvyrkin G. N.
Bauman Moscow State Technical University 105005, Moscow, Russian Federation
fn2@bmstu.ru
A mathematical model of heat transfer in a composite material with identically oriented anisotropic ellipsoidal inclusions was created. Formulas for calculating the effective thermal conductivity of composite materials, which are anisotropic with respect to the property of thermal conductivity, were obtained on the basis of developed mathematical model. To estimate the possible error of the results the dual formulation of the variational problem of stationary heat conduction in a heterogeneous solid body was applied. These results can be used to predict the effective thermal conductivity of the composite materials with ellipsoidal inclusions, which are among others nanostructured elements (including carbon nanotubes). Obtained formulas could be used for estimating electric conductivity of composite materials with identically oriented anisotropic ellipsoidal inclusions due to electro thermal analogy .
References
1. Kats E.A. Fullereny, uglerodnye nanotrubki i nanoklastery. Rodoslovnaiaform i idei [Fullerenes, carbon nanotubes and nanoclusters. Genealogy of forms and ideas]. Moscow, LKI Publ., 2008. 296 p.
2. Arzamasov B.N., Krasheninnikov A.I., Pastukhova Zh.P., Rakhshtadt A.G. Nauchnye osnovy materialovedeniia [Scientific foundations of Materials Science]. Moscow, Bauman MSTU Publ., 1994. 366 p.
3. Van Vlack L.H. Elements of Materials Science and Engineering. 2nd ed. Addison-Wesley Publishing Co., Massachusetts, 1964. 600 p. (Russ. ed.: VanFlekL. Teoreticheskoe iprikladnoe materialovedenie. Moscow, Atomizdat, 1975. 472 p.)
4. Carslaw H.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.)
5. 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.)
6. 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 science], 2012, no. 3, pp. 76-85.
7. Zarubin VS. Inzhenernye metody resheniia zadach teploprovodnosti [Engineering methods for solving problems of thermal conductivity]. Moscow, Energoatomizdat, 1983. 328 p.
8. 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.
