Математические модели потери устойчивости неоднородных цилиндрических оболочек от неравномерной радиальной нагрузки Текст научной статьи по специальности «Физика»

МЕХАНИКА

УДК 539.3

МАТЕМАТИЧЕСКИЕ МОДЕЛИ ПОТЕРИ УСТОЙЧИВОСТИ

НЕОДНОРОДНЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК ОТ НЕРАВНОМЕРНОЙ РАДИАЛЬНОЙ НАГРУЗКИ

Э.В. Антоненко, Т.Э. Шульга*

Саратовский государственный университет, кафедра вычислительного эксперимента в механике;

* Саратовский государственный социально-экономический университет, факультет информатики и информационных технологий

E-mail: vem@info.sgu.ru, *shulga@ssea.runnet.ru

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

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

Mathematical Models of Stability Loss of Nonuniform Cylindrical Shells Because of Nonuniform Radial Loading

E.V. Antonenko, T.E. Shulga*

Saratov State University

Chair of Mechanics Computational Experiment;

* Saratov State Socio-Economical University

Faculty of Computer Sciences and Informational Technologies E-mail: vem@info.sgu.ru, *shulga@ssea.runnet.ru

The circular cylindrical shell with variable thickness along the axis of elongation is considered. The axisymmetric radial pressure along the axis of shell is suggested. The one of values (for the law of pressure variation) which effects the stability loss of shell is determinated.

Key words: cylindrical shell, theory of shells, stability of shells, radial loading, critical pressure, shell thickness, the theory of elasticity.

Исследования устойчивости оболочек от действия неравномерного давления, по-видимому, начинаются с работ Б.О. Элмроса [1] и В.И. Моссаковского [2]. Решение некоторых задач для неоднородных оболочек с переменной радиальной нагрузкой имеется в работах [3, 4]. Однако в них рассматривались ситуации отдельно для неоднородности оболочки и отдельно от действия переменной нагрузки. На конференциях по механике деформированных тел, теории оболочек и в литературе [4-6] ставились задачи о получении зависимостей для расчета критического давления с учетом неоднородности оболочек и нагрузок, аналогичных зависимостям [7].

Задачу статической устойчивости будем решать в линейной постановке, используя классические зависимости теории оболочек [6, 8-10] и деформационных гипотез нерастяжимости поперечного сечения и отсутствия сдвигов в срединной поверхности. Применим энергетический метод в форме [9], используя две математические

© Э.В. Антоненко, Т.Э. Шульга, 2009

79

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

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

Составим функционал энергии единицы длины оболочки, включая в него наиболее весомые в энергетическом отношении составляющие: энергию изгиба поперечного сечения от кольцевых изгибающих моментов mр и энергию депланаций от осевых напряжений <x:

Г = ±/ [m^ Хр + <xó(x)ex - p(x)Rwx^]Rdp, (1)

где последнее слагаемое есть работа сил переменного вдоль оболочки давления p(x) на радиальных перемещениях, толщину оболочки 5(x) полагаем также переменной вдоль продольной оси оболочки x.

Все усилия и деформации в (1) выражаются через функцию радиальных перемещений w(x, p) и ее производные. Прогиб оболочки зададим в виде

w(x, -фп (x)cos np, (2)

n

где -фп(x) — неизвестные функции, меняющиеся по длине оболочки, n = 2,3,4,... — числа натурального ряда, соответствующие количеству волн в поперечном сечении оболочки при потере устойчивости. Далее индекс n опускается.

Кольцевые и осевые перемещения в силу принятых гипотез имеют вид

"0(x) • R ,,( ч

v = — > - sin np, u = — > -tW (x)cos np.

n n2

n n

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

1 (д2 w \ v-^ n2 - 1 . , , ^

= -R2 [ + w) = ^(x)cosnp, m = DX9,

\ г / n

du R2 ,,, .

£x = дХ = П2 ^ (x)cos n^, = Eex.

(3)

Подстановка (3) в (1) с учетом ортогональности тригонометрических функций дает

Г = 1 У ÍD(x)(n^ ^2(x) + ^^[^"(x)]2 - ^Р(Ф2(Х)\nR, (4)

^ ^ 1 v 7 R4 n4 Lr v 7J R

n 4

где учтено, что f cos2npdp = n, D(x) = ES3(x)/12(1 — д2).

Дифференциальное уравнения задачи потери устойчивости получим из условия минимума функционала (4), применяя уравнение Эйлера - Лагранжа вариационной задачи:

дГ d ( дГ \ d2 ( дГ \ п

+ = 0- (5)

д-0(ж) dx \дф'(x) J dx2 \дф''(x) Выполняя для (4) операции (5) получаем:

,JVfs , J'(x)ul„, , , 5''{x)u„, , , n4(n2 - 1) ÍE(n2 - 1)¿3(x)

,tv, s (x) ,/"/ \ ó''(x) ,,,, , n4(n2 - 1) ÍE(n2 - 1)ó3(x) , Л , , , 4,tv(x) +2_i_2^(x) + <(x) + _3e^ -фп(x) = 0. (6)

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

Для определения величины критического давления закон изменения давления в (6) необходимо представлять в форме, содержащей основной параметр этого давления, например,

Р(х) = Ро/(х). (7)

Величина р0, которую будем называть критическим давлением, зависит, как это видно из (6), от функций /(х) и граничных условий. Только для случая (7) можно говорить о задаче на собственные значения дифференциального уравнения.

Для оболочек с постоянной толщиной 5(х) = <50 при действии равномерного давления р(х) = р0 из (6) как частный случай следует известное в теории оболочек дифференциальное уравнение с постоянными коэффициентами. Оно имеет аналитическое решение, приводящее, например, к формуле П.Ф. Попковича для критического давления.

Аналитическое определение величины р0 из дифференциального уравнения с переменными коэффициентами (6) не представляется возможным, численные решения связаны со значительными трудностями.

Предлагается более простая модель рассматриваемой задачи, следующая из принципа минимума потенциальной энергии и, когда в состоянии безразличного равновесия оболочки при потере устойчивости работа внутренних усилий равна работе внешних сил

I

и = Jrdx = 0, (8)

0

где I - длина оболочки.

Условие (8) с учетом (4) и (7) позволяет получить формулу для критического давления:

2 1 3 1

12(Г-Йд» / (х)^2(х) dx + I S(x)[ф"(x)]2dx

| =-0-1-0-. (9)

/ /(х)ф2 (х) dx 0

Для оболочек средней длины, у которых число волн при потере устойчивости п > 4, в формуле (9)можно пренебречь единицей по сравнению с п2. Из условия минимума р0:

ф0 = 0, (10)

дп2 следует

4 УЗ/3/2 (11)

Р0 =о т-' (11)

где интегралы, входящие в (11), имеют вид

I I I

/1 = У 12(1 -3|х2)д3 (х) dx, 12 = I Е6(х)[ф"(х)]2dx, 1з = I /(х)ф2(х) dx.

0 0 0

При использовании (9) и (10) функции продольной координаты ф(х) должны задаваться и удовлетворять граничным условиям оболочки.

В настоящей статье представлены результаты расчетов по формуле (11) для четырех схем закрепления краев оболочки. Функции продольной координаты ф(х)двух видов представлены в табл.1. Численный эксперимент, выполненный авторами, показал, что для однородных оболочек при равномерном давлении результаты расчетов при рассмотренных видах функций ф(х) практически не отличаются, если удовлетворяются граничные условия задачи. В данной работе функции ф(х) задавались в виде

Граничные условия

Таблица 1

Номер схемы Граничные условия Виды функций ф(ж) Коэффициент а

1 ф(0) = 0 ф'' (0) = 0 ф(1) = 0 ф'' (1) = 0 A sin п

x 2lx I l x

2 ф(0) = 0 ф' (0) = 0 ф(1) = 0 ф' (1) = 0 A sin2 ^p 1, 5п

x 2lx I l x

3 ф(0) = 0 ф' (0) = 0 ф'' (1) = 0 ф'''(1) = 0 A (1 - cos nx) 0, 6п

x4 - 4lx3 + 6l2x2

4 ф(0) = 0 ф' (0) = 0 ф(1) = 0 ф'' (1) = 0 2x4 - 5lx3 + 3l2x2 1, 25п

полиномов. В этих случаях получаем формулу для расчета критического давления в форме [7]. При коэффициенте Пуассона ^ = 0,3 имеем

п2 = 2, 39а ■

R

т

Ро = 0, 293aER (R

2,5

(12)

где а — коэффициент критичности, учитывающий граничные условия. При граничных условиях Новье из (12) следует формула П.Ф. Папковича (а = п).

Рассмотрена устойчивость оболочки с линейно изменяющейся толщиной

д(х) = ¿о(1 - вХ), в = . (13)

I до

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

¡ ч (Л хч Ро

P(x) = Ро(1 - 7т), Y =-

1 Ро

Pi

(14)

В формулах (13) и (14) ви 7 — показатели изменяемости толщины оболочки и давления, д0 и р0 — толщина оболочки и давление в сечении х=0, д и р — эти же параметры на краю оболочки х = I.

Подстановка выражений д(х)и р(х)вида (13) и (14) в формулы (9) и (11) учетом (10) дают зависимости для расчета числа волн при потере устойчивости и величины критического давления в виде аналогичном (12):

Я (до

n2 = n0 f (в), Ро = 0, 293аЕ

1 V R

2,5

kg

(15)

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

Таблица 2

Выражения корректирующих функций

Номер схемы f (в) ks kp

1 J1 - 0, 5в V <Мв) <1 (в) = 1 - 1, 5в + 0, 85в2 - 0,174в2 4(1 - 0, 5в)<3(в) 1 1 - 0, 5Y

2 J1 - 0, 5в V <2 (в) <2 (в) = 1 - 1, 5в + 0, 818в2 - 0, 159в2 4(1 - 0, 5в)<3(в) 1 1 - 0, 5Y

3 J1 - 0,167в V ^з(в) <3(в) = 1 - 2, 415в + 2, 01в2 - 0, 573в2 4(1 - 0,167в)<3 (в) 1 1 - 0, 802Y

4 J1 - 0, 417в V <4 (в) <£>4 (в) = 1 - 0, 698в + 1, 04в2 - 0, 226в2 4(1 - 0, 417в)<3 (в) 1 1 - 0, 565y

Зависимости к$ = /(в) и кр = /(7) при линейном изменении толщины и давления для расчета удобнее представлять в виде функция от соотношения толщин и давлений на краях оболочки:

ß =1 -

öi

Y = 1 -

Po

Графические зависимости корректирующих функций к$ и кр от соотношения толщин и давления приведены на рис. 1 и 2, где номера кривых соответствуют номерам схем граничных условий в табл. 1 и 2.

Kp

4 чУ ч3

1-2

0,8 0,6 0,4 0,2

<"4

1 v

0

0,2 0,4

0,6

0,8 P,/p

Рис. 1. Зависимости корректирующих функций изменения толщины от соотношения толщины при разных граничных условиях оболочки

0 0,2 0,4 0,6 0,8 8[/80

Рис. 2. Зависимости корректирующих функций изменения давления от соотношения давлений при разных граничных условиях оболочки

Для однородной оболочки (¿о = ¿1) в = 0 при равномерном давлении (р0 = Р1) 7 = 0 значения всех корректирующих функций становятся равными единице. Из формул (15) как частный случай следуют известные зависимости [6-10] для однородных оболочек при равномерном давлении.

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

Математическая модель потери устойчивости оболочки, построенная на прямом методе, позволяет по формулам (9)-(11) находить критические нагрузи практически для любых законов изменения толщины и давления вдоль оболочки.

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

1. Элмрос Б.О. Выпучивание цилиндрической оболочки, подверженной действию неравномерного внешнего давления // Прикладная механика: Тр. Амер. общества инженеров - механиков. М.: Изд-во иностр. лит., 1962. № 4. С. 27-31.

2. Моссаковский В.И., Андреев Л.В., Зюзин В.А. Некоторые вопросы устойчивости цилиндрической оболочки под действием неравномерного давления // V Всесоюз. конф. по теории пластин и оболочек: Тез. докл. М.: Наука, 1965. С. 47-51.

3. Кан С.Н., Тодчук В.А. Устойчивость цилиндрической оболочки при действии переменной радиальной нагрузки // Избранные главы по строительной механике оболочек. Харьков, 1965. Вып. V. С. 182-187.

4. Антоненко Э.В., Хлопцева Н.С. Устойчивость цилиндрических оболочек переменной толщины при дей-

ствии наружного давления // Тр. XXI Междунар. конф. по теории оболочек и пластин. Саратов: Изд-во Сарат. ун-та, 2005. С. 7-13.

5. Актуальные проблемы механики оболочек: Тр. междунар. конф. Казань: Новое знание, 2000. 488 с.

6. Товстик П.Е. Устойчивость тонких оболочек. М.: Наука, 1995. 308 с.

7. Антоненко Э.В. Свободные колебания и устойчивость оболочек с упругими краевыми ребрами // Прикладная механика. Киев, 1975. Т. XI, вып. 6. С. 44-50.

8. Вольмир А.С. Устойчивость деформируемых систем. М.: Наука, 1974. 640 с.

9. Кан С.Н. Строительная механика оболочек М.: Машиностроение, 1966. 508 с.

10. Прочность. Устойчивость. Колебания: В 3 т. М.: Машиностроение, 1968. Т. 3. 568 с.

4

3