Computational analysis on mechanostructural properties of hydroxyapatite–alumina–titanium nanocomposite

In this research, Taguchi–grey relational analysis has been applied to mitigate the insufficient assumptions made on the optimization of mechanical and structural (mechanostructural) properties of synthesized hydroxyapatite (HAp)–alumina–titanium nanocomposite. This nanocomposite has already been developed and studied in the previous study. This paper employs the L9 (3**3) orthogonal array, including displaying factors and levels of 3, 5, 7 wt % for alumina, 5, 10, 15 wt % for titanium, and 1100, 1150, 1200 °C sintering temperature. The computational analysis presents the predicted mechanostructural grey relational response as 0.7271, higher than the highest response shown in the ninth experimental run. The optimal control factors are analyzed to be 7 wt % alumina, 15 wt % titanium, and 1200 °C sintering temperature. The obtained result elucidates the hypothesis that a singular response optimization is not enough in the fabrication of biomedical material, disproving the assumption made in the previous literature. Importantly, to fabricate a high clinical grade HAp–alumina–titanium nanocomposite, titanium is the most invaluable contributor with a contribution of 49.11%, followed by alumina (45.52%), and then sintering temperature (3.2%). Although the confidence level and probability distribution analysis show that all the experimental mechanostructural responses were within the 95% confidence level, the employment of the predicted optimal factors is strongly recommended for experimentation.


Introduction
Hydroxyapatite of biogenic or synthetic origin has extensive applications in biomedical repairs, substitutions, and augmentations and as scaffolds in tissue engineering for bone regeneration [1][2][3].HAp can be applied as abrasives to roughen metal implant surfaces and as a bioactive coating on orthopedic and dental implants [4][5][6][7][8].This material has also been used as transfection agents, drug carriers, and percutaneous devices.In the same vein, HAp is disadvantaged as it relates to its fracture toughness which somewhat restricts its use for load-bearing clinical applications.Therefore, extensive studies are geared toward methods that have the potential of enhancing the mechanical properties of HAp without overly compromising its biocompatibility.These methods stem from tailoring the processing of the HAp powders, carefully optimizing the sintering temperature, pressure, and reinforcements for better biomedical performance of HAp [9][10][11].
Alumina has received great attention as a promising material in improving the osteoconductivity of hydroxyapatite.Because of the bioinert nature of the material, it allows the bone and tissues surrounding it to re-grow and integrate without causing any negative immune response.In dentistry, bioinert materials are an essential component of implants, bone grafts, prostheses, and fillings.In addition to not harming the body, a bioinert substance in dentistry will also be osteoconductive and promote osteogenesis.These attributes show the fact that the material possesses unique properties like nontoxicity, excellent biocompatibility, and in-vivo biodegradability [12,13].Alumina-based materials find applications most often as bone substitute material, implant coats, and drug delivery systems [14][15][16][17].One of the major benefits offered by alumina-based biomaterials is their ability to produce the apatite layer when immersed in a simulated body fluid.This action plays a crucial part in inducing bonding between the implant and human bone [17].Many biological studies involving alumina-based materials have demonstrated that these materials have superior mechanical properties, and can enhance the rate and quality of bone tissue repair compared with pure HAp [18].Titanium and titanium alloys have been shown in several studies to have a distinctive combination of strength and biocompatibility, which makes them useful in medical applications.They have gained ground in biomedical engineering because of their exceptional properties, namely, fluid effects resistance, high tensile strength, high toughness, higher bending strength, high corrosion resistance, and high flexibility [19][20][21][22][23]. Hence, a mixture of alumina and titanium with an effective sintering technique tends to improve the mechanical and biological properties of HAp.Importantly, careful optimization of reinforcements and sintering temperature will give an excellent mechanically and biologically enhanced HAp suitable for tissue engineering.
Taguchi experimental design technique is an efficient technique that ensures excellent performance in the design of processes or products.The technique can reduce the number of experiments for better products or process performance, thereby reducing cost and mitigating errors [24].Despite the advantage of Taguchi experimental design, it can only optimize singular performance characteristic, which limits it for optimizing multiple performance characteristics [25].However, grey computational analysis can convert nonadditive multiple performance characteristics into a singular characteristic, thereby assisting Taguchi experimental design to optimize multiple responses of products or processes [26,27].In this recent time, researchers have started exploring the assistance of grey computational analysis with Taguchi design experimental design technique to optimize multiple response of a product or process [23,28,29].
Hence, in this study, Taguchi experimental design technique assisted by grey computational analysis was employed to optimize and to investigate the effect of alumina, titanium, and sintering temperature on mechanostructural properties of hydroxyapatite.

Experimental data curation and orthogonal array experiment
The synthesis, characterization, and mechanostructural properties evaluation of hydroxyapatite-alumina-titanium nanocomposite have been conducted in the study of Ebrahimi et al. [30].The L9 (3**3) Taguchi orthogonal array employed in this study agrees with the previous study, having the control factors and levels, shown in Table 1.This design resulted to only nine experimental runs, shown in Table 2.The mechanostructural experimental data employed for the Taguchi-grey computational analysis was obtained from the previous report, and it is displayed in Table 3.

Mechanostructural grey computational analysis
As established in the grey computational analysis technique, the five experimental mechanostructural data were first preprocessed to range within zero and one.The preprocessing data technique used for phase decomposition, open porosity, toughness, and bending strength was the higher-the-better (Eq.1), while for the grain size was the smaller-the-better (Eq.2).The reason for the choice is that a higher porosity, higher phase decomposition rate, and smaller grain size enhance cell proliferation and surface growth [2,[31][32][33][34], while a higher mechanical strength is desired for load-bearing biomedical application.The preprocessed data are displayed in Table 4, and it is usually called grey relational generation.A comparison was done with an ideal sequence, x 0 (k) (k = 1, 2,…,9) for the five mechanostructural data.
where x i (k) is the preprocessed data for the ith experiment, and y i (k) is the initial sequence of the mean of the responses.The deviation sequence (Eq.3) was then computed to process the grey relational coefficient.The deviation sequence of the five experimental data is displayed in Table 5. (1) where Δ oi (k) , x o (k), and x i (k) are the deviation, reference, and comparability sequences, respectively.Then the grey relational coefficient (GRC) was computed using Eq. 4. The grey relational coefficient shows the relationship between the expected and obtained experimental data.It also gives the possibility to relate and carry out statistical mean of the non-additive multiple experimental data.
where i (k) reflects the GRC value of the individual mecha- nostructural responses computed as a function of Δ min and Δ max , the minimum and the maximum deviations of each response variable.is the distinguishing coefficient (0-1), but a different weight of 0.5 is usually assigned to each parameter.
Next, the grey relational grade, i.e., GRG in Eq. 5 was computed by averaging the five grey relational coefficients.The GRC and GRG values are displayed in Table 6.The grey relational grade gives the overall response of the ( 4) five experimental mechanostructural data.In other words, the unfeasible optimization of the five complicated data becomes feasible by the singular response (GRG) generated by the grey computational analysis.Optimization of the materials design parameters is the level showing the highest GRG.
where i = the value of GRG determined for the ith experiment, n is the aggregate count of the performance characteristics.

Taguchi design of experiment (DOE)
The obtained GRG values were subjected to the Minitab software and analyzed according to the considered orthogonal array design.Then the signal-to-noise ratio corresponding to the respective GRG values is displayed in Fig. 1.Table 7 displays the response table for GRG.The respective delta value shown in Table 7 is the difference between the maximum and the minimum value of GRG of each factor, ( and it is ranked.The ranking shows the order of significance on the resultant GRG values. Figure 2 displays the effect of control factors on the mechanostructural response (GRG) of the composite.It is seen that the increase in factor level increased the mechanostructural response.This shows that the optimal mehanostructural  response is derived with factors at level 3, i.e., 7 wt % alumina, 15 wt % titanium, and 1200 °C sintering temperature are the optimal control factors for better mechanostructural properties of the composite.These findings disagree with the conclusive optimal control factors in the study of Ebrahimi et al. [30].The reason for the disagreement is that Ebrahimi et al. [30] assumed only one response (structural property) to be the most desired property, and then singularly optimizing it for composite production.Whereas all the properties are essential for a functional biomedical material (See Sect.3.0 for the importance of each factor).The grey computational analysis has been able to optimize all the five complicated mechanostructural data for better biomedical performance.

Analysis of variance
Analysis of variance (ANOVA) helps to investigate the level of significance of each control factor on the mechanostructural response of the composite.important to note that residual error is insignificant, having a contribution of 3.2%.

Validation analysis
Having known the optimal and the significance of the composite fabrication parameters, it is expedient to investigate the predicted response from the design analysis.This is shown in Eq. 6: [23,35].
0 is the highest average value of GRG at the optimal fabrication parameter levels and m is the average of GRG.q is the number of the fabrication parameters.
In Eq. 5 and Table 6, the value of GRG was predicted using the optimal fabrication parameter levels.The predicted GRG was gotten to be 0.7271.
Confidence interval (CI) was then computed using Eq. 7.This was done to determine the correlation between the experimental GRG value and the predicted GRG value [23,36]: ratio required for ; = risk; f e = DOF of error; V e = variance of error; eff = effective number of replications, which is Eq. 8 given below: R = number of replications when the experiment is carried out for confirmation; N = total number of experiments.
Therefore, V e = 0.000829; f e = 2.Total DOF of control factors = 6.R = 1, N = 9. = 0.5 (95% confidence interval).F 0.5 (1, 2) = 18.51 (tabulated values from the F-Tables ) 95% confidence interval of the predicted optimal grey relational grade is given in Eq. 9, as posited by Abifarin (2021): The optimal fabrication parameters predicted in this study happened to be outside the experimental run in the study.There is a need to fabricate a composite based on the optimal fabrication parameters and then conduct its mechanostructural properties, thereafter benchmarking it with the generated confidence interval having the minimum and maximum mechanostructural response boundaries of 0.5619 and 0.8923.Putting all the fabrication conditions, it is expected that the mechanostructural grey relational response should correlate with the predicted response.
Having recommended a need to conduct experiment based on the predicted optimal conditions, probability distribution is another way of validating the predicted result.This is shown in Fig. 3.The experimental mechanostructural grey relational responses displayed in Fig. 3 are within the 95% confidence level.In other words, the optimal mechanostructural grey relational response is expected to correlate with the predicted mechanostructural response under 95% confidence level.

Multi-mechanostructural modeling and interaction plot of the composite
The experimental and modeled mechanostructural response of the composite is shown in Fig. 4a.The graph shows a very similar and close pattern between the experimental and the modeled mechanostructural response as generated by regression analysis on the GRG values obtained.Also, the mathematical model of the mechanostructural response of the composite is shown in Eq. (10). Figure 4b-d displays the reactions between the considered factors relative to their corresponding mechanostructural response.It is shown from the plot, the combination of factors settings resulting to as high as possible the mechanostructural response of the composite.
where A, B, and C denote the responses generated by alumina, titanium, and sintering temperature, respectively.

Conclusion
Taguchi-grey relational analysis for the optimization of control factors of hydroxyapatite-alumina-titanium composite synthesis for complicated multiple mechanostructural properties has been reported in this paper.The highest GRG value, which was 0.7010 for the nine experimental run, was obtained at the ninth experimental run (7 wt % alumina, 15 wt% titanium, and 1150 °C sintering temperature), but was proven not to be the actual optimal conditions.The best optimal control factors are validated to be 7 wt % alumina, 15 wt % titanium, and 1200 °C sintering temperature.These findings validated the research hypothesis that a singular response optimization is not enough for a biomedical material, disproving the assumption of Ebrahimi et al. [30].Quantitatively, titanium displayed the highest significance of 49.11%, followed by alumina (45.52%), and then sintering temperature (3.2%).From the study, the mathematical and graphical model displaying similar pattern of the experimental mechanostructural response shows the efficacy of the design analysis.

M
Ma ai i n n E Ef ff fe ec ct ts s P P l l o ot t f fo or r M Me ea an ns s Data Means

Table 1
Control factors and levels for composite synthesis

Table 2
Experimental design for composite synthesis

Table 5
Deviation sequence Experimental run Table 8 displays the ANOVA of GRG.The table was generated by Minitab 16 software considering the resultant GRG values, but the percentage of contribution was calculated to quantitatively know the level of significance of each control factor.A higher F value and a lower P show a greater level of significance.The result shows titanium has the highest percentage of contribution of 49.11%, followed by alumina and then a very little contribution from sintering temperature.It is

Table 7
GRG response table

Table 8
ANOVA of GRG