CC BY
23
Конечно-элементное моделирование в ACELAN пьезоэлектрического элемента устройства накопления "зеленой" энергии
А.Н. Соловьев1, Д.А.Ермаков1, О.Н. Лесняк1, Ле В. З.2
1 Донской государственный технический университет, г.Ростов-на-Дону,
Российская Федерация 2Государственный технический университет имени Ле Куи Дона, г. Ханой, Вьетнам
Аннотация: Статья посвящена моделированию работы пьезоэлектрического генератора (ПЭГ) в виде круглой биморфной пластины с двумя пьезоактивными слоями и инерционной массой. Рассматриваемый ПЭГ может служить элементом устройства накопления энергии, в качестве источника электрической энергии, получаемой из вибраций элементов конструкций и машин.
Рассматриваются изгибные осесимметричные колебания составного упругого и электроупругого тело в линейной постановке, с учетом диссипации механической и электрической энергии. Задача решается с помощью метода конечных элементов (МКЭ), в ЛСБЬЛК. Постановка задачи моделирует закрепление ПЭГ в подвижной конструкции. Проведен модальный анализ устройства, найдены первые частоты резонанса, антирезонанса и коэффициента электромеханической связи. Далее решается задача вынужденных гармонических колебаний на первой частоте антирезонанса и определяется выходной электрический потенциал. Построены зависимости выходных характеристик ПЭГ от его геометрических параметров.
Ключевые слова: пьезоэлектрический генератор, накопление энергии, электрических потенциал, метод конечных элементов.
Требования подвижности энергетической независимости являются главными требованиями при проектировании сенсорных и независимых систем. Помимо этого существует область малой энергетики, где мощности потребляемой энергии малы: бытовая техника, мобильные телефоны, беспроводные средства диагностики технического состояния различных объектов. В современном мире имеются разнообразные решения этой проблемы. Преимущество источников питания телефонов и ноутбуков уделяется литий-ионным батареям. Однако, как показывает практика, их ресурс ограничен.
Последнее время появились попытки уменьшить энергопотребление датчиков и преобразователей и спроектировать батареи со сверхдлительным
сроком эксплуатации и передачу различных энергий (например, сверхчастотной, световой и магнитной) по воздушным каналам.
Упомянутые факторы побудили ученых и инженеров искать эффективные самовозобновляемые автономные источники энергии. Одним из видов таких устройств являются «пьезоэлектрические устройства накопления энергии» - устройства, собирающие энергию с пьезоэлементов, преобразующих рассеянную механическую энергию в электрическую с последующим накоплением и передачей приемному устройству. Анализируя требования, которые предъявляются к накопителям и преобразователям энергии окружающей среды, можно сказать, что они существенно зависят от того, для каких типов устройств они предназначены.
Пьезоэлектрические материалы дают возможность прямого преобразования механической энергии деформации в электрический заряд и наоборот. Пьезоэлектрические генераторы (ПЭГ) являются основным элементом устройств накопления энергии и выполняют функцию получения энергии из среды, и преобразования ее в электрическую. Концепция устройств накопления энергии стремится к тому, чтобы создать элементы, которые не нуждались в замене.
Пьезоэлектрические генераторы позволяют собирать энергию от маломощных источников (вибрация от качающихся мостов, движущихся механизмов и т.д). В настоящее время пьезоэлектрические устройства накопления энергии активно изучаются.
Возможность выработки энергии изучается из использования повседневных движений человека, таких как кровяное давление, тела, сжатие и расширение грудной клетки при дыхании, и движение ног при ходьбе [13].
В работах [4,5] рассматривается построение конечно-элементых (КЭ) моделей и решение задач оптимизации конструкции устройств накопления
энергии на основе пьезоэлектрических элементов для различных типов пьезоэлектрических генераторов, в работах [6,7] описывается построение моделей ПЭГ в пакетах АКБУБ, АСЕЬАК и БкхРБЕ и результаты расчетов, с помощью которых проведен анализ применимости прикладных теорий расчета ПЭГ.
Авторы работы [8] предложили устройство, которое может преобразовать механическую энергию в электрическую, если его встроить в вибрирующую среду и использовать для питания других устройств.
Проанализировав эффективность пьезоэлектрического вещества в устройстве стековой конфигурации на основе аналитической модели в работе был сделан вывод о том, что, что максимум выработки можно достигнуть за счет минимизации количества энергии [9].
Одним из наиболее результативных направлений проектирования пьезоэлектрических генераторов считается применение механических вибраций для появления энергии деформации в пьезокерамике. Источниками энергии могут быть колебания окружающей среде после их преобразования в электрическую энергию[10].
В данной работе рассматривается пьезоэлектрический генератор, представляющий собой круглую пластину, которая имеет конструкцию биморфа. На дисковидный центральный слой приклеены два симметрично расположенных пьезоактивных слоя меньшего диаметра из пьезокерамики. Также исследуются различные характеристики при гармонически изменяющейся во времени вибрации и рассмотрены пути оптимизации конструкции. Достижение наибольшего выходного потенциала производится путем изменения соотношения различных геометрических параметров
Материалы и методы. Пьезоэлектрический генератор (рис. 1) устройства накопления энергии это составная конструкция, состоящая из упругих и электроупругих элементов, совершающих малые колебания в
подвижном системе координат, которая совершает прямолинейное поступательное движение вдоль оси симметрии (кинематическое возбуждение)
Рис. 1. - Параметры модели В этих условиях математической моделью функционирования устройства является начально-краевая задача линейной теории электроупругости [7].
Для пьезоэлектрической среды: püt + арй t - ст.,= f; Di. = ^
а. = суы (£и + ) - etjA; Dt + QDt = etk¡ (ski+SdéM)+э,Л
£k¡ = (ükj + üi,k )/2; Ek = (, где p - плотность материала; üi - компоненты вектора - функция перемещений; а. - компоненты тензора механических напряжений; f -компоненты вектора плотности массовых сил; Di - компоненты вектора электрической индукции; cjjkl - компоненты тензора упругих модулей; Ei -компоненты вектора напряженности электрического поля; ejjk - компоненты
(1) (2) (3)
тензора пьезомодулей; sti - компоненты тензора деформаций; р -
электрическим потенциал; э } - компоненты тензора диэлектрических
проницаемостей; а,Р,да - неотрицательные коэффициенты демпфирования. Для упругой среды:
рйх + арйг - оц,} = £; (4)
^ = } ); (5)
о=(йк ,1+й л)/2; (6)
Механические условия задаются на объединении непересекающихся областей 5 = 8й и St и :
йА3и = й (7)
,, = } 15, = Р, (8)
йи = щщ^ = О, Т = о (9)
Электрические краевые условия задаются на совокупности областей
м
5 = и ^ и ^ :
m=1
-на электродах задается значение электрического потенциала
P\se = Й = const (10)
- на неэлектродированной поверхности
D„\Sd = Dini\sD = De (= 0) (11)
- если электрод подключен к внешней электрической цепи, то значение электрического потенциала на нем в (10) неизвестно и для его нахождения используется дополнительное условие
J DndS = Im (= 0) (12)
e
^m
где Im - ток в цепи.
В программном комплексе ACELAN для описания потерь механической и электрической энергии используются коэффициенты
затухания а,Р,да, которые выражаются через добротность б следующим образом:
2¥л!г 2 в = с= 1
б/ + / 2) 2Л()(/Л + /г2) (13)
где /г1, /г2 - частоты резонанса, которые находятся при использовании граничного условия (10) или частоты антирезонанса /а1, /а2 для определения
которых используются (10), (12).
В качестве инструмента решения краевой задачи (1)-(12) используется конечно-элементный комплекс ЛСБЬЛК [7-11]. Рассматриваются установившиеся колебания, находятся частоты резонанса и антирезонанса, определяется выходной электрический потенциал на частоте антирезонанса.
На границах областей (рис. 1) задаются следующие граничные условия: линия 1 - ось симметрии (9); на линию 2 действует нормальное давление (8); линия 3 - гладкий контакт, условие (9); линии 4 и 5 - электроды, 4 - условие (10) потенциал равен 0; 5 -условия (10), (12); точка 6 - шарнирное закрепление.
Результаты исследования. Конструкция ПЭГ имеет ряд геометрических параметров, на рис 1 представлена половина осевого сечения конструкции с обозначением геометрических параметров. В работе изучается влияние этих параметров на выходные характеристики. В частности будет рассмотрено влияние на выходной потенциал р следующих параметров: радиуса пьезоэлемента гр, высоты инерционной массы кт, толщины
пьезоэлемента Ър и радиуса стойки г.
Расчеты приводятся при неизменных в процессе исследования следующих значениях параметров: радиуса инерционной массы гт = 20 мм, высоты стойки гк = 2.5 мм, высоты пластины ^ = 0.3 мм и радиуса пластины г■ = 20 мм. При этом приняты ограничения на геометрические
размеры конструкции согласно рисунку 3: кр = 0,2 + 1 мм; гр = 13 ^19 мм; гг = 2 ^ 6 мм; кт = 5 + 10 мм.
Первым шагом численного эксперимента является определение частот резонанса, антирезонанса и коэффициента механической связи К (КЭМС).
К =
1 -
й
/2
с/ а
(14)
Рис.2. - График зависимости частот резонанса, антирезонанса и КЭМС от
радиуса пьезоэлемента На рис. 2 представлены зависимости частот резонанса, антирезонанса и КЭМС от радиуса пьезоэлемента, из которого видно, что частоты увеличиваются с увеличением радиуса, а КЭМС достигает максимума при радиусе пьезоэлемента 18,3 мм. Расчеты проводятся при следующих значениях параметров: кр = 0,6мм, гг = 2,5 мм, кт = 7,5 мм.
Таблица 1
Зависимость частот резонанса, антирезонанса и КЭМС от толщины
пьезоэлемента
^р (мм) й ^ а Гг К
0,2 271 295 0,395
0,3 356 390 0,4
0,4 441 480 0,41
0,5 518 558 0,42
0,6 582 619 0,43
0,7 638 671 0,44
0,8 681 710 0,45
0,9 712 740 0,46
1 730 761 0,47
В табл.1 представлены зависимости частот резонанса и антирезонанса, а также КЭМС от толщины пьезокерамического слоя, из которого видно, что частоты увеличиваются с увеличением толщины, а максимум достигается при высоте пьезокерамического слоя 1 мм. Расчеты проводятся при следующих значениях параметров: гр = 15 мм., гг = 2,5 мм., кт = 7,5 мм.
Таблица 2
Зависимость частот резонанса, антирезонанса и КЭМС от высоты инерционной массы
к / \ т (мм) /" а Гг К
5,5 672 715 0,33
6 646 684 0,36
6,5 623 663 0,37
7 602 640 0,39
7,5 582 619,5 0,4
8 565 601 0,41
8,5 549 584 0,42
9 535 570 0,43
9,5 522 555 0,44
10 509 545 0,44
В табл. 2 представлены зависимости частот резонанса и антирезонанса, а также КЭМС от высоты инерционной массы, из которых видно, что
частоты уменьшаются с увеличением высоты инерционной массы, а максимум КЭМС достигается при высоте инерционной массы 10 мм. Расчеты проводятся при следующих значениях параметров: гр = 15 мм.,
гг = 2,5 мм., кр = 0.6 мм.
Предполагается, что подвижная система координат, в которой закреплен ПЭГ движется по закону:
где г0 = 0.01 мм - амплитуда колебаний, с - круговая частота колебаний.
При этом инерционные силы, действующие на элементы ПЭГ, учитываются с помощью давления приложенного (рис. 1 линия 2) которое зависит от параметров движения системы координат (15) и массы элементов
На инерционную массу прикладывается равномерно распределенное давление ап
где тин - инерционная масса, тр - масса пьезоэлементов, тст - масса стойки, тпл - масса пластины, с = 2 -п- /а. При этом на рис. 3-5 и табл. 3 представлены зависимости выходного потенциала при изменении геометрических параметров модели. При этом колебания возбуждаются на частоте антирезонанса.
г (/) = 20 е
(16)
ПЭГ.
^п = (тии + тг + тст + тгл)20 ■ с
(17)
■£■ Инженерный вестник Дона. №1 (2019) ВЦ ivdon.ru/ru/magazine/arcliive/nly2019/5587
(Р 250 200 150 100 50 0
13 13,5 14 14,5 15 15,5 16 16,5 17,2 18 у (лш)
Рис. 3. - Зависимость выходного потенциала от радиуса пьезоэлемента На рис. 3 представлена зависимость выходного потенциала от радиуса пьезоэлемента, из которой видно, что потенциал увеличивается с увеличением радиуса. Расчеты проводятся при следующих значениях
параметров: кр = 0,6 мм., гг = 2,5 мм,. кт = 7,5 мм. (р ™>
600 500 400 300 200 100
0 Н-1-1-1-1-1-1-1-1-1
2 2,5 3 3,5 4 4,5 5 5,5 /• (ад ад)
Рис. 4. - Зависимость выходного потенциала от радиуса стойки На рис. 4 представлена зависимость выходного потенциала от радиуса стойки, из которой видно, что потенциал увеличивается с увеличением радиуса. Расчеты проводятся при следующих значениях параметров кр = 0,6 мм, гр = 15 мм, кт = 7,5 мм.
(Р 450 400 350 300 250 200 150 100 50
2 2,5 3 3,5 4 4,5 5 5,5 б 6,5 Ц ^ у
Рис. 5. - Зависимость выходного потенциала от высоты инерционной массы
На рис. 5 представлена зависимость выходного потенциала от высоты инерционной массы, из которой видно, что потенциал увеличивается с увеличением высоты. Расчеты проводятся при следующих значениях параметров: Расчеты проводятся при следующих значениях параметров: гр = 15 мм, гг = 2,5 мм, кр = 0.6 мм.
Таблица 3
Зависимость выходного потенциала от толщины пьезоэлемента
кр (мм) 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
40 87 120 171.5 208 238 257.6 268 270
В табл.3 представлена зависимость выходного потенциала от высоты пьезоэлемента, из которой видно, что потенциал увеличивается с увеличением высоты. Расчеты проводятся при следующих значениях параметров: гр = 15 мм, гг = 2,5 мм, кт = 7,5 мм.
Обсуждение и заключение. В данной работе исследована задача эффективности ПЭГ в зависимости от геометрических параметров и
характера возбуждения. Построена конечно-элементная модель в программном комплекса ACELAN, позволяющая менять геометрические параметры устройства, вычислять выходной электрический потенциал.
При кинематическом возбуждении выходной потенциал и частоты антирезонанса увеличиваются во всех рассмотренных случаях, таких как увеличение высоты и радиуса пьезоэлемента, высоты инерционной массы и стойки. Коэффициент электромеханической связи достигает максимума при значении радиуса пьезоэлемента rp = 18,3 мм, высоте инерционной массы
hm = 7,5 мм и высоте пьезокерамического слоя hp = 1мм.
Литература
1. Hasler E., Stein L., Harbauer G. Implantable physiological power supply with PVDF film // Ferroelectrics. - 1984. - V. 60. - № 1. - pp. 277-282.
2. Kymissis J., Kendall C., Paradiso J., Gershenfeld N. Parasitic power harvesting in shoes. // IEEE, 1998. pp. 132-139.
3. Starner T. Human-powered wearable computing // IBM systems Journal. - 1996. - V. 35. - № 3.4. - pp. 618-629.
4. Соловьев А.Н., Зыонг Ле В. Конечно-элементное моделирование пьезоэлектрического устройства накопления энергии на основе кантилевера // Вестник Донского государственного технического университета. - 2014. -№ 1. - C. 169-179.
5. Зыонг Л.В., Соловьев А.Н. Конечно-элементное моделирование и анализ пьезоэлектрического устройства накопления энергии в форме круглой пластины с пьезоэлементами // Экологический вестник научных центров Черноморского экономического сотрудничества. - 2013. - Т. 1. - № 4. - C. 112119.
6. Акопьян В.А., Паринов И.А., Рожков Е. В., Чебаненко В.А., Зыонг Л.В. Экспериментальное и конечно-элементное моделирование устройств накопления энергии // "Математическое моделирование и
биомеханика в современном университете». Сборник трудов IX всероссийской школы-семинара (г. Ростов-на-Дону, 26 мая - 30 мая 2014 года). Изд-во ЮФУ. - 2014. - C. 57.
7. Зыонг Л.В., Рожков Е. В. Оптимизация конструкции устройств накопления энергии на основе пьезоэлектрических элементов с помощью конечно-элементного моделирования // "Математическое моделирование и биомеханика в современном университете». Сборник трудов VIII всероссийской школы-семинара (г. Ростов-на-Дону, 27 мая - 31 мая 2013 года). Изд-во ЮФУ. - 2013. - C. 61.
8. Williams C.B., Yates R.B. Analysis of a micro-electric generator for microsystems // Sensors and Actuators A: Physical. - 1996. - V. 52. - № 1. - pp. 811.
9. Соловьев А.Н., Оганесян П. А. Конечноэлементное моделирование конструкций из неоднородных материалов с усложненными свойствами // Вестник ДГТУ. 2013. № 7-8, с. 21-27.
10. Skaliukh, A.S., Soloviev, A.N.,Oganesyan, P.A. Modeling of piezoelectric elements with inhomogeneous polarization in ACELAN // Ferroelectrics 2015 483 (1), pp. 95-101
11. Мартыненко А.В. Исследование параметров осесимметричного изгибного пьезоэлектрического преобразователя при различных граничных условиях // Инженерный вестник Дона, 2010, №3. URL: ivdon.ru/magazine/archive/n3y2010/206
12. Зыонг Ле В. Конечноэлементный анализ применимости прикладных теорий расчета пьезоэлектрического устройства накопления энергии стековой конфигурации // Инженерный вестник Дона, 2014, №2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2310/.
References
1. Hasler E. Ferroelectrics. International journal. 1984. V. 60. № 1. pp. 277282.
2. Kymissis J., Kendall C., Paradiso J., Gershenfeld N. IEEE Spectrum. International journal. 1998. pp. 132-139.
3. Starner T. IBM systems Journal. International journal. 1996. V. 35. № 3.4. pp. 618-629.
4. Solovev A.N., Zyong L. V. Vestnik Donskogo gosudarstvennogo texnicheskogo universiteta. Nauchny'j zhurnal 2014. № 1. pp. 169-179.
5. Soloviev A.N., Zyong L.V. E'kologicheskij vestnik nauchnyx centrov Chernomorskogo e'konomicheskogo sotrudnichestva. Nauchno-obrazovatelnyj i prikladnoj zhurnal. 2013. № 4. pp. 112-119.
6. Akopyan V.A., Parinov I.A., Rozhkov E. V., Chebanenko V.A., Zyong L.V. Vserossiyskaya shkola-seminar "Matematicheskoe modelirovanie i biomexanika v sovremennom universitete": trudy (Proc. All-Russian Workshop School. "Mathematical modeling and biomechanics in modern university"). Rostov-on-Don, 2014, pp. 57.
7. Zyong L.V. Rozhkov E. V. Vserossiyskaya shkola-seminar "Matematicheskoe modelirovanie i biomexanika v sovremennom universitete": trudy (Proc. All-Russian Workshop School. "Mathematical modeling and biomechanics in modern university"). Rostov-on-Don, 2013, pp. 61.
8. Williams C.B., Yates R.B. Sensors and Actuators A: Physical. International journal. 1996. № 1. pp. 8-11.
9. Soloviev A.N., Oganesyan P.A. Vestnik Donskogo gosudarstvennogo texnicheskogo universiteta. Nauchny'j zhurnal. 2013. № 7-8, pp. 21-27.
10. Skaliukh, A.S., Soloviev, A.N.,Oganesyan, P.A. Ferroelectrics. International journal. 2015 483 (1), pp. 95-101.
11. Martynenko A.V. Inzenernyj vestnik Dona (Rus), 2010, №3. URL: ivdon.ru/magazine/archive/n3y2010/206/.
12. Ье V. 1п2епегпу| уеБ1шк Бопа (Киб), 2014, №2. иКЬ:
ivdon.ru/ru/magazine/archive/n2y2014/2310/.
