CC BY
49
УДК 678.8
В.А. Гончаров1, А.Е. Раскутин1
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ИНФУЗИИ ПРИ ИЗГОТОВЛЕНИИ КОМПОЗИТНОГО АРОЧНОГО ЭЛЕМЕНТА
Описывается последовательность операций для моделирования технологического процесса изготовления композитного арочного элемента методом инфузии. Отражены феноменологические процессы, протекающие при формовании, а также характеристики, которые необходимо учитывать при имитации.
Ключевые слова: моделирование процесса инфузии, проницаемость, метод конечных элементов, вязкость.
The article contains the sequence of operations for modeling of manufacturing by infusion composite arched element. Described are phenomenological processes arising from the formation and critical characteristics that must be taken into account in the simulation.
Keywords: modeling of infusion, permeability, finite element method, viscosity.
Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт авиационных материалов» Государственный научный центр Российской Федерации [Federal state unitary enterprise «All-Russian scientific research institute of aviation materials» State research center of the Russian Federation] E-mail: admin@viam.ru
Введение
Одним из основных преимуществ ПКМ является то, что материал, технология и конструкция создаются одновременно - этим определяется высокая степень инноваций на всех этапах жизненного цикла материала, что требует применения новых подходов к проектированию и созданию конструкций из полимерных композиционных материалов (ПКМ) на основе сложных математических моделей, в результате использования которых появляется возможность спрогнозировать не только свойства ПКМ и конструкций из них, но и возможное образование дефектов формования [1-4].
Можно привести множество удачных примеров использования композиционных материалов в конструкциях современных изделий, таких как Boeing 787 и Space Shuttle NASA. Кроме авиационно-космической отрасли данные материалы достаточно давно используются в строительстве, автомобилестроении и судостроении.
В настоящее время ведутся разработки, направленные на снижение стоимости изготовления композиционных материалов, путем исключения необходимости использования автоклавного оборудования и изготовления препрега.
Альтернативной технологией формования ПКМ является технология вакуумной инфузии. Точность воспроизведения теоретического контура композиционных панелей, изготавливаемых методом инфузии, определяется формой оснастки, на которой осуществляется пропитка и выкладка. Основная проблема технологии вакуумной инфузии заключается в трудности получения детали со стабильными геометрическими и физико-механическими характеристиками [5-11].
Материалы и методы
С учетом процессов, протекающих при формовании инфузионными методами, определены требования, предъявляемые к исходным данным, для реализации компьютерного моделирования в программе PAM-RTM. К ним относятся характеристики связующего и наполнителя - физические, реологические и теплофизические:
- для наполнителя - градиент проницаемости, плотность, теплоемкость и теплопроводность, поверхностная площадь, а также начальная толщина преформы;
- для связующего - функция вязкости от температуры, теплоемкость и теплопроводность, плотность.
Необходимо отметить, что для наполнителя необходимо учитывать изменение проницаемости в результате процесса предварительной формовки и связанного с этим смещения его волокон относительно первоначального положения [12, 13].
Моделирование процесса инфузии сводится к решению задач гидродинамики и теплообмена. Решение задачи гидродинамики при моделировании описывается законом Дарси, определяющим расход однородной жидкости через пористую среду при ламинарном режиме потока. Основой для выполнения закона Дарси является обеспечение закона сохранения масс с учетом несжимаемости жидкости. Поскольку масса связующего в выделенном объеме со временем не меняется, то полный поток этого вектора через замкнутую поверхность обязательно равен нулю [14].
Параметрами, влияющими на процесс инфузии, являются проницаемость наполнителя, вязкость связующего, градиент давления. Чем больше проницаемость наполнителя и скорость подачи связующего, тем быстрее произойдет полная пропитка армирующего материала. Обратное влияние на процесс оказывает вязкость связующего [15, 16].
КДР
v=-
ц-ДЬ
где v - скорость пропитки, м/с; K - проницаемость армирующего наполнителя, м2; ц - вязкость эпоксидного связующего, Па с; Ь - габаритный размер зоны пропитки, м; ДР - градиент давления при пропитке, Па.
Результаты
Рассмотрим исходные данные, необходимые для моделирования процесса пропитки при изготовлении оболочки арочного элемента методом инфузии в программе СЛБ-анализа РАМ-ЯТМ:
По связующему
- Уш - объемное содержание в ПКМ;
- вязкость:
- ц=соп81 - постоянная вязкость;
- ц=А(т) - изменение вязкости по времени.
Следует отметить, что если вязкость связующего за время пропитки изменяется незначительно, то в программе РАМ-ЯТМ при вводе исходных данных она может задаваться константой. В ином случае данная характеристика описывается исключительно функциональной зависимостью от времени, полученной по результатам экспериментальных исследований с помощью вискозиметра (рис. 1).
Рис. 1. Изменения вязкости связующего ВСЭ-15 при нагреве в изотермическом режиме при температурах 80, 100 и 120°С
По армирующему наполнителю
- - объемное содержание в ПКМ;
- проницаемость пакета наполнителя:
- К1 - коэффициент проницаемости в направлении основы;
- К2 - коэффициент проницаемости в направлении утка;
- К3 - коэффициент проницаемости в трансверсальном направлении.
По распределительной сетке (РС) с повышенной филътруемостъю
- Крс - коэффициент проницаемости в направлении основы и утка.
По геометрической форме образца
- оболочечная или трехмерная САО-модель с разбивкой на конечные прямоугольные и/или треугольные элементы.
Требования к геометрической форме образцов при моделировании
2,5Б - для моделирования ПКМ на поверхностных моделях, которые представляют собой срединные поверхности с эквивалентными свойствами. Недостатком данного подхода является невозможность определения объемных эффектов - например, пористости по толщине композита.
Для моделирования полной объемной картины ПКМ рекомендуется использовать ЗБ-расчет [17]. Так, для моделирования процесса инфузии при изготовлении композитной оболочки арочного элемента определены значения характеристик связующего ВСЭ-38 (изменение вязкости по времени) и наполнителя БШех 4716 (коэффициент проницаемости преформы с учетом распределительной сетки).
По результатам исследования реологических свойств установлено, что вязкость связующего ВСЭ-38 при температуре 25°С изменяется с 0,37 до 1 Па с за 214 мин, что вводит ограничения по времени пропитки - не более 3,5 ч.
По анализу скорости течения жидкости через армирующий материал при подаче с одного торца образца, а отводе - с противоположной стороны, проводили исследование фильтрационных характеристик и определяли коэффициент проницаемости преформы на основе наполнителя БШех 4716 с использованием распределительной сетки: 2,Г10-9 м2 (рис. 2).
Силиконовый герметик
Воздушный зазор
Подача связующего
Ткань
Рис. 2. Схема заготовки для эксперименталь- отвод ного определения коэффициента проницаемости связующего
По результатам проведенных исследований скорости фильтрации через армирующий наполнитель установлено, что при градиенте давления в 0,3 ат (0,03 МПа) поло-
жение фронта жидкости относительно точки подачи сместилось на 210 мм за 60 с, что позволяет определить коэффициент проницаемости:
где х - координата распределения фронта связующего; ц - вязкость; Р - давление; т - продолжительность пропитки; ф - объемное содержание армирующего наполнителя.
Полученное высокое значение коэффициента проницаемости обусловлено наличием на поверхности распределительной сетки, применение которой увеличивает скорость фильтрации связующего в направлениях основы и утка для обеспечения полной пропитки изделия. Значения вязкости и проницаемости являются ключевыми характеристиками для использования в качестве входных данных для моделирования пропитки в программе РЛМ-ЯТМ [18-23].
Для проведения расчета геометрическая форма композитной оболочки арочного элемента с величиной пролета 8000 мм, диаметром в сечении 300 мм и высотой подъема 1000 мм построена в виде оболочечной САЛ-модели в программе ЗоННЕёре (рис. 3).
С целью минимизации влияния граничных условий на точность последующего расчета осуществлен переход от криволинейной геометрической к эквивалентной прямолинейной САЛ-модели.
После перевода САЛ-модели из формата в программу FisualMеsh сгенерирована сетка с разбивкой на конечные треугольные элементы. Далее в программе РАМ-RTM заданы граничные условия и положение точек инжекции (подачи связующего) и вакуумирования для первой схемы пропитки. Особенностью данной схемы является наличие точки подачи для инжекции связующего по центру конструкции и двух - для отвода связующего, расположенных по торцам моделируемого объекта.
По результатам моделирования определено положение фронта связующего для прогнозирования времени, затрачиваемого на пропитку при изготовлении композитной оболочки по первой схеме. Установлено, что через 4 ч после начала инжекции полной пропитки не произошло, что не соответствует технологической жизнеспособности связующего ВСЭ-38 (рис. 4).
2
к = Х -Ф-Ц
ехр 2Р-Х '
Рис. 3. САБ-модель композитной оболочки арочного элемента
Рис. 4. Продолжительность пропитки композитной оболочки по первой схеме
С учетом полученных результатов моделирования по первой схеме расположения точек подачи и отвода связующего в программе РАМ-КТМ скорректированы граничные условия и положение точек инжекции и вакуумирования для второй схемы пропитки. Особенностью данной схемы является наличие трех точек подачи связующего для инфузии и двух - для отвода, расположенных между точками подачи на одинаковом удалении.
По результатам имитационного компьютерного моделирования определено положение фронта связующего для прогнозирования времени, затрачиваемого на пропитку при изготовлении по второй схеме (рис. 5).
Рис. 5. Положение фронта связующего через 60 мин после начала пропитки
Продолжительность полной пропитки по второй схеме составила 2 ч, что соответствует технологической жизнеспособности связующего.
Обсуждение и заключения
По результатам сравнения данных моделирования процесса инфузии при изготовлении композитной оболочки с различными вариантами расположения точек подачи связующего установлено сокращение продолжительности полной пропитки в 3 раза путем оптимизации схемы при добавлении дополнительных двух точек подачи связующего и одной точки вакуумирования.
По результатам моделирования процесса инфузии композитной оболочки арочного элемента оптимизировано расположение точек инжекции для пропитки с учетом жизнеспособности связующего, а также установлено, что применение средств компьютерной имитации позволяет определить продолжительность пропитки, а также спрогнозировать характер распределения фронта связующего еще на стадии расчета.
ЛИТЕРАТУРА
1. Каблов E.H. Инновационные разработки ФГУП «ВИАМ» ГНЦ РФ по реализации «Стратегических направлений развития материалов и технологий их переработки на период до 2030 года» //Авиационные материалы и технологии. 2015. №1 (34). С. 3-33.
2. Каблов E.H. Материалы и химические технологии для авиационной техники //Вестник Российской академии наук. 2012. Т. 82. №6. С. 520-53G.
3. Каблов E.H., Каримова С.А., Семенова Л.В. Коррозионная активность углепластиков и защита металлических силовых конструкций в контакте с углепластиком //Коррозия: материалы, защита. 2011. №12. С. 1-7.
4. Каблов E.H., Сиваков Д.В., Гуляев И.Н., Сорокин КВ., Федотов М.Ю., Гончаров В.А. Методы исследования конструкционных композиционных материалов с интегрированной электромеханической системой //Авиационные материалы и технологии. 2010. №4. С. 17-2G.
5. Раскутин А.Е., Гончаров В.А. Компьютерное моделирование технологического процесса изготовления ПКМ методом вакуумной инфузии //Авиационные материалы и технологии. 2012. №S. С. 286-291.
6. Душин М.И., Мухаметов P.P., Платонов A.A., Меркулова Ю.И. Исследование фильтрационных характеристик армирующих наполнителей и связующих при разработке технологии безавтоклавного формования полимерных композиционных материалов //Авиационные материалы и технологии. 2013. №2. С. 22-25.
7. Григорьев М.М., Коган Д.И., Твердая О.Н., Панина H.H. Особенности изготовления ПКМ методом RFI //Труды ВИАМ. 2013. №4. Ст. 03 (viam-works.ru).
8. Соколов ИИ, Раскутин А.Е. Углепластики и стеклопластики нового поколения //Труды ВИАМ. 2013. №4 (viam-works.ru).
9. Душин М.И., Хрульков A.B., Мухаметов P.P., Чурсова Л.В. Особенности изготовления изделий из ПКМ методом пропитки под давлением //Авиационные материалы и технологии. 2012. №1. С. 18-26.
10. Михайлин Ю.А. Конструкционные полимерные композиционные материалы. СПб.: Научные основы и технологии. 2008. 822 с.
11. Чурсова Л.В., Душин М.И., Коган Д.И., Панина H.H., Ким М.А., Гуревич Я.М., Платонов A.A. Пленочные связующие для RFI-технологии //Российский химический журнал. 2010. T. LIV. №1. С. 63-66.
12. Гончаров В.А., Шиенок A.M., Сорокин К.В., Федотов М.Ю., Зеленский П.В. Моделирование инфузионных технологий изготовления слоистых полимерных композиционных материалов //Все материалы. Энциклопедический справочник. 2013. №1. С. 43-49.
13. Хрульков A.B., Душин М.И., Попов Ю.О., Коган Д.И. Исследования и разработка автоклавных и безавтоклавных технологий формования ПКМ //Авиационные материалы и технологии. 2012. №S. С. 292-3G1.
14. Гончаров В.А. Моделирование инфузионных технологий изготовления ПКМ на основе эпоксидных связующих и углеродных тканей Porcher //Новости материаловедения. Наука и техника. 2013. №5. Ст. 08 (materialsnews.ru).
15. Бунаков В.А., Головкина Г.С. Армированные пластики. М.: МАИ. 1997. 4G4 с.
16. Antonucci V., Esposito M., Ricciardi M.R., Raffone M., Zarrelli M., Giordano M. Permeability characterization of stitched carbon fiber preforms by fiber optic sensors //Express Polymer Letters. 2G11. V. 5. P. 1G75-1G84.
17. Porto J., Letzow M., Dos Santos E., Amico S., Souza J., Isoldi L. Computational modeling of RTM and LRTM processes applied to complex geometries //EngenhariaTérmica (Thermal Engineering). 2G12. V. 11. P. 93-99.
18. Котов B.B. Моделирование производства композиционных материалов - вызов времени //Перспективные материалы для авиакосмической промышленности. 2011. №1. С. 21-23.
19. Li J., Zhang C., Liang R., Wang B., Walsh S. Modeling and analysis of thickness gradient and variations in vacuum-assisted resin transfer molding process //Polymer Composites. 2GG8. V. 29. P. 473-482.
20. Гончаров B.A., Федотов М.Ю., Сорокин K.B., Раскутин А.Е. Моделирование стадии пропитки получения углепластиков на основе эпоксидных связующих для RTM и VARTM технологий //Справочник. Инженерный журнал. 2013. №12. С. 24-28.
21. Гончаров B.À., Федотов М.Ю., Шиенок A.M. Моделирование полимерных композиционных материалов /В сб. трудов конф. «Новые материалы и технологии глубокой переработки сырья - основа инновационного развития экономики России». М.: ВИАМ. 2010. С. 8.
22. Гончаров В.А., Федотов М.Ю. Моделирование физико-механических свойств интеллектуального углепластика с оптоволоконными сенсорами /В сб. материалов. V Международной конф. «Деформация и разрушение материалов и наноматериалов» DFMN 2013. М. 2013. С. 822-823.
23. Гончаров В.А., Сорокин КВ., Федотов М.Ю., Раскутин А.Е. Моделирование технологического процесса изготовления методом инфузии ПКМ на основе связующего ВСЭ-21 и ткани Tenax-e Style 450 //Труды ВИАМ. 2014. №10. Ст. 04 (viam-works.ru).
