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

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

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

Наблюдаемые в опытах начальные несовершенства в виде пор и трещин сильно удлинены вдоль волокон, поэтому в дальнейшем пред­ положим следующее:

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

б) рост трещины не изменяет исходного отношения полуосей:

в) сопротивление росту трещим в матрице и на межфазной грани­ це при нормальном отрыве и сдвиге определяется соответственно дву­ мя параметрами къ кг и к01г к02 В случаях, когда эти коэффициенты равны между собой, возможно смешанное разрушение композиционной

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

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

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

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

Ы = 2 2 й,„«т8„>, I. к, в, п = 1,2,3.

(8.1)

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

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

131

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

Рассмотрим волокнистый материал с двоякопериодической систе­ мой одинаковых по величине и ориентации межфазных трещин.

Пусть / + /0 — межфазная граница; участок /0 определяет разрыв •сплошности — трещину, имеющую бесконечно малую ширину; I

.плоскость с совершенным контактом фаз; п — вектор упругих смеще­

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

Тпа = Тп у На — И> , Тб/,

(8.2)

Т^а = Тп = О, Т б Здесь т = ае1и — координата точки на границе; величины, относящи­

еся к волокну, имеют индекс а, к матрице — без индекса; знаки «+» и «—» указывают, что предельное значение определяется вдоль поло­ жительного и отрицательного направлений нормали. Если волокна одинаковы, а их центры расположены в узлах двоякопериодической сетки, то условия (8.2) дополняются требованием инвариантности со­ стояния отдельной ячейки к операциям трансляционной и поворотной симметрии, допустимыми для выбранной структуры. Условия перио­ дичности требуют

»№(* + <•>)=Ощ{х),

(8.3)

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

В частности,

 

<°«> = -Г - I

(8.4)

132

где V— объем ячейки; 5* — площадь грани, перпендикулярной к оси хг\ ©! — компоненты вектора периода о.

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

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

 

 

^

=

0. °,г = о

^ .

°№ = О т -ш

(85)

Для

волокна

соответствующие

равенства

следуют из

приведенных

при введении

индекса а у всех величин. Краевые условия

 

 

иа =

и,

оа =

а1г,

т е /;

г = а;

|О| < Ф0,

(8.6)

 

оа1г = 0,

а1г =

0, т 6/; г = а\

й0< |$ |< :п ;.

 

 

При

однородном

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

ячеек

следует принять

 

 

 

 

 

 

 

 

и =

 

гсозО + - ^ г - г з т й ,

г >•оо.

(8.7)

Методом разделения переменных решение строится элементарно:

«„ = « « + 2 рп{РпеШ + Рпе-‘Л

п = ‘ ^ (8.8)

и = 11,,, + чре1* + яре~л

+ 2

Р* (Чпе'"* + Чпе~‘"%

 

л=1

 

где р = г/а; а — радиус волокна;

и0,

— несущественные постоян­

ные жесткого смещения компонентов;

 

гг®___«гг..

 

133

Строим формальное решение задачи. Подчиняя функцию иа краевому условию

^ = - т ~ ж = ° - р = 1 : ♦ •< | * | < я *

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

 

2 (Я„е'”в - Рпе~ш ) = О, „ < | ОI < я.

(8.10)

л> 0

 

Это функциональное уравнение выразим через условие аналитическо­ го продолжения за пределы круга вспомогательной функции X (г)

комплексной переменной г = ре1<>

 

|Х + (2) =

у

Х „г\

 

| г | <

1,

 

*(г) =

 

я>0

 

 

 

 

 

 

(8. 11)

 

У

Х

- ! р ,

| г| > 1.

 

* “ (*)=

 

 

 

 

я > 0

 

 

 

 

 

 

 

Уравнение (8.10)

заменяется соотношением

 

 

 

 

 

Х+ (т) — X- (т) =

0,

 

 

 

(8.12)

если Х п = Рп\ отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.13)

или

 

 

 

 

 

 

 

 

 

 

X (г) = 2

 

^ ( х ) = 2

*»*"*•

 

 

л>0

 

 

 

 

л>0

 

 

Подчиняя функцию смещения матрицы краевому условию

<7*’- = ' § " | г

= 0 - Р =

1;

» . < | * | < Я .

 

приходим к функциональному уравнению

 

 

 

 

 

2

(Чпем ~ ч пе~ш ) -

Че‘* + де~№ = 0.

 

п > О

 

 

 

 

 

 

 

 

 

Если принять уп = Уп,

то можно

ввести

вторую вспомогательную

функцию

 

 

 

 

 

 

 

 

 

 

 

К+ (г) = 2

^

п +

?

4 .

 

| г | <

1,

У(г) =

л > 0

 

 

 

 

 

 

(8.14)

 

 

 

 

 

 

 

 

 

К~(г) = л>05 Кп^ +<?г’

 

|г|>1‘

Условие аналитического

продолжения

функции

У (2) будет

К+ (т) — К-

(т) = 0,

р =

1;

#0< | ^ | < я ,

(8.15)

134

следовательно,

 

 

 

 

 

 

 

 

 

 

К(г) =

у ( 4 - ) .

 

 

 

(8.16)

На участке совершенного контакта компонентов имеем систему

функциональных уравнений

 

 

 

 

 

 

 

[Х (т)+аУ (т)]+ -[Х (т ) +

аУ(т)Г = 0, р =

1 ;

 

 

[X (т) -

У (т)]+ +

[X (т) — У (т)]~ =

0,

 

(8.17)

где а = ОЮа. Первое уравнение

решаем с учетом предельного усло­

вия (8.9) и равенств (8.13), (8.15)

 

 

 

 

 

 

Х(г) + аУ(г) = а ^ г

+ ? -^ -|.

 

 

(8.18)

Второе уравнение решаем с помощью функции Племеля

 

где

 

X (г) — У (г) =

<2 (г) р (г),

 

 

 

 

 

 

 

 

 

 

 

 

р (г) =

]^1 — 2 г с о з +

г2;

<2 (г) =

+ (}22] +

... Н— "

4-

далее

 

 

 

 

 

 

 

 

 

 

1*+ (т) + 1*

(т) =

0,

|'9 |< д 0,

 

 

 

 

 

- 2 К Л

|* |< 1 ,

 

 

 

 

ц(г) =

 

 

 

 

 

 

(8.19)

 

 

2

^

-

'*»'•

 

 

 

Здесь

 

л = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\> = 1;

X, = — соз

Хп = Рп (со$ 00) — 2 соз й0Рп_, (соз #0) +

 

 

+ Рп(соз1?0);

л > 2 ;

 

 

 

Рп (созО0) — полином

Лежандра

первого рода я-го

порядка.

Общее

решение системы (8.17) будет

 

 

 

 

 

 

 

 

X (г) =

(дг + д

 

Р (а) <2 (г),

 

 

 

 

 

 

 

 

 

(*)•

 

Для определения неизвестных постоянных полинома (}(г) вос­

пользуемся предельными условиями,

вытекающими из равенств (8.14)

и (8.16):

(г)—^оо, У(г)-*?2,

откуда

С}0= — д,

 

= д.

Искомое

135

решение, удовлетворяющее всем условиям внутри ячейки, имеет вид

х ^ = - 4 4 - [

* * + $

4 + 1и* *- >г - * ) ] •

 

у(г) --= - у ^ Г | а

(<?2 + 9 т )

— |Л<г) (^"Г — 9)]

(8,20)

где значение ц определено формулой (8.9).

В данном случае вспомогательные функции непосредственно свя­

заны с функциями перемещений компонентов

 

иа = Х(г) + Х (г), и - У (г) + У (г).

(8.21)

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

энергии (1.24) и установим связь

между напряжениями

а?з и

средними деформациями (у12), (ухз):

 

 

 

- 4 1ш [ Р ' (г) +

У~Щ\ (<т12 — |а1з) А г= -~ <о12> {у12) +

4

<ст1з> <?и>-

I

 

 

 

 

 

 

 

Заменяя под интегралом

напряжения на усредненные

(ог°2— ш°3),

получаем

 

 

 

 

 

 

 

.

г

. < У 1 а —

4‘Т 1 3 > ( 1 —

+ " о “ +

*°~с)

^ 1 2

— “У1а)

«°и- * й - с ( 1

+ - у

 

°

'

 

 

 

 

 

[ -

СХ1+ О + 0

* г -

 

 

Используя второе представление

энергии

 

 

 

а

[ иГ (г) Аг =

2(1^ 4

а|ЦгТ (0„ (Т12> + 013<713>2

4 1т

'— 2612^23№2)(?1з))

и заменяя смещения и на их среднее значение (1.17), запишем формулы для эффективных модулей

с181 “ 0 1 - &

+ (1 +

о о/оа±

а 2'

(8.22)

 

М-23 =

^32 =

О-

 

 

Здесь верхние знаки относятся к (312,

нижние — к 0 13; X, = — соз

к =

 

 

 

 

Для весьма малых трещин отслоения, когда # 0 -* л, значение уп­ ругих модулей согласуется с формулой (1.26), когда на всей межфаз­

136

ной границе соблюдаются условия совершенного контакта. Если во­ локно полностью отслоилось от матрицы, что соответствует значению 1>о = 0, получаем наинизшее значение модуля сдвига 612 = <?13 =

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

найденными результатами.

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

и = а р8созЗй + а-—2- р з т Зд.

Решение задачи и для этого случая строится аналогично. Оконча­ тельно имеем

X

= 1+ а [ ^ +

' ? Т +

И-(2 )< Э з(г)|.

У(2>=

Т П Г [“

+ Я- г )

- 1* й Яг (г)]•

Здесь <7 определено формулой (8.9);

<Эз (г) =

Я [Рг (соз #о) +

(со8 *о) +

22] —

1

,

?1 (С 05 ’&с)

,

Р 2 (С 05 Ос) 1

~~Ч \уг

-1

р

1-

г%

]•

Для произвольного закона возрастания N напряжений получаем решение методом индукции

х (г) =

1?г" +

Я -рг + ц (г) V

(сов &„) (4 - - 92'-') 1 .

 

 

/=1

 

У М = Т Т »

[ ( ^ +

^ т г ) “ - и (г)

(соз «о) (4 - - /гг'-')].

Суммируем все решения при различных я и находим достаточно общий вид взаимодействия данного волокна со смежными

* ( 2> = 7 Т 1 г 2 к 2" +

-7 - + И *) 2 р п ч « * »») (^Г -

 

л=Т1

/=1

 

 

=

+ - 7 г ) - ^ г> ^ рп

ч

х

л = 1

/= 1

 

 

(8-231

137

Строгое решение поставленной задачи ищем по методу сшивок [18], согласно которому полагаем, что поле смещений и напряжений по фор­ мулам (8.23) должно быть согласовано на некоторой окружности г= = Не™, расположенной внутри ячейки, с функцией, определяющей взаимодействие в двоякопериодической структуре. Последнюю для простой гексагональной или тетрагональной решетки берем в виде ряда по эллиптическим функциям Вейерштрасса С (г)

У (г) =

Аг + 2

АкЪ(г)'*-11.

(8.24)

 

А=1

 

 

Постоянные Ак и <7 согласовываем равенством смещений

и напряже­

ний по окружности сшивки

решений

т0 = Не™ и <р (т0) =

У (т0).

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

^

51 а °

- Р

+

йр+иьЧь) +

С 5^[ ( ~о~

 

— ^1,/»)Чр “I"

 

р>0,к^1

 

 

 

 

 

 

 

+

( - %

- “

< 0 Я,} =

 

.

(8-25)

 

Яп +

2

 

(Рр+1,кЯн +

йр+1. ^ ) =

0.

 

Здесь апф — коэффициенты разложения (р + 1)-й производной ^-функ­ ции Вейерштрасса:

 

а к + Р

1 - О Ц А р

(Д+01

с * , = т т щ

1 + в/Оа

аГ-*СГва

«.

_ а*-к(цаа

/0. кфр,

Т + Щ (1 + «*.);

**

1 + 0/0о

 

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

= 2 2 Ъ Л р +

=

<8 2 6)

Р>1

 

 

Сохраняя в системах (8.25), (8.26) первые доминирующие слагаемые, приходим к формулам (8.22). Последнее указывает на полную согла­ сованность результатов, найденных методом последовательной регу­ ляризации (см. гл. 5) с доминирующими членами в точном решении задачи.

138

Более строгие зависимости

 

 

(видимо, не отличающиеся

от

 

 

точных в первых четырех зна­

 

 

чащих

цифрах)

определены

 

 

численно.

Результаты

иссле­

 

 

дований

модулей

сдвига

от

 

 

угла

 

раскрытия

трещины

 

 

2 (л —#0) для стеклопластика

 

 

(С/Оа =

25)

при ^ =

0,5 при­

 

 

ведены

на рис. 59,

где кри­

 

 

вые 3, ^соответствуют гекса­

 

 

гональной,

а /, 2 — тетраго­

 

 

нальной

упаковкам.

Кри­

 

 

вые 1, 3 определяют

отно­

 

 

шение

С13/С,

а 2,

4 — Сг12/(7.

 

 

Как

видно,

 

влияние

вида

 

 

упаковки материала с трещи­

 

 

нами

на

его упругие

харак­

 

 

теристики

несущественно,

ес­

 

 

ли ? ^

0,5.

Это

следует

из

 

 

поведения

кривых на

рис. 60

 

 

(обозначения

см. на рис. 59).

 

 

Отметим появление значитель­

 

 

ной

анизотропии

в

упругих

 

 

свойствах первоначально изо­

 

 

тропного на

сдвиг материала

 

 

при

однотипной

ориентации

 

 

всех

трещин;

в

частности,

 

 

модуль С13 слабо изменяется,

 

 

если трещины ориентированы

Рис. 61

Рис. 62

нормально к плоскости сдви­

га

при

 

их

раскрытии

 

 

вплоть

до

углов

я/4;

превышение

этих пределов

ведет к почти

скачкообразному падению сдвиговой жесткости материала. Модуль С12 до значительных углов раскрытия трещин снижается практически по линейному закону. Некоторые отклонения от закона прямой свя­ заны с изменением ориентации площадки сдвига. Для более плотных упаковок при $ = 0,85 наблюдается локальное влияние ближайших волокон на искривление кривой (рис. 61) при гексагональной и тет­ рагональной (рис. 62) структурах. Кривые /, 2 соответствуют измене­ нию С13/С и С12/С. Для малых углов раскрытия трещин выполняется

условие Н т - А — (01кО,-А) = 0; 0% — модуль бездефектного мате­

риала.

Особый интерес представляет сравнение результатов, получен­ ных на основе первого приближения по методу последовательной регу­ ляризации и строгого решения задачи. Результаты расчетов приведе­ ны на рис. 63 при Е = 0,7 для рассматриваемого стеклопластика; сплошные кривые построены по точному решению, штриховые — по

139

формуле (8.22). Кривые 1 и 2 соответствуют поведению 013/0 и 01г/0 с ростом трещин. В целом наблюдается удовлетворительное сог­ ласование результатов при достаточно высо­ ком значении Очевидно, что погрешность будет расти для более плотных упаковок.

Рассмотрим условия хрупкого разруше­ ния волокнистой среды при продольном сдви­ ге и росте трещин в поперечном сечении ма­ териала в соответствии с линейной механи­ кой разрушения [5]. Скорость роста упругой энергии Д№ объема ячейки при действии на­ пряжений (о12) и (о13> будет

д №

(а»)2 д

( Р \ ^ (сг13)2 д

 

 

(8.27)

где Р — площадь ячейки. Пусть энергия

поверхностного натяжения

V распределена равномерно по контуру трещины:

V =

4ауЪ0

(8.28)

Здесь у — экспериментально определяемая постоянная. При этом при­ нято, что рост трещин не нарушает симметрию структуры материала, поэтому одновременно движутся оба конца трещины.

Из условия хрупкого разрушения Гриффитса следует выражение для соотношения между критическими средними напряжениями

ж ^

= 4 -

 

 

т ё +

 

ж ) '

это уравнение удобно записать так:

 

 

 

 

 

( СТ1г) а___ г

( а 1 з ) а

__

1

 

(8.29)

 

аг

Т-

 

 

где

 

М .

,2_

_8ау _

з

/ м

2 = _8уо__а_/

р

дЪ0 [ О ъ ) '

0 -

Р

д

[ О и ) '

Уравнение хтрочности (8.29) — начало

роста

трещин — определяет

на плоскости средних напряжений

(о12), (а13> эллипс. Как следует из

кривых (см. рис. 59, 60), при 2 (я — т?0) <

я/4 податливость 1/013 для

данной ориентации трещин слабо изменяется с ростом размера дуговой трещины, поэтому а2 > Ь2, и эллипс сильно вытянут вдоль оси (а12 >. Если (ог13) = 0, то связь между углом раскрытия трещины и крити­ ческими напряжениями

К 2> - л / '

_д_У1_

(8.30)

*

* ял г...

 

140