Next Article in Journal
Outcomes and Complications of 33 Soft-Milled Cobalt-Chromium-Ceramic Full-Arch Screw-Retained Implant-Supported Prostheses: A Retrospective Study with up to 10-Year Follow-Up
Next Article in Special Issue
Enhancing Polymethyl Methacrylate Prostheses for Cranioplasty with Ti mesh Inlays
Previous Article in Journal
TiNi-Based Material with Shape-Memory Effect for Surgical Treatment of Diseases of Small Intestine in Newborn and Young Children
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Machine-Learning-Based Approach for Predicting Mechanical Performance of Semi-Porous Hip Stems

1
Department of Engineering Management, Prince Sultan University, P.O. Box 66833, Riyadh 11586, Saudi Arabia
2
Department for Mechanical and Industrial Engineering, Qatar University, Doha 2713, Qatar
*
Author to whom correspondence should be addressed.
J. Funct. Biomater. 2023, 14(3), 156; https://doi.org/10.3390/jfb14030156
Submission received: 8 February 2023 / Revised: 7 March 2023 / Accepted: 13 March 2023 / Published: 15 March 2023

Abstract

:
Novel designs of porous and semi-porous hip stems attempt to alleviate complications such as aseptic loosening, stress shielding, and eventual implant failure. Various designs of hip stems are modeled to simulate biomechanical performance using finite element analysis; however, these models are computationally expensive. Therefore, the machine learning approach is incorporated with simulated data to predict the new biomechanical performance of new designs of hip stems. Six types of algorithms based on machine learning were employed to validate the simulated results of finite element analysis. Afterwards, new designs of semi-porous stems with outer dense layers of 2.5 and 3 mm and porosities of 10–80% were used to predict the stiffness of the stems, stresses in outer dense layers, stresses in porous sections, and factor of safety under physiological loads using machine learning algorithms. It was determined that decision tree regression is the top-performing machine learning algorithm as per the used simulation data in terms of the validation mean absolute percentage error which equals 19.62%. It was also found that ridge regression produces the most consistent test set trend as compared with the original simulated finite element analysis results despite relying on a relatively small data set. These predicted results employing trained algorithms provided the understanding that changing the design parameters of semi-porous stems affects the biomechanical performance without carrying out finite element analysis.

1. Introduction

Conventional hip stems are made of dense metals such as cobalt-chromium, stainless steel, and titanium. The stiffness of these dense stems is 5–15 fold higher than that of the cortical bone and 50–100 fold higher than that of the cancellous bone [1,2]. This high mismatch of stiffness causes several complications such as stress shielding, aseptic loosening, corrosion, and implant failure [3,4]. Thus, revision of surgeries is required to fix the failed stems, which is a costly and painful process. To overcome the aforementioned complications, researchers have developed fascinating materials and designs using the additive manufacturing technique to reduce the stiffness of stems that can last longer with excellent functionality [5,6]. The choices of materials are limited due to biocompatibility; however, the implant designs have been extensively studied in recent decades. Several studies have shown excellent biocompatibility of titanium alloys when tested in vivo. Titanium alloys are inert and showed excellent bone ingrowth into the porous surfaces [7,8,9]. Moreover, advanced manufacturing techniques such as additive-manufactured implants also showed excellent biocompatibility; however, the development of new materials and designs of medical devices does not contain biocompatibility testing until the prototype stage is reached [8]. Recently, different stem designs were based on the additive manufacturing concept that provides wide freedom to alter the design [10,11,12]. Porous stems with various architectures are manufactured with additive manufacturing to reduce the stiffness of the stems. These porous cellular architectures include circular, cubic, body-centered cubic, diamond, and gyroid shapes [13]. Moreover, the porosity of these cells can be easily controlled within stems. These porous structures are usually filled inside the stem with an outer dense layer; thus, these design parameters generate a lot of design alternatives that are costly and time-consuming, using experiments. Therefore, computer models [14,15] are attracting the attention of researchers to reduce the cost, time, and uncertainty in experimental work. Alkhatib et al. [14] investigated the biomechanical performance of hip stems under various physiological conditions such as walking and stair climbing, and found that hip stems showed more stresses in the climbing activity as compared with the walking activity. Moreover, porous hip stems alleviated stress shielding in surrounding bone as compared with dense stems under different activities. However, these simulations are also costly, as they require the running of hundreds of simulations for these design parameters. Thus, researchers have recently given more attention to incorporating machine learning (ML) techniques in the medical field, which are capable of predicting the results [16,17,18].
Machine learning extends to a large variety of fields and can assume different approaches. However, one simple and clear definition was presented by Murphy [19] stating that machine learning is “a set of methods that can automatically detect patterns in data, and then use the uncovered patterns to predict future data”. In other words, machine learning techniques are statistical models that can be autonomously trained to predict data using complex algorithms [20]. Since the last decade, the technique has grown tremendously and is now being used in many fields such as industrial [21], project management [22], finance [23,24], construction and materials [25,26], and medicine [27], among others, to predict multiple types of data. It is worth mentioning that machine learning methods can be divided into classifiers and predictors. Classifiers are used for categorical variables whereas predictors are used for numerical variables, the latter being the case in this paper.
In medicine, several researchers have turned to ML techniques as a computationally cheaper alternative to Finite Element (FE) model simulations. Villamor et al. [28] made comparisons between Support Vector Machine (SVM), Logistic Regression, Shallow Neural Networks, and Random Forest ML methods to determine the best-performing model to predict osteoporotic hip fracture in postmenopausal women based on FE analyses. Alastruey-López et al. [29] used Artificial Neural Networks (ANN) and a parametric FE simulation to predict impingement and dislocation in total hip arthroplasty. Their efforts were aimed at identifying the optimal prosthesis design to reduce the probability of dislocation. Similarly, Jun et al. [30] used results produced using an FE model analysis to train a machine learning method that combines both principal component analysis (PCA) and support vector regression (SVR) in an effort to predict the contact stress of the hip prosthesis acetabular lining. The prediction model performance was then compared with the ridge regression and lasso models for validation. Cilla et al. [20] also combined FE modeling, ANNs, and SVMs in an effort to optimize the commercial short-stem hip prosthesis design. Their work focused on predicting the optimal stem length, thickness in the lateral and medial, and the distance between the implant neck and the central stem surface.
Within this context, this study aims to address stress shielding in dense hip implants by introducing a porous hip implant. ML techniques were used to validate finite element results of porous hip stems with different designs. Six types of algorithms were used to investigate errors. Then, a new data set was created with semi-porous hip stems which have outer dense layers with different thicknesses and inner porous cellular structures with different porosities. The most efficient validated algorithm was used to predict the outcomes of new implant designs, which are stiffness of the implants, stresses in dense layers and porous cellular structures, and factor of safety of the implants.
There are multiple stem designs that use finite element analysis (FEA) to investigate biomechanical performance, as mentioned above. One of these unique designs that help alleviate stress shielding was proposed by Mehboob et al. [12]. Accordingly, this paper contributes to the existing literature as follows:
  • Using predictive machine learning techniques to validate the FEA-based models presented by Mehboob et al. [12] to reduce the in vivo experimental cost.
  • Comparing multiple machine learning algorithms to determine the best-performing method for the chosen model.
Quick prediction of the designs of the stems using machine learning will be readily available for printing as a personalized implant which will decrease the complications of revision of surgery and reduce the burden on the health system. In addition, the lower cost achieved through machine-learning-aided designs will potentially reduce the cost of the implant on patients. Moreover, designing the semi-porous implants and printing using additive manufacturing will further reduce the wastage of the material.
The remainder of this paper is organized as follows. Section 2 explains the methodology used in the development of the hip implant through a finite element analysis and the training data set, as well as the methodology for developing the predictive models using machine learning techniques. Section 3 presents and discusses the results obtained from all six machine learning models in addition to identifying the best-performing model as well as an interpretation of these results. The final section presents an overall conclusion in addition to some limitations and avenues for future research.

2. Methodology

2.1. Finite Element Analysis

Finite element models of various designs of hip implants and bone were constructed in SolidWorks in a previous study [12]. These models were assembled in simulation code, ABAQUS v6.17 (Dassault Systemes, Vélizy-Villacoublay, France), and physiological loads and boundary conditions were applied to mimic a realistic situation, as shown in Figure 1a. A parametric study of the influence of hip implant design on the stiffness of stems, stresses in implants, and fatigue life in terms of the factor of safety was investigated in a previous study [12]. In the previous study [12], the layer thickness and porosity were changed to investigate the factor of safety. The factor of safety was calculated using the Soderberg approach ( σ a S e   +   σ m S y = 1 n ), where, σa is stress amplitude, σs is mean stress, Se is endurance limit, Sy is the yield strength of the material and n is the factor of safety. The factor of safety is calculated to determine the safety of the structure under a certain load [31]. For instance, if the factor of safety is greater than 1 for a certain load, then the structure is considered safe under that load; otherwise, the structure is considered unsafe. In the Soderberg approach, the material properties and the results of simulations were used to calculate the factor of safety. However, these finite element simulations and post-processing calculations are time-consuming and computationally expansive compared with machine learning approaches. In this study, values of inputs of stem designs (layer thickness and porosity) and results (stresses, stiffness, and factor of safety), obtained from Mehboob et al. [12], were used to train the machine learning algorithms and a new test dataset was created based on layer thickness and porosities to predict the stiffness of stems, stresses in implants, and the factor of safety based on previously trained values.
In this study, the data was used to feed and train six types of algorithms. These algorithms were validated with the published data and errors were calculated. After training these algorithms, new design parameters of the hip stem were created with different thicknesses of outer dense layers (2.5 and 3 mm) and inner porosities (10–80%) as shown in Figure 1b. The validated algorithms were used to predict the stiffness of stems, stresses in outer dense layers and porous sections, and factor of safety of newly designed implants.

2.2. Machine Learning

The dataset used in this research is based on finite element analysis as mentioned in Section 2.1. Machine learning algorithms are used to predict the outputs of the simulations for various design parameter settings. The input variables are dense layer thickness (DLT) in mm and porosity of porous section (PPS) in percent: v1 and v2, respectively. The output variables, v3 to v6, are stiffness of stems (SS) in N/mm, maximum stresses in dense layer (MSDL) in MPa, maximum stresses in porous section (MSPS) in MPa, and factor of safety (FS). Table 1 shows the dataset of dimension 30 × 6.
The dataset requires minor preprocessing before machine learning algorithms are utilized for the purposes of predicting the simulated results. A MinMax normalization was applied to obtain a normalized dataset in the range [0.1, 1]. It is noted that a normalization of range [0, 1] is problematic due to the upcoming calculations of prediction accuracies. The MinMax normalization is applied using the following equation.
n o r m a l i z e d   i n p u t = v i j min ( v i j ) max ( v i j ) min ( v i j )
where i is the input variable number and j is the simulation record. No data imputation techniques were used due to the lack of missing values in the dataset. The entire dataset is then vertically split into x = [v1, v2] and y = [v3, v4, v5, v6]. A shuffle split is then horizontally implemented to divide the dataset into training which accounts for 70% and validation which accounts for the remaining 30% of the dataset. The resulting training and validation matrices are of dimensions 21 × 6 and 9 × 6, respectively. The vertical split is shuffled using a NumPy seed [32] to allow for a non-sequential sampling with replacement which is next implemented 10 times for each of the utilized machine learning algorithms.
As mentioned earlier, machine learning techniques are generally categorized into classification and prediction algorithms. The machine learning algorithms used in this paper are generally prediction based as the output to be predicted is numerical. Six different algorithms are used, namely decision tree regression (DTR), linear regression, ridge regression (RR), lasso regression (LSR), elastic net (EN), and multilayer perceptron (MLP) regression. Each of the previously mentioned techniques is implemented using its default hyperparameter settings, as shown in Table 2, as per the Scikit-learn library [33]. A brief about each of the machine learning algorithms is presented next.

2.2.1. Decision Tree Regression (DTR)

Decision trees can be used for both classification and regression tasks [34]. More specifically, decision trees can be referred to as classification trees and regression trees depending on the studied task. A decision tree regressor or a regression tree is used in this paper in order to predict numerical values. A decision tree is considered to be a collection of splits based on threshold values at the training set level. The information obtained from a trained decision tree is then applied to validation and test sets for the purposes of prediction. A decision tree contains root nodes, decision nodes for splitting, and leaf nodes where the final results are shown. A simple demonstration of the decision tree structure, inspired by Bulbul et al. [34], is presented in Figure 2. In Figure 2, decision nodes are where specific variable values are decided upon, to be assigned to one of the following two leaves based on a criterion established in the decision node. This criterion can be a threshold value for numerical variables or a voting system in categorical variables. Variables are ordered based on which variable is a better splitter of the data to produce more accurate predictions. It should be noted that decision nodes in Figure 2 are simultaneously used for the prediction of all four leaf nodes but can be subdivided as depicted in the figure.

2.2.2. Linear Regression (LR)

One of the most commonly used machine learning predictors is linear regression. It has been shown to aid in the prediction of output variables based on input variables by means of best-fit linear relationship navigation [35] which is done using least squares minimization. The following equation illustrates the mechanism of the linear regression predictor.
Y = β 0 + n = 1 N β n X n + ε n
where Y is the output vector, X n are the multiple input variables, β 0 is a constant, β n is the estimated linear parameter signified by the slope which is the regression coefficient, and ε n is the error. According to Ogutu et al. [36], and for the purposes of illustration, a basic linear regression model can be written as in the following equation.
y = µ 1 n + X β + e
where β is a vector of coefficients and e is the residual error vector.

2.2.3. Ridge Regression (RR)

Ridge regression is an extension of linear regression [36] where 2 regularization is used as shown in the following equation.
β ^ ( r i d g e ) = arg min β y X β 2 2 + λ β 2 2  
where y X β 2 2 = i = 1 n ( y i x i T β ) 2 is the loss function with 2 regularization.

2.2.4. Lasso Regression (LSR)

Another extension of linear regression is lasso regression [36] where 1 regularization is utilized as presented next.
β ^ ( l a s s o ) = arg min β y X β 2 2 + λ β 1
where β 1 = j = 1 p | β j | is the loss function with 1 regularization.

2.2.5. Elastic Nets (EN)

The third and final extension of linear regression that is utilized in this paper is elastic nets. It utilizes both 1 and 2 regularization [36]. It is useful for high-dimensional data, which is not the case in this paper. Nonetheless, EN was applied for the purposes of comparison. Elastic nets can be mathematically described as shown in the following equation.
β ^ ( E N ) = ( 1 + λ 2 n ) { arg min β y X β 2 2 + λ 2 β 2 2 + λ 1 β 1 }

2.2.6. Multilayer Perceptron (MLP)

An artificial neural network, also called a multilayer perceptron, can be utilized in its most basic form for the purposes of supervised learning. It nonlinearly maps inputs to outputs by utilizing nodes and their associated weights [37]. The weights connect the nodes to produce the output which is considered to be the sum of inputs. This mapping is implemented using an activation function which is the rectified linear unit (Relu) in this paper. More specifically, weights are updated through a backpropagation process performed by the algorithm to further refine the predicted output. This process is known as training the neural network. In this research, the weights are applied to the dense layer thickness and porosity of porous section as inputs, leading to stiffness of stems, stresses in dense layer and porous sections, and factor of safety as outputs. Neural networks can be and already are being used in almost all fields of scientific studies. Figure 3 shows a basic fully connected neural network schematic where w i represent the weights, i l is the input layer, h l i are the hidden layers, and o l is the output layer.
Each of the machine learning predictors is used to fit the data 10 times according to the NumPy seed in the range of [0, 9] which creates a sampling with a replacement scheme. The algorithms are trained on the 70% training portion and validated on the 30% validation portion of the dataset. This approach creates 60 different runs that were implemented on Python 3 [38]. The validation root mean squared error (RMSE) and mean absolute percentage error (MAPE) are reported in the results and discussion section for each of the 60 runs. RMSE and MAPE are calculated using Eqs. 7 and 8, respectively.
R M S E = 1 n t = 1 n ( A t P t ) 2
M A P E = ( 1 n t = 1 n | A t P t | A t ) 100 %
where A is the actual value and P is the predicted value. RMSE and MAPE results are reported in the results and discussion section where a comprehensive comparison is shown across the six utilized machine learning algorithms and their associated samples. An average, standard deviation, maximum, and minimum result from Equations (7) and (8) is also shown for the purposes of clearly determining the lowest validation prediction error.
Additionally, the trained algorithms are tested against an unseen test set. The test set inputs are shown in Table 3 where different design parameters are chosen to truly test the algorithms’ feasibility in accurately following the trend of the training and validation sets’ output variables.
After using the trained and validated models to predict the outputs of the test set, the trend is observed, and the best algorithm is chosen based on a voting system. The voting system utilizes the proportional trend validity for each algorithm and produces a score out of 100. The best algorithm is then used to showcase its trend results for each of the four output variables. The criteria for choosing the best algorithm depend on two factors. The first factor is the monotonicity of the prediction where no two consecutive predicted outputs are repeated. To clarify, monotonicity means always increasing or always decreasing. The second factor is the nonexistence of negative value predictions as some negative values are wrongly predicted by some algorithms where no negative values exist in the original data.

3. Results and Discussion

3.1. Finite Element Analysis

This study investigates the biomechanical performance of semi-porous hip stems to address issues such as stress shielding, which mainly occurs in Gruen Zone 7 due to a solid dense metallic stem. Designing porous and semi-porous implants can address the complications of stiff dense implants. However, in silico, in vitro, and in vivo investigations of the newly designed porous and semi-porous complex structure are time-consuming and costly; therefore, the results of computer simulations were validated using machine learning in this study. After validation, the algorithms were trained and employed to predict the biomechanical performance of newly designed semi-porous hip stems. Finite element analysis of implants showed that by increasing the thickness of the outer layer, the stresses in the layer were decreased, as discussed in a previous study [12]. In addition, by increasing the porosity, the stresses were decreased. These thicknesses of the layer and the change in porosities affected the stiffness of the stems [12]. Similar findings were observed when a stochastic open-cell porous structure was incorporated in the porous femoral stem [39]. The results of the study were validated experimentally, and showed a reduction of 31% in flexural stiffness when a 33% porous stem was simulated compared with the dense stem. In another study [40], a cobalt-chromium implant with a porous architecture was introduced with various porosities. The results showed that the stiffness was decreased by increasing the porosity, which is in agreement with the current study. Hazlehurst et al. [41] also modeled porous and dense stems simulated in finite element analysis. The results also showed that the porous stem takes more stress as compared with the dense stem. Further validating the finite element results in [12], the new designs in this study also predicted similar trends. Increasing the layer thickness and decreasing the porosity increased the stiffness of the stem. Even while keeping the same porosity, increasing the layer thickness caused an increase in implant stiffness which was predicted by the algorithms. Similarly, keeping the same layer thickness and decreasing the porosity also caused an increase in stiffness. Moreover, increasing the thickness of the outer layer and increasing the porosity improved the factor of safety, which is also consistent with the previous study [7]. These increases in stiffness and factor of safety are good in the view of longevity of any part; however, much higher stiffness may cause stress shielding and reduce bone density around the hip stem. Therefore, an optimal design is required to satisfy both simultaneously; this will be included in future research.

3.2. Machine Learning Predictions

Following the application of the machine learning methodology in Section 2.2, the results are showcased and discussed in this section. Table 4 shows the validation RMSE and MAPE results across all utilized algorithms and sample seeds.
The lowest average MAPE in terms of percentage is the decision tree regression algorithm. However, since observing a reasonable trend is the desired outcome, DTR is excluded from being the algorithm of choice. In fact, all algorithms, after training and validation, were used to predict the test set, and the trend was observed. Table 5 shows the scoring of trend validity on the test set using the trained and validated machine learning algorithms.
It is evident from Table 5 that the two algorithms with the best trend score are ridge and linear regression. Linear regression was excluded due to its prediction of negative values according to the second selection criterion mentioned in the methodology section. To reiterate, the criteria for choosing the best algorithm depended on two factors. The first factor is the monotonicity of the prediction, where no two consecutive predicted outputs are repeated. The second factor is the nonexistence of negative value predictions. In Table 5, a score of 0 or 1 is used to describe criteria satisfaction, where a score of 1 indicates that the criterion of monotonicity is met, whereas a score of 0 indicates that monotonicity is not met. Figure 4 shows the validation prediction of SS, MSDL, MSPS, and FS, using the ridge regression algorithm on the fourth sample seed.
As shown in Figure 4, ridge regression is able to reasonably predict the validation set outputs. However, ridge regression is not the best algorithm in terms of average validation accuracy. On the one hand, it is noted that the ridge regression algorithm under-predicts the values of the stiffness of stems and the factor of safety. On the other hand, the predictions of maximum stresses in the dense layer and maximum stresses in the porous section are each considered to be a mix of both over-prediction and under-prediction. For the purposes of comparison, and since the fourth sample seed was used to produce Figure 4, DTR is used to produce Figure 5 as follows.
Figure 5 is illustrated to showcase the validation prediction of DTR which is the highest-performing algorithm in terms of average validation MAPE. Nonetheless, the trend of test set predictions as a result of using DTR does not satisfy the criteria established in this paper. As shown in Figure 5, decision tree regression is able to more reasonably predict the validation set outputs as compared with ridge regression. In fact, it is the best-performing machine learning algorithm based on the average MAPE results. Decision tree regression mostly under-predicts the true values of stiffness of stems and the factor of safety. The algorithm of DTR strictly over-predicts maximum stresses in the dense layer while producing a mixture of over- and under-predictions for the maximum stresses in the porous section.
The following figure, Figure 6, shows the test set trend of predicted results based on input design variables of the hip stem using the ridge algorithm as trained on the fourth sample seed. The predicted results show that the stiffness of the stem, maximum stresses in the porous section, and factor of safety were increased when porosity was decreased while keeping the same outer layer thickness of 2.5 mm. A similar trend was observed when the outer layer thickness was 3 mm and porosity was decreased. Contrarily, the stresses in porous sections were decreased by decreasing the porosity of stems in both outer layer thicknesses of 2.5 and 3 mm. The trend of these predicted results agreed with the validated results of finite element simulations. The increase in stiffness of the stem by decreasing the porosity is logical because the amount of material in the stem was increased, which resists more against loads. Similarly, the stresses in the dense layer were increased by decreasing the porosity, which showed that the dense section started taking more load as the material of the porous section was decreased. In addition, the factor of safety usually increases when the volume of material is increased, and this was the case in the prediction. However, the stresses in the porous section were decreased by decreasing the porosity because the porous material was not able to take greater loads and yielded.
It is clear that the trends of the predicted test set, as shown in Figure 6, are logical and in agreement with the original dataset. For this reason, it can be concluded that ridge regression can be used for the purposes of future design considerations without the need for further simulation experiments. The necessity to conduct simulation experiments is time-consuming and may be completely replaced by machine learning prediction efforts when the algorithms are fully optimized and proven to perfectly perform predictions.
It is evident that the trend produced by the ridge regression predictions of the test set followed the original dataset’s trend in terms of the output variables v3 to v6, stiffness of stems, maximum stresses in a dense layer, maximum stresses in a porous section, and factor of safety. For the purposes of comparison, DTR was also tested on the same test set. As can be seen in Figure 7, the trend of the predicted test set using DTR is a poor representation of the original dataset’s trend in terms of the output variables v3 to v6. Figure 7 shows that almost every two successive predicted outputs have the same level, which does not show realistic results, as shown in Figure 6 for ridge regression. When the outer layer thickness or porosity is changed, the stiffness of the stem and stresses in the dense and porous sections should be changed, which is not accurately predicted by DTR.
Figure 8 shows the calculation times of FEA calculations, ML validation, testing, and total ML time. Typically, it takes hours to a few days to perform finite element calculations through modeling and finite element analysis, whereas machine learning models take significantly less time. Even if machine learning training time is also taken into account, the total time required is still significantly lower than that taken for FEA calculations. Hence, machine learning models are computationally more efficient and significantly accelerate the predictions of the biomechanical performance of new designs for hip stems. In addition to time reduction, it is worth noting that these results were also obtained using a relatively small data set.
The framework of this study has great potential to aid new stem designs, and thus may be applicable for clinical use. It allows for the rapid exploration of the biomechanical performance of various designs of hip stems to be used and applied to the personal conditions of patients. Thus, this study has reduced the computational cost and time for the design process of stems. In general, machine learning models show robust performance in predicting the biomechanical performance of different designs of hip stems; however, the prediction efficiency can always be improved by investigating a larger dataset.

4. Conclusions

Finite element simulations are computationally expensive as compared with machine learning algorithms, but give more accurate results. The validation of finite element results was successfully achieved using machine learning algorithms and this validation of finite element results would be far too expensive in in vivo experimentations.
Using the simulated data, a variety of machine learning algorithms were utilized for the purposes of predicting the stiffness of stems, maximum stresses in a dense layer, maximum stresses in a porous section, and factor of safety. The ridge regression algorithm was shown to produce the most accurate test set prediction in terms of trend as compared with the original dataset’s trend. Thus, the trend of biomechanical performance of new stem designs was successfully predicted using trained machine learning algorithms.
Future research potential depends on further optimizing machine learning algorithms to achieve a better prediction scheme for the biomechanical performance of different stem designs. An example of further optimization could depend on using grid search to optimize the selected hyperparameter values. Another avenue of future research could include utilizing different test sets to facilitate machine-learning-informed implant design without the need for running computer simulations. In addition, future research can also focus on using machine learning techniques in the designing of high-performing hip implants rather than simply predicting design parameters.
One of the limitations of this paper is the relatively small dataset that was used to train the machine learning models. Accordingly, the accuracy of the predictive models was affected. The use of a larger dataset may improve the accuracy of the prediction. Another limitation is that these results are predicted based on finite element analysis and machine learning approaches; however, the biomechanical performance should be investigated using in vivo experiments. Another limitation is the fact that ML model complexity was not considered when evaluating the models used in this research. More complex models are usually expected to perform better but they are less interpretable. The current algorithms are not trained to predict biochemical effects, including biocompatibility, stress shielding, and bone remodeling, which are gradual changes in the peri-stem bone changes and are important in assessing the secondary stability of implants. Hip stems can be optimized considering changes in bone density to maximize the life of total hip arthroplasty. Once the implant design is optimized, a biochemical study should be carried out to study the effect of material and design on biological response.

Author Contributions

Conceptualization, K.A., H.M., R.A. and F.T.; methodology, K.A. and H.M.; software, K.A., H.M. and F.T.; validation, K.A., H.M., R.A. and F.T.; formal analysis, K.A., H.M. and R.A.; investigation, K.A., H.M., R.A. and F.T.; data curation, K.A. and H.M.; writing—original draft preparation, K.A., H.M. and R.A.; writing—review and editing, K.A., H.M., R.A. and F.T.; visualization, K.A. and H.M.; Discussion, K.A., H.M., R.A. and F.T.; Conclusion, K.A., H.M., R.A. and F.T.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Prince Sultan University, for paying the Article Processing Charges (APC) of this publication. This paper was made possible by NPRP grant# NPRP 8-876-2-375 from the Qatar National Research Fund (a member of Qatar Foundation).

Data Availability Statement

Data is available in Table 1 and Table 3.

Acknowledgments

The authors would like to thank Prince Sultan University for their support. This paper was made possible by NPRP grant# NPRP 8-876-2-375 from the Qatar National Research Fund (a member of Qatar Foundation). The findings achieved herein are solely the responsibility of the authors.

Conflicts of Interest

The authors have declared no conflict of interest.

References

  1. Mi, Z.; Shuib, S.; Hassan, A.; Shorki, A.; Ibrahim, M. Problem of stress shielding and improvement to the hip Implat designs: A review. J. Med. Sci. 2007, 7, 460–467. [Google Scholar]
  2. Sumner, D. Long-term implant fixation and stress-shielding in total hip replacement. J. Biomech. 2015, 48, 797–800. [Google Scholar] [CrossRef]
  3. Savio, D.; Bagno, A. When the total hip replacement fails: A review on the stress-shielding effect. Processes 2022, 10, 612. [Google Scholar] [CrossRef]
  4. Thakur, A.; Kaya, S.; Abousalem, A.; Sharma, S.; Ganjoo, R.; Assad, H.; Kumar, A. Computational and experimental studies on the corrosion inhibition performance of an aerial extract of Cnicus Benedictus weed on the acidic corrosion of mild steel. Process Saf. Environ. Prot. 2022, 161, 801–818. [Google Scholar] [CrossRef]
  5. Arabnejad, S.; Johnston, B.; Tanzer, M.; Pasini, D. Fully porous 3D printed titanium femoral stem to reduce stress-shielding following total hip arthroplasty. J. Orthop. Res. 2017, 35, 1774–1783. [Google Scholar] [CrossRef] [PubMed]
  6. Cortis, G.; Mileti, I.; Nalli, F.; Palermo, E.; Cortese, L. Additive manufacturing structural redesign of hip prostheses for stress-shielding reduction and improved functionality and safety. Mech. Mater. 2022, 165, 104173. [Google Scholar] [CrossRef]
  7. Head, W.C.; Bauk, D.J.; Emerson, R.H., Jr. Titanium as the material of choice for cementless femoral components in total hip arthroplasty. Clin. Orthop. Relat. Res. 1995, 311, 85–90. [Google Scholar]
  8. Sidambe, A.T. Biocompatibility of advanced manufactured titanium implants—A review. Materials 2014, 7, 8168–8188. [Google Scholar] [CrossRef]
  9. Zheng, Q.; Mao, L.; Shi, Y.; Fu, W.; Hu, Y. Biocompatibility of Ti-6Al-4V titanium alloy implants with laser microgrooved surfaces. Mater. Technol. 2022, 37, 2039–2048. [Google Scholar] [CrossRef]
  10. Alkhatib, S.E.; Tarlochan, F.; Mehboob, H.; Singh, R.; Kadirgama, K.; Harun, W.S.B.W. Finite element study of functionally graded porous femoral stems incorporating body-centered cubic structure. Artif. Organs 2019, 43, E152–E164. [Google Scholar] [CrossRef]
  11. Mehboob, H.; Tarlochan, F.; Mehboob, A.; Chang, S.-H. Finite element modelling and characterization of 3D cellular microstructures for the design of a cementless biomimetic porous hip stem. Mater. Des. 2018, 149, 101–112. [Google Scholar] [CrossRef]
  12. Mehboob, H.; Tarlochan, F.; Mehboob, A.; Chang, S.-H.; Ramesh, S.; Harun, W.S.W.; Kadirgama, K. A novel design, analysis and 3D printing of Ti-6Al-4V alloy bio-inspired porous femoral stem. J. Mater. Sci. Mater. Med. 2020, 31, 78. [Google Scholar] [CrossRef]
  13. Du Plessis, A.; Razavi, S.M.J.; Benedetti, M.; Murchio, S.; Leary, M.; Watson, M.; Bhate, D.; Berto, F. Properties and applications of additively manufactured metallic cellular materials: A review. Prog. Mater. Sci. 2021, 125, 100918. [Google Scholar] [CrossRef]
  14. Alkhatib, S.E.; Mehboob, H.; Tarlochan, F. Finite element analysis of porous titanium alloy hip stem to evaluate the biomechanical performance during walking and stair climbing. J. Bionic Eng. 2019, 16, 1103–1115. [Google Scholar] [CrossRef]
  15. Tarlochan, F.; Mehboob, H.; Mehboob, A.; Chang, S.-H. Biomechanical design of a composite femoral prosthesis to investigate the effects of stiffness, coating length, and interference press fit. Compos. Struct. 2018, 204, 803–813. [Google Scholar] [CrossRef]
  16. Borjali, A.; Chen, A.F.; Muratoglu, O.K.; Morid, M.A.; Varadarajan, K.M. Detecting total hip replacement prosthesis design on plain radiographs using deep convolutional neural network. J. Orthop. Res. 2020, 38, 1465–1471. [Google Scholar] [CrossRef]
  17. Hambli, R.; Katerchi, H.; Benhamou, C.-L. Multiscale methodology for bone remodelling simulation using coupled finite element and neural network computation. Biomech. Model. Mechanobiol. 2011, 10, 133–145. [Google Scholar] [CrossRef] [PubMed]
  18. Li, Z. Predicting bone regeneration from machine learning. Nat. Comput. Sci. 2021, 1, 509–510. [Google Scholar] [CrossRef]
  19. Murphy, K.P. Machine Learning: A Probabilistic Perspective; MIT Press: Cambridge, MA, USA, 2012. [Google Scholar]
  20. Cilla, M.; Borgiani, E.; Martínez, J.; Duda, G.N.; Checa, S. Machine learning techniques for the optimization of joint replacements: Application to a short-stem hip implant. PLoS ONE 2017, 12, e0183755. [Google Scholar] [CrossRef]
  21. Bertolini, M.; Mezzogori, D.; Neroni, M.; Zammori, F. Machine Learning for industrial applications: A comprehensive literature review. Expert Syst. Appl. 2021, 175, 114820. [Google Scholar] [CrossRef]
  22. Kanakaris, N.; Karacapilidis, N.; Kournetas, G.; Lazanas, A. Combining machine learning and operations research methods to advance the project management practice. In International Conference on Operations Research and Enterprise Systems; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  23. Ashfaq, T.; Khalid, R.; Yahaya, A.S.; Aslam, S.; Azar, A.T.; Alsafari, S.; Hameed, I.A. A Machine Learning and Blockchain Based Efficient Fraud Detection Mechanism. Sensors 2022, 22, 7162. [Google Scholar] [CrossRef]
  24. Henrique, B.M.; Sobreiro, V.A.; Kimura, H. Literature review: Machine learning techniques applied to financial market prediction. Expert Syst. Appl. 2019, 124, 226–251. [Google Scholar] [CrossRef]
  25. Sharma, N.; Thakur, M.S.; Sihag, P.; Malik, M.A.; Kumar, R.; Abbas, M.; Saleel, C.A. Machine learning techniques for evaluating concrete strength with waste marble powder. Materials 2022, 15, 5811. [Google Scholar] [CrossRef]
  26. Upadhya, A.; Thakur, M.S.; Al Ansari, M.S.; Malik, M.A.; Alahmadi, A.A.; Alwetaishi, M.; Alzaed, A.N. Marshall Stability Prediction with Glass and Carbon Fiber Modified Asphalt Mix Using Machine Learning Techniques. Materials 2022, 15, 8944. [Google Scholar] [CrossRef] [PubMed]
  27. Khan, M.A.; Ahmad, S.; El-Affendi, M.A.; Zaka, R.; Mahmood, S.; Jehangir, M. Predicting Mental Health of Best Human Capital for Sustainable Organization through Psychological and Personality Health Issues: Shift from Traditional to Novel Machine Learning-Supervised Technique Approach. Biomed Res. Int. 2022, 2022, 5775640. [Google Scholar] [CrossRef] [PubMed]
  28. Villamor, E.; Monserrat, C.; Del Río, L.; Romero-Martín, J.; Rupérez, M.J. Prediction of osteoporotic hip fracture in postmenopausal women through patient-specific FE analyses and machine learning. Comput. Methods Programs Biomed. 2020, 193, 105484. [Google Scholar] [CrossRef] [PubMed]
  29. Alastruey-López, D.; Ezquerra, L.; Seral, B.; Pérez, M. Using artificial neural networks to predict impingement and dislocation in total hip arthroplasty. Comput. Methods Biomech. Biomed. Eng. 2020, 23, 649–657. [Google Scholar] [CrossRef]
  30. Jun, Z.; Youqiang, Z.; Wei, C.; Fu, C. Research on prediction of contact stress of acetabular lining based on principal component analysis and support vector regression. Biotechnol. Biotechnol. Equip. 2021, 35, 462–468. [Google Scholar] [CrossRef]
  31. Al Zoubi, N.F.; Tarlochan, F.; Mehboob, H. Mechanical and Fatigue Behavior of Cellular Structure Ti-6Al-4V Alloy Femoral Stems: A Finite Element Analysis. Appl. Sci. 2022, 12, 4197. [Google Scholar] [CrossRef]
  32. Harris, C.R.; Millman, K.J.; Van Der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J. Array programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef]
  33. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  34. Bulbul, A.M.R.; Khan, K.; Nafees, A.; Amin, M.N.; Ahmad, W.; Usman, M.; Nazar, S.; Arab, A.M.A. In-Depth Analysis of Cement-Based Material Incorporating Metakaolin Using Individual and Ensemble Machine Learning Approaches. Materials 2022, 15, 7764. [Google Scholar] [CrossRef]
  35. Shaban, M.; Alsharekh, M.F.; Alsunaydih, F.N.; Alateyah, A.I.; Alawad, M.O.; BaQais, A.; Kamel, M.; Nassef, A.; El-Hadek, M.A.; El-Garaihy, W.H. Investigation of the Effect of ECAP Parameters on Hardness, Tensile Properties, Impact Toughness, and Electrical the Conductivity of Pure Cu through Machine Learning Predictive Models. Materials 2022, 15, 9032. [Google Scholar] [CrossRef] [PubMed]
  36. Ogutu, J.O.; Schulz-Streeck, T.; Piepho, H.-P. Genomic Selection Using Regularized Linear Regression Models: Ridge Regression, Lasso, Elastic Net and Their Extensions. In BMC Proceedings; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  37. Gardner, M.W.; Dorling, S. Artificial neural networks (the multilayer perceptron)—A review of applications in the atmospheric sciences. Atmos. Environ. 1998, 32, 2627–2636. [Google Scholar] [CrossRef]
  38. Guido, V.R.; Drake, F., Jr. Python 3 Reference Manual; CreateSpace: Scotts Valley, CA, USA, 2009. [Google Scholar]
  39. Simoneau, C.; Terriault, P.; Jetté, B.; Dumas, M.; Brailovski, V. Development of a porous metallic femoral stem: Design, manufacturing, simulation and mechanical testing. Mater. Des. 2017, 114, 546–556. [Google Scholar] [CrossRef]
  40. Limmahakhun, S.; Oloyede, A.; Sitthiseripratip, K.; Xiao, Y.; Yan, C. Stiffness and strength tailoring of cobalt chromium graded cellular structures for stress-shielding reduction. Mater. Des. 2017, 114, 633–641. [Google Scholar] [CrossRef]
  41. Hazlehurst, K.B.; Wang, C.J.; Stanford, M. A numerical investigation into the influence of the properties of cobalt chrome cellular structures on the load transfer to the periprosthetic femur following total hip arthroplasty. Med. Eng. Phys. 2014, 36, 458–466. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Finite element models of semi–porous stems; (a) finite element models used in the previous study [12], (b) new design parameters for prediction of biomechanical performance using machine learning.
Figure 1. Finite element models of semi–porous stems; (a) finite element models used in the previous study [12], (b) new design parameters for prediction of biomechanical performance using machine learning.
Jfb 14 00156 g001
Figure 2. Decision tree structure showing root, decision, and leaf nodes as it relates to the implemented model.
Figure 2. Decision tree structure showing root, decision, and leaf nodes as it relates to the implemented model.
Jfb 14 00156 g002
Figure 3. MLP basic structure shows an input layer, two hidden layers, and an output layer.
Figure 3. MLP basic structure shows an input layer, two hidden layers, and an output layer.
Jfb 14 00156 g003
Figure 4. Ridge regression validation set predictions of (a) stiffness of stems; (b) maximum stresses in dense layer; (c) maximum stresses in porous section; and (d) factor of safety.
Figure 4. Ridge regression validation set predictions of (a) stiffness of stems; (b) maximum stresses in dense layer; (c) maximum stresses in porous section; and (d) factor of safety.
Jfb 14 00156 g004
Figure 5. Decision tree regression validation set predictions of (a) stiffness of stems; (b) maximum stresses in dense layer; (c) maximum stresses in porous section; and (d) factor of safety.
Figure 5. Decision tree regression validation set predictions of (a) stiffness of stems; (b) maximum stresses in dense layer; (c) maximum stresses in porous section; and (d) factor of safety.
Jfb 14 00156 g005
Figure 6. Ridge regression test set predictions of all four output variables.
Figure 6. Ridge regression test set predictions of all four output variables.
Jfb 14 00156 g006
Figure 7. Decision tree regression test set predictions of all four output variables.
Figure 7. Decision tree regression test set predictions of all four output variables.
Jfb 14 00156 g007
Figure 8. The computational time for the biomechanical performance of hip stems required by FEA calculations and ML.
Figure 8. The computational time for the biomechanical performance of hip stems required by FEA calculations and ML.
Jfb 14 00156 g008
Table 1. Dataset details.
Table 1. Dataset details.
SimulationDLTPPSSSMSDLMSPSFS
1090134901040.165116842
20772981731360.273558063
30635421892170.311399648
40478441893120.355395013
503012841892960.649761978
601816711892890.956662736
70.590900969310.439016544
80.577999848640.504168444
90.56311417091060.607191898
100.54713215831470.754651623
110.53015934521901.012067069
120.51818453682171.274545641
131901340578150.870252011
141771412540350.906265976
151631508493650.978832319
1614716264451011.09870317
1713018033831481.299523956
1811819673371841.503057065
191.5901639445101.277959118
201.5771691427251.288267624
211.5631757403471.371124031
221.5471835378761.459330144
231.53019493451171.644608986
241.51820553181521.820361118
25290184537981.696200284
262771883368191.6979466
272631928355371.740895669
282471978340621.787807737
292302050321991.943079537
3021821153051332.079133065
Table 2. Selected default hyperparameter values for utilized algorithms.
Table 2. Selected default hyperparameter values for utilized algorithms.
Algorithm HyperparameterHyperparameter Value
Decision tree regression Criterion Squared error
SplitterBest
Linear, ridge, lasso, elastic netAlpha1.0
Fit interceptTrue
Multilayer perceptronActivationRectified linear unit (Relu)
Hidden layer sizes100
SolverAdaptive momentum (Adam)
Table 3. Test set values.
Table 3. Test set values.
#DLTPPS
12.580
22.570
32.560
42.550
52.540
62.530
72.520
82.510
9380
10370
11360
12350
13340
14330
15320
16310
Table 4. RMSE and MAPE comparison.
Table 4. RMSE and MAPE comparison.
DTRLRRRLSRENMLP
SeedRMSEMAPERMSEMAPERMSEMAPERMSEMAPERMSEMAPERMSEMAPE
00.1116.910.1024.220.1325.880.2249.700.2249.700.1432.13
10.0813.870.1020.630.1326.590.2147.980.2147.980.1430.28
20.1128.100.1747.350.2276.060.30106.340.30106.340.2170.74
30.1016.990.0818.720.1227.440.2046.700.2046.700.1020.67
40.0915.290.0712.540.1118.210.2237.530.2237.530.0713.39
50.1116.600.0921.730.1327.340.2352.720.2352.720.1125.69
60.1022.750.1531.140.1952.620.2779.880.2779.880.1952.65
70.0917.100.1335.670.1642.260.2572.020.2572.020.2155.02
80.1126.500.1548.140.1967.980.2899.770.2899.770.1964.38
90.1522.040.1420.110.1934.350.2753.360.2753.360.1930.10
Average0.1019.620.1228.030.1639.870.2464.600.2464.600.1639.50
Sample σ0.024.900.0312.220.0419.580.0323.770.0323.770.0519.62
Max0.1528.100.1748.140.2276.060.30106.340.30106.340.2170.74
Min0.0813.870.0712.540.1118.210.2037.530.2037.530.0713.39
Table 5. Trend scores based on the test set.
Table 5. Trend scores based on the test set.
Trend Scores
SeedDTRLRRRLSRENMLP
0011001
1011001
2011000
3011001
4011001
5011001
6011000
7011001
8011000
9011000
Trend Score01010006
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

Akkad, K.; Mehboob, H.; Alyamani, R.; Tarlochan, F. A Machine-Learning-Based Approach for Predicting Mechanical Performance of Semi-Porous Hip Stems. J. Funct. Biomater. 2023, 14, 156. https://doi.org/10.3390/jfb14030156

AMA Style

Akkad K, Mehboob H, Alyamani R, Tarlochan F. A Machine-Learning-Based Approach for Predicting Mechanical Performance of Semi-Porous Hip Stems. Journal of Functional Biomaterials. 2023; 14(3):156. https://doi.org/10.3390/jfb14030156

Chicago/Turabian Style

Akkad, Khaled, Hassan Mehboob, Rakan Alyamani, and Faris Tarlochan. 2023. "A Machine-Learning-Based Approach for Predicting Mechanical Performance of Semi-Porous Hip Stems" Journal of Functional Biomaterials 14, no. 3: 156. https://doi.org/10.3390/jfb14030156

APA Style

Akkad, K., Mehboob, H., Alyamani, R., & Tarlochan, F. (2023). A Machine-Learning-Based Approach for Predicting Mechanical Performance of Semi-Porous Hip Stems. Journal of Functional Biomaterials, 14(3), 156. https://doi.org/10.3390/jfb14030156

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