CC BY
26Радиооптика
ХДК 517.1; 530.1
Оценки диэлектрической проницаемости сферопластика
Зарубин В. С.1'*, Кувыркин Г. Н.1, Савельева И. Ю.1 *fn2@bmstu.ru
1МГТУ им. Н.Э. Баумана, Москва, Россия
Ссылка на статью:
//Радиооптика. МГТУ им. Н.Э. Баумана. Электрон. журн. 2016. №3. С. 29-46.
Б01: 10.7463/^ор1.0316.0846170
Представлена в редакцию: 31.07.2016 © МГТУ им. Н.Э. Баумана
Двойственная вариационная формулировка задачи электростатики в неоднородной сплошной среде, содержащая два альтернативных функционала (минимизируемый и максимизируемый), достигающих на истинном решении задачи совпадающих экстремальных значений, использована в сочетании с четырехфазной моделью структуры сферопластика для построения гарантированных двусторонних границ области параметров, в которой расположены истинные значения эффективной диэлектрической проницаемости сферопластика. Путем перехода от исходной четырехфазной к более простой трехфазной модели структуры сферопластика получены расчетные зависимости для эквивалентной диэлектрической проницаемости микросфер, а затем и расчетные формулы для эффективной диэлектрической проницаемости сферопластика в целом. Количественный анализ полученных соотношений подтвердил адекватность построенных математических моделей. Эти соотношения могут быть применены для прогноза эффективной диэлектрической проницаемости сферопластика по известным характеристикам микросфер и связующего
Ключевые слова: сферопластик; микросфера; эквивалентная диэлектрическая проницаемость микросферы; эффективная диэлектрическая проницаемость сферопластика
Введение
Композиционные материалы с матрицей из полимерных и кремнийорганических связующих, армированные полыми шаровыми включениями (микросферами), принято называть сферопластиками. Распространенным материалом микросфер является стекло, но применяют и керамические, углеродные и полимерные микросферы, имеющие диаметр порядка миллиметра и толщину стенки в несколько микрометров [1, 2]. Малые плотность и теплопроводность и достаточно высокие механические свойства сферопластиков определили их широкое использование в качестве легких и прочных теплоизоляционных материалов. Но подбор материалов микросфер и связующего позволяет создавать сферопластики с необхо-димимыми свойствами применительно к различным областям техники [3, 4], в том числе с диэлектрическими характеристиками, которыми должны обладать конструкционные материалы разнообразных радиотехнических устройств и объектов [5, 6, 7].
Для создания сферопластиков с заданными диэлектрическими свойствами необходимо располагать математическими моделями, позволяющими описывать механическое и электрическое взаимодействие микросфер и связующего и строить расчетные зависимости для оценки диэлектрической проницаемости сферопластика, которую можно получить при использовании различных сочетаний материалов микросфер и связующего. Такие модели могут быть построены с привлечением основных положений механики и электродинамики [8, 9, 10, 11].
В данной работе построен вариант математической модели, представляющей сферо-пластик в виде сочетания четырех фаз: однородного материала с искомой эффективной диэлектрической проницаемостью и взимодействующего с представительным элементом структуры сферопластика, состоящим, в свою очередь, из шарового слоя связующего, полой микросферы и ее наполнения. Для гарантированной двусторонней оценки возможных значений эффективной диэлектрической проницаемости сферопластика к этой модели применена двойственная вариационная формулировка задачи электростатики в неоднородной сплошной среде, включающая два альтернативных функционала (минимизируемый и максимизируемый), достигающих на истинном решении задачи совпадающих экстремальных значений. Переходом от исходной четырехфазной модели структуры сферопластика к более простой трехфазной удалось вывести расчетные зависимости для эквивалентной диэлектрической проницаемости микросферы, а затем и расчетные формулы для эффективной диэлектрической проницаемости сферопластика в целом. Количественный анализ этих зависимостей подтвердил адекватность построенных математических моделей для описания механического и электрического взаимодействия элементов структуры сферопластика и достоверность определяемого полученными расчетными зависимостями прогноза эффективной диэлектрической проницаемости сферопластика.
1. Постановка задачи электростатики в неоднородной сплошной среде
Из полной системы уравнений Максвелла [9] вытекают уравнения электростатики в виде [12]
Ух Е = 0, V- Б = ре, (1)
где V — дифференциальный оператор Гамильтона; Е и Б — векторы напряженности электростатического поля и электрического смещения (электрической индукции) соответственно; 0 — нулевой вектор; ре — объемная плотность свободных электрических зарядов, которую в дальнейшем примем равной нулю. Первое уравнение (1) можно удовлетворить тождественно, если ввести соотношением
Е = -Уи (2)
скалярный электрический потенциал и. Для изотропной среды векторы Б и Е коллинеарны и связаны равенством [9]
Б = ££0Е, (3)
где е — относительная диэлектрическая проницаемость (для ваккума е = 1, а для диэлектриков е > 1), ео = 8,8542 ■ 10-12 ^м — электрическая постоянная. Далее слово
В
«относительная» для краткости будет опущено.
Изотропную сплошную среду в области V считаем неоднородной, т.е. диэлектрическая проницаемость е(М) этой среды зависит от положения точки М € V в этой области. Тогдаиз второго уравнения (1) и формул (2) и (3) при ре(М) = 0 (М € V) получим дифференциальное уравнение[12]
V ■ (е(М)Уи(М)) = 0. (4)
Для однозначного решения уравнения (4) в точках N € Б поверхности Б, ограничивающей область V, необходимо сформулировать граничные условия. На неконтактирующих между собой участках Б1 и Бо этой поверхности зададим значения
и(^ = и > 0, N € Б1; и(^ = 0, N € Бо (5)
электрического потенциала, а остальную часть Б* = Б \ (Б1 и Б0) поверхности области примем идеально электроизолированной, т.е.
е(^п^) ^и(^ = 0, N € Б*, (6)
где п^) — единичный вектор внешней нормали к поверхности Б* в точке N € Б*.
2. Двойственная вариационная формулировка задачи электростатики
Формулировке задачи электростатики в виде дифференциального уравнения (4) с граничными условиями (5) и (6) можно поставить в соответствие вариационную формулировку [8, 13], содержащую функционалы с определенными экстремальными свойствами. Сначала рассмотрим путь построения минимизируемого функционала [14].
Равенство (6) умножим на вариацию ¿и^) ^ € Б*) и проинтегрируем это произведение по идеально электроизолированному участку Б* поверхности Б, ограничивающей рассматриваемую область V, а произведение уравнения (4) и вариации ¿и(М) (М € V) проинтегрируем по этой области. Каждый из полученных интегралов равен нулю. Поэтому и их разность также будет равна нулю, т.е.
е(^п^) ■ VU¿и(^ ¿Б(^ - /(V ■ (е(М^и(М))) ¿и(М) ¿V(М) = 0.
Преобразуем второй интеграл в соответствии с первой формулой Грина [15] (обозначения точек опущены)
У (V ■ ^и^ ¿и ¿V = У е п ■ (VU) ¿и ¿Б - ^ еVU ■ V(¿U) ¿V. (7)
V 6 V
и учтем, что вариация би(Ж) = 0 (Ж е Би) на участках Би = Б1 и Б0 поверхности Б, поскольку на этих участках значения электрического потенциала заданы. Тогда в итоге получим равенство
б I (уи (м ))2 -V (М) = 0,
которое является условием стационарности функционала
3 [и] = ^ е(М )(уи (М ))2 -V (М). (8)
V
Функционал (8) допустимо рассматривать на множестве непрерывных в области V и кусочно дифференцируемых в ней функций и(М), М е V, принимающих на участках Бо и Б1 поверхности Б заданные значения электрического потенциала. Этот функционал на множестве допустимых функций является строго выпуклым (вниз) [8, 16], т.е. существует единственная функция и*(М), М е V, определенная в замкнутой области V = V и Б, которая является стационарной точкой функционала (8), соответствующей его наименьшему значению 3[и*]. В случае достаточной гладкости функции и*(М) она будет удовлетворять дифференциальному уравнению (4) и граничному условию (6).
Подынтегральная функция в соотношении (8) пропорциональна объемной плотности энергии электростатического поля в области V, а значение интеграла пропорционально суммарной энергии Ш условного электрического конденсатора с неоднородным диэлектриком в объеме области V, заключенным между двумя электродами на участках Б1 и Бо поверхности этой области с заданными электрическими потенциалами и1 и ио = 0 соответственно. Введя емкость С такого конденсатора, для стационарной точки функционала (8) можно записать
ш = у/ £(м )(уи *(м ))2 -V (М ) = Си?.
V
Для построения функционала, альтернативного по отношению к функционалу (8) и достигающего в стационарной точке наибольшего значения, равного 3[и*], расширим область определения функционала (8) путем введения векторной функции Ю, удовлетворяющей дополнительному условию
Ю(М)+ е(М(М) = 0, М е V. (9)
Отсюда получим Уи(М) = —Ю(М)/е(М) и после подстановки в функционал (8) запишем
3[и. 4 = 1 /(ЮМ?-V(М), (10)
а затем при помощи векторного множителя Лагранжа Ь введем условие (9) в соотношение
3[и, Ю, Ь] = 3[и, Ю] — I (Ю(М) + е(М)еоУи(М)) ■ Ь(М) -V(М). (11)
V
Равенство нулю вариации функционала (11) с учетом ее преобразования по первой формуле Грина вида (7) (обозначение точек опущено)
¿7[и, Б, Ь] = у ^(Б - Ь) ■ ¿Б - (Б + е^и) ■ ¿Ь + ееоУ ■ Ь¿^ ¿V - / (Ь ■ п)£и^Б = 0
приводит к условиям стационарности этого функционала, в которые помимо условия (9) входят равенства
Ь(М) = Б(М); V- Ь(М) = 0, М € V; ) ■ п(Ж) = 0, N е Б*. Используя эти равенства и соотношения (8), (10) и (11), приходим к функционалу
[и, Б] = 7[и] - ^ (Б(М) + е(М)^и(М))2 ¿V (12)
V
с дополнительными условиями (5) и условиями стационарности (9) и
V ■ Б(М) = 0 УМ е V; Б(^ ■ n(N) = 0, N е (13)
Из соотношения (12)) следует, что 7'[и, Б] ^ 7[и], а при выполнении условия стационарности (9) стационарные значения этих функционалов совпадают: 7'[и*, Б*] = 3[и*], причем
Б*(М) = -е(М)^и*(М), М е V. (14)
Если потребовать выполнения условий стационарности (13) в качестве дополнительных, то после преобразования функционала (12) с учетом функционала (8) и его дополнительного условия (5) получим
1 , (Б(М))2 и ,
1 [Б] = - 2 / ^ф^Т ^ (М) - Ц n(N) ■Б(N) (15)
V 4 У 0 0 51
Функционал (15) допустимо рассматривать на множестве непрерывно дифференцируемых в области V векторных функций Б(М), М € V, удовлетворяющих указанным выше дополнительным условиям, одно из которых совпадает со вторым уравнением (1), и интегрируемых на участках 5* поверхности Б. Этот функционал является строго выпуклым вверх [8, 16] и в стационарной точке Б*(М) (М е V и Б и) достигает наибольшего значения I [Б*], совпадающего со значением 7[и*].
Альтернативные функционалы (8) и (15) в сочетании с цепочкой неравенств
7[и] ^ 7[и*] = I[Б*] ^ I[Б] (16)
составляют двойственную вариационную формулировку задачи электростатики в неоднородном диэлектрике. Использование этой формулировки позволяет установить двусторонние границы, между которыми должны быть расположены эффективные значения характеристик такого диэлектрика, а также количественно оценить интегральную погрешность, возникающую при приближенном решении указанной задачи. В частности, соотношение (16) в сочетании с формулой (10) дает возможность получить двусторонние оценки емкости С условного конденсатора, соответствующего области V с неоднородным диэлектриком.
3. Представительный элемент структуры сферопластика
Структуру сферопластика представим как совокупность составных частиц шаровой формы с наружным радиусом Яу. Каждая из этих частиц включает шаровой слой связующего и микросферу с внутренним Я0 и наружным Я = Я^у3 радиусами, где Су- — объемная концентрация микросфер в сферопластике. Условно примем, что значение Яу может изменяться от некоторого конечного до бесконечно малого, что позволяет рассматривать сферопластик как сплошную среду и считать, что Су- е [0; 1].
Значение £° диэлектрической проницаемости для большинства используемых в сферо-пластиках связующих лежит в интервале 3,5 ... 5,0 [1, 6, 17, 18]. Для материалов, применяемых при изготовлении микросфер, значение электрической проницаемости £0 может изменяться в достаточно широких пределах. Например, для кварцевого стекла оно составляет 3,75, а для свинцового силиката, содержащего 20% РЬО — 16,2, причем £0 пропорционально плотности р в кг/м3 с коэффициентом около 2,4 [19]. С увеличением температуры диэлектрическая проницаемость стекла в интервале 293 ... 403 К возрастает лишь на 3 ... 10%, но с ростом температуры от 600 К до 1000 К удваивает свое значение [19, 20]. При кристаллизации стекол в силу упорядочивания структуры диэлектрическая проницаемость уменьшается (например, для боросиликатного стекла от 9,8 до 6,5).
Если в шаровой полости микросферы находится одноатомный или двухатомный газ, то его диэлектрическую проницаемость можно считать равной единице. В патентной литературе описаны способы получения микросфер на основе поливинилхлорида (£0 = 4) и полиметилметакрилата (£0 = 3,6), содержащих четыреххлористый углерод и другие низко-кипящие жидкости с диэлектрической проницаемостью £ = 2,26 — 2,78 [21].
Таким образом, модель структуры сферопластика включает четыре фазы: однородный материал с искомой эффективной диэлектрической проницаемостью £ сферопластика, связующее, материал полого шарового включения и наполнение микросферы с диэлектрической проницаемостью £;. Примем, что рассмотренная составная шаровая частица является представительным элементом структуры сферопластика. Основное свойство этого элемента состоит в предположении, что его замена равновеликим шаром с диэлектрической проницаемостью £ не вызывает возмущения электростатического поля в массиве однородного материала. Это допущение дает возможность использовать двойственную вариационную формулировку задачи электростатики в содержащей неоднородный материал замкнутой области для построения двусторонних оценок значения £.
4. Двусторонние оценки эффективной диэлектрической проницаемости сферопластика
Замкнутую область V представим в виде прямого цилиндра высотой Н ^ Яу с идеально теплоизолированной боковой поверхностью и заданными на каждом из оснований
площадью F ^ nf 2 значениями электрических потенциалов U0 = 0 и U > 0. Диаметральное поперечное сечение представительного элемента структуры сферопластика совместим с основанием цилиндра, имеющим нулевое значение электрического потенциала, т.е. область V будет содержать половину этого элемента. Остальная часть этой области будет занята однородным материалом с искомым значением е диэлектрической проницакмости. В качестве допустимого для функционала (8) примем линенйное по высоте цилиндра распределение электрического потенциала с постоянной составляющей G = U /Н вектора градиента. Тогда из формулы (8) получим
Ji[G] = eG*(HF - íf3) + „G- f3) + - f3) + . (17)
Для функционала (15) простейшим допустимым распределением вектора электрического смещения D является постоянное значение D = — ee0G единственной составляющей этого вектора, перпендикулярной основаниям цилиндра. В этом случае формула (15) примет вид
Ii [G] = —£G2( - — ^) — п ^ (+ ^ + R3) + eG2HF. (18)
Принятые допустимые распределения электрического потенциала и вектора электрического смещения для неоднородной области отличаются от действительных и поэтому значения Ji[G] и A[G] не будут совпадать, но будут удовлетворять соотношению (16), причем J1 [G] > I1 [G]. В промежутке между этими значениями должно быть расположено и значение J0 = eG2HF/2 функционала (8) для однородной области с диэлектрической проницаемостью е, поскольку замена половины представительного элемента структуры композита равновеликим объемом однородного материала не вызовет возмущения в принятом линейном распределении электрического потенциала, которое в этом случае будет соответствовать истинному распределению U*(M) (M £ V) в рассматриваемой замкнутой области. Тогда с учетом соотношения (17) из условия J1[G] ^ J0 получим для отношения ё = е/е° верхнюю оценку
е+ = ее° = 1 — Cv + cvе+ ^ ё, (19)
где = е+/е° = е0 + (ё — е0)С0, е0 = е0/е°, ё = е'/е° и C0 = (f0/f)3 — пористость
микросфер, а при использовании формулы (18) из условия 11[G] ^ J0 найдем нижнюю оценку
е_ 1
ё- = е- =1 — Cv + Cv/е- ^ё, (20) где . е* 1
- 1/ёс + (1/ё - 1/ёс)Со'
Отметим, что формулу для ё+ можно получить с применением теории смесей [22]. Однако из этой теории в отличие от использованной двойственной вариационной формулировки, содержащей функционалы (8) и (15), не следует, что ё+ является верхней оценкой значения ё.
Значения e+ и el можно рассматривать как соответственно верхнюю и нижнюю оценки некоторой эквивалентной диэлектрической проницаемости e* микросферы при ее замене равновеликим сплошным шаром радиусом R.
5. Эквивалентный коэффициент диэлектрической проницаемости микросферы
Рассмотрим взаимодействие отдельно взятой микросферы с внутренним R0 и наружным R радиусами и окружающего ее однородного материала с диэлектрической проницаемостью e° связующего сферопластика. Центр полого шара поместим в начале сферической системы координат. Примем, что на большом расстоянии r от начала координат задан модуль G вектора градиента электростатического поля в однородном материале. Этот вектор направим по оси сферической системы координат, от которой происходит отсчет угловой координаты 0. Тогда при r ^ то распределение электрического потенциала в этом материале будет описывать функция U^(r, 0) = U0 + Gr cos 0, которая удовлетворяет уравнению Лапласа, имеющего в сферических координатах вид [15]
1 д ( 2 dU\ 1 д / ^U\ 1 d2U _Q r2 dr V dr / r2 sin 0 д0 V д0 / r2 sin2 0 dp2 '
В данном случае благодаря параллельности заданного вектора градиента электростатического поля оси отсчета угловой координаты 0 распределение электрического потенциала симметрично относительно этой оси и не зависит от угловой координаты p, т.е. д2U^/dp2 = 0.
При приближении к микросфере возникает возмущение электростатического поля в однородном материале, описываемое также удовлетворяющим уравнению (21) дополнительным слагаемым AU(r, 0) = (B°/r2) cos 0 [23], где B° — подлежащий определению постоянный коэффициент. Таким образом, электростатическое поле в однородном материале, удовлетворяющее заданному условию при r ^ то и уравнению (21), описывает функция
/ B ° \
U°(r,0) = U<»(r,0) + AU(r,0) = Uo + ÍGr + —J cos0. (22)
Аналогичные зависимости описывают распределения электрического потенциала в шаровом слое микросферы
U (r, 0) = Uo + (Ar + B^j cos 0 (23)
и в шаровой полости
/ B' \
U'(r, 0) = Uo + í A'r + cos 0, (24)
причем в силу ограниченности значения электрического потенциала при г = 0 в соотношении (24) В' = 0.
В равенства (22)-(24) входят 4 неизвестных коэффициентов, которые необходимо найти из четырех граничных условий на сферических поверхностях с радиусами Л0 и Л, определяющих непрерывность на этих поверхностях распределения электрического потенциала
и радиальной составляющей вектора электрического смещения. Из граничных условий на сферической поверхности радиусом Ro с учетом равенств В' = 0, (23) и (24) получим
л, л В м ( л 2B N ёо
A' = А + , A' = А--т —.
До V R3 У ё'
Отсюда, исключая коэффициент A', находим
a = — —, (25)
R3 1 - А'
где ё' = ё'/ёо = е'/е0. Граничные условия на сферической поверхности радиусом R с учетом равенств (22) и (23) имеют вид
A + В = A + f- (A - !>о = А° - f. (26)
Заменим полую микросферу сплошным шаром радиусом R с эквивалентным значением ё* диэлектрической проницаемости, при котором коэффициенты A° и В° в формулах (26) сохраняют свои значения. При такой замене распределение электрического потенциала в сплошном шаре вместо равенства (23) будет определять зависимость U* = U0 + A*r cos 0, а
граничные условия (26) следует заменить равенствами
B° 2B°
A* = A° + —, А*ё* = A°---,
R3' R3 '
где ё* = е*/е°. Тогда с учетом равенств (26) запишем
А° + R= А° - £ = (А - RB)ё = (А + |)ё*.
Отсюда, используя формулу (25), для эквивалентной диэлектрической проницаемости микросферы находим
= е* = 2 + Г - 2(1 - ёрСО
ео 2 + ё' + (1 - ё')Со . ( )
В безразмерных параметрах ё+ и ё*, входящих в формулы (19) и (20), изменим нормировку, получив двусторонние оценки ё+ = 1 + (её' - 1)Со и её- = (1 + (1/ё - 1)Со)-1 для безразмерной величины ёё* эквивалентной диэлектрической проницаемости микросферы. На рис. 1 приведены графики зависимостей от Со верхней ё+ (штрихпунктирные линии) и нижней ё* (штриховые линии) оценок определяемой формулой (27) величины ё* (сплошные кривые) для значений ё', равных ординатам кривых при Со = 1. Видно, что все кривые для ёё* расположены в пределах границ, определяемых соответствующими двусторонними оценками. При Со = 0 и Со = 1 двусторонние оценки совпадают с соответствующим значением ёё**, но по мере отклонения параметра Со от этих значений разность Дё* = ё+ - ё* возрастает, достигая наибольшего значения (1 - \[ё')2 при Со = "\/fy(1 + Vf).
Величину в* = Дё*/(ё+ + ёё-) следует рассматривать как наибольшую возможную относительную погрешность, которая может возникнуть при выборе для микросферы в качестве ёё* полусуммы двусторонних оценок. На рис. 2 представлены графики зависимости в * от Со для некоторых значений ёё = 1 /ёё', причем эти графики симметричны относительно вертикальной прямой с абсциссой Со = 0,5.
Рис. 1. Двусторонние границы областей расположения значений эквивалентной диэлектрической проницаемости микросферы
Рис. 2. Зависимость наибольшей возможной относительной погрешности эквивален-той диэлектрической проницаемости микросферы от ее пористости при различных значениях параметра ё
6. Эффективный коэффициент диэлектрической проницаемости сферопластика
Проведенная выше замена полой микросферы сплошным шаром с эквивалентной диэлектрической проницаемостью е* позволяет от исходной четырехфазной модели структуры сферопластика перейти к трехфазной модели с представительным элементом в виде такого шара радиусом Я, покрытым шаровым слоем связующего с наружным радиусом . Если провести построение двусторонних оценок для трехфазной модели структуры сферопла-стика, то вместо формул (19) и (20) получим соответственно
ё+ = — = 1 + (е* - 1 )Су > ё, ё- = — = ----1——- ^ ё. (28)
+ е° 1 7 е° 1 + (1/е* - 1)Су 4 у
Эти оценки аналогичны двусторонним оценкам для величины ё* эквивалентной диэлектрической проницаемости микросферы. Поэтому для нахождения наибольшей возможной
относительной погрешности в = —I-, возникающей при выборе в качестве значения ё
+ е-
полусуммы двусторонних оценок, определяемых соотношениями (28), можно использовать рис. 2 с заменой в * на в, С0 на Су и ё на е*.
Для трехфазной модели известно решение задачи о тепловом взаимодействии аналогичного представительного элемента с неограниченным массивом однородного материала А [24], совпадающее по структуре с искомым решением для эффективного коэффициента диэлектрической проницаемости сферопластика и принимающее вид
ё = = (2 + е* - 2(1 - £*)Су)(2 + е* + (1 - £*)Су), (29)
аналогичный формуле (27). Это дает возможность после замен на рис. 1 £* на £, С0 на Су и £ на е* использовать данный рисунок для анализа влияния параметров Су и е* на величину Необходимо отметить, что формулу (29) можно преобразовать к виду, совпадающему с известной формулой Максвелла для среды с дисперсными шаровыми включениями, приведенной в работе [25] и полученной на основе двухфазной модели, состоящей из включения в виде сплошного шара с диэлектрической проницаемостью е* и окружающей его среды с диэлектрической проницаемостью е°, заполняющей неограниченную область.
Помимо соотношения (29) для вычисления эффективного значения диэлектрической проницаемости сферопластика можно использовать приближенную формулу [9, 26], получаемую осреднением свойств смеси из микросфер и частиц связующего,
е1 = е1/е° « (1 - Су + (е*)1/3)3, (30)
или соотношение
1п е2 = Су 1п £*, (31)
следующее из формулы Лихтенекера [6, 19]. Возможный подход к расчету эффективной диэлектрической проницаемости сферопластика связан также с применением метода самосогласования [27], приводящего к расчетной зависимости
£з = £3 = ^(1 - (1 - е*)(1 - 3Су) + ^(1 - (1 - ё*)(1 - 3Су))2 + 8£*). (32)
10
5
2 1
0,5
0,2 0,1
/>»/
8+, в—, 8,81,87,8^
•1, &2> &3
й-
/
/
у* у
/ ✓ 'А
\ ч
\ «...
4 4
\
1 ч^.- ¡^к
, \ \\ V
0
0,2 0,4 0,6 0,£
Су
Рис. 3. Двусторонние границы областей расположения значений эффективной диэлектрической проницаемости сферопластика, рассчитанных по формулам (29)-(32)
На рис. 3 в полулогарифмических координатах представлены построенные с использованием соотношений (28)-(32) графики зависимостей от Су двусторонних оценок £+
и е— (штрихпунктирные и штриховые линии соответственно), £, ¿Ti, £2 и £3 (соответственно сплошные кривые, пунктирные линии, сплошные кривые со светлыми квадратами и кружками) при различных значениях параметра £* (каждая кривая соответствует значению £*, равному ординате этой кривой при CV = 1). Все графики, определяемые расчетными зависимостями (29)-(32), не выходят за границы областей, определяемых соответствующими линиями для двусторонних оценок. Следовательно, математические модели, использованные для построения этих зависимостей, можно считать в первом приближении адекватно описывающими диэлектрические характеристики сферопластика.
Заключение
Двойственная вариационная формулировка задачи электростатики в неоднородной сплошной среде, содержащая два альтернативных функционала (минимизируемый и максимизируемый), достигающих на истинном решении задачи совпадающих экстремальных значений, позволила установить гарантированные двусторонние границы области параметров, в которой находятся истинные значения эффективной диэлектрической проницаемости сферопластика. Путем перехода от исходной четырехфазной модели структуры сферопла-стика к более простой трехфазной получены расчетные зависимости для эквивалентной диэлектрической проницаемости микросферы, а затем и расчетные формулы для эффективной диэлектрической проницаемости сферопластика, количественный анализ которых подтвердил адекватность построенных математических моделей. Полученные соотношения могут быть применены для прогноза эффективной диэлектрической проницаемости сферо-пластика по известным характеристикам микросфер и связующего.
Работа выполнена по гранту МК-6573.2015.8 программы Президента РФ государственной поддержки молодых кандидатов наук, а также в рамках проекта 1712 по государственному заданию № 2014/104 Минобрнауки РФ и государственного задания по проекту № 1.2640.2014.
Список литературы
1. Полимерные композиционные материалы: структура, свойства, технология / Под ред. A.A. Берлина. СПб.: Профессия, 2011. 560 с.
2. Наполнители для полимерных композиционных материалов: Справочное пособие / Пер. с англ. под ред. П.Г. Бабаевского. М.: Химия, 1981. 736 с.
3. Ушков С.С., Николаев Г.И., Михайлов В.И., Матвеев Г.В., ХесинЮ.Д. Конструкционные материалы для глубоководных аппаратов // Судостроение. 2004. № 5. С. 111-114.
4. Соколов И.И., Долматовский М.Г. Сферопластики // Полимерные материалы. 2005. № 9. С. 20-21.
5. Селиванов О.Г., Михайлов В.А. Теплоизоляционные синтактовые материалы на основе термостойкого кремнийорганического полимера // Международный журнал прикладных и фундаментальных исследований. 2014. №7. С. 12-13.
6. Михайлов В.А. Синтактные материалы с высокими диэлектрическими свойствами на основе кремнийорганического полимера // Успехи современного естествознания. 2015. № 12. С. 47-50.
7. Погосян М.А., Барковский А.Ф., Рожков А.И., Поляков Ю.Г., Господарский С.А. Антенный обтекатель, способ его изготовления и способ изготовления слоя антенного обтекателя: Патент РФ №2186444, 2002.
8. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 512 с.
9. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 8. Электродинамика сплошных сред. М.: Наука, 1992. 664 с.
10. Зарубин В.С. Моделирование. М.: Издательский центр «Академия», 2013. 336 с.
11. Зарубин В.С., Кувыркин Г.Н. Особенности математического моделирования технических устройств //Математическое моделирование и численные методы. 2014. № 1(1). С. 5-17. DOI: 10.18698/2309-3684-2014-1-517
12. Толмачев В.В., Головин А.М., Потапов В.С. Термодинамика и электродинамика сплошной среды. М.: Изд-во Моск. ун-та, 1988. 232 с.
13. Зарубин В.С., Селиванов В.В. Вариационные и численные методы механики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 1993. 360 с.
14. Зарубин В.С., Станкевич И.В. Расчет теплонапряженных конструкций. М.: Машиностроение, 2005. 352 с.
15. Власова Е.А., Зарубин В.С., Кувыркин Г.Н. Приближенные методы математической физики. 2-е изд. М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. 700 с.
16. Ванько В.И., Ермошина О.В., Кувыркин Г.Н. Вариационное исчисление и оптимальное управление. 2-е изд. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 488 с.
17. Медведев А.М. Печатные платы. Конструкции и материалы. М.: Техносфера, 2005. 304 с.
18. Чухланов В.Ю., Сысоев Э.П. Диэлектрические характеристики синтактных материалов на основе полых керамических микросфер с кремнийорганическими связующими // Стекло и керамика. 2004. № 6. С. 19-20.
19. Тареев Б.М. Физика диэлектрических материалов. М.: Энергоиздат, 1982. 320 с.
20. Применко В.И. Влияние состава на теплопроводность стекла / В сб. Вопросы химии и химической технологии. Вып. 62. Харьков: Вища школа. Изд-во при Харьк. ун-те, 1981. С. 72-74.
21. Физические величины: Справочник / Под ред. И.С. Григорьева, Е.З. Мейлихова. М.: Энергоатомиздат, 1991. 1232 с.
22. Головин Н.Н., Зарубин В.С., Кувыркин Г.Н. Смесевые модели механики композитов, Ч. 1. Термомеханика и термоупругость многокомпонентной смеси // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2009. №3. С. 36-49.
23. Зарубин В.С., Кувыркин Г.Н., Пугачев О.В. Оценки электрофизических характеристик композита с диэлектрической матрицей и дисперсными проводящими включениями // Радиооптика. МГТУ им. Н.Э.Баумана. Электрон. журн. 2015. №3. С. 51-67. DOI: 10.7463/rdopt.0315.0800066
24. Зарубин В.С., Кувыркин Г.Н., Савельева И.Ю. Эффективный коэффициент теплопроводности композита с шаровыми включениями // Тепловые процессы в технике. 2012. № 10. С. 470-474.
25. Maxwell C. Treatise on electricity and magnetism. In 2 vols. Oxford: Clarendon Press, 1873.
26. Зарубин В.С., Кувыркин Г.Н., Савельева И.Ю. Оценки диэлектрической проницаемости композита с дисперсными включениями // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2015. No3. С. 50-64. DOI: 10.18698/0236-3933-2015-3-50-64
27. Зарубин В.С., Кувыркин Г.Н., Савельева И.Ю. Оценка эффективной теплопроводности композита с шаровыми включениями методом самосогласования // Наука и образование. МГТУ им. Н.Э.Баумана. Электрон. журн. 2013. №9. С. 435-444. DOI: 10.7463/0913.0601512
Radiooptics of the Bauman MSTU, 2016, no. 3, pp. 29-46.
DOI: 10.7463/rdopt.0316.0846170
Received: 31.07.2016
© Bauman Moscow State Technical University
The Spheroplastic Permittivity Estimates
Zarubin V. S.1'*, Kuvyrkin G. N.1, Savel'eva I. Yu.1 * fti2@bmstu.ru
1 Bauman Moscow State Technical University, Russia
Keywords: spheroplastic, microsphere, the equivalent permittivity of the microsphere, effective permittivity of the spheroplastic
The composite materials with a matrix from polymer and silicone binders, reinforced by the hollow spherical inclusions (microspheres) are commonly referred to as spheroplastics. A widespred microsphere material is glass, but, it is also possible to use ceramic, carbon, and polymeric microspheres having a diameter within a millimeter and a wall thickness of several micrometers. Low density and thermal conductivity and sufficiently high mechanical properties of the spheroplastics decided their widespread use as the light and strong thermal insulation materials. But a selection of the microsphere and binder materials allows us to create the spheroplastics with appropriate properties to be used in variety of technology domains, including those with dielectric properties, which the structural materials of various wireless devices and objects should have.
To create the spheroplastic with specified dielectric properties it is necessary to have a mathematical model allowing us to describe the mechanical and electrical interaction of the microspheres and binder and build the calculated relationships to estimate a spheroplastic permittivity, which can be obtained when using the various combinations of microsphere materials and binder. Such models can be built involving the essentials of mechanics and electrodynamics
The paper offers a constructed mathematical model to represent the spheroplastic as a combination of four phases: a homogeneous material with the desired effective permittivity and a spheroplastic interacting with the representative structure element consisting, in turn, of a ball bonding layer, a hollow microsphere and and its filling. To ensure a two-sided estimate of the possible values of spheroplastic effective permittivity this model uses a dual variational formulation of electrostatic problems in nonhomogeneous continuum, including two alternative functionals (minimized and maximized), reaching the matching extreme values to true solution. The transition from the initial four-phase model of the spheroplastic structure to the simple three-phase one enabled us to derive the calculated relationships for the equivalent permittivity of microspheres, and then also the calculating formulas for effective permittivity of the spheroplastic as a whole. A quantitative analysis of these relationships proved the adequacy of the constructed mathematical
Radiooptics
Electronic journal of the Bauman MSTU
http://radiooptics.ru
models to describe the interaction between mechanical and electrical elements of the structure and the predicted accuracy of the effective spheroplastic permittivity determined owing to calculated dependences.
References
1. Berlin A.A., ed. Polimernye kompozitsionnye materialy: struktura, svoistva, tekhnologiia [Polymeric composites: structure, properties, technology]. Sankt-Peterburg, Professiia Publ., 2011. 560 p. (in Russian).
2. Babayevsky P.G., ed. Napolniteli dlya polimernykh kompozitsionnykh materialov [Fillers for polymeric composite materials]. Moscow, KhimiiaPubl., 1981. 736 p. (in Russian).
3. Ushkov S.S., Nikolaev G.I., Mikhailov V.I., Matveev G.V., Khesin Yu.D. Structural materials for deep-water submersibles. Sudostroenie = Shipbuilding, 2004, no. 5, pp. 111-114. (in Russian).
4. Sokolov I.I., Dolmatovskii M.G. Spheroplastics. Polimernye materialy = Polymer materials, 2005, no. 9, pp. 20-21. (in Russian).
5. Selivanov O.G., Mikhailov V.A. Insulation sintakt materials based on thermal stability of silicone polymers. Mezhdunarodnyi zhurnal prikladnykh i fundamental'nykh issledovanii = International Journal of Applied and Fundamental Research, 2014, no. 7, pp. 12-13. (in Russian).
6. Mikhailov V.A. Syntactic materials with high dielectric properties based on silicone polymer. Uspekhi sovremennogo estestvoznaniia = Advances in current natural sciences, 2015, no. 12, pp. 47-50. (in Russian).
7. Poghosyan M.A., Barkovskii A.F., Rozhkov A.I., Polyakov Yu.G., Gospodarskii S.A. Anten-nyy obtekatel', sposob ego izgotovleniya i sposob izgotovleniya sloya antennogo obtekatelya [Antenna radome: its manufacturing method and a method for manufacturing antenna radome layer]. Patent RF no. 2186444, 2002.
8. Zarubin V.S., Kuvyrkin G.N. Matematicheskie modeli mekhaniki i eleknhjdinamiki sploshnoi sredy [Mathematical models of continuum mechanics and electrodynamics]. Moscow, Bauman MSTUPubl., 2008. 512 p. (in Russian).
9. Landau L.D., Lifshitz E.M. Teoreticheskaya fizika. T.8. Electrodynamika sploshnykh sred [Theoretical physics. Vol. 8. Electrodynamics of Continuous Media]. Moscow, Nauka Publ., 1992. 664 p. (in Russian).
10. Zarubin V.S.Modelirovanie [Modeling]. Moscow, "Academia" Publishing Center, 2013. 336 p. (in Russian).
11. Zarubin V. S., Kuvyrkin G.N. Features of mathematical modeling of technical devices. Matem-aticheskoe modelirovanie i chislennye metody = Mathematical simulation and numerical methods, 2014, no. 1(1), pp. 5-17. DOI: 10.18698/2309-3684-2014-1-517 (in Russian).
12. Tolmachev V.V., Golovin A.M., Potapov V.S. Termodinamica i elektrodinamicasploshnoi sredy [Thermodynamics and electrodynamics continuous medium]. Moscow, MGU Publ., 1988. 232 p. (in Russian).
13. Zarubin V.S., Selivanov V.V. Variatsionnye i chislennye metody mekhaniki sploshnoi sredy [Variational and numerical methods of continuum mechanics]. Moscow, Bauman MSTU Publ., 1993. 360 p. (in Russian).
14. Zarubin V.S., Stankevich I.V. Raschet teplonapryazhennykh konstruktsiy [Calculation of heat-stressed designs]. Moscow, Mashinostroenie Publ., 2005. 352 p. (in Russian).
15. Vlasova E.A., Zarubin V.S., Kuvyrkin G.N. Priblizhennye metody matematicheskoi fiziki [Approximate methods of mathematical physics]. Moscow, Bauman MSTU Publ., 2004. 700 p. (in Russian).
16. Van'ko V.I., Ermoshina O.V., Kuvyrkin G.N. Variatsionnoe ischislenie i optimal'noe upravle-nie [Calculus of variations and optimal control]. Moscow, Bauman MSTU Publ., 2001. 488 p. (in Russian).
17. Medvedev A.M. Pechatnye platy. Konstruktsii i materialy [Printed circuit boards. Structures and materials]. Moscow, Tekhnosfera Publ., 2005. 304 p. (in Russian).
18. Chukhlanov V.Yu., Sysoev E.P. Dielectric characteristics of syntact materials based on hollow ceramic microspheres with organosilicon binders. Steklo i keramika, 2004, no. 6, pp. 19-20. (English version of journal: Glass and ceramics, 2004, vol.61, no. 5, pp. 190-191. DOI: 10.1023/B:GLAC.0000043090.07100.67).
19. Tareev B.M. Fizika dielektricheskikh materialov [Physics of dielectric materials]. Moscow, EnergoizdatPubl., 1982. 320 p. (in Russian).
20. Primenko V.I. Influence of structure on the thermal conductivity of glass. Voprosy khimii I khimicheskoi tekhnologii [Chemistry and chemical technology questions], 1981, vol.62, pp. 72-74. (in Russian).
21. Grigoryev I.S., Meilikhov E.Z., eds. Fizicheskie velichiny [Physical quantities]. Moscow, EnergoatomisdatPubl., 1991. 1232 p. (in Russian).
22. Golovin N.N., Zarubin V.S., Kuvyrkin G.N. Mixture Models of Composite Mechanics. P. 1. Thermal Mechanics and Thermoelasticity of Multicomponent Mixture. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki = Herald of the Bauman Moscow State Technical University. Ser. Natural sciences, 2009, no. 3, pp. 36-49. (in Russian).
23. Zarubin V.S., Kuvyrkin G.N., Pugachev O.V. Estimating Electro-physical Properties of the Composite with Dielectric Matrix and Dispersed Conductive Inclusions. Radiooptika. MGTU im. N.E.Baumana = Radiooptics of the Bauman MSTU, 2015, no. 3, pp. 51-67. DOI: 10.7463/rdopt.0315.0800066 (in Russian).
24. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Iu. The Effective Coefficients of Thermal Conductivity of Composites with Spherical Inclusions. Teplovye protsessy v tekhnike = Thermal Processes in Engineering, 2012, no. 10, pp. 470-474. (in Russian).
25. Maxwell C. Treatise on electricity and magnetism. In 2 vols. Oxford, Clarendon Press, 1873.
26. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Evaluation of Dielectric Permittivity of Composite with Dispersed Inclusions. Vestnik MGTU im. N.E. Baumana. Ser. Priborostroenie = Herald of the Bauman Moscow State Technical University. Ser. Instrument Engineering, 2015, no. 3, pp. 50-64. DOI: 10.18698/0236-3933-2015-3-50-64 (in Russian).
27. Zarubin V.S., Kuvyrkin G.N., Savel'eva I.Yu. Evaluation of composite effective thermal conductivity with spherical inclusions by the self-consistency method. 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).
