CC BY
47
УДК 539.3
МОМЕНТНЫЕ ФУНКЦИИ СТОХАСТИЧЕСКОЙ КРАЕВОЙ ЗАДАЧИ СТРУКТУРНОЙ МЕХАНИКИ МАТРИЧНЫХ
КОМПОЗИТОВ
Ю.О. Аристова, Н.В. Евлампиева, А.А. Ташкинов (г. Пермь)
Abstract
In this work a stochastic boundary-value elasticity problem is formulated for structural-inhomogeneous media taking into account the moment functions which are capable of describing the real material structure. The appropriate analytical solution has been given to determine the field moments of structural stresses and strains of a composite consisting of randomly located spherical inclusions and a matrix. Numerical modeling has been performed to obtain the part of the structure in terms of the Monte-Carlo method, to statistically analyze this structure and to derive its moment functions of the second and third order. It has been found that these functions can be approximated analytically, which is necessary to calculate the moments of random fields of structural stresses and strains.
В настоящее время для построения механических моделей композитов широкое распространение получил структурно-феноменологический подход, основанный на том, что феноменологические уравнения и критерии механики деформируемого твердого тела рассматриваются на двух уровнях: структурном, в рамках которого описывается поведение элементов структуры, и макроскопическом, описывающем композиционный материал как однородный с эффективными свойствами.
В связи с этим удобно вместо всего объема композита ограничиться рассмотрением представительного объема I-го порядка малости V с характерным размером l таким, что l << L, где L - характерный размер конструкции [8].
В пределах элементарного объема 2-го порядка малости Vk с характерным размером l k таким, что l k << l, композиционному материалу присваиваются свойства элементов структуры. На этом уровне представляется возможным исследовать процессы с помощью моделей и методов механики деформируемого твердого тела. При этом результаты исследований используются в континуальных уравнениях макроскопически однородной среды с помощью некоторых осредненных параметров.
Таким образом, в рамках структурно-феноменологической модели стохастическая краевая задача механики композитов в отсутствие массовых сил записывается в виде:
foj (г )
=о
4 2 У dxj dxt ^
О ij (r ) = Cjmn (r)S mn (r ) ,
(1)
на границе Г заданы перемещения
и, (0| Г = в,]г], (2)
где а у (г) - структурные напряжения, в у (г) - структурные деформации, Сутп (г )-
тензор модулей упругости неоднородной среды, ву - компоненты произвольного
заданного тензора малых макродеформаций.
Поля структурных напряжений а у (г) и деформаций в у (г) являются
случайными однородными всюду, за исключением области малой окрестности, прилегающей к границе Г. Для их описания используется статистический подход, связанный с нахождением моментов
^■3 =<а у(г)а у(Г) >, мй =<в \} (г)в(г) >.
Поле структурных модулей упругости Сутп (г) является статистически однородным и описывается следующей зависимостью:
СЦ„т (г ) = ЩЩ,,,,, + (1 - X(Г))Ст .
Систему уравнений (1) преобразуем к виду
с д 2ит' (г) =_ _
< 1}тП > дХп дХ} 11 ’Г
п , = д
((< (тп (Г) > + С'утп (Г))етп + С'утп (Г)и'т,п )•
11 ’1 дХ;
Полученное неоднородное дифференциальное уравнение будем решать методом функций Грина. При постоянном тензоре < С1;тп > функция Грина Отк(г,Г1) или
тензор Кельвина-Сомильяны вместе со своими производными обращается на бесконечности в нуль и удовлетворяет дифференциальным уравнениям
д 2 Стк (г, Г)
< с»» =_81к8(г _^
где 5 (г _ г1) - дельта-функция или функция Дирака.
Если < С11тп > - изотропный тензор, то
^ , 5 тк „ (гт _ Г1т )(гк _ Г1к )
Сутп (г , Г1) = АТ^^\ + В-----і----Ї3-----,
Г _гі\ \Г _ А\
где коэффициенты А и В связаны с постоянными Ламе тензора < С11тп >
соотношениями:
<Х > + < Зи, > < X > + < и, >
А =------------------ --------, В = •
8п < X >< (X + 2д) > 8п < X >< (X + 2д) >
Г
Возвращаясь к решению краевой задачи, имеем для неизвестного поля и, (г) интегро-дифференциальное уравнение
ди,' (г) гдОу (г, г ') д и \ \
=1~к—аХ7 ((< ( ( ( >+с]пк1 (г'))« + С) (г ')и'кЛ (г ')).
Полученное уравнение будем решать методом последовательных приближений. Приведем решение в первом корреляционном приближении, поскольку в дальнейшем ограничимся рассмотрением только этого приближения:
зи,'(г) гас,,(г,г•) д , г > _
= !—дХ-------(С((г )М' •
&а V дха дхп
Для пульсаций структурных напряжений и деформаций справедливы следующие соотношения:
^ га (Г) = - 21 {Ог] а (Г, Г' ) + ^ау,г (Г, Г')) ^ (СШ (Г ')еш )^' ,
2 V дхп
° ] (Г) Сутп (г )етп < Сг’тп (Г)б 'тп (Г) > +Сутп (Г)б 'тп (Г) •
Видно, что пульсации напряжений и деформаций являются случайными однородными полями, определяемыми геометрией элементов структуры, их упругими свойствами и характером взаимного расположения. На основе этого переходим к решению задачи в моментах:
<8у’ (Г*8ав' (Г) >= 4 е„пе - (С«тп - Ст<т„ )(С^„ - ) X (3)
х //(0,ку (Г,Г'') + 01К, (Г,Г "))(Са,(Г , Г ''') + Л (?, Г"’)) < А’(Г”)А'(?'") >,е dV,
<а(Паав (г) >= е е < С' (г)С' (г) > +е < С >< С' (г)в (г) > +
у\ / ар кі уд і}кі аруд кі аруд. і}Ш уд
+ е < С >< С' (г ) 8 ' (Г ) > + < Сік >< С >< 8 ’ (г ) 8 ' (г) > +
уд і]кІ аруд кі ікі аруд уд кі
+ е < С' (г)С' (Т)8’ (г) >+е < С’ (г)С’ (г)8 (г) > +
кі і}кІ ’ аруд ’ уд4' уд аруд 4 ’ укП ’ кІу ’
+ < С;кі >< (г)81д(г)8 к:(г) > + < Сар^д >< С’и (г)8 Ід(г)8 кі(г) > +
+ < С’кі(г )8 кі(г) >< СарТд(г )8 ід(г) > + < СаРтд(Г )С'кі(г )8 ід(г )8 кі(г) >• (4)
В качестве модели структурно-неоднородной среды рассматривается материал, который представляет собой однородную матрицу , где случайным, но равномерным образом расположены включения равного радиуса.
Для построения такой разреженной матричной структуры используется алгоритм, известный в литературе как метод Монте-Карло [4]. Суть его заключается в следующем.
В некоторую ограниченную область поочередно “выбрасываются” точки (центры сфер). Начиная со второй точки, производится вычисление расстояния между центрами сфер, которое не должно превышать значения 2d ^ - радиус сферы). В случае, когда данное условие не выполняется, точка “уничтожается” и “выбрасывание” повторяется. Процесс продолжается до тех пор, пока не будет
достигнута заданная объемная доля. При заданных объемной доле /, и числе
^ А
включений N вычисляется радиус сферы d = а
4
В противном случае, когда сделано предельное значение выбросов, но необходимое число точек не смогло разместиться в заданном объеме, вычисляется объемная доля, соответствующая вместившемуся количеству включений. Данный метод позволяет синтезировать разреженные структуры данного типа с объемным содержанием включений до 0,3.
Количественными характеристиками свойств синтезированной случайной структуры материала, необходимыми для решения статистических краевых задач механики структурно-неоднородных сред, являются моментные функции случайной индикаторной функции А(г) . Однако определение их вида для конкретных материалов требует ввода гипотез, нуждающихся в достаточном обосновании. Проблема может быть решена, если воспользоваться методами структурного моделирования на ЭВМ.
Вычислим моментную функцию второго порядка для синтезированной ранее случайной структуры:
К А (Г, Г) =< А '(Г) А '(Г) >=< А (Г) А (Г) > - < А ( г ) >2, (5)
где А (Г) есть индикаторная функция случайного поля, описывающая геометрию двухфазной структуры,
”(г) = <
1, г є сфе ре;
0, г є матрице,
здесь Г - радиус-вектор, А '(Г) - пульсация индикаторной функции, имеющая вид А '(Г) = А(г ) - (а (г ^, (..) - оператор статического осреднения.
При построении корреляционной функции в синтезированном фрагменте случайной структуры выделяется куб, который разбивается сеткой с шагом И = у 2 d . В узлах сетки проверяется наличие включения или матрицы и параметру А(Г) присваивается значение 1 или 0 соответственно. Для всех пар узловых точек,
г - г1
вычисляются произведения
отстоящих друг от друга на расстоянии А(Г)А(Г4) и определяется их среднее значение суммированием всех произведений А(Г)А(Г1) и делением их на количество рассматриваемых пар узловых точек. Для нулевого расстояния между точками находится значение (А (Г)) . Затем по формуле (5) вычисляется моментная функция к А (Г, Г), которая при Г = Г дает дисперсию величины А'(Г) : ВАк = ^(1-^) , где р - объемное содержание наполнителя.
В силу того, что А(Г) - статистически однородная и изотропная функция, К А зависит только от расстояний между рассматриваемыми узлами Г,Г4,Г2,Г3,... и инвариантна к их взаимному расположению, т. е. искомую моментную функцию п-го порядка можно представить в виде
К
(г , гі , г2,...)= я” Д (\г - гі ^ |г - |гі - г2 |,...) , (6)
где /” - нормированная корреляционная функция п-го порядка, - центральный
момент п-го порядка.
Рис. 1. Нормированная корреляционная функция 2-го порядка для случайных разреженных структур со сферическими включениями равного радиуса (г=0,39) с объемной долей р = 0,28; 0,24 при N = 70 (число разбиений ребра куба)
Рис. 2. Нормированная корреляционная функция 2-го порядка для случайных разреженных структур со сферическими включениями равного радиуса (г=0,39) с объемной долей р = 0,28 при N = 25, 35, 70
На представленных выше графиках корреляционная функция быстро убывает при увеличении аргумента. При переходе к меньшей объемной доле кривые становятся более пологими и сдвигаются в сторону увеличения аргумента (рис.1).Увеличение числа разбиений ребра куба N приводит к значительному затуханию нормированной
корреляционной функции относительно оси абсцисс, что подтверждает достоверность полученных результатов (рис.2).
Статистический подход механики структурно-неоднородных сред предполагает вычисление безусловных моментов структурных напряжений и деформаций, что непосредственно связано с поиском явного вида нормированных моментных функций /% и /I случайного параметра %(г). Для этого построенная ранее на ЭВМ первая
производная в нулевой точке [3]:
I- _ 12 Г I- _ 12 I _ |2
г - г
Л \г - А\ = exp(-с1_“Г
-)
с2 cos(c3 -
г - г
г - г
-) + с4 sin(c5 —
)
(7)
где ё - радиус включений. Показано, что константы с1 - с4 зависят от объемной доли включений и их численные значения соответственно равны (таблица):
константа р=0,28 р=0,24
с 1 1,03377243 1,04283448
с 2 0,9904881 0,993981383
с 3 1,8743548 1,859772621
с 4 0,7296027 0,720358433
Рис.3. Аппроксимированная нормированная корреляционная функция 2-го порядка для случайной разреженной структуры со сферическими включениями равного радиуса
(г = 0,39) с объемной долей р = 0,28
В соответствии с [2,3,7] может быть записано аналитическое выражение для моментной функции 3-го порядка случайного параметра %(г ) :
К %(0, г, Г) = £%3Л3(0, г, г-),
£>%= р(1 - р)(1 - 2 р),
где /% (0,г , г1)- нормированная корреляционная функция 3-го порядка, -
центральный момент 3-го порядка параметра Х(г ) ,
Л (0, г, г1) = ехР(-
с3
с1 2й 2
+ г11 + г - Г- )) х
с5
c2cos(—-(г + Г1 + Г - Г1 )) + c4sin(—— (г + ГМ + Г - гМ ))
(8)
С помощью представленных формулами (6), (7), (8) функций определяются моментные функции 3-го и 4-го порядков, необходимые в дальнейшем для определения моментов 2-го порядка полей структурных деформаций и напряжений:
К
;(0,0,г ) = (1- 2р) кх(0, г),
К % (°Д Г, г ) = (1 - 2р) к % (а Г, г ) + р(1 - р)К % (г, г).
Таким образом, в данной работе была рассмотрена задача определения моментов полей структурных напряжений и деформаций для синтезированных случайных структур со сферическими включениями равного радиуса. Дальнейшее вычисление моментов (3), (4) производится с помощью численного интегрирования по области статистической зависимости.
Библиографический список
1. Богачев И.Н., Вайнштейн А. А., Волков С. Д. Введение в статистическое металловедение. - М.: Металлургия, 1972. - 216 с.
2. Вильдеман В.Э., Соколкин Ю.В., Ташкинов А.А. Механика неупругого деформирования и разрушения композиционных материалов / Под ред. Ю.В. Соколкина. - М.: Наука, 1997. -288 с.
3. Волков С. Д., Ставров В.П. Статистическая механика композиционных материалов. -Минск: Изд-во БГУ им. В. И. Ленина, 1977. - 206 с.
4. Гаришин О.К. Геометрический синтез и исследование случайных структур // Структурные механизмы формирования механических свойств зернистых полимерных композитов. -Екатеринбург, 1997. - С. 48-81.
5. Кристенсен Р. Введение в механику композитов. - М.: Мир, 1984. - 336 с.
6. Победря Б..Е. Механика композиционных материалов. -М.: Изд-во Моск. ун-та, 1984. - 336 с.
7. Структурно-феноменологический подход к оценке прочности анизотропных композитных конструкций / Танкеева М.Г., Ташкинов А.А., Соколкин Ю.В., Постных А.М. // ДАН СССР. - Свердловск, 1989. - С. 3-24.
8. Шермергор Т.Д. Теория упругости микронеоднородных материалов. - М.: Наука, 1977. - 400 с.
Получено 15. 04.99.
