Анализ некоторых методов соединения ячеек в механическом тетрахиральном метаматериале | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/3

Анализ некоторых методов соединения ячеек в механическом тетрахиральном метаматериале

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

Analysis of some methods of integration of cells in a mechanical metamaterial.pdf Метаматериалы - это искусственно созданные материалы, уникальные свойства которых обусловлены их микроструктурой, а не химическим составом базового материала. Другими словами, механическое поведение метаматериалов ближе к поведению конструкций, а не материалов. В материалах элементарная ячейка представляет собой набор атомов, в то время как метаматериал строится из крупных ячеек, на порядки больше атомного значения. Одним из методов создания метаматериалов является механическое изготовление иерархической ячеистой структуры на мезо- или макромасштабе. Данный подход в основном применяется для создания механических метаматериалов. Механические метаматериалы характеризуются необычными механическими свойствами и классифицируются по константам упругости: (1) модуль Юнга, (2) сдвиговой модуль, (3) объемный модуль, (4) коэффициент Пуассона [1]. 28 Ахметшин Л.Р., Смолин И.Ю. Анализ некоторых методов соединения ячеек Среди прочих видов метаматериалов ячеистые структуры достигают наивысшей эффективности благодаря меньшему удельному весу [2]. В последние годы интерес к исследованиям в области ячеистых структур метаматериалов расширился от чисто прикладных до фундаментальных исследований свойств. Очень популярным типом структур метаматериалов является структура с хиральными элементами. Хиральность - это свойство объекта не накладываться на свое зеркальное отображение. Хиральность может быть как левосторонней, так и правосторонней [3]. Простой хиральный элемент имеет центральное кольцо и ребра, выходящие из него [4]. Количество ребер будет определять название хиральной структуры. Впервые трехмерный образец из метаматериала с хиральной структурой, состоящей из кольца и четырех связок, был продемонстрирован в работе [5]. Авторы этой работы показали, что при одноосном нагружении стержня из механического метаматериала он проявляет необычный эффект, заключающийся в скручивании стержня. Полученный результат является аналогом оптической активности и именуется «механическая активность». Метаматериалы активно развиваются в настоящее время и уже продемонстрировали исключительный потенциал для широкого спектра применений, в частности для теплоизоляции, электродов аккумуляторов, носителей катализаторов и гашения акустической, вибрационной или ударной энергии [6]. Вызывает интерес их применение для преобразования или демпфирования механических волн. В метаматериалах энергия может быть сконцентрирована в элементах структуры и эффективно поглощена [7]. В биомедицинской инженерии существует множество перспектив применения метаматериалов: в качестве протезов, имплантатов, стентов, скаффолдов, расширителей, швов, фиксаторов связок / мускулов, бинтов [8]. Особую актуальность разработке изделий из механических метаматериалов придают успехи в развитии современных технологий 3D-печати. Технологии производства трехмерных изделий являются перспективными и конкурентоспособными по сравнению с традиционными благодаря высокой производительности и возможности создавать детали со сложной геометрией для достижения ранее недоступных свойств [9]. Настоящая работа является следствием предыдущих работ [10, 11], по результатам которых было обнаружено, что отсутствуют исследования, связанные с изучением соединения элементарных ячеек в метаматериалах. Цель данной работы - анализ различных методов соединения элементарных ячеек в механическом тетрахиральном метаматериале. Структура метаматериала Для исследования способов соединения элементарных ячеек в механическом ячеистом метаматериале можно провести анализ поведения системы двух элементарных ячеек. Чтобы создать геометрическую модель соединения двух ячеек были выделены этапы построения геометрической модели одной элементарной ячейки механического тетрахирального метаматериала (рис. 1). Затем элементарные ячейки соединяются между собой. Тетрахиральность означает, что в структуре присутствуют кольцо и четыре связки (ребра), соединенные с кольцом по касательной и взаимодействующие с другими ячейками. 29 Механика / Mechanics a Рис. 1. Структура метаматериала: эскиз (a), хиральная структура (b), элементарная ячейка (с) Fig. 1. Structures of a metamaterial: (a) scheme, (b) chiral structure, and (с) unit cell После создания элементарной ячейки ее необходимо размножить для создания макротела. Как было выяснено и будет показано в данной работе, расположение элементарных ячеек играет важную роль в деформационном поведении образца из метаматериала. Создание геометрической модели производили в модуле Design Modeler программного комплекса Ansys WB 2020R2. Для создания математической модели были приняты следующие значения параметров: l = 50 мм - длина элементарной ячейки, t = 5 мм - ширина ребра, h = 5 мм - толщина ребра, r2 =17.5 мм - внешний радиус кольцевого элемента, r1 = 12.5 мм - внутренний радиус кольцевого элемента, Ѳ - угол наклона ребра. Угол Ѳ строится между ребром метаматериала и горизонтальной плоскостью, является зависимым от параметров кольца и вычисляется следующим образом [12]: ( 2 2L V / Ѳ = arccos Методы соединения элементарных ячеек Основу метаматериала составляют элементарные ячейки. При создании метаматериала следует различать способы их соединения (рис. 2). Наиболее очевидным является метод «присоединения», когда одна ячейка присоединяется к другой, как показано на рис. 2, а. Другой метод соединения элементарных ячеек представлен на рис. 2, b. В этом случае происходит наложение краев ребер. Будем называть такой способ методом «внахлест». Видно, что при этом грань одной ячейки является также гранью другой ячейки. Различия в методах соединения элементарных ячеек заметны даже на примере двух ячеек. Видно, что при использовании метода «присоединения» в области соединения двух ячеек возникает удвоение толщины. В методе «внахлест» не возникает утолщений, и вероятно, что такой способ будет более экономичным за счет отсутствия некоторых граней ячеек. 30 Ахметшин Л.Р., Смолин И.Ю. Анализ некоторых методов соединения ячеек a b Рис. 2. Методы соединения ячеек в метаматериале: присоединение (a), внахлест (b) Fig. 2. Methods of integration of cells in a metamaterial: (a) joining and (b) overlapping Математическая постановка и особенности численного моделирования Постановка краевой задачи выражается системой уравнений теории упругости для интересующих нас полей перемещений Ыі и напряжений Oj в трехмерной постановке. Система включает уравнения равновесия (1), соотношения Коши (2) для определения деформаций через перемещения и определяющие соотношения (закон Гука) (3): V, -а,, = 0; (1) 1 ( S,2 ди. ди, + dXj dxj a- =^8-+ 2 - (2) (3) где Oj - компоненты тензора напряжений, e j - компоненты тензора деформаций, xi - пространственные координаты, ыі - компоненты вектора перемещений, X, ц -коэффициенты Ламе, 5j - символ Кронекера, i, j = 1, 2, 3. Начальные условия в момент времени t0 = 0 соответствовали отсутствию начальных напряжений и деформаций и выглядят следующим образом: а- {to ) = s- {t0 ) = 0 . (4) В работе рассмотрена задача одноосного нагружения образца из механического метаматериала. Для этого на нижней грани образца из метаматериала были наложены условия жесткой заделки, а на верхней грани задавалось перемещение, соответствующее продольному сжатию на 3%. Перемещение может быть задано со знаком плюс или минус, что соответствует растяжению или сжатию. На остальных гранях образца заданы свободные граничные условия. В этом случае граничные условия можно записать следующим образом: ыі = 0 для xi е S\\, ы2 = -0.03/ для xi е S2, а-n- = 0 для xi е S3, (5) 31 Механика / Mechanics где S1 - нижняя грань, S2 - верхняя грань, S3 - оставшиеся грани образца, l -начальная длина образца вдоль оси 2(y), nj - компоненты вектора нормали к поверхности S3. В данной работе для упругих постоянных были приняты следующие значения: E = 2.6 ГПа - модуль Юнга, ѵ = 0.4 - коэффициент Пуассона. Связь модуля Юнга и коэффициента Пуассона с постоянными Ламе задается формулами: , _ ѵЕ _E (1 + ѵ)(і - 2ѵ) ’ Ц 2(і + ѵ) . Приятые значения констант соответствуют ABS-пластику. Свойства метаматериалов в большей степени зависят не от значений упругих модулей, а от геометрии макроструктуры метаматериала. Численное моделирование проводилось методом конечных элементов в программном комплексе Ansys WB 2020R2. Элементарная ячейка рассматривалась как система стержней (балок), которые при расчете методом конечных элементов моделировались как совокупность трехмерных твердотельных элементов. Деформирование образца происходит без контактных взаимодействий составляющих его элементов. Результаты и обсуждение Сеточная сходимость Для получения достоверных результатов математического моделирования, а также для уменьшения вычислительных затрат необходимо было провести исследование на сеточную сходимость. Это позволяет определить оптимальный размер элементов конечно-элементной сетки. Поскольку наибольший интерес в исследовании представляет описание закручивания образца из метаматериала, было решено оценивать сеточную сходимость по величине угла поворота крайней грани образца из метаматериала. Рис. 3. Сеточная сходимость по осредненной характеристике - повороту грани Fig. 3. Grid convergence with respect to an averaged face rotation 32 Ахметшин Л.Р., Смолин И.Ю. Анализ некоторых методов соединения ячеек Проведенные расчеты (рис. 3) показали, что оптимальными конечными элементами будут являться тетраэдральные элементы с характерным размером 0.6 мм. Диапазон проверяемых значений размеров элементов составил 5, 1, 0.8, 0.6, 0.4, 0.2, 0.1 мм. При выборе любого из размеров конечных элементов погрешность результатов составляла не более 3.6%. Система двух элементарных ячеек Тетрахиральность позволяет двумерным структурам скручиваться вдоль некоторых направлений. Это говорит о том, что при приложении одноосной нагрузки (сжатие / растяжение) появляется дополнительная степень свободы, а именно вращение кольца [13]. Соединение ячеек в метаматериале методом «присоединения» может привести к разнонаправленному вращению граней в области контакта ячеек, что является отрицательным эффектом для закручивания всего образца из такого метаматериала (рис. 4, а). Соединение ячеек методом «внахлест» не дает подобной реакции, поэтому закономерным будет предположить, что создание образца из метаматериала представленным методом позволит достичь лучших значений механического эффекта закручивания (рис. 4, с). a b с Рис. 4. Векторы смещений в образцах из двух элементарных ячеек Fig. 4. Displacement vectors in the samples build up from two cells d Наиболее интересный эффект проявляется, если рассматривать систему из двух ячеек целиком. Как видно из рис. 4, a, b, в методе «присоединения» центр вращения возникает в месте соединения граней. Это связано с тем, что в этой точке находится центр масс системы. Вращение направлено против часовой стрелки. При рассмотрении системы ячеек, соединенных методом «внахлест», обнаружено, что такая система вращается относительно центра одного из колец тетрахиральной структуры (отмечен крестиком на рис. 4, d). То есть наблюдается асимметричная картина механического поведения, которая не замечена в системе двух ячеек, соединенных методом «присоединения». 33 Механика / Mechanics Система девяти элементарных ячеек При рассмотрении результатов для системы двух ячеек получены интересные результаты, связанные с ассиметричным деформированием. Предполагается, что подобный эффект вызван несимметричным расположением ячеек по осям в плоскости XZ. Необходимо рассмотреть образец, в котором количество ячеек, расположенных как по оси X, так и по оси Z одинаково, например, 3 х 3 ячеек. Результаты для рассматриваемой системы представлены на рис. 5. a b Рис. 5. Механическое поведение в системе из девяти ячеек, соединенных методами присоединения (а), внахлест (b) Fig. 5. Mechanical behavior in a system of nine cells integrated by the methods of (a) joining and (b) overlapping of cells По соответствующим векторам движения узлов конечно-элементной сетки можно сделать вывод, что механическое поведение при обоих способах соединения ячеек симметричное, а по распределению векторов можно заключить, что эффект скручивания в методе «внахлест» выражен более сильно, чем в методе «присоединения». Отклонение от оси X при скручивании структуры составляет 2.38 мм и 2.26 мм для методов «внахлест» и «присоединение» соответственно. Заключение Рассмотрена квазистатическая задача нагружения системы элементарных ячеек механического тетрахирального метаматериала в рамках линейной теории упругости для случая трехмерных однородных хиральных ячеистых структур. Эта работа является продолжением исследований авторов по изучению необычного эффекта «нагружение-скручивание» [10]. Математическое моделирование методом конечных элементов является удобным средством для прогнозирования деформационного поведения образца из механического метаматериала. В общем смысле проектирование материалов и метаматериалов путем математического моделирования позволяет избежать проб и ошибок при проведении натурных экспериментов. 34 Ахметшин Л.Р., Смолин И.Ю. Анализ некоторых методов соединения ячеек В работе рассмотрены некоторые методы соединения элементарных ячеек в механическом тетрахиральном метаматериале. Описаны преимущества и недостатки каждого из методов соединения в системе двух ячеек. Обнаружено, что в одном из случаев вращение происходит вокруг центра масс системы, что благоприятно влияет на стабильную работу системы двух ячеек, но уменьшает значение угла поворота. При этом возникает симметричное деформационное поведение в обеих ячейках. В другом случае показано, как отсутствие одной из граней влияет на деформацию системы двух ячеек и приводит к асимметричной картине механического поведения. Нарушение симметрии вращения может вызвать эффекты изгиба, приводящие к нестабильности. При рассмотрении системы из ячеек, симметрично расположенных по осям плоскости, состоящей из 3 х 3 ячеек, отсутствует асимметрия механического поведения. Можно утверждать, что поворот осуществляется симметрично относительно осей X и Z. Данное утверждение дает нам понимание того, что центр вращения находится в центре образцов для обоих методов. Хотя оптимизация химического состава играет важную роль в разработке новых материалов, механические и функциональные свойства большинства материалов в значительной степени определяются их структурой на уровне выше атомного. Целенаправленное проектирование «архитектурных» микроструктур позволяет получать материалы с оптимальными механическими и новыми функциональными свойствами [14]. Мы предполагаем, что продолжение исследований в выбранном нами направлении позволит создать метаматериалы и конструкции с программируемым механическим поведением.

Ключевые слова

механический метаматериал, хиральность, скручивание, соединение ячеек, численное моделирование

Авторы

ФИООрганизацияДополнительноE-mail
Ахметшин Линар РишатовичТомский государственный университет; Институт физики прочности и материаловедения СО РАНаспирант кафедры механики деформируемого твердого тела физико-технического факультета; младший научный сотрудникakhmetshin.lr@gmail.com
Смолин Игорь ЮрьевичТомский государственный университет; Институт физики прочности и материаловедения СО РАНдоктор физико-математических наук, профессор кафедры прочности и проектирования физико-технического факультета; главный научный сотрудникsmolin@ispms.ru
Всего: 2

Ссылки

Yu X., Zhou J., Liang H., Jiang Zh., Wu L. Mechanical metamaterials associated with stiffness, rigidity and compressibility: a brief review // Progress in Materials Science. 2018. V. 94. P. 114-173. doi: 10.1016/j.pmatsci.2017.12.003
Cummer S.A., Christensen J., Alu A. Controlling sound with acoustic metamaterials // Nature Reviews Materials. 2016. V. 1. Art. 16001. doi: 10.1038/natrevmats.2016.1
Grima J.N., Gatt R., Farrugia P.-S. On the properties of auxetic meta-tetrachiral structures // Physica Status Solidi (B). 2008. V. 245. P. 511-520. doi: 10.1002/pssb.200777704
Prall D., Lakes R.S. Properties of a chiral honeycomb with a Poisson’s ratio of -1 // Intemational Journal of Mechanical Sciences. 1997. V. 39, No. 3. P. 305-314.
Frenzel T., Kadic M., Wegener M. Three-dimensional mechanical metamaterials with a twist // Science. 2017. V. 358. Art. 6366. P. 1072-1074. doi: 10.1126/science.aao4640
Schaedler T., Jacobsen A., Torrents A., Sorensen A., Lian J., Greer J., Valdevit L., Carter W.B. Ultralight metallic microlattices // Science. 2011. V. 334. Art. 6058. P. 962-965. doi: 10.1126/science.1211649
Tan T., Yan Zh., Zou H., Ma K., Liu F., Zhao L., Peng Zh., Zhang W. Renewable energy harvesting and absorbing via multi-scale metamaterial systems for Internet of things // Applied Energy. 2019. V. 254. Art. 113717. doi: 10.1016/j.apenergy.2019.113717
Bhullar S.K., Lala N.L., Ramkrishna S. Smart biomaterials - a review // Review on Advanced, Materials Science. 2015. V. 40, No. 3. P. 303-314.
Kweun J.M., Lee H.J., Oh J.H., Seung H.M., Kim Y.Y. Transmodal Fabry-Perot Resonance: Theory and Realization with Elastic Metamaterials // Physical Review Letters. 2017. V. 118. P. 205901-1-205901-6. doi: 10.1103/PhysRevLett.118.205901
Akhmetshin L.R., Smolin I.Yu. Influence of unit cell parameters of tetrachiral mechanical metamaterial on its effective properties // Nanoscience and Technology: An International Journal. 2020. V. 11, No. 3. P. 265-273. doi: 10.1615/NanoSciTechnolIntJ.2020033737
Akhmetshin L.R., Smolin I.Yu. The localization of deformations in mechanical metamaterial with a twist. Numerical investigation // AIP Conference Proceedings. 2020. V. 2310. Art. 020008. doi: 10.1063/5.0034083
Akhmetshin L.R., Iokhim K. V. Analytical and Numerical Analysis of the Deformation Behavior of Tetrachiral Metamaterial // AIP Conference Proceedings. 2020. V. 2310. Art. 20006. doi: 10.1063/5.0034072
Fu M.-H., Zheng B.B., Li W.-H. A novel chiral three-dimensional material with negative Poisson’s ratio and the equivalent elastic parameters // Composite Structures. 2017. V. 176. P. 442-448. doi: 10.1016/j.compstruct.2017.05.027
Bonfanti S., Guerra R., Zaiser M., Zapperi S. Digital strategies for structured and architected materials design // APL Materials. 2021. V. 9. Art. 020904. doi: 10.1063/5.0026817
 Анализ некоторых методов соединения ячеек в механическом тетрахиральном метаматериале | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/3

Анализ некоторых методов соединения ячеек в механическом тетрахиральном метаматериале | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/3