Piezo Bender

Пьезоэлектрический биморфный брус прямоугольного сечения

  • Библиотека:
  • Simscape/Электрический/Электромеханический/Мехатронный Приводы

Описание

Блок Piezo Bender моделирует пьезоэлементный пучок биморфа прямоугольного сечения.

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

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

Уравнения

Этот рисунок показывает декартовую систему координат

где:

  • L - длина балки.

  • w - ширина луча.

  • d - толщина балки.

Это конститутивные уравнения для пьезоэлектрического материала в композиции стресс-заряда,

T=[c]S[e]TE(1)
D=[e]S[ϵS]E(2)

где:

  • T - поле напряжения.

  • [c] - тензор податливости.

  • S - поле деформации.

  • [e] - тензор коэффициента пьезометрического напряжения.

  • E - электрическое поле.

  • D - поле электрического смещения.

  • S] - тензор диэлектрической проницаемости при постоянной или нулевой деформации.

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

Этот блок рассматривает только силы, приложенные в y направлении, и пьезоэлектрический материал поляризован, чтобы сгибаться только в x-y плоскости. По этой причине, чтобы описать кинематику, необходимо только задать вертикальное перемещение в y направлении центра тяжести каждого поперечного сечения, y(x) и вращение вокруг оси z каждого поперечного сечения, φz(x).

Из предыдущих предположений поле деформации в луче Эйлера-Бернулли, подлежащем изгибу, равно:

Sxx(x,y)=ydφzdx(x).(3)

Потому что электрическое поле постоянно между положительной и отрицательной пластинами, Ey=vdблок подставляет уравнение 3 в уравнение 1:

Txx(x,y)=Eydφzdx(x)e31vd.

В этом уравнении E = c11 является модулем Янга материала и e31 является (3,1) пьезоэлектрическим коэффициентом связи напряжение-заряд, σxx=e31Ey.

Это уравнение задает изгибающий момент из поля напряжения:

Mz(x)=yTxx(x,y)dS=Ey2dφzdx(x)+e31vdydS.

Потому что поляризация материала для y=[d2,0] является противоположным, чем поляризация для y=[0,d2](3,1) пьезоэлектрический коэффициент связи напряжение-заряд изменяет знак, и изгибающий момент определяется как

Mz(x)=E12wd3dφzdx(x)+wde31v=EIdφzdx(x)+wde31v,

где I=112wd3 - второй момент площади прямоугольного поперечного сечения.

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

Затем блок подставляет уравнение 3 в уравнение 2:

Dy(x,y)=e31ydφzdx(x)+ϵvd.

Электрический заряд внутри объема равен интегралу Гаусса электрического перемещения:

dq(x)=DydS=e31dwdxdφzdx=e31dwdφz.

Затем это уравнение определяет заряд, накопленный между двумя участками луча из-за пьезоэлектрического эффекта:

q=e31dw(φz(x2)φz(x1)).

Наконец, с механической точки зрения можно смоделировать пьезогенератор как балку Эйлера-Бернулли, нагруженную равномерным крутящим моментом, пропорциональным напряжению:

Mz(x)=EIdφzdx(x)+wde31v

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

q=e31dw(φzRφzC)+ϵwldv.

Формулировка конечного элемента

Чтобы дискретизировать и решить уравнения Эйлера-Бернулли с пьезоэлектрическим связыванием, Piezo Bender блок использует метод Конечного Элемента.

Блок дискретизирует пьезогенерирующую балку на ряд срезов в направлении длины с той же шириной, w и толщиной, d. Длина каждого элемента равна общей длине балки, разделенной на количество элементов, l=LNelements.

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

[FCTCFRTR]=[12EIl36EIl212EIl36EIl26EIl24EIl6EIl22EIl12EIl36EIl212EIl36EIl26EIl22EIl6EIl24EIl][yCφCyRφR].

Затем, чтобы получить уравнения для пьезогенерирующего элемента луча, добавьте условия связи и большую матрицу для инерции:

ρlwd420[15622l5413l22l4l213l3l25413l15622l13l3l222l4l2][d2yCdt2d2φCdt2d2yRdt2d2φRdt20]+[12EIl36EIl212EIl36EIl206EIl24EIl6EIl22EIle31wd12EIl36EIl212EIl36EIl206EIl22EIl6EIl24EIle31wd0e31wd0e31wdϵwld][yCφCyRφRv]=[FCTCFRTRq].

Наконец, это уравнение для пьезоизгибающего элемента балки с демпфированием:

ρlwd420[15622l5413l22l4l213l3l25413l15622l13l3l222l4l2][d2yCdt2d2φCdt2d2yRdt2d2φRdt20]+[[B]000][dyCdtdφCdtdyRdtdφRdt0]+[12EIl36EIl212EIl36EIl206EIl24EIl6EIl22EIle31wd12EIl36EIl212EIl36EIl206EIl22EIl6EIl24EIle31wd0e31wd0e31wdϵwld][yCφCyRφRv]=[FCTCFRTRq]

где:

  • l - длина элемента.

  • w - ширина элемента.

  • d - толщина элемента.

  • I=112wt3 - второй момент площади.

  • E - модуль Янга.

  • m = ρlwd - масса элемента, где ρ - плотность массы.

  • e31 - (3,1) пьезоэлектрический коэффициент связи напряжение-заряд, σxx=e31Ey.

  • ε - электрическая диэлектрическая проницаемость.

  • B=bkK+bmM - матрица демпфирования.

  • bm - коэффициент Релея демпфирования, пропорциональный массе.

  • bk - коэффициент Релея демпфирования, пропорциональный жесткости.

  • K=[12EIl36EIl212EIl36EIl26EIl24EIl6EIl22EIl12EIl36EIl212EIl36EIl26EIl22EIl6EIl24EIl] - матрица конечного элемента жесткости.

  • M=ρlwd420[15622l5413l22l4l213l3l25413l15622l13l3l222l4l2] - большая матрица.

  • yC - отклонение вдоль y оси левого конца элемента .

  • yR - отклонение вдоль оси y правого конца элемента .

  • φC - поворот вокруг z оси левого конца элемента .

  • φR - поворот вокруг z оси правого конца элемента .

  • FC - сила вдоль оси y левого конца элемента .

  • FR - сила вдоль оси y правого конца элемента .

  • TC - крутящий момент в z -оси левого конца элемента .

  • TR - крутящий момент в z -оси правого конца элемента .

  • v - напряжение на верхнем и нижнем электродах.

  • q - накопленный заряд между электродами и пьезоэлектрическим материалом.

Параметризация таблицы данных

Таблица данных пьезобендера обычно предоставляет следующие данные:

  • Размерности (l, w, d)

  • Масса, m

  • Номинальное напряжение, vrated

  • Свободный прогиб при номинальном напряжении, yfree

  • Блокирующая сила при номинальном напряжении, Fblock

  • Емкость, Cpiezo

  • Первая резонансная частота, f1

Можно вычислить основные параметры материала пьезогенератора с помощью параметров таблицы данных.

Во-первых, блок решает зависимости отклонения силы напряжения из установившихся уравнений без приложенного крутящего момента и без зажима строения:

[12EIl36EIl206EIl24EIle31wd0e31wdεwld][yRφRv]=[FR0q].

Эти уравнения определяют отношение между отклонением совета, напряжением и силой совета:

yR=l2(4FRl+6dwe31v)12EI=l2(4FRl+6dwe31v)Ewd3=FR4l3Ewd3+v6e31l2Ed2yfree=6l2dwe31vrated12EI=6e31l2vratedEd2Fblock=3dwe31vrated2l

Блок вычисляет емкость, принимая нулевую приложенную силу:

Cpiezo=εlwd+e312ld2w2EI=εlwd+e31212lwEd=(ε+12e312E)lwd.

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

m=ρlwd.

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

e31=2lFblock3dwvratedE=4Fblockl3yfreed3wε=dlw(Cpiezo+4Fblockyfreevrated2)

Затем подставьте эти уравнения в конститутивные уравнения:

vvrated=FRFblock+yRyfreeq=Cpiezov+yfreevratedFR

Таблица данных динамики

Можно вычислить первую резонансную частоту безжимного луча равномерного сечения при помощи этого уравнения:

2πf0=1.8552EIml3.

Затем блок параметрирует динамику непосредственно путем определения желаемой собственной частоты.

Граничные условия

Балки имеют различные граничные условия на левом и правом концах:

  • Свободно - и перемещение, и вращение равны любому значению.

  • Просто поддерживается - перемещение равно 0.

  • Зажим - и перемещение, и вращение равны 0.

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

ConfigurationМодель
Без зажима

Поддерживаемые

Зажимно-зажимной

Порты

Сохранение

расширить все

Электрический порт сопоставлен с положительной клеммой пьезоизгиба.

Электрический порт сопоставлен с отрицательной клеммой пьезоизгиба.

Порт механической передачи, связанная с делом.

Механический вращательный порт сопоставлен с корпусом.

Порт механической передачи, сопоставленный с ротором.

Механический вращательный порт сопоставлен с ротором.

Параметры

расширить все

Размерности

Количество элементов, в которых пьезогенерирующий луч дискретизирован с помощью уравнений конечных элементов Эйлера-Бернулли.

Общая длина балки.

Ширина балки.

Толщина балки.

Установившееся состояние

Параметризация параметров устойчивого состояния из таблицы данных или из свойств материала.

Емкость пьезогенератора.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from a datasheet for clamped-free configuration.

Номинальное напряжение привода.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from a datasheet for clamped-free configuration.

Отклонение совета, когда никакая сила не приложена. Знак этого параметра должен отличаться от знака параметра Blocking force at Vrated.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from a datasheet for clamped-free configuration.

Сила, которая обнуляет отклонение совета. Знак этого параметра должен отличаться от знака параметра Free deflection at Vrated.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from a datasheet for clamped-free configuration.

Молодой модуль, или модуль упругости в напряжении.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from material properties.

(3,1) пьезоэлектрический коэффициент связи напряжение-заряд.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from material properties.

Диэлектрическая проницаемость, которая представляет относительную диэлектрическую проницаемость материала. Этот параметр должен быть больше, чем вакуумная диэлектрическая проницаемость.

Зависимости

Чтобы включить этот параметр, установите Parameterization равным Specify from material properties.

Динамика

Моделировать ли инерционные эффекты.

Параметризация параметров инерции из таблицы данных или из свойств материала.

Зависимости

Чтобы включить этот параметр, установите Inertial effects равным On.

Масса балки.

Зависимости

Чтобы включить этот параметр, установите Inertial effects равным On и Parameterization к Specify from material properties.

Первая естественная частота.

Зависимости

Чтобы включить этот параметр, установите Inertial effects равным On и Parameterization к Specify from a datasheet for clamped-free configuration.

Коэффициент демпфирования, пропорциональный жесткости, для демпфирования Релея.

Коэффициент демпфирования, пропорциональный массе, для демпфирования Релея.

Зависимости

Чтобы включить этот параметр, установите Inertial effects равным On.

Ссылки

[1] Tadmor, E. B., and G. Kosa. «Электромеханическая коррекция связи для пьезоэлектрических слоистых пучков». Журнал микроэлектромеханических систем, том 12, № 6, декабрь 2003, стр. 899-906. DOI.org (Crossref), doi:10.1109/JMEMS.2003.820286.

[2] Benjeddou A, Trindade MA, Ohayon R. «Единая модель конечного элемента луча для механизмов продления и сдвига пьезоэлектрических приводов». Журнал интеллектуальных материальных систем и структур. 1997;8(12):1012-1025. doi:10.1177/1045389X9700801202

[3] Гэвин, Анри П. «Beam Element Stiffness Matrices». 421L CEE. Матричный структурный анализ. Университет Дьюка, 2014.

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ Simulink ®

.
Введенный в R2021a