CC BY
46
УЧЕНЫЕ ЗАПИСКИ ПЕТРОЗАВОДСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Июнь, № 4 Физико-математические науки 2014
УДК 539.3
АНДРЕЙ АНАТОЛЬЕВИЧ ПАНЬКОВ
доктор физико-математических наук, доцент, профессор кафедры механики композиционных материалов и конструкций аэрокосмического факультета, Пермский национальный исследовательский политехнический университет (Пермь, Российская Федерация) mkmk_pr@pstu. ru
ДИЭЛЕКТРИЧЕСКИЕ ПРОНИЦАЕМОСТИ КОМПОЗИТА С УЧЕТОМ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ СТРУКТУРЫ И МАКСВЕЛЛ-ВАГНЕРОВСКОЙ РЕЛАКСАЦИИ*
Получено уточненное аналитическое решение для тензора эффективных проницаемостей композита с разупорядоченной полидисперсной структурой в виде разложения через известные решения для базовых структур и приведенный коэффициент корреляции, рассчитываемый на основе анализа корреляционных функций структур. Точность решения проиллюстрирована на примере расчета эффективных диэлектрических проницаемостей композита с разупорядоченной полидисперсной структурой. Проведен численный расчет и анализ влияния с учетом максвелл-вагнеровской релаксации на действительную и мнимую части эффективной диэлектрической проницаемости композита полимер/феррит с различными полидисперсными структурами величины наполнения ферритовыми волокнами и частоты приложенного электрического поля.
Ключевые слова: композит, эффективные свойства, диэлектрическая проницаемость, максвелл-вагнеровская релаксация
ВВЕДЕНИЕ
Гетерогенные радиотехнические материалы благодаря своим уникальным физико-механическим свойствам находят широкое применение в различных областях науки и техники. Композиты находят применение в тех случаях, когда традиционные материалы: кристаллы, керамика, сплавы - не обеспечивают необходимого комплекса требуемых характеристик. Возможность оптимизации и управления структурой композитов открывает путь к созданию новых материалов с наперед заданными свойствами. В результате взаимодействия на микроуровне элементов структуры на макроуровне композита возникают качественно новые эффекты по сравнению с однородными материалами. Эффективные свойства композита определяются из условия, чтобы однородное тело с такими свойствами взаимодействовало с приложенным, например, электромагнитным полем на макроуровне так же, как микронеоднородное тело той же формы из композита. Нерегулярный характер реальных структур композитов обусловливает необходимость решения этой задачи в вероятностной постановке.
1. Корреляционные функции полидисперсных структур
Рассмотрим нормированные корреляционные функции для разупорядоченной полидисперсной структуры [3] на рис. 1б
^n(p) = y|-(i1/(rX/(r1)} (1)
Mi
и для вспомогательных структур на рис. 1а, в соответственно
*,» = '(r К '(r,)), (2)
М i
к, (Р) = Mi," '(r )S '(r,)), (3)
М
где начальные значения: к, ,(0) = к,, (0) =
= к,, (0) = 1 , дисперсии:
Mi=MP = MS = чО-ч) (4)
с учетом равенств V = v1 = V объемных долей волокон во всех трех структурах, ^..^ - оператор осреднения по представительной области композита V с учетом статистической однородности и эргодичности [2] структур, вектор разности
Р = ri- r; (5)
здесь и далее верхние индексы p, sc указывают на принадлежность к вспомогательным полидисперсным структурам на рис. 1а и рис. 1в соответственно.
Распределение размеров радиусов поперечных сечений однотипных цилиндрических ячеек варьируется в широком диапазоне, включающем и бесконечно малые значения. Считаем, что точка r случайным образом лежит внутри некоторой фиксированной цилиндрической ячейки с радиусом поперечного сечения R и радиус поперечного сечения содержащегося в ней волокна rv для полидисперсной структуры на рис. 1а
© Паньков А. А., 2014
108
А. А. Паньков
0,75
0,25 -
2 2,5 3
-0,25
Рис. 1. Фрагменты (а)-(в), нормированные корреляционные функции (г) полидисперсных структур
1'
и для структуры на рис.
16 Г
R
1+ S’
где заданная величина относительной толщины
прослойки матрицы 5 =--------. В произволь-
r-
ных ячейках полидисперсной структуры на рис. 1а волокна расположены в них с вероятностью р0 = 1, для структуры на рис. 1в - с вероятностью р0 = v1, для структуры на рис. 1б - с вероятностью
р0 е jV ; 1^ . Для структуры на рис. 1б величина относительной толщины прослойки матрицы 5
может быть рассчитана 5 = —- -1 через за-
данные значения величин вероятности р0 и относительной объемной доли волокон V е ^1; v1max ). Максимальное значение объемной доли волокон в композите v1 = р и для этого предельного
случая (рис. 1в) прослойка 5 = 0.
Численный расчет двухточечных корреляционных функций (1)—(3) проведем на основе алгоритма, изложенного в [10], [12]. На рис. 1г представлены результаты расчета нормированных корреляционных функций: &п(р) (•) для структуры на рис. 1б при р0 = 0,8, ^(р) (□) для структуры на рис. 1а, ^(р) (О) для структуры на рис. 1в, расстояние между точками р = |р| (9). Дополнительно на рис. 1г приведен результат
аппроксимации (о) нормированной корреляционной функции
£ц(р) * рЛп(р)+(i-р, Кт (р) (6)
через функции £*(р) (□), ^[(р) (О) и аппроксимирующий коэффициент р. - приведенный коэффициент корреляции.
Рассмотрим метод корреляционных составляющих на примере вычисления тензора эффективных диэлектрических проницаемостей к композита с разупорядоченной полидисперсной структурой (см. рис. 1б) через известные решения для базовых структур (см. рис. 1а, в) на основе разложения (6) нормированной корреляционной функции.
2. Метод корреляционных составляющих
Постановка задачи. Для микронеоднородной представительной области V уравнение непрерывности и краевые условия
[Д (Д)Фг (r)],; = ^ Ф\Т = -fyi (7)
относительно искомого поля электрического потенциала ф(г), где коэффициенты
ВД=! */(*•)>-/ (8)
м
дифференциального оператора - это быстро осциллирующие функции координат r; индикаторная функция if (r) = 1, r е Vf; if (r) = 0, r £ Vf; kf - тензор диэлектрических проницаемостей f-й фазы и Vf - область f-й фазы 2
в V = ^ Vf (считаем, что f = 1 для волокон и
f= / \
f = 2 для матрицы композита); {гЛ = vf - относительное объемное содержание фазы f в композите, E* = ^ E - заданный вектор однородной напряженности электрического поля или «макронапряженности» композита. На межфазных поверхностях выполняются условия идеального контакта: непрерывность потенциала и нормальной к поверхности составляющей индукции электрического поля.
Решение для поля перемещений
Ф(г ) = {-r, + К (r) )Щ (9)
и для тензора эффективных диэлектрических проницаемостей к композита
д=Д) (10)
сводится к поиску решения для поля h(r). Искомый тензор входит в определяющие соотношения
Д)=Д,Д) (11)
на макроуровне композита.
Диэлектрические проницаемости композита с учетом корреляционной функции структуры...
109
Эффективные диэлектрические проницаемости. Для наглядности ограничимся в традиционном решении [2], [4], [9] для пульсаций электрического потенциала
Ф>) = Иг (r) E (12)
лишь «корреляционным приближением» или первым членом ряда
К (r) = jG(r ri ) [• (ri )],/r: +...; (13)
V
пульсации напряженностей
£» = -ф/п (r) = К ,п (r) E*,
где
Mr) = JG,>’ri)K(0]/ri +••• (14)
у ’
с учетом (12), (13); функция Грина G(r,rj) является решением уравнения
^•G>y(rJri) = -^(r-ri) (15)
и вместе со своими производными обращается на бесконечности в ноль, <5(r — гД- обобщенная дельта-функция Дирака, g - тензор диэлектрических проницаемостей «среды сравнения» [9].
Уточненное решение для тензора X, получим на основе аппроксимации (6) нормированной корреляционной функции kn (р). Для этого домножим со сверткой по двум парам индексов левую и правую части уравнения (14) на пульсации X1 (г) и осредним
(адО=-/ «.(г."-, )(4> ад,)}*, щ+-(16)
V ’
Выражение в левой части (16) представим в виде
(•ад}=(• (17)
через искомый тензор эффективных диэлектрических проницаемостей X композита с разупо-рядоченной полидисперсной структурой. В полученном уравнении
fc -(к))е; =
= -v^l - v |G„(r,r, ) 8 ки(рЩ E* +...
(18)
(l)s
используем аппроксимацию (6) нормированной корреляционной функции *n(p) через корреляционные функции kpx (p), ^ц(Р) традиционных полидисперсных структур (рис. 1а, в). Таким образом, из (18) получим
= Р.
fc-(4)K*= .
-V,(1 - v.^CL. jG.Cr.r,)8 Kipyir, ё; +...
(i)i
+ (!-a)
|G„(r,r1)-^-A17(p)*1 %+...
i drm, .
или
X' - {X) = p, {xp* - {X))+ (l - p, \sc' - {X)). (19)
В результате из (19) получим искомое уточненное решение
X* =р,Хр*+(l-p,)xsc* (20)
для тензора X эффективных диэлектрических проницаемостей композита с разупорядоченной полидисперсной структурой (рис. 1б), приведенный коэффициент корреляции р. учитывает реальный вид корреляционных функций случайной структуры (6). Известные решения Xр , Xsc в (20) могут быть получены как частные случаи обобщенного сингулярного приближения [9], в котором тензор g (15) диэлектрических проницаемостей среды сравнения последовательно приравнен к тензорам Х2, Xsc соответственно. Учет проводимостей У / фаз f = 1,2 и частоты о приложенного электрического поля через комплексную форму записи [1], [5], [6], [7], [8], [11]
Xf=X'f-i^ (21)
со
тензоров диэлектрических проницаемостей Xj-фаз с действительными частями Х'^ приводит к комплексному виду
*
Х'=Х'*-1У- (22)
о
тензора эффективных диэлектрических проницаемостей композита X и, как следствие, к возникновению на макроуровне композита дисперсии и энергетических потерь в переменных электрических полях, известных под названием «максвелл-вагнеровская релаксация».
Отметим, что решение (20) является уточнением известной формулы
X* =рХр* +(l-p)Xsc*, (23)
полученной ранее [4] в обобщенном сингулярном приближении с учетом лишь одноточечных моментов различных порядков разупорядочен-ных (квазипериодических [4]) случайных структур, где коэффициент корреляции (периодичности)
рассчитывается через корреляционный момент vn = (i^ ip^j, геометрический смысл которого -
относительное объемное содержание области пересечения включений при мысленном наложении разупорядоченной структуры на неразу-порядоченную (периодическую) структуру.
110
А. А. Паньков
Рис. 2. Концентрационные (а) и частотные (б) зависимости действительной Л и мнимой ЛЦ частей эффективной диэлектрической проницаемости
3. Численный расчет
На рис. 2 представлены результаты расче-
- л f* iff*
та действительных Л\ и мнимых Л частей эффективной диэлектрической проницаемости Л в плоскости изотропии ГГ2 композита -полиэтилен с ферритовыми волокнами с раз-упорядоченной полидисперсной структурой (рис. 1в) (о) по формуле (20) с использованием приведенного коэффициента корреляции р. в зависимости от относительного объемного содержания Vi керамических волокон в полимерной матрице (рис. 2а) и в зависимости от частоты о электрического поля (рис. 2б) в сравнении с точным решением (•) по «трехфазной модели» [3]: волокно с прослойкой матрицы в эффективной среде для заданного численного значения параметра Ро = 0,8. На рис. 2б частотные зависимости действительной Л и мнимой ЛЦ частей обозначены соответственно сплошной и пунктирной линиями. Диэлектрическая проницаемость Л(2)11 = °,5Л0) и проводимость Г(2)11 = 10 (Омхм)1 полиэтилена [7], для
феррита: у0)и =10— (Омхм)-1 [5], Л'(Г)и =10Л0, где диэлектрическая проницаемость вакуума Ло ~ 8,85 *1012 Ф/м. Графики на рис. 2а построены для предельного случая о и здесь мнимая часть к'* ^ 0 (22), на рис. 2б - для объемной доли волокон V = 0,7. Дополнительно на рис. 2а для сравнения приведены результаты расчета (Л) (23) с использованием коэффициента корреляции p (24) и для базовых полидисперсных структур на рис. 1а (□) и рис. 1б (◊). Нижняя граница (—) Хашина - Штрикмана [1] совпадает с базовым решением (□), а верхняя граница (+ )
совпадает при инверсии свойств фаз композита с нижней границей ( — ).
ЗАКЛЮЧЕНИЕ
Представлен метод корреляционных составляющих для решения стохастической краевой задачи стационарного процесса протекания, например, электро- и магнитопроницаемости, тепло- и электропроводности, диффузии в композите со случайной структурой, заданной через разупорядочивание некоторой начальной «более простой» случайной структуры. На основе анализа корреляционных функций искомое решение краевой задачи представлено через расчет соответствующих отклонений от считающегося известным решения для неразупорядоченной, в частном случае периодической, случайной структуры. Получено уточненное аналитическое решение для тензора эффективных проницаемостей композита с разупорядоченной полидисперсной структурой в виде разложения через известные решения для базовых структур и приведенный коэффициент корреляции, рассчитываемый на основе анализа корреляционных функций случайной структуры. Точность решения проиллюстрирована на примере расчета эффективных диэлектрических проницаемостей композита с разупорядоченной полидисперсной структурой. Проведен численный расчет и анализ влияния с учетом максвелл-вагнеровской
релаксации на действительную Л и мнимую Л части эффективной диэлектрической проницаемости Л1*1 композита полимер/феррит с различными полидисперсными структурами величины наполнения ферритовыми волокнами
V и частоты о приложенного электрического поля. Мнимая часть ЛЛ эффективной диэлектрической проницаемости обусловлена учетом проводимостей у(f)п фаз композита и в пределе: ^ 0, когда f )и ^ 0 или о ^ 0.
Для частотных зависимостей Л (о) полидисперсного композита характерна глубокая нормальная релаксация, обусловленная тем, что на высоких частотах не успевает накапливаться объемный заряд вблизи межфазных границ [6], [7], [8]. При нормальной релаксации абсолютное
значение действительной части Л монотонно уменьшается с ростом частоты, тогда как для мнимой части характерен четко выраженный релаксационный экстремум (см. рис. 2б) и могут быть рассчитаны время релаксации тк и релак-
сационная частота ок = 1/тк , соответствующая экстремуму мнимой части ЛЛ эффективной диэлектрической проницаемости.
* Работа выполнена при финансовой поддержке гранта РФФИ № 14-01-96004 р_урал_а.
Диэлектрические проницаемости композита с учетом корреляционной функции структуры...
111
СПИСОК ЛИТЕРАТУРЫ
1. Виноградов А. П. Электродинамика композитных материалов. М.: Эдиториал УРСС, 2001. 208 с.
2. Волков С. Д., Ставров В. П. Статистическая механика композитных материалов. Минск: Изд-во Белорус. гос. ун-та, 1978. 208 с.
3. Кристенсен Р. Введение в механику композитов. М.: Мир, 1982. 334 с.
4. Паньков А. А. Механика пьезокомпозитов. Электро- и магнитоупругость неоднородных сред. LAP LAMBERT Academic Publishing, 2011. 476 c.
5. Петров В. М., Бичурин М. И., Srinivasan G. Максвелл-вагнеровская релаксация в магнитоэлектрических композиционных материалах // Письма в ЖТФ. 2004. Т. 30. № 8. С. 81-87.
6. Соцков В. А. Экспериментальная оценка концентрационной зависимости действительной части диэлектрической проницаемости в неупорядоченной макросистеме парафин-графит // Письма в ЖТФ. 2004. Т. 30. № 12. С. 1-5.
7. Турик А. В., Радченко Г. С. Гигантский пьезоэлектрический эффект в слоистых композитах сегнетоэлектрик-полимер // Физика твердого тела. 2003. Т. 45. № 9. С. 1676-1679.
8. Турик А. В., Радченко Г С., Чернобабов А. И., Турик С. А. Диэлектрическая проницаемость полимерных матриц, содержащих изолированные включения: гигантское диэлектрическое усиление вместо коллективного резонанса // Письма в ЖТФ. 2004. Т. 79. № 9. С. 512-514.
9. Шер мер гор Т. Д. Теория упругости микронеоднородных сред. М.: Наука, 1977. 399 с.
10. Monet to I., Drug an W. J. A micromechanics-based nonlocal constitutive equation for elastic composites containing randomly oriented spheroidal heterogeneities // Journal of the Mechanics and Physics of Solids. 2004. Vol. 52. P. 359-393.
11. Pan’kov A. A. Maxwell-wagner relaxation in fibrous polydisperse magnetoelectric piezocomposites // Mechanics of Composite Materials. 2013. Vol. 49. № 1. P. 45-50.
12. Torquato S., Stell G. Micro structure of two^hase random media. I. The nqjoint рп^ЬЕ^у functions // Journal of Chemical Physics. 1982. Vol. 77. № 4. P. 2071-2077.
Pan’kov A. A., State National Research Politechnical University of Perm (Perm, Russian Federation)
COMPOSITE’S DIELECTRIC PERMEABILITY INCLUSIVE OF ITS CORRELATED STRUCTURAL FUNCTION AND MAXWELL-WAGNER RELAXATION
The specified analytical decision for the tensor of effective dielectric permeability of a composite with a disordered polydisperse structure in the form of decomposition through known decisions for basic structures and the specified coefficient of correlation is received. The correlation coefficient was obtained on the basis of the analysis of structures’ correlation functions. The calculation accuracy is illustrated on the example of calculation of effective dielectric permeability of the composite with disordered polydisperse structure. Numerical calculations and the influence analysis taking into account Maxwell-Wagner relaxation on the real and imaginary parts of effective dielectric permeability of the polymer/ferrite composite with various polydisperse structures filled ^ by different volumes of ferrite fibers and frequencies of enclosed electric fields are carried out.
Key words: composite, effective properties, dielectric permeability, Maxwell-Wagner relaxation
REFERENCES
1. Vinogradov A. P Elektrodinamika kompozitnykh materialov [Electrodynamics of composite materials]. Moscow, Editorial URSS Publ., 2001. 208 р.
2. Volkov S. D., Stavrov V. P Statisticheskaya mekhanika kompozitnykh materialov [Statistical mechanics of composite materials]. Minsk, Izd-vo Belorus. gos. un-ta, 1978. 208 р.
3. Kristensen R. Vvedenie v mekhaniku kompozitov [Introduction into mechanics of composites]. Moscow, Mir Publ., 1982. 334 р.
4. Pan’kov A. A. Mekhanika p’ezokompozitov. Elektro- i magnitouprugost’ neodnorodnykh sred [Mechanics of piezocomposites. Electro- and magnetoelasticity of non-uniform environments]. LAP LAMBERT Academic Publishing, 2011. 476 р.
5. Petrov V. M., Bichurin M. I Srinivasan G. Maxwell-Wagner relaxation in magnetoelectric composite materials [Maksvell-vagnerovskaya relaksatsiya v magnitoelektricheskikh kompozitsionnykh materialakh]. Pis’ma v ZhTF [JETP Letters]. 2004. Vol. 30. № 8. P 81-87.
6. Sotskov V. A. Experimental assessment of concentration dependence of the real part of dielectric permeability in the disordered macrosystem of paraffin-graphite [Eksperimental’naya otsenka kontsentratsionnoy zavisimosti deystvitel’noy chasti dielektricheskoy pronitsaemosti v neuporyadochennoy makrosisteme parafin-grafit]. Pis’ma v ZhTF [JETP Letters]. 2004. Vol. 30. № 12. P 1-5.
7. Turik A. V, Radchenko G. S. Huge piezoelectric effect in layered composites of ferroelectric material^^mer [Gi-gantskiy p’ezoelektricheskiy effekt v sloistykh kompozitakh segnetoelektrik-polimer]. Fizika tverdogo tela [Physics of the Solid State]. 2003. Vol. 45. № 9. P 1676-1679.
8. Turik A. V., Radchenko G. S., Chernobabov A. I., Turik S. A. Dielectric permeability of the polymeric matrixes containing isolated inclusions: huge dielectric strengthening instead of collective resonance [Dielektricheskaya ргот-tsaemost’ polimernykh matrits, soderzhashchikh izolirovannye vklyucheniya: gigantskoe dielektricheskoe usilenie vmesto kollektivnogo rezonansa]. Pis’ma v ZhTF [JETP Letters]. 2004. Vol. 79. № 9. P. 512-514.
9. Shermergor T. D. Teoriya uprugosti mikroneodnorodnykh sred [Theory of elasticity of micro non-uniform media]. Moscow, Nauka Publ., 1977. 399 р.
10. Monet to I., Drug an W. J. A micromechanics-based nonlocal constitutive equation for elastic composites containing randomly oriented spheroidal heterogeneities // Journal of the Mechanics and Physics of Solids. 2004. Vol. 52. P 359-393.
11. Pan’kov A. A. Maxwell - Wagner relaxation in fibrous polydisperse magnetoelectric piezocomposites // Mechanics of Composite Materials. 2013. Vol. 49. № 1. P 45-50.
12. Torquato S., Stell G. Micro structure of two-phase random media. I. The n^oint рп^ЬЕ^у functions // Journal of Chemical Physics. 1982. Vol. 77. № 4. P 2071-2077.
Поступила в редакцию 30.10.2013
