Accurate in vitro identification of fracture onset in bones: failure mechanism of the proximal human femur

The Impact of Coronal Configuration of the Proximal Femur on its Mechanical Properties and the Validation of a New Theoretical Model: Finite Element Analysis and Biomechanical Examination

OBJECTIVE

This study aims to establish the coronal configuration of the proximal femur as an independent factor for its mechanical properties and provide validation for the theoretical model "fulcrum-balance-reconstruction."

METHODS

The digital 3D femur model constructed with the lower extremity high-resolution computed tomography of a senior subject was applied with the axial compression of 2100N under 5 different α angles of 10°, 5°, 0°, -5°, -10°. The equivalent stress distribution of the femoral geometric model under each angle were calculated. Under the same five α angles, fatigue test was performed on 15 composite artificial left femurs (three specimens in each angle group) to obtain the failure cycle and fracture site. The statistical analysis was accomplished using One-Way ANOVA.

RESULTS

The maximum stress of the entire femur in physiological angle (α = 10°) occurred below femoral neck with a value of 63.91 MPa. When the proximal femur is in extreme abducted angle (α = -10°), the maximum stress shift to the lower medial cortex of femoral shaft with a value of 105.2 MPa. As the α angle changed from 10° to -10°, the greater trochanteric region had the largest increment in maximum stress (2.78 times for cortex and 1.67 times for cancellous bone) locally at the proximal femur. The failure cycles of the artificial femurs with a variety of abduction angle were averagely 9126 ± 2453.87 (α = -10°), 58,112.33 ± 1293.84 (α = -5°), 92,879.67 ± 2398.54 (α = 0°), 172,045.3 ± 11011.11 (α = 5°), and 264,949.3 ± 35,067.26 (α = 10°), and the statistical analysis revealed that the α angle of the group of concern is proportional to the P value of the corresponding group compared to the 10° group(α = 5° & α = 10°, P = 0.01; α = 0 & α = 10°, P = 0.001; α = -5°, -10° & α = 10°, P < 0.001). In fatigue test, the fracture appeared on femoral neck for the α angles of 10° (three subcapital), 5° (two basal; one transcervical), and 0° (one transcervical). Fracture sites located at trochanteric region were observed with the more abducted angles including 0° (two subtrochanteric) and -5° (two intertrochanteric; one subtrochanteric). The fracture line was only found on femoral shaft in the -10° group.

CONCLUSION

With increasing hip abduction, the proximal femur shows declining mechanical properties, which suggests higher risk of hip fracture and increasement in the fraction of trochanteric fracture subtype. Furthermore, the hypothesis of "fulcrum-balance-reconstruction" was validated by our study to a certain extent.

Physical activity induced adaptation can increase proximal femur strength under loading from a fall onto the greater trochanter

Physical activity enhances proximal femur bone mass, but it remains unclear whether the benefits translate into an enhanced ability to resist fracture related loading. We recently used baseball pitchers as a within-subject controlled model to demonstrate physical activity induced proximal femur adaptation in regions associated with weight bearing and increased strength under single-leg stance loading. However, there was no measurable benefit to resisting common injurious loading (e.g. a fall onto the greater trochanter). A lack of power and a small physical activity effect size may have contributed to the latter null finding. Softball pitchers represent an alternative within-subject controlled model to explore adaptation of the proximal femur to physical activity, exhibiting greater dominant-to-nondominant (D-to-ND) leg differences than baseball pitchers. The current study used quantitative computed tomography, statistical parametric mapping, and subject-specific finite element (FE) modeling to explore adaptation of the proximal femur to physical activity in female softball pitchers (n = 25). Female cross-country runners (n = 15) were included as symmetrically loaded controls, showing very limited D-to-ND leg differences. Softball pitchers had D-to-ND leg differences in proximal femur, femoral neck, and trochanteric volumetric bone mineral density and content, and femoral neck volume. Voxel-based morphometry analyses and cortical bone mapping showed D-to-ND leg differences within a large region connecting the superior femoral head, inferior femoral neck and medial intertrochanteric region, and within the greater trochanter. FE modeling revealed pitchers had 19.4% (95%CI, 15.0 to 23.9%) and 4.9% (95%CI, 1.7 to 8.2%) D-to-ND leg differences in predicted ultimate strength under single-leg stance loading and a fall onto the greater trochanter, respectively. These data affirm the spatial and strength adaptation of the proximal femur to weight bearing directed loading and demonstrate that the changes can also have benefits, albeit smaller, on resisting loads associated with a sideways fall onto the greater trochanter.

Heterogeneous Spatial and Strength Adaptation of the Proximal Femur to Physical Activity: A Within-Subject Controlled Cross-Sectional Study

Physical activity (PA) enhances proximal femur bone mass, as assessed using projectional imaging techniques. However, these techniques average data over large volumes, obscuring spatially heterogeneous adaptations. The current study used quantitative computed tomography, statistical parameter mapping, and subject-specific finite element (FE) modeling to explore spatial adaptation of the proximal femur to PA. In particular, we were interested in adaptation occurring at the superior femoral neck and improving strength under loading from a fall onto the greater trochanter. High/long jump athletes (n = 16) and baseball pitchers (n = 16) were utilized as within-subject controlled models as they preferentially load their take-off leg and leg contralateral to their throwing arm, respectively. Controls (n = 15) were included but did not show any dominant-to-nondominant (D-to-ND) leg differences. Jumping athletes showed some D-to-ND leg differences but less than pitchers. Pitchers had 5.8% (95% confidence interval [CI] 3.9%-7.6%) D-to-ND leg differences in total hip volumetric bone mineral density (vBMD), with increased vBMD in the cortical compartment of the femoral neck and trochanteric cortical and trabecular compartments. Voxel-based morphometry analyses and cortical bone mapping showed pitchers had D-to-ND leg differences within the regions of the primary compressive trabeculae, inferior femoral neck, and greater trochanter but not the superior femoral neck. FE modeling revealed pitchers had 4.1% (95% CI 1.4%-6.7%) D-to-ND leg differences in ultimate strength under single-leg stance loading but no differences in ultimate strength to a fall onto the greater trochanter. These data indicate the asymmetrical loading associated with baseball pitching induces proximal femur adaptation in regions associated with weight bearing and muscle contractile forces and increases strength under single-leg stance loading. However, there were no benefits evident at the superior femoral neck and no measurable improvement in ultimate strength to common injurious loading during aging (ie, fall onto the greater trochanter), raising questions as to how to better target these variables with PA. © 2019 American Society for Bone and Mineral Research.

Numerical modelling of hip fracture patterns in human femur

BACKGROUND AND OBJECTIVE

Hip fracture morphology is an important factor determining the ulterior surgical repair and treatment, because of the dependence of the treatment on fracture morphology. Although numerical modelling can be a valuable tool for fracture prediction, the simulation of femur fracture is not simple due to the complexity of bone architecture and the numerical techniques required for simulation of crack propagation. Numerical models assuming homogeneous fracture mechanical properties commonly fail in the prediction of fracture patterns. This paper focuses on the prediction of femur fracture based on the development of a finite element model able to simulate the generation of long crack paths.

METHODS

The finite element model developed in this work demonstrates the capability of predicting fracture patterns under stance loading configuration, allowing the distinction between the main fracture paths: intracapsular and extracapsular fractures. It is worth noting the prediction of different fracture patterns for the same loading conditions, as observed during experimental tests.

RESULTS AND CONCLUSIONS

The internal distribution of bone mineral density and femur geometry strongly influences the femur fracture morphology and fracture load. Experimental fracture paths have been analysed by means of micro-computed tomography allowing the comparison of predicted and experimental crack surfaces, confirming the good accuracy of the numerical model.

Biomechanics of Femoral Neck Fractures and Implications for Fixation

Fractures of the femoral neck can occur in young healthy individuals due to high loads occurring during motor vehicle accidents, impacts, or falls. Failure forces are lower if impacts occur sideways onto the greater trochanter as compared with vertical loading of the hip. Bone density, bone geometry, and thickness of cortical bone at the femoral neck contribute to its mechanical strength. Femoral neck fractures in young adults require accurate reduction and stable internal fixation. The available techniques for fracture fixation at the femoral neck (cannulated screws, hip screw systems, proximal femur plates, and cephallomedullary nails) are reviewed with respect to their competence to provide biomechanical stability. Mechanically unstable fractures require a load-bearing implant, such as hip screws, with antirotational screws or intramedullary nails. Subcapital or transcervical fracture patterns and noncomminuted fractures enable load sharing and can be securely fixed with cannulated screws or solitary hip screw systems without compromising fixation stability.

Clinical hip fracture is accompanied by compression induced failure in the superior cortex of the femoral neck

Hip fractures pose a major health problem throughout the world due to their devastating impact. Current theories for why these injuries are so prevalent in the elderly point to an increased propensity to fall and decreases in bone mass with ageing. However, the fracture mechanisms, particularly the stress and strain conditions leading to bone failure at the hip remain unclear. Here, we directly examined the cortical bone from clinical intra-capsular hip fractures at a microscopic level, and found strong evidence of compression induced failure in the superior cortex. A total of 143 sections obtained from 24 femoral neck samples that were retrieved from 24 fracturing patients at surgery were examined using laser scanning confocal microscopy (LSCM) after fluorescein staining. The stained microcracks showed significantly higher density in the superior cortex than in the inferior cortex, indicating a greater magnitude of strain in the superior femoral neck during the failure-associated deformation and fracture process. The predominant stress state for each section was reconstructed based on the unique correlation between the microcrack pattern and the stress state. Specifically, we found clear evidence of longitudinal compression and buckling as the primary failure mechanisms in the superior cortex. These findings demonstrate the importance of microcrack analysis in studying clinical hip fractures, and point to the central role of the superior cortex failure as an important aspect of the failure initiation in clinical intra-capsular hip fractures.

Damage Identification on Vertebral Bodies During Compressive Loading Using Digital Image Correlation

MINI: Identifying fracture is important for understanding vertebral mechanics. Isolated cadaveric thoracolumbar vertebrae were compressed, and surface strains were measured using digital image correlation. Fracture locations from video analysis were qualitatively similar to the locations of high compressive strains and local damage occurred before the maximum force was reached.

STUDY DESIGN

Ex vivo compression experiments on isolated cadaveric vertebrae.

OBJECTIVE

To qualitatively compare the fracture locations identified in video analysis with the locations of high compressive strain measured with digital image correlation (DIC) on vertebral bodies and to evaluate the timing of local damage to the cortical shell relative to the global yield force.

SUMMARY OF BACKGROUND DATA

In previous ex vivo experiments, cortical bone fracture has been identified using various methods including acoustic emission sensors, strain gages, video analysis, or force signals. These methods are, however, limited in their ability to detect the location and timing of fracture. We propose use of DIC, a noncontact optical technique that measures surface displacement, to quantify variables related to damage.

METHODS

Isolated thoracolumbar human cadaveric vertebral bodies (n = 6) were tested in compression to failure at a quasi-static rate, and the force applied was measured using a load cell. The surface displacement and strain were measured using DIC. Video analysis was performed to identify fractures.

RESULTS

The location of fractures identified in the video corresponded well with the locations of high compressive strain on the bone. Before reaching the global yield force, more than 10% of the DIC measurements reached a minimum principal strain of 1.0%, a previously reported threshold for cortical bone damage.

CONCLUSION

DIC measurements provide an objective measure that can be used to identify the location and timing of fractures during ex vivo vertebral experiments. This is important for understanding fracture mechanics and for validating vertebral computational models that incorporate failure.

LEVEL OF EVIDENCE

N /A.

The inferomedial femoral neck is compromised by age but not disease: Fracture toughness and the multifactorial mechanisms comprising reference point microindentation

The influence of ageing on the fracture mechanics of cortical bone tissue is well documented, though little is known about if and how related material properties are further affected in two of the most prominent musculoskeletal diseases, osteoporosis and osteoarthritis (OA). The femoral neck, in close proximity to the most pertinent osteoporotic fracture site and near the hip joint affected by osteoarthritis, is a site of particular interest for investigation. We have recently shown that Reference Point micro-Indentation (RPI) detects differences between cortical bone from the femoral neck of healthy, osteoporotic fractured and osteoarthritic hip replacement patients. RPI is a new technique with potential for in vivo bone quality assessment. However, interpretation of RPI results is limited because the specific changes in bone properties with pathology are not well understood and, further, because it is not conclusive what properties are being assessed by RPI. Here, we investigate whether the differences previously detected between healthy and diseased cortical bone from the femoral neck might reflect changes in fracture toughness. Together with this, we investigate which additional properties are reflected in RPI measures. RPI (using the Biodent device) and fracture toughness tests were conducted on samples from the inferomedial neck of bone resected from donors with: OA (41 samples from 15 donors), osteoporosis (48 samples from 14 donors) and non age-matched cadaveric controls (37 samples from 10 donoros) with no history of bone disease. Further, a subset of indented samples were imaged using micro-computed tomography (3 osteoporotic and 4 control samples each from different donors) as well as fluorescence microscopy in combination with serial sectioning after basic fuchsin staining (7 osteoporotic and 5 control samples from 5 osteoporotic and 5 control donors). In this study, the bulk indentation and fracture resistance properties of the inferomedial femoral neck in osteoporotic fracture, severe OA and control bone were comparable (p > 0.05 for fracture properties and <10% difference for indentation) but fracture toughness reduced with advancing age (7.0% per decade, r = -0.36, p = 0.029). Further, RPI properties (in particular, the indentation distance increase, IDI) showed partial correlation with fracture toughness (r = -0.40, p = 0.023) or derived elastic modulus (r = -0.40, p = 0.023). Multimodal indent imaging revealed evidence of toughening mechanisms (i.e. crack deflection, bridging and microcracking), elastoplastic response (in terms of the non-conical imprint shape and presence of pile-up) and correlation of RPI with damage extent (up to r = 0.79, p = 0.034) and indent size (up to r = 0.82, p < 0.001). Therefore, crack resistance, deformation resistance and, additionally, micro-structure (porosity: r = 0.93, p = 0.002 as well as pore proximity: r = -0.55, p = 0.027 for correlation with IDI) are all contributory to RPI. Consequently, it becomes clear that RPI measures represent a multitude of properties, various aspects of bone quality, but are not necessarily strongly correlated to a single mechanical property. In addition, osteoporosis or osteoarthritis do not seem to further influence fracture toughness of the inferomedial femoral neck beyond natural ageing. Since bone is highly heterogeneous, whether this finding can be extended to the whole femoral neck or whether it also holds true for other femoral neck quadrants or other material properties remains to be shown.

Physical Activity for Strengthening Fracture Prone Regions of the Proximal Femur

PURPOSE OF REVIEW

Physical activity improves proximal femoral bone health; however, it remains unclear whether changes translate into a reduction in fracture risk. To enhance any fracture-protective effects of physical activity, fracture prone regions within the proximal femur need to be targeted.

RECENT FINDINGS

The proximal femur is designed to withstand forces in the weight-bearing direction, but less so forces associated with falls in a sideways direction. Sideways falls heighten femoral neck fracture risk by loading the relatively weak superolateral region of femoral neck. Recent studies exploring regional adaptation of the femoral neck to physical activity have identified heterogeneous adaptation, with adaptation principally occurring within inferomedial weight-bearing regions and little to no adaptation occurring in the superolateral femoral neck. There is a need to develop novel physical activities that better target and strengthen the superolateral femoral neck within the proximal femur. Design of these activities may be guided by subject-specific musculoskeletal modeling and finite-element modeling approaches.

Focal osteoporosis defects play a key role in hip fracture

BACKGROUND

Hip fractures are mainly caused by accidental falls and trips, which magnify forces in well-defined areas of the proximal femur. Unfortunately, the same areas are at risk of rapid bone loss with ageing, since they are relatively stress-shielded during walking and sitting. Focal osteoporosis in those areas may contribute to fracture, and targeted 3D measurements might enhance hip fracture prediction. In the FEMCO case-control clinical study, Cortical Bone Mapping (CBM) was applied to clinical computed tomography (CT) scans to define 3D cortical and trabecular bone defects in patients with acute hip fracture compared to controls. Direct measurements of trabecular bone volume were then made in biopsies of target regions removed at operation.

METHODS

The sample consisted of CT scans from 313 female and 40 male volunteers (158 with proximal femoral fracture, 145 age-matched controls and 50 fallers without hip fracture). Detailed Cortical Bone Maps (c.5580 measurement points on the unfractured hip) were created before registering each hip to an average femur shape to facilitate statistical parametric mapping (SPM). Areas where cortical and trabecular bone differed from controls were visualised in 3D for location, magnitude and statistical significance. Measures from the novel regions created by the SPM process were then tested for their ability to classify fracture versus control by comparison with traditional CT measures of areal Bone Mineral Density (aBMD). In women we used the surgical classification of fracture location ('femoral neck' or 'trochanteric') to discover whether focal osteoporosis was specific to fracture type. To explore whether the focal areas were osteoporotic by histological criteria, we used micro CT to measure trabecular bone parameters in targeted biopsies taken from the femoral heads of 14 cases.

RESULTS

Hip fracture patients had distinct patterns of focal osteoporosis that determined fracture type, and CBM measures classified fracture type better than aBMD parameters. CBM measures however improved only minimally on aBMD for predicting any hip fracture and depended on the inclusion of trabecular bone measures alongside cortical regions. Focal osteoporosis was confirmed on biopsy as reduced sub-cortical trabecular bone volume.

CONCLUSION

Using 3D imaging methods and targeted bone biopsy, we discovered focal osteoporosis affecting trabecular and cortical bone of the proximal femur, among men and women with hip fracture.

Can CT image deblurring improve finite element predictions at the proximal femur?

The aim of this study was to determine if a CT image deblurring algorithm can improve CT-based FE modelling accuracy at the proximal femur. Experimental data (CT scans of fourteen proximal fresh-frozen cadaveric femurs, non-destructive surface strain measurements in stance and sideways fall loading configurations on all femurs, and failure loads obtained in stance for seven specimens, in sideways fall for the other seven) were taken from a recent study (Schileo et al., 2014). An estimate of the 3D Point Spread Function for each CT scan was used within a deconvolution solver to perform the deblurring. The most proximal regions of three specimens were scanned using an HRpQCT scanner and compared to the original and deblurred CT images to quantify errors in bone contour estimates and determine correlation of intensity values within the bone contours. Subject-specific FE models of the proximal femur were generated. The accuracy of deblurred FE predictions against experimental measurements was compared to the published (non-deblurred) FE results. When compared to HRpQCT, CT deblurring led to lower mean surface distances (0.31 vs. 0.49mm) and higher CT intensity correlations with respect to the original CT. All indicators of strain prediction accuracy were significantly improved in deblurred FE models, more markedly at the femoral neck (peak error reduced by 38%). Failure load prediction, based on a simple elastic limit model, was also improved in deblurred FE models, although differently for stance and sideways fall loading conditions. In stance, correlation was unchanged, but specimen-wise errors were reduced (mean error 10% vs. 15%). In sideways fall, correlation notably increased (R(2)=0.95 vs. 0.81), despite a general overestimation of failure load. In summary, the proposed CT deblurring technique yielded moderate but significant improvements in FE predictions, and may thus be considered a first step toward the improvement of CT-based FE models of the human femur.

Strain distribution in the proximal Human femur during in vitro simulated sideways fall

This study assessed: (i) how the magnitude and direction of principal strains vary for different sideways fall loading directions; (ii) how the principal strains for a sideways fall differ from physiological loading directions; (iii) the fracture mechanism during a sideways fall. Eleven human femurs were instrumented with 16 triaxial strain gauges each. The femurs were non-destructively subjected to: (a) six loading configurations covering the range of physiological loading directions; (b) 12 configurations simulating sideways falls. The femurs were eventually fractured in a sideways fall configuration while high-speed cameras recorded the event. When the same force magnitude was applied, strains were significantly larger in a sideways fall than for physiological loading directions (principal compressive strain was 70% larger in a sideways fall). Also the compressive-to-tensile strain ratio was different: for physiological loading the largest compressive strain was only 30% larger than the largest tensile strain; but for the sideways fall, compressive strains were twice as large as the tensile strains. Principal strains during a sideways fall were nearly perpendicular to the direction of principal strains for physiological loading. In the most critical regions (medial part of the head-neck) the direction of principal strain varied by less than 9° between the different physiological loading conditions, whereas it varied by up to 17° between the sideways fall loading conditions. This was associated with a specific fracture mechanism during sideways fall, where failure initiated on the superior-lateral side (compression) followed by later failure of the medially (tension), often exhibiting a two-peak force-displacement curve.

To what extent can linear finite element models of human femora predict failure under stance and fall loading configurations?

Proximal femur strength estimates from computed tomography (CT)-based finite element (FE) models are finding clinical application. Published models reached a high in-vitro accuracy, yet many of them rely on nonlinear methodologies or internal best-fitting of parameters. The aim of the present study is to verify to what extent a linear FE modelling procedure, fully based on independently determined parameters, can predict the failure characteristics of the proximal femur in stance and sideways fall loading configurations. Fourteen fresh-frozen cadaver femora were CT-scanned. Seven femora were tested to failure in stance loading conditions, and seven in fall. Fracture was monitored with high-speed videos. Linear FE models were built from CT images according to a procedure already validated in the prediction of strains. An asymmetric maximum principal strain criterion (0.73% tensile, 1.04% compressive limit) was used to define a node-based risk factor (RF). FE-predicted failure load, mode (tensile/compressive) and location were determined from the first node reaching RF=1. FE-predicted and measured failure loads were highly correlated (R(2)=0.89, SEE=814N). In all specimens, FE models correctly identified the failure mode (tensile in stance, compressive in fall) and the femoral region where fracture started (supero-lateral neck aspect). The location of failure onset was accurately predicted in eight specimens. In summary, a simple FE model, adaptable in the future to multiple loads (e.g. including muscles), was highly correlated with experimental failure in two loading conditions on specimens ranging from normal to osteoporotic. Thus, it can be suitable for use in clinical studies.

DXA and pQCT predict pertrochanteric and not femoral neck fracture load in a human side-impact fracture model

The validity of dual energy X-ray absorptiometry (DXA) and peripheral quantitative computed tomography (pQCT) measurements as predictors of pertrochanteric and femoral neck fracture loads was compared in an experimental simulation of a fall on the greater trochanter. 65 proximal femora were harvested from patients at autopsy. All specimens were scanned with use of DXA for areal bone mineral density and pQCT for volumetric densities at selected sites of the proximal femur. A three-point bending test simulating a side-impact was performed to determine fracture load and resulted in 16 femoral neck and 49 pertrochanteric fractures. Regression analysis revealed that DXA BMD trochanter was the best variable at predicting fracture load of pertrochanteric fractures with an adjusted R(2) of 0.824 (p < 0.0001). There was no correlation between densitometric parameters and the fracture load of femoral neck fractures. A significant correlation further was found between body weight, height, femoral head diameter, and neck length on the one side and fracture load on the other side, irrespective of the fracture type. Clinically, the DXA BMD trochanter should be favored and integrated routinely as well as biometric and geometric parameters, particularly in elderly people with known osteoporosis at risk for falls.