Научная статья на тему 'Задача о продольной устойчивости подкреплённой стрингерами оболочки в разномодульной упругой среде'

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

CC BY
75
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УСТОЙЧИВОСТЬ / ШАРНИРНО-ОПЁРТАЯ ЦИЛИНДРИЧЕ­СКАЯ ОБОЛОЧКА / КОМБИНИРОВАННЫЙ АЛГОРИТМ ПЕРЕБОРА ВАРИАНТОВ

Аннотация научной статьи по математике, автор научной работы — Михайловский Евгений Ильич, Кораблев Анатолий Юрьевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Михайловский Евгений Ильич, Кораблев Анатолий Юрьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

Вестник Сыктывкарского университета. Сер. 1. Вып. 17.2013

УДК 539.3

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

Е.И. Михайловский. А.Ю. Кораблев

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

1. Введение

В докторской диссертации [1] (форм. (25)} построена так называемая деформационная теория ребристых оболочек, главная особенность которой заключается в том. что в ней впервые наряду с реактивной силой учтён реактивный момент от ребра жёсткости. Названная теория подробно изложена в монографии [2]. В частности, уравнения статики конструктивно ортотропной цилиндрической оболочки, получаемые путём "размазывания" регулярной системы стрингеров имеют вид:

с0Ги = —R2q + -

V

Куд^ Кг В? <Э£4

¿Иги

В? д^дф2 Кп д4и2

Кг д3и2 В? д£2д(р

(1.1)

Д2 <Э£4

Здесь первая часть уравнения представляет собой матричную запись уравнений статики цилиндрической оболочки в смещениях; Ки, Кп -

© Михайловский Е.И., Ксрабяев А.Ю., 2013.

жёсткости стрингера при изгибе соответственно в нормальной плоскости и из этой плоскости. - жёсткость при кручении. - жёсткость стрингера при ратяжении (сжатии); I - расстояние между соседними стрингерами по дуге поперечного сечения; £ = х/Я. (р = у/Я.

Примем следующие допущения:

для расчёта оболочек без рёбер допустимо использовать упрощённую теорию Доннела-Власова ([3 , табл. 13.2. форм. (13.18с));

конструктивно ортотропная оболочка испытывает осесимметрич-ную деформацию, т.е. df/d^p = О V/;

ni

стрингеры изгибаются в нормальной плоскости как гибкие стержни. т.е.:

(fui Кп dAU2

ct

de в? de

— 0;

¡у) на оболочку действует лишь нормальная нагрузка, т.е. ql = q^2 = 0. Тогда уравнению (1.1) можно придать вид

щ _ R2 "0"

U2 0

Со

_Qn_

, с0 =

Eh

где

(1.2)

А =

«зз =

• 0

д?

о 1^-Ц о ¿Ю

-р-^ 0 а33

У 1 ь? (1 - S)J\ э4()

\ 12 i?2 IhR2 ) д?^1'-

I - тождественный оператор, J - момент инерции поперечного сечения стержня. И - толщина неподкреплённой оболочки. Е - модуль Юнга, у - коэффициент Пуассона материала ребристой оболочки.

Для получения из системы (1.2) соответствующего уравнения относительно функции прогиба воспользуемся операторным методом [4].

"Детерминант" матрицы А имеет вид

Задача об устойчивости подкреплённой стрингерами оболочки... 29 I . * 1 -ЧУ (1-у*)Лс1\) 1-у ?<1\)

Заменяя в матрице А последний столбец правой частью уравнения (1,2} и вычисляя "детерминант" так составленной матрицы, получим

аеъАт - 2 ^ Применяя теперь формально правило Крамера, можно записать

/ Т А \ 1 А

ии — ——— или (ае£А)ги — ае£Ат аеЬ А

Последнее равенство после элементарных преобразований представляет искомое уравнение в виде

где

4Ь4 = 12(1 — ¿о— Ек3

к2' _ 12(1-1/2) ■

Представим удельную нормальную нагрузку в виде суммы

Цп = дп + Яп (1-4)

В условиях наличия внутри и вне оболочки винклеровых сред различной жёсткости (с^сг). можно записать [5]

(¡п = -С1-Ш+ - с2гу_ (1.5)

Кроме этого при рассмотрении продольной устойчивости цилиндрической оболочки от действия сжимаюших усилий Т0. следует положить

[3]

Подставив соотношения (1.4) в уравнение (1.3) и выполнив замену

= (1-6)

~ 7Г.К ^ 7Г-Й X ТГХ Гп ,

< = ^ = тгте.м

(Ь - длина оболочки)

и. сохранив за новой переменной £ прежнее обозначение окончательно получим

(1 + a)— + A^w + A— + klW+ + k2w_ = 0 (1.7)

di

Здесь

. ттх rn n EJ Л T0L2

4 _ 12(1 - _ k - — P ~ ' ° " 12(1 - i/2)' ** ~ 7r4d0

(1.8)

где L - длина оболочки, TQ - сжимающие усилия. Обратим внимание, что от сжимаюших усилий Т0 зависит параметр Л.

2. Постановка задачи

Рассматриваем дифференциальное уравнение (1.7):

(1 + a)wIV + A/34w + kxw+ + k2w_ = -Xw" (2.9)

и найдём такое Л, при котором имеет нетривиальное решение краевая задача с граничными условиями шарнирного опирания:

w( 0) = ги(тг) = 0; w"{0) = w"(ir) = 0

Проинтегрируем по частям функционал, образованный уравнением (2,9). умноженный на w(£). Имеем

/*7Г Л7Г /*7Г /*7Г Л7Г

(1+а) / и/" d£+4/34 / / w+2d£+A;2 / = Л / го'

JO JO JO JO J 0

(2.10)

Заменим формулу (2.10) приближённой с использованием дискретного представления функции w. задаваемой её значениями на равномерной сетке, т.е. Wi = w(xi),i = 0..п

Далее аппроксимируем производные конечно-разностными схемами:

wi+1 - 1 „ - 2wi + tUj-i —^—-^--(2.П)

Интегралы будем вычислять по квадратурной формуле трапеций

/

J о

Г /(Ode = ^ + /(6) + /(6) + - + Д£п-1) + Щ^) (2.12)

Значения срезок функции в узлах сетки представляем формулами

«>+(&) = Ъг^и ю-(6) = С1 ~ = { о ^ < О (2ЛЗ)

Граничные условия шарнирного опирания аппроксимируем формулами

мо = иоп = 0, и;_1 = — гУ1, го„+1 = -ги„_1 (2.14)

После преобразования уравнение (2,9) примет вид

Ли/ + Сиз = \Q-ib (2.15)

Матрицы А. <5 - пятидиагональные:

5 -4 1

-4 6 -4 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 -4 6 -4 1

1 -4 6 -4

1 -4 6 -

1 -4

3 0 -1

0 2 0 -1

-1 0 2 0 -1

-1 0 2 0 -1

-1 0 2 0

-10 3

С = (Над [кф{ + к2{1 - Ь) + 4/З4] к Необходимое условие минимума функционала имеет вид

Аш + Сю - \<Эгю = 0

3. Комбинированный алгоритм

Для решения задачи будем использовать комбинированный алгоритм. который включает в себя применение на первой стадии алгоритма полного перебора вариантов (ППВ). а на последующих - локального перебора вариантов ЛПВ [6].

Для построения части собственного спектра уравнения применяется алгоритм ППВ форм изгиба, который заключается в следующем:

[) перебираются все 2т~1 возможных представления вектора формы

6 = [61,62, ...,6Ш_1]Т:

и) для каждого варианта вектора формы решается задача на собственные значения;

ш) запоминается собственная пара (число и форма) для которой форма изгиба согласуется с выбранным вектором формы.

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

После этого будем применять алгоритм ЛПВ. используя в качестве приближения полученную качественно-адекватную форму:

,]) последовательно удваиваем число узлов сетки путём деления интервалов пополам;

|]) осуществляем перебор вариантов лишь вблизи корней последнего приближения к искомой собственной форме;

Д]) процесс продолжается до тех пор. пока соответствующее собственное значение не стабилизируется с требуемой точностью.

4. Численное решение

Применим комбинированный алгоритм "ППВ+ЛПВ" (ППВ при п = 6. ЛПВ при п = 24) при следующих значениях параметров: к\ = 16, = 18. Значения параметра а будем изменять в пределах от 0 до 90.

Из таблицы 1 можно заметить, что при прочих равных параметрах если увеличивается а(а — 0 - оболочка не подкреплённая стрингерами), то оболочка укрепляется, т.к. увеличивается значение А, а значит и значение первой критической силы,

Таблица 1

а 4/34 Л

0 43.68 17.54060017038568

10 43.68 61.73360008456593

30 43.68 91.95393728433567

50 43.68 112.1229929869527

70 43.68 132.2920486895459

90 43.68 152.4611043921551

Таблица 2

а 1/200 1/100 к/Я 1/50 1/30 1/20

10 573 254 125 77 62

30 885 433 208 150 93

50 1079 515 291 170 113

70 1274 598 373 190 133

90 1468 681 393 211 153

Зафиксируем значения Ь=200см и г/—0.3 и получим зависимость Л от И/Я при различных а.

Построим соответствующий график при а = 10.

Литература

1. Михайловский Е.И. Деформационная теория ребристых оболочек и её приложения: автореф. дисс.... докт, физ.-мат. наукр. Спеп. 01.02.04 - Мех, деф. тв. тела, Л.: Ленинградский гос. ун-т. 1989. 31с.

2. Новожилов В.В.. Черных К.Ф., Михайловский Е.И. Линейная теория тонких оболочек. Л.: Политехника, 1991, 656 с.

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

4. Филин А.П. Элементы теории оболочек. Л.: Стройиздат. 1987. 384 с.

5. Михайловский Е.И. Элементы конструктивно-нелинейной механики, Сыктывкар: Изд-во Сыктывкарского ун-та. 2011. 212 с,

6. Михайловский Е.И.. Тулубенская Е.В. Алгоритм локального перебора вариантов в одной существенно нелинейной спектральной задаче // ПММ. 2010. Т. 71 Вып. 2. С. 299-S10.

Summary

Mikhailovskii E.I., Korablev A.J. The longitudinal stability of a cylindrical cover supported by stringers in a multimodulus elastic surroundings

In the work the stability of a longitudinal compressed hinge-sup ported cylindrical cover stiffened by stringers and located on the border of two Winkler's ambiences is considered. The problem is solved using a combined exhaustive search algorithm.

Keywords: stability, hinge-supported cylindrical cover, combined exhaustive search algorithm.

Сыктывкарский университет

Поступила 09.06.2013

i Надоели баннеры? Вы всегда можете отключить рекламу.