Full Text

Введение В различных технических устройствах, использующих принцип прямого и обратного пьезоэффекта, в качестве преобразователя энергии используются тонкие биморфные пластины [1-3]. Наиболее эффективной конструкцией, обладающей высокой механической прочностью и чувствительностью, является круглая металлическая подложка с наклеенными с двух сторон пьезокерамическими пластинами. Для исследования их напряженно-деформированного состояния, как правило, используется прикладная теория для тонких пластин, в которой гипотезы для тонкостенных упругих элементов дополняются аналогичными допущениями о характере изменения электрического поля [4-12]. Данный подход позволяет с помощью простых расчетных соотношений описать работу конструкции и подобрать ее геометрические размеры, позволяющие повысить эффективность трансформации энергии. Однако вопрос о характере изменения физических полей в пьезокерамических пластинах, особенно при исследовании систем с разрезными электродами, до настоящего времени изучен недостаточно полно. Поэтому для решения данной проблемы исследуется биморфная пластина в рамках теории электроупругости в трехмерной постановке. Новое замкнутое решение построено методом конечных интегральных преобразований. 1. Постановка задачи Пусть круглая жестко закрепленная пластина, занимающая в цилиндрической системе координат область : , состоит из металлической подложки толщиной и двух наклеенных на нее пьезокерамических элементов высотой , выполненных из материала гексагональной системы класса 6mm. Изгибные осесимметричные колебания возбуждаются за счет подведения к лицевым электродам пьезопластин с противоположным направлением вектора аксиальной поляризации, электрического напряжения и заземления подложки (рис.1). Рис. 1. Биморфная пластина Дифференциальные уравнения движения и электростатики, а также краевые условия в цилиндрической системе координат и безразмерной форме имеют следующий вид [13,14]: (1.1) ; , , , (1.2) , , ; , (1.3) , ; , (1.4) , , ; , , (1.5) , ; где ; , радиальная и аксиальная компоненты вектора перемещений; нормальные и касательные механические напряжения; потенциал электрического поля; , пьезомодули и коэффициенты диэлектрической проницаемости электроупругого материала объемная плотность и модули упругости пьезокерамического и упругого материалов; известные в начальный момент времени перемещения, скорости перемещений; , . Соотношения (1.1)-(1.5) представляют математическую формулировку рассматриваемой краевой задачи электроупругости. 2. Построение общего решения Решение осуществляется методом интегральных преобразований с использованием последовательно преобразования Ханкеля [15] с конечными пределами по переменной и обобщенного конечного преобразования (КИП) [16,17] по радиальной координате . Первоначально соотношения (1.1)-(1.5) приводятся к стандартной форме, позволяющей провести данную процедуру разделения переменных. Для этого последнее равенство (1.2) заменяется условием отсутствия касательных напряжений: , (2.1) и на основании теоремы о суперпозиции решений вводится новая функция , связанная с соотношением , (2.2) где дважды дифференцируемая функция, определяемая в процессе решения задачи из условия отсутствия вертикальных перемещений срединной поверхности пластины при *. В результате подстановки (2.2) в (1.1)-(1.5), (2.1) получаем новую краевую задачу относительно функций . При этом правые части второго дифференциального уравнения (1.1) и начальных условий (1.5) принимают следующий вид: , а условия (1.2) с учетом (2.1) при определяются равенствами** , , (2.3) К краевой задаче относительно применяем преобразование Ханкеля с конечными пределами по переменной , используя следующие трансформанты: , (2.4) и формулы обращения , (2.5) , , где - положительные нули функции расположенные в порядке их возрастания . В пространстве изображений получаем новую краевую задачу: (2.6) ; , (2.7) , ; , (2.8) , (2.9) где , . На втором этапе решения процедура стандартизации связана с приведением граничных условий (2.7),(2.8) к однородным при использовании следующих разложений: , (2.10) , где ; единичная функция Хэвисайда [20]; непрерывно дифференцируемые функции, удовлетворяющие следующим условиям*: , , (2.11) , , , , . В результате подстановки (2.10) в (2.6)-(2.9) получаем новую краевую задачу относительно функций . При этом правые части дифференциальных уравнений (2.6) и начальные условия (2.9) следует заменить на : , , , . Начально-краевую задачу (2.6)-(2.9) относительно функций решаем, используя структурный алгоритм обобщенного метода конечных интегральных преобразований (КИП) [16]. Введем на сегменте вырожденное КИП с неизвестными компонентами вектора-функции ядра преобразования: , (2.12) , (2.13) , где - положительные параметры, образующие счетное множество . Круговые частоты осесимметричных колебаний цилиндра связаны с зависимостью . (2.14) Подвергаем систему уравнений (2.6) и начальные условия (2.9) относительно функций преобразованиям КИП в соответствии со структурным алгоритмом [16], представляя кусочно-гладкие функции в виде (2.15) В результате получаем счетное множество задач Коши для трансформанты , решение которых имеет вид , (2.16) и однородную краевую задачу с учетом граничных условий (2.7), (2.8) относительно компонент ядра преобразований : (2.17) (2.18) , , , (2.19) , . В соотношениях (2.16)-(2.19) приняты следующие обозначения: , , . При исследовании системы уравнений (2.17) необходимо отдельно рассмотреть задачи теории электроупругости и упругости . В работах автора [21,22] приведены общие решения данных задач, в которых получены выражения для функций . Их подстановка в граничные условия (2.18), (2.19) формирует однородную систему уравнений относительно постоянных . Разыскивая ее нетривиальные решения, получаем трансцендентное уравнение для вычисления собственных значений , а также выражения постоянных интегрирования. Функции , входящие в разложения (2.10), с учетом условий (2.11) принимаются в виде , , , , , , , . Окончательные выражения функций , получим, применяя к трансформанте (2.16) последовательно формулы обращения КИП (2.13) и метод конечных преобразований Ханкеля (2.5). В результате с учетом (2.2), (2.10) имеем , (2.20) , . Функция определяется с помощью второго равенства (2.20) при удовлетворении условия об отсутствии вертикальных перемещений срединной поверхности пластины при : . (2.21) Полученные расчетные соотношения (2.20) удовлетворяют дифференциальным уравнениям (1.1) и краевым условиям (1.2)-(1.5) (последнее граничное условие (1.2) только для срединной плоскости), т.е. представляют замкнутое решение рассматриваемой краевой задачи электроупругости. 3. Численные результаты. Выводы В качестве примера рассматривается биморфная пластина, имеющая следующие физические и геометрические характеристики аксиально поляризованных пьезокерамических пластин состава ЦТС-19 [23] и металлической стальной подложки: H/м, Кл/м, Ф/м, кг/м, H/м, кг/м, м, м. Анализ полученных расчетных соотношений подтверждает экспериментальные данные, а именно: при жестком закреплении биморфной конструкции изгибные колебания возбуждаются только при наличии системы разрезных кольцевых электродов, расположенных на лицевых поверхностях пьезокерамических пластин. Причем в этом случае электрический потенциал подводится к соседним электродам в противофазе. Рассмотрим работу электроупругой системы на первой резонансной частоте. В этом случае для наиболее эффективного преобразования внешнего электрического воздействия необходимо использовать два электрода (количество и размеры электродов определяют нулевые значения функции ) с радиусом их раздела м Представляем электрическую нагрузку в виде , где амплитуда и частота внешнего воздействия в безразмерной форме. На рис. 2-5 изображены графики, характеризующие изменение по пространственным координатам и времени составляющих механических и электрических полей напряжений На основании анализа результатов расчета можно сделать следующие выводы. 1. При высокочастотном внешнем воздействии наблюдается сложная зависимость изменения напряженно-деформированного состояния системы во времени (рис. 2). Поэтому обычно применяемое при исследовании подобных задач [2-4] допущение об установившемся режиме вынужденных колебаний можно использовать только при решении задач на собственные значения. Рис. 2. Изменение - 1; - 2 по времени (пунктирная линия показывает характер изменения по времени электрической нагрузки) Рис. 3. Изменение амплитудных значений - 1; - 2; - 3 по радиальной координате Рис. 4. Изменение амплитудных значений - 1; - 2; - 3 по аксиальной координате Рис. 5. Изменение амплитудных значений - 1; - 2; - 3 по аксиальной координате 2. Использование двух разрезных круговых электродов и отсутствие угла поворота в жестком закреплении (см. рис. 3, кривая 1) приводят к образованию в пьезокерамических пластинах в радиальной плоскости зон растяжения и сжатия (см. рис. 3, кривая 3). Кроме того, максимальное значение радиальной компоненты вектора перемещений наблюдается в области разреза электродов (см. рис. 3, кривая 2). 3. График «» (см. рис. 4, кривая 1) показывает, что гипотезу плоских сечений, применяемую в прикладной теории для тонких систем, можно использовать также при исследовании составных тонких пластин, имеющих различные физические характеристики по высоте сечения. 4. В зоне раздела электродов (см. рис. 4, кривая 2, ) наблюдаются максимальные значения касательных напряжений . Поэтому для расчета тонких электроупругих систем с разрезными электродами с помощью прикладных теорий необходимо использовать систему уравнений Тимошенко, учитывающую деформации поперечного сдвига. 5. При расчете тонких пластин используется допущение об отсутствии давления слоев друг на друга. Однако в пьезокерамических тонких пластинах нормальные напряжения существенны (см. рис. 4, кривая 3). 6. Потенциал электрического поля по высоте пьезокерамической пластины изменяется по параболической зависимости (см. рис. 5, кривая 1), соответственно, график, описывающий изменение аксиальной компоненты вектора напряженности , представляет собой прямую линию (см. рис. 5, кривая 2). 7. Радиальная компонента вектора напряженности электрического поля практически равна нулю, за исключением точек раздела двух электродов. Максимальное значение данной функции наблюдается на лицевых поверхностях пьезокерамических пластин . Однако даже незначительная величина электрического поля в радиальной плоскости вносит существенный вклад в график изменения аксиальной компоненты индукции электрического поля по высоте пьезокерамической пластины (в прикладных теориях функция принимается постоянной по высоте сечения [5, 6, 8, 12]). В данной задаче разница численных значений по координате (см. рис. 5, кривая 3) достигает 30 %. В заключение следует отметить, что данный подход можно использовать также при расчете многослойных упругих и электроупругих пологих сферических оболочек.

About the authors

D A Shlyakhin

Samara State University of Architecture and Civil Engineering

References

Statistics

Views

Abstract - 135

PDF (Russian) - 54