An upper extremity inverse dynamics model for pediatric Lofstrand crutch-assisted gait

The objective of this study was to develop an instrumented Lofstrand crutch system, which quantifies three-dimensional (3-D) upper extremity (UE) kinematics and kinetics using an inverse dynamics model. The model describes the dynamics of the shoulders, elbows, wrists, and crutches and is compliant with the International Society of Biomechanics (ISB) recommended standards. A custom designed Lofstrand crutch system with four, six-degree-of-freedom force transducers was implemented with the inverse dynamics model to obtain triaxial UE joint reaction forces and moments. The crutch system was validated statically and dynamically for accuracy of computing joint reaction forces and moments during gait. The root mean square (RMS) error of the system ranged from 0.84 to 5.20%. The system was demonstrated in children with diplegic cerebral palsy (CP), incomplete spinal cord injury (SCI), and type I osteogenesis imperfecta (OI). The greatest joint reaction forces were observed in the posterior direction of the wrist, while shoulder flexion moments were the greatest joint reaction moments. The subject with CP showed the highest forces and the subject with SCI demonstrated the highest moments. Dynamic quantification may help to elucidate UE joint demands in regard to pain and pathology in long-term assistive device users.     

Kinematics of lower limbs during walking are emulated by springy walking model with a compliantly connected,off-centered curvy foot

The dynamics of the center of mass (CoM) during walking and running at various gait conditions are well described by the mechanics of a simple passive spring loaded inverted pendulum (SLIP). Due to its simplicity, however, the current form of the SLIP model is limited at providing any further information about multi-segmental lower limbs that generate oscillatory CoM behaviors and their corresponding ground reaction forces. Considering that the dynamics of the CoM are simply achieved by mass-spring mechanics, we wondered whether any of the multi-joint motions could be demonstrated by simple mechanics. In this study, we expand a SLIP model of human locomotion with an off-centered curvy foot connected to the leg by a springy segment that emulates the asymmetric kinematics and kinetics of the ankle joint. The passive dynamics of the proposed expansion of the SLIP model demonstrated the empirical data of ground reaction forces, center of mass trajectories, ankle joint kinematics and corresponding ankle joint torque at various gait speeds. From the mechanically simulated trajectories of the ankle joint and CoM, the motion of lower-limb segments, such as thigh and shank angles, could be estimated from inverse kinematics. The estimation of lower limb kinematics showed a qualitative match with empirical data of walking at various speeds. The representability of passive compliant mechanics for the kinetics of the CoM and ankle joint and lower limb joint kinematics implies that the coordination of multi-joint lower limbs during gait can be understood with a mechanical framework.     

Effect of arm swinging on lumbar spine and hip joint forces

During level walking, arm swing plays a key role in improving dynamic stability. In vivo investigations with a telemeterized vertebral body replacement showed that spinal loads can be affected by differences in arm positions during sitting and standing. However, little is known about how arm swing could influence the lumbar spine and hip joint forces and motions during walking. The present study aims to provide better understanding of the contribution of the upper limbs to human gait, investigating ranges of motion and joint reaction forces.A three-dimensional motion analysis was carried out via a motion capturing system on six healthy males and five patients with hip instrumented implant. Each subject performed walking with different arm swing amplitudes (small, normal, and large) and arm positions (bound to the body, and folded across the chest). The motion data were imported in a commercial musculoskeletal analysis software for kinematic and inverse dynamic investigation.The range of motion of the thorax with respect to the pelvis and of the pelvis with respect to the ground in the transversal plane were significantly associated with arm position and swing amplitude during gait. The hip external-internal rotation range of motion statistically varied only for non-dominant limb. Unlike hip joint reaction forces, predicted peak spinal loads at T12-L1 and L5-S1 showed significant differences at approximately the time of contralateral toe off and contralateral heel strike.Therefore, arm position and swing amplitude have a relevant effect on kinematic variables and spinal loads, but not on hip loads during walking.     

A musculoskeletal foot model for clinical gait analysis

Several full body musculoskeletal models have been developed for research applications and these models may potentially be developed into useful clinical tools to assess gait pathologies. Existing full-body musculoskeletal models treat the foot as a single segment and ignore the motions of the intrinsic joints of the foot. This assumption limits the use of such models in clinical cases with significant foot deformities. Therefore, a three-segment musculoskeletal model of the foot was developed to match the segmentation of a recently developed multi-segment kinematic foot model. All the muscles and ligaments of the foot spanning the modeled joints were included. Muscle pathways were adjusted with an optimization routine to minimize the difference between the muscle flexion–extension moment arms from the model and moment arms reported in literature. The model was driven by walking data from five normal pediatric subjects (aged 10.6±1.57 years) and muscle forces and activation levels required to produce joint motions were calculated using an inverse dynamic analysis approach. Due to the close proximity of markers on the foot, small marker placement error during motion data collection may lead to significant differences in musculoskeletal model outcomes. Therefore, an optimization routine was developed to enforce joint constraints, optimally scale each segment length and adjust marker positions. To evaluate the model outcomes, the muscle activation patterns during walking were compared with electromyography (EMG) activation patterns reported in the literature. Model-generated muscle activation patterns were observed to be similar to the EMG activation patterns.     

Upper limb loadings of gait with crutches

Long-term crutch users and patients with arthritis are particularly susceptible to upper limb joint degeneration during aided gait. The function of the walking aid for stability, support, and restraint/propulsion must be optimized with the upper limb loadings caused by the aids. Post-operative total hip replacement (THR) patients, tibial fracture, and paraplegic subjects using sticks and elbow crutches were analyzed in this study. Elbow and shoulder joint centers and aid orientations were monitored simultaneously in three dimensions and combined with aid forces to determine upper limb moment loadings. Three loading effects were observed: tendency for the aids to cause 1) the elbow to flex and shoulder to extend, 2) the elbow and shoulder to extend, and 3) the shoulder to abduct. Moment values of up to 0.10 Nm per body weight (BW) causing the shoulder to extend were measured, i.e., of similar magnitude to the moments at the hip in unaided gait. A modification of the elbow crutch, designed to improve medial-lateral stability, was unsuccessful in use due to wrist instability. This reinforced the requirement that crutch designs integrate the aid's function in gait with the ability of the upper limb joints to balance the applied loads.     

Direct comparison of measured and calculated total knee replacement force envelopes during walking in the presence of normal and abnormal gait patterns

Knee joint forces measured from instrumented implants provide important information for testing the validity of computational models that predict knee joint forces. The purpose of this study was to validate a parametric numerical model for predicting knee joint contact forces against measurements from four subjects with instrumented TKRs during the stance phase of gait. Model sensitivity to abnormal gait patterns was also investigated. The results demonstrated good agreement for three subjects with relatively normal gait patterns, where the difference between the mean measured and calculated forces ranged from 0.05 to 0.45 body weights, and the envelopes of measured and calculated forces (from three walking trials) overlapped. The fourth subject, who had a "quadriceps avoidance" external moment pattern, initially had little overlap between the measured and calculated force envelopes. When additional constraints were added, tailored to the subject's gait pattern, the model predictions improved to complete force envelope overlap. Coefficient of multiple determination analysis indicated that the shape of the measured and calculated force waveforms were similar for all subjects (adjusted coefficient of multiple correlation values between 0.88 and 0.92). The parametric model was accurate in predicting both the magnitude and waveform of the contact force, and the accuracy of model predictions was affected by deviations from normal gait patterns. Equally important, the envelope of forces generated by the range of solutions substantially overlapped with the corresponding measured envelope from multiple gait trials for a given subject, suggesting that the variable strategic processes of in vivo force generation are covered by the solution range of this parametric model.     

Effect of body weight support variation on muscle activities during robot assisted gait: a dynamic simulation study

Background and Objectives: While body weight support (BWS) intonation is vital during conventional gait training of neurologically challenged subjects, it is important to evaluate its effect during robot assisted gait training. In the present research we have studied the effect of BWS intonation on muscle activities during robotic gait training using dynamic simulations. Methods: Two dimensional (2-D) musculoskeletal model of human gait was developed conjointly with another 2-D model of a robotic orthosis capable of actuating hip, knee and ankle joints simultaneously. The musculoskeletal model consists of eight major muscle groups namely; soleus (SOL), gastrocnemius (GAS), tibialis anterior (TA), hamstrings (HAM), vasti (VAS), gluteus maximus (GLU), uniarticular hip flexors (iliopsoas, IP), and Rectus Femoris (RF). BWS was provided at levels of 0, 20, 40 and 60% during the simulations. In order to obtain a feasible set of muscle activities during subsequent gait cycles, an inverse dynamics algorithm along with a quadratic minimization algorithm was implemented. Results: The dynamic parameters of the robot assisted human gait such as joint angle trajectories, ground contact force (GCF), human limb joint torques and robot induced torques at different levels of BWS were derived. The patterns of muscle activities at variable BWS were derived and analysed. For most part of the gait cycle (GC) the muscle activation patterns are quite similar for all levels of BWS as is apparent from the mean of muscle activities for the complete GC. Conclusions: Effect of BWS variation during robot assisted gait on muscle activities was studied by developing dynamic simulation. It is expected that the proposed dynamic simulation approach will provide important inferences and information about the muscle function variations consequent upon a change in BWS during robot assisted gait. This information shall be quite important while investigating the influence of BWS intonation on neuromuscular parameters of interest during robotic gait training.     

A Dynamic Finite Element Analysis of Human Foot Complex in the Sagittal Plane during Level Walking

The objective of this study is to develop a computational framework for investigating the dynamic behavior and the internal loading conditions of the human foot complex during locomotion. A subject-specific dynamic finite element model in the sagittal plane was constructed based on anatomical structures segmented from medical CT scan images. Three-dimensional gait measurements were conducted to support and validate the model. Ankle joint forces and moment derived from gait measurements were used to drive the model. Explicit finite element simulations were conducted, covering the entire stance phase from heel-strike impact to toe-off. The predicted ground reaction forces, center of pressure, foot bone motions and plantar surface pressure showed reasonably good agreement with the gait measurement data over most of the stance phase. The prediction discrepancies can be explained by the assumptions and limitations of the model. Our analysis showed that a dynamic FE simulation can improve the prediction accuracy in the peak plantar pressures at some parts of the foot complex by 10%–33% compared to a quasi-static FE simulation. However, to simplify the costly explicit FE simulation, the proposed model is confined only to the sagittal plane and has a simplified representation of foot structure. The dynamic finite element foot model proposed in this study would provide a useful tool for future extension to a fully muscle-driven dynamic three-dimensional model with detailed representation of all major anatomical structures, in order to investigate the structural dynamics of the human foot musculoskeletal system during normal or even pathological functioning.     

Static and dynamic optimization solutions for gait are practically equivalent

The proposition that dynamic optimization provides better estimates of muscle forces during gait than static optimization is examined by comparing a dynamic solution with two static solutions. A 23-degree-of-freedom musculoskeletal model actuated by 54 Hill-type musculotendon units was used to simulate one cycle of normal gait. The dynamic problem was to find the muscle excitations which minimized metabolic energy per unit distance traveled, and which produced a repeatable gait cycle. In the dynamic problem, activation dynamics was described by a first-order differential equation. The joint moments predicted by the dynamic solution were used as input to the static problems. In each static problem, the problem was to find the muscle activations which minimized the sum of muscle activations squared, and which generated the joint moments input from the dynamic solution. In the first static problem, muscles were treated as ideal force generators; in the second, they were constrained by their force-length-velocity properties; and in both, activation dynamics was neglected. In terms of predicted muscle forces and joint contact forces, the dynamic and static solutions were remarkably similar. Also, activation dynamics and the force-length-velocity properties of muscle had little influence on the static solutions. Thus, for normal gait, if one can accurately solve the inverse dynamics problem and if one seeks only to estimate muscle forces, the use of dynamic optimization rather than static optimization is currently not justified. Scenarios in which the use of dynamic optimization is justified are suggested.     

Development and validation of a 3-D model to predict knee joint loading during dynamic movement

The purpose of this study was to develop a subject-specific 3-D model of the lower extremity to predict neuromuscular control effects on 3-D knee joint loading during movements that can potentially cause injury to the anterior cruciate ligament (ACL) in the knee. The simulation consisted of a forward dynamic 3-D musculoskeletal model of the lower extremity, scaled to represent a specific subject. Inputs of the model were the initial position and velocity of the skeletal elements, and the muscle stimulation patterns. Outputs of the model were movement and ground reaction forces, as well as resultant 3-D forces and moments acting across the knee joint. An optimization method was established to find muscle stimulation patterns that best reproduced the subject's movement and ground reaction forces during a sidestepping task. The optimized model produced movements and forces that were generally within one standard deviation of the measured subject data. Resultant knee joint loading variables extracted from the optimized model were comparable to those reported in the literature. The ability of the model to successfully predict the subject's response to altered initial conditions was quantified and found acceptable for use of the model to investigate the effect of altered neuromuscular control on knee joint loading during sidestepping. Monte Carlo simulations (N = 100,000) using randomly perturbed initial kinematic conditions, based on the subject's variability, resulted in peak anterior force, valgus torque and internal torque values of 378 N, 94 Nm and 71 Nm, respectively, large enough to cause ACL rupture. We conclude that the procedures described in this paper were successful in creating valid simulations of normal movement, and in simulating injuries that are caused by perturbed neuromuscular control.     

A three-dimensional whole-body model to predict human walking on level ground

Predictive simulation of human walking has great potential in clinical motion analysis and rehabilitation engineering assessment, but large computational cost and reliance on measurement data to provide initial guess have limited its wide use. We developed a computationally efficient model combining optimization and inverse dynamics to predict three-dimensional whole-body motions and forces during human walking without relying on measurement data. Using the model, we explored two different optimization objectives, mechanical energy expenditure and the time integral of normalized joint torque. Of the two criteria, the sum of the time integrals of the normalized joint torques produced a more realistic walking gait. The reason for this difference is that most of the mechanical energy expenditure is in the sagittal plane (based on measurement data) and this leads to difficulty in prediction in the other two planes. We conclude that mechanical energy may only account for part of the complex performance criteria driving human walking in three dimensions.

     

Muscle and external load contribution to knee joint contact loads during normal gait

Large knee adduction moments during gait have been implicated as a mechanical factor related to the progression and severity of tibiofemoral osteoarthritis and it has been proposed that these moments increase the load on the medial compartment of the knee joint. However, this mechanism cannot be validated without taking into account the internal forces and moments generated by the muscles and ligaments, which cannot be easily measured. Previous musculoskeletal models suggest that the medial compartment of the tibiofemoral joint bears the majority of the tibiofemoral load, with the lateral compartment unloaded at times during stance. Yet these models did not utilise explicitly measured muscle activation patterns and measurements from an instrumented prosthesis which do not portray lateral compartment unloading. This paper utilised an EMG-driven model to estimate muscle forces and knee joint contact forces during healthy gait. Results indicate that while the medial compartment does bear the majority of the load during stance, muscles provide sufficient stability to counter the tendency of the external adduction moment to unload the lateral compartment. This stability was predominantly provided by the quadriceps, hamstrings, and gastrocnemii muscles, although the contribution from the tensor fascia latae was also significant. Lateral compartment unloading was not predicted by the EMG-driven model, suggesting that muscle activity patterns provide useful input to estimate muscle and joint contact forces.     

Fibre operating lengths of human lower limb muscles during walking

Muscles actuate movement by generating forces. The forces generated by muscles are highly dependent on their fibre lengths, yet it is difficult to measure the lengths over which muscle fibres operate during movement. We combined experimental measurements of joint angles and muscle activation patterns during walking with a musculoskeletal model that captures the relationships between muscle fibre lengths, joint angles and muscle activations for muscles of the lower limb. We used this musculoskeletal model to produce a simulation of muscle-tendon dynamics during walking and calculated fibre operating lengths (i.e. the length of muscle fibres relative to their optimal fibre length) for 17 lower limb muscles. Our results indicate that when musculotendon compliance is low, the muscle fibre operating length is determined predominantly by the joint angles and muscle moment arms. If musculotendon compliance is high, muscle fibre operating length is more dependent on activation level and force-length-velocity effects. We found that muscles operate on multiple limbs of the force-length curve (i.e. ascending, plateau and descending limbs) during the gait cycle, but are active within a smaller portion of their total operating range.     

Influence of altered gait patterns on the hip joint contact forces

Children who exhibit gait deviations often present a range of bone deformities, particularly at the proximal femur. Altered gait may affect bone growth and lead to deformities by exerting abnormal stresses on the developing bones. The objective of this study was to calculate variations in the hip joint contact forces with different gait patterns. Muscle and hip joint contact forces of four children with different walking characteristics were calculated using an inverse dynamic analysis and a static optimisation algorithm. Kinematic and kinetic analyses were based on a generic musculoskeletal model scaled down to accommodate the dimensions of each child. Results showed that for all the children with altered gaits both the orientation and magnitude of the hip joint contact force deviated from normal. The child with the most severe gait deviations had hip joint contact forces 30% greater than normal, most likely due to the increase in muscle forces required to sustain his crouched stance. Determining how altered gait affects joint loading may help in planning treatment strategies to preserve correct loading on the bone from a young age.     

A probabilistic method to estimate gait kinetics in the absence of ground reaction force measurements

Human joint torques during gait are usually computed using inverse dynamics. This method requires a skeletal model, kinematics and measured ground reaction forces and moments (GRFM). Measuring GRFM is however only possible in a controlled environment. This paper introduces a probabilistic method based on probabilistic principal component analysis to estimate the joint torques for healthy gait without measured GRFM. A gait dataset of 23 subjects was obtained containing kinematics, measured GRFM and joint torques from inverse dynamics in order to obtain a probabilistic model. This model was then used to estimate the joint torques of other subjects without measured GRFM. Only kinematics, a skeletal model and timing of gait events are needed. Estimation only takes 0.28 ms per time instant. Using cross-validation, the resulting root mean square estimation errors for the lower-limb joint torques are found to be approximately 0.1 Nm/kg, which is 6–18% of the range of the ground truth joint torques. Estimated joint torque and GRFM errors are up to two times smaller than model-based state-of-the-art methods. Model-free artificial neural networks can achieve lower errors than our method, but are less repeatable, do not contain uncertainty information on the estimates and are difficult to use in situations which are not in the learning set. In contrast, our method performs well in a new situation where the walking speed is higher than in the learning dataset. The method can for example be used to estimate the kinetics during overground walking without force plates, during treadmill walking without (separate) force plates and during ambulatory measurements.     

Predicting muscle forces in gait from EMG signals and musculotendon kinematics

An EMG-driven muscle model for determining muscle force-time histories during gait is presented. The model, based on Hill's equation (1938), incorporates morphological data and accounts for changes in musculotendon length, velocity, and the level of muscle excitation for both concentric and eccentric contractions. Musculotendon kinematics were calculated using three-dimensional cinematography with a model of the musculoskeletal system. Muscle force-length-EMG relations were established from slow isokinetic calibrations. Walking muscle force-time histories were determined for two subjects. Joint moments calculated from the predicted muscle forces were compared with moments calculated using a linked segment, inverse dynamics approach. Moment curve correlations ranged from r = 0.72 to R = 0.97 and the root mean square (RMS) differences were from 10 to 20 Nm. Expressed as a relative RMS, the moment differences ranged from a low of 23% at the ankle to a high of 72% at the hip. No single reason for the differences between the two moment curves could be identified. Possible explanations discussed include the linear EMG-to-force assumption and how well the EMG-to-force calibration represented excitation for the whole muscle during gait, assumptions incorporated in the muscle modeling procedure, and errors inherent in validating joint moments predicted from the model to moments calculated using linked segment, inverse dynamics. The closeness with which the joint moment curves matched in the present study supports using the modeling approach proposed to determine muscle forces in gait.     

Pattern of anterior cruciate ligament force in normal walking

The goal of this study was to calculate and explain the pattern of anterior cruciate ligament (ACL) loading during normal level walking. Knee-ligament forces were obtained by a two-step procedure. First, a three-dimensional (3D) model of the whole body was used together with dynamic optimization theory to calculate body-segmental motions, ground reaction forces, and leg-muscle forces for one cycle of gait. Joint angles, ground reaction forces, and muscle forces obtained from the gait simulation were then input into a musculoskeletal model of the lower limb that incorporated a 3D model of the knee. The relative positions of the femur, tibia, and patella and the forces induced in the knee ligaments were found by solving a static equilibrium problem at each instant during the simulated gait cycle. The model simulation predicted that the ACL bears load throughout stance. Peak force in the ACL (303 N) occurred at the beginning of single-leg stance (i.e., contralateral toe off). The pattern of ACL force was explained by the shear forces acting at the knee. The balance of muscle forces, ground reaction forces, and joint contact forces applied to the leg determined the magnitude and direction of the total shear force acting at the knee. The ACL was loaded whenever the total shear force pointed anteriorly. In early stance, the anterior shear force from the patellar tendon dominated the total shear force applied to the leg, and so maximum force was transmitted to the ACL at this time. ACL force was small in late stance because the anterior shear forces supplied by the patellar tendon, gastrocnemius, and tibiofemoral contact were nearly balanced by the posterior component of the ground reaction.     

A practical solution to reduce soft tissue artifact error at the knee using adaptive kinematic constraints

Musculoskeletal modeling and simulations have vast potential in clinical and research fields, but face various challenges in representing the complexities of the human body. Soft tissue artifact from skin-mounted markers may lead to non-physiological representation of joint motions being used as inputs to models in simulations. To address this, we have developed adaptive joint constraints on five of the six degree of freedom of the knee joint based on in vivo tibiofemoral joint motions recorded during walking, hopping and cutting motions from subjects instrumented with intra-cortical pins inserted into their tibia and femur. The constraint boundaries vary as a function of knee flexion angle and were tested on four whole-body models including four to six knee degrees of freedom. A musculoskeletal model developed in OpenSim simulation software was constrained to these in vivo boundaries during level gait and inverse kinematics and dynamics were then resolved. Statistical parametric mapping indicated significant differences (p < 0.05) in kinematics between bone pin constrained and unconstrained model conditions, notably in knee translations, while hip and ankle flexion/extension angles were also affected, indicating the error at the knee propagates to surrounding joints. These changes to hip, knee, and ankle kinematics led to measurable changes in hip and knee transverse plane moments, and knee frontal plane moments and forces. Since knee flexion angle can be validly represented using skin mounted markers, our tool uses this reliable measure to guide the five other degrees of freedom at the knee and provide a more valid representation of the kinematics for these degrees of freedom.