Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

PRECONDITIONS OF REGULATION OF THERMAL CONDUCTIVITY OF POLYMER COMPOSITES BASED ON MODELING OF HEAT FLOW AND ADDITIVE TECHNOLOGIES

Pavlov S.P. 1 Makarova V.M. 1 Zlobina I.V. 1 Bekrenev N.V. 1
1 Yuri Gagarin State Technical University of Saratov
The simulation of temperature fields distribution in the volume of polymer composite materials (PCM) in the presence of macrostructure and composition inhomogeneities is performed. It is shown that by introducing into the base material structures with high thermal conductivity, which are distributed in a certain way over the cross section of the sample on the basis of the developed computer model, it is possible to significantly increase the thermal conductivity of such a composite structure taking into account the features of the product design. Examples of the organization of introduction of heat-conducting structures in the bulk of the material for different values of the required thermal conductivity and structural solutions of the unit cell of the product, taking into account the share of the input heat-conducting material. It is stated that modern additive technologies, in particular FDM, allow to practically implement the proposed approach to the regulation of thermophysical characteristics of promising PCM, used products operating under high temperature gradients: aerospace engineering, construction and other engineering structuВres designed for operation in high mountains, Arctic and Antarctic. The use of these technologies is especially important given the existing solutions for the topological optimization of the strength characteristics of polymer composite materials by forming in their volume hardening structures in accordance with the established pattern of distribution of stress fields from external loads.
polymer composite materials
thermal characteristics
thermal conductivity
temperature gradients
thermal expansion
stresses

Согласно материалам публичных аналитических докладов Сколковского Института Науки и Технологий 2014–2015 гг. [1, 2] основное внимание в ведущих промышленных странах уделяется нескольким направлениям научных исследований в интересах совершенствования производства: CAD/САМ/CAE технологии компьютерного проектирования сложных технических систем и процессов, робототехника, аддитивные технологии, композиционные материалы с иерархической структурой. При этом можно сделать вывод, что последние два направления тесно взаимосвязаны, поскольку наиболее эффективным методом формирования иерархической структуры материалов в условиях многономенклатурного производства, согласно разнообразным требованиям заказчиков, является послойное нанесение моно- или поликомпонентных составов (аддитивные процессы). При этом большое внимание уделяется качеству материалов.

За последние 20 лет объем применения композиционных, в том числе керамических, материалов в аэрокосмической технике увеличился более чем в два раза и данная динамика продолжает усиливаться [3–5]. В перспективе ожидается расширенное применение полимерных композиционных материалов (ПКМ) в ракетно-космической технике. Эта мировая тенденция характерна и для российской космической промышленности. В России разрабатывается и создается новое поколение космических аппаратов, в которых может быть заметно увеличена доля и роль композиционных и керамических материалов: перспективная экранновакуумная теплоизоляция с функцией микрометеороидной защиты для «Научно-энергетического модуля» МКС, лобовой теплозащитный экран для многоразового пилотируемого космического корабля «Федерация» и др. Планирующееся возобновление работ по исследованию и освоению Луны, в том числе с созданием на ее поверхности долговременных обитаемых модулей, также предполагает применение композиционных материалов со специальными функциональными свойствами.

При этом направления исследований в области материаловедения для космических аппаратов определяются спецификой долговременного нахождения данных объектов в открытом космическом пространстве под действием высоких градиентов температур, ионизирующих излучений, корпускулярных потоков и динамических нагрузок вследствие столкновений с микрометеоритами и «космическим мусором». В ходе исследований, выполненных космонавтами и астронавтами на долговременных космических станциях («Салют», «Скайлэб», «Мир», МКС), выявлены следующие основные факторы, влияющие на прочность и надежность конструкций из композиционных материалов: микроэрозионное повреждение поверхности и доотверждение связующего. При этом глубина повреждений и степень доотверждения зависят от состава и структуры материала, технологии его изготовления, исходной степени отверждения, максимальной температуры термоциклов и продолжительности экспонирования. Из-за температурного градиента по толщине образцов в процессе экспонирования формируется градиент физико-механических показателей ПКМ по толщине пластины. Вследствие упрочнения слоев гибридных композитов под воздействием термоциклов изменяется характер разрушения при изгибном нагружении [6]. Исследование поведения конструкций из ПКМ в космосе выявило особую важность их термостабильности в интервале температур 100…450 K в условиях периодических теплосмен, вызванных движением через теневые участки орбиты. Термостабильные конструкции создаются применением ПКМ с малыми значениями коэффициента линейного термического расширения (КЛТР) и высокой теплопроводностью. Считается [7], что для этой цели наилучшим образом подходят углепластики при условии значительного увеличения их теплопроводности в направлении больших ожидаемых температурных градиентов. Однако проблема создания теплопроводных композитных конструкций осложнена тем, что действительные физико-механические свойства ПКМ можно определить лишь после изготовления конкретной детали, т.е. в отвержденном состоянии. Кроме того, ПКМ отличаются значительной анизотропией физико-механических и теплофизических свойств, определяемой как схемами армирования, так и различием на порядок и более данных свойств для матрицы и наполнителя [8, 9]. Часто высокие физико-механические характеристики ПКМ слабо коррелируют с его теплофизическими свойствами, что создает трудности при разработке внешних конструкций как космических аппаратов, так и инженерных сооружений, предназначенных для эксплуатации в условиях Крайнего Севера, Антарктиды и высокогорья.

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

Для расчета эффективной термической проводимости элементарной периодической ячейки применяется метод гомогенизации, в котором гомогенизированные выражения могут быть получены из уравнений теплопроводности. По [10, 11] эффективная теплопроводность может быть определена в следующей форме:

pavl01.wmf (1)

с локальной координатой y.

Если обозначить pavl02.wmf [11], то, аналогично схеме, используемой для интерполяции локального тензора жесткости, характеристическая функция Ti – является решением уравнения теплопроводности:

pavl03.wmf (2)

Тогда эффективный тензор теплопроводности (1) может быть записан в виде

pavl04.wmf (3)

Для двумерной ячейки тензор эффективной термической проводимости представляется как

pavl05.wmf (4)

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

pavl06.wmf (5)

Граничные условия для определения pavl07.wmf приведены на рис. 1.

pavl1a.tif pavl1b.tif pavl1c.tif

а) б) в)

Рис. 1. Элементарная периодическая ячейка композита: а) область топологической оптимизации Y и граничные условия для задач теплопроводности (б), (в)

Граничные условия для pavl08.wmf:

pavl09.wmf

Граничные условия для pavl10.wmf:

pavl11.wmf

Материалы и методы исследования

На сегодняшний день, метод гомогенизации [10] является эффективным подходом к вычислению глобальных физических свойств композита, таких как объемный модуль жесткости или термическая проводимость.

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

Предположим, что элементарная периодическая ячейка разбита на конечные элементы и каждому элементу присваивается переменная плотность, ρn(n = 1, N). Следуя закону SIMP [12], в котором в качестве переменной управления в задаче оптимизации выступает искусственно введённая плотность ρ(x), коэффициент теплопроводности:

pavl12.wmf. (6)

Здесь k0, k1 – коэффициенты термической проводимости материалов соответственно.

Показатель степени p ≥ 1 является фактором штрафа и увеличение р приводит к более четкому решению. Показатель степени p обычно выбирают равным p = 4 или p = 5 [12].

Обычно задачи оптимизации на максимум для удобства расчетов сводятся к задачам на минимум.

Целевую функцию определим в виде

pavl13.wmf. (7)

Здесь kb – коэффициент теплопроводности, используемый для нормировки. Второе слагаемое является функцией штрафа для исключения эффекта «шахматной доски» в процессе оптимизации, h0 первоначальный размер сетки и hmax – текущий размер сетки. Величина 0 ≤ q ≤ 1 – заданный коэффициент, позволяющий сбалансировать функцию цели и функцию штрафа друг с другом. Следует отметить, что поиск минимума функции (8) соответствует задаче поиска максимума функции, стоящей после множителя (q – 1) за счет того, что q – 1 ≤ 0 [13].

Ограничения для искусственно введенной плотности ρ(x) выберем в виде

pavl14.wmf (8)

pavl15.wmf (9)

В формуле (8) A обозначает общий объем материала оптимизируемой области Y в элементарной периодической ячейке при ρ(x) = 1 и γ обозначает процент экономии материала с параметрами.

Поскольку для получения численного решения жесткость не должна полностью исчезать, положим, в неравенстве (9) δ равной некоторой достаточно малой величине, чтобы избежать сингулярности начальной матрицы жесткости при оптимизации.

Результаты исследования и их обсуждение

Рассмотрим элементарную ячейку композита из двух материалов, коэффициенты термической проводимости которых соответственно равны k0 = 10, k1 = 100. Оптимальную микроструктуру композита будем строить для различных значений от γ = 0,2 до γ = 0,9, γ – объемная доля материала с коэффициентом теплопроводности k1 = 100.

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

pavl16.wmf. (10)

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

Эффективная теплопроводность композита в зависимости от доли компонента с повышенной теплопроводностью в общем объеме

γ

0,2

0,3

0,4

0,5

 

pavT1a.tif

pavT1b.tif

pavT1c.tif

pavT1d.tif

tr(k)/2

19,9672

26,2511

33,2693

41,0807

γ

0,6

0,7

0,8

0,9

 

pavT1e.tif

pavT1k.tif

pavT1l.tif

pavT1m.tif

tr(k)/2

49,8471

59,7396

70,202

84,2823

pavl2.tif

Рис. 2. Теоретическая зависимость теплопроводности ПКМ в зависимости от % % доли компонента с условной теплопроводностью γ = 100

Построенная по результатам моделирования зависимость представлена на рис. 2.

Анализ графика (рис. 2) позволяет сделать вывод, что введение в структуру основного материала дополнительного компонента с высокой теплопроводностью, распределенного с учетом рациональной топологии, позволяет существенно повысить данную теплофизическую характеристику. Например, размещение 20 % от общего объема дополнительного компонента приводит к увеличению теплопроводности в 2 раза, а замена 80 % материала на теплопроводный компонент вызывает семикратное увеличение теплопроводности. С учетом возможных различий в механических характеристиках компонентов использование топологического структурирования ПКМ материалом с отличающейся от основного теплопроводностью позволяет путем введения относительно малой его доли, не вызывающей значимого изменения механических характеристик, существенно изменить в требуемом направлении теплофизические свойства. Таким образом, появляется инструмент направленного регулирования теплофизических свойств ПКМ с учетом их прочности и других механических
свойств.

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

В последнее время получают все большее распространение в различных областях производства аддитивные технологии прямого формирования трехмерных объектов любой сложности формы на основе разработанной компьютерной твердотельной модели [10–12]. Среди различных методов данного технологического направления технология FDM позволяет формировать композиционные изделия из различных полимерных материалов в соответствии с заданной структурой.

Ранее нами были получены положительные результаты по формированию упрочняющей топологии в виде армированного углеродными волокнами композита в объектах трехмерной печати путем выкладки армирующего материала в соответствии с полученным на основе компьютерного моделирования распределением полей напряжений. При этом полости для заполнения, имеющие сложную форму, были получены в процессе печати основы из пластика ABS [13]. Изучены теплофизические и механические свойства таких объектов, в том числе после дополнительного электрофизического воздействия [14, 15].

На основании изложенного представляется возможным предложить следующий алгоритм формирования структуры ПКМ с оптимизированными прочностными и теплофизическими характеристиками:

– на основе известных требований к изделию проектируется его конструкция;

– на основе конструкции формируется твердотельная модель изделия;

– с учетом распределения и расчетных численных эксплуатационных нагрузок моделируется распределение полей внутренних напряжений для последующего топологического упрочнения материалов с высокими модулем упругости и пределом выносливости и статической прочности;

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

– разрабатывается обобщенная твердотельная модель изделия, учитывающая упрочняющую топологию и распределение теплопроводных структур;

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

– методом аддитивных технологий на основе многокомпонентной трехмерной печати формируется изделие.

Выводы

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

Теоретически установлено, что в зависимости от доли содержания дополнительного компонента: (20–90) % условная теплопроводность композита повышается соответственно от 2 до 8 раз, что позволяет путем незначительных изменений содержания теплопроводного компонента в широких пределах направленно регулировать теплопроводность изделия из ПКМ.

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

Сочетание упрочняющей и теплофизической оптимизации структуры ПКМ позволит сформировать заданный условиями эксплуатации объекта комплекс свойств, что особенно важно для динамичных объектов, находящихся в условиях высоких температурных градиентов, например, космических аппаратов, а также инженерных конструкций и транспортной техники для зон Крайнего Севера и Антарктиды.

Исследования выполнены при поддержке гранта РФФИ № 17-03-00720 «Методология оптимизационного микроконструирования композиционных материалов для объектов сложной формы повышенной динамической прочности, послойно формируемых электротехнологическими методами».