Orientation of axes in the elbow and forearm for biomechanical modelling

H.E.J. (DirkJan) Veeger1, B. Yu2

1Institute for Fundamental and Clinical Human Movement Sciences, Department of Human Movement Sciences, VU, 1081 BT Amsterdam, the Netherlands

2Orthopaedic Biomechanics Laboratory, Mayo Clinic, Rochester MN 55905


Published in:

Bajpaj P.K. (ed.) Proceedings of the Fifteenth Southern Biomedical Engineering Conference,

Dayton, Ohio, pp. 377-380


ACKNOWLEDGMENT

This study was supported by a NATO-Science Fellowhip awarded by the Netherlands Organisation for Scientific Research (NWO) and by NIH grants HD07447 and AR41171.
 

Abstract - To determine the three-dimensional positions and orientations of axes of rotation for elbow flexion and forearm pro-supination, the flexion-extension and pro-supination movements were measured for five arms. Four right and one left arm were taken from four fresh cadavers. Movement was measured with an electromagnetic tracking system. Axes were found to be almost perpendicular (88.9 ± 5.1°) and crossed at 3.3 ± 0.8 mm. It is suggested that arm movements can be modelled with a two Degree-Of-Freedom model.

I. INTRODUCTION

Biomechanical modelling of the upper extremity can play an important role in understanding the causes of upper extremity disorders and the corresponding treatment. Also, a biomechanical model can be used to get insight into more fundamental issues on motor control. Recently, a three-dimensional (3-D) model of the shoulder has been developed [1]. This model is based on morphological data that were collected in an extensive cadaver study [2]. However, this study did not embrace data for extension of the model to the elbow and the elbow function had to be modelled with a three Degree-of-Freedom (DOF) joint [1]. To improve the range of the model and allow for more accurate modelling of the function of the bi-articular arm muscles, the addition of a morphologically more realistic description of the arm was judged to be necessary. Since previous studies generally reported qualitative descriptions of elbow flexion axes [3,4,5,6] and forearm pro-supination axes [3,4,7], those data were not available from the literature. The aim of this project was to collect quantitative information on the 3-D orientations of rotation axes in the elbow and forearm.

II. METHODS

Five upper extremities were harvested from four fresh specimen at the level of the scapula (four right, one left). Subsequently, the extremities were fastened onto a measuring table such that the scapula was fixed and the arm was fully free to move. Following fixation of 3-D electromagnetic sensors (Isotrack, Polhemus) on scapula, humerus, ulna and radius, the positions of local anatomical landmarks were digitized with an extra sensor, mounted with a stylus. The anatomical landmarks were at a later stage used for the definition of local orientation axes. Subsequently, each arm was moved through a selection of standard directions (elbow flexion-extension, forearm pro-supination and glenohumeral movements), during which the positions and orientations of each sensor were collected. A typical movement consisted of two cycles in which a joint was moved through its full range of motion. The sampling period was 15 seconds with a sampling frequency of Fs=10Hz.

Sensor data were used to determine the rotation axes for elbow flexion and pro/supination of the forearm. Estimations of the rotation axes were based on the Instantaneous Helical Axes (IHA) algorithms [8]. Prior to calculation of the instantaneous helical axes, for each sensor the position and angle time series were filtered with a cut-off frequency of 1.5 Hz. The angular velocity of each sensor was calculated using the method described by Woltring [9]. The Instantaneous Helical Axis was calculated for each time sample where the angular velocity of the distal segment, relative to the proximal segment was larger than 0.25 rad.s -1. To determine a mean helical axis and pivot, the optimal pivot point was calculated as sopt, as the mean 'pivot', closest to all Instantaneous Helical Axes in a least-squared sense [3]:

BME_1.gif (1) where: bme_2.gif, and BME_3.gif (2)

v = unit direction vector;

s = projection of the sensor origin onto the IHA.

Fig. 1 Typical example for IHA calculations for elbow flexion, represented in the measurement co-ordinate system. The X-Z plane roughly corresponds with the frontal plane.

IHA are drawn with a length of 2 cm. Dots are data samples from the surfaces of humerus, ulna and radius. The two asterisks, connected with a dashed line, form the line EL-EM.

Fig. 2 Typical example for IHA calculations for forearm pronation, represented in the measurement co-ordinate system. The X-Z plane roughly corresponds with the frontal plane.

IHA are drawn with a length of 2 cm. Dots are data samples from the surfaces of humerus, ulna and radius. The two asterisks, connected with a dashed line, form the line EL-US.

Analogous to the position vector, the optimal direction vector vopt wascalculated as the vector with the smallest angle between axes. For each specimen, these axes were subsequently expressed relative to a local co-ordinate system with its origin in the Epicondylus Medialis (EM) and its X-axis in the direction of the Epicondylu Lateralis (EL), and processus styloideus ulnaris (US) in the X-Y plane.

To allow for calculation of a mean axis over all five specimen, each arm was mathematically set in a theoretical anatomical position in which the long axes of the humerus and ulna both lie within the same plane, perpendicular to the elbow axis.. This position was obtained through a rotation of the humerus around the individual elbow axes, to the same plane as the ulna (Z=0). In this position, the elbow flexion angle is zero. The rotation to the theoretical anatomical positions was performed on all landmarks on the humerus: the Acromion (AC), EL and EM.

TABLE I.

Optimal Axes for Flexion-Extension and Pro-Supination, estimated for five Specimen (four right Arms and one left Arm). Data are given in Cm and are expressed in their local Ulnar Co-Ordinate Systems., with EM as the Origin. e Is calculated as the mean absolute Error.

Flexion Soptx  Sopty Soptz e (cm) Voptx Vopty Voptz e (deg)
#1_r 3.58 0.73 0.99 0.65 -0.9990 0.0099 0.0434 1.61
#2_r 4.77 0.47 0.70 0.78 -0.9907 0.1326 0.0303 2.59
#3_r 4.31 0.79 0.15 0.63 -0.9975 0.0709 0.0081 1.86
#4_r 4.73 0.64 0.64 0.90 -0.9881 0.0332 -0.1505 2.12
#2_l 2.42 0.99 -0.72 1.27 -0.9935 0.0952 -0.0622 1.45

Pronation
#1_r 6.27 7.94 0.63 0.24 0.0579 0.9982 0.0136 0.40
#2_r 6.38 4.81 0.80 0.43 -0.0125 0.9998 -0.0127 1.43
#3_r 6.75 11.62 -0.14 0.59 0.0605 0.9981 0.0068 1.24
#4_r 6.87 10.87 0.56 0.40 0.1059 0.9942 0.0194 0.73
#2_l 6.15 4.91 -0.66 0.55 0.0289 0.9990 0.0335 0.91

III. RESULTS AND DISCUSSION

Figures 1 and 2 illustrate the helical axes as determined for elbow flexion and forearm pronation. The data are plotted in the measuring

co-ordinate system. IHA calculations showed fluctuations in position, but not as much in orientation. The mean smallest distances between the IHA axes and the optimal axes were generally well below 1 cm, while the mean absolute angle stayed below 2.6° (Table 1).

The individual optimal axes for flexion-extension and pro-supination are given in Table 1. These axes are expressed in the local ulnar system with the origin in EM and the local X-axis in the direction of EL, but with the arm in the originally measured (anatomical) position. In this position, the arm is in approximately 17.5° flexion. From the calculated error in Table 1, it can be seen that the pronation of the forearm appears to take place around a tighter axis than for elbow flexion. The kinematic flexion-extension axes deviate 6.0±2.6° from the axes through EM and EL. The pro-supination axis crosses the flexion-extension axis at 3.3 ± 0.8 mm and an angle of 88.9 ± 5.1°. This axis runs through the radial head and lies close to the anatomical landmarks EL (13.1±2.2 mm) and US (8.0±4.5 mm). The average distance to the center of the radial head was found to be 3.2±4.0 mm.

The mean elbow flexion-extension and forearm pro-supination axes have also been determined in the (theoretical) anatomical position where the elbow flexion angle is zero. Table 2 contains the average data for all five specimen. The data associated with the humerus (EL, AC, Huhe) have mathematically been transformed to a position in which the elbow angle is zero, which is the theoretical anatomical position. The data for AC and the estimated center of the humeral head indicate that the orientations of the humerus relative to the ulna was relatively constant over the five specimen. Obviously, this procedure still results in differences in position for the landmarks on humerus and ulna, since the specimen differ in build and in segment orientation.

The findings of this study confirm that the elbow-forearm complex can be modelled as a two DOF system. The orientations of the estimated elbow axes confirm previous qualitative observations that describe a flexion-extension axis as running through the centre of the trochlea and the capitulum humeri [3,4,5,6]. Also, the pro-supination axes measured in this study confirmed earlier descriptions [3,4,7] on an orientation of that axis through the center of the radial head and the most distal point on the ulna (which is approximated by US).

The results of this study are being implemented in a model of the shoulder and arm. First results indicate a good fit between in-vivo measured arm movements during wheelchair propulsion and kinematic model simulations.

TABLE II.

Mean (N=5) Axis Orientations for Elbow Flexion-Extension and Forearm Pro-Supination, after rotation around the Elbow Axis to align the Ulna and Radius. The Mean Axes are calculated over Five Specimen.

Sflexion = Position Vector for Flexion-Extension, Vflexion = Unit Direction Vector for Flexion-Extension. AC = Acromion, Hum. Head = Center of the Humeral Head, EL = Epicondylus Lateralis, US = Proc. Styloideus Ulnaris.

Mean (N=5) X (cm) (std) Y (cm) (std) Z (cm) (std)
Hum head 4.33 1.38 -29.52 1.53 0.00 0.00
AC 3.22 1.53 -34.58 1.73 --0.70 1.82
EL 7.13 0.11 0.00 0.00 0.00 0.00
US 8.79 1.84 26.79 0.93 0.00 0.00
Sflexion 3.96 0.99 0.52 0.24 0.44 0.82
Vflexion -0.9759 0.0173 0.1364 0.0860 -0.0392 0.1635
Spronation 6.97 0.67 7.62 3.25 0.31 0.61
Vpronation 0.1177 0.0511 0.9918 0.0069 0.0078 0.0176

REFERENCES

[1] F.C.T. Van der Helm "A finite element musculoskeletal model of the shoulder mechanism". J. Biomech. 27(5); 551-570.,1994

[2] H.E.J. Veeger, F.C.T. van der Helm, L.H.V. van der Woude, G.M. Pronk, R.H. Rozendal "Inertia and muscle contraction parameters for musculoskeletal modelling of the shoulder mechanism" .J. Biomech. 24(7); 615-629, 1991.

[3] R. Fick Handbuch der Anatomie und Mechanik der Gelenke. Jena: Gustav Fischer, 1911.

[4] B.F. Morrey, E.Y.S. Chao, "Passive motion of the elbow joint". I. Bone Joint Surg. 58-A; 501-508, 1976.

[5] J.T. Deland, A. Garg, P.S .Walker, "Biomechanical basis for elbow hinge-distractor design" Clin. Orthop. Rel. Res. 215; 303-312, 1987.

[6] Y. Youm, R.F. Dryer, K. Thambyrajah, A.E. Flatt, B.L. Sprague, "Biomechanical analyses of forearm pronation-supination and elbow flexion-extension." J. Biomech. 12; 245-255, 1979

[7] A.M. Hollister, H. Gellman, R.L. Waters "The relationship of the interosseous membrane to the axis of rotation of the forearm." Clin. Orth. Rel. Res. 298; 272-276, 1994.

[8] H.J. Woltring, "Data processing and error analysis." In: Biomechanics of human movement: Applications in Rehabilitation, Sport and Ergonomics. (Edited by Capozzo, A and Berme P.) 10.1, pp 203-237, 1990.

[9] H.J. Woltring, "Definition and Calculus of Attitude Angles and Instantaneous Helical Axes from noisy Position and Attitude data". In: Proc. on the Int. Symp. on 3-D Analysis of human Movement, Montreal, 59-62, 1991.