Simulation of a knee joint replacement during a gait cycle using explicit finite element analysis.

The stress distribution within the polyethylene insert of a total knee joint replacement is dependent on the kinematics, which in turn are dependent on the design of the articulating surfaces, the relative position of the components and the tension of the surrounding soft tissues. Implicit finite element analysis techniques have been used previously to examine the polyethylene stresses. However, these have essentially been static analyses and hence ignored the influence of the kinematics. The aim of this work was to use an explicit finite element approach to simulate both the kinematics and the internal stresses within a single analysis. A simulation of a total knee joint replacement subjected to a single gait cycle within a knee wear simulator was performed and the results were compared with experimental data.The predicted kinematics were in close agreement with the experimental data. Various solution-dependent parameters were found to have little influence on the predicted kinematics. The predicted stresses were found to be dependent on the mesh density. This study has shown that an explicit finite element approach is capable of predicting the kinematics and the stresses within a single analysis at relatively low computational cost.     

Computational analysis of knee joint stability following total knee arthroplasty

The overall objective of this study was to introduce knee joint power as a potential measure to investigate knee joint stability following total knee arthroplasty (TKA). Specific aims were to investigate whether weakened knee joint stabilizers cause abnormal kinematics and how it influences the knee joint kinetic (i.e., power) in response to perturbation.Patient-specific musculoskeletal models were simulated with experimental gait data from six TKA patients (baseline models). Muscle strength and ligament force parameter were reduced by up to 30% to simulate weak knee joint stabilizers (weak models). Two different muscle recruitment criteria were tested to examine whether altered muscle recruitment pattern can mask the influence of weakened stabilizers on the knee joint kinematics and kinetics. Level-walking knee joint kinematics and kinetics were calculated though force-dependent kinematic and inverse dynamic analyses. Bode analysis was then recruited to estimate the knee joint power in response to a simulated perturbation.Weak models resulted in larger anterior-posterior (A-P) displacement and internal-external (I-E) rotation compared to baseline (I-E: 18.4 ± 8.5 vs. 11.6 ± 5.7 (deg), A-P: 9.7 ± 5.6 vs. 5.5 ± 4.1 (mm)). Changes in muscle recruitment criterion however altered the results such that A-P and I-E were not notably different from baseline models. In response to the simulated perturbation, weak models versus baseline models generated a delayed power response with unbounded magnitudes. Perturbed power behavior of the knee remained unaltered regardless of the muscle recruitment criteria.In conclusion, impairment at the knee joint stabilizers may or may not lead to excessive joint motions but it notably affects the knee joint power in response to a perturbation. Whether perturbed knee joint power is associated with the patient-reported outcome requires further investigation.     

In vitro analysis of varus-valgus laxity of the knee joint: Comparison of clinical evaluation with measurements using a reference motion analysis system

IntroductionNumerous measurement devices can help clinicians during the knee examination. However, manual evaluation still remains routinely used to assess the knee laxities. The present study evaluated how accurate was a clinician for a varus-valgus stress test. We compared the clinician evaluation to the objective measurement of the knee movements during the same test session.MethodsWe studied six fresh-frozen anatomical lower limbs. The clinician performed a varus-valgus stress test in extension and at 25° flexion. The limbs were equipped with intracortical pins in femur and tibia, and spherical retro-reflecting markers were glued on the pins. Objective knees movements were measured by means of a Motion Analysis^® system (Motion Analysis Corporation, Santa Rosa, CA, USA). Two statistical analyses were performed. A single sample t-test was first used to verify the required 25? flexion. Then, a multivariate anova was performed to analyse the varus-valgus laxity under the fixed factors of measurement method and flexion of the knee.ResultsThe results for varus-valgus and total laxity of the clinician always exhibited a greater variability than objective measurements of the device. Test condition is a factor of grouping differences for Valgus and for global mediolateral laxity. Statistical analysis revealed that the objective measurement was able to show a difference between extension and 25° flexion for global mediolateral laxity, whereas the clinician was not.DiscussionThe clinician was relatively accurate in his manual evaluation. However, we demonstrated that a measurement device could clearly help clinician to exhibit differences in laxity. This can be very useful to compare a knee to itself in two successive conditions, e.g., before and after a surgery.     

Precision MRI-based joint surface and cartilage density analysis of the knee joint using rapid water-excitation sequence and semi-automatic segmentation algorithm]

The aim of this study was to analyse the precision of three-dimensional joint surface and cartilage thickness measurements in the knee, using a fast, high-resolution water-excitation sequence and a semiautomated segmentation algorithm. The knee joint of 8 healthy volunteers, aged 22 to 29 years, were examined at a resolution of 1.5 mm x 0.31 mm x 0.31 mm, with four sagittal data sets being acquired after repositioning the joint. After semiautomated segmentation with a B-spline Snake algorithm and 3D reconstruction of the patellar, femoral and tibial cartilages, the joint surface areas (triangulation), cartilage volume, and mean and maximum thickness (Euclidean distance transformation) were analysed, independently of the orientation of the sections. The precision (CV%) for the surface areas was 2.1 to 6.6%. The mean cartilage thickness and cartilage volume showed coefficients of 1.9 to 3.5% (except for the femoral condyles), the value for the medial femoral condyle being 9.1%, and for the lateral condyle 6.5%. For maximum thickness, coefficients of between 2.6 and 5.9% were found. In the present study we investigate for the first time the precision of MRI-based joint surface area measurements in the knee, and of cartilage thickness analyses in the femur. Using a selective water-excitation sequence, the acquisition time can be reduced by more than 50%. The poorer precision in the femoral condyles can be attributed to partial volume effects that occur at the edges of the joint surfaces with a sagittal image protocol. Since MRI is non-invasive, it is highly suitable for examination of healthy subjects (generation of individual finite element models, analysis of functional adaptation to mechanical stimulation, measurement of cartilage deformation in vivo) and as a diagnostic tool for follow-up, indication for therapy, and objective evaluation of new therapeutic agents in osteoarthritis.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

Repeatability of gait data using a functional hip joint centre and a mean helical knee axis

Repeatability of traditional kinematic and kinetic models is affected by the ability to accurately locate anatomical landmarks (ALs) to define joint centres and anatomical coordinate systems. Numerical methods that define joint centres and axes of rotation independent of ALs may also improve the repeatability of kinematic and kinetic data. The purpose of this paper was to compare the repeatability of gait data obtained from two models, one based on ALs (AL model), and the other incorporating a functional method to define hip joint centres and a mean helical axis to define knee joint flexion/extension axes (FUN model). A foot calibration rig was also developed to define the foot segment independent of ALs. The FUN model produced slightly more repeatable hip and knee joint kinematic and kinetic data than the AL model, with the advantage of not having to accurately locate ALs. Repeatability of the models was similar comparing within-tester sessions to between-tester sessions. The FUN model may also produce more repeatable data than the AL model in subject populations where location of ALs is difficult. The foot calibration rig employed in both the AL and FUN model provided an easy alternative to define the foot segment and obtain repeatable data, without accurately locating ALs on the foot.     

Kinematic and dynamic analysis of an anatomically based knee joint

This paper presents a knee-joint model to provide a better understanding on the interaction between natural joints and artificial mechanisms for design and control of rehabilitation exoskeletons. The anatomically based knee model relaxes several commonly made assumptions that approximate a human knee as engineering pin-joint in exoskeleton design. Based on published MRI data, we formulate the kinematics of a knee-joint and compare three mathematical approximations; one model bases on two sequential circles rolling a flat plane; and the other two are mathematically differentiable ellipses-based models with and without sliding at the contact. The ellipses-based model taking sliding contact into accounts shows that the rolling-sliding ratio of a knee-joint is not a constant but has an average value consistent with published measurements. This knee-joint kinematics leads to a physically more accurate contact-point trajectory than methods based on multiple circles or lines, and provides a basis to derive a knee-joint kinetic model upon which the effects of a planar exoskeleton mechanism on the internal joint forces and torque during flexion can be numerically investigated. Two different knee-joint kinetic models (pin-joint approximation and anatomically based model) are compared against a condition with no exoskeleton. The leg and exoskeleton form a closed kinematic chain that has a significant effect on the joint forces in the knee. Human knee is more tolerant than pin-joint in negotiating around a singularity but its internal forces increase with the exoskeleton mass-to-length ratio. An oversimplifying pin-joint approximation cannot capture the finite change in the knee forces due to the singularity effect.     

Measurement of joint capsule tissue loading in the cat knee using calibrated mechanoreceptors

The vertical loading in the posterior capsule of the cat knee has been measured while the knee is rotated into hyperextension. Tissue loading was determined using a previously verified model of the capsule that represents its upper edge as a catenary suspension cable. Tensile loads in the cable were measured using the discharge of mechanoreceptive sensory neurons that had been calibrated as load sensors. The results revealed that the capsule is very lightly loaded in extension rotations. Less than 4% of the applied moment is sustained by the capsule.     

Uncertainty of knee joint muscle activity during knee joint torque exertion: the significance of controlling adjacent joint torque.

In the single-joint torque exertion task, which has been widely used to control muscle activity, only the relevant joint torque is specified. However, the neglect of the neighboring joint could make the procedure unreliable, considering our previous result that even monoarticular muscle activity level is indefinite without specifying the adjacent joint torque. Here we examined the amount of hip joint torque generated with knee joint torque and its influence on the activity of the knee joint muscles. Twelve healthy subjects were requested to exert various levels of isometric knee joint torque. The knee and hip joint torques were obtained by using a custom-made device. Because no information about hip joint torque was provided to the subjects, the hip joint torque measured here was a secondary one associated with the task. The amount of hip joint torque varied among subjects, indicating that they adopted various strategies to achieve the task. In some subjects, there was a considerable internal variability in the hip joint torque. Such variability was not negligible, because the knee joint muscle activity level with respect to the knee joint torque, as quantified by surface electromyography (EMG), changed significantly when the subjects were requested to change the strategy. This change occurred in a very systematic manner: in the case of the knee extension, as the hip flexion torque was larger, the activity of mono- and biarticular knee extensors decreased and increased, respectively. These results indicate that the conventional single knee joint torque exertion has the drawback that the intersubject and/or intertrial variability is inevitable in the relative contribution among mono- and biarticular muscles because of the uncertainty of the hip joint torque. We discuss that the viewpoint that both joint torques need to be considered will bring insights into various controversial problems such as the shape of the EMG-force relationship, neural factors that help determine the effect of muscle strength training, and so on.     

Effect of knee joint flexion and femur rotation on retropatellar contact of the human knee joint]

The aim of the present study was to evaluate retropatellar contact characteristics at different angles of flexion of the knee joint. To this end, 6 cadaveric legs were examined using pressure sensitive film (Fuji Prescale type "super low") at angles of flexion of 45 degrees, 60 degrees, 90 degrees and 120 degrees both in neutral rotation and 10 degrees internal and external rotation of the femur in the same knee joints. A force of 140 N was applied to both the vastus medialis and lateralis, and a comparison made with a medially and a laterally dominating muscle force. The contact areas decreased with increasing angles of flexion. The medially dominating muscle traction increased the contact area. Comparison between internal and external rotation revealed a decrease in contact area on internal rotation. The pressure measurements were comparable in all loading situations. Comparison between neutral and medial traction revealed significant differences in contact area, pressure and force. The influence of femoral rotation showed no significant difference. A comparison of the different angles of flexion revealed only few significant differences. To prevent the development of retropatellar arthrosis, maximum contact areas are necessary. The study has shown an advantage for medially dominating muscle traction, and external rotation of the femur.     

Recruitment of knee joint ligaments

On the basis of earlier reported data on the in vitro kinematics of passive knee-joint motions of four knee specimens, the length changes of ligament fiber bundles were determined by using the points of insertion on the tibia and femur. The kinematic data and the insertions of the ligaments were obtained by using Roentgenstereophotogrammetry. Different fiber bundles of the anterior and posterior cruciate ligaments and the medial and lateral collateral ligaments were identified. On the basis of an assumption for the maximal strain of each ligament fiber bundle during the experiments, the minimal recruitment length and the probability of recruitment were defined and determined. The motions covered the range from extension to 95 degrees flexion and the loading conditions included internal or external moments of 3 Nm and anterior or posterior forces of 30 N. The ligament length and recruitment patterns were found to be consistent for some ligament bundles and less consistent for other ligament bundles. The most posterior bundle of each ligament was recruited in extension and the lower flexion angles, whereas the anterior bundle was recruited for the higher flexion angles. External rotation generally recruited the collateral ligaments, while internal rotation recruited the cruciate ligaments. However, the anterior bundle of the posterior cruciate ligament was recruited with external rotation at the higher flexion angles. At the lower flexion angles, the anterior cruciate and the lateral collateral ligaments were recruited with an anterior force. The recruitment of the posterior cruciate ligament with a posterior force showed that neither its most anterior nor its most posterior bundle was recruited at the lower flexion angles. Hence, the posterior restraint must have been provided by the intermediate fiber bundles, which were not considered in the experiment. At the higher flexion angles, the anterior bundles of the anterior cruciate ligament and the posterior cruciate ligament were found to be recruited with anterior and posterior forces, respectively. The minimal recruitment length and the recruitment probability of ligament fiber bundles are useful parameters for the evaluation of ligament length changes in those experiments where no other method can be used to determine the zero strain lengths, ligament strains and tensions.     

The finite helical axis of the knee joint (a non-invasive in vivo study using fast-PC MRI)

An understanding of the in vivo knee joint kinematics is critical for the further improvement and validation of knee joint models and for the development of better surgical and rehabilitative protocols. Unfortunately, most studies exploring the finite helical axis (FHA) tend to produce excellent qualitative results, but quantitative results are often lacking. Thus, the purpose of this study was to non-invasively and in vivo quantify the tibiofemoral FHA in a relatively large normal population during volitional knee extension using fast-PC MRI, to report the data relative to consistent coordinate systems (making it available for modeling input, experimental comparison and for device design), to determine the variability of the FHA, to investigate the screw home mechanism and to test the hypothesis that knee joint kinematics are independent of gender. Intra- and inter-subject repeatability was excellent. The intra- (inter-) subject repeatability of the FHA orientation in the frontal and axial planes was 1.8% (3.3%) and 3.7% (6.0%) of the average value, respectively. At the beginning of extension, the FHA was directed laterally and slightly superiorly and at the end of extension, it was directed in the lateral-inferior direction, indicative of the screw-home mechanism. The FHA location was not fixed during extension. There was small, but significant differences in all FHA parameters between genders and normalizing positional data relative to epicondylar width helped to reduce this difference. The data obtained in the current study forms an excellent base for future knee joint modeling and clinical studies.