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

УДК 539.3

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

ЕЛО. Макарова

Abstract

Using the boundary value micromechanical problem for heterogeneous solid new equations between deformation in stochastic composites and deformation in composites with periodical structure are achieved.

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

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

при макроскопически однородном деформированном состоянии. В формулах (1), (2),

(3) через а, е, а и 0 обозначены структурные поля напряжений, деформаций.

перемещений и модулей упругости соответственно, е - заданное поле макроскопических деформаций.

основы решения аналогичной краевой задачи для среды с регулярной микроструктурой:

V  ст - 0, є = def и, ст ” 8 • в ,

(О

с граничными условиями

v

которым, как известно [1], эквивалентны условия на поверхности Б тела V

Идея излагаемого ниже метода заключается в использовании в качестве

V-ст'1” = о, 8(р) = def u(p) , а<р) = С<р) -є(|,),

[8

где и(р), Б(р|, ст(р- - детерминированные периодические функции структурных

перемещений, деформаций и напряжений , С(р) -тензор структурных модулей упругости среды с регулярной структурой. Предположим, что решение краевой задачи

(4) нам извес тно [ 3 ] :

Е1р),(г) = ТЧ(р)(г)--е*< е(р)=е* +Е(Р)’,

с*(р> /с(р)(г)\ + ^с(р)(г)„м(р)(Г)\ 5 ст(р) =с*{р) -8*,

где С*(р) - эффективные модули упругости среды с регулярной структурой, Ы|р)(г)-структурные функции [3], () - оператор осреднения по представительному объему.

Ниже устанавливается важное свойсгво упругих микронеоднородных деформируемых сред. Если микронеоднородная среда макроскопически однородна и квази изотропна, перемещения границы 5 тела V, имеющего конечные размеры, детерминированы, дисперсии физических свойств среды конечны, структурные деформации - микроскопически гладкие функции координат, то существует случайный функционал Ф^р)(0), зависящий от граничных условий, такой, что пульсации

о

структурных деформаций е(г) связаны со структурными деформациями в регулярной среде 8<р) соотношением:

е

Для доказательства соотношения (5) рассмотрим решение краевой задачи (6), которая приводится к интегрированию системы дифференциальных уравнений при нулевых граничных условиях:

° с

\/.

о о о /о о\

г де Г1 - 0- -1^' и'Р}+ 0- с!е1' и- | 0- с1еГ и }.

/

Уравнения (6) можно рассматривать как уравнения краевой задачи теории упругости микронеоднородных сред с регулярной структурой С<р)(г) и

перемещениями и(г), обусловленными действием фиктивных случайных объемных сил V 11.

Вводя функцию Грина среды с регулярной структурой С(р)(г,г'), система

уравнений (6) преобразуется в систему интегро-дифференциальных уравнений:

и

V

Для определения полей структурных деформаций необходимо знать градиент пульсаций структурных перемещений, поэтому дифференцируем (7):

Уи

V

Уравнение (7) решаем методом последовательных приближений при

ограничениях, сформулированных в виде макроскопической однородности и квазиизотропности микронеоднородной среды и микроскопической гладкости (в в /-окрестности микронеоднородной среды) структурных деформаций в регулярной структуре.

В первом приближении полагаем:

°^1'1 г < \ 0

V и

V

Для макроскопически однородной среды интегралы в (9) фактически

распространяются на е2/-окрестность микронеоднородной среды, где е(р>-постоянны, поэтому соотношения (9) принимают вид

0П)

Уи

г

где Ур(р)(1) = |УС(р)(г>Г')-(У'-0 )с1У' ,

V

' (1) ^

а р|р> (в,г)- тензор-функционал третьего ранга относительно физических свойств среды.

Подставляя (9) в (8), с учетом (10) получаем второе приближение:

о<2>

У и ^(Ур<р,“' + Ур(р)(2')--8(р) ,

г

Ур(р)(2) |уо(р) • (V'  (0 • У'р(р)'кЧ) )ё\".

V

Окончательно получаем:

Уи(г)= Ур(р)--е(р),

Ур

где К Ь

к=1

Поскольку пульсации структурных деформаций определяются выражением

о о

е(г) = ёеГ и(г), то в силу (11) приходим к соотношению (5):

8^ф

где функционал Ф(р)(0) определяется соотношением:

Ф

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

Библиографический список

1 Соколкин Ю.В., Ташкинов А.А. Механика деформирования и разрушения структурно-неоднородных тел. -М.: Наука, 1984. -116с.

2 Соколкин Ю.В., Вотинов AM., Ташкинов А.А., Постных А.М., Чекалкин А.А Технология и проектирование углерод-углеродных композитов и конструкций - М: Наука, Физматлит, 1996. -239 с.

3 Победря Б Е Механика композиционных материалов. М.: МГУ, 1984, -336 с.