Next Article in Journal
Hybrid Zinc Coatings with Improved Corrosion Resistance Based on Chitosan Oligosaccharides
Previous Article in Journal
Rheological Behavior of TRIP600 Steel during Deformation
Previous Article in Special Issue
Aging Treatment Induces the Preferential Crystallographic Orientation of αs in the Near-α Titanium Alloy Ti60
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Machine Learning in Constitutive Relationship Prediction of Porous Titanium Materials for Artificial Bone

School of Mechanical Engineering, Jiangsu University of Science and Technology, Zhenjiang 212000, China
*
Author to whom correspondence should be addressed.
Metals 2024, 14(6), 634; https://doi.org/10.3390/met14060634
Submission received: 26 April 2024 / Revised: 18 May 2024 / Accepted: 22 May 2024 / Published: 27 May 2024
(This article belongs to the Special Issue Light Alloys and Composites)

Abstract

:
Artificial bone porous titanium materials are widely used in orthopedic implants. However, the traditional constitutive model is often limited by the complexity and accuracy of the model, and it is difficult to accurately and efficiently describe the constitutive relationship of porous titanium materials. In this study, structured data were established based on experimental data from published papers, and goodness of fit ( R 2 ), mean absolute error ( M A E ) and mean absolute percentage error ( M A P E ) were used to evaluate the model. The prediction effect of random forest (RF), multi-layer perceptron (MLPR) and support vector machine (SVR) on the constitutive relationship of porous titanium materials was discussed. Through comprehensive comparison, it can be seen that the RF model with max_depth of 24 and n_estimators of 160 has the best performance in prediction, and the average absolute percentage error is less than 4.4%, which means it can accurately predict the temperature sensitivity and strain rate sensitivity of porous titanium materials. And its predictive ability is better than that of the traditional constitutive model, which provides a new idea and method for the constitutive modeling of porous titanium materials.

1. Introduction

In today’s engineering field, the research and application of porous materials have been paid more and more attention. In particular, as an important functional material, porous titanium has excellent mechanical properties and biocompatibility, so it is widely used in orthopedic implants [1,2,3]. The research on the mechanical properties of porous titanium materials is the basis for further research and application of porous titanium. At present, the constitutive model of porous titanium materials is mainly based on the traditional theoretical model [4,5].
In recent years, with the rapid development of machine learning, it has been widely used in various fields, such as Reka using machine learning with machine vision to identify the degree of fruit decay [6], Zheng et al. using machine learning to optimize the installation position of hydraulic cylinders on medical robots [7] and so on. For material prediction, machine learning is also widely used in this field because of its characteristics of learning the laws contained in a wide range of data sets [8,9,10,11]. At present, scholars at home and abroad have launched a certain exploration in this field. In 2023, Lei et al. [12] established a fatigue test data set containing Hastelloy material characteristics, and then built a machine learning framework for fatigue life prediction based on the data set. After comparison, it was found that the prediction effect of the model built based on machine learning was greatly improved compared with traditional prediction methods. Some scholars have tested the prediction effect of common machine learning algorithms such as SVM, RF, AdaBoost and DNN on material prediction and found the most suitable machine learning algorithm corresponding to material properties and the optimal parameters of each algorithm through comprehensive comparison [13,14,15,16,17]. For hyperelastic porous materials, Liu Zhentao et al. [18] built a constitutive model based on artificial neural networks to predict their mechanical behavior. Meanwhile, some scholars have improved the algorithm model for different cases. For example, Zhao [19] and Ismail [20] used the genetic algorithm and particle swarm optimization algorithm, respectively, to optimize the neural network and predict the creep fatigue life of metal oxides and the wear rate of nanocomposites. In 2021, Guo et al. [21] combined convolutional neural networks and an improved long short-term memory neural network to propose a deep learning model that can quickly identify constitutive model parameters and has good denoising performance. Some scholars [22,23] combined intelligent optimization algorithms such as Gray Wolf optimizer and the particle swarm optimization algorithm with support vector machine and XGBoost algorithms to form a hybrid model, and obtained a better model by comparing the performance differences between different hybrid models.
To sum up, machine learning has been used to predict the properties of materials. However, the traditional constitutive model is often based on phenomenology or physics, which is limited by its own explicit mathematical description ability, and it is very difficult to find a constitutive model that can accurately describe the complex mechanical behavior of materials under the combined action of multiple effects [24]. Therefore, when describing some materials with complex properties, it is necessary to optimize and modify the classical constitutive model based on the current material properties, such as introducing dislocation dynamics, grain evolution and crystal plasticity theories [25]. In addition, the prediction range of the traditional constitutive model is limited, and it does not have enough predictive properties to describe the properties of materials under certain conditions of high temperature or high strain rate [4,26]. In contrast, the data set of the constitutive model based on machine learning can be continuously expanded, and the training and verification time of the new model is shorter [27]; so, the flexibility and prediction range of the model are more advantageous than those of the traditional constitutive model. However, the construction of a constitutive model using machine learning relies on a large amount of data to capture the nonlinear, complex and multivariate behavior of materials; so, the application of machine learning in this construction model relies heavily on the training database [26,28]. Moreover, the machine learning algorithm has difficulties in parameter adjustment and convergence, and poor algorithm parameter adjustment will lead to problems such as model convergence failure and overfitting. Although machine learning still has some problems in the field of material prediction, it has great prospects for development. At present, the research on porous materials, especially porous titanium materials, is still limited, and the adaptability and prediction accuracy of different models for porous titanium materials have not been systematically studied.
In this study, a structured data set containing four characteristic parameters was established through missing value processing, outlier processing and standardization processing, and the prediction effect of RF, MLPR and SVR algorithms on the constitutive relationship of porous titanium materials was discussed. Through comprehensive comparison, the RF model was found to be the best. Then, the RF model was compared with the traditional constitutive model, and the prediction accuracy of the RF model was found to be higher. The purpose of this study was to provide new ideas and methods for the establishment of the constitutive model of porous titanium materials, promote the development of porous materials, and promote the research of related fields to make new breakthroughs and more progress.

2. Materials and Methods

2.1. Model Principle

Three algorithms, random forest, multi-layer perceptron and support vector machine, are used in this study. The random forest model is built by integrating multiple decision trees, which has good robustness and prediction accuracy, but it is easy to overfit the noisy data. Multi-layer perceptrons capture complex relationships in data through nonlinear changes in activation functions in hidden layers. They can learn complex nonlinear relationships and perform well on large-scale data sets. However, they require a lot of data and computational resources to train the model, and the model structure is complex and difficult to interpret. SVM uses the kernel function to capture complex relationships in data. It performs well in high-dimensional space and is suitable for high-dimensional data. However, it has high computational complexity for large-scale data sets and high-dimensional data.

2.1.1. Random Forest Regression (RF)

The random forest algorithm is a supervised learning algorithm. Its main idea is to build a forest composed of decision trees in a random way, and the samples are classified by Information Gain, Gain Ratio and Gini Index at the nodes of each tree. The random forest algorithm achieves high prediction accuracy through random selection of samples, features and integrated prediction.
The implementation process of the random forest algorithm is as follows:
  • N′ (N > N′ > 2N/3) samples are selected each time through N random sampling for the full set D of samples with a capacity of N, so as to form n sample training sets. The samples that are not drawn are divided into test sets.
  • n weak learners are generated based on n sample training sets, but only S′ (S′ << S) are randomly selected from S feature attributes of the sample as feature variables.
  • Therefore, the test set tests n weak learners and obtains n predicted values. After averaging them, the final predicted value can be obtained.

2.1.2. Multi-Layer Perceptron Regression (MLPR)

The multi-layer perceptron algorithm is a supervised learning algorithm. Its main idea is to build an algorithm model by simulating the connection mechanism of the human brain neural network. The multi-layer perceptron algorithm consists of three layers: input layer, hidden layer and output layer. The input layer has the same number of nodes as the feature variable and does not include the activation function. The output layer has the same number of nodes as the output variable; each layer of the activation function in the hidden layer is nonlinearly transformed to capture complex relationships in the input data.
The implementation process of the multi-layer perceptron algorithm is as follows:
  • The weight parameters and threshold parameters of each layer are initialized. The data are passed through the input layer to the first hidden layer. The weighted input is calculated in the hidden layer, the activation function is applied to pass the result to the next layer, and then propagated forward in turn until the result is finally obtained in the output layer.
  • According to the gradient descent method, the connection weight parameters and threshold parameters between neurons are updated to minimize the total loss function.
  • Steps (1) and (2) are repeated until a preset stop condition such as network convergence or the maximum number of iterations is reached.

2.1.3. Support Vector Regression (SVR)

Support vector regression is a supervised learning algorithm. The main idea is to map the feature variables to the high-dimensional feature space by using the kernel function, and then perform linear regression in the high-dimensional feature space and nonlinear regression in the low-dimensional space.
The implementation process of the support vector machine algorithm is as follows [29]:
  • One must estimate the sample data x 1 y 1 ( x N , y N ) , x i , y i R which defines the regression function as follows:
    f ( x ) = w T x + b
    where f(x) is the predicted value, w is the weight vector and b is the bias term.
  • In order to find the minimum w, the relaxation factor ξ i ,   ξ i * is introduced to avoid underfitting the model. The optimization objectives are as follows:
    min 1 2 w 2 + C i = 1 N ( ξ i + ξ i * )
    f ( x i ) y i ε + ξ i y i f ( x i ) ε + ξ i * ξ i , ξ i * 0
    where C is the penalty coefficient. The Lagrange function is introduced, α , α * is introduced for each constraint, and then the kernel function k ( x i , x ) is introduced in the process of solving α , α * .
  • The regression function is as follows [30]:
    f ( x ) = i = 1 N ( α i * α i ) k ( x i , x ) + b

2.2. Model Building

In this work, a constitutive model of porous titanium materials was constructed according to Figure 1. The input characteristics in the model are strain, porosity, temperature and strain rate, and the response variable is stress.

2.2.1. Data Description

The sample data come from the experimental data in a previously published paper [4]. The experimental data are the stress–strain curves of porous titanium materials at different temperatures and strain rates, which are obtained by conducting orthogonal experiments using the split Hopkinson pressure bar (SHPB). The raw material of the experiment is high-purity titanium powder provided by Shijiazhuang Yida Co., Ltd. (Shijiazhuang, China), and the experimental material is made of cylindrical samples with a diameter of 8 mm and a height of 3 mm, 4 mm and 6 mm by using the powder sintering method. The sintering additive uses a 2 wt% polyvinyl alcohol aqueous solution. Considering the high chemical activity of titanium, it can react with oxygen, hydrogen, nitrogen and other elements, so the vacuum sintering method is used, and the vacuum degree is 10−4 Pa. The sintering temperature of titanium powder with a particle size less than or equal to 27 μ m is 1200 °C, the holding time is 2 h and natural cooling is achieved. The sintering temperature of titanium powder with a particle size less than or equal to 74 μm is 1100 °C, the holding time is 2 h and natural cooling is achieved. In the split Hopkinson pressure bar experiment, the sample size was adjusted according to different loading strain rates. A pressure bar with a diameter of 15 mm was used in the experiment, and the experimental conditions of different strain rates were achieved by adjusting the speed of the impact bar. The experiment was repeated 2 to 3 times for each group, and the average value of the effective data of the experimental results was taken to reduce the experimental error. The specific parameters of porous titanium materials are shown in Table 1, which shows the percentage of each element of the experimental material in the overall mass, particle size, pore size and porosity, and the sample data are shown in Table 2. A JSM-6480 scanning electron microscope (Nihon Technology Co., LTD., Osaka, Japan) was used to observe the surface morphology of the sample, and the results are shown in Figure 2.
Pearson correlation coefficient is used to view the degree of linear correlation between features. Its value ranges from −1 to +1. The larger the value, the stronger the correlation. Its expression [30] is as follows:
p = cov ( X , Y ) σ X σ Y = E ( X μ X ) ( Y μ Y ) σ X σ Y
where cov ( X , Y ) is the covariance of the data set X , Y ; μ X , μ Y are the average value of for X , Y ; σ X , σ Y are the deviation of X , Y , respectively. It can be seen from Figure 3 that most of the features are nonlinear. However, there is a certain positive correlation between stress and strain, which is caused by the proportional increase in stress in the plastic deformation stage and the nonlinear increase in the elastic stage. The porosity will reduce the mechanical strength of the material, so there is a weak negative correlation between porosity and stress.

2.2.2. Data Preprocessing

In order to improve the final prediction quality of the data quality assurance model, it is necessary to carry out data cleaning on the sample set, such as missing value processing, outlier processing and standardization processing. Among them, standardization processing can eliminate the dimension gap between the data and accelerate the convergence rate of the model. Here, the Z-score standardization method is used, and its expression [30] is as follows:
X = X μ σ
where X is the value after data standardization; X is the value before data standardization; μ is the average value of the sample; σ is the standard deviation of the sample. Because the activation function in the MLPR model has scale requirements on the data, if the scale difference between features is large, it will affect the convergence speed of the gradient descent method; so, the data used to construct the MLPR model need to be standardized. Because the calculation of the kernel function in the SVR model depends on the distance between features, the calculation accuracy of the kernel function will be affected if the feature scale is different, so the data used to construct the SVR model needs to be standardized. The prediction result of the RF model is mainly related to the distribution of variables and the conditional probability between variables, and has nothing to do with the size of the value itself; so, the data used to build RF model do not need to be standardized.

2.2.3. Data Set Partitioning and Model Evaluation

On the basis of data preprocessing, the sample set is divided into a training set and a test set. The training set is used to train the model, and the test set is used to test the model’s performance. In general, the more data used to train the model, the higher the prediction accuracy of the model, but excessive training will lead to a decline in the overfitting generalization ability of the model; so, the ratio of the training set to the test set is 8:2.
Goodness of fit R 2 , mean absolute error M A E and mean absolute percentage error M A P E were used to evaluate the performance of regression models. The best fit represents the degree to which the regression curve fits the experimental value, and the closer the goodness of fit is to 1, the better the regression curve fits the predicted value. The mean absolute error represents the mean of the absolute error between the predicted value and the experimental value. The calculation formula [31] is as follows:
R 2 = 1 i = 1 n ( Y i Y i ) 2 i = 1 n ( Y i Y ¯ ) 2
M A E = 1 n i = 1 n ( Y i Y i )
M A P E = 100 % n i = 1 n Y ^ i Y i Y i
where Y i represents the predicted value of the model, Y i represents the experimental value, Y ¯ represents the average of the experimental value, and n represents the number of samples.

2.2.4. Model Hyperparameter Tuning

Random Forest Model (RF)

The number of weak learners (n_estimators) and the maximum depth (max_depth) are tuned in the RF model. In this study, the number of weak learners ranged from 10 to 235, with steps of 15, and the maximum depth ranged from 10 to 38, with steps of 2, respectively; R 2 and M A E were used to evaluate the influence of parameter values on the predicted results. According to Figure 4a,b, it is found that R 2 increases with the increase in the number of weak learners and the maximum depth, and M A E decreases with the increase in the number of weak learners and the maximum depth. However, after the number of weak learners reaches 70, M A E fluctuation greatly occurs at 70 and 160, respectively, and 160 is taken as the minimum value after comparison. When the maximum depth is greater than 24, R 2 and· M A E converge gradually; so, the optimal parameters of the number of weak learners and the maximum depth are 160 and 24, respectively.

Multi-Layer Perceptron Model (MLPR)

We now discuss how to tune the solver, activation and hidden_layer_sizes in the MLPR model. The solvers used for weight optimization include lbfgs, which is a quasi-Newtonian optimization method; sgd, which is a stochastic gradient descent method; adam, which is an optimization method based on stochastic gradient. Excitation functions include Rectified Linear Unit function (relu), hyperbolic tangent function (tanh), logistic function and identity function. The identity function is suitable for tasks with linear relationships that lack nonlinearity, but it limits the representation of the model. Both tanh and logistic belong to generalized S-type activation functions, both of which have bidirectional saturation, but there are still some problems such as gradient disappearance and large computation in the derivation process. The relu function solves the gradient disappearance problem and is fast to compute, but it may result in the “death” of neurons. The formula for these four activation functions are as follows:
r e l u : f ( x ) = max ( 0 , x )
tanh : f ( x ) = e x e x e x + e x
log i s t i c : f ( x ) = 1 1 + e x
identity :   f ( x ) = x
where x is the input value.
In order to improve the computational efficiency and generalization ability of the model, the multi-layer perceptron here uses only one hidden layer. As shown in Figure 5a,b, in the combination of 12 solvers and excitation functions, the relu function and lbfgs solver have the best performance; so, the MLPR model should use the relu function and lbfgs solver. According to Figure 5c, it is found that R 2 increases with the increase in the number of hidden layer elements, and M A E decreases with the increase in the number of hidden layer elements. When the number of hidden layer elements is greater than 70, R 2 and· M A E gradually converge, so the optimal parameter of the number of hidden layer elements is 70.

Support Vector Machine Model (SVR)

In the SVR model, the error penalty coefficient (C) and kernel function are optimized, and the influence of parameter changes on the prediction results is evaluated by using R 2 and M A E , respectively. C is an important parameter used to balance the complexity of the model and the accuracy of the fitted data. The larger the value, the more accurately the model learns from the sample, but the easier it is to overfit. The kernel includes linear kernel functions (linear), polynomial kernel functions (poly), radial basis kernel functions (rbf) and S-type kernel functions (sigmoid). The linear kernel function is fast to calculate, but its ability to fit nonlinear data is limited. The polynomial kernel function can handle nonlinear data to some extent, but the computational complexity is high. The radial basis kernel function has a strong ability to fit nonlinear data, but it is easy to overfit. Because the error penalty coefficient is too large, the model is easy to overfit, so C is taken from 0.1 to 100. Figure 6a,b shows that the SVR model should use the radial basis kernel function (rbf), and C should take 100.

2.3. Comprehensive Analysis of Model Prediction Ability

Comparing R 2 and M A E of the optimized model, the results are shown in Figure 7. After tuning, R 2 of the three models is relatively different, and R 2 of the RF model is the highest. However, M A E of the three models is quite different, among which the RF model has the smallest M A E . Compared with the three optimized models, the RF model has excellent comprehensive performance, so the RF model should be used for the establishment of a porous titanium constitutive model.
Zhu Tianyu et al. [4] used the DP model to build a constitutive model of porous titanium, and verified the accuracy of the DP model by comparing the predicted values of the model with the experimental values. The DP model was proposed by Drucker and Prager in 1952. It was initially applied to the study of plastic mechanics of geological materials. Later, researchers found that the model parameters could be adjusted to describe the plastic deformation process of metals. In this paper, the prediction ability of the DP model and RF model for porous titanium materials is compared and analyzed.
Porous titanium is temperature-sensitive due to its substrate material. As shown in Figure 8, when the temperature is less than 300 °C, the yield limit and flow stress of porous titanium decrease with the increase in temperature, indicating that the material has a temperature-softening effect. The RF model predicted the data well, and MAPE did not exceed 3.5%. The predicted value of the DP model can reflect the temperature sensitivity, but the error is large, and the average MAPE is as high as 11.4%.
Affected by the microstructure strain efficiency and the matrix material, porous titanium has strain rate sensitivity due to the influence of its matrix material [5,32,33,34]. From the perspective of the microstructure, the pore wall of porous titanium has plastic bending during the compression process, so the strain rate effect will be produced. As shown in Figure 9, when the strain rate is less than 3000/s−1, the yield limit and flow stress of the porosity increase with the increase in the strain rate. When the strain rate reaches 3000/s−1, the distribution of the stress–strain curve is irregular. The predicted results of the RF model can reflect the strain sensitivity of porous titanium materials, and the average MAPE is 3.8%. However, the prediction results of RF show small errors in part of the stress–strain curve. The predicted value of the DP model can reflect the strain rate sensitivity of porous titanium materials when the strain rate is less than 3000/s−1, but the error is large, and the average MAPE is as high as 27.5%.
As shown in Figure 10, the presence of pores significantly reduces the mechanical strength of porous titanium materials. The RF model can predict the data well. The prediction results of the DP model can show that the existence of pores significantly reduces the mechanical strength of porous titanium materials, but the error is large.
As shown in Figure 8, Figure 9 and Figure 10, both the RF model and the traditional DP model can predict the strain rate effect and temperature sensitivity of porous titanium materials, but the prediction accuracy of the latter is lower than that of the former. Because the traditional constitutive model is usually based on physics and phenomenology and describes the structure and properties of materials through mathematical equations, its prediction accuracy is limited by its own explicit mathematical description ability. Machine learning models, on the other hand, predict the properties of materials by training large amounts of data to establish patterns and associations. As shown in Table 3, the average M A P E of the RF model is much smaller than that of the DP model, so the predictive performance of the constitutive model constructed using the RF model is better than that of the DP model.

3. Results

  • For the data set obtained by the split Hopkinson pressure bar experiment, the performance of the RF model is optimal when max_depth and n_estimators are 24 and 160, respectively. When the relu function is combined with the lbfgs solver and the number of hidden layer elements are 70, the performance of the MLPR model is optimal. When the rbf function is used and C is 100, the SVR model has the best performance.
  • Based on the data set obtained from the split Hopkinson pressure bar experiment, the prediction performance of three machine learning algorithms, RF, MLPR and SVR, on porous titanium materials is discussed. Through comparative analysis, it can be seen that the prediction performance of the RF model with 160 weak learners and a maximum depth of 24 is the best, and R 2 is 0.9998, while M A E is 1.413.
  • The traditional DP model can predict the strain rate effect and temperature sensitivity of porous titanium materials. However, the prediction accuracy of the latter was lower than that of the former, with an MAPE of 4.4% for the RF model and 25.5% for the DP model.
  • In general, the RF model has good predictive performance for the constitutive relationship of porous titanium materials.

Author Contributions

Conceptualization, C.T. and Z.L.; methodology, C.T.; software, C.T.; validation, C.T. and C.L.; formal analysis, C.T.; investigation, Z.L.; resources, Z.L.; data curation, C.T.; writing—original draft preparation, C.T.; writing—review and editing, C.T. and C.L.; visualization, C.T.; supervision, Z.L.; project administration, Z.L.; funding acquisition, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in this study are included in the article; further inquiries can be directed to the corresponding author.

Acknowledgments

The authors thank Liu Zhiqiang and Li Chunjin for providing the experimental data required for this project and guiding the writing of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
R 2 Goodness of fit
M A E Mean absolute error
M A P E Mean absolute percentage error
RFRandom forest
MLPRMulti-layer perceptron
SVRSupport vector machine

References

  1. Wang, T.; Yu, Z.; Song, Y.; Yang, L.; Chen, L.; Dang, J.; Zeng, X.; Guo, S.; Su, Y.; Song, B. Dura mater assists 3D printed porous titanium scaffolds in healing rabbit calvarial defects. Mater. Today Commun. 2023, 34, 105111. [Google Scholar] [CrossRef]
  2. Villa, T.; Tenconi, A.A.; Colombo, M.; Banci, L. Hip arthrodesis conversion to dual mobility total hip arthroplasty in a young adult using a custom-made highly-porous titanium monobloc cup: A case report. Int. J. Surg. Case Rep. 2024, 116, 109375. [Google Scholar] [CrossRef] [PubMed]
  3. Haroun, F.; Benmoussa, N.; Bidault, F.; Lassau, N.; Moya-Plana, A.; Leymarie, N.; Honart, J.-F.; Kolb, F.; Qassemyar, Q.; Gorphe, P. Outcomes of mandibular reconstruction using three-dimensional custom-made porous titanium prostheses. J. Stomatol. Oral Maxillofac. Surg. 2023, 124, 101281. [Google Scholar] [CrossRef]
  4. Zhu, T. Research on Performance Test and Micro-Cutting Simulation Experiment of Porous Titanium Material for Artificial Bone. Master’s Thesis, Jiangsu University of Science and Technology, Zhenjiang, China, 2016. [Google Scholar]
  5. Meng, S. Study on Mechanical Properties and Surface Microstructure of Porous Titanium in Artificial Bone. Master’s Thesis, Jiangsu University of Science and Technology, Zhenjiang, China, 2021. [Google Scholar] [CrossRef]
  6. Reka, S.S.; Bagelikar, A.; Venugopal, P.; Ravi, V.; Devarajan, H. Deep learning-based classification of rotten fruits and identification of shelf life. Comput. Mater. Contin. 2024, 78, 781–794. [Google Scholar] [CrossRef]
  7. Zheng, Y.; Wang, Y.; Liu, J. Research on structure optimization and motion characteristics of wearable medical robotics based on Improved Particle Swarm Optimization Algorithm. Future Gener. Comput. Syst. 2022, 129, 187–198. [Google Scholar] [CrossRef]
  8. Stergiou, K.; Ntakolia, C.; Varytis, P.; Koumoulos, E.; Karlsson, P.; Moustakidis, S. Enhancing property prediction and process optimization in building materials through machine learning: A review. Comput. Mater. Sci. 2023, 220, 112031. [Google Scholar] [CrossRef]
  9. Priyadharshini, K.V.; Vijay, A.; Swaminathan, K.; Avudaiappan, T.; Banupriya, V. Materials property prediction using feature selection based machine learning technique. Mater. Today Proc. 2022, 69, 710–715. [Google Scholar] [CrossRef]
  10. Nguyen, B.D.; Potapenko, P.; Demirci, A.; Govind, K.; Bompas, S.; Sandfeld, S. Efficient surrogate models for materials science simulations: Machine learning-based prediction of microstructure properties. Mach. Learn. Appl. 2024, 16, 100544. [Google Scholar] [CrossRef]
  11. Lin, Y.; Li, C.; Zhou, K.; Guo, Z.; Zang, C. A Constitutive Model Study of Chemical Corrosion Sandstone Based on Support Vector Machine and Artificial Bee Colony Algorithm. Sustainability 2023, 15, 13415. [Google Scholar] [CrossRef]
  12. Lei, L.; Li, B.; Wang, H.; Huang, G.; Xuan, F. High-temperature high-cycle fatigue performance and machine learning-based fatigue life prediction of additively manufactured Hastelloy X. Int. J. Fatigue 2024, 178, 108012. [Google Scholar] [CrossRef]
  13. Gao, Z.; Fan, X.; Gao, S.; Xve, W. Research on Mechanical Properties of Alloy Structural Steel Tempering Based on Multi-model Machine Learning. Mater. Rev. 2023, 37, 148–154. [Google Scholar]
  14. Barshan, D.; Ashikur, M.R.; Jahidul, M.I.; Rahman, M.Z.; Zhu, D. Properties prediction of composites based on machine learning models: A focus on statistical index approaches. Mater. Today Commun. 2024, 38, 107659. [Google Scholar]
  15. Joo, C.; Park, H.; Lim, J.; Cho, H.; Kim, J. Machine learning-based heat deflection temperature prediction and effect analysis in polypropylene composites using catboost and shapley additive explanations. Eng. Appl. Artif. Intell. 2023, 126, 106873. [Google Scholar] [CrossRef]
  16. Gao, S.; Yue, X.; Wang, H. Predictability of Different Machine Learning Approaches on the Fatigue Life of Additive-Manufactured Porous Titanium Structure. Metals 2024, 14, 320. [Google Scholar] [CrossRef]
  17. Suryawanshi, A.; Behera, N. Application of machine learning for prediction dental material wear. J. Polym. Mater. 2023, 40, 305–316. [Google Scholar] [CrossRef]
  18. Liu, Z. Structural Design of Hyperelastic Porous Materials Aided by Machine Learning. Master’s Thesis, Southwest University of Science and Technology, Mianyang, China, 2023. [Google Scholar]
  19. Zhao, P.; Zeng, X.; Kou, H.; Chen, H. A kind of numerical model combined with genetic algorithm and back propagation neural network for creep-fatigue life prediction and optimization of double-layered annulus metal hydride reactor and verification of ASME-NH code. Int. J. Hydrog. Energy 2024, 54, 1251–1263. [Google Scholar] [CrossRef]
  20. Ismail, N.; Ayman, S.; Alam, M.N.; Fathy, A. Prediction of wear rates of Al-TiO2 nanocomposites using artificial neural network modified with particle swarm optimization algorithm. Mater. Today Commun. 2023, 35, 105743. [Google Scholar]
  21. Guo, Z.; Bai, R.; Lei, Z.; Jiang, H.; Liu, D.; Zou, J.; Yan, C. CPINet: Parameter identification of path-dependent constitutive model with automatic denoising based on CNN-LSTM. Eur. J. Mech.—A/Solids 2021, 90, 104327. [Google Scholar] [CrossRef]
  22. Xia, Y.; Zhang, C.; Wang, C.; Liu, H.; Sang, X.; Liu, R.; Zhao, P.; An, G.; Fang, H.; Shi, M.; et al. Prediction of bending strength of glass fiber reinforced methacrylate-based pipeline UV-CIPP rehabilitation materials based on machine learning. Tunn. Undergr. Space Technol. 2023, 140, 105319. [Google Scholar] [CrossRef]
  23. Zhao, J.; Li, D.; Jiang, J.; Luo, P. Uniaxial compressive strength prediction for rock material in deep mine using boosting-based machine learning methods and optimization algorithms. Comput. Model. Eng. Sci. 2024, 140, 275–304. [Google Scholar] [CrossRef]
  24. Feng, Y. Constitutive Relationship Modeling Method and Experimental Verification of Advanced High-Strength Steel Based on Machine Learning. Master’s Thesis, Shanghai Jiao Tong University, Shanghai, China, 2021. [Google Scholar] [CrossRef]
  25. Li, D.-W.; Liu, J.-X.; Fan, Y.-S.; Yang, X.-G.; Huang, W.-Q. A preliminary discussion about the application of machine learning in the field of constitutive modeling focusing on alloys. J. Alloys Compd. 2023, 976, 173210. [Google Scholar] [CrossRef]
  26. Lin, Y.; Chen, X.-M. A critical review of experimental results and constitutive descriptions for metals and alloys in hot working. Mater. Des. 2011, 32, 1733–1759. [Google Scholar] [CrossRef]
  27. Brueckl, H.; Breth, L.; Fischbacher, J.; Schrefl, T.; Kuehrer, S.; Pachlhofer, J.; Schwarz, M.; Weirather, T.; Czettl, C. Machine learning based prediction of mechanical properties of WC-Co cemented carbides from magnetic data only. Int. J. Refract. Met. Hard Mater. 2024, 121, 106665. [Google Scholar] [CrossRef]
  28. Jia, X.; Hao, K.; Luo, Z.; Fan, Z. Plastic Deformation Behavior of Metal Materials: A Review of Constitutive Models. Metals 2022, 12, 2077. [Google Scholar] [CrossRef]
  29. Wang, D.; Fang, T.; Tang, Y.; Ma, Y. A Review of Support Vector Machine regression theory and control. Pattern Recognit. Artif. Intell. 2003, 16, 192–197. [Google Scholar]
  30. Xu, X.; Hu, Z.; Liu, J.; Li, W.; Liu, J. Prediction of concrete strength of Three Gorges Dam based on machine learning regression model. Mater. Rev. 2023, 37, 45–53. [Google Scholar]
  31. Dai, D.; Liu, Q.; Hu, R.; Wei, X.; Ding, G.; Xu, B.; Xu, T.; Zhang, J.; Xu, Y.; Zhang, H. Method construction of structure-property relationships from data by machine learning assisted mining for materials design applications. Mater. Des. 2020, 196, 109194. [Google Scholar] [CrossRef]
  32. Hu, Z.; Li, M.; Xu, W.; Chen, Q.; Hu, W.; Su, Y. Preparation and Mechanical properties of porous titanium by powder metallurgy. Prog. Titan. Ind. 2010, 27, 12–15. [Google Scholar] [CrossRef]
  33. Wang, X.; Li, J.; Hu, R.; Kou, H. Mechanical properties and pore structure deformation behavior of biomedical porous titanium. Trans. Nonferrous Met. Soc. China 2015, 25, 1543–1550. [Google Scholar] [CrossRef]
  34. Liu, Z.; Ji, F.; Wang, M.; Zhu, T. One-Dimensional Constitutive Model for Porous Titanium Alloy at Various Strain Rates and Temperatures. Metals 2017, 7, 24. [Google Scholar] [CrossRef]
Figure 1. Flow chart of machine learning algorithm model establishment and analysis.
Figure 1. Flow chart of machine learning algorithm model establishment and analysis.
Metals 14 00634 g001
Figure 2. Sample surface topography. Reprinted with permission from ref. [4]. 2016 Zhu, T.
Figure 2. Sample surface topography. Reprinted with permission from ref. [4]. 2016 Zhu, T.
Metals 14 00634 g002
Figure 3. Correlation matrix diagram.
Figure 3. Correlation matrix diagram.
Metals 14 00634 g003
Figure 4. Hyperparameter tuning of RF model: (a) tuning of n_estimators; (b) tuning of max_depth.
Figure 4. Hyperparameter tuning of RF model: (a) tuning of n_estimators; (b) tuning of max_depth.
Metals 14 00634 g004
Figure 5. Hyperparameter tuning of MLPR model: (a) R2 evaluation, (b) MAE evaluation and (c) hidden_layer_sizes optimization.
Figure 5. Hyperparameter tuning of MLPR model: (a) R2 evaluation, (b) MAE evaluation and (c) hidden_layer_sizes optimization.
Metals 14 00634 g005
Figure 6. Hyperparameter tuning of SVR model: (a) using R 2 for evaluation; (b) using MAE for evaluation.
Figure 6. Hyperparameter tuning of SVR model: (a) using R 2 for evaluation; (b) using MAE for evaluation.
Metals 14 00634 g006
Figure 7. Comparison between R2 and MAE after three models are optimized.
Figure 7. Comparison between R2 and MAE after three models are optimized.
Metals 14 00634 g007
Figure 8. Comparison of predictions of RF model and DP model with same porosity and strain rate at different temperatures: (a) RF-26%-1000/s−1; (b) DP-26%-1000/s−1.
Figure 8. Comparison of predictions of RF model and DP model with same porosity and strain rate at different temperatures: (a) RF-26%-1000/s−1; (b) DP-26%-1000/s−1.
Metals 14 00634 g008
Figure 9. Comparison of prediction of RF model and DP model at same porosity and temperatures and different strain rates: (a) RF-26%-200 °C, (b) DP-26%-200 °C, (c) RF-36%-100 °C and (d) DP-36%-100 °C.
Figure 9. Comparison of prediction of RF model and DP model at same porosity and temperatures and different strain rates: (a) RF-26%-200 °C, (b) DP-26%-200 °C, (c) RF-36%-100 °C and (d) DP-36%-100 °C.
Metals 14 00634 g009
Figure 10. Comparison of prediction of RF model and DP model at same strain rate and temperatures and different porosity: (a) RF-100 °C-1300/s−1; (b) DP-100 °C-1300/s−1.
Figure 10. Comparison of prediction of RF model and DP model at same strain rate and temperatures and different porosity: (a) RF-100 °C-1300/s−1; (b) DP-100 °C-1300/s−1.
Metals 14 00634 g010
Table 1. Specific parameters of materials.
Table 1. Specific parameters of materials.
PorosityApertureGrain SizeTiFeCuCON
26%15 μ m ≤27 μ m ≥99.7%≤0.15%≤0.005%≤0.05%≤0.2%≤0.03%
36%250 μ m ≤74 μ m ≥99.7%≤0.25%≤0.003%≤0.06%≤0.2%≤0.03%
Table 2. Sample data information.
Table 2. Sample data information.
Porosity/%Temperature/°CStrain Rate/s−1
26251200, 2000, 3000, 3600, 5200
100950, 1200, 2200, 3000, 4200
200800, 1500, 1950, 2750, 3800
3001200, 2000, 2900,3600, 3700
36251000, 2000, 3000
1001380, 2050, 2350, 3400, 3700
2001000,1800, 2000, 2400, 3000
3001200, 2000, 3000, 3400, 4500
Table 3. Error comparison between DP model and RF model.
Table 3. Error comparison between DP model and RF model.
Model Maximum   M A E Maximum   M A P E Mean   M A E Mean   M A P E
RF22.58.0%10.94.4%
DP185.944.9%94.425.5%
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

Tan, C.; Li, C.; Liu, Z. Application of Machine Learning in Constitutive Relationship Prediction of Porous Titanium Materials for Artificial Bone. Metals 2024, 14, 634. https://doi.org/10.3390/met14060634

AMA Style

Tan C, Li C, Liu Z. Application of Machine Learning in Constitutive Relationship Prediction of Porous Titanium Materials for Artificial Bone. Metals. 2024; 14(6):634. https://doi.org/10.3390/met14060634

Chicago/Turabian Style

Tan, Chengzhi, Chunjin Li, and Zhiqiang Liu. 2024. "Application of Machine Learning in Constitutive Relationship Prediction of Porous Titanium Materials for Artificial Bone" Metals 14, no. 6: 634. https://doi.org/10.3390/met14060634

APA Style

Tan, C., Li, C., & Liu, Z. (2024). Application of Machine Learning in Constitutive Relationship Prediction of Porous Titanium Materials for Artificial Bone. Metals, 14(6), 634. https://doi.org/10.3390/met14060634

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