КОЛЕБАНИЯ КОНСОЛЬНО-ЗАЩЕМЛЕННОЙ ТОЛСТОЙ ПЛАСТИНЫ Текст научной статьи по специальности «Физика»

Папков С.О. Колебания консольно-защемленной толстой пластины // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 2. С. 106-117. БО!: 10.15593/регш.шесЬ/2021.2.10

Papkov S.O. Vibrations of a cantilevered thick plate. PNRPU Mechanics Bulletin, 2021, no. 2, pp. 106-117. DOI: 10.15593/perm.mech/2021.2.10

ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2021 PNRPU MECHANICS BULLETIN

https://ered.pstu.ru/index.php/mechanics/index

■H

БО1: 10.15593/регш.шесЬ/2021.2.10 УДК 539.3

КОЛЕБАНИЯ КОНСОЛЬНО-ЗАЩЕМЛЕННОЙ ТОЛСТОЙ ПЛАСТИНЫ С.О. Папков

Севастопольский государственный университет, Севастополь, Россия

О СТАТЬЕ

АННОТАЦИЯ

Получена: 20 октября 2020 г. Принята: 12 мая 2021 г. Опубликована: 12 июля 2021 г.

Ключевые слова:

консольно-защемленная пластина, теория Рейснера - Миндлина, колебания, собственные частоты, собственные формы, аналитическое решение, бесконечная система линейных уравнений.

В работе впервые представлено аналитическое решение задачи о свободных колебаниях консольно-защемленной толстой ортотропной пластины. Данная задача является достаточно громоздкой для использования точных методов теории упругости, в связи с чем для ее решения разрабатывались методы на основе вариационного подхода. В статье предлагается использовать метод суперпозиции для построения общего решения уравнений колебаний пластины в виде ряда по частным решениям, полученным на основе разделения переменных. При этом частные решения по одной из координат выбираются в форме тригонометрических функций специального вида (модифицированная тригонометрическая система). Построенное решение, в отличие известных в литературе решений на основе вариационного подхода, точно удовлетворяет уравнениям колебаний. Использование модифицированной тригонометрической системы функций позволяет получить единообразные выкладки для четных и нечетных форм колебаний и уменьшить число граничных условий на сторонах пластины с 12 до 9, при этом 5 из 9 граничных условий также выполняются точно. Структура представленного решения такова, что на границе пластины каждая из кинематических или силовых характеристик пластины представляется в виде суммы двух рядов - тригонометрического ряда и ряда по гиперболическим функциям. Относительно неопределенных коэффициентов рядов, представляющих решение, из оставшихся граничных условий, получена бесконечная система линейных алгебраических уравнений. Сходимость решения бесконечной системы по методу редукции исследуется численно. Приведены примеры численной реализации, на основе полученного решения выполнены численные исследования спектра собственных частот консольно-защемленной толстой пластины, как при вариации упругих характеристик материала, так и при вариации геометрических параметров.

©ПНИПУ

© Папков Станислав Олегович - д.ф.-м.н., доц., зав. каф., e-mail: stanislav.papkov@gmail.com, : 0000-0001-8624-0846.

Stanislav О. Papkov - Doctor of Physical and Mathematical Sciences, Associate Professor, Head of Department, e-mail: stanislav.papkov@gmail.com, : 0000-0001 -8624-0846.

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

VIBRATIONS OF A CANTILEVERED THICK PLATE

S.O. Papkov

Sevastopol State University, Sevastopol, Russian Federation

ABSTRACT

It has been for the first time that an analytical solution to the problem of free vibrations of a cantilevered thick orthotropic plate is presented. This problem is quite cumbersome for using the exact methods of the theory of elasticity; therefore, methods based on the variational approach were developed to solve it. The paper suggests using the superposition method to construct a general solution of the vibration equations of a plate in the series form of particular solutions obtained with the help of a variables separation. The particular solutions of one of the coordinates are built in the form of trigonometric functions of a special type (modified trigonometric system). The constructed solution, in contrast to the solutions known in the literature on the basis of the variational approach, accurately satisfies the equations of vibrations. The use of a modified trigonometric system of functions makes it possible to obtain uniform formulas for even and odd vibration shapes and to reduce the quantity of boundary conditions on the plate sides from twelve to nine ones, while five of the nine boundary conditions are also accurately satisfied. The structure of the presented solution on the plate boundary is such that, each of the kinematic or force characteristics of the plate is represented as a sum of two series, i.e. a trigonometric series and a series in hyperbolic functions. Remaining boundary conditions make it possible to obtain an infinite system of linear algebraic equations with respect to the unknown coefficients of the series representing the solution. The convergence of the solution by the reduction method of the infinite system is investigated numerically. Examples of the numerical implementation are given; numerical studies of the spectrum of natural frequencies of the cantilevered thick plate were carried out based on the obtained solution, both with varying elastic characteristics of the material and with varying geometric parameters.

©PNRPU

ARTICLE INFO

Received: 20 October 2020 Accepted: 12 May 2021 Published: 12 July 2021

Keywords:

cantilevered plate, Reisner - Mindlin theory, vibrations, natural frequencies, natural shape modes, analytical solution, infinite system of linear equations.

Введение

Проблема колебаний консольно-защемленных пластин относится к одной из наиболее давних задач механики деформированного твердого тела. Элементы в форме консольно-защемленных по одному из краев прямоугольных пластин можно встретить во многих технических и строительных конструкциях, они являются неотъемлемыми элементами зданий, авиационного и автомобильного транспорта, различных приборов. По этой причине первые работы, посвященные проблеме колебаний подобных пластин, в рамках теории тонких пластин относятся к середине ХХ в. Young [1], используя произведения балочных функций в качестве базиса метода Релея - Ритца, одним из первых получил достаточно точные верхние оценки собственных частот колебаний изотропной квадратной консольно-защемленной пластины. В работах [2, 3] для решения той же проблемы было предложено использовать метод, основанный на вариационном принципе Рейснера. Заметим, что при этом подходе граничные условия на свободных сторонах пластины выполнялись не полностью, в частности условия на сдвиговые силы игнорировались. В работах [4-6] развиты модификации метода Релея - Ритца для анализа колебаний прямоугольных изотропных пластин в случае консольного защемления, которые позволили с достаточной точностью найти первые собственные частоты пластины и построить собственные формы колебаний. Можно заметить, что во всех упомянутых работах используются методы, основанные на вариационном подходе.

Однако, как известно, в ряде прикладных задач теория тонких пластин зачастую не позволяет адекватно описать структурный элемент, что привело к развитию теорий толстых пластин, в частности теории Рейснера -Миндлина [7]. Одно из первых теоретических и экспериментальных исследований колебаний толстых консольно-защемленных пластин представлено в [8]. В работе Lim и др. [9] был развит метод Ритца для анализа колебаний толстых изотропных пластин для ряда классических граничных условий, в том числе для консольного защемления. В качестве базиса метода Ритца авторами выбирается тригонометрическая система функций, обеспечивающая как выполнение граничных условий, так и достаточно быструю сходимость метода. Метод Ритца и его модификации остаются основой для теоретического анализа консольно-защемленных пластин и на сегодняшний день. Среди недавних работ можно отметить работы Liew и др. [10], где получены собственные частоты колебаний на основе pb-2 метода Релея - Ритца, в работе Cheung и Zhou [11] в качестве базиса метода Релея -Ритца предлагается взять статические балочные функции. Вариационный DSC - element method был предложен [12] для анализа свободных колебаний толстых пластин при различных граничных условиях. В работе [13] строится решение задачи о колебаниях консольно-защемленной пластины на основе вариационного подхода применительно к трехмерным уравнениям теории упругости.

На основе аналитических решений [14] в работах [15-17] развит Dynamic Stiffness method для колебаний

изотропных толстых пластин. В статье [18] развит асимптотический подход к анализу высокочастотных колебаний прямоугольных пластин при различных граничных условиях. Для получения численного решения исследуемой задачи в [19, 20] используется метод конечных разностей, развитый с использованием теории степенных рядов. Численный подход на основе метода конечных разностей и метода конечных элементов также можно найти в работах [21-24]. В статье [25] авторы фактически используют метод суперпозиции для сведения задачи о колебаниях консольно-защемленной изотропной тонкой пластины к бесконечной системе линейных алгебраических уравнений.

В представленной статье на основе подхода [26-28] предлагается аналитическое решение задачи о колебаниях консольно-защемленной толстой ортотропной пластины. В отличие от предыдущих работ, посвященных колебаниям консольно-защемленной толстой пластины, построенное решение точно удовлетворяет разрешающим уравнениям колебаний. Выполнение граничных условий приводит к бесконечной системе алгебраических уравнений относительно неопределенных коэффициентов в общем решении.

1. Основные уравнения и структура общего решения задачи

где О —

4 _ phm 4 h 4

Д

О4 - частотные параметры;

и _ А 7, _ Д(4 7, _ ^55 7 _ Аб й _ Аб + А2

Кг. --• /С. --• /Сс --. /С,- --• /С --.

д д д д д

Упругие константы Д, Д, Д2 Д66, и А связаны с техническими константами (модулями Юнга Е и Е, модулями сдвига С12 и Оп, С23, и коэффициентами Пуассона у12 и у21) посредством соотношений

А =

к Е,

Д =-

к Е,

12(1 -У12У21) ; 12(1 -у12у21) к3у,, Е,

Д12 =

12(1 -У12У21)

Дб =

12

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

х = -а : Ш = фх = фу = 0, х — а : М — М — О — 0,

х XV ^х '

у — ±Ь : М — М — О — 0,

у ху

(3)

Рассмотрим свободные колебания однородной тол- где моменты и сдвиговые силы имеют вид

стой ортотропной пластины {(х, у) е[-а; а]х[-Ь; Ь]} постоянной толщины h. Согласно гипотезе теории толстых пластин Рейснера - Миндлина [7] элемент пластины смещается в прямолинейном направлении относительно срединной поверхности пластины, что приводит к следующему представлению компонентов вектора смещений:

и (х, ) — и0 (х, у, г) + 2Ф0 (х, у, г), V(х, у, 2, /) — V0 (х, у, () + 2ф°, (х, у, /),

Ш (х, у, 2, /) — V0 (х, у, г )

(1)

В случае гармонических колебаний углы поворота ф0 (х, у, г) — Ф (х, у)е,т', ф0 (х, у, г) — Ф (х, у)е'т' и прогиб пластины w0( х, у, г) — Ш (х, у)е,т' должны удовлетворять уравнениям колебаний толстой ортотропной пластины:

•52ФУ

дх2

д2

ду2

ду2

д2 ф дх2

ймЭу

у + к^-к, + фу | + О4фу — 0, (2)

дхду

ду

к ¡дШ+дФх|+к

41 дх2 дх 1 5

(^2

д 2Ш дф.

Л

+ О4Ш — 0,

мх дф!

.ар0, м

5Ф:,, а?;

Д

дх

ду д

&

Мху А

— К

ду дх

(4)

О — к д к4

дх

+ф0

о

О — К Д

0

-+ Ф0

ду ф

Четное решение относительно оси у (индекс ] =0) или же нечетное решение по у (=1) краевой задачи (2)-(3)

Фх

чф у J

— I

фхк

строится в форме суммы рядов с неопределенными коэффициентами Хп и Уы по модифицированной тригонометрической системе функций {Тк (ал4: х} и

{Т (PnjУ)\. Используя подход, представленный в статье [28], можно получить

Ш — 8к0 (Х1,0кА,0кНЛ (РшУ) + Х2,0к4,0кН (Р2,0ку)) + + 8Л (¥1,0кЕ1,0Нк (?1,0х + ¥2,0кЕ2,0Нк (?2,0х) +

+ZZ X мАмН (Р, My)Tk (ankx) +

n= 1 1= 1

w 3

+ ZZ ^nAnHk (qunjX)Tj (Pnj-y);

n=1 1=1

Фхк = 8j0 (Y1,0kH'; (qi,0jx) + Y2,0kHk (?2,0jX)) +

w 3

+ ZZ V, nHj (p, ,nky)Tk'(«nkX) +

n=1 1=1

w 3

+ ZZ , „Л (q, ,njX)Tj j

n=1 1=1

Фу,к = 8k0 (X1,0kH_ j (P1,0ky) + X2 ,0k j (P2,0ky)) +

w 3

+ZZ x nkH; (Pin nky)Tk (ankX)+

n=1 1=1

w 3

+ ZZYKnG,.njHk (jT/j (5)

n=1 1=1

где T(z) = (1-k)cosz + ksinz ; H(z) = (1-k)coshz +

+ ksinhz ; dmn - символ Кронекера; a^ =—[ n-К|;

a ^ 2 )

в —(»- 2 )•

Величины q; ., pímt (n =1, 2, ...) являются корнями

ветви l следующих характеристических бикубических уравнений:

С0 Р6 + С1Р4 + С2 Р2 + С3 = 0, (6)

daqe + dqA + ^q2 + d3 = 0, (7)

где

C0 = k2k5k6 .

q = a2 (k5 (k2-k¡)~k2 (k5 + k4k6)) + + k^ (k + kg) + ^^ k^k^ k^k^k^,

c2 = a4 {k4 ('k2 + k2 - P) + ksk6) + + a2 (q4 (k2 -k¡)-k5Q.¡ (l + k6)-k2 (q4 + k4Q.¡) +

+k4(2k,(k6+k)-k6nl)) +

+ Q4 (Q4 (k2 + k6) - k^k4) + k5 (Q4 (Q4 - k4) - k6Q4) c3 = -a6k4k6 + a4 (-kk + kQ4 + kQ¿ (1 + k )) +

+a2 (q4 ( к + kk -Q4)+ +Q4 (k4k5 - k¡Q4 - K4Q4 ))+Q4 (q4 - k) (Qh - k5);

d0 = К4К6 '

rf, = P2 [k4 (k2 -k¡-k2)-k5k6) + D.X (1 + k6).+ QX -k4k5

с12 = р4 (к5 [к2 +к26-к2 ^ + к2к4к6) + + Р2 ¡2к4к5 (к6 +к)-аА (к2 +к26 -к2)-(к2к4 + к5 + кАк6 + к5к6)) + +о4 (п4 (1+к6)- к - кА к6)+к404 (04 - К),

= -в6к2к5к6 + в4 (-к2к4к5 + к2к6О4 + к04 (к + к)) +

+в2 (о4 (к2 к4 + к5к6 - к204)+ +04 (к4к5 - к6о4 - к5о4 ))+о4 (04 - к4) (04 - к5).

Для случая п = 0 уравнения (6), (7) вырождаются [28] в биквадратные уравнения вида

к2к5 р04 + (к2о4 + к5о4 4 р02 +о4 (04 - к5) = 0, (8)

к4^0 +(04 + к4П4 ) ^ +04 (04 - кА ) = 0. (9)

Заметим, что все представленные характеристические уравнения (6)-(9) допускают точное решение.

При этом значения величин А, м , Е, п , В, м и О, п, входящих в решение (5), могут быть выражены следующим образом для каждой ветви l характеристических уравнений:

Щр1„к - к4 {кгр21пк - к6а2пк -к5+01) Pi.uk (кк5р1„к ~кк4а2пк +к4к5 +Ш4) к\р\л + {к2Р1пк - К°4 - К + ) {к5Р1пк ~ Ка1к + )

апкР1,пк {кк5Р1пк ~кк4а1к +Кк5 +ка4)

кКч1т - к5 (д1„1 - КР1 -к4+П4 )

{кКч1т-ЩК +к4к5 +т4)

A1,nk = B1 ,nk =

E1 ,nj

Gl,nj =

kWi,„j +{qt„j -KPlj -К ){kAql„j -kJl +Q")

qi.njPnj {kkslnj -kkJl +k4k5 +Ш4)

. _ k5 P1,0k . p _ k4 q1,0 j A ,0k _ - " - ' _

к5р20к +0^ к49207 +04 '

2. Сведение к бесконечной системе линейных алгебраических уравнений

Прежде всего, можно заметить, что с учетом четности решения по переменной у, из 12 граничных условий (3) достаточно выполнить девять (граничные значения кинематических и силовых характеристик для сторон у = Ь и у = -Ъ совпадают с точностью до знака). Далее, из данных 9 условий 5 условий можно выполнить тождественно, так как в этом случае на границе пластины соответствующая характеристика пластины представляются полными тригонометрическими рядами с неопределенными коэффициентами. Например, из условия ф^ (-а, у) = 0 и формул (5) получаем

w3

Фх,0,] (-а, у) + Фх,1,] (-a, у ) — — -8]0 (11,00Н0 («1,0]а ) + 72,00Н0 («2,0]а )) +

+ 8]0 (У1,01Н («1,0]а) + 72,01Н1' («2,0]а)) -

ж 3

Н0( ]) Т] (Рп-у)+

п—1 1—1

ж 3

+££ ,щН;( ) Т] (р„у)—0.

п—1 1—1

Отсюда, очевидно, следуют соотношения между неизвестными коэффициентами У1пк

"í1,00Hó (^,0 ,a) - Y2,00Hc'^2,0 ,a) + YWHÍ(4i,0 Ja) + 72,ЮН/^2,0 ,a) = 0

—£ Y,,c H'^a) + £ YnMw) = 0.

(10)

Кроме упомянутого граничного условия, аналогичным образом можно выполнить также условия М — О — 0 при х — а , если положить

(«1 ,0 ]Е1,0 ] 1)(11,00 Н0(«1,0 ]а) + ¥1 ,01 Н/(«1,0]а)) + + («2,0]Е2,0] - 1)(^,00Нк («2 ,0 ]а) + ^2,01Н1' («2,0 ]а) ) — 0,

3

£(«, А] +Рп] К^На' (] + у1Мн;(ч, ])) — 0, (11)

I—1

3

£( Ч, пА ]-ЩпаН'М ,п]а) + У^НХЯ])) — 0.

Аналогично, на сторонах пластины у — ±Ь можно точно выполнить условия Мху — О — 0 , если для неизвестных Х выполняются следующие соотношения:

X1,00 (Л,00A1,00 - i) HJ (Л,00Ь) + + X2,00 (P2,004,00 -1) H'j (P2,00b) = 0

3

£X^ (plMBlM +ank) H'j (pl» = 0,

(12)

3

£ X, пк (Р, пкА пк -1) Н; (Р, МЬ) — 0.

I—1

Уравнения (10)-(12) позволяют явно выразить неизвестные коэффициенты 7Лп1 через коэффициенты ^,„о,

а неизвестные Х2,пк , Х3,пк - через Х1,пк .

Действительно, введем для удобства новые неизвестные

у,пк — (Ч,,п,а) и х,,пк — Х,пЛН';(Р,пкЬ), (13)

а также матрицы, соответствующие левым частям систем (11), (12):

А

T =

Л

1 1 1

qi,njE1,nj — 1 q2,njE2,nj — 1 q3,njE3,nj — 1

V^nAn + Р* q2,rjG2,rj + Pnj q3,njG3,nj + Pnj у

Snk =

1 1 1

P1,nkA1,nk — 1 P2,nkA2,nk — 1 P3,nkA3,nk — 1

P1,nkB1,nk + ani P2,nkB2,nk + ani P3,nkB3,nk + ani y

тогда из (10) и (11) получаем следующие соотношения:

Я ,n1 =

üíkLV -det T £ Я',г

(14)

(15)

где Тх1п и 5|; пк - алгебраические дополнения к элементам с индексами (1, I) соответственно матриц Тп и ^ . Аналогичные соотношения можно найти и для слу-

чая n = 0:

fv л

71,01

V V2,01

= T

v v1,00

V V2,00 У

X2,00 S0X1,00 '

где

T = T0

q1,0 jE1,0 j q2,0 jE2,0 j

«1,0 ]Е1,0 ] «2,0 ]Е2,0 ] + 2 2ч2,0 ]Е2,0 ] + 2 у 2«1,0 ]-Е1,0 ] — 2 «1,0 ]Е1,0 ] + «2,0 ]Е2,0 ] — 2

Р1,00 ^1,00 — 1

—--^—!Г .

Р2,00 ^2,00 1

Таким образом, для решения поставленной задачи остается выполнить 4 граничных условия на сторонах пластины, используя неопределенные коэффициенты у,п0 (l =1, 2, 3) и х1пк (Л = 0, 1). При этом структура

представленного решения такова, что на границе пластины каждая из кинематических или силовых характеристик пластины представляется в виде суммы двух рядов - тригонометрического ряда и ряда по гиперболическим функциям. Используя схожесть разложения [29] на отрезке г е [—г; т] гиперболических функций

соэЬ«г по тригонометрическим функциям соб

mt

. , , „ . ж(2н — 1)/ и sinh qt по системе функции sin-:- согласно

Hk (q/ ) =

2т

я (—1T+kTk£\m — k|t)

т

H'k (qt) Tq т m=1 л2 , k \ , 2

-2 ' ks 21m ^ I + q2

получаем из краевых условий (3) бесконечную систему линейных алгебраических уравнений относительно неизвестных коэффициентов решения следующего вида:

i=1

/=1

T

4,00 , A,00 ^ -+ s0-

Р1,<

PlA

\00 + b ( E1,0 j Coth0 ( 41,0 ja )-

T11,0 E1,0j Coth1 ( ?1,0 ja )-T21,0 E2,0j Coth1 ( ?2,0ja )) Д,00 +

+ Ь ( E2,0 j Coth0 ( 4l,0ja )-T12,0 E1,0j C0th1 ( 41,0 ja )-

1 » 3

"T22,0 E2,0 j C0th1 ( 42,0ja )) 72,00 + ZZX„t Z

k=0 n=1 i=1 pin

= 0,(16)

5j0 (k* (1 + S ) *1,00 + b (41,0 jCoth0 (41,0ja) + + 41,0 jT11,0Coth1 (41,0 ja) + 42,0 jT21,0Coth1 (42,0 ja)) 71,0 + Ь ( 42,0 j Coth0 ( 42,0ja )+ 41,0 jT12,0Coth1 ( 41,0 ja ) + + 42,0jT22,0Coth1 (42,0ja)) 72,C

- Z Z (-1)k ju Z

. BLnk U

k -

Pin

,00 \ \

s

V J

= 0, (17)

ak2 (P1,00Cothj (P1,00b) + S0P2,00Cothj (P2,00b)) X XX1,00 + k 5j0 (71,00 + 72,00 ) +

bZ (-1)"+j Z

A1,00 P1,00 , S0A2,00 P2,0

(

k - k

KG,

nj l,nj

4,

У, ,n0 = 0, (18)

nj У

x1,00 +2 ZZb* Z

1 l=1 Pmj + Pl,nk

у Pmj + P1,00 Pmj + P2,00

+ b(-1)m+j Z (Ei,m (Coth0 (4j) + Coth1 (qj)) -

1=1

2 -jlE^Cot\{qimJa)flim\ylm0=0, (19)

det T" tf

PL + P1200 Kj + P22,0

*%00 + 2ZZ^"k Z

k=0 n=1 i=1 k-+pi.nk

b Z(Gi,m (Coth0 (4i,mJa) + Coth (^a))-

pm 1=1

2 3 ^ -7Т7ГZ ( w):^ = 0, (20)

detT ti J

2k

f P2

s0 P2,00

VPmj + P1,00 Pmj + P2

mj 1 F2,00 J

„Х^Х^/ Pl,nk\k Pl.nk ^nk^l,nk)^\l,nk

2ZZH)

k=0 n=1 l=1

3

Pmj + Pl,nk

+ b (-1)m+j ZZ ((qhmj - k*PmjGi,mj ) E^ (Cot^ (4j)-

2

det T" 1=1

- Coth1 (4,,mja)) +

У fa. .-k* p G. ) Coth, (a. a)f.

/ ArLi.mi rmj i,mj f 1 \rLi,mj f \i,t.

x У,, m 0 = ° =

(21)

2k 5

j 0

41,0 j71,00 42,0 j72,00 ^ 2 2 + 2 2 v 41,0 j + Um0 42,0j + Um0

3

i=1

, TV/' l^'V* 4l,nj (k 4l,nj - k2P"jGl,nj )

+ 2Z(-1) Z-__2 . .2-у, ,n 0 = 0 (22)

n =1

1=1

+ 42

nj

2k 5

j 0

41,0 jT11,0 , 422,0 jT2

0 j 21,0

420 j +°-2m1 422,0 j +Um1

41,0 jT12,0 , 4^,0 jT2

0 j^ 22,0

420j+«m1 422,0j+«m1

71,00 +

л

72,0

- aZ( k2 Pi

,m1 k am1Bl,m1 ) Coth j (Pi

,m1 . ///1X ///

i=1

+ 2L{~4 la I^TL-_.2 . 2-lM,n~

n=1

l=1

V

det T„ t!

2 2 a , + 4. .

m1 ±i,nj

4l,nj ( k*4ln - k2PnjG/,nj )

am + 4i, n

У1 ,n0 = 0,

(23)

Я.(2) (-l)X,

где k' =k-k6, Coth. (z) = —i— , =

h; (z)

3. Численные результаты

Для исследования и решения полученной бесконечной системы использовался пакет Mathematica. Согласно методу редукции полученная бесконечная система (16)-(23) сводилась к конечной системе относительно первых N неизвестных %nk и у, п0, при этом в расчетах

удерживались первые (5N + 3) уравнения системы в случае симметричных относительно оси Oy колебаний (j = 0) и (5N + 1) уравнений в случае антисимметричных колебаний (j =1). Определитель конечной системы служил в качестве дисперсионного уравнения для определения собственных частот колебаний, а нетривиальное решение - для определения собственных форм колебаний. В табл. 1 представлена сходимость метода редукции при увеличении порядка аппроксимирующей системы.

В табл. 2 приведены первые собственные частоты квадратной изотропной пластины в сравнении с результатами [10] и [12]. Заметим, что в [10] используется модифицированный вариационный подход (DSC метод), в работе [12] классический метод - Релея - Ритца. Полученные результаты оказываются в хорошем соответствии с данными [10] и [12], что позволяет сделать вывод об адекватности предложенного подхода.

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

+

+

k=0 n=1

l=1

+

+

+

X

J

Таблица 1

Сходимость метода редукции на первой собственной частоте Ц = 0.9301 для изотропной квадратной пластины при v = 0,3, к = 5/6, h/a = 0,1

Table 1

Convergence of reduction method for the first fundamental natural frequency Ц = 0.9301 of the isotropic square plate with v = 0.3, к = 5/6, h/a = 0.1

колебаниям (/ = 1). Сверяя узловые линии построенных собственных форм с известными в литературе формами для тонких пластин [30], можно заметить, что отличия незначительны.

Таблица 3

Первые собственные частоты консольно-защемленной квадратной пластины при к — 0,8601 Е — 60,7 ГПа,

G12 = G13 = G23 = 12 ГПа,

v12 = 0,23,, v21 = 0,094

n 1 2 3 5

Л„0 при N = 10 -0,32742 0,195743 -0,152986 -0,114373

Л„0 при N = 20 -0,351037 0,209318 -0,164008 -0,124418

Л„0 при N = 40 -0,357335 0,212983 -0,167122 -0,127549

Таблица 2

4Q

Первые собственные частоты Q = —— изотропной

71"

квадратной пластины при v = 0,3, к = 5/6, h/a = 0,1

Table 2 4Q2

First fundamental natural frequencies Q = ——

71"

of the isotropic square plate with v = 0.3, к = 5/6, h /a = 0.1

1 2 3 4 5 6 7 8

Представленное 0,3505 0,8170 2,0350 2,5839 2,8621 4,8172 5,4788 5,7741

решение

[10] 0,3735 0,8403 2,0502 2,5997 2,8752 4,8285 5,4885 5,7860

Г121 0,3476 0,8168 2,0356 2,5836 2,8620 4,8162 5,4834 5,7769

В табл. 3 представлены первые десять собственных частот пластины из эпоксидного стекла при вариации ее толщины. Можно заметить, что при малой толщине пластины к /а — 0.01 собственные частоты колебаний в рамках теории Рейснера - Миндлина мало отличаются от соответствующих частот, вычисленных согласно теории тонких пластин. С увеличением толщины пластины различие между собственными значениями согласно теории толстых пластин и теории тонких пластин несколько увеличивается, но не превосходит 3 %.

Нетривиальное решение бесконечной системы (16)-(23) на собственной частоте колебаний позволяет построить согласно (5) собственные формы колебаний пластины. На рис. 1 представлены первые собственные формы колебаний квадратной консольно-защемленной пластины в случае ортотропного материала, представленного в табл. 3. Первая, третья и десятая формы колебаний соответствуют симметричным относительно оси Оу колебаниям (в приведенных формулах / = 0), вторая, пятая и седьмая формы колебаний - антисимметричным

Table 3

First fundamental natural frequencies of the cantilevered square plate with k = 0.8601 E = 60.7 GPa ,

G12 = G = G3 = 12 GPa ,

v12 = 0.23,

v21 = 0.094

1 2 3 4 5 6 7 8 9 10

h/a = 0,01 0,9363 1,3124 2,1443 2,3484 2,6147 3,2433 3,2674 3,9217 4,1035 4,4660

h/a = 0,05 0,9357 1,3080 2,1368 2,3417 2,6021 3,2202 3,2540 3,8955 4,0708 4,4362

h/a = 0,10 0,9346 1,2998 2,1212 2,3218 2,5725 3,1693 3,2197 3,8198 3,9877 4,3885

Заметим, что схожая ситуация наблюдается и для других граничных условий (например, полностью защемленной пластины и пластины со свободными краями [28]). Более того, соответствующие по номеру собственные формы в случае изотропного материала имеют тот же качественный характер, что и формы на рис. 1, количественные различия можно увидеть, например, на рис. 2, где представлены линии уровня первых двух собственных форм для ортотропного и изотропного материалов. Можно заметить, что в приведенном случае собственные формы колебаний почти совпадают.

Чтобы оценить, насколько ортотропные свойства материала влияют на спектр собственных значений консольно-защемленной пластины, рассмотрим пластину с защемлением по стороне у — —Ь из того же материала Е — 60,7 ГПа, — £13 — 623 —12 ГПа, [ПW1] у12 — 0,23, у21 — 0,094, что и в предыдущих примерах.

В этом случае материал уже оказывается более жестким в направлении, перпендикулярном защемленной стороне, что сказывается как на значениях собственных частот (табл. 4), так и на отвечающих им собственных формах колебаний (рис. 3[ПW2]).

Действительно, если первая частота колебаний в табл. 4 почти совпадает для обоих случаев защемления, то для второго типа защемления (по у — -Ь ) собственные частоты лежат существенно выше. Кроме этого, меняется и четность ряда форм колебаний; так, четвертая и пятая моды колебаний являются антисимметичной и симметричной, в то же время при защемлении пластины по стороне х — -а они соответственно симметричная и антисимметричная. На рис. 3 представлены вторая, третья, четвертая и пятая формы колебаний, соответствующих случаю защемления орто-тропной пластины по стороне у — -Ь .

Рис. 1. Собственные формы колебаний консольно-защемленной толстой пластины: а, б, в, г, д, е - первая, вторая,

третья, пятая, седьмая и десятая формы колебаний

Fig. 1. The natural modes of the cantilevered thick plate; (a), (b), (c), (d), (e), (f) are the first, second, third, fifth,

seventh and tenth shape modes, respectively [nW3]

■LO -0.? 0.0 0.3 10 -10 -0 5 0.0 0.5 1.0

а б

Рис. 2. Линии уровня собственных форм колебаний консольно-защемленной толстой пластины: а, б - первая и вторая формы колебаний (сплошная линия - ортотропный материал с данными из табл. 3, пунктирная линия - изотропный материал при v = 0.3)

Fig. 2. Level lines of natural vibration modes of the cantilevered thick plate; (a), (b) are the first and second shape modes (the solid line shows the orthotopic material with data from Table 3, the dotted line shows the isotropic material at v = 0.3)

Таблица 4

Первые собственные частоты консольно-защемленной толстой квадратной пластины (E = 60,7 ГШ, G = G = G23 = 12 ГШ, v12 = 0,23, v21 = 0,094; h/a = 0,10) в случае защемления

различных сторон пластины

Table 4

First fundamental natural frequencies of the cantilevered thick square plate for the case of different clamped sides (к = 0.8601 E = 60.7 GPa , G12 = G13 = G23 = 12 GPa , v12 = 0.23, v21 = 0.094)

1 2 3 4 5

0.9346 1.2998 2.1212 2.3218 2.5725

L V Ч V " 0.9348 1.5301 2.3255 2.8626 3.0903

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

Заключение

Таким образом, в работе впервые на основе метода суперпозиции построено аналитическое решение задачи о свободных колебаниях консольно-защемленной орто-тропной пластины. Данное решение, в отличие от вариационного подхода, точно удовлетворяет уравнениям колебаний, 5 из 9 граничных условий также выполнены

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

Ш =2.8626

в

Рис. 3. Собственные формы колебаний консольно-защемленной Fig. 3. The natural modes of the cantilevered thick plate; (a),

Благодарности

Исследование выполнено при финансовой поддержке РФФИ и города Севастополя в рамках научного проекта № 18-41-920001.

Библиографический список

1. Young D. Vibration of rectangular plates by the Ritz method // J. Appl. Mech. - 1950. Vol. 17, no 4. - P. 448-453.

2. Plass H.J., Gaines J.H., Newsom C.D. Application of Reissner's variational principle to cantilever plate deflection and vibration problems // J. Appl. Mech. - 1962. - Vol. 29, no 1. -P. 127-135.

3. Austin R.N., Caughfield D.A., Plass H.J. Application of Reissner's variational principle to the vibration analysis of square flat plates with various root support conditions // Developments in theoretical and applied mechanics. - N.Y.: Plenum press, 1963. -P. 1 -24.

4. Barton M.V. Vibration of rectangular and skew cantilever plates // J. Appl. Mech. - 1951. - Vol. 18, no 1. - P. 129-134.

5. Martin A.I. On the vibration of a cantilever plate // Quart. J. Appl. Math. - 1956. - Vol. 9. - P. 94-102.

б

Q5 = 3.0903

г

толстой пластины: а, б, в, г - вторая, третья, четвертая и пятая (b), (c), (d) are the second, third, fourth, fifth shape modes

Acknowledgments

The study was carried out with the financial support of the Russian Foundation for Basic Research and the city of Sevastopol within the framework of the scientific project No. 18-41-920001.

6. Claassen R.W., Thome C.J. Vibrations of a rectangular cantilever plate // J. Aerospace Sci. - 1962. - Vol. 29, no 11. -P. 1300-1305.

7. Mindlin R. Influence of rotatory inertia and shear on fleuxural motion of isotropic elastic plates // ASME Journal Applied Mechanic. - 1951. - Vol. 18. - P. 31-38.

8. Plunkett R. Natural frequencies of uniform and nonuniform rectangular cantilever plates // J. Mech. Eng. Sci. -1963. - Vol. 5, no 2. - P. 146-156.

9. Lim C. W., Liew K.M., Kitipornchai S. Numerical aspects for free vibration of thick Part I: Formulation and verification plates // Comput. Methods Appl. Mech. Eng. - 1998. - Vol. 156. -P. 15-29.

10. Liew K.M., Xiang Y., Kitipornchai S., Transverse vibration of thick rectangular plates - I. Comprehensive sets of bounda-

ry conditions //Computers& Structures - 1993. - Vol. 49(2). -P. 1-29.

11. Chung J.H., Zhou D. Vibration of moderately thick rectangular plates in terms of a set of static Timoshenko beam functions // Computers& Structures. - 2000. - Vol. 78(6). -P. 757-768.

12. Xiang Y., Lai S.K., Zhou L. DSC- element method for free vibration analysis of rectangular Mindlin plates // Int. J. of Mech. Science. - 2010. - Vol. 52. - P. 548-560.

13. Seok J., Tiersten H.F., Scarton H.A. Free vibrations of rectangular cantilever plates. Part 1: out-of-plane motion // Journal of Sound and Vibration - 2004. - Vol. 271, no. 22. - P. 131-146.

14. Gorman D. J. Free vibration analysis of Mindlin plates with uniform elastic edge support by the superposition method // Journal of Sound and Vibration. - 1997. - Vol. 207, no. 3. -P. 335-350.

15. Kolarevic N., Nefovska-Danilovic M., Petronijevic M. Dynamic stiffness elements for free vibration analysis of rectangular Mindlin plate assemblies // Journal of Sound and Vibration. -2015. - Vol. 359 - P. 84-106.

16. Shear deformable dynamic stiffness elements for a free vibration analysis of composite plate assemblies - Part I: Theory / M. Nefovska-Danilovic, N. Kolarevic, M. Marjanovic', M. Petronijevic // Composite Structures. - 2017. - Vol. 159. -P. 728-744.

17. Free vibration study of sandwich plates using a family of novel shear deformable dynamic stiffness elements: limitations and comparison with the finite element solutions / M. Marjanovic, N. Kolarevic, M. Nefovska-Danilovic, M. Petronijevic // Thin-Walled Structures. - 2016. - Vol. 107. - P. 678-694.

18. Иванова Е.А Асимптотический и численный анализ высокочастотных свободных колебаний прямоугольных пластин // МТТ. - 1998. - № 2. - C. 163-174.

19. Usarov М. Dynamic Design of Thick Orthotropic Cantilever Plates with Consideration of Bimoments // World Journal of Mechanics. - 2016. - No. 6. - P. 341-356.

20. On the Simple and Mixed First-Order Theories for Plates Resting on Elastic Foundations / A.M. Zenkour, M.N.M. Allam,

References

1. Young D. Vibration of rectangular plates by the Ritz method. J. Appl. Mech, 1950, vol. 17, no 4, pp. 448 - 453.

2. Plass H.J., Gaines J.H., Newsom C.D. Application of Reissner's variational principle to cantilever plate deflection and vibration problems. J. Appl. Mech., 1962, vol. 29, no 1, pp. 127-135.

3. Austin R.N., Caughfield D.A., Plass H.J. Application of Reissner's variational principle to the vibration analysis of square flat plates with various root support conditions. Developments in theoretical and applied mechanics. N.Y.: Plenum press, 1963. pp. 1-24.

4. Barton M.V. Vibration of rectangular and skew cantilever plates. J. Appl. Mech, 1951, vol. 18, no 1, pp. 129-134.

5. Martin A.I. On the vibration of a cantilever plate. Quart. J. Appl. Math., 1956, vol.9, pp. 94-102.

6. Claassen R.W., Thorne C.J. Vibrations of a rectangular cantilever plate. J. Aerospace Sci., 1962, vol. 29, no 11, pp. 1300 -1305.

7. Mindlin R. Influence of rotatory inertia and shear on fleuxural motion of isotropic elastic plates. ASME Journal Applied Mechanic, 1951, vol.18, pp. 31-38.

M.O. Shaker, A.F. Radwan // Acta Mechanica - 2011. -Vol. 220. - P. 33-46.

21. Torabi K., Afshari H. Vibration analysis of a cantile-vered trapezoidal moderately thick plate with variable thickness // Engineering Solid Mechanics. - 2017. - Vol. 5(1). - P. 71-92.

22. Huang C.S., Leissa A.W., Chang M.J. Vibrations of skewed cantilevered triangular, trapezoidal and parallelogram Mindlin plates with considering corner stress singularities // International Journal for Numerical Methods in Engineering. - 2005. -Vol. 62. - P. 1789-1806.

23. Free vibration analysis of plates using least-square-based on finite difference method / M. Huang, X.O. Ma, T. Sakiyama, M. Matuda, C. Morita // Journal of Sound and Vibration. - 2005. -Vol. 288. - P. 931-955.

24. Nguyen-Xuan H., Liu G.R., Thai-Hoang C. An edge-based smoothed finite element method (ES-FEM) with stabilized discrete shear gap technique for analysis of Reissner-Minslin // Computer Methods in Applied Mechanics and Engineering. -2010. - Vol. 199. - P. 471-489.

25. Sukhoterin M., Baryshnikov S., Aksenov D. Free Vibration Analysis of Rectangular Cantilever Plates Using the Hyperbolic-Trigonometric Series // American Journal of Applied Sciences. - 2016. - Vol. 13 (12). - P. 1442-1451.

26. Papkov S.O. A new method for analytical solution of inplane free vibration of rectangular orthotropic plates based on the analysis of infinite systems // Journal of Sound and Vibration. -2016. - Vol. 369. - P. 228-245.

27. Papkov S.O., Banerjee J.R. A new method for free vibration and bucking analysis of rectangular orthotropic plates // Journal of Sound and Vibration. - 2015. - Vol. 339. - P. 342-358.

28. Papkov S.O. New analytical solutions for vibration problem of thick plates // PNRPU Mechanics Bulletin. - 2019. -№ 4. - P. 145-156.

29. Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. Элементарные функции. - М.: Наука. Глав. ред. физ.-мат. лит-ры, 1981. - 800 с.

30. Leissa A.W. Vibration of Plates (NASA SP-160). -Washington, DC: Govement Printing office, 1969. - 353 p.

8. Plunkett R. Natural frequencies of uniform and nonuniform rectangular cantilever plates. J. Mech. Eng. Sci., 1963, vol. 5, no 2, pp. 146-156.

9. Lim C. W., Liew K.M., Kitipornchai S. Numerical aspects for free vibration of thick Part I: Formulation and verification plates. Comput. Methods Appl. Mech. Eng., 1998, 156, pp. 15-29

10. Liew K.M., Xiang Y., Kitipornchai S., Transverse vibration of thick rectangular plates - I. Comprehensive sets of boundary conditions. Computers& Structures, 1993, 49(2), pp. 1-29.

11. Chung J.H., Zhou D. Vibration of moderately thick rectangular plates in terms of a set of static Timoshenko beam functions. Computers& Structures, 2000, 78(6), pp. 757-768.

12. Xiang Y., Lai S.K., Zhou L. DSC- element method for free vibration analysis of rectangular Mindlin plates. Int. J. of Mech. Science, 2010, 52, pp. 548-560

13. Seok J., Tiersten H.F., Scarton H.A. Free vibrations of rectangular cantilever plates. Part 1: out-of-plane motion. Journal of Sound and Vibration, 2004, vol. 271, 22, pp. 131-146.

14. Gorman D. J. Free vibration analysis of Mindlin plates with uniform elastic edge support by the superposition method. Journal of Sound and Vibration, 1997, vol. 207, no.3, pp.335-350.

15. Kolarevic N., Nefovska-Danilovic M., Petronijevic M. Dynamic stiffness elements for free vibration analysis of rectangular Mindlin plate assemblies. Journal of Sound and Vibration, 2015, 359, pp. 84-106

16. Nefovska-Danilovic M., Kolarevic N., Matjanovic' M., Petronijevic M. Shear deformable dynamic stiffness elements for a free vibration analysis of composite plate assemblies - Part I: Theory. Composite Structures, 2017, 159, pp.728-744.

17. Marjanovic M., Kolarevic N., Nefovska-Danilovic M., Petronijevic M. Free vibration study of sandwich plates using a family of novel shear deformable dynamic stiffness elements: limitations and comparison with the finite element solutions. Thin-Walled Structures, 2016, 107, pp. 678-694.

18. Ivanova E. A. Asimptoticheskiy i chislennyy analiz vysokochastotnykh svobodnykh kolebaniy pryamougol'nykh plastin [Asymptotic and numerical analysis of high-frequency free vibrations of rectangular plates] Mechanics of Solids, 1998, no. 2, pp. 163-174

19. Usarov M. Dynamic Design of Thick Orthotropic Cantilever Plates with Consideration of Bimoments. World Journal of Mechanics, 2016, no 6, pp. 341-356.

20. Zenkour, A.M., Allam, M.N.M., Shaker, M.O. and Rad-wan, A.F. On the Simple and Mixed First-Order Theories for Plates Resting on Elastic Foundations. Acta Mechanica, 2011, 220, pp. 33-46

21. Torabi K., Afshari H. Vibration analysis of a cantile-vered trapezoidal moderately thick plate with variable thickness. Engineering Solid Mechanics, 2017, 5(1), pp. 71-92.

22. Huang C.S., Leissa A.W., Chang M.J. Vibrations of skewed cantilevered triangular, trapezoidal and parallelogram

Mindlin plates with considering corner stress singularities. International Journal for Numerical Methods in Engineering, 2005, 62, pp. 1789-1806.

23. Huang M., Ma X.O., Sakiyama T., Matuda M., Morita C. Free vibration analysis of plates using least-square-based on finite difference method. Journal of Sound and Vibration, 2005, 288, pp. 931-955.

24. Nguyen-Xuan H., Liu G.R., Thai-Hoang C. An edge-based smoothed finite element method (ES-FEM) with stabilized discrete shear gap technique for analysis of Reissner-Minslin. Computer Methods in Applied Mechanics and Engineering, 2010, 199, pp. 471-489.

25. Sukhoterin M., Baryshnikov S., Aksenov D. Free Vibration Analysis of Rectangular Cantilever Plates Using the Hyperbolic-Trigonometric Series. American Journal of Applied Sciences, 2016, 13 (12), pp. 1442.1451.

26. Papkov S.O. A new method for analytical solution of inplane free vibration of rectangular orthotropic plates based on the analysis of infinite systems. Journal of Sound and Vibration, 2016, 369, pp. 228-245.

27. Papkov S.O., Banerjee J.R. A new method for free vibration and bucking analysis of rectangular orthotropic plates. Journal of Sound and Vibration, 2015, 339, pp. 342 - 358.

28. Papkov S.O. New analytical solutions for vibration problem of thick plates. PNRPU Mechanics Bulletin, 2019, no 4, pp. 145-156.

29. Prudnikov A.P., Brychkov Yu. A., Marichev O.I. Integraly i ryady [Integrals and series], Moskow, Nauka, 1985, 800 p.

30. Leissa A.W. Vibration of Plates (NASA SP-160). Washington, DC: GovementPrinting office, 1969, 353 p.