Evaluating knee replacement mechanics during ADL with PID-controlled dynamic finite element analysis

Validated computational knee simulations are valuable tools for design phase development of knee replacement devices. Recently, a dynamic finite element (FE) model of the Kansas knee simulator was kinematically validated during gait and deep flexion cycles. In order to operate the computational simulator in the same manner as the experiment, a proportional–integral–derivative (PID) controller was interfaced with the FE model to control the quadriceps actuator excursion and produce a target flexion profile regardless of implant geometry or alignment conditions. The controller was also expanded to operate multiple actuators simultaneously in order to produce in vivo loading conditions at the joint during dynamic activities. Subsequently, the fidelity of the computational model was improved through additional muscle representation and inclusion of relative hip–ankle anterior–posterior (A–P) motion. The PID-controlled model was able to successfully recreate in vivo loading conditions (flexion angle, compressive joint load, medial–lateral load distribution or varus-valgus torque, internal–external torque, A–P force) for deep knee bend, chair rise, stance-phase gait and step-down activities.     

Gait symmetry - A valid parameter for pre and post planning for total knee arthroplasty

Objectives:We aimed to determine whether GS can help to plan and rearrange the treated side by using IMUs to measure the joint angle of the hip, knee, and ankle. We hypothesized that the kinematics in healthy individuals for both sides are approximately equal during walking.Methods:IMUs were used to measure the joint angles of 25 healthy participants during walking. The participants performed the 10-meter walk test. The normalized symmetry index (SI_norm) was used to calculate the symmetry of joint angles for the hip, knee, and ankle throughout the gait cycle.Results:The SI_norm demonstrated high symmetry between both legs; and the ranges were -1.5% and 1.1% for the hip, -3.0% and 3.1% for the knee, and -12% and 9.2% for the ankle joint angle throughout the gait cycle.Conclusion:The SI_norm provides strong information that can be helpful in the planning process for the surgeries. Further, the IMUs system gives the possibility to measure the patients before their surgeries and use their data to plan and rearrange for the operated side.     

Gait posture estimation using wearable acceleration and gyro sensors

A method for gait analysis using wearable acceleration sensors and gyro sensors is proposed in this work. The volunteers wore sensor units that included a tri-axis acceleration sensor and three single axis gyro sensors. The angular velocity data measured by the gyro sensors were used to estimate the translational acceleration in the gait analysis. The translational acceleration was then subtracted from the acceleration sensor measurements to obtain the gravitational acceleration, giving the orientation of the lower limb segments. Segment orientation along with body measurements were used to obtain the positions of hip, knee, and ankle joints to create stick figure models of the volunteers. This method can measure the three-dimensional positions of joint centers of the hip, knee, and ankle during movement. Experiments were carried out on the normal gait of three healthy volunteers. As a result, the flexion–extension (F–E) and the adduction–abduction (A–A) joint angles of the hips and the flexion–extension (F–E) joint angles of the knees were calculated and compared with a camera motion capture system. The correlation coefficients were above 0.88 for the hip F–E, higher than 0.72 for the hip A–A, better than 0.92 for the knee F–E. A moving stick figure model of each volunteer was created to visually confirm the walking posture. Further, the knee and ankle joint trajectories in the horizontal plane showed that the left and right legs were bilaterally symmetric.     

Influence of the 3D inverse dynamic method on the joint forces and moments during gait

The joint forces and moments are commonly used in gait analysis. They can be computed by four different 3D inverse dynamic methods proposed in the literature, either based on vectors and Euler angles, wrenches and quaternions, homogeneous matrices, or generalized coordinates and forces. In order to analyze the influence of the inverse dynamic method, the joint forces and moments were computed during gait on nine healthy subjects. A ratio was computed between the relative dispersions (due to the method) and the absolute amplitudes of the gait curves. The influence of the inverse dynamic method was negligible at the ankle (2%) but major at the knee and the hip joints (40%). This influence seems to be due to the dynamic computation rather than the kinematic computation. Compared to the influence of the joint center location, the body segment inertial parameter estimation, and more, the influence of the inverse dynamic method is at least of equivalent importance. This point should be confirmed with other subjects, possibly pathologic, and other movements.     

Invariant hip moment pattern while walking with a robotic hip exoskeleton

Robotic lower limb exoskeletons hold significant potential for gait assistance and rehabilitation; however, we have a limited understanding of how people adapt to walking with robotic devices. The purpose of this study was to test the hypothesis that people reduce net muscle moments about their joints when robotic assistance is provided. This reduction in muscle moment results in a total joint moment (muscle plus exoskeleton) that is the same as the moment without the robotic assistance despite potential differences in joint angles. To test this hypothesis, eight healthy subjects trained with the robotic hip exoskeleton while walking on a force-measuring treadmill. The exoskeleton provided hip flexion assistance from approximately 33% to 53% of the gait cycle. We calculated the root mean squared difference (RMSD) between the average of data from the last 15 min of the powered condition and the unpowered condition. After completing three 30-min training sessions, the hip exoskeleton provided 27% of the total peak hip flexion moment during gait. Despite this substantial contribution from the exoskeleton, subjects walked with a total hip moment pattern (muscle plus exoskeleton) that was almost identical and more similar to the unpowered condition than the hip angle pattern (hip moment RMSD 0.027, angle RMSD 0.134, p<0.001). The angle and moment RMSD were not different for the knee and ankle joints. These findings support the concept that people adopt walking patterns with similar joint moment patterns despite differences in hip joint angles for a given walking speed.     

Discriminating age and disability effects in locomotion: neuromuscular adaptations in musculoskeletal pathology.

We identified biomechanical variables indicative of lower extremity dysfunction, distinct from age-related gait adaptations, and examined interrelationships among these variables to better understand the neuromuscular adaptations in gait. Sagittal plane ankle, knee, and hip peak angles, moments, and powers and spatiotemporal parameters were acquired during preferred-speed gait in 120 subjects: 45 healthy young, 37 healthy elders, and 38 elders with functional limitations due to lower extremity musculoskeletal pathology, primarily arthritis. Multiple analysis of covariance with discriminate analysis, adjusted for gait speed, was used to identify the variables discriminating groups. Correlation analysis was used to explore interrelationships among these variables within each group. Healthy elders were discriminated (sensitivity 76%, specificity 82%) from young adults via decreased late-stance ankle plantar flexion angle, increased late-stance knee power absorption, and early-stance hip extensor power generation. Disabled elders were discriminated (sensitivity 74%, specificity 73%) from healthy elders via decreased late-stance ankle plantar flexor moment and power generation, increased early-stance ankle dorsiflexor moment, and late-stance hip flexor moment and power absorption. Relationships among variables showed a higher degree of coupling for the disabled elders compared with the healthy groups, suggesting a reduced ability to alter motor strategies. Our data suggest that, beyond age-related changes, elders with lower extremity dysfunction rely excessively on passive action of hip flexors to provide propulsion in late stance and contralateral ankle dorsiflexors to enhance stability. These findings support a growing body of evidence that gait changes with age and disablement have a neuromuscular basis, which may be informative in a motor control framework for physical therapy interventions.     

New assistive walker improved local dynamic stability in young healthy adults

In this study, we investigated the effect of walker type on gait pattern characteristics comparing normal gait (NG), gait with a regular walker (RW), and gait with a newly developed walker with vertical moveable handlebars, the Crosswalker (CW).Partial weight bearing (PWB) of the feet, peak joint angles and largest Lyapunov exponent (λ_max) of the lower extremities (hip, knee, ankle) in the sagittal plane, and gait parameters (gait velocity, stride length, cadence, stride duration) were determined for 18 healthy young adults performing 10 walking trials for each walking condition. Assistive gait with the CW improved local dynamic stability in the lower extremities (hip, knee, ankle) compared with RW and was not significantly different from NG. However, peak joint angles and stride characteristics in CW were different from NG. The PWB on the feet was lower with the RW (70.3%) compared to NG (82.8%) and CW (80.9%). This improved stability may be beneficial for the elderly and patients with impaired gait. However, increased PWB is not beneficial for patients during the early stages of rehabilitation.     

Effect of boot shaft stiffness on stability joint energy and muscular co-contraction during walking on uneven surface

Increased boot shaft stiffness may have a noticeable impact on the range of motion of the ankle joint. Therefore, the ability of the ankle joint to generate power for propulsion might be impaired. This might result in compensatory changes at the knee and hip joint. Besides, adaptability of the subtalar joint to uneven surface might be reduced, which could in turn affect stability. The aim of the study was therefore to investigate the influence of boot shaft stiffness on biomechanical gait parameters.Fifteen healthy young adults walked over coarse gravel wearing two different hiking boots that differed by 50% in passive shaft stiffness. Leg kinematics, kinetics and electromyography were measured. Gait velocity and indicators for stability were not different when walking with the hard and soft boot shaft over the gravel surface. However, the hard boot shaft decreased the ankle range of motion as well as the eccentric energy absorbed at the ankle joint. As a consequence, compensatory changes at the knee joint were observed. Co-contraction was increased, and greater eccentric energy was absorbed. Therefore, the efficiency of gait with hard boots might be decreased and joint loading at the knee might be increased, which might cause early fatigue of knee muscles during walking or hiking. The results of this study suggest that stiffness and blocking of joint motion at the ankle should not be equated with safety. A trade-off between lateral stiffness and free natural motion of the ankle joint complex might be preferable.     

Identification of passive elastic joint moment-angle relationships in the lower extremity

The purpose of this study was to develop a method for identifying subject-specific passive elastic joint moment-angle relationships in the lower extremity, which could subsequently be used to estimate passive contributions to joint kinetics during gait. Twenty healthy young adults participated in the study. Subjects were positioned side-lying with their dominant limb supported on a table via low-friction carts. A physical therapist slowly manipulated the limb through full sagittal hip, knee, and ankle ranges of motion using two hand-held 3D load cells. Lower extremity kinematics, measured with a passive marker motion capture system, and load cell readings were used to compute joint angles and associated passive joint moments. We formulated a passive joint moment-angle model that included eight exponential functions to account for forces generated via the passive stretch of uni-articular structures and bi-articular muscles. Model parameters were estimated for individual subjects by minimizing the sum of squared errors between model predicted and experimentally measured moments. The model predictions closely replicated measured joint moments with average root-mean-squared errors of 2.5, 1.4, and 0.7 Nm about the hip, knee, and ankle respectively. We show that the models can be coupled with gait kinematics to estimate passive joint moments during walking. Passive hip moments were substantial from terminal stance through initial swing, with energy being stored as the hip extended and subsequently returned during pre- and initial swing. We conclude that the proposed methodology could provide quantitative insights into the potentially important role that passive mechanisms play in both normal and abnormal gait.     

On the estimation of joint kinematics during gait.

In gait analysis, the concepts of Euler and helical (screw) angles are used to define the three-dimensional relative joint angular motion of lower extremities. Reliable estimation of joint angular motion depends on the accurate definition and construction of embedded axes within each body segment. In this paper, using sensitivity analysis, we quantify the effects of uncertainties in the definition and construction of embedded axes on the estimation of joint angular motion during gait. Using representative hip and knee motion data from normal subjects and cerebral palsy patients, the flexion-extension axis is analytically perturbed +/- 15 degrees in 5 degrees steps from a reference position, and the joint angles are recomputed for both Euler and helical angle definitions. For the Euler model, hip and knee flexion angles are relatively unaffected while the ab/adduction and rotation angles are significantly affected throughout the gait cycle. An error of 15 degrees in the definition of flexion-extension axis gives rise to maximum errors of 8 and 12 degrees for the ab/adduction angle, and 10-15 degrees for the rotation angles at the hip and knee, respectively. Furthermore, the magnitude of errors in ab/adduction and rotation angles are a function of the flexion angle. The errors for the ab/adduction angles increase with increasing flexion angle and for the rotation angle, decrease with increasing flexion angle. In cerebral palsy patients with flexed knee pattern of gait, this will result in distorted estimation of ab/adduction and rotation. For the helical model, similar results are obtained for the helical angle and associated direction cosines.(ABSTRACT TRUNCATED AT 250 WORDS)     

Influence of joint constraints on lower limb kinematics estimation from skin markers using global optimization

In order to obtain the lower limb kinematics from skin-based markers, the soft tissue artefact (STA) has to be compensated. Global optimization (GO) methods rely on a predefined kinematic model and attempt to limit STA by minimizing the differences between model predicted and skin-based marker positions. Thus, the reliability of GO methods depends directly on the chosen model, whose influence is not well known yet.This study develops a GO method that allows to easily implement different sets of joint constraints in order to assess their influence on the lower limb kinematics during gait. The segment definition was based on generalized coordinates giving only linear or quadratic joint constraints. Seven sets of joint constraints were assessed, corresponding to different kinematic models at the ankle, knee and hip: SSS, USS, PSS, SHS, SPS, UHS and PPS (where S, U and H stand for spherical, universal and hinge joints and P for parallel mechanism). GO was applied to gait data from five healthy males.Results showed that the lower limb kinematics, except hip kinematics, knee and ankle flexion–extension, significantly depend on the chosen ankle and knee constraints. The knee parallel mechanism generated some typical knee rotation patterns previously observed in lower limb kinematic studies. Furthermore, only the parallel mechanisms produced joint displacements.Thus, GO using parallel mechanism seems promising. It also offers some perspectives of subject-specific joint constraints.     

3D joint dynamics analysis of healthy children's gait

The 3D joint moments and 2D joint powers have been largely explored in the literature of healthy children's gait, in particular to compare them with pathologic subjects’ gait. However, no study reported on 3D joint power in children which could be due to the difficulties in interpreting the results. Recently, the analysis of the 3D angle between the joint moment and the joint angular velocity vectors has been proposed in order to help 3D joint power interpretation.Our hypothesis is that this 3D angle may help in characterizing the level of gait maturation. The present study explores 3D joint moments, 3D joint power and the proposed 3D angle for both children's and adults’ gaits to highlight differences in the strategies used. The results seem to confirm that children have an alternative strategy of mainly ankle stabilization and hip propulsion compared to the adults’ strategy of mainly ankle resistance and propulsion and hip stabilization.In the future, the same 3D angle analysis should be applied to different age groups for better describing the evolution of the 3D joint dynamic strategies during the growth.     

Evaluation of the kinematic parameters of normal-paced gait in subjects with gonarthrosis and the influence of gonarthrosis on the function of the ankle joint and hip joint

The aim of the study was to assess the variability of parameters characterising the gait of persons suffering from degenerative changes of the knee joint and their influence on the ankle and hip joints. The values of the angular changes in the knee, ankle and hip joints in the three planes of motion were assessed. Locomotion tests were performed on 27 persons, aged between 60 and 74, using Vicon 250, the three-dimensional analysis system. The sharpest deviations from the results of the control group were revealed in the transverse and frontal planes. Degenerative knee joint disease has changed the gait stereotype causing a reduction in the economy of gonarthrosis patients' locomotion, the influence of the disease on the function of the neighbouring joints is also distinctly marked.     

Comparison of hierarchical and six degrees-of-freedom marker sets in analyzing gait kinematics

The objective of this study was to determine how marker spacing, noise, and joint translations affect joint angle calculations using both a hierarchical and a six degrees-of-freedom (6DoF) marker set. A simple two-segment model demonstrates that a hierarchical marker set produces biased joint rotation estimates when sagittal joint translations occur whereas a 6DoF marker set mitigates these bias errors with precision improving with increased marker spacing. These effects were evident in gait simulations where the 6DoF marker set was shown to be more accurate at tracking axial rotation angles at the hip, knee, and ankle.     

A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb

Concurrent multiscale simulation strategies are required in computational biomechanics to study the interdependence between body scales. However, detailed finite element models rarely include muscle recruitment due to the computational burden of both the finite element method and the optimization strategies widely used to estimate muscle forces. The aim of this study was twofold: first, to develop a computationally efficient muscle force prediction strategy based on proportional-integral-derivative (PID) controllers to track gait and chair rise experimental joint motion with a finite element musculoskeletal model of the lower limb, including a deformable knee representation with 12 degrees of freedom; and, second, to demonstrate that the inclusion of joint-level deformability affects muscle force estimation by using two different knee models and comparing muscle forces between the two solutions. The PID control strategy tracked experimental hip, knee, and ankle flexion/extension with root mean square errors below 1°, and estimated muscle, contact and ligament forces in good agreement with previous results and electromyography signals. Differences up to 11% and 20% in the vasti and biceps femoris forces, respectively, were observed between the two knee models, which might be attributed to a combination of differing joint contact geometry, ligament behavior, joint kinematics, and muscle moment arms. The tracking strategy developed in this study addressed the inevitable tradeoff between computational cost and model detail in musculoskeletal simulations and can be used with finite element musculoskeletal models to efficiently estimate the interdependence between muscle forces and tissue deformation.     

Joint forces and torques when walking in shallow water

This study reports for the first time an estimation of the internal net joint forces and torques on adults' lower limbs and pelvis when walking in shallow water, taking into account the drag forces generated by the movement of their bodies in the water and the equivalent data when they walk on land. A force plate and a video camera were used to perform a two-dimensional gait analysis at the sagittal plane of 10 healthy young adults walking at comfortable speeds on land and in water at a chest-high level. We estimated the drag force on each body segment and the joint forces and torques at the ankle, knee, and hip of the right side of their bodies using inverse dynamics. The observed subjects' apparent weight in water was about 35% of their weight on land and they were about 2.7 times slower when walking in water. When the subjects walked in water compared with walking on land, there were no differences in the angular displacements but there was a significant reduction in the joint torques which was related to the water's depth. The greatest reduction was observed for the ankle and then the knee and no reduction was observed for the hip. All joint powers were significantly reduced in water. The compressive and shear joint forces were on average about three times lower during walking in water than on land. These quantitative results substantiate the use of water as a safe environment for practicing low-impact exercises, particularly walking.     

Gait and neuromuscular pattern changes are associated with differences in knee osteoarthritis severity levels

Knee osteoarthritis (OA) is a multifactoral, progressive disease process of the musculoskeletal system. Mechanical factors have been implicated in the progression of knee OA, but the role of altered joint mechanics and neuromuscular control strategies in progressive mechanisms of the disease have not been fully explored. Previous biomechanical studies of knee OA have characterized changes in joint kinematics and kinetics with the disease, but it has been difficult to determine if these biomechanical changes are involved in the development of disease, are in response to degenerative changes in the joint, or are compensatory mechanisms in response to these degenerative changes or other related factors as joint pain. The goal of this study was to explore the association between biomechanical changes and knee OA severity in an effort to understand the changing role of biomechanical factors in the progression of knee OA. A three-group cross-sectional model was used that included asymptomatic subjects, subjects clinically diagnosed with moderate knee OA and severe knee OA subjects just prior to total joint replacement surgery. Principal component analysis and discriminant analysis were used to determine the combinations of electromyography, kinematic and kinetic waveform pattern changes at the knee, hip and ankle joints during gait that optimally separated the three levels of severity. Different biomechanical mechanisms were important in discriminating between severity levels. Changes in knee and hip kinetic patterns and rectus femoris activation were important in separating the asymptomatic and moderate OA gait patterns. In contrast, changes in knee kinematics, hip and ankle kinetics and medial gastrocnemius activity were important in discriminating between the moderate and severe OA gait patterns.     

The passive elastic moment at the knee and its influence on human gait

The elastic component of the passive moment at the knee was measured in situ. The force needed to manually range the knee from approximately 90 degrees of flexion to full extension was measured. Hip and ankle angle were held fixed. The passive knee moment, computed from the force and knee angle data, was compared to the total knee moment required for normal gait. This comparison suggested that the passive moment can contribute a significant portion of the total joint moment during some phases of the gait cycle.     

Paired nonlinear behavior of active and passive joint torques associated with preparation for walk-to-run gait transition

This study investigated the lower extremity torque's active and passive features during the walk-to-run gait transition with continuously increased walking speed. Fourteen volunteers participated in the experiment. Kinematic and kinetic data were collected synchronously. Five strides leading up the gait transition were examined. Peaks of the passive (e.g., contact) and active (e.g., generalized muscle torques), along with net joint torque, and time to peak torques exhibited significant differences at the last stride before gait transition, compared to the first four strides, at the ankle, knee, and hip joints, respectively. Selected peak joint active and passive torques showed significant and opposite trends at critical events within a stride cycle: such ankle joint right after heel-contact, knee joint during weight acceptance, and both hip and knee joints right before toe-off. The magnitude and the corresponding time to active and passive peak torque changed in a nonlinear pattern before the transition from walk to run. The lower extremity segment-interaction during gait transition appeared to be an active reorganization exemplified by the interaction between the lower extremity's active and passive torque components.     

Empirical evaluation of gastrocnemius and soleus function during walking

Distinguishing gastrocnemius and soleus muscle function is relevant for treating gait disorders in which abnormal plantarflexor activity may contribute to pathological movement patterns. Our objective was to use experimental and computational analysis to determine the influence of gastrocnemius and soleus activity on lower limb movement, and determine if anatomical variability of the gastrocnemius affected its function. Our hypothesis was that these muscles exhibit distinct functions, with the gastrocnemius inducing limb flexion and the soleus inducing limb extension. To test this hypothesis, the gastrocnemius or soleus of 20 healthy participants was electrically stimulated for brief periods (90 ms) during mid- or terminal stance of a random gait cycle. Muscle function was characterized by the induced change in sagittal pelvis, hip, knee, and ankle angles occurring during the 200 ms after stimulation onset. Results were corroborated with computational forward dynamic gait models, by perturbing gastrocnemius or soleus activity during similar portions of the gait cycle. Mid- and terminal stance gastrocnemius stimulation induced posterior pelvic tilt, hip flexion and knee flexion. Mid-stance gastrocnemius stimulation also induced ankle dorsiflexion. In contrast mid-stance soleus stimulation induced anterior pelvic tilt, knee extension and plantarflexion, while late-stance soleus stimulation induced relatively little change in motion. Model predictions of induced hip, knee, and ankle motion were generally in the same direction as those of the experiments, though the gastrocnemius? results were shown to be quite sensitive to its knee-to-ankle moment arm ratio.