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

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

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

11. Вакуленко А. А., Смотрин Н. Т., Чебанов В. М.

О связи

параметров Лоде

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

1967, №

2, с. 373—374.

12.Баев Л. В., Малинин Н. И. Исследование ползучести органического стекла (полиметнлметакрилата) при совместном действии растяжения с кручением. — Меха­ ника полимеров, 1966, № 5, с. 671—677.

13.Вилкс У. К., Крегерс А. Ф. Сложное нагружение нелинейно ползучего полимер­ ного материала. — Механика полимеров, 1967, Nq 3, с. 421—426.

14.Васин Р. А., Громова Н. А., Никиточкин А. Н., Огибалов П. М. Эксперимен­

тальное исследование полиэтилена при сложном нагружении. — Механика полимеров, 1974, № 1, с. 10— 17.

15.Вилкс У. К., Крегер А. Ф. Маложесткий трехкомпонентный измеритель дефор­ маций трубчатых образцов. — Механика полимеров, 1975, № 5, с. 949—951.

16.Крегер А. Ф., Вилкс У. К. Аппроксимация семейств кривых физически-нели-

нейной ползучести полимерных материалов на ЭВМ. — Механика полимеров, 1976,

2, с. 220—229.

17.Дзене И. Я., Крегерс А. Ф., Вилкс У. К. Особенности процесса деформирова­ ния при ползучести и повторной ползучести полимеров в условиях одноосного растя­

жения. Часть 1.; ■— Механика полимеров, 1974, № 3, с. 399—404.

18.Дзене Й. Я., Крегерс А. Ф., Вилкс У. К. Особенности процесса деформирования при ползучести и повторной ползучести полимеров в условиях одноосного растяжения. Часть 2. — Механика полимеров, 1974, № 4, с. 589—598.

19.Дзене И. Я., Путане А. В. Коэффициент Пуассона при одномерной ползучести полиэтилена. — Механика полимеров, 1967, № 5, с. 947—949.

20.Дзене И. Я. Изменение коэффициента Пуассона при полном цикле одномерной ползучести. — Механика полимеров, 1968, № 2, с. 227—231.

21.Дзене И. Я-, Путане А. В. Коэффициент Пуассона полиэтилена при обратной ползучести. — Механика полимеров, 1968, № 3, с. 563—564.

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

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

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

 

УДК 539.4:678.067.4

Е. Н. Сахарова, А. С. Овчинский

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

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

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

Рассмотрение сил, действующих на движущийся элемент некоторого волокна в осевом направлении, приводит к уравнению

d 2Ui

(1)

d N ji-p fd z —^ -----Tidz = 0,

где t= 0, 1,2,3. — индекс волокон. Укладка волокон

принимается гек­

сагональной. Индекс 0 соответствует разрушившемуся

волокну; индекс

1 — волокнам первого пояса, непосредственно окружающим разрушен­ ное; индекс 2 — волокнам следующего пояса и т. д., индексы 3, 4М5 .. присваиваются волокнам по мере удаления их от разрушившегося! [3].

Сделав ряд предположений, характерных для одномерныхмоделей, исследующих перераспределение напряжений только вдоль одной осе­ вой координаты, а именно, предположив, что осевые напряжения рас­ пределены равномерно по сечению волокон (ofi = Nfi/Ff) , приняв во внимание дщ /дг = е_а и что волокна вплоть до разрушения деформиру­ ются упруго е/ = сг//£7, приведем уравнение (1) к виду

д2иг

рf

д2и{

Ti

(2)

~dz2

 

dt2

~ EfFf

 

 

где Ui — осевые перемещения

некоторого волокна; р/ — плотность;

Ff — площадь поперечного сечения волокна; Ef — модуль упругости волокон; Ti — касательные силы, действующие на элемент некоторого

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

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

 

Г(= | т,г,М,

(3)

О

 

где г/ — радиус волокна; 0 — угловая координата.

Касательные напряжения связаны со сдвиговыми деформациями матрицы, которая в общем случае может деформироваться упругоплас­ тически. Сдвиговые деформации матрицы между волокнами, yij согласно

известной гипотезе

[ 1—3] прямо пропорциональны разности перемеще­

ний этих волокон

(Ui — Uj) и обратно пропорциональны

расстоянию

между волокнами

[3] 6 (0):

 

 

yij={ui-U j)/b{Q ).

(4)

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

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

T = G my при у^ ут ‘, (5) т = С гу+ {Gm — GT)yr при у>ут, (6)

где ут=тг/б?т; хт — предел текучести матрицы на сдвиг; Gm — модуль сдвига матрицы. Б случае, когда вслед за пластическим деформирова­ нием матрицы на сдвиг возникнет процесс упругой разгрузки, имеет место соотношение

x = G my — (Gm—GT) {утш\ — ут)\

(7)

здесь ушах — соответствует величине

сдвиговой деформации

матрицы

в момент начала разгрузки.

 

 

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

(5) — (7) в

(3), а затем с учетом

(4) под­

ставив Тг в (2), получим

три системы

уравнений. Количество

уравне­

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

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

стадиях деформирования матрицы на сдвиг:

 

 

D 2u0—а у2и0+ (x y2U \= 0 ;

С£ 2

2

ф 2

сс ^

D 2U \-\— р —Uo— —- ocy2WiЧ— ^— u 2-\— -p Из = 0 ;

 

 

b

o

3

b

 

 

 

 

 

(8)

D2U2+

иi - <xy2u2+

u3 = 0;

D2u3 +

иx

u2 - а угиъ= 0;

3

3

b

3

Заметим, что полные деформации волокон e,i = Ef + dUi/dz, где е/ — уро­ вень деформаций в момент разрушения.

Если в процессе движения разрушившегося волокна касательные напряжения на его границе превысят значение сдвиговой прочности связи волокон и матрицы, то начнется процесс отслоения разрушивше­ гося волокна от матрицы. На этом участке z ^ z 0 движение будет опи­ сываться следующей системой уравнений:

 

D2u0 = KRa y2uT sign (й0) ;

 

 

 

D2u, + ^ - KRUT sign (йо) -

uv+ ^ ~

и2+

иг= 0;

(11)

 

О

1 6

 

о

 

2

2

2

 

2

 

D2UQ-\—

щ — <Ху2и,2-\— ih = 0;

D2u$-\

u\Н— ^— и.2— (Ху2и$= 0.

о

о

о

 

о

 

Здесь, как и в работе [2], предполагается, что участок волокна испы­ тывает сухое трение о матрицу с постоянной интенсивностью сил трения Тд= /СяТг, где KR — некоторый коэффициент.

Системы уравнений (8) — (11) принадлежат к классу линейных ги­ перболических систем уравнений второго порядка. Совместное их ре­ шение с учетом условий стыковки по границе упругих и пластических деформаций, а также зоны отслоения было получено численно с по­ мощью метода характеристик, который позволяет получить решение, не обладающее числовой дисперсией вблизи фронта распространения волны разгрузки, что наблюдается, например, при использовании ме­ тода сеток [2]. Расчеты проводились по программе, составленной на языке ГДР-АЛГОЛ на ЭВМ БЭСМ-6.

З.ДИсследовалось перераспределение напряжений при разрушении волокон в бороалюминии./Исходя из реальных свойств компонентов,

 

Рис. 3.

 

Рис.

4.

Рис.

3. Перераспределение перемещений (а) и напряжений (б) в разрушившемся

 

волокне для двух моментов времени.

 

Рис.

4. Перераспределение перемещений

(а) и напряжений

(б)

в волокнах первого

 

пояса для двух

моментов времени.

 

 

качестве исходных данных использовали параметры E = Ef/Gm= 14, Yr = 0,0037; q2= G m/GT = 20; KR = 0,5. Объемная доля волокон У/ (в дит через Ь), значение сдвиговой прочности связи волокон и матри Tjb/тг, параметр Т1 = е//ут, характеризующий уровень нагрузки, при торой происходит разрушение волокна, варьировались в широких п делах. Результатами расчетов являются эпюры распределения осев перемещений Ui(t,z)/uT и растягивающих напряжений Ofi(t, z )/а/00 длине волокон и во времени, для различных сочетаний парамет] модели.

На рис. 2—6 сплошными кривыми показано распределение пере: щений и напряжений при упругопластическом деформировании м рицы на сдвиг без отслоения; штриховые кривые характеризуют вл ние отслоения на процесс деформирования волокон. Штрихпунктирг кривые на рис. 2 и 3 соответствуют результатам решения для «од волоконной» модели, когда волокна, соседние с разрушенным, закр лены [2]. Кривые на рис. 2—6 построены для значений парамет; У/ = 0,5; Tjb/xr=l,2; т|= 1.

П а

рис. 2 представлены перемещения конца разрушившегося

 

,1локна.

Прямая

1 соответствует решению без учета динамических

:

II фектов

[3]. Из

рисунка видно, что кривые 2 и 3 близки между соб

п т. е. рассмотрение «многоволоконной» модели не вносит существенн поправок в поведение разрушившегося волокна, что свидетельств; о правомерности исследования процессов перераспределения напря] ний в разрушенном волокне в рамках «одноволоконной» модели [ Отслоение разрушившегося волокна от матрицы приводит к увеличен перемещений конца волокна (кривая 4 рис. 2), и, в результате дейсп сил трения, колебательного процесса не возникает. |

Т На рис. 3—6 приведены эпюры распределения осевых перемеще! Ui(t,z)/uT и соответствующих напряжений Ofi{t, z)/a/°° по длине раз шившегося волокна, по длине волокон первого, второго и треть j поясов для двух моментов времени t/tf = 4 и 16 (обозначены на рис. 3

!звездочкой). Здесь tf = dfj^Eflр/ — время прохождения волны разгру: на величину диаметра волокна df, a/°° — уровень напряжений в вол нах в момент разрушения.

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

.UJ/UT

 

 

Рис. 5.

 

Рис.

5.

Перераспределение перемещений

(а) и напряжений (б) в волокнах вто

Рис.

6.

пояса для двух моментов времени.

Перераспределение перемещений (а) и напряжений (б) в волокнах трет

 

 

пояса для двух

моментов времени.

вслед за волной разгрузки (см. рис. 3—б). Перегрузка волокон, сосед-/ них с разрушившимся, имеет также сложный, динамический характер.] Г~Во-первых, можно выделить перегрузку волокон вблизи плоскости раз­ рыва (0^г/с(/^2,5)/как это было показано при статическом анализе [3] jtBfTBTopbix, волны перегрузки, которые формируются в плоскости разрыва и пробегают по значительной длине волокон, соседних с раз­ рушившимся. При этом для волокон второго и третьего поясов пере­

грузки второго типа выше первых^ Т Анализ соответствующих зависимостей (штриховые кривые) в слу­

чае отслоения разрушившегося волокна от матрицы показывает, что для начальных моментов (t/tf = 4) наличие отслоения вносит несуще­ ственные различия. Но в последующие моменты времени характер про­ цесса перераспределения напряжений принципиально меняется. В ре­ зультате действия сил трения в волокнах второго и третьего поясов воз­ мущение полей перемещений и напряжений на участках, прилегающих к плоскости разрыва, быстро затухает. Для волокон первого пояса наблюдается снижение уровня перегрузки, в плоскости разрыва вместо пика напряжений образуется площадка, длина которой примерно равна зоне отслоения. Развитие процессов отслоения волокон от матрицы приводит к уменьшению вероятности их разрушения{ Эти результаты важны при исследовании процессов разрушения композитов в вероят­ ностном аспекте [4], а также при моделировании механизмов разру­ шения на ЭВМ [5, 6].

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

Параметром, необходимым для моделирования процессов разруше­ ния на ЭВМ, является коэффициент максимальных дополнительных перегрузок [6] A/C*= (a/i(Z=0) — а/00)/'^00. На рис. 8 приведены зависи­ мости АКг от объемных долей волокон с учетом перераспределения части нагрузки на матрицу [6]. Полученные значения динамических перегрузок примерно вдвое превышают значения, полученные при ста­ тическом анализе [3].

Рис. 7.

Рис. 8.

Рис. 7. Распределение максимальных перегрузок по длине волокон первого (/), вто­ рого (2) и третьего поясов (3).

Рис. 8. Зависимости коэффициентов дополнительных перегрузок от объемных долей

волокон Vf для волокон первого (1),

второго (2) и

третьего поясов (3).ri = 0,3 ( ---------

)

и

0,5 (-----------

).

 

 

Численное решение было получено [9] для циклически стабильной матрицы, обладающей линейным упрочнением, для которой (Ти=Ф (еи) имеет вид:

Г ЗСтпбц

еи^ег;

° и '•3G7ne„(l — X )+3G mXer

еи>Ег.

Из сравнения соответствующих эпюр работ

[2, 9] следует, что учет е™

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

т ~ у качественно не меняет

характера процесса пере­

распределения

напряжений, приводит к некоторому запаздыванию в

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

Выводы. 1."|Па основе одномерной модели приведен анализ динами­ ческих эффектов перераспределения напряжений как в .разрушившемся волокне, так и в волокнах, соседних с разрушившимся^

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

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

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

4. |Развитие процесса отслоения разрушившегося волокна от мат­ рицы приводит к уменьшению перегрузок соседних волокон. J

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

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

2.Сахарова Е. Н., Овчинский А. С. Динамика перераспределения напряжений в разрушившемся волокне композиционного материала. — Механика композитных ма­ териалов, 1979, № 1, с. 57—64.

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

деление напряжений при разрыве хрупких волокон в металлических композиционных материалах. — Механика полимеров, 1977, № 1, с. 19—29.

4.Тамуж В. П. Объемное разрушение однонаправленных композитов. — Меха­ ника композитных материалов, 1979, № 2, с. 260—267.

5.Копьев И. М., Овчинский А. С., Билсагаев Н. К. Моделирование на ЭЦВМ

процессов разрушения композитов с дефектами прочности связи между компонен­ тами. — Механика композитных материалов, 1979, № 2, с. 217—221.

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

7.Ильюшин А. А. Пластичность. М., 1948. 376 с.

8.

Москвитин

В. В. Пластичность при переменных нагружениях. М., 1965. 310 с.

9.

Сахарова Е. Н. Исследование динамического напряженного состояния при раз­

рушении волокна

композиционного материала. — В кн.: Упругость и неупругость,

1979,

вып. 6, с. 85—90 (М.).

Институт металлургии им. А. А. Байкова

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

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

 

УДК 539.4:678.067.4:539.2

В.И. Костиков, С. А. Колесников, Е. И. Холодилова, С. В. Вавилки

А.Н. Полилов, В. К. Хохлов, Г М. Бутырин, А. П. Набатников,

И.П. Шепилов, А. Н. Шуршаков

ВЛИЯНИЕ СТРУКТУРЫ НА СВОЙСТВА ВОЛОКНИСТЫХ УГЛЕРОДНЫХ МАТЕРИАЛОВ*

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

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

Ор = Omax@~bPi

 

 

 

ГДе Gmax

ПРОЧНОСТЬ

 

бб(

 

 

 

ристого материала;

Р

 

 

 

ристость матрицы; 6 = 0,405

 

 

 

+ 0,318^ + 1,22

[4]

(/,

d

 

 

 

длина и диаметр поры с(

 

 

 

ветственно).

 

 

 

 

 

 

 

 

Испытания

образцов

 

 

 

струе

плазмы воздуха

пок<

 

 

 

ли, что унос массы углер

 

 

 

углеродного

материала

 

та:

 

 

 

связан

экспоненциальной

 

 

 

висимостыо с пористостью

 

 

 

териала

[5].

Уровень

 

по|

Рис. 1. Интегральные кривые распределения пор

тости

углерод-углеродных

териалов

определяется

 

те;

по размерам для композитных материалов на

логическими

режимами

и

основе высокомодульного волокна (1), низкомо­

дульной ткани (2) и

нпзкомодульного

во­

товления,

типом

волокна

локна

(<?).

 

связующего.

В

литера!

Доклад, представленный на IV Всесоюзную конференцию

«Структура,

свой

и применение композиционных материалов»

(Москва, ноябрь 1978 г.).