The incidence of hip fracture in western society is becoming an increasingly serious concern as the population ages. While the economic and social costs are on the rise, the number of patients receiving treatment for bone fragility is in decline. Addressing this treatment gap should begin by investigating the relatively weak diagnostic technology currently being utilized in standard clinical practice. With problems including low sensitivity and a high rate of over diagnosis, it is clear that better tools are required for identifying individuals predisposed to hip fracture. Finite element (FE) models of the proximal femur derived from computed tomography have been proposed as a potential alternative. These simulations have outperformed areal bone mineral density (aBMD) measurements in predicting femoral strength within cadaver experiments, however FE models have not significantly outperformed conventional density measurements when predicting hip fracture from retrospective clinical data.

The objective of this thesis was to address the limitations of state-of-the-art FE simulations of hip fracture in order to develop an improved hip fracture prediction model. A hip fracture is the result of a fall onto the hip, leading to an impact force greater than the load bearing capacity of the proximal femur. The three components of hip fracture, i.e. fall risk, fall severity, and bone strength, have typically been investigated individually, and no study thus far has attempted to fully characterize risk of fracture using a holistic understanding of all three aspects. This thesis attempts to address this knowledge gap, and was divided into three aims, each designed to enhance the biofidelity of FE simulations:​

In the first study, the anisotropy of trabecular bone was measured within numerous sub-volumes extracted from high-resolution scans using the conventional mean intercept length measurement. These anisotropic properties were subsequently mapped to FE models of the proximal femur, loaded in a sideways fall configuration. For the first time, isotropic and anisotropic models were compared using different material mapping strategies representing the extremes of previously published modeling techniques. The resulting whole bone stiffness and surface strains were then compared to cadaver experiments tested with similar boundary conditions. The adv dition of anisotropy had very little effect on whole bone stiffness and only a small effect on the resulting surface strain. Differences in principal compressive strain were identified in the femoral head, neck, and greater trochanter. The study concluded that anisotropic material properties, mapped using morphological measurements of high-resolution CT scans, has little impact on macroscopic, organ-level properties, but anisotropy could still affect localized internal strains, and could therefore be relevant depending on the modeling objectives.

In the second study, a novel material mapping strategy was presented for explicit FE models of the proximal femur, validated with drop tower experiments. These non-linear material properties were designed to account for large deformations that could occur during an impact load on the hip. These properties included tensile damage, compressive densification, as well as tension/compression asymmetry and strain rate dependency. For the first time, the ultimate force predicted by the dynamic FE models was correlated with the ultimate force measured in dynamic experiments. Additionally, the simulated impulse response was strongly correlated with the force-time response measured in the drop tower experiments. Thus, these results represent the current benchmark in dynamic FE modeling of the proximal femur. Compressive strain rates over 100/s were observed in elements located in the femoral head, neck, and greater trochanter, suggesting that a better understanding of the strain rate dependency of bone tissue is still required for these loading rates. The dynamic models were able to simulate fracture patterns initiating in the subcapital region, but were not able to predict inter-trochanteric fractures, suggesting that the material mapping strategy could still be improved in this anatomical region.

The third study presented a biofidelic FE modeling technique that included models of the pelvis, soft tissue, and lower extremities that were morphed based on subject-specific biometrics. This approach was tested in a large, retrospective, clinical study to determine whether hip fracture classification models based on FE simulations provided a more accurate prediction of hip fracture incidence compared to clinical standard aBMD measurements. Logistic regression models based on simulated ultimate femur force, and surface strain in the femoral neck, marginally outperformed aBMD in terms of sensitivity and specificity, however the differences were not significant until subjects that did not report falling at baseline were excluded from the analysis. This resulted in a statistically significant difference between the simulated surface strain at the femoral neck and total femur aBMD, with AUC values of 0.85 and 0.74, respectively. These results indicate that subject-specific fall risk must be accurately estimated before mechanical assessments of fracture risk can be compared. This also suggested that fall risk could be a useful parameter for pre-screening hip fracture risk in the population. Additionally, the large number of biofidelic models tested in this study provided a novel measurement of subject-specific impact force transmitted to the proximal femur. A linear regression model found that the variance in ultimate femur force could be explained as a function of soft tissue thickness, pelvis width, and femoral head radius (R² = 0.79;​ RMSE = 0.46 kN).

This thesis has made several important contributions to the field of hip fracture biomechanics. It has demonstrated that simulating anisotropic tissue-level properties had minimal effect on organ-level mechanics, relative to changes in the material properties. A novel material mapping strategy for explicit FE models of the proximal femur was designed in the absence of model tuning, and subsequently validated against drop tower experiments. The importance of subject-specific soft tissue thickness and fall risk was demonstrated using biofidelic FE models tested in a large retrospective study, suggesting that femoral strength alone is insufficient for accurately characterizing hip fracture risk. Finally, this thesis has presented the strongest evidence, to date, that FE simulations can predict hip fracture more accurately than conventional aBMD measurements. However, this observation was statistically significant only when subjects with high risk of falling were tested. Looking forward, the retrospective study presented in this thesis could serve as the baseline for future in silico trials simulating prophylactic interventions using biofidelic FE models. This could enable a novel estimate of the potential reduction in hip fracture incidence resulting from a preventative treatment. In conclusion, FE simulations of have the potential to improve both diagnostic and treatment technology, which could lead to meaningful change in the current standard of care for individuals predisposed to hip fracture.

Hüftfrakturen in der westlichen Gesellschaft werden auf Grund der demographischen Alterung der Bevölkerung immer häufiger. Während die wirtschaftlichen und sozialen Kosten steigen, ist die Zahl der Patienten, die wegen Brüchigkeit ihrer Knochen behandelt werden, rückläufig. Ursache dieses Trends ist die Unzuverlässigkeit der klinischen Diagnosemittel, die verwendet werden um das Hüftfrakturrisko eines Patienten zu ermitteln. Kliniker verlassen sich Allgemeinen auf eine zweidimensionale Analyse der Knochenmineraldichte (aBMD), aber diese Messung ist mit einer geringen Sensitivität und einer hohen überdiagnose verbunden. Eine Schlieung der daraus resultierenden Behandlungslücke erfordert neue diagnostische Werkzeuge, die eine überlegene Sensitivität und Spezifität aufweisen. Eine alternative Diagnosemethode sind Festigkeitsberechnungen mit Finite Elemente (FE) Modellen, die auf Computertomographie (CT) basieren. Obwohl diese Modelle, wenn verglichen mit der Bruchfestigkeit des Femurs in Kadaverexperimenten, bessere Ergebnisse als Knochendichtemessungen erzielen, konnten mit diesem Ansatz keine verbesserte Diagnosemethode zur Vorhersage von Hüftfrakturen entwickelt werden.

Die Motivation dieser Doktorarbeit ist die Beurteilung der Grenzen der Vorhersagefähigkeit von Hüftfrakturen mit FE-Modellen, um letztendlich ein verbessertes Vorhersagemodell für Hüftfrakturen zu entwickeln. Hüftfrakturen sind das Ergebnis eines Falls auf die Hüfte, in Kombination mit einem Hüftschlag, der schwere als die Belastbarkeit des Femurs ist. Diese drei Komponenten der Hüftfraktur (Sturzrisiko, Schwere des Falls und Knochenstärke) wurden typischerweise einzeln untersucht, wobei in Wirklichkeit ein Verständnis aller drei Faktoren erforderlich ist, um die tatsächliche Risiko der Hüftfraktur zu messen, und Personen mit hohem Risiko zu identifizieren. Diese Arbeit ist in drei Ziele unterteilt, um die Biofidelity von FESimulationen zu verbessern:​

In der ersten Studie wurde die Anisotropie und Direktionalität von trabekulärer Knochen in zahlreichen Teilvolumen gemessen, die aus hochaufgelösten CT-Scans extrahiert wurden. Dabei wurde die etablierte Methode der ”mean intercept length” verwendet. Die Ergebnisse wurden dann auf FE-Modelle des proximalen Femurs übertragen, die die Bruchfestigkeit dieser Modelle, unter Randbedingung, die einen Fall zur Seite abbilden, simulieren. In dieser Studie wurden zum ersten Mal Modelle mit isotropischen und anisotropischen Materaleigenschaften und verschieden Mapping-Strategien verglichen, die die Extreme früher veröffentlichter Modellierungstechniken darstellen. Die resultierende Gesamtknochensteifigkeit und Dehnungen an der Oberfläche wurden dann mit Kadaverexperimenten verglichen, die unter ähnlichen Randbedingungen getestet wurden. Die Zugabe von Anisotropie hatte nur sehr geringe Auswirkungen auf die Knochensteifigkeit und nur eine geringe Auswirkung auf die resultierenden Dehnungen an der Oberfläche. Unterschiede in den Dehnungen entlang der Hauptspannungsrichtungen wurden im Femurkopf, Femurhals und im Trochanter Major festgestellt. Das Ergebnis der Studie war, dass anisotrope Materialeigenschaften, die mit morphologischen Messungen von hochauflösenden CT-Scans zugeordnet wurden, wenig Einfluss auf die makroskopischen mechanischen Eigenschaften auf Organebene haben. Wenn lokale interne Dehnungen von Interesse sind, könnte die Anisitropie jedoch relevant sein.

In der zweiten Studie wurde eine neuartige Materialzuordnungsverfahren für dynamische FE-Modelle des proximalen Femurs vorgestellt, die mit Fallturmversuchen validiert wurden. Nichtlinearen Materialeigenschaften wurden entwickelt, um großen Verformungen Rechnung zu tragen, die während eines dynamischen Aufpralls auftreten werden. Diese beinhalteten Zugschädigung, Druckversagen und Druckverdichtung als auch Zug-Druck-Asymmetrie und Dehnratenabhängigkeit. Die ermittelte Korrelation zwischen dynamischer Femurfestigkeiten im Versuch und in der Simulation, wurder noch nicht zuvor bestimmt und stellt damit den Benchmark der dynamischen Femurmodellierung für Fallaufprallsimulationen dar. Darüber hinaus korrelierte die von den dynamischen Modellen vorhergesagte Impulsantwort stark mit der Kraft-Zeit Antwort der Falltrumversuche. Druckdehnungsraten über 100/s traten in Elementen, im Femurkopf, Femurhals und Trochanter Major auf, was nahelegt, dass ein gutes Verständis der Materialeigenschaften von Knochengewebe für Belastungen bis zu diesen Dehnraten erforderlich ist. Die dynamischen Modelle waren in der Lage, Bruchbilder von Sub-Kapital-Frakturen waren jedoch nicht, intertrochantäre Frakturen vorherzusagen, deshlab ist anzunehmen, dass die Materialzuordnungsverfahren in dieser anatomischen Region noch verbessert werden könnte.

In der dritten Studie wurde eine neuartige biofidelic Modellierungstechnik vorgestellt, die Modelle des Beckens, des Weichgewebes und der untere Extremitäten beinhaltete. Die Massen, Geometrien und Trägheitmomente dieser Körperteile basierten auf patientenspezifischer Biometrie. Diese Modellierungstechnik wurde in einer groß angelegten, klinischen retrospektiven Studie getestet, um zu bestimmen, ob Frakturvorhersagemodelle, die auf Parametern aus biofidelic Modellen basieren, eine bessere Vorhersage des Erleidens einer Hüftfrakture, verglichen zum klinischen Standard liefern. Makroskopische Messungen der Femurfestigkeit und der Parameter für Dehnungen an der Oberfläche ergaben nur eine geringfügig bessere Sensitivität und Spezifität im Vergleich zu aBMD, was darauf hindeutet, dass wichtige patientenspezifische Informationen im Vorhersagemodell noch fehlen. Die FE-basierten Kriterien übertrafen aBMD deutlicher im hochspezifischen Teil der Recieve-OporatorKennlinie. Die Berücksichtigung von Eigenangaben zur Fallhäufigkeit verbesserte die Vorhersagekraft des Modelles nicht. Dies weist darauf hin, dass bessere Fallrisikoindikatoren für Kohortenstudien notwendig sind.

Die Dissertation leistet mehrere wichtige Beiträge zur Bewertung des Hüftbruchrisikos basierend auf FE-Simulationen. Der erste Teil veranschaulichte, dass gewisse Gewebeeigenschaften, wie z.B. Anisotropie, im Vergleich zur Methode der heterogenen Materialeigenschaftszuordnung, keine wesentlichen Rolle für die Mechanik auf Organebene spielen. Des Weiteren, konnten Modelle des proximalen Femurs auf Organebene mit einer sorgfältig entwickelten Strategie zur heterogenen Materialzuordnung, zum ersten Mal die Festigkeit und Impulsantwort unter Stoßbelastung vorhersagen. Darüber hinaus, diese Forschung präsentierte einen neuen Ansatz zur Vorhersage des Frakturrisikos, welcher das Sturzrisiko, die Heftigkeit des Aufpralls und die Knochenfestigkeit verbindet. Anschlieend wurde die Methode in einer gro angelegten retrospektiven Studie and klinischen Daten getested. Obwohl einige Einschränkungen bei der Verwendung dieser Methode berücksichtigt werden sollte, kann diese Modellierungs-Pipeline als Versuchsplattform dienen. Durch die Implementierung zusätzlicher patientenspezifischer Informationen kann deren Potential zur Verbessung der Vorhersage von Hüftfrakturen in zukünftigen Studien bewertet werden. Dies gilt ins Besondere für die Heftigkeit des Aufpralls und das Sturzrisiko. Obwohl einige Einschränkungen bei der Verwendung dieser Methode berücksichtigt werden sollte, kann diese Modellierungs-Pipeline als Versuchsplattform dienen. Durch die Implementierung zusätzlicher patientenspezifischer Informationen kann deren Potential zur Verbessung der Vorhersage von Hüftfrakturen in zukünftigen Studien bewertet werden. Dies gilt ins Besondere für die Heftigkeit des Aufpralls und das Sturzrisiko. Letztendlich, hat sich der derzeitige Ansatz zur biofidelic Modellierung als mindestens ebenso wirksam erwiesen, wie der derzeitige diagnostische Standard. Folglich können die Modelle verwendet werden, um vorbeugende Manahmen, die darauf abzielen die Häufigkeit von Hüftfrakturen in der Gesellschaft zu reduzieren, auf ihre Wirksamkeit zu untersuchen.