- Research article
- Open Access
- Open Peer Review
Prediction of fracture load and stiffness of the proximal femur by CT-based specimen specific finite element analysis: cadaveric validation study
BMC Musculoskeletal Disorders volume 18, Article number: 536 (2017)
Finite element analysis (FEA) of the proximal femur has been previously validated with large mesh size, but these were insufficient to simulate the model with small implants in recent studies. This study aimed to validate the proximal femoral computed tomography (CT)-based specimen-specific FEA model with smaller mesh size using fresh frozen cadavers.
Twenty proximal femora from 10 cadavers (mean age, 87.1 years) were examined. CT was performed on all specimens with a calibration phantom. Nonlinear FEA prediction with stance configuration was performed using Mechanical Finder (mesh,1.5 mm tetrahedral elements; shell thickness, 0.2 mm; Poisson’s coefficient, 0.3), in comparison with mechanical testing. Force was applied at a fixed vertical displacement rate, and the magnitude of the applied load and displacement were continuously recorded. The fracture load and stiffness were calculated from force–displacement curve, and the correlation between mechanical testing and FEA prediction was examined.
A pilot study with one femur revealed that the equations proposed by Keller for vertebra were the most reproducible for calculating Young’s modulus and the yield stress of elements of the proximal femur. There was a good linear correlation between fracture loads of mechanical testing and FEA prediction (R2 = 0.6187) and between the stiffness of mechanical testing and FEA prediction (R2 = 0.5499). There was a good linear correlation between fracture load and stiffness (R2 = 0.6345) in mechanical testing and an excellent correlation between these (R2 = 0.9240) in FEA prediction.
CT-based specimen-specific FEA model of the proximal femur with small element size was validated using fresh frozen cadavers. The equations proposed by Keller for vertebra were found to be the most reproducible for the proximal femur in elderly people.
Proximal femoral neck fracture is a major cause of high morbidity and mortality in elderly people with osteoporosis [1, 2]. A decline in the bone mineral density (BMD) is associated with fracture; hence, to predict fracture risk, dual energy x-ray absorptiometry and quantitative computed tomography (CT) have been widely used to measure BMD [3,4,5,6,7]. However, BMD is correlated only with bone strength and cannot indicate other mechanical properties of the proximal femur . Furthermore, its ability to predict bone strength varies [9,10,11]. CT-based finite element analysis (FEA) can account for various aspects of the bone, such as bone geometry, cortical and trabecular bone distribution, and loading direction, and can improve the predictive accuracy of bone strength. The first validation study was conducted by Keyak et al. , and several studies have reported the accuracy and usefulness of CT-based FEA of the proximal femur in the stance configuration [13,14,15,16,17].
Recently, a large number of studies using FEA of the proximal femur have been reported as a postoperative evaluation with small implants like screws, plates, and the femoral stem [18,19,20,21,22,23]. In future, these studies may be useful for preoperative planning to select the most suitable implant or fixation methods. However, FEA for the bone behavior with implants is controversial because mesh sizes of models that were used in previous validation studies were relatively larger than implant sizes. Thus, FEA models should be validated with a smaller mesh size to improve their reliability.
This study aimed to document a proximal femoral FEA model using fresh frozen cadavers with smaller mesh size and to verify the accuracy of FEA prediction using CT-based specimen-specific FEA compared with mechanical testing.
Twenty femora (right, 10 and left, 10) from 10 fresh frozen cadavers (males, 5 and females, 5) were obtained from the Clinical Anatomy Laboratory in our university. The mean age was 87.1 years [range, 74–101 years; standard deviation (SD), 9.13]. Subjects underwent no previous hip surgeries. The cause of death included pneumonia (n = 3); senility (n = 2); heart failure (n = 2); and chronic renal failure, breast cancer, and pancreatic cancer (n = 1 each). Cadavers were stored at −22 °C, and after thawing at room temperature, whole femurs were retrieved from the body and all soft tissues were removed. CT was performed using Aquilion ONE (Toshiba Medical Systems, Tokyo, Japan) with imaging parameters including 320-row detector; 120 kV; 200 mA; slice thickness, 0.5 mm; and pixel width, 0.3 mm. A calibration phantom (QRM-BDC, QRM, Möhrendorf, DE) containing three hydroxyapatite rods (0, 100, and 200 mg/cm3) was tested together with the specimen in water. The proximal femora were then sawed 12 cm distal from the tip of the greater trochanter, sloping 20° in the coronal plane to the shaft axis. Specimens were sprayed with a saline solution to maintain their moisture during the procedure and were not refrozen.
To verify fracture load assessment, quasi-static compression testing was conducted. Specimens were loaded using a universal testing machine (Autograph AG-20kN X Plus; Shimadzu, Kyoto, Japan). The distal 3-cm portion of each specimen parallel to the sawed surface was fixed using resin cement (resin box), and a resin cap was molded on the femoral head to apply a uniform compressive load (Fig. 1). The force was applied on the resin cap at a fixed vertical displacement rate of 5 mm/min until proximal femoral fracture occurred. The magnitude of the applied load and displacement were continuously recorded, and mechanical fracture was identified when the slope of the force–displacement curve (stiffness) rapidly decreased. Mechanical stiffness was calculated between 20% and 80% of the maximum fracture load using the force–displacement curve.
Nonlinear FEA prediction
Data from CT were transferred to an HP Z400 workstation (Hewlett-Packard, Palo Alto, CA, USA). The proximal femur model (12 cm distal from the tip of the greater trochanter) was made using FEA software (Mechanical Finder, Research Center for Computational Mechanics, Tokyo, Japan). All femoral trabecular bone and inner parts of the cortex were meshed using linear tetrahedral elements with 1.5-mm global edge length and overlaid with 1.5 × 1.5 × 0.2-mm triangular shell elements simulating the outer cortex.
CT value of each element was set as the average of the voxels contained in one element. Mechanical properties of each element were calculated in Hounsfield units (HU) . There are several equations for calculating the Young’s modulus and yield stress for proximal femur. Therefore, we first performed a pilot study using CT DICOM data from an 85-year-old male to determine the most reproducible equations for proximal femur at this mesh size, and the most reproducible equations were adopted in subsequent tests. Modulus values of <0.01 MPa were designated as 0.01 MPa, and those >20 GPa as 20 GPa . Young’s modulus and yield stress of the shell element were calculated, assuming its CT value was 1000 HU. Drucker–Prager equivalent criterion was adopted for the yield of the element . The tensile yield stress was assumed to be 0.8 times the compressive yield stress, in agreement with previous studies [15, 27, 28]. Poisson’s coefficient for each element was set at 0.3 .
To reproduce the real mechanical testing, the FEA model was sloped at 20° in the coronal plane to the shaft axis and fixed with a resin box 3 cm distally, and a resin cap was placed on the femoral head. A uniaxial compressive load with a uniform distribution was applied on the resin cap. The degree of displacement and the reaction force of each point were recorded, and the force–displacement curve was determined. FEA-predicted fracture load was defined as the load when stiffness had declined by >20% of the estimated values, and FEA-predicted stiffness was calculated in the same way as for the actual mechanical test.
Fracture load and stiffness were compared in mechanical testing and by FEA prediction using Student’s t tests and Pearson χ2 tests (SPSS 16.0, SPSS, Inc., Chicago, Illinois, USA). Age, gender, and side differences were calculated using multiple regression analysis as predictive factors for the fracture load in the proximal femur. A p-value of <0.05 was considered significant.
The force–displacement curves obtained from mechanical testing and FEA prediction using equations proposed by Keyak  (Fig. 2a), Carter  (Fig. 2b), Minamisawa  (Fig. 2c), and Keller for vertebra  (Fig. 2d) are shown. Considering the fracture load (the tip of the curve) and stiffness (the slope of the curve), we concluded that Keller’s equation for vertebra was the most reproducible for the proximal femur. The following specific equations were used:
Young’s modulus (E, MPa).
Yield stress (σ, MPa).
Accuracy of FEA prediction in fracture load
The mean fracture loads from mechanical testing and FEA prediction using Keller’s equations for vertebra were 3435.8 N (SD, 1802.1) and 4520.3 N (SD, 1879.0), respectively. There was a good linear correlation between these values (Fig. 3).
Accuracy of FEA prediction in stiffness
The mean stiffness from mechanical testing and FEA prediction were 1280 N/mm (SD, 847) and 1566 N/mm (SD, 664), respectively. There was also a good linear correlation between these values (Fig. 4).
Correlation between fracture load and stiffness
Laterality in mechanical testing and FEA prediction
Fracture load and stiffness were not significantly different between the right and left femora for both mechanical testing and FEA prediction (Table 1).
Effect of age on mechanical testing and FEA prediction
There was a significant linear correlation between fracture load and age from mechanical testing; the load declined 142.6 N per year (Fig. 6a). FEA prediction also produced a significant linear correlation between fracture load and age, with load declining by 153.8 N per year (Fig. 6b).
FEA prediction was validated with a CT-based subject-specific proximal femoral FEA model with a smaller element size using fresh frozen cadavers. This study adopted a finer size tetrahedral element (1.5 mm) for the entire model than that reported in previous studies [13,14,15,16,17]. FEA models in previous studies used a 3.0-mm cube-shaped element , a l.5-mm cube-shaped element only below the lesser trochanter , a 3.0-mm tetrahedral element , or a 1.5–4.0-mm tetrahedral element [16, 17]. These element sizes are no longer sufficient to match the progress of modern implant technology and are insufficient to analyze the proximal femur with small screws or a stem-inserted model of the femur. Cube-shaped elements can make a mechanically precise model, but the model cannot be automatically made and is time-consuming to construct. An adequate model can be made using tetrahedral elements. At least 1.5-mm tetrahedral elements of the entire model are essential to simulate the proximal femur with small implants.
Preliminary testing revealed that Keller’s equations for vertebra were the most reproducible for calculating Young’s modulus and yield stress for the proximal femur. Although Keller’s equations for vertebra predicted approximately 1.6 times greater stiffness than mechanical testing, Keyak’s equations, which were adopted in several previous studies, predicted approximately 10 times greater stiffness than mechanical testing. We consider the error of Keller’s equations for vertebra to be acceptable and better than that reported in other studies. The correspondence between the proximal femur and vertebra derives from the high proportion of sponge-like cancellous bone with a thin cortical bone shell in both. Moreover, both proximal femur and vertebra present comparable risks for fracture in people with osteoporosis.
Correlations between mechanical testing and FEA prediction were significant for both fracture load and stiffness in this study. However, previous studies have shown extremely high correlations (fracture load, R2 = 0.73–0.96 and stiffness, R2 = 0.62–0.82) [13,14,15,16,17]. These discrepancies in the results may probably be due to differences in the equations adopted for FEA models. Smaller element size might be another reason for the discrepancy. It has been shown that differences in element size influence the results of FEA predictions [15, 33]. We adopted a compressive displacement condition to reproduce the results from mechanical testing, not to simulate the actual fracture situation, which may also have influenced the results. Furthermore, bone strength is determined not only by BMD but also by bone quality. Bone quality includes bone turnover, microarchitecture, collagen quality, mineralization, microdamage, and bone matrix and mineral composition [34,35,36]. Approximately 70% of bone strength is dependent on BMD and 30% on bone quality . CT-based FEA can calculate the influence of BMD, but it does not consider all aspects of bone quality, which suggests that there is some systematic error inherent to FEA prediction. Nevertheless, we believe that our results are an acceptable estimate of mechanical behavior.
Good correlations between fracture load and stiffness have been reported previously . This suggests that the fracture load of each specimen depends on its stiffness rather than on the degree of displacement at the point of fracture. The correlation between load and stiffness was far higher with FEA prediction than with mechanical testing. Bone quality, which was not considered in CT-based FEA, might have influenced these results.
Our study had certain limitations. First, most subjects were aged >80 years; hence, these results cannot be applied to younger people. However, FEA prediction of proximal femoral mechanical properties is most useful for elderly people at the risk of proximal femoral fractures. Second, we did not examine the fall configuration that was the most suitable to simulate actual hip fracture with a high correlation between mechanical and FEA-predicted values [12, 16, 38,39,40]. We only adopted the one-legged stance configuration to examine the behavior of the proximal femur because our focus was to validate the model and make it possible to analyze the model with small implants. Recent studies using FEA prediction of the proximal femur with small implants have provided much important knowledge [18,19,20,21,22,23]; however, the validation was insufficient. Third, we did not compare mechanical fracture sites with FEA-predicted fracture sites because most mechanical fracture sites were not obvious macroscopically or with CT imaging. The fracture sites could be evaluated with micro CT, but it was not practical or necessary for this study.
FEA model of the proximal femur with small element size was validated using fresh frozen cadavers. The equations proposed by Keller for vertebra were found to be the most reproducible for the proximal femur in elderly people.
Bone mineral density
Finite element analysis
Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int. 2006 Dec;17(12):1726–33.
Cummings SR, Melton LJ. Epidemiology and outcomes of osteoporotic fractures. Lancet. 2002 May 18;359(9319):1761–7.
Hayes WC. Biomechanics of cortical and trabecular bone: implications for assessment of fracture risk. In: Hayes WC, editor. Basic orthopaedics biomechanics. Raven press; 1991. p. 93–142.
Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. Report of a WHO Study Group. World Health Organ Tech Rep Ser. 1994;843:1–129. Review.
Adams JE. Quantitative computed tomography. Eur J Radiol. 2009 Sep;71(3):415–24.
Black DM, Bouxsein ML, Marshall LM, Cummings SR, Lang TF, Cauley JA, Ensrud KE, Nielson CM, Orwoll ES. Osteoporotic fractures in men (MrOS) research group. Proximal femoral structure and the prediction of hip fracture in men: a large prospective study using QCT. J Bone Miner Res. 2008 Aug;23(8):1326–33.
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. Vivo discrimination of hip fracture with quantitative computed tomography: results from the prospective European femur fracture study (EFFECT). J Bone Miner Res. 2011 Apr;26(4):881–93.
Ammann P, Rizzoli R. Bone strength and its determinants. Osteoporos Int. 2003;14(Suppl 3):S13–8.
Bouxsein ML, Coan BS, Lee SC. Prediction of the strength of the elderly proximal femur by bone mineral density and quantitative ultrasound measurements of the heel and tibia. Bone. 1999 Jul;25(1):49–54.
Lochmüller EM, Müller R, Kuhn V, Lill CA, Eckstein F. Can novel clinical densitometric techniques replace or improve DXA in predicting bone strength in osteoporosis at the hip and other skeletal sites? J Bone Miner Res. 2003 May;18(5):906–12.
Faulkner KG, Cummings SR, Black D, Palermo L, Glüer CC, Genant HK. Simple measurement of femoral geometry predicts hip fracture: the study of osteoporotic fractures. J Bone Miner Res. 1993 Oct;8(10):1211–7.
Keyak JH, Rossi SA, Jones KA, Skinner HB. Prediction of femoral fracture load using automated finite element modeling. J Biomech. 1998 Feb;31(2):125–33.
Keyak JH. Improved prediction of proximal femoral fracture load using nonlinear finite element models. Med Eng Phys. 2001 Apr;23(3):165–73.
Cody DD, Gross GJ, Hou FJ, Spencer HJ, Goldstein SA, Fyhrie DP. Femoral strength is better predicted by finite element models than QCT and DXA. J Biomech. 1999 Oct;32(10):1013–20.
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.
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 Jan;52(1):27–38.
Dragomir-Daescu D, Salas C, Uthamaraj S, Rossman T. Quantitative computed tomography-based finite element analysis predictions of femoral strength and stiffness depend on computed tomography settings. J Biomech. 2015 Jan 2;48(1):153–61.
Noda M, Saegusa Y, Takahashi M, Tezuka D, Adachi K, Naoi K. Biomechanical study using the finite element method of internal fixation in Pauwels type III vertical femoral neck fractures. Arch Trauma Res. 2015 Aug 26;4(3):e23167.
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.
Chang CW, Chen YN, Li CT, Peng YT, Chang CH. Role of the compression screw in the dynamic hip-screw system: a finite-element study. Med Eng Phys. 2015 Dec;37(12):1174–9.
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 Feb;46(2):213–7.
Peng MJ, Chen HY, Hu Y, Ju X, Bai B. Finite element analysis of porously punched prosthetic short stem virtually designed for simulative uncemented hip arthroplasty. BMC Musculoskelet Disord. 2017 Jul 11;18(1):295.
Burchard R, Braas S, Soost C, Graw JA, Schmitt J. Bone preserving level of osteotomy in short-stem total hip arthroplasty does not influence stress shielding dimensions - a comparing finite elements analysis. BMC Musculoskelet Disord. 2017 Aug 7;18(1):343.
Les CM, Keyak JH, Stover SM, Taylor KT, Kaneps AJ. Estimation of material properties in the equine metacarpus with use of quantitative computed tomography. J Orthop Res. 1994;12:822–33.
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.
Drucker DC, Prager W. Soil mechanics and plastic analysis of limit design. Q Appl Math. 1952;10:157–65.
Kaneko TS, Pejcic MR, Tehranzadeh J, Keyak JH. Relationships between material properties and CT scan data of cortical bone with and without metastatic lesions. Med Eng Phys. 2003 Jul;25(6):445–54.
Keaveny TM, Wachtel EF, Ford CM, Hayes WC. Differences between the tensile and compressive strengths of bovine tibial trabecular bone depend on modulus. J Biomech. 1994;27:1137–46.
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.
Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure. J Bone Joint Surg Am. 1977 Oct;59(7):954–62.
Minamisawa T. Japanese only. 1981.
Keller TS. Predicting the compressive mechanical behavior of bone. J Biomech. 1994 Sep;27(9):1159–68.
Keyak JH, Skinner HB. Three-dimensional finite element modelling of bone: effects of element size. J Biomed Eng. 1992;14:483–9.
Fonseca H, Moreira-Gonçalves D, Coriolano HJ, Duarte JA. Bone quality: the determinants of bone strength and fragility. Sports Med. 2014 Jan;44(1):37–53.
Gourion-Arsiquaud S, Faibish D, Myers E, Spevak L, Compston J, Hodsman A, Shane E, Recker RR, Boskey ER, Boskey AL. Use of FTIR spectroscopic imaging to identify parameters associated with fragility fracture. J Bone Miner Res. 2009 Sep;24(9):1565–71.
Riggs BL, Hodgson SF, O'Fallon WM, Chao EY, Wahner HW, Muhs JM, Cedel SL, Melton LJ 3rd. Effect of fluoride treatment on the fracture rate in postmenopausal women with osteoporosis. N Engl J Med 1990 Mar 22;322(12):802–809.
NIH Consensus Development Panel on Osteoporosis Prevention, Diagnosis, and Therapy. Osteoporosis prevention, diagnosis, and therapy. JAMA 2001 Feb 14;285(6):785–795. Review.
Varga P, Schwiedrzik J, Zysset PK, Fliri-Hofmann L, Widmer D, Gueorguiev B, Blauth M, Windolf M. Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup. J Mech Behav Biomed Mater. 2016 Apr;57:116–27.
Dragomir-Daescu D, Op Den Buijs J, McEligot S, Dai Y, Entwistle RC, Salas C, Melton LJ 3rd, Bennet KE, Khosla S, Amin S. Robust QCT/FEA models of proximal femur stiffness and fracture load during a sideways fall on the hip. Ann Biomed Eng 2011 Feb;39(2):742–755.
Nishiyama KK, Gilchrist S, Guy P, Cripton P, Boyd SK. Proximal femur bone strength estimated by a computationally fast finite element analysis in a sideways fall configuration. J Biomech. 2013 Apr 26;46(7):1231–6.
Junichi Nakamura received JSPS KAKENHI Grant Number 25870125. Study sponsors had no involvement in study design, data collection or interpretation, writing of the manuscript. The other authors have nothing to disclose.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
The research protocol of this cadaveric study was in compliance with the Helsinki Declaration, approved by the Institutional Review Boards (Research Ethics Committees of Graduate School of Medicine, Chiba University; the reference number #1782) and written consent was obtained for the donation and participate the study before death.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.