Статистическая механика композитных материалов
..pdfФункция (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г |
||||||
статистической модели 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у детерминирован: |
||||||
0у = СП |
пусть |
макрогладкая |
неслучайная функция |
||||
Далее, |
|||||||
Н (х, х') имеет особенность г~2 |
или |
г-3 |
в точке М(х) |
||||
тела |
V*: |
|
|
|
|
|
|
Н (х, |
х') = |
г 3/г(х, |
х'); г = |х |
— х'|; |
|h(x, |
x ')|< Q < o o ; |
70