CC BY
71
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Математическое моделирование диэлектрических свойств полимер-керамических композиционных материалов методом асимптотического осреднения # 10, октябрь 2013 Б01: 10.7463/1013.0623343
Димитриенко Ю. И., Соколов А. П., Маркевич М. Н.
УДК 539.8
Россия, МГТУ им. Н.Э. Баумана dimit@serv.bmstu.ru а1 sokol o@bmstu. ги mamarkevi@gmai1.com
Введение
В связи с необходимостью создания изоляционных многокомпонентных гетерогенных систем, обладающих диэлектрическим усилением, возникает вопрос о моделировании эффективных свойств композиционных материалов [9, 10, 12]. Описание воздействия электромагнитного поля на сильно неоднородные среды осуществляется с помощью уравнений с быстроосциллирующими коэффициентами, характеризующими свойства отдельных компонент композиционного материала. Основная идея метода асимптотического осреднения дифференциальных уравнений (МАО), предложенного в работах Н. С. Бахвалова, E. Санчез-Паленсии, А. Бенсуссана [1, 11, 13], состоит в замене гетерогенного материала однородной средой, обладающей усредненными физико-механическими параметрами. Быстроосциллирующие коэффициенты уравнений представляются в виде периодических функций. Проводится асимптотическое исследование решений поставленной задачи. Однако, применимость МАО ограничивается поиском эффективных характеристик композитов, обладающих достаточно простой геометрической структурой. Для расчета эффективных свойств сложноструктурированных композиционных материалов, обладающих свойством центральной симметрии (рис.1), разработана модернизация МАО приводящая задачу с периодическими граничными условиями к серии задач, так называемых Ьр , классического типа [4-8]. Решается цикл обратных задач Ьр с адаптированными компонентами в
направлениях пространства для установления эффективных характеристик композита и тангенса угла диэлектрических потерь.
Решение серии Ьр задач и поиск эффективных характеристик композита
Для композиционного материала, характерный размер которого Ь, и существует
характерный размер ячейки периодичности (ЯП) I введем малый параметр к = — << 1, и
Ь
два типа безразмерных координат: — = — - глобальные, ^ = —- = —— локальные, где хг -
Ь к I
декартовы координаты.
Исходя из периодичности структуры композита, периодический характер присваивается диэлектрической проницаемости, комплексная амплитуда которой будет рассматриваться как функция координат :е*() = е*(хг + аг), где аг - длина сторон ЯП.
Относительно потенциала р* строится асимптотического разложения по параметру к :
> (—,) = р'(и)(— ) + к>(1)(—,) + к >"(2)(—,) +... + к>'(п)(—) =
ад
= р*(0)( —) + £кУ(п)( — ).
(1)
П=1
а] а}-
Для объема ЯП У^ = : —1 ,1 = 1,3} введены фазы Уаа,а = 1,N
2 31 2
принадлежащие разному типу компонент. С использованием разложения (1), свойств периодичной функции, правил дифференцирования сложной функции рассмотрим систему уравнений, описывающую воздействие переменного электрического поля на композит
А*а (0) = о
Б *а(0) = е*а£ *а(0)
— е У
Е*а(0) = е* +>а (1)
(
а(о) _т^*
а(1) _ **(1) /г\*а(0)
= >(1),(Д'а(0) - Б** (0))пг = о
а
— е Уа
— е Уа
— еХ
а
< >а(1) >= 0,[[Бг*(0)]]Пг = 0,[>а(1)]] = 0, х. е X
■а
(2),
*(1) - [[>*(1)]] где <>'>- операция осреднения по «ячейке периодичности», -|-1-условия
периодичности, Е* -средняя напряженность композита.
Задача (2) содержит нестандартные условия нормировки в виде интегральных уравнений Вольтера первого рода и условия периодичности на границе ЯП, что создает трудности при численном ее решении. Для упрощения поиска решения вводим функции псевдопотенциала
р=1
(3)
где^(Р) = -Е* <!;р + 5(0)(&к); ЛаО^)- новые неизвестные функции от &, уже не являющиеся периодическими.
Далее вводится предположение о центральной симметрии ЯП относительно центра локальных координат, что позволяет вместо решения на всей ЯП V& решение на 1/8 ЯП
V& : V, = ^п > 0).
Используя выражения (4) получаем задачу электродинамики композита на 1/8 ЯП
В а = 0
^г ( р )/г и
) = б*аЕ*(Хр), хг е Ег*р) =^(Т)/г , Х Е ~
-е. )
(4)
« = 0 В*р) - В-))п1 = 0, Хг е£
где Е = Е п V, - поверхности контакта компонентов внутри подобласти,
координатные плоскости - Е = = 0}, торцевые поверхности «ячейки периодичности»- Е = = 1/2},. = 1,2,3.
К задаче (4) присоединяются граничные условия на 1/8 ЯП
а р т^*
| = _р_ Е" | = 0
^(р) = „ Ер ,&(р) ^ = 0
Е
2 р, (р) Е
(5)
С/г 1Е = 0,^5)/г1Е/ = 0,р * г
Задачи Ьр имеют «классический» вид относительно функций 3(0 (^), т.е. являются «контактными» задачами электростатики с граничными условиями на поверхности области V^.
Компоненты тензора эффективной комплексной диэлектрической проницаемости среды определяется как
Д*
-* = (р)
гр е * ' р
__N
где Д*р) =< ДО,) >= X |ДО)аУ -осредненная по объему ЯП индукция электрического
а=1 V (а)
поля.
Метод конечных элементов для серии Ьр задач
На объеме 1/8 ЯП, аналогично работе [5], рассмотрим вариационную формулировку задачи электродинамики
Ш^(р)/г>>*р)/гV = .ОХ) Дг^)^ , (6)
~ X!
7~л *е
где Дг( ) - комплексная амплитуда вектора индукции внешнего поля, заданного на поверхности «ячейки периодичности» композита X , на поверхности Х2 - нулевые
граничные условия; 83*) и 8(3*р),г) - виртуальное изменение псевдопотенциала и
производная виртуального изменения.
Для решения вариационного уравнения (6) используется метод конечных элементов. Область V^ разбивается на четырехузловые элементы Уе: Уе . Для каждого
е
элемента Уе записывается уравнение (6) имеет вид
{{{8(3*р)/г) те5рР)1 ¿У = |{8(3(*р))т Д*^ (7)
Уе Хе
Псевдопотенциал 3*р) аппроксимируем как
3(*р) = [Ф]{^«}, (8)
где {ф(Р )}Т = {ф(Р )(^(1)),--,ф(*р)(^ (И))}Т = {ф('р)(^(1}),ф('Р )(й (1)Х"-,ф('р)(£ (т))ф)(й (т) )}
координатный столбец значений электрического потенциала
г(1)) = {^(Р)(й( 1 ))'^(Р)(^( 1))} 6 Узлах конечного элемента, £■(;)- координата 7"го узла конечного элемента,]=1...т, т- число узлов, [Ф] - матрица функций формы.
Производные от псевдопотенциала 3*Р) представим в виде координатного столбца
{^(Р)/1>^(Р )/2,^(р )/3} = [Ь]^(
(Р) ■
(9)
где [ Ь]Т =
д 0 д 0 д 0
дй дй2 дй
0 д 0 д 0 д
дй1 дй2 дй3
- матрица операторов дифференцирования.
Т
Вариационное уравнение (7) с учетом (8) и (9) имеет вид
Ш^{ф(Р)}Т [£]Т^[£]{ф(Р)}^ = [Иф(Р)}Т [Ф]Т А*Р) ^
где [В] = [Ь][Ф]- матрица градиентов потенциала. Вынося вариацию <5{Ф(Р)} за знак
интегралов, получаем разрешающую систему линейных алгебраических уравнений (СЛАУ):
[<]{ф(Р)} = {/*}. (10)
В СЛАУ (10) матрица [К*]является локальной матрицей «жесткости» для конечного элемента, представленной выражением:
[К*] = {{{[ В]У[ В]ё¥.
где тензор диэлектрической проницаемости е . Локальный вектор правой части {/*}Т СЛАУ (10) имеет вид
{/Л = Я [Ф]Т В Р^ .
Таким образом, если собрать СЛАУ (10) для всех конечных элементов, то получим глобальную СЛАУ относительно псевдопотенциала для всей рассматриваемой области У,
[К *]{ф*р)} = {/•}, (3.34)
где [К *]- глобальная матрица жесткости, {/ *}- глобальный столбец внешних нагрузок.
Решением (3.34) получаем распределение электрического поля внутри проектируемого наноструктурированного композита.
Численный расчет эффективных диэлектрических характеристик
В качестве исследуемого образца выбран композит на основе поливинилхлорида (ПВХ) - матрица и сегнетоэлектрической керамики ЦТС-19 - наполнитель. Для определения диэлектрических свойств композита при воздействии переменного электромагнитного поля проведен расчет серии Ьр задач с применением метода конечных элементов и метода бисопряженных градиентов, на базе программного комплекса, разработанного в НОЦ «СИМПЛЕКС» МГТУ имени Н.Э. Баумана.
Для двух типов геометрий композиционных материалов - с ячейкой периодичности (ЯП) гексагонального типа с различными радиусами армирующих включений (рис. 1), рассчитана эффективная диэлектрическая проницаемость при воздействии
2 4
электрического поля диапазона частот от 10" до10 Гц. Установлена частотная зависимость мнимой и действительной части комплексной диэлектрической проницаемости композитов на основе пластиката для мелкодисперсного наполнителя ЦТС-19 характерными размерами ~1 мкм и ~10 мкм при равной объемной доле включений (рис. 2 и рис. 3).
Рис. 1. Композиционный материал с мелкодисперсным наполнителем 2-х типов: с гексагональной формой ЯП и с кубической симметрией ЯП
Оценка адекватности математической модели и точности численных методов расчета эффективных диэлектрических характеристик получены посредством сравнения расчетных характеристик действительной и мнимой частей е' = Яее и е" = 1те* эффективной диэлектрической проницаемости, а также тангенса угла диэлектрических потерь tgд = е"/е' композита с результатами экспериментальных данных [2, 3].
Рис. 2. Частотные зависимости в' и в" для композита на основе ПВХ с керамическими
включениями ЦТС-19, мкм, ф/=40%
Рис. 3. Частотные зависимости в и в композита на основе ПВХ с керамическими включениями ЦТС-19, ^~10мкм, ф=40%
Проведен расчет тангенса угла диэлектрических потерь для каждого типа композита с коэффициентами армирования ^ равными 10% и 40 %, что позволяет проанализировать энергоемкость композита и в дальнейшем установить оптимальную концентрацию керамических включений в композите (рис. 4-5). Характер поведения тангенса угла диэлектрических потерь, при варьировании объемной доли включений, свидетельствует о хорошей точности используемых математических методов и соответствии физическому смыслу математической модели.
а) б)
Рис. 4. Расчетные и экспериментальные зависимости тангенс угла диэлектрических потерь tg8 компонентов композита: (ПВХ -РУС), керамики ЦТС-19 и самого композита, ^
=10% , ё~Ю мкм (а), мкм (б)
|д(0 1д(Я
а) б)
Рис. 5. Расчетные и экспериментальные зависимости тангенс угла диэлектрических потерь tg8 компонентов композита: (ПВХ -РУС), керамики ЦТС-19 и самого композита, ф/
=40% , ё~10 мкм (а), мкм (б) .
Заключение
В данной работе проведено асимптотическое исследование эффективных диэлектрических характеристик сложноструктурированных композитов с применением МАО. Сформулирована серия локальных задач на «ячейке периодичности» композита. Проведен
расчет эффективной диэлектрической проницаемости полимерного композита с мелкодисперсными включениями керамики ЦТС-19 с применением метода конечных элементов. Проанализирована частотная зависимость диэлектрических характеристик и тангенса угла диэлектрических потерь гетерогенных структур с различными концентрациями и размерами мелкодисперсных включений при воздействии
2 4
электрического поля в диапазоне частот от 10 до10 Гц. Получено хорошее совпадение численных расчетов и экспериментальных данных для диэлектрических свойств дисперсно-армированных полимер-керамических композитов[2, 3].
Поддержка
Работа проведена при финансовой поддержке Министерства образования и науки Российской Федерации по Соглашениям о предоставлении гранта в форме субсидий № 14.B37.21.0448, 14.B37.21.1869, НИР № 1.5433.2011, № 14.132.21.1699 в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы и РФФИ № 12-08-00998.
Список литературы
1. Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. Математические задачи механики композиционных материалов. М.: Наука, 1984. 352 с.
2. Гефле О.С., Лебедев С.М., Похолков Ю.П. Частотные спектры комплексной диэлектрической проницаемости композиционных диэлектриков на основе поливинилхлорида // Известия Томского политехнического университета. 2007. Т. 310, № 1. С. 87-91.
3. Гефле О.С., Лебедев С.М., Ткаченко С.Н. Поведение полимерных композиционных материалов с наполнителем из сегнетоэлектрической керамики в электрическом поле // Известия Томского политехнического университета. 2005. Т. 308, № 4. С. 64-68.
4. Димитриенко Ю.И., Кашкаров А.И. Расчет эффективных характеристик композитов с периодической структурой методом конечных элементов // Вестник МГТУ им.
Н.Э. Баумана. Сер. Естественные науки. 2002. № 2. С. 95-108.
5. Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов // Математическое моделирование. 2012. Т. 24, № 5. С. 3-20.
6. Димитриенко Ю.И., Соколов А.П. Разработка системы автоматизированного вычисления эффективных упругих характеристик композитов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2008. № 2. С. 57-67.
7. Димитриенко Ю.И., Соколов А.П. Система автоматизированного прогнозирования свойств композиционных материалов // Информационные технологии. 2008. № 8. С. 3138.
8. Димитриенко Ю.И., Соколов А.П., Маркевич М.Н. Моделирование диэлектрических характеристик композиционных материалов на основе метода асимптотического осреднения // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 1. DOI: 10.7463/0113.0531682
9. Емец Ю.П. Эффективная диэлектрическая проницаемость трехкомпонентных композиционных материалов с анизотропной структурой // Журнал технической физики. 2005. Т. 75, вып. 2. С. 67-72.
10. Шелухин В.В., Терентьев С.А. Гомогенизация уравнений Максвелла и дисперсия Максвелла-Вагнера // Доклады Академии наук. 2009. Т. 424, № 3. С. 402-406.
11. Bensoussan A., Lions JL., Papanicalaou G Asymptotic analysis for periodic structures. Amsterdam; New York: North-Holland Pub. Co.,1978. 396 p.
12. Kettunen H., Qi J., Wallen H., Sihvola A. Frequency dependence of effective permittivity of simple dielectric composites // Proceedings of ACES 2010. Tampere, Finland. 2010. P. 248-253.
13. Sanchez-Palencia E. Book chapters // In: Homogenization techniques for composite media / Sanchez-Palencia E., Zaoui A. (eds.). Berlin, Heidelberg: Springer, 1987. P. 121-192. DOI: 10.1007/3-540-17616-0 (Ser. Lecture notes in Physics; no. 272).
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Mathematical simulation of dielectric properties of polymer-ceramic
composite materials, using the asymptotical averaging method
# 10, October 2013
DOI: 10.7463/1013.0623343
Dimitrienko Yu.I., Sokolov A.P., Markevich M.N.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
dimit@serv.bmstu.ru al sokol o@bmstu. ru mamarkevi@gmail.com
In this paper, the authors presented a mathematical model for calculating effective dielectric characteristics of complex structured composite materials influenced by an external alternating electric field. An asymptotic investigation of electro-dynamic differential equations with fast oscillating coefficients was carried out. A local electro-dynamic problem on the composite's "periodicity cell" was formulated. An effective complex permittivity of composite materials with various volume fractions of the fine-dispersed ferroelectric ceramic inclusions was calculated by the finite element method and the biconjugate gradient stabilized method. The authors analyzed frequency dependence of dielectric characteristics and dielectric loss tangent of composite ceramics. The obtained results are good congruous to experimental data.
Publications with keywords: mathematical modeling, composites, asymptotic averaging method, effective permittivity, ferroelectric ceramics
Publications with words: mathematical modeling, composites, asymptotic averaging method, effective permittivity, ferroelectric ceramics
References
1. Bakhvalov N.S., Panasenko G.P. Osrednenieprotsessov vperiodicheskikh sredakh. Matematicheskie zadachi mekhaniki kompozitsionnykh materialov [Averaging processes in periodic media. Mathematical problems of mechanics of composite materials]. Moscow, Nauka, 1984. 352 p.
2. Gefle O.S., Lebedev S.M., Pokholkov Yu.P. Chastotnye spektry kompleksnoy dielektricheskoy pronitsaemosti kompozitsionnykh dielektrikov na osnove [Complex permittivity frequency spectra of composite dielectrics on the basis of polyvinylchloride]. Izvestiya Tomskogo politekhnicheskogo universiteta [Bulletin of the Tomsk Polytechnic University], 2007, vol. 310, no. 1, pp. 87-91.
3. Gefle O.S., Lebedev S.M., Tkachenko S.N. Povedenie polimernykh kompozitsionnykh materialov s napolnitelem iz segnetoelektricheskoy keramiki v elektricheskom pole [Behaviour of polymeric composite materials filled with ferroelectric ceramics in electric field]. Izvestiya Tomskogopolitekhnicheskogo universiteta [Bulletin of the Tomsk Polytechnic University], 2005, vol. 308, no. 4, pp. 64-68.
4. Dimitrienko Iu.I. Kashkarov A.I. Raschet effektivnykh kharakteristik kompozitov s periodicheskoi strukturoi metodom konechnykh elementov [The calculation of the effective characteristics of composite materials with periodic structure by finite element method]. Vestnik MGTU im. N.E.Baumana. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2002, no. 2, pp. 95-108.
5. Dimitrienko Yu.I., Sokolov A.P. Mnogomasshtabnoe modelirovanie uprugikh kompozitsionnykh materialov [Multiscale modeling of elastic composite materials]. Matematicheskoe modelirovanie, 2012, vol. 24, no. 5, pp. 3-20.
6. Dimitrienko Iu.I., Sokolov A.P. Razrabotka sistemy avtomatizirovannogo vychisleniia effektivnykh uprugikh kharakteristik kompozitov [Development of automated technology of calculation of effective elastic characteristics of composites by method of asymptotic averaging]. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2008, no. 2, pp. 56-67.
7. Dimitrienko Yu.I., Sokolov A.P. Sistema avtomatizirovannogo prognozirovaniya svoystv kompozitsionnykh materialov [Automated forecasting of composite material properties by means of homogenization method]. Informatsionnye tekhnologii, 2008, no. 8, pp. 31-38.
8. Dimitrienko Yu.I., Sokolov A.P., Markevich M.N. Modelirovanie dielektricheskikh kharakteristik kompozitsionnykh materialov na osnove metoda asimptoticheskogo osredneniya [Modeling of dielectric properties of composite materials on the basis of asymptotic averaging]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2013, no. 1. DOI: 10.7463/0113.0531682
9. Emets Yu.P. Effektivnaya dielektricheskaya pronitsaemost' trekhkomponentnykh kompozitsionnykh materialov s anizotropnoy strukturoy [Effective permittivity of three-phase composite materials with an anisotropic structure]. Zhurnal tekhnicheskoy fiziki, 2005, vol. 75, no. 2, pp. 67-72. (English Translation: Technical Physics, 2005, vol. 50, iss. 2, pp. 207-212. DOI: 10.1134/1.1866436 ).
10. Shelukhin V.V., Terent'ev S.A. Gomogenizatsiya uravneniy Maksvella i dispersiya Maksvella-Vagnera [Homogenization of Maxwell equations and the Maxwell-Wagner dispersion]. Doklady Akademii nauk, 2009, vol. 424, no. 3, pp. 402-406. (English Translation:
Doklady Earth Sciences, 2009, vol. 424, iss. 1, pp. 155-159.
DOI: 10.1134/S1028334X09010334 ).
11. Bensoussan A., Lions J.L., Papanicalaou G. Asymptotic analysis for periodic structures. Amsterdam; New York: North-Holland Pub. Co., 1978. 396 p.
12. Kettunen H., Qi J., Wallen H., Sihvola A. Frequency dependence of effective permittivity of simple dielectric composites. Proceedings of ACES 2010, Tampere, Finland, 2010, pp. 248-253.
13. Sanchez-Palencia E. Book chapters. In: Sanchez-Palencia E., Zaoui A., eds. Homogenization techniques for composite media. Berlin, Heidelberg: Springer, 1987, pp. 121-192. DOI: 10.1007/3-540-17616-0 (Ser. Lecture notes in Physics; no. 272).
