CC BY
51
УДК 539.3, 519.6
Вестник СПбГУ. Сер. 1. Т. 2(60). 2015. Вып. 2
НЕОСЕМИММЕТРИЧНАЯ ПОТЕРЯ УСТОЙЧИВОСТИ НЕОДНОРОДНЫХ КОЛЬЦЕВЫХ ПЛАСТИН*
Е. Б. Воронкова, К. А. Минина
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7/9
В работе рассматривается устойчивость осесимметричных форм равновесия неоднородных круглых и кольцевых пластин, загруженных нормальным давлением. Внешний край пластины закреплен от поворотов, но точки края свободно смещаются в радиальном и окружном направлении, а внутренний — свободно смещается в направлении оси пластины, но не поворачивается. Полагается, что модуль упругости пластины меняется при движении от центра к краю. Наименьшее значение нагрузки найдено численным методом в предположении, что несимметричная составляющая решения системы носит периодический характер. Исследовано влияние отверстия в центре пластины и степени неоднородности материала на величину критической нагрузки и форму потери устойчивости. Библиогр. 5 назв. Ил. 1. Табл. 2.
Ключевые слова: кольцевая пластина, потеря устойчивости, неоднородность.
1. Введение. При больших прогибах пологие пластины и оболочки, нагруженные симметричной нагрузкой, могут потерять свою осесимметричную форму, при этом по краю пластины образуются складки [1, 2]. Существование несимметричного решения для симметрично нагруженной пластины было доказано Н.Ф.Морозовым в [3], а единственность такого решения доказана В.Пихоки в [4]. В работе [5] были найдены значения критической нагрузки, при которой для неоднородной круглой пластины возможен переход в неосесимметричное состояние.
В настоящей работе рассматривается задача о потере устойчивости осесиммет-ричной формы равновесия неоднородной кольцевой пластины, модуль упругости которой изменяется при движении от центра пластины к ее краю. Исследуется влияние отверстия в центре пластины и неоднородности материала пластины на величину критической нагрузки.
2. Постановка задачи. Рассмотрим кольцевую пластину внешним радиусом К, внутренним радиусом Кп и толщиной Н, лицевая поверхность которой загружена нормальным давлением. Материал пластины полагается изотропным и неоднородным в радиальном направлении.
Уравнения равновесия пластины в безразмерных переменных имеет вид [5]
д1(т)ААю + д[ (т)Ь+(ю) + д"(т)Ь+ (ю) = р + Ь(ю,Е),
д2(т)ААЕ + д'2(т)Ь-(Е) + д'2(т)Ь- (Г) = -Ь(т,т)/2, (1)
()' = д()/дт, () = 8()/дв.
Здесь ю(т,в), Г(т,в) —неизвестные безразмерные нормальный прогиб и функция усилий; т, в — полярные координаты, введенные на срединной поверхности пластины; д\(т) = Е(т)/Еау (т) —достаточно гладкая функция, определяющая закон изменения модуля упругости в плоскости пластины, где Е(т) — модуль упругости, Еау = //Е(т)т Лгйв/Б — среднее значение модуля упругости, $ — площадь пластины; я
* Работа выполнена при финансовой поддержке РФФИ (гранты №15-01-06311, №13-01-00523.
д2(г) = 1/д^г); Д — оператор Лапласа, Ь± (г = 1, 2) —дифференциальные операторы
у) = ж'' (у'/г + у/г2) + у'' (ж'/г + ж/г2) — 2(ж/г)'(у/г)', Ь±(у) = 2у''' +(2 ± ^)у''/г + 2 (у)' /г2 — у'/г2 — 3у/г3, Ь±(у) = у'' ± ^(у'/г + у/г2) .
Безразмерные величины связаны с размерными соотношениями
4= Т"» 4= 4= Т"»4 77* 4=
где V — коэффициент Пуассона пластины.
Положим, что внешний край пластины закреплен от поворотов, но точки края свободно перемещаются в радиальном и окружном направлении, а внутренний — свободно смещается в направлении оси симметрии пластины, но не поворачивается. Тогда на внешнем крае пластины (г =1) предполагается равенство нулю растягивающих и сдвигающих усилий, на внутреннем (г = 6) обращается в нуль поперечное реактивное усилие, складывающееся из поперечной силы и производной от крутящего момента.
В терминах нормального перемещения т и функции усилий Р граничные условия запишутся в виде
г = 1 : т = т' = Р '/г + Р/г2 = — /г^ = 0;
г = 6 : т' = — ( Р/г ) = 0,
1 , (2)
w''' + w''/r - w'/r2 + (2 - г/) (W)' /r2 - (3 - v)W/r3 = 0,
Для сплошной пластины, учитывая ограниченность искомых решений, в центре пластины при r = 0 полагаем w' = F' = 0.
Схема решения задачи подробно описана в [2, 5]. Полагая, что для малых значений нагрузки p система (1)-(2) имеет только симметричное решение, а несимметричное решение этой системы появляется при возрастании нагрузки, представим решение в виде
w(r, 0) = ws(r) + wn(r)cos n0, F(r, 0) = Fs(r) + Fn(r)cos n0, (3)
где функции ws, Fs описывают докритическое симметричное решение, а wns(r, 0) = wn(r)cos n0, Fns (r, 0) = Fn (r) cos n0 — закритическое состояние пластины (n — в окружном направлении число волн, образовавшихся после потери устойчивости).
После разделения переменных (3) исходная система (1)-(2) распадается на две: нелинейную, для нахождения симметричного решения ws (r), Fs (r), и линейную систему уравнений относительно wn(r), Fn(r), так как функции wns, Fns являются малыми сразу после перехода пластины в неосесимметричное состояние. Для каждого числа волн в окружном направлении n будем искать такие значения нагрузки pn, при которых существуют отличные от нуля функции wn(r), Fn(r). Критической нагрузке соответствует pcr = min pn.
n
3. Результаты. Проведены серии расчетов для неоднородных сплошной и кольцевой пластин при изменении модуля упругости пластины по линейному, E = E(1)(1 +
д\г), и экспоненциальному, Е = Е(2)е92Г, законам. Параметры Е(г) и ^ (г = 1, 2) выбирались так, чтобы среднее значение модуля упругости пластины Еау оставалось постоянным. Полученное в [5] критическое значение нагрузки для однородной пластины равно р°сг = 62598, а соответствующее этой нагрузке волновое число п = 14. Для неоднородных сплошных пластин результаты расчетов приведены на рисунке. Значение 0 параметра д соответствует пластине с постоянным модулем упругости.
1.1
P/PO
0.9
0.8
-0.4 -0.2 0 0.2 0.4 0.6
q
Изменение критической нагрузки в зависимости от неоднородности пластины q, po — критическая нагрузка для однородной пластины.
Расчеты показывают, что если модуль упругости пластины уменьшается к ее краю, то потеря устойчивости осесимметричных форм равновесия происходит при более низкой нагрузке и с образованием большего числа складок в окружном направлении, чем для однородной пластины. Чем быстрее убывает модуль упругости, тем ниже становится величина критической нагрузки. При увеличении модуля упругости к центру пластины критическая нагрузка возрастала.
Безразмерные значения критической нагрузки p и число волн n, образующихся по краю пластины при переходе в неосесимметричное состояние, для однородной кольцевой пластины при различных радиусах центрального отверстия приведены в табл. 1, для неоднородной кольцевой пластины при изменении модуля упругости по закону E = Eoeqr —в табл. 2. В обоих случаях p°r соответствует критической нагрузке для сплошной пластины (S = 0)
Таблица 1. Критическая нагрузка для однородной кольцевой пластины
S = 0 й = 0.05 й = 0.1 й = 0.15
Рсг/Рсг 1 1.03 1.08 1.15
п 14 13 12 12
Видно, что потеря устойчивости осесимметричных форм равновесия при увеличении радиуса центрального отверстия происходит с образованием меньшего числа складок в окружном направлении.
Как и для сплошной пластины, при убывания модуля упругости от центра пластины к ее краю переход (д < 0) в неосеммитричное состояние происходит при более низкой нагрузке, чем для однородной пластины, а при увеличении (д > 0) — при существенно более высокой.
ж" у,'«»""
У
er" S и ✓
■O-E = E^V -B-E = E02)( 1 + qr)
Таблица 2. Критическая нагрузка для неоднородной кольцевой пластины при & = 0.05, Е = Еоедг)
q = -0.5 q = -0.1 q = 0 q = 0.1 q = 0.5 q= 1
Per/Per 0.75 1.002 1.03 1.06 1.17 1.34
n 13 13 13 13 13 13
Потеря устойчивости осемимметричных форм равновесия локализована вблизи внешнего края пластины и связана с появлением в окрестности внешнего края сжимающих окружных усилий [2, 5]. Этим объясняется влияние на критическую нагрузку свойств пластины в окрестности внешнего края.
5. Заключение. В работе численным методом исследована потеря устойчивости осесимметричных форм равновесия неоднородных кольцевых пластин. Найдено критическое давление, при котором возможна бифуркация пластины в неосесиммет-ричное состояние. Показано влияние степени неоднородности и размера центрального отверстия пластины на изменение критической нагрузки. Произведен численный анализ критической нагрузки для различных законов изменения модуля упругости пластины, размеров внутреннего края в случае кольцевой пластины.
Литература
1. Панов Д.Ю., Феодосьев В. И. О равновесии и потере устойчивости пологих оболочек при больших прогибах // ПММ. Т. XII. 1948. C. 389-406.
2. Cheo L.S., Reiss E.L. Unsymmetric wrinkling of circular plates // Quart. Appl. Math. N31. 1971. P. 75-91.
3. Морозов Н. Ф. К вопросу о существовании несимметричного решения в задаче о больших прогибах круглой пластины, загруженной симметричной нагрузкой // Изв. Высш. Уч. Заведений, Математика, 1961, №2. С. 126-129.
4. Piechocki W. On the non-linear theory of thin elastic spherical shells // Arch. Mech. Stos. №21, 1969. P. 81-101.
5. Бауэр С.М., Воронкова Е.Б., Романова А. А. О потере устойчивости симметричных форм равновесия круглых пластин под действием нормального давления // Вестн. С.-Петерб. ун-та. Сер. 1: Математика. Механика. Астрономия. 2012. Вып. 1. С. 80-85.
Статья поступила в редакцию 26 марта 2015 г. Сведения об авторах
Воронкова Ева Боруховна — кандидат физико-математических наук, доцент; e.voronkova@spbu.ru Минина Ксения Андреевна — аспирант, ksushik@mail.ru
BUCKLING OF AXISYMMETRIC EQUILIBRIUM STATES OF ANNULAR PLATES
Eva B. Voronkova, Kseniya A. Minina
St.Petersburg State University, Universitetskaya nab., 7/9, St.Petersburg, 199034, Russian Federation; e.voronkova@spbu.ru
The stability of axisymmetric equilibrium states of an isotropic non-homogeneous plate under uniform pressure is considered. The unsymmetric part of the solution is sought in terms of multiples of the harmonics of the angular coordinates. A numerical method is employed to obtain the lowest load value, which leads to the appearance of waves in the circumferential direction. It is shown that if the elasticity modulus decreases away from the center of a plate, the critical pressure for unsymmetric buckling is sufficiently lower than for a plate with constant mechanical properties. Refs 5. Figs 1. Tables 2. Keywords: circular plate, buckling, non-homogeneity.
References
1. Panov D.U., Feodos'ev V. I., "Equilibrum and loss of stability of shallow shells with large deflections", PMM 12, 389-406 (1948) [in Russian].
2. Cheo L. S., Reiss E. L., "Unsymmetric wrinkling of circular plates", Quart. Appl. Math. N31, 75-91 (1971).
3. Morozov N. F., "On the existence of a non-symmetric solution in the problem of large deflections of a circular plate with a symmetric load", Izv. Vyssh. Uchebn. Zaved. Mat. N2, 126-129 (1961) [in Russian].
4. Piechocki W., "On the non-linear theory of thin elastic spherical shells", Arch. Mech. Stos. N21, 81-101 (1969).
5. Bauer S. M., Voronkova E. B., Romanova A. A., "On the unsymmetric buckling of circular plates under normal pressure", Vestnik St.Petersburg University. Ser. 1 Issue 1, 80-85 (2012).
