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

книги / Технология и проектирование углерод-углеродных композитов и конструкций

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

 

 

 

 

 

 

 

Аналогичные

исследо­

J ,M .o tc/n z

 

 

 

 

 

вания проводились в

рабо­

 

 

 

 

 

 

те [98]. Поскольку в этом

 

 

 

 

 

 

случае

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

 

 

 

 

 

 

композит

 

структуры

2D

 

 

 

 

 

 

обладал

нехрупким

пове­

 

 

 

 

 

 

дением,

то было показано,

 

 

___________

что

трещиностойкость

ма-

о,2 0,4

териала

 

 

более

 

точно

о,в

о,в

£,мм

описывается

кривыми

со­

Рис. 4.29. Зависимость /-интеграла от пере­

противления,

которые

от­

ражают зависимость удель­

мещения

б-торцов образцов

 

ной

энергии

разрушения

 

 

 

 

 

 

от

продвижения

трещины.

Ф,кД&с/М‘

 

 

 

 

 

Рассматривался рост

 

тре­

 

 

 

 

 

щины для образцов с про­

 

 

 

 

 

 

 

 

 

 

 

 

дольной

 

и

 

поперечной

 

 

 

 

 

 

схемой

нанесения

надреза.

 

 

 

 

 

 

Результаты

экспериментов

 

 

 

 

 

 

показали, что в том слу­

 

 

 

 

 

 

чае,

когда

фронт

трещины

 

 

 

 

 

 

нормален к

плоскости

сло­

 

 

16 cL,MH

 

ев, энергия разрушения су­

 

 

 

щественно

больше за

 

счет

 

 

 

 

 

 

 

Рис. 4.30. Зависимость сопротивления от ве­

дополнительных

механиз­

личины трещиы, распространяющейся перпен­

мов

поглощения

при

 

раз­

дикулярно

к плоскости

слоев

 

рыве волокон. На рис. 4.30

 

 

 

 

 

 

приведена

 

кривая

 

со­

 

 

 

 

 

 

противления

 

Ф

 

для

 

 

 

 

 

 

alW = 0,3.

разрушаются,

 

 

 

 

 

 

УУКМ

 

 

 

 

 

 

 

как

правило,

хрупко

при

 

 

 

 

 

 

действии

 

 

циклических

 

 

 

 

 

 

нгрузок. При этом процесс

 

 

 

 

 

 

разрушения

идет сразу

же

 

 

 

 

 

 

на всех

структурных уров­

—-=------ -------------- -----

нях. В работе [99] иссле-

дована

усталостная

проч-

10

10

10

 

10

N

ность

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

Рис. 4.31. Кривые

Велера

при

цикличес-

материалов

 

различного

ком растяжении

для

двух

различных

строения

при

растяжении

омпозитов

 

 

 

 

образцов с частотой 2

и Ю

Гц. На рис. 4.31 приведены кривые долговечности для некоторых композитов (JV—число циклов растяжения).

Г Л А В А 5

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

5.1.Краевая задача механики структурно неоднородных сред

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

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

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

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

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

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

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

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

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

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

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

наряду с тензором структурных модулей упругости 0 „шп, тензора микроповреждений ш,ушп(ер>г)* зависящего от условий нагру­

жения,

 

 

 

а№

= e u r J rW m„„ ~ ‘"„,W AV > |£OT'

<5Л)

ГДе *ijmn

—единичный тензор четвертого ранга;

—единичный

тензор второго ранга; е^—инвариант

тензора

деформаций.

Таким образом, в качестве математической модели процессов квазистатического деформирования и разрушения в рамках такого подхода поставлена стохастическая краевая задача ме­ ханики композитов [79],

+

« 0 .

u t(r) I

-B ?W . °,/г)пА

=S»(r),

<J.3>

l ^ a

Ire Гa

 

 

ще [/—вектор структурных перемещений; X—вектор структурных объемных (массовых) сил; а—тензор структурных коэффициентов теплового расширения; АТ—однородный перепад температуры в

элементе конструкции; [/^—заданный на части поверхности

элемента конструкции вектор перемещений; 5°(г)—заданный на части поверхности вектор поверхностных сил. Материальные

величины X,® , со, а для композитов с однородными компонентами

(фазами) являются кусочно-постоянными функциями радиуса-век­ тора г и изменяются по закону

*(т) для г б Vm,

*(') =

для г б Vjу

где Vm и Vf—соответственно объемы, занимаемые матрицей и

волокном. Практически все реальные структуры удается описать с помощью аппарата обобщенных функций. Примеры такого описания для структур типа 3D и SD приведены в работах [31, 32]. Таким образом возможно аналитическое задание материаль­ ных функций, что особенно важно при разработке численных и аналитических методов решения краевых задач микромеханики композитов.

Основные уравнения (5.2) и граничные условия (5.3) краевой задачи для элемента конструкции из композита необходимо дополнить условиями на межфазной поверхности. Например, для идеального контакта матрицы и волокон

[{Л(г)]+=[£Л(г)1 [у..]+=(о.лу]-

Отличительной особенностью краевой задачи (5.2), (5.3) является содержащийся в физических уравнениях (5.1) тензор повреждаемости со четвертого ранга, описывающий структурное разрушение, стадию формирования и роста магистральной трещи­ ны, что отражается на равновесной нелинейной диаграмме деформирования наличием ниспадающего участка [15].

Обобщенное (в силу разрывности коэффициентов) решение задачи (5.2), (5.3) прямым путем (т. е. непосредственно) получить для практически важных случаев, как правило, не удается. Поэтому предлагается поэтапное решение исходной краевой задачи (табл. 5.1).

На первом этапе по заданным феноменологическим структур- - ным уравнениям состояния с учетом формы и взаимного расположения элементов структуры строится макроскопическая модель среды. Для этого последовательное решение уравнений системы (5.2) и (так как важно найти именно макроскопические

S5

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

Т в б л и ц а 5.1

Схема решения исходной краевой задачи

Исходная информация

Второй этап решения подразумевает отыскание «осредненного» решения, т. е. отыскание макроскопических напряжений, деформаций и перемещений. Краевая задача для макровеличин, соответствующая краевой задаче (5.2), (5.3), записывается в виде

да1 (г)

+ Х*.(г) = О,

 

 

-

в 1 ,п п ^ м -

(5.4)

 

W

. r = ^ ' лМ

(5.5)

| г 6 Г = 5! »

Поскольку на втором этапе микронеоднородная среда, заполняю­ щая область V, заменяется эквивалентной макроскопически однородной средой, для краевой задачи (5.4), (5.5) можно применять современные методы механики однородных тел [62].

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

точке М(г) тела V известны макронапряжения о*. =

Тоща в

элементарном макрообьеме AJvt выделенном вокруг точки Щг), структурные поля деформирования удовлетворяют уравнениям (5.2) и граничным условиям

V/') 'rA - W .

(5-6)

где Т.—вектор сил на поверхности Гд/и макрообъема.

Если заданы макродеформации е* = е.., то для определения полей <7„(г) и е^.(г) следует рассматривать уравнения (5.2) совместно с граничными условиями в перемещениях

Согласно [78] решения краевых задач (5.2), (5.6) и (5.2), (5.7) совпадают. Очевидно, эти решения не зависят от геометрии

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

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

Вероятность микроразрушения, в свою очередь, можно связать с вероятностью макроразрушения, т. е. вероятностью разрушения

элементарного объема AJv. Это составляет основную задачу заключительного, пятого этапа решения исходной краевой за­ дачи. Эта задача впервые рассмотрена в [78], где установлена связь между инвариантными мерами структурного разрушения (вероятностью микроразрушения) и инвариантными мерами мак­ роразрушения (вероятностью макроразрушения) в общем случае. В частном случае степенной зависимости распределения микро-

5.2.Метод осреднения стохастической краевой задачи микромеханики композитов на основе модели однородной среды

Необходимость прогнозирования упругих свойств углеродных матриц с учетом разупорядоченности их поликристаллической структуры и наличия пористости приводит к постановке сто­ хастических краевых задач механики композитов. К аналогичной постановке сводится задача о прогнозировании эффективных свойств углеродных волокон с аксиальной текстурой. Для иллю­ страции основных идей метода рассмотрим частный случай краевой задачи (5.2), (5.3)

Э*~°йг(Г) ~ ° ’ =

+

<«>

W|r(Г 'и 'г

«•»

где е*—заданный тензор макродеформации;

0 ^ и(г)—тензор мо­

дулей упругости в лабораторной системе, связанной с тензором моделей упругости С'утп в кристаллографической системе ко­

ординат формулами преобразования тензоров четвертого ранга,

e ijmn(r ) = Xia

Xf i Хту Хп2р'afr6'

(5 Л 0 )

где X..—косинусы углов между лабораторной и кристаллографиче­

ской системами

координат.

среды введем

Для учета

пористости микронеоднородной

случайную индикаторную функцию Д(г), определяемую соотно­ шением

1 ,

если

г е

Vj

 

А(г) = О,

если

г £

Vm,

(5.11)

где Vj и Vm—объемы, занимаемые соответственно порами и

матрицей.

Тогда соотношение (5.10) для микронеоднородных сред с учетом распределения пористости (5.11) записывается в виде

(5.12)

® ijm t№ ~ Xia Xjp Хщ ХпЗС 'аруб(1 Х^ '

В модели (5.10) предполагаются известными упругие свойства кристаллов (зерен) С'„тд, а в модели (5.12), кроме того,

известны моментные функции различных порядков распределения величины Х(г). Требуется определить эффективные характе­ ристики микронеоднородной среды и поля деформирования в указанных выше случаях на основании решения краевой задачи (5.8), (5.9).

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

 

 

 

 

 

 

 

3

 

 

 

 

 

С'..

=

а'6.6

 

+

Ь'(6. 6. +6 .6 . ) + 'V

ц'6. 6.6

6

 

, (5.13)

ijtnn

 

ij тп

 

4 im )п

in jm'

P=1

r

ip jp

mp np

 

 

 

 

 

 

 

 

 

 

 

 

 

где a '= CJJ, 5 '=

С4'4,/г '= CJ'J -

CJJ -

3C4'4.

В

случае

среды с

гексагональной структурой

 

 

 

 

 

 

 

С>ijmn

=

a '8 ij8 mn +

b'$ im 6jn + 8 in8jn )

+

 

 

 

 

 

 

 

+ У'51Ъ5ЦдтЪ8пЪ+ X'(8i38J38mn + V«As> +

 

 

+ р'(6. 6 J

 

+ 6.J6 ,6.

+ <5. 6J6 , +

<5..<5 ,6 . ),

 

(5.14)

 

r

v т J3

m3

(3 «3 jm

т

J3 пЗ

i3 m3

упу>

 

 

где а \ Ь \у ',к ',р '—материальные постоянные поликристалла гек­ сагональной структуры.

Подставляя (5.14) в (5.10) и осредняя полученное уравнение,

находим

 

 

 

 

 

 

 

 

С..

=<0..

) = а'6.6

тп

+ в'(6. 6.

+ 6 .6 . ),

(5.15)

ijmn

' ijmtv

IJ

“ ' im jn

in jmn

 

где для

рассматриваемого

 

случая

 

 

’ = (“'+ b

 

'+\ к)

 

г - (*■+& ' +1"')•

<5Л6>

Таким образом, тензор средних модулей упругости микронеодно­ родной среды определяется как тензор модулей упругости изотропной среды с модулями Ляме X' и /?'. Из уравнений (5.10), (5.14) и (5.16)

®ijmn ~ ®ijntn

^ ijm r ) ~ ^ ^/3^/3^тЗ^ пЗ

 

 

 

- T s V ifm n +8im8jn + W

+* W

m n +

+ 8 ijK i3 *пЗ ~

3 8 ij 8 m J + Р '^ а *пЗ *jm + 8 in 8j3 8 m3 +

 

+ * 3*тЗ 8]П+ 8im ^J3*j3*пЗ ~ 3 *to 8jm

3 ^toi

(5‘17>

Как видно из уравнений (5.17), для вычисления моментов связи различных порядков

<9,7т„ « в р ^ )> . < < W ')

и Т- Д'

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

^ /3 ' V ’ ( А/3 А>3 Ат 3 АлЗ^ И Т. д.

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

Из формул (5.18) нетрудно установить алгоритм, по которому можно вычислить средние значения произведения направляющих косинусов любого порядка. Так, средние значения произведения восьми направляющих косинусов определяются уравнением

 

 

 

 

+

+ S

im

8

+ 8 8

 

 

jnpgrs

in Ji

 

В общем случае

 

 

^ (3AI23 А«зз —Х1г> ~ (п + 1)!! а/,*2... г

.(5.19)

Таким образом, пользуясь формулами (5.18) и (5.19), нетрудно вычислить моменты любого порядка структурных модулей упру­ гости микронеоднородной среды с гексагональной структурой модели (5.10).

Для определения моментов' различных порядков с произволь­ ной анизотропией элементов структуры необходимо знать средние значения четных произведений направляющих косинусов. Сред­ нее значение произведения двух направляющих косинусов равно

^ V ” 3 дЦ6сф'

(5.20)

 

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

30 v ) Х