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

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

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

Функция (2.16) обладает следующим свойством. Если относительные макроперемещения V точек тела V* детер­ минированы (ц1 = и), то ф°(х)==0.

Для доказательства выберем начало координат (xt) в точке М — центре шара ДV При таком выборе системы

координат абсолютные перемещения точек сферы A'S сов­ падают с относительными перемещениями по отношению к точке М. Следовательно, составляющие %г вектора

тождественны флуктуациям составляющих т)1 вектора ц1

относительного макроперемещения точек сферы. По условию вектор TJ1 детерминирован: т]1О= 0. Но тогда (х') = О

под знаком интеграла (2.16). Следовательно, ср,(х) = 0, что и требовалось.

Другое свойство состоит в следующем. Пусть тело V* — односвязная область с конечными линейными раз­ мерами; М (х) — внутренняя точка тела У*, расположен­ ная на расстоянии Кг'^.1ф„ от границы 5 тела V*; отно­ сительные макроперемещения точек границы S детер-_ минированы; #*Дх, х ') — непрерывная при М (х')65 функция. Тогда при любом сколь угодно малом значении параметра е* справедливо равенство

ФI (х) = j я .ф(х> х') (*') dS' = °-

(2.17)

(S )

 

Действительно, в общем случае абсолютные переме­ щения у(х') точек М(х') границы 5 заданы в какой-ли­

бо произвольной системе координат (хс) с началом в точ­ ке О', не связанной с данным телом. Не изменяя суще­ ства задачи, можно выбрать новую систему координат (Xi) с началом в точке Л4(х) (осуществить параллельный перенос первоначальной системы координат). В результа­ те после переноса системы координат под знаком инте­

грала

(2.17) получим флуктуации относительных переме­

щений точек М{х')

границы 5 относительно точки М(х) .

Эти величины будут макроскопическими,

так как при

е*-^0

расстояния

от точки М(х)

до

любой точки

Af(x')G5 не меньше величин первого порядка малости.

Следовательно, вместо величин (х)

под знаком интег­

рала (2.17) будем иметь величины r)t?°(x). Но по условию относительные макроперемещения т)(х) детерминированы:

61

T)i°(x) = 0.

Тогда,

очевидно, ср°(х) =

0, что и требовалось.

Рассмотрим эргодическое свойство

перемещений. Пусть

О — начало

неподвижной системы координат; ОМ—отрезок

на оси хь

имеющий конечную длину

г; точки

М(х')

от­

резка ОМ

имеют

случайные перемещения

х* (х');

=

= < Xi ) — моменты первого порядка случайных перемеще­ ний, заданные на множестве реализаций в соответственных

точках М2(х) тел

Vz статистической модели V* (см. п. 3

гл. 1). Если составляющая

(х) вектора х(х) абсолютного пе­

ремещения точки

М (х)

детерминирована

(х* (х) = щ (х)),

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

 

 

 

(х) =

уг (х) = < х,-(х) >

(2.18)

Действительно, случайное абсолютное перемещение х,- (х) точки М(х) можно найти как предел суммы относительных случайных перемещений концов непересекающихся от­

резков Дкг длиной (/0е*)2, образующих в сумме отрезок ОМ:

 

Хг (х) = Иш У

=

и, (х).

 

 

 

т-*- оо

 

 

 

 

Применяем к этому равенству операцию осреднения по

множествам реализаций

в соответственных точках тел

статистической модели V*. Учитывая, что

(х/{) = <

(x/t) > ,

находим:

 

 

 

 

 

 

 

и, (х) =

lim

У [U; (x/t) vt (xk_t)]

Vi (x) v, (0).

 

m -*■oo

k~1

 

 

 

 

 

Здесь Vi(x)

— среднее

абсолютное

перемещение

точки

М (х); а,(0) — среднее

переносное

перемещение

той же

точки, равное перемещению точки О, так как в перенос­

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

отрезка ОМ остается

постоянной.

Точка О по условию неподвижна: Уг(0)=0. Поэтому

из последнего равенства следует равенство

(2.18).

Случайные деформации. Симметричная часть тензора

dy

 

 

 

 

дисторсий----

 

 

 

 

dx

,. =J_ (Ъи4-^i

(2.19)

 

11

2 ^ dxj

dxt

 

62

образует тензор случайных деформаций е = iapeaP. Анти­ симметричная часть

эквивалентна случайному вектору поворота элемента dV среды как абсолютно твердого тела.

В тензорных обозначениях геометрические соотноше­ ния (2.19) записывают в виде

8 = def х, или 8 = —

[ух + (VX)*b

(2.19')

где def — оператор деформации;

звездочка над тензором

означает транспонирование. Для симметричного тензора е, как известно, е* = е, для антисимметричного (о: со* = —а>.

В рассматриваемой далее линейной теории деформа­ ции величины еij с вероятностью Р= 1 малы. Под этим подразумевается малость почти всех реализаций этих ве­ личин в соответственных точках Mz тел Vz статистической модели У*.

Макро- и микродеформации. Макроскопическими, по определению, являются деформации элементов dl V тела V*£B2, микроскопическими — элементов dllV

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

Эргодичность случайных деформаций. Поскольку век­ тор перемещения как отмечено выше, задан на множе­ стве реализаций в соответственных точках Mz тел Vz ста­ тистической модели У*, аналогично задан тензор в слу­ чайных деформаций. Естественно возникает вопрос: при каких условиях те или иные статистические характери­ стики тензора деформаций можно определить по реализа­ циям в одном теле У*. Ответ на этот вопрос связан с по­ нятием эргодичности случайных функций (п. 3 гл. 1 и

п. 1 гл. 2).

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

63

= <eij>, то последние можно измерить эксперименталь­ но на одном образце композитного материала или вычи­ слить по реализациям в одном теле. В этом случае результаты расчета можно применить к отдельным кон­ струкциям или изделиям.

Пусть U— длина отрезка, первоначально параллель­ ного оси Xi, га — относительное удлинение элемента dxi. Тогда относительное удлинение отрезка в направлении оси Xi определится выражением

Если Vij при i=/=j есть половина относительного сдви­ га параллелепипеда, составленного из отрезков длиной

(2.20)

Из предыдущей формулы видно, что формула (2.20) справедлива как при f=й=/, так и при i= j.

Формула (2.20) аналогична формуле (2.4). Поэтому поле случайных деформаций имеет аналогичные эргодические свойства.

Положив vu (x) = ( еи (х) > и LWj (х, х') = < е*. (х)х Xe°if (хг) ) , нетрудно убедиться в том, что статистически

однородное в широком смысле поле ei; (х) эргодично по отношению к начальному моменту первого порядка v;j при конечном радиусе корреляции, если выполнено условие

lim -4 -

ГГ LiHj(x, х') dx^xl = 0.

(2.21)

h

J J

 

При радиусе корреляции г* = 0 имеем

Существенное различие между условиями эргодично­ сти (2.21) и (2.22) заключается в том, что первое из них может быть выполнено только в случае, если тело имеет

64

неограниченно большие размеры, а второе — для тел с

конечными

или

даже малыми

линейными размерами.

В отличие

от

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

условий эргодичности,

приведенных в теории напряжений, условия (2.21) и (2.22) можно проверить экспериментально на образцах композитных материалов, удовлетворяющих критерию статистической однородности в широком смысле. Соглас­ но этому критерию, поле средних деформаций иг:/(х) дол­ жно быть однородным (Vij = const относительно х) и моментные функции Lmj (х, х') должны зависеть только от разности (х'—х) векторов х 'и х (Lijij(x9x/)= L ijij(x/— —х)). Если поле деформаций изотропно, то Lijij(xi х') = = 1гЛ (|х'—х|).

Измерения геометрических характеристик расположе­ ния компонентов и случайных деформаций зерен квазиизотропных поликристаллов [5, 72, 84—87] показали, что радиус корреляции г* имеет сравнительно малую ве­ личину. В зернистой структуре он равен 3—4 средним линейным размерам зерен.

Пусть тело У* принадлежит к классу В2; Alh = li s* — линейные размеры элемента А1^.* Макродеформации е!/ и микродеформации г// связаны соотношениями

(Д1/;)

аналогичными формулам (2.8). Перефразировав соот­ ветствующие доказательства, нетрудно убедиться в том,

что поле микррдеформаций е!/(х), удовлетворяющей условиям локальной статистической однородности (усло­ виям макрогладкости по п. 1 гл. 2), имеет следующие свойства.

Чтобы статистически однородная в локальном смысле (макрогладкая) случайная функция ej} (х) была локально-

эргодической по отношению к функции vU (х) = < е^.1(х) > , необходимо и достаточно выполнить условие

lim

— р---- С Г Ц/' / (х, х') dx,dx: = 0.

(2.23)

«.-о

(А /г)2 J J

 

 

(д’/г)

 

5. Зак. 674

 

65

У словие (2.23) будет выполнено, если полож ить

'ФО,

\xi — x'l \<l°i г'+а

Lljijix, х')

(2.24)

= 0,

\xi— x 'i \ > l°tei+a

0<C a<l; 1°— характерные размеры области корреляции

функции е}/ (х). Предполагается,

что функции

(2.24)

ограничены по модулю.

 

имеет форму

шара

В случае, если область корреляции

(/] = /°), радиус шара

можно задать

соотношением

г* =

= /0еН-а. Величина г*

представляет

радиус корреляции слу­

чайной функции е” (х).

Если область статистической зависимости—эллипсоид или

имеет какую-либо иную форму, отличную

от шара,

то в

экспериментальных условиях обычно можно

задать

ее ха­

рактерные размеры величинами r%k =

 

аналогичными

радиусу статистический зависимости.

Постоянные

и а

можно подобрать в соответствии с экспериментом.

 

мик­

Например, по результатам эксперимента

[87]

для

родеформаций г\\ можно положить: /°

=

1 мм

(средний

размер зерна), /° ej+a = 4 мм (характерный размер области корреляции), 1\ = 50 мм (минимальный размер расчетной части образца). Решая два уравнения с двумя неизвестны­

ми, находим: е* = 0,14, а =

0,28.

 

 

Чтобы получить модель V* класса В2, нужно неограни­

ченно уменьшать параметр

е* при

неизменных

значениях

величин /° и а (п. 3 гл.

1).

При этом размеры тела будут

оставаться неизменными,

а

размеры элементов

структуры

в модели V* будут неограниченно

убывать с одновремен­

ным увеличением их количества в

одном данном элементе

AlV В пределе, при е*->0, моментные функции (2.24) будут удовлетворять условию (2.23). Следовательно, в мо­ дели V* класса В2 функция еЧ (х) будет локально-эргоди-

ческой по отношению к функции vjj (х).

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

66

Далее, пусть макроскопические деформации в|. данного тела V* детерминированы (е}. = еи). Тогда моментные функ­ ции микродеформаций порядка т

L}}(% 1,

 

xm) = < &IJ° (х,)

е '/* (хJ > ,

где е!Г = E].J vjJ,

удовлетворяют условию

Нш

1

I

IL ll” !l(x ‘-

п) Щ')

 

е*-0 (Л%Г

 

 

 

(Д1//)

 

 

 

 

dx\m) = О,

(2.25)

аналогичному условию (2.23)

 

совместности случай­

Дифференциальные уравнения

ных деформаций. В сплошных телах Vz статистической

модели

(п. 3 гл. 1) реализации ez тензора деформаций

е удовлетворяют известным

из

классической механики

сплошных сред дифференциальным уравнениям совмест­ ности деформаций. Пусть

 

0,

i -- j или / =--- т или i =

m;

 

в»*,

1,

ijm— четная перестановка

чисел 1,

2, 3;

 

1 ,

ijm— нечетная перестановка чисел

1 , 2, 3.

Полагая, что с вероятностью Я=1 тела Vz имеют совме­ стные деформации, можно представить условия совмест­ ности случайных деформаций тела V* уравнениями

d2e«p

= 0

(2.26)

бiay^m дхудх6

Случайный тензор второго ранга TJ = (riij),

 

— ®iav®iP6

дх/1х6

(2.27)

 

называют тензором несовместности случайных

дефор­

маций.

Переходя к тензорным обозначениям, уравнения

(2.26) и (2.27) запишем в виде

 

у X (у X е)* = 0;

(2.26')

5*

67

Т] =

— V X (v х е)*,

(2.27')

где у X означает ротор

(в данном случае от тензора вто­

рого ранга) [88].

Если статистическая модель У* принадлежит к классу В2, дифференциальные уравнения совместности (2.26) можно составить как для макро- (ц1), так и для микроде­ формаций (ц11).

При детерминированных макродеформациях тензор несовместности макродеформаций т]г детерминирован. Если, кроме того, макродеформации совместны, то т]1 = 0. При этом микродеформации могут быть случайными, совместными или несовместными.

3. ФИЗИЧЕСКИЕ УРАВНЕНИЯ

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

®ij =

^lj = ПиарОар,

(2.28)

записанного для составляющих тензоров о и е, или соот­ ношениями

а = 0

е, 8 = П

о

(2.28')

для самих тензоров. Смысл операции свертывания, обо­ значенной точками, расшифровывается выражениями (2.28).

В уравнениях (2.28') случайные тензоры четвертого ранга 0 и П представляют физические свойства среды. Поэтому в отличие от уравнений равновесия (2.2) и гео­ метрических соотношений (2.19) теории деформаций уравнения (2.28) принято называть физическими. Рас­ пределения постоянных упругости 0 и П исследованы в п. 5 гл. 1.

Макро- и микросвойства композитной среды. Соглас­ но основным свойствам статистической модели У* класса В2, модули упругости 0 1, экспериментально определяе­ мые на образцах из композитных материалов (на элемен­ тах ДХУ), являются макроскопическими. При е*-*0 эти величины — модули упругости элементов dlV статистиче­ ской модели У*. Соответственно величины 0 П, экспери*

68

ментально определяемые на образцах из компонентов (на элементах ДПУ), являются микроструктурными модуля­ ми упругости. При е*-Я) они становятся модулями упру­ гости элементов du V статистической модели К*.

Чтобы тензор модулей упругости 0 был детерминиро­ ван, необходимо и достаточно положить равными нулю дисперсии его составляющих. Свойства упругого тела 1/* класса В2назовем макродетерминированными, если свой­ ства элемента AlV с центром описанной сферы в точке М(х) тела V* детерминированы: 0 1(х)=С 1(х) независи­ мо от величины х при е*->0.

Тело К*, по определению, макроскопически однородно в широком смысле, если Cx(x)=const относительно х.

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

%В дальнейшем будем предполагать, что модули упру­

гости 0 п (х) — случайные функции

детерминированного

радиуса-вектора х точки М(х) тела

V* с конечными зна­

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

Тело V* будем называть макрооднородным в узком смысле, если поле тензора 0 п (х) статистически однород­ но в узком смысле (см. п. 3 гл. 1).

Тело К*, по определению, макроизотропно, если тен­ зор макромодулей упругости 0 1 — изотропный. Матрица составляющих изотропного тензора содержит три раз­ личных элемента, из которых два независимы вследствие соотношения

0 1 1 1 1 01122

= 202323

(2.29)

Компоненты макроизотропного композита могут

быть

как анизотропными, так и изотропными.

 

Локальная эргодичность

модулей упругости 0 п(х).

Среднее значение модуля упругости, вычисленное путем осреднения по координатам элемента AlV при е*-^0

0 „(x )= lim

ГQu (x)dV

(2.30)

е,-*0 Д V

.)

 

 

(4V)

 

69

микроструктурных модулей

упругости 0 п (х),

в общем

случае является случайным.

Оно не совпадает с макро­

скопическим модулем

упругости

0 х(х)

(подробнее об

этом см.

гл. 3) и

со средним

значением

Сп (х) =

= <0п (х)>,

определенным

на реализациях в телах Vz

статистической модели V*.

 

 

 

 

Случайная функция 0 п (х) эргодична по отношению к

функции Сп (х), если последнюю

можно

вычислить по

формуле (2.30). Поскольку эта формула аналогична фор­

мулам

(2.8) и (2.20),

условия

эргодичности

функции

0 й (х)

очевидны. В дальнейшем

будут

нужны следую­

щие свойства функции 0 п (х),

связанные с ее локальной

эргодичностью.

 

был детерминированным

(0у =

Чтобы тензор 0у(х)

= Су),

необходимо и достаточно

положить, что макро-

гладкая функция 0 п (х)

локально-эргодична по отноше­

нию к начальной

моментной

функции Сп (х) =<0п (х) >

первого порядка модулей упругости

0 п (х).

При

этом

0у = Су = С11 (х).

 

 

 

 

детерминирован и

Следствие. Если тензор

0 у ( х )

0Ц° = 0 П—Сп, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.31)

Действительно,

пусть 0^ = 0у — С11

Тогда по форму-

ле (2.30) находим

 

 

 

 

 

 

 

 

По условию тензор Qv

детерминирован, следовательно,

©у = 0. Из последнего

равенства получаем выражение

(2.31).

 

Справедливо также обратное утверждение: если функ­

ция 0 й (х)

локальио-эргодическая, имеет место равенст­

во (2.31) и, следовательно, тензор

детерминирован:

0у = СП

пусть

макрогладкая

неслучайная функция

Далее,

Н (х, х') имеет особенность г~2

или

г-3

в точке М(х)

тела

V*:

 

 

 

 

 

 

Н (х,

х') =

г 3/г(х,

х'); г = |х

— х'|;

|h(x,

x ')|< Q < o o ;

70