The purpose of this research is to analyze the mechanical behavior of auxetic re-entrant-based metamaterials with properties similar to Inconel 625 using homogenization techniques. Through a thorough analysis, this study investigates the displacement patterns exhibited in various materials throughout a range of thicknesses. The examination also includes analyzing how the Young’s modulus changes with varying strut thickness after homogenization. This detailed investigation provides information on the stiffness and deformation response of the material. This research advances our knowledge of the complex mechanical properties of re-entrant-based auxetic metamaterials that resemble Inconel 625 by interpreting these displacement and Young’s modulus patterns.

Architected materials, which have precisely defined architectures at the micro- and nanoscales, are a new and important class of advanced engineered materials.1–3 Unprecedented combinations of qualities, such as strength, stiffness, toughness, and multifunctionality, can be obtained by intentionally modifying the geometry of the internal material structure from the microscopic to the macroscopic level, overcoming limits inherent in traditional materials. Significant research has concentrated on lightweight lattice-based structures with extraordinary mechanical efficiency, such as microlattices and nanolattices. The micro- and nano-architectures’ high degree of control and tunability gives mechanisms for engineering materials with specific target attributes. Such complicated interior material patterns can be produced at scale using additive manufacturing methods.4–7 Because their possible property space far exceeds that of traditional materials, architected materials are poised to be a highly disruptive approach in materials science and engineering, opening up new frontiers in structural and multifunctional materials’ design.

Topology, which is concerned with the organization and geometric arrangement of these interior formations, acts as the foundational framework determining the operational properties of the material. The intricacies of topology also have an impact on qualities such as energy absorption capacities, customization potential, and the ability to fulfill many purposes concurrently. Importantly, the growing convergence of topology optimization and cutting-edge manufacturing processes, such as 3D printing, accelerates the realization of complicated interior structures. The capacity of auxetic materials—especially those in the form of lattice structures—to absorb impact energy more effectively than that of conventional materials is being studied. This feature is useful in situations where reducing the impact force is essential for user safety, such as sports equipment, helmets, and protective gear. The application of auxetic structures in tissue engineering and biomedical implants is being investigated. These materials may be appropriate for implants and scaffolds that closely resemble the mechanical properties of genuine tissues because the negative Poisson’s ratio can imitate some characteristics of biological tissues. Auxetic materials can be used to make elastic and flexible parts for wearables and electronics. These lattice structures are appropriate for applications such as stretchy conductors and sensors because they permit deformation without sacrificing the material’s integrity. It is possible to create auxetic lattice structures that function well as sieves or filters. Researchers can alter the size of the gaps in the structure by varying the lattice’s characteristics, which makes it beneficial for applications such as liquid and air filtration. Auxetic materials can be used to create lightweight structures with better energy absorption capabilities in aeronautical applications. This may aid in the creation of stronger and more effective aerospace components, such as impact-resistant materials and aircraft parts. Auxetic lattice structures can be included into mechanical joints and seals to improve deformation resistance and sealing effectiveness under different loading scenarios. In mechanical systems, this may result in increased longevity and dependability. Auxetic materials’ distinct mechanical properties can be used to reduce vibration in many applications. Researchers hope to improve the overall stability and performance of structures by lowering vibrations in mechanical systems through the use of auxetic lattice structures. Gaining knowledge of auxetic lattice structures enhances the comprehension of materials science and creative design concepts. Investigating materials with negative Poisson’s ratios creates new opportunities to modify material characteristics to suit particular uses, advancing materials science and engineering.8–11 

The importance of auxetic metamaterials stems from their extraordinary mechanical properties, which contradict conventional norms. These materials have the unusual capacity to expand under tension, bringing a variety of benefits with far-reaching ramifications. This feature correlates with exceptional toughness and resilience, making them particularly ideal for applications requiring good shock absorption, such as protective gear, sports equipment, and cushioning solutions. What distinguishes auxetic materials is their ability to provide customized mechanical reactions via intentional structural design, allowing for precise customization aligned with specific requirements. Shen et al.12 took the concept of pattern alteration triggered by buckling to an advanced level, applying it to the creation of a novel collection of three-dimensional metamaterials. These materials were engineered to exhibit a negative Poisson’s ratio across a wide range of strains. Lei et al.13 introduced an original auxetic metamaterial design. By harnessing the shape memory characteristics of its constituent elements, they achieved the ability to dynamically customize the in-plane moduli and Poisson’s ratios. Cheng et al.14 devised re-entrant unit cells featuring diverse variable stiffness factors (VSFs) with the aim of manipulating stiffness adjustability by tuning densification strain. This study revealed that the compaction points of these re-entrant structures could be precisely adjusted through the utilization of the specified VSF, thereby introducing a novel approach for the strategic design of unit cells with negative Poisson’s ratio.

The primary purpose of this research is to gain a full understanding of the mechanical behavior and properties presented by these complex materials at a macroscopic level. Because of their intricate microstructural arrangements, re-entrant-based auxetic materials exhibit strange and illogical deformation patterns. Homogenization procedures were generally employed to bridge the gap between the material’s overall macroscopic reactivity and its microscale behavior, which is produced from individual structural components. The goal of this investigation was to see how varying strut thicknesses affected the homogenized properties of these materials. These findings are expected to result in increased performance across a wide range of real-world applications and the advancement of optimal material design.

Topology optimization is a mathematical technique that optimizes the distribution of material inside a specific design area based on a number of loads, constraints, and performance requirements. The best material arrangement, or “topology,” must be determined in order to maximize a performance metric, such as stiffness or strength, while also conforming to production and weight constraints. It entails discretizing the design space into finite elements and iteratively altering each element’s density between solid and void until the best material arrangement is reached. This iterative process is powered by numerical optimization algorithms that simulate different material distributions and update densities to improve the objective function. The resulting topologies are usually characterized by holes and thin, organic features that improve structural efficiency. The mechanism of topology optimization is depicted in Fig. 1.

FIG. 1.

Mechanism of topology optimization.

FIG. 1.

Mechanism of topology optimization.

Close modal

Topology optimization generates compact, high-performance systems by adjusting the material architecture to the physics of the challenge. However, making the intricate forms it creates may be difficult. Geometric constraints and manufacturing considerations must be incorporated into the formulation in order to develop designs that are suitable for production. Because of advancements in additive manufacturing, topology optimization is becoming more prominent as an automated design process in engineering domains, such as aerospace, automotive, and biomechanics. The purpose of current topology optimization method research is to shorten solution durations, widen their scope, and make them more manufacturing-ready.

Topology has a direct impact on the bulk physical qualities and performance of these cutting-edge engineered materials, making it essential when designing architected materials. Rather than their fundamental chemical makeup, architect materials acquire their specific properties from the internal organization, connection, and geometry that have been painstakingly constructed. Architected materials can be divided into three different types according to the degree of control over geometric features and the reliance on additive manufacturing processes, as shown in Fig. 2 of this study. The two-dimensional honeycomb and three-dimensional foam arrangements stand out among these varieties due to their popularity and basic characteristics. The exceptional mechanical characteristics of two-dimensional honeycomb structures, characterized by regularly spaced hexagonal cells on a flat lattice, are particularly remarkable. These characteristics include high strength-to-weight ratios and efficient material consumption. These structures have a wide range of applications in aviation components, sandwich panels, and lightweight composites. They are similar to natural examples, such as beehives, which serve as sturdy yet mobile frames for honey storage and larval development. They are incredibly useful to engineers because of their capacity to distribute loads efficiently, absorb energy, and provide exact geometric control. However, the first class of cellular structures, which includes honeycombs and foams, has inherent restrictions that force the development of a second class, represented by lattices with nodes and struts. In honeycombs and foams, the fixed cell configuration limits the design flexibility and prevents customization to satisfy particular engineering requirements. Contrarily, lattices with nodes and struts provide for greater design flexibility, giving engineers the ability to optimize connections and geometry for improved mechanical properties. In addition, this second class may have better load-bearing capacities than honeycombs and foams, making it appropriate for applications requiring increased strength and load-bearing capacities. More design freedom is actually possible with the second class of cellular structures, which are characterized by lattices with struts and nodes. Lattice systems show stress concentration in response to external stresses, leading to concentrated areas of higher stress, particularly at strut and node connections. When there are acute angles or abrupt transitions in the lattice geometry, this phenomenon might result in a localized collapse and diminished structural integrity. The issue of stress concentration is addressed by the third class of cellular structures, known as Triply Periodic Minimal Surfaces (TPMS). Small surface areas and continuous, smooth curves of TPMS aid in more evenly dispersing stress throughout the material. The interconnected and regular structure of TPMS allows for gradual transitions between cells, as opposed to lattices with sharp edges, which lowers stress concentrations. The improved mechanical properties and material efficiency of TPMS are due to their increased surface area-to-volume ratio. TPMS are thus appropriate in circumstances needing an enhanced load-bearing capacity and a higher structural performance.

FIG. 2.

Classification of architected materials.

FIG. 2.

Classification of architected materials.

Close modal

A higher strength-to-weight ratio is frequently correlated with a larger surface area-to-volume ratio. Weight is a function of volume, and strength in materials and constructions is the capacity to bear an applied force without cracking or deforming. The strength of the structure will be enhanced in relation to its weight if you increase the surface area while maintaining the volume constant. This is because the applied forces are more evenly distributed over a greater area. Structures with less material can be created by increasing the surface area-to-volume ratio. This is crucial for a lightweight design when cutting down on material usage is a top priority. Because of their complex patterns, TPMS structures use less material to attain a large surface area, which increases their efficiency in using resources. An increased ratio of surface area to volume may be a factor in the increased stiffness. Increased stiffness and structural integrity are the results of reducing the susceptibility to deformation by distributing the forces over a wider surface area.

Auxetic metamaterials are a class of materials that, when stretched longitudinally, expand transversely due to a negative Poisson’s ratio as shown in Fig. 3. When stretched, most common materials, such as rubber or metals, have a positive Poisson’s ratio and normally compress laterally; meanwhile, auxetic materials, or those with a negative Poisson’s ratio, expand laterally. Because of their distinct mechanical characteristics, auxetic materials are frequently used in the design of structures, such as lattice systems. Auxetic metamaterials have special properties, such as increased fracture toughness, indentation resistance, shear resistance, and synclastic curvature, because of this counterintuitive behavior. Various geometries and structures, including re-entrant models, rotating polygonal models, chiral models, crumpled sheets, and perforated sheets, can be used to create auxetic metamaterials. They can also be found naturally in a few materials, such as cat skin, cancellous bone, and some cubic metals. The development of auxetic metamaterials—metals and composites—to enhance their mechanical properties is still being researched.

FIG. 3.

(a) Axial re-entrant unit-cell model, (b) axial re-entrant unit-cell porous auxetic-lattice structure, (c) transversely oriented re-entrant unit-cell model, (d) transversely oriented porous auxetic-lattice structure, and (e) two-phase auxetic-lattice structure15 [Tashkinov et al., Polymers 15(20), 4076 (2023), licensed under a Creative Commons Attribution (CC BY) license.].

FIG. 3.

(a) Axial re-entrant unit-cell model, (b) axial re-entrant unit-cell porous auxetic-lattice structure, (c) transversely oriented re-entrant unit-cell model, (d) transversely oriented porous auxetic-lattice structure, and (e) two-phase auxetic-lattice structure15 [Tashkinov et al., Polymers 15(20), 4076 (2023), licensed under a Creative Commons Attribution (CC BY) license.].

Close modal

In order to address the inherent conflicts between the strength and toughness of architected mechanical metamaterials, Wang et al.16 designed and produced a novel hybrid auxetic re-entrant metal–ceramic lattice. The synergistic effects of auxetic re-entrant metal honeycombs and infilled ceramic materials are studied experimentally and numerically, and auxetic deformation features and failure modes are characterized using the digital image correlation (DIC) technique as well as illustrated in Fig. 4.

FIG. 4.

Strain contour of the built-in ceramics obtained from DIC measurements along the horizontal direction: (a)–(c) conventional hex and (d)–(f) re-entrant hex. (a) and (d) The strain contour at the start of the compression test (0.5% in the vertical direction). (b) and (e) The strain of the compress test in the vertical direction is 2.88% for conventional hex and 3.70% for re-entrant hex. The strain contour at half time of the compression failure period is as follows: (c) and (f) The strain of the compress test in the vertical direction is 5.75% for conventional hex and 7.39% for re-entrant hex. The strain contour at the instant before compression failure16 [Wang et al., Appl. Sci. 13(13), 7564 (2023), licensed under a Creative Commons Attribution (CC BY) license.].

FIG. 4.

Strain contour of the built-in ceramics obtained from DIC measurements along the horizontal direction: (a)–(c) conventional hex and (d)–(f) re-entrant hex. (a) and (d) The strain contour at the start of the compression test (0.5% in the vertical direction). (b) and (e) The strain of the compress test in the vertical direction is 2.88% for conventional hex and 3.70% for re-entrant hex. The strain contour at half time of the compression failure period is as follows: (c) and (f) The strain of the compress test in the vertical direction is 5.75% for conventional hex and 7.39% for re-entrant hex. The strain contour at the instant before compression failure16 [Wang et al., Appl. Sci. 13(13), 7564 (2023), licensed under a Creative Commons Attribution (CC BY) license.].

Close modal

In terms of stress and micromotion distributions, the new porous femoral hip meta-implant with a graded Poisson’s ratio distribution by Ghavidelnia et al.17 was compared to three other femoral hip implants (one solid implant and two porous meta-implants, one with a positive distribution of Poisson’s ratio and the other with a negative distribution) shown in Fig. 5.

FIG. 5.

(a) The meta-implant’s CAD model and (b) its meshing finite element (FE) model, both based on the 3D re-entrant structure. (c) Complete implant and femur bone model assembly, complete with matching mesh. (d) The implant’s boundary condition and the applied load on the femur bone model 17 [Ghavidelnia et al., Materials 14(1), 114 (2021), licensed under a Creative Commons Attribution (CC BY) license].

FIG. 5.

(a) The meta-implant’s CAD model and (b) its meshing finite element (FE) model, both based on the 3D re-entrant structure. (c) Complete implant and femur bone model assembly, complete with matching mesh. (d) The implant’s boundary condition and the applied load on the femur bone model 17 [Ghavidelnia et al., Materials 14(1), 114 (2021), licensed under a Creative Commons Attribution (CC BY) license].

Close modal

Auxetic metamaterials are a distinct class of materials distinguished by the paradoxical property of becoming thinner (decreasing in the cross-sectional area) when compressed and thicker (increasing in the cross-sectional area) when stretched, in contrast to conventional materials that behave in the opposite way. Their unique internal structure leads to a negative Poisson’s ratio, which results in this peculiar behavior as shown in Fig. 6. The present work utilizes the Inconel 625 material for the development of 2 × 2 × 2 based periodic auxetic metamaterial using nTopology software. Figure 7 shows the cell map of the dimension 10 × 10 × 10 mm3. In the context of lattice systems, the term “cell map” refers to a graphic representation or a diagram that shows the configuration of individual cells within the lattice. Repeating unit cells make up lattice structures, which have patterns that fill the entire space. A cell map’s individual cells each stand in for a basic lattice repeating unit. The cell’s geometry is depicted on the map, including its size, shape, and connections to nearby cells. Understanding the mechanical behavior of the lattice and how it responds to stresses, pressures, and other stimuli requires knowledge of this information.

FIG. 6.

Unit cell structure of the re-entrant auxetic metamaterial.

FIG. 6.

Unit cell structure of the re-entrant auxetic metamaterial.

Close modal
FIG. 7.

Cell map used in the present work.

FIG. 7.

Cell map used in the present work.

Close modal

In the present work, four lattice structures having a unit cell dimension of 2 × 2 × 2 of strut thickness 0.2, 0.3, 0.4, and 0.4 mm are taken into account as shown in Fig. 8.

FIG. 8.

Representation of the 2 × 2 × 2 re-entrant auxetic metamaterial of strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 8.

Representation of the 2 × 2 × 2 re-entrant auxetic metamaterial of strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal

A combination of analytical models, numerical simulations (such as finite element analysis), and experimental testing is used to analyze the mechanical behavior of lattice systems. An analytical model includes honeycomb and truss models, periodic cell models, and homogenization techniques. A numerical simulation based on finite element analysis includes mesh based simulations, material modeling, and parametric studies. Experimental studies include compression and tensile testing, 3D printing and prototyping, and digital image correlation.

A powerful method for accurately mimicking the behavior of infinite periodic lattice structures with only a single repeating unit cell is periodic boundary conditions. This assumes that the lattice is flawlessly periodic, with no flaws or deviations. The displacements and strains must continue and repeat in all three dimensions across the unit cell’s opposing boundaries. Constraining adjacent nodes on opposing faces to have the same displacements ensures this continuity. The displacement field ux,y,z is said to be periodic if Eq. (1) is obeyed,
(1)
where Lx, Ly and Lz, are the unit cell dimensions.
Similarly, the strain field εx,y,z and stress field σ(x, y, z) are periodic if Eqs. (2) and (3) are obeyed,
(2)
(3)
This physically joins opposing borders without gaps, as though the unit cell was tiled endlessly in all directions. The effective homogeneous material response is represented by reactions and volume-averaged stresses and strains over the unit cell, which roughly correspond to the overall lattice behavior. One unit cell modeling is sufficient to reproduce the entire lattice response when periodic boundary conditions are used. When compared to examining a huge lattice, this significantly lowers the computing cost, enabling useful design and optimization research. Thus, periodic boundaries are a crucial enabler for lattice structure multiscale modeling and simulation.
Homogenization is an important technique in multiscale modeling of materials with microstructural heterogeneity. The core concept involves replacing a heterogeneous material with a fictitious homogeneous equivalent material that exhibits the same effective properties. Considering a heterogeneous material with a microstructure characterized by a periodic unit cell, the heterogeneous response is governed by equations relating stresses and strains as shown in the following equation:
(4)
where C(x) is the position dependent elasticity tensor.
The goal is to derive a homogeneous equivalent material with constant properties independent of x. This is accomplished by volume averaging the equations over a representative volume element (RVE) as shown in the following equation:
(5)
where V is the RVE volume, ⟨ ⟩ denotes volume averaging, and C* is known as the homogenized or effective elasticity tensor and relates the average stress to average strain. Boundary conditions on the RVE determine the imposed average strains based on the applied loading. C* can be extracted by solving boundary value problems on the RVE and computing average quantities. This homogenized constant tensor C* replaces the heterogeneous elasticity tensor C(x) in the larger-scale analysis.

Homogenization simulations were performed on a re-entrant auxetic lattice structure consisting of a 2 × 2 × 2 unit cell made of Inconel 625. The strut thickness was varied from 0.2 to 0.4 mm to study its effect on the overall auxetic behavior and stiffness as shown in Figs. 9(a)9(c). Periodic boundary conditions were applied to allow the unit cell response to represent the effective homogenized properties. The displacements in the principal X, Y, and Z directions along with XY, XZ, and YZ shear, respectively, were extracted for each thickness case as shown in Table I and are further illustrated in Figs. 1015. It was observed that increasing the strut thickness from 0.2 to 0.4 mm reduced displacements in all normal and shear directions. This indicates an enhancement in stiffness and auxetic effects, as the lattice expands laterally when loaded axially. The periodic boundary conditions enable the efficient prediction of the macroscale auxetic behavior by modeling a single repeating unit cell. Varying the geometry and thickness allows tuning the design for desired stiffness and Poisson’s ratio. The reduced displacements with thicker struts demonstrate that tuning strut thickness is an effective approach for tailoring the performance of re-entrant auxetic lattices.

FIG. 9.

Homogenization of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 9.

Homogenization of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 10.

Displacement in the X-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 10.

Displacement in the X-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 11.

Displacement in the Y-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 11.

Displacement in the Y-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 12.

Displacement in the Z-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 12.

Displacement in the Z-direction of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 13.

Displacement in the YZ-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 13.

Displacement in the YZ-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 14.

Displacement in the XZ-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 14.

Displacement in the XZ-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
FIG. 15.

Displacement in the XY-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 15.

Displacement in the XY-plane of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
TABLE I.

Displacement in mm obtained after homogenization.

ThicknessXYZYZXZXY
0.2 5.085 10 5.071 96 5.096 51 5.050 01 5.049 19 5.057 82 
0.3 5.215 64 5.116 16 5.055 45 5.066 40 5.062 12 5.080 76 
0.4 5.217 24 5.406 39 5.107 65 5.141 07 5.133 29 5.106 66 
ThicknessXYZYZXZXY
0.2 5.085 10 5.071 96 5.096 51 5.050 01 5.049 19 5.057 82 
0.3 5.215 64 5.116 16 5.055 45 5.066 40 5.062 12 5.080 76 
0.4 5.217 24 5.406 39 5.107 65 5.141 07 5.133 29 5.106 66 

The material’s behavior after homogenization is thoroughly explained in Table II, with special attention paid to the variation in Young’s modulus across various thicknesses. A crucial mechanical measure called Young’s modulus reveals a material’s stiffness shown in Fig. 16 and tendency to respond to outside forces. The results show a distinct pattern: the increase in material thickness is accompanied by a perceptible rise in Young’s modulus. This pattern suggests that bigger specimens exhibit a higher deformation resistance, highlighting their superior stiffness compared to slimmer specimens.

FIG. 16.

Stiffness plot of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

FIG. 16.

Stiffness plot of the auxetic metamaterial having strut thickness (a) 0.2 mm, (b) 0.3 mm, and (c) 0.4 mm.

Close modal
TABLE II.

Obtained Young’s modulus after homogenization.

Thickness (mm)Young’s modulus (MPa)
0.2 117.944 
0.3 593.775 
0.4 1699.95 
Thickness (mm)Young’s modulus (MPa)
0.2 117.944 
0.3 593.775 
0.4 1699.95 

The key findings of this study include that the rigidity of the structure tends to rise with increasing lattice strut thickness. This is so that larger loads can be supported without causing severe bending or buckling since thicker struts have a stronger resistance to deformation. The strut thickness is not the only geometric aspect that affects the lattice’s rigidity; the length and placement of the struts are also important factors. Generally speaking, though, thicker struts tend to be stiffer.

In conclusion, this study effectively used homogenization approaches to investigate the mechanical characteristics of auxetic re-entrant-based metamaterials that have characteristics similar to those of Inconel 625. A thorough examination of the displacement patterns and Young’s modulus values over a range of material thicknesses has provided important new information about how these materials react to external forces. In particular, the special qualities of Inconel 625 have improved our comprehension of their behavior. The promise of these metamaterials in applications needing specialized mechanical responses and increased stiffness is demonstrated by findings. From these findings, several intriguing directions for future research are apparent. Analyzing how microstructural changes affect homogenized qualities may help us comprehend the material’s behavior better. In addition, the reliability and legitimacy of the results would be improved by experimental validation of our findings. The range of potential applications would be increased by expanding this study to include dynamic loading situations and taking into account additional material properties in addition to Young’s modulus. Advanced computer methods could also make it possible to analyze complicated structures in greater detail.

The authors have no conflicts to disclose.

Ethics approval is not required.

All the authors have read and agreed to the published version of the manuscript.

Akshansh Mishra: Conceptualization (equal); Data curation (equal); Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Visualization (equal); Writing – original draft (equal). Vijaykumar S Jatti: Data curation (equal); Supervision (equal). Eyob Messele Sefene: Conceptualization (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). Emad Makki: Conceptualization (equal); Investigation (equal); Project administration (equal); Supervision (equal).

The data that support the findings of this study are available within the article.

1.
S.
Mora
,
N. M.
Pugno
, and
D.
Misseroni
, “
3D printed architected lattice structures by material jetting
,”
Mater. Today
59
,
107
132
(
2022
).
2.
A.
Du Plessis
,
N.
Razavi
,
M.
Benedetti
,
S.
Murchio
,
M.
Leary
,
M.
Watson
,
D.
Bhate
, and
F.
Berto
, “
Properties and applications of additively manufactured metallic cellular materials: A review
,”
Prog. Mater. Sci.
125
,
100918
(
2022
).
3.
R. D. K.
Misra
and
K. P.
Misra
, “
Process-structure-biofunctional paradigm in cellular structured implants: An overview and perspective on the synergy between additive manufacturing, bio-mechanical behaviour and biological functions
,”
Artif. Cells, Nanomed., Biotechnol.
51
(
1
),
630
640
(
2023
).
4.
C.
Bregoli
,
F.
Buccino
,
F.
Picca
,
S.
Bagherifard
,
C. A.
Biffi
,
A.
Tuissi
, and
L. M.
Vergani
, “
Additive manufacturing of AISI 316L specimens with distributed inner bone-type cavities: Processability and characterization
,”
IOP Conf. Ser.: Mater. Sci. Eng.
1275
(
1
),
012001
(
2023
).
5.
L.
Greco
,
F.
Buccino
,
Z.
Xu
,
L.
Vergani
,
F.
Berto
,
M.
Guagliano
,
N.
Razavi
, and
S.
Bagherifard
, “
Design and analysis of energy absorbent bioinspired lattice structures
,”
J. Bionic Eng.
20
,
1670
1686
(
2023
).
6.
Z.
Xu
,
I.
La Mendola
,
N.
Razavi
, and
S.
Bagherifard
, “
Additive manufactured triply periodical minimal surface lattice structures with modulated hybrid topology
,”
Eng. Struct.
289
,
116249
(
2023
).
7.
F.
Buccino
,
G.
Martinoia
, and
L. M.
Vergani
, “
Torsion—Resistant structures: A nature addressed solution
,”
Materials
14
(
18
),
5368
(
2021
).
8.
G. I.
Rozvany
, “
A critical review of established methods of structural topology optimization
,”
Struct. Multidiscip. Optim.
37
,
217
237
(
2009
).
9.
O.
Sigmund
and
K.
Maute
, “
Topology optimization approaches: A comparative review
,”
Struct. Multidiscip. Optim.
48
(
6
),
1031
1055
(
2013
).
10.
R. J.
Yang
and
C. J.
Chen
, “
Stress-based topology optimization
,”
Struct. Optim.
12
,
98
105
(
1996
).
11.
O.
Ibhadode
,
Z.
Zhang
,
J.
Sixt
,
K. M.
Nsiempba
,
J.
Orakwe
,
A.
Martinez-Marchese
,
O.
Ero
,
S. I.
Shahabad
,
A.
Bonakdar
, and
E.
Toyserkani
, “
Topology optimization for metal additive manufacturing: Current trends, challenges, and future outlook
,”
Virtual Phys. Prototyping
18
(
1
),
e2181192
(
2023
).
12.
J.
Shen
,
S.
Zhou
,
X.
Huang
, and
Y. M.
Xie
, “
Simple cubic three-dimensional auxetic metamaterials
,”
Phys. Stat. Solidi B
251
(
8
),
1515
1522
(
2014
).
13.
M.
Lei
,
W.
Hong
,
Z.
Zhao
,
C.
Hamel
,
M.
Chen
,
H.
Lu
, and
H. J.
Qi
, “
3D printing of auxetic metamaterials with digitally reprogrammable shape
,”
ACS Appl. Mater. Interfaces
11
(
25
),
22768
22776
(
2019
).
14.
X.
Cheng
,
Y.
Zhang
,
X.
Ren
,
D.
Han
,
W.
Jiang
,
X. G.
Zhang
,
H. C.
Luo
, and
Y. M.
Xie
, “
Design and mechanical characteristics of auxetic metamaterial with tunable stiffness
,”
Int. J. Mech. Sci.
223
,
107286
(
2022
).
15.
M.
Tashkinov
,
A.
Tarasova
,
I.
Vindokurov
, and
V. V.
Silberschmidt
, “
Composites with Re-entrant lattice: Effect of filler on auxetic behaviour
,”
Polymers
15
(
20
),
4076
(
2023
).
16.
R.
Wang
,
Y.
Chen
,
X.
Yan
,
N.
Cong
,
D.
Fang
,
P.
Zhang
,
X.
Liang
, and
W.
Wu
, “
Experimental investigations on the mechanical performances of auxetic metal-ceramic hybrid lattice under quasi-static compression and dynamic ballistic loading
,”
Appl. Sci.
13
(
13
),
7564
(
2023
).
17.
N.
Ghavidelnia
,
M.
Bodaghi
, and
R.
Hedayati
, “
Femur auxetic meta-implants with tuned micromotion distribution
,”
Materials
14
(
1
),
114
(
2021
).