CC BY
15
ВЯЗКОУПРУГИЕ ПЛАСТИНОВЫЕ МЕХАНИЧЕСКИЕ СИСТЕМЫ С ТОЧЕЧНЫМИ СВЯЗЯМИ И ИХ СОБСТВЕННЫЕ КОЛЕБАНИЯ
Салиева Олима Камаловна
канд. техн. наук, доцент, Бухарский инженерно-технологический институт,
Узбекистан, г. Бухара E-mail: saliyevaok@mail. ru
Кулдашева Фируза Салимовна
ассистент,
Бухарский инженерно-технологический институт,
Узбекистан, г Бухара E-mail: saliyevaok@mail. ru
VISCOELASTIC PLASTIC MECHANICAL SYSTEMS WITH SPOT LINKS AND THEIR OWN VIBRATIONS
Olima Saliyeva
candidate of technical Sciences, associate Professor Bukhara engineering and technological Institute,
Uzbekistan, Bukhara
Firuza Kuldasheva
Assistant
Bukhara engineering and technological Institute,
Uzbekistan, Bukhara
АННОТАЦИЯ
В данной работе рассматриваются собственные колебания вязкоупругих пластинчатых механических систем с точечными связями. Получены частотные уравнения и решено численно, методом Мюллера. Приведен параметрический анализ комплексных собственных частот в зависимости от геометрических параметров.
ABSTRACT
In this article, the natural oscillations of viscoelastic lamellar mechanical systems with point connections are considered. Frequency equations are obtained and solved numerically by the Muller method. A parametric analysis of complex eigenfrequencies depending on the geometric parameters is given
Ключевые слова: Свободные колебания, диссипативная система, собственные колебания, вязкоупругие системы.
Keywords: Free oscillations, dissipative system, vibrations, viscoelastic system.
Структурная неоднородность системы определяется наличием в ней вязкоупругих элементов с разными диссипативными свойствами (в противном случае это структурно - однородная вязкоупругая система). Под механической системой здесь понимается прямоугольная пластина, пакет прямоугольных пластин, оболочка вращения, система оболочек вращения, имеющие точечные связи.
Свободные колебания диссипативной системы носят затухающий характер. Амплитуды форм колебаний с течением времени уменьшаются, поэтому такой процесс, строго говоря, не является периодическим. Но частоты соответствующих форм при этом остаются постоянными [1, 2] и в этом смысле диссипативную систему можно исследовать как систему, обладающую собственными колебаниями.
Постановка задачи. Рассмотрим механическую систему, состоящую из N изотропных вязкоупругих тел, занимающих объем Ц1 и ограниченных поверхностями Пп (п = 1, ...,Ю. При этом предполагается, что один линейный размер каждого тела намного меньше двух остальных. При каждом п на части поверхности п-го тела заданы П^ однородные граничные условия, на остальной свободной поверхности Пв = Пи/Пгр в конечном числе точек наложены связи кинематического и динамического характера: точечные жесткие, упругие и (или) вязкоупругие шарнирного типа опоры (жесткие опоры могут быть защемлены), жесткие упругие и (или) вязкоупругие амортизаторы, соединяющие тела (при N > 1), сосредоточенные массы
Библиографическое описание: Салиева О.К., Кулдашева Ф.С. Вязкоупругие пластиновые механические системы с точечными связями и их собственные колебания //Universum: Технические науки: электрон. научн. журн. 2019. J№ 6(63). URL: http://7universum.com/ru/tech/archive/item/7530
Мчп (Я = 1,—^). Расположение связей и масс на поверхностях П.П произвольно.
В общем случае диссипативные свойства элементов системы различны. Частным случаем такой структурно неоднородной вязкоупругой системы является система с упругими и вязкоупругими элементами. Для последнего случая N = + Ип, где Иу-количество упругих элементов системы, Nп- количество вязкоупругих элементов. При N > 2 тела расположены параллельно друг другу свободными поверхностями (пакеты пластин оболочек). При N = 1 стойки отсутствуют. Требуется определить частоты собственных колебаний вязкоупругой системы, а также оценить ее демпфирующую способность. В математической постановке вязкоупругость выглядит следующим образом. Пусть все точки п-го тела подчиняются гармоническому закону колебаний, т.е.
UnJ(xn,t) = U0j(xn)e-iwt,n = 1,.
■N,j = l...../,(1)
где и^(хп)- ]-я компонента вектора перемещений п-го тела, I- число компонент вектора перемещений, хп = (хП, хЩ, хП)-радиус-вектор точки п-го тела, ш = шк + ш1 - искомая комплексная частота системы, причем шк - собственная частота, а ш1 -коэффициент демпфирования (ш1 < 0). Поскольку каждая компонента вектора перемещений уже имеет индекс п, то последний для обозначения компонент радиус-вектора в дальнейшем не используется. Для прямоугольных пластин = 1 и
и^1(х1,х2) = Щ0(х,у), для оболочек вращения / = 3 и
и01(х1,х2) = и0(х,у),и02(х,у) = Упо(х,у),и03(х,у) =Шп0(х,у),
где х,у— координаты. Исходя из принципа возможных перемещений, приравняем нулю сумму работ всех активных сил, включая силы инерции на возможных перемещениях 8ип^ (х, Ь):
где рп,Уп-плотность и объем п-го тела, Мцп- q-я присоединения масса п- го тела с координатами
хч = (хч хч 3) 1
- число пружин (амортизаторов) между п-м и (п+1)-м телами, Qп - число сосредоточенных масс на п-м теле, 1п - число упругих (вязкоупругих) опор на
п п
п-м теле, Отк,£тк,а1 , н , а1/,£1/ - компоненты тензоров напряжений и деформаций соответственно п-го тела, 1-й пружины (амортизатора) и 1/-й упругой (вяз-коупругой) опоры.
Физические и геометрические соотношения для упругого элемента или упругой связи системы запишем с помощью обобщенного закона Гука
(t) = ^ п® ( ^^т к + 2№п^тк (t),
где Лп, Дп-интегральные операторы Вольтерра, которые ниже заменяются на один оператор. Выражая Хп, р.п по известным формулам через Ёп,0>пи учитывая, что
= ип = const, где %cp)(t) = Ёп [p(t) - fair (t - T)cp(T)dx\ (4)
здесь Ёп- мгновенный модуль упругости, а 1п- ядро релаксации.
Учитывая (1), функцией времени в равенстве (4) будет p(t) = exp(-iMt) с медленно меняющейся амплитудой. Предполагая малость интеграла f™ R(r)dT, с помощью метода замораживания заменим соотношение (4) приближенным:
Ёп = Ёп[1 - Гс(шв) - is(coR)]cp, где
fiW^CnZV
о
Это позволяет исключить из вариационного уравнения интегральные члены и, в конечном итоге, время. В символическом виде его можно представить в виде
SAа + 8Аа + 8Ат = 0,
(2)
SG(u0j(X),v2) = 0.
(5)
где 5А„, 8Аа, ЗАт- виртуальные работы внутренних сил тел пружин, а также сил инерции с учетом сосредоточенных масс. Эти работы можно представить следующими соотношениями:
5Аа = -
п 8£^У,8Аа
п=1
N Ln
11
N Ln
п п п п
аГбЕГ-Ц^бъ,
п=1 ¡/=1
SA.
т = -Ж=1Рп fVnUrn (x,t^u^yy- -
Если п-я пластина, 1-я пружена и /-я опора вязко-упругие, то Ё>п, С1п, С/п представляются следующими формулами:
Ёп = Ёп1п(шл),С1п = С1п/ьпС./п = С1/п/1/пМ
где
f(^R) = 1- rc(^R) - iTs((0R)
- комплексная функция, числовые коэффициенты которой зависят от параметров ядра релаксации соответствующих вязкоупругих элементов,
D =
Evh-
(l-v%) , С1п Cl/'
п=1 =1
12
обобщенные мгновенные жесткости соответственно п-й пластины, 1-го амортизатора, 1-й опоры. В упругом случае Т)п = С1п = С1п, С1/п = С1/п, где Бп, С1П, С^п, - обобщенные жесткости соответственно п-й пластины, 1-го пружины, 1/-й опоры.
Необходимо найти спектр комплексных собственных частот
шк = + ш*,
где Шд - частоты, а - коэффициенты демпфирования собственных затуханий колебания.
Численные результаты. Рассмотрим конструкцию, представляющую собой пакет из двух параллельных квадратных упругих пластин с амортизатором и присоединенной массой. Ядро релаксации для амортизатора выбрано в виде
И^) = Аехр(-рс)са-1,
где А,р,а - параметры ядра [2].
Вязкость амортизатора принята такой, чтобы его деформация ползучести при квазистатическом процессе составляла малую долю от общей (~12%). Для этого случая параметры ядра следующие: А = 0.01, а = 0.1,/3 = 0.05 [2].
В отличие от упругой задачи, здесь исследовалась зависимость двух низких частот и соответствующих коэффициентов демпфирования от величины мгновенной жесткости амортизатора. Последняя изменялась от 10-4 до 10-1. Справа этот диапазон ограничен величиной, т.к. при С=С2 происходит смена второй формы. На рис. 1 показана зависимость первых двух частот и соответствующих коэффициентов
демпфирования от величины мгновенной
жесткости амортизатора С. Из анализа графиков следует, что диссипативные свойства этой системы в целом определяются не только реологией ее элементов, но существенно зависят от взаимодействия колебаний собственных форм.
Рисунок 1. Зависимость частот и коэффициентов демпфирования
Данный эффект выражается в том, что при некоторых условиях (о них ниже) и до некоторого значения жесткости амортизатора энергически более емкая (в данном случае вторая) форма диссипирует энергии меньше, чем форма менее энергоемкая. Затем, начиная с некоторого значения мгновенной жесткости амортизатора (в данном случае с С* = 5.4 • 10-3), процесс диссипации энергии собственными формами нормализуется и протекает согласно энергетической иерархии форм.
Практический вывод следующий: демпфирующие способности конструкции в основном определяет минимальный по абсолютной величине коэффициент демпфирования (в этом случае последними затухают колебания именно этой формы); глобальным (определяющим) коэффициентом демпфирования системы
1 2 является сначала ш1 до точки пересечения, а затем ш2.
Оптимальный в смысле затухания режим колебаний конструкции будет при С=С*, когда этот глобальный коэффициент демпфирования максимален.
^иток литературы:
1. Каюмов С.С., Сафаров И.И. Распространение и дифракция волн в диссипативно - неоднородных цилиндрических деформируемых механических системах. Ташкент: Фан, 2004 г. -250 с.
2. Колтунов М.А. Ползучесть и релаксация. - М.: Высшая школа, 1976. - 276 с.
