Optimizing the printing parameters for dimensional accuracy of distal femur bone by using Taguchi’s method

Background: Fused deposition modelling (FDM) is a popular additive manufacturing technique with capability of producing complex and integrate shapes. One of the critical aspects of FDM is the dimensional accuracy of 3D (three-dimension) printed model, especially in medical science applications, as proper fit and function with human body can prevent patient’s discomfort, complication or even harm. Objective: In this research work, the optimisation of print parameters: layer height, nozzle temperature, printing speed, infill pattern and infill density for improving the dimensional accuracy of distal femur bone, an irregular and complex shaped geometry is carried out using Taguchi’s method and to study its influence using ANOVA (analysis of variance). Methodology: 3D CAD (computer-aided design) model of the distal femur bone is generated from a CT (computerized tomography) scan using 3D slicer and its associated errors are corrected using Ansys SpaceClaim. The model is prepared for printing using Ultimaker Cura as per L 16 orthogonal array experimental layout where TEA (trans epicondylar axis), which is the distance between the most prominent point of the lateral and medial epicondyle, is set at 45° from X-axis in XY plane, i.e. diagonally on the plane of printing bed. It is then printed with PLA (polylactic acid) filament. Length along TEA is compared accordingly with 3D CAD model. Taguchi’s method of ‘smaller the better’ is applied for reducing deviation. Further, ANOVA analysis is done on the data set and a linear regression model is also developed. Result: Through Taguchi’s method, the optimum parameters were found to be triangle for infill pattern, 200 °C for nozzle temperature, 30 mm/s for nozzle speed, 0.1 mm for layer height and 40% for infill density. ANOVA analysis shows that all parameters contribute significantly with layer height being the most influential parameter, followed by infill pattern, nozzle speed, nozzle temperature and infill density. Mathematical model through multiple linear regression method was developed with determination of coefficient value of 96.91% and standard residual value is within the acceptable range of ± 2 indicating that there is no outliner in the data.


Introduction
FDM (fused deposition modelling) is a widely used additive manufacturing technique that has the capability of producing complex and integrated shapes [1].Its applicability can be found in various fields such as medicine, industry, agriculture, etc. [2][3][4].It works by melting thermoplastic filament and extruding it layer by layer through the nozzle to create the model as per 3D CAD (computer-aided drawing) model data [5,6].One of the critical aspects of FDM is the dimensional accuracy of the printed model especially in medical science applications [7,8].The dimensional accurateness of a model in general is to ensure proper fit and assembly, to maintain quality functional prototypes or products, especially in the field of medical science as FDM is used to create custom implants, prosthetics and anatomical models for surgical guides.A 3D print of custom anatomical model for implants or a surgical guide with accurate dimensions will ensure the safety of patient as a proper fit and function with the human body will avoid the patient's discomfort, complication or even harm [9][10][11].Our proposed research focuses on the dimensional accuracy of custom 3D printed distal femur as knee replacement surgery is a frequently performed procedure worldwide, with hundreds of thousands of surgeries conducted annually [12].With the help of the printed model of distal femur bone, customised femoral component of the knee implant can be manufactured which will be beneficial compared to the standard implants currently available in the markets by reducing the reduction of healthy bone.The printed model can also act as a surgical guide for complicated surgery [11].With the increase in the research of bioactive composite that can be 3D printed such as PLA-hydroxyapatite composite, PLA-chitosan composite, etc. [13] for the printing of customised implants, our research work might be used as a basis for printing with greater dimensional accuracy.
The dimensional deviation of the 3D printed model is affected by various factors such as shrinkage or wrapping, stepping effect, overhanging structure without support at a higher angle, build direction, etc. which are influenced by various 3D printing parameters [14][15][16][17][18].These parameters have been studied by various researchers to reduce the dimensional deviation of the 3D printed model.
Agarwal et al. [19] examined the print parameters' effect including print speed, wall thickness, infill density, layer thickness, build plate temperature and extrusion temperature on the dimensional accuracy of central composite design 3D printed model using ABS.Robles et al. [20] analysed the effect of layer thickness and infill density of FDM on dimensional accuracy of ABS parts by printing cubic shapes.Bolat et al. [21] analysed the layer height's effect on the dimensional accuracy of PLA, PET-G and ABS on 3D printed model of standard ASTM D638-14 type IV.Suaidi et al. [22] optimized the FDM printing parameters viz.layer thickness, infill density and build orientation for improving the surface roughness, dimensional accuracy and tensile strength using 3D printed ASTM D638-14 type IV model.Mohanty et al. [23] examined the effect of raster angle, part orientation, air gap, layer thickness and raster width on the dimensional accurateness of a 3D-printed geometrical model using ABS material.Akbas et al. [24] examined the influence of nozzle temperature and feed rates on the dimensional accuracy of PLA and ABS using strip as printed models.Zharylkassyn et al. [25] reviewed over 31 articles that examined the effect of layer thickness, part orientation and extrusion temperature in terms of the material used viz.ABS, PLA, etc. on the dimensional accurateness of printed model.
Various parameters and their significances concerning the dimensional accurateness of the printed model have been extensively examined through the use of standard geometrical 3D models such as cylinder, cuboid, ASTM D638-14 type IV model, etc., yet, no researchers have attempted to study an irregular and complex shape geometry such as human anatomy.Hence, our research focuses on improving the dimensional accuracy of 3D printed model of distal femur bone through the optimisation of five significant FDM printing parameters viz.layer height, nozzle temperature, infill density, printing speed and infill pattern with four levels each using Taguchi's method.Furthermore, the significance and its effect on dimensional accuracy are analysed using ANOVA.A linear regression model, using the experimental data, is developed to understand the correlation between the dimensional deviation and printing parameters.The linear regression model is validated by comparing it with the experimental data.Overall, this research will help in the fabrication of dimensionally accurate customised implants, implant mould and other complex and irregularly shaped geometry.

Methods
This research aims to optimize the FDM print parameters for dimensional accuracy of an irregular and complex shape using distal femur bone through Taguchi's method.A flowchart of the process employed for the research is shown in Fig. 1.The CT scan of a patient is imported into the 3D slicer software where 3D CAD model is exported in STL file format from DICOM file format.The errors associated with the STL file 3D model are corrected with the help of Ansys SpaceClaim.Once the correction is done to the STL file model, the model is then imported into the Ultimaker Cura, the model across which the length is to be measured i.e., trans epicondylar axis (TEA), which is the distance between the most prominent point of the lateral and medial epicondyle, is set at 45° from X-axis in XY plane i.e. diagonally on the plane of printing bed as shown in Fig. 2.
Polylactic acid (PLA) filament is considered as the printing material.PLA is one of the most commonly used materials because of its low cost, easily available, comparatively low toxicity and low print temperature [26,27].PLA is also a biocompatible material which researchers have been using PLA as base material for making bioactive composite material [13].The low thermal coefficient of PLA highly influences dimensional accuracy as it influences the expansion and shrinkage [28].Similarly, good intra-bonding of layers of PLA helps in maintaining the dimension of the model.The colour of the PLA filament and the fabrication process of PLA was also observed to be influencing the dimensional accuracy which is due to colour dependent on process crystallinity or difference in thermal properties due to the presence of colour additives [29][30][31].
A pilot test was conducted to determine the range of FDM printing parameters that affect the dimensional accuracy of the 3D-printed models.Accordingly, the printing parameters and their range that have been considered are infill pattern which is made to vary between lines, cubic, triangle and octet, infill density is varied between 20 to 80% with a step size of 20%, layer height is varied from 0.1 mm to 0.4 mm with a step size of    1.The parameters shown in Table 2 are fixed while the remaining parameters are set as per the default setting of the Ultimaker Cura.The experimental layout if used in full factorial design, the total experiment to be conducted will be 1024, i.e. 4 5 for 5 factors and 4 levels; however, through the use of Taguchi's method of orthogonal array, the total number of experiments can be reduced while giving statically same valid result.As the degree of freedom (DOF) in Taguchi's design is given by subtracting 1 from number of levels, the DOF for each factor with 4 level will be 3.Therefore, the total DOF for 5 factors with 4 level will be 15.The appropriate orthogonal array will be L 16 .Hence, an experimental layout of the print parameters as per L 16 orthogonal array is listed in Table 3.Additionally, to observe the improvement in accurateness of dimension, a model is prepared for printing as a baseline model with the default setting of ultimaker cura where infill pattern is cubic shaped, layer height is 0.2 mm, infill density of 20%, nozzle temperature of 200 °C and with a print speed of 50 mm/s.All these parameters are then set accordingly and slicing is done.The G code is then transferred to the FDM printer (Creality Ender 3 Pro) and models are then printed accordingly using PLA filament of 1.75 mm diameter as printing material.A digital vernier calliper with a least count of 0.01 mm is used for measuring the 3D printed distal femur bone across the trans epicondylar axis (TEA).The measurement is taken thrice to reduce human errors in the measurement.The measurement is then compared with the TEA length measured in 3D CAD model using Ansys SpaceClaim and absolute deviation is calculated using Eq. 1. ( Absolute Deviation = |TEA length of 3D printed model−TEA length of 3D CAD model|

Results and discussion
Modelling and printing of distal femur bone 3D model of the distal femur bone as shown in Fig. 3 is generated from CT scan to STL file format using 3D slicer.The final model generated from Ansys SpaceClaim after the correction of STL file format has about 974,096 faces and 487,050 vertices as shown in Fig. 4. The model is then sliced into individual layers and converted the STL file into G code using Ultimaker Cura.The G code makes the tools in the 3D printer move as per the generated layer information during the slicing process.A total of 16 models are printed accordingly as shown in Fig. 5.

Optimisation of the printing parameters for dimensional accuracy
The length of TEA in 3D CAD model, when measured using Ansys SpaceClaim, is 79.72 mm.The absolute deviation of the printed model from that of the 3D CAD model is calculated using Eq. 1 and is listed in Table 4.The model printed as baseline was found to deviate with an absolute mean value of 0.2933 mm.As lower deviation of the model is required to increase the dimensional accuracy of the model, Taguchi's  'smaller the better' characteristic is taken into consideration whose value is given by the formula: where 'Y' is the result value, i.e. absolute deviation and 'n' is the number of experiments.
The absolute deviation along with the value for signal-to-noise ratio (S/N) calculated using Eq. 2 is listed in Table 4.The S/N ratio for each level base on the smaller is better characteristic is shown in Table 5 which identifies the optimum parameter, i.e. parameter with the highest S/N ratio for lowering of absolute deviation and the influence of each factor has been ranked accordingly as per the Delta value.Therefore, as per the mentioned values in Table 5 and the graph shown in Fig. 6, the optimal parameter for higher dimensional accuracy is when the infill pattern is triangular, nozzle temperature is 200 °C, printing speed is 30 mm/s, layer height is 0.1 mm and infill density is 40%.Fig. 6 Signal-to-noise ratios graph for each parameter using smaller is better characteristic The data is further analysed using ANOVA whose result is shown in Table 6.It can be seen that all values of p are very close to 0; hence, the null hypothesis can be rejected.This indicates that all the printing parameters are of very significant in relation to the dimension accuracy of the 3D printed model.The layer height has the highest contribution of 55.47% in regards to improvement in the dimensional accuracy, followed by an infill pattern of 24.75%, then by printing speed, nozzle temperature and infill density.The infill density has the least contribution in regards to dimensional accuracy.Further, the linear regression analysis is performed on the data set which generates an empirical equation for a predictive mathematical model as shown below in Eq. 3 with determination of coefficient, R2 value of 96.91%.The empirical equation, i.e. linear regression model helps in establishing a relation between the dependent variable, i.e. absolute deviation and independent variables which are the 3D printing parameters.To validate the empirical equation, a comparison is made between the absolute mean deviation and calculated deviation (fits) as shown in a graph in Fig. 7 and Table 7 which shows that both are in agreement with a little difference between them.The normal probability graph shown in Fig. 8 shows that the standard residual value is within the acceptable range of ± 2 indicating that there is no outliner in the data.

Validation of the result
Layer height highly influences the dimension accuracy of the 3D printed distal femur bone as seen from ANOVA analysis in Table 6.The optimal dimension is achieved when the layer height is 0.1 mm.With an increase in the layer height, it is observed that the dimensional deviation increased as seen in the S/N graph shown in Fig. 6 as PLA is deposited from the nozzle in the molten state; due to its fluidity, it bulges out which results in dimensional deviation.The increase in layer height increases the volume of the material which in turn increases the bulging of material resulting in greater dimensional deviation.The dimensional deviation is further amplified by the staircase effect [32].This makes the layer height the most significant parameter among the other parameters in controlling dimensional accuracy.The result observed in this research is similar to the finding of Nidagundi et al. [33] and Polak et al. [26].Among the four infill pattern types, the optimal infill pattern is observed with triangle infill pattern which is similar to the finding of Alafaghani et al. [34] where the dimensional accuracy is observed with diamond infill pattern.The diamond infill pattern is similar to the triangular infill as combining the base of two triangles forms the diamond infill.The dimensional accuracy decreases with linear infill pattern followed by cubic infill pattern.The highest dimensional deviation is observed with octet infill pattern.
The dimensional deviation of the distal femur bone increases as the printing speed increases from 30 mm/s to 150 mm/s as the time available for proper deposition of the material decreases.Hence, the optimal printing speed for dimensional accuracy is observed at 30 mm/s.It is in sync with the findings of Polak et al. [26] and Alafaghani et al. [34].
Optimal nozzle temperature for dimensional accuracy of the printed distal femur bone is observed at 200 °C.The dimensional deviation increases as the nozzle temperature increases from 200 °C to 215 °C.This is due to an increase in the fluidity of PLA with increasing temperature thereby making the dimensional control difficult as PLA expands freely [34][35][36].The observation made in this research is in line with the findings of Akbas et al. [24], Kaveh et al. [37] and Frunzaverde et al. [27].
Even though infill density is significant concerning dimensional accuracy, its contribution among the other control print parameters is less.The optimal infill density for dimensional accuracy of the printed distal femur bone is observed at 40%.The accuracy increases as the infill density increases from 20 to 40% but the dimensional deviation increases as the density increases from 40 to 80% with a slight variation of dimensional difference between the printed distal femur bone with 40% and 60%.This observation is similar to the finding of Alafaghani et al. [34] and Suaidi et al. [22].The increase in accuracy with an increase in infill density is due to the presence of high volume of molten material and lesser space within the printed model part causing it to flow outwards which counteract the dimensional error due to the shrinkage [18].However, with further increase in the infill density, dimensional error increases as internal stress is created within the part due to uneven heat distribution during cooling and solidification with increases in density as the shrinkage cycle of the molten material increases resulting in wrapping and distortion making it difficult to maintain or control dimensional accuracy [19].
The observations made in this research are validated as it is in line with the finding of various researchers.

Optimal 3D printing parameter
The obtained optimum 3D printing parameter using Taguchi's method of smaller the better characteristic as shown in Table 8 is used for printing the model.The TEA of 3D printed distal femur bone is measured using a digital vernier calliper.Its absolute mean deviation is also compared with that of the deviation calculated (fits) using the empirical equation in the Eq. 3 as shown in Table 9.The value obtained is within the 95% confidence interval of between 0.241035 and 0.260465.Overall, it can be seen that the obtained absolute mean deviation with the optimum printing parameters is less than the other deviation.The printed model using the optimum parameter is shown in Fig. 9.

Overall evaluation of process parameters
Figure 10 shows a graph that compares the dimensional accuracy of 16 sets of printed models with that of model printed with optimum parameters and baseline model.When compared to the baseline model that is printed with the default setting of Ultimaker Cura with that of the optimum parameters, dimensional accuracy is improved by 17.24%, i.e. by 0.05 mm.The difference between them is less as the layer height of baseline model is 0.1 mm more than that of the optimum, the difference in speed is 20 mm/s and the density of the baseline model is less by 20% while the temperature is same in both cases.As layer height and infill pattern contribute significantly more than others, the dimensional difference in between the baseline and optimum model did not differ significantly.This can be seen from Table 10 that in model 1, when the infill pattern and layer height are maintained same with that of the optimum parameters, the mean dimensional deviation is of 0.25667 which is closer to the dimension deviation of the optimum model whereas in model 2, when the infill pattern and layer height is maintained same with that of the baseline parameters, the mean dimensional deviation is of 0.27333 which is closer to the dimension deviation of the baseline model.The difference in the dimensional deviation of model 1 and 2 also iterates the significant contribution of layer height and infill pattern.The 3D printed distal femur bone of baseline, model 1 and model 2 is shown in Fig. 11.
The limitation of our research is the unable to control the environmental factors such as ambient temperature, humidity, etc. which can affect the material, printing process and, consequently, dimensional accurateness of the model [35].

Conclusions
This research paper focuses on the optimisation of the important printing parameters viz.infill pattern, nozzle temperature, printing speed, layer height and infill density for improvement in dimensional accuracy of the irregular and complex shape of human anatomy (distal femur bone) when measured along the TEA which was aligned at a 45°   • Using the 'smaller the better' characteristic of Taguchi's method, the minimum dimensional deviation was observed when infill pattern is Triangle; nozzle temperature is at 200 °C; nozzle speed is 30 mm/s; layer height is 0.1 mm; and infill density is of 40%.• ANOVA analysis of the data has shown that all the parameters contribute significantly as the p-value is very close to 0. The dimensional accuracy of the distal femur bone was highly influenced by layer height with a contribution of 55.47% followed infill pattern of 24.75%, then by printing speed, nozzle temperature and infill density.• A close agreement was observed between the fits and the experimental absolute deviation in the developed linear regression model.This empirical equation can be used for the prediction of absolute deviation for the given printing parameters.

Fig. 2
Fig. 2 TEA of distal femur bone set at 45° angle from X-axis in Ultimaker Cura

Fig. 3 Fig. 4
Fig. 3 3D CAD model of distal femur bone in STL file format generated from CT scan using 3D slicer software

Fig. 5
Fig. 5 3D printed distal femur bone as per the experimental layout of L 16 orthogonal array

Fig. 9
Fig. 9 3D printed model of the distal femur bone using optimum parameter

Fig. 10
Fig. 10 Comparative graph of absolute mean deviation of the 3D printed model of distal femur

Table 1
Control FDM printing parameters with 4 levels each

Table 2
Fixed FDM printing parameters , nozzle temperature is varied from 200 °C to 215 °C with a step size of 5 °C and print speed is varied from 30 mm/s to 150 mm/s with a step size of 40 mm/s shown in Table

Table 3
Experimental layout of the print parameters using L 16 orthogonal array

Table 4
Absolute dimensional deviation of TEA of 3D printed distal femur bone with S/N ratio and mean

Table 5
Signal-to-noise ratio for each level of each parameter along with delta value

Table 6
ANOVA analysis of the result Fig. 7 Comparative graph between the measured absolute deviation and fits Singh et al.Journal of Engineering and Applied Science (2024) 71:10

Table 7
A comparative table between mean absolute deviation and fits calculated using linear regression model Normal probability graph plot for standardized residual

Table 8
Optimum 3D printing parameter with respect to dimensional accuracy

Table 9
Dimensional deviation of the 3D printed model using optimal value

Table 10
3D printing parameters for models 1 and 2