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

Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 11. С. 748-770.

Б01: 10.7463/1114.0734246

Представлена в редакцию: 09.08.2014 Исправлена: 27.10.2014

© МГТУ им. Н.Э. Баумана УДК 539.3

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

Димитриенко Ю. И.1' , Губарева Е. А.1, ацш.ьташаgmaii.com

1 7

Сборщиков С. В. , Федонюк Н. Н.

:МГТУ им. Н.Э. Баумана, Москва, Россия 2Крыловский государственный научный центр, Санкт-Петербург,

Россия

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

Ключевые слова: вязкоупругие характеристики, композиты, полимерные композиционные материалы, слоисто-волокнистые материалы, метод асимптотического осреднения, комплексные модули упругости, тангенс угла потерь

Введение

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

Наука и Образование

МГТУ им. Н.Э. Баумана

Сетевое научное издание

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

Среди работ, посвященных исследованию вязкоупругих характеристик ПКМ, отметим [1-8]. Диссипация энергии в материалах при вибрационном нагружении обусловлена внутренним трением, и с точки зрения термомеханики деформируемых твердых сред ("из первых принципов") это явление описывается моделями термовязкоупругих сред [2,5,9-11]. Разработке моделей вязкоупругого поведения композитов посвящены работы [12-18]. Целью настоящей работы является разработка метода расчета вязкоупругих характеристик ПКМ слоисто-волокнистой структуры, которые широко применяются при создании авиационных и судовых конструкций, исходя их диссипативных характеристик составляющих компонентов ПКМ - матрицы и армирующих волокон. Метод основан на применении теории асимптотического осреднения периодических структур [19-24] для вязкоупругих композитов при установившихся циклических колебаниях.

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

колебаниях

Рассмотрим задачу механики линейно вязкоупругой среды [1 1] при переменных нагрузках, при относительно невысоких частотах, когда инерционными силами можно пренебречь:

V-о = 0,

8 = и + V иг ),

о = 4 С "8 -[4 К&-ТО -е(т) —,

* а

о ав

о |2 -п = Б(0, и = ^(0,

здесь обозначены: V - набла-оператор [25], о - тензор напряжений, е - тензор малых деформаций, и - вектор перемещений, 8(?) - вектор внешних усилий на части границы ^ ,

ие (?) - вектор заданных перемещений на части границы Е2, п - вектор внешней нормали, 4С - тензор 4-го ранга - тензор модулей упругости, а 4К(^ - т, ) - тензор ядер релаксаций.

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

Т Т йт х = Г йт а = ехр

1 Л п (й(т\\ ' 1 ¡лШ-тЛЛ' в п + КЙ

о ав

где (0(г)) - функция температурного сдвига, а, а - константы, 0- температура, Л0 = в-в0, а 0О - начальное значение температуры.

Рассмотрим случай установившихся колебаний, когда изменения векторов внешних усилий 8(') и заданных перемещений ив (') происходит по моногармоническому закону

8

где 8*, ив* - комплексные амплитуды колебаний, о - заданная частота колебаний. Тогда будем искать перемещения и, являющиеся решением задачи (1) также в моногармоническом виде

и

Подставляя выражение (4) в (1) и (2), получаем, что деформации и напряжения также являются моногармоническими функциями, если пренебречь начальным периодом колебаний: £ = Яе(£ в10'), с = Яе(сеш'), где с , £ - комплексные амплитуды тензоров напряжений и деформаций, связанные следующими соотношениями:

г и о* = 4С*(й5)--£*. Здесь обозначен тензор комплексных модулей

упругости

4С*

о

а со = соав(6) - приведенная частота колебаний.

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

[V-с*=0

о* = 4С* (й5)-е*

1

£* = ^ (V и* +1

с* I -п = S*, и* = ив*.

Вязкоупругие характеристики изотропных компонентов композитов

Для изотропных вязкоупругих материалов тензор модулей упругости имеет 2 независимые константы: К - модуль объемного сжатия и О - модуль сдвига, тензор ядер релаксации 4 К (г) имеет 2 независимые функции - Кк ('), (') - ядра объемной и сдвиговой релаксации, причем для большинства твердых сред, объемной релаксацией можно пренебречь [5,20 ], тогда тензоры 4С и 4К(') можно представить в следующем виде [11]

I

2

4 С = (к - - а)Е  Е + 20А , 4 К(0 = Ка (0(Е  Е + 2А),

3

где А - единичный тензор 4-го ранга, Е - метрический тензор. После подстановки выражений (7) в (5), получаем

9

+

где

0*

0

-комплексный модуль сдвига.

Принимаем для ядра релаксации Ке (г) модель в виде суммы экспонент [5,11]

п

Ка ) = £ г,

где А,, г, - вязкоупругие константы материала (спектры релаксации и времен

релаксации). Тогда для комплексного модуля сдвига а* имеем следующее аналитическое выражение от приведенной частоты колебаний

N А п А зг

а' = а + га\ а = а+£ —, а" = £ ^^.

, % 1 +

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

матрицы Е* и комплексный коэффициент Пуассона у* могут быть вычислены по формулам, формально совпадающим с формулами для упругих констант [5,11]:

9Ш*

Е * =

*

у =

зК + а*' зк - 2а* 6к + 2а*'

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

композитов

Однонаправленные композиционные материалы (Ш композиты), состоящие из пучков волокон, собранных в нити, и связанных между собой полимерной матрицей, представляют собой трансверсально-изотропные материалы, даже, если матрица и волокна являются изотропными. Для вычисления компонент тензоров комплексных модулей упругости Ш композита С*и воспользуемся приближенной моделью смесевого типа [20,26,27], в которой каждая нить представляет собой систему большого числа

параллельно расположенных однонаправленных (Ш) элементов цилиндрической формы. Тогда комплексные модули упругости нити можно вычислить по следующим формулам:

Е* = ЕрРг + Е»<рг), Е*т = р +

у*=угРг рг), у* =у„

; 2

1 (1У) , (1 ~Р/)(1 + К,\-1

где Е* - продольный комплексный модуль упругости Ш композита в направлении

Г-,» *

ориентации волокон, Е» - поперечный комплексный модуль упругости нити, у* -продольный комплексный коэффициент Пуассона, у» - поперечный комплексный

Е*

коэффициент Пуассона, О* - продольный комплексный модуль сдвига, О* =-^—

2

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

продольный и поперечный коэффициенты Пуассона моноволокон, О ^ - продольный

модуль сдвига моноволокон, е* и у* - комплексные модуль упругости и коэффициент

Пуассона матрицы, р^- - относительное объемное содержание волокон.

По полученным значениям упругих констант (формулы (12)) составим компоненты

тензора комплексных упругих податливостей П*0 Ш композита:

ук1

(П*0) =

\п п п1111 П 1122 п п1133 0 0 0

П п 2222 П 1 2233 0 0 0

П *0 п3333 0 0 0

2П1313 0 0

сим. 2 П 2П 2323 0 2П*0 2П1212

1 Е; Е; Е; 0 0 0

Е; 1 0 0 0

е; 1 0 0 0

0 0 0 1 20» 0 0

0 0 0 0 1 20* 0

0 0 0 0 0 1

(13).

2О*

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

'*0 _ ( гт*0 4

уШ = (ПуШ1 ;

^ / ФА . _1

комплексных упругих податливостей: С*° = (Пш) .

Слоисто-волокнистые композиты представляют собой многослойные материалы, каждый слой которых представляет собой Ш композит, повернутый на угол фа вокруг

оси 0%3 (рис.1).

Рис.1.Ячейка периодичности слоисто-волокнистого композита с числом слоев N=3.

Тогда в повернутой (собственной) системе координат 0^г(а), ось 0Е,[а) которой совпадает с направлением ориентации волокон каждого слоя, компоненты С*0 = (П*0)— тензоров комплексных модулей упругости вычисляем по формулам (13), а в единой для всех слоев системе координат СВК О^ компоненты тензора модулей упругости а -го слоя вычисляются с помощью формул преобразования компонент тензора 4-го ранга [25]:

с:^)=с:лр^1т0а.,„0\0аи1, £ е ^,

здесь Qaim - элементы матрицы поворота слоя с номером а, эта матрица имеет следующий вид:

" COS Фа sin Фа О [Qajj ] = - sin фа COS Фа .

О 0 1

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

композитов

Рассмотрим теперь общий случай вязкоупругого композита, состоящего из N компонент, и обладающего периодической структурой. На границах раздела компонентов будем полагать выполненными условия идеального контакта

[с*] • п = 0, [и*] = 0,

где [ ] - скачок функций. Тогда для композита в целом можно сформулировать задачу

композита [21,22], тогда: 4С* = 4 С* (¿5,^). Введем малый параметр к = 1!Ь, где 1 -характерный размер ЯП, Ь - глобальный размер конструкции из ПКМ, и х - глобальные координаты, связанные с локальными координатами соотношением \ = х/к . Функции

и*(х,^) и 4С* = будем считать периодичными по \ . Введем также обозначения

д г д

для производных по х и ^ : ^х= е' -—; ^ = е1 -—, где е' - векторы базиса. Имеет место

следующее правило дифференцирования периодических функций: У = .

Обозначим (/(х, = — | /(х, - оператор осреднения по ЯП.

Тогда в соответствии с методом асимптотического осреднения [19-23 ] представим перемещения и *, являющиеся решением задачи (6), в виде асимптотического разложения по степеням малого параметра к

и*

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

£* = £*(0) +К£*(1) +К2..., О* = О*(0) +КС*(1) +К2...,

где

£*

причем в нулевом приближении имеют место следующие определяющие соотношения вязкоупругости: с*(и) = 4С*(й5,^)--е*(0).

Подставляя разложения (18) в уравнение системы (6), и собирая в каждом уравнении члены при одинаковых степенях малого параметра к, получаем из условия равенства нулю членов при наименьших степенях параметра к следующую задачу вязкоупругости "на ячейке периодичности"

V о*(0) = 0

О

= 4С*((Э, • •£*

*(0) _

1

=£+ ^

и

*

[о*(0)]• п = 0, [и*(1)] = 0, [[о*(0)]] • п = 0, [[и*(1)]] = 0,

< и*(1) >= 0.

и

*(1)Гч

здесь обозначены [[о*(0)]] • п = 0, [[и*(1)]] = 0- условия периодичности функций на границах ячейки периодичности, а < и*(1) >= 0 - условие нормировки. Формальное решение задачи (20) в силу ее линейности может быть записано в виде

1 2

3

и*

где N(5) - неизвестная тензор-функция 3-го ранга. Тогда, подставляя эти выражения в определяющее соотношение системы (22), и, осредняя его по ЯП, получаем искомые эффективные определяющие соотношения вязкоупругости

< с*

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

С» =< 4С'(а, 5) --(Д +0 3N + ^ 0 3N)) >.

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

СВК композитов

Рассмотрим теперь частный случай СВК композитов, обладающих слоистой структурой, периодической только по одному направлению Ох3. Соответствующую

локальную координату обозначим как . Характерный размер ЯП I для слоистого

композита представляет собой толщину повторяющегося пакета слоев, глобальный размер Ь - это толщина СВК, а малый параметр - к = I / Ь = 1/ N, где N - число повторяющихся пакетов слоев в составе СВК. Задача (20) на ЯП для СВК является одномерной, и в компонентной записи в декартовом базисе эта задача для слоистых композитов может быть записана следующим образом:

дс™

—= 0, дд '

*:г=с;к1(сдле;?\

йи *

2£*(0) = 2£ + 5 1 + 5 к 2ьк1 2йИ +5к3 - - + 513 - - ,

дд дд

*(1) *(1) _*(0) _*(0) /о/14

и 1 = и 1 , С 3

- 1д=-0.5 - 1д=0.5 -3 1д=-0.5 13 1д=0.5

[с*(0)] = 0, [и(1)] = 0, 5 ^

(и*(1)) = 0

л. - *(1) « о-* ди*(0) ди*(0) относительно функций и* ', здесь обозначены: 2е * = —---1-----компоненты тензора

1 дХ. дХг.

осредненных деформаций, являющиеся "входными данными задачи". Локальная

координата д принимает значения в диапазоне от -0.5 до 0.5. Подобно аналогичной задаче в теории упругих композитов [20], задача (24) допускает явное аналитическое решение, получим его несколько иным методом, отличным от предложенного в [20]. Интегрируя уравнение равновесия в системе (24), получаем = A * = const,

*(0)_r* *(0),r* *(0) *(0)_r* *(0),r* *(0) Г9<-ч

UIJ C IJk3 ck3 C IJKL CKL , ui3 C i3k3 ck3 C i3KL CKL ,

*(0)_тт* 9*(0^9F* ^UkX> | du3(1) e

bKL Ь KL, 2bk3 2fc k3 ^ ^ °k3,

од од

где А - константы решения, большие индексы 1ЦК,Ь принимают значения 2 и 3, а малые '/^/-значения 1,2,3. Выразим из второго уравнения этой системы деформации е*(0):

е/*3 ^ = Ск3г3

где Скз)3 - компоненты тензора комплексных упругих податливостей, обратного к С *ш3, и подставим в это выражение деформации £*(0) (5 -е уравнение системы (25)), тогда получим:

ди*

=

по индексу р здесь суммирования нет. Заметим, что условие периодичности (6-е уравнение системы (25)), можно записать в виде

А *

0=и„() - = —— ад =<—— >, следовательно, среднее значение от

р |д=зо_5 р \д=0_5 ^зд зд

выражения (27) должно обращаться в ноль, в результате получаем уравнение для нахождения константы А': < с*"1 ш(А 3С*лкь £*кь)-£*р1 >= 0. Отсюда получаем

А =< С^РЗ'З >з1

Подставляя это выражение в (28), находим

*

ак3 с к3'3((<с '3 р3 > <с р3 у3 с У3КЬ > 3с 13КЬ )аКЬ <с '3 р3 > ар1) '

Наконец, подставляя это выражение в определяющие соотношения (2-е и 3-е уравнения системы (25)), получаем:

— *(0) _ тт -* -*

ии Н11К1аК1 +нЦр 3ар 3,

*

и'3 11'3к3ак3 +н'3КЬаКЬ ,

тт = П* 4-С* С*1 С *з1 *з1 Г1* Г^* л^*1

НЦКЬ С ЦКЬ ^^ Цк3 С к3'3 < С ¿3р3 > < С р3у'3 С ]3КЬ > 3с Цк3 С к3'3 С '3КЬ

Н =С * С * 1 < С * 1 >з1

НЦк1 С 1Л3 С !'3'3 <с '3к3 > • У^^)

Осредняя получившееся выражение по ЯП, находим эффективные определяющие соотношения вязкоупругости слоистого композита:

< с*(0) >= С * Ё* +С * Ё*

<ии > С Цк 3 ек 3 С ЦКЬ еКЬ,

<с*(°) >= С * е* +С* Ё

<с3 > С -3к3 ек3 С -3КЬ еКЬ

а также компоненты эффективного тензора комплексных модулей упругости слоистого композита:

СЦКЬ =< НЦКЬ > = <СЦКЬ > + < СЦк3(Ск3-3) >< (С-3 р3) > < (Ср313) с,* 3КЬ > —< сщ 3(Ск3 -3) С--3КЬ >

Ц =< я*,* >=< Щ )-1 ><

Заметим, что < Ник3 >=< Нк31 > , но Н'Ц3 ^ Нк31 .

В явном виде формулы (32) имеют следующий вид:

с * =< >-1

С3333 < ^* >

С3333

- С * 1

у"|* 33Ц 1 -1 ( =< —33и ><-> 1

С3333 С3333

* * * * С * =< С * > - < С33цС33КЬ > + < С33Ц >< С33КЬ >< 1 >-1

СЦКЬ < СЦКЬ > < > + < ^* >< ^* >< ^* >

=<_С313и_ ^ _С2323_ ^^ _С1313_ ^ _

313 Ц х^* \2 \2 ^г* \ 2

С2323С1313 —

*

— <_С1323_ >2)-1

С2323С1313 -

здесь I, J, К, L = 1,2.

Результаты численных расчетов вязкоупругих характеристик матрицы,

^ композита и СВК-композита

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

20 1 + у

От = 0,345 ГПа; Кт = —тт—^, = 0.45

3

пт = 1, А = 1 ГПа, Тх = 0.01с , а = 20К, а2 = 300К ,

а также характеристики волокон, которые считались чисто упругими :

Ег = 200ГПа; ^ = 0.25 , А = 0 .

Коэффициент армирования Ш композита был принят равным 0.6. СВК-композит имел следующую структуру армирования: [0/+45/-45/90], с равным соотношением всех 4-х слоев - по 25% каждого слоя.

Были рассчитаны действительные и мнимые части компонент тензора комплексных модулей упругости: С*к1(а) = Яе(с*н(а)) +11т(С**ы(а)), а также тангенс угла потерь для модулей упругости

1т

1Еб(С;к1 (а)) =

Расчеты были проведены в диапазоне частот колебаний от 1 до 70 Гц.

На рисунках 2...4 представлены результаты расчетов вязкоупругих характеристик матрицы в зависимости от частоты колебаний. Результаты расчетов показывают, что на частоте а =9,5 Гц у матрицы существует пик (экстремум) тангенса угла потерь tgSG = tgS(Gm(а)) модуля сдвига. Значения тангенса угла потерь по сдвиговым свойствам составляют 0.75. Аналогичный пик имеется у графика тангенса угла потерь коэффициента Пуассона, но его максимальное значение составляет tgSv =0.12, которое достигается на частоте со =6,5 Гц.

Рис.2 Зависимость действительного модуля сдвига матрицы от частоты колебаний

/

/

/

!

/

12 24 36 48 во 72

ю> Гц

Рис. 3 Зависимость действительного части комплексного коэффициента Пуассона матрицы

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

;

\

\

<

/

12 Гц

Рис.4 Зависимость тангенса угла потерь комплексного модуля сдвига и комплексного коэффициента

Пуассона матрицы от частоты колебаний

На рисунках 5...8 представлены результаты расчетов вязкоупругих характеристик Ш композита в зависимости от частоты колебаний. Результаты расчетов показывают, что экстремум функции тангенса угла потерь Х%8Си= достигается при той же

частоте, что и для матрицы - при о =9,5. Однако, если для значения тангенса угла потерь в продольном направлении tg8c составляют всего 0.004, то в поперечном направлении

гё8Съъъъ составляют 0.52, а при сдвигах Щ8Сгъгъ, ^8С его значения близки к

соответствующим значениям матрицы.

12 24 36 48 60 72

Гц

Рис.5 Зависимость действительной части комплексного продольного модуля упругости С 0цц

Ш композита от частоты колебаний

'уЫ

\

\

\ С 2222 /

У \ / 2233

вР С*0

12

24

36

48

60

72

га, Гц

Рис.6 Зависимость действительной части комплексных модулей упругости С Ш композита от частоты

'Ук!

\\7 2323

х \

/ \ 4 V / улъи

1 3333 Xх

1 •—

12

24

36

48

60

72

ю, Гц

Рис.7 Зависимость тангенса угла потерь комплексных модулей упругости С ¿¿и Ш композита от частоты

колебаний

tgër,*o

уЫ

/Л

/ \ / \

\ \

1 1&5 *о 1212 Ч

о

- *

12 24 36 48 60 72

со, Гц

Рис.8 Зависимость тангенса угла потерь комплексных модулей упругости С 0уи Ю композита от частоты

колебаний

На рисунках 9...12 представлены результаты расчетов вязкоупругих характеристик СВК-композита. Результаты расчетов показывают, что на частоте а =9,5 Гц, так же, как и у матрицы, и у Ш _композита, у СВК-композита существует максимум значений тангенса угла потерь Х%8Си= tg8(C*и(«)), всех комплексных модулей сдвига. Минимальные

значения тангенса угла потерь в продольном направлении tg8c составляют всего 0.004,

в поперечном направлении tg8c составляют 0.52, а при межслоевых сдвигах tg8c ,

tg8c его значения близки к соответствующим значениям матрицы - 0.7. Тангенсы угла

потерь модуля сдвига в плоскости слоя tg8c значительно меньше - 0.06, но превышает

значения tg8c . Аналогичные пики имеется у графиков тангенса угла потерь tg8c ,

tg8c , но их абсолютные максимальные значения не превышают 0.025.

Таким образом, можно сделать вывод, что вязкоупругие свойства матрицы достаточно эффективно реализуются в составе композита СВК - но только при межслойных сдвигах и в поперечном направлении. Эта особенность накладывает специфические требования к расчету вязкоупругих свойств конструкций из ПКМ -традиционные методы расчета пластин и оболочек по теории Кирхгофа-Лява и Тимошенко не могут обеспечить точного расчета трансверсальных и межслоевых напряжений, играющих ключевую роль в механизме диссипации энергии композитов. Поэтому необходима разработка специализированных методик расчета напряжений в тонкостенных многослойных композитных конструкциях при расчете их упруго-диссипативных характеристик. Такая методика была предложена в работах [28-32] для

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

Рис. 9 Зависимость действительной части комплексного модуля упругости Ц ш СВК от частоты

колебаний

ГПа

к

оо

—----¿ Ке

*е

\

Яе

, Ке

и 12 24 36 48 60 72

о, Гц

Рис.10 Зависимость действительной части комплексных модулей упругости СВК от частоты

Рис.11 Зависимость тангенса угла потерь комплексных модулей упругости

ст свк

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

( -^ ш

/

1У

// *-3333

1/

12 24 36 48 60 72

(о, Гц

Рис.12 Сравнительные зависимости тангенса угла потерь комплексного модуля сдвига матрицы и комплексных модулей упругости СВК при межслойном сдвиге Ц313 и в трансверсальном направлении

С *

С3333

Заключение

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

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

Исследование выполнено за счет гранта Российского научного фонда (проект №1419-00847).

Список литературы

1. Hashin Z. Viscoelastic Behavior of Heterogeneous Media // J. Appl. Mech. Trans. ASME. 1965. Vol. 32, iss. 3. P. 630-636. DOI: 10.1115/1.3627270

2. Christensen R.M. Theory of Viscoelasticity. 2nd ed. New York: Academic Press, 1982. 356 p.

3. Christensen R.M. Mechanics of Composite Materials. New York: John Wiley & Sons, 1979. 324 p.

4. Ferry J.D. Viscoelastic properties of Polymers. 2nd ed. New York: John Wiley & Sons, 1979. 316 p.

5. Ильюшин А.А., Победря Б.Е. Основы математической теории термовязкоупругости. М.: Наука, 1970. 356 с.

6. Победря Б.Е., Димитриенко Ю.И. Связанные задачи линейной термомеханики деформируемых твердых тел // Успехи механики. 1987. № 2. С. 97-137.

7. Зиновьев П.А., Смердов А.А., Кулиш Г.Г. Экспериментальное исследование упруго-диссипативных характеристик углепластиков // Механика композитных материалов. 2003. Т. 39, № 5. С. 595-602.

8. Димитриенко Ю.И., Лимонов В.А. Влияние ориентации волокон на диссипативный разогрев и деформативность вязкоупругих композитов при циклическом нагружении // Механика композитных материалов. 1988. Т. 24, № 5. С. 797-805.

9. Dimitrienko Yu.I. Nonlinear Continuum Mechanics and Large Inelastic Deformations. Springer Netherlands, 2010. 722 p. (Ser. Solid Mechanics and Its Applications; vol. 174.). DOI: 10.1007/978-94-007-0034-5

10. Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 2. Универсальные законы механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 560 с.

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

12. Sheldon Imaoka. Analyzing Viscoelastic Materials // ANSYS Advantage. 2008. Vol. 2, no. 4. P. 46-47.

13. Matzenmiller A., Gerlach S. Micromechanical modeling of viscoelastic composites with compliant fiber-matrix bonding // Computational Materials Science. 2004. Vol. 29, iss. 3. P. 283-300. DOI: 10.1016/j.commatsci.2003.10.005

14. Zmindak M., Riecky D., Dudinsky M. Finite Element Analysis of Viscoelastic Composite Solids // Proc. of the 4th International conf. "Modelling of Mechanical and Mechatronic Systems 2011". 20-22 September, 2011. Herlany, Slovak Republic, 2011. P. 576-584.

15. Michel J.C., Moulinec H., Suquet P. Effective properties of composite materials with periodic microstructure: a computational approach // Comput. Methods Appl. Mech. Engrg. 1999. Vol. 172, iss. 1-4. P. 109-143. DOI: 10.1016/S0045-7825(98)00227-8

16. Shibuya Y. Evaluation of creep compliance of carbon-fiber-reinforced composites by ho-mogenization theory // JSME Int. J. Ser. A. 1997. Vol. 40. P. 313-319.

17. Haasemann G., Ulbricht V. Numerical evaluation of the viscoelastic and viscoplastic behavior of composites // Technische Mechanik. 2010. Vol. 30, no. 1-3. P. 122-135.

18. Masoumi S., Salehi M., Akhlaghi M. Nonlinear Viscoelastic Analysis of Laminated Composite Plates - A Multi Scale Approach // International Journal of Recent Advances in Mechanical Engineering. 2013. Vol. 2, no. 2. P. 11-18.

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

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

21. Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов // Математическое моделирование. 2012. Т. 24, № 5. С. 3-20.

22. Димитриенко Ю.И., Соколов А.П. Численное моделирование композиционных материалов с многоуровневой структурой // Известия Российской академии наук. Серия физическая. 2011. Т. 75, № 11. С. 1549-1554.

23. Димитриенко Ю.И., Соколов А.П. Разработка автоматизированной технологии вычисления эффективных упругих характеристик композитов методом асимптотического осреднения // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2008. № 2. С. 56-66.

24. Димитриенко Ю.И., Соколов А.П. Система автоматизированного прогнозирования свойств композиционных материалов // Информационные технологии. 2008. № 8. С. 31-38.

25. Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 1. Тензорный анализ. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 463 с.

26. Димитриенко Ю.И., Сборщиков С.В., Соколов А.П., Садовничий Д.Н., Гафаров Б.Р. Численнное и экспериментальное моделирование прочностных характеристик сферопластиков // Композиты и наноструктуры. 2013. № 3. С. 35-51.

27. Димитриенко Ю.И., Сборщиков С.В., Соколов А.П. Численное моделирование микроразрушения и прочностных характеристик пространственно-армированных

композитов // Механика композиционных материалов и конструкций. 2013. Т. 19, № 3. С. 365-383.

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

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

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

31. Димитриенко Ю.И., Федонюк Н.Н., Губарева Е.А., Сборщиков С.В., Прозоровский А.А. Многомасштабное конечно-элементное моделирование трехслойных сотовых композитных конструкций // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 10. С. 359-382. Б01: 10.7463/1014.0730105

32. Димитриенко Ю.И. Асимптотическая теория многослойных тонких пластин // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2012. № 3. С. 86-100.

Science and Education of the Bauman MSTU, 2014, no. 11, pp. 748-770.

DOI: 10.7463/1114.0734246

Received: Revised:

09.08.2014 27.10.2014

Science ^Education

of the Bauman MSTU

ISSN 1994-0448 © Bauman Moscow State Technical Unversity

Simulation of Viscoelastic Properties of Fibrous Laminated Polymer Composite Materials

Yu.I. Dimitrienko1'*, E.A. Gubareva1, dmut bmstu agmaü com

S.V. Sborschikov1, N.N. Fedonyuk2

:Bauman Moscow State Technical University, Moscow, Russia 2Krylov State Research Centre, St. Petersburg, Russia

Keywords: viscoelastic properties, composites, polymer composite materials, fibrous laminated composites, asymptotic homogenization method, complex elastic modules, tangent of loss angle

The paper presents a developed technique for calculating the viscoelastic properties of fibrous laminated polymer composites at steady cyclic vibrations. The technique uses the asymptotic homogenization procedure. It allows us to predict a complete set of complex elastic modules using rather simple analytical relations. The paper gives examples of computational modeling for viscoelastic unidirectional and fibrous laminated polymer composites. They have shown that maximum values of viscoelastic properties determined by the loss angle tangent of complex elastic modules are realized in the transverse direction towards the relative fibers orientation and in interlayer shears. It is shown that these peculiarities of viscoelastic composite properties impose some demands for stress calculation methods in thin-walled composite structures with the assumption that it is necessary to calculate transversal throw-thickness and interlayer stresses, which the classic theories of thin plates, usually, neglect.

References

1. Hashin Z. Viscoelastic Behavior of Heterogeneous Media. J. Appl. Mech. Trans. ASME, 1965, vol. 32, iss. 3, pp. 630-636. DOI: 10.1115/1.3627270

2. Christensen R.M. Theory of Viscoelasticity. 2nd ed. New York, Academic Press, 1982. 356 p.

3. Christensen R.M. Mechanics of Composite Materials. New York, John Wiley and Sons, 1979. 324 p.

4. Ferry J.D. Viscoelastic Properties of Polymers. 2nd ed. New York, John Wiley and Sons, 1979. 316 p.

5. Il'iushin A.A., Pobedria B.E. Osnovy matematicheskoi teorii termoviazkouprugosti [Bases of mathematical theory of thermoviscoelasticity]. Moscow, Nauka Publ., 1970. 356 p. (in Russian).

6. Pobedrya B.E., Dimitrienko Yu.I. Coupled problems of linear thermomechanics of deforma-ble solids. Uspekhi mekhaniki = Advances in Mechanics, 1987, no. 2, pp. 97-137. (in Russian).

7. Zinov'ev P.A., Smerdov A.A., Kulish G.G. Experimental Investigation of Elastodissipative Characteristics of Carbon-Fiber-Reinforced Plastics. Mekhanika kompozitnykh materialov, 2003, vol. 39, no. 5, pp. 595-602. (English translation: Mechanics of Composite Materials, 2003, vol. 39, iss. 5, pp. 393-398. DOI: 10.1023/B:MQCM.0000003289.12297.84 ).

8. Dimitrienko Yu.I., Limonov V.A. Effect of the orientation of the reinforcement in dissipative heating and deformability of viscoelastic composites on cyclic loading. Mekhanika kompozitnykh materialov, 1988, vol. 24, no. 5, pp. 797-805. (English translation: Mechanics of Composite Materials, 1988, vol. 24, iss. 5, pp. 595-602. DOI: 10.1007/BF00608843 ).

9. Dimitrienko Yu.I. Nonlinear Continuum Mechanics and Large Inelastic Deformations. Springer Netherlands, 2010. 722 p. (Ser. Solid Mechanics and Its Applications; vol. 174.). DOI: 10.1007/978-94-007-0034-5

10. Dimitrienko Iu.I. Mekhanika sploshnoi sredy. V 4 t. T. 2. Universal'nye zakony mekhaniki i elektrodinamiki sploshnoi sredy [The mechanics of a continuous medium. In 4 vols. Vol. 2. The universal laws of mechanics and electrodynamics of continua]. Moscow, Bauman MSTU Publ., 2011. 560 p. (in Russian).

11. Dimitrienko Iu.I. Mekhanika sploshnoi sredy. V 4 t. T. 4. Osnovy mekhaniki tverdykh sred [Continuum Mechanics. In 4 vols. Vol. 4. Bases of Mechanics of Solid Medium]. Moscow, Bauman MSTU Publ., 2013. 624 p. (in Russian).

12. Sheldon Imaoka. Analyzing Viscoelastic Materials. ANSYS Advantage, 2008, vol. 2, no. 4, pp. 46-47.

13. Matzenmiller A., Gerlach S. Micromechanical modeling of viscoelastic composites with compliant fiber-matrix bonding. Computational Materials Science, 2004, vol. 29, iss. 3, pp. 283-300. DOI: 10.1016/j.commatsci.2003.10.005

14. Zmindák M., Riecky D., Dudinsky M. Finite Element Analysis of Viscoelastic Composite Solids. Proc. of the 4th International conf. "Modelling of Mechanical and Mechatronic Systems 2011", 20-22 September, 2011, Herl'any, Slovak Republic, pp. 576-584.

15. Michel J.C., Moulinec H., Suquet P. Effective properties of composite materials with periodic microstructure: a computational approach. Comput. Methods Appl. Mech. Engrg., 1999, vol. 172, iss. 1-4, pp. 109-143. DOI: 10.1016/S0045-7825(98)00227-8

16. Shibuya Y. Evaluation of creep compliance of carbon-fiber-reinforced composites by homog-enization theory. JSMEInt. J. Ser. a, 1997, vol. 40, pp. 313-319.

17. Haasemann G., Ulbricht V. Numerical evaluation of the viscoelastic and viscoplastic behavior of composites. Technische Mechanik, 2010, vol. 30, no. 1-3, pp. 122-135.

18. Masoumi S., Salehi M., Akhlaghi M. Nonlinear Viscoelastic Analysis of Laminated Composite Plates - A Multi Scale Approach. International Journal of Recent Advances in Mechanical Engineering, 2013, vol. 2, no. 2, pp. 11-18.

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

20. Pobedria B.E. Mekhanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, MSU Publ., 1984. 336 p. (in Russian).

21. Dimitrienko Iu.I., Sokolov A.P. Multiscale modeling of elastic composite materials. Matematicheskoe modelirovanie, 2012, vol. 24, no. 5, pp. 3-20. (in Russian).

22. Dimitrienko Iu.I., Sokolov A.P. Numerical modeling of composites with multiscale microstructure. Izvestiia RAN. Seriia fizicheskaia, 2011, vol. 75, no. 11, pp. 1549-1554. (English translation: Bulletin of the Russian Academy of Sciences: Physics, 2011, vol. 75, no. 11, pp. 1457-1461. DOI: 10.3103/S1062873811110074 ).

23. Dimitrienko Iu.I., Sokolov A.P. Development of Automated Technology of Calculation of Effective Elastic Characteristics of Composites by Method of Asymptotic Averaging. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki = Herald of the Bauman MSTU. Ser. Natural science, 2008, no. 2, pp. 56-66. (in Russian).

24. Dimitrienko Iu.I., Sokolov A.P. Automated Forecasting of Composite Material Properties by Means of Homogenization Method. Informatsionnye tekhnologii = Information Technologies, 2008, no. 8, pp. 31-38. (in Russian).

25. Dimitrienko Iu.I. Mekhanika sploshnoi sredy. V 4 t. T. 1. Tenzornyi analiz [Continuum Mechanics. In 4 vols. Vol. 1. Tensor Analysis]. Moscow, Bauman MSTU Publ., 2011. 463 p. (in Russian).

26. Dimitrienko Yu.I., Sborshchikov S.V., Sokolov A.P., Sadovnichiy D.N., Gafarov B.R. Computer and experimental study modeling of failure of micro-sphere filled composite. Kompozity i nanostruktury = Composites andnanostructures, 2013, no. 3, pp. 35-51. (in Russian).

27. Dimitrienko Yu.I., Sborshchikov S.V., Sokolov A.P. Computational Modeling of Microdestruction and Strength of Multidimensional Reinforced Composites. Mekhanika kompozitsionnykh materialov i konstruktsii = Journal on Composite Mechanics and Design, 2013, vol. 19, no. 3, pp. 365-383. (in Russian).

28. Dimitrienko Yu.I., Gubareva E.A., Sborshchikov S.V. Asymptotic theory of constructive-orthotropic plates with two-periodic structures. Matematicheskoe modelirovanie i chislennye metody, 2014, no. 1, pp. 36-57. (in Russian).

29. Dimitrienko Iu.I., Iakovlev D.O. Asymptotic theory of thermoelasticity of multilayer composite plates. Mekhanika kompozitsionnykh materialov i konstruktsii = Journal on Composite Mechanics and Design, 2014, vol. 20, no 2, pp. 259-282. (in Russian).

30. Dimitrienko Iu.I., Iakovlev N.O., Erasov V.S., Fedoniuk N.N., Sborshchikov S.V., Gubareva E.A., Krylov V.D., Grigor'ev M.M., Prozorovskii A.A. Development of a multilayer polymer composite material with discrete structural-orthotropic fillers. Kompozity i nanostruktury = Composites andNanostructures, 2014, vol. 6, no. 1, pp. 32-48. (in Russian).

31. Dimitrienko Yu.I., Fedonyuk N.N., Gubareva E.A., Sborshchikov S.V., Prozorovskiy A.A. Asymptotic Theory of Viscoelastic Multilayer Thin Composite Plates. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 10, pp. 359-382. DOI: 10.7463/1014.0730105 (in Russian).

32. Dimitrienko Yu.I. Asymptotic Theory of Multilayer Thin Plates. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki = Herald of the Bauman MSTU. Ser. Natural science, 2012, no. 3, pp. 86-100. (in Russian).