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

книги / Оболочки и пластины

..pdf
Скачиваний:
21
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

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

*»«., = j y f [2 V 2 + e - * { V l - X +

V l + x )] . (4,24,19)

и потому коэффициент динамичности

(4,24,20)

Коэффициент динамичности кольцевого давления оказывается зна­ чительно больше коэффициента динамичности внутреннего давления. Например, при х=0,45 имеем /Сд5' = 1,35, /Сду= 1,035; для оболочки с от­ ношением h/R«0,2 это соответствует скорости движения нагрузки #«1200 м/сек.

§ 25. О РАСПРОСТРАНЕНИИ УПРУГИХ ВОЛН В ОБОЛОЧКЕ

Представляют интерес исследования распространения возмущений в тонкостенных конструкциях в связи с проблемой их прочности и жест­ кости [52]. Рассмотрим распространение упругих волн в оболочках [53].

Вначале выпишем систему дифференциальных уравнений [54] для осесимметричных колебаний круговой цилиндрической оболочки тол­ щины 2h и радиуса срединной поверхности г0, которые затем будут использованы при изучении распространения упругих волн в оболочке.

Исходя из общих уравнений теории упругости без каких-либо ги­ потез о характере деформирования, основываясь на алгоритме Н. А. Кильчевского [55], путем исключения из матричных операторов всех функций перемещений, кроме одной или нескольких, И. Т. Селезов построил обобщенные дифференциальные уравнения (с точностью до членов порядка | 3) колебаний [56]:

(4,25,1)

- { ( 2 + 1 4 ,1 -6 % - ~ + 1%

д3 I

Р\ — Рг

+

 

d t * 3 /

2 г

 

+ {[5'. + Р«

- 1*'.

д х *г3 + П

^d4t*2d xг}* I-4

2^- +

 

+ {[6 +

£•/.1 — Ъ31»

SP

д2

1 Pi +

Р2

(4,25,2)

+ S3/»

 

 

 

где х* — переменная,

отсчитываемая вдоль

оси

оболочки; ®о и «о —

радиальное и осевое перемещения; q\$.— радиальная нагрузка на внут-

тренней и внешней поверхностях оболочки;

/?*,ч— осевая

нагрузка на

тех же поверхностях; а3 и Ь, — постоянные,

зависящие от

пуассонова

числа 1/v. Здесь приняты такие безразмерные параметры:

 

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

Исследуем бесконечную цилиндрическую оболочку, к которой на расстоянии ** = 0 приложен сосредоточенный осесимметричный импульс

Qa(t*), где o(t*) — функция Хевисайда.

Задачу будем

решать, исклю­

чая точки приложения сосредоточенного

импульса. Тогда

приходим

к дифференциальным уравнениям вида:

 

 

 

|[(--Мо + а7^в) + ё2 (---аФо +

Мо + а8^в + аФ^\

дх*2----Ь

+ [ао + I2 (ai + b2aQ)\ —^

\о.фй

Ь3а0+ Д768 + афв]

— |-

+ I— ьо+ ¥■(— Ь0аз —<h— Mo + bt + а.Ья + а106в)1

 

+

+ I 2 [^0а 5 +

а 4 4" 6.1

 

 

 

 

d t*2d x **

 

 

 

+ I2 [fle -г Ьъ\

шо =

0;

- Ь0 +

 

+ [1 + £ 2Ьй

 

+ t*b3

(>-

д х *2

1

d t*2

дх**>

1 Метод!нуждается в строгом математическом обосновании-

-

е’6‘ ^

г

+

5‘6* - £ г }“: +

 

+ { [ - » . -

6 4 1

• £

6 4 - £ г -

6 4 - 4 ^ ) № = 0 .

Эти уравнения получены из (4,25,1) и (4,25,2), причем члены порядка выше | 3 отброшены.

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

Ы

дЦ

_____ 1_

 

 

дх*3

д х *3 - о

~

5* Х

2

 

(4,25,5)

 

«о 1+о =

“о |—о =

0.

 

 

 

 

 

 

Три условия удовлетворяются на бесконечности. Здесь

24 (2А. + 3|i)

Х =

X-ф- 2(1 '

где X и р, — постоянные

Ламе. В момент

t*=0,

если

принять

нулевые

начальные условия, выполняются такие начальные условия, также вы­ текающие из точной постановки задачи:

дЩ d t*

du0 d t*

0 1

d t*3

_ d3u0

d t* 2

_

~

1

d3w'0

d t*3

0 *

dt*3

= 0

(4,25,6)

= 0 .

Решение уравнений (4,25,3) и (4,25,4) в пространстве лапласовых изображений, удовлетворяющее условиям сопряжения (4,25,5), усло­ виям на бесконечности л начальным условиям (4,25,6), имеет вид

2 |2 w

( х .

ч

у

 

4 fee x p ( - n fe|** |)

 

%Q

°

£

2

pU r i - 4 ) A r + (nl-4)A*\

 

 

 

 

(A j

А г) exp (— n3j x * | )

 

(4,25,7)

 

 

pnз 1(4 — n23) A i + (n j n22) A«] ’

 

 

 

%LU(X'

p) = —

V

lkAk^ ( - n k\x*\)

 

XQ

 

 

 

2 Pgk [(«? — n2)

+ («2 — гф Л2]

 

 

 

_j_

<з (4i -fr A2) exp (— n31x* | )

(4,25,8)

 

 

Pga [(«i — Л3) 4 , + ( 4 — «2) A d ’

 

 

 

где n h ( k = 1,2,3)

корни

характеристического

уравнения Re«h>0,

W0(x',p) = Cje ->’‘,w ;(x',C)d r,

0

U0 ( x ' , p )

=

f e - P

‘- u ; ( x ' , t')d t'

 

 

 

О

 

 

преобразования Лапласа, A\=gi(l2g3hs^), А^

получается из A i кру­

говой перестановкой индексов:

 

 

 

 

I t =

1 2Ь вп *+

[ ( Ь в +

1 % +) ТО>2],

g. = i *b3ni - [(й„ -

Т О

+

Тор 2]"2* +1(1 +

Vbt) р2 + Тор 4]-

Т1о теореме обращения [57] получим решение задачи в виде контурных интегралов:

212

. (x .

n =

у

_ L1 _

Г

Akexp (pt*

nfenft1х*1x* II))

dp

xQ

°

'

k = f 2

2 Я /

L

pnk K n i ~ rt3> Ai +

( 4 — nl ) A il

 

—4

(Ax Лг) exp (pt*— я31** I)

dp;

(4,25,9)

 

i

J

J

 

/

рл3 [(n2 — n2) .4t — (ra| — n j) 42]

 

 

_2£

 

 

lkAkexp (pt*— ftfe 1x* I)

dp

xQ ■ w

. f l -

k^Sf 2 lj Pgk [(«1 - 4) Ax+ (n\ -

n\) Аг]

 

+

_ i _

Г

l3(Ax + A2)exp(pt*-na\x*\) dp

(4,25,10)

 

2ni

}

Рёз [{n\ — /ф Ax + (4 — n\) A2]

 

 

 

 

§ 26. О РАСПРОСТРАНЕНИИ УПРУГО-ПЛАСТИЧЕСКОЙ ВОЛНЫ НАГРУЖЕНИЯ В ОБОЛОЧКЕ

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

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

Пусть на цилиндрическую свободную (незакрепленную) оболочку длины L в момент t = 0 начинает действовать внешнее давление, сим­ метричное относительно оси вращения и распространяющееся по поверх­ ности оболочки с некоторой скоростью v, которая может быть меньше или больше скорости распространения упругих волн растяжения и из-

гиба в оболочке -V-р (1 — V2) .Пусть также ширина нагрузки £г^оо,

т. е. нагрузка, или сходит в момент te = l/V с оболочки, или остается на ней. При V<k\o деформации в оболочке распространяются со скоростью fei, а при V>ki деформированной окажется вся оболочка позади фронта

нагрузки, т. е. деформации будут распространяться со скоростью V Следовательно, при У<Аю надлежит решать вторую смешанную задачу, а при V>k\ — третью смешанную задачу.

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

оболочки поворачиваться и перемещаться в радиальном направле­ нии, т. е.

<р = 0

при * =

0,

х = L,

 

(4,26,1)

w = 0

при .г =

0,

х — L.

 

(4,26,2)

В осевом направлении концы оболочки могут иметь

жесткую

(непод­

вижную) заделку

 

 

 

 

 

ы=0

 

 

(4, 26, 3)

или скользящую (подвижную) заделку

 

 

 

 

7

^ = 0

 

 

(4,

26, 4)

Кратко воспроизведем ход решения

задачи для

различных

скоро­

стей движения нагрузки. В случае дозвуковой скорости движения на­

грузки V<k\o (рис.

4.15). На линии x = k i0t имеем

нулевые условия:

<р = ш= и= 0, (fx = (pt =

w x = W t — u x = u l = 0 . Дойдя до

противоположного

края оболочки x = L, упругая волна x = kl0t отразится от него. Отражен­

ная волна

пойдет

в обратном направлении

по характеристике

4 х = fix' dt,

отделяя

область бегущей волны (1)

от области отражен­

ной волны (рис. 4.15); начальная скорость отраженной волны будет равна, конечно, k\0.

Решая вторую смешанную задачу, определим сначала решение в области (1), ограниченной прямыми x = k\0t, х = 0 и характеристикой

dx = k{\ ) dt, включая границы этой области; подобно определим и об­

ласть (2), ограниченную прямой x = L, характеристикой dx kS? dt, на которой функции и, w, ср и их первые производные непрерывны, и ха­

рактеристикой dx = k®) dt, исходящей из точки встречи характеристики

dx = k\2)dt с прямой х = 0.

В случае сверхзвуковой скорости движения нагрузки v>k\0 (рис. 4.16). Здесь решение в области, ограниченной характеристиками первого семейства d x = ± k xdt, исходящими из точки О (0,0), (линия ОМ)

л точки (линия LM), не зависит от граничных условий на

концах оболочки, а чтобы найти его, надо решить задачу Коши по дан­ ным на отрезке прямой x = W ^ 0 < * < L , в треугольнике,

ограниченном этой прямой и характеристиками ОМ и LM (рис. 4.16). Из анализа уравнений движения цилиндрической оболочки видно,

что при нагрузке постоянной интенсивности решение задачи Коши вдоль прямых* параллельных фронту движения нагрузки x = vt, не ме­ няется. Это значительно уменьшает трудоемкость вычислений при ре­ шении задачи Коши, ибо потребуется определять лишь точки, лежащие непосредственно на характеристике ОМ.

Для начала интегрирования возьмем на прямой x = vt точку Оi, достаточно близкую к точке О. Теперь, решая задачу Коши по данным на отрезке ООь определим точку Мх на характеристике ОМх положи­ тельного направления первого семейства. В точке Nu лежащей на пере­ сечении характеристики первого семейства положительного направле­ ния, исходящей из точки 0\, с отрезком M\N\ прямой, параллельной линии x = vt, решение будет то же, что и в точке М\. Поэтому следую­ щую точку М2 на характеристике ОМ, а также все остальные точки на

ней можно определить или общим методом, решая задачу Гурса по данным на характеристиках М\0\ и 0\N U или, воспользовавшись по­ стоянством решения вдоль отрезка M\N\, решить задачу Коши по дан­ ным на отрезке M\NX. В обоих случаях основные исходные точки i и / одинаковы; вспомогательные точки, лежащие на характеристиках вто-

v

рого и третьего семейств, разные. Точкой i служит предыдущая искомая точка (т). Точкой I при решении задачи Гурса служит старая точка /. Новая точка ] определится из уравнений

■Xj — Xt = Ащ (*,— */),

по формулам

Vtj k jiti + x i — Xj

V - k Ui

(4,26,5)

x i = xi + V (tj — ii).

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

При определении точки т из решения задачи Коши координаты точки / при данном интервале Ax=Xjх{ определяются по формулам

Xj = х + Ал:,

(4,26,6)

tj= tt + — .

1 1 у

Вдоль характеристики LM решение определяется как и на характе­ ристике ОМ, причем решение на характеристике LM равно решению

на характеристике ОМ, но со смещением во времени на tL = — , т. е. для любой другой функции f будем иметь

/ьл1(0 = /ом ^ ---^ .

В частном случае мгновенного приложения равномерного давления по всей длине оболочки решение в углу между характеристиками ОМ

и LM й осью х в каждый данный момент времени t=t*

одинаково во

всех точках отрезка прямой

t = t*, заключенного между

характеристи­

ками ОМ и LM.

же типа для конической оболочки. Если

Рассмотрим задачу того

на свободную коническую оболочку действует внешняя

нагрузка, то

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

После того как возмущения дойдут до края оболочки x = L, она начнет перемещаться в направлении оси оболочки как абсолютно твер­ дое тело, имея, конечно, одновременно и перемещения вследствие де­ формации. Следовательно, чтобы иметь возможность применять урав­ нения к оболочке в. ее движении относительно центра массы, необхо­ димо остановить центр массы оболочки, для чего нужно к центру массы приложить силу инерции Даламбера, равную осевой (в рассматривае­ мом случае — горизонтальной) составляющей равнодействующей внеш­ ней нагрузки. Сила Даламбера действует на каждый элемент массы оболочки. Поэтому ускорение центра массы (перегрузка) равно

(4,26,7)

м

где Рх — горизонтальная составляющая равнодействующей внешних на­ грузок; М — масса, например, всей конструкции. На оболочку (в от­ дельности) приходится часть силы Даламбера, равная

Р' = дМ*

(4,26,8)

Здесь М* — масса оболочки. Если с оболочкой не связаны никакие со­ средоточенные массы, находящиеся внутри, то интенсивность этой на­ грузки можно записать в виде

р*

(4,26,9)

р* = S ’

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

диусом г0);

тогда для новой переменной граничные

условия при х = 0

и X — L для

конической оболочки будут иметь тот же

вид, что и у ци­

линдрической.

При решении задачи Коши для определения решения на характе­ ристике ОМ (рис. 4.16) мы уже не сможем в случае постоянной нагруз­ ки определить все величины на этой характеристике, не зная их значе­ ний во всем треугольнике OML, из-за того что решение вдоль прямых, параллельных линии x = vt, не будет неизменным, так как коэффициенты в уравнениях движения теперь будут зависеть от S(rQ=5sinp).

Однако для конической оболочки малой конусности радиус попе­ речного сечения г0 мало меняется по длине, и потому, без большой по-

грешности, для простоты расчета, когда длина оболочки не очень ве­

лика, этим изменением можно пренебречь.

Можно считать, что в узких полосах шириной Ал;', прилегающих к характеристикам ОМ и LM, решение вдоль отрезков, параллельных прямой x = S —sQ= vt, меняться не будет даже при большой конусности.

Таким образом, в случае конической оболочки, по которой бежит равномерная нагрузка со скоростью v>k\0i решение на характеристи­ ках ОМ и LM можно найти с большой степенью точности, не определяя его во всем треугольнике OML. Причем, в отличие от цилиндрической оболочки, в самом треугольнике OML решение вдоль прямых, парал­ лельных OL, меняется тем сильнее, чем больше угол конусности обо­ лочки и чем длиннее сама оболочка. Заметим также, что при малой конусности оболочки, когда на нее действует лишь внешнее давление, можно заменить его, без большой погрешности (существенно упрощает расчет), нагрузкой, направленной перпендикулярно оси конуса (не вво­ дить силы Даламбера).

Приведем пример числового расчета усеченной конической стальной (сталь 3) обо­

лочки с углом полураствора Р = 11°, для

^ 136. Расчет произведен на внешнее дав­

ление р ,х= \ \ 2 к г!см 2 и касательную нагрузку р 8 = 22 кг/см 2, распространяющиеся вдоль оболочки со скоростью \0 км /сек в сторону возрастания го полосой, ширины /= 200 /г,,

причем касательная нагрузка направлена тоже в сторону подрастания /*0. При расчете предполагалось, что материал оболочки не сжимаем.

Решение найдено

на

характеристике ОМ .

Предполагалось, что на

расстоянии

Да'= 0,85 h по направлению прямой x = vt состояние оболочки не

меняется. Результаты

расчета приведены на рис. 4.17 и в табл. 4.3.

 

 

 

Оказалось, что при

/ =

h

 

х=л'в=11,1

Л впервые

i ? s = l l i l ----------------------- в сечении

 

 

р (1 —

V2)

 

 

в пролете оболочки возникают пластические деформации. Предполагалось, что на конце

Jt=0 осуществлена

скользящая заделка.

Результаты расчета даны на том же рис. 4.1?

и в табл. 4.4. При

/0 = 0 ,5 6 5 ------ ^ -

— в заделке деформации были еще упругие..

Р (1 — V2)

Произведенные -расчеты показывают, что деформации поперечного сдвига по сво­

ей 1величине примерно равна

деформация

изгиба

(удлинения— сжатия

во

внешних

слоях);

деформация

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

меньше

деформаций

изгиба и поперечного сдвига. Скорость

распространения

упруго­

пластической волны

нагружения, -возникшей

в точке -(**, te)

характеристики

ОМ,

в мо­

мент t = t B

будет примерно

равна

скорости

движения нагрузки V.

Вдоль

характери­

стики L M

решение находится аналогично тому, как и на ОМ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

4, 3

 

X

 

 

 

е°

 

10*/t

 

10*

io* -^2L

 

. 0 = ^

точек

 

~

± Y

— E—

ю*__fi

 

 

sz

 

 

h У

P ( 1 -

V =)

**

 

 

 

es

 

 

 

es

 

 

 

 

 

 

 

 

 

 

 

 

0,0

 

0

 

0

 

0

 

 

0

0

0

 

 

0

1,1

 

0,85

 

0,85

—0,113

 

 

0,359

0,621

0,312

 

 

0,457

2,2

 

1,70

 

1,70

—0,231

 

0,267

1,90

1,10

 

 

1,17

3,3

 

2,55

 

2,55

-0,357

 

1,05

2,47

1,74

 

 

1,44

4,4

 

3,40

 

3,40

—0,491

 

3*,05

3,23

2,94

 

 

2,25

5,5

 

4,25

 

4,25

—0,629

 

5,92

4,28

4,73

 

 

3,74

6,6

 

5,10

—5,10

—0,785

 

9,98

5,74

7,30

 

 

6,03

7,7

 

5,95

 

5,95

—0,926

 

— 15,7

7,79

10,9

 

 

9,43

8,8

 

6,80

 

6,80

—1,06

 

— 23,9

10,7

16,0

 

14,3

9,9

 

7,65

 

7,65

-1 ,1 9

 

— 35,4

14,8

23,1

 

21,3

0,10

 

8,50

 

8,50

—1,30

 

— 51,7

20,6

33,2

 

31,2

1,11

 

9,35

 

9,35

—1,37

 

— 74,17

28,8

.47,4

 

45,2

2,12

 

10,20

10,20

—1,40

 

—107,0

40,4

67,6

 

65,1

3,13

 

11,14

11,14

—1,35

 

—159,0

62,1

100,0

 

97.2

4,14

 

11,60

11,58

—1,71

 

—200,0

109

134,0

 

130,0

5,15

 

12,04

12,07

—5.72

 

—287,0

130,0

174,0

 

167,0

6,16

 

12,93

13,13

-1 2 ,1

 

—441,0

33,3

270,0

 

241,0

 

 

 

 

 

 

 

 

 

 

 

 

Т аб л и ц а

4.4

 

X

‘ V

E

 

10*

 

10=Л —

 

юг _!H_

 

Ш’ - Ъ -

 

h

 

 

10= Ysz

 

точек

 

h V p(l-v)2

 

 

 

 

 

 

 

 

 

 

0.1

 

0

0,142

 

—0,0000382

—0,663

—1,16

0,773

 

0,773

0,2

 

0

0,565

 

—0,000139

—1,08

—1,77

1,20

 

1,20

§ 27. КОЛЕБАНИЯ ПЛАСТИНЫ, ОБУСЛОВЛЕННЫЕ ТЕПЛОВЫМ УДАРОМ

Рассмотрим шарнирно

опертую по всему контуру пластину

к поверхности z = h/2 которой, начиная с момента / = б,

внезапно подводится тепло,

причем количество тепла, подводимое за

единицу времени к единице поверхности, равно постоянной величине q. Нижняя поверхность пластины z = —Л/2 и края пластины х = ± . а , у = ± Ь предполагают идеально теплоизолированными.

Распределение температуры при этом зависит только от коорди­ наты z и времени t, следовательно, определяется уравнением

 

 

д*Т

_

дТ

 

 

(4,27,1)

 

ат дг2

 

dt

 

 

при начальном условии

 

 

 

 

 

 

 

 

 

Т == 0 при

t =

0

 

 

(4,27,2)

и граничных условиях

 

 

 

 

 

 

 

 

 

дТ

= q

при

 

 

h

(4,27,3)

X - дг

z =

 

2

дТ

0 при z

=

h

 

(4,27,4)

дг

- =

2

*

 

 

 

 

 

 

 

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

ном изменении температуры ее поверхности приведено в работе [59].

этой

Здесь мы остановится на варианте А. Д. Коваленко [60] изложения

задачи.

 

 

 

 

 

 

 

 

 

Решение уравнения (4,27,1) при условиях (4,27,2) легко может быть

найдено при помощи преобразования Лапласа и имеет вид [61]

 

Г

ЯЬ

Г

,

12г * + 1 2 h z - h ? ____ 2

_ у

( - 1 ) " с- п,р,

Я*

/ 2 2

Л '

 

X

[

я 2

'

24И?

я 2

ЛшЛ

п?

2

[ h

) \ '

 

 

 

 

 

 

 

 

 

 

 

(4,27,5)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о =

атДг

 

 

 

 

 

 

 

 

 

Н

 

Л*

'

 

 

 

Используя уравнение температурного изгиба пластины

V V » —

+

 

Eab?

v *

T

M=

 

о

 

(4,27,6)

 

 

 

( l - v ) D

 

 

 

 

 

 

 

и заменяя в соответствии

с принципом

Даламбера

интенсивность по-

 

 

 

 

1

d*w

, получаем уравнение дви­

верхностных сил силами инерции — рН -

 

жения

 

 

 

 

 

 

 

 

 

 

 

 

у 2у 2w +

Рл

d2w

+

Eah2

 

V

T M

=

O

,

(4,27,7)

 

D

 

(1 -

v) D

 

2

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h_

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

/i2

j

T (z) z d z .

 

 

 

 

 

(4,27,8)

 

 

 

 

 

 

 

 

 

 

 

 

h