CC BY
115
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Моделирование диэлектрических характеристик композиционных материалов на основе метода асимптотического осреднения # 01, январь 2013 Б01: 10.7463/0113.0531682
Димитриенко Ю. И., Соколов А. П., Маркевич М. Н.
УДК 539.8
Россия, МГТУ им. Н.Э. Баумана [email protected] al sokol o@bmstu. ш [email protected]
1. Введение
Проектирование композиционных материалов с заданными электромагнитными свойствами является важной технической проблемой. Данные материалы находят широкое применение в различных областях техники, в частности, используются в качестве современных пьезоэлементов, электроизоляционных конструкций [1, 2], радиопрозрачных и радиопоглощающих конструкций. Для математического моделирования характеристик композиционных материалов используюися различные приближенно-аналитические и численные методы. Наиболее перспективным является метод асимптотического осреднения (МАО) или «метод гомогенизации», основы которого разработаны в работах Э. Санчес-Паленсии [3], Бахвалова Н.С. [4], Победри Б.Е. [5] и других. Данный метод позволяет математически точно вычислять эффективные характеристики композитов с помощью решения специальных «локальных задач» на ячейках периодичности. Однако, аналитическое решение такого рода задач для композитов, обладающих сложной геометрической структурой, невозможно, а численное осложняется интегро-дифференциальной постановкой локальных задач с неклассическими краевыми условиями периодичности. Методика преобразования локальных задач к классическим краевым задачам с граничными условиями первого и второго рода предложена в работах [6-10]. В настоящей работе данная методика применяется для решения локальных задач электростатики, на основе которой расчитываются эффектиные тензоры диэлектической проницаемости композитов со сложными структурами армирования.
2. Метод асимптотического осреднения для задачи электростатики
Рассмотрим композиционный материал, занимающий в пространстве Я3 область V с поверхностью Рассматриваемая область V состоит из N фаз: Va,a = 1...N-1-наполнители (например, волокна), VN - матрица. Введем обозначения: Еа - поверхности областей Уа, и Е^ - поверхности контакта матрицы и волокон, волокна полагаем не контактирующими между собой, Еае - часть поверхности Е композита, занятая а -ым
n-1
компонентом (причем Еа = Е^ иЕае - для волокон и ЕN = ЕаN иЕш- для матрицы).
а=1
Компоненты композита полагаем изотропными. В каждой области Va,a = 1..^
рассмотрим следующую задачу электростатики, являющуюся следствием уравнений Максвелла [12] для установившихся электромагнитных колебаний:
д Дг а = 0, х,. е Va
Д а, X е Va и ^
Е*а=д1р*а, , х е Va (1)
ра=Р,(Да- Д7 )п* = 0, х* еXaV
*а | * т-1 *а | тл*а v"*
Р 1 = P*, Дг пг 2 = Д , х* е ^
где Di - компоненты комплексной амплитуды вектора электрической индукции, -компоненты комплексной амплитуды вектора напряженности электрического поля, р* -
*а
комплексная амплитуда электрического потенциала, е - комплексная амплитуда диэлектрической проницаемости. Символ * обозначает, что рассматриваются амплитуды комплексных величин (например, комплексная амплитуда электрического потенциала
*а ^ а • п а • \
р = р + гр , где г - мнимая единица), соответствующие электромагнитным колебаниям с известной частотой с. Все компоненты векторов отнесены к прямоугольной системе координат, дг - дифференцирование по 1-ой декартовой
координате х .
Для применения МАО введем предположение о том, что композиционный материал обладает периодической структурой, область ячейки периодичности (ЯП) ^ которого
состоит из фаз V^a,a = 1, N, матрица является связанной областью. Поверхности раздела компонентов композита в ЯП обозначим через ^^ .Тогда можно ввести малый
параметр к =1 << 1, и два типа безразмерных координат: хг - глобальные и ^ -
локальные, которые определяются через декартовые координаты Хг следующими соотношениями
— Хг £: Х1 Х1
Х = ~Т , Л=~ = -7-, Ь к I
где I - длина ребра ячейки периодичности, а Ь - характерный размер всего композита. Согласно общей концепции МАО функции / = {Д*, Ё*,ю*} для композита рассматриваются как квазипериодические: /(Хг) = /(Х1+ а,), Х1 е VЛ, е V^ , где а, - компонеты целочисленного вектора. Эти функции дифференцируются по правилу:
дю* * 1 * * дю* * дю* _, ,1Ч
--> ю ,г+—ю /г, где ю,г = , Ю/г =-• Тогда решение задачи (1) строится в виде
дХ К дХ1
асимптотического разложения:
ад
юа (Х , л) = юмо) (Х) + кюМ1) (Хг , Л) + куа(2) (Х Л,) + ••• = юа(0) (Х) + 2кп ю*а( п) (X, Л, )
п=1
1
Ё*а = юа(0),. +кю,*а(1) +к —ю*"(1)/г +к2ю,а(2) +••• = ^к" (ю*а("),1 +юМ"+1)/г) = 2кпЁ*а(п), (2)
К п=0 п=0
д*а = 2 кпд*а(п) д*а(п) = £*аё*а(") п > 0
2 ^^^ I ~> i I '
Векторы электической индукции и напряженности электрического поля «нулевого уровня» имеют вид
д *а(0) = £*а Ё *®(0)
г*а(0) _ тр* *а(1) г ~ ^г ^ У/,
-г ' V/г ~* _ л,*«(0)
(3)
Ёг =ю*
Подставив (2) и (3) в задачу (1), получим следующую систему уравнений в нулевом приближении:
' ) а(0) = 0, г/ г >
д *а(0) = Ё *«(0)
Х г е
Ё*а(0) = Ё* +ю*а(1)
ю
/
а(1) =ю*м(1), (Да(0) - (0) )п = 0, Х е 2, юа(1) >= 0,[[= 0,[[юМ1)]] = 0, Хг е 2
Хг е
Хг е
(4)
ЛаЫ
здесь Д г ',ю - компоненты вектора электрической индукции и электрический потенциал фаз композита Р^, а = N в рамках одной ЯП; Ё* =ю*а(0)- компоненты
эффективного вектора напряженности электрического поля; < р*(1) >= 0-условие нормировки, необходимое для единственности решения в классе периодических функций,
N . а
где < р*а(1) >= £ Г р*а(1)dV^ - операция осреднения по «ячейке периодичности»;
а=1
[[р ()]] = 0, [[Д ]]пг = 0- условия периодичности для электрического потенциала р и вектора электрической индукции Д*(0). В силу периодичности функции р*(1) имеет
7^"* 7—' * а (0) *а (0)
место соотношение: Е* =< Е* ^ ' >= р { '.
3. Преобразование локальной задачи к задачам «классического типа»
Решение задачи (4) будем искать в виде сумм
ра (1) = £ р*р) , (5)
р=1
где ррР) - функции следующего вида
р(Р) = -Ё^Р +$Р)({к ), (6)
здесь 5(*р)(£к )- новые неизвестные функции от , уже не являющиеся периодическими. Производные от функций (6) по локальным координатам имеют вид
Р* а (1)/г = £ р^П = -I Е;51р + £ = - Е* + £ ^)/г . (7)
Р=1 Р=1 р=1 р=1
Вычислим напряженность через потенциал, используя (7):
3 3
Ег*а(0) = Е* +р/?(1) = Е* - Е* + £%р)/г =£$
(р)/1 •
Р=1 Р=1
Отсюда получим, что напряженность Е * (0) , подобно потенциалу, представима в виде сумм
3
Е *а(0) = v Е*а г ~ Аа^Кр)' р=1
где функции Е*"), называемые псевдонапряженностями, имеют вид
Е*аР) = р«/г . (8)
Электрическую индукцию Д*а (0) можно так же представить в виде сумм
3
д *а(0) = д *а Дг - АиДг(р),
Р=1
где псевдоиндукции связаны с псевдонапряженностями линейными соотношениями
л* а = *а Е* а Л(Р)~Б е1(Р).
Д * а =*а Е * а (9)
г ( р) ~ь р)• v)
В силу линейности задачи псевдоиндукция будет удовлетворять соотношениям:
ДО)/, = 0, Хг е ^
Граничные условия для функций 3^ (Л ) :
п*а _ п*N
3(Р) = 3(Р) ■
(Дй - =о, Хг
Также должны выполняться условия нормировки:
(11)
<юа >= 2<ю*р) > = -2ё(*р) <Лр > +2<3(*Р) > = 0^
Р=1 р=1 р=1
Так как <ЛР > = 0, то функции 3(*Р)(Л) должны удовлетворять условиям нормировки:
<3(р,(4 )>=0-
Из условия периодичности псевдопотенциала
(12)
[[юа(1)]] = 2[[юа(р)]] = -2 ё(*р)[[Л ]]+2 [[О = 0^
Р=1 р=1 р=1
получаем условия для 3(р):
ко=Ё(*р)[Лр ]]• (13)
Так как величина [ЛР ]] удовлетворяет следующему соотношению: [ЛР ]]а = ар8ра, то условия (13) для функций $(*Р)(Л ) примут вид
[[3(*"]]. = Ё* )а 8
^ (Р)л -|г г(р) Р I
г(р) р гр
(14)
т.е функции 3(*Р) (Л ) уже не являются периодическими.
Из уравнений (8)-(12), (14) следует, что функции 3*0)(Л) являются решением следующих задач на «ячейке периодичности»
'Д 0 = 0
г(р)/г и
(р) = 8 Ёг(р)
тр*а _ п*а Ёг( ) ( )/г
п*а n*N , -т\*а Гл*N \ п.
3(р) = 3(р),(а(р) - а(р))пг = 0
'(р) ~ "(р)'^г(р)
<3(*Р) >= 0
(15)
[[3(р)]]г = Ё(*Р) а р81р [[[ Д*сР)]]п1 = 0
Задача (15) в отличие от задачи (4) не содержат входных данных в соотношения Гаусса, а имеют входные данные - функции Ё*, заданные на поверхности «ячейки
периодичности», но они по-прежнему содержат интегральное условие нормировки и граничные условия периодического типа. Подобные задачи по расчету электрических
полей внутри композита разрешимы на ЯП с достаточно простой геометрической структурой [8]. Для более сложных структур, обладающих свойством симметрии, введем модификацию МАО.
Далее предположим, что в ЯП V¡, является симметричной относительно координатных плоскостей (^ = 0) ,относительно поворота на угол п вокруг каждой оси координат О£ и при преобразовании центральной симметрии с центром в точке О. Тогда задаче (15) можно в соответствие поставить задачу 1/8 ЯП (в первом координатном октанте, составляющем подобласть : = п (^ > 0).
А а = 0
А") = ^ Е^ Хг е и£, )
Е:0Р) =«Хг е V
€) = 0 ДКР) - ))Пг = 0, Хг av
(16)
Поверхности контакта компонентов £ = £ п , координатные плоскости -
£ = = 0}, а торцевые поверхности ЯП - £ = = а!, /2}, ^ = 1,2,3. Граничные условия на торцевых и координатных плоскостях различны для всех р:
а _
п *а | =_РЕ * 1*а | = 0 *(Р) = „ ЕР Р) ^ = 0,
£р
(17)
С/г1£ = 0,^«/г и = 0, Р * г
2 ^ (Р) £ Р £
Задачи (16) на 1/8 «ячейке периодичности» У^ = ^ с граничными условиями
а=1
(17) будем называть задачам Ьр, р = 1,2,3.
4. Расчет эффективного тензора диэлектрической проницаемости композита
После решения серии задач Ьр дляр=1, 2, 3 и нахождения псевдопотенциалов 3*р) и
7~ч* а
псевдоиндукции Д р) во всех компонентах композита, вычислим среднюю
электрическую индукцию, используя интегрирование по областям, занятым компонентами композиционного материала:
< Д*(0) >=£ А
г (р)
р=1
где
— N
ДР) =< Д0р) >=2 I Дг*Р)dV • (19)
а=1 V ( а )
В силу линейности задач ЬР их решения Д*; линейно зависят от входных данных Ё*р • Следовательно, существуют тензоры 8*(р>), связывающие комплексные амплитуды
7—' * 7~\ * а
вектора напряженности Ёр и вектора индукции электрического поля Д р):
Д0р)=8Г Ё; • (20)
Следовательно, зная псевдоиндукции Д *0), можем вычислить компоненты тензора диэлектрической проницаемости 8г*Р(0) нулевого приближения:
Д *а Л )
*(0)
г; Ё * '
Тогда, подставив (20) в (19), получаем осредненную задачу с эффективным тензором диэлектрической проницаемости 8Р;, связывающим среднюю индукцию < Дг > и
среднюю напряженность Ё*р электрического поля внутри композиционного материала:
< Д*(0) >=8 *Ё* ,
г гр р >
где эффективный тензор комплексных амплитуд диэлектрической проницаемости 8*; вычисляется по формуле
Д
г; Ё* ' ЁР
5. Вариационная формулировка задач электростатики Ьр
Для произвольного объема V е дадим вариационную формулировку задач Ь; (16)
и (17). Для этого рассмотрим класс комплекснозначных функций Ю(*Р) (возможные
значения электрического потенциала), которые определены во всей области V и
являются гладкими в подобластях V0л, удовлетворяют условиям непрерывности на границе раздела фаз и граничным условиям первого рода на торцевых поверхностях
^^ р и ^^; :
ар^ *
[ю*р )]=о, юр) 2=о, Ю(*р)!2: (21)
Поверхность 2^ области подобласти VЛ представляем как совокупность
поверхностей: 2 = 2 ^ 2г . На поверхности 2 задано внешнее поле Д*
Д(р) = Д**р)п , а на поверхности - нулевые граничные условия. Обозначим также
3(р) - вариации возможных значений электрического потенциала, удовлетворяющие тем же условиям (21), но с нулевыми граничными условиями на торцевых поверхностях X' . Истинный электрический потенциал 3*р), удовлетворяющий всем уравнениям задачи (16)-(17), отличается от всех возможных значений потенциала (р*р) тем, что для него и только для него лагранжиан Ь имеет минимальное значение:
3Ь = 0,
Ь =| (р)Д(Р)^ - 21 ((р)/г УУ^р)/¿V.
(22)
I,
6. Метод конечных элементов для задач Ьр
Для решения вариационного уравнения (22) применим метод конечных элементов
(МКЭ), согласно которому всю область разбиваем на конечное, достаточно большое количество подобластей стандартного типа Ve в форме тетраэдров: Ve . Для
е
каждого элемента Ve записываем уравнение (22), при этом правая часть вариационного уравнения аппроксимирует поток электрической индукции Д*) = е*3*(р)]1через
поверхность отдельно взятого конечного элемента:
¡3(3^,гУеХ,/^ = | ¿Ц*р))т
(23)
В каждом конечном элементе аппроксимируем псевдопотенциал 3(р )£ г' ( 1 ) ) = 3р)£ г(1 )),3(р)(^( 1))} и его вариацию 33(р) как
3 = [Ф]{ф*р)Ь 3) = ms{Фlp)},
где {ф(*р)}т= {ф'р )(£ г (1)), ф('р) (£(1) Х-, ф('р) (т)), ффр) (т))}Т координатный столбец, составленный из значений электрического потенциала в узлах конечного элемента, ) -координата 1-го узла конечного элемента,1=1...т, т- число узлов, [Ф]2х8 - матрица функций формы, имеющая вид
[Ф] =
Производные от псевдопотенциала представим в виде координатного столбца {3(*р)/1,3(*р)/2,3(*р)/3} = [В]ф(*р), где [В] = [Ь][Ф]- матрица градиентов потенциала,
Ф1 0 Ф 2 0 Ф 3 0 Ф 4 0
0 Ф1 0 ф2 0 ф3 0 Ф 4 _
[¿]т =
д 0 д 0 д 0
дЛ дЛ дЛ
0 д 0 д 0 д
дЛ дЛ дЛ
- матрица операторов дифференцирования. С учетом введенных операторов, перепишем уравнение (23):
\ЩР)}т[В]т8*[Б]{Ф1 )}ё¥ = | 5{Ф1 )}т[Ф]тД¿И.. (24)
Вынося вариацию 5{ффр)}т за знак интегралов, из (24) получаем разрешающую систему линейных алгебраических уравнений (СЛАУ):
[<]{Ф(*р)} = {/.*}, (25)
где [Ке*]8х8 =|[В]та*[Б]ё¥ - матрица жесткости, {/*} = | [Ф]тВ*рёЪ- локальный вектор
Ее
правой части. Если просуммировать СЛАУ (25) для всех конечных элементов, то получим глобальную СЛАУ для всей рассматриваемой области VЛ
[К *]{Ф(*р)} = {/•}, (26)
где [К *]- глобальная матрица жесткости, {/ *}- глобальный столбец правой части.
7. Расчет для 3D ортогонально-армированного композиционного материала
Решение глобальной СЛАУ (26) осуществлялось путем разделения комплексных переменных на действительную и мнимую части с последующим привлечением QMR методов. С помощью такого подхода были решены Ьр задачи для 3Б ортогонально-армированного композита [6,10], число составляющих элементов которого N=4, где а = 1,2,3-волокна, а = 4 - матрица. Все волокна полагались одинаковыми и концентрация волкон по 3-м координатным направляниям также была одинакова. Такой композит
вцелом является материалом с кубическим типом симметрии [13], и тензор эффективной
_* _* _*
диэлектической проницаемости а рр имеет одну независимую компоненту а = а и
(суммирования здесь нет).
В численных расчетах число конечных элементов при решении Ьр задач составляло 6343, а число степеней свободы - 1502. При численной реализации учитывалось, что армирующие композита могут состоять из различных материалов и иметь различные радиусы.
т
Расчет действительной и мнимой частей электрического псевдопотенциала для задачи L3 (рис.1 и 2) производился в композите с коэффициентом армирования (pf = Vf /Vm =0.18, где Vf ,Vm- содержание волокон и матрицы соответственно.
Модельные значения коэффициентов диэлектрической проницаемости
Res Ims Ims f „ . удовлетворяют соотношениям: -— = 0.1,-— = 0.1,-— = 0.1.
Res f
Res
Res f
Рис.1. Распределение действительной части псевдопотенциала Re5* (безразмерное
значение) в ЯП.
Приближенно-аналитическое значение эффективного коэффициента диэлектрической проницаемости по методу Фойгта (линейная зависимость) с помощью следующих соотношений
а* = е( +а*т(1 ). (27)
По методу Рейса (обратно-линейная зависимость) данное значение ивычисляется следующим образом:
= ^ (1 -Pf)
_* * + * .
s s* sm
(28)
На рис. 3 изображена зависимость действительной части коэффициента диэлектрической проницаемости композита Res*, отнесенной к действительной части коэффициента проницаемости матрицы Res*m, от коэффициента армирования cpf, на рис.
4 - мнимой части Ims*, отнесенной к Res*m, а на рис. 5 - абсолютной величины компоненты модуля s* диэлектрической проницаемости композита.
m
Результаты численных расчетов показывают, что значения Res*,Ims*, s* для 3D
ортогонально-армированного композиционного материала укладываются в вилку Фойгта-Рейсса, что свидетельствует о хорошей точности предложенного метода расчета эффективных характеристик.
Рис. 2. Распределение мнимой части псевдопотенциала 1т 3* (безразмерное значение) в
ЯП.
В то же время сами соотношения (27)-(28) нельзя выбрать в качестве даже приближенных выражений для эффективного коэффициента диэлектрической проницаемости в виду того, что вилка Фойгта-Рейсса слишком широка. Различия между значениями |а*|, вычисленными по предложенному методу и по методу Фойгта-Рейсса
составляло примерно 60% для коэффициента армирования ср^- = 0.55, что является весьма плохим результатом.
9t
Рис. 3. Зависимость Res* от коэффициента для 3D композита, рассчитанная по методу
Фойгта(Е)-Рейсса^) и МАО (точки).
0.0 -г г г I I Г
0.0 0.1 0.2 0 3 0 4 0.5 О.в 0.7 0.8
Рис. 4. Зависимость 1т а * от коэффициента для 3Б композита, рассчитанная по методу
Фойгта(Р)-Рейсса^) и МАО (точки)
Таким образом, для расчета эффективных электрических характеристик композитов со сложными структурами армирования целесообразно приемнять метод асимптотического осреднения, рассмотренный в данной работе. Этот метод позволяется вычислять математически точные значения эффективных диэлектричских характеристик композиционных материалов. Возможные погрешности метода могут быть связаны
только с погрешностями численного метогда расчета, которые достаточно малы (не более 1%) и могут быть еще уменьшены за счет выбора более мелких конечно-элементных сеток, а также с погрешностями реальной геометрической структуры волокон, которая может отличаться от идеальной формы, использованной в расчетах. Однако и эти погрешности могут быть снижены, за счет более точного учета геометрической формы волокон в численном расчете.
Рис. 5. Зависимость | е* | от коэффициента для 3Б композита, рассчитанная по методу
Фойгта(Б)-Рейсса(К) и МАО (точки)
Заключение
Предложена математическая модель для расчета эффективного тензора диэлектрической проницаемости композиционных материалов на основе метода асимптотического осреднения периодических структур. Для численного решения локальной задачи электростатики использован метод конечных элементов. Проведены тестовые расчеты эффективных диэлектрических характеристик 3Б ортогонально-армированного композита. Сравнение результатов, полученных при численной реализации метода асимптотического осреднения и метода Фойгта-Рейсса, показало, что предложенная модель расчета обеспечивает высокую точность вычислений эффективных диэлектрических характеристик композиционных материалов со сложными структурами армирования. Данная математическая модель может быть применена для прогнозирования диэлектрических характеристик новых синтезируемых материалов.
Поддержка
Работа проведена при финансовой поддержке Министерства образования и науки Российской Федерации по Соглашениям о предоставлении гранта в форме субсидий №14.В37.21.0448 и №14.132.21.1699 в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на
2009-2013 годы.
Список литературы
1. Бычков И.В., Дубровских Д.В., Зотов И.С., Федий А.А. Исследование эффективной диэлектрической проницаемости композитного материала Са804-2Н20 - графит // Вестник Челябинского государственного университета. 2011. № 7 (222). Физика. Вып. 9. С. 7-15.
2. Най Дж. Физические свойства кристаллов : пер. с англ. М.: Мир, 1967. 385 с.
3. Санчес-Паленсия Э. Неоднородные среды и теория колебаний: пер. с англ. М.: Мир, 1984. 472 с.
4. Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. М.: Наука,1984. 352 с.
5. Победря Б.Е. Механика композиционных материалов. М.: Изд-во МГУ, 1984. 336 с.
6. Димитриенко Ю.И. Кашкаров А.И. Расчет эффективных характеристик композитов с периодической структурой методом конечных элементов // Вестник МГТУ им. Н.Э.Баумана. Естественные науки. 2002. № 2. С. 95-108.
7. Димитриенко Ю.И., Соколов А.П. Разработка системы автоматизированного вычисления эффективных упругих характеристик композитов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2008. № 2. С. 57-67.
Димитриенко Ю.И., Морозов А.Н., Соколов А.П., Ничеговский Е.С. Моделирование эффективных пьезоэлектроупругих свойств композиционных материалов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2010. № 3. С. 86-97.
9. Димитриенко Ю.И., Ничеговский Е.С. Численное моделирование магнитных свойств композиционных материалов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2010. № 1. С. 3-11.
10. Димитриенко Ю.И., Соколов А.П. Метод конечных элементов для решения локальных задач механики композиционных материалов. М.: Изд-во МГТУ им. Баумана, 2010. 67 с.
11. Шелухин В.В., Терентьев С.А. Гомогенизация уравнений Максвелла и дисперсия Максвелла-Вагнера // Доклады Академии Наук. 2009. Т. 424, № 3. С. 402-406.
12. Димитриенко Ю.И. Механика сплошной среды : учеб. пособие. В 4 т. Т. 2. Универсальные законы механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э.Баумана, 2011. 560 с.
13. Димитриенко Ю.И. Тензорное исчисление. М.: Высшая школа, 2001. 576 с.
SCIENTIFIC PERIODICAL OF THE RAIJMAN MS TU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-040S
electronic scientific and technical journal
Modeling of dielectric properties of composite materials on the basis
of asymptotic averaging
# 01, January 2013
DOI: 10.7463/0113.0531682
Dimitrienko Yu.I., Sokolov A.P., Markevich M.N.
Russia, Bauman Moscow State Technical University
[email protected] al sokol o@bmstu. ru [email protected]
In this paper, the authors consider calculation of effective dielectric parameters of composite materials with complex three-dimensional armoring. The method of asymptotic averaging of periodic structures is used for averaging of long-term quasi-static electric fields in composites with dielectric components. Series of local electrostatic tasks were formulated with the use of periodicity cells; electrostatic variational problems were also formulated. To solve spatial local problems for composite materials with three-dimensional complex armoring, the finite element numerical method was used. In this article, the authors proposed an algorithm for computation of effective tensor of dielectric permeability of three-dimensional armored composite materials. Test computations were carried out in order to obtain distribution of the local electric field, effective dielectric characteristics of three-dimensional orthogonal armored composite materials with different inclusion volume fractions.
Publications with keywords: composites, multiscale homogenization method, dielectric properties
Publications with words: composites, multiscale homogenization method, dielectric properties
References
1. Bychkov I.V., Dubrovskikh D.V., Zotov I.S., Fedii A.A. Issledovanie effektivnoi dielektricheskoi pronitsaemosti kompozitnogo materiala CaSO4-2H2O - grafit [Study of the effective permittivity of a composite material CaSO4-2H2O - graphite]. Vestnik Cheliabinskogo gosudarstvennogo universitet [Bulletin of the Chelyabinsk state University], 2011, no. 7 (222), Fizika [Physics], iss. 9, pp. 7-15.
2. Nye. J.F. Physical Properties of Crystals. London, Oxford University Press, 1957. (Russ. ed.: Nai Dzh. Fizicheskie svoistva kristallov. Moscow, Mir, 1967. 385 p.).
3. Sanches-Palencia E. Non-Homogeneous media and Vibration Theory. Berlin, Springer, 1980. (Lectures Notes in Physics, vol. 127). (Russ. ed.: Sanches-Palensiia E. Neodnorodnye sredy i teoriia kolebanii. Moscow, Mir, 1984. 472 p.).
4. Bakhvalov N.S., Panasenko G.P. Osrednenieprotsessov vperiodicheskikh sredakh [Averaging processes in periodic media]. Moscow, Nauka, 1984. 352 p.
5. Pobedria B.E. Mekhanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, MSU Publ., 1984. 336 p.
6. 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.
7. 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.
Dimitrienko Iu.I., Morozov A.N., Sokolov A.P., Nichegovskii E.S. Modelirovanie effektivnykh p'ezoelektrouprugikh svoistv kompozitsionnykh materialov [Simulation of efficient piezoelectric-elastic composite materials]. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2010, no. 3, pp. 86-97.
9. Dimitrienko Iu.I., Nichegovskii E.S. Chislennoe modelirovanie magnitnykh svoistv kompozitsionnykh materialov [Numerical simulation of magnetic properties of composites]. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2010, no. 1, pp. 3-11.
10. Dimitrienko Iu.I., Sokolov A.P. Metod konechnykh elementov dlia resheniia lokal'nykh zadach mekhaniki kompozitsionnykh materialov [The method of finite elements for the solution of local problems in the mechanics of composite materials]. Moscow, Bauman MSTU Publ., 2010. 67 p.
11. Shelukhin V.V., Terent'ev S.A. Gomogenizatsiia uravnenii Maksvella i dispersiia Maksvella-Vagnera [Homogenization of the equations of Maxwell and dispersion of the Maxwell-Wagner]. Doklady Akademii Nauk, 2009, vol. 424, no. 3, pp. 402-406.
12. Dimitrienko Iu.I. Mekhanika sploshnoi sredy : ucheb. posobie. V 4 t. T. 2. Universal'nye zakony mekhaniki i elektrodinamiki sploshnoi sredy [The mechanics of a continuous medium. In 4 vols. Vol. 2. The universal laws of mechanics and electrodynamics of continua]. Moscow, Bauman MSTU Publ., 2011. 560 p.
13. Dimitrienko Iu.I. Tenzornoe ischislenie [Tensor calculus]. Moscow, Vysshaia shkola, 2001. 576 p.
