- Research article
- Open Access
Finite element analysis of the femoral diaphysis of fresh-frozen cadavers with computed tomography and mechanical testing
© The Author(s). 2018
- Received: 17 May 2018
- Accepted: 23 July 2018
- Published: 31 July 2018
The purpose of this study was to validate a diaphyseal femoral fracture model using a finite element analysis (FEA) with mechanical testing in fresh-frozen cadavers.
We used 18 intact femora (9 right and 9 left) from 9 fresh-frozen cadavers. Specimens were obtained from 5 males and 4 females with a mean age of 85.6 years. We compared a computed tomography (CT)-based FEA model to diaphyseal femoral fracture loads and stiffness obtained by three-point bending. Four material characteristic conversion equations (the Keyak, Carter, and Keller equations plus Keller’s equation for the vertebra) with different shell thicknesses (0.3, 0.4, and 0.5 mm) were compared with the mechanical testing.
The average fracture load was 4582.8 N and the mean stiffness was 942.0 N/mm from actual mechanical testing. FEA prediction using Keller’s equation for the vertebra with a 0.4-mm shell thickness showed the best correlations with the fracture load (R2 = 0.76) and stiffness (R2 = 0.54). Shell thicknesses of 0.3 and 0.5 mm in Keller’s equation for the vertebra also showed a strong correlation with fracture load (R2 = 0.66 for both) and stiffness (R2 = 0.50 and 0.52, respectively). There were no significant correlations with the other equations.
We validated femoral diaphyseal fracture loads and stiffness using an FEA in a cadaveric study.
- Finite element analysis
- The femoral diaphysis fracture
- Validation study
- Fresh frozen cadaver
- Computed tomography
Although proximal femoral fractures are common in the elderly, femoral shaft fractures are seen in all generations. Because of the difference of fracture type and bone quality, we may need individual consideration of proximal femoral fracture and femoral diaphyseal fracture.
Recently, computed tomography-based finite element analysis (CT-based FEA) has been widely used for mechanical analysis of the femur. Many reports have described fracture models of the femur following traffic accidents and simulations of stress distribution after joint replacement [1–13]. When we consider a femoral diaphyseal fracture and the stress distribution for femoral diaphyseal after prosthesis replacement using FEA, the whole femoral shaft must be evaluated. However, there are no reports of FEA models of diaphyseal femoral fractures. Although experimental data can be used to validate CT-based FEA models, most of the fracture models in previous reports were of the proximal femur [1–6]. The proximal of femur have abundant cancellous bone, while the femoral diaphysis mainly consists of cortical bone. The bone architecture of the proximal femur and the diaphysis are different, and thus, they should be examined separately.
The aim of this study is to validate the newly constructed CT-based FEA models of the femur by comparing with the data obtained from the actual mechanical fracture tests using the original fresh-frozen cadaveric femurs.
The research protocol was in compliance with the Helsinki Declaration; it was approved by the Research Ethics Committee of our institution and registered with the University Hospital Medical Information Network. Written informed consent was obtained from all the donors before death.
The specimens were loaded using the Autograph AG-20000N X Plus Precision Universal Tensile Tester (Shimadzu, Kyoto, Japan). The proximal and the distal 2.5-cm portions of the specimen were embedded in resin cement boxes shaped square parallel to the cut surface. A three-point bending test was done without fixing the resin cement boxes to the jig (Fig. 1b). A compression load was applied at the midpoint of the specimen on the anterior surface of the femoral shaft with a crosshead speed of 5 mm/min. Load was applied until the yield point was confirmed on the load-deformation curve, which was continuously recorded in 100 Hz throughout the mechanical tests. Fracture was defined by a rapid decrease on the load-displacement curve. The yield point was considered as the time point where a fracture of the specimen occurred. The load recorded at the yield point was defined as the fracture load. The stiffness was calculated as the slope of the load-deformation curve between 20 and 80% of the fracture load reference to the previous study that was reported by Miura et al. .
Non-linear FEA study
CT-based morphological data were transmitted to an HP Z400 workstation (Hewlett-Packard, Palo Alto, CA, USA). A non-linear FEA model of the femur was constructed from the CT data using Mechanical Finder software (Research Center for Computational Mechanics, Tokyo, Japan). The FEA model was made to exactly match the specimen used for the mechanical test (proximal end, 12 cm distal to the tip of the greater trochanter; distal end, 20 cm distal to the proximal end). The cancellous bone and the inner cortical bone were simulated with 3-mm linear tetrahedral elements according to the previous reports , and the outer cortical bone was constructed with 3 × 3 × 0.3–0.5 mm triangular shell elements. The shell elements were to reinforce the margin of the bone, because the tetrahedral elements of the bone margin might show lower HU values as they calculated the average of HU values in the element. Three models were created for each specimen with different shell thicknesses of 0.3, 0.4, and 0.5 mm.
Equations proposed by Keyak, Carter, Keller, and Keller for vertebra
Young’s modulus (E: MPa)
Yield stress (σ: MPa)
E = 0.001 (ρ = 0)
σ = 1.0 × 1020 (ρ ≦ 0.2)
E = 33,900 ρ2.20 (0 < ρ ≦ 0.27)
E = 5307 ρ + 469 (0.27 < ρ < 0.6)
E = 10,200 ρ2.01 (0.6 ≦ ρ)
σ = 137 ρ1.88 (0.2 < ρ < 0.317), σ = 114ρ1.72 (0.317 ≦ ρ)
E = 0.001 (ρ = 0)
σ = 1.0 × 1020 (ρ ≦ 0.2)
E = 3790(0.01)0.05 ρ (ρ < 0)
σ = 68(0.01)0.06 ρ2 (0.2 < ρ)
E = 0.001 (ρ = 0)
σ = 1.0 × 1020 (ρ ≦ 0.2)
E = 10,500 ρ2.51 (ρ < 0)
σ = 117 ρ1.93 (0.2 < ρ)
Keller for vertebra
E = 0.001 (ρ = 0)
σ = 1.0 × 1020 (ρ ≦ 0.2)
E = 1890 ρ1.92 (ρ < 0)
σ = 284 ρ2.27 (0.2 < ρ)
FEA results were compared with the results of mechanical tests using Pearson’s correlation coefficients. Influence of the shell thickness was evaluated for the most highly correlated Keller-vertebra equation. A p value of less than 0.05 was considered statistically significant (BellCurve for Excel, Social Survey Research Information Co., Ltd. Tokyo, Japan).
All femoral fractures occurring in mechanical tests were transverse fractures at the center of the femoral shafts. The FEA model also reproduced the same fracture site.
Shell thickness effect
The FEA was developed in 1970, and three-dimensional analysis using CT started around 1990 . In recent years, the FEA has been widely used in the field of orthopedic surgery. There are several validation studies reported, including those for the proximal femur, the vertebra, and the distal radius [1–3, 5, 17, 21–23]. In validation studies of the proximal femur, Bessho et al. and Keyak reported a strong positive correlation between the FEA results and actual mechanical tests (R2 = 0.96 and 0.94, respectively) [2, 5]. For the vertebra and distal radius, Imai et al. and Matsuura et al. also reported positive correlations (R2 = 0.96 and 0.97, respectively) [17, 21]. However, to our knowledge, a validation study for the femoral diaphysis has not been reported. In the present study, the FEA showed significant correlations with the femoral shaft fracture load (R2 = 0.76) and stiffness (R2 = 0.54) calculated using the Keller-vertebra equation with a 0.4-mm shell thickness. The correlation coefficients in our study were not as high as those in the previous reports [2, 5, 17, 21]. It may be due to the fact that those studies excluded an analysis of the stiffness. In the studies that investigated stiffness, the values were similar to the current study. Dall'Ara et al.  showed strong similarity between the FEA and actual mechanical tests for fracture load (R2 = 0.72) and stiffness (R2 = 0.54) in a validation study of the proximal femur. In a vertebral study, Matsuura et al. also reported lower values for fracture load (R2 = 0.48) and stiffness (R2 = 0.79). We believe this is the first report on a successful validation of the FEA model of the femoral diaphysis.
Keyak’s equation has been widely used in the previous FEA reports [2, 3, 5, 9–11, 21]. However, we could not find any correlation between the mechanical tests and values produced by the Keyak equation, nor by Carter’s or Keller’s equations, except for the Keller-vertebra equation. To explain this, we calculated stiffness and fracture load using two equations, Keyak’s equation and the Keller-vertebra equation. The ash density of cortical bone of the femoral diaphysis was substituted with 0.935 g/cm3 by reference to Bousson’s report . As a result, stiffness showed about a fivefold difference between the Keyak and the Keller-vertebra equations (8911 and 1661 MPa, respectively). Moreover, the fracture load calculated using Keyak’s equation was less than half of that using Keller’s equation for the vertebra (101 and 243 MPa, respectively). Therefore, Keyak’s bone model seems too stiff to reproduce elasticity in bending. Using Keyak’s model, the fractures likely occur at an early stage due to the low yield stress. Although many studies have used Keyak’s equation, the present study suggests that it may provide an inaccurate estimate of the mechanical properties of the femoral diaphysis. We propose that the Keller-vertebra equation is a suitable material characteristic conversion equation for the FEA of the femoral diaphysis when an FEA of femoral diaphyseal fractures is applied to examine stress distribution.
The shell thickness has often been defined as 0.3–0.4 mm in previous validation studies [5, 17, 23]. However, the optimum shell thickness is unknown. The current study showed that a shell thickness of 0.4 mm had the best correlation coefficient for the femoral diaphysis both for fracture load and stiffness. Further study is needed to determine the optimum shell thickness for other bone sites.
The present study has several limitations. First, only fractures from three-point bending tests were evaluated. Because of the limited number of cadavers, we could not perform other fracture tests. Clinically, there are various types of femoral diaphyseal fractures. Further study is necessary to validate our FEA for other types of fractures, such as spiral fractures or those from axial compression forces. Secondly, the steel bar model might not be made to exactly match the bar used for the mechanical test. But, in both cases, the load was applied to the center point of the femoral diaphyseal, and the behavior of the bone was similar, so we thought that the research had no influence. Third, all cadavers were from elderly individuals. It is possible that the fracture load or stiffness would be different in younger individuals. However, it is difficult to obtain femurs from young people. We studied matched pairs of right and left femurs in nearly equal numbers of men and women, and we observed high reproducibility without differences of laterality (data not shown). Thus, we think that the FEA of the femoral diaphysis is predictive of mechanical properties regardless of laterality or gender.
In conclusion, we have confirmed that FEA can predict fracture loads of the femoral diaphysis. We propose that Keller’s equation for the vertebra with a 0.4-mm shell thickness allows prediction of fracture load and stiffness.
This study was supported by JSPS KAKENHI grant number 17K10954. Senior author, Junichi Nakamura, was supported by JSPS KAKENHI grant number 17K10954 and Takeda Science Foundation.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
YW performed all the testing and analysis, interpreted the data, and wrote manuscript as the corresponding author. MM supported all the testing and analysis. JN and YM conceived of this study, participated in its design and coordination, and helped to draft the manuscript. YW, MM, and TS participated in obtaining the specimens from the cadavers and performing mechanical testing. SH, SOrita, KY, YK, MS, KI, and KN participated in the design and performed the statistical analysis. SOhtori finally checked this manuscript. All authors have read and approved this final submitted manuscript.
Ethics approval and consent to participate
The research protocol was in compliance with the Helsinki Declaration, approved by the Research Ethics Committees of Graduate School of Medicine, Chiba University (reference number 670) and registered with the University Hospital Medical Information Network. Written informed consent was obtained from all study participants.
Consent for publication
Informed consent for publication was obtained from the families of all the participants.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Miura M, Nakamura J, Matsuura Y, Wako Y, Suzuki T, Hagiwara S, Orita S, Inage K, Kawarai Y, Sugano M, Nawata K, Ohtori S. Prediction of fracture load and stiffness of the proximal femur by CT-based specimen specific finite element analysis: cadaveric validation study. BMC Musculoskelet Disord. 2017;18(1):536.View ArticlePubMedPubMed CentralGoogle Scholar
- Keyak JH, Rossi SA, Jones KA, Skinner HB. Prediction of femoral fracture load using automated finite element modeling. J Biomech. 1998;31(2):125–33.View ArticlePubMedGoogle Scholar
- Keyak JH. Improved prediction of proximal femoral fracture load using nonlinear finite element models. Med Eng Phys. 2001;23(3):165–73.View ArticlePubMedGoogle Scholar
- Keyak JH, Rossi SA, Jones KA, Les CM, Skinner HB. Prediction of fracture location in the proximal femur using finite element models. Med Eng Phys. 2001;23:657–64.View ArticlePubMedGoogle Scholar
- Bessho M, Ohnishi I, Matsuyama J, Matsumoto T, Imai K, Nakamura K. Prediction of strength and strain of the proximal femur by a CT-based finite element method. J Biomech. 2007;40(8):1745–53.View ArticlePubMedGoogle Scholar
- Dall'Ara E, Luisier B, Schmidt R, Kainberger F, Zysset P, Pahr D. A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro. Bone. 2013;52(1):27–38.View ArticlePubMedGoogle Scholar
- Kennedy EA, Hurst WJ, Stitzel JD, Cormier JM, Hansen GA, Smith EP, Duma SM. Lateral and posterior dynamic bending of the mid-shaft femur: fracture risk curves for the adult population. Stapp Car Crash J. 2004;48:27–51.PubMedGoogle Scholar
- Li Z, Zou D, Liu N, Zhong L, Shao Y, Wan L, Huang P, Chen Y. Finite element analysis of pedestrian lower limb fractures by direct force: the result of being run over or impact? Forensic Sci Int. 2013;229(1–3):43–51.View ArticlePubMedGoogle Scholar
- Hirata Y, Inaba Y, Kobayashi N, Ike H, Fujimaki H, Saito T. Comparison of mechanical stress and change in bone mineral density between two types of femoral implant using finite element analysis. J Arthroplast. 2013;28(10):1731–5.View ArticleGoogle Scholar
- Ike H, Inaba Y, Kobayashi N, Hirata Y, Yukizawa Y, Aoki C, Choe H, Saito T. Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses. Comput Methods Biomech Biomed Engin. 2015;18(10):1056–65.View ArticlePubMedGoogle Scholar
- Oba M, Inaba Y, Kobayashi N, Ike H, Tezuka T, Saito T. Effect of femoral canal shape on mechanical stress distribution and adaptive bone remodelling around a cementless tapered-wedge stem. Bone Joint Res. 2016;5(9):362–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hrubina M, Horak Z, Skotak M, Letocha J, Baca V, Dzupa V. Assessment of complications depending on the sliding screw position - finite element method analysis. Bratisl Lek Listy. 2015;116(5):302–10.PubMedGoogle Scholar
- Leonidou A, Moazen M, Lepetsos P, Graham SM, Macheras GA, Tsiridis E. The biomechanical effect of bone quality and fracture topography on locking plate fixation in periprosthetic femoral fractures. Injury. 2015;46(2):213–7.View ArticlePubMedGoogle Scholar
- Keyak JH, Skinner HB. Three-dimensional finite element modelling of bone: effects of element size. J Biomed Eng. 1992;14:483–9.View ArticlePubMedGoogle Scholar
- Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure. J Bone Joint Surg Am. 1977;59(7):954–62.View ArticlePubMedGoogle Scholar
- Keller TS. Predicting the compressive mechanical behavior of bone. J Biomech. 1994;27(9):1159–68.View ArticlePubMedGoogle Scholar
- Matsuura Y, Kuniyoshi K, Suzuki T, Ogawa Y, Sukegawa K, Rokkaku T, Takahashi K. Accuracy of specimen-specific nonlinear finite element analysis for evaluation of distal radius strength in cadaver material. J Orthop Sci. 2014;19(6):1012–8.View ArticlePubMedGoogle Scholar
- Drucker DC, Prager W. Soil mechanics and plastic analysis of limit design. Q Appl Math. 1952;10:157–65.View ArticleGoogle Scholar
- Muller M, Mitton D, Moilanen P, Bousson V, Talmant M, Laugier P. Prediction of bone mechanical properties using QUS and pQCT: study of the human distal radius. Med Eng Phys. 2008;30:761–7.View ArticlePubMedGoogle Scholar
- Zysset PK, Dall'Ara E, Varga P, Pahr DH. Finite element analysis for prediction of bone strength. Bonekey Rep. 2013;2:386.View ArticlePubMedPubMed CentralGoogle Scholar
- Imai K, Ohnishi I, Bessho M, Nakamura K. Nonlinear finite element model predicts vertebral bone strength and fracture site. Spine (Phila Pa 1976). 2006;31(16):1789–94.View ArticleGoogle Scholar
- Matsuura Y, Giambini H, Ogawa Y, Fang Z, Thoreson AR, Yaszemski MJ, Lu L, An KN. Specimen-specific nonlinear finite element modeling to predict vertebrae fracture loads after vertebroplasty. Spine (Phila Pa 1976). 2014;39(22):E1291–6.View ArticleGoogle Scholar
- Matsuura Y, Kuniyoshi K, Suzuki T, Ogawa Y, Sukegawa K, Rokkaku T, Thoreson AR, An KN, Takahashi K. Accuracy of specimen-specific nonlinear finite element analysis for evaluation of radial diaphysis strength in cadaver material. Comput Methods Biomech Biomed Engin. 2015;18(16):1811–7.View ArticlePubMedGoogle Scholar
- Bousson VD, Adams J, Engelke K, Aout M, Cohen-Solal M, Bergot C, Haguenauer D, Goldberg D, Champion K, Aksouh R, Vicaut E, Laredo JD. In vivo discrimination of hip fracture with quantitative computed tomography: results from the prospective European Femur Fracture Study (EFFECT). J Bone Miner Res. 2011;26(4):881–93.View ArticlePubMedGoogle Scholar