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

книги / Релаксационные явления в полимерах

..pdf
Скачиваний:
3
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

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

кристалличности.

Действительно, как видно из

соотношений

(VI. 76) — (VI. 79),

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

объемного и

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

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

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

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

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

Исходное равенство в методе вириального разложения выра­ жает связь между средними напряжениями и деформациями во

всем материале и в отдельных фазах:

 

(<*т) = 2 °а ( О

(ет) = 2 °а <em>

(VI. 80)

Выражая в первом из равенств (VI. 80) напряжения в фазах через деформации, а во втором — деформации через напряжения и учитывая, что связь между средними напряжениями (деформация’ ми) в фазах со средним по материалу должна определяться по­ стоянным тензором

' « ) = Аатп (оп); <е“ > = „ <е„) (VI. 81)

приходим к следующим ' выражениям для эффективных коэффи­ циентов упругости и податливости:

стп = 2 v^m qBqn>

smn = 2 °a.smqAqn

(VI82)

a

a

 

Матрицы Amn и Bmn находим

из решения задачи об

упругих

напряжениях и деформациях внутри сферического зерна, заключен­ ного в неограниченную матрицу. Задавая средние по материалу на­ пряжения (Ощ) и проводя усреднение по зерну и матрице вы­

численных напряжений (стт)« находим матрицу А%п. Аналогично

находим матрицу BLr

Для рассмотренного случая сферических включений малой кон­ центрации соответствующий расчет был выполнен Кривоглазом и Черевко [13], которые нашли:

r = /T ,[l-t>2

3 ( l - v i ) ( f t i — /Q

1

2 /C ,(l-2 v ,) + ^ ( I + v , )

J

H* = p,[l - t>s

15 (1 — V|) (Ц, — Ц2)

1

Pi (7 — 5v,) + 2(»j (4 — 5vj) J

(VI. 83)

(VI. 84)

Здесь индексы 1 и 2 относятся к матрице и включению соответ­ ственно, а V — коэффициент Пауссона.

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

Применительно к композиционным материалам метод самосо­ гласования был использован Будянским [14] и Хиллом [15], хотя основы метода были разработаны раньше применительно к вычис­ лению модулей упругости поликристаллов Хирши [16] и Крёнером [17].

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

принимается сферической, но и в том,

что толщина

переходного

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

равной нулю.

внутри

сферического

включения

Пусть средние

напряжения

равны:

 

 

 

 

 

К ) = Ата К

- Не.

И*) (а .)

(VI. 85)

Здесь считается, что фазы композиционного материала являют­ ся изотропными и характеризуются модулями упругости Ка и ця.

Принимаяпервый раз, что сферическое включение представ­ ляет собой зерно фазы № 1, расположенной в матрице с эффектив­ ными свойствами, а второй раз, что включение соответствует фазе № 2, получим:

(<Тт> - М я ю № , H i, к \ Ц’ ) (ffrt) + M m e № . Ц * Г , Ц*) (<Г„> (V I . 8 6 )

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

Используя решение Эшелби [18] для задачи о сферическом включении в изотопной матрице, в случае W-фазного материала

найдем [14]:

JV

2 . . [ i

+

( • ( $

-

 

(VI. 87)

а

 

 

 

 

 

N

 

 

 

 

 

s а - [ ■ + * ■ - № )

 

(VI. 8 8 )

Г

-

2 (4 — 5v*) .

 

1+V*

(VI. 89)

1-5 (1 - v*) ’

а "

8 (1 — v*)

 

Уравнения (VI. 87) — (VI. 89)

вместе

с

обычным соотношением

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

 

 

 

 

 

,

3К ' - 2 ц’

 

 

(VI. 90)

v ”

6К* + 2ц*

 

 

 

 

 

полностью определяют зависимости К* и р* от модулей упругости Ка и ра отдельных фаз.

Эффективные модули упругости армированных композиционных материалов, состоящих из изотропных компонентов, вычислялись этим методом Хиллом [19.]

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

упругими включениями, что ограничивает возможность использо­ вания этого метода, для матричных смесей. В то же время формулы (VI. 83) и (VI. 84) для двух рассмотренных предельных случаев приводят качественно к правильным результатам.

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

Метод, как показал Хилл [21], в конечном счете базируется на принципе минимума дополнительной энергии. Однако ниже будут изложены основные идеи метода в форме, близкой к трактовке Ха­ шина и Штрикмана.

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

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

и i ~ u i ~ и% г ц —

Ъ ц

ъ^ц, c t j k i — c i i k t с ц ы

(VI. 91)

Определим теперь тензор упругой поляризаций рц\

 

Pi! =

ац —

c°m i4i

(VI. 92)

Отсюда видно, что если ац и гц — те напряжения и деформации, которые устанавливаются в неоднородном материале при заданных граничных условиях, то

P t l = c ljkfinl

(VI. 93)

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

Vp — Vо “ "2 I (Н ц м Р ц Р

ы “ * Ptle li ~~ 2P*y8? / ) d V

(V I - 94)

Uo = j S ° y u d V ,

Нцыс'ым- I t iM

(VI. 95)

взятого по объему неоднородного анизотропного тела при условии

( « W w + P l / L - 0

(VI-96)

« î ( S ) - 0

(VI. 97)

Функционал Up оказывается стационарным, если неоднородные поля напряжений и деформаций связаны соотношением (VI. 93).

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

 

 

a2D

Кг) ^ Г (Ю <

 

1 +

<i -

 

a2i>2 (K

 

 

 

 

bîDp

 

— <и> <

 

I +

btvt (Hi — Ц2)

 

ЗЯа 4- 8Ца

Да-

1

,

2

К — 2ц0

15

ji0 (Хо Ч- 2ца)

1 -

atv!(Ki -

Кг)

(VI. 98)

 

b,D,

 

 

 

lUIL

 

(VI. 99)

I — biVi (Ц, — Ц2)

 

А>а = Ка

g"На (VI. 100)

Здесь Dx = v\Vz(x\ — х2) 2— дисперсия величины х.

Аналогичные формулы могут быть получены для армированных материалов [22, 23].

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

ма решения задачи, развитая в работах [28—30].

форме:

Запишем уравнения равновесия

(VI. 49) в операторной

Lu — — F,

L = VCV

(VI. 101)

где С — тензор модулей упругости.

Найдем такой оператор D*, который приводит к уравнению рав­ новесия для усредненных смещений:

V (и) = - F, L* = VC*V

(VI. 102)

Тогда тензор С*, определяющий перенормированный оператор L* будет представлять собой искомый эффективный тензор моду­ лей упругости.

На базе дифференциального оператора L построим интеграль­ ный оператор М, а также регулярные операторы М° и L0:

ML = - /, М° (L ) = —/, (М) L0 = —/

(VI. 103)

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

и' = MF — (М ) F = - (M L0 + 1) (u)=sR (и)

(VI. 104)

При помощи равенств (VI. 103) операторы М и 1° могут быть выражены через известный оператор М° и случайную составляю­ щую оператора L:

М =* 2 ХпМ \

L° = (L) 2 ( - 2

(**>)

(VI. 105)

0

О V 1

/

 

X = M°U, M°F = J О (г — Го) F (Го) dro

(VI. 106)

Подставляя выражения (VI. 105) в (VI. 104)', находим разложе­ ние для случайного матричного оператора:

00

л—1

&+£<я—I

 

R - 2 у'я; Yn -

- 2

<**> + s

U*> U 1) -

1

ft=i

k,i= 1

(VI. 107)

Здесь ядро G интегрального оператора M° представляет собой тен­ зорную функцию Грина, явный вид которой дается равенством (VI. 53), a L' = VC'V.

Переходя от смещений к тензору деформаций и учитывая, что согласно (VI. 102) оператор дифференцирования будет стоять со­ множителем непосредственно перед щ для каждого слагаемого ря­ да, можно записать аналогичное разложение для оператора Р, вве­ денного первым из равенств (VI. 107). Таким образом, выражение (VI. 107) вместе с (VI. 68) дает формальное решение задачи о вы­ числении эффективного тензора модулей упругости.

Аналогично можно показать, что ряд (VL107) будет определять эффективный тензор податливостей. Явное значение оператора X в этом случае будет определяться случайной составляющей опера­ тора совместности £.,-Wm = eijPehngVj4nSpqim и тензорной функцией Грина этого оператора.

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

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

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

Согласно корреляционному приближению для двухфазных изо­ тропных композиционных материалов будем иметь [31]:

 

 

alJK

(VI. 108)

 

3 { К ) 1

( З К + .)2

6

 

( / ( +

2 ц ) £)ц 1

(VI. 109)

5

( З К + 4ц> (ц>2 J

 

Учет всех многочастичных взаимодействий в сингулярном при­ ближении путем суммирования рядов по относительным флюктуа­ циям модулей упругости и использование функции Грина уравне­ ния равновесия дает для верхнего значения эффективных модулей упругости [32]:

K “ =

( K ) - D K ( о , / с 2

+

+ а ) " 1

(VI. 110)

Ц“ =

(И) “ £>ц

+ ViH, + Ь)~1

(VI. Ill)

в

ft— g-diX e/C -H n)-1

(VI. 112)

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

 

 

 

 

 

 

 

 

(VI. ИЗ)

 

т Н

 

^ - Ч ^ Н

'1

 

 

(VI. 114)

 

 

 

 

 

«

- < £ >

* - ф

( *

+

*

х

*(VI.+115)| >

Если в равенствах (VI. ПО) и (VI. Ill)

или (VI. 113)

и

(VI. 114)

положить Оц" = 0, то получаем точное значение эффективного объ­

емного модуля

 

 

D

 

 

 

 

 

 

К ' шт (

/

о ---------------------------------£ -------------

 

 

 

( V I . 116)

 

 

 

viKi + V-Ai + -jj-u

 

 

 

 

Это решение

впервые

 

было получено

Хиллом

иным

методом

[33]. Оно представляет собой единственный случай, когда удается найти точное значение эффективного модуля упругости при произ­ вольном пространственном распределении фаз.

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

К О Р Р Е Л Я Ц И О Н Н Ы Е Ф У Н К Ц И И У П Р У Г О Г О п о л я

Деформирование неоднородных материалов приводит к появле­ нию как регулярной, так и случайной составляющей упругого поля.

Связь между этими

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

(V I.67), из которых

следует, что структура случайного упругого

поля зависит как от топологии распределения компонентов в не­ однородной среде, так и от координатной зависимости регулярной составляющей поля напряжений или деформаций. Ниже будут рас­ сматриваться лишь однородные макродеформации, когда выполне­ ноусловие (оп) = const или эквивалентное ему требование (en) — const. Такой подход соответствует стационарным задачам статистической радиотехники.

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

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

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

Ниже рассмотрены следующие корреляционные функции:

иц{г ~ го) =

<«/М « / ( Го)>

(VI. 117)

Оц (r -

ro)в

М И/ ы )

(VI. 118)

Еам (г -

г„)

(в'ц (г) е'к1(г0)>

(VI. 119)

2 j/« (г ~

го) =

(а'цW < 4 Ы >

(VI. 120)

 

= (сmi М

(VI. 121)

SpqrS {Г ~

Го) =

SpgrsWЫ )

(VI. 122)

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

В соответствии с определением, бинарные корреляционные фун­ кции представляют собой тензоры ранга 2N, где N — ранг тензора (вектора), для которого определяется данная корреляционная функция. Поскольку корреляции находятся для двух точек г и г0, то, вообще говоря, корреляционная функция второго порядка должна зависеть от двух переменных: г и Го или г — Го и г + г 0. Од­ нако в случае статистической однородности из условия инвариант­ ности относительно преобразования переноса следует независимость корреляционных функций от г + г0. При этом имеет место услозие эргодичности, когда пространственное среднее совпадение со сред­ ним по ансамблю эквивалентных систем.

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

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

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

Будем для определенности рассматривать корреляционные функции тензоров модулей упругости.

Для расчетов оказывается весьма эффективной гипотеза о раз­ делении тензорной и координатной зависимостей [3,37].

(Г) = A ф (г) (VI. 123)

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

Рассмотрим отдельно координатную и тензорную зависимости

корреляционного тензора Ср,Д. Поскольку корреляционные связи по мере удаления рассматриваемых точек друг от друга ослабляют­ ся, функция ф(т), описывающая координатную зависимость корре­ ляционного тензора модулей упругости, должна иметь следующие предельные значения: ф (0)= 1, ф(оо) = 0. Кроме того, производная функции ф(г) при переходе через нуль должна претерпевать раз­ рыв, если неоднородность связана с резкими границами раздела между фазами. Если же неоднородность связана с плавным изме­ нением свойств, например вследствие случайных изменений плот­ ности однофазной среды, то производная от ф(г) при г = 0 должна обращаться в нуль.

Исходя из этого, обычно используют следующие функции:

ф (г)= е

р,

р= -£-

cp (г) =s e“ ncos Рр

. и -

2V

T (v)

<P(r)=e_,,î

(VI. 124)

(VI. 125)

(VI. 126)

(VI. 127)

где r ( v ) — гамма-функция, a Kv(p) — функция Бесселя второго рода от мнимого аргумента (функция Макдональда). При v = Va

кармановская функция (VI. 126) вырождается

в

экспоненту

(VI. 124), поэтому

выражение (VI. 126), так же

как

и (VI. 125),

представляет собой

обобщение экспоненты (VI. 124).

 

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

Функции (VI. 124) — (VI. 127) могут быть использованы и для описания текстурированных композиционных материалов, анизо­ тропия которых обусловлена формой включений. В этом случае необходимо провести замену [42]

 

р ->(2*?аГу *

 

 

(VI. 128)

Отсюда при а! = 02

= 0, а3 ф 0 получаем

текстуру

типа

пла­

стин (однородность по

координатам Х\ и Хч),

а

при ai

= 0,

a<i

= о3 Ф 0 — волокнистую текстуру. Аналогично

могут

быть

рас­

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

Перейдем теперь к анализу тензорной зависимости корреля­ ционной функции второго' порядка модулей упругости. Из равен­

ства (VI. 124) видно, что она определяется сомножителем Apfrls. Рассмотрим механическую смесь двух изотропных компонентов. Подставим в тождество

 

A p q r s

( c i j k l c p q r s ) = ( p i 1 k l c p q r s ) i p i j k l ) ( p p q r s )

(VI. 129)

значения СцЫсогласно

(VI. 29)

 

 

 

 

A pqrs = 9DKV im V pqrs +

6DKll {V iJklDpqrs +

Dijk lV pqrs) + AD^Dijk lDpqrs

(VI. 1-30)

где

 

 

(X'y') = t»,o2 (*, — x2) (yx—

y-i)

(VI. 131)

 

D xy

 

 

D x ! ' D x x =

( * i — X i Y

 

 

Если в выражении (VI. 131) провести свертку четырех внутрен­

них индексов, то перекрестное слагаемое пропадает:

 

 

A Hrs -

W ïrs = Ы к У а г * + «W

*

(VI. 132)

Отсюда

с помощью

соотношений

(VI. 27)

и (VI. 28)

находим

объемную

и девиаторную свертки:

 

 

 

 

 

 

A U s s ~ 2 7 D rc; А 1

ш

- 20° р

 

(V I. 133)

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