Assessment of Hip Fracture Risk Using Cross-Section Strain Energy Determined by QCT-Based Finite Element Modeling

Kheirollahi, Hossein; Luo, Yunhua

doi:https://doi.org/10.1155/2015/413839

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion References Copyright Related Articles

Research Article Corrigendum

!

A Corrigendum for this article has been published. To view the article details, please click the ‘Corrigendum’ tab above.

Special Issue

Emerging Insights on Surgical Techniques and Biomaterials for Total Hip and Knee Arthroplasty

View this Special Issue

Research Article | Open Access

Volume 2015 | Article ID 413839 | https://doi.org/10.1155/2015/413839

Assessment of Hip Fracture Risk Using Cross-Section Strain Energy Determined by QCT-Based Finite Element Modeling

Hossein Kheirollahi¹and Yunhua Luo^1,2

Academic Editor: Kengo Yamamoto

Received24 Apr 2015

Revised28 Jun 2015

Accepted29 Jun 2015

Published25 Oct 2015

Abstract

Accurate assessment of hip fracture risk is very important to prevent hip fracture and to monitor the effect of a treatment. A subject-specific QCT-based finite element model was constructed to assess hip fracture risk at the critical locations of femur during the single-leg stance and the sideways fall. The aim of this study was to improve the prediction of hip fracture risk by introducing a novel failure criterion to more accurately describe bone failure mechanism. Hip fracture risk index was defined using cross-section strain energy, which is able to integrate information of stresses, strains, and material properties affecting bone failure. It was found that the femoral neck and the intertrochanteric region have higher fracture risk than other parts of the femur, probably owing to the larger content of cancellous bone in these regions. The study results also suggested that women are more prone to hip fracture than men. The findings in this study have a good agreement with those clinical observations reported in the literature. The proposed hip fracture risk index based on strain energy has the potential of more accurate assessment of hip fracture risk. However, experimental validation should be conducted before its clinical applications.

1. Introduction

Two of the major determinants of proximal femur fractures among the elderly are osteoporosis and sideways fall. The most common serious injury associated with the fall of an elderly person is hip fracture. Furthermore, hip fracture is associated with an up to 20% chance of death, a 25% chance of long term institutionalization, and less than a 50% chance of full recovery [1]. The total number of hip fractures in men and women in 1990 was estimated to be 338,000 and 917,000, respectively, over the world [2]. Assuming no change in the age- and sex-specific incidence, the number of hip fractures is estimated to approximately double to 2.6 million by the year 2025 and 4.5 million by the year 2050 over the world [2]. By the increasing trend in hip fractures because of the aging of the population, the worldwide annual costs of hip fractures in the year 2050 have been estimated to be $131.5 billion [3]. So, hip fracture risk should be assessed in individuals who are at a risk for an osteoporotic hip fracture to provide proper plans to prevent future probable fractures.

Statistical models and bone mineral density (BMD) captured by Dual-Energy X-ray Absorptiometry (DXA) are used for in vivo osteoporotic fracture risk assessment [4, 5]. However, their accuracy in assessment of fractures is limited; most osteoporotic fractures actually occur with BMD measurements that are above the conventional osteoporotic threshold [6]. Fracture Risk Assessment Tool (FRAX) is a tool to evaluate an individual’s fracture probability in the next 10 years, adopted by the WHO in 2008 [7]. FRAX does not take into account fall-induced impact force that is critically important in the hip fracture risk assessment [8, 9]. The main limitations of the FRAX include the following: it is a statistical model and fracture risk is not consistent within 10 years with some of the treatment results [10]. Hip structure analysis (HSA) program is now commercially available and is used to automatically assess the geometric and structural parameters of the femur. Whereas HSA is based on a beam model, the deformation of bone is oversimplified, especially in the femoral neck and the intertrochanteric region where osteoporotic femoral fractures most often occur; therefore the deformation is too complicated to be described as a beam model [11].

Finite element models constructed from quantitative computed tomography (QCT) are very helpful in assessing hip fracture risk, as they are based on well-established biomechanical principles and theories. In QCT-based finite element modeling, choosing a proper failure criterion is very important to accurate assessment of hip fracture risk. Commonly used bone failure criteria include von Mises stress and strain criteria [12–15], maximum principle stress and strain criteria [16–19], maximum shear stress criterion [20], and maximum distortion energy criterion [20]. The femur consists of inhomogeneous (porous) cancellous bone and nearly homogenous cortical bone, and their failure mechanisms are different due to their different microstructure. Failure mechanism of the cancellous bone is mostly in the form of spicule (trabecula) buckling, and the failure of denser cancellous bone and the cortical bone is mostly characterized by local cracking [21, 22]. Cortical and cancellous bone have different properties. Cortical bone is usually more brittle [23], while cancellous bone is more ductile. The current failure criteria are not convenient in describing their properties. Therefore, the total strain energy, which integrates all information of stress, strain, and material properties, is a better option [21]. Mirzaei et al. [21, 24] predicted failure strength and failure patterns of human proximal femur and human vertebrae using the strain energy criterion with a QCT-based FE model. Their predictions of the failure loads and failure locations were in a good agreement with the experimental findings. The strain energy criterion is widely used in fracture analysis of engineering materials. It is usually used in crack problems [25–27], composite laminates [28, 29], and bone cement analysis [30]. Therefore, computation of hip fracture risk index (FRI) over a region of interest (ROI) based on the strain energy criterion theoretically should be more accurate for assessing hip fracture risk. To the best of our knowledge, there are currently no published studies that use the strain energy criterion for the hip fracture risk assessment. The objective of this study is to improve the hip fracture risk assessment procedure previously developed by Luo et al. [11] by introducing the strain energy criterion.

2. Materials and Methods

The proposed methodology for assessment of hip fracture risk in the critical regions of femur using the strain energy criterion determined from QCT-based finite element model is shown in Figure 1. The procedure is explained in detail in the following.

2.1. QCT-Based Finite Element Model

2.1.1. QCT-Scan of Femur

The purpose of this study is to accurately assess hip fracture risk, so, a 3D finite element model of subject’s femur is required to achieve it. The 3D model can be constructed from the subject’s femur QCT images. QCT slices are produced using multiple scanners with a set of proper acquisition and reconstruction parameters (Figure 2(a)). Slice thickness of 1 mm is commonly used. The scanned QCT images are stored in the format of Digital Imaging and Communications in Medicine (DICOM), which can be used for the construction of a 3D FE model. A proper segmentation is done to separate the femur for constructing the 3D model. Each voxel in the QCT-scan has an intensity (or grey scale) that is expressed as Hounsfield Unit (HU), which is correlated to bone density [31, 32]. Threshold value for the cancellous bone is from 100 HU to 200 HU and for the cortical bone is between 200 HU and 2000 HU. QCT images of 60 clinical cases (30 females and 30 males) were acquired from the Winnipeg Health Science Centre in an anonymous way under a human research ethics approval. The subjects are in the age scope from 50 to 82 years (average of 65 years). Statistical information of the clinical cases is listed in Table 1.

(a)

(b)

(c)

(d)

(e)

(f)

2.1.2. Generation of Finite Element Mesh

In the first step, the geometrical model of the femur is generated from clinical QCT images using Mimics (Materialise, Leuven, Belgium). QCT images (in DICOM format) are imported to Mimics for segmentation (Figure 2(a)) and construction of 3D geometric model of the femur (Figure 2(b)). With the 3D geometric model, a FE mesh is generated using the 3-matic module in Mimics (Figure 2(c)). The 4-node linear tetrahedral element SOLID72 in ANSYS was used in this study. To investigate model convergence, FE models with different maximum element edge lengths were created. For each FE model, displacement was calculated under the same loading and boundary conditions. The maximum element edge length that produced converged finite element solutions was obtained and used in all the rest of FE simulations. The number of elements for each case is assigned based on the maximum element edge length that provided converged FE solution, meaning that the number of elements is different from case to case.

2.1.3. Assignment of Material Properties

To construct a more faithful FE model, bone material properties are considered inhomogeneous and isotropic in this study. Information on the inhomogeneous isotropic mechanical properties of the bone can be derived from the CT data using a mathematical relationship between the CT numbers and the mechanical properties of bone. The following empirical equation was used to determine bone ash density () from HU number [33, 34]:

Equations (2) and (3), derived by Keller [35], were, respectively, used to assign Young’s modulus () and yield stress () according to the bone ash density:

A constant Poisson’s ratio () was considered as suggested in the literature [13, 36, 37]. To assign material properties, elements are grouped into a number of discrete material bins using Mimics (Materialise, Leuven, Belgium), to approximately represent the continuous distribution of the inhomogeneous bone mechanical properties. To determine the maximum number of material bins, convergence study was performed. Models with different material bins were created for convergence study. For each FE model, displacement was calculated under the same loading and boundary conditions. The maximum number of material bins that generated converged finite element solutions was obtained. Figure 2(d) shows the isotropic inhomogeneous distribution of material properties.

2.2. Finite Element Analysis Using ANSYS

The finite element model of femur with the assigned material properties was output from Mimics and then imported to ANSYS for finite element analysis. For a precise assessment of hip fracture risk during the single-leg stance and the sideways fall, loading and boundary conditions simulating the single-leg stance and the sideways fall configurations are required in the FE model. To simulate the single-leg stance configuration, 2.5 times the patient’s body weight was applied as a distributed load on the femoral head in direction of femoral shaft axis [38] and femur was fixed at the distal end [13, 39]; see Figure 2(e). Considerwhere is the subject’s body weight in Newton (N). To simulate sideways fall, the distal end of femur was completely fixed and the surface of femoral head was fixed in the loading direction (Figure 2(f)) [40, 41]. The impact force during the sideways fall acting on the greater trochanter perpendicularly (Figure 2(f)) is given by [38, 42]:where is the height of the subject in centimeter (cm). Loading and boundary conditions on the greater trochanter, the femoral head, and the distal end of femur were applied to a group of nodes using APDL codes (Figures 2(e) and 2(f)). After importing the QCT-based FE model and applying the loading and boundary conditions, finite element analysis was performed and finite element solutions were obtained. In all the analysis, the nodal displacements, stresses, and strains were obtained for each subject.

2.3. Detection of the Three Critical Cross-Sections on the Femur

Hip fractures usually occur at one of the anatomical locations: the femoral neck, the intertrochanter, and the subtrochanter as illustrated in Figure 3. According to clinical observations, 49 percent of hip fractures are intertrochanteric, 37 percent are at femoral neck, and 14 percent are subtrochanteric [43]. Therefore, the smallest femoral neck cross-section (SFN CS), the intertrochanteric cross-section (IntT CS), and the subtrochanteric cross-section (SubT CS) are the three critical cross-sections of femur that usually have the highest fracture risk (Figure 3). To determine the smallest femoral neck cross-section and the intertrochanteric cross-section, neck-shaft angle is needed. The neck-shaft angle is the angle between the femoral neck axis and the femoral shaft axis. This angle traditionally is measured on conventional radiography images, or using 2D images projected from CT/MRI data. In spite of their popularity, these methods are based on oversimplification of the real 3D anatomy and may lead to large errors due to the inaccuracies in selection of the measurement plane [44–46]. In this study, the neck-shaft angle was measured using a 3D measurement technique based on fitting functions. In this technique, the shapes of particular parts of the femur are approximated using geometric entities such as circle, cylinder, and sphere, which are well fitted to the actual anatomy, and the geometrical relationships between these entities are obtained to estimate the neck-shaft angle.

First, a sphere was fitted to the femoral head, to obtain the position of the joint’s centre of rotation, which is also the femoral head centre. Then, the femoral neck axis and the femoral shaft axis are identified by applying the “fit ruled surface direction” function on the femoral neck and shaft. All fitting functions were applied using the 3-matic module in Mimics. The neck-shaft angle was also measured by 3-matic module of Mimics (Figure 4). With the femoral neck-shaft angle, the intertrochanteric cross-section and the smallest femoral neck cross-section were found using in-house computer codes. The smallest femoral neck cross-section is chosen as the cross-section with the smallest area in the neck region and the intertrochanteric cross-section is chosen as the cross-section that has the largest area in the intertrochanteric region [47]. By using APDL codes, perpendicular planes on the femoral neck axis were determined and then areas of the cross-sections were calculated. The planes with the smallest and the largest areas were chosen, respectively, as the smallest femoral neck cross-section and the intertrochanteric cross-section. The subtrochanteric cross-section is considered five centimeters below the lesser trochanter [48] (Figure 3).

2.4. Hip Fracture Risk Index Definition Using the Strain Energy Criterion

Based on the previous discussion of the bone failure mechanism and microstructure, the strain energy criterion is theoretically more favourable for hip fracture risk assessment. The strain energy at the three critical cross-sections of femur induced by the applied forces was computed using in-house developed MATLAB codes and the data extracted by APDL codes from the obtained finite element solutions. The plane boundaries of the three critical cross-sections, extracted from the finite element mesh, were imported to MATALB to generate a 2D mesh for calculating the cross-section strain energy. Figure 5 shows the generated triangle elements over the smallest femoral neck cross-section, the intertrochanteric cross-section, and the subtrochanteric cross-section.

(a)

(b)

(c)

The strain energy at the three critical cross-sections induced by the applied forces is the sum of strain energy in all the triangle elements; that is, where is the strain energy in Element induced by the applied forces and is the number of triangle elements created over the concerned cross-section. Gaussian integration method was used to calculate the strain energy in elements. Integration points in each triangle element were determined using in-house MATLAB codes. By using Gaussian integration method, the strain energy of Element induced by the applied forces is calculated aswhere is the strain energy density of Element ; is the strain energy density at the integration point of in Element ; is the weight at the integration point; is the determinant of the Jacobean matrix of the triangle element; and is the number of integration points over the triangle element. The strain energy density at an integration point () was determined from the finite element solutions obtained by the QCT-based FE model; that is, where and . The strain energy density at each integration point can be expressed by the finite element solutions aswhere is the displacement vector consisting of displacements at the nodes of Element ; matrix consists of the derivatives of shape functions of the element; and is the material property matrix. Considerwhere Poisson’s ratio is constant () and Young’s modulus is a function of the bone density as given in (2). For each integration point, its Young’s modulus is determined by the bone density at the point.

The maximum allowable strain energy over a critical cross-section of the femur was also computed from the obtained finite element solutions using in-house MATLAB codes. The maximum allowable strain energy (or the yield strain energy) over a critical cross-section is obtained aswhere is the yield strain energy in Element and is the number of triangle elements over the concerned cross-section. The Gaussian integration method was also used to calculate the maximum allowable strain energy in each triangle element. The maximum allowable strain energy that a triangle element () can sustain is given bywhere is the yield strain energy density in Element ; is the number of integration points; and is the yield strain energy density at integration point and is calculated aswhere and are, respectively, the yield stress and Young’s modulus at the integration point. Both of them are function of bone density, as given in (2) and (3).

Hip fracture risk index () at a critical cross-section is defined aswhere and are, respectively, obtained from (6) and (11).

3. Results

3.1. Convergence Studies

3.1.1. Element Size in Femur Finite Element Analysis

The convergence of finite element solutions in a representative case is shown in Figure 6. The convergence study showed that the finite element displacements converged with the maximum element edge length smaller than 8 mm. Therefore, in construction of the rest of femur FE models, the maximum element edge length was set to 8 mm.

3.1.2. Assignment of Inhomogeneous Material Properties

For convergence study in assigning the inhomogeneous material properties, 3D femur FE models with different material bins were created. For each FE model with different material bins, the maximum displacement at the smallest femoral neck cross-section was monitored under the same loading and boundary conditions. Data on the displacement were compared among the FE models with different material bins. The results of the convergence study showed that the displacement did not change significantly with the number of material bins larger than 50. Therefore, in the assignment of material properties for all the cases, 50 discrete material bins were considered.

3.1.3. Element Size in Calculating Cross-Section Strain Energy

Convergence study was also performed to determine the element size used in integrating cross-section strain energy, as it affects the calculated fracture risk index (FRI). The FRI at the concerned cross-section was calculated with different maximum element edge lengths. The results are plotted in Figure 7. The FRI did not change significantly with the maximum element edge length smaller than 5 mm. Therefore, the maximum element edge length was set to 5 mm in calculating cross-section strain energy.

3.1.4. The Number of Integration Points in Calculating Cross-Section Strain Energy

The effect of the number of integration points on the calculated FRI was investigated. FRI at the smallest femoral neck cross-section was computed for 5 clinical cases with different number of integration points. The relative errors between FRIs obtained with 3 and 7 integration points are shown in Table 2. As it can be seen, the errors are not significant. Therefore, the 3-point integration rule was used in this study to reduce computational time.

3.2. Stress and Strain Patterns at the Three Critical Cross-Sections

For the 10 clinical cases (5 females and 5 males, totally 20 right and left femurs), the maximum von Mises stress and strain at the three critical cross-sections of femur during both the single-leg stance and the sideways fall are shown in Figures 8 and 9. It can be observed that, during the sideways fall, the femoral neck and the intertrochanteric region received higher stresses than the subtrochanteric region (Table 4). But during the single-leg stance, the patterns in the stresses are different (Table 3); first, the differences between the stresses over the three regions are much smaller; for some cases, the stresses at the subtrochanteric region are higher than those in the other two regions (Figure 8). Strains in the three regions have similar patterns during both the single-leg stance and the sideways fall (Tables 5 and 6).

3.3. Comparison of Hip Fracture Risk at the Three Critical Cross-Sections

For the 60 clinical cases (30 females and 30 males), hip fracture risk indices based on the strain energy criterion were calculated for the smallest femoral neck, the intertrochanteric, and the subtrochanteric cross-section of femur during the single-leg stance and the sideways fall. The calculated fracture risk indices are shown in Figures 10 and 11.

As shown in Tables 7, 8, 9, and 10 and Figures 12 and 13, the average FRI at the smallest femoral neck was higher than those at the intertrochanteric and the subtrochanteric cross-section during both the single-leg stance and sideways fall.

For the single-leg stance, FRIs for all cases at the three critical cross-sections are much lower than one (FRI ≪ 1), indicating that the possibility of hip fracture incidence in the single-leg stance was low. For the sideways fall, FRIs of 8 right femurs and 7 left femurs at the smallest femoral neck cross-section and FRIs of 7 right femurs and 5 left femurs at the intertrochanteric cross-section were higher than one (FRI > 1), meaning that there is possibility for fracture occurring in these regions; but the FRIs at the subtrochanteric cross-section in all cases were much lower than one (FRI ≪ 1), indicating that there is lower possibility of fracture in this region. Figure 14 shows the number of possible fractures at the three critical cross-sections of femur during the sideways fall, that is, the cases that have FRI larger than one (FRI > 1) at one of the three critical cross-sections.

3.4. Comparison of Hip Fracture Risk in Women and Men

In this study, hip fracture risk for the 30 females and 30 males was assessed. The average hip fracture risk of females was generally higher than that of males. As it can be seen from Tables 11–13 and Figures 15 and 16, the average FRI at the smallest femoral neck cross-section, the intertrochanteric cross-section, and the subtrochanteric cross-section of the 30 females is higher than that of the 30 males for both the single-leg stance and the sideways fall. Figure 17 shows the number of possible hip fractures in women and men during the sideways fall for the studied cases.

4. Discussion

Hip fracture may occur anywhere from the articular cartilage of the hip joint to the femur shaft [48]. Not only the location of the fracture types but also the etiology differs. It was reported that women with the intertrochanteric fracture have significantly lower BMD than those with the femoral neck fracture [49–51]. On the other hand, the femoral neck fracture may not be mainly attributed to low BMD but may be related to external causes such as sideways fall [52]. Femoral neck and intertrochanteric fractures are often the result of falls from standing height and impact onto the greater trochanter, particularly for the elderly. The subtrochanteric fractures, on the other hand, are typically the result of high energy impacts such as motor vehicle accidents and falls from a height [53].

The selection of bone failure criterion is challenging. In literatures, stress and strain based failure criteria such as the von Mises stress and strain criteria and the maximum principle stress and strain criteria were commonly used to assess hip fracture risk. To the best of our knowledge, the strain energy based failure criterion has not been used yet for hip fracture risk assessment. Whereas the cancellous bone failure is in the form of buckling and deformation (strain intensity) and the cortical bone failure is related to its local cracking (stress intensity), strain energy failure criterion, which is a combination of both stress and strain intensities, is theoretically more reasonable than other failure criteria for hip fracture risk assessment. The differences between the strains in the three critical regions of femur during both the single-leg stance and the sideways fall (Figure 9, Tables 5 and 6) are much higher than the differences between the corresponding stresses (Figure 8, Tables 3 and 4), indicating that bone failure is more sensitive to the strains because of its fragility property and the effects of strains should also be considered in bone fracture risk assessment.

Results of this study show that the femoral neck and the intertrochanteric region have higher fracture risk than the subtrochanteric region (Figures 12 and 13 and Tables 7–10), which is consistent with the fact that the femoral neck and the intertrochanteric region have a larger proportion of cancellous bone than the subtrochanteric region; and the cancellous bone is generally weaker than the cortical bone. Therefore, hip fracture is most likely to initiate first at the femoral neck and then in the intertrochanteric region or in the subtrochanteric region. For all subjects, there is very low hip fracture risk during the single-leg stance (FRI < 1). Among all the cases, during the sideways fall, 15 femurs at the femoral neck and 12 femurs at the intertrochanteric region have FRI higher than one (FRI > 1), while there is very low fracture risk at the subtrochanteric region for all the subjects (FRI < 1) (Figure 14). Our findings have a good agreement with the previous clinical observations. According to the clinical observations, the femoral neck is the most common location for a hip fracture, accounting for 45% to 53% of hip fractures; the intertrochanteric fractures account for approximately 38% to 50% of all hip fractures; and the subtrochanteric fractures are less common than the femoral neck and intertrochanteric fractures, accounting for approximately 5% to 15% of hip fractures [54, 55]. Also it was observed that, in old people with age between 65 and 99, femoral neck and intertrochanteric fractures occurred with approximately the same frequency [56]. Data on 176 geriatric patients with hip fractures showed that 59% were the intertrochanteric fractures while the rest (41%) were the intracapsular neck fractures [57]. In another study by Michelson et al. [43], it was observed that 37% of hip fractures are in the femoral neck, 49% are intertrochanteric, and 14% are subtrochanteric. Therefore, the femoral neck and the intertrochanteric region are more disposed to fracture than the subtrochanteric region.

Both the fracture risk level and the potential fracture location are patient-dependent and depend on BMD and other conditions that have not been studied in this paper. We have found that women generally have lower bone strength than men and thus are exposed to higher hip fracture risk during both the single-leg stance and the sideways fall (Figures 15 and 16 and Tables 11–13). Based on the calculated FRI, in the 30 women potentially there are 12 femoral neck fractures and 7 intertrochanteric fractures, while for the 30 men there are only 3 potential femoral neck and 5 potential intertrochanteric fractures (Figure 17). Our finding is consistent with other studies that have found that hip fracture risk in women is generally higher than in men. The results of a worldwide study by Dhanwal et al. [58] on incidence and epidemiology of hip fractures in Asia, Africa, Europe, Latin America, North America, and Oceania show that women are more disposed to the hip fracture risk than men in different countries over the world. In the study by Jacobsen et al. [59], the age-adjusted hip fracture incidence rates in white US males and females were 4.3 and 8.1 per 1,000 per year, respectively. Higher hip fracture risk in women may be related to their lower BMD.

5. Conclusion

A reliable methodology to assess hip fracture risk in individuals is crucially important for preventing hip fracture and initiating a treatment. The purpose of this study is to propose a more effective hip fracture risk index that is based on the strain energy failure criterion, and it is able to better describe bone failure mechanism and microstructure. The proposed fracture risk index can predict not only the fracture risk level, but also the potential fracture location. The results of this study showed that there is a very low hip fracture risk during the single-leg stance, while, during the sideways fall, there is a high fracture risk at the femoral neck and the intertrochanteric region, compared to the subtrochanteric region. Based on the results obtained from this study, women are more prone to hip fracture than men. The procedure described in this study can be implemented into computer programming and used in hip fracture prevention and monitoring of osteoporosis treatments in the elderly. The method may also help design more effective hip protectors by providing feedback information such as stresses/strains in the hip for adjusting the design parameters. The main limitation of this study is that no experiment has been conducted to validate the predicted fracture risk levels and potential fracture locations, which has been set as future work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

N. M. Resnick and S. L. Greenspan, “‘Senile’ osteoporosis reconsidered,” Journal of the American Medical Association, vol. 261, no. 7, pp. 1025–1029, 1989.
View at: Publisher Site | Google Scholar
B. Gullberg, O. Johnell, and J. A. Kanis, “World-wide projections for hip fracture,” Osteoporosis International, vol. 7, no. 5, pp. 407–413, 1997.
View at: Publisher Site | Google Scholar
O. Johnell, “The socioeconomic burden of fractures: today and in the 21st century,” The American Journal of Medicine, vol. 103, no. 2, pp. 20S–26S, 1997.
View at: Google Scholar
J. S. Bauer and T. M. Link, “Advances in osteoporosis imaging,” European Journal of Radiology, vol. 71, no. 3, pp. 440–449, 2009.
View at: Publisher Site | Google Scholar
D. M. Albertsson, D. Mellström, C. Petersson, and R. Eggertsen, “Validation of a 4-item score predicting hip fracture and mortality risk among elderly women,” Annals of Family Medicine, vol. 5, no. 1, pp. 48–56, 2007.
View at: Publisher Site | Google Scholar
D. Marshall, O. Johnell, and H. Wedel, “Meta-analysis of how well measures of bone mineral density predict occurrence of osteoporotic fractures,” British Medical Journal, vol. 312, no. 7041, pp. 1254–1259, 1996.
View at: Publisher Site | Google Scholar
J. A. Kanis, H. Johansson, A. Oden, and E. V. McCloskey, “Assessment of fracture risk,” European Journal of Radiology, vol. 71, no. 3, pp. 392–397, 2009.
View at: Publisher Site | Google Scholar
S. Kaptoge, L. I. Benevolenskaya, A. K. Bhalla et al., “Low BMD is less predictive than reported falls for future limb fractures in women across Europe: results from the European Prospective Osteoporosis Study,” Bone, vol. 36, no. 3, pp. 387–398, 2005.
View at: Publisher Site | Google Scholar
T. L. N. Järvinen, H. Sievänen, K. M. Khan, A. Heinonen, and P. Kannus, “Shifting the focus in fracture prevention from osteoporosis to falls,” British Medical Journal, vol. 336, no. 7636, pp. 124–126, 2008.
View at: Publisher Site | Google Scholar
M. R. McClung, “To FRAX or not to FRAX,” Journal of Bone and Mineral Research, vol. 27, no. 6, pp. 1240–1242, 2012.
View at: Publisher Site | Google Scholar
Y. Luo, Z. Ferdous, and W. D. Leslie, “A preliminary dual-energy X-ray absorptiometry-based finite element model for assessing osteoporotic hip fracture risk,” Proceedings of the Institution of Mechanical Engineers Part H: Journal of Engineering in Medicine, vol. 225, no. 12, pp. 1188–1195, 2011.
View at: Publisher Site | Google Scholar
Y. Luo, Z. Ferdous, and W. D. Leslie, “Precision study of DXA-based patient-specific finite element modeling for assessing hip fracture risk,” International Journal for Numerical Methods in Biomedical Engineering, vol. 29, no. 5, pp. 615–629, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
J. H. Keyak, S. A. Rossi, K. A. Jones, and H. B. Skinner, “Prediction of femoral fracture load using automated finite element modeling,” Journal of Biomechanics, vol. 31, no. 2, pp. 125–133, 1997.
View at: Publisher Site | Google Scholar
J. C. Lotz, E. J. Cheal, and W. C. Hayes, “Fracture prediction for the proximal femur using finite element models: part I—linear analysis,” Journal of Biomechanical Engineering, vol. 113, no. 4, pp. 353–360, 1991.
View at: Publisher Site | Google Scholar
J. C. Lotz, E. J. Cheal, and W. C. Hayes, “Fracture prediction for the proximal femur using finite element models: part II—nonlinear analysis,” Journal of Biomechanical Engineering, vol. 113, no. 4, pp. 361–365, 1991.
View at: Publisher Site | Google Scholar
D. Testi, M. Viceconti, F. Baruffaldi, and A. Cappello, “Risk of fracture in elderly patients: a new predictive index based on bone mineral density and finite element analysis,” Computer Methods and Programs in Biomedicine, vol. 60, no. 1, pp. 23–33, 1999.
View at: Publisher Site | Google Scholar
T. Ota, I. Yamamoto, and R. Morita, “Fracture simulation of the femoral bone using the finite-element method: how a fracture initiates and proceeds,” Journal of Bone and Mineral Metabolism, vol. 17, no. 2, pp. 108–112, 1999.
View at: Publisher Site | Google Scholar
E. Schileo, F. Taddei, L. Cristofolini, and M. Viceconti, “Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro,” Journal of Biomechanics, vol. 41, no. 2, pp. 356–367, 2008.
View at: Publisher Site | Google Scholar
H. Gong, M. Zhang, Y. Fan, W. L. Kwok, and P. C. Leung, “Relationships between femoral strength evaluated by nonlinear finite element analysis and BMD, material distribution and geometric morphology,” Annals of Biomedical Engineering, vol. 40, no. 7, pp. 1575–1585, 2012.
View at: Publisher Site | Google Scholar
J. H. Keyak and S. A. Rossi, “Prediction of femoral fracture load using finite element models: an examination of stress- and strain-based failure theories,” Journal of Biomechanics, vol. 33, no. 2, pp. 209–214, 2000.
View at: Publisher Site | Google Scholar
M. Mirzaei, M. Keshavarzian, and V. Naeini, “Analysis of strength and failure pattern of human proximal femur using quantitative computed tomography (QCT)-based finite element method,” Bone, vol. 64, pp. 108–114, 2014.
View at: Publisher Site | Google Scholar
J. S. Stölken and J. H. Kinney, “On the importance of geometric nonlinearity in finite-element simulations of trabecular bone failure,” Bone, vol. 33, no. 4, pp. 494–504, 2003.
View at: Publisher Site | Google Scholar
J. Cordey and E. Gautier, “Strain gauges used in the mechanical testing of bones. Part I: theoretical and technical aspects,” Injury, vol. 30, no. 1, pp. SA7–SA13, 1999.
View at: Google Scholar
M. Mirzaei, A. Zeinali, A. Razmjoo, and M. Nazemi, “On prediction of the strength levels and failure patterns of human vertebrae using quantitative computed tomography (QCT)-based finite element method,” Journal of Biomechanics, vol. 42, no. 11, pp. 1584–1591, 2009.
View at: Publisher Site | Google Scholar
G. C. Sih, “Strain-energy-density factor applied to mixed mode crack problems,” International Journal of Fracture, vol. 10, no. 3, pp. 305–321, 1974.
View at: Publisher Site | Google Scholar
G. C. Sih, Mechanics of Fracture Initiation and Propagation, Springer, Dordrecht, The Netherlands, 1991.
View at: Publisher Site
Y. Wei, “An extended strain energy density failure criterion by differentiating volumetric and distortional deformation,” International Journal of Solids and Structures, vol. 49, no. 9, pp. 1117–1126, 2012.
View at: Google Scholar
W. E. Wolfe and T. S. Butalia, “A strain-energy based failure criterion for non-linear analysis of composite laminates subjected to biaxial loading,” Composites Science and Technology, vol. 58, no. 7, pp. 1107–1124, 1998.
View at: Publisher Site | Google Scholar
T. S. Butalia and W. E. Wolfe, “A strain-energy-based non-linear failure criterion: comparison of numerical predictions and experimental observations for symmetric composite laminates,” Composites Science and Technology, vol. 62, no. 12-13, pp. 1697–1710, 2002.
View at: Publisher Site | Google Scholar
R. Vasu, D. R. Carter, and W. H. Harris, “Evaluation of bone cement failure criteria with applications to the acetabular region,” Journal of Biomechanical Engineering, vol. 105, no. 4, pp. 332–337, 1983.
View at: Publisher Site | Google Scholar
J. H. Keyak, J. M. Meagher, H. B. Skinner, and C. D. Mote Jr., “Automated three-dimensional finite element modelling of bone: a new method,” Journal of Biomedical Engineering, vol. 12, no. 5, pp. 389–397, 1990.
View at: Publisher Site | Google Scholar
T. M. Keaveny, R. E. Borchers, L. J. Gibson, and W. C. Hayes, “Trabecular bone modulus and strength can depend on specimen geometry,” Journal of Biomechanics, vol. 26, no. 8, pp. 991–1000, 1993.
View at: Publisher Site | Google Scholar
C. M. Les, J. H. Keyak, S. M. Stover, K. T. Taylor, and A. J. Kaneps, “Estimation of material properties in the equine metacarpus with use of quantitative computed tomography,” Journal of Orthopaedic Research, vol. 12, no. 6, pp. 822–833, 1994.
View at: Publisher Site | Google Scholar
D. Dragomir-Daescu, J. O. D. Buijs, S. McEligot et al., “Robust QCT/FEA models of proximal femur stiffness and fracture load during a sideways fall on the hip,” Annals of Biomedical Engineering, vol. 39, no. 2, pp. 742–755, 2011.
View at: Publisher Site | Google Scholar
T. S. Keller, “Predicting the compressive mechanical behavior of bone,” Journal of Biomechanics, vol. 27, no. 9, pp. 1159–1168, 1994.
View at: Publisher Site | Google Scholar
D. T. Reilly and A. H. Burstein, “The elastic and ultimate properties of compact bone tissue,” Journal of Biomechanics, vol. 8, no. 6, pp. 393–405, 1975.
View at: Publisher Site | Google Scholar
W. C. van Buskirk and R. B. Ashman, “The elastic moduli of bone,” Transaction of American Society of Mechanical Engineers, vol. 45, pp. 131–143, 1981.
View at: Google Scholar
T. Yoshikawa, C. H. Turner, M. Peacock et al., “Geometric structure of the femoral neck measured using dual-energy x-ray absorptiometry,” Journal of Bone and Mineral Research, vol. 9, no. 7, pp. 1053–1064, 1994.
View at: Google Scholar
M. Bessho, I. Ohnishi, T. Matsumoto et al., “Prediction of proximal femur strength using a CT-based nonlinear finite element method: differences in predicted fracture load and site with changing load and boundary conditions,” Bone, vol. 45, no. 2, pp. 226–231, 2009.
View at: Publisher Site | Google Scholar
J. E. M. Koivumäki, J. Thevenot, P. Pulkkinen et al., “Ct-based finite element models can be used to estimate experimentally measured failure loads in the proximal femur,” Bone, vol. 50, no. 4, pp. 824–829, 2012.
View at: Publisher Site | Google Scholar
K. K. Nishiyama, S. Gilchrist, P. Guy, P. Cripton, and S. K. Boyd, “Proximal femur bone strength estimated by a computationally fast finite element analysis in a sideways fall configuration,” Journal of Biomechanics, vol. 46, no. 7, pp. 1231–1236, 2013.
View at: Publisher Site | Google Scholar
S. N. Robinovitch, W. C. Hayes, and T. A. McMahon, “Prediction of femoral impact forces in falls on the hip,” Journal of Biomechanical Engineering, vol. 113, no. 4, pp. 366–374, 1991.
View at: Publisher Site | Google Scholar
J. D. Michelson, A. Myers, R. Jinnah, Q. Cox, and M. Van Natta, “Epidemiology of hip fractures among the elderly: risk factors for fracture type,” Clinical Orthopaedics and Related Research, no. 311, pp. 129–135, 1995.
View at: Google Scholar
J. S. Kim, T. S. Park, S. B. Park, J. S. Kim, I. Y. Kim, and S. I. Kim, “Measurement of femoral neck anteversion in 3D. Part 1: 3D imaging method,” Medical and Biological Engineering and Computing, vol. 38, no. 6, pp. 603–609, 2000.
View at: Publisher Site | Google Scholar
B. Atilla, A. Oznur, O. Caǧlar, M. Tokgözoǧlu, and M. Alpaslan, “Osteometry of the femora in Turkish individuals: a morphometric study in 114 cadaveric femora as an anatomic basis of femoral component design,” Acta orthopaedica et Traumatologica Turcica, vol. 41, no. 1, pp. 64–68, 2007.
View at: Google Scholar
E. Sariali, A. Mouttet, G. Pasquier, and E. Durante, “Three-dimensional hip anatomy in osteoarthritis: analysis of the femoral offset,” Journal of Arthroplasty, vol. 24, no. 6, pp. 990–997, 2009.
View at: Publisher Site | Google Scholar
R. Nikander, P. Kannus, P. Dastidar et al., “Targeted exercises against hip fragility,” Osteoporosis International, vol. 20, no. 8, pp. 1321–1328, 2009.
View at: Publisher Site | Google Scholar
B. Abrahamsen, T. van Staa, R. Ariely, M. Olson, and C. Cooper, “Excess mortality following hip fracture: a systematic epidemiological review,” Osteoporosis International, vol. 20, no. 10, pp. 1633–1650, 2009.
View at: Publisher Site | Google Scholar
E. Vega, C. Mautalen, H. Gómez, A. Garrido, L. Melo, and A. O. Sahores, “Bone mineral density in patients with cervical and trochanteric fractures of the proximal femur,” Osteoporosis International, vol. 1, no. 2, pp. 81–86, 1991.
View at: Publisher Site | Google Scholar
N. Nakamura, T. Kyou, K. Takaoka, K. Ohzono, and K. Ono, “Bone mineral density in the proximal femur and hip fracture type in the elderly,” Journal of Bone and Mineral Research, vol. 7, no. 7, pp. 755–759, 1992.
View at: Google Scholar
S. L. Greenspan, E. R. Myers, L. A. Maitland, T. H. Kido, M. B. Krasnow, and W. C. Hayes, “Trochanteric bone mineral density is associated with type of hip fracture in the elderly,” Journal of Bone and Mineral Research, vol. 9, no. 12, pp. 1889–1894, 1994.
View at: Google Scholar
C. A. Mautalen, E. M. Vega, and T. A. Einhorn, “Are the etiologies of cervical and trochanteric hip fractures different?” Bone, vol. 18, no. 3, supplement, pp. 133S–137S, 1996.
View at: Google Scholar
R. B. Birrer, Field Guide to Fracture Management, Lippincott Williams & Wilkins, 2005.
K. J. Koval and J. D. Zuckerman, “Hip fractures: I. Overview and evaluation and treatment of femoral-neck fractures,” Journal of the American Academy of Orthopaedic Surgeons, vol. 2, no. 3, pp. 141–149, 1994.
View at: Google Scholar
Fixing Hip Fractures, March 2015, http://www.hopkinsmedicine.org/gec/series/fixing_hip_fractures.
M. R. Karagas, G. L. Lu-Yao, J. A. Barrett, M. L. Beach, and J. A. Baron, “Heterogeneity of hip fracture: age, race, sex, and geographic patterns of femoral neck and trochanteric fractures among the US elderly,” American Journal of Epidemiology, vol. 143, no. 7, pp. 677–682, 1996.
View at: Publisher Site | Google Scholar
V. Trikha and S. Rastogi, “Epidemiology and rehabilitation of hip fractures in the geriatric population,” Indian Journal of Physical Medicine and Rehabilitation, vol. 16, no. 1, pp. 16–19, 2005.
View at: Google Scholar
D. K. Dhanwal, E. M. Dennison, N. C. Harvey, and C. Cooper, “Epidemiology of hip fracture: worldwide geographic variation,” Indian Journal of Orthopaedics, vol. 45, no. 1, pp. 15–22, 2011.
View at: Publisher Site | Google Scholar
S. J. Jacobsen, J. Goldberg, T. P. Miles, J. A. Brody, W. Stiers, and A. A. Rimm, “Hip fracture incidence among the old and very old: a population-based study of 745,435 cases,” The American Journal of Public Health, vol. 80, no. 7, pp. 871–873, 1990.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2015 Hossein Kheirollahi and Yunhua Luo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2208

Downloads

1186

Citations