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

книги из ГПНТБ / Рикардс Р.Б. Устойчивость оболочек из композитных материалов

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

4.2.

Устойчивость под действием

внешнего давления

123;

оболочки

h/R 1/10, 1/25

 

 

и 1/100;

график построен

 

 

для

 

величины

RjL= 1.

 

 

Параметр h-*~1, когда мо­

 

 

дуль межслойного сдвига

 

 

стремится к бесконечнос­

 

 

ти. Для сравнительно тон­

 

 

ких оболочек (кривая 3)

 

 

влияние поперечных сдви­

 

 

гов при реальных величи­

 

 

нах

модуля

межслойного

 

 

сдвига

мало

ощутимо. В

 

 

то же время, у сравни­

 

 

тельно

толстых

оболочек

 

 

(кривые 1 и 2 )

это влия­

Рис. 4.2.2. Кривые изменения

параметра k~

ние

значительно. В слу­

в зависимости от величины модуля меж­

чае

внешнего

давления

слойного сдвига.

 

на

критическое состояние

 

 

конструкции и на форму потери устойчивости влияет относи-- тельная длина оболочки. На рис. 4.2.3 представлен график из­ менения параметра k в зависимости от длины оболочки. Кривые- 1—4 соответствуют параметрам | 4 = ^5 = 0,1; 0,05; 0,02 и 0,01; для всех кривых принято h/R= 1/10. С уменьшением длины увеличи­ вается разница между результатами, полученными на основе классической теории и найденными с помощью расчетов по фор--

муле (2. 15).

Из рис. 4.2.3 видно, что у коротких оболочек, например,, когда L/R = 0,5, погрешность, вносимая при расчете с использо­ ванием гипотезы Кирхгофа—Лява, даже в случае большой;

Рис. 4.2.3. Зависимость критического внешнего давления от длины оболочки. Обозначения те же, что на.

рис. 4.2.1.

Глава IV . Устойчивость оболочэк при статических нагрузках

124

■сдвиговой жесткости (кривая ]), значительна. В случае длинных оболочек расчет критических нагрузок по классической теории и теории типа Тимошенко дает близкие результаты. У коротких оболочек большее влияние поперечных сдвигов на критическую нагрузку можно объяснить тем, что с уменьшением длины увели­ чивается число волн, по которым' выпучивается оболочка, т. е. локальные кривизны больше как в поперечном, так и в продоль­ ном направлениях. Следовательно, при выпучивании коротких оболочек деформации поперечных сдвигов больше, чем длин­ ных. Эксперимент, проведенный на оболочках из стеклопластика, показал, что под действием внешнего давления критическая на­ грузка падает с уменьшением длины по сравнению с теорети­ ческим расчетом по классической теории. Отметим также, что с уменьшением модуля межслойного сдвига изменяется форма потери устойчивости. Так, например, для оболочки с парамет­ рами h/R = 1/25; R/L = 1; £4 = £5 = оо; 0,05; 0,01 число волн по ок­ ружности, при котором значение критического давления (2.15) имеет минимум, /г*= 5,7 и 8 соответственно. С уменьшением мо­ дуля межслойного сдвига увеличивается число волн, по которым

стремится

выпучиваться оболочка. Для иллюстрации на

рис. 4.2.4

приведены кривые, которые характеризуют зависи­

мость внешнего давления от параметра волнообразования. Гра­ фик построен для h/R= 1/25 и R/L= 1; кривые 1—5 соответст­ вуют параметрам £4 = | 5 = оо, 0,3; 0,1; 0,05 и 0,01. Параметр волнообразования может принимать только дискретные целочис­ ленные значения вблизи точки минимума (на рисунке обозна­

чены кружками).

устремив £4 = ^5->-0, можно получить

Отметим, что из (2. 15),

формулу

для

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

усилия, когда модуль

 

 

 

межслойного сдвига стремится

 

 

 

к нулю. Интересно отметить,

 

 

 

что

критическое

давление

в

 

 

 

этом

случае

не равно нулю, а

 

 

 

составляет

 

примерно

5—10%

 

 

 

от

величины

классического

 

 

 

давления. При этом оболочка

 

 

 

стремится

выпучиваться с об­

 

 

 

разованием

большого

количе­

 

 

 

ства волн по окружности.

 

 

 

 

Под нагрузкой оболочка из

 

 

 

стеклопластика до потери ус­

 

 

 

тойчивости

 

может

расслаи­

 

 

 

ваться

[66],

следовательно,

в

 

 

 

некоторых зонах модуль меж­

разования.

1—5

соответствуют

слойного

сдвига оказывается

54=15=°°;

0,3; 0,1; 0,05; 0,01.

равным нулю и нагрузка, кото­

4.3. Геометрически нелинейная задача

12S

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

Из приведенных расчетов можно сделать вывод о пределах применимости гипотезы Кирхгофа—Лява при расчете на мгно­ венную критическую нагрузку цилиндрической оболочки под

внешним давлением. При относительной

толщине

h / R ^ 0,04,

длине

L/R^ . 2 и

при величине параметра

модуля межслойного

сдвига

^4 = ^5 —0,1

расчет следует вести

на базе

уточненных

теорий.

 

 

 

 

4.3. ГЕОМЕТРИЧЕСКИ НЕЛИНЕЙНАЯ ЗАДАЧА

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

Известно [44], что более полную картину потери устойчивости и послекритического поведения дает решение в геометрически нелинейной постановке. В экспериментальной работе [66], посвя­ щенной исследованию стеклопластиковых оболочек, установ­ лено, что после потери устойчивости у оболочек с параметром R/h = 20—150 наблюдается сдвиг слоев. Это свидетельствует о превышении предела прочности на срез. Классическая теория оболочек (гипотеза Кирхгофа—Лява) и уточненная геометри­ чески линейная теория не дают возможности достаточно убеди­ тельно объяснить этот факт.

Рассмотрим задачу об устойчивости в «большом» цилин­ дрической ортотропной стеклопластиковой оболочки под дейст­ вием внешнего давления*.

Оси х, у, z направим вдоль образующей, по окружности и нормали. Примем, что оси ортотропии материала совпадают с главными направлениями оболочки. Считаем, что оболочка имеет начальные неправильности в виде геометрических откло­ нений срединной поверхности от идеальной формы. Используем систему уравнений устойчивости (3.5.25) для оболочек боль­ шого прогиба с начальными несовершенствами, заменив упруго­ вязкие операторы упругими:

(Vy—Yv(0). Y*~Y*(0)}

+

 

+ НАыЦ\2{уу — yv{0),

®—“>«»} =0;

(4-3. 1)

* Подобная задача с использованием классической теории решена в

{31, 225 и др.].

Глава IV . Устойчивость оболочек при статических нагрузках

126

У= 4 ;1I{’YX—Yx(0), Yy —Yy(0)} +

(4.3.2)

+ ^55'Пп {Yx —y.t<0). o' —o>(o)} =0;

z ^ X (Y x - Yx(0), YyYy(0)} -

 

 

 

(4.3.3)

1

d2 { w - w (0)) _

 

V=V*F+ —

dx2

 

HR

(4.3.4)

- k {w , w} +Х{Ш(0)>аУ(о)} =0.

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

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

Принимаем, что форма начального волнообразования нахо­ дится «в резонансе» с волнообразованием оболочки в процессе деформации. Такое допущение несколько увеличивает влияние начальных неправильностей в сторону запаса устойчивости [44], но значительно облегчает решение. Тогда для начальных проги­ бов и углов поворота нормального волокна

-

w0

 

^

w0= - = *i0sinoasin By+ х2Эsin2 ах;

.

 

h

...

•у**0*= *3ocos ax sin $y + xw sin 2oa;

 

 

Yv(0) = x50 sin <xxcos $y.

 

(4.3.6)

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

.4-3. Геометрически нелинейная задача

127

Подставим (3.5) и (3.6) в уравнение (3.4). После интегри­ рования этого уравнения получим выражение для функции на­ пряжений, в граничных условиях которой учтено, что на обо­ лочку действует внешнее давление (напряжение 022° = —qR/h):

F=-

=[

х102) ~

h2s

--^гЬрл-Ярг^г' 1 2 — х20) ]cos2ax + -r^-npiX

Хп -^ ц ^ !2 —Л'ю2) cos 2 ^y+ — (xlx2 - x l0x20) х

У/2j

Xsin3axsin -{—^-[«“2р2- 1(*1 —х10) —

 

, k 2

 

dR

(4.3.7)

- (x1x2- x Io^2o)]sin ах sin fiy— ^ j - x 2.

Здесь безразмерные параметры материала и геометрии оболочки определены по (2. 12), а коэффициенты ki и k2 следующим спо­ собом:

^1 = 9ц % 2 {n~2K2pi2 --2 U l{-%2- x\xrx+ U -{+ ^ n 2n~2p r 2l r lP~l :

k 2 = ц~1\2~1п~2я2рI2—2 |3£1~112-1 + |б~1+ п2л~2р\-21 г 1ц_!•

Далее решаем систему (3.1) —(3.3) методом Бубнова—Га- леркина, подставив в правые части этих уравнений разложения (3.5) и (3.6), а также выражение (3.7) для функции напряже­ ний. Согласно этому методу, интегрированием по поверхности оболочки необходимо найти следующие уравнения:

2ЯЯ L

I I Z sin ах sin fiydxdy = 0;

2яН Ь

J J Z sin2 axdxdy = 0 \

оо

2ЯЯ L

II Y cos ах sin fiydxdy = 0;

оо

2ЯЯ L

2яЯ

L

 

J

J X sin ах cos $ydxdy = 0.

I

S

У sin 2аxdxdy = 0;

(4. 3. 8)

О

0

0

0

 

 

Глава IV . Устойчивость оболочек при статических нагрузках

12а

После интегрирования получаем систему, состоящую из двух, нелинейных и трех линейных алгебраических уравнений относи­ тельно пяти неизвестных х4 (i 1, 2,... , 5):

У 1 — a i-'6 3 + О 2 Х 1 Х 2 2 — 03X 1X 2 + Q4X1 + 05X2 + 03X3 + О7Х5 +

Og — 0;

У 2 — Ь \Х \ 2Х 2 — Ь 2Х 12 - \ - Ь з Х { + 6 4 X 2 + ^5X4 + Ь 6 = 0 ]

 

Уз— С\Х\ — с2х3- с3х5 + сА= 0;

 

У*= d[X2 — d2Xj + d3 = 0;

 

г/s = / 1-^1 — / 2-^3 — / 3-^5 + / 4 = 0 .

( 4 . 3 . 9 )

Здесь коэффициенты a, b, с, d, f выражаются через геометричес­ кие и физические параметры материала, величину нагрузки, ам­ плитуд начальных прогибов и углы поворота нормального во­ локна. Представим коэффициенты а* для первого уравнения системы:

Oi=—-l (Zin2pl2n - 2 +%2n - 2p r 2n2);

“2“

х

ж

+ ~ й г ) ;

 

 

«з=4- Ь ^ _2р г 2Р2_! + 4 - «~2р2“';

 

 

О

 

 

/22

 

 

« 4 =

0 ,

• 0

- 4 р 2 - 2 ------— q n ~ 2p

r 2 - a i X w 2 -\— — ^ 2 l t _ 2 p + 2p 2 _ I ^ 2 0 ;

 

 

Z&2

 

 

£

О

 

0,5 =

1

 

-2

-1

 

 

" о и ~ П

2Р 2

Хх\0 ~ а 2х \Ъх 2Ъ,

 

 

1 к<1

 

 

 

 

« 6 = -2 ^ " [liJtp i« _2P2_ 1+ (|з +

2 5 б ) я -,р г 1р2- 1 ];

 

 

 

[Е2П- 2р Г 2Р2- 1 « + ( | з + 2|б) « “ ’рг- 1 ];

 

« 8=

-« б Х з 0 -О 7Х50 + - 4 — (M- 2p2_ 1X l0 X 20 -n _4P2_2JClo).

(4. 3. 10а)

 

 

 

 

2к2

 

 

Таким же образом получены коэффициенты для второго уравне­ ния системы:

4.3. Геометрически нелинейная задача

129

Ьз= - ~ п - 2р2 - 1х 1й\ J.K2

Ь^=~%2П-2р Г 2П~2^2~2\

Ьь= - ^ 1 \Щ\Р2~1П-2\

6б= -у^-|2Л'-2р1~2р2-1-Х:К>2 — ^5 X4 3 — Ь^Х2 0 -

(4.3. 106)

Коэффициенты di и fi содержат параметры межслойной сдви­ говой жесткости:

С1=|5Яр1Р2,

С2 ——^(^,\Л2р12р22 + ^6р22П2) +£5;

 

С з ~ ~ ^ 2 (^3+

1б) Л Р 1Р 2 2ГГ, с 4 = С2Х 30 +

С3Х50 — С\ХГю;

 

^1 = Е 5Я р1р2,

^ 2 = _T 'E l I1t2p i 2p22 +

i5 ,

d 3 = d 2X^Q

й [Х 2о1

 

о

 

 

 

/l=^4«p2,

/2=-у2'(|з + 5б)яр1Цр22;

 

 

[з= _^2'(^бЯ2р12р22+ ^2р22я2) + |4;

fi = f2^30“Ь/3^50

fiMo* (4. 3. 10в)

В выражениях (3. 10а) — (3. 10в) безразмерные параметры попрежнему определяются формулами (2. 12).

Решение системы (3. 9) можно найти итерационным методом Ньютона—Рафсона на ЭЦВМ, для чего необходимо задавать начальные приближения к решению и использовать способ «по­ следовательного нагружения» [145], т. е. создать картину, подоб­ ную нагружению в реальной конструкции.

На рис. 4.3.1 сплошной линией показано нагружение в мо­ дели, прерывистой — нагружение в реальной оболочке. Сначала находим верхнюю критическую нагрузку qB по формуле (2. 15), затем шагом, например 0,01 qB, изменяем нагрузку, причем за первое начальное приближение неизвестных Xi принимаем их значения при q = 0, т. е. величины начальных амплитуд прогибов и углов поворота нормального волокна. В итоге мы получаем

9 - 1744

Глава IV . Устойчивость оболочек при статических нагрузках

130

Рис. 4.3.1. К определению верхней предель­ ной и нижней критических нагрузок.

кривую нагрузка—амплитуда прогиба. Сходимость при решении системы (3.9) была хорошей только в докритической стадии. Для закритической стадии (после «хлопка» в модели) по­ лучить кривую нагрузка—амплитуда прогиба не удалось. Для улучшения сходимости система (3.9), состоящая из пяти урав­ нений, была преобразована в систему двух уравнений путем ис­ ключения из первого и второго уравнений системы (3.9) неиз­ вестных х3, х4 и х5. Из четвертого уравнения системы получаем выражения для неизвестного

х4=

d\X2 + d3

(4.3.11)

 

d2

Решая совместно третье и пятое уравнения системы (3.9), полу­ чаем выражения для неизвестных х3 и х5:

/4С3 —сДз_______Cjf3 C3fi ^

f2c3 — c2f3 / 2 ^ 0 £2 / з

^2/4—£4/2_____Cif2 — fjC2 ^

(4.3.

12)

Czfz — C'Ji

c2f3 — c3f2

 

 

Подставляя (3.11) и (3.12) в первое и второе уравнения сис­ темы (3.9), получаем систему двух нелинейных алгебраических уравнений относительно двух неизвестных Х\ и х2, сумма кото­ рых f= x i+ x 2 представляет амплитуду прогибов:

4.3. Геометрически нелинейная задача

131

у 1—а\Х1ъ + а2Хух22 — агх хХ2+ у а 4—

 

а б(с \ ! з ~ сзЫ -\-a 7 {c \ h ~ flC2)

Ху+

 

 

/гсз — сз}з

 

 

 

 

 

a7{c2fi~cif2) - =

 

+

+ G8 +' аб((4с3 — с4з)

0;

 

 

сз!2~с2!з

 

 

 

 

M l

 

~}— —

0. (4.3.13)

У2 — ЬуХу2Х2Ь2Ху2-}-ЬууХу (ь^-\------ ----+

 

'

«о

'

do

 

Сходимость итерационного процесса для системы

(3. 13) значи­

тельно лучше, чем для (3.9). Используя предложенный способ «последовательных нагружений» в сочетании с методом Нью­ тона—Рафсона, мы получили кривые нагрузка—амплитуда про­ гиба для закритического состояния. Характерная диаграмма на­ грузка—амплитуда прогиба представлена на рис. 4.3.1. На­ грузку, при которой в модели происходит «хлопок», будем счи­

тать предельной

(точка А на рис. 4.3. 1),

а нагрузку, при кото­

рой имеет место

«выхлоп», примем за

нижнюю критическую

(точка С на рис.

4.3.1). Величины qn и qH полностью определя­

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

Графики на рис. 4.3.1 и 4.3.2 построены при следующих значениях параметров: | 1= | 2=1; £з= £б = 0,15; g4 = E5 = 0,l; р! = = 0,55; р2 = 0,04; п* = 4. На рис. 4.3.2 принята амплитуда началь­ ного прогиба /o= ^io + -^2o = 0,l; 0,5; 0,9 для кривых 13 соответ-

U_1.11_______________ 1

J

0

I

2

3

4

5

Рис. 4.3.2. Характерные диаграммы нагрузка—амплитуда прогиба для различных значений начальной ам­ плитуды. 1—3 соответствуют ампли­ туде начального прогиба f0=0,l; 0,5; 0,9.

Рис. 4.3.3. Развертка поверх­ ности цилиндрической оболоч­ ки (заштрихованы зоны пов­ реждения).

9*

Глава IV . Устойчивость оболочек при статических нагрузках

132

ственно. В вычислениях принято допущение (2.1.30).

В этом

случае величины х30, х40 и х50 выражаются через амплитуду на­ чальных прогибов Хю и х20 следующим образом:

-Гзо = х40 = лрфг^го, -*5o = ?lP2^io- (4.3.14)

Для приведенного примера значение верхней критической на­ грузки qB~ 8,2210~2, значение верхней предельной и нижней критической нагрузок для кривой 1 <р= 7,85-10~2 и ^н = 6,98-10-2 соответственно. Из рис. 4.3.2 видно, что в случае больших на­ чальных искривлений с увеличением нагрузки амплитуда моно­ тонно возрастает, т. е. отсутствует «хлопок». Разница между верхней и верхней предельной критической нагрузками неболь­ шая (меньше, чем у изотропной оболочки с такими же парамет­ рами толщины и длины). Этот факт установлен эксперимен­ тально и исследован теоретически [66, 205]. Показано, что для анизотропных оболочек (например, из боропластика) влияние начальных неправильностей в виде отклонений срединной по­ верхности от идеальной формы играет меньшую роль, чем у изо­ тропных. Тем не менее, в оболочках из армированных пластиков существенную роль может играть неоднородность материала. Это обстоятельство можно учитывать, если отнести случайную неоднородность к отклонению поверхности оболочки от идеаль­ ной формы и тем самым приближенно оценивать возможность снижения критической нагрузки.

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

d(w — wо)

823= -

( У у ~ У у Ф ) ) +

ду

 

 

 

 

выражения для уу, уут и w, w0 (3.5) и (3.6),

получаем

823 = [ -

(*5-*5о) + (X i-x10)ttp2]sinaxcos $У-

(4-3- 15)

Соответствующие сдвиговые напряжения огз определяются по

закону

Гука 023= 823623

(здесь 6 2з —

модуль межслойного

сдвига).

Максимальными

напряжениями

огз характеризуются

точки оболочки, где sin ax cos рг/= 1, что соответствует х — Ц 2 и

y = nR/m

(т= 1 ,2 ,...,/г*). Примем А ц = А 22 = 2,5-105 кгс/см2.

Тогда для

нашего примера максимальные сдвиговые напряжения

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