Механика композитных материалов 6 1980
..pdf
|
|
|
Рп(Х) имеют четыре локаль |
|||||||
|
|
|
ных экстремума, а если 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 |
|
|
/ |
|
|
7Г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 г.).