Development of a balanced experimental–computational approach to understanding the mechanics of proximal femur fractures

Method for computer tomography voxel-based finite element analysis and validation with digital image correlation system

Understanding the mechanical behavior of heterogeneous materials is becoming increasingly crucial across various fields, including aerospace engineering, composite materials development, geology, and biomechanics. While substantial literature exists on this topic, conventional methods often rely on commercial software packages. This study presents a framework for computed tomography (CT) scan-based finite element (FE) analysis of such materials using open-source software in most of the workflow. Our work focuses on three key aspects:1.Mesh generation that incorporates spatially varying mechanical properties and well-defined boundary conditions.2.Validation of the FE results through comparison with digital image correlation (DIC) system measurements.3.Open-source software utilization throughout the entire process, making it more accessible and cost-effective.This work aims to demonstrate the effectiveness of this framework for analyzing heterogeneous materials in various fields, offering a more accessible and affordable approach.

Prognostic Factors of Hip Fracture in Elderly: A Systematic Review

The hip fracture causes significant disabilities in many elderly people. Many studies around the world have identified various risk factors for the hip fracture. The aim of this study was to systematically investigate the risk factors of hip fractures. This study is a systematic review of risk factors for hip fractures. All published papers in English and Persian languages on patients in Iran and other countries between 2002 - 2022 were examined. The search strategy used keywords matching the mesh, including : predictors, hip fracture, and disability. Articles were selected from international databases (PubMed, Proquest ,Web of Sience, Scopus, Google scholar and Persian(Sid,Magiran), and the Newcastle Ottawa Scale was used to assess the risk of bias. The study has identified several factors that were significantly correlated with the risk of hip fracture, including age, cigarette and alcohol consumption, visual and hearing problems, low BMI levels, history of falling, weakness, and diseases such as stroke, cardiovascular disease, high blood pressure, arthritis, diabetes, dementia, Alzheimer's, Parkinson's, liver and kidney diseases, bone density, osteoporosis, vertebral fracture, and hyperthyroidism. However, the study did not find any significant correlations between the consumption of calcium and vitamin D, history of fractures, cognitive disorders, schizophrenia, and household income, and the risk of hip fracture. The results of this study reveal the determining role of some risk factors in hip fracture in older persons. Therefore, it is recommended that health policy makers provide the possibility of early intervention for some changeable factors.

Fracture in porous bone analysed with a numerical phase-field dynamical model

A dynamic phase-field fracture finite element model is applied to discretized high-resolution three-dimensional computed tomography images of human trabecular bone to analyse rapid bone fracture. The model is contrasted to quasi-static experimental results and a quasi-static phase-field finite element model. The experiment revealed complex stepwise crack evolution with multiple crack fronts, and crack arrests, as the global tensile displacement load was incrementally increased. The quasi-static phase-field fracture model captures the fractures in the experiment reasonably well, and the dynamic model converges towards the quasi-static model when mechanically loaded at low rates. At higher load rates, i.e., at larger impulses, inertia effects significantly contribute to an increased initial global stiffness, higher peak forces and a larger number of cracks spread over a larger volume. Since the fracture process clearly is different at large impulses compared to small impulses, it is concluded that dynamic fracture models are necessary when simulating rapid bone fracture.

Experimental validation of a voxel-based finite element model simulating femoroplasty of lytic lesions in the proximal femur

Femoroplasty is a procedure where bone cement is injected percutaneously into a weakened proximal femur. Uncertainty exists whether femoroplasty provides sufficient mechanical strengthening to prevent fractures in patients with femoral bone metastases. Finite element models are promising tools to evaluate the mechanical effectiveness of femoroplasty, but a thorough validation is required. This study validated a voxel-based finite element model against experimental data from eight pairs of human cadaver femurs with artificial metastatic lesions. One femur from each pair was left untreated, while the contralateral femur was augmented with bone cement. Finite element models accurately predicted the femoral strength in the defect (R² = 0.96) and augmented (R² = 0.93) femurs. The modelled surface strain distributions showed a good qualitative match with results from digital image correlation; yet, quantitatively, only moderate correlation coefficients were found for the defect (mean R² = 0.78) and augmented (mean R² = 0.76) femurs. This was attributed to the presence of vessel holes in the femurs and the jagged surface representation of our voxel-based models. Despite some inaccuracies in the surface measurements, the FE models accurately predicted the global bone strength and qualitative deformation behavior, both before and after femoroplasty. Hence, they can offer a useful biomechanical tool to assist clinicians in assessing the need for prophylactic augmentation in patients with metastatic bone disease, as well as in identifying suitable patients for femoroplasty.

Effect of fall direction on the lower hip fracture risk in athletes with different loading histories: A finite element modeling study in multiple sideways fall configurations

Physical loading makes bones stronger through structural adaptation. Finding effective modes of exercise to improve proximal femur strength has the potential to decrease hip fracture risk. Previous proximal femur finite element (FE) modeling studies have indicated that the loading history comprising impact exercises is associated with substantially higher fracture load. However, those results were limited only to one specified fall direction. It remains thus unclear whether exercise-induced higher fracture load depends on the fall direction. To address this, using magnetic resonance images of proximal femora from 91 female athletes (mean age 24.7 years with >8 years competitive career) and their 20 non-athletic but physically active controls (mean age 23.7 years), proximal femur FE models were created in 12 different sideways fall configurations. The athletes were divided into five groups by typical loading patterns of their sports: high-impact (H-I: 9 triple- and 10 high-jumpers), odd-impact (O-I: 9 soccer and 10 squash players), high-magnitude (H-M: 17 powerlifters), repetitive-impact (R-I: 18 endurance runners), and repetitive non-impact (R-NI: 18 swimmers). Compared to the controls, the FE models showed that the H-I and R-I groups had significantly (p < 0.05) higher fracture loads, 11-17% and 22-28% respectively, in all fall directions while the O-I group had significantly 10-11% higher fracture loads in four fall directions. The H-M and R-NI groups did not show significant benefit in any direction. Also, the analyses of the minimum fall strength (MFS) among these multiple fall configurations confirmed significantly 15%, 11%, and 14% higher MFSs in these impact groups, respectively, compared to the controls. These results suggest that the lower hip fracture risk indicated by higher fracture loads in athletes engaged in high impact or repetitive impact sports is independent of fall direction whereas the lower fracture risk attributed to odd-impact exercise is more modest and specific to the fall direction. Moreover, in concordance with the literature, the present study also confirmed that the fracture risk increases if the impact is imposed on the more posterolateral aspect of the hip. The present results highlight the importance of engaging in the impact exercises to prevent hip fractures and call for retrospective studies to investigate whether specific impact exercise history in adolescence and young adulthood is also associated with lower incidence of hip fractures in later life.

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.

Mechanism and microstructure based concept to predict skull fracture using a hybrid-experimental-modeling-computational approach

Cellular and tissue-scale indent/impact thresholds for different mechanisms of functional impairments to the brain would be the preferred method to predict head injuries, but a comprehensive understanding of the dominant possible injury mechanisms under multiaxial stress-states and rates is currently not available. Until then, skull fracture could serve as an indication of head injury. Therefore the ability to predict the initiation of skull fracture through finite element simulation can serve as an in silico tool for assessing the effectiveness of various head protection scenarios. For this objective, skull fracture initiation was represented with a microstructurally-inspired, mechanism-based (MIMB) failure surface assuming three different dominant mechanisms of skull failure: each element, with deformation and failure properties selected based on its microstructure, was allowed to fail either in tension, compression, or shear, corresponding to clinical linear, depressed or penetrating shear-plug failure (fracture), respectively. Microstructure-inspired a priori values for the initiation threshold of each mechanism, obtained previously from uniaxial and simple-shear experiments, were iterated and optimized for the predicted load-displacement to represent that of the corresponding indentation experiment. Element-level failure enabled the visualization of the evolution of fracture by different mechanisms. The final crack pattern at the time of macroscopic (clinically-identifiable) injury was compared between the simulation and experiment obtained through 3D tomography. Even though the timing was slightly different, the simulated prediction represented remarkably well the experimental crack pattern before the appearance of the catastrophic unstable fast crack in the experiment, thus validating the implemented hybrid-experimental-modeling-computational (HEMC) concept as a tool to predict skull fracture initiation.

Femoral strength and strains in sideways fall: Validation of finite element models against bilateral strain measurements

Low impact falls to the side are the main cause of hip fractures in elderly. Finite element (FE) models of the proximal femur may help in the assessment of patients at high risk for a hip fracture. However, extensive validation is essential before these models can be used in a clinical setting. This study aims to use strain measurements from bilateral digital image correlation to validate an FE model against ex vivo experimental data of proximal femora under a sideways fall loading condition. For twelve subjects, full-field strain measurements were available on the medial and lateral side of the femoral neck. In this study, subject-specific FE models were generated based on a consolidated procedure previously validated for stance loading. The material description included strain rate dependency and separate yield and fracture strain limits in tension and compression. FE predicted fracture force and experimentally measured peak forces showed a strong correlation (R² = 0.92). The FE simulations predicted the fracture initiation within 3 mm distance of the experimental fracture line for 8/12 subjects. The predicted and measured strains correlated well on both the medial side (R² = 0.87) and the lateral side (R² = 0.74). The lower correlation on the lateral side is attributed to the irregularity of the cortex and presence of vessel holes in this region. The combined validation against bilateral full-field strain measurements and peak forces has opened the door for a more elaborate qualitative and quantitative validation of FE models of femora under sideways fall loading.

Digital image correlation in dental materials and related research: A review

OBJECTIVE

Digital image correlation (DIC) is a non-contact image processing technique for full-field strain measurement. Although DIC has been widely used in engineering and biomechanical fields, it is in the spotlight only recently in dental materials. Therefore, the purpose of this review paper is introducing the working principle of the DIC technique with some modifications and providing further potential applications in various dental materials and related fields.

METHODS

The accuracy of the algorithm depending on the environmental characteristics of the DIC technique, as well as the advantages and disadvantages of strain measurement using optical measurements, have been elaborated in dental materials and related fields. Applications to those researches have been classified into the following categories: shrinkage behavior of light-cured resin composite, resin-tooth interface, mechanical properties of tooth structure, crack extension and elastic properties of dental materials, and deformation of dental restoration and prosthesis. This classification and discussion were performed using literature survey and review based on numerous papers in the international journals published over the past 20 years. The future directions for predicting the precise deformation of dental materials under various environments, as well as limitations of the DIC technique, was presented in this review.

RESULTS

The DIC technique was demonstrated as a more effective tool to measure full-field polymerization shrinkage of composite resin, even in a simulated clinical condition over the existing methods. Moreover, the DIC combined with other technologies can be useful to evaluate the mechanical behavior of material-tooth interface, dentine structure and restorative and prosthetic materials with high accuracy. Three-dimensional DIC using two cameras extended the measurement range in-plane to out-of-plane, enabling measure of the strain directly on the surface of dental restorations or prosthesis.

SIGNIFICANCE

DIC technique is a potential tool for measuring and predicting the full-field deformation/strain of dental materials and actual prostheses in diverse clinical conditions. The versatility of DIC can replace the existing complex sensor devices in those studies.

Damage tolerance and toughness of elderly human femora

Observations of elastic instability of trabecular bone cores supported the analysis of cortical thickness for predicting bone fragility of the hip in people over 60 years of age. Here, we falsified the hypothesis that elastic instability causes minimal energy fracture by analyzing, with a micrometric resolution, the deformation and fracture behavior of entire femora. Femur specimens were obtained from elderly women aged between 66 - 80 years. Microstructural images of the proximal femur were obtained under 3 - 5 progressively increased loading steps and after fracture. Bone displacements, strain, load bearing and energy absorption capacity were analyzed. Elastic instability of the cortex appeared at early loading stages in regions of peak compression. No elastic instability of trabecular bone was observed. The subchondral bone displayed local crushing in compression at early loading steps and progressed to 8 - 16% compression before fracture. The energy absorption capacity was proportional to the displacement. Stiffness decreased to near-zero values before fracture. Three-fourth of the fracture energy (10.2 - 20.2 J) was dissipated in the final 25% force increment. Fracture occurred in regions of peak tension and shear, adjacent to the location of peak compression, appearing immediately before fracture. Minimal permanent deformation was visible along the fracture surface. Elastic instability modulates the interaction between cortical and trabecular bone promoting an elastically stable fracture behavior of the femur organ, load bearing capacity, toughness, and damage tolerance. These findings will advance current methods for predicting hip fragility.

Biomechanical Comparison of a Notched Head Locking T-Plate and a Straight Locking Compression Plate in a Juxta-Articular Fracture Model

OBJECTIVE

 This investigation compared the biomechanical properties of a 2.0 mm locking compression notched head T-plate (NHTP) and 2.0 mm straight locking compression plate (LCP), in a simple transverse juxta-articular fracture model.

STUDY DESIGN

 Two different screw configurations were compared for the NHTP and LCP, modelling short (configuration 1) and long working length (configuration 2). Constructs were tested in compression, perpendicular and tension non-destructive four point bending and torsion. Plate surface strain was measured at 12 regions of interest (ROI) using three-dimensional digital image correlation. Stiffness and strain were compared between screw configurations within and between each plate.

RESULTS

 The LCP was stiffer than the NHTP in all three planes of bending and torsion (p < 0.05). The NHTP had greater strain than the LCP during compression bending and torsion at all ROI (p < 0.0005). The short working length was stiffer in all three planes of bending and in torsion (p < 0.05) than the longer working length for both plates. The long working length showed greater strain than the short working length at most ROI.

CONCLUSION

 In this experimental model, a 2.0 mm LCP with two screws in the short fragment was significantly stiffer and had lower plate strain than a 2.0 mm NHTP with three screws in the short fragment. Extending the working length significantly reduced construct stiffness and increased plate strain. These findings may guide construct selection.

Body Anthropometry and Bone Strength Conjointly Determine the Risk of Hip Fracture in a Sideways Fall

We hypothesize that variations of body anthropometry, conjointly with the bone strength, determine the risk of hip fracture. To test the hypothesis, we compared, in a simulated sideways fall, the hip impact energy to the energy needed to fracture the femur. Ten femurs from elderly donors were tested using a novel drop-tower protocol for replicating the hip fracture dynamics during a fall on the side. The impact energy was varied for each femur according to the donor’s body weight, height and soft-tissue thickness, by adjusting the drop height and mass. The fracture pattern, force, energy, strain in the superior femoral neck, bone morphology and microarchitecture were evaluated. Fracture patterns were consistent with clinically relevant hip fractures, and the superior neck strains and timings were comparable with the literature. The hip impact energy (11 – 95 J) and the fracture energy (11 – 39 J) ranges overlapped and showed comparable variance (CV = 69 and 61%, respectively). The aBMD-based definition of osteoporosis correctly classified 7 (70%) fracture/non-fracture cases. The incorrectly classified cases presented large impact energy variations, morphology variations and large subcortical voids as seen in microcomputed tomography. In conclusion, the risk of osteoporotic hip fracture in a sideways fall depends on both body anthropometry and bone strength.

Comparison of femur strain under different loading scenarios: Experimental testing

Bone fracture, formation and adaptation are related to mechanical strains in bone. Assessing bone stiffness and strain distribution under different loading conditions may help predict diseases and improve surgical results by determining the best conditions for long-term functioning of bone-implant systems. In this study, an experimentally wide range of loading conditions (56) was used to cover the directional range spanned by the hip joint force. Loads for different stance configurations were applied to composite femurs and assessed in a material testing machine. The experimental analysis provides a better understanding of the influence of the bone inclination angle in the frontal and sagittal planes on strain distribution and stiffness. The results show that the surface strain magnitude and stiffness vary significantly under different loading conditions. For the axial compression, maximal bending is observed at the mid-shaft, and bone stiffness is also maximal. The increased inclination leads to decreased stiffness and increased magnitude of maximum strain at the distal end of the femur. For comparative analysis of results, a three-dimensional, finite element model of the femur was used. To validate the finite element model, strain gauges and digital image correlation system were employed. During validation of the model, regression analysis indicated robust agreement between the measured and predicted strains, with high correlation coefficient and low root-mean-square error of the estimate. The results of stiffnesses obtained from multi-loading conditions experiments were qualitatively compared with results obtained from a finite element analysis of the validated model of femur with the same multi-loading conditions. When the obtained numerical results are qualitatively compared with experimental ones, similarities can be noted. The developed finite element model of femur may be used as a promising tool to estimate proximal femur strength and identify the best conditions for long-term functioning of the bone-implant system in future study.

Predicting population level hip fracture risk: a novel hierarchical model incorporating probabilistic approaches and factor of risk principles

Fall-related hip fractures are a major public health issue. While individual-level risk assessment tools exist, population-level predictive models could catalyze innovation in large-scale interventions. This study presents a hierarchical probabilistic model that predicts population-level hip fracture risk based on Factor of Risk (FOR) principles. Model validation demonstrated that FOR output aligned with a published dataset categorized by sex and hip fracture status. The model predicted normalized FOR for 100000 individuals simulating the Canadian older-adult population. Predicted hip fracture risk was higher for females (by an average of 38%), and increased with age (by15% per decade). Potential applications are discussed.

Individualized prediction of pedicle screw fixation strength with a finite element model

Pedicle screws are used for the treatment of a wide variety of spinal pathologies. A good screw holding power in bone is required for treatment success, but has so far not been predictable computationally. The goal of this study was to develop an automated tool able to predict patient-specific screw fixation strength through finite element simulation. We compared the simulation results with results from biomechanical pull-out tests performed on animal lumbar specimens. Experimental and simulation pull-out strengths were highly correlated [Formula: see text] and the mean error was 20.25%. The fixation strength was also associated to great extent with pull-out stiffness and strain energy, as well as the screw size and mean vertebral density.

EFFECT OF MATERIAL MODEL SELECTION ON LATERAL IMPACT SIMULATIONS OF PELVIC COMPLEX

Material mechanical behavior plays an important role in pelvic complex simulation under lateral impact. Aiming to investigate effects of material model selection on the responses of lateral impact simulations, a seating pelvic complex model was constructed. The model was subjected to a series of impacts at velocity of 3–10[Formula: see text]m/s, and two material models were, respectively, assigned to the pelvic bone to evaluate the accuracy of the simulation. The results showed that the pelvic response and fracture pattern with plastic–elastic material model agreed well with the literature, while linear elastic material model was dissatisfied factory, especially the pelvic response at low velocity deviated from most cadaveric test data. In addition, drastic change of arterial pressure was responsible for hemorrhages associated with pelvic fracture. Ligament loading sequence verified that the posterior pelvic ring bore the greatest amount of load during the impact. Based on the above findings, we concluded that a plastic–elastic with strain rate effect material model can improve the simulation accuracy of pelvic complex under lateral impact, and pelvic fracture pattern may help to estimate the parameters’ selection in impact simulation.

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.

Biofidelic finite element models for accurately classifying hip fracture in a retrospective clinical study of elderly women from the AGES Reykjavik cohort

Clinical retrospective studies have only reported limited improvements in hip fracture classification accuracy using finite element (FE) models compared to conventional areal bone mineral density (aBMD) measurements. A possible explanation is that state-of-the-art quasi-static models do not estimate patient-specific loads. A novel FE modeling technique was developed to improve the biofidelity of simulated impact loading from sideways falling. This included surrogate models of the pelvis, lower extremities, and soft tissue that were morphed based on subject anthropometrics. Hip fracture prediction models based on aBMD and FE measurements were compared in a retrospective study of 254 elderly female subjects from the AGES-Reykjavik study. Subject fragility ratio (FR) was defined as the ratio between the ultimate forces of paired biofidelic models, one with linear elastic and the other with non-linear stress-strain relationships in the proximal femur. The expected end-point value (EEV) was defined as the FR weighted by the probability of one sideways fall over five years, based on self-reported fall frequency at baseline. The change in maximum volumetric strain (ΔMVS) on the surface of the femoral neck was calculated between time of ultimate femur force and 90% post-ultimate force in order to assess the extent of tensile tissue damage present in non-linear models. After age-adjusted logistic regression, the area under the receiver-operator curve (AUC) was highest for ΔMVS (0.72), followed by FR (0.71), aBMD (0.70), and EEV (0.67), however the differences between FEA and aBMD based prediction models were not deemed statistically significant. When subjects with no history of falling were excluded from the analysis, thus artificially assuming that falls were known a priori with no uncertainty, a statistically significant difference in AUC was detected between ΔMVS (0.85), and aBMD (0.74). Multivariable linear regression suggested that the variance in maximum elastic femur force was best explained by femoral head radius, pelvis width, and soft tissue thickness (R2 = 0.79; RMSE = 0.46 kN; p < 0.005). Weighting the hip fracture prediction models based on self-reported fall frequency did not improve the models' sensitivity, however excluding non-fallers lead to significant differences between aBMD and FE based models. These findings suggest that an accurate assessment of fall probability is necessary for accurately identifying individuals predisposed to hip fracture.

Dynamic response and material sensitivity analysis of pelvic complex numerical model under side impact

BACKGROUND

The surrogate design and clinical diagnostic suggest that the pelvic dynamic response should be the basis of bone fracture mechanism study under side impact. Pelvic response indicators are the impact force, compression (C), viscous criterion (VC), bone stress, and bone strain. However, no evaluation of these indicators has been conducted.

OBJECTIVE

To evaluate pelvic response indicators under side impact.

METHODS

A sitting pelvic finite element (FE) complex model comprising bone, artery, ligaments, and soft tissue was constructed. The dynamic response of the model under side impact with initial velocity of 3 m/s was investigated and material sensitivity analysis was complemented by changing bone elastic modulus.

RESULTS

The pelvic FE model could predict response under side impact. Specifically, the indicators such as artery pressure and strain, together with the ligaments axial force and strain were provided. The sensitivity analysis showed the impact force, bone stress, and axial force were sensitive to the elastic modulus, whereas, C, VC, bone strain, and artery pressure were not.

CONCLUSIONS

The sitting FE model in this study can predict pelvic dynamic response, and C, VC, bone strain and artery pressure are proposed for pelvic tolerance instead of impact force under side impact.

Dynamic Simulation of Biomechanical Behaviour of the Pelvis in the Lateral Impact Loads

The objective of this study was to develop and validate a novel 3D dynamic model of a pelvic side-impactor system. The biomechanical responses of a pelvic flexible model (having .mnf file suffix) under the lateral impact load for predicting the bone fracture mechanism are investigated as well. The 3D solid model of the side-impactor system was imported into MSC/ADAMS software for analyzing the dynamic model, and the pelvic flexible model was extracted from the CT images of a Chinese female volunteer. The flexible model of the pelvis system was developed considering a wide range of mechanical properties in the bone complex and soft tissue to achieve a realistic biomechanical response during a lateral impact. Good agreements were achieved between the dynamic simulations and the experimental results of pelvic side impacts, in terms of the biomechanical criteria. The dynamic model of impactor system could be employed to investigate the hip protector effectiveness, improving the vehicle safety, and biomechanical response of the other human organs.

Impact loading history modulates hip fracture load and location: A finite element simulation study of the proximal femur in female athletes

Sideways falls impose high stress on the thin superolateral cortical bone of the femoral neck, the region regarded as a fracture-prone region of the hip. Exercise training is a natural mode of mechanical loading to make bone more robust. Exercise-induced adaptation of cortical bone along the femoral neck has been previously demonstrated. However, it is unknown whether this adaption modulates hip fracture behavior. The purpose of this study was to investigate the influence of specific exercise loading history on fall-induced hip fracture behavior by estimating fracture load and location with proximal femur finite element (FE) models created from magnetic resonance images (MRI) of 111 women with distinct exercise histories: 91 athletes (aged 24.7 ± 6.1 years, >8 years competitive career) and 20 women as controls (aged 23.7 ± 3.8 years). The athletes were divided into five groups based on typical loading patterns of their sports: high-impact (H-I: 9 triple-jumpers and 10 high jumpers), odd-impact (O-I: 9 soccer and 10 squash players), high-magnitude (H-M: 17 power-lifters), repetitive-impact (R-I: 18 endurance runners), and repetitive non-impact (R-NI: 18 swimmers). Compared to the controls, the H-I, O-I, and R-I groups had significantly higher (11-26%, p < 0.05) fracture loads. Also, the fracture location in the H-I and O-I groups was significantly more proximal (7-10%) compared to the controls. These results suggest that an exercise loading history of high impacts, impacts from unusual directions, or repetitive impacts increases the fracture load and may lower the risk of fall-induced hip fracture.

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.

DYNAMIC RESPONSE OF PELVIC COMPLEX FINITE ELEMENT STUDY AND VALIDATION UNDER SIDE IMPACT

Objective: To investigate and validate dynamic response of the pelvis, a finite element model of seated pelvic complex comprising of bone, ligaments, abdominal artery and soft tissue was developed and concurrently, a cadaver experiment was set up. Materials and Methods: Based on supine scanned CT images, we first developed an FE pelvic complex model and modified it to construct a seated pelvic model by anteriorly rotating the proximal femur to 90[Formula: see text]. For the cadaver experiment, a customized pelvic impact apparatus was designed and optical devices, strain gauges and pressure detectors were used to measure the pelvic response. Results: The results of the FE analysis and the cadaver tests were congruent in terms of impact force and fracture sites. Dynamic arterial response to the lateral impact showed hemodynamic instability that was displayed in pressure variation. The response of ligaments indicated that the posterior ligaments of pelvic ring experienced a larger amount of load. Conclusion: FE results provided the impact of ligaments and arteries besides impact force, compression (C) and viscous criterion (VC). Accordingly, the cadaver experiment measured arterial pressure, impact force, bone strain and compression. The compatibility between the FE and cadaver analyses demonstrates the high bio-fidelity of our pelvic complex model.

Spatial Differences in the Distribution of Bone Between Femoral Neck and Trochanteric Fractures

There is little knowledge about the spatial distribution differences in volumetric bone mineral density and cortical bone structure at the proximal femur between femoral neck fractures and trochanteric fractures. In this case-control study, a total of 93 women with fragility hip fractures, 72 with femoral neck fractures (mean ± SD age: 70.6 ± 12.7 years) and 21 with trochanteric fractures (75.6 ± 9.3 years), and 50 control subjects (63.7 ± 7.0 years) were included for the comparisons. Differences in the spatial distributions of volumetric bone mineral density, cortical bone thickness, cortical volumetric bone mineral density, and volumetric bone mineral density in a layer adjacent to the endosteal surface were investigated using voxel-based morphometry (VBM) and surface-based statistical parametric mapping (SPM). We compared these spatial distributions between controls and both types of fracture, and between the two types of fracture. Using VBM, we found spatially heterogeneous volumetric bone mineral density differences between control subjects and subjects with hip fracture that varied by fracture type. Interestingly, femoral neck fracture subjects, but not subjects with trochanteric fracture, showed significantly lower volumetric bone mineral density in the superior aspect of the femoral neck compared with controls. Using surface-based SPM, we found that compared with controls, both fracture types showed thinner cortices in regions in agreement with the type of fracture. Most outcomes of cortical and endocortical volumetric bone mineral density comparisons were consistent with VBM results. Our results suggest: 1) that the spatial distribution of trabecular volumetric bone mineral density might play a significant role in hip fracture; 2) that focal cortical bone thinning might be more relevant in femoral neck fractures; and 3) that areas of reduced cortical and endocortical volumetric bone mineral density might be more relevant for trochanteric fractures in Chinese women. © 2017 American Society for Bone and Mineral Research.

Exercise loading history and femoral neck strength in a sideways fall: A three-dimensional finite element modeling study

Over 90% of hip fractures are caused by falls. Due to a fall-induced impact on the greater trochanter, the posterior part of the thin superolateral cortex of the femoral neck is known to experience the highest stress, making it a fracture-prone region. Cortical geometry of the proximal femur, in turn, reflects a mechanically appropriate form with respect to habitual exercise loading. In this finite element (FE) modeling study, we investigated whether specific exercise loading history is associated with femoral neck structural strength and estimated fall-induced stresses along the femoral neck. One hundred and eleven three-dimensional (3D) proximal femur FE models for a sideways falling situation were constructed from magnetic resonance (MR) images of 91 female athletes (aged 24.7±6.1years, >8years competitive career) and 20 non-competitive habitually active women (aged 23.7±3.8years) that served as a control group. The athletes were divided into five distinct groups based on the typical loading pattern of their sports: high-impact (H-I: triple-jumpers and high-jumpers), odd-impact (O-I: soccer and squash players), high-magnitude (H-M: power-lifters), repetitive-impact (R-I: endurance runners), and repetitive non-impact (R-NI: swimmers). The von Mises stresses obtained from the FE models were used to estimate mean fall-induced stresses in eight anatomical octants of the cortical bone cross-sections at the proximal, middle, and distal sites along the femoral neck axis. Significantly (p<0.05) lower stresses compared to the control group were observed: the H-I group - in the superoposterior (10%) and posterior (19%) octants at the middle site, and in the superoposterior (13%) and posterior (22%) octants at the distal site; the O-I group - in the superior (16%), superoposterior (16%), and posterior (12%) octants at the middle site, and in the superoposterior (14%) octant at the distal site; the H-M group - in the superior (13%) and superoposterior (15%) octants at the middle site, and a trend (p=0.07, 9%) in the superoposterior octant at the distal site; the R-I group - in the superior (14%), superoposterior (23%) and posterior (22%) octants at the middle site, and in the superoposterior (19%) and posterior (20%) octants at the distal site. The R-NI group did not differ significantly from the control group. These results suggest that exercise loading history comprising various impacts in particular is associated with a stronger femoral neck in a falling situation and may have potential to reduce hip fragility.

In silico investigation of vertebroplasty as a stand-alone treatment for vertebral burst fractures

BACKGROUND

The use of percutaneous vertebroplasty as a stand-alone treatment for stable vertebral burst fractures has been investigated in vitro and in clinical studies. These studies present inconsistent results on the mechanical response of vertebroplasty-treated burst fractures. In addition, observations of the loss of sagittal alignment after vertebroplasty raise questions on the applicability of vertebroplasty for burst fractures. Therefore, the aim of this study was to investigate the mechanical stability of burst fractures after stand-alone treatment by vertebroplasty.

METHODS

Finite element simulations were performed with models generated from two laboratory-induced burst fractures in human thoracolumbar specimens. The burst fracture models were virtually injected with various cement volumes using a unipedicular or bipedicular approach. The models were subjected to four individual loads (compression, lateral bending, extension and torsion) and a multi-axial load case in the physiological range.

FINDINGS

All treated burst fractures showed improvements in stiffness and a reduction in inter-fragmentary displacements, thus potentially providing a suitable mechanical environment for fracture healing. However, large volumes of the trabecular bone (<43%), cement (<53%) and bone-cement composite (<58%) were predicted to experience strain levels exceeding the yield point. While damage was not specifically modeled, this implies a potential collapse of the treated vertebra due to local failure.

INTERPRETATION

To improve the primary stability and to prevent the collapse of treated burst fractures, the use of posterior instrumentation is suggested as an adjunct to vertebroplasty.

Full-field in vitro measurements and in silico predictions of strain shielding in the implanted femur after total hip arthroplasty

Alterations in bone strain as a result of implantation may contribute towards periprosthetic bone density changes after total hip arthroplasty. Computational models provide full-field strain predictions in implant–bone constructs; however, these predictions should be verified using experimental models wherever it is possible. In this work, finite element predictions of surface strains in intact and implanted composite femurs were verified using digital image correlation. Relationships were sought between post-implantation strain states across seven defined Gruen zones and clinically observed longer-term bone density changes. Computational predictions of strain distributions in intact and implanted femurs were compared to digital image correlation measurements in two regions of interest. Regression analyses indicated a strong linear correlation between measurements and predictions (R = 0.927 intact, 0.926 implanted) with low standard error (standard error = 38 µε intact, 26 µε implanted). Pre- to post-operative changes in measured and predicted surface strains were found to relate qualitatively to clinically observed volumetric bone density changes across seven Gruen zones: marked proximal bone density loss corresponded with a 50%−64% drop in surface strain, and slight distal density changes corresponded with 4%−14% strain increase. These results support the use of digital image correlation as a pre-clinical tool for predicting post-implantation strain shielding, indicative of long-term bone adaptations.

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.