The influence and biomechanical role of cartilage split line pattern on tibiofemoral cartilage stress distribution during the stance phase of gait

This study presents an evaluation of the role that cartilage fibre ‘split line’ orientation plays in informing femoral cartilage stress patterns. A two-stage model is presented consisting of a whole knee joint coupled to a tissue-level cartilage model for computational efficiency. The whole joint model may be easily customised to any MRI or CT geometry using free-form deformation. Three ‘split line’ patterns (medial–lateral, anterior–posterior and random) were implemented in a finite element model with constitutive properties referring to this ‘split line’ orientation as a finite element fibre field. The medial–lateral orientation was similar to anatomy and was derived from imaging studies. Model predictions showed that ‘split lines’ are formed along the line of maximum principal strains and may have a biomechanical role of protecting the cartilage by limiting the cartilage deformation to the area of higher cartilage thickness.     

In vivo biomechanical behavior of the human heel pad during the stance phase of gait

A technique is introduced for simultaneous measurements of the heel pad tissue deformation and the heel–ground contact stresses developing during the stance phase of gait. Subjects walked upon a gait platform integrating the contact pressure display optical method for plantar pressure measurements and a digital radiographic fluoroscopy system for skeletal and soft tissue motion recording. Clear images of the posterior-plantar aspect of the calcaneus and enveloping soft tissues were obtained simultaneously with the pressure distribution under the heel region throughout the stance phase of gait. The heel pad was shown to undergo a rapid compression during initial contact and heel strike, reaching a strain of 0.39±0.05 in about 150 ms. The stress–strain relation of the heel pad was shown to be highly non-linear, with a compression modulus of 105±11 kPa initially and 306±16 kPa at 30% strain. The energy dissipation during heel strike was evaluated to be 17.8±0.8%. The present technique is useful for biomechanical as well as clinical evaluation of the stress–strain and energy absorption characteristics of the heel pad in vivo, during natural gait.     

Tibiofemoral kinematics and condylar motion during the stance phase of gait

Accurate knowledge of the dynamic knee motion in-vivo is instrumental for understanding normal and pathological function of the knee joint. However, interpreting motion of the knee joint during gait in other than the sagittal plane remains controversial. In this study, we utilized the dual fluoroscopic imaging technique to investigate the six-degree-of-freedom kinematics and condylar motion of the knee during the stance phase of treadmill gait in eight healthy volunteers at a speed of 0.67 m/s. We hypothesized that the 6DOF knee kinematics measured during gait will be different from those reported for non-weightbearing activities, especially with regards to the phenomenon of femoral rollback. In addition, we hypothesized that motion of the medial femoral condyle in the transverse plane is greater than that of the lateral femoral condyle during the stance phase of treadmill gait. The rotational motion and the anterior–posterior translation of the femur with respect to the tibia showed a clear relationship with the flexion–extension path of the knee during the stance phase. Additionally, we observed that the phenomenon of femoral rollback was reversed, with the femur noted to move posteriorly with extension and anteriorly with flexion. Furthermore, we noted that motion of the medial femoral condyle in the transverse plane was greater than that of the lateral femoral condyle during the stance phase of gait (17.4±2.0 mm vs. 7.4±6.1 mm, respectively; p<0.01). The trend was opposite to what has been observed during non-weightbearing flexion or single-leg lunge in previous studies. These data provide baseline knowledge for the understanding of normal physiology and for the analysis of pathological function of the knee joint during walking. These findings further demonstrate that knee kinematics is activity-dependent and motion patterns of one activity (non-weightbearing flexion or lunge) cannot be generalized to interpret a different one (gait).     

Finite element analysis of a total ankle replacement during the stance phase of gait

Total ankle replacement (TAR) designs have still several important issues to be addressed before the treatment becomes fully acceptable clinically. Very little is known about the performance, in terms of the contact pressures and kinematics of TAR when subjected to daily activities such as level gait. For this purpose, an explicit finite element model of a novel 3-component TAR was developed, which incorporated a previously validated mechanical model of the ankle ligament apparatus. The intermediate mobile polyethylene meniscal bearing was modelled as an elastic-plastic continuum while the articulating surfaces of the tibial and talar metal components as rigid bodies. Overall kinematics, contact pressures and ligament forces were analysed during passive, i.e. virtually unloaded, and active, i.e. stance phase of gait, conditions. Simulation of passive motion predicted similar kinematics as reported previously in an analytical four-bar linkage model. The meniscal bearing was observed to move 5.6 mm posteriorly during the simulated stance and the corresponding antero-posterior displacement of the talar component was 8.3 mm. The predicted pattern and the amount (10.6 degrees ) of internal-external rotation of the ankle complex were found to be in good agreement with corresponding in vivo measurements on normal ankles. A peak contact pressure of 16.8 MPa was observed, with majority of contact pressures below 10 MPa. For most ligaments, reaction forces remain within corresponding physiological ranges. A first realistic representation of the biomechanical behaviour of the human ankle when replaced by prosthetic joints is provided. The applied methodology can potentially be applied to other TAR designs.     

A finite element model of the lower limb during stance phase of gait cycle including the muscle forces

Background

Results of finite element (FE) analyses can give insight into musculoskeletal diseases if physiological boundary conditions, which include the muscle forces during specific activities of daily life, are considered in the FE modelling. So far, many simplifications of the boundary conditions are currently made. This study presents an approach for FE modelling of the lower limb for which muscle forces were included.

Methods

The stance phase of normal gait was simulated. Muscle forces were calculated using a musculoskeletal rigid body (RB) model of the human body, and were subsequently applied to a FE model of the lower limb. It was shown that the inertial forces are negligible during the stance phase of normal gait. The contact surfaces between the parts within the knee were modelled as bonded. Weak springs were attached to the distal tibia for numerical reasons.

Results

Hip joint reaction forces from the RB model and those from the FE model were similar in magnitude with relative differences less than 16%. The forces of the weak spring were negligible compared to the applied muscle forces. The maximal strain was 0.23% in the proximal region of the femoral diaphysis and 1.7% in the contact zone between the tibia and the fibula.

Conclusions

The presented approach based on FE modelling by including muscle forces from inverse dynamic analysis of musculoskeletal RB model can be used to perform analyses of the lower limb with very realistic boundary conditions. In the present form, this model can be used to better understand the loading, stresses and strains of bones in the knee area and hence to analyse osteotomy fixation devices.

     

Crouched postures reduce the capacity of muscles to extend the hip and knee during the single-limb stance phase of gait

Many children with cerebral palsy walk in a crouch gait that progressively worsens over time, decreasing walking efficiency and leading to joint degeneration. This study examined the effect of crouched postures on the capacity of muscles to extend the hip and knee joints and the joint flexions induced by gravity during the single-limb stance phase of gait. We first characterized representative mild, moderate, and severe crouch gait kinematics based on a large group of subjects with cerebral palsy (N=316). We then used a three-dimensional model of the musculoskeletal system and its associated equations of motion to determine the effect of these crouched gait postures on (1) the capacity of individual muscles to extend the hip and knee joints, which we defined as the angular accelerations of the joints, towards extension, that resulted from applying a 1N muscle force to the model, and (2) the angular acceleration of the joints induced by gravity. Our analysis showed that the capacities of almost all the major hip and knee extensors were markedly reduced in a crouched gait posture, with the exception of the hamstrings muscle group, whose extension capacity was maintained in a crouched posture. Crouch gait also increased the flexion accelerations induced by gravity at the hip and knee throughout single-limb stance. These findings help explain the increased energy requirements and progressive nature of crouch gait in patients with cerebral palsy.     

Muscular contributions to hip and knee extension during the single limb stance phase of normal gait: a framework for investigating the causes of crouch gait

Crouch gait, a troublesome movement abnormality among persons with cerebral palsy, is characterized by excessive flexion of the hips and knees during stance. Treatment of crouch gait is challenging, at present, because the factors that contribute to hip and knee extension during normal gait are not well understood, and because the potential of individual muscles to produce flexion or extension of the joints during stance is unknown. This study analyzed a three-dimensional, muscle-actuated dynamic simulation of walking to quantify the angular accelerations of the hip and knee induced by muscles during normal gait, and to rank the potential of the muscles to alter motions of these joints. Examination of the muscle actions during single limb stance showed that the gluteus maximus, vasti, and soleus make substantial contributions to hip and knee extension during normal gait. Per unit force, the gluteus maximus had greater potential than the vasti to accelerate the knee toward extension. These data suggest that weak hip extensors, knee extensors, or ankle plantar flexors may contribute to crouch gait, and strengthening these muscles--particularly gluteus maximus--may improve hip and knee extension. Abnormal forces generated by the iliopsoas or adductors may also contribute to crouch gait, as our analysis showed that these muscles have the potential to accelerate the hip and knee toward flexion. This work emphasizes the need to consider how muscular forces contribute to multijoint movements when attempting to identify the causes of abnormal gait.     

Quantitative assessment of gait determinants during single stance via a three-dimensional model--Part 1. Normal gait

In this two-part paper, a variety of three-dimensional, dynamical models are constructed for simulating the single support phases of normal and pathological human gait. A major objective of this work is to quantify the influence of individual gait determinants on the ground reaction forces generated during normal, level walking. To this end, Part 1 presents a three-dimensional, seven degree-of-freedom model incorporating five of the six fundamental determinants of gait. On the basis of crude muscle-force and/or joint-moment trajectories, body-segmental motions and ground reaction forces are synthesized open loop. Through a quantitative comparison with experimental gait data, the model's predictions are evaluated. Our simulation results suggest that pelvic list is not as dominant a dynamical determinant as either stance knee flexion-extension or foot and knee interaction. Transverse pelvic rotation, however, makes an important contribution by limiting the magnitude of the horizontal ground reaction prior to opposite heel-strike.     

Quantitative assessment of gait determinants during single stance via a three-dimensional model--Part 2. Pathological gait

A three-dimensional model for normal gait formulated in Part 1 is now altered to simulate the dynamics of pathological walking. Mechanisms fundamental to the production of a normal gait pattern are systematically removed, in order to assess contributions from individual gait determinants. Four separate pathological cases are studied: a model neglecting ankle plantarflexor activity; absence of stance knee flexion-extension and foot and knee interaction; both pelvic list and transverse pelvic rotation removed; and finally, a model with all major gait determinants missing. These are used collectively to show that stance knee flexion-extension and foot and knee interaction successively dominate lower-extremity dynamical response during the single support phase of normal gait. The hip abductor muscles, while effecting pelvic list, serve to stabilize this limb, rather than actively determine whole-body vertical acceleration. Mechanisms compensating for a loss in joint motion are also explored. Complete ankle loss may be successfully compensated with increased hip abductor muscle activity; the loss of both ankle and knee, however, demand unacceptable levels of vertical pelvic displacement.     

In vivo tibiofemoral contact analysis using 3D MRI-based knee models

This paper quantified the motion of the tibiofemoral contact points during in vivo weight bearing flexion using MRI- based 3D knee models and two orthogonal fluoroscopic images. The contact points on the medial and lateral tibial plateau were calculated by finding the centroid of the intersection of the tibial and femoral cartilage layers and by using the bony geometry alone. Our results indicate that the medial femoral condyle remains in the central portion of the tibial plateau and the lateral condyle translates posteriorly with increasing flexion. Using the bony contact model increased the total translation of the medial and lateral condyles by 250 and 55%, respectively, compared to the cartilage contact model. These results suggest that using the bony geometry alone may not accurately represent the articular surfaces of the knee. Articular cartilage geometry may have to be used to accurately quantify tibiofemoral contact.     

In vivo kinematics and articular surface congruency of total ankle arthroplasty during gait

Relatively high rates of loosening and implant failure have been reported after total ankle arthroplasty. Abnormal kinematics and incongruency of the articular surface may cause increased contact pressure and rotational torque applied to the implant, leading to loosening and implant failure. We measured in vivo kinematics of two-component total ankle arthroplasty (TNK ankle), and assessed congruency of the articular surface during the stance phase of gait. Eighteen ankles of 15 patients with a mean age of 75±6 years (mean±standard deviation) and follow-up of 44±38 months were enrolled. Lateral fluoroscopic images were taken during the stance phase of gait. 3D-2D model-image registration was performed using the fluoroscopic image and the implant models, and three-dimensional kinematics of the implant and incongruency of the articular surface were determined. The mean ranges of motion were 11.1±4.6°, 0.8±0.4°, and 2.6±1.5° for dorsi-/plantarflexion, inversion/eversion, and internal/external rotation, respectively. At least one type of incongruency of the articular surface occurred in eight of 18 ankles, including anterior hinging in one ankle, medial or lateral lift off in four ankles, and excessive axial rotation in five ankles. Among the four ankles in which lift off occurred during gait, only one ankle showed lift off in the static weightbearing radiograph. Our observations will provide useful data against which kinematics of other implant designs, such as three-component total ankle arthroplasty, can be compared. Our results also showed that evaluation of lift off in the standard weightbearing radiograph may not predict its occurrence during gait.     

Statistical shape model-based prediction of tibiofemoral cartilage

Abstract

Computed tomography is used more routinely to design patient-specific instrumentation for knee replacement surgery. Its moderate imaging cost and simplified segmentation reduce design costs compared with magnetic resonance (MR) imaging, but it cannot provide the necessary cartilage information. Our method based on statistical shape modelling proved to be successful in predicting tibiofemoral cartilage in leave-one-out experiments. The obtained accuracy of 0.54?mm for femur and 0.49?mm for tibia outperforms the average cartilage thickness distribution and reported inter-observer MR segmentation variability. These results suggest that shape modelling is able to predict tibiofemoral cartilage with sufficient accuracy to design patient-specific instrumentation.     

Muscle contributions to support and progression during single-limb stance in crouch gait

Pathological movement patterns like crouch gait are characterized by abnormal kinematics and muscle activations that alter how muscles support the body weight during walking. Individual muscles are often the target of interventions to improve crouch gait, yet the roles of individual muscles during crouch gait remain unknown. The goal of this study was to examine how muscles contribute to mass center accelerations and joint angular accelerations during single-limb stance in crouch gait, and compare these contributions to unimpaired gait. Subject-specific dynamic simulations were created for ten children who walked in a mild crouch gait and had no previous surgeries. The simulations were analyzed to determine the acceleration of the mass center and angular accelerations of the hip, knee, and ankle generated by individual muscles. The results of this analysis indicate that children walking in crouch gait have less passive skeletal support of body weight and utilize substantially higher muscle forces to walk than unimpaired individuals. Crouch gait relies on the same muscles as unimpaired gait to accelerate the mass center upward, including the soleus, vasti, gastrocnemius, gluteus medius, rectus femoris, and gluteus maximus. However, during crouch gait, these muscles are active throughout single-limb stance, in contrast to the modulation of muscle forces seen during single-limb stance in an unimpaired gait. Subjects walking in crouch gait rely more on proximal muscles, including the gluteus medius and hamstrings, to accelerate the mass center forward during single-limb stance than subjects with an unimpaired gait.     

A procedure to estimate normal and friction contact parameters in the stance phase of the human gait

A new approach to estimate normal and tangential contact parameters in the foot-ground contact during human gait was proposed. A correct estimation of the contact parameters would be very important in the resolution of predictive forward dynamic problems. The normal contact forces have been well estimated in the literature. But accurate estimation of tangential forces has not been reached yet. This work proposed a new procedure to accurately estimate friction forces. The approach has been based on the consideration of the modulus of the tangential force instead of its components. This modulus was introduced together with the modulus of the normal contact force and its two associated moments in an optimization algorithm to fit the contact forces provided by the model to the experimental data obtained with a force plate. An inverse dynamics problem was solved as a step previous to the optimization algorithm. The results showed that both the normal and tangential forces and the moments in the horizontal plane were in agreement with the experimental measurements. This work also analyzed the influence on the results of the friction law. The results obtained with the general friction law, which considered dry (static and dynamic) and viscous friction, were compared with results provided by simpler laws. The analysis of the components of the friction forces pointed out the importance of the Stribeck component in the resultant force instead of the viscous friction which played a minimal role. But for modelling the stick-slip transition, the implementation of a general friction law is necessary.     

In vivo patellar tracking and patellofemoral cartilage contacts during dynamic stair ascending

The knowledge of normal patellar tracking is essential for understanding the knee joint function and for diagnosis of patellar instabilities. This paper investigated the patellar tracking and patellofemoral joint contact locations during a stair ascending activity using a validated dual-fluoroscopic imaging system. The results showed that the patellar flexion angle decreased from 41.9° to 7.5° with knee extension during stair ascending. During first 80% of the activity, the patella shifted medially about 3.9mm and then slightly shifted laterally during the last 20% of the ascending activity. Anterior translation of 13mm of the patella was measured at the early 80% of the activity and the patella slightly moved posteriorly by about 2mm at the last 20% of the activity. The path of cartilage contact points was slightly lateral on the cartilage surfaces of patella and femur. On the patellar cartilage surface, the cartilage contact locations were about 2mm laterally from heel strike to 60% of the stair ascending activity and moved laterally and reached 5.3mm at full extension. However, the cartilage contact locations were relatively constant on the femoral cartilage surface (～5mm lateral). The patellar tracking pattern was consistent with the patellofemoral cartilage contact location pattern. These data could provide baseline knowledge for understanding of normal physiology of the patellofemoral joint and can be used as a reference for clinical evaluation of patellofemoral disorders.     

Biomechanical analysis of the stance phase during barefoot and shod running

This study investigated spatio-temporal variables, ground reaction forces and sagittal and frontal plane kinematics during the stance phase of nine trained subjects running barefoot and shod at three different velocities (3.5, 4.5, 5.5 m s(-1)). Differences between conditions were detected with the general linear method (factorial model). Barefoot running is characterized by a significantly larger external loading rate than the shod condition. The flatter foot placement at touchdown is prepared in free flight, implying an actively induced adaptation strategy. In the barefoot condition, plantar pressure measurements reveal a flatter foot placement to correlate with lower peak heel pressures. Therefore, it is assumed that runners adopt this different touchdown geometry in barefoot running in an attempt to limit the local pressure underneath the heel. A significantly higher leg stiffness during the stance phase was found for the barefoot condition. The sagittal kinematic adaptations between conditions were found in the same way for all subjects and at the three running velocities. However, large individual variations were observed between the runners for the rearfoot kinematics.     

Toe-out gait in patients with knee osteoarthritis partially transforms external knee adduction moment into flexion moment during early stance phase of gait: a tri-planar kinetic mechanism

Altered gait kinematics and kinetics are observed in patients with medial compartment knee osteoarthritis. Although various kinematic adaptations are proposed to be compensatory mechanisms that unload the knee, the nature of these mechanisms is presently unclear. We hypothesized that an increased toe-out angle during early stance phase of gait shifts load away from the knee medial compartment, quantified as the external adduction moment about the knee. Specifically, we hypothesized that by externally rotating the lower limb anatomy, primarily about the hip joint, toe-out gait alters the lengths of ground reaction force lever arms acting about the knee joint in the frontal and sagittal planes and transforms a portion of knee adduction moment into flexion moment. To test this hypothesis, gait data from 180 subjects diagnosed with medial compartment knee osteoarthritis were examined using two frames of reference. The first frame was attached to the tibia (reporting actual toe-out) and the second frame was attached to the laboratory (simulating no-toe-out). Four measures were compared within subjects in both frames of reference: the lengths of ground reaction force lever arms acting about the knee joint in the frontal and sagittal planes, and the adduction and flexion components of the external knee moment. The mean toe-out angle was 11.4 degrees (S.D. 7.8 degrees , range -2.2 degrees to 28.4 degrees ). Toe-out resulted in significant reductions in the frontal plane lever arm (-6.7%) and the adduction moment (-11.7%) in early stance phase when compared to the simulated no-toe-out values. These reductions were coincident with significant increases in the sagittal plane lever arm (+33.7%) and flexion moment (+25.0%). Peak adduction lever arm and moment were also reduced significantly in late stance phase (by -22.9% and -34.4%, respectively) without a corresponding increase in sagittal plane lever arm or flexion moment. These results indicate that toe-out gait in patients with medial compartment knee osteoarthritis transforms a portion of the adduction moment into flexion moment in early stance phase, suggesting that load is partially shifted away from the medial compartment to other structures.     

In vivo measurement of surface skin strain during human gait to improve the design of rehabilitation devices

Abstract

When designing any rehabilitation, sportswear or exoskeleton device the mechanical behaviour of the body segment must be known, specifically the skin, because an excessive tissue strain may lead to ulceration and bedsores. To date, it is not known if the kinematic variability between subjects have an effect on the skin strain field, and therefore, in the design and manufacturing of rehabilitation products, such as orthoses. Several studies have analysed the skin deformation during human motion, nevertheless, the comparison between the skin strain field in different subjects during normal or pathological gait has not been reported yet. This work presents a comparison of skin strain analysis for different gait patterns to study the differences between people and, specifically, if it is possible to standardize the orthotic design between subjects with the same gait disorder. Moreover, the areas with relatively minimum strain during the ankle-foot motion are compared to improve the design of structural parts of rehabilitation devices. In this case, a validated 3D digital image correlation system has been used for this purpose combined with strain ellipse theory. The results demonstrate variations in the skin strain field between subjects with the same pathology and similarities between subjects with normal gait. However, more studies and experiments are necessaries to validate this hypothesis and also to test it between different gait pathologies.     

The natural frequency of the foot-surface cushion during the stance phase of running

Researchers have reported on the stiffness of running in holistic terms, i.e. for the structures that are undergoing deformation as a whole rather than in terms of specific locations. This study aimed to estimate both the natural frequency and the viscous damping coefficient of the human foot-surface cushion, during the period between the heel strike and the mid-stance phase of running, using a purposely developed one degree-of-freedom inverted pendulum state space model of the leg. The model, which was validated via a comparison of measured and estimated ground reaction forces, incorporated a novel use of linearized and extended Kalman filter estimators. Investigation of the effect of variation of the natural frequency and/or the damping of the cushioning mechanism during running, using the said model, revealed the natural frequency of running on said foot-surface cushion, during the stance phase, to lie between 5 and 11 Hz. The "extended Kalman filter (EKF)" approach, that was used here for the first time to directly apply measured ground forces, may be widely applicable to the identification process of combined estimation of both unknown physiological state and mechanical characteristics of the environment in an inverse dynamic model.