MODEL OF VISCOELASTIC THERMOMECHANICAL BEHAVIOR OF FIBROUS COMPOSITE AND ITS EXPERIMENTAL IDENTIFICATION
- Authors: Smetannikov O.I.1, Il’inykh G.V1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 4 (2017)
- Pages: 51-72
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/amcs/article/view/2202
- DOI: https://doi.org/10.15593/прикладная%20математика%20и%20вопросы%20управления%20/%20applied%20mathematics%20and%20control%20sciences.v0i4.2202
- Cite item
Abstract
Phenomenological determining relationships for fibrous composite materials under the conditions of thermo-relaxation (vitrification) transitions are created. The model is based on the approach used earlier by the authors to describe the thermomechanical behavior of glassy amorphous polymers with a generalization to the class of anisotropic glass materials. A technique for the experimental identification of the obtained physical model was developed and implemented. On the example of two materials - epoxy binder and fabric glass fiber reinforced plastic on its basis thermomechanical constants included in the model were obtained. In the testing experiments, the adequacy of the proposed approach was confirmed.
Full Text
Проблеме построения определяющих соотношений для волокнистых композиционных полимерных материалов в широком диапазоне температур, включающем диапазон температур релаксационного перехода (стеклования) в связующем, уделяется большое внимание исследователей в связи с ее практической и научной значимостью, особенно с точки зрения описания поведения композитов в условиях их изготовления, что зачастую сопровождается процессами отверждения связующего. Среди многочисленных исследований большое внимание уделяется особенностям деформационного поведения аморфных полимеров, использующихся в качестве связующих в композитах. В отечественной и зарубежной литературе хорошо известны модели поведения отверждающихся полимерных связующих материалов, предложенные в разное время в работах И.И. Бугакова [1, 2], В.В. Болотина [3], Р.А. Турусова [4], В.Т. Томашевского [5], B.W. Shaffer, M. Lewitsky [6, 7] и др. Общая черта этих моделей - стремление получить возможность простой оценки остаточных напряжений, без точного учета влияния релаксационных свойств и наследственных механизмов. В результате теряется существенный объем важной информации, например информации об эволюции технологических напряжений в процессе изготовления. В работах последних лет (M.C. Boyce, R.B. Dupaix [8], A.G. Varghese, R.C. Batra [9] и других авторов) предложен ряд моделей для описания механического поведения аморфных полимеров как в застеклованном состоянии, так и в температурном диапазоне, включающем температуру стеклования. Различным аспектам деформационного поведения полимеров, связанным с исследованием эффектов памяти формы, посвящено большое число публикаций, в том числе обзорного характера (работы C. Liu, H. Qin, P.T. Mather [10, 11], B. Dietsch, T.A. Tong [12], I.S. Gunes, S.C. Jana [13] и др.). Анализ публикаций показал, что большинство предлагаемых в них моделей предназначены для описания поведения полимера в достаточно узком температурном диапазоне и моделирования конкретного механизма памяти формы и не могут описывать другие важные термомеханические эффекты. В статье приводится модель термоханического поведения композиционных материалов, основанная на разработанном ранее подходе к построению определяющих соотношений при релаксационных и фазовых переходах [14-17]. Представлены результаты нескольких видов экспериментов, предназначенных для идентификации и верификации полученной модели применительно как к полимерному связующему, так и к композитному материалу. 1. Построение феноменологических определяющих соотношений для волокнистых композитных материалов в условиях терморелаксационных переходов Алгоритм построения основан на подходе, примененном ранее для описания термомеханического поведения стеклующихся аморфных полимеров с обобщением на класс анизотропных стеклующихся материалов [14-17]. 1.1. Определяющие соотношения упругого приближения для композита в режиме релаксационного перехода В композите благодаря взаимодействию волокна и связующего с различными термомеханическими характеристиками переход связующего в твердую фазу (стеклование) приводит к достаточно резкому изменению ЛКТР. Необходимо также учитывать анизотропию изменения эффективных свойств композита при переходе через интервал стеклования. Вывод определяющих соотношений основан на введении следующих упрощающих гипотез [1]: 1) характерные времена внешних воздействий намного меньше времен релаксации связующего в стеклообразном состоянии 2) характерные времена внешних воздействий значительно превышают пределы релаксационного спектра высокоэластического связующего Введение гипотез 1, 2 позволяет считать композитный материал за пределами интервала стеклования связующего упругим. Будем также полагать, что поведение полимера в интервале стеклования при уменьшении температуры обусловлено нарастанием величины механической жесткости вследствие постепенного уменьшения сегментальной подвижности и увеличения энергии межмолекулярного взаимодействия. Для количественного описания этих процессов в дальнейшем используется удельная свободная энергия упругого материала, завершенность процесса характеризуется скалярным параметром - степенью стеклования вида (1) где - параметр, определяющий ширину интервала стеклования; - температура стеклования. Обозначим через удельную свободную энергию композита с высокоэластическим связующим, через - удельную свободную энергию дополнительных связей полностью застеклованного композита. На основе сделанных предположений получим выражение удельной свободной энергии композиционного материала в виде , (2) где (3) - компоненты эффективных тензоров жесткости и коэффициентов линейного температурного расширения композита при ; - те же величины для образующихся в процессе стеклования связующего связей; - удельная энтропия; ; . После подстановки (3) в (2) получаем (4) Для определения рассмотрим поведение композита в стеклообразном состоянии . При этом степень стеклования = 1 и не меняется с изменением температуры, поэтому интегралы в правой части выражения (4) становятся постоянными величинами. Пусть при фиксированной температуре материал деформируется, в результате чего компоненты тензора деформаций получают приращение Напряжения в этом случае изменятся следующим образом: . Поскольку композит по условию полностью отвержден, коэффициент пропорциональности между и есть компонент тензора упругих характеристик стеклообразного материала Следовательно, . (5) Пусть далее при той же температуре произошел подогрев на . Считая, что материал остался недеформированным, запишем выражение для приращения напряжений из формулы (4): (6) где . Определим ту же величину через эффективные характеристики застеклованного композита: (7) где Приравнивая правые части (6), (7), получим систему в общем случае шести уравнений относительно неизвестных независимых компонент тензора . Для трансверсально-изотропных и ортотропных материалов число независимых компонентов тензора деформаций сокращается до трех: (8) Разрешая систему трех линейных уравнений (8), получим искомые значения . Если = 1, то , и свободную деформацию можно определить по формуле . Для приближенных вычислений можно принимать . Методика численного расчета входящих в формулу (4) констант по известным свойствам волокна и связующего, а также численный эксперимент по проверке адекватности данной методики описаны ниже. 1.2. Вязкоупругие определяющие соотношения для композита в режиме релаксационного перехода Будем полагать, что, в отличие от рассмотренной в п. 1.1 модели, композиционный материал в стеклообразном состоянии является вязкоупругим. При этом переход из высокоэластического в стеклообразное состояние сопровождается наращиванием межмолекулярных связей в связующем, каждая из которых в момент возникновения обладает вязкоупругими свойствами. После применения подхода, описанного в п. 1.1, с учетом термореологически простого поведения композита при температурах ниже получим следующие определяющие соотношения: • при охлаждении: , (9) где , - тензор функций релаксации композита со связующим в стеклообразном состоянии; , , - приведенные с учетом температурно-временной аналогии времена; • при нагреве: . (10) При выводе выражения (10) полагалось, что к моменту начала нагрева деформационная история «забыта» (макронапряжения полностью отрелаксировали). Объединяя формулы (9) и (10), получим , (11) где ; - время окончания охлаждения и начала нагрева. Соотношения (11) допускают обобщение на случай многократного повторения циклов охлаждения-нагрева. 2. Термомеханические испытания цилиндрических одноосных образцов и одноосная релаксация эпоксидного связующего Механические свойства композита в существенной мере определяются характеристиками связующего. Соотношения (11) в частном случае изотропного материала пригодны для описания термомеханического поведения полимера в условиях релаксационного перехода, поэтому на первом этапе была проведена серия идентификационных экспериментов для эпоксидного связующего ARALDITE LY8615/ARADUR 8615. Испытания проводились на динамическом механическом анализаторе DMA Q800 TA Instruments, который предназначен для исследования механических свойств различных материалов и имеет следующие возможности: диапазон температур - от -145 до 600 °C при скорости нагревания до 20 °C/мин. Изменение свойств образца определяется по семи экспериментальным параметрам: температуре, времени, частоте, механическому напряжению, силе, смещению и деформации. При проведении квазистатических термомеханических испытаний (TMA) используется цилиндрический образец (D = 6 мм, L = 7,8 мм). Образец нагружается при комнатной температуре постоянным весом 4,9 Н и нагревается до температуры 250 °C со скоростью 2 °C/мин. На рис. 1 показана экспериментальная температурная зависимость деформации образца под постоянной нагрузкой. Зона падения деформации соответствует релаксационному переходу в высокоэластическое состояние. Рис. 1. Экспериментальная зависимость деформации нагруженного образца от температуры, °С На рис. 2 приведен результат дополнительного эксперимента по нагреву свободного (ненагруженного) образца, результаты которого использованы для определения параметров коэффициента температурного расширения. Рис. 2. Температурная деформация образца. Точки - эксперимент, линия - аппроксимация полиномом Аппроксимация ЛКТР осуществлялась полиномом первого порядка Коэффициенты определялись методом наименьших квадратов. Найденные значения в числе прочих материальных констант модели представлены в табл. 1. Следующий вид испытаний связующего ARALDITE LY8615/ARADUR 8615 - релаксация при комнатной температуре. Эксперимент проводился на электромеханической универсальной машине с термокриокамерой Instron 5882. Результат эксперимента приведен на рис. 3. Начальное напряжение в образце - 20 МПа. Время испытания - 12 часов. Рис. 3. Функция релаксации. Точки - эксперимент, линия - расчет по модели (11) 2.1. Идентификация вязкоупругой модели связующего по результатам эксперимента Определялись материальные константы модели (11), преобразованной для одноосного напряженного состояния с функцией релаксации в виде суммы экспонент , (12) где - модуль Юнга связующего в высокоэластическом состоянии; , - мгновенный модуль застеклованного связующего; , - равновесный модуль застеклованного связующего; ; - количество экспонент функции релаксации; ; - времена релаксации; - коэффициент температурно-временного сдвига по модели Вильямса-Ланделла-Ферри, ; ; ; . Исходными данными служили результаты испытаний на релаксацию (точки на рис. 4) и термомеханические кривые (см. рис. 4). Температурная деформация, аппроксимированная ранее, исключалась из рассмотрения путем вычета из измеренной в TMA эксперименте деформации (см. рис. 4). Рис. 4. Активная (за вычетом температурной) деформация нагруженного образца при нагревании. Точки - эксперимент; линии - расчет по модели (1); скорость нагрева - 1 К/с (тонкая кривая), 2 К/с (средняя кривая) и 4 К/с (жирная кривая) Функция релаксации представлялась суммой из экспонент. Времена релаксации задавались распределенными равномерно на логарифмической шкале времени в интервале от 1 до 1018 с. Константы считались кусочно-линейной функцией от , количество вершин ломаной, распределенных равномерно на интервале [1, 18], при этом составило = 7. Определялись следующие константы: • , , , , , , и семь ординат ломаной описывающей функцию При этом вычислялся из термомеханической зависимости (см. рис. 4) по формуле • Модуль где определялся из экспериментальной кривой функции релаксации (см. рис. 3): • Аналогично по данным рис. 3. Таким образом, остается 11 неизвестных коэффициентов (вектор варьируемых параметров ): и Поиск оставшихся неизвестных констант проводился методом нелинейного программирования в пакете MatLab с использованием функции fminsearch. Решалась задача поиска минимума комбинированной невязки между экспериментальными и рассчитанными по модели (1), (12) характеристиками: где , , индексы и относятся к термомеханическому и релаксационному экспериментам соответственно; и - экспериментальные и расчетные значения в каждой точке измерения. Дополнительно накладывались ограничения на неотрицательность коэффициентов В табл. 1 представлены найденные значения материальных констант. Таблица 1 Материальные константы вязкоупругих определяющих соотношений Е1, МПа Е2, МПа Е2¥, МПа Тg, °С gL, °С Aa, °С-1 Ba, °С-2 D1 D2, °С 21 3380 119 202 8,63 17,1 0,481 10,2 39,4 На рис. 5 показана зависимость от температуры степени стеклования N(t), вычисленная по формуле (1) со значениями констант из табл. 1. На рис. 6 представлены значения найденных констант разложения функции релаксации Ci. Рис. 5. Температурная зависимость степени стеклования s Рис. 6. Константы разложения функции релаксации eст На рис. 4 расчетные кривые построены для трех различных скоростей нагрева. Видно, что рост скорости изменения температуры приводит к смещению зоны падения деформации на графиках вправо. Это соответствует экспериментальным данным, наблюдаемым для сшитых полимеров, к которым относятся эпоксидные смолы [1, 2]. В дальнейшем предполагается подтвердить полученные данные расчета в аналогичных натурных испытаниях. 3. DMA-испытания пластинчатых образцов из трехслойного тканевого стеклопластика на трехточечный изгиб с постоянной частотой и переменной температурой 3.1. Программа эксперимента Эксперимент проводился на образцах стеклопластика на основе эпоксидного связующего ARALDITE LY8615/ARADUR 8615, армированного равнопрочной стеклотканью Porcher 7781 с сатиновым 8-ремизным плетением. В качестве испытательной машины был использован динамический механический анализатор DMA Q800 TA Instruments. Образцы - 65´12,55´1,85 мм, расстояние между опорами - 50 мм. Применялось предварительное нагружение F0 = 5 Н для устранения несовершенств формы с последующим знакопеременным нагружением амплитудой Fа = 2 Н с частотой w = 1 Гц. Интервал изменения температуры - от 25 до 100 °С. Результаты измерения компонентов комплексного модуля представлены на рис. 7, а, б (точки). а б Рис. 7. Действительная (а) и мнимая (б) части комплексного модуля стеклопластика. Точки - эксперимент; сплошные линии - расчет 3.2. Идентификация термомеханических констант композита Целью DMA-испытаний в применении к разработанной физической модели является построение температурной зависимости комплексного модуля материала в диапазоне ниже интервала стеклования С учетом того что степень стеклования в указанном интервале , соотношения (12) принимают вид линейных наследственных соотношений. , (13) где , , . В дальнейших выкладках учтено, что температурная нагрузка квазистационарна (температура на протяжении нескольких циклов внешнего периодического воздействия постоянна). При гармоническом кинематическом нагружении деформацией после подстановки в (17) и математических преобразований получаем на следующий вид зависимости напряжения от времени . (14) Использовав к (14) известные соотношения для расчета компонент комплексного модуля из работы [18], получаем следующие выражения: ; (15) . (16) Для поиска входящих в выражения для модулей констант применялась методика поиска минимума комбинированной невязки между экспериментальными и расчетными значениями, аналогичная описанной в разделе 2. Функция релаксации представлялась суммой из экспонент. Времена релаксации задавались распределенными равномерно на логарифмической шкале времени в интервале от 1 до 107 с. Константы считались кусочно-линейной функцией от i, количество вершин ломаной, распределенных равномерно на интервале [1, 13], при этом составило = 5. Определялись следующие константы: , , D1, и пять ординат ломаной (k = 1, 2, …, 5), описывающей функцию , -всего девять неизвестных коэффициентов. Результаты поиска показаны в табл. 2 и на рис. 8. Таблица 2 Материальные константы вязкоупругости стеклопластика ГПа ГПа D1 °K 17,48 14,70 18,7 151 Рис. 8. Зависимость коэффициентов ядра релаксации стеклопластика Сi от времен релаксации с Из анализа рис. 7 следует, что максимальная относительная погрешность расчетных результатов не превышает 6 %. 3.3. Верификационный эксперимент. Ступенчатый трехточечный изгиб с переменной температурой Эксперимент проводился на описанных выше пластинчатых образцах. Всего было проведено два опыта. Графики изменения температуры и нагрузки во времени представлены на рис. 9. Там же показаны кривые измеренного перемещения точки приложения изгибающей силы. Рис. 9. Программа ступенчатого температурно-силового нагружения и измеренного перемещения при трехточечном изгибе пластинчатого стеклопластикового образца (результаты двух опытов) Для численного моделирования данной краевой задачи использовался пакет ANSYS Mechanical APDL. На прямоугольную область продольного сечения пластины нанесена сетка изопараметрических элементов Plane183, работающих по схеме обобщенного плоского деформированного состояния. Материал считается изотропным, линейно-вязкоупругим (модель Prony). Вводится также гипотеза о постоянстве объемного модуля, коэффициент Пуассона при нормальных условиях = 0,4. Константы модели Prony с учетом температурно-временного сдвига взяты из эксперимента. Температурная деформация не учитывается. Программа температурно-силового нагружения соответствует экспериментальной. На рис. 10 красная кривая отражает расчетную зависимость перемещений в точке приложения силы от времени. Как видно на рисунке, отличие от экспериментальных данных довольно существенно, особенно на заключительных этапах нагружения. Это объясняется накоплением погрешности эксперимента, связанной с периодическим изменением знака нагрузки, а также температурной деформацией элементов трансмиссии испытательного оборудования. Для исключения данных факторов построены кривые перемещения контрольной точки по отдельности для каждого цикла ступенчатого увеличения нагрузки относительно начального для цикла значения перемещения (рис. 11). Рис. 10. Зависимость перемещений в контрольной точке, мкм, от времени, с. Жирная линия - эксперимент; тонкая - расчет в ANSYS Для наглядности начало каждого цикла нагрузки смещено относительно предыдущего на 1000 с. Рис. 11. Поэтапные относительные перемещения контрольной точки, мкм, от времени, с. Черные линии - эксперимент; красные - расчет в ANSYS Рис. 11 можно интерпретировать как последовательность кривых ползучести при изгибе, построенных для различных постоянных во времени температур. Как видно на рисунке, в этом случае максимальная погрешность не превышает 7 %, что подтверждает адекватность предложенной физической модели. Таким образом, в представленной работе построены феноменологические определяющие соотношения для волокнистых композитных материалов в условиях терморелаксационных (стеклование) переходов. Разработана и реализована методика экспериментальной идентификации полученной физической модели. На примере двух материалов - эпоксидного связующего и тканевого стеклопластика - получены термомеханические константы, входящие в модель. Экспериментальная верификация показала адекватность предложенного подхода. В дальнейшем планируется разработка численных алгоритмов многоуровневого прогнозирования эффективных свойств композитного стеклующегося материала по свойствам компонентов, основанная на представленной физической модели.About the authors
O. Iu Smetannikov
Perm National Research Polytechnic University
G. V Il’inykh
Perm National Research Polytechnic University
References
- Бугаков И.И. Способ оценки остаточных температурных напряжений в полимерных телах // Изв. АН СССР. МТТ. - 1978. - № 3. - С. 68-74.
- Бугаков И.И. Определяющие уравнения для материалов с фазовым переходом // Механика твердого тела. - 1989. - № 3. - С. 111-117.
- Болотин В.В., Воронцов А.Н., Мурзаханов Р.Х. Анализ технологических напряжений в намоточных изделиях из композитов на протяжении всего процесса изготовления // Механика композитных материалов. - 1990. - № 3. - С. 500-508.
- Турусов Р.А. Механические явления в полимерах и композитах (в процессах формирования): дис. … д-ра физ.-мат. наук. - М., 1983.
- Tomashevskii V.T., Yakovlev V.S. Models in the engineering mechanics of polymer-matrix composite systems // International Applied Mechanics. - 2004. - Vol. 40, № 6. - P. 601-621.
- Shaffer B.W., Lewitsky M. Termoelastic constitutive equation for chemically hardering materials // Journal of Appl. Mech. - 1974. - Vol. 41, № 3. - Р. 652-657.
- Lewitsky M., Shaffer B.W. Residual thermal stresses in a solid sphere cast from a termosetting material // Journal of Appl. Mech. - 1975. - Vol. 42, № 9. - Р. 651-655.
- Dupaix R.B., Boyce M.C. Constitutive modeling of the finite strain behavior of amorphous polymers in and above the glass transition // Mechanics of Materials. - 2006. - Vol. 39, № 1. - Р. 39-52.
- Varghese A.G., Batra R.C. Constitutive equations for thermo-mechanical deformation of glassy polymers // International Journal of Solids and Structures. - 2009. - № 46. - Р. 4079-4094.
- Liu C., Qin H., Mather P.T. Review of progress in shape-memory polymers // J. Mater. Chem. - 2007. - Vol. 17. - P. 1543-1558.
- Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling / Y.P. Liu, K. Gall, M.L. Dunn, A.R. Greenberg, J. Diani // Int. J. Plast. - 2006. - Vol. 22. - P. 279-313.
- Dietsch B., Tong T. A review: features and benefits of shape memory polymers (SMPs) // J. Adv. Mater. - 2007. - Vol. 39. - P. 3-12.
- Gunes I.S., Jana S.C. Shape memory polymers and their nanocomposites: a review of science and technology of new multifunctional materials // J. Nanosci. Nanotechnol. - 2008. - Vol. 8. - P. 1616-1637.
- Труфанов Н.А., Сметанников О.Ю., Труфанов А.Н. Модели формирования полей технологических и остаточных напряжений в условиях релаксационного перехода (стеклования) // Вестник Нижегород. ун-та им. Н.И. Лобачевского. - 2011. - № 4 (5). - С. 2534-2536.
- Models of thermomechanical behavior of polymeric materials undergoing glass transition / V.P. Matveenko, O.Yu. Smetannikov, N.A. Trufanov, I.N. Shardakov // Acta Mechanica. - 2012. - Vol. 223, № 6. - P. 1261-1284.
- Сметанников О.Ю., Самусев И.В. Экспериментальная идентификация параметров определяющих соотношений для фотополимерного композита // Механика композиционных материалов и конструкций. - 2013. - Т. 19, № 1. - С. 105-116.
- Constitutive relations for viscoelastic materials under thermorelaxation transition / V.P. Matveenko, O.Yu. Smetannikov, N.A. Trufanov, I.N. Shardakov // Acta Mechanica. - 2015. - Vol. 226, № 7. - P. 2177-2194.
- ГОСТ ИСО 10112-2002. Материалы демпфирующие. Графическое представление комплексных модулей упругости. - М., 2002.
Statistics
Views
Abstract - 61
PDF (Russian) - 28
Refbacks
- There are currently no refbacks.