In vivo determination of the anatomical axes of the ankle joint complex: An optimization approach

This study investigates the feasibility of a subject-specific three-dimensional model of the ankle joint complex for kinematic and dynamic analysis of movement. The ankle joint complex was modelled as a three-segment system, connected by two ideal highe joints: the talocrural and the subtalar joint. A mathematical formulation was developed to express the three-dimensional translation and rotation between the foot and shank segments as a function of the two joint angles, and 12 model parameters describing the locations of the joint axes. An optimization method was used to fit the model parameters to three-dimensional kinematic data of foot and shank markers, obtained during test movements throughout the entire physiological range of motion of the ankle joint. The movement of the talus segment, which cannot be measured non-invasively, is not necessary for the analysis.

This optimization method was used to determine the position and orientation of the joint axes in 14 normal subjects. After optimization, the discrepancy between the best fitting model and actual marker kinematics was between 1 and 3 mm for all subjects. The predicted inclination of the subtalar joint axis from the horizontal plane was 37.4±2.7°, and the medial deviation was 18.0±16.2°. The lateral side of the talucrural axis was directed slightly posteriorly (6.8±8.1°), and inclined downward by 7.0±5.4°. These results are similar to previously reported typical results from anatomical, in vitro, studies. Reproducibility was evaluated by repeated testing of one subject, which resulted in variations of about one-fifth of the standard deviation within the group, the inclination of the subtalar joint axis was significantly correlated to the arch height and a radiographic ‘tarsal index’. It is concluded that this optimization method provides the opportunity to incorporate inter-individual anatomical differences into kinematic and dynamic analysis of the ankle joint complex. This allows a more functional interpretation of kinematic data, and more realistic estimates of internal forces.     

Minimizing errors during in vitro testing of multisegmental spine specimens: considerations for component selection and kinematic measurement

Apparatus-induced artifacts may invalidate standard spine testing protocols. Kinematic measurements may be compromised by the configuration of motion capture equipment. This study has determined: (1) the influence of machine design (component friction) on in vitro spinal kinetics; (2) the sensitivity of kinematic measurements to variations in the placement of motion capture markers. A spinal loading simulator has been developed to dynamically apply pure bending moments (three axes) with or without a simultaneous compressive preload. Two linear slider types with different friction coefficients, one with caged ball bearings and one with high-precision roller bearings on rails, were mounted and specimen response compared in sequential tests. Three different optoelectronic marker cluster configurations were mounted on the specimen and motion data was captured simultaneously from all clusters during testing. A polymer tube with a uniform bending stiffness approximately equivalent to a polysegmental lumbar spine specimen was selected to allow reproducible behavior over multiple tests. The selection of sliders for linear degrees of freedom had a marked influence on parasitic shear forces. Higher shear forces were recorded with the caged-bearing design than with the high-precision rollers and consequently a higher moment was required to achieve a given rotation. Kinematic accuracy varied with each marker configuration, but in general higher accuracy was achieved with larger marker spacings and situations where markers moved predominantly parallel to the camera's imaging plane. Relatively common alternatives in the mechanical components used in an apparatus for in vitro spine testing can have a significant influence on the measured kinematic and kinetics. Low-magnitude parasitic shear forces due to friction in sliders induces a linearly increasing moment along the length of the specimen, precluding the ideal of pure moment application. This effect is compounded in polysegmental specimens. Kinematic measurements are highly sensitive to marker design and placement, despite equivalent absolute precision of individual marker measurements, however marker configurations can be designed to minimize errors related to spatial distribution and system bias.     

Predictive modelling of cervical disc implant wear

This study presents a chain of simulations aimed at estimating the wear in a cervical disc implant and providing insight into the in vivo biomechanical performance of the implant. The simulation chain can start with determining a representative maximum range of motion (ROM) of a person's head. The ROM is used as motion input to a kinematic simulation of the cervical spine containing a disc implant. The cervical spine geometry is obtained from computed tomography (CT) scans and converted to STL format using reverse engineering software. The time histories of the loads imposed by the adjacent vertebrae on the implant, as well as the vertebral relative rotations can be extracted from the kinematic simulation. Alternatively, force and motion profiles prescribed by wear test protocols (e.g. ISO 18192-1 and ASTM F2423-05) can be used. The force and motion profiles are applied as boundary conditions to a non-linear finite element model (FEM) of the implant to determine the time-varying contact stress and slip velocity distributions at the interface between the two halves of the implant. The stresses and slip velocities are used in a linear wear model to estimate the wear rate distribution at the FEM's nodal points where contact occurs. Reverse engineering software is used to triangulate the contact surface so that the total wear volume can be calculated. The simulation chain's predicted wear rate shows good agreement with in vitro results in the literature. The simulation chain is thereby demonstrated to be suitable for comparative pre-experimental studies of spinal implant designs.     

A computational model for dynamic analysis of the human gait

Biomechanical models are important tools in the study of human motion. This work proposes a computational model to analyse the dynamics of lower limb motion using a kinematic chain to represent the body segments and rotational joints linked by viscoelastic elements. The model uses anthropometric parameters, ground reaction forces and joint Cardan angles from subjects to analyse lower limb motion during the gait. The model allows evaluating these data in each body plane. Six healthy subjects walked on a treadmill to record the kinematic and kinetic data. In addition, anthropometric parameters were recorded to construct the model. The viscoelastic parameter values were fitted for the model joints (hip, knee and ankle). The proposed model demonstrated that manipulating the viscoelastic parameters between the body segments could fit the amplitudes and frequencies of motion. The data collected in this work have viscoelastic parameter values that follow a normal distribution, indicating that these values are directly related to the gait pattern. To validate the model, we used the values of the joint angles to perform a comparison between the model results and previously published data. The model results show a same pattern and range of values found in the literature for the human gait motion.     

Animation of in vitro biomechanical tests.

Interdisciplinary communication of three-dimensional kinematic data arising from in vitro biomechanical tests is challenging. Complex kinematic representations such as the helical axes of motion (HAM) add to the challenge. The difficulty increases further when other quantities (i.e. load or tissue strain data) are combined with the kinematic data. The objectives of this study were to develop a method to graphically replay and animate in vitro biomechanical tests including HAM data. This will allow intuitive interpretation of kinematic and other data independent of the viewer's area of expertise. The value of this method was verified with a biomechanical test investigating load-sharing of the cervical spine. Three 3.0 mm aluminium spheres were glued to each of the two vertebrae from a C2-3 segment of a human cervical spine. Before the biomechanical tests, CT scans were made of the specimen (slice thickness=1.0 mm and slice spacing=1.5 mm). The specimens were subjected to right axial torsion moments (2.0 Nm). Strain rosettes mounted to the anterior surface of the C3 vertebral body and bilaterally beneath the facet joints on C3 were used to estimate the force flow through the specimen. The locations of the aluminium spheres were digitised using a space pointer and the motion analysis system. Kinematics were measured using an optoelectronic motion analysis system. HAMs were calculated to describe the specimen kinematics. The digitised aluminium sphere locations were used to match the CT and biomechanical test data (RMS errors between the CT and experimental points were less than 1.0 mm). The biomechanical tests were "replayed" by animating reconstructed CT models in accordance with the recorded experimental kinematics, using custom software. The animated test replays allowed intuitive analysis of the kinematic data in relation to the strain data. This technique improves the ability of experts from disparate backgrounds to interpret and discuss this type of biomechanical data.     

Kinematical models to reduce the effect of skin artifacts on marker-based human motion estimation

The estimation of the skeletal motion obtained from marker-based motion capture systems is known to be affected by significant bias caused by skin movement artifacts, which affects joint center and rotation axis estimation. Among different techniques proposed in the literature, that based on rigid body model, still the most used by commercial motion capture systems, can smooth only part of the above effects without eliminating their main components. In order to sensibly improve the accuracy of the motion estimation, a novel technique, named local motion estimation (LME), is proposed. This rests on a recently described approach that, using virtual humans and extended Kalman filters, estimates the kinematical variables directly from 2D measurements without requiring the 3D marker reconstruction. In this paper, we show how such method can be extended to include the computation of the local marker displacement due to skin artifacts. The 3D marker coordinates, expressed in the corresponding local reference coordinate frames, are inserted into the state vector of the filter and their dynamics is automatically estimated, with adequate accuracy, without assuming any particular deformation function. Simulated experiments of lower limb motion, involving systematic mislocations (5, 10, 20 mm) and random errors of the marker coordinates and joint center locations (+/-5, +/-10, +/-15 mm), have shown that artifact motion can be substantially decoupled from the global skeletal motion with an effective increase of the accuracy wrt standard techniques. In particular, the comparison between the nominal kinematical variables and the one recovered from markers attached to the skin surface proved LME to be sensibly superior (50% in the worse condition) to the methods imposing marker-bone rigidity. In conclusion, while requiring further validation on real movement data, we argue that the proposed method can constitute an appropriate approach toward the improvement of the human motion estimation.     

Extrinsic flexor muscles generate concurrent flexion of all three finger joints

The role of the forearm (extrinsic) finger flexor muscles in initiating rotation of the metacarpophalangeal (MCP) joint and in coordinating flexion at the MCP, the proximal interphalangeal (PIP), and distal interphalangeal (DIP) joints remains a matter of some debate. To address the biomechanical feasibility of the extrinsic flexors performing these actions, a computer simulation of the index finger was created. The model consisted of a planar open-link chain comprised of three revolute joints and four links, driven by the change in length of the flexor muscles. Passive joint characteristics, included in the model, were obtained from system identification experiments involving the application of angular perturbations to the joint of interest. Simulation results reveal that in the absence of passive joint torque, shortening of the extrinsic flexors results in PIP flexion (80°), but DIP (8°) and MCP (7°) joint extension. The inclusion of normal physiological levels of passive joint torque, however, results in simultaneous flexion of all three joints (63° for DIP, 75° for PIP, and 43° for MCP). Applicability of the simulation results was confirmed by recording finger motion produced by electrical stimulation of the extrinsic flexor muscles for the index finger. These findings support the view that the extrinsic flexor muscles can initiate MCP flexion, and produce simultaneous motion at the MCP, PIP, and DIP joints.     

A videofluoroscopy-based tracking algorithm for quantifying the time course of human intervertebral displacements

The motions of individual intervertebral joints can affect spine motion, injury risk, deterioration, pain, treatment strategies, and clinical outcomes. Since standard kinematic methods do not provide precise time-course details about individual vertebrae and intervertebral motions, information that could be useful for scientific advancement and clinical assessment, we developed an iterative template matching algorithm to obtain this data from videofluoroscopy images. To assess the bias of our approach, vertebrae in an intact porcine spine were tracked and compared to the motions of high-contrast markers. To estimate precision under clinical conditions, motions of three human cervical spines were tracked independently ten times and vertebral and intervertebral motions associated with individual trials were compared to corresponding averages. Both tests produced errors in intervertebral angular and shear displacements no greater than 0.4° and 0.055 mm, respectively. When applied to two patient cases, aberrant intervertebral motions in the cervical spine were typically found to correlate with patient-specific anatomical features such as disc height loss and osteophytes. The case studies suggest that intervertebral kinematic time-course data could have value in clinical assessments, lead to broader understanding of how specific anatomical features influence joint motions, and in due course inform clinical treatments.     

In situ comparison of A-mode ultrasound tracking system and skin-mounted markers for measuring kinematics of the lower extremity

Skin-mounted marker based motion capture systems are widely used in measuring the movement of human joints. Kinematic measurements associated with skin-mounted markers are subject to soft tissue artifacts (STA), since the markers follow skin movement, thus generating errors when used to represent motions of underlying bone segments. We present a novel ultrasound tracking system that is capable of directly measuring tibial and femoral bone surfaces during dynamic motions, and subsequently measuring six-degree-of-freedom (6-DOF) tibiofemoral kinematics. The aim of this study is to quantitatively compare the accuracy of tibiofemoral kinematics estimated by the ultrasound tracking system and by a conventional skin-mounted marker based motion capture system in a cadaveric experimental scenario. Two typical tibiofemoral joint models (spherical and hinge models) were used to derive relevant kinematic outcomes. Intra-cortical bone pins equipped with optical markers were inserted in the tibial and femoral bones to serve as a reference to provide ground truth kinematics. The ultrasound tracking system resulted in lower kinematic errors than the skin-mounted markers (the ultrasound tracking system: maximum root-mean-square (RMS) error 3.44° for rotations and 4.88 mm for translations, skin-mounted markers with the spherical joint model: 6.32° and 6.26 mm, the hinge model: 6.38° and 6.52 mm). Our proposed ultrasound tracking system has the potential of measuring direct bone kinematics, thereby mitigating the influence and propagation of STA. Consequently, this technique could be considered as an alternative method for measuring 6-DOF tibiofemoral kinematics, which may be adopted in gait analysis and clinical practice.     

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.     

Kinematic Analysis and Experimental Verification on the Locomotion of Gecko

This paper presents a kinematic analysis of the locomotion of a gecko,and experimental verification of the kinematicmodel.Kinematic analysis is important for parameter design,dynamic analysis,and optimization in biomimetic robot research.The proposed kinematic analysis can simulate,without iteration,the locomotion of gecko satisfying the constraint conditionsthat maintain the position of the contacted feet on the surface.So the method has an advantage for analyzing the climbing motionof the quadruped mechanism in a real time application.The kinematic model of a gecko consists of four legs based on 7-degreesof freedom spherical-revolute-spherical joints and two revolute joints in the waist.The motion of the kinematic model issimulated based on measurement data of each joint.The motion of the kinematic model simulates the investigated real gecko’smotion by using the experimental results.The analysis solves the forward kinematics by considering the model as a combinationof closed and open serial mechanisms under the condition that maintains the contact positions of the attached feet on the ground.The motions of each joint are validated by comparing with the experimental results.In addition to the measured gait,three othergaits are simulated based on the kinematic model.The maximum strides of each gait are calculated by workspace analysis.Theresult can be used in biomimetic robot design and motion planning.     

Kinematic correlates of walking cadence in the foot

Evidence has frequently been reported of modifications in gait patterns within the lower limb related to the cadence of walking. Most reports have concerned relationships between cadence and kinematic and the kinetic changes occurring in the main joints and muscles of the lower limb as a whole. The aim of the present study was to assess whether significant changes are also measurable in kinematics of the foot segments. An existing 15 marker-set protocol allowed a four-segment foot and shank model to be defined for relative rotations between the segments to be calculated. Stereophotogrammetry was employed to record marker position data from ten subjects walking at three cadences. The slow- and normal cadence datasets showed similar profiles of joint rotation in three anatomical planes, but significant differences were found between these and the fast cadence. At all joints, frame-by-frame statistical analysis revealed increased dorsiflexion from heel-strike to midstance (p<0.05) and increased plantarflexion from midstance to toe-off (p<0.05) with increasing cadence. From foot-flat to heel-rise, the fast cadence kinematic data showed a decreased range of motion in the sagittal-plane between forefoot and rearfoot (3.2°±1.2° at slow cadence; 2.0°±0.8° at fast cadence; p<0.05). The cadences imposed and the multisegment protocol revealed significant kinematic changes in the joints of the foot during barefoot walking.     

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.     

Robust estimation of dominant axis of rotation

A simple method is developed for robustly estimating a fixed dominant axis of rotation (AoR) of anatomical joints from surface marker data. Previous approaches which assume a model of circular marker trajectories use plane fitting to estimate the direction of the AoR. However, when there is limited joint range of motion and rotation due to a second degree of freedom, minimizing only the planar error can give poor estimates of the AoR direction. Optimizing a cost function which includes the error component within a plane, instead of only the component orthogonal to a plane, leads to improved estimates of the AoR direction for joints which exhibit additional rotational motion from a second degree of freedom. Results from synthetic data validation show the ranges of motion where the new method has lower estimation error compared to plane-fitting techniques. Estimates of the flexion-extension AoR from empirical motion capture data of the knee and index finger joints were also more anatomically plausible.     

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.     

An hypothesis for the role of the spine in human locomotion: A challenge to current thinking

Locomotion was first achieved by the motion of the spine. The limbs came after, as an improvement, not as a substitute; and yet, analysis of bipedal gait concentrates almost exclusively on the motion of the limbs. The requirements for land locomotion are examined from a general point of view and the evolution of the vertebrate spine is presented as a mechanism designed to move the animal. The necessary spinal movements are also analysed; the role of the musculoskeletal system is discussed and it is shown that the lumbar spine is a key structure in land locomotion, the pelvis being driven by the spine. The optimum control of motion demands that the stress at all the intervertebral joints should be minimized and equalized. This theory of locomotion requires the central nervous system to control the torque at those intervertebral joints and suggests that a breakdown of the control system would result in torsional failure of the spine. The theory is supported by EMG, force and torque data collected from several sources.     

The three-dimensional kinematics and flexibility characteristics of the human ankle and subtalar joints--Part I: Kinematics

The in-vitro, three dimensional kinematic characteristics of the human ankle and subtalar joint were investigated in this study. The main goals of this investigation were: 1) To determine the range of motion of the foot-shank complex and the associated range of motion of the ankle and subtalar joints; 2) To determine the kinematic coupling characteristics of the foot-shank complex, and 3) To identify the relationship between movements at the ankle and subtalar joints and the resulting motion produced between the foot and the shank. The tests were conducted on fifteen fresh amputated lower limbs and consisted of incrementally displacing the foot with respect to the shank while the motion of the articulating bones was measured through a three dimensional position data acquisition system. The kinematic analysis was based on the helical axis parameters describing the incremental displacements between any two of the three articulating bones and on a joint coordinate system used to describe the relative position between the bones. From the results of this investigation it was concluded that: 1) The range of motion of the foot-shank complex in any direction (dorsiflexion/plantarflexion, inversion/eversion and internal rotation/external rotation) is larger than that of either the ankle joint or the subtalar joint.; 2) Large kinematic coupling values are present at the foot-shank complex in inversion/eversion and in internal rotation/external rotation. However, only a slight amount of coupling was observed to occur in dorsiflexion/plantarflexion.; 3) Neither the ankle joint nor the subtalar joint are acting as ideal hinge joints with a fixed axis of rotation.; 4) Motion of the foot-shank complex in any direction is the result of rotations at both the ankle and the subtalar joints. However, the contribution of the ankle joint to dorsiflexion/plantarflexion of the foot-shank complex is larger than that of the subtalar joint and the contribution of the subtalar joint to inversion/eversion is larger than that of the ankle joint.; 5) The ankle and the subtalar joints have an approximately equal contribution to internal rotation/external rotation movements of the foot-shank complex.