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.     

Secondary hip dysplasia increases risk for early coxarthritis after Legg-Calve-Perthes disease. A study of 255 hips

Abstract

The biomechanical parameters of the hip joint articular surface were analysed in 141 adult hips after Legg-Calve Perthes Disease, and 114 contralateral unaffected hips (controls), by using HIPSTRESS mathematical models. Geometrical parameters, assessed from anteroposterior and axial radiograms, were used as input to models for resultant hip force and contact hip stress. Results confirm previous indications that head enlargement after the Legg-Calve-Perthes Disease compensates the values of hip stress. Furthermore, it was found that an increased risk for coxarthritis development after the disease is secondary to concomitant hip dysplasia, with considerable and statistically significantly lower centre-edge angle and unfavourable distribution of stress.     

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.     

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.     

Joint contact stresses calculated for acetabular dysplasia patients using discrete element analysis are significantly influenced by the applied gait pattern

Gait modifications in acetabular dysplasia patients may influence cartilage contact stress patterns within the hip joint, with serious implications for clinical outcomes and the risk of developing osteoarthritis. The objective of this study was to understand how the gait pattern used to load computational models of dysplastic hips influences computed joint mechanics. Three-dimensional pre- and post-operative hip models of thirty patients previously treated for hip dysplasia with periacetabular osteotomy (PAO) were developed for performing discrete element analysis (DEA). Using DEA, contact stress patterns were calculated for each pre- and post-operative hip model when loaded with an instrumented total hip, a dysplastic, a matched control, and a normal gait pattern. DEA models loaded with the dysplastic and matched control gait patterns had significantly higher (p = 0.012 and p < 0.001) average pre-operative maximum contact stress than models loaded with the normal gait. Models loaded with the dysplastic and matched control gait patterns had nearly significantly higher (p = 0.051) and significantly higher (p = 0.008) average pre-operative contact stress, respectively, than models loaded with the instrumented hip gait. Following PAO, the average maximum contact stress for DEA models loaded with the dysplastic and matched control patterns decreased, which was significantly different (p < 0.001) from observed increases in maximum contact stress calculated when utilizing the instrumented hip and normal gait patterns. The correlation between change in DEA-computed maximum contact stress and the change in radiographic measurements of lateral center-edge angle were greatest (R² = 0.330) when utilizing the dysplastic gait pattern. These results indicate that utilizing a dysplastic gait pattern to load DEA models may be a crucial element to capturing contact stress patterns most representative of this patient population.     

Biomechanical and clinical alterations of the hip joint following femoral neck fracture and implantation of bipolar hip endoprosthesis

The implantation of a bipolar partial hip endoprosthesis is a treatment of choice for displaced medial femoral neck fracture. We present an experimental study which asses and compare biomechanical and clinical status through period before and after hip fracture and implantation of bipolar partial hip endoprosthesis. This study encompassed 75 patients who suffered from an acute medial femoral neck fracture and were treated with the implantation of a bipolar partial hip endoprosthesis. Their biomechanical status (stress distribution on the hip joint weight bearing area) and clinical status (Harris Hip Score) were estimated for the time prior to the injury and assessed at the follow-up examination that was, on average, carried out 40 months after the operation. Despite ageing, the observed Harris Hip Score at the follow-up examination was higher than that estimated prior to the injury (77.9 > 69.6; p = 0.006). Similarly, the hip stress distribution was reduced (2.7 MPa < 2.3 MPa; p = 0.001). While this reduction can be attributed to a loss of weight due to late ageing, the principal improvement came from the operative treatment and corresponding restoration of the biomechanical properties of the hip joint. The implantation of a bipolar partial hip endoprosthesis for patients with displaced medial femoral neck fractures improves the biomechanical and clinical features of the hip, what should have on mind during making decision about treatment.     

Normal hip joint contact pressure distribution in single-leg standing--effect of gender and anatomic parameters.

A practical and easy-to-use analysis technique that can study the patient's hip joint contact force/pressure distribution would be useful to assess the effect of abnormal biomechanical conditions and anatomical deformities on joint contact stress for treatment planning purpose. This technique can also help to establish the normative database on hip joint contact pressure distribution in men and women in different age groups. Twelve anatomic parameters and seven biomechanical parameters of the hip joint in a normal population (41 females, 15 males) were calculated. The inter-parameter correlations were investigated. The pressure distribution in the hip joint was calculated using a three-dimensional discrete element analysis (DEA) technique. The 3D contact geometry of the hip joint was estimated from a 2D radiograph by assuming that the femoral head and the acetabular surface were spherical in shape. The head-trochanter ratio (HT), femoral head radius, pelvic height, the joint contact area, the normalized peak contact pressure, abductor force, and the joint contact force were significantly different between men and women. The normalized peak contact pressure was correlated both with acetabular coverage and head-trochanter ratio. Change of abductor force direction within normal variation did not affect the joint peak contact pressure. However, in simulated dysplastic conditions when the CE angle is small or negative, abductor muscle direction becomes very sensitive in joint contact pressure estimation. The models and the results presented can be used as the reference base in computer simulation for preoperative planning in pelvic or femoral osteotomy.     

Effect of hip joint angle at seat-off on hip joint contact force during sit-to-stand movement: a computer simulation study

Background

Sit-to-stand movements are a necessary part of daily life, and excessive mechanical stress on the articular cartilage has been reported to encourage the progression of osteoarthritis. Although a change in hip joint angle at seat-off may affect hip joint contact force during a sit-to-stand movement, the effect is unclear. This study aimed to examine the effect of the hip joint angle at seat-off on the hip joint contact force during a sit-to-stand movement by using a computer simulation.

Methods

A musculoskeletal model was created for the computer simulation, and eight muscles were attached to each lower limb. Various sit-to-stand movements were generated using parameters (e.g., seat height and time from seat-off to standing posture) reported by previous studies. The hip joint contact force for each sit-to-stand movement was calculated. Furthermore, the effect of the hip joint angle at seat-off on the hip joint contact force during the sit-to-stand movement was examined. In this study, as the changes to the musculoskeletal model parameters affect the hip joint contact force, a sensitivity analysis was conducted.

Results and conclusions

The hip joint contact force during the sit-to-stand movement increased approximately linearly as the hip flexion angle at the seat-off increased. Moreover, the normal sit-to-stand movement and the sit-to-stand movement yielding a minimum hip joint contact force were approximately equivalent. The effect of the changes to the musculoskeletal model parameters on the main findings of this study was minimal. Thus, the main findings are robust and may help prevent the progression of hip osteoarthritis by decreasing mechanical stress, which will be explored in future studies.

     

Effect of cup inclination on predicted contact stress-induced volumetric wear in total hip replacement

In order to increase the lifetime of the total hip endoprosthesis, it is necessary to understand mechanisms leading to its failure. In this work, we address volumetric wear of the artificial cup, in particular the effect of its inclination with respect to the vertical. Volumetric wear was calculated by using mathematical models for resultant hip force, contact stress and penetration of the prosthesis head into the cup. Relevance of the dependence of volumetric wear on inclination of the cup (its abduction angle ?_A) was assessed by the results of 95 hips with implanted endoprosthesis. Geometrical parameters obtained from standard antero-posterior radiographs were taken as input data. Volumetric wear decreases with increasing cup abduction angle ?_A. The correlation within the population of 95 hips was statistically significant (P = 0.006). Large cup abduction angle minimises predicted volumetric wear but may increase the risk for dislocation of the artificial head from the cup in the one-legged stance. Cup abduction angle and direction of the resultant hip force may compensate each other to achieve optimal position of the cup with respect to wear and dislocation in the one-legged stance for a particular patient.     

基于三维CT的模拟手术在大龄DDH患儿外科治疗中的应用研究

目的:通过对大龄发育性髋关节脱位(developmental dysplasia of the hip,DDH)患儿进行术前模拟手术,实现术中精确截骨及旋转角度,从而达到个体化治疗,改善患儿预后的目的。方法:本研究按照术前规划方式分为两组,一组为传统手术组;另一组为模拟手术组。共20例患儿均采用骨盆三联截骨术+股骨截骨术治疗,传统手术组10例,模拟手术组10例,手术时平均年龄为11.3岁,平均随访时间24.2个月。所有患儿均于术前行骨盆三维重建CT检查,测量CE角、股骨前倾角及髋臼指数,在mimics软件中,模拟手术方案,确定术中股骨截骨需要旋转的角度及骨盆截骨的位置,术中按照模拟手术的结果进行操作。术前评价指标使用Tonnis分级,术后评价指标使用改进的Trevor评分系统。结果:模拟手术组Tonnis分级3级4髋,Tonnis分级4级8髋;传统手术组Tonnis分级3级4髋,Tonnis分级4级9髋,两组患儿术前严重程度无显著性差异。依据Trevor评分,模拟手术组8髋(67%)优秀,3髋(25%)良好,1髋(8%)一般。传统手术组5髋(38%)优秀,5髋(38%)良好,3髋(23%)一般。两组有显著性差异。并发症:术后传统手术组3例患儿有不同程度的股骨头坏死。结论:大龄DDH患儿术前模拟手术,可以达到术中精确截骨及旋转角度,可改善患儿预后,实现该类患者的个体化治疗。     

Discrete element analysis is a valid method for computing joint contact stress in the hip before and after acetabular fracture

Evaluation of abnormalities in joint contact stress that develop after inaccurate reduction of an acetabular fracture may provide a potential means for predicting the risk of developing post-traumatic osteoarthritis. Discrete element analysis (DEA) is a computational technique for calculating intra-articular contact stress distributions in a fraction of the time required to obtain the same information using the more commonly employed finite element analysis technique. The goal of this work was to validate the accuracy of DEA-computed contact stress against physical measurements of contact stress made in cadaveric hips using Tekscan sensors. Four static loading tests in a variety of poses from heel-strike to toe-off were performed in two different cadaveric hip specimens with the acetabulum intact and again with an intentionally malreduced posterior wall acetabular fracture. DEA-computed contact stress was compared on a point-by-point basis to stress measured from the physical experiments. There was good agreement between computed and measured contact stress over the entire contact area (correlation coefficients ranged from 0.88 to 0.99). DEA-computed peak contact stress was within an average of 0.5 MPa (range 0.2–0.8 MPa) of the Tekscan peak stress for intact hips, and within an average of 0.6 MPa (range 0–1.6 MPa) for fractured cases. DEA-computed contact areas were within an average of 33% of the Tekscan-measured areas (range: 1.4–60%). These results indicate that the DEA methodology is a valid method for accurately estimating contact stress in both intact and fractured hips.     

内固定与关节置换手术治疗骨质疏松性髋部骨折的临床效果及术后并发症的影响因素分析

目的:分析内固定与关节置换手术治疗骨质疏松性髋部骨折的临床效果及其术后并发症的影响因素。方法:将2017年4月至2018年5月因骨质疏松性髋部骨折于我院进行手术治疗的78例患者作为研究对象,参考患者自身意愿按照采取手术方案的不同将所有患者分为内固定组与关节置换组。内固定组主要采用动力髋螺钉内固定治疗;关节置换组采用全髋关节置换术治疗。对比分析两组治疗后并发症的发生情况及关节功能恢复程度。结果:关节置换组手术治疗后关节功能的恢复情况显著优于内固定组(P0.05),术后并发症发生率低于内固定组(P0.05);患者术后并发症发生的主要影响因素包括:术前存在合并症、手术时机≥2h、采用内固定手术。结论:与内固定术相比,关节置换手术治疗骨质疏松性髋部骨折患者的疗效和安全性均更高,但采用该手术治疗时需注意手术操作的规范性及手术时机。     

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 ～0.5 mm). 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.     

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.     

Changes in Parkinsonian gait kinematics with self-generated and externally-generated cues: a comparison of responders and non-responders

Abstract

Purpose: Rhythmic auditory stimulation such as listening to music can alleviate gait bradykinesia in people with Parkinson disease (PD) by increasing spatiotemporal gait features. However, evidence about what specific kinematic alterations lead to these improvements is limited, and differences in responsiveness to cueing likely affect individual motor strategies. Self-generated cueing techniques, such as singing or mental singing, provide similar benefits but no evidence exists about how these techniques affect lower limb joint movement. In this study, we assessed immediate effects of external and self-generated cueing on lower limb movement trajectories during gait.

Methods: Using 3D motion capture, we assessed sagittal plane joint angles at the hip, knee, and ankle across 35 participants with PD, divided into responders (n?=?23) and non-responders (n?=?12) based on a clinically meaningful change in gait speed. Joint motion was assessed as overall range of motion as well as at two key time points during the gait cycle: initial contact and toe-off.

Results: Responders used both cue types to increase gait speed and induce increases in overall joint ROM at the hip while only self-generated cues also increased ROM at the ankle. Increased joint excursions for responders were also evident at initial contact and toe-off.

Conclusions: Our results indicate that self-generated rhythmic cues can induce similar increases in joint excursions as externally-generated cues and that some people may respond more positively than others. These results provide important insight into how self-generated cueing techniques may be tailored to meet the varied individual needs of people with PD.     

A Voxel-based Formulation for Contact Finite Element Analysis

To date, voxel-based finite element models have not been feasible for contact problems, owing to the inherent stair-step boundary discontinuities. New preprocessing techniques are reported herein to mesh these boundaries smoothly, for purposes of contact stress analysis. Further, new techniques are reported to concentrate the mesh resolution automatically near the articular surface, thus reducing the problem size to levels compatible with executing nonlinear problems on contemporary engineering workstations. Close approximations to Hertzian analytical solutions were obtained for spherical and cylindrical geometries meshed in this manner, and an illustrative anatomical contact problem of the human hip joint is presented.     

Subject-specific geometry affects acetabular contact pressure during gait more than subject-specific loading patterns

Abstract

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics.     

Effect of knee position on hip and knee torques during the barbell squat

Some recommendations suggest keeping the shank as vertical as possible during the barbell squat, thus keeping the knees from moving past the toes. This study examined joint kinetics occurring when forward displacement of the knees is restricted vs. when such movement is not restricted. Seven weight-trained men (mean +/- SD; age = 27.9 +/- 5.2 years) were videotaped while performing 2 variations of parallel barbell squats (barbell load = body weight). Either the knees were permitted to move anteriorly past the toes (unrestricted) or a wooden barrier prevented the knees from moving anteriorly past the toes (restricted). Differences resulted between static knee and hip torques for both types of squat as well as when both squat variations were compared with each other (p < 0.05). For the unrestricted squat, knee torque (N.m; mean +/- SD) = 150.1 +/- 50.8 and hip torque = 28.2 +/- 65.0. For the restricted squat, knee torque = 117.3 +/- 34.2 and hip torque = 302.7 +/- 71.2. Restricted squats also produced more anterior lean of the trunk and shank and a greater internal angle at the knees and ankles. The squat technique used can affect the distribution of forces between the knees and hips and on the kinematic properties of the exercise. PRACTICAL APPLICATIONS: Although restricting forward movement of the knees may minimize stress on the knees, it is likely that forces are inappropriately transferred to the hips and low-back region. Thus, appropriate joint loading during this exercise may require the knees to move slightly past the toes.     

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.     

Assessment of the applicability of the Hertzian contact theory to edge-loaded prosthetic hip bearings

The components of prosthetic hip bearings may experience in-vivo subluxation and edge loading on the acetabular socket as a result of joint laxity, causing abnormally high, damaging contact stresses. In this research, edge-loaded contact of prosthetic hips is examined analytically and experimentally in the most commonly used categories of material pairs. In edge-loaded ceramic-on-ceramic hips, the Hertzian contact theory yields accurate (conservatively, <10% error) predictions of the contact dimensions. Moreover, the Hertzian theory successfully captures slope and curvature trends in the dependence of contact patch geometry on the applied load. In an edge-loaded ceramic-on-metal pair, a similar degree of accuracy is observed in the contact patch length; however, the contact width is less accurately predicted due to the onset of subsurface plasticity, which is predicted for loads >400N. The Hertzian contact theory is shown to be ill-suited to edge-loaded ceramic-on-polyethylene pairs due to polyethylene's nonlinear material behavior. This work elucidates the methods and the accuracy of applying classical contact theory to edge-loaded hip bearings. The results help to define the applicability of the Hertzian theory to the design of new components and materials to better resist severe edge loading contact stresses.