Bone creep-fatigue damage accumulation

Creep and fatigue tests were performed on human femoral cortical bone and the results were compared to a cumulative damage model for bone fracture. Fatigue tests in tension, compression, and reversed loading with a tensile mean stress were conducted at 2 Hz and 0.02 Hz. Load frequency had a strong influence on the number of cycles to failure but did not influence the total time to failure. Bone displayed poor creep-fracture properties in both tension and compression. The fracture surfaces of the tensile creep specimens are distinctly different than those of the compressive specimens. The results suggest that tensile cyclic loading creates primarily time-dependent damage and compressive cyclic loading creates primarily cycle-dependent damage. However, data for load histories involving both tensile and compressive loading indicate lower time to failure than predicted by a simple summation of time-dependent and cycle-dependent damage.     

Development of a surrogate model based on patient weight,bone mass and geometry to predict femoral neck strains and fracture loads

Osteoporosis and related bone fractures are an increasing global burden in our ageing society. Areal bone mineral density assessed through dual energy X-ray absorptiometry (DEXA), the clinically accepted and most used method, is not sufficient to assess fracture risk individually. Finite element (FE) modelling has shown improvements in prediction of fracture risk, better than aBMD from DEXA, but is not practical for widespread clinical use. The aim of this study was to develop an adaptive neural network (ANN)-based surrogate model to predict femoral neck strains and fracture loads obtained from a previously developed population-based FE model. The surrogate model performance was assessed in simulating two loading conditions: the stance phase of gait and a fall.The surrogate model successfully predicted strains estimated by FE (r² = 0.90–0.98 for level gait load case, r² = 0.92–0.96 for the fall load case). Moreover, an ANN model based on three measurements obtainable in clinics (femoral neck length (level gait) or maximum femoral neck diameter (fall), femoral neck bone mass, body weight) was able to give reasonable predictions (r² = 0.84–0.94) for all of the strain metrics and the estimated femoral neck fracture load. Overall, the surrogate model has potential for clinical applications as they are based on simple measures of geometry and bone mass which can be derived from DEXA images, accurately predicting FE model outcomes, with advantages over FE models as they are quicker and easier to perform.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.     

Effect of boundary conditions,impact loading and hydraulic stiffening on femoral fracture strength

Patient specific quantitative CT (QCT) imaging data together with the finite element (FE) method may provide an accurate prediction of a patient's femoral strength and fracture risk. Although numerous FE models investigating femoral fracture strength have been published, there is little consent on the effect of boundary conditions, dynamic loading and hydraulic strengthening due to intra-medullary pressure on the predicted fracture strength. We developed a QCT-derived FE model of a proximal femur that included node-specific modulus assigned based on the local bone density. The effect of three commonly used boundary conditions published in literature were investigated by comparing the resulting strain field due to an applied fracture load. The models were also augmented with viscoelastic material properties and subject to a realistic impact load profile to determine the effect of dynamic loads on the strain field. Finally, the effect of hydraulic strengthening was investigated by including node specific permeability and performing a coupled pore diffusion and stress analysis of the FE model. Results showed that all boundary conditions yield the same strain field patterns, but peak strains were 22% lower and fracture load was 18% higher when loaded at the greater trochanter than when loaded at the femoral head. Comparison of the dynamic models showed that material viscoelasticity was important, but inertial effects (vibration and shock) were not. Finally, pore pressure changes did not cause significant hydraulic strengthening of bone under fall impact loading.     

Impact response and simulation of damaged ulna with internal fixation

The objectives of this work were to explore a methodology that combines static and dynamic finite element (FE) analysis, linear elastic fracture mechanics (LEFM) and experimental methods to investigate a worst-case scenario in which a previously damaged bone plate system is subjected to an impact load. Cadaver ulnas with and without midshaft dynamic compression plates are subjected to a static three-point bend test and loaded such that subcritical crack growth occurs as predicted by a hybrid method that couples LEFM and static FE. The plated and unplated bones are then unloaded and subsequently subjected to a midshaft transverse impact test. A dynamic strain-based FE model is also developed to model the midshaft transverse impact test. The average value of the impact energy required for failure was observed to be 10.53% greater for the plated set. There appears to be a trade-off between impact damage and impact resistance when ulnas are supported by fixation devices. Predictions from the dynamic FE model are shown to corroborate inferences from the experimental approach.     

Fracture mechanics of the femoral neck in a composite bone model: Effects of platen geometry

Load applicator (platen) geometry used for axial load to failure testing of the femoral neck varies between studies and the biomechanical consequences are unknown. The purpose of this study was to determine if load application with a flat versus a conical platen results in differing fracture mechanics. Femurs were aligned in 25° of adduction and an axial compressive force was applied to the femoral heads at a rate of 6 mm/min until failure. Load application with the conical platen resulted in an average ultimate failure load, stiffness, and energy to failure of 9067 N, 4033 N/mm, and 12.12 J, respectively. Load application with the flat platen resulted in a significant (p<0.05) reduction in ultimate failure load (7620 N) and stiffness (2924 N/mm). Energy to failure (12.30 J) was not significantly different (p=0.893). Different fracture patterns were observed for the two platens and the conical platen produced fractures more similar to clinical observations. Use of a flat platen underestimates the strength and stiffness of the femoral neck and inaccurately predicts the associated fracture pattern. These findings must be considered when interpreting the results of prior biomechanical studies on femoral neck fracture and for the development of future femoral neck fracture models.     

Mechanical behaviour of femoral bones in bending loading

A series of 33 human femoral bones have been subjected to a four point bending test at high strain rates. Two different failure modes were recognized. A Y shaped fracture at the middle region induced by a pure bending moment yielded a zone of non-linearity at the load vs deformation curve and a higher bending force, more deformation of the structure and higher strain energy to fracture compared with the less frequently occurring oblique fracture at the distal third of the structure resulting in a failure without a 'plastic' portion at the load-deformation curve. Estimated values of bending modulus and maximum bending moment based upon a simple uniform beam model showed high correlation coefficients with the experimentally determined values. Scanning electron microscopic examination of the Y fracture showed distortion and void formation of the material at the structural level. This could explain the extensive non-elastic deformation prior to failure.     

Mechanobiological prediction of proximal femoral deformities in children with cerebral palsy

Children with cerebral palsy (CP) walk with altered gait and frequently exhibit proximal femoral deformities, such as anteversion and coxa valga. The objective of this research was to investigate the effect of specific gait patterns on the femoral morphology in CP.

In this study, the mechanobiological principles were implemented on a 3D finite element (FE) model of the proximal femur in order to predict changes in morphology over time in healthy and CP children. This model relies on the assumption that cyclic octahedral shear stress promotes growth and cyclic hydrostatic compressive stress inhibits growth. Growth was simulated over 16 iterations, representing approximately 5 months of growth.

The FE model predicts an increase in the femoral anteversion and coxa valga for CP loading conditions when compared with healthy ones. Understanding the role of loading in skeletal morphogenesis may help prevent bone deformities and improve function in children with gait abnormalities.     

Parametric finite element analysis of vertebral bodies affected by tumors

The vertebral column is the most frequent site of metastatic involvement of the skeleton. Due to the proximity to the spinal cord, from 5% to 10% of all cancer patients develop neurologic manifestations. As a consequence, fracture risk prediction has significant clinical importance. In this study, we model the metastatically involved vertebra so as to parametrically investigate the effects of tumor size, material properties and compressive loading rate on vertebral strength. A two-dimensional axisymmetric finite element model of a spinal motion segment consisting of the first lumbar vertebral body (no posterior elements) and adjacent intervertebral disc was developed to allow the inclusion of a centrally located tumor in the vertebral body. After evaluating elastic, mixed, and poroelastic formulations, we concluded that the poroelastic representation was most suitable for modeling the metastatically involved vertebra's response to compressive load. Maximum principal strains were used to localize regions of potential vertebral trabecular bone failure. Radial and axial vertebral body displacements were used as relative indicators of spinal canal encroachment and endplate failure. Increased tumor size and loading rate, and reduced trabecular bone density all elevated axial and radial displacements and maximum tensile strains. The results of this parametric study suggest that vertebral tumor size and bone density contribute significantly to a patients risk for vertebral fracture and should be incorporated in clinical assessment paradigms.     

During sideways falls proximal femur fractures initiate in the superolateral cortex: Evidence from high-speed video of simulated fractures

Results of recent imaging studies and theoretical models suggest that the superior femoral neck is a location of local weakness due to an age-related thinning of the cortex, and thus the site of hip fracture initiation. The purpose of this study was to experimentally determine the spatial and temporal characteristics of the macroscopic failure process during a simulated hip fracture that would occur as a result of a sideways fall. Twelve fresh frozen human cadaveric femora were used in this study. The femora were fractured in an apparatus designed to simulate a fall on the greater trochanter. Image sequences of the surface events related to the fractures were captured using two high-speed video cameras at 9111 Hz. The videos were analyzed with respect to time and load to determine the location and sequence of these events occurring in the proximal femur. The mean failure load was 4032 N (SD 370 N). The first surface events were identified in the superior femoral neck in eleven of the twelve specimens. Nine of these specimens fractured in a clear two-step process that initiated with a failure in the superior femoral neck, followed by a failure in the inferior femoral neck. This cadaveric model of hip fracture empirically confirms hypotheses that suggested that hip fractures initiate with a failure in the superior femoral neck where stresses are primarily compressive during a sideways fall impact, followed by a failure in the inferior neck where stresses are primarily tensile. Our results confirm the superolateral neck of the femur as an important region of interest for future hip fracture screening, prevention and treatment research.     

Finite element analysis of anterior tooth root stresses developed during endodontic treatment

Vertical tooth root fractures are diagnostically challenging, frustrating, and an increasingly common cause of failure of tooth restoration. These vertical root fractures have been associated with many causes, including the endodontic process itself. To investigate these endodontic causes, various phases of crown replacement for an anterior tooth were modeled using nonlinear, plane strain finite element (FE) analysis. Stresses developed during obturation, post positioning, crown placement, and masticatory and occlusal loading of the restored tooth were estimated using this analysis method. The minimum (compressive) principal stress was greatest during obturation of cones 1 and 2, and during mastication of the restored tooth. Although these stresses were significant (-150 to -280 MPa), they did not exceed the compressive strength of dentin. The maximum (tensile) principal stresses, 160 to 300 MPa, were also observed during obturation of cones 1 and 2. These peak stresses exceed the dentin tensile strength.     

Calibration of hyperelastic material properties of the human lumbar intervertebral disc under fast dynamic compressive loads

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.     

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

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.     

Critical evaluation of known bone material properties to realize anisotropic FE-simulation of the proximal femur

PURPOSE: In a meta-analysis of the literature we evaluated the present knowledge of the material properties of cortical and cancellous bone to answer the question whether the available data are sufficient to realize anisotropic finite element (FE)-models of the proximal femur. MATERIAL AND METHOD: All studies that met the following criteria were analyzed: Young's modulus, tensile, compressive and torsional strengths, Poisson's ratio, the shear modulus and the viscoelastic properties had to be determined experimentally. The experiments had to be carried out in a moist environment and at room temperature with freshly removed and untreated human cadaverous femurs. All material properties had to be determined in defined load directions (axial, transverse) and should have been correlated to apparent density (g/cm(3)), reflecting the individually variable and age-dependent changes of bone material properties. RESULTS: Differences in Young's modulus of cortical [cancellous] bone at a rate of between 33% (58%) (at low apparent density) and 62% (80%) (at high apparent density), are higher in the axial than in the transverse load direction. Similar results have been seen for the compressive strength of femoral bone. For the tensile and torsional strengths, Poisson's ratio and the shear modulus, only ultimate values have been found without a correlation to apparent density. For the viscoelastic behaviour of bone only data of cortical bone and in axial load direction have been described up to now. CONCLUSIONS: Anisotropic FE-models of the femur could be realized for most part with the summarized material properties of bone if characterized by apparent density and load directions. Because several mechanical properties have not been correlated to these main criteria, further experimental investigations will be necessary in future.     

Predicting the failure response of cement-bone constructs using a non-linear fracture mechanics approach

A non-linear fracture mechanics approach was used to predict the failure response of complex cement-bone constructs. A series of eight mechanical tests with a combination of tensile and shear loading along the cement-bone interface was performed. Each experiment was modeled using the finite element method with non-linear constitutive models at the cement-bone interface. Interface constitutive parameters were assigned based on the quantity of bone interdigitated with the cement. There was a strong correlation (r2 = 0.80) between experimentally measured and finite element predicted ultimate loads. The average error in predicted ultimate load was 23.9 percent. In comparison to the ultimate load predictions, correlations and errors for total energy to failure (r2 = 0.24, avg. error = 38.2 percent) and displacement at 50 percent of the ultimate load (r2 = 0.27, avg. error = 52.2 percent) were poor The results indicate that the non-linear constitutive laws could be useful in predicting the initiation and progression of interface failure of cemented bone-implant systems. However improvements in the estimation of post-yield interface properties from the quantity of bone interdigitated with cement are needed to enhance predictions of the overall failure response.     

Effect of displacement rate on the tensile mechanics of pediatric cervical functional spinal units

This study examined the effect of loading (displacement) rate on the tensile mechanics of cervical spine functional spinal units. A total of 40 isolated functional spinal units (two vertebrae and the adjoining soft tissues) from juvenile male baboons (10+/-0.6-human equivalent years old) were subjected to tensile loading spanning four orders of magnitude from 0.5 to 5000 mm/s. The stiffness, ultimate failure load, and corresponding displacement at failure were measured for each specimen and normalized by spinal geometry to examine the material properties as well as the structural properties. The tensile stiffness, failure load, normalized stiffness, and normalized failure load significantly increased (ANOVA, p<0.001) with increasing displacement rate. From the slowest to fastest loading rate, a two-fold increase in stiffness and four-fold increase in failure load were observed. The tensile failure strains (1.07+/-0.31 mm/mm strain) were not significantly correlated with loading rate (ANOVA, p=0.146). Both the functional (non-destructive stiffness and normalized stiffness) and failure mechanics of isolated functional spinal units exhibited a power-law relationship with displacement rate. Modeling efforts utilizing these rate-dependent characteristics will enhance our understanding of the tensile viscoelastic response of the spine and enable improved dynamic injury prevention schemes.     

Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur

Functional adaptation of the femur has been investigated in several studies by embedding bone remodelling algorithms in finite element (FE) models, with simplifications often made to the representation of bone’s material symmetry and mechanical environment. An orthotropic strain-driven adaptation algorithm is proposed in order to predict the femur’s volumetric material property distribution and directionality of its internal structures within a continuum. The algorithm was applied to a FE model of the femur, with muscles, ligaments and joints included explicitly. Multiple load cases representing distinct frames of two activities of daily living (walking and stair climbing) were considered. It is hypothesised that low shear moduli occur in areas of bone that are simply loaded and high shear moduli in areas subjected to complex loading conditions. In addition, it is investigated whether material properties of different femoral regions are stimulated by different activities. The loading and boundary conditions were considered to provide a physiological mechanical environment. The resulting volumetric material property distribution and directionalities agreed with ex vivo imaging data for the whole femur. Regions where non-orthogonal trabecular crossing has been documented coincided with higher values of predicted shear moduli. The topological influence of the different activities modelled was analysed. The influence of stair climbing on the properties of the femoral neck region is highlighted. It is recommended that multiple load cases should be considered when modelling bone adaptation. The orthotropic model of the complete femur is released with this study.     

Femoral strength is better predicted by finite element models than QCT and DXA.

Clinicians and patients would benefit if accurate methods of predicting and monitoring bone strength in-vivo were available. A group of 51 human femurs (age range 21-93; 23 females, 28 males) were evaluated for bone density and geometry using quantitative computed tomography (QCT) and dual energy X-ray absorptiometry (DXA). Regional bone density and dimensions obtained from QCT and DXA were used to develop statistical models to predict femoral strength ex vivo. The QCT data also formed the basis of a three-dimensional finite element (FE) models to predict structural stiffness. The femurs were separated into two groups; a model training set (n = 25) was used to develop statistical models to predict ultimate load, and a test set (n = 26) was used to validate these models. The main goal of this study was to test the ability of DXA, QCT and FE techniques to predict fracture load non-invasively, in a simple load configuration which produces predominantly femoral neck fractures. The load configuration simulated the single stance phase portion of normal gait; in 87% of the specimens, clinical appearing sub-capital fractures were produced. The training/test study design provided a tool to validate that the predictive models were reliable when used on specimens with "unknown" strength characteristics. The FE method explained at least 20% more of the variance in strength than the DXA models. Planned refinements of the FE technique are expected to further improve these results. Three-dimensional FE models are a promising method for predicting fracture load, and may be useful in monitoring strength changes in vivo.     

3D real-time micromechanical compressive behaviour of bone-cement interface: experimental and finite element studies

The integrity of bone-cement interface is essential for the long-term stability of cemented total joint arthroplasty. Although several studies have been carried out on bone-cement interface at continuum level, micromechanics of the interface has been studied only recently for tensile and shear loading cases. Fundamental studies of bone-cement interface at microstructural level are critical to the understanding of the failure processes of the interface, where multiple factors may contribute to failure. Here we present a micromechanical study of bone-cement interface under compression, which utilised in situ mechanical testing, time-lapsed microcomputed tomography (CT) and finite element (FE) modelling. Bovine trabecular bone was used to interdigitate with bone cement to obtain bone-cement interface samples, which were tested in step-wise compression using a custom-made loading stage within the μCT chamber. A finite element model was built from the CT images of one of the tested samples and loaded similarly as in the experiment. The simulated stress-displacement response fell within the range of the experimental responses, and the predicted local strain distribution correlated well with the failure pattern in the subject-specific experimental model. Damage evolution with load in the samples was monitored both experimentally and numerically. The results from the FE simulations further revealed the development of damage in the regions of interest during compression, which may be useful towards a micromechanics understanding of the failure processes at bone-cement interface.     

Use of a statistical model of the whole femur in a large scale,multi-model study of femoral neck fracture risk

Interpatient variability is often overlooked in orthopaedic computational studies due to the substantial challenges involved in sourcing and generating large numbers of bone models. A statistical model of the whole femur incorporating both geometric and material property variation was developed as a potential solution to this problem. The statistical model was constructed using principal component analysis, applied to 21 individual computer tomography scans. To test the ability of the statistical model to generate realistic, unique, finite element (FE) femur models it was used as a source of 1000 femurs to drive a study on femoral neck fracture risk. The study simulated the impact of an oblique fall to the side, a scenario known to account for a large proportion of hip fractures in the elderly and have a lower fracture load than alternative loading approaches. FE model generation, application of subject specific loading and boundary conditions, FE processing and post processing of the solutions were completed automatically. The generated models were within the bounds of the training data used to create the statistical model with a high mesh quality, able to be used directly by the FE solver without remeshing. The results indicated that 28 of the 1000 femurs were at highest risk of fracture. Closer analysis revealed the percentage of cortical bone in the proximal femur to be a crucial differentiator between the failed and non-failed groups. The likely fracture location was indicated to be intertrochantic. Comparison to previous computational, clinical and experimental work revealed support for these findings.