Three-dimensional dynamic hip contact area and pressure distribution during activities of daily living

Estimation of the hip joint contact area and pressure distribution during activities of daily living is important in predicting joint degeneration mechanism, prosthetic implant wear, providing biomechanical rationales for preoperative planning and postoperative rehabilitation. These biomechanical data were estimated utilizing a generic hip model, the Discrete Element Analysis technique, and the in vivo hip joint contact force data. The three-dimensional joint potential contact area was obtained from the anteroposterior radiograph of a subject and the actual joint contact area and pressure distribution in eight activities of daily living were calculated. During fast, normal, and slow walking, the peak pressure of moderate magnitude was located at the lateral roof of the acetabulum during mid-stance. In standing up and sitting down, and during knee bending, the peak pressures were located at the edge of the posterior horn and the magnitude of the peak pressure during sitting down was 2.8 times that of normal walking. The peak pressure was found at the lateral roof in climbing up stairs which was higher than that in going down stairs. These results can be used to rationalize rehabilitation protocols, functional restrictions after complex acetabular reconstructions, and prosthetic component wear and fatigue test set up. The same model and analysis can provide further insight to soft tissue loading and pathology such as labral injury. When the pressure distribution on the acetabulum is inverted onto the femoral head, prediction of subchondral bone collapse associated with avascular necrosis can be achieved with improved accuracy.     

A Three-Dimensional Finite Element Model for Biomechanical Analysis of the Hip

The objective of this study was to construct a three-dimensional (3D) finite element model of the hip. The images of the hip were obtained from Chinese visible human dataset. The hip model includes acetabular bone, cartilage, labrum, and bone. The cartilage of femoral head was constructed using the AutoCAD and Solidworks software. The hip model was imported into ABAQUS analysis system. The contact surface of the hip joint was meshed. To verify the model, the single leg peak force was loaded, and contact area of the cartilage and labrum of the hip and pressure distribution in these structures were observed. The constructed 3D hip model reflected the real hip anatomy. Further, this model reflected biomechanical behavior similar to previous studies. In conclusion, this 3D finite element hip model avoids the disadvantages of other construction methods, such as imprecision of cartilage construction and the absence of labrum. Further, it provides basic data critical for accurately modeling normal and abnormal loads, and the effects of abnormal loads on the hip.     

Biomechanics of the hip and birth in early Homo

A complex of traits in the femur and pelvis of Homo ereclus and early “erectus-like” specimens has been described, but never satisfactorily explained. Here the functional relationships between pelvic and femoral structure in humans are explored using both theoretical biomechanical models and empirical tests within modern samples of diverse body form (Pecos Amerindians, East Africans). Results indicate that a long femoral neck increases mediolateral bending of the femoral diaphysis and decreases gluteal abductor and hip joint reaction forces. Increasing biacetabular breadth along with femoral neck length further increases M-L bending of the femoral shaft and maintains abductor and joint reaction forces at near “normal” levels. When compared to modern humans, Homo erectus and early “erectus-like” specimens are characterized by a long femoral neck and greatly increased M-L relative to A-P bending strength of the femoral shaft, coupled with no decrease in hip joint size and a probable increase in abductor force relative to body size. All of this strongly suggests that biacetabular breadth as well as femoral neck length was relatively large in early Homo. Several features preserved in early Homo partial hip bones also indicate that the true (lower) pelvis was very M-L broad, as well as A-P narrow. This is similar to the lower pelvic shape of australopithecines and suggests that nonrotational birth, in which the newborn's head is oriented transversely through the pelvic outlet, characterized early Homo as well as Australopithecus. Because M-L breadth of the pelvis is constrained by other factors, this may have limited increases in cranial capacity within Homo until rotational birth was established during the late Middle Pleistocene. During or after the transition to rotational birth biacetabular breadth decreased, reducing the body weight moment arm about the hip and allowing femoral neck length (abductor moment arm) to also decrease, both of which reduced M-L bending of the proximal femoral shaft. Variation in femoral structural properties within early Homo and other East African Early Pleistocene specimens has several taxonomic and phylogenetic implications. © 1995 Wiley-Liss, Inc.     

How should dysplastic human hips be evaluated?

Dysplasia of the human hip is characterised by insufficient anterolateral covering of femoral head by the acetabulum. In our study, we evaluated dysplastic human hip joints using biomechanical parameters (the peak contact stress in the weight-bearing area of the hip - pmax) and X-ray image parameters (the centre-edge angle - thetaCE, the transverse acetabular inclination angle - thetaUS, the acetabular index of the weight-bearing zone - thetaAC, the ACM angle - thetaACM, and the hip value - HV). The purpose of this study is to make use of X-ray and biomechanical parameters to evaluate hips diagnosed with "hip dysplasia", and to establish whether or not there is a correlation between the two. Our results show a statistically significant correlation between pmax and thetaCE, thetaUS and thetaAC. The correlation between pmax and thetaACM and HV is not statistically significant.     

Association between contact hip stress and RSA-measured wear rates in total hip arthroplasties of 31 patients

BACKGROUND: The main concern in the long run of total hip replacements is aseptic loosening of the prosthesis. Optimization of the biomechanics of the hip joint is necessary for optimization of long-term success. A widely implementable tool to predict biomechanical consequences of preoperatively planned reconstructions still has to be developed. A potentially useful model to this purpose has been developed previously. The aim of this study is to quantify the association between the estimated hip joint contact force by this biomechanical model and RSA-measured wear rates in a clinical setting. METHODS: Thirty-one patients with a total hip replacement were measured with RSA, the gold standard for clinical wear measurements. The reference examination was done within 1 week of the operation and the follow-up examinations were done at 1, 2 and 5 years. Conventional pelvic X-rays were taken on the same day. The contact stress distribution in the hip joint was determined by the computer program HIPSTRESS. The procedure for the determination of the hip joint contact stress distribution is based on the mathematical model of the resultant hip force in the one-legged stance and the mathematical model of the contact stress distribution. The model for the force requires as input data, several geometrical parameters of the hip and the body weight, while the model for stress requires as input data, the magnitude and direction of the resultant hip force. The stress distribution is presented by the peak stress-the maximal value of stress on the weight-bearing area (p(max)) and also by the peak stress calculated with respect to the body weight (p(max)/W(B)) which gives the effect of hip geometry. Visualization of the relations between predicted values by the model and the wear at different points in the follow-up was done using scatterplots. Correlations were expressed as Pearson r values. RESULTS: The predicted p(max) and wear were clearly correlated in the first year post-operatively (r = 0.58, p = 0.002), while this correlation is weaker after 2 years (r = 0.19, p = 0.337) and 5 years (r = 0.24, p = 0.235). The wear values at 1, 2 and 5 years post-operatively correlate with each other in the way that is expected considering the wear velocity curve of the whole group. The correlation between the predicted p(max) values of two observers who were blinded for each other's results was very good (r = 0.93, p < 0.001). CONCLUSION: We conclude that the biomechanical model used in this paper provides a scientific foundation for the development of a new way of constructing preoperative biomechanical plans for total hip replacements.     

Subject-specific hip geometry affects predicted hip joint contact forces during gait

Hip loading affects bone remodeling and implant fixation. In this study, we have analyzed the effect of subject-specific modeling of hip geometry on muscle activation patterns and hip contact forces during gait, using musculoskeletal modeling, inverse dynamic analysis and static optimization. We first used sensitivity analysis to analyze the effect of isolated changes in femoral neck-length (NL) and neck-shaft angle (NSA) on calculated muscle activations and hip contact force during the stance phase of gait. A deformable generic musculoskeletal model was adjusted incrementally to adopt a physiological range of NL and NSA. In a second similar analysis, we adjusted hip geometry to the measurements from digitized radiographs of 20 subjects with primary hip osteoarthrosis. Finally, we studied the effect of hip abductor weakness on muscle activation patterns and hip contact force. This analysis showed that differences in NL (41-74 mm) and NSA (113-140 degrees ) affect the muscle activation of the hip abductors during stance phase and hence hip contact force by up to three times body weight. In conclusion, the results from both the sensitivity and subject-specific analysis showed that at the moment of peak contact force, altered NSA has only a minor effect on the loading configuration of the hip. Increased NL, however, results in an increase of the three hip contact-force components and a reduced vertical loading. The results of these analyses are essential to understand modified hip joint loading, and for planning hip surgery for patients with osteoarthrosis.     

Hip joint degeneration due to cam impingement: a finite element analysis

The goal of this study was to investigate the impact of cam impingement, a biomechanical risk factor, on hip joint degeneration and ultimately coxarthrosis. 3D finite element solid models of a healthy and a pathologic hip were developed based on clinical reports. The biphasic characteristics of cartilaginous tissues were considered to identify localised solid matrix overloading during normal walking and sitting down (SD). Localised femoral intrusion at the anterior-superior pelvic horn was revealed in the pathologic hip during SD, where the radial and meridional solid stresses in the acetabular cartilage and circumferential solid stresses within the acetabular labrum increased by 3.7, 1.5 and 2.7 times, respectively. The increased solid-on-solid stresses, reduction in fluid-load support and associated higher friction during articulation may result in joint wear and other degenerative changes in the hip.     

Hip contact stress during normal and staircase walking: the influence of acetabular anteversion angle and lateral coverage of the acetabulum

Hip contact stress is considered to be an important biomechanical factor related to development of coxarthrosis. The effect of the lateral coverage of the acetabulum on the hip contact stress has been demonstrated in several studies of hip dysplasia, whereas the effect of the anterior anteversion remains unclear. Therefore, the joint hip contact stress during normal level walking and staircase walking, in normal and dysplastic hips, for small and large acetabular anteversion angle was computed. For small acetabular anteversion angle, the hip contact stress is slightly increased (less than 15%) in staircase walking when compared with normal walking. In hips with large angle of acetabular anteversion, walking downstairs significantly increases the maximal peak contact stress (70% in normal hips and 115% in dysplastic hips) whereas walking upstairs decreases the peak contact stress (4% in normal hips and 34% in dysplastic hips) in comparison to normal walking. Based on the presented results, we suggest that the acetabular anteversion should be considered in biomechanical evaluation of the hips, especially when the lateral coverage of the acetabulum is small.     

"In vivo" determination of hip joint separation and the forces generated due to impact loading conditions

Numerous supporting structures assist in the retention of the femoral head within the acetabulum of the normal hip joint including the capsule, labrum, and ligament of the femoral head (LHF). During total hip arthroplasty (THA), the LHF is often disrupted or degenerative and is surgically removed. In addition, a portion of the remaining supporting structures is transected or resected to facilitate surgical exposure. The present study analyzes the effects of LHF absence and surgical dissection in THA patients. Twenty subjects (5 normal hip joints, 10 nonconstrained THA, and 5 constrained THA) were evaluated using fluoroscopy while performing active hip abduction. All THA subjects were considered clinically successful. Fluoroscopic videos of the normal hips were analyzed using digitization, while those with THA were assessed using a computerized interactive model-fitting technique. The distance between the femoral head and acetabulum was measured to determine if femoral head separation occurred. Error analysis revealed measurements to be accurate within 0.75mm. No separation was observed in normal hips or those subjects implanted with constrained THA, while all 10 (100%) with unconstrained THA demonstrated femoral head separation, averaging 3.3mm (range 1.9-5.2mm). This study has shown that separation of the prosthetic femoral head from the acetabular component can occur. The normal hip joint has surrounding capsuloligamentous structures and a ligament attaching the femoral head to the acetabulum. We hypothesize that these soft tissue supports create a passive, resistant force at the hip, preventing femoral head separation. The absence of these supporting structures after THA may allow increased hip joint forces, which may play a role in premature polyethylene wear or prosthetic loosening.     

A Wider Pelvis Does Not Increase Locomotor Cost in Humans,with Implications for the Evolution of Childbirth

The shape of the human female pelvis is thought to reflect an evolutionary trade-off between two competing demands: a pelvis wide enough to permit the birth of large-brained infants, and narrow enough for efficient bipedal locomotion. This trade-off, known as the obstetrical dilemma, is invoked to explain the relative difficulty of human childbirth and differences in locomotor performance between men and women. The basis for the obstetrical dilemma is a standard static biomechanical model that predicts wider pelves in females increase the metabolic cost of locomotion by decreasing the effective mechanical advantage of the hip abductor muscles for pelvic stabilization during the single-leg support phase of walking and running, requiring these muscles to produce more force. Here we experimentally test this model against a more accurate dynamic model of hip abductor mechanics in men and women. The results show that pelvic width does not predict hip abductor mechanics or locomotor cost in either women or men, and that women and men are equally efficient at both walking and running. Since a wider birth canal does not increase a woman’s locomotor cost, and because selection for successful birthing must be strong, other factors affecting maternal pelvic and fetal size should be investigated in order to help explain the prevalence of birth complications caused by a neonate too large to fit through the birth canal.     

Frontal joint dynamics when initiating stair ascent from a walk versus a stand

Ascending stairs is a challenging activity of daily living for many populations. Frontal plane joint dynamics are critical to understand the mechanisms involved in stair ascension as they contribute to both propulsion and medio-lateral stability. However, previous research is limited to understanding these dynamics while initiating stair ascent from a stand. We investigated if initiating stair ascent from a walk with a comfortable self-selected speed could affect the frontal plane lower-extremity joint moments and powers as compared to initiating stair ascent from a stand and if this difference would exist at consecutive ipsilateral steps on the stairs. Kinematics data using a 3-D motion capture system and kinetics data using two force platforms on the first and third stair treads were recorded simultaneously as ten healthy young adults ascended a custom-built staircase. Data were collected from two starting conditions of stair ascent, from a walk (speed: 1.42 ± 0.21 m/s) and from a stand. Results showed that subjects generated greater peak knee abductor moment and greater peak hip abductor moment when initiating stair ascent from a walk. Greater peak joint moments and powers at all joints were also seen while ascending the second ipsilateral step. Particularly, greater peak hip abductor moment was needed to avoid contact of the contralateral limb with the intermediate step by counteracting the pelvic drop on the contralateral side. This could be important for therapists using stair climbing as a testing/training tool to evaluate hip strength in individuals with documented frontal plane abnormalities (i.e. knee and hip osteoarthritis, ACL injury).     

The influence of the biomechanical parameters of the hip on the outcome of treatment of hips subject to avascular necrosis of the femoral head

We were interested in whether or not the biomechanical status of the hip influences the course of avascular necrosis of the femoral head. To investigate this, we used a computer aided system based on a three dimensional mathematical model for determining the stress distribution in the hip joint from standard anteroposterior rentgenographs (X-ray images) of both hips and pelvis. Based on the results of our study, we suggest that the biomechanical parameters of the hip play an important role in the outcome of treatment of hips affected by avascular necrosis of the femoral head.     

Influence of weak hip abductor muscles on joint contact forces during normal walking: probabilistic modeling analysis

The weakness of hip abductor muscles is related to lower-limb joint osteoarthritis, and joint overloading may increase the risk for disease progression. The relationship between muscle strength, structural joint deterioration and joint loading makes the latter an important parameter in the study of onset and follow-up of the disease. Since the relationship between hip abductor weakness and joint loading still remains an open question, the purpose of this study was to adopt a probabilistic modeling approach to give insights into how the weakness of hip abductor muscles, in the extent to which normal gait could be unaltered, affects ipsilateral joint contact forces. A generic musculoskeletal model was scaled to each healthy subject included in the study, and the maximum force-generating capacity of each hip abductor muscle in the model was perturbed to evaluate how all physiologically possible configurations of hip abductor weakness affected the joint contact forces during walking. In general, the muscular system was able to compensate for abductor weakness. The reduced force-generating capacity of the abductor muscles affected joint contact forces to a mild extent, with 50th percentile mean differences up to 0.5 BW (maximum 1.7 BW). There were greater increases in the peak knee joint loads than in loads at the hip or ankle. Gluteus medius, particularly the anterior compartment, was the abductor muscle with the most influence on hip and knee loads. Further studies should assess if these increases in joint loading may affect initiation and progression of osteoarthritis.     

Computer determination of contact stress distribution and size of weight bearing area in the human hip joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

A 42-year-old patient presenting with femoral head migration after hemiarthroplasty performed 22 years earlier: a case report

Introduction

Treatment of femoral neck fractures in young adults may require total hip arthroplasty or hip hemiarthroplasty using a bipolar cup. The latter can, however, result in migration of the femoral head and poor long-term results.

Case presentation

We report a case of femoral head migration after hemiarthroplasty performed for femoral neck fracture that had occurred 22 years earlier, when the patient (a Japanese man) was 20 years old. He experienced peri-prosthetic fracture of the femur, subsequent migration of the prosthesis, and a massive bone defect of the pelvic side acetabular roof. After bone union of the femoral shaft fracture, the patient was referred to our hospital for reconstruction of the acetabular roof. Intra-operatively, we placed two alloimplants of bone from around the transplanted femoral head into the weight-bearing region of the acetabular roof using an impaction bone graft method. We then implanted an acetabular roof reinforcement plate and a cemented polyethylene cup in the position of the original acetabular cup. Eighteen months post-operatively, X-rays showed union of the transplanted bone.

Conclusions

Treatment of femoral neck fractures in young adults is usually accomplished by osteosynthesis, but it may be complicated by femoral head avascular necrosis or by infection or osteomyelitis. In such cases, once an infection has subsided, either hip hemiarthroplasty using a bipolar cup or total hip arthroplasty may be required. However, if the acetabular side articular cartilage is damaged, a bipolar cup should not be used. Total hip arthroplasty should be performed to prevent migration of the implant.     

Computer Determination of Contact Stress Distribution and Size of Weight Bearing Area in the Human Hip Joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

The influence of the acetabular labrum on hip joint cartilage consolidation: a poroelastic finite element model

The goal of this study was to investigate the influence of the acetabular labrum on the consolidation, and hence the solid matrix strains and stresses, of the cartilage layers of the hip joint. A plane-strain finite element model was developed, which represented a coronal slice through the acetabular and femoral cartilage layers and the acetabular labrum. Elements with poroelastic properties were used to account for the biphasic solid/fluid nature of the cartilage and labrum. The response of the joint over an extended period of loading (10,000s) was examined to simulate the nominal compressive load that the joint is subjected to throughout the day. The model demonstrated that the labrum adds an important resistance in the flow path of the fluid being expressed from the cartilage layers of the joint. Cartilage layer consolidation was up to 40% quicker in the absence of the labrum. Following removal of the labrum from the model, the solid-on-solid contact stresses between the femoral and acetabular cartilage layers were greatly increased (up to 92% higher), which would increase the friction between the joint surfaces. In the absence of the labrum, the centre of contact shifted towards the acetabular rim. Subsurface strains and stresses were much higher without the labrum, which could contribute to fatigue damage of the cartilage layers. Finally, the labrum provided some structural resistance to lateral motion of the femoral head within the acetabulum, enhancing joint stability and preserving joint congruity.     

Evaluation of Constant Thickness Cartilage Models vs. Patient Specific Cartilage Models for an Optimized Computer-Assisted Planning of Periacetabular Osteotomy

Modern computerized planning tools for periacetabular osteotomy (PAO) use either morphology-based or biomechanics-based methods. The latter relies on estimation of peak contact pressures and contact areas using either patient specific or constant thickness cartilage models. We performed a finite element analysis investigating the optimal reorientation of the acetabulum in PAO surgery based on simulated joint contact pressures and contact areas using patient specific cartilage model. Furthermore we investigated the influences of using patient specific cartilage model or constant thickness cartilage model on the biomechanical simulation results. Ten specimens with hip dysplasia were used in this study. Image data were available from CT arthrography studies. Bone models were reconstructed. Mesh models for the patient specific cartilage were defined and subsequently loaded under previously reported boundary and loading conditions. Peak contact pressures and contact areas were estimated in the original position. Afterwards we used a validated preoperative planning software to change the acetabular inclination by an increment of 5° and measured the lateral center edge angle (LCE) at each reorientation position. The position with the largest contact area and the lowest peak contact pressure was defined as the optimal position. In order to investigate the influence of using patient specific cartilage model or constant thickness cartilage model on the biomechanical simulation results, the same procedure was repeated with the same bone models but with a cartilage mesh of constant thickness. Comparison of the peak contact pressures and the contact areas between these two different cartilage models showed that good correlation between these two cartilage models for peak contact pressures (r = 0.634 ∈ [0.6, 0.8], p < 0.001) and contact areas (r = 0.872 > 0.8, p < 0.001). For both cartilage models, the largest contact areas and the lowest peak pressures were found at the same position. Our study is the first study comparing peak contact pressures and contact areas between patient specific and constant thickness cartilage models during PAO planning. Good correlation for these two models was detected. Computer assisted planning with FE modeling using constant thickness cartilage models might be a promising PAO planning tool when a conventional CT is available.     

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%.     

Hemiarthroplasty of hip joint: An experimental validation using porcine acetabulum

Biphasic properties of articular cartilage allow it to be an excellent bearing material and have been studied through several simplified experiments as well as finite element modelling. However, three-dimensional biphasic finite element (FE) models of the whole joint are rare. The current study was carried out to experimentally validate FE methodology for modelling hemiarthroplasty. Material properties such as equilibrium elastic modulus and permeability of porcine acetabular cartilage were initially derived by curve-fitting an experimental deformation curve with that obtained using FE. These properties were then used in the hemiarthroplasty hip joint modelling. Each porcine acetabular cup was loaded with 400N using a 34mm diameter CoCr femoral head. A specimen-specific FE model of each acetabular cup was created using μCT and a series of software processes. Each model was analysed under conditions similar to those tested experimentally. Contact stresses and contact areas predicted by the model, immediately after loading, were then compared with the corresponding experimentally measured values. Very high peak contact stresses (maximum experimental: 14.09MPa) were recorded. A maximum difference of 12.42% was found in peak contact stresses. The corresponding error for contact area was 20.69%. Due to a fairly good agreement in predicted and measured values of contact stresses and contact areas, the integrated methodology developed in this study can be used as a basis for future work. In addition, FE predicted total fluid load support was around 80% immediately after loading. This was lower than that observed in conforming contact problems involving biphasic cartilage and was due to a smaller local contact area and variable clearance making fluid exudation easier.