ТРАНСФОРМАЦИЯ КОЛЕБАНИЙ НЕУСТОЙЧИВОЙ СИСТЕМЫ В ПРЕОБРАЗОВАТЕЛЕ-НАКОПИТЕЛЕ ЭНЕРГИИ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ MATHEMATICAL MODELLING

Научная статья DOI: 10.18287/2541-7525-2023-29-2-7-18

УДК 517; 519.6; 519.8 Дата: поступления статьи: 01.03.2023

после рецензирования: 11.04.2023 принятия статьи: 30.06.2023

С.В. Борзунов

Воронежский государственный университет, г. Воронеж, Российская Федерация E-mail: sborzunov@gmail.com. ORCID: https://orcid.org/0000-0002-5099-9655

ТРАНСФОРМАЦИЯ КОЛЕБАНИЙ НЕУСТОЙЧИВОЙ СИСТЕМЫ В ПРЕОБРАЗОВАТЕЛЕ-НАКОПИТЕЛЕ ЭНЕРГИИ

АННОТАЦИЯ

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

Ключевые слова: преобразователь-запасатель энергии; перевернутый маятник; устойчивость системы

Цитирование. Борзунов С.В. Трансформация колебаний неустойчивой системы в преобразователе-накопителе энергии // Вестник Самарского университета. Естественнонаучная серия. 2023. Т. 29, № 2. С. 7-18. DOI: http://doi.org/10.18287/2541-7525-2023-29-2-7-18.

Информация о конфликте интересов: авторы и рецензенты заявляют об отсутствии конфликта интересов.

© Борзунов С.В., 2023

Сергей Викторович Борзунов — кандидат физико-математических наук, доцент кафедры цифровых технологий, Воронежский государственный университет, 394018, Российская Федерация, г. Воронеж, Университетская пл., 1.

Введение

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

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

Различные механизмы поддержания резонансных режимов посредством вибрационных колебательных систем с настраиваемыми характеристиками рассматривались в работах [1—6]. Однако рассматриваемые в этих работах системы требует прецизионной настройки параметров, не всегда реализуемой в реальных технических системах. В частности, большое внимание уделяется величинам электрического напряжения [7], температурным условиям [8], шумовым воздействиям [9].

В частности, ряд работ посвящен электрическим цепям, включающим в свой состав пьезоэлектрические материалы, т. е. вещества, в которых при упругих деформациях возникает электрическая поляризация. Они представляют собой кристаллические вещества без центра симметрии и характеризуются сложной, нелинейной зависимостью между приложенным механическим напряжением и создаваемым электрическим полем. Среди пьезоэлектриков выделяют класс сегнетоэлектрических материалов, которые обладают в определенном диапазоне температур ненулевой поляризацией, изменяющейся за счет внешних воздействий [10; 11]. Именно сегнетоэлектрики служат материалом для построения микроэлектромеханических систем (microelectromechanical systems, MEMS): излучателей звука, акселерометров, прецизионных датчиков микроперемещений [8; 12; 13].

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

Экспериментальные данные показывают, что известные модели с достаточной точностью описывают пьезокерамический преобразователь, когда он подвергается воздействию низкочастотного электрического напряжения. Точность моделей ухудшается, если диапазон частот напряжения становится шире, а также в случае относительно больших механических нагрузок [14].

1. Модель преобразователя-запасателя энергии в линейном приближении

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

Модель изучаемой системы состоит из перевернутого математического маятника, закрепленного на легкой горизонтальной платформе и соединенного механической связью с сегнетоэлектрическим конденсатором, который включен в замкнутую электрическую цепь. Платформа P может перемещаться в горизонтальном направлении. Обозначим угол отклонения маятника относительно вертикали через <^(t), координату платформы — через u(t) (рис. 1.1). Длина маятника равна l, масса его груза равна m. Связь между механической и электрической подсистемами содержит звено, подчиняющееся закону вязкого трения с коэффициентом диссипации с [соответствующая сила трения равна Ffr = cv(t) = сjt(1ф/2 + u)].

Пьезоэлектрический материал, образующий конденсатор с емкостью C, включен в электрическую цепь с внешней омической нагрузкой R. Напряжение на нагрузке обозначим через V(t).

Рассматриваемая динамическая система описывается системой обыкновенных дифференциальных уравнений

{т1ф + с(1ф/2 + U) — mg sin ф + AV = —mU,

(1.1)

CV+ R = B (1ф/2 + U).

Первое из уравнений системы (1.1) представляет собой уравнение движения груза m в условиях действия

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

Точкой над символом обозначена производная по времени t.

В приближении небольших отклонений маятника от положения равновесия sin ф ~ ф получаем (далее

2

0

введены обозначения yo = c/(2m), w2 = g/l):

f ф + 7оФ - ^0ф + mv = -1it- 2f u,

I ^ + 11av - %ф = Cu '

В этой системе A и B — параметры связи (coupling) механической и электрической подсистем, которые выводятся на основе следующих рассуждений. Уравнение относительно неизвестной функции ф(Ь] представляет собой уравнение движение осциллятора, находящегося под воздействием, помимо механических сил, внешней силы немеханического происхождения [-Fext(t) = AV(t)], а также силы инерции Finer(t) = -mu(t).

Рис. 1.1. Схематическое описание математического маятника, связанного с пьезоэлектрическим

генератором

Fig. 1.1. Scheme of a mathematical pendulum associated with a piezoelectric generator

Закон упругих деформаций (закон Гука) связывает натяжение а и относительную деформацию (для однородных и изотропных материалов в случае умеренных абсолютных величин упругих деформаций):

а = Ее, (1.3)

где Е — модуль Юнга, зависящий только от материала пьезоэлектрика и его физического состояния [в нелинейном случае а = Ее + Е2е2 + 0(е2)].

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

Px = dn(rx - Ту), (1.4)

где тх и ту — механические натяжения, действующие параллельно осям Ox и Oy, a du — постоянная, называемая пьезоэлектрическим модулем. Выразим заряд Q = CV(t), сформированный на гранях образца (примем за S площадь одной грани, C — электрическая емкость конденсатора, образованного пьезоэлектрической пластинкой):

Q = Px S, (1.5)

Q = dnF. (1.6)

Следовательно, между упругой силой, действующей со стороны пьезоэлектрика на груз m, и смещением S(t) = 1ф/2 + u существует связь в виде линейной пропорциональности:

F = ^ V = AV, (1.7)

C

где A — константа, зависящая только от материала и его диэлектрика термодинамических свойств. Другими словами, напряжение линейно зависит от смещения по формуле V(S) = ELS(t), где L — длина пьезоэлектрической пластинки [т. к. закон Гука можно представить в виде F/S = E х (S/L)].

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

В электрической цепи, состоящей из внешнего источника тока (за счет поляризации пластин диэлектрика) и параллельно соединенных конденсатора C и резистора R (выполняющего роль активной нагрузки), действует закон Ома:

CV+ V = BS, (1.8)

R

где B = CL = dEL =const.

Заметим, что в литературе, посвященной моделям накопителей энергии, описаны электрические подсистемы, построенные на тех же принципах, но удовлетворяющие другим дифференциальным уравнениям, в частности, в работе [15]:

в работах [16; 17]:

V V

BV + CV + v + v =0, R L

RQ - C x + Q =0

(1.9)

(1.10)

(здесь d(x) — зависящий от механического параметра х коэффициент связи между механической и электрической подсистемами).

Среди многих вариантов механической колебательной системы на основе маятника, предлагаемых для использования в накопителях энергии, отметим маятник на перевернутой балке [18], перевернутый маятник с ограничителями амплитуды (inverted pendulum amplitude limiters) [19], перевернутая консольная балка с массой на конце (cantilever beam with a tip mass) [20].

1.1. Приведение системы к безразмерному виду

Перейдем от переменных t и V к безразмерным переменным по следующему правилу:

t = tc • т, V = Vc • v, (1.11)

где характерные величины времени tc и напряжения Vc будут определены ниже. Угол отклонения маятника ф(Ь) является безразмерной величиной. Система (1.1) принимает следующий вид:

1 Y0 ф - ф + m Vcv = - ^ U(tcT) - 2 Y0 utr),

£v + ReVcv - ОСЬф = Ct¿fa).

(1.12)

Здесь и далее точкой над символом обозначена производная по безразмерному времени т (двумя точками — вторая производная по т). Выполним алгебраические преобразования:

ф + YotcФ - -О^ф + AVctOv = -1 V(tcT) - 2^U(tcT),

v + Retcv - Oevkфр = cku(tcT)-

(1.13)

Пусть функция, определяющая движение платформы Р, равна п{Ь) = и(Ьст), а значения размерных коэффициентов: Ьс = ЕС, Ус = А(пС)2. Тогда получим основную систему в безразмерном виде:

V + Yo(RC)Ф - -O(RC)2ф + v = -1 U(t) - 2

,7o (RC) ~

(т),

(1.14)

* + - - АВЕС= АВЕ2Си(т).

С целью упрощения и ясности дальнейших вычислений введем обозначения для безразмерных величин: коэффициента затухания, характеристической частоты маятника и внешнего воздействия: 7 = 70 • (ЕС), ш = шо • (ЕС), и> = и(т)/1. Также удобно обозначить коэффициент связи электрической и механической подсистем через а = Е2С. Окончательно

¡V + Yфр - —2ф + v = -w(r) - 2yw(t), v + v - афр = 2aw(r),

или в более удобной матричной форме:

d dr

V 0 1 0 V

ф = -О -Y -1 ф +

v 0 а -1 v

0

-w - 2yw 2aw

(1.15)

(1.16)

где ф(т) = ф(т) — угловая скорость маятника. Заметим, что согласно физической постановке задачи выполняются неравенства ш > 0, 7 ^ 0.

1.2. Исследование устойчивости линеаризованной системы

Исследуем вопрос об устойчивости решений системы обыкновенных дифференциальных уравнений (1.16). Для этого выпишем характеристическое уравнение системы (1.16) в автономном случае:

-Л

ш2 0

1 0 -7 - Л -1 а —1— Л

0,

(1.17)

оно приводится к алгебраическому уравнению третьей степени с постоянными коэффициентами:

А3 + (7 + 1)А2 + (7 + а - ш2)А- ш2 =0. (1.18)

Как известно, необходимым и достаточным условием устойчивости системы дифференциальных уравнений является положительность всех трех главных диагональных миноров Д^, Д2, Д3 определителя Гурвица

ai a3 0

Г = ao a2 0 , (1.19)

0 ai аз

где щ, г € {0,1, 2, 3} — коэффициенты алгебраического уравнения (1.18), причем ао = 1 (см., например, [21]). Поскольку а3 = —ш2 < 0, то условия критерия Гурвица Д1 > 0, Д2 > 0, Д3 > 0:

а1 > 0,

а1а2 — а0а3 > 0, (1.20)

аз(а1 а2 — аоаз) > 0,

как легко видеть, несовместны ни при каких допустимых значениях параметров ш, 7, а. Следовательно, система (1.16) является неустойчивой. Тем не менее введение управления с помощью внешнего воздействия и>(т) позволяет в ряде случаев перевести динамическую систему в устойчивый режим колебаний. В следующем разделе рассмотрим, при каких значениях параметров задачи это возможно.

1.3. Управление амплитудой колебаний маятника

В настоящем разделе покажем, что введение управления с помощью внешнего воздействия и>(т) позволяет стабилизировать движение маятника и предоставляет возможность управления амплитудой ф(т) его колебаний.

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

Запишем систему уравнений (1.16) в условиях управления вида «;(£) = аф(Ь) + ¿вгПг, где й € С:

или

d dr

d dr

ф 0 1 0 ф

Ф = ш2 -7 -1 ф +

v 0 а -1 - v

ф г 0 1 0

ф = и2 - 27а -(7 + а) -1

v 2аа а -1

0

-аф - 27аф - вШт(itt + 27)d 2ааф + 2aiQeT d

ф

ф

v

+

0

-(in + 27)eiQT

2aiQeinT

d.

(1.21)

(1.22)

Как видно из (1.22), коэффициенты матрицы системы в результате действия управления изменились. Проверим условия критерия Гурвица в этом случае. Выпишем характеристическое уравнение

Л3 + (1 + 7 + а)Л2 + [7 + а - и2 + (1 + 27)а]Л + [-ш2 + 2(7 + а)а] = 0 и матрицу Гурвица

,2

Г(а) =

1+ 7 + а -ш2 + 2(7 + а)а 1 0

0

7 + а - и2 + (1 + 27) а 0

1+ 7 + а -ш2 + 2(7 + а)а

(1.23)

(1.24)

которая приводит к следующим условиям:

1 + 7 + a> 0,

(1 + 27)a2 + (1 + 27 + 272 - a - ш2)а + (1 + 7X7 + a) - 7J2 > 0, -J1 + 2(7 + a)a > 0.

(1.25)

Итак, система (1.22) обладает устойчивостью тогда и только тогда, когда ее параметры удовлетворяют условиям (1.25). Ее решение относительно параметра управления а не представляет сложностей:

а >

2(т+('

если V > 0,

а > max

2(7+^)'

,аъ) , если D ^ 0,

(1.26)

где аь — больший из двух корней квадратного уравнения (1.25), V — его дискриминант:

V = (1 + 27 + 272 - а - ш2)2 - 4(1 + 27)(7 + 72 + а + 7а - 7ш2). (1.27)

Запишем решение системы (1.22), представленное в виде суммы общего решения однородной системы и частного решения неоднородной:

Ф

Ф

Re (C\eXlт/ + С2вХ2Т/ + С3вХзТ f3 +

Аф Аф

Av

.

(1.28)

В формуле (1.28) Л1, Л2, Лэ — собственные значения матрицы системы (1.22), а /1, /2, /э отвечающие им собственные векторы.

В установившемся режиме, т. е. при достаточно больших значениях времени т ^ » тах(1/|Ие Л1|, 1/|Ке Л2|, 1/|Ке Лэ|)

' Ф ' Аф

Ф — Re Аф

v Av

Шт

(1.29)

Как видно из (1.29), решение в установившемся режиме содержит слагаемые, изменяющиеся по гармоническому закону ~ вгПт. Их амплитуды Аф, Аф, Ау, как несложно показать непосредственной подстановкой {^(т) = АфвгПт, ф(т) = АфвгПт, у(т) = АувгПт} в (1.22), удовлетворяют системе алгебраических уравнений:

-in 1 0 ■ Ф 0

J2 - 27а -(7 + а + in) -1 Ф + -(in+ 27) d — 0. (1.30)

2аа a -1 - in v 2ain

Таким образом, для амплитуды колебаний маятника в зависимости от амплитуды гармонического члена с частотой П в управлении и>(т) получаем формулу:

Аф

1

[27 + i(1 + 27 + 2a)n - n2] d,

(1.31)

ф Т(П)

где Т(П) = ш2 -2а(7 + а) -[(1 + 27)а + 7 + а - ш2] + П2(а + 7 +1) + ,0э — определитель системы (1.30) (заметим, что этот определитель также может быть получен из правой части (1.23) путем выполнения формальной замены Л ^ 1,0,). Установившиеся колебания возможны при условии Т(О) =0.

Аналогичным образом вычисляются амплитуды колебаний угловой скорости маятника и напряжения:

Аф —

in

FW)

[27 + i(1 + 27 + 2а)П - n2] d,

Av —

F (П)

[4а7 + 2i(a + 7 - 2aj + ш2)П + (2а + 27 - 1)П2 + 2in3] d.

(1.32)

(1.33)

Результаты численных расчетов области устойчивости, корней характеристического уравнения и решения системы при различных значениях параметров представлены на рис. 1.2—1.3.

Графики зависимости амплитуды вынужденных колебаний маятника А'А в пространстве параметров {(с1, О)} представлены на рис. 1.4—1.5. На этих рисунках видны убывание максимальной амплитуды с увеличением 7, а также смещение резонансной области в сторону больших частот О в этом случае. В силу линейности системы (1.30) величина АА пропорциональна параметру d в управляющей функции и>(£).

2

2

v

e

a

3.0 2.5 2.0 ш 1.5 1.0 0.5 0.0

3.0 2.5 2.0 ш 1.5 1.0 0.5

0.0

3.0 2.5 2.0 ш 1.5 1.0 0.5 0.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

Ф(т)

ФМ

10

20

30

40

50

2.0 1.5 1.0 0.5 0.0 -0.5

2.5 2.0 1.5 1.0 0. 5 0.0 -0.5

Рис. 1.2. Результаты численных расчетов области устойчивости (1.22), корней характеристического уравнения и решения (1.22) при y = 0,2, а = 0,4. Параметры и начальные условия для угла отклонения ш = 0,4, d = 1,0, Q = 5,0, ф(0) = 1.0, ф(0) = 0.0, -у(0) = 1.0: a — область устойчивости для а = 0,05; b — область устойчивости для а = 0,5; с — область устойчивости для а = 1,5; d — вещественные (толстые линии) и мнимые части (тонкие линии) корней характеристического уравнения для а = 0,05; e — вещественные (толстые линии) и мнимые части (тонкие линии) корней

характеристического уравнения для а = 0,5; f — вещественные (толстые линии) и мнимые части (тонкие линии) корней характеристического уравнения для а = 1,5; g — угол отклонения маятника для а = 0,05; h — угол отклонения маятника для а = 0,5; k — угол отклонения маятника для а = 1,5 Fig. 2.2. Results of numerical calculations of the stability region (1.22), roots of the characteristic equation and the solution (1.22) for y = 0,2, a = 0,4. Parameters and initial conditions for the deflection angle are ш = 0,4, d = 1,0, О = 5,0, ф(0) = 1,0, ф(0) = 0,0, v(0) = 1,0: a — region of stability for a = 0,05; b — region of stability for a = 0,5; с — region of stability for 7 = 0,2; d — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 0,05; e — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 0,5; f — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 1,5; g — pendulum deflection angle for a = 0,05; h — pendulum deflection angle for a = 0,5; k — pendulum deflection angle for a = 1,5

b

a

с

1,2,3

1,2,3

1,2,3

3

d

e

0

h

k

g

3.0 2.5 2.0 ш 1.5 1.0 0.5 0.0

3.0 2.5 2.0 ш 1.5 1.0 0.5 0.0

3.0 2.5 2.0 ш 1.5 1.0 0.5 0.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.0 0.5 1.0 1.5 2.0 2.5 3.0

b

a

с

0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50

г г г

g h k

Рис. 1.3. Результаты численных расчетов области устойчивости (1.22), корней характеристического уравнения и решения (1.22) при 7 = 1,1, а = 1,5. Параметры и начальные условия для угла отклонения ш = 1,2, d = 1,0, l = 5,0, ф(0) = 1,0, ф(0) = 0,0, -у(0) = 1,0: a — область устойчивости системы для а = 0,05; b — область устойчивости системы для а = 0,5; с — область устойчивости системы для а = 1,5; d — вещественные (толстые линии) и мнимые части (тонкие линии) корней характеристического уравнения для а = 0,05; e — вещественные (толстые линии) и мнимые части (тонкие линии) корней характеристического уравнения для а = 0,5; f — вещественные (толстые линии) и мнимые части (тонкие линии) корней характеристического уравнения для а = 1,5; g — угол отклонения маятника для а = 0,05; h — угол отклонения маятника для а = 0,5; k — угол отклонения

маятника для а =1,5

Fig. 2.3. Results of numerical calculations of the stability region (1.22), roots of the characteristic equation and the solution (1.22) for 7 =1,1, a = 1,5. Parameters and initial conditions for the deflection angle are ш = 1,2, d = 1,0, О = 5,0, ф(0) = 1,0, ф(0) = 0,0, v(0) = 1,0: a — region of stability for a = 0,05; b — region of stability for a = 0,5; с — region of stability for a = 1,5; d — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 0,05; e — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 0,5; f — real (thick lines) and imaginary parts (thin lines) of the roots of the characteristic equation for a = 1,5; g — pendulum deflection angle for a = 0,05; h — pendulum deflection angle for a = 0,5; k — pendulum deflection angle for a = 1,5

Рис. 1.4. Амплитуда вынужденных колебаний маятника \Av(d, согласно (1.31) для значений параметра 7, характеризующего трение: a — 7 = 0.2; b — 7 = 1.1; c — 7 = 1.9. Расчет выполнен

для а = 0.5, а = 1.5, ш = 1.2 Fig. 1.4. Amplitude of forced oscillations of the pendulum \Av(d, Q)| according to (1.31) for the values of the 7 parameter characterizing the friction: a — 7 = 0.2; b — 7 = 1.1; c — 7 = 1.9. The calculation was made for

a = 0.5, a = 1.5, и = 1.2

Рис. 1.5. Амплитуда вынужденных колебаний маятника \Av(d, П)\ согласно (1.31) для значений параметра 7, характеризующего трение: a — 7 = 0.2; b — 7 = 1.1; c — 7 = 1.9. Расчет выполнен

для а = 1.5, а =1.5, ш =1.2 Fig. 1.5. Amplitude of forced oscillations of the pendulum \Av(d, Q)| according to (1.31) for the values of the 7 parameter characterizing the friction: a — 7 = 0.2; b — 7 = 1.1; c — 7 = 1.9. The calculation was made for

a = 1.5, a = 1.5, и = 1.2

Заключение

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

Литература

[1] Зайцев В.В., Федюнин Э.Ю. Модель автоколебаний без гармоник основной частоты // Вестник Самарского университета. Естественнонаучная серия. 2018. Т. 24, № 3. С. 53-59. DOI: http://doi.org/10.18287/2541-7525-2018-24-3-53-59. EDN: https://elibrary.ru/ytbbet.

[2] Ibrahim P., Arafa M., Anis Y. An Electromagnetic Vibration Energy Harvester with a Tunable Mass Moment of Inertia // Sensors. 2021. Vol. 21, Issue 16. Article number 5611. DOI: https://doi.org/10.3390/s21165611.

[3] Le Scornec J., Guiffard B., Seveno R., Le Camb V. Frequency tunable, flexible and low cost piezoelectric micro-generator for energy harvesting // Sensors and Actuators A: Physical. 2020. Vol. 312, Article number 112148. DOI: https://doi.org/10.1016/j.sna.2020.112148.

[4] Wang Y., Li L., Hofmann D., Andrade J.E., Daraio C. Structured fabrics with tunable mechanical properties // Nature. 2021. Vol. 596, pp. 238-243. DOI: https://doi.org/10.1038/s41586-021-03698-7.

[5] Semenov M.E., Solovyov A.M., Meleshenko P.A. Stabilization of coupled inverted pendula: From discrete to continuous case // Journal of Vibration and Control. 2021. Vol. 27, Issue 1-2, Pp. 43-56. DOI: https://doi.org/10.1177/1077546320923436. EDN: https://elibrary.ru/uztjrc.

[6] Semenov M.E., Meleshenko P.A., Borzunov S.V., Reshetova O.O., Barsukov A.I. A Simple Model of the Energy Harvester within a Linear and Hysteresis Approach // Micromachines. 2023. Vol. 14, Issue 2, Article number 310. DOI: https://doi.org/10.3390/mi14020310.

[7] Wang Z., He L., Zhang Z., Zhou Z., Zhou J., Cheng G. Research on a Piezoelectric Energy Harvester with Rotating Magnetic Excitation // Journal of Electronic Materials. 2021. Vol. 50, pp. 3228-3240. DOI: https://doi.org/10.1007/s11664-021-08910-y.

[8] Kovacova V., Glinsek S., Girod S., Defay E. High Electrocaloric Effect in Lead Scandium Tantalate Thin Films with Interdigitated Electrodes // Sensors. 2022. Vol. 22, Issue 11, article number 4049. DOI: https://doi.org/10.3390/s22114049.

[9] Bouhedma S., Hu S., Schutz A., Lange F., Bechtold T., Ouali M., Hohlfeld D. Analysis and Characterization of Optimized Dual-Frequency Vibration Energy Harvesters for Low-Power Industrial Applications // Micromachines. 2022. Vol. 13, Issue 7. Article number 1078. DOI: https://doi.org/10.3390/mi13071078.

10] Damjanovic D. Ferroelectric, dielectric and piezoelectric properties of ferroelectric thin films and ceramics // Reports of Progress in Physics. 1998. Vol. 61, Number 9, Pp. 1267-1324. DOI: https://doi.org/10.1088/0034-4885/61/9/002.

11] Dawber M. Physics of thin-film ferroelectric oxides // Review of Modern Physics. 2005. Vol. 77, Issue 4. Article number 1083. DOI: https://doi.org/10.1103/RevModPhys.77.1083.

12] Cheng C., Peters T., Dangi A., Agrawal S., Chen H., Kothapalli S.-R., Trolier-McKinstry S. Improving PMUT Receive Sensitivity via DC Bias and Piezoelectric Composition // Sensors. 2022. Vol. 22, Issue 15. Article number 5614. DOI: https://doi.org/10.3390/s22155614.

13] Ge C., Cretu E. Simple and Robust Microfabrication of Polymeric Piezoelectric Resonating MEMS Mass Sensors // Sensors. 2022. Vol. 22, Issue 8, Article number 2994. DOI: https://doi.org/10.3390/s22082994.

14] Hu H., Ben Mrad R. On the classical Preisach model for hysteresis in piezoceramic actuators // Mechatronics. 2003. Vol. 13, Issue 2. Pp. 85-94. DOI: https://doi.org/10.1016/S0957-4158(01)00043-5.

15] Renno J.M., Daqaq M.F., Inman D.J. On the optimal energy harvesting from a vibration source // Journal of Sound and Vibration. 2009. Vol. 320, Issues 1-2, pp. 386-405. DOI: https://doi.org/10.1016/jjsv.2008.07.029.

16] Triplett A., Quinn D.D. The effect of nonlinear piezoelectric coupling on vibration-based energy harvesting // Proceedings of the IMECE2008 ASME International Mechanical Engineering Congress and Exposition, Boston, MA, USA, 31 October — 6 November 2008; IMECE2008-66393. Pp. 1-6.

17] Tusset A.M., Rocha R.T., Iliuk I., Balthazar J.M., Litak G. Dynamics and Control of Energy Harvesting from a Non-ideally Excited Portal Frame System with Fractional Damping // Balthazar J.M. (eds) Vibration Engineering and Technology of Machinery. Mechanisms and Machine Science, vol 95. Cham: Springer, 2021, pp. 383-395. DOI: https://doi.org/10.1007/978-3-030-60694-7_25.

18] Pan J., Qin W., Deng W., Zhang P., Zhou Z. Harvesting weak vibration energy by integrating piezoelectric inverted beam and pendulum // Energy. 2021. Vol. 227. Article number 120374. DOI: https://doi.org/10.1016/j.energy.2021.120374.

19] Lygas K., Wolszczak P., Litak G. Broadband frequency response of a nonlinear resonator with clearance for energy harvesting // MATEC Web Conference. 2018. Vol. 148, Article number 12003. DOI: https://doi.org/10.1051/matecconf/201814812003.

20] Friswell M.I., Ali S.F., Bilgen O., Adhikari S., Lees A.W., Litak G. Non-linear piezoelectric vibration energy harvesting from a vertical cantilever beam with tip mass // Journal of Intelligent Material Systems and Structures. 2012. Vol. 23, Issue 13. Pp. 1505-1521. DOI: https://doi.org/10.1177/1045389X12455722.

21] Гантмахер Ф.Р. Теория матриц. Москва: ФИЗМАТЛИТ, 2010. 560 с. ISBN 978-5-9221-0524-8. URL: https://djvu.online/file/9RBsN3mYxlxaq.

DOI: 10.18287/2541-7525-2023-29-2-7-18 Submited: 01.03.2023

Revised: 11.04.2023 Accepted: 30.06.2023

S.V. Borzunov

Voronezh State University, Voronezh, Russian Federation E-mail: sborzunov@gmail.com. ORCID: https://orcid.org/0000-0002-5099-9655

TRANSFORMATION OF OSCILLATIONS OF AN UNSTABLE SYSTEM

IN AN ENERGY HARVESTER

ABSTRACT

In this work, a simple model of energy harvester transforming the energy of unstable mechanical vibrations into useful electric energy is proposed. The mechanical part of the system is presented in the form of an inverted pendulum. The electric part consists of a converter of mechanical energy based on a direct piezoelectric effect, and a payload. The dynamics of the system is considered in the linearized production, the areas of stability are determined, various dynamics modes in the parameter space are identified. It was established that the stabilizing pendulum of management based on the principles of feedback allows you to transfer the system into a stable mode of functioning. The optimal characteristics of the system that meets maximum power were found.

Key words: energy harvester; inverted pendulum; stability of the system.

Citation. Borzunov S.V. Transformation of oscillations of an unstable system in an energy harvester. Vestnik Samarskogo universiteta. Estestvennonauchnaia seriia = Vestnik of Samara University. Natural Science Series, 2023, vol. 29, no. 2, pp. 7-18. DOI: http://doi.org/10.18287/2541-7525-2023-29-2-7-18. (In Russ.)

Information about the conflict of interests: authors and reviewers declare no conflict of interests.

© Borzunov S.V., 2023

Sergei V. Borzunov — Candidate of Physical and Mathematical Sciences, associate professor of the Department of Digital Technologies, Voronezh State University, 1, Universitetskaya Square, Voronezh, 394018, Russian Federation.

References

[1] Zaitsev V.V., Fedyunin E.Yu. Model of self-oscillations without harmonicas of the base frequency. Vestnik Samarskogo universiteta. Estestvennonauchnaia seriia = Vestnik of Samara University. Natural Science Series, 2018, vol. 24, no. 3, pp. 53-59. DOI: http://doi.org/10.18287/2541-7525-2018-24-3-53-59. EDN: https://elibrary.ru/ytbbet. (In Russ.)

[2] Ibrahim P., Arafa M., Anis Y. An Electromagnetic Vibration Energy Harvester with a Tunable Mass Moment of Inertia. Sensors, 2021, vol. 21, issue 16, article number 5611. DOI: https://doi.org/10.3390/s21165611.

[3] Le Scornec J., Guiffard B., Seveno R., Le Camb V. Frequency tunable, flexible and low cost piezoelectric micro-generator for energy harvesting. Sensors and Actuators A: Physical, 2020, vol. 312, p. 112148. DOI: https://doi.org/10.1016/j.sna.2020.112148.

[4] Wang Y., Li L., Hofmann D., Andrade J.E., Daraio C. Structured fabrics with tunable mechanical properties. Nature, 2021, vol. 596, pp. 238-243. DOI: https://doi.org/10.1038/s41586-021-03698-7.

[5] Semenov M.E., Solovyov A.M., Meleshenko P.A. Stabilization of coupled inverted pendula: From discrete to continuous case. Journal of Vibration and Control, 2021, vol. 27, issue 1-2, pp. 43-56. https://doi.org/10.1177/1077546320923436

[6] Semenov M.E., Meleshenko P.A., Borzunov S.V., Reshetova O.O., Barsukov A.I. A Simple Model of the Energy Harvester within a Linear and Hysteresis Approach. Micromachines, 2023, vol. 14, issue 2, p. 310. DOI: https://doi.org/10.3390/mi14020310.

[7] Wang Z., He L., Zhang Z., Zhou Z., Zhou J., Cheng G. Research on a Piezoelectric Energy Harvester with Rotating Magnetic Excitation. Journal of Electronic Materials, 2021, vol. 50, pp. 3228-3240. DOI: https://doi.org/10.1007/s11664-021-08910-y.

[8] Kovacova V., Glinsek S., Girod S., Defay E. High Electrocaloric Effect in Lead Scandium Tantalate Thin Films with Interdigitated Electrodes. Sensors, 2022, vol. 22, issue 11, article number 4049. DOI: https://doi.org/10.3390/s22114049.

[9] Bouhedma S., Hu S., Schutz A., Lange F., Bechtold T., Ouali M., Hohlfeld D. Analysis and Characterization of Optimized Dual-Frequency Vibration Energy Harvesters for Low-Power Industrial Applications. Micromachines, 2022, vol. 13, issue 7, article number 1078. DOI: https://doi.org/10.3390/mi13071078.

10

11

12

13

14

15

16

17

18

19

20

21

Damjanovic D. Ferroelectric, dielectric and piezoelectric properties of ferroelectric thin films and ceramics. Reports of Progress in Physics, 1998, vol. 61, number 9, pp. 1267-1324. DOI: https://doi.org/10.1088/0034-4885/61/9/002.

Dawber M. Physics of thin-film ferroelectric oxides. Reviews of Modern Physics, 2005, vol. 77, issue 4, article number 1083. DOI: https://doi.org/10.1103/RevModPhys.77.1083.

Cheng C., Peters T., Dangi A., Agrawal S., Chen H., Kothapalli S.-R., Trolier-McKinstry S. Improving PMUT Receive Sensitivity via DC Bias and Piezoelectric Composition. Sensors, 2022, vol. 22, issue 15, article number 5614. DOI: https://doi.org/10.3390/s22155614.

Ge C., Cretu E. Simple and Robust Microfabrication of Polymeric Piezoelectric Resonating MEMS Mass Sensors. Sensors, 2022, vol. 22, issue 8, article number 2994. DOI: https://doi.org/10.3390/s22082994.

Hu H., Ben Mrad R. On the classical Preisach model for hysteresis in piezoceramic actuators. Mechatronics. 2003. Vol. 13, Issue 2. P. 85-94. DOI: http://doi.org/10.1016/S0957-4158(01)00043-5.

Renno J.M., Daqaq M.F., Inman D.J. On the optimal energy harvesting from a vibration source. Journal of Sound and Vibration, 2009, vol. 320, issues 1-2, pp. 386-405. DOI: https://doi.org/10.1016/jjsv.2008.07.029.

Triplett A., Quinn D.D. The effect of nonlinear piezoelectric coupling on vibration-based energy harvesting. In: Proceedings of the IMECE2008 ASME International Mechanical Engineering Congress and Exposition, Boston, MA, USA, 31 October-6 November 2008; IMECE2008-66393. pp. 1-6.

Tusset A.M., Rocha R.T., Iliuk I., Balthazar J.M., Litak G. Dynamics and Control of Energy Harvesting from a Non-ideally Excited Portal Frame System with Fractional Damping. In: Balthazar J.M. (eds) Vibration Engineering and Technology of Machinery. Mechanisms and Machine Science, vol. 95. Cham: Springer, 2021, pp. 383-395. DOI: https://doi.org/10.1007/978-3-030-60694-7_25.

Pan J., Qin W., Deng W., Zhang P., Zhou Z. Harvesting weak vibration energy by integrating piezoelectric inverted beam and pendulum. Energy, 2021, vol. 227, article number 120374. DOI: https://doi.org/10.1016/j.energy.2021.120374.

Lygas K., Wolszczak P., Litak G. Broadband frequency response of a nonlinear resonator with clearance for energy harvesting. In: MATEC Web Conference, 2018, vol. 148, article number 12003. DOI: https://doi.org/10.1051/matecconf/201814812003.

Friswell M.I., Ali S.F., Bilgen O., Adhikari S., Lees A.W., Litak G. Non-linear piezoelectric vibration energy harvesting from a vertical cantilever beam with tip mass. Journal of Intelligent Material Systems and Structures, 2012, vol. 23, issue 13, pp. 1505-1521. DOI: https://doi.org/10.1177/1045389X12455722.

Gantmakher F.R. Theory of matrices. Moscow: FIZMATLIT, 2010, 560 p. ISBN 978-5-9221-0524-8. Available at: https://djvu.online/file/9RBsN3mYxlxaq. (In Russ.)