Surface SH-Waves in Pre-Stressed Piezoelectrics with Functionally Graded

Abstract


A model of a ferroelectric structure consisting of a homogeneous piezoactive half-space with an inhomogeneous coating which is either a layer or a packet of homogeneous or functionally graded piezoactive layers is suggested. It is assumed that the half-space, as well as the coating, being piezoelectrics of the hexagonal syngony of a 6 mm class in their intrinsic state, are in the conditions of the action of the initial mechanical stresses. Investigations of dynamic properties of the functionally oriented pre-stressed structures are carried out in Lagrangian rectangular coordinate system. Linearized constitutive relations and motion equations are used. The boundary problem for the system of differential equations in partial derivatives is reduced to the system of ordinary differential equations by means of the operational calculus. In case of the structure homogeneous components, Green’s function is constructed analytically in a closed form on the basis of solving the system of differential equations with constant coefficients. In case of inhomogeneous (functionally graded) components, the system of differential equations with variable coefficients is reduced to the system of Cauchy initial problems by means of a special substitution. In this case, Green’s function is constructed numerically on using Runge-Kutta numerical methods with Merson modification, which allows us to control the error of calculations in an effective way. When constructing Green’s function of the ferroelectric structure with an inhomogeneous coating, we used the matrix approach that allows us to combine analytical and numerical methods of constructing its isolated components. The influence of the type and the value of the initial stresses on the peculiarities of the surface wave propagation in heterostructures is studied. It became possible to determine the conditions under which the action of the initial mechanical stresses leads to the increase of Bleustein-Gulyaev wave velocity with respect to the velocity of the original material, as well as the conditions under which the piezoelectric structure ceases to be a weakly inhomogeneous one.

Full Text

Введение Проблема распространения поверхностных акустических волн является одной из актуальных проблем современной механики. С одной стороны, это связано с широким спектром приложений в акустике, геофизике, сейсмологии, дефектоскопии и др., с другой стороны, с бурным развитием технологий создания новых пьезоэлектрических материалов с уникальными, ранее не наблюдавшимися в природе свойствами. Возможность создания на их основе принципиально новых типов электронных приборов и устройств обусловило значительный рост работ, посвященных математическому моделированию процессов распространения волн в функционально-градиентных средах [1-6]. В [1] построена общая теория электродинамики сплошной среды, приведены основные типы термодинамического потенциала, предложены подходы к построению линеаризованных соотношений динамики преднапряженной электроупругой среды. В [2-4] приведены определяющие соотношения динамики сплошной пьезоактивной среды при наличии начальных механических или электростатических полей. В [5-7] дан вывод определяющих соотношений динамики упругой и пьезоактивной среды в произвольной, в общем случае криволинейной системе координат. Особенности распространения сдвиговых волн в слоистых преднапряженных средах рассмотрены в [8-12]. В [8-11] изучено влияние начальных напряжений на скорости распространения волн Лява и Гуляева-Блюштейна. В [12] изучены вопросы распространения волн в слабо неоднородных средах, неоднородность которых обусловлена наличием начальных напряжений или поляризации. Показано, что при специальных условиях в слабо неоднородном полупространстве могут возникать вторые моды поверхностных волн типа волн Гуляева-Блюштейна. В [13-16] исследуются процессы распространения волн в функционально-градиентных пьезоактивных средах, параметры неоднородности которых допускают построение аналитического решения в той или иной форме. В [13] исследованы особенности распространения волн Гуляева-Блюштейна в среде, все характеристики которой изменяются по одному и тому же закону. Рассмотрены случаи полиномиальной и экспоненциальной зависимостей. Исследовано влияние параметра неоднородности на скорость, дисперсию, коэффициенты затухания, глубину профилей, электромеханический коэффициент связи. В [14] исследованы особенности распространения волн Гуляева-Блюштейна в среде, характеристики которой изменяются по специальному, допускающему аналитическое решение закону. Изучено влияние параметра неоднородности на характеристики волнового поля. В [15] изучаются волны Лява в среде из поляризованной керамики с функционально-градиентным покрытием. Для построения аналитического решения используется асимптотический подход. Исследовано влияние параметра неоднородности на дисперсионные соотношения, в частности на фазовые скорости поверхностных волн, а также на электромеханический коэффициент связи. В [16] изучается влияние параметра неоднородности и начальных напряжений на характеристики волнового поля. В настоящей работе предложена модель сегнетоэлектрической структуры, состоящей из однородного пьезоактивного полупространства с неоднородным покрытием, представляющим собой либо слой, либо пакет однородных или функционально-градиентных пьезоактивных слоев. Предполагается, что полупространство, равно как и покрытие, находятся в условиях воздействия начальных механических напряжений, свойства покрытия описываются произвольными функциями. Для этого случая система дифференциальных уравнений с переменными коэффициентами специальной заменой сводится к системе начальных задач Коши. Функция Грина строится численно на основе использования процедур Рунге-Кутты с модификацией Мерсона, которая позволяет эффективно контролировать погрешность вычислений. При построении функции Грина сегнетоэлектрической структуры с неоднородным покрытием использован матричный подход, позволяющий сочетать аналитические и численные методы построения отдельных ее составляющих. В рамках линеаризованной теории электроупругости в работе исследованы особенности распространения поверхностных волн в сегнетоэлектрических гетероструктурах с учетом различного характера механических воздействий на составляющие структуры. Изучено влияние величины наведенной деформации, вида и характера преднапряжений в составляющих структуры на особенности поведения и скорость распространения поверхностных акустических волн (ПАВ) для ряда наиболее востребованных сегнетоэлектриков. Установлены условия, при которых действие начальных механических напряжений приводит к увеличению скорости волны Гуляева-Блюштейна относительно скорости исходного материала, а также условия, при которых пьезоэлектрическая структура перестает быть слабо неоднородной. 1. Постановка задачи Рассматривается задача о распространении сдвиговых горизонтально поляризованных волн, движущихся в направлении по поверхности составной преднапряженной пьезоактивной среды. Полагаем, что колебания среды вызваны действием удаленного источника гармонических колебаний, среда представляет собой однородное преднапряженное полупространство , с преднапряженным покрытием (рис. 1 а, б). Рис. 1. Геометрия задачи Fig. 1. Geometry of the problem Покрытие моделируется либо одним слоем (см. рис. 1, б), либо пакетом , как однородных, так и неоднородных слоев из функционально-градиентного материала. В качестве исходного материала структуры используется пьезоэлектрик класса 6 mm гексагональной сингонии, ось симметрии которого в естественном состоянии (ЕС) совпадает с осью , векторы поляризации составляющих покрытия и полупространства совпадают либо противоположны по направлению. Начально-деформированное состояние (НДС) каждой из составляющих структуры однородно и наводится за счет действия начальных механических напряжений, внешние начальные электрические воздействия отсутствуют [1-7, 17]: (1) Здесь - радиус-векторы точки среды в начально-деформированном и естественном состоянии соответственно; - относительные удлинения волокон, направленных в естественной конфигурации вдоль осей, совпадающих с декартовыми координатами; - символ Кронекера. Неоднородность начального напряженного состояния составляющих структуры вызвана только неоднородностью физических свойств. Исследования проводятся в лагранжевой системе координат, совпадающей с прямоугольной декартовой системой, режим колебаний полагается установившимся, динамический процесс удовлетворяет условиям (2) При этом для слабо неоднородной пьезоэлектрической структуры (см. рис. 1, а, б) в ЕС выполняются соотношения: для однородных составляющих [6, 7] либо (3) для неоднородных составляющих [12] (4) Далее используем безразмерные параметры: линейные параметры отнесены к характерной толщине k-го слоя , плотность - к плотности полупространства , упругие параметры - к модулю сдвига подстилающего полупространства . При переходе к безразмерным параметрам пьезоэлектрических и диэлектрических констант используется множитель , при этом , ( - диэлектрическая проницаемость вакуума). Далее штрихи опускаем. В качестве безразмерной частоты используется либо параметр , либо (где - скорость сдвиговой волны полупространства и - скорость сдвиговой волны полупространства с учетом пьезоэлектрических свойств в ЕС [18]). В рамках принятых предположений краевая задача о колебаниях преднапряженной электроупругой среды описывается линеаризованными уравнениями [6-7, 12]: (5) (6) для вакуума (7) с граничными условиями на поверхности: отсутствие механических воздействий: (8) электрически свободная поверхность: (9) металлизированная поверхность: (10) на границе раздела сред ( ) (11) на бесконечности , (12) Здесь - оператор Гамильтона; - расширенный вектор перемещений; - вектор внешней нормали к поверхности среды определены в системе координат, связанной с естественным состоянием (ЕС); - плотность материала n-составляющей в ЕС; - оператор Лапласа. Верхний индекс « » соответствует номеру среды, составляющей структуру ( - вакуум, - поверхностный слой, - полупространство). Линеаризованные тензор напряжений и вектор индукции с учетом (1) представляются в виде [12] (13) где (14) Участвующие в (14) компоненты тензора Кирхгофа , а также упругие константы зависят как от свойств материала, так и от вида начального напряженного состояния среды [1, 6, 7, 12]: (15) Здесь использована упрощенная запись термодинамического потенциала [1, 6, 7]: где - компоненты тензора деформаций Коши в случае однородной деформации (1) определены формулами где - компоненты материального вектора электрической индукции; и - упругие константы II и III порядка; - константы диэлектрической восприимчивости, связанные с константами диэлектрической проницаемости соотношениями Для удобства дальнейшего изложения введем обозначения (16) Подставляя выражения компонент тензора напряжений и вектора индукции (13) с учетом (2), свойств (3), (4), обозначений (16) в (5)-(12) для преднапряженных однородных составляющих структуры, получаем (индексами после запятой отмечены производные по соответствующим координатам) (17) для преднапряженных функционально градиентных составляющих (18) для вакуума (19) граничные условия (20) (21) (22) (23) (24) В рамках настоящей работы исследуем две краевых задачи о распространении поверхностных горизонтально поляризованных волн в преднапряженных пьезоэлектриках с покрытием [1, 6-12]: задача I - со свободной поверхностью, описывается уравнениями движения (17)-(19) с граничными условиями (20), (21), (23), (24); задача II - с металлизированной поверхностью, описывается уравнениями движения (17)-(19) с граничными условиями (20), (22)-(23). 2. Решение краевых задач. Дисперсионное уравнение для преднапряженного пьезоактивного полупространства с преднапряженным покрытием При решении краевых задач I и II используем преобразование Фурье по координате , - параметр преобразования. Движение пьезоэлектрической структуры в трансформантах Фурье имеет вид: для преднапряженных однородных составляющих (25) для преднапряженных функционально градиентных составляющих (26) для вакуума (27) с граничными условиями (28) (29) (30) (31) , (32) Решение задач (25)-(27) с граничными условиями (28)-(32) в трансформантах Фурье для однородных составляющих покрытия ищем в виде ( ) [6, 7, 12] (33) для неоднородных составляющих (34) для полупространства и вакуума (35) Участвующие в представлениях (33), (35) удовлетворяют характеристическому уравнению (36) Коэффициенты определяются из решения однородной системы линейных уравнений с матрицей коэффициентов . в (34) - линейно независимые решения задачи Коши с начальными условиями для уравнения (37) (38) Следует отметить, что при построении функции Грина для электроупругой среды с неоднородным покрытием могут быть использованы различные и весьма эффективные подходы [19-23]. В настоящей работе использован предложенный в [24] подход, сочетающий аналитические решения (33), (35) для однородных с численными методами для неоднородных составляющих структуры, в частности, при решении систем (37), (38) применен модифицированный метод Рунге-Кутты. Неизвестные (33)-(35) определяются из удовлетворения граничным условиям (28), (29), (31), (32) для задачи I или (28), (30)-(32) для задачи II. Следуя [24] представим дисперсионное уравнение задач в виде (39) (40) Размеры матрицы и матриц ее составляющих определены геометрией задачи и условиями на поверхности. Задача I: матрицы и - квадратные, размер определяется геометрией задачи и равен соответственно и . , и - матрицы размера , и , вид которых полностью определяется свойствами верхнего слоя, подстилающего полупространства и диэлектрической проницаемостью вакуума соответственно и не зависит от свойств и количества возможных промежуточных слоев. (41) для преднапряженного однородного поверхностного слоя (42) для преднапряженного функционально градиентного поверхностного слоя (43) использованы обозначения (44) В случае - двухслойное полупространство (см. рис. 1, б) вид матрицы определен свойствами слоя, при этом для преднапряженного однородного слоя (45) для преднапряженного функционально градиентного слоя совпадает с диагональной единичной матрицей размера . В случае - полупространство со слоистым покрытием матрицы и принимают вид (46) (47) Участвующие в (46) матрицы и определены свойствами слоев с общей границей , имеют размер , верхний индекс отвечает номеру слоя, аргумент - границе раздела (см. рис. 1, а). Для преднапряженного однородного слоя матрица с учетом обозначений (44) определена формулой (48) для преднапряженного функционально градиентного слоя (49) Задача II: матрица имеет размер , в зависимости от геометрии задачи может быть представлена выражениями (40) и (47). Размер и вид матрицы совпадает с соответствующей матрицей задачи I, определяется формулами (45) и (46). Матрицы , и , участвующие в (40), имеют размер , и . Вид матриц определяется из представлений (41)-(43), путем вычеркивания в выражении 2-й строки и 3-го столбца - матрица становится нулевой, в представлении 3-го столбца и в представлениях (42), (43) 2-й строки. Размер и вид матриц и полностью совпадает с соответствующими матрицами задачи I (48), (49). 3. Определение НДС В рамках сделанных предположений начально деформированное состояние в каждой из составляющих структуры однородно (1) и может различаться по характеру и величине начальных деформаций. Неоднородность начального напряженного состояния составляющих вызвана только неоднородностью их физических свойств. Характер начального напряженного состояния определяется тензором Кирхгофа, компоненты которого (15) с учетом свойств материала и предположений об однородности НДС, наведенного за счет механических воздействий в отсутствие внешних начальных электростатических полей определены формулами [6, 7, 12] (50) Здесь коэффициенты позволяют учитывать влияние модулей упругости III порядка (сумма по немым индексам): . В силу предположения об отсутствии начальных внешних электрических воздействий и в зависимости от способа задания НДС неизвестными могут быть либо , либо, например, определяемые из решения системы (50). В работе исследовано влияние одноосных НДС: ; двухосных НДС: и гидростатического НДС: . Исследования проводились для структур из монокристаллических ZnO [25], CdSe [26] и поликристаллических сегнетоэлектриков, таких как пьезокерамики Soft - PZT DL-61HD, Hard PZT DL-40 [27]. Физические параметры сред приведены в таблице, = 8,85·10-12 Ф/м, - коэффициент электромеханической связи. Относительные скорости волны Гуляева-Блюштейна (ВГБ) для полупространства со свободной и металлизированной поверхностью рассчитаны без учета модулей упругости III порядка. Физические свойства материалов в ЕС Physical properties of materials in natural state (NS) ZnO CdSe PZT DL61 PZT DL40 ρ [кг/м3] 5680 5504 8200 7700 с11 [н/м2] 2,09718·1011 7,4057·1010 14,6000·1010 17,8000·1010 с12 [н/м2] 1,2114·1011 4,5155·1010 9,6000·1010 10,1000·1010 с13 [н/м2] 1,0513·1011 3,9302·1010 10,0000·1010 9,2000·1010 с33 [н/м2] 2,10941·1011 8,3551·1010 13,0000·1010 12,4000·1010 с44 [н/м2] 4,2449·1010 1,3167·1010 3,9000·1010 2,3000·1010 е15 [Кл/м2] -0,59 -0,138 33,10 6,20 е31 [Кл/м2] -0,61 -0,159 -15,80 -0,10 е33 [Кл/м2] 1,14 0,347 25,30 9,00 ε11/ε(0) 8,3 9,33 2810,0 290,0 ε33/ε(0) 8,8 10,2 2520,0 210,0 VP [м/с] 6097 3668 4219,58 4808 VS [м/с] 2743 1547 2180,85 1728,3 VSe [м/с] 2892 1560 3182,58 2220,8 0,9999417 0,9999986 0,999999982 0,999999082 0,994944441 0,99985181 0,847723696 0,918947544 0,01 0,0003 0,28 0,1555 4. Численные результаты Для детального изучения влияния величины наведенной деформации, вида и характера преднапряжений в составляющих структуры на особенности поведения и скорость распространения поверхностных акустических волн (ПАВ) рассмотрим простейший случай преднапряженного пьезоэлектрического полупространства с преднапряженным однослойным покрытием (M = 2, рис. 1, б). В качестве материала структуры (покрытия и основания) используем гексагональные моно- и поликристаллические сегнетопьезоэлектрики (таблица) с различным коэффициентом электромеханической связи, без учета модулей упругости III порядка. Полагаем, что ось симметрии материала в ЕС совпадает с осью , векторы поляризации составляющих покрытия и полупространства совпадают, либо противоположны по направлению (3). На рис. 2-8 показано влияние характера и величины начальных механических воздействий на относительные фазовые скорости поверхностных волн ( , где , - решение дисперсионного уравнения (39) с матрицей (40), (47) в зависимости от геометрии задачи в обозначениях (41)-(46), (48)-(49)). Цифрами 0, 1, 2, …, 10 на рисунках отмечено сочетание наведенных деформаций - 1/1, 0,97/1, 1,03/1, Рис. 2. Влияние величины, характера и типа начальных воздействий на относительные скорости ВГБ для структуры ZnO/ZnO Fig. 2. Influence of value, nature and type of the initial impacts on relative rates of Bleustein-Gulyaev wave for ZnO / ZnO structure 0,95/0,97, 0,97/0,97, 1/0,97, 1,03/0,97, 0,97/1,03, 1/1,03, 1,03/1,03, 1,06/1,03 соответственно, при этом нижние индексы «f» и «m» обозначают свободную и металлизированную поверхность, верхний индекс «e» случай . На рис. 2-4 показано влияние величины и характера начальных напряжений, действующих на составляющие структуры в рамках одного типа НДС, на относительные фазовые скорости ПАВ. Цифрами на рисунках отмечено сочетание деформаций, наведенных действием начальных напряжений в структурах из монокристаллических сегнетоэлектриков ZnO/ZnO (см. рис. 2), CdSe/CdSe (см. рис. 3) и поликристаллических сегнетопьезокерамик PZT DL-61/PSZ DL-61, PZT DL-40/PZT DL40 (см. рис.4). Поляризация в слое и полупространстве совпадает, , сплошными линиями на рисунках обозначены скорости ПАВ задачи I - со свободной поверхностью, штриховыми линиями - задачи II с металлизированной поверхностью. На рис. 2 в, 3 в и 4 а, б показано сравнение влияния одноосных 1х1 (сплошные линии) и 1х3 (пунктирные линии) начальных напряжений на ПАВ в рассматриваемых структурах с металлизированной поверхностью. Рис. 3. Влияние величины, характера и типа начальных воздействий на относительные скорости ВГБ для структуры CdSe/CdSe Fig. 3. Influence of value, nature and type of the initial impacts on relative rates of Bleustein-Gulyaev wave for CdSe/CdSe structure Из сравнения рис. 2 и 3 видно, что при малом значении КЭМС ( ) различие в зависимости фазовых скоростей ПАВ от частоты для задач со свободной и металлизированной поверхностью незначительно (см. рис. 3, а, в). Влияние начальных напряжений, наводящих одноосные НДС в составляющих структуры, минимально в случае 1х2. При 1х1 и 1х3 - влияние соизмеримо и зависит от свойств материала и величины наведенных деформаций. Следует отметить, что воздействие начальных механических напряжений может существенно менять свойства структуры. Так, при определенном соотношении наведенных деформаций (кривые 1, 3, 8 рис. 2, а и 2, 5 рис. 2, г, д) сегнетоэлектрическая структура остается слабо неоднородной: по поверхности распространяется волна Гуляева-Блюштейна (ВГБ), скорость которой в зависимости от характера воздействий может быть больше или меньше значения скорости в ЕС, диапазон изменения незначительный. Соотношение наведенных деформаций 0,97/1,03 при НДС 1х1 (кривые 7 на рис. 2 а, в и рис. 3 а, в) приводит к увеличению диапазона изменения скорости, появляются вторые, а при определенных типах НДС и более высокие моды. Структура перестает быть слабо неоднородной, по поверхности распространяются сдвиговые горизонтально поляризованные волны Лява. Соотношение деформаций 1,03/1, 1/0,97, 1,03/0,97, 1,06/1,03 в рамках Рис. 4. Влияние НДС на относительные скорости ВГБ для структур PZT DL-61/PZT DL-61 и PZT DL-40/PZT DL-40 Fig. 4. Influence of IDS on relative rates of Bleustein-Gulyaev for PZT DL-61/PZT DL-61 and PZT DL-40/PZT DL-40 structures НДС 1х1 приводит к увеличению скорости ВГБ до значения , волна существует в ограниченном частотном диапазоне, который определяется величиной КЭМС, соотношением деформаций и характером начальных воздействий (см. рис. 2, б, е, рис. 3, б). На рис. 4, а-е показано влияние величины, характера и типа наведенных НДС на фазовые скорости ПАВ в структурах из сегнетопьезокерамик PZT DL-61 (см. рис. 4, а, в, д) и PZT DL-40 (см. рис. 4, б, г, е) с большим значением КЭМС (таблица). Сравнение одноосных НДС 1х1 и 1х3 дано на рис. 4, а, б для первых мод ПАВ. В случае структуры из PZT DL-40 в рассматриваемом частотном диапазоне, начиная с = 6,7, появляются вторые моды ПАВ для всех приведенных соотношений деформаций как при НДС 1х1, так и при НДС 1х3. В случае PZT DL-61 вторые моды в рассматриваемом диапазоне появляются только при НДС 1х3. Из сравнения рис. 2, 3 и 4 видно, что действие начальных напряжений, наводящих в структуре одноосные НДС 1х1 или 1х3 соизмеримы по влиянию на скорости ПАВ. Для многих пьезоэлектриков преобладающее влияние имеет НДС 1х1, однако для ряда сегнетопьезокерамик с высоким значением пьезоэлектрических и диэлектрических модулей в ЕС влияние НДС 1х3 может значительно превзойти влияние НДС 1х1. Для пьезокерамик с большим значением КЭМС в случае задачи с металлизированной поверхностью характерно существование поверхностной волны во всем частотном диапазоне даже при соотношении наведенных деформаций, приводящих к увеличению фазовой скорости относительно значения скорости ВГБ соответствующего полупространства (рис. 4, д, е). Рис. 5. Влияние величины, характера и типа НДС на относительные скорости ВГБ для структуры ZnO/ZnO Fig. 5. Influence of value, nature and type of IDS on relative rates of Bleustein-Gulyaev for ZnO/ZnO structure На рис. 5, а-г показано влияние направления поляризации составляющих структуры ZnO/ZnO (см. рис. 5, а), величины, характера наведенных деформаций и типа начальных напряжений (см. рис. 5, б-г) на изменение скорости ВГБ для задач со свободной и металлизированной поверхностью. Предполагается, что наведенные в составляющих структуры деформации равны. На рис. 5, а приведено характерное поведение скорости ВГБ в ЕС при и . В случае НДС 1х1 при соотношении наведенных деформаций 1,03/1,03 (кривые 9 рис. 5, б) значение относительных скоростей и достигает значений 1,0356454 и 1,30791675 соответственно. Из рисунков видно, что изменение направления вектора поляризации при всех прочих равных параметрах приводит к изменению характера частотой зависимости скорости ПАВ, среда перестает быть однородной, наблюдается дисперсия свойств (см. рис. 5, а, б). В рамках действия различных типов НДС в структуре ZnO/ZnO при условии наибольшее значение скорости и достигают при НДС 1х1 (кривые 9, рис. 5, б, в) и НДС 2х2 (см. рис. 5, г). На рис. 6, а-г и 7, а-г показано влияние на относительные фазовые скорости ПАВ величины и характера начальных напряжений, действующих на составляющие структуры ZnO/ZnO в рамках 1х1 (см. рис. 6) и 2х1 (см. рис. 7). Рис. а, б иллюстрируют влияние направления вектора поляризации в случае преднапряженного покрытия (рис. а) и в случае Рис. 6. Влияние на относительные фазовые скорости ПАВ соотношения величины и характера наведенных деформаций в структуре ZnO/ZnO в рамках НДС: 1х1 Fig. 6. Influence on relative phase rates of SAW of relation between value and nature of the induced deformations in structure ZnO/ZnO within IDS: 1х1 преднапряженной структуры (рис. б). Рис. в, г отражают влияние величины и характера преднапряжений на частотную зависимость фазовых скоростей ПАВ в случае . Цифрами на рисунках отмечено соотношение деформаций, наведенных действием начальных напряжений. Из сравнения рис. 2, б, е с рис. 6 и 7 видно, что эффекты, вызванные действием различного рода начальных механических напряжений, в пьезоэлектрической структуре существенно зависят от направленности векторов поляризации ее составляющих: в случае относительно меняется характер частотной зависимости скорости ПАВ, усиливается влияние преднапряжений. Рис. 7. Влияние на относительные фазовые скорости ПАВ соотношения величины и характера наведенных деформаций в структуре ZnO/ZnO в рамках НДС: 2х1 Fig. 7. Influence on relative phase rates of SAW of relation between value and nature of the induced deformations in structure ZnO/ZnO within IDS: 2х1 На рис. 8, а, б приведено влияние на фазовые скорости ПАВ начальных напряжений, наводящих в покрытии и основании НДС различных типов. Для удобства сравнения в качестве НДС основания использовано НДС 2х2, в качестве НДС покрытия одноосные НДС - 1х1, 1х2, 1х3 (см. рис. 8, а). На рис. 8, б дано сравнение одноосного 1х1 и двухосного 2х1 НДС покрытия. Соотношение наведенных деформаций на рис. 8, а, б отвечает случаю 1.03/0.97, векторы поляризации покрытия и основания противоположны по направлению. Из рис. 8, а, б видно, что при рассматриваемом соотношении деформаций независимо от НДС покрытия в рассматриваемом частотном диапазоне существуют вторые моды ПАВ, минимальное изменение скорости достигается при НДС 1х1/2х2 (см. рис. 8, а), максимальное - при НДС 2х1/2х2 (см. рис. 8, б). Сравнение рис. 6, а, б с рис. 3 в, г показывает Рис. 8. Влияние типов НДС покрытия и основания на фазовые скорости ПАВ Fig. 8. Influence of IDS types and foundations on SAW phase rates возможность увеличения скорости ПАВ за счет различий в начальных механических воздействиях на составляющие пьезоэлектрической структуры.

About the authors

T I Belyankova

Southern Scientific Center of Russian Academy of Sciences; Southern Federal University

Email: tbelen415@mail.ru

V V Kalinchuk

Southern Scientific Center of Russian Academy of Sciences; Southern Federal University

Email: kalin@ssc-ras.ru

References

  1. Maugin G.A. Continuum Mechanics of Electromagnetic Solids. - Amsterdam, Elsevier Science Publishers, 1991. - 621 p.
  2. Thurston R.N., Brugger K. Third-order elastic constants and the velocity of small amplitude elastic waves in homogeneously stressed media // Phys. Rev. - 1964 - Vol. 133. - No. 6A - P. A1604-A1610.
  3. Tiersten H.F. Electroelastic equations for electrode thin plates subject to large driving voltages // J. Appl. Phys. - 1993. - Vol. 74. - No. 5 - P. 3389-3393.
  4. Chai J.F., Wu T.T. Propagation of surface waves in a prestressed piezoelectric material // J. Acoust. Soc. Amer. - 1996. - Vol. 100. - No. 4. - Pt. 1 - P. 2112-2122.
  5. Калинчук В.В., Белянкова Т.И. Динамические контактные задачи для предварительно напряженных полуограниченных тел. - М.: Физматлит, 2008. - 240 с.
  6. Калинчук В.В., Белянкова Т.И., Евдокимова О.В. Определяющие соотношения динамики преднапряженной пьезоактивной среды в отсутствие внешних электрических полей // Вестн. Южного научного центра РАН. - 2006. - Т. 2, № 1. - С. 16-23.
  7. Евдокимова О.В., Белянкова Т.И., Калинчук В.В. Уравнения динамики преднапряженной пьезоактивной среды при наличии внешнего электростатического поля // Вестн. Южного научного центра РАН. - 2007. - Т. 3, № 4. - С. 19-25.
  8. Liu H., Wang Z.K., Wang T.J. Effect of initial stress on the propagation behavior of Love waves in a layered piezoelectric structure // Int. J.Eng Sci. - 2001. - Vol. 38. - P. 37-51.
  9. Jin F., Wang Z., Wang T. The Bleustein-Gulyaev (B-G) wave in a piezoelectric layered half-space // Int. J.Eng Sci. - 2001. - Vol. 39. - P. 1271-1285.
  10. Liu H., Kuang Z.B., Cai Z.M. Propagation of Bleustein-Gulyaev waves in a prestressed layered piezoelectric structure // Ultrasonics. - 2003. - Vol. 41. - P. 397-405.
  11. Love waves propagation in a piezoelectric layered structure with initial stresses / Z. Qian, F. Jin, Z. Wang, Xi’an China, K. Kishimoto // Acta Mechanica. - 2004. - Vol. 171. - P. 41-57.
  12. Белянкова Т.И., Лыжов В.А. Некоторые особенности динамики слабо неоднородных пьезоактивных структур // Вестн. Южного научного центра РАН. - 2010. - Т. 6, № 2. - С. 3-10.
  13. Collet B., Destrade M., Maugin G.A. Bleustein-Gulyaev waves in some functionally graded materials // European Journal of Mechanics A/Solids. - 2006. - Vol. 25. - Р. 695-706.
  14. Bleustein-Gulyaev waves in a functionally graded piezoelectric material layered structure / C. Xiaoshan, J. Feng, W. ZiKun, L. TianJia // Science in China Series G: Physics, Mechanics & Astronomy. - 2009. - Vol. 52. - No. 4. - Р. 613-625.
  15. Transverse surface waves on a piezoelectric material carrying a functionally graded layer of finite thickness / Z. Qian, F. Jin, Z. Wang, K. Kishimoto // International Journal of Engineering Science. - 2007. - Vol. 45 - P. 455-466.
  16. Effect of initial stress on Love waves in a piezoelectric structure carrying a functionally graded material layer / Z.-H. Qian, F. Jin, T. Lu, K. Kishimoto, S. Hirose // Ultrasonics. - 2010. - Vol. 50 - P. 84-90.
  17. Лурье А.И. Нелинейная теория упругости. - М.: Наука, 1980. - 512 с.
  18. Surface Acoustic Waves in Inhomogeneous Media / S.V. Biryukov, Y.V. Gulyaev, V.V. Krylov, V.P. Plessky. - New York, Springer-Verlag, 1995. - 287 p.
  19. Калинчук В.В., Белянкова Т.И., Богомолов А.С. К проблеме моделирования неоднородных материалов с заданными свойствами // Эколог. вестн. науч. центров Черноморского экономического сотрудничества. - 2006. - № 2. - С. 26-32.
  20. Численно-аналитическое построение матриц Грина трехмерных теорий упругости и электроупругости / Л.А. Игумнов, С.Ю. Литвинчук, В.П. Пазин, А.Н. Петров // Вестн. Нижегород. ун-та им. Н.И. Лобачевского. - 2010. - № 3-1. - С. 134-140.
  21. Igumnov L.A., Markov I.P., Rataushko Y.Y. Modeling the dynamics of 3-d elastic anisotropic solids using boundary element method // Advanced Materials Research. - 2014. - Vol. 1040. - P. 633-637.
  22. Balogun S., Achenbach J.D. Surface waves on a half-space with depth dependent properties // J. Acoust. Soc. Am. - 2012. - Vol. 132. - No. 3. - P. 1336-1345.
  23. Balogun S., Achenbach J.D. Surface waves generated by a line load on a half-space with depth-dependent properties // Wave Motion. - 2013. - Vol. 50 - P. 1063-1072.
  24. Белянкова Т.И., Калинчук В.В. К проблеме анализа динамических свойств слоистого полупространства // Акустический журнал. - 2014. - Т. 60, № 5. - С. 492-504.
  25. Akusticheskie kristally: Spravochnik. [Acoustic crystals: a Handbook]. Eds. Shaskolskaya M.P. - Moscow: Nauka Publishers, 1982. - 632 p.
  26. Sharma J.N., Pal M., Chand D. Propagation characteristics of Rayleigh waves in transversely isotropic piezothermoelastic materials // Journal of Sound and Vibration. - 2005. - Vol. 284. - P. 227-248.
  27. Material Specification Sheet. Available at: www.delpiezo.com/products.

Statistics

Views

Abstract - 302

PDF (Russian) - 61

Cited-By


PlumX


Copyright (c) 2016 Belyankova T.I., Kalinchuk V.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies