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

книги из ГПНТБ / Амиро И.Я. Ребристые цилиндрические оболочки

.pdf
Скачиваний:
14
Добавлен:
24.10.2023
Размер:
11.82 Mб
Скачать

 

 

 

2 blk

(Д4тЛ)т

Д3т^0т)] •

Qzlm

q 7.0m

T + i 'oi

L

lk,m

=; —

 

 

m

1 1 — V

 

£)

 

 

 

 

lk,m .

 

 

 

 

 

 

 

 

^m

/i2f 2

«2 \

 

 

tt ,m — Dlk,m

(U

— vd ),

 

 

 

 

 

г . * - =

ф

^ =

T

5 r - | *

J (2 + v)<4l,

Ift.m

Dlk,m

 

 

Jlk,m

 

 

 

1 + V

2

 

 

 

 

= "ВДГ (’ +

“’M

D“ .- =

(1 ~

V „ - ( < £ + W .

? * » = 4

' г 5.

-( 1 - v V

q ~ ~ w ~ q-

После подстановки (IV.2) в (IV. 1) получаем1

^ lk .m ~ q z0т {^Om^Oi 1 -j- Öq; ^ l b .m (A 2m^0m Д 3 т ^ 0 т ) "Ь

+

*Л т (Д4 п /от~ Д3т^0т)1 ] '

 

^ lk .m ~

^ iflzOm l^fft.m

( Д 2;гДот Д 3 т ^ 0 т ) "I-

 

+

* '1» .» (А ,Л « - А 3„Ч„)1.

(IV.3>

^1к,т ^ zGrn Orrfi01

2 blh

'■v*lk,mI (' Д 2mr J ’ Om ~~ A 3mf

Ofm) +

] 4 - 6

01

 

 

1

 

 

“I" F'lk.m (Д4т^0т

Д3т^0т)))‘

 

Анализируя общие члены рядов (IV. 1), представленные в виде (IV.3), можно сделать следующие общие выводы.

Соотношения Аш,т, Вш,т, Cik,m можно представить в виде суммы двух слагаемых, одно из которых описывает деформацию неподкрепленной оболочки под действием заданной нагрузки, второе ■— учи­ тывает влияние ребер на оболочку. Сходимость рядов по т, общие члены которых равны первым слагаемым указанных сумм, такаяlI1*

1 Д2 т, Аз™, Дчт определяются по формулам, приведенным на стр. 97, куда вместо dim нужно подставить d m.

Ill

же, как сходимость рядов, описывающих поведение неподкрепленной оболочки. Сходимость рядов по I (при заданном т), общие члены которых равны вторым слагаемым указанных сумм, такая же, как сходимость рядов по I, описывающих поведение неподкрепленной оболочки под действием периодической нагрузки, распределенной вдоль поверхности контакта обшивки и ребер; если ширина ребра

не равна нулю, то ряды для и сходятся при больших I как —t

ряды для V — как 1 , а ряды для w — как 1 ; при ширине ребра,

равной нулю, ряды для и, v, w сходятся соответственно как

, JJ-; ряды по I, определяющие нормальные усилия и изгибаю­

щие моменты в обшивке и ребрах, сходятся как

; ряды

по

I,

описывающие крутящие моменты, сходятся как

-і- . Ряды

по

т,

общие члены которых определяются вторыми слагаемыми указан­ ных выше сумм, сходятся не медленнее, чем соответствующие ряды для неподкрепленной оболочки, находящейся под действием перио­ дической в окружном направлении нагрузки; сходимость рядов, в которые раскладывается эта нагрузка по т, будет такой же, как сходимость рядов, составленных из A2m<g0m—A3mFem и Д4mF0m—

—A3mü0ro, а эти ряды сходятся не медленнее, чем ряды, описываю­ щие заданную нагрузку; следовательно, ряды по т сходятся не медленнее, чем сходятся ряды по т для неподкрепленной оболоч­ ки, нагруженной периодической в окружном направлении и посто­ янной в меридиональном направлении нагрузкой.

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

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

будет зависеть от - и k. Чем больше у и /г, тем меньше то число

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

112

Действительно, характер изменения напряженно-деформирован­ ного состояния в окружном направлении определяется скоростью

убывания е т' при удалении от ребра і— вещественная часть корня характеристического уравнения для панели [48], заключен­ ной между соседними ребрами и подверженной действию нагрузки,

которая

изменяется в меридиональном

направлении по закону

. яягЕ.

и зависит как от

к

так и от

т

т. е. от числа полуволн

sin -Z-2)

 

 

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

Рассмотрим теперь сходимость рядов по I на числовом примере. Будем исследовать сходимость рядов по /, ограничиваясь в решении (IV. 1) одним членом ряда по т {т = 1), и сопоставим результаты с данными, полученными при табулировании точного решения за­ дачи, представленного в виде одинарного тригонометрического ряда. Рассмотрим оболочку с такими параметрами: k — 8, ѵ = 0,3, ус=

= 6с = Лс = 4,6 а2Л 2 = 1, Дс = 0, у = 2gQ. Результаты вы­ числений амплитудных значений прогибов и изгибающих моментов

при т = 1 представлены в табл.

17. Пользуясь табл.

17 (где "Ѳ= &Ѳ),

вычисляем:

 

 

 

 

 

г г (1

— -V2) .

6

 

w = w* —

— - q sin jtg,

 

j =

qr2а k2G\ sin я |,

 

(IV.'4)

M22 =

qr2a2k2G‘sin я£,

 

 

Анализируя данные табл. 17, можно сделать следующие вы­ воды.

При увеличении числа членов ряда по I решение, полученное в виде ряда, стремится к точному решению; для получения трех верных знаков при определении прогибов можно ограничиться пятью членами ряда по I (Zmax= 4), но уже два члена ряда по I вполне удовлетворительно описывают прогибы оболочки; для опре­ деления изгибающих моментов в поперечных сечениях обшивки

(Gl) с погрешностью до 10% необходимо удержать не менее девяти членов ряда по /; такую же погрешность при определении моментов

в продольных сечениях обшивки (Gl) можно получить, удерживая не менее 17 членов ряда; погрешность решения в рядах становится

8—39

113

Т а б л и ц а 17

 

 

 

о»* -10—1

 

 

 

 

 

«

 

 

 

 

 

 

С2

 

 

Решения

 

 

 

 

 

 

 

 

 

С1

 

 

 

 

 

 

 

 

 

 

_

я

Зл

 

0=0

 

1

 

_

Зя

 

 

_

я

 

_ Зя

 

 

Ѳ=0

ѳ=-

Ѳ =Я

 

я "

Ѳ=я

0=0

- Я

Ѳ—я

 

ѳ

= т

; 4

 

 

 

 

 

 

 

 

 

 

Ѳ=2

ѳ

= -

 

 

ѳ

= г

ѳ = т

Ѳ=т

 

Точное

0,1309

0,1354

0,1414

0,1444

0,1452

0,1281

0,2065

 

.0,2340

0,2336

0,2301

—0,1850

0,0550

0,1189

0,1035

0,0882

В виде ряда

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

ІІП Н Х -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,1331

0,1351

0,1398

0,1446

0,1466

0,1849

0,1939

 

0,2156

0,2373

0,2463

—0,0060

0,0147

0,0647

0,1147

0,1354

4

0,1311

0,1355

0,1413

0,1443

0,1451

0,1480

0,2066

 

0,2324

0,2355

0,2280

—0,1199

0,0549

0,1134

0,1097

0,0813

8

0,1310

0,1355

0,1414

0,1444

0,1451

0,1390

0,2064

 

0,2334

0,2330

0,2294

—0,1493

0,0544

0,1168

0,1013

0,0861

16

0,1310

0,1355

0,1414

0,1444

0,1451

0,1341

0,2065

 

0,2338

0,2334

0,2299

—0,1650

0,0547

0,1180

0,1027

0,0875

32

0,1310

0,1355

0,1414

0,1444

0,1451

0,1316

0,2066

 

0,2339

0,2335

0,2306

—0,1741

0,0550

0,1184

0,1031

0,0879

40

0,1310

0,1355

0,1414

0,1444

0,1451

0,1310

0,2066

 

0,2339

0,2335

0,2300

—0,1758

0,0550

0,1184

0,1031

0,0880

несущественной (меньше 5%), если при выполнении вычислений удерживать не менее 41 члена ряда.

При определении напряженно-деформированного состояния реб­ ристых оболочек мы ограничивались 41 членом ряда по I. При этом принималось во внимание, что рассматриваются оболочки, усилен­ ные достаточно большим числом ребер {(fk4 > 1), а это, как ука­ зывалось, должно привести к улучшению сходимости рядов по I. Учитывалось также, что с возрастанием т коэффициенты раз­ ложения нагрузки в тригонометрические ряды уменьшаются и это ослабляет потерю точности расчета в результате ухудшения схо­ димости рядов по /.

Исследование влияния числа удерживаемых членов в двойных тригонометрических рядах (И.4) на величины характеристик напря­ женно-деформированного состояния ребристой оболочки выполнено также на примере рассмотренной выше (§ 2 главы III) оболочки с перекрестной системой ребер \ подверженной действию локальной нагрузки, ребра которой размещены симметрично относительно сре­ динной поверхности обшивки.

Характеристики напряженно-деформированного состояния, вы-

£

численные в точке £ = -у-, Ѳ = 0 в седьмом приближении при

различном числе удерживаемых членов, приведены в табл. 18, где т!— число сохраненных членов в рядах по |, а п' — по Ѳ. Нетруд­ но заметить, что при определении прогибов даже сравнительно небольшое число членов ряда дает хорошую точность. Так, по­ правка к значению прогиба, вычисленному при сохранении m' X X n' = 5 X 5 членов ряда при дальнейшем увеличении числа

1 Здесь использованы результаты, полученные В. Н. Хитровым.

членов оказывается меньше 9%. Для изгибающих моментов полу­ чение достоверных результатов требует сохранения значительно большего числа членов ряда. В табл. 19 представлены величины

Т а б л и ц а 18

Прогибы, усилия.

 

 

 

 

m 'X n '

 

 

 

 

моменты

5X5

11X11

15X15

21X21

31X31

41X41

51X51

75X75

100X100

 

шЕг-102

26453

28639

28963

29064

29121

29136

29142

29144

29145

Pz

 

 

 

 

 

 

 

 

 

Мг-10б

1083

1350

1421

1466

1509

1526

1537

1549

1551

рг

 

 

 

 

 

 

 

 

 

Ми-105

1241

1700

1866

1953

2040

2083

2108

2135

2147

р2

 

 

 

 

 

 

 

 

 

NyrAW

14087

14213

14185

14179

14189

14191

14195

14197

14199

Pz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 19

 

Относи-

 

 

 

m'X n'

 

 

 

 

тельная по­

 

 

 

 

 

 

 

 

 

грешность

5X5

11X11

15X15

21X21

31X31

41X41

51X51

75X75

 

моментов

 

ш х

30,17

12,95

8,38

5,48

2,71

1,61

0,90

0,13

 

т у

42,19

20,83

13,09

9,04

4,98

2,98

1.81

0,55

 

114

8’

115

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

§ 2. Продольно подкрепленные оболочки

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

с =

1,

т|* = ІО2,

£' =

££2 =

5л,

V = 0,3,

(IV.5)

 

 

 

 

 

 

 

 

 

 

 

 

8С=

ѵс =

0,

Ас =

0,05;

 

 

 

с =

6,75,

т)' =

2,32 • ІО2,

£ ' = 3 0 ,

ѵ =

0,3,

(ІѴ.6)

 

 

 

 

 

 

 

 

 

 

k \ =

6,69,

ѵс =

0,238,

Дс =

0,005.

 

 

Более подробно изучена оболочка, имеющая параметры (ІѴ.5).

Для нее построены эпюры перемещений п и ш и

эпюры изгибающих

моментов в обшивке М х и М у (рис. 12—15). В табл.

20

приведены

величины усилий N x и Ny в обшивке и усилий Л/с и моментов Мс в ребрах. В табл. 21 и на рис. 16— 19 приведены данные, иллюстриру-

рующие изменение перемещений п и

w, изгибающих моментов М х,

М у в оболочке при удалении от ребра.

 

 

 

 

В табл. 20, 21

и на рис.

12—19 приняты такие обозначения:

 

( 1 - ѵ Ѵ

*

(1

V1 г

 

 

(1 — Va) r2

t

и =

 

 

 

 

2.

. 2

 

 

 

 

 

Ehk

qir,

V = — Ehk

■qv*, w =

 

ЕІх

q W '

NK=qrT\,

Ny = qrT’v

N xy = ~ N

 

1 — V qrS*,

 

 

 

 

 

 

 

V х

 

2

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 20

 

 

 

 

*

 

 

*

 

 

 

 

 

 

 

10 T1

 

 

 

 

*

 

 

Б/Ь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тс

-О с

 

 

 

Ѳ'=0

Ѳ'=л

Ѳ'=0

 

Ѳ'=л

 

 

 

 

 

0,05

0,3404

—0,2760

0,3196

0,3702

0,1426

0,0096

 

0,10

0,4945

—0,4330

0,5678

0,6234

0,2414

0,0192

 

0,15

0,4730

—0,4386

0,7330

0,7670

0,2934

0,0226

 

0,20

0,3484

—0,3388

0,8330

0,8432

0,3128

0,0264

 

0,25

0,1966

—0,2002

0,8904

0,8856

0,3150

0,0290

 

0,30

0,0620

—0,0702

0,9236

0,9126

0,3112

0,0308

 

0,35

—0,0382

0,0296

0,9434

0,9320

0,3068

0,0320

 

0,40

—0,1030

0,0952

0,9554

0,9458

0,3036

0,0328

 

0,45

—0,1384

0,1312

0,9620

0,9542

0,3020

0,0330

 

0,50

0,1493

0,1426

0,9642

0,9572

0,1516

0,0332

116

В =П

О

0,10

0,20

0,30

0,40 £ /£ 2

 

 

Рис.

12.

Рис. 13.

Рис. 16.

Рис. 19.

Mx = - q r W G ] ,

My = - q r V k % ,

Mxy = Myx = qrWk2H*,

EI

E J s

( l - v 2)q rc, Ѳ' = kQ.

M : ^ ^ ( \ - v

2) q k % Nc = ^ -

 

Eh

 

 

 

(IV.7)

Анализируя данные, приведенные в табл. 20, 21 и на рис. 12— 19, нетрудно заметить, что наиболее важные (с точки зрения ра­ счета) пиковые значения характеристик напряженно-деформиро- ванного состояния обшивки лежат в ее сечениях, проведенных че­ рез Ѳ' = 0 и Ѳ' = я. Можно выбрать также две характерные точки по меридиональной координате: одной из них обязательно должна

быть £ =

0,5

£2, вторая должна лежать у края

оболочки.

В даль­

нейшем

в качестве второй

выбиралась либо

£ =

0,625

£2, либо

É = 0,1

£2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 21

 

 

 

1=0 5 6,

 

6=0.1 Ь

 

 

 

«V

—ЮГ)

*

*

*

—и*

*

 

 

л

- г 2

10

- га

— I0S]

 

 

0

0,1493

0,9642

0,4965

0,5677

1,9031

0

 

 

1

0,1488

0,9642

0,4945

0,5678

1,9028

0,0056

 

 

2

0,1473

0,9642

0,4887

0,5681

1,9020

0,0112

 

 

3

0,1449

0,9641

0,4790

0,5686

1,9001

0,0168

 

 

4

0,1414

0,9641

0,4657

0,5692

1,8988

0,0224

 

 

5

0,1370

0,9640

0,4490

0,5700

1,8965

0,0281

 

 

6

0,1318

0,9639

0,4290

0,5710

1,8937

0,0337

 

 

7

0,1257

0,9637

0,4061

0,5721

1,8904

0,0394

 

 

8

0,1188

0,9636

0,3805

0,5734

1,8866

0,0450

 

 

9

0,1111

0,9634

0,3525

0,5749

1,8825

0,0505

 

 

10

0,1028

0,9632

0,3224

0,5764

1,8779

0,0559

 

 

И

0,0938

0,9630

0,2905

0,5781

1,8731

0,0611

 

 

12

0,0843

0,9628

0,2571

0,5799

1,8679

0,0660

 

 

13

0,0743

0,9626

0,2226

0,5819

1,8624

0,0707

 

 

14

0,0639

0,9623

0,1871

0,5838

1,8567

0,0750

 

 

15

0,0531

0,9621

0,1510

0,5859

1,8568

0,0789

 

 

16

0,0421

0,9618

0,1145

0,5880

1,8448

0,0824

 

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

119

Анализируя кривые, приведенные на рис. 16—19, нетрудно за­ метить, что V, до, Мх, Му в рассмотренном интервале изменения Ѳ— монотонные функции 1 Ѳ, при этом изгибающие моменты на ребре максимальны, прогибы на ребре минимальны. Однако по­ следний вывод не является общим. На рис. 20 приведены эпюры прогибов для оболочки с параметрами (IV.6). Нетрудно заметить, что в сечении £ = 0,5 прогиб обшивки на ребре больше прогиба

/fff'

О

2

4

6

іб £ '

 

 

Рис.

20.

 

обшивки между ребрами. Этот факт можно объяснить тем, что про­ гиб обшивки в достаточно длинной тонкой оболочке должен быть близким к прогибу бесконечно длинной оболочки, а прогиб в центре жесткого ребра может оказаться больше прогиба бесконечно длин­ ной оболочки. Эпюры изгибающих моментов Ми для первой (IV.5) и второй (IV.6) оболочек аналогичны (см. рис. 15 и 21).

Отметим, что данные табл. 20 и 21 и рис. 12 — 21 достаточно убедительно свидетельствуют о необходимости учета дискретного размещения ребер при определении напряженно-деформированного состояния подкрепленных оболочек.

Было изучено влияние параметров ребристой оболочки на про­ гибы до и изгибающие моменты в обшивке. Введение координат

=

kl,

Ѳ' = АѲ позволило записать формулы для определения

и*,

ѵ*,

до*, Т\, Т*т S*, G*, G*, Я*, Т* и G* в виде, где указанные

безразмерные перемещения и усилия зависят только от шести па­ раметров (ѵ, с, т|*, 6*, ус. Q- На числовых примерах исследова­ лось влияние изменения параметров ребер (ц*, б*, ус) при фикси­

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

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

1 Сглаженные вершины кривых для Мх и Му — следствие учета ширины

7L

ребер (напомним, что угловая ширина ребра в приведенном расчете —).

120

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

женно-деформированного состояния, кроме и* и Т\. При увеличе­

нии 6Сот 0 до 10 они изменяются не более чем на 10% (для сравнения использовались максимальные по абсолютной величине значения безразмерных функций). При этом вблизи ребра указанное влияние более существенно (точки / и II), чем между ребрами (точки III

и IV).

Безразмерные функции, как правило, зависят от 5Смонотонно.

Так, с увеличением бё Т\ увеличивается, Gl уменьшается, ТІ на ребре (Ѳ' = 0) уменьшается, а между ребрами (Ѳ ' = я) увеличи­

вается и т. д. Влияние бс более существенно у края оболочки-

Т а б л и ц а 22

*

 

Чс

0,3285

0

5

0,2556

10

0,2380

15

0,2271

20

0,2196

25

0,2140

30

0,2097

35

0,2062

40

0,2035

45

0,2011

50

0,1992

*

ѵс=о

Чс

5

0,0082

10

0,0199

15

0,0279

20

0,0337

25

0,0380

30

0,0414

35

0,0440

40

0,0462

45

50

0,0496

 

 

 

*

 

П

 

*2

 

 

і и

IV

0,3285

1,1206

1,1206

0,2840

1,0536

1,0317

0,2685

0,9806

0,9608

0,2586

0,9310

0,9131

0,2517

0,8955

0,8791

0,2465

0,8688

0,8536

0,2425

0,8481

0,8338

0,2393

0,8315

0,8179

0,2368

0,8180

0,8050

0,2346

0,2327

0,7971

0,7851

*

 

 

Т а б л и*ц а 23

 

 

 

*2

П

О

со

ѴС=0

ѵс=о,з

 

II?

 

1,0536

1,0321

—0,1108

—0,0953

0,9806

0,9604

—0,0847

0,9310

0,9117

—0,0771

0,8955

0,8767

—0,0714

0,8688

0,8565

—0,0670

0,8481

0,8301

—0,0634

0,8315

0,8138

—0,0605

0,8180

0,8005

—0,0581

0,7893

—0,0560

0,7971

0,7799

І2І

Соседние файлы в папке книги из ГПНТБ