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

Механика композитных материалов 6 1980

..pdf
Скачиваний:
4
Добавлен:
15.11.2022
Размер:
10.98 Mб
Скачать

 

 

 

Рп(Х) имеют четыре локаль­

 

 

 

ных экстремума, а если 6^

 

 

 

=^/г^15 — один экстремум,

 

 

 

так же как и в работе [1],

 

 

 

Наименьшие

значения рн в

 

 

 

зависимости

от

отношения

 

 

 

Е\/Е2 при

Е0= Е[

даны

на

 

 

 

рис.

4.

Для EJE2 = 4, 6,

Ю,

 

 

 

20,

60,

200

соответственно

 

 

 

при Х^4,5, 4,5, 5,0, 5,5, 6,5,

 

 

 

8,5

имеет

место

осесиммет­

 

 

 

ричная потеря устойчивости.

 

 

 

Как видно из рис. 4, для

 

 

 

EJE2= Q0,

200

с

увеличе­

Рис. 4. Кривые неосесимметричной потери устой­

нием параметра

пологости),

чивости

однослойной

ортотропной сферической

критические

нагрузки

не­

оболочки

в зависимости

от отношения EI/E2 (зна­

осесимметричного

выпучи­

чения указаны справа) при G = £ 2/3, V i = 0,25.

вания рв уменьшаются. Для

 

1—15 — номера гармоники п.

Рис. 5. Кривая потери устойчивости для Х=\2,

всех

приведенных

графиков

hi/h0=l при комбинированном действии двух на­

с увеличением X возрастает

 

грузок.

номер гармоники п, которой

 

 

 

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

минимальное

 

 

 

значение Рп.

 

 

 

 

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

х

 

х

Ь= J [p ig i - a ig 2-

z

B2h\) ] dxl2pajx % d x .

0

о

Здесь Рь сц, Вь В2 определяются из следующих линейных краевых

задач:

U n'(Bu B2)=h\(x)-,

U n4Bu B2)= h \ (x)-t

 

 

Bi(k) =Bil {l) = Т2п'В2= Т2п2В2= 0;

 

 

Si (•^PiII+ PiI) — S2— Pi —ФР1+ {х—P)ai = gi (х );

 

 

S 5 ( x a i n + a i 1 ) — 5 7 — a i — (л: — p ) P i — ^ 2 ( л : ) ;

(7)

 

pi (A.) =>.«,! (X) +

4 ^ - a, (X) =0,

 

 

 

Л22

 

 

где p, ф являются осесимметричным

решением,

найденным

при р=рь

а функции

g 1, g2, h*1, /1*2, входящие

в системы

(6), (7), определяются

по формулам работы [3]

(h*i(x) =hi(x), h*2 (x) =h 2 (x)).

Для физических параметров (5)

при h\/h0= 0,2 и Е0= Еi1 имеем:

Х = 9;

Рп= 0,0550;

п = 3;

Ъ= -0,184;

X = 12;

рв= 0,0459;

п = 5;

Ь= -0,241;

Я= 15;

ри= 0,0453;

п= 7;

Ъ= -0,287;

Х = 20;

рн = 0,0407;

п= 6;

Ь = -0,278.

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

дМх

дМху

 

у

<Э2фх

 

дх

+

ду

^

G

12

дР

;

дМу

 

дМху _

~

 

hz ^

_

_ п

ду

+ '

дх

Qv

G

12

dt2

- 0 ‘

Уравнения состояния для цилиндрической ортотропной оболочки в декартовой системе'" координат х, у, z на ее поверхности имеют вид [1] (Л 1=Л2 = 1; *i = 0; * 2 = -1 /Я)

Nx>=h(An1!ЕХ+ Л1122ЕУ) ; Ny=h (Л2222ЕУ+ Л221 ‘Е*);

Мь

12

Nxy = Ny X= h A ' 2'2 (Exy + E Vx)\

 

ax

Qv= hA2323\

 

 

ag

R

(2)

M

 

l

* +

 

h 3

( A 1111

дх

- I л 1122

 

 

\

 

 

 

12 И

+ л

dy

I

 

 

 

(?фх

1

h 3

A1212

/

 

 

•+Л2211-

1 II A К

 

{

ду

дх

ду

дх

 

12 A

 

\

Физические компоненты деформаций выражаются через компоненты вектора перемещений срединной поверхности:

ди 1 /

dw \2

 

dv

1

/

dw \2

w

Ез:=='д х + Т \ 1

х~ )

 

v = ~d^+ ~2

\ ~ d f I

+ т

ди

dv

dw

dw

 

Еху = Еух = ~ду

дх

дх

ду

( 3 )

Ехг —фх‘

dw

 

 

1

dw

v

 

д х

Еyz— фу-! R

ду

R '

 

Физические компоненты тензора жесткости выражаются через техни­ ческие характеристики ортотропного материала:

Л11и=-

Д 2 2 2 2 - ,

Д 1 1 2 2 _ Д 2211 _

£lV21

I - V 12V2I

I - V 12V21

1 — V12V21 ’

 

 

 

 

(4)

£Уу12

. A1212= G l 2 ; A1313= G l 3 ; A2323=G23.

1 - V 12V21

Будем считать, что на поверхности оболочки действует подвижная нагрузка k=kcos (fix—at), где к — параметр нагрузки, k — волновое число, со — круговая частота. Подставляя в уравнения (1) соотношения

(2) — (4), получим уравнения движения оболочки в перемещениях:

k

d2u

 

dw d2w

 

 

 

 

d2v

 

d2w

dw

n

( +

дх

дх2

)

+ л

1122

( дхду

дхду

ду

+ т г ж ) ]

+hA 12,22 (

д2и

 

d2v

 

 

d2w

dw

dw

d2w

ду2

1

дхду

 

 

дхду

ду

 

дх

ду2 ■ ) -

/

 

dw

\

d2w

М -

(

/

 

dw

_

v \

I

d2w

- М13,з ( ^ + _

)

 

_

_

 

^ + _ _

)

(

dxgy

 

 

 

 

 

dv

\

 

 

 

d2u

 

 

 

 

 

(5 )

Rdx ) ~ p h

i [ л 2222

(

d2v

 

dw

d2w

+

1

dw \ +42211 (

 

 

 

 

dw

d2w

 

dy2

+ —------ —

R

dy

 

 

 

 

 

 

dxdy

 

dx

dxdy •)

 

 

'

dy

dy2

 

 

 

 

 

 

 

 

 

 

 

 

dw

\

I 1

dv

 

d2w

 

 

 

 

 

 

/

 

 

dw

 

v \

,+/L41313( Ipx

dx

/

'

R

dx

 

dxdy

)+ h A ^ { ^ y+ —

- - ) x

 

 

 

 

X l iR~

 

w

 

d2w

 

 

 

 

 

d2v

 

 

 

 

 

 

 

 

 

 

( \ R 2

 

dy2 )

1 ~ (,h ~ d f

*=0;

 

 

 

 

 

 

 

(

дфх

 

d2w

 

 

 

 

 

 

 

 

 

 

d2w

 

1

dv \

 

 

hA1313

dx

 

dx2

)

+ M2323 (

4

dy

 

'

dy2

~ ITdiT / +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

2

\

dy

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

d2w

 

 

„ „

/

da

 

 

dv

-+■

dw

\

/

d2w

 

 

 

+ ~R

] }

 

dx2

■+4hA1212 у

dy

 

 

dx

dxdy

) (

V dxdy

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

1

 

dv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

dx

I

 

 

1“

 

L dy

 

 

 

d y!

 

'

R

 

 

 

 

 

 

 

 

1

/

 

 

\21

\

Г

1

 

(

w

 

 

2

 

dv

 

d2w

 

 

 

 

 

+ T

\~dx~ /

J

/

 

 

 

\ Я2~+ Я~

%

 

 

dy<

 

 

 

 

 

 

 

 

 

 

 

 

 

d2w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

 

 

 

 

 

 

ph—— — bp cos (kx — mt) =0;

 

 

 

 

 

h3

/ Лч п -

^

 

 

 

 

d2tyy

\

 

h3 ^41212

/

д 2фх

.

d2,tyy

)

-

~T2

V

 

dx2

 

 

 

dxdy

1

+ _ 12

 

 

'

dy2

 

dxdy

12

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dw У

 

v

 

h3

 

(52фх

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

:+

dx

1

 

 

12

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

if

 

G

 

 

 

 

 

 

 

 

 

h?_ ^/[2222. d ^1/

 

| ^2211

д 2фх

\

 

л3 .41212 /

д 2фх

 

d2^y

)

-

1 2

 

 

 

 

 

 

 

 

dxdy

)

+ 1 2

 

 

 

'

dxdy

 

dx2

 

 

 

— hA2323

(фу +

dw

 

 

\

 

 

уу

h3h3

d2tyy

= 0 .

 

 

 

 

 

 

dy

 

 

/

 

 

G T2~12

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При

действии

подвижной

нагрузки

 

Я cos (kx—cot)

 

на

поверхности

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

Аппроксимирующее выражение для разрешающих функций выберем в следующей форме:

Щх,у) — ^ о(х )"Ь Un(x) cos

пу_

 

пу

R

V(x,y)— ^n(x) Sin R

 

 

 

ny

w{x,y)— wQ{x)+ Wn(x) COS

 

 

 

R

фх(х,1/) — фхО(х) “b фхп(х) COS

ny

 

ny

R

tyvix.y) — ф)/п(х) S in R

На каждом шаге приращения параметра подвижной нагрузки Д), строится периодическое решение А«о, Аип, Аw0, Awn, Avn, Дфх(о)> А'фх(п), А'Фу(п) системы уравнений вида (6), составленной в виде суперпозиции нормированных фундаментальных решений [2]. Запишем на t-ом шаге систему уравнений движения (6) в форме

АХа)=В'х(Хм) AX(i)+fo)(t),

(7)

где В'х — линеаризованный оператор нелинейных уравнений движения

оболочки;

AA'(i) —

решение

системы

уравнений

(7);

=

= AX(i)COS (kx—at) + /’(i)(/); r^)(t)

— вектор

невязок.

иТ-периоди-

 

Система

(7)имеет

Т-периодические коэффициенты

 

ческую правую часть. Необходимо найти ее Г-периодическое решение.

Как известно, для того чтобы решение AX^)(t) такой

системы было

Г-периодическим, необходимо и достаточно, чтобы

 

ЛА'(г)(0) =AX(i)(T) .

(8)

Тогда решение системы (7) будем искать в виде

 

Д^(г)(0 =y(i)°{t) +K (i)(f)a(i),

 

где t/(i)°(0 — решение задачи Коши для системы (7) при нулевых на­ чальных условиях; Y(i)(t) — нормированная фундаментальная матрица решений системы t/(i)= В'х №i))T(i); ct(i)= {ai^, т.. , an(i)} — вектор постоянных, которые подбираются так, чтобы выполнялись условия (8). Таким образом, вектор Щг) должен удовлетворять системе линейных уравнений

(У№(Г )- £ )а « = № > 0(Г).

где Y(Т) — матрица монодромии; Е — единичная матрица. Если мат­ рица Y(T)—E невырождена, то задача построения Т-периодического решения линеаризованной системы динамики оболочки решается одно­ значно. В случае, когда определитель матрицы равен нулю, задача про­ должения решения по параметру становится некорректной. Ее регу­ ляризация достигается построением и исследованием уравнения развет­ вления Ляпунова—Шмидта.

Изложенным методом было исследовано явление распространения бегущих волн в тонкой цилиндрической оболочке при действии подвиж­ ной нагрузки со следующими значениями геометрических и механи­ ческих параметров: # = 50 см; /i = 0,1 см; 6 = 0,01 1/см; G = 981 см/с2; (о= 5000 1/с; Е\ = 1,2-105 кгс/см2; Е2= 1,9-105 кгс/см2; GI2 = 0,25-105 кгс/см2; Gi3 = 0,l-105 кгс/см2; G 23= 0,l-105 кгс/см2; vi2 = 0,16; v2i = 0,129.

Рис. 2.

В результате решения постав­

ленной

задачи

найдены

формы

волн,

установлены

критические

точки,

от

отделены

устойчивые

ветви

неустойчивых.

Анализ

результатов

показал,

что при

действии

подвижной

нагрузки

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

оболочка

обла­

дает

высокой

несущей

 

способ­

ностью.

рис.

1

показана

кривая

На

0

1 —2

зависимости

амплитуд­

ного перемещения w0 от пара­ метра нагрузки X. Она имеет не­ линейный характер. При значе­ нии нагрузки 10 689 кгс/см2 (точка 1 ) определитель матрицы монодромии выродился, и от по­ лученной кривой ответвилось ре­ шение с отличными от нуля зна­ чениями функции w3 (рис. 2). Анализ значений мультипликато­ ров показал, что ветвь 0— 1 рис. 1

устойчива, а

ветви 1 —2 рис.

1 и

1— 3 рис.

2 неустойчивы.

На

рис. 3 представлена форма бегу-

п у

щей волны w0+ w3 cos-^ - для со­

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

потери устойчивости. Можно заметить, что вмятины и выпучины распо­ лагаются в шахматном порядке. По мере продолжения по ветви 1—3 на рис. 2 амплитуда волны возрастает.

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

1.Рикарде Р. Б., Тетере Г. А. Устойчивость оболочек из композитных материалов. Рига, 1974. 310 с.

2.Гуляев В. И., Дехтярюк Е. С. Построение периодических решений нелинейных дифференциальных уравнений. — В кн.: Сопротивление материалов и теория соору­ жений. Киев, 1978, с. 106— 110.

Киевский инженерно-строительный институт

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

УДК 678.067:536.6:539.373

А. И. Бейль, Г Г Портнов, И. В. Санина, В. А. Якушин

УСТРАНЕНИЕ НАЧАЛЬНЫХ ТЕРМИЧЕСКИХ НАПРЯЖЕНИИ В НАМОТОЧНЫХ ИЗДЕЛИЯХ ИЗ к о м п о з и т о в ИЗМЕНЕНИЕМ УГЛА НАМОТКИ ПО ТОЛЩИНЕ*

1. Изготовление изделий типа тел вращения намоткой композитов сопровождается созданием в них системы начальных технологических напряжений. При намотке с малыми усилиями натяжения основную часть их образуют термические напряжения, возникающие при охлажде­ нии изделия после его термообработки и полимеризации связующего. Особую опасность представляют радиальные растягивающие напряже­ ния, вызывающие расслоение намоточных изделий. Теория образования начальных напряжений исследована достаточно полно и изложена в ра­ ботах [1—4]. К настоящему времени предложено несколько способов борьбы с этими напряжениями. Одним из них является намотка с изме­ нением усилия натяжения по рационально выбранной программе [5, 6]. В данном случае напряжения, возникающие вследствие неравномерности натяжения витков, полностью или частично компенсируют температур­ ные. Диапазон относительных толщин, в котором этот метод оказыва­ ется эффективным, определяется степенью радиальной податливости на­ матываемого полуфабриката, от которой зависит падение натяжения в наматываемых витках, и предельно возможной величиной натяжения, т. е. прочностью полуфабриката на растяжение. Этот диапазон может быть несколько расширен сочетанием программированной намотки с по­ слойным отверждением, ужесточающим изделие в радиальном направ­ лении [7, 8].

Другой подход связан с воздействием на процесс охлаждения изде­ лия. Искусственно удлинняя его, добиваются частичной релаксации на­ чальных напряжений [9]. Представляет интерес также метод темпера­ турных градиентов [10, 11], основная идея которого состоит в создании температурного поля с определенной неоднородностью, обеспечивающей подчинение температурных деформаций уравнению совместности. При таком условии можно безопасно миновать в процессе охлаждения тем­ пературу, при которой прочность в радиальном направлении мини­ мальна. Осуществляя процесс полимеризации в условиях неоднородного температурного поля, можно несколько понизить и уровень напряжений в готовом изделии [12]. Однако этот метод имеет ряд ограничений, не позволяющих достаточно точно реализовать наперед заданный закон из­ менения температуры по толщине изделия. Одним из способов уменьше­ ния опасности трещинообразования является опрессовка изделия при охлаждении, также позволяющая без трещинообразования пройти при охлаждении опасное сочетание изменяющихся во времени радиальной прочности и радиальных напряжений [13].

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

* Доклад, представленный на IV Всесоюзную конференцию по механике полимер­ ных и композитных материалов (Рига, октябрь 1980 г.).