CC BY
34
УДК 539.3:534.1
В.П. Матвеенко, E.II. Клигман, Н.А. Юрлова, Д.В. Грачев Институт механики сплошных сред УрО РАН
УПРАВЛЕНИЕ ДИНАМИЧЕСКИМИ СВОЙСТВАМИ МЕХАНИЧЕСКИХ СИСТЕМ, ВЫПОЛНЕННЫХ ИЗ SMART-МАТЕРИАЛОВ НА ОСНОВЕ ПЬЕЗОЭЛЕКТРИКОВ
Abstract
Application of composites materials involving piezoelectric elements in engineering systems offers considerable scope for optimization of their mechanical characteristics. WV propose to solve this problem by means of embedding piezoelectric elements into the body of shell to achieve vibration suppression through dissipation of energy in external RLC-cin uih. Complex eigenfrequencies obtained consist of the reed part, which is vibration frequeni i. and the imaginary part— damping index. Complex eigenvector contains the ■.nbration mode shape and phase. The approach of piecewise linear normal with the use of kinematic shell theor: for each layer were used to handle the essentially different physical properties of layers Computations performed for shells of arbitrary spatial geometry showed effectiveness of the proposed method of control of dissipative properties.
Новые технологии в производстве композитов породили такое понятие, как SMART-материалы. Основной особенностью таких материалов является их способность целенаправленно изменять свои свойства в зависимости от внешних условий. Как правило, в состав SMART-композитов входят датчики, фиксирующие изменение параметров окружающей среды или термомеханическое состояние конструкции, и активные элементы, гак называемые актуаторы, которые в нужном направлении изменяют механические свойства конструкции. Такой тип SMAR Г-материалов является активным и подразумевает наличие цепей обратной связи и систем управления актуаторами.
Наиболее часто для изготовления SMART-ком позитов используются пьезоэлектрические материалы. В частности, это объясняется наличием у них прямого и обратного пьезоэффекта, что позволяет использовать пьезоэлементы как в качеаве датчиков, так и в качестве актуаторов, В настоящее время наибольшее распространение получили пьезокерамические материалы на основе титаната бария и циркоьата -титаната свинца.
Для решения проблем, связанных с вибрациями конструкций, нами предлагается использовать полуактивные SMART-материалы на базе пьезокерамики, в которых функции датчиков и актуаторов совмещены. Для этого электроднровянные поверхности пьезоэлементов, расположенных в определенных местах конструкции, соединяются пассивными ^¿С-цепями с точкой нулевого потенциала. Принцип гашения вибраций основан на преобразовании механической энергии в электрическую с ее последующем рассеиванием во внешних С-цепях в виде тепла и электромагнитного излучения.
Наибольший эффект от использования полуактивных SMART-материалов достигается при решении задач оптимизации динамических свойств оболочечных элементов конструкций.
Диссипативные свойства конструкции мотуг быть определены из решения спектральной задачи теории злекгровязкоупругости. Для большинства конструкционных материалов спектральная задача теории электровязкоупругости может быть сформулирована в форме с комплексными динамическими модулями, являющимися аналогами дифференциальных и интегральных операторов в физических соотношениях (работы A.C. Кравчука, В.П. Матвеенко, Е.П. Клигмана, В.В. Болотина, Ю.Н. Новичкова), и комплексными сопротивлениями внешних электрических RLC-цепей.
В этой постановке найденные комплексные собственные значения будут определять частоты собственных колебаний и показатели демпфирования, Комплексные частоты и формы собственных колебаний зависят не только от конструктивных и технологических особенностей оболочечных систем, но и от типа и способа электрического соединения (RLC- параметры) пьезоэлементов конструкции. Целенаправленное изменение этих параметров позволит оптимизировать динамические характеристики (резонансные частоты, формы собственных колебаний и показатели демпфирования) системы
Вариационное уравнение движения тела, обладающего пьезоэлектрическими свойствами
Вариационное уравнение движения тела, состоящего из упругого F, и пьезоэлектрического V2 элементов, может быть получено на основе соотношений линейной теории упругости и квазистатических уравнений Максвелла. В случае изотермического процесса и при отсутствии внешних сил вариационное уравнение может быть записано следующим образом:
|(&} 8е>;- + ри 8н,- У V + Jö£,; - D' ЪЕ, + ри Sw. 'jtlV - j qeЩс1й. = 0, (1)
п ’ Vi
где а. г - тензоры напряжений и деформаций; и - вектор перемещений; D. Е -векторы электрической индукции и напряженности электрического поля; р -плотность материала; £1р - поверхность, ограничивающая пьезоэлектрический элемент; qs и ф - поверхностная плотность зарядов и электрический потенциал.
Компоненты тензора деформаций удовлетворят соотношениям Коши
е;/ =^-(н +u,f). Электрическое поле считается потенциальным, то есть выполняется
Зф Г п
условие —- = £,■ Для изотермических процессов в линеиных электроупругих средах
Эл-1
справедливы следующие физические соотношения:
<fJ = С&'% - для F,, (2а)
а1' =Сик1єк1 -Р «кЕ,
где С°к> - тензор упругих констант; Р'/А и ек' - тензоры пьезоэлектрических и
диэлектрических коэффициентов.
Если элемент тела V, обладает вязкоупругими свойствами, то согласно
Применение пьезоматериалов для создания оболочечных конструкций требует использования моделей оболочек, учитывающих все компоненты тензоров напряжении и деформаций. Это объясняется тем, что такие материалы, как правило, поляризованные в направлении толщины, проявляют свои пьезоэлектрические свойства при трансверсальном деформировании. Учет всех компонент тензоров напряжении и деформаций может быть осуществлен в рамках шестимодальной теории оболочек Тимошенко.
Обобщая кинематггческие гипотезы шестимодальной теории Тимошенко, приходим к следующим соотношениям для компонент вектора перемещений и электрического потенциала точки оболочки с Гауссовыми координатами .г'1 (о = 1.2)
и нормальной координатой з"‘ = г:
Здесь иа- перемещения в направлении базисных векторов Гауссовых координат, и, -перемещения в направлении нормали, (р - электрический потенциал в точке: Vг. и . /
- перемещения и потенциал точек поверхности приведения,, - углы поворота нормали, г) - обжатие по толщине оболочки. Ф - изменение потенциала по толщине оболочки.
При получении физических соотношений мы получили обобщенные векторы напряжений и деформаций, которые помимо механических величин содержат и электрические. Введем теперь векторы обобщенных перемещений и нагрузок.
Вариационное уравнение движения пьезоэлектрической оболочки в матричной форме будет иметь следующий вид:
принципу Вольтерра тензор упругих констант С"А/ в зависимости (2а) должен быть заменен соответствующим вязкоупругим оператором.
Уравнения электроупругосги оболочек с внешними Ш( пенями
Я8{Е}' ММ + 8{и}7 [р]{«}) ¿п = ]5{и}' {/} СІП.
г>
где {м}={м1(и2,ф}г, [р] = (1іа£(р,р,р,0], {г} = ,г2,/-3,‘ : -
вектор
обобщенных деформаций. Ь - дифференциальный оператор.
Квазигармонические колебания оболочек
Рассмотрим процессы, которые описываются временной функцией
/(') = /<> *"**■
где СО - комплексная частота;,/« - комплексная амплитуда.
Тогда вариационное уравнение (4) для амплитудных значений компонент состояния примет вид (временной множитель опущен):
|(5{е}г [2>]{е} -0)28{«}г [р]{и})сШ = ]8{«}г .
п а
Интеграл по поверхности в уравнении (4) определяет электрические граничные условия для пьезоэлектрических участков оболочки.
Рассмотрим случай, когда один из пьезоэлементов с злектродированной поверхностью 0{| соединен с точкой нулевого потенциала проводником сопротивлением Я, емкостью С и индуктивностью Для математического описания этого типа граничных условий воспользуемся законом Ома для переменного тока:
Ф=^ + д/ + 1/=^/гё + 1ё. (5)
Здесь <р - потенциал и Q~j дс!£1 - заряд на электроде; / = <2 - ток в проводнике.
«с
Так как мы рассматриваем квазигармонический процесс, то есть ф{7) = ф(1 е'"'1', то дифференциальное уравнение (5) можно разрешить относительно
Фп ■ £ к"'
0(0=-., г"~......*■
С -со ■ I - г • со- К
С учетом эквипотенциальности поверхности П, и того, что КЬС является внешней цепью, будет выполняться условие
I С -со I- ¿сой
При отсутствии иных поверхностных сил. окончательно уравнение движения примет вид
}(8{е}г [д]{е}-шг5{ы}Г мм) ¿п+ 5ф» о, (6,
’ С - со ь -1 ож
где ф0 - потенциал на злектродированной части пьезоэлемента П Уравнение (6)
является однородным и может рассматриваться как вариационная задача на собственные значения.
Конечно - элементная реализация поставленной задачи
Вариационная задача (6) с помощью метода конечных элементов (МКО) сводится к алгебраической проблеме собственных значений:
([к]~со2[Лф[е(а))]){Х}:=0, (7)
где [А'1 - матрица жесткости, [М] - матрица масс. [С] ~ матрица "жесткости'
внешних ^¿С-цепей,
Условием существования нетривиального решения является
О(со) = ёег ([А'] - со‘ [А/] + [(7 (со)]) = 0 .
Уравнение (7) существенно отличается от обобщенной проблемы на
собственные значения наличием матрицы [С(со) |, которая содержит нули и элементы
вида
С“ - п) А - /О)/?
Для расчета слоистых оболочек использован следующий подход, позволяющий учитывать изменение свойств материала по толщине.
Представим рассматриваемую конструкцию совокупностью п независимых однослойных оболочек. Выберем в качестве единой поверхности приведения (координатной поверхности) внутреннюю поверхность первого слоя Б’'. Построим на поверхности конечно-элементную сетку с одинаковой аппроксимацией всех
компонент вектора перемещений.
Для каждого к-то слоя вводится свой вектор узловых параметров
{ак щ ,\|/£,%к ,Фа)Г={{^Л)>Ф*}" ^ определяется соответствующая ему
матрица жесткости [АГ^]. Здесь {{/*} включает перемещения координатной поверхности ик, и 'к , мк, а компонента {Ч/*}— углы поворота нормали и<т, и обжагие по толщине Хк • Преобразуем вектор узловых перемещений ¿-го узла для А-ю слоя к новым компонентам:
ик
Ф*
[1] О! [I] 0'
1
•І А1>. Ф і
где и к - вектор перемещений внутренней поверхности слоя; Л(Л - приращение вектора перемещений в к-м слое; [1] и [0] - единичная и нулевая матрицы соответствующей размерности.
В результате перехода к вектору неизвестных в приращениях матрица жесткости будет определяться следующим образом:
к;]=кГку]И 1-
Условием совместности деформации первого и второго слоя является соотношение:
Ы={ЦМД£А}-
В матричном виде оно приводится ниже. Здесь учтено, что первый слой не обладает пьезоэлектрическими свойствами (узловые переменные: Ъ\, АС/|), а второй является пьезоэлектриком (С/2, АСА, ф):
и, '[1] 0 0 о"
шх 0 [1] 0 0
Щ [1] [1] 0 0
д и2 0 0 [1] 0
ф о 0 0 1_
' их '
АС/,
А и2
ф
= [т2] {Щ},
вектор узловых параметров
К'з]
где [Т2] - матрица «склейки» 1-го и 2-го слоев.
Данное уравнение вводит новый
описывающий совместную деформацию двух первых слоев. При этом матрицы жесткости этих слоев объединяются следующим образом:
О
[4] = ^]
Т\Ц 10
КІ
[Тг]-
Теперь - матрица жесткости двух первых объединенных слоев. Аналогично
вычисляются матрицы жесткости для любого числа слоев. Матрица масс слоистой оболочки и векторы узловых нагрузок получаются но этой же схеме.
Численное исследование и оптимизация динамических свойств конструкции
Для иллюстрации возможностей предложенного подхода был выполнен расчсч н проведена оптимизация демпфирующих свойств стальной квадратной пластинки размером 120x120x0,5 мм, защемленной по контуру. Предполагалось найти место расположения пьезоэлемента и величину сопротивления шунтирующего резистора (емкость и индуктивность отсутствуют) для наилучшего гашения третьей моды колебаний (рис. 1). Пьезоэлемент был выполнен из керамики J,Z^5 и имел размер 15x15x0,5 мм.
Производился поиск оптимальных областей расположения и величин сопротивления для третьей моды колебаний. График зависимости показателя демпфирования (мнимой части
комплексных собственных частот) от места расположения пьезоэлемента на пластинке
Рис. 1, Третья форма колебаний пластинки
для третьей моды показан на рис.2. Зависимость показателя демпфирования от
величины сопротивления шунтирующего резистора Л в случае установки
пьезоэлемента в оптимальном месте приведена на рис.З.
да._
А?
Рис.2. Зависимость показателя демпфирования от места расположения пьезоэлемента на плас тинке для третьей моды
Рис.З. Зависимость показателя демпфирования от величины шунтирующего сопротивления Л'
При расположении пьезоэлемента в оптимальной зоне показатель демпфирования может быт ь увеличен более чем в 2 раза.
Для оценки снижения резонансных амплитуд третьей моды колебаний были построены амплитудно-частотные характеристики (АЧХ) различных вариантов пластинки с оптимальным расположением пьезоэлемента при возбуждении колебании сосредоточенной силой Р - С05{р?), приложенной на поверхности пластинки. Участок АЧХ в окрестности третьего резонанса показан на рис. 4.
Рис. 4. Амплитудно-частотная характеристика системы с оптимальным демпфированием третьей моды колебаний
Кривые на графике соответствуют следующим вариантам:
- 1 - колебания упругой пластинки;
- 2 - колебания пластинки, обладающей вязкоупругими свойствами с
тангенсом потерь 5 = 0^001;
- 3 - колебания пластинки при оптимальном значении сопротивления
шунтирующего резистора R=R0;
- 4 - колебания пластинки при значении сопротивления шунтирующего
резистора R=R0/W;
- 5 - колебания пластинки при значении сопротивления шунтирующего
резистора R=10Ro,
Полученные результаты позволяют сделать вывод о том, что оптимальное значение сопротивления шунтирующего резистора не зависит от места установки пьезоэлемента, а определяется резонансной частотой со, на которую настроен демпфер. При этом выполняется условие равенства сопротивления внешней цепи и емкостного сопротивления пьезоэлемента R0 =1/(сйС), где С = cc()S¡el - емкость пьезоэлемента; S - плошадь электродированной поверхности; d - расстояние между электродами; е -относительная диэлектрическая проницаемость пьезокерамики; е0 - диэлектрическая проницаемость вакуума.
Этот факт хорошо известен в электротехнике и говорит о том, что максимальная отдача энергии реализуется при равенстве внутреннего сопротивления источника и нагрузки.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант РФФИ № 00-01-640).
Библиографический список
1. Matveyenko V.P., Kligman Е.Р. Natura! Vibration Problem of Viscoelastic Solids as Applied to Optimization of Dissipative Properties of Constructions // Journal of Vibration and Control - 1997. - 3. - P. 87-102.
2. Болотин В.В., Новичков Ю.Н. Механика многослойных конструкций. - М.: Машиностроение, 1980. - 376 с.
3. Партон В.З., Кудрявцев Б.А.. Электромагнитоупругость пьезоэлектрических и электропроводных тел. - М.; Наука, 1988. - 472 с.
4. Васидзу К. Вариационные методы в теории упругости и пластичности. - М.; Мир. 1987. - 542 с.
5. Кравчук A.C., Майборода В.П., Уржумцев Ю.С. Механика полимерных и композиционных материалов. - М.: - Наука, 1985. - 304 с.
Получено 20.03.2001
