Full Text

Введение На протяжении долгого времени люди использовали относительно сплошные материалы. Так, например, в строительстве активно применяются бетон, сталь, стекло и другие материалы. Между тем в природе для тех же целей используются также пенообразные, клеточные, ячеистые и сотовые структуры (кости, дерево, кораллы и т.п.) Их преимущество состоит в относительно высокой жесткости при малой плотности, в пористости и низкой теплопроводности [1, 2]. В настоящее время клеточные материалы могут эффективно производиться искусственно из различных видов пластика или металлов [3-14]. В связи с этим развитие методов предсказательного моделирования для определения относительной плотности, эффективных физических свойств и оптимальной геометрической структуры пенных материалов является важной задачей. Современные обзоры работ по данным направлениям можно найти в [15-17]. Для моделирования пенообразных и ячеистых структур предложено множество решений, но самым популярным, по-видимому, является ячейка Гибсона-Эшби [18], которая вместе со своими различными модификациями была подробно описана в работах авторов [2, 18-21]. В последнее время было проведено большое количество экспериментов, которые показали, что модели Гибсона-Эшби дают удовлетворительные предсказания физических свойств высокопористых пен с регулярной структурой [9, 11, 12, 22-24]. Нерегулярные пенообразные и ячеистые структуры с клетками, отличными от ячеек Гибсона-Эшби, исследовались теоретическими, экспериментальными и численными методами в [25-32]. Нерегулярные решетки из ячеек Гибсона-Эшби изучались в [33, 34]. Однако в этих работах нерегулярность понималась как различие в размерах базовой ячейки Гибсона-Эшби по осям координат, которая далее при построении представительного объема транслировалась вдоль осей уже без изменений. При расчетах нерегулярных структур в указанных работах отмечалась анизотропия эффективных свойств, вызванная нарушением свойств симметрии геометрий решеток. Отметим, что для определения эффективных свойств пенообразных и ячеистых структур помимо моделей Гибсона-Эшби можно использовать многие известные методы усреднения, применимые для композитов различных типов [15, 35, 36], а также методы вычислительной гомогенизации, учитывающие особенности пенообразных и ячеистых материалов [15]. Так, при использовании в задачах гомогенизации метода конечных элементов важны модели представительных объемов, которые специфичны для различных пенообразных, ячеистых и сотовых структур [6, 9, 15, 23, 26-28, 30-34, 37-40]. Для клеточных структур были разработаны различные модификации метода конечных элементов, например метод конечных ячеек [41, 42], а также специальное программное обеспечение: FoamGeo, GeoDict (Math2Market); Digimat FE (MSC Software Corporation, Part of Hexagon); MAVI - Modular Algorithms for Volume Images (Fraunhofer Institute for Industrial Mathematics, ITWM); Rhinoceros (Robert McNeel & Associates) и др. Тем не менее анализ проведенных исследований показывает, что модели пенообразных материалов, составленные из нерегулярных решеток с ячейками Гибсона-Эшби случайной структуры, ранее не изучались. Такие модели и составляют основной предмет исследований в настоящей статье. В работе исследуются высокопористые упругие материалы, моделируемые массивами ячеек Гибсона-Эшби как регулярной, так и нерегулярной структуры. Вначале дается математическая постановка задачи гомогенизации по методу эффективных модулей для общего случая анизотропного двухфазного материала. Обсуждаются задачи гомогенизации для материалов определенных классов физической и геометрической анизотропии, а также отмечается возможность предельного перехода от задачи гомогенизации двухфазного композита к соответствующей задаче для пористого материала. В следующем разделе подробно описываются алгоритмы создания твердотельных моделей регулярных и нерегулярных массивов ячеек Гибсона-Эшби средствами геометрического моделирования конечно-элементного пакета ANSYS. Два варианта алгоритмов позволяют создавать нерегулярные структуры с открытой пористостью как для высокопористых, так и для низкопористых пенообразных материалов. Построенные твердотельные структуры позволяют далее легко генерировать качественные гексаэдральные конечно-элементные сетки с управляемыми размерами элементов. В качестве численного примера рассмотрены регулярные и нерегулярные структуры пористого поликарбоната, рассматриваемого как изотропный материал. Проведено сравнение результатов компьютерных расчетов в ANSYS для модели Гибсона-Эшби, для регулярных и нерегулярных массивов ячеек Гибсона-Эшби в зависимости от пористости. Анализировались эффективные модули Юнга, модули жесткости и коэффициенты Пуассона, а также влияние нерегулярности геометрических структур на анизотропию свойств результирующих гомогенных материалов. 1. Теоретические аспекты метода эффективных модулей Пусть - представительный объем, состоящий из двух частей и . Объем представляет собой связанную конструкцию , состоящую из соединенных между собой трехмерных ячеек Гибсона-Эшби . Вторая область заполнена гипотетическим упругим материалом с очень малыми упругими жесткостями и моделирует пустоты в пористом материале. Примем, что материал в области является упругим, в общем случае анизотропным, и характеризуется плотностью и упругими жесткостями . В области для фиктивного материала примем , , . Таким образом, область рассматривается как упругий композит с кусочно-неоднородными плотностью и упругими модулями : , , , . При этом в случае пористого композита параметр контраста между модулями различных фаз имеет экстремально малое значение. Задача гомогенизации состоит в том, чтобы, зная материальные модули и геометрические характеристики фаз композита, определить эффективные свойства эквивалентной гомогенной среды. В качестве критерия эквивалентности будем принимать равенство энергии композита и энергии гомогенной среды в объеме при некоторых одинаковых внешних воздействиях. Рассмотрим в представительном объеме статическую задачу теории упругости для вектора перемещений со специальными граничными условиями на внешней границе . Уравнения линейной теории упругости запишем в векторно-матричном виде: (1) (2) (3) (4) где - массив, состоящий из компонент тензора напряжений ; - массив, состоящий из компонент тензора деформаций ; - матрица упругих жесткостей ; - операция транспонирования. Здесь и далее латинские индексы изменяются в пределах от 1 до 3, а греческие - в пределах от 1 до 6. При переходе от тензорных обозначений , , к векторно-матричным обозначениям Фойхта , , используются стандартные формулы соответствия между индексами , : , , , , , . Как часто принимается в теории композитов [35, 36], на внешней границе зададим линейные по пространственным координатам краевые условия для перемещений: , (5) где ; - некоторые не зависящие от постоянные величины. Положим в (5) лишь одну из компонент массива отличной от нуля: (6) где - символ Кронекера, индекс принимает какое-то целое число от 1 до 6; - некоторая постоянная величина. В соответствии с общей теорией, после решения задачи (1)-(6) и нахождения полей напряжений , можно определить столбец матрицы эффективных модулей с номером : (7) где при , при , а через обозначено осреднение по представительному объему, (8) В общем случае анизотропии требуется решить шесть задач (1)-(6) при , которые в совокупности позволяют найти эффективные модули для всех шести столбцов матрицы . Как несложно заметить, первые три задачи с являются задачами о растяжении объема вдоль осей , , соответственно, а задачи с являются задачами о сдвигах представительного объема в различных плоскостях. Описанный выше метод эффективных модулей основывается на том, что для задач (1)-(6) для композитного тела с модулями и для гомогенного тела с модулями при одинаковых внешних воздействиях (5), (6) оказываются одинаковыми средние деформации и выполняется тождество Хилла [35, 36] . Поэтому принятие условия равенства потенциальных энергий , , влечет за собой равенство средних напряжений , что и обосновывает формулы (7), (8). В случае материала ромбической системы матрица упругих жесткостей имеет девять независимых коэффициентов. Например, матрицу податливостей можно представить в виде (9) причем в силу симметрии матрицы жесткости выполняются равенства , , . Здесь - модули Юнга; - коэффициенты Пуассона; - модули сдвига. Для изотропного материала имеется только два независимых модуля, поскольку , , , . Отметим, что если исходные модули двух фаз композита обладали каким-то одним и тем же типом анизотропии, то его эффективные модули могут иметь более сложный тип анизотропии, поскольку представительный объем может не иметь тех свойств симметрии, которые имеются у материальных свойств фаз. Поэтому, решая задачи (1)-(6) при различных и определяя по (7), (8) эффективные модули, надо проверять, каким типом анизотропии будет обладать гомогенный материал. При решении задач (1)-(6) численно методом конечных элементов равенства некоторых эффективных модулей могут выполняться не точно, а с некоторой погрешностью. Естественно, для оценки типа анизотропии эффективных модулей этой вычислительной погрешностью можно пренебречь. В приводимом ниже численном примере рассматривался поликарбонат [43], который является изотропным материалом. Поэтому для регулярного представительного объема с пустотами эквивалентный гомогенный материал получается изотропным, и для нахождения его эффективных констант можно решить всего одну задачу (1)-(6), например, с . Если же представительный объем не имеет геометрической изотропии, то надо решать все шесть задач (1)-(6) и оценивать степень различия модулей из решений разных задач и получаемых модулей. При этом, как следует из доказательства для более общего случая электроупругого композита [44], матрица эффективных модулей должна быть симметричной, т.е. должны выполняться равенства . Описанная выше методология применима к двухфазным упругим композитам с условиями полного контакта на границе раздела фаз , , (10) где - скачок функции через интерфейсную границу ; - вектор единичной нормали к , внешней, например, по отношению к области . Условие (10) выражает непрерывность перемещений и вектора напряжений на межфазной границе . Для рассматриваемого случая, когда область заполнена материалом с очень малыми модулями жесткости, можно решать задачу (1)-(6) только для области с границей , , задавая условие (6) только на границе и считая границу свободной от напряжений. . (11) При малых модулях жесткости , , эффективные модули, рассчитанные по (7) из решения задачи (1)-(4) в , (5), (6) на , (10) и из решения задачи (1)-(4) в , (5), (6) на и (11), будут практически совпадать. При этом осреднение по (8) в (7) для второй задачи надо заменить на осреднение по объему функций , заданных только на : (12) Для высокопористых материалов первый вариант задачи гомогенизации более удобен для теоретического обоснования метода, а второй вариант более удобен для практического решения. В связи с этим далее будем использовать второй вариант с решением задачи гомогенизации только в . 2. Описание алгоритма построения представительного объема Согласно подходу Гибсона-Эшби [2, 18-21], на простейшем уровне трехмерный высокопористый материал или пену можно смоделировать как массив кубических элементов, состоящий из тонких балок. Примем, что отдельные ячейки имеют вдоль всех осей декартовой системы координат длину , все балки имеют одинаковое квадратное поперечное сечение со стороной , а размер внутреннего куба, измеренный как расстояние между серединами его противоположных балок, равен . На рис. 1, а изображен один из вариантов твердотельной модели ячейки Гибсона-Эшби, сгенерированной в пакете ANSYS, а на рис. 1, б показана та же модель в плоскости ( , , ). При построении ячейки Гибсона-Эшби вначале строится основная часть ячейки, состоящая из внутреннего кубического остова, а затем все ее соединения. На рис. 2, а приведен пример построения нижних балок основной части ячейки. Отметим, что структура состоит из 16 объемов двух видов. Это необходимо для того, чтобы в дальнейшем не было проблем с генерацией регулярных конечно-элементных сеток при разбиении твердотельной модели на гексаэдральные конечные элементы. Именно для этого были размещены дополнительные кубические объемы на местах контактов с соединительными балками. Затем аналогичным образом достраивается вся основная (внутренняя) часть ячейки (рис. 2, б). На последнем этапе построения геометрии отдельной ячейки строятся ее соединительные балки, и в итоге получается ячейка Гибсона-Эшби (см. рис. 1). а б Рис. 1. Вариант ячейки Гибсона-Эшби: а - общий вид в трехмерном пространстве; б - вид в плоскости Fig. 1. Variant of the Gibson-Ashby cell: a a general view in three-dimensional space; b a view in plane а б Рис. 2. Построение кубического каркаса ячейки: а - нижняя часть; б - полный кубический каркас Fig. 2. Construction of the cubic skeleton for a cell: a - a bottom part; b - a full cubic skeleton Отметим, что при варьировании размера стороны поперечного сечения балки изменяется и пористость всей структуры. Если постепенно увеличивать толщину балки, то можно прийти к случаю, когда геометрические объемы, из которых состоит ячейка, начнут пересекаться друг с другом. Поэтому для случаев больших толщин была предложена другая модель. Предположим, что . В таком случае кубики, которые расположены в углах и в серединах балок основной части ячейки, будут настолько большими, что будут пересекаться друг с другом. Чтобы избежать таких перекрытий, дополнительно выделяются линии предполагаемых пересечений больших кубов и затем по всем граничным линиям строятся непересекающиеся объемы в форме параллелепипедов. Затем аналогичным образом достраивается основная часть ячейки, и добавляются соединительные балки. На рис. 3, а приведен пример построения нижних балок основной части ячейки для случая , а один из вариантов сгенерированной итоговой ячейки показан на рис. 3, б. а б Рис. 3. Построение ячейки при больших толщинах балок: а - нижняя часть кубического каркаса; б - полная ячейка Fig. 3. Construction of a cell for large thicknesses of beams: a - a lower part of the cubic skeleton; b - a full cell Как известно [3, 9, 22], при производстве искусственных пен размеры отдельных пор могут достаточно сильно варьироваться. Поэтому далее будут рассмотрены алгоритмы для построения решеток двух видов: регулярной и нерегулярной. Регулярная решетка - это структура, которая состоит из ячеек Гибсона-Эшби одинаковой внутренней структуры. Нерегулярная решетка строится из ячеек, у которых размеры внутренних кубических каркасов выбираются случайным образом. В случае регулярной решетки ее сборка из отдельных ячеек Гибсона-Эшби не вызывает никаких трудностей. Так как все ячейки идентичны, то они идеально стыкуются между собой соединительными балками. В случае нерегулярной решетки алгоритм становится более сложным. Примем, что толщина балок одинакова для всех ячеек, это обеспечивает наиболее простой контакт соединительных балок различных ячеек. Заметим, что при больших значениях размеры основной внутренней части ячейки могут меняться в меньших пределах, чем при малых . Поэтому при толстых балках решетка становится более регулярной. В связи с этим для нерегулярной решетки мы будем рассматривать только случай, когда . Обратим внимание, что в представленной модели все ячейки имеют одинаковый размер . Но для каждой ячейки размер внутренней части выбирается случайным образом. При уменьшении размера внутренний кубик ячейки также уменьшается, но увеличиваются длины соединительных балок. Напротив, при возрастании внутренний кубик ячейки увеличивается, но уменьшаются длины балок. На рис. 4 изображены две ячейки с существенно различными размерами их внутренних частей. Рис. 4. Цепочка из двух ячеек Гибсона-Эшби с существенно различными размерами кубических каркасов Fig. 4. A chain of two Gibson-Ashby cells with significantly different sizes of cube skeletons Перед началом геометрического построения всей структуры с использованием генератора случайных чисел создается массив, который содержит значения для всех ячеек решетки. На этом этапе необходимо задать ограничения на . Очевидно, что . Но если размер будет близок к крайним значениям, то при построении появятся объемы, у которых размеры в одном измерении будут много меньше размеров в других измерениях. В такой ситуации нужно будет генерировать или очень много конечных элементов с пропорциональными размерами, или использовать элементы с непропорциональными размерами по осям. Оба варианта являются плохими для расчетов. В первом варианте размеры конечно-элементных матриц будут очень большими, и расчеты будут проводиться достаточно долго. Второй вариант может приводить к большим ошибкам в конечно-элементных решениях. В связи с этим ограничения на размер были изменены следующим образом: , После того как массив со значениями сгенерирован, в цикле строятся ячейки вдоль координатных осей по аналогии с регулярной решеткой. Так как размеры основных частей каждой ячейки различные, то их необходимо соединить дополнительными балками. Эти дополнительные балки могут быть двух видов. В случае если расстояние между балками двух ячеек небольшое (меньше ), то дополнительные балки строятся под наклоном, как показано на рис. 5, между крайней правой и центральной ячейками. В случае если расстояние между балками больше толщины балки, то дополнительные элементы строятся перпендикулярно, как показано на рис. 5, между центральной и крайней левой ячейками. а б Рис. 5. Нерегулярная цепочка из трех ячеек Гибсона-Эшби: а - общий вид в трехмерном пространстве; б - вид в плоскости Fig. 5. An irregular chain of three Gibson-Ashby cells: a - a general view in a three-dimensional space; b - a view in plane На рис. 6 для наглядности приведены примеры регулярной (а) и нерегулярной (б) решеток, включающих по три ячейки вдоль каждой оси. Как можно заметить, в нерегулярном случае ячейки с большими каркасами могут генерироваться в случайных местах решетки. Естественно, при этом ячейки с маленькими каркасами концентрируются в других местах. Такое расположение, очевидно, может приводить к геометрической анизотропии каркаса, что повлечет за собой и анизотропию эффективной среды композита. Далее это будет подтверждено численными экспериментами. Так как задачи гомогенизации решались численно в ANSYS, то твердотельный представительный объем необходимо было разбить на конечные элементы. Для этих целей был выбран стандартный конечный элемент SOLID185 из библиотеки ANSYS, который имеет восемь узлов с тремя степенями свободы перемещений в каждом узле. Средний размер конечных элементов был принят равным L/20. В результате из твердотельной модели в ANSYS была создана конечно-элементная сетка, включающая только гексаэдральные элементы. Отметим, что в силу линейности задач гомогенизации значения эффективных модулей, вычисленные по (7), (8) или по (7), (12), будут одинаковыми при геометрических структурах представительных объемов с подобными значениями , , и . Поэтому можно принять, например, (м), и , и задавать как относительные размеры (отнесенные к ). а б Рис. 6. Регулярная (а) и нерегулярная (б) решетки Fig. 6. Regular (a) and irregular (b) lattices 3. Численные результаты Как отмечалось выше, при вычислительных экспериментах рассматривался поликарбонат, который считался линейным изотропным упругим материалом. Для расчетов использовались следующие материальные константы [43]: плотность кг/м3; модуль Юнга ГПа, коэффициент Пуассона Представительный объем как для регулярной, так и для нерегулярной решетки был составлен из двадцати семи ячеек Гибсона-Эшби (по три ячейки вдоль каждой оси). Расчеты проводились при различной пористости , которая варьировалась при задании толщины балок . Для регулярных решеток были проведены численные эксперименты при пористости от 5 до 99 %, а для нерегулярных решеток - от 62 до 99 %. Все результаты далее приводятся в относительных величинах: - отношение компонент матриц жесткостей; - отношение модулей Юнга, - отношение коэффициентов Пуассона. Здесь , , - значения эффективных модулей пористого тела, а , , - соответствующие значения этих же модулей для материала матрицы композита (т.е. для тела с нулевой пористостью). Согласно одной из моделей Гибсона-Эшби для изотропного материала, относительное значение модуля Юнга равно квадрату относительной плотности, умноженному на некоторый коэффициент. Для открытых ячеек Гибсона-Эшби при преобладающих деформациях растяжения-сжатия этот коэффициент можно принять за единицу, и формулу для относительного модуля Юнга можно представить в виде [2, 18-23] (13) На рис. 7 изображены графики зависимостей от пористости относительных модулей Юнга для одной ячейки Гибсона-Эшби и для регулярной решетки из двадцати семи ячеек Гибсона-Эшби. Зеленая линия соответствует значениям относительных модулей Юнга для модели Гибсона-Эшби, рассчитанным по формуле (13). Синяя линия соответствует значениям, полученным численными экспериментами в ANSYS из решений задач гомогенизации (1)-(4) в , (5), (6) на , (11), (7), (12). Из графиков, приведенных на рис. 7, видно, что формула (13) дает значения, близкие к полученным для решетки ячеек Гибсона-Эшби, с размерами только для высокопористых пен. Пористость, % Рис. 7. Зависимости относительного модуля Юнга от пористости для модели Гибсона-Эшби и регулярной решетки Fig. 7. Dependences of the relative Young's modulus on porosity for the Gibson-Ashby model and for a regular lattice На рис. 8-11 изображены графики зависимостей от пористости относительных модулей Юнга , коэффициента Пуассона и модулей жесткости и для регулярных и нерегулярных решеток из ячеек Гибсона-Эшби. При этом на рис. 8 приведена также зеленая кривая, соответствующая значениям относительных модулей Юнга, полученным по формуле (13). Синие кривые на рис. 8-11 построены по значениям, полученным для регулярной решетки, а красные области содержат значения, полученные для нерегулярных решеток. Здесь следует отметить, что для нерегулярных решеток размеры ячеек выбирались случайным образом, и поэтому получаемые результаты изменялись от одного вычислительного эксперимента к другому. Для уменьшения влияния случайности для каждой фиксированной толщины балок проводилось по пять расчетов и осуществлялось дополнительное усреднение коэффициентов матрицы эффективных упругих модулей (9) до модулей изотропного материала. Для дальнейшей обработки результатов использовались статистические методы с вычислением доверительных интервалов по следующей методологии [45]. Пусть , , …, - значения для одной выборки относительных эффективных модулей с числом степеней свободы при малом числе испытаний . Тогда среднее арифметическое и выборочное стандартное отклонение определялись в виде , (14) а для расчета доверительного интервала среднего значения при t-распределении с доверительной вероятностью использовалась формула (15) с критическим значением t-критерия Стьюдента при , , найденным из таблицы A.3 [45]. Проводя вычисления по (14), (15) для значений , можно построить зависимости от пористости верхних и нижних значений относительных эффективных модулей с доверительными интервалами . Сглаженные кривые красного цвета получены аппроксимацией значений доверительных интервалов. Они определяют на рис. 8-11 границы областей, в которых с доверительной вероятностью 95% находятся значения эффективных модулей для нерегулярных решеток из ячеек Гибсона-Эшби. Из приведенных графиков видно, что значения, полученные для нерегулярных решеток, могут иметь достаточно большой разброс. Однако они сохраняют общую тенденцию и за исключением значений относительных коэффициентов Пуассона близки к значениям, полученным для регулярных решеток. Но, как известно [31, 46], эффективные коэффициенты Пуассона пористых материалов наиболее существенно зависят от микроструктуры, и даже тенденции их изменения от пористости могут меняться в зависимости от значений коэффициентов Пуассона соответствующих сплошных материалов и внутренней структуры представительных объемов. Пористость, % Рис. 8. Зависимости относительного модуля Юнга от пористости для модели Гибсона-Эшби и для регулярных и нерегулярных решеток Fig. 8. Dependences of the relative Young's modulus on porosity for the Gibson-Ashby model and for regular and irregular lattices Пористость, % Рис. 9. Зависимости относительного модуля жесткости от пористости для регулярных и нерегулярных решеток Fig. 9. Dependences of the relative stiffness modulus on porosity for regular and irregular lattices Пористость, % Рис. 10. Зависимости относительного модуля жесткости от пористости для регулярных и нерегулярных решеток Fig. 10. Dependences of the relative stiffness modulus on porosity for regular and irregular lattices Пористость, % Рис. 11. Зависимости относительного коэффициента Пуассона от пористости для регулярных и нерегулярных решеток Fig. 11. Dependences of the relative Poisson's ratio on porosity for regular and irregular lattices Другим интересным результатом является то, что нерегулярная решетка может не сохранять изотропные свойства материала матрицы композита. Для анализа этого явления были рассчитаны относительные отклонения эффективных модулей для композитов с нерегулярными решетками. Так, были вычислены следующие величины: % - относительные отклонения модулей Юнга; % - относительные отклонения модулей жесткости и ; % - относительные отклонения коэффициентов Пуассона. На рис. 12 приведены графики средних относительных отклонений для модулей жесткости и (красная кривая), модулей и (зеленая кривая), модулей Юнга (синяя кривая) и коэффициентов Пуассона (желтая кривая). Интересно отметить, что зависимости от пористости относительных отклонений и очень близки, а поэтому синяя кривая практические совпадает с красной и не видна на рис. 12. Пористость, % Рис. 12. Разброс относительных модулей нерегулярных решеток, характеризующий их анизотропию Fig. 12. Scatter of the relative moduli for irregular lattices characterizing their anisotropy На рис. 12 изображены аппроксимации средних относительных отклонений. Видно, что они не превышают 30 %, однако их неосредненные значения в некоторых случаях могут быть гораздо больше, особенно для высокопористых структур. Приведенные результаты подтверждают мысль, что пористые пены с нерегулярной решеткой могут иметь достаточно выраженные анизотропные свойства, даже если материалы их каркасов являются изотропными. Эта анизотропия обусловлена возможной несимметричностью геометрий и большими напряжениями в нерегулярных решетках. На рис. 13 приведены примеры расчетов двух задач гомогенизации для различных структур с нерегулярными решетками, в которых наблюдаются большие концентрации напряжений в ячейках вдоль одной оси. Обратим внимание, что каркасы этих ячеек имеют небольшие размеры относительно других ячеек. Чем меньше каркасы, тем длиннее их соединительные балки и, соответственно, наблюдаются большие напряжения в этих балках и вблизи данных ячеек [32, 33, 47]. а б Рис. 13. Концентрация напряжения на соединительных балках при деформации вдоль оси x (а) и оси y (б) Fig. 13. Stress concentration on the connecting beams during deformation along the x-axis (a) and y-axis (b) Так, на рис. 13, а изображено распределение x-компоненты напряжения после растяжения вдоль оси Ox. Из данного рисунка видно, что непосредственно в части представительного объема, где длинные соединительные балки выстроены вдоль оси Ox, наблюдается наибольшее напряжение. Аналогичная тенденция видна на рис. 13, б, но, в отличие от рис. 13, а, на нем изображено распределение y-компоненты напряжения после растяжения вдоль оси Oy другой нерегулярной решетки. Отметим, что здесь представлены результаты для регулярных и нерегулярных решеток, составленных из трех ячеек Гибсона-Эшби вдоль каждой оси. При большем числе ячеек в регулярных и нерегулярных решетках статистически получаются более близкие результаты для эффективных упругих модулей, а эффекты концентрации напряжений нивелируются. Однако для отдельных вариантов сильно нерегулярных решеток эффективные модули по-прежнему могут принимать экстремальные значения. Разработанные модели нерегулярных структур из ячеек Гибсона-Эшби могут использоваться в специализированном программном обеспечении для компьютерного дизайна пенообразных материалов, как в рамках приближений линейной теории упругости, так и для исследования их нелинейного поведения. Заключение Было проведено исследование эффективных упругих свойств пористого пенного материала, моделируемого массивами ячеек Гибсона-Эшби с регулярной и нерегулярной структурой. Разработан программный инструментарий на языке APDL ANSYS, позволяющий строить такие массивы ячеек, получать из них конечно-элементные модели с управляемыми характеристиками, решать численно задачи гомогенизации и определять полный набор эффективных модулей жесткости для материалов произвольного типа анизотропии. Вычислительные эксперименты для поликарбонатных пен показали, что для нерегулярных решеток может наблюдаться достаточно существенный разброс значений упругих модулей, особенно для эффективных коэффициентов Пуассона. Отмечено, что для нерегулярных решеток характерно появление геометрической анизотропии эффективных свойств композита, которая может достаточно сильно проявляться для высокопористых пен. С использованием разработанных конечно-элементных технологий в дальнейшем предполагается исследовать свойства представительных объемов из ячеек Гибсона-Эшби большой размерности, а также наноструктурированные пены [48] по теории поверхностных напряжений Гуртина-Мурдоха путем добавления в разработанные модели поверхностных конечных элементов [49-51].

About the authors

A. S Kornievsky

Southern Federal University

A. V Nasedkin

Southern Federal University

References

Statistics

Views

Abstract - 567

PDF (Russian) - 262