CC BY
122
УДК 539.376:678
ЧИСЛЕННОЕ ПРОГНОЗИРОВАНИЕ ЭФФЕКТИВНЫХ ТЕРМОВЯЗКОУПРУГИХ ХАРАКТЕРИСТИК ОДНОНАПРАВЛЕННОГО ВОЛОКНИСТОГО КОМПОЗИТА С ВЯЗКОУПРУГИМИ КОМПОНЕНТАМИ
© 2009 Е.В. Куимова, Н.А. Труфанов1
Рассматриваются вопросы численного определения эффективных термовязкоупругих характеристик однонаправленно армированных полимерных композитов по свойствам компонентов. Дана общая постановка краевой задачи термовязкоупругости на ячейке периодичности, описаны частные вычислительные эксперименты для отыскания эффективных свойств композита. Рассмотрено влияние объемной релаксации связующего на свойства однонаправленного композита. Получены эффективные обобщенные функции релаксации однонаправленного органопластика и эффективные функции температурно-временного сдвига. Подтверждено простое термореологическое поведение органопластика при наличии такого же поведения волокон и связующего.
Ключевые слова: полимерные материалы, однонаправленный волокнистый композит, эффективные свойства, термовязкоупругость, численное прогнозирование, обобщенные функции релаксации, функции темпера-турно-временного сдвига.
Введение
Волокнистые полимерные композиционные материалы (ПКМ) широко используются в машино-, авиа- и судостроении, а также электро- и радиотехнике и других отраслях. Основу волокнистых ПКМ составляют армирующие волокнистые наполнители, объединенные в монолитный композиционный материал полимерной матрицей. Одной из основных задач механики композиционных материалов является определение механических характеристик композита по свойствам его компонентов. Элементы конструкций из композиционных материалов составляются обычно из различно ориентированных однонаправленно армированных слоев, уложенных в определенной
хКуимова Елена Владимировна (ekuimova@yandex.ru), Труфанов Николай Александрович (truf2@perm.ru), кафедра вычислительной математики и механики Пермского государственного технического университета, 614000, Россия, г. Пермь, пр. Комсомольский, 29а.
последовательности, поэтому в основе методов расчета и проектирования таких конструкций лежат механические характеристики однонаправленно армированного слоя, которые требуется определить экспериментально или расчетным путем. Механическое поведение таких композитов значительно зависит от времени, что обусловлено ярко выраженными вязкоупругими свойствами полимерных связующих и некоторых типов волокон.
Проблеме прогнозирования вязкоупругих характеристик однонаправленных волокнистых ПКМ по известным вязкоупругим свойствам волокон и связующего уделено значительное число работ [1—5]. В общем случае должны быть известны из независимых экспериментов все вязкоупругие характеристики волокон и связующего. Из наиболее распространенных волокон существенными деформациями ползучести обладают органические волокна, ползучестью других типов волокон (стеклянных, углеродных, борных) можно в большинстве случаев пренебречь [1, 2]. Органическое волокно обладает трансверсальной изотропией вязкоупругих свойств, и для него требуется определение пяти независимых функций ползучести, однако на практике удается найти только характеристику продольной ползучести [3, 6], а определение трансверсальных характеристик встречает, в силу особенностей строения волокна, значительные трудности и производится в ряде работ косвенным способом [7] по найденным из эксперимента свойствам однонаправленного композита и известным свойствам связующего, или данные характеристики считаются упругими [2]. Полимерные связующие чаще всего являются изотропными вязкоупругими материалами, свойства которых определяются двумя независимыми вязкоупругими характеристиками. В ряде работ изучены свойства продольной ползучести полимерных связующих [3, 5, 8, 9], имеются данные о наличии существенной объемной ползучести [10—12].
С математической точки зрения решение задачи прогнозирования вяз-коупругих свойств однонаправленного композита осложнено наличием в постановке задачи вязкоупругости нескольких независимых вязкоупругих операторов, характеризующих свойства волокна и связующего, что делает неприменимыми традиционные методы решения краевых задач вязкоупру-гости, ориентированные на наличие одного оператора [5, 13, 14]. Во многих работах вводятся различные гипотезы, уменьшающие число независимых операторов, одной из самых распространенных является гипотеза об отсутствии объемной релаксации у полимерного связующего [5, 15]. В этом случае поведение связующего определяется одной вязкоупругой характеристикой (например, функцией продольной ползучести). Если рассматривать композиты в рамках этой гипотезы и только с упругими волокнами, то решение задачи возможно осуществить известными хорошо разработанными методами. По терминологии Б.Е. Победри [16], такие вязкоупругие композиты называют простыми. Однако влияние применения гипотезы об отсутствии объемных релаксационных свойств связующего на количественные результаты по прогнозированию эффективных вязкоупругих характеристик
не изучено. В рамках данной работы выполнены сравнительные исследования численного прогнозирования эффективных вязкоупругих свойств однонаправленных композитов с упругими волокнами с использованием и без использования гипотезы о чисто упругом поведении материала связующего в условиях действия гидростатического давления.
В случае, если в качестве наполнителя используются органические волокна, задача осложняется присутствием вязкоупругих операторов, определяющих свойства волокна. Для решения этой задачи в работе предложен вариант метода квазиконстантных операторов с частичными аппроксимациями [13, 17]. Известно, что свойства полимерных композитов существенно зависят от температуры [3, 4, 15, 18, 19]. Экспериментальные исследования свойств компонентов композита в широком температурном диапазоне приведены в работах [3, 8, 15, 20]. Во многих работах показано соблюдение принципа температурно-временной аналогии для композита и его компонентов [3, 15, 20, 21, 22], однако при решении задачи термовязкоупругости для органопластика возникает ситуация, когда функции температурно-вре-менной редукции компонент композита значительно различаются [15, 22], в связи с этим возникает проблема проверки гипотезы о термореологиче-ски простом поведении композита в целом и определения функции темпера-турно-временной редукции композита. В работе приведены результаты прогнозирования термовязкоупругих свойств органопластика в широком температурном диапазоне, в том числе и функции температурно-временного сдвига, построенные по функциям такого же сдвига его компонентов.
1. Постановка задачи
Рассматриваются однонаправленно армированные композиционные материалы, волокна которых являются упругими изотропными (стекловолокно), упругими монотропными (борное и углеродное волокно) или вязкоупру-гими монотропными (органическое волокно) материалами, а связующее является изотропным материалом, обладающим вязкоупругими свойствами (например, эпоксидное связующее ЭДТ-10).
Общая постановка краевой задачи теории термовязкоупругости для тела, занимающего область V = Vв и Vе с границей £, где Vв — область, принадлежащая волокну, Vе — область, принадлежащая связующему, включает в себя для Ь > 0 : уравнения равновесия
diva (x, t) = 0, x € V,
геометрические соотношения теории малых деформаций
(1.1)
2¿(x,t)= Vu (x, t) + (Vu (x,t))T , x € V,
(1.2)
граничные условия
и (х,*)= и (х,£), х € (1.3)
а(х,£) ■ п (х) = 0, х € Ба. (1.4)
Для замыкания системы уравнений (1.1)—(1.4) необходимо записать определяющие соотношения. Для описания поведения изотропного связующего воспользуемся линейной теорией термовязкоупругости с учетом объемной релаксации и гипотезой о термореологически простом поведении:
а (х, ¿) = И^* (в (х, ¿) - 3ас (Т (х, ¿) - Т£)) = = / & - (в (х, Т) - 3ас (Т (х, Т) - т£)) ,
4
в (х, ¿) = И§*е (х, ¿) = / (¿2 - т^е (х, т)
о
___ г ^
= /;
= /;
«т (Т (х,5))'
х € Vс, х € Vс,
(1.5)
Физические соотношения для волокна в рамках линейной теории вязко-упругости при условии термореологически простого поведения имеют вид:
а (х, ¿) = 4ИБ* ■ ■ (I (х, ¿) - аБ (Т (х, ¿) - ТБ)) =
= 0 4ЯБ (¿3 - т3) ■ ^ (I (х, т) - аБ (Т (х, т) - ТБ)) , х € VБ, (16)
= Г_£_
13 = 0 «В (Т(х,0) .
В соотношениях (1.1)—(1.6): а, I — тензоры напряжений и деформаций с компонентами а^ и соответственно; и, п — векторы перемещений и внешней единичной нормали к поверхности 5 с компонентами и^, п^; х — радиус-вектор произвольной точки тела, имеющий компоненты ж^; И — заданный на части границы 5 вектор перемещений с компонентами Ц^;
и = Б
; г, j = 1, 2, 3; ¿1, ¿2, ¿3 — время и модифицированные времена связующего и волокна соответственно; ^Б (¿), (¿) — функции объемной и сдвиговой релаксации связующего, аТ, — соответствующие функции температурно-временного сдвига связующего; 4_ЙБ (¿) — тензор функций релаксации волокна; аБ — функция температурно-временного сдвига волокна; Тя — начальная температура; То — температура приведения; индекс С характеризует параметры связующего, индекс В — параметры волокна.
Выберем декартову систему координат таким образом, чтобы ее ось О^ была направлена вдоль направления укладки волокон.
Тензор функций релаксации 4_ЙБ (¿) трансверсально изотропного волокна содержит пять независимых функций релаксации и, соответственно, должно быть пять разных модифицированных времен волокна. Однако в
литературе удается найти экспериментальные данные только для продольной вязкоупругой характеристики волокон, поэтому остальные параметры волокна будем считать упругими, а модифицированное время ¿3 в соотношениях (1.6) имеет отношение только к функции продольной релаксации вдоль оси 02.
Для упругих волокон (углеродное, стеклянное и борное) тензор 4Кв (Ь) = 4Кв не зависит от времени, соотношения (1.6) принимают вид закона термоупругого деформирования
а (х,Ь) = 4Кв ■■ (е (х,Ь) - ав (Т (х,Ь) - Т§)) , х € V, (1.7)
упругие константы, определяющие поведение материала, приведены в табл. 1.1.
Таблица 1.1
Упругие характеристики компонентов композита [15]
Волокно Е1, ГПа Е2, ГПа ^12 ^23 С12, ГПа
ЭДТ-10 3,1 - 0,42 - -
Углеродное 226 12,9 0,31 0,2 60
Стеклянное 93,2 - 0,24 - -
Борное 370 - 0,15 - -
Органическое 121 3,35 0,27 0,17 2,6
При рассмотрении материала композита в некотором малом объеме можно считать его неоднородным. Но в этом случае задача описания его свойств становится очень сложной. Для преодоления этой трудности вводится гипотеза континуума [23], которая включает процедуру осреднения, посредством нее структура и состояние материала идеализируются таким образом, что материал считается однородным, для которого присущие однородной среде характерные свойства одинаковы во всех точках. Основная задача состоит в использовании процедуры осреднения для предсказания эффективных свойств идеализированной однородной среды через свойства фаз и геометрические характеристики.
Рассматриваем однонаправленно армированный слой как неограниченную непрерывно армированную среду со структурой, образованной выпрямленными и одинаково ориентированными цилиндрическими волокнами с равными круговыми поперечными сечениями. Пространство между волокнами заполнено связующей средой — матрицей, характеристики которой отличны от характеристик волокон.
Выделим из данной среды элемент симметрии (будем рассматривать гексагональную укладку волокон) — ячейку периодичности (рис. 1.1). Предположим, что в границах ячейки градиент внешнего поля воздействий изменяется незначительно.
Рис. 1.1. Ячейка периодичности и принятая система координат (заштрихованные области — волокно): а, 6, с — ширина, высота и глубина ячейки; К — радиус волокна; Sf, Бь, Б;, Бг, Бд, — передняя, задняя, левая, правая, нижняя и верхняя грани ячейки периодичности соответственно
Относительное объемное содержание волокна в композите через параметры ячейки периодичности можно выразить соотношением: у = = уе2++2Ув = пЮОь • Для гексагональной упаковки а = -—3 •
Исходя из структуры однонаправленного композита с гексагональной укладкой волокон и согласно гипотезе континуума, материал можно считать однородной средой, которая имеет ось упругой симметрии, совпадающую с направлением волокон, т. е. является трансверсально-изотропным материалом. Таким образом, композит характеризуется пятью операторами: операторными модулями упругости Е1 и Е2, соответствующими направлениям вдоль и поперек волокон, операторными коэффициентами Пуассона и ^23 и операторным модулем сдвига С^:
Е\<р (*) = / Е\ (т) ^ (£ - т) (т, Е2^ (¿)= / Е2 (т) ^ (£ - т) (т,
—те
J 4
V 12^ (¿) = / ^12 (т) (£ - т) (т, (¿) = / ^23 (т) <£' (£ - т) (т,
— те —те
С^Р (*) = / С 12 (т) ^ (* - т) (т,
ядра которых требуется определить из решения задачи (1.1)—(1-6).
Опишем вычислительные эксперименты, которые необходимо поставить на трехмерной ячейке периодичности для нахождения вязкоупругих харак-
— те
теристик композита: продольная релаксация вдоль направления волокон, поперечная релаксация в плоскости, перпендикулярной направлению волокон, и сдвиговая релаксация вдоль направления волокон.
1.1. Продольная релаксация ячейки периодичности
Для данного случая деформирования необходимо реализовать следующие условия на границах ячейки периодичности:
х €
х € Б
я'
и (х, ¿) = 0, и2 (х, ¿) = 0, и3 (х, ¿) = 0, х € 5; а12 (х, ¿) = а1з (х, ¿) = 0 а12 (х, ¿) = а2з (х, ¿) = 0 а1з (х, ¿) = а2з (х, ¿) = 0
и1 (х, ¿) = Ц0Л (¿), х € 5/;
и>2 (х, ¿) = Ц (¿), х € ;
из (х, ¿) = Ц (¿), х € 5Г; х € 5/; х €
х € ,
где Л, (¿) — единичная ступенчатая функция Хевисайда.
Величина и0 задается, а Ц и Ц — неизвестные функции времени. Для описания механических свойств композита введем понятие средних напряжений и средних деформаций:
<ау (¿)) = и а^ (х, ¿) ¿V, ¿> 0;
у
<%• (¿)) = V / ^ (х,^) ¿V,
* > 0.
Деформация композита в продольном направлении при £ ^ 0 постоянна и запишется следующим образом:
и о
111 = <111 (¿)) = —.
о
Из решения задачи на ячейке периодичности необходимо найти следующие величины: релаксацию осредненных продольных напряжений <ац (¿)) и изменение осредненных поперечных деформаций во времени
<122 (¿)) =
<133 (¿)) =
из (¿)
Ь..... а
При этом непосредственно из решения должно следовать
<а22) = <а33) = <а12) = <а13) = <а23) = 0. Эффективные характеристики композита: ядра операторного модуля упругости Е1 (¿) и операторного коэффициента Пуассона ^12 (¿) определяются из следующих соотношений:
Е (¿) = <а11 101; ^12 (¿) = <122 1О1; ^13 (¿) = <133 (¿))/ 101В силу трансверсальной изотропии свойств непосредственно из решения должно получиться выполнение равенства ^12 (¿) = ^13 (¿).
1.2. Поперечная релаксация ячейки периодичности
Для данного случая деформирования необходимо реализовать следующие условия на границах ячейки периодичности:
х €
х €
91
Щ (х, £) = и (£) , и2 (х, £) = (£) из (х, £) = из (£), х € Sf; х € St; х € Sr.
х € Sf; х € St; х€ Sr ;
и1 (х, £) = 0, и2 (х, £) = 0, и3 (х, £) = 0, х € S1; СТ12 (х, £) = СТ13 (х, £) = 0, СТ12 (х, £) = СТ23 (х, £) = 0, СТ13 (х, £) = СТ23 (х, £) = 0,
Величина П° задается, а и и П3 — неизвестные функции времени.
Деформация композита в поперечном направлении при £ ^ 0 постоянна и запишется следующим образом:
П °
= (^22 (£)> = -
Из решения задачи на ячейке периодичности необходимо найти следующие величины: релаксацию осредненных поперечных нормальных напряжений (ст22 (£)> и изменение осредненных деформаций во времени
(ец (£)> =
П1 (*).
(е33 (£)> =
П3 (£)
с а
При этом непосредственно из решения должно следовать
(СТ11> = (СТ33 > = (СТ12> = (СТ13> = (СТ23> = 0.
Эффективные характеристики композита: ядра операторного модуля упругости Е2 (£) и операторного коэффициента Пуассона ^23 (£) определяются из следующих соотношений:
Е (£) = (СТ22 (£)>/е°2; ^23 (£) = (е33 (£)>/е°2; (£) = (еп (£))/е°2.
В силу трансверсальной изотропии свойств непосредственно из решения должно получаться выполнение операторного равенства ^2^2 = ^ЕЦ или
г г
J ^21 (£ - т) Е (т) ^т = J ^12 (£ - Т) Е (т) ^т.
1.3. Сдвиговая продольная релаксация ячейки периодичности
Для реализации условий продольного сдвига на ячейке периодичности необходимо задать следующие граничные условия:
«1 (х, = в1уЛ (¿), «2 (х, ¿) = Ц (¿) , и1 (х, ¿) = и2 (х, ¿) = 0, «3 (х, = Ц3 (¿), и3 (х, ¿) = 0,
а23 (х, ¿) = 0;
= 0; х € 5/, ,
= 0; х €
х € 5Я,
= 0; х € 5Г,
= 0; х € 5.
направлении при t
постоянна и запишется следующим образом:
^12 = <^12 (¿)) = в1/2-Из решения задачи на ячейке периодичности необходимо найти следующие величины: релаксацию осредненных касательных напряжений <а12 (¿)) и изменение осредненных поперечных деформаций во времени <^22 (¿)) = ^; <^33 (¿)) = ^.
При этом непосредственно из решения должно следовать
<ац) = <а22) = <а33) = <а13) = <а23) = 0.
Ядро эффективного операторного модуля сдвига ^12 (¿) определяется из следующего соотношения:
«12 (¿) = <а12 (¿))/в1-Непосредственно из решения должно получиться выполнение равенства и2 (¿) = и3 (¿) = 0, так как в силу трансверсальной изотропии свойств нет связи между сдвиговыми и линейными деформациями.
2. Учет объемной релаксации связующего
Во многих работах (например, [5, 15]) вводится гипотеза об отсутствии объемной релаксации связующего, упрощающая постановку задачи (1.1)—(1.6). Исследуем, какой вклад вносит учет объемной релаксации связующего в вязкоупругие характеристики композита.
Будем рассматривать однонаправленно армированные композиционные материалы, волокна которых являются упругими изотропными (стекловолокно) или упругими монотропными (углеродное волокно) материалами, а связующее является изотропным материалом, обладающим вязкоупруги-ми свойствами (эпоксидная смола) (задача (1.1)—(1.5),(1.7)). Для описания вязкоупругих свойств связующего использовались результаты работы [11], в которой приведены данные для кривых объемной податливости при различных температурах. Для расчета выбрана температура Т = 303 К. Данные по ползучести в продольном направлении были взяты из работы [3] (рис. 2.1). На рисунке сплошной линией показана функция объемной релаксации В (¿), пунктирной линией — функция продольной релаксации Е (¿).
Е(Ч, Па ВО), Па 7.Е+09
6.Е+09
5.Е+09
4.Е+09
З.Е+09
2.Е+09
1.Е+09
0.Е+00
1 10 100 1000 10000 100000
Рис. 2.1. Вязкоупругие свойства эпоксидного связующего: — функция объемной релаксации В (£); - - функция продольной релаксации Е (£))
Постановка данной краевой задачи (1.1)—(1.5), (1.7) на ячейке периодичности содержит два независимых вязкоупругих оператора. Для ее решения был использован приближенный метод — метод квазиконстантных операторов, основные положения которого для кусочно-однородных вязко-упругих сред подробно изложены в работе [13]. В соответствии с теорией метода при постоянных во времени внешних воздействиях решение задачи вязкоупругости с достаточно высокой точностью можно получать, заменяя в постановке задачи вязкоупругие операторы соответствующими функциями релаксации или ползучести, что сводится, по сути, в нашем случае к построению зависимости решения трех краевых задач теории упругости (одноосное растяжение вдоль направления волокон, одноосное растяжение в плоскости, перпендикулярной направлению волокон, и сдвиг вдоль направления волокон) от переменных во времени упругих характеристик компонентов.
Будем использовать методику определения эффективных свойств композита по свойствам его компонентов на ячейке периодичности методом конечных элементов (подробное описание этой методики приведено в работе [5]). Предварительно эта методика была оттестирована на задачах по прогнозированию упругих характеристик композитов. Результаты решения таких задач представлены в табл. 2.1.
Погрешность численного прогнозирования, определяемая как относительная к среднестатистическому значению соответствующей экспериментальной характеристики, приводилась только в случаях, когда результат расчета выходил за рамки доверительного интервала результатов эксперимента. Наблюдается достаточно хорошее соответствие экспериментальных и расчетных упругих характеристик за исключением коэффициента Пуассона
ч
ч
^23 углепластика, что, возможно, связано с недостаточно точным определением трансверсальных упругих характеристик углеродного волокна.
Таблица 2.1
Упругие характеристики композитов. Сравнение с данными экспериментов
Е1, ГПа Е2, ГПа ^12 ^23 С12, ГПа
Исходные данные [23]
Органи-
ческое 127 3,35 0,26 0,17 3
волокно
Углеродное
волокно 226 12,9 0,31 0,2 60
Связующее
ЭДТ-10 3,1 0,42
Данные эксперимента [23]
Органо- 50,5 3,63 0,37 0,46 1,58
пластик 0,4 ±5, 7 ±0, 44 ±0, 06 ±0, 02 ±0,19
Угле- 87,3 5,88 0,38 0,46 3,82
пластик 0,4 ±8, 6 ±0, 76 ±0, 03 ±0, 03 ±0,11
Результаты расчета
Органо-
пластик 0,4 52,70442 3,64932 0,36645 0,50215 1,66667
Погрешность
расчетов 4,8 %
Угле-
пластик 0,4 92,31620 5,84279 0,37536 0,58269 3,20837
Погрешность
расчетов 20,1 % 13,1 %
По описанным параметрам компонентов композита проведено два типа расчета: с учетом и без учета объемной релаксации связующего. Результаты расчетов для углепластика приведены на рис. 2.2. Рассчитана относительная погрешность которая определяется следующим образом:
¿X =
Xь - X
X "
100%,
где X — расчетное свойство композита, полученное при расчете без использования гипотезы об отсутствии объемной релаксации; X'" — расчетное свойство композита, полученное при расчете с учетом объемной релаксации связующего.
ЕЩ Па 1.15Е+11
1.14Е+11
1.14Е+11
10000 100000
С12(Ц, Па
5.0Е+09 т
4.5Е+09
4.0Е+09
3.5Е+09
3.0Е+09
2.5Е+09
2.0Е+09
1.5Е+09
1.0Е+09
5.0Е+08
0.0Е+00
I, Ч
10000 100000
Рис. 2.2. Вязкоупругие характеристики углепластика при температуре Т = 303 К (— с учетом объемной релаксации; - - без учета объемной релаксации): а - функция продольной релаксации (¿), Па; б - функция поперечной релаксации Е2 (¿), Па; в - ядро операторного коэффициента Пуассона в плоскости вдоль направления волокон ^12 (¿); г - ядро операторного коэффициента Пуассона в плоскости перпендикулярной направлению укладки волокон ^23 (¿); д - функция сдвиговой релаксации 012 (¿), Па
По проведенным расчетам можно сделать следующие выводы. Расчет функций продольной, поперечной и сдвиговой релаксации простых композитов с использованием гипотезы об отсутствии объемной релаксации связующего приводит к незначительной погрешности, которая не превышает 1 %. При определении ядер операторных коэффициентов Пуассона следует учитывать, что введение данной гипотезы качественно меняет картину данных характеристик (при отсутствии объемной релаксации ядра операторных коэффициентов монотонно растут во времени, при наличии — имеется зона убывания функций времени с последующим более быстрым ростом), что, однако, не сильно влияет на их численные значения, и максимальная погрешность составляет порядка 3 %.
з. Эффективные термовязкоупругие характеристики органопластика
Рассмотрим решение краевой задачи вязкоупругости (1.1)—(1.6) применительно к ячейке периодичности органопластика, считая волокно вязко-упругим при растяжении в продольном направлении, принимая остальные характеристики волокна упругими. Связующее (эпоксидная смола ЭДТ-10) является вязкоупругим изотропным материалом, при этом используется гипотеза об отсутствии объемной релаксации.
Таким образом, в задаче имеется два различных вязкоупругих оператора. Для решения воспользуемся приближенным методом квазиконстантных операторов [24], погрешность которого мала, если для ядер вязкоупругих операторов, входящих в решение, выполняются условия квазиконстантности. Для ядра Ев (t) оператора продольной релаксации Ев* органоволокна показатель квазиконстантности равен вв = 0,13, а для ядра Eс (t) оператора продольной релаксации EC* связующего ЭДТ-10 показатель квазиконстантности равен вв = 0,324 [24]. Существующие теоретические оценки [24] позволяют утверждать, что относительная погрешность решения по методу квазиконстантных операторов в этом случае гарантированно не превысит 10 % (величина порядка квадрата максимального по величине показателя квазиконстантности). При реализации описанных краевых задач на ячейке периодичности внешние воздействия представляют собой мгновенно приложенные в начальный момент времени и неизменяющиеся далее величины, в этом случае решение по методу квазиконстантных операторов по форме соответствует так называемому методу переменных параметров упругости
и, например, для компонент тензора напряжений имеет вид [24]:
ay (x, t) = aj [x, EB (t), EB, Gf2, vB2, vB3, ec (t), BC] +
+O (вВ) +0(вввс) + о (во,
(e)
где под aj понимается решение соответствующей упругой задачи.
Рассмотрим также применение метода квазиконстантных операторов с частичными аппроксимациями [17], который позволяет увеличить точность при наличии в решении оператора с относительно большим показателем квазиконстантности (в данном случае показатель квазиконстантности связующего).
Решение упругой задачи на ячейке периодичности имеет вид:
(х, *) = 4е) (х,Е?, ЕВ, с?2> ^3, Вс, ) Н (г), (3.1)
где для удобства использована упругая константа
шс = 2Сс!3Вс = (9Вс/(2Ес) - 1/2)-1.
Введем в рассмотрение систему дробно-рациональных функций аргумен-
с
та ис:
= (^с) 1, ^2 = 1, ^3 = , ^п = 1+131пШа , п = 4,Ж, (3.2)
где вп — набор положительных чисел.
Выбор системы базисных функций в виде (3.2) соответствует методу аппроксимаций [14, 16]. Аппроксимируем каждую компоненту тензорной функции ¿т^ (х, Е^, Е^, С12, иВ2, ивз, Вс, шс) конечным рядом по системе функций (3.2), тогда
N
4е) (х, г) * ^ ВЩ7 (х, ЕВ, ЕВ, С?2, < *23, Вс) ^п (ис), (3.3)
п=1
где для каждой фиксированной точки х = хо € V, Вщ7 — подлежащие определению коэффициенты аппроксимации.
Используя представление аналитической функции В^7 в виде интеграла Коши по аргументу Е^ и применяя принцип Вольтерры к соотношениям (3.3) с учетом (3.1), запишем вязкоупругое решение в форме
N
£
п=1
¿7 к (х, г) =
(2пг) У (г-Е**)
^п (ис) н (г),
(3.4)
где г — мнимая единица; 7 — замкнутый контур в комплексной плоскости, представляющий собой, например, окружность радиуса г с центром ЕВ. Контур выбирается таким образом, чтобы он располагался в области аналитичности функций В7 по соответствующему переменному Е^ и охватывал точку г = ЕВ . Функции параметра г в (3.4) являются аналитическими продолжениями по г упругих решений (3.3).
После проведения операций с квазиконстантным оператором (г - ЕВ*)-1 [13] получим
1
N t
El
n=1 0
j (x, t)
= Zf FT [x, Ef (t), Ef> Gf2> vf v£> Bp] | [p„ (wp*) h (t)] dr + (3.5)
+ o (ef [ + 0(вв вс) + O (ep [ .
Для реализации описанного метода в виде (3.3) необходимо решить два вопроса. Во-первых, требуется определить коэффициенты разложения F^. Точное представление упругого решения в форме (3.3) возможно лишь для узкого круга простейших задач. В остальных случаях приходится прибегать к методу численной реализации упругого решения [14], который в данном случае позволяет для некоторого фиксированного r = rm отыскать коэффициенты F^ из системы алгебраических уравнений
j) [x, Ef (Tm) > E22> Gf2, Vf2, V2|> BC> WlC] =
N _
= £ FjmPn (wp), l = 1>L,
n=1
wp — набор чисел.
с с где ш2, ш2 ,
Вычисление интеграла в (3.5) удобно вести численно с использованием квадратурных формул Гаусса, в которых тт — узлы интегрирования, тогда система (3.4) должна быть решена столько раз, сколько узлов содержит квадратурная формула.
Во-вторых, для реализации (3.5) необходимо расшифровать выражения (ш2*) Н (*) и (ш2*) Н (¿). При
(t) = f wp (t - r) p (r) dr, wp (t) = wp 1 - J rg (r) dr
выражение (3.5) приобретает вид:
<ту (x> t) = ^ + F2j (x> t) + Fj (x> t)wp + E +
t
+ /
0
(x,T) Kp (
Kg(t - r) - Fj(x>r)wprg(t - r) + E
N
n=4 Fn (x,t)
n=4
Kp (t - r)
dr,
где К,2 (¿) — ядро сдвиговой ползучести материала связующего; Г, (¿) — ядро сдвиговой релаксации материала связующего; (¿) — резольвента
ядра Гр впШ/(в„ + ш).
Приведем результаты расчетов методом квазиконстантных операторов
и методом квазиконстантных операторов с частичными аппроксимациями (рассматриваются моменты времени 1250, 2500, 5000, 10000 и 20000 ч) в сравнении с данными эксперимента [2].
Рассматривается однонаправленный органопластик с объемным содержанием волокон 50 % и их гексагональной укладкой. Необходимые для рас-
t
t
с
ш
n
чета характеристики упругости приведены в табл. 2.1. Для описания вязко-упругих свойств органических волокон использовались результаты работы [15], для ЭДТ-10 использованы данные по ползучести работы [3]. Результаты расчетов методом квазиконстантных операторов и методом квазиконстантных операторов с частичными аппроксимациями в сравнении с данными эксперимента приведены на рис. 3.1.
Е1(Ц, Па
1.Е+11 9.Е+10 8.Е+10 7.Е+10 6.Е+10 5.Е+10 4.Е+10 З.Е+10
2.Е+10 1.Е+10 0.Е+00
Ми12(Ц 0,39 0,38 0,37 0,36 0,35 0,34 0,33 0,32 0,31 0,30
<— С А- *->
5000 10000 15000
I, Ч 20000
5000
- Ъ Ч
10000 15000 20000
в С12(Ч, Па 2.0Е+09 1.8Е+09 1.6Е+09 1.4Е+09 1.2Е+09 1.0Е+09 8.0Е+08 6.0Е+08 4.0Е+08 2.0Е+08 0.0Е+00
Е2(Ц, Па 4.0Е+09 3.5Е+09 3.0Е+09 2.5Е+09 2.0Е+09 1.5Е+09 1.0Е+09 5.0Е+08 0.0Е+00
№23(1) 0,80 0,70 0,60 0,50 0,40 0,30 0,20 0,10 0,00
Тг-
А, ' АА А
5000
1,4
10000 15000 20000
5000 10000 15000
20006 4
5000 10000
д
15000
Рис. 3.1. Обобщенные вязкоупругие характеристики органопластика:---
метод квазиконстантных операторов; х - метод квазиконстантных операторов с частичными аппроксимациями; ▲ - результаты эксперимента [2]; а - функция продольной релаксации Е1 , Па; б - функция поперечной релаксации Е2 (¿), Па; в - ядро операторного коэффициента Пуассона в плоскости вдоль направления волокон ^12 (¿); г - ядро операторного коэффициента Пуассона в плоскости, перпендикулярной направлению укладки волокон ^23 (¿); д - функция сдвиговой релаксации 012 (£), Па
1п ат
Рис. 3.2. Функция температурно-временной редукции для вязкоупругих характеристик органопластика 1п (ау) и для Е1 (£) композита — о; для Е2 (£) и ^12 (£) композита — □; для ^12 (£) и ^23 (£) композита — Д; для Е (£) связующего — -; для Е1 (£) волокна — - - -
В качестве исходных данных для расчетов использовались обобщенные функции продольной релаксации волокон и связующего, функции темпера-турно-временной редукции ат которых различны. В результате получены пять эффективных обобщенных функций композита, для которых значения ат неизвестны. Для определения проведена серия расчетов для различных значений температуры и построены зависимости свойств композита в полулогарифмической системе координат, затем путем горизонтальных сдвигов определены значения 1п ат для каждой вязкоупругой характеристики композита. Полученные данные приведены на графике вместе с функциями температурно-временной редукции связующего и волокон (рис. 3.2).
Заключение
Совместное применение метода квазиконстантных операторов и метода конечных элементов удобно в практическом использовании и дает достаточно точный численный алгоритм для прогнозирования эффективных вязко-упругих характеристик однонаправленных композитов.
По результатам вычислительных экспериментов можно сделать следующие выводы: однонаправленный органопластик проявляет термореологи-чески простое поведение, при этом функция температурно-временной редукции для обобщенной функции продольной релаксации органопластика практически совпадает с функцией для органических волокон, а функции температурно-временной редукции для остальных обобщенных характеристик практически совпадают с функцией связующего.
Литература
[1] Плуме Э.З. Сравнительный анализ ползучести однонаправленных композитов, армированных волокнами различного типа // Механика композитных материалов. 1992. № 4. С. 557—566.
[2] Максимов Р.Д., Плуме Э.З. Ползучесть однонаправленно армированных полимерных композитов // Механика композитных материалов. 1984. № 2. С. 215—223.
[3] Максимов Р.Д., Плуме Э.З. Прогнозирование ползучести однонаправ-ленно армированного пластика с термореологически простыми структурными компонентами // Механика композитных материалов. 1982. № 6. С. 1081—1089.
[4] Кочетков В.А. Прогнозирование термического деформирования слоистых гибридных композитов с учетом термовязкоупругих свойств связующего и волокон // Механика композитных материалов. 1993. № 3. С. 317—323.
[5] Термомеханика вязкоупругих материалов в условиях релаксационных и фазовых переходов / Н.А. Труфанов [и др.] // Региональный конкурс РФФИ-Урал. Результаты научных исследований, полученные за 2005 г. Аннотационные отчеты. Екатеринбург: ПНЦ УрО РАН. 2006. С. 93—97.
[6] Зелин В.И., Янсон Ю.О. Определение ядер ползучести по результатам кратковременных испытаний // Механика полимеров. 1977. № 6. С. 972—975.
[7] Труфанов Н.А., Сметанников О.Ю. Приближенное определение транс-версальных вязкоупругих свойств органоволокна в составе однонаправленного органопластика // Численное моделирование статического и динамического деформирования конструкций. Свердловск: УрО АН СССР. 1990. С. 114—118.
[8] Максимов Р.Д., Плуме Э.З. Длительная ползучесть органопластика // Механика композитных материалов. 2001. № 4. С. 435—450.
[9] Уржумцев Ю.С. Прогнозирование длительного сопротивления полимерных материалов М.: Наука, 1982. 222 с.
[10] Дзене И.Я., Ректиньш М.Ф. Объемные деформации полимерных смол в условиях кратковременной ползучести // Механика композитных материалов. 1986. № 6. С. 1114—1117.
[11] Гольдман А.Я. Прогнозирование деформационно-прочностных свойств полимерных и композиционных материалов Л.: Химия, 1988. 272 с.
[12] Малмейстер А.А., Янсон Ю.О. Прогнозирование релаксационных свойств эпоксидного связующего ЭДТ-10 при сложном напряженном состоянии // Механика композитных материалов. 1983. № 5. С. 889-894.
[13] Методы прикладной вязкоупругости / А.А. Адамов [и др.]. Екатеринбург: УрО РАН, 2003. 411 с.
[14] Победря Б.Е. Численные методы в теории ползучести и пластичности. М.: Изд-во МГУ, 1981. 343 с.
[15] Максимов Р.Д., Кочетков В.А. Прогнозирование термического деформирования гибридных композитов с вязкоупругими компонентами // Механика композитных материалов. 1989. № 6. С. 969—979.
[16] Победря Б.Е. Механика композиционных материалов. М.: Изд-во МГУ, 1984. 336 с.
[17] Труфанов Н.А. Применение частичных аппроксимаций в методе квазиконстантных операторов // Вестник ПГТУ. Полимерные материалы. 1997. № 3. С. 86—89.
[18] Ташкинов А.А., Шавшуков В.Е. Тепловое расширение однонаправленных и пространственных ортогонально армированных волокнистых композитов // Механика композиционных материалов и конструкций. 2003. № 2. С. 133—141.
[19] Федоровский Г.Д. Исследование термовязкоупругих характеристик намоточных стеклопластиков в поперечном направлении // Механика композитных материалов. 1983. № 4. С. 713—718.
[20] Киселев А.Д., Кочетков В.А. Ползучесть эпоксидного связующего при термоциклировании // Механика композитных материалов. 1992. № 4. С. 557—566.
[21] Янсон Ю.О., Дмитриенко И.П., Зелин В.И. Прогнозирование деформаций ползучести однонаправлено армированного органопластика по результатам квазистатических испытаний // Механика композитных материалов. 1983. № 4. С. 610—613.
[22] Анискевич К.К. Оценка функции температурно-временной редукции ЭДТ-10 по результатам дилатометрических испытаний // Механика композитных материалов. 1989. № 4. С. 737—753.
[23] Ванин Г.А. Микромеханика композиционных материалов. Киев: Наук. думка, 1985. 304 с.
[24] Малый В.И., Труфанов Н.А. Метод квазиконстантных операторов в теории вязкоупругости анизотропных нестареющих материалов // Изв. АН СССР. МТТ. 1987. № 6. С. 148—154.
Поступила в редакцию 19/11/2009.
в окончательном варианте — 19/11/2009.
THE NUMERICAL PREDICTION OF EFFECTIVE THERMOVISCOELASTIC PROPERTIES OF UNIDIRECTIONAL FIBER COMPOSITE WITH THE VISCOELASTIC COMPONENTS
© 2009 E.V. Kuimova, N.A. Trufanov2
The problems of numerical identification of effective thermoviscoelas-tic properties of the unidirectional reinforced polymeric composite on the base of components properties are considered. General statement of a boundary value problem of thermoviscoelasticity for a cell of periodicity is yielded; special computing experiments for searching effective properties of a composite are featured. The effect of a volumetric relaxation of the binder on the properties of the unidirectional composite is considered. Effective generalized functions of a relaxation of unidirectional organoplasty and effective functions of temperature-time shift are gained. It is confirmed that organoplasty shows simple thermorheological behaviour at presence of simple thermorheological behaviour of fibers and the binder. Simple thermorheological behaviour of organoplasty at presence of simple thermorheological behaviour of fibers and the binder is confirmed.
Key words and phrases: polymeric materials, unidirectional fiber, composite, effective properties, thermoviscoelasticity, numerical prediction, generalized functions of a relaxation, functions of temperature-time shift.
Paper received 19/77/2009. Paper accepted 19/77/2009.
2Kuimova Elena Vladimirovna (ekuimova@yandex.ru), Trufanov Nikolay Aleksandrovich (truf2@perm.ru), Dept. of Numerical Mathematics and Mechanics, Perm State Technical University, Perm, 614000, Russia.
