Minsk, Belarus
employee
Metal–Polymer Research Institute of National Academy of Sciences of Belarus (the laboratory 2.1 “Mechanics of Composites and Biopolymers”, Leading researcher)
employee
Gomel, Belarus
Minsk, Belarus
GRNTI 20.01 Общие вопросы информатики
GRNTI 20.15 Организация информационной деятельности
GRNTI 20.17 Документальные источники информации
GRNTI 20.19 Аналитико-синтетическая переработка документальных источников информации
GRNTI 55.01 Общие вопросы машиностроения
GRNTI 55.03 Машиноведение и детали машин
GRNTI 55.13 Технология машиностроения
GRNTI 55.15 Литейное производство
GRNTI 55.16 Кузнечно-штамповочное производство
GRNTI 55.18 Сборочное производство
GRNTI 55.19 Резание материалов
GRNTI 55.20 Электрофизикохимическая обработка
GRNTI 55.21 Термическая и упрочняющая обработка
GRNTI 55.22 Отделка поверхностей и нанесение покрытий
GRNTI 55.23 Производство изделий из порошковых материалов
GRNTI 55.24 Производство неметаллических изделий
GRNTI 55.29 Станкостроение
GRNTI 55.30 Робототехника
GRNTI 55.31 Инструментальное производство
GRNTI 55.33 Горное машиностроение
GRNTI 55.35 Металлургическое машиностроение
GRNTI 55.37 Турбостроение
GRNTI 55.41 Локомотивостроение и вагоностроение
GRNTI 55.42 Двигателестроение
GRNTI 55.43 Автомобилестроение
GRNTI 55.45 Судостроение
GRNTI 55.47 Авиастроение
GRNTI 55.51 Подъемно-транспортное машиностроение
GRNTI 55.53 Строительное и дорожное машиностроение
GRNTI 55.55 Коммунальное машиностроение
GRNTI 55.57 Тракторное и сельскохозяйственное машиностроение
GRNTI 55.69 Прочие отрасли машиностроения
GRNTI 73.01 Общие вопросы транспорта
GRNTI 73.29 Железнодорожный транспорт
GRNTI 73.31 Автомобильный транспорт
GRNTI 73.39 Трубопроводный транспорт
GRNTI 73.41 Промышленный транспорт
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 polyvinylidene 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 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 con-nected 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.
pressure sensor, aluminum anodic oxide, polymeric piezoelectric, fiber-reinforced composite, electroelasticity
Введение
При создании сенсорных элементов микроэлектронных и микроэлектромеханических систем активно используются высокоупорядоченные тонкие керамические покрытия с наноразмерными порами, получаемые анодным оксидированием металлов [1; 2]. Перспективным является использование этих покрытий в датчиках давления матричного типа, а также в качестве компонентов умных (интеллектуальных) материалов. С этой целью поры оксида, как керамической основы, импрегнируются полимером, реализующим прямой пьезоэлектрический эффект.
Ранее авторами был предложен способ импрегнирования нанопористого анодного оксида алюминия (АОА) антифрикционными фторполимерами [3] и разработана методика прогнозирования деформационно-прочностных и триботехнических параметров получаемого полимер-керамического покрытия, как трансверсально-изотропного композита [4]. Вместе с тем функциональные показатели подобных покрытий, характеризующие эффективность их использования в качестве датчиков давления, изучены недостаточно полно.
В связи с вышесказанным, целью настоящего исследования является совершенствование ранее разработанной механико-математической модели импрегнированного полимером нанопористого АОА и описание прямого пьезоэффекта в полимерном наполнителе.
Упрощенная формулировка задачи электроупругости
Физические уравнения для электроупругого материала без учета намагничиваемости и влияния температуры можно записать в виде [5]
(1)
Здесь sij, eij – компоненты тензоров напряжений и деформаций соответственно; Di, Ei – компоненты векторов индукции и напряженности электрического поля соответственно; Cijnm, eijr, lij – компоненты тензоров модулей упругости, пьезоэлектрических постоянных и диэлектрической проницаемости соответственно; индексы i, j, n, m, r = 1, 2, 3.
При пренебрежении действием магнитного поля и отсутствии некомпенсированных объемных электрических зарядов уравнения Максвелла [6] сводятся к соотношению
. (2)
Запятая перед индексом означает производную по соответствующей координате.
Компоненты тензора sij должны удовлетворять условию равновесия элементарного объема материала:
. (3)
В системе уравнений (3) не учитываются объемные силы. Компоненты тензора деформаций связаны с компонентами вектора упругих смещений u [7]:
. (4)
Компоненты вектора напряженности электрического поля связаны с электрическим потенциалом y [8]:
. (5)
Используя соотношения (4) и (5) в физических уравнениях (1) и подставив полученные выражения в (2) и (3), получим систему четырех дифференциальных уравнений в частных производных для четырех функций координат – компонент вектора смещений и электрического потенциала.
Сформулированная связанная задача электроупругости существенно упрощается, если имеется возможность принять решение дифференциального уравнения (2) в виде
. (6)
Исходя из первых трех уравнений системы (1), выражения для компонент вектора напряженности электрического поля подставим в оставшиеся шесть физических уравнений. Таким образом, допущение (6) позволяет свести задачу электроупругости к задаче теории упругости при измененных значениях модулей упругости с последующим вычислением компонент напряженности электрического поля.
В качестве материала наполнителя нанопористого АОА был взят поливинилиденфторид (ПВДФ), характеристики которого приведены, в частности, в работах [8; 9]. Значения ненулевых компонент тензоров Cijnm, eijr, lij для рассматриваемого трансверсально-изотропного композита следующие:
C1111 = C2222 = 8,6 ГПа, C3333 = 7,1 ГПа, C1122 = 5,6 ГПа, C1133 = C2233 = 5,4 ГПа; C1313 = C2323 = 1 ГПа, C1212 = (C1111 – С1122)/2 = 1,5 ГПа;
e311 = e322 = -1,1 Кл/м2, e333 = 2,9 Кл/м2, e113 = e223 = 2,3 Кл/м2;
l11 = l22 = 14,7l0, l33 = 12,4l0.
Здесь l0 = 8,85×10-12 Ф/м – диэлектрическая проницаемость вакуума.
При использовании для данного материала допущения (6) связь компонент вектора напряженности электрического поля с компонентами тензора деформаций имеет вид
. (7)
Компоненты тензора Cijnm0 в задаче теории упругости, к которой при выполнении (6) сводится исходная задача электроупругости, задаются соотношениями
(8)
Напряженно-деформированное состояние пьезоактивного наполнителя
Рассматриваемый пьезоэлектрический пленочный датчик давления может быть представлен в виде двухслойного покрытия, размещенного на условно недеформируемом основании. Применительно к наноструктурированным покрытиям на основе анодных оксидов металлов более тонкий поверхностный композитный слой толщиной h1 образован нанопористым АОА, импрегнированным полимером (рис. 1). Метод газофазной полимеризации позволяет достичь различной степени заполнения пор в зависимости от длительности данного процесса.
а) б)
в) г)
Рис. 1. Структура композитного покрытия на основе анодного оксида алюминия
после 1 ч (а, б) и 2 ч (в, г) газофазной полимеризации октафторциклобутана [3; 4]
После соответствующей обработки ряд фторполимеров, например ПВДФ, приобретают пьезоэлектрические свойства, что используется в датчиках описанной конструкции.
Материал поверхностного слоя будем рассматривать как одноосно армированный сплошными цилиндрическими волокнами композит. Волокна ориентированы по нормали к поверхности пластины. Материалом волокон является ПВДФ, материалом матрицы – анодный оксид алюминия. Нижний слой толщиной h2 представляет собой подложку из алюминия. На верхнюю поверхность покрытия действует равномерно распределенное давление p. Нижний слой покрытия адгезионно связан с недеформируемым основанием. Общая толщина покрытия много меньше характерного размера горизонтального поперечного сечения. Следовательно, с достаточной точностью можно принять [10], что при заданном режиме нагружения и наложенных граничных условиях покрытие находится в состоянии одноосной деформации, при которой ezz ¹ 0, exx = eyy = exy = exz = eyz = 0.
При анализе деформаций верхнего слоя покрытия будем использовать трехфазную модель волоконно-армированного композита [11; 12]. Основанная на данной модели методика расчетного определения механических характеристик импрегнированного полимером покрытия АОА подробно изложена в [4]. Трехфазная модель в рассматриваемом случае представляет собой составной цилиндр, помещенный в деформируемое пространство (рис. 2), характеристики которого соответствуют эффективным характеристикам композита. Внутренний цилиндр радиусом r1 образован материалом волокна (ПВДФ), внешний слой с внутренним радиусом r1 и внешним радиусом r2 образован материалом матрицы (АОА). Значения радиусов r1 и r2 связаны с объемной долей полимерного наполнителя Vf соотношением 
.
Рис. 2. Трехфазная модель волоконно-
армированного композита
Описание напряженно-деформированного состояния модели осуществляется в цилиндрической системе координат z, r, j. Ось z совпадает с осью составного цилиндра, направлена по нормали к поверхности покрытия и совпадает с осью симметрии электромеханических характеристик ПВДФ. При действии на поверхность пластины давления p осевая деформация ezz для различных элементов трехфазной модели (волокна, матрицы и композита) будет одинаковой и совпадающей с макроскопической деформацией ezz1 верхнего (композитного) слоя:
.
Верхние индексы f, m, с означают, что данная величина относится к волокну, матрице или композиту соответственно.
В силу осевой симметрии модели и рассматриваемого режима нагружения компонента uj вектора упругих смещений точек модели будет равна нулю. С учетом вышесказанного общее решение уравнений теории упругости в перемещениях для элементов трехфазной модели будет иметь вид [11]
. (9)
Здесь Ai, Bi (i = f, m, c) – не зависящие от координат величины, определяемые из граничных условий.
Требование ограниченности смещений при r = 0 будет выполняться при Bf = 0. Макроскопическая радиальная деформация композитного слоя определяется по формуле
.
Как было сказано выше, покрытие находится в состоянии одноосной деформации. Следовательно, err1 = 0 и Ac = 0. Энергетический критерий гомогенизации Эшелби [13] будет выполняться при Bk = 0. Оставшиеся три константы определяются из условий неразрывности радиального смещения ur и радиальной компоненты тензора напряжений на границе раздела волокно – матрица (r = r1), а также условия неразрывности радиального смещения на границе раздела матрица – композит (r = r2). После выполнения математических преобразований получим:
(10)
Здесь C1111m, C1122m, C1133m – компоненты тензора модулей упругости матричного материала. Не снижая общности дальнейших выкладок, будем рассматривать АОА как изотропный линейно упругий материал, механические свойства которого характеризуются модулем Юнга EАОА = 140 ГПа и коэффициентом Пуассона nАОА = 0,32 [4]. Тогда
Нижний индекс e в (10) означает, что константы af, am, bm определены при условии адгезионной связи датчика с основанием. Для сравнения рассмотрим датчик, свободно лежащий на основании. Тогда состояние одноосной деформации пластины сменится одноосным напряженным состоянием. Вместо условия равенства нулю макроскопической радиальной деформации err1 = 0 (Ac = 0) вводится требование равенства нулю макроскопического радиального напряжения srr1 = 0. Как при наличии адгезионной связи покрытия с основанием, так и при ее отсутствии энергетический критерий гомогенизации подразумевает Bc = 0. При выполнении данного условия радиальная компонента напряжения srrc в «композитном» элементе трехфазной модели не будет зависеть от координаты r. Следовательно, при одноосном напряженном состоянии покрытия для определения констант Af, Am, Bm можно использовать условия неразрывности радиального смещения ur и радиальной компоненты тензора напряжений на границе раздела волокно – матрица (r = r1), а также условие равенства нулю радиального напряжения в матрице при r = r2 (на границе с композитом). После выполнения математических преобразований получим:
(11)
Промежуточные коэффициенты K1 и K2 определены в (10). Нижний индекс s в (11) означает, что константы af, am, bm определены при отсутствии сцепления покрытия с основанием.
Макроскопическая осевая компонента тензора напряжений определяется усреднением компоненты szz по площади поперечного сечения составного цилиндра в трехфазной модели:
. (12)
Константы af, am, bm в (12) определяются либо по формулам (10), либо из соотношений (11) – в зависимости от условий закрепления.
Макроскопическое осевое напряжение связано с приложенным к пластине внешним давлением: szz1 = -p. Используя равенство (12), определим макроскопическую осевую деформацию композитного слоя, равную осевой деформации пьезоактивного полимерного наполнителя:
. (13)
Радиальная errf и окружная ejjf деформации наполнителя совпадают с константой Af:
. (14)
Зная компоненты тензора деформаций наполнителя, по формулам (7) определим компоненты вектора напряженности электрического поля. В рассматриваемой упрощенной постановке задачи электроупругости отличной от нуля будет только осевая компонента:
. (15)
Здесь деформация ezz1 определяется соотношением (13).
Компонента напряженности электрического поля не зависит от координаты. Поэтому разность потенциалов определяется произведением осевой компоненты напряженности на толщину композитного пьезоактивного слоя. Электрический потенциал распределен только по поверхности наполнителя. Для оценки эффективного значения потенциала следует умножить его на отношение площади поверхности наполнителя к общей площади поверхности тензочувствительного покрытия, равное объемной доле наполнителя. Следовательно, разность потенциалов между поверхностью покрытия и подложкой можно определить следующим образом: 
.
Эффективность датчика будем характеризовать удельной чувствительностью, определяемой как отношение разности потенциалов U к приложенному давлению p и h1:
. (16)
На рис. 3 представлены расчетные зависимости данной характеристики от объемной доли полимерного наполнителя. Можно отметить, что при отсутствии сцепления датчика с основанием чувствительность Q практически линейно возрастает с ростом Vf вплоть до Vf » 0,8. При увеличении объемной доли наполнителя до единицы скорость роста чувствительности несколько возрастает. Для датчика, адгезионно связанного с основанием, зависимость Q(Vf) существенно нелинейна и близка к квадратичной. При Vf < 0,58 величина Q для датчика, адгезионно связанного с основанием, меньше, чем при отсутствии сцепления. При Vf > 0,58 датчик, жестко связанный с основанием, обладает большей чувствительностью, чем датчик, свободно лежащий на основании.
Рис. 3. Зависимость удельной чувствительности датчика от объемной доли
полимерного пьезоактивного наполнителя (сплошная кривая – при отсутствии
сцепления датчика с основанием; пунктирная – при адгезионной связи с основанием)
Заключение
Разработана методика расчетного определения чувствительности пьезоэлектрического датчика давления, представляющего собой тонкое покрытие из нанопористого анодного оксида алюминия, поры которого наполнены поливинилиденфторидом. Методика основана на трехфазной модели волоконно-армированного композита и упрощенной постановке задачи электроупругости. Установлено, что для рассматриваемого покрытия, свободно лежащего на жестком основании, чувствительность прямо пропорциональна объемной доле полимерного наполнителя. Для адгезионно связанного с основанием покрытия зависимость чувствительности от доли наполнителя близка к квадратичной. Расчетные оценки исследуемой характеристики датчика при различных способах крепления совпадают при 58 % объемного содержания полимерного пьезоэлектрика.
Работа выполнена при поддержке БРФФИ (проект Т18КОРГ-004).
1. Nanoporous Alumina: Fabrication, Structure, Properties and Applications / eds.: D. Losic, A. Santos. - Springer, 2015. - XIV. - 362 p.
2. Muhurov, N.I. Elektromehanicheskie mikroustroystva / N.I. Muhurov, G.I. Efremov. - Minsk: Belaruskaya navuka, 2012. - 257 s.
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. Shil'ko, S.V. Opredelenie parametrov deformirovaniya i treniya gradientnyh kompozitnyh pokrytiy, obrazovannyh chastichnym zapolneniem nanopor anodnogo oksida alyuminiya polimerom / S.V. Shil'ko, D.A. Chernous, T.V. Ryabchenko // Aktual'nye voprosy mashinovedeniya: sb. nauch. tr. OIM. - Minsk, 2017. - Vyp. 6. - S. 141-144.
5. Parton, V.Z. Elektromagnitouprugost' p'ezoelektricheskih i elektroprovodnyh tel / V.Z. Parton, B.A. Kudryavcev. - M.: Nauka, 1988. - 472 s.
6. Savel'ev, I.V. Osnovy teoreticheskoy fiziki. V 2 t. T. 1. Mehanika i elektrodinamika / I.V. Savel'ev. - 2-e izd., ispr. - M.: Nauka, Gl. red. fiz.-mat. lit., 1991. - 496 s.
7. Starovoytov, E.I. Osnovy teorii uprugosti, plastichnosti i vyazkouprugosti / E.I. Starovoytov. - Gomel': BelGUT, 2001. - 344 s.
8. Horoshun, L.P. Prognozirovanie effektivnyh svoystv p'ezoaktivnyh kompozitnyh materialov / L.P. Horoshun, B.P. Maslov, P.V. Leschenko. - Kiev: Nauk. dumka, 1989. - 208 s.
9. Sessler, G.M. Piezoelectricity in polyvinylideneflyoride / G.M. Sessler // J. Acoust. Soc. Amer. - 1981. - V. 70. - № 6. - P. 1596-1608.
10. Mozharovskiy, V.V. Prikladnaya mehanika sloistyh tel iz kompozitov / V.V. Mozharovskiy, V.E. Starzhinskiy. - Minsk: Nauka i tehnika, 1988. - 271 s.
11. Kristensen, R. Vvedenie v mehaniku kompozitov / R. Kristensen. - M.: Mir, 1982. - 334 s.
12. Dispersno-napolnennye polimernye kompozity tehnicheskogo i medicinskogo naznacheniya / B.A. Lyukshin [i dr.]. - Novosibirsk: Izd-vo SO RAN, 2017. - 311 s.
13. Eshelbi, Dzh. Kontinual'naya teoriya dislokaciy / Dzh. Eshelbi. - M.: Inostr. lit., 1963. - 247 s.
