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

книги / Устойчивость и колебания пластинок и оболочек с отверстиями

..pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
17.03 Mб
Скачать

Общее решение уравнения (7.14) представляется в виде интеграль­ ных уравнений:

•йУ:

г т { - [ / 2 (*) z" ~ Л (л) z lv] -

 

/о (■*)

— 2

I / 2 W ГУ4" /4 W гу] Л/з £Tix

U 2 (*) г) +

/ 4

(*) г)\ ( -

l )y+1

t rJ* ^ (ле

; (7.22)

 

У-1

 

 

 

 

 

 

4

 

 

 

 

 

wl= z l +

^ Ayoz-yerrv+

^

( - 1У+1

o e r^ ^ jxe г^ х ;

 

 

ySt

7 S1

 

о

(7.23)

 

 

 

 

 

 

 

 

 

( - 1 ) * 1-D ГУ ]Х\ **ГЧХС1Х

Здесь A jo— 'произвольные постоянные интегрирования

( / = 1, 2,

3, 4); функция z определяется линейным дифференциальным урав­ нением с постоянными коэффициентами:

aAzlv -\-(izZu -{~aQZ = 0.

(7.24)

Параметр D есть определитель Вандермонда:

 

1

1

1

1

 

 

г х

г2

 

г 3

г4

 

D = _2

2

 

2

2 = const ^-0,

(7. 25)

Гх

Г 2

 

Г 3

Г 4

 

-3

_3

 

3

з

 

Гх

Г 2

 

Гз

Г 4

 

где dj — определитель третьего порядка, получающийся из матри­ цы определителя D путем вычеркивания первой строки и /-го столбца.^

Функция [х определяется зависимостью вида

=

б0 -j-

—[—02^50^^ —f—бз'гед*11 —|—

(7 .2 6 )

гд е

в0= 0 ;

в1=

/ „ ( Л ) - * 1

_ 1 ;

 

 

 

 

а0

 

 

9 _ _ / » (* )/о М ~ / о

{X) /2(х)

0L _ a L .

 

 

/о (•*)

 

Л0

д0

 

 

 

 

 

(7 .2 7 )

a

fi( x ) a 4— f 4{x)a2 . п

fo(*)

1

до /4(*)

/о (*) во

Исследуемую полоску берем по ширине, рав­

р

 

ной расстоянию между центрами вырезов вдоль

 

дуги. В этом случае уравнение

изгиба полоски

 

запишется в форме (7.14). Его -решение ищем по

 

методу последовательных приближений Н. К. Ку­

-й -тз-й -

ликова. Оно записывается в форме (7.22).

 

h It

Чтобы решить нашу задачу в нулевом

приб­

 

лижении,

следует в (7.22) предположить, что

ч Т г т Г П

р = 0 . Так как общее решение должно удовлетво­

 

рять граничным условиям задачи, то, подставляя

Рнс. 7.1

в них (7.22)

и (7.29), .получим:

 

 

 

 

 

 

r i Л ю - { - Г 2Л 20

г з Л 3 о / ' 4 ^ 4 о = 0 ; |

 

 

гie * Лю-Ьггв * Лго-Ь^зС^ Л3о-|-,‘4вг<*Л4о= 0;]3

Л [/2 (0) ~ r \]iA\o-j- Л [/2 (0) — г1| Л2о+

гз [/2 (0) — rl]Лзо+З

 

+ г4 [/2 (0) — Л] Л40= 0 ;

(7.28)

Л[f2 [I)r'\\ Л 10 4~ Л[/2 (/) rl]ёг‘*Л2о + r3j[/2 (7)—r| j er,f Л 30-f-

^4] er*^ 4 o = 0 .

Эта система имеет решение, отличное от нуля лишь в том слу­ чае, если определитель, составленный из коэффициентов при неиз­ вестных Aj[0| равен нулю. Раскрывая его и учитывая значения кор­ ней гь г2, г3 и г4, определяемые соотношениями (7.21), получим два -алгебраических уравнения:

г2 — г ? = 0 ;

 

(7 .29)

| Q—(ri4Ta)/

__Q

(7 .30)

 

 

Рассмотрим второе из них. Оно легко преобразуется, в результате имеем

 

cos {rt-f- г2) / — COS (rj — r 2) l =

0

(7 .3 1 )

или

sin r x/ =

0;

sin r2l = 0 .

 

(7 .3 2 )

Отсюда

r it = n m

m = 1,

2,

3 ,...

 

(7 .3 3 )

 

г21 = л т

m = 1,

2,

3 ,...

 

(7 .3 4 )

Освобождаясь от корней, из зависимости

(7.33)

найдем

 

 

^ 2 = —r

flo"fa2»

 

(7 .3 5 )

 

 

a*

 

 

 

 

 

*где

 

а==щ л//.

 

 

 

 

Если число полуволн т достаточно велико, то можно записать условие минимума коэффициента a2t включающего внеш­ нюю нагрузку Р, приравнивая нулю производную от а2 по а:

-

а 0+ 2 а = 0 ,

 

da

аз

0 1

 

откуда

 

4 —

(7 .36)

 

а = у а0.

Подставляя (7.36) в (7.35), найдем выражение для критической силы

р _ 2 Уао

(7 .37)

«2

Такое же выражение получается из уравнения (7.34) и (7.2Э). Таким образом, мы получили приближенное выражение для критической нагрузки, действующей на полоску, выделенную из оболочки, ослабленной вырезами произвольной формы. Переходя от него к критическому напряжению для перфорированной оболоч­

ки в целом, найдем

Эта формула получается из (7.37) путем деления на Л.

Решения (7.37) и (7.38) представляют собой лишь нулевое при­ ближение окончательного результата. В -случае необходимости за­ дача может -быть решена с более высокой точностью. Как показы­ вают решения других задач устойчивости [44], выполненные с по­ мощью данного метода, погрешность нулевого приближения, в сравнении с точным, получается небольшая, в пределах 10% . Для нашей задачи погрешность будет зависеть от характера изменения жесткости полоски от участка к участку, а также геометрических размеров оболочки и вырезов.

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

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

выражение для функции р путем подстановки выражений до' и w"y найденных ранее, в отношение (7.26). При этом учитываются за­ висимости (7.21). Кроме того, определяются коэффициенты функ­

ции р по (7.27) с использованием (7.19). Имея выражение для р,

можем записать общее решение уравнения (7.14) в первом при­ ближении. Удовлетворяя граничным условиям, вновь приходим к системе однородных уравнений относительно Аю, А 20, Лзо и Л40.

После приравнивания определителя системы, составленного из ко­ эффициентов при неизвестных А#, нулю, получим алгебраическое уравнение, из решения которого совместно с (7.21) находятся не­ известные величины г: и Р. По окончании расчета первого прибли­ жения вместо (7.38) получим новое выражение для критического напряжения:

ак р

( 7 . 3 9 )

где | — коэффициент ( £ < 1), изменяющийся от приближения к при­

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

Если предположить, что диаметр кругового выреза значительно меньше длины и диаметра оболочки, то справедливо будет счи­ тать действие реактивной силы отброшенной части оболочки оди­ наковым по всей длине выделенной полоски.

Далее рассмотрим некоторые частные случаи. Вначале иссле­ дуем устойчивость нагруженных осевой силой оболочек с опреде­ ленным числом вырезов квадратной формы. Прежде всего, подстав­ ляя в (7.38) зависимости (7.15) и (7.5), получим

( 7 . 4 0 )

где сткр определяется формулой (7.11) или (7.13),

да у

( 7 . 4 1 )

i J и * )

f

О

Судя по зависимости (7.41), коэффициент к зависит от числа, фор­ мы и места расположения вырезов.

Рассмотрим устойчивость цилиндрической оболочки, ослаблен­ ной четырьмя вырезами, расположенными симметрично относитель­ но оси на средней параллели. Длина стороны выреза равна а. Тогда расстояние между центрами вырезов будет равно /i= 2jtR /4. Момент

инерции полоски шириною 1\, мысленно выделенной из оболочки,

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

В пределах каждого участка он постоянен. На участках от нуля до

(I а) 12, а также от

(1 + а )/2 до I жесткость полоски единичной

ширины определяется

по зависимости (7.16). На участке с коорди­

натами ( / — а )/2 . . . (

1 + а ) /2 жесткость полагается

постоянной по

всей ширине полоски,

равной среднему значению,

получающемуся

к

О

0,7 О,U

0,в

0,6 ф

 

Рис.

7.2

Рис. 7.3

делением истинной жесткости на ширину участка 1\. Тогда для по­ лоски единичной ширины она равна

(7.42).

Подставляя (7.16) и (7.42) в (7.41) и произведя интегрирование, получим

На рис. 7.2 приведены кривые, показывающие характер изме­ нения величины параметра k в функции от отношения длины квад­ ратного выреза к диаметру и к длине оболочки. О степени влия­ ния на устойчивость числа вырезов можно судить по графику, при­ веденному на рис. 7.3. Расчет был выполнен для оболочки с пара­ метрами a l(2 R )= 0,2 и а /1= 0,4 и различным числом квадратных вырезов, расположенных в среднем сечении.

Предлагаемый метод может быть применен для исследования устойчивости оболочек с подкрепленными вырезами. Подкрепле­ ние подбирается в соответствии с решением задачи прочности тако­ вым, чтобы оболочка получилась равножесткой во всех направле­ ниях. Если же подкрепление задано, то площадь его поперечного сечения можно равномерно распределить («размазать») до толщи­ ны стенки оболочки, уменьшая действительное отверстие. Форма фиктивного эквивалентного выреза аналогична действительной, а его контур эквидистантен.

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

7.3.ЦИЛИНДРИЧЕСКИЕ ОБОЛОЧКИ ПРИ ВНЕШНЕМ ДАВЛЕНИИ

Рассмотрим устойчивость длинной цилиндрической оболочки, 'ослабленной системой прямоугольных вырезов, расположенных на равном расстоянии друг от друга вдоль образующей. Предполага­ ется, что при потере устойчивости кривизна оболочки изменяется лишь в окружном направлении, а образующие остаются прямоли­ нейными (рис. 7 .4). Это предположение позволяет приближенно судить об устойчивости оболочки по устойчивости кольца, мыслен­ но выделенного двумя перпендикулярными оси сечениями. Анало­ гичные предпосылки были приняты в работе [68]. В отличие от нее здесь избран другой метол решения, позволивший получить окон­ чательные простые зависимости для определения величины крити­ ческого внешнего давления.

Кольцевой элемент, мысленно выделяемый из оболочки, при­ нимался по ширине равным расстоянию между центрами вырезов по образующей. Действующий на него изгибающий момент М, как известно [93], связан с изменением кривизны следующим соотно­ шением:

£ /у.= М ,

(7.44)

тде Е — модуль упругости; I — момент

инерции поперечного сече­

ния кольца, который для исследуемого случая является функцией угла 9.

Решение поставленной задачи ищем с учетом малых дефор­ маций.

Путем простых операций дифференциальное уравнение (7.44) преобразуется к виду

d2w

W = MR2

(7.45)

Пт

El

 

где Я — средний начальный радиус кривизны кольца.

Изгибающий момент появляется вследствие радиального упру­ гого смещения. Он равен M = N w , где N — нормальное усилие, возникающее в поперечном сечении кольца. Величину нормально­

го усилия принимаем

так

же,

как в работе

«а

]68], постоянной и равной 1среднвму значению:

 

N = p R l (*

2^яГ )'

<7 - 46>

 

Здесь п — число прямоугольных вырезов в по­

 

перечном сечении оболочки;

5 — площадь се­

 

чения одного выреза.

 

 

 

Рис. 7.4

Подставляя (7.46) в выражение для момента

M — Nw, а по­

следний в уравнение (7.45), получим

 

 

 

d2w

k w = 0,

(7 .47)

 

 

d02

 

 

 

 

где

k = k (B )

p& L Г i _

(7 .48)

 

 

£ /(0 )l

 

Функция прогиба, найденная из решения уравнения (7 .47), дол­ жна удовлетворять условиям симметрии изогнутой формы, выра­ жающимся в том, что

———= 0 при 0— 0

и 0 = — .

(7.49)

dd

2

 

Линейное дифференциальное уравнение (7.48) с переменным коэффициентом k может быть решено, как и в предыдущем разде­ ле, методом последовательных приближений. Сравнивая уравнение (7.47) с (7.16), найдем, что

п = 2; / 2 (0 )= 1 ; / 2 (0)=О ; / 1(6 )= 0 ;

' • №- ' ^

[ 1- w

]

+

, - ' F ( e ,+ li

(7- 50>

fi( 0 )= p F '( e ) ;

/

(0)— 0.

 

Характеристическое уравнение (7.18)

будет иметь вид

 

r 2-f-a1r - }- a 0= 0 .

 

 

(7 .5 1 )

Решая его относительно г, получим, что

 

 

гх— a J 2 -\-k] г2=

— a J 2 — k,

 

если

k2= a i /4 — % > ();

(7 .52)

Г\= ~aJ2-\-qi\ г 2 =

—aJ2—qi, если

д2—а0—а ? / 4 > 0 .

 

Коэффициенты уравнения (7.51) а0 и а\ пока неопределенные,

т\ФН Оот^О. Уравнение -(7.23) для нашей задачи имеет вид Кроме того, согласно (7.24) имеем

D = 1 1

— г2 f i = co n st ф 0; ах= г 2; а2— гл.

г2

 

Подставим далее (7.50) в зависимости (7.25) и (7 .27):

 

|L=A1-a// -j-X2'ZE//;

 

xl = p m . - i ; x2 = - f i w _ i i

(7 .5 3 )

«о

Оо^(в) Ч

 

Общее решение дифференциального уравнения (7.47) в соответст­ вии с (7.23) имеет вил

w

0 Чъ)

+ Г^ 2°е''’ +

^

J ^ " >м ~

= ~

 

 

e r»8

^ (ie

r*eof0l;

 

 

 

ГЧ.

 

 

J

 

 

 

Оо

 

 

 

 

а I

 

*У2

£ er«e ^ jxe_ r‘flrf0 —

 

■ш' = ЛюГ1ег,в+ ^ 2 о ^ е г»0 -|

 

 

 

 

 

П — гх\

 

 

 

— ег»в|‘ (1.е

 

Га°^/вJ;

(7.54)

.^"=Л 10г1ег,°+ Л20г|ег , в

Гrxe^e ^ [Te_ri8d6—

1

о0

— г2ег»еf {Ге~Г1°^б1.

Оо

J

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

Рассмотрим вначале нулевое приближение. В этом случае jx=Of тогда решение преобразуется таким образом:

® =

j (г1-4 1оег‘*+''2Ла)е'>1);

 

-а!>'=Л10г1е'««+Л 20г2е'*«;

(7.55)

=Л10г?ег'° + Л2ог1ег»0.

Вполученное решение входят неопределенные параметры Сь До

инеопределенные корни Г\ и г2. Для нахождения а<> и а\ используем метол, применяемый при решении интегральных уравнений, суть которого состоит в осреднении соответствующих функций. На ос­ новании этого неопределенные параметры могут быть вычислены по следующей формуле:

/,(« )< » . 7 = 0 > ] > 2 . - м П.

(7 .56)

Осредняя функции /0(6) и f i(0) на рассматриваемом интервале по зависимостям (7.55), найдем

a' = 1

7

=

^ /,(e )d e ;

 

О»

 

 

Оо

0,

 

 

 

 

 

а 0= I— \ f ( Q ) d B = — £— ^ ( 6 ) ^ 0 .

(7 .5 7 )

aiао J

 

 

®i— во

«)

 

во

 

 

 

О0

 

Затем по зависимостям (7.57)

и

(7.52)

определим

неизвестны^

корни

 

 

 

 

 

r x=qi\ r2 = — qi\ Ч=

л

/

^ ( е)^ е +

(7 .5 8 )

 

 

Y

 

Оо

 

Решение (7.55) содержит

произвольные постоянные А \0 и Л20.

Их можно найти, удовлетворяя дополнительным условиям, которые следуют из условия симметрии выпученного кольца. Эти условия достаточно выписать для двух .поперечных сечений А к Б (см. рис.

7.4) с координатами 0 i= O и 02=

jt/2:

 

700,-0= 0;

7О0,-«/2=О.

(7 .5 9 )

Объединяя (7.5Э) и (7.55), получим систему алгебраических уравнений, с помощью которой можно определить постоянные Лю и Л 20*'

r 1 еГ1®' Л 10+ г 2ег»°» Л20= 0 ; гхe rifl* Л 10-f г2 ег*в» Л20= 0 . (7 .6 0 )

Эта система будет иметь решение, отличное от нудя, в том слу­ чае, если определитель, составленный из коэффициентов при неиз­

вестных, равен нулю:

 

 

гхег,0>

г2 ег*в‘ = 0.

(7 .6 1 )

гхег,в>

г2ег*в»

 

Учитывая соотношения (7.58) и раскрывая определитель, най­

дем, что sm«/i(02 — 0 i ) = 0 . Отсюда легко определим

наименьший

положительный корень, отличный от нуля: <7 (02 — 0 i)= jt или, под­

ставляя вместо 02 и 0i их значения, найдем q = 2 . Найденное зна­ чение q подставим в (7.58), в результате получим формулу для вы­ числения критического давления в нулевом приближении:

3 (Bi — 0О)

(7 .6 2 )

0>

|>(б)<го

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

P i с р — 5

EI

(7 .6 3 )