Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Механика композитных материалов 5 1982..pdf
Скачиваний:
14
Добавлен:
15.11.2022
Размер:
10.34 Mб
Скачать

У Д К 6 2 4 .0 7 3 .0 0 1 :6 7 8 .0 6 7

В. А. Поляков, В. В. Хитрое

УСТОЙЧИВОСТЬ ПЛАСТИН ИЗ КОМПОЗИТОВ ПРИ СДВИГЕ В ПЛОСКОСТИ ДЛЯ СМЕШАННЫХ ГРАНИЧНЫХ УСЛОВИЙ

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

Большая часть работ, относящихся к расчетной методике по устойчи­ вости тонкостенных анизотропных пластин, выполнена в предположении однородности свойств материала, причем преобладают исследования на случай нормальных нагрузок в плоскости пластины. Расчет при каса­ тельных усилиях, распределенных по контуру, проводился энергетиче­ ским методом в работе [1] и статическим методом в [2]. Обзор результа­ тов этих работ и анализ решения, полученного для полосы при сущест­ венной анизотропии материала, производится в работе [3]. Центральным вопросом в решении этих задач является выбор функции для прогиба при потере устойчивости пластины. Это обусловлено необходимостью рассмотрения ряда гармоник [1] или введения в аргумент тригономет­ рических функций угла наклона косых волн [2] с тем, чтобы обеспечить нетривиальную форму прогиба при изгибе от сдвига.

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

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

Пусть прямоугольная ортотропная пластинка, толщиной h и разме­ рами в плане а и Ъ (рис. 1—б), деформируется касательными усилиями t = hxxy, распределенными равномерно по всем четырем сторонам. Глав­ ные оси упругости параллельны сторонам; стороны х = 0 и х = а жестко

Рис. 1. Кручение заделанного по торцу коробчатого стержня

(а) и схема касательных уси­ лий и размеры пластинки — грани стержня (б).

заделаны, а стороны у = 0, у=Ь свободно оперты. Критическое значение /1ф, при котором пластинка теряет устойчивость, будем определять энер­ гетическим методом. Согласно этому методу [3], следует приравнять по­ тенциальную энергию пластинки при изгибной форме равновесия П ра­ боте касательных усилий А при деформации сдвига, обусловленного из­ гибом пластины. При представлении прогиба в виде функционального ряда критическое значение нагрузки находится из условия вырожден­

ное™ детерминанта следующей системы уравнений [3]:

 

д с г = о;

w

где £/=П —A; Ci — коэффициент при i-й функции прогиба в его разло­ жении в ряд. Для смешанных граничных условий

№ = 0,

dW

=0 при х = 0,х = а; W=0 при у = 0, у=Ь (2)

примем следующую двухчленную форму представления прогиба:

W'=Cn(p„(.v') sin^|- + Cn+icpn+1 (х) sin » (3)

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

dAфп

—Яп4фп-,

л — 1 , 2 , . . .

(4)

dx4

 

 

 

и нормированы в промежутке [0, а]:

 

 

а

 

 

 

| (pn2dx=l;

л = 1 , 2 ,

(5)

о

 

 

 

Запись их имеет следующий вид:

 

 

(cos Хпх — ch Хпх) (sin pin—sh р„) —(sin Хпх — sh Хпх) X

 

Ф„ (x)

X (cos р„ —ch pn)

(6)

] a (sin pn—sh pn)

 

 

 

л = 1,2, 3,.

 

где собственные числа Xn и p„ = aXn определяются из уравнения

 

 

cos pin ch p„= l.

(7)

Условия (2) при выборе прогиба в форме (3) с учетом (6) выполняются автоматически*.

При построении решения укажем [4], что первый корень уравнения 3

(7) имеет значение pi=^ л+0,0176, а для последующих его значений

справедлива с точностью более 0,1% следующая приближенная фор­ мула:

2л+1

(о)

р « = — ^— л; п=2>3’

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

считывается с учетом (3) по формуле

 

 

а Ь

4

а

 

Г f d\V dW

Г

(9)

A = t J J

dxdy = — tCnCn+i J (cpncp^+i —fp/n(pn+i) dx.

* Функции фп имеют размерность [Z.-1/:], так что коэффициенты С„, Сп-и имеют со­ ответственно размерность [UU\.

Прежде чем вычислить интеграл в правой части (9), установим, опи­ раясь на условия (2), (4), (7), некоторые соотношения между интегра­ лами от произведений собственных функций ф„(х) и их производных. Используя уравнение (4) и условие фп|*=о,а = ф/п|*=о,а=0, которое вы­ полняется в силу граничных условий (2), получим после трехкратного интегрирования по частям следующее интегральное соотношение:

а

 

а

 

V J фпф'гЖ^ + Атж4 |

ф»+1ф'п<*Х= [ф,,яф/'п+1]® .

(Ю)

0

0

 

 

Составим также очевидное в силу условий (2) соотношение

аа

J ф п ф ' п + l ^ + J

ф п + 1 ф ' п ^ Л ; = 0 .

( 1 1 )

0 0

Решая совместно (10), (11) относительно интегралов, получим

Г

,

,

-

Г ,

.

[ф//пф//п+1]а

(12)

J фпф n + \ d x =

J

ф n * f n + l d x =

------ т— .

«

 

 

 

п

 

 

Ап

+1

 

Расчет граничных значений вторых производных собственных функ­ ций (6) с учетом (7), (8) сводится к соотношениям

//

2JW2

„ . .

2А.П2 cos

sh p „ -sin

р„ ch рп

 

Ф п (0) = ------

; ф „ (а) =

---у ---------------------a

;--------------------------------

 

 

 

Уа

 

sin |.iri—sn цп

 

-

2Ъп2 ( - l) n sin p „ - tg р„

= ( - I )t .n2 K

2

(13)

 

уa

sin p„ —(—l) n tg p„

Уа,

 

При записи последнего выражения в (13) было использовано следующее из анализа (7) соотношение sh р„= ( —1)" tg р„; л = 1 ,2 ,... Таким обра­ зом, окончательная запись выражения (12) имеет вид

о

и работа касательных усилий согласно (9) выразится формулой

64 Яп2Ап+12

с пс 71+1 t.

За Кп2—Я.п+12

Потенциальная энергия изгиба при потере устойчивости пластины имеет запись

аb

П = - М

J Ф (x,y)dxdy,

(14)

z о

о

 

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

Ф= £>1

d2W дПУ

+ 4Dh

(15)

+ 2DiVi2 дх2 ду2 + 7)2

где цилиндрические жесткости определяются через упругие константы материала в главных осях упругой симметрии:

'

/г3

Ег о

/23

(16)

 

D\,2= г* ----- :— ;

7)/1=— -G 12; D^=D\v\2 + 2Dk.

 

12

1—V12V21

12

 

Здесь h — толщина пластины. В целях упрощения последующих резуль­ татов последним выражением в (16) приведена комбинированная жест­

кость, которая в случае изотропии материала вырождается, так что

^3 = ^1,2- Подставив (3) в (15), при интегрировании по (14) используем соот­

ношение (5) и следующие из него и (2) интегральные зависимости для собственных функций:

 

 

 

 

а

 

 

 

 

 

а

<pn2dx=Xn4;

п=1 ,2 ,...

 

 

 

 

 

 

 

J

(p"2ndx= Xn4 j

 

 

 

(17)

 

 

 

 

о

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

а

ф'2ndx= -

а

ф„ф,/„^=J

 

 

 

 

 

 

 

 

 

J

 

Л „2Хп2-

;

п= 1 ,2 ,...

 

 

(18)

 

о

 

о

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

г,

о

 

моч

параметр

 

 

COS(X„-ch(X„

 

sin (An

 

 

от-

Введенный в (18)

хп= —--------- ;----- =

—;— гг--------------

 

 

 

 

 

 

 

 

 

 

 

 

sin pn -sh р.„

( —1)“ —cos

 

 

 

личается от единицы менее чем на 0,1% при

2.

 

 

 

 

В результате, учитывая соотношения (17), (18), получим интегриро­

ванием по формуле (14) следующее выражение для потенциальной энер­

гии пластины:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

I f ^

 

D i b

 

 

It D 2n 4

, / ) 3 ( р п 2Х п 2 - 2 р , п > С п ) п 2

 

 

 

п =

2 VС4

1

?

•*" + 1 5 Г

+ -----------М,------------

 

 

 

 

 

 

 

Цп+14 +

16Z)2jt4

 

4£>з(|хп-нп+г

2|in+i^n+i) л-

 

 

 

 

 

 

 

~ 2 Ь ~ +

 

аЧ>

 

 

 

 

Приравнивая детерминант системы двух уравнений (1)

при i=n,

л+1

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

сдвигающего усилия при потере устойчивости пластины:

 

 

 

 

 

 

 

 

 

 

tкр —

 

3

Цп+14 — Цп4 l / D l D j

 

 

 

 

 

 

 

 

 

 

 

128

 

р„2р„+12

Ь'2

 

 

 

 

 

X VivD2 С3

 

 

 

 

 

 

 

 

 

П

 

' Di

 

 

С

 

yDiD2 -I

х

 

 

 

 

D

1 Цп4

|

 

 

у

Д 2 |

2 я 2(р п2Хп2- 2 р п К п )

D i

1

 

X

ГDl ]

/

Ц

п1+ ‘41

6

л

'

° 21

/ |

8

n 2

( P n +

l2K n + l2

-

2

p n + l X

'

D2

С3

Л

 

Di

 

 

С

T/nI/D1.DnJ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19)

где для относительной длины принят параметр С—а/Ь.

Следует иметь в виду, что в зависимости от порядка собственных функций п критическая нагрузка, найденная по формуле (19), будет со­ ответствовать различающимся между собой изгибным формам потери устойчивости пластины. При я=1 прогиб W изменяется вдоль длины пластины согласно двум функциям — ф] и ф2, первая из которых описы­ вает одну полуволну, вторая — две с нулем посередине длины. С увели­ чением п число полуволн пропорционально возрастает.

Для одной и той же нагрузки ^кр при соответствующем соотношении размеров пластины в плане могут существовать две смежные формы по­ верхности при потере устойчивости. Значение параметра С=Сппр, при котором происходит изменение формы поверхности, соответствующей функциям ф„, фп+ь к форме поверхности, соответствующей функциям

фп+ь фп+2, определяется из уравнения

 

^ t

^

4)2 Рп (С) Qn+l (С) = - - n+2V n+;4)2 Рп+i (С) Q„+2(С),

(20)

Р п

JIn+1*

(ln + lV n + 24

 

где полиномы имеют вид

Р п { С ) = п * ~ \ 1 — - С * +

2 л 2( ц п2х п2 — 2 ц пХп ) D 3

д а т

1 Dt

Qn(C) = 16я4 ] / - ~ с ^ + 8л2 (цп2хп2 —2Ц п Х п )D ifD\D2

Между двумя значениями С„,_I"P и Сппр находится такое значение С*„, при котором критическая нагрузка, приводящая к потере устойчивости пластины по форме, описываемой двумя собственными функциями срп и фп+1, будет минимальной*. Это значение относительно параметра опреде­

ляется как С*п= Ух, где х — наименьший вещественный корень урав­ нения

4 / D \ 1 6j.ln^ “1“ P'n+l^ Z?32 (Ц?1^%712 2(Ll7lX?l) ( Ц7? + 1“%7? + 1 2 (ln + lX77-f 1) \ 2

* _ \ ~D 2

+ z v

^

/ * "

D3DI

 

 

 

4 n W [ (

2|XnX?i) [In-f-14 ”T 4 (|l7i + l“%?i+ l“

2(.ln + iX n + l ) |Я ^ ]Х

 

3

ZV

 

 

16я8

Z)22Цп4Р'?г+14 —0.

 

Сравнение расчетных значений критической нагрузки с результатами, полученными при использовании приближенных формул в [1], показы­ вает, что выражение (19) дает несущественное ее отличие (порядка 3—5%) в случае поперечного армирования пластинки, т. е. при укладке всех волокон в направлении оси у. Этот случай соответствует предпо­ сылке о существенной анизотропии материала, принятой в [1]. Вырож­ дение (19) на случай изотропии материала дает завышенные оценки для ^кр по сравнению с приведенными в [5] (отличие порядка 25%). Для по­ следнего случая, очевидно, недостаточно двухчленного представления прогиба в форме (3). Таким образом, выражение (19) можно рекомендо­ вать для качественного анализа изменения критической нагрузки при потере устойчивости ортотропной пластинки от сдвига.

На рис. 2 представлены кривые изменения приведенных критических напряжений пластин из углепластика в зависимости от относительной длины С= а/Ь для различной анизотропии свойств материала. Построе­ ние расчетных кривых приводилось при минимизации тКр по порядку соб­ ственных функций п. В случае продольно-поперечной укладки волокон в углепластике (см. рис. 2—а) изменение анизотропии упругих свойств характеризуется параметром й2, равным относительному содержанию волокон, уложенных в направлении оси у. В этом случае монотонное из­ менение цилиндрических жесткостей таково, что между их граничными значениями имеется следующая связь: D\(й2 = 0) =£>2(й2 = 1) и £>3(й2) =

= const. В случае косоугольной

укладки волокон под углами ±0 (см.

рис. 2—б) для цилиндрических

жесткостей справедливо соотношение

Z)I (0 = O)=D2 |0 = ту"| = £ I № = 0); н о D3(0) ^const. С увеличением доли

волокон в поперечном направлении у или с увеличением угла 0 критиче­ ская нагрузка возрастает в несколько раз. Так, изменение Ткр при С < 3 от смены укладок с продольной на поперечную составляет ~400% . Ана­ логичная тенденция имеет место и для случая свободного опирания по всему периметру пластины (см. [2]). Отметим, что в обоих случаях (см. рис. 2) минимизация нагрузки при переборе порядка функций п обнару-

* В случае n= 1 левая граница С0пр условно принимается равной нулю, а не нахо­ дится из уравнения (20).

Рис. 2. Влияние относительной длины с=а)Ь пластин и схемы укладки однонаправлен­ ного углепластика (£ц= 13,2 • 105 кгс/см2, £ ± =1,22*105 кгс/см2, “ 0,25 • 105 кгс/см2,

v llJL = 0,33) на приведенные критические напряжения: а — ортогонально армированный материал с различной долей й2 поперечных волокон; б — косоармироваиный материал с

различным углом армирования ±0 по отношению к оси х.

1—7 — значения п.

Л2 = 0 (I);

0,2 (II); 0,4 (III); 0,8 (IV); 1,0 (V). 0 = 15 (I); 30 (И);

45 (III); 60 (IV) и

75° (V).

живает немонотонное изменение Ткр по степени анизотропии упругих

свойств материала

(см. кривые на рис. 2 при йг= 0,8 и

0= ±6О°).

Сопоставление

критических напряжений

типовой

пластины (h/b =

= 0,01) из косоармированного углепластика

с приведенными в [5] дан­

ными по прочности при сдвиге подобного материала в плоскости арми­ рования 700—1200 кгс/см2 показывает, что в широком диапазоне 0 не­ сущая способность пластин будет определяться их жесткостью.

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

СП И С О К Л И Т Е Р А Т У Р Ы

1.Балабух Л. И. Устойчивость фанерных пластинок. — Техника воздуш. флота, 1937, Ко 9.

2.Секерж-Зснькович Я. И. К расчету на устойчивость листа фанеры, как анизотроп­ ной пластинки. — В кп.: Тр. ЦАГИ, 1931, № 76.

3.Лехницкий С. Г. Анизотропные пластинки. М., 1957. 463 с.

4.Канторович Л. В., Крылов В. И. Приближенные методы высшего анализа. М.; Л. 1962. 708 с.

5. Работное Ю. Н., Данилова Н. Н., Полилов А. Н., Соколова Г. В., Карпсйкин И. С., Вайнбсрг М. В. Исследование прочности намоточных эпоксидных угле- и стек­ лопластиков при кручении, растяжении и поперечном изгибе. — Механика полимеров, 1978, № 2, с. 219—225.

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

Поступило в редакцию 22.01.82

АН Латвийской ССР, Рига

 

М Е Х А Н И К А К О М П О З И Т Н Ы Х М А Т Е Р И А Л О В , 1982, Лг° 5, с. 839— 843

У ДК 624.074.001:678.067

В. П. Трошин

ВЛИЯНИЕ ПРОДОЛЬНОГО РАССЛОЕНИЯ в СЛОИСТОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКЕ НА ВЕЛИЧИНУ КРИТИЧЕСКОГО ВНЕШНЕГО ДАВЛЕНИЯ

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

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

рим половину поперечного сечения оболочки

(рис. 1).

 

Уравнения

нейтрального равновесия цилиндрических панелей I, II,

III имеют вид

[4]

 

 

 

 

TUA + T2U2 + N ,R -' = Q-

1

+

 

 

МиЛ+М21,2- ^ = 0;

/

чГ

U

N \11+ N2,2 Т22^?2_I + *7n —0.

Здесь (.. ,) { = д/д1й gb g2 — продольная и окружная координаты; R 2 =

= /?, R1= 00; R — радиус цилиндрической оболочки. Символ 1,2 обозна­

чает, что невыписанные формулы получаются из приведенных переста­ новкой индексов 1 и 2. Подставляя значение нормальной составляющей внешней поверхностной нагрузки qn— q^2R (где к2 — изменение кри­ визны нейтральной поверхности в окружном направлении) в последнее уравнение системы (1), получим уравнения устойчивости оболочки при

действии внешнего давления q.

 

 

 

В смежной форме равновесия связь между усилиями и моментами

Tij,

N u Mij, с одной стороны, и перемещениями щ,

w, с другой, задается

с помощью соотношений теории тонких ортотропных оболочек [4]

 

^ 1 1 = ^ П [ ^ 1 , 1 -\~ w R i 1+ V2 1 (^2 , 2 + w R 2 !) ] i

 

 

Т\2=*В\2[u2,\ + ^1,2+ "7Г h2Rz~l (-W, 12+^2,1^?2_1 + ^1,2^1_I) ] i

(2)

 

О

 

11

~ ^11 [^,11 —^l.l^l-1 + V21 (w,22~ U2>2R 2~ l) ] ’>

 

 

12=^ i2 (ДО,12~~ U2,\R2~X— U\i2R l_1) •

Здесь Bi^EihKl-VijVji); B ij=Gijh; Dij=— B ijh2; Eh Gijt v0- — модули

Интегрирование системы дифференциальных уравнений (7) с гра­ ничными условиями (8) проводим методом Кутта—Мерсоиа с промежу­ точным ортонормированием векторов решения [5, 6]. Значение крити­ ческой нагрузки находится из системы алгебраических уравнений (5). Алгоритм, реализующий программу поиска критической нагрузки, раз­ работан для БЭСМ-6.

Параметрический анализ проведем на примере круговой цилиндри­ ческой оболочки (jR/Am = 100) из композитного материала GI2/£' = 0,25 (Е\ = Е2, vi2 = V2i =v = 0,25). Величину критического внешнего давле­

ния для оболочки длиной I определим по формуле [7]

 

#о=

^ п^223(1—V12V21).

(9)

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

Вначале рассмотрим расслоение, расположенное на нейтральной по­ верхности оболочки. Параметр h* = hI/hm i определяющий место распо­ ложения трещины по толщине оболочки, в этом случае равен А* = 0,5. На рис. 2 сплошными линиями обозначены зависимости критического параметра q*= qlq<n от параметра ширины трещины l*2 = h ! ^ R при раз­ личных значениях параметра /*= ///?, определяющего длину оболочки и расслоения. Здесь q01— величина критического внешнего давления обо­ лочки длиной /*= 1, определяемая по формуле (9); /2, 1\ — ширина и длина расслоения = Штриховыми линиями на рис. 2 обозначены критические нагрузки нереализующихся форм потери устойчивости.

Как видно из рис. 2, изменение длины оболочки существенно влияет на абсолютные значения критического внешнего давления. Однако об­ щий характер зависимостей q*= q*{l*2) при варьировании параметра длины Г сохраняется. Все эти зависимости при изменении /*2 от нуля до /*20 = 0,015—0,02 имеют горизонтальный участок, в пределах которого критический параметр q*= const. Следовательно, значение 1*2о опреде­ ляет предельную ширину допустимого расслоения. Меньшие по разме­ рам расслоения не снижают первоначального уровня критической на­ грузки q. Большие по размерам расслоения приводят к снижению критического внешнего давления и к локализации общей формы волно­ образования.

Общая потеря устойчивости происходит по одной из двух локальных форм: симметричной (рис. 3—а) и кососимметричной (рис. 3—б) (име­ ется в виду симметрия относительно центра расслоения). Изменение от­ носительных прогибов w* = w/w0 на рис. 3 показано для четверти попе­ речного сечения оболочки (w0 — наибольший прогиб). Исследование форм потери устойчивости показывает, что при симметричной форме мо-

Рис. 3. Общие формы потери устойчивости:

а — симметричная (/*2= 0,12; /* = 2);

б — кососиммстрнчная

(/*2 = 0,05; /* = 2).

Рис. 4. Зависимость критического параметра от размеров трещины, различно распо­ ложенной по толщине оболочки: /г*= 0,05 (/); 0,1 (2); 0,15 (5). /*=1.

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

Для рассмотренных ранее оболочек длиной /* = 1 на рис. 6 представ­ лена зависимость параметра предельно допустимых расслоений /*20 = = l2o/2 nR от параметра Л*, характеризующего место расположения тре­ щины по толщине оболочки. Сплошная линия соответствует расслое­ ниям, приводящим к общей кососимметричной форме. Можно считать, что эта линия симметрична относительно вертикальной оси, проходящей через значение /г* = 0,5. Следовательно, расслоения, расположенные сим­ метрично относительно нейтральной поверхности оболочки в интервале 0,35—0,65, равноопасны. В общем случае на критический параметр и размер предельных расслоений влияет кривизна поверхности. По этой причине расслоение, расположенное ближе к внутренней поверхности оболочки, в большей степени снижает величину критического внешнего давления, чем расслоение, находящееся ближе к наружной поверхности. В соответствии с этим различаются и предельные размеры расслоений, обозначенные на рис. 6 штриховыми линиями. Эти размеры определены на основе смешанных форм потери устойчивости. Штрихпунктирные ли­ нии на рис. 6 соответствуют расслоениям, расположенным вблизи на­ ружной или внутренней поверхностей оболочки. Предельный размер та­ ких расслоений установлен исходя из возможной потери устойчивости по местной форме. Так как при местной форме следует ожидать неустой­ чивости только одного тонкого слоя, то потери несущей способности всей оболочки в этом случае не происходит (при условии нераспространения трещины). Следовательно, найденные на основе местной формы потери устойчивости предельные размеры расслоений (штрихпунктирные ли­ нии) не являются ограничениями по несущей способности.

Напомним, что все результаты получены при значении параметра толщины оболочки R/hJ1I= 100. Как показывают дополнительные рас­ четы, изменение этого параметра в широких пределах (R/hm = 30-^300) не меняет общего характера представленных зависимостей.

СП И С О К Л И Т Е Р А Т У Р Ы

1.Болотин В. В., Зебельян 3. X., Курзин А. А. Устойчивость сжатых элементов с дефектами типа расслоений. — Проблемы прочности, 1980, № 7, с. 3—8.

2.Трошин В. П. К устойчивости цилиндрических оболочек с расслоениями. — Ме­

ханика композит, материалов, 1981, № 4, с. 729—731.

3. Андреев JI. В., Ободан Н. И. Задачи устойчивости цилиндрической оболочки с переменной жесткостью при внешнем давлении. — Прикл. механика, 1968, т. 4, вып. 12, с. 82—88.

4.Новожилов В. В. Теория тонких оболочек. Л., 1962. 430 с.

5.Годунов С. К. О численном решении краевых задач для систем обыкновенных

линейных дифференциальных уравнений. — Успехи мат. наук, 1961, т. 16, вып. 3(39),

с.171— 174.

6.Кармиишн А. В., Лясковец В. А., Мяченков В. И., Фролов А. Н. Статика и ди­

намика тонкостенных оболочечных конструкций. М., 1975. 376 с.

7. Королев В. И. Слоистые анизотропные пластинки и оболочки из армированных пластмасс. М., 1965. 272 с.

Поступило в редакцию 28.09.81