Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

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

Степаненко И.С. 1 Хаймович И.Н. 1, 2
1 ФГАОУ ВО «Самарский национальный исследовательский университет имени академика С.П. Королева»
2 ЧОУ ВО «Международный институт рынка»
Рассмотрен метод подготовки модели структуры композиционного материала, армированного дискретной фазой, для проведения прочностного расчёта в CAE-системе. Метод основан на алгоритме гомогенизации по среднему полю напряжений и деформаций и учитывает влияние распределения и направленности дисперсной фазы. Современные CAE-системы моделирования процесса литья могут предложить параметры процесса лишь исходя из рекомендуемых параметров переработки материала. При этом выбор оптимальных параметров остаётся за технологом, и такая оптимизация параметров должна проводиться посредством соответствующих методов, в том числе встроенных в CAE-систему (например, метод Тагути). Однако для верификации конкретной численной модели требуется проведение её калибровки для устранения ошибки при расхождении граничных условий моделирования и реальных параметров технологического процесса, т.к. модель всегда является упрощением. В частности, можно контролировать теоретический и экспериментальный профили течения и оценивать разницу между ними. Апробация метода была проведена при имитационном конечно-элементном моделировании инжекционного литья лопатки из полиэфирэфиркетона, армированного коротким углеволокном. Полученная картина распределения и ориентации армирующей фазы в рамках концепции мультимасштабного подхода является основой для расчета анизотропных свойств материала для ячейки представительного объёма материала, переноса этих свойств к сетке конечных элементов и экспорта данных в CAE-систему прочностного анализа.
композиционный материал
имитационное моделирование
мультимасштабный подход
алгоритмы гомогенизации
статорная лопатка
полиэфирэфиркетон
анизотропия свойств
1. Гречников Ф.В., Хаймович И.Н. Разработка информационных систем управления конструкторско-технологической подготовкой производства как интегрированной базы информационных и функциональных структур // Кузнечно-штамповочное производство. Обработка материалов давлением. – 2008. – № 3. – С. 34–41.
2. Хаймович И.Н. Автоматизация проектирования объектов заготовительно-штамповочного производства компрессорных лопаток авиационных двигателей // Известия вузов. Авиационная техника. – 2014 – № 2. – С. 44–48.
3. Хаймович И.Н. Методология организации согласованных механизмов управления процессом конструкторско-технологической подготовки производства на основе информационно-технологических моделей: Автореф. дис. док. техн. наук. – Самара, 2008. – 26 с.
4. Хаймович И.Н., Клентак Л.С. Усовершенствование методов сглаживания сложных поверхностей с использованием интерполяционных сплайнов //Фундаментальные исследования. – 2013. – № 10–12. – С. 2634–2638.
5. Хаймович И.Н., Хаймович А.И. Проектирование и реализация системы автоматизированного проектирования штамповки компрессорных лопаток из титановых сплавов // Известия высших учебных заведений. Цветная металлургия. – 2015. – № 2. – С.37–43.
6. Хаймович И.Н., Хаймович А.И. Процедурные правила разработки и согласования бизнес-процессов кузнечно-штамповочного производства // Вестник Самарского государственного аэрокосмического университета им. академика С.П. Королева (национального исследовательского университета). – 2008. – № 1(14). – С. 248–252.
7. Doghri I., Ouaar A. Homogenization of two-phase elasto-plastic composite materials and structures: study of tangent operators, cyclic plasticity and numerical algorithms //International Journal of Solids and structures. – 2003. – Vol. 40. – № 7. – P. 1681–1712.
8. Fu S.Y., Lauke B. Effects of fiber length and fiber orientation distributions on the tensile strength of short-fiber-reinforced polymers //Composites Science and Technology. – 1996. – Vol. 56. – № 10. – P. 1179–1190.
9. Ogierman W., Kokot G. A study on fiber orientation influence on the mechanical response of a short fiber composite structure //Acta Mechanica. – 2016. – Vol. 227. – № 1. – P. 173–183.
10. Shokri P., Bhatnagar N. Effect of packing pressure on fiber orientation in injection molding of fiber reinforced thermoplastics //Polymer Composites. – 2007. – Vol. 28. – № 2. – P. 214–223.
11. Sun X. et al. Measurement and quantitative analysis of fiber orientation distribution in long fiber reinforced part by injection molding //Polymer Testing. – 2015. – Vol. 42. – P. 168–174.
12. Swain S.S. et al. Investigation of fibre orientation using SEM micrograph and prediction of mechanical properties through micromechanical modelling // Bulletin of Materials Science. – 2016. – Т. 39. – № 3. – P. 837–846.
13. Wu L. et al. A multiscale mean-field homogenization method for fiber-reinforced composites with gradient-enhanced damage models //Computer Methods in Applied Mechanics and Engineering. – 2012. – Vol. 233. – P. 164–179.

Для литья под давлением (инжекционного литья) полимерных композиционных материалов (ПКМ), армированных дискретным высокомодульным волокном, характерно то, что распределение волокна (его направленность и длина) в готовом изделии будет значительно варьироваться. Для моделирования данного процесса можно использовать математические и компьютерные модели [1–6]. Таким образом, заданные условия процесса литья влияют на конечные механические свойства детали [8–11]. При этом композитный материал является анизотропным и гетерогенным, что сильно осложняет проведение достоверной симуляции классическим подходом, основанным на макроскопической модели. Прогнозирование свойств такого материала возможно посредством мультимасштабного подхода [7, 13].

Гетерогенное твёрдое тело состоит из фазы материала-матрицы и множества фаз-включений (например, волокон, микрополостей, частиц или микротрещин). Можно выделить 2 уровня: микроскопический и макроскопический. Гетерогенность структуры проявляется на микроуровне, в то время как на макроуровне твёрдое тело можно рассматривать как локально гомогенное.

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

1. Экспериментально де факто либо теоретически численными методами определяются ключевые параметры микроструктуры. Для коротковолокнистых ПКМ – это распределение и направленность армирующего волокна в полимерной матрице. Длина волокна известна априори.

2. На основе параметров микроструктуры (микромодель структуры материала) с применением методов гомогенизации механики сплошных сред рассчитывают анизотропные механические свойства материала в элементарном макрообъеме (макромодель структуры материала).

3. С использованием макромодели, когда в любой точке материала рассчитаны его анизотропные механические свойства, производят необходимый прочностный анализ в зависимости от граничных и начальных условий.

На практике численными методами невозможно получить механические свойства на уровне достаточно малого масштаба. Поэтому берут в рассмотрение макроуровень и считают, что каждая материальная точка является центром ячейки представительского объёма (ЯПО) с соответствующей гетерогенной микроструктурой. Классическая механика твёрдых тел применяется для макромасштаба, при этом в каждой расчётной точке значения напряжений и деформаций является граничным условием (ГУ) соответствующей ЯПО. Решением каждой ЯПО являются значения напряжений и жёсткости, которые действуют на макроуровне.

Для ЯПО с классическими граничными условиями макронапряжения и деформации равны среднеобъёмным ЯПО неизвестных полей микронапряжений и деформаций внутри ЯПО. Согласно теории линейной эластичности, соотнесение двух средних значений даёт эффективные и общие жёсткости композита на макроуровне.

Для решения проблемы ЯПО, можно использовать конечно-элементный (КЭ) метод. Несмотря на то что метод является очень распространённым и точным, создание достаточно реалистичной сетки микроструктуры является сложной задачей. Кроме того, при нелинейной постановке задачи (например, неэластичное поведение материала) возникают большие процессорные времязатраты.

Совершенно другой подход предлагает метод гомогенизации по среднему полю (ГСП), который основан на предположении о взаимосвязях средних полей напряжений и деформаций в каждой фазе ЯПО (рис. 1). По сравнению с прямым КЭ-методом и другими существующими методами перехода между макро- и микроуровнями, ГСП является самым простым в использовании и наиболее быстрым. Тем не менее, его недостатком является невозможность предоставить детальную информацию о полях напряжений и деформаций в каждой фазе, и он применим только для эллипсоидной формы включений.

step1.tif

Рис. 1. Ячейка представительского объёма композиционного материала: фаза-включение (1), фаза материала-матрицы (2)

Тензор ориентации второго порядка (aij) является эффективным способом описания ориентации волокна при инжекционном литье. Тензор состоит из девяти компонент с индексами, обозначающими направление потока, поперечное направление и направление вдоль толщины детали. Вследствие свойства симметричности aij = aji и условия нормализации (a11 + a22 + a33 = 1) девять компонент сокращаются до 5 независимых переменных:

step01.wmf

Наибольший интерес представляют следующие три компоненты:

a11 – ориентация волокна в направлении течения (значение изменяется от 0 до 1);

a22 – ориентация волокна в поперечном потоку направлении (0…1);

a13 – угол наклона в плоскости 1–3 (– 0,5…0,5).

step2.tif

Рис. 2. Схема алгоритма метода среднего поля

Одним из получивших распространение в системах численного моделирования алгоритмов гомогенизации свойств материала является метод эффективного поля [1]. Схема алгоритма представлена на рис. 2.

1. На временном интервале [tn, tn+1] исходными данными являются общие деформации на макроуровне step02.wmf и приращения деформации step03.wmf, а также значения остальных параметров к моменту времени tn.

2. Начало расчёта приращения усреднённой деформации в фазе включения.

3–5. Определение матрицы расчётных коэффициентов жёсткости материала включений step04.wmf и материала матрицы step05.wmf.

6. Расчёт усреднённой деформации в матрице.

7. Выделение изотропной части матрицы жёсткости step06.wmf.

8. Вычисление тензора Эшелби S(I, step07.wmf) для единичного включения матрицы модулей жёсткости.

9. Вычисление значений дискретных расчетных модулей жесткости step08.wmf и step09.wmf на временном шаге tn+α:

step3a.tif а) step3b.tif б)

Рис. 3. Исходная сетка инжекционного литья (а) и структурная сетка прочностного расчёта (б)

step4.tif

Рис. 4. Результат процедуры маппинга

10. Вычисление тензора концентраций деформаций.

11. Проверка сходимости расчётной и исходной (для шага tn+α) усреднённой деформаций во включении. Если |R| в пределах допуска на усреднённую деформацию в фазе наполнителя, то выходим из цикла. Иначе осуществляется новая итерация (переход на первый шаг) с уменьшенным значением step10.wmf.

12. После завершения цикла вычисляются значения макроскопических матриц жёсткости step11.wmf и тензора напряжений step12.wmf.

Рассмотренный алгоритм был апробирован при численном моделировании формирования структуры ПКМ и расчета прочности литой лопатки из ПКМ в связке CAE-систем Moldex3D-Digimat-ANSYS.

Объект моделирования – статорная лопатка компрессора ГТД. Полимерный композиционный материал – полиэфиэфиркетон, на 40 % армированный углеволокном (PEEK90HMF40).

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

Как уже было отмечено выше, для определения прочностных характеристик неоднородного материала необходимо определение направления и распределения волокон [5, 6], являющихся исходными данными вышеизложенного алгоритма гомогенизации. Эти параметры описываются тензором ориентации второго порядка. Для их графического представления генерируется собственный вектор. Собственный вектор показывает главные направления ориентации волокна, и собственные значения дают статистическое распределение волокна в направлении главных осей (от 0 до 1). Исходя из этих данных, определяется направленность эллипсоида, которая полностью описывает распределение волокна по направлениям для каждого конечного элемента.

Объёмное распределение и направленность волокон в форме тензора ориентации есть результат моделирования инжекционного литья в CAE-системе. При этом КЭ-сетка, используемая в расчёте в Moldex3D, не обязательно будет совпадать с сеткой ANSYS, поскольку алгоритмы прочностного анализа отличаются от алгоритмов моделирования течения расплава (рис. 3). Поэтому необходимо провести перенос тензоров ориентации между сетками (маппинг), что осуществляется с помощью DIGIMAT-MAP (рис. 4).

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


Библиографическая ссылка

Степаненко И.С., Хаймович И.Н. ПРИМЕНЕНИЕ МУЛЬТИМАСШТАБНОГО ПОДХОДА ДЛЯ ОПРЕДЕЛЕНИЯ СВОЙСТВ КОМПОЗИЦИОННОГО МАТЕРИАЛА ПРИ ИНЖЕКЦИОННОМ ЛИТЬЕ // Фундаментальные исследования. – 2016. – № 12-3. – С. 536-541;
URL: https://fundamental-research.ru/ru/article/view?id=41128 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674