СВЯЗАННЫЕ КОЛЕБАНИЯ ВЯЗКОУПРУГИХ ТРЕХСЛОЙНЫХ КОМПОЗИТНЫХ ПЛАСТИН. 2. ЧИСЛЕННЫЙ ЭКСПЕРИМЕНТ

Вестник СПбГУ. Математика. Механика. Астрономия. 2021. Т. 8 (66). Вып. 1 УДК 534.121.1:678.067 МБС 74Е30

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

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

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

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

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

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

Для цитирования: Рябов В. М., Ярцев Б. А., Паршина Л. В. Связанные колебания вязко-упругих трехслойных композитных пластин. 2. Численный эксперимент // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2021. Т. 8 (66). Вып. 1.

C. 88-100. https://doi.org/10.21638/spbu01.2021.108

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

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

1. Математическая модель и численный метод решения дифференциальных уравнений затухающих колебаний трехслойных пластин, образованных двумя жесткими моноклинными слоями и мягким средним изотропным слоем из вязкоупругого полимера [1] позволяют определять параметры динамического отклика таких структур — собственные частоты, собственные формы колебаний, коэффициенты механических потерь. Очевидно, что эти параметры зависят от относительной толщины мягкого слоя вязкоупругого полимера Н, структуры армирования жестких слоев и температуры окружающей среды Тс. Многообразие структур армирования жестких слоев трехслойной пластины не оставляет надежд на рассмотрение бесконечного

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

Рис. 1. Схема (а) симметричной и (б) асимметричной трехслойных пластин.

множества их возможных вариантов. Поэтому далее рассматриваются лишь две специальные структуры армирования жестких слоев, компонуемых из однонаправленных слоев углепластика HMS/DX-209 (h(i) = h(3) = h = а/100 = 0.0016 м = const): симметричная и асимметричная. Симметричная структура армирования жестких слоев 1 и 3 образована совокупностью однонаправленных слоев углепластика, ориентированных в направлениях 0^) = 0(3) = 0 (рис. 1, а). В асимметричной структуре армирования однонаправленные слои HMS/DX-209 ориентированы в направлениях 0(1) = — 0(3) = 0 (рис. 1, б). В качестве материала мягкого слоя используется изотропный вязкоупругий полимер ВДМ-2.

2. Вначале рассмотрим влияние ориентации армирования на величины собственных частот и коэффициентов механических потерь квазиоднородной пластины из углепластика HMS/DX-209 (а = 0.16 м, b = а/4 = 0.04 м, h = 0.0016 м). В этом случае r =1, B1m = 0, E = G = p[2] =0, I1 = 0 и связанная система дифференциальных уравнений колебаний трехслойных пластин [1] распадается на две системы, первая из которых описывает мембранно-сдвиговые колебания квазиоднородной моноклинной пластины

Aii + Ai2 v,xy + Aie(û,xy + v,xx ) + Ai6 U,xy+

+A26v,yy + Абб

Ai6 U ,xx + A26 3,ху + A66(U,xy + 3,xx ) + Ai2 U,xy + +A22v,yy + A26 (U,yy + v,xy ) + W21oV = 0.

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

A44 (в ,y + w,yy) + A45 (à,y + w,xy) + A45 (в ,x + W,xy ) + A55 (<ä,x + W,xx ) + w21o w = 0,

Dii + Di2 ß ,xy + Di6(à ,xy + /?,xx) + Di6 à,xy + D26 ß,yy + D66(à ,yy + ß ,xy —

— A45 (/3 + w,y ) — A55 (à + w,x ) + w2/2 à = 0,

Di6 à,xx + D26 ß,xy + D66 (à ,xy + /5,xx ) + Di2à,xy + D22 Àyy + D26 (à ,yy + ß,xy ) —

— A44 (/3 + w,y ) — A45 (à + w,x ) + w2 /2/3 = 0.

а б

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

Если квазиоднородная пластина ортотропна, то A16 = A26 = A45 = =

D26 = 0 и система дифференциальных уравнений (1) принимает вид

A11 + A12 + Аб6 + ) + ш2 /с u = 0, Аб6 (u,xy + ) + A12 + A22 + ш2 /с v = 0.

Также упрощается и система дифференциальных уравнений (2):

A44 (Ду + Wyy) + A55 (а+ w3,xx) + ш2/сw3 = 0, D11 a+ D12¡3,Ху + D66(а,уу + /3,xy) - A55(а + w3,x) + ш2/2а = 0, D66(а,xy + /3,xx) + D12а,xy + D22Дуу - A44(/3 + w3,y) + ш2/2/3 = 0.

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

При проведении численных исследований рассматривались семь низших тонов связанных изгибно-крутильных колебаний квазиоднородной пластины из однонаправленного углепластика HMS/DX-209 (рис. 2). Неклассический вид графиков зависимостей собственных частот /г = /г (в) и коэффициентов механических потерь Пг = П (в) от угла ориентации армирования однонаправленного углепластика в Е [0°, 90°] (г = 1, 7), представленных на рис. 2, свидетельствует о фейерверке эффектов, порождаемых связанностью мод колебаний пластины. Анализ полученных результатов позволил установить, что в зависимости от значения угла в каждой моде собственных колебаний соответствуют несколько различных собственных форм, порождаемых чередой взаимных трансформаций собственных форм связанных мод колебаний.

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

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

Для демонстрации описанных эффектов рассмотрим последовательность взаимных трансформаций первой и второй связанных мод колебаний, графики зависимостей собственных частот / = /¿(0) и коэффициентов механических потерь П = Пг(^) ( =1, 2) которых приведены на рис. 2 (кривые 1, 2). На рис. 2, а также нанесены графики зависимостей парциальных изгибных / = /"¿(0) и крутильных / = /¿(0) собственных частот (кривые 1, 1). Величины парциальных собственных частот определялись из решения систем дифференциальных уравнений (3), (4) по «свободным» комплексным модулям упругости Ехх = Ехх(0) и «свободным» комплексным модулям межслойного сдвига Ох2 = Ох2 (0), Оу2 = Оу2 (0), в то время как «свободный» комплексный модуль сдвига в плоскости армирования Оху = Оху (0) заменялся «чистым» комплексным модулем сдвига в плоскости армирования Оху(0) = Оху(0)/[1 - Мху,х(0) • Мх,ху(0)]. Здесь ^х,ху(0), Мху,х(0) — коэффициенты взаимного влияния первого и второго рода. Понятия «свободных» и «чистых» модулей введены Фойгтом [2]. Изменения собственных форм низших тонов связанных колебаний квазиоднородной пластины из углепластика ИМ8/БХ-209 в зависимости от угла ориентации армирования приведены на рис. 3.

Совместный анализ приведенных на рис. 2, 3 зависимостей позволяет утверждать, что на отрезке 0 € [0°, 30о] происходят взаимные трансформации собственных форм первой и второй мод связанных колебаний. Об этом, в частности, свидетельствует появление точки пересечения графиков зависимостей их парциальных изгибных / = /"¿(0) и крутильных / = /"¿(0) собственных частот от угла ориентации армирования 0 (кривые 11, 11 на рис. 2, а) внутри области, ограниченной функциями /1 = /1(0) и /2 = /2(0) связанных мод колебаний (кривые 1, 2 на рис. 2, а). Пере-

сечение графиков функций / = /¿(0) и / = /¿(в) в точке в « 8° свидетельствует о том, что на примыкающем к ней отрезке с нечеткими границами по оси в начинается взаимная трансформация собственных форм первой и второй связанных мод колебаний. Собственная форма первой связанной моды, идентифицируемая при в < 8° как первая квазикрутильная, при в > 8° трансформируется в первую ква-зиизгибную собственную форму. Собственная форма второй связанной моды, идентифицируемая при в < 8° как первая квазиизгибная, при в > 8° трансформируется в первую квазикрутильную собственную форму. Аналогичная картина наблюдается и для более высоких тонов колебаний. Например, пересечение графиков функций /2 = /2 (в) и /2 = /2(в) в точке в « 11с (кривые 2, 2 на рис. 2, а) свидетельствует о том, что на примыкающем к ней отрезке изменения в происходит взаимная трансформация собственных форм третьей и пятой связанных мод колебаний. При в > 11° собственная форма третьей связанной моды трансформируется в собственную форму пятой связанной моды, а собственная форма пятой связанной моды — в собственную форму третьей связанной моды. В свою очередь, по мере увеличения угла в функция парциальной изгибной частоты /5 = /5 (в) при в « 33° пересекается с функцией парциальной крутильной частоты /3 = /з(в) (кривые 2, 2 на рис. 2, а), что также свидетельствует о появлении области взаимной трансформаций четвертой связанной (при в G [0°, 33°] — первой квазиизгибной относительно оси x) и пятой связанной (при в G [0°, 33°] — третьей квазикрутильной относительно оси x) мод колебаний. Результатом этой трансформации является тождественность собственной формы пятой связанной моды при в = 90° собственной форме седьмой связанной моды при в = 0°. Собственная форма седьмой связанной моды в результате взаимной трансформации с более высокой связанной модой на отрезке с центром в точке в « 13° (кривые 2, 2 на рис. 2, а) при в = 90° вырождается в четвертую изгибную моду колебаний относительно оси x.

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

3. Перейдем к рассмотрению двух специальных видов трехслойных пластин — симметричной и асимметричной. Квазиоднородные жесткие слои этих пластин имеют одинаковую толщину h(1) = h(3) = h, одинаковые инерционные характеристики

Р(1) = Р(3) = Р, /С1) = /с3) = /с = ph, /11) = /13) = h =0, = /23) = /2 = ph3/12 и отличаются лишь ориентацией армирования углепластика HMS/DX-209. В симметричной пластине все слои углепластика ориентированы под углами в(1) = в(3) = в к направлениям глобальных осей x(1) и x(3) жестких слоев 1 и 3 (рис. 1, а), а в асимметричной — под углами в(1) = — в(3) = в (рис. 1, б) к направлениям тех же осей. В этом случае справедливы равенства u(1) = u(3) = u, v(1) = v(3) = v, w(1) = w(3) = w, а(1) = а(3) = а, в(1) = в(3) = в

Для симметричной трехслойной пластины элементы матриц мембранных, смешанных и изгибных жесткостей жестких слоев связаны соотношениями A(m =

Am = Aim, Bm = Bm = 0, Dm = Dm = Am, подстановка которых в связанную систему дифференциальных уравнений колебаний трехслойных пластин приводит к двум независимым системам дифференциальных уравнений. Первая из этих систем описывает мембранно-сдвиговые колебания глобально моноклинной трехслой-

ной пластины:

А11и,хх + ,ху + А16(и,ху + V ,хх) + А16и,ху + А26 3 ,уу + А66 (и,уу + 3 ,ху ) +

+ ш2 /с +

2

и = 0,

А16Й,хх + А26«,ху + Аз6(Й,ху + Й,хх) + А12Й,ху + А22 У,уу + А26 (И,уу + £,ху ) +

уу ,ху

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

А44 (Ду + й,уу ) + А45(а,у + й,ху ) + А4б(/?,х + й,ху ) + А55 (а,х + й,хх) +

+

а

^[2]

+ ш2 /с +

2

й = 0,

^11й,хх + ^12/?,ху + ^16(а,ху + Дхх) + ^16«,ху + ^26/?,уу + ^66(а,уу + /3,ху)-

- А45(/3 + №„) - А55(а + гйж) - А ( т—а - ) + I /2 +

2 \Н

[2]

Р[ 2]Ь2Ь[2] 24

а = 0,

^16«,хх + ^26/3,ху + ^66(а,ху + Дхх) + ^12а,ху + ^22 Дуу + ^26(а,уу + /?,ху)-

2 и[2]

- А44(у8+ «;,„) -^45(а Н-йх) - ^ ( --/3 - йу ) + со2 [ 12 + ' ) /3 = 0.

24

Если жесткие слои симметричной трехслойной пластины ортотропны (А16 = А26 = А45 = ^16 = ^26 = 0), то системы дифференциальных уравнений (5) и (6) упрощаются:

А11и,хх + А12^,ху + А66(и,уу + ^,ху ) + ш2 /с + А66(и,ху + £,хх) + А12«,ху + А22£,уу + ш2 ( /с +

р[ 2]^[2] 2

2

и = 0,

V = 0,

а

А44 (/3

+

+ ш2 /с +

2

й = 0,

+ 012Р,хУ + 066(алу + ¡3^ху) - А55(а + юх)--— & - ] +

.2 I т , Р[2]^[2]

+ ш2 /2 +

24

а = 0,

ал. / л.

£>66+ Джж) + -012«,+ £>22Дуу - ^44(/3 + гй у)--— ( —¡3 - ] +

+ ш2 /2 +

[2] 24

/3 = 0.

2

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

4. Для асимметричной трехслойной пластины элементы матриц мембранных,

смешанных и изгибных жесткостей жестких слоев связаны соотношениями аЦ =

А(3) = Ли А(1) = А(3) = А(1) = - А(3) = А(1) = А(3) = Иоо А(1) ^ А13) =

а11 = а11 , а12 = А12 = а12, а16 = а16 = А16, а22 = А22 = А22, а26 = А26 =

Л(1) = Л3) = а(1) = а(3) = а(1) = -Л(3) = а(1) = Л(3) = Акк

а26, а66 = а66 = а66, а44 = А44 = а44, а45 = А45 = а45 , а55 = а55 = а55 ,

в(1) = в(3) = 0 (/ т = 1 2 6) п(1) = п(3) = П п(1) = п(3) = п п(1) = п(3) =

= = 0 (/ , т = 1 2 , b), п11 = п11 = П1Ь п12 = п12 = П12, П16 = -П16 =

П П(1) = П(3) = П п(1) = п(3) = п п(1) = п(3) = п п16, п22 = п22 = п26 = п26 = П26, П66 = п66 = п66.

Подстановка этих соотношений в связанную систему дифференциальных уравнений колебаний трехслойных пластин приводит к ее разделению на две независимые системы. Первая система, описывающая мембранные и сдвиговые колебания асимметричной трехслойной пластины, совпадает с системой дифференциальных уравнений (7), описывающей мембранные и сдвиговые колебания глобально орто-тропной трехслойной пластины. Вторая система, описывающая изгибные и крутильные колебания асимметричной трехслойной пластины, совпадает с системой дифференциальных уравнений (8), описывающей изгибные и крутильные колебания глобально ортотропной трехслойной пластины.

Таким образом, связанные затухающие колебания симметричной трехслойной пластины описываются двумя системами дифференциальных уравнений (5), (6) структуры которых близки структурам систем соответствующих дифференциальных уравнений, описывающих затухающие колебания квазиоднородной моноклинной пластины (1), (2). В то же время связанные затухающие колебания асимметричной трехслойной пластины описываются двумя системами дифференциальных уравнений (7), (8), совпадающими с системами соответствующих дифференциальных уравнений, описывающих затухающие колебания глобально ортотропной трехслойной пластины.

5. Влияние ориентации армирования жестких слоев симметричной и асимметричной трехслойных пластин на величины их собственных частот и коэффициентов механических потерь исследовалось путем варьирования угла 0 при фиксированных значениях относительной толщины мягкого слоя вязкоупругого полимера Н = ^[2]/Н = 0.2 и температуры окружающей среды ТС = 30°.

Результаты численных исследований зависимостей собственных частот / = /¿(0), коэффициентов механических потерь п = ^¿(0) и собственных форм шести низших тонов колебаний симметричной трехслойной пластины приведены на рис. 4, 5. Сопоставление графиков зависимостей / = /¿(0) и собственных форм колебаний для глобально моноклинной симметричной трехслойной пластины с аналогичными зависимостями для моноклинной квазиоднородной пластины из однонаправленного углепластика ИМ8/БХ-209 (рис. 2, 3) позволяет отметить подобие их изменений, обусловленное изгибно-крутильной связанностью. Об этом свидетельствует появление множества областей взаимных трансформаций собственных форм связанных мод колебаний. Последовательность взаимных трансформаций собственных форм симметричной трехслойной пластины определяется ориентацией армирования ее жестких слоев, а также частотой колебаний и температурой окружающей среды, поскольку упруго-диссипативные характеристики мягкого слоя вязкоупру-гого полимера демонстрируют существенную зависимость от этих параметров.

О 15 30 45 60 75 6,° 0 15 30 45 60 75 8 а б

Рис. 4- Зависимости собственных частот (а) и коэффициентов механических потерь (б) низших тонов колебаний симметричной трехслойной пластины от угла ориентации армирования жестких слоев.

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

Рис. 5. Зависимость собственных форм связанных мод колебаний симметричной трехслойной пластины от угла ориентации армирования жестких слоев (Тс = 30°О).

Из сопоставления представленных на рис. 4, 5 и на рис. 2, 3 результатов следует, что для всех 0 Е [0°, 90°] последовательности взаимных трансформаций собственных форм трех низших тонов связанных колебаний симметричной трехслойной пластины совпадают с аналогичными последовательностями для связанных колебаний квазиоднородной пластины. Для иллюстрации сказанного на графике за-

висимостей собственных частот связанных колебаний симметричной трехслойной пластины / = /¿(0) дополнительно нанесены графики зависимостей парциальных изгибных / = /"¿(0) и крутильных / = /"¿(0) собственных частот (соответственно точечные кривые 1~4 и пунктирные кривые 1-3 на рис. 4, а). Величины парциальных собственных частот определялись из решения системы дифференциальных уравнений (8) по «свободным» комплексным модулям упругости = Ехх(0), «свободным» комплексным модулям межслойного сдвига = (0), = (0) и «чистому» комплексному модулю сдвига в плоскости армирования (0) =

«ху (0)/[1 (0) (0)].

Из совместного анализа приведенных на рис. 4, 5 зависимостей следует, что на отрезке 0 € [0°, 30°] происходят взаимные трансформации собственных форм первой и второй, третьей и четвертой, а также пятой и седьмой мод связанных колебаний. Об этом свидетельствует появление точек пересечения графиков зависимостей соответствующих парциальных изгибных f¿ = /"¿(0) и крутильных / = /"¿(0) собственных частот от угла ориентации армирующих слоев 0 внутри областей, ограниченных функциями /т = /т(0) и /п = /п(0) связанных мод колебаний (пересечения кривых 1, 1 при 0 « 10°, кривых 1 3 при 0 « 8°, кривых 1 3 при 0 « 12°).

6. Результаты численных исследований зависимостей собственных частот / = /¿(0), коэффициентов механических потерь п = ^¿(0) и собственных форм шести низших тонов колебаний асимметричной трехслойной пластины приведены на рис. 6, 7. В безопорной асимметричной трехслойной пластине возникает взаимодействие изгибных мод колебаний в двух взаимно ортогональных плоскостях (изгибно-изгибное взаимодействие). Для иллюстрации этого взаимодействия на графике зависимостей собственных частот колебаний / = /¿(0) нанесены графики зависимостей парциальных изгибных собственных частот / = /"¿(0) (точечные кривые 1~1 на рис. 6, а). Величины парциальных изгибных собственных частот определялись из решения системы дифференциальных уравнений (8) по «чистым» комплексным модулям упругости Ехх(0) = Ехх(0)/[1 — • ^х,ху(0)] и «свободным» комплексным модулям межслойного сдвига = (0), = (0).

Совместный анализ зависимостей, приведенных на рис. 6, 7, позволяет утверждать, что на отрезке 0 € [15°, 30°] происходит взаимная трансформация собственных форм шестой и седьмой мод колебаний. Об этом свидетельствует появление точки пересечения графиков зависимостей соответствующих парциальных собственных частот при 0 « 23° (кривые 3, 3). В результате взаимной трансформации собственная форма шестой связанной моды, при 0 = 0° идентифицируемая как первая изгибная относительно оси х, при 0 = 90° идентифицируется как третья изгибная относительно оси у (см. рис. 7). Представленные результаты позволяют заключить, что для безопорной асимметричной трехслойной пластины взаимодействие изгиб-ных мод колебаний в двух взаимно ортогональных плоскостях возникает в том случае, если в направлениях пластины х или у обе собственные формы характеризуются либо четным, либо нечетным числом четвертей волны.

7. Численные исследования влияния относительной толщины мягкого слоя изотропного вязкоупругого полимера Н на величины собственных частот / = /¿(Н) и коэффициентов механических потерь п = ^¿(Н) шести низших тонов колебаний (¿ = 1,6) проводились для фиксированных углов ориентации армирования жестких слоев трехслойных пластин. Рассматривались глобально моноклинная симметричная (0(1) = 0(3) = 0 = 15°) и глобально ортотропная асимметричная (0(1) = — 0(3) = 0 = 15°) трехслойные пластины. Полученные результаты приведены на рис. 8, 9 в

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

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

Рис. 7. Зависимость собственных форм связанных мод колебаний асимметричной трехслойной пластины от угла ориентации армирования жестких слоев (Тс = 30°С).

виде графиков зависимостей / = /(Н) и п = П (Н). Судя по приведенным на рис. 8, 9 данным, согласующимся с результатами работы [3], каждой моде собственных колебаний как симметричных, так и асимметричных трехслойных пластин соответствует эффективная относительная толщина мягкого слоя изотропного вязкоупругого полимера Не, при которой функция п = П (Н) достигает максимума.

0.7

Л/

0.5

0.3

0.1

6 ±1

^^ 3

2 А]

1/ Л

0.2

0.4 0.6

б

0.8

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

0.65

0.45

0.25

0.05

0

11

1

(4

6 Л

0.2

0.4

0.6

0.8

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

а

Дальнейшее увеличение величины Н > Ке не только не приводит к увеличению значений коэффициентов механических потерь пг = Пг ( Н), но часто сопровождается снижением уровней демпфирования. Последнее объясняется температурно-частотной зависимостью упруго-диссипативных характеристик вязкоупругого полимера, порождающей аналогичную зависимость величины Не = Не (/, Тс). При

h = he функции n = достигают глобального максимума. Эти же области

значений аргумента h соответствуют максимальным значениям функций д/¿(h)/dh. Сопоставление представленных на рис. 8, 9 функций n = n«(h) (h G [0, 1]) позволяет отметить, что изгибно-крутильная связанность колебаний симметричных трехслойных пластин позволяет реализовать уровни демпфирования, превышающие уровни демпфирования асимметричных трехслойных пластин.

Литература

1. Рябов В.М., Ярцев Б. А., Паршина Л. В. Связанные колебания вязкоупругих трехслойных композитных пластин. 1. Постановка задачи. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 7 (65), вып. 3, 469—480 (2020). https://doi.org/10.21638/spbu01.2020.309

2. Voigt W. Lehrbuch der Kristallphysik. Leipzig und Berlin, Teubner (1928).

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

Статья поступила в редакцию 2 февраля 2020 г.;

после доработки 29 февраля 2020 г.; рекомендована в печать 17 сентября 2020 г.

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

Рябов Виктор Михайлович — д-р физ.-мат. наук, проф.; v.ryabov@spbu.ru

Ярцев Борис Александрович — д-р техн. наук, нач. сектора; boris_yar@mail.ru

Паршина Людмила Валентиновна — канд. техн. наук, ст. науч. сотр.; ludaparshina@mail.ru

Coupled vibrations of viscoelastic three-layer composite plates.

2. Numerical experiments

V. M. Ryabov1, B. A. Yartsev1'2, L. V. Parshina2

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

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

For citation: Ryabov V. M., Yartsev B. A., Parshina L. V. Coupled vibrations of viscoelastic

three-layer composite plates. 2. Numerical experiments. Vestnik of Saint Petersburg University.

Mathematics. Mechanics. Astronomy, 2021, vol. 8(66), issue 1, pp. 88-100.

https://doi.org/10.21638/spbu01.2021.108

Results of numerical studies of the influence of reinforcement orientation of hard layers and the relative thickness of the soft layer of an isotropic viscoelastic polymer on the values of natural frequencies and coefficients of mechanical losses are discussed. In a symmetrical three-layer plate, a bending-torsion interaction arises, generating mutual transformations of the eigenmodes of the coupled vibration modes, if at least one of its own forms in plate directions is characterized by an even number quarters of a wave, and another proper form — an odd number of quarters the waves. In an unsymmetrical asymmetric three-layer plate, bending-bending interaction arises, generating mutual transformations eigenmodes of coupled vibration modes in two mutually orthogonal planes if in the main directions of the plate both of their own forms characterized by either an even or an odd number of quarters of a wave. It was found that each mode of natural vibrations of three-layer

plates corresponds to the effective relative thickness of the soft isotropic layer viscoelastic polymer.

Keywords : plate, composite, anisotropy, viscoelastic polymer, temperature-frequency dependence, coupled vibrations, natural frequency, loss factor.

References

1. Ryabov V.M., Yartsev B.A., Parshina L. V. Coupled vibrations of viscoelastic three-layer composite plates. 1. Formulation of Problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 7 (65), iss. 3, 469-480 (2020). https://doi.org/10.21638/spbu01.2020.309 (In Russian) [Engl. transl.: Vestnik St. Petersb. Univ. Math. 53, iss. 3, 320-328 (2020). https://doi.org/10.1134/S1063454120030127].

2. Voigt W. Lehrbuch der Kristallphysik. Leipzig und Berlin, Teubner (1928).

3. Parshina L. V., Ryabov V. M., Yartsev B. A. Energy dissipation during vibrations of heterogeneous composite structures: 3. Numerical Experiments. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 6 (65), iss. 1, 144-156 (2019). https://doi.org/10.21638/11701/spbu01.2019.111 (In Russian) [Engl. transl.: Vestnik St. Petersb. Univ. Math. 52, iss. 1, 102-111 (2019). https://doi.org/10.3103/S1063454119010114].

Received: February 2, 2020 Revised: February 29, 2020 Accepted: September 17, 2020

Authors' information:

Victor M. Ryabov — v.ryabov@spbu.ru Boris A. Yartsev — boris_yar@mail.ru Ludmila V. Parshina — ludaparshina@mail.ru