Quadrupedal locomotion in squirrel monkeys (Cebidae: Saimiri sciureus): a cineradiographic study of limb kinematics and related substrate reaction forces

Quadrupedal locomotion of squirrel monkeys on small-diameter support was analyzed kinematically and kinetically to specify the timing between limb movements and substrate reaction forces. Limb kinematics was studied cineradiographically, and substrate reaction forces were synchronously recorded. Squirrel monkeys resemble most other quadrupedal primates in that they utilize a diagonal sequence/diagonal couplets gait when walking on small branches. This gait pattern and the ratio between limb length and trunk length influence limb kinematics. Proximal pivots of forelimbs and hindlimbs are on the same horizontal plane, thus giving both limbs the same functional length. However, the hindlimbs are anatomically longer than the forelimbs. Therefore, hindlimb joints must be more strongly flexed during the step cycle compared to the forelimb joints. Because the timing of ipsilateral limb movements prevents an increasing amount of forelimb retraction, the forelimb must be more protracted during the initial stance phase. At this posture, gravity acts with long moment arms at proximal forelimb joints. Squirrel monkeys support most of their weight on their hindlimbs. The timing of limb movements and substrate reaction forces shows that the shift of support to the hindlimbs is mainly done to reduce the compressive load on the forelimb. The hypothesis of the posterior weight shift as a dynamic strategy to reduce load on forelimbs, proposed by Reynolds ([1985]) Am. J. Phys. Anthropol. 67:335-349; [1985] Am. J. Phys. Anthropol. 67:351-362), is supported by the high correlation of ratios between forelimb and hindlimb peak vertical forces and the range of motion of shoulder joint and scapula in primates.     

Direct kinematic modeling of the upper limb during trunk-assisted reaching

The study proposes a rigid-body biomechanical model of the trunk and whole upper limb including scapula and the test of this model with a kinematic method using a six-dimensional (6-D) electromagnetic motion capture (mocap) device. Large unconstrained natural trunk-assisted reaching movements were recorded in 7 healthy subjects. The 3-D positions of anatomical landmarks were measured and then compared to their estimation given by the biomechanical chain fed with joint angles (the direct kinematics). Thus, the prediction errors was attributed to the different joints and to the different simplifications introduced in the model. Large (approx. 4 cm) end-point prediction errors at the level of the hand were reduced (to approx. 2 cm) if translations of the scapula were taken into account. As a whole, the 6-D mocap seems to give accurate results, except for prono-supination. The direct kinematic model could be used as a virtual mannequin for other applications, such as computer animation or clinical and ergonomical evaluations.     

The role of intersegmental dynamics during rapid limb oscillations

The interactive dynamic effects of muscular, inertial and gravitational moments on rapid, multi-segmented limb oscillations were studied. Using three-segment, rigid-body equations of motion, hip, knee and ankle intersegmental dynamics were calculated for the steady-state cycles of the paw-shake response in adult spinal cats. Hindlimb trajectories were filmed to obtain segmental kinematics, and myopotentials of flexors and extensors at each of the three joints were recorded synchronously with the ciné film. The segmental oscillations that emerged during the paw-shake response were a consequence of an interplay between active and passive musculotendinous forces, inertial forces, and gravity. During steady-state oscillations, the amplitudes of joint excursions, peak angular velocities, and peak angular accelerations increased monotonically and significantly in magnitude from the proximal joint (hip) to the most distal joint (ankle). In contrast to these kinematic relationships, the maximal values of net moments at the hip and knee were equal in magnitude, but of significantly lower magnitude than the large net moment at the ankle joint. At both the ankle and the knee, the flexor and extensor muscle moments were equal, but at the hip the magnitude of the peak flexor muscle moment was significantly greater than the extensor muscle moment. Muscle moments at the hip not only acted to counterbalance accelerations of the more distal segments, but also acted to maintain the postural orientation of the hindlimb. Large muscle moments at the knee functioned to counterbalance the large inertial moments generated by the large angular accelerations of the paw. At the ankle, the muscle moments dominated the generation of the paw accelerations. At the ankle and the knee, muscle moments controlled limb dynamics by slowing and reversing joint motions, and the active muscle forces contributing to ankle and knee moments were derived from lengthening of active musculotendinous units. In contrast to the more distal joints, the active muscles crossing the hip predominantly shortened as a result of the interplay among inertial forces and gravitational moments. The muscle function and kinetic data explain key features of the complex interactions that occur between central control mechanisms and multi-segmented, oscillating limb segments during the paw-shake response.     

Agreement of measured and calculated muscle activity during highly dynamic movements modelled with a spherical knee joint

The inclusion of muscle forces into the analysis of joint contact forces has improved their accuracy. But it has not been validated if such force and activity calculations are valid during highly dynamic multidirectional movements. The purpose of this study was to validate calculated muscle activation of a lower extremity model with a spherical knee joint for running, sprinting and 90°-cutting. Kinematics, kinetics and lower limb muscle activation of ten participants were investigated in a 3D motion capture setup including EMG. A lower extremity rigid body model was used to calculate the activation of these muscles with an inverse dynamics approach and a cubic cost function. Correlation coefficients were calculated to compare measured and calculated activation. The results showed good correlation of the modelled and calculated data with a few exceptions. The highest average correlations were found during walking (r = 0.81) and the lowest during cutting (r = 0.57). Tibialis anterior had the lowest average correlation (r = 0.33) over all movements while gastrocnemius medius had the highest correlation (r = 0.9). The implementation of a spherical knee joint increased the agreement between measured and modelled activation compared to studies using a hinge joint knee. Although some stabilizing muscles showed low correlations during dynamic movements, the investigated model calculates muscle activity sufficiently.     

A computer-based biomechanical analysis of the three-dimensional motion of cementless hip prostheses.

A computer-based mathematical technique was developed to measure and completely describe the migration and micromotion of a femoral hip prosthesis relative to the femur. This technique utilized the mechanics of rigid-body motion analysis and apparatus of seven linear displacement transducers to measure and describe the complete three-dimensional motion of the prosthesis during cyclic loading. Computer acquisition of the data and custom analysis software allowed one to calculate the magnitude and direction of the motion of any point of interest on the prostheses from information about the motion of two points on the device. The data were also used to replay the tests using a computer animation technique, which allowed a magnified view of the three-dimensional motion of the prosthesis. This paper describes the mathematical development of the rigid-body motion analysis, the experimental method and apparatus for data collection, the technique used to animate the motion, the sources of error and the effect of the assumptions (rigid bodies) on the results. Selected results of individual test runs of uncemented and cemented prostheses are presented to demonstrate the efficacy of the method. The combined effect of the vibration and electrical noise resulted in a resolution of the system of about 3-5 microns motion for each transducer. Deformation effects appear to contribute about 3-15 microns to the measurement error. This measurement and analysis technique is a very sensitive and powerful means of assessing the effects of different design parameters on the migration and micromotion of total joint prostheses and can be applied to any other case (knee, dental implant) where three-dimensional relative motion between two bodies is important.     

An improved inverse dynamics formulation for estimation of external and internal loads during human sagittal plane movements

Planar musculoskeletal models are common in the inverse dynamics analysis of human movements such as walking, running and jumping. The continued interest in such models is justified by their simplicity and computational efficiency. Related to a human planar model, a unified formulation for both the flying and support phases of the sagittal plane movements is developed. The actuation involves muscle forces in the lower limbs and the resultant muscle torques in the other body joints. The dynamic equations, introduced in absolute coordinates of the segments, are converted into useful compact forms using the projective technique. The solution to a determinate inverse dynamics problem allows for the explicit determination of the external reactions (presumed to vanish during the flying phases) and the resultant muscle torques in all the model joints. The indeterminate inverse dynamics problem is then focused on the assessment of muscle forces and joint reaction forces selectively in the supporting lower limb. Numerical results of the inverse dynamics simulation of sample sagittal plane movements are reported to illustrate the validity and effectiveness of the improved formulation.     

Critical analysis of musculoskeletal modelling complexity in multibody biomechanical models of the upper limb

The inverse dynamics technique applied to musculoskeletal models, and supported by optimisation techniques, is used extensively to estimate muscle and joint reaction forces. However, the solutions of the redundant muscle force sharing problem are sensitive to the detail and modelling assumptions of the models used. This study presents four alternative biomechanical models of the upper limb with different levels of discretisation of muscles by bundles and muscle paths, and their consequences on the estimation of the muscle and joint reaction forces. The muscle force sharing problem is solved for the motions of abduction and anterior flexion, acquired using video imaging, through the minimisation of an objective function describing muscle metabolic energy consumption. While looking for the optimal solution, not only the equations of motion are satisfied but also the stability of the glenohumeral and scapulothoracic joints is preserved. The results show that a lower level of muscle discretisation provides worse estimations regarding the muscle forces. Moreover, the poor discretisation of muscles relevant to the joint in analysis limits the applicability of the biomechanical model. In this study, the biomechanical model of the upper limb describing the infraspinatus by a single bundle could not solve the complete motion of anterior flexion. Despite the small differences in the magnitude of the forces predicted by the biomechanical models with more complex muscular systems, in general, there are no significant variations in the muscular activity of equivalent muscles.     

A 3-D dynamic model of human finger for studying free movements.

The purpose of this work is to develop a 3D inverse dynamic model of the human finger for estimating the muscular forces involved during free finger movements. A review of the existing 3D models of the fingers is presented, and an alternative one is proposed. The validity of the model has been proved by means of two simulations: free flexion-extension motion of all joints, and free metacarpophalangeal (MCP) adduction motion. The simulation shows the need for a dynamic model including inertial effects when studying fast movements and the relevance of modelling passive forces generated by the structures studying free movements, such as the force exerted by the muscles when they are stretched and the passive action of the ligaments over the MCP joint in order to reproduce the muscular force pattern during the simulation of the free MCP abduction-adduction movements.     

Skeletal kinematics of the midtarsal joint during walking: Midtarsal joint locking revisited

The kinematics of the human foot complex have been investigated to understand the weight bearing mechanism of the foot. This study aims to investigate midtarsal joint locking during walking by noninvasively measuring the movements of foot bones using a high-speed bi-planar fluoroscopic system. Eighteen healthy subjects volunteered for the study; the subjects underwent computed tomography imaging and bi-planar radiographs of the foot in order to measure the three-dimensional (3D) midtarsal joint kinematics using a 2D-to-3D registration method and anatomical coordinate system in each bone. The relative movements on bone surfaces were also calculated in the talonavicular and calcaneocuboid joints and quantified as surface relative velocity vectors on articular surfaces to understand the kinematic interactions in the midtarsal joint. The midtarsal joint performed a coupled motion in the early stance to pronate the foot to extreme pose in the range of motion during walking and maintained this pose during the mid-stance. In the terminal stance, the talonavicular joint performed plantar-flexion, inversion, and internal rotation while the calcaneocuboid joint performed mainly inversion. The midtarsal joint moved towards an extreme supinated pose, rather than a minimum motion in the terminal stance. The study provides a new perspective to understand the kinematics and kinetics of the movement of foot bones and so-called midtarsal joint locking, during walking. The midtarsal joint continuously moved towards extreme poses together with the activation of muscle forces, which would support the foot for more effective force transfer during push-off in the terminal stance.     

Incorporating the length-dependent passive-force generating muscle properties of the extrinsic finger muscles into a wrist and finger biomechanical musculoskeletal model

Dynamic movement trajectories of low mass systems have been shown to be predominantly influenced by passive viscoelastic joint forces and torques compared to momentum and inertia. The hand is comprised of 27 small mass segments. Because of the influence of the extrinsic finger muscles, the passive torques about each finger joint become a complex function dependent on the posture of multiple joints of the distal upper limb. However, biomechanical models implemented for the dynamic simulation of hand movements generally don’t extend proximally to include the wrist and distal upper limb. Thus, they cannot accurately represent these complex passive torques. The purpose of this short communication is to both describe a method to incorporate the length-dependent passive properties of the extrinsic index finger muscles into a biomechanical model of the upper limb and to demonstrate their influence on combined movement of the wrist and fingers. Leveraging a unique set of experimental data, that describes the net passive torque contributed by the extrinsic finger muscles about the metacarpophalangeal joint of the index finger as a function of both metacarpophalangeal and wrist postures, we simulated the length-dependent passive properties of the extrinsic finger muscles. Dynamic forward simulations demonstrate that a model including these properties passively exhibits coordinated movement between the wrist and finger joints, mimicking tenodesis, a behavior that is absent when the length-dependent properties are removed. This work emphasizes the importance of incorporating the length-dependent properties of the extrinsic finger muscles into biomechanical models to study healthy and impaired hand movements.     

Hand function: peripheral and central constraints on performance.

The hand is one of the most fascinating and sophisticated biological motor systems. The complex biomechanical and neural architecture of the hand poses challenging questions for understanding the control strategies that underlie the coordination of finger movements and forces required for a wide variety of behavioral tasks, ranging from multidigit grasping to the individuated movements of single digits. Hence, a number of experimental approaches, from studies of finger movement kinematics to the recording of electromyographic and cortical activities, have been used to extend our knowledge of neural control of the hand. Experimental evidence indicates that the simultaneous motion and force of the fingers are characterized by coordination patterns that reduce the number of independent degrees of freedom to be controlled. Peripheral and central constraints in the neuromuscular apparatus have been identified that may in part underlie these coordination patterns, simplifying the control of multi-digit grasping while placing certain limitations on individuation of finger movements. We review this evidence, with a particular emphasis on how these constraints extend through the neuromuscular system from the behavioral aspects of finger movements and forces to the control of the hand from the motor cortex.     

Kinematic models of lower limb joints for musculo-skeletal modelling and optimization in gait analysis

Kinematic models of lower limb joints have several potential applications in musculoskeletal modelling of the locomotion apparatus, including the reproduction of the natural joint motion. These models have recently revealed their value also for in vivo motion analysis experiments, where the soft-tissue artefact is a critical known problem. This arises at the interface between the skin markers and the underlying bone, and can be reduced by defining multibody kinematic models of the lower limb and by running optimization processes aimed at obtaining estimates of position and orientation of relevant bones. With respect to standard methods based on the separate optimization of each single body segment, this technique makes it also possible to respect joint kinematic constraints. Whereas the hip joint is traditionally assumed as a 3 degrees of freedom ball and socket articulation, many previous studies have proposed a number of different kinematic models for the knee and ankle joints. Some of these are rigid, while others have compliant elements. Some models have clear anatomical correspondences and include real joint constraints; other models are more kinematically oriented, these being mainly aimed at reproducing joint kinematics. This paper provides a critical review of the kinematic models reported in literature for the major lower limb joints and used for the reduction of soft-tissue artefact. Advantages and disadvantages of these models are discussed, considering their anatomical significance, accuracy of predictions, computational costs, feasibility of personalization, and other features. Their use in the optimization process is also addressed, both in normal and pathological subjects.     

Similarities in the Neural Control of the Shoulder and Elbow Joints Belie Their Structural Differences

Movement of the hand in three dimensional space is primarily controlled by the orientation of the shoulder and elbow complexes. Due to discrepancies in proprioceptive acuity, overlap in motor cortex representation and grossly different anatomies between these joints, we hypothesized that there would be differences in the accuracy of aimed movements between the two joints. Fifteen healthy young adults were tested under four conditions – shoulder motion with the elbow constrained and unconstrained, and elbow motion with the shoulder constrained and unconstrained. End point target locations for each joint were set to coincide with joint excursions of 10, 20 or 30 degrees of either the shoulder or elbow joint. Targets were presented in a virtual reality environment. For the constrained condition, there were no significant differences in angular errors between the two joints, suggesting that the central nervous system represents linked segment models of the limb in planning and controlling movements. For the unconstrained condition, although angle errors were higher, hand position errors remained the same as those of the constrained trials. These results support the idea that the CNS utilizes abundant degrees of freedom to compensate for the potentially different contributions to end-point errors introduced by each joint.     

Motion analysis of leg joints associated with escape turns of the cockroach,Periplaneta americana

Considerable information is now available on the neural organization of the escape system of the American cockroach. To relate these data to the behavior, we need detailed information on the movements made at the principle leg joints that produce the turn. We used motion analysis of high speed video records to acquire such information. Records from both free ranging and tethered animals were analyzed. 1. We analyzed individual joint movements using a tethered preparation. Stimuli from 4 different angles around the animal were used. For all wind angles, the femur-tibia (FT) joint on the mesothoracic leg that is ipsilateral to the wind source extended while the contralateral mesothoracic FT joint flexed. This moved both of these legs laterally toward the wind source. In freely moving animals the FT movements provide forces that turn the animal away from the wind source. 2. The ipsilateral mesothoracic coxa-femur (CF) joint extended for all wind angles. The contralateral mesothoracic CF joint extended in response to most winds from the rear, but switched to flexion in response to wind from the side and front. As a result of these joint movements, rear wind resulted in rearward movements of the contralateral mesothoracic leg, while side and front wind resulted in more forward movements of that leg. 3. The CF and FT joints for both ipsilateral and contralateral metathoracic legs extended to wind from the rear and switched to flexion as the wind was placed at more anterior positions around the animal. In freely moving animals, extension of these joints would push the animal forward. Flexion would pull the animal backward. 4. Several of the joints showed correlations between rate of movement and initial joint angle. That is, joints that were already flexed at the onset of stimulation tended to move at a faster rate to a final position than joints that started at a more extended position. 5. Metathoracic FT and CF joints showed a high degree of positive correlation during the escape movements. Indeed, many curves showing movement of metathoracic FT and CF joints with time were virtually identical.     

Effect of a pedicle-screw-based motion preservation system on lumbar spine biomechanics: A probabilistic finite element study with subsequent sensitivity analysis

Pedicle-screw-based motion preservation systems are often used to support a slightly degenerated disc. Such implants are intended to reduce intradiscal pressure and facet joints forces, while having a minimal effect on the motion patterns.In a probabilistic finite element study with subsequent sensitivity analysis, the effects of 10 input parameters, such as elastic modulus and diameter of the elastic rod, distraction of the segment, level of bridged segments, etc. on the output parameters intervertebral rotations, intradiscal pressures, and facet joint forces were determined. A validated finite element model of the lumbar spine was employed. Probabilistic studies were performed for seven loading cases: upright standing, flexion, extension, left and right lateral bending and left and right axial rotation.The simulations show that intervertebral rotation angles, intradiscal pressures and facet joint forces are in most cases reduced by a motion preservation system. The influence on intradiscal pressure is small, except in extension. For many input parameter combinations, the values for intervertebral rotations and facet joint forces are very low, which indicates that the implant is too stiff in these cases. The output parameters are affected most by the following input parameters: loading case, elastic modulus and diameter of the elastic rod, distraction of the segment, and angular rigidity of the connection between screws and rod.The designated functions of a motion preservation system can best be achieved when the longitudinal rod has a low stiffness, and when the connection between rod and pedicle screws is rigid.     

Peculiarities of Activation of Muscles of the Shoulder Belt in Voluntary Two-Joint Movements of the Upper Limb

We studied coordination of central motor commands (СMCs) coming to muscles of the shoulder and shoulder belt in the course of single-joint and two-joint movements including flexion and extension of the elbow and shoulder joints. Characteristics of rectified and averaged EMGs recorded from a few muscles of the upper limb were considered correlates of the CMC parameters. Special attention was paid to coordination of CMCs coming to two-joint muscles that are able to function as common flexors (m. biceps brachii, caput breve, BBcb) and common extensors (m. triceps brachii, caput longum, TBcl) of the elbow and shoulder joints. Upper limb movements used in the tests included planar shifts of the arm from one spatial point to another resulting from either simultaneous changes in the angles of the shoulder and elbow joints or isolated sequential (two-stage) changes in these joint angles. As was found, shoulder muscles providing movements of the elbow with changes in the angle of the elbow joint, i.e., BBcb and TBcl, were also intensely involved in the performance of single-joint movements in the shoulder joint. The CMCs coming to two-joint muscles in the course of two-joint movements appeared, in the first approximation, as sums of the commands received by these muscles in the course of corresponding single-joint movements in the elbow and shoulder joints. Therefore, if we interpret the isolated forearm movement performed due to a change in the angle of the elbow joint as the main motor event, while the shoulder movement is considered the accessory one, we can conclude that realization of a two-joint movement of the upper-limb distal part is based on superposition of CMCs related to basic movements (main and accessory). Neirofiziologiya/Neurophysiology, Vol. 41, No. 1, pp. 48–56, January–February, 2009.     

An explicit expression for the moment in multibody systems.

The equations of motion are formulated for a set of interacting rigid bodies. An explicit formulation could be derived for the total moment of the forces on such a system. This is applied to (1) the intersegmental moment in a multisegment rigid-body model, and (2) the interpretation of the 'centre of pressure' as defined in the studies of human posture.     

Joint moments and contact forces in the foot during walking

The net force and moment of a joint have been widely used to understand joint disease in the foot. Meanwhile, it does not reflect the physiological forces on muscles and contact surfaces. The objective of the study is to estimate active moments by muscles, passive moments by connective tissues and joint contact forces in the foot joints during walking. Joint kinematics and external forces of ten healthy subjects (all males, 24.7 ± 1.2 years) were acquired during walking. The data were entered into the five-segment musculoskeletal foot model to calculate muscle forces and joint contact forces of the foot joints using an inverse dynamics-based optimization. Joint reaction forces and active, passive and net moments of each joint were calculated from muscle and ligament forces. The maximum joint reaction forces were 8.72, 4.31, 2.65, and 3.41 body weight (BW) for the ankle, Chopart’s, Lisfranc and metatarsophalangeal joints, respectively. Active and passive moments along with net moments were also obtained. The maximum net moments were 8.6, 8.4, 5.4 and 0.8%BW∙HT, respectively. While the trend of net moment was very similar between the four joints, the magnitudes and directions of the active and passive moments varied between joints. The active and passive moments during walking could reveal the roles of muscles and ligaments in each of the foot joints, which was not obvious in the net moment. This method may help narrow down the source of joint problems if applied to clinical studies.     

Development and validation of a multi-body model of the canine stifle joint

Multi-body musculoskeletal models that can be used concurrently to predict joint contact pressures and muscle forces would be extremely valuable in studying the mechanics of joint injury. The purpose of this study was to develop an anatomically correct canine stifle joint model and validate it against experimental data. A cadaver pelvic limb from one adult dog was used in this study. The femoral head was subjected to axial motion in a mechanical tester. Kinematic and force data were used to validate the computational model. The maximum RMS error between the predicted and measured kinematics during the complete testing cycle was 11.9 mm translational motion between the tibia and the femur and 4.3° rotation between patella and femur. This model is the first step in the development of a musculoskeletal model of the hind limb with anatomically correct joints to study cartilage loading under dynamic conditions.