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

Деформирование и разрушение композитов

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

АКАДЕМИЯ НАУК СССР

УРАЛЬСКИЙ НАУЧНЫЙ ЦЕНТР

ДЕФОРМИРОВАНИЕ И РАЗРУШЕНИЕ КОМПОЗИТОВ

/

СВЕРДЛОВСК 1985

УДК 539.2 : 539.3 : 678.0

Деформирование и разрушение композитов: [Сб. статей]. Сверд­ ловск: УНЦ АН СССР, 1985.

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

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

О т в е т с т в е н н ы е р е д а к т о р ы докт. физ.-мат. наук Ю. В. Соколкин, канд. физ.-мат. наук А. А. Ташкинов, канд. физ-мат. наук В. Н. Аптуков.

Р е ц е н з е н т докт. техн. наук Г. К. Ибраев

„20304—1330—101 83

_

© УНЦ АН СССР, 1985

Д

------------------------- 5— - 11—1985

 

055(02)7 .

 

 

АКАДЕМИЯ НАУК СССР

УРАЛЬСКИЙ НАУЧНЫЙ ЦЕНТР

ДЕФОРМИРОВАНИЕ И РАЗРУШЕНИЕ КОМПОЗИТОВ

1985

Ю. В. СОКОЛКИН, А. А. ТАШКИНОВ

ПРИНЦИП ЛОКАЛЬНОСТИ И МЕТОДЫ РЕШЕНИЯ СТАТИЧЕСКИХ ЗАДАЧ МЕХАНИКИ СТРУКТУРНО-НЕОДНОРОДНЫХ ТЕЛ

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

Постановка задачи.

Рассмотрим статическую задачу о дефор­

мировании

области

V

с

границей Г. Будем считать, что мате­

риал

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

собой композит матричного типа,

причем

характерный размер

неоднородностей

(волоков частиц

наполнителя и т. д.)

намного меньше характерного/^азмера об­

ласти.

Предположим,

что

форма

и свойства

однородных компо­

нентов структуры детерминированы и заданы,

а/взаимное рас­

положение таково, что кусочно-постоянное

поле структурных

модулей упругости

С (г)

есть

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

 

 

 

 

 

< А imn{r) )

=COnst,

 

 

 

(1)

•где

( . . . )

— оператор

осреднения по

случайным

однородным

полям. Если случайное поле С (г) еще и эргодичное,

то его сред­

нее значене не зависит от реализации.

 

 

 

 

Пусть в некоторой

декартовой

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

щие

соотношения,

связывающие

тензор

деформаций е и тензор

напряжений

а, задаются в виде

 

 

 

 

 

 

 

 

 

<*ii(r) = Ciimn{r)Emn(r).

 

 

(2)

.Деформации будем считать малыми, так что выполняются соот­ ношения Коши

М г) = 7 а [“ t , j ( r) + u i,t (г)1-

(3)

*

Пусть также заданы уравнения равновесия среды

<г//./(г) = 0-

(4>

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

o\pnj= aif)n}.

(5>

Тогда стохастической краевой задачей I теории упругости структурно-неоднородных сред будем называть краевую задачу, описываемую замкнутой системой уравнений (2) — (4), условиями контакта (5) и краевыми условиями Дирихле

Щ= еиг},

(6>

где е — произвольный симметричный тензор второго

ранга,

имеющий смысл тензора макродеформаций.Действительно

при

принятых условиях справедливо

 

<е(г)>5ЕЕе*=е

(7)

(звездочка вверху означает макроскопические физические вели­ чины).

Стохастическую краевую задачу II можно получить заменой в задаче I краевых условий (6) краевыми условиями Неймана

T l {r)=pijnj {r) |г ,

(8)

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

<з(г)> =<т*=р.

(9>

Принцип локальности. «В расположении и взаимодействии струк­ турных неоднородностей имеет место ближний порядок».

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

полей

деформирования в

произвольном включении и в матрице

вокруг

него — следствие

ближнего порядка во взаимодействии.

Метод периодических

составляющих. Локальные моментные

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

дических

составляющих

основан

на разложении поля

С (г),

а также искомых полей

и (г)

е (г),

а (г) в стохастических задачах

I и II на

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

составляющие, известные из

реше­

4

ния периодической

задачи для области Vp с регулярной струк­

турой, и соответствующие случайные отклонения.

 

Рассмотрим детерминированную задачу для области Vp, отли­

чающуюся от аналогичной ей стохастической задачи

I только тем,

что упругие свойства

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

периодические

кусочно-непрерывные функции. Предположим, что

решение этой

задачи, которое

может

быть получено, например,

методом Бах­

валова— Победри

[4], известно:

 

 

 

 

efi(r)=F?lmn(r)em„;

(10)

индексом

р вверху

обозначены физические величины детермини­

рованной

задачи для

Vp,

(е р ) =е. Под оператором осреднения

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

Условие, с помощью которого случайной структуре

области V

ставится в соответствие регулярная структура

области

Vp, выбе­

рем в виде

 

 

<С1/тЛ г)) = <С?1тп(г)).

 

(11)

Отметим, что условие (11) не является, вообще говоря, одно­ значным.

Используя разложение случайных полей задачи I

 

C (r)= 0 (r)+ C (r),

(12)

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

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

полей деформирования:

 

 

[Cfjm n {?) Ит,п 7)] ,} =

[Сц'тп (r ) етп (7 + С ijmn (r ) и т,п (f)] ,/>

e«j (г)= 7* Щ,/ (г)+ И/,*(г)],

(13)

Oii{r) = Ctjmn {г)^тч {г)+С1Ып (г) б£п(г), и |г =0.

Решение задачи (13) будем строить с помощью функции Грина для неограниченной среды с регулярной структурой. Тогда всюду, за исключением малой окрестности, прилегающей к гра­

нице Г,

для

составляющей

и (г) случайного

поля и (г) спра­

ведливо

 

 

 

 

 

 

= J

 

Pi) \Cqkmn (pi) етп{ri)~\~Cqkmn (^l)

,n (^1)] ,k drlt

(14)

V

 

 

 

 

 

 

где функция

Грина

Gp(r, rx) обращается вместе со своими

про­

изводными в ноль на

бесконечности и удовлетворяет уравнению

 

[ей™, (г) с и ( г ,

г.)]. ) = - 6 (г - г.)

 

(15)

(6 r j

— функция Дирака, р — символ Кронекера).

 

5

Решая интегро-дифференциальное уравнение методами после­ довательных приближений или малого параметра [6], получаем

«(,/(/•)= 2

1М№п(г)еа„

(16)

Мцтп(г) —V Giq.j

» rl)

ifl) Fstmn if1)] ,k drlt

(17)

M lfU r )- J G?q,j(r,

r,)

 

(/"1)].лdrx.

(18)

В итоге, например, для поля микродеформаций будем иметь

е„ (г)= \р!1ап(г)+ Д M jjL {г)

(19)

К преимуществам метода периодических составляющих сле­ дует отнести: а) установление связи между решениями задач для тел с регулярной и случайной структурой; б) новое разложение полей деформирования в виде (11) позволяет уже в первом при­ ближении решения стохастической задачи учесть все упомянутые выше тонкие эффекты микромеханики.

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

 

V - а (г)=0,

 

 

 

Ink [DP (г).. <у(г)]= - Ink [D (г)., а (г)+ & (г)..аР (г)],

(20)

 

ст(г)|г =0,

 

 

 

где CijmnDmnpq = I Upq,

/ — единичный тензор

4-го ранга,

 

Ink (в) = ЕЕ jftjE jmn^

rfihm'

(21)

Метод локального

приближения.

Идея

метода заключается

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

А

соседних включений и однородного поля напряжений aiJt задан­

ного на

достаточном удалении от ансамбля включений. В силу

наличия

ближнего порядка во взаимодействии неоднородностей,

 

А

можно подобрать такие значения комнонент тензора <J , при кото­ рых поля деформирования в центральной области ансамбля включений совпадут с искомыми напряжениями и деформациями стохастической задачи II для области V. Естественно, что реше­

ние стохастической

краевой

задачи таким путем будет получено

в реализациях. В

рамках

теории упругости структурно-неодно­

6

родных тел

связь между заданными макронапряжениями Ъ*—р

и тензором

Л

 

 

 

 

а линейная:

 

 

где А — симметричный

j = АцтпРтп*

 

(22)

тензор 4-го ранга

(Aijmn= A jimn= A ijnm).

Таким образом, метод позволяет свести стохастическую (детер­

минированную)

задачу для области V (Vp) к

краевой задаче для

области, содержащей малое число включений.

 

В

работах

[3, 6]

содержатся обоснование и примеры реали­

зации

метода

локального приближения

для

тел с регулярной

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

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

Л

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

 

 

Л И Т Е Р А Т У Р А

1.

В е н т ц е л ь Е.

С. Теория вероятностей. М.: Физ-матлетгиз, 1962.

415 с.

В о л к о в С. Д., С т а в р о в В. П. Статистическая механика композит­

2.

ных материалов. Минск: Изд-во БГУ, 1978. 206 с.

3.

И в а н о в В. Н.,

Т а ш к и н о в А. А. Расчет полей структурных на­

пряжений в микронеоднородных упругих средах с регулярной структурой.—

В сб.: Структурные превращения в

полимерных и жидких кристаллах.

Свердловск: УНЦ АН СССР, 1981, с. 120—124.

4.

П о б е д р я Б. Е. Численные

методы теории упругости и пластич­

ности. М.: Изд-во МГУ, 1981, с. 269—280.

5.

Со ко л кин Ю. В., И в а н о в

В. Н., Т а ш к и н о в А. А. Структур­

ный анализ упругопластического деформирования однонаправленных во­ локнистых композитов.—Тезисы докладов 4-й Всесоюзной конференции по композиционным материалам. М.: Изд-во МГУ, 1981, ч. 2, с. 174.

6. К решению стохастических краевых задач теории упругости струк­ турно неоднородных тел А. А. Ташкинов; Ин-т toex. сплошн. сред УНЦ АН

СССР. Пермь. 1981, 20 с., библиогр. 9 назв. Рукопись деп. в ВИНИТИ 26.08.81, № 4329-81. Деп.

АКАДЕМИЯ НАУК СССР

УРАЛЬСКИЙ НАУЧНЫЙ

ЦЕНТР

ДЕФОРМИРОВАНИЕ И РАЗРУШЕНИЕ КОМПОЗИТОВ

1985

С. Е. ЕВЛАМПИЕВА

НАПРЯЖЕННОЕ СОСТОЯНИЕ УПРУГОЙ МАТРИЦЫ ПРИ ХАОТИЧНОМ И РЕГУЛЯРНОМ ЗАПОЛНЕНИИ ОБЪЕМА

КОМПОЗИТА КРУГЛЫМИ ЖЕСТКИМИ ВКЛЮЧЕНИЯМИ

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

Нами было показано \ что напряжения в матрице, содер­ жащей регулярно расположенные жесткие включения, распре­ делены неравномерно. Актуальным является изучение характе­ ра этой неоднородности. Так как в реальных системах части­ цы наполнителя расположены нерегулярно, необходимо выяс­ нить, в какой мере нарушение регулярности (при том же объ­ емном наполнении) изменяет картину НДС, присущую регу­ лярным системам.

На рис. 1 показаны две системы из семи близких друг дру­ гу регулярных и хаотично расположенных жестких включений. Указанная группа включений размещена в неограниченной эла­ стомерной матрице и находится в состоянии плоской дефор­ мации. На удалении от включений системы нагружены единич­ ным растягивающим напряжением в направлении оси х. Объ­ емная доля включений в области их расположения составляет 50%.

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

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

1 Е в л а м п и е в а

С. Е. Напряженное состояние упругой

матрицы при

регулярном заполнении

объема композита круглыми жесткими

включениями

(плоская задача).— В сб.: Структурная механика неоднородных сред. Сверд­ ловск: УНЦ АН СССРГ1982, с. 69—71.

2 М у с х е л и ш в и л и Н. И. Некоторые основные задачи математической теории упругости. 5-е изд., неправ, и доп. М.: Наука, 1966. 707 с.

8

Рис. 1. Схема нагруже­ ния регулярных и хао­ тично расположенных жестких включений.

Пояснения см. в тексте.

-© Q

©т*1© з © I;

е > е -- е е

Рис. 2. Распределение полей гидростатических (а) и интенсивностных (б) напряжений вокруг включений.

Рис. 3. Функция плотно­ сти распределения ох.

Пояснения см. в тексте.

нагрузки и иллюстрирующих внутреннее сопротивление внеш­ ней (со стороны матрицы) нагрузке.

Из распределения поля средних гидростатических напряже­

ний для

регулярно и хаотично расположенных включений

(рис. 2а)

видно, что максимальные напряжения сосредоточены

впромежутках между включениями в направлении действия внешней нагрузки. Количественная разница существенна: для регулярных систем максимальные напряжения равны 3,35, для хаотично расположенных 7,6. Вероятность разрушения от дей­ ствия гидростатической составляющей напряженного состояния

вслучае нерегулярности выражена сильнее.

Изолинии главных деформаций представлены на рис. 26. Для регулярной системы включений максимальные их значе­ ния равны 2,02, для хаотично расположенных 3,9. Разница в деформациях почти двухкратная.

На рис. 4 представлены поля напряжений а*. Интенсивность напряженности характеризуется толщиной темных линий. И в этом случае очевидно существенное усиление неоднородности по на­ пряжениям ох, а именно, они равны 8,6 в нерегулярных систе­ мах и 3,5 — в регулярных. График функции плотности распре­ деления Ох для обоих вариантов представлен на рис. 3. Штри­ ховая линия относится к регулярным структурам, сплошная — к нерегулярным. Из рис. 4 можно сделать вывод, что нерегу­ лярность приводит к резкому усилению напряженности в не­ значительном объеме матрицы, а также к существенному увеличению объема матрицы, который характеризуется низки­ ми и отрицательными значениями ох. В формировании внут­ реннего противодействия активно участвует около 40 % матри­ цы, 60 % вообще не нагружены. Напряженные участки матри­ цы в совокупности с частицами образуют сетеобразную внутренюю структуру.

10