Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

Experimental validation of intact and implanted distal femur finite element models

Four finite element (FE) models of intact and distal femur of knee replacements were validated relative to measured bone strains. FE models of linear tetrahedrons were used. Femoral replacements with cemented stemless, cemented and noncemented femoral stems of the PFC Sigma Modular Knee System were analyzed. Bone strains were recorded at ten locations on the cortex. The magnitude of the FE bone strains corresponded to the mean measured strains, with an overall agreement of 10%. Linear regression between the FE and mean experimental strains produced slopes between 0.94 and 1.06 and R(2) values between 0.92 and 0.99. RSME values were less than 12%. The FE models were able to adequately replicate the mechanical behavior of distal femur reconstructions.     

Modelling the fibrous tissue layer in cemented hip replacements: experimental and finite element methods

The long-term fixation of cemented femoral components may be jeopardised by the presence of a fibrous tissue layer at the bone-cement interface. This study used both experimental and finite element (FE) methods to investigate the load transfer characteristics of two types of cemented hip replacements (Lubinus SPII and Müller-Curved) with a fibrous tissue layer.The experimental part investigated six stems of each type, where these were implanted in composite femurs with a specially selected silicone elastomer modelling the soft interfacial layer. Two fibrous tissue conditions were examined: a layer covering the full cement mantle, representing a revision condition; and a layer covering the proximal portion of the cement mantle, representing a non-revised implant with partial debonding and fibrous tissue formation. The FE method was used to model the full fibrous tissue layer condition, for both implants. The layer was modelled as a homogeneous, linearly isotropic material. A cross-comparison was performed of the experimental and FE findings.Agreement between experimental and FE models was verified to be within 15%. Varying the stiffness parameter of the FE soft tissue layer had little influence on the cortical bone strains, though had considerable effect on the cement strains. Stress shielding occurred for both stems under both fibrous tissue conditions, with the greatest reduction around the calcar. However, the cortical bone strains were generally larger than those for the equivalent well-fixed stems. The fibrous tissue layer was not found to increase the general strain pattern of the cement mantle, though localised regions of high stress were detected.     

Cement mantle fatigue failure in total hip replacement: experimental and computational testing

One possible loosening mechanism of the femoral component in total hip replacement is fatigue cracking of the cement mantle. A computational method capable of simulating this process may therefore be a useful tool in the preclinical evaluation of prospective implants. In this study, we investigated the ability of a computational method to predict fatigue cracking in experimental models of the implanted femur construct. Experimental specimens were fabricated such that cement mantle visualisation was possible throughout the test. Two different implant surface finishes were considered: grit blasted and polished. Loading was applied to represent level gait for two million cycles. Computational (finite element) models were generated to the same geometry as the experimental specimens, with residual stress and porosity simulated in the cement mantle. Cement fatigue and creep were modelled over a simulated two million cycles. For the polished stem surface finish, the predicted fracture locations in the finite element models closely matched those on the experimental specimens, and the recorded stem displacements were also comparable. For the grit blasted stem surface finish, no cement mantle fractures were predicted by the computational method, which was again in agreement with the experimental results. It was concluded that the computational method was capable of predicting cement mantle fracture and subsequent stem displacement for the structure considered.     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Mechanical response comparison in an implant overdenture retained by ball attachments on conventional regular and mini dental implants: a finite element analysis

This study investigates the bone/implant mechanical responses in an implant overdenture retained by ball attachments on two conventional regular dental implants (RDI) and four mini dental implants (MDI) using finite element (FE) analysis. Two FE models of overdentures retained by RDIs and MDIs for a mandibular edentulous patient with validation within 6% variation errors were constructed by integrating CT images and CAD system. Bone grafting resulted in 2 mm thickness at the buccal side constructed for the RDIs-supported model to mimic the bone augmentation condition for the atrophic alveolar ridge. Nonlinear hyperelastic material and frictional contact element were used to simulate characteristic of the ball attachment-retained overdentures. The results showed that a denture supported by MDIs presented higher surrounding bone strains than those supported by RDIs under different load conditions. Maximum bone micro strains were up to 6437/2987 and 13323/5856 for MDIs/RDIs under single centric and lateral contacts, respectively. Corresponding values were 4429/2579 and 9557/5774 under multi- centric and lateral contacts, respectively. Bone micro strains increased 2.06 and 1.96-folds under single contact, 2.16 and 2.24-folds under multiple contacts for MDIs and RDIs when lateral to axial loads were compared. The maximum RDIs and MDIs implant stresses in all simulated cases were found by far lower than their yield strength. Overdentures retained using ball attachments on MDIs in poor edentulous bone structure increase the surrounding bone strain over the critical value, thereby damaging the bone when compared to the RDIs. Eliminating the occlusal single contact and oblique load of an implant-retained overdenture reduces the risk for failure.     

Sensitivity of periprosthetic stress-shielding to load and the bone density-modulus relationship in subject-specific finite element models

Subject-specific finite element (FE) computer models of the proximal femur in hip replacement could potentially predict stress-shielding and subsequent bone loss in individual patients. Before such predictions can be made, it is important first to determine if between subject differences in stress-shielding are sensitive to poorly defined parameters such as the load and the bone material properties. In this study we investigate if subject-specific FE models provide consistent stress-shielding patterns in the bone, independent of the choice of the loading conditions and the bone density-modulus relationship used in the computer model. FE models of two right canine femurs with and without implants were constructed based on contiguous computed tomography (CT) scans so that subject-specific estimates of stress-shielding could be calculated. Four different loading conditions and two bone density-modulus relationships were tested. Stress-shielding was defined as the decrease of strain energy per gram bone mass in the femur with the implant in place relative to the intact femur.The analyses showed that for the four loading conditions and two bone density-modulus relationships the difference in stress-shielding between the two subjects was essentially constant (1% variation) when the same loading condition and density-modulus relationship was used for both subjects. The severity of stress-shielding within a subject was sensitive to these input parameters, varying up to 20% in specific regions with a change in loading conditions and up to 10% for a change in the assumed density-modulus relationship. We conclude that although the choice of input parameters can substantially affect stress-shielding in an individual, this choice had virtually no effect on the relative differences in femoral periprosthetic stress-shielding between individuals. Thus, while care should be taken in the interpretation of the absolute value of stress-shielding calculated with these type of models, subject-specific FE models may be useful for explaining the variation in bone adaptation responsiveness between different subjects in experimental or clinical studies.     

Validation of an efficient method of assigning material properties in finite element analysis of pelvic bone

Bone in the pelvis is a composite material with a complex anatomical structure that is difficult to model computationally. Rather than assigning material properties to increasingly smaller elements to capture detail in three-dimensional finite element (FE) models, properties can be assigned to Gauss points within larger elements. As part of a validation process, we compared experimental and analytical results from a composite beam under four-point load to FE models with material properties assigned to refined elements and Gauss points within larger elements. Both FE models accurately predicted deformation and the analytical predictions of internal shear stress.     

Strain and micromotion in intact and resurfaced composite femurs: Experimental and numerical investigations

Understanding the load transfer within a resurfaced femur is necessary to determine the influence of mechanical factors on potential failure mechanisms such as early femoral neck fractures and stress shielding. In this study, an attempt has been made to measure the stem-bone micromotion and implant cup-bone relative displacements (along medial-lateral and anterior-posterior direction), in addition to surface strains at different locations and orientations on the proximal femur and to compare these measurements with those predicted by equivalent FE models. The loading and the support conditions of the experiment were closely replicated in the FE models. A new experimental set-up has been developed, with specially designed fixtures and load application mechanism, which can effectively impose bending and deflection of the tested femurs, almost in any direction. High correlation coefficient (0.92–0.95), low standard error of the estimate (170–379 με) and low percentage error in regression slope (12.8–17.5%), suggested good agreement between the numerical and measured strains. The effect of strain shielding was observed in two (out of eight) strain gauges located on the posterior side. A pronounced strain increase occurred in strain gauges located on the anterior head and neck regions after implantation. Experimentally measured stem-bone micromotion and implant cup-bone relative displacements (0–13.7 μm) were small and similar in trends predicted by the FE models (0–25 μm). Despite quantitative deviations in the measured and numerical results, it appears that the FE model can be used as a valid predictor of the actual strain and stem-bone micromotion.     

Comparative in vitro study on the long term performance of cemented hip stems: validation of a protocol to discriminate between "good" and "bad" designs

The long-term clinical success of cemented hip stems is influenced both by the implant design, and by the surgical procedure. A methodology is proposed for discriminating between implant designs with different clinical outcomes. The protocol was designed with industrial pre-clinical validation in mind.Two cemented stem types were tested, one (Lubinus SPII) having good and the other (Müller Curved) having poor clinical outcomes. Three implants for each type were subjected to a mechanical in vitro test of one million loading cycles. Each cycle reproduced the load components of stair climbing. Interface shear micromotion was measured during the test in the direction of rotation and along the stem axis. The stem roughness before and after the test was compared. After the test, the cement mantles were retrieved and inspected through dye penetrants to detect evidences of micro-damage. For each specimen, the events of the loosening process were examined, based on the in vitro data available, so as to analyze the whole failure mechanism.The protocol developed was sensitive to the implant design, with significantly different results being found for the two stem types, both in terms of stem-cement micromotions, surface roughness alteration, and cement mantle damage. The information yielded by the three different investigation techniques was consistent for each of the two groups of specimens tested, allowing a better understanding of the failure process. In vitro inducible micromotion and permanent migration measurements, together with cement-stem interface fretting damage and cement fatigue damage, can help predicting the clinical performance of cemented stems.     

On the importance of considering porosity when simulating the fatigue of bone cement

Fatigue cracking in the cement mantle of total hip replacement has been identified as a possible cause of implant loosening. Retrieval studies and in vitro tests have found porosity in the cement may facilitate fatigue cracking of the mantle. The fatigue process has been simulated computationally using a finite element/continuum damage mechanics (FE/CDM) method and used as a preclinical testing tool, but has not considered the effects of porosity. In this study, experimental tensile and four-point bend fatigue tests were performed. The tensile fatigue S-N data were used to drive the computational simulation (FE/CDM) of fatigue in finite element models of the tensile and four-point bend specimens. Porosity was simulated in the finite element models according to the theory of elasticity and using Monte Carlo methods. The computational fatigue simulations generated variability in the fatigue life at any given stress level, due to each model having a unique porosity distribution. The fracture site also varied between specimens. Experimental validation was achieved for four-point bend loading, but only when porosity was included. This demonstrates that the computational simulation of fatigue, driven by uniaxial S-N data can be used to simulate nonuniaxial loadcases. Further simulations of bone cement fatigue should include porosity to better represent the realities of experimental models.     

Experimental validation of finite element models of intact and implanted composite hemipelvises using digital image correlation

A detailed understanding of the changes in load transfer due to implantation is necessary to identify potential failure mechanisms of orthopedic implants. Computational finite element (FE) models provide full field data on intact and implanted bone structures, but their validity must be assessed for clinical relevance. The aim of this study was to test the validity of FE predicted strain distributions for the intact and implanted pelvis using the digital image correlation (DIC) strain measurement technique. FE models of an in vitro hemipelvis test setup were produced, both intact and implanted with an acetabular cup. Strain predictions were compared to DIC and strain rosette measurements. Regression analysis indicated a strong linear relationship between the measured and predicted strains, with a high correlation coefficient (R?=?0.956 intact, 0.938 implanted) and a low standard error of the estimate (SE?=?69.53?με, 75.09?με). Moreover, close agreement between the strain rosette and DIC measurements improved confidence in the validity of the DIC technique. The FE model therefore was supported as a valid predictor of the measured strain distribution in the intact and implanted composite pelvis models, confirming its suitability for further computational investigations.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.     

Polymerisation stress modelling in acrylic bone cement

Fatigue failure of the cement mantle has been proposed as one of the failure processes contributing to aseptic loosening of cemented joint replacements. It has also been suggested that fatigue failure is dramatically accelerated by residual stress generated during the cement polymerisation process. Previous computational models of the polymerisation process have investigated only the latter part of polymerisation by assuming both instantaneous hardening of the material (a stress locking point) and that all residual stress results from thermal shrinkage after this stress locking point. In this study, finite element models which use the local degree of polymerisation to calculate material properties and shrinkage have been used to predict residual stresses in two models of total hip replacement cement mantles. Results indicate that the final value of cement mantle stress may not be the highest stresses that the cement is subjected to during the polymerisation process. Two models are presented, a 2-dimensional model, which was adapted from a similar model in the literature (Lennon and Prendergast, 2002) and a 3-dimensional concentric-cylinders model. In both cases a chemical kinetics model was used to predict the progress of the polymerisation reaction and a second linear model used to predict cement mechanical properties and density, and so stress generation and volume change, over time. There was good agreement of the results of the 2D model with its counterpart in the literature. For the 3D model, the final residual stress magnitudes and patterns showed good agreement with similar physical and computational models in the literature.     

3D characterization of bone strains in the rat tibia loading model

Bone strain is considered one of the factors inducing bone tissue response to loading. Nevertheless, where animal studies can provide detailed data on bone response, they only offer limited information on experimental bone strains. Including micro-CT-based finite element (micro FE) models in the analysis represents a potent methodology for quantifying strains in bone. Therefore, the main objective of this study was to develop and validate specimen-specific micro FE models for the assessment of bone strains in the rat tibia compression model. Eight rat limbs were subjected to axial compression loading; strain at the medio-proximal site of the tibiae was measured by means of strain gauges. Specimen-specific micro FE models were created and analyzed. Repeated measurements on each limb indicated that the effect of limb positioning was small (COV?= 6.45 ± 2.27 %). Instead, the difference in the measured strains between the animals was high (54.2%). The computational strains calculated at the strain gauge site highly correlated to the measured strains (R ²?=?0.95). Maximum peak strains calculated at exactly 25% of the tibia length for all specimens were equal to 435.11 ± 77.88 microstrains (COV?=?17.19%). In conclusion, we showed that strain gauge measurements are very sensitive to the exact strain gauge location on the bone; hence, the use of strain gauge data only is not recommended for studies that address at identifying reliable relationships between tissue response and local strains. Instead, specimen-specific micro FE models of rat tibiae provide accurate estimates of tissue-level strains.     

A 1-year study of osteoinduction in hydroxyapatite-derived biomaterials in an adult sheep model: part II. Bioengineering implants to optimize bone replacement in reconstruction of cranial defects

The present study investigated hydroxyapatite biomaterials implanted in critical-size defects in the calvaria of adult sheep to determine the optimal bioengineering of hydroxyapatite composites to facilitate bone ingrowth into these materials. Five calvarial defects measuring 16.8 mm in diameter were made in each of 10 adult sheep. Three defects were filled with cement paste composites of hydroxyapatite and beta-tricalcium phosphate as follows: (1) 100 percent hydroxyapatite-cement paste, (2) 60 percent hydroxyapatite-cement paste, and (3) 20 percent hydroxyapatite-cement paste. One defect was filled with a ceramic composite containing 60 percent hydroxyapatite-ceramic, and the fifth defect remained unfilled. One year after implantation, the volume of all biomaterials was determined by computed tomography, and porosity and bone replacement were determined using backscatter electron microscopy. Computed tomography-based volumetric assessment 1 year after implantation demonstrated that none of the unfilled cranial defects closed over the 1-year period, confirming that these were critical-size defects. There was a significant increase in volume in both the cement paste and ceramic implants containing 60 percent hydroxyapatite (p < 0.01). There was no significant change in volume of the remaining cement paste biomaterials. Analysis of specimens by backscatter electron microscopy demonstrated mean bone replacement of 4.8 +/- 1.4 percent (mean +/- SEM) in 100 percent hydroxyapatite-cement paste, 11.2 +/- 2.3 percent in 60 percent hydroxyapatite-cement paste, and 28.5 +/- 4.5 percent in 20 percent hydroxyapatite-cement paste. There was an inverse correlation between the concentration of hydroxyapatite and the amount of bone replacement in the cement paste for each composite tested (p < 0.01). Bone replacement in 60 percent hydroxyapatite-ceramic composite (13.6 +/- 2.0 percent) was not significantly different from that in 60 percent hydroxyapatite-cement paste. Of note is that the ceramic composite contained macropores (200 to 300 microm) that did not change in size over the 1-year period. All cement paste composites initially contained micropores (3 to 5 nm), which remained unchanged in 100 percent hydroxyapatite-cement paste. Cement paste implants containing increased tricalcium phosphate demonstrated a corresponding increase in macropores following resorption of the tricalcium phosphate component. Bone replacement occurred within the macropores of these implants. In conclusion, there was no significant bone ingrowth into pure hydroxyapatite-cement paste (Bone Source, Stryker-Leibinger Inc., Dallas, Texas) in the present study. The introduction of macropores in a biomaterial can optimize bone ingrowth for reconstruction of critical-size defects in calvaria. This was demonstrated in both the ceramic composite of hydroxyapatite tested and the cement paste composites of hydroxyapatite by increasing the composition of a rapidly resorbing component such as beta-tricalcium phosphate.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

Biomechanical validation of finite element models for two silicone metacarpophalangeal joint implants

Silicone implants are used for prosthetic arthroplasty of metacarpophalangeal (MCP) joints severely damaged by rheumatoid arthritis. Different silicone elastomer MCP implant designs have been developed, including the Swanson and the NeuFlex implants. The goal of this study was to compare the in vitro mechanical behavior of Swanson and NeuFlex MCP joint implants. Three-dimensional (3D) finite element (FE) models of the silicone implants were modeled using the commercial software ANSYS and subjected to angular displacement from 0 deg to 90 deg. FE models were validated using mechanical tests of implants incrementally bent from 0 deg to 90 deg in a joint simulator. Swanson size 2 and 4 implants were compared with NeuFlex size 10 and 30 implants, respectively. Good agreement was observed throughout the range of motion for the flexion bending moment derived from 3D FE models and mechanical tests. From 30 deg to 90 deg, the Swanson 2 demonstrated a greater resistance to deformation than the NeuFlex 10 and required a greater bending moment for joint flexion. For larger implant sizes, the NeuFlex 30 had a steeper moment-displacement curve, but required a lower moment than the Swanson 4, due to implant preflexion. On average, the stress generated at the implant hinge from 30 deg to 90 deg was lower in the NeuFlex than in the Swanson. On average, starting from the neutral position of 30 deg for the preflexed NeuFlex implant, higher moments were required to extend the NeuFlex implants to 0 deg compared with the Swanson implants, which returned spontaneously to resting position. Implant toggling within the medullary canals was less in the NeuFlex than in the Swanson. The differential performance of these implants may be useful in implant selection based on the preoperative condition(s) of the joint and specific patient functional needs.     

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.     

Patient-specific finite-element analyses of the proximal femur with orthotropic material properties validated by experiments

Patient-specific high order finite-element (FE) models of human femurs based on quantitative computer tomography (QCT) with inhomogeneous orthotropic and isotropic material properties are addressed. The point-wise orthotropic properties are determined by a micromechanics (MM) based approach in conjunction with experimental observations at the osteon level, and two methods for determining the material trajectories are proposed (along organs outer surface, or along principal strains). QCT scans on four fresh-frozen human femurs were performed and high-order FE models were generated with either inhomogeneous MM-based orthotropic or empirically determined isotropic properties. In vitro experiments were conducted on the femurs by applying a simple stance position load on their head, recording strains on femurs' surface and head's displacements. After verifying the FE linear elastic analyses that mimic the experimental setting for numerical accuracy, we compared the FE results to the experimental observations to identify the influence of material properties on models' predictions. The strains and displacements computed by FE models having MM-based inhomogeneous orthotropic properties match the FE-results having empirically based isotropic properties well, and both are in close agreement with the experimental results. When only the strains in the femoral neck are being compared a more pronounced difference is noticed between the isotropic and orthotropic FE result. These results lay the foundation for applying more realistic inhomogeneous orthotropic material properties in FEA of femurs.