Next Article in Journal
Diatom-Inspired Structural Adaptation According to Mode Shapes: A Study on 3D Structures and Software Tools
Previous Article in Journal
A Mechanical Approach for Comparing Jaws in Fishes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Simulation of the Compressive Strength of Bioinspired Lightweight Structures Manufactured by a Stereolithography 3D Printer

by
Cristina Alía García
,
Álvaro Rodríguez Ortiz
,
José Manuel Arenas Reina
*,
Juan David Cano-Moreno
and
Manuel Gómez Gómez
Escuela Técnica Superior de Ingeniería y Diseño Industrial, Universidad Politécnica de Madrid, 28012 Madrid, Spain
*
Author to whom correspondence should be addressed.
Biomimetics 2024, 9(4), 240; https://doi.org/10.3390/biomimetics9040240
Submission received: 7 March 2024 / Revised: 6 April 2024 / Accepted: 16 April 2024 / Published: 17 April 2024
(This article belongs to the Section Biomimetics of Materials and Structures)

Abstract

:
The use of metamaterials is a good alternative when looking for structures that can withstand compression forces without increasing their weight. In this sense, using nature as a reference can be an appropriate option to design this type of material. Therefore, in this work, a comparative study of a selection of eight representative models of a wide variety of existing solutions, both bioinspired and proposed by various researchers, is presented. These models have been manufactured using stereolithography (SLA) printing, which allows complex geometries to be obtained in a simple way that would be more complicated to achieve by other procedures. Additionally, the manufacturing cost of each model has been determined. The compression tests of the different models have made it possible to evaluate the breaking force and its corresponding deformation. Likewise, a finite element analysis of the manufactured models has been carried out to simulate their behavior under compression, achieving results very similar to those obtained in the experimental tests. In this way, it has been concluded that, among the three-dimensional patterns, the structure called “3D auxetic” is the one that supports the greatest breaking force due to the topographic characteristics of its bar structure. Similarly, among the two-dimensional patterns, the structure called “Auxetic 1”, with a topography based on curves, is capable of supporting the greatest deformation in the compression direction before breaking. Moreover, the highest resistance-force-to-cost ratio has been obtained with a “3D auxetic” structure.

1. Introduction

In the current context, protection against compressive forces has become a key area of study and development in various fields, such as civil engineering, aerospace, and the automotive industry, among others. The goal is to safeguard the lives and integrity of individuals, prevent damage to infrastructures, and minimize the consequences of these forces [1,2].
Furthermore, technological advancements and the constant pursuit of innovative solutions have led to the creation and discovery of new materials with extraordinary properties. Among these materials are metamaterials. Mechanical metamaterials are artificially designed and architected structures that derive their mechanical performance from the unique design of their unit cells [3]. These materials are designed to exhibit characteristics that conventional materials cannot. Due to their versatile behavior, metamaterials are divided into various main groups, including positive Poisson’s ratio (PPR) structures [4], zero Poisson’s ratio (ZPR) structures [5,6,7,8], and negative mechanical metamaterial indices such as auxetics [9,10,11], negative stiffness, and thermal expansion structures [12,13,14,15]. As a result, some of the features they can possess include ultra-high specific strength and a high volumetric modulus-to-shear modulus ratio, among others. These properties make metamaterials an ideal choice for creating structures that absorb and attenuate the effects of compressive forces [16,17,18]. Hence, they are employed in various applications such as sports equipment [19], damping vibration [20], actuators [21] or sound absorption [22,23], providing good compressive strength while maintaining a low weight.
On the other hand, observing and imitating nature is another method for scientists and engineers to find solutions to problems. Animals, plants and microorganisms have evolved over millions of years, developing efficient mechanisms to survive environmental challenges. The strategies organisms use to resist and absorb compressive forces can serve as inspiration for improving or developing new structures. For example, bamboo vascular bundles have been idealized to create honeycomb tubular structures to enhance the energy absorption capacity of vehicle bumpers [24]. The luffa plant inspired the creation of nickel-plated luffa sponge structures with higher energy absorption capacity and toughness compared to other engineering materials like aluminum foam and metal lattice materials [25]. Butterfly wings, with their excellent bending properties, have also been reproduced for engineering applications [26].
In addition, a fundamental tool for development in various scientific and technological fields is finite element analysis (FEA). By using calculations, models, and simulations, FEA predicts the behavior of an object under different physical conditions. Its use can reduce the number of physical prototypes created and experiments conducted, evaluate alternative materials, optimize designs, analyze different solutions quickly and decrease the time and costs associated with design and manufacturing [27].
The advent of 3D printing has revolutionized the manufacturing field, enabling the direct and efficient creation of prototypes and products. This technology allows for the design and manufacture of complex geometries and lightweight structures with high strength, which are challenging to produce using conventional machining and processing techniques [28]. It also enables cost-effective mass customization and reduces the amount of material used during production, as the material is applied only where needed [29,30,31]. In this regard, the use of 3D printing to manufacture metamaterials is widely documented in the literature. Li et al. [32] use SLS technology to fabricate beam lattice structures (BLS) with a negative Poisson’s ratio (NPR) effect by applying sinusoid-shaped beams, conducting experimental compression tests and subsequently validating with finite element analysis. Alomarah et al. [33] employ Multi-Jet Fusion (MJF) to manufacture auxetic structures for compression tests. Hamzehei et al. [34] present new 3D metamaterials with zero Poisson’s ratio for energy absorption applications using SLA 3D printing. Sahin et al. [35] use a new 4D-printing method by rotating the printed axis of bilayers inspired by the large-flowered butterwort (Pinguicula grandiflora). However, when using this type of 3D printing, the associated manufacturing costs are often not considered despite their importance as another variable to consider.
Therefore, the objective of this article is to carry out a technical-economic analysis of bioinspired light structures manufactured by an SLA 3D printer that integrates the results of the compression tests of each support with the evaluation of the cost associated with its manufacture to facilitate the choice of the alternative that best combines the “mechanical resistance force to cost” relationship.
As the number of bioinspired light structures that can be considered is very large, in the present work, eight representative structures have been selected that present a wide functional and formal diversity (five with a two-dimensional pattern and three with a three-dimensional pattern). For each of these structures, their manufacturing cost has been calculated and their mechanical resistance to compression has been determined both experimentally and through finite element simulation using the ABAQUS software (2023 version) [36].
In the selection of the eight representative models, efforts have been made to ensure that there is functional and formal diversity of structures, both with two-dimensional and three-dimensional patterns, that are adequate to resist compression forces. Additionally, the selected models are especially appropriate for manufacturing using 3D printing by stereolithography (SLA).
In this way, for two-dimensional patterns, a first conventional rhombus structure is proposed (model A) that allows behavior with a certain elasticity. Along these lines, the hexagonal structure (model B) presents a somewhat more complex geometry, although with better mechanical behavior. The selection of auxetic geometries allows considering non-conventional structures with a specific mechanical behavior (auxetic models are compressed in their horizontal direction when subjected to a vertical compression stress). In this way, two representative models of auxetic structures have been chosen: models C (cells with smooth curves) and D (stackable cells). Additionally, in order to evaluate whether there is an influence of changes in direction on the propagation of the breakage on the failure stress, a structure based on the Voronoi geometry (model E) characterized by its high degree of irregularity has been selected.
On the other hand, in the case of three-dimensional patterns, a model with a conventional cubic structure and diagonal reinforcements (F models), a three-dimensional auxetic model (G model) and a bioinspired model with a very irregular structure and unpredictable mechanical behavior (model H).

2. Materials and Methods

2.1. Manufacturing of Bioinspired Models

The criterion followed for the design of the structures has been the combination of existing models in the literature (models A, C, D, F, G [37,38,39]) and, in some cases, bio-inspiration (models B, E, and H). Thus, 8 different models have been created, aiming to represent the functional and formal diversity that structural elements subjected to compression forces may exhibit. Furthermore, the various designed models have been grouped into two patterns. On one hand, the structures with a two-dimensional pattern (2D) developed by repeating the unit cells along with the two principal in-plane directions and then extruding in a third direction [40]. On the other hand, the structures with a three-dimensional pattern (3D) are designed by repeating the unit cells along with the three principal directions [41]. The following is an explanation of each of the 8 models divided into 2D and 3D patterns:
  • Two-Dimensional Patterns:
    • Rhombuses. It presents a grid with ribs in two directions perpendicular to each other and rounded corners (Figure 1). This model has a relatively conventional geometry that may be reminiscent of mobile accordion-type structures such as folding fences or architectural elements such as the facade of the Hearst Tower.
    • Hexagonal. It presents a grid similar to the structure of honeycombs or basalt formations (Figure 2). Although the advantages of this type of hexagonal structure are successfully applied to honeycomb boards, in this study, they are subjected to stress in their most unfavorable direction.
    • Auxetic 1. This geometry was chosen following an auxetic pattern typical of literature due to its simplicity, regularity and sense of continuity given by the curved shapes. Said model will be compressed in its horizontal direction when subjected to a vertical compression stress (Figure 3).
    • Auxetic 2. Given the variety of authentic patterns that exist, a second representative model of this type of metamaterials has been selected (Figure 4). The chosen geometry is inspired by stacking systems such as glasses or chairs.
    • Organic Voronoi. This model follows a Voronoi geometry and is characterized by irregularity (Figure 5). The chosen geometry is inspired by the fur of giraffes or the wings of some insects.
  • Three-Dimensional Patterns:
    f.
    Cubic with diagonals. This model follows a conventional geometry following the repetition of cells formed by cylindrical bars in a cubic arrangement with reinforcements on their diagonals (Figure 6). The chosen geometry is inspired by the BCC crystalline atomic structure.
    g.
    Three-dimensional (3D) auxetic. For this model, a design from the literature [42] has been adapted, modifying the bar-shaped structure to a structure with more continuity that gives the sensation of being made up of two-dimensional elements. In this proposal, axial compression generates compression in the other two directions of space (Figure 7), while the model (auxetic 2) only compresses one direction indirectly.
    h.
    Organic (irregular mesh). To design this model, the base pattern created for model F, which is slightly elongated in the vertical direction, was used as a starting point. Said structure would be figuratively wrapped in an organic and irregular mesh, inspired by how a spider wraps prey that falls into its web or bird nests (Figure 8).
Once the models are selected by combining literature and bioinspiration, the 8 designs are modeled through different parametric design tools: Auto-desk Inventor, Autodesk Fusion 360, Autodesk Meshmixer and Chitubox. Figure 9 shows the modeling result of each of the designs. Furthermore, in order to compare the models with each other, all structures share maximum dimensions of 50 × 50 × 102.67 [mm] and a net volume of 77.88 mL ± 1%.

2.2. Equipment

Once the bioinspired structures are selected, they are manufactured using a first-generation stereolithography (SLA) 3D printer (ELEGOO, Shenzhen, China), which uses a 4.72 × 2.68 square inch LCD screen and a workspace height of 6.1 inches. Figure 10 shows the main components of the SLA printer.
The printing parameters used to manufacture the different structures have been kept the same for all designs in order to allow for a comparison of all models. Table 1 shows the printing parameters.
For the tests, a model TN-MD machine (HOYTOM, S.L., Bilbao, Spain), motorized with automatic control via a computer, was used. Its capacity is 200 kN, the piston stroke length is 125 mm, and the displacement rate was fixed at 2 mm/min. Three tests have been carried out on each of the designed models.

2.3. Materials

The material used was a standard photopolymer resin (ELEGOO, Shenzhen, China) for SLA printing [43]. This resin has high fluidity and low viscosity, fast photocuring ability, high accuracy printing, good prototyping effect, easy cleaning and curing and wide applications. Once the specimen is printed it must be washed using an isopropyl alcohol solution for a short period of time, as described by the provider. This resin is commonly used to prototype parts that would be made from polypropylene (PP). The use of this resin is due to the fact that its mechanical properties [44] seem to be suitable for use in elements that will be subjected to large deformations or impacts, which is the potential application. Table 2 shows the technical characteristics of the resin used.
Table 3 specifies, for each of the eight structure models manufactured, the main dimensions of each lattice, the diameter of the reinforcement, the mass and the density.

2.4. Manufacturing Costs

The determination of the manufacturing cost in 3D printing manufacturing processes is a relevant element in evaluating the technical-economic viability of the parts obtained by these procedures [45]. In this way, the model that best combines mechanical performance and adaptation to the SLA manufacturing process can be selected.
The manufacturing cost has been divided into the cost of the material used in the model, cost of labor and cost of operating the workplace. The SLA printer software (Elegoo Mar 5 version) itself allows you to determine the mass of resin used and the printing time. This software also provides the material cost (including the material used as a support, which, in this case, is the same resin) when the type of resin and the cost per unit of mass (EUR/kg) are indicated. The cost of labor is calculated as a product of the printing time times the weighted salary of the operator for this type of process. The weighted day considers that the operator only dedicates a small fraction of the model processing time to attending to the 3D printer (mainly for starting up the machine, periodic supervision and performing post-curing), dedicating the rest of the time to other processes and activities. In this way, the weighted salary is obtained as a product of the standard salary for the aforementioned fraction of time dedicated to the model development process (in our case, this weighted salary is 5.30 EUR/hour). The most relevant factors of the operating cost of the 3D printer are the electrical energy consumption and the amortization of the machine. Considering the energy consumed per hour of the machine and that the expected amortization period is four years, the operating cost per hour of the machine has been estimated at 0.035 EUR/hour. The operating cost for processing each model is obtained as a product of the operating cost/hour times the corresponding printing time. Table 4 shows the detailed manufacturing costs (material, labor and operating costs) for each of the eight models manufactured by SLA.

2.5. Methodology of Finite Element Analysis (FEA)

The finite element analysis of the most representative models (A, B, C, D, E, F, G) has been carried out. To do this, the experimental compression test was considered, where the specimen was supported on its lower base while the piston of the testing machine came into contact with a downward movement on its upper face. To replicate these boundary conditions in the FEM model, an embedment condition has been placed on the lower base, and pressure has been applied to the upper face (Figure 11).
The pressure value that has been considered for each model is the equivalent of the average failure force recorded in the experimental tests.
The mesh used is the result of a compromise between the available computing power and the maximum quality of the mesh. The quality of the mesh is determined by the size, distribution and type. In terms of size, we chose one that was equivalent to the minimum thickness that existed in the model. The distribution of the elements can be structured or free. ABAQUS allows a structured distribution when the models have regular shapes. In this case, this technique was available for some models, so for the analysis, it was decided to use the free distribution in all models. When using a free distribution, the elements have to be tetrahedral, leaving it to be decided whether they have to be linear or quadratic. The difference between the two is the number of existing nodes per element, being 4 for the linear ones and 10 for the quadratic ones. The decision was made to use quadratic elements, which, although they require more calculation resources, result in greater precision. In Figure 12, you can see the result of the meshing for model A, where C3D10-type elements have been used.

3. Results

3.1. Experiment Results

Once the tests carried out and the models used have been described, this section will show both the results obtained and their analysis. The different models used will allow us to study how each of the patterns behaves when subjected to compression stresses.
Figure 13 shows the force–displacement curves of the specimens tested in compression. It can be seen that in all cases, except in model D (auxetic 2), the force increases to a maximum value where breakage occurs. In the case of model D, due to its layer configuration, it makes the breakage different, having an intermediate partial breakage before the total breakage. Furthermore, differences in the maximum supported force are observed between the different models and an increase in rigidity is observed as a function of deformation. This aspect is apparently linked to local geometric nonlinearity effects and is very different from that usually found in compact materials. Thus, in structures in which the applied loads produce deformations of great magnitude, such that the hypothesis that the final deformed position coincides with the initial position cannot be accepted. This fact makes the response of the structure highly nonlinear, and it is common to call this phenomenon geometric nonlinearity [46]. Furthermore, this nonlinearity affects the mode of failure that occurs and the speed at which it occurs, causing total failure, as shown in the graphs.
After carrying out the mechanical compression tests, the experimental results shown in Table 5 have been obtained.
Figure 14 shows a correlation diagram that relates the failure stress (in %) to the compression failure load for the eight models considered. In the case of two-dimensional patterns, which offer anisotropic properties intended for unidirectional loads, it is observed that model E (organic Voronoi) has the highest failure load while the auxetic models (D and C) have the lowest load they support. Regarding the breaking deformation, all models with a two-dimensional pattern present a deformation between 3.18% and 4.83%.
In the case of three-dimensional patterns intended to support 3D loads, it is observed that model G (3D Auxetic) stands out for its high breaking load (13.39 kN), while model H (organic and irreg. mesh) is the one that presents the least load. The fracture stress values for the three-dimensional pattern models are between 4.02% and 5.01%.
An analysis of the breaking force and the manufacturing cost of each of the models has also been conducted (Table 6). In this way, it is observed that model G (Three-dimensional auxetic) is the one with the highest manufacturing cost. The rest of the models, except for the H model (organic with irregular mesh), whose cost is also high, have a similar cost, although there is a difference in the breaking forces that they are capable of withstanding. This comparison is very practical when the economic part is relevant in decision-making. If we compare the force/cost ratio, model G presents the best value, eight times higher than the D model and five times higher than the C model.
On the other hand, the type of breakage observed in each model is analyzed. The fracture line followed by each model is shown in red. In the case of models A (rhombuses) and B (hexagonal), the angle of failure is 45°, as are materials with ductile behavior (Figure 15).
However, models C (auxetic 1), D (auxetic 2), E (organic Voronoi) and F (cubic with diagonals) have a fracture angle of 90° typical of brittle materials (Figure 16).
On the other hand, in models G (three-dimensional auxetic) and H (organic with irregular mesh), the breakage begins in one corner of the model, and then the propagation begins (Figure 17).
Table 7 shows a summary of the breakage types obtained. Therefore, there are very different fracture behaviors depending on the model design.

3.2. Results of Finite Element Analysis (FEA)

This section describes the results obtained when applying finite element analysis (FEA) to the structural models considered in Section 2.5.
With these conditions, the results of the finite element analysis are shown. The stress field corresponding to the Von Misses stresses (units in Pascals) is shown for both the complete specimen and the middle section for each of the analyzed models (Figure 18). According to the manufacturer, in its technical sheet, the compressive strength is between 36 and 53 MPa, and the bending strength is between 59 and 70 MPa. The stress state of the test specimen is a mix of compression and traction stress states depending on the areas. The failure criterion is an average value of the minimum and maximum failure stresses provided by the manufacturer ((70 + 36)/2 = 53 ~ 50 MPa). With this criterion, the test specimen collapses when any area reaches 50 MPa von Misses stresses. The visualization of the stresses in the ABAQUS models was set, placing the visualization between 0 MPa (blue color) and 50 MPa (red color).
The models have been loaded by using a linear increase in force. Figure 18 shows the moment in which the areas that reach a stress of 50 MPa are relevant. That is, from that moment on, it is considered that the material has reached its limit of structural resistance, and consequently, the collapse of the structure occurs. Additionally, in the figure, there are also small grey areas that indicate that in these areas, the maximum stress (50 MPa) has been exceeded, which corroborates that this instant is considered the moment in which the collapse of the structure.
Table 8 shows the comparison between the effective failure stress values of each model and those obtained by the finite element analysis. With these data, it can be verified that there is a coherence and high correlation between the experimental tests and the finite element models. This will allow us to make more complex models, taking the failure criterion as the moment at which a stress of 50 MPa is reached.

4. Discussion

In the compression test of the analyzed models, it was observed that the force–displacement curves are linear with a single loading cycle and abrupt failure in all cases except in model D (auxetic 2), which presents two loading cycles and one more gradual breakage.
In view of the experimental tests carried out, it can be stated that design G (three-dimensional auxetic) is the model that supports the greatest breaking force and is also the most expensive to manufacture, as well as that design C (auxetic 1) is capable of withstanding the greatest deformation in the compression direction before breaking, having an average manufacturing cost in relation to the rest of the models. Additionally, design D (auxetic 2) has shown that in its breaking process, the fragments that fractured were not projected toward the outside of the piece, but rather an implosion phenomenon took place. The regions where a fracture occurred were not a starting point for the propagation of said fracture. However, in other designs such as A (rhombuses), B (hexagonal), E (Voronoi organic) or F (cubic with diagonals), a direction of fracture propagation is clearly seen, design D (auxetic 2) tends to show different decentralized fracture points.
Furthermore, model D (auxetic 2) maintains its integrity and stability despite suffering numerous partial fractures even when there is a continuous fracture surface that splits the model in two. In this case, the piece is divided into several parts and has zero resistance to traction efforts.
On the other hand, the H design pieces (organic with irregular mesh) show a very interesting behavior that can be highly useful in structures that need to ensure their integrity for safety reasons. This model is capable of resisting compression forces of more than 4000 N until it compresses approximately 5%. From this point on, a partial fracture allows the structure to continue compressing while resisting stresses of up to 2000 N, according to experimental results. In this way, the change in behavior warns of a partial rupture in a non-catastrophic way.
The finite element analysis of the models analyzed in the simulation of the compression test has made it possible to determine the breaking stresses and the failure mode of the different models. The high degree of similarity between the results of the experimental test and those obtained in the FEM analysis allows the simulation carried out with the considered hypotheses to be validated. In particular, the appropriateness of using the moment at which a stress of 50 MPa is reached as the model’s failure criterion has been confirmed.
Additionally, it has been observed that the initial failure of the models varies depending on the design. Similar behaviors are observed for each design, both in the model and in the experimental tests. Once the test failure occurs, a rapid progression of the breakage appears, while these dynamic effects have not been introduced in the model.

5. Conclusions

The following points provide the main conclusions of the presented research:
  • Up to eight different models have been designed and manufactured, with two-dimensional and three-dimensional patterns of structures that could work as compression stress dissipative layers, with some being based on models found in the literature and others being bioinspired. The objective has been to achieve lightweight structures resistant to compressive stresses in critical service conditions where it is necessary to guarantee the safety of people and property.
  • The manufacturing of the models has been carried out using SLA 3D printing technology. This system has been very useful in the creation of prototypes, allowing the manufacturing of complex geometries quickly and efficiently. It has also made it possible to evaluate and verify the functionality of the designed models.
  • The experimental results of the compression tests have shown that the choice of one model or another will depend on the specific needs of each situation. Geometry is a determining factor in the behavior of structures. In situations where greater deformation is sought, a more flexible structure, such as model C (auxetic 1), will be required. On the other hand, given a greater need for stability and load, a more rigid and resistant structure, such as the G model (three-dimensional auxetic), will be preferable. Moreover, analyzing costs, the G model (three-dimensional auxetic) is the most expensive but also presents the best profitability since it has the highest force/cost ratio.
  • On the other hand, it has been experimentally proven that the types of breakage are also different depending on the configuration they present. An assessment of the manufacturing cost of each of the models has also been made, considering all the factors that affect them. In this way, the relationship between strength at breakage and manufacturing cost is analyzed, which can be very relevant when deciding on one model or another for a given application.
  • Furthermore, through finite element analysis, the behavior of the designed structures has also been evaluated, allowing their behavior and stress distribution to be predicted, which are key characteristics for the optimization of designs. The results obtained have allowed us to validate the simulation carried out with the considered hypotheses. With this, the ability to create and analyze more complex models in less time and optimize the selection of the most appropriate alternative has been achieved.

Author Contributions

Conceptualization, C.A.G. and M.G.G.; methodology, C.A.G., M.G.G. and J.M.A.R.; software, Á.R.O.; validation, C.A.G., Á.R.O. and J.D.C.-M.; formal analysis, all authors; investigation, all authors; resources, J.M.A.R.; data curation, M.G.G., C.A.G. and Á.R.O.; writing—original draft preparation, all authors; writing—review and editing, J.D.C.-M.; supervision, J.M.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Dataset is available upon request from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mansoori, H.; Hamzehei, R.; Dariushi, S. Crashworthiness Analysis of Cylindrical Tubes with Coupling Effects under Quasi-Static Axial Loading: An Experimental and Numerical Study. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2022, 236, 647–662. [Google Scholar] [CrossRef]
  2. Bejjani, R.; Odelros, S.; Öhman, S.; Collin, M. Shift of Wear Balance Acting on CVD Textured Coatings and Relation to Workpiece Materials. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2021, 235, 114–128. [Google Scholar] [CrossRef]
  3. Rezaei, S.; Kadkhodapour, J.; Hamzehei, R.; Taherkhani, B.; Anaraki, A.P.; Dariushi, S. Design and Modeling of the 2D Auxetic Metamaterials with Hyperelastic Properties Using Topology Optimization Approach. Photonics Nanostr. Fundam. Appl. 2021, 43, 100868. [Google Scholar] [CrossRef]
  4. Ha, N.S.; Pham, T.M.; Tran, T.T.; Hao, H.; Lu, G. Mechanical Properties and Energy Absorption of Bio-Inspired Hierarchical Circular Honeycomb. Compos. Part B Eng. 2022, 236, 109818. [Google Scholar] [CrossRef]
  5. Chen, Y.; Zheng, B.-B.; Fu, M.-H.; Lan, L.-H.; Zhang, W.-Z. Doubly Unusual 3D Lattice Honeycomb Displaying Simultaneous Negative and Zero Poisson’s Ratio Properties. Smart Mater. Struct. 2018, 27, 045003. [Google Scholar] [CrossRef]
  6. Zolfagharian, A.; Bodaghi, M.; Hamzehei, R.; Parr, L.; Fard, M.; Rolfe, B.F. 3D-Printed Programmable Mechanical Metamaterials for Vibration Isolation and Buckling Control. Sustainability 2022, 14, 6831. [Google Scholar] [CrossRef]
  7. Huang, J.; Gong, X.; Zhang, Q.; Scarpa, F.; Liu, Y.; Leng, J. In-Plane Mechanics of a Novel Zero Poisson’s Ratio Honeycomb Core. Compos. Part B Eng. 2016, 89, 67–76. [Google Scholar] [CrossRef]
  8. Hamzehei, R.; Zolfagharian, A.; Dariushi, S.; Bodaghi, M. 3D-Printed Bio-Inspired Zero Poisson’s Ratio Graded Metamaterials with High Energy Absorption Performance. Smart Mater. Struct. 2022, 31, 035001. [Google Scholar] [CrossRef]
  9. Li, D.; Yin, J.; Dong, L.; Lakes, R.S. Strong Re-Entrant Cellular Structures with Negative Poisson’s Ratio. J. Mater Sci. 2018, 53, 3493–3499. [Google Scholar] [CrossRef]
  10. Lv, W.; Li, D.; Dong, L. Study on Blast Resistance of a Composite Sandwich Panel with Isotropic Foam Core with Negative Poisson’s Ratio. Int. J. Mech. Sci. 2021, 191, 106105. [Google Scholar] [CrossRef]
  11. Li, D.; Shen, G. Study on Mechanical Properties of an Isotropic Negative Poisson’s Ratio Voronoi Foam and Its Foam-Filled Tube. Smart Mater. Struct. 2022, 31, 065017. [Google Scholar] [CrossRef]
  12. Morris, C.; Bekker, L.; Spadaccini, C.; Haberman, M.; Seepersad, C. Tunable Mechanical Metamaterial with Constrained Negative Stiffness for Improved Quasi-Static and Dynamic Energy Dissipation. Adv. Eng. Mater. 2019, 21, 1900163. [Google Scholar] [CrossRef]
  13. Wang, Z.; Luan, C.; Liao, G.; Liu, J.; Yao, X.; Fu, J. Progress in Auxetic Mechanical Metamaterials: Structures, Characteristics, Manufacturing Methods, and Applications. Adv. Eng. Mater. 2020, 22, 2000312. [Google Scholar] [CrossRef]
  14. Wang, H.; Xiao, S.-H.; Zhang, C. Novel Planar Auxetic Metamaterial Perforated with Orthogonally Aligned Oval-Shaped Holes and Machine Learning Solutions. Adv. Eng. Mater. 2021, 23, 2100102. [Google Scholar] [CrossRef]
  15. Hamzehei, R.; Kadkhodapour, J.; Anaraki, A.P.; Rezaei, S.; Dariushi, S.; Rezadoust, A.M. Octagonal Auxetic Metamaterials with Hyperelastic Properties for Large Compressive Deformation. Int. J. Mech. Sci. 2018, 145, 96–105. [Google Scholar] [CrossRef]
  16. Surjadi, J.U.; Gao, L.; Du, H.; Li, X.; Xiong, X.; Fang, N.X.; Lu, Y. Mechanical Metamaterials and Their Engineering Applications. Adv. Eng. Mater. 2019, 21, 1800864. [Google Scholar] [CrossRef]
  17. Hamzehei, R.; Rezaei, S.; Kadkhodapour, J.; Anaraki, A.P.; Mahmoudi, A. 2D Triangular Anti-Trichiral Structures and Auxetic Stents with Symmetric Shrinkage Behavior and High Energy Absorption. Mech. Mater. 2020, 142, 103291. [Google Scholar] [CrossRef]
  18. Imbalzano, G.; Tran, P.; Ngo, T.D.; Lee, P.V.S. A Numerical Study of Auxetic Composite Panels under Blast Loadings. Compos. Struct. 2016, 135, 339–352. [Google Scholar] [CrossRef]
  19. Sanami, M.; Ravirala, N.; Alderson, K.; Alderson, A. Auxetic Materials for Sports Applications. Procedia Eng. 2014, 72, 453–458. [Google Scholar] [CrossRef]
  20. Ma, Y.; Scarpa, F.; Zhang, D.; Zhu, B.; Chen, L.; Hong, J. A Nonlinear Auxetic Structural Vibration Damper with Metal Rubber Particles. Smart Mater. Struct. 2013, 22, 084012. [Google Scholar] [CrossRef]
  21. Imbalzano, G.; Tran, P.; Ngo, T.D.; Lee, P.V. Three-Dimensional Modelling of Auxetic Sandwich Panels for Localised Impact Resistance. J. Sandw. Struct. Mater. 2017, 19, 291–316. [Google Scholar] [CrossRef]
  22. Rossiter, J.; Takashima, K.; Scarpa, F.; Walters, P.; Mukai, T. Shape Memory Polymer Hexachiral Auxetic Structures with Tunable Stiffness. Smart Mater. Struct. 2014, 23, 045007. [Google Scholar] [CrossRef]
  23. Scarpa, F.; Bullough, W.A.; Lumley, P. Trends in Acoustic Properties of Iron Particle Seeded Auxetic Polyurethane Foam. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2004, 218, 241–244. [Google Scholar] [CrossRef]
  24. Hu, D.; Wang, Y.; Song, B.; Dang, L.; Zhang, Z. Energy-Absorption Characteristics of a Bionic Honeycomb Tubular Nested Structure Inspired by Bamboo under Axial Crushing. Compos. Part B Eng. 2019, 162, 21–32. [Google Scholar] [CrossRef]
  25. Yin, S.; Wang, H.; Li, J.; Ritchie, R.O.; Xu, J. Light but Tough Bio-Inherited Materials: Luffa Sponge Based Nickel-Plated Composites. J. Mech. Behav. Biomed. Mater. 2019, 94, 10–18. [Google Scholar] [CrossRef] [PubMed]
  26. Pelanconi, M.; Ortona, A. Nature-Inspired, Ultra-Lightweight Structures with Gyroid Cores Produced by Additive Manufacturing and Reinforced by Unidirectional Carbon Fiber Ribs. Materials 2019, 12, 4134. [Google Scholar] [CrossRef] [PubMed]
  27. Chakraverty, S.; Pradhan, K.K. Computational Structural Mechanics: Static and Dynamic Behaviors, 1st ed.; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  28. Paris, H.; Mokhtarian, H.; Coatanéa, E.; Museau, M.; Ituarte, I.F. Comparative Environmental Impacts of Additive and Subtractive Manufacturing Technologies. CIRP Ann. 2016, 65, 29–32. [Google Scholar] [CrossRef]
  29. Håkansson, K.M.O.; Henriksson, I.C.; de la Peña Vázquez, C.; Kuzmenko, V.; Markstedt, K.; Enoksson, P.; Gatenholm, P. Solidification of 3D Printed Nanofibril Hydrogels into Functional 3D Cellulose Structures. Adv. Mater. Technol. 2016, 1, 1600096. [Google Scholar] [CrossRef]
  30. Guvendiren, M.; Molde, J.; Soares, R.M.D.; Kohn, J. Designing Biomaterials for 3D Printing. ACS Biomater. Sci. Eng. 2016, 2, 1679–1693. [Google Scholar] [CrossRef]
  31. Ngo, T.D.; Kashani, A.; Imbalzano, G.; Nguyen, K.T.Q.; Hui, D. Additive Manufacturing (3D Printing): A Review of Materials, Methods, Applications and Challenges. Compos. Part B Eng. 2018, 143, 172–196. [Google Scholar] [CrossRef]
  32. Li, S.; Hassanin, H.; Attallah, M.M.; Adkins, N.J.E.; Essa, K. The Development of TiNi-Based Negative Poisson’s Ratio Structure Using Selective Laser Melting. Acta Mater. 2016, 105, 75–83. [Google Scholar] [CrossRef]
  33. Alomarah, A.; Ruan, D.; Masood, S.; Gao, Z. Compressive Properties of a Novel Additively Manufactured 3D Auxetic Structure. Smart Mater. Struct. 2019, 28, 085019. [Google Scholar] [CrossRef]
  34. Hamzehei, R.; Serjouei, A.; Wu, N.; Zolfagharian, A.; Bodaghi, M. 4D Metamaterials with Zero Poisson’s Ratio, Shape Recovery, and Energy Absorption Features. Adv. Eng. Mater. 2022, 24, 2200656. [Google Scholar] [CrossRef]
  35. Sahin, E.S.; Cheng, T.; Wood, D.; Tahouni, Y.; Poppinga, S.; Thielen, M.; Speck, T.; Menges, A. Cross-Sectional 4D-Printing: Upscaling Self-Shaping Structures with Differentiated Material Properties Inspired by the Large-Flowered Butterwort (Pinguicula Grandiflora). Biomimetics 2023, 8, 233. [Google Scholar] [CrossRef] [PubMed]
  36. ABAQUS HTML Documentation 6.7.0.0. Available online: https://abaqus-html-documentation.updatestar.com (accessed on 6 March 2024).
  37. Kolken, H.M.A.; Zadpoor, A.A. Auxetic Mechanical Metamaterials. RSC Adv. 2017, 7, 5111–5129. [Google Scholar] [CrossRef]
  38. Kim, Y.; Son, K.H.; Lee, J.W. Auxetic Structures for Tissue Engineering Scaffolds and Biomedical Devices. Materials 2021, 14, 6821. [Google Scholar] [CrossRef] [PubMed]
  39. Wu, W.; Hu, W.; Qian, G.; Liao, H.; Xu, X.; Berto, F. Mechanical Design and Multifunctional Applications of Chiral Mechanical Metamaterials: A Review. Mater. Des. 2019, 180, 107950. [Google Scholar] [CrossRef]
  40. Zhang, W.; Zhao, S.; Scarpa, F.; Wang, J.; Sun, R. In-Plane Mechanical Behavior of Novel Auxetic Hybrid Metamaterials. Thin-Walled Struct. 2021, 159, 107191. [Google Scholar] [CrossRef]
  41. Valle, R.; Pincheira, G.; Tuninetti, V. Design of an Auxetic Cellular Structure with Different Elastic Properties in Its Three Orthogonal Directions. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2021, 235, 1341–1350. [Google Scholar] [CrossRef]
  42. Balabanov, S.V.; Makogon, A.I.; Sychov, M.M.; Gravit, M.V.; Kurakin, M.K. Mechanical Properties of 3D Printed Cellular Structures with Topology of Triply Periodic Minimal Surfaces. Mater. Today Proc. 2020, 30, 439–442. [Google Scholar] [CrossRef]
  43. Palucci Rosa, R.; Rosace, G. Nanomaterials for 3D Printing of Polymers via Stereolithography: Concept, Technologies, and Applications. Macromol. Mater. Eng. 2021, 306, 2100345. [Google Scholar] [CrossRef]
  44. Martín-Montal, J.; Pernas-Sánchez, J.; Varas, D. Experimental Characterization Framework for SLA Additive Manufacturing Materials. Polymers 2021, 13, 1147. [Google Scholar] [CrossRef] [PubMed]
  45. Arenas, J.M.; Alía, C.; Blaya, F.; Sanz, A. Multi-Criteria Selection of Structural Adhesives to Bond ABS Parts Obtained by Rapid Prototyping. Int. J. Adhes. Adhes. 2012, 33, 67–74. [Google Scholar] [CrossRef]
  46. Celigüeta Lizarza, J.T. Análisis de Estructuras Con No Linealidad Geométrica; Universidad Pública de Navarra = Nafarroako Unibertsitate Publikoa: Pamplona, Spain, 2009; ISBN 978-84-921970-3-3. [Google Scholar]
Figure 1. Rhombuses model (A) and Hearst Tower (B).
Figure 1. Rhombuses model (A) and Hearst Tower (B).
Biomimetics 09 00240 g001
Figure 2. Hexagonal model (A) and basalt formations (B).
Figure 2. Hexagonal model (A) and basalt formations (B).
Biomimetics 09 00240 g002
Figure 3. Auxetic 1 model (A) and auxetic pattern study (B).
Figure 3. Auxetic 1 model (A) and auxetic pattern study (B).
Biomimetics 09 00240 g003
Figure 4. Auxetic 2 model (A) and stacked chair system (B).
Figure 4. Auxetic 2 model (A) and stacked chair system (B).
Biomimetics 09 00240 g004
Figure 5. Organic Voronoi model (A) and insect wing (B).
Figure 5. Organic Voronoi model (A) and insect wing (B).
Biomimetics 09 00240 g005
Figure 6. Cubic with diagonals model (A) and BCC structure (B).
Figure 6. Cubic with diagonals model (A) and BCC structure (B).
Biomimetics 09 00240 g006
Figure 7. Three-dimensional auxetic model (A) and mechanical properties of 3D double-U structures (B) [39].
Figure 7. Three-dimensional auxetic model (A) and mechanical properties of 3D double-U structures (B) [39].
Biomimetics 09 00240 g007
Figure 8. Organic model with irregular mesh (A) and bird nest (B).
Figure 8. Organic model with irregular mesh (A) and bird nest (B).
Biomimetics 09 00240 g008
Figure 9. Designed models: (A) rhombuses, (B) hexagonal, (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi, (F) cubic with diagonals, (G) three-dimensional auxetic and (H) organic with irregular mesh.
Figure 9. Designed models: (A) rhombuses, (B) hexagonal, (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi, (F) cubic with diagonals, (G) three-dimensional auxetic and (H) organic with irregular mesh.
Biomimetics 09 00240 g009
Figure 10. Main components of the printer: (1) main body, (2) protective cover with UV filter, (3) bucket and (4) manufacturing platform.
Figure 10. Main components of the printer: (1) main body, (2) protective cover with UV filter, (3) bucket and (4) manufacturing platform.
Biomimetics 09 00240 g010
Figure 11. Boundary conditions in the FEM model A (rhombuses).
Figure 11. Boundary conditions in the FEM model A (rhombuses).
Biomimetics 09 00240 g011
Figure 12. Mesh in the FEM model A (rhombuses).
Figure 12. Mesh in the FEM model A (rhombuses).
Biomimetics 09 00240 g012
Figure 13. Force–displacement curves for proposed structure models: (A) rhombuses, (B) hexagonal, (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi, (F) cubic with diagonals, (G) three-dimensional auxetic and (H) organic with irregular mesh.
Figure 13. Force–displacement curves for proposed structure models: (A) rhombuses, (B) hexagonal, (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi, (F) cubic with diagonals, (G) three-dimensional auxetic and (H) organic with irregular mesh.
Biomimetics 09 00240 g013aBiomimetics 09 00240 g013b
Figure 14. Comparative diagram between breaking force and breaking deformation for the eight models considered (A to G).
Figure 14. Comparative diagram between breaking force and breaking deformation for the eight models considered (A to G).
Biomimetics 09 00240 g014
Figure 15. Breakage type of designs (A) (rhombuses) and (B) (hexagonal).
Figure 15. Breakage type of designs (A) (rhombuses) and (B) (hexagonal).
Biomimetics 09 00240 g015
Figure 16. Breakage type of designs: (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi and (F) cubic with diagonals.
Figure 16. Breakage type of designs: (C) auxetic 1, (D) auxetic 2, (E) organic Voronoi and (F) cubic with diagonals.
Biomimetics 09 00240 g016
Figure 17. Breakage type of designs: (G) 3D auxetic and (H) organic with irregular mesh.
Figure 17. Breakage type of designs: (G) 3D auxetic and (H) organic with irregular mesh.
Biomimetics 09 00240 g017
Figure 18. Stress field of finite element models simulated with ABAQUS.
Figure 18. Stress field of finite element models simulated with ABAQUS.
Biomimetics 09 00240 g018
Table 1. SLA 3D printing parameters.
Table 1. SLA 3D printing parameters.
Layer Height0.05 mm
Bottom Layer Count8
Transition TypeLinear
Exposure Time6 s
Bottom Exposure Time50 s
Lifting Distance5 mm
Lifting Speed65 mm/min
Retract Speed150 mm/min
Table 2. Resin technical properties.
Table 2. Resin technical properties.
Hardness/Contraction84D/7.1%
Viscosity (25 °C)150–200 mPa.s
Density (liquid)1.100 g/cm3
Density (solid)1.195 g/cm3
Bending resistance59–70 MPa
Compressive strength36–53 MPa
Elongation at break14.2
Price36.99 EUR/kg
Table 3. Characteristics of the eight manufactured structure models.
Table 3. Characteristics of the eight manufactured structure models.
Pattern TypeModelLattice Dimensions (mm)Truss Diameter
(mm)
Mass (g)Density (g/cm3)
2D(A) Rhombuses12.5 × 12.51.6930.372
(B) Hexagonal10.2 × 9.11.7930.372
(C) Auxetic 119.5 × 13.91.4940.376
(D) Auxetic 2 18.9 × 12.10.9950.380
(E) Organic VoronoiIrregularIrregular1020.408
3D(F) Cubic with diagonals12.4 × 12.42.7920.368
(G) 3D auxetic17.7 × 9.23.2940.376
(H) Organic–irregular meshIrregularIrregular930.372
Table 4. Manufacturing cost of the eight models printed by SLA.
Table 4. Manufacturing cost of the eight models printed by SLA.
Pattern TypeModelMass (g)Printing Time (h)Material Cost
(EUR)
Labor Cost
(EUR)
Operating Cost
(EUR)
Manufacturing Cost (EUR)
2D(A) Rhombuses9310.23.4454.340.3658.14
(B) Hexagonal939.823.4452.370.3456.15
(C) Auxetic 1949.823.4852.390.3456.21
(D) Auxetic 2 9510.33.5154.90.3658.77
(E) Organic Voronoi10210.533.7756.190.3760.33
3D(F) Cubic with diagonals9210.153.4054.080.3557.83
(G) 3D auxetic9415.323.4880.970.5484.99
(H) Organic with irregular mesh9314.83.4478.250.5282.21
Table 5. Experimental outputs of 8 proposed models.
Table 5. Experimental outputs of 8 proposed models.
PatternsModelsMean Breaking Force
(kN)
Mean Break Deformation
(%)
Mean
Effective Breaking Stress
(MPa)
Mean Young’s
Modulus
(MPa)
Absolute Difference
(kN)
Relative
Difference
(%)
2D(A) Rhombuses2.803.811.12796.011.45107
(B) Hexagonal4.583.871.831279.433.23239
(C) Auxetic 11.494.830.60401.950.1410
(D) Auxetic 2 1.353.200.54221--
(E) Organic Voronoi4.953.181.981588.473.6267
3D(F) Cubic and diagonals5.774.052.3116674.42327
(G) 3D auxetic13.394.025.363389.312.04892
(H) Organic and irreg. mesh4.345.011.741410.632.99221
Table 6. Comparative analysis between breaking force and manufacturing costs.
Table 6. Comparative analysis between breaking force and manufacturing costs.
PatternsModelsMean Breaking Force
(kN)
Manufacturing
Cost
(EUR)
Force/Cost
Ratio
(kN/EUR)
Absolute Difference
(kN/EUR)
Relative Difference
(%)
2D(A) Rhombuses2.8058.140.050.03150
(B) Hexagonal4.5856.150.080.06300
(C) Auxetic 11.4956.210.030.0150
(D) Auxetic 2 1.3558.770.02--
(E) Organic Voronoi4.9560.330.080.06300
3D(F) Cubic and diagonals5.7757.830.10.08400
(G) 3D auxetic13.3984.990.160.14700
(H) Organic and irreg. mesh4.3482.210.050.03150
Table 7. Type of breakage of each model.
Table 7. Type of breakage of each model.
PatternsModelsType of Breakage
2D(A) RhombusesDuctile
(B) HexagonalDuctile
(C) Auxetic 1Fragile
(D) Auxetic 2 Fragile
(E) Organic VoronoiFragile
3D(F) Cubic and diagonalsFragile
(G) 3D auxeticOther
(H) Organic and irreg. meshOther
Table 8. Comparison of results between experimental tests and finite element model (FEM) outputs.
Table 8. Comparison of results between experimental tests and finite element model (FEM) outputs.
PatternModelsEffective
Breaking
Stress
(MPa)
Breaking Stress FEM
(MPa)
Absolut Value
(MPa)
Relative Value
(%)
2D(A) Rhombuses1.121.060.065.36
(B) Hexagonal1.831.680.158.19
(C) Auxetic 10.600.650.058.33
(D) Auxetic 2 0.540.60.0611.11
(E) Organic Voronoi1.981.890.094.54
3D(F) Cubic and diagonals2.312.210.14.33
(G) 3D auxetic5.364.840.529.70
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alía García, C.; Rodríguez Ortiz, Á.; Arenas Reina, J.M.; Cano-Moreno, J.D.; Gómez Gómez, M. Analysis and Simulation of the Compressive Strength of Bioinspired Lightweight Structures Manufactured by a Stereolithography 3D Printer. Biomimetics 2024, 9, 240. https://doi.org/10.3390/biomimetics9040240

AMA Style

Alía García C, Rodríguez Ortiz Á, Arenas Reina JM, Cano-Moreno JD, Gómez Gómez M. Analysis and Simulation of the Compressive Strength of Bioinspired Lightweight Structures Manufactured by a Stereolithography 3D Printer. Biomimetics. 2024; 9(4):240. https://doi.org/10.3390/biomimetics9040240

Chicago/Turabian Style

Alía García, Cristina, Álvaro Rodríguez Ortiz, José Manuel Arenas Reina, Juan David Cano-Moreno, and Manuel Gómez Gómez. 2024. "Analysis and Simulation of the Compressive Strength of Bioinspired Lightweight Structures Manufactured by a Stereolithography 3D Printer" Biomimetics 9, no. 4: 240. https://doi.org/10.3390/biomimetics9040240

APA Style

Alía García, C., Rodríguez Ortiz, Á., Arenas Reina, J. M., Cano-Moreno, J. D., & Gómez Gómez, M. (2024). Analysis and Simulation of the Compressive Strength of Bioinspired Lightweight Structures Manufactured by a Stereolithography 3D Printer. Biomimetics, 9(4), 240. https://doi.org/10.3390/biomimetics9040240

Article Metrics

Back to TopTop