exponenta event banner

Пьезо Бендер

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

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

Описание

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

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

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

Уравнения

На этом рисунке показана декартова опорная рамка

где:

  • L - длина луча.

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

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

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

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

где:

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

  • [c] - тензор соответствия.

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

  • [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) коэффициент связи пьезоэлектрического напряжения с зарядом, startxx = e31Ey.

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

Mz (x) =−∬yTxx (x, y) dS=∬Ey2dφzdx (x) + e31vdy dS.

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

Mz (x) = E12w d3doftzdx (x) + w d e31v = EIdoftzdx (x) + w d e31v,

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

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

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

Dy (x, y) = e31ydoftzdx (x) +ϵvd.

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

dq (x) =∯Dy dS = e31dw dxdoftzdx = e31dw díz.

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

q = e31dw (фз (х2) фз (х1)).

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

Mz (x) = EIdízdx (x) + w d e31v

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

q = e31d w (фzR фzC) +ϵwldv.

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

Для дискретизации и решения уравнений Эйлера - Бернулли с пьезоэлектрической муфтой блок Пьезо Бендера использует метод конечных элементов.

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

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

[FCTCFRTR] = [12IEl36EIl2 12EIl36EIl26EIl24EIl 6EIl22EIl 12EIl212EIl3 6EIl26EIl22EIL L 6

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

αlwd420 [15622l54 13l22l4l213l 3l25413l156 22l − 13l 3l2 22l4l2] [d2yCdt2d2φCdt2d2yRdt2d2φRdt20]+[12EIl36EIl2−12EIl36EIl206EIl24EIl−6EIl22EIle31wd−12EIl3−6EIl212EIl3−6EIl206EIl22EIl−6EIl24EIl−e31wd0−e31wd0e31wdϵwld] [yCФРRv] = [FCTCFRTRq].

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

αlwd420 [15622l54 13l22l4l213l 3l25413l156 22l − 13l 3l2 22l4l2] [d2yCdt2d2φCdt2d2yRdt2d2φRdt20] + [[B]⋯0⋮⋱⋮0⋯0] [dyCdtdφCdtdyRdtdφRdt0]+[12EIl36EIl2−12EIl36EIl206EIl24EIl−6EIl22EIle31wd−12EIl3−6EIl212EIl3−6EIl206EIl22EIl−6EIl24EIl−e31wd0−e31wd0e31wdϵwld] [yCФРRv] = [FCTCFRTRq]

где:

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

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

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

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

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

  • m = αlwd - масса элемента, где "" "" "" "" "" "" "" "" "" "" "" "" "" "

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

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

  • B = bkK + bmM - демпфирующая матрица.

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

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

  • K = [12EIl36EIl2 12EIl36EIl26EIl24EIl 6EIl22EIl 12EIl3 6EIl212EIl3 6EIl26EIl22EIl − 6EIl24EIl] - матрица конечных элементов жёсткости.

  • M = startlwd420 [15622l54 13l22l4l213l 3l25413l156 22l 13l − 3l2 − 22l4l2] - матрица масс.

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

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

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

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

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

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

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

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

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

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

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

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

  • Размеры (l, w, d)

  • Масса, м

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

  • Свободное отклонение при номинальном напряжении, без

  • Сила блокировки при номинальном напряжении, Fблок

  • Емкость, Cpiezo

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

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

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

[12EIl3 6EIl20 6EIl24EIl e31w d0e31w d

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

yR=l2 (4FR l+6 d w e31 v) 12EI=l2 (4FR l+6 d w e31 v) E w d3=FR4l3Ewd3+v6 e31l2E d2yfree=6l2 d w e31 vrated12EI=6 e31l2 vratedE d2Fblock =−3 d w e31 vrated2l

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

Cpiezo = αl wd + e312l d2w2EI = αl wd + e31212l wEd = (start+ 12e312E) l wd.

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

m = startlwd.

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

e31 =−2 l Fblock3 d w vratedE =−4 Fblockl3yfreed3wε = dl w (Cpiezo+4Fblockyfreevrated2)

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

vvrated=FRFblock+yRyfreeq=Cpiezov+yfreevratedFR

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

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

2πf0=1.8552EIm l3.

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

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

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

  • Свободно (Free) - как смещение, так и вращение равны любому значению.

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

  • Зажатый - как смещение, так и вращение равны 0.

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

КонфигурацияМодель
Без зажима

Поддерживается-поддерживается

Зажатый-зажатый

Порты

Сохранение

развернуть все

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

Электрический консервационный порт, связанный с отрицательным выводом пьезо-бендера.

Механический консервационный порт, связанный с корпусом.

Механическое поворотное защитное отверстие, связанное с корпусом.

Механический консервационный порт, связанный с ротором.

Механическое отверстие для сохранения вращения, связанное с ротором.

Параметры

развернуть все

Размеры

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

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

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

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

Стационарный режим

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

Емкость пьезобогатителя.

Зависимости

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

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

Зависимости

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

Отклонение наконечника при отсутствии приложения силы. Знак этого параметра должен отличаться от знака силы блокировки на параметре Vrated.

Зависимости

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

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

Зависимости

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

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

Зависимости

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

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

Зависимости

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

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

Зависимости

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

Динамика

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

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

Зависимости

Чтобы включить этот параметр, установите для параметра Инерционные эффекты значение On.

Масса балки.

Зависимости

Чтобы включить этот параметр, установите для параметра Инерционные эффекты значение On и параметризация для Specify from material properties.

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

Зависимости

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

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

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

Зависимости

Чтобы включить этот параметр, установите для параметра Инерционные эффекты значение On.

Ссылки

[1] Тадмор, Э. Б. и Г. Коса. «Коррекция электромеханической муфты для пьезоэлектрических многослойных балок». Журнал микроэлектромеханических систем, том 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] Гэвин, Анри П. «Балочный элемент Stiffness Матрицы». 421L ЦВЕ. Матричный структурный анализ. Университет Дьюка, 2014.

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

Создание кода C/C + +
Создайте код C и C++ с помощью Simulink ® Coder™

.
Представлен в R2021a