An anatomical subject-specific FE-model for hip fracture load prediction

In order to reduce the socio-economic burden induced by osteoporotic hip fractures, finite element models have been evaluated as an additional diagnostic tool for fracture prediction. For a future clinical application, the challenge is to reach the best compromise between model relevance and computing time. Based on this consideration, the current study focused on the development and validation of a subject-specific FE-model using an original parameterised generic model and a specific personalization method. A total of 39 human femurs were tested to failure under a quasi-static compression in stance configuration. The corresponding FE- models were generated and for each specimen the numerical fracture load (FFEM) was compared with the experimental value (FEXP), resulting in a significant correlation (FEXP = 1.006 FFEM with r2 = 0.87 and SEE = 1220 N, p < 0.05) obtained with a reasonable computing time (30 mn). Further in vivo study should confirm the ability of this FE-model to improve the fracture risk prediction.     

Specimen-specific modeling of hip fracture pattern and repair

Hip fracture remains a major health problem for the elderly. Clinical studies have assessed fracture risk based on bone quality in the aging population and cadaveric testing has quantified bone strength and fracture loads. Prior modeling has primarily focused on quantifying the strain distribution in bone as an indicator of fracture risk. Recent advances in the extended finite element method (XFEM) enable prediction of the initiation and propagation of cracks without requiring a priori knowledge of the crack path. Accordingly, the objectives of this study were to predict femoral fracture in specimen-specific models using the XFEM approach, to perform one-to-one comparisons of predicted and in vitro fracture patterns, and to develop a framework to assess the mechanics and load transfer in the fractured femur when it is repaired with an osteosynthesis implant. Five specimen-specific femur models were developed from in vitro experiments under a simulated stance loading condition. Predicted fracture patterns closely matched the in vitro patterns; however, predictions of fracture load differed by approximately 50% due to sensitivity to local material properties. Specimen-specific intertrochanteric fractures were induced by subjecting the femur models to a sideways fall and repaired with a contemporary implant. Under a post-surgical stance loading, model-predicted load sharing between the implant and bone across the fracture surface varied from 59%:41% to 89%:11%, underscoring the importance of considering anatomic and fracture variability in the evaluation of implants. XFEM modeling shows potential as a macro-level analysis enabling fracture investigations of clinical cohorts, including at-risk groups, and the design of robust implants.     

Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement

In the prediction of bone remodelling processes after total hip replacement (THR), modelling of the subject-specific geometry is now state-of-the-art. In this study, we demonstrate that inclusion of subject-specific loading conditions drastically influences the calculated stress distribution, and hence influences the correlation between calculated stress distributions and changes in bone mineral density (BMD) after THR.For two patients who received cementless THR, personalized finite element (FE) models of the proximal femur were generated representing the pre- and post-operative geometry. FE analyses were performed by imposing subject-specific three-dimensional hip joint contact forces as well as muscle forces calculated based on gait analysis data. Average values of the von Mises stress were calculated for relevant zones of the proximal femur. Subsequently, the load cases were interchanged and the effect on the stress distribution was evaluated. Finally, the subject-specific stress distribution was correlated to the changes in BMD at 3 and 6 months after THR.We found subject-specific differences in the stress distribution induced by specific loading conditions, as interchanging of the loading also interchanged the patterns of the stress distribution. The correlation between the calculated stress distribution and the changes in BMD were affected by the two-dimensional nature of the BMD measurement.Our results confirm the hypothesis that inclusion of subject-specific hip contact forces and muscle forces drastically influences the stress distribution in the proximal femur. In addition to patient-specific geometry, inclusion of patient-specific loading is, therefore, essential to obtain accurate input for the analysis of stress distribution after THR.     

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

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

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2–5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R² = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R² = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem.     

股骨骨折术后1 年随访骨愈合模型的有限元分析

目的：对股骨骨折髓内钉术后1 年骨愈合模型快速建模,通过有限元分析研究对比术前术后模型,通过术前判定内固定取出后骨折断端是否断裂。方法：运用Mimics、Geomagic Studio、Abaqus等软件采用快速个体化建模方法对股骨骨折髓内钉术后1年内固定取出术前后的多层螺旋CT 数据进行快速建立模型,术前模型模拟剥除钢板后进行有限元分析,施加人体单腿站立时的静力载荷和约束,并将分析结果与术后模型进行对比,观察米塞斯应力分布情况、最大值及其所处部位。结果：按照材料属性进行区别显示米赛斯应力的最大值及最小值,在不同应力载荷下,手术前后各类型材料的米赛斯应力最大值及最小值部位相同,各类型材料中,最大值均没有位于骨折断端,不同方法的最大应力值部位相近,均在股骨中远端1/4 交界处,手术前后应力分布基本相同。结论：采用个体化建模方法可以对骨折内固定取出前的骨愈合模型进行运算分析,快速预判术后是否导致骨折断端断裂。     

股骨骨折术后1年随访骨愈合模型的有限元分析

目的：对股骨骨折髓内钉术后1年骨愈合模型快速建模,通过有限元分析研究对比术前术后模型,通过术前判定内固定取出后骨折断端是否断裂。方法：运用Mimics、Geomagic Studio、Abaqus等软件采用快速个体化建模方法对股骨骨折髓内钉术后1年内固定取出术前后的多层螺旋CT数据进行快速建立模型,术前模型模拟剥除钢板后进行有限元分析,施加人体单腿站立时的静力载荷和约束,并将分析结果与术后模型进行对比,观察米塞斯应力分布情况、最大值及其所处部位。结果：按照材料属性进行区别显示米赛斯应力的最大值及最小值,在不同应力载荷下,手术前后各类型材料的米赛斯应力最大值及最小值部位相同,各类型材料中,最大值均没有位于骨折断端,不同方法的最大应力值部位相近,均在股骨中远端1/4交界处,手术前后应力分布基本相同。结论：采用个体化建模方法可以对骨折内固定取出前的骨愈合模型进行运算分析,快速预判术后是否导致骨折断端断裂。     

Computed tomography and automated image analysis of prehistoric femora

Non-invasive characterization of limb bone cross-sectional geometry would be useful for biomechanical analyses of skeletal collections. Computed tomography (CT) is potentially the method of choice. Additionally, CT images are suitable for automated analysis. CT is here shown to be both accurate and precise in the analysis of cross-sectional geometry of prehistoric femora. Beam hardening artifacts can be reduced by using a water bath. As the availability of CT for research increases, both bone density and geometry could be determined simultaneously with this method.     

Scanner influence on the mechanical response of QCT-based finite element analysis of long bones

Patient-specific QCT-based finite element (QCTFE) analyses enable highly accurate quantification of bone strength. We evaluated CT scanner influence on QCTFE models of long bones.A femur, humerus, and proximal femur without the head were scanned with K₂HPO₄ phantoms by seven CT scanners (four models) using typical clinical protocols. QCTFE models were constructed. The geometrical dimensions, as well as the QCT-values expressed in Hounsfield unit (HU) distribution was compared. Principal strains at representative regions of interest (ROIs), and maximum principal strains (associated with fracture risk) were compared. Intraclass correlation coefficients (ICCs) were calculated to evaluate strain prediction reliability for different scanners. Repeatability was examined by scanning the femur twice and comparing resulting QCTFE models.Maximum difference in geometry was 2.3%. HU histograms before phantom calibration showed wide variation between QCT scans; however, bone density histogram variability was reduced after calibration and algorithmic manipulation. Relative standard deviation (RSD) in principal strains at ROIs was <10.7%. ICC estimates between scanners were >0.9. Fracture-associated strain had 6.7%, 8.1%, and 13.3% maximum RSD for the femur, humerus, and proximal femur, respectively. The difference in maximum strain location was <2 mm. The average difference with repeat scans was 2.7%.Quantification of strain differences showed mean RSD bounded by ∼6% in ROIs. Fracture-associated strains in “regular” bones showed a mean RSD bounded by ∼8%. Strains were obtained within a ±10% difference relative to the mean; thus, in a longitudinal study only changes larger than 20% in the principal strains may be significant. ICCs indicated high reliability of QCTFE models derived from different scanners.     

Discrete element analysis is a valid method for computing joint contact stress in the hip before and after acetabular fracture

Evaluation of abnormalities in joint contact stress that develop after inaccurate reduction of an acetabular fracture may provide a potential means for predicting the risk of developing post-traumatic osteoarthritis. Discrete element analysis (DEA) is a computational technique for calculating intra-articular contact stress distributions in a fraction of the time required to obtain the same information using the more commonly employed finite element analysis technique. The goal of this work was to validate the accuracy of DEA-computed contact stress against physical measurements of contact stress made in cadaveric hips using Tekscan sensors. Four static loading tests in a variety of poses from heel-strike to toe-off were performed in two different cadaveric hip specimens with the acetabulum intact and again with an intentionally malreduced posterior wall acetabular fracture. DEA-computed contact stress was compared on a point-by-point basis to stress measured from the physical experiments. There was good agreement between computed and measured contact stress over the entire contact area (correlation coefficients ranged from 0.88 to 0.99). DEA-computed peak contact stress was within an average of 0.5 MPa (range 0.2–0.8 MPa) of the Tekscan peak stress for intact hips, and within an average of 0.6 MPa (range 0–1.6 MPa) for fractured cases. DEA-computed contact areas were within an average of 33% of the Tekscan-measured areas (range: 1.4–60%). These results indicate that the DEA methodology is a valid method for accurately estimating contact stress in both intact and fractured hips.     

Virtual structural analysis of tibial fracture healing from low-dose clinical CT scans

Quantitative assessment of bone fracture healing remains a significant challenge in orthopaedic trauma research. Accordingly, we developed a new technique for assessing bone healing using virtual mechano-structural analysis of computed tomography (CT) scans. CT scans from 19 fractured human tibiae at 12 weeks after surgery were segmented and prepared for finite element analysis (FEA). Boundary conditions were applied to the models to simulate a torsion test that is commonly used to access the structural integrity of long bones in animal models of fracture healing. The output of each model was the virtual torsional rigidity (VTR) of the healing zone, normalized to the torsional rigidity of each patient’s virtually reconstructed tibia. This provided a structural measure to track the percentage of healing each patient had undergone. Callus morphometric measurements were also collected from the CT scans. Results showed that at 12 weeks post-op, more than 75% of patients achieved a normalized VTR (torsional rigidity relative to uninjured bone) of 85% or above. The predicted intact torsional rigidities compared well with published cadaveric data. Across all patients, callus volume and density were weakly and non-significantly correlated with normalized VTR and time to clinical union. Conversely, normalized VTR was significantly correlated with time to union (R² = 0.383, p = 0.005). This suggests that fracture scoring methods based on the visual appearance of callus may not accurately predict mechanical integrity. The image-based structural analysis presented here may be a useful technique for assessment of bone healing in orthopaedic trauma research.     

A normative database of hip and knee joint biomechanics during dynamic tasks using anatomical regression prediction methods

Many methodologies exist to predict the hip joint center (HJC), of which regression based on anatomical landmarks appear most common. Despite the fact that predicted HJC locations vary depending upon chosen method, inter-study comparisons and inferences about populations are commonly made. The purpose of this study was to create a normative database of hip and knee biomechanics during walking, running, and single leg landings based on five commonly utilized HJC methods to serve as a reference for inter-study comparisons. Secondarily, we devised to provide comparisons of peak knee angles and hip angles, moments, and powers from the five HJC methods. Thirty healthy young adults performed walking, running, and single leg landing tasks at self-selected speeds (walking/running) and at 90% of their maximum jump height (landing). Three-dimensional motion capture and ground reaction forces were collected during all tasks. Five different HJC prediction methods: Bell, Davis, Hara, Harrington, and Greater Trochanter were implemented separately in a 6 degree of freedom model. Predicted HJC locations, direct kinematics, and inverse dynamics were computed for all tasks. Predicted HJC mediolateral, anteroposterior, and superior-inferior locations differed between methods by an average of 1.3, 2.9, and 1.4 cm, respectively. A database was created using the mean of all subjects for all five methods. In addition, one-way ANOVAs were used to compare triplanar peak angles, moments, and powers between the methods. The database of hip and knee biomechanics illustrates (1) variability between methods increases with more dynamic tasks (running/landing vs. walking) and (2) frontal and transverse plane hip and knee biomechanics are more variable between methods. Comparisons between methods found 38 and 16 main effect differences in hip and knee biomechanics, respectively. The Greater Trochanter method provided the most differences compared with other methods, while the Davis method provided the least differences. The database constructed provides an important reference for inter-study comparisons and details the impact of anatomical regression methods for predicting the HJC.     

Accounting for spatial variation of trabecular anisotropy with subject-specific finite element modeling moderately improves predictions of local subchondral bone stiffness at the proximal tibia

IntroductionPreviously, a finite element (FE) model of the proximal tibia was developed and validated against experimentally measured local subchondral stiffness. This model indicated modest predictions of stiffness (R² = 0.77, normalized root mean squared error (RMSE%) = 16.6%). Trabecular bone though was modeled with isotropic material properties despite its orthotropic anisotropy. The objective of this study was to identify the anisotropic FE modeling approach which best predicted (with largest explained variance and least amount of error) local subchondral bone stiffness at the proximal tibia.MethodsLocal stiffness was measured at the subchondral surface of 13 medial/lateral tibial compartments using in situ macro indentation testing. An FE model of each specimen was generated assuming uniform anisotropy with 14 different combinations of cortical- and tibial-specific density-modulus relationships taken from the literature. Two FE models of each specimen were also generated which accounted for the spatial variation of trabecular bone anisotropy directly from clinical CT images using grey-level structure tensor and Cowin’s fabric-elasticity equations. Stiffness was calculated using FE and compared to measured stiffness in terms of R² and RMSE%.ResultsThe uniform anisotropic FE model explained 53–74% of the measured stiffness variance, with RMSE% ranging from 12.4 to 245.3%. The models which accounted for spatial variation of trabecular bone anisotropy predicted 76–79% of the variance in stiffness with RMSE% being 11.2–11.5%.ConclusionsOf the 16 evaluated finite element models in this study, the combination of Synder and Schneider (for cortical bone) and Cowin’s fabric-elasticity equations (for trabecular bone) best predicted local subchondral bone stiffness.     

Soft tissue stiffness over the hip increases with age and its implication in hip fracture risk in older adults

Risk of hip fracture depends on the bone strength as well as the impact force delivered to the proximal femur during falls, and femoral soft tissue may help to reduce the hip fracture risk by attenuating the impact force. Femoral soft tissue stiffness was measured from a large sample, and compared how this was affected by age, gender and site.One hundred fifty healthy individuals (fifty-two young (aged between 19 and 29), forty-eight middle-aged (30–64), and fifty old (over 65)) participated. Each age group included an equal number of males and females. Using an automated hand-held indentation device, soft tissue stiffness was measured over twelve sites with respect to the greater trochanter (GT).For both left and right hips, the stiffness was associated with age (p < 0.0005), gender (p < 0.0005), and site (p < 0.0005). On average, the stiffness was 26% greater in older than young adults (321.5 versus 254.3 N/m). On average across twelve sites, the regression analysis indicated that the stiffness increases 1.33 N/m every year (“soft tissue stiffness over the hip = 1.33*age + 221.8”; R = 0.518, p < 0.0005). Furthermore, the stiffness was 18% greater in male than female (308.8 versus 262.6 N/m), and differed across twelve sites over the hip, being greatest (424.2 N/m) at the GT, and least (206.3 N/m) at the superior gluteal region.The results provide insights into the shock absorbing property of soft tissue over the hip, and inform the improvement of fall-related injury prevention interventions (i.e., hip protector, safe landing strategies) in older adults.     

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a ‘direct-mechanics’ method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.     

The mechanics of focal chondral defects in the hip

There is a mean incidence of osteoarthritis (OA) of the hip in 8% of the overall population. In the presence of focal chondral defects, defined as localized damage to the articular cartilage, there is an increased risk of symptomatic progression toward OA. This relationship between chondral defects and subsequent development of OA has led to substantial efforts to develop effective procedures for surgical cartilage repair. This study examined the effects of chondral defects and labral delamination on cartilage mechanics in the dysplastic hip during the gait cycle using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with simulated chondral defects at the chondrolabral junction on the posterior acetabulum during five distinct points in the gait cycle. Focal chondral defects increased maximum shear stress on the osteochondral surface of the acetabular cartilage, when compared to the intact case. This effect was amplified with labral delamination. Additionally, chondral defects increased the first principal Lagrange strain on the articular surface of the acetabular cartilage and labrum. Labral delamination relieved some of this tensile strain. As defect size was increased, contact stress increased in the medial zone of the acetabulum, while it decreased anteriorly. The results suggest that in the presence of chondral defects and labral delamination the cartilage experiences elevated tensile strains and shear and contact stress, which could lead to further damage of the cartilage, and subsequent arthritic progression. The framework presented here will serve as the procedure for future finite element studies on cartilage mechanics in hips with varying disease states with simulated chondral defects and labral tears.     

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.     

A biomechanical comparison of composite femurs and cadaver femurs used in experiments on operated hip fractures

Fourth generation composite femurs (4GCFs, models #3406 and #3403) simulate femurs of males <80 years with good bone quality. Since most hip fractures occur in old women with fragile bones, concern is raised regarding the use of standard 4GCFs in biomechanical experiments. In this study the stability of hip fracture fixations in 4GCFs was compared to human cadaver femurs (HCFs) selected to represent patients with hip fractures.Ten 4GCFs (Sawbones, Pacific Research Laboratories, Inc., Vashon, WA, USA) were compared to 24 HCFs from seven females and five males >60 years. Proximal femur anthropometric measurements were noted. Strain gauge rosettes were attached and femurs were mounted in a hip simulator applying a combined subject-specific axial load and torque. Baseline measurements of resistance to deformation were recorded. Standardized femoral neck fractures were surgically stabilized before the constructs were subjected to 20,000 load-cycles. An optical motion tracking system measured relative movements.Median (95% CI) head fragment migration was 0.8 mm (0.4 to 1.1) in the 4GCF group versus 2.2 mm (1.5 to 4.6) in the cadaver group (p=0.001). This difference in fracture stability could not be explained by observed differences in femoral anthropometry or potential overloading of 4GCFs. 4GCFs failed with fracture-patterns different from those observed in cadavers.To conclude, standard 4GCFs provide unrealistically stable bone-implant constructs and fail with fractures not observed in cadavers. Until a validated osteopenic or osteoporotic composite femur model is provided, standard 4GCFs should only be used when representing the biomechanical properties of young healthy femurs.