Конечно-элементный расчет эффективных упругопластических характеристик композитов на основе метода асимптотического осреднения

ПРИКЛАДНАЯ МАТЕМАТИКА И МЕТОДЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

УДК 539.3

Ю. И. Димитриенко, А. И. Кашкаров, А. А. Макашов

КОНЕЧНО-ЭЛЕМЕНТНЫЙ РАСЧЕТ ЭФФЕКТИВНЫХ УПРУГОПЛАСТИЧЕСКИХ ХАРАКТЕРИСТИК КОМПОЗИТОВ НА ОСНОВЕ МЕТОДА АСИМПТОТИЧЕСКОГО ОСРЕДНЕНИЯ

Разработан метод решения задач теории малых упруго-пластических деформаций А.А.Ильюшина на ячейке периодичности для композитов с пространственной структурой армирования, основанный на сведении исходной задачи к задаче "классического" типа, для которой применен метод конечных элементов. Разработана численная процедура, и решена задача расчета эффективных упругопластических характеристик для пространственно-армированного композита.

Метод асимптотического осреднения периодических структур, идея которого была предложена в 1974 г. Н.С. Бахваловым и независимо Э. Санчес-Паленсией (см. список литературы к работам [1], [2]), позволяет математически точно вычислять так называемые эффективные характеристики композиционных материалов, как линейные, так и нелинейные. Применение этого метода для расчета упругих модулей композитов, а также для расчета нелинейно-упругих и упруго-пластических свойств было впервые предложено Б.Е. Победрей [3]. Было показано, что для этих целей необходимо решение специальных задач на ячейке периодичности композита. Однако методы решения этих задач до настоящего времени удалось разработать только для очень ограниченного типа структур — слоистых и однонаправленно-армированных. В недавних работах [4-7] был предложен новый метод решения задач на ячейке периодичности, основанный на сведении их к задачам классического типа, к которым уже может быть применен метод конечного элемента. С помощью этого метода были рассчитаны упругие модули пространственно-армированных композитов. Цель настоящей работы — дальнейшее развитие этого метода для расчета упруго-пластических характеристик композитов со сложными структурами армирования.

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

Рис. 1. Композиционный материал с пространственной (3-0) схемой армирования

которому в пространстве К3 соответствует область V с поверхностью £. Композит состоит из N компонентов. Компоненты с индексами а = 1,..., N—1 представляют собой волокна различных типов, ориентированные по а-му направлению в М3, а компонент с индексом а = N — это матрица. Обозначим Va, а = 1... N, — области в М3, соответствующие а-му компоненту композита; £а — поверхности областей Va и — поверхности контакта матрицы и волокон (волокна полагаем не контактирующими между собой); £ае — часть поверхности £ композита, занятая а-м компонентом (причем £а = £aN и £ае —

N-1

для волокон и £N = У и £Ne — для матрицы). Волокна и ма-

а=1

трицу полагаем изотропными упруго-пластическими, соответствующими деформационной теории пластичности А.А. Ильюшина (теория малых упруго-пластических деформаций) [8]. Тогда в каждой области Уа, а = 1... N, можно рассмотреть следующую задачу малых упруго-пластических деформаций:

Здесь м",

о

_ 0,

в Va;

_a _ Ipa (ca N

<ij _ Fij \bkl) ,

1

2

_a _ _

bij _ -

aa i,j _Г j,i

В VaU Ea

Va

a N

u'a — uN,

ua ui

<jj- oN

ij ij

nj _0

на El

<ij nj Sie,

на EaN;

на EL.

<

— перемещения, напряжения и деформации в a-м

a

в

компоненте, последние два условия в системе (1) — это условия идеального контакта матрицы и волокон, и условия на а.

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

В этом случае для такой структуры может быть применен метод асимптотического осреднения [1-5], согласно которому решение задачи (1) для матрицы и волокон строится в виде асимптотических разложений [5]:

— иг

^ = <(0)(хк 1) + ка°(1)(хк 1) + к2 ...,

причем по аргументу эти функции полагаются периодическими. Деформации и напряжения "нулевого уровня" имеют следующий вид:

е?0) = % + !

1 / .

% = 2

_а

сгз _ Ач I е

V

Здесь ,1 = д/дх1 и / = д/д£1 — производные по двум типам координат. При выводе формул (3)-(5) и далее используется правило дифференцирования асимптотических разложений:

д/ ^ д/ (хг+ д/ (хг,^г) д6 = д/ ^ ^ + 1 д/ ^ ^

дх^ дх^ д£к дх^ дх^ к д£к

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

<

a(0) ij/j

_ 0, В V*,

<a(0) _ Fa /_a(0) <<ij _ Fij \fckl

В V* и Es,

_a

"ij

- , 1 f a

_ + 2 (V/ + jV

в V*,

a

ui _ ui

<

u

ij

a

на E

*aN,

) nj

_ 0,

u

a

_ 0.

Здесь

(< > _ E / uadV*, a=vi

j _

N

e

a=l

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

В выражении (6) условие [и"] = 0 — это условие периодичности функций на границе ячейки периодичности, а условие (и") = 0 вызвано требованием единственности решения локальной задачи [3]. В силу периодичности функций иа( 1 ^ имеет место следующее соотношение:

£ij

_a

ij

_ 1

2 Vui>j + uji

Метод упругих решений для локальных задач Жр<1 на ячейке периодичности.. Будем полагать далее, что волокна и матрица — изотропные, а ЯП является симметричной при зеркальном отражении относительно трех координатных плоскостей £ = в = 0}, а также симметричной относительно поворота на угол п вокруг каждой оси координат О^и симметричной при преобразовании центральной симметрии с центром в точке О (подробнее об этих преобразованиях см. работу [9]). Это ограничение назовем основным допущением симметрии композита. Согласно этому допущению, вместо решения локальной задачи (6) на всей ЯП V, можно перейти к решению задачи на области V, представляющей 1/8 ЯП.

Поскольку задача (6) нелинейная, то для ее решения применим итерационный метод, являющийся разновидностью метода упругих решений [8]. Согласно этому методу определяющие соотношения в системе (6) для компонентов композита линеаризуются следующим

0

образом:

a{m} _ fYa(m> a{m} тга (ca{m~1}\ _ a{m-1}

где а^ и е ^ — значения напряжений а^ и деформаций е ^ на т-м шаге итерационного цикла, а Сак1 — тензоры модулей упругости компонентов композита на т-м шаге итерации. В простейшем варианте метода эти модули являются константами и совпадают с Сак1.

г^г а{ш| « а

Обозначим также — значения перемещений на т-м шаге

итерационного цикла. Тогда на т-м шаге итерации вместо задачи (6) получаем следующую линеаризованную задачу:

f a{m} n j г

аЛ. =0 в Vt,

J/J 5'

a{m} _ r,a{m> a{m^ раЛa{m-l}\ ^aim} a{m-l^ ^ . . I I V

ег/ } = + 2 (Vj } + V» J в ,

ua{m} = Mf{m}, ^

на S taN.

(af} - aN {m})nj = 0

Согласно предложенному в работе [5] варианту метода асимптотического осреднения, перемещения первого уровня и напряжения нулевого уровня при каждом значении m представляются в виде следующих сумм :

3 3 3 a{m} _ \ Л a{m} a{m} _ \ Л a{m} a{m} _ \ Л a{m} Ai 1\

= ui

p,q=1 Pi9=1 Pi9=1

если £k E Vai, a = l,..., N, причем для функций в каждой

комбинации (pq) выделяется линейная часть по локальным координатам:

u

am=-p.*

где — символ Кронекера, а (6) — некоторые функции, назы-

ваемые псевдоперемещениями, для которых при каждом фиксированном наборе индексов (рд) получаем следующую задачу, называемую "линеаризованной локальной задачей Жрд на ячейке периодичности":

a{m} т/

jq)}j B '

a{m} _ r,a{m> a{m} r,a{m> a{m-l} 1 ~a{m-1} ^ .. v' I I V

¿j(pq) ijki ¿j(pq) ijki j(pq) g

ea{(m> - 1 fUa{m}. + Ua({m},) в К, j(pq) 2 V i(pq)/j j(pq)}y

Ua{m} — иN{m} ¿(pq) ¿(pq) ,

f a{m} N {m}\ n f На S «"N•

(j} - — 0 I

Здесь обозначено:

^ = ^ ^

\Р,5=1 /

= 1 + . (15) Кроме того, к системе (13) присоединяются условия на координатных плоскостях £ = = 0} и на торцевых поверхностях ЯП £'3 = = 1/2}, 8 = 1, 2, 3, которые записываются следующим образом:

пРи p — q ^pq) — 1/2epq ^ jpq) — 0, Sfc(pq) — 0 На ,

i — j — k — i

пРи p — q^pq) — j ) — 0 Ua(pq) — 0

на ,i,j — {p,q};

Si(pq) — 0, Sj(pq) — 0, Uk(pq) — 0 На , г — j — k — ^

а также

пРи p — q, Uapq) — 0, S(pq) — 0 Sa(pq) — 0 На

г — j — k — i

пРи p—q ^pq) — 0, jpq) — 0, Ua(pq) — 0 На ,i,j — { Sa(pq) — 0, jpq) — 0, Ua(pq) — 0 На , i — j — k —

Вывод этих граничных условий для линеаризованных задач (13) такой же, как и полученный в работе [5], для линейных задач на ячейке периодичности. Здесь введены обозначения для векторов усилий б1^).'

^

Решение задачи (13)-(17) разыскивается в области V;, представляющей 1/8 ЯП: V; = V П > 0), здесь также обозначена Ё— поверхность контакта компонентов внутри V;: Ё= Ё;а(з П V.

Появление условий (16), (17) связано с требованием периодичности всех функций (1) на границе ЯП, а также наличием симметрии ЯП относительно трех координатных плоскостей для рассматриваемого типа композита. Более подробно эти условия обсуждаются в работе [5]. Там же на рисунках приведены граничные условия для задач Ж33, Ж13 и Ж31.

Вариационная формулировка локальной задачи Жрд. Для произвольного конечного объема У"с V вариационная формулировка задачи Жрд (13) при фиксированных р и д имеет вид

У ¿ета^ = У Б^Е + У ¿ета{т-1}^.

V Е V

Здесь обозначены координатные столбцы псевдоперемещений и, напряжений а, деформаций е и поверхностных усилий Б:

U =

тт-ajm} j j-a{m} Tj-a{m} U1

а =

a{m} a{m} a{m} a{m} ; ri1(pq)' a22(pq)' a33(pq)' 713(pq)/V 2

7.

a{m} 23

/72

7

a{m} 12

/72

e =

a{m} a{m} a{m} a{m}

11

)/T2,

/А Ä/Vä ]т ,

-23

т

т

r< _ Г na na na 1т

S = LS1(pq),S2(pq),S3(pq^ ,

7{m-1} =

^a{m-1} ^a{m-1} ^a{m-1} ^a{m-1} ^a{m-1} ^a{m-1} 711

^ajm— 1} ^-yajmj a-fm— 1} ~ajm—1| <7 / ч = 0-1, £,,/ ч--Oi

j(pq) ki(pq) g

Линеаризованные определяющие соотношения, следующие из второй группы уравнений в системе (13), с использованием координатных столбцов записываются следующим образом:

а = С е,

где С — матрица упругости размером [6 х 6], составленная из компонентов тензора С,а{т} стандартным образом [9]:

т

C =

' C1111 C1112 C1113 0 0 0

/~t a{m} C2222 C1123 0 0 0

C a{m} C3333 0 0 0

сим. 2C a{m} 2C1212 0 2C1313 0 0

V 2C2323

/

Соотношения Коши (третья группа уравнений в системе (13)) в матричном виде записываются следующим образом:

е = DU,

где Б — матрица линейных дифференциальных операторов дифференцирования (д = 1):

D=

Ö1 0 0

0 d2 0

0 0 d3

^/78 0 Ö1/V8

0 дз/78 02/78 _ ^/78 дх/^8 0

С учетом определяющих отношений (23) и (25) вариационное уравнение (21) можно представить в виде

[

./У Уе

Метод конечного элемента для задач Жрд. Аппроксимируя псевдоперемещения и в КЭ линейными функциями

и = Фд,

где д — координатный столбец псевдоперемещений в узлах КЭ, а ф (^) — матрица функции формы, зависящая от типа КЭ, получаем итоговую разрешающую систему линейных алгебраических уравнений:

Кд = /,

где

К = у Вт СВ^

V

— локальная матрица жесткости, В = ДФ, а

/ = I Фт^£ + у Вта¿V

Е V

— столбец нагрузок.

Глобальная матрица жесткости задачи составляется из локальной матрицы жесткости стандартным образом [10], после ее формирования к системе линейных алгебраических уравнений (СЛАУ) применяются граничные условия (16) и (17). Граничные условия идеального контакта (последняя группа соотношений в системе (13)) не требуют специального учета, так как при данном варианте МКЭ они удовлетворяются автоматически.

Решая СЛАУ, находим перемещения д в узлах, по которым вычисляем псевдоперемещения и = Фд, деформации е = Вд и напряжения а{т} = сВд — <г{т-1} в КЭ. Для решения СЛАУ применялся метод сопряженных градиентов.

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

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

Е

I

3 X -1

E|_(a){m-1}

Гвт

^ 8,

,3,7=1 ' ' \3,Т=1 ' ' '

где 8 = 0,01... 0,001. Число итераций т = М, обеспечивающее выполнение данного условия, было различным для разных задач Жрд и для разных значений ерд, но не превышало 10... 15. Напряжения при максимальном значении номера итерации т = М обозначались как а^а^РМ)}. Для повышения устойчивости счета задачи Ж11, Ж22 и Ж33 решались совместно — как общая задача с входными данными е11, е22и е33. Также совместно решались задачи Ж13 и Ж31, Ж12 и Ж21, Ж23 и Ж32. Далее совершался еще один цикл решения задач Жрд — по значениям входных данных .

Расчет эффективных упруго-пластических характеристик композиционного материала. После решения серии задач Жрд (13), (16), (19) указанным методом для всех pq проинтегрируем напряжения ав7(РМ)} по областям, занятым волокнами и матрицей,

а, =

p,q=i

В результате получим осредненные напряжения в локальной упруго-пластической задаче (6). Тогда эффективные упруго-пластические соотношения, связывающие осредненные напряжения а, и деформации можно записать в следующем символическом операторном виде:

а, = Fj

Поскольку ранее указан алгоритм вычисления средних напряжений а, (формулы (33), в которых } вычисляются по методу, описанному

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

Данный символический оператор (34) может быть конкретизирован, если имеется информация о типе геометрической симметрии ЯП композита и типе анизотропии волокон и матрицы. Поскольку принято основное допущение симметрии, то решение всех задач Жрд будет обладать указанными типами симметрии (т.е. не изменится во всей ЯП при указанных выше преобразованиях). Это означает, что и оператор (34) должен обладать данным типом симметрии. Но перечисленные выше преобразования вместе с тождественным преобразованием образуют группу ортотропии (8 преобразований [9]), следовательно по терминологии, введенной в работе [9], оператор (34) будет тензорной функцией, индифферентной относительно группы ортотропии. Но тогда для такой функции можно использовать представление ее в тензорном базисе группы ортотропии, которое имеет вид

3 1

а, = Fj (£pq) = Е j + ц(^4ё23 + ^бё12ё1з)(^г2^,3 + ¿г3£,2) +

7=1

+ ^(^5 £13 + ^бё12ё2з)(£г1£,3 + ¿¿3^1) +

+ 2

где

=

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

= еШ = а22, ^ = е33, ^ = е;23,

£ _2 2___ £ _2

= е 13, = е 13е 12е 23, = е 12.

Эти функции <7 фактически и представляют собой искомые эффективные упругопластические характеристики композита. Зная значения и , функции <7 могут быть вычислены следующим образом. Запишем формулы (35) в явном виде:

а 11 = <£1(/1,...,/1),

а22 = Ы/ъ ... ,/г),

а33 = ^1

а 23 = <4е23 + <6е 12е 13,

а 13 = <5е13 + <6а12е23,

а 12 = <7^12 + ^бё13ё23,

а затем образуем из этих соотношений семь инвариантов тензора средних напряжений , 7 = 1,..., 7, аналогичных инвариантам 17, 7 = 1,..., 7 (37), в результате получим систему семи нелинейных уравнений относительно семи функций <7:

/Г = <1, /2 = <2, /3 = <3,

/Г = <4/7 + 2^6^4/17/7 + <2/б7, /Г = ^5/7 + 2^6^5/4/7 + ^2/7, /Г = ^7/7 + 2^7^6/7/7 + <6 /7,

/6г = (<4 <5<7 + <2(< + <5 + ^г))/| + <5<6<7 /52/72 +

+ <4<6<7/472/72 + <4<5<6/|2/|2 + <6/| 2, решая которую находим искомые эффективные упругопластические функции композита (34).

Теорема. Пусть для композита с периодической структурой выполнены следующие условия: 1) принято основное допущение симметрии, 2) волокна и матрица являются квазилинейными упруго-пластическими изотропными материалами, т.е. их операторы (5) ^ (е^) = ^ /7<а)^ зависят только от двух инвариантов тензо-

ра деформаций относительно группы изотропии:

т£(а) __7"е(а) —

= е11 + е22 + е33, ^2 =

=

и не зависят от третьего инварианта в этой группе /|(а) = det(—а), тогда эффективные упругопластические функции (33), (34) композита являются: а) индифферентными относительно группы ортотропии, б) квазилинейными: зависящими только от линейных и квадратичных инвариантов тензора осредненных деформаций /|, /|, /|, /|, /| Т£ "

'2> ^З' -Ч» ±Ъ' 7

I6

вид:

не зависящими от кубического инварианта т.е. имеют следующий

- _ 3 1 СТ^ = ^ (—) = ^ ^Т+ о^4—23(^2^3 + ¿¿3^2) +

1 ^=1 1 2

+ 2 ^5-13

а = 0.

Доказательство. Утверждение а) было нами доказано ранее. Эффективные упругопластические функции (34) композита полностью определяются видом символического оператора (34) эффективных определяющих соотношений, который в свою очередь определяется решением локальной задачи (13), (16), (17). Решение же этой задачи зависит от только через граничные условия (16), в которых деформации входят в виде комбинаций /| = - 11, /| = -22, /| = -33, = |-2з|, = | — 131, = | — 121, совпадающих с шестью

указанными инвариантами тензора —. Зависимости от кубического инварианта /| у данного оператора нет. Таким образом и утверждение б) теоремы, а следовательно, и сама теорема доказаны.

Если выполнены условия теоремы, то соотношения (39) для нахождения эффективных упруго-пластических функций упрощаются и имеют вид

= /", ^2 = /", ^3 = Ц,

^2 = //|, ^2 = //|, ^7 =

Для того чтобы воспользоваться этими соотношениями, необходимо: 1) в шестимерном пространстве инвариантов /|, /|, /|, /|, /|, /7 ввести область допустимых значений, например, шестимерный куб

= {атт < /| < атах, 7 = 1, 2, 3, 4, 5, 7}; ввести сетку = г7/, 7 = 1,..., 5, 7, где /^ — значения инварианта в узле сетки с шестимерным номером г7, 7 = 1,..., 5, 7, а / = (атах — атт)/^ — шаг сетки; 2) для каждого значения /^, 7 = 1,..., 5, 7 в узле сетки решить серию задач Жрд и, получив значения осредненных напряжений (33) СТ^, найти значения их инвариантов ; 3) по формулам (43) найти значения ¿7 функций .

Данный способ приводит к значительному числу вычислений, так, если введена сетка с числом узлов на ребре куба N=10, то общее

число узлов сетки составит 106, а число необходимых для решения задач Жрдсоставит 6 • 106. Однако, если в расчетах используется модель малых упруго-пластических деформаций А.А. Ильюшина, то указанное число задач Жрд может быть существенно сокращено.

Случай малых упруго-пластических деформаций А.А.Ильюшина. В случае малых упруго-пластических деформаций пластичности А.А.Ильюшина [8] операторы Fj (е^) в системе (1) для волокон и матрицы имеют вид

2Ga

где Ga — модуль сдвига, Ka - модуль объемного сжатия, — функция пластичности А.А. Ильюшина, причем линеаризованные определяющие соотношения (9) на m-м шаге будут иметь следующий вид:

<{m}= (3Ga (1 - <{m}) - K«) e^f}¿у+2Ga (1 - <{m}) e^-

-(2Ga (1 - <{m}) - Ke) eOf-% +

+ 2Ga

Соотношения (42) не зависят от траектории в пространстве деформаций, и часто решения прикладных задач ограничиваются случаем простого (пропорционального) деформирования, при котором (t) = ^(t), здесь t — параметр кривой в пространстве деформаций £pq. Тогда в шестимерном пространстве инвариантов этой кривой будет соответствовать луч /^(t) = /0Ф^), где /0 выражаются через по формулам (37). Вводя сетку вдоль этого луча /^ = гФ/0, Y = 1,..., 5, 7, где Ф = Фтах/N — приращение параметра деформирования (шаг сетки), получим, что число задач Жрд, необходимое для расчета эффективных упруго-пластических функций составит 6N.

Случай кубической симметрии. Для композита, у которого армирование осуществлено одинаковым образом относительно всех трех осей координат, ячейка периодичности будет иметь кубическую симметрию [9], а число независимых инвариантов тензоров /£ и равно трем. В качестве таких инвариантов могут быть выбраны спектральные инварианты [9], которые образуются с помощью спектрального разложения симметричных тензоров второго ранга относительно кубической группы симметрии [9]:

1 3 1

Pij(1) = 3akk, Pij(2) 'У ] ^aaAiaAja — Ц^kk,

a=l

3

Pj(3) = ^ ] ^aaAiaAja,

a=1

Q

1

3

Pij(1) — 1 ekk, Pij'(2) — ^ ^ eaaAiaAja „ ekk,

a=1

3

Pij(3) ej ^ ] eaaAia¿ja.

a=1

Спектральные инварианты и /f представляют собой свертки

)<г ре :

aa

соответствующих ортопроекторов Pj(a), Pi

ja — . /P.<j

a у j(a) j(a)' a у j(a) ij(a)' ' ' ' v '

и имеют следующий вид:

j. — 1 0"Ü, j. — ^ 7(an - 022 ) 2 + (^22 - ^33)2 + (^33 - ^ll)2,

J. — \[2\/a22 + ^3 + ^23,

Ii — 3J2 — ^ 7(e 11 - £22)2 + (e 22 - £ 33)2 + (e 33 - e 11)2,

/| = ^—1 + —13 + —223.

Тогда эффективные определяющие соотношения (34) принимают следующий вид:

- (- ) — + /Л А " --7"еЛ 1 +

ij Fij (epq) ^ Aij + J / ^ ^öiaöjae aa 3 j1 J +

2 a=1

у/3 j /е / у ^"i^j^aa

3

+ I3 ^ij - E ¿ia ¿ja-aa) ,

3 a=1

где функции связывают спектральные инварианты

/а =

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

Функции (/|, /|, /|) находились численным образом, процедура их вычисления была следующей. По формуле (33) вычислялись

средние напряжения с помощью решения соответствующих задач Жрд при фиксированных значениях , затем вычислялись значения инвариантов и по формулам (47), тогда значение = и является значением функции при фиксированных значениях . Меняя затем значения , получаем новые значения функции <^а.

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

/г = А^, /: = ^

где Аа — эффективные константы упругости композита, а (/^) — эффективные функции пластичности А.А. Ильюшина. Ниже приведены результаты проверки применимости простейшей модели пластичности для композитов с кубической симметрией.

Результаты численного моделирования упруго-пластического деформирования композита. При вычислительном эксперименте рассматривался пространственно-армированный углеалюминиевый композит, с алюминиевой матрицей и углеродными волокнами, которые для простоты предполагались изотропными и упругими. Константы упругости волокна были взяты следующими: Е = 250 ГПа, V = 0,35. Матрица полагалась упруго-пластической, функция пластичности выбиралась в виде

Г 0, ^ < /^

<

( (1 - к(Л0)(1 - /^Д^), /^ > /^ ,

который соответствует кусочно-линейной зависимости инвариантов ) = ) ^. Характеристики матрицы следующие:

Е = 70 ГПа, V = 0,25, /2 ) = а^)/Е), где а^) = 300 МПа —

Е<*)

предел текучести матрицы, к(м) = ) = 0,1 — отношение модулей упругости в пластической и упругой зонах.

Результаты моделирования приведены на рис. 2-9.

На рис. 2 и 3 показаны эффективные диаграммы пластичности: напряжение а33-деформация ё33 и напряжение ¿^-деформация ё13 для композита с коэффициентом армирования 0,29, а также для волокна и матрицы. Из этой диаграммы видно, что даже в направлении армирования композит проявляет пластические свойства, причем предел

"33

МПа 1200

800

400

/

/

У с

111

0.5

1.0

1.5

Езз °°

Рис. 2. Эффективная диаграмма пластичности:

напряжение ст33-деформация е33 для упруго-пластического композита с коэффициентом армирования 0,35

Рис. 3. Эффективная диаграмма пластичности:

напряжение ст13-деформация е13 для упруго-пластического композита с коэффициентом армирования 0,35

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

На рис. 4 приведена эффективная диаграмма (Г\3-€13 для упруго-пластического композита, полученная разработанным методом, а также аналогичные диаграммы, построенные приближенно-аналитическими методами по Фойгту (деформации сдвига у матрицы и волокна равны, а напряжения сдвига суммируются) и по Рейссу (наоборот, напряжения сдвига у матрицы и волокна равны, а деформации сдвига суммируются). Эффективные определяющие соотношения по Фойгту и по Рейссу записываются следующим образом:

Рис. 4. Эффективные диаграммы:

напряжение ст^-деформация £13 для упруго-пластического композита

Рис. 5. Распределение напряжений в 1/8 ЯП упруго-пластического композита в задаче Ж13+Ж31, коэффициент армирования 0,35

1 - <Pf

(1 - Шт(ст1з)) + 2Gt(1 - wf (ai3)))ai3 (по Рейссу),

V 2Ст

^13 = 2(Ст(1 - Рг)(1 -(ё1э)) + Gf (1 -Шf (¿Г1з)))-13 (по Фойгту), где Ст, шт(ё13), Gf, Шf (ё13) — модули сдвига и функции пластичности А.А. Ильюшина для матрицы и волокна. Проведенные расчеты показали, что вычисления эффективных диаграмм деформирования композита по этим приближенно-аналитическим методам приводят к значительной погрешности (см. рис. 4); отличие от результатов, получаемых разработанным методом, достигает 50 % по напряжениям.

На рис. 5 показано распределение напряжений а13 в 1/8 ЯП упруго-пластического композита в задаче Ж13+Ж31, а на рис.6 — распределение напряжений а33 в задаче Ж33.

Эффективная зависимость инварианта Ц от Ц (функция (49) вида = 0, 0)) для упруго-пластического композита с коэффициен-

том армирования 0,35 показана на рис. 7. Зависимости инвариантов

Рис. 6. Распределение напряжений СТ33 в 1/8 ЯП упруго-пластического композита в задаче Ж33, коэффициент армирования 0,35

0.5 1.0 1.5

Рис. 7. Эффективная зависимость инварианта от (функция

Ту = Ц, Ц)) для упруго-пластического композита с коэффициен-

том армирования 0,35

от при различных фиксированных значениях инварианта (функция /2 = (0,/|, /|) приведены на рис.8, а зависимости инварианта /д от /3 при различных фиксированных значениях инварианта (функция /3 = ^з(0,/|,/|) показаны на рис.9. Этими вычислениями проверялась возможность принятия простейшей модели (49) для композита. Они показали, что такая модель может быть принята, так как третий инвариант /3 практически не влияет на зависимость между вторыми инвариантами /2 и /|, и наоборот, второй инвариант /| практически не оказывает влияния на эффективную зависимость между третьими инвариантами. Этот вывод, по крайней мере, справедлив для исследованного диапазона изменений характеристик матрицы, волокна и коэффициента армирования.

Рис. 8. Эффективные зависимости инварианта I% от /| при различных значениях инварианта Ц для упруго-пластического композита с коэффициентом армирования 0,35

т° BOO х3

мш

700 BOO 500 400

эоо|-.................i.............../...............í....................................;..................j................-

200 ..................уГ.......... ..................j....................................................... .................-

100 -............................................;.................................... ...................................-

_I_¡_I_i_i_I_

0 0.002 0.004 0.Ü0B 0.000 0.01 0.012 .£

'з

Рис. 9. Эффективные зависимости инварианта If от If при различных значениях инварианта If для упруго-пластического композита с коэффициентом армирования 0,35

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

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

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

С помощью разработанного численного метода проведена проверка применимости простейшей модели анизотропной пластичности Б.Е. Победри, показавшая адекватность этой модели для исследованного типа композитов. Это позволяет ее рекомендовать для проведения расчетов конструкций из упруго-пластических композитов.

СПИСОК ЛИТЕРАТУРЫ

1.Бахвалов Н. С., Панасенко Г. П. Осреднение процессов в периодических средах. - М.: Наука, 1984. - 352 с.

2. Санчес-Паленсия Э. Неоднородные среды и теория колебаний. - М.: Мир, 1984. -472 с

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

4. Димитриенко Ю. И. Механика композиционных материалов при высоких температурах. - М.: Машиностроение, 1997. - 366 с.

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

6. Димитриенко Ю. И., Кашкаров А. И., Макашов А. А. Конечно-элементное моделирование процесса разрушения пространственно-армированных композитов с периодической структурой / Современные естественно-научные и гуманитарные проблемы. - М.: Логос, 2005. - С. 485498.

7. Д и м и т р и е н к о Ю. И., К а ш к а р о в А. И., М а к а ш о в А. А. Разработка конечно-элементного метода решения локальных задач теории упругости "на ячейке периодичности" для композитов с периодической пространственной структурой / Математика в современном мире: Под ред. Ю.А. Дробышева. -Калуга: Изд-во КГПУ, 2004. - С. 177-191.

8. И л ь ю ш и н А. А. Пластичность. - М.: ГИИТЛ. - 1948.

9. Димитриенко Ю. И. Тензорное исчисление. - М.: Высшая школа, 2001. - 575 с.

10. П о п о в Б. Г. Расчет конструкций вариационно-матричными методами. -М.:Изд-во МГТУ им. Н.Э. Баумана, 1993. - 294 с.

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

Юрий Иванович Димитриенко родился в 1962 г., окончил в 1984г. МГУ им. М.В.Ломоносова. Д-р физ.-мат. наук, профессор, заведующий кафедрой "Вычислительная математика и математическая физика" МГТУ им. Н.Э.Баумана, действительный член Академии инженерных наук. Автор более 100 научных работ в области вычислительной механики, нелинейного тензорного анализа, термомеханики композитов, математического моделирования в материаловедении.

Yu.I. Dimitrienko (b.1962) graduated from the Lomonosov Moscow State University in 1984. D.Sc.(Phys.-Math.), professor, head of "Computational Mathematics and Mathematical Physics" department of the Bauman Moscow State Technical University. Full member of the Russian Academy of Engineering Sciences. Author of more than 100 publications in the field of computational mechanics, nonlinear tensor analysis, thermomechanics of composite materials, mathematical simulations in science of materials.

Александр Игоревич Кашкаров родился в 1978 г., окончил в 2001г. МГТУ им. Н.Э. Баумана. Ассистент кафедры "Вычислительная математика и математическая физика" МГТУ им. Н.Э. Баумана. Автор ряда научных работ по численным методам в механике композитов.

A. I. Kashkarov (b.1978) graduated from the Bauman Moscow State Technical University in 2001. Assistant of "Computational Mathematics and Mathematical Physics" department of the Bauman Moscow State Technical University. Author of some publications in the field of numerical methods for mechanics of composites.

Алексей Анатольевич Макашов родился в 1983 г., окончил в 2005 г. МГТУ им. Н.Э. Баумана. Аспирант кафедры "Вычислительная математика и математическая физика" МГТУ им. Н.Э. Баумана. Имеет ряд научных работ по применению метода конечных элементов для решения задач механики композитов.

A. A. Makashov (b.1983) graduated from Bauman Moscow State Technical University in 2005. Post-graduate of "Computational Mathematics and Mathematical Physics" department of the Bauman Moscow State Technical University. Author of some publications in the field of finite element method applications for solving problems of mechanics of composites.