CC BY
42
Федулов Б.Н., Федоренко А.Н., Антонов Ф.К., Ломакин Е.В. Алгоритм топологической оптимизации конструкции, выполненной из анизотропного материала с учетом параметров ориентации армирования // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 3. С. 182-189. DOI: 10.15593/perm.mech/2021.3.17
Fedulov B.N., Fedorenko A.N., Antonov F.K., Lomakin E.V. Algorithm for Topology Optimization of a Structure Made of Anisotropic Material with Consideration of the Reinforcement Orientation Parameters. PNRPU Mechanics Bulletin, 2021, no. 3, pp. 182-189. DOI: 10.15593/perm.mech/2021.3.17
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 3,2021 PNRPU MECHANICS BULLETIN
https://ered.pstu.ru/index.php/mechanics/index
Б01: 10.15593/регш.шееЬ/2021.3.17 УДК 531
АЛГОРИТМ ТОПОЛОГИЧЕСКОЙ ОПТИМИЗАЦИИ КОНСТРУКЦИИ, ВЫПОЛНЕННОЙ ИЗ АНИЗОТРОПНОГО МАТЕРИАЛА С УЧЕТОМ ПАРАМЕТРОВ ОРИЕНТАЦИИ АРМИРОВАНИЯ
19 1 о
Б.Н. Федулов1'2, А.Н. Федоренко1, Ф.К. Антонов3, Е.В. Ломакин2
1
Сколковский институт науки и технологий, Москва, Россия
2
Московский государственный университет им. М.В. Ломоносова, Москва, Россия
3
ООО «Анизопринт», Москва, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 9 ноября 2020 г. Принята: 9 июля 2021 г. Опубликована: 22 октября 2021 г.
Ключевые слова:
топологическая оптимизация, анизотропная топологическая оптимизация, выбор ориентации волокна, полимерные композиты, аддитивные технологии, термопластичное связующее, оптимизация армирования.
Современные методы проектирования включают в себя разнообразные подходы к оптимизации конструкции. Одним из хорошо себя зарекомендовавших методов является топологическая оптимизация. Существуют различные формы реализации такого вида анализа конструкций, но наиболее распространен вариант использования плотности материала как обобщенного параметра. Развитие технологий 3D-печати усиливает интерес к алгоритмам такого типа. При этом, как правило, речь идет о печати металлом или пластиком, что подразумевает изотропию материала. Появление систем печати непрерывным волокном, развитие выкладочных машин, работающих с композитными пленками, позволяет управлять не только позицией материала, но и выбирать локальную ориентацию армирующих элементов. Такие системы требуют развития алгоритмов оптимизации с учетом дополнительных параметров, характеризующих анизотропию материала. Традиционная задача топологической оптимизации в этом случае должна быть прежде всего модифицирована для модели ортотропного материала с ограничениями на величину угла поворота линий армирования вдоль траектории печати. В данной работе предложен вариант такого алгоритма топологической оптимизации, реализованный на основе комплекса конечно-элементного моделирования Abaqus, и рассмотрены примеры решения типовых задач.
©ПНИПУ
© Федулов Борис Никитович - д.ф.-м.н., в.н.с., e-mail: fedulov.b@mail.ru, : 0000-0002-6797-2424. Федоренко Алексей Николаевич - к.т.н., с.н.с., e-mail: alexey.n.fedorenko@gmail.com, : 0000-0002-3260-7531 Антонов Федор Константинович - к.ф.-м.н., дир.
Ломакин Евгений Викторович - д.ф.-м.н., зав. каф., e-mail: evlomakin@yandex.ru, : 0000-0002-8716-5363.
Boris N. Fedulov - Doctor of Physical and Mathematical Sciences, Leading Research, e-mail: fedulov.b@mail.ru, : 0000-0002-6797-2424.
Alexey N. Fedorenko - CSc in Technical Sciences, Senior Research, e-mail: alexey.n.fedorenko@gmail.com,
: 0000-0002-3260-7531. Fedor K. Antonov - CSc in Physical and Mathematical Sciences, Director Evgeny V. Lomakin - Doctor of Physical and Mathematical Sciences, Head of Department, e-mail: evlomakin@yandex.ru, : 0000-0002-8716-5363.
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
ALGORITHM FOR TOPOLOGY OPTIMIZATION OF A STRUCTURE MADE OF ANISOTROPIC MATERIAL WITH CONSIDERATION OF THE REINFORCEMENT ORIENTATION PARAMETERS
B.N. Fedulov1, 2, A.N. Fedorenko1, F.K. Antonov3, E.V. Lomakin2
Skolkovo institute of science and technology, Moscow, Russian Federation Lomonosov Moscow State University, Moscow, Russian Federation
3
Anisoprint LLC, Moscow, Russian Federation
ARTICLE INFO ABSTRACT
Advanced design methods offer not only a number of methods of modeling structures but also quite profound optimization approaches. The topology optimization is one of such advantageous methods. There are several implementations of this type of structural analysis, but the most common option is to deal with the material density as a generalized parameter. The development of 3D printing technologies increases the interest in algorithms of this type. However, it mostly refers to printing with metal or plastic, which implies the material isotropy. Advances in continuous fiber printing systems and automatic placement machines using composite tapes make it possible for engineers to control not only a material's position, but also to choose the local orientation of reinforcing elements. Such systems require optimization algorithms taking into account additional parameters characterizing the material anisotropy. In this case the traditional problem of the topology optimization should first be modified for the model of an orthotropic material with restrictions on the value of the angle of rotation of the reinforcement lines along the print path. In this paper, an idea of such a topology optimization algorithm is proposed, which is implemented using the Abaqus finite element modeling facilities, and examples of solving typical problems are considered.
©PNRPU
Received: 9 November 2020 Accepted: 9 July 2021 Published: 22 October 2021
Keywords:
topology optimization, anisotropic topology optimization, fiber orientation choice, polymer composites, additive technologies, thermoplastic matrix, reinforcement optimization.
Введение
В последние годы активно развиваются технологии 3Б-печати полимерами с добавлением армирующих компонентов [1-11]. Одним из способов является печать расплавленным термопластичным связующим совместно с непрерывным волокном на основе стекло-или углеволокна, предварительно пропитанным связующим для обеспечения необходимой жесткости. Такое связующее обеспечивает связь между монофила-ментами в жгуте волокна, а термопластичная матрица обеспечивает хорошие межслойные механические свойства и высокую ударную вязкость композита [5-6]. Объемное содержание волокна при такой технологии составляет до 30 %. Данный метод является хорошей заменой громоздким и дорогостоящим роботизированным комплексам выкладки в тех случаях, когда требуется изготовить новые детали относительно небольших размеров (детали БПЛА, кронштейны, лопасти) в малом количестве. Среди производителей оборудования для 3Б-печати непрерывным волокном стоит выделить компании «Анизопринт» (Россия) и Markforged (США).
По сравнению с выкладочным комплексом головка 3Б-принтера способна как плавно изменять направление ориентации волокна, так и достаточно резко повернуть траекторию выкладки. Это позволяет технически реализовывать оптимальные ориентации слоев, которые могут быть определены расчетными методами. Например, реализовать укладки с небольшим варьированием
направления волокна внутри одного слоя для повышения устойчивости композитной панели [12].
С развитием технологии 3D-печати геометрия изделий менее подвержена производственным ограничениям. Это в свою очередь делает актуальной задачу оптимизации геометрической формы изделия в очень разнообразных постановках. Многие исследования при решении таких задач используют подход, известный как топологическая оптимизация. В работах [13-28] предложены общие подходы и алгоритмы решения задачи топологической оптимизации с целью максимизации жесткости изделия (или минимизации энергии упругих деформаций) материала. Численная реализация алгоритмов обычно осуществляется в совокупности с решением задачи деформирования методом конечных элементов.
При рассмотрении анизотропного материала увеличивается количество параметров, входящих в постановку задачи, а также ограничений на них. Наиболее важным из таких параметров является угол, задающий направление печати волокном. Одним из способов является комбинирование топологической оптимизации с дискретным выбором ориентации волокна внутри конечного элемента оптимизируемого объема. В данном направлении можно выделить работы [29-30]. Тем не менее они имеют ряд недостатков: в первой рассматриваются только дискретные непроизвольные углы ориентации, во второй работе авторы анализируют балочные конструкции; также в обоих случаях рассматриваются
плоские задачи. В предлагаемой работе рассматривается более общий алгоритм поиска оптимального решения, основанный на градиентном методе. Предлагаемый подход наиболее эффективно можно использовать при работе с 3Б-печатью композитным материалом, при создании изделий на основе выкладочных машин, а также более абстрактных конструкционных элементов, если их механические свойства можно улучшить анизотропной линейно-упругой моделью материала.
1. Постановка задачи оптимизации
Топологическая оптимизация конструкций позволяет выбрать оптимальную геометрию изделия, при этом возможно появление полостей, новых свободных поверхностей, изменение топологии оптимизируемой детали. Для математической постановки оптимизации, следуя [8, 9], предположим, что оптимизируемое тело О разбито на большое количество малых в сравнении с О подобластей Ои (рис. 1). При работе алгоритма в сочетании с методом конечных элементов в качестве Оп удобно рассматривать пространство, занимаемое самими элементами, используемыми в моделировании. Теперь, если в такой дискретной постановке соотнести каждую подобласть - элементу Оп параметр ри, так, что если рп = 0 , то материал в данной области отсутствует, а если ри = 1, то материал присутствует, задача формализуется в некотором приближении. Степень точности, как и в методе конечных элементов, определяется количеством таких подобластей и соответствующих им параметров р . Таким образом, задача оптимизации состоит в выборе такого набора нулей и единиц { р1, р2, р3 ... ри...}, чтобы достичь необходимого
оптимума. В расслабленной постановке предполагается, что значения параметров меняются непрерывно в диапазоне 0 < ря < 1, это позволяет применять градиентные методы оптимизации, а после получения решения по степени близости параметра р к 0 или 1 интерпретировать результат для выбора дискретного значения [1315]. При этом можно положить, что р - это функция координат, и если точка среды имеет координаты x внутри элемента с индексом п, то р (х) = р (х сПл),
в противном случае ря( х) = 0 (х
Определяющие соотношения для материала, свойства которого зависят от параметра ри, записываются
в виде а = рр„Щиеу, где p - безразмерный параметр,
как правило, больший единицы. Тогда можно рассмотреть, например, задачу минимизации энергии деформаций внутри О по всем рп:
mm
р
in \Em (р)6-
где р = {р1,р2,...,ря},Еви = рр£°,; в„ - компоненты
тензора малых деформаций.
Если добавить ограничения на общую массу используемого материала
Е^ <М0'
п
и ограничения на параметры
0<р,. < 1,
то получим наиболее известную постановку задачи топологической оптимизации [1-28]. Популярность такой задачи объясняется возможностью выразить градиент для целевой функции (1), для любого р .
Рис. 1. Разбиение области оптимизации на подобласти
Fig. 1. Dividing the optimization area into sub-areas
д
дРП
I j zudQ.
Pn ¿П
(2)
При этом градиент имеет существенно локальное выражение. В случае нагружения тела сосредоточенной силой минимизация энергии деформации строго эквивалентна максимизации жесткости конструкции, что придает данной задаче понятный физический смысл. Доказательство можно найти, например, в работе [13].
В случае работы с анизотропным материалом и возможностью выбирать наравне с плотностью ориентацию материала необходимо дополнить искомые параметры. Если рассматривать выкладочные машины или системы с печатью непрерывным волокном, то на данный момент достаточно рассмотреть только один угол ориентации а, указывающий направление выкладываемого волокна.
Таким образом, с учетом ориентации материала или угла направления волокна определяющие соотношения примут следующий вид:
= Ет (р, а)6, =РР„Е°И (а)6,.
(3)
(1)
где (ап) - упругие модули анизотропного материала в зависимости от угла ориентации волокна в подоб-
ст
Q
ласти Ои. Такие формулы широко применяются и известны в теории слоистых композитов. Аналогично параметру плотности примем, что а - это функция координат х, такая что аи(х)=ая (х сОл), в противном
случае ая(х) = 0 (х . Здесь параметр а можно
рассматривать как угол выкладки волокна экструдером или композитной пленки выкладочной машины. В таком случае в качестве конкретного вида линейного
упругих модулей Е н имеет смысл рассмотреть транс-
версально изотропное тело.
Таким образом, задача оптимизации в случае изменяемого угла ориентации материала будет выглядеть следующим образом:
Ш1И
р,а
т \Ет (р,аКеиё°
(4)
где р = {р1;р2,р„}, а = {а1, а2, .... а„} ,
ХП < М0' (5)
п
0 <р. < 1.
При этом ограничения на величины углов можно не рассматривать в силу периодичности зависимостей модулей от ориентации армирования
Ет (р. а) = (р. а ± п) .
2. Градиент энергии деформации
Для решения поставленной задачи необходимо получить выражение для градиента целевой функции (4), аналогичное (2). Следуя работам [13, 14], рассмотрим принцип виртуальных работ:
{СТ (и)еу (8у)ёО = |/ЬуёО + |тЬг^,
О О -о
где и - перемещения, отвечающие равновесию; ст#. (и) -
напряжения, полученные при помощи поля перемещений и; / - массовая сила; т - поверхностные нагрузки; 8у - произвольное поле перемещений, удовлетворяющее граничным условиям; е (8у) - деформации,
полученные на поле перемещений 8у .
Теперь можно записать следующие выражения:
{Еде (р,а)еу (и)ек! (и)ёО = |/иёО + { тий8
О -о
{/ЬуёО + {тЬуЛ-{ст,.. (и)еу (8у)ёО
_ п -О п .
где е (и) - тензор деформаций, полученный на поле перемещений и.
Последняя скобка тождественно равна нулю, что делает данное выражение верным. Теперь, дифференцируя по параметру а, принимая в расчет, что 8у произвольное, но постоянное поле перемещений = 0), можно получить следующее выражение:
(р,аКе«ёО =
аа _
{/^ёп + {т*л-г~ (8у)ёо.
• да ^ да ^ да 1
Далее
дъц (и^д(р% (а)вИ (и))
да
да
д(рРЕк (а)) да
д1 (а)
е« ( и) + рРЕ'
РГ 0 ^ы
^ы (и)
'ук!
да
=рр
да
ек!(и ) + рРЕ"к! вк! }
Теперь можем записать
и
^{Е1к, (р,а)е1ек!ёО =
{/ — ёО + {т — ds да я™
да
дЕ0 (а)
-рРе! (иН(8у)ёО-{сту(и)е„ (8у)ёО.
И окончательно, подставляя вместо 8у перемещения,
полученные для равновесия
ди
да
можно заметить, что
{/^ёО +{т — ё* - Гст..(«1 ёО - 0.
О да до да О'1 ( ) 1 ^
Тогда для градиента энергии деформации можно получить следующее выражение:
да
дЕ0 (а)
/Е^ (р,а)еуенёО = рР{-Е^ен (и)е? (и)ёО. (6)
„ о да
Соотношение (6) дополняет соотношение (2) для получения градиентов по всем используемым параметрам.
Далее, используя метод градиентного спуска, можно задавшись начальными параметрами р = {р, р2,..., р„},
а = {а, а, . ., а„} , получать оптимальные решения или, более строго, точки локальных минимумов системы.
3. Численная реализация
Имея полный набор градиентов от целевой функции, возможно решать задачу оптимизации методом градиентного спуска.
Для параметра плотности р можно воспользоваться хорошо зарекомендовавшим себя вариантом поиска минимума в задачах топологической оптимизации [13]:
п
„к +1 p„
max{(1 -C)pf,Pmin} если p^ <max{(l"C)pf,Pm,„}, mm{(1 + C)pf,1} если mm{(1 + C)pf,l}<P^, pK-B^ в остальных случаях,
Вк = Л--^ JEjjUEj,s„dO
pn П„
Здесь Лк - множитель, который выбирается на каждом шаге из условия ограничения массы (6), степень параметра плотности р регулирует контрастность результата, типовое значение 3, параметры q и ^ регулируют скорость сходимости, типовые значения 0,5 и 0,2 соответственно.
Для параметра угла ориентации никаких ограничений в нашей задаче нет, поэтому используем классический градиентный метод:
к О____ „,к л ^ к о
ак+1 = аи 1
ак - а° если ак - AK < аК - а°
ак+а0 ----- ~к
ак - A„
если ак - AK >а^ + а° в остальных случаях,
Ак = X
(pf )p J ^ 8н8я dO
9а
^IdiJ
где X, у и а° - параметры, контролирующие процесс сходимости. Типовые значения параметров X = 1,
у = 0,5 и а° = 5°.
4. Примеры решения конкретных задач
Предложенные алгоритмы были всторены в систему конечно-элементного анализа Abaqus, что позволяет рассмотреть конкретные примеры анализа.
В качестве первой задачи рассмотрим изгиб консоли. Ограничения по массе выберем 30 % от начальной массы. Свойства материала, получаемые при изготовлении на системах 3Б-печати компании «Анизопринт», приведены в таблице.
Механические свойства материала
Material mechanical properties
E1, ГПа E2, ГПа G12, ГПа V12
145 8,34 4,71 0,233
Конкретные значения упругих констант не столь важны для рассматриваемого анализа, важны только относительные величины. При изменении величины нагрузки, масштаба модели и пропорциональном увеличении жесткостей результат анализа не изменится.
На рис. 2 показана схема нагружения и результат оптимизации: распределение параметров плотности р и распределение углов ориентации материала а. Угол ориентации отсчитывается от нулевого направления, которое схематично обозначено на рисунке и совпадает с ребром оптимизируемого объема. Начальными параметрами р и а были 0,3 и 0 соответственно для всех элементов анализируемого объема.
¡к
I !й
в
Рис. 2. Изгиб консоли: а - постановка задачи, плоскость выкладки волокна; б - распределение плотностей р; в - распределение параметров угла ориентации а°
Fig. 2. Bending of the cantilever: a - statement of the problem, a plane of the fiber placement; b - density distribution р, c - distribution of the orientation angle parameters a°
Видно, что значения плотности не везде равны строго единице, это связано с тем, что в анализе используется расслабленная постановка, где параметр может изменяться в пределах 0 < p < 1. На рис. 2 скрыт материал, плотность которого меньше p < 0,3. В соответствии с уравнением (3) жесткость материала падает при уменьшении параметра плотности, причем в данном примере параметр p=3, что снижает жесткость до ~3 % от изначальной.
i
I
б
Рис. 3. Изгиб консоли в трехмерной постановке: а -распределение плотностей р; б - распределение параметров угла ориентации а°
Fig. 3. Bending of the console in a three-dimensional setting a - distribution of densities р; b - distribution of parameters of the orientation angle a°
у
На рис. 3 показан результат оптимизации в полной трехмерной постановке. Видно, что общая форма повторяет двумерную задачу, а наибольшая разница проявляется в местах приложения нагрузки и закреплений.
Рис. 4. Изгиб шарнирно опертой балки: а - постановка задачи, плоскость выкладки волокна; б - распределение плотностей р; в - распределение параметров угла ориентации а°
Fig. 4. Bending of a hingedly supported beam: a - the problem statement, the plane of the fiber placement; b - density distribution р, c - distribution of orientation angle parameters a°
На рис. 4 показан результат оптимизации для на-гружения шарнирно опертой балки с центральным приложение сосредоточенной нагрузки.
В качестве ограничения по массе выбрано значение 15 % от изначальной массы. Видно, что конструкция
Библиографический список
1. Manufacturing and 3D printing of continuous carbon fiber prepreg filament / Q. Hu [et al.] // Journal of Materials Science. -2018. - Vol. 53, no. 3. - P. 1887-1898.
2. UV-assisted 3D printing of glass and carbon fiber-reinforced dual-cure polymer composites / M. Invernizzi [et al.] // Materials. - 2016. - Vol. 9, no. 7. - 583 p.
3. Build orientation determination for multi-material deposition additive manufacturing with continuous fibers / Y. Zhang [et al.] // Procedia Cirp. - 2016. - Vol. 50, no. 2016. -P. 414-419.
4. 3D printing for continuous fiber reinforced thermoplastic composites: mechanism and performance / Chuncheng Yang, Xiaoyong Tian, Tengfei Liu, Yi Cao, Dichen Li // Rapid Prototyp J. - 2017. - Vol. 23.
5. Development of a two-matrix composite material fabricated by 3D printing / A.V. Azarov [et al.] // Polymer Science, Series D. - 2017. - Vol. 10, no. 1. - P. 87-90.
6. Composite 3D printing for the small size unmanned aerial vehicle structure / A.V. Azarov [et al.] // Composites Part B: Engineering. - 2019. - Vol. 169. - P. 157-163.
имеет явно вьграженный трехмерный вид, несмотря на то, что печать производится в плоскости по слоям. Интересен горизонтальный элемент, подкрепляющий основные силовые опоры. Можно сделать вывод, что в оптимальной конструкции задействуется изгибная жесткость материала.
Имеет смысл отметить, что во всех представленных примерах результатом являются воспроизводимые технологически конструкции. Видно, что изменение всех параметров довольно гладкое, но тем не менее перенос геометрии в специализированные программы подготовки к SD-печати требуют определенной работы конструктора. Таким образом, результаты такого типа оптимизации имеет смысл рассматривать как прообраз конструкции, позволяющий выбрать оптимальную силовую схему изделия.
Заключение
В работе представлен алгоритм оптимизации конструкций, получаемых посредством технологии 3D-печати с использованием непрерывного волокна. Приведены выражения для градиентов энергии деформации в зависимости от всех используемых параметров, то есть плотностей и углов ориентации материала. Рассмотрены конкретные задачи оптимизации.
Благодарность
Работа выполнена при финансовой поддержке гранта РНФ (грант № 20-11-20230).
Acknowledgements
This work was supported by a grant from the Russian Science Foundation (grant no. 20-11-20230.
7. Additively manufactured carbon fiber-reinforced composites: State of the art and perspective / N. van de Werken [et al.] // Additive Manufacturing. - 2020. - Vol. 31. - P. 100962.
8. Novel continuous fiber bi-matrix composite 3D-printing technology / A. Adumitroaie [et al.] // Materials. - 2019. -Vol. 12, no. 18. - P. 3011.
9. 3D printing of optimized composites with variable fiber volume fraction and stiffness using continuous fiber / K. Sugiyama [et al.] // Composites Science and Technology. - 2020. -Vol. 186. - P. 107905.
10. Anisogrid thermoplastic composite lattice structure by innovative out-of-autoclave process / D. Santoro [et al.] // The International Journal of Advanced Manufacturing Technology. -2020. - Vol. 109, no. 7. - P. 1941-1952.
11. Yasa E., Ersoy K.A. Review on the additive manufacturing of fiber reinforced polymer matrix composite // Proceedings of the 29th Annual International Solid Freeform Fabrication Symposium. - Austin, TX, USA, 2018. - P. 13-15.
12. Van Campen J.M.J.F., Kassapoglou C., Gurdal Z. Generating realistic laminate fiber angle distributions for optimal vari-
able stiffness laminates // Composites Part B: Engineering. - 2012. -Vol. 43, no. 2. - P. 354-360.
13. Bendsoe M.P., Sigmund O. Topology optimization: theory, methods, and applications. - Springer Science & Business Media, 2013.
14. Bendsoe M.P., Kikuchi N. Generating optimal topologies in structural design using a homogenization method. - 1988.
15. Bends0e M.P. Optimal shape design as a material distribution problem // Structural Optimization. - 1989. - Vol. 1, no. 4. - P. 193-202.
16. Bends0e M.P., Sigmund O. Material interpolation schemes in topology optimization // Archive of Applied Mechanics. - 1999. - Vol. 69, no. 9-10. - P. 635-654.
17. Sigmund O., Petersson J. Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima // Structural Optimization. - 1998. - Vol. 16, no. 1. - P. 68-75.
18. Sigmund O. A 99 line topology optimization code written in Matlab // Structural and Multidisciplinary Optimization. -2001. - Vol. 21, no. 2. - P. 120-127.
19. Sigmund O. On the design of compliant mechanisms using topology optimization // Journal of Structural Mechanics. -1997. - Vol. 25, no. 4. - P. 493-524.
20. Sigmund O., Maute K. Topology optimization approaches // Structural and Multidisciplinary Optimization. - 2013. -Vol. 48, no. 6. - P. 1031-1055.
21. Rozvany G.I.N. A critical review of established methods of structural topology optimization // Structural and Multi-disciplinary Optimization. - 2009. - Vol. 37, no. 3. -P. 217-237.
22. Sigmund O., Torquato S. Design of materials with extreme thermal expansion using a three-phase topology optimiza-
References
1. Hu Q. et al. Manufacturing and 3D printing of continuous carbon fiber prepreg filament Journal of materials science, 2018, vol. 53, no. 3, pp. 1887-1898.
2. Invernizzi M. et al. UV-assisted 3D printing of glass and carbon fiber-reinforced dual-cure polymer composites., Materials,
2016, vol. 9, no. 7, pp. 583.
3. Zhang Y. et al. Build orientation determination for multimaterial deposition additive manufacturing with continuous fibers., Procedia Cirp, 2016, vol. 50, no. 2016, pp. 414-419.
4. Yang Chuncheng, Tian Xiaoyong, Liu Tengfei, Cao Yi, Li Dichen. 3D printing for continuous fiber reinforced thermoplastic composites: mechanism and performance. Rapid Prototyp J 2017;23.
5. Azarov A. V. et al. Development of a two-matrix composite material fabricated by 3D printing., Polymer Science, Series D,
2017, vol. 10, no. 1, pp. 87-90.
6. Azarov A. V. et al. Composite 3D printing for the small size unmanned aerial vehicle structure., Composites Part B: Engineering, 2019, vol. 169, pp. 157-163.
7. van de Werken N. et al. Additively manufactured carbon fiber-reinforced composites: State of the art and perspective., Additive Manufacturing, 2020, vol. 31, pp. 100962.
8. Adumitroaie A. et al. Novel Continuous Fiber Bi-Matrix Composite 3-D Printing Technology., Materials, 2019, vol. 12, no. 18, pp. 3011.
9. Sugiyama K. et al. 3D printing of optimized composites with variable fiber volume fraction and stiffness using continuous fiber., Composites Science and Technology, 2020, vol. 186, pp. 107905.
tion method // Journal of the Mechanics and Physics of Solids. -1997. - Vol. 45, no. 6. - P. 1037-1067.
23. Bourdin B. Filters in topology optimization // International Journal for Numerical Methods in Engineering. - 2001. -Vol. 50, no. 9. - P. 2143-2158.
24. Duysinx P., Bends0e M.P. Topology optimization of continuum structures with local stress constraints // International Journal for Numerical Methods in Engineering. - 1998. - Vol. 43, no. 8. - P. 1453-1478.
25. Bruggi M. On an alternative approach to stress constraints relaxation in topology optimization // Structural and Mul-tidisciplinary Optimization. - 2008. - Vol. 36, no. 2. - P. 125-141.
26. Fedulov B.N., Fedorenko A.N. Residual strength estimation of a laminated composite with barely visible impact damage based on topology optimization // Structural and Multidisciplinary Optimization. - 2020. - P. 1-19.
27. Fedulov B.N., Fedorenko A.N., Lomakin E.V. Evaluation of the residual strength of structures made of composite materials based on a conservative distribution of damage parameters // IOP Conference Series: Materials Science and Engineering. - IOP Publishing, 2019. - Vol. 581, no. 1. - P. 012023.
28. Fedulov B., Fedorenko A. The analysis of the worst-case distribution of the damage in composite material imposed by a low velocity impact // Procedia Structural Integrity. - 2019. -Vol. 18. - P. 399-405.
29. Stegmann J., Lund E. Discrete material optimization of general composite shell structures // International Journal for Numerical Methods in Engineering. - 2005. - Vol. 62, no. 14. -P. 2009-2027.
30. Zhou K., Li X. Topology optimization of structures under multiple load cases using a fiber-reinforced composite material model // Computational Mechanics. - 2006. - Vol. 38, no. 2. - P. 163-170.
10. Santoro D. et al. Anisogrid thermoplastic composite lattice structure by innovative out-of-autoclave process., The International Journal of Advanced Manufacturing Technology, 2020, vol. 109, no. 7, pp. 1941-1952.
11. Yasa E., Ersoy K. A. Review on the Additive Manufacturing of Fiber Reinforced Polymer Matrix Composite., Proceedings of the 29th Annual International Solid Freeform Fabrication Symposium, Austin, TX, USA, 2018, pp. 13-15.
12. Van Campen J. M. J. F., Kassapoglou C., Gurdal Z. Generating realistic laminate fiber angle distributions for optimal variable stiffness laminates., Composites Part B: Engineering, 2012, vol. 43, no. 2, pp. 354-360.
13. Bendsoe M. P., Sigmund O. Topology optimization: theory, methods, and applications, Springer Science & Business Media, 2013.
14. Bendsoe M.P., Kikuchi N. Generating optimal topologies in structural design using a homogenization method, 1988.
15. Bends0e M. P. Optimal shape design as a material distribution problem., Structural optimization, 1989, vol. 1, no. 4, pp. 193-202.
16. Bends0e M. P., Sigmund O. Material interpolation schemes in topology optimization., Archive of applied mechanics, 1999, vol. 69, no. 9-10, pp. 635-654.
17. Sigmund O., Petersson J. Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima., Structural optimization, 1998, vol. 16, no. 1, pp. 68-75.
18. Sigmund O. A 99 line topology optimization code written in Matlab., Structural and multidisciplinary optimization, 2001, vol. 21, no. 2, pp. 120-127.
19. Sigmund O. On the design of compliant mechanisms using topology optimization., Journal of Structural Mechanics, 1997, vol. 25, no. 4, pp. 493-524.
20. Sigmund O., Maute K. Topology optimization approaches., Structural and Multidisciplinary Optimization, 2013, vol. 48, no. 6, pp. 1031-1055.
21. Rozvany G. I. N. A critical review of established methods of structural topology optimization., Structural and Multidis-ciplinary Optimization, 2009, vol. 37, no. 3, pp. 217-237.
22. Sigmund O., Torquato S. Design of materials with extreme thermal expansion using a three-phase topology optimization method., Journal of the Mechanics and Physics of Solids, 1997, vol. 45, no. 6, pp. 1037-1067.
23. Bourdin B. Filters in topology optimization., International journal for numerical methods in engineering, 2001, vol. 50, no. 9, pp. 2143-2158.
24. Duysinx P., Bends0e M. P. Topology optimization of continuum structures with local stress constraints., International journal for numerical methods in engineering, 1998, vol. 43, no. 8, pp. 1453-1478.
25. Bruggi M. On an alternative approach to stress constraints relaxation in topology optimization., Structural and multi-disciplinary optimization, 2008, vol. 36, no. 2, pp. 125-141.
26. Fedulov B. N., Fedorenko A. N. Residual strength estimation of a laminated composite with barely visible impact damage based on topology optimization., Structural and Multidisciplinary Optimization, 2020, pp. 1-19.
27. Fedulov B. N., Fedorenko A. N., Lomakin E. V. Evaluation of the Residual Strength of Structures Made of Composite Materials Based on a Conservative Distribution of Damage Parameters., IOP Conference Series: Materials Science and Engineering, IOP Publishing, 2019, vol. 581, no. 1, pp. 012023.
28. Fedulov B., Fedorenko A. The analysis of the worst-case distribution of the damage in composite material imposed by a low velocity impact., Procedia Structural Integrity, 2019, vol. 18, pp. 399-405.
29. Stegmann J., Lund E. Discrete material optimization of general composite shell structures., International Journal for Numerical Methods in Engineering, 2005, vol. 62, no. 14, pp. 2009-2027.
30. Zhou K., Li X. Topology optimization of structures under multiple load cases using a fiber-reinforced composite material model., Computational Mechanics, 2006, vol. 38, no. 2, pp. 163-170.
