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

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

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

стоянии b от края оболочки). Выполняются краевые условия свободного опирания.

Экспериментальной обработке подвергается формула (3,19,18), т. е.

(3,20,1)

где Ма и То — критическая комбинация действующих нагрузок; М*кр — критическая сумма четырех локальных моментов, действующих на обо­ лочку с подкрепляющим поясом, когда осевое сжатие отсутствует;

Т кр — критическая сумма четырех локальных осевых усилий, когда мо-

ментная нагрузка отсутсувует. По результатам первой части этого па­ раграфа для четырех осевых сил имеем Т кр =0,4 Гкр. где Тщ> — крити­ ческое осевое усилие, когда на оболочку действует равномерно распре­ деленная нагрузка. Естественно предположить, что жесткость подкреп­ ляющего пояса оказывает существенное влияние на величину критиче­ ского изгибающего момента. Наличие подкрепляющего пояса повлияет и на величину критического осевого усилия, но, как показывают экспе­ рименты,, когда жесткости оболочки и пояса сравнимы друг с другом, это влияние несущественно и оно идет в запас устойчивости.

Представим (3,20,1) в виде

/

Мр у

/_л_у

+

То

 

(3,20,2)

\ м

кр;

\ н )

0,4Т, =

1 ,

 

где Л1Кр — критическая сумма четырех локальных моментов, действую­

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

m кр —

cEh3

(3,20,3)

(1 — V2) 1?"

 

 

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

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

^кр

 

c-jEhi3

(3,20,4)

(1

- v 3)R'

 

 

где коэффициент С\ подлежит определению из эксперимента. Уравнение (3,20,2) теперь запишется так:

Мо_ ( l - v a)Ky

/_Л_у ,

Г»

(3,20,5)

С!

Eh3 ]

\ Н /

0,4Гкр

 

После проведения серии экспериментов, в которых изменялись гео­ метрические параметры подкрепляющего пояса 26 и Я, а также крити­ ческие комбинации нагрузок М0 и Т0 при фиксированных размерах поя­ са, были определены неизвестные коэффициенты С\, р и ц. Окончатель­ ный вид формулы для определения критической комбинации М0 и Т0:

:М0

( 1 - у 3) R 1P

/ h

у

Гр

(3,20,6)

0,267 '

Eh3 J

\ н

)

0,4Гкр

 

где р находится по графику рис. 3.60, показывающему зависимость р от относительной ширины подкрепляющего пояса 2 b/R.

f

3.Устойчивость оболочки, ослабленной круговым отверстием

В§ 19 гл. 3 получена формула (3,19,20) для критического усилйя сжатия, содеожащая неизвестную функцию /, зависящую от параметра р2IRh.

Экспериментальное исследование, в ходе которого испытывались

оболочки с различными отверстиями, выявило зависимость f от р2/Rh. Всего было испытано 30 оболочек. Зависимость / от р2/Rh показана на графике рис. 3.61, где кружками и треугольниками отмечены экспери­ ментальные точки *. Экспериментальная кривая проведена так, что все точки лежат на кривой или ниже. Это сделано для того, чтобы при на­ хождении Т0щ) по формуле (3,19,20) и графику на рис. 3.61 результат получался наименьшим, т. е. чтобы запас был в сторону устойчивости1

1 Экспериментальная кривая на рис. 3.61 соответствует локальной потере устой­ чивости оболочки с образованием первых волн в зоне концентрации напряжений около’ отверстия.

оболочки. Следует отметить, что график построен для оболочек из ма­ териала, имеющего коэффициент Пуассона v = 0,3. Для материалов, имеющих другой коэффициент Пуассона, результаты несколько отлича­ ются от полученных.

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

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

 

(3,19,20), в

которой значение функции f(p2/Rh) находится из графика

на рис. 3.61,

а Г0„р вычисляется для оболочки с такими же параметрами,

но б.ез от­

верстия, с учетом начальных несовершенств.

в

[148—150].

Изложенные в этом параграфе результаты получены

§ 21. ТЕОРЕТИКО-ЭКСПЕРИМЕНТАЛЬНЫЙ МЕТОД ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ КОНИЧЕСКИХ ОБОЛОЧЕК

1. Линейная задача

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

Обратимся к линейным уравнениям нейтрального равновесия ко­ нической оболочки, учитывающим безмоментность ее первоначального состояния. Последние, как и уравнения (3,19,6), можно получить из общих нелинейных уравнений (3,19,2) путем варьирования. Соответ­ ствующие линеаризованные уравнения будут

(3,21,1)

где Гю, Т2о, S\o — мембранные усилия до потери устойчивости.

В безмоментом состоянии в случае равномерного внешнего давле­ ния эти усилия определяются по формулам

Вводя их в уравнения (3,21,1) и сохраняя члены, соответствующие полубезмоментной теории, получим уравнения

\

(3,21,2)

1 d*6w г]2

г дг2 ‘ гз

где

/ = 5ф/£ ctg Y;

Ehr\ c t g 2 у

vl = qr\\gyiD.

т)2 =

;

 

 

D

 

Задаваясь решением уравнений (3,21,2) в форме

6ф = Д (г) cos пхcpi*

6w =

(г) cos

/гх = м/sin у

н вводя переменную х = \п

— , эту систему уравнений приводим к одно-

му уравнению вида

 

ri

 

 

dwi

d2^

+ ту-2 (п\е-2х o*0n6{ex) w1 = О,

 

dx4

dx2

(3,21,3)

 

— / < х < 0, t = \п ——.

 

 

Го

Отметим весьма

интересный факт. Несмотря на различие общих

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

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

жет быть получено решение и для цилиндра.

Различие двух

систем

уравнений (переменные коэффициенты и т. д.)

обусловливает

лишь

степень трудности математического решения соответствующих

задач

для конуса и цилиндра. Однако это различие не существенно для экс­ периментального решения.

В самом деле, уравнения для конуса и цилиндра в смысле наличия -производных по дуговой координате не отличаются друг от друга и имеют одинаковые решения по этой переменной. В силу этого уравнение (3,19,8) при а*=т* = 0 совпадает по своему типу с уравнением (3,21,3).

Допустим, что уравнение (3,21,3) решается любым* из приближен­ ных методов для какого-то вида граничных условий (например, свобод­ ное опирание).

После минимизации критической нагрузки и исключения волново­

го числа П\ получим формулу

 

 

 

 

 

 

 

Як = / ( 0 ,1

%

l t /

— ^

У

J

/2;

* =

(3,21,4)

 

(1

V»)

‘ ‘fit V

rx t g

у

 

Л0

 

где f ( i ) — функция,

подлежащая

экспериментальному

определению.

Для замкнутой в вершине конической

оболочки In —— »оо и f ( t ) / t =

= С = const. В этом

случае формула (3,21,4)

го

 

примет вид

 

 

Як =

с

Eh

 

 

 

 

 

(3,21,5)

 

( l _ v 2 )V V i

V

 

t g

Y

 

 

 

 

 

\

T i t g

у /

 

По результатам

экспериментальных

 

исследований

защемленных

конических оболочек из алюминия [135] константа С в формуле (3,21,5) изменяется от 2,18 до 2,74.

Однако в идеальном случае для определения константы С в этой формуле требуется эксперимент с одним образцом.

Для цилиндрической оболочки длины I в формуле (3,21,4) следует выполнить операцию предельного перехода, а именно

/(О =

/

[ ln

+ у - ] ]

=

/ ( f )

=

/(0 ) =

С,

так как f

-* 0,

 

 

 

r xt =

r x In

=

rx In (

1 +

— ^ =

г il/г о =&/,

 

 

 

 

 

/■о

 

V

 

Г0 /

 

 

 

 

 

 

 

 

 

 

M g y =

tf

 

 

 

 

и формула

(3,21,4)

при этом

будет

совпадать

с

первой

из формул

( 3 , 1

9 , 8

"

) .

формулы

 

(3,21,4)

следует,

что

для эксперименталь­

Из

структуры

 

ного определения функции f ( t ) в опыте необходимо изменить лишь от­ ношение ri/ro. При этом вид этой функции не зависит от угла конусно­

сти. Значения функции

f ( t )

для усе­

ченной защемленной конической

обо­

лочки из алюминия представлены на

рис. 3.62. Кривая 1 построена по ре­

зультатам

испытаний

 

защемленных

усеченных конических оболочек с уг­

лом конусности 2у = 45° при изменении

отношения г\!г0 от 1,143 до 4,0. Кри­

вая 2 представляет теоретическое ре­

шение для функции f{ t)

[135].

 

Выполненное здесь

 

решение част­

ной задачи

позволяет

сделать следу­

ющий вывод. Если деформация кони­

ческой оболочки в дуговом направле­

нии превалирующая,

как

это

имеет

место при внешнем давлении, то

фор­

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

цилиндрической оболочки путем введения в них вместо длины I и ра­ диуса R соответственные величины для конуса r xt и ritgy. При этом в формулы для конической оболочки следует также ввести в качестве множителя функцию вида f { i ) , подлежащую определению из экспери­ мента. Легко видеть, что для определения вида этой функции требуется минимальное количество экспериментов.

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

2. Математическая аналогия в экспериментальных исследованиях

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

;

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

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

h = h0r\ K = Eh0r/(l — ч*) = Ко'П D = Е/г^г3/ 12 (1 — v2) = D0r.

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

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

Учитывая это, уравнения равновесия для оболочки переменной тол­ щины можно получить из уравнения (2,15,23), вводя вместо h значение Л = Л0г. Преобразуем уравнения с помощью подстановок:

 

т

tp =

 

*

w = rx ctg yezwx,

wx =

Ь2

w0;

 

г - In— ;

e2zf;

 

ro

 

 

 

 

 

 

 

 

 

 

 

 

r?

 

 

/о = f/EhjPb2 ctg2 Y;

2 = — ta;

фх = — 0;

/ = In

r0

 

 

 

 

 

 

 

 

 

 

r\

 

 

 

0 <

6 < 1;

1 > a > 0;

hx = h0rx,

 

 

получим уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2w0

 

d2w0

 

d*w0 __ /

d2w0 \ 2 _

Q

 

 

 

 

 

da2

 

da2

 

 

d№

’ V

dadO

)

~~

 

 

D*ViV^o

0Vo

Wo.

 

 

_

d2f0

 

d*w0

,

 

(3,21,6)

 

da2

 

da2

d№

 

 

d№

 

da2

 

 

 

 

 

 

 

 

 

 

 

 

+ 2 - ^ - J ^ + q* = 0,

 

 

 

 

 

где

 

 

dadd

 

dadO

 

 

 

 

 

 

 

 

 

 

 

 

* (

)

frjf_ у

d*(

) .

 

 

 

 

 

V? =

 

 

 

 

 

V ь

)

дв2

 

 

 

 

 

 

 

да2

 

 

 

 

D

_____1

 

( - L ) 4.

<T =

_ j _ ( ± Y-

 

 

62

 

12(1—v2)ft*3

Vrxt

)

 

Ek*3 \ It

)

 

 

V i tgy ’

Ь— дуговая ширина панели при r= ru h\=hori— толщина оболочки

ПрИ Г = Г\.

Из сопоставления уравнений (.3,21,6) с аналогичным образом пре­ образованными уравнениями для цилиндрической панели устанавли­ ваем следующие соответствия:

цилиндрическая

панель

и* _

Ь2

 

Rh

1

а

К = —

ь

q'

 

коническая

 

коническая

 

панель (усеченная)

 

панель

(замкнутая)

 

k' -

62

 

/г* -

62

 

 

rji tg у

 

 

rih tg Y

«-*

 

b

 

Ji — ri

 

 

 

%1~ T -

~

q = q’/%(К, i)

~

? =

<77X M-

Заметим, что из сопоставления уравнений типа (3,21,6) для цилин­ дра и конуса вытекает лишь соответствие (цилиндр—усеченный конус). Соответствие (цилиндр—замкнутый конус) установлено путем дополни­ тельных исследований.

Приближенная аналогия позволяет трансформировать решение для цилиндрической оболочки в решение для конической. Поскольку при выводе уравнений (3,21,6) мы пренебрегли всеми низшими производны­ ми по координате z по сравнению со старшими по той же координате, то возникающие при этом погрешности в решении для конической обо­ лочки, полученном из решения для цилиндра, следует компенсировать приведенной поверхностной нагрузкой q = q * / % ( K t ) , где % ( k , t ) или %(Х)— функции, подлежащие определению из эксперимента. Вид этих функций не зависит от угла конусности у. Другими словами, как и выше при экспериментальном решении задачи, выполнение условий геомет­ рического подобия не обязательно. Установленная приближенная анало­ гия полезна при решении задач с помощью экспериментальных мето­ дов.

В рассматриваемой здесь задаче значения функций х(^, i) опреде­ лялись по результатам испытания на устойчивость жестко защем­ ленных по всем кромкам конических панелей под действием равномер­ ного внешнего давления [135]. Конические панели изготавливались из алюминия и константана. Были испытаны панели с углом конусности 2у = 45°, 90, 130°, отношением ri/r0 =1,2; 1,5; 3,0 и углом раствора а в плоскости нижнего основания, изменяющимся в пределах от 10 до 60° Для каждой испытанной панели замерялось критическое давление хлопка q x л, соответствующее появлению первой волны, и критическое давление выхлопа qBых (нижнее критическое давление, соответствую­ щее выпученной форме равновесия панели).

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

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

ния. Так как функция %(Х, t) в формуле (3,21,8) не зависит от угла ко­ нусности, то значения ее определялись по результатам замера крити­ ческого давления-^хл для панелей 2у = 90° при переменных X и t.

На рис. 3.63 и 3.64 представлены значения функции %(Х, t) ДЛЯ алюминия с коэффициентом Пуассона v = 0,26.

На рис. 3.65 и 3.66 представлены значения функции х(Х, t) для константана с коэффициентом Пуассона v=0,31.

Zttj

Значение функции x(X, t) зависит от материала, что обусловленовлиянием коэффициента Пуассона. Функцию х(^> О для панелей с п/г0г^;1,2 можно считать постоянной для всех значений X.

Для панелей

с — = 1,2; х(^>

^) = 1,61 (алюминий) и x(?t,

t) = 1,35

(константан).

r0

 

 

давление qh для

конической панели любого

угла ко­

Критическое

нусности находим по формуле (3,21,8), определяя для данной панели с параметрами X и. t значение х(^> 0 по приводимым графикам. Крити­ ческое давление qh по формуле (3,21,8) получается завышенным для большинства панелей на 10—15% по сравнению с экспериментальным значением qxЛ. Можно получить нижнюю границу для критических дав­ лений конической панели, если при определении функции %{Х, t) в фор­

муле (3,21,8) исходить

из экспериментальных значений критического

давления выхлопа qBых.

'

3.Локальная устойчивость конических оболочек

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

очень короткие волны, длина которых порядка V Rh, для исследования устойчивости можно применить уравнения локальной устойчивости. Эти

уравнения при допущении в них внутренней погрешности порядка / I

по сравнению с единицей являются уравнениями с постоянными коэф­ фициентами и содержат в себе лишь высшие производные. Эта идея упрощения уравнений, первоначально использованная в линейной тео­ рии краевого эффекта, в теории устойчивости впервые, по-видимому,, была применена И. Я. Штаерманом [138]. Позднее Ю. Н. Работнов раз­ вил эту теорию, назвав ее теорией локальной устойчивости оболочек [140]. В монографии [139] В. 3. Власов обобщил эту теорию на случай нелинейных задач. В статье [137] А. В. Саченков уточнил процедуру применения уравнений локальной устойчивости, показав, что последниеможно использовать не только в тех случаях, когда область выпучива­ ния незначительна по сравнению с характерными размерами оболочки, но и когда наблюдается общая потеря устойчивости с образованием одной полуволны по длине оболочки, например, в случае конической оболочки при внешнем давлении. Применительно к конической оболоч­ ке ниже излагаются результаты, принадлежащие этому автору.

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

Преобразуем уравнения равновесия (3,19,2'), вводя следующие подстановки:

z = In— ;

ср = ezf; w = rx ctg yezwx.

(3,21,9)

При этом получим систему

двух нелинейных уравнений

относительно

f и wx. Эта преобразованная система уравнений в свою очередь может быть значительно упрощена.

Уравнения равновесия, которым функция усилий ф удовлетворяет тождественно, приближенные. Они не учитывают влияния перерезыва­ ющих сил. Неучет этих сил эквивалентен пренебрежению единицей по сравнению с квадратом числа окружных волн п2. Если оболочка не является длинной, то п2^ 1, и пренебрежение единицей по сравнению с п2 не ведет к значительным погрешностям. Математически это озна­ чает пренебрежение функцией по сравнению с ее второй производной' по фЬ т. е. функции при дифференцировании считаются существенно

возрастающими.

'cPwi

d2w1

Допуская, что

и пренебрегая в преобразованных с по-

d z a

дер2

мощью (3,19,9) уравнениях функцией по сравнению с ее второй произ­ водной, приходим к следующей системе:

v V

/ =

&

d2Wx

d2Wx

dwx

d2Wi

? c t g * v « *

( ( ^ ) *

 

 

 

 

 

~д&

дер2

dz

дер2

 

 

 

dw±

d2W\

 

d2Wi

 

dwx

 

 

 

 

 

 

 

dz

 

dz2

 

 

dz2

 

dz

J ’

 

 

 

D y 2y

2W j_ -

 

 

- +

-

d f

I

/

d2wx

1

 

\

a»/

+

l

dz2

dz

4

 

dz2

dz

 

 

 

 

 

 

П

)

a<p2

 

 

 

d2wx x J

 

 

df_

dWl \

2

d2f

 

d2wt

Л з г

d f

\

-

f .

 

1 mUл 5=- 0V7.

dz

J

'

dz

V,

dz

 

 

dz

)

 

dzd(pi

dzd(pi

c tg у

(3 ,2 1 ,1 0 )

Эти уравнения будем считать исходными для несимметричной дефор­ мации конической оболочки при /г2>1. Интегрирование их представляет большие трудности. Однако они могут быть сильно упрощены в случае весьма распространенных видов нагружения, сопровождающихся появ­ лением мелких волн во взаимно перпендикулярных направлениях. При этом функция усилий и функция прогибов изменяются как быстро воз­ растающие функции и отношение двух последовательных производ­

ных величин порядка л /

JL. Поэтому в случае тонких оболочек, когда

величина

мала

по

сравнению

с единицей,

будем

пренебрегать

функциями f и Ш] по сравнению с их первыми производными.

Таким

образом,

уравнения

(3,19,10)

заменяем

приближенными

уравнениями

 

 

 

 

 

 

 

 

 

 

 

у2у2/ = Er\ ctg2 yez

 

 

d2w1

d2Wi

 

 

 

 

 

~~d#

~dfi

 

 

 

 

 

 

 

 

 

 

 

 

D y 2v 2ffi'j — he2°

hez° j

5%!

 

 

 

(3,21,11)

 

0 *

'

а Ф2

+

 

3<p2

&f_ _

2

 

d2f )

r]ezc tg yq = 0.

 

L&2

 

dzdipi

dzdrfi J

 

 

 

 

 

Так как в этих уравнениях сохраняются лишь старшие производные, а переменные коэффициенты при них, являющиеся геометрическими ха­ рактеристиками оболочки, изменяются плавно, то последние можно рассматривать как постоянные параметры, подлежащие определению из условия минимума критической нагрузки. Это соответствует рассмот­ рению устойчивости малой части оболочки вблизи z= zc, метрику кото­ рой можно считать евклидовой. Поэтому в уравнениях (3,21,11) мы по­ ложили

ez — eZf = const.

Выполняя в (3,21,11) дополнительные преобразования

/ = F/ezc\ wx = wlrx ctg уег< ; z = xfrxe \ cpx = у/гхег',