Numerical Modeling and Experimental Investigation of Effective Numerical Modeling and Experimental Investigation of Effective Elastic Properties of the 3D Printed Gyroid Infill Elastic Properties of the 3D Printed Gyroid Infill

: A numerical homogenization approach is presented for the effective elastic moduli of 3D printed cellular inﬁlls. A representative volume element of the inﬁll geometry is discretized using either shell or solid elements and analyzed using the ﬁnite element method. The elastic moduli of the bulk cellular material are obtained through longitudinal and shear deformations of a representative volume element under periodic boundary conditions. The method is used to analyze the elastic behavior of gyroid inﬁlls for varying inﬁll densities. The approach is validated by comparing the Young’s modulus and Poisson’s ratio with those obtained from compression experiments. Results indicate that although the gyroid inﬁll exhibits cubic symmetry, it is nearly isotropic with a low anisotropy index. The numerical predictions are used to develop semi-empirical equations of the effective elastic moduli of gyroid inﬁlls as a function of inﬁll density in order to inform design and topology optimization workﬂows.


Introduction
Additive manufacturing (AM) is a method whereby parts are manufactured by progressively depositing material layer by layer to achieve a desired shape. This differs from traditional (subtractive) manufacturing, where material is progressively removed from a parent block until it reaches the desired shape. Additive manufacturing can be performed in a variety of ways (including material-extrusion, powder-bed-fusion, binder-jetting, and stereolithography [1,2]), each of which has unique benefits and drawbacks. Additive manufacturing in general, though, has a number of advantages over traditional (subtractive) manufacturing methods [3][4][5][6][7][8]. One of these is the ability to manufacture lightweight parts which have an external skin and sparse internal infill. This decreases the weight of a part, which also decreases the amount of material required to manufacture the part. This, in turn, decreases the cost of the part in terms of machine-time, operator-time, and energyconsumption [9]. In order to design a part which takes advantage of this, however, it is first necessary to characterize the structural properties of the infill.
There are a variety of infills commonly implemented in AM softwares, all of which have both benefits and drawbacks. The most common class of infills are the 2D infills (including the particularly common rectilinear infill) which are generated by printing the same pattern for every layer. This results in a structure that has a constant cross-section vertically (as if extruded in this direction). These infills have the advantage of simplicity, which makes them quick and reliable to print, but the infill geometry results in highly anisotropic structural properties. Bhandari and Lopez-Anido [10][11][12], have developed methods of homogenizing the behavior of sparse rectilinear infills into orthotropic effective properties, which could be very useful for designing components using this infill.
The next class of common infills is the 3D infills, such as cubic infill, which consist of a structure that spatially repeats in all three directions. Many of these exhibit cubic symmetry, which more closely approximates isotropic behavior. This benefit, however, comes at the cost of added complexity and more difficulty in printing. One infill that is particularly promising for future design work is the gyroid, a 3D infill based on a mathematical construct known as a triply periodic minimal surface. This infill is nearly isotropic, prints relatively quickly, and is quite strong [13][14][15]. All of these properties are very useful for simplifying design, but very little specific elastic property data is available in the literature.
The final category of infills is the least well defined. It consists of those infills which do not fall into either of the previous categories due to their non-periodic nature, which often arises from optimization processes. These can take several forms, including using graded infills [16], internal structures containing voids [17][18][19], or internal structures with repeating optimized lattice structures which can be tuned to have specific elastic properties or graded to have spatially varying properties [19][20][21][22][23][24][25][26][27]. These repeating structures often take the form of standard rectangular, rhombic or hexagonal infill patterns whose properties have been studied in [28][29][30][31]. Research has also been performed on strut-based lattices with spatiallyvarying strut-thicknesses [32][33][34][35]. Several researchers have used topology optimization in order to directly tailor the microstructures to exhibit specific properties, thus allowing tuning of elastic moduli [24], Poisson's ratio [20] , and general orthotropic properties [36]. Martinez et al. [36], developed a systematic analysis of stochastic foam-like infill materials, which can be tuned to nearly any set of orthotropic material properties. Most of these researchers have generated unique infill forms which are not readily implemented in existing slicing software. For the purposes of this research, we focus on the common infill patterns that exist in current slicing software, so as to make this work relevant without a purpose-built slicer.
One thing that is missing in order to utilize shell-infill structures in everyday part design is data on how these infills behave as their relative density varies. "Relative density" or "infill density" of a 3D printed infill represents the fraction of space occupied by material, ranging from 0 (no material) to 1 (solid material). This value is determined as the density of infill divided by density of neat material, and is often converted to percentage for readability. While this data can be acquired experimentally, the process will need to be conducted again for every new infill type before implementing it in future designs. At the opposite end of the design spectrum, internal structural optimization could be done for each individual part, thus generating a specialized non-periodic infill from scratch for every component designed. Presented in this research is the intermediate approach, where finite element analysis (FEA) is used to analyze the behavior of a common, periodic infill at the level of the unit-cell, thus developing a method for future analysis of other infills. To this end, a method has been developed to predict the effective elastic properties of the gyroid infill at various relative densities. Experiments are conducted to verify the FEA predictions, and a set of empirical equations have been developed to make these predictions convenient for use in making design decisions. This analysis methodology is straightforward, and could be readily adapted in order to make similar behavioral predictions of other infill types.

Modeling and Finite Element Analysis
In this section, the gyroid infill geometry is introduced. This geometry is used to develop methods of model mesh generation, which are then used in conjunction with periodic boundary conditions to analyze the gyroid's structural behavior.

Gyroid Infill Geometry
The gyroid surface ( see Figure 1) is one member in the mathematical family known as "triply-periodic minimal surfaces", and was first discovered by Schoen [37]   Upon its recent addition as an infill option in some slicers, the gyroid has produced much interest among hobbyist printers due to its beneficial properties such as near isotropy, quick print time, and strength [13][14][15]. For this research, the gyroid has been selected as the primary focus, and will be analyzed using the finite element method in order to develop predictive models of its behavior as the relative density varies, while the focus of this work is on the gyroid, the same method will be useful as additional 3D infills are implemented. In particular, infills with complex-shapes (such as other triply-periodic minimal surfaces) may warrant detailed investigation.
It is important to note that there are two types of gyroid infill. The first, which is the focus of this work, is that defined by thickening the gyroid minimal surface. The second is a separate but related infill that has been investigated by Zhang et al. [38] and others. Rather than a thickened minimal surface, this second type is a lattice of intersecting curved struts defined by filling in the empty space that the gyroid surface encloses. This "gyroid-strut" model has great potential regarding continuously variable relative density [38], however, it has not yet been implemented in common slicers.
In order to determine the elastic response of the gyroid infill, finite element analysis (FEA) was implemented to model a unit cell, similar to that done by Bhandari and Lopez-Anido [10][11][12]. Unlike the rectangular infill, however, the gyroid is a complexly curved surface and cannot be reasonably simplified to a space-frame model. The logical first step in analyzing this shape, then, was to use a shell model.

Shell Element Discretization
In a shell model, the midsurface of the part is represented by a series of thin elements. The formulation is based on the plane stress assumption wherein the only significant stress components are the in-plane normal and shear stresses due to bending moments and in-plane forces. The out-of-plane normal and shear stresses are typically small for a thin shell element and can be neglected. For the purpose of modeling the gyroid, it was assumed that, for low relative densities, these through-thickness effects are small enough as to be insignificant.
The first step in modeling the gyroid is to generate a mesh of the surface using triangular shell elements. Since tools for generating such a mesh based on an implicit equation were not readily available, this needed to be done using a purpose-built script, which built up the mesh using implicit function plotting tools to locate nodes along the surface. These nodes were then triangulated into a mesh. Since the implicit plotting tools are optimized for properly displaying the geometry rather than for numerical analysis, the resulting mesh contained a number of elements whose shapes were unacceptable for FEA modeling. To remedy this, a series of mesh cleanup steps were performed using mesh-smoothing, and remeshing tools included in an FEA software. These cleanup tools rearrange nodes within the surface with minimal out-of-plane deviation in order to produce a mesh which matches the geometry, but whose elements are closer to the ideal equilateral triangles for FEA. The final mesh can be seen in Figure 2. The meshed gyroid was then analyzed using FEA with periodic boundary conditions (discussed in Section 2.4), and assuming a linear elastic response. Due to the nature of the periodic boundary conditions, each node on a boundary of the RVE must correspond exactly with one node on the opposite boundary. This adjustment is made by re-seeding the boundaries with nodes that correspond with each other, and re-aligning the mesh to these new bounding-nodes. In order to simulate different relative densities, several different shell thicknesses were modeled, and a curve of stiffness vs. relative density was generated. These showed clear deviation from expected behavior, particularly at high relative densities where the gyroid is expected to approach neat-material behavior. At these high relative densities, the structure is no longer thin-walled, which means that the shell assumptions are invalid and through-thickness stress effects are not negligible. Additionally, in some cases the thickness in the shell creates material overlap, thus double counting certain regions.

Solid Element Discretization
In order to ameliorate the shell issues and accurately predict the high relative density gyroid's behavior, a solid-model was developed.
To generate this model, the shell midsurface was converted to a solid model by first taking a 3 × 3 × 3 unit-cell shell mesh as an .stl file, and thickening it by offsetting it along the normal. This geometry is then trimmed down to a single unit-cell to remove the edge distortion caused by the offsetting tool. The result was an .stl of the outer surface of the thickened gyroid. This mesh is cleaned up using the same remeshing tools as the shell model, and the boundary nodes adjusted for compatibility with the periodic boundary constraints. Here, the boundary adjustment was accomplished by completely replacing a boundary face with the mesh from the opposite face. Finally, the thickened surface mesh was converted to a solid using an automatic tetramesher which builds a mesh of solid tetrahedral elements to fill the space enclosed by the triangular surface mesh. The resulting solid model (see Figure 3) was analyzed using the same process as the shell model (linear response under periodic boundary conditions). In addition, a mesh convergence study was performed, the results of which can be found in Section 4.

Periodic Boundary Conditions and Elastic Property Calculations
Periodic boundary conditions allow for the modeling of a bulk material using a representative volume element (RVE). In order to do this, the nodes on the boundary of the RVE are constrained so that each boundary node's displacements are matched to those of the corresponding node on the opposite boundary. An RVE constrained in this way has been shown to deform as if it were a single member in a bulk material composed of identical RVEs repeating in all directions [39]. The requirement that nodes be matched with corresponding nodes on opposite faces adds one additional step to the mesh generation, wherein the boundaries are re-seeded with nodes that correspond with each other, and the mesh realigned to these new bounding-nodes.
Once the boundary constraints are in place, a displacement is applied to one face, with the remaining faces free of traction. The resulting forces and displacements of the nodes are then calculated using FEA. Sun and Vaidya [39] demonstrated using Gauss' divergence theorem that the average strain through the RVE can be determined directly from the displacement of the surface (2), where V is RVE volume, S is a boundary surface, u i is displacement in the i-direction, n i is the i-component of the surface normal, and¯ ij is the averaged i, j component of the strain tensor. (3) shows how this reduces to the conventional strain definition for longitudinal strains in the x-direction,¯ where¯ xx is the average normal component of strain in the x-direction, L x and ∆L x are the length and elongation in the x-direction, respectively. Likewise, using the principle of strain energy and external work equivalence, the average stress is shown to be by the applied force divided by the area of the surface [39] (4), whereσ xx is the average normal component of stress in the x-direction, F x is the normal force resultant summed across the boundary (x = L x ) of the RVE, and A x is the area of this face. While the local stresses and strains vary throughout the model, these "average" properties are the behavior observed in a bulk material composed of these RVEs. The applied strain and predicted stress can then be used to calculate the effective elastic properties of this bulk-structure. Several similar analyses are needed to determine all of the properties of the bulk structure in different directions. Normal strains applied in any of the cardinal directions allow for calculation of the corresponding longitudinal modulus and the two corresponding Poisson's ratios. For example, applying a normal displacement ∆L x to the boundary (x = L x ) of the RVE causes it to deform such that there is also a length change in the y, and z-directions (∆L y and ∆L z ), and nodal resultant forces along the face x = L x which sum to F x . Dividing F x by the area of the face yields a resultant stress. These deformations/resultants are used in the standard definition of elastic modulus and Poisson's ratio (5)-(7).
Likewise, to calculate shear modulus G xy = G yx , a pure-shear state is imposed on the RVE by displacing the face x = L x in the y-direction by an amount ∆B xy , and the face y = L y face in the x-direction by an amount ∆B yx . This generates resultant shear forces on these faces (F xy and F yx , respectively). As before, these deformations and resultants are used in the standard formula for shear modulus (8). Figure 4 demonstrates the deformation in these two cases. For convenience, the constraints were generated, and requisite analyses were performed using an Abaqus plugin developed by Omairey et al. [40]. This results in 12 separate properties for the gyroid, however, there is little difference in the longitudinal moduli (E x , Ey, and E z ), or the shear moduli (G xy , G xz , and G yz ), or Poisson's ratios (ν xy , ν yx , ν xz , ν zx , ν yz , and ν zy ). The gyroid exhibits cubic symmetry, meaning that the elastic behavior is fully defined using three independent constants: E, ν, and G. Any difference between the different directional properties arises from mesh imperfection. For this reason, the several properties are averaged to one constant each (see (9)-(11)).
Additional discussion on the theory behind periodic boundary conditions is found in [41][42][43]. Figure 5 shows the deformation of a representative gyroid solid model under the different loading modes. It is clear that the stress is not distributed uniformly through the structure, but rather concentrates in certain locations. The bulk elastic properties, however, are calculated from the average stresses and strains of the entire RVE.

Compression Experiments
In order to validate the FEA predictions, a series of compression experiments were conducted. The samples tested are prismatic rectangular blocks 152.4 mm × 50.8 mm × 50.8 mm (6 × 2 × 2 ) printed of Ultimaker toughened PLA on an Ultimaker 5 machine using a 0.4 mm nozzle with layer height 0.2 mm, extruder-temperature 205 • C, bed-temperature 60 • C, and print-speed 70 mm/s. Density was adjusted using the "Infill Density" and "Infill Line Multiplier" parameters in the Cura slicer. Wall thickness and top/bottom thickness parameters set to 0 to print bare infill. All other parameters remain default for an Ultimaker 5 machine. The size was chosen in order to capture the bulk behavior of the gyroid infill better than the standard compression specimen (which is much smaller) due to the printed unit-cell size. Cell-size and wall-thickness are the two driving factors for gyroid relative density. Since wall-thickness is fixed by the printer's bead-width, the slicer adjusts relative density by changing the cell-size. The double-bead printed specimens at 25%, for example, are only 6 unit-cells wide. A smaller specimen would likely have significant edge-effects at this relative density. While the majority of these specimens are printed so that the testing direction corresponded with the z-axis of the printer, others are printed in the x-, and y-directions, respectively, in order to capture any directional effects. These specimens are labeled separately in the results. All of the specimens were tested using an Instron 100 kN servohydraulic test frame in quasi-static compression. Compression was chosen so as to capture longitudinal stresses while avoiding potential issues with gripping a porous structure. This test was based on ASTM D695 [44], but with the aforementioned deviation in specimen size. These experiments are displacement-controlled with a feed rate of 1.27 mm/min (0.05 in/min). In addition, any extraneous moments due to fixture misalignment are eliminated by using a swiveling top platen (see Figure 6).
Load data is captured using a 100 kN load cell and strain data is captured using digital image correlation (DIC). DIC data is captured in the form of a pointwise displacement field based on a series of dots adhered to the surface of each specimen. Analyzing this data using the GOM Correlate software, a surface is generated wherein the dots become the corner nodes of a mesh. The longitudinal and transverse strains are then calculated at each node using the built-in analysis tools. It is noteworthy that neither longitudinal nor transverse strain vary noticeably across the measured region (Figure 7 demonstrates how little the transverse strain varies), which indicates that the gauge region is sufficiently distant from the testing fixtures as to be outside the edge-effected area.  This test was performed with several specimens each at 25%, 50% , 75%, and 100% nominal relative density. For each relative density (except 100%) two different geometries are used: one-bead and two-bead wall thicknesses. This is in order to determine whether the wall-thicknesses are a significant factor in the response. Nominal density of each specimen is defined as a parameter in the slicing software, while the actual density is determined by measuring its size and mass. This is then divided by the published density of PLA to convert to relative density.

Mesh Convergence
In order to guarantee that mesh-size effects are not significant in the final results, a convergence study has been performed using the 50% relative density 3D model. Based on the results of this study (Figure 8), all subsequent analyses were performed using 10-noded quadratic tetrahedral elements. These elements have an average size that is 2% of the cell size. Average size is specified here because the models for different relative densities all have different numbers of elements, but the average size remains constant in all of the models.

Shell vs. Solid Models
In this section, the results obtained from the shell model are compared with those obtained from the solid model. Figure 9 shows that the shell model matches the solid for the lower relative density portion of the plot, but the two models rapidly diverge above 40% relative density. At these higher relative densities, the solid model more accurately matches the expected result by approaching 100% of the neat stiffness as it approaches 100% relative density. As such, the solid model results have been used for all comparisons to experimental data and for curve fitting. This additional accuracy does come at the cost of increased complexity, and computation time. In addition, to simulate different relative densities using solid models requires that the thickening process be performed to generate a new mesh for each relative density to be analyzed, where the shell-model's relative density can be adjusted by simply altering the shell-thickness parameter. For these reasons, it would be best to use shell models in any instance where their accuracy is sufficient such as on the low end of the relative density range.

Compression Results
Stress-strain curves are shown in Figure 10. These samples are clearly seen to exhibit elastic-plastic behavior similar to those samples of ULTEM 9085 tested by Bhandari and Lopez-Anido [12]. Based on qualitative observation of the many simulation results, this failure seems to be a combination of the local elastic instability/buckling and plastic failure. As a bulk property, however, this can be captured simply as an elastic-plastic response.  Table 1 shows a summary of results from the compression tests. From this data, we can see that the actual measured relative densities of 2-bead thickness specimens are much closer to the nominal relative density, especially for higher relative densities. For example, the single-bead 75% infill specimens are actually closer to 50% relative density, whereas the double-bead variant is exactly the nominal 75% relative density. For this reason, it is assumed that the printer cannot properly resolve the gyroid geometry using a single-bead. It is important to note that the relative density of the gyroid is determined by two factors: the unit-cell size, and the wall thickness. Since wall thickness is fixed here (by the beadwidth of the print), the way in which relative density is altered is by varying unit-cell size. This means that of the specimens printed, the 75% single-bead specimens have the smallest cell size, which explains why the resolution issue is most pronounced in these specimens. Elastic modulus was determined for each of the specimens by fitting a line to the initial linear portion of the stress-strain response. Figure 11 shows a comparison between experimental data and predicted values. On this plot, it is clear that, while there is some variation from specimen to specimen, the experimental data agrees well with that predicted. Note that due to the variation in results, these data were normalized against a reference material properties drawn from the PLA filament manufacturer's published modulus and density data (rather than against measured values from the 100% infill test-specimens).
It is also worth noting that for both modulus and Poisson's ratio, the z-orientation (default) specimens vary only slightly from their xand y-orientation counterparts at 50% relative density.

Poisson's Ratio
Poisson's ratio of the test specimens was determined by first averaging the strains across the gauge region of the specimen, then taking the ratio of average transverse strain and average longitudinal strain. In order to minimize data noise and capture the overall behavior through the linear portion of the test, this was calculated using a linear regression for¯ transverse vs.¯ longitudinal . Figure 12 shows the measured Poisson's ratio data points along with the FEA predicted ranges for the gyroid's Poisson response. Since published data for the Poisson's ratio of neat-PLA ranges from 0.3 to 0.35, it is anticipated that these gyroid specimens, being printed from PLA, should have a Poisson's ratio that falls within the 0.3 < ν m < 0.35 window. This data, however, shows a significant amount of scatter; enough that the variation between specimens is larger in some cases than the variation that was predicted for 0 < ν m < 0.45. In addition, it is clear that the double-bead thickness specimens match the FEA predictions much better than do the single-bead specimens. This, in addition to the issue with matching desired relative density, indicates that it is best to use a double-bead infill for structural predictability.
From this data, we can see that the modulus predictions match reality quite well, and the Poisson's ratio predictions match for specimens with 2-bead wall thicknesses.

Semi-Empirical Models
In Section 2, a method was presented for generating predictive solid models of the gyroid infill. These solid-model analyses were performed several times at each of several relative densities over a range of test materials with modulus values (E m ) from 1000 to 4000 MPa, and Poisson's ratio values (ν m ) ranging from 0 to 0.45. (Since the base material is isotropic, the shear modulus G m is calculated from them using the standard isotropic relationship.) This range of properties was chosen to encompass the behavior of a variety of common 3D printing polymers, so as to allow for predictions with all of these materials. In order to best use the resulting predictions for design, a relatively simple equation is defined which fits the data. This is similar to how simple mathematical models have been used to predict the behavior of other cellular structures such as honeycombs or foams [45].

Longitudinal Modulus
The starting point for this effort is to find a model that properly predicts the longitudinal modulus of the gyroid at various relative densities. By looking at a plot of the FEA results, it is clear that the gyroid's longitudinal modulus (E g ), when normalized against the neat material modulus (E m ), is nearly independent of all parameters except relative density. A variety of models have been found in the literature relating to cellular structures such as foams and honeycombs. For example, a power law is used to represent the modulus of closed-cell foams in [46]. In addition, a variety of polynomial ratios have been used for various foam and honeycomb predictions in [47][48][49][50][51]. None of these, however, fit the data satisfactorily. The model found to best fit the data was a modified power law, where E n is the normalized longitudinal modulus, E m is the neat-material elastic modulus, C i are curve fitting coefficients, and D is the relative density of infill (0 ≤ D ≤ 1 which corresponds to 0% ≤ D% ≤ 100%). The best value of fitting parameters C 1 and C 2 were determined to be 0.77 and 1.2, respectively. This results in an average error of 0.44% (maximum 1.13%). Figure 13 shows the FEA and semi-emperical predictions.

Shear Modulus
While shear modulus has not yet been verified experimentally, it was calculated in the FEA analysis so the results have also been curve fit for simplicity of use. Unlike the longitudinal modulus, the normalized shear modulus is seen to be a function of both relative density and Poisson's ratio. As such, very few different types of model have been found in cellular structures or composites literature which are applicable. Gong et al. [50] and Ghezal et al. [48] predict different responses for shear response, but neither of these results in a satisfactory fit for the gyroid. Likewise, simple modifications of the basic mathematical models yield unsatisfactory results. It is apparent, however, that the normalized shear modulus follows a similar curve to the longitudinal modulus. In fact, simply using the E/E m curve results in a max error of only 9.3%. This is improved significantly with the addition of a Poisson's ratio dependent term, yielding (13). The form of this term was determined by curve fitting various functional forms and choosing the one which resulted in the smallest error.
where G n is the normalized shear modulus, G m is the neat-material shear modulus, and ν m is the neat-material Poisson's ratio. The best value of fitting parameter C 3 was determined to be 0.5. This results in an average error of 0.35% (maximum 1.36%). Figure 14 shows the FEA and semi-emperical predictions.

Poisson's Ratio
Similar to the shear modulus, there are very few models in the literature to predict the bulk Poisson's ratio of cellular microstructures. By removing the purely geometric component of the response (Poisson's ratio resulting from a gyroid of zero-Poisson material), it was found that this also follows a similar curve to the longitudinal modulus. Figure 15 shows the FEA and semi-emperical predictions.
where the geometric component is represented by the polynomial: The best value of fitting parameters C 4 , C 5 , and C 6 were determined to be 0.28, 0.088, and 0.16, respectively. This produces an average error of 1.59% (maximum 6.54%).

Degree of Anisotropy
A common metric for determining the anisotropy of a material with cubic symmetry is the Zener index [52] (or Zener ratio), which is defined by (16).
Isotropic materials have Zener index of one. The further from one, the more anisotropic the material. The Zener index has been calculated for the gyroid using the semi-empirical equations developed previously, and the results plotted in Figure 16. The figure shows that the gyroid's Zener index varies with varying density, but that it remains close to one (0.95 ≤ A ≤ 1.25). By comparison, some other cubically symmetric infills have been shown to exhibit Zener ratios above two [53]. It is clear, then, that the gyroid is very-nearly isotropic. This near isotropy greatly simplifies implementation of gyroids into design and optimization workflows.

Conclusions
A novel homogenization method has been developed whereby FEA can be utilized to determine effective elastic properties of 3D printed infills. By discretizing the infill's unit cell and applying periodic boundaries, a prediction of the bulk elastic properties is made for a given relative density.
This method has been applied to the gyroid infill by performing a full suite of simulations. This allowed for generation of predictive data for the effective elastic properties at the full range of relative densities, using various base-materials. The FEA predictions have been verified using compression experiments whose results correspond well to the predicted behavior. At 50% nominal density, the error in longitudinal modulus is only 4% of neat modulus, which is less than the variation within the test results (5%). Finally, the FEA predictions have been used in order to inform semi-empirical equations, which can be used directly in designing parts using gyroid infill. These equations match the FEA predictions do a high degree of accuracy (0.44% error for the longitudinal modulus).
It is clear that a method of predicting homogeneous bulk properties of infills is an important component in optimal design using AM technology. The method presented here will likely be very useful in designing components using any of the complex 3D infills.