На правах рукописи

ЛЫЖОВ ВЯЧЕСЛАВ АЛЕКСАНДРОВИЧ

ДИНАМИЧЕСКИЕ СМЕШАННЫЕ ЗАДАЧИ

ДЛЯ ПЬЕЗОАКТИВНЫХ НЕОДНОРОДНЫХ

СТРУКТУР

Специальность 01.02.04 – Механика деформируемого

твердого тела

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Ростов-на-Дону 2013

Работа выполнена в комплексном отделе механики, химии, физики и нанотехнологий Федерального государственного бюджетного учреждения науки Южного научного центра Российской академии наук Научный доктор физико-математических наук, профессор руководитель: Калинчук Валерий Владимирович Официальные Пряхина Ольга Донатовна оппоненты: доктор физико-математических наук, профессор Суворова Татьяна Виссарионовна доктор физико-математических наук, доцент ФГБОУ ВПО РГУПС, профессор Ведущая ФГБОУ ВПО "Нижегородский государственный университет организация: им. Н.И. Лобачевского" (НИИ механики)

Защита состоится 24 декабря 2013 г. в 12-00 на заседании диссертационного совета 212.101.07 при ФГОУ ВПО "Кубанский государственный университет" по адресу: 350040, г. Краснодар, ул. Ставропольская 149, ауд 231.

С диссертацией можно ознакомиться в научной библиотеке Кубанского государственного университета по адресу: 350040, г. Краснодар, ул. Ставропольская 149.

Автореферат разослан «23» ноября 2013 г.

Ученый секретарь диссертационного совета Зарецкая М.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования.

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

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

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

Научная новизна результатов работы Для небольшого количества электродов (менее 20) предложен и реализован строгий метод, позволяющий исследовать тонкие вопросы процесса взаимодействия электродов с предварительно напряженной пьезоактивной средой. Исследованы особенности распределения полей при взаимодействии ПАВ с многоэлектродными структурами, изучено влияние начальных напряжений на динамические характеристики волнового процесса на поверхности среды.

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

феноменологической теории, учитывающей возникновение напряжений в тонкой пленке за счет деформации несоответствия (misfit strain) из-за разницы параметров решетки пленки и подложки. Метод реализован на структуре BST0.8/MgO.

Исследовано влияние толщины пленки на область существования и фазовые скорости распространяющихся поверхностных акустических.

акустических волн для ряда неоднородных пьезоактивных структур, находящихся в начальном напряженном состоянии, разработана методика определения электрических, механических полей и динамических характеристик для систем планарных электродов на неоднородной пьезоактивной подложке. В процессе выполнения исследований показано, что в зависимости от вида и величины начальных напряжений изменяются скорости поверхностных волн и частотные диапазоны их существования. Изучены особенности распределения физических полей в зависимости от геометрических полупространстве. Установлено, что при уменьшении толщины слоя существенно увеличивается неоднородность распределения электрического поля и возникает новый максимум нормальной компоненты электрической индукции на границе раздела пьезоэлектрик – диэлектрик, обусловленный наличием механических напряжений.

Исследовано влияние начальных напряжений на фазовые скорости волн ГуляеваБлюштейна, а также на коэффициент электромеханической связи среды.

Научное и практическое значение результатов работы.

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

Важной практической задачей является развитие методов моделирования акустоэлектронных устройств, выполненных из неоднородных пьезоактивных материалов и использующих в качестве возбуждающего или принимающего элемента многоэлектродные структуры. Используемые методы исследований и полученные результаты нашли свое применение при выполнении работ по плановой НИР №г/р 01201354242, грантам РФФИ 05-08-65499, 06-01-00726, 09-01-00695, 10-08-01082, 11-08-00884-а, 12-08-31219 мол_а.

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

Достоверность результатов Достоверность результатов диссертационной работы основана на использовании строгих математических методов электроупругости; на использовании определяющих соотношений линеаризованной теории распространения упругих волн; на детальном исследовании свойств ядер интегральных представлений для разработки устойчивых численных алгоритмов; на сравнении результатов решения рассмотренных задач с точными решениями и результатами других исследователей.

На защиту выносится:

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

многоэлектродной структуры (количество электродов более 20) с неоднородной предварительно напряженной пьезоактивной средой, основанная на использовании приближенного метода решения интегрального уравнения на системе отрезков;

исследования закономерностей взаимодействия многоэлектродных структур с неоднородными предварительно напряженными пьезоактивными средами, изучения особенностей поверхностных волновых полей;

4. математическая модель, описывающая динамику гетероструктуры с покрытием из тонкой сегнетоэлектрической пленки, основанная на использовании феноменологической теории, учитывающей возникновение напряжений за счет деформации несоответствия (misfit strain) из-за разницы параметров решетки пленки и подложки и линейной деформации подложки при изменении температуры;

5. результаты анализа влияния начальных напряжений на характеристики волнового поля, инициированного многоэлектродными структурами;

неоднородных пьезоактивных структур, находящихся в начальном напряженном состоянии;

7. методика определения электрических, механических полей и динамических характеристик для систем планарных электродов на неоднородной пьезоактивной подложке;

8. закономерности влияния вида и величины начальных напряжений в неоднородных пьезоактивных средах на скорости поверхностных волн и частотные диапазоны их существования;

геометрических параметров задачи для пьезоэлектрического слоя на диэлектрическом полупространстве;

10. закономерности влияния начальных напряжений на фазовые скорости волн Гуляева-Блюштейна, а также на коэффициент электромеханической связи среды;

11. закономерности влияния параметров пленки на область существования и фазовые скорости распространяющихся поверхностных акустических волн.

Апробация работы Результаты диссертационной работы докладывались на XI международной конференции «Современные проблемы механики сплошной среды» (Ростов-на-Дону, 2007), Всероссийской конференции "Проблемы нелинейной механики деформируемого твердого тела" (Пермь, 2008), XIV международной конференции «Современные проблемы механики сплошной среды» (Ростов-на-Дону, Азов, 2010), X Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики (Н.Новгород, 2011), международной конференции «Современные проблемы механики», посвященной 100-летию Л.А Галина (Москва, 2012), VII Всероссийской (с международным участием) конференции по механике деформируемого твердого тела (Ростов-на-Дону, 2013), а также на семинарах отдела волновых процессов НИИ механики и прикладной математики им. Воровича И.И. и на семинарах комплексного отдела механики, химии, физики и нанотехнологии ЮНЦ РАН.

Публикации Основные результаты исследований, выполненных по теме диссертации, содержатся в 13 публикациях, из них 5 работ – в рецензируемых изданиях из списка ВАК РФ.

Личный вклад автора Выбор темы, цели и задачи, а также возможные пути проведения исследований проводились совместно с научным руководителем В.В. Калинчуком. Обсуждение полученных результатов, их анализ и формулировка выводов, выявление особенностей и общих закономерностей осуществлялись автором диссертационной работы совместно с научным руководителем В.В. Калинчуком и ведущим научным сотрудником Т.И.

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

Структура и объем работы Диссертационная работа состоит из введения, четырех глав, заключения, списка литературы из 121 наименования. Общий объем работы составляет 106 страниц машинописного текста, включая 42 рисунка и 3 таблицы.

Содержание работы Во введении приводится обзор современного состояния вопроса в литературе и обоснование актуальности учета начальных напряжений и решения задач в строгой электроупругих средах. Существенный вклад в разработку различных вопросов по теме исследования внесли В.А. Бабешко, Т.И. Белянкова, М.А. Балакирев, Д.И. Бардзокас, А.В. Белоконь, А.О. Ватульян, И.А. Викторов, И.И. Ворович, Е.В. Глушков, Н.В.

Н.Ф. Морозов, А.В. Наседкин, В. Новацкий, В.З. Партон, О.Д. Пряхина, Д. Руайе, Н.А.

Сеник, А.Н. Соловьев, А.В. Смирнова, А.Ф. Улитко, Ю.А. Устинов, Дж. Фарнелл, М.Л. Фильштинский, Н.А. Шульга, J.C. Baumhauer, L. Bergman, K.Bragger, J. Helsing, R.B. Hetnarski, Q.H. Qin, S. Shen, H.F. Tiersten, R.N. Thurston и другие исследователи.

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

Актуальная конфигурация (АК) рассматривается в виде наложения малых возмущений на некоторое начальное конечное деформированное состояние (НДС), которое предполагается однородным. Используются линеаризованные уравнения движения электроупругости и квазиэлектростатики компоненты линеаризованных тензора напряжений и вектора индукции имеют вид В выражениях (2) коэффициенты при частных производных играют такую же роль, что и обычные упругие, пьезоэлектрические и диэлектрические модули, но зависят как от свойств материала подложки, так и от вида и величины начальных напряжений:

Здесь Clksp, e plk элементарного объема в направлении i-ой оси, Plp – компоненты тензора напряжений напряженного состояния определяются из решения задачи для начальноPlp деформированного состояния.

Из выражений (3) следует, что различные виды начального напряженного состояния могут приводить к изменению класса симметрии кристалла, поэтому при построении математической модели рассматривается материал класса 2mm, в рамках которого удается описать кристаллы классов 6mm, 4mm и 2mm, а также пьезокерамику в условиях воздействия начальных напряжений.

Рис.1. Геометрия задачи о колебаниях слоисто неоднородного пьезоактивного полупространства под действием системы поверхностных планарных электродов.

Модель акустоэлектронного устройства представляет собой систему из N тонких протяженных электродов с заданными потенциалами на поверхности слоисто неоднородного пьезоактивного полупространства (рис. 1). Поверхность внешнего слоя свободна от механических напряжений. На внутренних границах между слоями предполагается непрерывность полей смещений и напряжений, потенциала и нормальной компоненты электрической индукции. На бесконечности учитываются условия излучения Зоммерфельда.

В большинстве случаев можно считать продольные размеры электродов в значительной степени превосходящими поперечные. Если к тому же свойства подложки однородны вдоль оси x3, можно считать все параметры независящими от координаты x3 и рассматривать плоскую задачу для волн типа Рэлея или антиплоскую для волн типа Лява и Гуляева-Блюштейна.

Во второй главе получены элементы матрицы функции Грина и представление решений сформулированных краевых задач для основных типов пьезоактивных подложек – однородного с учетом вакуума, слоисто-неоднородного полупространства, пьезоэлектрического слоя на диэлектрическом полупространстве. С помощью преобразования Фурье задача сводится к системе обыкновенных дифференциальных уравнений (ОДУ). Для получения аналитического представления решения системы, рассматривается аналогичная вспомогательная задача с известным распределением нормальной компоненты электрической индукции D2*(x1). Решение полученной системы ОДУ строится традиционным способом. Возвращаясь в пространство оригиналов, можно выписать выражения для полей смещения и потенциала в любой точке подложки через распределение напряжений и индукции на поверхности.

Например, для сдвиговой задачи они имеют вид:

Здесь T32*(x1) и D2*(x1) – распределение механических напряжений и электрической индукции на поверхности, Knm – компоненты матрицы функции Грина, описывающие структуру и свойства материала подложки. Выражения для всех рассмотренных случаев получены в аналитическом виде.

Выражение для потенциала из (5) можно использовать в качестве интегрального уравнения (ИУ) для определения D2*(x1). В случае отсутствия механических напряжений на поверхности, ИУ будет иметь вид:

Однозначная разрешимость уравнения (6) определяется свойствами функций Knm() – они имеют счетное количество чередующихся нулей и полюсов на вещественной оси и убывают при как O.

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

представляется в виде:

K соответственно. Остальные z k k n1 1,..., N и k k n2 1,..., N ), N maxn1, n2 – комплексные полюсы и нули функции K, лежащие в полосе Im 0. Далее используется техника, предложенная в цикле работ В.А. Бабешко и О.Д. Пряхиной, позволяющая свести уравнение (7) к уравнению относительно неизвестной функции, t m, которая связана с решением исходного уравнения (7) соотношением Для решения уравнения (10) в диссертации используется численный метод БубноваГалеркина. Сочетание метода фиктивного поглощения с численными методами позволило обобщить его на среды типа многослойного полупространства и исследовать тонкие вопросы влияния начальных напряжений на динамику пьезоактивной неоднородной среды.

При использовании приближенного метода функция распределения заряда qn ( x1 ) под электродом с номером n представляется в виде разложения по ортонормированной системе функций k ( x1 ), k 1,2,..., M, определенных на отрезке an, bn, который совпадает с положением самого электрода:

Используя представление (11), интегральное уравнение (6) сводится к решению системы линейных алгебраических уравнений:

Характерный вид распределения D2*(x1) для многоэлектродной структуры при однополярном режиме возбуждения приведен на рис. 2. Влияние начальных напряжений на решение интегрального уравнения (6) носит сложный характер и проиллюстрировано графиками на рис. 3 и 4.

0. 0. -0. -0. -0. Рис. 2. Распределение плотности заряда для системы 9-ти электродов на поверхности пьезоэлектрического полупространства. Сплошная линия – вещественная часть, пунктирная – мнимая.

Рис. 3. Влияние начальных напряжений на распределение плотности заряда для системы 9-ти электродов (абсолютное изменение). Сплошная линия – вещественная часть, пунктирная – мнимая.

Рис. 4. Влияние начальных напряжений на распределение плотности заряда для системы 9-ти электродов (относительное изменение). Сплошная линия – вещественная часть, пунктирная – мнимая.

В третьей главе рассматриваются результаты решения уравнения (6) для распределения полей многоэлектродных систем и влияния краевых эффектов на характеристики устройств были решены задачи для одного электрода, 3-х, 5-ти, 7-ми, 9-ти электродных структур и так далее в сторону увеличения количества электродов.

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

Влияние краевых эффектов, приводящее к существенному изменению формы распределения зарядов, было выявлено только для 3-4 крайних электродов.

Прямое решение уравнения (6) для многоэлектродных (более 30) систем связано подлежащих определению неизвестных, при построении системы уравнений (12) было сделано допущение, что распределение зарядов под внутренними электродами совпадает с распределением зарядов под центральным электродом:

Учитывая равенства (13) систему уравнений (12) можно переписать в виде:

Анализ возможности замены решения строгой задачи для большого количества электродов на решение задачи по упрощенной схеме выявил следующий недостаток: в месте стыковки решения для крайних электродов с периодической структурой из внутренних электродов возникала ошибка, растущая с увеличением количества электродов в системе и достигающая 40% уже для 25 электродов. Указанная ошибка электродами и с некоторым изменением абсолютных значений qn ( x1 ) по мере приближения к краю электродной структуры. Для уменьшения разницы между точным нецентральными электродами были дополнительно домножены на весовую функцию w( x1 ). В качестве весовой функции была выбрана сплайновая аппроксимация зависимости отклонения распределения заряда под электродом от распределения Рис. 5. Отклонение распределения заряда под электродами по сравнению с центальным электродом, построенные в относительных координатах для некоторых электродных структур.

заряда под центральным электродом координатах x структурой подложки и не зависит от количества электродов.

Окончательно, предлагаемая схема расчета многоэлектродных систем сводится к последовательному решению двух вспомогательных задач – для одного и семи электродов. Учет 3-х крайних электродов позволяют хорошо описать краевые эффекты, а периодическая система между ними взвешивается согласно функции распределения заряда w( x1 ) под всей электродной структурой. Вид функции w( x1 ) определяется из решения контактной задачи для одного электрода. Ошибка расчета электрического поля по такой упрощенной схеме не превышает 3–4 % и обусловлена в основном несимметричностью влияния электродной системы на нецентральные электроды. При этом вычислительные затраты на построение решения могут быть уменьшены в десятки раз в зависимости от количества электродов в моделируемой системе.

Рис. 6. Ошибка расчета электрического поля электродной структуры по упрощенной первоначальной модели, сплошная – для модифицированной.

использования приближенного метода с использованием весовой функции и без нее.

Нетрудно видеть, что использование весовой функции в значительной степени повышает эффективность метода.

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

Так как диэлектрическая проницаемость пленки намного выше, чем материала подложки, электрическое поле в основном локализуется в слое между электродами. Это приводит к увеличению электрической индукции в зазоре по мере уменьшения толщины слоя, структура распределения становится более однородной, что иллюстрируется графиками на рис. 7.

На границе раздела «пьезоэлектрик-диэлектрик» при уменьшении толщины происходит увеличение касательных механических напряжений T32, что влечет за собой появление нового максимума нормальной компоненты электрической индукции D2.

Рис. 7. Изменение структуры распределения нормальной компоненты электрической индукции D2 при уменьшении толщины пьезоэлектрического слоя.

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

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

Рис. 8. Изменение скорости ПАВ при учете деформаций, характерных для тонких сегнетоэлектрических пленок, в рамках линеаризованной теории электроупругости и в рамках феноменологической теории.

С другой стороны, фиксированные деформации среды для сегнетоэлектрической пленки позволяют использовать линеаризованную около конечного начального деформированного состояния теорию электроупругости. Деформации в среде определялись по изменению параметров решетки в пленке по сравнению с объемным образцом. Проведено сравнение результатов расчета фазовых скоростей ПАВ с использованием констант, полученных из феноменологической теории и констант, рассчитанных исходя из линеаризованной теории динамики преднапряженных электроупругих сред.

На рис. 8 приведена зависимость изменения фазовых скоростей ПАВ V f (V f V f ) *1000 от безразмерной частоты для различных наборов материальных параметров: a) объемный образец BST08 на подложке MgO; b) пленка BST08/MgO, h < 40 nm (линеаризация около НДС); c) пленка BST08/MgO, h > 105 nm (линеаризация около НДС); d) пленка BST08/MgO (феноменологическая теория с учетом misfit strain).

Константы для объемного образца твердого раствора BST08 в естественном состоянии также были взяты из феноменологической теории. При наличии сжимающих в плоскости подложки напряжений (c, d) фазовая скорость ПАВ увеличивается, при наличии растягивающих напряжений (b) – уменьшается. В области низких частот эффект изменения скорости незначителен, но с ростом частоты изменение скорости становится существенным. Точка выхода второй моды * смещается в сторону высоких частот по мере увеличения сжимающих напряжений.

Полученные результаты позволяют сделать вывод об адекватности предложенных в работе подходов к описанию динамических процессов в пьезоэлектрических структурах с покрытием из тонких сегнетоэлектрических пленок. Изменение динамических характеристик в обоих моделях под действием сжимающих напряжений качественно совпадает. Окончательно выбор расчетной модели можно будет сделать только после проверки на лучшее совпадение с экспериментальными данными.

В заключении формулируются основные результаты, полученные в диссертации, обсуждаются возможные направления продолжения исследований.

Результаты, выносимые на защиту

:

1) Разработан численно-аналитический метод решения связанных краевых задач, позволяющий с высокой точностью учитывать динамические свойства среды и рассматривать структурно неоднородные полуограниченные структуры в условиях начального напряженного состояния. Разработаны соответствующие численные алгоритмы и программное обеспечение, реализующие названный 2) Исследовано влияние начальных напряжений на область существования и скорости ПАВ, на структуру механических и электрических полей, на интегральные динамические характеристики для пьезоактивных подложек различной конфигурации.

3) Исследованы особенности распределения электрических полей многоэлектродных систем на поверхности неоднородного пьезоактивного полупространства в различных режимах возбуждения.

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

5) Исследовано влияние толщины пьезоэлектрического слоя на структуру механических и электрических полей в подложке и на интегральные характеристики устройства.

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

1. Белянкова Т.И., Лыжов В.А. Некоторые особенности динамики слабо неоднородных пьезоактивных структур // Вестник ЮНЦ РАН. 2010. Т. 6. № 2.

С. 3-10.

2. Белянкова Т.И., Калинчук В.В., Лыжов В.А. Связанная смешанная задача для системы электродов на поверхности преднапряженного электроупругого структурно неоднородного полупространства // ПММ. 2010. Т. 74.Вып. 6. С. 897– 3. Шейдаков Д.Н., Лыжов В.А., Михайлова И.Б. Потеря устойчивости сжатого пористого стержня с однородным покрытием // Вестник ЮНЦ РАН. 2011, Т.7, № 1, С.5-12.

4. Белянкова Т.И., Калинчук В.В., Лыжов В.А. Роль размерных параметров в формировании волновых полей в неоднородных пьезоактивных структурах // Вестник ЮНЦ РАН. 2011, Т.7, № 3, С.3-12.

5. Белянкова Т.И., Лыжов В.А. Моделирование преднапряженных функционально градиентных материалов и покрытий // Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2019– пьезоактивного полупространства в условиях неоднородного начального напряженного состояния // Труды X международной конференции «Современные проблемы механики сплошной среды». Ростов-на-Дону, 5-9 декабря 2006 г.

Ростов-на-Дону: Издательство ООО «ЦВВР», 2007. Т. 2. С. 59-63.

2. Белянкова Т.И., Лыжов В.А. Об учете констант III порядка в задачах динамики электроупругих тел // Труды XI международной конференции «Современные проблемы механики сплошной среды». Ростов-на-Дону, 26-29 ноября 2007 г.

Ростов-на-Дону: Издательство ООО «ЦВВР», 2008. Т. 2. С. 27-31.

3. Лыжов В.А., Тукодова О.М., Ворович Е.И. Взаимодействие многоэлектродных структур с пьезоактивными средами // Современные проблемы механики сплошной среды. Труды XIII международной конференции, г. Ростов-на-Дону, 13–15 октября 2009 г. Ростов-на-Дону: Издательство ООО «ЦВВР», 2009. Т. 2. С.

138–142.

4. Лыжов В.А., Тукодова О.М., Ворович Е.И. Взаимодействие многоэлектродных структур с предварительно напряженными электроупругими средами // Современные проблемы механики сплошной среды. Труды XIV международной конференции, Азов, 19–24 июня 2010 г. Ростов-на-Дону: Издательство ЮФУ, 2010. Т. 1. С. 211– преднапряженной электроупругой средой // VI научная конференция студентов и аспирантов базовых кафедр ЮНЦ РАН: тезисы докладов, Ростов-на-Дону, 19– апреля 2010 г. Ростов-на-Дону: Издательство ЮНЦ РАН, 2010. С. 389.

многоэлектродных структур // Тезисы докладов международной конференции «Современные проблемы механики», посвященной 100-летию Л.А. Галина.

Москва, 2012 г. С. 51-52.

7. Лыжов В. А., Тукодова О. М., Ворович Е. И., Агаян К.Л. Эффективный приближенный метод построения связанных полей многоэлектродных структур // Современные проблемы механики сплошной среды. Труды XVI международной конференции, г. Ростов-на-Дону, 16-19 октября 2012 г. Ростов-на-Дону:

Издательство ЮФУ, 2012. Т. 2. С. 129– 8. Лыжов В.А. Динамические свойства наноразмерных сегнетоэлектрических пленок // Труды VII Всероссийской (с международным участием) конференции по механике деформируемого твердого тела (г. Ростов-на-Дону, 15–18 октября 2013 г.). Ростов-на-Дону: Издательство ЮФУ, 2013. Т. 2. С. 76-80.