Home | About | Journals | Submit | Contact Us | Français |

**|**Biomed Eng Online**|**v.11; 2012**|**PMC3549900

Formats

Article sections

- Abstract
- Background
- Methods
- Results
- Discussion
- Conclusions
- Abbreviations
- Competing interests
- Authors' contributions
- References

Authors

Related links

Biomed Eng Online. 2012; 11: 84.

Published online 2012 November 14. doi: 10.1186/1475-925X-11-84

PMCID: PMC3549900

Enrique López,^{1} Elena Ibarz,^{2} Antonio Herrera,^{3,}^{4,}^{5} Jesús Mateo,^{3,}^{4,}^{5} Antonio Lobo-Escolar,^{3,}^{4,}^{5} Sergio Puértolas,^{2} and Luis Gracia^{}^{2,}^{6}

Enrique López: se.razinu@czepole; Elena Ibarz: se.razinu@zrabie; Antonio Herrera: se.nogara.dulas@arerreha; Jesús Mateo: moc.ono@aoetamj; Antonio Lobo-Escolar: se.oohay@ralocseobola; Sergio Puértolas: se.razinu@bps; Luis Gracia: se.razinu@ivargul

Received 2012 August 31; Accepted 2012 November 7.

Copyright ©2012 López et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Osteoporotic hip fractures represent major cause of disability, loss of quality of life and even mortality among the elderly population. Decisions on drug therapy are based on the assessment of risk factors for fracture, from BMD measurements. The combination of biomechanical models with clinical studies could better estimate bone strength and supporting the specialists in their decision.

A model to assess the probability of fracture, based on the Damage and Fracture Mechanics has been developed, evaluating the mechanical magnitudes involved in the fracture process from clinical BMD measurements. The model is intended for simulating the degenerative process in the skeleton, with the consequent lost of bone mass and hence the decrease of its mechanical resistance which enables the fracture due to different traumatisms. Clinical studies were chosen, both in non-treatment conditions and receiving drug therapy, and fitted to specific patients according their actual BMD measures. The predictive model is applied in a FE simulation of the proximal femur. The fracture zone would be determined according loading scenario (sideway fall, impact, accidental loads, etc.), using the mechanical properties of bone obtained from the evolutionary model corresponding to the considered time.

BMD evolution in untreated patients and in those under different treatments was analyzed. Evolutionary curves of fracture probability were obtained from the evolution of mechanical damage. The evolutionary curve of the untreated group of patients presented a marked increase of the fracture probability, while the curves of patients under drug treatment showed variable decreased risks, depending on the therapy type.

The FE model allowed to obtain detailed maps of damage and fracture probability, identifying high-risk local zones at femoral neck and intertrochanteric and subtrochanteric areas, which are the typical locations of osteoporotic hip fractures.

The developed model is suitable for being used in individualized cases. The model might better identify at-risk individuals in early stages of osteoporosis and might be helpful for treatment decisions.

Currently, osteoporotic fractures represent major cause of disability, loss of quality of life and even death among the elderly population [1]. Osteoporosis is caused by a skeletal involution linked to aging, which is more prevalent in women: the lifetime risk for a fragility fracture at the age of 50 lies within the range of 40% in women [2].

Hip fracture is considered the most devastating osteoporotic fracture. It generally occurs around age 80 [3] and affects women more often than men in a ratio 3:1 [4]. The worldwide incidence of hip fractures in 2000 was estimated at 1.6 million [5] and, if present population forecasts are borne out, there will be 6.26 million in 2050 [6]. Hip fracture is now a major cause of morbidity, mortality and disability, and represents a significant economic cost [4,5,7,8].

Decisions on drug therapy in osteoporotic patients without previous fractures are mainly based on the analysis of risk factors which predispose to fracture. A risk assessment tool called FRAX® (Fracture Risk Assessment Tool) has been developed by the World Health Organization (WHO) for this purpose [9-12]. The risk of major osteoporotic fractures (hip, vertebrae, humerus and wrist), or specifically hip fracture over the next 10 years can be estimated with the FRAX® tool. The probability of fracture is calculated on the basis of age, body mass index and several dichotomized variables (previous fracture, smoking, rheumatoid arthritis, etc.). Optionally, bone mineral density (BMD) of the femoral neck can be included for risk calculation. Other studies have questioned the effectiveness of FRAX as a tool for predicting fracture risk [13,14].

Several previous surveys have assessed the risk of fracture using various methodologies, but mostly based on BMD measurements [15,16]. BMD measurements have also been used for determining the mechanical strength [17] or to develop statistical models for predicting the risk of fracture [18]. Another study used morphological measurements of the femoral head and neck to determine the risk of fracture [19].

Regarding simulation by means of finite element method, both micro and macro-mechanical models have been suggested, with different characteristics and methodologies. These models can be used for prediction of bone strength at different ages, or to prediction of fracture risk. Thus, Lee [20] provides a micromechanical model of bone behavior under different densities, which is difficult to extrapolate to the scale required to get realistic predictions. Boccacio, Zhang [21,22] provide more advanced macro-mechanical models, which analyze a complete functional unit of the spine (two vertebrae with their intervertebral disk) in terms of mechanical behavior depending on bone density. Macneil [23] sets a 2D model in the sagittal plane (vertebrae L1-L4), using bone geometry and BMD measurements obtained from radiographs and DXA. In this model, developed for the analysis of vertebral fractures, the stiffness is calculated based on the patient's age, taking into account an exponential decline. In the case of proximal femoral fractures, Kaneko [24] develops a model based on imaging techniques (Quantitative Computed Tomography (QCT)), focused on the estimation of bone strength as a function of age in normal populations. The model aimed to establish a statistical correlation between the prediction of bone strength and the risk of osteoporotic fracture. A different methodology is used by Bryan [25], who suggests a parametric model incorporating both the geometry and the properties of bone, which allows a range of results reflecting the statistical variation of the model parameters. A similar methodology is used by Bessho [26], although the parametric analysis, in this case, refers to load and support conditions of the model. Some authors have begun to incorporate yield criteria for fracture risk prediction. In this respect, Derikx [27] applies the Drucker-Prager criterion on a model made from QCT, with asymmetric yielding in tension and compression. Similarly Tellache [28] applies an anisotropic yield criterion on a model constructed from imaging (CT scan) for prediction of fracture risk. With a different approach, Amin [29] performed a comparative analysis of fracture risk predictions based on BMD measurements against those ones obtained from an FE model developed from QCT, correlating bone strength with fracture risk. Finally, on the matter of drug treatments, Keaveny [30] discusses the influence on bone strength of PTH and alendronate, using a FE model developed from QCT scans of osteoporotic patients. Currently, the most popular tool to assess the fracture risk is FRAX, combining different clinical factors and BMD measurements. However, all the above methods have some limitations concerning their approaches, including clinical or mechanical magnitudes related to bone fracture in an independent way, but without consider their mutual influence as actually happens. So, clinical BMD values allow evaluating the mechanical properties of bone (stiffness and resistance), which condition the stresses patterns causing cumulative mechanical damage. That cumulative damage would cause fractures in different load scenarios.

In view of the current models limitations, we have developed a model for predicting the risk of osteoporotic hip fractures based on the Damage Mechanics and Fracture Mechanics. This model will be incorporated into a finite element code to simulate their evolution over time. Thus, we will estimate the probability of hip fracture in the mid and long terms. Bone mineral density (BMD) measurements, from dual energy X-ray absorptiometry (DXA) scans, will be incorporated into the model in order to fit it to clinical conditions. The model is not intended for simulating the bone fracture, but to predict the degenerative process in the skeleton, with the consequent lost of bone mass and hence the decrease of its mechanical resistance which enables the fracture due to different traumatisms.

In order to develop a predictive model which takes into account the mechanical parameters involved in the fracture process, a correlation between these magnitudes and those ones measured in clinical terms is firstly required. To this effect, Carter and Hayes [31] established a direct relationship between Young's modulus and bone density for low strain rates (0.01):

$$E=2875\phantom{\rule{0.12em}{0ex}}{\rho}^{3}$$

(1)

On the other hand, the relationship between the BMD value and the apparent density is adjusted, according to experimental results, as:

$$\rho ={\rho}_{\mathit{m\xe1x}}{\left(\frac{BMD}{BM{D}_{\mathit{\text{m\xe1x}}}}\right)}^{\lambda}$$

(2)

being λ a parameter that depends on the sample data. For the present study a value of 9/25 was assigned for λ, in order to fit the actual data presented in [32].

BMD is the current standard for diagnosis of osteoporosis. Over 100 worldwide published papers assessing hip BMD evolution, both in natural conditions and in patients under drug therapy were selected for analysis. Among them, three treatments have been selected for the comparative study: alendronate 10 mg per day [33], oral ibandronate 2.5 mg per day [34] and PTH 1–84 100 mg per day [35]. These drugs have proven to be effective, and their BMD evolution curves show the required regularity in the analyzed time period.

Bisphosphonates settle in the bone tissue and its effect persists during some time after its administration, and hence it is completely accepted in the clinical practice the use of intermittent or discontinuous treatments. In most of the studies a maximum five years follow-up is done, although a treatment continuous for ten years is also accepted. Moreover, there are several studies concerning the safety of long term treatments (ten years) [36-42]. With respect to the use of combined therapies (teriparatide and bisphosphonates), it is a usual clinical practice in the osteoporosis treatment [43-47].

Regarding the natural history of BMD, the average curve published by Mazess [48] was chosen as a reference. Briefly, the natural curve of BMD linked with age was compared to the curves of BMD in patients treated with bisphosphonates or derivatives of parathyroid hormone.

Since the standard adjustment techniques do not provide enough accuracy and reliability to be applied to the predictive model, higher level, more complex, adjustments have been set out to obtain continuous curves of regression. These make it possible to obtain continuous curves for the BMD evolution. The evolution trends of the different selected cases in the study can be extrapolated from those continuous curves, ensuring a consistent behavior. The following adjustments were proposed:

Polynomial:

$$\rho ={a}_{n}{t}^{n}+{a}_{n-1}{t}^{n-1}+\dots +{a}_{1}t+{a}_{0}$$

(3)

Exponential:

$$\rho =k\left(1-{e}^{-\left({a}_{n}{t}^{n}+{a}_{n-1}{t}^{n-1}+\dots +{a}_{1}t+{a}_{0}\right)}\right)$$

(4)

Exponential asymptotic:

$$\rho ={\rho}_{0}+\left({\rho}_{\mathit{l\xedm}}-{\rho}_{0}\right)\left(1-{e}^{-\frac{at}{{t}_{m}-t}}\right)$$

(5)

These adjustments (Equations 3, 4 and 5) have been applied to the four examined curves (natural evolution and the three therapies). Several choices have been made during this adjustment: the lowest mean square error, the closest to unity correlation coefficient R^{2}, and a 10-year standardized follow-up period with the possibility of extrapolation to longer periods up to 15 years (Figure (Figure1).1). That extrapolation is intended to be consistent in mathematical terms, independently of the actual duration of the considered study. This condition will permit in the future an easy incorporation of more long term studies when available. Standardization of BMD measurements included in each published paper was required in order to make comparisons among them [49], because de data were obtained from different densitometric equipment depending on the study (Hologic, Lunar and Norland).

The Karganovin^{′}s Damage Mechanics model [50] was selected for the simulation of the degenerative process. This model defines the mechanical damage, *D*, as a function of the equivalent strain, *ε*_{c}:

$$D=1-\kappa {\left({\epsilon}_{c}\right)}^{\gamma}$$

(6)

being:

$${\epsilon}_{c}=\sqrt{{\epsilon}_{I}^{2}+{\epsilon}_{\mathit{II}}^{2}+{\epsilon}_{\mathit{III}}^{2}-{\epsilon}_{I}{\epsilon}_{\mathit{II}}-{\epsilon}_{I}{\epsilon}_{\mathit{III}}-{\epsilon}_{\mathit{II}}{\epsilon}_{\mathit{III}}}$$

(7)

the equivalent strain, where *ε*_{I}, *ε*_{II}, *ε*_{III} are the principal components of the strain tensor. This mechanical damage leads to a decrease of the bone stiffness, according with the relationship:

$$D=1-\sqrt{\frac{E}{{E}_{0}}}$$

(8)

where *E*_{0} corresponds to the Young’s modulus of healthy bone and *E* is the actual value of Young’s modulus of the bone with cumulative mechanical damage. The *κ* and *γ* constants depend on the critical damage value, the critical strain and the strain threshold, *ε*_{0}, below which no damage is produced (*ε*_{c}=*ε*_{0} implies *D*=*0*):

$$\begin{array}{cccc}\hfill \gamma =\frac{ln\left(\frac{1}{1-{D}_{\mathit{cri}}}\right)}{ln\left(\frac{{\epsilon}_{0}}{{\epsilon}_{\mathit{cri}}}\right)}\hfill & \hfill \hfill & \hfill \hfill & \hfill \kappa =\frac{1}{{\left({\epsilon}_{0}\right)}^{\gamma}}\hfill \end{array}$$

(9)

where *ε*_{cri} corresponds to the equivalent strain value that produces the critical damage, *D*_{cri}. Once the damage model is defined, a relationship between the damage level and the probability of fracture must be set. Obviously several factors are involved in the equation and may difficult the calculations: damage location, amount of damage, range and type of load cycles, and so on. However, since osteoporosis is a generalized disease affecting bone mass extensively, and it occurs in older people with little variation in their life habits, a simplified model could be used. In this work, a model of fracture probability based on the law of Paris [51], which explains the stable crack growth under monotonic loading, was developed. The number of loading cycles needed to increase damage from the *D*_{i} value to the *D*_{f} value was expressed as follows:

$$N=\frac{1}{\alpha \phantom{\rule{0.12em}{0ex}}{\left(\Delta \sigma \right)}^{\beta}{\gamma}^{\beta /2}\left(\frac{\beta}{2}-1\right)}\left[{\left(\frac{1}{\omega {D}_{i}}\right)}^{\frac{\beta}{2}-1}-{\left(\frac{1}{\omega {D}_{f}}\right)}^{\frac{\beta}{2}-1}\right]$$

(10)

where *α* and *β* are material parameters associated with the law of Paris, *γ* is a parameter related to the stress intensity factor defined in Linear Elastic Fracture Mechanics, which depends on geometric aspects and stress distribution, Δ*σ* is the variation of stress in each loading cycle and *ω* is a parameter which relates the mechanical damage to the crack size *a* (*a*=*ωD*).

By normalizing the probability of fracture, assigning value 1 to critical damage (*D*=*D*_{cri})_{,} and value 0 to no damage (*D*=*0*), we finally obtain the probability of fracture as a function of both the number of cycles and the damage:

$$P\left(fracture\right)=1-\frac{N\left(D\right)}{{N}_{max}}={\left(\frac{D}{{D}_{\mathit{cr\xed}}}\right)}^{\frac{\beta}{2}-1}$$

(11)

In Equation (11), the magnitude *N*_{max} represents the maximum number of load cycles necessary for a mechanical damage equivalent to the critical value, *D*_{cri}. A value *5* has been given to the *β* coefficient for cortical bone, according to Taylor [52], and, in accordance with Kargarnovin [50], a critical damage, *D*_{cri}, of *0*.*38* and a critical strain, *ε*_{cri}, of *0*.*0174* have been considered, with a strain threshold, *ε*_{0}, of *0*.*0015* (no damage is produced by strains below *0*.*0015*).

Progression curves are basic references and provide only general information. In order to apply the model to specific patients we must consider, in addition to the trend, the reference density value of the patient. According to the law of interpolated natural evolution, the density matching with the age of the patient is given by:

$${\rho}_{N}\left({t}_{0}\right)\ne {\rho}_{0}$$

(12)

being *ρ*_{N}(*t*_{0}) the average density measured in the patient and *ρ*_{0} the matching value from the reference progression curve, that is, the actual patient’s density doesn’t necessarily be equal to the value corresponding to the average curve for the considered population. There is an offset that should be added to the reference progression curve of density, providing a translation of the average curve, allowing an adaptation for each individual patient, so:

$${\rho}_{p}^{N}\left(t\right)={\rho}_{N}\left(t\right)+\left[{\rho}_{0}-{\rho}_{N}\left({t}_{0}\right)\right]$$

(13)

When a drug therapy is applied to the same patient, a similar correction is required because a new offset arises. So, the progression curve for this patient under a treatment would be:

$${\rho}_{p}^{T}\left(t\right)={\rho}_{T}\left(t\right)+\left[{\rho}_{0}-{\rho}_{T}\left({t}_{0}\right)\right]$$

(14)

In the Equations (12) to (14), the subscripts or superscripts *N* and *T* represent natural evolution or evolution with treatment, respectively. All these adjustments provide the estimated BMD value for any type of patient, at any age, and under any prescribed therapy. From this value, mechanical properties of bone can be calculated. It must be noticed that the considered curves represent the mean evolutionary curves for the population, and an individual patient could not follow the curve exactly but in an approximate way.

Densitometric data of the healthy femur have been taken as the starting point for this study, based on a previous work of our group on the biomechanical behavior of a femoral stem [32]. In that survey, densitometric data were correlated both with the apparent volumetric density, and with the Young moduli of each of the Gruen zones. Table Table11 depicts the BMD, the apparent volumetric density, and the Young moduli correlation corresponding to a healthy femur.

From all the previous calculations, an evolutionary algorithm has been implemented (Figure (Figure2),2), which has been used combined with a finite element model of the femur [32] using the Abaqus software [53]. The procedure was based on a finite element model of the upper half of the femur, made up of tetrahedral elements with quadratic approximation (C3D10 in Abaqus nomenclature), as shown in Figure Figure3A3A and b. To generate the model a 3D laser scanner Roland PIZCA was used for the outer geometry and thirty transverse direction tomographic cross-sections and eight longitudinal direction cross-sections were taken using CAT (General Electric Brightspeed Elite) to determine the geometry of the cancellous bone. Moreover a 3D reconstruction was made to obtain the shape of the medullar cavity [32].

3D FE model of proximal femur: a) Adaptation of Gruen zones to the healthy femur; b) Boundary conditions on the proximal femur.

In the final mesh, consisting of 408094 elements and 75223 nodes, cortical bone, cancellous bone and bone marrow have been differentiated (229931 elements for cortical bone, 166220 elements for cancellous bone and 11943 elements for bone marrow). To guarantee the accuracy in the FE results, a sensitivity analysis was performed with a mesh refinement in order to achieve a convergence towards a minimum of the potential energy, with a tolerance of 1% between consecutive meshes.

In addition, the femur has been divided into 7 zones of Gruen [54], since the available information on BMD was related to these zones (Figure (Figure3A).3A). The mechanical properties obtained from the evolutionary algorithm have been assigned to each of the Gruen zones. For this purpose, the apparent density is calculated for each Gruen zone from its BMD measure according Equation (2), and then, the Young’ modulus is obtained from Equation (1). The same process was followed for cancellous bone. Any case, cortical bone stresses are slightly influenced by the mechanical properties of cancellous bone and bone marrow since their stiffness is much lower than the corresponding to cortical bone. The function of cancellous bone, in mechanical terms, consists in providing stability to the thin walled cortical bone. A linear elastic behavior has been defined for all the materials. The final collection of results was focused on the proximal femur, mainly on the area located within the edge of the femoral head and the subtrochanteric region.

As a boundary condition, the middle third of the femoral diaphysis was clamped (Figure (Figure3B),3B), since this area is far enough from the proximal femur to avoid significant perturbations in the stress distribution. Thus, computational cost was reduced if compared with the entire femur model.

In the finite element model construction, as important as the reaction strength on the femoral head due to the body weight, is the inclusion of the muscle forces to be considered in the simulation. In our model only the abductor muscle forces were included, in line with several authors [55,56]. As a rule of thumb, abductor muscles cause a reaction on the femoral head of 2.75 times the body weight. However, load increases as much as 4 times the body weight when the heel impacts the ground, and during the double support stage of the gait [57]. The latter situation was considered, as the worst case, in order to set the boundary conditions. According to data from the study [32] a 79.3 kg body-weight was set as a reference. So, two load conditions were imposed (Figure (Figure3B3B):

• Reaction strength on the femoral head due to the body weight (3110 N).

• Load due to the abductor muscles, applied to the proximal area of the greater trochanter (1360 N).

The above values are slightly higher than those included in recent studies [58,59]. This model has been used in predicting the evolution of femoral fracture probability, by comparing the natural history and the expected evolution under different therapies.

Firstly, adjustment models were applied both to BMD physiological curve (weighted average from [48]) and the curves of patients under three different therapies: alendronate (10 mg per day), oral ibandronate (2.5 mg per day) and PTH 1–84 (100 mg per day). Figure Figure11 shows the different evolutionary curves obtained in each case, together with the interpolated curves. Despite an extreme irregularity of some values, the correlation coefficients have been 0.992, 0.980, 0.991 and 0.982, respectively.

Comparisons among various treatments can be drawn regarding different simulations. Figure Figure4A4A shows BMD evolution curves under natural conditions and under three therapies (alendronate, ibandronate, and PTH 1–84 plus oral alendronate). As can be seen, patients treated with alendronate showed a significant initial increase in bone density, which remained stable until the end of the study period. Oral ibandronate also leads to a remarkable increase in BMD during the first stage; a progressive decline after treatment interruption (similar to the natural physiological curve), and a new increase with the therapy reintroduction. In opposition, the physiological curve was characterized by a substantial and progressive decline in BMD linked to an increased risk of fracture.

Evolution of different magnitudes at femoral neck: a) BMD; b) Young’s modulus; c) Equivalent strain; d) Mechanical damage; e) Evolution of the fracture probability at femoral neck in different conditions; f) Evolution of the fracture probability **...**

Equivalent comparisons can be made with different parameters. Thus, Figure Figure4B4B shows the evolution of Young's modulus in untreated patients and in those under drug therapy (alendronate, ibandronate, and PTH 1–84 plus oral alendronate). Similar patterns to the observed in BMD curves were obtained, according with the direct relationship between the two parameters. Figure Figure4C4C shows the evolution of the equivalent strain in subjects under natural conditions and under the same three therapies (alendronate, ibandronate, and PTH 1–84 plus oral alendronate). This equivalent deformation was estimated as a weighted average of the considered bone area. Contrary to previous trends, as bone stiffness decreases the average strain showed a significant increase (as it occurs in the natural history of BMD). In the case of therapies which preserve BMD, and therefore bone stiffness, the average strain value remained stable or even decreases. Finally, Figure Figure4D4D illustrates the evolution of mechanical damage under natural conditions and under the three proposed therapies (alendronate, ibandronate, and PTH 1–84 plus oral alendronate). Similar trends to those obtained for strain can be seen. In fact, when the average strains value increases, mechanical damage increases too.

As a final result, the evolutionary curves of fracture probability were obtained from the evolution of mechanical damage. The estimated probability, according to mechanical damage caused by strains, is calculated for the initial patient’s state (Figure (Figure4E).4E). Evolutionary curves of fracture probability increase can be obtained by referencing all results to the initial state (Figure (Figure4F).4F). As can be seen, the fracture probability showed a marked increase in the natural evolution curve, while the curves of the treated patients showed lower degrees of risk, depending on the therapy type.

In addition to previous results, programmed subroutines make it possible to obtain damage and fracture probability maps and to identify high-risk zones of the femoral bone (Figure (Figure5A5A and b). Femoral neck and intertrochanteric and subtrochanteric areas are the zones at greatest risk, in coincidence with the typical locations of osteoporotic hip fractures. It must be pointed out that Figure Figure55 doesn’t represent actual fracture zones, but zones with poorer bone mechanical strength due to cumulative damage. The actual fracture zone would be determined according loading scenario (sideway fall, impact, accidental loads, etc.).

A novel method for estimating the risk in osteoporotic patients has been developed. Clinical data (DXA measures) and mechanical magnitudes related to bone strength were combined in this tool. The mechanical properties of bone are updated from BMD values obtained from clinical data of untreated patients and in those under different treatments. The model uses Damage and Fracture Mechanics concepts to evaluate the fracture probability in an evolutionary algorithm.

The model can be used in a personalized way from BMD measurements in each case. The model can contribute to the development of diagnostic tools for detection of early stages of osteoporosis. It may also be helpful for treatment decisions in selected patients. Many studies have been carried out, both in the clinical [9-19] and the simulation fields [20-30]. But a simple and reliable model, useful as a tool for diagnosis and prevention in our daily practice, has not yet been achieved.

Several predictive models can be found in the literature, but statistical models are currently the most reliable [9-12], regardless of mechanical issues involving bone strength estimated for different conditions and ages. The development of new techniques for measuring BMD has focused much of the recent research in the clinical setting, but the mechanical aspects have not been adequately studied [13-16]. In other cases, the improvements have been applied on statistical models previously [18,19]. Only in one published paper [17] a new method derived from DXA measures was developed for bone strength assessment.

Concerning the finite element simulation, and based on previous micromechanical models [20], various methods have been developed, but they are difficult to extrapolate to the scale required to get realistic predictions from different approaches [21,22]. The incorporation of the latest imaging techniques (QCT) and BMD measurement (DXA) [24,27-29] has allowed the improvement of these models. However, most models use standard yield criteria for estimating the risk of fracture [27,28], without considering essential aspects in fracture analysis. Nonetheless, all models assume that bone mineral density is the basic measurement, and it should therefore be used as a benchmark in predicting fracture risk.

Recent works [59] establish a correlation between BMD at femoral neck and fracture risk considering FE analysis and experimental data. Their obtain a von Mises stress distribution with maximum values at femoral neck and subtrochanteric area, just the same zones where the present model predicts the maximum fracture probability.

From the mechanical point of view, the exposition of the bone to cyclic loads of high value in a damaged bone, once the degenerative process is started, decreases its strength over the time and produces a cumulative damage which can lead to a final fracture. It seems apparent that Damage Mechanics and Fracture Mechanics criteria should be incorporated in any model intending to obtain reliable results. In this regard, our model combines all these requirements, and might be useful as a basis for future more sophisticated models.

Moreover, this model enables to incorporate future developments with the same methodology. In the first term, a more accurate bone density distribution, that is, not by Gruen zones but by individual elements in the mesh, could be used. That requires a planned collection of BMD data, by means of DXA or CT scan images. More complex damage models can be added, including mechanical behavior of anisotropic or mixed models, based on both strains and stresses. It would also be possible to include crack growth models fitting to the results of in vitro bone fracture. Finally, a parametric finite element model of the femoral head could be performed including both the loads produced on the bone and the shape and dimensions of a specific patient.

Despite DXA measurements just quantify bone mass and not bone quality, it is widely accepted as a macroscopic indicator of bone strength and stiffness and also that micro-fractures exert an important influence on the mechanical strength of the bone.

Finally, clinical trials are needed to validate the proposed model in order to apply it to the clinical practice helping for treatment decisions.

A mechanical model based on Damage and Fracture Mechanics and DXA measurements, for predicting the probability of fracture in osteoporotic patients has been carried out. The model represents a first step towards the development of new tools for diagnosis and prevention of osteoporosis. The incorporation of clinical measurements and simulation results will be useful for an individualized monitoring and treatment in specific patients.

BMD: Bone mineral density; PTH: Parathyroid hormone; FE: Finite elements; FRAX: Fracture risk assessment tool; WHO: World health organization; DXA: Dual-emission x-ray absorptiometry; NHANES: National health and nutrition examination survey; QCT: Quantitative computed tomography; CT: Computed tomography.

The authors declare that they have no conflict of interest.

AH and LG conceived the approach of this work. EL and LG conceived and developed the predictive model for fracture probability. AH, ALE and JM contribute with clinical measurements and experience with osteoporotic patients. EL, EI and SP conceived and developed the finite element model and carried out all the simulations. AH and LG coordinated the work between surgeons and engineers. All authors participated in the drawing up of the manuscript, and read and approved the final manuscript.

This work has been partially financed through the research project “Development and Clinical Application of a Mechanical Model for the Analysis of Spine Fractures due to Low Energy Trauma, based on Bone Mineral Density Measurements”, supported by Foundation Mutua Madrileña.

- National Osteoporosis Foundation. Clinician’s Guide to Prevention and Treatment of Osteoporosis. Washington, DC: National Osteoporosis Foundation; 2010. ISBN: 978-0-9798989-9-0.
- Johnell O, Kanis J. Epidemiology of osteoporotic fractures. Osteoporos Int. 2005;16(Suppl 2):S3–S7. doi: 10.1007/s00198-004-1702-6. [PubMed] [Cross Ref]
- Instituto de Información Sanitaria. Estadísticas Comentadas: La Atención a La Fractura De Cadera En Los Hospitales Del SNS (on line) Ministerio de Sanidad y Política Social, Madrid; 2010. http:es/estadEstudios/estadisticas/cmbdhome.htm. Accessed 7 September 2010.
- Herrera A, Martínez AA, Ferrández L, Moreno A. Epidemiology of osteoporotic hip fractures in Spain. Int Orthop. 2006;30(1):11–14. doi: 10.1007/s00264-005-0026-2. [PMC free article] [PubMed] [Cross Ref]
- Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int. 2006;17(12):1726–1733. doi: 10.1007/s00198-006-0172-4. [PubMed] [Cross Ref]
- Kannus P, Parkkari J, Sievanen H, Heinonen A, Vuori I, Järvinen M. Epidemiology of hip fractures. Bone. 1996;18(Suppl 1):S57–S63. doi: 10.1016/8756-3282(95)00381-9. [PubMed] [Cross Ref]
- Bessette L, Jean S, Lapointe-Garant MP, El B, Davison KS, Ste-Marie LG, Brown JP. Direct medical costs attributable to peripheral fractures in Canadian post-menopausal women. Osteoporos Int. 2011. (Epub ahead of print) [PubMed] [Cross Ref]
- Bass E, French DD, Bradham DD, Rubestein LZ. Risk-adjusted mortality rates of elderly veterans with hip fractures. Ann Epidemiol. 2007;17(7):514–519. doi: 10.1016/j.annepidem.2006.12.004. [PubMed] [Cross Ref]
- Kanis JA, Borgstrom F, De Laet C, Johansson H, Johnell O, Jonsson B, Oden A, Zethraeus N, Pfleger B, Khaltaev N. Assessment of fracture risk. Osteoporos Int. 2005;16(6):581–589. doi: 10.1007/s00198-004-1780-5. [PubMed] [Cross Ref]
- Kanis JA, Black D, Cooper C, Dargent P, Dawson-Hughes B, De Laet C, Delmas P, Eisman J, Johnell O, Johnsson B, Melton L, Oden A, Papapoulos S, Pols H, Rizzoli R, Silman A, Tenenhouse A. A new approach to the development of assessment guidelines for osteoporosis. Osteoporos Int. 2002;13(7):527–536. doi: 10.1007/s001980200069. [PubMed] [Cross Ref]
- Kanis JA, Oden A, Johnell O, De Laet C, Brown J, Burckhardt P, Cooper C, Christiensen C, Cummings S, Eisman JA, Fujiwara S, Glüer C, Goltzman D, Hans D, Krieg MA, La Croix A, McCloskey E, Mellstrom D, Melton LJ, Pols H, Reeve J, Sanders K, Schott AM, Silman A, Torgerson D, van Staa T, Watts NB, Yoshimura N. The use of clinical risk factors enhances the performance of BMD in the prediction of hip and osteoporotic fractures in men and women. Osteoporos Int. 2007;18(8):1033–1046. doi: 10.1007/s00198-007-0343-y. [PubMed] [Cross Ref]
- Kanis JA, Johnell O, Oden AJH, McCloskey E. FRAX and the assessment of fracture probablility in men and women from the UK. Osteoporos Int. 2008;19(4):385–397. doi: 10.1007/s00198-007-0543-5. [PMC free article] [PubMed] [Cross Ref]
- Van Geel TA, van den Bergh JP, Dinant GJ, Geusens PP. Individualizing fracture risk prediction. Maturitas. 2010;65(2):143–148. doi: 10.1016/j.maturitas.2009.12.007. [PubMed] [Cross Ref]
- Ensrud KE, Lui LY, Taylor BC, Schousboe JT, Donaldson MG, Fink HA, Cauley JA, Hillier TA, Browner WS, Cummings SR. A comparison of prediction models for fractures in older women: is more better? Arch Intern Med. 2009;169(22):2087–2094. doi: 10.1001/archinternmed.2009.404. [PMC free article] [PubMed] [Cross Ref]
- Moayyeri A, Kaptoge S, Dalzell N, Bingham S, Luben RN, Wareham NJ, Reeve J, Khaw KT. Is QUS or DXA better for predicting the 10-year absolute risk of fracture? J Bone Miner Res. 2009;24(7):1319–1325. doi: 10.1359/jbmr.090212. [PubMed] [Cross Ref]
- Langsetmo L, Leslie WD, Zhou W, Goltzman D, Kovacs CS, Prior J, Josse R, Olszynski WP, Davison KS, Anastassiades T, Towheed T, Hanley DA, Kaiser S, Kreiger N. Using the same bone density reference database for men and women provides a simpler estimation of fracture risk. J Bone Miner Res. 2010;25(10):2108–2114. doi: 10.1002/jbmr.112. [PMC free article] [PubMed] [Cross Ref]
- Boehm HF, Horng A, Notohamiprodjo M, Eckstein F, Burklein D, Panteleon A, Lutz J, Reiser M. Prediction of the fracture load of whole proximal femur specimens by topological analysis of the mineral distribution in DXA-scan images. Bone. 2008;43(5):826–831. doi: 10.1016/j.bone.2008.07.244. [PubMed] [Cross Ref]
- Schechner Z, Luo G, Kaufman JJ, Siffert RS. A Poisson process model for hip fracture risk. Med Biol Eng Comput. 2010;48(8):799–810. doi: 10.1007/s11517-010-0638-6. [PubMed] [Cross Ref]
- Baker-LePain JC, Luker KR, Lynch JA, Parimi N, Nevitt MC, Lane NE. Active shape modeling of the hip in the prediction of incident hip fracture. J Bone Miner Res. 2011;26(3):468–474. doi: 10.1002/jbmr.254. [PMC free article] [PubMed] [Cross Ref]
- Lee TC, McHugh PE, O'Brien FJ, O’Mahoney D, Taylor D, Bruzzi M, Rackard SM, Kennedy OD, Mahony NJ, Harrison N, Lohfield S, Brennan O, Gleeson J, Hazenberg JG, Mullins L, Tyndyk M, McNamara LM, O’Kelly KU, Prendergast PJ. In: Topics in Bio-Mechanical Engineering. Prendergast PJ, McHugh PE, editor. Dublin: Trinity Centre for Bio-Engineering & National Centre for Biomedical Engineering Science; 2004. Bone for life: osteoporosis, bone remodelling and computer simulation; pp. 58–93.
- Boccaccio A, Vena P, Gastaldi D, Franzoso G, Pietrabissa R, Pappalettere C. Finite element analysis of cancellous bone failure in the vertebral body of healthy and osteoporotic subjects. Proc Inst Mech Eng H. 2008;222(7):1023–1036. doi: 10.1243/09544119JEIM296. [PubMed] [Cross Ref]
- Zhang L, Yang G, Wu L, Yu B. The biomechanical effects of osteoporosis vertebral augmentation with cancellous bone granules or bone cement on treated and adjacent non-treated vertebral bodies: a finite element evaluation. Clin Biomech. 2010;25(2):166–172. doi: 10.1016/j.clinbiomech.2009.10.006. [PubMed] [Cross Ref]
- Macneil JA, Adachi JD, Goltzman D, Josse RG, Kovacs CS, Prior JC, Olszynski W, Davison KS, Kaiser SM. Predicting fracture using 2D finite element modelling. Med Eng Phys. 2011. (Epub ahead of print) [PMC free article] [PubMed] [Cross Ref]
- Kaneko M, Ohnishi I, Bessho M, Matsumoto T, Ohashi S, Tobita K, Nakamura K. Prediction of proximal femur strength by a quantitative computed tomography-based finite element method – creation of predicted strength data of the proximal femur according to age range in a normal population and analysis of risk factors for hip fracture. J Bone Joint Surg Br. 2011;93-B no(SUPP II):216. [PubMed]
- Bryan R, Nair PB, Taylor M. Use of a statistical model of the whole femur in a large scale, multi-model study of femoral neck fracture risk. J Biomech. 2009;42(13):2171–2176. doi: 10.1016/j.jbiomech.2009.05.038. [PubMed] [Cross Ref]
- Bessho M, Ohnishi I, Matsumoto T, Ohashi S, Matsuyama J, Tobita K, Kaneko M, Nakamura K. 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. 2009;45(2):226–231. doi: 10.1016/j.bone.2009.04.241. [PubMed] [Cross Ref]
- Derikx LC, Vis R, Meinders T, Verdonschot N, Tanck E. Implementation of asymmetric yielding in case-specific finite element models improves the prediction of femoral fractures. Comput Methods Biomech Biomed Engin. 2011;14(2):183–193. doi: 10.1080/10255842.2010.542463. [PubMed] [Cross Ref]
- Tellache M, Pithioux M, Chabrand P, Hochard C. Femoral neck fracture prediction by anisotropic yield criteria. Eur J Comp Mech. 2009;18(1):33–41.
- Amin S, Kopperdhal DL, Melton LJ 3rd, Achenbach SJ, Therneau TM, Riggs BL, Keaveny TM, Khosla S. Association of hip strength estimates by finite‐element analysis with fractures in women and men. J Bone Miner Res. 2011;26(7):1593–1600. doi: 10.1002/jbmr.347. [PMC free article] [PubMed] [Cross Ref]
- Keaveny TM, Hoffmann PF, Singh M, Palermo L, Bilezikian JP, Greenspan SL, Balck DM. Femoral Bone strength and its relation to cortical and trabecular changes after treatment with PTH, Alendronate, and their combination as assessed by finite element analysis of quantitative CT scans. J Bone Miner Res. 2008;23(12):1974–1982. doi: 10.1359/jbmr.080805. [PubMed] [Cross Ref]
- Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure. J Bone Joint Surg Am. 1977;59(7):954–962. [PubMed]
- Herrera A, Panisello JJ, Ibarz E, Cegoñino J, Puértolas JA, Gracia L. Long-term study of bone remodelling after femoral stem: a comparison between DEXA and finite element simulation. J Biomech. 2007;40(16):3615–3625. doi: 10.1016/j.jbiomech.2007.06.008. [PubMed] [Cross Ref]
- Bone HG, Hosking D, Devogelaer JP, Tucci JR, Emkey RD, Tonino RP, Rodríguez-Portales JA, Downs RW, Gupta J, Santora AC, Liberman UA. Ten years’ experience with alendronate for osteoporosis in postmenopausal women. N Engl J Med. 2004;350(12):1189–1199. doi: 10.1056/NEJMoa030897. [PubMed] [Cross Ref]
- Riis BJ, Ise J, von Stein T, Bagger Y, Christiansen C. Ibandronate: a comparison of oral daily dosing versus intermittent dosing in postmenopausal osteoporosis. J Bone Miner Res. 2001;16(10):1871–1878. doi: 10.1359/jbmr.2001.16.10.1871. [PubMed] [Cross Ref]
- Greenspan SL, Bone HG, Ettinger MP, Hanley DA, Lindsay R, Zanchetta JR, Blosch CM, Mathisen AL, Morris SA, Marriott TB. Effect of recombinant human parathyroid hormone (1–84) on vertebral fracture and bone mineral density in postmenopausal women with osteoporosis: a randomized trial. Ann Intern Med. 2007;146(5):326–339. [PubMed]
- Watts NB, Diab DL. Long-term use of bisphosphonates in osteoporosis. J Clin Endocrinol Metab Apr. 2010;95(4):1555–1565. doi: 10.1210/jc.2009-1947. [PubMed] [Cross Ref]
- Felsenberg D, Czerwinski E, Stakkestad J, Christiansen C, Civitelli R, Drezner MK. et al. Efficacy of monthly oral ibandronate is maintained over 5 years: the MOBILE LTE study. Osteoporos Int. 2009;20(Suppl. 1):S15. Abstract 0C32.
- Di Munno O, Delle Sedie A. Efficacy of ibandronate: a long term confirmation. Clin Cases Miner Bone Metab Jan. 2010;7(1):23–26. [PMC free article] [PubMed]
- Pazianas M, Cooper C, Ebetino FH, Russell RG. Long-term treatment with bisphosphonates and their safety in postmenopausal osteoporosis. Ther Clin Risk Manag. 2010;6:325–343. PMID: 20668715. [PMC free article] [PubMed]
- Rossini M, Idolazzi L, Adami S. Evidence of sustained vertebral and nonvertebral antifracture efficacy with ibandronate therapy: a systematic review. Ther Adv Musculoskelet Dis Apr. 2011;3(2):67–79. doi: 10.1177/ 1759720X10395651. [PMC free article] [PubMed] [Cross Ref]
- Pazianas M, Abrahamsen B. Safety of bisphosphonates. Bone Jul. 2011;49(1):103–110. doi: 10.1016/j.bone.2011.01.003. [PubMed] [Cross Ref]
- Bianchi G, Czerwinski E, Kenwright A, Burdeska A, Recker RR, Felsenberg D. Long-term administration of quarterly IV ibandronate is effective and well tolerated in postmenopausal osteoporosis: 5-year data from the DIVA study long-term extensión. Osteoporos Int. 2012;23:1769–1778. doi: 10.1007/s00198-011-1793-9. [PubMed] [Cross Ref]
- Kanis JA, Burlet N, Cooper CP, Delmas D, Reginster JY, Borgstrom F, Rizzoli R. European guidance for the diagnosis and management of osteoporosis in postmenopausal women. Osteoporos Int. 2008;19:399–428. doi: 10.1007/s00198-008-0560-z. [PMC free article] [PubMed] [Cross Ref]
- Compston J. The use of combination therapy in the treatment of postmenopausal osteoporosis. Endocrine. 2012;41:11–18. doi: 10.1007/s12020-011-9554-2. [PubMed] [Cross Ref]
- Schafer AL, Sellmeyer DE, Palermo L, Hietpas J, Eastell R, Shoback DM, Black DM. Six months of parathyroid hormone (1–84) administered concurrently versus sequentially with monthly Ibandronate over two years: the PTH and Ibandronate combination study (PICS) randomized trial. J Clin Endocrinol Metab Oct. 2012;97(10):3522–3529. doi: 10.1210/jc.2012-1844. [PubMed] [Cross Ref]
- Walker MD, Cusano NE, Sliney J, Jr, Romano M, Zhang C, McMahon DJ, Bilezikian JP. Combination therapy with risedronate and teriparatide in male osteoporosis. Endocrine. 2012. Oct 26. [Epub ahead of print]. PMID:23099796. [PubMed]
- Cusano NE, Bilezikian JP. Combination anabolic and antiresorptive therapy for osteoporosis. Endocrinol Metab Clin North Am Sep. 2012;41(3):643–654. doi: 10.1016/j.ecl.2012.04.005. Epub 2012 May 19. [PubMed] [Cross Ref]
- Mazess RB, Barden H. Bone density of the spine and femur in adult white females. Calcif Tissue Int. 1999;65(2):91–99. doi: 10.1007/s002239900663. [PubMed] [Cross Ref]
- Hanson J. Standardization of femur BMD. J Bone Miner Res. 1997;12(8):1316–1317. doi: 10.1359/jbmr.1997.12.8.1316. [PubMed] [Cross Ref]
- Kargarnovin MH, Bagher-Ebrahimi M, Katoozian HR. Damage Initiation and Growth in a Long Bone Under increasing Monotonic Loading Using the Continuum Damage Mechanics Principle. 2006. (Proc. Fifth International Conf on Eng Computational Technology, paper 191).
- Paris P, Erdogan F. A critical analysis of crack propagation laws. Jour Basic Engineering. 1963;85(4):528–534. doi: 10.1016/0753-3322(92)90002-O. [Cross Ref]
- Taylor D. Microcrack growth parameters for compact bone deduced from stiffness variations. J Biomech. 1998;31(7):587–592. doi: 10.1016/S0021-9290(98)00050-5. http://dx.doi.org/10.1016/S0021-9290(98)00050-5. [PubMed] [Cross Ref]
- Dassault Systèmes. [ http://www.3ds.com/]
- Gruen TA, McNeice GM, Amstutz HC. Modes of failure of cemented stem-type femoral components: a radiographic analysis of loosening. Clin Orthop. 1979;141:17–27. [PubMed]
- Kerner J, Huiskes R, Van Lenthe GH, Weinans H, van Rietbergen B, Engh CA, Amis AA. Correlation between pre-operative periprosthetic bone density and post-operative bone loss in THA can be explained by strain-adaptative remodelling. J Biomech. 1999;32(7):695–703. doi: 10.1016/S0021-9290(99)00041-X. [PubMed] [Cross Ref]
- Rietbergen B, Huiskes R. Load transfer and stress shielding of the hydroxyapatite ABG hip. A Study of stem length and proximal fixation. J Arthroplasty. 2001;16(8 suppl 1):55–63. doi: 10.1054/arth.2001.28369. [PubMed] [Cross Ref]
- Gutiérrez P, Doménech P, Roca J. In: Patología De La Cadera En El Adulto. López-Prats F, editor. Madrid: Sociedad Española de Cirugía Ortopédica y Traumatología; 2004. Biomecánica de la cadera (in Spanish) pp. 11–19.
- Sarikanat M, Yildiz H. Determination of bone density distribution in proximal femur by using the 3D orthotropic bone adaptation model. J of Engineering in Medicine. 2011;225:365–375.
- Tsouknidas A, Anagnostidis K, Maliaris G, Michailidis N. Fracture risk in the femoral hip region: a finite element analysis supported experimental approach. J Biomech. 2012;45:1959–1964. doi: 10.1016/j.jbiomech.2012.05.011. [PubMed] [Cross Ref]

Articles from BioMedical Engineering OnLine are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |