Andrew E. Anderson

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.

Validation of finite element predictions of cartilage contact pressure in the human hip joint.

Methods to predict contact stresses in the hip can provide an improved understanding of load distribution in the normal and pathologic joint. The objectives of this study were to develop and validate a three-dimensional finite element (FE) model for predicting cartilage contact stresses in the human hip using subject-specific geometry from computed tomography image data, and to assess the sensitivity of model predictions to boundary conditions, cartilage geometry, and cartilage material properties. Loads based on in vivo data were applied to a cadaveric hip joint to simulate walking, descending stairs, and stair-climbing. Contact pressures and areas were measured using pressure sensitive film. CT image data were segmented and discretized into FE meshes of bone and cartilage. FE boundary and loading conditions mimicked the experimental testing. Fair to good qualitative correspondence was obtained between FE predictions and experimental measurements for simulated walking and descending stairs, while excellent agreement was obtained for stair-climbing. Experimental peak pressures, average pressures, and contact areas were 10.0 MPa (limit of film detection), 4.4-5.0 MPa, and 321.9-425.1 mm(2), respectively, while FE-predicted peak pressures, average pressures, and contact areas were 10.8-12.7 MPa, 5.1-6.2 MPa, and 304.2-366.1 mm(2), respectively. Misalignment errors, determined as the difference in root mean squared error before and after alignment of FE results, were less than 10%. Magnitude errors, determined as the residual error following alignment, were approximately 30% but decreased to 10-15% when the regions of highest pressure were compared. Alterations to the cartilage shear modulus, bulk modulus, or thickness resulted in +/-25% change in peak pressures, while changes in average pressures and contact areas were minor (+/-10%). When the pelvis and proximal femur were represented as rigid, there were large changes, but the effect depended on the particular loading scenario. Overall, the subject-specific FE predictions compared favorably with pressure film measurements and were in good agreement with published experimental data. The validated modeling framework provides a foundation for development of patient-specific FE models to investigate the mechanics of normal and pathological hips.

Radiographic prevalence of femoroacetabular impingement in collegiate football players: AAOS Exhibit Selection.

BACKGROUND The prevalence of femoroacetabular impingement may be greater in athletes than in the general population because of increased loading of the hip during sports. This study evaluated the radiographs of collegiate football players in order to quantify the prevalence of femoroacetabular impingement in asymptomatic athletes. METHODS Sixty-seven male collegiate football players (age, 21 ± 1.9 years) participated in this prospective study. Both hips (n = 134) were evaluated independently by two orthopaedic surgeons for radiographic signs of femoroacetabular impingement. The alpha angle and femoral head-neck offset were measured on frog-leg lateral radiographs. The lateral center-edge angle, acetabular index, crossover sign, and alpha angle were measured on anteroposterior radiographs. Data for continuous variables were averaged between observers prior to assessing prevalence. Cam femoroacetabular impingement was considered to be present if the femoral head-neck offset was <8 mm and/or the alpha angle was >50° on either radiograph. Pincer femoroacetabular impingement was considered to be present if the lateral center-edge angle was >40°, the acetabular index was <0°, and/or a positive crossover sign was detected by both observers. RESULTS Ninety-five percent of the 134 hips had at least one sign of cam or pincer impingement, and 77% had more than one sign. Twenty-one percent had only one sign of cam femoroacetabular impingement and 57% had both signs. Fifty-two percent had only one sign of pincer femoroacetabular impingement, 10% had two, and 4% had all three signs. Specifically, 72% had an abnormal alpha angle, 64% had a decreased femoral head-neck offset, 61% had a positive crossover sign, 16% had a decreased acetabular index, and 7% had an increased lateral center-edge angle. Fifty percent of all hips had at least one sign of pincer femoroacetabular impingement and at least one sign of cam impingement. Interobserver and intraobserver repeatability was moderate or better for each measure (range, 0.59 to 0.85). CONCLUSIONS Morphologic abnormalities associated with cam and pincer femoroacetabular impingement were common in these collegiate football players. The prevalence of cam and pincer femoroacetabular impingement was substantially higher than the previously reported prevalence in the general population.

Verification, validation and sensitivity studies in computational biomechanics

Computational techniques and software for the analysis of problems in mechanics have naturally moved from their origins in the traditional engineering disciplines to the study of cell, tissue and organ biomechanics. Increasingly complex models have been developed to describe and predict the mechanical behavior of such biological systems. While the availability of advanced computational tools has led to exciting research advances in the field, the utility of these models is often the subject of criticism due to inadequate model verification and validation (V&V). The objective of this review is to present the concepts of verification, validation and sensitivity studies with regard to the construction, analysis and interpretation of models in computational biomechanics. Specific examples from the field are discussed. It is hoped that this review will serve as a guide to the use of V&V principles in the field of computational biomechanics, thereby improving the peer acceptance of studies that use computational modeling techniques.

Effects of Idealized Joint Geometry on Finite Element Predictions of Cartilage Contact Stresses in the Hip

Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error approximately 0.5mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.

Role of the acetabular labrum in load support across the hip joint

The relatively high incidence of labral tears among patients presenting with hip pain suggests that the acetabular labrum is often subjected to injurious loading in vivo. However, it is unclear whether the labrum participates in load transfer across the joint during activities of daily living. This study examined the role of the acetabular labrum in load transfer for hips with normal acetabular geometry and acetabular dysplasia using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with and without the labrum during activities of daily living. The labrum in the dysplastic model supported 4-11% of the total load transferred across the joint, while the labrum in the normal model supported only 1-2% of the total load. Despite the increased load transferred to the acetabular cartilage in simulations without the labrum, there were minimal differences in cartilage contact stresses. This was because the load supported by the cartilage correlated with the cartilage contact area. A higher percentage of load was transferred to the labrum in the dysplastic model because the femoral head achieved equilibrium near the lateral edge of the acetabulum. The results of this study suggest that the labrum plays a larger role in load transfer and joint stability in hips with acetabular dysplasia than in hips with normal acetabular geometry.

Finite element prediction of cartilage contact stresses in normal human hips

Our objectives were to determine cartilage contact stress during walking, stair climbing, and descending stairs in a well‐defined group of normal volunteers and to assess variations in contact stress and area among subjects and across loading scenarios. Ten volunteers without history of hip pain or disease with normal lateral center‐edge angle and acetabular index were selected. Computed tomography imaging with contrast was performed on one hip. Bone and cartilage surfaces were segmented from volumetric image data, and subject‐specific finite element models were constructed and analyzed using a validated protocol. Acetabular contact stress and area were determined for seven activities. Peak stress ranged from 7.52 ± 2.11 MPa for heel‐strike during walking (233% BW) to 8.66 ± 3.01 MPa for heel‐strike during descending stairs (261% BW). Average contact area across all activities was 34% of the surface area of the acetabular cartilage. The distribution of contact stress was highly non‐uniform, and more variability occurred among subjects for a given activity than among activities for a single subject. The magnitude and area of contact stress were consistent between activities, although inter‐activity shifts in contact pattern were found as the direction of loading changed. Relatively small incongruencies between the femoral and acetabular cartilage had a large effect on the contact stresses. These effects tended to persist across all simulated activities. These results demonstrate the diversity and trends in cartilage contact stress in healthy hips during activities of daily living and provide a basis for future comparisons between normal and pathologic hips.

Validation of computational models in biomechanics

Abstract The topics of verification and validation have increasingly been discussed in the field of computational biomechanics, and many recent articles have applied these concepts in an attempt to build credibility for models of complex biological systems. Verification and validation are evolving techniques that, if used improperly, can lead to false conclusions about a system under study. In basic science, these erroneous conclusions may lead to failure of a subsequent hypothesis, but they can have more profound effects if the model is designed to predict patient outcomes. While several authors have reviewed verification and validation as they pertain to traditional solid and fluid mechanics, it is the intent of this paper to present them in the context of computational biomechanics. Specifically, the task of model validation will be discussed, with a focus on current techniques. It is hoped that this review will encourage investigators to engage and adopt the verification and validation process in an effort to increase peer acceptance of computational biomechanics models.

The effect of three-component total ankle replacement malalignment on clinical outcome: pain relief and functional outcome in 317 consecutive patients.

BACKGROUND Total ankle replacement has become an increasingly popular treatment for patients with end-stage ankle osteoarthritis. The surgery is technically demanding and generally performed by only experienced foot and ankle surgeons. An important complication of total ankle replacement is malposition of the talar component. The biomechanical effect of malposition has been reported; however, the functional outcomes of patients with varying degrees of talar component malposition have not. The purpose of this study was to assess the influence of talar component malposition on postoperative pain relief and functional outcome. METHODS This retrospective cohort study included 317 total ankle replacements in 317 patients. The anteroposterior offset ratio was measured with use of lateral ankle radiographs made with the patient in a standing, weight-bearing position. Patients were classified into one of three groups: those with an anteroposterior offset ratio of 0 (127 ankles), those with a ratio of >0 (103 ankles), and those with a ratio of <0 (eighty-seven ankles). Postoperative pain relief was assessed with use of a visual analogue scale. Functional outcome was assessed with the American Orthopaedic Foot & Ankle Society (AOFAS) hindfoot score and measurement of ankle range of motion. The mean duration of follow-up (and standard deviation) was 53.2 ± 18.4 months (range, twenty-four to ninety-eight months). RESULTS The postoperative pain level in the group with an anteroposterior offset ratio of 0 was significantly lower than that for both the group with a ratio of >0 (p < 0.001) and the group with a ratio of <0 (p = 0.017). Also, the functional outcome, measured with use of the AOFAS hindfoot score, was significantly higher, and ankle motion was significantly greater, in the group with an anteroposterior offset ratio of 0 than in the group with a ratio of >0 (p = 0.003 and p < 0.001, respectively) and the group with a ratio of <0 (p = 0.007 and p = 0.080). CONCLUSIONS The anteroposterior offset ratio may be a useful predictor of outcome in patients with total ankle replacement with regard to both pain and function. LEVEL OF EVIDENCE Prognostic Level II. See Instructions for Authors for a complete description of levels of evidence.

Patient-specific Analysis of Cartilage and Labrum Mechanics in Human Hips with Acetabular Dysplasia

BACKGROUND Acetabular dysplasia is a major predisposing factor for development of hip osteoarthritis (OA), and may result from alterations to chondrolabral loading. Subject-specific finite element (FE) modeling can be used to evaluate chondrolabral mechanics in the dysplastic hip, thereby providing insight into mechanics that precede OA. OBJECTIVE To evaluate chondrolabral contact mechanics and congruency in dysplastic hips and normal hips using a validated approach to subject-specific FE modeling. METHODS FE models of ten subjects with normal acetabula and ten subjects with dysplasia were constructed using a previously validated protocol. Labrum load support, and labrum and acetabular cartilage contact stress and contact area were compared between groups. Local congruency was determined at the articular surface for two simulated activities. RESULTS The labrum in dysplastic hips supported 2.8-4.0 times more of the load transferred across the joint than in normal hips. Dysplastic hips did not have significantly different congruency in the primary load-bearing regions than normal hips, but were less congruent in some unloaded regions. Normal hips had larger cartilage contact stress than dysplastic hips in the few regions that had significant differences. CONCLUSIONS The labrum in dysplastic hips has a far more significant role in hip mechanics than it does in normal hips. The dysplastic hip is neither less congruent than the normal hip, nor subjected to elevated cartilage contact stresses. This study supports the concept of an outside-in pathogenesis of OA in dysplastic hips and that the labrum in dysplastic hips should be preserved during surgery.