Numerical simulation of 3D printed FFF parts through G-code derived geometry
Simulating 3D printed parts obtained through Fused Filament Fabrication (FFF) technologies presents a significant challenge due to the lack of tools to accurately reconstruct the submillimeter inner geometry resulting from both infill and perimeter strategies. Researchers at ProtoLab, F. Rotini and S. Marconi addressed this challenge by developing an in-house MATLAB code.
This code generates precise interior (infill) and exterior (perimeters) geometries for any FFF components based on G-code instructions.
The underlying concept originates from the original idea of G. Alaimo and A. Cattenone, and the project received supervision and support from Prof. F. Auricchio.
Methods
Ten beam samples were 3D printed and tested by B. Rossetti and F. Rotini using the Litem BYC 200 machine, available at the Department of Civil Engineering and Architecture (DICAR) of Pavia University.
The experimental setup is shown in Figure 1a. To run the test, a total displacement of 4 mm was applied to the specimen midpoint at a deformation rate of 1.7 mm/min according to the ASTM D790 standard. The deformed configuration of a beam sample is shown in Figure 1b.
![litem setup fff litem setup fff](https://compmech.unipv.it/wp-content/uploads/2024/07/ImportedPhoto.743592843.95764-scaled-500x500.jpeg)
Figure 1a: three point bending test experimental setup
![3D printed giroide Bending test 3D printed giroide Bending test](https://compmech.unipv.it/wp-content/uploads/2024/07/ImportedPhoto.743592863.884121-scaled-500x500.jpeg)
Figure 1b: beam deformed configuration
Figure 2 shows the averaged load-displacement curve of the experimental test, along with the 98% confidence interval.
![experimental data](https://compmech.unipv.it/wp-content/uploads/2024/07/c98dbb41ccb3f001db2eb3b28b3671e5.png)
Figure 2 Average load-displacement curve, along with the standard deviation
The geometry of the samples was reconstructed through the proposed Matlab code, whereas the G-code visualization through slicing software is shown in figure.3
![](https://compmech.unipv.it/wp-content/uploads/2024/07/06242e7f68b5f676e51f42abd6f4086b.png)
Figura 3: G-code geometry reconstruction
The experimental three-point bending test has been replicated through the numerical simulation carried out in Ansys, with frictionless contact and large deformations. The material model is linear elastic, with material parameters as per manufacturer’s data sheet. Set-up and meshed geometry are shown in figure 4 and figure 5 whereas total deformation and Y – normal stress resulting from the FE analysis are shown in figure 6a and 6b, respectively.
![](https://compmech.unipv.it/wp-content/uploads/2024/08/con-logo-2.jpg)
Figure 4 Numerical simulation set-up
![Ansys meshed geometry](https://compmech.unipv.it/wp-content/uploads/2024/08/235.gif)
Figure 5 Meshed geometry with details
![](https://compmech.unipv.it/wp-content/uploads/2024/07/fff2.jpg)
Figure 6a
![](https://compmech.unipv.it/wp-content/uploads/2024/07/fig-45.jpg)
Figure 6b
Results
Experimental and numerical data comparison is shown in Figure 7. The analysis of the figure shows that the numerical data falls within the confidence interval of the experimental data. Such an evidence allows us to conclude that the geometry reconstruction is accurate enough for a reliable prediction of the mechanical response of the gyroid beam samples.
![](https://compmech.unipv.it/wp-content/uploads/2024/07/9b4d369f8876e44f93d91cffdb3e5801.png)
Figure 7 Comparison between experimental and numerical data