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

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

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

новленным фактом при введении наполнителя в мо-

 

номер является повышение начальной скорости по­

 

лимеризации, наблюдающееся в том случае, когда

 

между ними имеет

место взаимодействие,

например,

 

с образованием водородных или химических связей.

 

Такое взаимодействие

может быть связано также и

 

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

 

полненной системе.

 

 

 

 

Нами были поставлены опыты по изучению ки­

Влияние модификации поверхности

нетики полимеризации

бутилметакрилата

(БМА) в

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

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

ции бутилметакрилата (инициатор

действующего с мономером. В качестве наполнителя

ПБ 0.060 моль/л, 7'=75°); / — без

наполнителя; 2 — 5% А-ТМС; 3 —

использовали аэросил

калушский — А -175, немоди-

5% А-175; 4 — 5% А-МТЭС; 5 — и

фнцировапный и модифицированный триметилхлор-

присутствии метакрилатметилтри-

этоксисилана.

силаном (А-ТМС)

и метакрилатметилтриэтоксисила-

 

ном (А-МТЭС). Аэросил А-ТМС получали по методике [3]. Аэросил А-МТЭС приго­ тавливали обработкой А-175 6% водно-спиртовым раствором метакрилатметилтриэтоксисилаиа, взятым в трехкратном количестве к стереохимическому на гидроксильные группы поверхности А-175. Обработку проводили в течение 3 ч при 20° С, 1 ч при 60° С, после чего аэросил отделяли и высушивали.

В случае аэросила А-ТМС силанольные группы поверхности замещали триметнлснлоксигруппами, в А-МТЭС — метакрилатметилсилоксигруппами. В аэросиле А-МТЭС мо­ дифицирующие группы химически связаны с поверхностью. Методом инфракрасной спектроскопии было установлено, что при наполнении БМА аэросилами А-175 и А-ТМС образуются водородные связи [4]. Элементный анализ образцов аэросилов показал, что в случае А-ТМС только 50% силанольных групп аэросила были замещены, а в А-МТЭС были замещены все силанольные группы. Характеристики использованных в полимери­ зации аэросилов приведены в табл. 1 .

На рисунке приведены кинетические кривые полимеризации ненаполненного БМА и БМА, наполненного 5% аэросила разных видов — А-ТМС, А-175 и А-МТЭС. Из кине­ тических данных видно, что введение-в мономер модифицированных и немодифнцированного аэросилов приводит к увеличению скорости полимеризации БМА. При этом наибольший эффект достигается в случае аэросилов А-МТЭС и А-175; в меньшей сте­ пени влияет аэросил А-ТМС, в котором лишь 50% силанольных групп остались свобод­ ными, т. е. способными к взаимодействию с мономером с образованием водородных связей.

Контрольный опыт изучения кинетики полимеризации БМА при введении в него

метакрилатметилтриэтокснсплапа в соотношении, равном тому, которое вводилось вместе с 5% А-МТЭС, показал отсутствие какого-либо влияния этой добавки на ско­ рость процесса (см. рис. — кривая 5, совпадающая с кривой 1). Таким образом, уско­ ряющее влияние этого мономера на скорость полимеризации БМА имеет место лишь в том случае, когда он химически связан с силанольными группами поверхности аэросила.

 

 

 

 

 

 

Табл. 1

 

 

Характеристика поверхности аэросилов

 

 

 

 

 

Количество

Содержание углерода, %

Аэросил

 

Поверхностные группы

 

 

 

групп,

%

рассчитано

найдено

 

 

 

 

 

А-175

 

—ОН

100

 

 

 

 

 

 

А-ТМС

 

—О—Si— (СН3)з

50

 

3,50

1,72

 

 

—ОН

50

 

А-МТЭС

- ° \

 

100

 

8,95

8,50

—О—Si—СН2—О—С—С =СН 2

 

 

0 /

1

 

 

 

 

 

U

СН3

 

 

 

 

 

 

Табл.

2

 

 

 

Табл. S

Изменение

количества

привитого

 

Влияние химии поверхности

наполнителя

полимера в

зависимости от степени

(5% ) на механические свойства

превращения при полимеризации

 

 

полибутилметакрилата

бутилметакрилата с А-МТЭС

 

 

 

 

 

 

 

Количество

привитого

 

 

Проч­

Относи­

 

Степень

пзли-БМА, % в расчете

 

 

ность

тельное

Внешний вид

иа исходный аэроспл

 

Наполнитель

на раз­

превращения

 

 

 

удлинение,

пленок

%

мономер

аэросил

 

 

рыв,

%

 

 

 

 

кгс/см2

 

 

23

0,11

2

 

Без напол-

70

420

Прозрачная

35

0,45

8

 

нителя

85

380

 

 

А-175

Опалесцн-

60

2,00

40

 

 

А-МТЭС

400

102

рующаи

90

26,00

500

 

Мутная

 

А-ТМС

100

450

Прозрачная

 

 

 

 

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

Были исследованы свойства композитных полимеров, полученных полимеризацией в присутствии описанных наполнителей. При полимеризации БМА в присутствии А-175 на поверхности наполнителя не образуется привитого полимера. Экстракция бензолом приводит практически к полному удалению полимера с поверхности аэросила. Полимер же, полученный в присутствии А-МТЭС, содержит нерастворимую часть, не экстраги­ руемую бензолом и днметнлформамндом с поверхности наполнителя. Как видно из табл. 2, количество привитого полимера увеличивается с глубиной полимеризации. Та­ ким образом, использование метакрилатметнлтриэтоксисилана для модификации аэросила позволяет получать полимерные композиции, в которых часть полимера химиче­ ски связана с поверхностью модифицированного наполнителя вследствие сополнмерпзацпп основного мономера с привитым к поверхности мономером-модификатором.

Все исследованные наполнители изменяют свойства исходного полимера. Механи­ ческие свойства и внешний вид пленок из полученных композитных полимеров приве­ дены в табл. 3. Из таблицы видно, что наибольшее влияние на прочность, увеличивая ее в четыре раза, оказывает наполнитель, химически прививающийся к полимерной матрице. Вызывает интерес также влияние на свойства композитных пленок наполни­ теля А-ТМС, в котором часть силанольных групп замещена триметилсилокснгруппамп, увеличивающими совместимость аэросила с полимерной матрицей за счет образования переходного слоя у поверхности наполнителя. Введение А-ТМС приводит к возрастанию прочности и эластичности композитного полимера, а также прозрачности по сравнению с полимером, наполненным ^модифицированным аэросилом А-175.

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

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

1.Морозова Е. М., Елисеева В. И. О структурных превращениях при полимериза­ ции наполненных мономерных систем. — В кн.: Гетерогенные полимерные материалы. Киев, 1973, с. 76—83.

2.Морозова Е. М., Асламазова Т. Р., Елисеева В. И. О влиянии наполнителя па стадию инициирования радикальной полимеризации. — Высокомолекуляр. соединения. Сер. Б, 1978, т. 20, № 11, с. 859—862.

3. Никитин 10. С. Способы приготовления некоторых адсорбентов и носителей для молекулярной хроматографии. Химическое модифицирование поверхности кремния. — В кн.: Экспериментальные методы в адсорбции и молекулярной хроматографии. М., 1973, с. 446—447.

4. Schulz D., Kunath D., Morozova E. M. Infrared spectroscopic study of methacry­ lates on aerosil. — React. Kinet. Catal. Lett., 1979, vol. 12, N 2, p. 113— 117.

Институт физической химии АН СССР, Москва

Поступило в редакцию I3.02.8C

Механика композитных материалов. 1980, № 5, с. 924—927

УДК 624.074:678.067:539.3

И. Г Стрельченко

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

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

Методами, изложенными в работах [1, 2], были изготовлены три конструкции пере­ секающихся под прямым углом цилиндрических оболочек п плоские образцы для опре­ деления упругих постоянных материала. Стеклопластпковые составные оболочки и об­ разцы выполнены на основе стеклонити алюмобороснликатного состава БС-6-18Х1Х6 на парафиновом замаслнвателе, пропитанной эпокспфенольной композицией ЭФ32-301-16. Технология изготовления образцов и сопрягающихся оболочек сохранялась неизменной с целью получения сопоставимых результатов.

Сопрягающиеся оболочки с соотношением слоев армирования в кольцевом и про­

дольном

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

параметры:

радиус

сре­

динной

поверхности, толщина п длина основной оболочки II

были равны

0,119,

0,003

и 0,6 м соответственно, а оболочки-патрубка I — 0,051, 0,003 н 0,17 м.

 

 

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

Методы электротензометрии позволяют провести такие измерения с сохранением стабильной чувствительности во всем диапазоне изменения деформаций и при повтор­ ном приложении нагрузки. Тензодатчики для измерения деформаций в окрестности со­ пряжения оболочек должны обладать малой базой и низкой поперечной чувствитель­ ностью при достаточно высоком сопротивлении (не менее 40 Ом) для использования современных регистрирующих приборов. Всем этим требованиям вполне отвечают тен­ зодатчики сопротивления типа 2ПКБ-10-100ХВ, электрические сигналы которых выво­ дились на статический измеритель деформаций СД-1.

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

Эта симметрия и обусловила наклеивание эпоксидной смолой холодного отвержде­ ния тензодатчиков на внешнюю поверхность основной оболочки. Датчики располагались Т-образно на расстоянии 0,001 м от линии пересечения при ф= 0, 30, 45, 60, 9 0 ° ,..., причем их направления совпадали с введенными на основной оболочке полярными коор­ динатами (р, ф). На рис. 1 приведена схема расположения тензодатчиков (вид со стороны оси патрубка). Показание каждой пары датчиков дублировалось другой парой, симметричной ей относительно продольной и поперечной плоскостей. Для обнаружения

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

также были наклеены тензодатчики. При проведении экспериментов на

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

крышками и опорами. Испытания проводились с помощью машины RH-30 «Шимадзу»

(Япония).

Оболочканагружалась

ступенями осевым сжимающим усилием р=

= 1800-9,81

Н. Накаждой ступени

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

усилие—деформация.

Упругие характеристики материала конструкции, определенные с помощью плоских образцов, приведены в работах [3, 4]. Напряжения подсчитывались по формулам

Ер

(ep+|i(pe<p); Оф=

Е<р

op—

(бф~Ь|лрбр),

1 |-1р|-1ф

1 [^р Р-ф

где Ер, £ ф, |.ip, рф — модули

упругости и коэффициенты Пуассона материала конструк­

ции пересекающихся оболочек в полярной системе координат, определяемые формулами преобразования упругих постоянных [5].

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

ходит однородное напряжение 0 = —80,5-9,81 • 104 Н/м2.

 

 

 

 

 

 

Табл.

1

ф°

ер • ю»

Еф - 10»

(Тр •

(9,81 X

пф •

(9,81 X

 

ХЮ4

Н/м2)-'

Х104

Н/м2)-'

 

 

 

 

 

0

0,20

0,41

 

65

 

104

 

30

0,19

-0 ,2 8

 

7

 

- 3 2

 

45

-0 ,6 5

-0 ,4 8

-1 6 1

- 1 5 3

 

60

-0 ,4 6

- 1,10

-1 9 4

- 2 4 8

 

90

-0 ,9 8

-1 ,6 0

- 3 0 6

- 4 4 5

 

 

 

 

 

 

 

Табл.

2

Расстояние от линии

ст0 • (0,81 •

104 Н/м2)-'

сгф • (9,81 104 Н/м2)-1

 

 

 

 

 

 

 

пересечения, м

LL*

 

 

LL*

 

 

 

 

К К **

К К **

 

 

0,005

65

- 3 0 6

104

- 4 4 5

 

 

0,021

31

- 2 6 2

49

- 2 5 9

 

 

0,037

15

- 2 2 4

46

13

 

 

0,053

14

- 1 8 3

54

19

 

 

0,069

- 5 7

- 1 8 0

51

29

 

 

0,075

- 8 2

 

 

40

 

 

Вдоль образующей LL.

Вдоль окружного направления КК.

Ё табл. 2 приведены Напряжения (?р И аФвДоЛь направлений LL при ф=б° и КК при

ф=90°. Согласно экспериментальным данным табл. 2, наибольшие напряжения наблю­ даются на линии пересечения оболочек. При удалении от линии пересечения возмущен­ ное напряженное состояние затухает, и это затухание вдоль направления LL происхо­ дит несколько быстрее, чем в окружном направлении КК-

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

1.Стрельченко И. Г Изготовление и прочность пересекающихся под прямым углом ортотропных стеклопластиковых цилиндрических оболочек. — Механика полимеров, 1976, N° 4, с. 743—744.

2.Игнатов И. В., Малицкая И. Г., Попенко В. Г., Работное Н. А., Стрельченко И. Г., Юрьев С. В. Исследование механических характеристик образцов из стеклопластика.—

Вкн.: Устойчивость и деформативность элементов конструкций из композиционных ма­ териалов. Киев, 1972, с. 108— 112.

3.Игнатов И. В., Стрельченко И. Г., Юрьев С. В. Статистические характеристики

механических констант стеклопластиков. — Механика полимеров, 1972 N° 6 с. 1125— 1128.

4.

Игнатов И. В., Работное Н. А. Вероятностные характеристики модуля сдвига

стеклопластиков. — Механика полимеров,

1973, № 1, с. 160— 162.

5.

Амбарцумян С. А. Теория анизотропных оболочек. М., 1961. 384 с.

Институт механики АН Украинской ССР,

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

Киев

 

 

Механика композитных материалов, 1980, № 5, с. 927—929

УДК 624.074.4.001:539.4:678.067

В. И. Кобелев, В. В. Кобелев, В. А. Потопахин

ОБ ОДНОМ ВАРИАНТЕ УРАВНЕНИЙ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ МНОГОСЛОЙНЫХ ПЛАСТИН И ОБОЛОЧЕК

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

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

Рассмотрим многослойную оболочку, состоящую из чередующихся «жестких» слоев — наполнителя — толщиной t и мягких слоев — связующего — толщиной с.

Введем

триортогональную

систему координат

х, у, z. Слои

с координатой z > 0 будем

считать

положительными

(i= 1, 2, 3 , . . . , п),

с координатой

z < 0 — отрицательными

( / = 1 , 2 , 3 , . . . , л ) .

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

определенному закону (линейному, кубической параболы и т. д.).

Запишем перемещения любого произвольного слоя оболочки через перемещения со­ седних слоев. Перемещения i слоя связующего выражаются через перемещения слоев наполнителя i+ 1 и / —1 :

(

t dwi+i

t dWi-i

 

 

( 1)

/

t d w i + i

\ r

/

if

^ ,

“,= \ “i+,+7 —д Г

)

v

T

a» / “'

W i = W i + i F 2i + W i - lF l i .

Перемещения слоя наполнителя i 1 определим через перемещения слоев связующего

1—2 и i:

t

 

C

/

dfli

 

 

\ 1

 

 

2 sinl ^ 7 “ ^ 7 /

 

 

Г

C

/

d V i —2

 

d W i - 2 \ 1

+

 

 

 

 

 

ax

 

 

 

 

c

/

---------

do>i

\

1

Ui-----sin

I

 

I ^21-1 +

[

 

2

\

dz

 

ду

I

J

Г

c

l

du.i-2

dwi-z

\

1

+ Г - 2+Т 5Ш1

 

*

-

ay Л

( 2)

Здесь u, u, w соответствуют перемещениям слоев в направлении осей х, у, z. Функции Р\и Pii описывают характер деформации слоев. Они должны удовлетворять условиям

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

 

/

г н

1 Г

-

i t - 2 z

1

;-) -

Ч

 

Рц\

=

1

 

 

=

- -(

-1 - +-J

i-

 

 

для слоев наполнителя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

i

1

Г

 

ic—2z

1

)

+ Pi2 1 =

1 P ii’-

 

J

;

 

i =

 

(

1

+

i

 

^

в выражениях

(2)

углы

dwi

dwi

dvj

дщ

 

такими,

что для них

Будем считать

-----, ----- ,

——, —

 

выполняются условия.

 

 

 

 

 

 

дх

 

ду

dz

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvi

dWi

 

dv i

dwi

 

 

 

/

dui

dWi

>\

dui

dWi

dz

дх

1

dz

dx

 

 

 

\

dz

 

dy

J1

'

dz

dy

 

cos

|

dvi

 

dW;i

\

 

(

 

дщ

dwi

\

s i .

 

 

 

dz

 

dx

1

- C0о 1

dz

 

dy

)

 

 

 

 

,

 

/

 

V

 

"

 

 

 

Используя формулы (1)— (3), а также условия равенства перемещений на границах слоев, получим системы уравнений для определения перемещений любого внутреннего слоя i через перемещения двух граничных слоев п и п :

 

/

t

dwi

 

\

 

 

 

 

 

 

 

 

 

t

dWi-2

\

 

 

 

 

( Vi +

 

dx

 

F 2 i

 

 

 

 

 

 

 

 

2

dx

 

) Fu-1=0;

 

 

 

2

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv i_i

 

dW i-1 \

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

dx

) ] Pu-,-V '-,+ [<-«-,+!■ (

 

 

 

 

 

 

 

 

)

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[ °7+*- 7

('

dvn+3

 

dwn+3

 

) I F2n+2—Vn+2+

 

 

 

 

VZ

 

OX

 

 

 

 

 

 

J

J

 

 

 

 

 

 

 

 

,

Г -

 

 

c l

l

dvZ+i

 

dWn+i

\

1

J

_

 

 

 

 

 

+ L”*+1+r

 

- -ai- - -i

P

 

)

f ‘”+2=0'

 

 

/

_

t

dw~+2 \

 

 

_

 

_

 

[

 

_

 

t

dw~ \

_

 

 

Г "+2+7

 

a*

 

 

 

 

 

 

 

 

 

 

7 -1 7 - j fi

 

 

/

/

dz&)n

\

F2n- 1

 

 

/

un -

 

 

t

dw n - 2

\

Fin -i—0;

 

I

Wn H —

 

-

 

I

— Un- 1 + I

2 ----------------- ■)

 

\

2

ду

 

I

 

 

 

 

 

\

 

 

 

 

2

 

dy

 

I

 

 

 

 

Г

 

 

 

c l dun- 1

-

dwn- 1 \ 1

 

 

 

 

 

 

 

 

L“”-‘“7 l

 

a2

 

ay

 

I J f2"-2- “- 2+

 

 

 

 

 

 

 

 

c / dun- 3

дшп-з

\

I

 

 

 

 

 

 

 

 

 

 

“' - +7 l ~ i ;

 

 

 

i p

) J

 

‘”_!=0;

 

 

 

t

 

dwi

\

 

 

 

 

 

(

 

 

 

t

dWi- 2

\

 

 

(5)

 

U i+ —

 

-—

 

J F 2 i - \ — « 1 - 1+ ( U i —2 — — ■— ;

 

‘ ) ^ 11- 1 = 0 ;

 

 

 

 

 

2

ду

1

 

 

 

 

 

 

2

 

 

 

2

du

J

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

 

 

Г

 

 

 

c

(

 

dui-i

 

dWi-i

 

\

1

 

 

 

 

 

 

 

 

l “‘ 1- 7

\

 

 

dz

 

 

 

 

 

 

 

 

 

 

U>-2+

 

 

 

 

 

Г

 

 

c

 

(

dUi-з

 

 

dWi-з

\

1

 

 

 

 

 

 

 

+ l “i_3+7 l

ai

 

 

i p

)

J

f "“2"0;

 

 

 

 

 

 

 

 

 

 

dun+3

 

dwn+г

 

 

 

 

 

 

 

 

 

 

 

[- + 3 7c l V J z

 

 

i p )

\ l r _

 

 

 

 

 

 

 

Г

 

 

C

 

l

dUn+x

 

 

dwZ+x \

1

 

_

 

:

 

 

 

(+l “”+1+7 1~iz

 

 

 

ay

 

)J

 

‘"+2

 

 

l

_

t

dWn+2 \

 

_

W71 + I +

P

 

i

P

£

)

. F l n + x — 0;

I Un+---- ----- г------

I Г2П+ 1-

 

\

 

2

 

duy

J

 

 

 

 

 

\

 

 

 

2 dudy

'

 

 

W n F 2 n 1CO n -| + ffln2 F 1 n —1 = 0j

W n

—XF 2 n —2

W

n -

2 + W

n - 3 F x n —2

0,

W i F 2 i - [X — W i - l + W i - 2 F

i i - l =

Q-,

W

i - \ F

2 i - 2 —

W i - 2

+

W i - z F x i - 2 —

® \

(6)

COn+zF2n +2—®n+2"l"®n+|f ln+2—0; Wn +2^*2n + I COn +1 + W n F 1 n +1 0.

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

например, в работе [7].

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

16 порядка; она имеет общий вид

 

A n V i= P u

(7)

где i= 1,2, 3 , . . . , 6; j = N lt N2, N3, NT, W2T^з; индексы

1,2, 3,1,2, 3 относятся соответст­

венно к перемещениям наружных слоев; v n, ип, w n, v„, «й, иуй; Aij — некоторые диф­ ференциальные операторы, являющиеся функциями физико-механических характерис­ тик материалов и геометрических размеров слоев и оболочки, а также некоторых интег­ ралов сложных функций от z; Vj — компоненты тензоров перемещений двух наружных слоев; Pi — компоненты вектора внешней нагрузки, приложенной к наружным слоям. Необходимо отметить, что число слоев влияет на вид выражений интегралов, входящих в операторы Aij, но не влияет на порядок системы уравнений (7).

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

зоров

перемещений, деформаций,

напряжений

всех слоев

с помощью

выражений

(4) — (6), а также названных формул работы [7]. Особенностью полученной

системы

уравнений (7) является ее аналогия системе

уравнений

трехслойных

пластин и

оболочек.

 

 

 

 

 

В

качестве примера решена

задача изгиба

пятислойной

балки в двух

случаях.

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

лись

следующими: £ 2= £ 0= £Г=5,6 • 103 кгс/мм2; £ i = £ f = 3 ’ 102 кгс/мм2;

GXZi = G V2- =

= 8,5

кгс/мм2; V2= v o= v 2"=0,3; Vi = vT= 0,35; длина балки а = 120 мм;

^ = 0,2 мм;

с=0,05 мм; параметр нагрузки «70=0,02 кгс/мм2. В первом случае справедливым счита­

лось допущение об отсутствии сближения слоев,

т. е. юг = Wo=W2 :=wi = w1'.

На рис. 1 показано изменение нормальных

напряжений ох в наружном несущем

слое 2 и касательных напряжений i xz в слое связующего 1 по длине балки.

Рис.

1.

Изменение нормальных напряжений ах в слое наполнителя 2 и

касательных напряжений

xxz в слое связующего

I по длине балки нагруженной поперечной синусоидальной нагрузкой, шар­

 

 

 

нирно закрепленной за элемент среднего слоя 0.

 

Рис.

2.

Изменение нормальных напряжений ах в слое наполнителя 2 и

касательных напряжений

т хг

в слое связующего

1 по длине свободно опертой балки, нагруженной

сосредоточенной нагруз­

 

 

 

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

 

Согласно расчетам, максимальные касательные напряжения в слое связующего зна­ чительно ниже максимальных напряжений в слое наполнителя и, следовательно, рас­ слоение невозможно. На рис. 2 приведены расчеты тех же напряжений для случая на­ гружения балки сосредоточенной силой в середине пролета. Здесь в связи с ярко выра­ женной деформацией по нормали принимать допущение w2= w 0=w'2 = w l = wi нельзя.

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

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

1.Григолюк Э. И., Коган Ф. А. Современное состояние теории многослойных обо­ лочек. — Прикл. механика, 1972, т. 8, вып. 6, с. 3— 17.

2.Болотин В. В. К теории слоистых плит. — Изв. АН СССР. Механика и машино­ строение, 1963, № 3, с. 65—72.

3.Либреску Л. К уточненной линейной теории упругих анизотропных многослойных

оболочек. — Механика полимеров, 1975, № 6, с. 1038— 1050.

4.Новичков Ю. Н. Нелинейная теория и устойчивость толстых многослойных обо­ лочек. '— Прикл. математика и механика, 1973, т. 37, вып. 3„ с. 532—543.

5.Рябов А. Ф. Основные уравнения теории многослойных оболочек с учетом де­ формации поперечного сдвига. — В кн.: Сопротивление материалов и теория сооруже­ ний, 1968, № 3, с. 17—27 (Киев).

6. Прусаков А. П. Нелинейные уравнения изгиба пологих многослойных оболочек. Прикл. механика, 1971, т. 7, вып. 3, с. 3—6.

7. Власов В. 3. Общая моментная теория

оболочек. М., 1962 (Избранные

труды; Т. 1). 440 с.

 

Ростовское высшее военное командное училище

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

им. Главного маршала артиллерии Неделина М. И.

 

Московский физико-технический институт

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

 

1980, Ms 5, С. 929—933

УДК 624.074.4:539.3:678.067

В.Н. Бакулин, Д. А. Мысык

КРАСЧЕТУ ТРЕХСЛОЙНЫХ ОБОЛОЧЕК С ПЕРЕМЕННОЙ ТОЛЩИНОЙ

ЗАПОЛНИТЕЛЯ

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

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

Рассмотрим трехслойную цилиндрическую оболочку с несоосными несущими слоями разной толщины (рис. 1). Несоосность несущих слоев приводит к переменной толщине среднего слоя в окружном направлении. Оболочка нагружена внутренним р\ и внешним Рг давлением.

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

при таком нагружении и при таких характеристиках оболочки может привести к зна­ чительной погрешности, так как для несжимаемого в радиальном направлении запол­ нителя расчет будет проводиться на результирующую нагрузку, приложенную к одному из несущих слоев трехслойной оболочки. Но, как показано в работе [1], напряженно-де­ формированное состояние существенно зависит от того, к какому несущему слою при­ ложена нагрузка.

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

Несущие слои трехслойного прямоугольного конечного элемента [2], с помощью которого проводится расчет рассматриваемой трехслойной оболочки с переменной тол­ щиной заполнителя, подчиняются гипотезе Кирхгофа—Лява; в заполнителе учитыва­ ются деформации поперечного сдвига и нормальная деформация.

Расчетная схема приведена на рис. 1. Трехслойная оболочка оперта на краях каждым несущим слоем. Вследствие симметрии оболочки в окружном и меридиональном на­ правлениях решение проводится для четвертой части оболочки ABCD. Кинематические

граничные условия для узлов, расположенных

на линиях АВ и CD,

выглядят так:

u/= '6'i/ = 0; для узлов на торце

оболочки

(ВС)

ыу/ = 0; для

узлов на линии AD —

ut=ftVf=0.

 

 

 

 

 

 

 

Оболочка имеет

следующие

характеристики: —

= 12 ; —

= 5,7 • 10~3; ^ ша* =6,13;

 

 

 

 

R.3

Яз

 

^2

2/imin= ^2', — =4,33;

E\f=E2f\ Gi3=G23‘,

=6,67• 103; -—^-=2,8; ——=0,5; — =375

r2

 

^i3

 

G13

G13

p2

В данной работе

применяются обозначения,

использованные в [2].

По окружности

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

вом направлении

применяется равномерная

разбивка

на пять элементов.

 

 

 

 

На рис. 2 сплошными линиями показано

 

 

изменение радиальных

перемещений

несущих

 

 

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

 

 

нии

оболочки / х=

при разбиении

в коль­

 

 

цевом

направлении на семь элементов. Как

 

 

 

К)

 

 

 

 

видно из графиков рис. 2, радиальные пере­

 

 

мещения несущих слоев существенно разли­

 

 

чаются. Это объясняется тем, что заполни­

Рис. 2. Расчетные кривые прогибов. 1

тель сжимаем в радиальном направлении.

так как

внутреннее

давление превышает

*mln: ^

^ср> ^ *тах-

А