Force track analysis of contemporary hip simulators

In an earlier paper, the authors presented the first verified method of computation of slide tracks in the relative motion between femoral head and acetabular cup of total hip prostheses. The method was applied for gait and for two hip simulator designs, and in a subsequent paper, for another eight designs. In the present paper, the track drawn by the resultant contact force, the so-called force track, was studied in depth. The variations of sliding distance, sliding velocity and direction of sliding during a cycle, all of which are important with respect to wear, were computed for gait and for 11 hip simulator designs. Moreover, the product of the instantaneous load and increment of sliding distance was numerically integrated over a cycle. This integral makes it possible to compare clinical wear rates with those produced by hip simulators in terms of a wear factor. For the majority of contemporary hip simulators, the integral has so far been unknown. The computations revealed considerable differences, which are likely to explain the substantial differences in wear produced by the simulators. With the most common head diameter, 28 mm, the ranges for sliding distance per cycle, mean sliding velocity, total change of direction of sliding and integral were: 19.7-34.3 mm, 19.7-49.0 mm/s, 360-1513 degrees, and 17.4-43.5 Nm, respectively.     

Enhanced computational prediction of polyethylene wear in hip joints by incorporating cross-shear and contact pressure in additional to load and sliding distance: Effect of head diameter

A new definition of the experimental wear factor was established and reported as a function of cross-shear motion and contact pressure using a multi-directional pin-on-plate wear testing machine for conventional polyethylene in the present study. An independent computational wear model was developed by incorporating the cross-shear motion and contact pressure-dependent wear factor into the Archard's law, in additional to load and sliding distance. The computational prediction of wear volume was directly compared with a simulator testing of a polyethylene hip joint with a 28 mm diameter. The effect of increasing the femoral head size was subsequently considered and was shown to increase wear, as a result of increased sliding distance and reduced contact pressure.     

Comparison of global and joint-to-joint methods for estimating the hip joint load and the muscle forces during walking

A three-dimensional musculoskeletal model of the lower limb was developed to study the influence of biarticular muscles on the muscle force distribution and joint loads during walking. A complete walking cycle was recorded for 9 healthy subjects using the standard optoelectronic motion tracking system. Ground contact forces were also measured using a 6-axes force plate. Inverse dynamics was used to compute net joint reactions (forces and torques) in the lower limb. A static optimization method was then used to estimate muscle forces. Two different approaches were used: in the first one named global method, the biarticular muscles exerted a torque on the two joints they spanned at the same time, and in the second one called joint-by-joint method, these biarticular muscles were divided into two mono-articular muscles with geometrical (insertion, origin, via points) and physiological properties remained unchanged. The hip joint load during the gait cycle was then calculated taking into account the effect of muscle contractions. The two approaches resulted in different muscle force repartition: the biarticular muscles were favoured over any set of single-joint muscles with the same physiological function when using the global method. While the two approaches yielded only little difference in the resultant hip load, the examination of muscle power showed that biarticular muscles could produce positive work at one joint and negative work at the other, transferring energy between body segments and thus decreasing the metabolic cost of movement.     

Contributions of individual muscles to hip joint contact force in normal walking

The human hip joint withstands high contact forces during daily activity and is therefore susceptible to injury and structural deterioration over time. Knowledge of muscle-force contributions to hip joint loading may assist in the development of strategies to prevent and manage conditions such as osteoarthritis, femoro-acetabular impingement and fracture. The main aim of this study was to determine the contributions of individual muscles to hip contact force in normal walking. Muscle contributions to hip contact force were calculated based on a previously published dynamic optimization solution for normal walking, which provided the time histories of joint motion, ground reaction forces, and muscle forces during the stance and swing phases of gait. The force developed by each muscle plus its contribution to the ground reaction force were used to determine the muscle’s contribution to hip contact force. Muscles were the major contributors to hip contact force, with gravitational and centrifugal forces combined contributing less than 5% of the total contact force. Four muscles that span the hip – gluteus medius, gluteus maximus, iliopsoas, and hamstrings – contributed most significantly to the three components of the hip contact force and hip contact impulse (integral of hip contact force over time). Three muscles that do not span the hip – vasti, soleus, and gastrocnemius – also contributed substantially to hip joint loading. These results provide additional insight into lower-limb muscle function during walking and may also be relevant to studies of cartilage degeneration and bone remodelling at the hip.     

Slide track analysis of eight contemporary hip simulator designs

In an earlier paper, the authors presented a new method of computation of slide tracks in the relative motion between femoral head and acetabular cup of total hip prostheses. For the first time, computed tracks were verified experimentally and with an alternative method of computation. Besides being an efficient way to illustrate hip kinematics, the shapes of the slide tracks are known to be of fundamental importance regarding the wear behaviour of prostheses. The verified method was now applied to eight contemporary hip simulator designs. The use of correct motion waveforms and an Euler sequence of rotations in each case was again found to be essential. Considerable differences were found between the simulators. For instance, the shapes of the tracks drawn by the resultant contact force included a circle, ellipse, irregular oval, leaf, twig, and straight line. Computation of tracks correctly for the most widely used hip simulator, known as biaxial, was made possible by the insight that the device is actually three-axial. Slide track patterns have now been computed for virtually all contemporary hip simulators, and both for the heads and for the cups. This comparative analysis forms a valuable basis for studies on the relationship between the type of multidirectional motion and wear. These studies can produce useful information for the design of joint simulators, and improve the understanding of wear phenomena in prosthetic joints.     

A dynamic model of the forearm including fatigue

A biomechanical model of the forearm, consisting of 61 muscle-tendon systems or tendons and 8 sections, is presented. The model can be used to calculate the muscle forces when resultant of the external forces and the motion is known. Calculations are based on constraints of muscle forces, joint forces, contact forces, and tendon junctions, and a load sharing principle telling which of the feasible solutions are likely and which are not. Fatigue is accounted for by updating the upper limits of the muscle forces according to the loading history. As an example, the model is used to predict the load sharing between the fingers when they are pressed against a table with a given total force.     

Tibiofemoral joint contact force in deep knee flexion and its consideration in knee osteoarthritis and joint replacement

The aim of the study was to estimate the tibiofemoral joint force in deep flexion to consider how the mechanical load affects the knee. We hypothesize that the joint force should not become sufficiently large to damage the joint under normal contact area, but should become deleterious to the joint under the limited contact area. Sixteen healthy knees were analyzed using a motion capture system, a force plate, a surface electromyography, and a knee model, and then tibiofemoral joint contact forces were calculated. Also, a contact stress simulation using the contact areas from the literature was performed. The peak joint contact forces (M +/- SD) were 4566 +/- 1932 N at 140 degrees in rising from full squat and 4479 +/- 1478 N at 90 degrees in rising from kneeling. Under normal contact area, the tibiofemoral contact stresses in deep flexion were less than 5 MPa and did not exceed the stress to damage the cartilage. The contact stress simulation suggests that knee prosthesis having the contact area smaller than 200 mm2 may be problematic since the contact stress in deep flexion would become larger than 21 MPa, and it would lead damage or wear of the polyethylene.     

Determination of muscle loading at the hip joint for use in pre-clinical testing

The stability of joint endoprostheses depends on the loading conditions to which the implant-bone complex is exposed. Due to a lack of appropriate muscle force data, less complex loading conditions tend to be considered in vitro. The goal of this study was to develop a load profile that better simulates the in vivo loading conditions of a "typical" total hip replacement patient and considers the interdependence of muscle and joint forces. The development of the load profile was based on a computer model of the lower extremities that has been validated against in vivo data. This model was simplified by grouping functionally similar hip muscles. Muscle and joint contact forces were computed for an average data set of up to four patients throughout walking and stair climbing. The calculated hip contact forces were compared to the average of the in vivo measured forces. The final derived load profile included the forces of up to four muscles at the instances of maximum in vivo hip joint loading during both walking and stair climbing. The hip contact forces differed by less than 10% from the peak in vivo value for a "typical" patient. The derived load profile presented here is the first that is based on validated musculoskeletal analyses and seems achievable in an in vitro test set-up. It should therefore form the basis for further standardisation of pre-clinical testing by providing a more realistic approximation of physiological loading conditions.     

Hip joint contact force in the emu (Dromaius novaehollandiae) during normal level walking

The emu is a large, (bipedal) flightless bird that potentially can be used to study various orthopaedic disorders in which load protection of the experimental limb is a limitation of quadrupedal models. An anatomy-based analysis of normal emu walking gait was undertaken to determine hip contact forces for comparison with human data. Kinematic and kinetic data captured for two laboratory-habituated emus were used to drive the model. Muscle attachment data were obtained by dissection, and bony geometries were obtained by CT scan. Inverse dynamics calculations at all major lower-limb joints were used in conjunction with optimization of muscle forces to determine hip contact forces. Like human walking gait, emu ground reaction forces showed a bimodal distribution over the course of the stance phase. Two-bird averaged maximum hip contact force was approximately 5.5 times body weight, directed nominally axially along the femur. This value is only modestly larger than optimization-based hip contact forces reported in literature for humans. The interspecies similarity in hip contact forces makes the emu a biomechanically attractive animal in which to model loading-dependent human orthopaedic hip disorders.     

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.     

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.     

Improved single- and multi-contact life-time testing of dental restorative materials using key characteristics of the human masticatory system and a force/position-controlled robotic dental wear simulator

This paper presents a new in vitro wear simulator based on spatial parallel kinematics and a biologically inspired implicit force/position hybrid controller to replicate chewing movements and dental wear formations on dental components, such as crowns, bridges or a full set of teeth. The human mandible, guided by passive structures such as posterior teeth and the two temporomandibular joints, moves with up to 6 degrees of freedom (DOF) in Cartesian space. The currently available wear simulators lack the ability to perform these chewing movements. In many cases, their lack of sufficient DOF enables them only to replicate the sliding motion of a single occlusal contact point by neglecting rotational movements and the motion along one Cartesian axis. The motion and forces of more than one occlusal contact points cannot accurately be replicated by these instruments. Furthermore, the majority of wear simulators are unable to control simultaneously the main wear-affecting parameters, considering abrasive mechanical wear, which are the occlusal sliding motion and bite forces in the constraint contact phase of the human chewing cycle. It has been shown that such discrepancies between the true in vivo and the simulated in vitro condition influence the outcome and the quality of wear studies. This can be improved by implementing biological features of the human masticatory system such as tooth compliance realized through the passive action of the periodontal ligament and active bite force control realized though the central nervous system using feedback from periodontal preceptors. The simulator described in this paper can be used for single- and multi-occlusal contact testing due to its kinematics and ability to exactly replicate human translational and rotational mandibular movements with up to 6 DOF without neglecting movements along or around the three Cartesian axes. Recorded human mandibular motion and occlusal force data are the reference inputs of the simulator. Experimental studies of wear using this simulator demonstrate that integrating the biological feature of combined force/position hybrid control in dental material testing improves the linearity and reduces the variability of results. In addition, it has been shown that present biaxially operated dental wear simulators are likely to provide misleading results in comparative in vitro/in vivo one-contact studies due to neglecting the occlusal sliding motion in one plane which could introduce an error of up to 49% since occlusal sliding motion D and volumetric wear loss V(loss) are proportional.     

Study of micromotion in modular acetabular components during gait and subluxation: a finite element investigation

Polyethylene wear after total hip arthroplasty may occur as a result of normal gait and as a result of subluxation and relocation with impact. Relocation of a subluxed hip may impart a moment to the cup creating sliding as well as compression at the cup liner interface. The purpose of the current study is to quantify, by a validated finite element model, the forces generated in a hip arthroplasty as a result of subluxation relocation and compare them to the forces generated during normal gait. The micromotion between the liner and acetabular shell was quantified by computing the sliding track and the deformation at several points of the interface. A finite element analysis of polyethylene liner stress and liner/cup micromotion in total hip arthroplasty was performed under two dynamic profiles. The first profile was a gait loading profile simulating the force vectors developed in the hip arthroplasty during normal gait. The second profile is generated during subluxation and subsequent relocation of the femoral head. The forces generated by subluxation relocation of a total hip arthroplasty can exceed those forces generated during normal gait. The induced micromotion at the cup polyethylene interface as a result of subluxation can exceed micromotion as a result of the normal gait cycle. This may play a significant role in the generation of backsided wear. Minimizing joint subluxation by restoring balance to the hip joint after arthroplasty should be explored as a strategy to minimize backsided wear.     

Postoperative Changes in In Vivo Measured Friction in Total Hip Joint Prosthesis during Walking

Loosening of the artificial cup and inlay is the most common reasons for total hip replacement failures. Polyethylene wear and aseptic loosening are frequent reasons. Furthermore, over the past few decades, the population of patients receiving total hip replacements has become younger and more active. Hence, a higher level of activity may include an increased risk of implant loosening as a result of friction-induced wear. In this study, an instrumented hip implant was used to measure the contact forces and friction moments in vivo during walking. Subsequently, the three-dimensional coefficient of friction in vivo was calculated over the whole gait cycle. Measurements were collected from ten subjects at several time points between three and twelve months postoperative. No significant change in the average resultant contact force was observed between three and twelve months postoperative. In contrast, a significant decrease of up to 47% was observed in the friction moment. The coefficient of friction also decreased over postoperative time on average. These changes may be caused by ‘running-in’ effects of the gliding components or by the improved lubricating properties of the synovia. Because the walking velocity and contact forces were found to be nearly constant during the observed period, the decrease in friction moment suggests an increase in fluid viscosity. The peak values of the contact force individually varied by 32%-44%. The friction moment individually differed much more, by 110%-129% at three and up to 451% at twelve months postoperative. The maximum coefficient of friction showed the highest individual variability, about 100% at three and up to 914% at twelve months after surgery. These individual variations in the friction parameters were most likely due to different ‘running-in’ effects that were influenced by the individual activity levels and synovia properties.     

Stair climbing is more critical than walking in pre-clinical assessment of primary stability in cementless THA in vitro

Pre-clinical testing of hip endoprostheses is a mandatory requirement before clinical release. Inadequate loading conditions may lead to lower elastic and plastic interface movements than those occurring post-operatively in vivo. This study investigated the influence of patient activity on the primary stability of cementless prostheses with a special emphasis on active simulation of muscle forces. A loading setup, based on validated musculo-skeletal analyses, was used to generate the hip contact force during walking and stair climbing by transmitting muscle forces through the femur. In addition, a loading configuration which only generated the hip contact force occurring during stair climbing at the prosthesis head was simulated. CLS prostheses were implanted in 18 composite femora and subjected to cyclical loading. The relative micro-movements at the bone-prosthesis interface were determined and appeared to be extremely sensitive to the specific patient activity. Compared to walking, stair climbing generated higher micro-movements, with pronounced axial and rotational components. Stair climbing with the femur loaded by the resultant hip contact force only exhibited a characteristic valgus tilt of the stem with significantly lower interface micro-movements than under active simulation of muscle forces. The analyses suggest that stair climbing induced the highest mechanical instability at the bone-prosthesis interface, a level which may compromise the necessary osseointegration process. Active simulation of muscle forces considerably affects the primary stability of cementless hip endoprostheses. Pre-clinical in vitro tests should therefore simulate stair climbing and include muscle activity in the assessment of initial implant stability, otherwise micro-movements may be underestimated and the primary stability overestimated.     

Quantification of the effect of cross-shear on the wear of conventional and highly cross-linked UHMWPE

A computational model has been developed to quantify the degree of cross-shear of a polyethylene pin articulating against a metallic plate, based on the direct simulation of a multidirectional pin-on-plate wear machine. The principal molecular orientation (PMO) was determined for each polymer site. The frictional work in the direction perpendicular to the PMO was assumed to produce the greatest orientation softening [Wang et al., 1997. Orientation softening in the deformation and wear of ultra-high molecular weight polyethylene. Wear 203-204, 230-241]. The cross-shear ratio (CS) was defined as the frictional work perpendicular to the PMO direction, divided by the total frictional work. Cross-shear on the pin contact surface was location specific, and of continuously changing magnitude because the direction of frictional force continuously changed due to pin rotation. The polymer pin motion was varied from a purely linear track (CS=0) up to a maximum rotation of +/-55 degrees (CS=0.254). The relationship between wear factors (K) measured experimentally and theoretically predicted CS was defined using logarithmic functions for both conventional and highly cross-linked ultra-high molecular weight polyethylene (UHMWPE). Cross-shear increased the apparent wear factor for both polyethylenes by more than fivefold compared to unidirectional wear.     

Quantification of the Contact Area at the Head-Stem Taper Interface of Modular Hip Prostheses

Corrosion of modular taper junctions of hip implants may be associated with clinical failure. Taper design parameters, as well as the intraoperatively applied assembly forces, have been proposed to affect corrosion. Fretting corrosion is related to relative interface shear motion and fluid ingress, which may vary with contact force and area. It was hypothesised in this study that assembly forces modify the extent and distribution of the surface contact area at the taper interface between a cobalt chrome head and titanium stem taper with a standard threaded surface profile. Local abrasion of a thin gold coating applied to the stem taper prior to assembly was used to determine the contact area after disassembly. Profilometry was then used to assess permanent deformation of the stem taper surface profile. With increasing assembly force (500 N, 2000 N, 4000 N and 8000 N) the number of stem taper surface profile ridges in contact with the head taper was found to increase (9.2±9.3%, 65.4±10.8%, 92.8±6.0% and 100%) and the overall taper area in contact was also found to increase (0.6±0.7%, 5.5±1.0%, 9.9±1.1% and 16.1±0.9%). Contact was inconsistently distributed over the length of the taper. An increase in plastic radial deformation of the surface ridges (-0.05±0.14 μm, 0.1±0.14 μm, 0.21±0.22 μm and 0.96±0.25 μm) was also observed with increasing assembly force. The limited contact of the taper surface ridges at lower assembly forces may influence corrosion rates, suggesting that the magnitude of the assembly force may affect clinical outcome. The method presented provides a simple and practical assessment of the contact area at the taper interface.     

Effect of heel height on in-shoe localized triaxial stresses

Abnormal and excessive plantar pressure and shear are potential risk factors for high-heeled related foot problems, such as forefoot pain, hallux valgus deformity and calluses. Plantar shear stresses could be of particular importance with an inclined supporting surface of high-heeled shoe. This study aimed to investigate the contact pressures and shear stresses simultaneously between plantar foot and high-heeled shoe over five major weightbearing regions: hallux, heel, first, second and fourth metatarsal heads, using in-shoe triaxial force transducers. During both standing and walking, peak pressure and shear stress shifted from the lateral to the medial forefoot as the heel height increased from 30 to 70mm. Heel height elevation had a greater influence on peak shear than peak pressure. The increase in peak shear was up to 119% during walking, which was about five times that of peak pressure. With increasing heel height, peak posterolateral shear over the hallux at midstance increased, whereas peak pressure at push-off decreased. The increased posterolateral shear could be a contributing factor to hallux deformity. It was found that there were differences in the location and time of occurrence between in-shoe peak pressure and peak shear. In addition, there were significant differences in time of occurrence for the double-peak loading pattern between the resultant horizontal ground reaction force peaks and in-shoe localized peak shears. The abnormal and drastic increase of in-shoe shear stresses might be a critical risk factor for shoe-related foot disorders. In-shoe triaxial stresses should therefore be considered to help in designing proper footwear.     

One-legged stance as a representative static body position for calculation of hip contact stress distribution in clinical studies

It was shown in several clinical studies that static one-legged stance may be a relevant body position to describe the loads acting at the hip. However, the stress distribution averaged during movement may better describe hip load than hip contact stress distribution in the static body position. Using data on the resultant hip force during walking taken from the measurements of Bergmann (2001), spatial distribution of contact stress over the articular surface was calculated by the HIPSTRESS method and compared with the stress distribution in one-legged stance. It is shown, that the shape of the contact stress distribution during one-legged stance closely resembled the averaged contact stress distribution during the walking cycle (Pearson's correlation coefficient R2 equals; .986; p < .001). This finding presents a link between the hypothesis that the averaged contact stress distribution during a walking cycle is crucial for cartilage development and the results of clinical studies in which the calculated distribution of contact stress in one-legged stance was successfully used to predict the clinical status of the hip.     

In vivo contact kinematics and contact forces of the knee after total knee arthroplasty during dynamic weight-bearing activities

Analysis of polyethylene component wear and implant loosening in total knee arthroplasty (TKA) requires precise knowledge of in vivo articular motion and loading conditions. This study presents a simultaneous in vivo measurement of tibiofemoral articular contact forces and contact kinematics in three TKA patients. These measurements were accomplished via a dual fluoroscopic imaging system and instrumented tibial implants, during dynamic single leg lunge and chair rising-sitting. The measured forces and contact locations were also used to determine mediolateral distribution of axial contact forces. Contact kinematics data showed a medial pivot during flexion of the knee, for all patients in the study. Average axial forces were higher for lunge compared to chair rising-sitting (224% vs. 187% body weight). In this study, we measured peak anteroposterior and mediolateral forces averaging 13.3% BW during lunge and 18.5% BW during chair rising-sitting. Mediolateral distributions of axial contact force were both patient and activity specific. All patients showed equitable medial-lateral loading during lunge but greater loads at the lateral compartment during chair rising-sitting. The results of this study may enable more accurate reproduction of in vivo loads and articular motion patterns in wear simulators and finite element models. This in turn may help advance our understanding of factors limiting longevity of TKA implants, such as aseptic loosening and polyethylene component wear, and enable improved TKA designs.