CC BY
65
УДК 539.3, 519.6
Вестник СПбГУ. Сер. 1. Т. 3(61). 2016. Вып. 1
О СОБСТВЕННЫХ ЧАСТОТАХ КОЛЕБАНИЙ ТРАНСВЕРСАЛЬНО-ИЗОТРОПНЫХ КРУГЛЫХ ПЛАСТИН*
С. М. Бауэр, Е. В. Воронкова
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
В работе обсуждается влияние свойств материала круглой трансверсально-изотропной пластины на частоты ее собственных колебаний. Частоты собственных колебаний однородной пластины найдены с помощью двух неклассических теорий пластин и оболочек с учетом нормальных сил инерции и инерции вращения. Для исследования влияния неоднородных по радиусу свойств пластины проведены расчеты частот колебаний в конечно-элементном пакете Сотэо1 МиШрЬуБЮБ 5.0. Проведенные расчеты показывают, что неоднородность пластины оказывает наибольшее влияние на первую частоту колебаний, в то время как различие тангенциальных и нормальных модулей упругости оказывает большее влияние на высокочастотные колебания пластины. Библиогр. 2 назв. Ил. 1. Табл. 1.
Ключевые слова: круглая пластина, трансверсальная изотропия, колебания.
1. Введение. В работе рассматривается задача о колебаниях трансверсально-изотропной круглой пластины и исследуется влияние свойств материала пластины на частоты ее собственных колебаний. Классическая теория пластин Кирхгофа—Лява (КЛ) учитывает только свойства пластины в срединной поверхности, поэтому частоты, получаемые для изотропных и трансверсально-изотропных пластин по классической теории, одинаковы. Теория анизотропных оболочек С. А. Амбарцумяна [1] позволяет учесть влияние деформаций сдвига в направлении толщины пластины на ее напряженно-деформированное состояние. Теория анизотропных пластин и оболочек, разработанная В. А. Родионовой, Б. Ф. Титаевым и К. Ф. Черныхом [2], в общем случае позволяет учитывать не только поперечные сдвиги, но и деформируемость нормали к срединной плоскости. В данной статье задача о нахождении частот собственных колебаний решается с помощью теорий [1] и [2], уточняющих теорию КЛ. Для исследования влияния неоднородных по радиусу свойств пластины проводятся расчеты частот колебаний в конечно-элементом пакете Comsol Multiphysics 5.0.
2. Постановка задачи. Рассматривается задача о нахождении частот собственных колебаний трансверсально-изотропной круглой пластины радиуса Д, толщины к. Материал пластины подчиняется закону Гука [2], при этом плоскость пластины является плоскостью изотропии:
ац = Ецвы + Е^в^ + Егквии, % = 3 = к, а^ = Сцвц, г = 3, (1)
где
Е1Е1 - Е3Е1(Е1и12 + Ез^2з) Е^з^з(1 + Ей --р-, Е12 ----, Ь13 --р-,
е,СП е,СП е,СП
Е33 = д^за-^ц Еаеп = (1+_ Щ2) _ 2Ез^
Е,-
den
* Работа выполнена при финансовой поддержке гранта РФФИ № 15-01-06311 и с использованием оборудования ресурсного центра Научного парка СПбГУ «Обсерватория экологической безопасности».
(¡5 Санкт-Петербургский государственный университет, 2016
Здесь Ег{г = г, в, г) — модули упругости в направлении координаты г, (г, в, г) —введенные цилиндрические координаты, О^ —модули сдвига в плоскости (г — ]), — коэффициенты Пуассона. Для трансверсально-изотропной пластины О12 = Е\/2(1 + ^12), 013 = О23.
Система уравнений колебаний трансверсально-изотропной однородной пластины по уточненной теории С. А. Амбарцумяна [1] записывается в виде
д2 цп
ВААи> + рЫ1 -кА)—г =0, = 0, (2)
дЬ2
где
Ь2 Л О12 \ Л2 10 О13 2т;--^13 , <5
10(1 — ^Н О13 У' ь2 О12'
V — нормальный прогиб, Е — функция, связывающая прогиб и перерезывающие усилия, Д — оператор Лапласа, р — плотность материала, О = Е1Ь3/12(1 — ^2) —цилиндрическая жесткость пластины.
Характеристическое уравнение для определения частот собственных колебаний жестко заделанной по контуру пластины имеет вид
4-1Ы _ д-Лг-1(А) ^
1п{р) Лг(А)
2<52Д2(1 -г/12) ' Ам
+ - „12) + -Тйт) = (3)
Частоты колебаний этой пластины связаны с корнями А„ т уравнения (3) соотношениями
А __А п,т I __А„.т
<т ) 1 / , , .1 -
-А / и о4 ' гп,т ,-
Здесь 1п — функции Бесселя, п, т — число узловых диаметров и окружностей соответственно.
Уравнение, описывающее колебания пластины по классической теории пластин КЛ, и соответствующее ему характеристическое уравнение могут быть получены из (2) и (3), если положить к = 0 и 6
3. Колебания по теории Родионовой—Титаева—Черныха. Теория анизотропных пластин и оболочек, разработанная В. А. Родионовой, Б.Ф.Титаевым и К.Ф.Черныхом (PTЧ) [2], в общем случае позволяет учитывать деформируемость нормали к срединной плоскости. При построении теории полагается, что касательные напряжения распределены по толщине по закону полинома третьей степени, а нормальное напряжение — полинома второй степени от координаты В данной статье рассматривается задача о колебаниях на основе теории, изложенной в [2], с учетом только поперечных сдвигов.
При отсутствии внешних поверхностных сил векторно-матричное дифференциальное уравнение движения с учетом нормальных сил инерции и инерции вращения имеет вид
д2у „ д2У „ д 2У „ дУ „ дУ „ ,г „ д 2У
ТУ + А12 -7Г7Ш + 22 + ~я~ + 02 -яз + ¿оо V + Аа -— = 0, 4
дг2 гдгдв г2 дв2 дг гдв дЬ2
в котором V = (ад* (г, в, Ь), 7* (г, в, ¿), 1*> (г, в, ¿)) — искомая вектор-функция. В квадратных матрицах А^- следующие элементы отличны от нуля:
<> = A&V =
*22
Aif = а255) = dlb Aif = A233) = d66,
A^3 = Al322) = di2 + d66,
.(И) _ 5
^io = —G^h, A\X> = -G
(12) _ 5
4(22) _ dn Aw - —>
432X) = -Icisfc.
4(33) _ de6
6r
4(13) _ 5r Л)2 —
10
3
^13,
4оХ) =
13,
л (23) _ л (32) _ ¿12 + ¿66
Л02 — ~A02 — > r
4(11) = -Pk, A(22)
00 3r 13.
h2 Am
4 (22) _ ~dn
^oo — —
P-
6
— /
dii =
E(k2
6(1 - ^12)'
di2 = V(2d(i,
a.
1 в '
1 - Vi2
- з«13,
iii.
A
(33) _
00
=-dee/r --G13,
r
2
После разделения переменных для функций V* = Шп(г)cos(nв)cos(wt), -у* = д1п(г)сов(пв)сов(со1), 7* = д2п(г) вт(пв) вт(^) система (4) может быть сведена к виду
У' = А(ш2)У, (5)
где У = {\Уп,д1п,д2п,'№';1,д'1п,д'2п), А(ш2) — матрица размера 6 х 6, ( )' = ■£( ). Решение системы (5) представимо виде
Y(r) = L
еЛ6 r
L-iC.
(6)
л
1 r
0
e
0
Здесь Xi — собственные числа, а L — матрица, составленная из собственных векторов матрицы A(w2), C — вектор произвольных постоянных. После подстановки (6) в соответствующие граничные условия уравнение для определения частот wnm выводится из условия существования нетривиального решения системы (5).
4. Результаты. В таблице приведены нормированные параметры частот собственных колебаний однородной пластины Xn,m = Xn,m/XKm, найденные по двум уточненным теориям оболочек [1, 2], а также методом конечных элементов (МКЭ) в пакете Comsol Multhiphysics 5.0. Предполагалось, что k/R = 1/10, а отношение модулей упругости равно Ei/E3 = Gn/Gi3 = 1 для изотропной пластины и Ei/E3 = Gn/Gi3 = 10 для трансверсально-изотропной пластины. Индексы KL, A, RTC, FEM соответствуют значениям, найденным по классической теории КЛ, теориям Амбарцумяна, уточненной теории РТЧ и с помощью МКЭ соответственно. Видно, что трансверсальная изотропия оказывает большее влияние на высшие частоты.
Расчеты частот собственных колебаний неоднородной пластины были проведены в пакете Comsol Multhiphysics 5.0. При расчетах предполагалось, что модуль упругости в плоскости изотропии пластины изменяется при движении от центра к краю пластины Ei = Eiof (r), где f (r) —достаточно гладкая положительная на отрезке [0,R] функция. При проведении расчетов параметр Ei0 выбирался так, чтобы среднее зна-
Нормированные частоты колебаний
при h/R = 1/10
Ei/Es = G12/G13 = 1
Ei/Es = G12/G13 = Ю
0.97 0.95 0.92 0.92 0.90 0.89
0.97 0.96 0.92 0.92 0.90 0.89
0.98 0.96 0.92 0.92 0.89 0.89
0.77
0.67
0.6
0.58
0.54
0.51
0.79
0.68
0.6
0.59
0.55
0.52
0.8
0.7
0.61
0.6
0.56
0.53
1 hom
• (0.0)
- (1.0)
-ш- - (2.0)
-0.5 -0.3 -0.1 0.1 0.3 0.5 Параметр неоднородности, (<|)(|
Влияние неоднородности пластины на частоты колебаний
n,m
n,m
n
m
чение модуля упругости пластины Eav оставалось постоянным:
2п R
Eav = —^ J J Ei(r)rdrd,e = const.
0 0
На рисунке на примере трех низших частот колебаний показано влияние свойств материала пластины на частоты ее собственных колебаний. Предполагалось, что модуль Юнга в плоскости изотропии изменяется по закону Ei = Ею(1 + qr/R), и его среднее значение в 10 раз превосходит модуль упругости пластины в направлении толщины. Частоты Л10 и Л2,о неоднородной пластины отличаются от соответствующих частот однородной пластины (q = 0) не более чем на 7%, а для первой (наименьшей) частоты Ло,о это изменение превышает 10%.
5. Заключение. Проведенные расчеты показывают, что различие тангенциальных и нормальных модулей упругости оказывает большее влияние на высокочастотные колебания пластины. Неоднородность пластины больше влияет на низшие частоты.
Литература
1. Амбарцумян С. А. Теория анизотропных пластин. М., 1967. 226 с.
2. Родионова В. А., Титаев Б. Ф., Черных К. Ф. Прикладная теория анизотропных пластин и оболочек. СПб.: Изд-во С.-Петерб. ун-та, 1996. 280 c.
Статья поступила в редакцию 22 октября 2015 г.
Сведения об авторах
Бауэр Светлана Михайловна — доктор физико-математических наук, профессор; s.bauer@spbu.ru Воронкова Ева Боруховна — кандидат физико-математических наук, доцент; e.voronkova@spbu.ru
ON NATURAL FREQUENCIES OF
TRANSVERSALLY ISOTROPIC CIRCULAR PLATES
Svetlana M. Bauer, Eva B. Voronkova
St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation; s.bauer@spbu.ru, e.voronkova@spbu.ru
Impact of material properties of transversely isotropic circular plates on its natural frequencies is discussed in the present paper. Two refined theories of plates are used to analyze the free vibration behavior of homogeneous plates. Both theories take into account normal and rotary inertias. Fundamental frequencies for plates with radial inhomogeneity are obtained with the help of finite element package Comsol Multiphysics 5.0. It is shown that the inhomogeneity of the plate have a profound impact on the first (lowest) frequency of the plate while the plate orthotropic influences more on the second and higher vibration modes. Refs 2. Figs 1. Tables 1.
Keywords: circular plate, transversely isotropic plate, vibrations.
References
1. Ambartsumyan S.A., Theory of Anisotropic Plates (Technomic, Stamford, 1970).
2. Rodionova V.A., Titaev B.F., Chernykh K. F., Applied Theory of Anisotropic Plates and Shells (St. Peterb. Univ., St. Petersburg, 1996) [in Russian].
