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.     

Scaling of form and function in the xenarthran femur: a 100-fold increase in body mass is mitigated by repositioning of the third trochanter

How animals cope with increases in body size is a key issue in biology. Here, we consider scaling of xenarthrans, particularly how femoral form and function varies to accommodate the size range between the 3 kg armadillo and its giant relative the 300 kg glyptodont. It has already been noted that femoral morphology differs between these animals and suggested that this reflects a novel adaptation to size increase in glyptodont. We test this idea by applying a finite element analysis of coronal plane forces to femoral models of these animals, simulating the stance phase in the hind limb; where the femur is subject to bending owing to longitudinal compressive as well as abduction loads on the greater trochanter. We use these models to examine the hypothesis that muscles attaching on the third trochanter (T3) can reduce this bending in the loaded femur and that the T3 forces are more effective at reducing bending in glyptodont where the T3 is situated at the level of the knee. The analysis uses traditional finite element methods to produce strain maps and examine strains at 200 points on the femur. The coordinates of these points before and after loading are also used to carry out geometric morphometric (GM) analyses of the gross deformation of the model in different loading scenarios. The results show that longitudinal compressive and abductor muscle loading increases bending in the coronal plane, and that loads applied to the T3 reduce that bending. In the glyptodont model, the T3 loads are more effective and can more readily compensate for the bending owing to longitudinal and abductor loads. This study also demonstrates the usefulness of GM methods in interpreting the results of finite element analyses.     

Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis.

In Part I we reported the results of linear finite element models of the proximal femur generated using geometric and constitutive data collected with quantitative computed tomography. These models demonstrated excellent agreement with in vitro studies when used to predict ultimate failure loads. In Part II, we report our extension of those finite element models to include nonlinear behavior of the trabecular and cortical bone. A highly nonlinear material law, originally designed for representing concrete, was used for trabecular bone, while a bilinear material law was used for cortical bone. We found excellent agreement between the model predictions and in vitro fracture data for both the onset of bone yielding and bone fracture. For bone yielding, the model predictions were within 2 percent for a load which simulated one-legged stance and 1 percent for a load which simulated a fall. For bone fracture, the model predictions were within 1 percent and 17 percent, respectively. The models also demonstrated different fracture mechanisms for the two different loading configurations. For one-legged stance, failure within the primary compressive trabeculae at the subcapital region occurred first, leading to load transfer and, ultimately, failure of the surrounding cortical shell. However, for a fall, failure of the cortical and trabecular bone occurred simultaneously within the intertrochanteric region. These results support our previous findings that the strength of the subcapital region is primarily due to trabecular bone whereas the strength of the intertrochanteric region is primarily due to cortical bone.     

Fracture prediction for the proximal femur using finite element models: Part I--Linear analysis.

Over 90 percent of the more than 250,000 hip fractures that occur annually in the United States are the result of falls from standing height. Despite this, the stresses associated with femoral fracture from a fall have not been investigated previously. Our objectives were to use three-dimensional finite element models of the proximal femur (with geometries and material properties based directly on quantitative computed tomography) to compare predicted stress distributions for one-legged stance and for a fall to the lateral greater trochanter. We also wished to test the correspondence between model predictions and in vitro strain gage data and failure loads for cadaveric femora subjected to these loading conditions. An additional goal was to use the model predictions to compare the sensitivity of several imaging sites in the proximal femur which are used for the in vivo prediction of hip fracture risk. In this first of two parts, linear finite element models of two unpaired human cadaveric femora were generated. In Part II, the models were extended to include nonlinear material properties for the cortical and trabecular bone. While there was poor correspondence between strain gage data and model predictions, there was excellent agreement between the in vitro failure data and the linear model, especially using a von Mises effective strain failure criterion. Both the onset of structural yielding (within 22 and 4 percent) and the load at fracture (within 8 and 5 percent) were predicted accurately for the two femora tested. For the simulation of one-legged stance, the peak stresses occurred in the primary compressive trabeculae of the subcapital region.(ABSTRACT TRUNCATED AT 250 WORDS)     

Biomechanical consequences of subtalar joint arthroereisis in treating posterior tibial tendon dysfunction: a theoretical analysis using finite element analysis

Subtalar joint arthroereisis (SJA) has been introduced to control the hyperpronation in cases of flatfoot. The objective of this study is to evaluate the biomechanical consequence of SJA to restore the internal stress and load transfer to the intact state from the attenuated biomechanical condition induced by posterior tibial tendon dysfunction (PTTD). A three-dimensional finite element model of the foot and ankle complex was constructed based on clinical images of a healthy female (age 28 years, height 165 cm, body mass 54 kg). The boundary and loading condition during walking was acquired from the gait experiment of the model subject. Five sets of simulations (conditions) were completed: intact condition, mild PTTD, severe PTTD, mild PTTD with SJA, severe PTTD with SJA. The maximum von Mises stress of the metatarsal shafts and the load transfer along the midfoot during stance were analyzed. Generally, SJA deteriorated the joint force of the medial cuneonavicular and calcaneocuboid joints during late stance, while that of the metatarsocuneiform joints during early stance were over-corrected. Only the calcaneocuboid joint force at 45% stance demonstrated a trend of improvement. Besides, SJA exaggerated the increased stress of the metatarsals compared to the PTTD conditions, except that of the first metatarsal. Our study did not support the hypothesis that SJA can restore the internal load transfer and midfoot stress. SJA cannot compensate the salvage of midfoot stability attributed by PTTD and could be biomechanically insufficient to restore the biomechanical environment. Additional procedures such as orthotic intervention may be necessary.     

Combined finite element model of human proximal femur behaviour considering remodeling and fracture

The purpose of this work was to develop a combined remodeling-to-fracture finite element model allowing for the combined simulation of human proximal femur remodeling under a given boundary conditions followed by the simulation of its fracture behaviour under quasi-static load. The combination of remodeling and fracture simulation into one unified model consists in considering that the femur properties resulting from the remodeling simulation correspond to the initial state for the fracture prediction. The remodeling model is based on a coupled strain and fatigue damage stimulus approach. The fracture model is based on continuum damage mechanics in order to predict the progressive fracturing process, which allows to predict the fracture pattern and the complete force-displacement curve under quasi-static load. To investigate the potential of the proposed unified remodeling-to-fracture model, we performed remodeling simulations on a 3D proximal femur model for a duration of 365 days followed by a side fall fracture simulation reproducing.     

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

Proximal femur strength estimates from computed tomography (CT)-based finite element (FE) models are finding clinical application. Published models reached a high in-vitro accuracy, yet many of them rely on nonlinear methodologies or internal best-fitting of parameters. The aim of the present study is to verify to what extent a linear FE modelling procedure, fully based on independently determined parameters, can predict the failure characteristics of the proximal femur in stance and sideways fall loading configurations.     

Prediction of strength and strain of the proximal femur by a CT-based finite element method

Hip fractures are the most serious complication of osteoporosis and have been recognized as a major public health problem. In elderly persons, hip fractures occur as a result of increased fragility of the proximal femur due to osteoporosis. It is essential to precisely quantify the strength of the proximal femur in order to estimate the fracture risk and plan preventive interventions. CT-based finite element analysis could possibly achieve precise assessment of the strength of the proximal femur. The purpose of this study was to create a simulation model that could accurately predict the strength and surface strains of the proximal femur using a CT-based finite element method and to verify the accuracy of our model by load testing using fresh frozen cadaver specimens. Eleven right femora were collected. The axial CT scans of the proximal femora were obtained with a calibration phantom, from which the 3D finite element models were constructed. Materially nonlinear finite element analyses were performed. The yield and fracture loads were calculated, while the sites where elements failed and the distributions of the principal strains were determined. The strain gauges were attached to the proximal femoral surfaces. A quasi-static compression test of each femur was conducted. The yield loads, fracture loads and principal strains of the prediction significantly correlated with those measured (r=0.941, 0.979, 0.963). Finite element analysis showed that the solid elements and shell elements in undergoing compressive failure were at the same subcapital region as the experimental fracture site.     

A model of bone adaptation as an optimization process

The internal bone adaptation of the proximal femur is considered. A three-dimensional finite element model of the proximal femur is used. The bone remodeling in this work is numerically described by an evolutionary remodeling scheme with anisotropic material parameters and time-dependent loading. The memory of past loading is included in the model to account for the delay in the bone response from the load changes. The remodeling rate equation is derived from the structural optimization task of maximizing the stiffness in each time step. Additional information can be extracted from the optimization process; the remodeling equilibrium parameter where no apposition or resorption takes place, is defined as the element optimality conditions and the optimal design is used as an initial design for the onset of the remodeling simulation. Two examples of bone adaptation resulting from load changes are given, and the irreversible nature of the model is illustrated.     

Automated three-dimensional finite element modelling of bone: a new method

Three-dimensional finite element stress analysis of bone is a key to understanding bone remodelling, assessing fracture risk, and designing prostheses; however, the cost and complexity of predicting the stress field in bone with accuracy has precluded the routine use of this method. A new, automated method of generating patient-specific three-dimensional finite element models of bone is presented — it uses digital computed tomographic (CT) scan data to derive the geometry of the bone and to estimate its inhomogeneous material properties. Cubic elements of a user-specified size are automatically defined and then individually assigned the CT scan-derived material properties. The method is demonstrated by predicting the stress, strain, and strain energy in a human proximal femur in vivo. Three-dimensional loading conditions corresponding to the stance phase of gait were taken from the literature and applied to the model. Maximum principal compressive stresses of 8–23 MPa were computed for the medial femoral neck. Automated generation of additional finite element models with larger numbers of elements was used to verify convergence in strain energy.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

Effect of varus/valgus malalignment on bone strains in the proximal tibia after TKR: an explicit finite element study

Malalignment is the main cause of tibial component loosening. Implants that migrate rapidly in the first two post-operative years are likely to present aseptic loosening. It has been suggested that cancellous bone stresses can be correlated with tibial component migration. A recent study has shown that patient-specific finite element (FE) models have the power to predict the short-term behavior of tibial trays. The stresses generated within the implanted tibia are dependent on the kinematics of the joint; however, previous studies have ignored the kinematics and only applied static loads. Using explicit FE, it is possible to simultaneously predict the kinematics and stresses during a gait cycle. The aim of this study was to examine the cancellous bone strains during the stance phase of the gait cycle, for varying degrees of varus/valgus eccentric loading using explicit FE. A patient-specific model of a proximal tibia was created from CT scan images, including heterogeneous bone properties. The proximal tibia was implanted with a commercial total knee replacement (TKR) model. The stance phase of gait was simulated and the applied loads and boundary conditions were based on those used for the Stanmore knee simulator. Eccentric loading was simulated. As well as examining the tibial bone strains (minimum and maximum principal strain), the kinematics of the bone-implant construct are also reported. The maximum anterior-posterior displacements and internal-external rotations were produced by the model with 20 mm offset. The peak minimum and maximum principal strain values increased as the load was shifted laterally, reaching a maximum magnitude for -20 mm offset. This suggests that when in varus, the load transferred to the bone is shifted medially, and as the bone supporting this load is stiffer, the resulting peak bone strains are lower than when the load is shifted laterally (valgus). For this particular patient, the TKR design analyzed produced the highest cancellous bone strains when in valgus. This study has provided an insight in the variations produced in bone strain distribution when the axial load is applied eccentrically. To the authors' knowledge, this is the first time that the bone strain distribution of a proximal implanted tibia has been examined, also accounting for the kinematics of the tibio-femoral joint as part of the simulation. This approach gives greater insight into the overall performance of TKR.     

The importance of femur/acetabulum cartilage in the biomechanics of the intact hip: experimental and numerical assessment

Experimental studies have been made to study and validate the biomechanics of the pair femur/acetabulum considering both structures without the presence of cartilage. The main goal of this study was to validate a numerical model of the intact hip. Numerical and experimental models of the hip joint were developed with respect to the anatomical restrictions. Both iliac and femur bones were replicated based on composite replicas. Additionally, a thin layer of silicon rubber was used for the cartilage. A three-dimensional finite element model was developed and the boundary conditions of the models were applied according to the natural physiological constrains of the joint. The loads used in both models were used just for comparison purposes. The biomechanical behaviour of the models was assessed considering the maximum and minimum principal bone strains and von Mises stress. We analysed specific biomechanical parameters in the interior of the acetabular cavity and on femur's surface head to determine the role of the cartilage of the hip joint within the load transfer mechanism. The results of the study show that the stress observed in acetabular cavity was 8.3 to 9.2 MPa. When the cartilage is considered in the joint model, the absolute values of the maximum and minimum peak strains on the femur's head surface decrease simultaneously, and the strains are more uniformly distributed on both femur and iliac surfaces. With cartilage, the cortex strains increase in the medial side of the femur. We prove that finite element models of the intact hip joint can faithfully reproduce experimental models with a small difference of 7%.     

Prospects of implant with locking plate in fixation of subtrochanteric fracture: experimental demonstration of its potential benefits on synthetic femur model with supportive hierarchical nonlinear hyperelastic finite element analysis

Background

Effective fixation of fracture requires careful selection of a suitable implant to provide stability and durability. Implant with a feature of locking plate (LP) has been used widely for treating distal fractures in femur because of its favourable clinical outcome, but its potential in fixing proximal fractures in the subtrochancteric region has yet to be explored. Therefore, this comparative study was undertaken to demonstrate the merits of the LP implant in treating the subtrochancteric fracture by comparing its performance limits against those obtained with the more traditional implants; angle blade plate (ABP) and dynamic condylar screw plate (DCSP).

Materials and Methods

Nine standard composite femurs were acquired, divided into three groups and fixed with LP (n?=?3), ABP (n?=?3) and DCSP (n?=?3). The fracture was modeled by a 20?mm gap created at the subtrochanteric region to experimentally study the biomechanical response of each implant under both static and dynamic axial loading paradigms. To confirm the experimental findings and to understand the critical interactions at the boundaries, the synthetic femur/implant systems were numerically analyzed by constructing hierarchical finite element models with nonlinear hyperelastic properties. The predictions from the analyses were then compared against the experimental measurements to demonstrate the validity of each numeric model, and to characterize the internal load distribution in the femur and load bearing properties of each implant.

Results

The average measurements indicated that the constructs with ABP, DCPS and LP respectively had overall stiffness values of 70.9, 110.2 and 131.4?N/mm, and exhibited reversible deformations of 12.4, 4.9 and 4.1?mm when the applied dynamic load was 400?N and plastic deformations of 11.3, 2.4 and 1.4?mm when the load was 1000?N. The corresponding peak cyclic loads to failure were 1100, 1167 and 1600?N. The errors between the displacements measured experimentally or predicted by the nonlinear hierarchical hyperelastic model were less than 18?%. In the implanted femur heads, the principal stresses were spatially heterogeneous for ABP and DCSP but more homogenous for LP, meaning LP had lower stress concentrations.

Conclusion

When fixed with the LP implant, the synthetic femur model of the subtrochancteric fracture consistently exceeds in the key biomechanical measures of stability and durability. These capabilities suggest increased resistance to fatigue and failure, which are highly desirable features expected of functional implants and hence make the LP implant potentially a viable alternative to the conventional ABP or DCSP in the treatment of subtrochancteric femur fractures for the betterment of clinical outcome.     

A finite element model of the human knee joint for the study of tibio-femoral contact

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (ie., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.     

Numerical simulation of load-induced bone structural remodelling using stress-limit criterion

A simple and efficient numerical method for predicting the remodelling of adaptive materials and structures under applied loading was presented and implemented within a finite element framework. The model uses the trajectorial architecture theory of optimisation to predict the remodelling of material microstructure and structural organisation under mechanical loading. We used the proposed model to calculate the density distribution of proximal femur in the frontal plane. The loading considered was the hip joint contact forces and muscular forces at the attachment sites of the muscles to the bone. These forces were estimated from a separate finite element calculation using a heterogeneous three-dimensional model of the proximal femur. The density distributions obtained by this procedure has a qualitative similarity with in vivo observations. Solutions displayed the characteristic high-density channels that are evident in the Dual X-ray Absorptiometry scan. There is also evidence of the intramedullary canal, as well as low-density regions in the femoral neck. Several parametric studies were carried out to highlight the advantages of the proposed method, which includes fast convergence and low-computational cost. The potential applications of the proposed method in predicting bone structural remodelling in cancer are also briefly discussed.     

Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur

Functional adaptation of the femur has been investigated in several studies by embedding bone remodelling algorithms in finite element (FE) models, with simplifications often made to the representation of bone’s material symmetry and mechanical environment. An orthotropic strain-driven adaptation algorithm is proposed in order to predict the femur’s volumetric material property distribution and directionality of its internal structures within a continuum. The algorithm was applied to a FE model of the femur, with muscles, ligaments and joints included explicitly. Multiple load cases representing distinct frames of two activities of daily living (walking and stair climbing) were considered. It is hypothesised that low shear moduli occur in areas of bone that are simply loaded and high shear moduli in areas subjected to complex loading conditions. In addition, it is investigated whether material properties of different femoral regions are stimulated by different activities. The loading and boundary conditions were considered to provide a physiological mechanical environment. The resulting volumetric material property distribution and directionalities agreed with ex vivo imaging data for the whole femur. Regions where non-orthogonal trabecular crossing has been documented coincided with higher values of predicted shear moduli. The topological influence of the different activities modelled was analysed. The influence of stair climbing on the properties of the femoral neck region is highlighted. It is recommended that multiple load cases should be considered when modelling bone adaptation. The orthotropic model of the complete femur is released with this study.     

The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model

The finite element (FE) model of the pelvic joint is helpful for clinical diagnosis and treatment of pelvic injuries. However, the effect of an FE model boundary condition on the biomechanical behavior of a pelvic joint has not been well studied. The objective of this study was to study the effect of boundary condition on the pelvic biomechanics predictions. A 3D FE model of a pelvis using subject-specific estimates of intact bone structures, main ligaments and bone material anisotropy by computed tomography (CT) gray value was developed and validated by bone surface strains obtained from rosette strain gauges in an in vitro pelvic experiment. Then three FE pelvic models were constructed to analyze the effect of boundary condition, corresponding to an intact pelvic joint, a pelvic joint without sacroiliac ligaments and a pelvic joint without proximal femurs, respectively. Vertical load was applied to the same pelvis with a fixed prosthetic femoral stem and the same load was simulated in the FE model. A strong correlation coefficient (R(2)=0.9657) was calculated, which indicated a strong correlation between the FE analysis and experimental results. The effect of boundary condition changes on the biomechanical response depended on the anatomical location and structure of the pelvic joint. It was found that acetabulum fixed in all directions with the femur removed can increase the stress distribution on the acetabular inner plate (approximately double the original values) and decrease that on the superior of pubis (from 7 MPa to 0.6 MPa). Taking sacrum and ilium as a whole, instead of sacroiliac and iliolumber ligaments, can influence the stress distribution on ilium and pubis bone vastly. These findings suggest pelvic biomechanics is very dependent on the boundary condition in the FE model.     

Use of a statistical model of the whole femur in a large scale,multi-model study of femoral neck fracture risk

Interpatient variability is often overlooked in orthopaedic computational studies due to the substantial challenges involved in sourcing and generating large numbers of bone models. A statistical model of the whole femur incorporating both geometric and material property variation was developed as a potential solution to this problem. The statistical model was constructed using principal component analysis, applied to 21 individual computer tomography scans. To test the ability of the statistical model to generate realistic, unique, finite element (FE) femur models it was used as a source of 1000 femurs to drive a study on femoral neck fracture risk. The study simulated the impact of an oblique fall to the side, a scenario known to account for a large proportion of hip fractures in the elderly and have a lower fracture load than alternative loading approaches. FE model generation, application of subject specific loading and boundary conditions, FE processing and post processing of the solutions were completed automatically. The generated models were within the bounds of the training data used to create the statistical model with a high mesh quality, able to be used directly by the FE solver without remeshing. The results indicated that 28 of the 1000 femurs were at highest risk of fracture. Closer analysis revealed the percentage of cortical bone in the proximal femur to be a crucial differentiator between the failed and non-failed groups. The likely fracture location was indicated to be intertrochantic. Comparison to previous computational, clinical and experimental work revealed support for these findings.