Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методическое пособие 787

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
12.92 Mб
Скачать

b h3

h

2

 

 

I 2

1

b h

 

1 0.25

 

 

(4.1)

 

 

12

1

 

2

 

 

 

 

 

 

 

 

 

 

 

Момент сопротивления составного сечения равен:

2 I

 

 

 

 

W

 

 

 

(4.2)

 

 

 

 

 

 

H 0.25

 

 

5. Проверка принятых размеров поперечного сечения арки по I и II предельным состояниям:

I.

 

N

 

Mдеф.

R

(5.1)

F

Kw W

 

 

 

сж

 

.

где N - продольная сила; F - площадь поперечного сечения; Mдеф. - изгибающий момент от

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

Mдеф.

M g M N

,

(5.2)

 

 

 

 

 

 

 

 

 

 

 

H h

 

 

 

 

 

 

где: M

g

M

0,25l

, M

N

N e N

T

 

, h

- высота опорной диафрагмы h = h

0.2

h ;

 

 

 

 

 

2

T

T

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

e- эксцентриситет, e H hT ; - коэффициент, учитывающий деформацию приращения,

2

определяемый по формуле:

 

=1

 

2 N

 

 

(5.3)

 

3000 R

 

F

 

 

 

 

 

 

с.

бр.

 

II:

fn

 

fц.

fu

 

1

l ,

(5.4)

 

kI

200

 

 

 

 

 

 

 

 

где fц. - прогиб как цельного элемента сечения равен:

fц.

 

M H l2

,

(5.5)

 

 

1

10

E I

 

 

 

где l1 2l ; M H -максимальный изгибающий момент от нормативной нагрузки; kI - коэффи-

циент, учитывающий влияние податливости связей на момент инерции /1/; fu - предельный нормативный прогиб [2].

6. Определение сдвигающих сил на полудлине плоскости сечения арки: 1) от поперечной нагрузки [1]:

T

 

M g S

,

(6.1)

 

g

I

 

 

 

 

 

 

 

 

где S - статический момент, равный:

 

 

 

 

 

S b h h1

 

(6.2)

 

 

1

2

 

 

2) от изгиба моментом:

80

T

 

M N S

(6.3)

 

M N

I

 

 

 

 

7. Определение количества связей на полудлине арки:

n

k

 

Tg TM N

,

(7.1)

T

c

T

 

 

 

 

 

 

НП

 

 

где kT - коэффициент, учитывающий неравномерность распределения сдвигающих усилий, kT 1; TНП - расчетная несущая способность нагельной пластины, определяемая по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TНП Tна m. nН ,

(7.2)

где Tна - расчетная несущая способность на один срез нагеля, [1]:Tна

= 1,4 (кН); m. - коли-

чество «срезов», m. = 1; nн - количество нагелей в пластине, nн = 6 шт.

Количество связей на длине арке определяется:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nc 2 nc 1.

(7.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1. Трехшарнирная арка из деревянных брусьев на нагельных пластинах: а) схема арки; б) сечение арки; в) нагельная пластина (НП)

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

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

1.СНиП II-25-83. Деревянные конструкции. Нормы проектирования. - М.: Стройиз-

дат, 1983.

2.СНиП 2.01.07-85* Нагрузки и воздействия. Госстройиздат. - М., 2003.

References

1.SNIP II – 25 – 83. Structural timber. Desing standards. – M.: Stroyizdut, 1983

2.SNIP 2.01.07 –85. External forces and effect. State building publishing house. Gosstroyizdut – M.:, 2003.

81

УДК 624.071.2

 

Воронежский государственный

Voronezh State University of Architecture

архитектурно-строительный университет

and Civil Engineering

Д-р техн. наук, проф. С.Н. Булатов

Dr. tehn. sciences professor S.N. Bulatov

Россия, г. Воронеж, тел. 8(4732)71-52-02

Russia, Voronezh, ph. 8(4732)71-52-02

С.Н. Булатов

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

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

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

S.N. Bulatov

TO THE CALCULATION ON STABILITY OF CYLINDRICAL SHELLS

OF DIFFICULT STRUCTURE

Within the framework of technical theory of thin shells a power method is get the exact decision of task on determination of stability of circular cylindrical shell, supported stringer and structural ring of regular structure. Inflexibility of shell is assumed variable along formative. In a maximum transition calculation formulas coincide with known from literature.

Keywords: barrel shell, stability, rigidity.

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

Врамках теории [1] ниже рассмотрена цилиндрическая оболочка, подкреплённая стрингерами и шпангоутами регулярной структуры. Такие оболочки принято называть кон- структивно-ортотропными. Область их применения в основном – авиационная и ракетнокосмическая техника.

Всоответствии с [2] задачу об устойчивости этой оболочки, нагруженной в плоскости симметрии погонной кольцевой нагрузкой интенсивности q, решим аналитически. В целях повышения несущей способности конструкции, в плоскости действия нагрузки q предполагается наличие силового шпангоута с изгибной жёсткостью EI. Здесь и ниже обозначения и

рисунок приняты из [2], но толщина h1 срединной усиленной секции собственно оболочки на длине -ℓ≤x≤ℓ, предполагается симметричной относительно срединной поверхности, что дает возможность исключить ранее принятую гипотезу о допустимости её разрыва вследствие

©Булатов С.Н., 2009

82

перепада толщины h2 в концевых и центральной секциях. Допущения об отсутствии сдвигов в срединной поверхности и нерастяжимости контура, принимаемые в технической теории оболочек, остаются в силе. Следуя [1], радиальные прогибы W представим через функциональную неизвестную Ψ (x), угловую координату φ и число волн в окружном направлении

 

 

 

 

 

 

 

 

W (x)cosn .

 

 

 

 

 

(1)

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

 

 

 

 

 

 

 

 

 

V 1

(x)sin n

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

осевые U перемещения вдоль оси x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

R

d (x) cosn

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

 

dx

 

 

 

 

 

 

 

 

Тогда относительные линейные 1 ,

1

и угловые ,

x

деформации определяются

выражениями

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 n2

 

 

 

 

 

 

 

 

 

 

 

 

(x)cosn ,

 

H (x)cosn ,

 

 

 

H (n

2

1)R

2

(x)cos n ,

 

Hkp (n

2

1)(nR)

1

 

(4)

 

 

 

 

 

 

 

(x)sin n .

Статические силовые факторы с учётом (4) принимают вид

 

 

1 ERn

2

 

 

 

 

 

 

2

h(x)n

3

 

 

 

 

 

 

 

(x)cosn ,

 

 

N ER

 

(x)sin n ,

 

M1 D(x)

(x)cos n ,

 

 

 

 

 

 

M D(x)(n2

1)R 2 (x)cos n .

 

 

(5)

Потенциальная энергия деформации

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1dx 2 2dy 0 ,

 

 

 

(6)

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

где Гi – энергия на единицу длины оболочки толщины h1(i=1) и концевых секций толщиной h2(i=2), П0 – энергия силового шпангоута, расположенного в плоскости действия погонной нагрузки q.Раскрывая Гi в соответствии с [2], получим

 

Di (n

2

1)

2

2

 

 

2

 

 

2

 

 

 

R

 

 

i

 

 

 

ER

hi d

 

i

, (i=1,2)

(7)

R

4

 

i

n

4

 

dx

2

 

2

 

 

 

 

 

 

 

 

 

 

 

где Di-цилиндрическая жёсткость соответствующих секций оболочки, Ψii(x), hi – приведенная толщина, R – радиус срединной поверхности, Е – модуль Юнга. Потенциальная энергия деформации шпангоута и внутренних сил оболочки в плоскости её симметрии х=0 равна

 

 

 

EI

(n2 1)2

2

q

 

n2 1 2

 

ER2hi d 3 1

 

 

R .

(8)

0

 

4

kp

 

 

 

 

2R

1

 

2R

1

 

n

4

dx

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

Из условия минимума потенциальной энергии деформации (6)– вариационного уравнения Эйлера-Пуассона

i

 

d

 

i

 

 

 

d

2

 

 

 

 

 

 

 

 

 

 

 

 

i

0 , (i=1,2)

(9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

i

 

 

 

 

 

dx

 

 

 

 

 

 

dx i

 

 

i

 

 

вытекают два однородных дифференциальных уравнения четвертого порядка с постоянными коэффициентами

d 4

 

 

 

0

, (i=1,2)

(10)

i 4K

dx4

 

 

i

i

 

 

 

 

где

 

 

 

 

 

 

 

 

K

 

 

D n4 (n2

1)

2

.

(11)

i

 

i

 

 

 

4ER6h

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

83

 

 

 

 

Из характеристического уравнения: r4

4K 4

0 ,

(i=1,2)

i

i

 

 

вытекают корни

 

 

 

ri Ki 2 i0 , (i=1,2);

ri

Ki

2 i0 , (i=3,4),

где i0 –мнимая единица.

 

 

 

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

1

Ñ1 cosk1x chk1x C2 sin k1x shk1x C3 cosk1x shk1x C4 sin k1x chk1x ,

(12)

2

Ñ5 cosk2 x chk2 x C6 sin k2 x shk2 x C7 cosk2 x shk2 x C8 sin k2 x chk2 x ,

(13)

Выполняя условия сопряжения двух концевых секций оболочек с центральной

 

 

1 2 ,

1 2 ,

h1 1 h2 2 ,

h1 1 h2 2 ,

(14)

приходим к системе четырёх алгебраических уравнений, решение которой возможно совместно с уравнениями, вытекающими из граничных условий при x=0. Из (14) следует, что радиальные прогибы, углы поворота, нормальные и касательные усилия в сечениях x=±ℓ, равны между собой.

В плоскости симметрии конструкции при x=0 отсутствует депланация сечения, или

d 1

K1(C3 Ñ4 ) 0 ,

Ñ3 Ñ4

(15)

dx x 0

Вторым граничным условием в плоскости действия погонной нагрузки q является уравнение смешанной вариационной задачи

 

1

 

d

 

 

 

 

0

 

 

 

 

 

 

1

 

 

0 ,

(16)

 

 

 

 

 

 

 

1

 

dx

1

 

1 x 0

 

 

что после подстановки входящих величин, даёт дополнительно алгебраическое уравнение вида

C (EI

n2 1

q

)

n2 1

C

8ER2h K 3

 

(17)

 

R3

 

 

 

R3

 

i

1 0 .

1

 

 

 

 

kp

 

 

 

3

n4

 

 

 

Определитель системы четырёх уравнений, вытекающих из (14), равен

=-2, то есть

система имеет одно решение. Тогда постоянные интегрирования C5 ,C6 ,C7 ,C8

определяются

из формул Крамера

1

 

 

 

 

2

 

 

 

 

 

3 ,

 

4 .

 

C

,

 

C

,

 

C

C

(18)

5

 

 

 

 

6

 

 

 

7

 

 

 

8

 

 

В концевых сечениях оболочки x=±L отсутствуют радиальные прогибы и изгибающие

моменты, что равносильно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

x L 0 ,

 

x L 0 .

 

 

(19)

 

 

 

 

 

 

 

Удовлетворяя (19), нетрудно получить ещё два алгебраических уравнения и с уравнением (17) имеем систему трёх однородных алгебраических уравнений. Её определитель равен нулю в случае нейтрального равновесия оболочки, то есть, когда форма поперечных сечений слегка изогнута. Опуская сравнительно простые, но громоздкие выкладки, получим

искомую формулу для определения величины критической нагрузки qkp , приложенной в плоскости силового шпангоута при x=±L

q

 

EI

n2 1

 

8ER3h K 3

A,

(20)

kp

R3

i 1

 

 

 

n4 (n2 1)

 

 

где A - коэффициент, учитывающий геометрию оболочки и условия её опирания в концевых сечениях. Первое слагаемое правой части – нагрузка, воспринимаемая силовым шпангоутом, второе – собственно оболочкой.

84

Следует отметить, что при любых граничных концевых условиях, структура формулы (20) сохраняется; изменяется лишь содержание коэффициента A. Так, например, для жёстко заделанной по концам оболочки граничные условия имеют вид

2

 

x L 0 ,

 

 

x L 0 ,

(21)

 

 

что равносильно отсутствию радиальных прогибов и углов поворота концевых сечений. Для оболочки со свободными краями концевые граничные условия соответствуют ра-

венствам

2

 

x L 0 ,

2

x L 0 ,

(22)

 

то есть нормальные и касательные усилия должны быть равны нулю.

Достоверность полученных результатов подтверждается предельным переходом к оболочкам постоянной толщины известным из литературы [1]:

 

q

kp

EI

n2 1

 

8ER3hK 3

A, K

4

 

Dn4 (n2 1)2

,

 

(23)

 

 

R3

 

n4 (n2 1)

 

 

4ER6h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где для оболочки со свободно опёртыми, заделанными и свободными краями

 

 

A

2(sh2 KL cos2 KL)

,

 

A

 

 

shKL

sin 2KL

 

 

,

 

A

shKL sin 2KL

.

sh2KL sin 2KL

 

 

 

2(sh2 KL sin2 KL)

 

2(ch2 KL cos2

KL)

 

 

 

 

 

 

 

 

 

 

 

 

 

Для длинной оболочки n=2, A=1. Тогда из (20) имеем формулу,

 

 

 

 

 

 

 

q

kp

 

EI

 

3

 

 

0,38

 

Eh2

 

h

,

 

 

 

 

 

 

 

 

 

R3

(1 2 )0,75

R

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полученную Л.Н. Балабухом и В.М. Марченко [3] на основе интегрального уравнения Фредгольма с симметричным ядром (EI=0)

q

 

12

 

Eh

h 1,5

 

 

 

 

 

.

9

 

(1 2 )0,75

 

kp

 

 

 

R

 

Численные расчёты, выполненные на ЭВМ, показывают, что коэффициент A в значи-

тельной

степени

зависит от безразмерного параметра

K2 L , относительной толщины

h h1 h2

и длины

L усиленного пояса. Для оболочки с жёстко защемлёнными и сво-

бодными краями коэффициент A →1 при K2 L →∞. В то время, когда K2 L →0, величина ко-

эффициента A →∞. Это факт свидетельствует о том, что концевые диафрагмы сближаются с

уменьшением длины оболочки и дляK2 L 0 совпадают с силовым шпангоутом, поэтому

q

.

 

 

 

kp

Для оболочки со свободными краями коэффициент A →0 при K2 L →0, так как в этом

случае концевые диафрагмы отсутствуют. Выявлена также зависимость A=A( h ) при фикси-

рованных значениях безразмерных параметров =0,2;0,6

и K2L 0,6;0,9;1,8 Безразмерная

толщина изменялась в пределах 0≤h ≤2,05. Анализ графических зависимостей показывает,

что для защемлённых и шарнирно-опёртых оболочек при прочных равных условиях, коэффициент A →1 с увеличением h . Для оболочек со свободными краями при отсутствии диафрагм 0<A≤1. Графики A=A( ) и K2L 0,6;0,9;1,8 при h =1,9 показывают рост A с увеличением параметра и для случая жёсткой заделки A= 4,25 при =0,5 и K2 L 0,6 .

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

85

зультатов опирается на предельный переход к формулам, известным из литературы и сравнительный анализ численных расчётов.

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

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

2.Булатов С.Н. Устойчивость ступенчатой цилиндрической оболочки под действием радиальных давлений и осевых усилий. – М.: Строительная механика и расчет сооружений.

М., 1970, №3.

3.Балабух Л.Н., Марченко В.М. Устойчивость тонкой оболочки, нагруженной местными усилиями. Труды ЦАГИ, 1952.

References

1.Caen S.N. Building mechanics of shells - Moscow: Machinery construction (Mashinostroenie), 1966. - 508p.

2.Bulatov S.N. Stability of stepped barrel shell under radial pressure and axial thrusts. - Moscow: Building mechanics and analysis of structure, 1970, №3

3.Balabukh L.N., Marchenko V.M. Stability of thin shell loaded with local load. - Researches of TsAGI, 1952

УДК 624.071.2

 

Воронежский государственный

Voronezh State University of Architecture

архитектурно-строительный университет

and Civil Engineering

Д-р техн. наук, проф. С.Н. Булатов

Dr. tehn. sciences professor S.N. Bulatov

Россия, г. Воронеж, тел. 8(4732)71-52-02

Russia, Voronezh, ph. 8(4732)71-52-02

С.Н. Булатов

ОБ ОДНОМ ПОДХОДЕ К ИССЛЕДОВАНИЮ ДЛИННЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК ПЕРЕМЕННОЙ ЖЁСТКОСТИ

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

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

S.N. Bulatov

ABOUT ONE GOING NEAR RESEARCH OF CYLINDRICAL SHELLS

OF VARIABLE INFLEXIBILITY

Stability long structurally orthotropic a cylindrical cover of a variable thickness is investigated at local loading by lateral pressure. Simple enough settlement dependences for definition of bearing ability of a thin-walled design are received.

Keywords: barrel shells, side pressure, load-carrying ability.

В работе предложен один из возможных подходов к определению критической радиальной нагрузки, приложенной к боковой поверхности конструктивно-ортотропной оболоч-

© Булатов С.Н., 2009

86

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

На рис. 1 приведена расчётная схема и геометрические размеры оболочки, нагруженной равномерно распределённым по поясу толщины h1 и длинны 2 , радиальным давлением интенсивности P. При P=Pкр конструкция разрушается – теряет устойчивость с образованием регулярно расположенных по окружности волн. Задача по определению Pкр решается энергетическим методом.

Рис. 1. Расчетная схема оболочки, нагруженной равномерно распределенным по поясу радиальным давлением

Следуя [1], потенциальную энергию деформации длинной оболочки запишем в виде

2 1

dx 2 2 dx 0,

(1)

0

 

 

где потенциальная энергия деформации первой Г1 и Г2 второй секций на единицу длины со-

ответственно равна

D1 n2

1 2

 

2

 

ER2 h1

 

" 2

 

n2

1

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PKP

 

R,

(2)

2R4

 

2n4

2R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ER h2

 

 

 

 

 

 

 

 

 

D1 n

 

1

2

 

 

" 2 R.

 

 

(3)

 

 

 

 

 

 

2n4

 

 

 

 

2

 

2R4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь и ниже обозначения приняты по [2]. Функциональная неизвестная (x) не

терпит разрыв в силу симметричного расположения h1 относительно срединной поверхности оболочки (рис.1). В плоскости симметрии хоу предполагается отсутствие силового шпангоута (EI=0) и депланации контура поперечного сечения. Так что угол поворота при х=0

d

 

x 0

0.

(4)

 

dx

 

 

 

 

Условию (4) удовлетворяет функциональная неизвестная

 

fe kx sin kx cos kx ,

(5)

где f - амплитудное значение функции ,k -

неизвестный коэффициент затухания.

Про-

дольные U, окружные V и радиальные W перемещения на краях длинной оболочки равны нулю. Поэтому Pкр не зависит от краевых граничных условий. С учётом (2) ÷ (5) выражение потенциальной энергии деформации (1) получает вид

6 1 t

4 C 1 t

2

 

2

1

P

0 ,

(6)

 

 

 

 

1

1

 

3

 

 

 

 

 

 

 

 

1

êð

 

 

 

откуда искомая критическая нагрузка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

 

 

 

в 1 t

 

 

 

 

C 1 t

 

 

 

 

1

 

,

(7)

кр

 

1

1

 

2

2

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

87

 

D1 n2 1

 

 

ER2 h1 K 4

 

 

 

 

1

 

 

 

1

2k

2

cos 2k sin 2k ,

 

 

 

 

 

 

 

 

 

 

где в

R3

,

C

 

,

t1 1,

t2

 

 

1

, 1

 

3 e

 

n4 n2 1

 

 

 

 

 

 

2 e 2k cos 2k sin 2k 2 .

Для изотропной оболочки переменной толщины t1 t2 t

P

(1 t ) 4 C

(1 t

2

)

(1 ) .

(8)

kp

 

 

 

1

 

3

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда, как частные случаи, вытекают формулы для изотропной оболочки постоян-

ной толщины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D(n2 1

 

 

4 ER3hK

4

 

 

Pkp

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

(1 1) .

(9)

 

 

R

 

 

3 n

4

(n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

 

С увеличением ширины нагружаемого пояса ( L) параметр

k , а 1 0 , и

формула (9) принимает вид

 

D n2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pкр

 

 

4 ER3hK 4

 

 

 

 

R3

 

 

 

 

3

n4

n2

1

.

(10)

Для длинных оболочек n=2, поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

3D

 

ER3hK 4

 

 

 

(11)

 

 

 

 

R3

 

 

 

36

 

 

 

 

 

 

 

кр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Минимальное значение Pкр имеет место при k=0. Тогда из (11) получим известную

формулу Брайана [4] для нагруженной оболочки по всей поверхности

 

 

 

Pкр

 

 

 

 

Eh3

 

 

 

 

 

 

(12)

 

 

4 1 2 R3

 

 

 

Рассмотрим случай нагружения длинной оболочки погонной кольцевой нагрузкой q в плоскости симметрии xoy при отсутствии силового шпангоута (рис. 2).

Рис. 2. Длинная оболочка, нагруженная погонной кольцевой нагрузкой

Условие равенства работ внутренних и внешних сил для оболочки со ступенчатым изменением толщины записывается в виде

L

2 1

dx 2 2 dx 0

(13)

0

 

 

где первые два слагаемых – потенциальная энергия деформации нагруженного кольцевого пояса ширины 2 и концевых секций оболочки, T-работа внешних сил для координаты x=0:

 

D

2

 

2

hi

 

 

 

 

 

i n2

1

2

ER

" 2

R

, (i=1,2)

(14)

2R4

 

 

i

 

 

2n4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

88

 

 

 

 

 

T

1

qкр R 0 H 0 Rd qкр

n2 1

 

2

0 R .

 

 

 

(15)

 

 

 

 

 

2

2R

 

 

 

 

 

 

С учётом (14) и (15) уравнение (13) получает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

L

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Di n

 

 

1

 

 

 

2

2

 

 

2ER

h1

 

 

2

 

 

 

 

2

 

 

 

2

 

qкр

 

 

 

 

dx h

 

 

 

 

 

 

 

"

dx h "

 

 

 

0 (16)

 

R

3

 

 

dx

n

4

n

2

 

 

 

 

dx

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

Функцию деформации после потери устойчивости оболочки, удовлетворяющую граничному условию (4), примем в форме (5). Выполняя интегрирование в (16), получим расчётную формулу

q

 

в 1 t

C 1 t

 

 

,

в 1.5D n2

1 KR3 ,

C 2ER3

 

K 3

n4 n2 1 , (17)

кр

2

h

 

1

1

 

 

2

1

 

 

 

 

 

 

 

 

 

 

h1

 

h2 ,

i

и ti – коэффициенты;

(i=1,2).

 

 

 

 

h

 

 

Из формулы (17) вытекает ряд частных случаев. Например, для изотропной оболочки постоянной толщины имеем расчётную формулу

 

3 D1 n2 1

 

2ER3

 

K 3

 

 

 

h

 

 

qкр 2 KR3

 

 

n4 n2 1

.

(18)

 

Коэффициент затухания K найдём

из

условия минимума

критической нагруз-

ки qк

K 0

 

 

 

 

 

 

 

 

K=[Dn4(n2-1)/4ER6h]0.5 .

(19)

Подстановка (19) в (18) с учётом n=2 для длинных оболочек приводится к известной [3] формуле

q

 

 

12

 

Eh

h 1.5

 

 

 

 

 

 

 

 

.

(20)

кр

9

1 2 0.75

 

 

 

 

 

R

 

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

Решим теперь задачу устойчивости оболочки, представленной на рис.3, когда толщина оболочки h=h(x) изменяется вдоль образующей по линейному закону. Условие равенства работ представимы в виде

 

 

2 dx T 0

(21)

0

 

89