Моделирование несжимаемых слоистых композитов с конечными деформациями на основе метода асимптотического осреднения

УДК 539.3

БОТ 10.18698/2309-3684-2017-1-3254

Моделирование несжимаемых слоистых композитов с конечными деформациями на основе метода асимптотического осреднения

© Ю.И. Димитриенко, Е.А. Губарева, Д.Ю. Кольжанова, С.Б. Каримов

МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

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

Введение. Значительный интерес для промышленности представляют композиционные материалы, способные выдерживать большие деформации (до 100 % и более) и состоящие из резиноподобных и эла-стомерных матриц и высокопрочных армирующих наполнителей (тканей, волокон и др.). Для расчета прочности таких композитов с конечными деформациями в зависимости от их микроструктуры необходимы методики, позволяющие связать эффективные механические характеристики — диаграммы деформирования в области больших деформаций с нелинейно-упругими характеристиками отдельных фаз композита. Для малых деформаций существует достаточно много методов решения этой задачи [1-7], однако для композитов с конечными деформациями большинство из них неприменимы в силу существенной нелинейности задачи.

Наиболее перспективен для расчета эффективных диаграмм деформирования метод асимптотического осреднения [1-3]. В работах [8, 9] этот метод был применен для расчета нелинейно-упругих свойств слоистых композитов с большими деформациями. В настоящей работе предложено дальнейшее развитие этого метода для расчета диаграмм деформирования несжимаемых слоистых композитов с конечными деформациями. Несжимаемыми средами с конечными деформациями являются практически все резиноподобные материалы [10-15]. Известно, что несжимаемые материалы некорректно рассматривать как предельный случай слабосжимаемых материалов, для них используется специальная постановка задач теории упругости с дополнительным условием несжимаемости, требующая особых методов решения.

Для определяющих соотношений фаз композита были использованы так называемые универсальные представления моделей несжимаемых нелинейно-упругих сред с конечными деформациями [9], предложенные Ю.И. Димитриенко [15, 16]. Они позволяют в рамках единого методического подхода решать задачи нелинейной теории упругости сразу для нескольких классов моделей сред с конечными деформациями.

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

к 0

точек в отсчетной и актуальной конфигурациях как х и Хк, лагран-

жевы координаты как X1 (их будем полагать совпадающими с де-о .

картовыми X1 = х1). Рассмотрим слоистый композит, в отсчетной

0

конфигурации представляющий собой область V, являющуюся системой параллельных периодически повторяющихся N слоев, ортогональных к направлению ОХ . Тогда можно ввести ячейку периодичности (ЯП) — набор из N слоев суммарной толщиной I. Введем малый параметр к = I / Ь«1, как отношение толщины ЯП I

о

к общей толщине композита Ь в отсчетной конфигурации К, а так-

о

же локальную лагранжеву координату 2, в К, которая связана с X1 соотношениями

X3 — X' Ъ= —, X' = ^.

к Ь

Поместим нулевое положение локальной лагранжевой координаты в середину ЯП, тогда £ будет изменяться в диапазоне -0,5 < £ < 0,5.

Границы раздела слоев в ЯП обозначим как £ = £а, а = 1, ..., N -1.

Рассматривая слоистый композит как неоднородную среду, сформулируем для него задачу нелинейной теории упругости в ла-гранжевом описании в общей формулировке [14] с использованием универсальных моделей сред с конечными деформациями, предложенных Ю.И. Димитриенко, для несжимаемых упругих сред

0 .. 0 . . о

ViP1J +р fJ = 0, X1 eV;

k k 0 k .oo F\ =5f +Viuk, X1 eV<^Z;

det

0 .. . . 0

m [P ] = 0, [u1 ] = 0, X1 eZap;

0 . . 0 0

mP. = t., X1 eZ1, u1 = u., X1 eZ2.

Здесь: (2) — уравнения равновесия; (3) — определяющие соотношения нелинейно-упругой среды; (4) — кинематическое соотношение; (5) — условие несжимаемости; (6) — условия идеального контакта

0

на поверхностях раздела Zap а и Р-й компоненты композита; (7) —

00

граничные условия на частях Z1 и Z2 внешней поверхности компо-0 0 0 0 зита ( Z1 ^Z2 = 5 V ). Скачок функций на границе раздела Zap компонент композита обозначен [PlJ ]. Все компоненты векторов и тензоров отнесены к неподвижному ортонормированному базису ek —

0

отсчетной конфигурации K.

0

Кроме того, введены следующие обозначения: р — плотность, p — гидростатическое давление, а также компоненты тензора напряжений Пиолы — Кирхгофа P1J, тензора градиента деформации Fk, обратного градиента деформации F-1J, вектора перемещений

k00 u , вектора нормали m к поверхности в K, вектора поверхностных

усилий / вектора заданных перемещений поверхности игг, вектора

плотности массовых сил набла-оператор V; =

а

ах1

В соотношении (3) — компоненты тензора, описывающего

нелинейно-упругие свойства фаз композита согласно моделям Вп упругих несжимаемых сред с конечными деформациями [15, 16]

г ((, *)=р (к ) ^

о

Л

рг

V

sq

орг, £

а

а ощ

V

((п) \

орг, £

)

( п) Е01М

/ , \ 3 л 0 . .00 (п) 1

(к )= х < бу б /е е, орг = 1

у ,ю=1

3 0 0

—3 ^ п"3 ере;-

п 3 у=1

10у.Щ

Здесь Е Jq — компоненты тензора энергетической эквивалент-

((п) Л

ности [16]; V

о

рг

упругий потенциал, различный для каждой

фазы композита, и поэтому зависящий явно от О рг — компоненты симметричного тензора энергетических мер деформации [16]; ¿УсД^у) — функции собственных значений тензоров искажений,

собственные значения являются также функциями только от ком-

0

понент тензора // [16]; еУ (к ) и е» ) — матрицы собственных

векторов левого и правого тензора искажений, являющиеся неявными функциями от .

Асимптотические разложения в нелинейной задаче механики для слоистого композита. Решение нелинейной задачи механики (1)-(7) находим относительно поля вектора перемещений ик = = ик (х1, £) и гидростатического давления р (X1, £) в виде квазипериодических функций координат, т. е. периодических функций по локальной лагранжевой координате и зависящих от глобальных

лагранжевых координат X1:

0

ик

[X1,£) = ик [X1,£ + а), р[X1,£) = р[[,£ + а),

где а — произвольное целое число.

Дифференцирование квазипериодических функций /к [ X1, осуществляем по формальным правилам дифференцирования сложной функции, тогда с учетом (9) имеем

0 1

VI ик = ик +-и,к 513,

,1 К I3

где ик -X,ик(X'.5), и3 = |ик[X',5).

Будем искать решение задачи (1)-(5) для композита периодической структуры в виде асимптотических разложений по параметру к :

ик [X1, £) = ик(0) [X1) + кик(1) [[, £) + к2...,

р [, ф р

Подобные асимптотические разложения использованы в задачах о тонких упругих многослойных пластинах [17-19]. Подставляя (11) в соотношения (4) с учетом правила дифференцирования (10), находим асимптотическое разложение для градиента деформации:

рк = рк (0) [X1, + (1) [[, £) + к2...,

РЩ) [X1, ^) = 5к + ик

' \ 'V ' ,' /3 '3, ' ,' /з '3 4 '

Подставляя разложение (11) в соотношения (3) и используя формулу Тейлора, находим асимптотические разложения определяющих соотношений и тензора Пиолы — Кирхгофа:

р1 = р1 (0) [xk,£) + кР1 (1) [[к,£) + к2Рт) [[к,£) +...,

р/

1 = _р(0) [[-1 )(1)1 _р(1) [[-1 )(0)1 [к(0), (1).

Подставляя разложение (12) в уравнение (5), с использованием формулы Тейлора находим асимптотическое разложение уравнения несжимаемости:

к (1) +к2... = 1.

После подстановки разложения (13) в уравнения равновесия (2) с учетом (10), получаем асимптотическое разложение уравнений равновесия и граничных условий:

1

-РЬ(0)5;3 +

(

к

13

рЬ

1(0) + Р/3(1)^3 + Р Г ] + к ((1) + р/3 (2)8Й ) + к2... = 0,

0

п

и

3 к

р3Ь (0)'

+ к

0

+ кп3

к

р3] (1)'

+ к2... = 0,

и

+ к2... = 0, ХЬ еЕар,

Пзру(0) +к^зр3'(1) +к2... = г{, ХЬ еЕь

ик (0) +кик (1) + к2... = и\, X1 еЕ2 .

Локальные задачи нелинейной теории упругости для несжимаемых сред. Приравнивая к нулю в уравнениях (13)—(15) члены при одинаковых степенях к, получаем рекуррентную последовательность Ьк локальных задач нелинейной теории упругости для несжимаемых сред. Задача Ь0 имеет следующий вид:

р3Ь(0) = 0,

рЬ

Рк(0) (, ^) = Рк + ик3(1)5/3, (((0)к ) = 1,

р

3Д0)"

= 0,

к

= 0, ^ = , а = 1,...,п_1,

ик (1Л = 0

к

= 0, (р(0) = р.

Здесь введена операция осреднения по ЯП V ^

0,5

< ик(1) >=

| ик(1)^,

-0,5

обозначены средний градиент деформации композита

Ёк =8? + ик <0>

и среднее гидростатическое давление в композите р.

Наличие условий (ик(1)^ = 0 и = Р вызвано требованием

единственности задачи на ячейке и делает ее интегро- дифференциальной. Условие ик (1) = 0 означает периодичность неизвестных функций на

границе ЯП

V(1)" = ик (1) _ ик (1) = °

_ - - _ з

Задачу Ь° рассмотрим относительно неизвестных перемещений

к

-г? —

при этом средний тензор-градиент г 1 и среднее давление р

являются «входными данными» локальной задачи Ь°. Локальная задача Ь1 имеет вид:

рШ + р/(1)8з + £ ^ =

р,(1) =_р(°)(-1 )(1)* -р(1)(г-1)1

.лШ' д

дрк (°) г? (1) = и} + икз(2) 8/3,

^ ** (р))(1) = o,

? °* (( (°), (1),

р3, (1)'

= °,

к (2)"

= °, £ = £а, а = 1,...,п-1,

ик (2Л =

и

к (2)'

= °, (Р(1)) = °.

Эта система уравнений рассматривается относительно перемещений ик (2) и давления р(1) и является линейной относительно данных функций.

Осредненная задача нелинейной теории упругости для композита. В силу периодичности функций на ЯП имеет место соотношение

= °.

Осредняя систему уравнений (15), в том числе уравнения (13), с учетом (20), получаем осредненную задачу нелинейной упругости для композита

<Pij(0) >f +<р fj >= 0,

< Pj(0) >=-< p(0) (f-1 )(0)j > + < F0j (Fk(0), é) < F/(0) >=5k + u\(0),

0 0 0 m < Pij(0) >=< tj >, Xi eSl, uk(0) =< uj >, Xi e Z2.

Покажем, что осредненный градиент деформации < Fk > совпа-

дает со средним градиентом деформаций ¥к :

< Fk >= Fk.

Действительно, используя третью формулу из системы (16) и осредняя ее по ЯП, в силу периодичности функций ик (1), получаем:

0,5

<

Fk(0) >=< Fk + u/k3(1)5/3 >= Fk + 5/3 J u^:

k

-0,5

= FK +5/3luk

^=0,5

- u

4=-0,5

=Fk.

Вычислим теперь детерминант градиента деформации ¥/(0)к. С учетом (12) и (16), получаем

det (((0) k) = det (Fk + uk3(1) 5/3) = det

Раскрывая детерминант матрицы по третьему столбцу, получаем линейную дифференциальную форму

F11 F1 F31 +«m» 1

det F12 F22 F32 + "S1'

F13 F3 F33 +•3!1> J

det (((0)k ) = det

+ det

F12 F221

F F23 J

F11 F1

F2 F2 J

(( + u/31))-det

(( + u/33(1)).

f F11 F21 1 VF13 F23J

( + u/23(1)) +

Осредним теперь соотношение (25) по ЯП и примем во внимание, что детерминанты матриц размером 2 в этом выражении от £ не за-

висят, а < и/3(1) >= ° в силу периодичности и] (1). Тогда

,

<

ёе!

Г^2

>= ёе!

г2 к

2

V К

к

< гз1 + и

1(1) /з

> _

2

_ ёе!

щ щ1 >

V г>з Щзу

< Кз2 + и

2

/з

> +

+ ёе!

_ ёе!

щ ц ^

А К?,

щ ЩА

V Кз Щз у

Г^2

< гзз + и?з(1) >= ёе!

г2 ц

2

Кз Кз

V11

Кз _

2 У

Кз2 + ёе!

щ

г2

К

1

ц

2 2

2 У

Цз = ё*Ц).

В силу того, что все компоненты композита являются несжимаемыми, ёе!

ёе! Ц)

= 1.

Иначе говоря, композит в целом — несжимаемый материал.

Из соотношения (27) вытекает, что в системе из N +1 уравнения

ёе!(((( ,£)) = 1, _£а<£<£а+1, ёе! (( (х* )) = 1

а = °, 1, ...,N_ 1,

независимыми являются только N уравнений.

Перемещение ик(1), являющееся решением задачи Ь° (16), может быть представлено как функция локальных координат и «входных данных задачи»: среднего градиента деформаций Цк и среднего гидростатического давления р. Тогда в виде такой же зависимости

Рт = рк

может быть представлен и градиент Ц^ (°) нулевого приближения.

Осредненные определяющие соотношения для композита с несжимаемыми фазами. Подставляя (29) во вторую группу соотношений (21), получаем осредненные определяющие соотношения композита, записанные в неявной форме:

< Р*(0) >= - < р(0) (¥-1 ) > + у ((к, р),

П /_ ч (") I- /_ ч \

(0) / -1 \(0)У

Вычислим компоненты тензора < р (¥ ) >. Поскольку

ёй (¥(0)у ) = 1, компоненты обратного тензора (¥-1)* , вычисленные через алгебраические дополнения элементов ¥(0)у, в явном виде могут быть записаны как

(¥-1 )(0,и = ¥2 (3+и33'>)-¥? (( + и2з<"),

(¥-1 )<0)22=¥ (+и33">)-¥3 ((+ и>™),

(-1 У0'33 = ¥.¥ - ¥',

(¥-1 )<0'13 = ¥> (( + и.™)-¥? (( + и)™), (¥),0,3, = ¥,2¥ - ¥¥22, (¥)(0)23 = ¥2 ( + и'<»)-¥' (( + и3("), (¥),0,32 = ¥^3 - ¥*¥}.

Часть этих компонент (¥ ) не зависят от координаты £

и совпадают с соответствующими компонентами обратного осреднен-

/—_1\У I -л(0)у

ного тензора (¥ ) . Остальные компоненты (¥ ) являются линейными функциями от частных производных ик(1). Тогда при осреднении компонент из первой группы, умноженных на р(0) ((к, £), получаем:

< Р

При осреднении компонент из второй группы, умноженных на р(°) ((к, £), получаем выражения

< р(°) (К_1 )(0)23 >=< р(°) (щ2 (Кз1 + и/з1))_ Ц (( + и/зз(1)) >=

=< p(0) >

Очевидно, что (р22Щ1 _ рЩ3 ) = (р 1) . Положим далее, что ЯП обладает симметрией относительно плоскости £ = °. Тогда в силу периодичности функции р(°), ик(1) являются симметричными по координате р(°)((,£) = р(°)(хк,_£), ик(1)(хк,£) = ик(1)(хк,_£). Производные и/к,(1) от симметричной функции являются антисимметричными: ик(1) ((к, £) = _и/3(1) (хк, _£). Антисимметричными будут и произведения функций р (°) ((к, £)икз(1) (хк, £) = _ р (°) (хк, _£)х

х ик(1) ((к, _£). Тогда интегралы по всей ЯП от антисимметричных функций обращаются в нуль:

< р

_°,5 °

Следовательно, в выражениях (33) все слагаемые, содержащие множители < р(°)ик(1) >, равны нулю. Исходя из этого, получаем

< Р

Аналогично преобразуем все остальные компоненты тензора (0)/ г7-Л(0)*

< Р IF ) > второй группы компонент:

<р

Объединяя формулы (32) и (36), получаем, что для всех компонент обратного тензора < p(0) (f-1J > имеет место соотношение

(36). Тогда осредненные определяющие соотношения (30) для композита можно записать:

P = -p

- (?) ,_ 4 (") /_ /_ 4 \

Pij =< Pj(0) >, F0j ((, p)=< F0j (Fk + uk(1) ((, p, £m )/3 5,3 ), £> .

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

F0j (F , p) определяющих соотношений композита. Для фаз композита имеет место более простое определяющее соотношение (3).

Тензор истинных напряжений Коши композита вычисляем по формуле Tij = Fimpmj\ Тогда из формулы (37) получаем

_ (-) /_ ч n /_ ч _ (n) /_ ч

Tij = -p5j +F0ij

т. е. представление определяющих соотношений композита в виде связи тензора напряжений Коши T1J, осредненного градиента деформации Fk и среднего давления p.

Решение локальной задачи нулевого приближения. Задача (16) является нелинейной, но одномерной, в ней все функции зависят только от что позволяет найти ее формальное решение.

Интегрируя уравнения равновесия в системе (16), получаем, что напряжения P3j (0) постоянны в ЯП

P3j

где CJ — постоянные интегрирования.

Из третьего уравнения системы (16) следует, что из девяти компонент градиента деформаций Ff(0) от координаты £ зависят только три, а остальные шесть совпадают с компонентами среднего градиента

Fk

Тогда, подставляя (37) во второе соотношение в системе (16) для напряжений Pij(0) вместе с условием несжимаемости (четвертое со-

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

и р(0):

4

с" =-р(0) (-1 )(0)3/ +^03/ (( (0), П, %),

ёе1

р(0Л = р, -£а<£<£а+1, а = 0, 1, ..., N-1.

Средний градиент в системе (41) считается известным и удовлетворяющим осредненному условию несжимаемости (27). Тогда с учетом (28) только для N -1 слоя система первых двух уравнений в системе (41) является независимой. Для всех N слоев эта система будет зависимой, поэтому условие несжимаемости для одного из

слоев в системе (41) заменена на условие нормировки = р.

Запишем формальное решение системы (41):

з(0) =

^(0) =^к (,рк, р, %),

Тогда, подставляя первую формулу системы (42) в выражение (40), получаем систему трех обыкновенных линейных дифференциальных

уравнений относительно перемещений ик

- (п) , _ ч

^ + ик

которую легко интегрируем:

%

ик(1) =

-0,5

Здесь Вк — константы интегрирования.

Условия идеального контакта в системе

р3, (0)"

= 0,

ик

= 0 выполняются автоматически. Подставляя выражение

в условие нормировки в системе (16) (и?(1)) = °, находим константы

интегрирования:

£ (п) , _ ч Вк =< | £к ((, Цк, р, £)>.

_°,5

После подстановки выражения (44) в условие периодичности системы (16), получаем

<

>= Кзк,

представляющее собой систему трех нелинейных алгебраических уравнений относительно трех констант С].

Запишем формальное решение этой системы уравнений в виде

(п) , _ ч

С = (рк, р).

Подставляя (47) в системы (42) и (43), находим:

(п) ( (п) . , _ Л

Кзк(°) = ^к У' ((,р),К/,£

,(°)=р

((,р),К/,£

икз(1) = ^к

У1 ((, р), К/, £

_ ц.

После подстановки третьего выражения

г^к

и осредненным градиентом г ,:

( (п) ( (п) Л Л

рк (°)=ц? + р (, р), цт, £

V V

_ Кз?

у у

8/з.

Подставляя выражения (48), (49) для р (°) в определяющие соотношения системы (16), находим компоненты тензора напряжений Пиолы — Кирхгофа в нулевом приближении

рУ (0) = _ Р

_а(0)У

(р/,р),р',% (р~*) +

(и) ( ( (и) ( (иК ч ^ > \

+^ Рк + SJ ((, р), р/, % _ р3к 8/3, %

V V V ) ) )

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

Г)1} ((, р )=< Г)1]

г

V V

Sj ((,р),р,%

\ \ >

_ П 8/3, % >

) ) )

> .

Алгоритм численной реализации предложенного метода решения локальной задачи и вычисления функции осредненных определяющих соотношений (51) аналогичен алгоритму, описанному для сжимаемых нелинейно-упругих композитов [8, 9].

Задача об одноосном растяжении слоистой композитной пластины. Рассмотрим задачу об одноосном растяжении пластины из слоистого композита. Для этой задачи средний тензор деформаций р I композита имеет три ненулевые компоненты, не зависящие от координат X1,

р = 1 + м1

Отсюда находим закон движения пластины из осредненного композита при одноосном растяжении

X7 =Р7}Хх, у = 1, 2, 3.

Из условия несжимаемости (27) следует, что независимыми являются только две компоненты р77:

,2

3

гк

р3 = 1.

Будем искать решение локальной задачи нулевого приближения (39)-(47) для случая одноосного растяжения в виде

рк

7=1

Тогда из системы (40) и условия несжимаемости слоев следует, что все градиенты деформации слоев совпадают со средним градиентом

Ру

Тогда система (41) принимает вид

С3 =_р (°)/Щ +/°33 ((, К1, К22, £,

р (°0 = р

и допускает решение в явном виде относительно С и р( ):

С3=_р (°)/ Кз3 + Гъъ ((, р1, р22, £

= р + Р 3

^°33

Уравнение (46) с учетом (56) выполняется тождественно. Формулы (58) вместе с и/3з(1) = ° представляют собой явный вид решения (48) для одноосного растяжения. Выражения (56) — это явная форма уравнений (49) для данного случая. Формула (5°) для напряжений в слоях композита принимает следующий вид:

рШ = _

Осредненное определяющее соотношение

_ _ / (п) _ _ _ \

р* =_

Выражение (37) для компонент осредненного тензора напряжений Коши принимает вид

Туу = _р

+ Ру^^?°)уу

Используем граничные условия на боковых поверхностях компо-

—22 —33

зита Т = Т = ° при одноосном растяжении. В силу того, что напряжения Туу в пластине однородны, из граничных условий и формул

Т11=_р+Р^ Г^11 ((, Р1, Р22, £ ° = _р + Ру^^((3,Ц1,Р22,£)), У = 1, 2,

Р1Щ2 Рз3 = 1

получаем систему четырех уравнений относительно четырех неизвестных: Т11, Р22, Рз3, р, являющихся функциями от р/. Решая эту систему, находим диаграммы деформирования

Т11 = Т11

композита при одноосном растяжении.

Результаты численного моделирования. В качестве конкретных моделей нелинейно-упругих сред была рассмотрена модель несжимаемых сред класса Бу [ 16] (модель Муни, для которой индекс п

принимает значение 5). Тензорная функция Г(Р/ , £) для нее имеет явное выражение через градиент деформации

где ¡1 = р^р? — первый инвариант меры деформации Ог/ = ркрк1; д(£) и Р(£) — упругие константы слоев. Вычислять собственные

°  / \

векторы бУ (Р/ ) и собственные значения для этой модели не требуется.

Для численного решения системы (62) был использован метод минимизации численной ошибки путем покоординатного спуска [8, 9]. Ячейка периодичности композита состояла из трех слоев. Были выбраны следующие упругие константы слоев:

_ слой 1: = °,8 МПа, Р1 = 1, °, относительная толщина слоя

¿1 = °, 2;

_ слой 2: ^2 = 1,4 МПа, Р2 = 1, °, относительная толщина слоя

И = °,6;

_ слой 3: = °, 7 МПа, Рз = °, 9, относительная толщина слоя

Из = °, 2.

Графики функций (63), построенные с помощью описанного выше алгоритма для отдельных слоев композита и для самого композита, отражены на рис. 1.

Рп> МПа

Рис. 1. Диаграммы деформирования Р11 ) для слоистого композита (К) и отдельных составляющих его слоев (1, 2 и 3), построенные с помощью разработанного метода

Функция Р1 (р11 для композита располагается между соответствующими функциями для отдельных его слоев: слоя 2 с более высокой жесткостью и слоев 1 и 3 с меньшей жесткостью, чем сам композит.

Были проведены расчеты для трех значений относительных толщин слоев композита: Ъ^ = 0,8, ^ = 0,6 и Л2 = 0,4, где ^ — толщина «жесткого» слоя, толщины слоев 1 и 2 были одинаковы. Результаты расчетов приведены на рис. 2. С увеличением толщины жесткого слоя диаграмма деформирования композита приближается к диаграмме деформирования жесткого слоя.

На рис. 3 представлены результаты расчетов напряжений Р1 в слоях композита с толщиной ^ = 0,6 при фиксированном значении р1 = 2,0 . Напряжения постоянны по каждому слою, но меняются при переходе от слоя к слою; наибольшие их значения возникают в жестком слое, превышая среднее напряжение Рп ) в композите

л, ().

МПа

Рис. 2. Диаграммы деформирования ^ (ц1 ) для слоистого композита, построенные с помощью разработанного метода, при разных относительных толщинах Н2 второго (жесткого) слоя:

1 — 0,8; 2 — 0,6; 3 — 0,4

Рп,МПа

8,0

6,0

4,0

-0,5 -0,4 -0,3 -0,2 -0,1 0 0,1 0,2 0,3 0,4 Е,

Рис. 3. Распределение напряжения ^ по толщине слоистого композита с ЯП из трех слоев при продольном растяжении

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

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

Приведенный пример построения диаграммы деформирования композита с конечными деформациями при одноосном растяжении

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

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

ЛИТЕРАТУРА

[1] Christensen R.M. Mechanics of composite materials. New York, John Wiley & Sons, 1979, 324 p.

[2] Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. Москва, Наука, 1984, 356 с.

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

[4] Dimitrienko Yu.I. A Structural thermos-mechanics model of textile composite materials at high temperatures. Composites Science and Technology, 1999, vol. 59, issue 7, pp. 1041-1053.

[5] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Конечно-элементное моделирование эффективных вязкоупругих свойств однонаправленных композиционных материалов. Математическое моделирование и численные методы, 2014, № 2, с. 28-48.

[6] Димитриенко Ю.И., Кашкаров А.И. Расчет эффективных характеристик композитов с периодической структурой методом конечного элемента. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2002, № 2, с. 95-108.

[7] Димитриенко Ю.И., Кашкаров А.И., Макашов А.А. Конечно-элементный расчет эффективных упругопластических характеристик композитов на основе метода асимптотического осреднения. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2007, № 1, с. 102-116.

[8] Димитриенко Ю.И., Губарева Е.А., Кольжанова Д.Ю. Моделирование слоистых композитов с конечными деформациями методом асимптотической гомогенизации. Инженерный журнал: наука и инновации, 2015, вып. 5. URL: http://engjournal.ru/catalog/msm/pmcm/1405.html (дата обращения 24.05.2017). DOI 10.18698/2308-6033-2015-5-1405

[9] Димитриенко Ю.И. Моделирование нелинейно-упругих характеристик композитов с конечными деформациями методом асимптотического осреднения. Известия высших учебных заведений. Машиностроение, 2015, № 11, с. 68-77. DOI 10.18698/0536-1044-2015-11-68-77

[10] Yang Q., Xu F. Numerical modeling of nonlinear deformation of polymer composites based on hyperelastic constitutive law. Frontiers of Mechanical Engineering in China, 2009, vol. 4, issue 3, pp. 284-288.

[11] Aboudi J. Finite strain micromechanical modeling of multiphase composites. International Journal Multiscale Computational Engineering, 2008, vol. 6, no. 5, pp. 411-434.

[12] Zhang B., Yu X., Gu B. Micromechanical modeling of large deformation in se-piolite reinforced rubber sealing composites under transverse tension. Polymer Composites, 2015. DOI 10.1002/pc.23596

[13] Ge Q., Luo X., Iversen C.B., Nejad H.B., Mather P.T., Dunn M.L., Qi H.J. A finite deformation thermomechanical constitutive model for triple shape

polymeric composites based on dual thermal transitions. International Journal of Solids and Structures, 2014, vol. 51, pp. 2777-2790.

[14] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4. Основы механики твердых сред. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2013, 624 с.

[15] Димитриенко Ю.И., Даштиев И.З. Модели вязкоупругого поведения эластомеров при конечных деформациях. Вестник МГТУ им. Н. Э. Баумана. Сер. Естественные науки, 2001, № 1, с. 21-41.

[16] Димитриенко Ю.И. Нелинейная механика сплошной среды. Москва, Физ-матлит, 2009, 610 с.

[17] Димитриенко Ю.И., Яковлев Н.О., Ерасов В.С., Федонюк Н.Н., Сборщиков С.В., Губарева Е.А., Крылов В.Д., Григорьев М.М., Прозоровский А.А. Разработка многослойного полимерного композиционного материала с дискретным конструктивно-ортотропным заполнителем. Композиты и наноструктуры, 2014, т. 6, № 1, с. 32-48.

[18] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой. Математическое моделирование и численные методы, 2014, № 1, с. 36-57.

[19] Димитриенко Ю.И., Яковлев Д.О. Асимптотическая теория термоупругости многослойных композитных пластин. Механика композиционных материалов и конструкций, 2014, т. 20, № 2, с. 260-282.

Статья поступила в редакцию 31.01.2017

Ссылку на эту статью просим оформлять следующим образом:

Димитриенко Ю.И., Губарева Е.А., Кольжанова Д.Ю., Каримов С.Б. Моделирование слоистых несжимаемых композитов с конечными деформациями методом асимптотического осреднения. Математическое моделирование и численные методы, 2017, № 1, с. 32-54.

Димитриенко Юрий Иванович — д-р физ.-мат. наук, профессор, директор Научно-образовательного центра «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана, заведующий кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 350 научных работ в области механики сплошных сред, вычислительной механики, механики и термомеханики композитов, математического моделирования в науке о материалах, вычислительной газодинамики. e-mail: dimit.bmstu@gmail.com

Губарева Елена Александровна — канд. физ.-мат. наук, доцент, заместитель заведующего кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 50 научных работ в области механики сплошных сред, механики контактного взаимодействия, математического моделирования, механики композитов. e-mail: gubareva_ea@pochta.ru

Кольжанова Дарья Юрьевна — ассистент кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. e-mail: darya.kolzhanova@yandex.ru

Каримов Суннат Бахтиярович — студент кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. e-mail: sunnat49@yandex.ru

Incompressible layered composites with finite deformations on the basis of the asymptotic averaging method

© Yu.I. Dimitrienko, E.A. Gubareva, D.Yu. Kolzhanova, S.B. Karimov Bauman Moscow State Technical University, Moscow, 105005, Russia

The article considers the modeling results of incompressible layered composites with finite strains deformation according to the individual layers characteristics. The article proposes an asymptotic averaging method version for layered nonlinearly elastic incompressible composites with finite deformations and periodic structure. We are using a universal representation of the defining relations for incompressible composite layers, proposed by Yu.I. Dimitrienko, which allows us to simulate simultaneously for a complex of various nonlinear elastic media models characterizedby the choice of a pair of energy tensors. We proved that if all composite layers are incompressible, the composite as a whole is also an incompressible, but anisotropic, medium. The article considers the problem of laminated plate uniaxial stretching from incompressible layers with finite deformations. Using the developed method, we calculated the effective deformation diagrams connecting the averaged Piola — Kirchhoff stress tensors components and the strain gradient, as well as the stress distribution in the composite layers. The developed method for calculating effective deformation diagrams and stresses in composite layers can be used in the design of elastomeric composites with specified properties.

Keywords: layered composites, finite deformations, incompressible media, asymptotic averaging method, Piola — Kirchhoff stress tensor, strain gradient, universal models of determining relations, deformation diagrams

REFERENCES

[1] Christensen R.M. Mechanics of composite materials. New York, John Wiley & Sons, 1979, 324 p.

[2] Bakhvalov N.S., Panasenko G.P. Osrednenie protsessov v periodicheskikh sredakh [Averaging processes in periodic media]. Moscow, Nauka Publ., 1984, 356 p.

[3] Pobedrya B.E. Mekhanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, Lomonosov Moscow State University Publ., 1984, 336 p.

[4] Dimitrienko Yu.I. Composites Science and Technology, 1999, vol. 59, issue 7, pp. 1041-1053.

[5] Dimitrienko Yu.I., Gubareva E.A., Sborschikov S.V. Matematicheskoe mode-lirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 2, pp. 28-48.

[6] Dimitrienko Yu.I., Kashkarov A.I. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki — Herald of Bauman Moscow State Technical University. Series Natural Sciences, 2002, no. 2, pp. 95-108.

[7] Dimitrienko Yu. I., Kashkarov A.I., Makashov A.A. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki — Herald of Bauman Moscow State Technical University. Series Natural Sciences, 2007, no.1, pp. 102-116.

[8] Dimitrienko Yu.I., Gubareva E.A., Kolzhanova D.Yu. Inzhenerny zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2015, no. 5. Available at: http://engjournal.ru/catalog/msm/pmcm/1405.html (accessed May 24, 2017). DOI 10.18698/2308-6033-2015-5-1405

[9] Dimitrienko Yu.I. Izvestiya vyshikh uchebnykh zavedeny. Mashinostroenie — Proceedings of Higher Educational Institutions. Machine Building, 2015, no. 11, pp. 68-77. DOI 10.18698/0536-1044-2015-11-68-77

[10] Yang Q., Xu F. Frontiers of Mechanical Engineering in China, 2009, vol. 4, no. 3, pp. 284-288.

[11] Aboudi J. International Journal Multiscale Computational Engineering, 2008, vol. 6, no. 5, pp. 411-434.

[12] Zhang B., Yu X., Gu B. Polymer Composites, 2015. DOI 10.1002/pc.23596

[13] Ge Q., Luo X., Iversen C.B., Nejad H.B., Mather P.T., Dunn M.L., Qi H.J. International Journal of Solids and Structures, 2014, vol. 51, pp. 2777-2790.

[14] Dimitrienko Yu.I. Mekhanika sploshnoy sredy. V 4 tomakh. Tom 4. Osnovy mekaniki tverdykh sred [Continuum mechanics. In 4 vols. Vol. 4. Fundamentals of mechanics of solid mechanics]. Moscow, BMSTU Publ., 2013, 624 p.

[15] Dimitrienko Yu.I., Dashtiev I.Z. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki — Herald of Bauman Moscow State Technical University. Series. Natural Sciences, 2001, no. 1, pp. 21-41.

[16] Dimitrienko Yu.I. Nelineynaya mekhanika sploshnoy sredy [Nonlinear mechanics of continuous medium]. Moscow, Fizmatlit Publ., 2009, 610 p.

[17] Dimitrienko Yu.I., Yakovlev N.O., Erasov V.S., Fedonyuk N.N., Sborschikov S.V., Gubareva E.A., Krylov V.D., Grigorev M.M., Prozorovsky A.A. Kompozity i nano-struktury — Composites andNanostructures, 2014, vol. 6, no. 1, pp. 32-48.

[18] Dimitrienko Yu.I., Gubareva E.A., Sborschikov S.V. Matematicheskoe mode-lirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 1, pp. 36-57.

[19] Dimitrienko Yu.I., Yakovlev D.O. Mekhanika kompozitsionnykh materialov i konstruktsiy — Mathematical modeling and Computational Methods, 2014, vol. 20, no. 2, pp. 260-282.

Dimitrienko Yu.I., Dr. Sc. (Phys.-Math.), Professor, Head of the Scientific and Educational Center for Supercomputer Engineering Modeling and Development of Software Complexes, Bauman Moscow State Technical University, Head of Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University. Author of over 350 research papers in the field of continuum mechanics, computational mechanics, mechanics and thermomechanics of composites, mathematical modeling in the science of materials, computational gas dynamics. e-mail: dimit.bmstu@gmail.com

Gubareva E.A., Cand. Sci. (Phys.-Math.), Assoc. Professor, Deputy Director of Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University. Author of over 50 research papers in the field of continuum mechanics, the mechanics of contact interaction, mathematical modeling, and the mechanics of composites. e-mail: gubareva_ea@pochta.ru

Kolzhanova D.Yu., assistant of Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University. e-mail: darya.kolzhanova@yandex.ru

Karimov S.B., student of Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University. e-mail: sunnat49@yandex.ru