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

УДК 539.4:678.067

Ю. Н. Работное, В. П. Когаев, А. Н. Полипов, В. Б. Стрекалов

КРИТЕРИЙ МЕЖСЛОЙНОЙ ПРОЧНОСТИ УГЛЕПЛАСТИКОВ ПРИ ЦИКЛИЧЕСКИХ НАГРУЗКАХ

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

В настоящей работе критерий прочности композитов при статических нагрузках [2, 3], позволяющий объяснить зависимость межслойной сдви­ говой прочности от отношения длины пролета балки I к ее толщине А, развивается на случай действия циклических нагрузок.

1. Усталостные испытания образцов на межслойный сдвиг по методу трехточечного изгиба короткой балки проводили на машине, разработанной в Институте машиноведе­ ния им. А. А. Благонравова АН СССР. Образцы из однонаправленного углепластика на основе эпоксифенолформальдегидной смолы имели толщину 6 = 2 мм, ширину 6 = 10 мм, отношение //6 составляло 4, 5, 7,5, 9, 12,5. Испытания проводили на воздухе при комнат­ ной температуре с коэффициентом асимметрии цикла R= 0 и постоянной амплитудой перемещения нажимного ролика. Частота нагружения была равна 25 Гц. Сила Р, дейст­ вующая на образец, определялась по динамометру. Наибольшие касательные напряже­ ния при разрушении, именуемые обычно условной прочностью межслоевого сдвига, вы-

3Р числяли по формуле т==4^/1 ’ где нагРУзка Р измерялась в начальный момент нагружения

при заданной амплитуде смещения нажимного ролика. База испытаний была принята 106 циклов. За разрушение образца принимали момент падения его жесткости на 20—22%, причем в большинстве случаев межслойные трещины усталости визуально не обнаруживались. Пределы ограниченной выносливости т0 на базе 10б циклов при пульса-

ционном цикле (рис. 1) характеризуются величинами:

при

//6 = 4 то=22,8 МПа,

при

//6 = 5 — 22 МПа, при //6=7,5 — 20 МПа, при Г/6=9

18,5 МПа, при //6=12,5

15,5 МПа.

 

 

 

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

Рис. 1. Кривые усталости образцов из однонаправленного углепластика при трехточеч­ ном изгибе и //6 = 4 (У); 5 (2)\ 7,5 (5); 9 (4); 12,5 (5). Штриховыми линиями показаны средние значения пределов прочности при межслойном сдвиге ть, полученные по резуль­ татам статических испытаний. 0 — единая кривая усталости (5) для однонаправленного углепластика при растяжении.

касательные тху напряжения, что при некоторых допущениях [2, 3] мо#но

учесть в критерии прочности: ay+ mTjcy = Cf

(1)

где /п, с — экспериментально определяемые параметры материала; ось у направлена вдоль длины балки.

Согласно критерию (1) условная прочность при межслойном сдвиге

3Р

о

т = —-

оказывается не постоянной, как предсказывает элементарная оа~

лочная теория, а зависящей от отношения пролета балки к ее толщине

cm

(2)

т=- т2+ (Цh) 2 »

что подтверждается в экспериментах при квазистатическом нагружении {2, 3].

В случае действия циклических нагрузок постоянной амплитуды па­ раметры т и с можно представить в виде функций от числа циклов N. Для аналитического описания этих функций экспериментальные данные по циклической прочности были нанесены в координатах 1/т—(///i)2» в которых согласно (2) они должны ложиться на прямые линии 1/т = = ll[c{N)m(N)](l/h)2 + m{N)/c(N) при любом фиксированном числе цик­ лов. На рис. 2 все экспериментальные точки для статики N= \ и для фик­ сированных баз, взятых для удобства Af=104, 105, 106 циклов, получен­ ных интерполированием данных на рис. 1, хорошо ложатся на парал­ лельные прямые, что подтверждает справедливость зависимости (2) при циклических нагрузках. Значение т согласно [3] равно критическому значению l/h при смене механизма разрушения; при l/h<.m должно про­ исходить расслоение по критерию (1), при l/h> m — разрушение от нор­ мальных напряжений. Определенные по рис. 2 через коэффициенты пря­ мых линий значения c{N) и m{N) в зависимости от логарифма числа циклов представлены на рис. 3. Как видно, m{N) >13,5 и, следовательно, для исследованных отношений l / h ^ 12,5 должен осуществляться единый механизм разрушения, описываемый зависимостью (2). Согласие экспе­ риментальных точек с прямыми линиями на рис. 2 подтверждает этот вывод.

Зависимости параметров c(N) и m(N) на рис. 3 удобно описать ана­ литически с помощью двух простейших квадратичных функций, вклю­ чающих кроме статических значений параметров с0 = 502 МПа, /п0=13,5

всего лишь два экспериментальных безразмерных

параметра а\

и

с(N) = с0 (1 —ai lg2 N );

 

(3)

m(N)=m0(l + a2lg2N),

 

(4)

где а.\ = 5,60• 10-3; «2= 6,38’ 10-3. Зависимость (2)

с подстановкой

(3),

(4) принципиально не сложнее, чем обычная аппроксимация участков усталостных кривых прямыми линиями, но в отличие от прямых линий

Рис. 2. Зависимость обратной условной прочности при межслойном сдвиге 1/т от отно­ шения l/h при JV=10 (/); 10° (2)\ 105 (5); 10б циклов (4).

Рис. 3. Изменение параметров с (!) и т (2) в зависимости от \gN. Штриховая линия — прямая (6).

(3) и (4) обеспечивает единое описание зависимостей т {N) в пределах от статических нагрузок, по крайней мере, до 106 циклов. В то же время для определения а\ и а2 в дополнение к статическим данным при описа­ нии зависимостей т {N) при любых l/h необходимы всего две независи­ мые экспериментальные точки при произвольных l/h и N. В настоящей работе были обработаны все экспериментальные точки для подтвержде­ ния справедливости предложенного метода. Как видно из рис. 1, кривые, рассчитанные по формулам (2) —(4), удовлетворительно согласуются с результатами эксперимента.

Кривые усталости углепластиков при растяжении имеют небольшой угол наклона к оси абсцисс [4]. Согласно (2) с ростом l/h влияние на межслойное разрушение нормальных напряжений, характеризуемых па­ раметром с, растет по сравнению с влиянием касательных напряжений, характеризуемых параметром т. Из (3), (4) видно, что влияние числа циклов на параметр т несколько больше, чем на с, поэтому с ростом l/h, когда уменьшается вклад касательных напряжений, угол наклона ка­ сательных к кривым усталости в исследуемой области снижается и при­ ближается (см. рис. 1) к наклону кривой (0) при растяжении.

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

При обработке данных, приведенных в [4], оказалось, что для широ­ кого класса углепластиков (гладких и с отверстиями) участки усталост­ ных кривых в диапазоне 103—106 циклов очень хорошо аппроксимиру­ ются в безразмерном виде следующей линейной зависимостью в полуло­ гарифмических координатах:

атах/(Тв= 1,04-0,035 lg N.

(5)

Однако следует, отметить, что зависимость (5)

не обеспечивает точную

экстраполяцию на статическую прочность ав.

 

 

Были проведены также усталостные испытания на

растяжение образцов

из того

же однонаправленного углепластика, из которого вырезали балки для изгиба. Испыта­ ния образцов прямоугольного поперечного сечения 300X20X2 мм проводили при асим­

метричном растяжении (коэффициент асимметрии цикла R= 0,2) в режиме мягкого на­ гружения. Испытания проводили на электрогидравлической машине (фирмы MTS) с записью диаграммы деформирования. Частота нагружения в зависимости от уровня нап­ ряжений была выбрана в диапазоне 10—20 Гц. Образцы испытывали до полного разру­ шения. База испытаний была принята 106 циклов. Среднее значение статической проч­ ности ов = 760 МПа.

На рис. 1 приведены результаты испытаний на циклическое растяже­ ние образцов из однонаправленного углепластика, а также кривая 0, по­ строенная по зависимости (5). Как видно, и в этом случае единая кривая усталости хорошо согласуется с экспериментом.

Параметр с в критерии (1) по сути — изгибная прочность достаточно длинных образцов (///i>m ), когда можно пренебречь влиянием сдви­ гов. Хотя изгибная прочность в общем случае не совпадает с прочностью на растяжение, естественно предположить, что характер изменения изгибной прочности с ростом числа циклов в безразмерной форме такой же, что и для прочности на растяжение. Тогда для с считаем справедли­ вой зависимость (5), известную из независимых экспериментов:

с/с0= 1,04-0,035 \gN.

(6)

Прямая (6) изображена штриховой линией на рис. 3.

Учитывая, что для исследованного класса углепластиков данной тол­

щины наклон прямых на рис. 2 постоянен и

 

ШоСо

 

m(N)c(N)

 

независимо от числа циклов,

получаем возможность определить из (6)

зависимость m(N) и согласно

(2)

получить

 

т(Л0 = т 02 (1,04 -

с0т0

(7)

0,035 lg Л0- 2+ (///г)2

Сопоставление расчетных т(А)

Расчет по зависимости

(7) приведен

 

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

(7)

 

в таблице и хорошо согласуется со

с экспериментальными значениями,

средними

экспериментальными

ре­

полученными по средним долговечностям

зультатами

как

при

статическом,

 

для всех уровней испытаний

 

 

 

 

 

так и при циклическом нагружениях

Uh

 

lgN

 

Статика

(расхождение расчетных x{N) с экс­

 

 

периментальными

значениями

со­

 

 

 

 

 

 

 

 

ставляет менее 7%).

 

 

4

3,8855

4,9886

5,8721

В результате предложенного под­

 

28,4

26,1

24,4

35,4

хода по статическим испытаниям на

 

30

26

23

34,2

трехточечный изгиб (в принципе, до­

5

3,5900

4,9259

6,0000

статочно двух испытаний при раз­

 

27,9

25,4

23,4

33,8

ных l/h [3]) можно определить то и

 

28

25

22

32,4

с0, и, используя в качестве c{N)

не­

7,5

3,3464

4,4398

5,8250

зависимо определенную единую кри­

 

25,1

23,4

21,4

29,3

вую (5),

построить

усталостные

 

26

24

20

28,2

кривые в диапазоне по крайней мере

9

4,0934

4,3432

5,8238

до 10е циклов при любых ijh.

 

 

22

21,7

19,8

26,4

Можно осуществить и обратную

 

22

21

19

25,8

процедуру: проведя усталостные ис­

12,5

4,3032

5,1045

5,6977

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

 

17,5

16,9

16,4

20,4

какого-либо композита, по зависи­

 

18

17

16

20

мостям c(N)

типа представленных

Примечание.

В числителе

показаны

на рис. 3, оценить характер уста­

лостной

кривой

при

растяжении.

расчетные значения т(А ), МПа, в знаме­

Это приведет к экономии материала

нателе — экспериментальные.

 

и к снижению мощности необходи­ мых испытательных средств.

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

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

1. Rabotnov Yu. N., Polilov A. N. Strength criteria for In: Adv. Res. on the Strength and Fracture of Materials. June, 1977. Vol. 3B. Applications and non-metals. Sess. 3.

fibre-reinforced plastics. __ ICF 4, Waterloo, Canada, Composites, p. 1059— 1065*

N.Y., 1978.

2.Полилов A. H., Хохлов В. К. Критерий межслойной прочности композитов при поперечном изгибе. — Машиноведение, 1977, № 3, с. 56—59.

6.Полилов А. //., Хохлов В. К. Расчетный критерий прочности композитных балок при изгибе. — Машиноведение, 1979, N° 2, с. 53—57.

4. Гусенков А. П., Когаев В. П., Березин А. В., Стрекалов В. Б., Пешехонов Б. А., Звягин JI. /С. Сопротивление усталости углепластиков в связи с конструктивно-техноло­ гическими факторами. — Механика композит, материалов, 1981, N° 3, с. 437—442.

Институт машиноведения им. А. А. Благонравова

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

АН СССР, Москва

’ ’

УДК 539.4.001:678.067

А.Н. Гузь, М. А. Черевко

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

ПРИ СЖАТИИ

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

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

упругопластическими. Отметим, что часть результатов вошла в сообщо-

ние [10].

Основные соотношения. Рассмотрим двухкомпонентный волокнистый композитный материал. Полагаем, что содержание армирующих волокон небольшое, и взаимодействие между ними мало сказывается на значении нагрузки, при которой происходит потеря прочности композита, т. е. при­ мем модель матрицы с одним круговым волокном радиуса R. Волокно будем считать либо сплошным, либо полым с концентрической по­ лостью радиуса R\. Свяжем с осью волокна систему координат (г, <р, z) так, чтобы ось z совпадала с осью волокна. Нагружение вдоль волокна осуществляется путем одинакового деформирования как матрицы, так и волокна, задавая деформацию в направлении оси волокна e2z° = e. Для связующего используем соотношения теории малых упругопластических деформаций [11], полагая, что между интенсивностью напряжений и ин­ тенсивностью деформаций существует степенная зависимость

ои°=Аеп0к

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

Err0+ £фф°+ е= 0,

а для возмущений [9]

g-nmg-^jV

nUj = 0 .

(i)

В большинстве случаев достаточно учитывать влияние пластических деформаций только в матрице, так как в качестве волокон в основном используются достаточно хрупкие материалы, хотя в некоторых случаях необходимо учитывать и пластические свойства наполнителя. Для реше­ ния задачи воспользуемся линеаризованными уравнениями теории устой­ чивости при малых докритических деформациях, когда докритическое состояние определяется по геометрически линейной теории. Поскольку нагружение в большинстве случаев реализуется «мертвыми» нагрузками, то выполняются достаточные условия применимости статического метода исследования упругопластической устойчивости [9]. Поэтому как для матрицы, так и для волокна в докритическом состоянии должно выпол­ няться уравнение равновесия в силу осесимметричности нагружения

doгг® 1

— ^ — Ь— (<Тгг°—<Тфф°) = 0 ,

(2 )

а для возмущений уравнение

V ^ = 0 (*,/ = 1,2,3),

(3)

где для матрицы tij = Ki№ V $иа+ gingnjP, а для волокна /^ = со^аРУрИ*а. Предположим, что в композитном материале отсутствуют дефекты, т. е. считаем, что как матрица, так и волокно однородны, изотропны, а

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

иг°=и*°г;

и для возмущений

Ui= u*i;

I

•ч II

*»= <*»; r=R.

(4)

(5)

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

Рассмотрим два подхода к решению поставленной задачи: в первом случае полагаем, что докритическое состояние однородное, во втором —

неоднородное.

Однородное докритическое состояние. Д ля

волокна и матрицы при [9]

hi°= 'k20¥=h°\ п0=СТ22°;

8ц0= 6220

перемещения и деформации имеют вид

 

^n0 = 6in (Я*0 l)^ij 8nm®= 6nm (Ят°—1) .

В силу принятых предположений для матрицы потребуем выполнения условия несжимаемости

AI0+A.2°+A,3°= 0.

Возмущения перемещений в матрице представим в виде

1

д

д2

uv = ~

_<?d

i

d

2

uz= Д Х ;

Ur = — ~

drdz

W

r

— - X ;

г

дер

v

~дгdr

 

depdz

P = [ ( а ц - О 1з - р . ) А + ( ц + а з з 0)-^ -7 - ] ^ - - Х ,

а для волокна — в виде

 

J__c>

 

 

д 2

X *;

 

 

 

d

i

d

2

 

 

и*г = ~ ~ ' ¥ * ~

drdz

« % = - - —

--------д - т Х * ;

 

 

г

д ф

*

 

 

 

д г

 

г

depdz

 

 

 

 

а

ц

( A + £ ± ^ J L ) x .,

 

 

 

U*z—-

 

 

 

где

 

 

а*1з +

И-*

\

 

а*и

 

d z2 >

 

 

 

 

 

 

д2

1

д

1

 

д2

 

 

 

 

 

 

 

А = -

 

 

 

 

 

 

 

 

дг2

г

дг

г2

 

дер2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а ц =

£ — + —

(£ — 1) ] Ес; a33= — Ec ( 2 k + l );

 

 

a i3= —2a i2= — -^ -E c(k— 1);

р = - ^ Е с;

Ec= A ( l —X)h-';

 

 

 

a*u = a*33= V + 2 n * ;

о *12 =

а *1з =

Я*;

 

 

 

a33°=Ec(l - X);

а*«зз=(Я-1)-

 

и*(ЗЯ* + 2ц*)

 

 

 

V

+ ц*

-

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

Функции 'F, X, 'F* и X* определяются из уравнений

 

 

гд е

(4+s,si&-) ч,- 0;

( 4+ь!-^-)

 

( 4+ь!^ ) .х - о .

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{l. _ J S ± ^ ;

и , . , с ± у с* - * ± ? !

 

 

 

 

 

IX

 

 

 

г

 

 

 

 

 

 

 

 

2 с =

(Дп +

#зз +

СГ330— 2 # i3

2\i)/\i\

 

 

л

ц*+ а*°зз

 

 

ь

т/

с*?-

(#*зз+ а*°зз) {\1* + о*°гз)

 

£*»,=■

Г \з = с * ±

У

а-------- -------------------

 

 

LL

 

 

 

 

' J

 

 

 

 

# ПИ*

 

 

 

о*зз+ <т*°зз

[X* + Р*°зз

 

(а*1з+ц*)2 .

 

 

2с* =

*

 

 

а*п

 

а

п и

 

 

 

 

 

 

 

 

 

 

 

 

Я* и ц* — постоянные Ламе материала волокна.

 

 

 

 

Решение уравнений (6) представим для матрицы в виде

 

Т = Л 1Я , ( 1?1a )

sin

ф sin

a z ;

Х =

[Л 2Я , ( Ы

 

+ A 3H l (%3a ) ] c o s ф c o s <xz;

(7)

а Для воЛбкйа в форме

sin ф sinia2;

 

^ * = A J i

 

Х*=[Л5/, (tr2a) + A 6Ji (i£*3a)] cos ф cos аг.

(8)

Удовлетворяя условия контакта (5), получаем систему алгебраических уравнений относительно i4s:

6

 

% A 6bsh=0-, А=1,...,6,

(9)

S = 1

где коэффициенты bSk находятся из выражений для перемещений и соот­ ветствующих напряжений при подстановке в них (7) и (8).

Из условия существования нетривиальных решений (9) получаем характеристическое уравнение

| Щ = 0

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

Неоднородное докритическое состояние. Из условия несжимаемости находим напряженно-деформированное состояние матрицы в виде

Cx_^_ re

( Cl

e \

о C ‘

e

T

8rr°=

1

r2

+ 2 /

Byip-----p

T *

~ 2 ~ ’

l

 

II о

3 CO

l / 4

Cl2

E2-

2

, „

h

+

Grr0 = Po +

Л е г Л и 0''1- 1) •

r4

6

О

 

 

 

 

2

 

О^фф = Ро

О

Значение ро определяется из (2) с учетом условия «на бесконечности». Для волокна напряженно-деформированное состояние из (2) находим

в виде

и * С 2г+- у- \ а*°гг= (V + ц*) С22ц* ~ '+ Я *е.

Для сплошного волокна полагаем Сз=0. Неопределенные постоянные С\ и С2 находим из условия на поверхности контакта волокна и матрицы (4), а для полого С3 из условия на полости. В результате получаем

 

 

 

W

(k - Ъ ) (6 -5 ) ... (A -2S -1)

 

 

 

Cf_

 

X

Л ел- 3(6 -

•1) R 6 Н

+ Е

28sl(2s+3)

У

 

 

 

Х

 

- 2 ( ^ + I A* ) C 2 + V § ' - ^ 6 = 0 ;

 

Z\l

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

(3) для матрицы запишем в виде

ur=cpi (0 cos Ф cos az\ tt(p = 92(r)sin ф cos a z\

wz=cp3(r)cos ф sin a z\

l ч

(Ю)

p = p(r)cos ф cos az,

4

!гдё a = nkil -=* параметр волнообразования; / — длина полуволны во­ локна после потери устойчивости. Решение для волокна представим ана­ логично:

и*г= TJ?I (г) COS ф COS а г\

w%=^2(^)sin ф cos а z\ u*z=ty3(r)cos ф sin аг.

( И )

Подставляя (10) и (11) в уравнения (3) и учитывая для матрицы усло­ вие несжимаемости (1), получаем для определения ф* и ф* однородные системы обыкновенных дифференциальных уравнений

-^р-=А(г)у(г);

 

 

 

(12)

-^1=Л *(г)у*(г);

 

R > r ^ R h

(13)

. .

Г

 

dq>2

dcp3 \

 

у(0 = 1 ф1;ч»;фз; р;

dr

;

dr } ;

 

. . .

/ u

d^i

dty

d$з

\

У (0 -

\ Фь Фг, фз,

dr

dr

' dr

 

Системы дифференциальных уравнений (12) и (13) решали численно, воспользовавшись программой, приведенной в [12]. Переход от решений для матрицы к решениям для волокна осуществлялся с помощью усло­ вий контакта (5). На правом конце интервала интегрирования получаем при решении однородную алгебраическую систему уравнений. Из усло­ вия нетривиальности решений этой системы, приравнивая ее определи­ тель нулю, получаем относительно е характеристическое уравнение. Зна­ чение критической деформации еКр выбираем как наименьшее из всех полученных е.

Численная реализация и анализ результатов. Для получения число­ вых результатов для однородного и для неоднородного докритических состояний были составлены вычислительные программы. При численном интегрировании систем обыкновенных дифференциальных уравнений бесконечное тело заменялось цилиндром радиуса /?«>. Радиус волокна разделялся на N равных промежутков, а в матрице шаг интегрирования выбирался переменным — увеличивался с радиусом по формуле Ап =

R - R 1

. Для сплошного волокна не включался в область чис­

N

 

ленного интегрирования его центр, т. е. исключался цилиндр радиуса 0,02/?, причем увеличение его радиуса до 0,05R практически не отража­ лось на полученных результатах, как об этом свидетельствует сравнение,

приведенное в табл.

1 для /?оо =

 

 

 

 

Табл. 1

= 10R и N=3. Точность вычисле­

 

 

 

 

 

 

 

 

 

ний

контролировалась

сравне­

 

 

 

E/A

 

нием

результатов,

полученных

k

R iIR

2000

250

50

для различных /?«> и /V. В табл. 2

 

 

 

 

 

приведены значения екр для /?«> =

1

0,05

0,0160

0,0493

0,1161

= 10R и N = 3, а также для /?«,=

 

0,02

0,0158

0,0489

0,1160

= 30

и N = 6. Как видно,

числен­

0,75

0,05

0,0270

0,0704

0,1433

ное решение хорошо сходится, и

 

0,02

0.0260

0,0704

0,1431

поэтому

все представленные

ре­

0,5

0,05

0,0396

0,0900

0,1642

зультаты

были

вычислены

для

 

0,02

0,0395

0,0899

0,1634

Roo= 10R

и N = 3.

В дальнейшем

0,25

0,05

0,0527

0,1060

0,1741

результаты решения задачи с уче­

 

0,02

0,0525

0,1062

0,1748

том

неоднородности

докритиче-

0,1

0,05

0,0594

0,1111

 

ского

состояния

показаны

на

 

0,02

0,0591

0,1117

 

 

 

 

 

Табл. 2

 

 

 

екр

при

б,

 

k

EjA

 

 

 

R o o - m

я*, =зоя

%

 

 

 

N = 3

N=6

 

 

 

 

При R i / R = 0,2

 

 

1

4000

- 0,0112

-0,0116

3,5

 

 

250

-0,0510

-0,0521

2,1

 

 

20

-0,1835

-0,1853

1,0

 

0,1

4000

-0,0493

-0,0470

4,9

 

 

 

При R , I R =0,75

4,7

 

1

4000

-0,0203

0,0213

 

 

250

-0,0873

-0,0887

1,6

 

 

20

 

0,1

4000

-0,0690

-0,0670

3,0

Рис. 1.

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

На рисунках показаны графики зависимости енр от отношения EjA

для сплошного волокна

(рис. 1) и для волокна с радиусом полости R i=

= 0,7R (рис. 2), когда

коэффициент Пуассона волокна v = 0,3; от v во­

локна для случаев однородного и неоднородного докритического состоя­

ния при £/Л = 250

(рис. 3); от радиуса полости,для 6 = 0,25 при различ­

ных EjA (рис. 4);

а также зависимости екр

от отношения показателей

скоростей деформационного упрочнения AJA

волокна и матрицы, если

волокно упругопластическое для 6 = 0,1 и различных значений коэффи­ циента деформационного упрочнения волокна 6В (рис. 5). При 6=1 при­ ведены результаты для расчетной схемы, когда волокно и матрица принимались упругими, что соответствует решению [13]. В работе [14] экс­ периментально определялось значение критической деформации однонап­ равленного волокнистого композита чистый алюминий — нержавеющая

сталь, и для малой концентрации наполнителя

(4,1%) найдено, что

екр~ 6%. Приняв для чистого алюминия Л = 108 Па

и 6 = 0,1 [14,

15] и

считая волокно упругим с модулем Юнга Е = 2 - 1011

Па, по рис. 1

полу­

чаем значение е1ф ~ 5,9%. Если принять волокно

упругопластическим,

то согласно [16] для нержавеющей стали АМ-355, содержащей те же ком­ поненты, что и сталь волокон в работе [14], приближенно имеем Ав= = 2,0-109 Па, 6В= 0,12, и тогда екр«7% (см. рис. 5).

Полученные

результаты

показы­

-кр

 

 

/V

вают,

что

увеличение полости в

во­

 

 

 

 

 

 

0 5 / /0'35

локне

или

учет

упругопластических

0.15 I-

 

 

 

 

 

свойств волокон увеличивают значения

 

 

 

 

критического укорочения е«р, при этом

0,125 k

 

 

 

кривые для различных коэффициентов

 

/ / / / / / / Л К

деформационного

укорочения

имеют

VK

/ / / / ' / / / / V1

качественно

одинаковый характер

из­

/ / / /

/ / / / Л/кв=°п

менения (см. рис. 1, 2, 5) в зависи­

У'

/ / /' / / /

мости от отношения жесткостных пара­

0,075b

/ / / /

У

метров волокна и матрицы. Сравнение

"' /'

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

0,05

 

 

 

ношений EJA от 4000 до 20 результатов

 

 

 

 

с учетом и без учета неоднородности

0.025

 

 

 

докритического

состояния

показало,

 

 

 

Ав/А

что эти два подхода дают близкие ре­

 

 

 

 

 

 

 

зультаты

(различающиеся

менее чем

500 250 125

62,5

31

16

на 5%) для всего диапазона изменения

 

Рис.

5.

 

коэффициента Пуассона

волокна

и

 

 

лишь для Е/А<^20 возможно и необ­ ходимо учитывать неоднородность докритического состояния. Это позво­

ляет сделать заключение о возможности проводить исследование устой­ чивости арматуры в упругопластической матрице, не учитывая неодно­ родности докритического состояния (см. рис. 3). Исследование зависи­ мости критической деформации от радиуса полости в упругом волокне показало, что для малых радиусов полости, когда R \ /R < .0 ,2-^-0,3, ре­ зультаты менее чем на 5% отличаются от результатов для сплошного во­ локна. Поэтому, возможно, используемые для изготовления, например, борных волокон нити вольфрама, поскольку £ ВольфРама <Ебора, по-види­ мому оказывают очень малое влияние на прочность композита. Отметим также, что чем больше полость в волокне, тем меньше влияние неодно­

родности.

Выводы. 1. Привлечение трехмерной линеаризованной теории устой­ чивости для решения задач о разрушении однонаправленных^ волокнис­ тых композитных материалов с упругопластической матрицей при сжа­ тии является эффективным и позволяет достаточно точно предсказать разрушающие нагрузки, о чем свидетельствует сравнение с эксперимен­ тальными результатами работы [14], для которых eKp^6% , тогда как в настоящей работе для упругих волокон еКр~5,9%, а для упругопласти­

ческих 7% 2. Для отношений Е/А в пределах от 20 до 4000, как показали резуль­

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

получить, только предполагая однородность докритического состс^а

3. Полость в волокнах, имеющая радиус меньше 0,2 0,3 радиус^

 

локна, оказывает небольшое влияние на критическое укорочени

на_

можно пренебречь. Для больших полостей необходим их учет, та

увеличением радиуса полости происходит уменьшение критичес

 

грузки.

 

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

материалы3 М ^1% 7*сН5 4 -9 6 Р°ЧНеНИЯ К0МП03ИЦ|1Й- — В км.: Волокнистые композитные

зитя n ± T Z

А" К°3ак М- Дж- Влияние поверхности раздела на характеристики компо-

К мпозн хР У К

°1 Л паС Т И Ч е С К 0 Й 0бласти- -

В

™ : Поверх' ости раздела в металлических

композн^х Композиционные материалы. М„ 1978, т. 1, с. 231-265 .

жур,1. прим.

i

 

 

-

m.i«x c„s,P

„"а,«,р1й |с1 % ^

г к

» 1‘Ж

Р:п, й й 1 ьогоп ,iЬгe'mda,

63 — 1316

5.Гузь О. М. Про виЗиаЧення Теоретично? границ! мщиост! при стиску армо&ЗиИК матер1ал!в. — Доп. АН УРСР. Сер. А, 1969, № 3, с. 236—238.

6.Гузь А. И. О построении теории устойчивости однонаправленных волокнистых

материалов. — Прикл. механика, 1969, т. 5, № 2, с. 62—70.

7. Бабич И. Ю., Гузь А. Н., Щульга И. А. Исследование динамики и устойчивости композитных материалов в трехмерной постановке. — Прикл. механика, 1982, т. 18,

1, с. 3—27.

8.Гузь А. Н. Устойчивость трехмерных деформируемых тел. Киев, 1971. 276 с.

9.Гузь А. И. Основы теории устойчивости горных выработок. Киев, 1977. 203 с.

10.Гузь А. Н., Черевко М. А. К механике разрушения волокнистого композитного материала с упругопластической матрицей при сжатии. — Докл. АН СССР, 1981, т. 256,

4, с. 806—809.

11.Ильюшин А. А. Пластичность. Ч. 1. Упругопластические деформации. М.; Л., 1948. 376 с.

12.Шинкарь А. И., Китайгородский А. Б., Борщевская С. К. Решение линейных диф­

ференциальных уравнений (программа на языке ФОРТРАН). — В кн.: Алгоритмы и программы решения задач механики твердого деформируемого тела. Киев, 1976,

с.157— 170.

13.Бабич И. Ю. Об устойчивости волокна в матрице при малых деформациях. — Прикл. механика, 1973, т. 9, № 4, с. 29—35.

14.Pinnel М. R., Lawley A. Correlation of unianial yielding and substructure in

aluminium-stainless steel composite. — Met. Trans., 1970, vol. 1, N 5, p. 1337— 1348.

15.Промышленные деформируемые, спеченные и литейные алюминиевые сплавы. Справочное пособие. М., 1972. 552 с.

16.Бабаков А. А., Приданиев М. В. Коррозионностойкие стали и сплавы. М., 1971. 320 с.

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

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

Киев

 

УДК 539.375.001:678.067

М.Б. Победря

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

Вшироком классе композитных материалов особое место занимают волокнистые композиты. Актуальной проблемой является прочностной расчет деталей и конструкций из них. В настоящее время существуют различные подходы к решению данной проблемы [1—4]. Во многих из этих подходов, в том числе и статистическом, требующем определения

«неэффективной» длины волокна в композите [5], необходимо знание пе­ рераспределения напряжений в композитном материале при разрушении одного волокна. При этом информацию о перераспределении чаще всего надо иметь в «аналитическом» виде. При некоторых указанных ниже предположениях перераспределение напряжений в композите при разру­ шившемся упругом волокне с учетом упругих и упругопластических свойств матрицы в статической постановке было найдено в работе [6]. Динамика перераспределения напряжений при разрушившемся упругом волокне в упругой, матрице в предположении, что соседние (неразрушен­ ные) волокна закреплены, была исследована в работе [7]. Решение такой же задачи при отказе от гипотезы о закреплении соседних волокон было дано в [8].

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

1. Постановка задачи. Рассмотрим волокнистый композит, содержа­ щий разрушенное волокно. Все величины, относящиеся к волокну, будем снабжать нижним индексом /, а величины, относящиеся к матрице, ин­ дексом m [4]. Материал нагружается в направлении армирования (по оси z). После того, как однородная деформация некоторого волокна до­ стигла значения е/, произошло его разрушение. Обозначим перемещение разрушенного волокна вдоль оси z через и0, а соответствующую дефор­ мацию через 80 и вообще присвоим разрушенному волокну лндекс 0, во­ локнам, непосредственно окружающим его, индекс 1, волокнам следую­ щих поясов индексы 2, 3, 4 и т. д. Например, перемещения волокон соот­ ветствующих поясов вдоль оси z *отсчитываемые от плоскости разрыва, будут обозначаться ии и2, и3 и т. д. Тогда, очевидно, что

6 i= e /+ ^ T ; i=Q' *’ ■

С1-1).

Будем считать справедливой следующую гипотезу [4, 6—8]. Сдвиговые деформации на границе некоторого волокна прямо пропорциональны перемещению этого волокна относительно соседних и обратно пропор­ циональны расстоянию от его поверхности до поверхности соседних во­ локон. Для гексагональной упаковки волокон расстояние между волок­ нами является функцией угла 0, периодически изменяющегося с полупериодом л/6. Обозначим это расстояние через 6(0) [6].

Распределение Сдвиговых дефорЫСцПй иб йерйметру разрушенного йолокна и волокон, окружающих разрушенное, будет иметь вид

Уа=

'» i,/ = 0 ,1 ,2 ,...; 1Ф'].

(1-2)

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

тij(t)=Gm [ у .;( 0 - 1

0(£-т)у»Ит)^т] ,

(1.3)

о

 

 

где Gm — мгновенный модуль сдвига матрицы; G {t) — ядро сдвиговой релаксации матрицы. Точно так же для связи между растягивающими

напряжениями в волокнах a,, t'=0,1,2,... и соответствующими деформа­

циями (1.1) можно записать

 

t

 

Oi{t) =Ef [ei(/) - J E(t-x)zi{x)dx ]

(1.4 )

о

 

где Ef — мгновенный модуль Юнга волокна, £ (/) — ядро релаксации волокна при растяжении.

Рассматривая силы, действующие на элемент волокна, получим урав­

нения движения

 

 

 

 

да

Ti

_

д2щ

(1.5)

dz

nrf2

^

dt2

 

где rf — радиус волокна; Ti — касательные силы, действующие по еди­ ничной длине волокна t-ro пояса; р/ — плотность волокна. Касательная сила, действующая на элемент разрушенного волокна, равна

я/6

Т0= J

XoirjdQ= 12 J Xoir/dQ.

оо

Суммарная касательная сила, действующая на единичной длине волокна 1-го пояса, равна

Я /6

 

Я /6

 

Я /6

 

 

Т\ = 2

J t i o / * / d 0 +

4 J* T i 2 r / d 0

+XxtffdQ,2 J*

0

 

0

 

0

 

 

на единичной длине волокна 2-го пояса —

 

 

 

я/6

 

я/6

я/6

 

Г2 = 4 J t2iT/d0 + 4 J

Х23r/d.0 + 4 J

T24ffdQ,

3-го пояса —

о

 

о

%

0

 

 

 

 

 

 

 

 

 

Я/6

 

 

я/6

 

я/6

 

я/6

T3=2 J x3irfdQ + i

j T32r/d0 + 4

J Тз4гД0 + 2

J x3srjdQ

0

0

 

 

0

0

 

 

и т. д. для Ti, Tb, ... —

Р/

я/6

 

 

 

 

 

 

 

 

 

Ti=

^ 1

J

rfAihXihdQ,

i=4, 5 ,...

(16)

 

k = a i

0

 

 

 

 

 

где

a ,^ 2 ,

(a < p ) — некоторые целые числа { k ^ i) .

Подставляя

(1.6) в уравнения (1.5) и учитывая физические соотно­

шения

(1.3), (1.4), а также геометрические соотношения (1.1) и (12),

получим, вообще говоря, бесконечную систему интегродифференциальных уравнений относительно перемещений щ (i= 0 ,1 ,2 ,...):

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

Л-а.

 

 

t

 

 

 

 

 

 

 

X [ Ui (z, / ) — I G (/—T ) «

, ( Z , т ) е / т

] + - ^ - ( X ^ i *

) [uk(z,t)~

 

 

О

 

 

 

 

 

 

 

 

 

t

 

 

 

1

d2^ (z, t)

 

 

 

- i

G (t- x)u h(z,x)dx ]

 

 

(1.7)

 

 

~a2

dt2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

где а2=

 

2Gm

6

f

cB

 

К системе

уравнений (1.7)

 

Efrjbc

я

J0

ft(0)'

 

 

следует добавить граничные условия

 

 

 

 

 

 

д«о(0, t)

= -er,

Ui(0,t)=0

({=1,2,.

.)■,

Ui(z,i)

О при z-+oo

dz

 

 

 

 

 

 

 

 

( 1.8)

 

 

 

(i = 0, 1,...)

 

 

 

и начальные данные

 

 

 

 

 

 

 

 

Ui(z, 0)=0; dUif ’-0)~ =Q

(i = 0 ,l , ... ).

 

(1.9)

Таким образом, задача заключается в решении бесконечной системы уравнений (1.7) при выполнении граничных условий (1.8) и начальных данных (1.9). Если такая задача будет каким-либо образом решена, т. е. найдены перемещения щ(г, t) всех поясов, то по формулам (1.1) и (1.2) могут быть найдены соответственно деформации растяжения волокон и деформации сдвига на границе волокна и матрицы, а по формулам (1.3)

и(1.4) можно будет определить растягивающие напряжения в волокнах

икасательные напряжения в матрице.

2.Частные случаи. Поставленная в предыдущем пункте задача

(1.7) — (1.9) является динамической задачей теории вязкоупругости. Если в уравнениях (1.7) пренебречь силами инерции, то получим урав­ нения равновесия

Ьг(и)= 0; £= 0, 1 ,2 ,...

(2.1)

Тогда задача, заключающаяся в решении уравнений (2.1) при удовлет­ ворении граничным условиям (1.8), является квазистатической задачей теории вязкоупругости. Если волокна являются упругими, в уравнениях (1.7) или (2.1) следует положить

£ ( 0 = 0 .

(2.2)

Если матрица также окажется упругой, т. е.

G ( t ) ^ 0,

(2.3)

то задача (1.7) —(1.9) будет динамической задачей теории упругости, а задача (1.8), (2.1) — статической задачей теории упругости.

Чтобы оборвать бесконечную цепочку уравнений (1.7) или (2.1), не­ обходимо ввести дополнительную гипотезу. Назовем Л/'-волоконной мо­

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

 

 

«л-= 0;

MJY+I = 0;

«JV+2 = 0 ...

 

(2-4)

Например,

для

одноволоконной модели (т. е. для

модели

при «1= 0,

ы2 = 0,...)

имеем из (1.7) одно уравнение движения

 

 

т ,

ч

d2u0(z,t)

t

Btt ^ d 2u0(z,x)

J

^

f

LQ(U) —r

 

J

T).

Qz2 dX

V [ W(2. 0

J G(t-x)up {z>T)dT] = 4

d2«o(z, t)

(2.5)

dt2

с граничными условиями

 

 

 

 

du0(0, t)

= -e/;

lim «0(z, 0 = 0

(2.6)

dz

 

Z-^OO

' '■

 

 

 

 

 

и начальными данными

 

 

 

 

«o(z, 0) =0;

dup{z, 0)

=Q

(2.7)

 

 

dt

 

 

Соответствующая квазистатическая задача будет заключаться в реше­ нии уравнения равновесия

 

 

 

 

 

Lp{u)= 0

 

 

 

 

 

(2.8)

при выполнении граничных условий (2.6).

 

 

 

 

(2.2) и

(2.3)

Решение задачи

(2.5) —(2.7)

при выполнении условий

дано в [7], а задачи

(2.8), (2.6)

при тех же условиях

(2.2) и (2.3)

— в

работе [4].

 

 

 

 

 

 

 

 

 

 

 

Для четырехволоконной модели, т. е. при условии

 

 

 

 

 

 

«4= 0;

«5= 0;

«6=0

 

 

 

 

 

(2.9)

получим из (1.7) систему уравнений

 

 

 

 

 

 

 

Mp{u) S

-

J E ( t - x ) ~ f

z; T)

 

d x - y 2{

Up(z, t) -« i(z , 0 -

 

t

 

 

 

 

 

)

1

d2Uo(z, t)

 

 

f

 

 

 

 

 

 

 

~ J G(t-x)[u0{z,x) - u l{z,x)]dxf = -^------^ ---- ;

 

M■(“) =

d2Ugz2’ t] ~

j

E (t- x)

 

dx

 

{ 4«!(z, t)-Up(z, 0 -

—2«2 (z, t) —«3(%, 0 — J G(t-x)[4ul(z,x)-Up(z,x)-2u2(z ,x ) -

 

 

 

 

o

 

 

 

 

 

 

 

 

 

 

\

1

d2ui(z,t)

;

 

( - ) ..■

d2u2 (z, 0

 

 

 

 

 

 

dp

 

 

a - — ~

 

 

t

 

 

 

 

 

 

 

 

 

 

 

-

J E (t - x ) - 2ц0

’т)

d x -

{б«2(z, t) —2«i (Z, 0 - 2 a 3(z, 0 -

 

t

 

 

 

 

 

 

1

 

1

d2Uo(z,t)

 

f

 

 

 

 

 

 

 

- J

G (/—t)[6«2(z,T)-2Mi(z,T)-2«3(z,T)]e?T J =

a2

dt2

 

M / ,л _

d2u3(z,t)

t

. v2 f

x d°-u3(z,x)

Мз{и\—

} E ( t - x ) -----—

d x - —- {бы3(2, 0 - и ,( 2> <>-

- 2 « 2(Z, 0 - J 0(/-т)[6и8(г,т)-и1(г1т)-2и2(г1т ) ] * } = —

о

 

 

 

 

i

a2 df2

 

 

 

 

 

 

(2.10)

Заметим, что уравнения (2.10) получаются из (1.7), если

в последних

положить

 

 

 

 

 

 

i= 0, 1,2,3; о0 = Ро = 1; Л0х= 12;

 

«1= 0, Pi = 3; Лю=2;

^4I2 = 4;

А\3=2',

 

« 2 = 1 ,

Рг= 4 ;

Л 21 = 4 ;

Л 2з = 4 ;

Л 24= 4 ;

 

«з=1;

Рз= 5;

Л31= 2;

Лз2 =4;

Лз4 = 4; Лз5 = 2

 

и учесть соотношения

(2.9). К уравнениям

(2.10) следует добавить гра­

ничные условия (1.8)

и начальные данные

(1.9), в которых также сле­

дует учесть условия

(2.9).

 

 

 

 

Квазистатическая задача для четырехволоконной модели заключа­

ется в решении четырех уравнений

 

Mi (и) =0; i= 0, 1,2,3

(2.11)

при удовлетворении граничным условиям (1.8) с учетом (2.9). Стати­ ческая задача (2.11), (1.8), (2.9) при выполнении условий (2.2) и (2.3) решена в работе (4], а динамическая задача для трехволоконной модели при выполнении условий (2.2) и (2.3) — в работе [8].

Заметим, что для обрыва бесконечной цепочки уравнений (1.7) или (2.1) могут быть приняты и другие условия, отличные от условий (2.4). Например, в отличие от «жесткой» четырехволоконной модели (2.9) мо­ жет быть рассмотрена «податливая» четырехволоконная модель, для ко­ торой вместо условий (2.9) можно записать [4]

T i3 = 0 ; т 2з = 0 ; т 24 = 0 ,

В этом случае получим систему уравнений, отличную от (2.10) или (2.11). Для случая упругих волокон, т. е. при выполнении условий (2.2)

и(2.3), для статической задачи в [4] показано, что решения для жесткой

иподатливой моделей различаются незначительно.

3.Решение поставленной задачи. Применим к системе уравнений

(1.7), записанной для

^-волоконной

модели, преобразование Лап­

ласа [10]:

оо

 

 

 

u*i(z,s) =

J e~stUi(z,t)dt

(i= 0 ,1 ,2 ,...);

 

о

 

Е*(s) = J

e~stE{t)dt\ G*(s)z= J er«G(t)dt,

о

 

о

где s — комплексная переменная. Тогда с учетом начальных данных (1.9) получим систему обыкновенных дифференциальных уравнений для

изображений и*i

d2u*j{z,_s)

+ V B i}(s)U*)(z,s)=0 \

(3.1)

dz2

. j=0

 

Здесь

£ <«-e-> ( iU) ] * . +

 

A . W - i - V t

 

 

"=<Xf

(3.2)'

где величины Aijf входящие в соотношения (1.6) и (1.7), обладают сле­ дующим свойством: Aij=0 при i=j и j<.au /> Р ь а бij — символы Кро-

некера [111-

уравнений (3.1) при удовлетворении граничным

Решение системы

условиям

(1.8) имеют вид

 

 

 

 

 

N

 

 

 

 

u*i(z, s)= ^

Cii(s)e-Pj<s>z;

i= 0, 1... ,N,

(3.3)

где Pj(s);

 

j-o

 

 

 

/=0,1

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

корни характеристического

уравнения

(предположим для простоты, что все они различны)

 

 

det

| Bn+p28ij \ =0; t,/ = 0 ,1 ,.... TV,

(3.4)

а величины Сц (б) определяются из уравнений

 

 

N

 

N

 

 

 

PiCoj (s) = —

’, ^ jC ,j(s)= 0 ; t = 1, 2 , ... , Л/,

 

 

j=o

s

j=l

 

 

которые являются следствием граничных условий (1.8).

(3.2) счи­

Для решения квазистатической задачи нужно в выражении

тать, что а->-оо. Если в композите волокна являются упругими, следует

в (3.2) положить £* = 0, а если

матрица является упругой, то G*s=0.

Чтобы по заданному изображению

(3.3) найти оригинал, можно вос­

пользоваться обратным преобразованием Лапласа [10]

 

1

а + г о о

Ui(z, t) =

|

u*i(z, s)eatds.

 

2ni, — / о о

 

В ряде случаев эффективным методом приближенного обращения преоб­ разования Лапласа может оказаться метод Шепери [9], с помощью ко­ торого оригинал записывается в виде

j=o

Для одноволоконной модели из (3.3) имеем

 

u*0(z,s)

£-Р0(S)Z^

(3.5)

где

spo(s)

 

с2

1/2

 

1-G*(s) + -V

OiL

po(s) =Y

1 -E*(s)

так что приближенное решение по методу Шепери имеет вид

«о (г, 0 =

Z

 

Для четырехволоконной модели

характеристическое уравнение (3.4)

может быть решено только приближенно:

, S 2

1/2

- ^ + M 2[l-G * (s)]

Pi(s)

i= 0,1,2,3,

где b0= 1,12848; 6, = 4,73509; b2 = 7,26491; b3 = 7,62155.

В заключение заметим, что для описания поведения оригинала при t-*~оо можно воспользоваться предельной теоремой [10]

Пт щ (t) = Пт su*i (s).

t—►oo s-*0

Например, для одноволоконной модели из (3.5) имеем

Пт «o(z, t) =

 

е/ g - P o ( 0 ) Z j

t—►оо

Ро(0)

где

 

1/2

 

 

1 -

J

G(t)dt

Ро(0)=у

О_______

оо

 

1 -

j

E(t)dt

 

о

 

Автор благодарен В. В. Москвитину, под руководством которого вы­ полнена данная работа.

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

1.Малмейстер А. К., Тамуж В. П., Тетере Г. А. Сопротивление полимерных и, композитных материалов. 3-е изд. Рига, 1980. 572 с.

2. Волокнистые и дисперсно-упрочненные композиционные материалы. М.,

1976. 213 с.

3.Механика композиционных материалов. М., 1978. 556 с.

4.Копьев И. М., Овчинский А. С. Разрушение металлов, армированных волокнами.

М., 1977. 240 с.

5.Работное Ю. Н. Механика деформируемого твердого тела. М., 1979. 744 с.

6.Овчинский А. С., Копьев И. М., Сахарова Е. Н., Москвитин В. В. Перераспреде­

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

териалах. — Механика полимеров, 1977, № 1, с. 19—29.

7. Москвитин В. В., Овчинский А. С. Динамика перераспределения напряжений в разрушившемся волокне при упругом деформировании компонентов композиционного материала. — Изв. АН СССР. Механика твердого тела, 1979, № 1, с. 120— 124.

8. Москвитин В. В. Динамика перераспределения напряжений в волокнах компози­ ционного материала при упругом деформировании компонентов. — Вести. Москов. ун-та.

Сер. 1. Математика, механика, 1980, № 6, с. 95—98.

 

9.

Москвитин В. В. Сопротивление вязкоупругих материалов. М., 1972. 328 с.

10.

Дёч Г. Руководство к практическому применению

преобразования Лапласа и

Z-преобразования. М., 1971. 288 с.

 

11.

Победря Б. Е. Лекции по тензорному анализу. М.,

1979. 224 с.

Московский государственный университет

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

им. М. В. Ломоносова

 

УДК 620.178:678.067

М. Я. Микельсон, А. X. Курземниекс

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

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

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

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

1231-У 10 при посто­

янной амплитуде напряжения с асимметрией цикла R= 0,02 при

частоте 2 цикла в

минуту и базе испытаний N = \03 циклов. Для измерения силы использовали динамометр, входящий в комплект испытательной машины 1231-У 10; для измерения перемещений ис­ пользовали измеритель деформации, описанный в [2]. Во время усталостных испытаний сигнал, поступающий с динамометра и измерителя деформации, регистрировался и на­ капливался в памяти ЭВМ с частотой отсчета два раза в секунду.

После определенного числа циклов нагружения образец вместе с захватами снимали с испытательной машины 1231-У 10 и вместе со специальной приставкой устанавливали на рентгеновский аппарат ДРОН-2. Образцы просвечивались излучением трубки с мед­ ным анодом, и путем шагового сканирования регистрировались меридиональные реф­ лексы, обусловленные структурой армирующих волокон материала [3]. По изменениям характеристик рефлексов определяли соответствующие изменения в структуре волокон, которые прямым образом связаны с механическими характеристиками материала в це­ лом. Для испытаний использовали однонаправленный органопластик с коэффициентом армирования ф =0,6. Размеры рабочей части образцов 1X 4X60 мм.

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

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

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

напряжения

 

 

 

 

Касательные

модули

упругости

 

 

 

 

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

в процессе

циклического нагружения

минимальный

 

 

перед

Сечение

исходный,

 

 

£тах * К)'3 МПа

 

 

разрушением

при о,

 

 

 

 

 

 

Е • 10_3 МПа

МПа

Е • ЮЛ

 

 

 

 

 

 

 

 

МПа

 

 

При О /V/ ., МПа

 

 

 

 

700

1050

1370

1470

1500

 

1570

300

62,1

63,8

67,4

67,7

71,1

66,6

63,5

57,3

500

59,4

66,8

70,7

72,3

73,7

70.0

65,1

59,4

700

•49,0

72,9

74,5

76,7

74,0

67,3

61,6

900

39,5

74,3

76,7

80,6

77,2

70,6

64,8

1100

42,3

78,1

83,0

79,3

73,2

67,6

1300

48,2

79,0

74,8

65,1

изменения механических характеристик однонаправленного органоплас­ тика во время малоциклового нагружения. Средняя прочность мате­ риала при статическом растяжении составляет 1670±50 МПа. На рис. 1 показана зависимость прочности от количества циклов нагружения. Из рисунка видно, что прочность органопластика слабо зависит от коли­ чества циклов нагружения.

Для выяснения влияния малоциклового нагружения на остаточную прочность однонаправленного органопластика образцы подвергались предварительному нагружению 1000 циклов в диапазоне напряжений от 700 до 1500 МПа и после 20-ч отдыха испытывались на статическое рас­ тяжение до разрушения. Средняя остаточная прочность предварительно нагруженных образцов составила 1690±60 МПа. Разница между исход­ ной и остаточной прочностями статистически незначима.

Изменение механических свойств однонаправленного органопластика можно проиллюстрировать изменением зависимости а —е одного образца в разных циклах нагружения. На рис. 2 в качестве примера приведена зависимость о —г в 1-м, 100-м циклах нагружения и перед разрушением при о= 1490 МПа. Из рисунка видно, что кривая о —г в 100-м цикле идет более круто по сравнению с соответствующей кривой в 1-м цикле нагру-

б МПа

Рис. 1. Зависимость разрушающего напряжения от количества циклов нагружения.

Рис. 2. Зависимость а —е в разных циклах нагружения: 1 — 1-й цикл; 2 — 100-й, 3 — перед разрушением.

жения; уменьшилась также площадь

Суммарная удельная рассеянная энергия

петли гистерезиса. На рис. 3 в ка­

при разных величинах циклического

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

 

напряжения

 

касательного

модуля

упругости в

Максимальное

Количество

Суммарная

зависимости

от количества

циклов

нагружения в сечении кривой а —е

циклическое

циклов

удельная

напряжение

нагружения

рассеянная

при а = 300 и 900 МПа для разных

ст ^ МПа

до разрушения

энергия W,

 

N, циклы

МПа

величин

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

напряжения.

 

 

30,0

Как видно из рисунка, во время

•1670

1

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

нагружения при базе

1570

13

35,5

испытаний N= 103 циклов и напря­

1560

18

44,9

1500

420

328,5

жениях а„=<1470 МПа значения ка­

1370

1000

620,0

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

 

(не разру­

 

чиваются. При циклических

напря­

 

шился)

 

жениях

с~;>1470

МПа,

когда

происходит усталостное разрушение образцов на базе N=103 циклов, касательные модули упругости в на­ чальной стадии нагружения увеличиваются, но не достигают таких зна­ чений, как при а^ = 1470 МПа. После определенного числа циклов на­ гружения модуль упругости начинает снижаться вплоть до разрушения образца. Для иллюстрации в табл. 1 приводятся численные значения модулей упругости и их изменение в результате циклического нагру­ жения. На рис. 4 кривые 1, 2 показывают изменение площади петли гис­ терезиса в зависимости от количества циклов нагружения. При напря­ жении г^1470 МПа площадь петли гистерезиса в процессе нагруже­ ния уменьшается. При а~>1470 МПа после определенного числа цик­ лов нагружения площадь петли гистерезиса начинает увеличиваться. В табл. 2 приведены численные значения суммы рассеянной энергии до разрушения образца. Как следует из таблицы, при уменьшении вели­ чины циклического напряжения суммарная рассеянная энергия увели­ чивается.

На рис. 4 приведена также зависимость изменения минимальной и максимальной деформации от количества циклов нагружения при раз­ ных уровнях циклического напряжения (кривые 3,4). Минимальная де­ формация перед разрушением может достигнуть величины 1,5%. Из рис. 4 видно, что при N<100 циклов emin выходит на линейный участок.

R ГРЧРП1Ш кпипп?,С^Ь касатель”ого М°ДУЛЯ упругости от количества циклов нагружения

Ро е. (7—300 (а) и[ 900 МПа (б) при величине циклического напряже­

 

 

ния 1050

(/); 1470 (2);

1500 МПа

(3).

жрння гш

нр1пмпмТИ величииы пл°Щади петли гистерезиса

от количества циклов нагру-

Ы

:

еп Г " 14™

иапРяжения

1370 (1) и

1500 МПа (2) и максималь­

ной и минимальной деформации от количества циклов нагружения при величине цикли­ ческого напряжения 1370 (3) и 1500 МПа (4).

Рис. 5. Зависимости минимальной деформации (I),

остаточной деформации (2) и остаточной микро­ деформации (3) от напряжения циклической на­ грузки после 200 циклов нагружения.

На рис. 5 показана зависимость emin после 200 циклов нагружения от вели­ чины циклического напряжения (кри­ вая /). После разгрузки образца остаточ­ ная деформация уменьшается. Более рез­

кое снижение наблюдается

за первые

15 мин, но после 20-ч отдыха

при комнат­

ной температуре она уже практически не изменяется. На рис. 5 показана также е0Ст после 200 циклов нагружения и отдыха в течение 20 ч (кри­ вая 2). Остаточная деформация увеличивается при увеличении цикличе­ ского напряжения.

При рентгенографическом изучении образцов, подвергнутых цикличе­ скому нагружению, были обнаружены изменения по сравнению с исход­ ным состоянием (рис. 6), которые, как нам представляется, имеют пря­ мую связь с изменениями механических свойств материала в целом. Во-первых, было установлено наличие остаточной микродеформации 8ц0СТ (остаточная деформация в масштабе межплоскостных расстояний в мезоморфных областях волокон) при напряжении циклической на­ грузки а ^2^900 МПа. При меньших напряжениях ецост не наблюдалось. Как и е0ст образца в целом, ецост после разгрузки уменьшается, выходя на уровень 0,3—0,4%, после чего дальнейшего уменьшения практически не наблюдается. Зависимости е цост от напряжения циклического нагру­ жения за оч =900 МПа практически не обнаружено (кривая 5, рис. 5). Как и в [3], появление ецост мы связываем с конформационными перехо­ дами в макромолекулах волокон в мезоморфных областях при достиже­ нии определенного напряжения, что прямым образом влияет на характер деформирования материала, т. е. частично определяет описанное выше изменение зависимости о е материала и является одной из составляю­

щих е 0ст образца в целом.

При термообработке (100° С в течение 2 ч) циклически нагруженных образцов наблюдается частичное восстановление механических характе­ ристик, а также исчезновение ецост.

Рис. 7. Зависимость относительных из­ менений модуля упругости в сечении

кривой

а —е при

а= 900 U) и

700 МПа (2) и относительных изме­

нений

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

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

рефлексов (5) от циклического напря­ жения после 200 циклов нагружения.

Определенное значение в изменении зависимости а —е при цикличе­ ском нагружении, на наш взгляд, оказывает и вязкоупругий характер де­ формирования материала. Как видно из рис. 4, при данной скорости цик­ лического нагружения растет максимальная деформация, а минималь­ ная перед следующим циклом не успевает релаксировать до е 0ст- Следовательно, на молекулярном уровне материал все более ориентиру­

ется и модуль упругости растет.

Как видно из рис. 6, кроме сдвига дифракционных максимумов (изза появления едост) наблюдается изменение их интенсивности, что свиде­ тельствует об изменении степени упорядоченности структуры волокон. На рис. 7 показано относительное изменение интенсивности и относитель­ ное изменение касательного модуля упругости в зависимости от цикличе­ ского напряжения после 200 циклов нагружения. Как видно, в некоторой области напряжений имеет место благоприятное воздействие цикличе­ ской нагрузки на структуру волокна — степень упорядоченности растет. Очевидно, этот процесс связан с конформационными переходами в це­ пях, так как тоже начинается при достижении о ^~900 МПа. Выше напряжения 1100—1150 МПа, что при данной базе испытаний является самым благоприятным, эффект уменьшается, а при о > 1300—1400 МПа наблюдаем снижение степени мезоморфности. Последнее можно объяс­ нить разрушением в упорядоченных областях структуры — внутри фиб­ рилл. Очевидно, при напряжениях сг= 1100—1400 МПа над разрушением еще преобладает ориентирующий эффект циклического нагружения, а потом определяющую роль играет разрушение.

Изменения степени мезоморфности должны повлечь за собой и неко­ торое изменение механических свойств, например, модуля упругости, так как модуль упругости мезоморфных областей на 10—15% выше модуля волокна в целом [3].

Из рис. 7 можно видеть, что величина Е/Е0 следует за 1/10 с некото­ рым опозданием. Очевидно, дефекты, которые мы обнаруживаем в во­ локнах по спаду ///о, малы и оказывают малоощутимое воздействие на модуль упругости, что полностью перекрывается противоположным влия­ нием других факторов. Лишь при укрупнении дефектов при больших напряжениях а>0,88ар наблюдаем спад отношения Е/Е0.

Зависимость 1/10 от при напряжениях выше 1400 МПа указывает, что действительно происходит укрупнение дефектов, так как ///0 стано­ вится слабо зависимым от а~. Следовательно, новых дефектов появля­

ется мало и разрушение локализуется — увеличиваются наиболее опасные дефекты.

Выводы. 1. Малоцикловое нагружение вызывает значительные изме­ нения механических свойств однонаправленного органопластика _ из­ менение характера зависимости а —е и соответствующие изменения мо­

дуля упругости, площади петли гистерезиса, накопление остаточной деформации.

2. Рентгенографическими исследованиями установлены изменения в структуре армирующих волокон — появление остаточных микродефор-

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

фектов в структуре волокон.

^

3. Малоцикловое нагружение органопластика напряжениями

^0,88сгр вызывает увеличение модуля упругости. При напряжениях <7~ > 0,88сгр в начальном интервале модуль упругости также растет, но не достигает уровня, наблюдающегося при меньших напряжениях, и после определенного числа циклов начинается снижение модуля. Харак­ тер изменения модуля упругости может быть объяснен влиянием двух факторов: увеличение связано с конформационными переходами в цепях макромолекул и упорядочением структуры, снижение — с накоплением повреждений.

4.Величина остаточной деформации во время усталостного нагруже­ ния повышается с увеличением циклического напряжения и перед разру­ шением может достигнуть 1,5%.

5.Площадь петли гистерезиса при циклическом нагружении умень­ шается. Перед разрушением площадь петли увеличивается. Суммарная рассеянная энергия увеличивается с уменьшением циклического напря­ жения.

6.Предварительное циклическое нагружение на базе N= 103 циклов не оказывает заметного влияния на прочность однонаправленного орга­ нопластика.

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

1.Тамуж В. П., Куксенко В. С. Микромеханика разрушения полимерных материа­ лов. Рига, 1978. 294 с.

2.Вилкс У. /С. Устройство для измерения деформаций. Авт. свидетельство СССР

№ 355486. — Открытия. Изобретения. Пром. образцы. Товар, знаки, 1972, № 31. 3. Курземниекс А. X. Деформативные свойства структуры органических волокон на

основе параполиамидов. — Механика композит, материалов, 1979, № 1, с. 10— 14.

4.Перов Б. В., Скудра А. М., Машинская Г. П., Булаве Ф. Я. Особенности разру­ шения органопластиков и их влияние на прочность. — Механика композит, материа­ лов, 1979, № 2, с. 317—321.

5.Morgan R. Mones Е. Т., Stelle W. Deutscher S. В. The failure modes and

durability of Kevlar/49 epoxy composites. — SAMPE Quart., 1981, vol. 12, N 3, p. 26—31.

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

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

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

 

УДК 539.4:678.2

В.И. Соломатов, А. Н. Бобрышев, А. П. Прошин

ОВЛИЯНИИ РАЗМЕРНЫХ ФАКТОРОВ ДИСПЕРСНОГО НАПОЛНИТЕЛЯ НА ПРОЧНОСТЬ эпоксидных композитов

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

Настоящая статья посвящена изучению прочностных свойств ЭК и выявлению оптимальных значений размерных факторов наполнителя.

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

С целью выяснения влияния степени наполнения проводили прочностные испытания на образцах ЭК, изготовленных литьем при нормальном давлении из эпоксидной смолы

ЭД-20, наполненной молотым кварцевым песком, средний размер частиц которого рав­ нялся 2,38 -10-2 мм. Образцы отверждали при комнатной температуре в течение 24 ч и затем подвергали термической обработке при температуре 80° С в течение 8 ч. Вели­ чина степени наполнения / (процентное отношение массы наполнителя к массе смолы) варьировалась в диапазоне 0—600%.

Образцы испытывались на сжатие, изгиб и ударную циклическую прочность. В испы­ таниях на сжатие и изгиб использовались образцы с размерами ( 1 x 1 X1 ) • 10~2 м и (1X1X5) • 10-2 м соответственно. Скорость деформирования в этих испытаниях остава­ лась постоянной и равнялась 4 • 10_3 м/мин. В испытаниях на ударную циклическую прочность использовались образцы типа «грибок», форма которых была заимствована из работы [1].

Ввиду того, что частицы песка имели большую прочность, чем поли­ мерная матрица ЭК, разрушение носило смешанный адгезионно-когези­ онный характер. Результаты экспериментов представлены на рис. 1, 2. Как видно из графиков рис. 1, изменение прочности при изменении сте­ пени наполнения ЭК идентично (с незначительными отклонениями) во всех видах испытаний и соответствует данным работ {2—4].

Анализ полученных данных выявляет характерный спад прочности при изгибе с увеличением / в сравнении с прочностью при сжатии (рис. 1), что отмечалось в работе [2]. Такое несоответствие в прочности, а также наличие сдвига максимума выпуклой части (кривой 1 рис. 1) в сравнении с кривой 2 в область более низких значений I вызвано чувст­ вительностью изгиба к концентрациям напряжений, источниками кото­ рых являются локальные зоны у поверхности частиц наполнителя [3]. Из сказанного следует, что конструкции из эпоксидных композиций с дис­

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

Тем не менее, идентичность характера зависимостей (см. рис. 1, 2) свидетельствует о едином физическом механизме, ответственном за из­ менение прочности ЭК в зависимости от варьирования степени наполне­ ния. Основными параметрами, которые могут существенно повлиять на прочность при изменении /, являются величина общей площади поверх-

Моста $ частиц найбЛнйтеля й ёреднёё значенйё расстояния / между смежными частицами наполнителя. Эти параметры непосредственно свя­ заны между собой через величину числа частиц наполнителя и оказы­ вают на прочность как совместное, так и самостоятельное влияние. Действительно, изменение величины S оказывает влияние на изменение адгезионной прочности, тогда как значение I в основном влияет на коге­ зионную прочность полимерной матрицы.

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

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

/ = (CSp/M) 100% ,

( 1)

где С — постоянная величина; р — плотность наполнителя; М — масса связующего.

Изменение удельной свободной энергии ЭК, обусловленное смачива­ нием [5], равно

Y c == Y2.0 — Y i ,о + YU 2,

2)

где V2,o и 7i,o — поверхностное натяжение полимерной матрицы и напол­ нителя по отношению к воздуху; 71,2 — поверхностное натяжение между наполнителем и полимерной матрицей. Величины 72,0 и 71,0 являются по­ стоянными для определенной полимерной композиции. При больших значениях степени наполнения y'c = f(S ) ¥ zyc {5]. Согласно (2) опреде­ лим удельную работу адгезии в виде

иа = ии-~

(3)

Рис. 1. Зависимость прочности ЭК при сжатии сгс (/) и изгибе Он (2) от степени напол­ нения.

Рис. 2. Зависимости прочности (1—3) и относительной деформации (4—6) ЭК при ударном циклическом сжатии от степени наполнения. Данные получены при разрушаю­ щем числе циклов 106 (/, 4) \ 3,17106 (2, 5) и 107 (3, 6).

где Uk=272,0 — удельная работа когезии материала матрицы. Общую

удельную работу и0 адгезии и когезии ввиду их аддитивности найдем из уравнения

ио= 2ик- у с.

(4)

Подставляя ус из (2) в (3) и заменяя yi,o в соответствии с уравнением Неймана [6] значением 71,0 = 71,2+ 72,0 cos <р, получим известную зависи­ мость Юнга [7]: ^a= 72,o(l+coscp), где <р — краевой угол смачивания наполнителя связующим. Соотношение (4) является уравнением адгези­ онно-когезионной прочности ЭК. Оно показывает, что уменьшение по­ верхностного натяжения 71,2 путем введения в полимерный раствор раз­ личных пластификаторов и поверхностно-активных веществ [8—10] или путем улучшения технологии получения компаундов [11] однозначно при­ водит к увеличению общей работы wo> затрачиваемой на разрушение композита.

Условие термодинамической устойчивости ЭК при любых / с учетом

(1) и (3) запишется в виде [5]

dyc/dS>0\ Адус/д1>0\ Адиа/д1<0,

где Л = Ср • 100%]М — согласно (1) постоянная величина. Из (4) сле­ дует, что с ростом I прочность ЭК должна падать, что и установлено в результате обработки экспериментальных данных (левые ветви графи­ ков рис. 1). Однако, как видно из рис. 1, 2, некоторой части кривой (участок увеличения прочности) отвечает условие Адиа/д1> 0. Это усло­ вие соответствует существованию термодинамически метастабильных со­ стояний в ряде композитов, хотя механически они вполне устойчивы. Согласно [5] существует аналогия между характером изменения проч­ ности гетерофазных систем, подобных тем, что рассматриваются нами, и ван-дер-ваальсовыми изотермами. Действительно, меняя степень на­ полнения, можно заметить, что малонаполненные ЭК более склонны к вязкому разрушению, которое сопровождается пластической деформа­ цией и переориентацией макромолекул [12], тогда как высоконаполненные ЭК разрушаются при меньших деформациях хрупким образом. Су­ ществование этих двух, резко различающихся характеров разрушения, очевидно, должно предопределять наличие переходных состояний, кото­ рые являются метастабильными.

В практике использования ЭК большее значение имеет усиливающее действие наполнителя, т. е. зона метастабильных состояний, где Адиа/д1 = 0 и Ад2иа/д21<.0, так как при этом можно достичь значений

Рис. 3. Зависимости упругих постоянных ЭК от степени наполнения / и дисперсности d

частиц песка-наполнителя: динамического (/) и статического (2) модулей упругости от степени наполнения; коэффициента внутреннего трения от степени наполнения (3) и дисперсности (4) наполнителя соответственно; динамического модуля упругости от дис­ персности наполнителя (5).

Рис. 4. Зависимости прочности ЭК при изгибе (1) и сжатии (2) от дисперсности напол­ нителя. Точкам на графиках соответствуют данные экспериментов, сплошным линиям — расчетные значения, полученные по формуле (11).

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

Результаты испытаний на ударную циклическую прочность (кривые 13) свидетельствуют об увеличении прочности в области повышенных значений /, что было обнаружено в работе [13]. Этот эффект обязан тому, чтр при небольших величинах деформаций (кривые 4—6 рис. 2) цикли­ чески напряженные связи макромолекул поперечных сшивок в высоконаполненных ЭК деформируются главным образом в области упругости (кривая 1 рис. 3) ввиду их ограниченности в подвижности, не претерпе­ вая разрывов. Поэтому в целом такие композиты ведут себя как упругие тела. Описанный эффект имеет прямое практическое значение при экс­ плуатации циклически малонапряженных конструкций, при изготовле­ нии которых следует отдавать предпочтение более дешевым, высоконаполненным композитам.

Общее снижение циклической прочности (см. рис. 2) по сравнению со статической прочностью при сжатии (кривая 2 рис. 1) вызвано как разо­ гревом композита [14] и запаздыванием релаксационных процессов [15, 16], так и разупрочняющим взаимодействием микронапряжений и дест­ рукций предыдущих циклов нагружения с микронапряжениями текущих циклов [17].

2. Влияние дисперсности наполнителя на прочность ЭК. Известно, что увеличение дисперсности частиц наполнителя повышает прочность полимерных композитов [6]. Однако для практики большее значение имеет количественная связь между этими параметрами.

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

s = 4tt0,

где 0 = 20102/(01 + 02) — средняя кривизна поверхности частиц напол­ нителя; 0ь 02 — главные радиусы кривизны поверхности.

С целью установления количественной связи между 0 и прочностью ЭК необходимо иметь начальное краевое решение. Это решение, являю­ щееся условием адгезионного контакта плоской поверхности материала наполнителя с материалом полимерной матрицы, запишем в виде 0 = оо; Ua = Uoo. Условие, при котором 0 = 0, относится к чистому полимеру и ли­ шает смысла решение по отношению к композиту. Для всех остальных значений О < 0 < о о работу адгезии иа найдем из уравнения

ua = Uoo + Up,

(5)

где Uoo — удельная работа адгезии при условии контакта наполнителя и матрицы по плоскости; ир — удельная поверхностная энергия.

Отмеченный выше адгезионно-когезионный характер разрушения ЭК свидетельствует о разрыве у поверхности частиц наполнителя адгезион­ ного слоя, переходящего в структуру полимерной матрицы ЭК. С учетом

этого поверхностную энергию ир определим в виде

 

UP = WS4JI[®2—(0 —Л)2],

(6)

где Ws _ доля энергии ир, приходящаяся на одну молекулу адгезионного слоя; h — толщина адгезионного слоя. Общую работу адгезии и когезии согласно (4) —(6) определим из соотношения

w0N = WcoN+ ws\n [02 —(0 —h)2] + whNy

(7)

где WoN= u0\ WooN-u^-, WhN— Uk',w0, w„, Wk удельные значения общей работы, работы адгезии при 0 = 0 и работы когезии, приходящиеся на одну молекулу; N — число молекул в адгезионном слое. Дифференци­ руя (7) по N, получим

w0 = Wco + 2wsvld + Wh,

(8)

где d = 20 — средний размер частиц наполнителя;

о = 4/31/Л7лс{©3 —

—(©—/г)3] — объем, занимаемый одной молекулой адгезионного слоя. Во втором слагаемом в (8) определяли dQJdN, учитывая то, что 0 » /t.

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

к

dV = v d N = M 2 ^ Id (s'i —Si), i= 1

где l — средняя высота пирамид, составляющих многогранник частицы наполнителя; s'i >— площадь внешней поверхности адгезионного слоя на i-й грани; s* — площадь поверхности i-й грани многогранника. В ре­ зультате подстановки 2 (s' —s*) в (6) и (7) и выполнения преобразова­ ний, аналогичных тем, что производились для округлых частиц наполни­ теля, найдем соотношение для Wo>соответствующее уравнению (8).

Напряжение разрушения композита а связано с величиной w со­ гласно закону Больцмана [5] равенством

1/а = ехр( —wjkT),

(9)

где k — постоянная Больцмана; Т — абсолютная температура. Из начального краевого решения известно значение о<х>напряжения разру­ шения в зоне адгезионного контакта материала матрицы ЭК с плоской поверхностью (0 = оо, ua = Uoo) материала наполнителя. Исходя из этого определим напряжение разрушения ЭК, наполненного частицами со средним размером d, из условия нормировки (9):

оо

J exp (-wfkT)dw

1 / с = 1 / л . ---------------------

.

(Ю)

J exp(-wlkT)dw

Wоо

В результате интегрирования (10) с учетом (8) получим

а=ОооС exp (2wsv/dkT),

(11)

где C= exp(Wk/kT) — постоянная величина, характеризующая когези­ онное разрушение композита. Зависимость (11), как можно убедиться [19], аналогична по своей форме уравнению Томсона, выражающему из­

менение давления пара над каплей жидкости при изменении размеров последней.

Прочностные испытания ЭК с различной степенью дисперсности на­ полнителя проводились при сжатии и изгибе на образцах, размеры кото­ рых указывались ранее. Степень наполнения во всех ЭК оставалась по­ стоянной и равнялась / = 200%. Дисперсность наполнителя изменялась в интервале^ (0,0238—1,25) мм. В качестве наполнителя использовали кварцевый песок, который предварительно подвергался дисперсионному рассеву на размерные фракции. Скорость деформирования в испытаниях равнялась 4 - 10-3 м/мин.

Как видно из кривых 4, 5 рис. 3, модуль упругости и коэффициент внутреннего трения композитов не претерпевают значительных измене­ ний в диапазоне исследованных значений дисперсности песка-наполни­ теля. Результаты экспериментальных данных (рис. 4) в достаточной сте­ пени коррелируют с расчетными, полученными из уравнения (11). Из (11) следует, что при возрастании величины, обратной среднему значе­ нию размера частиц наполнителя d, на единицу разрушающее напряже­ ние а увеличивается в е раз. Температурная зависимость разрушающего напряжения по формуле (11) находится в соответствии с данными [20, 21], проверенными при 'различных температурах и скоростях деформиро­ вания.

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

1.Левандовский М. Я., Топоров Г. В. Методика исследования сопротивления оста­ точной деформации и контактному разрушению сталей при ударных циклических на­

грузках. — Завод, лаб., 1968, т. 34, № 12, с. 1498— 1501.

2. Ленг Ф. Ф. Разрушение композитов с дисперсными частицами в хрупкой мат­ рице. — В кн.: Композиционные материалы. Т. 5. Разрушение и усталость. М., 1978,

с.11—57.

3.Эйрих Ф. Р., Смит Т. Л. Молекулярно-механические аспекты изотермического разрушения эластомеров. — В кн.: Разрушение. Т. 7. Разрушение неметаллических и композиционных материалов. М., 1976, ч. 2, с. 104—390.

4.Гуль В. Е. Взаимосвязь между структурой и свойствами полимеров. — В кн.: Взаимосвязь структуры и свойств полимеров. М., 1975, с. 29—57.

5.Френкель Я. И. Поверхностные явления. — В кн.: Кинетическая теория жидкос­

тей. Л., 1975, с. 351—408.

аутогезии и когезии тел. —

6. Москвитин Н. И. Молекулярный механизм адгезии,

В кн.: Физико-химические основы процессов склеивания

и прилипания. М., 1974,

с.6—29.

7.Берлин А. А., Басин В. Е. Термодинамическая трактовка явлений адгезии. —

Вкн.: Основы адгезии полимеров. М., 1974, с. 52—85.

8.Малощук Ю. С., Киселев В. Я-, Воюцкий С. С. Влияние пластификаторов на

адгезию кристаллических полимеров к эластомерам. — Пласт, массы, 1967, 12,

с.68—71.

9.Хозин В. Г., Султанаев Р. М., Прохорова Н. С. Особенности модификации эпок­

сидных полимеров полихлордифенилом. — Механика полимеров, 1975, № 2, с. 218—223. 10. Маския Л. Добавки для модификации механических свойств. — В кн.: Д о­

бавки для пластических масс. М., 1978, с. 53— 117.

11.Ениколопов И. С. Композиционные полимерные материалы. — Природа, 1980, 8, с. 62—67.

12.Шишкин Н. И., Милагин М. Ф. Изучение предразрывных состояний твердых по­

лимеров. 2. Предел упругости и

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

1976,

№ 2,

с. 195—201.

 

 

 

13. Роузен Б. Феноменологическое исследование процесса разрушения

эластомеров

в стеклообразном состоянии. —

В кн.: Разрушение твердых полимеров. М.,

1971,

с.235—285.

14.Feltner С. Е. Cycle-dependent fracture of РММА. — J. Appl. Phys., 1967, vol. 9,

p.3576—3584.

15.Бугло С. T., Ратнер С. Б. Влияние релаксационных явлений на выносливость пластмасс при гармоническом и ударном нагружениях. — Механика полимеров, 1972, 1, с. 165— 168.

16. Bueche F. Tensile strength of plastics: effects of flaws and chain relaxation. —

J.Appl. Phys., 1958, N 8, p. 1231— 1234.

17.Степанов В. А., Никонов Ю. А., Беляева Л. А., Власов А. С. О причинах сни­ жения долговечности полимеров при циклическом нагружении. — Механика полиме­

ров, 1976, 2, с. 279—283.

18.Постников В. С. Физика и химия поверхности. — В кн.: Физика и химия твер­ дого состояния. М., 1978, с. 440—478.

19.Беккер Р. Давление пара над небольшой каплей. — В кн.: Теория теплоты.

М., 1974, с. 95— 100.

20. Гуль В. Е. К вопросу о разрушении полимерных материалов. — Механика по­

лимеров, 1975, 2, с. 195— 199.

21. Шамраевская Т. В., Щукин М. В., Гуль В. Е. К вопросу об оценке энергии активации по температурно-скоростным зависимостям прочности. — Механика полиме­ ров, 1976, 4, с. 591—596.

Пензенский инженерно-строительный институт

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

УДК 624.073.001:678.067

Б. А. Киладзе, И. Я. Преображенский, А. Ш. Цхведиани

КОЛЕБАНИЯ МНОГОСЛОЙНОЙ ЦИЛИНДРИЧЕСКОЙ ПАНЕЛИ

САНИЗОТРОПНЫМИ СЛОЯМИ ПРИ БОЛЬШИХ ПРОГИБАХ

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

нографиях [1—3].

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

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

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

Связующие слои считаются упругими, анизотропными, податливыми на сдвиг, а деформации в направлении нормали полагаются равными нулю. Толщина несущего слоя — ft; толщина связующего слоя — hc. Ко­ эффициент объемного армирования

ft

(i )

ft + ftc

Плотность материала несущего слоя — рн, материала связующего — рс. Оболочка считается тонкой и пологой, радиус срединной поверхности ра­ вен R ; полная толщина оболочки H<^R (рис. 1). На рис. 2 показано се­ чение в плоскости х —z элемента оболочки в процессе деформирования.

Рис. 1.

Рис. 2.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]