Biomechanical evaluation of different strain judging criteria on the prediction precision of cortical bone fracture simulation under compression
Ruoxun Fan
Jie Liu2 - 1Department of Traffic Engineering, Yangzhou Polytechnic Institute, Yangzhou, China
- 2Department of Aerospace Engineering, Jilin Institute of Chemical Technology, Jilin, China
- 3Department of Mechanical and Aerospace Engineering, Jilin University, Changchun, China
Introduction: The principal strain or equivalent strain is mainly used in current numerical studies to determine the mechanical state of the element in the cortical bone finite element model and then perform fracture simulation. However, it is unclear which strain is more suitable for judging the element mechanical state under different loading conditions due to the lack of a general strain judging criterion for simulating the cortical bone fracture.
Methods: This study aims to explore a suitable strain judging criterion to perform compressive fracture simulation on the rat femoral cortical bone based on continuum damage mechanics. The mechanical state of the element in the cortical bone finite element model was primarily assessed using the principal strain and equivalent strain separately to carry out fracture simulation. The prediction accuracy was then evaluated by comparing the simulated findings with different strain judging criteria to the corresponding experimental data.
Results: The results showed that the fracture parameters predicted using the principal strain were closer to the experimental values than those predicted using the equivalent strain.
Discussion: Therefore, the fracture simulation under compression was more accurate when the principal strain was applied to control the damage and failure state in the element. This finding has the potential to improve prediction accuracy in the cortical bone fracture simulation.
Introduction
With a load-bearing role, the cortical bone structure bears the majority of the external load and protects the soft tissue from damage (Kumar and Ghosh, 2022). Both the instantaneous impact and continuous static loads may lead to structural damage and crack occurrence in the cortical bone structure, and even direct fracture in extreme cases (Cui et al., 2022; Rosa et al., 2022). Therefore, investigating the relationships among the cortical bone mechanical properties, damage and failure mechanisms, and different loading conditions is necessary to explore preventive measures to reduce the possibility of cortical bone damage (Nourisa et al., 2019). However, the achievement of these first requires accurate observation of the complete failure process in the cortical bone structure.
The failure process in the cortical bone structure includes a series of mechanical behaviors from elastic deformation to crack initiation, propagation, and final fracture (Salem et al., 2022; Snow et al., 2022). Due to the limitation of image development, observing these mechanical behaviors at the micro-level by experiment is difficult (Ural, 2021). Therefore, most studies simulate the failure process of the cortical bone structure using finite element (FE) analysis. Early numerical research predicted the fracture load in the cortical bone structure by developing a linear elastic fracture model (Stitzel et al., 2004; Ural and Vashishth, 2006). With the development of fracture mechanics, the elastic-plastic fracture model was used to perform cortical bone fracture, which simulated the failure process by describing the change in the mechanical behavior during crack propagation (Prasannavenkadesan and Pandithevan, 2021; Remache et al., 2020). The cortical bone fracture can also be simulated using the continuum damage mechanics theory and the extended finite element method. Both methods can simulate the complete failure process in the cortical bone structure by defining the crack initiation, softening model, and element failure, and can simulate both brittle and ductile fractures by setting different softening models (Ng et al., 2017; Soni et al., 2021; Yadav et al., 2021; Kumar et al., 2022).
The above numerical fracture models differed but shared one feature: the strain was regarded as an indicator to control the damage and failure state in the element. Many studies have considered that strain is more accurate in determining the mechanical state of bone elements than stress (Cohen et al., 2017; Lovrenic-Jugovic et al., 2020). However, different strain judging criteria were used among the fracture models discussed above. One part applied the principal strain as the element damage and failure criterion in the cortical bone FE model, while the other part used the equivalent strain (Ng et al., 2017; Remache et al., 2020; Soni et al., 2021; Kumar et al., 2022). Because no standardized strain judging criterion for cortical bone fracture simulation was provided, it is unclear which strain is more suitable for judging the element mechanical state in the FE models under different loading conditions.
Cortical bone is a transversely isotropic material with mechanical properties that vary depending on the orientation of the osteon (Li et al., 2013; Szabo and Rimnac, 2022). Meanwhile, the main difference between the equivalent strain and principal strain is whether the shear strain is considered, so the numerical variation process in the two types of strain is different under various loading conditions (Xia and Wang, 2001). Therefore, the applicable range of the two strains in judging the element mechanical state should differ in cortical bone fracture simulation under different loading conditions.
This study aims to explore a suitable strain judging criterion for simulating cortical bone fracture under compression. The compressive experiment on the rat femoral cortical bone specimen was first performed. Then, the cortical bone FE model was established based on the femur micro images, and the equivalent and principal strains were separately set as the judging criteria to simulate the failure process in the cortical bone structure. Finally, the predicted fracture parameters were compared with the corresponding experimental data, and the effects of applying different strain judging criteria on the prediction precision of cortical bone fracture were discussed, the outcome of which has the potential to improve prediction accuracy in cortical bone fracture simulation under compression.
Materials and methods
Compressive experiment
Six 3-month-old male Wistar rats were selected to obtain six right femurs. The location of the femoral mid-diaphysis was set as a benchmark, and the cortical bone specimen with a height of 5 mm was cut along the femoral shaft axis. All the cortical bone specimens were subjected to axial compression on the electronic testing machine, as shown in Figure 1A. The pressure head was compressed at a uniform speed of 1 mm/min to perform quasi-static loading, and the test started with an increase to 30 N, followed by a decrease to 0 N to ensure no slide on the specimen (Fang et al., 2019).
FIGURE 1. The schematic diagrams of the experiment and simulation on the rat femoral cortical bone structure (A) Compressive experiment; (B) Fracture simulation.
Establishment of the cortical bone finite element model
The cortical bone micro images were obtained based on micro-CT scanning of the rat femurs. The images were first imported into the MIMICS software to reconstruct the geometry model of the femoral cortical bone structure, then the geometry model was imported into the ABAQUS software to establish the cortical bone FE model using the C3D8 element. To simulate the boundary condition in the compressive experiment, rigid circular plates were created above and below the cortical bone structure. The frictionless interaction was established between the lower rigid plate and the cortical bone structure, and the tangential behavior friction with a penalty coefficient of 0.2 was created between the upper rigid plate and the cortical bone structure (Hambli and Allaoui, 2013; Harrison et al., 2013). A reference point coupled to the upper rigid plate was set up above the model. An axial compression displacement of 0.5 mm was applied to the reference point, and all the degrees of freedom except the loading direction on the reference point were constrained. Meanwhile, all the degrees of freedom in the lower rigid plate were constrained. The boundary condition in the FE models can be seen in Figure 1B. Assigning the longitudinal and transverse elastic moduli for the FE model is necessary because the cortical bone is a transversely homogeneous material. After a previous study on this batch of rat femurs, i.e., the nanoindentation test on these cortical bone specimens, the average longitudinal and transverse elastic moduli of the cortical bone tissue were measured as 13,498 and 11,025 MPa, respectively, and the Poisson’s ratio was set to 0.3 (Zhang et al., 2015).
Fracture simulation for the cortical bone finite element model
The fracture method based on continuum damage mechanics theory was adopted in this study. The implemented ABAQUS user’s subroutine (UMAT) was related to the damage variable D and used to perform fracture simulation. The stress–strain relationship after the onset of the damage can be expressed as (Cornin et al., 2022):
where
The damage variable D was set to depend on the principal strain or equivalent strain. The change in the damage variable can be expressed as (Gaziano et al., 2022):
Where ε is the principal strain or equivalent strain in the element of the FE model, εf is the critical failure strain in the cortical bone material.
The critical failure strain in the femoral cortical bone of 3-month-old Wistar rats was obtained in our previous study, with a compressive failure strain threshold of 4.3% (Fan et al., 2016). As the compression condition was considered in this study, the compressive strain in the element was dominant, so the minimum principal strain was compared with the equivalent strain to determine the suitable strain judging criterion. In the process of compression simulation, the minimum principal strain and equivalent strain were separately compared with the critical compressive failure strain threshold of the cortical bone material to determine the element damage initiation. Element damage occurred when the minimum principal or equivalent strain exceeded 4.3%. In the damaged element, the damage variable D increased as expression 3), and the elastic modulus decreased as expression 2). The damage elasticity matrix was adopted to update the stress and the Jacobian matrix in the element. The Jacobian matrix, after the onset of the damage, could be described as expression 4) (Gao et al., 2018). The maximum value of the damage variable was set to 0.9 (Hambli et al., 2012). The element failed and lost its load-bearing capacity when the elastic modulus reduced to 10% of the original. The apparent failure occurred in the cortical bone structure when the failed elements accumulated to a certain degree. The flow chart of the fracture simulation can be seen in Figure 2. By comparing the numerical differences between the simulation results obtained by applying the minimum principal and equivalent strains, the strain judging criterion that is more suitable for simulating the fracture of the rat femoral cortical bone under compression can be determined.
FIGURE 2. Detailed flow chart for simulating the progressive failure in the cortical bone structure using UMAT.
Results
Mesh sensitivity analysis
Mesh sensitivity analysis was first conducted to determine the suitable element size for the FE models. Five sizes (60, 80, 100, 120, and 140 μm) were used for creating five cortical bone FE models, and the fracture simulation was performed by applying the equivalent strain to observe the effects of different sizes on the failure process. Figure 3 shows that the fracture load increased with the fine mesh, indicating that the damage variable D rose faster for the large element size and led to a faster decrease in the structural stiffness. Meanwhile, when the element size was in the range of 60–100 μm, the shapes of the load–displacement curves did not differ much. The fine mesh was needed because the fracture method adopted in this study cannot make a crack in the element. Therefore, the element size was selected to 80 μm for creating the cortical bone FE models. The element number of the six established cortical bone FE models was in the range of 15,230–16,843 in this study.
FIGURE 3. Mesh sensitivity analysis on the cortical bone FE models.
Comparison of the fracture patterns in the simulation and experiment
Figure 4 compares the fracture patterns in the simulations and experiment. The simulation using the equivalent strain expressed a transverse slide-open fracture at the end of compression. The crack was produced in the middle of the structure and propagated in the direction of perpendicular to the compression load. In contrast, the simulation predicted by applying the minimum principal strain exhibited a longitudinal slide-open fracture, and the crack propagated along a certain angle to the loading direction. The comparisons of the fracture patterns stated that the simulation predicted by applying the minimum principal strain was in better agreement with the experimental results.
FIGURE 4. Comparison of the fracture patterns in the simulations and compressive experiment (A) Fracture simulation predicted by applying the equivalent strain; (B) Fracture simulation predicted by applying the minimum principal strain; (C) Compressive experiment.
Comparison of the predicted curves between the simulations and experiments
Figure 5 shows the load–displacement curves obtained from the experiments and the corresponding simulations separately by applying the minimum principal and equivalent strains. The experimental curves exhibited that the quasi-brittle fracture occurred in the rat femoral cortical bone structure under compression. The yielding behavior was not appear after the elastic phase, and the apparent stiffness also hardly changed before the complete fracture. The simulation curves also exhibited quasi-brittle fracture in the cortical bone FE models. On the premise that the elastic moduli in the cortical bone FE models were all acquired from the nanoindentation test, the apparent stiffness predicted by the simulations using different strain judging criteria agreed with the experimental results.
FIGURE 5. Comparison of the load–displacement curves in the six specimens and corresponding FE models obtained from the experiments and simulations separately using the minimum principal strain and equivalent strain.
The shapes of the predicted curves obtained by applying the two strain judging criteria were essentially similar, but the fracture parameters differed. At the beginning of loading, the curves coincided because the assigned elastic moduli were the same. With increasing compression, the simulation performed by applying the minimum principal strain first entered the fracture stage, while the simulation performed by applying the equivalent strain underwent a relatively long elastic phase before the complete fracture. Therefore, the fracture time in the cortical bone FE model predicted by applying the equivalent strain was remarkably later than that predicted by applying the minimum principal strain, which directly led to great differences in the fracture load. Furthermore, the comparison with the experimental curves showed that the simulated results obtained by applying the minimum principal strain were close to the experimental results, while the results predicted by applying the equivalent strain were different from the experimental data.
Discussion
Adopting different strain criteria in judging the element mechanical state in the FE model may lead to differences in the predicted results, even for the same structure (Russell et al., 2012; Yavari et al., 2013). To explore the suitable strain judging criterion in compression conditions, this study simulated the cortical bone fracture under compression load based on the continuous damage mechanics theory. The prediction accuracy of the fracture simulation using the minimum principal and equivalent strains was first determined, then the reasons for the differences in the mechanical responses of the same cortical bone structure under the two simulations were exhibited and explained.
It was observed experimentally that quasi-brittle fracture occurred in the femoral cortical bone structure under compression in this study. Based on the literature, for brittle fracture, the consistency of the failure process can be verified by comparing the strength limit, the apparent failure strain, and the apparent stiffness in the elastic phase (Dall'Ara et al., 2013; Huang et al., 2010). Due to the lack of a significant plastic yielding phase in the quasi-brittle failure process, the load–displacement relationships were more suited to conducting the comparison among the fracture parameters. In this study, the predicted fracture load, fracture time, and apparent stiffness in the elastic phase were all compared to the experimental data to discuss the prediction accuracy between the two simulations. Furthermore, the comparisons of the fracture patterns between the simulations and experiments were also used to verify the conclusions obtained from the match of the load–displacement curves. The comparisons of the load–displacement curves and the fracture patterns stated that the simulation predicted by applying the minimum principal strain was in better agreement with the experimental results. The fracture parameters and the direction of the crack propagation in the simulation using the minimum principal strain were similar to the experimental data. Therefore, compared with the equivalent strain, the application of the minimum principal strain to control the element mechanical state was more suitable for simulating the cortical bone fracture under compression conditions.
The discrepancy in the predicted results must be caused by adopting different strain judging criteria because the structural and material parameters in the two simulations were identical. The cortical bone structures were subjected to an axial compression load. The minimum principal strain in the element should be the unidirectional strain similar to the loading direction, while the equivalent strain was calculated by considering the transverse, longitudinal, and shear strains (Onaka, 2010; Hambli and Allaoui, 2013). Thus, the intrinsic reason for the differences in the results may depend on whether the shear strain was considered and which strain was dominant during the compression process.
When element damage occurred, the differences in the crack initiation location appeared in the two simulations, as shown in Figure 6. The oblique crack, which was expressed at a certain angle to the loading direction, appeared in the simulation predicted by applying the minimum principal strain. This phenomenon was mainly because the minimum principal strain was unidirectional in the load direction, and the load direction may transmit at a certain angle to the osteon in the cortical bone specimen due to friction during compression (Hambli and Allaoui, 2013). However, the transverse crack, which was perpendicular to the loading direction, appeared in the simulation predicted by applying the equivalent strain, indicating that the shear strain in the element increased faster and dominated at the crack initiation stage (Feerick et al., 2013). Both the cracks then propagated with increasing compression. Once the crack initiation location is confirmed, the driving force will expand in the direction of the crack and will not easily change direction due to the limitation of the boundary condition (Zimmermann et al., 2009). Thus, one propagated along a certain angle to the loading direction, and the other propagated perpendicular to the loading direction. Because transverse propagation required crossing many osteons, which needed more energy than longitudinal propagation, the predicted results expressed a greater fracture load in the simulation using the equivalent strain than that of using the minimum principal strain.
FIGURE 6. The schematic diagrams of the failure processes in the two simulations (A) Fracture simulation predicted by applying the minimum principal strain; (B) Fracture simulation predicted by applying the equivalent strain.
It can be inferred that the equivalent strain was not suitable for judging the element mechanical status and performing the compressive fracture simulation for the cortical bone structure because the shear strain in the element increased faster than the longitudinal strain at the beginning of the loading and dominated the driving force at the crack initiation stage. Based on these, the innovation of this paper was to further improve the simulation accuracy of the cortical bone fracture under compression. The principal strain or equivalent strain was always selected as the criterion to judge the crack initiation and propagation in the element of the FE model (Dall’Ara and Tozzi, 2022; Hambli and Allaoui, 2013). However, various types of strain rise differently, even under the same loading condition (Dall’Ara et al., 2017). The inappropriate application of the strain judging criterion for performing fracture simulation may result in differences in the damage initiation and propagation processes, thus affecting the whole failure process. Therefore, the strain judging criterion suitable for simulating compressive fracture of the cortical bone structure was explored in this study, which may further improve the prediction accuracy of fracture simulation on the basis of previous simulations.
The fracture simulations in this study had several limitations. First, only compressive load was discussed, and the conclusions obtained in this study were only applicable to this loading condition. The cortical bone may also be subjected to other loads, e.g., bending and torsion in daily activity, and the change in the element strain in the cortical bone FE model may differ in various loading conditions. Therefore, the three-point bending simulation and experiment on the rat femur are planned to be conducted in the future to explore the suitable strain judging criterion in fracture simulation under bending load. Second, only the minimum principal strain in the element was considered because the cortical bone was subjected to a compressive load. Although the cortical bone structure was relatively simple and the bone unit aligned in the same direction as the loading direction, generating tensile strain in the element was also possible (Li et al., 2013; Dai et al., 2022). However, the accuracy of the simulation results in this study can be illustrated by comparing them with the experimental results. These indicated that the tensile strain generated in the cortical bone FE model under compressive load was relatively low, and the failure to consider the tensile strain did not have a substantial effect on the simulation results.
This study evaluated the prediction accuracy of the fracture simulations predicted by applying the minimum principal and equivalent strains. Experimental verification showed that the application of the minimum principal strain to control the damage and failure state in the element of the FE model could accurately simulate the failure process in the cortical bone structure under compression. Meanwhile, based on the comparisons of the load–displacement curves and fracture patterns, the reasons for the differences between the simulation using the equivalent strain and the experimental results were also revealed. Thus, this study provides a foundation for exploring the suitable strain judging criterion for cortical bone fracture under compression, which has the potential to improve prediction accuracy in the bone biomechanics field.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by The Medical Ethics Committee of the First Hospital of Jilin University (No. 2018-238).
Author contributions
RF implemented the material model in the finite element code and simulated the various boundary conditions. JL and ZJ developed models, assisted in data analysis, and generated figures for the manuscript. All authors contributed to this manuscript and approved the submitted version.
Funding
The work was supported by the Natural Science Foundation of Jilin Province (YDZJ202301ZYTS250).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Cohen, H., Kugel, C., May, H., Medlej, B., Stein, D., Slon, V., et al. (2017). The influence of impact direction and axial loading on the bone fracture pattern. Forensic Sci. Int. 277, 197–206. doi:10.1016/j.forsciint.2017.05.015
Cronin, D. S., Watson, B., Khor, F., Gierczycka, D., and Malcolm, S. (2022). Cortical bone continuum damage mechanics constitutive model with stress triaxiality criterion to predict fracture initiation and pattern. Front.Bioeng. Biotech. 56, 1022506–1023374. doi:10.3389/fbioe.2022.1022506
Cui, Y. Y., Xiang, D. D., Shu, L. M., Duan, Z., Liao, Z., Wang, S., et al. (2022). Incremental element deletion-based finite element analysis of the effects of impact speeds, fall postures, and cortical thicknesses on femur fracture. Materials 15 (8), 2878. doi:10.3390/ma15082878
Dai, X. F., Fan, R. X., Wu, H. J., and Jia, Z. B. (2022). Investigation on the differences in the failure processes of the cortical bone under different loading conditions. Appl. Bionics. Biomech. 2022, 1–9. doi:10.1155/2022/3406984
Dall'Ara, E., Luisier, B., Schmidt, R., Kainberger, F., Zysset, P., and Pahr, D. (2013). A nonlinear QCT-Based finite element model validation study for the human femur tested in two configurations in vitro. Bone 52 (1), 27–38. doi:10.1016/j.bone.2012.09.006
Dall’Ara, E., Peña-Fernández, M., Palanca, M., Giorgi, M., Cristofolini, L., and Tozzi, G. (2017). Precision of digital volume correlation approaches for strain analysis in bone imaged with micro-computed to mography at different dimensional levels. Front. Bioeng. Biotech. 4, 31. doi:10.3389/fmats.2017.00031
Dall’Ara, E., and Tozzi, G. (2022). Digital volume correlation for the characterization of musculoskeletal tissues: Current challenges and future developments. Front. Bioeng. Biotech. 10, 1010056. doi:10.3389/fbioe.2022.1010056
Fan, R. X., Gong, H., Zhang, R., Gao, J. Z., Jia, Z. B., and Hu, Y. J. (2016). Quantification of age-related tissue-level failure strains of rat femoral cortical bones using an approach combining macrocompressive test and microfinite element analysis. J. Biomech. Eng-T. ASME 138 (4), 041006. doi:10.1115/1.4032798
Fang, J., Gao, J. Z., Gong, H., Zhang, T. L., Zhang, R., and Zhan, B. C. (2019). Multiscale experimental study on the effects of different weight-bearing levels during moderate treadmill exercise on bone quality in growing female rats. Biomed. Eng. Online 18, 33. doi:10.1186/s12938-019-0654-1
Feerick, E. M., Liu, X. Y., and Mcgarry, P. (2013). Anisotropic mode-dependent damage of cortical bone using the extended finite element method (XFEM)finite element method (XFEM). J. Mech. Behav. Biomed. Mat. 20, 77–89. doi:10.1016/j.jmbbm.2012.12.004
Gao, X., Chen, M. H., Yang, X. G., Harper, L., Ahmed, I., and Lu, J. W. (2018). Simulating damage onset and evolution in fully bio-resorbable composite under three-point bending. J. Mech. Behav. Biomed. Mat. 81, 72–82. doi:10.1016/j.jmbbm.2018.02.022
Gaziano, P., Falcinelli, C., and Vairo, G. (2022). A computational insight on damage-based constitutive modelling in femur mechanics. Eur. J. Mech. A-Solids 93, 104538. doi:10.1016/j.euromechsol.2022.104538
Hambli, R., and Allaoui, S. (2013). A robust 3D finite element simulation of human proximal femur progressive fracture under stance load with experimental validation. Ann. Biomed. Eng. 41 (12), 2515–2527. doi:10.1007/s10439-013-0864-9
Hambli, R., Bettamer, A., and Allaoui, S. (2012). Finite element prediction of proximal femur fracture pattern based on orthotropic behaviour law coupled to quasi-brittle damage. Med. Eng. Phys. 34 (2), 202–210. doi:10.1016/j.medengphy.2011.07.011
Harrison, N. M., McDonnell, P., Mullins, L., Wilson, N., O’Mahoney, D., and McHugh, P. E. (2013). Failure modelling of trabecular bone using a non-linear combined damage and fracture voxel finite element approachfinite element approach. Biomech. Model. Mechan. 12 (2), 225–241. doi:10.1007/s10237-012-0394-7
Huang, H. L., Tsai, M. T., Lin, D. J., Chien, C. S., and Hsu, J. T. (2010). A new method to evaluate the elastic modulus of cortical bone by using a combined computed tomography and finite element approach. Comput. Biol. Med. 40 (4), 464–468. doi:10.1016/j.compbiomed.2010.02.011
Kumar, A., and Ghosh, R. (2022). A review on experimental and numerical investigations of cortical bone fracture. Proc. Inst. Mech.Eng. H. 236 (3), 297–319. doi:10.1177/09544119211070347
Kumar, A., Shitole, P., Ghosh, R., Kumar, R., and Gupta, A. (2022). Experimental and numerical comparisons between finite element method, element-free Galerkin method, and extended finite element method predicted stress intensity factor and energy release rate of cortical bone considering anisotropic bone modelling. Proc. Inst. Mech. Eng. H. 233 (8), 823–838. doi:10.1177/0954411919853918
Li, S., Demirci, E., and Silberschmidt, V. V. (2013). Variability and anisotropy of mechanical behavior of cortical bone in tension and compression. J. Mech. Behav. Biomed. Mat. 21, 109–120. doi:10.1016/j.jmbbm.2013.02.021
Lovrenic-Jugovic, M., Tonkovic, Z., and Skozrit, I. (2020). Experimental and numerical investigation of cyclic creep and recovery behavior of bovine cortical bone. Mech. Mat. 146, 103407. doi:10.1016/j.mechmat.2020.103407
Ng, T. P., Koloor, S. S. R., Djuansjah, J. R. P., and Abdul Kadir, M. (2017). Assessment of compressive failure process of cortical bone materials using damage-based model. J. Mech. Behav. Biomed. Mat. 66, 1–11. doi:10.1016/j.jmbbm.2016.10.014
Nourisa, J., and Rouhi, G. (2019). Prediction of the trend of bone fracture healing based on the results of the early stages simulations: A finite element study. J. Mech. Med. Biol. 19 (5), 1950021. doi:10.1142/S0219519419500210
Onaka, S. (2010). Equivalent strain in simple shear deformation described by using the Hencky strain. Phil. Mag. Lett. 90 (9), 633–639. doi:10.1080/09500839.2010.489030
Prasannavenkadesan, V., and Pandithevan, P. (2021). Johnson-cook model combined with cowper-symonds model for bone cutting simulation with experimental validation. J. Mech. Med. Biol. 21 (2), 2150010. doi:10.1142/S021951942150010X
Remache, D., Semaan, M., Rossi, J., Pithioux, M., and Milan, J. (2020). Application of the Johnson-Cook plasticity model in the finite element simulations of the nanoindentation of the cortical bone. J. Mech. Behav. Biomed. Mat. 101, 103426. doi:10.1016/j.jmbbm.2019.103426
Rosa, N., Moura, M. F. S. F., Olhero, S., Simoes, R., Magalhães, F. D., Marques, A. T., et al. (2022). Bone: An outstanding composite material. Appl. Sci-Basel. 12 (7), 3381. doi:10.3390/app12073381(
Russell, C. M., Choo, A. M., Tetzlaff, W., Chung, T. E., and Oxland, T. R. (2012). Maximum principal strain correlates with spinal cord tissue damage in contusion and dislocation injuries in the rat cervical spine. J. Neurotraum 29 (8), 1574–1585. doi:10.1089/neu.2011.2225
Salem, M., Westover, L., Adeeb, S., and Duke, K. (2022). Prediction of fracture initiation and propagation in pelvic bones. Comput. Method. Biomec. 25 (7), 808–820. doi:10.1080/10255842.2021.1981883
Snow, T., Woolley, W., Metcalf, R. M., Rosenberg, J., Acevedo, C., and Kingstedt, O. T. (2022). Effect of collagen damage induced by heat treatment on the mixed-mode fracture behavior of bovine cortical bone under elevated loading rates. Int. J. Fract. 233 (1), 85–101. doi:10.1007/s10704-021-00612-0
Soni, A., Negi, A., Kumar, S., and Kumar, N. (2021). An IGA based nonlocal gradient-enhanced damage model for failure analysis of cortical bone. Eng. Fract. Mech. 255, 107976. doi:10.1016/j.engfracmech.2021.107976
Stitzel, J. D., Barretta, J. T., and Duma, S. M. (2004). Predicting fractures due to blunt impact: A sensitivity analysis of the effects of altering failure strain of human rib cortical bone. Int. J. Crashworthiness 9 (6), 633–642. doi:10.1533/ijcr.2004.0317
Szabo, E., and Rimnac, C. (2022). Biomechanics of immature human cortical bone: A systematic review. J. Mech. Behav. Biomed. Mat. 125, 104889. doi:10.1016/j.jmbbm.2021.104889
Ural, A. (2021). Biomechanical mechanisms of atypical femoral fracture. J. Mech. Behav. Biomed. Mat. 124, 104803. doi:10.1016/j.jmbbm.2021.104803
Ural, A., and Vashishth, D. (2006). Cohesive finite element modeling of age-related toughness loss in human cortical bone. J. Biomech. 39 (16), 2974–2982. doi:10.1016/j.jbiomech.2005.10.018
Xia, K., and Wang, J. (2001). Shear, principal, and equivalent strains in equal-channel angular deformation. Metal. Mater.Trans. A 32 (10), 2639–2647. doi:10.1007/s11661-001-0054-5
Yadav, R. N., Uniyal, P., Sihota, P., Kumar, S., Dhiman, V., Goni, V. G., et al. (2021). Effect of ageing on microstructure and fracture behavior of cortical bone as determined by experiment and extended finite element method (XFEM). Med. Eng. Phys. 93, 100–112. doi:10.1016/j.medengphy.2021.05.021
Yavari, S. A., Weinans, H., and Zadpoor, A. (2013). Full-field strain measurement and fracture analysis of rat femora in compression test. J. Biomech. 46 (7), 1282–1292. doi:10.1016/j.jbiomech.2013.02.007
Zhang, R., Gong, H., Zhu, D., Ma, R. S., Fang, J., and Fan, Y. B. (2015). Multi-level femoral morphology and mechanical properties of rats of different ages. Bone 76, 76–87. doi:10.1016/j.bone.2015.03.022
Keywords: cortical bone, compression, fracture, principal strain, equivalent strain
Citation: Fan R, Liu J and Jia Z (2023) Biomechanical evaluation of different strain judging criteria on the prediction precision of cortical bone fracture simulation under compression. Front. Bioeng. Biotechnol. 11:1168783. doi: 10.3389/fbioe.2023.1168783
Received: 18 February 2023; Accepted: 03 April 2023;
Published: 13 April 2023.
Edited by:
Lizhen Wang, Beihang University, ChinaReviewed by:
Yongtao Lu, Dalian University of Technology (DUT), ChinaAbdelwahed Barkaoui, International University of Rabat, Morocco
Copyright © 2023 Fan, Liu and Jia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ruoxun Fan, fanruoxun@126.com