Next Article in Journal
Applications of Smart Material Sensors and Soft Electronics in Healthcare Wearables for Better User Compliance
Next Article in Special Issue
Structural Design of a Special Machine Tool for Internal Cylindrical Ultrasonic-Assisted Electrochemical Grinding
Previous Article in Journal
Functional Surface Generation by EDM—A Review
Previous Article in Special Issue
Study on the Fatigue Strength of Welding Line in Injection Molding Products under Different Tensile Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Behavior of Honeycomb, Grid and Triangular PLA Structures under Symmetric and Asymmetric Bending

by
Vasile Cojocaru
,
Doina Frunzaverde
* and
Calin-Octavian Miclosina
Department of Engineering Science, Babeș-Bolyai University, Traian Vuia Square 1-4, 320085 Reșița, Romania
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(1), 120; https://doi.org/10.3390/mi14010120
Submission received: 31 October 2022 / Revised: 30 November 2022 / Accepted: 28 December 2022 / Published: 31 December 2022
(This article belongs to the Special Issue Manufacturing Methods or Processing Methods in Micromachines)

Abstract

:
Additive manufacturing technologies enable the production of components with lightweight cores, by means of infills with various patterns and densities. Together with reduced mass and material consumption, infill geometries must ensure that strength and stiffness conditions are fulfilled. For the proper correlation of the infill type with the loading case of the part, the mechanical behavior of the infill along all three principal axes of inertia has to be known. In this paper, the behavior in symmetric and asymmetric bending of three infill geometries, commonly used in 3D printing processes (honeycomb, grid and triangles) is analyzed. The variations of deflections as a function of force orientation are presented, showing that honeycomb and triangular structures exhibit similar behaviors along the Y and Z principal axes of inertia. Furthermore, the displacements obtained for the three types of structures are compared, in relation to the consumed volume of material. The larger displacements of the grid structure compared to the honeycomb and triangular structures are highlighted.

1. Introduction

Additive manufacturing technologies have developed intensively in recent decades. Specific to these technologies is the building of parts layer by layer, allowing the core of the parts to be generated with different infill geometries (alternative to full core parts). A decrease in mass and material consumption is achieved using an infill pattern. Furthermore, the infill patterns should be selected in relation to the loads applied to the part, aiming to optimize the stress distribution and static/dynamic behavior [1,2,3,4].
Infill geometries are defined by two main parameters: the infill pattern and the infill density [5]. The various slicer software used in FFF (fused filament fabrication) technology have predefined various infill patterns, the number of which has been increasing in recent years. Engineers can also build optimized infill geometries directly in computer-aided design (CAD) software according to the stress and strain states of the parts. Several categories of infill geometries can be defined:
-
Infill geometries defined by a 2D model, identical for two successive layers (the result is a 2.5D structure with prismatic cells oriented along the Z axis of the printer);
-
Infill geometries defined by 3D models (the infill pattern is defined by several successive layers with different configurations; this pattern is repeated after a specific number of layers; the cells can be open or closed) [6,7];
-
Lattice-type infill geometries [8];
-
Infill geometries with variable infill, changed from layer to layer [9] or changed between specified areas in the part [10,11,12];
-
Infill geometries with variable structure resulting from topological optimization.
Figure 1 shows examples of infill geometries defined by 2D models. The four examples (grid, triangles, concentric, cross) were generated with an infill density of 40%, for a prismatic volume 40 × 40 × 10 mm3, using the Ultimaker Cura 5.0.0 slicer software (Ultimaker B.V., Utrecht, Netherlands). For a better display, the infill geometries were generated without top/bottom layers and without peripheral contours (without shell).
Figure 2 shows four examples of infill geometries defined by 3D models, generated under the same conditions presented above.
The density of the infill indicates the percentage of the inner volume that is occupied by material (100% infill density indicates a solid volume; 0% infill density indicates a hollow part). Figure 3 shows a 40 × 40 × 10 mm3 prismatic volume with grid infill at four different densities: 20%, 40%, 60% and 80%. The infill was generated by the Ultimaker Cura 5.0.0 slicer software, using three perimeter contours and without upper and lower layers.
The peripheral contours, the lower layers and the upper layers constitute a shell that closes the filling volume and gives the shape of the piece. The ratio of the shell material volume to the infill material volume increases under certain conditions: small parts; low infill densities; high number of top/bottom layers and peripheral contours. In such cases, the contribution of the shell to the load bearing capacity of the piece increases as well [13]. An eloquent example in this regard is a tensile specimen with a small cross-section. For instance, the area occupied by the infill of an ASTM D638-14 type IV tensile specimen [14] with a 3.2 × 6 mm2 cross-section, printed in the ZX or ZY build orientation [15], with two peripheral layers of 0.4 mm (Figure 4), represents only 37% of the cross-section. To conclude, for an accurate analysis of the results obtained from mechanical tests on specimens with small cross-sections, the following parameters should also be considered: the number and the thickness of the lower and upper layers, the wall line number, and the wall thickness.
In the fused filament fabrication (FFF) technology, polymeric materials (polylactic acid—PLA, acrylonitrile butadiene styrene—ABS, polyethylene terephthalate—PET, etc.) and composite materials are used. For the components made by FFF printing, the mechanical properties are influenced by several process parameters, but also by pre-process conditions (filament state) and post-process conditions (storage conditions, post-process treatments, working environment) [16,17]. The orientation of the part on the printer build plate and the type of raster can determine the anisotropic behavior of the material (even if the material is isotropic in conventional manufacturing). Differences are observed mainly between the behavior on the two axes of the printer build plate (X and Y, according to ISO/ASTM 52921:2013) and the Z axis on which successive layers are built [16]. The pattern and the density of infill can increase the anisotropic behavior. Under these conditions, it is necessary to carry out the characterization of the mechanical behavior of components with different types of infill following each of the three principal axes of inertia.
Research carried out so far on polymeric materials printed with different infills has mainly focused on the analysis of tensile, compressive or bending behavior. A synthesis of the literature is presented below, including only papers in which specimens are printed from polylactic acid using the FFF technology.
Doboș et al. [18] analyze the tensile behavior for four types of infill (rectilinear, Hilbert curve, concentric and honeycomb) with 40–80% densities. The force–elongation curves show that at 80% density the highest breaking force is obtained for honeycomb infill, while at 40% and 60% densities the highest breaking force is obtained for concentric infill.
Guan [19] compares the tensile behavior of honeycomb, triangular and rectilinear infill specimens and shows that both UTS and Young’s modulus (E) are higher for honeycomb and triangular infill.
Gonzalez-Rebanque et al. [20] analyze the tensile behavior of PLA specimens made with rectilinear, triangle and honeycomb infill and show that at 80% infill density the highest value of UTS (ultimate tensile strength) is obtained for rectilinear infill specimens, and at 50% and 30% densities the highest UTS value is obtained for honeycomb infill.
In [21] the tensile behavior of six types of infill is analyzed: hexagonal, triangular, square, diamond, diamond angle and square angle. These infill structures were disposed only in the calibrated area of the ASTM D638 Type I specimen directly by CAD modeling. No top layers were defined, but peripheral contours were used. The stress–strain curves show that the highest tensile strengths were obtained for the square and triangle structures and the lowest tensile strength was obtained for the diamond structure.
Comparative analysis of tensile behavior for the ASTM D638-type I specimens with three peripheral contours [22] indicates that at 75% infill density, concentric-type structures have higher UTS compared to grid and tri-hexagonal structures.
The comparison of the tensile behavior of ISO 527-2 [23] specimens made with seven infill patterns, 60% infill density and two peripheral contours shows close values of the mechanical strength [24]. Analyzing the specimens, it can be concluded that the peripheral contours determine the mechanical behavior while the infill influence is low.
Research by Alafaghani et al. [25] indicates that the UTS obtained for ASTM D638 type IV specimens made with hexagonal infill and diamond infill is slightly higher than the UTS obtained for specimens with linear infill. The small differences obtained can be correlated with the small dimensions of the specimen cross-section (6 × 3.5 mm2) and the influence of the shell walls.
Tensile tests performed on ASTM D638 specimens [26] showed that the UTS is higher for grid infill compared to triangle and gyroid. This conclusion was obtained for all three infill densities analyzed (60%, 75%, 90%). The authors explain this behavior by the alignment of the grid infill with the axis of the tensile specimen.
Pernet and co-authors [27] test 14 types of infill in compression and plot the ratios of maximum compressive force to weight of specimens. The specimens used are cylindrical with a diameter of 12.7 mm and a length of 25.4 mm, and have three peripheral contours. The highest force/mass ratios are obtained for grid and cross infills.
In [28] the compression behavior of prismatic specimens (12.7 × 12.7 × 25.4 mm3) with six types of infill (octagram spiral, Hilbert curve, rectilinear, line, Archimedean curve and honeycomb) and four infill densities (ranging from 20% to 80%) are analyzed. It is shown that the Hilbert curve pattern specimens have the highest compressive strength (approximately double compared to the other infills).
Analyzing four types of infill (triangle, grid, quarter-cubic and tri-hexagon) in low-velocity impact tests and a compression test, Aloyaydi et al. [29] show that the grid type infill has the highest compressive strength and the triangle type infill has the highest absorbed energy in impact tests.
Birosz et al. [30] compare the bending behavior for three types of infill: grid, honeycomb and gyroid. The specimens with square cross-section of 10 × 10 mm2 were made with shells of constant thickness—0.8 mm (the section occupied by the infill is 8.4 × 8.4 mm2). The tests were performed with the force orientated along two axes, Y and Z. The results obtained by the authors highlight the higher breaking force of specimens with hexagonal and gyroid infill. It is noted that the breaking force and modulus of elasticity obtained along the X-axis are close to the breaking force and modulus of elasticity obtained along the Z-axis (this is evidenced for all three types of infill).
Honeycomb infill specimens (10×4 mm2 cross-section) show higher modulus of elasticity compared to linear and rectilinear infill specimens in four-point bending tests [31].
In [32] the bending properties of grid and honeycomb infills made at three different densities (10%, 25% and 50%) are analyzed. For the 50% density the honeycomb infill shows flexural strength higher by approx. 13% compared to the grid infill. The ISO 178:2011 [33] type specimens (10 × 4 mm2 section) are made with three peripheral contours. At low infill densities, the influence of the shell on the bending behavior can be higher than the influence of the infill.
In [34] the tensile and bending behavior of ASTM D638-14 [14] and ASTM D790-10 [35] specimens with rectilinear, concentric, honeycomb and Hilbert-curve infill are analyzed. No details are given on printing parameters, number of peripheral contours or number of top layers. In both bending and tensile tests, the highest mechanical strengths were obtained for the rectilinear and honeycomb infills.
The analysis of the impact behavior (Izod test) for 12 infill patterns at different infill densities reveals that the linear, zig-zag and concentric infills exhibit the highest energy absorbed on impact, at both 50% and 75% density [36].
Bonada et al. [37] analyze the type of infill correlated with the build orientation of the specimens and the type of the raster. It should be noted that the build orientation causes significant differences in infill orientation and specimen mass. Thus, at XY and YX orientations the infill is built on the thickness of the specimen, at YZ and XZ orientations the infill is built on the width of the specimen, and at ZY and ZX orientations the infill is built on the length of the specimen.
Most of the research reviewed was conducted on standardized specimens. The reduced thickness of these specimens leads to a low ratio of the infill material volume to the shell material volume. Under these conditions, the relevance of analyses of the mechanical behavior of different infill patterns is low. Furthermore, the analysis of the variation of the mechanical properties of the cellular structures in relation to the orientation of the force to the main axes of inertia is less studied in the literature.
The objective of the current research was to analyze the symmetric and asymmetric bending behavior for three types of infill structures widely used in additive manufacturing: honeycomb, grid and triangular. In particular, it aimed to determine the variation of the deflections as a function of the orientation of the force to the principal axes of inertia Y and Z. Furthermore, considering the volume of each type of structure, a comparative analysis of the behavior of the three structures was also carried out. The characterization of asymmetric bending behavior of infill structures is of major importance for mechanical components subjected to loads with variable direction (e.g., turbine blades, propellers).

2. Materials and Methods

The geometries of the three types of structures analyzed in this paper (honeycomb, grid and triangles) were modelled in the SolidWorks 2020-2021 software (Dassault Systèmes, Vélizy-Villacoublay, France). The structures were fitted into a prismatic volume of 120 mm length and 20 × 20 mm2 cross-section (Figure 5). A square cross-section was chosen so that a comparative analysis of Z- and Y-axis behavior could be performed. The cross-section dimensions of the specimen for the three-point bending tests were increased relative to standard dimensions in order to increase the number of cell rows arranged across the width of the specimen. The length of the specimen was correlated with the cross-section.
No peripheral contours or upper/lower layers were generated, as their dimensions introduce additional variables into the analysis. In the reference system attached to the three types of structures, the X-axis defines the length, and the Y- and Z-axes define the cross-section.
The structures defined by 2D patterns with regular polygonal cells (hexagon, square, triangle) are characterized by the specific dimensions of the polygons and the distances between two neighboring polygons (wall thickness) [38]. In this research the cell wall thicknesses were identical in all three structures, and the characteristic dimensions of the polygonal cells were chosen assuming that the specimen volumes were of the same scale order. A specific dimension of 3 mm was used to define the inner of the cells and the wall thickness was kept constant at 1 mm (Figure 6). The layout of the hexagonal/ grid/ triangular cells was made so that the structures were symmetrical relative to the median planes. Therefore, the X, Y and Z axes were the main axes of inertia.
To compare the mechanical behavior of the three structures, the differences in actual volume (reflecting material consumption) were calculated, based on the 3D models from SolidWorks (Table 1). The volume of the grid structure is 1.82% smaller compared to the hexagonal structure, while the volume of the triangular structure is 23.96% higher.
The 3D printing was done by fused filament fabrication (FFF), using an Ultimaker 2+ Connect printer (Ultimaker B.V., Utrecht, Netherlands). The geometries generated in SolidWorks were imported as stl file into the Ultimaker Cura 5.0.0 slicer software. The material used was Ultimaker Silver PLA (Ultimaker B.V., Utrecht, Netherlands) with a diameter of 2.85 mm. The main properties of the material (according to the manufacturer’s specifications) were: density—1.24 g/cm3, melting point—145–160 °C, thermal decomposition—250 °C, tensile modulus of elasticity—2890 MPa up to 3393 MPa, tensile stress at break—30.3 MPa up to 57.5 MPa, flexural modulus of elasticity—2693 MPa up to 3106 MPa and flexural strength—52 MPa up to 101.3 MPa. No pre-process treatments were applied to the filament. The print parameters (Table 2) were maintained constant for all three structures.
The selection of the printing parameters was done in accordance with the manufacturer’s recommendations and previous research [16,39]. The specimens were placed in the center of the build plate, following the XY orientation. No support material was used. Each specimen was printed individually.
Figure 7 shows the aspect of the three types of test specimens. No post-processing treatments were applied after 3D printing.
Two testing phases were completed to investigate the symmetric and asymmetric bending behavior of the structures:
(a)
The analysis of the three-point bending behavior (symmetric bending);
(b)
The analysis of the asymmetric bending behavior for cantilever beam loaded by a force disposed inclined to the main axes of inertia Y and Z.
The experimental three-point bending tests were carried out for two loading cases: force applied along the Y principal axis of inertia and force applied along the Z principal axis of inertia (Figure 8, Figure 9 and Figure 10). The tests were performed on a Mecmesin MultiTest 2.5 dV (PPT Group UK Ltd., Slinfold, UK) mechanical testing machine. The span between the 3-point bending device supporting pins was 100 mm and the force was applied in the middle section (Figure 8c). Displacements were measured for seven force values: 10 N, 20 N, 30 N, 40 N, 50 N, 60 N and 70 N. For the Z-direction force analysis, the specimens were positioned rotated (with the axis of the prismatic cells placed horizontally—Figure 8b, Figure 9b and Figure 10b). Three experimental tests were performed for each type of structure. The results are presented as force-deflection variations. Assuming that the force-deflection variations are linear (as the deflections are obtained in the elastic region of the material), the tangents of these variations were calculated.
Finite element analysis of the 3-point bending deflections was performed for loading layouts similar to the experimental tests (Figure 11 and Figure 12—example for the hexagonal structure; similarly loading layouts were used for the grid and the triangle structures). The force was varied from 10 N to 70 N, determining the force-deflection variations and the tangents of these variations. The set-up of the mesh was based on the convergence of the Y-axis maximum displacements. Subsequently, the mesh parameters were maintained identical for both bending cases of the same structure. Furthermore, the overall size of the finite elements was set to the same value (2 mm) for the three structures, resulting in: 52,444 finite elements with 95,468 nodes for the hexagonal structure, 50,663 finite elements with 95,739 nodes for the grid structure and 52,115 finite elements with 98,174 nodes for the triangle structure.
The deviations between the tangents of the force-deflection variations obtained experimentally and those obtained by finite element analysis were analyzed. After the validation of the FEM simulations for 3-point bending, the analysis models for asymmetric bending were settled. Asymmetric bending assumes that the resultant force acts along an axis inclined to the principal axes of inertia Y and Z.
In three-point bending the variation of the force direction implies the change of the fixing surfaces, which leads to additional variables introduced into the analysis. Therefore, the cantilever beam setup was used for the analysis of the displacement variation in asymmetric bending. The fixing conditions were applied on the left cross-section and the force was applied on the right cross-section of the beam for all loading cases (Figure 13—illustration for the hexagonal structure). An angle α was used to define the direction of the force in the (YZ) plane, relative to the Y-axis (Figure 13b). For α = 0° the force is applied along the Y-axis (Figure 13a) and for α = 90° the force is applied along the Z-axis (Figure 13c). The value of the applied force was kept constant at F = 5 N and the direction of the force was varied by 5° in the range 0°–90° (19 directions analyzed for each structure). In SolidWorks, the inclined force was applied as force projections (Figure 13b). For each type of structure, the variation of the displacement ( U r e s ) versus angle α was plotted.

3. Results and Discussion

This chapter presents the results obtained in the two stages of the analysis: three-point bending along the Y- and Z-axes and asymmetric bending of a cantilever beam.

3.1. Three-Point Bending

The three-point bending layout corresponds to a simply supported beam with force applied in the middle section. In this case, the maximum deflection U 3 PB   max of a beam with constant cross-section can be calculated with Equation (1):
U 3 PB   max = F l 3 E I
where l is the beam length (the length between the two supporting pins of the three-point bending device), E is Young’s modulus and I is the axial moment of inertia ( I z at Y-axis loading, respectively I y at Z-axis loading).
For a square cross-section the moments of inertia along the Y and Z axes are equal I z / I y = 1 , resulting in equal deflections on bending along the two axes U y = U z . For the structures with a cellular infill, the axial moments of inertia vary with the position of the cross-section along the X-axis. In addition, the axial moments of inertia I z and I y in any cross-section are not equal to each other and the ratio I z / I y depends on the position of the cross-section. The ratio of maximum deflection U y / U z will thus be dependent on the variation of the I z / I y ratios.
The experimental tests carried out at three-point bending aimed to determine the ratios of the maximum deflections U y / U z and to calculate the tangents of the force-deflection variations (these tangents reveal the elastic behavior of the structures).
The deflections obtained in the three tests carried out for each type of structure were represented graphically, plotting the average force-deflection variation. The tests were carried out in the elastic domain; thus, the variations were linear. Figure 14, Figure 15 and Figure 16 show the deflection variations U y = f ( F y ) (force applied on the Y-axis of the specimen) and U z = f ( F z ) (force applied on the Z-direction) resulting from the experimental tests for the three types of structures. For the hexagonal structure the variations U y = f ( F y ) and U z = f ( F z ) have similar slopes. The same type of behavior is observed for the triangle structure. In the grid structure the tangent of the variation U z = f ( F z ) is substantially higher than the tangent of the variation U y = f ( F y ) .
Figure 17 shows a comparison of the variations U z = f ( F z ) for the three types of structure. It can be seen that the triangular structure has the smallest deflections, followed by the hexagonal structure. The differences between the deflections of these two types of structure and the deflections of the hexagonal structure are large.
The comparative analysis must also include the volume of each structure. As shown in Table 1, hexagonal cell structures and grid cell structures have similar volumes V g r i d / V h e x = 0.98 . The ratio of the tangents of the deflection functions   θ z = tan ( U z = f ( F z ) ) is θ z g r i d / θ z h e x = 3.94 . Thereby, at the same material consumption, the hexagonal structure has smaller deflections compared to the grid structure.
For the comparative analysis between the hexagonal structure and the triangle structure the volume ratio is V t r i / V h e x = 1.23 , while the ratio of the tangents of the deflection functions is θ z t r i / θ z h e x = 0.70 . It follows that a 23% higher material consumption in the triangle structure leads to a 30% decrease of deflection compared to the hexagonal structure.
Figure 18 shows a comparison of the variations U y = f ( F y ) for the three types of structures. The tangent ratios of the deflection variations θ y = tan ( U y = f ( F y ) ) are: θ y g r i d / θ y h e x = 3.18 and θ y t r i / θ y h e x = 0.69 . The tendencies are similar with the U z = f ( F z ) variations.
FEM analyses of three-point bending were performed for seven force values, ranging from 10 N to 70 N, determining the functions U y = f ( F y ) and U z = f ( F z ) and the tangents of these functions. Figure 19, Figure 20 and Figure 21 show the charts of displacements for F = 70   N .
Table 3 summarizes the values of the tangents θ z and θ y obtained experimentally and by numerical simulation of three-point bending. The values of the coefficient of determination R2 are also presented. It is observed that for the variations determined by experimental measurements, the values of this coefficient are situated between 0.975 and 0.997, indicating small values of the deviations.
The conclusions of this analysis of the three-point bending displacements should be correlated with the orientation of the structures’ cells. The effect of changing the orientation of hexagonal/triangular/grid cells through rotation in the (XZ) plane was not included in this analysis. The cell orientation influences the consumption of the material [40] and the mechanical behavior of the structures [41]. Furthermore, changing the position of the force relative to the center of the polygonal cells determines the changes in the mechanical behavior of the structure [42].

3.2. The Variation of Resultant Displacement in Unsymmetrical Bending

For a cantilever beam loaded with asymmetric force F, oriented at angle α to the y-axis, the maximum resultant displacement can be calculated by the equation:
U r e s = ( ( F cos α l 3 E I z ) 2 + ( F sin α l 3 E I y ) 2 ) 1 / 2
where: l is the length of the bar; E is Young’s modulus; and I y and I z are the axial moments of inertia. From Equation (2) it is observed that if I z = I y (e.g., square cross-section), then U _ r e s is constant for any value of angle α.
For the three geometries considered here the cross-section is fitted into a 20 × 20 mm2 square, the Y and Z axes are the principal axes of inertia, but the axial moments of inertia I y and I z are not equal and are variable, depending on the position of the cross-section on the X axis.
To establish the variation of the displacements as a function of the orientation of the force (angle α ), the value of the resultant force was kept constant and the angle α was varied by 5°, in the range 0°–90°. Figure 22, Figure 23 and Figure 24 show the distribution of displacements for the three structures at angles α = 0° and α = 90°. The similar values of the maximum displacements for the hexagonal structure are observed, with the ratio U h e x ( α = 0 ) / U h e x ( α = 90 ) = 1.02 . For the triangular cell structure, the displacement ratio is U t r i ( α = 0 ) / U t r i ( α = 90 ) = 0.95 . The grid cell structure shows a significant variation in the resulting displacement U g r i d ( α = 0 ) / U g r i d ( α = 90 ) = 0.64 and larger displacements relative to the hexagonal and triangular structures.
Figure 25 shows the variation of the resultant displacements as a function of the angle α of inclination of the force. It can be seen that the extreme values are obtained for α = 0° and α = 90°. For the grid structure the largest displacement is obtained for the orientation of the force along the Z-axis (α = 90°), a conclusion also drawn from the experimental analysis at three-point bending.
The low variation of the resultant displacement in hexagonal and triangular geometries shows that these structures can be used in beam-like mechanical components loaded with force of variable orientation.
The conclusions resulting from the analysis of the displacements under asymmetric bending should be correlated with the cell orientations used for the three types of structure. Changing the orientation of the cells by rotation in the (XZ) plane can cause the change of the position of the principal axes of inertia.

4. Conclusions

The use of cellular structures for the core of mechanical components leads to a reduction in mass and material consumption and may conduce, at the same time, to an optimized stress distribution and an improved elastic behavior. To achieve these advantages, the mechanical behavior of these structures, as well as the variation of the mechanical properties in relation to the main axes of inertia, has to be known.
The number and thickness of top/bottom layers and the number of peripheral contours constitute variables that influence the analysis of infill types, especially when using low infill densities or test specimens with a small cross-section area. An alternative is the testing of structures without shells, as presented in this paper.
The analysis of symmetric and asymmetric bending displacements for the three structures considered in this paper (hexagonal structure, grid structure and triangular structure) led to the following conclusions:
-
The grid cell structure has significantly larger displacements than the hexagonal cell structure at the same material consumption.
-
At the same material consumption, the triangular cell structure has slightly smaller displacements than the hexagonal cell structure.
-
For hexagonal/triangular cell structures the displacements do not vary significantly with the orientation of the force relative to the principal axes of inertia. These types of structure are suitable for mechanical components where the force changes direction in the YZ plane.
-
For the grid structure the variation of displacements with the orientation of the force relative to the principal axes of inertia is large; displacements are higher when the force acts perpendicular to the axis of the prismatic cells of the grid geometry.
The analyzed structures show symmetry with respect to the median planes. The rotation of the cells in the XZ plane can cause the change of the principal axes of inertia and the modification of the bending behavior. These issues will be analyzed in future research.

Author Contributions

Conceptualization, V.C. and D.F.; methodology, V.C. and D.F.; software and validation, V.C. and C.-O.M.; investigation and data curation, V.C. and C.-O.M.; writing—original draft preparation, V.C.; writing—review and editing, V.C., D.F. and C.-O.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goh, G.D.; Yap, Y.L.; Tan, H.K.J.; Sing, S.L.; Goh, G.L.; Yeong, W.Y. Process–Structure–Properties in Polymer Additive Manufacturing via Material Extrusion: A Review. Crit. Rev. Solid State Mater. Sci. 2019, 45, 113–133. [Google Scholar] [CrossRef]
  2. Liu, J.; Gaynor, A.T.; Chen, S.; Kang, Z.; Suresh, K.; Takezawa, A.; Li, L.; Kato, J.; Tang, J.; Wang, C.C.L.; et al. Current and Future Trends in Topology Optimization for Additive Manufacturing. Struct. Multidiscip. Optim. 2018, 57, 2457–2483. [Google Scholar] [CrossRef] [Green Version]
  3. Bergonzi, L.; Vettori, M.; Stefanini, L.; D’Alcamo, L. Different Infill Geometry Influence on Mechanical Properties of FDM Produced PLA. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1038, 012071. [Google Scholar] [CrossRef]
  4. Kaidouchi, H.; Kebdani, S.; Slimane, S.A. Vibro-Acoustic Analysis of the Sound Transmission through Aerospace Composite Structures. Mech. Adv. Mater. Struct. 2022, 29, 1–11. [Google Scholar] [CrossRef]
  5. Qamar Tanveer, M.; Mishra, G.; Mishra, S.; Sharma, R. Effect of Infill Pattern and Infill Density on Mechanical Behaviour of FDM 3D Printed Parts—A Current Review. Mater. Today Proc. 2022, 62, 100–108. [Google Scholar] [CrossRef]
  6. Ma, R.; Li, M.; Xu, Y.; Meloni, M.; Feng, J.; Cai, J. Geometry Design and In-Plane Compression Performance of Novel Origami Honeycomb Material. Thin-Walled Struct. 2022, 181, 110111. [Google Scholar] [CrossRef]
  7. Zhang, K.; Xie, X.; Wang, C.; Wang, H.; Xu, F.; Wang, H.; Zhang, X.; Guan, H.; Qu, H.; Zhang, J. Optomechanical Performances of Advanced Lightweight Mirrors Based on Additive Manufacturing. Micromachines 2022, 13, 1334. [Google Scholar] [CrossRef]
  8. Qin, D.; Sang, L.; Zhang, Z.; Lai, S.; Zhao, Y. Compression Performance and Deformation Behavior of 3D-Printed PLA-Based Lattice Structures. Polymers 2022, 14, 1062. [Google Scholar] [CrossRef]
  9. Kim, S.; Andreu, A.; Kim, I.; Kim, J.-H.; Lee, J.; Yoon, Y.-J. Continuously Varied Infill Pattern (ConVIP): Improvement of Mechanical Properties and Printing Speed of Fused Filament Fabrication (FFF) 3D Printing. J. Mater. Res. Technol. 2022, 18, 1055–1069. [Google Scholar] [CrossRef]
  10. Paczos, P.; Wichniarek, R.; Magnucki, K. Three-Point Bending of Sandwich Beam with Special Structure of the Core. Compos. Struct. 2018, 201, 676–682. [Google Scholar] [CrossRef]
  11. Qi, C.; Jiang, F.; Yang, S. Advanced Honeycomb Designs for Improving Mechanical Properties: A Review. Compos. Part B Eng. 2021, 227, 109393. [Google Scholar] [CrossRef]
  12. Maszybrocka, J.; Dworak, M.; Nowakowska, G.; Osak, P.; Łosiewicz, B. The Influence of the Gradient Infill of PLA Samples Produced with the FDM Technique on Their Mechanical Properties. Materials 2022, 15, 1304. [Google Scholar] [CrossRef] [PubMed]
  13. Lubombo, C.; Huneault, M.A. Effect of Infill Patterns on the Mechanical Performance of Lightweight 3D-Printed Cellular PLA Parts. Mater. Today Commun. 2018, 17, 214–228. [Google Scholar] [CrossRef]
  14. ASTM D638-14; Standard Test Method for Tensile Properties of Plastics. ASTM International: West Conshohocken, PA, USA, 2014.
  15. ISO/ASTM 52921:2013; Standard Terminology for Additive Manufacturing—Coordinate Systems and Test Methodologies. International Organization for Standardization: Geneva, Switzerland, 2013.
  16. Cojocaru, V.; Frunzaverde, D.; Miclosina, C.-O.; Marginean, G. The Influence of the Process Parameters on the Mechanical Properties of PLA Specimens Produced by Fused Filament Fabrication—A Review. Polymers 2022, 14, 886. [Google Scholar] [CrossRef] [PubMed]
  17. Dey, A.; Yodo, N. A Systematic Survey of FDM Process Parameter Optimization and Their Influence on Part Characteristics. J. Manuf. Mater. Process. 2019, 3, 64. [Google Scholar] [CrossRef] [Green Version]
  18. Dobos, J.; Hanon, M.; Oldal, I. Effect of infill density and pattern on the specific load capacity of FDM 3D-printed PLA multi-layer sandwich. J. Polym. Eng. 2022, 42, 118–128. [Google Scholar] [CrossRef]
  19. Guan, R.; Smith, D. Influence of Infill Parameters on the Tensile Mechanical Properties of 3D Printed Parts. J. Emerg. Investig. 2020, 2, 1–5. Available online: https://emerginginvestgators.org (accessed on 15 September 2022).
  20. Rebenaque, A.G.; González-Requena, I. Study of bending test of specimens obtained through FDM processes of additive manufacturing. Procedia Manuf. 2019, 41, 859–866. [Google Scholar] [CrossRef]
  21. John, J.; Devjani, D.; Ali, S.; Abdallah, S.; Pervaiz, S. Optimization of 3D Printed Polylactic Acid Structures with Different Infill Patterns Using Taguchi-Grey Relational Analysis. Adv. Ind. Eng. Polym. Res. 2022, in press. [Google Scholar] [CrossRef]
  22. Rismalia, M.; Hidajat, S.C.; Permana, I.G.R.; Hadisujoto, B.; Muslimin, M.; Triawan, F. Infill Pattern and Density Effects on the Tensile Properties of 3D Printed PLA Material. J. Phys. Conf. Ser. 2019, 1402, 044041. [Google Scholar] [CrossRef] [Green Version]
  23. ISO 527-2:2012; Plastics—Determination of Tensile Properties—Part 2: Test Conditions for Molding and Extrusion Plastics. International Organization for Standardization: Geneva, Switzerland, 2012.
  24. Pandzic, A.; Hodzic, D. Mechanical Properties Comparison of PLA, Tough PLA and PC 3D Printed Materials with Infill Structure—Influence of Infill Pattern on Tensile Mechanical Properties. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1208, 012019. [Google Scholar] [CrossRef]
  25. Alafaghani, A.; Qattawi, A.; Alrawi, B.; Guzman, A. Experimental optimization of fused deposition modelling processing parameters: A design-for-manufacturing approach. Procedia Manuf. 2017, 10, 791–803. [Google Scholar] [CrossRef]
  26. Ambati, S.S.; Ambatipudi, R. Effect of infill density and infill pattern on the mechanical properties of 3D printed PLA parts. Mater. Today Proc. 2022, 64, 804–807. [Google Scholar] [CrossRef]
  27. Pernet, B.; Nagel, J.K.; Zhang, H. Compressive Strength Assessment of 3D Printing Infill Patterns. Procedia CIRP 2022, 105, 682–687. [Google Scholar] [CrossRef]
  28. Yadav, P.; Sahai, A.; Sharma, R.S. Strength and Surface Characteristics of FDM-Based 3D Printed PLA Parts for Multiple Infill Design Patterns. J. Inst. Eng. India Ser. C 2020, 102, 197–207. [Google Scholar] [CrossRef]
  29. Aloyaydi, B.; Sivasankaran, S.; Mustafa, A. Investigation of infill-patterns on mechanical response of 3D printed poly-lactic-acid. Polym. Test. 2020, 87, 106557. [Google Scholar] [CrossRef]
  30. Birosz, M.T.; Ledenyak, D.; Ando, M. Effect of FDM infill patterns on mechanical properties. Polym. Test. 2022, 113, 107654. [Google Scholar] [CrossRef]
  31. Travieso-Rodriguez, J.A.; Jerez-Mesa, R.; Llumà, J.; Traver-Ramos, O.; Gomez-Gras, G.; Roa Rovira, J.J. Mechanical Properties of 3D-Printing Polylactic Acid Parts subjected to Bending Stress and Fatigue Testing. Materials 2019, 12, 3859. [Google Scholar] [CrossRef] [Green Version]
  32. Kołodziej, A.; Żur, P.; Baier, A.; Borek, W. Influence of 3D-Printing Parameters on Mechanical Properties of PLA Defined in the Static Bending Test. Eur. J. Eng. Sci. Technol. 2019, 2, 65–70. [Google Scholar] [CrossRef] [Green Version]
  33. ISO 178:2019; Plastics—Determination of Flexural Properties. International Organization for Standardization: Geneva, Switzerland, 2019.
  34. Khan, S.F.; Zakaria, H.; Chong, Y.L.; Saad, M.A.M.; Basaruddin, K. Effect of Infill on Tensile and Flexural Strength of 3D Printed PLA Parts. IOP Conf. Ser. Mater. Sci Eng. 2018, 429, 012101. [Google Scholar] [CrossRef]
  35. ASTM D790-10; Standard Test Methods for Flexural Properties of Unreinforced and Rein-forced Plastics and Electrical Insulating Materials. ASTM International: West Conshohocken, PA, USA, 2010.
  36. Mishra, P.K.; Senthil, P.; Adarsh, S.; Anoop, M.S. An Investigation to Study the Combined Effect of Different Infill Pattern and Infill Density on the Impact Strength of 3D Printed Polylactic Acid Parts. Compos. Commun. 2021, 24, 100605. [Google Scholar] [CrossRef]
  37. Bonada, J.; Pastor, M.M.; Buj-Corral, I. Influence of Infill Pattern on the Elastic Mechanical Properties of Fused Filament Fabrication (FFF) Parts through Experimental Tests and Numerical Analyses. Materials 2021, 14, 5459. [Google Scholar] [CrossRef] [PubMed]
  38. Harpool, T.D.; Alarifi, I.M.; Alshammari, B.A.; Aabid, A.; Baig, M.; Malik, R.A.; Mohamed Sayed, A.; Asmatulu, R.; EL-Bagory, T.M.A.A. Evaluation of the Infill Design on the Tensile Response of 3D Printed Polylactic Acid Polymer. Materials 2021, 14, 2195. [Google Scholar] [CrossRef] [PubMed]
  39. Frunzaverde, D.; Cojocaru, V.; Ciubotariu, C.-R.; Miclosina, C.-O.; Ardeljan, D.D.; Ignat, E.F.; Marginean, G. The Influence of the Printing Temperature and the Filament Color on the Dimensional Accuracy, Tensile Strength, and Friction Performance of FFF-Printed PLA Specimens. Polymers 2022, 14, 1978. [Google Scholar] [CrossRef] [PubMed]
  40. Suteja, J. Effect of Infill Pattern, Infill Density, and Infill Angle on the Printing Time and Filament Length of 3D Printing. J. Rekayasa Mesin 2021, 12, 145. [Google Scholar] [CrossRef]
  41. Rajpal, R.; Lijesh, K.P.; Gangadharan, K.V. Parametric Studies on Bending Stiffness and Damping Ratio of Sandwich Structures. Addit. Manuf. 2018, 22, 583–591. [Google Scholar] [CrossRef]
  42. Slimane, S.; Kebdani, S.; Boudjemai, A.; Slimane, A. Effect of Position of Tension-Loaded Inserts on Honeycomb Panels Used for Space Applications. Int. J. Interact. Des. Manuf. 2017, 12, 393–408. [Google Scholar] [CrossRef]
Figure 1. The top views of a prismatic volume (40 × 40 × 10 mm3) with various infill patterns defined by 2D models: (a) grid, (b) triangles, (c) concentric, (d) cross.
Figure 1. The top views of a prismatic volume (40 × 40 × 10 mm3) with various infill patterns defined by 2D models: (a) grid, (b) triangles, (c) concentric, (d) cross.
Micromachines 14 00120 g001
Figure 2. The perspective views of a prismatic volume (40 × 40 × 10 mm3) with various infill patterns defined by 3D models: (a) cross 3D, (b) gyroid, (c) quarter cubic, (d) octet subdivision.
Figure 2. The perspective views of a prismatic volume (40 × 40 × 10 mm3) with various infill patterns defined by 3D models: (a) cross 3D, (b) gyroid, (c) quarter cubic, (d) octet subdivision.
Micromachines 14 00120 g002
Figure 3. The top views of a prismatic volume (40 × 40 × 10 mm3) with grid infill at various densities: (a) 20%, (b) 40%, (c) 60%, (d) 80%.
Figure 3. The top views of a prismatic volume (40 × 40 × 10 mm3) with grid infill at various densities: (a) 20%, (b) 40%, (c) 60%, (d) 80%.
Micromachines 14 00120 g003
Figure 4. Cross-section of an ASTM D638-14 type IV specimen with 90% grid infill, two wall lines and ZY build orientation.
Figure 4. Cross-section of an ASTM D638-14 type IV specimen with 90% grid infill, two wall lines and ZY build orientation.
Micromachines 14 00120 g004
Figure 5. The specimens’ overall dimensions—illustration for the hexagonal structure.
Figure 5. The specimens’ overall dimensions—illustration for the hexagonal structure.
Micromachines 14 00120 g005
Figure 6. The detailed views of structures’ cells: (a) hexagonal, (b) grid, (c) triangles.
Figure 6. The detailed views of structures’ cells: (a) hexagonal, (b) grid, (c) triangles.
Micromachines 14 00120 g006
Figure 7. The test specimens: (a) hexagonal structure, (b) grid structure, (c) triangle structure.
Figure 7. The test specimens: (a) hexagonal structure, (b) grid structure, (c) triangle structure.
Micromachines 14 00120 g007
Figure 8. The 3-point bending setup for the hexagonal cells specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen, (c) device layout.
Figure 8. The 3-point bending setup for the hexagonal cells specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen, (c) device layout.
Micromachines 14 00120 g008aMicromachines 14 00120 g008b
Figure 9. The 3-point bending setup for the grid cells specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen.
Figure 9. The 3-point bending setup for the grid cells specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen.
Micromachines 14 00120 g009
Figure 10. The 3-point bending setup for the triangle cell specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen.
Figure 10. The 3-point bending setup for the triangle cell specimen: (a) force along the Y-axis of the specimen, (b) force along the Z-axis of the specimen.
Micromachines 14 00120 g010
Figure 11. The FEM analysis setup for 3-point bending of the hexagonal structure—force applied along the Y-axis of the specimen: (a) general view, (b) detailed view of the loading force, (c) detailed view on fixing conditions (identically for both fixing areas).
Figure 11. The FEM analysis setup for 3-point bending of the hexagonal structure—force applied along the Y-axis of the specimen: (a) general view, (b) detailed view of the loading force, (c) detailed view on fixing conditions (identically for both fixing areas).
Micromachines 14 00120 g011
Figure 12. The FEM analysis setup for 3-point bending of the hexagonal structure—force applied along the Z-axis of the specimen: (a) general view, (b) detailed view of the loading force, (c) detailed view on fixing conditions (identically for both fixing areas).
Figure 12. The FEM analysis setup for 3-point bending of the hexagonal structure—force applied along the Z-axis of the specimen: (a) general view, (b) detailed view of the loading force, (c) detailed view on fixing conditions (identically for both fixing areas).
Micromachines 14 00120 g012
Figure 13. The FEM analysis setup for asymmetric bending of the hexagonal structure: (a) α = 0°, (b) α = 45°, (c) α = 90°.
Figure 13. The FEM analysis setup for asymmetric bending of the hexagonal structure: (a) α = 0°, (b) α = 45°, (c) α = 90°.
Micromachines 14 00120 g013
Figure 14. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the hexagonal structure, in experimental three-point bending.
Figure 14. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the hexagonal structure, in experimental three-point bending.
Micromachines 14 00120 g014
Figure 15. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the grid structure, in experimental three-point bending.
Figure 15. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the grid structure, in experimental three-point bending.
Micromachines 14 00120 g015
Figure 16. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the triangle structure, in experimental three-point bending.
Figure 16. The variations of deflections Uy = f(Fy) and Uz = f(Fz) for the triangle structure, in experimental three-point bending.
Micromachines 14 00120 g016
Figure 17. The comparative analysis of the variation of deflection U z = f ( F z ) for the hexagonal, triangle and grid structures, in experimental three-point bending.
Figure 17. The comparative analysis of the variation of deflection U z = f ( F z ) for the hexagonal, triangle and grid structures, in experimental three-point bending.
Micromachines 14 00120 g017
Figure 18. The comparative analysis of the variation of deflection U y = f ( F y ) for the hexagonal, triangle and grid structures, in experimental three-point bending.
Figure 18. The comparative analysis of the variation of deflection U y = f ( F y ) for the hexagonal, triangle and grid structures, in experimental three-point bending.
Micromachines 14 00120 g018
Figure 19. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the hexagonal structure: (a) Fy = 70 N, (b) Fz = 70 N.
Figure 19. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the hexagonal structure: (a) Fy = 70 N, (b) Fz = 70 N.
Micromachines 14 00120 g019
Figure 20. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the grid structure: (a) Fy = 70 N, (b) Fz = 70 N.
Figure 20. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the grid structure: (a) Fy = 70 N, (b) Fz = 70 N.
Micromachines 14 00120 g020
Figure 21. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the triangle structure: (a) Fy = 70 N, (b) Fz = 70 N.
Figure 21. The Uy and Uz displacements charts obtained from FEM analysis of three-point bending of the triangle structure: (a) Fy = 70 N, (b) Fz = 70 N.
Micromachines 14 00120 g021
Figure 22. Cantilever beam. The displacement charts for the hexagonal structure: (a) α = 0°, (b) α = 90°.
Figure 22. Cantilever beam. The displacement charts for the hexagonal structure: (a) α = 0°, (b) α = 90°.
Micromachines 14 00120 g022
Figure 23. Cantilever beam. The displacement charts for the grid structure: (a) α = 0°, (b) α = 90°.
Figure 23. Cantilever beam. The displacement charts for the grid structure: (a) α = 0°, (b) α = 90°.
Micromachines 14 00120 g023
Figure 24. Cantilever beam. The displacement charts for the triangle structure: (a) α = 0°, (b) α = 90°.
Figure 24. Cantilever beam. The displacement charts for the triangle structure: (a) α = 0°, (b) α = 90°.
Micromachines 14 00120 g024
Figure 25. The influence of force orientation on the resultant displacement at asymmetric bending, U r e s = f ( α ) .
Figure 25. The influence of force orientation on the resultant displacement at asymmetric bending, U r e s = f ( α ) .
Micromachines 14 00120 g025
Table 1. The comparative analysis of volume for the three specimens.
Table 1. The comparative analysis of volume for the three specimens.
SampleVolume
[mm3]
Volume Variation rel.
to Hexagonal
[%]
Hexagonal21,353.070.00
Grid20,964.76−1.82
Triangular26,468.8723.96
Table 2. The 3D printing parameters.
Table 2. The 3D printing parameters.
ParametersValues
Processfused filament fabrication (FFF)
Layer thickness, t 0.1 mm
Printing head temperature, TH210 °C
Build plate temperature, TB60 °C
Printing speed, sp50 mm/s
Nozzle diameter, dn0.40 mm
Filament diameter, df2.85 mm
Build orientation (acc. to [15])XY
MaterialUltimaker PLA Silver
Table 3. The comparative analysis the tangents of the U y = f ( F y ) and U z = f ( F z ) functions obtained by experiment and the FEM simulation of three-point bending.
Table 3. The comparative analysis the tangents of the U y = f ( F y ) and U z = f ( F z ) functions obtained by experiment and the FEM simulation of three-point bending.
GeometryUy = f(Fy)Uz = f(Fz)
ExperimentalFEMExperimentalFEM
θ y R2 θ y R2 θ z R2 θ z R2
Hexagonal0.00950.9820.009510.01140.9890.01131
Grid0.03030.9970.018310.04500.9970.04781
Triangular0.00660.9750.005810.00800.9870.00881
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

Cojocaru, V.; Frunzaverde, D.; Miclosina, C.-O. On the Behavior of Honeycomb, Grid and Triangular PLA Structures under Symmetric and Asymmetric Bending. Micromachines 2023, 14, 120. https://doi.org/10.3390/mi14010120

AMA Style

Cojocaru V, Frunzaverde D, Miclosina C-O. On the Behavior of Honeycomb, Grid and Triangular PLA Structures under Symmetric and Asymmetric Bending. Micromachines. 2023; 14(1):120. https://doi.org/10.3390/mi14010120

Chicago/Turabian Style

Cojocaru, Vasile, Doina Frunzaverde, and Calin-Octavian Miclosina. 2023. "On the Behavior of Honeycomb, Grid and Triangular PLA Structures under Symmetric and Asymmetric Bending" Micromachines 14, no. 1: 120. https://doi.org/10.3390/mi14010120

APA Style

Cojocaru, V., Frunzaverde, D., & Miclosina, C. -O. (2023). On the Behavior of Honeycomb, Grid and Triangular PLA Structures under Symmetric and Asymmetric Bending. Micromachines, 14(1), 120. https://doi.org/10.3390/mi14010120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop