Рассеяние энергии при колебаниях неоднородных композитных структур. 3. Численный эксперимент Текст научной статьи по специальности «Физика»

УДК 534.222 Вестник СПбГУ. Математика. Механика. Астрономия. 2019. Т. 6 (64). Вып. 1 MSC 74E30

Рассеяние энергии при колебаниях неоднородных композитных структур. 3. Численный эксперимент

Л. В. Паршина1, В. М. Рябов2, Б. А. Ярцев1'2

1 Крыловский государственный научный центр,

Российская Федерация, 196158, Санкт-Петербург, Московское шоссе, 44

2 Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

Для цитирования: Паршина Л. В., Рябов В. М., Ярцев Б. А. Рассеяние энергии при колебаниях неоднородных композитных структур. 3. Численный эксперимент // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2019. Т. 6(64). Вып. 1. С. 144-156. https://doi.org/10.21638/11701/spbu01.2019.111

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

Ключевые слова: анизотропия, колебания, собственная частота, коэффициент механических потерь.

1. Влияние относительной толщины «жесткого» слоя изотропного вязкоупругого полимера на величины собственных частот f и коэффициентов механических потерь r/i низших тонов колебаний (г = 1,5) изучалось на примере двухслойной безопорной прямоугольной пластины (a = 0,16 м, b = a/2 = 0, 08 м). Силовой слой двухслойной пластины толщиной hi = a/100 = 0, 0016 м = const представляет собой квазиоднородную структуру, образованную совокупностью конечного числа армирующих слоев углепластика HMS/DX-209. Ориентация всех армирующих слоев

(¡5 Санкт-Петербургский государственный университет, 2019

силового слоя относительно связанной с направлениями осей пластины глобальной системы координат хуг одинакова. Предполагалось, что в рассматриваемом диапазоне температур упруго-диссипативные характеристики однонаправленного углепластика HMS/DX-209 частотно независимы. На одну из наружных поверхностей силового слоя нанесен слой частотно зависимого «жесткого» изотропного вязкоупру-гого полимера. Величины упруго-диссипативных характеристик всех компонуемых в слоистую структуру материалов приведены в [1].

Численные исследования проводились путем варьирования относительной толщины слоя «жесткого» изотропного вязкоупругого полимера Н € [0, 3] при фиксированных значениях углов ориентации армирующих слоев силового слоя пластины в = 0° и в = 15°. Полученные результаты представлены на рис. 1 и 2 в виде графиков зависимостей собственных частот ]г = fi(h) и коэффициентов механических потерь Пг = Пг(Н) первых пяти тонов колебаний рассматриваемой двухслойной безопорной прямоугольной пластины для трех значений температуры окружающей среды: Тс = 0° С, 20° С, 40° а

Прежде чем перейти к анализу полученных результатов, отметим, что в случае ориентации армирующих слоев силового слоя пластины относительно связанной с направлениями осей пластины глобальной системы координат хуг в направлениях в = 0°, 90° элементы матрицы жесткости Л^ = Л26 = Л45 = В16 = В26 = ^16 = ^26 = 0. Тогда описывающая связанные затухающие колебания пластины система дифференциальных уравнений движения (см. [2]) упрощается и записывается в виде

Л1 1 ио,хх + Л66и0,уу + (Л12 + Л66) го,ху + В1 1а0,хх +

+ В66ао,уу + (В12 + В66) во,ху + и2 (1оио + Дао) = 0 , (Л12 + Л66) ио,ху + Л66«о,хх + Л22ъо,уу + (В12 + В66) ао,ху+

+ В66во,хх + В22во,уу + и2 (/го + /во) = 0 , Л44 (во,у + ыо,уу) + Л55 (ао,х + ыо,хх) + и2/ЧЩ = 0 , (1)

ВцПо,хх + В66«о,уу + (В12 + В66) Уо,ху + ^11«о,хх +

+ ^66«о,уу + (^12 + ^66) во,ху - Л55 (ао + ^о,х) + и2 (/«о + /ао) = 0 , (В12 + В66) ио,ху + В66«о,хх + В22«о,уу + (^12 + ^66) ао,ху +

+ В66во,хх + ^22во,уу - Л44 (во + ^о,у) + и2 (/го + /во) = 0 ,

где и — круговая собственная частота колебаний.

Если армирующие слои силового слоя пластины произвольно ориентированы относительно осей пластины (в = 0°, 90°), то система дифференциальных уравнений движения такова:

Л1 1 ио,хх + 2Л16ио,ху + Л66ио,уу + Л16го,хх + (Л12 + Л66) го,ху +

+ Л26 го,уу + В1 1ао,хх + 2В16ао,ху + В66ао,уу + В16 вз,хх +

+ (В12 + В66) во,ху + В26во,уу + и2 (1оио + Дао) = 0 , Л16ио,хх + (Л12 + Л66) ио,ху + Л26ио,уу + Л66го,хх + 2Л26«о,ху +

+ Л22го,уу + В16ао,хх + (В12 + В66) ао,ху + В26ао,уу + В66во,хх+

+ 2В26во,ху + В22 во,уу + и2 (/ого + А во) = 0 , Л44 (во,у + Ыо,уу) + Л45 (ао,у + во,х + 2юо,ху) + Л55 (ао,х + ^о,хх) + и2/о^о = 0 ,

,xx + (B12 + B6e) v0,xy + + B26v0,yy + D 11a0,xx + 2Di6ao,xy + D66ao,yy + D^ß0,xx + + (D12 + D66) ß0,xy + D26ß0,yy - A45 (ß0 + W0,y) -

- A55 (a0 + W0,x) + w2 (I1U0 + ^0) =0 , (2) B16U0,xx + (B12 + B66) U0,xy + B26U0,yy + B66V0,xx +

+ 2B26V0,xy + B22V0,yy + D16a0,xx + (D12 + D66) a0,xy + D26a0,yy+ + D66ß0,xx + 2D26ß0,xy + D22ß0,yy - A44 ß + W0,y) —

- A45 (a0 + W0,x) + w2 (/1V0 + hß0) = 0 .

Анализ системы (1) позволяет заключить, что в случае двухслойной неоднородной пластины, силовой слой которой представляет собой ортотропную или трансверсально-изотропную структуру, возникает лишь мембранно-изгибное взаимодействие.

Система дифференциальных уравнений движения (2) описывает мембранно-изгибно-сдвиговые связанные затухающие колебания анизотропной слоистой композитной пластины. Связанность колебаний, обусловленная как несимметрией по толщине структуры, так и произвольностью углов ориентации армирующих слоев силового слоя пластины относительно глобальной системы координат xyz, порождает взаимные трансформации собственных форм.

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

Судя по приведенным на рис. 1 и 2 данным, согласующимся с результатами работы [3], каждой моде собственных колебаний двухслойной пластины соответствует эффективная относительная толщина слоя «жесткого» изотропного вязкоупругого полимера he, отвечающая максимуму функции щ = Vi (h). В зависимости от температуры окружающей среды Тс дальнейшее увеличение h двухслойной пластины с ортотропным (трансверсально-изотропным) силовым слоем может приводить как к возрастанию, так и к убыванию величин собственных частот fi без значимого изменения величин щ (рис. 1). Если же силовой слой двухслойной пластины представляет собой моноклинную структуру, то для некоторых мод колебаний при h > he происходит быстрое убывание функций щ = ni(h) (рис. 2, кривые 2 для Тс = 0° C, 20° C, 40° C и кривая 5 для Тс = 400C на графиках зависимостей щ = ni(h)). При этом в зависимости от Тс происходит как возрастание, так и убывание функций fi = fi (h). Отмеченное различие характера изменения функций щ = ni (h) объясняется тем, что в двухслойной пластине, силовой слой которой либо ортотропен, либо трансверсально-изотропен, возникает лишь мембранно-изгибное взаимодействие, не сопровождающееся изменением собственных форм колебаний по мере увеличения относительной толщины слоя «жесткого» изотропного вязкоупругого полимера на рассматриваемом отрезке изменения h G [0, 3].

Если же структура армирования силового слоя моноклинная, то в двухслойной пластине на отрезке h G [0, 3] возникает мембранно-изгибно-сдвиговое взаимодействие, приводящее к изменению собственных форм колебаний (рис. 3 демонстрирует изменения собственных форм второй моды колебаний).

2. Влияние углов ориентации армирующих слоев силового слоя пластины в на величины собственных частот fi и коэффициентов механических потерь щ низших

Рис. 1. Зависимости собственных частот и коэффициентов механических потерь низших тонов колебаний двухслойной безопорной прямоугольной пластины от относительной толщины слоя вязкоупругого материала (в = 0°).

Рис. 2. Зависимости собственных частот и коэффициентов механических потерь низших тонов колебаний двухслойной безопорной прямоугольной пластины от относительной толщины слоя вязкоупругого материала (в = 15°).

Рис. 3. Зависимость собственных форм низших тонов колебаний двухслойной безопорной прямоугольной пластины от относительной толщины слоя вязкоупругого материала (в = 15°, Т = 20° С).

тонов колебаний (г = 1,7) изучалось на примере двухслойной безопорной прямоугольной пластины (a = 0,16 м, b = a/2 = 0, 08 м, hi = h2 = a/100 = 0,0016 м, h =1), описанной в п. 1.

Вначале рассматривалась квазиоднородная прямоугольная пластина, образованная совокупностью одинаково ориентированных армирующих слоев углепластика HMS/DX-209 при отсутствии слоя «жесткого» изотропного вязкоупругого полимера. Такая пластина симметрична относительно срединной плоскости и, следовательно, Bim = 0 (l, m =1, 2, 6), Ii = 0. Тогда описывающая затухающие колебания квазиоднородной анизотропной пластины связанная система пяти дифференциальных уравнений движения второго порядка (2) упрощается и разделяется на две независимые системы, первая из которых характеризует связанные мембранные колебания:

Ai 1 u0,xx + 2A16u0,xy +A66u0,yy+A16V0,xx +

+ (A12+A66 ) v0,xy +A26v0,yy I0u0 — 0 , (3)

Ai6U0,xx+(Ai2+A66) U0,xy +A26U0,yy+A66V0,xx+2A26V0,xy +A22V0,yy +W2I0V0 = 0 .

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

A44 (@0,y + W0,yy) + A45 («0,y + fi0,x + 2w0,xy) + A55 (a0,x + w0,xx) + W2I0W0 = 0 , ^iia0,xx + 2Di6a0,xy + D66a0,yy + Di6^0,xx + (D12 + D66) @0,xy+

+ D26@0,yy - A45 (в0 + W0,y) - A55 (a0 + W0,x) + w2I2a0 = 0 , (4) Di6a0,xx + (D12 + D66) a0,xy + D26a0,yy + D66@0,xx + 2D26@0,xy +

+ D22@0,yy - A44 (в0 + W0,y) - A45 (a0 + W0,x) + w2I2^0 = 0 .

\\ 9 Г 7 \ 11 V \\

\1 а/ у!

д\ / /

ч ■■•«С* £ \з

Гч\Т /

'я;----- у*

ЦТ V

"г-- чи. -С

и

ТУ

3 2

5 /\ 6

4 у

75 е..

Рис. 4. Зависимости (а) собственных частот и (б) коэффициентов механических потерь низших тонов колебаний квазиоднородной пластины от угла ориентации армирующих слоев углепластика HMS/DX-209.

Неклассический вид графиков зависимостей собственных частот / = ^(в) и коэффициентов механических потерь гц = Цг(в) от угла ориентации слоев однонаправленного углепластика в € [0°, 90°] (г = 1, 7), представленных на рис. 4, свидетельствует о фейерверке эффектов, порождаемых связанностью мод колебаний пластины. Анализ полученных результатов позволил установить, что в зависимости от значения угла в (в € [0°, 90°]) каждой моде собственных колебаний соответствуют несколько различных собственных форм, порождаемых чередой взаимных трансформаций собственных форм взаимодействующих мод колебаний.

Для демонстрации описанных эффектов рассмотрим последовательность взаимных трансформаций взаимодействующих первой, второй и четвертой связанных мод колебаний, графики зависимостей собственных частот / = }г(в) и коэффициентов механических потерь гц = Цг(в) (г =1, 2, 4) которых приведены на рис. 4 (кривые 1, 2, 4). На рис. 4, а также нанесены графики зависимостей парциальных собственных частот / = $г(в) (кривые 1, 2, 4). Расчет значений парциальных собственных частот производился при выполнении условия Ql6 = Q26 = ^45 =0 в комплексной матрице жесткости слоя углепластика HMS/DX-209 [1]. Величины парциальных собственных частот определялись по «свободным» комплексным модулям упругости Ехх = Ехх(в) и «свободным» комплексным модулям межслойного сдвига Охг = Охг(в), Оуг = Суг(в), в то время как «свободный» комплексный модуль сдвига в плоскости армирования Оху = Оху (в) заменялся «чистым» комплексным модулем сдвига в плоскости армирования Охе(в) = Оху(в)/[1 — цху,х(в) • цх,ху(в)], соответствующим полному «запрещению» изгибных деформаций при сдвиге пластины в плоскости армирования. Здесь цхху(в), цху,х(в) —коэффициенты взаимного влияния первого и второго рода слоя ПКМ в локальной системе координат. Понятия «свободных» и «чистых» модулей введены Фойгтом [4]. Изменения собственных форм низших тонов связанных колебаний квазиоднородной пластины в зависимости от угла ориентации армирующих слоев углепластика HMS/DX-209 приведены на рис. 5. Собственные формы парциальных мод колебаний, не зависящие от ориентации армирующих слоев, представлены на рис. 6.

о = о0 в = 15° 0 = 30° в = 45° 9 = 60° в = 75° в = 90°

Рис. 5. Зависимость собственных форм связанных мод колебаний квазиоднородной пластины от угла ориентации армирующих слоев углепластика HMS/DX-209.

г = 1 г = 2 г = 3 г = 4 г = 5 г = 6

г = 7 г = 8 г = 9 г = 10 г = 11

1111111 И! ршш

1111111 ШИШ/Шё 11111 НИ наш

Рис. 6. Собственные формы парциальных мод колебаний квазиоднородной пластины из углепластика HMS/DX-209.

Совместный анализ приведенных на рис. 4-6 зависимостей позволяет утверждать, что на отрезке в € [0°, 30°] происходят взаимные трансформации собственных форм первой, второй и четвертой мод связанных колебаний. Об этом свидетельствует появление точек пересечения графиков зависимостей их парциальных собственных частот от угла ориентации армирующих слоев /- = /-¿(в) (кривые 1, 2, 4 на рис. 4, а) внутри областей, ограничиваемых функциями /т = /т(в) и /п = /п(в) связанных мод колебаний (кривые 1, 2, 4 на рис. 4, а). Пересечение графиков функций /2 = /2(в) и /4 = /(в) в точке в « 8° свидетельствует о том, что на примыкающем к ней отрезке с нечеткими границами по оси в начинается взаимная трансформация собственных форм второй и четвертой связанных мод колебаний. Поэтому при в > 8° собственная форма второй связанной моды трансформируется в собственную форму четвертой связанной моды, а собственная форма четвертой связанной моды — в собственную форму второй связанной моды. На отрезке с условным центром в точке в « 14° процесс взаимной трансформации собственных форм второй и четвертой мод дополняется процессом взаимной трансформации собственных форм первой и второй мод связанных колебаний. Это подтверждается пересечением графиков функций / 1 = /1 (в) и /4 = /4(в). Следовательно, при в > 14° собственная форма первой связанной моды трансформируется в собственную форму второй связанной моды, а собственная форма второй связанной моды трансформируется

в собственную форму первой связанной моды. Процесс взаимных трансформаций окончательно завершается при в = 90°. Поэтому собственные формы первой и четвертой мод связанных колебаний при в = 0° тождественны собственным формам второй и первой мод связанных колебаний при в = 90°. Следствием множества взаимных трансформаций на отрезке в € [0°, 90°] является тождественность собственной формы второй моды связанных колебаний при в = 0° собственной форме четырнадцатой моды связанных колебаний при в = 90°.

Аналогичная картина наблюдается и для более высоких тонов колебаний. Например, пересечение графиков функций /3 = /з(в) и /5 = /5(в) в точке в « 8° свидетельствует о том, что на примыкающем к ней отрезке изменения в происходит взаимная трансформация собственных форм третьей и пятой связанных мод колебаний. При в > 8° собственная форма третьей связанной моды трансформируется в собственную форму пятой связанной моды, а собственная форма пятой связанной моды — в собственную форму третьей связанной моды. В свою очередь, по мере увеличения угла ориентации армирующих слоев при в « 16° и в « 18° графики парциальных частот /3 = /3(в) и /5 = /5(в) пересекаются графиком парциальной частоты /10 = /10(в), что также свидетельствует о появлении областей взаимных трансформаций третьей и пятой связанных мод колебаний. Результатом этих трансформаций является тождественность собственных форм пятой моды связанных колебаний при в = 0° и в = 90°. Собственная форма третьей связанной моды колебаний при в = 90°, тождественная собственной форме десятой связанной моды колебаний при в = 0°, приведена на рис. 5.

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

3. Перейдем теперь к исследованию влияния угла ориентации армирующих слоев силового слоя в на величины собственных частот и коэффициентов механических потерь низших тонов колебаний двухслойной безопорной пластины. Поскольку в € [0°, 90°], то двухслойная пластина совершает связанные мембранно-изгибно-сдвиговые колебания, описываемые системой дифференциальных уравнений (2). Численные исследования проводились путем варьирования угла в при фиксированном значении относительной толщины слоя вязкоупругого полимера Н = 1.

Результаты исследования представлены на рис. 7, 8 в виде графиков зависимостей собственных частот /- = /¿(в), коэффициентов механических потерь щ = П-(в) и собственных форм первых семи тонов колебаний рассматриваемой пластины для трех значений температуры окружающей среды: Тс = 0° С, 20° С, 40° С. Приведенные на рис. 7, 8 данные качественно согласуются с полученными ранее результатами. Свидетельством этого является множество взаимных трансформаций собственных форм связанных мод колебаний безопорной двухслойной прямоугольной пластины. Последовательность взаимных трансформаций собственных форм помимо ориентации армирующих слоев силового слоя пластины также определяется частотой колебаний и температурой окружающей среды, поскольку упруго-диссипативные характеристики слоя «жесткого» изотропного вязкоупругого полимера демонстрируют существенную зависимость от этих параметров. Из сопоставления результатов, представленных на рис. 7, 8, следует, что при Тс =0° Си Тс = 20° С для всех в € [0°, 90°] последовательность собственных форм первых семи тонов связанных

Рис. 7. Зависимости собственных частот и коэффициентов механических потерь низших тонов колебаний двухслойной безопорной прямоугольной пластины от угла ориентации армирующих слоев углепластика HMS/DX-209.

Рис. 8. Зависимость собственных форм низших тонов колебаний двухслойной безопорной прямоугольной пластины от угла ориентации армирующих слоев углепластика HMS/DX-209 (Н = 1,

Тс = 0° С, Тс = 20° С.

колебаний двухслойной пластины совпадает. В то же время слой «жесткого» изотропного вязкоупругого полимера приводит к отличию этой последовательности от последовательности собственных форм квазиоднородной пластины из углепластика HMS/DX-209. В частности, в отличие от квазиоднородной пластины, в процесс трансформации собственных форм низших тонов колебаний которой вовлечены четыре связанные моды — первая, вторая и четвертая (рис. 4, 5), в двухслойной пластине в аналогичном процессе участвуют лишь две связанные моды (рис. 8) —первая и вторая. Взаимодействие первой и второй связанных мод колебаний приводит к взаимной трансформации их собственных форм. Результатом этой трансформации является тождественность собственных форм первой и второй мод связанных колебаний при в = 0° собственным формам второй и первой мод связанных колебаний при в = 90°.

Аналогичная картина наблюдается и для более высоких тонов колебаний.

Дальнейшее повышение температуры окружающей среды до Тс = 40° С сопровождается очередными изменениями последовательности взаимных трансформаций собственных форм колебаний.

Кстати, стоит отметить, что вследствие температурно-частотной зависимости вещественной части комплексного модуля упругости Ие Е(/, Тс) и коэффициента механических потерь п(/, Тс) «жесткого» вязкоупругого полимера повышение температуры окружающей среды на отрезке Тс € [0° С, 40° С] влечет за собой понижение собственных частот и повышение коэффициентов механических потерь всех рассматриваемых мод колебаний при в € [0°, 90°].

Литература

1. Паршина Л. В., Рябов В.М., Ярцев Б. А. Рассеяние энергии при колебаниях неоднородных композитных структур. 2. Метод решения // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2018. Т. 5(63). Вып. 4. С. 678—688.

2. Паршина Л. В., Рябов В. М., Ярцев Б. А. Рассеяние энергии при колебаниях неоднородных композитных структур. 1. Постановка задачи // Вестник Санкт-Петербургского университета. Ма-

тематика. Механика. Астрономия. 2018. Т. 5(63). Вып. 2. С. 300-309. https://doi.org/10.21638/11701 /spbu01.2018.210

3. Чернышев В. М. Демпфирование колебаний механических систем покрытиями из полимерных материалов. М.: Наука, 2004. 288 с.

4. Voigt W. Lehrbuch der Kristallphysik. Leipzig und Berlin: Teubner, 1928. 978 s.

Статья поступила в редакцию 20 ноября 2017 г.;

после доработки 29 декабря 2017 г.; рекомендована в печать 27 сентября 2018 г.

Контактная информация:

Паршина Людмила Валентиновна — канд. техн. наук, ст. науч. сотр.; [email protected]

Рябов Виктор Михайлович — д-р физ.-мат. наук, проф.; [email protected]

Ярцев Борис Александрович — д-р техн. наук, начальник сектора; [email protected]

Energy dissipation during vibrations of non-uniform composite structures.

3. Numerical experiments

L. V. Parshina1, V. M. Ryabov2, B. A. Yartsev1'2

1 Krylov State Research Center, Moskovskoye shosse, 44, St. Petersburg, 196158, Russian Federation

2 St. Petersburg State University, Universitetskaya nab., 7—9, St. Petersburg, 199034, Russian Federation

For citation: Parshina L.V., Ryabov V. M., Yartsev B. A. Energy dissipation during vibrations of non-uniform composite structures. 3. Numerical experiments. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2019, vol. 6(64), issue 1, pp. 144-156.

https://doi.org/10.21638/11701/spbu01.2019.111 (In Russian)

The influence of the relative thickness of the "stiff" viscoelastic layer polymer and the orientation of the reinforcing layers of the power layer on the values of their natural frequencies and loss factors of damped oscillations unsupported two-layered composite plate is investigated. It has been established that each natural vibration mode of two-layered plate corresponds to a certain effective relative thickness of isotropic viscoelastic polymer layer. Further increase in relative thickness for plates with orthotropic or transversally isotropic bearing layer depending on ambient temperature may result in both growth and decrease of natural frequencies without any considerable change of their energy dissipation. For the plates with monoclinic bearing layer, increase in relative thickness of isotropic viscoelastic polymer layer is accompanied by both decrease and increase of loss factor for different vibration modes because these structure have membrane-bending-shearing interaction that results in the change of natural vibration shapes. It is shown that change in stacking sequence of reinforcing layers in the bearing layer with respect to global axes of two-layered plate brings about mutual transformation domains of coupled vibration modes. It has been established that mutual transformation of natural coupled vibration mode shapes of unsupported quasi-uniform rectangular plate appears in case if at least in one of the plate directions one of the natural shapes has an even number of quarter-wavelengths, whereas another natural shape has odd. It has been demonstrated that temperature-frequency relationship of elastic-dissipative characteristics of the "stiff" viscoelastic polymer has a considerable effect upon natural frequencies and loss factors of all the investigated vibration modes for unsupported rectangular two-layered plate. Keywords: anisotropy, vibrations, natural frequency, loss factor.

References

1. Parshina L.V., Ryabov V.M., Yartsev B.A., "Energy dissipation during vibrations of nonuniform composite structures. 2. Method of solution", Vestnik St. Petersburg University. Mathematics. Mechanics. Astronomy 5 (63), issue 4, 678-688 (2018). (In Russian)

2. Parshina L. V., Ryabov V. M., Yartsev B. A., "Energy dissipation during vibrations of nonuniform composite structures. 1. Formulation of problem", Vestnik St. Petersburg University. Mathematics. Mechanics. Astronomy 5(63), issue 2, 300-309 (2018). (In Russian) https://doi.org/10.21638 /11701/spbu01.2018.210

3. Chernyshev V.M., Vibration damping of mechanical systems by means of polymeric coatings (Nauka Publ., Moscow, 2004, 288 p.). (In Russian)

4. Voigt W., Lehrbuch der Kristallphysik (Teubner, Leipzig, Berlin, 1928, 978 s.).

Received: November 20, 2017 Revised: December 29, 2017 Accepted: September 27, 2018

Author's information:

Ludmila V. Parshina — [email protected]

Victor M. Ryabov — [email protected]

Boris A. Yartsev — [email protected]