Расчет чувствительности пьезоэлектрического датчика давления на основе импрегнированного полимером анодного оксида алюминия Текст научной статьи по специальности «Физика»

УДК 681.586.773

DOI: 10.30987/article_5d2d923295da08.42628593

С.В. Шилько, Д.А. Черноус, Т.В. Рябченко

РАСЧЕТ ЧУВСТВИТЕЛЬНОСТИ ПЬЕЗОЭЛЕКТРИЧЕСКОГО ДАТЧИКА ДАВЛЕНИЯ НА ОСНОВЕ ИМПРЕГНИРОВАННОГО ПОЛИМЕРОМ АНОДНОГО ОКСИДА АЛЮМИНИЯ

Разработана методика расчетного определения чувствительности датчика давления, представляющего собой тонкое покрытие из нанопористого анодного оксида алюминия, поры которого наполнены полимерным пьезоэлектриком. Датчик моделируется тонким двухслойным покрытием, размещенным на условно недеформируемом основании. Методика основана на трехфазной модели воло-конно-армированного композита и упрощенной

постановке задачи электроупругости. При моделировании прямого пьезоэффекта на примере поли-винилиденфторида в качестве наполнителя показано, что зависимость чувствительности от доли наполнителя близка к квадратичной.

Ключевые слова: датчик давления, анодный оксид алюминия, полимерный пьезоэлектрик, волоконно-армированный композит, электроупругость.

S.V. Shilko, D A. Chernous, T.V. Ryabchenko

SENSITIVITY CALCULATION OF PIEZOELECTRIC PRESSURE SENSOR BASED ON ALUMINUM ANODIC OXIDE IMPREGNATED WITH POLYMER

The application of highly ordered thin ceramic coatings with nano-dimensional pores obtained by metal anodic oxidation in pressure sensors of a matrix type is promising. For this purpose the pores of a ceramic base are impregnated with polymer realizing a direct piezoelectric effect.

There is developed a procedure for the calculated definition of piezoelectric pressure sensor sensitivity which is a thin coating made of nano-porous aluminum anodic oxide the pores of which are filled with polyvi-nylidene fluoride.

The procedure is based on a three-phase model of fiber-reinforced composite and a simplified problem setting of electroelasticity. The sensor is modeled by a

Введение

При создании сенсорных элементов микроэлектронных и микроэлектромеханических систем активно используются высокоупорядоченные тонкие керамические покрытия с наноразмерными порами, получаемые анодным оксидированием металлов [1; 2]. Перспективным является использование этих покрытий в датчиках давления матричного типа, а также в качестве компонентов умных (интеллектуальных) материалов. С этой целью поры оксида, как керамической основы, импрегни-руются полимером, реализующим прямой пьезоэлектрический эффект.

Ранее авторами был предложен способ импрегнирования нанопористого

thin two-layer coating located on a conditionally non-deformable base.

It is defined that for a coating lying freely on a rigid foundation the sensitivity is directly proportional to volume fraction of the polymer filler. For a coating connected in an adhesion way with the base the sensitivity dependence upon a filler share is close to a quadratic one. The estimated assessments of sensor characteristics under analysis at different ways of fastening coincide at 58% volume content of polymer piezoelectric.

Key words: pressure sensor, aluminum anodic oxide, polymeric piezoelectric, fiber-reinforced composite, electroelasticity.

анодного оксида алюминия (АОА) антифрикционными фторполимерами [3] и разработана методика прогнозирования деформационно-прочностных и триботехни-ческих параметров получаемого полимер-керамического покрытия, как трансвер-сально-изотропного композита [4]. Вместе с тем функциональные показатели подобных покрытий, характеризующие эффективность их использования в качестве датчиков давления, изучены недостаточно полно.

В связи с вышесказанным, целью настоящего исследования является совершенствование ранее разработанной механико-математической модели импрегниро-

ванного полимером нанопористого АОА и описание прямого пьезоэффекта в поли-

мерном наполнителе.

Упрощенная формулировка задачи электроупругости

Физические уравнения для электроупругого материала без учета намагничи-ваемости и влияния температуры можно записать в виде [5]

П — ' Е + Р Р

Пг ' у Е] + Р Пт Р пт

СТ,

— С р — р Е

гупт пт РУ Р

(1)

Здесь Сту, Ру - компоненты тензоров напряжений и деформаций соответственно; Пу, Еу - компоненты векторов индукции и напряженности электрического поля соответственно; Супт, Рур, 'у - компоненты тензоров модулей упругости, пьезоэлектрических постоянных и диэлектрической проницаемости соответственно; индексы , у, п, т, р = 1, 2, 3.

При пренебрежении действием магнитного поля и отсутствии некомпенсированных объемных электрических зарядов уравнения Максвелла [6] сводятся к соотношению

ёпЪ — Ц г — 0 . (2)

Запятая перед индексом означает производную по соответствующей координате.

Компоненты тензора Сту должны удовлетворять условию равновесия элементарного объема материала:

ст., - 0. (3)

В системе уравнений (3) не учитываются объемные силы. Компоненты тензора деформаций связаны с компонентами вектора упругих смещений и [7]:

1

8'- 2

дн ды

Л

-+-

дх, дх

V ] 1

(4)

Компоненты вектора напряженности электрического поля связаны с электрическим потенциалом у [8]:

Е — — ^га^у ^ Е — —У г. (5)

Используя соотношения (4) и (5) в физических уравнениях (1) и подставив полученные выражения в (2) и (3), получим систему четырех дифференциальных уравнений в частных производных для четырех функций координат - компонент вектора смещений и электрического потенциала.

Сформулированная связанная задача электроупругости существенно упрощается, если имеется возможность принять решение дифференциального уравнения (2) в виде

Б, — 0. (6)

Исходя из первых трех уравнений системы (1), выражения для компонент вектора напряженности электрического поля подставим в оставшиеся шесть физических уравнений. Таким образом, допущение (6) позволяет свести задачу электроупругости к задаче теории упругости при измененных значениях модулей упругости с последующим вычислением компонент напряженности электрического поля.

В качестве материала наполнителя нанопористого АОА был взят поливини-лиденфторид (ПВДФ), характеристики которого приведены, в частности, в работах [8; 9]. Значения ненулевых компонент тензоров Сгупт, Рур, Ау для рассматриваемого трансверсально-изотропного композита следующие:

С1111 = С2222 = 8,6 ГПа, С3333 = 7,1 ГПа, С1122 = 5,6 ГПа, С1133 = С2233 = 5,4 ГПа; С1313 = С2323 = 1 ГПа, С1212 = (С1111 - Сш2)/2 = 1,5 ГПа; Р311 = Р322 = -1,1 Кл/м2, Р333 = 2,9 Кл/м2, Р113 = Р223 = 2,3 Кл/м2; '11 = '22 = 14,7'о, '33 = 12,4'о.

Здесь 'о = 8,85-10-12 Ф/м - диэлектрическая проницаемость вакуума.

При использовании для данного материала допущения (6) связь компонент

вектора напряженности электрического поля с компонентами тензора деформаций имеет вид

Е — —

М13

'11

'13

Е — — Е 2 —

'223

'23

Е — —

Е 3 —

22

^ [

Р311 ( Р11 + Р 22

) + Р

333 Р 33

(7)

33

Компоненты тензора Сцпш0 в задаче теории упругости, к которой при выполне-

нии (6) сводится исходная задача электроупругости, задаются соотношениями

е2 е2

с = с = с + — с = с +е311 с11110 с 22220 с1111 + - , с1 1220 с1122 +

Л33

Л33

с =с =с +е31е333 с =с +

с1 1330 с22330 с1133 + п , с33330 с3333 +

333

Л33

Л

(8)

33

с = с = с +-211 с = с

с13130 с23230 с1313 + ~ , с12120 с1212

Л11

Напряженно-деформированное состояние пьезоактивного наполнителя

Рассматриваемый пьезоэлектрический пленочный датчик давления может быть представлен в виде двухслойного покрытия, размещенного на условно неде-формируемом основании. Применительно к наноструктурированным покрытиям на основе анодных оксидов металлов более

тонкий поверхностный композитный слой толщиной к\ образован нанопористым АОА, импрегнированным полимером (рис. 1). Метод газофазной полимеризации позволяет достичь различной степени заполнения пор в зависимости от длительности данного процесса.

в) г)

Рис. 1. Структура композитного покрытия на основе анодного оксида алюминия после 1 ч (а, б) и 2 ч (в, г) газофазной полимеризации октафторциклобутана [3; 4]

После соответствующей обработки ряд фторполимеров, например ПВДФ, приобретают пьезоэлектрические свойства, что используется в датчиках описанной конструкции.

Материал поверхностного слоя будем рассматривать как одноосно армированный сплошными цилиндрическими волокнами композит. Волокна ориентирова-

ны по нормали к поверхности пластины. Материалом волокон является ПВДФ, материалом матрицы - анодный оксид алюминия. Нижний слой толщиной к2 представляет собой подложку из алюминия. На верхнюю поверхность покрытия действует равномерно распределенное давление р. Нижний слой покрытия адгезионно связан с недеформируемым основанием. Общая

2

2

толщина покрытия много меньше характерного размера горизонтального поперечного сечения. Следовательно, с достаточной точностью можно принять [10], что при заданном режиме нагружения и наложенных граничных условиях покрытие находится в состоянии одноосной деформации, при которой Ргг ^ 0, Рхх Руу Ржу Рхг Рyz 0.

При анализе деформаций верхнего слоя покрытия будем использовать трехфазную модель волоконно-армированного композита [11; 12]. Основанная на данной модели методика расчетного определения механических характеристик импрегниро-ванного полимером покрытия АОА подробно изложена в [4]. Трехфазная модель в рассматриваемом случае представляет собой составной цилиндр, помещенный в деформируемое пространство (рис. 2), характеристики которого соответствуют эффективным характеристикам композита. Внутренний цилиндр радиусом п образован материалом волокна (ПВДФ), внешний слой с внутренним радиусом п и внешним радиусом Г2 образован материалом матрицы (АОА). Значения радиусов Г1 и Г2 связаны с объемной долей полимерного

( „ V

наполнителя ¥{ соотношением —

V Г2 У

Описание напряженно-деформированного состояния модели осуществляется в цилиндрической системе координат г, г, ф. Ось г совпадает с осью составного цилиндра, направлена по нормали к поверхности покрытия и совпадает с осью симметрии электромеханических характеристик

ы{ — ыт — ыС — е„л 2; 2 2 2 221 ?

Здесь Л1, Б1 (У = / т, с) - не зависящие от координат величины, определяемые из граничных условий.

Требование ограниченности смещений при г = 0 будет выполняться при Б = 0. Макроскопическая радиальная деформация композитного слоя определяется по формуле

ёы

ёт

— АС

Как было сказано выше, покрытие находится в состоянии одноосной дефор-

ПВДФ. При действии на поверхность пластины давления р осевая деформация Ргг для различных элементов трехфазной модели (волокна, матрицы и композита) будет одинаковой и совпадающей с макроскопической деформацией Ргг1 верхнего (композитного) слоя:

Р 22 Р 22 Р 22 Р 221 .

Рис. 2. Трехфазная модель волоконно-армированного композита

Верхние индексы /, т, с означают, что данная величина относится к волокну, матрице или композиту соответственно.

В силу осевой симметрии модели и рассматриваемого режима нагружения компонента % вектора упругих смещений точек модели будет равна нулю. С учетом вышесказанного общее решение уравнений теории упругости в перемещениях для элементов трехфазной модели будет иметь вид [11]

г

ы;„ —

ыф — ыф — 0; ы'т — Ат +

В'

(9)

мации. Следовательно, Ргг1 = 0 и Лс = 0. Энергетический критерий гомогенизации Эшелби [13] будет выполняться при Бк = 0. Оставшиеся три константы определяются из условий неразрывности радиального смещения иг и радиальной компоненты тензора напряжений на границе раздела волокно - матрица (г = п), а также условия неразрывности радиального смещения на границе раздела матрица-композит (г = г2). После выполнения математических преобразований получим:

с

т

Ртт1

Af =Е , a

f

zz1^ S 5

Am = s „tf

Bm = r2S bm

B r1 S zz1bS

am

Sm

m m m

133

m m m m

^—1111 + *1122 C(K 2 1)(*1111 *1122)

af=K + K2 aSm,

bm = af - am

K1

m

*1133 —

11330

nf -i-Cf 4-Cm —Cm *11110 + *11220 + *—1111 *1122

(10)

m

TS — TS 2^1111 K = K-т

2 1 si m _ /~if

*1133 *11

1330

Здесь Сш1т, Сш2т, Спззт - компоненты тензора модулей упругости матричного материала. Не снижая общности дальнейших выкладок, будем рассматривать АОА как изотропный линейно упру-

/im _ /im _ /im

*—1111 = * 2222 = *333

гий материал, механические свойства которого характеризуются модулем Юнга Еаоа = 140 ГПа и коэффициентом Пуассона УАОА = 0,32 [4]. Тогда

ЕАОА — ^АОА )

mmm * 1122 = *1133 = * 2233 =

Нижний индекс в в (10) означает, что константы О, ат, Ьт определены при условии адгезионной связи датчика с основанием. Для сравнения рассмотрим датчик, свободно лежащий на основании. Тогда состояние одноосной деформации пластины сменится одноосным напряженным состоянием. Вместо условия равенства нулю макроскопической радиальной деформации вГг1 = 0 (Ас = 0) вводится требование равенства нулю макроскопического радиального напряжения аГг1 = 0. Как при наличии адгезионной связи покрытия с основанием, так и при ее отсутствии энергетический критерий гомогенизации подра-

Лг =е,л а,, Ли =в

am,

zz1 а '

a„ = —

VfK1

1 + V,(К2 -1)'

Промежуточные коэффициенты К1 и К2 определены в (10). Нижний индекс а в (11) означает, что константы О, ат, Ьт определены при отсутствии сцепления покрытия с основанием.

af =

(1 +^аоа)(1 - 2VAOA )' E v

_еАОА v АОА_

(1 +VAOA)(1 - 2VAOA )'

зумевает Bc = 0. При выполнении данного условия радиальная компонента напряжения Orrc в «композитном» элементе трехфазной модели не будет зависеть от координаты r. Следовательно, при одноосном напряженном состоянии покрытия для определения констант Af, Am, Bm можно использовать условия неразрывности радиального смещения Ur и радиальной компоненты тензора напряжений на границе раздела волокно - матрица (r = ri), а также условие равенства нулю радиального напряжения в матрице при r = r2 (на границе с композитом). После выполнения математических преобразований получим:

Bm = rSbm,

(ii)

K + KaZ-

1 2 а -

bm = af — am.

а а а

а zz 1 =

2

Л

Ja fzrdr + Ja

Макроскопическая осевая компонента тензора напряжений определяется усреднением компоненты Охх по площади поперечного сечения составного цилиндра в трехфазной модели:

Г2 >

rdr

(12)

bm в (12) определя-

Константы а , ат ются либо по формулам (10), либо из соотношений (11) - в зависимости от условий закрепления.

Макроскопическое осевое напряжение связано с приложенным к пластине

Г1 У

внешним давлением: а2г1 = -р. Используя равенство (1 2), определим макроскопическую осевую деформацию композитного слоя, равную осевой деформации пьезоак-тивного полимерного наполнителя:

2

r

2

0

р 22 р 221

С ( С33330 + 2 а^С1У1330 ) + (1 С )(С3333 + 2 а С1133

Радиальная Рг/ и окружная Рфф^ деформации наполнителя совпадают с константой Л:

. (14)

Зная компоненты тензора деформаций наполнителя, по формулам (7) определим компоненты вектора напряженности электрического поля. В рассматриваемой

'т + 2 атСт33)

(13)

а'

22 1

упрощенной постановке задачи электроупругости отличной от нуля будет только осевая компонента:

Е2 [Р311(

' 33

Р / + Р ) + Р™ ТР.

' гг

фф

333 22

■ ]

'22 1

(2Р311 + Р333 ) .

(15)

33

Здесь деформация Р221 определяется соотношением (13).

Компонента напряженности электрического поля не зависит от координаты. Поэтому разность потенциалов определяется произведением осевой компоненты напряженности на толщину композитного пьезоактивного слоя. Электрический потенциал распределен только по поверхности наполнителя. Для оценки эффективного значения потенциала следует умножить его на отношение площади поверхности наполнителя к общей площади поверхно-

Q-U

сти тензочувствительного покрытия, равное объемной доле наполнителя. Следовательно, разность потенциалов между поверхностью покрытия и подложкой можно определить следующим образом:

и — 1у,Ег.

Эффективность датчика будем характеризовать удельной чувствительностью, определяемой как отношение разности потенциалов и к приложенному давлению р и М:

— V,

1 + e'

3333

¡311-1 с(С33330 + C1133o) + (1 С)(С3333 + Cllзз)

(16)

На рис. 3 представлены расчетные зависимости данной характеристики от объемной доли полимерного наполнителя. Можно отметить, что при отсутствии сцепления датчика с основанием чувствительность Q практически линейно возрастает с ростом V/ вплоть до V/« 0,8. При увеличении объемной доли наполнителя до единицы скорость роста чувствительности несколько возрастает. Для датчика, ад-

гезионно связанного с основанием, зависимость Q(Vf) существенно нелинейна и близка к квадратичной. При V/ < 0,58 величина Q для датчика, адгезионно связанного с основанием, меньше, чем при отсутствии сцепления. При V/> 0,58 датчик, жестко связанный с основанием, обладает большей чувствительностью, чем датчик, свободно лежащий на основании.

Q,

В/Па-м

Рис. 3. Зависимость удельной чувствительности датчика от объемной доли полимерного пьезоактивного наполнителя (сплошная кривая - при отсутствии сцепления датчика с основанием; пунктирная - при адгезионной связи с основанием)

0,05

0

0

0,2

0,4

0,6

0,8

с

Заключение

Разработана методика расчетного определения чувствительности пьезоэлектрического датчика давления, представляющего собой тонкое покрытие из нанопо-ристого анодного оксида алюминия, поры которого наполнены поливинилиденфто-ридом. Методика основана на трехфазной модели волоконно-армированного композита и упрощенной постановке задачи электроупругости. Установлено, что для рассматриваемого покрытия, свободно ле-

Работа выполнена при поддержке БРФФИ СПИСОК ЛИТЕРАТУРЫ

1. Nanoporous Alumina: Fabrication, Structure, Properties and Applications / еds.: D. Losic, A. Santos. -Springer, 2015. - XIV. - 362 p.

2. Мухуров, Н.И. Электромеханические микроустройства / Н.И. Мухуров, Г.И. Ефремов. -Минск: Беларуская навука, 2012. - 257 с.

3. Shil'ko, S.V. Estimation of the friction coefficient of a nanostructured composite coating / S.V. Shil'ko, D. Chernous, T.V. Ryabchenko, V.V. Hat'ko // Mechanics of Composite Materials. -2017. - Vol. 53. - № 5. - P. 579-588.

4. Шилько, С.В. Определение параметров деформирования и трения градиентных композитных покрытий, образованных частичным заполнением нанопор анодного оксида алюминия полимером / С.В. Шилько, Д.А. Черноус, Т.В. Рябченко // Актуальные вопросы машиноведения: сб. науч. тр. ОИМ. - Минск, 2017. - Вып. 6. - С. 141-144.

5. Партон, В.З. Электромагнитоупругость пьезоэлектрических и электропроводных тел / В.З. Партон, Б.А. Кудрявцев. - М.: Наука, 1988. -472 с.

6. Савельев, И.В. Основы теоретической физики. В 2 т. Т. 1. Механика и электродинамика / И.В.

жащего на жестком основании, чувствительность прямо пропорциональна объемной доле полимерного наполнителя. Для адгезионно связанного с основанием покрытия зависимость чувствительности от доли наполнителя близка к квадратичной. Расчетные оценки исследуемой характеристики датчика при различных способах крепления совпадают при 58 % объемного содержания полимерного пьезоэлектрика.

ч Т18КОРГ-004).

Савельев. - 2-е изд., испр. - М.: Наука, Гл. ред. физ.-мат. лит., 1991. - 496 с.

7. Старовойтов, Э.И. Основы теории упругости, пластичности и вязкоупругости / Э.И. Старовойтов. - Гомель: БелГУТ, 2001. - 344 с.

8. Хорошун, Л.П. Прогнозирование эффективных свойств пьезоактивных композитных материалов / Л.П. Хорошун, Б.П. Маслов, П.В. Лещен-ко. - Киев: Наук. думка, 1989. - 208 с.

9. Sessler, G.M. Piezoelectricity in polyvinylidene-flyoride / G.M. Sessler // J. Acoust. Soc. Amer. -1981. - V. 70. - № 6. - P. 1596-1608.

10. Можаровский, В.В. Прикладная механика слоистых тел из композитов / В.В. Можаровский, В.Е. Старжинский. - Минск: Наука и техника, 1988. - 271 с.

11. Кристенсен, Р. Введение в механику композитов / Р. Кристенсен. - М.: Мир, 1982. - 334 с.

12. Дисперсно-наполненные полимерные композиты технического и медицинского назначения / Б.А. Люкшин [и др.]. - Новосибирск: Изд-во СО РАН, 2017. - 311 с.

13. Эшелби, Дж. Континуальная теория дислокаций / Дж. Эшелби. - М.: Иностр. лит., 1963. - 247 с.

1. Nanoporous Alumina: Fabrication, Structure, Properties and Applications / eds.: D. Losic, A. Santos. -Springer, 2015. - XIV. - 362 p.

2. Mukhurov, N.I. Electromechanical Micro-devices / N.I. Mukhurov, G.I. Efremov. - Minsk: belorussian Science, 2012. - pp. 257.

3. Shil'ko, S.V. Estimation of the friction coefficient of a nanostructured composite coating / S.V. Shil'ko, D. Chernous, T.V. Ryabchenko, V.V. Hat'ko // Mechanics of Composite Materials. -2017. - Vol. 53. - № 5. - P. 579-588.

4. Shilko, S.V. Parameter definition of deformation and friction of gradient composite coatings formed with partial filling of aluminum anodic oxide nano-pores by polymer / S.V. Shilko, D.A. Chernous, T.V. Ryabchenko // Actual Problems of Machine

Science: OIM Proceedings. - Minsk, 2017. - edition 6. - pp. 141-144.

5. Parton, V.Z. Electromagnetoelasticity of Piezoelectric and Electroconductive Bodies / V.Z. Parton, B.A. Kudryavtsev. - M.: Science, 1988. - pp. 472.

6. Saveliev, I.V. Fundamentals of Theoretical Physics. In 2 Vol., Vol.1. Mechanics and Electrodynamics / I.V. Saveliev. - 2-d Edition revised. - M.: Science: Main editorial board of phys.-math. Lit., 1991. -pp. 496.

7. Starovoitov, E.I. Fundamentals of Theory of Elasticity, Plasticity, Viscous Elasticity / E.I, Starovoi-tov. - Gomel: BelSUT, 2001. - pp. 344.

8. Khoroshun, L.P. Prediction of Piezoactive Composite Efficient Properties / L.P. Khoroshun, B.P. maslov, P.V. Leshchenko. - Kiev: Scientific Thought, 1989. - pp. 208.

9. Sessler, G.M. Piezoelectricity in polyvinylidene-flyoride / G.M. Sessler // J. Acoust. Soc. Amer. -1981. - V. 70. - № 6. - P. 1596-1608.

10. Mozharovsky, V.V. Applied Mechanics of Laminated Bodies and Their Composites / V.V. Mozhrovsky, V.E. Starzhinsky. - Minsk: Science and Engineering, 1988. - pp. 271.

11. Cristensen, R. Introduction into Composite Mechanics / R. Cristensen. - M.: World, 1982. - pp. 334.

12. Dispersed-filled Polymeric Composites of Engineering and Medical Purpose / B.A. Lyukshin [et al.]. - Novosibirsk: Publishing house of SB RAS, 2017. - pp. 311.

13. Eshelby, J. Continuum Theory of Dislocations / J. Eshelby. - M.: Foreign Lit., 1963. - pp. 247.

Статья поступила в редакцию 18.04.19 Рецензент: д.т.н., профессор Белорусского государственного

университета транспорта Тариков Г.П.

Статья принята к публикации 7. 06. 19.

Сведения об авторах:

Шилько Сергей Викторович, к.т.н., доцент, зав. лабораторией «Механика композитов и биополимеров» ИММС НАН Беларуси, e-mail: shilko mpri@mail.ru.

Черноус Дмитрий Анатольевич, к.т.н., доцент, вед. науч. сотрудник ИММС НАН Беларуси, тел.: +375(232)34-06-31.

Рябченко Татьяна Владимировна, науч. сотрудник ИММС НАН Беларуси, e-mail: rits@tut.by.

Shilko Sergey Victorovich, Can. Sc. Tech., Assistant Prof., Head of the Lab. "Mechanics of Composites and Biopolymers", IMMS NAS of Belarus, e-mail: shilko mpri@mail.ru.

Chernous Dmitry Anatolievich, Can. Sc. Tech., Assistant Prof., Leading researcher, IMMS NAS of Belarus, phone: +375(232)34-06-31. Рябченко Татьяна Владимировна, науч. сотрудник ИММС НАН Беларуси, e-mail: rits@tut.by.