A proposal for a new definition of the axial rotation angle of the shoulder joint.

The Euler/Cardan angles are commonly used to define the motions of the upper arm with respect to the trunk. This definition, however, has a problem in that the angles of both the horizontal flexion/extension and the axial rotation of the shoulder joint become unstable at the gimbal-lock positions. In this paper, a new definition of the axial rotation angle was proposed. The proposed angle was stable over the entire range of the shoulder motion. With the new definition, the neutral position of the axial rotation agreed with that in the conventional anatomy. The advantage of the new definition was demonstrated by measuring actual complex motions of the shoulder with a three-dimensional motion capture system.     

On the estimation of joint kinematics during gait.

In gait analysis, the concepts of Euler and helical (screw) angles are used to define the three-dimensional relative joint angular motion of lower extremities. Reliable estimation of joint angular motion depends on the accurate definition and construction of embedded axes within each body segment. In this paper, using sensitivity analysis, we quantify the effects of uncertainties in the definition and construction of embedded axes on the estimation of joint angular motion during gait. Using representative hip and knee motion data from normal subjects and cerebral palsy patients, the flexion-extension axis is analytically perturbed +/- 15 degrees in 5 degrees steps from a reference position, and the joint angles are recomputed for both Euler and helical angle definitions. For the Euler model, hip and knee flexion angles are relatively unaffected while the ab/adduction and rotation angles are significantly affected throughout the gait cycle. An error of 15 degrees in the definition of flexion-extension axis gives rise to maximum errors of 8 and 12 degrees for the ab/adduction angle, and 10-15 degrees for the rotation angles at the hip and knee, respectively. Furthermore, the magnitude of errors in ab/adduction and rotation angles are a function of the flexion angle. The errors for the ab/adduction angles increase with increasing flexion angle and for the rotation angle, decrease with increasing flexion angle. In cerebral palsy patients with flexed knee pattern of gait, this will result in distorted estimation of ab/adduction and rotation. For the helical model, similar results are obtained for the helical angle and associated direction cosines.(ABSTRACT TRUNCATED AT 250 WORDS)     

A simple method to obtain consistent and clinically meaningful pelvic angles from euler angles during gait analysis

Clinical gait analysis usually describes joint kinematics using Euler angles, which depend on the sequence of rotation. Studies have shown that pelvic obliquity angles from the traditional tilt-obliquity-rotation (TOR) Euler angle sequence can deviate considerably from clinical expectations and have suggested that a rotation-obliquity-tilt (ROT) Euler angle sequence be used instead. We propose a simple alternate approach in which clinical joint angles are defined and exactly calculated in terms of Euler angles from any rotation sequence. Equations were derived to calculate clinical pelvic elevation, progression, and lean angles from TOR and ROT Euler angles. For the ROT Euler angles, obliquity was exactly the same as the clinical elevation angle, rotation was similar to the clinical progression angle, and tilt was similar to the clinical lean angle. Greater differences were observed for TOR. These results support previous findings that ROT is preferable to TOR for calculating pelvic Euler angles for clinical interpretation. However, we suggest that exact clinical angles can and should be obtained through a few extra calculations as demonstrated in this technical note.     

Globographic visualisation of three dimensional joint angles

Three different methods for describing three dimensional joint angles are commonly used in biomechanics. The joint coordinate system and Cardan/Euler angles are conceptually quite different but are known to represent the same underlying mathematics. More recently the globographic method has been suggested as an alternative and this has proved particularly attractive for the shoulder joint. All three methods can be implemented in a number of ways leading to a choice of angle definitions. Very recently Rab has demonstrated that the globographic method is equivalent to one implementation of the joint coordinate system. This paper presents a rigorous analysis of the three different methods and proves their mathematical equivalence. The well known sequence dependence of Cardan/Euler is presented as equivalent to configuration dependence of the joint coordinate system and orientation dependence of globographic angles. The precise definition of different angle sets can be easily visualised using the globographic method using analogues of longitude, latitude and surface bearings with which most users will already be familiar. The method implicitly requires one axis of the moving segment to be identified as its principal axis and this can be extremely useful in helping define the most appropriate angle set to describe the orientation of any particular joint. Using this technique different angle sets are considered to be most appropriate for different joints and examples of this for the hip, knee, ankle, pelvis and axial skeleton are outlined.     

Determination of 3D spinal kinematics without defining a local vertebral coordinate system

In this paper a method is presented to calculate Euler's angles of rotation of a body segment during locomotion without a priori defining the location of the center of rotation, and without defining a local vertebral coordinate system. The method was applied to in vivo spinal kinematics. In this method, the orientation of each segment is identified by a set of three markers. The orientation of the axes of rotation is calculated based on the average position of the markers during one stride cycle. Some restrictions and assumptions should be made. The approach is viable only when the average orientation of the anatomical axes of rotation of each spinal segment during a stride cycle coincides with the three axes of the laboratory coordinate system. Furthermore, the rotations should be symmetrical with respect to both sides of the plane of symmetry of the spinal segment, and the subject should move parallel to one axis of the laboratory coordinate system. Since in experimental conditions these assumptions will only be met approximately, errors will be introduced in the calculated angles of rotation. The magnitude of the introduced errors was investigated in a computer simulation experiment. Since the maximal errors did not exceed 0.7° in a range of misalignments up to 10° between the two coordinate systems, the approach proved to be a valid method for the estimation of spinal kinematics.     

A 3D Generic Inverse Dynamic Method using Wrench Notation and Quaternion Algebra

In the literature, conventional 3D inverse dynamic models are limited in three aspects related to inverse dynamic notation, body segment parameters and kinematic formalism. First, conventional notation yields separate computations of the forces and moments with successive coordinate system transformations. Secondly, the way conventional body segment parameters are defined is based on the assumption that the inertia tensor is principal and the centre of mass is located between the proximal and distal ends. Thirdly, the conventional kinematic formalism uses Euler or Cardanic angles that are sequence-dependent and suffer from singularities.

In order to overcome these limitations, this paper presents a new generic method for inverse dynamics. This generic method is based on wrench notation for inverse dynamics, a general definition of body segment parameters and quaternion algebra for the kinematic formalism.     

A 3D generic inverse dynamic method using wrench notation and quaternion algebra

In the literature, conventional 3D inverse dynamic models are limited in three aspects related to inverse dynamic notation, body segment parameters and kinematic formalism. First, conventional notation yields separate computations of the forces and moments with successive coordinate system transformations. Secondly, the way conventional body segment parameters are defined is based on the assumption that the inertia tensor is principal and the centre of mass is located between the proximal and distal ends. Thirdly, the conventional kinematic formalism uses Euler or Cardanic angles that are sequence-dependent and suffer from singularities. In order to overcome these limitations, this paper presents a new generic method for inverse dynamics. This generic method is based on wrench notation for inverse dynamics, a general definition of body segment parameters and quaternion algebra for the kinematic formalism.     

Use of dual Euler angles to quantify the three-dimensional joint motion and its application to the ankle joint complex

This paper presents a modified Euler angles method, dual Euler angles approach, to describe general spatial human joint motions. In dual Euler angles approach, the three-dimensional joint motion is considered as three successive screw motions with respect to the axes of the moving segment coordinate system; accordingly, the screw motion displacements are represented by dual Euler angles. The algorithm for calculating dual Euler angles from coordinates of markers on the moving segment is also provided in this study. As an example, the proposed method is applied to describe motions of ankle joint complex during dorsiflexion–plantarflexion. A Flock of Birds electromagnetic tracking device (FOB) was used to measure joint motion in vivo. Preliminary accuracy tests on a gimbal structure demonstrate that the mean errors of dual Euler angles evaluated by using source data from FOB are less than 1° for rotations and 1 mm for translations, respectively. Based on the pilot study, FOB is feasible for quantifying human joint motions using dual Euler angles approach.     

The analysis of golf swing as a kinematic chain using dual Euler angle algorithm

The manner in which anatomical rotation from an individual segment contributes to the position and velocity of the endpoint can be informative in the arena of many athletic events whose goals are to attain the maximal velocity of the most distal segment. This study presents a new method of velocity analysis using dual Euler angles and its application in studying rotational contribution from upper extremity segments to club head speed during a golf swing. Dual Euler angle describes 3D movement as a series of ordered screw motions about each orthogonal axis in a streamlined matrix form-the dual transformation matrix- and allows the translation and rotation component to be described in the same moving frame. Applying this method in biomechanics is a novel idea and the authors have previously applied the methodology to clinical studies on its use in displacement analysis. The focus of this paper is velocity analysis and applications in sports biomechanics. In this study, electrogoniometers (Biometrics, UK) with a frequency of 1000 Hz were attached to a subject during the execution of the swing to obtain the joint angles throughout the motion. The velocity of the club head was then analyzed using the dual velocity which specifies the velocity distribution of a rigid body in screw motion at any point in time as the dual vector. The contributions of each segment to the club-head velocity were also compared. In order to evaluate this method, the calculated position and velocity of the club head were compared to the values obtained from video image analysis. The results indicated that there is good agreement between calculated values and video data, suggesting the suitability of using the Dual Euler method in analyzing a kinematic chain motion.     

The envelope of passive knee joint motion

The purpose of this study is to create an accurate experimental database for the passive (in vitro) freedom-of-motion characteristics of the human knee joint on a subject to subject basis, suitable for the verification and enhancement of mathematical knee-joint models. Knee-joint specimens in a six degree-of-freedom motion rig are moved through flexion under several combinations of external loads, including tibial torques, axial forces and AP-forces. Euler rotation angles and translation vectors, describing the relative, spatial motions of the joint are measured using an accurate Roentgen Stereo Photogrammetric system. Conceptually the joint is considered as a two degrees-of-freedom of motion mechanism (flexion-tibial rotation), whereby the limits of internal and external tibial rotation are defined at torques of +/- 3 Nm. The motion pathways along these limits are defined as the envelopes of passive knee joint motion. It is found that these envelope pathways are consistent and hardly influenced by additional axial forces up to 300 N and AP-forces of 30 N. Within the envelope of motion, however, the motion patterns are highly susceptible to small changes in the external load configuration. It is shown that the external tibial rotation during extension ('screw-home mechanism') is not an obligatory effect of the passive joint characteristics, but a direct result of the external loads. Anatomical differences notwithstanding, the inter-individual discrepancies in the motion patterns of the four specimens tested, showed to be relatively small in a qualitative sense. Quantitative differences can be explained by small differences in the alignment of the coordinate systems relative to the joint anatomy and by differences in rotatory laxity.     

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.     

Rationalization of kinematic descriptors for three-dimensional hand and finger motion

Modelling joint motion in three dimensions is often based on techniques taken from classical dynamics, each analysis resulting in a set of six parameters describing the relative motion betwen two body segments. The literature on joint kinematics has been difficult to compare due to use of different anatomical landmarks, axis nomenclature, and analytical methods. It is here shown that with care in sequence definition, the three alignment-based systems (Euler, Cardan, floating axis) give identical results for angular parameters. While the equivalent screw displacement axis system can be related simply to the other methods only if the functional axis of motion is aligned with a coordinate axis, the basic matrix for relating rigid body positions before and after a motion can always be reconstructed. Therefore the changes in alignment angles may be obtained from screw displacement parameters, permitting the results of different analyses to be compared. Translation parameters are most difficult to interpret in any system. Examples of the way in which simple planar motions are characterized by the various analytical methods are given.     

A novel method for assessing the 3-D orientation accuracy of inertial/magnetic sensors

A novel method for assessing the accuracy of inertial/magnetic sensors is presented. The method, referred to as the “residual matrix” method, is advantageous because it decouples the sensor's error with respect to Earth's gravity vector (attitude residual error: pitch and roll) from the sensor's error with respect to magnetic north (heading residual error), while remaining insensitive to singularity problems when the second Euler rotation is close to ±90°. As a demonstration, the accuracy of an inertial/magnetic sensor mounted to a participant's forearm was evaluated during a reaching task in a laboratory. Sensor orientation was measured internally (by the inertial/magnetic sensor) and externally using an optoelectronic measurement system with a marker cluster rigidly attached to the sensor's enclosure. Roll, pitch and heading residuals were calculated using the proposed novel method, as well as using a common orientation assessment method where the residuals are defined as the difference between the Euler angles measured by the inertial sensor and those measured by the optoelectronic system. Using the proposed residual matrix method, the roll and pitch residuals remained less than 1° and, as expected, no statistically significant difference between these two measures of attitude accuracy was found; the heading residuals were significantly larger than the attitude residuals but remained below 2°. Using the direct Euler angle comparison method, the residuals were in general larger due to singularity issues, and the expected significant difference between inertial/magnetic sensor attitude and heading accuracy was not present.     

Interpolation of segment Euler angles can provide a robust estimation of segment angular trajectories during asymmetric lifting tasks

Video-based field methods that estimate L5/S1 net joint moments from kinematics based on interpolation in the sagittal plane of joint angles alone can introduce a significant error on the interpolated joint angular trajectory when applied to asymmetric dynamic lifts. Our goal was to evaluate interpolation of segment Euler angles for a wide range of dynamic asymmetric lifting tasks using cubic spline methods by comparing the interpolated values with the continuous measured ones. For most body segments, the estimated trajectories of segment Euler angles have less than 5° RMSE (in each dimension) with 5-point cubic spline interpolation when there is no measurement error of interpolation points. Sensitivity analysis indicates that when the measurement error exists, the root mean square error (RMSE) of estimated trajectories increases. Comparison among different lifting conditions showed that lifting a load from a high initial position yielded a smaller RMSE than lifting from a low initial position. In conclusion, interpolation of segment Euler angles can provide a robust estimation of segment angular trajectories during asymmetric lifting when measurement error of interpolation points can be controlled at a low level.     

A kinematic model of the human ankle

The spatial gross motion of the foot with respect to the shank is modelled as rotations about two fixed ankle axes: the upper ankle rotation axis (plantarflexion/dorsiflexion) and the subtalar rotation axis (inversion/eversion). The positions of the axes are determined by externally visible bony landmarks of the lower leg and are measured for a living subject. The model input data are the plantarflexion/dorsiflexion and inversion/eversion rotation angles; the model output is a 4 × 4 transformation matrix which quantitatively describes the relative position of a foot coordinate system with respect to a shank coordinate system.     

Comparison of three local frame definitions for the kinematic analysis of the fingers and the wrist

Because the hand is a complex poly-articular limb, numerous methods have been proposed to investigate its kinematics therefore complicating the comparison between studies and the methodological choices. With the objective of overcoming such issues, the present study compared the effect of three local frame definitions on local axis orientations and joint angles of the fingers and the wrist. Three local frames were implemented for each segment. The “Reference” frames were aligned with global axes during a static neutral posture. The “Landmark” frames were computed using palpated bony landmarks. The “Functional” frames included a flexion–extension axis estimated during functional movements. These definitions were compared with regard to the deviations between obtained local segment axes and the evolution of joint (Cardan) angles during two test motions. Each definition resulted in specific local frame orientations with deviations of 15° in average for a given local axis. Interestingly, these deviations produced only slight differences (below 7°) regarding flexion–extension Cardan angles indicating that there is no preferred method when only interested in finger flexion–extension movements. In this case, the Reference method was the easiest to implement, but did not provide physiological results for the thumb. Using the Functional frames reduced the kinematic cross-talk on the secondary and tertiary Cardan angles by up to 20° indicating that the Functional definition is useful when investigating complex three-dimensional movements. Globally, the Landmark definition provides valuable results and, contrary to the other definitions, is applicable for finger deformities or compromised joint rotations.     

Simulation of Codman's paradox reveals a general law of motion

Codman's paradox refers to a specific pattern of motion at the shoulder joint. It asked how a mysterious axial rotation about the longitudinal axis of the arm occurred during two or three sequential arm rotations that did not involve rotation about the long-axis. The objective of this paper was to find how the mysterious axial rotation occurred in the Codman's paradox. First, Codman's paradox and Codman's rotation were defined in general situations that involved arm rotations about orthogonal axes starting from the neutral attitude as well as rotations about non-orthogonal axes and starting from non-neutral attitudes. Then a general law of motion was proposed to answer the question of the Codman's paradox, which is stated as when the long-axis of the arm performs a closed-loop motion by three sequential rotations defined as Codman's rotation, it produces an equivalent axial rotation angle about the long-axis. The equivalent axial rotation angle equals the angle of swing-the second rotation in the three sequential long-axis rotations. Validity of the proposed law of motion is demonstrated by computer simulation of various Codman's rotations. Clinical relevance of the proposed law of motion is also discussed in the paper.     

A correction for axis misalignment in the joint angle curves representing knee movement in gait analysis

This paper derives a simple mathematical model relating changes in the orientations of the two Cartesian coordinate systems involved in recording knee movement and the varus-valgus and the internal-external rotation angles for describing the knee's motion. Rotation matrices are given for changing the orientations of the two Cartesian coordinate systems in such a way that the quadratic variations in the varus-valgus and in the external-internal angles are minimal. These estimated rotation matrices are used to correct for axis misalignment. The correction is calibrated by considering the impact of the new orientation of the thigh Cartesian coordinate system on the hip joint angles. The procedure is applied to kinematic data collected on normal subjects. The uncertainty about the specification of the thigh Cartesian coordinate system is shown to explain some of the between subject variability in the varus-valgus and in the internal-external rotation angles curves.