Математическое моделирование процесса неизотермического отверждения полимерных композитных конструкций РКТ Текст научной статьи по специальности «Химические технологии»

ТЕХНОЛОГИЯ И ТЕХНОЛОГИЧЕСКИЕ МАШИНЫ

УДК 678.5

В. А. Тарасов, Е. В. Беляков

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

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

E-mail: [email protected]; [email protected]

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

Уровень энергомассовых характеристик ракет-носителей (РН) определяется эффективностью применения полимерных композиционных материалов преимущественно в конструкциях третьей ступени и головного обтекателя РН [1].

Наименьшая масса головного обтекателя РН достигается при использовании обшивок из полимерных композитных материалов (ПКМ) на основе углеродных волокон вследствие их высоких жесткостных характеристик при низкой плотности.

До недавнего времени в производстве конструкционных углепластиков для изделий аэрокосмического назначения широко использовались дефицитные углеродные ленты типа ЛУ-П и ЭЛУР [2] и связующее марки ЭНФБ (ТУ 1-596-36-2005.)

В последние несколько лет в производстве головных обтекателей российских РН в рамках коммерческих запусков находят применение импортные углеродные наполнители, в частности производства фирмы Porcher Industries.

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

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

Достижение поставленных перед производством целей потребовало построения ряда математических моделей [4, 7].

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

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

Процессы химической кинетики, как правило, моделируются кинетическим уравнением Аррениуса с постоянными коэффициентами:

где Еа — энергия активации отверждения, кДж/моль; Я — универсальная газовая постоянная, кДж/(моль-К); Т — абсолютная температура, К; Ф(а) — кинетическая функция; а — степень отверждения.

Кинетическая функция Ф(а) может аппроксимироваться уравнениями ви-

Из работ [8-11] следует, что эти зависимости приемлемы для описания процесса отверждения чистых связующих или наполненных инертным наполнителем [12], оказывающим незначительное влияние на кинетику и реологию процесса. Использование других армирующих наполнителей существенно влияет на кинетику и реологию процесса отверждения [13-15].

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

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

Цель настоящей работы — построение математической модели отверждения препрега на основе углеродной ткани и связующего марки ЭНФБ, позволяющей по данным динамических экспериментов прогнозировать такие параметры препрега, как степень отверждения, комплексную вязкость системы и точку гелеобразования.

Исследование кинетики отверждения препрега проводится методом дифференциальной сканирующей калориметрии с использованием дифференциального сканирующего калориметра ДСК-Д. При исследовании образцов в условиях нагрева с постоянной скоростью регистрируется тепловыделение. По зафиксированному в экспериментах значению пика тепловыделений рассчитывается степень отверждения и скорость реакции отверждения. Исследование кинетики отверждения препрега будет описываться кинетической функцией вида Ф(а) = K(1 — а)п.

Исследование вязкости препрега в процессе отверждения проводится методом динамического механического анализа в условиях сдвига с использованием вискоанализатора VA 4000 фирмы Metravib (Франция).

В рамках исследования используется препрег одной партии на основе связующего ЭНФБ и углеродной ткани арт. 3692 компании Porcher Industries. Эксперименты проводятся при нескольких различающихся температурных режимах, перекрывающих область использования модели в реальном производстве. Нагрев образцов проводится в условиях очень близких скоростей нагрева на скоростях 1, 2, 3 ° С/мин.

да [8-11]

Ф(а) = K (1 — а)п ; Ф(а) = K (1 — а)п am; Ф(а) = K (1 — а)(1 + коа).

Кинетика отверждения препрега в условиях динамического эксперимента с постоянной скоростью нагрева описывается уравнением

(1(1 Еа

— = Ае-Ш(1 - (1)

где а — степень отверждения; Ь, мин, — время отверждения; А, мин-1, — предэкспоненциальный коэффициент; Еа, кДж/моль, — энергия активации; Я = 8,31 • 10-3, кДж/(моль^К) — универсальная газовая постоянная; Т, К — температура; п — порядок реакции.

При нахождении показателя п левая часть уравнения (1) представляется в виде

(а (а (Т (а (Т (а

( = Иь • (Т = (Т • Иь = (Тв ( )

где в, °С/мин — скорость нагрева; Т, — температура.

Логарифмируя левую и правую части уравнения (1), с учетом соотношения (2) получаем

1п( =1п(А)+ п 1п(1 - а)+ (-§) Т. (3)

Для произвольного числа к динамических экспериментов с постоянной скоростью нагрева вк составляется система

'(а1 „ )\ , , „ „ ( Еа) / 1

Еа ) / 1

(К jf )) k='n<A)+»с-*1 - «'»k+(-R К TO k;

(K^))=ln(A)+"(ln(1 -<*»»+(-R) (T2)k; (4)

(Чвк))к = 1п(А)+ п(1п(1 - ар))к + (-Я)(^,

где р — число точек разбиения температурного интервала проведения экспериментов.

Из системы (4) параметры А, п и Еа определяются с помощью линейной множественной регрессии, в которой зависимой переменной является

1п^(ав), а независимыми переменными — соответственно 1п(а), 1п(1 — а) 1

и —. Т

В рамках настоящей работы для расчета параметров А, п и Еа системы (4) эксперименты по определению кинетики отверждения препрега проведены при скоростях, °С/мин: в1 = 1, в2 = 2, в3 = 3.

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

Линейная модель принимает следующий вид:

14(^в) = 22,9904 + 1,411п(1 - а) + (-10805);!. (5)

Из табл. 1 следует, что зависимость между откликом и предикторами очень сильная (Я2 > 0,75). Построенная линейная регрессия адекватно описывает взаимосвязь между откликом и предикторами, свободный член статистически значим [16].

Таблица 1

Результаты расчета по модели (4) на основе экспериментальных данных

Число точек разбиения температурного интервала проведения экспериментов

1000

Коэффициент детерминации Я^иагеё

0,98

Пред экспоненциальный коэффициент А, мин

-1

9,95 • 109

Показатель реакции п

1,41

Энергия активации Еа, кДж/моль

89,79

Анализ кинетики 36 различных партий препрегов на основе связующего ЭНФБ и трех марок углеродного наполнителя арт. 3692, арт. 4510, арт.4500 на скоростях 1... 100С/мин позволяет утверждать, что линейная модель (5) приемлема для моделирования кинетики отверждения препрега в условиях нагрева с постоянной скоростью. Отклонение расчетных значений степени отверждения от экспериментальных в области скоростей нагрева 1 0С/мин^ в ^3°С/мин не превышает 3 %, в области скоростей нагрева 3 °С/мин< в ^10 °С/мин не превышает 7 %.

Сравнительные экспериментальные данные по исследованию кинетики отверждения препрега на трех скоростях нагрева 1, 2, 3 °С/мин и расчета по модели (5) приведены на рис. 1.

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

Ж = Ф[Т (4)'аЬ.

П = П [Т, а]

(6)

Кинетика процесса отверждения препрега описывается уравнением вида

(7)

da -Ea И

— = (1 - a)n,

где А(а), мин-1, — предэкспоненциальный коэффициент; Еа(а), кДж/моль, — энергия активации.

Реология процесса отверждения препрега описывается уравнением типа

Рис. 1. Зависимость степени отверждения а от температуры Т, 0 С, при отверждении препрега на трех скоростях нагрева 1, 2, 3 0 С/мин (кривые 1, 2, 3 соответственно). Точки на графике — экспериментальные данные, линии — расчетные по предложенной модели (5)

Аррениуса-Френкеля-Эйринга [17]:

ln [n(T,a)] = ln [V(X> (a)] +

Er, (a) RT(t),

где Пга (а), Па^с, — предельное значение вязкости при Т ^ то; Еп(а), кДж/моль, — энергия активации вязкого течения.

С учетом уравнений (1), (2), (7), (8) система (6) принимает следующий

вид:

da -Ea И

—в = A(a)e"RT"(1 - a)n;

ln [n(T, a)] = ln [nx,(a)] +

En (a) RT (t)'

(9)

Большинство параметров, входящих в систему уравнений (9), не являются константами, а представляют собой функции температуры Т, степени отверждения а или времени £. Решить такую систему уравнений можно только численно.

Неизвестные функции А(а) и Еа(а) системы (9) находятся с помощью метода наименьших квадратов. Дифференциальное уравнение кинетики решается методом Рунге-Кутты четвертого порядка точности. Неизвестные функции Пга (а) и Ем (п) системы (9) находятся методом наименьших квадратов, совмещением вязкостного эксперимента, проведенного на определенной скорости нагрева, с кривой изменения степени отверждения а на той же скорости нагрева, рассчитанной по макрокинетической модели (5).

Сравнительные данные эксперимента по исследованию реологии отверждения препрега на трех скоростях нагрева и расчета значений вязкости по системе уравнений (9) приведены на рис. 2.

Точка гелеобразования Тде\ отверждающейся вязкой системы определяется при динамических вязкостных испытаниях из условия равенства единице тангенса угла механических потерь в соответствии с работой [9]. Тангенс

п'

угла механических потерь определяется как tg 5 = —- [18]. Для построения

п

Рис. 2. Зависимость комплексной вязкости п °т температуры Т, °С, при отверждении препрега на трех скоростях нагрева 1, 2, 3 ° С/мин (кривые 1, 2, 3 соответственно). Точки на графике — экспериментальные данные, линии — расчетные, по предложенной модели (9)

Ъд "

3,0 -

3,0 -2,5 -

2,0 "

1,5 -

1,0 -

0,5 -

0 -1-1-1-1-1-1- -г I I в

30 45 60 75 90 105 120 135 150 165 180 Г, °С

Рис.3. Зависимость тангенса механических потерь tg5 от температуры Т при нагревании с постоянной скоростью 1, 2, 3 0 Омин (кривые 1, 2, 3 соответственно). Точки на графике — экспериментальные данные, линии — рассчитанные по предложенной модели (10)

функциональной зависимости tg 5(а) используется математическая модель

•&(а)

tg(T,а) = С(а)ект№ , (10)

где С (а) и 'д(а) определяются аналогично функциям п^ (а) и Ем (п).

Сравнительные данные эксперимента по исследованию тангенса угла механических потерь в процессе отверждения препрега на трех скоростях нагрева и расчета значений тангенса угла механических потерь по уравнению (10) приведены на рис. 3.

Сравнительные результаты проведения экспериментов по исследованию реологии процесса отверждения препрега методом ДМА и расчета по уравнению (10) показаны в табл. 2.

Таблица 2

Анализ результатов проведения вязкостного эксперимента и расчета координаты точки гелеобразования

Параметры отверждения Экспериментальные данные Расчетные данные

Скорость нагрева,0С/мин 1 2 3 1 2 3

Температура точки гелеобразования Тде1, 0 С 121 132 137 123 133 141

Степень отверждения а 0,163 0,182 0,179 0,183 0,177 0,190

В реальном производстве композитов используются ступенчатые темпе-ратурно-временные режимы отверждения. Один из них предусматривает нагрев со скоростью 2°С/мин до температуры 60 0С, выдержку при заданной температуре 1ч, нагрев со скоростью 2°С/мин до температуры 165 °С, выдержку при заданной температуре 6 ч.

Проведен эксперимент по определению кинетики и реологии отверждения препрега по рассмотренному режиму. Сравнительные результаты проведения эксперимента и расчета по системе уравнений (9) степени отверждения а и комплексной вязкости п приведены на рис. 4.

Реализация алгоритмов решения поставленных задач организована на базе математического пакета Mathematica 7.0 компании Wolfram Research.

Рис. 4. Зависимость комплексной вязкости п и степени отверждения а от времени £ при отверждении по стандартному температурно-временному режиму. Точки — экспериментальные данные, линии — расчетные, по предложенной модели (9)

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

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

1. Конструктивн о-технологические решения силовых соединений трехслойных углепластиковых оболочек головного обтекателя ракеты-носителя / И.Г. Оленин и др. // Теория и практика технологий производства изделий из композиционных материалов и новых металлических сплавов: Тр. междунар. конф. - М.: МГУ, 2005. - С. 726.

2. Эффективность применения углеродных однонаправленных и армированных тканей при замене углеродных лент типа ЛУ- П и ЭЛУР для изделий аэрокосмического назначения / E.B. Беляков, H.B. Bымoркoв, А.К. Хмельницкий и др. // Теория и практика технологий производства изделий из композиционных материалов и новых металлических сплавов: Тр. междунар. конф. - М.: МГУ, 2007. - С. 725.

3. Autoclave vs. non-autoclave composite processing / R.Dave, J.L.Kardos, S.J. Choi // 32nd Int. SAMPE Symp. Covina, California, 1987. - Vol. 32. - P. 325337.

4. Loos A. C., S p r i n g e r G. S. Curing of epoxy matrix composites // Journal of Composite Materials. - 1983. - Vol. 17. No. 3. - P. 135-169.

5. М а т е м а т и ч е с к о е моделирование процесса отверждения изделий из полимерных композиционных материалов методом вакуумного автоклавного формования в технологическом пакете / О.С. Дмитриев, CB. Мищенко, A.B. Шаповалов и др. // Beстник ТГТУ. - 2001. - Т. 7, № 1. - С. 7-19.

6. Метод исследования параметров течения связующего при отверждении композитов / О.С. Дмитриев, CB. Мищенко, А.О. Дмитриев // Beстник ТГТУ. -2005.-Т. 11.-С. 53-60.

7. Хемовязкостная модель неизотермического отверждения термореактивных связующих / А.Г. Савин, М.А. Бокий, Н.С. Обухова // Механика композиционных материалов. - 1988. - № 5. - С. 884-887.

8. Соколов А. Д. Расчет температурно-временных режимов формования изделий из фенопластов // Методы расчета и контроля энерготехнологических параметров изготовления и переработки полимерных материалов: Материалы семинара. - М., 1982. - С. 137-142.

9. М а л к и н А. Я., К у л и ч и х и н С. Г. Реология в процессах образования и превращения полимеров. - М.: Химия, 1985. - 240 с.

10. М а л к и н А. Я., Б е г и ш е в В. П. Химическое формование полимеров. - М.: Химия, 1991.-С. 239.

11. Денисов Е. Т. Кинетика гомогенных химических реакций: Учеб. пособие для вузов. - 2-е изд., перераб. и доп. - М.: Химия, 1980. - 312 с.

12. Р о з е н б е р г Б. А., ЕниколопянН. С. Проблемы технологической монолитности изделий из композиционных материалов // Журн. Всесоюз. хим. пр-ва. - 1978. - Т. 23. - № 3. - С. 298-304.

13. З а р и н А. В., А н д р е е в А. С., В а й х а н с к и й Л. Э. Влияние армирующих химических волокон на кинетику отверждения эпоксидных связующих // Композиционные полимерные материалы. - 1985. - Вып. 24. - С. 7-10.

14. А р т е м е н к о С. В., К а ш д а ш М. М., Мальков Ю. Е. Кинетика отверждения термореактивных связующих в присутствии химических волокон // Пластич. массы. - 1988. - № 6. - С. 51-53.

15. Исследования физических и физико-химических взаимодействий в наполненных эпоксидных смолах / Яковлева Р.А. и др. // Теория и практика технологий производства изделий из композиционных материалов и новых металлических сплавов: Тр. межд. конф. - М.: МГУ, 2007. - С. 210-216.

16. Х а л а ф я н А. А. STATISTICA 6. Статистический анализ данных: Учебник. -3-е изд.: - М.: ООО "Бином-Пресс", 2008 г. - 512 с.

17. В и н о г р а д о в Г. В., М а л к и н А. Я. Реология полимеров. - М.: Химия, 1977. - 440 с.

18. М а л к и н А. Я., И с а е в А. И. Реология: концепции, методы, приложения / Пер. с англ. - СПб.: ЦОП "Профессия", 2010. - 560 с.

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

Владимир Алексеевич Тарасов родился в 1946 г., окончил МВТУ им. Н.Э. Баумана в 1969 г. Д-р техн. наук, профессор, заведующий кафедрой "Технологии ракетно-космического машиностроения" МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области технологии машиностроения, контроля и диагностики.

V.A. Tarasov (b. 1946) graduated from the Bauman Moscow Higher Technical School in 1969. D. Sc. (Eng.), professor, head of "Technologies of Rocket and Space Mechanical Engineering" department of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of technology of mechanical engineering, check and diagnostics. Евгений Владимирович Беляков родился в 1983 г., окончил МГТУ им. Н.Э. Баумана в 2006 г. Аспирант кафедры "Технологии ракетно-космического машиностроения" МГТУ им. Н.Э. Баумана. Специализируется в области производства изделий из полимерных композиционных материалов. Ye.V. Belyakov (b. 1983) graduated from the Bauman Moscow State Technical University in 2006. Post-graduate of "Technologies of Rocket and Space Mechanical Engineering" department of the Bauman Moscow State Technical University. Specializes in the field of fabrication of products from polymeric composite materials.