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

In vitro mechanical analysis of X-shaped femoroplasty with polymethyl methacrylate boundary a fall on the greater trochanter. Injury 2023;54 Suppl 6:110747. [PMID: 38143120 DOI: 10.1016/j.injury.2023.04.034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

To evaluate with mechanical testing (MT) using synthetic femurs, an X-shaped femoroplasty technique with polymethyl methacrylate (PMMA), analyzing the results applied to the prophylaxis of proximal femur (PF) fractures caused by low-energy trauma. MT was performed simulating a fall on the greater trochanter, using fifteen Sawbones™ models. They were divided into three experimental groups (n = 5): control (DP) group, drilled without augmentation (DWA) group, and X-shaped augmentation (DX) group. Maximum load, stiffness, absorbed energy and displacement were analyzed primarily in all groups; and secondarily then, morphology and fracture type were verified in all groups while PMMA volume, temperature and time polymerization were analyzed only in the DX group. The MT results obtained for synthetic models respectively in the DP, DWA, and DX groups were: mean maximum load (5562.0 ± 464.8) N, (4798.0 ± 121.2) N, and (7132.0 ± 206.9) N; mean stiffness values (673 ± 64.34) N/mm, (636 ± 8.7) N/mm, and (738 ± 17.13) N/mm, and mean absorbed energy values (36,203 ± 3819) N.mm, (27,617 ± 3011) N.mm, (44,762 ± 3219) N.mm; mean displacement values (13.6 ± 1.45) N, (11.1 ± 0.5) N, and (13.2 ± 0.69) N. The mean volume, temperature reached during filling in the DX group were 9.8 mL, 42.54ºC with 1' 56" of polymerization. The fracture types were similar between the DP and DWA groups, affecting the trochanteric region, as distinctly to those in the DX group, which were restricted to the femoral neck. The values obtained in MT showed statistical significance when analyzed by one-way ANOVA (5%) for maximum load, stiffness, and absorbed energy between groups. In conclusion, X-shaped PMMA augmentation presents a protective biomechanical characteristic against PF fractures generated in synthetic models by boundary a fall on the greater trochanter.

The role of bone marrow on the mechanical properties of trabecular bone: a systematic review

Background

Accurate evaluation of the mechanical properties of trabecular bone is important, in which the internal bone marrow plays an important role. The aim of this systematic review is to investigate the roles of bone marrow on the mechanical properties of trabecular bone to better support clinical work and laboratory research.

Methods

A systematic review of the literature published up to June 2022 regarding the role of bone marrow on the mechanical properties of trabecular bone was performed, using PubMed and Web of Science databases. The journal language was limited to English. A total of 431 articles were selected from PubMed (n = 186), Web of Science (n = 244) databases, and other sources (n = 1).

Results

After checking, 38 articles were finally included in this study. Among them, 27 articles discussed the subject regarding the hydraulic stiffening of trabecular bone due to the presence of bone marrow. Nine of them investigated the effects of bone marrow on compression tests with different settings, i.e., in vitro experiments under unconfined and confined conditions, and computer model simulations. Relatively few controlled studies reported the influence of bone marrow on the shear properties of trabecular bone.

Conclusion

Bone marrow plays a non-neglectable role in the mechanical properties of trabecular bone, its contribution varies depending on the different loading types and test settings. To obtain the mechanical properties of trabecular bone comprehensively and accurately, the solid matrix (trabeculae) and fluid-like component (bone marrow) should be considered in parallel rather than tested separately.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-022-01051-1.

Which experimental procedures influence the apparent proximal femoral stiffness? A parametric study

Background

Experimental validation is the gold standard for the development of FE predictive models of bone. Employing multiple loading directions could improve this process. To capture the correct directional response of a sample, the effect of all influential parameters should be systematically considered. This study aims to determine the impact of common experimental parameters on the proximal femur’s apparent stiffness.

Methods

To that end, a parametric approach was taken to study the effects of: repetition, pre-loading, re-adjustment, re-fixation, storage, and μCT scanning as random sources of uncertainties, and loading direction as the controlled source of variation in both stand and side-fall configurations. Ten fresh-frozen proximal femoral specimens were prepared and tested with a novel setup in three consecutive sets of experiments. The neutral state and 15-degree abduction and adduction angles in both stance and fall configurations were tested for all samples and parameters. The apparent stiffness of the samples was measured using load-displacement data from the testing machine and validated against marker displacement data tracked by DIC cameras.

Results

Among the sources of uncertainties, only the storage cycle affected the proximal femoral apparent stiffness significantly. The random effects of setup manipulation and intermittent μCT scanning were negligible. The 15^∘ deviation in loading direction had a significant effect comparable in size to that of switching the loading configuration from neutral stance to neutral side-fall.

Conclusion

According to these results, comparisons between the stiffness of the samples under various loading scenarios can be made if there are no storage intervals between the different load cases on the same samples. These outcomes could be used as guidance in defining a highly repeatable and multi-directional experimental validation study protocol.

Fluid-structure interaction (FSI) modeling of bone marrow through trabecular bone structure under compression

The present study has sought to investigate the fluid characteristic and mechanical properties of trabecular bone using fluid-structure interaction (FSI) approach under different trabecular bone orientations. This method imposed on trabecular bone structure at both longitudinal and transverse orientations to identify effects on shear stress, permeability, stiffness and stress regarded to the trabeculae. Sixteen FSI models were performed on different range trabecular cubes of 27 mm³ with eight models developed for each longitudinal and transverse direction. Results show that there was a moderate correlation between permeability and porosity, and surface area in the longitudinal and transverse orientations. For the longitudinal orientation, the permeability values varied between 3.66 × 10^-8 and 1.9 × 10^-7 and the sheer stress values varied between 0.05 and 1.8 Pa, whilst for the transverse orientation, the permeability values varied between 5.95 × 10^-10 and 1.78 × 10^-8 and the shear stress values varied between 0.04 and 3.1 Pa. Here, transverse orientation limits the fluid flow from passing through the trabeculae due to high shear stress disturbance generated within the trabecular bone region. Compared to physiological loading direction (longitudinal orientation), permeability is higher within the range known to trigger a response in bone cells. Additionally, shear stresses also increase with bone surface area. This study suggests the shear stress within bone marrow in real trabecular architecture could provide the mechanical signal to marrow cells that leads to bone anabolism and can depend on trabecular orientation.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Biomechanical Computed Tomography analysis (BCT) for clinical assessment of osteoporosis

The surgeon general of the USA defines osteoporosis as "a skeletal disorder characterized by compromised bone strength, predisposing to an increased risk of fracture." Measuring bone strength, Biomechanical Computed Tomography analysis (BCT), namely, finite element analysis of a patient's clinical-resolution computed tomography (CT) scan, is now available in the USA as a Medicare screening benefit for osteoporosis diagnostic testing. Helping to address under-diagnosis of osteoporosis, BCT can be applied "opportunistically" to most existing CT scans that include the spine or hip regions and were previously obtained for an unrelated medical indication. For the BCT test, no modifications are required to standard clinical CT imaging protocols. The analysis provides measurements of bone strength as well as a dual-energy X-ray absorptiometry (DXA)-equivalent bone mineral density (BMD) T-score at the hip and a volumetric BMD of trabecular bone at the spine. Based on both the bone strength and BMD measurements, a physician can identify osteoporosis and assess fracture risk (high, increased, not increased), without needing confirmation by DXA. To help introduce BCT to clinicians and health care professionals, we describe in this review the currently available clinical implementation of the test (VirtuOst), its application for managing patients, and the underlying supporting evidence; we also discuss its main limitations and how its results can be interpreted clinically. Together, this body of evidence supports BCT as an accurate and convenient diagnostic test for osteoporosis in both sexes, particularly when used opportunistically for patients already with CT. Biomechanical Computed Tomography analysis (BCT) uses a patient's CT scan to measure both bone strength and bone mineral density at the hip or spine. Performing at least as well as DXA for both diagnosing osteoporosis and assessing fracture risk, BCT is particularly well-suited to "opportunistic" use for the patient without a recent DXA who is undergoing or has previously undergone CT testing (including hip or spine regions) for an unrelated medical condition.

Methods for Post Hoc Quantitative Computed Tomography Bone Density Calibration: Phantom-Only and Regression

Quantitative computed tomography (qCT) relies on calibrated bone mineral density data. If a calibration phantom is absent from the CT scan, post hoc calibration becomes necessary. Scanning a calibration phantom after-the-fact and applying that calibration to uncalibrated scans has been used previously. Alternatively, the estimated density is known to vary with CT settings, suggesting that it may be possible to predict the calibration terms using CT settings. This study compares a novel CT setting regression method for post hoc calibration to standard and post hoc phantom-only calibrations. Five cadaveric upper limbs were scanned at 11 combinations of peak tube voltage and current (80-140 kV and 100-300 mA) with two calibration phantoms. Density calibrations were performed for the cadaver scans, and scans of the phantoms alone. Stepwise linear regression determined if the calibration equation terms were predictable using peak tube voltage and current. Peak tube voltage, but not current, was significantly correlated with regression calibration terms. Calibration equation slope was significantly related to the type of phantom (p < 0.001), calibration method (p = 0.026), and peak tube voltage (p < 0.001), but not current (p = 1.000). The calibration equation vertical intercept was significantly related to the type of phantom (p < 0.001), and peak tube voltage (p = 0.006), but not calibration method (p = 0.682), or current (p = 0.822). Accordingly, regression can correlate peak tube voltage with density calibration terms. Suggesting that, while standard qCT calibration is preferable, regression calibration may be an acceptable post hoc method when necessary.

Previous Damage Accumulation Can Influence Femoral Fracture Strength: A Finite Element Study

To manage osteoporotic hip fracture risk, it is necessary to understand failure mechanisms of bone at both the material and organ level. The structural response of bone is dependent on load history. Repeated loading causes progressive microstructural cracking, resulting in reduced apparent-level stiffness and, if damage is significant, reductions to peak load bearing capability. However, the effect of previous damage accumulation has not been well explored at the organ level. It was hypothesized that femoral fracture load and fracture pattern may be sensitive to damage accumulation from previous loading events. Six cadaveric specimens were used to develop patient specific finite element (FE) models from quantitative tomographic (qCT) scans. Material properties were assigned from qCT intensity at each element location, and damage evolution was predicted using a previously validated quasi-brittle FE model. Three scenarios were investigated: stumble followed by another stumble (S-S), fall followed by another fall (F-F), and stumble followed by a fall (S-F). Fracture load and pattern were compared to FE predictions for a single stumble (S) or single fall (F) loading event. Most specimens were resilient to accumulated damage, showing little (<5%) change in fracture load from the multiple-load scenarios (S-S, F-F, and S-F) compared to an equivalent single load scenario (S or F). Only one specimen demonstrated moderate (5-15%) reductions in strength from all three multiple-load scenarios. However, two specimens experienced moderate (20-30%) increase in fracture load in some load cases. In these cases, initial damage caused the load to be more evenly distributed upon subsequent loading events. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:2197-2203, 2019.

The effect of boundary and loading conditions on patient classification using finite element predicted risk of fracture

BACKGROUND

Osteoporotic proximal femoral fractures associated to falls are a major health burden in the ageing society. Recently, bone strength estimated by finite element models emerged as a feasible alternative to areal bone mineral density as a predictor of fracture risk. However, previous studies showed that the accuracy of patients' classification under their risk of fracture using finite element strength when simulating posterolateral falls is only marginally better than that of areal bone mineral density. Patients tend to fall in various directions: since the predicted strength is sensitive to the fall direction, a prediction based on certain fall directions might not be fully representative of the physical event. Hence, side fall boundary conditions may not be completely representing the physical event.

METHODS

The effect of different side fall boundary and loading conditions on a retrospective cohort of 98 postmenopausal women was evaluated to test models' ability to discriminate fracture and control cases. Three different boundary conditions (Linear, Multi-point constraints and Contact model) were investigated under various anterolateral and posterolateral falls.

FINDINGS

The stratification power estimated by the area under the receiver operating characteristic curve was highest for Contact model (0.82), followed by Multi-point constraints and Linear models with 0.80. Both Contact and MPC models predicted high strains in various locations of the proximal femur including the greater trochanter, which has rarely reported previously.

INTERPRETATION

A full range of fall directions and less restrictive displacement constraints can improve the finite element strength ability to classify patients under their risk of fracture.

Influence of geometry on proximal femoral shaft strains: Implications for atypical femoral fracture

INTRODUCTION

Atypical femoral fractures (AFF) are characterized as low-energy fractures of the femoral shaft or subtrochanteric region. Femoral geometry is known to play a role in AFF risk; it is hypothesized that high-risk geometries are associated with elevated femoral shaft strain. However, it is not well known which geometric parameters have the greatest effect on strain, or whether interaction between parameters is significant. The purpose of this study was to thoroughly quantify the relationship between femoral geometry and diaphyseal strain, using patient specific finite element (FE) modelling in concert with parametric mesh morphing.

METHODS

Ten FE models were generated from computed tomography (CT) images of cadaveric femora. Heterogeneous material properties were assigned based on average CT intensities at element locations and models were subject to loads and boundary conditions representing the stance phase of gait. Mesh morphing was used to manipulate 8 geometric parameters: neck shaft angle (NSA), neck version angle (NV), neck length (NL), femoral length (FL), lateral bowing angle (L.Bow), anterior bowing angle (A.Bow), shaft diameter (S.Dia), and cortical bone thickness (C·Th). A 2-Level full factorial analysis was used to explore the effect of different combinations of physiologically realistic minimum and maximum values for each parameter. Statistical analysis (Generalized Estimating Equations) was used to assess main effects and first order interactions of each parameter.

RESULTS

Six independent parameters and seven interaction terms had statistically significant (p<0.05) effects on peak strain and strained volume. For both measures, the greatest changes were caused by S.Dia, L.Bow, and A.Bow, and/or first order interactions involving two of these variables.

CONCLUSIONS

As hypothesized, a large number of geometric measures (six) and first order interactions (seven) are associated with changes in femoral shaft strain. These measures can be evaluated radiographically, which may have important implications for future studies investigating AFF risk in clinical populations.

Validation of an alignment method using motion tracking system for in-vitro orientation of cadaveric hip joints with reduced set of anatomical landmarks

Accurate in-vitro orientation of cadaveric hip joints is challenging due to limited available anatomical landmarks. Published hip joint in-vitro investigations commonly lack details on methods used to achieve reported orientations and the accuracy with which the desired orientation has been achieved. The aim of this study was to develop an accurate method for orienting hip joints with limited anatomical landmarks for in-vitro investigations, and to compare this method against orientation using guiding axes and by visual approximation. The proposed orientation method resulted in orientation angles achieved to within one degree (SD ± 0.58°). For most specimens, orientation using physical tools resulted in errors of ±8° and ±12° in at least one of three orientation angles used to place the femur and pelvis in neutral orientation, respectively. Precision was also worse, with SDs ranging from ±1° to ±5° for orientation angles of femoral specimens and SDs ranging from ±1° to ±8° for pelvic specimens. The error in the orientation angles was worse for orientation by visual approximation and the range of SDs were greater for both the femur and pelvis. Finite element modeling was used to assess the effects of observed orientation errors, on prediction of fracture load. In most cases, the largest error in fracture load among all trials exceeded 30%, relative to a femur oriented without any error in the orientation angles.

Sideways fall-induced impact force and its effect on hip fracture risk: a review

Osteoporotic hip fracture, mostly induced in falls among the elderly, is a major health burden over the world. The impact force applied to the hip is an important factor in determining the risk of hip fracture. However, biomechanical researches have yielded conflicting conclusions about whether the fall-induced impact force can be accurately predicted by the available models. It also has been debated whether or not the effect of impact force has been considered appropriately in hip fracture risk assessment tools. This study aimed to provide a state-of-the-art review of the available methods for predicting the impact force, investigate their strengths/limitations, and suggest further improvements in modeling of human body falling.

METHODS

We divided the effective parameters on impact force to two categories: (1) the parameters that can be determined subject-specifically and (2) the parameters that may significantly vary from fall to fall for an individual and cannot be considered subject-specifically.

RESULTS

The parameters in the first category can be investigated in human body fall experiments. Video capture of real-life falls was reported as a valuable method to investigate the parameters in the second category that significantly affect the impact force and cannot be determined in human body fall experiments.

CONCLUSIONS

The analysis of the gathered data revealed that there is a need to develop modified biomechanical models for more accurate prediction of the impact force and appropriately adopt them in hip fracture risk assessment tools in order to achieve a better precision in identifying high-risk patients. Graphical abstract Impact force to the hip induced in sideways falls is affected by many parameters and may remarkably vary from subject to subject.

Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature

BACKGROUND

Finite element modeling of human bone provides a powerful tool to evaluate a wide variety of outcomes in a highly repeatable and parametric manner. These models are most often derived from computed tomography data, with mechanical properties related to bone mineral density (BMD) from the x-ray energy attenuation provided from this data. To increase accuracy, many researchers report the use of quantitative computed tomography (QCT), in which a calibration phantom is used during image acquisition to improve the estimation of BMD. Since model accuracy is dependent on the methods used in the calculation of BMD and density-mechanical property relationships, it is important to use relationships developed for the same anatomical location and using the same scanner settings, as these may impact model accuracy. The purpose of this literature review is to report the relationships used in the conversion of QCT equivalent density measures to ash, apparent, and/or tissue densities in recent finite element (FE) studies used in common density-modulus relationships. For studies reporting experimental validation, the validation metrics and results are presented.

RESULTS

Of the studies reviewed, 29% reported the use of a dipotassium phosphate (K2HPO4) phantom, 47% a hydroxyapatite (HA) phantom, 13% did not report phantom type, 7% reported use of both K2HPO4 and HA phantoms, and 4% alternate phantom types. Scanner type and/or settings were omitted or partially reported in 31% of studies. The majority of studies used densitometric and/or density-modulus relationships derived from different anatomical locations scanned in different scanners with different scanner settings. The methods used to derive various densitometric relationships are reported and recommendations are provided toward the standardization of reporting metrics.

CONCLUSIONS

This review assessed the current state of QCT-based FE modeling with use of clinical scanners. It was found that previously developed densitometric relationships vary by anatomical location, scanner type and settings. Reporting of all parameters used when referring to previously developed relationships, or in the development of new relationships, may increase the accuracy and repeatability of future FE models.

A two-level subject-specific biomechanical model for improving prediction of hip fracture risk

BACKGROUND

Sideways fall-induced hip fracture is a major worldwide health problem among the elderly population. However, all existing biomechanical models for predicting hip fracture mainly consider the femur related parameters. Their accuracy is limited as hip fracture is significantly affected by loading conditions as well. The objective of this study was to develop a biomechanical model for improving assessment of hip fracture risk by subject-specific prediction of fall-induced loading conditions.

METHOD

All information required to construct the models was extracted from the subject's whole-body and hip medical image in order to make the models subject-specific. Fall-induced hip fracture risk for eighty clinical cases was calculated under two sets of loading conditions: subject-specific determined by the proposed model, and non-subject-specific obtained from empirical functions. The predicted hip fracture risk indices were then compared with clinical observations.

FINDINGS

It was found that the subject-specific prediction of fall-induced loading conditions significantly improves the hip fracture risk assessment. Consistent to the clinical observations, the fracture risk predicted by the proposed model suggested that obesity is a protective factor for hip fracture and underweight subjects are more likely to experience a hip fracture.

INTERPRETATIONS

This study shows that hip fracture risk is affected by a number of factors, including body weight, body height, impact force, body mass index, hip soft tissue thickness, and bone quality. The proposed model provides a comprehensive, fast, accurate, and non-expensive method for prediction of hip fracture risk which should lead to more effective prevention of hip fractures.

Effects of hip abductor muscle forces and knee boundary conditions on femoral neck stresses during simulated falls

Through experiments that simulated sideways falls with a mechanical hip impact simulator, we demonstrated the protective effect of hip abductor muscle forces in reducing peak stresses at the femoral neck and the corresponding risk for hip fracture.

INTRODUCTION

Over 90% of hip fractures are due to falls, and an improved understanding the factors that separate injurious and non-injurious falls (via their influence on the peak stress generated at the femoral neck) may lead to improved risk assessment and prevention strategies. The purpose of this study was to measure the effect of muscle forces spanning the hip, and knee boundary conditions, on peak forces and estimated stresses at the femoral neck during simulated falls with a mechanical system.

METHODS

We simulated hip abductor muscle forces and knee boundary conditions with a mechanical hip impact simulator and measured forces and stresses at the femoral neck during sideways falls.

RESULTS

Peak compressive and tensile stresses, shear force, bending moment, and axial force are each associated with hip abductor muscle forces and knee boundary conditions (p < 0.0005). When muscle force increased from 400 to 1,200 N, peak compressive and tensile stresses decreased 24 and 56%, respectively. These effects were similar to the magnitude of decline in fracture strength associated with osteoporosis and arose from the tension-band effect of the muscle in reducing the bending moment by 37%. Furthermore, peak compressive and tensile stresses averaged 40 and 51% lower, respectively, in the free knee than fixed knee condition.

CONCLUSIONS

Contraction of the hip abductor muscles at the moment of impact during a fall, and landing with the knee free of constraints, substantially reduced peak compressive and tensile stresses at the femoral neck and risk for femoral fracture in a sideways fall.