Spasticity is an important complication after stroke, especially in the anti-gravity muscles, i.e. lower limb extensors. However the contribution of hyperexcitable muscle spindle reflex loops to gait ...impairments after stroke is often disputed. In this study a neuro-musculoskeletal model was developed to investigate the contribution of an increased length and velocity feedback and altered reflex modulation patterns to hemiparetic gait deficits.
A musculoskeletal model was extended with a muscle spindle model providing real-time length and velocity feedback of gastrocnemius, soleus, vasti and rectus femoris during a forward dynamic simulation (neural control model). By using a healthy subject's base muscle excitations, in combination with increased feedback gains and altered reflex modulation patterns, the effect on kinematics was simulated. A foot-ground contact model was added to account for the interaction effect between the changed kinematics and the ground. The qualitative effect i.e. the directional effect and the specific gait phases where the effect is present, on the joint kinematics was then compared with hemiparetic gait deviations reported in the literature.
Our results show that increased feedback in combination with altered reflex modulation patterns of soleus, vasti and rectus femoris muscle can contribute to excessive ankle plantarflexion/inadequate dorsiflexion, knee hyperextension/inadequate flexion and increased hip extension/inadequate flexion during dedicated gait cycle phases. Increased feedback of gastrocnemius can also contribute to excessive plantarflexion/inadequate dorsiflexion, however in combination with excessive knee and hip flexion. Increased length/velocity feedback can therefore contribute to two types of gait deviations, which are both in accordance with previously reported gait deviations in hemiparetic patients. Furthermore altered modulation patterns, in particular the reduced suppression of the muscle spindle feedback during swing, can contribute largely to an increased plantarflexion and knee extension during the swing phase and consequently to hampered toe clearance.
Our results support the idea that hyperexcitability of length and velocity feedback pathways, especially in combination with altered reflex modulation patterns, can contribute to deviations in hemiparetic gait. Surprisingly, our results showed only subtle temporal differences between length and velocity feedback. Therefore, we cannot attribute the effects seen in kinematics to one specific type of feedback.
Abstract The present study aimed to re-examine the influence of the isometric plantarflexors contraction on the Achilles tendon moment arm (ATMA) and the factors influencing the ATMA in ...three-dimensions. A series of coronal magnetic resonance images of the right ankle were recorded at foot positions of 10° of dorsiflexion, neutral position, and 10° of plantarflexion for the rest condition and the plantarflexors contraction condition at 30% maximal voluntary effort. The shortest distance between the talocrural joint axis and the line of action of the Achilles tendon force projected to the orthogonal plane of the talocrural joint axis was determined as the ATMA. The ATMA determined in the contraction condition was significantly greater by 8 mm than that determined in the rest condition. The talocrural joint axis was displaced anteriorly by 3 mm and distally by 2 mm due to the muscle contraction. As the same time, the line of action of the Achilles tendon force was displaced posteriorly by 5 mm and medially by 2 mm. These linear displacements of the talocrural joint axis and the line of action of the Achilles tendon force accounted for the difference in the ATMAs between the two conditions by 35.9 and 62.4%, respectively. These angular displacements accounted for the total of 0.4% increase in the ATMA. These results confirm the previous findings reported in two-dimensional studies and found that the linear displacement of the line of action of the Achilles tendon force is the primary source of the contraction-induced increase in the ATMA.
Mechanical and metabolic energy conservation is considered to be a defining characteristic in many common motor tasks. During human gait, the storage and return of elastic energy in compliant ...structures is an important energy saving mechanism that may reduce the necessary muscle fiber work and be an important determinant of the preferred gait mode (i.e., walk or run) at a given speed. In the present study, the mechanical work done by individual muscle fibers and series-elastic elements (SEE) was quantified using a musculoskeletal model and forward dynamical simulations that emulated a group of young healthy adults walking and running above and below the preferred walk-run transition speed (PTS), and potential advantages associated with the muscle fiber-SEE interactions during these gait modes at each speed were assessed. The simulations revealed that: (1) running below the PTS required more muscle fiber work than walking, and inversely, walking above the PTS required more muscle fiber work than running, and (2) SEE utilization in running was greater above than below the PTS. These results support previous suggestions that muscle mechanical energy expenditure is an important determinant for the preferred gait mode at a given speed.
Abstract Background The manual wheelchair user population experiences a high prevalence of upper-limb injuries, which are related to a high load on the shoulder joint during activities of daily ...living, such as handrim wheelchair propulsion. An alternative mode of propulsion is handcycling, where lower external forces are suggested to be applied to reach the same power output as in handrim wheelchair propulsion. This study aimed to quantify glenohumeral contact forces and muscle forces during handcycling and compare them to previous results of handrim wheelchair propulsion. Methods Ten able-bodied men propelled the handbike on a treadmill at two inclines (1% and 4% with a velocity of 1.66 m/s) and two speed conditions (1.39 and 1.94 m/s with fixed power output). Three-dimensional kinematics and kinetics were obtained and used as input for a musculoskeletal model of the arm and shoulder. Output variables were glenohumeral contact forces and forces of important shoulder muscles. Findings The highest mean and peak glenohumeral contact forces occurred at 4% incline (420 N, 890 N respectively). The scapular part of the deltoideus, the triceps and the trapezius produced the highest force. Interpretation Due to the circular movement and the continuous force application during handcycling, the glenohumeral contact forces, as well as the muscle forces were clearly lower compared to the results in the existing literature on wheelchair propulsion. These findings prove the assumption that handcycling is mechanically less straining than handrim wheelchair propulsion, which may help preventing overuse to the shoulder complex.
Abstract The rat is of increasing importance for experimental studies on fracture healing. The healing outcome of long bone fractures is strongly influenced by mechanical factors, such as the ...interfragmentary movement. This movement depends on the stability of the fracture fixation and the musculoskeletal loads. However, little is known about these loads in rats. The musculoskeletal loads during gait were estimated using an inverse-dynamic musculoskeletal model of the right hindlimb of the rat. This model was based on a micro-CT scan of the lower extremities and an anatomical study using 15 rat cadavers. Kinematics were reconstructed from X-ray movies, taken simultaneously from two perpendicular directions during a gait cycle. The ground reaction forces were taken from the literature. The muscle forces were calculated using an optimization procedure. The internal forces and moments varied over the gait cycle and along the femoral axis. The greatest internal force (up to 7 times bodyweight) acted in the longitudinal direction. The greatest internal moment (up to 13.8 bodyweight times millimeter) acted in the sagittal plane of the femur. The validity of the model was corroborated by comparing the estimated strains caused by the calculated loads on the surface of the femoral mid-shaft with those from the literature. Knowledge of the internal loads in the femur of the rat allows adjustment of the biomechanical properties of fixation devices in fracture healing studies to the desired interfragmentary movement.
An accurate, dynamic, functional model of the skull that can be used to predict muscle forces, bite forces, and joint reaction forces would have many uses across a broad range of disciplines. One ...major issue however with musculoskeletal analyses is that of muscle activation pattern indeterminacy. A very large number of possible muscle force combinations will satisfy a particular functional task. This makes predicting physiological muscle recruitment patterns difficult. Here we describe in detail the process of development of a complex multibody computer model of a primate skull (Macaca fascicularis), that aims to predict muscle recruitment patterns during biting. Using optimisation criteria based on minimisation of muscle stress we predict working to balancing side muscle force ratios, peak bite forces, and joint reaction forces during unilateral biting. Validation of such models is problematic; however we have shown comparable working to balancing muscle activity and TMJ reaction ratios during biting to those observed in vivo and that peak predicted bite forces compare well to published experimental data. To our knowledge the complexity of the musculoskeletal model is greater than any previously reported for a primate. This complexity, when compared to more simple representations provides more nuanced insights into the functioning of masticatory muscles. Thus, we have shown muscle activity to vary throughout individual muscle groups, which enables them to function optimally during specific masticatory tasks. This model will be utilised in future studies into the functioning of the masticatory apparatus.
► We describe a new musculoskeletal computer model of a primate (Macaque) skull. ► It predicts bite force, TMJ forces and muscle recruitment patterns during biting. ► The results compare very well with experimental in vivo EMG and bite force data. ► They reveal the complexity of muscle recruitment even during simple biting.
Accurate muscle geometry (muscle length and moment arm) is required to estimate muscle function when using musculoskeletal modelling. In shoulder, muscles are often modelled as a collection of ...independent line segments, leading to non-physiological muscles trajectory, especially for the rotator cuff muscles. To prevent this, a surface mesh model was developed and validated against 7 MRI positions in one participant. Mean moment arm errors was 11.4% for the line vs. 8.8% for the mesh model. While the model with independent lines led to some non-physiological trajectories, the mesh model gave lower misestimations of muscle lengths and moment arms.
The purpose of this study was to determine the influence of lumbar spine extension and erector spinae muscle activation on vertical jump height during maximal squat jumping. Eight male athletes ...performed maximal squat jumps. Electromyograms of the erector spinae were recorded during these jumps. A simulation model of the musculoskeletal system was used to simulate maximal squat jumping with and without spine extension. The effect on vertical jump height of changing erector spinae strength was also tested through the simulated jumps. Concerning the participant jumps, the kinematics indicated a spine extension and erector spinae activation. Concerning the simulated jumps, vertical jump height was about 5.4 cm lower during squat jump without trunk extension compared to squat jump
.
These results were explained by greater total muscle work during squat jump, more especially by the erector spinae work (+119.5 J). The erector spinae may contribute to spine extension during maximal squat jumping. The simulated jumps confirmed this hypothesis showing that vertical jumping was decreased if this muscle was not taken into consideration in the model. Therefore it is concluded that the erector spinae should be considered as a trunk extensor, which enables to enhance total muscle work and consequently vertical jump height.
A novel open-source biomechanical model of the index finger with an electromyography (EMG)-constrained static optimization solution method are developed with the goal of improving co-contraction ...estimates and providing means to assess tendon tension distribution through the finger. The Intrinsic model has four degrees of freedom and seven muscles (with a 14 component extensor mechanism). A novel plugin developed for the OpenSim modelling software applied the EMG-constrained static optimization solution method. Ten participants performed static pressing in three finger postures and five dynamic free motion tasks. Index finger 3D kinematics, force (5, 15, 30 N), and EMG (4 extrinsic muscles and first dorsal interosseous) were used in the analysis. The Intrinsic model predicted co-contraction increased by 29% during static pressing over the existing model. Further, tendon tension distribution patterns and forces, known to be essential to produce finger action, were determined by the model across all postures. The Intrinsic model and custom solution method improved co-contraction estimates to facilitate force propagation through the finger. These tools improve our interpretation of loads in the finger to develop better rehabilitation and workplace injury risk reduction strategies.
•Multiple balance strategies could be identified by a single optimization criterion.•Mixed ankle/hip strategies with straight knees were found when the toes were fixed.•Without fixing the toes, the ...hop strategy was identified under severe perturbation.•Simulated movements mimicked spontaneous recruiting or suppressing DOF in the CNS.•DOF number in balance depended on perturbation intensity and movement constraints.
In human balance recovery, different strategies have been proposed with generally overlooked knee motions but extensive focus on the ankle, hip, and step strategies. It is not well understood whether maintaining balance is regulated at the lower “muscular–articular” level of coordinating segment joints or at a higher level of controlling whole body dynamics. Whether balance control is to minimize joint degrees of freedom (DOF) or utilize all the available DOF also remains unclear. This study aimed to use a realistic musculoskeletal human model to identify multiple balance recovery strategies with a single optimization criterion. Movements were driven by neural excitations (which activated muscle force generation) and were assumed to be symmetric. Balance recoveries were simulated with forward-inclined straight body postures as the initial conditions. When the position of the toes was fixed, balance was regained with virtually straight knees and mixed ankle/hip strategies. Under a severely perturbed condition, use of the forward hop strategy after releasing the fixed-toes constraint indicated spontaneous recruitment or suppression of DOF, which mimicked functions of optimally computed CNS commands in humans. The results also indicated that increase/decrease in the number of DOF depends on the imposed perturbation intensity and movement constraints.